Despite the high mechanical demands associated with climbing, the ability to ascend vertically has evolved independently in most major animal lineages. However, little is known about the kinetics, mechanical energy profiles or spatiotemporal gait characteristics of this locomotor mode. In this study, we explored the dynamics of horizontal locomotion and vertical climbing on both flat substrates and narrow poles in five Australian green tree frogs (Ranoidea caerulea). Vertical climbing is associated with slow, deliberate movements (i.e. reduced speed and stride frequency and increased duty factors) and propulsive fore–aft impulses in both the forelimb and hindlimb. By comparison, horizontal walking was characterized by a braking forelimb and a propulsive hindlimb. In the normal plane, tree frogs mirrored other taxa in exhibiting a net pulling forelimb and a net pushing hindlimb during vertical climbing. In terms of mechanical energy, tree frogs matched theoretical predictions of climbing dynamics (i.e. the total mechanical energetic cost of vertical climbing was predominantly driven by potential energy, with negligible kinetic contributions). Utilizing power as a means of estimating efficiency, we also demonstrate that Australian green tree frogs show total mechanical power costs only slightly above the minimum mechanical power necessary to climb, highlighting their highly effective locomotor mechanics. This study provides new data on climbing dynamics in a slow-moving arboreal tetrapod and raises new testable hypotheses about how natural selection can act upon a locomotor behavior that is notably constrained by external physical forces.

Navigating an arboreal environment requires extreme locomotor flexibility to successfully traverse substrates of varying diameter, compliance, orientation and continuity (e.g. Grand, 1972; Jenkins, 1974; Cartmill, 1985; Preuschoft, 2002). Extensive efforts have been made to describe the biomechanics of arboreal locomotion in tetrapods, typically focusing on horizontal movement in mammals (e.g. Jenkins, 1974; Larson, 1998; Lemelin and Cartmill, 2010; Lemelin et al., 2003; Schmitt, 1999), squamates (e.g. Higham and Jayne, 2004; Fischer et al., 2010) or birds (e.g. Reader et al., 2022; Young et al., 2022a,b). Less understood are the mechanisms associated with ascending vertically oriented substrates (i.e. vertical climbing; but see Autumn et al., 2006; Granatosky et al., 2019b; Hanna et al., 2017; Wang et al., 2015; Young et al., 2022a). Unlike horizontal movement, in which the animal's body mass is supported by the substrate, vertical climbing necessitates the consistent opposition against gravitational acceleration to prevent falling (Cartmill, 1985; Preuschoft, 2002), resulting in high energetic costs. In spite of this, the ability to climb has evolved independently in numerous lineages of limbed animals (e.g. Cartmill, 1985; Goldman et al., 2006; Labonte and Federle, 2015, 2016).

Mechanics of climbing

Although behavioral modes of climbing have been broadly described across myriad taxa (e.g. Cartmill, 1985; Granatosky et al., 2022; Hill et al., 2018; Hunt et al., 1996; Preuschoft, 2002), quantitative data regarding either the spatiotemporal gait characteristics, kinetics or mechanical energy profiles of climbing are lacking (e.g. Autumn et al., 2006; Bock and Winkler, 1978; Goldman et al., 2006; Theunissen et al., 2014). Broadly, vertical climbing is associated with slow, deliberate movements (i.e. reduced speed and stride frequency and increased stride lengths, contact times, duty factor and overall stride duration; see Granatosky et al., 2019b; Hanna et al., 2022; Herrel et al., 2013; Karantanis et al., 2018; Theunissen et al., 2014). Further, animals are thought to prioritize long strides to regulate speed as a mechanism for decreasing substrate oscillations on thin arboreal supports (Granatosky and McElroy, 2022; Granatosky et al., 2019b; Hanna et al., 2022; Karantanis et al., 2018). These strategies are thought to increase stability and reduce substrate oscillations on thin arboreal substrates (Granatosky and McElroy, 2022; Granatosky et al., 2019b; Hanna et al., 2022; Herrel et al., 2013; Karantanis et al., 2018).

Single limb forces (SLFs) collected during vertical climbing are reported for several primate taxa (Hanna et al., 2017; Hirasaki et al., 1993), parrots (Reader et al., 2022; Young et al., 2022a), geckos (Autumn et al., 2006; Wang et al., 2015) and cockroaches (Clemente and Federle, 2008; Goldman et al., 2006; Watson et al., 2002). It is important to note that owing to the reorientation of the body axis associated with vertical climbing, these forces are not directly comparable to horizontal locomotion (Hanna et al., 2017). Fore–aft forces, which traditionally move an animal forward and backward along a horizontal plane, become upward and downward forces, respectively. Meanwhile, normal forces, which are typically imparted into the ground during horizontal walking, become pushing and pulling forces directed toward and away from the support. These normal forces are of particular importance during vertical ascent as climbing inherently creates a gravitational torque pitching the animal away from the support (Norberg, 1986; Preuschoft, 2002; Spring, 1965). Climbing animals are therefore required to generate and balance pulling forces generated by a rostrally positioned limb loaded in tension and pushing forces generated by a caudally positioned limb loaded in compression, such that the sum of the normal forces is zero (Hanna et al., 2017; Norberg, 1986; Preuschoft, 2002; Spring, 1965). Hexapedal insects broadly match the SLF patterns described above (Cruse, 1976; Goldman et al., 2006). Goldman et al. (2006) provided the most in-depth analysis of climbing cockroaches, and determined that all limbs produce propulsive fore–aft forces. Like quadrupedal climbers (Autumn et al., 2006; Hanna et al., 2017; Wang et al., 2015), cockroaches countered pitchback moments due to gravity by using the front-leg normal forces to pull the body toward the substrate while the hindlimbs push away from the substrate. Middle limb-pairs in cockroaches generate small normal forces (Goldman et al., 2006).

Movements of the center of mass (COM) are also altered by climbing. During horizontal locomotion in walking gaits, the limbs act as struts of an inverted pendulum and the COM sways from posterior to the limb at touchdown before assuming an anterior position at lift-off (Ahn et al., 2004; Bertram, 2016; Cavagna et al., 1977; Granatosky et al., 2018; Reilly et al., 2006; Riskin et al., 2016). Through conversion of potential energy to kinetic energy, mechanical energy can be recovered between strides, with the magnitude of energy recovered dependent upon the phasic relationship between potential and kinetic energies (Ahn et al., 2004; Blickhan and Full, 1987; Cavagna et al., 1977; O'Neill and Schmitt, 2012). During climbing, however, this phasic relationship ceases to exist. As animals gain height while climbing, potential energy steadily increases (Autumn et al., 2006; Goldman et al., 2006). While ascent continues, it is therefore impossible for the animal to recover energy within the stride by converting the gained potential energy into kinetic energy. Thus, measurements of mechanical efficiency during climbing are instead derived by comparing total mechanical power with the product of gravity and speed: the theoretical minimum power needed for ascent (Autumn et al., 2006; Goldman et al., 2006). Such patterns appear to hold true regardless of the number of limbs or adhesion mechanism (Autumn et al., 2006; Goldman et al., 2006), and these similarities led Goldman et al. (2006) to propose that the physical constraints of climbing are so great that animals are restricted to a singular global template to effectively ascend vertical surfaces.

Climbing in tree frogs

Tree frogs have evolved across multiple anuran families (Chan et al., 2018; Peloso et al., 2018; Richards and Moore, 1996; Wiens et al., 2005). As the name implies, these frogs are typically found in trees or other high-growing vegetation and do not normally descend to the ground, except to mate and spawn. Some species eschew even this, instead building foam nests on leaves and thus rarely leaving the trees at all as adults (Granatosky and Krysko, 2011; Kasuya et al., 1996). Tree frogs are typically smaller and more slender than terrestrial frogs, given that their weight must be supported by the branches and twigs in their habitats (Herrel et al., 2013). These small branches are typically traversed via digital grasping, and tree frogs have been shown to generate strong intrinsic grips around thin supports (Herrel et al., 2013; Hill et al., 2018; Manzano et al., 2008). However, larger diameter substrates are typically negotiated via the use of wet adhesion clinging (Federle et al., 2006). Although numerous studies have explored characteristics of arboreal locomotion in tree frogs (Endlein et al., 2017; Herrel et al., 2013; Hill et al., 2018; Kamada et al., 2017; Manzano et al., 2008), there remains no information about patterns of SLFs or overall movements of the COM. Here, we present impulse data (as derived from single limb force kinetics) from the Australian tree frog across four locomotor conditions (horizontal flat, horizontal pole, vertical flat and vertical pole). We elected to present impulses rather than individual single limb forces, as impulses better reflect the function of the limb during locomotion (Dickinson et al., 2022; Larson and Demes, 2011; Raichlen et al., 2009). We also present data on COM mechanical energy across horizontal flat and vertical flat conditions. Thus, we test the extent to which climbing in tree frogs mirrors patterns seen in other climbing species (Autumn et al., 2006; Goldman et al., 2006; Hanna et al., 2017; Hirasaki et al., 1993; Wang et al., 2015). Specifically, we assess the following three predictions.

Prediction 1

Vertical climbing will be associated with slow, deliberate movements (i.e. reduced speed and stride frequency and increased stride lengths, contact times, duty factor and overall stride duration). Further, animals will increase stride length, as opposed to stride frequency, when increasing speed as a mechanism for decreasing substrate oscillations on thin arboreal supports (Granatosky and McElroy, 2022; Granatosky et al., 2019b; Hanna et al., 2022; Karantanis et al., 2018).

Prediction 2

The functional role of each limb pair (e.g. forelimb versus hindlimb) will change as animals shift from horizontal walking to vertical climbing. Namely, both forelimbs and hindlimbs will adopt propulsive fore–aft roles during climbing, and they will differentiate into a pulling forelimb and pushing hindlimb (Autumn et al., 2006; Goldman et al., 2006; Hanna et al., 2017; Wang et al., 2015). Mediolateral forces will not be altered in response to vertical climbing, but will differ from primarily laterally directed on flat substrates to primarily medially directed on poles (Carlson et al., 2005; Demes et al., 2006; Franz et al., 2005; Schmitt, 2003).

Prediction 3

Total mechanical power will be greater during vertical climbing versus horizontal walking owing to the unavoidable accumulation of height and subsequently potential energy oscillations (Autumn et al., 2006).

Spatiotemporal gait variables, SLFs and COM movements were collected from five Australian green tree frogs [Ranoidea caerulea (White 1790); body mass range 26.5–48.0 g, average body mass=34.1±7.6 g] at the New York Institute of Technology College of Osteopathic Medicine (Old Westbury, NY, USA). All data collection protocols were approved by the Institutional Animal Care and Use Committee (IACUC protocol: MCG-2019-01). Animal subjects were adults and free from any visible gait pathologies.

Spatiotemporal gait data collection and processing

Spatiotemporal variables (i.e. contact time, total stride time, duty factor, stride length and stride frequency) were measured from lateral video footage of locomotor trials. Touchdown, lift-off and subsequent touchdown frames of each stride were recorded and divided by the recording frame rate to obtain time. Contact time was calculated by taking the difference between initial touchdown and lift-off. Duty factor was calculated by dividing contact time by total stride time. Stride length was calculated by multiplying speed by stride time. Stride frequency is the inverse of stride duration.

SLF data collection and processing

SLFs were collected across four experimental conditions, namely a flat runway and a pole each oriented at either a 0 deg (flat) or 90 deg (vertical) angle relative to the ground. Both the flat and pole runways measured 0.508 m in length, with widths of 0.152 m (flat) and 0.012 m (pole), respectively. Each runway featured a centrally positioned cut-out measuring 0.045 m in length, large enough to accommodate the palm length and sole length of the frog's hand or foot, while limiting the potential for overlapping handfalls or footfalls. These cut-outs were instrumented with a small-load force plate (model HE6X6; Advanced Mechanical Technology, Inc., Watertown, MA, USA) using a 3D printed platform (see Fig. 1). All trials began by zeroing out the force plate to remove any drift or offsets from previous recordings. In most cases, minimal motivation, such as a gentle tap, was provided to the frogs to encourage movement across both flat and cylindrical platforms. Only trials in which the frogs moved with no visual acceleration or deceleration, and included a clear handfall or footfall, were considered successful. All trials were recorded using two high-speed cameras (XC-2; Xcitex Inc., Woburn, MA, USA) mounted to capture frog movement at 125 Hz from lateral and posterior views. Force data (sampled at 1250 Hz) and video data were synchronized using Procapture (Xcitex). Stride length, stride duration and speed were calculated within ImageJ (Schneider et al., 2012) using the position of the frog's eye and the square-grid background of known distance (0.1 m) to calibrate the space. A stride was a defined as forelimb touchdown to subsequent touchdown of the same forelimb.

Fig. 1.

Experimental setup for single limb force (top) and center of mass (bottom) data collection. Schematic of (A) flat runway oriented horizontally and (B) cylindrical runway oriented vertically instrumented with an AMTI small-load force plate (FP; model HE6X6; load capacity=4 N; Watertown, MA, USA). The majority of the runway is mounted separately from the force plate (grey regions), although small segments (dark regions) are instrumented with force plate. Schematic of flat runway oriented (C) horizontally and (D) vertically instrumented with an AMTI small-load force plate (same as above). The runway is large enough to allow multiple strides. Each panel is accompanied by snapshot of Australian green tree frog (Ranoidea caerulea) traversing that experimental setup.

Fig. 1.

Experimental setup for single limb force (top) and center of mass (bottom) data collection. Schematic of (A) flat runway oriented horizontally and (B) cylindrical runway oriented vertically instrumented with an AMTI small-load force plate (FP; model HE6X6; load capacity=4 N; Watertown, MA, USA). The majority of the runway is mounted separately from the force plate (grey regions), although small segments (dark regions) are instrumented with force plate. Schematic of flat runway oriented (C) horizontally and (D) vertically instrumented with an AMTI small-load force plate (same as above). The runway is large enough to allow multiple strides. Each panel is accompanied by snapshot of Australian green tree frog (Ranoidea caerulea) traversing that experimental setup.

Single limb forces were processed using a custom-written MATLAB script to correct for direction of travel, orientation and contacting limb (right or left). Forces were standardized such that fore–aft forces were split into negative braking and positive propulsive forces, normal forces were divided into positive pushing and negative pulling forces, and mediolateral forces were divided into negative medially directed forces and positive laterally directed forces. All forces were filtered through a low-pass Fourier filter at 15 Hz (Hanna et al., 2017; Nyakatura et al., 2014; Ogihara et al., 2014) and normalized to a percentage of the animal's body weight (% BW) for peak forces and body weight seconds (% BWS) for impulses to allow statistical comparison between individuals (see below). All net impulses were calculated as the area under the force–time curve in the respective orientation. Although peak force data are reported in Table S1, for ease of readability and to reduce the number of statistical tests, we have elected to only describe single limb force data in terms of net impulse. Net impulse was preferable in the present study as it provides a meaningful summary of the limb's overall role during locomotion (for review, see Demes et al., 1994; Granatosky et al., 2018; Larson and Demes, 2011; Raichlen et al., 2009). For example, a limb's functional braking or propulsive role has been best represented by net fore–aft impulse in numerous studies (e.g. Byron et al., 2017; Demes et al., 1994; Granatosky et al., 2018). Further, weight support distribution between the limbs during horizontal walking has often been analyzed in terms of net normal impulse rather than peaks (Dickinson et al., 2022; Larson and Demes, 2011; Raichlen et al., 2009).

COM dynamics data collection and processing

To collect COM movements, animals walked or climbed across a fully instrumented (HE6X6 force plate as above) flat runway 0.508 m in length and 0.152 m in width. The pole substrate was not used for COM data because the frogs often twisted around the pole, affecting directionality of force readings. Data were sampled, recorded and synchronized using the same methods as above.

Energy fluctuations of the COM were calculated following previously published methodology (Autumn et al., 2006; Blickhan and Full, 1987; Cavagna et al., 1977; Heglund et al., 1982; Reilly et al., 2006). In summary, COM data were collected as the animal moved across a horizontal and vertical flat instrumented runway long enough to ensure multiple strides (Fig. 1). From the resulting kinetic data, the product of mass and acceleration was subtracted from forces in the plane of gravity (from normal forces in horizontal trials and from fore–aft forces in vertical trials) and subsequently divided by the animal's body mass to arrive at the instantaneous acceleration of the COM. Detrended acceleration (see Bishop et al., 2008; Cole et al., 2014; Gutierrez-Farewik et al., 2006; Ihlen et al., 2013; Walker, 1998) was integrated to calculate instantaneous velocity and secondarily integrated to identify the position of the COM. To obtain constants for integration of velocity (initial fore–aft velocity), we digitized a consistent point on the animal's head over a known distance (see Riskin et al., 2016). Once positional data were acquired, kinetic energy (EK) in the fore–aft (EKx), normal (EKz) and mediolateral (EKy) planes was calculated as:
formula
(1)
where m is the mass of the animal and v is the velocity of the COM in each respective axis. Potential energy of the COM (EP) was calculated using:
formula
(2)
where m is the mass of the animal, g is gravitational force (9.81 m s−1) and h is the vertical displacement in each plane COM in each respective axis. Total fore–aft, mediolateral and normal kinetic (EK) and potential energies (EP) were summed during each stride to obtain total energy (ET) of the COM. The power output of each individual component [power of kinetic energy in each of the fore–aft (PEKx), normal (PEKz) and mediolateral (PEKy) planes as well as total kinetic energy (PEK), total potential energy (PEP) and total energy (PET)] was determined by summing the positive increments over a stride and dividing by the duration of the stride.

Statistical analysis

Data on SLFs, mechanical power and spatiotemporal variables (contact time, stride time, duty factor, stride length and stride frequency) were analyzed in R (https://www.r-project.org/) using the packages ‘lmerTest’ (Kuznetsova et al., 2017) and ‘lme4’ (Bates et al., 2015). Normality of datasets was assessed with Shapiro–Wilk and Levene's tests (Sokal and Rohlf, 2012) before rank transforming all data to satisfy assumptions of subsequent statistical analyses (Sokal and Braumann, 1980). Speed, which has been suggested to influence the magnitude of force data (for statistical discussion, see Granatosky et al., 2020), was included in all linear mixed models discussed below. Additionally, individual behavioral idiosyncrasies were accounted for by including individual as a random effect following Winter (2013 preprint) and Bates et al. (2015) in all models.

Spatiotemporal gait characteristics

Slow walking gaits were used on all substrates (i.e. pole versus flat) and orientations (i.e. vertical versus horizontal). Both substrate and orientation had a significant (all P<0.001; Table 1, Fig. 2, see gait sequence in Fig. 3; prediction 1) influence on locomotor speed, duty factor, stride time, contact time and stride frequency. Stride length was significantly influenced by orientation (P<0.001), but not substrate (P=0.136) (see Fig. 2).

Fig. 2.

Bivariate plot demonstrating the relationship between stride frequency and speed (left) and stride length and speed (right). Horizontal (orange) and vertical (purple) orientations are depicted by color. Flat (circle) and pole (triangle) subtrates are denoted using different shapes. Dashed and solid lines represent regression lines for the pole and flat substrates, respectively. Regression in stride frequency for horizontal flat: y=0.05x+0.05, R2=0.56; horizontal pole: y=0.02x+0.03, R2=0.14; vertical flat: y=0.08x–0, R2=0.37; vertical pole: y=0.03x+0.04, R2=0.08. Regression in stride length for horizontal flat: y=0.83x+0.09, R2=0.09; horizontal pole: y=0.72x+0.02, R2=0.82; vertical flat: y=0.46x+0.05, R2=0.07; vertical pole: y=0.55x+0.02, R2=0.49.

Fig. 2.

Bivariate plot demonstrating the relationship between stride frequency and speed (left) and stride length and speed (right). Horizontal (orange) and vertical (purple) orientations are depicted by color. Flat (circle) and pole (triangle) subtrates are denoted using different shapes. Dashed and solid lines represent regression lines for the pole and flat substrates, respectively. Regression in stride frequency for horizontal flat: y=0.05x+0.05, R2=0.56; horizontal pole: y=0.02x+0.03, R2=0.14; vertical flat: y=0.08x–0, R2=0.37; vertical pole: y=0.03x+0.04, R2=0.08. Regression in stride length for horizontal flat: y=0.83x+0.09, R2=0.09; horizontal pole: y=0.72x+0.02, R2=0.82; vertical flat: y=0.46x+0.05, R2=0.07; vertical pole: y=0.55x+0.02, R2=0.49.

Fig. 3.

Average COM force, speed, COM position, energy profiles, and a representative gait sequence during a stride cycle in horizontal walking (left) and vertical climbing (right). A stride was defined from forelimb touchdown to subsequent forelimb touchdown. Solid lines represent the fore–aft direction (braking/propulsive), dashed lines represent mediolateral, and dotted lines represent normal (pulling and pushing in/from substrate) in force and speed traces. Energy components are represented as a solid line for potential energy (EP), dashed for kinetic energy (EK) and dotted for total energy (ET). LH, left hindlimb; RH, right hindlimb; LF, left forelimb; RF, right forelimb.

Fig. 3.

Average COM force, speed, COM position, energy profiles, and a representative gait sequence during a stride cycle in horizontal walking (left) and vertical climbing (right). A stride was defined from forelimb touchdown to subsequent forelimb touchdown. Solid lines represent the fore–aft direction (braking/propulsive), dashed lines represent mediolateral, and dotted lines represent normal (pulling and pushing in/from substrate) in force and speed traces. Energy components are represented as a solid line for potential energy (EP), dashed for kinetic energy (EK) and dotted for total energy (ET). LH, left hindlimb; RH, right hindlimb; LF, left forelimb; RF, right forelimb.

Table 1.

Spatiotemporal variables [duty factor, stride time, stride frequency and stride length (means±s.d.)] collected from five Australian tree frogs (Ranoidea caerulea) across each of the four experimental setups (flat and pole substrates oriented horizontally and vertically)

Spatiotemporal variables [duty factor, stride time, stride frequency and stride length (means±s.d.)] collected from five Australian tree frogs (Ranoidea caerulea) across each of the four experimental setups (flat and pole substrates oriented horizontally and vertically)
Spatiotemporal variables [duty factor, stride time, stride frequency and stride length (means±s.d.)] collected from five Australian tree frogs (Ranoidea caerulea) across each of the four experimental setups (flat and pole substrates oriented horizontally and vertically)

Single limb forces

Speed did not significantly impact impulses in any plane (Table 2) during either horizontal or vertical locomotion (all P>0.232), likely reflecting the limited range of locomotor speeds exhibited by the frogs. Frogs were found to produce significantly greater net fore–aft impulses on flat compared with the pole substrate (P=0.033; Tables 2 and 3, Fig. 4; prediction 2); however, inter-substrate effects were not observed in mediolateral and normal impulses. The hindlimbs of Australian tree frogs were found to produce larger net fore–aft and normal impulses compared with the forelimbs (all P≤0.001). Net fore–aft impulses were significantly greater during vertical climbing compared with horizontal walking (P<0.001; Tables 2, 3, Fig. 4). Conversely, net mediolateral and net normal impulses were significantly greater during horizontal walking (P<0.001; Tables 2, 3, Fig. 4), but did not differ between pole and flat conditions.

Fig. 4.

Impulse (% BWS) differences in the fore–aft (braking and propulsive), mediolateral and normal (pulling and pushing) axes from three Australian green tree frogs (Ranoidea caerulea) between horizontal (left of dashed line) and vertical (right of dashed line) experimental conditions. All forelimb impulses are denoted in blue and hindlimb boxes are in yellow. Substrate differences are denoted by circles for flat and triangles for pole.

Fig. 4.

Impulse (% BWS) differences in the fore–aft (braking and propulsive), mediolateral and normal (pulling and pushing) axes from three Australian green tree frogs (Ranoidea caerulea) between horizontal (left of dashed line) and vertical (right of dashed line) experimental conditions. All forelimb impulses are denoted in blue and hindlimb boxes are in yellow. Substrate differences are denoted by circles for flat and triangles for pole.

Table 2.

Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects (speed, orientation, limb, substrate) on net impulses in the fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes

Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects (speed, orientation, limb, substrate) on net impulses in the fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes
Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects (speed, orientation, limb, substrate) on net impulses in the fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes
Table 3.

Means±s.d. of impulses (% BWS) in fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes from five Australian green tree frogs (Ranoidea caerulea) during horizontal walking and vertical climbing on pole or flat substrates and stratified by forelimb versus hindlimb

Means±s.d. of impulses (% BWS) in fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes from five Australian green tree frogs (Ranoidea caerulea) during horizontal walking and vertical climbing on pole or flat substrates and stratified by forelimb versus hindlimb
Means±s.d. of impulses (% BWS) in fore–aft (braking/propulsive), normal (pull/push) and mediolateral axes from five Australian green tree frogs (Ranoidea caerulea) during horizontal walking and vertical climbing on pole or flat substrates and stratified by forelimb versus hindlimb

Center of mass dynamics

Compared with horizontal walking, the total mechanical energetic cost of vertical climbing was predominantly driven by potential energy, with negligible kinetic contributions (Fig. 3). Although phasic relationships between potential and kinetic energy profiles are present in some trials, these are absent during vertical climbing (Fig. 3). Total mechanical power was significantly greater during vertical climbing than during horizontal walking (P<0.001; Tables 4 and 5, Fig. 5A), driven by a major increase in potential power during vertical climbing. Kinetic energy in all planes (fore–aft, mediolateral and normal) was greater during vertical climbing than during horizontal walking (all P<0.001; Tables 4, 5, Fig. 5A).

Fig. 5.

Costs of horizontal walking and vertical climbing. (A) Mass-specific positive power (as calculated by the summation of positive increments in each component of energy over the stride divided by stride period) between horizontal (orange) and vertical (purple) orientations. Power of: kinetic energy in the fore–aft direction (PEKx), kinetic energy in the mediolateral direction (PEKy), kinetic energy in the normal direction (PEKz), total kinetic energy (PEK), total potential energy (PEP) and total positive mechanical power (PET). (B) Scatterplot of mechanical power as a function of speed. Orientation is differentiated by shape, with circles representing horizontal walking and triangles depicting vertical climbing. The red-legged running frog (Kassina maculata; Ahn et al., 2004) is represented in red, the Australian green tree frog (Ranoidea caerulea; present study) in light green and the gecko (Hemidactylus garnotii; Autumn et al., 2006; Chen et al., 2006) in dark green. Data from Ahn et al. (2004), Autumn et al. (2006) and Chen et al. (2006) were digitized using WebPlotDigitizer (Drevon et al., 2017; Rohatgi, 2017). The dashed line represents the minimum possible power production necessary to climb, as calculated by speed times gravity (∼9.81 m s−2).

Fig. 5.

Costs of horizontal walking and vertical climbing. (A) Mass-specific positive power (as calculated by the summation of positive increments in each component of energy over the stride divided by stride period) between horizontal (orange) and vertical (purple) orientations. Power of: kinetic energy in the fore–aft direction (PEKx), kinetic energy in the mediolateral direction (PEKy), kinetic energy in the normal direction (PEKz), total kinetic energy (PEK), total potential energy (PEP) and total positive mechanical power (PET). (B) Scatterplot of mechanical power as a function of speed. Orientation is differentiated by shape, with circles representing horizontal walking and triangles depicting vertical climbing. The red-legged running frog (Kassina maculata; Ahn et al., 2004) is represented in red, the Australian green tree frog (Ranoidea caerulea; present study) in light green and the gecko (Hemidactylus garnotii; Autumn et al., 2006; Chen et al., 2006) in dark green. Data from Ahn et al. (2004), Autumn et al. (2006) and Chen et al. (2006) were digitized using WebPlotDigitizer (Drevon et al., 2017; Rohatgi, 2017). The dashed line represents the minimum possible power production necessary to climb, as calculated by speed times gravity (∼9.81 m s−2).

Table 4.

Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects on potential, kinetic and total power

Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects on potential, kinetic and total power
Statistical parameters derived from least squares regressions demonstrating the statistical importance of various fixed effects on potential, kinetic and total power
Table 5.

Sum of positive mechanical power (calculated as the sum of the positive increments of power over a stride divided by the stride duration) represented as means±s.d., collected from three Australian tree frogs (Ranoidea caerulea) while navigating horizontally and vertically oriented flat substrates

Sum of positive mechanical power (calculated as the sum of the positive increments of power over a stride divided by the stride duration) represented as means±s.d., collected from three Australian tree frogs (Ranoidea caerulea) while navigating horizontally and vertically oriented flat substrates
Sum of positive mechanical power (calculated as the sum of the positive increments of power over a stride divided by the stride duration) represented as means±s.d., collected from three Australian tree frogs (Ranoidea caerulea) while navigating horizontally and vertically oriented flat substrates

Gait characteristics

The Australian green tree frogs in this study primarily adopted slow, walking (i.e. duty factors above 50%) gaits (Table 5). Consistent with studies of arboreal movement in other taxa (Granatosky et al., 2019b; Hanna et al., 2022; Herrel et al., 2013; Manzano et al., 2008; Nyakatura and Heymann, 2010), horizontal pole and vertical locomotion resulted in a decrease in locomotor speed and an increase in both duty factor and contact time compared with flat horizontal walking. Thus, prediction 1 was supported (Fig. 2).

Fore–aft forces

Fore–aft force profiles for Australian green tree frogs during horizontal locomotion largely match global patterns observed across tetrapods (Bertram, 2016; Granatosky et al., 2018, 2020). These include a braking to propulsive transition, and clear braking and propulsive functional differentiation between the forelimb and the hindlimb (Bertram, 2016; Granatosky et al., 2018), supporting prediction 2 (Fig. 4). Indeed, across all horizontal substrate conditions, Australian green tree frogs demonstrated a tendency for the forelimb to serve a net braking function whereas the hindlimb acts in a primarily propulsive role (Table 3, Fig. 4), a pattern consistent across non-belly-dragging quadrupeds (Nyakatura et al., 2014). Peak fore–aft magnitudes were generally similar in Australian green tree frogs during horizontal locomotion, once corrected for directionality and irrespective of substrate (i.e. flat or pole); however, on flat substrates, forelimb peak braking force was higher than any other fore–aft force magnitude (Table S1). Although not quantified, we noticed a tendency for the forelimb to ‘slap’ the substrate at touchdown, which may explain the higher braking forces. Otherwise, impulse patterns between the forelimbs and hindlimbs are generally similar across substrates. Specifically, the forelimb has relatively larger braking impulses, whereas the hindlimb shows greater propulsive impulses (Table 3, Fig. 4). Accordingly, fore–aft impulse is net braking in the forelimb and net propulsive in the hindlimb. Fore–aft impulse is generally greater while animals move on the horizontally oriented pole compared with the horizontal flat support (Table 3, Fig. 4). Such a finding can be attributed to the relatively longer contact times and slower velocities observed while Australian green tree frogs moved on simulated arboreal supports (Table 1).

In contrast, during vertical climbing, both the forelimb and hindlimb show exclusively propulsive forces (Table 3, Fig. 4; Table S1), a finding supported across climbing insects and tetrapods (Autumn et al., 2006; Cruse, 1976; Goldman et al., 2006; Wang et al., 2015) and in accordance with prediction 2.

Normal forces

Across tetrapods there is a general trend, with the exception of primates (e.g. Kimura, 1979; Reynolds, 1985) and arboreal opossums (e.g. Schmitt and Lemelin, 2002), to support most of the body weight on the forelimbs (Cieri et al., 2021; Fish et al., 2021; Gray, 1944). The Australian green tree frog is no exception to this pattern despite the substantial size differentiation between the smaller forelimbs and large hindlimbs. Although there was a tendency for more equal contributions of weight support between the limbs when Australian green tree frogs moved on the horizontally oriented pole (Table 3, Fig. 4), substrate effects in linear mixed-effect models were insignificant.

Australian green tree frogs follow generalized animal patterns of functionally differentiating the cranially positioned limb (pulling) and caudually positioned limb (pushing) (Autumn et al., 2006; Goldman et al., 2006; Hanna et al., 2017) when climbing. This pattern appears to be true irrespective of the clinging mechanism (Autumn et al., 2006; Clemente and Federle, 2008; Cruse, 1976; Hanna et al., 2017; Wang et al., 2015). To our knowledge, no study has yet reported limb-loading patterns during vertical climbing of a claw-clinging species (e.g. squirrels, felines, arboreal varanids, etc.), but data from Young et al. (2022a) showed that the claw-clinging hindlimb of parrots follows the same principles described above.

Hanna et al. (2017) reported that during ascent the forelimbs and hindlimbs of primates each experience a pull/push transition, meaning that as the forelimbs and hindlimbs touch down there is an initial pulling normal force subsequently followed by a pushing force at the end of stance phase. Such a pull/push transition during climbing is likely attributable to a limb moving from a cranial to caudal position relative to the COM (Granatosky et al., 2018; Hanna et al., 2017; Norberg, 1986). Although the hindlimb normal force profiles of Australian green tree frogs did indeed demonstrate a pull/push transition (Table 3, Fig. 4; Table S1), we observed no pushing forces in the forelimb. This pattern is consistent with cockroaches (Goldman et al., 2006), geckos (Autumn et al., 2006; Wang et al., 2015) and the beak and hindlimbs of climbing parrots (Young et al., 2022a). Such a taxonomic difference between primates and other climbing animals requires further investigation and may be attributable to several factors including the tendency for primates to demonstrate large limb excursions (Granatosky et al., 2019a; Larson et al., 2001), overstriding (Cartmill et al., 2020; Demes et al., 1994) or perhaps some unknown phylogenetic differentiation.

Mediolateral forces

The analysis and interpretation of single limb mediolateral forces has always been challenging in animal locomotion studies. Accordingly, such data is often not reported (Hanna et al., 2017; Hirasaki et al., 1993; Young et al., 2022a) and overall patterns remain elusive for climbing. Mediolateral forces were significantly (P<0.001) influenced by substrate orientation, but not substrate type (i.e. pole versus flat), contrary to prediction 2. During horizontal walking on a flat substrate, net forelimb mediolateral impulses are approximately equally distributed between medial and lateral directions, but, when traversing a vertical flat substrate, forelimb impulses resolved as medially oriented (Fig. 4). A similar pattern is also observed in the hindlimb (Fig. 4). Mediolateral forces were also highly variable and low in magnitude compared with the normal and fore–aft force components. Thus, perhaps it is incorrect to think about mediolateral force patterns as some preordained motor ‘goal’, and such forces should instead be considered as step-to-step adjustments for maintaining balance throughout a stride (Bishop et al., 2018).

Center of mass mechanics

Although COM mechanics during horizontal walking have been widely discussed (Ahn et al., 2004; Blickhan and Full, 1987; Cavagna et al., 1977; Heglund et al., 1982; Reilly et al., 2006; Rubenson et al., 2004), COM mechanics for vertical climbing remain poorly studied (Autumn et al., 2006; Goldman et al., 2006). During horizontal walking, a phasic relationship between kinetic and potential energy drives energy recovery. However, during climbing, the phasic relationship between kinetic and potential energy is completely uncoupled owing to the continuous accumulation of potential energy associated with gaining height (prediction 3; Fig. 3). Instead, kinetic components of energy should drive efficiency (Autumn et al., 2006). It can be inferred that efficient climbers should minimize oscillations in their mediolateral and normal axes, thereby reducing total kinetic energetic costs to the bare minimum of fore–aft kinetic energy, lying in their sole axis of travel (Autumn et al., 2006; Goldman et al., 2006; prediction 3). However, this was not observed within our data; instead, both normal and mediolateral kinetic energies were greater during vertical climbing than during horizontal walking (Table 3).

One potential explanation for this trend is that climbing frogs do not utilize particularly efficient locomotion. We can assess climbing efficiency by comparing the total external mechanical power exerted during ascent with the minimum power requirement to lift one's COM against the force of gravity (i.e. the product of gravity and speed; Autumn et al., 2006). Comparing this metric of efficiency between our data and climbing geckos published by Autumn et al. (2006), we observe that frogs are closer to this theoretical minimum of power production than geckos [residual sum of squares (RSS) of frogs=0.023; RSS of geckos=3.22; see Fig. 5]. Therefore, this lack of concordance between our data and theory cannot be ascribed to inefficient climbing.

Instead, we must consider other possible mechanisms to explain efficient climbing. One such explanation is that out-of-plane oscillations of the COM may be less deleterious than originally theorized. In support of this idea, it is noteworthy that mediolateral oscillations of the COM are reported as typical among expert human rock climbers, though no mechanistic explanation has been proposed for these movements (Saul et al., 2019; Zampagni et al., 2011). It is plausible that the pendular nature of left–right swinging in an axis perpendicular to that of the direction of travel may serve as analogous to normal oscillations of the COM during horizontal walking, allowing animals to recover greater energy while climbing. This theory, however, must be further explored through close examination of kinetic energy profiles in future studies.

Finally, we would be remiss to discount the fact that Australian green tree frogs may simply be an inappropriate model to assess mechanical energy between horizontal and vertical locomotion, owing to their small body size and the start–stop nature of their locomotion (Hooper, 2012; Reilly et al., 2007). This is because energy-saving mechanisms associated with greater pendular oscillations of the COM would reduce the relative cost of horizontal walking at large body sizes (Hooper, 2012; Reilly et al., 2007). Such mechanisms are not available during climbing. This disparity between the relative cost of horizontal locomotion and vertical climbing as body size increases has been demonstrated based on metabolic energetic expenditure (Hanna and Schmitt, 2011; Hanna et al., 2008; Taylor et al., 1972). Taken together, we predict that it is at large body sizes that true mechanical energy savings during climbing, most likely driven by modulation of kinetic energy (Autumn et al., 2006), may come into play and provide the necessary variation for natural selection to act upon and drive changes in anatomy and behavior in arboreal species. It is also possible that there is a critical mass at which the relative mechanical costs of climbing no longer represent an energetically feasible form of locomotion (Cant, 1992). However, the presence of climbing behaviors in large-bodied taxa such as bears (Amanat et al., 2020) and gorillas (Remis, 1998) brings such speculation into question. To address these possibilities would require a sample across a broad range of body sizes to assess allometry of climbing mechanical costs.

Conclusions

During horizontal walking, Australian green tree frogs demonstrate the typical quadruped pattern wherein the forelimb serves a net braking function while the hindlimb acts in a primarily propulsive role. Conversely, during vertical climbing both the forelimb and hindlimb show exclusively propulsive forces. In the normal plane, the forelimb serves a net pulling function while the hindlimb is net pushing (Table 2, Fig. 2), in line with observations in other climbing taxa. Finally, mediolateral forces remain enigmatic to interpret and may reflect adaptations towards maintaining stability during a stride (Bishop et al., 2018).

This study is one of the few (Autumn et al., 2006; Goldman et al., 2006) to report COM movements and energy profiles during non-human climbing. Overall, Australian green tree frogs match theoretical predictions of climbing dynamics (i.e. total energetic requirements are primarily driven by the accumulation of potential energy as the animal ascends; Autumn et al., 2006; Goldman et al., 2006). We suggest that attention should be directed towards fluctuations in kinetic energy, the only means by which natural selection can act to reduce overall mechanical costs. Utilizing power as a means for estimating efficiency, we observed that Australian green tree frogs are highly mechanically efficient, barely exceeding the minimum mechanical power necessary to climb at a certain speed (Fig. 5). Although it is often assumed that gaits maximizing stability are in direct conflict with speed-related decreases in overall cost of transport (Heglund and Taylor, 1988; Heglund et al., 1982; O'Neill, 2012), it would appear that Australian green tree frogs move both slowly and efficiently. Indeed, a tendency for careful movement during vertical climbing seems to be a broad pattern observed across climbing species (Goldman et al., 2006; Granatosky et al., 2019b; Hanna et al., 2022; Herrel et al., 2013; Karantanis et al., 2018), and may represent a viable means to achieve both stability and mechanical efficiency. All climbing species are not slow, however, and assessing how such species balance stability and energetic considerations would be of great interest. The findings in this study raise exciting considerations about what makes a good climber from a mechanical perspective, and future work is required to quantify variation in climbing efficiency across species of differing anatomical configurations and body sizes.

We are grateful to two anonymous reviewers for their thoughtful and constructive feedback on an earlier version of the paper. We thank Jeanne Quidore-Jermann and Sandra Kahler, and the Core Facilities team at NYITCOM's animal care facility, for their help with animal housing and care. We also thank Daniel Tanis for initial animal procurement and care, Aleksandra Ratkiewicz for assistance in figure generation, and Todd Green for photographing our experimental animals.

Author contributions

Conceptualization: M.W.Y., N.D.F., E.D., M.C.G.; Methodology: M.W.Y., N.D.F., J.R.Y., A.R., N.G., M.C.G.; Software: M.W.Y., N.D.F., M.C.G.; Validation: M.W.Y., N.D.F., J.R.Y., E.D., M.C.G.; Formal analysis: M.W.Y., N.D.F., J.R.Y., A.R., N.G., M.C.G.; Investigation: M.W.Y., N.D.F., J.R.Y., A.R., N.G., E.D., M.C.G.; Resources: M.C.G.; Data curation: M.W.Y., J.R.Y., A.R., N.G., E.D., M.C.G.; Writing - original draft: M.W.Y., J.R.Y., A.R., N.G., E.D., M.C.G.; Writing - review & editing: M.W.Y., J.R.Y., A.R., N.G., E.D., M.C.G.; Visualization: M.W.Y., N.D.F., E.D., M.C.G.; Supervision: M.C.G.; Project administration: M.W.Y., M.C.G.; Funding acquisition: M.C.G.

Funding

This study was funded in part by the Center for Biomedical Innovation at New York Institute of Technology.

Data availability

All data used for statistical analysis within this study are available from GitHub: https://github.com/melodyoungg/TreeFrogs

Ahn
,
A. N.
,
Furrow
,
E.
and
Biewener
,
A. A.
(
2004
).
Walking and running in the red-legged running frog, Kassina maculata
.
J. Exp. Biol.
207
,
399
-
410
.
Amanat
,
S.
,
Mayer
,
J.
,
Paracha
,
H.
,
Ali
,
Z.
and
Granatosky
,
M. C.
(
2020
).
Bear locomotion
. In
Encyclopedia of Animal Cognition and Behavior
(ed.
J.
Vonk
and
T.
Shackelford
), pp.
1
-
6
.
Cham
:
Springer International Publishing
.
Autumn
,
K.
,
Hsieh
,
S. T.
,
Dudek
,
D. M.
,
Chen
,
J.
,
Chitaphan
,
C.
and
Full
,
R. J.
(
2006
).
Dynamics of geckos running vertically
.
J. Exp. Biol.
209
,
260
-
272
.
Bates
,
D.
,
Mächler
,
M.
,
Bolker
,
B.
and
Walker
,
S
. (
2015
).
Fitting linear mixed-effects models using lme4
.
J. Stat. Soft.
67
,
1-48
. doi:10.18637/jss.v067.i01
Beck
,
H. K.
,
Schultz
,
J. T.
and
Clemente
,
C. J.
(
2021
).
A bio-inspired robotic climbing robot to understand kinematic and morphological determinants for an optimal climbing gait
.
Bioinspir. Biomim.
17
,
016005
.
Bertram
,
J. E.
(
2016
).
Concepts in locomotion: levers, struts, pendula and springs
. In
Understanding Mammalian Locomotion: Concepts and Applications
, pp.
79
-
110
.
Hoboken, NJ
:
John Wiley & Sons, Inc
.
Bishop
,
K. L.
,
Pai
,
A. K.
and
Schmitt
,
D.
(
2008
).
Whole body mechanics of stealthy walking in cats
.
PLoS One
3
,
e3808
.
Bishop
,
P. J.
,
Graham
,
D. F.
,
Lamas
,
L. P.
,
Hutchinson
,
J. R.
,
Rubenson
,
J.
,
Hancock
,
J. A.
,
Wilson
,
R. S.
,
Hocknull
,
S. A.
,
Barrett
,
R. S.
,
Lloyd
,
D. G.
et al. 
(
2018
).
The influence of speed and size on avian terrestrial locomotor biomechanics: predicting locomotion in extinct theropod dinosaurs
.
PLoS One
13
,
e0192172
.
Blickhan
,
R.
,
Full
,
R.
(
1987
).
Locomotion energetics of the ghost crab: II. Mechanics of the centre of mass during walking and running
.
J. Exp. Biol.
130
,
155
-
174
.
Bock
,
W. J.
and
Winkler
,
H.
(
1978
).
Mechanical analysis of the external forces on climbing mammals
.
Zoomorphologie
91
,
49
-
61
.
Byron
,
C. D.
,
Granatosky
,
M. C.
and
Covert
,
H. H.
(
2017
).
An anatomical and mechanical analysis of the douc monkey (genus Pygathrix), and its role in understanding the evolution of brachiation
.
Am. J. Phys. Anthropol.
164
,
801
-
820
.
Cant
,
J. G. H.
(
1992
).
Positional behavior and body size of arboreal primates: a theoretical framework for field studies and an illustration of its application
.
Am. J. Phys. Anthropol.
88
,
273
-
283
.
Carlson
,
K. J.
,
Demes
,
B.
and
Franz
,
T. M.
(
2005
).
Mediolateral forces associated with quadrupedal gaits of lemurids
.
J. Zool.
266
,
261
-
273
.
Cartmill
,
M.
(
1985
).
Climbing
. In
Functional Vertebrate Morphology
(ed.
M.
Hildebrand
,
D. M.
Bramble
,
K. F.
Liem
and
D. B.
Wake
), pp.
73
-
88
.
Cambridge, MA
:
Belknap
.
Cartmill
,
M.
,
Brown
,
K.
,
Atkinson
,
C.
,
Cartmill
,
E. A.
,
Findley
,
E.
,
Gonzalez-Socoloske
,
D.
,
Hartstone-Rose
,
A.
and
Mueller
,
J.
(
2020
).
The gaits of marsupials and the evolution of diagonal-sequence walking in primates
.
Am. J. Phys. Anthropol.
171
,
182
-
197
.
Cavagna
,
G. A.
,
Heglund
,
N. C.
and
Taylor
,
C. R.
(
1977
).
Mechanical work in terrestrial locomotion: two basic mechanisms for minimizing energy expenditure
.
Am. J. Physiol.
233
,
R243
-
R261
.
Chan
,
K. O.
,
Grismer
,
L. L.
and
Brown
,
R. M.
(
2018
).
Comprehensive multi-locus phylogeny of Old World tree frogs (Anura: Rhacophoridae) reveals taxonomic uncertainties and potential cases of over- and underestimation of species diversity
.
Mol. Phylogenet. Evol.
127
,
1010
-
1019
.
Chen
,
J. J.
,
Peattie
,
A. M.
,
Autumn
,
K.
and
Full
,
R. J.
(
2006
).
Differential leg function in a sprawled-posture quadrupedal trotter
.
J. Exp. Biol.
209
,
249
-
259
.
Cieri
,
R. L.
,
Dick
,
T. J. M.
,
Irwin
,
R.
,
Rumsey
,
D.
and
Clemente
,
C. J.
(
2021
).
The scaling of ground reaction forces and duty factor in monitor lizards: implications for locomotion in sprawling tetrapods
.
Biol. Lett.
17
,
20200612
.
Clemente
,
C. J.
and
Federle
,
W.
(
2008
).
Pushing versus pulling: division of labour between tarsal attachment pads in cockroaches
.
Proc. R. Soc. B Biol. Sci.
275
,
1329
-
1336
.
Cole
,
M. H.
,
Van Den Hoorn
,
W.
,
Kavanagh
,
J. K.
,
Morrison
,
S.
,
Hodges
,
P. W.
,
Smeathers
,
J. E.
and
Kerr
,
G. K.
(
2014
).
Concurrent validity of accelerations measured using a tri-axial inertial measurement unit while walking on firm, compliant and uneven surfaces
.
PLoS One
9
,
e98395
.
Cruse
,
H.
(
1976
).
The function of the legs in the free walking stick insect, Carausius morosus
.
J. Comp. Physiol.
112
,
235
-
262
.
Demes
,
B.
,
Larson
,
S. G.
,
Stern
,
J. T.
,
Jungers
,
W. L.
,
Biknevicius
,
A. R.
and
Schmitt
,
D.
(
1994
).
The kinetics of primate quadrupedalism: ‘hindlimb drive’ reconsidered
.
J. Hum. Evol.
26
,
353
-
374
.
Demes
,
B.
,
Carlson
,
K. J.
and
Franz
,
T. M.
(
2006
).
Cutting corners: the dynamics of turning behaviors in two primate species
.
J. Exp. Biol.
209
,
927
-
937
.
Dickinson
,
E.
,
Young
,
M. W.
and
Granatosky
,
M. C.
(
2022
).
Testing mechanisms for weight support distribution during inverted quadrupedalism in primates
.
J. Exp. Zool. Part Ecol. Integr. Physiol.
337
,
699
-
708
.
Drevon
,
D.
,
Fursa
,
S. R.
and
Malcolm
,
A. L.
(
2017
).
Intercoder reliability and validity of WebPlotDigitizer in extracting graphed data
.
Behav. Modif.
41
,
323
-
339
.
Endlein
,
T.
,
Ji
,
A.
,
Yuan
,
S.
,
Hill
,
I.
,
Wang
,
H.
,
Barnes
,
W. J. P.
,
Dai
,
Z.
and
Sitti
,
M.
(
2017
).
The use of clamping grips and friction pads by tree frogs for climbing curved surfaces
.
Proc. R. Soc. B Biol. Sci.
284
,
20162867
.
Faltings
,
L.
,
Young
,
M. W.
,
Ross
,
C. F.
and
Granatosky
,
M. C.
(
2022
).
Got rhythm? Rhythmicity differences reflect different optimality criteria in feeding and locomotor systems
.
Evolution
76
,
2181
-
2190
.
Federle
,
W.
,
Barnes
,
W. J. P.
,
Baumgartner
,
W.
,
Drechsler
,
P.
and
Smith
,
J. M.
(
2006
).
Wet but not slippery: boundary friction in tree frog adhesive toe pads
.
J. R. Soc. Interface
3
,
689
-
697
.
Fischer
,
M. S.
,
Krause
,
C.
and
Lilje
,
K. E.
(
2010
).
Evolution of chameleon locomotion, or how to become arboreal as a reptile
.
Zoology
113
,
67
-
74
.
Fish
,
F. E.
,
Sheehan
,
M. J.
,
Adams
,
D. S.
,
Tennett
,
K. A.
and
Gough
,
W. T.
(
2021
).
A 60:40 split: differential mass support in dogs
.
Anat. Rec.
304
,
78
-
89
.
Franz
,
T. M.
,
Demes
,
B.
and
Carlson
,
K. J.
(
2005
).
Gait mechanics of lemurid primates on terrestrial and arboreal substrates
.
J. Hum. Evol.
48
,
199
-
217
.
Goldman
,
D. I.
,
Chen
,
T. S.
,
Dudek
,
D. M.
and
Full
,
R. J.
(
2006
).
Dynamics of rapid vertical climbing in cockroaches reveals a template
.
J. Exp. Biol.
209
,
2990
-
3000
.
Granatosky
,
M.
and
Krysko
,
K.
(
2011
).
Ontogenetic behavioral shifts in habitat utilization of treefrogs (Hylidae) in north-central Florida
.
IRCF Reptil. Amphib.
18
,
195
-
201
.
Granatosky
,
M. C.
and
Mcelroy
,
E. J.
(
2022
).
Stride frequency or length? A phylogenetic approach to understand how animals regulate locomotor speed
.
J. Exp. Biol.
225
,
jeb243231
.
Granatosky
,
M. C.
,
Fitzsimons
,
A.
,
Zeininger
,
A.
and
Schmitt
,
D.
(
2018
).
Mechanisms for the functional differentiation of the propulsive and braking roles of the forelimbs and hindlimbs during quadrupedal walking in primates and felines
.
J. Exp. Biol.
221
,
jeb162917
.
Granatosky
,
M. C.
,
Mcelroy
,
E. J.
,
Laird
,
M. F.
,
Iriarte-Diaz
,
J.
,
Reilly
,
S. M.
,
Taylor
,
A. B.
and
Ross
,
C. F.
(
2019a
).
Joint angular excursions during cyclical behaviors differ between tetrapod feeding and locomotor systems
.
J. Exp. Biol.
222
,
jeb200451
.
Granatosky
,
M. C.
,
Schmitt
,
D.
and
Hanna
,
J.
(
2019b
).
Comparison of spatiotemporal gait characteristics between vertical climbing and horizontal walking in primates
.
J. Exp. Biol.
222
,
jeb185702
.
Granatosky
,
M. C.
,
Mcelroy
,
E. J.
,
Lemelin
,
P.
,
Reilly
,
S. M.
,
Nyakatura
,
J. A.
,
Andrada
,
E.
,
Kilbourne
,
B. M.
,
Allen
,
V. R.
,
Butcher
,
M. T.
,
Blob
,
R. W.
et al. 
(
2020
).
Variation in limb loading magnitude and timing in tetrapods
.
J. Exp. Biol.
223
,
jeb201525
.
Granatosky
,
M. C.
,
Young
,
M. W.
,
Herr
,
V.
,
Chai
,
C.
,
Raidah
,
A.
,
Kairo
,
J. N.
,
Anaekwe
,
A.
,
Havens
,
A.
,
Zou
,
B.
,
Ding
,
B.
et al. 
(
2022
).
Positional behavior of introduced monk parakeets (Myiopsitta monachus) in an urban landscape
.
Animals
12
,
2372
.
Grand
,
T. I.
(
1972
).
A mechanical interpretation of terminal branch feeding
.
J. Mammal.
53
,
198
.
Gray
,
J.
(
1944
).
Studies in the mechanics of the tetrapod skeleton
.
J. Exp. Biol.
20
,
88
-
116
.
Gutierrez-Farewik
,
E. M.
,
Bartonek
,
Å.
and
Saraste
,
H.
(
2006
).
Comparison and evaluation of two common methods to measure center of mass displacement in three dimensions during gait
.
Hum. Mov. Sci.
25
,
238
-
256
.
Hanna
,
J. B.
and
Schmitt
,
D.
(
2011
).
Locomotor energetics in primates: gait mechanics and their relationship to the energetics of vertical and horizontal locomotion
.
Am. J. Phys. Anthropol.
145
,
43
-
54
.
Hanna
,
J. B.
,
Schmitt
,
D.
and
Griffin
,
T. M.
(
2008
).
The energetic cost of climbing in primates
.
Science
320
,
898
.
Hanna
,
J. B.
,
Granatosky
,
M. C.
,
Rana
,
P.
and
Schmitt
,
D.
(
2017
).
The evolution of vertical climbing in primates: evidence from reaction forces
.
J. Exp. Biol.
220
,
3039
-
3052
.
Hanna
,
C. S.
,
Alihosseini
,
C.
,
Fischer
,
H. M.
,
Davoli
,
E. C.
and
Granatosky
,
M. C.
(
2022
).
Are they arboreal? Climbing abilities and mechanics in the red-backed salamander (Plethodon cinereus)
.
J. Exp. Zool. Ecol. Integr. Physiol.
337
,
238
-
249
.
Heglund
,
N. C.
and
Taylor
,
C. R.
(
1988
).
Speed, stride frequency and energy cost per stride: how do they change with body size and gait?
J. Exp. Biol.
138
,
301
-
318
.
Heglund
,
N. C.
,
Fedak
,
M. A.
,
Taylor
,
C. R.
and
Cavagna
,
G. A.
(
1982
).
Energetics and mechanics of terrestrial locomotion. IV. Total mechanical energy changes as a function of speed and body size in birds and mammals
.
J. Exp. Biol.
97
,
57
-
66
.
Herrel
,
A.
,
Perrenoud
,
M.
,
Decamps
,
T.
,
Abdala
,
V.
,
Manzano
,
A.
and
Pouydebat
,
E.
(
2013
).
The effect of substrate diameter and incline on locomotion in an arboreal frog
.
J. Exp. Biol.
216
,
3599
-
3605
.
Higham
,
T. E.
and
Jayne
,
B. C.
(
2004
).
In vivo muscle activity in the hindlimb of the arboreal lizard, Chamaeleo calyptratus: general patterns and the effects of incline
.
J. Exp. Biol.
207
,
249
-
261
.
Hill
,
I. D. C.
,
Dong
,
B.
,
Barnes
,
W.
,
Jon
,
P.
,
Ji
,
A.
and
Endlein
,
T.
(
2018
).
The biomechanics of tree frogs climbing curved surfaces: a gripping problem
.
J. Exp. Biol.
221
,
jeb168179
.
Hirasaki
,
E.
,
Kumakura
,
H.
and
Matano
,
S.
(
1993
).
Kinesiological characteristics of vertical climbing in Ateles geoffroyi and Macaca fuscata
.
Folia Primatol.
61
,
148
-
156
.
Hooper
,
S. L.
(
2012
).
Body size and the neural control of movement
.
Curr. Biol.
22
,
R318
-
R322
.
Hunt
,
K. D.
,
Cant
,
J. G. H.
,
Gebo
,
D. L.
,
Rose
,
M. D.
,
Walker
,
S. E.
and
Youlatos
,
D.
(
1996
).
Standardized descriptions of primate locomotor and postural modes
.
Primates
37
,
363
-
387
.
Ihlen
,
E. A. F.
,
Skjæret
,
N.
and
Vereijken
,
B.
(
2013
).
The influence of center-of-mass movements on the variation in the structure of human postural sway
.
J. Biomech.
46
,
484
-
490
.
Jenkins
,
F. A.
(
1974
).
Tree shrew locomotion and the origins of primate arborealism
. In
Primate Locomotion
(ed.
F. A.
Jenkins
), pp.
85
-
115
.
New York
:
Academic Press
.
Kamada
,
K.
,
Tachibanagi
,
R.
and
Nakagawa
,
H.
(
2017
).
Strategy of landing behavior of the tree frog Hyla japonica
.
J. Exp. Zool. Ecol. Integr. Physiol.
327
,
600
-
610
.
Karantanis
,
N.-E.
,
Rychlik
,
L.
,
Herrel
,
A.
and
Youlatos
,
D.
(
2018
).
Vertical locomotion in Micromys minutus (Rodentia: Muridae): insights into the evolution of eutherian climbing
.
J. Mamm. Evol
.
25
,
277
-
289
.
Kasuya
,
E.
,
Hirota
,
M.
and
Shigehara
,
H.
(
1996
).
Reproductive behavior of the Japanese treefrog, Rhacophorus arboreus (Anura: Rhacophoridae)
.
Popul. Ecol.
38
,
1
-
10
.
Kimura
,
T.
(
1979
).
Kinesiological characteristics of primate walking
. In
Environment, Behavior, and Morphology: Dynamic Interactions in Primates
(ed.
H
Morbeck
,
H.
Preuschoft
and
N.
Gomberg
), pp.
297
-
311
.
New York, NY
:
G. Fischer
.
Kuznetsova
,
A.
,
Brockhoff
,
P. B.
and
Christensen
,
R. HB
. (
2017
).
lmerTest package: tests in linear mixed effects models
.
J. Stat. Softw.
82
,
1
-
26
.
Labonte
,
D.
and
Federle
,
W.
(
2015
).
Scaling and biomechanics of surface attachment in climbing animals
.
Philos. Trans. R. Soc. B Biol. Sci.
370
,
20140027
.
Labonte
,
D.
and
Federle
,
W.
(
2016
).
Biomechanics of shear-sensitive adhesion in climbing animals: peeling, pre-tension and sliding-induced changes in interface strength
.
J. R. Soc. Interface
13
,
20160373
.
Larson
,
S. G.
(
1998
).
Unique aspects of quadrupedal locomotion in nonhuman primates
. In
Primate Locomotion
(ed.
E.
Strasser
,
J. G.
Fleagle
,
A. L.
Rosenberger
and
H. M.
McHenry
), pp.
157
-
173
.
Boston, MA
:
Springer US
.
Larson
,
S. G.
and
Demes
,
B.
(
2011
).
Weight support distribution during quadrupedal walking in Ateles and Cebus
.
Am. J. Phys. Anthropol.
144
,
633
-
642
.
Larson
,
S. G.
,
Schmitt
,
D.
,
Lemelin
,
P.
and
Hamrick
,
M.
(
2001
).
Limb excursion during quadrupedal walking: how do primates compare to other mammals?
J. Zool.
255
,
353
-
365
.
Lemelin
,
P.
and
Cartmill
,
M.
(
2010
).
The effect of substrate size on the locomotion and gait patterns of the kinkajou (Potos flavus)
.
J. Exp. Zool. A Ecol. Genet. Physiol.
313
,
157
-
168
.
Lemelin
,
P.
,
Schmitt
,
D.
and
Cartmill
,
M.
(
2003
).
Footfall patterns and interlimb co-ordination in opossums (Family Didelphidae): evidence for the evolution of diagonal-sequence walking gaits in primates
.
J. Zool.
260
,
423
-
429
.
Manzano
,
A. S.
,
Abdala
,
V.
and
Herrel
,
A.
(
2008
).
Morphology and function of the forelimb in arboreal frogs: specializations for grasping ability?
J. Anat.
213
,
296
-
307
.
Norberg
,
R. Å.
(
1986
).
Treecreeper climbing; mechanics, energetics, and structural adaptations
.
Scand. J. Ornithol.
17
,
191
-
209
.
Nyakatura
,
J. A.
and
Heymann
,
E. W.
(
2010
).
Effects of support size and orientation on symmetric gaits in free-ranging tamarins of Amazonian Peru: implications for the functional significance of primate gait sequence patterns
.
J. Hum. Evol.
58
,
242
-
251
.
Nyakatura
,
J. A.
,
Andrada
,
E.
,
Curth
,
S.
and
Fischer
,
M. S.
(
2014
).
Bridging ‘Romer's gap’: limb mechanics of an extant belly-dragging lizard inform debate on tetrapod locomotion during the early carboniferous
.
Evol. Biol.
41
,
175
-
190
.
Ogihara
,
N.
,
Oku
,
T.
,
Andrada
,
E.
,
Blickhan
,
R.
,
Nyakatura
,
J. A.
and
Fischer
,
M. S.
(
2014
).
Planar covariation of limb elevation angles during bipedal locomotion in common quails (Coturnix coturnix)
.
J. Exp. Biol.
217
,
3968
-
3973
.
O'Neill
,
M. C.
(
2012
).
Gait-specific metabolic costs and preferred speeds in ring-tailed lemurs (Lemur catta), with implications for the scaling of locomotor costs
.
Am. J. Phys. Anthropol.
149
,
356
-
364
.
O'Neill
,
M. C.
and
Schmitt
,
D.
(
2012
).
The gaits of primates: center of mass mechanics in walking, cantering and galloping ring-tailed lemurs, Lemur catta
.
J. Exp. Biol.
215
,
1728
-
1739
.
Peloso
,
P. L. V.
,
Oliveira
,
R. M. D.
,
Sturaro
,
M. J.
,
Rodrigues
,
M. T.
,
Lima-Filho
,
G. R.
,
Bitar
,
Y. O. C.
,
Wheeler
,
W. C.
and
Aleixo
,
A.
(
2018
).
Phylogeny of map tree frogs, Boana semilineata species group, with a new Amazonian species (Anura: Hylidae)
.
South Am. J. Herpetol.
13
,
150
-
169
.
Preuschoft
,
H.
(
2002
).
What does ‘arboreal locomotion’ mean exactly and what are the relationships between ‘climbing’, environment and morphology?
Z. Für Morphol. Anthropol
83
,
171
-
188
.
Proske
,
U.
(
1979
).
The Golgi tendon organ
.
Trends Neurosci.
2
,
7
-
8
.
Proske
,
U.
(
1997
).
The mammalian muscle spindle
.
Physiology
12
,
37
-
42
.
Raichlen
,
D. A.
,
Pontzer
,
H.
,
Shapiro
,
L. J.
and
Sockol
,
M. D.
(
2009
).
Understanding hind limb weight support in chimpanzees with implications for the evolution of primate locomotion
.
Am. J. Phys. Anthropol.
138
,
395
-
402
.
Reader
,
L. L.
,
Carrier
,
D. R.
,
Goller
,
F.
,
Isaacs
,
M. R.
,
Moore Crisp
,
A.
,
Barnes
,
C. J.
and
Lee
,
D. V.
(
2022
).
Climbing parrots achieve pitch stability using forces and free moments produced by axial-appendicular couples
.
J. Exp. Biol.
225
,
jeb242305
.
Reilly
,
S. M.
,
Mcelroy
,
E. J.
,
Odum
,
R. A.
and
Hornyak
,
V. A.
(
2006
).
Tuataras and salamanders show that walking and running mechanics are ancient features of tetrapod locomotion
.
Proc. R. Soc. Lond. B Biol. Sci.
273
,
1563
-
1568
.
Reilly
,
S. M.
,
Mcelroy
,
E. J.
and
Biknevicius
,
A. R.
(
2007
).
Posture, gait and the ecological relevance of locomotor costs and energy-saving mechanisms in tetrapods
.
Zoology
110
,
271
-
289
.
Remis
,
M. J.
(
1998
).
The gorilla paradox
. In
Primate Locomotion: Recent Advances
(ed.
E.
Strasser
,
J. G.
Fleagle
,
A. L.
Rosenberger
and
H. M.
McHenry
), pp.
95
-
106
.
Boston, MA
:
Springer US
.
Reynolds
,
T.
(
1985
).
Stresses on the limbs of quadrupedal primates
.
Am. J. Phys. Anthropol.
67
,
351
-
362
.
Richards
,
C. M.
and
Moore
,
W. S.
(
1996
).
A phylogeny for the African treefrog family hyperoliidae based on mitochondrial rDNA
.
Mol. Phylogenet. Evol.
5
,
522
-
532
.
Riskin
,
D. K.
,
Kendall
,
C. J.
and
Hermanson
,
J. W.
(
2016
).
The crouching of the shrew: mechanical consequences of limb posture in small mammals
.
PeerJ
4
,
e2131
.
Rohatgi
,
A.
(
2017
).
WebPlotDigitizer
.
Austin, TX
,
USA
.
Ross
,
C. F.
,
Blob
,
R. W.
,
Carrier
,
D. R.
,
Daley
,
M. A.
,
Deban
,
S. M.
,
Demes
,
B.
,
Gripper
,
J. L.
,
Iriarte-Diaz
,
J.
,
Kilbourne
,
B. M.
,
Landberg
,
T.
et al. 
(
2013
).
The evolution of locomotor rhythmicity in tetrapods
.
Evolution
67
,
1209
-
1217
.
Rubenson
,
J.
,
Heliams
,
D. B.
,
Lloyd
,
D. G.
and
Fournier
,
P. A.
(
2004
).
Gait selection in the ostrich: mechanical and metabolic characteristics of walking and running with and without an aerial phase
.
Proc. R. Soc. B Biol. Sci.
271
,
1091
-
1099
.
Saul
,
D.
,
Steinmetz
,
G.
,
Lehmann
,
W.
and
Schilling
,
A. F.
(
2019
).
Determinants for success in climbing: a systematic review
.
J. Exerc. Sci. Fit.
17
,
91
-
100
.
Schmitt
,
D.
(
1999
).
Compliant walking in primates
.
J. Zool.
248
,
149
-
160
.
Schmitt
,
D.
(
2003
).
Mediolateral reaction forces and forelimb anatomy in quadrupedal primates: implications for interpreting locomotor behavior in fossil primates
.
J. Hum. Evol.
44
,
47
-
58
.
Schmitt
,
D.
and
Hanna
,
J. B.
(
2004
).
Substrate alters forelimb to hindlimb peak force ratios in primates
.
J. Hum. Evol.
46
,
239
-
254
.
Schmitt
,
D.
and
Lemelin
,
P.
(
2002
).
Origins of primate locomotion: gait mechanics of the woolly opossum
.
Am. J. Phys. Anthropol.
118
,
231
-
238
.
Schneider
,
C. A.
,
Rasband
,
W. S.
and
Eliceiri
,
K. W.
(
2012
).
NIH Image to ImageJ: 25 years of image analysis
.
Nat. Methods
9
,
671
-
675
.
Schultz
,
J. T.
,
Beck
,
H. K.
,
Haagensen
,
T.
,
Proost
,
T.
and
Clemente
,
C. J.
(
2021
).
Using a biologically mimicking climbing robot to explore the performance landscape of climbing in lizards
.
Proc. R. Soc. B Biol. Sci.
288
,
20202576
.
Sokal
,
R. R.
and
Braumann
,
C. A
. (
1980
).
Significance tests for coefficients of variation and variability profiles
.
Syst. Biol.
29
,
50
-
66
.
Sokal
,
R. R.
and
Rohlf
,
F. J.
(
2012
).
Biometry: The Principles and Practice of Statistics in Biological Research
, 6th edn.
New York
:
Freeman & Company
.
Spring
,
L. W.
(
1965
).
Climbing and pecking adaptations in some north american woodpeckers
.
Condor
67
,
457
-
488
.
Taylor
,
C. R.
,
Caldwell
,
S. L.
and
Rowntree
,
V. J.
(
1972
).
Running up and down hills: some consequences of size
.
Science
178
,
1096
-
1097
.
Theunissen
,
L. M.
,
Vikram
,
S.
and
Dürr
,
V.
(
2014
).
Spatial co-ordination of foot contacts in unrestrained climbing insects
.
J. Exp. Biol.
217
,
3242
-
3253
.
Walker
,
J. A.
(
1998
).
Estimating velocities and accelerations of animal locomotion: a simulation experiment comparing numerical differentiation algorithms
.
J. Exp. Biol.
201
,
981
-
995
.
Wang
,
Z.
,
Dai
,
Z.
,
Ji
,
A.
,
Ren
,
L.
,
Xing
,
Q.
and
Dai
,
L.
(
2015
).
Biomechanics of gecko locomotion: the patterns of reaction forces on inverted, vertical and horizontal substrates
.
Bioinspir. Biomim.
10
,
016019
.
Watson
,
J. T.
,
Ritzmann
,
R. E.
,
Zill
,
S. N.
and
Pollack
,
A. J.
(
2002
).
Control of obstacle climbing in the cockroach, Blaberus discoidalis
.
I. Kinematics. J. Comp. Physiol. A
188
,
39
-
53
.
Wiens
,
J. J.
,
Fetzner
,
J. W.
, Jr
,
Parkinson
,
C. L.
and
Reeder
,
T. W.
(
2005
).
Hylid frog phylogeny and sampling strategies for speciose clades
.
Syst. Biol.
54
,
778
-
807
.
Winter
,
B
. (
2013
).
Linear models and linear mixed effects models in R with linguistic applications
.
arXiv:13085499
.
Young
,
M. W.
,
Dickinson
,
E.
,
Flaim
,
N. D.
and
Granatosky
,
M. C.
(
2022a
).
Overcoming a ‘forbidden phenotype’: the parrot's head supports, propels and powers tripedal locomotion
.
Proc. Biol. Sci.
289
,
20220245
.
Young
,
M. W.
,
Lynch
,
S. K.
,
Dickinson
,
E.
,
Currier
,
A. A.
, III
,
Davoli
,
E. C.
,
Hanna
,
C. S.
,
Fischer
,
H. M.
,
Diubaldi
,
G. A.
and
Granatosky
,
M. C.
(
2022b
).
Patterns of single limb forces during terrestrial and arboreal locomotion in rosy-faced lovebirds (Psittaciformes: Agapornis roseicollis)
.
J. Exp. Biol.
225
,
jeb244571
.
Zampagni
,
M. L.
,
Brigadoi
,
S.
,
Schena
,
F.
,
Tosi
,
P.
and
Ivanenko
,
Y. P.
(
2011
).
Idiosyncratic control of the center of mass in expert climbers
.
Scand. J. Med. Sci. Sports
21
,
688
-
699
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information