Distinct latitudinal patterns and drivers of topsoil nitrogen and phosphorus across urban forests in eastern China.

Nitrogen (N) and phosphorus (P) are the two most important macronutrients supporting forest growth. Unprecedented urbanization has created growing areas of urban forests that provide key ecosystem services for city dwellers. However, the large-scale patterns of soil N and P content remain poorly understood in urban forests. Based on a systematic soil survey in urban forests from nine large cities across eastern China, we examined the spatial patterns and key drivers of topsoil (0-20 cm) total N content, total P content, and N:P ratio. Topsoil total N content was found to change significantly with latitude in the form of an inverted parabolic curve, while total P content showed an opposite latitudinal pattern. Variance partition analysis indicated that regional-scale patterns of topsoil total N and P contents were dominated by climatic drivers and partially regulated by time and pedogenic drivers. Conditional regression analyses showed a significant increase in topsoil total N content with lower mean annual temperature (MAT) and higher mean annual precipitation (MAP), while topsoil total P content decreased significantly with higher MAP. Topsoil total N content also increased significantly with the age of urban park and varied with pre-urban soil type, while no such effects were found for topsoil total P content. Moreover, topsoil N:P ratio showed a latitudinal pattern similar to that of topsoil total N content and also increased significantly with lower MAT and higher MAP. Our findings demonstrate distinct latitudinal trends of topsoil N and P contents and highlight a dominant role of climatic drivers in shaping the large-scale patterns of topsoil nutrients in urban forests.

Can ethylenediurea (EDU) alter the effects of ozone on the source-sink regulation of nitrogen uptake and remobilization during grain filling period in rice?

Increased surface ozone (O3) pollution seriously threatens crop production, and ethylenediurea (EDU) can alleviate crop yield reduction caused by O3. However, the reason for the decrease in grain nitrogen (N) accumulation caused by O3 and whether EDU serves as N fertilizer remain unclear. An experiment was conducted to investigate the impacts of factorial combinations of O3 enrichment (ambient air plus 60 ppb) and EDU (foliage spray with 450 ppm solutions) on N concentration, accumulation and remobilization in hybrid rice seedlings. Compared to ambient condition, elevated O3 significantly inhibited the N accumulation in vegetative organs during anthesis and grain N accumulation during the maturity stage. Elevated O3 significantly decreased the total N accumulation during anthesis and maturity stages, with a greater impact at the latter stage. The decrease in grain N accumulation caused by O3 was attributed to a decrease in N remobilization of vegetative organs during the grain filling period as well as to a decrease in post-anthesis N uptake. However, there was no significant change in the proportion of N remobilization and N uptake in grain N accumulation. The inhibitory effect of O3 on N remobilization in the upper canopy leaves was greater than that in the lower canopy leaves. In addition, elevated O3 increased the N accumulation of panicles at the anthesis stage, mainly by resulting in earlier heading of rice. EDU only increased N accumulation at the maturity stage, which was mainly attributed to an increase in rice biomass by EDU. EDU had no significant effect on N concentration, N remobilization process, and N harvest index. The findings are helpful to better understand the utilization of N fertilizer by rice under O3 pollution, and can also provide a theoretical basis for sustainable nutrient management to alleviate the negative impact of O3 on crop yield and quality.

Global distribution of surface soil organic carbon in urban greenspaces.

Urban greenspaces continue to grow with global urbanization. The global distribution and stock of soil organic carbon (SOC) in urban greenspaces remain largely undescribed and missing in global carbon (C) budgets. Here, we synthesize data of 420 observations from 257 cities in 52 countries to evaluate the global pattern of surface SOC density (0-20 cm depth) in urban greenspaces. Surface SOC density in urban greenspaces increases significantly at higher latitudes and decreases significantly with higher mean annual temperature, stronger temperature and precipitation seasonality, as well as lower urban greenness index. By mapping surface SOC density using a random forest model, we estimate an average SOC density of 55.2 (51.9-58.6) Mg C ha-1 and a SOC stock of 1.46 (1.37-1.54) Pg C in global urban greenspaces. Our findings present a comprehensive assessment of SOC in global urban greenspaces and provide a baseline for future urban soil C assessment under continuing urbanization.

Rapid migration of Mongolian oak into the southern Asian boreal forest.

The migration of trees induced by climatic warming has been observed at many alpine treelines and boreal-tundra ecotones, but the migration of temperate trees into southern boreal forest remains less well documented. We conducted a field investigation across an ecotone of temperate and boreal forests in northern Greater Khingan Mountains of northeast China. Our analysis demonstrates that Mongolian oak (Quercus mongolica), an important temperate tree species, has migrated rapidly into southern boreal forest in synchrony with significant climatic warming over the past century. The average rate of migration is estimated to be 12.0 ± 1.0 km decade-1 , being slightly slower than the movement of isotherms (14.7 ± 6.4 km decade-1 ). The migration rate of Mongolian oak is the highest observed among migratory temperate trees (average rate 4.0 ± 1.0 km decade-1 ) and significantly higher than the rates of tree migration at boreal-tundra ecotones (0.9 ± 0.4 km decade-1 ) and alpine treelines (0.004 ± 0.003 km decade-1 ). Compared with the coexisting dominant boreal tree species, Dahurian larch (Larix gmelinii), temperate Mongolian oak is observed to have significantly lower capacity for light acquisition, comparable water-use efficiency but stronger capacity to utilize nutrients especially the most limiting nutrient, nitrogen. In the context of climatic warming, and in addition to a high seed dispersal capacity and potential thermal niche differences, the advantage of nutrient utilization, reflected by foliar elementomes and stable nitrogen isotope ratios, is also likely a key mechanism for Mongolian oak to coexist with Dahurian larch and facilitate its migration toward boreal forest. These findings highlight a rapid deborealization of southern Asian boreal forest in response to climatic warming.

Attributing Atmospheric Phosphorus in the Himalayas: Biomass Burning vs Mineral Dust.

Atmospheric phosphorus is a vital nutrient for ecosystems whose sources and fate are still debated in the fragile Himalayan region, hindering our comprehension of its local ecological impact. This study provides novel insights into atmospheric phosphorus based on the study of total suspended particulate matter at the Qomolangma station. Contrary to the prevailing assumptions, we show that biomass burning (BB), not mineral dust, dominates total dissolved phosphorus (TDP, bioavailable) deposition in this arid region, especially during spring. While total phosphorus is mainly derived from dust (77% annually), TDP is largely affected by the transport of regional biomass-burning plumes from South Asia. During BB pollution episodes, TDP causing springtime TDP fluxes alone accounts for 43% of the annual budget. This suggests that BB outweighs dust in supplying bioavailable phosphorus, a critical nutrient, required to sustain Himalayas' ecological functions. Overall, this first-hand field evidence refines the regional and global phosphorus budget by demonstrating that BB emission, while still unrecognized, is a significant source of P, even in the remote mountains of the Himalayas. It also reveals the heterogeneity of atmospheric phosphorus deposition in that region, which will help predict changes in the impacted ecosystems as the deposition patterns vary.

A distinctive latitudinal trend of nitrogen isotope signature across urban forests in eastern China.

Rapid urbanization has greatly altered nitrogen (N) cycling from regional to global scales. Compared to natural forests, urban forests receive much more external N inputs with distinctive abundances of stable N isotope (δ15 N). However, the large-scale pattern of soil δ15 N and its imprint on plant δ15 N remain less well understood in urban forests. By collecting topsoil (0-20 cm) and leaf samples from urban forest patches in nine large cities across a north-south transect in eastern China, we analyzed the latitudinal trends of topsoil C:N ratio and δ15 N as well as the correlations between tree leaf δ15 N and topsoil δ15 N. We further explored the spatial variation of topsoil δ15 N explained by corresponding climatic, edaphic, vegetation-associated, and anthropogenic drivers. Our results showed a significant increase of topsoil C:N ratio towards higher latitudes, suggesting lower N availability at higher latitudes. Topsoil δ15 N also increased significantly at higher latitudes, being opposite to the latitudinal trend of soil N availability. The latitudinal trend of topsoil δ15 N was mainly explained by mean annual temperature, mean annual precipitation, and atmospheric deposition of both ammonium and nitrate. Consequently, tree leaf δ15 N showed significant positive correlations with topsoil δ15 N across all sampled plant species and functional types. Our findings reveal a distinctive latitudinal trend of δ15 N in urban forests and highlight an important role of anthropogenic N sources in shaping the large-scale pattern of urban forest 15 N signature.

Reassessment of growth-climate relations indicates the potential for decline across Eurasian boreal larch forests.

Larch, a widely distributed tree in boreal Eurasia, is experiencing rapid warming across much of its distribution. A comprehensive assessment of growth on warming is needed to comprehend the potential impact of climate change. Most studies, relying on rigid calendar-based temperature series, have detected monotonic responses at the margins of boreal Eurasia, but not across the region. Here, we developed a method for constructing temporally flexible and physiologically relevant temperature series to reassess growth-temperature relations of larch across boreal Eurasia. Our method appears more effective in assessing the impact of warming on growth than previous methods. Our approach indicates widespread and spatially heterogeneous growth-temperature responses that are driven by local climate. Models quantifying these results project that the negative responses of growth to temperature will spread northward and upward throughout this century. If true, the risks of warming to boreal Eurasia could be more widespread than conveyed from previous works.

Distinct urban-rural gradients of air NO2 and SO2 concentrations in response to emission reductions during 2015-2022 in Beijing, China.

Nitrogen dioxide (NO2) and sulfur dioxide (SO2) are two major air pollutants in urban environment. Emission reduction policies have thus been implemented to improve urban air quality, especially in the metropolises. However, it remains unclear whether the air concentrations of NO2 and SO2 in and around large cities follow a same spatial pattern and how their characteristics change over time in response to the emission reductions. Using ground-based monitoring datasets of air NO2 and SO2 concentrations in Beijing, China, we tested the hypothesis of urban air pollutant islands and evaluated their seasonal and inter-annual variations during 2015-2022. The results showed that air NO2 concentrations increased significantly towards the urban core, being in line with the hypothesis of urban air pollutant island, while air SO2 concentrations showed no such spatial patterns. The urban air NO2 island varied seasonally, with larger radius and higher air NO2 concentrations in spring and winter. In response to the emission reduction, the annual mean radius of the urban air NO2 island showed a rapid decrease from 45.8 km to zero km during the study period. The annual mean air NO2 concentration at the urban core showed a linear decrease at a rate of 4.5 µg m-3 yr-1. In contrast, air SO2 concentration decreased nonlinearly over time and showed a legacy in comparison to the emission reduction. Our findings suggest different urban-rural gradients of air NO2 and SO2 concentrations and highlight their distinct responses to the regional reductions of anthropogenic emissions.

Urban CO2 imprints on carbon isotope and growth of Chinese pine in the Beijing metropolitan region.

Rapid urbanization has occurred globally and resulted in increasing CO2 emissions from urban areas. Compared to natural forests, urban forests are subject to higher atmospheric CO2 concentrations in view of strong urban-periurban-rural gradients of CO2 emissions. However, relevant insights in the CO2-associated urban imprints on the physiology and growth of regional forests remain lacking. By sampling foliage and tree rings of Chinese pine (Pinus tabuliformis) in the Beijing metropolitan region, China, we explored whether and how urban CO2 emissions affect stable carbon isotope ratios (δ13C) and tree growth spatially and/or temporally. The results indicate a significant decrease in foliar δ13C values towards the urban center and this pattern was mainly explained by the urban-periurban-rural gradients of CO2 emissions as surrogated by trunk road density. Tree-ring δ13C values showed a significant decrease over last four decades and this trend was mainly explained by rising levels of CO2 and secondarily mediated by the variations of aridity index during growing season. Moreover, annual basal area increment of Chinese pine was significantly accelerated during last two decades, being mainly driven by increasing CO2 emissions and secondarily mediated by climate variations. These findings reveal significant CO2-associated imprints of urbanization on plant growth and provide empirical evidences of significant CO2-induced alteration of carbon cycles in urban forests.

Hierarchical change-point regression models including random effects to estimate empirical critical loads for nitrogen using Bayesian Regression Models (brms) and JAGS.

The concept of critical loads is used in the framework of the Convention on Long-range Transboundary Air Pollution (UNECE) to define thresholds below which no damaging effects on habitats occur based on the latest scientific knowledge. Change-point regression models applied in a Bayesian framework are useful statistical tools to estimate critical empirical loads. While hierarchical study designs are common in ecological research, previous methods to estimate critical loads using change-point regression did not allow to analyse data collected under such a design. This method update provides an implementation of hierarchical data structure by including random effects such as study sites or as in this example tree species within the Bayesian approach of change-point regression models using two different approaches. The example data set is an European wide gradient study of the impact of climate change and air pollution on forest tree health assessed by foliar nutrient status of nitrogen (N) to phosphorus (P) from 10 different conifer tree species originated from 88 forest sites and 9 countries covering 22 years (1995-2017). Both modelling approaches using JAGS and Bayesian Regression Models using 'Stan' (brms) resulted in reasonable and similar estimations of the critical empirical load for nitrogen (CLempN) for temperate forests. These methodological examples of using different approaches of Bayesian change-point regression models dealing with random effects could prove useful to infer CLempN for other ecosystems and long-term data sets.•Hierarchical change-point regression models are suitable for estimating critical empirical loads.•The Bayesian framework of these models provides the inclusion of the current critical load and various confounding or modifying variables.•Here we present two ways of implementing hierarchical data sets in Bayesian change-point regression models using JAGS and brms.

Strategic roadmap to assess forest vulnerability under air pollution and climate change.

Although it is an integral part of global change, most of the research addressing the effects of climate change on forests have overlooked the role of environmental pollution. Similarly, most studies investigating the effects of air pollutants on forests have generally neglected the impacts of climate change. We review the current knowledge on combined air pollution and climate change effects on global forest ecosystems and identify several key research priorities as a roadmap for the future. Specifically, we recommend (1) the establishment of much denser array of monitoring sites, particularly in the South Hemisphere; (2) further integration of ground and satellite monitoring; (3) generation of flux-based standards and critical levels taking into account the sensitivity of dominant forest tree species; (4) long-term monitoring of N, S, P cycles and base cations deposition together at global scale; (5) intensification of experimental studies, addressing the combined effects of different abiotic factors on forests by assuring a better representation of taxonomic and functional diversity across the ~73,000 tree species on Earth; (6) more experimental focus on phenomics and genomics; (7) improved knowledge on key processes regulating the dynamics of radionuclides in forest systems; and (8) development of models integrating air pollution and climate change data from long-term monitoring programs.

High-level nitrogen additions accelerate soil respiration reduction over time in a boreal forest.

Increased nitrogen (N) inputs are widely recognised to reduce soil respiration (Rs), but how N deposition affects the temporal dynamics of Rs remains unclear. Using a decade-long fertilisation experiment in a boreal larch forest (Larix gmelini) in northeast China, we found that the effects of N additions on Rs showed a temporal shift from a positive effect in the short-term (increased by 8% on average in the first year) to a negative effect over the longer term (decreased by 21% on average in the 11th year). The rates of decrease in Rs for the higher N levels were almost twice as high as those of the low N level. Our results suggest that the reduction in Rs in response to increased N input is accelerated by high-level N additions, and experimental high N applications are likely to overestimate the contribution of N deposition to soil carbon sequestration in a boreal forest.

Plant sizes and shapes above and belowground and their interactions with climate.

Although the above and belowground sizes and shapes of plants strongly influence plant competition, community structure, and plant-environment interactions, plant sizes and shapes remain poorly characterized across climate regimes. We investigated relationships among shoot and root system size and climate. We assembled and analyzed, to our knowledge, the largest global database describing the maximum rooting depth, lateral spread, and shoot size of terrestrial plants - more than doubling the Root Systems of Individual Plants database to 5647 observations. Water availability and growth form greatly influence shoot size, and rooting depth is primarily influenced by temperature seasonality. Shoot size is the strongest predictor of lateral spread, with root system diameter being two times wider than shoot width on average for woody plants. Shoot size covaries strongly with rooting system size; however, the geometries of plants differ considerably across climates, with woody plants in more arid climates having shorter shoots, but deeper, narrower root systems. Additionally, estimates of the depth and lateral spread of plant root systems are likely underestimated at the global scale.

Retention of deposited ammonium and nitrate and its impact on the global forest carbon sink.

The impacts of enhanced nitrogen (N) deposition on the global forest carbon (C) sink and other ecosystem services may depend on whether N is deposited in reduced (mainly as ammonium) or oxidized forms (mainly as nitrate) and the subsequent fate of each. However, the fates of the two key reactive N forms and their contributions to forest C sinks are unclear. Here, we analyze results from 13 ecosystem-scale paired 15N-labelling experiments in temperate, subtropical, and tropical forests. Results show that total ecosystem N retention is similar for ammonium and nitrate, but plants take up more labelled nitrate ([Formula: see text]%) ([Formula: see text]) than ammonium ([Formula: see text]%) while soils retain more ammonium ([Formula: see text]%) than nitrate ([Formula: see text]%). We estimate that the N deposition-induced C sink in forests in the 2010s  is [Formula: see text] Pg C yr-1, higher than previous estimates because of a larger role for oxidized N and greater rates of global N deposition.

Nonlinear responses of ecosystem carbon fluxes to nitrogen deposition in an old-growth boreal forest.

Nitrogen (N) deposition is known to increase carbon (C) sequestration in N-limited boreal forests. However, the long-term effects of N deposition on ecosystem carbon fluxes have been rarely investigated in old-growth boreal forests. Here we show that decade-long experimental N additions significantly stimulated net primary production (NPP) but the effect decreased with increasing N loads. The effect on soil heterotrophic respiration (Rh) shifted from a stimulation at low-level N additions to an inhibition at higher levels of N additions. Consequently, low-level N additions resulted in a neutral effect on net ecosystem productivity (NEP), due to a comparable stimulating effect on NPP and Rh, while NEP was increased by high-level N additions. Moreover, we found nonlinear temporal responses of NPP, Rh and NEP to low-level N additions. Our findings imply that actual N deposition in boreal forests likely exerts a minor contribution to their soil C storage.

The diurnal cycle of summer tropospheric ozone concentrations across Chinese cities: Spatial patterns and main drivers.

China is experiencing severe tropospheric ozone pollution, especially during the summer period in cities. Previous studies have assessed the role of meteorological conditions and anthropogenic precursors in shaping the diurnal variation of ozone concentration in some Chinese cities or the spatial patterns of daytime ozone concentration, but less is known about the spatial variation and main regulators of the diurnal cycle of summer ozone concentrations in Chinese cities. Using monitoring data from 367 cities, we analyzed the spatial patterns and main regulators of daytime maximum, nighttime minimum and diurnal difference of summer (June-August) ozone concentration during 2015-2019. National mean values and standard deviations of daytime maximum and nighttime minimum of summer surface ozone concentration were 124.1 ± 27.5 and 33.4 ± 13.0 µg m-3, resulting in a diurnal difference of 90.7 ± 25.2 µg m-3. High values of daytime maximum, nighttime minimum, and diurnal difference of summer ozone concentration occurred in cities in northern China, especially in the North China Plain, and several city agglomerations in southern China. Daytime maximum ozone concentration was higher in cities with higher daytime PM2.5 and NO2 concentrations, lower daytime precipitation and lower elevation. Nighttime minimum ozone concentration increased with lower nighttime precipitation, lower NO2 concentration and CO concentration, higher nighttime maximum PM2.5 concentration and higher elevation. Diurnal difference of ozone concentration increased with lower elevation, lower daytime precipitation, and higher diurnal difference of CO and NO2 concentrations. Our findings highlight different regulators for daytime and nighttime ozone and imply the need of joint regulation of PM2.5 and NO2 emissions to control ozone pollution.

Spatially divergent trends of nitrogen versus phosphorus limitation across European forests.

Nitrogen (N) and phosphorus (P) are essential nutrients that widely limit plant growth in global terrestrial ecosystems. Rising atmospheric CO2 concentration generally stimulates terrestrial net primary productivity and consequently may cause or aggravate N and P limitation due to a dilution effect, but the spatial variation of temporal trends in N versus P limitation and its key regulating factors is poorly understood. Using the leaf N:P ratio of 15 dominant tree species as an indicator, we analysed the spatial variation of plot-level shift towards N or P limitation across 163 European forest plots during 1995-2017. Phosphorus limitation increased from 25% to 33% of the studied plots between 1995-1997 and 2015-2017, while N limitation occurred in a negligible number of plots. A major proportion (56%) of the plots showed no significant trend in leaf N:P ratio, implying no shifts in N versus P limitation status. In the remaining plots, 38% of the plots showed a significant increase of leaf N:P ratio and only 6% of the plots showed a significant decrease of leaf N:P ratio. The spatial variation in the rate of decrease in leaf N:P ratio was associated with a significant decrease in leaf N concentration and mainly explained by the rate of decrease in N deposition. In contrast, the spatial variation in the rate of increase in leaf N:P ratio was associated with a significant decrease in leaf P concentration and mainly explained by forest category (broadleaf vs. conifer), mean annual temperature and soil C:N ratio. Our findings highlight a remarkable spatial divergence in temporal trends of nutrient limitation status across European forests over the past two decades, but overall, P is becoming more limiting versus N, especially in broadleaved forests.

Effects of nitrogen addition on microbial residues and their contribution to soil organic carbon in China's forests from tropical to boreal zone.

Atmospheric nitrogen (N) deposition has a significant influence on soil organic carbon (SOC) accumulation in forest ecosystems. Microbial residues, as by-products of microbial anabolism, account for a significant fraction of soil C pools. However, how N deposition affects the accumulation of soil microbial residues in different forest biomes remains unclear. Here, we investigated the effects of six/seven-year N additions on microbial residues (amino sugar biomarkers) in eight forests from tropical to boreal zone in eastern China. Our results showed a minor change in the soil microbial residue concentrations but a significant change in the contribution of microbial residue-C to SOC after N addition. The contribution of fungal residue-C to SOC decreased under low N addition (50 kg N ha-1 yr-1) in the tropical secondary forest (-19%), but increased under high N addition (100 kg N ha-1 yr-1) in the temperate Korean pine mixed forest (+21%). The contribution of bacterial residue-C to SOC increased under the high N addition in the subtropical Castanopsis carlesii forest (+26%) and under the low N addition in the temperate birch forest (+38%), respectively. The responses of microbial residue-C in SOC to N addition depended on the changes in soil total N concentration and fungi to bacteria ratio under N addition and climate. Taken together, these findings provide the experimental evidence that N addition diversely regulates the formation and composition of microbial-derived C in SOC in forest ecosystems.

Effects of nitrogen addition on soil methane uptake in global forest biomes.

Nitrogen (N) deposition has been conventionally thought to decrease forest soil methane (CH4) uptake, while the biome specific and dose dependent effect is poorly understood. Based on a meta-analysis of 63 N addition trials from 7 boreal forests, 8 temperate forests, 13 subtropical and 4 tropical forests, we evaluated the effects of N addition on soil CH4 uptake fluxes across global forest biomes. When combining all N addition levels, soil CH4 uptake was insignificantly decreased by 7% in boreal forests, while N addition significantly decreased soil CH4 uptake by 39% in temperate forests and by 21% in subtropical and tropical forests, respectively. Meta-regression analyses, however, indicated a shift from a positive to a negative effect on soil CH4 uptake with increasing N additions both in boreal forests (threshold = 48 kg N ha-1 yr-1) and temperate forests (threshold = 27 kg N ha-1 yr-1), while no such shift was found in subtropical and tropical forests. Considering that current N deposition to most boreal and temperate forests is below the abovementioned thresholds, N deposition likely exerts a positive to neutral effect on soil CH4 uptake in both forest biomes. Our results provide new insights on the biome specific and dose dependent effect of N addition on soil CH4 sink in global forests and suggest that the current understanding that N deposition decreases forest soil CH4 uptake is flawed by high levels of experimental N addition.

Family-level leaf nitrogen and phosphorus stoichiometry of global terrestrial plants.

Leaf nitrogen (N) and phosphorus (P) concentrations are critical for photosynthesis, growth, reproduction and other ecological processes of plants. Previous studies on large-scale biogeographic patterns of leaf N and P stoichiometric relationships were mostly conducted using data pooled across taxa, while family/genus-level analyses are rarely reported. Here, we examined global patterns of family-specific leaf N and P stoichiometry using a global data set of 12,716 paired leaf N and P records which includes 204 families, 1,305 genera, and 3,420 species. After determining the minimum size of samples (i.e., 35 records), we analyzed leaf N and P concentrations, N:P ratios and N∼P scaling relationships of plants for 62 families with 11,440 records. The numeric values of leaf N and P stoichiometry varied significantly across families and showed diverse trends along gradients of mean annual temperature (MAT) and mean annual precipitation (MAP). The leaf N and P concentrations and N:P ratios of 62 families ranged from 6.11 to 30.30 mg g-1, 0.27 to 2.17 mg g-1, and 10.20 to 35.40, respectively. Approximately 1/3-1/2 of the families (22-35 of 62) showed a decrease in leaf N and P concentrations and N:P ratios with increasing MAT or MAP, while the remainder either did not show a significant trend or presented the opposite pattern. Family-specific leaf N∼P scaling exponents did not converge to a certain empirical value, with a range of 0.307-0.991 for 54 out of 62 families which indicated a significant N∼P scaling relationship. Our results for the first time revealed large variation in the family-level leaf N and P stoichiometry of global terrestrial plants and that the stoichiometric relationships for at least one-third of the families were not consistent with the global trends reported previously. The numeric values of the family-specific leaf N and P stoichiometry documented in the current study provide critical synthetic parameters for biogeographic modeling and for further studies on the physiological and ecological mechanisms underlying the nutrient use strategies of plants from different phylogenetic taxa.