Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1 clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H2SO4(H2O)n where n = 1–6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller–Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)n clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H2SO4(H2O)n clusters are favorable in colder regions of the troposphere (T = 216.65–273.15 K) for n = 1–6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H2SO4–H2O system must contain more than one H2SO4 and are in concert with recent findings(1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Observations of middle atmospheric H2O and O3 during the 2010 major sudden stratospheric warming by a network of microwave radiometers

In this study, we present middle atmospheric water vapor (H2O) and ozone (O3) measurements obtained by ground-based microwave radiometers at three European locations in Bern (47° N), Onsala (57° N) and Sodankylä (67° N) during Northern winter 2009/2010. In January 2010, a major sudden stratospheric warming (SSW) occurred in the Northern Hemisphere whose signatures are evident in the ground-based observations of H2O and O3. The observed anomalies in H2O and O3 are mostly explained by the relative location of the polar vortex with respect to the measurement locations. The SSW started on 26 January 2010 and was most pronounced by the end of January. The zonal mean temperature in the middle stratosphere (10 hPa) increased by approximately 25 Kelvin within a few days. The stratospheric vortex weakened during the SSW and shifted towards Europe. In the mesosphere, the vortex broke down, which lead to large scale mixing of polar and midlatitudinal air. After the warming, the polar vortex in the stratosphere split into two weaker vortices and in the mesosphere, a new, pole-centered vortex formed with maximum wind speed of 70 m s−1 at approximately 40° N. The shift of the stratospheric vortex towards Europe was observed in Bern as an increase in stratospheric H2O and a decrease in O3. The breakdown of the mesospheric vortex during the SSW was observed at Onsala and Sodankylä as a sudden increase in mesospheric H2O. The following large-scale descent inside the newly formed mesospheric vortex was well captured by the H2O observations in Sodankylä. In order to combine the H2O observations from the three different locations, we applied the trajectory mapping technique on our H2O observations to derive synoptic scale maps of the H2O distribution. Based on our observations and the 3-D wind field, this method allows determining the approximate development of the stratospheric and mesospheric polar vortex and demonstrates the potential of a network of ground-based instruments.

Forcing of stratospheric chemistry and dynamics during the Dalton Minimum

The response of atmospheric chemistry and dynamics to volcanic eruptions and to a decrease in solar activity during the Dalton Minimum is investigated with the fully coupled atmosphere–ocean chemistry general circulation model SOCOL-MPIOM (modeling tools for studies of SOlar Climate Ozone Links-Max Planck Institute Ocean Model) covering the time period 1780 to 1840 AD. We carried out several sensitivity ensemble experiments to separate the effects of (i) reduced solar ultra-violet (UV) irradiance, (ii) reduced solar visible and near infrared irradiance, (iii) enhanced galactic cosmic ray intensity as well as less intensive solar energetic proton events and auroral electron precipitation, and (iv) volcanic aerosols. The introduced changes of UV irradiance and volcanic aerosols significantly influence stratospheric dynamics in the early 19th century, whereas changes in the visible part of the spectrum and energetic particles have smaller effects. A reduction of UV irradiance by 15%, which represents the presently discussed highest estimate of UV irradiance change caused by solar activity changes, causes global ozone decrease below the stratopause reaching as much as 8% in the midlatitudes at 5 hPa and a significant stratospheric cooling of up to 2 °C in the mid-stratosphere and to 6 °C in the lower mesosphere. Changes in energetic particle precipitation lead only to minor changes in the yearly averaged temperature fields in the stratosphere. Volcanic aerosols heat the tropical lower stratosphere, allowing more water vapour to enter the tropical stratosphere, which, via HOx reactions, decreases upper stratospheric and mesospheric ozone by roughly 4%. Conversely, heterogeneous chemistry on aerosols reduces stratospheric NOx, leading to a 12% ozone increase in the tropics, whereas a decrease in ozone of up to 5% is found over Antarctica in boreal winter. The linear superposition of the different contributions is not equivalent to the response obtained in a simulation when all forcing factors are applied during the Dalton Minimum (DM) – this effect is especially well visible for NOx/NOy. Thus, this study also shows the non-linear behaviour of the coupled chemistry-climate system. Finally, we conclude that especially UV and volcanic eruptions dominate the changes in the ozone, temperature and dynamics while the NOx field is dominated by the energetic particle precipitation. Visible radiation changes have only very minor effects on both stratospheric dynamics and chemistry.

Climate and chemistry effects of a regional scale nuclear conflict

Previous studies have highlighted the severity of detrimental effects for life on earth after an assumed regionally limited nuclear war. These effects are caused by climatic, chemical and radiative changes persisting for up to one decade. However, so far only a very limited number of climate model simulations have been performed, giving rise to the question how realistic previous computations have been. This study uses the coupled chemistry climate model (CCM) SOCOL, which belongs to a different family of CCMs than previously used, to investigate the consequences of such a hypothetical nuclear conflict. In accordance with previous studies, the present work assumes a scenario of a nuclear conflict between India and Pakistan, each applying 50 warheads with an individual blasting power of 15 kt ("Hiroshima size") against the major population centers, resulting in the emission of tiny soot particles, which are generated in the firestorms expected in the aftermath of the detonations. Substantial uncertainties related to the calculation of likely soot emissions, particularly concerning assumptions of target fuel loading and targeting of weapons, have been addressed by simulating several scenarios, with soot emissions ranging from 1 to 12 Tg. Their high absorptivity with respect to solar radiation leads to a rapid self-lofting of the soot particles into the strato- and mesosphere within a few days after emission, where they remain for several years. Consequently, the model suggests earth's surface temperatures to drop by several degrees Celsius due to the shielding of solar irradiance by the soot, indicating a major global cooling. In addition, there is a substantial reduction of precipitation lasting 5 to 10 yr after the conflict, depending on the magnitude of the initial soot release. Extreme cold spells associated with an increase in sea ice formation are found during Northern Hemisphere winter, which expose the continental land masses of North America and Eurasia to a cooling of several degrees. In the stratosphere, the strong heating leads to an acceleration of catalytic ozone loss and, consequently, to enhancements of UV radiation at the ground. In contrast to surface temperature and precipitation changes, which show a linear dependence to the soot burden, there is a saturation effect with respect to stratospheric ozone chemistry. Soot emissions of 5 Tg lead to an ozone column reduction of almost 50% in northern high latitudes, while emitting 12 Tg only increases ozone loss by a further 10%. In summary, this study, though using a different chemistry climate model, corroborates the previous investigations with respect to the atmospheric impacts. In addition to these persistent effects, the present study draws attention to episodically cold phases, which would likely add to the severity of human harm worldwide. The best insurance against such a catastrophic development would be the delegitimization of nuclear weapons.

Diurnal variations in middle-atmospheric water vapor by ground-based microwave radiometry

Abstract. In this paper, we compare the diurnal variations in middle-atmospheric water vapor as measured by two ground-based microwave radiometers in the Alpine region near Bern, Switzerland. The observational data set is also compared to data from the chemistry–climate model WACCM. Due to the small diurnal variations of usually less than 1%, averages over extended time periods are required. Therefore, two time periods of five months each, December to April and June to October, were taken for the comparison. The diurnal variations from the observational data agree well with each other in amplitude and phase. The linear correlation coefficients range from 0.8 in the upper stratosphere to 0.5 in the upper mesosphere. The observed diurnal variability is significant at all pressure levels within the sensitivity of the instruments. Comparing our observations with WACCM, we find that the agreement of the phase of the diurnal cycle between observations and model is better from December to April than from June to October. The amplitudes of the diurnal variations for both time periods increase with altitude in WACCM, but remain approximately constant at 0.05 ppm in the observations. The WACCM data are used to separate the processes that lead to diurnal variations in middle-atmospheric water vapor above Bern. The dominating processes were found to be meridional advection below 0.1 hPa, vertical advection between 0.1 and 0.02 hPa and (photo-)chemistry above 0.02 hPa. The contribution of zonal advection is small. The highest diurnal variations in water vapor as seen in the WACCM data are found in the mesopause region during the time period from June to October with diurnal amplitudes of 0.2 ppm (approximately 5% in relative units).

Emerging Areas in Atmospheric Photochemistry

Sunlight is a major driving force of atmospheric processes. A detailed knowledge of atmospheric photochemistry is therefore required in order to understand atmospheric chemistry and climate. Considerable progress has been made in this field in recent decades. This contribution will highlight a set of new and emerging ideas (and will therefore not provide a complete review of the field) mainly dealing with long wavelength photochemistry both in the gas phase and on a wide range of environmental surfaces. Besides this, some interesting bulk photochemistry processes are discussed. Altogether these processes have the potential to introduce new chemical pathways into tropospheric chemistry and may impact atmospheric radical formation.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

The impacts of volcanic aerosol on stratospheric ozone and the Northern Hemisphere polar vortex: separating radiative-dynamical changes from direct effects due to enhanced aerosol heterogeneous chemistry

After major volcanic eruptions the enhanced aerosol causes ozone changes due to greater heterogeneous chemistry on the particle surfaces (HET-AER) and from dynamical effects related to the radiative heating of the lower stratosphere (RAD-DYN). We carry out a series of experiments with an atmosphere–ocean–chemistry–climate model to assess how these two processes change stratospheric ozone and Northern Hemispheric (NH) polar vortex dynamics. Ensemble simulations are performed under present day and preindustrial conditions, and with aerosol forcings representative of different eruption strength, to investigate changes in the response behaviour. We show that the halogen component of the HET-AER effect dominates under present-day conditions with a global reduction of ozone (−21 DU for the strongest eruption) particularly at high latitudes, whereas the HET-AER effect increases stratospheric ozone due to N2O5 hydrolysis in a preindustrial atmosphere (maximum anomalies +4 DU). The halogen-induced ozone changes in the present-day atmosphere offset part of the strengthening of the NH polar vortex during mid-winter (reduction of up to −16 m s-1 in January) and slightly amplify the dynamical changes in the polar stratosphere in late winter (+11 m s-1 in March). The RAD-DYN mechanism leads to positive column ozone anomalies which are reduced in a present-day atmosphere by amplified polar ozone depletion (maximum anomalies +12 and +18 DU for present day and preindustrial, respectively). For preindustrial conditions, the ozone response is consequently dominated by RAD-DYN processes, while under present-day conditions, HET-AER effects dominate. The dynamical response of the stratosphere is dominated by the RAD-DYN mechanism showing an intensification of the NH polar vortex in winter (up to +10 m s-1 in January). Ozone changes due to the RAD-DYN mechanism slightly reduce the response of the polar vortex after the eruption under present-day conditions.

Evaluating atmospheric methane inversion model results for Pallas, northern Finland

A state-of-the-art inverse model, CarbonTracker Data Assimilation Shell (CTDAS), was used to optimize estimates of methane (CH4) surface fluxes using atmospheric observations of CH4 as a constraint. The model consists of the latest version of the TM5 atmospheric chemistry-transport model and an ensemble Kalman filter based data assimilation system. The model was constrained by atmospheric methane surface concentrations, obtained from the World Data Centre for Greenhouse Gases (WDCGG). Prior methane emissions were specified for five sources: biosphere, anthropogenic, fire, termites and ocean, of which bio-sphere and anthropogenic emissions were optimized. Atmospheric CH 4 mole fractions for 2007 from northern Finland calculated from prior and optimized emissions were compared with observations. It was found that the root mean squared errors of the posterior esti - mates were more than halved. Furthermore, inclusion of NOAA observations of CH 4 from weekly discrete air samples collected at Pallas improved agreement between posterior CH 4 mole fraction estimates and continuous observations, and resulted in reducing optimized biosphere emissions and their uncertainties in northern Finland.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.