12 Search Results

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects andmore » that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.« less

Abstract Land use change (LUC) alters the global carbon (C) stock, but our estimation of the alteration remains uncertain and is a major impediment to predicting the global C cycle. The uncertainty is partly due to the limited number and geographical bias of observations, and limited exploration of its predictors. Here we generated a comprehensive global database of 5,980 observations from 790 articles. The number of sites evaluated is at least seven times larger than in previous meta‐analyses. Our constrained estimates of different LUC's effects on soil organic C (SOC) and their variations across global climates reveal underestimation/overestimation in previousmore » estimates. Converting forests and grasslands to croplands reduced SOC by 24.5% ± 1.53% (−11.03 ± 1.06 Mg ha⁻¹) and 22.7% ± 1.22% (−8.09 ± 0.67 Mg ha⁻¹), while 28.0% ± 1.56% (4.46 ± 0.42 Mg ha⁻¹) and 33.5% ± 1.68% (5.8 ± 0.38 Mg ha⁻¹) increases, respectively, were obtained in the reverse processes. Converting forests to grasslands decreased SOC by 2.1% ± 1.22% (−1.13 ± 0.44 Mg ha⁻¹), while the reverse process increased SOC by 18.6% ± 1.73% (3.31 ± 0.51 Mg ha⁻¹). Modeled relative importance of 10 drivers of LUC's impact on SOC revealed that higher initial SOC (iSOC) does not solely determine SOC loss in SOC‐negative LUC scenarios as previously proposed. Across four decades, reconverting croplands to forests and grasslands recovered only 49.5% (6.1 ± 0.51 Mg ha⁻¹) and 75.3% (7.0 ± 0.38 Mg ha⁻¹) of the iSOC, respectively, indicating the need for protecting C‐rich ecosystems. Our global data set advances information on LUC's effect on SOC and can be valuable to constrain Earth system models to reliably estimate global SOC stocks and plan climate change mitigation strategies.« less

Abstract Soils store the largest amount of carbon (C) in the biosphere, and the C pool in soil is critical to the global C balance. Numerous microbial models have been developed over the last few decades to represent microbial processes that regulate the responses of soil organic carbon (SOC) to climate change. However, the representation of microbial processes varies, and how microbial processes are incorporated into SOC models has not been well explored. Here, we reviewed 71 microbial models to characterize the microbial processes incorporated into SOC models and analyzed variations in mechanistic complexity. We revealed that (a) four processesmore » (microbial‐mediated decomposition, mineral interaction, microbial necromass recycling, and active and dormant microbial dynamics) are commonly incorporated in microbial models, (b) ∼48% of models simulate only one microbial process (i.e., microbial‐mediated decomposition) and 35% of models simulate two microbial processes: for example, microbial‐mediated decomposition and mineral interaction, (c) more than 80% microbial models use nonlinear equations, such as forward Michaelis‐Menten kinetics, to represent SOC decomposition, (d) the concept of persistence of SOC due to its intrinsic properties has been replaced by organo‐mineral interaction (∼39% of microbial models) that protects SOC from decomposition, and (e) various temperature and moisture modifiers and pH effects have been used to explain the environmental effect on microbial processes. In the future, to realistically incorporate microbial processes into Earth System Models, it is imperative to identify experimental evidence on rate limitation processes and firmly ground model structure on the field and laboratory data.« less

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchangesmore » among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.« less

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. Themore » model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.« less

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140more » balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.« less

Abstract Soil microbial respiration is an important source of uncertainty in projecting future climate and carbon (C) cycle feedbacks. However, its feedbacks to climate warming and underlying microbial mechanisms are still poorly understood. Here we show that the temperature sensitivity of soil microbial respiration ( Q ₁₀ ) in a temperate grassland ecosystem persistently decreases by 12.0 ± 3.7% across 7 years of warming. Also, the shifts of microbial communities play critical roles in regulating thermal adaptation of soil respiration. Incorporating microbial functional gene abundance data into a microbially-enabled ecosystem model significantly improves the modeling performance of soil microbial respiration by 5–19%,more » and reduces model parametric uncertainty by 55–71%. In addition, modeling analyses show that the microbial thermal adaptation can lead to considerably less heterotrophic respiration (11.6 ± 7.5%), and hence less soil C loss. If such microbially mediated dampening effects occur generally across different spatial and temporal scales, the potential positive feedback of soil microbial respiration in response to climate warming may be less than previously predicted.« less

The dynamics of soil phosphorus (P) control its bioavailability. Yet it remains a challenge to quantify soil P dynamics. In this work we developed a soil P dynamics (SPD) model. We then assimilated eight data sets of 426-day changes in Hedley P fractions into the SPD model, to quantify the dynamics of six major P pools in eight soil samples that are representative of a wide type of soils. The performance of our SPD model was better for labile P, secondary mineral P, and occluded P than for nonoccluded organic P (Po) and primary mineral P. All parameters describing soilmore » P dynamics were approximately constrained by the data sets. The average turnover rates were labile P 0.040 g g^-1 day-1, nonoccluded Po 0.051 g g^-1 day^-1, secondary mineral P 0.023 g g^-1 day^-1, primary mineral P 0.00088 g g^-1 day^-1, occluded Po 0.0066 g g^-1 day^-1, and occluded inorganic P 0.0065 g g^-1 day^-1, in the greenhouse environment studied. Labile P was transferred on average more to nonoccluded Po (transfer coefficient of 0.42) and secondary mineral P (0.38) than to plants (0.20). Soil pH and organic C concentration were the key soil properties regulating the competition for P between plants and soil secondary minerals. The turnover rate of labile P was positively correlated with that of nonoccluded Po and secondary mineral P. The pool size of labile P was most sensitive to its turnover rate. Overall, we suggest data assimilation can contribute significantly to an improved understanding of soil P dynamics.« less

Global and regional projections of climate change by Earth system models are limited by their uncertain estimates of terrestrial ecosystem productivity. At the middle to low latitudes, the East Asian monsoon region has higher productivity than forests in Europe–Africa and North America, but its estimate by current generation of terrestrial biosphere models (TBMs) has seldom been systematically evaluated. Here, we developed a traceability framework to evaluate the simulated gross primary productivity (GPP) by 15 TBMs in the East Asian monsoon region. The framework links GPP to net primary productivity, biomass, leaf area and back to GPP via incorporating multiple vegetationmore » functional properties of carbon–use efficiency (CUE), vegetation C turnover time (τ_veg), leaf C fraction (F_leaf), specific leaf area (SLA), and leaf area index (LAI)–level photosynthesis (P_LAI), respectively. We then applied a relative importance algorithm to attribute intermodel variation at each node. The results showed that large intermodel variation in GPP over 1901–2010 were mainly propagated from their different representation of vegetation functional properties. For example, SLA explained 77% of the intermodel difference in leaf area, which contributed 90% to the simulated GPP differences. In addition, the models simulated higher CUE (18.1 ± 21.3%), τ_veg (18.2 ± 26.9%), and SLA (27.4±36.5%) than observations, leading to the overestimation of simulated GPP across the East Asian monsoon region. Furthermore, these results suggest the large uncertainty of current TBMs in simulating GPP is largely propagated from their poor representation of the vegetation functional properties and call for a better understanding of the covariations between plant functional properties in terrestrial ecosystems.« less

The ability to forecast ecological carbon cycling is imperative to land management in a world where past carbon fluxes are no longer a clear guide in the Anthropocene. However, carbon–flux forecasting has not been practiced routinely like numerical weather prediction. This study explored (1) the relative contributions of model forcing data and parameters to uncertainty in forecasting flux– versus pool–based carbon cycle variables and (2) the time points when temperature and CO₂ treatments may cause statistically detectable differences in those variables. We developed an online forecasting workflow (Ecological Platform for Assimilation of Data (EcoPAD)), which facilitates iterative data–model integration. EcoPADmore » automates data transfer from sensor networks, data assimilation, and ecological forecasting. We used the Spruce and Peatland Responses Under Changing Experiments data collected from 2011 to 2014 to constrain the parameters in the Terrestrial Ecosystem Model, forecast carbon cycle responses to elevated CO₂ and a gradient of warming from 2015 to 2024, and specify uncertainties in the model output. Our results showed that data assimilation substantially reduces forecasting uncertainties. Interestingly, we found that the stochasticity of future external forcing contributed more to the uncertainty of forecasting future dynamics of C flux–related variables than model parameters. However, the parameter uncertainty primarily contributes to the uncertainty in forecasting C pool–related response variables. Given the uncertainties in forecasting carbon fluxes and pools, our analysis showed that statistically different responses of fast–turnover pools to various CO₂ and warming treatments were observed sooner than slow–turnover pools. In conclusion, our study has identified the sources of uncertainties in model prediction and thus leads to improve ecological carbon cycling forecasts in the future.« less