Increasingly negative tropical water-interannual CO2 growth rate coupling.

Terrestrial ecosystems have taken up about 32% of the total anthropogenic CO2 emissions in the past six decades1. Large uncertainties in terrestrial carbon-climate feedbacks, however, make it difficult to predict how the land carbon sink will respond to future climate change2. Interannual variations in the atmospheric CO2 growth rate (CGR) are dominated by land-atmosphere carbon fluxes in the tropics, providing an opportunity to explore land carbon-climate interactions3-6. It is thought that variations in CGR are largely controlled by temperature7-10 but there is also evidence for a tight coupling between water availability and CGR11. Here, we use a record of global atmospheric CO2, terrestrial water storage and precipitation data to investigate changes in the interannual relationship between tropical land climate conditions and CGR under a changing climate. We find that the interannual relationship between tropical water availability and CGR became increasingly negative during 1989-2018 compared to 1960-1989. This could be related to spatiotemporal changes in tropical water availability anomalies driven by shifts in El Niño/Southern Oscillation teleconnections, including declining spatial compensatory water effects9. We also demonstrate that most state-of-the-art coupled Earth System and Land Surface models do not reproduce the intensifying water-carbon coupling. Our results indicate that tropical water availability is increasingly controlling the interannual variability of the terrestrial carbon cycle and modulating tropical terrestrial carbon-climate 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.

Accelerating the energy transition towards photovoltaic and wind in China.

China's goal to achieve carbon (C) neutrality by 2060 requires scaling up photovoltaic (PV) and wind power from 1 to 10-15 PWh year-1 (refs. 1-5). Following the historical rates of renewable installation1, a recent high-resolution energy-system model6 and forecasts based on China's 14th Five-year Energy Development (CFED)7, however, only indicate that the capacity will reach 5-9.5 PWh year-1 by 2060. Here we show that, by individually optimizing the deployment of 3,844 new utility-scale PV and wind power plants coordinated with ultra-high-voltage (UHV) transmission and energy storage and accounting for power-load flexibility and learning dynamics, the capacity of PV and wind power can be increased from 9 PWh year-1 (corresponding to the CFED path) to 15 PWh year-1, accompanied by a reduction in the average abatement cost from US$97 to US$6 per tonne of carbon dioxide (tCO2). To achieve this, annualized investment in PV and wind power should ramp up from US$77 billion in 2020 (current level) to US$127 billion in the 2020s and further to US$426 billion year-1 in the 2050s. The large-scale deployment of PV and wind power increases income for residents in the poorest regions as co-benefits. Our results highlight the importance of upgrading power systems by building energy storage, expanding transmission capacity and adjusting power load at the demand side to reduce the economic cost of deploying PV and wind power to achieve carbon neutrality in China.

Sub-continental-scale carbon stocks of individual trees in African drylands.

The distribution of dryland trees and their density, cover, size, mass and carbon content are not well known at sub-continental to continental scales1-14. This information is important for ecological protection, carbon accounting, climate mitigation and restoration efforts of dryland ecosystems15-18. We assessed more than 9.9 billion trees derived from more than 300,000 satellite images, covering semi-arid sub-Saharan Africa north of the Equator. We attributed wood, foliage and root carbon to every tree in the 0-1,000 mm year-1 rainfall zone by coupling field data19, machine learning20-22, satellite data and high-performance computing. Average carbon stocks of individual trees ranged from 0.54 Mg C ha-1 and 63 kg C tree-1 in the arid zone to 3.7 Mg C ha-1 and 98 kg tree-1 in the sub-humid zone. Overall, we estimated the total carbon for our study area to be 0.84 (±19.8%) Pg C. Comparisons with 14 previous TRENDY numerical simulation studies23 for our area found that the density and carbon stocks of scattered trees have been underestimated by three models and overestimated by 11 models, respectively. This benchmarking can help understand the carbon cycle and address concerns about land degradation24-29. We make available a linked database of wood mass, foliage mass, root mass and carbon stock of each tree for scientists, policymakers, dryland-restoration practitioners and farmers, who can use it to estimate farmland tree carbon stocks from tablets or laptops.

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.

The land-to-ocean loops of the global carbon cycle.

Carbon storage by the ocean and by the land is usually quantified separately, and does not fully take into account the land-to-ocean transport of carbon through inland waters, estuaries, tidal wetlands and continental shelf waters-the 'land-to-ocean aquatic continuum' (LOAC). Here we assess LOAC carbon cycling before the industrial period and perturbed by direct human interventions, including climate change. In our view of the global carbon cycle, the traditional 'long-range loop', which carries carbon from terrestrial ecosystems to the open ocean through rivers, is reinforced by two 'short-range loops' that carry carbon from terrestrial ecosystems to inland waters and from tidal wetlands to the open ocean. Using a mass-balance approach, we find that the pre-industrial uptake of atmospheric carbon dioxide by terrestrial ecosystems transferred to the ocean and outgassed back to the atmosphere amounts to 0.65 ± 0.30 petagrams of carbon per year (±2 sigma). Humans have accelerated the cycling of carbon between terrestrial ecosystems, inland waters and the atmosphere, and decreased the uptake of atmospheric carbon dioxide from tidal wetlands and submerged vegetation. Ignoring these changing LOAC carbon fluxes results in an overestimation of carbon storage in terrestrial ecosystems by 0.6 ± 0.4 petagrams of carbon per year, and an underestimation of sedimentary and oceanic carbon storage. We identify knowledge gaps that are key to reduce uncertainties in future assessments of LOAC fluxes.

Wetland emission and atmospheric sink changes explain methane growth in 2020.

Atmospheric methane growth reached an exceptionally high rate of 15.1 ± 0.4 parts per billion per year in 2020 despite a probable decrease in anthropogenic methane emissions during COVID-19 lockdowns1. Here we quantify changes in methane sources and in its atmospheric sink in 2020 compared with 2019. We find that, globally, total anthropogenic emissions decreased by 1.2 ± 0.1 teragrams of methane per year (Tg CH4 yr-1), fire emissions decreased by 6.5 ± 0.1 Tg CH4 yr-1 and wetland emissions increased by 6.0 ± 2.3 Tg CH4 yr-1. Tropospheric OH concentration decreased by 1.6 ± 0.2 per cent relative to 2019, mainly as a result of lower anthropogenic nitrogen oxide (NOx) emissions and associated lower free tropospheric ozone during pandemic lockdowns2. From atmospheric inversions, we also infer that global net emissions increased by 6.9 ± 2.1 Tg CH4 yr-1 in 2020 relative to 2019, and global methane removal from reaction with OH decreased by 7.5 ± 0.8 Tg CH4 yr-1. Therefore, we attribute the methane growth rate anomaly in 2020 relative to 2019 to lower OH sink (53 ± 10 per cent) and higher natural emissions (47 ± 16 per cent), mostly from wetlands. In line with previous findings3,4, our results imply that wetland methane emissions are sensitive to a warmer and wetter climate and could act as a positive feedback mechanism in the future. Our study also suggests that nitrogen oxide emission trends need to be taken into account when implementing the global anthropogenic methane emissions reduction pledge5.

Delayed use of bioenergy crops might threaten climate and food security.

The potential of mitigation actions to limit global warming within 2 °C (ref. 1) might rely on the abundant supply of biomass for large-scale bioenergy with carbon capture and storage (BECCS) that is assumed to scale up markedly in the future2-5. However, the detrimental effects of climate change on crop yields may reduce the capacity of BECCS and threaten food security6-8, thus creating an unrecognized positive feedback loop on global warming. We quantified the strength of this feedback by implementing the responses of crop yields to increases in growing-season temperature, atmospheric CO2 concentration and intensity of nitrogen (N) fertilization in a compact Earth system model9. Exceeding a threshold of climate change would cause transformative changes in social-ecological systems by jeopardizing climate stability and threatening food security. If global mitigation alongside large-scale BECCS is delayed to 2060 when global warming exceeds about 2.5 °C, then the yields of agricultural residues for BECCS would be too low to meet the Paris goal of 2 °C by 2200. This risk of failure is amplified by the sustained demand for food, leading to an expansion of cropland or intensification of N fertilization to compensate for climate-induced yield losses. Our findings thereby reinforce the urgency of early mitigation, preferably by 2040, to avoid irreversible climate change and serious food crises unless other negative-emission technologies become available in the near future to compensate for the reduced capacity of BECCS.

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.

Changes in above- versus belowground biomass distribution in permafrost regions in response to climate warming.

Permafrost regions contain approximately half of the carbon stored in land ecosystems and have warmed at least twice as much as any other biome. This warming has influenced vegetation activity, leading to changes in plant composition, physiology, and biomass storage in aboveground and belowground components, ultimately impacting ecosystem carbon balance. Yet, little is known about the causes and magnitude of long-term changes in the above- to belowground biomass ratio of plants (η). Here, we analyzed η values using 3,013 plots and 26,337 species-specific measurements across eight sites on the Tibetan Plateau from 1995 to 2021. Our analysis revealed distinct temporal trends in η for three vegetation types: a 17% increase in alpine wetlands, and a decrease of 26% and 48% in alpine meadows and alpine steppes, respectively. These trends were primarily driven by temperature-induced growth preferences rather than shifts in plant species composition. Our findings indicate that in wetter ecosystems, climate warming promotes aboveground plant growth, while in drier ecosystems, such as alpine meadows and alpine steppes, plants allocate more biomass belowground. Furthermore, we observed a threefold strengthening of the warming effect on η over the past 27 y. Soil moisture was found to modulate the sensitivity of η to soil temperature in alpine meadows and alpine steppes, but not in alpine wetlands. Our results contribute to a better understanding of the processes driving the response of biomass distribution to climate warming, which is crucial for predicting the future carbon trajectory of permafrost ecosystems and climate feedback.

A comprehensive quantification of global nitrous oxide sources and sinks.

Nitrous oxide (N2O), like carbon dioxide, is a long-lived greenhouse gas that accumulates in the atmosphere. Over the past 150 years, increasing atmospheric N2O concentrations have contributed to stratospheric ozone depletion1 and climate change2, with the current rate of increase estimated at 2 per cent per decade. Existing national inventories do not provide a full picture of N2O emissions, owing to their omission of natural sources and limitations in methodology for attributing anthropogenic sources. Here we present a global N2O inventory that incorporates both natural and anthropogenic sources and accounts for the interaction between nitrogen additions and the biochemical processes that control N2O emissions. We use bottom-up (inventory, statistical extrapolation of flux measurements, process-based land and ocean modelling) and top-down (atmospheric inversion) approaches to provide a comprehensive quantification of global N2O sources and sinks resulting from 21 natural and human sectors between 1980 and 2016. Global N2O emissions were 17.0 (minimum-maximum estimates: 12.2-23.5) teragrams of nitrogen per year (bottom-up) and 16.9 (15.9-17.7) teragrams of nitrogen per year (top-down) between 2007 and 2016. Global human-induced emissions, which are dominated by nitrogen additions to croplands, increased by 30% over the past four decades to 7.3 (4.2-11.4) teragrams of nitrogen per year. This increase was mainly responsible for the growth in the atmospheric burden. Our findings point to growing N2O emissions in emerging economies-particularly Brazil, China and India. Analysis of process-based model estimates reveals an emerging N2O-climate feedback resulting from interactions between nitrogen additions and climate change. The recent growth in N2O emissions exceeds some of the highest projected emission scenarios3,4, underscoring the urgency to mitigate N2O emissions.

Managing nitrogen to restore water quality in China.

The nitrogen cycle has been radically changed by human activities1. China consumes nearly one third of the world's nitrogen fertilizers. The excessive application of fertilizers2,3 and increased nitrogen discharge from livestock, domestic and industrial sources have resulted in pervasive water pollution. Quantifying a nitrogen 'boundary'4 in heterogeneous environments is important for the effective management of local water quality. Here we use a combination of water-quality observations and simulated nitrogen discharge from agricultural and other sources to estimate spatial patterns of nitrogen discharge into water bodies across China from 1955 to 2014. We find that the critical surface-water quality standard (1.0 milligrams of nitrogen per litre) was being exceeded in most provinces by the mid-1980s, and that current rates of anthropogenic nitrogen discharge (14.5 ± 3.1 megatonnes of nitrogen per year) to fresh water are about 2.7 times the estimated 'safe' nitrogen discharge threshold (5.2 ± 0.7 megatonnes of nitrogen per year). Current efforts to reduce pollution through wastewater treatment and by improving cropland nitrogen management can partially remedy this situation. Domestic wastewater treatment has helped to reduce net discharge by 0.7 ± 0.1 megatonnes in 2014, but at high monetary and energy costs. Improved cropland nitrogen management could remove another 2.3 ± 0.3 megatonnes of nitrogen per year-about 25 per cent of the excess discharge to fresh water. Successfully restoring a clean water environment in China will further require transformational changes to boost the national nutrient recycling rate from its current average of 36 per cent to about 87 per cent, which is a level typical of traditional Chinese agriculture. Although ambitious, such a high level of nitrogen recycling is technologically achievable at an estimated capital cost of approximately 100 billion US dollars and operating costs of 18-29 billion US dollars per year, and could provide co-benefits such as recycled wastewater for crop irrigation and improved environmental quality and ecosystem services.

Climatic and biotic factors influencing regional declines and recovery of tropical forest biomass from the 2015/16 El Niño.

The 2015/16 El Niño brought severe drought and record-breaking temperatures in the tropics. Here, using satellite-based L-band microwave vegetation optical depth, we mapped changes of above-ground biomass (AGB) during the drought and in subsequent years up to 2019. Over more than 60% of drought-affected intact forests, AGB reduced during the drought, except in the wettest part of the central Amazon, where it declined 1 y later. By the end of 2019, only 40% of AGB reduced intact forests had fully recovered to the predrought level. Using random-forest models, we found that the magnitude of AGB losses during the drought was mainly associated with regionally distinct patterns of soil water deficits and soil clay content. For the AGB recovery, we found strong influences of AGB losses during the drought and of [Formula: see text]. [Formula: see text] is a parameter related to canopy structure and is defined as the ratio of two relative height (RH) metrics of Geoscience Laser Altimeter System (GLAS) waveform data-RH25 (25% energy return height) and RH100 (100% energy return height; i.e., top canopy height). A high [Formula: see text] may reflect forests with a tall understory, thick and closed canopy, and/or without degradation. Such forests with a high [Formula: see text] ([Formula: see text] ≥ 0.3) appear to have a stronger capacity to recover than low-[Formula: see text] ones. Our results highlight the importance of forest structure when predicting the consequences of future drought stress in the tropics.

Increasing and widespread vulnerability of intact tropical rainforests to repeated droughts.

Intact tropical rainforests have been exposed to severe droughts in recent decades, which may threaten their integrity, their ability to sequester carbon, and their capacity to provide shelter for biodiversity. However, their response to droughts remains uncertain due to limited high-quality, long-term observations covering extensive areas. Here, we examined how the upper canopy of intact tropical rainforests has responded to drought events globally and during the past 3 decades. By developing a long pantropical time series (1992 to 2018) of monthly radar satellite observations, we show that repeated droughts caused a sustained decline in radar signal in 93%, 84%, and 88% of intact tropical rainforests in the Americas, Africa, and Asia, respectively. Sudden decreases in radar signal were detected around the 1997-1998, 2005, 2010, and 2015 droughts in tropical Americas; 1999-2000, 2004-2005, 2010-2011, and 2015 droughts in tropical Africa; and 1997-1998, 2006, and 2015 droughts in tropical Asia. Rainforests showed similar low resistance (the ability to maintain predrought condition when drought occurs) to severe droughts across continents, but American rainforests consistently showed the lowest resilience (the ability to return to predrought condition after the drought event). Moreover, while the resistance of intact tropical rainforests to drought is decreasing, albeit weakly in tropical Africa and Asia, forest resilience has not increased significantly. Our results therefore suggest the capacity of intact rainforests to withstand future droughts is limited. This has negative implications for climate change mitigation through forest-based climate solutions and the associated pledges made by countries under the Paris Agreement.

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.

Evaluation of effects of heat released from SOC decomposition on soil carbon stock and temperature.

Heat released from soil organic carbon (SOC) decomposition (referred to as microbial heat hereafter) could alter the soil's thermal and hydrological conditions, subsequently modulate SOC decomposition and its feedback with climate. While understanding this feedback is crucial for shaping policy to achieve specific climate goal, it has not been comprehensively assessed. This study employs the ORCHIDEE-MICT model to investigate the effects of microbial heat, referred to as heating effect, focusing on their impacts on SOC accumulation, soil temperature and net primary productivity (NPP), as well as implication on land-climate feedback under two CO2 emissions scenarios (RCP2.6 and RCP8.5). The findings reveal that the microbial heat decreases soil carbon stock, predominantly in upper layers, and elevates soil temperatures, especially in deeper layers. This results in a marginal reduction in global SOC stocks due to accelerated SOC decomposition. Altered seasonal cycles of SOC decomposition and soil temperature are simulated, with the most significant temperature increase per unit of microbial heat (0.31 K J-1) occurring at around 273.15 K (median value of all grid cells where air temperature is around 273.15 K). The heating effect leads to the earlier loss of permafrost area under RCP8.5 and hinders its restoration under RCP2.6 after peak warming. Although elevated soil temperature under climate warming aligns with expectation, the anticipated accelerated SOC decomposition and large amplifying feedback on climate warming were not observed, mainly because of reduced modeled initial SOC stock and limited NPP with heating effect. These underscores the multifaceted impacts of microbial heat. Comprehensive understanding of these effects would be vital for devising effective climate change mitigation strategies in a warming world.

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.

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.

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.