Global evidence for joint effects of multiple natural and anthropogenic drivers on soil nitrogen cycling.

Global soil nitrogen (N) cycling remains poorly understood due to its complex driving mechanisms. Here, we present a comprehensive analysis of global soil δ15N, a stable isotopic signature indicative of the N input-output balance, using a machine-learning approach on 10,676 observations from 2670 sites. Our findings reveal prevalent joint effects of climatic conditions, plant N-use strategies, soil properties, and other natural and anthropogenic forcings on global soil δ15N. The joint effects of multiple drivers govern the latitudinal distribution of soil δ15N, with more rapid N cycling at lower latitudes than at higher latitudes. In contrast to previous climate-focused models, our data-driven model more accurately simulates spatial changes in global soil δ15N, highlighting the need to consider the joint effects of multiple drivers to estimate the Earth's N budget. These insights contribute to the reconciliation of discordances among empirical, theoretical, and modeling studies on soil N cycling, as well as sustainable N management.

Microbial evolution-An under-appreciated driver of soil carbon cycling.

Although substantial advances in predicting the ecological impacts of global change have been made, predictions of the evolutionary impacts have lagged behind. In soil ecosystems, microbes act as the primary energetic drivers of carbon cycling; however, microbes are also capable of evolving on timescales comparable to rates of global change. Given the importance of soil ecosystems in global carbon cycling, we assess the potential impact of microbial evolution on carbon-climate feedbacks in this system. We begin by reviewing the current state of knowledge concerning microbial evolution in response to global change and its specific effect on soil carbon dynamics. Through this integration, we synthesize a roadmap detailing how to integrate microbial evolution into ecosystem biogeochemical models. Specifically, we highlight the importance of microscale mechanistic soil carbon models, including choosing an appropriate evolutionary model (e.g., adaptive dynamics, quantitative genetics), validating model predictions with 'omics' and experimental data, scaling microbial adaptations to ecosystem level processes, and validating with ecosystem-scale measurements. The proposed steps will require significant investment of scientific resources and might require 10-20 years to be fully implemented. However, through the application of multi-scale integrated approaches, we will advance the integration of microbial evolution into predictive understanding of ecosystems, providing clarity on its role and impact within the broader context of environmental change.

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.

Drivers of natural gas use in U.S. residential buildings.

Natural gas is the primary fuel used in U.S. residences, yet little is known about its consumption patterns and drivers. We use daily county-level gas consumption data to assess the spatial patterns of the relationships and the sensitivities of gas consumption to outdoor air temperature across U.S. households. We fitted linear-plus-plateau functions to daily gas consumption data in 1000 counties, and derived two key coefficients: the heating temperature threshold (Tcrit) and the gas consumption rate change per 1°C temperature drop (Slope). We identified the main predictors of Tcrit and Slope (like income, employment rate, and building type) using interpretable machine learning models built on census data. Finally, we estimated a potential 2.47 million MtCO2 annual emission reduction in U.S. residences by gas savings due to household insulation improvements and hypothetical behavioral change toward reduced consumption by adopting a 1°C lower Tcrit than the current value.

Forest demography and biomass accumulation rates are associated with transient mean tree size vs. density scaling relations.

Linking individual and stand-level dynamics during forest development reveals a scaling relationship between mean tree size and tree density in forest stands, which integrates forest structure and function. However, the nature of this so-called scaling law and its variation across broad spatial scales remain unquantified, and its linkage with forest demographic processes and carbon dynamics remains elusive. In this study, we develop a theoretical framework and compile a broad-scale dataset of long-term sample forest stands (n = 1,433) from largely undisturbed forests to examine the association of temporal mean tree size vs. density scaling trajectories (slopes) with biomass accumulation rates and the sensitivity of scaling slopes to environmental and demographic drivers. The results empirically demonstrate a large variation of scaling slopes, ranging from -4 to -0.2, across forest stands in tropical, temperate, and boreal forest biomes. Steeper scaling slopes are associated with higher rates of biomass accumulation, resulting from a lower offset of forest growth by biomass loss from mortality. In North America, scaling slopes are positively correlated with forest stand age and rainfall seasonality, thus suggesting a higher rate of biomass accumulation in younger forests with lower rainfall seasonality. These results demonstrate the strong association of the transient mean tree size vs. density scaling trajectories with forest demography and biomass accumulation rates, thus highlighting the potential of leveraging forest structure properties to predict forest demography, carbon fluxes, and dynamics at broad spatial scales.

Mitigating nitrogen losses with almost no crop yield penalty during extremely wet years.

Climate change-induced precipitation anomalies during extremely wet years (EWYs) result in substantial nitrogen losses to aquatic ecosystems (Nw). Still, the extent and drivers of these losses, and effective mitigation strategies have remained unclear. By integrating global datasets with well-established crop modeling and machine learning techniques, we reveal notable increases in Nw, ranging from 22 to 56%, during historical EWYs. These pulses are projected to amplify under the SSP126 (SSP370) scenario to 29 to 80% (61 to 120%) due to the projected increases in EWYs and higher nitrogen input. We identify the relative precipitation difference between two consecutive years (diffPr) as the primary driver of extreme Nw. This finding forms the basis of the CLimate Extreme Adaptive Nitrogen Strategy (CLEANS), which scales down nitrogen input adaptively to diffPr, leading to a substantial reduction in extreme Nw with nearly zero yield penalty. Our results have important implications for global environmental sustainability and while safeguarding food security.

Coarse-resolution burned area datasets severely underestimate fire-related forest loss.

Global coarse-resolution (≥250 m) burned area (BA) products have been used to estimate fire related forest loss, but we hypothesised that a significant part of fire impacts might be undetected because of the underestimation of small fires (<100 ha), especially in the tropics. In this paper, we analysed fire-related forest cover loss in sub-Saharan Africa (SSA) for 2016 and 2019 based on a BA product generated from Sentinel-2 data (20 m), which was observed to have significantly lower omission errors than the coarse-resolution BA products. Using these higher resolution BA datasets, we found that fires contribute to >46 % of total forest losses over SSA, more than twice the estimates from coarse-resolution BA products. In addition, burned forest areas showed more than twofold likelihood of subsequent loss compared to unburned ones. In moist tropical forests, the most fire-vulnerable biome, burning had even six times more chance to precede forest loss than unburned areas. We also found that fire-related characteristics, such as fire size and season, and forest fragmentation play a major role in the determination of tree cover fate. Our results reveal that medium-resolution BA detects more fires in late fire season, which tend to have higher impact on forests than early season ones. On the other hand, small fires represented the major driver of forest loss after fires and the vast majority of these losses occur in fragmented landscapes near forest edge (<260 m). Therefore medium-resolution BA products are required to obtain a more accurate evaluation of fire impacts in tropical ecosystems.

Critical role of water conditions in the responses of autumn phenology of marsh wetlands to climate change on the Tibetan Plateau.

The Tibetan Plateau, housing 20% of China's wetlands, plays a vital role in the regional carbon cycle. Examining the phenological dynamics of wetland vegetation in response to climate change is crucial for understanding its impact on the ecosystem. Despite this importance, the specific effects of climate change on wetland vegetation phenology in this region remain uncertain. In this study, we investigated the influence of climate change on the end of the growing season (EOS) of marsh wetland vegetation across the Tibetan Plateau, utilizing satellite-derived Normalized Difference Vegetation Index (NDVI) data and observational climate data. We observed that the regionally averaged EOS of marsh vegetation across the Tibetan Plateau was significantly (p < .05) delayed by 4.10 days/decade from 2001 to 2020. Warming preseason temperatures were found to be the primary driver behind the delay in the EOS of marsh vegetation, whereas preseason cumulative precipitation showed no significant impact. Interestingly, the responses of EOS to climate change varied spatially across the plateau, indicating a regulatory role for hydrological conditions in marsh phenology. In the humid and cold central regions, preseason daytime warming significantly delayed the EOS. However, areas with lower soil moisture exhibited a weaker or reversed delay effect, suggesting complex interplays between temperature, soil moisture, and EOS. Notably, in the arid southwestern regions of the plateau, increased preseason rainfall directly delayed the EOS, while higher daytime temperatures advanced it. Our results emphasize the critical role of hydrological conditions, specifically soil moisture, in shaping marsh EOS responses in different regions. Our findings underscore the need to incorporate hydrological factors into terrestrial ecosystem models, particularly in cold and dry regions, for accurate predictions of marsh vegetation phenological responses to climate change. This understanding is vital for informed conservation and management strategies in the face of current and future climate challenges.

Anthropogenic disturbance exacerbates resilience loss in the Amazon rainforests.

Uncovering the mechanisms that lead to Amazon forest resilience variations is crucial to predict the impact of future climatic and anthropogenic disturbances. Here, we apply a previously used empirical resilience metrics, lag-1 month temporal autocorrelation (TAC), to vegetation optical depth data in C-band (a good proxy of the whole canopy water content) in order to explore how forest resilience variations are impacted by human disturbances and environmental drivers in the Brazilian Amazon. We found that human disturbances significantly increase the risk of critical transitions, and that the median TAC value is ~2.4 times higher in human-disturbed forests than that in intact forests, suggesting a much lower resilience in disturbed forests. Additionally, human-disturbed forests are less resilient to land surface heat stress and atmospheric water stress than intact forests. Among human-disturbed forests, forests with a more closed and thicker canopy structure, which is linked to a higher forest cover and a lower disturbance fraction, are comparably more resilient. These results further emphasize the urgent need to limit deforestation and degradation through policy intervention to maintain the resilience of the Amazon rainforests.

Temperature dependence of spring carbon uptake in northern high latitudes during the past four decades.

In the northern high latitudes, warmer spring temperatures generally lead to earlier leaf onsets, higher vegetation production, and enhanced spring carbon uptake. Yet, whether this positive linkage has diminished under climate change remains debated. Here, we used atmospheric CO2 measurements at Barrow (Alaska) during 1979-2020 to investigate the strength of temperature dependence of spring carbon uptake reflected by two indicators, spring zero-crossing date (SZC) and CO2 drawdown (SCC). We found a fall and rise in the interannual correlation of temperature with SZC and SCC (RSZC-T and RSCC-T ), showing a recent reversal of the previously reported weakening trend of RSZC-T and RSCC-T . We used a terrestrial biosphere model coupled with an atmospheric transport model to reproduce this fall and rise phenomenon and conducted factorial simulations to explore its potential causes. We found that a strong-weak-strong spatial synchrony of spring temperature anomalies per se has contributed to the fall and rise trend in RSZC-T and RSCC-T , despite an overall unbroken temperature control on net ecosystem CO2 fluxes at local scale. Our results provide an alternative explanation for the apparent drop of RSZC-T and RSCC-T during the late 1990s and 2000s, and suggest a continued positive linkage between spring carbon uptake and temperature during the past four decades. We thus caution the interpretation of apparent climate sensitivities of carbon cycle retrieved from spatially aggregated signals.

A meta-analysis highlights globally widespread potassium limitation in terrestrial ecosystems.

Potassium (K+ ) is the most abundant inorganic cation in plant cells, playing a critical role in various plant functions. However, the impacts of K on natural terrestrial ecosystems have been less studied compared with nitrogen (N) and phosphorus (P). Here, we present a global meta-analysis aimed at quantifying the response of aboveground production to K addition. This analysis is based on 144 field K fertilization experiments. We also investigate the influences of climate, soil properties, ecosystem types, and fertilizer regimes on the responses of aboveground production. We find that: K addition significantly increases aboveground production by 12.3% (95% CI: 7.4-17.5%), suggesting a widespread occurrence of K limitation across terrestrial ecosystems; K limitation is more prominent in regions with humid climates, acidic soils, or weathered soils; the effect size of K addition varies among climate zones/regions, and is influenced by multiple factors; and previous N : K and K : P thresholds utilized to detect K limitation in wetlands cannot be applied to other biomes. Our findings emphasize the role of K in limiting terrestrial productivity, which should be integrated into future terrestrial ecosystems models.

A Gradient-Descent Optimization of CO2-CO-NOx Emissions over the Paris Megacity─The Case of the First SARS-CoV-2 Lockdown.

Urban greenhouse gas emissions monitoring is essential to assessing the impact of climate mitigation actions. Using atmospheric continuous measurements of air quality and carbon dioxide (CO2), we developed a gradient-descent optimization system to estimate emissions of the city of Paris. We evaluated our joint CO2-CO-NOx optimization over the first SARS-CoV-2 related lockdown period, resulting in a decrease in emissions by 40% for NOx and 30% for CO2, in agreement with preliminary estimates using bottom-up activity data yet lower than the decrease estimates from Bayesian atmospheric inversions (50%). Before evaluating the model, we first provide an in-depth analysis of three emission data sets. A general agreement in the totals is observed over the region surrounding Paris (known as Île-de-France) since all the data sets are constrained by the reported national and regional totals. However, the data sets show disagreements in their sector distributions as well as in the interspecies ratios. The seasonality also shows disagreements among emission products related to nonindustrial stationary combustion (residential and tertiary combustion). The results presented in this paper show that a multispecies approach has the potential to provide sectoral information to monitor CO2 emissions over urban areas enabled by the deployment of collocated atmospheric greenhouse gases and air quality monitoring stations.

Including land management in a European carbon model with lateral transfer to the oceans.

The use of cover crops (CCs) is a promising cropland management practice with multiple benefits, notably in reducing soil erosion and increasing soil organic carbon (SOC) storage. However, the current ability to represent these factors in land surface models remains limited to small scales or simplified and lumped approaches due to the lack of a sediment-carbon erosion displacement scheme. This precludes a thorough understanding of the consequences of introducing a CC into agricultural systems. In this work, this problem was addressed in two steps with the spatially distributed CE-DYNAM model. First, the historical effect of soil erosion, transport, and deposition on the soil carbon budget at a continental scale in Europe was characterized since the early industrial era, using reconstructed climate and land use forcings. Then, the impact of two distinct policy-oriented scenarios for the introduction of CCs were evaluated, covering the European cropping systems where surface erosion rates or nitrate susceptibility are critical. The evaluation focused on the increase in SOC storage and the export of particulate organic carbon (POC) to the oceans, compiling a continental-scale carbon budget. The results indicated that Europe exported 1.95 TgC/year of POC to the oceans in the last decade, and that CCs can contribute to reducing this amount while increasing SOC storage. Compared to the simulation without CCs, the additional rate of SOC storage induced by CCs peaked after 10 years of their adoption, followed by a decrease, and the cumulative POC export reduction stabilized after around 13 years. The findings indicate that the impacts of CCs on SOC and reduced POC export are persistent regardless of their spatial allocation adopted in the scenarios. Together, the results highlight the importance of taking the temporal aspect of CC adoption into account and indicate that CCs alone are not sufficient to meet the targets of the 4‰ initiative. Despite some known model limitations, which include the lack of feedback of erosion on the net primary productivity and the representation of carbon fluxes with an emulator, the current work constitutes the first approach to successfully couple a distributed routing scheme of eroded carbon to a land carbon model emulator at a reasonably high resolution and continental scale. SHORT ABSTRACT: A spatially distributed model coupling erosion, transport, and deposition to the carbon cycle was developed. Then, it was used to simulate the impact of cover crops on both erosion and carbon, to show that cover crops can simultaneously increase organic carbon storage and reduce particulate organic carbon export to the oceans. The results seemed persistent regardless of the spatial distribution of cover crops.

High-resolution assessment of coal mining methane emissions by satellite in Shanxi, China.

Accurate assessment of coal mine methane (CMM) emissions is a prerequisite for defining baselines and assessing the effectiveness of mitigation measures. Such an endeavor is jeopardized, however, by large uncertainties in current CMM estimates. Here, we assimilated atmospheric methane column concentrations observed by the TROPOMI space borne instrument in a high-resolution regional inversion to estimate CMM emissions in Shanxi, a province representing 15% of the global coal production. The emissions are estimated to be 8.5 ± 0.6 and 8.6 ± 0.6 Tg CH4 yr-1 in 2019 and 2020, respectively, close to upper bound of current bottom-up estimates. Data from more than a thousand of individual mines indicate that our estimated emission factors increase significantly with coal mining depth at prefecture level, suggesting that ongoing deeper mining will increase CMM emission intensity. Our results show robustness of estimating CMM emissions utilizing TROPOMI images and highlight potential of monitoring methane leakages and emissions from satellites.

Temperature extremes of 2022 reduced carbon uptake by forests in Europe.

The year 2022 saw record breaking temperatures in Europe during both summer and fall. Similar to the recent 2018 drought, close to 30% (3.0 million km2) of the European continent was under severe summer drought. In 2022, the drought was located in central and southeastern Europe, contrasting the Northern-centered 2018 drought. We show, using multiple sets of observations, a reduction of net biospheric carbon uptake in summer (56-62 TgC) over the drought area. Specific sites in France even showed a widespread summertime carbon release by forests, additional to wildfires. Partial compensation (32%) for the decreased carbon uptake due to drought was offered by a warm autumn with prolonged biospheric carbon uptake. The severity of this second drought event in 5 years suggests drought-induced reduced carbon uptake to no longer be exceptional, and important to factor into Europe's developing plans for net-zero greenhouse gas emissions that rely on carbon uptake by forests.

Non-symmetric responses of leaf onset date to natural warming and cooling in northern ecosystems.

The northern hemisphere has experienced regional cooling, especially during the global warming hiatus (1998-2012) due to ocean energy redistribution. However, the lack of studies about the natural cooling effects hampers our understanding of vegetation responses to climate change. Using 15,125 ground phenological time series at 3,620 sites since the 1950s and 31-year satellite greenness observations (1982-2012) covering the warming hiatus period, we show a stronger response of leaf onset date (LOD) to natural cooling than to warming, i.e. the delay of LOD caused by 1°C cooling is larger than the advance of LOD with 1°C warming. This might be because cooling leads to larger chilling accumulation and heating requirements for leaf onset, but this non-symmetric LOD response is partially offset by warming-related drying. Moreover, spring greening magnitude, in terms of satellite-based greenness and productivity, is more sensitive to LOD changes in the warming area than in the cooling. These results highlight the importance of considering non-symmetric responses of spring greening to warming and cooling when predicting vegetation-climate feedbacks.