Abstract. Atmospheric inverse modeling describes the process of estimating greenhouse gas fluxes or air pollution emissions at the Earth's surface using observations of these gases collected in the atmosphere. The launch of new satellites, the expansion of surface observation networks, and a desire for more detailed maps of surface fluxes have yielded numerous computational and statistical challenges for standard inverse modeling frameworks that were often originally designed with much smaller data sets in mind. In this article, we discuss computationally efficient methods for large-scale atmospheric inverse modeling and focus on addressing some of the main computational and practical challenges. We develop generalized hybrid projection methods, which are iterative methods for solving large-scale inverse problems, and specifically we focus on the case of estimating surface fluxes. These algorithms confer several advantages. They are efficient, in part because they converge quickly, they exploit efficient matrix–vector multiplications, and they do not require inversion of any matrices. These methods are also robust because they can accurately reconstruct surface fluxes, they are automatic since regularization or covariance matrix parameters and stopping criteria can be determined as part of the iterative algorithm, and they are flexible because they can be paired with many different types of atmospheric models. We demonstrate the benefits of generalized hybrid methods with a case study from NASA's Orbiting Carbon Observatory 2 (OCO-2) satellite. We then address the more challenging problem of solving the inverse model when the mean of the surface fluxes is not known a priori; we do so by reformulating the problem, thereby extending the applicability of hybrid projection methods to include hierarchical priors. We further show that by exploiting mathematical relations provided by the generalized hybrid method, we can efficiently calculate an approximate posterior variance, thereby providing uncertainty information.
more »
« less
Geostatistical inverse modeling with very large datasets: an example from the Orbiting Carbon Observatory 2 (OCO-2) satellite
Abstract. Geostatistical inverse modeling (GIM) has become a common approach to estimating greenhouse gas fluxes at the Earth's surface using atmospheric observations. GIMs are unique relative to other commonly used approaches because they do not require a single emissions inventory or a bottom–up model to serve as an initial guess of the fluxes. Instead, a modeler can incorporate a wide range of environmental, economic, and/or land use data to estimate the fluxes. Traditionally, GIMs have been paired with in situ observations that number in the thousands or tens of thousands. However, the number of available atmospheric greenhouse gas observations has been increasing enormously as the number of satellites, airborne measurement campaigns, and in situ monitoring stations continues to increase. This era of prolific greenhouse gas observations presents computational and statistical challenges for inverse modeling frameworks that have traditionally been paired with a limited number of in situ monitoring sites. In this article, we discuss the challenges of estimating greenhouse gas fluxes using large atmospheric datasets with a particular focus on GIMs. We subsequently discuss several strategies for estimating the fluxes and quantifying uncertainties, strategies that are adapted from hydrology, applied math, or other academic fields and are compatible with a wide variety of atmospheric models. We further evaluate the accuracy and computational burden of each strategy using a synthetic CO2 case study based upon NASA's Orbiting Carbon Observatory 2 (OCO-2) satellite. Specifically, we simultaneously estimate a full year of 3-hourly CO2 fluxes across North America in one case study – a total of 9.4×106 unknown fluxes using 9.9×104 synthetic observations. The strategies discussed here provide accurate estimates of CO2 fluxes that are comparable to fluxes calculated directly or analytically. We are also able to approximate posterior uncertainties in the fluxes, but these approximations are, typically, an over- or underestimate depending upon the strategy employed and the degree of approximation required to make the calculations manageable.
more »
« less
- Award ID(s):
- 1720398
- PAR ID:
- 10165258
- Date Published:
- Journal Name:
- Geoscientific Model Development
- Volume:
- 13
- Issue:
- 3
- ISSN:
- 1991-9603
- Page Range / eLocation ID:
- 1771 to 1785
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Eddy covariance measurement systems provide direct observation of the exchange of greenhouse gases between ecosystems and the atmosphere, but have only occasionally been intentionally applied to quantify the carbon dynamics associated with specific climate mitigation strategies. Natural climate solutions (NCS) harness the photosynthetic power of ecosystems to avoid emissions and remove atmospheric carbon dioxide (CO2), sequestering it in biological carbon pools. In this perspective, we aim to determine which kinds of NCS strategies are most suitable for ecosystem-scale flux measurements and how these measurements should be deployed for diverse NCS scales and goals. We find that ecosystem-scale flux measurements bring unique value when assessing NCS strategies characterized by inaccessible and hard-to-observe carbon pool changes, important non-CO2 greenhouse gas fluxes, the potential for biophysical impacts, or dynamic successional changes. We propose three deployment types for ecosystem-scale flux measurements at various NCS scales to constrain wide uncertainties and chart a workable path forward: “pilot”, “upscale”, and “monitor”. Together, the integration of ecosystem-scale flux measurements by the NCS community and the prioritization of NCS measurements by the flux community, have the potential to improve accounting in ways that capture the net impacts, unintended feedbacks, and on-the-ground specifics of a wide range of emerging NCS strategies.more » « less
-
Continuous greenhouse gas monitoring at sub-zero temperatures is needed for monitoring greenhouse gas emission in cold environments such as the Arctic tundra. This work reports a single-frequency electrochemical impedance sensing (SF-EIS) method for real-time continuous monitoring of carbon dioxide (CO2) at a wide range of temperatures (−15 to 40 °C) by using robust ionic liquid (IL) sensing materials and noninvasive, low-power, and low-cost impedance readout mechanisms since they cause minimal changes in the sensing interface, avoiding the baseline change for long-term continuous sensing. In addition, a miniaturized planar electrochemical sensor was fabricated that incorporates a hydrophobic 1-butyl-1-methylpyrrolidinium bis(trifluromethylsulfonyl)imide ([Bmpy][NTf2]) IL electrolyte and Pt black electrode materials. The high viscosity of the ILs facilitates the formation of thin, ordered, and concentrated layers of ionic charges, and the inverse relationship of IL viscosity with temperature makes them especially suited for impedance sensing at low temperatures. The unique low-temperature properties of ILs together with EIS transduction mechanisms are shown to be sensitive and selective for continuously monitoring CO2 at a −15 to 40 °C temperature range via impedance changes at a specifically selected frequency at the open circuit potential (OCP). Molecular dynamics simulations revealed insights into the structure and dynamics of the IL at varying temperatures in the presence of methane and CO2 and provided potential explanations for the observed sensing results. The miniaturized and flexible planar electrochemical sensor with the [Bmpy][NTf2] electrolyte was tested repeatedly at subzero temperatures over a 58-day period, during which good stability and repeatability were obtained. The CO2 impedance sensor was further tested for sensing CO2 from soil samples and shows promising results for their use in real-time monitoring of greenhouse gas emissions in cold temperatures such as permafrost soils.more » « less
-
null (Ed.)Herbivory can have strong impacts on greenhouse gas fluxes in high-latitude ecosystems. For example, in the Yukon-Kuskokwim (Y-K) Delta in western Alaska, migratory goose grazing affects the magnitude of soil carbon dioxide (CO2) and methane (CH4) fluxes. However, the underlying drivers of this relationship are unclear, as few studies systematically tease apart the processes by which herbivores influences soil biogeochemistry. To examine these mechanisms in detail, we conducted a laboratory incubation experiment to quantify changes in greenhouse gas fluxes in response to three parameters altered by herbivores in situ: temperature, soil moisture content, and nutrient inputs. These treatments were applied to soils collected in grazing lawns and nearby ungrazed habitat, allowing us to assess how variation in microbial community structure influenced observed responses. We found pronounced differences in both fungal and prokaryotic community composition between grazed and ungrazed areas. In the laboratory incubation experiment, CO2 and CH4 fluxes increased with temperature, soil moisture, and goose fecal addition, suggesting that grazing-related changes in the soil abiotic environment may enhance soil C losses. Yet, these abiotic drivers were insufficient to explain variation in fluxes between soils with and without prior grazing. Differences in trace gas fluxes between grazed and ungrazed areas may result both from herbivore-induced shifts in abiotic parameters and grazing-related alterations in microbial community structure. Our findings suggest that relationships among herbivores and soil microbial communities could mediate carbon-climate feedbacks in rapidly changing high-latitude ecosystems.more » « less
-
Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Emissions and atmospheric concentrations of CH4 continue to increase, maintaining CH4 as the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the factors explaining the well-observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in-situ and greenhouse gas observing satellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land-surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full datasets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this edition benefits from important progress in estimating inland freshwater emissions, with better accounting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double accounting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double accounting that still exists (average of 23 Tg CH4 yr-1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr-1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr-1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr-1 or ~65 % are attributed to direct anthropogenic sources in the fossil, agriculture and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr-1 or 63–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr-1 (range 9–40). Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr-1) larger global emissions (669 Tg CH4 yr-1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr-1 in Saunois et al. (2016, 2020), respectively), and for the first time uncertainty in bottom-up and top-down budgets overlap. The latitudinal distribution from atmospheric inversion-based emissions indicates a predominance of tropical and southern hemisphere emissions (~65 % of the global budget, <30° N) compared to mid (30° N–60° N, ~30 % of emissions) and high-northern latitudes (60° N–90° N, ~4 % of global emissions). This latitudinal distribution is similar in the bottom-up budget though the bottom-up budget estimates slightly larger contributions for the mid and high-northern latitudes, and slightly smaller contributions from the tropics and southern hemisphere than the inversions. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. We identify five major priorities for improving the CH4 budget: i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; ii) further development of process-based models for inland-water emissions; iii) intensification of CH4 observations at local (e.g., FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; iv) improvements of transport models and the representation of photochemical sinks in top-down inversions, and v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).more » « less
An official website of the United States government

