The large theropod dinosaur Tyrannosaurus rex underwent remarkable changes during its growth from <10 kg hatchlings to >6000 kg adults in <20 years. These changes raise fascinating questions about the morphological transformations involved, peak growth rates, and scaling of limb muscle sizes as well as the body's centre of mass that could have influenced ontogenetic changes of locomotion in T. rex. Here we address these questions using three-dimensionally scanned computer models of four large, well-preserved fossil specimens as well as a putative juvenile individual. Furthermore we quantify the variations of estimated body mass, centre of mass and segment dimensions, to characterize inaccuracies in our reconstructions. These inaccuracies include not only subjectivity but also incomplete preservation and inconsistent articulations of museum skeletons. Although those problems cause ambiguity, we conclude that adult T. rex had body masses around 6000–8000 kg, with the largest known specimen (“Sue”) perhaps ∼9500 kg. Our results show that during T. rex ontogeny, the torso became longer and heavier whereas the limbs became proportionately shorter and lighter. Our estimates of peak growth rates are about twice as rapid as previous ones but generally support previous methods, despite biases caused by the usage of scale models and equations that underestimate body masses. We tentatively infer that the hindlimb extensor muscles masses, including the large tail muscle M. caudofemoralis longus, may have decreased in their relative size as the centre of mass shifted craniodorsally during T. rex ontogeny. Such ontogenetic changes would have worsened any relative or absolute decline of maximal locomotor performance. Regardless, T. rex probably had hip and thigh muscles relatively larger than any extant animal's. Overall, the limb “antigravity” muscles may have been as large as or even larger than those of ratite birds, which themselves have the most muscular limbs of any living animal.

Also pertaining to the third aim of our study (reconstructing locomotor ontogeny), we focus on a second hypothesis. We propose that the aforementioned changes of body proportions in tyrannosaur ontogeny would have shifted the centre of mass (COM) cranially (and dorsally). Such a shift would have further worsened locomotor performance because of increasing demands placed on limb muscles by lowered effective mechanical advantage in the hindlimbs [20] , [21] , [28] , [47] – [49] . Few studies have examined ontogenetic changes of COM (or inertia) in any taxa, extant or extinct, so this is a novel feature of our study. Initial examinations of extant crocodiles (Crocodylus johnstoni) and junglefowl Gallus gallus (ancestral wild chickens) indicate that they may shift their COMs craniodorsally during ontogeny [31] , as we propose for tyrannosaurs, which may thus be a common (or even homologous) archosaurian pattern. Alternatively, the COM might maintain a conservative position or even shift caudally, potentially forestalling decreases in running performance. By combining our estimates of hindlimb extensor muscle masses with COM position we will re-evaluate previous estimates of locomotor ontogeny (and locomotor abilities at adult sizes) in tyrannosaurs while still considering the ambiguities inherent to these estimates.

These methods are then used to test our hypothesis that juvenile tyrannosaurids, using the “Jane” specimen as an exemplar, had relatively larger extensor (e.g., antigravity) muscles than adults. This hypothesis agrees with observations of the proportionately smaller heads, longer limbs, more lightly-built presacral regions and possibly shorter tails in juvenile tyrannosaurs [5] , [13] , [17] , [66] . It also fits well with speculations about their increased relative locomotor performance and differing predatory ecologies [21] , [39] , [56] as well as general trends in the ontogenetic scaling of extant taxa (references cited above). Contrary to this hypothesis, Persons and Currie [19] concluded the opposite; that leg muscles scaled with positive ontogenetic allometry in tyrannosaurs to mitigate the relative decline in locomotor muscle force and power output as body size increases interspecifically (e.g., [47] – [55] , [72] ), although other species may do the opposite during ontogeny (e.g., [49] – [55] ). Alternatively, it is also important to consider that perhaps unknown parameters are too numerous to falsify either hypothesis or the null hypothesis of no difference in leg muscle mass across ontogeny in extinct dinosaurs.

We develop two methods here for analyzing limb muscle masses. We use a methodology (detailed below) similar to that of Persons and Currie [19] to investigate tyrannosaur ontogenetic changes in the mass of M. caudofemoralis longus, certainly the largest extensor muscle in the hindlimb. As a simple but critical validation step, we apply the same method blindly to a crocodile skeleton and compare the muscle mass estimate to that obtained via dissection of the same specimen. A second approach that we implement here is to subtract actual bone volumes from estimated segment volumes (as suggested by [30] ) to estimate soft tissue volume, then use those volumes to estimate how large the extensor muscles acting about each hindlimb joint may have been.

Our study's third aim is to re-examine locomotor ontogeny in Tyrannosaurus. We do this first by estimating hindlimb muscle masses. The sizes of those muscles are pivotal information for reconstructing how large and small tyrannosaurs moved [20] , [21] . The hindlimb extensor (antigravity) muscles would have been critical determinants of maximal locomotor performance in tyrannosaurs [20] – [22] , [24] , [25] , [67] , [68] . This notion fits prevailing evidence for speed limitations in other bipeds (reviewed in [21] , [67] ; also [69] – [71] ). Hutchinson et al. [28] provided simple geometric estimates of how much muscle mass different limb segments could contain (ranging from 2–10% body mass from the shank to the thigh, depending on assumptions about the bulk of soft tissues). An improved reconstruction methodology for estimating the mass of M. caudofemoralis, an important hip extensor muscle group, was presented by Persons and Currie [19] .

However, a reasonably complete juvenile specimen probably pertaining to Tyrannosaurus rex has recently been discovered (BMR P2002.4.1; [58] – [59] ). This specimen (still awaiting formal description) was originally referred to the small Maastrichtian tyrannosaurid Nanotyrannus [60] , which many workers consider to represent juvenile individuals of T. rex [16] , [61] – [62] . Discovery of this specimen, combined with past studies of tyrannosaur ontogeny and individual variation [16] , [61] – [66] , opens new potential to quantify ontogenetic changes of 3D body proportions in T. rex using modern computational tools, and to compare these results to data from previous studies of T. rex and other taxa (extant and extinct). This would illuminate more precise details of how the biology of this famous carnivore may have changed during its lifetime or otherwise varied between individuals, related to the second aim of our study noted above.

The remarkable limb and body proportions of tyrannosaurs changed appreciably as the animals grew from relatively more slender and long-legged (i.e., cursorial; [43] , [44] ) juveniles into more robust adults [5] , [17] . A similar pattern seems to hold for other lineages of giant theropods (e.g., Allosauroidea: [45] ). Tyrannosaurus rex is unusual relative to other dinosaurs in that its femur scales according to a static stress model (i.e., length scales half as fast as diameter) through ontogeny, the only one of 24 examined non-avian dinosaurs species to unambiguously exhibit that growth pattern [37] . Tyrannosaur growth itself is fascinating from a biomechanical perspective, because T. rex took less than 20 years to develop from <10kg hatchlings into >6000 kg giants [39] , [46] . The fast growth of larger tyrannosaurs should have caused them to face rapidly shifting size-related biomechanical constraints involving ontogenetic allometry [31] , [47] – [55] . Hutchinson and colleagues [20] , [21] , [25] used an integration of limb musculoskeletal biomechanics and estimates of inertial properties to assess Currie's [17] , [56] hypothesis that young tyrannosaurs were at least relatively, if not absolutely, more athletic than larger, older tyrannosaurs and thus that tyrannosaur ecology changed drastically across ontogeny (as in some other species; e.g., Alligator [57] ). Juvenile Tyrannosaurus specimens were unknown at the time so they used a juvenile Gorgosaurus libratus specimen for comparison, and found support for Currie's [17] , [56] inference.

Body mass is a key parameter for understanding the biology of any animal, let alone taxa that represent physical extremes such as Tyrannosaurus rex and other large dinosaurs. Most studies that use body mass estimates for dinosaurs rely on equations that are extrapolated from multi-taxon datasets of presumed adult individuals [32] , [36] . These approaches approximate an average model value based on the sampled taxa and thus minimize taxonomic differences in body shape and volume. They also do not address individual and ontogenetic variation, which is critical because intraspecific scaling patterns do not follow interspecific scaling in non-avian dinosaurs [37] . Studies of dinosaurs that incorporate ontogenetic series [38] – [41] often utilize Developmental Mass Extrapolation (DME) [42] to estimate mass for immature specimens using the cube of femoral length and an assumption of isometric scaling. While this method derives support from observations on scaling in extant archosaurs, it remains untested for fossil taxa. The application of advanced 3D scanning equipment in palaeontology allows estimation of body mass for individual specimens, and thereby a potentially more accurate understanding of both ontogenetic changes in body mass as well as individual variation in this critical parameter. Such insights should increase the precision in analyses of dinosaurian biology that employ body mass estimates and allow for independent evaluation of the DME approach. Our study's second aim is to apply these techniques to Tyrannosaurus to revisit estimates of its growth rate and the reliability of DME.

The use of computer models to make these estimates of inertial properties avoids some major problems of alternative approaches, such as use of scaling equations (e.g., [32] – [33] ). Such models have the advantage of producing inertial property estimates specific to the taxon (and specimen) of interest and its entire anatomy. The accuracy of these estimates is presumably no worse than that of alternative methods and quite possibly better, but still has a substantial margin of error of 50% or more (e.g., [31] ), which always demands careful consideration. Here our primary aim is to use three-dimensional (3D) scanning of a putative juvenile T. rex specimen (see below and [34] ) as well as four of the best-preserved and most well known adult specimens ( Table 1 ); including the famously well-preserved skeleton “Sue” (FMNH PR 2081); to estimate individual and ontogenetic variation in limb and body dimensions in T. rex ( Figure 1 ) and their repercussions for tyrannosaur biology. As part of this aim, we seek to further understand the degree of errors intrinsic to such 3D reconstructions; as well as alternative methods such as scaling equations, following up on previous studies [26] – [31] , [35] .

Most recently, studies of tyrannosaur hindlimbs have emphasized biomechanical procedures to quantitatively test functional hypotheses about locomotion [20] – [28] . Furthermore as part of a renaissance in musculoskeletal biomechanics that was initiated by Alexander [9] – [10] and hastened by technological developments, tyrannosaur morphology has been a major subject for computational analyses of three-dimensional (3D) whole-body shape. In particular, studies have focused on inertial properties such as mass, centre of mass and inertia [26] – [31] that are critical for biomechanical analyses.

Tyrannosaurus rex Osborn 1906 (Dinosauria: Theropoda: Tyrannosauridae) [1] has become an exemplar taxon for palaeobiological studies because of its large body size, suite of near-complete skeletons, and popular appeal [2] . It is immediately recognizable to many people worldwide for not only its giant size but also its unusual anatomy: large head, small forelimbs, robust torso and tail, and long hindlimbs. The anatomy of its hindlimbs in particular has been the focus of comparative studies for almost a century [3] – [19] .

Filled circles represent our mass estimates derived from digital modelling; squares represent data points generated using Developmental Mass Extrapolation [42] . The red line represents the best fit curve for modelled data; the black one is for the DME estimate. Stippled lines represent the 95% confidence intervals for the model-based estimates. The MOR specimen is treated as being 16 years old in this plot (see text).

Given the overall larger specimen masses derived from our modelling of scans of mounted specimens, it comes as no surprise that the adjusted growth curve for Tyrannosaurus rex finds significantly higher growth rates during the exponential phase ( Fig. 6 ), than previously reported [39] . Using the model based mass estimates, maximum growth rates for T. rex during the exponential stage are 1790 kg/year-- a more than twofold increase from the reported value [39] . Comparisons between model-based and DME mass estimates are relatively favourable. The DME-derived estimates fall within the 95% confidence limit when the MOR specimen is assumed to be 14 years of age, though not if it is considered to be 16 years old ( Fig. 6 ). Under the latter scenario, its mass and age are almost identical to that of Stan, essentially removing any scatter in a plot with so few data points. Both the AICc and F-tests overwhelmingly favour a single model (i.e., one growth curve) to fit both data series.

By subtracting bone volumes from thigh and shank segment volumes, then adjusting these volumes to estimate extensor muscle masses, we aimed to derive more accurate muscle mass estimates for Tyrannosaurus. Generally, the lower end range of our estimates ( Table 9 ) compares reasonably well with others' [28] (total extensor masses per limb of 12–14% body mass). Total extensor muscle masses estimated for each of the three major joints of one pelvic limb varied in similar ways for each of our models ( Table 9 ). The hip joint extensor muscle masses, including our CFL estimates, ranged from 8.8–15.1% body mass for the adult tyrannosaurs and 11.6–15.9% for the Jane specimen. The knee extensor masses ranged from 3.3–7.1% body mass in adults vs. 4.9–7.3% for Jane (squarely overlapping the results for other models made by KB). Finally, the ankle extensor masses ranged from 1.9–3.2% body mass in adults vs. 2.5–3.9% for Jane. Thus the total extensor muscle masses for one hindlimb were about 14–25% body mass for adults vs. 19–27% body mass for the Jane specimen. These estimates naturally reflected similarities and differences noted above for CFL mass and segment mass ranges for the adults vs. the juvenile. They are heavily dominated by the hip/thigh musculature including the CFL (proportion of total limb extensor mass roughly 60% hip vs. 40% knee and ankle). Our estimates do not consider the small (∼6%) difference between the density of our models' homogeneous flesh (1000 kg m −3 not including air cavities, which are negligible for the limbs and tail) and the density of vertebrate muscle.

Our validation test for the extant Crocodylus johnstoni specimen indicated reasonable accuracy for our methodology, concurring with a similar reconstruction technique [19] . The estimated muscle mass was 0.311 kg and the measured mass was 0.316 kg, although the assumed density (1000 kgm −3 ) was lower than the likely actual density (∼1060 kgm −3 ); a more comparable estimate would be 0.330 kg (4.4% overestimate of volume or mass vs. actual). This corresponds to a CFL mass of 1.6% body mass (per hindlimb) for the adult crocodile; identical to values in adult Alligator [54] .

Dorsal and right lateral views are shown on top, and in the bottom row are caudal views of the right femur and then caudal vertebrae (8 th and 17 th ). Red shaded volumes are the M. caudofemoralis longus reconstruction. Note a small space for M. caudofemoralis brevis (not reconstructed) is left around the ilium/sacrum and lateral to the CFL insertion.

Our models of four adult Tyrannosaurus had remarkably consistent CFL muscle mass estimates ( Figure 5 ; Table 8 ) of 162–189 kg for minimal models (the smaller value surely caused by the smaller tail segment of MOR). The relative masses ranged from 1.95–3.9% body mass for all the different size (Min and Max) models. This contrasts with the potentially slightly larger CFL muscle mass in Jane (at 3.4–3.8% body mass).

Absolute COM positions are shown in Table 6 , but in Table 7 we focus on relative COM positions, using a variety of normalizing metrics (following [31] ). Here we focus on the results for femur length as a normalizing factor because a femur is well preserved for all specimens whereas the other metrics are more likely to have been influenced by skeletal preservation and mounting biases, although they are all presented in Table 7 . Moderate variation in COM position was found along the y (dorsoventral) axis between the different iterations of each specimen, but the dorsoventral COM position varied more widely between specimens. Expressed as % femur length ventral to the hip, the COM was estimated to be 22–35% for the Carnegie and MOR specimens, a surprisingly conservative but otherwise similar 26–28% for Sue, 38–50% for Stan, and 35–54% for Jane. As expected from previous studies [28] , [31] , we found greater variation in craniocaudal (x-axis) COM position for each of our specimens. Again using femur length to normalize the results, the COM distance from the hip joint ranged from −1.9% (i.e., behind the hip; this highly implausible result was only found for the Carnegie specimen) to 88.8% (in Sue specimen) femur length in “Most Cranial” and “Most Caudal” extreme models. These were intentionally designed (see Methods ) to maximize the possible range of COM positions given a set value for morphological variability (here, ±20% radial dimensions for body segments), rather than to satisfy the valid but hard to quantify criteria of appearing ‘natural’ (subjectively plausible to an observing anatomist). Hence they should not be assumed to be equally as plausible as less circumspect, ‘natural looking’ models. Instead, the range of ∼30–45% femur length for most minimal and maximal (also most dorsal and ventral) models; albeit with Sue at 44–61%; would cover our qualitative assessment of the most plausible COM range for adult Tyrannosaurus. Jane's COM estimates (30–35% femur length) fall at the lower boundary of this range, but as these ranges overlap appreciably it is not conclusive whether the COM experienced a craniad ontogenetic shift (0–15% femur length).

Estimated COM position varied relatively more among the four adult Tyrannosaurus models than among different iterations of individual specimens. However, the similar-sized Carnegie, Stan and MOR specimens, despite apparent investigator bias resulting in different body mass estimates, had roughly comparable cranial positions of the COM (∼0.55 m from hip) for the minimal models. The pronounced difference of these models' COM cranial positions from the Sue model (0.80 m) can be partly explained by the greater size of the FMNH skeleton -- with contributions from errors such as the artefactually wide, tall ribcage (see above and Discussion ). Likewise, the cause of the ∼50% more ventral COM position in the Stan model (0.495 m vs. ∼0.34 m below the hips in the other three large models) is unclear but likely a result of modelling errors.

Table 6 shows our results for our five whole-body models. For the four adult Tyrannosaurus models we obtained a wide range of body masses and COM positions because our methodology purposefully produced wide “error bars.” We expanded our range of models to extremes that we qualitatively consider to be too skinny, too fat, or too disproportionate. We did this with the expectation that the actual (but effectively unknowable) body mass and COMs would fall within this range. Body mass (for minimal and maximal models; here rounded to nearest 100 kg) ranged from ∼5800–10800 kg for the comparably-sized Stan and MOR models made by KB, whereas the Carnegie and Sue models made by VA+JM produced proportionately larger masses: ∼7400–14600 kg for the former and ∼9500–18500 for the latter. For the Jane specimen we obtained masses of 639–1269 kg. Despite the occasional differences in individual body segment dimensions noted above, our maximal models were consistently around 1.9 times as massive as our minimal models.

There were several obvious patterns of differences between our five models, which is most clearly summarized by considering the total masses dedicated to limb vs. body segments in the models by VA+JM and KB. Investigators VA+JM tended to produce models having 17.3–20.4% body mass in all four limbs, whereas KB's two adult Tyrannosaurus models did not overlap with this range (in minimal or maximal models) at 25.6–31.3% body mass, and his Jane model was larger (26.8–36.2% body mass). Conversely, the main body segments were relatively larger when constructed by VA+JM vs. KB. We consider these and other errors or biases in the Discussion .

Pelvic limb masses, however, tended to be proportionately quite similar among the four adult Tyrannosaurus models. The thigh segment was the main exception – it ranged from 12.6%–17.4% body mass (with larger relative values for KB's models). Jane (with the relatively shorter femur; Table 2 ) had a relatively smaller thigh mass than KB's other models at 15.5%. These differences, and wide variation in scaling the thigh segment to maximal dimensions, are likely mainly caused by subjective investigator biases but also by large differences in how the torso was assembled in the various mounts, which affects overall mass and its proportional distribution. The shank mass estimates fell within a relatively narrow range of 4.3–5.9% body mass. Jane (with the relatively longest “tibiotarsus”; Table 2 ) was at the upper end of this range, again as expected. Finally, the foot (“tarsometatarsus” plus digits) segment masses varied from 1.3%–2.9% body mass, with the Sue specimen standing out as having a relatively small total segment. This was probably caused by differing subjective reconstructions including the inflated size of the torso, but also slightly by the absence of digit 1 and metatarsal 5 in the Sue -- and Carnegie-- specimens.

Our analysis of main (axial) body segments raised the question: Do these apparent errors and biases persist in the limb segments? Pectoral limb masses of the adult specimens ranged from 0.1–0.4% body mass which we anticipated due to great subjectivity in assigning masses to the proximal forelimb vs. torso. Jane's forelimb mass was strikingly larger at 2.8% body mass, which did not seem to be an artefact of our methodology and fits well with ontogenetic changes discussed in the Introduction .

The entire tail segment volume also showed wide variation corresponding to that observed for skeletal tail lengths ( Table 3 ); 10–21% body mass, with Jane at 11.5% body mass. The models by VA+JM had similar relative masses (11.6–12.3%) and were scaled similarly for the maximal models (2.8–3 times more massive), whereas the models by KB had relatively smaller (9.9%; Stan) and larger (20.7%; MOR) masses and were scaled to maximal sizes by 2–2.2 times (2.3 times for Jane). Again the MOR estimates seem to have been skewed by different skeleton dimensions in particular segments, for the same reasons noted above.

The body (i.e., torso; base of neck to caudal end of sacrum) segment masses were remarkably similar for most specimens in terms of relative sizes, ranging from 54.2–58.5% body mass for the Carnegie, Sue and Stan specimens. Jane's body segment mass estimate was below this range (51.2%). However, clearly relating to its implausibly short torso (see above), the MOR specimen had a very small relative torso mass of 35.3% body mass. Almost 10% of the body mass that should have been apportioned to this segment instead ended up in the head and neck, plus more in the tail (below). We discuss this more later. Finally, even though they agreed on and used the same methods VA+JM and KB scaled their body segments differently: the maximal models were ∼1.9 times larger for the former but only ∼1.7 times larger for the latter (see Discussion ).

Neck segment masses also varied widely, from 3.2–7.3% of body mass; again with Jane lying in the middle of this range (4.6%). The largest values are for the Sue and MOR specimens. The much larger value for the MOR neck segment compared to Stan and Carnegie again relates to their drastically different skeletal proportions, which determine the main body:neck segment lengths in our models. This difference (i.e., a relatively long neck, but shorter main body segment in MOR) is reflected in the relative masses of the main body and neck segments for MOR and Stan in the study of Bates et al. [30] . The caudal limit of the main body segment was defined by the caudal tip of the ischium, and the cranial limit by the cranial margin of the pectoral girdle (i.e., coracoids). The boundaries of the neck segment were defined by the latter boundary and the back of the skull. Thus any caudal shift of the pectoral girdle's boundary would increase the neck segment's volume. Therefore differences in mounted orientations of skeletal elements as well as missing/reconstructed elements, on top of real individual variation, contribute to these wide discrepancies.

The heads of all four adult tyrannosaurs constituted 4.1–11.1% of the total (minimal) body mass, with the Jane specimen falling into the middle of that range (5.8%). The smaller values are for models made by VA+JM whereas the larger values come from KB's models and thus there may be some investigator bias in these models (see Discussion ). Furthermore the MOR model by KB has an anomalously large head which is an artefact of the lower scan quality. The lower resolution of the raw scans results in a poor quality triangular mesh in the skeletal model of the MOR specimen [30] , giving the digital skull bones an irregular topology. This error artificially inflates the mediolateral dimensions of the skull, creating both a larger external outline and smaller internal chambers, hence lowering the size zero-density respiratory structures. This has led to a substantially greater mass estimate—e.g., the density of the MOR reconstruction is 985 kg m −3 vs. 870 kg m −3 for Stan and 850 kg m −3 for Jane (cf. 791 and 807 kg m −3 for the Sue and Carnegie specimens by VA+JM). The head boundaries for Sue were derived from a CT scan of the mounted reconstructed skull, which attempts to correct for crushing and warping in the original. Although we acknowledge the potential for such reconstruction to affect overall volume, the ratio of head to body mass in Sue is roughly comparable to those of our Carnegie and Jane (and to a lesser degree, Stan) models.

Results for our five models are shown, segment by segment, in Table 4 for body (head to tail) and Table 5 for limb (pectoral limb and individual pelvic limb segments). Here we focus only on the “minimal” models, for brevity, but maximal model results are also included in Tables 4 and 5 . Figures 2 , 3 , 4 show the resulting models.

These similarities and differences are reflected in the skeletal lengths and proportions: all four large tyrannosaurs had femur lengths around 1.3 m, body lengths 11–12 m and limb lengths 3–3.3 m, with Jane's lengths about 50–65% of these. The MOR specimen has an anomalously short body or torso (GAD), about 70% the length of the other three adults, whereas the fragmentary Carnegie specimen's extensive reconstruction is evident in its short tail length (<6 m and <50% of body length; all other specimens are closer to 7 m and >50% of body length). Again, the Sue specimen tends to have the greatest lengths of the four adult specimens and the other three are markedly similar in most dimensions. However Sue's tail is reconstructed beyond the 27 th caudal and could be artificially foreshortened as a result (cf. its tail length vs. that of Stan). Yet an alternative explanation, based on the number of actually preserved caudals (none in the CM 9380, 15 in MOR 555, 31 in BHI 3033, 36 in FMNH PR 2081), is that the seemingly foreshortened tail in the Sue specimen is actually more representative of the actual tail length, whereas other reconstructions have overestimated tail lengths. More complete Tyrannosaurus discoveries would resolve this issue more conclusively (see also discussion in [16] ).

The estimated volumes of each major pelvic limb bone (or set thereof) are shown in Table 2 with additional data on femur/tibiotarsus volume (as a gauge of proximal vs. distal limb size as well as preservation quality or other biases). Table 3 shows skeleton-only measurements of femur length, total head-tail tip body length, gleno-acetabular distance (GAD; between the approximate centres of those joints)/body length, limb length, and tail length (from the first free caudal vertebra)/body length. The bone volumes of our four adult tyrannosaur specimens are roughly similar. The largest specimen (Sue) has larger absolute bone volumes than all but the Carnegie specimen's, whose dimensions are partly inflated (especially for the lower limb) by the presence of a supportive metal mounting framework ( Figs. 2 , 3 , 4 ) that was too tightly integrated with the bones to easily be separated. However, because the measurements for Sue were derived from CT rather than coarser-resolution optical scans, they also represent the most accurate volume measurements in the dataset. The MOR specimen exhibits a disproportionately large tibiotarsus and small femur. The latter could be partly attributable to crushing ( [11] ; although Hutchinson et al. [22] corrected for much of this) whereas the cause of the former difference is unclear but could easily be due to scanning methodologies and preservational distortion. Hence all values must be viewed with some caution and not read literally. The “Jane” specimen's bone volumes are around one-tenth to one-fifth those of the adult specimens. Figures 2 , 3 , 4 show the 3D scans of the skeletons used.

We begin by outlining our results for models of the skeletons alone, which we expect to have the highest precision (adding flesh to the skeletons inevitably increases subjective errors). We then compare the masses of individual fleshed-out body segments among our five main specimens. Third, we examine mass and COM differences for the whole body models (and sensitivity analysis thereof) for all five specimens. Fourth, we detail our results for key extensor muscle mass estimations for each specimen. Finally, we present revised growth rate estimates for Tyrannosaurus and compare these with DME methods.

Discussion

What were the body segment dimensions of an adult Tyrannosaurus rex and how might they have changed during ontogeny? This was the primary question of our study. Although our results are complicated by specimen variation, incomplete preservation, mounting errors and investigator biases (discussed below), there are some general patterns that we infer to be indicative of potential ontogenetic changes in the morphology of T. rex. It is not clear how much (if at all) head, neck or tail mass changed, but the forelimbs clearly became proportionately lighter from juveniles to adults (by as much as a factor of 10), the torso (despite uncertainties about gastral basket shape) became longer and thus heavier (by perhaps 5% of body mass) and the limbs became relatively lighter overall (but presumably heavier for the longer thigh/femur). These trends were suggested by skeletal data (including our Table 3) but reinforced by our “fleshed-out” models (Tables 4,5), across an order of magnitude change in body mass (Table 6). Resulting whole-body density values for all five minimal models overlap those of previous studies [28], [30]: ∼840–890 kg m−3; vs. less plausible values of ∼970 kg m−3 for the maximal models. These changes of proportional masses may have altered the body's COM position. Discounting the anomalous results for the Stan specimen (noted above), our results indicate that a dorsal shift of the COM position in tyrannosaur ontogeny (by <33% femur length; Table 7) is the most plausible conclusion. The COM also may have shifted slightly cranially; possibly by as much as 15% femur length; although the latter shift is more ambiguous and deserves closer examination in the future. However it would be interesting whatever result was supported-- the two extant archosaurs examined to date seem to show a craniodorsal shift of COM position during ontogeny [31] so T. rex might be unusual if it did not experience such a cranial shift. Contrarily, if there was a cranial shift of COM position it could have altered locomotor abilities, as discussed below.

Sources of error in body segment reconstructions of tyrannosaurs The error in any reconstructions of extinct organisms must be addressed as part of interpretations of their biology that rely on those reconstructions, so we consider this problem in relation to the primary aim of our study (above). We discuss two major sources of error here; neither was discussed in any detail by our preceding studies [28], [30]–[31], [35]. First, while our approach offers more precision in terms of gauging individual and ontogenetic variation than many previous methods of mass estimation (see below), it certainly is sensitive to how specimens were restored and reconstructed. Second, investigator biases in digitally reconstructing our specimens to create volumetric models are also evident. Other sources of error such as whole body density (similar for all specimens here although varying for some segments such as the head; see above) and body segment orientations (generally standardized here) were previously addressed [31], so we avoid repetition here. Although relatively little difference exists between the four adult specimens in most linear dimensions, they do differ considerably in the shape and volume of their torso, due to conflated biases of incomplete preservation (and thus varying reconstructions), distortion (of preserved or reconstructed elements) and different mounting techniques (articulations of preserved or reconstructed skeletal elements). This is particularly germane to Sue, and to a lesser degree the Carnegie specimen. Both are mounted with a proportionately longer, wider, and taller torso than the others (Fig. 7). This barrel-chested reconstruction in Sue is an artefact of the dorsal displacement of the transverse processes on the trunk vertebrae, which forced a dorsal displacement of the tubercular articulations and a lateral expansion of the rib cage as a whole. Better understanding of the anatomy of the partial gastral basket (a problem for all our specimens) and potential span of the probable, but pathological, furcula [73] may help constrain the span of the front end of the rib cage and correct the effects of the inflated chest dimensions. Because the torso is by far the largest body segment, overall body mass and the mass distribution between segments should be strongly affected by changes in its dimensions. However, one study of Allosaurus [35] found that this effect was smaller than one might expect, largely because the zero-density lungs and air sacs were expanded with the ribcage in their models. Indeed, the error caused by the large ribcage in our study may have more to do with our coarse distribution of octagonal hoops used to loft flesh around skeletal landmarks (see Methods and Figure 1), but a much broader sensitivity analysis and validation study would be needed to test this. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 7. Comparison of torso/body dimensions for the four large Tyrannosaurus rex specimens. Linear measurements from our digital models show how the gleno-acetabular distance (GAD) is anomalously short in the MOR specimen and the chest is anomalously wider in the Sue and Carnegie specimens. Otherwise all four specimens compare fairly well, considering that the Sue specimen is known to be somewhat larger. https://doi.org/10.1371/journal.pone.0026037.g007 We also commented above on other problems with our various specimens, such as the short torso and large head, neck, and tail in the MOR specimen (e.g., >100% difference of relative masses of head and neck segments in MOR vs. Stan specimens). These problems of incomplete preservation and unstandardized mounting techniques plague all anatomical reconstructions and thus body mass estimates of tyrannosaurs. They thereby limit our understanding of intraspecific variation, such as the likely range of body segment dimensions in adult Tyrannosaurus or ontogenetic changes of morphology [5], [16], [66]. The best-preserved specimen for our study's purposes appears to be the Stan specimen but even it displays some anomalous results; e.g., its ventral COM position. Subjectivity is unavoidable in our reconstructions, so different investigators might produce quite different models. Investigators VA+JM and KB each reconstructed two adult tyrannosaur specimens. The former investigators tended to produce more “fleshy” reconstructions (except for the head segment), partly caused by specimens with more barrel-chested torsos (see [26], [30], [35]), but also influenced by subjective judgements about segment dimensions that are difficult to estimate from skeletal dimensions. In some cases (especially the tail and thigh segments; less problematically the “foot” segment), the difference between adult tyrannosaur models made by different investigators was considerable and could not be ascribed to the above preservational/mounting biases. Furthermore, although our two investigator teams agreed on similar methods in advance there were clearly reproduced differences in how they executed them. For example, KB tended to use more “hoops” than VA+JM, producing more skeleton-hugging body outlines, and had a pronounced narrowed “waist” between the pelvis/sacrum/limbs and the chest (Figs. 2, 3, 4), which is likely a partial explanation for KB's tendency to produce smaller main body segments and differently proportioned leg segments (including min/max scaling differences; Tables 4,5) from VA+JM. Importantly, our minimal models are much more plausible reconstructions than our maximal models, especially for the Carnegie and Sue specimens, considering the realistic density (see above; similar to extant archosaurs'; [31]), less ridiculously rotund appearance (Figs. 2, 3, 4) and smaller masses (closer to previous studies' estimates) of the minimal models (Table 6). We certainly do not contend that our maximal models are results that should be assumed to be equally as valid as our minimal models. Rather, they form extreme upper end values on a spectrum of “fleshiness.” The actual mass of individual Tyrannosaurus specimens should have fallen well below these extremes and much closer to our minimal models', but how much below may never be known. Likewise, as noted above our maximal cranial, caudal, dorsal and ventral models (Table 7) are not as plausible as our minimal models, but rather form extreme boundaries that fairly fully explore the impact of reconstructions on mass, mass distribution and COM position. Regardless of errors in the reconstruction, even the large mass estimate (9502 kg) of the minimal model of the Sue specimen remains a plausible value. Yet considering those errors, the actual mass might have been smaller, whereas the minimal reconstructions by KB could be viewed as lower-end estimates. At present there is little that can be done about such subjectivity except to be aware of it and explicitly acknowledge it, but there is no reason that improvements could not be made in the future. One potential corrective method is that there could be consistent relationships between bone and “fleshed-out” outlines of skeletons in extant taxa that could be quantified and applied to extinct taxa. Allen et al. [31] proposed one fairly simple approach for saurian tails that led to more realistic reconstructions, even though our two teams of investigators still produced somewhat different tail dimensions due to some lingering subjectivity in that method.

What body mass estimate methods are most reliable? Our methods again raise this question that has been recurring in recent studies. Even with consideration afforded to artefacts of reconstruction and investigator biases in fleshing out skeletons (above), it is clear that even minimal body mass values for Tyrannosaurus rex estimated directly from skeletons tend to be larger than published values extrapolated from scale models or limb bone dimensions [9]–[12], [32], [36], [74]. This suggests that researchers have tended to favour ‘skinnier’ reconstructions in the production of scale models, with scaling perhaps magnifying the effects of such subjective choices. Our results (focusing on our minimal models which are more plausible as noted above and below) overlap better with those that employ graphic integrations [26], [75], though they are overall still greater in mass. They compare well with those of other computer volumetric studies using similar methods [28], [30]. As others have noted (see discussion in Hutchinson et al. [28]:p.674), 3D digital reconstruction-based estimates give more plausible results than scaling equations because they are specimen-specific and take into account the unique whole-body morphology and density (as well as explicit sensitivity analysis thereof) rather than assume a very tight relationship between one skeletal dimension and body mass for many species. They thus sidestep the “everyanimal” problem adduced by Pagel [76]. Furthermore the mass estimate equations such as Anderson et al. [32] have other errors and assumptions that leave their merits highly questionable. For example, a classic study [32] adjusted its mass estimates for bipeds using scale models of bipedal dinosaurs whose methodological accuracy is uncertain at best; dubious at worst. The statistical inaccuracies of the same study have also been a subject of recent debate [33], [77]–[78]. Of course, there are no error- or assumption-free methods for estimating dinosaur body dimensions. Although in some cases our results are biased by excessively bulky reconstructions, such as in the case of the Sue specimen, our mass estimates for T.rex derived from actual skeletal data are markedly heavier than proposed by several commonly used equations. These equations have been employed in understanding growth and other life history parameters for this taxon, so next we consider how their usage has influenced estimates of tyrannosaur biology.

Implications for growth rates The second aim of our study was to revisit published Tyrannosaurus rex growth rates with consideration of the validity of the DME approach as well as the accuracy of body masses used in such approaches. Erickson et al. [39] noted that in spite of its impressive magnitude, their estimate of T. rex maximal growth rate was only between a third and half of the expected value for a non-avian dinosaur of that size, when compared to regressions relating peak growth rate to body mass. Our new value for peak growth rate (Fig. 6) would largely erase the reported difference between observed T. rex peak growth rates and the expected theoretical value. Because the latter relationship was calculated using mass extrapolations as a function of limb bone diameters [32] rather than modelling, however, we cannot be sure how much confirmation of our mass estimation this seemingly better-fitting result offers. Considering that the upper end of the published theoretical regression [38] was defined by quadrupedal sauropodomorphs, whose masses are likely underestimated by the Anderson et al. equation ([32];see above), our greater estimate is still in keeping with what is to be expected of growth patterns in non-avian dinosaurs. One effect of finding a faster peak growth rate is that the age at which 50% of adult body mass is reached is shifted earlier into ontogeny. Compared with Erickson et al. [39] who found that Tyrannosaurus rex reached this rate around age 16, our results indicate that this threshold is passed between the ages of 12.9 and 14.2 years depending on how the age of the MOR specimen is set. The lower estimates seems unrealistic given the calculated mass of the Jane specimen, which is a 12 year old individual, but the upper bound appears reasonable if an animal of Jane's size is given two years to grow at or near the peak rate. Our results support the inference of an earlier attainment of somatic maturity (growth asymptotes) at around 16–17 years of age compared to the 18–19 years previously estimated from Sue's histology [39]. Given the paucity of data points, it is prudent to interpret the extrapolated growth curve in light of actual specimen data, which would slightly decrease peak growth rates and increase the age at which 50% of mass is reached. We did not implement such ‘biological’ corrections here, as our focus was on estimating the raw effects of change due to new mass estimates and also to compare our values with those estimated using DME. Likewise, we have not investigated the impact of our results on other inferences made based on these mass estimates and growth rates, such as the correlation of ageing rates with mortality [79]–[80]. With respect to the latter, our results, limited as they are, indicate that the DME method represents a robust first approximation for deriving masses for somatically immature ontogenetic stages of bipedal dinosaurs. Statistical evaluation of model fitting indicates that there is no significant difference between the constants for growth curves fitted to model-based data and to data generated by applying DME to the minimum body mass of Sue. The DME-generated results also fall within the 95% confidence interval for one of the possible age resolutions for the MOR specimen, though not the other, further bolstering our conclusion. Nevertheless, some caution is warranted as we were forced to use different estimates (i.e., minimum mass value) for Sue as opposed to the other specimens, for which we plotted average values. Although this fact was dictated by the exaggerated reconstruction of Sue's ribcage, we have no means to evaluate how much of an effect that would have. In our favour, however, is the observation that using the minimum value for Sue does result in a growth curve with an asymptote that reflects histological evidence (e.g., the external fundamental system) for a prolonged somatic stasis in Sue. Though more data are required to fully test the precision of DME for T. rex and other taxa, our results are generally encouraging.

Did locomotor performance change during Tyrannosaurus ontogeny? The third aim of our study was to test two related hypotheses about tyrannosaur ontogeny (Was there a relative decrease in hindlimb muscle mass? Did COM position shift cranially?), inspired by ontogenetic changes of relative segment dimensions that are indicated by gross skeletal shapes. Unfortunately, our reconstructed CFL masses and other hindlimb muscle mass data are somewhat ambiguous for juvenile vs. adult tyrannosaurs. Although the upper end of the range of muscle mass estimates for the Jane specimen tended to be relatively larger (as % body mass) than for the adult specimens (especially for the more rotund-bodied models of VA+JM), the ranges overlapped (CFL: ∼2–4% body mass in adults vs. 3–4% in Jane; all limb extensors: 14–25% body mass in adults vs. 19–27% in Jane). Thus, although there is some tentative support for the hypothesis that limb muscle mass decreased during ontogeny (by as much as double across a tenfold change of body mass), the null hypothesis of no difference is not excluded. Because the distal end of the tail is the most frequently missing part of the tail in our specimens, into which the CFL muscle would not have extended, the poor preservation of the tail (discussed above) is not a severe problem for our estimates of CFL muscle volume. Nonetheless, estimates of CFL masses in Tyrannosaurus deserve reconsideration as better specimens become available. However our results cast doubt on the inference of Persons and Currie [19], who used one small Gorgosaurus and one Tyrannosaurus specimen with no sensitivity analysis of error in their estimates, to infer that CFL muscle relative masses increased across tyrannosaurid ontogeny. Their hypothesis is not conclusively falsified either, but is less supported than an ontogenetic increase or lack of change, for the following reasons (other issues with their study were raised elsewhere [81]) in addition to our larger dataset. First, Persons and Currie [19] used smaller body mass estimates (3800–4500 kg) for Tyrannosaurus [81], which inflates the percentage of body mass that their reconstructed CFL muscle seems to represent. For example, if one assumes a body mass of 4500 kg for the Stan specimen then a 261 kg CFL muscle [19] makes up 5.8% body mass. If one instead assumes 7600 kg ([30]; in the middle range of this study's estimates from 5934–8385 kg) then that same 261 kg CFL muscle mass only makes up 3.4% of body mass, or 2.5% if this study's smaller estimate of 192 kg is used. Thus a factor of over 100% error can be introduced into relative muscle mass comparisons by error in body mass estimates. Second, our CFL reconstruction method obtained generally smaller mass estimates for CFL (see example above and Table 8) than Persons and Currie [19] did, for several reasons. These include that we did not also reconstruct the smaller M. caudofemoralis brevis (Fig. 5), and our reconstruction was not constrained to semicircular shapes but rather involved smooth ellipses. Our minimal estimates have a closer match to values obtained with the simpler (25% tail base segment mass) approach of Hutchinson et al. ([28]; ∼141 kg) but are superior in being more anatomically realistic and reasonably well validated. Likewise, the results of our COM position estimates for adult and juvenile Tyrannosaurus are not conclusive, but hint at the potential for a craniad (and dorsal) shift during ontogeny, as noted at the beginning of the Discussion above. If the latter was the case, then this would require greater active muscle volumes around the hip to stabilize the COM and could increase functional demands for muscles acting around other joints [21], [25], just as those muscle masses are potentially decreasing (as above). That would likely result in a steeper decrease of maximal locomotor performance (e.g., running speed) across tyrannosaur ontogeny. Overall, however, the consequences of our muscle and COM estimates for the ontogenetic scaling of locomotion in tyrannosaurids remain ambiguous. Most animals maintain, or slightly decrease, absolute maximal locomotor performance during growth while reducing relative maximal performance (e.g., maximal ground reaction forces sustainable, or maximal Froude number; or minimal duty factor; attainable) [22], [49]–[55], [82]–[83]. Our results are very consistent with these general patterns, and are bolstered by the biomechanical analyses of Hutchinson [21] and the scaling of the limb skeleton [17], [37] across tyrannosaur ontogeny. There is however no strong evidence for whether absolute maximal performance decreased or not, but there is little reason to doubt that relative maximal performance declined. To better resolve this mystery we suggest at least the following are needed: (1) sufficiently large sample sizes of 3D models of tyrannosaurs to properly characterize ontogenetic scaling patterns (once sufficient specimens exist; see [17], [37] for initial skeletal studies); (2) more rigorous validation techniques and sensitivity analysis for estimating muscle and body dimensions (including statistical treatments of these linked to point #1); and (3) improved, well-validated biomechanical models of tyrannosaur locomotion across ontogeny that test the influences of these dimensional estimates on quantitative locomotor performance.