A dynamic global vegetation model (DGVM) is applied in a probabilistic framework and benchmarking system to constrain uncertain model parameters by observations and to quantify carbon emissions from ...land-use and land-cover change (LULCC). Processes featured in DGVMs include parameters which are prone to substantial uncertainty. To cope with these uncertainties Latin hypercube sampling (LHS) is used to create a 1000-member perturbed parameter ensemble, which is then evaluated with a diverse set of global and spatiotemporally resolved observational constraints. We discuss the performance of the constrained ensemble and use it to formulate a new best-guess version of the model (LPX-Bern v1.4). The observationally constrained ensemble is used to investigate historical emissions due to LULCC (ELUC) and their sensitivity to model parametrization. We find a global ELUC estimate of 158 (108, 211) PgC (median and 90 % confidence interval) between 1800 and 2016. We compare ELUC to other estimates both globally and regionally. Spatial patterns are investigated and estimates of ELUC of the 10 countries with the largest contribution to the flux over the historical period are reported. We consider model versions with and without additional land-use processes (shifting cultivation and wood harvest) and find that the difference in global ELUC is on the same order of magnitude as parameter-induced uncertainty and in some cases could potentially even be offset with appropriate parameter choice.
Our understanding and quantification of global soil nitrous oxide (N2O) emissions and the underlying processes remain largely uncertain. Here, we assessed the effects of multiple anthropogenic and ...natural factors, including nitrogen fertilizer (N) application, atmospheric N deposition, manure N application, land cover change, climate change, and rising atmospheric CO2 concentration, on global soil N2O emissions for the period 1861–2016 using a standard simulation protocol with seven process‐based terrestrial biosphere models. Results suggest global soil N2O emissions have increased from 6.3 ± 1.1 Tg N2O‐N/year in the preindustrial period (the 1860s) to 10.0 ± 2.0 Tg N2O‐N/year in the recent decade (2007–2016). Cropland soil emissions increased from 0.3 Tg N2O‐N/year to 3.3 Tg N2O‐N/year over the same period, accounting for 82% of the total increase. Regionally, China, South Asia, and Southeast Asia underwent rapid increases in cropland N2O emissions since the 1970s. However, US cropland N2O emissions had been relatively flat in magnitude since the 1980s, and EU cropland N2O emissions appear to have decreased by 14%. Soil N2O emissions from predominantly natural ecosystems accounted for 67% of the global soil emissions in the recent decade but showed only a relatively small increase of 0.7 ± 0.5 Tg N2O‐N/year (11%) since the 1860s. In the recent decade, N fertilizer application, N deposition, manure N application, and climate change contributed 54%, 26%, 15%, and 24%, respectively, to the total increase. Rising atmospheric CO2 concentration reduced soil N2O emissions by 10% through the enhanced plant N uptake, while land cover change played a minor role. Our estimation here does not account for indirect emissions from soils and the directed emissions from excreta of grazing livestock. To address uncertainties in estimating regional and global soil N2O emissions, this study recommends several critical strategies for improving the process‐based simulations.
The ensemble of terrestrial biosphere models indicates that global soil N2O emissions have increased from 6.3 ± 1.1 Tg N2O‐N/year in the preindustrial period (the 1860s) to 10.0 ± 2.0 Tg N2O‐N/year in the recent decade (2007–2016). Cropland soil emissions increased from 0.3 Tg N2O‐N/year to 3.3 Tg N2O‐N/year over the same period, accounting for 82% of the total increase, among which 54% attributes to nitrogen fertilizer application. Regionally, China, South Asia, and Southeast Asia underwent rapid increases in cropland N2O emissions since the 1970s. However, European cropland N2O emissions appear to have decreased by 14%.
Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory ...and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm/yr (6.56×10^4 cu.km/yr) to 617.1 mm/yr (6.87×10^4 cu.km/yr). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm/sq.yr with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm/sq.yr, p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm/yr. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.
Atmospheric vapor pressure deficit (VPD) is a critical variable in determining plant photosynthesis. Synthesis of four global climate datasets reveals a sharp increase of VPD after the late 1990s. In ...response, the vegetation greening trend indicated by a satellite-derived vegetation index (GIMMS3g), which was evident before the late 1990s, was subsequently stalled or reversed. Terrestrial gross primary production derived from two satellite-based models (revised EC-LUE and MODIS) exhibits persistent and widespread decreases after the late 1990s due to increased VPD, which offset the positive CO
fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.
Abstract Carbon dioxide removal (CDR) is discussed for offsetting residual greenhouse gas emissions or even reversing climate change. All emissions scenarios of the Intergovernmental Panel on Climate ...Change that meet the ‘well below 2 °C’ warming target of the Paris Agreement include CDR. Ocean alkalinity enhancement (OAE) may be one possible CDR where the carbon uptake of the ocean is increased by artificial alkalinity addition. Here, we investigate the effect of OAE on modelled carbon reservoirs and fluxes in two observationally-constrained large perturbed parameter ensembles. OAE is assumed to be technically successful and deployed as an additional CDR in the SSP5-3.4 temperature overshoot scenario. Tradeoffs involving feedbacks with atmospheric CO 2 result in a low efficiency of an alkalinity-driven atmospheric CO 2 reduction of −0.35 −0.37 to −0.33 mol C per mol alkalinity addition (skill-weighted mean and 68% c.i.). The realized atmospheric CO 2 reduction, and correspondingly the efficiency, is more than two times smaller than the direct alkalinity-driven enhancement of ocean uptake. The alkalinity-driven ocean carbon uptake is partly offset by the release of carbon from the land biosphere and a reduced ocean carbon sink in response to lowered atmospheric CO 2 under OAE. In a second step we use the Bern3D-LPX model in CO 2 peak-decline simulations to address hysteresis and temporal lags of surface air temperature change (ΔSAT) in an idealized scenario where ΔSAT increases to ~2 °C and then declines to ~1.5 °C as result of CDR. ΔSAT lags the decline in CO 2 -forcing by 18 14–22 years, depending close to linearly on the equilibrium climate sensitivity of the respective ensemble member. These tradeoffs and lags are an inherent feature of the Earth system response to changes in atmospheric CO 2 and will therefore be equally important for other CDR methods.
Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic ...CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: “top‐down” atmospheric inversions and “bottom‐up” biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere–biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink–source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top‐down and bottom‐up approaches and yield more robust knowledge of regional CO2 budgets.
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.
THE GLOBAL N₂O MODEL INTERCOMPARISON PROJECT Tian, Hanqin; Yang, Jia; Lu, Chaoqun ...
Bulletin of the American Meteorological Society,
06/2018, Volume:
99, Issue:
6
Journal Article
Peer reviewed
Open access
Nitrous oxide (N₂O) is an important greenhouse gas and also an ozone-depleting substance that has both natural and anthropogenic sources. Large estimation uncertainty remains on the magnitude and ...spatiotemporal patterns of N₂O fluxes and the key drivers of N₂O production in the terrestrial biosphere. Some terrestrial biosphere models have been evolved to account for nitrogen processes and to show the capability to simulate N₂O emissions from land ecosystems at the global scale, but large discrepancies exist among their estimates primarily because of inconsistent input datasets, simulation protocol, and model structure and parameterization schemes. Based on the consistent model input data and simulation protocol, the global N₂O Model Intercomparison Project (NMIP) was initialized with 10 state-of-the-art terrestrial biosphere models that include nitrogen (N) cycling. Specific objectives of NMIP are to 1) unravel the major N cycling processes controlling N₂O fluxes in each model and identify the uncertainty sources from model structure, input data, and parameters; 2) quantify the magnitude and spatial and temporal patterns of global and regional N₂O fluxes from the preindustrial period (1860) to present and attribute the relative contributions of multiple environmental factors to N₂O dynamics; and 3) provide a benchmarking estimate of N₂O fluxes through synthesizing the multimodel simulation results and existing estimates from ground-based observations, inventories, and statistical and empirical extrapolations. This study provides detailed descriptions for the NMIP protocol, input data, model structure, and key parameters, along with preliminary simulation results. The global and regional N₂O estimation derived from the NMIP is a key component of the global N₂O budget synthesis activity jointly led by the Global Carbon Project and the International Nitrogen Initiative.
Measurements of the stable carbon isotope ratio (δ13C) on annual tree rings offer new opportunities to evaluate mechanisms of variations in photosynthesis and stomatal conductance under changing CO2 ...and climate conditions, especially in conjunction with process-based biogeochemical model simulations. The isotopic discrimination is indicative of the ratio between the CO2 partial pressure in the intercellular cavities and the atmosphere (ci∕ca) and of the ratio of assimilation to stomatal conductance, termed intrinsic water-use efficiency (iWUE). We performed isotope-enabled simulations over the industrial period with the land biosphere module (CLM4.5) of the Community Earth System Model and the Land Surface Processes and Exchanges (LPX-Bern) dynamic global vegetation model. Results for C3 tree species show good agreement with a global compilation of δ13C measurements on leaves, though modeled 13C discrimination by C3 trees is smaller in arid regions than measured. A compilation of 76 tree-ring records, mainly from Europe, boreal Asia, and western North America, suggests on average small 20th century changes in isotopic discrimination and in ci∕ca and an increase in iWUE of about 27 % since 1900. LPX-Bern results match these century-scale reconstructions, supporting the idea that the physiology of stomata has evolved to optimize trade-offs between carbon gain by assimilation and water loss by transpiration. In contrast, CLM4.5 simulates an increase in discrimination and in turn a change in iWUE that is almost twice as large as that revealed by the tree-ring data. Factorial simulations show that these changes are mainly in response to rising atmospheric CO2. The results suggest that the downregulation of ci∕ca and of photosynthesis by nitrogen limitation is possibly too strong in the standard setup of CLM4.5 or that there may be problems associated with the implementation of conductance, assimilation, and related adjustment processes on long-term environmental changes.
The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one‐third of anthropogenic CO2 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.
Plain Language Summary
Earth's natural vegetation absorbs about one‐third of CO2 emissions caused by human activities. This value is produced by a group of models rather than through direct observations. Our study assesses how well models reproduce the processes that drive the CO2 exchange between land and atmosphere using a wide range of data sets that are mainly derived from field measurements and satellite images. These reference data sets are prone to errors that are not quantified in a consistent manner. To account for such errors, we first compare different reference data sets against each other. We then compare model output against reference data and assess whether the differences are comparable to the differences among the reference data sets. We conclude that the performance of models is encouraging given how uncertain reference data are, but that ample potential for improvements remains.
Key Points
Differences between model and observations are often similar compared to differences between independently derived observation‐based data
We quantify differences between independently derived observations to disentangle model deficiencies from observational uncertainties
Future work should address biases in soil organic carbon, leaf area index, and the large spread of gross primary productivity among models