Finite-element analysis was used to investigate the extent of bias in the ichnological fossil record attributable to body mass. Virtual tracks were simulated for four dinosaur taxa of different sizes (Struthiomimus, Tyrannosaurus, Brachiosaurus and Edmontosaurus), in a range of substrate conditions. Outlines of autopodia were generated based upon osteology and published soft-tissue reconstructions. Loads were applied vertically to the feet equivalent to the weight of the animal, and distributed accordingly to fore- and hindlimbs where relevant. Ideal, semi-infinite elastic–plastic substrates displayed a ‘Goldilocks’ quality where only a narrow range of loads could produce tracks, given that small animals failed to indent the substrate, and larger animals would be unable to traverse the area without becoming mired. If a firm subsurface layer is assumed, a more complete assemblage is possible, though there is a strong bias towards larger, heavier animals. The depths of fossil tracks within an assemblage may indicate thicknesses of mechanically distinct substrate layers at the time of track formation, even when the lithified strata appear compositionally homogeneous. This work increases the effectiveness of using vertebrate tracks as palaeoenvironmental indicators in terms of inferring substrate conditions at the time of track formation. Additionally, simulated undertracks are examined, and it is shown that complex deformation beneath the foot may not be indicative of limb kinematics as has been previously interpreted, but instead ridges and undulations at the base of a track may be a function of sediment displacement vectors and pedal morphology.
The sample we have of the body fossil record is notoriously incomplete [1–4] and may be fundamentally biased by environmental and taxon-specific factors that potentially hamper our interpretation of ecological and evolutionary dynamics through deep time [5–8]. Interdependent environmental and taxon-specific biases are equally likely to affect the ichnological, or trace fossil, record. The potential for systematic bias towards ichnofossils produced by larger animals has previously been recognized in the field of vertebrate palaeoichnology, and particularly in the dinosaur track record . Sites preserving only the tracks of very large saurischian dinosaurs (i.e. track lengths greater than 0.5 m) are generally recognized as size-biased assemblages and presumed to represent sedimentary conditions in which only animals above a certain threshold of body mass were capable of producing recognizable tracks . However, beyond this supposition, the influence of body size on the recorded diversity of vertebrate ichnofossil assemblages is poorly understood both qualitatively and quantitatively. It is therefore imperative that the process of track formation and the variables associated with environment and animal biology are investigated.
The relationship between the size of an animal and the load applied to the substrate is not straightforward. Given that pressure is a measure of force over area, the resultant pressure exerted on the sediment surface is a function not only of the animal's mass (as weight), but also the geometry of the autopodia. Quadrupedal animals benefit from more feet in contact with the ground, further reducing the load on the substrate beneath any single foot, as compared with a similar-sized biped. In addition to size, foot morphology also plays an important role in determining the magnitude of deformation expressed as the depth of a track. Differing shapes present different paths for sediment movement, resulting in variable distributions of force that affect the extent to which any given foot may indent a substrate [10,11]. These considerations are not trivial if data on fossil track occurrences and abundances are to be used in ‘higher level’  palaeobiological and palaeoecological inferences. Allen  noted over a decade ago that a widespread understanding of track formation lagged behind knowledge of anatomical aspects and distributions of fossil tracks, and despite a number of rigorous experimental studies in the intervening years [14–20], this still remains the case.
Among vertebrates, the Dinosauria represent a useful model system for studying track formation; their high taxonomic diversity and long evolutionary history yield an array of disparate foot morphologies and a huge range in body mass with which to test for possible biological factors underpinning preservational bias. The group contains small and large obligate bipeds, quadrupeds and supposed intermediate locomotor strategists (e.g. facultative bipedalism [21–23]) that may have exerted different underfoot pressures according to foot geometry and body shape (i.e. mass distribution). Coupled with a vast quantity of research describing dinosaur tracks spanning more than a century and a half [24–26], dinosaurs provide the ideal basis on which to further our understanding of fossil track formation, and the size-related biases associated therewith.
In this paper, information on foot anatomy and mass distribution from osteological evidence and soft tissue reconstructions are integrated with geotechnical theory and computer simulation to explore the potential for size bias in the vertebrate track record. In addition, features related to underfoot pressure and foot morphology were examined in surface and subsurface planes (true tracks and undertracks). Previous work on track formation using computer simulation has explored independently the effects of substrate consistency , foot anatomy  and force . This paper aims to present a combined study in which the quantifiable variables of track formation are considered as a whole system, in the hope of elucidating aspects of preservational bias inherent in the fossil track record.
The following experiments used parallel finite-element analysis (FEA) software developed by authors Margetts and Falkingham, using the freely available ParaFEM libraries (www.parafem.org.uk) to model track formation [11,27–29]. A number of dinosaur tracks were simulated over a range of substrates in order to explore bias in their formation resulting from substrate- or taxon-specific factors.
2.1. The virtual foot
Four dinosaurs (Struthiomimus, Edmontosaurus, Tyrannosaurus and Brachiosaurus) were chosen to create a varied virtual track assemblage on a cohesive substrate, representing a range of body masses, and including obligate bipeds, an obligate quadruped and a facultative biped (table 1). These particular taxa were also chosen because they all have published data on body mass and centre of mass (CM) position [30,31], and represent a wide range in size, mass and autopodial morphology. The taxa were not selected in order to create some geotemporally correct track assemblage.
A track is formed through the interaction of three factors; force, foot anatomy and substrate [32,33]. Force applied and foot anatomy are both dependent upon the track maker. To apply a reasonable force, the body mass for each dinosaur was taken from the literature ([30,31]; see table 1). The reader is directed to Bates et al.  for a comprehensive discussion on the confidence of CM and body mass reconstructions. Animals spend a very small proportion of their time moving at anything more than a walking speed, and it would therefore be expected that most tracks are made by walking animals. Indeed, this is corroborated by the numbers of trackways showing walking, rather than running gaits [34,35]. Stride length is positively correlated with speed [36,37], meaning that at low speeds, the hip joint and CM will move a shorter distance horizontally from the contact between the foot and the ground, resulting in a smaller angle of ground reaction force (GRF) . As such, for the purposes of this paper, a purely vertical component to the applied force was assumed. Force distributed through feet in contact with the ground was taken as the weight of the animal, calculated as mass × gravity (9.81 m s−2). An animal of 100 kg would therefore exert a vertical force upon the ground of 981 N.
For a biped, maximum force is transmitted through a single foot when the opposite foot is raised, so the pressure applied in this case was equal to the weight of the animal divided by the surface area of a single foot. This is to approximate the peak force at any one time during limb support. In the case of quadrupedalism, CM plays a role in determining how much of the animal's weight is distributed to the fore- and hindlimbs, after which the two sets of limbs can be treated separately as bipeds . This is a simplification of the loads experienced by the autopodia of a quadruped during locomotion, but provides reasonable input values for the purposes of this paper. CM estimates for Edmontosaurus and Brachiosaurus were taken from the literature (see  for Edmontosaurus CM and  for Brachiosaurus) and used to apportion force between fore- and hindlimbs. The amount of the animal's weight given to each pair of limbs was equal to the relative position of the CM between the pelvic and pectoral girdles, i.e. a CM 60 per cent of the way from the pectoral girdle to the pelvic girdle would imply a weight distribution of 60 per cent to the hindlimbs, and 40 per cent to the forelimbs  (figure 1). Treating the animal as two linked bipeds with appropriate weights was sufficient for the purposes of these experiments (see  for walking models of quadrupeds represented as two bipeds in tandem). While this may be a very simplified solution that ignores the effects of complex gaits, walking velocities and limb kinematics of dinosaurs are unknown and employing the loading regime outlined above avoids incorporating additional unfounded assumptions into the simulations. Consideration is given to the effects of duty factor and locomotion later in the discussion.
Hadrosauridae have been interpreted as primarily bipedal with facultative quadrupedalism at either low [21,23] or high  speeds, based on anatomical features in the forelimbs suggestive of either mode of locomotion, and trackway evidence also supporting both gait reconstructions [35,39,40]. Edmontosaurus tracks were therefore simulated as being made by both a bipedal animal and a quadrupedal animal.
The indenters, or ‘virtual feet’, were created by producing outlines around ventral views of reconstructed skeletal autopodia (figure 2). Skeletal geometry was scaled to the same size as the specimens used by Bates et al.  and Henderson  so as to remain consistent with mass estimates. The outlines were then increased in size to account for soft tissue. The outline of the Edmontosaurus manus does not follow the osteology as closely as the other indenters, instead being based on the exceptionally preserved hadrosaur body fossil MRF 03 (though scaled to the specimen used by ), as figured in Sellers et al. [38, fig. 5], where the manus soft tissue takes a ‘mitten’-like form over the skeleton. This is supported by hadrosaur manus tracks illustrated by Lockley & Wright , and those described as ‘crescent shaped’ by Currie . The indenters representing the manus and pes of Brachiosaurus were generated as in Falkingham et al.  from reconstructions by Wright . These outlines defined the nodes and elements that would be loaded on the FE substrate volume (figure 2). For each animal, a volume of substrate was created for each foot to be indented into. Only one pes needed to be indented for each bipedal condition, and only one manus and one pes for the Brachiosaurus and quadrupedal Edmontosaurus.
2.2. The virtual substrate
An elastic-perfectly plastic von Mises model was applied in order to model a cohesive clay-like substrate. The mechanical properties of the substrate were defined by the undrained shear strength (Cu), Young's modulus (E) and Poisson's ratio (v). These parameters relate, respectively, to
— The strength of the substrate, that is, how much stress is needed before failure of the sediment (permanent deformation). Essentially a measure of cohesion between grains, shear strength is most strongly affected by water content [10,43]. Typical values of Cu in sediments located on the tidal banks of the Bahia Blanca Estuary, Argentina, were shown to range between 50 and 150 kN m−2 in the surface 1 m . Values of Cu according to the British Standards for Geotechnical Engineering are summarized along with field testing methods in table 2.
— The stiffness of the substrate—how much deformation is recoverable through elastic behaviour before (and after) plastic deformation takes place. The value of Young's modulus is typically 1000× the value of Cu in cohesive substrates .
— The compressibility of the substrate. In an entirely incompressible substrate, v = 0.5. Such a substrate will not change in volume when deformed, resulting in expansion equal to compression along an axis perpendicular to that of the primary stress . An incompressible substrate could be considered to be a fully saturated sediment, in which void space air has been completely replaced by water (note: though water is technically compressible to some extent, at the magnitude of forces dealt with here, it can safely be considered incompressible). Typical values for saturated clay or mud would be 0.4–0.5 .
Many palaeontological FEA studies concerning stress within bone use elastic models, in which there is a linear relationship between stress and strain (figure 3a), determined by E. The introduction of a failure criterion (Cu), however, produces an elastic-perfectly plastic model, whereby initial loading deforms the material in a recoverable elastic manner (line O–Y′ in figure 3b) until the load is sufficient to plastically deform the substrate. Further loading equals or exceeds the bearing capacity and results in failure, where the substrate can no longer support the load (line Y′–P in figure 3b). When the load is removed, recovery occurs along a line parallel to the original elastic deformation (line Y″–U in figure 3b). On the scale of individual elements, this relationship is clear and well defined, but over an entire mesh, where some elements may be in a plastic state and others in an elastic one, the relationship becomes less defined, with a curved portion where plastic deformation occurs (figure 3c).
A soft clay-like substrate was generated in the FEA simulations using 20-node hexahedral elements. The 20-node element is required in this case because of the nature of the deformation; indenting into soft substrate causes a large gradient of deformation from negative vertical displacement beneath the edge of the indenter, to positive vertical displacement adjacent to the indenter. Eight-node elements lack the numerical flexibility to deal with such a gradient, and so by increasing the number of nodes defining the element, a more accurate solution can be found. The volume of substrate modelled was equal to four times the foot length in all dimensions in order to avoid boundary effects .
2.3. The process of indenting
While the foot of an animal can move at joints and the soft tissue is deformable to an extent, as a whole, the foot can be considered rigid compared with the non-rigid substrate. In order to create a rigid loaded area, rigid-body interface elements were generated above the area on the mesh that would be loaded , essentially creating a solid meshed foot on the surface of the virtual substrate. The weight of the animal was then applied to the interface elements to generate a uniform load over the foot. The ‘foot’ was loaded vertically as a static analysis (i.e. independent of loading rate), and then removed vertically in order to allow the substrate to recover the elastic part of the deformation, as would be the case in the formation of a real track. The effects of loading rate are considered later in §4.
For each indenter, substrates were generated with a high Cu, and this was incrementally lowered until the substrate could no longer support the load (i.e. bearing capacity was exceeded). In all cases, E was equal to 1000× Cu, and v = 0.4. Maximum depth of indentation beneath the virtual foot was recorded in each experiment, as this value is a fair indication of the degree to which a track is observable. Additionally, surface tracks and undertracks were visualized and qualitatively observed.
If the body mass and total foot surface area of the models are logged, it can be shown that foot surface area is proportional to mass0.7 (see electronic supplementary material, S1). This is close to the relationship predicted by isometric scaling, where surface area is proportional to mass2/3 [48,49], and suggests that as animals increase in size, the pressure exerted on the substrate (discounting the effects of locomotion) increases at a proportionally greater rate. The small sample size used in this study means that any difference in this relationship between quadrupeds and bipeds cannot be observed, nor can the effects of allometric scaling be explored.
In homogeneous substrates, there is a very narrow range of Cu values for any given pressure that allow the formation of observable surface tracks (figure 4). If Cu is higher than this value, indenters fail to deform the substrate to an appreciable degree, attaining maximum track depths of less than a millimetre. Given that many of the autopodia used were tens of centimetres in length, such deformation cannot realistically be considered to be an observable track. Lower values of Cu than this narrow range cannot support the applied load, and the substrate fails. The maximum load a substrate can support beneath an indenter can be approximated by calculating the bearing capacity beneath a circular indenter under the specified load using the following equation : 3.1where S is a shape factor equal to 1 + 0.2 × (breadth/length).
Using this equation, it can be seen that a substrate for which Cu = 100 kN m−2 will fail when the load on a circular indenter (S = 1.2) reaches 616.8 kN m−2. The approximated failure point for circular indenters at any given load is plotted as a line in figure 4. This prediction is not the true value of bearing capacity for any specific track, however, owing to variations in foot morphology. However, as can be seen from figure 4, this approximation is sufficiently close as to highlight the relationship.
The range of Cu in which tracks of significant depth can be generated is very small regardless of foot morphology or load (figure 5). This limits tracks made in truly homogeneous substrates to an extremely narrow range of pressures and subsequently producer sizes and foot morphology. Below the minimum value of Cu, a track will be formed providing there is a firmer substrate layer beneath. If there is no firmer subsurface layer, the substrate cannot support the load, and the animal in question will be unable to traverse the area without becoming mired.
A key observation is that in simulating track formation in a homogeneous semi-infinite elastic-perfectly plastic substrate, generation of tracks to any significant depth is difficult to achieve—many tracks were so shallow that for real tracks of a similar depth, it would be unreasonable to expect discovery in the field. This is because the load required to plastically deform the substrate and the maximum load that the substrate can support are very close, implying that a very specific pressure is required to generate a track in a homogeneous substrate. There is therefore a ‘Goldilocks’ quality to homogeneous substrates regarding possible track formation. A faunal assemblage represented by tracks at a given tracksite will be strongly biased towards the largest animals the substrate can support, resulting in a very low diversity of recorded body sizes. Taxa exerting more pressure beneath their feet than the substrate can support will avoid the area or become mired, while animals producing less pressure than is required to create a track will not leave observable impressions.
More commonly, substrates are polyphasic, with heterogeneous mechanical properties varying vertically and laterally. If a substrate is underlain by a firmer layer (e.g. compacted sediment or rock), then tracks will be formed if the surface layer fails. If we consider the scenario of a series of substrate layers becoming progressively firmer with depth (figure 6), it is observed that any animal creating sufficient load as to deform the uppermost substrate will generate a track. Refining this stratification such that Cu increases gradually with depth results in the intuitive case that heavier animals generate deeper tracks. It can be seen from figure 4b that this being the case, there is a much larger range of possible track-bearing substrates for animals exerting a greater pressure (i.e. by being larger or moving faster). Given the shallow nature of the tracks modelled in homogeneous substrates, it becomes apparent that most real tracks must therefore be formed in mechanically heterogeneous substrates, or in relatively shallow homogeneous substrates underlain by rock.
If a substrate is stratified with mechanically distinct layers, then the depths and surface areas of present tracks can be used to infer the depth and mechanical properties of these layers at the time of track formation. A track-bearing surface on which small and medium tracks are impressed to a similar depth, yet on which tracks made by larger animals appear deeper, will indicate a mechanically homogeneous surface layer as deep as the small and medium tracks (after accounting for subsequent weathering/erosion). Such consideration of tracks as palaeo-penetrometers may prove useful in interpretations of palaeoenvironment, at least in so far as determining substrate conditions at the time of track formation.
Examining the biases inherent in track formation as a consequence of animal size permits discussion of genuine and artificial signals regarding diversity as interpreted from tracksites. A track site limited to large producers, e.g. a purely sauropod track assemblage, is likely to be a preservational artefact, or indistinguishable from such. Smaller animals may have been abundant at such a site, but unable to produce tracks in the substrate. Allen  noted in his discussion of the Flandrian deposits of the inner Bristol Channel and Severn Estuary that the fauna represented by tracks lacked records of smaller mammals such as foxes and dogs, arguing this was a preservational rather than an ecological issue. The results presented here support this hypothesis, and may be applicable to other fossil tracksites dominated by large fauna such as Fumanya, Spain [52,53], or the tracksites at Paluxy River, Glen Rose, TX, USA . In the case of the Glen Rose tracks, the site is dominated by large sauropod and medium to large theropod tracks. One trackway has been interpreted as showing the interaction of a theropod and sauropod , implying that the trackways were contemporary. Given the considerable depth of both the theropod and sauropod tracks, and the lack of tracks from smaller animals, the trackways appear to be consistent with the Goldilocks effect. Both the theropod and sauropod exceeded the bearing capacity of the surface mud, and indented deep tracks until supported by firmer substrate layers beneath the surface. However, the depth of the soft surface mud may have been too great for smaller animals to safely traverse while leaving tracks, resulting in the formation and subsequent preservation only of the largest animals present. A similar case can be made for the sauropod trackways at Fumanya, where deep sauropod tracks dominate. Shallow theropod tracks were previously reported from the site, but have since been subjected to weathering and are no longer present . In this case, only larger animals were able to produce tracks, resulting in an impoverished track assemblage, whose low diversity has been exacerbated further by weathering and the complete removal of the shallower tracks left by smaller animals.
Given that there is a strong bias towards greater underfoot pressures, the preservation potential of track assemblages representing mixed age groups (herd behaviour) is greatly reduced; there will be a strong bias towards preserving only the largest members of the group. If the adults within a group are particularly large, as in the case of sauropods for instance, the range of substrates traversable by the group will be constrained by the minimum substrate strength that can support the largest animals. As such, tracks from smaller individuals become far less likely to form and subsequently preserve, because the substrates over which the animals move may not be soft enough to record the passage of smaller, juvenile forms. This reduction in preservation potential of mixed-age herds is supported by the fossil record; Myers & Fiorillo  noted that of 13 sauropod trackway associations indicating gregarious behaviour, only three sites contained tracks from multiple age groups.
The presence on a single track-bearing surface of both small and large true tracks (rather than transmitted or undertracks), indented to approximately the same depth (evidently halted by a firmer subsurface layer), is likely to be more indicative of true diversity in the area at the time of track formation (providing effects of time averaging can be removed). Presence of small, shallow tracks and large, deep tracks may not be indicative of true diversity, however, if the large deep tracks are particularly deep. In such a scenario, it is possible that medium-sized animals produce too great a pressure underfoot to be supported by the soft surface layer, but sink too far before reaching a supportive layer as to be able to traverse the area. This highlights the importance of considering not just the tracks present at a tracksite, but their total three-dimensional morphology, including foot anatomy and track depth, in order to make interpretations about faunal diversity.
The Goldilocks effect described in this study is mitigated when a substrate is exposed for a period of time during which the mechanical properties alter, such as when a substrate is drying out. Changes in mechanical properties will undoubtedly be the rule, rather than the exception, but the rate at which these changes occur will determine the overall applicability of the Goldilocks effect to the fossilized tracksite. In cases where the substrate dries out over a relatively prolonged period, the Goldilocks effect will be applicable over short time spans, but will not be evident over the recording life of the substrate, or in the preserved track surface. The presence of sedimentary features such as drying cracks or displacement rims that are unique to some trackways and not others may shed some light on the preservational context, and as to whether the Goldilocks effect noted here applies to a given track assemblage.
The experiments carried out in this study have used body mass to apply a force through the autopodia in a number of taxa. Loading in this way assumes a direct relationship between body mass and force, and was carried out as a rate-independent (i.e. static) analysis as the focus of the work was to explore bias relating to size, not necessarily locomotor mode. This was done to avoid incorporating unfounded assumptions into the simulations, given that the habitual gaits of dinosaurs are unknown. Nevertheless, consideration must be given to the effects of locomotion, duty factor and limb kinetics and kinematics. As an animal begins to move, the GRF gains a horizontal (forward–backward) component in order to move the animal forwards . This force vector may also incorporate a lateral component depending on the animal's gait. As speed increases, the magnitude of the GRF also increases. In terms of pressure applied, the pressure beneath an animal's foot will increase as speed increases. As an animal increases in speed, the minimum Cu required to support the load also increases, such that a substrate that previously would be incapable of recording an animal standing or moving slowly may fail beneath the foot of a running animal.
As an animal traverses a substrate, the rate of loading is intrinsically linked to the speed and duty factor of the animal, with loading rate increasing as duty factor decreases. An increased loading rate results in more resistance from the substrate, and the result is that deformation occurs to a lesser extent (see electronic supplementary material, S2). However, at higher speeds, an animal exerts a greater force upon the substrate, as noted above. As such, although the loading rate increases with speed, the effects on substrate displacement will be mitigated by the increased load applied by the foot. It is important to note that the analyses carried out in this study were static, and thus did not account for the effects of rate-dependent loading and thixotropy. Exploring this complex interplay between loading magnitude and loading rate is beyond the scope of this study, and is impossible without comprehensive locomotor reconstructions of the animals in question. Instead, the Goldilocks effect can be considered to be the base mechanic around which other confounding factors such as limb dynamics, locomotion and substrate thixotropy have an effect. To provide some insight into these complex issues, and their relationship to the Goldilocks effect, the results from a series of simple dynamic simulations are provided and discussed in electronic supplementary material, S2.
Given the static vertical loading conditions employed in this study, the relationship between size and track-forming potential could be predicted to a reasonable degree with simple mechanics and geotechnical theory, as evidenced by the close correlation between equation (3.1) and the results (figure 4a). However, FEA provides benefits over simple mechanics. Simulating track formation allows for differing foot morphologies to be tested, which was shown to be a potentially important factor by Falkingham et al. . FEA also allows us to explore the full three-dimensional volume of simulated tracks. The importance of understanding the relationship between foot morphology, three-dimensional deformation and undertrack production will be demonstrated by the discussion of features observed in the models.
5. Discussion of individual track features
The simulations undertaken for this study present an opportunity to investigate track features at the original track surface, and in subsurface undertracks. Specific features related to autopodia morphology and undertrack depth are discussed here. In order to generate deeper tracks and associated undertracks, failure was allowed to occur, but was halted when maximum depth reached 0.05 m, approximating a firm subsurface layer.
The values of mass, foot morphology and CM position used for the Edmontosaurus produce differing pressures between manus and pes (151.93 kN m−2 and 103.31 kN m−2, respectively). These differing pressures imply that substrates of different Cu are required to support the loads, which in turn creates a range of substrates (Cu = 20–40 kN m−2), where the manus causes the substrate to fail, but the pes does not. In such substrates, if underlain by a firmer layer, only the manus will generate tracks. This is the same mechanism as described in detail by Falkingham et al.  for sauropod manus-only trackways. It is interesting to note that the resulting pressure beneath the pes when bipedal locomotion is assumed is less than for the manus in a quadrupedal mode of locomotion. Depending on the validity of the mass, CM and foot outline input parameters used here, this may support the hypothesis that differing modes of locomotion may have been used for traversing different substrates. The values employed in this study would suggest quadrupedalism to be potentially more advantageous on firmer substrates, where the manus will not sink, while a bipedal mode of locomotion would allow traversal of softer substrates. However, this suggestion is based only on underfoot pressures, and further factors such as stability will ultimately determine gait. It may therefore be unwise to infer otherwise unknown locomotor styles from trackways in which the substrate conditions at the time of track formation cannot be constrained. The reader is directed to Wilson et al.  for discussion of a trackway in which an ornithopod trackmaker transitions between bipedal and quadrupedal gaits as substrate changes.
As described above for the Edmontosaurus track simulations, and as described by Falkingham et al. , there are a range of substrates in which only the manus, and not the pes, of Brachiosaurus produce enough pressure to deform the substrate. When the subsurface undertracks generated by manus and pes are visualized, important features can be observed. The bowl-like form of the Brachiosaurus pes (figure 7) is reminiscent of a number of reported sauropod tracks (e.g. [59,60]). The simulations here indicate that such a bowl-like form is potentially characteristic of undertracks. This is consistent with the assumption that the plantar surface of the foot was approximately flat, given that the shape of foot required to form bowl-like tracks would prove unstable on firm ground.
When the Brachiosaurus manus track is observed as a subsurface undertrack at approximately 0.2 m depth, a ridge running transversely across the track can be seen (figure 8a,b). This ridge appears superficially similar to the undulating track surface hypothesized to result from three-phase movement of the foot [16,61]. However, with full control of all input variables, it is known that in this case the loading was carried out in an entirely vertical manner, evenly distributed through a flat foot, and so the ridge cannot be a function of limb kinematics or foot anatomy. Instead, this ridge is produced through the displacement of sediment according to Prandtl theory . As substrate is deformed by a load, it is pushed down and out from beneath the indenter (figure 8c). The base of the actively deforming zone of substrate undulates against the rigid, non-moving zone . A cross section through this area results in a subsurface track containing a ridge of non-deformed substrate. This effect is seen in the Brachiosaurus manus track because of the round shape of the indenter. This ridge is very subtle, and its absence from the other simulated tracks implies that its formation is closely linked with indenter morphology. Dynamic limb motion consisting of more complex, non-vertical loading will deform this ridge accordingly (e.g. see theropod track in ); however, fossil track evidence indicates that, as is the case here, sauropods placed the manus vertically, at least when traversing soft substrates . Further study is required to fully understand this phenomenon and to avoid erroneous interpretation of tracks with undulating bases. Note that some previously described ridges in sauropod tracks (e.g. ) appear more defined and morphologically different to those outlined here, and we do not suggest this mode of ridge formation for those cases.
Both theropod tracks (Struthiomimus and Tyrannosaurus) indented to a considerably greater degree at the posterior of the virtual foot (figure 9). The pes of the Edmontosaurus also exhibited this feature, albeit to a lesser extent. This effect is a function of the shapes of the indenters as seen in Falkingham et al. . The appearance of a deeper posterior track portion under uniform loading of a flat indenter has important consequences for interpretations of limb kinematics from fossil tracks. Commonly, the morphology of real tridactyl tracks is deeper beneath the distal areas of the digits as a result of the increased pressure as the animal kicks off, but in the simulated case there is no such loading regime. The development of a deeper track beneath a larger, more compact part of the foot would mean that the ‘two-phase’ interaction of the foot (weight bearing and toe off) described by Thulborn & Wade  could potentially produce a track with the appearance of a ‘three-phase’ foot–substrate interaction (which precedes the above phases with touch-down), where the heel and toes are deeper than the centre portion of the track [16,64]. It has been proposed that the two-phase and three-phase modes of locomotion represent knee-based and hip-based retraction of the limb, respectively ( and references therein), and that the associated pressures across the foot differ accordingly. However, if vertical loading, as has been used here, produces a deeper ‘heel’ in tridactyl tracks without a ‘heel-down’ kinematic phase, then attempting to differentiate between the locomotor modes of theropod dinosaurs and birds from fossil tracks may be considerably more difficult than has previously been assumed without ongoing experimental studies.
The resultant tracks from these simulations are relevant for studies of tracks where interpretations of locomotion have been made based on the ‘pitch’ of the track [57,65]. Such interpretations must consider as an alternative, or at least confounding factor, varying shear strength (as a function of water content) throughout the total substrate layer at the time of track formation, resulting in track pitch altering systematically along a trackway. Alternatively, variations in track ‘pitch’ may be influenced by grain size or compositional differences, given that sand responds in the opposite manner to mud , allowing greater deformation beneath digits.
The simulation of tracks from a series of dinosaur taxa ranging in size from 400 to 25 000 kg shows a linear relationship between body mass and substrate shear strength required to produce observable tracks. The point of failure for a given track and subsequently the shear strength of the substrate at the time of track formation can be approximated by calculating the bearing capacity required for a circular indenter of equal size and load. Variations around this approximation are due to the effects of foot shape.
Tracks of significant depth are not possible in homogeneous, cohesive substrates without the presence of a firmer subsurface layer, because failure of the substrate will result in the animal being unable to traverse the area. A homogeneous cohesive substrate will only record tracks from the largest animals that substrate can support without failing. There is, however, a strong bias towards tracks made by larger animals if there is a firmer substrate beneath a softer layer. This Goldilocks effect means that for a homogeneous substrate, loading conditions (that is, the animal size, locomotion and foot morphology) must be ‘just right’ in order for the animal to be able to traverse the area but still form tracks. This has wide-ranging implications for interpretations of palaeodiversity and palaeoecology based on vertebrate track assemblages preserved in lithified muds and silts.
Presence of small and large tracks indented to the same depth on a single track-bearing surface (assuming time-averaging/transmitted tracks can be accounted for) offer the highest possibility of presenting a true representation of faunal diversity in the area at the time of track formation. Caution is strongly advised in making any interpretations of faunal diversity or population dynamics from track assemblages where all tracks have been produced by similar-sized producers. Such assemblages most likely represent a strongly biased preservation, or an ‘instantaneous’ event. Substrates which have dried out over relatively long time periods will provide a fuller record of faunal diversity in the area, but will be subject to biases and other erroneous data associated with time-averaging.
Specific features regarding track and undertrack formation have been noted for this range of indenters based on dinosaur taxa. Bowl-like sauropod pedal impressions may be indicative of being undertracks, potentially of significant depth. Internal ridges may form in tracks due to the vectors of displacement beneath a uniform load (as in the Brachiosaurus manus), or as a result of autopodia morphology causing non-uniform displacement under uniform loading (as in the tridactyl tracks). That these features can be formed independent of limb kinematics is of great importance, and highlights the need for further experimental work to clarify the specifics of their formation.
The approach used here, of computer simulation using FEA, has allowed the generation of tracks and associated undertracks for a range of animal sizes that would be difficult to replicate using physical modelling. Employing computational methods has also catered for constancy in input variables between experiments, and has provided the ability to easily and systematically manipulate those variables. We recognize that this study makes a number of assumptions and simplifications in terms of loading, and expect subsequent research to build on the methods used here to produce more complex models. Many of the conclusions and observations recorded here are related to the mechanics of substrates under load, and we hope that this will encourage further research into the effects of complex limb kinematics and kinetics on track formation, in light of the confounding geotechnical effects described here.
P.L.F. and K.T.B. were funded by the Natural Environment Research Council (NER/S/A/2006/14033 and NER/S/A/2006/14101, respectively). FEA simulations were run on the HPCx supercomputing service, using Engineering and Physical Sciences Research Council grant EP/F055595/1, awarded to L.M. We would also like to thank James Jepson for commenting on an early draft of the manuscript, and also Jeff Wilson and one anonymous reviewer for their constructive, helpful comments.
- Received November 15, 2010.
- Accepted December 21, 2010.
- This Journal is © 2011 The Royal Society