Carbon-phosphorus cycle models overestimate CO2 enrichment response in a mature Eucalyptus forest.

The importance of phosphorus (P) in regulating ecosystem responses to climate change has fostered P-cycle implementation in land surface models, but their CO2 effects predictions have not been evaluated against measurements. Here, we perform a data-driven model evaluation where simulations of eight widely used P-enabled models were confronted with observations from a long-term free-air CO2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO2 effect on ecosystem carbon (C) sequestration, but they generally overestimated the effects on plant C uptake and growth. We identify leaf-to-canopy scaling of photosynthesis, plant tissue stoichiometry, plant belowground C allocation, and the subsequent consequences for plant-microbial interaction as key areas in which models of ecosystem C-P interaction can be improved. Together, this data-model intercomparison reveals data-driven insights into the performance and functionality of P-enabled models and adds to the existing evidence that the global CO2-driven carbon sink is overestimated by models.

Global critical soil moisture thresholds of plant water stress.

During extensive periods without rain, known as dry-downs, decreasing soil moisture (SM) induces plant water stress at the point when it limits evapotranspiration, defining a critical SM threshold (θcrit). Better quantification of θcrit is needed for improving future projections of climate and water resources, food production, and ecosystem vulnerability. Here, we combine systematic satellite observations of the diurnal amplitude of land surface temperature (dLST) and SM during dry-downs, corroborated by in-situ data from flux towers, to generate the observation-based global map of θcrit. We find an average global θcrit of 0.19 m3/m3, varying from 0.12 m3/m3 in arid ecosystems to 0.26 m3/m3 in humid ecosystems. θcrit simulated by Earth System Models is overestimated in dry areas and underestimated in wet areas. The global observed pattern of θcrit reflects plant adaptation to soil available water and atmospheric demand. Using explainable machine learning, we show that aridity index, leaf area and soil texture are the most influential drivers. Moreover, we show that the annual fraction of days with water stress, when SM stays below θcrit, has increased in the past four decades. Our results have important implications for understanding the inception of water stress in models and identifying SM tipping points.

Enhanced rock weathering increased soil phosphorus availability and altered root phosphorus-acquisition strategies.

Enhanced rock weathering (ERW) has been proposed as a measure to enhance the carbon (C)-sequestration potential and fertility of soils. The effects of this practice on the soil phosphorus (P) pools and the general mechanisms affecting microbial P cycling, as well as plant P uptake are not well understood. Here, the impact of ERW on soil P availability and microbial P cycling functional groups and root P-acquisition traits were explored through a 2-year wollastonite field addition experiment in a tropical rubber plantation. The results show that ERW significantly increased soil microbial carbon-use efficiency and total P concentrations and indirectly increased soil P availability by enhancing organic P mobilization and mineralization of rhizosheath carboxylates and phosphatase, respectively. Also, ERW stimulated the activities of P-solubilizing (gcd, ppa and ppx) and mineralizing enzymes (phoADN and phnAPHLFXIM), thus contributing to the inorganic P solubilization and organic P mineralization. Accompanying the increase in soil P availability, the P-acquisition strategy of the rubber fine roots changed from do-it-yourself acquisition by roots to dependence on mycorrhizal collaboration and the release of root exudates. In addition, the direct effects of ERW on root P-acquisition traits (such as root diameter, specific root length, and mycorrhizal colonization rate) may also be related to changes in the pattern of belowground carbon investments in plants. Our study provides a new insight that ERW increases carbon-sequestration potential and P availability in tropical forests and profoundly affects belowground plant resource-use strategies.

Size, distribution, and vulnerability of the global soil inorganic carbon.

Global estimates of the size, distribution, and vulnerability of soil inorganic carbon (SIC) remain largely unquantified. By compiling 223,593 field-based measurements and developing machine-learning models, we report that global soils store 2305 ± 636 (±1 SD) billion tonnes of carbon as SIC over the top 2-meter depth. Under future scenarios, soil acidification associated with nitrogen additions to terrestrial ecosystems will reduce global SIC (0.3 meters) up to 23 billion tonnes of carbon over the next 30 years, with India and China being the most affected. Our synthesis of present-day land-water carbon inventories and inland-water carbonate chemistry reveals that at least 1.13 ± 0.33 billion tonnes of inorganic carbon is lost to inland-waters through soils annually, resulting in large but overlooked impacts on atmospheric and hydrospheric carbon dynamics.

A global dataset on phosphorus in agricultural soils.

Numerous drivers such as farming practices, erosion, land-use change, and soil biogeochemical background, determine the global spatial distribution of phosphorus (P) in agricultural soils. Here, we revised an approach published earlier (called here GPASOIL-v0), in which several global datasets describing these drivers were combined with a process model for soil P dynamics to reconstruct the past and current distribution of P in cropland and grassland soils. The objective of the present update, called GPASOIL-v1, is to incorporate recent advances in process understanding about soil inorganic P dynamics, in datasets to describe the different drivers, and in regional soil P measurements for benchmarking. We trace the impact of the update on the reconstructed soil P. After the update we estimate a global averaged inorganic labile P of 187 kgP ha-1 for cropland and 91 kgP ha-1 for grassland in 2018 for the top 0-0.3 m soil layer, but these values are sensitive to the mineralization rates chosen for the organic P pools. Uncertainty in the driver estimates lead to coefficients of variation of 0.22 and 0.54 for cropland and grassland, respectively. This work makes the methods for simulating the agricultural soil P maps more transparent and reproducible than previous estimates, and increases the confidence in the new estimates, while the evaluation against regional dataset still suggests rooms for further improvement.

Increased precipitation over land due to climate feedback of large-scale bioenergy cultivation.

Bioenergy with carbon capture and storage (BECCS) is considered to be a key technology for removing carbon dioxide from the atmosphere. However, large-scale bioenergy crop cultivation results in land cover changes and activates biophysical effects on climate, with earth's water recycling altered and energy budget re-adjusted. Here, we use a coupled atmosphere-land model with explicit representations of high-transpiration woody (i.e., eucalypt) and low-transpiration herbaceous (i.e., switchgrass) bioenergy crops to investigate the range of impact of large-scale rainfed bioenergy crop cultivation on the global water cycle and atmospheric water recycling. We find that global land precipitation increases under BECCS scenarios, due to enhanced evapotranspiration and inland moisture advection. Despite enhanced evapotranspiration, soil moisture decreases only slightly, due to increased precipitation and reduced runoff. Our results indicate that, at the global scale, the water consumption by bioenergy crop growth would be partially compensated by atmospheric feedbacks. Thus, to support more effective climate mitigation policies, a more comprehensive assessment, including the biophysical effects of bioenergy cultivation, is highly recommended.

Overestimated nitrogen loss from denitrification for natural terrestrial ecosystems in CMIP6 Earth System Models.

Denitrification and leaching nitrogen (N) losses are poorly constrained in Earth System Models (ESMs). Here, we produce a global map of natural soil 15N abundance and quantify soil denitrification N loss for global natural ecosystems using an isotope-benchmarking method. We show an overestimation of denitrification by almost two times in the 13 ESMs of the Sixth Phase Coupled Model Intercomparison Project (CMIP6, 73 ± 31 Tg N yr-1), compared with our estimate of 38 ± 11 Tg N yr-1, which is rooted in isotope mass balance. Moreover, we find a negative correlation between the sensitivity of plant production to rising carbon dioxide (CO2) concentration and denitrification in boreal regions, revealing that overestimated denitrification in ESMs would translate to an exaggeration of N limitation on the responses of plant growth to elevated CO2. Our study highlights the need of improving the representation of the denitrification in ESMs and better assessing the effects of terrestrial ecosystems on CO2 mitigation.

Soil respiration-driven CO2 pulses dominate Australia's flux variability.

The Australian continent contributes substantially to the year-to-year variability of the global terrestrial carbon dioxide (CO2) sink. However, the scarcity of in situ observations in remote areas prevents the deciphering of processes that force the CO2 flux variability. In this study, by examining atmospheric CO2 measurements from satellites in the period 2009-2018, we find recurrent end-of-dry-season CO2 pulses over the Australian continent. These pulses largely control the year-to-year variability of Australia's CO2 balance. They cause two to three times larger seasonal variations compared with previous top-down inversions and bottom-up estimates. The pulses occur shortly after the onset of rainfall and are driven by enhanced soil respiration preceding photosynthetic uptake in Australia's semiarid regions. The suggested continental-scale relevance of soil-rewetting processes has substantial implications for our understanding and modeling of global climate-carbon cycle feedbacks.

Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

Machine learning for accelerating process-based computation of land biogeochemical cycles.

Global change ecology nowadays embraces ever-growing large observational datasets (big-data) and complex mathematical models that track hundreds of ecological processes (big-model). The rapid advancement of the big-data-big-model has reached its bottleneck: high computational requirements prevent further development of models that need to be integrated over long time-scales to simulate the distribution of ecosystems carbon and nutrient pools and fluxes. Here, we introduce a machine-learning acceleration (MLA) tool to tackle this grand challenge. We focus on the most resource-consuming step in terrestrial biosphere models (TBMs): the equilibration of biogeochemical cycles (spin-up), a prerequisite that can take up to 98% of the computational time. Through three members of the ORCHIDEE TBM family part of the IPSL Earth System Model, including versions that describe the complex interactions between nitrogen, phosphorus and carbon that do not have any analytical solution for the spin-up, we show that an unoptimized MLA reduced the computation demand by 77%-80% for global studies via interpolating the equilibrated state of biogeochemical variables for a subset of model pixels. Despite small biases in the MLA-derived equilibrium, the resulting impact on the predicted regional carbon balance over recent decades is minor. We expect a one-order of magnitude lower computation demand by optimizing the choices of machine learning algorithms, their settings, and balancing the trade-off between quality of MLA predictions and need for TBM simulations for training data generation and bias reduction. Our tool is agnostic to gridded models (beyond TBMs), compatible with existing spin-up acceleration procedures, and opens the door to a wide variety of future applications, with complex non-linear models benefit most from the computational efficiency.

Atmospheric phosphorus deposition amplifies carbon sinks in simulations of a tropical forest in Central Africa.

Spatial redistribution of nutrients by atmospheric transport and deposition could theoretically act as a continental-scale mechanism which counteracts declines in soil fertility caused by nutrient lock-up in accumulating biomass in tropical forests in Central Africa. However, to what extent it affects carbon sinks in forests remains elusive. Here we use a terrestrial biosphere model to quantify the impact of changes in atmospheric nitrogen and phosphorus deposition on plant nutrition and biomass carbon sink at a typical lowland forest site in Central Africa. We find that the increase in nutrient deposition since the 1980s could have contributed to the carbon sink over the past four decades up to an extent which is similar to that from the combined effects of increasing atmospheric carbon dioxide and climate change. Furthermore, we find that the modelled carbon sink responds to changes in phosphorus deposition, but less so to nitrogen deposition. The pronounced response of ecosystem productivity to changes in nutrient deposition illustrates a potential mechanism that could control carbon sinks in Central Africa. Monitoring the quantity and quality of nutrient deposition is needed in this region, given the changes in nutrient deposition due to human land use.

Forest expansion dominates China's land carbon sink since 1980.

Carbon budget accounting relies heavily on Food and Agriculture Organization land-use data reported by governments. Here we develop a new land-use and cover-change database for China, finding that differing historical survey methods biased China's reported data causing large errors in Food and Agriculture Organization databases. Land ecosystem model simulations driven with the new data reveal a strong carbon sink of 8.9 ± 0.8 Pg carbon from 1980 to 2019 in China, which was not captured in Food and Agriculture Organization data-based estimations due to biased land-use and cover-change signals. The land-use and cover-change in China, characterized by a rapid forest expansion from 1980 to 2019, contributed to nearly 44% of the national terrestrial carbon sink. In contrast, climate changes (22.3%), increasing nitrogen deposition (12.9%), and rising carbon dioxide (8.1%) are less important contributors. This indicates that previous studies have greatly underestimated the impact of land-use and cover-change on the terrestrial carbon balance of China. This study underlines the importance of reliable land-use and cover-change databases in global carbon budget accounting.

Process-oriented analysis of dominant sources of uncertainty in the land carbon sink.

The observed global net land carbon sink is captured by current land models. All models agree that atmospheric CO2 and nitrogen deposition driven gains in carbon stocks are partially offset by climate and land-use and land-cover change (LULCC) losses. However, there is a lack of consensus in the partitioning of the sink between vegetation and soil, where models do not even agree on the direction of change in carbon stocks over the past 60 years. This uncertainty is driven by plant productivity, allocation, and turnover response to atmospheric CO2 (and to a smaller extent to LULCC), and the response of soil to LULCC (and to a lesser extent climate). Overall, differences in turnover explain ~70% of model spread in both vegetation and soil carbon changes. Further analysis of internal plant and soil (individual pools) cycling is needed to reduce uncertainty in the controlling processes behind the global land carbon sink.

Convergence in phosphorus constraints to photosynthesis in forests around the world.

Tropical forests take up more carbon (C) from the atmosphere per annum by photosynthesis than any other type of vegetation. Phosphorus (P) limitations to C uptake are paramount for tropical and subtropical forests around the globe. Yet the generality of photosynthesis-P relationships underlying these limitations are in question, and hence are not represented well in terrestrial biosphere models. Here we demonstrate the dependence of photosynthesis and underlying processes on both leaf N and P concentrations. The regulation of photosynthetic capacity by P was similar across four continents. Implementing P constraints in the ORCHIDEE-CNP model, gross photosynthesis was reduced by 36% across the tropics and subtropics relative to traditional N constraints and unlimiting leaf P. Our results provide a quantitative relationship for the P dependence for photosynthesis for the front-end of global terrestrial C models that is consistent with canopy leaf measurements.

Toward a Global Model for Soil Inorganic Phosphorus Dynamics: Dependence of Exchange Kinetics and Soil Bioavailability on Soil Physicochemical Properties.

The representation of phosphorus (P) cycling in global land models remains quite simplistic, particularly on soil inorganic phosphorus. For example, sorption and desorption remain unresolved and their dependence on soil physical and chemical properties is ignored. Empirical parameter values are usually based on expert knowledge or data from few sites with debatable global representativeness in most global land models. To overcome these issues, we compiled from data of inorganic soil P fractions and calculated the fraction of added P remaining in soil solution over time of 147 soil samples to optimize three parameters in a model of soil inorganic P dynamics. The calibrated model performed well (r 2 > 0.7 for 122 soil samples). Model parameters vary by several orders of magnitude, and correlate with soil P fractions of different inorganic pools, soil organic carbon and oxalate extractable metal oxide concentrations among the soil samples. The modeled bioavailability of soil P depends on, not only, the desorption rates of labile and sorbed pool, inorganic phosphorus fractions, the slope of P sorbed against solution P concentration, but also on the ability of biological uptake to deplete solution P concentration and the time scale. The model together with the empirical relationships of model parameters on soil properties can be used to quantify bioavailability of soil inorganic P on various timescale especially when coupled within global land models.

Global maps and factors driving forest foliar elemental composition: the importance of evolutionary history.

Consistent information on the current elemental composition of vegetation at global scale and the variables that determine it is lacking. To fill this gap, we gathered a total of 30 912 georeferenced records on woody plants foliar concentrations of nitrogen (N), phosphorus (P) and potassium (K) from published databases, and produced global maps of foliar N, P and K concentrations for woody plants using neural networks at a resolution of 1 km2 . We used data for climate, atmospheric deposition, soil and morphoclimatic groups to train the neural networks. Foliar N, P and K do not follow clear global latitudinal patterns but are consistent with the hypothesis of soil substrate age. We additionally built generalized linear mixed models to investigate the evolutionary history effect together with the effects of environmental effects. In this comparison, evolutionary history effects explained most of the variability in all cases (mostly > 60%). These results emphasize the determinant role of evolutionary history in foliar elemental composition, which should be incorporated in upcoming dynamic global vegetation models.

Natural forests promote phosphorus retention in soil.

Soil phosphorus (P) availability often limits plant productivity. Classical theories suggest that total P content declines at the temporal scale of pedogenesis, and ecosystems develop toward the efficient use of scarce P during succession. However, the trajectory of ecosystem P within shorter time scales of succession remains unclear. We analyzed changes to P pools at the early (I), middle (II), and late (III) stages of growth of plantation forests (PFs) and the successional stages of natural forests (NFs) at 1969 sites in China. We found significantly lower P contents at later growth stages compared to earlier ones in the PF (p < .05), but higher contents at late successional stages than in earlier stages in the NF (p < .05). Our results indicate that increasing P demand of natural vegetation during succession, may raise, retain, and accumulate P from deeper soil layers. In contrast, ecosystem P in PF was depleted by the more rapidly increasing demand outpacing the development of a P-efficient system. We advocate for more studies to illuminate the mechanisms for determining the divergent changes, which would improve forest management and avoid the vast degradation of PF ecosystems suffering from the ongoing depletion of P.

Is the climate change mitigation effect of enhanced silicate weathering governed by biological processes?

A number of negative emission technologies (NETs) have been proposed to actively remove CO2 from the atmosphere, with enhanced silicate weathering (ESW) as a relatively new NET with considerable climate change mitigation potential. Models calibrated to ESW rates in lab experiments estimate the global potential for inorganic carbon sequestration by ESW at about 0.5-5 Gt CO2  year-1 , suggesting ESW could be an important component of the future NETs mix. In real soils, however, weathering rates may differ strongly from lab conditions. Research on natural weathering has shown that biota such as plants, microbes, and macro-invertebrates can strongly affect weathering rates, but biotic effects were excluded from most ESW lab assessments. Moreover, ESW may alter soil organic carbon sequestration and greenhouse gas emissions by influencing physicochemical and biological processes, which holds the potential to perpetuate even larger negative emissions. Here, we argue that it is likely that the climate change mitigation effect of ESW will be governed by biological processes, emphasizing the need to put these processes on the agenda of this emerging research field.