Graded substrates require legged animals to modulate their limb mechanics to meet locomotor demands. Previous work has elucidated strategies used by cursorial animals with upright limb posture, but it remains unclear how sprawling species such as alligators transition between grades. We measured individual limb forces and 3D kinematics as alligators walked steadily across level, 15 deg incline and 15 deg decline conditions. We compared our results with the literature to determine how limb posture alters strategies for managing the energetic variation that accompanies shifts in grade. We found that juvenile alligators maintain spatiotemporal characteristics of gait and locomotor speed while selectively modulating craniocaudal impulses (relative to level) when transitioning between grades. Alligators seem to accomplish this using a variety of kinematic strategies, but consistently sprawl both limb pairs outside of the parasagittal plane during decline walking. This latter result suggests alligators and other sprawling species may use movements outside of the parasagittal plane as an axis of variation to modulate limb mechanics when transitioning between graded substrates. We conclude that limb mechanics during graded locomotion are fairly predictable across quadrupedal species, regardless of body plan and limb posture, with hindlimbs playing a more propulsive role and forelimbs functioning to dissipate energy. Future work will elucidate how shifts in muscle properties or function underlie such shifts in limb kinematics.

Legged locomotion on land often requires animals to transition between level, incline and decline substrates. Graded substrates change the mechanical energy requirement of steady-state locomotion and animals must either perform more work (inclined surfaces) or dissipate more work (declined surfaces) with their limbs to compensate. Most studies of legged animals have focused on steady-state locomotion on level ground (Biewener and Daley, 2007; Birn-Jeffery and Higham, 2014), which establishes an important baseline for relating mechanics and energetics, but does not fully represent the diversity of locomotor conditions encountered in nature.

Data on individual limb forces across grades are available in a small subset of animals, most of which are cursorial species adapted for running. These include horses (Dutto et al., 2004), cats (Gregor et al., 2006), dogs (Lee, 2011), goats (Arnold et al., 2013) and primates (Gottschall and Kram, 2005; Schmidt, 2005). While the strategies that cursorial quadrupeds use to navigate grade changes are well understood, our understanding of how quadrupeds with sprawled postures modulate their limb mechanics when transitioning between grades remains limited. Recent studies in lizards have advanced our understanding of how limb mechanics are used to navigate the challenges of arboreal locomotion (Foster and Higham, 2012; Krause and Fischer, 2013; Wang et al., 2015).

Legged animals use some general strategies to transition between grades. Animals tend to reduce forward speed, increase duty factor and flex their limbs on graded conditions when compared with level (Birn-Jeffery and Higham, 2014). On inclines, most animals decrease stride length and increase stance phase duration relative to level (Carlson-Kuhta et al., 1998; Irschick and Jayne, 1998; Foster and Higham, 2012; Williams et al., 2009), usually resulting in increased vertical or craniocaudal (in the direction of travel) ground reaction forces (GRFs) (Dutto et al., 2004; Gottschall and Kram, 2005; Lammers et al., 2006). On declines, locomotor strategies are more variable, but common changes made by legged animals suggest that stability may be more important than energy modulation when moving downhill. Spatiotemporal characteristics of gait do not follow a general pattern on declines (Birn-Jeffery and Higham, 2014), but limb forces shift to net braking (Lammers et al., 2006; Lee, 2011) and limb joints typically flex more during stance phase relative to level (Jayne and Irschick, 1999; Krause and Fischer, 2013; Leroux et al., 2002; Smith et al., 1998). This increased limb flexion not only allows energy dissipation by large extensor muscles but also is thought to bring the center of mass (COM) closer to the ground and counteract the ‘toppling torques’ encountered when moving across non-level substrate. Overall, strategies to transition between grades depend on context and may be shaped by divergent demands for stability and efficiency.

American alligators, Alligator mississippiensis, are quadrupeds that walk primarily using a ‘walking trot’ and sprawled posture termed ‘semi-upright’ (Bakker, 1971), where proximal limb segments are abducted away from the body and hands, and feet do not lie directly underneath the body. While these general categories of posture serve as a convenient shorthand, the degrees of freedom in jointed systems limit their universal utility (Gatesy and Biewener, 1991). Alligators also have relatively shorter limbs than most terrestrial vertebrates (Blob, 2000) and drag around a large tail that accounts for ∼28% of their total body mass (Willey et al., 2004). These unique characteristics influence limb kinetics during level locomotion compared with other terrestrial quadrupeds. Specifically, mediolateral forces (forces directed toward the body midline) in alligator limbs are 100–140% the magnitude of craniocaudal forces (forces directed in the direction of travel) and hindlimbs support more body weight than forelimbs during walking (Willey et al., 2004), a feature shared broadly with primates (Pontzer et al., 2014), lizards (McElroy et al., 2014) and salamanders (Nyakatura et al., 2019). However, mediolateral forces in cursorial species are often negligible (e.g. Corbee et al., 2014) and forelimbs typically support more body weight in most quadrupeds (Arnold et al., 2013). Alligators and other sprawlers also use non-parasagittal limb movements during walking, which offers another axis of variation to modulate during locomotion across varying grades.

Here, we quantified the kinetic and kinematic changes used to transition between level, incline and decline walking in alligators. Specifically, we predicted that, similar to other quadrupeds, alligators will (1) support more body weight with their hindlimbs during incline and their forelimbs during decline walking, (2) use their hindlimbs to produce larger propulsive forces during incline walking compared with level walking, (3) use their forelimbs to produce larger braking forces during decline walking compared with level walking. Additionally, given the observation that proximal muscles generate more mechanical work during legged locomotion, we predicted that alligators will primarily modulate their proximal limb joint kinematics to transition across conditions. Finally, we predicted that because of their sprawled posture, alligators will modulate mediolateral kinetics and kinematics to navigate sloped terrain. We tested these predictions by measuring individual limb GRFs as alligators walk at a steady speed across a trackway and 3D high-speed videography to quantify kinematic changes across conditions.

Animals and experimental design

Three female juvenile American alligators, Alligator mississippiensis Daudin 1801 (body mass 524.4±78.4 g, snout–vent length 29.4±0.5 cm; mean±s.d.) were obtained from California State University, San Bernardino, USA, and transferred to an indoor vivarium at University of California Irvine, USA. The vivarium was kept at 28°C with a 12 h:12 h light:dark cycle with the humidity of the vivarium set at 75%. Large tanks within the vivarium allowed animals to fully submerge into water and emerge onto basking platforms ad libitum. Animals were fed weekly and trained to walk across a custom trackway. All animal husbandry and experimental protocols were approved by the University of California Irvine Institutional Animal Care and Use Committee (IACUC).

The trackway (2.5 m×0.3 m) was instrumented with two high-speed video cameras and a small circular force plate (∼0.018 m2) flush-mounted with the ground (Fig. 1A). The trackway was kept in a temperature-controlled room where the environmental temperature was maintained at 28.0±1.0°C (TC1 Controller, Environmental Growth Chambers, Chagrin Falls, OH, USA) with a humidity of 65–70%. This allowed all experiments to be conducted within alligators' preferred temperature range (Smith, 1975; Brattstrom, 1965). Level, incline and decline data were collected over a period of 4 weeks and animal body mass, total length and snout–vent length were measured on each day of data collection to account for minor changes across the span of this study. The substrate was made of acrylic lined with either 100-grit (level) or 120-grit sandpaper (incline and decline).

Fig. 1.

Schematic diagram of the of experimental set-up and kinematic analyses. (A) Dual high-speed video cameras recorded alligators as they walked steadily across a trackway. A small circular force plate flush-mounted with the ground measured individual limb forces during stance phase for either one forelimb or one hindlimb step per trial. The trackway was modified to a 15 deg incline and 15 deg decline following level data collection. Alligator image adapted from Reilly and Elias (1998). (B) 3D joint angles were calculated based on kinematic markers. Protraction/retraction angles for both limbs were based on markers along the vertebral column, right proximal joint, and the elbow (for shoulder protraction/retraction) or knee (for hip protraction/retraction). Abduction/adduction angles were based on markers on the left proximal joint, right proximal joint, and the elbow (for shoulder abduction/adduction) or knee (for hip abduction/adduction). (C) Sprawl area is defined as the area (in cm2) under the triangle formed by proximal limb joint height (hshoulder shown), effective limb length (ELL) and step width (wstep) and was used as a proxy for degree of sprawled posture.

Fig. 1.

Schematic diagram of the of experimental set-up and kinematic analyses. (A) Dual high-speed video cameras recorded alligators as they walked steadily across a trackway. A small circular force plate flush-mounted with the ground measured individual limb forces during stance phase for either one forelimb or one hindlimb step per trial. The trackway was modified to a 15 deg incline and 15 deg decline following level data collection. Alligator image adapted from Reilly and Elias (1998). (B) 3D joint angles were calculated based on kinematic markers. Protraction/retraction angles for both limbs were based on markers along the vertebral column, right proximal joint, and the elbow (for shoulder protraction/retraction) or knee (for hip protraction/retraction). Abduction/adduction angles were based on markers on the left proximal joint, right proximal joint, and the elbow (for shoulder abduction/adduction) or knee (for hip abduction/adduction). (C) Sprawl area is defined as the area (in cm2) under the triangle formed by proximal limb joint height (hshoulder shown), effective limb length (ELL) and step width (wstep) and was used as a proxy for degree of sprawled posture.

Kinematics

Joint positions were determined based on kinematic markers placed near each major limb joint's center of rotation. Kinematic markers (OptiTrack Facial Markers 3 mm, NaturalPoint Inc., Corvallis, OR, USA) were painted white and black to facilitate digitization and were placed on the right side of the animal. Forelimb (FL) markers were placed at the third metacarpophalangeal joint, the styloid process of the radius, the lateral epicondyle of the humerus, and directly above the right and left glenoid cavities. Hindlimb (HL) markers were placed at the third metatarsophalangeal joint, the lateral malleolus of the fibula, the lateral epicondyle of the femur, and directly above the right and left acetabula (see Reilly and Elias, 1998). Additional markers were placed along the vertebral column at 30% and 70% gleno-acetabular length. We assumed the marker at 70% along the craniocaudal axis represents COM location in juvenile alligators (Willey et al., 2004) and its displacement over time was used to calculate instantaneous forward walking velocity.

High-speed video cameras with a maximum resolution of 1280×1024 pixels (SC1, Edgertronic, Sanstreak Corp., San Jose, CA, USA) recorded dorsal and lateral views at 200 Hz as animals walked the length of the trackway. Each trial included several steps before and after the force plate, in addition to the step of interest on the force plate. Videographic data were used to determine stance phase duration (SP), swing phase duration (SW) of the pre-force plate step, and SW of the post-force plate step. The average between the two SW values was used to calculate stride duration for each trial, which was then used to calculate stride frequency (=1/stride duration), stride length (=stride frequency×speed) and duty factor (=stance phase duration/stride duration). Coordinates of kinematic markers overlying joint positions were determined by calibrating a 96-point calibration cube with direct-linear transformation software (Hedrick, 2008) in MATLAB (The MathWorks, Natick, MA, USA). 3D position data were used to calculate joint angles of each major limb joint (Fig. 1B), proximal limb joint height (in cm), effective limb length (in cm; the average 3D distance between the proximal limb joint and limb point of contact), step width (in cm; the average mediolateral distance between the proximal limb joint and limb point of contact), and sprawl area [in cm2; the area within the triangle formed by proximal joint height, effective limb length (ELL) and step width] during stance phase (Fig. 1C). Given the limitation of our markers, our analysis is best described as quasi-3D and quasi-static as we are unable to distinguish between the rotational axes used to arrive at a given joint position (Gatesy et al., 2010).

GRFs

A custom force plate (Gamma F/T Sensor ∼75 mm in diameter, ATI Industrial Automation Inc., Apex, NC, USA) fitted with an acrylic mask of the same diameter was used to record GRFs in vertical, craniocaudal and mediolateral directions of either a FL or HL step during a given walking trial. The dimensions of our sensor allowed us to accurately measure single foot forces (Movies 1–3); however, approximately 70% of walking trials were discarded because of partial contact with the force plate. The force plate was calibrated immediately before each experiment at the temperature and humidity used during data collection. Positive values in the vertical direction represent forces against gravity, positive values in the craniocaudal direction represent propulsive forces whereas negative values are braking forces, and positive values in the mediolateral direction represent lateral pushes from the limb. GRF data were collected at 1000 Hz using an A/D converter (National Instruments, Austin, TX, USA) and analyzed using Igor Pro software (Wavemetrics Inc., Portland, OR, USA). GRF and kinematics data were synchronized with a common external trigger. Raw forces were normalized to body weight units (BWUs) for each day of data collection. Normalized force traces were smoothed using a cubic spline function before calculation of impulses for each GRF component. Average GRF angle was calculated for both parasagittal and frontal planes using the following equations (Lee et al., 2004):
(1)
(2)
where iCC is net craniocaudal impulse, iV is net vertical impulse and iML is net mediolateral impulse. Average GRF angles based on impulse magnitudes best describe a limb's braking and propulsive efforts (θparasagittal), as well as laterally directed forces (θfrontal).

Level data were collected over the course of 3–4 days and once a sufficient number of trials had been collected from each animal, the trackway was modified to a 15 deg incline, and then a 15 deg decline (Fig. 1A). Animals usually walked on a level, incline or decline trackway for any given day of data collection. We did experiments with steeper slopes, but found less reliable GRF data during decline locomotion as steps became more variable and were often impacted by minor slips. Animals were encouraged to walk at their preferred speed via tail taps and allowed to rest after 4–5 consecutive walking bouts. Only trials where animals walked at a constant forward velocity (i.e. steady state), contacted the force plate with either the right FL or right HL in isolation, and remained in view of both cameras were included in our analyses. We considered animals to be walking steadily based on average speed calculated during two different time intervals and the quotient between them. Speed during stance phase (speedSP) was determined by taking the average instantaneous forward velocity of the COM during limb contact with the force plate. Speed during the total trial (speedtot) was determined by taking the average instantaneous forward velocity as animals walked the entire length of the trackway. In this study, trials were considered steady state if the speedSP/speedtot quotient fell between 0.9 and 1.1. A total of 96 steps from all individuals were analyzed during stance phase of walking across level, incline and decline conditions (n=48 forelimb steps, n=48 hindlimb steps). Alligator steps can be sloppy and often involve toe and finger dragging; therefore, we define stance phase in this study as the time between hand or foot flat on the substrate until the onset of significant wrist palmarflexion (for FL) or when the foot segment was perpendicular to the floor (for HL).

Statistical analyses

A linear mixed effects (LME) model was used to test for effects of condition (level, incline or decline) on the following dependent variables for each limb: stance phase duration, stride length, duty factor, GRF impulses, average GRF angles (parasagittal and frontal), proximal limb joint height, ELL, sprawl area and net excursion (=start angle−end angle) of each major limb joint during stance phase. Condition was treated as a fixed effect whereas individual alligator was treated as a random effect. We used a sequential Bonferroni correction to account for multiple statistical tests for each limb (LME α=0.006). All variables found to be significantly affected by condition in our LME model were further analyzed using post hoc pairwise comparison Tukey's tests (α=0.05). All statistics were performed using R (lme4 and emmeans packages; https://CRAN.R-project.org/package=lme4, https://CRAN.R-project.org/package=emmeans).

Speed and stride parameters

Condition (i.e. level, incline, decline) had no significant effect on speedSP, stance phase duration, stride length or duty factor for either limb (Table 1). Pooled (FL and HL trials) speedSP was 0.118±0.026 m s−1 during level walking, 0.116±0.041 m s−1 during incline walking and 0.116±0.028 m s−1 during decline walking, indicating no difference in average locomotor speed as alligators transition across grades.

Table 1.

Summary of linear mixed effects model results for spatiotemporal characteristics across grades

Summary of linear mixed effects model results for spatiotemporal characteristics across grades
Summary of linear mixed effects model results for spatiotemporal characteristics across grades

GRFs

We observed no significant differences in FL vertical impulse (LME, P=0.968), HL vertical impulse (P=0.110), FL mediolateral impulse (P=0.711) or HL mediolateral impulse (P=0.112) due to condition. However, craniocaudal impulse decreased significantly during decline walking and increased during incline walking (relative to level) for both limbs (FL LME, P<0.0001; HL LME, P<0.0001) (Fig. 2, Table 2). Consistent with these data, we found that average parasagittal GRF angle (θparasagittal) in both limbs became more negative during decline walking (indicating an increased braking effort) and more positive during incline walking (indicating a more propulsive effort) for each limb (FL LME, P<0.0001; HL LME, P<0.0001; Fig. 3A). We also found no significant differences in average frontal GRF angle (θfrontal) due to condition for either limb (FL LME, P=0.416; HL LME, P=0.919; Fig. 3B).

Fig. 2.

Representative ground reaction force (GRF) profiles under the three conditions. GRF is shown (in body weight units, BWUs) for the forelimb (left) and hindlimb (right) from (A) incline, (B) level and (C) decline trials normalized to percentage stance phase. Vertical (black) and mediolateral (gray) impulses did not change noticeably with condition; however, craniocaudal impulse (red) increased during incline and decreased during decline (relative to level) for each limb. Each profile is from the same alligator walking a steady-state speed.

Fig. 2.

Representative ground reaction force (GRF) profiles under the three conditions. GRF is shown (in body weight units, BWUs) for the forelimb (left) and hindlimb (right) from (A) incline, (B) level and (C) decline trials normalized to percentage stance phase. Vertical (black) and mediolateral (gray) impulses did not change noticeably with condition; however, craniocaudal impulse (red) increased during incline and decreased during decline (relative to level) for each limb. Each profile is from the same alligator walking a steady-state speed.

Fig. 3.

Box plots of average parasagittal GRF angle and average frontal GRF angle during stance phase. Each limb shifted to a more negative (or braking) parasagittal GRF angle during decline, and a more positive (or propulsive) angle during incline [linear mixed effect (LME) model, P<0.0001]. Average frontal GRF angle in the hindlimb decreased slightly during incline, but was not significantly different from the other two conditions (LME, P=0.92). Box plots show median, upper and lower quartiles and 1.5× interquartile range; circles are individual data points. Different lowercase letters denote significant differences between groups.

Fig. 3.

Box plots of average parasagittal GRF angle and average frontal GRF angle during stance phase. Each limb shifted to a more negative (or braking) parasagittal GRF angle during decline, and a more positive (or propulsive) angle during incline [linear mixed effect (LME) model, P<0.0001]. Average frontal GRF angle in the hindlimb decreased slightly during incline, but was not significantly different from the other two conditions (LME, P=0.92). Box plots show median, upper and lower quartiles and 1.5× interquartile range; circles are individual data points. Different lowercase letters denote significant differences between groups.

Table 2.

Summary of LME model results for force-derived variables during stance phase across grades

Summary of LME model results for force-derived variables during stance phase across grades
Summary of LME model results for force-derived variables during stance phase across grades

Whole-limb kinematics

Condition had no significant effect on shoulder height (LME, P=0.22) or hip height (P=0.84) (Fig. 4). ELL increased significantly during decline locomotion relative to the other two conditions in both limbs (FL P=0.002; HL P=0.005). Consistent with these data, we found that sprawl area significantly increased during decline relative to the other two conditions in both limbs (FL P=0.002; HL P=0.004), indicating alligator limbs are more abducted during decline walking.

Fig. 4.

Box plots of average proximal limb joint height, ELL and sprawl area during stance phase. Shoulder and hip height did not vary across locomotor condition (LME, P=0.22 and P=0.84, respectively). Increased ELL (forelimb P=0.002, hindlimb P=0.005) and sprawl area (forelimb P=0.005, hindlimb P=0.004) indicate both limbs were more sprawled out (i.e. laterally displaced) during decline locomotion. Box plots as in Fig. 3. Different lowercase letters denote significant differences between groups.

Fig. 4.

Box plots of average proximal limb joint height, ELL and sprawl area during stance phase. Shoulder and hip height did not vary across locomotor condition (LME, P=0.22 and P=0.84, respectively). Increased ELL (forelimb P=0.002, hindlimb P=0.005) and sprawl area (forelimb P=0.005, hindlimb P=0.004) indicate both limbs were more sprawled out (i.e. laterally displaced) during decline locomotion. Box plots as in Fig. 3. Different lowercase letters denote significant differences between groups.

Limb joint kinematics

We observed significant inter- and intra-individual variation in joint kinematics during stance phase of walking in alligators (Figs S1 and S2). With the exception of wrist net excursion during decline (Fig. 5, Table 3), most joint kinematic variables tested (e.g. start/end angles, joint flexion/extension) yielded non-significant results when run through our LME model. We found that the wrist has a net dorsiflexion excursion during level and incline walking (−8.8±7.5 deg and −10.8±7.7 deg, respectively), but shifts to a net palmarflexion excursion (1.4±14.0 deg) during decline walking (LME, P=0.005).

Fig. 5.

Average joint angle trajectories during each condition. Forelimb joint kinematics were not significantly affected by condition, apart from increased wrist plantarflexion during decline (LME P=0.005). Most changes in hindlimb kinematics due to condition were not significantly different from level (LME, P>0.006). Each line represents the pooled mean across all alligators for a given condition. The 95% confidence intervals for the level condition are denoted by gray shading; confidence intervals for the incline and decline conditions were omitted for clarity (see Figs S1 and S2 for the complete dataset).

Fig. 5.

Average joint angle trajectories during each condition. Forelimb joint kinematics were not significantly affected by condition, apart from increased wrist plantarflexion during decline (LME P=0.005). Most changes in hindlimb kinematics due to condition were not significantly different from level (LME, P>0.006). Each line represents the pooled mean across all alligators for a given condition. The 95% confidence intervals for the level condition are denoted by gray shading; confidence intervals for the incline and decline conditions were omitted for clarity (see Figs S1 and S2 for the complete dataset).

Table 3.

Summary of LME model results for net joint excursion during stance phase across grades

Summary of LME model results for net joint excursion during stance phase across grades
Summary of LME model results for net joint excursion during stance phase across grades

Body weight support across grades

We predicted that alligators would support a larger proportion of their body weight with their hindlimbs during incline walking and their forelimbs during decline walking. This prediction was partly based on the assumption that juvenile alligators had a COM closer to their hindlimbs (Willey et al., 2004). However, we found that vertical impulses were statistically similar between forelimbs and hindlimbs during level walking in the alligators used in this study (Table 2). Recent work shows that COM position shifts its craniocaudal position across the size range in alligators (Iijima et al., 2021); therefore, it is possible that the alligators used in this study (body mass range 0.47–0.61 kg) have a COM closer to 50% of their gleno-acetabular length.

Regardless of COM distribution, other quadrupeds such as horses (Dutto et al., 2004), dogs (Lee, 2011) and goats (Arnold et al., 2013) shift body weight support between limb pairs during graded locomotion when compared with level locomotion. This pattern is also observed in primates when comparing the relative role of limb pairs during walking versus climbing (Hanna et al., 2017). However, we did not observe significant differences in vertical impulses for either limb in alligators due to condition (Table 2). This suggests a 15 deg locomotor grade does not shift the relative contributions of the forelimb and hindlimb in juvenile alligators. It remains unclear whether differences in body weight support across grades between alligators and other quadrupeds are due to differences in body plan (i.e. relative body mass, COM distribution), limb posture or a combination of the two. Future work manipulating COM position with the addition of extrinsic loads in quadrupeds of varying sizes is needed to properly address the effect of COM position on limb mechanics during locomotion. In general, alligators seem to maintain consistent vertical impulses across conditions while selectively modulating craniocaudal impulses to maintain forward speed as they transition between grades (see below). This allows alligators to manage the increase in total GRF magnitude necessary for graded locomotion while prioritizing speed maintenance when transitioning between grades.

Braking and propulsive forces across grades

We predicted alligators would use their hindlimbs to produce larger propulsive forces during incline walking and their forelimbs to produce larger braking forces during decline walking (when compared with level walking). As predicted for the hindlimb, we found that θparasagittal increased from 7 deg for level walking to 22 deg for incline walking, indicating a stronger propulsive effort by the hindlimb when walking uphill. As predicted for the forelimb, we found that θparasagittal decreased from −5 deg for level walking to −22 deg for decline walking, indicating a stronger braking effort by the forelimb when walking downhill (Fig. 3). These results are consistent with individual limb data from dogs (Lee, 2011), horses (Dutto et al., 2004) and opossums (Lammers et al., 2006) moving across varying grades. In general, it seems that the hindlimbs of quadrupeds are better suited for propulsive force generation, whereas the forelimbs are better suited for braking. This is consistent with anatomical data in mammalian quadrupeds showing significantly more extensor muscle mass in the hindlimbs when compared with the forelimbs (Kilbourne and Hoffman, 2013), a pattern also observed in alligators (Allen et al., 2010, 2014). Alligators may modulate limb muscle recruitment or function (length trajectories) to maintain similar limb forces despite the dimensional disparity among forelimbs and hindlimbs. Future work measuring in vivo muscle function in forelimb and hindlimb muscles during quadrupedal locomotion will likely reveal how mechanical function shapes the physiology and control of limb muscles.

Limb joint kinematics across grades

We predicted alligators would primarily modulate proximal limb joint kinematics to transition across conditions. Most legged animals increase hip retraction (i.e. extension) on inclines (Carlson-Kuhta et al., 1998; Hoyt et al., 2002; Jayne and Irschick, 1999; Reznick et al., 2021) and protract (i.e. flex) limbs more on decline relative to level walking (Jayne and Irschick, 1999; Reznick et al., 2021; Smith et al., 1998). We found some support for alligators increasing hip retraction during incline when compared with level walking (Fig. 5), but this difference was not statistically significant (Table 3). We did however observe that the range of joint angles used tended to increase during incline walking compared with level locomotion (Fig. 5, Table 3). This pattern is consistent with the need to increase positive mechanical work production when walking up an incline and is most noticeable for the joints of the hindlimb. A similar pattern has also been observed in geckos (Birn-Jeffery and Higham, 2016) and goats (Arnold et al., 2013).

Our observation that alligators used a more sprawled limb posture during decline walking suggests that alligators may use movements outside of the parasagittal plane to improve static stability by increasing the base of support when walking downhill. Most legged vertebrates studied to date have parasagittal limb motions during locomotion and increase stability on grades by decreasing proximal joint height and bringing their COM closer to the ground (Birn-Jeffery and Higham, 2014). Decreased proximal joint height in upright animals causes increased limb joint flexion and decreased effective mechanical advantage across limb joints (Arnold et al., 2005; Biewener, 1989). Alligators and other sprawlers likely benefit from non-parasagittal limb movements as a mechanism for increasing stability without flexing their limbs. This interpretation is nicely exemplified by the limb posture in polypedal organisms (Martinez, 1996) or robots (Krummel et al., 2014), where non-parasagittal changes in posture are often associated with increased stability. It is also possible that the large tail of alligators, which drags on the substrate during locomotion, contributes to passive stability, which may alleviate the need for kinematic modulation of limbs across grades. Future work on the non-parasagittal limb movements used by sprawling species transitioning between grades is needed to support this notion.

Conclusion

In this study, we examined how alligators transition between level, incline and decline walking. We measured individual limb forces and 3D kinematics as alligators walked steadily across level, 15 deg incline and 15 deg decline conditions. We compared our results with strategies used by other legged animals in an effort to identify which features of an animal's body plan and locomotor mechanics effect overall limb function. We found that juvenile alligators maintain spatiotemporal characteristics of gait and locomotor speed by selectively modulating craniocaudal impulses (relative to level locomotion) when transitioning between grades. Alligators seem to accomplish this using a variety of kinematic strategies, but consistently sprawl both limb pairs outside of the parasagittal plane during decline walking. This latter result suggests alligators and other sprawling species likely benefit from movements outside of the parasagittal plane as another axis of variation to modulate when transitioning between graded substrates.

The authors would like to thank T. Owerkowicz for providing the alligators for the study. We thank members of the Azizi Lab for helpful comments on an early draft of the manuscript. Thank you to J. R. Hutchinson and one other anonymous reviewer for their comments that significantly improved the paper.

Author contributions

Conceptualization: A.A.A., E.A.; Methodology: A.A.A., E.A.; Validation: A.A.A.; Formal analysis: A.A.A., E.A.; Investigation: A.A.A., E.A.; Resources: E.A.; Data curation: A.A.A., E.A.; Writing - original draft: A.A.A.; Writing - review & editing: A.A.A., E.A.; Visualization: A.A.A., E.A.; Supervision: E.A.; Project administration: E.A.

Funding

This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

Data availability

All relevant data can be found within the article and its supplementary information.

Allen
,
V.
,
Elsey
,
R. M.
,
Jones
,
N.
,
Wright
,
J.
and
Hutchinson
,
J. R.
(
2010
).
Functional specialization and ontogenetic scaling of limb anatomy in Alligator mississippiensis
.
J. Anat.
216
,
423
-
445
.
Allen
,
V.
,
Molnar
,
J.
,
Parker
,
W.
,
Pollard
,
A.
,
Nolan
,
G.
and
Hutchinson
,
J. R.
(
2014
).
Comparative architectural properties of limb muscles in Crocodylidae and Alligatoridae and their relevance to divergent use of asymmetrical gaits in extant Crocodylia
.
J. Anat.
225
,
569
-
582
.
Arnold
,
A. S.
,
Anderson
,
F. C.
,
Pandy
,
M. G.
and
Delp
,
S. L.
(
2005
).
Muscular contributions to hip and knee extension during the single limb stance phase of normal gait: a framework for investigating the causes of crouch gait
.
J. Biomech.
38
,
2181
-
2189
.
Arnold
,
A. S.
,
Lee
,
D. V.
and
Biewener
,
A. A.
(
2013
).
Modulation of joint moments and work in the goat hindlimb with locomotor speed and surface grade
.
J. Exp. Bio.
216
,
2201
-
2212
.
Bakker
,
R. T.
(
1971
).
Dinosaur physiology and the origin of mammals
.
Evolution
25
,
636
-
658
.
Biewener
,
A. A.
(
1989
).
Scaling body support in mammals: limb posture and muscle mechanics
.
Science
245
,
45
-
48
.
Biewener
,
A. A.
and
Daley
,
M. A.
(
2007
).
Unsteady locomotion: integrating muscle function with whole body dynamics and neuromuscular control
.
J. Exp. Biol.
210
,
2949
-
2960
.
Birn-Jeffery
,
A. V.
and
Higham
,
T. E.
(
2014
).
The scaling of uphill and downhill locomotion in legged animals
.
Int. Org. Biol.
54
,
1159
-
1172
.
Birn-Jeffery
,
A. V.
and
Higham
,
T. E.
(
2016
).
Geckos decouple fore-and hind limb kinematics in response to changes in incline
.
Front. Zool.
13
,
11
.
Blob
,
R. W.
(
2000
).
Interspecific scaling of the hindlimb skeleton in lizards, crocodilians, felids and canids: does limb bone shape correlate with limb posture?
.
J. Zool.
250
,
507
-
531
. .
Brattstrom
,
B. H.
(
1965
).
Body temperatures of reptiles
.
Am. Midl. Nat.
73
,
376
-
422
.
Carlson-Kuhta
,
P.
,
Trank
,
T. V.
and
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
.
Corbee
,
R. J.
,
Maas
,
H.
,
Doornenbal
,
A.
and
Hazewinkel
,
H. A. W.
(
2014
).
Forelimb and hindlimb ground reaction forces of walking cats: assessment and comparison with walking dogs
.
Vet. J.
202
,
116
-
127
.
Dutto
,
D. J.
,
Hoyt
,
D. F.
,
Cogger
,
E. A.
and
Wickler
,
S. J.
(
2004
).
Ground reaction forces in horses trotting up an incline and on the level over a range of speeds
.
J. Exp. Biol.
207
,
3507
-
3514
.
Foster
,
K. L.
and
Higham
,
T. E.
(
2012
).
How forelimb and hindlimb function changes with incline and perch diameter in the green anole, Anolis carolinensis
.
J. Exp. Biol.
215
,
2288
-
2300
.
Gatesy
,
S. M.
and
Biewener
,
A. A.
(
1991
).
Bipedal locomotion: effects of speed, size and limb posture in birds and humans
.
J. Zool.
224
,
127
-
147
.
Gatesy
,
S. M.
,
Baier
,
D. B.
,
Jenkins
,
F. A.
and
Dial
,
K. P.
(
2010
).
Scientific rotoscoping: a morphology-based method of 3-D motion analysis and visualization
.
J. Exp. Zool.
313A
,
244
-
261
.
Gottschall
,
J. S.
and
Kram
,
R.
(
2005
).
Ground reaction forces during downhill and uphill running
.
J. Biomech.
38
,
445
-
452
.
Gregor
,
R. J.
,
Smith
,
D. W.
and
Prilutsky
,
B. I.
(
2006
).
Mechanics of slope walking in the cat: quantification of muscle load, length change, and ankle extensor EMG patterns
.
J. Neurophysiol.
3
,
1397
-
1409
.
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
.
Hedrick
,
T. L.
(
2008
).
Software techniques for two-and three-dimensional kinematic measurements of biological and biomimetic systems
.
Bioinspir. Biomimet.
3
,
034001
.
Hoyt
,
D. F.
,
Molinari
,
M.
,
Wickler
,
S. J.
and
Cogger
,
E. A.
(
2002
).
Effect of trotting speed, load and incline on hindlimb stance–phase kinematics
.
Equine Vet. J.
34
,
330
-
336
.
Iijima
,
M.
,
Munteanu
,
V. D.
,
Elsey
,
R. M.
and
Blob
,
R. W.
(
2021
).
Ontogenetic changes in limb posture, kinematics, forces and joint moments in American alligators (Alligator mississippiensis)
.
J. Exp. Biol.
224
,
242990
.
Irschick
,
D. J.
and
Jayne
,
B. C.
(
1998
).
Effects of incline on speed, acceleration, body posture and hindlimb kinematics in two species of lizard Callisaurus draconoides and Uma scoparia
.
J. Exp. Biol.
201
,
273
-
287
.
Jayne
,
B. C.
and
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
.
Kilbourne
,
B. M.
and
Hoffman
,
L. C.
(
2013
).
Scale effects between body size and limb design in quadrupedal mammals
.
PLoS ONE
8
,
e78392
.
Krause
,
C.
and
Fischer
,
M. S.
(
2013
).
Biodynamics of climbing: effects of substrate orientation on the locomotion of a highly arboreal lizard (Chamaeleo calyptratus)
.
J. Exp. Biol.
216
,
1448
-
1457
.
Krummel
,
G.
,
Kaipa
,
K. N.
and
Gupta
,
S. K.
(
2014
).
A Horseshoe Crab Inspired Surf Zone Robot With Righting Capabilities
.
Proceedings of the ASME 2014 Volume 5A: 38th Mechanisms and Robotics Conference
.
Buffalo, New York
,
USA
.
August 17–20, 2014. V05AT08A010. ASME
.
Lammers
,
A. R.
,
Earls
,
K. D.
and
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
.
Lee
,
D. V.
(
2011
).
Effects of grade and mass distribution on the mechanics of trotting in dogs
.
J. Exp. Biol.
3
,
402
-
411
.
Lee
,
D. V.
,
Stakebake
,
E. F.
,
Walter
,
R. M.
and
Carrier
,
D. R.
(
2004
).
Effects of mass distribution on the mechanics of level trotting in dogs
.
J. Exp. Biol.
207
,
1715
-
1728
.
Leroux
,
A.
,
Fung
,
J.
and
Barbeau
,
H.
(
2002
).
Postural adaptation to walking on inclined surfaces: I. Normal strategies
.
Gait Posture
15
,
64
-
74
.
Martinez
,
M. M.
(
1996
).
Issues for aquatic pedestrian locomotion
.
Am. Zool.
36
,
619
-
627
.
McElroy
,
E. J.
,
Wilson
,
R.
,
Biknevicius
,
A. R.
and
Reilly
,
S. M.
(
2014
).
A comparative study of single-leg ground reaction forces in running lizards
.
J. Exp. Biol.
217
,
735
-
742
.
Epub 2013 Nov 6. PMID: 24198262
.
Nyakatura
,
J. A.
,
Melo
,
K.
,
Horvat
,
T.
,
Karakasiliotis
,
K.
,
Allen
,
V. R.
,
Andikfar
,
A.
,
Andrada
,
E.
,
Arnold
,
P.
,
Lauströer
,
J.
,
Hutchinson
,
J. R.
et al. 
(
2019
).
Reverse-engineering the locomotion of a stem amniote
.
Nature
565
,
351
-
355
.
Pontzer
,
H.
,
Raichlen
,
D. A.
and
Rodman
,
P. S.
(
2014
).
Bipedal and quadrupedal locomotion in chimpanzees
.
J. Hum. Evol.
66
,
64
-
82
.
Reilly
,
S. M.
and
Elias
,
J. A.
(
1998
).
Locomotion in Alligator mississippiensis: kinematic effects of speed and posture and their relevance to the sprawling-to-erect paradigm
.
J. Exp. Biol.
201
,
2559
-
2574
.
Reznick
,
E.
,
Embry
,
K. R.
,
Neuman
,
R.
,
Bolívar-Nieto
,
E.
,
Fey
,
N. P.
and
Gregg
,
R. D.
(
2021
).
Lower-limb kinematics and kinetics during continuously varying human locomotion
.
Sci. Data
8
,
282
.
Schmidt
,
M.
(
2005
).
Quadrupedal locomotion in squirrel monkeys (Cebidae: Saimiri sciureus): a cineradiographic study of limb kinematics and related substrate reaction forces
.
Am. J. Phys. Anthropol.
128
,
359
-
370
.
Smith
,
E. N.
(
1975
).
Thermoregulation of the American alligator, Alligator mississippiensis
.
Physiol. Zool.
48
,
177
-
194
.
Smith
,
J. L.
,
Carlson-Kuhta
,
P.
and
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
.
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
.
Willey
,
J. S.
,
Biknevicius
,
A. R.
,
Reilly
,
S. M.
and
Earls
,
K. D.
(
2004
).
The tale of the tail: limb function and locomotor mechanics in Alligator mississippiensis
.
J. Exp. Biol.
207
,
553
-
563
.
Williams
,
R. J.
,
Nankervis
,
K. J.
,
Colborne
,
G. R.
,
Marlin
,
D. J.
and
Schroter
,
R. C.
(
2009
).
Heart rate, net transport cost and stride characteristics of horses exercising at walk and trot on positive and negative gradients
.
Comp. Exerc. Physiol.
6
,
113
-
119
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information