- Open Access
Can metrics of acceleration provide accurate estimates of energy costs of locomotion on uneven terrain? Using domestic sheep (Ovis aries) as an example
Animal Biotelemetry volume 10, Article number: 8 (2022)
Locomotion is often a necessity for animal survival and can account for a large proportion of an individual’s energy budget. Therefore, determining the energy costs of locomotion is an important part of understanding the interaction between an animal and its environment. Measures of animal acceleration, specifically ‘dynamic body acceleration’ (DBA) has proved to be a useful proxy of the energy cost of locomotion. However, few studies have considered the effects of interacting factors such as the animal’s speed or changes to the terrain slope on the putative acceleration versus energy expenditure relationship and how this may affect the relationship between DBA and energy expenditure.
Here we conducted a methodological study to evaluate the ability of the metric ‘vectorial dynamic body acceleration’, VeDBA, obtained from tri-axial accelerometer data loggers, to act as a proxy for energy expenditure in non-uniform environments. We used indirect calorimetry to measure the oxygen consumption (V̇O2) of domestic sheep (Ovis aries) that were exposed to different ambient temperatures when immobile (resting) and that walked at various speeds (0.8 to 2.9 km h−1) and slope angles (− 6° to 6°) on a treadmill while simultaneously measuring tri-axial acceleration recorded at 40 Hz by body-mounted tags.
The lower critical temperature of sheep was identified as 18 °C, and V̇O2 when they were immobile was 3.67 mL O2 kg−1 min−1. There were positive relationships between V̇O2, VeDBA, and speed of walking. However, VeDBA correlated less well with V̇O2 when the terrain slope either inclined or declined.
We advocate caution when using DBA metrics for establishing energy use in animals moving over uneven terrain and suggest that each study species or location must be examined on a case-by-case basis. Reliance upon the relationship described between acceleration and energy expenditure on horizontal-surface treadmills can lead to potential under- or over-estimates of energy expenditure when animals walk on uneven or inclined ground.
Locomotion occupies a pivotal position in animal ecology, because it is necessary in many species for, inter alia, food-finding, searching for mates and avoiding predators . However, translocation of the body is energetically costly, which explains why locomotion accounts for a large proportion of the energy budget of many animals [2,3,4]. Studies on the energetics of animal locomotion often seek to determine how movement costs affect fitness and survival [5, 6] to understand or predict how individuals and populations can survive .
Measurement of body movement via tri-axial accelerometery was first considered in relation to energy expenditure by Cavagna, Saibene, and Margaria . Meijer, Westerterp and Koper  proposed that acceleration-based metrics might relate to power use, while Terrier, Aminian and Schutz  related the mean of the integral of the vector magnitude of tri-axial accelerations measured by devices placed on human backs to oxygen consumption (V̇O2) during locomotion. However, it was not until work by Wilson et al.  that an acceleration-based metric, ‘dynamic body acceleration’ (DBA), was formalized and quantified for moving animals, initially cormorants (Phalacrocorax carbo) walking on a treadmill. The positive linear relationship between DBA and V̇O2 found by Wilson et al.  was echoed subsequently by a suite of authors working on several different taxa [12,13,14,15,16,17,18,19,20].
Despite the ease with which livestock can sometimes be studied compared to wild or non-domesticated animals, applications of tri-axial accelerometers to estimate the cost of locomotion of livestock are relatively rare (e.g., ). A study by Miwa et al.  measured DBA and heart rate in cattle, sheep and goats during locomotion of various intensities. While strong correlations were found, as the authors suggest, direct measurements of V̇O2 would provide a better test of the relationship in ruminants between DBA and energy expenditure. The relationships between DBA and V̇O2 derived from animals walking at various speeds on treadmills (e.g., [12, 13, 18, 23]) have faced criticism, as level treadmill-based estimates of energy expenditure may not reflect energy expenditure in the field . Specifically, the natural environment is rarely uniform, and changes in terrain substrate have been shown to influence the energetic cost associated with locomotion. Soft substrates, such as snow  and sand , incur a greater energetic cost than harder substrates. Importantly, variation in terrain angle can affect energy expenditure; e.g., increased costs of locomotion with increasing steepness of incline have been reported in a number of animals [27,28,29,30,31] due to the mechanical work required to move the body against gravity . Moreover, when the angle of the terrain decreases beyond level past a certain point, work is required to stabilise the body and control velocity of movement, which requires more energy than would be required when walking on a level surface or on shallow declines [29, 30, 33].
Much of the available data that defines relationships between V̇O2 and acceleration in terrestrial animals was conducted on a flat surface, with fewer studies considering the effect of incline [19, 31] and even less considering decline . In humans, the V̇O2-acceleration relationship can be markedly different between flat and angled surfaces. Terrier, Aminian and Schutz  found that energy expenditure increased with increasing incline, and the relationships between acceleration and energy expenditure were only apparent when data from each incline were considered independently. Similarly, Halsey et al.  found that DBA measured using tri-axial accelerometry could be used to predict the energy costs of walking at different inclines on a treadmill, but the V̇O2–DBA relationships were markedly different in each case. These different relationships need to be quantified, particularly for species that live typically in rugged terrains, such as sheep. Such terrains might vary in the steepness of slope from low (1–12°), to medium (13–25°) to high (> 25°) [34, 35] and also in evenness (where ‘uneven terrain’ refers to regions, where changes in terrain angle are frequent and localised). In the current study we explore the effect of slope. Once terrain slope is known and an appropriate V̇O2–DBA relationship for this slope has been determined, we can calculate the energetic cost of movement at this angle. Body pitch (\(\theta\)) of the animal (the gradient of the spine in a quadruped)  could provide a useful continuous measurement of terrain slope. In other words, the slope of the terrain can be expected to be reflected in the pitch of the animal. As a baseline against which to compare the energy costs of specific movements or individual behaviours, we need to determine energy expenditure, while the animal is stationary. Species such as sheep spend a lot of time standing, i.e., in a locomotor stance but travelling at zero speed and effectively resting, and thus this behaviour will represent a considerable part of their overall energy expenditure.
In the present study, we aimed to examine the relationship between V̇O2 and DBA using trained domestic sheep (Ovis aries) as a model quadruped, walking at different speeds and on different slope angles (both inclines and declines), to explore how V̇O2 relates to acceleration under differing circumstances. Ultimately, this work can enable us to understand if estimates of oxygen consumption and a measure of acceleration determined under a controlled environment can be used with confidence to predict how the topography of the environment affects the energy costs of locomotion of free ranging animals. Sheep are common as livestock species across the world, with economic importance both in developed and developing countries. In the UK (population 32.7 million sheep, ) the industry was worth c. $1.3 billion in 2020 . In Australia, the sheep meat industry was worth an estimated $4.5 billion in 2019–2020  and in China this was an estimated $106.7 in 2020 . Sheep are kept in a variety of environments from lowland farms to rugged mountainous terrain and exposed to different farming practices, such as rotational grazing during which sheep are moved from one pasture to another and fleece shearing. Clearly fleece removal will affect body insulation, heat loss and thermoregulatory capabilities, including, presumably, an increase in lower critical temperature (LCT) - i.e., the temperature below which the animal expends additional energy to keep warm . A secondary objective, therefore, was to determine the relationship between ambient temperature and V̇O2 when sheep were stationary. Consequently, as well as contributing to a general framework examining the validity of DBA metrics as a proxy for energy expenditure, this work will introduce terrain as an energy modulator into perspective for a species of commercial importance. The specific aims were to: (i) investigate the relationship between resting V̇O2 and ambient temperature; (ii) examine the effects of speed and angle of incline/decline upon V̇O2; (iii) investigate the effects of speed, angle and pitch on VeDBA; (iv) define a relationship between V̇O2 and VeDBA for sheep across various slope angles; (v) determine if the measurement of animal pitch (\(\theta\)) can be used as a suitable proxy of terrain angle; and (vi) examine the effect of fleece shearing on V̇O2.
Resting V̇O2 and the thermoneutral zone
Piecewise linear regression analysis indicated an inflection in V̇O2 at 18 °C for sheep A, representing the lower critical temperature (LCT). Two inflection points were identified for sheep B at 12.3 °C and 18 °C (Fig. 1B). However, above 12.3 °C, V̇O2 continued to decrease, by 22% to 3.33 mL kg−1 min−1 at 18 °C, beyond which V̇O2 stabilised. The LCT was, therefore, identified as 18 °C for sheep B, which was consistent with sheep A (Fig. 1B). The mean resting V̇O2 for the two sheep was 3.67 mL O2 kg−1 min−1 when they were un-sheared (i.e., with a full fleece). When sheared, V̇O2 decreased by 15% from a mean of 4.78 to 4.05 mL O2 kg−1 min−1 as ambient temperature increased from 23.5 to 27 °C (Fig. 1B). There was no significant difference in the resting V̇O2 of the sheep when they were not sheared compared with when they were sheared, at 27 °C (X2 = 2.106, df = 1, P = 0.147, Fig. 1B). As per the V̇O2 of non-sheared sheep, V̇O2 did not change when temperature increased above 18 °C for sheared sheep.
V̇O2 and locomotion
There was a significant positive relationship between animal pitch and slope angle (X2 = 4.949, df = 1, P = 0.026). V̇O2 increased significantly with both increasing slope angle (Fig. 2C) and treadmill speed (X2 = 6.299, df = 1, P = 0.012 and X2 = 4.534, df = 1, P = 0.033, respectively) (Fig. 2D). Inclusion of a quadratic rather than a linear function improved the fitted relationship (r2 = 0.33 versus r2 = 0.17, respectively) (Fig. 2C).
Factors affecting VeDBA
VeDBA increased significantly with increasing treadmill speed (X2 = 67.75, df = 1, P < 0.001) (Fig. 2B). However, there was no significant change in VeDBA with treadmill slope angle at any of the measured speeds (X2 = 0.348, df = 1, P = 0.55) or pitch of the animal’s body (X2 = 2.3067, df = 1, P = 0.129).
Relationship between V̇O2 and VeDBA
The GLMM model (Table 1, Model 1a), which incorporated data obtained from all slope angles and all speeds, revealed no significant relationship between V̇O2 and VeDBA (X2 = 2.06, df = 1, P = 0.107). Therefore, in subsequent analyses each slope angle was considered independently. The only significant relationship between V̇O2 and VeDBA was noted when the treadmill was level (0°), Spearman rs2 = 0.88, P = 0.007, (Fig. 3). The least-squares relationship was defined:
(Spearman rs2 = 0.77). No significant relationships were noted between V̇O2 and VeDBA for treadmill gradients of 4°, 6°, − 4° and − 6° (rs = 1, 0.8, 0.6, 0.8, P = 0.08, 0.33, 0.42, 0.33, respectively). Small sample size effectively precluded statistical power for significance (cf. rs = 1, above). Grouping data into incline walking (4° and 6°), level walking (0°) and decline walking (− 4° and − 6°) resulted in significant relationships between V̇O2 and VeDBA for both incline walking and level walking (rs = 0.86, P = 0.01, rs = 0.88, P = 0.007, respectively), although there was variation observed around the relationship defined for incline walking (Fig. 3). The relationship defined for incline walking was:
(Least squares regression, Spearman rs2 = 0.74). There was no significant relationship between V̇O2 and VeDBA for decline walking (rs2 = 0.64, P = 0.096).
This study sought to understand the relationships between DBA recorded by body-mounted accelerometer tags, and simultaneously measured energy expenditure using domestic sheep as a model quadruped. It also sought to determine whether DBA–V̇O2 relationships differed according to speed of travel and incline of the terrain. Furthermore, the study aimed to evaluate the potential use of animal pitch as a useful proxy for monitoring the incline of the terrain.
Relationship between DBA and V̇O2
A positive relationship between DBA and energy expenditure has been reported previously in a range of species (e.g., [11, 13, 16, 17, 19, 22]). Results of the current work are in broad agreement with these studies, with a positive relationship noted between DBA and V̇O2 when subjects walked on a level surface. However, no significant relationship emerged when data from level, positive and negative treadmill angles were included together in the analysis. Terrier, Aminian and Schutz  also found no relationship between acceleration metrics and energy expenditure for humans when multiple gradients were examined. They did, however, note that significant relationships existed when gradients were analysed separately. Similarly, when data from the current study were grouped into “incline”, “level” or “decline” walking, relationships between DBA and V̇O2 showed that “incline” walking was more energetically costly than “level” walking and “decline” walking. The fact that there were fewer data in the current study for sheep walking on slopes, compared to sheep walking on the level, may have inhibited our ability to determine DBA-energy expenditure relationships for each angle or indeed determine an overall relationship for V̇O2 which incorporates VeDBA and terrain slope.
The variation observed in the V̇O2–VeDBA relationship defined for “incline” walking does not provide confident estimates of V̇O2 with slope. Importantly, the use of VeDBA as a proxy for energy expenditure based upon a single relationship between V̇O2 and VeDBA, using Eq. 1 can result in large variation in predictions of the associated energy expenditure, ranging in this case from underestimates of 22.4 ± 12.7% for the 4° incline, 46.7 ± 2.3% for a 6° incline, and 16.2 ± 7.11% for a 6° decline when using VeDBA measurements on these incline. Where the decline angle is 4°, VeDBA could overestimate the cost of locomotion, depending on the speed of movement exhibited (6.8 ± 9% error, Eq. 1).
Terrier, Aminian and Schutz  suggested that a solution for poor predictive capacity of power (watts) from acceleration metrics in the case of incline movement was to assess slopes by independent methods (such as altimeters or differential GPS [42, 43]) so that slope-specific acceleration proxies could be applied. Certainly, the rapid development of sensors within animal-attached loggers such as altimeters (e.g., ) or logger pitch angle is making this increasingly feasible.
Factors affecting V̇O2
The positive relationship observed in the current study between V̇O2 and walking speed is in agreement with Clapperton,  who examined V̇O2 of sheep walking at similar inclines (Fig. 4). Dailey and Hobbs  examined the V̇O2 of sheep and goats walking up much steeper inclines (25°). Their V̇O2 values were higher, which is expected as increased energy is required to move their mass against gravity . Interestingly, in the current study which incorporated decline walking, once slope angle decreased to less than − 4°, energy expenditure increased (Fig. 2C). A similar trend was found by Minetti et al.  in humans and in a number of animals by Halsey and White . This suggests that steep downhill walking incurs greater costs than walking on a level surface, which may be due, for example, to the energy required for controlling velocity and for gait stabilization . In support of this, the use of a quadratic rather than a linear function to examine the relationship between V̇O2 and treadmill angle supports the notion that energy expenditure may increase with increasing downhill gradient (Fig. 2D). Nevertheless, the efficacy of using either a linear or quadratic relationship should be interpreted with caution for n = 2 subjects. We only measured two animals was for the following two reasons: The first was logistical. It took 6 weeks to tame and train sheep from lambs to adults that would walk on the treadmill at will. It also took a further 11 months to collect data during which only one measurement could be collected per animal per day. Second, we were interested in exploring the concept of DBA–V̇O2 relationship variation with terrain slope and walking speed, rather than determining species-specific values per se. Therefore, while inclusion of data from additional subjects would provide more information yielding more accurate energy expenditure values for sheep of that breed in general, it would probably not change the qualitative observation that V̇O2 and DBA vary with terrain slope and speed.
The lower critical temperature (LCT) of sheep with a full fleece in the current study was determined as 18 °C. This is similar to the 15–20 °C previously described [46, 47]. Above the LCT, resting V̇O2 of non-shorn sheep was determined at 3.67 mL kg−1 min−1. While the current study did not seek to determine the upper limits of the thermoneutral zone, previous studies have suggested that the upper critical temperature of sheep lies between 25 and 40 °C [48, 49]. The highest temperature to which the sheep were exposed in the current study was 32 °C. In terms of the V̇O2 of shorn sheep, although we did not identify the lower critical temperature, work by Blaxter et al.  reported a value of c. 23 °C. Here the V̇O2 values of shorn sheep were not different to the V̇O2 values of the sheep with a full fleece at temperatures above 18 °C but V̇O2 was higher at lower temperatures suggesting that there are indeed thermoregulatory costs associated with fleece removal.
Animal pitch as a proxy for terrain angle
Although a positive relationship was apparent between measurements of pitch and angle of the treadmill (equivalent to terrain slope), large variations in animal pitch (28%) occurred, predominantly during decline locomotion. Thus, the pitch of an animal’s body may not necessarily be an accurate representation of the terrain slope. Differences in terrain slope and animal pitch may be due to postural changes by the animal attempting to stabilise and control downhill movement. Presumably this is akin to a human maintaining an upright posture when walking up or downhill . Biewener  reported that animals may adopt a more upright posture of the limbs to reduce the stresses of downhill walking on muscles and joints, which may lead to a reduced downward body angle. As the current data were obtained across multiple locomotion speeds, it is possible that the variation in animal pitch may also be due to such changes in gait . Thus, although a measure of terrain slope can be estimated using animal pitch angle, we also suggest that additional methods, such as altimeters or GPS data , can provide additional information on whether an animal travels up or down an incline and subsequently be used to inform an estimate of energy expenditure in the field (e.g., ).
While current results describe associations between DBA and V̇O2, we understand that these are not representative of the species. However, results highlight significant variation between individuals and under different conditions (speed, terrain slope). Hence, we demonstrate that DBA–V̇O2 relationships depend on a variety of environmental conditions. Field studies using DBA as a proxy for energy expenditure should also consider other aspects of the terrain, such as the penetrability of the substrate on which the animal is moving, and terrain evenness (i.e., regions, where changes in terrain angle are frequent and localised, to the extent that perhaps the forelimbs might even be on ground that is at a different angle than the ground on which the hind limbs stand) and explore how these may influence the relationship between acceleration metrics and energy expenditure.
We recorded simultaneous measurements of oxygen consumption (V̇O2) by indirect calorimetry and body movement by tri-axial accelerometry to determine whether the metric ‘vectorial dynamic body acceleration’ (VeDBA) could be used as a proxy for instantaneous energy expenditure. Results showed that VeDBA increased linearly with V̇O2 as walking speed increased, and that this relationship varied with both incline and level terrain slope. Therefore, we recommend that future attempts to determine V̇O2–VeDBA relationships should examine a variety of incline and decline terrain slopes according to individual circumstance. We also show that although animal pitch could be used as a broad measure of terrain slope, this relationship was variable, especially for locomotion on a decline, presumably because subjects make adjustments in their posture.
Animals and study site
The experiments were conducted on a working farm near the town of Carnlough, Co. Antrim, Northern Ireland. The apparatus was set up inside a barn (c. 5 × 10 m floor space, 4 m height), which was exposed to natural lighting and ambient temperature with one whole side of the building remaining open. Two non-reproductive adult (1–2 years) Border-Cheviot ewes (66.47 ± 2.32 kg, sheep A and sheep B), which had been habituated to human contact since birth, were used. Prior to measurements of V̇O2 (see below), sheep initially underwent a 6-week training period which included becoming accustomed to a metabolic chamber and being exercised on a treadmill for 15 min at a time at 0.8, 1.5, 2.2 and 2.9 km h−1 twice per day, three times a week. During this period, V̇O2 measurements were not taken. Sheep were provided with c. 5 g of concentrate pellets (Thompsons Feed Innovation, Belfast) as a food reward when they entered the chamber and walked on the treadmill. Between training sessions, sheep were held in a field (0.12 ha) adjacent to the barn, where grass and water were available ad libitum.
Indirect calorimetry setup
Measurements of V̇O2 were determined using indirect calorimetry. A metabolic chamber (100 × 130 × 50 cm) made from clear Perspex® sheeting was mounted upon a treadmill (running area: 33 × 114 cm; Vfit, DOG Jogger™ Canine Rehabilitation Treadmill, West Yorkshire, UK). Ambient air was drawn from outside through plastic tubing (3 cm diameter, 100 cm length), to a pump (3.15 CFM 120 VAC, THOMAS, Munich, Germany). The air was then pushed from the pump through plastic tubing (1.5 cm diameter, 70 cm length), through a flowmeter (± 1.25% FSD, CT Platon NG series Flow Meter 10–100 L min−1 N2, Hampshire, England) and then through a copper coil (1.5 cm diameter, 150 cm length) which was placed inside a cabinet (PTC-1 Cabinet, Sable Systems Europe©, Berlin, Germany) controlled by a Pelt-5 temperature controller (Sable Systems Europe©, Berlin, Germany) which either warmed or cooled the air flowing through the coil to the desired temperature before entering the chamber at positive pressure (Fig. 5A). Within the metabolic chamber, air was mixed by five inbuilt electric fans. Temperature and pressure were measured inside the chamber using a digital thermometer and barometer (resolution ± 0.1°/± 0.1 hpa, Sunroad FR500, EverTrust™ UK). A subsample of air was drawn from inside the chamber at 300 mL min−1 and passed through a drying column (Drierite, Fisher Scientific, Leicestershire, England) before entering a gas analyser (Foxbox Field Gas Analysis System, Sable Systems Europe©, Berlin, Germany) which measured O2 concentration (O2%). Prior to experiments, the gas analyser was factory-calibrated by Sable Systems Europe© (Berlin, Germany).
The time taken for gases to reach equilibration within the system was determined by a series of Nitrogen (N2) gas injections [13, 53], which were undertaken prior to animal measurements. For these experiments, N2 was injected at various rates (L min−1) into the chamber and the time for O2% within the chamber to reach equilibration was measured at flow rates of air entering the chamber ranging from 20 to 90 L min−1 (Fig. 5B).
Measurements of resting V̇O2 and variation in V̇O2 with ambient temperature
Experiments were conducted between October 2016 and August 2017 during this time ambient temperatures ranged between − 3 and 24 °C. Sheep were weighed (large animal scales, Tree LC–VS 180, ± 0.05 kg, LW measurements LLC©, California) prior to measurements of oxygen concentration, which were determined between 11:00 and 16:00, during the normal period of activity for these animals . Before an individual was placed inside the chamber, the chamber conditions (temperature, air flow) were allowed to stabilise by running the pumps and the analyser for 60 min. Thereafter, the subject animal was encouraged to enter the chamber via a sliding door at the rear (Fig. 5A). This was then closed, secured, and the time the animal entered the chamber was recorded. Measurements of O2%, ambient pressure and temperature commenced as soon as the animal was observed to be stationary inside the chamber and were then taken every 5 min until the concentration of oxygen within the chamber was deemed to have ‘settled’. This was determined as a change of less than 0.1% O2 over a 5-min period . This process took approximately 20 min. Once stable conditions had been reached, the O2% within the chamber, temperature and pressure were recorded every 60 s for a further 5 min to enable 5 min of stable recordings to be taken with the animal at rest. If the animal became agitated at any time during the proceedings, or did not rest, measurements were abandoned and attempted on another day. Once measurements were complete, the animal was promptly removed from the respirometry chamber and the door of the chamber closed and secured. Measurements of O2% within the chamber were once again measured every 5 min until the O2% inside the chamber increased to a constant reading (e.g., to an ambient concentration of 20.95%), which took 15–20 min. At this point, measurements increased in frequency to one every 60 s. When O2% changed < 0.1% over a further 5-min period, the system was deemed to have stabilised. An average of these five measurements was then taken to determine the drift in O2% from the original span of 20.95%. This protocol was repeated for a range of temperatures between 8 and 32 °C for both subjects. One measurement was collected per temperature per animal per day. Towards the end of the study (July–August 2017), sheep were shorn (normal practice on the working farm) and V̇O2 measurements over a range of ambient temperatures were repeated. Note that we did not conduct experiments in which the subjects were walking on shorn sheep. As these latter measurements were conducted during the summer, the higher ambient temperatures (21–24 °C) precluded lower temperatures being achieved within the chamber.
Calculation of V̇O2
To correct for drift in the oxygen analyser, the final O2% was subtracted from the initial span of 20.95% to determine the difference in oxygen percentages. This was then divided by the time difference between the initial span and final recorded O2% to determine drift per minute, assuming a linear drift over time . The corrected O2% for each reading was determined as follows:
where O2% is the percentage O2 recorded, t is time between initial span and the reading obtained.
Rate of oxygen consumption (V̇O2)
To determine V̇O2, first flow rate of air through the chamber was corrected to standard volume (V̇std), (standard temperature (273.15 K) and pressure (760 mmHg) as in Winberry, ), as follows:
where V̇i is the volume of gas (L min−1) sampled at the measured atmospheric pressure (Pa) and temperature (Ta). Tstd (K) and Pstd (mmHg) are standard temperature (273.15 K (0 °C, 32F)) and standard pressure (100 kPa, 1 bar), respectively . V̇std was then converted to mL min−1.
V̇O2 (mL min−1) was then calculated following Withers  as
where Fi is the fractional percentage of oxygen flowing into the chamber and Fe is the fractional percentage of oxygen exiting the chamber. V̇O2 was then divided by animal body mass to return mass-specific rate of oxygen consumption (V̇O2; mL kg−1 min−1).
Tri-axial accelerometer loggers
Daily diary ‘DD’ loggers (Wildbyte Technologies Ltd, Swansea, UK) were secured to the sheep via a harness (detailed below) prior to them entering the metabolic chamber. These devices comprised a 12-bit resolution, tri-axial accelerometer (± 16 g), which recorded acceleration on three orthogonal axes (‘surge’—forwards and backwards; ‘heave’—up and down and ‘sway’—side to side, ) at 40 Hz. Devices were housed inside vacuum-formed plastic cases (4.5 × 1.8 × 4 cm) and were powered using a 3.7 V lithium-ion rechargeable battery. Data recorded by the device were saved to an onboard 32 GB MicroSD card, (SanDisk®, Western Digital Technologies, Inc). The device and battery were wrapped in waterproof, 19 mm PVC insulation tape to prevent water ingress.
Attachment of accelerometer loggers to sheep
Accelerometers were secured to a nylon ram harness (Nettex©; Kent, England) via Tesa®; tape (No. 4651; Tesa AG, Hamburg, Germany) and bound by four cable ties (300 × 4.5 mm) to minimise rotation or movement of the device not caused by animal activity. The harness was then fitted to the subject (Fig. 1A). Sheep remained stationary while being instrumented with the equipped harness and did not require restraint. The total weight of the device, batteries and harness was 460 g. Attachment of the logger and harness took less than 5 min per animal. Devices were positioned on the dorsal mid-line of each animal, behind the shoulders. Care was taken to ensure each accelerometer was orientated with the Y axis representing forward and backward movement, the Z axis representing up and down movement and the X axis representing side to side movement (Fig. 1A).
Measurements were conducted on one animal at a time. The temperature inside the chamber for these measurements was above the lower critical point for each animal. During this experiment, each sheep was initially placed on a level treadmill and subjected to four different speeds: 0.8, 1.5, 2.2 and 2.9 km h−1. As above, only one speed was measured on any one occasion, i.e., each session on the treadmill resulted in one data point. Speeds were chosen based on previous data of sheep walking with no changes in gait [27, 59]. Thereafter, experiments were conducted with the treadmill set at inclines of 4° and 6°, and declines of 4° and 6°, at two speeds, which were 0.8 and 2.2 km h−1. These gradients, at 4° and 6°, are equivalent to slopes of 7% or 1 in 14.3, and 10.5% or 1 in 9.5, respectively (Fig. 5A) and are the equivalents of walking on a gently sloping field. Once a sheep was inside the chamber, the treadmill was switched on and set to the desired speed. As per the description of the measurements of resting V̇O2 above, measurement of O2%, ambient pressure and chamber temperature were conducted every 5 min until the O2% within the chamber changed less than 0.1% over 5 min. At this point, the O2%, temperature and pressure were recorded every 60 s for a further 5 min. After this, the treadmill was safely stopped, and the animal was removed from the chamber. Again, as per the resting V̇O2 determination described above, measurements of O2% within the chamber were measured every 5 min until the O2% inside the chamber increased to a constant reading (ambient O2% of c. 20.95%). After this, measurements increased to every 60 s. When O2% changed < 0.1% over a further 5-min period, the system was deemed to have stabilised. An average of these five measurements was then taken as described above to determine the drift in O2% from the original span.
Calculation of vectorial dynamic body acceleration (VeDBA)
Acceleration data were measured during the 5-min steady-state period when the sheep was inside the treadmill, from which the measure vectorial dynamic body acceleration (VeDBA) values were derived. We used VeDBA as a proxy for activity-related energy expenditure, rather than other measures, such as overall dynamic body acceleration (ODBA) . The effectiveness of different measures of DBA has been debated previously (e.g., ). ODBA and VeDBA, are highly correlated, and both provide reliable estimates of V̇O2. However, we used VeDBA, because it has been suggested to be the preferred metric if device orientation cannot be maintained or be consistent , and therefore, the results should be applicable to a wide range of studies. For the calculation of VeDBA, the element of measured acceleration that is static acceleration was estimated by taking a running mean with a 2-s window of the raw acceleration data for each axis . Dynamic acceleration was then calculated by subtracting the static acceleration from the raw acceleration for each axis . VeDBA was then calculated according to Qasem et al.  as
where Ax is dynamic acceleration on the X axis (‘sway’), and Ay and Az are the dynamic accelerations of the Y (‘surge’) and Z (‘heave’) axes, respectively.
Calculation of animal pitch
The pitch of the animal, i.e., the angle of the tag placed on dorsal mid-line, behind the shoulders relative to the horizontal plane, \(\theta\) (Fig. 1a) was calculated according to Collins et al. :
where statAy is static acceleration recorded on the Y (‘surge’) axis and statAx and statAz are static acceleration recorded on the X (‘sway’) and Z (‘heave’) axes, respectively.
Analyses were conducted using R version 3.4.4 . To define the lower critical point, when sheep where stationary and exposed to different ambient temperatures, broken-stick linear regression models were fitted using the “segmented” package . This enabled inflection points for the relationship between V̇O2 and chamber temperature to be identified  for each sheep. Differences in V̇O2 before and after shearing were assessed with a general linear mixed model (GLMM) that included individual identity (i.e., sheep A and sheep B) as a random factor, using the “lme4” package . This package was also used to conduct another GLMM (Table 1), with log transformed V̇O2 as the dependent variable and VeDBA as the main effect and identity as a random factor. Further separate GLMMs with log transformed V̇O2 and VeDBA as the dependent variables were conducted to determine how V̇O2 and VeDBA varied depending on interactions between speed and pitch, and speed and slope. Speed, pitch, and slope were included in these two models as main effects, with individual identity as a random factor. Stepwise deletion refined each model, and the best model was selected based on the Akaike information criterion (AIC) (Table 1). Residuals and fitted values of the final models were checked via histograms and scatterplots to ensure model assumptions were met. V̇O2 was log-transformed to generate normal residuals for the models. The relationship between V̇O2 and treadmill angle was also explored using a quadratic relationship to examine whether there was any evidence for V̇O2 increasing at steeper gradients , both for inclines and declines. Due to limited data for each slope, Spearman’s rank correlations were used to assess relationships between V̇O2 and VeDBA at each individual treadmill angle using data obtained from walking at all speeds. Results were deemed to be significant if P < 0.05. Figures were created using the “ggplot2” package .
Availability of data and materials
The data sets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Hole CV, Goyens J, Prims S, Fransen E, Hernando MA, Van Cruchten S, et al. How innate is locomotion in precocial animals? A study on the early development of spatio-temporal gait variables and gait symmetry in piglets. J Exp Biol. 2017;220:2706–16.
Brosh A, Henkin Z, Ungar ED, Dolev A, Shabtay A, Orlov A, et al. Energy cost of activities and locomotion of grazing cows: a repeated study in larger plots. J Anim Sci. 2010;88:315–23.
Scantlebury DM, Mills MG, Wilson RP, Wilson JW, Mills ME, Durant SM, et al. Flexible energetics of cheetah hunting strategies provide resistance against kleptoparasitism. Science. 2014;346:79–81.
Bertram A. Understanding mammalian locomotion: concepts and applications. Hoboken: Wiley; 2016. p. 386.
Grémillet D, Lescroël A, Ballard G, Dugger KM, Massaro M, Porzig EL, Ainley DG. Energetic fitness: field metabolic rates assessed via 3D accelerometry complement conventional fitness metrics. Funct Ecol. 2018;32:1203–13.
Pagano AM, Atwood TC, Durner GM, Williams TM. The seasonal energetic landscape of an apex marine carnivore, the polar bear. Ecology. 2020;101:3.
Williams TM, Peter-Heide Jørgensen M, Pagano AM, Bryce CM. Hunters versus hunted: new perspectives on the energetic costs of survival at the top of the food chain. Funct Ecol. 2020;34:2015–29.
Cavagna GA, Saibene FP, Margaria R. External work in walking. J Appl Phys. 1963;18:1–9.
Meijer GA, Westerterp KR, Koper HANS. Assessment of energy expenditure by recording heart rate and body acceleration. Med Sci Sports Exerc. 1989;21:343–7.
Terrier P, Aminian K, Schutz Y. Can accelerometry accurately predict the energy cost of uphill/downhill walking? Ergonomics. 2001;44:48–62.
Wilson RP, White CR, Quintana F, Halsey LG, Liebsch N, Martin GR, Butler PJ. Moving towards acceleration for estimates of activity-specific metabolic rate in free-living animals: the case of the cormorant. J Anim Ecol. 2006;75:1081–90.
Halsey LG, Shepard ELC, Hulston CJ, Venables MC, White CR, Jeukendrup AE, Wilson RP. Acceleration versus heart rate for estimating energy expenditure and speed during locomotion in animals: tests with an easy model species, Homo sapiens. Zoology. 2008;111:231–41.
Halsey LG, Shepard ELC, Quintana F, Laich AG, Green JA, Wilson RP. The relationship between oxygen consumption and body acceleration in a range of species. Comp Biochem Phys Part A Mol Int Phys. 2009;152:197–202.
Gleiss AC, Dale JJ, Holland KN, Wilson RP. Accelerating estimates of activity-specific metabolic rate in fishes: testing the applicability of acceleration data-loggers. J Exp Mar Biol Ecol. 2010;385:85–91.
Halsey LG, White CR. Measuring energetics and behaviour using accelerometry in cane toads Bufo marinus. PLoS ONE. 2010;5:4.
Laich AG, Wilson RP, Gleiss AC, Shepard EL, Quintana F. Use of overall dynamic body acceleration for estimating energy expenditure in cormorants: does locomotion in different media affect relationships? J Exp Mar Biol Ecol. 2011;399:151–5.
Qasem L, Cardew A, Wilson A, Griffiths I, Halsey LG, Shepard EL, et al. Tri-axial dynamic acceleration as a proxy for animal energy expenditure; should we be summing values or calculating the vector? PLoS ONE. 2012;7:2.
Pagano AM, Carnahan AM, Robbins CT, Owen MA, Batson T, Wagner N, et al. Energetic costs of locomotion in bears: is plantigrade locomotion energetically economical? J Exp Biol. 2018;221:12.
Dunford CE, Marks NJ, Wilmers CC, Bryce CM, Nickel B, Wolfe LL, Scantlebury DM, Williams TM. Surviving in steep terrain: a lab-to-field assessment of locomotor costs for wild mountain lions (Puma concolor). Mov Ecol. 2020;8:1–12.
Halsey LG, Bryce CM. Proxy problems: why a calibration is essential for interpreting quantified changes in energy expenditure from biologging data. Funct Ecol. 2021;35:627–34.
di Virgilio A, Morales JM, Lambertucci SA, Shepard EL, Wilson RP. Multi-dimensional Precision Livestock Farming: a potential toolbox for sustainable rangeland management. PeerJ. 2018;6:e4867.
Miwa M, Oishi K, Nakagawa Y, Maeno H, Anzai H, Kumagai H, et al. Application of overall dynamic body acceleration as a proxy for estimating the energy expenditure of grazing farm animals: relationship with heart rate. PLoS ONE. 2015;10:6.
Green JA, Halsey LG, Wilson RP, Frappell PB. Estimating energy expenditure of animals using the accelerometry technique: activity, inactivity and comparison with the heart-rate technique. J Exp Biol. 2009;212:471–82.
Bidder OR, Goulding C, Toledo A, van Walsum TA, Siebert U, Halsey LG. Does the treadmill support valid energetics estimates of field locomotion? Int Comp Biol. 2017;57:301–19.
White RG, Yousef MK. Energy expenditure in reindeer walking on roads and on tundra. Can J Zool. 1978;56:215–23.
Davies SEH, Mackinnon SN. The energetics of walking on sand and grass at various speeds. Ergonomics. 2006;49:651–60.
Clapperton JL. The energy metabolism of sheep walking on the level and on gradients. Br J Nutr. 1964;18:47–54.
Brockway JM, Gessaman JA. The energy cost of locomotion on the level and on gradients for the red deer (Cervus elaphus). Q J Exp Physiol Cogn Med Sci Trans Int. 1977;62:333–9.
Dailey TV, Hobbs NT. Travel in alpine terrain: energy expenditures for locomotion by mountain goats and bighorn sheep. Can J Zool. 1989;67:2368–75.
Halsey LG, White CR. A different angle: comparative analyses of whole-animal transport costs when running uphill. J Exp Biol. 2017;220:161–6.
Carnahan AM, van Manen FT, Haroldson MA, Stenhouse GB, Robbins CT. Quantifying energetic costs and defining energy landscapes experienced by grizzly bears. J Exp Biol. 2021;224:6.
Wickler SJ, Hoyt DF, Biewener AA, Cogger EA, Kristin L. In vivo muscle function vs speed IIl Muscle function trotting up an incline. J Exp Biol. 2005;208:1191–200.
Minetti AE, Moia C, Roi GS, Susta D, Ferretti G. Energy cost of walking and running at extreme uphill and downhill slopes. J Appl Phys. 2002;93:1039–46.
Lambert MG, Clark DA, Grant DA, Costall DA, Fletcher RH. Influence of fertiliser and grazing management on North Island moist hill country: 1. Herbage accumulation. N Z J Agric Res. 1983;26:95–108.
López IF, Hodgson J, Hedderley DI, Valentine I, Lambert MG. Selective defoliation by sheep according to slope and plant species in the hill country of New Zealand. Grass Forage Sci. 2003;58:339–49.
Shepard EL, Wilson RP, Quintana F, Laich AG, Liebsch N, et al. Identification of animal movement patterns using tri-axial accelerometry. Endanger Species Res. 2008;10:47–60.
Department of Environment, Food and Rural Affairs (DEFRA), Farming Statistics—final crop areas, yields, livestock populations and agricultural workforce at 1 June 2020 United Kingdom (data to 1st June 2020). Last updated 22 December 2020. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/946161/structure-jun2020final-uk-22dec20.pdf. Accessed 26 Sept 2021.
Department of Environment, Food and Rural Affairs (DEFRA), Total Income from Farming in the United Kingdom, first estimate for 2020. Last updated 27 May 2021. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/989701/agricaccounts-tiffstatsnotice-27may21.pdf. Accessed 26 Sept 2021.
Australian Bureau of Agricultural and Resource Economics and Sciences (ABARES), Agricultural Commodities. Australia. September 2019. https://www.awe.gov.au/sites/default/files/sitecollectiondocuments/abares/agriculture-commodities/AgCommodities201909_SheepOutlook_v1.0.0.pdf. Accessed 30 Dec 2021.
IBISWorld statistics, Sheep and Cattle Farming Industry in China—Market Research Report. Last updated 4 January 2022. https://www.ibisworld.com/china/market-research-reports/sheep-cattle-farming-industry/#:~:text=The%20Sheep%20and%20Cattle%20Farming,through%202020%2C%20to%20%24106.7%20billion. Accessed 21 Jan 2022.
Blaxter KL, Graham NM, Wainman FW. Environmental temperature, energy metabolism and heat regulation in sheep III. The metabolism and thermal exchanges of sheep with fleeces. J Agric Sci. 1959;52:41–9.
Perrin O, Terrier P, Ladetto Q, Merminod B, Schutz Y. Improvement of walking speed prediction by accelerometry and altimetry, validated by satellite positioning. Med Biol Eng Comput. 2000;38:164–8.
Duncan GE, Lester J, Migotsky S, Higgins L, Borriello G. Measuring slope to improve energy expenditure estimates during field-based activities. Appl Phys Nutr Metab. 2012;38:352–6.
Shepard EL, Lambertucci SA, Vallmitjana D, Wilson RP. Energy beyond food: foraging theory informs time spent in thermals by a large soaring bird. PLoS ONE. 2011;6:11.
Hunter LC, Hendrix EC, Dean JC. The cost of walking downhill: is the preferred gait energetically optimal? J Biomech. 2010;43:1910–5.
Blaxter LL. The energy metabolism of ruminants. London: Hutchinson and Co.; 1962.
Khalifa HH, Shalaby T, Abdel-Khalek TMM. An approach to develop a biometeorological thermal discomfort index for sheep and goats under Egyptian conditions. In: Proceeding of the 17th international congress of biometeorology. Int Soc Biomet. 2005;118–122.
Alexander G. Heat loss from sheep. In: Heat loss from animals and man: assessment and control. 1974; 173–203.
Randall JM. Environmental parameters necessary to define comfort for pigs, cattle and sheep in livestock transporters. Anim Sci. 1993;57:299–307.
Leroux A, Fung J, Barbeau H. Postural adaptation to walking on inclined surfaces: I. Normal strategies. Gait Posture. 2002;15:64–74.
Biewener AA. Scaling body support in mammals: limb posture and muscle mechanics. Science. 1989;245:45–8.
Fukuoka Y, Kimura H, Cohen AH. Adaptive dynamic walking of a quadruped robot on irregular terrain based on biological concepts. Int J Robot Res. 2003;22:187–202.
Fedak MA, Rome L, Seeherman HJ. One-step N2-dilution technique for calibrating open-circuit V̇O2 measuring systems. J Appl Phys. 1981;51:772–6.
Fogarty ES, Manning JK, Trotter MG, Schneider DA, Thomson PC, et al. GNSS technology and its application for improved reproductive management in extensive sheep systems. Anim Prod Sci. 2015;55:1272–80.
McClune DW, Kostka B, Delahay RJ, Montgomery WI, Marks NJ, Scantlebury DM. Winter is coming: seasonal variation in resting metabolic rate of the European badger (Meles meles). PLoS ONE. 2015;10:9.
Lighton JRB. Chapter 3: Metabolic measurement techniques: baselining, mathematical correction of water vapour dilution and response correction. In: Gerrits W, Labussière E, editors. Indirect calorimetry: techniques, computations and applications. Wageningen: Wageningen Academic Publishers; 2015. p. 2993–3000.
Winberry W. Compendium of methods for the determination of inorganic compounds in Ambient Air:(Chapter IO-2.4) Calculations for Standard Volume, US Environmental Protection Agency EPA. 1999.
Withers PC. Design, calibration and calculation for flow-through respirometry systems. Aust J Zool. 2001;49:445–61.
Squires VR, Wilson AD, Daws GT. Comparisons of the walking activity of some Australian sheep. Proc Aust Soc Anim Prod. 1972;9:376–80.
Collins PM, Green JA, Warwick-Evans V, Dodd S, Shaw PJ, Arnould JP, Halsey LG. Interpreting behaviors from accelerometry: a method combining simplicity and objectivity. Ecol Evol. 2015;5:4642–54.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. 2018. https://www.R-project.org/.
Muggeo VM. Segmented: an R Package to fit regression models with broken-line relationships. R News. 2008;8:20–5.
McKechnie AE, Smit B, Whitfield MC, Noakes MJ, Talbot WA, Garcia M, et al. Avian thermoregulation in the heat: evaporative cooling capacity in an archetypal desert specialist, Burchell’s sandgrouse (Pterocles burchelli). J Exp Biol. 2016;219:2137–44.
Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, et al. Package ‘lme4.’ Convergence. 2015;12:1–113.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.
We would like to thank the farmer, Mr. H. McKay, for his assistance allowing access to his farm and animals for this study. Thank you to Dr. Natasha McGowan for her advice and assistance. We also thank Nick Decalmer, sadly lost to us before his time, who built the bespoke respirometer chamber used in this study.
We are grateful to the Department of Agriculture, Environment and Rural Affairs (DAERA), Northern Ireland, for providing a studentship to CCM, through which this research was conducted.
Ethics approval and consent to participate
Ethical approval was granted by the School of Biological Sciences Ethical Review Committee, Queen’s University Belfast and the UK Home Office. The owner of the animals, Mr H. McKay, consented to the research.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mulvenna, C.C., Marks, N.J., Wilson, R.P. et al. Can metrics of acceleration provide accurate estimates of energy costs of locomotion on uneven terrain? Using domestic sheep (Ovis aries) as an example. Anim Biotelemetry 10, 8 (2022). https://doi.org/10.1186/s40317-022-00281-3
- Dynamic body acceleration
- Energy expenditure
- Oxygen consumption
- Indirect calorimetry