scholarly journals Order of magnitude wall time improvement of variational methane inversions by physical parallelization: a demonstration using TM5-4DVAR

2021 ◽  
Author(s):  
Sudhanshu Pandey ◽  
Sander Houweling ◽  
Arjo Segers

Abstract. Atmospheric inversions are used to constrain the emissions of trace gases from atmospheric mole fraction measurements. The variational (4DVAR) inversion approach allows optimization of the emissions at a much higher temporal and spatial resolution than the ensemble or analytical approaches but provides limited opportunities for scalable parallelization as the optimization is performed iteratively. Multidecadal variational inversions are used to optimally extract information from the long measurement records of long-lived atmospheric trace gases like carbon dioxide and methane. However, the wall clock time needed—up to months— complicates these multidecadal inversions. The physical parallelization method introduced by Chevallier (2013) addresses this problem for CO2 inversions by splitting the time period of the chemical transport model into blocks that are run in parallel. Here we present a new implementation of the physical parallelization for variational inversion (PPVI) approach that is suitable for methane inversions as it accounts for methane’s atmospheric lifetime. The performance of PPVI is tested in an 11-year inversion using a TM5-4DVAR inversion setup that assimilates surface observations to optimize methane emissions at grid-scale. We find that the PPVI inversion approach improves the wall clock time performance by a factor of 5 and shows excellent agreement with the posterior emissions of a full serial inversion with identical configuration (global mean emissions difference = 0.06 % with an interannual variation correlation R = 99 %; regional mean emission difference < 5 % and interannual variation R > 0.95). The wall clock time improvement using the PPVI method increases with the size of the inversion period. The PPVI approach is planned to be used in future releases of the CAMS (Copernicus Atmosphere Monitoring Service) multidecadal methane reanalysis.

2011 ◽  
Vol 11 (17) ◽  
pp. 9253-9269 ◽  
Author(s):  
J. Angelbratt ◽  
J. Mellqvist ◽  
D. Simpson ◽  
J. E. Jonson ◽  
T. Blumenstock ◽  
...  

Abstract. Trends in the CO andC2H6 partial columns ~0–15 km) have been estimated from four European ground-based solar FTIR (Fourier Transform InfraRed) stations for the 1996–2006 time period. The CO trends from the four stations Jungfraujoch, Zugspitze, Harestua and Kiruna have been estimated to −0.45 ± 0.16% yr−1, −1.00 ± 0.24% yr−1, −0.62 ± 0.19 % yr−1 and −0.61 ± 0.16% yr−1, respectively. The corresponding trends for C2H6 are −1.51 ± 0.23% yr−1, −2.11 ± 0.30% yr−1, −1.09 ± 0.25% yr−1 and −1.14 ± 0.18% yr−1. All trends are presented with their 2-σ confidence intervals. To find possible reasons for the CO trends, the global-scale EMEP MSC-W chemical transport model has been used in a series of sensitivity scenarios. It is shown that the trends are consistent with the combination of a 20% decrease in the anthropogenic CO emissions seen in Europe and North America during the 1996–2006 period and a 20% increase in the anthropogenic CO emissions in East Asia, during the same time period. The possible impacts of CH4 and biogenic volatile organic compounds (BVOCs) are also considered. The European and global-scale EMEP models have been evaluated against the measured CO and C2H6 partial columns from Jungfraujoch, Zugspitze, Bremen, Harestua, Kiruna and Ny-Ålesund. The European model reproduces, on average the measurements at the different sites fairly well and within 10–22% deviation for CO and 14–31% deviation for C2H6. Their seasonal amplitude is captured within 6–35% and 9–124% for CO and C2H6, respectively. However, 61–98% of the CO and C2H6 partial columns in the European model are shown to arise from the boundary conditions, making the global-scale model a more suitable alternative when modeling these two species. In the evaluation of the global model the average partial columns for 2006 are shown to be within 1–9% and 37–50% of the measurements for CO and C2H6, respectively. The global model sensitivity for assumptions made in this paper is also analyzed.


2016 ◽  
Vol 9 (8) ◽  
pp. 2741-2754 ◽  
Author(s):  
Elham Baranizadeh ◽  
Benjamin N. Murphy ◽  
Jan Julin ◽  
Saeed Falahat ◽  
Carly L. Reddington ◽  
...  

Abstract. The particle formation scheme within PMCAMx-UF, a three-dimensional chemical transport model, was updated with particle formation rates for the ternary H2SO4–NH3–H2O pathway simulated by the Atmospheric Cluster Dynamics Code (ACDC) using quantum chemical input data. The model was applied over Europe for May 2008, during which the EUCAARI-LONGREX (European Aerosol Cloud Climate and Air Quality Interactions–Long-Range Experiment) campaign was carried out, providing aircraft vertical profiles of aerosol number concentrations. The updated model reproduces the observed number concentrations of particles larger than 4 nm within 1 order of magnitude throughout the atmospheric column. This agreement is encouraging considering the fact that no semi-empirical fitting was needed to obtain realistic particle formation rates. The cloud adjustment scheme for modifying the photolysis rate profiles within PMCAMx-UF was also updated with the TUV (Tropospheric Ultraviolet and Visible) radiative-transfer model. Results show that, although the effect of the new cloud adjustment scheme on total number concentrations is small, enhanced new-particle formation is predicted near cloudy regions. This is due to the enhanced radiation above and in the vicinity of the clouds, which in turn leads to higher production of sulfuric acid. The sensitivity of the results to including emissions from natural sources is also discussed.


2010 ◽  
Vol 10 (1) ◽  
pp. 977-1004
Author(s):  
C. Paton-Walsh ◽  
L. K. Emmons ◽  
S. R. Wilson

Abstract. In this paper we describe a new method for estimating trace gas emissions from large vegetation fires using measurements of aerosol optical depth from the MODIS instruments onboard NASA's Terra and Aqua satellites, combined with the atmospheric chemical transport model MOZART. The model allows for an estimate of double counting of enhanced levels of aerosol optical depth in consecutive satellite overpasses. Using this method we infer an estimated total emission of 10±3 Tg of carbon monoxide from the Canberra fires of 2003. Emissions estimates for several other trace gases are also given. An assessment of the uncertainties in the new method is made and we show that our estimate agrees (within expected uncertainties) with estimates made using current conventional methods of multiplying together factors for the area burned, fuel load, the combustion efficiency and the emission factor for carbon monoxide. The new method for estimating emissions from large vegetation fires described in this paper has some significant uncertainties, but these are mainly quantifiable and largely independent of the uncertainties inherent in conventional techniques. Thus we conclude that the new method is a useful additional tool for characterising emissions from vegetation fires.


2019 ◽  
Author(s):  
Erin Evoy ◽  
Adrian M. Maclean ◽  
Grazia Rovelli ◽  
Ying Li ◽  
Alexandra P. Tsimpidi ◽  
...  

Abstract. Information on the rate of diffusion of organic molecules within secondary organic aerosol (SOA) is needed to accurately predict the effects of SOA on climate and air quality. Often, researchers have predicted diffusion rates of organic molecules within SOA using measurements of viscosity and the Stokes-Einstein relation (D ∝ 1/η where D is the diffusion coefficient and η is viscosity). However, the accuracy of this relation for predicting diffusion in SOA remains uncertain. We measured diffusion coefficients over eight orders in magnitude in proxies of SOA including citric acid, sorbitol, and a sucrose-citric acid mixture. These results were combined with literature data to evaluate the Stokes-Einstein relation for predicting diffusion of organic molecules in SOA. Although almost all the data agrees with the Stokes-Einstein relation within a factor of ten, a fractional Stokes-Einstein relation (D ∝ C/ηt) with t = 0.93 and C = 1.66 is a better model for predicting diffusion of organic molecules in the SOA proxies studied. In addition, based on the output from a chemical transport model, the Stokes-Einstein relation can over predict mixing times of organic molecules within SOA by as much as one order of magnitude at an altitude ~ 3 km, compared to the fractional Stokes-Einstein relation with t = 0.93 and C = 1.66. These differences can be important for predicting growth, evaporation, and reaction rates of SOA in the middle and upper part of the troposphere. These results also have implications for other areas where diffusion of organic molecules within organic-water matrices is important.


2012 ◽  
Vol 12 (1) ◽  
pp. 371-380 ◽  
Author(s):  
R. Hossaini ◽  
M. P. Chipperfield ◽  
W. Feng ◽  
T. J. Breider ◽  
E. Atlas ◽  
...  

Abstract. We have used a global three-dimensional chemical transport model to quantify the impact of the very short-lived substances (VSLS) CHBr3, CH2Br2, CHBr2Cl, CHBrCl2, CH2BrCl and C2H5-Br on the bromine budget of the stratosphere. Atmospheric observations of these gases allow constraints on surface mixing ratios that, when incorporated into our model, contribute ~4.9–5.2 parts per trillion (ppt) of inorganic bromine (Bry) to the stratosphere. Of this total, ~76 % comes from naturally-emitted CHBr3 and CH2Br2. The remaining species individually contribute modest amounts. However, their accumulated total accounts for up to ~1.2 ppt of the supply and thus should not be ignored. We have compared modelled tropical profiles of a range of VSLS with observations from the recent 2009 NSF HIPPO-1 aircraft campaign. Modelled profiles agree reasonably well with observations from the surface to the lower tropical tropopause layer. We have also considered the poorly studied anthropogenic VSLS, C2H5Br, CH2BrCH2Br, n-C3H7Br and i-C3H7Br. We find the local atmospheric lifetime of these species in the tropical tropopause layer are ~183, 603, 39 and 49 days, respectively. These species, particularly C2H5Br and CH2BrCH2Br, would thus be important carriers of bromine to the stratosphere if emissions were to increase substantially. Our model shows ~70–73 % and ~80–85 % of bromine from these species in the tropical boundary layer can reach the lower stratosphere.


2010 ◽  
Vol 10 (1) ◽  
pp. 925-945 ◽  
Author(s):  
V. Naik ◽  
A. M. Fiore ◽  
L. W. Horowitz ◽  
H. B. Singh ◽  
C. Wiedinmyer ◽  
...  

Abstract. Energy security and climate change concerns have led to the promotion of biomass-derived ethanol, an oxygenated volatile organic compound (OVOC), as a substitute for fossil fuels. Although ethanol is ubiquitous in the troposphere, our knowledge of its current atmospheric budget and distribution is limited. Here, for the first time we use a global chemical transport model in conjunction with atmospheric observations to place constraints on the ethanol budget, noting that additional measurements of ethanol (and its precursors) are still needed to enhance confidence in our estimated budget. Global sources of ethanol in the model include 5.0 Tg yr−1 from industrial sources and biofuels, 9.2 Tg yr−1 from terrestrial plants, ~0.5 Tg yr−1 from biomass burning, and 0.05 Tg yr−1 from atmospheric reactions of the ethyl peroxide radical (C2H5O2) with itself and with the methyl peroxide radical (CH3O2). The resulting atmospheric lifetime of ethanol in the model is 2.8 days. Gas-phase oxidation by hydroxyl radical (OH) is the primary global sink of ethanol in the model (65%), followed by dry deposition to land (25%), and wet deposition (10%). Over continental areas, ethanol concentrations predominantly reflect direct anthropogenic and biogenic emission sources. Uncertainty in the biogenic ethanol emissions estimated at a factor of three may contribute to the 50% model underestimate of observations in the North American boundary layer. Furthermore, current levels of ethanol measured in remote atmospheres are an order of magnitude larger than those explained by surface sources or by in-situ atmospheric production from observed precursor hydrocarbons in the model, suggesting a major gap in understanding. Stronger constraints on the budget and distribution of ethanol and other VOCs are a critical step towards assessing the impacts of increasing use of ethanol as a fuel.


2016 ◽  
Author(s):  
Whitney Bader ◽  
Benoît Bovy ◽  
Stephanie Conway ◽  
Kimberly Strong ◽  
Dan Smale ◽  
...  

Abstract. Changes of atmospheric methane (CH4) since 2005 have been evaluated using Fourier Transform Infrared (FTIR) solar observations performed at ten ground-based sites, all members of the Network for Detection of Atmospheric Composition Change (NDACC). From this, we find an increase of atmospheric methane total columns that amounts to 0.31 ± 0.03 % year−1 (2-sigma level of uncertainty) for the 2005–2014 period. Comparisons with in situ methane measurements at both local and global scales show good agreement. We used the GEOS-Chem Chemical Transport Model tagged simulation that accounts for the contribution of each emission source and one sink in the total methane, simulated over the 2005–2012 time period and based on emissions inventories and transport. After regridding according to NDACC vertical layering using a conservative regridding scheme and smoothing by convolving with respective FTIR seasonal averaging kernels, the GEOS-Chem simulation shows an increase of atmospheric methane of 0.35 ± 0.03 % year−1 between 2005 and 2012, which is in agreement with NDACC measurements over the same time period (0.30 ± 0.04 % year−1, averaged over ten stations). Analysis of the GEOS-Chem tagged simulation allows us to quantify the contribution of each tracer to the global methane change since 2005. We find that natural sources such as wetlands and biomass burning contribute to the inter-annual variability of methane. However, anthropogenic emissions such as coal mining, and gas and oil transport and exploration, which are mainly emitted in the Northern Hemisphere and act as secondary contributors to the global budget of methane, have played a major role in the increase of atmospheric methane observed since 2005. Based on the GEOS-Chem tagged simulation, we discuss possible cause(s) for the increase of methane since 2005, which is still unexplained.


2016 ◽  
Author(s):  
Elham Baranizadeh ◽  
Benjamin N. Murphy ◽  
Jan Julin ◽  
Saeed Falahat ◽  
Carly L. Reddington ◽  
...  

Abstract. The particle formation scheme within PMCAMx-UF, a three dimensional chemical transport model, was updated with particle formation rates for the ternary H2SO4-NH3-H2O pathway simulated by the Atmospheric Cluster Dynamics Code (ACDC) using quantum chemical input data. The model was applied over Europe for May 2008, during which the EUCAARI-LONGREX campaign was carried out providing aircraft vertical profiles of aerosol number concentrations. The updated model reproduces the observed number concentrations of particles larger than 4 nm within one order of magnitude throughout the atmospheric column. This reasonable agreement is very encouraging considering the fact that no semi-empirical fitting was needed to obtain realistic particle formation rates. The cloud adjustment scheme for modifying the photolysis rate profiles within PMCAMx-UF was also updated with the TUV (Tropospheric Ultraviolet and Visible) radiative transfer model. Results show that although the effect of the new cloud adjustment scheme on total number concentrations is small, enhanced new particle formation is predicted near cloudy regions. This is due to the enhanced radiation above and in the vicinity of the clouds, which in turn leads to higher production of sulfuric acid. The sensitivity of the results to including emissions from natural sources is also discussed.


2020 ◽  
Author(s):  
Lei Zhu ◽  
Gonzalo González Abad ◽  
Caroline R. Nowlan ◽  
Christopher Chan Miller ◽  
Kelly Chance ◽  
...  

Abstract. Formaldehyde (HCHO) has been measured from space for more than two decades. Owing to its short atmospheric lifetime, satellite HCHO data are used widely as a proxy of volatile organic compounds (VOCs; please refer to Appendix A for abbreviations and acronyms), providing constraints on underlying emissions and chemistry. However, satellite HCHO products from different satellite sensors using different algorithms have received little validation so far. The accuracy and consistency of HCHO retrievals remain largely unclear. Here we develop a global validation platform for satellite HCHO retrievals using in situ observations from 12 aircraft campaigns with a chemical transport model (GEOS-Chem) as the intercomparison method. Application to the NASA operational OMI HCHO product indicates slight biases (−30.9 % to +16.0 %) under high-HCHO conditions partially caused by a priori shape factors used in the retrievals, while high biases (+113.9 % to +194.6 %) under low-HCHO conditions due mainly to slant column fitting and radiance reference sector correction. By providing quick assessment to systematic biases in satellite products over large domains, the platform facilitates, in an iterative process, optimization of retrieval settings and the minimization of retrieval biases. It is also complementary to localized validation efforts based on ground observations and aircraft spirals.


Sign in / Sign up

Export Citation Format

Share Document