Benefits of satellite XCO2 and newly proposed atmospheric CO2 observation network over India in constraining regional CO2 fluxes

Top-down modeling estimates are among the most reliable information available on the CO₂ fluxes of the earth system. The inadequate coverage of CO₂ observing stations over the tropical regions adds a limitation to this estimate, especially when the satellite XCO₂ is strictly screened for cloud contamination, aerosol, dust, etc. In this study, we investigated the potential benefit of a global ground-based observing station network, 17 newly proposed stations over India, and global satellite XCO₂ in reducing the uncertainty of terrestrial biospheric fluxes of Tropical Asia-Eurasia in TransCom cyclo-stationary inversion. The data from selected 80 global ground-based CO₂ observation stations, together with two additional stations from India (i.e., Cape Rama and Sinhagad) and satellite XCO₂, helps to reduce the temperate Eurasian terrestrial flux uncertainty by 23.8%, 26.4%, and 36.2%, respectively. This further improved to 54.7% by adding the newly proposed stations over India into the inversion. By separating the Indian sub-continent from temperate Eurasia (as inspired by the heterogeneity in the terrestrial ecosystems, prevailing meteorological conditions, and the orography of this vast region), the inversion evinces the capacity of existing CO₂ observations to reduce the Indian terrestrial flux uncertainty by 20.5%. The largest benefit (70% reduction of annual mean uncertainty) for estimating Indian terrestrial fluxes could be achieved by combining these global observations with data from the newly proposed stations over India. The existing two stations from India suggest Temperate Eurasia as a mild source of CO₂ (0.33 ± 0.57 Pg C yr^-1), albeit with prominent anthropogenic influences visible in these two stations during the dry seasons. This implies that the proposed new stations should be cautiously placed to avoid such effects. The study also finds that the newly proposed stations over India also have an impact in constraining nearby oceanic CO₂ fluxes.

The Terrestrial Biosphere Model Farm

Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901-2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.

Atmospheric Carbon and Transport - America (ACT-America) Data Sets: Description, Management, and Delivery

The ACT-America project is a NASA Earth Venture Suborbital-2 mission designed to study the transport and fluxes of greenhouse gases. The open and freely available ACT-America data sets provide airborne in situ measurements of atmospheric carbon dioxide, methane, trace gases, aerosols, clouds, and meteorological properties, airborne remote sensing measurements of aerosol backscatter, atmospheric boundary layer height and columnar content of atmospheric carbon dioxide, tower-based measurements, and modeled atmospheric mole fractions and regional carbon fluxes of greenhouse gases over the Central and Eastern United States. We conducted 121 research flights during five campaigns in four seasons during 2016-2019 over three regions of the US (Mid-Atlantic, Midwest and South) using two NASA research aircraft (B-200 and C-130). We performed three flight patterns (fair weather, frontal crossings, and OCO-2 underflights) and collected more than 1,140 h of airborne measurements via level-leg flights in the atmospheric boundary layer, lower, and upper free troposphere and vertical profiles spanning these altitudes. We also merged various airborne in situ measurements onto a common standard sampling interval, which brings coherence to the data, creates geolocated data products, and makes it much easier for the users to perform holistic analysis of the ACT-America data products. Here, we report on detailed information of data sets collected, the workflow for data sets including storage and processing of the quality controlled and quality assured harmonized observations, and their archival and formatting for users. Finally, we provide some important information on the dissemination of data products including metadata and highlights of applications of ACT-America data sets.

Parameter Optimization of the 3PG Model Based on Sensitivity Analysis and a Bayesian Method

Sensitivity analysis and parameter optimization of stand models can improve their efficiency and accuracy, and increase their applicability. In this study, the sensitivity analysis, screening, and optimization of 63 model parameters of the Physiological Principles in Predicting Growth (3PG) model were performed by combining a sensitivity analysis method and the Markov chain Monte Carlo (MCMC) method of Bayesian posterior estimation theory. Additionally, a nine-year observational dataset of Chinese fir trees felled in the Shunchang Forest Farm, Nanping, was used to analyze, screen, and optimize the 63 model parameters of the 3PG model. The results showed the following: (1) The parameters that are most sensitive to stand stocking and diameter at breast height (DBH) are nWs(power in stem mass vs. diameter relationship), aWs(constant in stem mass vs. diameter relationship), alphaCx(maximum canopy quantum efficiency), k(extinction coefficient for PAR absorption by canopy), pRx(maximum fraction of NPP to roots), pRn(minimum fraction of NPP to roots), and CoeffCond(defines stomatal response to VPD); (2) MCMC can be used to optimize the parameters of the 3PG model, in which the posterior probability distributions of nWs, aWs, alphaCx, pRx, pRn, and CoeffCond conform to approximately normal or skewed distributions, and the peak value is prominent; and (3) compared with the accuracy before sensitivity analysis and a Bayesian method, the biomass simulation accuracy of the stand model was increased by 13.92%, and all indicators show that the accuracy of the improved model is superior. This method can be used to calibrate the parameters and analyze the uncertainty of multi-parameter complex stand growth models, which are important for the improvement of parameter estimation and simulation accuracy.

Quantifying the value of surveillance data for improving model predictions of lymphatic filariasis elimination

Background

Mathematical models are increasingly being used to evaluate strategies aiming to achieve the control or elimination of parasitic diseases. Recently, owing to growing realization that process-oriented models are useful for ecological forecasts only if the biological processes are well defined, attention has focused on data assimilation as a means to improve the predictive performance of these models.

Methodology and principal findings

We report on the development of an analytical framework to quantify the relative values of various longitudinal infection surveillance data collected in field sites undergoing mass drug administrations (MDAs) for calibrating three lymphatic filariasis (LF) models (EPIFIL, LYMFASIM, and TRANSFIL), and for improving their predictions of the required durations of drug interventions to achieve parasite elimination in endemic populations. The relative information contribution of site-specific data collected at the time points proposed by the WHO monitoring framework was evaluated using model-data updating procedures, and via calculations of the Shannon information index and weighted variances from the probability distributions of the estimated timelines to parasite extinction made by each model. Results show that data-informed models provided more precise forecasts of elimination timelines in each site compared to model-only simulations. Data streams that included year 5 post-MDA microfilariae (mf) survey data, however, reduced each model’s uncertainty most compared to data streams containing only baseline and/or post-MDA 3 or longer-term mf survey data irrespective of MDA coverage, suggesting that data up to this monitoring point may be optimal for informing the present LF models. We show that the improvements observed in the predictive performance of the best data-informed models may be a function of temporal changes in inter-parameter interactions. Such best data-informed models may also produce more accurate predictions of the durations of drug interventions required to achieve parasite elimination.

Significance

Knowledge of relative information contributions of model only versus data-informed models is valuable for improving the usefulness of LF model predictions in management decision making, learning system dynamics, and for supporting the design of parasite monitoring programmes. The present results further pinpoint the crucial need for longitudinal infection surveillance data for enhancing the precision and accuracy of model predictions of the intervention durations required to achieve parasite elimination in an endemic location.

Although parasite transmission models offer powerful tools for predicting the impacts of interventions, there is growing realization that these models can be useful for this purpose only if their governing biological processes are well defined. Recently, model-data assimilation has been applied to address this problem and improve the performance of process-oriented models for ecological forecasting. Here, we developed an analytical framework that allowed the sequential coupling of the three existing lymphatic filariasis (LF) models with longitudinal infection monitoring data collected in field sites undergoing mass drug administrations (MDAs) to examine the relative value of such data for parameterizing these models and for improving their predictions of the required durations of drug interventions to break parasite transmission. We found that data-informed models provided more precise and reliable forecasts of elimination timelines in the study sites compared to model-only predictions, and that data collected up to 5 years post-MDA reduced each model’s predictive uncertainty most. We also found that this improved performance may be intriguingly related to temporal changes in system dynamics. Our results underscore the significance of sequential model-data fusion for enhancing the understanding of LF transmission dynamics, design of surveillance, and generation of reliable model predictions for management decision making.

Rate my data: quantifying the value of ecological data for the development of models of the terrestrial carbon cycle

Primarily driven by concern about rising levels of atmospheric CO2, ecologists and earth system scientists are collecting vast amounts of data related to the carbon cycle. These measurements are generally time consuming and expensive to make, and, unfortunately, we live in an era where research funding is increasingly hard to come by. Thus, important questions are: "Which data streams provide the most valuable information?" and "How much data do we need?" These questions are relevant not only for model developers, who need observational data to improve, constrain, and test their models, but also for experimentalists and those designing ecological observation networks. Here we address these questions using a model-data fusion approach. We constrain a process-oriented, forest ecosystem C cycle model with 17 different data streams from the Harvard Forest (Massachusetts, USA). We iteratively rank each data source according to its contribution to reducing model uncertainty. Results show the importance of some measurements commonly unavailable to carbon-cycle modelers, such as estimates of turnover times from different carbon pools. Surprisingly, many data sources are relatively redundant in the presence of others and do not lead to a significant improvement in model performance. A few select data sources lead to the largest reduction in parameter-based model uncertainty. Projections of future carbon cycling were poorly constrained when only hourly net-ecosystem-exchange measurements were used to inform the model. They were well constrained, however, with only 5 of the 17 data streams, even though many individual parameters are not constrained. The approach taken here should stimulate further cooperation between modelers and measurement teams and may be useful in the context of setting research priorities and allocating research funds.