ALBERT

Description: Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Description: Terrestrial ecosystems have absorbed roughly 30 % of anthropogenic CO2 emissions over the past decades, but it is unclear whether this carbon (C) sink will endure into the future. Despite extensive modeling and experimental and observational studies, what fundamentally determines transient dynamics of terrestrial C storage under global change is still not very clear. Here we develop a new framework for understanding transient dynamics of terrestrial C storage through mathematical analysis and numerical experiments. Our analysis indicates that the ultimate force driving ecosystem C storage change is the C storage capacity, which is jointly determined by ecosystem C input (e.g., net primary production, NPP) and residence time. Since both C input and residence time vary with time, the C storage capacity is time-dependent and acts as a moving attractor that actual C storage chases. The rate of change in C storage is proportional to the C storage potential, which is the difference between the current storage and the storage capacity. The C storage capacity represents instantaneous responses of the land C cycle to external forcing, whereas the C storage potential represents the internal capability of the land C cycle to influence the C change trajectory in the next time step. The influence happens through redistribution of net C pool changes in a network of pools with different residence times. Moreover, this and our other studies have demonstrated that one matrix equation can replicate simulations of most land C cycle models (i.e., physical emulators). As a result, simulation outputs of those models can be placed into a three-dimensional (3-D) parameter space to measure their differences. The latter can be decomposed into traceable components to track the origins of model uncertainty. In addition, the physical emulators make data assimilation computationally feasible so that both C flux- and pool-related datasets can be used to better constrain model predictions of land C sequestration. Overall, this new mathematical framework offers new approaches to understanding, evaluating, diagnosing, and improving land C cycle models.

Description: Carbon (C) turnover time is a key factor in determining C storage capacity in various plant and soil pools and the magnitude of terrestrial C sink in a changing climate. However, the effects of C turnover time on C storage have not been well quantified for previous researches. Here, we first analyzed the difference among different definition of mean turnover time (MTT) including ecosystem MTT(MTTEC) and soil MTT (MTTsoil) and its variability in MTT to climate changes, and then evaluated the changes of ecosystem C storage driven by MTT changes. Our results showed that total GPP-based ecosystem MTT (MTTEC_GPP : 25.0 ± 2.7 years) was shorter than soil MTT (35.5 ± 1.2 years) and NPP-based ecosystem MTT (MTTEC_NPP:50.8 ± 3 years) MTTEC_GPP = Cpool/GPP & MTTsoil = Csoil/NPP & MTTEC_NPP = Cpool/NPP, Cpool and Csoil referring as the ecosystem or soil carbon storage, respectively). At the biome scale, temperature is still the predictor for MTTEC (R2 = 0.77, p 

Description: Carbon (C) turnover time is a key factor in determining C storage capacity in various plant and soil pools as well as terrestrial C sink in a changing climate. However, the effects of C turnover time on ecosystem C storage have not been well explored. In this study, we compared mean C turnover times (MTTs) of ecosystem and soil, examined their variability to climate, and then quantified the spatial variation in ecosystem C storage over time from changes in C turnover time and/or net primary production (NPP). Our results showed that mean ecosystem MTT based on gross primary production (GPP; MTTEC_GPP =  Cpool/GPP, 25.0 ± 2.7 years) was shorter than soil MTT (MTTsoil =  Csoil/NPP, 35.5 ± 1.2 years) and NPP-based ecosystem MTT (MTTEC_NPP =  Cpool/NPP, 50.8 ± 3 years; Cpool and Csoil referred to ecosystem or soil C storage, respectively). On the biome scale, temperature is the best predictor for MTTEC (R2 =  0.77, p 

Description: Ecosystem carbon (C) transit time is a critical diagnostic parameter to characterize land C sequestration. This parameter has different variants in the literature, including a commonly used turnover time. However, we know little about how different transit time and turnover time are in representing carbon cycling through multiple compartments under a non-steady state. In this study, we estimate both C turnover time as defined by the conventional stock over flux and mean C transit time as defined by the mean age of C mass leaving the system. We incorporate them into the Community Atmosphere Biosphere Land Exchange (CABLE) model to estimate C turnover time and transit time in response to climate warming and rising atmospheric [CO2]. Modelling analysis shows that both C turnover time and transit time increase with climate warming but decrease with rising atmospheric [CO2]. Warming increases C turnover time by 2.4 years and transit time by 11.8 years in 2100 relative to that at steady state in 1901. During the same period, rising atmospheric [CO2] decreases C turnover time by 3.8 years and transit time by 5.5 years. Our analysis shows that 65 % of the increase in global mean C transit time with climate warming results from the depletion of fast-turnover C pool. The remaining 35 % increase results from accompanied changes in compartment C age structures. Similarly, the decrease in mean C transit time with rising atmospheric [CO2] results approximately equally from replenishment of C into fast-turnover C pool and subsequent decrease in compartment C age structure. Greatly different from the transit time, the turnover time, which does not account for changes in either C age structure or composition of respired C, underestimated impacts of warming and rising atmospheric [CO2] on C diagnostic time and potentially led to deviations in estimating land C sequestration in multi-compartmental ecosystems.

Description: The concentration–carbon feedback (β), also called the CO2 fertilization effect, is a key unknown in climate–carbon-cycle projections. A better understanding of model mechanisms that govern terrestrial ecosystem responses to elevated CO2 is urgently needed to enable a more accurate prediction of future terrestrial carbon sink. We conducted C-only, carbon–nitrogen (C–N) and carbon–nitrogen–phosphorus (C–N–P) simulations of the Community Atmosphere Biosphere Land Exchange model (CABLE) from 1901 to 2100 with fixed climate to identify the most critical model process that causes divergence in β. We calculated CO2 fertilization effects at various hierarchical levels from leaf biochemical reaction and leaf photosynthesis to canopy gross primary production (GPP), net primary production (NPP), and ecosystem carbon storage (cpool) for seven C3 plant functional types (PFTs) in response to increasing CO2 under the RCP 8.5 scenario. Our results show that β values at biochemical and leaf photosynthesis levels vary little across the seven PFTs, but greatly diverge at canopy and ecosystem levels in all simulations. The low variation of the leaf-level β is consistent with a theoretical analysis that leaf photosynthetic sensitivity to increasing CO2 concentration is almost an invariant function. In the CABLE model, the major jump in variation of β values from leaf levels to canopy and ecosystem levels results from divergence in modeled leaf area index (LAI) within and among PFTs. The correlation of βGPP, βNPP, or βcpool each with βLAI is very high in all simulations. Overall, our results indicate that modeled LAI is a key factor causing the divergence in β in the CABLE model. It is therefore urgent to constrain processes that regulate LAI dynamics in order to better represent the response of ecosystem productivity to increasing CO2 in Earth system models.

Description: Field measurements of aboveground net primary productivity (ANPP) in temperate grasslands suggest that both positive and negative asymmetric responses to changes in precipitation (P) may occur. Under normal range of precipitation variability, wet years typically result in ANPP gains being larger than ANPP declines in dry years (positive asymmetry), whereas increases in ANPP are lower in magnitude in extreme wet years compared to reductions during extreme drought (negative asymmetry). Whether the current generation of ecosystem models with a coupled carbon–water system in grasslands are capable of simulating these asymmetric ANPP responses is an unresolved question. In this study, we evaluated the simulated responses of temperate grassland primary productivity to scenarios of altered precipitation with 14 ecosystem models at three sites: Shortgrass steppe (SGS), Konza Prairie (KNZ) and Stubai Valley meadow (STU), spanning a rainfall gradient from dry to moist. We found that (1) the spatial slopes derived from modeled primary productivity and precipitation across sites were steeper than the temporal slopes obtained from inter-annual variations, which was consistent with empirical data; (2) the asymmetry of the responses of modeled primary productivity under normal inter-annual precipitation variability differed among models, and the mean of the model ensemble suggested a negative asymmetry across the three sites, which was contrary to empirical evidence based on filed observations; (3) the mean sensitivity of modeled productivity to rainfall suggested greater negative response with reduced precipitation than positive response to an increased precipitation under extreme conditions at the three sites; and (4) gross primary productivity (GPP), net primary productivity (NPP), aboveground NPP (ANPP) and belowground NPP (BNPP) all showed concave-down nonlinear responses to altered precipitation in all the models, but with different curvatures and mean values. Our results indicated that most models overestimate the negative drought effects and/or underestimate the positive effects of increased precipitation on primary productivity under normal climate conditions, highlighting the need for improving eco-hydrological processes in those models in the future.

Description: Representations of the terrestrial carbon cycle in land models are becoming increasingly complex. It is crucial to develop approaches for critical assessment of the complex model properties in order to understand key factors contributing to models' performance. In this study, we applied a traceability analysis which decomposes carbon cycle models into traceable components, for two global land models (CABLE and CLM-CASA′) to diagnose the causes of their differences in simulating ecosystem carbon storage capacity. Driven with similar forcing data, CLM-CASA′ predicted  ∼ 31 % larger carbon storage capacity than CABLE. Since ecosystem carbon storage capacity is a product of net primary productivity (NPP) and ecosystem residence time (τE), the predicted difference in the storage capacity between the two models results from differences in either NPP or τE or both. Our analysis showed that CLM-CASA′ simulated 37 % higher NPP than CABLE. On the other hand, τE, which was a function of the baseline carbon residence time (τ′E) and environmental effect on carbon residence time, was on average 11 years longer in CABLE than CLM-CASA′. This difference in τE was mainly caused by longer τ′E of woody biomass (23 vs. 14 years in CLM-CASA′), and higher proportion of NPP allocated to woody biomass (23 vs. 16 %). Differences in environmental effects on carbon residence times had smaller influences on differences in ecosystem carbon storage capacities compared to differences in NPP and τ′E. Overall, the traceability analysis showed that the major causes of different carbon storage estimations were found to be parameters setting related to carbon input and baseline carbon residence times between two models.

Description: Soil organic carbon (SOC) has a significant effect on the carbon emission and climate change. However, current SOC prediction accuracy of most models is very low. Most evaluation studies indicate that the prediction error mainly comes from parameter uncertainties, which can be obviously improved by parameter calibration. Data assimilation technique has been successfully employed for parameter calibration of SOC models. However, data assimilation algorithms such as Bayesian Markov Chain Monte Carlo (MCMC) generally require a large amount of computation cost and are not appropriate for complex global land models. This study proposes a new parameter calibration method based on surrogate optimization techniques for improving the prediction of SOC. Experiments on three types of land carbon cycle models, including Community Land Model with Carnegie-Ames-Stanford Approach biogeochemistry sub-model (CLM-CASA’) and two microbial models show that surrogate-based optimization method is more effective and efficient than MCMC on both accuracy and cost. The root mean squared errors (RMSE) between predictions of different SOC models calibrated by surrogate-base optimization and observations can be reduced up to 12% compared to the results by using Bayesian MCMC. Meanwhile, the corresponding computation cost required is only one thousandth of that with Bayesian MCMC.

Description: Terrestrial ecosystems absorb roughly 30 % of anthropogenic CO2 emissions since preindustrial era, but it is unclear whether this carbon (C) sink will endure into the future. Despite extensive modeling, experimental, and observational studies, what fundamentally determines transient dynamics of terrestrial C storage under climate change is still not very clear. Here we develop a new framework for understanding transient dynamics of terrestrial C storage through mathematical analysis and numerical experiments. Our analysis indicates that the ultimate force driving ecosystem C storage change is the C storage capacity, which is jointly determined by ecosystem C input (e.g., net primary production, NPP) and residence time. Since both C input and residence time vary with time, the C storage capacity is time-dependent and acts as a moving attractor that actual C storage chases. The rate of change in C storage is proportional to the C storage potential, the difference between the current storage and the storage capacity. The C storage capacity represents instantaneous responses of the land C cycle to external forcing, whereas the C storage potential represents the internal capability of the land C cycle to influence the C change trajectory in the next time step. The influence happens through redistribution of net C pool changes in a network of pools with different residence times. Moreover, this and our other studies have demonstrated that one matrix equation can exactly replicate simulations of most land C cycle models (i.e., physical emulators). As a result, simulation outputs of those models can be placed into a three-dimensional (3D) parameter space to measure their differences. The latter can be decomposed into traceable components to track the origins of model uncertainty. Moreover, the emulators make data assimilation computationally feasible so that both C flux- and pool-related datasets can be used to better constrain model predictions of land C sequestration. We also propose that the C storage potential be the targeted variable for research, market trading, and government negotiation for C credits.