Drought shortens subtropical understory growing season by advancing leaf senescence.

Subtropical forests, recognized for their intricate vertical canopy stratification, exhibit high resistance to extreme drought. However, the response of leaf phenology to drought in the species-rich understory remains poorly understood. In this study, we constructed a digital camera system, amassing over 360,000 images through a 70% throughfall exclusion experiment, to explore the drought response of understory leaf phenology. The results revealed a significant advancement in understory leaf senescence phenology under drought, with 11.75 and 15.76 days for the start and end of the leaf-falling event, respectively. Pre-season temperature primarily regulated leaf development phenology, whereas soil water dominated the variability in leaf senescence phenology. Under drought conditions, temperature sensitivities for the end of leaf emergence decreased from -13.72 to -11.06 days °C-1, with insignificance observed for the start of leaf emergence. Consequently, drought treatment shortened both the length of the growing season (15.69 days) and the peak growth season (9.80 days) for understory plants. Moreover, this study identified diverse responses among intraspecies and interspecies to drought, particularly during the leaf development phase. These findings underscore the pivotal role of water availability in shaping understory phenology patterns, especially in subtropical forests.

Enhanced causal effect of ecosystem photosynthesis on respiration during heatwaves.

Because of global warming, Earth's ecosystems have been experiencing more frequent and severe heatwaves. Heatwaves are expected to tip terrestrial carbon sequestration by elevating ecosystem respiration and suppressing gross primary productivity (GPP). Here, using the convergent cross-mapping technique, this study detected positive bidirectional causal effects between GPP and respiration in two unprecedented European heatwaves. Heatwaves enhanced the causal effect strength of GPP on respiration rather than respiration on GPP across 40 site-years of observations. Further analyses and global simulations revealed spatial heterogeneity in the heatwave response of the causal link strength between GPP and respiration, which was jointly driven by the local climate and vegetation properties. However, the causal effect strength of GPP on respiration showed considerable uncertainties in CMIP6 models. This study reveals an enhanced causal link strength between GPP and respiration during heatwaves, shedding light on improving projections for terrestrial carbon sink dynamics under future climate extremes.

Evidence for widespread thermal optimality of ecosystem respiration.

Ecosystem respiration (ER) is among the largest carbon fluxes between the biosphere and the atmosphere. Understanding the temperature response of ER is crucial for predicting the climate change-carbon cycle feedback. However, whether there is an apparent optimum temperature of ER ([Formula: see text]) and how it changes with temperature remain poorly understood. Here we analyse the temperature response curves of ER at 212 sites from global FLUXNET. We find that ER at 183 sites shows parabolic temperature response curves and [Formula: see text] at which ER reaches the maximum exists widely across biomes around the globe. Among the 15 biotic and abiotic variables examined, [Formula: see text] is mostly related to the optimum temperature of gross primary production (GPP, [Formula: see text]) and annual maximum daily temperature (Tmax). In addition, [Formula: see text] linearly increases with Tmax across sites and over vegetation types, suggesting its thermal adaptation. The adaptation magnitude of [Formula: see text], which is measured by the change in [Formula: see text] per unit change in Tmax, is positively correlated with the adaptation magnitude of [Formula: see text]. This study provides evidence of the widespread existence of [Formula: see text] and its thermal adaptation with Tmax across different biomes around the globe. Our findings suggest that carbon cycle models that consider the existence of [Formula: see text] and its adaptation have the potential to more realistically predict terrestrial carbon sequestration in a world with changing climate.

Nitrogen use strategy drives interspecific differences in plant photosynthetic CO2 acclimation.

Rising atmospheric CO2 concentration triggers an emergent phenomenon called plant photosynthetic acclimation to elevated CO2 (PAC). PAC is often characterized by a reduction in leaf photosynthetic capacity (Asat ), which varies dramatically along the continuum of plant phylogeny. However, it remains unclear whether the mechanisms responsible for PAC are also different across plant phylogeny, especially between gymnosperms and angiosperms. Here, by compiling a dataset of 73 species, we found that although leaf Asat increased significantly from gymnosperms to angiosperms, there was no phylogenetic signal in the PAC magnitude along the phylogenetic continuum. Physio-morphologically, leaf nitrogen concentration (Nm ), photosynthetic nitrogen-use efficiency (PNUE), and leaf mass per area (LMA) dominated PAC for 36, 29, and 8 species, respectively. However, there was no apparent difference in PAC mechanisms across major evolutionary clades, with 75% of gymnosperms and 92% of angiosperms regulated by the combination of Nm and PNUE. There was a trade-off between Nm and PNUE in driving PAC across species, and PNUE dominated the long-term changes and inter-specific differences in Asat under elevated CO2 . These findings indicate that nitrogen-use strategy drives the acclimation of leaf photosynthetic capacity to elevated CO2 across terrestrial plant species.

Across-model spread and shrinking in predicting peatland carbon dynamics under global change.

Large across-model spread in simulating land carbon (C) dynamics has been ubiquitously demonstrated in model intercomparison projects (MIPs), and became a major impediment in advancing climate change prediction. Thus, it is imperative to identify underlying sources of the spread. Here, we used a novel matrix approach to analytically pin down the sources of across-model spread in transient peatland C dynamics in response to a factorial combination of two atmospheric CO2 levels and five temperature levels. We developed a matrix-based MIP by converting the C cycle module of eight land models (i.e., TEM, CENTURY4, DALEC2, TECO, FBDC, CASA, CLM4.5 and ORCHIDEE) into eight matrix models. While the model average of ecosystem C storage was comparable to the measurement, the simulation differed largely among models, mainly due to inter-model difference in baseline C residence time. Models generally overestimated net ecosystem production (NEP), with a large spread that was mainly attributed to inter-model difference in environmental scalar. Based on the sources of spreads identified, we sequentially standardized model parameters to shrink simulated ecosystem C storage and NEP to almost none. Models generally captured the observed negative response of NEP to warming, but differed largely in the magnitude of response, due to differences in baseline C residence time and temperature sensitivity of decomposition. While there was a lack of response of NEP to elevated CO2 (eCO2 ) concentrations in the measurements, simulated NEP responded positively to eCO2 concentrations in most models, due to the positive responses of simulated net primary production. Our study used one case study in Minnesota peatland to demonstrate that the sources of across-model spreads in simulating transient C dynamics can be precisely traced to model structures and parameters, regardless of their complexity, given the protocol that all the matrix models were driven by the same gross primary production and environmental variables.

Nitrogen enrichment alters multiple dimensions of grassland functional stability via changing compositional stability.

Anthropogenic nutrient enrichment is known to alter the composition and functioning of plant communities. However, how nutrient enrichment influences multiple dimensions of community- and ecosystem-level stability remains poorly understood. Using data from a nitrogen (N) and phosphorus (P) addition experiment in a temperate semi-arid grassland that experienced a natural drought, we show that N enrichment, not P enrichment, decreased grassland functional and compositional temporal stability, resistance and recovery but increased functional and compositional resilience. Compositional stability and species asynchrony, rather than species diversity, were identified as key determinants of all dimensions of grassland functional stability, except for recovery. Whereas grassland functional recovery was decoupled from compositional recovery, N enrichment altered other dimensions of functional stability primarily through changing their corresponding compositional stability dimensions. Our findings highlight the need to examine ecological stability at the community level for a more mechanistic understanding of ecosystem dynamics in the face of environmental change.

Effects of Vegetation Phenology on Ecosystem Water Use Efficiency in a Semiarid Region of Northern China.

Water use efficiency (WUE) is an important ecosystem functional property for measuring coupled relationships of the carbon-water cycle. Both biotic and environmental factors significantly impact WUE in terrestrial ecosystems. Relationships between environmental factors and WUE have been well discussed in previous studies. Although the crucial role of vegetation phenology, a common indicator of climate-vegetation interactions, in regulating the WUE has been widely reported, the underlying mechanism has not yet to be fully elucidated. Here, we utilized multiple long-term remote sensing datasets to analyze the interannual changes in seasonal WUE, and discussed the potential associations between phenology and WUE in the Luanhe River basin, which is a typical semiarid region of China, from 1988 to 2015. Most of the pixels across Luanhe River basin showed increasing spring WUE and decreasing autumn WUE. The start of the growing season (SOS) was slightly advanced by an average of 0.06 days per year (d/y) in the whole study area, with a delayed trend (0.04 d/y) in the upper reaches region (UR) and an advanced trend (0.20 d/y) in the middle-lower reaches region (MLR). The end of the growing season (EOS) showed a generally delayed trend (0.11 d/y) across the basin. Furthermore, negative correlations were detected between spring WUE and SOS in the UR, while positive correlations were detected in the MLR. The opposite patterns of the relationship of autumn WUE and EOS were also found between UR and MLR. The results were attributed to the balance and compensation of biotic and abiotic factors in the regulation of ecosystem WUE. Our findings provide new insights into the interaction between vegetation dynamics and carbon-water cycle coupling.

Plant photosynthetic overcompensation under nocturnal warming: lack of evidence in subtropical evergreen trees.

BACKGROUND AND AIMS: Increased plant photosynthesis under nocturnal warming is a negative feedback mechanism to overcompensate for night-time carbon loss to mitigate climate warming. This photosynthetic overcompensation effect has been observed in dry deciduous ecosystems but whether it exists in subtropical wet forest trees is unclear. METHODS: Two subtropical evergreen tree species (Schima superba and Castanopsis sclerophylla) were grown in a greenhouse and exposed to ambient and elevated night-time temperature. The occurrence of the photosynthetic overcompensation effect was determined by measuring daytime and night-time leaf gas exchange and non-structural carbohydrate (NSC) concentration. KEY RESULTS: A reduction in leaf photosynthesis for both species and an absence of persistent photosynthetic overcompensation were observed. The photosynthetic overcompensation effect was transient in S. superba due to respiratory acclimation and stomatal limitation. For S. superba, nocturnal warming resulted in insufficient changes in night-time respiration and NSC concentration to stimulate overcompensation and inhibited leaf stomatal conductance by increasing the leaf-to-air vapour pressure deficit. CONCLUSIONS: The results indicate that leaf stomatal conductance is important for the photosynthetic overcompensation effect in different tree species. The photosynthetic overcompensation effect under nocturnal warming may be a transient occurrence rather than a persistent mechanism in subtropical forest ecosystems.

Can evolutionary history predict plant plastic responses to climate change?

Plant plastic responses are critical to the adaptation and survival of species under climate change, but whether they are constrained by evolutionary history (phylogeny) is largely unclear. Plant leaf traits are key in determining plants' performance in different environments, and if these traits and their variation are phylogenetically dependent, predictions could be made to identify species vulnerable to climate change. We compiled data on three leaf traits (photosynthetic rate, specific leaf area, and leaf nitrogen content) and their variation under four environmental change scenarios (warming, drought, elevated CO2 , or nitrogen addition) for 434 species, from 210 manipulation experiments. We found phylogenetic signal in the three traits but not in their variation under the four scenarios. This indicates that closely related species show similar traits but that their plastic responses could not be predicted from species relatedness under environmental change. Meanwhile, phylogeny weakened the slopes but did not change the directions of conventional pairwise trait relationships, suggesting that co-evolved leaf trait pairs have consistent responses under contrasting environmental conditions. Phylogeny can identify lineages rich in species showing similar traits and predict their relationships under climate change, but the degree of plant phenotypic variation does not vary consistently across evolutionary clades.

Soil phosphorus drives plant trait variations in a mature subtropical forest.

Earth system models are implementing soil phosphorus dynamic and plant functional traits to predict functional changes in global forests. However, the linkage between soil phosphorus and plant traits lacks empirical evidence, especially in mature forests. Here, we examined the soil phosphorus constraint on plant functional traits in a mature subtropical forest based on observations of 9943 individuals from 90 species in a 5-ha forest dynamic plot and 405 individuals from 15 species in an adjacent 10-year nutrient-addition experiment. We first confirmed a pervasive phosphorus limitation on subtropical tree growth based on leaf N:P ratios. Then, we found that soil phosphorus dominated multidimensional trait variations in the 5-ha forest dynamic plot. Soil phosphorus content explained 44% and 53% of the variance in the traits defining the main functional space across species and communities, respectively. Lastly, we found much stronger phosphorus effects on most plant functional traits than nitrogen at both species and community levels in the 10-year nutrient-addition experiment. This study provides evidence for the consistent pattern of soil phosphorus constraint on plant trait variations between the species and community levels in a mature evergreen broadleaf forest in the East Asian monsoon region. These findings shed light on the predominant role of soil phosphorus on plant functional trait variations in mature subtropical forests, providing new insights for models to incorporate soil phosphorus constraint in predicting future vegetation dynamics.

Consistent temperature-dependent patterns of leaf lifespan across spatial and temporal gradients for deciduous trees in Europe.

Temperature affects leaf lifespan (LL) across either space or time, driving long-term adaptation and short-term thermal acclimation, respectively. However, a comprehensive understanding of the phenomenon and the underlying phenological mechanisms remain poorly understood. The present study investigated the relationship between LL and temperature in six common deciduous trees across both spatial and temporal gradients, then explained the LL variation patterns based on phenological shifts. Using long-term (1971-2000) phenological records of six deciduous tree species at 54 sites across central Europe, we analyzed spatial and temporal variations of LL and leaf phenology along temperature gradients. We assessed the relative contribution of phenological shifts to LL variations by comparing absolute changes in leaf-out and leaf fall. We reported positive LL-temperature relationships across all observations along both spatial (+3.32 days/°C) and temporal (+4.43 days/°C) gradients. The paired t-test of the six deciduous tree species showed no significant difference in regression slopes of LL- temperature between the two gradients (t = -1.50, df = 5, P = 0.194). Prolonged LL can be explained mainly by earlier leaf-out induced by warmer temperatures both spatially (-3.22 days/°C) and temporally (-4.08 days/°C). The converging temperature-dependent patterns of LL across time and space indicate that short-term thermal acclimation keeps pace with long-term genetic adaptation for deciduous trees in Europe. Earlier leaf-out is the key force shaping the LL-temperature relationship. These results provide insights for predicting future vegetation dynamics under global warming.

Experimental warming reduces ecosystem resistance and resilience to severe flooding in a wetland.

Climate warming and extreme hydrological events are threatening the sustainability of wetlands across the globe. However, whether climate warming will amplify or diminish the impact of extreme flooding on wetland ecosystems is unknown. Here, we show that climate warming significantly reduced wetland resistance and resilience to a severe flooding event via a 6-year warming experiment. We first found that warming rapidly altered plant community structure by increasing the dominance of low-canopy species. Then, we showed that warming reduced the resistance and resilience of vegetation productivity to a 72-cm flooding event. Last, we detected slower postflooding carbon processes, such as gross ecosystem productivity, soil respiration, and soil methane emission, under the warming treatment. Our results demonstrate how severe flooding can destabilize wetland vegetation structure and ecosystem function under climate warming. These findings indicate an enhanced footprint of extreme hydrological events in wetland ecosystems in a warmer climate.

Nitrogen addition promotes soil microbial beta diversity and the stochastic assembly.

Nitrogen (N) deposition is one of major environmental concerns and alters the microbial communities in the pedosphere. A central debate in governing microbial community is on the relative importance of deterministic (ecological selection) vs. stochastic processes (dispersal, drift, diversification or speciation), which consequently limited our understanding of microbial assembly in response to N addition. Here, we conducted a global analysis of high-throughput sequencing data to reveal the mechanisms of N-addition effects on soil microbial communities. The results show that N addition significantly shifted the microbial community structure and promoted microbial beta diversity, particularly in the N-limited ecosystems. Changes in microbial structure and beta diversity increased significantly as the N addition rate, study duration, and the degree of soil acidification increased. The stochastic processes are more important than the deterministic processes for microbial community assembly, while N addition significantly increase the importance of stochastic processes whether the phylogenetic relationship is considered or not. Overall, the current study highlights the important of ecological stochasticity in regulating microbial assembly under N deposition scenarios.

Nocturnal warming accelerates drought-induced seedling mortality of two evergreen tree species.

Extreme drought is one of the key climatic drivers of tree mortality on a global scale. However, it remains unclear whether the drought-induced tree mortality will increase under nocturnal climate warming. Here we exposed seedlings of two wide-ranging subtropical tree species, Castanopsis sclerophylla and Schima superba, with contrasting stomatal regulation strategies to prolonged drought under ambient and elevated night-time temperature by 2 °C. We quantified the seedling survival time since drought treatment by measuring multiple leaf traits such as leaf gas exchange, predawn leaf water potential and water-use efficiency. The results showed that all seedlings in the ambient temperature died within 180 days and 167 days of drought for C. sclerophylla and S. superba, respectively. Night warming significantly shortened the survival time of C. sclerophylla, by 31 days, and S. superba by 28 days, under the drought treatment. A survival analysis further showed that seedlings under night warming suffered a 1.6 times greater mortality risk than those under ambient temperature. Further analyses revealed that night warming suppressed net leaf carbon gain in both species by increasing the nocturnal respiratory rate of S. superba across the first 120 days of drought and decreasing the photosynthetic rate of both species generally after 46 days of drought. These effects on net carbon gain were more pronounced in S. superba than C. sclerophylla. After 60 days of drought, night warming decreased the predawn leaf water potential and leaf water-use efficiency of C. sclerophylla but not S. superba. These contrasting responses are partially due to variations in stomatal control between the two species. These findings suggest that stomatal traits can regulate the response of leaf gas exchange and plant water-use to nocturnal warming during drought. This study indicates that nocturnal warming can accelerate tree mortality during drought. Night warming accelerates the mortality of two subtropical seedlings under drought.Night warming differently affects the drought response of leaf gas exchange and plant water-use between the two species due to species-specific stomatal morphological traits.Carbon metabolism changes and hydraulic damage play differential roles in driving night-warming impacts on the drought-induced mortality between the two species.

Warming reshaped the microbial hierarchical interactions.

Global warming may alter microbially mediated ecosystem functions through reshaping of microbial diversity and modified microbial interactions. Here, we examined the effects of 5-year experimental warming on different microbial hierarchical groups in a coastal nontidal soil ecosystem, including prokaryotes (i.e., bacteria and archaea), fungi, and Cercozoa, which is a widespread phylum of protists. Warming significantly altered the diversity and structure of prokaryotic and fungal communities in soil and additionally decreased the complexity of the prokaryotic network and fragmented the cercozoan network. By using the Inter-Domain Ecological Network approach, the cross-trophic interactions among prokaryotes, fungi, and Cercozoa were further investigated. Under warming, cercozoan-prokaryotic and fungal-prokaryotic bipartite networks were simplified, whereas the cercozoan-fungal network became slightly more complex. Despite simplification of the fungal-prokaryotic network, the strengthened synergistic interactions between saprotrophic fungi and certain prokaryotic groups, such as the Bacteroidetes, retained these phyla within the network under warming. In addition, the interactions within the fungal community were quite stable under warming conditions, which stabilized the interactions between fungi and prokaryotes or protists. Additionally, we found the microbial hierarchical interactions were affected by environmental stress (i.e., salinity and pH) and soil nutrients. Interestingly, the relevant microbial groups could respond to different soil properties under ambient conditions, whereas under warming these two groups tended to respond to similar soil properties, suggesting network hub species responded to certain environmental changes related to warming, and then transferred this response to their partners through trophic interactions. Finally, warming strengthened the network modules' negative association with soil organic matters through some fungal hub species, which might trigger soil carbon loss in this ecosystem. Our study provides new insights into the response and feedback of microbial hierarchical interactions under warming scenario.

A small climate-amplifying effect of climate-carbon cycle feedback.

The climate-carbon cycle feedback is one of the most important climate-amplifying feedbacks of the Earth system, and is quantified as a function of carbon-concentration feedback parameter (ß) and carbon-climate feedback parameter (γ). However, the global climate-amplifying effect from this feedback loop (determined by the gain factor, g) has not been quantified from observations. Here we apply a Fourier analysis-based carbon cycle feedback framework to the reconstructed records from 1850 to 2017 and 1000 to 1850 to estimate ß and γ. We show that the ß-feedback varies by less than 10% with an average of 3.22 ± 0.32 GtC ppm-1 for 1880-2017, whereas the γ-feedback increases from -33 ± 14 GtC K-1 on a decadal scale to -122 ± 60 GtC K-1 on a centennial scale for 1000-1850. Feedback analysis further reveals that the current amplification effect from the carbon cycle feedback is small (g is 0.01 ± 0.05), which is much lower than the estimates by the advanced Earth system models (g is 0.09 ± 0.04 for the historical period and is 0.15 ± 0.08 for the RCP8.5 scenario), implying that the future allowable CO2 emissions could be 9 ± 7% more. Therefore, our findings provide new insights about the strength of climate-carbon cycle feedback and about observational constraints on models for projecting future climate.

Reduced magnitude and shifted seasonality of CO2 sink by experimental warming in a coastal wetland.

Coastal wetlands have the highest carbon sequestration rate per unit area among all unmanaged natural ecosystems. However, how the magnitude and seasonality of the CO2 sink in coastal wetlands will respond to future climate warming remains unclear. Here, based on measurements of ecosystem CO2 fluxes in a field experiment in the Yellow River Delta, we found that experimental warming (i.e., a 2.4°C increase in soil temperature) reduced net ecosystem productivity (NEP) by 23.7% across two growing seasons of 2017-2018. Such a reduction in NEP resulted from the greater decrease in gross primary productivity (GPP) than ecosystem respiration (ER) under warming. The negative warming effect on NEP mainly occurred in summer (-43.9%) but not in autumn (+61.3%), leading to a shifted NEP seasonality under warming. Further analyses showed that the warming effects on ecosystem CO2 exchange were mainly controlled by soil salinity and its corresponding impacts on species composition. For example, warming increased soil salinity (+35.0%), reduced total aboveground biomass (-9.9%), and benefited the growth of plant species with high salt tolerance and late peak growth. To the best of our knowledge, this study provides the first experimental evidence on the reduced magnitude and shifted seasonality of CO2 exchange under climate warming in coastal wetlands. These findings underscore the high vulnerability of wetland CO2 sink in coastal regions under future climate change.

Warming-induced global soil carbon loss attenuated by downward carbon movement.

The fate of soil organic carbon (SOC) under warming is poorly understood, particularly across large extents and in the whole-soil profile. Using a data-model integration approach applied across the globe, we find that downward movement of SOC along the soil profile reduces SOC loss under warming. We predict that global SOC stocks (down to 2 m) will decline by 4% (~80 Pg) on average when SOC reaches the steady state under 2°C warming, assuming no changes in net primary productivity (NPP). To compensate such decline (i.e. maintain current SOC stocks), a 3% increase of NPP is required. Without the downward SOC movement, global SOC declines by 15%, while a 20% increase in NPP is needed to compensate that loss. This vital role of downward SOC movement in controlling whole-soil profile SOC dynamics in response to warming is due to the protection afforded to downward-moving SOC by depth, indicated by much longer residence times of SOC in deeper layers. Additionally, we find that this protection could not be counteracted by promoted decomposition due to the priming of downward-moving new SOC from upper layers on native old SOC in deeper layers. This study provides the first estimation of whole-soil SOC changes under warming and additional NPP required to compensate such changes across the globe, and reveals the vital role of downward movement of SOC in reducing SOC loss under global warming.

Research challenges and opportunities for using big data in global change biology.

Global change biology has been entering a big data era due to the vast increase in availability of both environmental and biological data. Big data refers to large data volume, complex data sets, and multiple data sources. The recent use of such big data is improving our understanding of interactions between biological systems and global environmental changes. In this review, we first explore how big data has been analyzed to identify the general patterns of biological responses to global changes at scales from gene to ecosystem. After that, we investigate how observational networks and space-based big data have facilitated the discovery of emergent mechanisms and phenomena on the regional and global scales. Then, we evaluate the predictions of terrestrial biosphere under global changes by big modeling data. Finally, we introduce some methods to extract knowledge from big data, such as meta-analysis, machine learning, traceability analysis, and data assimilation. The big data has opened new research opportunities, especially for developing new data-driven theories for improving biological predictions in Earth system models, tracing global change impacts across different organismic levels, and constructing cyberinfrastructure tools to accelerate the pace of model-data integrations. These efforts will uncork the bottleneck of using big data to understand biological responses and adaptations to future global changes.

Robust leaf trait relationships across species under global environmental changes.

Recent studies show coordinated relationships between plant leaf traits and their capacity to predict ecosystem functions. However, how leaf traits will change within species and whether interspecific trait relationships will shift under future environmental changes both remain unclear. Here, we examine the bivariate correlations between leaf economic traits of 515 species in 210 experiments which mimic climate warming, drought, elevated CO2, and nitrogen deposition. We find divergent directions of changes in trait-pairs between species, and the directions mostly do not follow the interspecific trait relationships. However, the slopes in the logarithmic transformed interspecific trait relationships hold stable under environmental changes, while only their elevations vary. The elevation changes of trait relationship are mainly driven by asymmetrically interspecific responses contrary to the direction of the leaf economic spectrum. These findings suggest robust interspecific trait relationships under global changes, and call for linking within-species responses to interspecific coordination of plant traits.