Introduction

The asymmetrical geographic distribution of life on Earth is often described by the latitudinal gradient of biodiversity1,2,3,4,5,6,7,8,9,10,11. Across many taxa, including mammals, birds, reptiles, amphibians12,13, plants14,15,16, marine bivalves17, ants18, liverworts19, and diseases20, species richness is highest near the equator and decays exponentially toward the poles. Moreover, these biodiversity gradients extend deep into the geological past8,9,21,22. It is now well-established that the biogeographic distribution of human cultural, linguistic, and economic diversity shows similar gradients1,3,4,23,24. While there is no unified mechanistic theory of biogeography to explain the origins of these biodiversity gradients, proposed hypotheses tend to focus on ecological (i.e., niche diversity), energetic (i.e., ecosystem kinetics), evolutionary (i.e., phylogenetic history), and/or historical mechanisms (i.e., dispersal)12,14,25,26,27,28,29. For cultural diversity similar mechanisms have been proposed, including latitudinal gradients in the length of the growing season23,24, environmental gradients1,4,30, and the timing of initial human colonization of particular continents2.

While latitudinal gradients may be empirically pervasive, latitude itself cannot be an explanatory variable as it is simply a geographic binning procedure that captures the steady loss biodiversity with increasing distance from the equator. Latitude indexes large-scale variation of Earth systems at many scales including the intensity of solar radiation, the length of the growing season, the distribution of precipitation, and the statistics of environmental temperature, all of which interact to constrain the distribution of environmental productivity around the planet31. Exactly how the productivity of local environments influences biodiversity has been central to ecological and evolutionary theory since Darwin and Wallace32.

Few studies have proposed mechanistic theory to explain why the global distributions of biodiversity and cultural diversity are so similar3,4,6. Ultimately, diversity must be correlated with the flux of energy through ecosystems as all organisms compete for environmental production to fuel metabolic processes of growth, maintenance and reproduction. Because environmental and climatic variation constrains the metabolic fluxes of free energy available to support biomass ecosystems differ in their capacity to support life26,31 and so the biogeographic distribution of biodiversity necessarily covaries with energy availability. The biogeographic distribution of human biomass is similarly constrained by geographic variation in energy availability despite technological innovations that enhance the ability of humans to compete for resources by reshaping their environments and redistributing materials, energy and information at multiple scales33,34,35,36,37.

To explain the ecological link between environmental production and species diversity the “more-individuals-hypothesis” proposed that increasingly productive environments support proportionally more individuals that belong to more species resulting in higher local diversity25,38,39. A similar hypothesis could be proposed for cultural diversity as human populations tend to be denser in more productive environments, which may lead to increased diversity. It is increasingly apparent, however, that species diversity is not simply a linear function of environmental productivity as the link between diversity and environmental productivity is more complex26. For example, while more productive environments house more species, tropical ecosystems house more diversity at all scales than would be expected by environmental production alone15,26. Similarly, studies of human cultures show that while diversity tends to be higher in more productive environments1,3,4,6,23,24,40 there are many other environmental, climatic, geological, and cultural processes that impact diversity at various scales40,41,42,43,44,45. For example, while global variation in hunter-gatherer population density, space use, and mobility is well-predicted by environmental productivity44,46,47, the spatial distribution of ethnolinguistic groups in general across the socioeconomic spectrum is better predicted by their level of sociopolitical complexity45,48.

In this paper, we use these studies as motivation to examine the macroecological relationship between mammal species diversity and human cultural diversity through the lens of the metabolic theory of ecology. Our goal in this paper is to explain why human cultural and mammalian diversity are so similarly distributed at the global scale and why both forms of diversity increase faster than would be predicted from environmental production alone. First, we derive predictive theory from the metabolic theory of ecology, which describes explicitly how the temperature-dependence of metabolic processes scales up to drive the temperature-dependence of species and cultural diversity. We then use this theory to develop a set of hypotheses that we test with global databases.

Theory development

Individual metabolism, temperature, and body size

We begin by considering the well-established link between individual and ecosystem level metabolism49. As first observed by Kleiber50, the basal metabolic rate of an organism, B, can be expressed in terms of body mass, M, as

$$ B = b_{0} M^{\alpha } $$
(1)

where b0 is a taxon-specific constant, independent of mass and temperature, and \(\alpha \approx 3/4\) ( see51,52,53,54 for further discussions of the empirical form of Kleiber’s Law). The metabolic theory of ecology derives theory that explains how the scaling of mass and metabolism arises from the fractal-like, space-filling nature of optimized internal networks that distribute the energy and resources that sustain organisms, such as the vascular system51,55,56,57,58,59. The scale parameter, b0, is derived from the underlying biochemical kinetics of metabolism and, as such, depends on the temperature, T, at which an organism operates57. This is given by the exponential Arrhenius-Boltzmann factor, \({\text{exp}}\left( { - E_{B} /kT} \right)\), where EB (~ 0.65 eV) is the average activation energy of the biochemical reactions contributing to metabolic processes of respiration, k is Boltzmann’s constant (8.62 × 10–5 eV K−1), and T is the absolute temperature at which the organism operates (°K)60. For endotherms, T is the internal body temperature, and for ectotherms, T is the ambient environmental temperature. Incorporating this temperature dependence into Eq. (1) yields the expression

$$ B = b_{0} e^{{ - E_{B} /kT}} M^{\alpha } . $$
(2)

which is the central equation of metabolic theory of ecology57,60,61.

Scaling up from body size to ecosystem metabolism

As ecosystems are composed of individual organisms, it follows that the respiration of an entire ecosystem is the sum over all individuals, the majority of which are microbes and plants49. Temperature is one of the primary drivers of energy flux in ecosystems62 the largest of which is gross primary production, GPP, the total mass of carbon entering an ecosystem through the conversion of solar radiation into biomass by photosynthesis over a given period of time63,64. At steady state, about half of GPP is respired by plants to support production, maintenance, and ion uptake65. The remainder is used by plants to fuel growth, and this is terrestrial net primary production (NPP), the net carbon gain in an ecosystem over a given period of time measured in units of g C m−2 yr−1. Following mass-energy equivalence, NPP is the total amount of ecosystem energy available for work (i.e., Gibbs free energy). It then follows that the flux of free energy in ecosystems is captured by the Arrhenius-Boltzmann temperature-dependence \(NPP \propto e^{{ - E_{NPP} /kT}}\)49,66. If \(E_{NPP}\) is assumed to be the activation energy of photosynthesis, then holding all other ecological and evolutionary processes constant \(E_{NPP} \approx 0.3 \) eV. However, all is not constant as the kinetics of NPP responds to other rate limiting constraints such as stoichiometry, climatic constraints, and water-availability26,31,62,67. Here we measure the temperature-dependence of NPP from data: introducing annual precipitation, R (mm yr−1), as a rate-limiting constraint on the kinetics of NPP, we have

$$ NPP = c_{1} R^{{\beta_{NPP} }} e^{{ - E_{NPP} /kT}} , $$
(3)

where \(\beta_{NPP}\) is an exponent capturing the effect of precipitation on NPP and c1 is a normalization constant independent of precipitation or temperature. We fit a multiple regression model to the logarithmic-transform of Eq. (3) and present the results in Table 1. Results show global scale variation in NPP is well-predicted by temperature and precipitation (\(R^{2} = 92\%\)), where the kinetics have a slope \(E_{NPP} = 0.52\) eV and precipitation is a sublinear constraint on NPP, \({ }\beta_{NPP} = 0.72\). Figure 1 illustrates these relationships in a 3D surface plot of NPP as a function of temperature and precipitation.

Table 1 Regression coefficients for the multiple regression model of NPP as a function of inverse temperature 1/kT and precipitation R.
Figure 1
figure 1

A 3D surface plot of net primary production (z) as a function of inverse temperature (x) and precipitation (y) in our data. The surface is a LOESS fit.

The kinetics of ecosystem complexity and biodiversity

Although energy availability within ecosystems is ultimately determined by temperature-dependent fluxes of carbon, increasingly productive environments are also more structurally complex. That is to say, a polar desert is not simply an energy-poor tropical rainforest, and a tropical rainforest is not simply an energy rich grassland. Increasingly productive environments not only have greater energy flux per unit area but also house more species26. Evolutionary explanations suggest speciation and extinction dynamics vary consistently across environments affecting the standing stock of diversity, and data suggest the tropics are more taxonomically and genetically diverse as a result68,69. Niche-based explanations suggest that in productive environments resources are apportioned in increasingly different ways, thus increasing the number of species that may co-exist in an ecosystem, and phylogenetically-related species have a tendency to inhabit similar niches in similar environments70,71,72. Ecological speciation explanations suggest increasing interactions between species in more productive environments promote enhanced diversity through increased competition and constraints on dispersal. The metabolic theory of ecology suggests a more general kinetic explanation, not necessarily exclusive of the preceding theories; because the pace of life is driven by metabolic processes, all the ecological and evolutionary dynamics that constitute the structure of ecosystems are governed by biological kinetics where biotic interactions of all kinds occur faster when the environment is warmer26,57. Diversity increases both as a function of more energy and faster biotic interactions in more productive environments.

Because human societies are also necessarily embedded within ecosystems, the abundance, density, and the intraspecific diversity of the human species are similarly subject to ecological kinetics. Groups of human compete with each other for finite sources of space and energy, as well as with other species, and cultural diversity is the result of intraspecific group competition and other evolutionary processes playing out over multiple generations73,74,75. In this paper we use ethnolinguistic diversity as our measure of cultural diversity. The term ethnolinguistic refers to a spatially and culturally discrete population of language speakers, which may or may not share a common language with other groups around the planet, such as English speakers in the British Isles, Australia, or the Caribbean, or Diné speakers in Alaska or the American Southwest. In both of these cases, each language is composed of multiple ethnolinguistic populations.

We model the kinetics of diversity by considering that the number of species increases with environmental productivity, and as the number of species increases there are more interactions between those species. We write the kinetics of diversity as a function of both the kinetics of energy flux and all other biotic interactions:

$$ S = c_{2} e^{{ - E_{S} /kT}} \propto e^{{ - \left( {E_{NPP} + E_{I} } \right)/kT}} , $$
(4)

where \(E_{S}\) is the overall kinetics of diversity and \(E_{I}\) is the kinetics of ecological interactions. Here, the kinetics of diversity, \(E_{S}\), will be faster than the kinetics of ecosystem energy flux, \(E_{NPP}\), when the kinetics of biological interactions \(E_{I} > 0\). From Eqs. (3) and (4) we can then express species diversity, S, as a function of NPP yielding

$$ S = c_{2} NPP^{{\beta_{S} }} $$
(5)

where \(\beta_{S} = E_{S} /E_{NPP} = \left( {E_{NPP} + E_{I} } \right)/E_{NPP} = 1 + \delta_{S}\), and so \(\delta_{S} = E_{I} /E_{NPP}\). Therefore, diversity will increase superlinearly (\(\beta_{S} > 1\)) with environmental productivity whenever the kinetics of diversity are faster than the kinetics of environmental productivity (\(E_{S} > E_{NPP}\)). If there are no interaction kinetics (\(E_{I} = 0\)), diversity will increase linearly with environmental productivity, and if interaction kinetics are negative (\(E_{I} < 0\))—from competitive exclusion, for example—diversity will increase sublinearly.

Four hypotheses

To establish whether this simple model of ecological kinetics captures variation in mammal and cultural diversity we use data on the global distribution of mammal species and ethnolinguistic populations to test four hypotheses using predictions derived from this model:

H1:

Diversity is a function of kinetics. We first test the hypothesis that diversity increases exponentially with temperature, as predicted by ecological theory at a rate consistent with the kinetics of environmental metabolism.

H2:

Diversity kinetics is faster than productivity kinetics. The kinetics of diversity, Es, should be faster than the kinetics of environmental productivity, \(E_{NPP}\), as species interaction rates, \(E_{I}\), are predicted to be both non-zero and increase with environmental productivity, \( NPP\).

H3:

Diversity is more than environmental productivity. Diversity, S, should increase with environmental productivity, P, at a superlinear rate given by \(\beta_{S} = E_{S} /E_{NPP} > 1\).

H4:

Diversity kinetics are faster in more productive environments. Holding NPP constant, diversity should increase with temperature, i.e., for a given range of NPP relatively warmer environments will be more diverse, and the slope of this response should increase in more productive environments.

Results

Hypotheses 1 and 2: Latitude and exponential temperature gradients

Figure 2 shows that the spatial distribution of mammal species (red points) and languages (blue points) across the surface of the planet are similar. In Fig. 2A, each data point is the centroid of a geographic distribution. Inset along the Y-axis of Fig. 2A are distributions of the relative frequency of ethnolinguistic groups and the number of mammal species by latitude. Figure 2B shows that the average size of mammal species ranges is significantly greater than language ranges, but Fig. 2C show that geographic ranges do not vary systematically with absolute latitude. As such, the centroids of ranges shown in Fig. 2A accurately capture the biogeographic distribution of both forms of diversity around the planet.

Figure 2
figure 2

The biogeographic distribution of mammal species and human cultural diversity. (A) A map of the global distribution of human ethnolinguistic groups (blue) and mammal species (red). Each data point is the centroid of the geographic range of a species or a language (see “Methods” for further discussion). The underlying map is composed of 56,598 cells (0.5*0.5-degree latitude grids) comprising the global sampling procedure. Inset along the y-axis are distributions of relative diversity controlling for available landmass per 1-degree latitudinal bin (i.e., centroid density). Both forms of diversity are highest near the equator and decay toward the poles. The Figure was generated in Sigmaplot v. 13 (www.systatsoftware.com). (B) Relative frequency distributions of geographic ranges for languages and mammals, where the solid vertical lines are the means. The geographic ranges of mammals are on average about 250-times larger than those of languages. (C) A plot of geographic range of as a function of absolute latitude for both forms of diversity shows no systematic relationship indicating that the centroids of geographic ranges capture the biogeographic distribution. However, mammal geographic ranges are 250-times larger than ethnolinguistic ranges, as shown in (B).

Figure 3A shows that after controlling for available landmass cultural diversity and mammal diversity exhibit similar responses to absolute latitude. Here we plot relative diversity, \(S^{\prime} = S_{i} /S_{{{\text{max}}}}\), for both mammal and cultural diversity and fit models of the form \({\text{ln}}S^{\prime} = {\text{ln}}S^{\prime}_{0} - \alpha L\), where \(\alpha\) is the slope of the response and L is latitude. Results show cultural diversity decreases with absolute latitude at a rate \(\alpha_{E} = 0.05 \pm 0.004\) (95% C.I.) (OLS regression: d.f. = 73, r2 = 83%, p < 0.0001) and mammalian diversity as \(\alpha_{M} = 0.05 \pm 0.01\) (OLS regression: d.f. = 68, r2 = 76%, p < 0.0001); both forms of diversity decrease at a rate of ~ 5% per degree latitude, and the correlation coefficient between the two data sets is 0.84 (Fig. 3A inset).

Figure 3
figure 3

Latitudinal gradients (A) and temperature gradients (B) in relative ethnolinguistic diversity (blue points), and relative mammal species diversity (red points), corrected by available landmass (i.e., the area per 0.5*0.5 sampling bin). In both cases both forms of diversity decrease exponentially with latitude and temperature at statistically indistinguishable rates. Inset figures are the Pearson’s r correlations between the two forms of diversity. Solid colored lines are OLS regression fits.

In support of hypothesis 1, Fig. 3B shows the exponential relationship between relative diversity and environmental temperature (in °C (mean annual temperature) along the upper x-axis and inverse temperature 1/kT along the lower x-axis) for both forms of diversity. We fit regression models of the form \({\text{ln}}S^{\prime} = {\text{ln}}S^{\prime}_{0} + E/kT\) is fit to both data sets. Cultural diversity decreases with temperature at a rate \(E_{L} = - 0.70\) eV (0.65–0.76 95% C.I.) (OLS regression: d.f. = 65, r2 = 91%, p < 0.0001) and mammal species diversity decreases with temperature at a rate \(E_{M} = - 0.71 \) eV (0.64–0.78) (OLS regression: d.f. = 60, r2 = 87%, p < 0.0001: for both data sets, statistical fits excluded the first three leftmost data points as they are outliers and have undue influence on the fit of the slope). The mammal and cultural slopes are not significantly different from each other and the two forms of diversity are highly correlated in space; \(r = 0.92\) (Fig. 3B inset). As such, the latitudinal gradients of diversity in Fig. 3A are proxies of the temperature-dependent diversity shown in Fig. 3B. Importantly, the slopes of both the mammal and cultural fits are significantly steeper than the slope of NPP shown in Fig. 1 (EP = 0.52 eV), therefore providing support for hypothesis 2: The kinetics of diversity are faster than the kinetics of productivity.

Hypothesis 3: Superlinear scaling of diversity and net primary production

Figure 4 shows that both forms of diversity scale superlinearly with net primary production, thus providing support for hypothesis 3. Cultural diversity increases superlinearly with NPP at a rate \(\beta_{L} = 1.40\) (1.25–1.55) (OLS regression: d.f. = 40, r2 = 91%, p < 0.0001) and mammalian diversity increases as \(\beta_{M} = 1.23\) (1.10–1.35) (OLS regression: d.f. = 36, r2 = 91%, p < 0.0001) (fits exclude the rightmost data point as it is an outlier and has a large influence of the fit of the slope). Both slopes are significantly greater than 1 and the data are highly correlated; \(r = 0.97\) (inset figure). Moreover, the confidence intervals around the slopes both encompass their theoretically-predicted values, thus providing support for hypothesis 4. For mammals, the predicted slope based on the ratio of the kinetic exponents is \(\lambda_{M} = E_{M} /E_{P} = 0.70/0.52 = 1.35\), which falls within the observed confidence interval 1.10–1.35. For cultural groups the predicted slope is \(\lambda_{L} = E_{L} /E_{P} = 0.71/0.52 = 1.37\), which falls within the observed confidence interval 1.25–1.55.

Figure 4
figure 4

Ethnolinguistic and mammalian diversity per thousand square kms by net primary production. Inset figure is the Pearson’s r correlation of ethnolingistic and mammalian diversity. Solid colored lines are OLS regression fits; blue = ethnolinguistic and red = mammal.

Hypothesis 4: Diversity kinetics and net primary production

Figure 5A and B show that both forms of diversity increase exponentially with temperature at faster rates in more productive environments. Here the diversity data are aggregated into bins of width 1lnNPP from the 0.5*0.5-degree latitude sampling cells, and linear functions are fit to the data within each bin. Details of the regression models and the statistical results are given in the Supplementary Information. These figures demonstrate that for both cultural groups and mammal species diversity begets diversity as the interaction rates in more productive environments increase with increasing diversity, thus providing support for hypothesis 4.

Figure 5
figure 5

Ethnolinguistic (A) and mammal species (B) diversity by inverse temperature in bins of 1lnNPP. Colored lines are OLS regression fits. The inset figures plot the slopes of each bin against the lnNPP bin showing the general trend of faster kinetics with higher NPP. The data show the temperature dependent kinetics of diversity tend to be faster in more productive environments.

Discussion

Our results show that diversity begets diversity for both mammal species and cultural diversity because the kinetics of diversity are faster than the kinetics of environmental productivity resulting in superlinearity between diversity and environmental productivity. Increasingly productive environments are also more interactive meaning that not only is more energy available to support individuals, species, and communities, but there are increasingly different ways in which that energy is apportioned in the environment and thus competed over. More productive environments are also more complex, allowing more ways for diversity to be housed per unit area. To test this, we derived a simple model of diversity kinetics from well-established ecological theory from which we proposed four hypotheses whose predictions could be tested with data. Our results provide support for the four hypotheses posed in this paper, suggesting both mammal and ethnolinguistic diversity are positively density-dependent.

Our first hypothesis proposed that diversity is a function of kinetics. For both mammal species and cultural populations diversity increases with temperature exponentially at rates consistent with ecological kinetics (Fig. 3B)62. As such, a linear increase in temperature leads to a multiplicative increase in diversity and the rate of this increase is well-predicted by ecological theory indicating that the rate of increase in diversity is consistent with the overall kinetics of ecosystem metabolism. This means that for human cultures richer interactive environments are increasingly fragmentary and thus more culturally diverse and so increasing interaction results in more competition, not cooperation. Moreover, the kinetics of cultural diversity are mathematically and statistically similar to the kinetics of mammalian diversity, a process which happens across a much broader range of taxonomies, body sizes, life histories, and ecologies. Perhaps increasingly complex environments provide an increasingly diverse range of niches for humans to exploit economically by adopting a range of different lifestyles and technological innovations, such as foraging, farming, or horticulture. If so, the results shown here suggest that ecological kinetics constrain the availability of ecological niches for humans and mammals in similar ways.

Our second hypothesis predicted that kinetics of diversity is faster than the kinetics of environmental productivity. This is because the kinetics of diversity is predicted to be a combination of the kinetics of productivity and interaction. Results show that this prediction is supported by data where the kinetics of both mammal and cultural diversity (\(E_{M} = 0.70\) and \(E_{L} = 0.71\), respectively) are both significantly faster than the kinetics of NPP, \({ }E_{NPP} = 0.52\). Therefore, cooler less productive environments lose ethnolinguistic groups and mammalian species at a faster rate than would be predicted by the decrease in productivity alone. Thus, the change in diversity across environments is not only a function of energy availability, but niche complexity.

Our third hypothesis predicted the scaling of diversity and productivity to be superlinear. Data show that for both mammal species and human cultures diversity scales superlinearly with environmental productivity (for mammals \(\beta_{M} = 1.23\) and for human cultures \(\beta_{L} = 1.40\)), and the scaling exponent is consistent with the predicted ratio of kinetic terms (Fig. 3B). This superlinearity describes increasing returns to scale in diversity with environmental productivity; for mammalian species a doubling of environmental productivity leads to ~ 23% more diversity than would be predicted by environmental energy availability alone, and in human cultures this increase is ~ 40%. The exact mechanisms that drive these increases are the subject of ongoing work.

Our fourth hypothesis predicted that as environmental productivity increases the temperature-dependence of the diversity response should also increase (Fig. 4). This result may be the strongest demonstration of the kinetic hypothesis as it shows holding environmental productivity constant, warmer environments are not only more diverse than less productive environments, but the rate of increase in diversity is faster in more productive environments. Richer environments house more species that interact with each other at increasingly faster rates. Therefore, mammalian and cultural diversity are both consistent with Red Queen dynamics driven by the temperature-dependence of interaction rates76.

This biogeographic convergence of biological and cultural diversity is striking given the differences in the evolutionary mechanisms and time scales of biological and cultural evolution77,78. Perhaps most obviously, by definition, mammal species diversity is interspecific whereas cultural diversity is intraspecific. Mammalian speciation occurs on time scales of millions of years while human languages speciate orders of magnitude faster at time scales of hundreds to thousands of years. Moreover, the pathways and currencies of biological and cultural evolution differ widely79,80,81, and so an a priori explanation of why these different evolutionary mechanisms should converge on similar biogeographic distributions is not obvious from the dynamics of the evolutionary mechanisms themselves.

However, it is interesting to note that while the centroid density of diversity is similar (Fig. 2A), the geographic ranges of mammal species and ethnolinguistic populations are quite different. Figure 2B and C show that on average mammalian geographic ranges are nearly 250-times larger than ethnolinguistic ranges. This suggests that while the centroids of their distributions converge, there are different dispersal mechanisms that impact the magnitude of their ranges. The larger geographic ranges of mammals could reflect the greater opportunities for species to co-occur in space across a broader range of body sizes and niches than is possible for humans. It may also be the case that faster rates of diversification in languages compared to species results in more restricted geographic ranges. As a result, the spatial turnover of mammalian and cultural diversity will differ at different spatial scales (i.e., differences in alpha, beta, and gamma-diversity82).

While biological and cultural evolutionary processes might be similar in principle35,73,83,84 and even share similar mathematical representations, the detailed mechanisms of the two evolutionary processes are entirely different. Biological evolution is a statistical consequence of the differential transmission of genetic and epigenetic information over generations in response to environmental selective constraints that impact the ability of individuals to compete for finite resources and produce viable offspring, and where speciation is the macroevolutionary consequence over large time scales85,86. Cultural evolution is a statistical consequence of the differential transmission of social information at many scales of social organization via multiple transmission pathways in both space and time that similarly impacts the ability of individuals to compete for finite resources and produce viable offspring73,74,84. Languages, cultures, and economies turnover at much faster time scales the biological species, from the near instantaneous to many thousands of years77,87. One axiomatic implication of the results shown here is that the gradients of diversity that emerge in response to ecological kinetics must be independent of specific evolutionary mechanism (i.e., genes vs. culture), otherwise they would be very differently distributed in space. It may be the case that large-scale geological, climatic, and environmental constraints impact speciation and extinction processes in both systems in similar ways across gradients of environmental productivity, despite the different evolutionary inheritance mechanisms and time scales involved. As such, the standing stock of diversity at any one location on the planet may be quite different between mammal species and cultures, but the gradients are similar.

Whatever the specific similarities or differences between biological and cultural evolution the general mechanism of convergence proposed here is kinetic. In general, more productive environments are warmer than less productive environments, and with increasing temperature not only is more energy available but there are increasing rates of interaction, including competition, specialization, and mutualisms, all of which combine to increase diversity above levels predicted by environmental productivity alone. Diversity increases superlinearly with environmental productivity because more productive environments are also more complex.

Methods

Temperature and environmental productivity data

We used the WorldClim database for mean annual temperature in °C and mean annual precipitation in mm yr−1 at 10 arc min resolution (version 1.4; https://www.worldclim.org/bioclim). WorldClim provides interpolated global climate data for stations with 10 to 30 years of data88. We used estimates of annual net primary productivity (g C m−2 yr−1) from an average of 17 global models at a spatial resolution of 0.5 degrees89. All data were first projected with the WGS 1984 Cylindrical Equal Area projection and then processed across our language polygons using the Zonal Statistics tool in ArcGIS. The Global Self-consistent, Hierarchical, High-resolution Shoreline dataset was used as the basemap90.

Mammal and ethnolinguistic data

Ethnolinguistic data were downloaded from the Ethnologue, and centroids (in decimal degrees of latitude and longitude) were recorded for each individual language polygon (N = 7627). Terrestrial mammal data were downloaded from PanTHERIA and centroids (in decimal degrees of latitude and longitude) were recorded for each species (N = 5089). Average temperature, annual precipitation, and net primary productivity were the extracted for each ethnolinguistic population and mammal species. Data layers were compiled in ArcGIS and then analyzed statistically in R91. We compiled a GIS of data layers for each of the variables used in the paper.

In this paper we choose to analyze the centroids of geographic ranges rather than the entire geographic range for several reasons. First, not all the geographic ranges for all languages and mammal species are available. Second, it is not clear whether the ranges that are available are directly comparable. This is because the size of geographic ranges can be estimated many ways, and the distribution of species or languages within the geographic ranges can be complex. Population densities are not constant across ranges and the entire geographic range within any available boundary will not be sampled evenly. Therefore, having confidence in the complete geographic distribution of species or language speakers within their ranges (and the lacunae within them) requires additional levels of sampling not currently available. As such, while imperfect, we chose the centroid of the species and language geographic distributions as these metrics are more robust to measurement error while capturing the approximate central tendency of spatial diversity.

Sampling procedure and statistics

Measuring diversity requires a sampling strategy. We use a statistical sampling technique designed to isolate the independent variable and reduce spatial autocorrelation. First, we divided the surface of the planet into a grid of 0.5 degrees of latitude by 0.5 degrees of longitude, extracted the non-oceanic cells (for a total of n = 56,598 cells) and calculated the available landmass per cell. For each cell we then counted the number of cultural populations and mammal species centroids that occur in that cell, and then divided that occurrence data by the available landmass. We also calculated the average temperature and net primary production for each cell. To sample global diversity, we then created bins of 0.1 lnNPP and 0.1 1/kT, and summed the total occurrence of ethnolinguistic and mammalian species densities, thus providing measures of the diversity per bin of independent variable controlling for the available landmass within that bin. This sampling approach has several important statistical advantages. First, we aggregate the data into bins by small increments of the independent variable, and then take statistics over the bin. This technique has the advantage of using the degrees of freedom available from the data set to estimate statistical moments within bins, and therefore isolates the effect of the independent variable on the dependent variable. Statistical error is thus minimized when estimating scaling exponents. Second, the binning procedure minimizes the effects of spatial autocorrelation as data points are aggregated from non-local populations collected from around the planet and averaged. Thus, data used in the scaling analyses are aggregates and are statistically independent. Third, the binning method minimizes the potential weighting biases of heteroscedasticity and unequal sampling across the data set.

Aside from the binning procedure, we do not attempt control for phylogenetic history in either of the datasets in this paper. This is because there is no phylogenetic supertree of human languages and so it is not statistically possible to control for the evolutionary structure of all languages in a single global database. Moreover, within the field of historical linguistics the concept of a supertree is highly controversial given that the speed at which languages evolve limit the ability to detect common ancestries beyond a few thousand years. Individual phylogenetic trees for most major language families of the world are now available, but not all, and some major geographic coverage is lacking, most notably for the Americas and Native American languages as a whole. While there is a phylogenetic supertree available for mammals we do not attempt to conduct phylogenetic comparisons on one database because this is not possible for the other. Therefore, comparing the two would be highly problematic and inconsistent. However, in the text we note that our binning procedure reshuffles data points distributed in space and therefore partially controls for the spatial autocorrelation inherent in evolutionary relationships, though we make no claim this is complete, nor equivalent to formal phylogenetic controls.

The data used in this paper are available in the Supplementary Material.