Validation of the Swiss methane emission inventory by atmospheric observations and inverse modelling

Atmospheric inverse modelling has the potential to provide observation-based estimates of greenhouse gas emissions at the country scale, thereby allowing for an independent validation of national emission inventories. Here, we present a regional-scale inverse modelling study to quantify the emissions of methane (CH₄) from Switzerland, making use of the newly established CarboCount-CH measurement network and a high-resolution Lagrangian transport model. In our reference inversion, prior emissions were taken from the "bottom-up" Swiss Greenhouse Gas Inventory (SGHGI) as published by the Swiss Federal Office for the Environment in 2014 for the year 2012. Overall we estimate national CH₄ emissions to be 196 ± 18 Gg yr⁻¹ for the year 2013 (1σ uncertainty). This result is in close agreement with the recently revised SGHGI estimate of 206 ± 33 Gg yr⁻¹ as reported in 2015 for the year 2012. Results from sensitivity inversions using alternative prior emissions, uncertainty covariance settings, large-scale background mole fractions, two different inverse algorithms (Bayesian and extended Kalman filter), and two different transport models confirm the robustness and independent character of our estimate. According to the latest SGHGI estimate the main CH₄ source categories in Switzerland are agriculture (78 %), waste handling (15 %) and natural gas distribution and combustion (6 %). The spatial distribution and seasonal variability of our posterior emissions suggest an overestimation of agricultural CH₄ emissions by 10 to 20 % in the most recent SGHGI, which is likely due to an overestimation of emissions from manure handling. Urban areas do not appear as emission hotspots in our posterior results, suggesting that leakages from natural gas distribution are only a minor source of CH₄ in Switzerland. This is consistent with rather low emissions of 8.4 Gg yr⁻¹ reported by the SGHGI but inconsistent with the much higher value of 32 Gg yr⁻¹ implied by the EDGARv4.2 inventory for this sector. Increased CH₄ emissions (up to 30 % compared to the prior) were deduced for the north-eastern parts of Switzerland. This feature was common to most sensitivity inversions, which is a strong indicator that it is a real feature and not an artefact of the transport model and the inversion system. However, it was not possible to assign an unambiguous source process to the region. The observations of the CarboCount-CH network provided invaluable and independent information for the validation of the national bottom-up inventory. Similar systems need to be sustained to provide independent monitoring of future climate agreements.

Sensitivity of atmospheric CO2 and climate to explosive volcanic eruptions

Impacts of low-latitude, explosive volcanic eruptions on climate and the carbon cycle are quantified by forcing a comprehensive, fully coupled carbon cycle-climate model with pulse-like stratospheric aerosol optical depth changes. The model represents the radiative and dynamical response of the climate system to volcanic eruptions and simulates a decrease of global and regional atmospheric surface temperature, regionally distinct changes in precipitation, a positive phase of the North Atlantic Oscillation, and a decrease in atmospheric CO2 after volcanic eruptions. The volcanic-induced cooling reduces overturning rates in tropical soils, which dominates over reduced litter input due to soil moisture decrease, resulting in higher land carbon inventories for several decades. The perturbation in the ocean carbon inventory changes sign from an initial weak carbon sink to a carbon source. Positive carbon and negative temperature anomalies in subsurface waters last up to several decades. The multi-decadal decrease in atmospheric CO2 yields a small additional radiative forcing that amplifies the cooling and perturbs the Earth System on longer time scales than the atmospheric residence time of volcanic aerosols. In addition, century-scale global warming simulations with and without volcanic eruptions over the historical period show that the ocean integrates volcanic radiative cooling and responds for different physical and biogeochemical parameters such as steric sea level or dissolved oxygen. Results from a suite of sensitivity simulations with different magnitudes of stratospheric aerosol optical depth changes and from global warming simulations show that the carbon cycle-climate sensitivity γ, expressed as change in atmospheric CO2 per unit change in global mean surface temperature, depends on the magnitude and temporal evolution of the perturbation, and time scale of interest. On decadal time scales, modeled γ is several times larger for a Pinatubo-like eruption than for the industrial period and for a high emission, 21st century scenario.

First middle-atmospheric zonal wind profile measurements with a new ground-based microwave Doppler-spectro-radiometer

We report on the wind radiometer WIRA, a new ground-based microwave Doppler-spectro-radiometer specifically designed for the measurement of middle-atmospheric horizontal wind by observing ozone emission spectra at 142.17504 GHz. Currently, wind speeds in five levels between 30 and 79 km can be retrieved which makes WIRA the first instrument able to continuously measure horizontal wind in this altitude range. For an integration time of one day the measurement error on each level lies at around 25 m s−1. With a planned upgrade this value is expected to be reduced by a factor of 2 in the near future. On the altitude levels where our measurement can be compared to wind data from the European Centre for Medium-Range Weather Forecasts (ECMWF) very good agreement in the long-term statistics as well as in short time structures with a duration of a few days has been found. WIRA uses a passive double sideband heterodyne receiver together with a digital Fourier transform spectrometer for the data acquisition. A big advantage of the radiometric approach is that such instruments can also operate under adverse weather conditions and thus provide a continuous time series for the given location. The optics enables the instrument to scan a wide range of azimuth angles including the directions east, west, north, and south for zonal and meridional wind measurements. The design of the radiometer is fairly compact and its calibration does not rely on liquid nitrogen which makes it transportable and suitable for campaign use. WIRA is conceived in a way that it can be operated remotely and does hardly require any maintenance. In the present paper, a description of the instrument is given, and the techniques used for the wind retrieval based on the determination of the Doppler shift of the measured atmospheric ozone emission spectra are outlined. Their reliability was tested using Monte Carlo simulations. Finally, a time series of 11 months of zonal wind measurements over Bern (46°57′ N, 7°26′ E) is presented and compared to ECMWF wind data.

Natural and anthropogenic methane fluxes in Eurasia: a mesoscale quantification by generalized atmospheric inversion

Eight surface observation sites providing quasi-continuous measurements of atmospheric methane mixingratios have been operated since the mid-2000’s in Siberia. For the first time in a single work, we assimilate 1 year of these in situ observations in an atmospheric inversion. Our objective is to quantify methane surface fluxes from anthropogenic and wetland sources at the mesoscale in the Siberian lowlands for the year 2010. To do so, we first inquire about the way the inversion uses the observations and the way the fluxes are constrained by the observation sites. As atmospheric inver- sions at the mesoscale suffer from mis-quantified sources of uncertainties, we follow recent innovations in inversion techniques and use a new inversion approach which quantifies the uncertainties more objectively than the previous inversion systems. We find that, due to errors in the representation of the atmospheric transport and redundant pieces of information, only one observation every few days is found valuable by the inversion. The remaining high-resolution quasi-continuous signal is representative of very local emission patterns difficult to analyse with a mesoscale system. An analysis of the use of information by the inversion also reveals that the observation sites constrain methane emissions within a radius of 500 km. More observation sites than the ones currently in operation are then necessary to constrain the whole Siberian lowlands. Still, the fluxes within the constrained areas are quantified with objectified uncertainties. Finally, the tolerance intervals for posterior methane fluxes are of roughly 20 % (resp. 50 %) of the fluxes for anthropogenic (resp. wetland) sources. About 50–70 % of Siberian lowlands emissions are constrained by the inversion on average on an annual basis. Extrapolating the figures on the constrained areas to the whole Siberian lowlands, we find a regional methane budget of 5–28 TgCH4 for the year 2010, i.e. 1–5 % of the global methane emissions. As very few in situ observations are available in the region of interest, observations of methane total columns from the Greenhouse Gas Observing SATellite (GOSAT) are tentatively used for the evaluation of the inversion results, but they exhibit only a marginal signal from the fluxes within the region of interest.

Quantifying differences in land use emission estimates implied by definition discrepancies

The quantification of CO2 emissions from anthropogenic land use and land use change (eLUC) is essential to understand the drivers of the atmospheric CO2 increase and to inform climate change mitigation policy. Reported values in synthesis reports are commonly derived from different approaches (observation-driven bookkeeping and process-modelling) but recent work has emphasized that inconsistencies between methods may imply substantial differences in eLUC estimates. However, a consistent quantification is lacking and no concise modelling protocol for the separation of primary and secondary components of eLUC has been established. Here, we review differences of eLUC quantification methods and apply an Earth System Model (ESM) of Intermediate Complexity to quantify them. We find that the magnitude of effects due to merely conceptual differences between ESM and offline vegetation model-based quantifications is ~ 20 % for today. Under a future business-as-usual scenario, differences tend to increase further due to slowing land conversion rates and an increasing impact of altered environmental conditions on land-atmosphere fluxes. We establish how coupled Earth System Models may be applied to separate secondary component fluxes of eLUC arising from the replacement of potential C sinks/sources and the land use feedback and show that secondary fluxes derived from offline vegetation models are conceptually and quantitatively not identical to either, nor their sum. Therefore, we argue that synthesis studies should resort to the "least common denominator" of different methods, following the bookkeeping approach where only primary land use emissions are quantified under the assumption of constant environmental boundary conditions.

Zeeman effect in atmospheric O₂ measured by ground-based microwave radiometry

In this work we study the Zeeman effect on stratospheric O₂ using ground-based microwave radiometer measurements. The interaction of the Earth magnetic field with the oxygen dipole leads to a splitting of O₂ energy states, which polarizes the emission spectra. A special campaign was carried out in order to measure this effect in the oxygen emission line centered at 53.07 GHz. Both a fixed and a rotating mirror were incorporated into the TEMPERA (TEMPERature RAdiometer) in order to be able to measure under different observational angles. This new configuration allowed us to change the angle between the observational path and the Earth magnetic field direction. Moreover, a high-resolution spectrometer (1 kHz) was used in order to measure for the first time the polarization state of the radiation due to the Zeeman effect in the main isotopologue of oxygen from ground-based microwave measurements. The measured spectra showed a clear polarized signature when the observational angles were changed, evidencing the Zeeman effect in the oxygen molecule. In addition, simulations carried out with the Atmospheric Radiative Transfer Simulator (ARTS) allowed us to verify the microwave measurements showing a very good agreement between model and measurements. The results suggest some interesting new aspects for research of the upper atmosphere.

Atmospheric histories and global emissions of halons H-1211 (CBrClF 2 ), H-1301 (CBrF 3 ), and H-2402 (CBrF 2 CBrF 2 )

We report ground-based atmospheric measurements and emission estimates for the halons H-1211 (CBrClF₂), H-1301 (CBrF₃), and H-2402 (CBrF₂CBrF₂) from the AGAGE (Advanced Global Atmospheric Gases Experiment) and the National Oceanic and Atmospheric Administration global networks. We also include results from archived air samples in canisters and from polar firn in both hemispheres, thereby deriving an atmospheric record of nearly nine decades (1930s to present). All three halons were absent from the atmosphere until ~1970, when their atmospheric burdens started to increase rapidly. In recent years H-1211 and H-2402 mole fractions have been declining, but H-1301 has continued to grow. High-frequency observations show continuing emissions of H-1211 and H-1301 near most AGAGE sites. For H-2402 the only emissions detected were derived from the region surrounding the Sea of Japan/East Sea. Based on our observations, we derive global emissions using two different inversion approaches. Emissions for H-1211 declined from a peak of 11 kt yr⁻¹ (late 1990s) to 3.9 kt yr⁻¹ at the end of our record (mean of 2013–2015), for H-1301 from 5.4 kt yr⁻¹ (late 1980s) to 1.6 kt yr⁻¹, and for H-2402 from 1.8 kt yr⁻¹ (late 1980s) to 0.38 kt yr⁻¹. Yearly summed halon emissions have decreased substantially; nevertheless, since 2000 they have accounted for ~30% of the emissions of all major anthropogenic ozone depletion substances, when weighted by ozone depletion potentials.