A global historical ozone data set and prominent features of stratospheric variability prior to 1979

We present a vertically resolved zonal mean monthly mean global ozone data set spanning the period 1901 to 2007, called HISTOZ.1.0. It is based on a new approach that combines information from an ensemble of chemistry climate model (CCM) simulations with historical total column ozone information. The CCM simulations incorporate important external drivers of stratospheric chemistry and dynamics (in particular solar and volcanic effects, greenhouse gases and ozone depleting substances, sea surface temperatures, and the quasi-biennial oscillation). The historical total column ozone observations include ground-based measurements from the 1920s onward and satellite observations from 1970 to 1976. An off-line data assimilation approach is used to combine model simulations, observations, and information on the observation error. The period starting in 1979 was used for validation with existing ozone data sets and therefore only ground-based measurements were assimilated. Results demonstrate considerable skill from the CCM simulations alone. Assimilating observations provides additional skill for total column ozone. With respect to the vertical ozone distribution, assimilating observations increases on average the correlation with a reference data set, but does not decrease the mean squared error. Analyses of HISTOZ.1.0 with respect to the effects of El Niño–Southern Oscillation (ENSO) and of the 11 yr solar cycle on stratospheric ozone from 1934 to 1979 qualitatively confirm previous studies that focussed on the post-1979 period. The ENSO signature exhibits a much clearer imprint of a change in strength of the Brewer–Dobson circulation compared to the post-1979 period. The imprint of the 11 yr solar cycle is slightly weaker in the earlier period. Furthermore, the total column ozone increase from the 1950s to around 1970 at northern mid-latitudes is briefly discussed. Indications for contributions of a tropospheric ozone increase, greenhouse gases, and changes in atmospheric circulation are found. Finally, the paper points at several possible future improvements of HISTOZ.1.0.

Influence of the sunspot cycle on the Northern Hemisphere wintertime circulation from long upper-air data sets

Here we present a study of the 11 yr sunspot cycle's imprint on the Northern Hemisphere atmospheric circulation, using three recently developed gridded upper-air data sets that extend back to the early twentieth century. We find a robust response of the tropospheric late-wintertime circulation to the sunspot cycle, independent from the data set. This response is particularly significant over Europe, although results show that it is not directly related to a North Atlantic Oscillation (NAO) modulation; instead, it reveals a significant connection to the more meridional Eurasian pattern (EU). The magnitude of mean seasonal temperature changes over the European land areas locally exceeds 1 K in the lower troposphere over a sunspot cycle. We also analyse surface data to address the question whether the solar signal over Europe is temporally stable for a longer 250 yr period. The results increase our confidence in the existence of an influence of the 11 yr cycle on the European climate, but the signal is much weaker in the first half of the period compared to the second half. The last solar minimum (2005 to 2010), which was not included in our analysis, shows anomalies that are consistent with our statistical results for earlier solar minima.

Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

Method for the determination of specific molecular markers of biomass burning in lake sediments

Fire has an influence on regional to global atmospheric chemistry and climate. Molecular markers of biomass burning archived in lake sediments are becoming increasingly important in paleoenvironmental reconstruction and may help determine the interaction between climate and fire activity. Here, we present a high performance anion exchange chromatography–mass spectrometry method to allow separation and analysis of levoglucosan, mannosan and galactosan in lake sediments, with implications for reconstructing past biomass burning events. Determining mannosan and galactosan in Lake Kirkpatrick, New Zealand (45.03°S, 168.57°E) sediment cores and comparing these isomers with the more abundant biomass burning markers levoglucosan and charcoal represents a significant advancement in our ability to analyze past fire activity. Levoglucosan, mannosan and galactosan concentrations correlated significantly with macroscopic charcoal concentration. Levoglucosan/mannosan and levoglucosan/(mannosan + galactosan) ratios may help determine not only when fires occurred, but also if changes in the primary burned vegetation occurred.

The quasi 16-day wave in mesospheric water vapor during boreal winter 2011/2012

This study investigates the characteristics of the quasi 16-day wave in the mesosphere during boreal winter 2011/2012 using observations of water vapor from ground-based microwave radiometers and satellite data. The ground-based microwave radiometers are located in Seoul (South Korea, 37° N), Bern (Switzerland, 47° N) and Sodankylä (Finland, 67° N). The quasi 16-day wave is observed in the mesosphere at all three locations, while the dominant period increases with latitude from 15 days at Seoul to 20 days at Sodankylä. The observed evolution of the quasi 16-day wave confirms that the wave activity is strongly decreased during a sudden stratospheric warming that occurred in mid-January 2012. Using satellite data from the Microwave Limb Sounder on the Aura satellite, we examine the zonal characteristics of the quasi 16-day wave and conclude that the observed waves above the mid-latitudinal stations Seoul and Bern are eastward-propagating s=−1 planetary waves with periods of 15 to 16 days, while the observed oscillation above the polar station Sodankylä is a standing oscillation with a period of approximately 20 days. The strongest relative wave amplitudes in water vapor during the investigated time period are approximately 15%. The wave activity varies strongly along a latitude circle. The activity of the quasi 16-day wave in mesospheric water vapor during boreal winter 2011/2012 is strongest over Northern Europe, the North Atlantic ocean and North-West Canada. The region of highest wave activity seems to be related to the position of the polar vortex. We conclude that the classic approach to characterize planetary waves zonally averaged along a latitude circle is not sufficient to explain the local observations because of the strong longitudinal dependence of the wave activity.

A climatology of the diurnal variations of stratospheric and mesospheric ozone over Bern, Switzerland

The ground-based radiometer GROMOS, stationed in Bern (47.95° N, 7.44° E), Switzerland, has a unique dataset: it obtains ozone profiles from November 1994 to present with a time resolution of 30 min and equal quality during night- and daytime. Here, we derive a monthly climatology of the daily ozone cycle from 17 yr of GROMOS observation. We present the diurnal ozone variation of the stratosphere and mesosphere. Characterizing the diurnal cycle of stratospheric ozone is important for correct trend estimates of the ozone layer derived from satellite observations. The diurnal ozone cycle from GROMOS is compared to two models: The Whole Atmosphere Community Climate Model (WACCM) and the Hamburg Model of Neutral and Ionized Atmosphere (HAMMONIA). Aura Microwave Limb Sounder (Aura/MLS) ozone data, from night- and daytime overpasses over Bern, have also been included in the comparison. Generally, observation and models show good qualitative agreement: in the lower mesosphere, daytime ozone is for both GROMOS and models around 25% less than nighttime ozone (reference is 22:30–01:30). In the stratosphere, ozone reaches its maximum in the afternoon showing values several percent larger than the midnight value. It is important that diurnal ozone variations of this order are taken into account when merging different data sets for the derivation of long-term ozone trends in the stratosphere. Further, GROMOS and models indicate a seasonal behavior of daily ozone variations in the stratosphere with a larger afternoon maximum during daytime in summer than in winter. At 0.35 hPa, observations from GROMOS and Aura/MLS show a seasonal pattern in diurnal ozone variations with larger relative amplitudes during daytime in winter (−25 ± 5%) than in summer (−18 ± 4%) (compared to mean values around midnight). For the first time, a time series of the diurnal variations in ozone is presented: 17 yr of GROMOS data show strong interannual variations in the diurnal ozone cycle for both the stratosphere and the mesosphere. There are some indications that strong temperature tides can suppress the diurnal variation of stratospheric ozone via the anticorrelation of temperature and ozone. That means the spatio-temporal variability of solar thermal tides seems to affect the diurnal cycle of stratospheric ozone.

Precipitation isoscape of high reliefs: interpolation scheme designed and tested for monthly resolved precipitation oxygen isotope records of an Alpine domain

Stable oxygen isotope composition of atmospheric precipitation (δ18Op) was scrutinized from 39 stations distributed over Switzerland and its border zone. Monthly amount-weighted δ18Op values averaged over the 1995–2000 period showed the expected strong linear altitude dependence (−0.15 to −0.22‰ per 100 m) only during the summer season (May–September). Steeper gradients (~ −0.56 to −0.60‰ per 100 m) were observed for winter months over a low elevation belt, while hardly any altitudinal difference was seen for high elevation stations. This dichotomous pattern could be explained by the characteristically shallower vertical atmospheric mixing height during winter season and provides empirical evidence for recently simulated effects of stratified atmospheric flow on orographic precipitation isotopic ratios. This helps explain "anomalous" deflected altitudinal water isotope profiles reported from many other high relief regions. Grids and isotope distribution maps of the monthly δ18Op have been calculated over the study region for 1995–1996. The adopted interpolation method took into account both the variable mixing heights and the seasonal difference in the isotopic lapse rate and combined them with residual kriging. The presented data set allows a point estimation of δ18Op with monthly resolution. According to the test calculations executed on subsets, this biannual data set can be extended back to 1992 with maintained fidelity and, with a reduced station subset, even back to 1983 at the expense of faded reliability of the derived δ18Op estimates, mainly in the eastern part of Switzerland. Before 1983, reliable results can only be expected for the Swiss Plateau since important stations representing eastern and south-western Switzerland were not yet in operation.

Daily ozone cycle in the stratosphere: global, regional and seasonal behaviour modelled with the Whole Atmosphere Community Climate Model

The Whole Atmosphere Community Climate Model (WACCM) is utilised to study the daily ozone cycle and underlying photochemical and dynamical processes. The analysis is focused on the daily ozone cycle in the middle stratosphere at 5 hPa where satellite-based trend estimates of stratospheric ozone are most biased by diurnal sampling effects and drifting satellite orbits. The simulated ozone cycle shows a minimum after sunrise and a maximum in the late afternoon. Further, a seasonal variation of the daily ozone cycle in the stratosphere was found. Depending on season and latitude, the peak-to-valley difference of the daily ozone cycle varies mostly between 3 and 5% (0.4 ppmv) with respect to the midnight ozone volume mixing ratio. The maximal variation of 15% (0.8 ppmv) is found at the polar circle in summer. The global pattern of the strength of the daily ozone cycle is mainly governed by the solar zenith angle and the sunshine duration. In addition, we find synoptic-scale variations in the strength of the daily ozone cycle. These variations are often anti-correlated to regional temperature anomalies and are due to the temperature dependence of the rate coefficients k2 and k3 of the Chapman cycle reactions. Further, the NOx catalytic cycle counteracts the accumulation of ozone during daytime and leads to an anti-correlation between anomalies in NOx and the strength of the daily ozone cycle. Similarly, ozone recombines with atomic oxygen which leads to an anti-correlation between anomalies in ozone abundance and the strength of the daily ozone cycle. At higher latitudes, an increase of the westerly (easterly) wind cause a decrease (increase) in the sunshine duration of an air parcel leading to a weaker (stronger) daily ozone cycle.

The quasi 16-day wave in mesospheric water vapor during boreal winter 2011/2012

This study investigates the characteristics of the quasi 16-day wave in the mesosphere during boreal winter 2011/2012 using observations of water vapor from ground-based microwave radiometers and satellite data. The ground-based microwave radiometers are located in Seoul (South Korea, 37° N), Bern (Switzerland, 47° N) and Sodankylä (Finland, 67° N). The quasi 16-day wave is observed in the mesosphere at all three locations, while the dominant period increases with latitude from 15 days at Seoul to 20 days at Sodankylä. The observed evolution of the quasi 16-day wave confirms that the wave activity is strongly decreased during a sudden stratospheric warming that occurred in mid-January 2012. Using satellite data from the Microwave Limb Sounder on the Aura satellite, we examine the zonal characteristics of the quasi 16-day wave and conclude that the observed waves above the midlatitudinal stations Seoul and Bern are eastward-propagating s = −1 planetary waves with periods of 15 to 16 days, while the observed oscillation above the polar station Sodankylä is a standing wave with a period of approximately 20 days. The strongest relative wave amplitudes in water vapor during the investigated time period are approximately 15%. The wave activity varies strongly along a latitude circle. The activity of the quasi 16-day wave in mesospheric water vapor during boreal winter 2011/2012 is strongest over northern Europe, the North Atlantic Ocean and northwestern Canada. The region of highest wave activity seems to be related to the position of the polar vortex. We conclude that the classic approach to characterize planetary waves zonally averaged along a latitude circle is not sufficient to explain the local observations because of the strong longitudinal dependence of the wave activity.

A climatology of the diurnal variations of stratospheric and mesospheric ozone over Bern, Switzerland

The ground-based radiometer GROMOS, stationed in Bern (47.95° N, 7.44° E), Switzerland, has a unique dataset: it obtains ozone profiles from November 1994 to present with a time resolution of 30 min and equal quality during night- and daytime. Here, we derive a monthly climatology of the daily ozone cycle from 17 yr of GROMOS observation. We present the diurnal ozone variation of the stratosphere and mesosphere. Characterizing the diurnal cycle of stratospheric ozone is important for correct trend estimates of the ozone layer derived from satellite observations. The diurnal ozone cycle from GROMOS is compared to two models: The Whole Atmosphere Community Climate Model (WACCM) and the Hamburg Model of Neutral and Ionized Atmosphere (HAMMONIA). Aura Microwave Limb Sounder (Aura/MLS) ozone data, from night- and daytime overpasses over Bern, have also been included in the comparison. Generally, observation and models show good qualitative agreement: in the lower mesosphere, daytime ozone is for both GROMOS and models around 25% less than nighttime ozone (reference is 22:30–01:30). In the stratosphere, ozone reaches its maximum in the afternoon showing values several percent larger than the midnight value. It is important that diurnal ozone variations of this order are taken into account when merging different data sets for the derivation of long-term ozone trends in the stratosphere. Further, GROMOS and models indicate a seasonal behavior of daily ozone variations in the stratosphere with a larger afternoon maximum during daytime in summer than in winter. At 0.35 hPa, observations from GROMOS and Aura/MLS show a seasonal pattern in diurnal ozone variations with larger relative amplitudes during daytime in winter (−25 ± 5%) than in summer (−18 ± 4%) (compared to mean values around midnight). For the first time, a time series of the diurnal variations in ozone is presented: 17 yr of GROMOS data show strong interannual variations in the diurnal ozone cycle for both the stratosphere and the mesosphere. There are some indications that strong temperature tides can suppress the diurnal variation of stratospheric ozone via the anticorrelation of temperature and ozone. That means the spatio-temporal variability of solar thermal tides seems to affect the diurnal cycle of stratospheric ozone.

The diurnal variation in stratospheric ozone from the MACC reanalysis, the ERA-Interim reanalysis, WACCM and Earth observations: characteristics and intercomparison

In this study we compare the diurnal variation in stratospheric ozone derived from free-running simulations of the Whole Atmosphere Community Climate Model (WACCM) and from reanalysis data of the atmospheric service MACC (Monitoring Atmospheric Composition and Climate) which both use a similar stratospheric chemistry module. We find good agreement between WACCM and the MACC reanalysis for the diurnal ozone variation in the high-latitude summer stratosphere based on photochemistry. In addition, we consult the ozone data product of the ERA-Interim reanalysis. The ERA-Interim reanalysis ozone system with its long-term ozone parametrization can not capture these diurnal variations in the upper stratosphere that are due to photochemistry. The good dynamics representations, however, reflects well dynamically induced ozone variations in the lower stratosphere. For the high-latitude winter stratosphere we describe a novel feature of diurnal variation in ozone where changes of up to 46.6% (3.3 ppmv) occur in monthly mean data. For this effect good agreement between the ERA-Interim reanalysis and the MACC reanalysis suggest quite similar diurnal advection processes of ozone. The free-running WACCM model seriously underestimates the role of diurnal advection processes at the polar vortex at the two tested resolutions. The intercomparison of the MACC reanalysis and the ERA-Interim reanalysis demonstrates how global reanalyses can benefit from a chemical representation held by a chemical transport model. The MACC reanalysis provides an unprecedented description of the dynamics and photochemistry of the diurnal variation of stratospheric ozone which is of high interest for ozone trend analysis and research on atmospheric tides. We confirm the diurnal variation in ozone at 5 hPa by observations of the Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) experiment and selected sites of the Network for Detection of Atmospheric Composition Change (NDACC). The latter give valuable insight even to diurnal variation of ozone in the polar winter stratosphere.

Radiocarbon Analysis of Elemental and Organic Carbon in Switzerland during Winter-Smog Episodes from 2008 to 2012 – Part I: Source Apportionment and Spatial Variability

While several studies have investigated winter-time air pollution with a wide range of concentration levels, hardly any results are available for longer time periods covering several winter-smog episodes at various locations; e.g., often only a few weeks from a single winter are investigated. Here, we present source apportionment results of winter-smog episodes from 16 air pollution monitoring stations across Switzerland from five consecutive winters. Radiocarbon (14C) analyses of the elemental (EC) and organic (OC) carbon fractions, as well as levoglucosan, major water-soluble ionic species and gas-phase pollutant measurements were used to characterize the different sources of PM10. The most important contributions to PM10 during winter-smog episodes in Switzerland were on average the secondary inorganic constituents (sum of nitrate, sulfate and ammonium = 41 ± 15%) followed by organic matter (OM) (34 ± 13%) and EC (5 ± 2%). The non-fossil fractions of OC (fNF,OC) ranged on average from 69 to 85 and 80 to 95% for stations north and south of the Alps, respectively, showing that traffic contributes on average only up to ~ 30% to OC. The non-fossil fraction of EC (fNF,EC), entirely attributable to primary wood burning, was on average 42 ± 13 and 49 ± 15% for north and south of the Alps, respectively. While a high correlation was observed between fossil EC and nitrogen oxides, both primarily emitted by traffic, these species did not significantly correlate with fossil OC (OCF), which seems to suggest that a considerable amount of OCF is secondary, from fossil precursors. Elevated fNF,EC and fNF,OC values and the high correlation of the latter with other wood burning markers, including levoglucosan and water soluble potassium (K+) indicate that residential wood burning is the major source of carbonaceous aerosols during winter-smog episodes in Switzerland. The inspection of the non-fossil OC and EC levels and the relation with levoglucosan and water-soluble K+ shows different ratios for stations north and south of the Alps (most likely because of differences in burning technologies) for these two regions in Switzerland.

Fossil vs. non-fossil sources of fine carbonaceous aerosols in four Chinese cities during the extreme winter haze episode of 2013

During winter 2013, extremely high concentrations (i.e., 4–20 times higher than the World Health Organization guideline) of PM2.5 (particulate matter with an aerodynamic diameter < 2.5 μm) mass concentrations (24 h samples) were found in four major cities in China including Xi'an, Beijing, Shanghai and Guangzhou. Statistical analysis of a combined data set from elemental carbon (EC), organic carbon (OC), 14C and biomass-burning marker measurements using Latin hypercube sampling allowed a quantitative source apportionment of carbonaceous aerosols. Based on 14C measurements of EC fractions (six samples each city), we found that fossil emissions from coal combustion and vehicle exhaust dominated EC with a mean contribution of 75 ± 8% across all sites. The remaining 25 ± 8% was exclusively attributed to biomass combustion, consistent with the measurements of biomass-burning markers such as anhydrosugars (levoglucosan and mannosan) and water-soluble potassium (K+). With a combination of the levoglucosan-to-mannosan and levoglucosan-to-K+ ratios, the major source of biomass burning in winter in China is suggested to be combustion of crop residues. The contribution of fossil sources to OC was highest in Beijing (58 ± 5%) and decreased from Shanghai (49 ± 2%) to Xi'an (38 ± 3%) and Guangzhou (35 ± 7%). Generally, a larger fraction of fossil OC was from secondary origins than primary sources for all sites. Non-fossil sources accounted on average for 55 ± 10 and 48 ± 9% of OC and total carbon (TC), respectively, which suggests that non-fossil emissions were very important contributors of urban carbonaceous aerosols in China. The primary biomass-burning emissions accounted for 40 ± 8, 48 ± 18, 53 ± 4 and 65 ± 26% of non-fossil OC for Xi'an, Beijing, Shanghai and Guangzhou, respectively. Other non-fossil sources excluding primary biomass burning were mainly attributed to formation of secondary organic carbon (SOC) from non-fossil precursors such as biomass-burning emissions. For each site, we also compared samples from moderately to heavily polluted days according to particulate matter mass. Despite a significant increase of the absolute mass concentrations of primary emissions from both fossil and non-fossil sources during the heavily polluted events, their relative contribution to TC was even decreased, whereas the portion of SOC was consistently increased at all sites. This observation indicates that SOC was an important fraction in the increment of carbonaceous aerosols during the haze episode in China.