Intercomparison of global foliar trait maps reveals fundamental differences and limitations of upscaling approaches

Foliar traits such as specific leaf area (SLA), leaf nitrogen (N), and phosphorus (P) concentrations play important roles in plant economic strategies and ecosystem functioning. Various global maps of these foliar traits have been generated using statistical upscaling approaches based on in-situ trait observations. Here, we intercompare such global upscaled foliar trait maps at 0.5° spatial resolution (six maps for SLA, five for N, three for P), categorize the upscaling approaches used to generate them, and evaluate the maps with trait estimates from a global database of vegetation plots (sPlotOpen). We disentangled the contributions from different plant functional types (PFTs) to the upscaled maps and quantified the impacts of using different plot-level trait metrics on the evaluation with sPlotOpen: community weighted mean (CWM) and top-of-canopy weighted mean (TWM). We found that the global foliar trait maps of SLA and N differ drastically and fall into two groups that are almost uncorrelated (for P only maps from one group were available). The primary factor explaining the differences between these groups is the use of PFT information combined with remote sensing-derived land cover products in one group while the other group mostly relied on environmental predictors alone. The maps that used PFT and corresponding land cover information exhibit considerable similarities in spatial patterns that are strongly driven by land cover. The maps not using PFTs show a lower level of similarity and tend to be strongly driven by individual environmental variables. Upscaled maps of both groups were moderately correlated to sPlotOpen data aggregated to the grid-cell level (R = 0.2–0.6) when processing sPlotOpen in a way that is consistent with the respective trait upscaling approaches, including the plot-level trait metric (CWM or TWM) and the scaling to the grid cells with or without accounting for fractional land cover. The impact of using TWM or CWM was relevant, but considerably smaller than that of the PFT and land cover information. The maps using PFT and land cover information better reproduce the between-PFT trait differences of sPlotOpen data, while the two groups performed similarly in capturing within-PFT trait variation. Our findings highlight the importance of explicitly accounting for within-grid-cell trait variation, which has important implications for applications using existing maps and future upscaling efforts. Remote sensing information has great potential to reduce uncertainties related to scaling from in-situ observations to grid cells and the regression-based mapping steps involved in the upscaling.