Are there absolute limits to the speed at which animals can run? If so, how close are present-day individuals to these limits? I approach these questions by using three statistical models and data from competitive races to estimate maximum running speeds for greyhounds, thoroughbred horses and elite human athletes. In each case, an absolute speed limit is definable, and the current record approaches that predicted maximum. While all such extrapolations must be used cautiously, these data suggest that there are limits to the ability of either natural or artificial selection to produce ever faster dogs, horses and humans. Quantification of the limits to running speed may aid in formulating and testing models of locomotion.

Legged locomotion is a complicated process (for reviews, see Alexander, 2003; Biewener, 2003). As it walks or runs, an animal periodically accelerates both its limbs and its center of gravity. These accelerations require the coordinated application of forces by muscles and skeletal `springs', and the mechanical and neural coordination of these forces can be complex. In turn, the acceleration of the body's various masses and the contraction of muscles place stresses on an organism's skeleton that can be potentially harmful. The metabolic demands of locomotion vary with the morphology, size, speed and gait of the animal.

The extent of our understanding of the complex process of legged locomotion can be assessed through a variety of metrics. One common index is a comparison between measured and predicted maximum speeds: if we understand the physiology and mechanics of locomotion in a particular animal, we should be able to accurately predict how fast that animal can run. Several approaches have been applied to this task. Maximum running speeds have been predicted based on (1)the mass of the body or of its locomotory musculature (e.g. Hutchison and Garcia, 2002; Weyand and Davis, 2005), (2)the rate at which energy can be provided to the limbs(Keller, 1973), (3) the ground force muscles can produce (e.g. Weyand et al., 2000), (4) the stiffness of the `spring' formed by the muscles, ligaments and skeleton (e.g. Farley, 1997), (5) the aerobic capacity of the lungs and circulatory system (e.g. Jones and Lindstedt, 1993; Weyand et al., 1994), and (6)the strength of bones, ligaments and tendons (e.g. Biewener, 1989; Biewener, 1990; Blanco et al., 2003; Iriarté-Diaz, 2002). All of these factors vary with body size, limb morphology and the distance over which speed is measured.

To assess the predictive accuracy of these models, we need empirical standards with which their predictions can be compared. The more precise the models become, the more precise these standards need to be. Therein lies a problem. For extinct species [e.g. Tyrannosaurus(Hutchison and Garcia, 2002)]it is impossible to measure speed directly. For many extant species, maximum running speed has been measured for only a few individuals and often under less than ideal conditions. Even under the best of circumstances the accuracy of speed measurements is often highly questionable (e.g. Alexander, 2003).

Human beings provide an illustrative case study. In many respects, humans are ideal experimental animals for the measurement of maximum speeds. They are intelligent and highly motivated to accomplish a given task, and a vast number of speed trials (competitive races) have been conducted over the years at a wide variety of distances. However, despite this wealth of experimental data,it has proven difficult to quantify the maximum running speeds of humans. In large part, this difficulty is due to the fact that measured speeds have changed through time. For example, the world record speed for men running 1500 m is 14% faster today than it was a century ago(Quercetani and Magnusson,1988; Lawson,1997) (International Association of Athletics Federations, http://www.iaaf.org/statistics/records/inout=O/index.html),and speed in the marathon (42.2 km) is nearly 23% faster than it was in 1920. Increases in speed among women are even more dramatic: 21% faster in the 1500 m race since 1944 and more than 60% faster in the marathon since 1963. If maximum running speed changes through time, it is difficult to use it as a standard for comparison.

The same problem applies to other well-studied species. Horses and dogs have raced competitively for centuries, and one might suppose that their maximum speeds would be well established. But, as with humans, the speed of horses and dogs has increased through time. The winning speed in the Kentucky Derby (a race for thoroughbred horses) increased by more than 13% in the 65 years from 1908 to 1973. The winning speed in the greyhound English Grand National has increased by nearly 15% in the 80 years since its inception in 1927.

The temporal variability in maximum running speeds of dogs, horses and humans raises two central questions. (1) Is there a definable maximum speed for a species running a given distance? The upward trend through time in race speeds for dogs, horses and humans demonstrates that advances in training and equipment, and evolution of the species itself (through either natural selection or selective breeding), can increase running performance. But improvements of the magnitude observed over the last century cannot continue indefinitely: for any given distance, any species will eventually reach its limits. However, it remains to be seen whether this limit can be reliably and accurately measured. (2) If there is a limit to speed, how does it compare with the speed of extant animals? Greyhounds and thoroughbreds have been the subject of intensive selective breeding. How successful has this breeding been in producing fleet animals? Humans are not bred for speed (at least not in the formal fashion of dogs and horses), but what they lack in breeding they have attempted to make up for with improvements in training, nutrition and equipment, and through the use of performance-enhancing drugs. How successful have these efforts been in producing faster humans?

Past attempts at predicting progress in running performance (primarily in humans) have been less than satisfying. Early attempts documented historical trends and extrapolated these trends linearly into the future (e.g. Whipp and Ward, 1992; Tatem et al., 2004). However,these analyses provide no hint of the absolute limits that must exist. Taken to their logical extremes they make absurd predictions: negative race times,speeds in excess of Mach 1. More recent efforts fit historical trends to exponential and logistic equations, models capable of explicitly estimating maximum speeds (e.g. Nevill and Whyte,2005). However, these models have been applied only to humans running a few distances and, even then, only to world records, a small subset of the available data.

Here I use the statistics of extremes and three statistical modeling approaches to estimate maximum running speeds for greyhounds, thoroughbred horses and elite human athletes.

I define `speed' as the average speed that an organism maintains over a fixed distance on level ground (=event distance/total elapsed time). This definition avoids the complications inherent in measuring instantaneous speed(e.g. Tibshirani, 1997). For thoroughbreds and greyhounds, extensive reliable data are available for only narrow ranges of distance (1911–2414 m for horses and 460–500 m for greyhounds), but for humans, data are available for distances varying by more than two orders of magnitude (100–42,195 m).

Historical records

Horses

Data were obtained for the US Triple Crown races (the Kentucky Derby,Preakness Stakes and Belmont Stakes), contested by 2 year olds. The Kentucky Derby has been run annually since 1875, but the current race distance (1.25 miles, 2012 m) was not set until 1896. I obtained winning race times for the years 1896–2008 from the race's official web site(www.kentuckyderby.com/2008/history/statistics). The Preakness Stakes has been run annually since 1873, and at the current race distance (1.1875 miles, 1911 m) since 1925. Data for 1925–2008 were obtained from the race's official web site(www.preakness-stakes.info/winners.php). The Belmont Stakes has been contested annually since 1867. The current race distance (1.5 miles, 2414 m) was set in 1926, and I obtained data for 1926–2008 from the official web site of the New York Racing Association(www.nyra.com/Belmont/Stakes/Belmont.shtml). Unlike racing dogs and humans, racing horses carry a jockey. The weight of the jockey and saddle in the Triple Crown races is typically 55–58 kg.

Dogs

Professional greyhound racing was established in Great Britain in 1927, and three races have been contested annually since that date (with a gap during World War II). Winning race times for three premier dog races (the English Derby, English Grand National and English Oaks, named after horse races) were obtained from www.greyhound-data.com. The length of each race has varied occasionally; only data for races of 460–480 m were used.

Humans

The ATFS (Association of Track and Field Statisticians, 1951–2006), Quercetani and Magnusson (Quercetani and Magnusson,1988), Magnusson and colleagues(Magnusson et al., 1991),Kök and colleagues (Kök et al.,1995; Kök et al.,1999) and the International Association of Athletics Federations(http://www.iaaf.org/statistics/records/inout=O/index.html;accessed June 2008) document annual world's best race times for both men and women for races varying from 100 m to the marathon (42,195 m). Unlike thoroughbreds (for which races are typically among individuals of a single year class), humans often race competitively for several years. On occasion, a single individual has recorded the world's best time in more than one year. This year-to-year connection is a problem for the methods applied here to the analysis of maximum speeds (extreme-value analysis, see below), which assume that individual data points are independent. To minimize this factor, only an individual's best time was used; other years in which that individual recorded the world's best time were deleted from the record. In cases where an individual recorded identical best times in separate years, the first occurrence of the time was used.

Data for men's races are available for some distances beginning in 1900,and for all distances from 1921. Except in the 100 m and 200 m races and the marathon (which include data from 2008), all records extend to 2007. Data for women's races from 100 m to 800 m are available from 1921 to 2007, but data for longer distances are more constrained. Data for the 1500 m race begin in 1944, with a gap from 1949 to 1966, sufficient for the present purposes. Data for the 3000 m, 5000 m and 10,000 m races are too scarce to be useful for this analysis. Records for the women's marathon are available from 1963 to 2007,and are sufficient for analysis.

Several minor adjustments and corrections were made to the human historical record (see Appendix 1).

Analytical approaches

Extreme-value analysis

Each data point in these historical race records is the annual maximum speed recorded from a select group of highly trained athletes in a given race. As such, each record is a sample of the extreme abilities of the species in question. The statistics of extremes(Gaines and Denny, 1993; Denny and Gaines, 2000; Coles, 2001; Katz et al., 2005) asserts that the distribution of these extreme values should conform asymptotically to a generalized extreme-value (GEV) distribution:
(1)
Here, P(V) is the probability that an annual maximum speed chosen at random is ≤V. The shape of this cumulative probability curve is set by three parameters: a, a shape parameter; b, a location parameter; and c, a scale parameter. Parameters aand b can take on any value; c is constrained to be >0. If a is ≥0, the shape of the distribution of extreme values is such that there is no defined limit to the extremes that can potentially be reached (Coles, 2001). In contrast, if a<0, P=1 when V=b–(c/a). In such a case, the distribution of extreme values has a defined absolute maximal value:
(2)

Because of this ability to define and quantify absolute maxima, the statistics of extremes is a promising approach for the study of maximum running speeds. All extreme-value analyses were carried out using extRemes software (E. Gilleland and R. W. Katz, NCAR Research Applications Laboratory,Boulder, CO, USA), an implementation in the R language of the routines devised by Coles (Coles, 2001). When fitting Eqn 1 to data, there is inevitably some uncertainty in the estimated value of each parameter and of Vmax. Confidence limits for these values were determined using the profile likelihood method (Coles,2001).

Method 1: no trend

Appropriate application of Eqn 1 to the estimate of maximum speeds depends on whether or not there is a trend in race data, either with time or with population size. For some races, speeds appear to have plateaued in recent years. The existence and extent of such a plateau was determined by sequentially calculating the correlation between year and maximum speed starting with data for the most recent 30 years and then extending point by point back in time. If there was no statistically significant trend in speed, a plateau was deemed to exist,and the beginning of the plateau was taken as the year associated with the lowest correlation coefficient. For the data in such a plateau, application of Eqn 1 is straightforward:parameters a, b and c are chosen to provide the best fit to the raw speed data in the plateau, the degree of fit being judged by a maximum likelihood criterion. If a is significantly less than 0 (that is, if its upper 95% confidence limit is <0), the absolute maximum is calculated according to Eqn 2 and confidence limits for this estimate are obtained. In cases where a≥0, no absolute maximum value can be determined. In these cases, I arbitrarily use the estimated maximum value for a return time of 100 years as a practical substitute for the absolute maximum.

Method 2: the logistic model

For cases in which trends are present in race data, the trend must first be modeled before the analysis of extremes can begin. Two models were used. The first model is the logistic equation mentioned in the Introduction:
(3)
which provides a flexible means to model values that, through time, approach an upper limit (Nevill and Whyte,2005). Eqn 3 is a natural fit to a basic assumption made here: that an absolute limit must exist to the speed at which animals can run. In Eqn 3, V(y) is the fastest speed recorded in year y, mn is the model's minimum fastest speed and mx is the model's maximum fastest speed (the parameter of most interest in the current context). k is a shape parameter determining how rapidly values transition from minimum to maximum, and t is a location parameter that sets the year at which the rate of increase is most rapid. Note that although the logistic equation incorporates the assumption that a maximum speed exists, it does not assume that speeds measured to date are anywhere near that maximum.

Historical race data were fitted to Eqn 3 using a non-linear fitting routine with a least-squares criterion of fit (Systat, SPSS, Chicago, IL, USA), providing an estimate(±95% confidence limit) for mx. Note that the confidence limits for mx indicate the range in which we expect to find this parameter of the model, but that this range does not necessarily encompass the variation of the data around the model. A heuristic example is shown in Fig. 1. For this example, I created a hypothetical set of race data by specifying a logistic curve (mn=14,mx=17, k=0.1, t=1940) for y=1900 to 2020 and adding to each deterministically modeled annual maximum speed a random speed selected uniformly from the range –0.35 to 0.35 m s–1. When provided with these hypothetical data, the fitting routine very accurately estimated the true values for the parameters of the model, but the estimated mx (17.00 m s–1, with 95% confidence limits ±0.06 m s–1) does not include the maximum speed in the data set(17.32 m s–1). To accurately estimate the overall maximum(rather than the maximum of the trend), it is necessary to characterize the effect of variation around the fitted logistic model.

To do so, the temporal trend in extreme values is first incorporated into the analysis by subtracting the best-fit logistic model from the measured annual maximum race speeds (Gaines and Denny, 1993; Denny and Gaines,2000). The distribution of the resulting deviations –trend-adjusted extreme values – is then analyzed using Eqn 1. If a for the distribution of trend-adjusted extremes is <0, there is a defined absolute maximum deviation from the logistic model(Eqn 2), and this absolute maximum deviation is added to the estimated upper limit of the logistic model(the best-fit mx) to yield an estimate of the overall absolute maximum speed. As before, if a≥0, I estimate a practical maximum by adding the 100 year return value to mx. In both cases, 95% confidence limits on the predicted maximum deviation give an indication of the statistical confidence in the overall estimate of maximum speed.

Unlike the logistic analysis of Nevill and Whyte(Nevill and Whyte, 2005),which used only world-record speeds, the analysis here uses the much more extensive measurements of annual maximum speed. The variation of these annual maxima around the underlying trend provides insight into the random variation in maximum speed present in the population of animals under consideration. And unlike the analysis of data in just the plateau portion of a record (method 1), analysis using deviations from the logistic equation utilizes all the data in the historical record.

Method 3: population-driven analysis

The larger the population from which racing dogs, horses or humans are selected each year, the higher the probability that an exceptional runner will be found by chance alone (Yang,1975). Thus, if population size grows through time, maximum recorded speed may increase. To assess the effect of this interaction between speed and population, I used the following analysis, illustrated here using human females as an example.

I begin by assuming that there exists an idealized distribution of the speeds at which individual women run a certain distance. Every time a woman runs a race at that distance, her speed is a sample from this distribution. I then suppose that in a given year S0 women run the race,and we record V, the fastest of these speeds. V is thus one sample from a different distribution – the distribution of maximum speeds for women running this distance. A different random sample of S0 individual speeds would probably yield a different V. In this fashion, repeated sampling provides information about the distribution of maximum speeds. Our job is to analyze the measured values of V to ascertain whether the distribution of these maximum speeds has a defined upper bound. Statistical theory of extremes(Coles, 2001) suggests that if we confine our exploration to maximum speeds above a sufficiently high threshold u, the cumulative distribution of our sampled maxima asymptotes to the generalized Pareto family of equations (GPE):
(4)
Here G(V) is the probability that the maximum speed of a sample of S0 women is >V.
Parameters ϵ and u can take on any value, σ is constrained to be >0. Note, if ϵ<0, G(V)=0 when V=u–(σ/ϵ). In other words, if ϵ is negative:
(5)
Thus, if we have empirical data for V and G(V), we can estimate u, σ and ϵ, and potentially calculate Vmax (Fig. 2A).
We have historical data for annual maximum speeds, but how can we estimate the corresponding probabilities, G? Here population size comes into play. The larger the population of running women available to be sampled in a given year, the greater the chance of having a woman run at an exceptionally high speed. Put another way, the larger the S, the population of female runners, the faster the V that we are likely to record. The faster the V we record in a large sample, the lower the probability that we would have exceeded that high V in a small sample. [Recall that G(V) is defined specifically for sample size S0.] Thus, V recorded at large S should have a relatively low G. The sense of this relationship between population size and G(V) is depicted in Fig. 2A (note the ordinate on the right). The theoretical relationship between population size and annual maximum velocity is described by a modification of the GPE:
(6)
Here N(V) is world population size in the year in which V was measured and N0 is world population size at the beginning of the historical record. For a derivation of this equation, see Appendix 2.

Note that total world population is used in Eqn 6 solely as an index of the population of runners. The actual number of runners is some unknown fraction f of the total population. Both as a practical matter and for the sake of simplicity, I assume that f is constant through time, in which case it cancels out of the equation when the ratio of population sizes is taken. I recognize the likelihood that, in fact, f has varied through time, and the potential effects of variation in f are addressed in the Discussion.

I apply this approach to the historical records of running speeds for races in which speed is correlated with population size. From records of population size and annual maximum speed, I construct an exceedance distribution of annual maximum speeds as described by Eqn 6, for which I then calculate the best fit values of ϵ, u, and σ using a least-squares criterion of fit. (Note that in this analysis, the value of the threshold u is chosen to give the best fit to the GPE rather than being chosen a priori.) The 95%confidence limits, for both individual parameters and the overall distribution, are calculated using 2000 iterations of a bootstrap sampling of the data with accelerated bias correction(Efron and Tibshirani, 1993). If ϵ≤0, the annual maximum speed estimated for an infinite population(G=0) provides an estimate of absolute maximum speed(Eqn 5), similar in principle(although not necessarily in magnitude) to the maximum speed estimated from the logistic equation (method 2).

To account for random variation about this best-fit population-driven model, deviations in speed from the best-fit GPE model are analyzed using Eqn 1. If the best-fit afor this distribution is <0, an absolute maximum deviation exists(Eqn 2), and this maximum deviation is added to the expected value for the population-driven model.

The similarity between the population-based model of maximum speed and the time-based logistic model is highlighted in Fig. 2B where the information of Fig. 2A has been replotted to show how modeled speed (the red line) varies as a function of population size.

As an alternative to the approach used here, the effect of population size on running speed could be addressed by incorporating population size as a covariate in a GEV analysis similar to that used in method 1 (see Coles, 2001; Katz et al., 2005). For many human races, speed increases approximately linearly with the logarithm of population size. Using log population size as a linear covariate produces results essentially similar to those obtained with method 1 described above. However, this alternative method has not been fully explored, and its results are not reported here.

Population size

To implement the population-driven model, information is required regarding the year-by-year population of potential contestants. The combined number of thoroughbred foals born each year in the United States, Canada and Puerto Rico has been recorded by the US Jockey Club(www.jockeyclub.com/factbook08/foalcrop-nabd.html),and I assume that this represents a good approximation of the potential population from which Triple Crown runners are chosen. Horses racing the Triple Crown are 2 yearolds, so the foal crop from a given year represents the potential racing population of the following year.

A request to the English Stud Book for records of racing greyhounds born in the UK each year was not answered, so an estimate of the trend in the potential population of racing greyhounds was obtained from records of the Irish Stud Book. This substitution was deemed acceptable for two reasons. First, many (perhaps most) greyhounds racing in Britain are born in Ireland. Second, Eqn 6 uses the ratio of initial population size to that in a given year. Thus, as long as the number of Irish greyhound births is proportional to that in the UK, the use of the Irish data is valid. Greyhounds begin racing at age 18 months to 2 years; I assume here that the number of dogs registered in a given year is approximately the number available to race in the following year.

Estimates of the world's human population were garnered from Cohen(Cohen, 1995) and the US Census Bureau(http://www.census.gov/ipc/www/idb/worldpop.html;`Total midyear population for the world: 1950–2050'; accessed June,2008). Population size N for a given year from 1850 to 2008 was estimated as:
(7)
where N is measured in billions and X is centuries since 1800 (r2>0.999). I assume a 1:1 gender ratio for humans; the potential runners' population for men and women is thus each half the total world population. Note that this equation yields spurious values if used outside the years 1850–2008.

Thoroughbreds

Temporal patterns of winning speeds for the US Triple Crown are shown in Fig. 3. There is no significant correlation between year and winning speed in the Kentucky Derby for the period 1949 to 2008. An apparent plateau was reached later in the Preakness Stakes (1971) and Belmont Stakes (1973). Significance levels for the correlation of speed with time in these plateau years are given in Appendix 3, Table A1.

Extreme-value analysis of race speeds during each plateau suggests that there is an absolute upper limit to speeds in each of these races (ais significantly less than 0 in each case), and the predictions are shown in Table 1 (no-trend model).

The temporal pattern of speed for each race is closely modeled by the logistic equation (Appendix 2 and Table A2) and, in each case, horses appear to have reached a plateau in speed (Fig. 3). Predictions of maximum speed made using a logistic model fitted to the entire data set for each race (Table 1, logistic model) are statistically indistinguishable from those obtained from the subset of plateau data.

For each race, the predicted absolute maximum running speed (averaged across methods) is only slightly (0.52% to 1.05%) faster than the current record.

The potential population of Triple Crown runners increased dramatically from the 1880s until the mid-1980s, but has decreased since(Fig. 4A). Plots of speed as a function of population size (e.g. Fig. 5A) demonstrate that changes in population size are not the controlling factor in winning speeds in these horse races. Speed in the Kentucky Derby is not correlated with population size when the population is above 8400 (P>0.882, Fig. 5A), in the Preakness Stakes when the population is greater than 24,300 (P>0.867), and in the Belmont Stakes when the population is greater than 25,700 (P>0.253). Because of the lack of correlation between population size and speed above certain population limits, the population-driven model of speeds was not applied to horses.

Greyhounds

In a pattern similar to that seen with horses, race speeds for greyhounds appear to have plateaued (Fig. 6). The plateau in the English Oaks began in approximately 1966 and in the English Grand National and English Derby in approximately 1971. The significance levels of the regression of speed on time during these plateaus are given in Appendix 3 (Table A3). The temporal pattern of speeds for each race is closely modeled by the logistic equation (details are given in Appendix 3, Table A4), and predicted maximum speeds calculated using this method(Table 2, logistic model) are very similar to those calculated from the plateaus alone.

Averaged across methods, predictions of maximum running speed for each race are only 0.29% to 0.92% higher than existing records(Table 2).

The estimated population of racing dogs increased gradually from 1950 (the earliest year in which records are available) to 2007(Fig. 4B), with substantial year-to-year variation. Speed in the English Oaks is not correlated with population size when the population is above 14,000 (P>0.354), in the English Grand National when the population is greater than 19,300(P>0.332), and in the English Derby when the population is greater than 19,565 (P>0.232). A representative example of the relationship between Irish greyhound population size and speed is shown in Fig. 5B. The fact that race speeds apparently plateaued while the population increased suggests that population size is not a substantial factor in the control of maximal speed in greyhounds, and the population-driven method of analysis was not applied to dogs.

Humans

Temporal patterns of human running speed are shown in Figs 7, 8, 9. For women running 100 m to 1500 m, speeds appear to have plateaued during the 1970s. Approximate onset years for each plateau and the corresponding probability level are given in Appendix 3 (Table A5). For these races, I applied extreme-value analysis directly to the data in each plateau. In the 200 m and 800 m races, an absolute maximum speed could be calculated (a was significantly <0, Table A5). In the 100 m, 400 m and 1500 m races, no absolute limit is defined for women's speeds; 100 year return values are given here (Table 3) and absolute maxima (if they exist) will be somewhat higher.

Data for all human races could be accurately fitted with a logistic model(for details see Appendix 3, Table A6). Results from the logistic models of human running speed are given in Table 3. In all cases except the women's 400 m and 1500 m races, a in the GEV fit was<0, and absolute maximum speeds could be predicted (Appendix 3, Table A6). For the women's 400 m and 1500 m races, values for 100 year return times are used. For the races in which speeds appear to have plateaued, predictions made using the logistic equation and the entire historical record are slightly higher than those from the analysis of the plateaus alone (Table 3).

The human population has exploded over the last century(Fig. 4C). In those races in which women's speeds have reached a plateau in recent years, a similar plateau is present in the relationship between speed and population (see Appendix 3, Table A7), and therefore it is unlikely that population is driving speed in these races. In all non-plateau races, however, a plot of speed versus population size shows a correlation throughout the record (a representative example is shown in Fig. 5C), and the population-driven model was applied. Due to large year-to-year variation in speeds recorded early in the twentieth century, the GPE fitted to data from men's 100 m and 200 m races had exceptionally large confidence intervals (e.g. the 95% confidence interval for predicted absolute maximum speed included 0 m s–1), and these questionable results are not included here. In the other races analyzed, the GPE provided an acceptable fit. Results from a representative example are shown in Fig. 10, and all results from this model are given in Table 3. In all non-plateau races, a was <0 (Appendix 3, Table A8), and an absolute maximum deviation from the fitted trend was calculated. Estimates from the population-driven model for non-plateau races closely match those obtained from other analytical approaches (Table 3), suggesting that increasing human population size will not be a major factor in future track records.

The results from all human races are summarized in Fig. 11 and Table 4. Speeds for which 100 year maxima are used (rather than absolute maxima) are shown as open symbols. Average predicted maximum speeds for men and women are only modestly faster than current world records (1.06% to 5.09% for men, 0.36% to 2.38% for women). The predicted potential for an increase in speed is not significantly correlated with race distance for men (P>0.93). There is a marginally significant negative correlation between the potential for increase and race distance in women (P=0.037), but this correlation is driven solely by the low predicted increase in speed in the marathon. Predicted maximum speeds for women are 9.3% to 13.4% slower than those for men, and in all but one instance (the 400 m race) the predicted maximum speed for women(including the confidence intervals) is less than the current record speed for men. There is a significant difference in the mean scope for increase(predicted maximum speed divided by current world record speed, minus 1)between men and women. For data pooled across all speeds, the mean scope for increase in predicted speed is 3.17% for men and 1.55% for women (Student's t-test, unequal variances, d.f.=12, P=0.008).

These results provide tentative answers to the questions posed in the Introduction. For greyhounds, thoroughbreds and humans, there appear to be definable limits to the speed at which they can cover a given distance, and current record speeds approach these predicted limits. If present-day dogs,horses and humans are indeed near their locomotory limits, these animals (and the limits they approach) can serve as appropriate standards against which to compare predictions from mechanics and physiology.

The case for defined limits in horses and dogs is particularly strong. Despite intensive programs to breed faster thoroughbreds and greyhounds,despite increasing populations from which to choose exceptional individuals,and despite the use of any undetected performance-enhancing drugs, race speeds in these animals have not increased in the last 40–60 years. Thus, for horses and dogs, a limit appears to have been reached, subject only to a slight (and bounded) further increase due to random sampling. The situation is less clear cut for humans, in particular for men. Logistic and population-driven models of the historical data suggest that a limit to male human speed exists, and that this speed is only a few per cent greater than that observed to date. But unlike speeds in horses and dogs, and sprint speeds for women, speeds for men have not yet reached a plateau.

An excellent example of the potential for a continued increase in men's speeds is provided by the recent world records set in the 100 m and 200 m races by Usain Bolt of Jamaica. Over a span of 3 days in the Olympic games of 2008, Bolt `shattered' the then existing records, lowering the record in the 100 m from 9.72 to 9.69 s and in the 200 m from 19.32 to 19.30 s. Because Bolt is exceptionally tall for a sprinter (6′5″, 1.96 m), he was hailed by the press as a physical `freak' and the harbinger of a new era of sprinting.

Should Bolt's records cast doubt on the predictions made here? The answer is no. Bolt's records are only small improvements on the existing records for the 100 m and 200 m races, 0.3% and 0.1%, respectively, and Bolt's records are not out of line with the logistic fit to the historical data (Figs 7 and 8, pink dots). Furthermore,there have previously been similar jumps in record speed. Thus, as admirable as they are, there is nothing in Bolt's records to suggest that the predictions made here are inaccurate or that human speeds in the 100 m and 200 m races are limitless.

For distances of 100 m to 1500 m, women's speeds appear to have plateaued(Figs 7 and 8), superficially giving added confidence in the logistic model of the data. These plateaus (and this confidence) should be viewed with some skepticism, however. In each of these races, the current world record was set in the early to mid-1980s, a time when performance-enhancing drugs were becoming prevalent in women track athletes but before reliable mechanisms were in place to detect these drugs(Holden, 2004; Vogel, 2004). If speeds were artificially high in the 1980s due to drug use, and drug use was absent in subsequent years, one might suspect that the apparent plateau in the historical record could be an artifact. However, removing the annual maximum speeds from the 1980s does not substantially alter the results of the logistic analyses. Thus, the temporal plateaus in speeds at these distances appear to be real. There is an interesting corollary to this conclusion: if performance-enhancing drugs are still being used by women, the effect of the drugs has itself reached an apparent plateau.

In contrast to times in the shorter races, women's speeds in the marathon have continued to increase in recent years; like men, women in this race have not reached a demonstrable plateau. In this case, however, the current world record (5.19 m s–1 set by Paula Radcliffe in 2005) is very close to the average predicted absolute maximum speed (5.21 m s–1); indeed, the current world record exceeds the maximum predicted from the population-driven model (although it lies within the confidence interval of this estimate). Given the upward trend in recent marathon speeds and the small difference between the current record and the predicted limit, this race is likely to provide the first test of the methods and predictions used here.

Maximum speeds predicted here are on average 1.63% higher than those predicted by Nevill and Whyte (Nevill and Whyte, 2005). The difference is probably due to the fact that the implementation of the logistic method here takes into account sampling variation in the maximum speeds.

My results for humans bolster the conclusion reached by Sparling and colleagues (Sparling et al.,1998) and Holden (Holden,2004) that the present gender gap between men and women will never be closed for race distances between 100 m and the marathon. Note that none of the data presented here speak to the possibility that women may someday out-run men at longer distances (Bam et al., 1997). Similarly, my results are consistent with the conclusion reached by Yang (Yang,1975) that increases in human population will have only a minor effect on speed.

Evolution

Is it reasonable to suppose that the evolution of speed in horses has reached its limits? In a restricted sense, the answer is yes. The equipment used in horse racing, and the surfaces of the tracks on which these races are contested, did not change appreciably during the years when speeds were increasing in the Triple Crown races, and they have not changed since. Nor were there any apparent breakthroughs in training or nutrition that led to the increases in speed in thoroughbreds in the first half of the twentieth century. It seems likely, then, that the initial increase in speed in horses was due primarily to selective breeding. If this is true, evidence from the Triple Crown races suggests that the process of selective breeding of thoroughbreds (as practiced in the US) is incapable of producing a substantially faster horse: despite the efforts of the breeders, speeds are not increasing, and current attempts to breed faster horses may instead be producing horses that are more fragile(Drape, 2008). The fastest speed in two of the Triple Crown races was set in 1973 by the same horse,Secretariat, and he was initially credited with a speed equal to the record in the third race (the Preakness Stakes) as well. (The timer malfunctioned in that race, however, and Secretariat's subsequently established official speed is slightly slower.) Thus, Secretariat approached the predicted absolute maximum speeds in all three of his Triple Crown races and therefore may represent a good approximation of the ultimate individual thoroughbred in races 1.25–1.5 miles long.

In a larger sense, however, the equine data presented here are preliminary at best. It may well be possible that different criteria for selective breeding of horses could produce a faster animal. Thoroughbreds have been recognized as a separate breed since the 1700s, and regulation of the breed has constrained its gene pool: thoroughbreds are less genetically diverse than other breeds of horses (Cunningham et al.,2001). The breed is effectively a closed lineage descended from as few as 12–29 individuals (Cunningham et al., 2001; Hill et al.,2002), and 95% of the paternal lineages in present-day thoroughbreds can be traced to a single stallion, The Darley Arabian. Selective breeding starting with different equine stock could perhaps yield faster horses. In this sense, then, the results presented here do not necessarily address the question of the maximum speed for the species Equus cabillus.

The same arguments apply to greyhounds. Greyhounds have been bred for speed since antiquity, and the results presented here show a clear plateau in the ability of current selective breeding to produce a faster greyhound. However,given the extraordinary malleability of canine morphology (which includes everything from Chihuahuas to Great Danes), it is quite possible that different breeding strategies (perhaps starting with a different breed of dog)could produce a faster Canus familiaris.

Once again, the situation is less clear for humans. Unlike the apparent case in horses and dogs, human runners have recently benefited from substantial improvements in training, equipment and nutrition. In some cases(women's sprints), these benefits may have reached their own limits. But in other races, continued improvement in training, equipment and nutrition may well be contributing to the continued increase in race speeds. Because these effects are inextricably entwined with the historical race data, the predictions I make here may be biased. In essence, these predictions assume that historical trends in training, equipment and nutrition (whatever they are) will continue into the future. It is always dangerous to make such assumptions. Competitive swimming provides an example of this potential effect: improvement in the design of full-body swimsuits (a breakthrough not contemplated 10 years ago) contributed to a rash of recent world records, 25 in the Olympics of 2008 alone. Until we know more about the mechanisms of improvement in training, equipment and nutrition, and more about their actual role in the historical running record, the magnitude of their effects on future running speeds will remain uncertain, and the predictions made here must be used cautiously.

And then there is the subject of artificial performance enhancement, which inevitably leads to philosophical questions pertaining to the definition of absolute maximum speed. For example, how should we define `male', `female' and even `human' for the purposes of this study? If a woman artificially enhances the concentration of testosterone in her body, a large number of changes accrue that make her physiologically more like a man and capable of higher speeds (e.g. Holden, 2004). In this altered state, she may exceed any limits that might exist for unaltered individuals. At what point should such a hormonally enhanced woman no longer count as a woman in the analysis of maximum female human speed? Stanislawa Walasiewicz (Stella Walsh) provides an intriguing example. She was the preeminent female sprinter of the 1930s, posting times that were not matched until the 1950s. She married in 1956 and was subsequently shot and killed during a supermarket robbery, where she was an innocent bystander. An autopsy revealed that (unbeknownst to her) she was a hermaphrodite, possessing both ovaries and testes. Presumably the anabolic steroids produced by her testes contributed to her athletic success(Lawson, 1997), but had she not been murdered, we would never have known. Where should she be placed in the record books? Even more vexing questions await us in the future. For example, the potential exists to genetically engineer human athletes for enhanced performance (e.g. Vogel,2004). At what point does a genetically altered person no longer count as human? (Similar questions can be raised regarding drugs and gender in greyhounds and horses.)

For scientific purposes, these sticky questions can in large part be circumvented through the use of arbitrary definitions. As long as one defines practical criteria for `female', `male' or `human' in formulating a mechanical/physiological model of locomotion, the predictions of that model can be compared with an appropriate set of test organisms. For present purposes, let us define a greyhound, thoroughbred or human (male or female) as an individual performing without drug or genetic enhancement. If drugs have contributed to the winning speeds in the races used here, speeds in the absence of these drugs would presumably have been slower. Thus, if we could definitively account for the effect of drugs in the historical record, the estimated maximum speeds I predict for unadulterated animals would, if anything, be slower, and my estimates are in this sense conservatively fast.

A similar conclusion applies to possible variation in the fraction of the overall human population that participates in running races. Recall that the analysis of the effect of population size on maximum speed assumed that f, a fixed (although unspecified) fraction of all men and women, runs each year. In reality, it seems likely that f has increased through time as awareness of the sport has spread and more prize money has been made available. In particular, the fraction of women competing in running races has probably increased in recent decades. However, if we could take the presumed increase in f into account in this analysis, the effect would be to reduce the predicted maximum speed. (See Appendix 2 for an explanation.) In this respect, my results are again conservatively fast.

Some cynics have suggested that the problem of artificial enhancement in sports should be `solved' by simply making drug and genetic enhancement legal. If such enhancement is allowed, the question of maximum human running speed becomes much more difficult to answer. On the one hand, it seems likely that humans have been, and still are, clandestinely employing performance-enhancing drugs despite the ban on their use. If so, the efficacy of these drugs is questionable: e.g. women's speeds in sprint races appear to have plateaued. On the other hand, it is impossible to rule out the possibility that new drugs or genetic enhancement could do for running what full-body swim suits have done for swimming: provide the means for dramatic improvement. In that case, the maximum speeds estimated here could be low.

What limits speed?

The analysis presented here deals solely with the results of competitive races, not with the factors that caused a certain individual to win or lose. In this respect, my results are as unsatisfying as those of previous statistical analyses: they tell us that speed has limits, but not what accounts for these limits. Nonetheless, the pattern of estimated maximal speeds provides information of potential value to physiologists and biomechanicians. It seems unlikely that a single mechanical or physiological factor could account for the limit to speed at all distances. The height and mass of elite runners differs among race distances(Weyand and Davis, 2005) as does the ability of aerobic capacity to predict speeds(Weyand et al., 1994; Weyand et al., 1999). It is striking, then, that the predicted scope for increased speed in humans is similar across distances ranging from 100 m to 42,195 m(Fig. 11; Table 4). This distance-independent scope for increase suggests that some sort of higher order constraint may act on the suite of physiological and mechanical factors to limit speed.

Context

The likelihood that there are limits to speed should not diminish the awe with which we view the performance of dogs, horses and humans. For example, a women running the estimated absolute fastest speed for 100 m would have beaten the world's fastest male in 1955, a feat that would have astounded contemporary spectators. The predicted maximum speed (5.83 m s–1) for a man running a marathon (42.2 km) would have been fast enough to beat the great Emil Zatopek in his world's best 10 km race in 1954. Those in the stands watching that race could not have imagined someone besting Zatopek by 16 s, and then simply continuing at that winning pace for another 32.2 km.

Notes on the historical records

Horses

Until 1991, Triple Crown races were timed by hand to the nearest 0.2 s;since then they have been timed electronically to the nearest 0.01 s. I have made no correction for the shift from hand to electronic timing.

Humans

In some women's races, the earliest data are extremely variable, perhaps because so few women participated in the sport at that time. For the 100 m race, women's data prior to 1931 were not used. For the 400 m race, data prior to 1928 were not used.

World's best times by a individual later proven to be using performance-enhancing drugs were deleted from the record; in those years, the second-best time was used. Times posted by individuals of uncertain gender (as noted by the ATFS) were not used; second-best times were substituted. The current women's 100 m world record (10.49 s) was set by Florence Griffith Joyner in 1984, but there is compelling evidence that that race was wind aided(Pritchard and Pritchard,1994). I have replaced this world's best for 1984 with a time of 10.61 s (Griffith Joyner's second-best time for that year).

Until the 1970s, human races were timed by hand. Because of the slight delay in starting a watch (due to human response time) and the potential for an early stop as an official anticipates the finish, hand times are slightly shorter than corresponding electronic times by approximately 0.165 s. This difference is negligible for races of 400 m and longer, but can be a substantial factor in 100 m and 200 m races. Here, I have added 0.165 s to hand-timed 100 m and 200 m times. For 2–3 years during the switch-over from hand to electronic timing, world's best times were recorded for both methods, and both hand-timed results (plus 0.165 s) and electronic results have been included here.

A race of 220 yards is very similar in length to a race of 200 m (220 yards=201.168 m). The ATFS has reconciled the two races by recording 220 yard times minus 0.1 s as equivalent to 200 m times for both men and women, and I follow this convention here. Similarly, I accept 440 yard times minus 0.3 s as equivalent to 400 m times.

Relating population size to exceedance probability

The generalized Pareto equation (GPE) is traditionally expressed as a cumulative probability distribution:
(A1)
where P(V) is the probability that the maximum speed in a randomly chosen sample of size S0 is ≤V(Coles, 2001). For the present purposes, it is more convenient to work with 1–P(V)=G(V), where G(V) is the chance of getting a sample maximum speed> V:
(A2)
How many samples of size S0 would one have to take (on average) to get a maximum speed >V? From the laws of probability(Feller, 1968; Denny and Gaines, 2000) it can be shown that this exceedance number E is:
(A3)
If each sample has S0 individuals and we need Esamples (on average) to get a maximum speed >V, the overall number of individuals we need to sample to exceed V is:
(A4)
Rearranging, we solve Eqn A4 for G(V):
(A5)
Upon inserting Eqn A5 into Eqn A2, we see that:
(A6)
If, in a given year, the world population of (for example) women is Nand fraction f are runners, S=fN women run the race that year. When measured in units of S0, the number of`samples' of maximum speed we take in a given year is therefore S/S0. If f is constant across years, S0/S=N0/N where N0 is the world population in the year in which S0 was measured. Thus:
(A7)
In this fashion, N0/N (and hence G) can be estimated from records of world population for each measured V. Plotting G as a function of V provides an estimate of the exceedance function described by Eqn 6 in the text. This is the recipe for estimating u,σ and ϵ discussed in the text.

It is possible that the fraction of the human population sampled by competition (especially the fraction of women) has increased over the last century. If so, maximum speeds predicted using the population-driven model are too high: estimates of cumulative probability incorporating an increasing fraction of runners would yield lower G values for recent years than those shown here, thereby lowering the estimates for extrapolated maximum speed (see Fig. A1).

Details of the statistical results

Details of the statistical results are given in Tables A1, A2, A3, A4, A5, A6, A7, A8.

I thank S. Davis, Secretary of the ATFS; S. Ito of the LA84 Foundation Library; C. Leavy and D. J. Histon of the Irish Stud Book; D. Kohrs and the Minneapolis Public Library for assistance in compiling race data; and J. M. Gosline for suggesting horses and greyhounds as experimental subjects.

Alexander, R. M. (
2003
).
Principles of Animal Locomotion
. Princeton, NJ: Princeton University Press.
Association of Track and Field Statisticians(
1951
-2006).
ATFS Annual
. Worcester, UK:SportsBooks.
Bam, J. T., Noakes, D., Juritz, J. and Steven, S. C.(
1997
). Could women outrun men in ultramarthon races?
Med. Sci. Sports Exerc.
29
,
244
-247.
Biewener, A. A. (
1989
). Support in mammals:limb posture and muscle mechanics.
Science
245
,
45
-48.
Biewener, A. A. (
1990
). Biomechanics of mammalian terrestrial locomotion.
Science
250
,
1097
-1103.
Biewener, A. A. (
2003
).
Animal Locomotion
. Oxford: Oxford University Press.
Blanco, R. E., Gambini, R. and Fariña, R. A.(
2003
). Mechanical model for theoretical determination of maximum running speed in mammals.
J. Theor. Biol.
222
,
117
-125.
Cohen, J. (
1995
).
How Many People Can the Earth Support?
New York: Norton.
Coles, S. (
2001
).
An Introduction to Statistical Modeling of Extreme Values
. London:Springer-Verlag.
Cunningham, E. P., Dooley, J. J., Splan, R. K. and Bradley, D. G. (
2001
). Microsatellite diversity, pedigree relatedness and the contributions of founder lineages to thoroughbred horses.
Anim. Genet.
32
,
360
-364.
Denny, M. W. and Gaines, S. (
2000
).
Chance in Biology
. Princeton, NJ: Princeton University Press.
Drape, J. (
2008
). Breeding for speed, ignoring durability.
New York Times
, May 11.
Efron, B. and Tibshirani, R. J. (
1993
).
An Introduction to the Bootstrap
. Boca Raton, FL: CRC Press.
Farley, C. T. (
1997
). Maximum speed and mechanical power output in lizards.
J. Exp. Biol.
200
,
2189
-2195.
Feller, W. (
1968
).
An Introduction to Probability Theory and its Applications. 3rd Edn.
New York:John Wiley.
Gaines, S. and Denny, M. W. (
1993
). The largest, smallest, highest, lowest, longest and shortest: extremes in ecology.
Ecology
74
,
1677
-1692.
Hill, E. W., Bradley, D. G., Al-Barody, M., Ertugrul, O., Splan,R. K., Zakharov, I. and Cunningham, E. P. (
2002
). History and integrity of the thoroughbred dam lines revealed in equine mtDNA variation.
Anim. Genet.
33
,
287
-294.
Holden, C. (
2004
). An everlasting gender gap?
Science
305
,
639
-640.
Hutchison, J. R. and Garcia, M. (
2002
). Tyrannosaurus was not a fast runner.
Nature
415
,
1018
-1021.
Iriarté-Diaz, J. (
2002
). Differential scaling of locomotor performance in small and large terrestrial animals.
J. Exp. Biol.
205
,
2897
-2908.
Jones, J. H. and Lindstedt, S. L. (
1993
). Limits to maximal performance.
Ann. Rev. Physiol.
55
,
547
-569.
Katz, R. W., Brush, G. S. and Parlange, M. B.(
2005
). Statistics of extremes: modeling ecological disturbances.
Ecology
86
,
1124
-1134.
Keller, J. B. (
1973
). A theory of competitive running.
Phys. Today
26
,
42
-46.
Kök, N., Magnusson, R., Potts, D. H. and Quercetani, R. L. (
1995
).
Track and Field Performances Through the Years
. Vol.
3
,
1945
-1950. Trowbridge, UK: Association of Track and Field Statisticians.
Kök, N., Magnusson, R., Potts, D. H. and Quercetani, R. L. (
1999
).
Track and Field Performances Through the Years
. Vol.
4
,
1921
-1928. Trowbridge, UK: Association of Track and Field Statisticians.
Magnusson, R., Potts, D. H. and Quercetani, R. L.(
1991
).
Track and Field Performances Through the Years
. Vol.
2
,
1937
-1944. Trowbridge, UK: Association of Track and Field Statisticians.
Lawson, G. (
1997
).
World Record Breakers in Track and Field Athletics
. Champaign, IL: Human Kinetics.
Nevill, A. M. and Whyte, G. (
2005
). Are there limits to running world records?
Med. Sci. Sports Exerc.
37
,
1785
-1788.
Pritchard, W. G. and Pritchard, J. K. (
1994
). Mathematical models of running.
Am. Sci.
82
,
546
-553.
Quercetani, R. L. and Magnusson, R. (
1988
).
Track and Field Performances Through the Years
. Vol.
1
,
1929
-1936. Trowbridge, UK:Association of Track and Field Statisticians.
Sparling, P. B., O'Donnell, E. M. and Snow, T. K.(
1998
). The gender difference in distance running performance has plateaued: an analysis of world rankings from 1980 to 1996.
Med. Sci. Sports Exerc.
30
,
1725
-1729.
Tatem, A. J., Guerra, C. A., Atkinson, P. M. and Hay, S. I.(
2004
). Momentous sprint at the 2156 Olympics?
Nature
432
,
525
.
Tibshirani, R. (
1997
). Who is the fastest man in the world?
Am. Stat.
51
,
106
-111.
Vogel, G. (
2004
). A race to the starting line.
Science
305
,
632
-635.
Weyand, P. G. and Davis, J. A. (
2005
). Running performance has a structural basis.
J. Exp. Biol.
208
,
2625
-2631.
Weyand, P. G., Cureton, K., Conley, D. S., Sloniger, M. A. and Liu, Y. L. (
1994
). Peak oxygen deficit predicts sprint and middle-distance performance.
Med. Sci. Sports Exerc.
26
,
1174
-1180.
Weyand, P. G., Lee, C. S., Martinez-Ruiz, R., Bundle, M. W.,Bellizzi, M. J. and Wright, S. (
1999
). High-speed running performance is largely unaffected by hypoxic reductions in aerobic power.
J. Appl. Physiol.
68
,
2059
-2064.
Weyand, P. G., Sternlight, D. B., Bellizzi, M. J. and Wright,S. (
2000
). Faster top running speeds are achieved with greater ground forces not more rapid leg movements.
J. Appl. Physiol.
89
,
1991
-1999.
Whipp, B. J. and Ward, S. A. (
1992
). Will women soon outrun men?
Nature
355
,
25
.
Yang, M. C. K. (
1975
). On the distribution of the inter-record times in an increasing population.
J. Appl. Probab.
12
,
148
-154.