Ontogeny of joint mechanics in squirrel monkeys (Saimiri boliviensis): functional implications for mammalian limb growth and locomotor development

Juvenile animals must often compete against adults for common resources,keep pace during group travel and evade common predators, despite ontogenetic limits on locomotor performance, such as reduced body size and an immature musculoskeletal system (Carrier,1996). Because, by definition, juveniles have yet to reproduce,there should be strong selection for compensatory mechanisms, such as allometric musculoskeletal growth trajectories or developmental changes in behavior, that could enhance locomotor performance despite size- and growth-related limitations on performance(Pennycuick, 1975; Carrier, 1983; Werner and Gillam, 1984; Carrier and Leon, 1990; Irschick, 2000; Trillmich et al., 2003; Main and Biewener, 2004; Irschick et al., 2005; Young, 2005; Herrel and Gibb, 2006; Lawler, 2006). This study tests the hypothesis that allometric changes in joint mechanics allow growing squirrel monkeys to compensate for ontogenetic limits on locomotor performance.

Carrier (Carrier, 1983)documented growth trajectories in black-tailed jackrabbits (Lepus californicus) and domestic cats (Felis domesticus) that theoretically permit these species to overcome ontogenetic limits on locomotor performance. In these animals, the length of the olecranon process at the elbow and the calcaneal tuberosity at the ankle, proportional to the average moment arms of m. triceps brachii and m. triceps surae, respectively, scale with negative allometry relative to the length of the anatomical segments distal to these joints (i.e. these muscles act on a relatively longer lever arm in younger animals). In a separate study of limb growth and locomotor development in domestic cats, Peters(Peters, 1983) corroborated Carrier's (Carrier, 1983)findings. Assuming that SRF load arms are proportional to limb length and SRF magnitudes are proportional to body weight, negative allometry of anatomical lever arms relative to limb length suggests that muscle mechanical advantage(i.e. r/R) should be greater in young animals, whereas requisite muscle should be reduced, perhaps compensating for limited muscle mass and smaller body size. In a subsequent study, Carrier(Carrier, 1995) found that young jackrabbits were indeed able to achieve adult-like jumping velocities at only 30% of adult body size, thus confirming the link between allometric growth trajectories and enhanced locomotor performance at young ages. Carrier(Carrier, 1983) suggested that because identical scaling trends were found in distantly related jackrabbits and cats, negative allometry of muscle mechanical advantage might represent a generalized mammalian solution to ontogenetic limits on locomotor performance. In support of Carrier's (Carrier,1983) hypothesis, Young(Young, 2005) recently documented a similar pattern of negative ontogenetic scaling of m. triceps brachii mechanical advantage in arboreal capuchin monkeys (Cebus apella and Cebus albifrons).

Carrier's (Carrier, 1983)and Young's (Young, 2005)analyses assumed that joint postures did not change during development and that SRF magnitudes remained a constant multiple of body weight. In other words, for the negative allometry of bony muscle lever arms to have the assumed effects, joint load arm lengths and joint moments should scale isometrically to body size. [For simplicity, external (i.e. SRF) joint moments will be referred to simply as joint moments throughout this paper. When referring to internal moments, the term muscle moment will be used.] If,however, developmental changes in joint postures and SRF magnitudes cause joint moments to scale allometrically, muscle lever growth trajectories may be more difficult to interpret. If joint moments are relatively greater in young animals, relatively longer bony muscle lever arms may not be sufficient to achieve adult-like levels of locomotor performance. Conversely, if joint moments are relatively smaller at early ages and increase over development,young mammals may not require relatively longer bony levers to increase performance.

This study focused on locomotor development in Bolivian (i.e. black-capped)squirrel monkeys [Platyrrhini: Saimiri boliviensis (Geoffroy and Blainville 1834)]. Squirrel monkeys are among the smallest anthropoid primates, with an average adult body mass of 811 g(Smith and Jungers, 1997). As a result, most squirrel monkey populations are under intense predation pressure (Fedigan et al.,1996). Predation risk has profoundly affected many aspects of squirrel monkey biology, including growth and development(Boinski, 1987; Boinski, 1999; Boinski et al., 2003). In a study of red squirrel monkey (Saimiri oerstedi) behavioral ecology in Costa Rica, Boinski (Boinski,1987) found that infants experienced higher rates of predation than any other age class, with more than 50% of infants dying within the first six months of life. Similar rates of predation have been reported for S. boliviensis in Peru (Boinski et al.,2002). As a means of coping with predation risk early in life,squirrel monkey behavioral development is markedly precocial relative to that of other primates (Elias,1977; Hartwig,1995; Garber and Leigh,1997). In captivity, infants engage in independent locomotion within the first month of life (Elias,1977; Kaack et al.,1979; Fragaszy et al.,1991). In the wild, infants begin foraging independently within the first or second month of life, are traveling primarily independently by three months of age and are weaned between four and eight months of age(Boinski and Fragaszy, 1989; Mitchell, 1990; Stone, 2006). By six months of age, when juveniles are only 40–50% of adult size, foraging activity and locomotor repertoires are generally indistinguishable from adults(Boinski, 1989; Boinski and Fragaszy, 1989). Moreover, once juvenile squirrel monkeys have become independent, they must travel an average of 2–4 km day^–1 to remain with the group and gain access to distributed foraging resources, such as fruit and invertebrate prey (Terborgh,1983; Mitchell,1990; Boinski,1999).

In summary, squirrel monkeys are an excellent group in which to investigate how developmental changes in limb mechanics might impact locomotor performance in young mammals. In this study, I used standard kinematic and kinetic techniques to investigate the ontogenetic scaling of joint postures, SRF magnitudes and SRF joint moments in a longitudinal sample of growing Bolivian squirrel monkeys. Because squirrel monkeys use a variety of gaits across a diversity of substrates (Fontaine,1990), I examined both symmetrical gaits (i.e. walking and running) and asymmetrical gaits (i.e. galloping and bounding) on terrestrial and simulated arboreal substrates (i.e. a flat runway and an 2.5 cm diameter elevated pole, respectively). Flat runway locomotion was sampled in order to remain consistent with previous studies of mammalian joint mechanics (e.g. Biewener, 1983; Schmitt, 1998; Polk, 2002; Witte et al., 2002) whereas locomotion on the pole was intended to represent a more naturalistic environment for arboreal squirrel monkeys. Additionally, in their natural environment, squirrel monkeys travel upon a variety of substrates that vary widely in diameter (Terborgh,1983; Boinski,1989; Fontaine,1990; Mitchell,1990). Sampling locomotion on the ground and a relatively narrow pole should represent the extremes of this variation in preferred substrate diameter. Because previous studies have assumed that limb growth is the primary determinant of ontogenetic variation in joint load arm lengths (e.g. Carrier, 1983; Young, 2005), I also measured limb segments (i.e. arm, forearm, thigh and leg lengths) at regular intervals and used path analysis (Li,1975) and hierarchical partitioning(Chevan and Sutherland, 1991)to investigate the independent influence of limb length and posture on developmental changes in joint load arm lengths.

I tested two null hypotheses, based on the morphometric findings of Carrier(Carrier, 1983) and Young(Young, 2005). First, when growing squirrel monkeys are traveling at the same absolute speed, joint angles, SRF magnitudes, SRF angles, joint load arm lengths and joint moments should scale isometrically to body size. Specifically, angular variables should remain constant throughout ontogeny, SRF magnitudes should scale to body weight, joint load arm lengths should scale to the cube root of body mass, and joint moments should scale to body mass. Second, I predicted that ontogenetic variation in joint load arm lengths should be primarily determined by ontogenetic variation in limb length distal to the joint in question. To examine the independent and overall influence of limb length on joint load arm length, I constructed path models (Li,1975) specifying the hypothesized influence of limb length, joint posture, limb posture and SRF orientation on joint load arm length. The general organization and predicted form of the path models are presented in Fig. 1.

Five female infant–juvenile squirrel monkeys constituted the sample for this study. An additional squirrel monkey participated in two experiments but had to be withdrawn because of an unrelated tail injury. Collectively, the squirrel monkeys ranged in age from 74 to 302 days, and in body mass from 218 to 535 g (29–71% of adult size; mean: 367 g). Participation duration for individual monkeys ranged from 151 to 185 days (mean: 172.2 days; median: 179 days), beginning at ages of 74–120 days (mean: 104.6 days; median: 118 days). The age range sampled in this study (i.e. 74–302 days) therefore encompasses the period during which wild squirrel monkeys begin locomoting independently, gain foraging competence and become independent juveniles. Monkeys lived in social groups of 15–30 individuals that were housed in large 1.5 m×2.1 m×4.5 m enclosures with 35–40 linear meters of perches and substrates for free-ranging locomotor activity. Research was performed at the Center for Neotropical Primate Research and Resources (CNPRR:Mobile, AL, USA). Institutional Animal Care and Use Committees (IACUC) at Stony Brook University and the CNPRR approved all procedures before the start of this research.

Before each experiment, individuals were weighed and the skin over the approximate centers of rotation of the shoulder, elbow, wrist, hip, knee and ankle were shaved and marked with retro-reflective tape (3 M Corporation, St Paul, MN, USA), a procedure that did not require the use of anesthetic. Anatomical landmarks used to identify approximate joint centers of rotation are listed in Table 1. Arm,forearm, thigh and leg lengths, measured as the distance between adjacent landmarks, were also recorded at this time. Monkeys were then filmed as they traversed a 2.75 m enclosed linear runway using a high-speed digital video camera (MotionMeter 1000, Redlake MASD, San Diego, CA, USA) operating at 250 Hz. Depending on experimental condition (e.g. ground or pole), monkeys traversed either the flat runway floor or a 2.5 cm diameter PVC pipe elevated 10.7 cm above the surface of the runway. Speed was not controlled during locomotion. Rather, at all ages, animals were permitted to self-select preferred speeds.

Two custom-built 30.5 cm×30.5 cm triaxial force platforms(Heglund, 1981; Biewener and Full, 1992),placed in series in the center of the runway, were used to measure locomotor kinetics. During the simulated arboreal trials, PVC segments (30.5 cm long and 2.5 cm in diameter) were attached to each force platform by bolts secured directly to the platform frame. Voltage outputs from channels corresponding to each force axis were routed through a National Instruments (Austin, TX, USA)SC-2345 chassis and recorded using a LabView virtual instrument. Cross talk between force channels was generally low, ranging between 0.3% and 3.5%without the pole segment and 0.3% and 10.5% when the pole segment was attached. Force platforms were calibrated daily following the recommendations of Biewener and Full (Biewener and Full,1992).

Following the method of Riskin et al.(Riskin et al., 2005), kinetic and kinematic data were synchronized using a 3.3 V square-wave pulse generated by the video camera and routed separately to a bank of LEDs positioned on the back wall of the runway and to the SC-2345 chassis. This circuit was normally interrupted by means of a handheld switch. During each trial, the switch was briefly closed, simultaneously illuminating the LEDs and changing the shape of the square wave in the data file. Using this procedure, it was possible to synchronize video and kinetic data to a resolution of four milliseconds.

Video files from each experiment were imported into the MATLAB DLT Dataviewer 2 digitizing platform (Hedrick,2007) for coding on a trial by trial basis. Individual strides were identified based on the cyclic touchdowns of a reference limb. For steps(i.e. stance phases) in which single-limb kinetic data were available, the two dimensional position of all limb landmarks and the fifth metapodial head were digitized at peak vertical force. Owing to the small size of the young monkeys, reflective markers could not be attached to the skin overlying the metacarpal and metatarsal heads. Therefore, the position of these landmarks was estimated to be at the base of the fifth manual/pedal ray.

Force data from each trial were imported into MATLAB (Mathworks, Natick,MA, USA), transformed into Newtons and corrected for cross talk. Force traces from each channel were smoothed using a zero-lag fourth-order Butterworth low-pass filter with a cut-off frequency of 25 Hz. Baseline drift during and between trials was corrected by sampling the average values of unloaded periods immediately before and following platform contact and subtracting these values from the force traces. Overall, baseline drift was extremely low,averaging less than 3 mN s^–1 across force channels and substrate conditions (ground: vertical=2.7±2.84 mN s^–1, fore–aft=2.4±6.10 mN s^–1;mediolateral=1.1±0.97 mN s^–1; pole:vertical=2.0±1.52 mN s^–1,fore–aft=1.1±1.28 mN s^–1;mediolateral=0.9±1.80 mN s^–1). On average, these values equated to less than 0.3% of the maximum force recorded during a given trial. Because step durations averaged 0.13±0.066 s across the dataset,the maximum average amount of drift during a single contact period with the force plates would have been 0.35 mN. Baseline drift thus exerted a negligible effect on force platform accuracy and precision.

Single-limb contacts were recognized when a forelimb or hind limb contacted the force platform in isolation and ended when the limb either left the platform or another limb touched down. Because the monkeys frequently placed limbs in close proximity during locomotion, obtaining isolated single-limb contacts was often difficult. Trials with overlapping limb contacts on a single force platform were retained only when peak vertical force was clearly identifiable and the vertical force trace had returned to 50% of the peak value prior to rising again in the case of forelimb contacts, or began its rise at no more than 50% of the former vertical force peak in the case of hind limb contacts. An exemplar stride illustrating the experimental apparatus used for kinematic and kinetic data collection is shown in supplementary movie 1(see supplementary movie).

Average locomotor speed was calculated from the sagittal displacement of either the shoulder or the hip, depending on marker visibility. Examination of trials in which both hip and shoulder data were available revealed that hip and shoulder speed varied on average by 0.02±0.054 m s^–1, or approximately 1.1% of average forward speed. Digitizing noise was corrected using a zero-lag fourth-order Butterworth low-pass filter with a cut-off frequency of 10 Hz. This cutoff frequency was selected as optimal using a residual analysis procedure described by Winter(Winter, 2005). Piecewise cubic spline interpolation was used to interpolate over gaps of missing data≤10 frames (i.e. 40 ms). After transforming raw pixel coordinates into meters by using the force platforms as calibration objects, linear least-squares regressions of corrected displacement data on time were used to calculate overall speed across each stride. Trials in which the coefficient of determination (i.e. R²) of reference marker position against time was less than 0.99 were designated unsteady and discarded.

To distinguish between symmetrical and asymmetrical strides, stride symmetry was calculated as the absolute duration between the touchdowns of contralateral forelimb and hind limb pairs (e.g. right and left forelimbs)expressed as a percentage of total stride duration. Following Hildebrand(Hildebrand, 1976), strides in which both forelimb and hind limb symmetry was between 43.75% and 56.25% were designated symmetrical; all other strides were designated asymmetrical. Data from leading and trailing limbs within asymmetrical gaits were treated separately for allometric analyses of joint mechanics.

Segment angles (i.e. forelimb, arm, hind limb and thigh angles) were calculated as the two-dimensional vector angle between the relevant limb segment and the vertical axis. Forelimb and hind limb segments were defined by a line joining the shoulder or hip to the metapodial head. Joint angles (e.g. elbow or knee angle) were calculated as the two-dimensional vector angle between the relevant limbs segments (e.g. arm and forearm or thigh and leg).

Because only sagittal (i.e. two-dimensional) kinematic data were recorded,the mediolateral component of SRF was not included in the above calculations. However, peak mediolateral force magnitudes were consistently minor (i.e. approximately 5% of peak vertical force magnitudes) and therefore had little effect on the magnitude or resultant orientation of the SRF. Following previous studies of mammalian joint mechanics (e.g. Biewener, 1983; Schmitt, 1998; Polk, 2002; Witte et al., 2002), SRFs were assumed to pass through the metapodial heads at the moment of peak vertical SRF. Joint moments were calculated as the cross (i.e. vector) product of the SRF and the two-dimensional position vector connecting the metacarpal/metatarsal head to the joint in question(Ozkaya and Nordin, 1999). Load arm lengths were computed by dividing joint moments by the SRF magnitude(Polk, 2001). Positive values signify flexing moments arms (dorsiflexing at the wrist and ankle), whereas negative values signify extending moments (palmar- and plantarflexing at the wrist and ankle).

It should be noted that the methods used here did not account for the effects of limb segment inertia or gravity on calculated joint moments. However, the present analysis was explicitly concerned with the effects of ontogenetic variation in limb posture and SRF on joint mechanics. The effects of ontogenetic changes in limb inertial properties is a related, but separate,issue that should certainly be examined in additional studies (e.g. Raichlen, 2006). Moreover, in the small-bodied quadrupeds examined here, inertial forces are likely to be minor relative to SRF (Biewener and Full,1992). Indeed, Witte et al.(Witte et al., 2002) have recently shown that inertial forces account for no more than 10% of total limb joint moments among small quadrupedal mammals (i.e. 150–400 g in body mass).

To adjust for ontogenetic differences in body size, SRF magnitudes were divided by body weight, load arm lengths were divided by the cube root of body mass and joint moments were divided by body mass, following Witte et al.(Witte et al., 2002). Scaling joint moments to body mass follows from the isometric expectation that muscle lever arm lengths, as linear dimensions, scale to the cube root of body mass(i.e. ∞M_b^0.33), whereas muscle force,proportional to physiological cross-sectional area, scales to the two-thirds power of body mass (i.e. ∞M_b^0.67). Muscle moments, equal to the product of force and length, should therefore scale directly to body mass (i.e. M_b^0.33 × M_b^0.67). Mass-adjusted joint moments therefore provide a metric of load magnitude relative to theoretically available muscle moments. This scaling procedure necessarily assumes that muscle moments scale isometrically to body mass in growing squirrel monkeys, a proposition that remains to be tested. Nevertheless, this assumption is consistent with previous morphometric studies of developmental joint mechanics (i.e. Carrier, 1996; Young, 2005), where it is argued that increased muscle mechanical advantage would be required to compensate for reduced muscularity at young ages and isometrically maintain similar mass-adjusted muscle moments throughout ontogeny.

To increase statistical power, data from individual monkeys were combined to create mixed longitudinal samples for all analyses. Body mass, which was strongly positively correlated with age across individuals(Fig. 2), was used as the primary independent variable in all ontogenetic analyses. Moreover, because of the precocial nature of squirrel monkey behavioral development(Elias, 1977; Hartwig, 1995), size should have a stronger effect on locomotor mechanics than age per se (cf. Schilling, 2005).

Morphometric data were tested for allometric growth by fitting log-transformed segment lengths and body masses to the standard allometric power function (Huxley, 1932)using Model II reduced major axis (RMA) regression(Ricker, 1984). Confidence intervals on RMA slopes were calculated following Pitman(Pitman, 1939). Allometry was identified when the 95% confidence intervals for the calculated slopes did not include 0.333. Regressions were performed using the (S)MATR software package(Falster et al., 2003).

To characterize intra-limb variability in joint loading across the ontogenetic sample, forelimb and hind limb joint moments were specified as within-subjects factors in repeated-measures analyses of variance (ANOVA). In cases where limb joint moments varied as a linear function of body mass and/or speed, repeated-measures analyses of covariance (ANCOVA) were used in lieu of ANOVA, specifying mass and/or speed as the covariate(s). Paired t-tests between group means (least-squares adjusted means following ANCOVA) were used to test for significant post-hoc differences between joints. Analyses were performed separately by substrate, gait type and limb order (in the case of asymmetrical gaits) using SPSS 11.0.4 (SPSS,Chicago, IL, USA).

Partial correlations were used to investigate the influence of body mass on each variable while controlling for speed. The null hypothesis for each test was that body mass had no association with the dependent variable. Because all variables were either dimensionless (i.e. angles) or were adjusted for body size differences prior to analysis, any significant correlation indicated allometry (Mosimann and James,1979). Correlation analyses were again performed separately by gait type, substrate and limb order (in the case of asymmetrical strides)using SPSS 11.0.4.

Hierarchical partitioning (Chevan and Sutherland, 1991) and path analysis(Li, 1975; Sokal and Rohlf, 1995) were used to dissect the independent influence of limb length and posture on observed variation in absolute joint load arm lengths. Both hierarchical partitioning and path analysis can be thought of as extensions of multiple regression. Given a linear model with one dependent variable (Y) and k independent variables (X₁, X₂, X₃...X_k),hierarchical partitioning quantifies the independent contribution of each X variable to the total coefficient of determination for Yas the average change in R² produced by adding the variable to a hierarchical series of increasingly complex models. The significance of each variable's independent contribution to the total coefficient of determination can be tested using a randomization procedure introduced by MacNally (MacNally,2002). Hierarchical partitioning was performed using the hier.part package (MacNally and Walsh,2004) of the R statistical platform(R Development Core Team,2008).

Path analysis (Li, 1975; Sokal and Rohlf, 1995) was used to clarify the independent and associative influence of limb lengths on load arm lengths. In a path analysis, each hypothesized causal link (i.e. direct association) between variables is represented diagrammatically by a single-headed arrow and is associated with a path coefficient –equivalent to a standardized partial regression coefficient in a multiple regression. Additionally, predictor variables can be joined viaundirected double-headed arrows (equivalent to a standard bivariate correlation), signifying that the variables are related to one another without specifying a necessary causal structure. Once a path diagram has been constructed, and all associated path coefficients have been specified, the`indirect' association between any predictor and a criterion variable can be quantified as the summation of all the products of coefficients linking the two variables (Sokal and Rohlf,1995).

Path models were fit to the data using iterative maximum likelihood estimation (Schermelleh-Engel et al.,2003). Each path model was initially fit with all of the hypothesized paths included in the model. Non-significant paths were then removed from the model specification and the model was refit, resulting in the most parsimonious model for the data. The overall fit of each path model was evaluated by calculating χ² tests that compared the empirically observed correlation matrix to the correlation matrix implied by the model. Non-significant χ² tests (i.e. P>0.05) indicate a good fit between the model and the data(Schermelleh-Engel et al.,2003). All path analyses and associated tests of fit were performed using the AMOS 16 software package (AMOS Development Corporation,Spring House, PA, USA). Only trials without missing data for any of the parameters were included in the path analyses.

Because the goal of hierarchical partitioning and path analyses was to test the hypothesis that joint load arms were directly proportional to anatomical limb length distal to the joint in question, only those joints for which the relevant limb lengths were available (i.e. shoulder, elbow, hip and knee joints) were included in these analyses. Additionally, because path analysis requires large sample sizes to obtain effective power (typically ≥10 times the number of parameters included in the model)(Schermelleh-Engel et al.,2003), data were pooled across gaits within substrates. Previous analyses of this data set indicated that the path models were statistically similar across gait types within substrates. Specifically, significant differences between path coefficients from symmetrical versusasymmetrical path models within substrates were found in only three out of 72 possible cases (i.e. nine parameters per path model × four joints× two substrates).

A total of 664 strides were analyzed. The tabulation of strides and forelimb and hind limb steps by substrate, gait type and limb order (for asymmetrical strides) is shown in Table 2.

In the sample of squirrel monkeys measured here, all segment lengths scaled to body mass with strong positive allometry(Fig. 3), indicating that smaller and younger monkeys were relatively short-limbed for their body size. Forelimb segments tended to scale with greater positive allometry than hind limb segments, although only the forearm scaled significantly faster than hind limb segments (P≤0.01 for all comparisons). Scaling exponents describing overall forelimb and hind limb growth indicated strong positive allometry for both limbs. Total forelimb length scaled with significantly greater positive allometry than total hind limb length(P<0.05).

Across the ontogenetic sample, most forelimb and hind limb joints primarily experienced flexing (i.e. positive) joint moments at peak vertical force(Fig. 4), regardless of gait type or substrate. The only consistent exception was the knee joint, which experienced a combination of weakly flexing and extending moments. The hip and shoulder joints also occasionally experienced extending moments, particularly during asymmetrical gaits. Forelimb joint moment magnitudes were consistently greatest at the elbow, regardless of substrate, gait or limb order (all P<0.001). Wrist moments were significantly greater than shoulder moments during all locomotion on the ground (all P<0.05), whereas shoulder and wrist moments were statistically similar during locomotion on the pole. Hind limb joint moment magnitudes were consistently smallest at the knee(all P<0.05). Hip moments were significantly greater than ankle moments during symmetrical gaits (all P<0.001), regardless of substrate, and for leading limbs during asymmetrical gaits on the ground(P<0.01). In all other conditions, hip and ankle moments were statistically similar.

Absolute locomotor speed and body mass were significantly positively correlated during asymmetrical gaits on the pole, but were unrelated during all other substrate-by-gait conditions(Table 3). Relative speed [i.e. Froude number: u(gh)^–0.5,where u is speed, g is gravitational acceleration(9.81 m s^–2) and h is the average of hip and shoulder heights at peak vertical force](Biewener, 2003) also changed little with speed (Table 3). As with absolute speed, relative speed was significantly positively correlated with body mass during asymmetrical gaits on the pole, but was unrelated to body mass during all other conditions. Overall, these data indicate that squirrel monkeys used a wide variety of speeds across the ontogenetic series and that the youngest monkeys frequently moved as fast as older and larger monkeys.

Partial correlations describing the independent association between body mass and all other kinematic/kinetic parameters, controlling for absolute speed, are shown in Table 3. Previous analyses of this dataset (Young, 2008a) indicated that the patterns shown in Table 3 remained consistent when partial correlations were calculated controlling for relative,rather than absolute, speed. Specifically, out of the 144 correlations displayed in Table 3 (i.e. 24 parameters × six substrate-by-gait conditions), significant differences between absolute speed-adjusted and relative speed-adjusted correlations were only found in seven cases.

Forelimb angle changed minimally with body mass, becoming significantly more protracted during symmetrical gaits on the ground and for leading limbs during asymmetrical gaits on the pole, but showing no other significant changes as body size increased. By contrast, elbow angle became significantly more flexed with increasing mass in all conditions except for symmetrical gaits on the ground and for leading limbs during asymmetrical gaits on the ground. Changes in arm angle principally tracked changes in elbow angle,becoming significantly more retracted during symmetrical gaits on the pole and for trailing limbs during asymmetrical gaits on the pole. Wrist angle did not change with size during symmetrical gaits, but became significantly more flexed for leading limbs during asymmetrical gaits on the ground and significantly more extended in all other asymmetrical gait conditions.

As with forelimb angle, hind limb angle changed little with increasing body mass, becoming significantly more protracted during symmetrical gaits on the ground and for leading limbs during asymmetrical gaits on the pole, but remaining static in all other conditions. Knee angle became significantly more flexed during symmetrical gaits on the pole but did not significantly change in any other condition. Size-related changes in thigh angle tracked changes in knee and hind limb angles, becoming significantly more flexed (i.e. protracted) during symmetrical gaits on the pole and for leading limbs during asymmetrical gaits on the pole. Ankle angle did not significantly vary with body size in any condition.

Relative forelimb force magnitudes significantly decreased with increasing body size across all conditions except for leading limbs during asymmetrical gaits on the ground. By contrast, relative hind limb peak force magnitudes significantly increased during symmetrical gaits on both substrates and for leading limbs during asymmetrical gaits on the pole, and showed non-significant trends to increase in all other conditions except for trailing limbs during asymmetrical gaits on the pole. Forelimb force angles significantly decreased with size (i.e. became significantly more caudally oriented) during symmetrical gaits on the ground and for leading limbs during asymmetrical gaits on the pole. Similarly, hind limb force angles significantly decreased for leading limbs during asymmetrical gaits on the pole. In each of these cases, caudal reorientation of the SRF vector corresponded to significant increases in limb protraction with increasing size. Similarly, forelimb force angles for leading limbs during asymmetrical gaits on the ground became significantly more cranially oriented as size increased, corresponding to a nearly significant increase in forelimb retraction (P=0.066). Altogether, these results suggest that both the fore- and hind limbs were operating at least partly as struts (i.e. the SRF vector was oriented along the effective axis of the limb)(Barclay, 1953; Gray, 1968).

Variation in relative joint load arm lengths principally tracked variation in joint angles. Increases in elbow flexion corresponded to significant size-related increases in relative elbow load arm lengths during symmetrical gaits on the pole, for leading limbs during asymmetrical gaits on the pole,and for trailing limbs during asymmetrical gaits on both substrates. Similarly, increases in wrist extension corresponded to decreases in relative wrist load arm lengths during asymmetrical gaits on the pole, regardless of limb order. Among the hind limb joints, significant changes were only observed during symmetrical gaits on the ground, where hip and ankle load arms became longer with increasing body size. By contrast, knee load arms decreased during symmetrical gaits on the ground, in fact passing from the flexing to the extending side of the joint (Fig. 4).

Significant size-related increases in relative moments were observed at several joints in multiple conditions. In the forelimb, relative elbow moments increased during symmetrical gaits on the pole, in leading forelimbs during asymmetrical gaits on the pole, and in trailing forelimbs during asymmetrical gaits on both substrates. In the hind limb, relative hip and ankle moments increased during all conditions except for trailing limbs during asymmetrical gaits on the ground. Size-related increases in relative knee moments were observed for trailing limbs during asymmetrical gaits on the ground. The only significant size-related decrease in relative moments was at the wrist joint for trailing limbs during asymmetrical gaits on the pole. Relative moments at all other joints and conditions remained static with increasing body size.

The results of hierarchical partitioning and path analyses of the influence of limb length and posture on shoulder, elbow, hip and knee load arm lengths are presented graphically in Fig. 5 for strides on the ground and Fig. 6 for strides on the pole. The direct, indirect and total effects of limb lengths on load arm lengths, as implied by the path models, are shown in Table 4. Overall, the path models described here predicted patterns of correlation among limb lengths,joint postures, limb postures, SRF angles and load arm lengths very well, as indicated by non-significant χ² statistics in all cases (all P≥0.23).

Variation in forelimb length, arm angle, forelimb angle and forelimb SRF angle collectively explained 97.6–97.8% of the variance in shoulder load arm lengths across the ontogenetic sample. Forelimb angle and SRF angle were consistently the best predictors of shoulder load arm length, irrespective of substrate. Together, these two variables explained more than 90% of the variation in shoulder load arm length on both substrates. Forelimb and SRF angles were negatively correlated, as expected if the limbs were at least partly operating as struts (Barclay,1953; Gray, 1968). Increases in forelimb angle (i.e. greater protraction), were associated with shorter shoulder joint load arms, even when controlling for SRF orientation,supporting the prediction that a more protracted forelimb places the shoulder joint more in line with the caudally inclined SRF vector, thereby mitigating load arms and moments.

Forelimb length explained just 2.2–3.6% of the variance in shoulder load arm length. Across substrates, most of the total correlation between forelimb length and shoulder load arm length was due to the direct effect of forelimb length (Table 4). Negative indirect correlations between forelimb length and shoulder load arm length, operating via arm and forelimb angle, slightly mitigated the total correlation between forelimb length and shoulder load arm lengths on both substrates. In sum, longer forelimbs were weakly associated with increases in shoulder joint load arms lengths, and this relationship was made even weaker when the effects of limb length on joint posture were included in the model.

Variation in forearm lengths, elbow angles, forelimb angles and forelimb SRF angles explained 85.7–87.5% of the variance in elbow load arm lengths. As predicted (Fig. 1),increases in forearm length, elbow flexion, forelimb protraction and forelimb SRF angle were all associated with longer elbow load arms. Across substrates,elbow joint angle was the best predictor, explaining 32.6–44.7% of the variance. In both conditions, forearm length followed elbow posture as the next best predictor of elbow load arm length, explaining 18.2–24.4% of the load arm length variance. Together, forearm length and elbow angle explained more than 50% of the variance in elbow load arm lengths on both substrates.

On the ground, forearm length was unassociated with elbow posture and only weakly associated with forelimb posture. Therefore, most of the total correlation between forearm length and elbow load arm length was due to the direct effects of forearm length (Table 4). By contrast, because forearm length was relatively strongly associated with elbow posture during locomotion on the pole, forearm length exerted a positive direct effect on elbow load arm length as well as a equally pronounced positive indirect effect because of the tendency of longer-limbed monkeys to use more flexed elbow postures. Overall, the total association between limb length and load arm length was greater at the elbow than at any other joint.

Together, hind limb length, thigh angle, hind limb angle and hind limb SRF angle explained 98.4–98.8% of the variation in hip load arm lengths across the ontogenetic sample. As predicted(Fig. 1), longer hip load arms were associated with greater hind limb protraction at peak vertical force and a more cranial inclination of the hind limb SRF. Additionally, the predicted relationship between thigh angles and hip load arms was upheld, with more extended thigh positions exerting a negative effect on hip load arm length. Hind limb angles were consistently the best predictors of hip load arm length,explaining approximately 50% of the variation in hip load arm length across substrates. On the ground, hind limb SRF angle was the second best predictor of hip load arm length, explaining 21.6% of the variance. On the pole, the second best predictors were thigh angle and hind limb SRF angle, explaining 28.3% and 23.1% of the variance, respectively.

As predicted, increases in total hind limb length were associated with longer hip joint load arms. However, the explanatory power of hind limb length was generally weak. On the pole, hind limb length was unable to explain a significant portion of load arm variance. Hind limb length performed better on the ground, explaining 9% of the variation in hip load arm length. On the ground, most of the total correlation between hind limb length and hip load arm length was attributable to the indirect influence of limb length on hind limb posture. On the pole, longer-limbed individuals tended to use more extended thigh postures, thus shortening hip load arms, mitigating the slight positive direct effect of hind limb length and reducing the total correlation between hind limb length and hip load arm length(Table 4).

Leg length, knee angle, hind limb angle and hind limb SRF angle collectively explained 96.8–97.6% of the variance in knee load arm lengths. On the ground, hind limb and knee angles were the best predictors of knee load arm lengths, explaining 41.7% and 39.1% of the variance,respectively. On the pole, hind limb SRF angle was the best predictor of knee load arm length, explaining 59% of the variance. Knee joint posture was the second best predictor, explaining 20.2% of the variance. As predicted(Fig. 1), increases in knee load arm lengths were associated with greater hind limb retraction, a more caudal orientation of the hind limb SRF and more flexed knee postures.

Leg length was consistently a poor predictor of knee load arm length,explaining just 1.4–4.2% of the variance. On the ground, a negative indirect correlation operating via hind limb angle negated the direct effects of leg length on knee joint load arm length, leading to a weak total correlation. Longer-legged individuals tended to use more flexed knee postures on the pole, and most of the total correlation between leg length and knee load arm length could be attributed to this indirect effect operating via knee joint posture.

Squirrel monkey limb growth during the first 10 months of life was characterized by extremely strong positive allometry, indicating that younger monkeys have relatively short limbs for their size(Fig. 3). All else being equal,relatively short limbs should reduce the distance between joint centers and reaction force vectors, shortening load arms and mitigating joint loading in younger animals. However, hierarchical partitioning and path analyses of the current dataset indicated that variation in anatomical limb length determined only a minor proportion of the variance in load arm lengths at most joints(i.e. less than 10% of the variance at all joints save the elbow; see Figs 5 and 6). Rather, across substrates,load arm lengths primarily varied as a function of posture and SRF angle. Together, joint angle, limb angle and SRF angle accounted for 89.5–97.7%of load arm length variance across all joints except the elbow. The tight correlation between load arm lengths and limb/SRF angles clarifies the predominant lack of ontogenetic variation in most load arm lengths(Table 3), despite significant increases in relative limb length. With the exception of persistent size-related changes in elbow and wrist postures, joint kinematics remained mostly static across the ontogenetic series. Owing to the strut-like behavior of the fore- and hind limbs, SRF angles were strongly associated with limb postures, and therefore ontogenetic variation in SRF orientation was similarly restricted. As a result, size-related changes in relative load arm lengths were primarily restricted to the elbow and, to a lesser extent, the wrist. Only during symmetrical gaits on the ground were ontogenetic changes in limb angles and SRF angles, specifically in the hind limb, sufficient to impact the relative length of load arms during growth.

In contrast to the other joints examined, forearm length alone accounted for 18.2–24.4% of elbow load arm length variance, and forearm length in combination with elbow angle accounted for 57–62.9% of the variance. The relatively strong dependency of elbow load arm length on forearm length is probably the result of two factors: limb geometry and hand posture at peak vertical force. First, the habitually pronated position of the mammalian forearm during quadrupedal locomotion ensures that the palm, and therefore the estimated hand center of pressure (COP), will be relatively distant from the elbow joint at peak vertical force (Fig. 1E). By contrast, typical hind limb postures ensure that the knee joint and the estimated foot COP could very well be collinear along the axis of the hind limb SRF vector. Second, the metatarsals were held in a significantly more elevated (i.e. digitigrade) position than the metacarpals at peak vertical force (Fig. 7)[P<0.001: Watson's two-sample test for circular data(Batschelet, 1981)], further increasing the distance between the estimated hand COP and the elbow joint. Altogether, the increased distance between the estimated hand COP and the elbow permitted more variation in limb length and elbow posture to be translated into load arm length. As such, ontogenetic increases in relative forearm length (Fig. 3) and elbow flexion caused significant size-related increases in relative elbow load arm lengths during most substrate-by-gait conditions(Table 3).

Whereas elbow postures became increasingly flexed as size increased ontogenetically, the wrist joint tended to become more extended (although significant increases in flexion were observed for leading limbs during asymmetrical gaits on the ground). As with the elbow joint, forearm pronation and sub-digitigrade postures at peak vertical force probably increased the distance between the forelimb SRF vector and the wrist joint, permitting more of the variation in wrist angle to be translated into load arm length variation. As a result, ontogenetic increases in wrist extension led to significant decreases in relative wrist load arm lengths during asymmetrical gaits on the pole. Additionally, although hand growth was not measured in the current study, negative ontogenetic allometry of hand length during development, as has been observed in numerous other mammalian taxa (see below; Fig. 8), may have contributed to the observed decrease in wrist load arm lengths.

Relative forelimb SRF magnitudes declined as body size increased in almost all substrate and gait conditions. By contrast, relative hind limb SRF magnitudes tended to increase with body size, although this trend was only significant in half of the substrate-by-gait conditions. Ontogenetic shifts in forelimb–hind limb peak force distribution were associated with a significant caudal translation in whole-body COM position(Young, 2008), a growth-related phenomenon has also been observed in other primates(Turnquist and Wells, 1994; Crompton et al., 1996; Shapiro and Raichlen,2007).

Ontogenetic shifts in posture and limb force distribution led to several significant size-related increases in relative moments at the elbow, hip and ankle joints. Because relative forelimb SRF magnitudes largely declined over development, increases in relative elbow moments were necessarily due to increases in relative load arm lengths alone. By contrast,relative hip and ankle load arm lengths increased significantly only during symmetrical gaits on the ground, although non-significant increasing trends were observed in most other conditions. Increases in relative hip and ankle moments must therefore be primarily due to relative increases in hind limb SRF magnitudes. Additional ontogenetic changes in joint moments included a relative increase in knee joint moments for trailing limbs during locomotion on the ground and a relative decrease in wrist moments for trailing limbs during asymmetrical gaits on the pole.

Despite well-established differences in the mechanics of symmetrical and asymmetrical gaits (Cavagna et al.,1977; Biewener,2003), size-related changes in joint mechanics were relatively consistent across gait types within substrates(Table 3). By contrast, across variation in gait type, size-related changes in joint kinematics were consistently observed, more frequently during pole locomotion than during ground locomotion (Table 3). Greater size-related adjustments during pole locomotion probably reflect the stability constraints of moving on narrow, fixed-diameter supports,particularly as body size increases(Cartmill, 1985; Schmitt, 1994; Schmitt, 1998; Schmitt, 1999; Franz et al., 2005; Wallace and Demes, 2008; Young, 2009).

The data presented here demonstrate that a combination of relatively short limbs, more extended elbow postures and changes in limb force distribution allowed infant squirrel monkeys to frequently reduce elbow, hip and ankle joint loading relative to older and larger individuals, even when traveling at the same absolute speed. Such changes should reduce the muscle force needed to maintain posture or effect movement at a joint, although data on the ontogenetic scaling of muscle moment arms in squirrel monkeys would be required to fully corroborate this prediction (i.e. Carrier, 1983; Young, 2005). In the following sections there is a review of what is currently known about the ontogenetic scaling of limb length, joint posture and SRF in other mammalian taxa in order to place the results presented here in a broader comparative context.

Data on ontogenetic allometry of limb growth in primates and other mammalian taxa were collated from the literature and are presented in Fig. 8. These data suggest that positive allometry of limb growth is typical of most mammals. Among primates,long bones grow with positive allometry in almost every species yet studied. Primates showing consistent isometry or negative allometry are either adapted for vertical clinging and leaping (i.e. sifakas, Propithecus spp.), a very peculiar form of arboreal locomotion that requires young individuals to be relatively long-limbed in order to produce the necessary accelerations(Ravosa et al., 1993; Demes et al., 1999), or are more terrestrial than most other primates (i.e. baboons, Papio cynocephalus). Among other mammalian orders, the only species showing consistent negative allometry of long bone growth are highly precocial and are required to stand and locomote with adults soon after birth (i.e. cuis, Galea musteloides, and domestic pigs, Sus scrofa). In these species, relatively longer limbs allow perinatal animals to keep pace with older and larger conspecifics during travel and effectively evade predation(Howell, 1944; Pennycuick, 1975; Trillmich et al., 2003).

Whereas the proximal and middle limb segments typically grow with positive allometry, autopodia (i.e. hands and feet) scale to body size with negative allometry or isometry in almost all mammals yet studied. Negative allometry of hand and foot size has been interpreted as an adaptation allowing young primates to cling to their mothers and young quadrupeds in general to negotiate `adult-sized' substrates and maintain a larger, more secure base of support (Jungers and Fleagle,1980; Raichlen,2005; Lawler,2006). There may, however, be a cost associated with relatively larger hands and feet – by lengthening the load arm of the SRF,relatively longer autopodia may increase joint moments at the wrist and ankle,particularly if these segments are positioned roughly perpendicular to the SRF vector, as in palmigrade or sub-digitigrade animals. As such, mechanical compensation for ontogenetic limits on force production might be particularly necessary at the wrist and ankle. In fact, strong negative allometry of triceps surae anatomical mechanical advantage at the ankle has previously been documented among black-tailed jackrabbits(Carrier, 1983; Carrier, 1995). Although I was unable to measure hand and foot lengths in the current study, future research should combine allometric analysis of growth in these segments with data on ontogenetic changes in wrist and ankle postures.

Previous studies of postural development in mammals have indicated some variability between taxa, in contrast to the relative uniformity observed in patterns of limb growth across mammals. In most species yet studied, including vervet monkeys, domestic cats, tree shrews and humans, early locomotor efforts are characterized by increased joint flexion(Peters, 1983; Vilensky and Gankiewicz,1989; Howland et al.,1995; Schilling,2005; Hallemans et al.,2006). Increased flexion during early locomotion has also been qualitatively noted in tufted capuchin monkeys, rhesus macaques, Japanese macaques, and chimpanzees (Hildebrand,1967; Kimura,1987; Fragaszy,1990; Nakano,1996). By contrast, young rats have been shown to use more extended limb postures than older individuals and adults(Westerga and Gramsbergen,1990). Finally, cuis and baboons show a combination of patterns:increasing flexion at some joints and extension at others(Schilling, 2005; Zeininger, 2007).

Patterns of ontogenetic change in SRF magnitudes are characterized by similar amounts of variability. Kimura(Kimura, 1987; Kimura, 2000) found that relative forelimb and hind peak vertical SRF magnitudes were significantly greater among infant Japanese macaques and chimpanzees than among juveniles and adults. Studies of human locomotor development have also documented patterns of decreasing relative peak vertical SRF magnitudes with increasing age (Beck et al., 1981; Diop et al., 2005). Similarly,using an approximation of SRF based on average limb contact time, Pennycuick(Pennycuick, 1975) estimated that free-ranging infant gnus (Connochaetes taurinus) experience significantly greater peak forces than adults traveling at identical absolute speeds. By contrast, other studies of a small, but diverse group of animals have documented isometric variation in SRF magnitudes during ontogeny. Main and Biewener (Main and Biewener,2004; Main and Biewener,2007) found that SRF magnitudes scaled in direct proportion to body mass when growing goats (Capra hircus) and emus (Dromaius novaehollandiae) were traveling at constant relative (i.e. dynamically similar) speeds. Hallemans et al.(Hallemans et al., 2006) also found that after controlling for relative speed, size-adjusted peak vertical and propulsive forces did not change with walking experience (i.e. age) among newly walking humans. As indicated by the level of variation between species and between different studies of the same species (i.e. humans), greater research on the patterns and causes of ontogenetic variation in SRF magnitudes is clearly needed. Additionally, in order to address both the ecological and physiological implications of ontogenetic changes in limb loading, future studies should make comparisons at matched absolute speeds as well as matched relative speeds.

In summary, long bone growth is positively allometric in most mammals,indicating that young mammals are relatively short-limbed for their size. However, it should be noted that the available growth data are strongly biased towards primates, indicating the need for additional ontogenetic study of other mammalian taxa. Patterns of postural development are more variable. Most animals thus far studied, are characterized by increased flexion during early locomotor efforts, although others taxa increased extension or a combination of increased flexion and extension, depending on the joint being examined. Patterns of ontogenetic change in SRF magnitudes appear equally variable, with some studies demonstrating relatively greater limb loading in younger animals whereas others document isometry of force magnitudes.

Pervasive positive allometry of proximal and middle limb segment growth may, by itself, constitute an independent solution to ontogenetic limits on locomotor performance. As illustrated in the current study (Figs 5 and 6) and in previous studies comparing the effects of intra- and interspecific variation in limb length on patterns of joint loading (Polk,2002; Gruss,2007), relatively shorter limbs have the potential to shorten joint load arms, thereby mitigating joint moments. The assumption that limb growth determines ontogenetic variation in joint loading motivated Carrier's(Carrier, 1983) and Young's(Young, 2005) morphometric analyses of developmental joint mechanics. However, as reviewed above, the data presented here failed to show a strong direct link between limb length and joint loading in growing squirrel monkeys, despite strong positive allometry of limb growth. Rather, variation in joint and limb posture were consistently the strongest predictors of variation in joint load arm lengths. It is probable that the flexed joint postures that appear to be typical of most animals during early locomotion are sufficient to lengthen joint load arms, although the precise interaction between developmental changes in posture and limb growth should be evaluated on a species by species basis. Provided SRF magnitudes scale isometrically or with negative allometry, as appears to be the case in the few animals studied thus far, joint moments are likely to be greater as well in animals characterized by increased flexion at early ages. As such, the available data generally support Carrier's(Carrier, 1983) and Young's(Young, 2005) assertion that compensatory growth patterns, such as negative ontogenetic allometry of bony muscle levers, would be required for young mammals to overcome ontogenetic limits on performance, particularly if relative whole-body muscle mass, and available muscle force, were also reduced at early ages.

Why, then, do young squirrel monkeys differ from most other animals in using more extended or similar joint postures relative to older and larger individuals? Among the taxa reviewed above, animals showing increased joint extension early in life are also relatively precocial in their locomotor behavior. Cuis must locomote independently soon after birth(Schilling, 2005). Although rats are born in an altricial state, mature gait patterns are evident at only 20 days after birth (Westerga and Gramsbergen, 1990). Ectothermic vertebrates required to move independently at early ages, such as desert iguanas (Dipsosaurus dorsalis), also show greater joint extension as juveniles(Irschick and Jayne, 2000). Path analyses and hierarchical partitioning of the current dataset emphatically demonstrated that walking with erect, strut-like limbs –where the SRF vector is aligned with the limb's axis – was the most effective way of shortening SRF load arms among developing squirrel monkeys. Postural adjustments to SRF load arm lengths have also been cited as the primary means of mitigating limb muscle force requirements between gaits within individuals (Biewener et al.,2004), between differently sized individuals within the same species (Polk, 2002; Gruss, 2007), and between differently sized species (Biewener,1983; Biewener,1989) [but see Day and Jayne(Day and Jayne, 2007)]. In sum, extended joints may constitute an effective behavioral means for relatively precocial infants to shorten joint load arms, reduce joint loading,and limit the muscle force required to maintain joint postures. To address the generality of this hypothesis, future studies should include animals encompassing a diverse array of life history strategies, sample a wider, more continuous range of ages (i.e. not just young infants and adults) and combine the performance-related measures examined here with morphometric and experimental data on the ontogenetic scaling of available muscle force and muscle moments arms. Additionally, in order to fully incorporate ontogenetic studies of joint mechanics into an adaptive evolutionary framework (e.g. Arnold, 1983), future research should combine morphometric and lab-based data on musculoskeletal growth and locomotor performance with field-based assessments of survivorship.

This study sought to answer two primary questions. First, among growing squirrel monkeys, how do limb lengths, joint postures, SRF magnitudes, SRF angles, joint load arm lengths and joint moments vary ontogenetically? Second,what proportion of the ontogenetic variation in load arm length can be attributed to limb growth? Previous morphometric studies of mammalian developmental joint mechanics have assumed that joint postures and SRF scale isometrically during growth and that load arm length is proportional to distal limb length (e.g. Carrier,1983; Carrier,1996; Young,2005). Isometry of joint loading across development would imply that infants, with reduced muscle mass and less available muscle force, would require compensatory musculoskeletal adaptations, such as negative allometry of anatomical mechanical advantage, if they were to match adult-like levels of locomotor performance (Pennycuick,1975; Carrier,1983; Carrier,1996; Trillmich et al.,2003; Main and Biewener,2004; Young,2005; Herrel and Gibb,2006; Lawler,2006). Contrary to predictions, joint postures, SRF vectors, joint load arm lengths and joint moments frequently did not scale isometrically to body size in growing squirrel monkeys. Rather, ontogenetic changes in joint postures and body mass distribution mitigated joint loading at the elbow, hip and ankle during several of the substrate-by-gait conditions examined here. A full assessment of the realized benefits of such ontogenetic changes in joint loading requires further data on the scaling of muscle force and muscle moment arms during growth.

In summary, the size-related kinematic and kinetic changes documented in this study constitute an example of `behavioral' compensation for ontogenetic limits on locomotor performance, a phenomenon previously identified only in ectothermic vertebrates (Jayne and Bennett, 1990; Irschick,2000; Miles, 2004; Irschick et al., 2005; Dial et al., 2008), suggesting that growing mammals may use a combination of anatomical and behavioral means to overcome ontogenetic limits on locomotor performance. The existence of such behavioral mechanisms may challenge the status of allometric musculoskeletal growth trajectories as `necessary' mammalian solutions to ontogenetic limits on locomotor performance (cf. Carrier,1983), at least among species characterized by patterns of locomotor development similar to the squirrel monkeys observed here. Morphological adaptations to ontogenetic limits on locomotion, such as negative ontogenetic allometry of anatomical mechanical advantage, may be more critical among animals characterized by more flexed limb postures early in life or in species where young are under extreme selective pressure to accelerate quickly to evade predation, such as the jackrabbits studied by Carrier (Carrier, 1983).

LIST OF ABBREVIATIONS

 
• SRF

  substrate reaction force

 
• r

  average muscle moment arm length (weighted by individual muscle force)

 
• R

  SRF load arm length

 
• F_S

  SRF magnitude

 
• F_M

  total muscle force

 
• F_x(Vpk)

  magnitude of the fore–aft component of the SRF at peak vertical force

 
• F_y(Vpk)

  magnitude of the vertical component of the SRF at peak vertical force

 
• M_b

  body mass