Microbial carbon use efficiency promotes global soil carbon storage.

Soils store more carbon than other terrestrial ecosystems1,2. How soil organic carbon (SOC) forms and persists remains uncertain1,3, which makes it challenging to understand how it will respond to climatic change3,4. It has been suggested that soil microorganisms play an important role in SOC formation, preservation and loss5-7. Although microorganisms affect the accumulation and loss of soil organic matter through many pathways4,6,8-11, microbial carbon use efficiency (CUE) is an integrative metric that can capture the balance of these processes12,13. Although CUE has the potential to act as a predictor of variation in SOC storage, the role of CUE in SOC persistence remains unresolved7,14,15. Here we examine the relationship between CUE and the preservation of SOC, and interactions with climate, vegetation and edaphic properties, using a combination of global-scale datasets, a microbial-process explicit model, data assimilation, deep learning and meta-analysis. We find that CUE is at least four times as important as other evaluated factors, such as carbon input, decomposition or vertical transport, in determining SOC storage and its spatial variation across the globe. In addition, CUE shows a positive correlation with SOC content. Our findings point to microbial CUE as a major determinant of global SOC storage. Understanding the microbial processes underlying CUE and their environmental dependence may help the prediction of SOC feedback to a changing climate.

Soil moisture-atmosphere feedback dominates land carbon uptake variability.

Year-to-year changes in carbon uptake by terrestrial ecosystems have an essential role in determining atmospheric carbon dioxide concentrations1. It remains uncertain to what extent temperature and water availability can explain these variations at the global scale2-5. Here we use factorial climate model simulations6 and show that variability in soil moisture drives 90 per cent of the inter-annual variability in global land carbon uptake, mainly through its impact on photosynthesis. We find that most of this ecosystem response occurs indirectly as soil moisture-atmosphere feedback amplifies temperature and humidity anomalies and enhances the direct effects of soil water stress. The strength of this feedback mechanism explains why coupled climate models indicate that soil moisture has a dominant role4, which is not readily apparent from land surface model simulations and observational analyses2,5. These findings highlight the need to account for feedback between soil and atmospheric dryness when estimating the response of the carbon cycle to climatic change globally5,7, as well as when conducting field-scale investigations of the response of the ecosystem to droughts8,9. Our results show that most of the global variability in modelled land carbon uptake is driven by temperature and vapour pressure deficit effects that are controlled by soil moisture.

Deep learning and process understanding for data-driven Earth system science.

Machine learning approaches are increasingly used to extract patterns and insights from the ever-increasing stream of geospatial data, but current approaches may not be optimal when system behaviour is dominated by spatial or temporal context. Here, rather than amending classical machine learning, we argue that these contextual cues should be used as part of deep learning (an approach that is able to extract spatio-temporal features automatically) to gain further process understanding of Earth system science problems, improving the predictive ability of seasonal forecasting and modelling of long-range spatial connections across multiple timescales, for example. The next step will be a hybrid modelling approach, coupling physical process models with the versatility of data-driven machine learning.

Nitrogen availability and summer drought, but not N:P imbalance, drive carbon use efficiency of a Mediterranean tree-grass ecosystem.

All ecosystems contain both sources and sinks for atmospheric carbon (C). A change in their balance of net and gross ecosystem carbon uptake, ecosystem-scale carbon use efficiency (CUEECO), is a change in their ability to buffer climate change. However, anthropogenic nitrogen (N) deposition is increasing N availability, potentially shifting terrestrial ecosystem stoichiometry towards phosphorus (P) limitation. Depending on how gross primary production (GPP, plants alone) and ecosystem respiration (RECO, plants and heterotrophs) are limited by N, P or associated changes in other biogeochemical cycles, CUEECO may change. Seasonally, CUEECO also varies as the multiple processes that control GPP and respiration and their limitations shift in time. We worked in a Mediterranean tree-grass ecosystem (locally called 'dehesa') characterized by mild, wet winters and summer droughts. We examined CUEECO from eddy covariance fluxes over 6 years under control, +N and + NP fertilized treatments on three timescales: annual, seasonal (determined by vegetation phenological phases) and 14-day aggregations. Finer aggregation allowed consideration of responses to specific patterns in vegetation activity and meteorological conditions. We predicted that CUEECO should be increased by wetter conditions, and successively by N and NP fertilization. Milder and wetter years with proportionally longer growing seasons increased CUEECO, as did N fertilization, regardless of whether P was added. Using a generalized additive model, whole ecosystem phenological status and water deficit indicators, which both varied with treatment, were the main determinants of 14-day differences in CUEECO. The direction of water effects depended on the timescale considered and occurred alongside treatment-dependent water depletion. Overall, future regional trends of longer dry summers may push these systems towards lower CUEECO.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Heatwave breaks down the linearity between sun-induced fluorescence and gross primary production.

Sun-induced fluorescence in the far-red region (SIF) is increasingly used as a remote and proximal-sensing tool capable of tracking vegetation gross primary production (GPP). However, the use of SIF to probe changes in GPP is challenged during extreme climatic events, such as heatwaves. Here, we examined how the 2018 European heatwave (HW) affected the GPP-SIF relationship in evergreen broadleaved trees with a relatively invariant canopy structure. To do so, we combined canopy-scale SIF measurements, GPP estimated from an eddy covariance tower, and active pulse amplitude modulation fluorescence. The HW caused an inversion of the photosynthesis-fluorescence relationship at both the canopy and leaf scales. The highly nonlinear relationship was strongly shaped by nonphotochemical quenching (NPQ), that is, a dissipation mechanism to protect from the adverse effects of high light intensity. During the extreme heat stress, plants experienced a saturation of NPQ, causing a change in the allocation of energy dissipation pathways towards SIF. Our results show the complex modulation of the NPQ-SIF-GPP relationship at an extreme level of heat stress, which is not completely represented in state-of-the-art coupled radiative transfer and photosynthesis models.

Emergence of the physiological effects of elevated CO2 on land-atmosphere exchange of carbon and water.

Elevated atmospheric CO2 (eCO2 ) influences the carbon assimilation rate and stomatal conductance of plants, thereby affecting the global cycles of carbon and water. Yet, the detection of these physiological effects of eCO2 in observational data remains challenging, because natural variations and confounding factors (e.g., warming) can overshadow the eCO2 effects in observational data of real-world ecosystems. In this study, we aim at developing a method to detect the emergence of the physiological CO2 effects on various variables related to carbon and water fluxes. We mimic the observational setting in ecosystems using a comprehensive process-based land surface model QUINCY to simulate the leaf-level effects of increasing atmospheric CO2 concentrations and their century-long propagation through the terrestrial carbon and water cycles across different climate regimes and biomes. We then develop a statistical method based on the signal-to-noise ratio to detect the emergence of the eCO2 effects. The eCO2 effect on gross primary productivity (GPP) emerges at relatively low CO2 increase (∆[CO2 ] ~ 20 ppm) where the leaf area index is relatively high. Compared to GPP, the eCO2 effect causing reduced transpiration water flux (normalized to leaf area) emerges only at relatively high CO2 increase (∆[CO2 ] >> 40 ppm), due to the high sensitivity to climate variability and thus lower signal-to-noise ratio. In general, the response to eCO2 is detectable earlier for variables related to the carbon cycle than the water cycle, when plant productivity is not limited by climatic constraints, and stronger in forest-dominated rather than in grass-dominated ecosystems. Our results provide a step toward when and where we expect to detect physiological CO2 effects in in-situ flux measurements, how to detect them and encourage future efforts to improve the understanding and quantification of these effects in observations of terrestrial carbon and water dynamics.

Mesophyll conductance in land surface models: effects on photosynthesis and transpiration.

The CO2 transfer conductance within plant leaves (mesophyll conductance, gm ) is currently not considered explicitly in most land surface models (LSMs), but instead treated implicitly as an intrinsic property of the photosynthetic machinery. Here, we review approaches to overcome this model deficiency by explicitly accounting for gm , which comprises the re-adjustment of photosynthetic parameters and a model describing the variation of gm in dependence of environmental conditions. An explicit representation of gm causes changes in the response of photosynthesis to environmental factors, foremost leaf temperature, and ambient CO2 concentration, which are most pronounced when gm is small. These changes in leaf-level photosynthesis translate into a stronger climate and CO2 response of gross primary productivity (GPP) and transpiration at the global scale. The results from two independent studies show consistent latitudinal patterns of these effects with biggest differences in GPP in the boreal zone (up to ~15%). Transpiration and evapotranspiration show spatially similar, but attenuated, changes compared with GPP. These changes are indirect effects of gm caused by the assumed strong coupling between stomatal conductance and photosynthesis in current LSMs. Key uncertainties in these simulations are the variation of gm with light and the robustness of its temperature response across plant types and growth conditions. Future research activities focusing on the response of gm to environmental factors and its relation to other plant traits have the potential to improve the representation of photosynthesis in LSMs and to better understand its present and future role in the Earth system.

Partitioning net carbon dioxide fluxes into photosynthesis and respiration using neural networks.

The eddy covariance (EC) technique is used to measure the net ecosystem exchange (NEE) of CO2 between ecosystems and the atmosphere, offering a unique opportunity to study ecosystem responses to climate change. NEE is the difference between the total CO2 release due to all respiration processes (RECO), and the gross carbon uptake by photosynthesis (GPP). These two gross CO2 fluxes are derived from EC measurements by applying partitioning methods that rely on physiologically based functional relationships with a limited number of environmental drivers. However, the partitioning methods applied in the global FLUXNET network of EC observations do not account for the multiple co-acting factors that modulate GPP and RECO flux dynamics. To overcome this limitation, we developed a hybrid data-driven approach based on combined neural networks (NNC-part ). NNC-part incorporates process knowledge by introducing a photosynthetic response based on the light-use efficiency (LUE) concept, and uses a comprehensive dataset of soil and micrometeorological variables as fluxes drivers. We applied the method to 36 sites from the FLUXNET2015 dataset and found a high consistency in the results with those derived from other standard partitioning methods for both GPP (R2  > .94) and RECO (R2  > .8). High consistency was also found for (a) the diurnal and seasonal patterns of fluxes and (b) the ecosystem functional responses. NNC-part performed more realistic than the traditional methods for predicting additional patterns of gross CO2 fluxes, such as: (a) the GPP response to VPD, (b) direct effects of air temperature on GPP dynamics, (c) hysteresis in the diel cycle of gross CO2 fluxes, (d) the sensitivity of LUE to the diffuse to direct radiation ratio, and (e) the post rain respiration pulse after a long dry period. In conclusion, NNC-part is a valid data-driven approach to provide GPP and RECO estimates and complementary to the existing partitioning methods.

Asymmetric responses of ecosystem productivity to rainfall anomalies vary inversely with mean annual rainfall over the conterminous United States.

The CONterminous United States (CONUS) presents a large range of climate conditions and biomes where terrestrial primary productivity and its inter-annual variability are controlled regionally by rainfall and/or temperature. Here, the response of ecosystem productivity to those climate variables was investigated across different biomes from 2010 to 2018 using three climate datasets of precipitation, air temperature or drought severity, combined with several proxies of ecosystem productivity: a remote sensing product of aboveground biomass, an net primary productivity (NPP) remote sensing product, an NPP model-based product and four gross primary productivity products. We used an asymmetry index (AI) where positive AI indicates a greater increase of ecosystem productivity in wet years compared to the decline in dry years, and negative AI indicates a greater decline of ecosystem productivity in dry years compared to the increase in wet years. We found consistent spatial patterns of AI across the CONUS for the different products, with negative asymmetries over the Great Plains and positive asymmetries over the southwestern CONUS. Shrubs and, to a lesser extent, evergreen forests show a persistent positive asymmetry, whilst (natural) grasslands appear to have transitioned from positive to negative anomalies during the last decade. The general tendency of dominant negative asymmetry response for ecosystem productivity across the CONUS appears to be influenced by the negative asymmetry of precipitation anomalies. AI was found to be a function of mean rainfall: more positive AIs were found in dry areas where plants are adapted to drought and take advantage of rainfall pulses, and more negative AIs were found in wet areas, with a threshold delineating the two regimes corresponding to a mean annual rainfall of 200-400 mm/year.

Nutrients and water availability constrain the seasonality of vegetation activity in a Mediterranean ecosystem.

Anthropogenic nitrogen (N) deposition and resulting differences in ecosystem N and phosphorus (P) ratios are expected to impact photosynthetic capacity, that is, maximum gross primary productivity (GPPmax ). However, the interplay between N and P availability with other critical resources on seasonal dynamics of ecosystem productivity remains largely unknown. In a Mediterranean tree-grass ecosystem, we established three landscape-level (24 ha) nutrient addition treatments: N addition (NT), N and P addition (NPT), and a control site (CT). We analyzed the response of ecosystem to altered nutrient stoichiometry using eddy covariance fluxes measurements, satellite observations, and digital repeat photography. A set of metrics, including phenological transition dates (PTDs; timing of green-up and dry-down), slopes during green-up and dry-down period, and seasonal amplitude, were extracted from time series of GPPmax and used to represent the seasonality of vegetation activity. The seasonal amplitude of GPPmax was higher for NT and NPT than CT, which was attributed to changes in structure and physiology induced by fertilization. PTDs were mainly driven by rainfall and exhibited no significant differences among treatments during the green-up period. Yet, both fertilized sites senesced earlier during the dry-down period (17-19 days), which was more pronounced in the NT due to larger evapotranspiration and water usage. Fertilization also resulted in a faster increase in GPPmax during the green-up period and a sharper decline in GPPmax during the dry-down period, with less prominent decline response in NPT. Overall, we demonstrated seasonality of vegetation activity was altered after fertilization and the importance of nutrient-water interaction in such water-limited ecosystems. With the projected warming-drying trend, the positive effects of N fertilization induced by N deposition on GPPmax may be counteracted by an earlier and faster dry-down in particular in areas where the N:P ratio increases, with potential impact on the carbon cycle of water-limited ecosystems.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites.

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate.

Global covariation of carbon turnover times with climate in terrestrial ecosystems.

The response of the terrestrial carbon cycle to climate change is among the largest uncertainties affecting future climate change projections. The feedback between the terrestrial carbon cycle and climate is partly determined by changes in the turnover time of carbon in land ecosystems, which in turn is an ecosystem property that emerges from the interplay between climate, soil and vegetation type. Here we present a global, spatially explicit and observation-based assessment of whole-ecosystem carbon turnover times that combines new estimates of vegetation and soil organic carbon stocks and fluxes. We find that the overall mean global carbon turnover time is 23(+7)(-4) years (95 per cent confidence interval). On average, carbon resides in the vegetation and soil near the Equator for a shorter time than at latitudes north of 75° north (mean turnover times of 15 and 255 years, respectively). We identify a clear dependence of the turnover time on temperature, as expected from our present understanding of temperature controls on ecosystem dynamics. Surprisingly, our analysis also reveals a similarly strong association between turnover time and precipitation. Moreover, we find that the ecosystem carbon turnover times simulated by state-of-the-art coupled climate/carbon-cycle models vary widely and that numerical simulations, on average, tend to underestimate the global carbon turnover time by 36 per cent. The models show stronger spatial relationships with temperature than do observation-based estimates, but generally do not reproduce the strong relationships with precipitation and predict faster carbon turnover in many semi-arid regions. Our findings suggest that future climate/carbon-cycle feedbacks may depend more strongly on changes in the hydrological cycle than is expected at present and is considered in Earth system models.

Effects of mesophyll conductance on vegetation responses to elevated CO2 concentrations in a land surface model.

Mesophyll conductance (gm ) is known to affect plant photosynthesis. However, gm is rarely explicitly considered in land surface models (LSMs), with the consequence that its role in ecosystem and large-scale carbon and water fluxes is poorly understood. In particular, the different magnitudes of gm across plant functional types (PFTs) are expected to cause spatially divergent vegetation responses to elevated CO2 concentrations. Here, an extensive literature compilation of gm across major vegetation types is used to parameterize an empirical model of gm in the LSM JSBACH and to adjust photosynthetic parameters based on simulated An  - Ci curves. We demonstrate that an explicit representation of gm changes the response of photosynthesis to environmental factors, which cannot be entirely compensated by adjusting photosynthetic parameters. These altered responses lead to changes in the photosynthetic sensitivity to atmospheric CO2 concentrations which depend both on the magnitude of gm and the climatic conditions, particularly temperature. We then conducted simulations under ambient and elevated (ambient + 200 µmol/mol) CO2 concentrations for contrasting ecosystems and for historical and anticipated future climate conditions (representative concentration pathways; RCPs) globally. The gm -explicit simulations using the RCP8.5 scenario resulted in significantly higher increases in gross primary productivity (GPP) in high latitudes (+10% to + 25%), intermediate increases in temperate regions (+5% to + 15%), and slightly lower to moderately higher responses in tropical regions (-2% to +5%), which summed up to moderate GPP increases globally. Similar patterns were found for transpiration, but with a lower magnitude. Our results suggest that the effect of an explicit representation of gm is most important for simulated carbon and water fluxes in the boreal zone, where a cold climate coincides with evergreen vegetation.

Climate extremes and the carbon cycle.

The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.

Towards physiologically meaningful water-use efficiency estimates from eddy covariance data.

Intrinsic water-use efficiency (iWUE) characterizes the physiological control on the simultaneous exchange of water and carbon dioxide in terrestrial ecosystems. Knowledge of iWUE is commonly gained from leaf-level gas exchange measurements, which are inevitably restricted in their spatial and temporal coverage. Flux measurements based on the eddy covariance (EC) technique can overcome these limitations, as they provide continuous and long-term records of carbon and water fluxes at the ecosystem scale. However, vegetation gas exchange parameters derived from EC data are subject to scale-dependent and method-specific uncertainties that compromise their ecophysiological interpretation as well as their comparability among ecosystems and across spatial scales. Here, we use estimates of canopy conductance and gross primary productivity (GPP) derived from EC data to calculate a measure of iWUE (G1 , "stomatal slope") at the ecosystem level at six sites comprising tropical, Mediterranean, temperate, and boreal forests. We assess the following six mechanisms potentially causing discrepancies between leaf and ecosystem-level estimates of G1 : (i) non-transpirational water fluxes; (ii) aerodynamic conductance; (iii) meteorological deviations between measurement height and canopy surface; (iv) energy balance non-closure; (v) uncertainties in net ecosystem exchange partitioning; and (vi) physiological within-canopy gradients. Our results demonstrate that an unclosed energy balance caused the largest uncertainties, in particular if it was associated with erroneous latent heat flux estimates. The effect of aerodynamic conductance on G1 was sufficiently captured with a simple representation. G1 was found to be less sensitive to meteorological deviations between canopy surface and measurement height and, given that data are appropriately filtered, to non-transpirational water fluxes. Uncertainties in the derived GPP and physiological within-canopy gradients and their implications for parameter estimates at leaf and ecosystem level are discussed. Our results highlight the importance of adequately considering the sources of uncertainty outlined here when EC-derived water-use efficiency is interpreted in an ecophysiological context.

The response of ecosystem water-use efficiency to rising atmospheric CO2 concentrations: sensitivity and large-scale biogeochemical implications.

Ecosystem water-use efficiency (WUE) is an important metric linking the global land carbon and water cycles. Eddy covariance-based estimates of WUE in temperate/boreal forests have recently been found to show a strong and unexpected increase over the 1992-2010 period, which has been attributed to the effects of rising atmospheric CO2 concentrations on plant physiology. To test this hypothesis, we forced the observed trend in the process-based land surface model JSBACH by increasing the sensitivity of stomatal conductance (gs ) to atmospheric CO2 concentration. We compared the simulated continental discharge, evapotranspiration (ET), and the seasonal CO2 exchange with observations across the extratropical northern hemisphere. The increased simulated WUE led to substantial changes in surface hydrology at the continental scale, including a significant decrease in ET and a significant increase in continental runoff, both of which are inconsistent with large-scale observations. The simulated seasonal amplitude of atmospheric CO2 decreased over time, in contrast to the observed upward trend across ground-based measurement sites. Our results provide strong indications that the recent, large-scale WUE trend is considerably smaller than that estimated for these forest ecosystems. They emphasize the decreasing CO2 sensitivity of WUE with increasing scale, which affects the physiological interpretation of changes in ecosystem WUE.

Plant functional traits and canopy structure control the relationship between photosynthetic CO2 uptake and far-red sun-induced fluorescence in a Mediterranean grassland under different nutrient availability.

Sun-induced fluorescence (SIF) in the far-red region provides a new noninvasive measurement approach that has the potential to quantify dynamic changes in light-use efficiency and gross primary production (GPP). However, the mechanistic link between GPP and SIF is not completely understood. We analyzed the structural and functional factors controlling the emission of SIF at 760 nm (F760 ) in a Mediterranean grassland manipulated with nutrient addition of nitrogen (N), phosphorous (P) or nitrogen-phosphorous (NP). Using the soil-canopy observation of photosynthesis and energy (SCOPE) model, we investigated how nutrient-induced changes in canopy structure (i.e. changes in plant forms abundance that influence leaf inclination distribution function, LIDF) and functional traits (e.g. N content in dry mass of leaves, N%, Chlorophyll a+b concentration (Cab) and maximum carboxylation capacity (Vcmax )) affected the observed linear relationship between F760 and GPP. We conclude that the addition of nutrients imposed a change in the abundance of different plant forms and biochemistry of the canopy that controls F760 . Changes in canopy structure mainly control the GPP-F760 relationship, with a secondary effect of Cab and Vcmax . In order to exploit F760 data to model GPP at the global/regional scale, canopy structural variability, biodiversity and functional traits are important factors that have to be considered.