Carbon monoxide (CO) plays an important role in controlling the oxidizing capacity of the atmosphere by reacting with OH radicals that affect atmospheric methane (CH4) dynamics. We develop a process-based biogeochemistry model to quantify the CO exchange between soils and the atmosphere with a 5 min internal time step at the global scale. The model is parameterized using the CO flux data from the field and laboratory experiments for 11 representative ecosystem types. The model is then extrapolated to global terrestrial ecosystems using monthly climate forcing data. Global soil gross consumption, gross production, and net flux of the atmospheric CO are estimated to be from-197 to-180, 34 to 36, and-163 to-145 Tg CO yr-1 (1 Tg Combining double low line 1012 g), respectively, when the model is driven with satellite-based atmospheric CO concentration data during 2000-2013. Tropical evergreen forest, savanna and deciduous forest areas are the largest sinks at 123 Tg CO yr-1. The soil CO gross consumption is sensitive to air temperature and atmospheric CO concentration, while the gross production is sensitive to soil organic carbon (SOC) stock and air temperature. By assuming that the spatially distributed atmospheric CO concentrations (∼ 128 ppbv) are not changing over time, the global mean CO net deposition velocity is estimated to be 0.16-0.19 mm s-1 during the 20th century. Under the future climate scenarios, the CO deposition velocity will increase at a rate of 0.0002-0.0013 mm s-1 yr-1 during 2014-2100, reaching 0.20-0.30 mm s-1 by the end of the 21st century, primarily due to the increasing temperature. Areas near the Equator, the eastern US, Europe and eastern Asia will be the largest sinks due to optimum soil moisture and high temperature. The annual global soil net flux of atmospheric CO is primarily controlled by air temperature, soil temperature, SOC and atmospheric CO concentrations, while its monthly variation is mainly determined by air temperature, precipitation, soil temperature and soil moisture.

Cold region ecosystems store vast amounts of soil organic carbon (C), which upon warming and decomposition can affect the C balance and potentially change these ecosystems from C sinks to carbon dioxide (CO2) sources. We quantified the decadal year-round CO2 flux from an alpine steppe-ecosystem on the Tibetan Plateau using eddy covariance and automatic chamber approaches during a period of significant warming (0.13°C per 10 years; and 0.18°C in the non-growing season alone: 1st October to next 30th April). The results showed that ongoing climate change, mainly warming within the topsoil layers, is the main reason for the site’s change from a sink for to a source of CO2 in the atmosphere. Non-growing-season ecosystem respiration accounted for 51% of the annual ecosystem respiration and has increased significantly. The growing seasons (1st May to 30th September) were consistent CO2 sink periods without significant changes over the study period. Observations revealed high-emission events from the end of the non-growing season to early in the growing season (1st March to fifteenth May), which significantly (p < 0.01) increased at a rate of 22.6 g C m−2 decade−1, ranging from 14.6 ± 10.7 g C m−2 yr−1 in 2012 to 35.3 ± 12.1 g C m−2 yr−1 in 2017. Structural equation modeling suggested that active layer warming was the key factor in explaining changes in ecosystem respiration, leading to significant changes in net ecosystem exchange over the period 2011–2020 and indicated that these changes have already transformed the ecosystem from a CO2 sink into a source. These results can be used to improve our understanding of the sensitivity of ecosystem respiration to increased warming during the non-growing period.

Accurate and cost-effective quantification of the carbon cycle for agroecosystems at decision-relevant scales is critical to mitigating climate change and ensuring sustainable food production. However, conventional process-based or data-driven modeling approaches alone have large prediction uncertainties due to the complex biogeochemical processes to model and the lack of observations to constrain many key state and flux variables. Here we propose a Knowledge-Guided Machine Learning (KGML) framework that addresses the above challenges by integrating knowledge embedded in a process-based model, high-resolution remote sensing observations, and machine learning (ML) techniques. Using the U.S. Corn Belt as a testbed, we demonstrate that KGML can outperform conventional process-based and black-box ML models in quantifying carbon cycle dynamics. Our high-resolution approach quantitatively reveals 86% more spatial detail of soil organic carbon changes than conventional coarse-resolution approaches. Moreover, we outline a protocol for improving KGML via various paths, which can be generalized to develop hybrid models to better predict complex earth system dynamics.

Agricultural nitrous oxide (N2O) emission accounts for a non-trivial fraction of global greenhouse gas (GHG) budget. To date, estimating N2O fluxes from cropland remains a challenging task because the related microbial processes (e.g., nitrification and denitrification) are controlled by complex interactions among climate, soil, plant and human activities. Existing approaches such as process-based (PB) models have well-known limitations due to insufficient representations of the processes or uncertainties of model parameters, and due to leverage recent advances in machine learning (ML) a new method is needed to unlock the "black box"to overcome its limitations such as low interpretability, out-of-sample failure and massive data demand. In this study, we developed a first-of-its-kind knowledge-guided machine learning model for agroecosystems (KGML-ag) by incorporating biogeophysical and chemical domain knowledge from an advanced PB model, ecosys, and tested it by comparing simulating daily N2O fluxes with real observed data from mesocosm experiments. The gated recurrent unit (GRU) was used as the basis to build the model structure. To optimize the model performance, we have investigated a range of ideas, including (1) using initial values of intermediate variables (IMVs) instead of time series as model input to reduce data demand; (2) building hierarchical structures to explicitly estimate IMVs for further N2O prediction; (3) using multi-task learning to balance the simultaneous training on multiple variables; and (4) pre-training with millions of synthetic data generated from ecosys and fine-tuning with mesocosm observations. Six other pure ML models were developed using the same mesocosm data to serve as the benchmark for the KGML-ag model. Results show that KGML-ag did an excellent job in reproducing the mesocosm N2O fluxes (overall r2Combining double low line0.81, and RMSECombining double low line3.6g€¯mgNm-2d-1 from cross validation). Importantly, KGML-ag always outperforms the PB model and ML models in predicting N2O fluxes, especially for complex temporal dynamics and emission peaks. Besides, KGML-ag goes beyond the pure ML models by providing more interpretable predictions as well as pinpointing desired new knowledge and data to further empower the current KGML-ag. We believe the KGML-ag development in this study will stimulate a new body of research on interpretable ML for biogeochemistry and other related geoscience processes.

Process-based land surface models are important tools for estimating global wetland methane (CH4) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH4 fluxes (FCH4) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH4 production). Our evaluation suggests the need to accurately replicate FCH4 variability, especially at short time scales, in future wetland CH4 model developments.

Agriculture contributes nearly a quarter of global greenhouse gas (GHG) emissions, which is motivating interest in adopting certain farming practices that have the potential to reduce GHG emissions or sequester carbon in soil. The related GHG emission (including N2O and CH4) and changes in soil carbon stock are defined here as “agricultural carbon outcomes”. Accurate quantification of agricultural carbon outcomes is the basis for achieving emission reductions for agriculture, but existing approaches for measuring carbon outcomes (including direct measurements, emission factors, and process-based modeling) fall short of achieving the required accuracy and scalability necessary to support credible, verifiable, and cost-effective measurement and improvement of these carbon outcomes. Here we propose a foundational and scalable framework to quantify field-level carbon outcomes for farmland, which is based on the holistic carbon balance of the agroecosystem: Agroecosystem Carbon Outcomes = Environment (E) × Management (M) × Crop (C). Following a comprehensive review of the scientific challenges associated with existing approaches, as well as their tradeoffs between cost and accuracy, we propose that the most viable path for the quantification of field-level carbon outcomes in agricultural land is through an effective integration of various approaches (e.g. diverse observations, sensor/in-situ data, and modeling), defined as the “System-of-Systems” solution. Such a “System-of-Systems” solution should simultaneously comprise the following components: (1) scalable collection of ground truth data and cross-scale sensing of environment variables (E), management practices (M), and crop conditions (C) at the local field level; (2) advanced modeling with necessary processes to support the quantification of carbon outcomes; (3) systematic Model-Data Fusion (MDF), i.e. robust and efficient methods to integrate sensing data and models at each local farmland level; (4) high computation efficiency and artificial intelligence (AI) to scale to millions of individual fields with low cost; and (5) robust and multi-tier validation systems and infrastructures to ensure solution fidelity and true scalability, i.e. the ability of a solution to perform robustly with accepted accuracy on all targeted fields. In this regard, we provide here the detailed scientific rationale, current progress, and future research and development (R&D) priorities to achieve different components of the “System-of-Systems” solution, thus accomplishing the Environment×Management×Crop framework to quantify field-level agricultural carbon outcomes.

Wetland CH₄ emissions are among the most uncertain components of the global CH₄ budget. The complex nature of wetland CH₄ processes makes it challenging to identify causal relationships for improving our understanding and predictability of CH₄ emissions. In this study, we used the flux measurements of CH₄ from eddy covariance towers (30 sites from 4 wetlands types: bog, fen, marsh, and wet tundra) to construct a causality-constrained machine learning (ML) framework to explain the regulative factors and to capture CH₄ emissions at sub-seasonal scale. We found that soil temperature is the dominant factor for CH₄ emissions in all studied wetland types. Ecosystem respiration (CO₂) and gross primary productivity exert controls at bog, fen, and marsh sites with lagged responses of days to weeks. Integrating these asynchronous environmental and biological causal relationships in predictive models significantly improved model performance. More importantly, modeled CH₄ emissions differed by up to a factor of 4 under a +1°C warming scenario when causality constraints were considered. These results highlight the significant role of causality in modeling wetland CH₄ emissions especially under future warming conditions, while traditional data-driven ML models may reproduce observations for the wrong reasons. Our proposed causality-guided model could benefit predictive modeling, large-scale upscaling, data gap-filling, and surrogate modeling of wetland CH₄ emissions within earth system land models.