In the present study, we show that the fastest runners and swimmers are becoming not only faster but also heavier, taller and more slender. During the past century, the world record speeds for 100 m-freestyle and 100 m-dash have increased with body mass (M) raised to the power 1/6, in accordance with the constructal scaling of animal locomotion. The world records also show that the speeds have increased in proportion with body heights (H)raised to the power 1/2, in accordance with animal locomotion scaling. If the athlete's body is modeled with two length scales (H, body width L), the (M, H) data can be used to calculate the slenderness of the body, H/L. The world records show that the body slenderness is increasing very slowly over time.

Animal locomotion scaling

Broadly speaking, larger animals move faster on earth than smaller animals(Bejan, 2000; Hoppeler and Weibel, 2005; Weibel, 2000). This aspect of animal scaling is often overlooked because large animals appear sluggish. Geese flap their wings infrequently. Whales swing their tails as if not in a hurry. The scaling relations of animal locomotion have been measured and studied empirically as three separate locomotion mechanisms: flying(Bartholomew and Casey, 1978; Greenewalt, 1975; Lighthill, 1974; Marden et al., 1997; May, 1995; Tennekes, 1997; Wakeling, 1997), running(Heglund et al., 1974; Iriarte-Diaz, 2002; Marsh, 1988; Pennycuick, 1975) and swimming(Arnott et al., 1998; Brett, 1965; Childress and Dudley, 2004; Drucker and Jensen, 1996; Kiceniuk and Jones, 1977; Peake and Farrell, 2004; Rohr and Fish, 2004; Videler, 1993).

More recently, a broader view of the commonality of animal locomotion has been emerging (Ahlborn, 2004; Bejan, 2000; Bejan, 2005; Bejan and Marden, 2006; Bejan and Lorente, 2008; Hoppeler and Weibel, 2005; Marden and Allen, 2002; Muller and van Leeuwen, 2004; Taylor et al., 2003; Weibel, 2000). For example,according to constructal theory, animal locomotion is a rhythm of body motion constructed such that the animal achieves a balance between two expenditures of useful energy: lifting weight on the vertical, and overcoming drag while progressing on the horizontal (this analysis is reviewed in the Appendix). The sum of the two efforts is minimal when the two efforts are of the same size. In this regime the body-mass scaling relations (slope and intercept) of animal locomotion are predicted and are in agreement with the measurements of swimmers, runners and fliers over the body mass (M) range 10–6–103 kg:

  • travel speeds (V) proportional to M1/6, where M is the body mass, ρ is the body density and gis gravitational acceleration, for example, for running on soft ground and swimming:
    \[\ V{\sim}\mathbf{g}^{1{/}2}{\rho}^{-1{/}6}M^{1{/}6},\]
  • body frequencies (t–1) (flapping, stride, fish tailing) proportional to M–1/6:
    \[\ t^{-1}{\sim}\mathbf{g}^{1{/}2}{\rho}^{1{/}6}M^{-1{/}6},\]
  • body force (F) scale equal to M:
    \[\ F{\sim}M\mathbf{g},\]
  • food requirement (useful energy, W) per distance traveled(Lx), which is proportional to M:
    \[\ \frac{W}{L_{\mathrm{x}}}{\sim}M\mathbf{g}.\]

The corresponding scaling laws for flying and running with air drag are similar to Eqns 1, 2, 3, 4. These relations are accurate within a dimensionless factor of order 1, as they were derived based on scale analysis. In spite of this built-in approximation, they agree well with the large body of experimental data available(Bejan and Marden, 2006).

Speed and body mass

We used this constructal framework to examine the evolution of speeds in modern athletics: the evolution of the sport (winning speeds and body metrics), not the evolution of the athletes. We focused on the two most documented probes for men, the 100 m-freestyle in swimming and the 100 m-dash in track. These are sprint probes, not endurance events. Sprint probes require intense expenditure of work during a relatively short period of time.

In Fig. 1A and Table 1 we see the evolution of the world speed record (V) for male 100 m-freestyle swimming since 1912. Because of the theoretical scaling (i), we also researched the evolution of the body masses of the record-breaking athletes(Fig. 1B). Both V and M have been increasing in time (t). By eliminating t between Fig. 1A,B,we found Fig. 1C, which shows the evolution of V vs M.

There is scatter in Fig. 1Cbecause the span of the V and M data is short, much shorter than the span of all biological cases of animal locomotion correlated in Bejan and Marden (Bejan and Marden,2006), where the M range was 10–6–103 kg. The shortness of the contemporary timeframe has the effect of magnifying the scatter of the data. Several additional factors also contribute to the scatter, for example,technology (space age swimming suits, running shoes and chronometry),competition environment (state of the art aquatic and track and field venues)and changes in the rules of competition. In spite of these random variables,the evolutionary direction of animal locomotion(Eqn 1) is respected: as an average, the faster swimmers are bigger. The best fit of the (V, M)data of Fig. 1C according to the theoretical proportionality between V and M1/6 (cf. Eqn 1) is:
\[\ V{\sim}0.72M^{1{/}6},\]
where V and M are expressed in m s–1 and kg, respectively. Eqn 5a was obtained by power law regression, with R2=0.171. The P-value is 0.028, and because it is less than 0.05, the correlation shown in Eqn 5a is statistically significant (Soong, 2004; Vogt, 2005).
The P-value is even smaller if we correlate the VM data of Fig. 1C with a more general power law V=aMb, in which the constants a and b can be optimized. The best fit of this kind is:
\[\ V{\sim}0.68M^{0.23},\]
for which the P-value is 0.023 and R2=0.19. Note the slight difference between the exponent 0.23 and the theoretical exponent 1/60.17. The scaling (Eqn 5b)represents a steeper increase in V with M than in the broad-range animal scaling (Eqn 1). Numerically, the two formulas(Eqn 5b and Eqn 1) agree in the Mrange of humans, and for this reason the VM1/6 scaling of Eqn 5a is sufficient for concluding that the animal scaling (Eqn 1) manifests itself in the evolution of swimming speeds vs mass among record holders.
For 100 m dash (Fig. 2), the trend, scatter and conclusion are the same. Eliminating t between Fig. 2A,B, we arrived at Fig. 2C. The power law regression equation for the data in Fig. 2C is:
\[\ V{\sim}4.85{\ }M^{1{/}6},\]
with R2=0.364. The P-value is 0.007, again indicating statistical significance.

Noteworthy are the factors 0.72 and 4.85 in Eqns 5a and 6, respectively. These factors of order 1 agree with the theoretical scaling(Eqn 1), which after substituting the values for g and ρ yields for swimming and running V∼1×M1/6. Here `1' is the intercept of the line plotted as log V vs log M. Also noteworthy is that the factor for running (4.85) is greater than the factor for swimming (0.72). This also agrees with the manner in which the empirical factor (not shown in Eqn 1 but reported in the Appendix) differentiates between the power-law correlations of animal speed data for runners and swimmers (Bejan and Marden, 2006).

Body height

The conclusion that body size has an effect on speed, Eqns 5a and 6, agrees fully with the doctrine of animal scaling (Bejan,2000; Hoppeler and Weibel,2005; Weibel,2000). Size can be expressed not only as M but also as body height (H). In the scale analysis that led to Eqns 1, 2, 3, 4, the body was modeled in the simplest possible way: with one length scale, which meant that M∼ρLb3, where Lb is the lone length scale. Accordingly, the theoretical speed V of Eqn 1 should be proportional to Lb1/2. In Fig. 3A and Fig. 4A we plotted the data of Tables 1 and 2 by using the H of the athletes as the length scale Lb. We determined the best correlations of type VLb1/2,based on power law regression:
\[\ V{\sim}1.37{\ }L_{\mathrm{b}}^{1{/}2}(100{\ }\mathrm{m}{\ }\mathrm{freestyle}),\]
\[\ V{\sim}7.45{\ }L_{\mathrm{b}}^{1{/}2}(100\mathrm{m}{\ }\mathrm{dash}),\]
where Lb is expressed in meters, R2=0.248 and 0.433 for Eqns 7 and 8, respectively. The corresponding P-values are 0.009 and 0.001, respectively; thus,indicating statistical significance for both correlations. The H data of Tables 1 and 2 are plotted vs t in Fig. 3B and Fig. 4B, respectively.

The proportionality between speeds and body length raised to the power 1/2(Eqns 7 and 8) suggests a simpler way to derive the speed–mass scaling rule (Eqn 1), much simpler than the analysis shown in the Appendix. During each cycle of locomotion the body falls from a height of order Lb. The time scale of the fall is of order t∼(Lb/g)1/2. The body falls forward to a distance of order Lb; therefore,the horizontal velocity scale is VLb/t∼(gLb)1/2. Combining this V scale with Lb∼(M/ρ)1/3 we arrive at Eqn 1.

Body slenderness

A body model that is more realistic than the single-scale model is a cylinder of height H and diameter (width) L. The Min this model is M=ρ(π/4)L2H. We used the body density (ρ∼1000 kg m–3), the recorded M and the recorded H to calculate the athlete's width scale L=(4M/ρΠH)1/2. We then used the recorded H and the calculated L to define the slenderness(S) of the body:
\[\ S=\frac{H}{L}=\left(\frac{{\pi}{\rho}H^{3}}{4{\ }M}\right)^{1{/}2}.\]
The S values calculated in this manner are reported for each athlete in Tables 1 and 2. Their significance becomes apparent if we combine this two-scale body model with the locomotion model proposed at the end of the preceding section.
In swimming, the vertical length scale is L, the time scale of the fall is (L/g)1/2 and the forward speed is of order V∼(gL)1/2. Omitting factors of order 1 (such as π/4), we combine the mass scale(M∼ρL2H) with S=H/L and obtain L∼(M/ρ)1/3S–1/3and:
\[\ V_{\mathrm{sim}}{\sim}\mathbf{g}^{1{/}2}{\rho}^{-1{/}6}M^{1{/}6}S^{-1{/}6}.\]
In running, the vertical length scale is H, and the corresponding scales are t∼(H/g)1/2, V∼(gH)1/2 and H∼(M/ρ)1/3S2/3. The speed–mass relation that replaces Eqn 1 is:
\[\ V_{\mathrm{run}}{\sim}\mathbf{g}^{1{/}2}{\rho}^{-1{/}6}M^{1{/}6}S^{1{/}3}.\]

The S effect differentiates between running and swimming. Dividing Eqns 11 and 10 we anticipate Vrun/VswimS1/2,which is a number of the same order as the ratio between Eqns 6 and 5a. The two-scale model also suggests that from among athletes with the same mass, the ones with larger S values are more likely to run fast. In swimming, the Seffect is the opposite but weaker: swimmers would be slightly faster if more robust (smaller S).

Fig. 5A,B show the S data plotted vs t, and indicate a weak progress toward larger S values for both swimming and running. The best linear fits for the two sets of data are:
\[\ S=0.48+0.0038t(100{\ }\mathrm{m}{\ }\mathrm{freestyle}),\]
\[\ S=1.7+0.0031t(100{\ }\mathrm{m}{\ }\mathrm{dash}),q\]
where t is the year, R2=0.31 and 0.13, and P=0.001 and 0.07, respectively.

The same S data are plotted against M in Fig. 6A,B. The data are too sparse to yield statistically significant correlations; however, qualitatively they suggest a slight increase in S vs M for running and a slight decrease in S vs M for swimming.

The scaling trends revealed by the speed data suggest that speed records will continue to be dominated by heavier and taller athletes. This trend is due to the scaling rules of animal locomotion, not to the contemporary increase in the average body size of humans. The mean height of humans has increased by roughly 5 cm from 1900 to 2002(Plastic Soldier Review,2002). During the same century, the mean height of champion swimmers and runners has increased by 11.4 cm and 16.2 cm, respectively(Fig. 3C, Fig. 4C).

The insight gained in this paper allows us to speculate what the running speeds might have been in ancient Greece and the Roman Empire. There is no record of what the winning speeds were then, because the competition was for winning the race, not for breaking a time record. Chronometry did not exist. In antiquity body masses were roughly 70% of what they are today(Plastic Soldier Review, 2002; Hpathy 2009; National Health and Nutrition Examination Survey, 1999). According to Eqn 6, this means that speeds were lower by a factor of roughly(0.7)1/6=0.94. In other words, if the 100 m dash in military training today is won in 13 s, 2000 years ago it would have been won in∼14 s.

This insight also teaches us why certain training techniques are successful in high-performance sports. For example, in modern speed swimming, the doctrine holds that the swimmer must raise his body to the highest level possible above the water. Two explanations are given for this swimming doctrine: air drag is much smaller than water friction, and the water wave generated by the body propels the body better(Collela, 2009). The doctrine is correct but for a different reason, which is evident in Eqn 7. When the body is high above the water it falls faster (and forward) when it reaches the water line. For the same reason, the speeds of all water waves exhibit the same scale as in Eqn 7, in which Lb is the length scale of the wave(Prandtl, 1969). The crest of the wave falls with a speed of order(gLb)1/2, which becomes visible as the forward speed of the traveling wave.

In the future, the fastest athletes can be expected to be heavier and taller. If the winners' podium is to include athletes of all sizes, then speed competitions might have to be divided into weight categories. This is not at all unrealistic in view of the body force scaling(Eqn 3), which was recognized from the beginning in the structuring of modern athletics. Larger athletes lift, push and punch harder than smaller athletes, and this led to the establishment of weight classes for weight lifting, wrestling and boxing. Larger athletes also run and swim faster.

Here is a brief summary of the scale analysis of animal locomotion, which leads to Eqns 1, 2, 3, 4. It was first done for flying(Bejan, 2000) and then generalized to all locomotion: running, flying and swimming(Bejan and Marden, 2006).

The animal body has a single length scale (Lb). Its mass scale is MρLb3. Locomotion is a rhythm – a sequence of cycles. During each cycle the body must perform work in two ways, in the vertical direction (W1) and in the horizontal direction (W2). Both W1and W2 are destroyed. The vertical work is necessary in order to lift the body to a height of order Lb:
\[\ W_{1}{\sim}M\mathbf{g}L_{\mathrm{b}}.\]
The horizontal work is necessary in order to push through and penetrate the surrounding medium. If the body length and speed scales (i.e. the Reynolds number) are large enough, then the horizontal work is of order:
\[\ W_{2}{\sim}F_{\mathrm{drag}}L_{\mathrm{x}},q\]
\[\ F_{\mathrm{drag}}{\sim}{\rho}_{\mathrm{m}}V^{2}L_{\mathrm{b}}^{2}C_{\mathrm{D}},\]
where Fdrag is the drag force, ρm is the density of the medium, CD∼1 is the drag coefficient and Lx is the distance traveled during the cycle. Together, these scales allow us to estimate the total work per distance travelled:
\[\ \frac{W_{1}+W_{2}}{L_{\mathrm{x}}}{\sim}\frac{M\mathbf{g}L_{\mathrm{b}}}{L_{\mathrm{x}}}+{\rho}_{\mathrm{m}}V^{2}L_{\mathrm{b}}^{2}.\]
The time scale of the cycle is the Galilean time of free fall from the height Lb:
\[\ t{\sim}\left(\frac{L_{\mathrm{b}}}{\mathbf{g}}\right)^{1{/}2}.\]
The horizontal travel during the cycle is LxVt, and Eqn A4 becomes:
\[\ \frac{W_{1}+W_{2}}{L_{\mathrm{x}}}{\sim}M\mathbf{g}^{3{/}2}\frac{L_{\mathrm{b}}^{1{/}2}}{V}+{\rho}_{\mathrm{m}}L_{\mathrm{b}}^{2}V^{2}.\]
The right side is a sum of two terms combined as A/V+BV2, where A and B are two constants and V may vary. This sum is minimal when V∼(A/B)1/3, which yields:
\[\ V{\sim}\left(\frac{{\rho}}{{\rho}_{\mathrm{m}}}\right)^{1{/}3}\mathbf{g}^{1{/}2}{\rho}^{-1{/}6}M^{1{/}6}.\]
The frequency associated with this cycle is t–1∼(g/Lb)1/2or:
\[\ t^{-1}{\sim}\mathbf{g}^{1{/}2}{\rho}^{1{/}6}M^{-1{/}6}.\]
The necessary body force scale F is dictated by the lifting work W1FLbMgLbtherefore:
\[\ F{\sim}M\mathbf{g}.\]
The minimum work per distance traveled is obtained by substituting Eqn A7 and Lb∼(M/ρ)1/3 into Eqn A6:
\[\ \left(\frac{W_{1}+W_{2}}{L_{\mathrm{x}}}\right)_{\mathrm{min}}{\sim}\left(\frac{{\rho}_{\mathrm{m}}}{{\rho}}\right)^{1{/}3}M\mathbf{g}.\]

In conclusion, the scaling relations (Eqns 1, 2, 3, 4) have been derived here in Eqns A7, A8, A9, A10. The modifying factor(ρm/ρ)1/3 depends on the medium. In flying, theρ m (air) is roughly equal to ρ/103, and the factor (ρm/ρ)1/3 is close to 1/10. In swimming,the ρm (water) is the same as the body density, and the factor(ρm/ρ)1/3 is 1. In running, the modifying factor is between 1/10 and 1, and depends on the nature of the running surface and the speed. For example, running through snow, mud and sand is represented by a(ρm/ρ)1/3 value close to 1. Running at high speed on a dry surface is represented more closely by a(ρm/ρ)1/3 factor similar to the one that represents flying.

In summary, the effect of the factor(ρm/ρ)1/3 is weak and of order 1, and for this reason it was left out of Eqns 1and 4. Important to note is that(ρm/ρ)1/3 differentiates between locomotion media in an unmistakable direction: if M is fixed, speeds increase in the direction sea → land → air (cf. Eqn A7); the work requirement decreases in the same direction (cf. Eqn A10). The animal speeds collected over the M=10–3–103 kg range in fig. 2 of Bejan and Marden(Bejan and Marden, 2006)confirm the differentiating effect of the surrounding medium. Each cloud of data is approximated by:
\[\ V{\sim}10M^{1{/}6}(\mathrm{flyers}),\]
\[\ V{\sim}4M^{1{/}6}(\mathrm{runners}),\]
\[\ V{\sim}1M^{1{/}6}(\mathrm{swimmers}),\]
It can be verified that Eqn A11, A12, A13 agree with Eqn 1 within a factor of order 1. They also agree with Eqns 5 and 6.

The data collected in fig. 2of Bejan and Marden (Bejan and Marden,2006) also confirm that the medium factor(ρm/ρ)1/3 does not have an effect on the frequencies (t–1) and forces (F) of flyers,runners and swimmers. This is in accordance with Eqns A8 and A9, in which(ρm/ρ)1/3 is not present.


  • CD

    drag coefficient

  • F

    force, N

  • Fdrag

    drag force, N

  • g

    gravitational acceleration, m s–2

  • H

    body height, m

  • L

    body width, m

  • Lb

    single body length scale, m

  • Lx

    distance traveled, m

  • M

    body mass, kg

  • P

    probability of true null hypothesis

  • R2

    coefficient of determination

  • S

    body slenderness, H/L

  • t

    period, s

  • t

    time, years

  • t–1

    frequency, s–1

  • V

    speed, m s–1

  • W

    ork, useful energy, J

  • W1

    work done vertically, J

  • W2

    work done horizontally, J

  • ρ

    body density, kg m–3

  • ρm

    density of the medium, kg m–3

J.C. is the starting 100 m breaststroke swimmer on Duke University's NCAA swimming team. A.B.'s research on the constructal law of design in nature is supported by grants from the US Air Force Office of Scientific Research and the National Science Foundation.

Ahlborn, B. K. (
Zoological Physics
. Berlin: Springer.
Arnott, S. A., Neil, D. M. and Ansell, A. D.(
). Tail-flip mechanism and size-dependent kinematics of escape swimming in the brown shrimp Cangon crangon.
J. Exp. Biol.
Bartholomew, G. A. and Casey, T. M. (
). Oxygen consumption of moths during rest, pre-flight warm-up, and flight in relation to body size and wing morphology.
J. Exp. Biol.
Bejan, A. (
Shape and Structure: From Engineering to Nature
. Cambridge: Cambridge University Press.
Bejan, A. (
). The constructal law of organization in nature: tree-shaped flows and body size.
J. Exp. Biol.
Bejan, A. and Lorente, S. (
Design with Constructal Theory
. Hoboken, NJ:Wiley.
Bejan, A. and Marden, J. H. (
). Unifying constructal theory for scale effects in running, swimming and flying.
J. Exp. Biol.
Brett, J. R. (
). The relation of size to rate of oxygen consumption and sustained swimming speed of sockeye salmon(Onocorhyncus nerka).
J. Fish. Res. Bd. Canada
Childress, S. and Dudley, R. (
). Transition from ciliary to flapping mode in swimming mollusc: flapping flight as a bifurcation in Reω.
J. Fluid Mech.
Collela, D. (
). Duke University Swimming Coach, Interview.
Drucker, E. G. and Jensen, J. S. (
). Pectoral fin locomotion in striped surfperch. II. Scaling swimming kinematics and performance at gait transition.
J. Exp. Biol.
Greenewalt, C. H. (
). The flight of birds:the significant dimensions, their departure from the requirements of geometrical similarity, and the effect on flight aerodynamics of that departure.
Trans. Am. Philos. Soc.
Heglund, N. C., Taylor, C. R. and McMahon, T. A.(
). Scaling stride frequency and gait to animal size: mice to horses.
Hoppeler, H. and Weibel, E. R. (
). Scaling functions to body size: theories and facts, special issue.
J. Exp. Biol.
Hpathy (
Ideal body weight for male and female adults
, Retrieved March 28, 2009.
Iriarte-Diaz, J. (
). Differential scaling of locomotor performance in small and large terrestrial mammals.
J. Exp. Biol.
Kiceniuk, J. W. and Jones, D. R. (
). The oxygen transport system in trout (Salmo gairdneri) during sustained exercise.
J. Exp. Biol.
Lighthill, J. (
Aeordynamic Aspects of Animal Flight
. British Hydrodynamics Research Association, 5th Fluid Science Lecture, June 1974, p.
Marden, J. H. and Allen, L. R. (
). Molecules, muscles, and machines: universal characteristics of motors.
Proc. Natl. Acad. Sci. USA
Marden, J. H., Wolf, M. R. and Weber, K. E.(
). Aerial performance of Drosophila melanogasterfrom populations selected for upwind flight ability.
J. Exp. Biol.
Marsh, R. L. (
). Ontogenesis of contractile properties of skeletal muscle and sprint performance in the lizard Dipsosaurus dorsalis.
J. Exp. Biol.
May, M. L. (
). Dependence of flight behavior and heat production on air temperature in the green darner dragonfly, Anax junius (Odonata: Aeshnidae).
J. Exp. Biol.
Muller, U. K. and van Leeuwen, J. L. (
). Swimming of larval zebrafish: ontogeny of body waves and implications for locomotory development.
J. Exp. Biol.
NHANES Study (
National Health and Nutrition Examination Survey
, Average weight for an adult man,1999-2002.
Peake, S. J. and Farrell, A. P. (
). Locomotory behaviour and post-exercise physiology in relation to swimming speed, gait transition and metabolism in free-swimming smallmouth bass(Micropterus dolomieu).
J. Exp. Biol.
Pennycuick, C. J. (
). On the running of the Gnu (Connochaetes taurinus) and other animals.
J. Exp. Biol.
Plastic Soldier Review (
Human height
,,Retrieved May 5, 2008.
Prandtl, L. (
Essentials of Fluid Dynamics
. London: Blackie.
Rohr, J. J. and Fish, F. E. (
). Strouhal numbers and optimization of swimming by odontocete cetaceans.
J. Exp. Biol.
Soong, T. T. (
Fundamentals of Probability and Statistics for Engineers
. Hoboken, NJ:Wiley.
Taylor, G. K., Nudds, R. L. and Thomas, A. L. R.(
). Flying and swimming animals cruise at a Strouhal number tuned for high power efficiency.
Tennekes, H. (
The Simple Science of Flight: From Insects to Jumbo Jets
. Cambridge, MA: MIT Press.
Videler, J. J. (
Fish Swimming
. London: Chapman & Hall.
Vogt, P. (
Dictionary of Statistics and Methodology
. Thousand Oaks, CA: SAGE Publications.
Wakeling, J. M. (
). Dragonfly flight III. Lift and power requirements.
J. Exp. Biol.
Weibel, E. R. (
Symmorphosis:On Form and Function in Shaping Life
. Cambridge, MA: Harvard University Press.