Search Results
Columnar modelling of nucleation burst evolution in the convective boundary layer - First results from a feasibility study, Part III: Preliminary results on physicochemical model performance using two "clean air mass" reference scenarios
2006, Hellmuth, O.
In Paper I of four papers, a revised columnar high-order model to investigate gas-aerosol-turbulence interactions in the convective boundary layer (CBL) was proposed. In Paper II, the model capability to predict first-, second- and third-order moments of meteorological variables in the CBL was demonstrated using available observational data. In the present Paper III, the high-order modelling concept is extended to sulphur and ammonia chemistry as well as to aerosol dynamics. Based on the previous CBL simulation, a feasibility study is performed using two "clean air mass" scenarios with an emission source at the ground but low aerosol background concentration. Such scenarios synoptically correspond to the advection of fresh post-frontal air in an anthropogenically influenced region. The aim is to evaluate the time-height evolution of ultrafine condensation nuclei (UCNs) and to elucidate the interactions between meteorological and physicochemical variables in a CBL column. The scenarios differ in the treatment of new particle formation (NPF), whereas homogeneous nucleation according to the classical nucleation theory (CNT) is considered. The first scenario considers nucleation of a binary system consisting of water vapour and sulphuric acid (H2SO4) vapour, the second one nucleation of a ternary system additionally involving ammonia (NH3). Here, the two synthetic scenarios are discussed in detail, whereas special attention is payed to the role of turbulence in the formation of the typical UCN burst behaviour, that can often be observed in the surface layer. The intercomparison of the two scenarios reveals large differences in the evolution of the UCN number concentration in the surface layer as well as in the time-height cross-sections of first-order moments and double correlation terms. Although in both cases the occurrence of NPF bursts could be simulated, the burst characteristics and genesis of the bursts are completely different. It is demonstrated, that observations from the surface layer alone are not conclusive to elucidate the origin of newly formed particles. This is also true with respect to the interpretation of box modelling studies. The binary and ternary NPF bursts observed in the surface layer differ with respect to burst amplitude and phase. New particles simulated in the binary scenario are formed in the forenoon in the upper part of the growing CBL, followed by turbulence-induced top-down transport. Hence, with respect to the burst observation site in the surface layer, new particles are formed ex situ. In opposite to this, the ternary case reveals a much more complex pattern. Here, NPF is initiated in the early morning hours in the surface layer, when temperature (T) is low and relative humidity (RH), sulphur dioxide (SO2) and NH3 concentrations are high, hence new particles are formed in situ. Shortly after that, ex situ NPF in the free troposphere sets in, followed by entrainment and top-down diffusion of newly formed particles into the surface layer. Altogether, these processes mainly contribute to the formation of a strong burst in the morning hours in the ternary scenario. While the time-height cross-section of the binary nucleation rate resembles a "blob"-like evolution pattern, the ternary one resembles a "sucking tube"-like pattern. The time-height cross-sections of the flux pattern and double correlations could be plausibly interpreted in terms of CBL turbulence and entrainment/detrainment processes both in the binary and in the ternary case. Although the present approach is a pure conceptual one, it shows the feasibility to simulate gas-aerosol-turbulence interactions in the CBL. Prior to a dedicated verification/validation study, further attempts are necessary to consider a more advanced description of the formation and activation of thermodynamically stable clusters according to modern concepts proposed by Kulmala et al. (2000), Kulmala (2003) and Kulmala et al. (2004a).
Homogeneous and heterogeneous ice nucleation at LACIS: Operating principle and theoretical studies
2011, Hartmann, S., Niedermeier, D., Voigtländer, J., Clauss, T., Shaw, R.A., Wex, H., Kiselev, A., Stratmann, F.
At the Leipzig Aerosol Cloud Interaction Simulator (LACIS) experiments investigating homogeneous and heterogeneous nucleation of ice (particularly immersion freezing in the latter case) have been carried out. Here both the physical LACIS setup and the numerical model developed to design experiments at LACIS and interpret their results are presented in detail. Combining results from the numerical model with experimental data, it was found that for the experimental parameter space considered, classical homogeneous ice nucleation theory is able to predict the freezing behavior of highly diluted ammonium sulfate solution droplets, while classical heterogeneous ice nucleation theory, together with the assumption of a constant contact angle, fails to predict the immersion freezing behavior of surrogate mineral dust particles (Arizona Test Dust, ATD). The main reason for this failure is the compared to experimental data apparently overly strong temperature dependence of the nucleation rate coefficient. Assuming, in the numerical model, Classical Nucleation Theory (CNT) for homogeneous ice nucleation and a CNT-based parameterization for the nucleation rate coefficient in the immersion freezing mode, recently published by our group, it was found that even for a relatively effective ice nucleating agent such as pure ATD, there is a temperature range where homogeneous ice nucleation is dominant. The main explanation is the apparently different temperature dependencies of the two freezing mechanisms. Finally, reviewing the assumptions made during the derivation of the CNT-based parameterization for immersion freezing, it was found that the assumption of constant temperature during ice nucleation and the chosen ice nucleation time were justified, underlining the applicability of the method to determine the fitting coefficients in the parameterization equation.
Retrievals of chlorine chemistry kinetic parameters from Antarctic ClO microwave radiometer measurements
2011, Kremser, S., Schofield, R., Bodeker, G.E., Connor, B.J., Rex, M., Barret, J., Mooney, T., Salawitch, R.J., Canty, T., Frieler, K., Chipperfield, M.P., Langematz, U., Feng, W.
Key kinetic parameters governing the partitioning of chlorine species in the Antarctic polar stratosphere were retrieved from 28 days of chlorine monoxide (ClO) microwave radiometer measurements made during the late winter/early spring of 2005 at Scott Base (77.85° S, 166.75° E). During day-time the loss of the ClO dimer chlorine peroxide (ClOOCl) occurs mainly by photolysis. Some time after sunrise, a photochemical equilibrium is established and the ClO/ClOOCl partitioning is determined by the ratio of the photolysis frequency, J, and the dimer formation rate, kf. The values of J and kf from laboratory studies remain uncertain to a considerable extent, and as a complement to these ongoing studies, the goal of this work is to provide a constraint on that uncertainty based on observations of ClO profiles in the Antarctic. First an optimal estimation technique was used to derive J/kf ratios for a range of Keq values. The optimal estimation forward model was a photochemical box model that takes J, kf, and Keq as inputs, together with a priori profiles of activated chlorine (ClOx = ClO+2×ClOOCl), profiles of ozone, temperature, and pressure. JPL06 kinetics are used as a priori in the optimal estimation and for all other chemistry in the forward model. Using the more recent JPL09 kinetics results in insignificant differences in the retrieved value of J/kf. A complementary approach was used to derive the optimal kinetic parameters; the full parameter space of J, kf, Keq and ClOx was sampled to find the minimum in differences between measured and modelled ClO profiles. Furthermore, values of Keq up to 2.0 times larger than recommended by JPL06 were explored to test the sensitivity of the J/kf ratio to changes in Keq. The results show that the retrieved J/kf ratios bracket the range of 1.23 to 1.97 times the J/kf value recommended by JPL06 over the range of Keq values considered. The retrieved J/kf ratios lie in the lower half of the large uncertainty range of J/kf recommended by JPL06 and towards the upper portion of the smaller uncertainty range recommended by JPL09.
Zonal asymmetries in middle atmospheric ozone and water vapour derived from Odin satellite data 2001-2010
2011, Gabriel, A., Körnich, H., Lossow, S., Peters, D.H.W., Urban, J., Murtagh, D.
Stationary wave patterns in middle atmospheric ozone (O3) and water vapour (H2O) are an important factor in the atmospheric circulation, but there is a strong gap in diagnosing and understanding their configuration and origin. Based on Odin satellite data from 2001 to 2010 we investigate the stationary wave patterns in O3 and H2O as indicated by the seasonal long-term means of the zonally asymmetric components O3* Combining double low line O3-[O3] and H2O* Combining double low line H2O-[H2O] ([O3], [H2O]: zonal means). At mid-and polar latitudes we find a pronounced wave one pattern in both constituents. In the Northern Hemisphere, the wave patterns increase during autumn, maintain their strength during winter and decay during spring, with maximum amplitudes of about 10-20 % of the zonal mean values. During winter, the wave one in O3* shows a maximum over the North Pacific/Aleutians and a minimum over the North Atlantic/Northern Europe and a double-peak structure with enhanced amplitude in the lower and in the upper stratosphere. The wave one in H2O* extends from the lower stratosphere to the upper mesosphere with a westward shift in phase with increasing height including a jump in phase at upper stratosphere altitudes. In the Southern Hemisphere, similar wave patterns occur mainly during southern spring. By comparing the observed wave patterns in O 3* and H2O3* with a linear solution of a steady-state transport equation for a zonally asymmetric tracer component we find that these wave patterns are primarily due to zonally asymmetric transport by geostrophically balanced winds, which are derived from observed temperature profiles. In addition temperature-dependent photochemistry contributes substantially to the spatial structure of the wave pattern in O 3* . Further influences, e.g., zonal asymmetries in eddy mixing processes, are discussed.
The contribution of sulphuric acid to atmospheric particle formation and growth: A comparison between boundary layers in Northern and Central Europe
2005, Fiedler, V., Dal Maso, M., Boy, M., Aufmhoff, H., Hoffmann, J., Schuck, T., Birmili, W., Hanke, M., Uecker, J., Arnold, F., Kulmala, M.
Atmospheric gaseous sulphuric acid was measured and its influence on particle formation and growth was investigated building on aerosol data. The measurements were part of the EU-project QUEST and took place at two different measurement sites in Northern and Central Europe (Hyytiälä, Finland, March-April 2003 and Heidelberg, Germany, March-April 2004). From a comprehensive data set including sulphuric acid, particle number size distributions and meteorological data, particle growth rates, particle formation rates and source rates of condensable vapors were inferred. Growth rates were determined in two different ways, from particle size distributions as well as from a so-called timeshift analysis. Moreover, correlations between sulphuric acid and particle number concentration between 3 and 6 nm were examined and the influence of air masses of different origin was investigated. Measured maximum concentrations of sulphuric acid were in the range from 1x106 to 16x106cm-3. The gaseous sulphuric acid lifetime with respect to condensation on aerosol particles ranged from 2 to 33min in Hyytiälä and from 0.5 to 8 min in Heidelberg. Most calculated values (growth rates, formation rates, vapor source rates) were considerably higher in Central Europe (Heidelberg), due to the more polluted air and higher preexistent aerosol concentrations. Close correlations between H2SO4 and nucleation mode particles (size range: 3-6 nm) were found on most days at both sites. The percentage contribution of sulphuric acid to particle growth was below 10% at both places and to initial growth below 20%. An air mass analysis indicated that at Heidelberg new particles were formed predominantly in air advected from southwesterly directions.
Long-term behavior of the concentration of the minor constituents in the mesosphere-a model study
2009, Grygalashvyly, M., Sonnemann, G.R., Hartogh, P.
We investigate the influence the rising concentrations of methane, nitrous oxide and carbon dioxide which have occurred since the pre-industrial era, have had on the chemistry of the mesosphere. For this investigation we use our global 3-D-model COMMA-IAP which was designed for the exploration of the MLT-region and in particular the extended mesopause region. Assumptions and approximations for the trends in the Lyman-flux (needed for the water vapor dissociation rate), methane and the water vapor mixing ratio at the hygropause are necessary to accomplish this study. To approximate the solar Lyman-α flux back to the pre-industrial time, we derived a quadratic fit using the sunspot number record which extends back to 1749 and is the only solar proxy available for the Lyman-α flux prior to 1947. We assume that methane increases with a constant growth rate from the pre-industrial era to the present. An unsolved problem for the model calculations consists of how the water vapor mixing ratio at the hygropause should be specified during this period. We assume that the hygropause was dryer during pre-industrial times than the present. As a consequence of methane oxidation, the model simulation indicates that the middle atmosphere has become more humid as a result of the rising methane concentration, but with some dependence on height and with a small time delay of few years. The solar influence on the water vapor mixing ratio is insignificant below about 80 km in summer high latitudes, but becomes increasingly more important above this altitude. The enhanced water vapor concentration increasesthe hydrogen radical concentration and reduces the mesospheric ozone. A second region of stronger ozone decrease is located in the vicinity of the stratopause. Increases in CO2 concentration enhance slightly the concentration of CO in the mesosphere. However, its influence upon the chemistry is small and its main effect is connected with a cooling of the upper atmosphere. The long-term behavior of water vapor is discussed in particular with respect to its impact on the NLC region.
The chemistry of OH and HO2 radicals in the boundary layer over the tropical Atlantic Ocean
2010, Whalley, L.K., Furneaux, K.L., Goddard, A., Lee, J.D., Mahajan, A., Oetjen, H., Read, K.A., Kaaden, N., Carpenter, L.J., Lewis, A.C., Plane, J.M.C., Saltzman, E.S., Wiedensohler, A., Heard, D.E.
Fluorescence Assay by Gas Expansion (FAGE) has been used to detect ambient levels of OH and HO2 radicals at the Cape Verde Atmospheric Observatory, located in the tropical Atlantic marine boundary layer, during May and June 2007. Midday radical concentrations were high, with maximum concentrations of 9 ×106 molecule cm−3 and 6×108 molecule cm−3 observed for OH and HO2, respectively. A box model incorporating the detailed Master Chemical Mechanism, extended to include halogen chemistry, heterogeneous loss processes and constrained by all available measurements including halogen and nitrogen oxides, has been used to assess the chemical and physical parameters controlling the radical chemistry. The model was able to reproduce the daytime radical concentrations to within the 1 σ measurement uncertainty of 20% during the latter half of the measurement period but significantly under-predicted [HO2] by 39% during the first half of the project. Sensitivity analyses demonstrate that elevated [HCHO] (~2 ppbv) on specific days during the early part of the project, which were much greater than the mean [HCHO] (328 pptv) used to constrain the model, could account for a large portion of the discrepancy between modelled and measured [HO2] at this time. IO and BrO, although present only at a few pptv, constituted ~19% of the instantaneous sinks for HO2, whilst aerosol uptake and surface deposition to the ocean accounted for a further 23% of the HO2 loss at noon. Photolysis of HOI and HOBr accounted for ~13% of the instantaneous OH formation. Taking into account that halogen oxides increase the oxidation of NOx (NO → NO2), and in turn reduce the rate of formation of OH from the reaction of HO2 with NO, OH concentrations were estimated to be 9% higher overall due to the presence of halogens. The increase in modelled OH from halogen chemistry gives an estimated 9% shorter lifetime for methane in this region, and the inclusion of halogen chemistry is necessary to model the observed daily cycle of O3 destruction that is observed at the surface. Due to surface losses, we hypothesise that HO2 concentrations increase with height and therefore contribute a larger fraction of the O3 destruction than at the surface.
Reactive Halogens in the Marine Boundary Layer (RHaMBLe): The tropical North Atlantic experiments
2010, Lee, J.D., McFiggans, G., Allan, J.D., Baker, A.R., Ball, S.M., Benton, A.K., Carpenter, L.J., Commane, R., Finley, B.D., Evans, M., Fuentes, E., Furneaux, K., Goddard, A., Good, N., Hamilton, J.F., Heard, D.E., Herrmann, H., Hollingsworth, A., Hopkins, J.R., Ingham, T., Irwin, M., Jones, C.E., Jones, R.L., Keene, W.C., Lawler, M.J., Lehmann, S., Lewis, A.C., Long, M.S., Mahajan, A., Methven, J., Moller, S.J., Müller, K., Müller, T., Niedermeier, N., O'Doherty, S., Oetjen, H., Plane, J.M.C., Pszenny, A.A.P., Read, K.A., Saiz-Lopez, A., Saltzman, E.S., Sander, R., von Glasow, R., Whalley, L., Wiedensohler, A., Young, D.
The NERC UK SOLAS-funded Reactive Halogens in the Marine Boundary Layer (RHaMBLe) programme comprised three field experiments. This manuscript presents an overview of the measurements made within the two simultaneous remote experiments conducted in the tropical North Atlantic in May and June 2007. Measurements were made from two mobile and one ground-based platforms. The heavily instrumented cruise D319 on the RRS Discovery from Lisbon, Portugal to São Vicente, Cape Verde and back to Falmouth, UK was used to characterise the spatial distribution of boundary layer components likely to play a role in reactive halogen chemistry. Measurements onboard the ARSF Dornier aircraft were used to allow the observations to be interpreted in the context of their vertical distribution and to confirm the interpretation of atmospheric structure in the vicinity of the Cape Verde islands. Long-term ground-based measurements at the Cape Verde Atmospheric Observatory (CVAO) on São Vicente were supplemented by long-term measurements of reactive halogen species and characterisation of additional trace gas and aerosol species during the intensive experimental period. This paper presents a summary of the measurements made within the RHaMBLe remote experiments and discusses them in their meteorological and chemical context as determined from these three platforms and from additional meteorological analyses. Air always arrived at the CVAO from the North East with a range of air mass origins (European, Atlantic and North American continental). Trace gases were present at stable and fairly low concentrations with the exception of a slight increase in some anthropogenic components in air of North American origin, though NOx mixing ratios during this period remained below 20 pptv (note the non-IUPAC adoption in this manuscript of pptv and ppbv, equivalent to pmol mol−1 and nmol mol−1 to reflect common practice). Consistency with these air mass classifications is observed in the time series of soluble gas and aerosol composition measurements, with additional identification of periods of slightly elevated dust concentrations consistent with the trajectories passing over the African continent. The CVAO is shown to be broadly representative of the wider North Atlantic marine boundary layer; measurements of NO, O3 and black carbon from the ship are consistent with a clean Northern Hemisphere marine background. Aerosol composition measurements do not indicate elevated organic material associated with clean marine air. Closer to the African coast, black carbon and NO levels start to increase, indicating greater anthropogenic influence. Lower ozone in this region is possibly associated with the increased levels of measured halocarbons, associated with the nutrient rich waters of the Mauritanian upwelling. Bromide and chloride deficits in coarse mode aerosol at both the CVAO and on D319 and the continuous abundance of inorganic gaseous halogen species at CVAO indicate significant reactive cycling of halogens. Aircraft measurements of O3 and CO show that surface measurements are representative of the entire boundary layer in the vicinity both in diurnal variability and absolute levels. Above the inversion layer similar diurnal behaviour in O3 and CO is observed at lower mixing ratios in the air that had originated from south of Cape Verde, possibly from within the ITCZ. ECMWF calculations on two days indicate very different boundary layer depths and aircraft flights over the ship replicate this, giving confidence in the calculated boundary layer depth.
A statistical proxy for sulphuric acid concentration
2011, Mikkonen, S., Romakkaniemi, S., Smith, J.N., Korhonen, H., Petäjä, T., Plass-Duelmer, C., Boy, M., McMurry, P.H., Lehtinen, K.E.J., Joutsensaari, J., Hamed, A., Mauldin III, R.L., Birmili, W., Spindler, G., Arnold, F., Kulmala, M., Laaksonen, A.
Gaseous sulphuric acid is a key precursor for new particle formation in the atmosphere. Previous experimental studies have confirmed a strong correlation between the number concentrations of freshly formed particles and the ambient concentrations of sulphuric acid. This study evaluates a body of experimental gas phase sulphuric acid concentrations, as measured by Chemical Ionization Mass Spectrometry (CIMS) during six intensive measurement campaigns and one long-term observational period. The campaign datasets were measured in Hyytiälä, Finland, in 2003 and 2007, in San Pietro Capofiume, Italy, in 2009, in Melpitz, Germany, in 2008, in Atlanta, Georgia, USA, in 2002, and in Niwot Ridge, Colorado, USA, in 2007. The long term data were obtained in Hohenpeissenberg, Germany, during 1998 to 2000. The measured time series were used to construct proximity measures ("proxies") for sulphuric acid concentration by using statistical analysis methods. The objective of this study is to find a proxy for sulfuric acid that is valid in as many different atmospheric environments as possible. Our most accurate and universal formulation of the sulphuric acid concentration proxy uses global solar radiation, SO2 concentration, condensation sink and relative humidity as predictor variables, yielding a correlation measure (R) of 0.87 between observed concentration and the proxy predictions. Interestingly, the role of the condensation sink in the proxy was only minor, since similarly accurate proxies could be constructed with global solar radiation and SO2 concentration alone. This could be attributed to SO2 being an indicator for anthropogenic pollution, including particulate and gaseous emissions which represent sinks for the OH radical that, in turn, is needed for the formation of sulphuric acid.
Sensitivity of polar stratospheric ozone loss to uncertainties in chemical reaction kinetics
2009, Kawa, S.R., Stolarski, R.S., Newman, P.A., Douglass, A.R., Rex, M., Hofmann, D.J., Santee, M.L., Frieler, K.
The impact and significance of uncertainties in model calculations of stratospheric ozone loss resulting from known uncertainty in chemical kinetics parameters is evaluated in trajectory chemistry simulations for the Antarctic and Arctic polar vortices. The uncertainty in modeled ozone loss is derived from Monte Carlo scenario simulations varying the kinetic (reaction and photolysis rate) parameters within their estimated uncertainty bounds. Simulations of a typical winter/spring Antarctic vortex scenario and Match scenarios in the Arctic produce large uncertainty in ozone loss rates and integrated seasonal loss. The simulations clearly indicate that the dominant source of model uncertainty in polar ozone loss is uncertainty in the Cl2O 2 photolysis reaction, which arises from uncertainty in laboratory-measured molecular cross sections at atmospherically important wavelengths. This estimated uncertainty in JCl 2O2 from laboratory measurements seriously hinders our ability to model polar ozone loss within useful quantitative error limits. Atmospheric observations, however, suggest that the Cl2O2 photolysis uncertainty may be less than that derived from the lab data. Comparisons to Match, South Pole ozonesonde, and Aura Microwave Limb Sounder (MLS) data all show that the nominal recommended rate simulations agree with data within uncertainties when the Cl2O2 photolysis error is reduced by a factor of two, in line with previous in situ ClOx measurements. Comparisons to simulations using recent cross sections from Pope et al. (2007) are outside the constrained error bounds in each case. Other reactions producing significant sensitivity in polar ozone loss include BrO + ClO and its branching ratios. These uncertainties challenge our confidence in modeling polar ozone depletion and projecting future changes in response to changing halogen emissions and climate. Further laboratory, theoretical, and possibly atmospheric studies are needed.