Introduction

Energy and life history traits are connected. Food energy, once processed by digestion, fuels biomass production, enabling growth and reproduction. The energy budgets of most species are constrained primarily by body size because body size sets fundamental limits on how much energy an organism can acquire and physiologically process1,2. This in turn sets a fundamental limit on the amount of energy that can be allocated to different components of the life history. Understanding variation in life histories requires understanding how constraints on the use and consumption of energy affect demographic traits. Figure 1 shows the relationship between two measures of biological energy use, basal metabolic rate and field metabolic rate and body mass, shown with per-capita industrial energy use across several nations for comparison3,4. Basal metabolic rate is the minimum energy required to support bodily functions in an inactive state, while field metabolic rate captures the energy demands of active individuals under natural environmental conditions. Both exhibit strong statistical relationships with body mass2,3,4.

Figure 1
figure 1

Three kinds of energy use as a function of body mass.

Black circles are basal metabolic rate for mammals4. Grey circles are field metabolic rate for mammals3 The open blue circles are per capita energy use among industrial human populations36 with a star at the median value. The arrows illustrate the linking of per-capita energy use to a corresponding mass by way of following the regression between metabolic rate and body size for mammals. Body mass for all human populations is set to 50kg.

Humans in industrialized nations have a very different relationship with energy, one that is independent of body mass. Energy use among the lowest energy-consuming countries is barely above the expected field metabolic rate of a mammal of our size (Figure 1). In contrast, the median level of energy use among industrial populations (i.e., human populations living in contemporary industrialized nations) is about 1,580 watts, indicated by the star in Figure 1, about 10 times greater than expected for a mammal of our size. Individuals living in countries near the top of the energy use spectrum would have basal metabolic rates equivalent to mammals weighing over 92,500 kg, about twice the size of an adult male Sperm Whale (Figure 1). The magnitude and pace of the increase in the human energy budget is likely unprecedented amongst biological species, in part because it occurred without an evolutionary increase to larger body mass.

Extra-metabolic energy can be used for many purposes that change demographic rates in ways that entirely bypass an organism's internal physiology. Most organisms must burn calories derived from food resources in order to lower their mortality rates by, for example, increasing immune function, working harder to stay warm, or avoiding predators. Industrial human environments provide the opportunity to alter mortality rates (and other demographic traits) via the distribution of goods and services that are largely independent of the calories obtained by digesting food. Inputs of extra-metabolic energy can alter demographic patterns both directly, via birth control or medication, for example, or indirectly by improving infrastructure, such as water quality, roads and schools5,6. In spite of wide interest in resource constraints on the growth and structure of populations, little is known about the basic statistical responses of human demographic and life history traits to variation in energy use. Understanding such responses is necessary for developing a complete theoretical picture of human demographic processes.

Here we analyze variability in key life history traits across a spectrum of populations that vary widely in per capita energy use. Specifically, we quantify the relationship between human life history traits and extra-metabolic energy use and ask whether different traits vary in the strength of their response. We take energy as the independent variable in order to test the hypothesis that increases in energy use, independent of body size, influence life history traits consistently, regardless of whether energy is used internally or externally (metabolically or extra-metabolically, respectively). The resulting relationships yield insights into the demography of contemporary human populations, including the perplexing issue of below-replacement fertility rates observed among many wealthy nations today7. Furthermore, we show that many of the life history allocation patterns of industrial human populations are simply extensions of the ecological energetics of biological populations.

We investigate human demographic traits using an allometric approach. To introduce this approach we first describe conventional allometry and its role in life history studies, which are generally based on body mass8. We convert these relationships from mass to energy. Allometric relationships are power functions of the form y = c0mz where y is some dependent variable of interest that scales with body mass m to the exponent z and c0 is a normalization constant. Functions of this form become linear on logarithmic axes, where z is the slope and log(c0) the intercept. The value of the slope captures the response of a change in y due to a change in m and when the coefficients are fit to data one says that y scales with mass to the z power.

Essentially all life history traits can be characterized as allometric functions of mass8, making body mass perhaps the most fundamental attribute of a species9. Biological rates such as birth rate and mortality rate tend to scale negatively with mass as z close to −1/4, while biological durations such as the time to maturity and the total lifespan tend to scale positively as z about +1/48,9,10,11. In the following analysis we are concerned primarily with the direction (sign) and steepness of the scaling exponent (slope) z, rather than the exact value. Body size is the key constraint on life history traits because, by setting the field metabolic rate, it strongly influences the total amount of energy that can be acquired and allocated among life history traits1 (Figure 1). From studies of the role of body size in life history it is clear that there are consistent selective pressures regulating how an organism's size-specific metabolic energy budget influences tradeoffs among life history traits. Figure 1 illustrates how fundamental this size-based constraint is across mammals3,4, as well as how different the industrial human energy budget is compared to energy budgets that are constrained primarily by mass. Because of the strong correlation between size and metabolic rate, any trait that scales with body size can also be expressed as a function of metabolic rate. By rewriting allometries in terms of metabolic rate rather than mass, we can examine the change in life history traits in response to changes in energy use. Considering energy instead of mass makes it possible to compare the vital rates of industrial human populations to metabolically constrained species of mammals and primates.

We take a macroscopic approach to assessing how life history traits respond to energy use in order to answer a general question about whether the source of that energy matters. To develop null predictions for the relationships between life history traits and per capita energy use we employ a combination of existing cross-species allometries and the mammal life history model developed in Charnov1,10, which predicts the slopes of many life history traits as functions of mass and after our extension, energy, that are grounded in explicit evolutionary tradeoffs. Given known patterns in cross-species allometries and Charnov's mammal life history model, the basic null predictions are that 1) rates of mortality, fertility and population growth should be negatively related to energy use; and 2) the time to a life history event, including to sexual maturity, first reproduction and lifespan, should be positively related to energy use.

Results

Birthrate and Juvenile Mortality

Moses and Brown12 established the approach of energy-based allometries by demonstrating that birthrate scales with energy use as an extension of the body mass relationships observed across mammals (Figure 2A). They showed that the scaling of birthrate with mass to the −1/4 power can be re-expressed as a scaling of birthrate with energy use to the −1/3 power (i.e., b M−1/4 and B M3/4 so b B−1/3 where b is birthrate, M is mass and B is metabolic rate). The predicted −1/3 power scaling between birthrate and energy consumption was found not only across biological species but also among industrialized human populations. Further, Moses and Brown (2003) demonstrated the remarkable robustness of this relationship for many different years across countries and over time within the United States between 1850 and the present.

Figure 2
figure 2

Life history traits as functions of per capita energy use.

(a) Birthrate. The slopes linking birthrate to per capita energy use are similar for mammals and industrial populations and both are near the predicted value of −1/312. (b) Mortality under age five is a steeply negative function of per capita energy use (slope near −1) among industrial nations. (c) Lifespan. Curved solid line is an asymptotic fit to the industrial population data: log(lifespan) = 4.56+(4.84)*exp(−0.36*log(energy)). The dotted line is the linear fit to the nonhuman primate data. Note that the primate data are maximum lifespans whereas the industrial nations are life expectancy at birth. Exactly comparable data for primates would cause the primate data points to shift downward noticeably. (d) Age at first reproduction for industrial nations, nonhuman primates and an average hunter-gatherer.

Figure 2B indicates that juvenile mortality, measured here as mortality under age five, is a steeply negative function of energy use with an exponent close to −1 (Table 1). Interestingly, this suggests that rates of juvenile mortality drop proportionally with increasing energy consumption, an exceptionally steep benefit to energy use, as twice the energy consumption leads to half the mortality, which contrasts markedly with adult lifespan (see below).

Table 1 Summary statistics for OLS regressions.

Adult Lifespan and its Divisions

While birthrate and under five mortality exhibit strong responses to per-capita energy use, adult lifespan exhibits an interesting curvature (Figure 2C). Initially, lifespan rises steeply with per capita energy use and then flattens toward an asymptote at high values13. This pattern follows the qualitative prediction in that the slope is positive, but the curvature deviates from the basic log-linear form of an allometric relationship resulting in a slope shallower than the theoretically-predicted 1/3. The observed energy-based allometry for lifespan also differs from the empirical trend across primate species. Primates generally deviate from the overall mammal life history in having slow growth, slow reproduction and long lifespans14,15. But industrial human lifespans respond to energy differently than either primates or non-primate mammals. The regression line for the primate data suggests that lifespans among industrial humans are actually shorter than expected based on simply extending the cross-primate allometry to higher levels of energy use. In other words, humans have shorter lifespans than a hypothetical primate species with biological metabolic rates equivalent to the current high levels of per capita energy use (dotted line, Figure 2C).

Age at first reproduction also scales positively, but with a shallower slope than the predicted +1/3. The scaling of age at first reproduction for humans of 0.1, although positive as predicted, is shallower than the scaling for primates, 0.45 (Table 1).

Figure 3A examines average lifespan alongside some of its major components. All of the pictured life history traits relevant to the divisions of the adult life-course (such as ages at menarche, first birth and menopause) are significant allometric functions of per capita energy use, but they are not all directionally consistent with the null life history predictions (Table 1). For example, while age at menopause scales positively with energy use, the shallow slope of 0.02 (95% CI: 0.012 – 0.038) demonstrates that age at menopause is less responsive to energy consumption than other life history traits. Because the age at first reproduction increases more steeply with energy use than age at menopause, the length of the reproductive lifespan, the time between age at first reproduction and the age at menopause, decreases with energy use (Table 1). As most organisms reproduce for the majority of their adult lives and lifespan should scale positively with body mass and energy, the negative scaling of reproductive lifespan with energy is contrary to the theoretical expectation and suggests an important constraint on the life history that emerges as a function of increasing energy inputs.

Figure 3
figure 3

Life history traits as functions of per capita energy use.

(a) The lifespan and its major divisions for industrial populations (nations). (b) Reproductive delay (age at menarche - age at first reproduction) varies similarly with per capita energy use among nonhuman primates and industrial human populations. The dashed line is a linear fit to the combined datasets, nonhuman primates and industrial nations. (c) The net reproductive rate for industrial populations. Horizontal dotted line drawn at log(R0) = 0. Growing populations are above the line and declining populations are below it (not accounting for migration or other dynamics). (d) Probability of survival to age five for industrial populations. Median value of energy use indicated with the dashed line.

Reproductive delay is the period of time between menarche and age at first reproduction. The relationships in Figure 3A indicate that this period of time is an increasing function of energy use among industrial human populations. Plotting reproductive delay across primates to facilitate comparison yields a somewhat surprising result (Figure 3B); reproductive delay increases with body size with an exponent near 1/4 and with metabolic rate to an exponent near 1/3. Interestingly, the slope for reproductive delay as a function of energy use among industrial populations is also near 1/3, almost identical to the primate value (Table 1). While the strength of the relationship is somewhat weak in primates, it is highly significant and the low variance could be attributed to the difficulty of observing the parameters on actual populations. The possibility that reproductive delay may be a general 1/4 power function of body mass across species as well as an underappreciated allometric measure for the pace of the life history warrants further attention.

The Net Reproductive Rate and Below Replacement Growth

Vital rates such as birthrate and mortality do not exist in isolation but are components of aggregate metrics, such as the net reproductive rate. Given the findings above, we may also wonder how the net reproductive rate responds to changes in per capita energy use. The net reproductive rate, R0, combines birthrate and survival (lifespan) into a single quantity in units of daughters per mother16. Non-growing populations have, by definition, lifetime reproductive rates such that each mother exactly replaces herself; i.e., R0 = 1 at replacement. If R0 is greater than 1 a population is growing and if less than 1 it is declining.

Across industrial populations R0 drops with per capita energy use (Figure 3C), consistent with the commonly noted trend for mortality and fertility to be lower among higher energy consuming nations or social strata17,18. Figure 3C shows that R0 declines with per capita energy use with a slope near −1/3 and that 51 of the 107 (48%) nations represented are at or below replacement levels of growth in terms of R0. If these were biological populations in the wild, we would predict that these populations were headed toward extinction (ignoring immigration). Because the component parts of R0 are energy-based allometries, our results shed light on why net reproductive rates frequently drop to such exceptionally low levels in human populations.

R0 is the product of birthrate b, lifespan E and the probability of survival to reproductive age S, so R0 = S * b * E10. It is well known that b and E are allometric functions of mass M, E = a0M1/4, b = b0M−1/4 8, where a0 and b0 are constants. Figures 2A and 2C demonstrate that they are functions of energy as well. As a result, R0 can be written as a function of mass:

Or equivalently, as a function of energy:

where B is either metabolic rate, or extra-metabolic energy use12,19. In equation 1, the exponents of 1/4 for lifespan and −1/4 for birthrate cancel and so R0 must equal one for any non-growing population regardless of mass. However, equation 2 provides a means of evaluating the effects of per capita energy use on R0 via the energy-based allometries for b and E, by writing it generally as a function of energy with unknown coefficients:

Because the balance between birth and mortality (inverse lifespan) drives changes in R0 we define δ = a3 * a4 and z = z1 * z2. Taking the logarithm of both sides yields

where the balance between fertility and mortality is captured by the value of z. Otherwise typical biological populations will have z = 0 because the empirical relationships between B and the vital rates would cancel each other out, as in equations 1 and 2. If S is held constant in equation 4 then it is apparent that if either b or E scaled differently than predicted (i.e., other than +/− 1/3) across populations, then R0 would vary with energy use. This means that z captures the balance between reproduction and lifespan (avoiding mortality) and demonstrates how population growth rates might be affected by varying allometric constraints on b and E. Holding S constant here is quite reasonable because S varies only slightly as a function of per capita energy use for high values, those above the median and in fact varies relatively little across the entire dataset (Figure 3D). As such, juvenile survival varies little among wealthier nations, which are generally those that have either gone through, or are going through the demographic transition to below-replacement fertility.

To investigate the effect of z on R0 we examine the derivative of log(R0) with respect to B (or M):

or as a function of log(B):

When z is negative the net reproductive rate declines and any species with a negative z should have fewer offspring per reproductive lifetime. The relationships in Figures 2A and 2C can be used to test this prediction. A linear fit to the lifespan relationship in Figure 2C gives a slope of 0.12, which can be combined with the slope for birthrate to give , which is just within the 95% confidence intervals for slope of log(R) vs. log(B) in Figure 3C (95% CI: −0.358 to −0.246), but this calculation underestimates z because it does not account for the curvature in the relationship between lifespan and per-capita energy use. If we subset the data to values of energy use above the asymptote in the lifespan data the slope falls more squarely within intervals. Nonetheless, it is clear that lifespan responds to energy-based allometric constraints differently and less plastically than does birthrate. Thus the imbalance between the energetic dependence of birth and death rates may be related to the drop in net reproductive rate to below replacement levels.

Discussion

Our results show that many human life history traits can be described as allometric functions of per capita energy use. In some cases these allometries are broadly consistent with the variation in metabolic rates across species, but in others the relationship differs in direction or shape. The ability to harvest large sums of extra-metabolic energy from the environment raises the total per capita energy budget to levels far beyond the field metabolic rate predicted for a primate of our size (Figure 1). This increase changes the constraints on human life history traits. Some life history traits respond to energy with steeper allometric exponents than others (Table 1), which indicates that some traits respond more readily to extra-metabolic investments of infrastructure, goods and services, while others are more closely tied to physiological constraints and hence are not as easily adjusted by technological development. The differences in response and their underlying constraints, can drive high-energy populations to below-replacement levels of net reproductive rate, which in turn suggests a route to a novel energy-based explanation for the phenomenon widely referred to as the demographic transition17,20.

The traits that are consistent with the null predictions, in that they scale with energy with the predicted sign (direction), are birthrate, juvenile mortality, reproductive delay, age at first reproduction and age at menopause. Traits that respond in the opposite directions and thus run contrary to a life history or allometric expectation, are the length of the reproductive lifespan and the age at menarche (Table 1). Birthrate is predicted especially well, quantitatively matching the slope for mammalian patterns12, adding to the idea that processes of industrial development, like urbanization, may have particularly consistent influences on offspring investment tradeoffs21. Lifespan increases, which is consistent with the direction of the cross-species pattern, but differs in shape from the standard log-linear allometries by bending toward an asymptote at high levels of per capita energy use13.

These traits can also be described as varying along a continuum with respect to how steeply each responds to per capita energy use (Table 1). Juvenile mortality is the most plastic, followed by birthrate, reproductive delay, age at first reproduction, lifespan and age at menopause is the least plastic. Not surprisingly, early life mortality and late life mortality are on opposite ends of this spectrum; mortality at young ages can be more efficiently reduced than mortality at late ages. Age at menopause and life expectancy exhibit the least-plastic responses to energy use of any of the traits examined. Both of these may be more closely tied to physiological processes than the traits that exhibit highly plastic (steep) responses to changes in energy use. For example, extra-metabolic energy can have a very strong effect on lowering infant mortality but menopause responds with a very shallow slope, perhaps reflecting the hard physiological constraint of follicle depletion that leads to the onset of menopause22. Menopause and lifespan are also influenced by extra-metabolic energy, but their shallower slopes demonstrate weaker responses and likewise imply that biology and actual physiological metabolism constrain these traits more than others. Moreover, this suite of statistical responses implies that the drivers of these patterns are closely related to the energetic properties of human environments, primarily extrinsic variables and not uniquely based on social or cultural phenomena that might be posited to exist independent of energy use.

Physiological constraints are the most likely cause of the asymptotic relationship between lifespan and energy use. While average lifespans can be extended using extra-metabolic energy (by improved sanitation and increased access to healthcare, for example), the biophysical requirement of having to metabolize energy for somatic maintenance sets a constraint that cannot be indefinitely stretched by inputs of extra-metabolic energy. Thus far, however, human populations have continued to lower mortality in all age groups and the limits of this trend have been notably difficult to estimate23,24. Nonetheless, while industrialized humans are exceptionally long-lived for mammals23, they are not as long-lived as predicted for a primate with an equivalent body size given our extra-metabolic energy consumption. That is, a hypothetical primate species following primate allometric patterns but with the metabolic rate equivalent to the energy use of a wealthy industrial population would be very large and have a very slow birth rate, but would also produce offspring for longer periods of time, keeping the net reproductive rate near replacement. Without the allometrically appropriate mass and physiology, however, contemporary human life histories face novel constraints that can drive the net reproductive rate to below replacement levels.

The actual length of the reproductive period decreases with increasing energy use, implying that the factors accompanying increased energy consumption may increase, rather than decrease, constraints on the length of the reproductive lifespan. For example, Korperlainen20 showed that reproductive periods reduced to very short intervals (less than five years) accompanying socio-economic development in Finland between 1870 and 1949. A more constrained, shortened, reproductive lifespan must affect tradeoffs between reproduction and investments in the accumulation of wealth and or/or status, or “maintenance costs” in life history terms, because the amount of time during which offspring can be produced directly constrains the number that can be produced in a lifetime25. Therefore, the shortening of the reproductive lifespan is likely an important part of the story leading to below replacement reproduction. While not specifically examined, Figure 3A also suggests that the maximum potential reproductive lifespan is increasing with energy use as the time between menarche and age at menopause seems to increase. It seems that maintaining a high-quality industrial niche requires that increasingly large fractions of the potential reproductive period be converted to forms of growth and maintenance (economic productivity) that compete with or constrain maximum potential reproductive output21,26.

While human energy budgets have increased dramatically with the advent of the industrial age, it is not necessarily the case that subsidizing the life history budget with additional resources occurs only in industrial contexts. Many anthropologists have noted that the evolved human life history pattern consists of costly and dependent offspring. Among hunter-gatherers it is clear that to raise an offspring from weaning to adulthood requires substantial subsidies from the foraging efforts of individuals other than the parents, including unrelated males, older siblings and grandmothers27,28,29,30. Indeed, a more remarkable feature of the evolved human life history pattern is that the relative cost of offspring seems to assume the presence of a social network that will increase food supplies for the offspring beyond that which the parents would be capable of supplying on their own27,28. As such, the contemporary industrial energy budget should be seen as the most recent manifestation of an evolving human energy budget, an amplification of a pattern that can be traced back to our hunter-gatherer roots. In this light the industrial energy budget could be seen as increasing the energy budget by large quantities, but it would not be strictly novel to industrial peoples (although the form of the energy is of course quite different). If so, the fact that some scaling patterns are preserved across such complex fuel sources and distribution networks is especially interesting and implies some emergent structure and/or deep evolutionary roots to how human societies allocate resources.

A note about causality

Our analysis begins with a factual premise. Humans have increased their energy budgets with the addition of energy from industrial sources. As with all organisms, our physiology and life histories evolved with an energy budget set primarily by body size. We find that demographic traits respond to this extra-metabolic energy in a variety of ways. We do not investigate the various possible proximate mechanisms for each trait we consider. There are numerous and varied routes by which energy may link goods, services and infrastructure to health, demographic rates and economic investment strategies. However, a necessary step preceding the explanation of any macroecological pattern is its empirical recognition31. We have introduced several statistically significant allometric relationships that beg explanation and our central objective was simply to ask whether such patterns exist and to identify their basic characteristics. Our analysis highlights the need to understand the ways that energy can limit the mechanisms generating these relationships. As energy is ultimately responsible for all biochemical work, it must play a causal role in the processes linking societal infrastructure to the growth and maintenance of the humans depending on the resources these infrastructures provide.

In sum, all of the life history traits we examined are significant functions of per capita energy use. Many of them vary in ways similar to the variation of metabolic rate across species. The balance of these traits provides important insights for the relationship of human demographics to the exceptionally large industrial energy supply. The patterns examined here take initial steps toward understanding how extra-metabolic energy use impacts human life history traits. While humans have a life history mediated through technology, energy and myriad social factors, we demonstrate that there is a consistent structure to the demographic patterns of industrialized populations, at least at the large aggregated scales examined here. On the one hand, given the flexibility and rapid changes of technology and energy, combined with the complex social interactions of modern societies, it may seem surprising that allometric relationships describe the life history traits of industrial humans at all. On the other hand, given that energy is fundamental to all ecological and evolutionary processes, including population dynamics and life histories, perhaps we should not be surprised to find that it affects human demographics as well.

Methods

Here we present the sources and relevant calculations for the variables depicted in the Figures. Basal and field metabolic rates (BMR and FMR) for mammals were obtained from the two most recent compilations we are aware of4,32. BMR was originally reported in units of milliliters of oxygen consumption per hour and FMR in units of kilojoules per day, each of which we converted into watts (joules per second). Per capita energy use for industrial populations is from the Earthtrends database at the World Resource Institute (WRI). The WRI provides methodological details, the most germane of which we repeat here: “Total energy consumption per capita measures the amount of primary energy consumed, on average, by each person living in a particular country or region for the year indicated. All primary sources of energy, including coal and coal products, oil and petroleum products, natural gas, nuclear, hydroelectric, etc., are included here. Data are reported in kilograms of oil equivalent (kgoe) per person. Consumption equals indigenous production + imports - exports - energy delivered to international marine bunkers +/− stock changes…Energy losses from transportation, friction, heat and other inefficiencies are included in these totals.” These data were also converted into Watts12.

A limitation of this measure of energy use is that we assume that the average level of energy consumption in an aggregated and coarse dataset is a proxy for variation in individual-level access to energy across countries. Within-country variation in access to energy is not included and this has to be an important source of variation in the patterns we depict here. An open and interesting question is the degree to which these patterns would be similar across social strata within countries given the asymmetry in the distribution of individual access to energy sources.

To derive allometric relationships linking body mass to metabolic rate in Watts for mammals and primates we use empirical fits to the data in3 and4. The final estimates were: FMR = 8.89*Mass0.734 for mammals and FMR = 7.56*Mass0.765 for primates (mass in kg). The primate intercept had to be estimated based on the difference between the heights of the functions for FMR and BMR in mammals (because of the strongly parallel nature of the fits of BMR and FMR to mass, Figure 1).

All non-human primate data are from33. Birthrate and body size data for mammals are from12; for hunter-gatherers from27; and for industrial populations from WRI. Annual birthrate for industrial populations is the total fertility rate divided by the difference of the best fit allometric functions between age at menopause and age at first reproduction vs. per capita energy use, with the difference between the first reproduction and menopause taken as the length of the reproductive lifespan. Moses and Brown use an average reproductive lifespan of 20 years and get the same basic result12. Mortality under age five and life expectancy data for industrial populations are from WRI. Because mortality at ages under five is available in units of deaths per 1000, the corresponding probability of survival was calculated as (1000-number of deaths under age five)/1000 (Figure 3C). Age at first reproduction for industrial populations were obtained as Selected Indicators from the United Nations Population Division34. Ages at menarche and menopause for industrial humans are from35. The reproductive lifespan was crudely estimated as the time between first birth and menopause in human populations and the time between first birth and lifespan in nonhuman primates. Both are certainly overestimates and the values in Table 1 should be seen as only rough approximations for this variable. For the reproductive delay regression we excluded negative values. Net reproductive rate data are national averages from 2000 to 2005 and were obtained from the UN website34.

All variables were natural-log-transformed before analysis, except for length of reproductive lifespan and best-fit lines obtained by ordinary least squares regression. The asymptotic fit applied to the lifespan data was implemented in R version 2.9.1. Outliers are removed from the industrial populations (all of which are oil producing nations), but including them does not appreciably alter the results and their relevance is discussed in Moses and Brown12 and DeLong et al.19.