Large size in aquatic tetrapods compensates for high drag caused by extreme body proportions

Aug 11, 2023

Communications Biology volume 5, Article number: 380 (2022) Cite this article

5635 Accesses

4 Citations

278 Altmetric

Metrics details

An Author Correction to this article was published on 19 May 2022

This article has been updated

Various Mesozoic marine reptile lineages evolved streamlined bodies and efficient lift-based swimming, as seen in modern aquatic mammals. Ichthyosaurs had low-drag bodies, akin to modern dolphins, but plesiosaurs were strikingly different, with long hydrofoil-like limbs and greatly variable neck and trunk proportions. Using computational fluid dynamics, we explore the effect of this extreme morphological variation. We find that, independently of their body fineness ratio, plesiosaurs produced more drag than ichthyosaurs and modern cetaceans of equal mass due to their large limbs, but these differences were not significant when body size was accounted for. Additionally, necks longer than twice the trunk length can substantially increase the cost of forward swimming, but this effect was cancelled out by the evolution of big trunks. Moreover, fast rates in the evolution of neck proportions in the long-necked elasmosaurs suggest that large trunks might have released the hydrodynamic constraints on necks thus allowing their extreme enlargement.

Tetrapods have undergone multiple independent transitions from land to sea1,2, each associated with major body plan innovations1,3,4. During the Mesozoic, numerous reptile clades took to the seas, becoming secondarily aquatic5, of which the ichthyosaurs and sauropterygians, which emerged following the end-Permian mass extinction, were the longest-persisting and most diverse6. Derived members of both groups became pelagic, a shift in ecology coupled with the evolution of medium to large body sizes7 (i.e. body lengths above 2 m), fast metabolic rates and the ability to thermoregulate6,8,9,10. Additionally, they became specialised, lift-based swimmers, but achieved this in different ways: ichthyosaurs as caudal oscillators, with stiff, deep bodies and broad lunate caudal tails11,12; and derived sauropterygians (i.e. plesiosaurs) as quadrupedal underwater fliers, with expanded girdles and large and rigid hydrofoil-like flippers13,14. The bodies of modern cetaceans are adapted for endurance and speed15,16 and have often been used as a functional reference for ichthyosaurs7,17,18. However, the swimming mode of plesiosaurs is unique among tetrapods, and no living analogues possess similar body plans. Consequently, although in recent years the biomechanics of plesiosaur swimming has been the subject of extensive research14,19,20, important aspects of their locomotory biology remain enigmatic.

Two distinct plesiosaurian body plans, the short-necked pliosauromorphs and the long-necked plesiosauromorphs, emerged independently in various clades21. The most extreme body proportions are found in the Cretaceous elasmosaurs, some of which had necks up to 6 m long22,23. Previous research suggested that ichthyosaurs could reach faster cruising speeds than plesiosaurs of the same body length7,24, and that long-necked plesiosauromorphs were slower than the short-necked morphotypes7,25. These differences were based largely on assumptions of a less efficient swimming mode in plesiosaurs and, to a lesser extent, on differences in their fineness ratios (FR, the proportion of maximum body depth to total length). However, the effect of whole-body morphology on drag is not yet fully understood.

Certain streamlined axisymmetric geometries produce the lowest drag for a given volume at a fineness ratio of 4.526,27, which led to an assumption that FR = 4.5 is the optimal proportion for low drag in aquatic animals7,28. Consequently, long-necked plesiosaurs and other elongated aquatic reptiles possessing suboptimal fineness ratios, have been classified as ‘slow’7,25. However, recent computational fluid dynamics (CFD) analysis in ichthyosaurs showed that drag is not correlated with fineness ratio for bodies of equal volume or mass29. Therefore, here we question this association between FR and drag also for plesiosaurs. Long necks have also been argued to add extra viscous drag due to their large surface area as well to increase pressure drag7,25,30. A recent CFD-based study of plesiosaurs concluded that drag was not affected by neck length during forward motion20. However, measures of skin friction (i.e. viscous) and pressure drag, the two main components of drag in fully submerged swimming, were not assessed, nor was the impact of the neck on the balance of drag to body mass20.

Here we take derived ichthyosaurs and plesiosaurs as paradigm models for high axial and appendicular locomotory specialisations in marine tetrapods and compare them to modern cetaceans using digital modelling and CFD. Our CFD protocol informs on skin friction and pressure drag components29, essential when assessing the flow over slender bodies for which drag is mostly frictional31,32. This allows us to address several questions. Did the bodies of both groups reach a similar level of low-drag form? Did the limbs contribute differently to drag in these groups? How do they compare to modern cetaceans? And importantly, how does body size influence the effects of shape? Then, we focus on plesiosaurs, exploring differences between the two morphotypes and the interplay between body size and body proportions for derived plesiosaurs, which show a great spread in fineness ratios due to high neck length plasticity23. Lastly, we explore the evolution of trunk length (used here as a proxy for body size) and neck proportions in Sauropterygia (i.e., plesiosaurs and their closest Triassic relatives) and analyse the effect of neck plasticity (i.e. the variability of neck to trunk ratio) on the drag-related costs of steady swimming, discussing functional and ecological implications.

At equal Reynolds numbers (same body length and same flow velocity), the total drag coefficients of plesiosaurs (Cd) are higher than the estimated values for ichthyosaurs and modern cetaceans (Fig. 1a). The limbless bodies, however, display similar Cd in all three groups and are even lower-than-average in the long-necked plesiosaurs, indicating that the limbs are responsible for the observed high Cd. The limbs of plesiosaurs contribute to more than 20% of their total drag coefficient: up to 32.2% in the basal Meyerasaurus and averaging 25% in derived plesiosaurs, with no major differences between plesiosaur morphotypes. In parvipelvian ichthyosaurs the contribution of the limbs to Cd is 11.2–15.6%, compared to 8.7–14.3% in modern cetaceans. Some of the living taxa we include provide a functional reference for this analysis. Our computed drag coefficient for the bottlenose dolphin model (Cd = 0.00413 at Re = 107) for example, is consistent with the estimates from a gliding living dolphin33 (Cd = 0.0034 at Re = 9.1 × 106) and other static CFD simulations34 (Cd = 0.00413 at Re = 107). It is worth noting that these values are, as expected, lower than estimates obtained from kinematic models, as motion is not accounted for35. In a former study, drag coefficients for a plesiosaur (Cryptoclidus), two ichthyosaurs and various cetaceans were obtained from rigid models in water tanks36. However, the pressure drag component (Cp) was likely overestimated due to the proximity of the models to the air–water interface, and thus are not directly comparable to ours.

a Total drag coefficient computed for the full models including the limbs (‘body + limbs’, circles) and the limbless models (‘body’, squares). Average (point) and range (bar) shown for calculations at Re = 5 × 106–107. The derived short-necked plesiosaurs are highlighted in orange; the parvipelvian ichthyosaurs in blue and the extant cetaceans in red. A basal plesiosaur included as a reference is highlighted in purple. b Representative two-dimensional plots of the flow velocity magnitude at Re = 5 × 106 (inlet velocity of 5 ms−1) in lateral view. For dorsal view see Supplementary Fig. 1. Images of Tursiops and the three ichthyosaurs modified from Gutarra et al.29.

In all models across the various clades, velocity plots display a stagnation point at the anterior tip of the model, a thin velocity gradient along the body corresponding to the boundary layer, an area of higher velocity around the greatest diameter and a low velocity wake behind the body, characteristic features of a fully developed external flow (Fig. 1b, Supplementary Fig. 1). The acceleration of flow results in areas of low pressure (Supplementary Fig. 2), while high pressure areas are observed where stagnation occurs. Our CFD methodology has been previously validated against experimental data from slender torpedo-like shapes26 and has been shown to provide a reliable distribution of internal drag components29 essential when dealing with streamlined bodies35. In all our simulations, the proportion of frictional and pressure drag was consistent with the expected values for slender geometries31: most of the drag originated from skin friction with a minor pressure drag component (Supplementary Fig. 2). The relatively larger limbs of plesiosaurs (Supplementary Table 1) produce a small increase in skin friction (Supplementary Fig. 2a), but a large increase in the pressure drag coefficient (Supplementary Fig. 2b), indicating that the latter largely explains differences in total drag coefficient between the groups. These effects might be explained by the low local Reynolds number of the flippers (resulting from a small chord length) producing high local Cd relative to the rest of the body31, alongside interference drag (i.e. drag caused by the interaction of flow fields where limbs and body meet), which might be higher for larger flippers.

When comparing morphologies at the same volume (proxy for body mass) and the same velocity, to focus on the effect of shape alone, derived plesiosaurs produce on average 30% more drag than parvipelvian ichthyosaurs and modern cetaceans (Fig. 2a, Supplementary Table 3; two-sample t-tests p < 0.001). Drag-per-unit-volume represents the contribution of drag to the cost of transport (COTdrag), with COT being the mass-normalised effort required for sustained forward swimming37. As for the drag coefficient, these differences are observed only when the full morphology is considered and not in the limbless models, indicating that the differences are caused by the relatively larger limb sizes in plesiosaurs. The model with the lowest absolute value of COTdrag is Tursiops, against which all other taxa were normalised. The highest COTdrag was estimated for the basal plesiosaur Meyerasaurus, which generates about 69% more drag than a bottlenose dolphin of the same mass. Among derived plesiosaurs, drag values increase from 29.2% in Thalassomedon to 42.6% in Dolichorhynchops relative to an equal-mass Tursiops, and no substantial differences are observed between the short-necked and long-necked morphotypes. The estimates of COTdrag in parvipelvian ichthyosaurs are about 4% to 15% higher than for the Tursiops model, very close to our estimates for the modern cetaceans Orcinus and Megaptera, which have relatively large fins. Overall differences between parvipelvian ichthyosaurs and cetaceans are non-significant (two-sample t-test p = 0.63; Supplementary Table 3).

a Relative drag per unit of volume (a proxy for the drag-related cost of steady locomotion or COTdrag) calculated for models scaled to the same total volume and compared at the same inlet velocity of 1 ms−1. Results are shown for the full models including the limbs (circles) and the limbless models (squares). Average of calculations performed with two different volumes (see Supplementary Data). b Relative drag per unit of volume for life-size scaled models compared at the same inlet velocity of 1 ms−1. Error bars represent minimum and maximum values accounting for taxon body size variation (see Supplementary Data). For an alternative set of calculations at 2 ms−1, see Supplementary Fig. 3. c–f Relative values of drag per unit of volume for models scaled to the same volume and measured at the same inlet velocity of 1 ms−1, corresponding to results in a, plotted against the fineness ratio, FR (c, e) and the surface area-to-volume ratio (d, f). Results are shown for limbless (c, d) and full (e, f) models. All values are normalised to the results for the Tursiops model. Derived short-necked plesiosaurs are highlighted in orange; the parvipelvian ichthyosaurs in blue and the extant cetaceans in red. A basal plesiosaur included as a reference is highlighted in purple.

Our CFD-based analysis thus shows that the overall morphology of plesiosaurs produced higher drag than parvipelvian ichthyosaurs and modern cetaceans, meaning that all other things being equal, an ichthyosaur should endure longer swims at a given speed or cruise at a faster velocity than a plesiosaur of the same mass. It is, however, uncertain to what extent all other things were equal. Propulsive efficiency estimates from living caudal oscillators such as cetaceans are generally higher than those of underwater fliers such as penguins, turtles and sea lions (for data of efficiency in extant animals and their sources see Fish, 200638). However, plesiosaurs were quadrupedal swimmers, with no functional reference among living tetrapods, and recent work suggested that their propulsive efficiency was enhanced by fine-tuning of the fore and hind flippers14. Further, plesiosaurs have more surface area dedicated to producing thrust. Thus, whether a more efficient propulsion compensated the extra drag of the flippers in plesiosaurs is not yet known.

Our results show there is no correlation between the COTdrag and the body fineness ratio (FR), regardless of whether the limbs are included or not (Fig. 2c, e). These observations, although opposite to what is generally assumed for aquatic animals7, are consistent with previous analyses in ichthyosaurs29. The widely extended concept of a FR range for minimum drag comes from the study of aerodynamically engineered forms, and only applies to certain shapes when all other geometric parameters are kept constant26,27, but cannot be extended to all complex streamlined forms. Instead, COTdrag displays a strong positive correlation with the ratio of surface area to volume only in simulations with full morphology (Pearson’s product-moment correlation, r2 = 0.89, p = 4.11 × 10−9; Fig. 2f). This is consistent with the expected behaviour of flow over streamlined forms for which drag is mainly frictional31. The large hydrofoil-shaped limbs in plesiosaurs, necessary for their lift-based quadrupedal appendicular swimming13,14,39, contribute to a large fraction of the surface area without adding much volume (Supplementary Table 1). In contrast, parvipelvian ichthyosaurs and modern cetaceans, both caudal oscillators, have lower proportions of body surface dedicated to the limbs (Supplementary Table 1), which add very little drag relative to a limbless body.

When body size is incorporated into the analysis (i.e. assessing the combined effect of shape and size by simulating the flow of life-size models for a constant velocity of 1 ms−1), the group differences detected in the volume-scaled simulations disappear (Fig. 2b, Supplementary Table 3; all two-sample t-tests p > 0.05). In these conditions, the drag-related costs of steady swimming of plesiosaurs fall within the range observed in both modern cetaceans and ichthyosaurs. Normalised against a 2.85 m-long Tursiops, the COTdrag for derived plesiosaurs ranges from 0.42, estimated for the large elasmosaur Thalassomedon, to 1.41 in the medium-sized Dolichorhynchops. In the parvipelvians, COTdrag spans from 0.33 estimated for the large Temnodontosaurus, to 1.76 in a 2.5 m-long Stenopterygius. Cetaceans show a smaller lower limit, because they include the largest animal in our sample, a 16 m-long humpback whale, with a COTdrag of 0.13 compared to Tursiops. The estimated cetacean upper COTdrag limit is 1.54 for a 1.9 m Tursiops. On the other hand, comparisons of the total drag power (Pdrag, i.e., the non-mass normalised version of COTdrag) for the same speed of 1 ms−1 (Fig. 3), show a different trend. Pdrag is highest for Megaptera, higher than in any fossil taxa included in this study, and is lowest in Tursiops. Thalassomedon is comparable both in total drag power and COTdrag to the killer whale. Similarly, the thalassophonean pliosaurid Liopleurodon matches the elasmosaurian Hydrotherosaurus in having a similarly low mass-normalised COTdrag but requiring about 4× more total drag power than Tursiops. Smaller forms like the polycotylid Dolichorhynchops and the thunnosaurian Ophthalmosaurus resemble the extant bottlenose dolphin in having a relatively high COTdrag and low total power.

Values of mass-normalised drag power (i.e., drag per unit of volume or COTdrag calculated as in Fig. 2b) in grey, and non-mass-normalised total drag power, in black, for an array of derived plesiosaurs, parvipelvian ichthyosaurs and modern cetaceans compared at the same inlet velocity of 1 ms−1. Error bars represent minimum and maximum values accounting for taxon body size variation (see Supplementary Data). Values are normalised to the results for Tursiops.

Thus, in contrast to the volume-normalised simulations, differences between animals at their life-size scale are mainly influenced by size. For example, medium-sized plesiosaurs and ichthyosaurs, such as Dolichorhynchops and Ophthalmosaurus, have values of COTdrag close to that of a dolphin, while large plesiosaurs like Thalassomedon are more like the parvipelvian ichthyosaur Temnodontosaurus and a modern Orcinus. It is worth noting that the inflow velocity of 1 ms−1, is a reference velocity used for comparative purposes, and is not equivalent to the optimal cruising speed (i.e. speed at which COT is minimum16). This parameter is known to vary little in nature, with most vertebrates displaying values of preferred speed between 1–2 ms−1 regardless of body size40,41,42, which means it is reasonable to assume all tested taxa, regardless of their size, were able to swim at this velocity. Using a different reference velocity (2 ms−1) has no effect on the relative values of drag per unit of volume and the mass-normalised drag power (Supplementary Fig. 3; Supplementary Data). A reduction of mass-normalised drag-related costs of cruising as body size increases is selectively advantageous, as energy savings can be used to extend foraging and mating range, increase swimming speed and fuel other activities42,43.

Our analysis shows that for highly aquatic tetrapods, size dominates over shape in affecting the drag-related costs of steady locomotion. This is because COTdrag (i.e., the balance of drag to volume) is highly sensitive to surface/volume proportion (Fig. 2f), and so is much influenced by isometry in streamlined animals.

Simulations at constant Reynolds number (i.e., comparing models at same total length and same flow velocity), show that necks up to 5× the length of the trunk do not increase substantially the total drag coefficient. Longer neck ratios up to 7× were found to impact the drag coefficient by as little as 3% (Fig. 4a). We estimated a 4–10% increase in skin friction drag coefficient for neck ratios of 3–7×, but also a comparable reduction in pressure drag resulting in almost no change in the total drag coefficient. A previous CFD-based study also found no differences in drag coefficient between plesiosaur models with variable neck proportions20, but further comparison is not possible because of great differences in the order of magnitude of Cd, the use of a different scaling reference area and the lack of information on skin and pressure drag20. Here, we have shown that long necks produce only a small increase in skin friction, although not as great as previously speculated25,30, and this is nullified by reduced pressure drag.

a Total drag coefficient and skin friction drag coefficient for an array of hypothetical plesiosaurs with varying neck ratios computed at Re = 5 × 106 (same total length and inflow velocity). b Drag per unit of trunk volume computed for the same array of models scaled at the same trunk length and tested at the same speed of 1 ms−1. The hypothetical models were created by modifying the length in the model of the basal plesiosaur Meyerasaurus victor which has a neck ratio of 0.87×. The limits of the trunk (which extends along the torso and includes the edges of the pectoral and pelvic girdles) are shown in red in the rendered models. c Three-dimensional models of a wide array of plesiosaurs, in dorsal view, at their life-size dimensions, showing the differences in body proportions and sizes. The limits of the trunk in the models (defined as in b) are coloured by group. Basal plesiosaurs are highlighted in purple. Among the derived groups, thalassophonean plesiosaurs (derived pliosaurid plesiosaurs) are highlighted in light orange, polycotylid plesiosaurs in dark orange and elasmosaurid plesiosaurs in green. d Scatterplot of trunk length (cm) and neck ratio showing the relative drag per unit of trunk volume as a gradient of colour for each taxon analysed and for the plot area in between (contour lines represent the interpolated values of drag per unit of volume). e Plot of the relative drag per unit of trunk volume versus the trunk length showing results highlighted by group. Line plots at the right-hand side show the range for each group. The D/Vtr and the trunk length show a significant negative correlation (Pearson’s correlation coefficient calculated with log-transformed variables, p = 2.28 × 10−7, R2 = −0.92). A small version of the fitted power curve (regression equation \(y=69.76{x}^{-0.94}\)) is shown on the right upper corner. The grey area around the curve represents a confidence interval of 95%. All values in b, d and e are normalized to the results for the Meyerasaurus model.

Next, we explored the impact of neck proportions on drag-related costs of swimming in simulations where the size factor is removed. We found that if trunk dimensions are kept constant while the neck is enlarged, the drag per unit of trunk volume does not change appreciably for neck ratios up to 2×. However, longer neck proportions did impact resistive forces. This was moderate for a 3× ratio, with 12% more drag per unit of trunk volume, but became more substantial for longer necks, with 22%, 35% and 59% excess drag for necks of 4×, 5× and 7× respectively (Fig. 4b). This means that elasmosaurine elasmosaurs, with necks commonly 3–4× the length of the trunk23 might have experienced higher drag than other plesiosaurs of similar trunk dimensions.

To test if the ‘long neck effect’ remains when body size is accounted for, we compared the relative amount of drag-per-unit-trunk-volume (D/Vtr) in a wide sample of plesiosaurs (Fig. 4c) at life-size scale for a constant velocity of 1 ms−1, including three species with neck ratios above 2×: Styxosaurus (2.76×), Hydrotherosaurus (3.18×) and Albertonectes (3.72×), the last being the elasmosaur with the longest reported neck44. Our results show great variability in D/Vtr. Small-bodied plesiosaurs such as Plesiosaurus, Meyerasaurus and Dolichorhynchops generated up to six times more D/Vtr than the largest plesiosaurs, Kronosaurus and Aristonectes (Fig. 4d, e). Comparisons per group show that both basal plesiosaurs and derived polycotylids, the groups with the smallest specimens, produced generally higher D/Vtr. Moreover, we did not find substantial differences between elasmosaurs and thalassophonean pliosauroids (Fig. 4e, Supplementary Table 4; all two-sample t-tests p > 0.05). Both groups had similarly low ranges of D/Vtr regardless of neck length, lower on average than in polycotylids. These results stand even if we exclude Aristonectes, which belongs to the aristonectines, an elasmosaur subfamily with reduced neck length23,45. Further comparisons by morphotype show no significant differences between short-necked pliosauromorphs (here arbitrarily including plesiosaurs with neck ratios below 2×) and long-necked plesiosauromorphs (Supplementary Table 4, all two-sample t-tests p > 0.05). The highest values of D/Vtr occur in animals with trunk lengths of 100 cm or less, followed by a steep decrease between 100–150 cm and a steadier decrease in longer trunks. This indicates a strong negative correlation between trunk dimensions and D/Vtr (Pearson’s product-moment correlation between the log-transformed variables, adjusted r2 = −0.92, p = 2.28 × 10−7). The curve that best describes this relationship is the power equation, D/Vtr = 69.76 × Trunk length−0.944 (Fig. 4e), an almost inversely proportional relationship, consistent with the streamlined nature of these animals for which skin friction drag is dominant.

Polycotylids and thalassophonean pliosaurs, both derived pliosauromorph plesiosaurs9,21, share the same general body proportions9,21,46, but the latter had larger bodies and therefore needed less power in relation to their muscles to move at the same speed. Elasmosaurs on the other hand, despite their disparate morphologies, were no different from thalassophonean pliosaurs in their drag-related costs of forward swimming (Fig. 4c–e) and therefore they were likely to have been equally efficient cruisers.

Earlier research suggested that, even if long necks did not add extra drag during forward swimming, speed in elasmosaurs would have been limited to avoid added drag when their necks bent20. However, when the neck is bent in living forms, the course of swimming changes, as does the flow direction, but the body remains streamlined in the direction of incoming flow. For example, sea lions perform non-powered turns initiated by the head in which the body glides smoothly in a curved position, limiting deceleration47. Further biomechanical research is needed to understand the role of plesiosaur necks in manoeuvrability and other aspects of swimming performance, as well as how these were influenced by shape and flexibility. The well-established idea that long-necked plesiosaurs were sluggish, slow swimmers7,30 is thus not supported here, not because long necks did not increase drag20, but because body size overrode this drag excess.

We analysed trends of body size and neck proportion in a wider sample of sauropterygians, including plesiosaurian and non-plesiosaurian Triassic sauropterygians. Long necks (neck ratio > 3×) occur in taxa with trunk lengths > 150 cm, whereas most sauropterygians had neck ratios of ≤ 2× (Fig. 5a). The great plasticity of body proportions of sauropterygians before and after their transition to a pelagic lifestyle after the Triassic has been well documented21,23,46, but this is the first time that neck and body size have been explored in the context of swimming performance for such a wide sample. We show that overall, sauropterygians and particularly plesiosaurs, mainly explored neck morphologies with little or no effect on drag costs and did not enter morphospaces that were suboptimal for aquatic locomotion (i.e., corresponding to small trunks with long necks; Fig. 5a). In fact, ancestral state reconstruction for trunk length shows that the ancestor of elasmosaurs was likely around 180 cm long and had a relatively short neck with a ratio smaller than 2× (Fig. 5b, c). This indicates that large trunks preceded neck elongation in elasmosaurs and suggests that extreme proportions might have been favoured by a release of hydrodynamic constraints.

a Bivariate plot of the length of trunk and the neck ratio of 79 sauropterygian taxa. Polygons in different colours show area occupied by the main sauropterygian groups. The functional trends describing the effect of each axis are based on results from flow simulations. On the top of this graph, a univariate plot shows the distribution and mean values of trunk length for each group. b, c Phenograms showing the disparity of trunk length (b) and neck ratio (c) in sauropterygians through time. The branches corresponding to basal Plesiosauria (including Rhomaleosauridae and Plesiosauridae), thalassophonean pliosaurs, polycotylids and elasmosaurs are highlighted (colour coding as in a). d, e Sauropterygian trees showing the evolutionary rates for trunk length (d) and neck ratio (e) represented by colour gradient (see Supplementary Fig. 5 for an alternative analysis to 5d using the log10-transformed trunk length). Consensus trees show average results from analyses of 20 cal3-dated trees (see Supplementary Figs. 4 and 6 for analysis on Hedman-dated trees). Rates are based on the mean scalar evolutionary rate parameter.

We next explored evolutionary rates of relative neck length and trunk length in sauropterygians. The pattern of trunk length evolution is consistent with a heterogeneous rates model, not a homogeneous Brownian motion model (log Bayes Factor48 (BF) > 5 in 100% of the sampled trees and > 10 in 92.5%, Supplementary Table 5). Analysis of non-transformed trunk data shows that through the evolution of Sauropterygia, there was a general increase in trunk length with some higher rates, in Triassic nothosauroids, Jurassic rhomaleosaurids and Cretaceous aristonectine elasmosaurs (Fig. 5d; Supplementary Fig. 4a). Additionally, analysis of the log10-transformed trunk data highlights variation in the small-to-medium size ranges and reveals high rates in Triassic eosauropterygians (Supplementary Figs. 5 and 6). The largest trunks evolved independently in two groups, thalassophonean pliosaurids and elasmosaurid plesiosauroids, with no evidence of high rates in the former. In the plesiosauroids, rates are not particularly high in the basal branches, but they are very high in derived aristonectines, and rates for the whole clade were significantly higher than the background rate in 40% of randomisation tests (Supplementary Fig. 7 and Table 6). A progressive increase in body mass over evolutionary time has been described for various clades of aquatic mammals49 and seems to be a common hallmark of the aquatic adaptation to marine pelagic lifestyles in secondarily aquatic tetrapods44. Whether body size reaches a plateau as is the case in cetaceans49 and what constraints influence the evolutionary patterns of size in plesiosaurs remains unexplored. Against this general trend, some derived plesiosaurs, such as polycotylids, saw a reduction in body size, which might have been related to pressures on niche selection, such as adaptation to specific prey, the need for higher manoeuvrability or other ecological factors. As shown earlier, small sizes require lower amounts of total power for a given speed, and therefore would be favoured if for example food resources were limited. This suggests that, in spite of the energy advantages of large size in terms of reduced mass-specific drag29 and metabolic rates49,50, which make it a common adaptation to the pelagic mode of life, other constraints limiting very large sizes were also at work50,51.

A heterogeneous evolutionary rates model for neck proportion is also strongly supported (log BF > 5 in 100% of the sampled trees and > 10 in 45%, Supplementary Table 5). Fast rates are consistently seen at the base of Pistosauroidea (including some Triassic forms and plesiosaurs) and, interestingly, also within elasmosaurs (Fig. 5e; Supplementary Fig. 4b). The neck proportions of elasmosaurs were found to evolve at a faster pace than the background rate in 90% of analyses (randomisation test p-value < 0.001 in 80% and < 0.01 in 10% of the sampled trees; Supplementary Fig. 7 and Table 6). Very fast rates in elasmosaurs are concentrated in the most derived branches (i.e., Euelasmosauridia from the late Upper Cretaceous52) and represent both rapid neck elongation in elasmosaurines and rapid neck shortening in weddellonectians (i.e., aristonectines and closely related taxa52). Additionally, various other independent instances of relative shortening of the neck occurred during the evolution of Sauropterygia, most notably in placodonts, pliosaurs and polycotylids, but these are not associated with high rates.

Our findings contrast with a previous study23 which did not identify any significant evolutionary rate shifts in the neck ratio across Sauropterygia. Here we use a larger number of taxa and a different model fitting approach, which might account for these discrepancies. The association between very long necks and large trunks, along with our flow simulations results and the evidence of high rates in the elongation of necks in elasmosaurines (Fig. 5e), suggests that neck elongation was facilitated by large body sizes. The question remains why neck ratios did not evolve longer than 4×. According to our data, hydrodynamic constraints might have operated against the selection of such long necks. However, it is possible that the primary function for which they were selected, which is still debated30,53, did not require necks with those characteristics. Neck anatomy is likely to be the result of a compromise between different functions/constraints, one of them being hydrodynamic, as shown by the results presented herein.

Six very complete plesiosaur specimens were selected for reconstruction as full-body, three-dimensional models: one basal plesiosaur, (i) the rhomaleosaurid Meyerasaurus victor (specimen exposed in ventral view, SMNS 12478) from the Early Jurassic of Germany; three derived short-necked plesiosaurs, including two thalassophonean pliosaurids, (ii) Peloneustes phylarchus (3D mounted specimen GPIT-RE-3182, previously GPIT 1754/3) from the Middle Jurassic of Germany and (iii) Liopleurodon ferox (3D mounted specimen GPIT-RE-3184, previously GPIT 1754/2) from the Middle Jurassic of the UK, and (iv) the polycotylid Dolychorhynchops osborni (3D mounted specimen KUVP 1300) from the Late Cretaceous of North America; and two long-necked elasmosaurid plesiosaurs, (v) Thalassomedon hanningtoni (3D mounted specimen DMNH 1588) and (vi) Hydrotherosaurus alexandrae (UCMP 33912, figured and reconstructed by Welles54) from the Late Cretaceous of North America. Specimens SMNS 12478, GPIT-RE-3182 and GPIT-RE-3184, belonging to Meyerasaurus, Peloneustes and Liopleurodon, respectively, were examined first-hand and measurements and photographs were taken to inform the modelling. The digital models of Dolichorhynchops55,56, Thalassomedon54,57 and Hydrotherosaurus54 were based on measurements, photographs and two-dimensional reconstructions from the literature (Supplementary Methods). The digital models of three derived ichthyosaurs, previously published in Gutarra et al.29, are based on almost complete specimens: Temnodontosaurus platyodon (NHMUK 2003), a neoichthyosaurian from the Early Jurassic of the U.K.; and two thunnosaurian ichthyosaurs, Stenopterygius quadriscissus (NHMUK R4086) and Ophthalmosaurus icenicus (NHMUK PV R3702, R3898, R4124) from the Early Jurassic of Germany and the Middle–Late Jurassic of the UK, respectively. Additionally, three extant cetaceans were included in this study, the odontocetes Tursiops truncatus (model previously described in Gutarra et al.29) and Orcinus orca, and the mysticete Megaptera novaeangliae. Life reconstructions in lateral and dorsal views, as well as photographs from live specimens, provided the information to model the body and appendages of the bottlenose dolphin Tursiops truncatus58 and the killer whale Orcinus orca58,59. The body and appendages of the humpback whale Megaptera novaeangliae were digitally modelled using an aerial photograph of the dorsal aspect of a wild specimen60, as well as published information on the planform and cross-sectional shape of the flippers61. Digital models for all taxa were created as NURBS geometries using Rhinoceros v. 5 (Supplementary Methods, Supplementary Fig. 8) and are available for download at https://doi.org/10.5281/zenodo.5979631. The geometric parameters of the models (i.e. surface area and volume) were calculated using Rhinoceros measuring tools.

Computer flow simulations were carried out using the commercial software ANSYS-Fluent (v. 18.1 Academic). Our protocol has been validated (i.e., computed drag compared to data from water tank experiments) for external flow over slender bodies29, and proved to replicate with an accuracy ≥ 95% the drag of standard rotational bodies of varying fineness ratios, within a broad range of Reynolds numbers (Re)29. We used the fully-turbulent shear stress transport (SST) model to solve the Reynolds-averaged Navier–Stokes (RANS) equations, as the Reynolds numbers of our analyses fall within the turbulent flow regime (> 106).

Models were imported into ANSYS, where a cylindrical enclosure and refinement box around the wake area were created using the geometry tools. As all taxa are bilaterally symmetrical, only half of the models and the enclosing domain were used in the simulations to economise on computational resources. The virtual flow domain was then meshed with the ANSYS meshing tool, producing grids of 5–15 million elements, depending on the geometry, that combined tetrahedral elements in the region of free flow, and 20–25 layers of prismatic elements in the boundary layer region (i.e., the area adjacent to the non-slip wall surface). Flow was simulated using a double precision, stationary pressure-based solver and a second-order discretization method. Convergence (i.e., the point where the simulation reaches a stable solution) and mesh independence (the influence of mesh size on results) were tested. The total drag coefficient (Cd), as well as the coefficients for its internal components, the viscous drag (Cf) and the pressure drag (Cp) were calculated using the formula:

where D is the drag force in N (total, viscous or pressure drag respectively), ρ is the density of water, 998.2 kg m−3 at 20 °C; u is the inlet velocity in ms−1 and S is the wetted surface area of the model in m2 (see Supplementary Table 2 for a sensitivity test on the flow physical parameters). Because the drag coefficient decreases with increasing Re, the comparison of Cd between taxa was done in conditions of dynamic similarity (i.e., same Re, that is same length and same velocity). The Reynolds numbers used here, 5 × 106–107, covers the range of Re values at which our selected animals likely moved, from about 2 × 106 for Stenopterygius, to 1.6 × 107 for Megaptera, considering a conservative velocity of 1 ms−1. Additionally, the Cd for 5 × 106–107 is numerically close to the average Cd of the much wider range 106–5 × 107 used in a previous study29. The drag results presented here correspond to conditions of zero lift, to eliminate potential variability in the results caused by induced drag. When required, small adjustments were made to the orientation of the models relative to the incoming flow so that the lift remained close to zero.

Our protocol using 3D static CFD simulations provides an objective assessment of the influence of morphology and size on drag forces independent of motion, as well as allowing for comparisons of wide arrays of taxa29. Moreover, aquatic animals commonly use inertial displacement in the absence of movement (i.e. gliding) during submerged swimming to economise energy62,63,64. Computed drag coefficients of dolphins obtained with this method are consistent with estimates obtained from gliding dolphins in water tanks33.

The drag per unit volume represents the contribution of drag to the cost of locomotion (i.e. the energy spent transporting a unit of mass a unit of distance37) in steady swimming (i.e. constant speed, when thrust equals drag force), also called here COTdrag. The total metabolic cost of transport results from dividing the total power (Pin) by the mass (\(m\)) and the velocity (\(u\)):

and therefore, COTdrag can be obtained by dividing the drag power \(({P}_{{{{{{\rm{drag}}}}}}})\) by the mass and the velocity,

where the volume (V), can be considered a proxy for body mass (assuming similar body density). For pelagic swimming animals it is reasonable to assume an approximate body density close to that of sea water65, consistent with measurements from living cetaceans66,67 and estimates of density in extinct marine reptiles57. Note that our calculations consider only the mechanical expenses of locomotion and do not account for the power invested in maintaining the basal metabolism or losses due to muscle efficiency43.

The drag per unit of volume (COTdrag) was calculated for an inflow velocity of 1 ms−1, for models scaled to equal total volume, in order to estimate potential differences in drag-associated costs from body shape alone. Simulations excluding the limbs were added to evaluate the contribution to drag of limbs and body separately. To test the effect of body size on COTdrag, simulations on full models at life-size scale were carried out for an inflow velocity of 1 ms−1. The total body length in the fossil taxa is the average of all available specimens for each genus, with values obtained from personal observations or from the literature (Supplementary Data). The range of sizes for adults of living cetaceans were taken from the literature (Supplementary Data). This study is not concerned with absolute values of the drag power, as it is known that the drag estimated from rigid bodies is smaller than dynamic drag. Hence results are normalised to the values obtained for the bottlenose dolphin Tursiops, here used as a reference.

Comparisons presented herein do not account for propulsive efficiency (\(\eta\)). The reason for this is that there are no clear differences in this parameter between highly specialised caudal oscillation and underwater flying. Estimates from large extant aquatic tetrapods displaying these two swimming styles, such as cetaceans and sea lions, have produced similarly high values of \(\eta\), 0.8 and 0.75–0.9 respectively16. Potential differences in performance due to the kinematics and the shape of propulsive elements should not be dismissed, however, they cannot be included in this model based on current knowledge.

Flow simulations were performed for a set of plesiosaur models with varying neck lengths at a constant Re = 5 × 106 to compute the total drag coefficient (Cd) as well as the coefficients for skin friction (Cf) and pressure drag (Cp). These models were built in Rhinoceros v. 5.0 by enlarging the neck of a basal plesiosaur, the rhomaleosaurid Meyerasaurus victor (Supplementary Fig. 8d), in which the ratio of neck length-to-trunk length (hereafter neck ratio) is 0.87, to encompass neck ratios of 2×, 3×, 4×, 5× and 7×. Meyerasaurus was chosen because of its plesiomorphic characteristics among plesiosaurs and relatively short neck. As shown in the Results, the drag of the limbless bodies does not differ significantly between plesiosaur models (Fig. 1a, Fig. 2a), thus making this model representative of a general plesiosaur morphology. We measured the neck as the distance from the base of the head to the edge of the pectoral girdle and the trunk as the distance between the acetabulum and the glenoid (i.e., inter-girdle distance). Elasmosaurinae, a subfamily of elasmosaurs52, are the plesiosaurs with the longest necks described so far, with neck ratios from 2 to 3.723 (Supplementary Data). Therefore, neck ratios up to 4× correspond to proportions observed in nature, while neck ratios above 4× represent hypothetical body shapes. Previous work suggested that thicker neck contours provide a hydrodynamic advantage in plesiosaurs by reducing the drag coefficient20. To control for the impact of neck thickness on our results, we performed sensitivity tests accounting for this parameter (Supplementary Fig. 9).

Next, the drag at a constant inflow velocity of 1 ms−1 was estimated for the same models of variable neck proportions, this time scaled to a constant trunk length. The total computed drag was then divided by the volume of the trunk. This analysis was aimed to test whether the enlargement of the neck length while maintaining a constant size of trunk would at some point become energetically costly. The trunk contains the muscles involved in locomotion and has previously been considered a better proxy for total body size than total length46 because of the enormous variation in body proportions in sauropterygians.

Finally, to test the interplay of neck ratio and body size a final set of simulations was carried out for various models of plesiosaurs at life-size dimensions, including seven extra limbless plesiosaur models to ensure a better representation of trunk sizes and neck proportions in derived plesiosaurs. We used photographs or reconstructions from the literature and where possible, material obtained from personal observation, of well-preserved adult specimens of the basal plesiosaurs Rhomaleosaurus thorntoni (reconstruction by Smith & Benson68) and Plesiosaurus sp. (3D mounted skeleton in the National Museum of Wales, personal observation); the pliosauroid Kronosaurus boyacensis (MJACM169); the polycotylid Polycotylus latippinus (reconstruction by O’Keefe and Chiappe70); the elasmosaurine elasmosaurs Styxosaurus sp. (SDSM 45171) and Albertonectes vanderveldei44; and finally the short-necked aristonectine elasmosaur Aristonectes quiriquinensis (SGO.PV.95745).

The drag force was computed from simulations at the same inlet velocity of 1 ms−1, then divided by the volume of the trunk for each model and finally normalised to the results of the Meyerasaurus model. The relative values of drag-per-unit-trunk-volume (D/Vtr) were visualised with a colour gradient over a two-dimensional plot of neck ratio and trunk length. Contour lines showing interpolated values for the rest of the plot area were added using the function geom_contour of the package ggplot2 v.3.3.272 in R v.3.6.273. Additionally, the drag-per-unit-trunk-volume was plotted against the trunk length, and the correlation between the log-transformed variables was calculated using Pearson’s product moment correlation coefficient. A univariate linear regression model was fitted using the R package moonBook74. The resulting linear equation \({{\log }}(y)=a{{\log}}\left(x\right)+b\), was finally transformed into a power equation to represent the relationship between the non-transformed variables, \(y={e}^{b}+{x}^{n}\).

We compiled a dataset of neck and trunk lengths for 79 sauropterygian species through the Mesozoic (Supplementary Data). Neck ratios and trunk lengths were visualised in a bivariate morphospace showing areas occupied by different sauropterygian clades. Functional explanations of the two axes on drag-per-unit-volume are based on previous CFD analyses on real and hypothetical models. The trunk length (cm) and neck proportions were plotted on a phylogeny of sauropterygians and values were inferred at ancestral nodes using the R package phytools v0.7-4775. For this, an informal composite sauropterygian tree was assembled based on a published phylogeny6 that combines phylogenetic relationships of Plesiosauria76 and Triassic sauropterygians77 (Supplementary Methods).

Rates of evolution were estimated for trunk length (using raw and log10-transformed data) and neck ratio on 40 time-calibrated phylogenies in a Bayesian framework, using the variable-rates model of BayesTraits v.2.0.278 and R v.3.6.2, with code from Stubbs et al.79. The sauropterygian tree was time-calibrated 20 times, using both the cal380,81 and Hedman82 scaling methods to account for uncertainties of fossil occurrences and differences in dating models (Supplementary Methods). Evolutionary rate heterogeneity was evaluated for the 20 trees from each dating approach with a reversible jump Markov Chain Monte Carlo algorithm (rjMCMC) using default prior distributions. The analysis used 2 billion iterations, of which the first 400 million were discarded as burn-in, and parameters were sampled every 80,000 iterations. The method detects shifts in rates of evolution by incorporating branch-specific scalars and rescales branch lengths that deviate from expectations of a homogeneous Brownian motion (BM) model. Rates results were summarised as consensus trees derived from the 20 iterations of each dating approach using phytools75, in which colours of branches represent mean rate scalars from all 20 trees. Convergence was tested using the minimum effective sample size function of the R package CODA v.0.19.383. The fit of a heterogeneous variable-rates model was tested against a null homogeneous random walk model (BM) using log Bayes factors48 (BF), calculated from the marginal likelihoods of these two models, obtained using the stepping-stone sampling method, with 100 stones per run for 1000 iterations84. Lastly, a randomisation test was applied to test for significantly different rates of evolution in Elasmosauridae, Polycotylidae and Thalassophonea. This analysis performs multiple random samplings of mean rate values and calculates differences between sampled and then non-sampled branches over 9999 replicates85.

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Supplementary Information includes Supplementary Figures, Supplementary Tables and Supplementary Methods; Supplementary Data includes calculations supporting the present results and datasets used in the evolutionary rates analyses. These files, as well as the digital models created for the computer flow simulations performed in this study, can be found in the GitHub repository (https://github.com/SusanaGutarra/Plesiosaur-hydrodynamics-evolution) and at https://doi.org/10.5281/zenodo.597963186.

The authors declare that the code supporting the findings of this study is available at https://doi.org/10.5281/zenodo.597963186.

A Correction to this paper has been published: https://doi.org/10.1038/s42003-022-03458-x

Kelley, N. P. & Pyenson, N. D. Evolutionary innovation and ecology in marine tetrapods from the Triassic to the Anthropocene. Science 348, aaa3716 (2015).

Article PubMed Google Scholar

Vermeij, G. J. & Motani, R. Land to sea transitions in vertebrates: the dynamics of colonization. Paleobiology 44, 237–250 (2018).

Article Google Scholar

Pyenson, N. D., Kelley, N. P. & Parham, J. F. Marine tetrapod macroevolution: physical and biological drivers on 250 Ma of invasions and evolution in ocean ecosystems. Palaeogeogr. Palaeoclimatol. Palaeoecol. 400, 1–8 (2014).

Article Google Scholar

Motani, R. & Vermeij, G. J. Ecophysiological steps of marine adaptation in extant and extinct non‐avian tetrapods. Biol. Rev. 96, 1769–1798 (2021).

Article PubMed Google Scholar

Motani, R. The evolution of marine reptiles. Evolution: Educ. Outreach 2, 224–235 (2009).

Google Scholar

Stubbs, T. L. & Benton, M. J. Ecomorphological diversifications of Mesozoic marine reptiles: the roles of ecological opportunity and extinction. Paleobiology 42, 547–573 (2016).

Article Google Scholar

Massare, J. A. Swimming capabilities of Mesozoic marine reptiles: implications for method of predation. Paleobiology 14, 187–205 (1988).

Article Google Scholar

Bernard, A. et al. Regulation of body temperature by some Mesozoic marine reptiles. Science 328, 1379–1382 (2010).

Article CAS PubMed Google Scholar

Benson, R. B. J., Evans, M. & Druckenmiller, P. S. High diversity, low disparity and small body size in plesiosaurs (Reptilia, Sauropterygia) from the Triassic–Jurassic boundary. PLoS ONE 7, e31838 (2012).

Article CAS PubMed PubMed Central Google Scholar

Fleischle, C. V., Sander, P. M., Wintrich, T. & Caspar, K. R. Hematological convergence between Mesozoic marine reptiles (Sauropterygia) and extant aquatic amniotes elucidates diving adaptations in plesiosaurs. PeerJ 7, e8022 (2019).

Article PubMed PubMed Central Google Scholar

Motani, R. Evolution of fish-shaped reptiles (Reptilia: Ichthyopterygia) in their physical environments and constraints. Annu. Rev. Earth Planet. Sci. 33, 395–420 (2005).

Article CAS Google Scholar

Motani, R., You, H. & McGowan, C. Eel-like swimming in the earliest ichthyosaurs. Nature 382, 347–348 (1996).

Article Google Scholar

O’Keefe, F. R. Ecomorphology of plesiosaur flipper geometry. J. Evolut. Biol. 14, 987–991 (2001).

Article Google Scholar

Muscutt, L. E. et al. The four-flipper swimming method of plesiosaurs enabled efficient and effective locomotion. Proc. R. Soc. B: Biol. Sci. 284, 20170951 (2017).

Article Google Scholar

Fish, F. E. & Rohr, J. J. Review of dolphin hydrodynamics and swimming performance. (1999).

Fish, F. E. Biomechanics and energetics in aquatic and semiaquatic mammals: platypus to whale. Physiol. Biochem. Zool. 73, 683–698 (2000).

Article CAS PubMed Google Scholar

Motani, R. Scaling effects in caudal fin propulsion and the speed of ichthyosaurs. Nature 415, 309–312 (2002).

Article PubMed Google Scholar

Motani, R. et al. Absence of suction feeding ichthyosaurs and its implications for Triassic mesopelagic paleoecology. PLoS ONE 8, e66075 (2013).

Article PubMed PubMed Central Google Scholar

Liu, S. et al. Computer simulations imply forelimb-dominated underwater flight in plesiosaurs. PLOS Computational Biol. 11, e1004605 (2015).

Article Google Scholar

Tröelsen, P. V., Wilkinson, D. M., Seddighi, M., Allanson, D. R. & Falkingham, P. L. Functional morphology and hydrodynamics of plesiosaur necks: does size matter? J. Vertebrate Paleontol. 39, e1594850 (2019).

Article Google Scholar

O’Keefe, F. R. The evolution of plesiosaur and pliosaur morphotypes in the Plesiosauria (Reptilia: Sauropterygia). Paleobiology 28, 101–112 (2002).

2.0.CO;2" data-track-action="article reference" href="https://doi.org/10.1666%2F0094-8373%282002%29028%3C0101%3ATEOPAP%3E2.0.CO%3B2" aria-label="Article reference 21" data-doi="10.1666/0094-8373(2002)0282.0.CO;2">Article Google Scholar

Otero, R. A. Taxonomic reassessment of Hydralmosaurus as Styxosaurus: new insights on the elasmosaurid neck evolution throughout the Cretaceous. PeerJ 4, e1777 (2016).

Article PubMed PubMed Central Google Scholar

Soul, L. C. & Benson, R. B. J. Developmental mechanisms of macroevolutionary change in the tetrapod axis: a case study of Sauropterygia. Evolution 71, 1164–1177 (2017).

Article PubMed PubMed Central Google Scholar

Motani, R. Swimming speed estimation of extinct marine reptiles-energetic approach revisited. Paleobiology 28, 251–262 (2002).

2.0.CO;2" data-track-action="article reference" href="https://doi.org/10.1666%2F0094-8373%282002%29028%3C0251%3ASSEOEM%3E2.0.CO%3B2" aria-label="Article reference 24" data-doi="10.1666/0094-8373(2002)0282.0.CO;2">Article Google Scholar

Massare, J. A. Swimming capabilities of Mesozoic marine reptiles: a review. in Mechanics and Physiology of Animal Swimming (eds. Maddock, L., Bone, Q. & Rayner, J. M. V.) 133–149 (Cambridge University Press, 1994).

Young, B. A. The calculation of the total and skin friction drags of bodies of revolution at zero incidence. Aeronautical Res. Comm. 1874, 51–70 (1939).

Google Scholar

Von Mises, R. Theory of Flight. (McGraw-Hill, 1945).

Fish, F. E. Imaginative solutions by marine organisms for drag reduction. In Proceedings of the International Symposium on Seawater Drag Reduction 443–450 (1998).

Gutarra, S. et al. Effects of body plan evolution on the hydrodynamic drag and energy requirements of swimming in ichthyosaurs. Proc. R. Soc. B: Biol. Sci. 286, 20182786 (2019).

Article Google Scholar

Noe, L., Taylor, M. & Gómez-Pérez, M. An integrated approach to understanding the role of the long neck in plesiosaurs. Acta Palaeontol. Polonica 62, 137–162 (2017).

Article Google Scholar

Hoerner, S. F. Fluid-Dynamic Drag: Practical Information on Aerodynamic Drag and Hydrodynamic Resistance. (Published by S. F. Hoerner, 1965).

Schlichting, H. & Gersten, K. Boundary-Layer Theory. (Springer, 2017).

Lang, T. G. & Daybell, D. A. Porpoise performance tests in a sea-water tank. Naval Ordnance Test Station. California Institute of Technology (1963).

Riedeberger, D. & Rist, U. Numerical simulation of laminar-turbulent transition on a dolphin using the γ-Re θ model. in High Performance Computing in Science and Engineering (eds. Nagel, W. E., Kröner, D. B. & Resch, M. M.) 379–391 (Springer, 2012).

Gutarra, S. & Rahman, I. A. The locomotion of extinct secondarily aquatic tetrapods. Biol. Rev. 97, 67–98 (2021).

Article PubMed Google Scholar

Aleyev, Y. G. Reducing Resistance to Movement. in Nekton (ed. Aleyev, Y. G.) 185–283 (Springer Netherlands, 1977).

Schmidt-Nielsen, K. Locomotion: energy cost of swimming, flying, and running. Science 177, 222–228 (1972).

Article CAS PubMed Google Scholar

Fish, F. E. Diversity, mechanics and performance of natural aquatic propulsors. In WIT Transactions on State of the Art in Science and Engineering (ed. Liebe, R.) 57–87 (WIT Press, 2006).

Smith, A. S. Reassessment of ‘Plesiosaurus’ megacephalus (Sauropterygia: Plesiosauria) from the Triassic-Jurassic boundary, UK. Palaeontologia Electron. 18, 1–20 (2015).

CAS Google Scholar

Sato, K. et al. Stroke frequency, but not swimming speed, is related to body size in free-ranging seabirds, pinnipeds and cetaceans. Proc. R. Soc. B: Biol. Sci. 274, 471–477 (2007).

Article Google Scholar

Williams, T. M. et al. Swimming and diving energetics in dolphins: a stroke-by-stroke analysis for predicting the cost of flight responses in wild odontocetes. J. Exp. Biol. 220, 1135–1145 (2017).

Article PubMed Google Scholar

Gough, W. T. et al. Scaling of swimming performance in baleen whales. J. Exp. Biol. 222, jeb204172 (2019).

Article PubMed Google Scholar

Williams, T. M. The evolution of cost efficient swimming in marine mammals: limits to energetic optimization. Philos. Trans. R. Soc. Lond. B: Biol. Sci. 354, 193–201 (1999).

Article Google Scholar

Kubo, T., Mitchell, M. T. & Henderson, D. M. Albertonectes vanderveldei, a new elasmosaur (Reptilia, Sauropterygia) from the Upper Cretaceous of Alberta. J. Vertebrate Paleontol. 32, 557–572 (2012).

Article Google Scholar

Otero, R. A., Soto-Acuña, S. & O’keefe, F. R. Osteology of Aristonectes quiriquinensis (Elasmosauridae, Aristonectinae) from the upper Maastrichtian of central Chile. J. Vertebrate Paleontol. 38, e1408638 (2018).

Article Google Scholar

O’Keefe, F. R. & Carrano, M. T. Correlated trends in the evolution of the plesiosaur locomotor system. Paleobiology 31, 656–675 (2005).

Article Google Scholar

Fish, F. E. Maneuverability by the sea lion Zalophus californianus: turning performance of an unstable body design. J. Exp. Biol. 206, 667–674 (2003).

Article PubMed Google Scholar

Kass, R. E. & Raftery, A. E. Bayes factors. J. Am. Stat. Assoc. 90, 773–795 (1995).

Gearty, W., McClain, C. R. & Payne, J. L. Energetic tradeoffs control the size distribution of aquatic mammals. Proc. Natl Acad. Sci. 115, 4194–4199 (2018).

Article CAS PubMed PubMed Central Google Scholar

Goldbogen, J. A. Physiological constraints on marine mammal body size. Proc. Natl Acad. Sci. USA 115, 3995–3997 (2018).

Article CAS PubMed PubMed Central Google Scholar

Gough, W. T. et al. Scaling of oscillatory kinematics and Froude efficiency in baleen whales. J. Exp. Biol. 224, jeb237586 (2021).

Article PubMed PubMed Central Google Scholar

O’Gorman, J. P. Elasmosaurid phylogeny and paleobiogeography, with a reappraisal of Aphrosaurus furlongi from the Maastrichtian of the Moreno Formation. J. Vertebrate Paleontol. 39, e1692025 (2019).

Article Google Scholar

Nagesan, R. S., Henderson, D. M. & Anderson, J. S. A method for deducing neck mobility in plesiosaurs, using the exceptionally preserved Nichollssaura borealis. R. Soc. Open Sci. 5, 172307 (2018).

Article PubMed PubMed Central Google Scholar

Welles, S. P. Elasmosaurid plesiosaurs with description of new material from California and Colorado. Mem. Univ. Calif. 13, 125–254 (1943).

Google Scholar

Bonner, O. W. An osteological study of Nyctosaurus and Trinacromerum with a description of a new species of Nyctosaurus. (Fort Hays Kansas State College, 1964).

Carpenter, K., Sanders, F., Reed, B., Reed, J. & Larson, P. Plesiosaur swimming as interpreted from skeletal analysis and experimental results. Trans. Kans. Acad. Sci. 113, 1–34 (2010).

Article Google Scholar

Henderson, D. Floating point: a computational study of buoyancy, equilibrium, and gastroliths in plesiosaurs. Lethaia 39, 227–244 (2006).

Article Google Scholar

Jefferson, T. A., Webber, M. A. & Pitman, R. Marine Mammals of the World. (Academic Press, 2015).

Heyning, J. E. & Dahlheim, M. E. Orcinus orca. Mamm. Species 304, 1–9 (1988).

Article Google Scholar

Christiansen, F., Dujon, A. M., Sprogis, K. R., Arnould, J. P. Y. & Bejder, L. Noninvasive unmanned aerial vehicle provides estimates of the energetic cost of reproduction in humpback whales. Ecosphere 7, e01468 (2016).

Article Google Scholar

Fish, F. E. & Battle, J. M. Hydrodynamic design of the humpback whale flipper. J. Morphol. 225, 51–60 (1995).

Article CAS PubMed Google Scholar

Weihs, D. Energetic advantages of burst swimming of fish. J. Theor. Biol. 48, 215–229 (1974).

Article CAS PubMed Google Scholar

Williams, T. M. et al. Sink or swim: strategies for cost-efficient diving by marine mammals. Science 288, 5 (2000).

Gleiss, A. C. et al. Convergent evolution in locomotory patterns of flying and swimming animals. Nat. Commun. 2, 352 (2011).

Article PubMed Google Scholar

Vogel, S. Life in Moving Fluids: The Physical Biology of Flow. (Princeton University Press, 1994).

Miller, P. et al. Body density and diving gas volume of the northern bottlenose whale (Hyperoodon ampullatus). J. Exp. Biol. 219, 2458–2468 (2016).

Article PubMed PubMed Central Google Scholar

Narazaki, T. et al. Body density of humpback whales (Megaptera novaeangliae) in feeding aggregations estimated from hydrodynamic gliding performance. PLoS ONE 13, e0200287 (2018).

Article PubMed PubMed Central Google Scholar

Smith, A. S. & Benson, R. B. J. Osteology of Rhomaleosaurus thorntoni (Sauropterygia: Rhomaleosauridae) from the Lower Jurassic (Toarcian) of Northamptonshire, England. Monogr. Palaeontographical Soc. 168, 1–40 (2014).

Article Google Scholar

Hampe, O. Ein großwüchsiger pliosauride (Reptilia: Plesiosauria) aus der Unterkreide (oberes Aptium) von Kolumbien. Cour. Forschungsinstitut Senckenberg 145, 1–32 (1992).

Google Scholar

O’Keefe, F. R. & Chiappe, L. M. Viviparity and K-selected life history in a Mesozoic marine plesiosaur (Reptilia, Sauropterygia). Science 333, 870–873 (2011).

Article PubMed Google Scholar

Welles, S. P. A review of the North American Cretaceous elasmosaurs. Univ. Calif. Publ. Geol. Sci. 29, 47–144 (1952).

Google Scholar

Wickham, H. ggplot2: Elegant Graphics for Data Analysis. (Springer International Publishing, 2016). https://doi.org/10.1007/978-3-319-24277-4.

R Core Team. R: a language and environment for statistical computing (R Foundation for Statistical Computing). (R Core Team, 2019).

Moon, K.-W. R statistics and graphs for medical papers. (Hannarae, 2015).

Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things): phytools: R package. Methods Ecol. Evolution 3, 217–223 (2012).

Article Google Scholar

Benson, R. B. J. & Druckenmiller, P. S. Faunal turnover of marine tetrapods during the Jurassic-Cretaceous transition: faunal turnover of marine. Biol. Rev. 89, 1–23 (2014).

Article PubMed Google Scholar

Cheng, Y.-N., Wu, X.-C., Sato, T. & Shan, H.-Y. A new eosauropterygian (Diapsida, Sauropterygia) from the Triassic of China. J. Vertebrate Paleontol. 32, 1335–1349 (2012).

Article Google Scholar

Baker, J., Meade, A., Pagel, M. & Venditti, C. Positive phenotypic selection inferred from phylogenies. Biol. J. Linn. Soc. 118, 95–115 (2016).

Article Google Scholar

Stubbs, T. L. et al. Ecological opportunity and the rise and fall of crocodylomorph evolutionary innovation. Proc. R. Soc. B. 288, 20210069 (2021).

Article PubMed PubMed Central Google Scholar

Bapst, D. W. paleotree: an R package for paleontological and phylogenetic analyses of evolution: Analyses of Paleo-Trees in R. Methods Ecol. Evolution 3, 803–807 (2012).

Article Google Scholar

Bapst, D. W. A stochastic rate-calibrated method for time-scaling phylogenies of fossil taxa. Methods Ecol. Evolution 4, 724–733 (2013).

Article Google Scholar

Hedman, M. M. Constraints on clade ages from fossil outgroups. Paleobiology 36, 16–31 (2010).

Article Google Scholar

Plummer, M., Best, N., Cowles, K. & Vines, K. CODA: convergence diagnosis and output analysis for MCMC. R. N. 6, 7–11 (2006).

Google Scholar

Xie, W., Lewis, P. O., Fan, Y., Kuo, L. & Chen, M.-H. Improving marginal likelihood estimation for Bayesian phylogenetic model selection. Syst. Biol. 60, 150–160 (2011).

Article PubMed Google Scholar

Ballell, A., Moon, B. C., Porro, L. B., Benton, M. J. & Rayfield, E. J. Convergence and functional evolution of longirostry in crocodylomorphs. Palaeontology 62, 867–887 (2019).

Article Google Scholar

Gutarra, S., Stubbs, T. L., Moon, B. C., Palmer, C. & Benton, M. J. Large size in aquatic tetrapods compensates for high drag caused by extreme body proportions (data and code). https://github.com/SusanaGutarra/Plesiosaur-hydrodynamics-evolution and https://doi.org/10.5281/zenodo.5979631.

Download references

We thank Imran Rahman (Natural History Museum, London), Stephan Lautenschlager (University of Birmingham) and Emily Rayfield (University of Bristol) for comments and discussions on early versions of this manuscript. We are grateful to Armin Elsler (University of Bristol) for sharing custom R-code to extract and visualise the results of the BayesTraits evolutionary rates analysis. We thank Ingmar Werneburg (Palaeontology Collection, Tübingen, Germany), Erin Maxwell (Staatliches Museum für Naturkunde Stuttgart, Germany) and Caroline Buttler (National Museum of Wales, UK) for facilitating access to fossil marine reptile specimens. We would also like to thank Neil Kelly, and an anonymous reviewer for constructive comments that helped improve this manuscript. S.G. was supported by the Natural Environment Research Council and the GW4+DTP (grant no. NE/L002434/1), with additional support from CASE partner National Museum of Wales; B.C.M. and T.L.S. are funded by NERC grant NE/P013724/1 and M.J.B by ERC grant 788203 (INNOVATION).

School of Earth Sciences, University of Bristol, Life Sciences Building, 24 Tyndall Avenue, Bristol, BS8 1TQ, UK

Susana Gutarra, Thomas L. Stubbs, Benjamin C. Moon, Colin Palmer & Michael J. Benton

Natural History Museum, Cromwell Road, London, SW7 5BD, UK

Susana Gutarra

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

S.G. designed and executed the CFD simulations, compiled the phenotypic data, assembled the sauropterygian tree and wrote the main manuscript. C.P. helped design and supervised simulation analyses; M.J.B, B.C.M. and T.L.S. helped design and supervised the evolutionary rates analyses. T.L.S performed the dating and rates analyses. All authors contributed to the final manuscript.

Correspondence to Susana Gutarra.

The authors declare no competing interests.

Communications Biology thanks Neil Kelley and the other, anonymous, reviewer for their contribution to the peer review of this work. Primary Handling Editor: Luke R. Grinham. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Gutarra, S., Stubbs, T.L., Moon, B.C. et al. Large size in aquatic tetrapods compensates for high drag caused by extreme body proportions. Commun Biol 5, 380 (2022). https://doi.org/10.1038/s42003-022-03322-y

Download citation

Received: 08 September 2021

Accepted: 25 March 2022

Published: 28 April 2022

DOI: https://doi.org/10.1038/s42003-022-03322-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

BMC Ecology and Evolution (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.