DOI: 10.5061/dryad.rn8pk0p82
title: Raw data accompanying: Ground reaction forces in monitor lizards (Varanidae) and the scaling of locomotion in sprawling tetrapods

description: Geometric scaling predicts a major challenge for legged, terrestrial
 locomotion. Locomotor support requirements scale identically with body
 mass (α M1), while force generation capacity should scale α M2/3 as it
 depends on muscle cross-sectional area. Mammals compensate with more
 upright limb postures at larger sizes, but it remains unknown how
 sprawling tetrapods deal with this challenge. Varanid lizards are an ideal
 group to address this question because they cover an enormous body size
 range while maintaining a similar bent-limb posture and body proportions.
 This study reports the scaling of ground reaction forces and duty factor
 for varanid lizards ranging from 7 g 37 kg. Impulses (force x time) scaled
 roughly as predicted by the inverted pendulum model (α M0.99-1.34) while
 peak forces (α M0.73-1.00) scaled higher than expected. Duty factor scaled
 α M0.04 and was higher for the hindlimb than the forelimb. The proportion
 of vertical impulse to total impulse increased with body size, and
 impulses decreased while peak forces increased with speed. These results
 provide valuable data into how locomotor forces vary with body size and
 suggest how other, extinct, sprawling tetrapods may have dealt with the
 biomechanical challenges associated with generating sufficient locomotor
 forces at larger body sizes.
description: Ground reaction forces
 (GRFs) were collected for 30 varanid lizards from 12 species (table 1)
 nearly covering the entire possible range of varanid body sizes. Data were
 recorded from 266 individual stance phases across speeds ranging from 0.96
 m s<sup>-1</sup> to 4.05 m s<sup>-1</sup>,
 comprising 105 forelimb (figure 1<i>b</i>) and 112 hindlimb
 (figure 1<i>c</i>) trials. Forces were collected at 1,000 or
 10,000 Hz using either a custom-built force-plate[12], a force-plate based
 on a Nano17 load cell (ATI Industrial Automation, Apex, NC, USA), or an
 IP65 gamma force-plate (ATI Industrial Automation, Apex, NC, USA). All
 animals (except <i>V. komodoensis, </i>which were zoo animals)
 were wild caught. All force-plates were calibrated using equal weights.
 Trials with aberrant gaits, or where the entire foot did not contact the
 force plate were excluded. Methods were approved under ethics
 SBS/195/12/ARC (QLD), ANA16104 (QLD), ANE1934 (QLD), ANE2054 (QLD), A2450
 (QLD), ANE2054 (QLD), and RA/3/100/1188 (WA), and lizards were collected
 under permits WISP11435612 (QLD) and SF009075 (WA), 08-001092-5 (WA), and
 WA0001919 (QLD). Dorsal and lateral video
 recordings were captured using two high-speed cameras (Fastec IL3-100;
 1280 x 1024 pixels, Fastec Imaging Corporation, San Diego, CA, USA)
 recording at 250 frames per second. Custom MATLAB (Mathworks, Natick, MA,
 USA) programs and an internal trigger were used to synchronize forces and
 videos and identify strides. GRFs were analyzed during the stance phase,
 from footfall to toe-off, and data were smoothed using a fourth-order,
 zero-lag Butterworth low-pass filter. The polarity of forces was
 transformed so that positive values indicate ground support
 (GRF<sub>Z</sub>), cranially-directed forces
 (GRF<sub>X</sub>), and medially-directed forces
 (GRF<sub>Y</sub>) relative to the body (figure
 1<i>e</i>). Negative (GRF<sub>X</sub>) forces
 indicate caudally-directed forces, and negative
 (GRF<sub>Y</sub>) values indicate laterally-directed forces.
 Total force magnitude [1,25] was calculated as in Eq. 1. Impulses were calculated as time integrals of the relevant force component. For cranial, caudal, medial, and lateral impulses, forces were first separated into positive and negative regions before integration of the absolute value.  Average speed was estimated by quantifying the 2D motion of a single marker on the dorsal aspect of the pectoral or pelvic girdle using DeepLabCut, a machine-learning video tracking software [26]. Speeds were converted to relative speed based on the maximum sprint speed for each species (table 1). Relative speed [27] was used instead of dynamic speed (Froude number), because although dynamic speed is useful for identifying similar gaits among diverse species [23] and was well-correlated with relative speed in our dataset (figure S4), it has been shown to be of limited utility in predicting similar gaits in closely-related lizards [28] and it is not obvious which characteristic length would best capture the dynamics of sprawling locomotion. All analyses were conducted using R 3.6.3 (R Development Core Team, 2020). The relationships between force parameters with body mass, relative speed, and limb (forelimb versus hindlimb) were assessed using mixed effects models using the <i>lmer.R </i> function of the lme4 and lmerTest packages including subject as a random factor [29,30]. All continuous variables were log-transformed. In all cases, a full-interaction model was initially constructed. If the interaction terms were not significant, a reduced model with no interaction terms was used instead, provided that an ANOVA test found no significant difference between these models. When full models were used, non-significant variables were successively dropped from the model using the <i>update.R </i>function as long as each, successive model was not found to be significantly different from the previous one using an ANOVA test (table S1). Individual regressions for body mass and relative speed (table 2) were calculated using <i>sim_slopes.R</i>, and plotted using <i>interact_plot</i>.<i>R </i>[31]<i>.</i> Percentage differences of peak forces and impulses between the fore and hindlimbs were calculated as the ratio of the inverse logs from the regression intercept values. To check for an influence of phylogenetic patterns, the slopes and confidence intervals of linear regressions for species means of each stride parameter were compared with phylogenetic independent contrasts of the same parameters (figure S1), calculated using the <i>pic.R </i>function [32]<i> </i>and the maximum-likelihood varanid tree built from 1030 bp of the NADH-2 gene [33]<i>. </i>In order to investigate the scaling of the position of the horizontal centre of mass, previously-published data from the two-scale method [34] (figure S5) were analysed (figure 1<i>a</i>). References 1.        Clemente CJ, Bishop PJ, Newman N, Hocknull SA. 2018 Steady bipedal locomotion with a forward situated whole-body centre of mass: the potential importance of temporally asymmetric ground reaction forces. <i>J. Zool.</i> <b>304</b>, 193–201. (doi:10.1111/jzo.12521) 2.        Blob RW, Biewener AA. 2001 Mechanics of limb bone loading during terrestrial locomotion in the green iguana (<i>Iguana iguana</i>) and American alligator (<i>Alligator mississippiensis</i>). <i>J. Exp. Biol.</i> <b>204</b>, 1099–1122. 3.        Biewener AA. 1989 Scaling body support in mammals: limb posture and muscle mechanics. <i>Science</i> <b>245</b>, 45–48. (doi:10.1126/science.2740914) 4.        Nath T, Mathis A, Chen AC, Patel A, Bethge M, Mathis MW. 2019 Using DeepLabCut for 3D markerless pose estimation across species and behaviors. <i>Nat. Protoc.</i> <b>14</b>, 2152–2176. (doi:10.1038/s41596-019-0176-0) 5.        Van Damme R, Aerts P, Vanhooydonck B. 1998 Variation in morphology, gait characteristics and speed of locomotion in two populations of lizards. <i>Biol. J. Linn. Soc.</i> <b>63</b>, 409–427. (doi:10.1006/bijl.1997.0202) 6.        Alexander RMN, Jayes AS. 1983 A dynamic similarity hypothesis for the gaits of quadrupedal mammals. <i>J. Zool.</i> <b>201</b>, 135–152. (doi:10.1111/j.1469-7998.1983.tb04266.x) 7.        Irschick DJ, Jayne BC. 1999 Comparative three-dimensional kinematics of the hindlimb for high-speed bipedal and quadrupedal locomotion of lizards. <i>J. Exp. Biol.</i> <b>202</b>, 1047–1065. 8.        Bates D, Mächler M, Bolker BM, Walker SC. 2015 Fitting linear mixed-effects models using lme4. <i>J. Stat. Softw.</i> <b>67</b>. (doi:10.18637/jss.v067.i01) 9.        Kuznetsova A, Brockhoff PB, Christensen RHB. 2017 lmerTest Package: Tests in Linear Mixed Effects Models . <i>J. Stat. Softw.</i> <b>82</b>. (doi:10.18637/jss.v082.i13) 10.      Long J. 2019 interactions: Comprehensive, user-friendly toolkit for probing interactions. 11.      Paradis E, Schliep K. 2019 Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. <i>Bioinformatics</i> <b>35</b>, 526–528. (doi:10.1093/bioinformatics/bty633) 12.      Clemente CJ. 2014 The evolution of bipedal running in lizards suggests a consequential origin may be exploited in later lineages. <i>Evolution (N. Y).</i> <b>68</b>, 2171–2183. (doi:10.1111/evo.12447)