Productivity, biodiversity, and pathogens influence the global hunter-gatherer population density

Miikka Tallavaara, Jussi T. Eronen, and Miska Luoto
PNAS February 6, 2018 115 (6) 1232-1237; published ahead of print December 27, 2017 https://doi.org/10.1073/pnas.1715638115

  1. Edited by Robert L. Kelly, University of Wyoming, Laramie, WY, and accepted by Editorial Board Member Richard G. Klein December 6, 2017 (received for review September 4, 2017)

Significance

Because of complex cumulative culture, human populations are often considered to be divorced from the environment and not be under the same ecological forcing as other species. However, this study shows that key environmental parameters net primary productivity, biodiversity, and environmental pathogen stress have strong influence on the global pattern of hunter-gatherer population density. Productivity and biodiversity exert the strongest influence in high and midlatitudes, whereas pathogens become more important in tropics. The most suitable conditions for preagricultural humans are found in temperate and subtropical biomes. Our results show that cultural evolution has not freed human hunter-gatherers from strong biotic and abiotic forcing.

Abstract

The environmental drivers of species distributions and abundances are at the core of ecological research. However, the effects of these drivers on human abundance are not well-known. Here, we report how net primary productivity, biodiversity, and pathogen stress affect human population density using global ethnographic hunter-gatherer data. Our results show that productivity has significant effects on population density globally. The most important direct drivers, however, depend on environmental conditions: biodiversity influences population density exclusively in low-productivity regions, whereas pathogen stress does so in high-productivity regions. Our results also indicate that subtropical and temperate forest biomes provide the highest carrying capacity for hunter-gatherer populations. These findings document that environmental factors play a key role in shaping global population density patterns of preagricultural humans.

Questions related to the environmental drivers of species distributions and abundances have been at the core of ecology throughout its history. However, these questions have rarely been addressed for preindustrial human populations, despite the fact that these populations provide a rare opportunity to investigate drivers of the distribution and abundance of a single species along global environmental gradients. Recent archaeological studies have correlated long-term human population dynamics and major dispersal and range shift events with climate changes (1⇓⇓⇓–5), but their results are mixed regarding the significance of climate in shaping human population processes. Ethnographic data provide finer-scale information about preindustrial human populations and their environments, potentially allowing more conclusive analyses. However, even for hunter-gatherers, ethnographic studies about the ecological determinants of population patterns are relatively few (6⇓⇓–9), despite hunter-gatherers’ close dependence on their immediate environment and wild resources.

Here, we use global ethnographic hunter-gatherer data (10, 11) to explore the effects of key environmental variables on human population densities. These variables are net primary productivity (NPP) calculated as a function of the annual mean temperature and annual precipitation; biodiversity as surrogated by a combination of mammal, bird, and vascular plant richness; and environmental pathogen stress as a combination of pathogen richness and the severity of exposure (Materials and Methods and Fig. S1). In previous work, only the effects of productivity on human population densities have been assessed. These studies suggest that both primary and secondary productivity have, at least regionally, positive effects on hunter-gatherer population density as well as on population home ranges (6⇓⇓–9). Such positive effects are expected, because hunter-gatherers access food directly from their surroundings, which vary widely in energy availability (12). Although hunter-gatherers depend on the productivity of wild plant and animal species, they appropriate only a small fraction of the production and in this regard, more resemble other land mammals of their size than they do humans of modern industrialized societies (13).

However, the amount of available food may not be all that matters in controlling hunter-gatherer abundance, as biodiversity can also play a role. Biodiversity influences ecosystem stability, such that with increasing levels of diversity, the temporal variance of ecosystem productivity tends to decrease and the temporal mean of ecosystem productivity tends to increase (14, 15). A current consensus is that higher biodiversity can enhance the temporal stability of aggregate ecosystem properties, such as biomass and productivity, but its effect on stability of individual populations remains controversial (15, 16). The stabilizing effect of biodiversity may come through processes like asynchrony in the responses of species to environmental variation (17). For hunter-gatherers, increased stability of ecosystem-level biomass production caused by higher biodiversity decreases subsistence-related risk, which can positively affect hunter-gatherer population densities.

In contrast to the effects of productivity and biodiversity, the effects of pathogens on hunter-gatherer abundance are negative. A well-known example is the demographic shock among Native Americans caused by infectious diseases after European contact (18, 19). Beyond the effects of such dramatic epidemiological events, we postulate that the background level of pathogen stress specific to each environment also influences human population density. Studies in animal ecology have shown that pathogens can have negative effects on hosts’ abundance (20⇓–22), and we expect that this applies to human populations as well.

Previously, Binford (10) has suggested that resource abundance and stability as well as pathogens would have played an important role in constraining hunter-gatherer population growth. However, Binford (10) does not analyze this hypothesis with the data, whereas our study presents a test of the joint effect of NPP, biodiversity, and pathogens on hunter-gatherer populations.

We use structural equation modeling to reveal the potentially hierarchical effects of NPP, biodiversity, and environmental pathogen stress (23) (Materials and Methods). This method allows us to explore the network of variables and thus, the indirect effects of the explanatory variables. In addition to its direct effect on hunter-gatherer population density, NPP may have an indirect effect mediated by biodiversity. Although the actual mechanism is still debated, plant and animal richness is related to productivity, and these relationships are predominantly positive and monotonic (24, 25). Furthermore, pathogen richness follows the same latitudinal pattern as biodiversity in general, and this pattern is most likely related to temperature and precipitation (26). Thus, climate-based NPP may also have an indirect effect on population density through pathogen stress. However, human pathogen richness is also related to the richness of alternative hosts (27). Consequently, the effects of climate-based NPP on population density can be simultaneously mediated by both biodiversity and pathogen stress. This hypothesizedx network of variables is illustrated in Fig. 1.

Fig. 1.
Fig. 1.

Structural equation metamodel showing hypothesized expectations about the relationships among NPP, biodiversity, environmental pathogen stress, and hunter-gatherer population density. Blue arrows indicate positive effects, and red arrows indicate negative effects.

In addition, we provide a geographical interpretation of environmental effects on population density by modeling a distribution map of the relative contributions of the limiting effects of NPP, biodiversity, and pathogen stress on global hunter-gatherer population density. We also model a potential distribution of global hunter-gatherer population density itself as a function of the three environmental variables. These maps are based on the combination of two regression modeling techniques and are projections to current environmental conditions assuming that the globe is populated by hunter-gatherers only (Materials and Methods).

Results

Visual inspections and breakpoint analysis of the bivariate relationship between climate-based NPP and hunter-gatherer population density indicate an NPP threshold of ∼1,360 g/m2 per year, above which the population density becomes largely independent of additional increases in annual productivity (Fig. 2A). However, above the same threshold, NPP exhibits a strong predictive power of pathogen stress (Fig. 2B). This raises the question as to whether the leveling off of hunter-gatherer population density in high-productivity environments is related to increased pathogen stress.

Fig. 2.
Fig. 2.

Contrasting effects of NPP on hunter-gatherer population density and pathogen stress. (A) Relationship between NPP and hunter-gatherer population density. Bubble size reflects the magnitude of pathogen stress. A bivariate piecewise regression model (black curve) explains 44% of the variation in population density. (B) Relationship between NPP and pathogen stress. Bubble size reflects hunter-gatherer population density. A bivariate piecewise regression model (black curve) explains 29% of the variation in pathogen stress. A and B also show linear model breakpoints, their SEs, and the mean of these breakpoints (1,360 g/m2 per year), which is used to define the low- (n = 234) and high-productivity (n = 66) environments. In addition, A and B show coefficients of NPP and their P values in regressions between NPP, population density, and pathogen stress in low- and high-productivity environments.

To further analyze these potentially complex relationships between variables, we used structural equation modeling. To ease interpretability of the direct and indirect effects, we built separate models with identical structure (Fig. 1) for the low- and high-productivity environments delineated by the NPP threshold identified in the breakpoint analysis.

The structural equation modeling results strongly indicate different constraints between the low- and high-productivity environments. In the low-productivity environments, biodiversity has the strongest direct effect on hunter-gatherer population density followed by NPP, whereas the effect of pathogen stress is not statistically significant (Fig. 3A, Fig. S2A, and Table S1). In the high-productivity regions, however, the importance of the explanatory variables is reversed. Here, pathogen stress has the strongest effect, and NPP still has relatively high direct and indirect effects, but biodiversity is no longer significant (Fig. 3B, Fig. S2B, and Table S1).

Fig. 3.
Fig. 3.

Relationships between variables. (A) Structural equation modeling of low-productivity environments (NPP ≤ 1,360 g/m2 per year). (B) Structural equation modeling of high-productivity environments (NPP > 1,360 g/m2 per year). In A and B, NPP is considered an exogenous variable. Arrow widths and accompanying numbers indicate the relative effects (standardized regression coefficients) of different variables. Unstandardized coefficients are shown in parentheses. Gray arrows indicate statistically nonsignificant relationships, and the red arrow indicates a negative effect. Effects are based on the component models including only statistically significant relationships (α = 0.05). Regression planes show how population density responses to the joint effects of NPP and biodiversity (C) and of NPP and pathogen stress (D) in the linear model fit to the whole data (i.e., data are not split into low- and high-productivity environments as in the structural equation modeling). These model results correspond to the structural equation modeling results about direct effects of biodiversity and pathogens. In C and D, red lines indicate the NPP threshold separating low- and high-productivity environments in the structural equation modeling.

Because of its indirect effect through biodiversity, NPP has the strongest total effect (combined direct and indirect effects) on population density in low-productivity regions (Table S1). In high-productivity regions, the total effect of NPP is weaker than the direct effect of pathogens (Table S1).

As a cross-check to the structural equation modeling results regarding the direct controls of hunter-gatherer population density, we fit a linear model to the whole data (i.e., the data were not split into low- and high-productivity environments as in the structural equation modeling). In this model, we included interaction terms between biodiversity and NPP and between pathogen stress and NPP. The model shows the same pattern of direct effects of biodiversity and pathogen stress on population density as the structural equation modeling (Fig. 3 C and D, Fig. S3, and Table S2): the effect of biodiversity decreases as NPP increases, and the effect of pathogens becomes stronger as NPP increases. This strongly supports our structural equation modeling results regarding the different controls of population density in low- and high-productivity environments.

The modeled geographical distribution of the relative limiting effects of NPP, biodiversity, and pathogens on hunter-gatherer population density suggests that, among the three variables, environmental pathogen stress has the strongest limiting effect on population density in the tropics and subtropics (Fig. 4 A and E). In contrast, NPP and biodiversity constrain hunter-gatherer population density in mid- and high latitudes (Fig. 4 A, C, and D). These two factors tend to be colimiting, especially in high latitudes. NPP assumes a more dominant role in North America and in northernmost Eurasia, whereas biodiversity assumes a more dominant role in dry regions and in midlatitude Eurasia. (Fig. 4 A, C, and D).

Fig. 4.
Fig. 4.

Global projection of hunter-gatherer population density and environmental constraints thereof under current conditions. (A) Geographical distribution of the limiting effects of NPP, biodiversity, and pathogen stress on population density. Red, green, and blue indicate dominant limiting effects of pathogens, NPP, and biodiversity, respectively. The higher the limiting effect of a variable, the farther away that it is from its optimal value. For example, pathogen stress is limiting in the Mediterranean region, although the pathogen stress there is lower than in sub-Saharan Africa. This is because the NPP and biodiversity are close to their optimum values in the Mediterranean region, leaving the constraining role to pathogens. (B) Map of predicted population densities. (C–F) Relative limiting effects of NPP (C), biodiversity (D), and pathogen stress (E) and predicted population density (F) in different ecological zones. Global ecological zones are according to the FAO (28, 29). In the box-and-whisker plots (C–F), the black horizontal lines represent the medians, the filled circles represent the means, the boxes cover the interquartile ranges, and the whiskers represent the minimums and maximums of the distributions.

Given the effects of NPP, biodiversity, and pathogen stress, the potential global distribution of hunter-gatherers highlights previously known areas of high hunter-gatherer abundance, such as the western coast of North America and the eastern and northeastern coasts of Australia (10, 11) (Fig. 4B and Fig. S4). However, the map also suggests a number of areas with potentially high hunter-gatherer abundances that have been occupied by agricultural populations from at least the late mid-Holocene onward and that are not well-represented in the historical ethnographic data on hunter-gatherers (Fig. 4B and Fig. S4). These areas include southwestern Mexico, southeastern North America, Western Europe, Southeast Asia, and Papua New Guinea. The highest predicted mean population densities are in temperate oceanic forest biome followed by subtropical humid forests, temperate continental forests, and tropical rainforests (Fig. 4F). However, the differences between these highly suitable biomes are small. When interpreting these results, one has to keep in mind that the global projection takes into account only the effects of three variables that reflect the influence of terrestrial environments on hunter-gatherers. Many coastal hunter-gatherer populations have been dependent on aquatic resources, especially in high latitudes, where terrestrial productivity is relatively low (10, 11). Thus, it is possible that the global projection of population density would change if productivity of aquatic environment is taken into account.

Discussion and Conclusion

Our results show that NPP has significant positive effects on hunter-gatherer abundance globally. However, there is no single determinant of the large-scale variation in abundance. Rather, the most influential direct drivers of population density are context-dependent. Biodiversity affects population density mostly in the low-productivity environments, whereas pathogen stress, which has insignificant effects in areas of low productivity, is a crucial constraint on population density in the high-productivity settings, especially the tropics. The insignificant effect of biodiversity on hunter-gatherer population density in high-productivity environments is not totally unexpected, as ecosystem properties tend to have a saturating response to increasing species richness (14, 15). In high-productivity settings, which often coincide with high biodiversity, increasing diversity does not increase ecosystem stability; consequently, the response of hunter-gatherer population density to biodiversity seems to saturate as well.

As such, our results indicate that the factors affecting hunter-gatherer abundance on a global scale are similar to those controlling global-scale biodiversity: macroscale studies on biodiversity indicate that, in high-latitude regions, productivity is limiting diversity (30⇓–32), whereas in the tropics, biotic interactions become a more important driver (33). In the case of hunter-gatherers, abundance and stability of resources limit population density in the low-productivity environments, whereas species interactions, as represented in this study by pathogen stress, become more important in controlling population density in the high-productivity environments. Similarities in the patterns of limiting factors of hunter-gatherer population density and global biodiversity might suggest similarities in the underlying mechanisms as well.

However, because we do not have direct measures of pathogen prevalence among individual study populations, we cannot completely rule out spurious relationships between pathogen stress and hunter-gatherer population density. Due to generally low population densities, it has been suggested that the prevalence of pathogens, such as influenza and measles viruses, and Plasmodium and Trypanosoma parasites is low among hunter-gatherers (34). In addition, it has been argued that high rates of mobility prevent the accumulation of fecal and other waste at hunter-gatherer dwelling sites, leading to a low prevalence of soil-transmitted helminths (35, 36). Even if hunter-gatherers are less susceptible to pathogens than farmers who live under similar conditions of environmental pathogen stress, our findings nevertheless indicate that the interaction between pathogens and hunter-gatherers is sufficient to create a significant negative effect on hunter-gatherer abundance.

By suggesting that subtropical and temperate biomes can sustain the highest hunter-gatherer population densities, our results corroborate the earlier untested hypothesis of Binford (10) regarding the role of resources and pathogens in creating optimal conditions for hunter-gatherers in temperate and subtropical biomes. Based on our results, it seems that, in subtropical and temperate biomes, there is an optimum tradeoff between the opposing effects of resource availability and pathogens: at the higher-productivity end, negative biotic interactions are too limiting, whereas at the lower-productive end, resources are too limiting. However, a large fraction of the area of these highly suitable biomes has been occupied by agricultural populations since at least the late mid-Holocene. It has been argued that stable environments with high resource abundance and richness provide favorable conditions for the innovation and adoption of agriculture (37, 38). Thus, it is understandable that the most suitable environments for hunter-gatherers are also among the first ones to support agriculture. Furthermore, potentially high hunter-gatherer population densities in these areas might have facilitated the innovation, spread, and development of the new subsistence technology. Alternatively, high hunter-gatherer population density itself, resulting in population packing, and reduced options for mobility might have forced hunter-gatherer populations toward agriculture (10, 39). Regardless of the favored explanation, by showing that hunter-gatherer population density and the abundance and stability of their resources are strongly linked, our results highlight the difficulty in distinguishing between environmental and demographic causes of the subsistence change. Given hunter-gatherers’ well-documented ability to affect species communities (40), we can nevertheless assume that the higher the hunter-gatherer abundance that a given biome can potentially sustain, the stronger the impact of anthropogenic drivers on the biome has been even before the transition to agriculture.

In summary, our analyses revealed that NPP, biodiversity, and environmental pathogen stress interact to impose complex and varying limitations on hunter-gatherer population density in different parts of the world. These findings highlight the key role that environmental factors play in shaping the global population density patterns of preagricultural humans. Although cultural evolution has made humans enormously potent ecosystem engineers and has enabled humans to survive under a variety conditions (for example, by extracting energy in low-productivity environments and compiling knowledge on medical plants to fight against pathogens), it has not freed hunter-gatherer populations from biotic and abiotic forcing. We can, therefore, assume that spatial and temporal variability in this forcing has been important for our biological and cultural evolution as well.

Materials and Methods

Data.

All data manipulation and analyses were conducted using R statistical software (41). The R script and the datasets are freely available at the Zenodo repository (https://doi.org/10.5281/zenodo.1069786), which makes it possible to reproduce this study. SI Materials and Methods gives more details regarding datasets and methods.

We extracted the annual mean temperature and annual precipitation values from the WorldClim.org database (42), and based on these values, we calculated climatic NPP using the empirical Miami model (43) (Fig. S1A). The Miami model, although relatively simple, has been shown to give reasonable estimates of NPP in current climatic and vegetation conditions (44) (SI Materials and Methods). Instead of using satellite-based values, we decided to use climate-based NPP estimates to avoid the potential problems caused by recent human appropriations of NPP and to reach an estimate that reflects more natural levels of NPP.

We used combined mammal, bird, and vascular plant richness as a surrogate of overall biodiversity. Global mammal and bird richness data were extracted from the BiodiversityMapping.org database (45). Global native vascular plant richness data are model based, because native plant species’ ranges are insufficiently documented in most regions of the world (46, 47). We interpolated plant richness data to cover the same global grid with the same spatial resolution as the NPP and mammal and bird richness data using inverse distance weighting interpolation with the “gstat” package in R (48⇓–50) (SI Materials and Methods and Fig. S5 A and B). The biodiversity variable was calculated as an average of the scaled animal (mammal + bird) and plant richness (Fig. S1B).

The estimations of environmental pathogen stress are based on observations of the local pathogen conditions prevailing among 186 nonindustrial societies around the globe (51⇓–53). These societies form the Standard Cross-Cultural Sample (www.worldcultures.org/) (54). We used the pathogen prevalence index, which is an average of the z scores of the prevalence of 10 pathogens (malaria, dengue, filariae, typhus, trypanosomes, leishmaniasis, schistosomes, plague, leprosy, and spirochetes). We interpolated the pathogen data to cover the same global grid and the same spatial resolution as the NPP and biodiversity data using kriging interpolation and an exponential variogram model with the gstat package in R (48⇓–50) (SI Materials and Methods and Figs. S1C and S5 C and D). We assume that the global interpolation of the pathogen data reflects the relative pathogen stress levels in the local environments of the study populations.

The hunter-gatherer population density data were combined from two ethnographic datasets compiled by Binford (10) and Kelly (11), with the addition of Saami in northern Europe from another source (55). To reduce the uncertainty related to ethnographic population density estimates, if a population was present in both datasets, we averaged the population density values from both sources. When a hunter-gatherer population was present in only one of the datasets, we used the value given by that dataset (SI Materials and Methods and Fig. S6). For each population in the data, we extracted the values of the environmental explanatory variables from the global environmental layers of NPP, biodiversity, and environmental pathogen stress.

To analyze biomewise differences in population density, we used the Food and Agriculture Organization of the United Nations (FAO) Global Ecological Zones data on the distribution of 20 natural vegetation zones (28, 29) available at www.fao.org/geonetwork/srv/en/main.home?uuid=2fb209d0-fd34-4e5e-a3d8-a13c241eb61b.

Methods.

Visual inspections of the relationships between NPP and the loge population density of hunter-gatherers and between NPP and environmental pathogen stress indicated shifts in the effects of NPP. This impression was confirmed using breakpoint analysis with the segmented package in R, which allows for the fitting of piecewise linear regression models and provides estimates of breakpoints (56, 57). This package uses an iterative model-fitting algorithm. We did not provide any starting values for the breakpoint parameters, which means that the algorithm used the median of the predictor variable as a starting value for breakpoint estimation. The Davies test (58) included in the segmented package was performed to test whether the difference between the slopes before and after the estimated breakpoint was different from zero. We calculated the mean of the two NPP breakpoints (for population density and pathogen stress) and used this value to distinguish between low- and high-productivity environments.

Structural equation modeling is an approach that interprets information about the observed correlations between system components to evaluate potentially complex causal relationships (23). The approach begins with hypothesizing the underlying structure of causal pathways in the system (a conceptual model or metamodel). Then, this structure is translated into regression equations, and finally, these equations are evaluated against data to support or refute the hypothesized structure. Through this process, structural equation modeling provides an understanding of the nature and magnitude of direct and indirect effects in a system (23), which in our case, is composed of three environmental factors and human population density.

Here, we performed structural equation modeling as a piecewise estimation of local relationships, in which the parameters are estimated separately for each regression equation (59). The fit of the structure of the piecewise model is evaluated using a test of directed separation, with Fisher’s C as the test statistic (59, 60). This test indicates whether missing paths between variables in a system should be included in the model. We performed modeling using the R package piecewiseSEM, which allows for the automated fitting of lists of structured equations and goodness-of-fit tests for full and component models (59). For the sake of interpretability, we fit separate models for low- and high-productivity environments.

We began with the models (for low- and high-productivity environments) that omitted the direct effects of NPP on loge population density. For low- and high-productivity structural equation models, the goodness-of-fit test strongly indicated that the direct effect of NPP should be included. In such saturated models, however, the effect of pathogen stress on population density (low productivity) and the effect of biodiversity on both pathogen stress and population density (high productivity) were nonsignificant. After omitting these effects, the goodness-of-fit test indicated that the structural equation models fit the data well in both low- and high-productivity environments.

However, after adjusting the P values for spatial autocorrelation (61) (SI Materials and Methods), the effect of NPP on pathogen stress turned out to not be statistically significant in the low-productivity environments. Therefore, we omitted this effect from the final structural equation model, although Fisher’s C indicated a lower fit than for the model that included the effect of NPP on pathogen stress.

For global projections of hunter-gatherer population density, we used two regression methods: ordinary least squares linear regression model (LM) (Fig. S3 and Table S2) and support vector regression (SVR) (62⇓–64). Both of these models included interaction terms between biodiversity and NPP and between pathogen stress and NPP. LM and SVR produce highly similar predictions (Fig. S7). However, because fitting algorithms are different (62⇓–64), their predictions are not identical, and therefore, these techniques are complementary (Fig. S7). Neither of the techniques seem to extrapolate, and the predicted values are within the boundaries of the original hunter-gatherer data (Fig. S7). In the context of global prediction and extrapolation, we nevertheless feel that it safer to use a combination of two techniques instead of single technique.

We fit global models using NPP, biodiversity, and pathogen stress as predictors of loge population density. SVR was implemented with the e1071 package in R (65). The LM and SVR (with linear kernel) models were trained without splitting the ethnographic data into low- and high-productivity subsets. After the training, we used model algorithms to predict the loge population density using global environmental data. We calculated the mean of the predictions of the LM and SVR models and used this average as our global projection of hunter-gatherer population densities as shown in Fig. 3B. The performance of the averaged (LM and SVR) model was evaluated using h-block cross-validation (SI Materials and Methods) (66, 67); h-block cross-validated (h = 500 km) prediction root mean squared error is 1.26, and R2 is 0.46. These metrics indicate that the model is performing well with the unseen data.

The limiting effect of a predictor variable is calculated as the difference between the population density of the prediction, in which the given predictor is set to its optimum value, and the “true” prediction (see above), in which all of the predictor variables are allowed to vary. The larger the difference, the stronger the limiting effect of the predictor. The map of the limiting effects of NPP, biodiversity, and pathogen stress was constructed by determining the optimum values of each predictor in the averaged (LM and SVR) model. The optimum value of a predictor is the value where the loge population density reaches its highest value in the ethnographic data, as seen in the response shapes of the averaged (LM and SVR) model. We then created three global datasets, in which one predictor variable at a time was set to its optimum and the other two predictors were allowed to vary. Next, we produced three global projections of population density using these datasets. From these projections, we subtracted the true prediction. The resulting differences were scaled to vary between 0 and 255. The final map was produced using the plotRGB function in the “raster” package in R (68). As humans have had profound impact on biodiversity, the limiting effect of biodiversity may be artificially high in certain regions. This is because biodiversity in these regions is further away from its optimum value than it would be without the recent human influence.

For biomewise comparisons, we rasterized the biome dataset (see above) and aligned it with the global projected limiting effects and the population density data. This allowed us to extract information about the distributions of the limiting effects and the loge population density in different biomes. For the biomewise comparisons, the limiting effects were scaled to sum to unity.

Acknowledgments

We thank M. Fortelius, J. Kankaanpää, R. L. Kelly, T. Rankama, and three anonymous reviewers for their constructive critique. M.T. and J.T.E. acknowledge financial support from the Kone Foundation. J.T.E. acknowledges this as a contribution to the Integrative Climate Change Biology Programme under International Union of Biological Sciences.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: miikka.tallavaara{at}helsinki.fi.

  • Author contributions: M.T. and M.L. designed research; M.T. performed research; M.T. and M.L. contributed new reagents/analytic tools; M.T., J.T.E., and M.L. analyzed data; and M.T., J.T.E., and M.L. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission. R.L.K. is a guest editor invited by the Editorial Board.

  • Data deposition: The data and scripts reported in this paper have been deposited in the Zenodo repository, https://doi.org/10.5281/zenodo.1069786.

  • See Commentary on page 1137.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1715638115/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Kelly RL,
    2. Surovell TA,
    3. Shuman BN,
    4. Smith GM
    (2013) A continuous climatic impact on Holocene human population in the Rocky Mountains. Proc Natl Acad Sci USA 110:443–447.
    .
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Shennan S, et al.
    (2013) Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat Commun 4:2486.
    .
    OpenUrlPubMed
  3. ↵
    1. Tallavaara M,
    2. Luoto M,
    3. Korhonen N,
    4. Järvinen H,
    5. Seppä H
    (2015) Human population dynamics in Europe over the last glacial maximum. Proc Natl Acad Sci USA 112:8232–8237.
    .
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Goldberg A,
    2. Mychajliw AM,
    3. Hadly EA
    (2016) Post-invasion demography of prehistoric humans in South America. Nature 532:232–235.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Timmermann A,
    2. Friedrich T
    (2016) Late pleistocene climate drivers of early human migration. Nature 538:92–95.
    .
    OpenUrlPubMed
  6. ↵
    1. Birdsell JB
    (1953) Some environmental and cultural factors influencing the structuring of Australian aboriginal populations. Am Nat 87:171–207.
    .
    OpenUrlCrossRef
  7. ↵
    1. Baumhoff MA
    (1963) Ecological Determinants of Aboriginal California Populations (Univ of California Press, Berkeley, CA).
    .
  8. ↵
    1. Keeley LH
    (1988) Hunter-gatherer economic complexity and “population pressure”: A cross-cultural analysis. J Anthropol Archaeol 7:373–411.
    .
    OpenUrlCrossRef
  9. ↵
    1. Hamilton MJ,
    2. Milne BT,
    3. Walker RS,
    4. Brown JH
    (2007) Nonlinear scaling of space use in human hunter-gatherers. Proc Natl Acad Sci USA 104:4765–4769.
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Binford LR
    (2001) Constructing Frames of Reference: An Analytical Method for Archaeological Theory Building Using Ethnographic and Environmental Data Sets (Univ of California Press, Berkeley, CA).
    .
  11. ↵
    1. Kelly RL
    (2013) The Lifeways of Hunter-Gatherers: The Foraging Spectrum (Cambridge Univ Press, Cambridge, UK), 2nd Ed.
    .
  12. ↵
    1. Hamilton MJ,
    2. Burger O,
    3. Walker RS
    (2012) Human ecology. Metabolic Ecology, eds Sibly RM, Brown JH, Kodric-Brown A (Wiley, Chichester, UK), pp 248–257.
    .
  13. ↵
    1. Burger JR,
    2. Weinberger VP,
    3. Marquet PA
    (2017) Extra-metabolic energy use and the rise in human hyper-density. Sci Rep 7:43869.
    .
    OpenUrl
  14. ↵
    1. Loreau M
    (2000) Biodiversity and ecosystem functioning: Recent theoretical advances. Oikos 91:3–17.
    .
    OpenUrlCrossRef
  15. ↵
    1. Hooper DU, et al.
    (2005) Effects of biodiversity on ecosystem functioning: A consensus of current knowledge. Ecol Monogr 75:3–35.
    .
    OpenUrlCrossRef
  16. ↵
    1. Tilman D,
    2. Isbell F,
    3. Cowles JM
    (2014) Biodiversity and ecosystem functioning. Annu Rev Ecol Evol Syst 45:471–493.
    .
    OpenUrlCrossRef
  17. ↵
    1. de Mazancourt C, et al.
    (2013) Predicting ecosystem stability from community composition and biodiversity. Ecol Lett 16:617–625.
    .
    OpenUrl
  18. ↵
    1. Llamas B, et al.
    (2016) Ancient mitochondrial DNA provides high-resolution Time scale of the peopling of the Americas. Sci Adv 2:e1501385.
    .
    OpenUrlFREE Full Text
  19. ↵
    1. O’Fallon BD,
    2. Fehren-Schmitz L
    (2011) Native Americans experienced a strong population bottleneck coincident with European contact. Proc Natl Acad Sci USA 108:20444–20448.
    .
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Patterson JEH,
    2. Neuhaus P,
    3. Kutz SJ,
    4. Ruckstuhl KE
    (2013) Parasite removal improves reproductive success of female North American red squirrels (Tamiasciurus hudsonicus). PLoS One 8:e55779.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hochachka WM,
    2. Dhondt AA
    (2000) Density-dependent decline of host abundance resulting from a new infectious disease. Proc Natl Acad Sci USA 97:5303–5306.
    .
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Hudson PJ,
    2. Dobson AP,
    3. Newborn D
    (1998) Prevention of population cycles by parasite removal. Science 282:2256–2258.
    .
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Grace JB
    (2006) Structural Equation Modeling and Natural Systems (Cambridge Univ Press, Cambridge, UK).
    .
  24. ↵
    1. Gillman LN,
    2. Wright SD
    (2006) The influence of productivity on the species richness of plants: A critical assessment. Ecology 87:1234–1243.
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. Cusens J,
    2. Wright SD,
    3. McBride PD,
    4. Gillman LN
    (2012) What is the form of the productivity-animal-species-richness relationship? A critical review and meta-analysis. Ecology 93:2241–2252.
    .
    OpenUrlPubMed
  26. ↵
    1. Guernier V,
    2. Hochberg ME,
    3. Guégan J-F
    (2004) Ecology drives the worldwide distribution of human diseases. PLoS Biol 2:e141.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Dunn RR,
    2. Davies TJ,
    3. Harris NC,
    4. Gavin MC
    (2010) Global drivers of human pathogen richness and prevalence. Proc Biol Sci 277:2587–2595.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. FAO
    (2001) Global Ecological Zoning of the Global Forest Resources Assessment 2000. Final Report (Food and Agriculture Organization of the United Nations, Rome).
    .
  29. ↵
    1. FAO
    (2012) Global Ecological Zones for FAO Forest Reporting: 2010 Update (Food and Agriculture Organization of the United Nations, Rome).
    .
  30. ↵
    1. Currie DJ
    (1991) Energy and large-scale patterns of animal- and plant-species richness. Am Nat 137:27–49.
    .
    OpenUrlCrossRef
  31. ↵
    1. Badgley C,
    2. Fox DL
    (2000) Ecological biogeography of North American mammals: Species density and ecological structure in relation to environmental gradients. J Biogeogr 27:1437–1467.
    .
    OpenUrlCrossRef
  32. ↵
    1. Hortal J,
    2. Rodríguez J,
    3. Nieto-Díaz M,
    4. Lobo JM
    (2008) Regional and environmental effects on the species richness of mammal assemblages. J Biogeogr 35:1202–1214.
    .
    OpenUrlCrossRef
  33. ↵
    1. Schemske DW,
    2. Mittelbach GG,
    3. Cornell HV,
    4. Sobel JM,
    5. Roy K
    (2009) Is there a latitudinal gradient in the importance of biotic interactions? Annu Rev Ecol Evol Syst 40:245–269.
    .
    OpenUrlCrossRef
  34. ↵
    1. Froment A
    (2001) Evolutionary biology and health of hunter-gatherer populations. Hunter-Gatherers: An Interdisciplinary Perspective, Biological Society Symposium Series, eds Panter-Brick C, Layton RH, Rowley-Conwy P (Cambridge Univ Press, Cambridge, UK), pp 236–266.
    .
  35. ↵
    1. London D,
    2. Hruschka D
    (2014) Helminths and human ancestral immune ecology: What is the evidence for high helminth loads among foragers? Am J Hum Biol 26:124–129.
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Page AE, et al.
    (2016) Reproductive trade-offs in extant hunter-gatherers suggest adaptive mechanism for the neolithic expansion. Proc Natl Acad Sci USA 113:4694–4699.
    .
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Smith BD
    (2011) A cultural niche construction theory of initial domestication. Biol Theory 6:260–271.
    .
    OpenUrlCrossRef
  38. ↵
    1. Zeder MA
    (2012) The broad spectrum revolution at 40: Resource diversity, intensification, and an alternative to optimal foraging explanations. J Anthropol Archaeol 31:241–264.
    .
    OpenUrlCrossRef
  39. ↵
    1. Rosenberg M
    (1990) The mother of invention: Evolutionary theory, territoriality, and the origins of agriculture. Am Anthropol 92:399–415.
    .
    OpenUrlCrossRef
  40. ↵
    1. Boivin NL, et al.
    (2016) Ecological consequences of human niche construction: Examining long-term anthropogenic shaping of global species distributions. Proc Natl Acad Sci USA 113:6388–6396.
    .
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. R Core Team
    (2017) R: A Language and Environment for Statistical Computing. Version 3.3.3. Available at https://www.R-project.org/. Accessed March 6, 2017.
    .
  42. ↵
    1. Hijmans RJ,
    2. Cameron SE,
    3. Parra JL,
    4. Jones PG,
    5. Jarvis A
    (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25:1965–1978.
    .
    OpenUrlCrossRef
  43. ↵
    1. Lieth H
    (1973) Primary production: Terrestrial ecosystems. Hum Ecol 1:303–332.
    .
    OpenUrlCrossRef
  44. ↵
    1. Adams B,
    2. White A,
    3. Lenton TM
    (2004) An analysis of some diverse approaches to modelling terrestrial net primary productivity. Ecol Modell 177:353–391.
    .
    OpenUrl
  45. ↵
    1. Jenkins CN,
    2. Pimm SL,
    3. Joppa LN
    (2013) Global patterns of terrestrial vertebrate diversity and conservation. Proc Natl Acad Sci USA 110:E2602–E2610.
    .
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Kreft H,
    2. Jetz W
    (2007) Global patterns and determinants of vascular plant diversity. Proc Natl Acad Sci USA 104:5925–5930.
    .
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Ellis EC,
    2. Antill EC,
    3. Kreft H
    (2012) All is not loss: Plant biodiversity in the anthropocene. PLoS One 7:e30535.
    .
    OpenUrlCrossRefPubMed
  48. ↵
    1. Zimmerman D,
    2. Pavlik C,
    3. Ruggles A,
    4. Armstrong MP
    (1999) An experimental comparison of ordinary and universal kriging and inverse distance weighting. Math Geol 31:375–390.
    .
    OpenUrl
  49. ↵
    1. Pebesma EJ
    (2004) Multivariable geostatistics in S: The gstat package. Comput Geosci 30:683–691.
    .
    OpenUrl
  50. ↵
    1. Gräler B,
    2. Pebesma EJ,
    3. Heuvelink G
    (2016) Spatio-temporal interpolation using gstat. R J 8:204–218.
    .
    OpenUrl
  51. ↵
    1. Low BS
    (1994) Pathogen intensity cross-culturally. World Cult 8:24–34.
    .
    OpenUrl
  52. ↵
    1. Cashdan E,
    2. Steele M
    (2013) Pathogen prevalence, group bias, and collectivism in the standard cross-cultural sample. Hum Nat 24:59–75.
    .
    OpenUrlCrossRefPubMed
  53. ↵
    1. Cashdan E
    (2014) Biogeography of human infectious diseases: A global historical analysis. PLoS One 9:e106752.
    .
    OpenUrl
  54. ↵
    1. Murdock GP,
    2. White DR
    (1969) Standard cross-cultural sample. Ethnology 8:329–369.
    .
    OpenUrlCrossRef
  55. ↵
    1. Itkonen TI
    (1948) Suomen lappalaiset vuoteen 1945. Ensimmäinen osa. (WSOY, Porvoo, Finland).
    .
  56. ↵
    1. Muggeo VMR
    (2003) Estimating regression models with unknown break-points. Stat Med 22:3055–3071.
    .
    OpenUrlCrossRefPubMed
  57. ↵
    1. Muggeo VMR
    (2008) Segmented: An R package to fit regression models with broken-line relationships. R News 8:20–25.
    .
    OpenUrl
  58. ↵
    1. Davies RB
    (1987) Hypothesis testing when a nuisance parameter is present only under the alternatives. Biometrika 74:33–43.
    .
    OpenUrlCrossRef
  59. ↵
    1. Lefcheck JS
    (2016) piecewiseSEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol Evol 7:573–579.
    .
    OpenUrlCrossRef
  60. ↵
    1. Shipley B
    (2000) A new inferential test for path models based on directed acyclic graphs. Struct Equ Model Multidiscip J 7:206–218.
    .
    OpenUrl
  61. ↵
    1. Moran PAP
    (1950) Notes on continuous stochastic phenomena. Biometrika 37:17–23.
    .
    OpenUrlCrossRefPubMed
  62. ↵
    1. Drucker H,
    2. Burges CJC,
    3. Kaufman L,
    4. Smola AJ,
    5. Vapnik V
    (1997) Support vector regression machines. Advances in Neural Information Processing Systems 9, eds Mozer MC, Jordan MI, Petsche T (MIT Press, Cambridge, MA), pp 155–161.
    .
  63. ↵
    1. Vapnik VN
    (1998) Statistical Learning Theory (Wiley, New York).
    .
  64. ↵
    1. Smola AJ,
    2. Schölkopf B
    (2004) A tutorial on support vector regression. Stat Comput 14:199–222.
    .
    OpenUrlCrossRef
  65. ↵
    1. Meyer D,
    2. Dimitriadou E,
    3. Hornik K,
    4. Weingessel A,
    5. Leisch F
    (2015) e1071: Misc Functions of the Department of Statistics (e1071), TU Wien. R Package Version 1.6-7. Available at CRAN.R-project.org/package=e1071. Accessed January 1, 2017.
    .
  66. ↵
    1. Telford RJ,
    2. Birks HJB
    (2009) Evaluation of transfer functions in spatially structured environments. Quat Sci Rev 28:1309–1316.
    .
    OpenUrl
  67. ↵
    1. Salonen JS, et al.
    (2016) Calibrating aquatic microfossil proxies with regression-tree ensembles: Cross-validation with modern chironomid and diatom data. Holocene 26:1040–1048.
    .
    OpenUrlCrossRef
  68. ↵
    1. Hijmans RJ
    (2016) raster: Geographic Data Analysis and Modeling. R Package Version 2.5-8. Available at https://CRAN.R-project.org/package=raster. Accessed January 1, 2017.
    .
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Drivers of the hunter-gatherer population density
Miikka Tallavaara, Jussi T. Eronen, Miska Luoto
Proceedings of the National Academy of Sciences Feb 2018, 115 (6) 1232-1237; DOI: 10.1073/pnas.1715638115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

You May Also be Interested in

A study uncovers genetic evidence of early colonization of Jamaica by primates and suggests that adaptation can shape species morphology in novel environments, even when species come from morphologically conservative families. Image courtesy of Lorraine Meeker (American Museum of Natural History, New York).
Primate colonization of the Caribbean
A study uncovers genetic evidence of early colonization of Jamaica by primates and suggests that adaptation can shape species morphology in novel environments, even when species come from morphologically conservative families.
Image courtesy of Lorraine Meeker (American Museum of Natural History, New York).
Researchers report that some moth species use sound absorbing-scales to avoid detection by bats and suggest that the findings could aid the design of biologically inspired sound absorbers for an array of applications, including noise mitigation and building acoustics. Image courtesy of Thomas R. Neil.
How moth scales help avoid detection by bats
Researchers report that some moth species use sound absorbing-scales to avoid detection by bats and suggest that the findings could aid the design of biologically inspired sound absorbers for an array of applications, including noise mitigation and building acoustics.
Image courtesy of Thomas R. Neil.
Proceedings of the National Academy of Sciences: 115 (50)
Current Issue

Submit

Jump to section