Mapping the global distribution of C4 vegetation using observations and optimality theory.

Plants with the C4 photosynthesis pathway typically respond to climate change differently from more common C3-type plants, due to their distinct anatomical and biochemical characteristics. These different responses are expected to drive changes in global C4 and C3 vegetation distributions. However, current C4 vegetation distribution models may not predict this response as they do not capture multiple interacting factors and often lack observational constraints. Here, we used global observations of plant photosynthetic pathways, satellite remote sensing, and photosynthetic optimality theory to produce an observation-constrained global map of C4 vegetation. We find that global C4 vegetation coverage decreased from 17.7% to 17.1% of the land surface during 2001 to 2019. This was the net result of a reduction in C4 natural grass cover due to elevated CO2 favoring C3-type photosynthesis, and an increase in C4 crop cover, mainly from corn (maize) expansion. Using an emergent constraint approach, we estimated that C4 vegetation contributed 19.5% of global photosynthetic carbon assimilation, a value within the range of previous estimates (18-23%) but higher than the ensemble mean of dynamic global vegetation models (14 ± 13%; mean ± one standard deviation). Our study sheds insight on the critical and underappreciated role of C4 plants in the contemporary global carbon cycle.

Diminishing carryover benefits of earlier spring vegetation growth.

Spring vegetation growth can benefit summer growth by increasing foliage area and carbon sequestration potential, or impair it by consuming additional resources needed for sustaining subsequent growth. However, the prevalent driving mechanism and its temporal changes remain unknown. Using satellite observations and long-term atmospheric CO2 records, here we show a weakening trend of the linkage between spring and summer vegetation growth/productivity in the Northern Hemisphere during 1982-2021. This weakening is driven by warmer and more extreme hot weather that becomes unfavourable for peak-season growth, shifting peak plant functioning away from earlier periods. This is further exacerbated by seasonally growing ecosystem water stress due to reduced water supply and enhanced water demand. Our finding suggests that beneficial carryover effects of spring growth on summer growth are diminishing or even reversing, acting as an early warning sign of the ongoing shift of climatic effects from stimulating to suppressing plant photosynthesis during the early to peak seasons.

Increased photosynthesis during spring drought in energy-limited ecosystems.

Drought is often thought to reduce ecosystem photosynthesis. However, theory suggests there is potential for increased photosynthesis during meteorological drought, especially in energy-limited ecosystems. Here, we examine the response of photosynthesis (gross primary productivity, GPP) to meteorological drought across the water-energy limitation spectrum. We find a consistent increase in eddy covariance GPP during spring drought in energy-limited ecosystems (83% of the energy-limited sites). Half of spring GPP sensitivity to precipitation was predicted solely from the wetness index (R2 = 0.47, p < 0.001), with weaker relationships in summer and fall. Our results suggest GPP increases during spring drought for 55% of vegetated Northern Hemisphere lands ( >30° N). We then compare these results to terrestrial biosphere model outputs and remote sensing products. In contrast to trends detected in eddy covariance data, model mean GPP always declined under spring precipitation deficits after controlling for air temperature and light availability. While remote sensing products captured the observed negative spring GPP sensitivity in energy-limited ecosystems, terrestrial biosphere models proved insufficiently sensitive to spring precipitation deficits.

AmeriFlux BASE data pipeline to support network growth and data sharing.

AmeriFlux is a network of research sites that measure carbon, water, and energy fluxes between ecosystems and the atmosphere using the eddy covariance technique to study a variety of Earth science questions. AmeriFlux's diversity of ecosystems, instruments, and data-processing routines create challenges for data standardization, quality assurance, and sharing across the network. To address these challenges, the AmeriFlux Management Project (AMP) designed and implemented the BASE data-processing pipeline. The pipeline begins with data uploaded by the site teams, followed by the AMP team's quality assurance and quality control (QA/QC), ingestion of site metadata, and publication of the BASE data product. The semi-automated pipeline enables us to keep pace with the rapid growth of the network. As of 2022, the AmeriFlux BASE data product contains 3,130 site years of data from 444 sites, with standardized units and variable names of more than 60 common variables, representing the largest long-term data repository for flux-met data in the world. The standardized, quality-ensured data product facilitates multisite comparisons, model evaluations, and data syntheses.

Biome-scale temperature sensitivity of ecosystem respiration revealed by atmospheric CO2 observations.

The temperature sensitivity of ecosystem respiration regulates how the terrestrial carbon sink responds to a warming climate but has been difficult to constrain observationally beyond the plot scale. Here we use observations of atmospheric CO2 concentrations from a network of towers together with carbon flux estimates from state-of-the-art terrestrial biosphere models to characterize the temperature sensitivity of ecosystem respiration, as represented by the Arrhenius activation energy, over various North American biomes. We infer activation energies of 0.43 eV for North America and 0.38 eV to 0.53 eV for major biomes therein, which are substantially below those reported for plot-scale studies (approximately 0.65 eV). This discrepancy suggests that sparse plot-scale observations do not capture the spatial-scale dependence and biome specificity of the temperature sensitivity. We further show that adjusting the apparent temperature sensitivity in model estimates markedly improves their ability to represent observed atmospheric CO2 variability. This study provides observationally constrained estimates of the temperature sensitivity of ecosystem respiration directly at the biome scale and reveals that temperature sensitivities at this scale are lower than those based on earlier plot-scale studies. These findings call for additional work to assess the resilience of large-scale carbon sinks to warming.

Global variations in critical drought thresholds that impact vegetation.

Identifying the thresholds of drought that, if crossed, suppress vegetation functioning is vital for accurate quantification of how land ecosystems respond to climate variability and change. We present a globally applicable framework to identify drought thresholds for vegetation responses to different levels of known soil-moisture deficits using four remotely sensed vegetation proxies spanning 2001-2018. The thresholds identified represent critical inflection points for changing vegetation responses from highly resistant to highly vulnerable in response to drought stress, and as a warning signal for substantial vegetation impacts. Drought thresholds varied geographically, with much lower percentiles of soil-moisture anomalies in vegetated areas covered by more forests, corresponding to a comparably stronger capacity to mitigate soil water deficit stress in forested ecosystems. Generally, those lower thresholds are detected in more humid climates. State-of-the-art land models, however, overestimated thresholds of soil moisture (i.e. overestimating drought impacts), especially in more humid areas with higher forest covers and arid areas with few forest covers. Based on climate model projections, we predict that the risk of vegetation damage will increase by the end of the twenty-first century in some hotspots like East Asia, Europe, Amazon, southern Australia and eastern and southern Africa. Our data-based results will inform projections on future drought impacts on terrestrial ecosystems and provide an effective tool for drought management.

Environmental controls on the light use efficiency of terrestrial gross primary production.

Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf-level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy-level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C3 species; however, additional assumptions are required to "scale up" from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed-effects model. The same model design was fitted to site-based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple "goodness of fit" comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

Diminishing seasonality of subtropical water availability in a warmer world dominated by soil moisture-atmosphere feedbacks.

Global warming is expected to cause wet seasons to get wetter and dry seasons to get drier, which would have broad social and ecological implications. However, the extent to which this seasonal paradigm holds over land remains unclear. Here we examine seasonal changes in surface water availability (precipitation minus evaporation, P-E) from CMIP5 and CMIP6 projections. While the P-E seasonal cycle does broadly intensify over much of the land surface, ~20% of land area experiences a diminished seasonal cycle, mostly over subtropical regions and the Amazon. Using land-atmosphere coupling experiments, we demonstrate that 63% of the seasonality reduction is driven by seasonally varying soil moisture (SM) feedbacks on P-E. Declining SM reduces evapotranspiration and modulates circulation to enhance moisture convergence and increase P-E in the dry season but not in the wet season. Our results underscore the importance of SM-atmosphere feedbacks for seasonal water availability changes in a warmer climate.

Increasing sensitivity of dryland vegetation greenness to precipitation due to rising atmospheric CO2.

Water availability plays a critical role in shaping terrestrial ecosystems, particularly in low- and mid-latitude regions. The sensitivity of vegetation growth to precipitation strongly regulates global vegetation dynamics and their responses to drought, yet sensitivity changes in response to climate change remain poorly understood. Here we use long-term satellite observations combined with a dynamic statistical learning approach to examine changes in the sensitivity of vegetation greenness to precipitation over the past four decades. We observe a robust increase in precipitation sensitivity (0.624% yr-1) for drylands, and a decrease (-0.618% yr-1) for wet regions. Using model simulations, we show that the contrasting trends between dry and wet regions are caused by elevated atmospheric CO2 (eCO2). eCO2 universally decreases the precipitation sensitivity by reducing leaf-level transpiration, particularly in wet regions. However, in drylands, this leaf-level transpiration reduction is overridden at the canopy scale by a large proportional increase in leaf area. The increased sensitivity for global drylands implies a potential decrease in ecosystem stability and greater impacts of droughts in these vulnerable ecosystems under continued global change.

Dispersal and fire limit Arctic shrub expansion.

Arctic shrub expansion alters carbon budgets, albedo, and warming rates in high latitudes but remains challenging to predict due to unclear underlying controls. Observational studies and models typically use relationships between observed shrub presence and current environmental suitability (bioclimate and topography) to predict shrub expansion, while omitting shrub demographic processes and non-stationary response to changing climate. Here, we use high-resolution satellite imagery across Alaska and western Canada to show that observed shrub expansion has not been controlled by environmental suitability during 1984-2014, but can only be explained by considering seed dispersal and fire. These findings provide the impetus for better observations of recruitment and for incorporating currently underrepresented processes of seed dispersal and fire in land models to project shrub expansion and climate feedbacks. Integrating these dynamic processes with projected fire extent and climate, we estimate shrubs will expand into 25% of the non-shrub tundra by 2100, in contrast to 39% predicted based on increasing environmental suitability alone. Thus, using environmental suitability alone likely overestimates and misrepresents shrub expansion pattern and its associated carbon sink.

The limits of forest carbon sequestration.

Current models may be overestimating the sequestration potential of forests.

Rising CO2 and warming reduce global canopy demand for nitrogen.

Nitrogen (N) limitation has been considered as a constraint on terrestrial carbon uptake in response to rising CO2 and climate change. By extension, it has been suggested that declining carboxylation capacity (Vcmax ) and leaf N content in enhanced-CO2 experiments and satellite records signify increasing N limitation of primary production. We predicted Vcmax using the coordination hypothesis and estimated changes in leaf-level photosynthetic N for 1982-2016 assuming proportionality with leaf-level Vcmax at 25°C. The whole-canopy photosynthetic N was derived using satellite-based leaf area index (LAI) data and an empirical extinction coefficient for Vcmax , and converted to annual N demand using estimated leaf turnover times. The predicted spatial pattern of Vcmax shares key features with an independent reconstruction from remotely sensed leaf chlorophyll content. Predicted leaf photosynthetic N declined by 0.27% yr-1 , while observed leaf (total) N declined by 0.2-0.25% yr-1 . Predicted global canopy N (and N demand) declined from 1996 onwards, despite increasing LAI. Leaf-level responses to rising CO2 , and to a lesser extent temperature, may have reduced the canopy requirement for N by more than rising LAI has increased it. This finding provides an alternative explanation for declining leaf N that does not depend on increasing N limitation.

Tropical extreme droughts drive long-term increase in atmospheric CO2 growth rate variability.

The terrestrial carbon sink slows the accumulation of carbon dioxide (CO2) in the atmosphere by absorbing roughly 30% of anthropogenic CO2 emissions, but varies greatly from year to year. The resulting variations in the atmospheric CO2 growth rate (CGR) have been related to tropical temperature and water availability. The apparent sensitivity of CGR to tropical temperature ([Formula: see text]) has changed markedly over the past six decades, however, the drivers of the observation to date remains unidentified. Here, we use atmospheric observations, multiple global vegetation models and machine learning products to analyze the cause of the sensitivity change. We found that a threefold increase in [Formula: see text] emerged due to the long-term changes in the magnitude of CGR variability (i.e., indicated by one standard deviation of CGR; STDCGR), which increased 34.7% from 1960-1979 to 1985-2004 and subsequently decreased 14.4% in 1997-2016. We found a close relationship (r2 = 0.75, p < 0.01) between STDCGR and the tropical vegetated area (23°S - 23°N) affected by extreme droughts, which influenced 6-9% of the tropical vegetated surface. A 1% increase in the tropical area affected by extreme droughts led to about 0.14 Pg C yr-1 increase in STDCGR. The historical changes in STDCGR were dominated by extreme drought-affected areas in tropical Africa and Asia, and semi-arid ecosystems. The outsized influence of extreme droughts over a small fraction of vegetated surface amplified the interannual variability in CGR and explained the observed long-term dynamics of [Formula: see text].

CO2 fertilization of terrestrial photosynthesis inferred from site to global scales.

SignificanceThe magnitude of the CO2 fertilization effect on terrestrial photosynthesis is uncertain because it is not directly observed and is subject to confounding effects of climatic variability. We apply three well-established eco-evolutionary optimality theories of gas exchange and photosynthesis, constraining the main processes of CO2 fertilization using measurable variables. Using this framework, we provide robust observationally inferred evidence that a strong CO2 fertilization effect is detectable in globally distributed eddy covariance networks. Applying our method to upscale photosynthesis globally, we find that the magnitude of the CO2 fertilization effect is comparable to its in situ counterpart but highlight the potential for substantial underestimation of this effect in tropical forests for many reflectance-based satellite photosynthesis products.

No evidence for a negative effect of growing season photosynthesis on leaf senescence timing.

The length of the growing season has a large influence on the carbon, water, and energy fluxes of global terrestrial ecosystems. While there has been mounting evidence of an advanced start of the growing season mostly due to elevated spring air temperatures, the mechanisms that control the end of the growing season (EOS) in most ecosystems remain relatively less well understood. Recently, a strong lagged control of EOS by growing season photosynthesis has been proposed, suggesting that more productive growing seasons lead to an earlier EOS. However, this relationship has not been extensively tested with in-situ observations across a variety of ecosystems. Here, we use observations from 40 eddy-covariance flux tower sites in temperate and boreal ecosystems in the northern hemisphere with more than 10 years of observations (594 site-years), ground observations of phenology, satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS), and three leaf senescence models to test the extent of a relationship between growing season photosynthesis and end of season senescence. The results suggest that there is no significant negative relationship between growing season photosynthesis and observed leaf senescence, flux-inferred EOS estimates, or remotely sensed phenological metrics, in most ecosystems. On the contrary, while we found negative effects of summer air temperatures and autumn vapor pressure deficit on EOS, more productive growing seasons were typically related to a later, not earlier, EOS. Our results challenge recent reports of carry-over effects of photosynthesis on EOS timing, and suggest those results may not hold over a large range of ecosystems.

Elevated CO2 moderates the impact of climate change on future bamboo distribution in Madagascar.

The distribution of bamboo is sensitive to climate change and is also potentially affected by increasing atmospheric CO2 concentrations due to its C3 photosynthetic pathway. Yet the effect of CO2 in climate impact assessments of potential changes in bamboo distribution has to date been overlooked. In this study, we proposed a simple and quantitative method to incorporate the impact of atmospheric CO2 concentration into a species distribution modeling framework. To do so, we implemented 10 niche modeling algorithms with regionally downscaled climatic variables and combined field campaign observations. We assessed future climate impacts on the distribution of an economically and ecologically important and widely distributed bamboo species in Madagascar, and examined the effect of increasing CO2 on future projections. Our results suggested that future climatic changes negatively impact potential bamboo distribution in Madagascar, leading to a decline of 34.8% of climatic suitability and a decline of 63.6 ± 3.2% in suitable areas towards 2100 under RCP 8.5. However, increasing atmosphere CO2 offsets the climate impact for bamboo, and led to a smaller reduction of 19.8% in suitability and a potential distribution expansion of +111.6 ± 9.8% in newly suitable areas. We also found that the decline in climatic suitability for bamboo was related to increasing monthly potential evapotranspiration of the warmest quarter and minimum temperature of the warmest month. Conversely, the decreasing isothermality and increasing precipitation of the warmest quarter contributed to projected increase in bamboo-suitable areas. Our study suggested that elevated CO2 may mitigate the decrease in climatic suitability and increase bamboo-suitable areas, through enhancing water use efficiency and decreasing potential evapotranspiration. Our results highlight the importance of accounting for the CO2 effect on future plant species distributions, and provide a mechanistic approach to do so for ecosystems constrained by water.

Exacerbated drought impacts on global ecosystems due to structural overshoot.

Vegetation dynamics are affected not only by the concurrent climate but also by memory-induced lagged responses. For example, favourable climate in the past could stimulate vegetation growth to surpass the ecosystem carrying capacity, leaving an ecosystem vulnerable to climate stresses. This phenomenon, known as structural overshoot, could potentially contribute to worldwide drought stress and forest mortality but the magnitude of the impact is poorly known due to the dynamic nature of overshoot and complex influencing timescales. Here, we use a dynamic statistical learning approach to identify and characterize ecosystem structural overshoot globally and quantify the associated drought impacts. We find that structural overshoot contributed to around 11% of drought events during 1981-2015 and is often associated with compound extreme drought and heat, causing faster vegetation declines and greater drought impacts compared to non-overshoot related droughts. The fraction of droughts related to overshoot is strongly related to mean annual temperature, with biodiversity, aridity and land cover as secondary factors. These results highlight the large role vegetation dynamics play in drought development and suggest that soil water depletion due to warming-induced future increases in vegetation could cause more frequent and stronger overshoot droughts.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.