Mapping local and global variability in plant trait distributions

Edited by William H. Schlesinger, Cary Institute of Ecosystem Studies, Millbrook, NY, and approved October 18, 2017 (received for review May 31, 2017)
December 1, 2017
114 (51) E10937-E10946

Significance

Currently, Earth system models (ESMs) represent variation in plant life through the presence of a small set of plant functional types (PFTs), each of which accounts for hundreds or thousands of species across thousands of vegetated grid cells on land. By expanding plant traits from a single mean value per PFT to a full distribution per PFT that varies among grid cells, the trait variation present in nature is restored and may be propagated to estimates of ecosystem processes. Indeed, critical ecosystem processes tend to depend on the full trait distribution, which therefore needs to be represented accurately. These maps reintroduce substantial local variation and will allow for a more accurate representation of the land surface in ESMs.

Abstract

Our ability to understand and predict the response of ecosystems to a changing environment depends on quantifying vegetation functional diversity. However, representing this diversity at the global scale is challenging. Typically, in Earth system models, characterization of plant diversity has been limited to grouping related species into plant functional types (PFTs), with all trait variation in a PFT collapsed into a single mean value that is applied globally. Using the largest global plant trait database and state of the art Bayesian modeling, we created fine-grained global maps of plant trait distributions that can be applied to Earth system models. Focusing on a set of plant traits closely coupled to photosynthesis and foliar respiration—specific leaf area (SLA) and dry mass-based concentrations of leaf nitrogen (Nm) and phosphorus (Pm), we characterize how traits vary within and among over 50,000 50×50-km cells across the entire vegetated land surface. We do this in several ways—without defining the PFT of each grid cell and using 4 or 14 PFTs; each model’s predictions are evaluated against out-of-sample data. This endeavor advances prior trait mapping by generating global maps that preserve variability across scales by using modern Bayesian spatial statistical modeling in combination with a database over three times larger than that in previous analyses. Our maps reveal that the most diverse grid cells possess trait variability close to the range of global PFT means.
Modeling global climate and the carbon cycle with Earth system models (ESMs) requires maps of plant traits that play key roles in leaf- and ecosystem-level metabolic processes (14). Multiple traits are critical to both photosynthesis and respiration, foremost leaf nitrogen concentration (Nm) and specific leaf area (SLA) (57). More recently, variation in leaf phosphorus concentration (Pm) has also been linked to variation in photosynthesis and foliar respiration (712). Estimating detailed global geographic patterns of these traits and corresponding trait–environment relationships has been hampered by limited measurements (13), but recent improvements in data coverage (14) allow for greater detail in spatial estimates of these key traits.
Previous work has extrapolated trait measurements across continental or larger regions through three methodologies: (i) grouping measurements of individuals into larger categories that share a set of properties [a working definition of plant functional types (PFTs)] (4, 15), (ii) exploiting trait–environment relationships (e.g., leaf Nm and mean annual temperature) (1, 1620), or (iii) restricting the analysis to species whose presence has been widely estimated on the ground (2124). Each of these methods has limitations—for example, trait–environment relationships do not well explain observed trait spatial patterns (1, 25), while species-based approaches limit the scope of extrapolation to only areas with well-measured species abundance. More critically, the first two global methodologies emphasized estimating a single trait value per PFT at every location, whereas both ground-based (5, 14) and remotely sensed (26) observations suggest that at ecosystem or landscape scales traits would be better represented by distributions. Here, we use an updated version of the largest global database of plant traits (14) coupled with modern Bayesian spatial statistical modeling techniques (27) to capture local and global variability in plant traits. This combination allows the representation of trait variation both within pixels on a gridded land surface and across global environmental gradients.
Information is lost when the range of measured trait values is compressed into a single PFT (Fig. 1). We observe that the global range of site-level SLA values for a single PFT such as broadleaf evergreen tropical trees (Fig. 1 A and C) is quite large (2.7–65.2 m2⋅kg−1). Even after limiting the scope to a single well-measured 0.5° × 0.5° pixel within Panama (Fig. 1 B and C), there is still a wide range of SLA values (4.7–37.7 m2⋅kg−1) with a local mean of 15.7 m2⋅kg−1 and a local standard deviation of 5.4 m2⋅kg−1—over one-third of the local mean. By contrast, the mean SLA value of all species associated with broadleaf evergreen tropical trees is 13.9 m2⋅kg−1, over 10% lower than the local average (Fig. 1C). Thus, single trait values per PFT fail to capture variability in trait values within or among grid cells, i.e., over a wide range of spatial scales.
Fig. 1.
Trait data. (A) Global locations and values of specific leaf area measurements for the PFT tropical broadleaf evergreen trees. (B) Locations and values of specific leaf area measurements for the tropical broadleaf evergreen trees in Panama. The square in the center indicates a 0.5°×0.5° pixel containing the Barro Colorado Island sites (see Fig. 5). These points have been jittered up to 0.05° to highlight the density of measurements. (C) The full distribution of specific leaf area values for all species classified as evergreen broadleaf tropical trees. The blue line is the global data while black is the local pixel, and the dashed vertical lines are the respective means.
Transitioning from a single trait value per PFT (within or among grid cells) to a distribution may lead to significantly different modeling results (20) as critical plant processes, such as photosynthesis, are nonlinear with respect to these traits (28). This is reinforced by recent modeling studies that have begun to incorporate distributions of traits at regional (29, 30) and global (31) scales. It has been shown that using trait distributions leads to different estimates of carbon dynamics (32) and that higher-order moments of trait distributions contribute to sustaining multiple ecosystem functions (33). While species-level mapping (21, 23, 24) does capture trait distributions, it has been limited geographically and restricted to subsets of functional groups.
Even the largest plant trait database offers only partial coverage across the globe in terms of site-level measurements. Hence, gap-filling approaches need to be adopted to extrapolate trait values at regions with no data coverage. Here, we overcome data limitations through PFT classification, trait–environment relationships, and additional location information to develop a suite of models capable of estimating trait distributions across the entire vegetated globe. The simplest one is a categorical model, which assigns traits to maps of remotely sensed PFTs. Every species, with its corresponding trait values, is associated with a PFT and these trait distributions are extrapolated to the satellite-estimated range of the PFT (SI Appendix, Figs. S1 and S2). The second one is a Bayesian linear model that complements the PFT information with trait–environment relationships. The third one is a Bayesian spatial model that, in addition to PFTs and the trait–environment relationships, leverages additional location information via Gaussian processes (Materials and Methods). The use of a spatial Gaussian process in this context is unique and model evaluation reveals the superior predictive performance of this model.
Each of these methods interpolates (and extrapolates) both mean trait values and entire trait distributions across space (i.e., across grid cells on a global map). These models are further stratified by three different levels of PFT categorization: (i) PFT-free, all plants in a single group (i.e., no PFTs); (ii) broad, 4 groups based on growth form and leaf type; and (iii) narrow, 14 groups based on further environmental, phenological, and photosynthetic categories (Materials and Methods). The PFT-free categorization groups all plants into a single class, while the broad grouping (4-PFT) is similar to the vegetation classification used in the Joint UK Land Environment Simulator land surface (34), and the narrow (14-PFT) category is equivalent to the classification used in the Community Land Model (CLM) (4, 15, 35).
The abovementioned methods allow for a representation of global vegetation that enables a more accurate formulation of functional diversity than the single-trait value per PFT paradigm that is widely used (4). The traits studied here—SLA, Nm, and Pm—are central to predicting variation in rates of plant photosynthesis (5, 6, 9, 11) and foliar respiration (10, 36). The importance of these traits and the more advanced representation of functional diversity developed here may be used to better capture the response of the land surface component of the Earth system to environmental change.

Results and Discussion

Model Evaluation.

Given the full suite of nine models proposed, we conducted extensive model evaluation (Table 1) to determine the trade-offs associated with each methodology and resolution of PFT. We assessed the predictive capability of the models, using the root-mean-square predictive error (RMSPE) based on out-of-sample data (SI Appendix, section S6). Among the nine models, the spatial narrow 14-PFT model emerged as the best predictor of mean trait values for SLA and Nm and the second best for Pm (Table 1). However, the spatial broad 4-PFT model performed nearly as well (Table 1). The models’ abilities to correctly estimate the spread of the trait distributions were assessed using the out-of-sample coverage probabilities (CPs)—the proportion of instances the model-predicted 95% confidence intervals contained the observed trait values. Most of the models provided adequate coverage (CP of around 90% or more). See SI Appendix, section S4, for more detailed definitions of the model comparison metrics.
Table 1.
Model evaluation
Modelps-R2, %RMSPECP, %
SLA   
 CfNA8.1391.2
 Cb16.97.1394.7
 Cn26.06.6695.8
 Lf4.67.9991.3
 Lb23.46.9394.0
 Ln30.76.5395.2
 Sf45.57.5493.6
 Sb58.56.3197.7
 Sn60.26.1397.7
Nm   
 CfNA7.1693.3
 Cb12.56.9593.2
 Cn19.46.4792.7
 Lf5.27.2893.2
 Lb16.76.7194.3
 Ln24.16.4294.6
 Sf44.27.1993.6
 Sb53.76.3696.1
 Sn54.86.1896.1
Pm   
 CfNA0.8690.5
 Cb5.30.8690.5
 Cn28.10.7891.1
 Lf25.60.8487.2
 Lb32.80.8585.3
 Ln35.40.8287.0
 Sf62.00.8390.7
 Sb66.70.8192.0
 Sn67.60.8091.3
Shown are the pseudo-R2 (ps-R2), RMSPE, and CP statistics for all nine models, for each of the three traits. The entries in boldface type correspond to the model producing highest ps-R2, lowest RMSPE, or CP closest to 0.95. The categorical PFT-free model (Cf) produces a constant estimate and hence ps-R2 is not defined. Each model is indicated by a two-letter abbreviation: C, categorical (no regression); L, linear (linear regression); and S, spatial (linear regression with spatial term) and the accompanying PFT resolution: f, PFT-free (no PFT information); b, broad (4-PFT); and n, narrow (14-PFT).
The improvement in prediction afforded by the inclusion of (i) a spatial term and (ii) PFT information (Table 1) invites further examination. First, the spatial term in our model likely incorporates some of the finer-scale variation that is unavailable given the relatively large grid cell size of the environmental covariates used in global studies. Thus, the spatial term allows for adjustment of trait values among neighboring or regional grid cells that the relatively coarse environmental metrics are not able to capture. Finer-scale studies that can evaluate local variations in climate, soil, or other relevant abiotic or biotic covariates may see less improvement from the inclusion of a spatial term, as they may directly measure local sources of variation. Second, the use of PFTs greatly improves the models, perhaps for similar reasons involving the degree of variation the raw data fail to incorporate. The greatest decrease in RMSPE occurs between the PFT-free grouping (a single category for all plants) and the broad (4-PFT) grouping across each of the models tested. If our trait data were perfectly predicted by environment, there would be no usefulness to including PFTs in mapping traits. That this is not so implies that the broad PFTs, based primarily on growth form and leaf type, offer superior predictive skill than environmental covariates on their own (19). However, the extra information in the narrow (14-PFT) grouping does further improve the fit and produces the most accurate predicted trait surface.

Global Maps.

We selected two sets of maps to describe, in broad strokes, how trait distributions vary across the land surface: the narrow 14-PFT spatial model and its categorical counterpart. The narrow 14-PFT spatial model is the best predictor of mean trait values and provided adequate coverage probability (Figs. 2 A and B, 3 A and B, and 4 A and B). For comparison, we also include the 14-PFT categorical model, which is most similar to maps currently used in ESMs (Figs. 2 C and D, 3 C and D, and 4 C and D). Maps for the other models can be found in SI Appendix, Figs. S8–S16. The mean and SD are presented as a summary of the full log-normal distribution within each pixel, but there are full distributions estimated in each pixel (Case Studies).
Fig. 2.
SLA maps. (A and B) Narrow (14-PFT) Bayesian spatial model pixel mean and SD estimates, respectively. (C and D) Narrow (14-PFT) categorical model pixel mean estimates and SD estimates, respectively. For clarity, the color bars have been truncated at the compound 5th and 95th percentiles of both models. Latitude tick marks indicate the equator, tropics, and Arctic Circle and longitude is marked at 100°W, 0°, and 100°E.
Fig. 3.
Nitrogen (mass) maps. (A and B) Narrow (14-PFT) Bayesian spatial model pixel mean and SD estimates, respectively. (C and D) Narrow (14-PFT) categorical model pixel mean estimates and SD estimates, respectively. For clarity, the color bars have been truncated at the compound 5th and 95th percentiles of both models. Latitude tick marks indicate the equator, tropics, and Arctic Circle and longitude is marked at 100°W, 0°, and 100°E.
Fig. 4.
Phosphorus (mass) maps. (A and B) Narrow (14-PFT) Bayesian spatial model pixel mean and SD estimates, respectively. (C and D) Narrow (14-PFT) categorical model pixel mean estimates and SD estimates, respectively. For clarity, the color bars have been truncated at the compound 5th and 95th percentiles of both models. Latitude tick marks indicate the equator, tropics, and Arctic Circle and longitude is marked at 100°W, 0°, and 100°E.
The SD maps (Figs. 2 B and D, 3 B and D, and 4 B and D) compared with the mean maps (Figs. 2 A and C, 3 A and C, and 4 A and C) highlight one of the central results of this analysis—the local SDs of trait values are of similar magnitudes to their respective means. Generally, we observed that the local SD is close to half the local mean value but can approach the global range of the trait mean values; e.g., Nm (Fig. 3) has a maximum local SD of 9 mg N/g, and the global mean range is only ∼10 mg N/g. The maps of the trait SDs follow similar patterns to the means, although there are several regions where the mean varies more markedly than the SD, such as SLA in the southeast United States and China in the categorical model (Fig. 2 C and D) and similarly for Nm in the spatial model across the Sahel in sub-Saharan Africa (Fig. 3 A and C). The lack of variation in the SD is most clear in the categorical model for Nm while both models show relatively modest variation in Pm.
For each of the three traits, the broad features of both the categorical and spatial models are similar, but there are numerous marked differences across regional and fine spatial scales (Figs. 24). The shared broad features of the maps from both models include SLA (Fig. 2) and Pm (Fig. 4) increasing from the tropics to the poles, while Nm (Fig. 3) has more modest variation, except that it tends to be lower in regions dominated by needle-leaved trees. Some of the notable differences between the models include the spatial model’s greater range and more marked variability of SLA within equatorial regimes (e.g., Brazil or central Africa); it also captures the low SLA of most of arid Australia better than the categorical model (Fig. 2A) and more strongly highlights the gradient of Pm from the tropics to the Arctic (16) (Fig. 4A).
The most consistent estimates between the categorical and spatial models are in the boreal regions dominated by needle-leaved trees; the measurements in this region are relatively sparse, which may have limited the ability of the spatial model to capture differences. On the other hand, broad-leaved trees span a wide range of environments, but a large portion of the measurements come from the tropics (66%), where there is a limited range of values among the climate covariates and therefore little variation with which to estimate a correlation. The grasses and shrubs have the largest SDs of the four broad PFTs (SI Appendix, Table S4) and dominate wide swaths of the land surface, but have fewer measurements—shrubs are the least measured of the broad PFTs in the database, and this appears to reduce the accuracy of the categorical model more than that of the spatial model (Table 1). The fact that shrubs are assumed to dominate in arid and boreal environments, which also tend to be undersampled, also likely contributes to these differences.
Our results also suggest that the breadth of functional niche space is reduced in both boreal and tropical biogeographic regions. The low variation across all three traits within the boreal forest implies that there is strong filtering and smaller niche space available in this relatively harsh environment. Surprisingly, despite the high species diversity in tropical forests, we also find that SLA and Pm have relatively low variation in these forests—suggesting that in this environment the trait space is reduced. This could be, in part, an artifact of the Earth system model PFT classification omitting herbaceous species. Conversely, grasslands and savannahs exhibit large variation in total trait space, suggesting these environments permit a wider range of strategies than in both the boreal and tropical regions. Most broadly, both the data and the spatial model suggest (SI Appendix, Figs. S24 and S25) lowest leaf nitrogen values in temperate climates that increase in both cooler and warmer regions; this may indicate a more complicated leaf biochemistry–temperature relationship than has previously been suggested (16).

Case Studies.

We conducted two regional case studies to provide a more in-depth analysis of the true and predicted shapes of trait distributions than can be provided by the SD maps and coverage probability. In these case studies trait data were pooled over an area to construct full trait distributions and then formally compared with the model predicted distributions.
We considered two areas with substantially different environmental conditions to evaluate the trait distributions obtained from the spatial and categorical models. We chose a single pixel that contained a highly studied site with numerous measurements of tropical trees, Barro Colorado Island (BCI), Panama; and a collection of pixels in an arid environment in which the mean estimates for SLA of the spatial and categorical models substantially disagreed, the southwestern United States. These areas were in the training data, and this analysis constituted a more detailed analysis of the models’ fit to the observed distribution of these locations. Here, the focus was on the structure of the full distribution of traits predicted at these sites; SI Appendix, Fig. S17 is a map of the measurements that comprised these locations and other sites included in this analysis. Both areas offer further insight into the structure of the distributions estimated by the categorical and spatial models.
In the pixel containing BCI, the categorical and spatial models broadly agreed for all three traits (Fig. 5 A, C, and E), although the spatial model means were only half as distant from the observed means for SLA and Nm (4% vs. 8% and 5% vs. 10%, respectively). There were only two PFTs present in this pixel: tropical broadleaf evergreen and deciduous trees. Despite the general similarity of the shapes of the distribution, the spatial model appears capable of capturing some subtle features. This is clearest for leaf nitrogen, where the peak of the distribution was quite broad. This is neatly captured in the narrow PFT model, and the pattern was detectable through the Kolmogorov–Smirnov (K-S) statistic, which evaluates the similarity of two full distributions. Indeed, the superiority of the spatial model was reinforced by a closer match for the Bayesian spatial model across all traits at BCI, although for Pm it was the PFT-free spatial model that fitted best (SI Appendix, Table S6).
Fig. 5.
Empirical trait distributions. Barro Colorado Island (A, C, and E) and the US Southwest (B, D, and F). A and B show SLA, C and D show leaf nitrogen, and E and F show leaf phosphorus. Each panel depicts the distribution of the data in solid black, the categorical model in blue, and the Bayesian spatial model in red. The dashed vertical lines indicate mean values.
The differences between the trait distributions of the categorical and Bayesian spatial models were stark in the southwestern United States, although the mean estimates for Nm and Pm were close (Fig. 5 B, D, and F). This may be a result of the topographic complexity of this region and the resulting difficulty of aggregating climate and soil covariates at the 0.5° pixel scale and the sparser sampling than at BCI. To get enough data to approximate a distribution, we aggregated 18 pixels with nine PFTs including every temperate category, although many of them are only marginally present. The inclusion of so many PFTs produced a noisier distribution in the categorical model than suggested by the data and estimated by the spatial model. Neither of the models produced distributions that matched as well with the observations; however, it is notable how close the mean values for both models matched the observations for Nm and Pm, and the spatial model did well for the mean SLA.

Environmental Covariates and the Spatial Term.

The improvement in prediction from the linear model to the spatial model is partially explained by weak trait–environment relationships (SI Appendix, Tables S1–S3). The magnitude of spatial variation explained by the Gaussian process model is comparable to that of the unexplained trait variation. For most of the spatial models, the estimated spatial range was around 300 km; this suggests a strong spatial effect and implies that the spatial model can provide more precise information about the trait distribution near the locations where we have data. This was largely borne out in the case studies and is illustrated more explicitly in Fig. 6 where the predicted trait SD for the spatial model was up to 50% lower than for the linear nonspatial model near locations with trait measurements. The spatial model leverages local information to reduce the uncertainty of trait estimation near data locations and may provide guidance for future data collection by identifying high-uncertainty regions.
Fig. 6.
Spatial learning. (A) The spatial model SD of Nm. The predicted variation near the data locations (black dots) is much lower than variation at locations away from any data point. (B) The linear model SD, which does not account for local spatial information, has no such pattern.

Applications for Trait Distributions.

Plant traits vary across a range of spatial scales, and the spatial model best captures changes across large spatial gradients (such as in Amazonia and Australia) as well as the subtleties within pixels. Maps for all of the models highlight how much information about local variability is lost when representing plant traits with a single value and suggest that a first application of these maps will be for ESMs to incorporate these scales of variability. For process-based ESMs, the simplest model to incorporate will likely be the categorical model as it is closest to the current PFT approach, but this model is also the least flexible. The more sophisticated models developed here provide more accurate large-scale variation and may be used to infer new trait values in a novel climate by perturbing the climatic covariates (37). However, given the likelihood of nonlinear trait–environment relationships, the spatial sparsity of the data, and the possibility of alternate strategies within a PFT that may alter the trait–environment relationship in a future climate some caution is called for when using these models for extrapolation. Future ecosystem models could also integrate the leaf-level variation in these maps with canopy-scale changes in leaf display traits—leaf angle, azimuth, and total area.
We have emphasized the quality of the Bayesian spatial model with narrow PFTs, but there is an intriguing possibility opened by the PFT-free model (SI Appendix, Figs. S8, S11, and S14)—that being the representation of vegetation without reference to PFTs (1). In this case the representation of vegetation would rely entirely on the structure of trait distributions at various landscape scales (1). Such a representation eliminates the need to separately model the future locations of PFTs (or species) when inferring the future distribution of traits; hence, the output of a model like that developed here could be updated with future environmental covariates, with the caveats that “out of sample prediction” may entail. At the same time, this method would allow for greater functional diversity than multiple PFTs with single-trait values, as is currently used in most ESMs. Adopting this approach does, however, raise the issue of how to deal with the paucity of surface observations in some regions, as evidenced by the greater errors associated with estimating out of sample values with this model (Table 1). Complementary work has retrieved leaf trait maps from a global carbon cycle model fused with Earth observations (38), providing another method that could be used for direct comparison against the trait maps produced here. While the methodology outlined in our analysis brings the possibility of a PFT-free land surface closer, we remain several steps away from being able to make such maps as accurately as we do using PFT characterizations for trait prediction. Several actions can bring us closer to that goal. First, incorporation of additional information (such as phylogenetic relatedness and trait–trait covariance) will likely improve trait maps, even using existing observations. Second, as the current level of observations is extremely sparse in some regions and sparse in most regions, expanded trait databases will also aid in development of PFT-free trait maps.

Conclusions

SLA and Nm are essential inputs into the land surface components of Earth system models, and while phosphorus has not yet been as widely incorporated into ESMs, it has been shown—particularly across the tropics—to be important to photosynthesis (9, 11, 3942) and respiration (11, 12, 36). The maps and trait–environment relationships presented here may be used by existing land surface models that use similar categories to classify vegetation. However, it should be noted that PFT-dependent models often have many other parameters that have been calibrated to historical estimates of particular trait values (4). Thus, the values developed here, while likely drawing from a larger pool of measurements than has been done previously, cannot necessarily be adopted without further modification of other model elements (37, 43). Nonetheless, these results can be incorporated into a wide class of models with relative ease. We can now provide global trait distributions at the pixel scale.
The global land surface is perhaps the most heterogeneous component of the Earth system. Reducing vegetation to a collection of PFTs with fixed trait values has been the preferred method to constrain this heterogeneity and group similar biochemical and biophysical properties; however, this has been at the expense of functional diversity. This analysis quantifies the substantial magnitude of this ignored trait variation. The approach and methods presented here retain the simplicity of the PFT representation, but capture a wider range of functional diversity.

Materials and Methods

Data.

The TRY database (www.try-db.org) (14) provided all data for leaf traits and the categorical traits to aggregate PFTs (TRY–Categorical Traits Dataset, https://www.try-db.org/TryWeb/Data.php#3, January 2016) used in the analysis. The TRY data may be requested from the TRY database custodians. See SI Appendix, section S10 for a complete list of the original publications associated with this subset of TRY. The extract from TRY used here has just under 45,000 measurements of individuals from 3,680 species with measurements of at least one of SLA, leaf nitrogen per dry leaf mass (Nm), and/or leaf phosphorus per leaf dry mass (Pm). The number of individual measurements varies from 32,315 for SLA on 2,953 species to 19,282 for Nm on 3,053 species down to 8,052 for Pm on 1,810 species; see SI Appendix, Table S4 for the number of unique measurements and species found in all categorizations used in the analysis. The species taxonomy was standardized using The Plant List (www.theplantlist.org/). Measurements were associated with environmental categories through Köppen–Geiger climate zones (44). All environmental variables are on a 0.5°×0.5° grid. Climate variables use 30-y climatologies from 1961 to 1990 as estimated by the Climate Research Unit (45, 46). Soil variables are from the International Soil Reference and Information Center–World Inventory of Soil Emission Potentials (ISRIC-WISE) (47). The spatial extent of PFTs has been previously estimated through satellite estimates of land cover around the year 2005 (48), and these estimates have been refined into climatic categories (15, 35). While TRY, and thus the data used here, represents the largest collection of plant traits in the world, most of the measurements come from a subset of global regions: North America, Europe, Australia, China, Japan, and Brazil. There are still large sections of the planet with extremely sparse measurements, notably much of the tropics outside of the Americas, large swaths of Central Asia, the Russian Federation, South Asia, and much of the Arctic (SI Appendix, Fig. S17). Improving data collection in these regions will greatly improve future modeling efforts. Until observations are more complete there remains the possibility of spurious patterns, although we have found little evidence to suggest their presence in this analysis, even in comparison with detailed regional studies (SI Appendix, Fig. S26) (49).

Classification of PFTs and Categorical Model.

We used three nested levels of PFT classification. In the first level, all plants are categorized into a single group (“PFT-free”). In the second level (“broad”), all plants are categorized into PFTs based on categorical traits associated with growth form (grass, shrub, tree) and leaf type (broad and needle-leaved), leading to the following four PFTs: grasses, shrubs, broad-leaved trees, and needle-leaved trees (SI Appendix, Fig. S1). In the third level (“narrow”), the broad PFTs are further refined by their climatic region—tropical, temperate, boreal—as well as leaf phenology and, for the grasses, photosynthetic pathway (C3 or C4). This produces 14 PFTs (SI Appendix, Fig. S2), which correspond exactly to those found in the CLM (4). Note that these PFT classifications exclude nonwoody eudicots (“herbs”), which were excluded from the analysis, on account of their lack of dominance within these PFT categories (50) and therefore, on account of being widely measured could overly influence the structure of the trait distributions if they were included. Satellite estimates of the PFT abundance that correspond to the narrow PFT categories defined above have already been calculated (15, 48) and we used these to assign a percentage of each 0.5°×0.5° pixel to each PFT present according to the fraction of the land surface within that pixel occupied by the PFT. The broad PFT fractions are calculated by summing the narrow PFT categories within each broad classification.
The categorical model uses the PFT categories and averages trait values for each species across individual measurements at each measured location. This defines the PFT as the interspecies range of trait values and ignores all local environmental factors. The results of the categorical model are summarized by the mean and SD of each PFT’s trait values (SI Appendix, Table S4) for all three resolutions of the model. Note that in the PFT-free case where no PFT information is used, the categorical model produces a constant trait distribution across the entire vegetated world. The categorical model and the Bayesian models described in the following section all use location-specific species mean values to estimate trait distributions. We assume no intraspecific variation in trait values. However, in regions dominated by a small number of species this may lead to biased predictions. The hyperdominance of a small group of species in the Amazon has recently been demonstrated (51) and thus serves as a case study to evaluate our assumption of equal species weighting (SI Appendix, section S8, Fig. S23). We found that equal weights (species means) produced trait distribution estimates closest to those of the hyperdominant trait abundances and this reinforces the use of this assumption globally. Further, as noted above, the omission of herbaceous species from tropical regions in this analysis (and ref. 51) may unduly limit trait diversity and calls for further research.

Bayesian Models.

A more fine-tuned depiction of geographical or spatial variation of plant trait values within each PFT can be achieved by leveraging environmental and location information, which allows trait values to adjust based on local conditions. Data for 17 climate- (45, 46) and soil-based (47) environmental predictors were available at the 0.5°×0.5°-pixel resolution used to create the trait maps. To avoid overfitting and collinearity issues, these 17 predictors were screened (SI Appendix, section S7) based on correlations among predictors, based on their individual correlation with the traits, and to include climate covariates along different axes of environmental stress and both chemical and physical soil covariates. We finally selected 5 predictors—mean annual temperature (MAT), total annual radiation (RAD), moisture index (precipitation/evapotranspiration) (MI), percentage of hydrogen (aqueous) (pH), and percentage of clay content (CLY). Remote-sensing data products, such as Normalized Difference Vegetation Index (52), are not used as covariates, to allow for inference outside of the historical observation period through perturbations of environmental covariates.
We used environment–trait relationships to obtain predictions of trait values (1, 1618, 37, 43) in a linear regression setup. The formal details of the initial model are as follows. We denote log-transformed trait values at a geographical location s as ytrait(s). This set of five predictors at a location s is denoted by the vector x(s)=(x1(s),x2(s),,x5(s)). A linear regression model relating the trait to the environmental predictors is specified as
ytrait(s)=b0+b1x1(s)+b2x2(s)++b5x5(s)+ϵ(s),
[1]
where bi are the regression coefficients and ϵ(s) is the error term explaining residual variation. Estimation of model parameters and prediction were achieved with a fully Bayesian hierarchical model. This enables inclusion of prior information and prediction of full trait distributions instead of representative values (like mean or median), thereby ensuring that the uncertainty associated with the estimation of model parameters is fully propagated into the predictive trait distributions.
We then generalized the above model into a Bayesian spatial linear regression model that borrows information from geographically proximal regions to capture residual spatial patterns beyond what is explained by environmental predictors. A customary specification of a spatial regression model is obtained by splitting up the error term ϵ(s) in Eq. 1 into the sum of a spatial process w(s) and an error term η(s) that accounts for the residual variation after adjusting for the spatial effects w(s). The underlying latent process w(s) accounts for local nuances beyond what is captured by the environmental predictors and is often interpreted as the net contribution from unobserved or unusable predictors. Gaussian processes (GPs) are widely used for modeling unknown spatial surfaces such as w(s), due to their convenient formulation as a multivariate Gaussian prior for the spatial random effect, unparalleled predictive performance (53), and ease of generating uncertainty-quantified predictions at unobserved locations. We use the computationally effective nearest-neighbor GP (27), which nicely embeds into the Bayesian hierarchical setup as a prior for w(s) in the second stage of the model specification. All technical specifications of the Bayesian spatial model are provided in SI Appendix, section S1.
The linear regression models used in previous studies (1, 1618) and both the spatial and nonspatial Bayesian models described above assume a global relationship between the traits and environment. Given the goal of predicting trait values for the entire land surface, the assumption of a universal trait–environment relationship may be an oversimplification (54). Moreover, if there is significant variation in plant trait values among different PFTs, the estimated parameters will be skewed toward values from abundantly sampled PFTs, such as broad-leaved trees. Additional information about plant characteristics at a specific location, if available, can potentially be used to improve predictions. As mentioned earlier, we have PFT classifications for each observation of the dataset used here and satellite estimates of PFT abundance at all pixels. The global regression approaches described above ignore this information and can yield biased predictions at locations dominated by PFTs poorly represented in the data, such as shrubs. Hence, we also incorporate the PFT information in these regression models by allowing the trait–environment relationship to vary between different PFTs. Finally, the PFT-specific distributions from the Bayesian models were weighted by the satellite-based PFT abundances to create a landscape-scale trait distribution, thereby enabling straightforward comparison between all three categorizations of PFT. Details of the PFT-based Bayesian models are provided in SI Appendix, section S2. The use of a GP-based spatial model as well as the Bayesian implementation of the regression models was unique to this application of plant trait mapping and, as results indicated, were critical to improving model predictions as well as properly quantifying trait distributions.

Data Availability

Data deposition: The code and data necessary to run the models are available at https://github.com/abhirupdatta/global_maps_of_plant_traits.

Acknowledgments

The authors appreciate the improvements suggested by two anonymous referees, which improved the clarity and depth of the manuscript. This research was supported as part of the Energy Exascale Earth System Model (E3SM) project, funded by the US Department of Energy, Office of Science, Office of Biological and Environmental Research (Grant DE-SC0012677 to E.E.B, H.F.M., M.C., K.R.W., A.B., and P.B.R.). O.K.A. acknowledges the support of the Australian Research Council (CE140100008). This research was also funded by programs from the NSF Long-Term Ecological Research (Grant DEB-1234162) and Long-Term Research in Environmental Biology (Grant DEB-1242531). A.B., F.F., and P.B.R. acknowledge funding from NSF Grant IIS-1563950. P.B.R. also acknowledges support from two University of Minnesota Institute on the Environment discovery grants. This study has been supported by the TRY initiative on plant traits (www.try-db.org). The TRY database is hosted at the Max Planck Institute for Biogeochemistry (Jena, Germany) and supported by DIVERSITAS/Future Earth, the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, and BACI (Towards a Biosphere Atmosphere Change Index) (Grant 640176). B.B. acknowledges a Natural Environment Research Council (NERC) independent research fellowship NE/M019160/1. J.P. acknowledges the financial support from the European Research Council Synergy Grant ERC-SyG-2013-610028 IMBALANCE-P, the Spanish Government Grant CGL2013-48074-P, and the Catalan Government Grant SGR 2014-274. B.B.-L. was supported by the Earth System Modeling program of the US Department of Energy, Office of Science, Office of Biological and Environmental Research. K.K. acknowledges the contribution of the Wageningen University and Research Investment theme Resilience for the project Resilient Forest (KB-29-009-003). P.M. acknowledges support from ARC Grant FT110100457 and NERC Grant NE/F002149/1. W.H. acknowledges support from the National Natural Science Foundation of China (Grant 41473068) and the “Light of West China” Program of the Chinese Academy of Sciences.

Supporting Information

Appendix (PDF)

References

1
PM Van Bodegom, JC Douma, LM Verheijen, A fully traits-based approach to modeling global vegetation distribution. Proc Natl Acad Sci USA 111, 13733–13738 (2014).
2
V Maire, et al., Global effects of soil and climate on leaf photosynthetic traits and rates. Glob Ecol Biogeogr 24, 706–717 (2015).
3
RS DeFries, et al., Mapping the land surface for global atmosphere-biosphere models: Toward continuous distributions of vegetation’s functional properties. J Geophys Res 100, 20867 (1995).
4
GB Bonan, et al., Improving canopy processes in the Community Land Model version 4 (CLM4) using global flux fields empirically inferred from FLUXNET data. J Geophys Res 116, 1–22 (2011).
5
PB Reich, DS Ellsworth, MB Walters, Leaf structure (specific leaf area) modulates photosynthesis–nitrogen relations : Evidence from within and across species and functional groups. Funct Ecol 12, 948–958 (1998).
6
J Kattge, W Knorr, T Raddatz, C Wirth, Quantifying photosynthetic capacity and its relationship to leaf nitrogen content for global-scale terrestrial biosphere models. Glob Change Biol 15, 976–991 (2009).
7
KY Crous, et al., Nitrogen and phosphorus availabilities interact to modulate leaf trait scaling relationships across six plant functional types in a controlled-environment study. New Phytol 215, 992–1008 (2017).
8
IJ Wright, et al., The worldwide leaf economics spectrum. Nature 428, 821–827 (2004).
9
PB Reich, J Oleksyn, IJ Wright, Leaf phosphorus influences the photosynthesis-nitrogen relation: A cross-biome analysis of 314 species. Oecologia 160, 207–212 (2009).
10
OK Atkin, et al., Global variability in leaf respiration in relation to climate, plant functional types and leaf traits. New Phytol 206, 614–636 (2015).
11
N Bahar, et al., Leaf-level photosynthetic capacity in lowland Amazonian and high elevation, Andean tropical moist forests of Peru. New Phytol 214, 1002–1018 (2016).
12
L Rowland, et al., Scaling leaf respiration with nitrogen and phosphorus in tropical forests across two continents. New Phytol 214, 1064–1077 (2016).
13
PB Reich, Global biogeography of plant chemistry: Filling in the blanks. New Phytol 168, 263–266 (2005).
14
J Kattge, et al., TRY - a global database of plant traits. Glob Change Biol 17, 2905–2935 (2011).
15
KW Oleson, et al., Technical description of version 4.5 of the Community Land Model (CLM) (National Center for Atmospheric Research, Boulder, CO), Technical Report NCAR/TN-503+STR. (2013).
16
PB Reich, J Oleksyn, Global patterns of plant leaf N and P in relation to temperature and latitude. Proc Natl Acad Sci USA 101, 11001–11006 (2004).
17
JC Ordoñez, et al., A global study of relationships between leaf traits, climate and soil measures of nutrient fertility. Glob Ecol Biogeogr 18, 137–149 (2009).
18
AH Simpson, SJ Richardson, DC Laughlin, Soil-climate interactions explain variation in foliar, stem, root and reproductive traits across temperate forests. Glob Ecol Biogeogr 25, 964–978 (2016).
19
PB Reich, IJ Wright, CH Lusk, Predicting leaf physiology from simple plant and climate attributes: A global GLOPNET analysis. Ecol Appl 17, 1982–1988 (2007).
20
PB Reich, RL Rich, X Lu, YP Wang, J Oleksyn, Biogeographic variation in evergreen conifer needle longevity and impacts on boreal forest carbon cycle projections. Proc Natl Acad Sci USA 111, 13703–13708 (2014).
21
NG Swenson, et al., The biogeography and filtering of woody plant functional diversity in North and South America. Glob Ecol Biogeogr 21, 798–808 (2012).
22
BA Hawkins, M Rueda, TF Rangel, R Field, JAF Diniz-Filho, Community phylogenetics at the biogeographical scale: Cold tolerance, niche conservatism and the structure of North American forests. J Biogeogr 41, 23–38 (2014).
23
I Šímová, et al., Shifts in trait means and variances in North American tree assemblages: Species richness patterns are loosely related to the functional space. Ecography 38, 649–658 (2015).
24
NG Swenson, et al., Phylogeny and the prediction of tree functional diversity across novel continental settings. Glob Ecol Biogeogr 26, 553–562 (2017).
25
JC Douma, MWA de Haan, R Aerts, JPM Witte, PM van Bodegom, Succession-induced trait shifts across a wide range of NW European ecosystems are driven by light and modulated by initial abiotic conditions. J Ecol 100, 366–380 (2012).
26
GP Asner, DE Knapp, CB Anderson, RE Martin, N Vaughn, Large-scale climatic and geophysical controls on the leaf economics spectrum. Proc Natl Acad Sci USA 113, E4043–E4051 (2016).
27
A Datta, S Banerjee, A Finley, A Gelfand, Hierarchical nearest-neighbor Gaussian process models for large geostatistical datasets. J Am Stat Assoc 111, 800–812 (2016).
28
GD Farquhar, S von Caemmerer, JA Berry, A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149, 78–90 (1980).
29
S Scheiter, SI Higgins, Impacts of climate change on the vegetation of Africa: An adaptive dynamic vegetation modelling approach. Glob Change Biol 15, 2224–2246 (2009).
30
S Scheiter, L Langan, SI Higgins, Next-generation dynamic global vegetation models: Learning from community ecology. New Phytol 198, 957–969 (2013).
31
R Pavlick, DT Drewry, K Bohn, B Reu, A Kleidon, The Jena Diversity-Dynamic Global Vegetation Model (JeDi-DGVM): A diverse approach to representing terrestrial biogeography and biogeochemistry based on plant functional trade-offs. Biogeosciences Discussions 9, 4627–4726 (2012).
32
C Pappas, S Fatichi, P Burlando, Terrestrial water and carbon fluxes across climatic gradients: Does plant diversity matter? New Phytol 16, 3663 (2014).
33
N Gross, et al., Functional trait diversity maximizes ecosystem multifunctionality. Nat Ecol Evol 1, 0132 (2017).
34
DB Clark, et al., The Joint UK Land Environment Simulator (JULES), model description – Part 2: Carbon fluxes and vegetation dynamics. Geosci Model Dev 4, 701–722 (2011).
35
GB Bonan, Landscapes as patches of plant functional types: An integrating concept for climate and ecosystem models. Glob Biogeochem Cycles 16, 5.1–5.18 (2002).
36
P Meir, J Grace, AC Miranda, Leaf respiration in two tropical rainforests: Constraints on physiology by phosphorus, nitrogen and temperature. Funct Ecol 15, 378–387 (2001).
37
LM Verheijen, et al., Inclusion of ecologically based trait variation in plant functional types reduces the projected land carbon sink in an earth system model. Glob Chang Biol 21, 3074–3086 (2015).
38
AA Bloom, JF Exbrayat, IR van der Velde, L Feng, M Williams, The decadal state of the terrestrial carbon cycle: Global retrievals of terrestrial carbon allocation, pools, and residence times. Proc Natl Acad Sci USA 113, 1285–1290 (2016).
39
P Meir, PE Levy, J Grace, PG Jarvis, Photosynthetic parameters from two contrasting woody vegetation types in West Africa. Plant Ecol 192, 277–287 (2007).
40
TF Domingues, et al., Co-limitation of photosynthetic capacity by nitrogen and phosphorus in West Africa woodlands. Plant Cell Environ 33, 959–980 (2010).
41
Q Zhang, YP Wang, AJ Pitman, YJ Dai, Limitations of nitrogen and phosphorous on the terrestrial carbon uptake in the 20th century. Geophys Res Lett 38, 1–5 (2011).
42
B Medlyn, et al., Using models to guide field experiments: A priori predictions for the CO2 response of a nutrient- and water-limited native Eucalypt woodland. Glob Change Biol 22, 2834–2851 (2016).
43
LM Verheijen, et al., Impacts of trait variation through observed trait-climate relationships on performance of an earth system model: A conceptual analysis. Biogeosciences 10, 5497–5515 (2013).
44
B Peel, BL Finlayson, TA McMahon, Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci 11, 1633–1644 (2007).
45
M New, M Hulme, P Jones, Representing twentieth-century space–time climate variability. Part I: Development of a 1961–90 mean monthly terrestrial climatology. J Clim 12, 829–856 (1999).
46
I Harris, PD Jones, TJ Osborn, DH Lister, Updated high-resolution grids of monthly climatic observations - the CRU TS3.10 Dataset. Int J Climatol 34, 623–642 (2014).
47
NH Batjes, ISRIC-WISE global data set of derived soil properties on a 0.5 by 0.5 degree grid (Version 3.0) (World Soil Information, Wageningen, The Netherlands), Report 2005/08. (2005).
48
PJ Lawrence, TN Chase, Representing a new MODIS consistent land surface in the Community Land Model (CLM 3.0). J Geophys Res 112, G01023 (2007).
49
GP Asner, et al., Airborne laser-guided imaging spectroscopy to map forest trait diversity and guide conservation. Science 355, 385–389 (2017).
50
DJ Gibson Grasses & Grassland Ecology (Oxford Univ Press, New York, 2009).
51
H ter Steege, et al., Hyperdominance in the Amazonian tree flora. Science 342, 1243092 (2013).
52
SV Ollinger, et al., Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests: Functional relations and potential climate feedbacks. Proc Natl Acad Sci USA 105, 19336–19341 (2008).
53
C Rasmussen, Evaluation of Gaussian processes and other methods for non-linear regression. PhD thesis (University of Toronto, Toronto). (1996).
54
LM Verheijen, R Aerts, G Bönisch, J Kattge, PM Van Bodegom, Variation in trait trade-offs allows differentiation among predefined plant functional types: Implications for predictive ecology. New Phytol 209, 563–575 (2016).

Information & Authors

Information

Published in

The cover image for PNAS Vol.114; No.51
Proceedings of the National Academy of Sciences
Vol. 114 | No. 51
December 19, 2017
PubMed: 29196525

Classifications

Data Availability

Data deposition: The code and data necessary to run the models are available at https://github.com/abhirupdatta/global_maps_of_plant_traits.

Submission history

Published online: December 1, 2017
Published in issue: December 19, 2017

Keywords

  1. plant traits
  2. Bayesian modeling
  3. spatial statistics
  4. global
  5. climate

Acknowledgments

The authors appreciate the improvements suggested by two anonymous referees, which improved the clarity and depth of the manuscript. This research was supported as part of the Energy Exascale Earth System Model (E3SM) project, funded by the US Department of Energy, Office of Science, Office of Biological and Environmental Research (Grant DE-SC0012677 to E.E.B, H.F.M., M.C., K.R.W., A.B., and P.B.R.). O.K.A. acknowledges the support of the Australian Research Council (CE140100008). This research was also funded by programs from the NSF Long-Term Ecological Research (Grant DEB-1234162) and Long-Term Research in Environmental Biology (Grant DEB-1242531). A.B., F.F., and P.B.R. acknowledge funding from NSF Grant IIS-1563950. P.B.R. also acknowledges support from two University of Minnesota Institute on the Environment discovery grants. This study has been supported by the TRY initiative on plant traits (www.try-db.org). The TRY database is hosted at the Max Planck Institute for Biogeochemistry (Jena, Germany) and supported by DIVERSITAS/Future Earth, the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, and BACI (Towards a Biosphere Atmosphere Change Index) (Grant 640176). B.B. acknowledges a Natural Environment Research Council (NERC) independent research fellowship NE/M019160/1. J.P. acknowledges the financial support from the European Research Council Synergy Grant ERC-SyG-2013-610028 IMBALANCE-P, the Spanish Government Grant CGL2013-48074-P, and the Catalan Government Grant SGR 2014-274. B.B.-L. was supported by the Earth System Modeling program of the US Department of Energy, Office of Science, Office of Biological and Environmental Research. K.K. acknowledges the contribution of the Wageningen University and Research Investment theme Resilience for the project Resilient Forest (KB-29-009-003). P.M. acknowledges support from ARC Grant FT110100457 and NERC Grant NE/F002149/1. W.H. acknowledges support from the National Natural Science Foundation of China (Grant 41473068) and the “Light of West China” Program of the Chinese Academy of Sciences.

Notes

This article is a PNAS Direct Submission.

Authors

Affiliations

Department of Forest Resources, University of Minnesota, St. Paul, MN 55108;
Abhirup Datta2,1 [email protected]
Department of Biostatistics, Johns Hopkins University, Baltimore, MD 21205;
Habacuc Flores-Moreno
Department of Forest Resources, University of Minnesota, St. Paul, MN 55108;
Department of Ecology, Evolution, and Behavior, University of Minnesota, St. Paul, MN 55108;
Ming Chen
Department of Forest Resources, University of Minnesota, St. Paul, MN 55108;
Kirk R. Wythers
Department of Forest Resources, University of Minnesota, St. Paul, MN 55108;
Farideh Fazayeli
Department of Computer Science and Engineering, University of Minnesota, Minneapolis, MN 55455;
Arindam Banerjee
Department of Computer Science and Engineering, University of Minnesota, Minneapolis, MN 55455;
Owen K. Atkin
Australian Research Council Centre of Excellence in Plant Energy, Research School of Biology, The Australian National University, Canberra, ACT 2601, Australia;
Division of Plant Sciences, Research School of Biology, The Australian National University, Canberra, ACT 2601, Australia;
Jens Kattge
Max Planck Institute for Biogeochemistry, 07745 Jena, Germany;
German Centre for Integrative Biodiversity Research Halle-Jena-Leipzig, 04103 Leipzig, Germany;
Bernard Amiaud
UMR 1137 Ecologie et Ecophysiologie Forestières, Université de Lorraine–Institut National de la Recherche Agronomique, 54506 Vandoeuvre-lès-Nancy, France;
Benjamin Blonder
Environmental Change Institute, University of Oxford, Oxford OX1 3BJ, United Kingdom;
Gerhard Boenisch
Max Planck Institute for Biogeochemistry, 07745 Jena, Germany;
Joint Global Change Research Institute, Department of Energy Pacific Northwest National Laboratory, College Park, MD 20740;
Kerry A. Brown
Department of Geography and Geology, Kingston University London, Surrey KT1 2EE, United Kingdom;
Chaeho Byun
School of Biological Sciences, Seoul National University, Seoul 08826, South Korea;
Giandiego Campetella
Plant Diversity and Ecosystems Management Unit, School of Biosciences & Veterinary Medicine, University of Camerino, 62032 Camerino, Italy;
Bruno E. L. Cerabolini
Department of Theoretical and Applied Sciences, University of Insubria, I-21100 Varese, Italy;
Johannes H. C. Cornelissen
Systems Ecology, Department of Ecological Science, Vrije Universiteit, 1081 HV Amsterdam, The Netherlands;
Jonah Ventures, Manhattan, KS 66502;
German Centre for Integrative Biodiversity Research Halle-Jena-Leipzig, 04103 Leipzig, Germany;
Department of Community Ecology, Helmholtz Centre for Environmental Research–UFZ, 06120 Halle (Saale), Germany;
Franciska T. de Vries
School of Earth and Environmental Sciences, The University of Manchester, Manchester M13 9PT, United Kingdom;
Sandra Díaz
Instituto Multidisciplinario de Biología Vegetal (Consejo Nacional de Invetigaciones Cientificas y Técnicas), Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba, CC 495 Córdoba, Argentina;
Departamento de Diversidad Biológica y Ecología, Facultad de Ciencias Exactas, Físicas y Naturales, Universidad Nacional de Córdoba, CC 495 Córdoba, Argentina;
Tomas F. Domingues
Faculdade de Filosofia Ciencias e Letras de Ribeirao Preto, Universidade de Sao Paulo, CEP 14040-901 Bairro Monte Alegre, Ribeirao Preto, Sao Paulo, Brazil;
Estelle Forey
Laboratory of Ecology Ecodiv, Institut National de Recherche en Sciences et Technologies pour l’Environnement et l’Agriculture, Normandie Université, 76821 Mont-Saint-Aignan, France;
Andrés González-Melo
Facultad de Ciencias Naturales y Matematicas, Universidad del Rosario, Bogota 110111, Colombia;
Departamento de Biología y Geología, Física y Química Inorgánica, Escuela Superior de Ciencias Experimentales y Tecnología, Universidad Rey Juan Carlos, 28933 Móstoles, Spain;
Institut National de la Recherche Agronomique, Unité Sous Contrat 1339, Centre d’Etude Biologique de Chizé, F 79360 Villiers en Bois, France;
Centre D’Étude Biologique de Chizé, CNRS–Université La Rochelle (UMR 7372), F-79360 Villiers en Bois, France;
Wenxuan Han
College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China;
Key Laboratory of Biogeography and Bio-Resource in Arid Land, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, Xinjiang, China;
Wesley N. Hattingh
School of Animal, Plant and Environmental Sciences, University of the Witwatersrand, WITS 2050, Johannesburg, South Africa;
Thomas Hickler
Senckenberg Biodiversity and Climate Research Centre (BiK-F), 60325 Frankfurt/Main, Germany;
Department of Physical Geography, Goethe-University, 60438 Frankfurt/Main, Germany;
Steven Jansen
Institute of Systematic Botany and Ecology, Ulm University, 89081 Ulm, Germany;
Koen Kramer
Team Vegetation, Forest and Landscape Ecology, Wageningen Environmental Research, 6708 PB Wageningen, The Netherlands;
Chairgroup Forest Ecology and Management, Wageningen University, 6708 PB Wageningen, The Netherlands;
Nathan J. B. Kraft
Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA 90095;
Hiroko Kurokawa
Department of Forest Vegetation, Forestry and Forest Products Research Institute, Tsukuba 305-8687, Japan;
Daniel C. Laughlin
Department of Botany, University of Wyoming, Laramie, WY 82071;
Patrick Meir
Division of Plant Sciences, Research School of Biology, The Australian National University, Canberra, ACT 2601, Australia;
School of Geosciences, University of Edinburgh, Edinburgh EH9 3FF, United Kingdom;
Vanessa Minden
Institute of Biology and Environmental Science, University of Oldenburg, 26111 Oldenburg, Germany;
Department of Plant Physiology, Estonian University of Life Sciences, 51014 Tartu, Estonia;
Yusuke Onoda
Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan;
CSIC, Unitat d’Ecologia Global CREAF-CSIC-UAB, Bellaterra 08193, Barcelona, Catalonia, Spain;
CREAF, Cerdanyola del Vallès 08193, Barcelona, Catalonia, Spain;
Quentin Read
Department of Forestry, Michigan State University, East Lansing, MI 48824;
Lawren Sack
Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA 90095;
Brandon Schamp
Department of Biology, Algoma University, Sault Ste. Marie, ON P6A 2G4, Canada;
Nadejda A. Soudzilovskaia
Institute of Environmental Sciences, Leiden University, 2333 CC Leiden, The Netherlands;
Marko J. Spasojevic
Department of Evolution, Ecology, and Organismal Biology, University of California, Riverside, CA 92521;
Enio Sosinski
Laboratório de Planejamento Ambiental, Embrapa Clima Temperado, Pelotas, RS, Brazil 96010-971;
Peter E. Thornton
Environmental Sciences Division and Climate Change Science Institute, Oak Ridge National Laboratory, Oak Ridge, TN 37831;
Fernando Valladares
Museo Nacional de Ciencias Naturales, CSIC, E-28006 Madrid Spain;
Peter M. van Bodegom
Institute of Environmental Sciences, Leiden University, 2333 CC Leiden, The Netherlands;
Mathew Williams
School of Geosciences, University of Edinburgh, Edinburgh EH9 3FF, United Kingdom;
Christian Wirth
Max Planck Institute for Biogeochemistry, 07745 Jena, Germany;
German Centre for Integrative Biodiversity Research Halle-Jena-Leipzig, 04103 Leipzig, Germany;
Department of Systematic Botany and Functional Biodiversity, University of Leipzig, 04103 Leipzig, Germany;
Peter B. Reich
Department of Forest Resources, University of Minnesota, St. Paul, MN 55108;
Hawkesbury Institute for the Environment, Western Sydney University, Penrith NSW 2751, Australia

Notes

2
To whom correspondence may be addressed. Email: [email protected] or [email protected].
Author contributions: E.E.B., A.D., H.F.-M., M.C., K.R.W., F.F., A.B., O.K.A., J.K., and P.B.R. designed research; E.E.B. and A.D. performed research; E.E.B., A.D., H.F.-M., and J.K. analyzed data; and E.E.B., A.D., H.F.-M., M.C., K.R.W., A.B., O.K.A., J.K., B.A., B.B., G.B., B.B.-L., K.A.B., C.B., G.C., B.E.L.C., J.H.C.C., J.M.C., D.C., F.T.d.V., S.D., T.F.D., E.F., A.G.-M., N.G., W.H., W.N.H., T.H., S.J., K.K., N.J.B.K., H.K., D.C.L., P.M., V.M., Ü.N., Y.O., J.P., Q.R., L.S., B.S., N.A.S., M.J.S., E.S., P.E.T., F.V., P.M.v.B., M.W., C.W., and P.B.R. wrote the paper.
1
E.E.B. and A.D. contributed equally to this work.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

Export the article citation data by selecting a format from the list below and clicking Export.

Cited by

    Loading...

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian

    Purchase options

    Purchase this article to access the full text.

    Single Article Purchase

    Mapping local and global variability in plant trait distributions
    Proceedings of the National Academy of Sciences
    • Vol. 114
    • No. 51
    • pp. 13303-E11060

    Figures

    Tables

    Media

    Share

    Share

    Share article link

    Share on social media