Engineering topics

Zaehle, Sönke

Publications and source records attributed to Zaehle, Sönke.

16 records

DOE OSTI2024

The key role of forest disturbance in reconciling estimates of the northern carbon sink

Northern forests are an important carbon sink, but our understanding of the driving factors is limited due to discrepancies between dynamic global vegetation models (DGVMs) and atmospheric inversions. We show that DGVMs simulate a 50% lower sink (1.1 ± 0.5 PgC yr –1 over 2001–2021) across North America, Europe, Russia, and China compared to atmospheric inversions (2.2 ± 0.6 PgC yr –1 ). We explain why DGVMs underestimate the carbon sink by considering how they represent disturbance processes, specifically the overestimation of fire emissions, and the lack of robust forest demography resulting in lower forest regrowth rates than observed. We reconcile net sink estimates by using alternative disturbance-related fluxes. We estimate carbon uptake through forest regrowth by combining satellite-derived forest age and biomass maps. We calculate a regrowth flux of 1.1 ± 0.1 PgC yr –1 , and combine this with satellite-derived estimates of fire emissions (0.4 ± 0.1 PgC yr –1 ), land-use change emissions from bookkeeping models (0.9 ± 0.2 PgC yr –1 ), and the DGVM-estimated sink from CO 2 fertilisation, nitrogen deposition, and climate change (2.2 ± 0.9 PgC yr –1 ). The resulting ‘bottom-up’ net flux of 2.1 ± 0.9 PgC yr –1 agrees with atmospheric inversions. The reconciliation holds at regional scales, increasing confidence in our results.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Trends and Drivers of Terrestrial Sources and Sinks of Carbon Dioxide: An Overview of the TRENDY Project

The terrestrial biosphere plays a major role in the global carbon cycle, and there is a recognized need for regularly updated estimates of land-atmosphere exchange at regional and global scales. An international ensemble of Dynamic Global Vegetation Models (DGVMs), known as the “Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide” (TRENDY) project, quantifies land biophysical exchange processes and biogeochemistry cycles in support of the annual Global Carbon Budget assessments and the REgional Carbon Cycle Assessment and Processes, phase 2 project. DGVMs use a common protocol and set of driving data sets. A set of factorial simulations allows attribution of spatio-temporal changes in land surface processes to three primary global change drivers: changes in atmospheric CO 2 , climate change and variability, and Land Use and Land Cover Changes (LULCC). Here, we describe the TRENDY project, benchmark DGVM performance using remote-sensing and other observational data, and present results for the contemporary period. Simulation results show a large global carbon sink in natural vegetation over 2012–2021, attributed to the CO 2 fertilization effect (3.8 ± 0.8 PgC/yr) and climate (–0.58 ± 0.54 PgC/yr). Forests and semi-arid ecosystems contribute approximately equally to the mean and trend in the natural land sink, and semi-arid ecosystems continue to dominate interannual variability. The natural sink is offset by net emissions from LULCC (–1.6 ± 0.5 PgC/yr), with a net land sink of 1.7 ± 0.6 PgC/yr. Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions.

58 GEOSCIENCES↗

DOE OSTI2024

The Modeled Seasonal Cycles of Surface N 2 O Fluxes and Atmospheric N 2 O

Nitrous oxide (N 2 O) is a greenhouse gas and stratospheric ozone-depleting substance with large and growing anthropogenic emissions. Previous studies identified the influx of N 2 O-depleted air from the stratosphere to partly cause the seasonality in tropospheric N 2 O (aN 2 O), but other contributions remain unclear. Here, we combine surface fluxes from eight land and four ocean models from phase 2 of the Nitrogen/N 2 O Model Intercomparison Project with tropospheric transport modeling to simulate aN 2 O at eight remote air sampling sites for modern and pre-industrial periods. Models show general agreement on the seasonal phasing of zonal-average N 2 O fluxes for most sites, but seasonal peak-to-peak amplitudes differ several-fold across models. The modeled seasonal amplitude of surface aN 2 O ranges from 0.25 to 0.80 ppb (interquartile ranges 21%–52% of median) for land, 0.14–0.25 ppb (17%–68%) for ocean, and 0.28–0.77 ppb (23%–52%) for combined flux contributions. The observed seasonal amplitude ranges from 0.34 to 1.08 ppb for these sites. The stratospheric contributions to aN 2 O, inferred by the difference between the surface-troposphere model and observations, show 16%–126% larger amplitudes and minima delayed by ~1 month compared to Northern Hemisphere site observations. Land fluxes and their seasonal amplitude have increased since the pre-industrial era and are projected to grow further under anthropogenic activities. Our results demonstrate the increasing importance of land fluxes for aN 2 O seasonality. Considering the large model spread, in situ aN 2 O observations and atmospheric transport-chemistry models will provide opportunities for constraining terrestrial and oceanic biosphere models, critical for projecting carbon-nitrogen cycles under ongoing global warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Carbon-phosphorus cycle models overestimate CO 2 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 CO 2 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 CO 2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO 2 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 CO 2 -driven carbon sink is overestimated by models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Estimating the CO 2 Fertilization Effect on Extratropical Forest Productivity From Flux‐Tower Observations

Abstract The land sink of anthropogenic carbon emissions, a crucial component of mitigating climate change, is primarily attributed to the CO 2 fertilization effect on global gross primary productivity (GPP). However, direct observational evidence of this effect remains scarce, hampered by challenges in disentangling the CO 2 fertilization effect from other long‐term confounding drivers, particularly climatic changes. Here, we introduce a novel statistical approach to separate the CO 2 fertilization effect on photosynthetic carbon uptake using eddy covariance (EC) records across 38 extratropical forest sites. We find the median stimulation rate of GPP to be 3.2 ± 0.9 gC m −2 yr −1 ppm −1 (or 16.4 ± 4.2% per 100 ppm) under increasing atmospheric CO 2 across these sites, respectively. To validate the robustness of our findings, we test our statistical method using factorial simulations of an ensemble of process‐based land surface models. We address additional factors, including nitrogen deposition and land management, that may impact plant productivity, potentially confounding the attribution to the CO 2 fertilization effect. Assuming these site‐specific effects offset to some extent across sites as random factors, the estimated median value still reflects the strength of the CO 2 fertilization effect. However, disentanglement of these long‐term effects, often inseparable by timescale, requires further causal research. Our study provides direct evidence that the photosynthetic stimulation is maintained under long‐term CO 2 fertilization across multiple EC sites. Such observation‐based quantification is key to constraining the long‐standing uncertainties in the land carbon cycle under rising CO 2 concentrations.

Environmental Sciences & Ecology↗

DOE OSTI2024

Synthesis of the land carbon fluxes of the Amazon region between 2010 and 2020

The Amazon is the largest continuous tropical forest in the world and plays a key role in the global carbon cycle. Human-induced disturbances and climate change have impacted the Amazon carbon balance. Here we conduct a comprehensive synthesis of existing state-of-the-art estimates of the contemporary land carbon fluxes in the Amazon using a set of bottom-up methods (i.e., dynamic vegetation models and bookkeeping models) and a top-down inversion (atmospheric inversion model) over the Brazilian Amazon and the whole Biogeographical Amazon domain. Over the whole biogeographical Amazon region bottom-up methodologies suggest a small average carbon sink over 2010-2020, in contrast to a small carbon source simulated by top-down inversion (2010-2018). However, these estimates are not significantly different from one another when accounting for their large individual uncertainties, highlighting remaining knowledge gaps, and the urgent need to reduce such uncertainties. Nevertheless, both methodologies agreed that the Brazilian Amazon has been a net carbon source during recent climate extremes and that the south-eastern Amazon was a net land carbon source over the whole study period (2010-2020). Overall, our results point to increasing human-induced disturbances (deforestation and forest degradation by wildfires) and reduction in the old-growth forest sink during drought.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Global Carbon Budget 1800-2018 Loss of Additional Sink Capacity and Present Transient Difference

A dataset of two global carbon fluxes: the loss of additional sink capacity (LASC) and present versus transient difference (PTD). The loss of additional sink capacity is the loss of indirect anthropogenic sink capacity in land ecosystems caused by the loss of ecosystems that can increase carbon in response to changing environments, primarily forests. The present transient difference refers to the direct anthropogenic flux in land ecosystems (from land cover and land use change) and the difference in that flux when calculated assuming "present day" (actually more like 1980) ecosystem carbon stocks or when calculated assuming transient carbon stocks (i.e. caused by environmental change). In this dataset, the LASC and PTD fluxes were calculated using the TRENDY ensemble of Dynamic Global Vegetation Models (DGVMs; Sitch et al., 2024) following the methods described therein and in Obermeier et al. (2021). For further details see Obermeier et al. (2021), Friedlingstein et al. (2019), Sitch et al. (2024), and Walker et al. (2025). The dataset includes ensemble means and standard deviations for annual fluxes, fluxes cumulated annually over the whole time period, and cumulated since 1959. These data were originally collected to compare different approaches (DGVMs versus bookkeeping models) for estimating land use and land cover change (LULCC) emissions (Obermeier et al. 2021). More recently, and similarly, these data were used to harmonize the estimates of the direct and indirect anthropogenic fluxes in the global carbon cycle and to recalculate the global carbon budget (Walker et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A Process‐Model Perspective on Recent Changes in the Carbon Cycle of North America

Continental North America has been found to be a carbon (C) sink over recent decades by multiple studies employing a variety of estimation approaches. However, several key questions and uncertainties remain with these assessments. Here we used results from an ensemble of 19 state-of-the-art dynamic global vegetation models from the TRENDYv9 project to improve these estimates and study the drivers of its interannual variability. Our results show that North America has been a C sink with a magnitude of 0.37 ± 0.38 (mean and one standard deviation) PgC year —1 for the period 2000–2019 (0.31 and 0.44 PgC year —1 in each decade); split into 0.18 ± 0.12 PgC year —1 in Canada (0.15 and 0.20), 0.16 ± 0.17 in the United States (0.14 and 0.17), 0.02 ± 0.05 PgC year —1 in Mexico (0.02 and 0.02) and 0.01 ± 0.02 in Central America and the Caribbean (0.01 and 0.01). About 57% of the new C assimilated by terrestrial ecosystems is allocated into vegetation, 30% into soils, and 13% into litter. Losses of C due to fire account for 41% of the interannual variability of the mean net biome productivity for all North America in the model ensemble. Finally, we show that drought years (e.g., 2002) have the potential to shift the region to a small net C source in the simulations (—0.02 ± 0.46 PgC year —1 ). Our results highlight the importance of identifying the major drivers of the interannual variability of the continental-scale land C cycle along with the spatial distribution of local sink-source dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

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.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

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 CO 2 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 CO 2 (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.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗

DOE OSTI2022

Are Land‐Use Change Emissions in Southeast Asia Decreasing or Increasing?

Abstract Southeast Asia is a region known for active land‐use changes (LUC) over the past 60 years; yet, how trends in net CO 2 uptake and release resulting from LUC activities (net LUC flux) have changed through past decades remains uncertain. The level of uncertainty in net LUC flux from process‐based models is so high that it cannot be concluded that newer estimates are necessarily more reliable than older ones. Here, we examined net LUC flux estimates of Southeast Asia for the 1980s−2010s from older and newer sets of Dynamic Global Vegetation Model simulations (TRENDY v2 and v7, respectively), and forcing data used for running those simulations, along with two book‐keeping estimates (H&N and BLUE). These estimates yielded two contrasting historical LUC transitions, such that TRENDY v2 and H&N showed a transition from increased emissions from the 1980s to 1990s to declining emissions in the 2000s, while TRENDY v7 and BLUE showed the opposite transition. We found that these contrasting transitions originated in the update of LUC forcing data, which reduced the loss of forest area during the 1990s. Further evaluation of remote sensing studies, atmospheric inversions, and the history of forestry and environmental policies in Southeast Asia supported the occurrence of peak emissions in the 1990s and declining thereafter. However, whether LUC emissions continue to decline in Southeast Asia remains uncertain as key processes in recent years, such as conversion of peat forest to oil‐palm plantation, are yet to be represented in the forcing data, suggesting a need for further revision.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Assessing the representation of the Australian carbon cycle in global vegetation models

Australia plays an important role in the global terrestrial carbon cycle on inter-annual timescales. While the Australian continent is included in global assessments of the carbon cycle such as the global carbon budget, the performance of dynamic global vegetation models (DGVMs) over Australia has rarely been evaluated. We assessed simulations of net biome production (NBP) and the carbon stored in vegetation between 1901 to 2018 from 13 DGVMs (TRENDY v8 ensemble). We focused our analysis on Australia's short-term (inter-annual) and long-term (decadal to centennial) terrestrial carbon dynamics. The TRENDY models simulated differing magnitudes of NBP on inter-annual timescales, and these differences resulted in significant differences in long-term vegetation carbon accumulation (–4.7 to 9.5 Pg C). We compared the TRENDY ensemble to several satellite-derived datasets and showed that the spread in the models' simulated carbon storage resulted from varying changes in carbon residence time rather than differences in net carbon uptake. Differences in simulated long-term accumulated NBP between models were mostly due to model responses to land-use change. The DGVMs also simulated different sensitivities to atmospheric carbon dioxide (CO 2 ) concentration, although notably, the models with nutrient cycles did not simulate the smallest NBP response to CO 2 . Our results suggest that a change in the climate forcing did not have a large impact on the carbon cycle on long timescales. However, the inter-annual variability in precipitation drives the year-to-year variability in NBP. We analysed the impact of key modes of climate variability, including the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), on NBP. While the DGVMs agreed on sign of the response of NBP to El Niño and La Niña and to positive and negative IOD events, the magnitude of inter-annual variability in NBP differed strongly between models. In addition, we find that differences in the timing of simulated phenology and fire dynamics are associated with differences in simulated or prescribed vegetation cover and process representation. We further find model disagreement in simulated vegetation carbon, phenology, and apparent carbon residence time, indicating that the models have different types and coverage of vegetation across Australia (whether prescribed or emergent). Our study highlights the need to evaluate parameter assumptions and the key processes that drive vegetation dynamics, such as phenology, mortality, and fire, in an Australian context to reduce uncertainty across models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Dynamic global vegetation models underestimate net CO 2 flux mean and inter-annual variability in dryland ecosystems

Despite their sparse vegetation, dryland regions exert a huge influence over global biogeochemical cycles because they cover more than 40% of the world surface (Schimel 2010 Science 327 418–9). It is thought that drylands dominate the inter-annual variability (IAV) and long-term trend in the global carbon (C) cycle (Poulter et al 2014 Nature 509 600–3, Ahlstrom et al 2015 Science 348 895–9, Zhang et al 2018 Glob. Change Biol. 24 3954–68). Projections of the global land C sink therefore rely on accurate representation of dryland C cycle processes; however, the dynamic global vegetation models (DGVMs) used in future projections have rarely been evaluated against dryland C flux data. Here, we carried out an evaluation of 14 DGVMs (TRENDY v7) against net ecosystem exchange (NEE) data from 12 dryland flux sites in the southwestern US encompassing a range of ecosystem types (forests, shrub- and grasslands). We find that all the models underestimate both mean annual C uptake/release as well as the magnitude of NEE IAV, suggesting that improvements in representing dryland regions may improve global C cycle projections. Across all models, the sensitivity and timing of ecosystem C uptake to plant available moisture was at fault. Spring biases in gross primary production (GPP) dominate the underestimate of mean annual NEE, whereas models' lack of GPP response to water availability in both spring and summer monsoon are responsible for inability to capture NEE IAV. Errors in GPP moisture sensitivity at high elevation forested sites were more prominent during the spring, while errors at the low elevation shrub and grass-dominated sites were more important during the monsoon. We propose a range of hypotheses for why model GPP does not respond sufficiently to changing water availability that can serve as a guide for future dryland DGVM developments. Our analysis suggests that improvements in modeling C cycle processes across more than a quarter of the Earth's land surface could be achieved by addressing the moisture sensitivity of dryland C uptake.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Integrating the evidence for a terrestrial carbon sink caused by increasing atmospheric CO 2

Atmospheric carbon dioxide concentration ([CO 2 ]) is increasing, which increases leaf-scale photosynthesis and intrinsic water-use efficiency. These direct responses have the potential to increase plant growth, vegetation biomass, and soil organic matter; transferring carbon from the atmosphere into terrestrial ecosystems (a carbon sink). A substantial global terrestrial carbon sink would slow the rate of [CO 2 ] increase and thus climate change. However, ecosystem CO 2 -responses are complex or confounded by concurrent changes in multiple agents of global change and evidence for a [CO 2 ]-driven terrestrial carbon sink can appear contradictory. In this work, we synthesise theory and broad, multi-disciplinary evidence for the effects of increasing [CO 2 ] (iCO 2 ) on the global terrestrial carbon sink. Evidence suggests a substantial increase in global photosynthesis since pre-industry. Established theory, supported by experiments, indicates that iCO 2 is likely responsible for about half of the increase. Global carbon budgeting, atmospheric data, and forest inventories indicate a historical carbon sink, and these apparent iCO 2 -responses are high in comparison with experiments and theory. Plant mortality and soil carbon iCO 2 -responses are highly uncertain. In conclusion, a range of evidence supports a positive terrestrial carbon sink in response to iCO 2 , albeit with uncertain magnitude and strong suggestion of a role for additional agents of global change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Organizing principles for vegetation dynamics

Understanding vegetation dynamics is very challenging because of the multitude of contributing processes at 64widely different spatial and temporal scales. In this Perspective we propose that understanding of vegetation dynamics can be improved, permitting better predictions, based on organizing principles that constrain plant and ecosystem be haviour: natural selection, self-organization, and entropy maximization. Although these ideas are increasingly used,a limited common understanding of their theoretical basis has prevented their full potential to be realized. We explain the power of natural selection-based optimality to predict photosynthesis and carbon allocation responses to multiple environmental drivers, and how individual plasticity leads to the predictable self-organization of forest canopies. We show how models of natural selection acting on a few key traits can generate realistic plant communities, and how entropy maximization can distinguish between stochastic and deterministic drivers of vegetation patterns. In combination with empirical exploration of patterns in plant functional variation, these principles can accelerate the development of dynamic vegetation models as well as trait-based ecology resting on strengthened theoretical and empirical foundations.

Geosciences↗