Oxidation of oleic acid at the air-water interface and its potential effects on cloud critical supersaturations

The oxidation of organic films on cloud condensation nuclei has the potential to affect climate and precipitation events. In this work we present a study of the oxidation of a monolayer of deuterated oleic acid (cis-9-octadecenoic acid) at the air-water interface by ozone to determine if oxidation removes the organic film or replaces it with a product film. A range of different aqueous sub-phases were studied. The surface excess of deuterated material was followed by neutron reflection whilst the surface pressure was followed using a Wilhelmy plate. The neutron reflection data reveal that approximately half the organic material remains at the air-water interface following the oxidation of oleic acid by ozone, thus cleavage of the double bond by ozone creates one surface active species and one species that partitions to the bulk (or gas) phase. The most probable products, produced with a yield of similar to(87 +/- 14)%, are nonanoic acid, which remains at the interface, and azelaic acid (nonanedioic acid), which dissolves into the bulk solution. We also report a surface bimolecular rate constant for the reaction between ozone and oleic acid of (7.3 +/- 0.9) x 10(-11) cm(2) molecule s(-1). The rate constant and product yield are not affected by the solution sub-phase. An uptake coefficient of ozone on the oleic acid monolayer of similar to 4 x 10(-6) is estimated from our results. A simple Kohler analysis demonstrates that the oxidation of oleic acid by ozone on an atmospheric aerosol will lower the critical supersaturation needed for cloud droplet formation. We calculate an atmospheric chemical lifetime of oleic acid of 1.3 hours, significantly longer than laboratory studies on pure oleic acid particles suggest, but more consistent with field studies reporting oleic acid present in aged atmospheric aerosol.

Oxidation of biogenic and water-soluble compounds in aqueous and organic aerosol droplets by ozone: a kinetic and product analysis approach using laser Raman tweezers

The results of an experimental study into the oxidative degradation of proxies for atmospheric aerosol are presented. We demonstrate that the laser Raman tweezers method can be used successfully to obtain uptake coeffcients for gaseous oxidants on individual aqueous and organic droplets, whilst the size and composition of the droplets is simultaneously followed. A laser tweezers system was used to trap individual droplets containing an unsaturated organic compound in either an aqueous or organic ( alkane) solvent. The droplet was exposed to gas- phase ozone and the reaction kinetics and products followed using Raman spectroscopy. The reactions of three different organic compounds with ozone were studied: fumarate anions, benzoate anions and alpha pinene. The fumarate and benzoate anions in aqueous solution were used to represent components of humic- like substances, HULIS; a alpha- pinene in an alkane solvent was studied as a proxy for biogenic aerosol. The kinetic analysis shows that for these systems the diffusive transport and mass accommodation of ozone is relatively fast, and that liquid- phase di. ffusion and reaction are the rate determining steps. Uptake coe. ffcients, g, were found to be ( 1.1 +/- 0.7) x 10(-5), ( 1.5 +/- 0.7) x 10 (-5) and ( 3.0 - 7.5) x 10 (-3) for the reactions of ozone with the fumarate, benzoate and a- pinene containing droplets, respectively. Liquid- phase bimolecular rate coe. cients for reactions of dissolved ozone molecules with fumarate, benzoate and a- pinene were also obtained: k(fumarate) = ( 2.7 +/- 2) x 10 (5), k(benzoate) = ( 3.5 +/- 3) x 10 (5) and k(alpha-pinene) = ( 1-3) x 10(7) dm(3) mol (-1) s (- 1). The droplet size was found to remain stable over the course of the oxidation process for the HULIS- proxies and for the oxidation of a- pinene in pentadecane. The study of the alpha- pinene/ ozone system is the first using organic seed particles to show that the hygroscopicity of the particle does not increase dramatically over the course of the oxidation. No products were detected by Raman spectroscopy for the reaction of benzoate ions with ozone. One product peak, consistent with aqueous carbonate anions, was observed when following the oxidation of fumarate ions by ozone. Product peaks observed in the reaction of ozone with alpha- pinene suggest the formation of new species containing carbonyl groups.

Coordination-driven self-assembly of a novel carbonato-bridged heteromolecular neutral nickel(II) triangle by atmospheric CO2 fixation

Formation of a quasi-symmetrical mu(3)-carbonato-bridged self-assembled heteromolecular triangle of Ni(II), [(mu(3)-CO3){Ni-2(salmeNH)(2)(NCS)(2)}[Ni(salmeNH(2))(2)]center dot Et2O center dot H2O (HsalmeNH = 2-[(3-methylamino-propylimino)-methyl]-phenol) involves atmospheric CO2 uptake in a neutral medium, by spontaneous self-reorganization of the starting mononuclear Ni(II)-Schiff-base complex, [Ni(salmeNH)(2)]. The environment around Ni(II) in two of the subunits is different from the third one. The starting complex, (Ni(salmeNH)(2)], and one of the possible intermediate species, [Ni(salmeNH(2))(2)(NCS)(2)], which has a very similar coordination environment to that in the third Ni(II) center, have been characterized structurally. A plausible mechanism for the formation of such a triangle has also been proposed. The compound shows a very strong antiferromagnetic coupling. Fit as a regular triangular arrangement gave J = -53.1, g = 2.24, and R = 1.5 x 10(-4).

Homogeneous and heterogeneous reactions of atmospheric mercury(II) with sulfur(IV)

Atmospheric models suggest that the reduction of Hg(II) to Hg(O) by S(W) prolongs the residence time of mercury. The redox reaction was investigated both in the aqueous phase (where the reductant is sulfite) and on particulate matter (where the reductant in SO2(g)). In both cases, one of the ultimate products is HgS. A mechanism is proposed involving formation of Hg(O) followed by mercury-induced disproportionation of SO2.

Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO+NO2 -> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)(1) and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.: R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data. 2000, 29, 167) 2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a. value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100,. 19398)4 of the bond dissociation energy, D-298(BrO-NO2) = 118 kJ mol(-1), corresponding to Delta H-0(circle) = 114.3 kJ mol(-1) at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing Delta H-0(circle) to an unrealistically high value (149.3 kJ mol(-1)) or by increasing , the average energy transferred in a downward collision, to an unusually large value (> 5000 cm(-1)). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.

Multiple-scattering in radar systems: a review

Although extensively studied within the lidar community, the multiple scattering phenomenon has always been considered a rare curiosity by radar meteorologists. Up to few years ago its appearance has only been associated with two- or three-body-scattering features (e.g. hail flares and mirror images) involving highly reflective surfaces. Recent atmospheric research aimed at better understanding of the water cycle and the role played by clouds and precipitation in affecting the Earth's climate has driven the deployment of high frequency radars in space. Examples are the TRMM 13.5 GHz, the CloudSat 94 GHz, the upcoming EarthCARE 94 GHz, and the GPM dual 13-35 GHz radars. These systems are able to detect the vertical distribution of hydrometeors and thus provide crucial feedbacks for radiation and climate studies. The shift towards higher frequencies increases the sensitivity to hydrometeors, improves the spatial resolution and reduces the size and weight of the radar systems. On the other hand, higher frequency radars are affected by stronger extinction, especially in the presence of large precipitating particles (e.g. raindrops or hail particles), which may eventually drive the signal below the minimum detection threshold. In such circumstances the interpretation of the radar equation via the single scattering approximation may be problematic. Errors will be large when the radiation emitted from the radar after interacting more than once with the medium still contributes substantially to the received power. This is the case if the transport mean-free-path becomes comparable with the instrument footprint (determined by the antenna beam-width and the platform altitude). This situation resembles to what has already been experienced in lidar observations, but with a predominance of wide- versus small-angle scattering events. At millimeter wavelengths, hydrometeors diffuse radiation rather isotropically compared to the visible or near infrared region where scattering is predominantly in the forward direction. A complete understanding of radiation transport modeling and data analysis methods under wide-angle multiple scattering conditions is mandatory for a correct interpretation of echoes observed by space-borne millimeter radars. This paper reviews the status of research in this field. Different numerical techniques currently implemented to account for higher order scattering are reviewed and their weaknesses and strengths highlighted. Examples of simulated radar backscattering profiles are provided with particular emphasis given to situations in which the multiple scattering contributions become comparable or overwhelm the single scattering signal. We show evidences of multiple scattering effects from air-borne and from CloudSat observations, i.e. unique signatures which cannot be explained by single scattering theory. Ideas how to identify and tackle the multiple scattering effects are discussed. Finally perspectives and suggestions for future work are outlined. This work represents a reference-guide for studies focused at modeling the radiation transport and at interpreting data from high frequency space-borne radar systems that probe highly opaque scattering media such as thick ice clouds or precipitating clouds.

The full-spectrum correlated-k method for longwave atmospheric radiative transfer using an effective Planck function

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas and performing a monochromatic radiation calculation for each point. In this presentation it is shown that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K/day due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such that they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide, and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K/day can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K/day for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases.

Twentieth century secular decrease in the atmospheric potential gradient

Current flowing in the global atmospheric electrical circuit (AEC) substantially decreased during the twentieth century. Fair-weather potential gradient (PG) observations in Scotland and Shetland show a previously unreported annual decline from 1920 to 1980, when the measurements ceased. A 25% reduction in PG occurred in Scotland 1920–50, with the maximum decline during the winter months. This is quantitatively explained by a decrease in cosmic rays (CR) increasing the thunderstorm-electrosphere coupling resistance, reducing the ionospheric potential VI. Independent measurements of VI also suggest a reduction of 27% from 1920–50. The secular decrease will influence fair weather atmospheric electrical parameters, including ion concentrations and aerosol electrification. Between 1920–50, the PG showed a negative correlation with global temperature, despite the positive correlation found recently between surface temperature and VI. The 1980s stabilisation in VI may arise from compensation of the continuing CR-induced decline by increases in global temperature and convective electrification.

The global atmospheric electrical circuit and climate

Evidence is emerging for physical links among clouds, global temperatures, the global atmospheric electrical circuit and cosmic ray ionisation. The global circuit extends throughout the atmosphere from the planetary surface to the lower layers of the ionosphere. Cosmic rays are the principal source of atmospheric ions away from the continental boundary layer: the ions formed permit a vertical conduction current to flow in the fair weather part of the global circuit. Through the (inverse) solar modulation of cosmic rays, the resulting columnar ionisation changes may allow the global circuit to convey a solar influence to meteorological phenomena of the lower atmosphere. Electrical effects on non-thunderstorm clouds have been proposed to occur via the ion-assisted formation of ultra-fine aerosol, which can grow to sizes able to act as cloud condensation nuclei, or through the increased ice nucleation capability of charged aerosols. Even small atmospheric electrical modulations on the aerosol size distribution can affect cloud properties and modify the radiative balance of the atmosphere, through changes communicated globally by the atmospheric electrical circuit. Despite a long history of work in related areas of geophysics, the direct and inverse relationships between the global circuit and global climate remain largely quantitatively unexplored. From reviewing atmospheric electrical measurements made over two centuries and possible paleoclimate proxies, global atmospheric electrical circuit variability should be expected on many timescales

Variability in surface atmospheric electric field measurements

An electrical current of the order one picoamp per metre squared flows vertically in the Earth's atmosphere, between the ionosphere at approximately 50km altitude and the surface. This current is generated by global thunderstorm activity and is modulated by galactic cosmic rays and atmospheric aerosol. In fair weather conditions, this current cause a vertical atmospheric electric field, commonly measured as a potential gradient. For circumstances other than fair weather conditions, the potential gradient varies, from small steady enhancements in fog to large fluctuations in thunderstorms. The atmospheric potential gradient is continuously monitored at the Reading University Atmospheric Observatory. An account of the variability of the potential gradient on a variety of time scales will be presented.

Retrieval of global atmospheric electrical activity at a polluted urban site

The atmospheric electrical Potential Gradient (PG) arises from global thunderstorm activity, but surface measurements of the atmospheric Potential Gradient (PG) are influenced by global thunderstorms and local aerosol concentration changes. The local aerosol change can be monitored independently, and in some cases the concentration changes are closely related to PG changes. For these circumstances, a general theory to remove the local aerosol influence on PG measurements has been developed. Continuous measurements of PG and aerosol mass concentration were made during 24–31 Dec, 2005 within an urban environment at Reading, UK. The average diurnal variation of PG showed a double diurnal cycle, with maxima in the early morning and evening hours. The aerosol concentration has similar double maxima. Removing the aerosol using from the PG and aerosol correlation returns a single diurnal cycle, suggestive of the more global PG diurnal cycle.

Carbon fractions associated with silt-size particles in surface and subsurface soil horizons

Information on the distribution and behavior of C fractions in soil particle sizes is crucial for understanding C dynamics in soil. At present little is known about the behavior of the C associated with silt-size particles. We quantified the concentrations, distribution, and enrichment of total C (TC), readily oxidizable C (ROC), hotwater- extractable C (HWC), and cold-water-extractable C (CWC) fractions in coarse (63–20-mm), medium (20–6.3-mm), and fine (6.3–2-mm) silt-size subfractions and in coarse (2000–250 mm) and fine (250–63 mm) sand and clay (<2-mm) soil fractions isolated from bulk soil (<2 mm), and 2- to 4-mm aggregate-size fraction of surface (0–25 cm) and subsurface (25–55 cm) soils under different land uses. All measured C fractions varied significantly across all soil particle-size fractions. The highest C concentrations were associated with the <20-mm soil fractions and peaked in the medium (20–6.3-mm) and fine (6.3–2-mm) silt subfractions in most treatments. Carbon enrichment ratios (ERC) revealed the dual behavior of the C fractions associated with the medium silt-size fraction, demonstrating the simultaneous enrichment of TC and ROC, and the depletion of HWC and CWC fractions. The medium silt (20–6.3-mm) subfraction was identified in this study as a zone where the associated C fractions exhibit transitory qualities. Our results show that investigating subfractions within the silt-size particle fraction provides better understanding of the behavior of C fractions in this soil fraction.

The full-spectrum correlated k method for longwave atmospheric radiative transfer: treatment of gaseous overlap

The correlated k-distribution (CKD) method is widely used in the radiative transfer schemes of atmospheric models, and involves dividing the spectrum into a number of bands and then reordering the gaseous absorption coefficients within each one. The fluxes and heating rates for each band may then be computed by discretizing the reordered spectrum into of order 10 quadrature points per major gas, and performing a pseudo-monochromatic radiation calculation for each point. In this paper it is first argued that for clear-sky longwave calculations, sufficient accuracy for most applications can be achieved without the need for bands: reordering may be performed on the entire longwave spectrum. The resulting full-spectrum correlated k (FSCK) method requires significantly fewer pseudo-monochromatic calculations than standard CKD to achieve a given accuracy. The concept is first demonstrated by comparing with line-by-line calculations for an atmosphere containing only water vapor, in which it is shown that the accuracy of heating-rate calculations improves approximately in proportion to the square of the number of quadrature points. For more than around 20 points, the root-mean-squared error flattens out at around 0.015 K d−1 due to the imperfect rank correlation of absorption spectra at different pressures in the profile. The spectral overlap of m different gases is treated by considering an m-dimensional hypercube where each axis corresponds to the reordered spectrum of one of the gases. This hypercube is then divided up into a number of volumes, each approximated by a single quadrature point, such that the total number of quadrature points is slightly fewer than the sum of the number that would be required to treat each of the gases separately. The gaseous absorptions for each quadrature point are optimized such they minimize a cost function expressing the deviation of the heating rates and fluxes calculated by the FSCK method from line-by-line calculations for a number of training profiles. This approach is validated for atmospheres containing water vapor, carbon dioxide and ozone, in which it is found that in the troposphere and most of the stratosphere, heating-rate errors of less than 0.2 K d−1 can be achieved using a total of 23 quadrature points, decreasing to less than 0.1 K d−1 for 32 quadrature points. It would be relatively straightforward to extend the method to include other gases.