Abstract. We use 2010–2015 observations of atmospheric methane columns from the GOSAT satellite instrument in a global inverse analysis to improve estimates of methane emissions and their trends over the period, as well as the global concentration of tropospheric OH (the hydroxyl radical, methane's main sink) and its trend. Our inversion solves the Bayesian optimization problem analytically including closed-form characterization of errors. This allows us to (1) quantify the information content from the inversion towards optimizing methane emissions and its trends, (2) diagnose error correlations between constraints on emissions and OH concentrations, and (3) generate a large ensemble of solutions testing different assumptions in the inversion. We show how the analytical approach can be used, even when prior error standard deviation distributions are lognormal. Inversion results show large overestimates of Chinese coal emissions and Middle East oil and gas emissions in the EDGAR v4.3.2 inventory but little error in the United States where we use a new gridded version of the EPA national greenhouse gas inventory as prior estimate. Oil and gas emissions in the EDGAR v4.3.2 inventory show large differences with national totals reported to the United Nations Framework Convention on Climate Change (UNFCCC), and our inversion is generally more consistent with the UNFCCC data. The observed 2010–2015 growth in atmospheric methane is attributed mostly to an increase in emissions from India, China, and areas with large tropical wetlands. The contribution from OH trends is small in comparison. We find that the inversion provides strong independent constraints on global methane emissions (546 Tg a−1) and global mean OH concentrations (atmospheric methane lifetime against oxidation by tropospheric OH of 10.8±0.4 years), indicating that satellite observations of atmospheric methane could provide a proxy for OH concentrations in the future.

Abstract. Methane is a greenhouse gas emitted by a range of natural and anthropogenic sources. Atmospheric methane has been measured continuously from space since 2003, and new instruments are planned for launch in the near future that will greatly expand the capabilities of space-based observations. We review the value of current, future, and proposed satellite observations to better quantify and understand methane emissions through inverse analyses, from the global scale down to the scale of point sources and in combination with suborbital (surface and aircraft) data. Current global observations from Greenhouse Gases Observing Satellite (GOSAT) are of high quality but have sparse spatial coverage. They can quantify methane emissions on a regional scale (100-1000 km) through multiyear averaging. The Tropospheric Monitoring Instrument (TROPOMI), to be launched in 2017, is expected to quantify daily emissions on the regional scale and will also effectively detect large point sources. A different observing strategy by GHGSat (launched in June 2016) is to target limited viewing domains with very fine pixel resolution in order to detect a wide range of methane point sources. Geostationary observation of methane, still in the proposal stage, will have the unique capability of mapping source regions with high resolution, detecting transient "super-emitter" point sources and resolving diurnal variation of emissions from sources such as wetlands and manure. Exploiting these rapidly expanding satellite measurement capabilities to quantify methane emissions requires a parallel effort to construct high-quality spatially and sectorally resolved emission inventories. Partnership between top-down inverse analyses of atmospheric data and bottom-up construction of emission inventories is crucial to better understanding methane emission processes and subsequently informing climate policy.

The global atmospheric sulfur budget and its emission dependence have been investigated using the coupled aerosol-chemistry-climate model SOCOL-AER. The aerosol module comprises gaseous and aqueous sulfur chemistry and comprehensive microphysics. The particle distribution is resolved by 40 size bins spanning radii from 0.39 nm to 3.2 μm, including size-dependent particle composition. Aerosol radiative properties required by the climate model are calculated online from the aerosol module. The model successfully reproduces main features of stratospheric aerosols under nonvolcanic conditions, including aerosol extinctions compared to Stratospheric Aerosol and Gas Experiment II (SAGE II) and Halogen Occultation Experiment, and size distributions compared to in situ measurements. The calculated stratospheric aerosol burden is 109 Gg of sulfur, matching the SAGE II-based estimate (112 Gg). In terms of fluxes through the tropopause, the stratospheric aerosol layer is due to about 43% primary tropospheric aerosol, 28% SO 2 , 23% carbonyl sulfide (OCS), 4% H 2 S, and 2% dimethyl sulfide (DMS). Turning off emissions of the short-lived species SO 2 , H 2 S, and DMS shows that OCS alone still establishes about 56% of the original stratospheric aerosol burden. Further sensitivity simulations reveal that anticipated increases in anthropogenic SO 2 emissions in China and India have a larger influence on stratospheric aerosols than the same increase in Western Europe or the U.S., due to deep convection in the western Pacific region. However, even a doubling of Chinese and Indian emissions is predicted to increase the stratospheric background aerosol burden only by 9%. In contrast, small to moderate volcanic eruptions, such as that of Nabro in 2011, may easily double the stratospheric aerosol loading.

Abstract. Individual countries report national emissions of methane, a potent greenhouse gas, in accordance with the United Nations Framework Convention on Climate Change (UNFCCC). We present a global inventory of methane emissions from oil, gas, and coal exploitation that spatially allocates the national emissions reported to the UNFCCC (Scarpelli et al., 2019). Our inventory is at 0.1∘×0.1∘ resolution and resolves the subsectors of oil and gas exploitation, from upstream to downstream, and the different emission processes (leakage, venting, flaring). Global emissions for 2016 are 41.5 Tg a−1 for oil, 24.4 Tg a−1 for gas, and 31.3 Tg a−1 for coal. An array of databases is used to spatially allocate national emissions to infrastructure, including wells, pipelines, oil refineries, gas processing plants, gas compressor stations, gas storage facilities, and coal mines. Gridded error estimates are provided in normal and lognormal forms based on emission factor uncertainties from the IPCC. Our inventory shows large differences with the EDGAR v4.3.2 global gridded inventory both at the national scale and in finer-scale spatial allocation. It shows good agreement with the gridded version of the United Kingdom's National Atmospheric Emissions Inventory (NAEI). There are significant errors on the 0.1∘×0.1∘ grid associated with the location and magnitude of large point sources, but these are smoothed out when averaging the inventory over a coarser grid. Use of our inventory as prior estimate in inverse analyses of atmospheric methane observations allows investigation of individual subsector contributions and can serve policy needs by evaluating the national emissions totals reported to the UNFCCC. Gridded data sets can be accessed at https://doi.org/10.7910/DVN/HH4EUM (Scarpelli et al., 2019).

Abstract. In terms of atmospheric impact, the volcanic eruption of Mt. Pinatubo (1991) is the best characterized large eruption on record. We investigate here the model-derived stratospheric warming following the Pinatubo eruption as derived from SAGE II extinction data including recent improvements in the processing algorithm. This method, termed SAGE_4λ, makes use of the four wavelengths (385, 452, 525 and 1024 nm) of the SAGE II data when available, and uses a data-filling procedure in the opacity-induced "gap" regions. Using SAGE_4λ, we derived aerosol size distributions that properly reproduce extinction coefficients also at much longer wavelengths. This provides a good basis for calculating the absorption of terrestrial infrared radiation and the resulting stratospheric heating. However, we also show that the use of this data set in a global chemistry–climate model (CCM) still leads to stronger aerosol-induced stratospheric heating than observed, with temperatures partly even higher than the already too high values found by many models in recent general circulation model (GCM) and CCM intercomparisons. This suggests that the overestimation of the stratospheric warming after the Pinatubo eruption may not be ascribed to an insufficient observational database but instead to using outdated data sets, to deficiencies in the implementation of the forcing data, or to radiative or dynamical model artifacts. Conversely, the SAGE_4λ approach reduces the infrared absorption in the tropical tropopause region, resulting in a significantly better agreement with the post-volcanic temperature record at these altitudes.