Using open-flow respirometry and video footage (25 frames s–1), the energy expenditure and hindlimb kinematics of barnacle geese, Branta leucopsis, were measured whilst they were exercising on a treadmill at gradients of +7 and –7 deg, and on a level surface. In agreement with previous studies, ascending a gradient incurred metabolic costs higher than those experienced on level ground at comparable speeds. The geese, however, are the first species to show an increased duty factor when ascending a gradient. This increased duty factor was accompanied by a longer stance time, which was probably to enable the additional force required for ascending to be generated. Contrary to previous findings, the geese did not experience decreased metabolic costs when descending a gradient. For a given speed, the geese took relatively shorter and quicker strides when walking downhill. This ‘choppy’ stride and perhaps a lack of postural plasticity (an inability to adopt a more crouched posture) may negate any energy savings gained from gravity's assistance in moving the centre of mass downhill. Also contrary to previous studies, the incremental increase in metabolic cost with increasing speed was similar for each gradient, indicating that the efficiency of locomotion (mechanical work done/chemical energy consumed) is not constant across all walking speeds. The data here suggest that there are species-specific metabolic responses to locomotion on slopes, as well as the established kinematics differences. It is likely that a suite of factors, such as ecology, posture, gait, leggedness and foot morphology, will subtly affect an organism's ability to negotiate gradients.

Most terrestrial animals are required, at some point, to locomote up or down gradients. Ascending an incline should be more mechanically costly than locomotion on level ground, because of the requirement for additional work to increase the potential energy of the centre of mass (CoM) as it ascends the slope (Parsons et al., 2008). To produce this extra mechanical work, a greater volume of limb muscle may be recruited (Carlson-Kuhta et al., 1998; de Leon et al., 1994; Gardiner et al., 1982; Leroux et al., 1999; Pierotti et al., 1989; Robert et al., 2000; Roy et al., 1991) or muscle-shortening strains increased (Gillis and Biewener, 2002). Additional muscle activity means that whole-animal metabolic costs (rate of O2 consumption, O2) should also increase (Ellerby et al., 2003), assuming that additional costs are met aerobically. Accordingly, increased metabolic costs were found in birds (Bamford and Maloiy, 1980; Ellerby et al., 2003; Rubenson et al., 2006; Snyder and Carello, 2008; Warncke et al., 1988), humans (Dill, 1965; Margaria et al., 1963; Minetti et al., 2002; Yousef et al., 1972), a range of quadrupeds (Armstrong et al., 1983; Bedford et al., 1979; Brooks and White, 1978; Chassin et al., 1976; Clapperton, 1964; Cohen et al., 1978; Eaton et al., 1995; Farley and Emshwiller, 1996; McDonough et al., 2002; Raab et al., 1976; Robbins et al., 1979; Snyder and Carello, 2008; Taylor et al., 1972; White and Yousef, 1978; Wickler et al., 2000; Wickler et al., 2003; Wunder and Morrison, 1974; Yousef et al., 1972) and invertebrates (Full and Tullis, 1990; Tullis and Andrus, 2011) when ascending gradients. Although there are notable exceptions – the marabou stork Leptoptilos crumeniferus (Bamford and Maloiy, 1980), cockroach Gromphadorhina portento (Herreid et al., 1981) and ant Camponotus sp. (Lipp et al., 2005) – generally, for most species, as gradient increases so does O2 at all locomotory speeds (Bedford et al., 1979; Brooks and White, 1978; Chassin et al., 1976; Clapperton, 1964; Cohen et al., 1978; Eaton et al., 1995; Ellerby et al., 2003; Full and Tullis, 1990; Raab et al., 1976; Robbins et al., 1979; Wunder and Morrison, 1974; Yousef et al., 1972). The relative cost of ascending or descending a gradient may also be size dependent, with a lower percentage increase or decrease (compared with locomotion on the level) of locomoting uphill or downhill, respectively, for smaller animals (Taylor et al., 1972).

Based upon similar reasoning, locomoting down a decline is predicted to be less metabolically costly than locomotion on the level, because here the CoM is moving ‘with’ and not ‘against’ gravity. Descending a gradient is not metabolically free, however, because the muscles must absorb mechanical energy to decrease the potential energy of the body. The absorption of potential energy is facilitated by leg muscles, which appear to shift function in response to gradient. In both bipeds (Roberts et al., 2007) and quadrupeds (Lammers et al., 2006) the leg muscles appear to generate additional energy (an increase in positive work) to propel the CoM uphill (against gravity), whereas on a decline they shift function to energy absorption (negative work) and operate as braking struts (Lammers et al., 2006; Roberts et al., 2007; Smith et al., 1998). Opposite to the situation during incline locomotion, less limb muscle is recruited relative to that during level locomotion when travelling downhill (Kai et al., 1997; Smith and Carlson-Kuhta, 1995).

Most empirical studies have recorded reduced metabolic costs for animals descending declines (Armstrong et al., 1983; Byrnes et al., 1985; Cohen et al., 1978; Margaria et al., 1963; Minetti et al., 2002; Raab et al., 1976; Robbins et al., 1979; Taylor et al., 1972; White and Yousef, 1978; Yousef et al., 1972). Ants are an exception, but this may be due to the resolution of the respirometry equipment not allowing detection of the small differences expected in very small animals (Lipp et al., 2005). Although there are exceptions (Yousef et al., 1972), O2 also generally decreases with a steepening decline (Cohen et al., 1978). The saving in metabolic costs when descending relative to locomotion on a horizontal surface tends not to be as great in absolute terms as the increased costs associated with ascending a slope of a similar angle (Armstrong et al., 1983; Cohen et al., 1978; Raab et al., 1976; Robbins et al., 1979; Taylor et al., 1972; White and Yousef, 1978). This implies that the whole-animal mechanical efficiency (mechanical work done/chemical energy consumed) is lower on declines.

Limb and body posture, and kinematics may also change in response to gradients (Carlson-Kuhta et al., 1998; Jayne and Irschick, 1999; Leroux et al., 1999; Smith et al., 1998). Although there are some broad trends in the responses of animals to locomoting on gradients, the findings of previous studies suggest that there are subtle species-specific differences in the kinematic responses (Gillis and Biewener, 2002). Hence, more comparative data are required before the constraints acting upon the locomotor system when negotiating slopes can be fully understood. Most research examining terrestrial locomotion with respect to gradients, including bird studies, has focused on species whose movements are primarily facilitated through terrestrial locomotion. These species are presumably highly adapted for efficient walking or running. Although capable of terrestrial locomotion, the barnacle goose has conflicting demands for flight and swimming, which compromise its morphology in terms of terrestrial locomotion. One obvious consequence of these competing demands is the barnacle goose's waddling walking gait. During their annual moult, barnacle geese are grounded and during the rearing period they spend much of their time walking with their goslings in search of food. Hence, effective terrestrial locomotion is still fundamental to a barnacle goose's survival, and as a generalist species with competing locomotor demands they may give different insight into the constraints operating on the vertebrate locomotor system than species more focused on cursoriality (Nudds et al., 2010).

In this study, barnacle geese had their metabolic rates measured whilst exercising at a range of speeds upon a treadmill on 7 deg positive and negative gradients, and on the level. It was predicted that, because barnacle geese are not specifically adapted to a terrestrial locomotory lifestyle, the metabolic cost of their uphill and downhill locomoting would be higher and not as low, respectively, relative to their level walking when compared with specialist terrestrial species.

Study species

Five barnacle geese Branta leucopsis (Bechstein 1803) as used in a previous study (Nudds et al., 2010) were raised from eggs within the BSF Animal Unit at the University of Manchester. Geese were housed within a 12 m2 room and provided with an area of turf and a 30 cm deep raised pond. They were maintained on a 13 h:11 h light:dark cycle and were given food ad libitum (Poultry Grower Pellets, Small Holder Range, Norfolk, UK), and were not fasted prior to the experiments. The geese were 6 months old when used in the experiments described here, at which point their body masses (mean 1.79 kg, range 1.5–2.1 kg) had stabilised and were comparable to those of an older 17 month old cohort (mean 1.75 kg, range 1.5–2.1 kg), which were previously reared under the same conditions. The geese were imprinted from hatching and from the age of 3 months were trained regularly to walk inside a respirometry chamber upon a Tunturi T60 (Turku, Finland) treadmill. All experimental procedures were conducted in accordance with the Animals (Scientific Procedures) Act 1986 under a UK Home Office Licence held by Dr Codd (40/3001) and were approved by the University of Manchester Ethics Committee.

Indirect calorimetry

A standard open-flow respirometry system was used to measure rates (ml min–1) of O2 consumption (O2) and CO2 production (CO2) (Lighton, 2008; Withers, 2001). The geese walked within a transparent walled chamber constructed of Perspex with internal dimensions of 66×48.5×46.5 cm and a volume of 148 l. The chamber was sealed upon the treadmill with rubber and brush seals. Air was pushed through the chamber by two vacuum pumps (models 2750CGH160 and 2688CHI44, Thomas, Sheboygan, WI, USA) at 255 l min–1. A high flow rate relative to the chamber size was used to ensure quick stabilisation of the measured gas volumes, reducing the amount of time the geese had to walk. Excurrent air was sub-sampled at 0.110 l min–1 downstream of the chamber and passed through an RH300 (Sable Systems International, Las Vegas, NV, USA) for measurement of water vapour pressure and relative humidity. The air was then scrubbed of H2O (using magnesium perchlorate) and passed through a CA-10A CO2 analyser (Sable Systems International), then scrubbed of CO2 (using Ascarite II) prior to passing through a dual absolute and differential O2 analyser (Oxzilla II, Sable Systems International). Therefore, the air entering the O2 analyser was both H2O and CO2 free. A second O2 analyser channel continuously measured O2 levels in ambient air. Before analysis by the second O2 channel, Ascarite II and magnesium perchlorate were used to scrub CO2 and H2O from the ambient air. The sub-sampling pumps (TR-SS3 mass flow meters, Sable Systems International) were positioned downstream of the gas analysers. Voltage outputs were recorded with a UI2 interface and ExpeData version 1.25 (both Sable Systems International). Scrubbing chemicals were supplied by Acros organics (Thermo Fisher Scientific, Fair Lawn, NJ, USA).

Because water was scrubbed from the sub-sampled air stream prior to O2 and CO2 analysis, the primary flow rate (FR) was adjusted to a dry-corrected flow rate (FRc) using equation 8.6 in Lighton (Lighton, 2008):
(1)
where Pb is barometric pressure and PH2O is water vapour pressure. O2 was calculated from equation 10.1 of Lighton (Lighton, 2008) as:
(2)
where FIO2 is the concentration of O2 flowing into the respirometry chamber and FEO2 is that measured in CO2-free air by the O2 analyser after leaving the chamber. The accuracy of the respirometry system was tested by leaching N2 into the respirometry chamber at 2 l min–1 (Fedak et al., 1981). The error in O2 measured was <1.96% with the treadmill belt running at both 1.00 and 2.00 m s–1 and the time to reach a stable reading was 163 s.

The metabolic rates of the geese were measured at five walking velocities (U=0.25, 0.50, 0.75, 1.00 and 1.25 m s–1) and at three gradients: 7 deg decline, 7 deg incline and on the level (i.e. 0 deg gradient). The data for walking on the level were obtained from a previous publication (Nudds et al., 2010). During the walking trials a companion bird was housed in a mesh-fronted holding crate and was visible to the walking goose. A mirror was also attached to the front of the respirometry chamber facing the walking goose. The companion goose and mirror reduced anxiety for the walking goose and ensured a steady walking performance. Each goose was walked at all of the five velocities on one of the gradients during a single experiment and the order of velocities and gradients used was randomised. The geese were made to walk until measured gas concentrations had stabilised (typically between 5 and 10 min). The maximum U sustainable by all the geese for a period long enough to obtain stable gas concentration readings on each gradient was 1.25 m s–1. At this speed the geese were still using a walking gait (Nudds et al., 2010). Between each walking U, before the first U measured and after the final U measured, the geese were allowed to rest until the gas concentrations recorded were stable (again 5–10 min). A stable reading when resting was taken to indicate recovery from the previous bout of walking. For each speed the O2 and CO2 concentrations were derived as the mean of the last 60 s (60 samples) of the period of stable readings.

Kinematics

A Sony HDR-SR8E video camera (25 frames s–1) was used to film the geese at each of the five walking speeds. The video footage was analysed using Tracker 2.51 (Open Source Physics). The position of the right foot (nearest to the camera) was tracked over five strides, to calculate the gait kinematics parameters: stride length (lstride), stride frequency (fstride), stance phase (tstance), swing phase (tswing) and duty factor.

Data analyses

Differences in the relationship between the dependent variables O2, respiratory exchange ratio (RER), kinematics variables and U for different gradients were investigated using general linear models (GLMs). Gradient was included as a factor, goose as a random factor and U as a covariate. Goose was included in the GLM, because of the repeated measures taken for each individual, i.e. each goose was measured once at each U and for each gradient (5 speeds × 3 gradients=15 measurements per goose). Initially, both the slopes and intercepts of the relationship between the dependent variable and U were tested for differences between gradients. If the interaction term (gradient × U) was non-significant, indicating similar slopes (i.e. the incremental change in the dependent variable with U was similar for each gradient), it was removed from the GLM and the GLM was rerun assuming parallel lines (similar slopes) and hence testing for a difference in the intercepts only. Where differences were found, post hoc Fisher's least significant difference (LSD) tests (α=0.05) were used to indicate specific differences between the gradients.

All statistical analyses were performed using IBM SPSS Statistics v.19. Results are displayed as means ± s.e.m. for O2 and RER, and as means ± s.d. for kinematics variables.

Rate of oxygen consumption

A GLM showed the slope of the relationship between O2 and U was similar for each gradient (Table 1). Simplifying the GLM by specifying a common slope (Fig. 1) indicated a difference in the intercepts of the relationship between O2 and U for the three gradients (Table 1). During incline walking, the geese had a higher O2 than when both decline and level walking across all speeds (Table 1, Fig. 1). Previous studies have not conducted GLM analyses. Instead, separate linear regression lines were used to determine the relationship between O2 and U for each gradient independently. Differences in the slopes of the O2 against U regressions were then taken as such without statistical validation. Similar separate regression analyses here resulted in lines of best fit for the three gradients tested of y=36.85+65.51x, r2=0.98, P<0.001 (level), y=34.40+70.68x, r2=0.98, P<0.001 (decline) and y=40.17+80.17x, r2=0.99, P<0.001 (incline).

The RERs were similar on each gradient (mean 0.69±0.01, 0.72±0.01 and 0.73±0.02 for decline, level and incline, respectively). RER also increased predictably with U and the incremental increase in RER with U was similar for each gradient (Table 1). The mean RER for each U in ascending order (0.25, 0.50, 0.75, 1.00 and 1.25 m s–1) was 0.67±0.03, 0.70±0.01, 0.71±0.02, 0.72±0.02 and 0.77±0.02.

Kinematics

With the exception of duty factor, the incremental change in each kinematics variable with U was similar among gradients (Table 1, Fig. 2). On all gradients, duty factor decreased with increasing U, but was higher at all speeds during incline walking than during decline and level walking, which were similar at all speeds (Table 1, Fig. 2A).

tstance and tswing also decreased with U at all gradients. tswing was of shorter duration during incline walking than when both decline and level walking (Table 1, Fig. 2B). Although a Fisher's LSD post hoc test did not detect a difference, tswing also appeared to be of shorter duration on a decline than when walking on the level. In contrast, tstance was of longer duration at all speeds during incline walking than when walking on level ground, which in turn showed a longer tstance than when walking on a decline (Table 1, Fig. 2B).

lstride was shorter and fstride quicker across all speeds during decline than during both incline and level walking, which were similar at all speeds (Table 1, Fig. 2C,D).

For each GLM performed there was a significant effect of goose upon the dependent variable (Table 1). The absolute values of kinematics variables and O2 differed between individuals. The overall (direction of) between-gradient effects and changes with U, however, were similar for each goose.

Across all speeds the metabolic cost (O2) of walking uphill in the barnacle geese was higher than during downhill and level walking, which had similar costs. The increase in metabolic cost when moving up an incline agrees with previous studies (Armstrong et al., 1983; Bamford and Maloiy, 1980; Bedford et al., 1979; Brooks and White, 1978; Chassin et al., 1976; Clapperton, 1964; Cohen et al., 1978; Dill, 1965; Eaton et al., 1995; Ellerby et al., 2003; Farley and Emshwiller, 1996; Margaria et al., 1963; McDonough et al., 2002; Minetti et al., 2002; Raab et al., 1976; Robbins et al., 1979; Rubenson et al., 2006; Snyder and Carello, 2008; Taylor et al., 1972; Warncke et al., 1988; White and Yousef, 1978; Wickler et al., 2000; Wickler et al., 2003; Wunder and Morrison, 1974; Yousef et al., 1972). The finding that descending a decline is not metabolically cheaper than locomoting on the level, however, is contrary to previous findings (Armstrong et al., 1983; Byrnes et al., 1985; Cohen et al., 1978; Margaria et al., 1963; Minetti et al., 2002; Raab et al., 1976; Robbins et al., 1979; Taylor et al., 1972; White and Yousef, 1978; Yousef et al., 1972). Likewise, the incremental cost of increasing speed was similar at each gradient, which is inconsistent with most previous studies that have found the slope of the relationship between O2 and U increases with incline angle and decreases with decline angle (e.g. Chassin et al., 1976; Clapperton, 1964; Cohen et al., 1978; Margaria et al., 1963; Raab et al., 1976; Taylor et al., 1972; White and Yousef, 1978; Wunder and Morrison, 1974; Yousef et al., 1972). The findings of the current study, then, pose two significant questions. First, why do the barnacle geese not save metabolic energy when walking downhill, and second, why is the incremental increase in O2 with U independent of gradient?

Table 1.

The results of the GLMs and associated Fisher's LSD statistical tests used to determine whether gradient (grad) differentially affected VO2, RER and each of the five kinematics variables

The results of the GLMs and associated Fisher's LSD statistical tests used to determine whether gradient (grad) differentially affected VO2, RER and each of the five kinematics variables
The results of the GLMs and associated Fisher's LSD statistical tests used to determine whether gradient (grad) differentially affected VO2, RER and each of the five kinematics variables

As highlighted by Gillis and Biewener (Gillis and Biewener, 2002), there is a dearth of studies providing comparative whole-limb kinematics data where animals are encouraged to walk on gradients (less for declines than inclines) and on level ground at similar speeds, and even fewer studies where metabolic data were collected simultaneously (Table 2). Consequently, making specific statements about the effects of gradients on limb kinematics may be premature. Nevertheless, it is possible to identify some broad trends (Table 2) and the data do hint at why in barnacle geese the metabolic costs are the same for decline as for level walking. The changes in barnacle goose kinematics when on a decline, compared with locomoting on level ground, are similar to those of the squirrel monkey (Vilensky et al., 1994) and the Sprague–Dawley rat (Gillis and Biewener, 2002) with no change in duty factor; increasing fstride and decreasing lstride and tswing. In contrast, however, the barnacle goose also showed reduced tstance. In humans, where both kinematics and metabolic data were collected, a negative gradient did not affect fstride or lstride, but O2 was reduced (Snyder and Farley, 2011). Similar kinematic changes were induced by declines in the cat (Carlson-Kuhta et al., 1998; Smith et al., 1998), turkey Meleagris gallopavo (Gabaldón et al., 2004) and iguanian lizard Dipsosaurus dorsalis (Jayne and Irschick, 1999), although concomitant metabolic data were not collected. A shorter, faster (‘choppy’) stride suggests that the barnacle geese were trying to get their feet down quickly to brake and avoid falling. Losing kinematic energy at each footfall requires negative work from the muscles, and the more energy absorbed (increased braking effort) the more metabolically costly it is to the organism. The assistance of gravity in helping move the CoM is probably negated by this need to brake and the increase in fstride (O2 tends to increase with fstride).

These geese are the only species so far to show an increase in duty factor and a decrease in tswing during incline locomotion. An animal must produce more force to ascend an incline at a similar speed to that on the level and it appears that how this is achieved may vary both inter- and intra-specifically (Table 2). The barnacle goose, rat (Gillis and Biewener, 2002) and some horses (van Oldruitenborgh-Ooste et al., 1997) require an increased tstance to generate the additional force, whereas other species appear to generate more force within the same tstance (Table 2). The barnacle goose maintains lstride and fstride at level walking intensities when ascending and so compensates for the increase in tstance by decreasing tswing, which contributes to an increase in duty factor. The muscles responsible for swinging the limb are accountable for 26% of the energy used in guinea fowl Numida meleagris (Marsh et al., 2004). Swinging the leg through more rapidly must increase the work done by the swing phase muscles and hence contributes to the increased metabolic cost. Although perhaps indicating a change in fuel source composition with increasing speed, the RER values suggest that, at the whole-animal level, the geese were always operating aerobically and gradient did not affect this.

Fig. 1.

Rate of oxygen consumption (VO2) plotted against velocity (U) for level, decline and incline walking. The data points represent the values averaged across all five geese and the lines of best fit were derived from the general linear model (GLM) and are O2=46.21+72.12U (incline, blue symbols and line), O2=31.89+72.12U (level, black symbols and line) and O2=33.32+72.12U (decline, red symbols and line).

Fig. 1.

Rate of oxygen consumption (VO2) plotted against velocity (U) for level, decline and incline walking. The data points represent the values averaged across all five geese and the lines of best fit were derived from the general linear model (GLM) and are O2=46.21+72.12U (incline, blue symbols and line), O2=31.89+72.12U (level, black symbols and line) and O2=33.32+72.12U (decline, red symbols and line).

Data available on the effects of gradient on leg kinematics where speed was also controlled are few, but suggest that differences between species persist and are probably driven by leggedness (bipedal or quadrupedal), ecological niche and postural effects (Table 2). The relative proportions of different leg muscles and their architecture are also likely to lead to species-specific responses to locomoting on gradients. Cats adopt a more crouched posture when negotiating both negative and positive gradients (Carlson-Kuhta et al., 1998; Smith et al., 1998) and lizard limbs are flexed more when travelling downhill (Jayne and Irschick, 1999) than when walking on the level. This reduction in effective limb length reduces hip height. Consequently, the lever arm between the foot and the CoM is reduced when going uphill. Similarly, crouching on a decline reduces the lever arm between the braking foot and CoM, which reduces the moment around the foot and presumably reduces the costs of braking. Although not measured directly, the geese did not appear to adjust their effective limb lengths on the gradients and a lack of plasticity in posture may be a contributor to the lack of a reduction in O2 for the geese walking downhill. The coefficient of friction may also be important in determining running performance uphill and downhill (Lammers et al., 2006). This must be dependent upon foot morphology and a webbed foot in geese may not provide the same level of traction as the clawed foot of a galliform. Posture will also determine how the foot is weighted and, consequently, influence traction. The timing and direction of weighting relative to the ground, and how much of the foot is in contact with the ground throughout stance, will affect traction.

Fig. 2.

Kinematics variables (A) duty factor, (B) stance time (tstance, filled circles) and swing time (tswing, crosses), (C) stride length (lstride) and (D) stride frequency (fstride) plotted against velocity (U) for incline (blue symbols), level (black symbols) and decline (red symbols) walking. The data points represent the values averaged across all five geese and the lines of best fit were derived from the GLMs, and are (A) duty factor: y=0.82–0.21U (incline, blue line), y=0.77–0.19U (level, black line) and y=0.79–0.23U (decline, red line), (B) stance: y=0.26U–0.88 (incline, blue line), y=0.25U–0.88 (level, black line) and y=0.24U–0.88 (decline, red line), swing: y=0.17U–0.30 (incline, blue line), y=0.19U–0.30 (level, black line) and 0.18U–0.30 (decline, red line), (C) lstride: y=0.44U0.33 (incline, black dashed line), y=0.44U0.33 (level, black dashed line) and 0.42U0.33 (decline, red line), and (D) fstride: y=0.49+1.79U (incline, black dashed line), y=0.48+1.79U (level, black dashed line) and y=0.57+1.79U (decline, red line). Because of the non-linear nature of the relationship between tstance, tswing and lstride and U, statistical analyses were conducted on log10-transformed data (Table 1). Note that where lines of best fit are not discernable either statistically or graphically between two different gradients a single black dashed line is plotted to represent both.

Fig. 2.

Kinematics variables (A) duty factor, (B) stance time (tstance, filled circles) and swing time (tswing, crosses), (C) stride length (lstride) and (D) stride frequency (fstride) plotted against velocity (U) for incline (blue symbols), level (black symbols) and decline (red symbols) walking. The data points represent the values averaged across all five geese and the lines of best fit were derived from the GLMs, and are (A) duty factor: y=0.82–0.21U (incline, blue line), y=0.77–0.19U (level, black line) and y=0.79–0.23U (decline, red line), (B) stance: y=0.26U–0.88 (incline, blue line), y=0.25U–0.88 (level, black line) and y=0.24U–0.88 (decline, red line), swing: y=0.17U–0.30 (incline, blue line), y=0.19U–0.30 (level, black line) and 0.18U–0.30 (decline, red line), (C) lstride: y=0.44U0.33 (incline, black dashed line), y=0.44U0.33 (level, black dashed line) and 0.42U0.33 (decline, red line), and (D) fstride: y=0.49+1.79U (incline, black dashed line), y=0.48+1.79U (level, black dashed line) and y=0.57+1.79U (decline, red line). Because of the non-linear nature of the relationship between tstance, tswing and lstride and U, statistical analyses were conducted on log10-transformed data (Table 1). Note that where lines of best fit are not discernable either statistically or graphically between two different gradients a single black dashed line is plotted to represent both.

Previous studies have found the slope of the line of best fit between O2 and U to increase with gradient (Armstrong et al., 1983; Bamford and Maloiy, 1980; Chassin et al., 1976; Clapperton, 1964; Cohen et al., 1978; Eaton et al., 1995; Farley and Emshwiller, 1996; Full and Tullis, 1990; Margaria et al., 1963; McDonough et al., 2002; Raab et al., 1976; Robbins et al., 1979; Taylor et al., 1972; Warncke et al., 1988; White and Yousef, 1978; Wunder and Morrison, 1974; Yousef et al., 1972) and the established way of calculating the efficiency of moving along gradients for comparison between different species is dependent upon there being a difference (Taylor et al., 1972).

The net metabolic cost of lifting 1 kg of body mass 1 m in J is calculated using:
(3)
where S is the slope of the line of best fit between O2 and U for gradient exercising, S0 is the same for locomoting on a level substrate, θ is the angle of incline and 20.10 is the thermal equivalent for conversion of O2 into J (Taylor et al., 1972). Efficiency (η) is mechanical work/energy expended or energy recovered/mechanical energy stored for uphill and downhill locomoting, respectively.
Table 2.

Trends in the metabolic cost and kinematic responses to locomoting on a decline or incline compared with on a level surface for a range of species

Trends in the metabolic cost and kinematic responses to locomoting on a decline or incline compared with on a level surface for a range of species
Trends in the metabolic cost and kinematic responses to locomoting on a decline or incline compared with on a level surface for a range of species

The incremental cost (O2) of walking on gradients in the barnacle goose, however, was independent of gradient angle. An increasing slope of O2 against U with increasing gradient is expected, if the change in cost is due solely to the external mechanical work associated with increasing the animal's potential energy (Ep=MgUsinθ, where M is body mass, g is the acceleration due to gravity and θ is the angle of the gradient with the horizontal). Only if η is constant across all speeds would a similar pattern in O2 variation with U to that predicted for Ep be expected. Alternatively, if η increased with U on inclines and decreased with U on declines then similar slopes between O2 and U for each gradient could result. Such speed and gradient variations in η were previously found for sheep (Clapperton, 1964).

The previous findings of increasing incremental costs with increasing gradient may be a real effect or alternatively a consequence of fitting separate regression lines to the data for each gradient and not using ANCOVA to test whether the slopes of these lines really (statistically) differ. For example, fitting independent linear regression analyses to the O2/U relationships for each gradient here indicated an increased slope for incline walking. Another confounding factor could stem from the fact that most studies have not determined the gaits used by the study animals and have included all the data together, which may have skewed the line(s) of best fit. The relationship between O2 and U is not necessarily linear (Nudds et al., 2009) and the nature of it may also differ between gaits (Nudds et al., 2011). There is also evidence that the effect of gradient upon limb kinematics is gait dependent. When galloping on an incline, horses increase their fstride and increase the work per stride for a given U when compared with similar speeds during galloping on level ground (Parsons et al., 2008). When trotting, however, they only increase the work per cycle, leaving fstride unchanged (Wickler et al., 2005). Alterations in joint angles and joint articulation angles are also dependent upon U as well as gradient (Carlson-Kuhta et al., 1998; Jayne and Irschick, 1999; Lammers et al., 2006; Leroux et al., 1999; Schmidt and Fischer, 2011; Smith et al., 1998; Vilensky et al., 1994). Perhaps, if only walking gaits were analysed in previous studies the slope of the relationship between O2 and U would also be independent of gradient, as found for the barnacle goose.

The lack of variation in the incremental increase in O2 with U between gradients means a comparison with other studies to determine whether the increased costs experienced by the geese are relatively higher is not possible. Therefore, only one of the two original predictions is supported, i.e. the geese saved less energy descending when compared with walking on the level than previous species studied. In addition, they showed a novel kinematic response to gradients. Raab and colleagues argued that there is little difference in the mechanical aspects of grade locomotion between bipeds and quadrupeds as indicated by similar net metabolism values at similar grades for humans and dogs (Raab et al., 1976). In subsequent work involving a broader range of species, Gillis and Biewener began to identify general patterns among species, but also postulate that species-specific differences in kinematic responses to gradients were likely (Gillis and Biewener, 2002). For the first time, the present study has identified no energy saving for decline running and an increase in duty factor when ascending a gradient. The barnacle goose data reported here therefore add important novel information to existing data and reinforce the notion that a suite of factors (such as ecological niche, posture, gait, leggedness and foot morphology) are likely to affect an organism's ability to negotiate gradients.

FUNDING

This work was supported by a research grant from the Leverhulme Trust [grant no. F/00 120/BH].

We thank J. D. Gardiner and P. G. Tickle for assistance during the experiments. We thank the two anonymous referees for their helpful comments.

Armstrong
R. B.
,
Laughlin
M. H.
,
Rome
L.
,
Taylor
C. R.
(
1983
).
Metabolism of rats running up and down an incline
.
J. Appl. Physiol.
55
,
518
-
521
.
Bamford
O. S.
,
Maloiy
G. M. O.
(
1980
).
Energy metabolism and heart rate during treadmill exercise in the Marabou stork
.
J. Appl. Physiol.
49
,
491
-
496
.
Bedford
T. G.
,
Tipton
C. M.
,
Wilson
N. C.
,
Oppliger
R. A.
,
Gisolfi
C. V.
(
1979
).
Maximum oxygen consumption of rats and its changes with various experimental procedures
.
J. Appl. Physiol.
47
,
1278
-
1283
.
Brooks
G. A.
,
White
T. P.
(
1978
).
Determination of metabolic and heart rate responses of rats to treadmill exercise
.
J. Appl. Physiol.
45
,
1009
-
1015
.
Byrnes
W. C.
,
Clarkson
P. M.
,
White
J. S.
,
Hsieh
S. S.
,
Frykman
P. N.
,
Maughan
R. J.
(
1985
).
Delayed onset muscle soreness following repeated bouts of downhill running
.
J. Appl. Physiol.
59
,
710
-
715
.
Carlson-Kuhta
P.
,
Trank
T. V.
,
Smith
J. L.
(
1998
).
Forms of forward quadrupedal locomotion. II. A comparison of posture, hindlimb kinematics, and motor patterns for upslope and level walking
.
J. Neurophysiol.
79
,
1687
-
1701
.
Chassin
P. S.
,
Taylor
C. R.
,
Heglund
N. C.
,
Seeherman
H. J.
(
1976
).
Locomotion in lions – energetic cost and maximum aerobic capacity
.
Physiol. Zool.
49
,
1
-
10
.
Clapperton
J. L.
(
1964
).
The energy metabolism of sheep walking on the level and on gradients
.
Br. J. Nutr.
18
,
47
-
54
.
Cohen
Y.
,
Robbins
C. T.
,
Davitt
B. B.
(
1978
).
Oxygen utilization by elk calves during horizontal and vertical locomotion compared to other species
.
Comp. Biochem. Physiol.
61A
,
43
-
48
.
de Leon
R.
,
Hodgson
J. A.
,
Roy
R. R.
,
Edgerton
V. R.
(
1994
).
Extensor- and flexor-like modulation within motor pools of the rat hindlimb during treadmill locomotion and swimming
.
Brain Res.
654
,
241
-
250
.
Dill
D. B.
(
1965
).
Oxygen used in horizontal and grade walking and running on the treadmill
.
J. Appl. Physiol.
20
,
19
-
22
.
Eaton
M. D.
,
Evans
D. L.
,
Hodgson
D. R.
,
Rose
R. J.
(
1995
).
Effect of treadmill incline and speed on metabolic rate during exercise in thoroughbred horses
.
J. Appl. Physiol.
79
,
951
-
957
.
Ellerby
D. J.
,
Cleary
M. E.
,
Marsh
R. L.
,
Buchanan
C. I.
(
2003
).
Measurement of maximum oxygen consumption in Guinea fowl Numida meleagris indicates that birds and mammals display a similar diversity of aerobic scopes during running
.
Physiol. Biochem. Zool.
76
,
695
-
703
.
Farley
C.
,
Emshwiller
M.
(
1996
).
Efficiency of uphill locomotion in nocturnal and diurnal lizards
.
J. Exp. Biol.
199
,
587
-
592
.
Fedak
M. A.
,
Rome
L.
,
Seeherman
H. J.
(
1981
).
One-step N2-dilution technique for calibrating open-circuit VO2 measuring systems
.
J. Appl. Physiol.
51
,
772
-
776
.
Full
R. J.
,
Tullis
A.
(
1990
).
Energetics of ascent: insects on inclines
.
J. Exp. Biol.
149
,
307
-
317
.
Gabaldón
A. M.
,
Nelson
F. E.
,
Roberts
T. J.
(
2004
).
Mechanical function of two ankle extensors in wild turkeys: shifts from energy production to energy absorption during incline versus decline running
.
J. Exp. Biol.
207
,
2277
-
2288
.
Gardiner
K. R.
,
Gardiner
P. F.
,
Edgerton
V. R.
(
1982
).
Guinea pig soleus and gastrocnemius electromyograms at varying speeds, grades, and loads
.
J. Appl. Physiol.
52
,
451
-
457
.
Gillis
G. B.
,
Biewener
A. A.
(
2002
).
Effects of surface grade on proximal hindlimb muscle strain and activation during rat locomotion
.
J. Appl. Physiol.
93
,
1731
-
1743
.
Herreid
C. F.
,
Full
R. J.
,
Prawel
D. A.
(
1981
).
Energetics of cockroach locomotion
.
J. Exp. Biol.
94
,
189
-
202
.
Jayne
B. C.
,
Irschick
D. J.
(
1999
).
Effects of incline and speed on the three-dimensional hindlimb kinematics of a generalized iguanian lizard (Dipsosaurus dorsalis)
.
J. Exp. Biol.
202
,
143
-
159
.
Kai
M.
,
Hiraga
A.
,
Kubo
K.
,
Tokurik
M.
(
1997
).
Comparison of stride characteristics in a cantering horse on a flat and inclined treadmill
.
Equine Vet. J.
29
Suppl.
,
76
-
79
.
Lammers
A. R.
,
Earls
K. D.
,
Biknevicius
A. R.
(
2006
).
Locomotor kinetics and kinematics on inclines and declines in the gray short-tailed opossum Monodelphis domestica
.
J. Exp. Biol.
209
,
4154
-
4166
.
Leroux
A.
,
Fung
J.
,
Barbeau
H.
(
1999
).
Adaptation of the walking pattern to uphill walking in normal and spinal-cord injured subjects
.
Exp. Brain Res.
126
,
359
-
368
.
Lighton
J. R. B.
(
2008
).
Measuring Metabolic Rates: a Manual for Scientists
.
New York
:
Oxford University Press
.
Lipp
A.
,
Wolf
H.
,
Lehmann
F. O.
(
2005
).
Walking on inclines: energetics of locomotion in the ant Camponotus
.
J. Exp. Biol.
208
,
707
-
719
.
Margaria
R.
,
Cerretelli
P.
,
Aghemo
P.
,
Sassi
G.
(
1963
).
Energy cost of running
.
J. Appl. Physiol.
18
,
367
-
370
.
Marsh
R. L.
,
Ellerby
D. J.
,
Carr
J. A.
,
Henry
H. T.
,
Buchanan
C. I.
(
2004
).
Partitioning the energetics of walking and running: swinging the limbs is expensive
.
Science
303
,
80
-
83
.
McDonough
P.
,
Kindig
C. A.
,
Ramsel
C.
,
Poole
D. C.
,
Erickson
H. H.
(
2002
).
The effect of treadmill incline on maximal oxygen uptake, gas exchange and the metabolic response to exercise in the horse
.
Exp. Physiol.
87
,
499
-
506
.
Minetti
A. E.
,
Moia
C.
,
Roi
G. S.
,
Susta
D.
,
Ferretti
G.
(
2002
).
Energy cost of walking and running at extreme uphill and downhill slopes
.
J. Appl. Physiol.
93
,
1039
-
1046
.
Nudds
R. L.
,
Codd
J. R.
,
Sellers
W. I.
(
2009
).
Evidence for a mass dependent step-change in the scaling of efficiency in terrestrial locomotion
.
PLoS ONE
4
,
e6927
.
Nudds
R. L.
,
Gardiner
J. D.
,
Tickle
P. G.
,
Codd
J. R.
(
2010
).
Energetics and kinematics of walking in the barnacle goose (Branta leucopsis)
.
Comp. Biochem. Physiol.
156A
,
318
-
324
.
Nudds
R. L.
,
Folkow
L. P.
,
Lees
J. J.
,
Tickle
P. G.
,
Stokkan
K. A.
,
Codd
J. R.
(
2011
).
Evidence for energy savings from aerial running in the Svalbard rock ptarmigan (Lagopus muta hyperborea)
.
Proc. R. Soc. London. B
278
,
2654
-
2661
.
Parsons
K. J.
,
Pfau
T.
,
Ferrari
M.
,
Wilson
A. M.
(
2008
).
High-speed gallop locomotion in the Thoroughbred racehorse. II. The effect of incline on centre of mass movement and mechanical energy fluctuation
.
J. Exp. Biol.
211
,
945
-
956
.
Pierotti
D. J.
,
Roy
R. R.
,
Gregor
R. J.
,
Edgerton
V. R.
(
1989
).
Electromyographic activity of cat hindlimb flexors and extensors during locomotion at varying speeds and inclines
.
Brain Res.
481
,
57
-
66
.
Raab
J. L.
,
Eng
P.
,
Waschler
R. A.
(
1976
).
Metabolic cost of grade running in dogs
.
J. Appl. Physiol.
41
,
532
-
535
.
Robbins
C. T.
,
Cohen
Y.
,
Davitt
B. B.
(
1979
).
Energy expenditure by elk calves
.
J. Wildl. Manage.
43
,
445
-
453
.
Robert
C.
,
Valette
J. P.
,
Denoix
J. M.
(
2000
).
The effects of treadmill inclination and speed on the activity of two hindlimb muscles in the trotting horse
.
Equine Vet. J.
32
,
312
-
317
.
Roberts
T. J.
,
Higginson
B. K.
,
Nelson
F. E.
,
Gabaldón
A. M.
(
2007
).
Muscle strain is modulated more with running slope than speed in wild turkey knee and hip extensors
.
J. Exp. Biol.
210
,
2510
-
2517
.
Roy
R. R.
,
Hutchison
D. L.
,
Pierotti
D. J.
,
Hodgson
J. A.
,
Edgerton
V. R.
(
1991
).
EMG patterns of rat ankle extensors and flexors during treadmill locomotion and swimming
.
J. Appl. Physiol.
70
,
2522
-
2529
.
Rubenson
J.
,
Henry
H. T.
,
Dimoulas
P. M.
,
Marsh
R. L.
(
2006
).
The cost of running uphill: linking organismal and muscle energy use in guinea fowl (Numida meleagris)
.
J. Exp. Biol.
209
,
2395
-
2408
.
Schmidt
A.
,
Fischer
M. S.
(
2011
).
The kinematic consequences of locomotion on sloped arboreal substrates in a generalized (Rattus norvegicus) and a specialized (Sciurus vulgaris) rodent
.
J. Exp. Biol.
214
,
2544
-
2559
.
Smith
J. L.
,
Carlson-Kuhta
P.
(
1995
).
Unexpected motor patterns for hindlimb muscles during slope walking in the cat
.
J. Neurophysiol.
74
,
2211
-
2215
.
Smith
J. L.
,
Carlson-Kuhta
P.
,
Trank
T. V.
(
1998
).
Forms of forward quadrupedal locomotion. III. A comparison of posture, hindlimb kinematics, and motor patterns for downslope and level walking
.
J. Neurophysiol.
79
,
1702
-
1716
.
Snyder
G. K.
,
Carello
C. A.
(
2008
).
Body mass and the energy efficiency of locomotion: lessons from incline running
.
Comp. Biochem. Physiol.
150A
,
144
-
150
.
Snyder
K. L.
,
Farley
C. T.
(
2011
).
Energetically optimal stride frequency in running: the effects of incline and decline
.
J. Exp. Biol.
214
,
2089
-
2095
.
Taylor
C. R.
,
Caldwell
S. L.
,
Rowntree
V. J.
(
1972
).
Running up and down hills: some consequences of size
.
Science
178
,
1096
-
1097
.
Tullis
A.
,
Andrus
S.
(
2011
).
The cost of incline locomotion in ghost crabs (Ocypode quadrata) of different sizes
.
J. Comp. Physiol. B
181
,
873
-
881
.
van Oldruitenborgh-Ooste
S.
,
Barneveld
A.
,
Schamhardt
H. C.
(
1997
).
Effects of treadmill inclination on kinematics of the trot in Dutch warmblood horses
.
Equine Vet. J. Suppl.
29
,
71
-
75
.
Vilensky
J. A.
,
Moore
A. M.
,
Libii
J. N.
(
1994
).
Squirrel monkey locomotion on an inclined treadmill: Implications for the evolution of gaits
.
J. Hum. Evol.
26
,
375
-
386
.
Warncke
G.
,
Bandholtz
J.
,
Schultzemotel
P.
(
1988
).
Metabolic cost and body temperatures during grade running in quail (Coturnix coturnix)
.
Comp. Biochem. Physiol.
89A
,
93
-
96
.
White
R. G.
,
Yousef
M. K.
(
1978
).
Energy-expenditure in reindeer walking on roads and on tundra
.
Can. J. Zool.
56
,
215
-
223
.
Wickler
S. J.
,
Hoyt
D. F.
,
Cogger
E. A.
,
Hirschbein
M. H.
(
2000
).
Preferred speed and cost of transport: the effect of incline
.
J. Exp. Biol.
203
,
2195
-
2200
.
Wickler
S. J.
,
Hoyt
D. F.
,
Cogger
E. A.
,
Myers
G.
(
2003
).
The energetics of the trot-gallop transition
.
J. Exp. Biol.
206
,
1557
-
1564
.
Wickler
S. J.
,
Hoyt
D. F.
,
Biewener
A. A.
,
Cogger
E. A.
,
De La Paz
K. L.
(
2005
).
In vivo muscle function vs speed. II. Muscle function trotting up an incline
.
J. Exp. Biol.
208
,
1191
-
1200
.
Withers
P. C.
(
2001
).
Design, calibration and calculation for flow-through respirometry systems
.
Aust. J. Zool.
49
,
445
-
461
.
Wunder
B. A.
,
Morrison
P. R.
(
1974
).
Red squirrel metabolism during incline running
.
Comp. Biochem. Physiol.
48A
,
153
-
161
.
Yousef
M. K.
,
Dill
D. B.
,
Freeland
D. V.
(
1972
).
Energetic cost of grade walking in man and burro, Equus asinus: desert and mountain
.
J. Appl. Physiol.
33
,
337
-
340
.