Diversity begets diversity in mammal species and human cultures

Across the planet the biogeographic distribution of human cultural diversity tends to correlate positively with biodiversity. In this paper we focus on the biogeographic distribution of mammal species and human cultural diversity. We show that not only are these forms of diversity similarly distributed in space, but they both scale superlinearly with environmental production. We develop theory that explains that as environmental productivity increases the ecological kinetics of diversity increases faster than expected because more complex environments are also more interactive. Using biogeographic databases of the global distributions of mammal species and human cultures we test a series of hypotheses derived from this theory and find support for each. For both mammals and cultures, we show that (1) both forms of diversity increase exponentially with ecological kinetics; (2) the kinetics of diversity is faster than the kinetics of productivity; (3) diversity scales superlinearly with environmental productivity; and (4) the kinetics of diversity is faster in increasingly productive environments. This biogeographic convergence is particularly striking because while the dynamics of biological and cultural evolution may be similar in principle the underlying mechanisms and time scales are very different. However, a common currency underlying all forms of diversity is ecological kinetics; the temperature-dependent fluxes of energy and biotic interactions that sustain all forms of life at all levels of organization. Diversity begets diversity in mammal species and human cultures because ecological kinetics drives superlinear scaling with environmental productivity.

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 planet 31 . Exactly how the productivity of local environments influences biodiversity has been central to ecological and evolutionary theory since Darwin and Wallace 32 .
Few studies have proposed mechanistic theory to explain why the global distributions of biodiversity and cultural diversity are so similar 3,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 life 26,31 and so the biogeographic distribution of biodiversity necessarily covaries with energy availability. The biogeographic distribution of human Scientific Reports | (2020) 10:19654 | https://doi.org/10.1038/s41598-020-76658-2 www.nature.com/scientificreports/ 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 scales [33][34][35][36][37] .
To explain the ecological link between environmental production and species diversity the "more-individualshypothesis" proposed that increasingly productive environments support proportionally more individuals that belong to more species resulting in higher local diversity 25,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 complex 26 . For example, while more productive environments house more species, tropical ecosystems house more diversity at all scales than would be expected by environmental production alone 15,26 . Similarly, studies of human cultures show that while diversity tends to be higher in more productive environments 1,3,4,6,23,24,40 there are many other environmental, climatic, geological, and cultural processes that impact diversity at various scales [40][41][42][43][44][45] . For example, while global variation in hunter-gatherer population density, space use, and mobility is well-predicted by environmental productivity 44,46,47 , the spatial distribution of ethnolinguistic groups in general across the socioeconomic spectrum is better predicted by their level of sociopolitical complexity 45,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 metabolism 49 . As first observed by Kleiber 50 , the basal metabolic rate of an organism, B, can be expressed in terms of body mass, M, as where b 0 is a taxon-specific constant, independent of mass and temperature, and α ≈ 3/4 ( see [51][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 system 51,55-59 . The scale parameter, b 0 , is derived from the underlying biochemical kinetics of metabolism and, as such, depends on the temperature, T, at which an organism operates 57 . This is given by the exponential Arrhenius-Boltzmann factor, exp(−E B /kT) , where E B (~ 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 which is the central equation of metabolic theory of ecology 57,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 plants 49 . Temperature is one of the primary drivers of energy flux in ecosystems 62 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 time 63,64 . At steady state, about half of GPP is respired by plants to support production, maintenance, and ion uptake 65 . 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 ∝ 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 ≈ 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-availability 26,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 where β NPP is an exponent capturing the effect of precipitation on NPP and c 1 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 www.nature.com/scientificreports/ a sublinear constraint on NPP, β NPP = 0.72 . Figure 1 illustrates these relationships in a 3D surface plot of NPP as a function of temperature and precipitation.
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 species 26 . 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 result 68,69 . Nichebased 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 environments [70][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 warmer 26,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 generations [73][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.  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.
Scientific Reports | (2020) 10:19654 | https://doi.org/10.1038/s41598-020-76658-2 www.nature.com/scientificreports/ 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: 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 where β S = E S /E NPP = (E NPP + E I )/E NPP = 1 + δ S , and so δ S = E I /E NPP . Therefore, diversity will increase superlinearly ( β 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: 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. H 2 Diversity kinetics is faster than productivity kinetics. The kinetics of diversity, E s , should be faster than the kinetics of environmental productivity, E NPP , as species interaction rates, E I , are predicted to be both nonzero and increase with environmental productivity, NPP. H 3 Diversity is more than environmental productivity. Diversity, S, should increase with environmental productivity, P, at a superlinear rate given by β S = E S /E NPP > 1. H 4 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 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 ′ = S i /S max , for both mammal and cultural diversity and fit models of the form lnS ′ = lnS ′ 0 − αL , where α is the slope of the response and L is latitude. Results show cultural diversity decreases with absolute latitude at a rate α E = 0.05 ± 0.004 (95% C.I.) (OLS regression: d.f. = 73, r 2 = 83%, p < 0.0001) and mammalian diversity as α M = 0.05 ± 0.01 (OLS regression: d.f. = 68, r 2 = 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).
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 lnS ′ = lnS ′ 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, r 2 = 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, r 2 = 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 (E P = 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 β L = 1.   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. www.nature.com/scientificreports/

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   www.nature.com/scientificreports/ 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 β M = 1.23 and for human cultures β 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 rates 76 .
This biogeographic convergence of biological and cultural diversity is striking given the differences in the evolutionary mechanisms and time scales of biological and cultural evolution 77,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 widely [79][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-diversity 82 ).
While biological and cultural evolutionary processes might be similar in principle 35,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 scales 85,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 offspring 73,74,84 . Languages, cultures, and economies turnover at much faster time scales the biological species, from the near instantaneous to many thousands of years 77,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 − 88 . 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 degrees 89 . 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 basemap 90 .
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 R 91 . 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.