Processes influencing ozone levels in Alaskan forest fire plumes during long-range transport over the North Atlantic

A case of long-range transport of a biomass burning plume from Alaska to Europe is analyzed using a Lagrangian approach. This plume was sampled several times in the free troposphere over North America, the North Atlantic and Europe by three different aircraft during the IGAC Lagrangian 2K4 experiment which was part of the ICARTT/ITOP measurement intensive in summer 2004. Measurements in the plume showed enhanced values of CO, VOCs and NOy, mainly in form of PAN. Observed O3 levels increased by 17 ppbv over 5 days. A photochemical trajectory model, CiTTyCAT, was used to examine processes responsible for the chemical evolution of the plume. The model was initialized with upwind data and compared with downwind measurements. The influence of high aerosol loading on photolysis rates in the plume was investigated using in situ aerosol measurements in the plume and lidar retrievals of optical depth as input into a photolysis code (Fast-J), run in the model. Significant impacts on photochemistry are found with a decrease of 18% in O3 production and 24% in O3 destruction over 5 days when including aerosols. The plume is found to be chemically active with large O3 increases attributed primarily to PAN decomposition during descent of the plume toward Europe. The predicted O3 changes are very dependent on temperature changes during transport and also on water vapor levels in the lower troposphere which can lead to O3 destruction. Simulation of mixing/dilution was necessary to reproduce observed pollutant levels in the plume. Mixing was simulated using background concentrations from measurements in air masses in close proximity to the plume, and mixing timescales (averaging 6.25 days) were derived from CO changes. Observed and simulated O3/CO correlations in the plume were also compared in order to evaluate the photochemistry in the model. Observed slopes change from negative to positive over 5 days. This change, which can be attributed largely to photochemistry, is well reproduced by multiple model runs even if slope values are slightly underestimated suggesting a small underestimation in modeled photochemical O3 production. The possible impact of this biomass burning plume on O3 levels in the European boundary layer was also examined by running the model for a further 5 days and comparing with data collected at surface sites, such as Jungfraujoch, which showed small O3 increases and elevated CO levels. The model predicts significant changes in O3 over the entire 10 day period due to photochemistry but the signal is largely lost because of the effects of dilution. However, measurements in several other BB plumes over Europe show that O3 impact of Alaskan fires can be potentially significant over Europe.

Relationships between lightning and properties of convective cloud clusters

Satellite observations of convective system properties and lightning flash rate are used to investigate the ability of potential lightning parameterizations to capture both the dominant land-ocean contrast in lightning occurrence and regional differences between Africa, the Amazon and the islands of the maritime continent. As found in previous studies, the radar storm height is tightly correlated with the lightning flash rate. A roughly second order power-law fit to the mean radar echo top height above the 0C isotherm is shown to capture both regional and land-ocean contrasts in lightning occurrence and flash rate using a single set of parameters. Recent developments should soon make it possible to implement a parameterization of this kind in global models. Parameterizations based on cloud top height, convective rain rate and convective rain fraction all require the use of separate fits over land and ocean and fail to capture observed differences between continental regions.

A review of the potential impacts of climate change on surface water quality

It is now accepted that some human-induced climate change is unavoidable. Potential impacts on water supply have received much attention, but relatively little is known about the concomitant changes in water quality. Projected changes in air temperature and rainfall could affect river flows and, hence, the mobility and dilution of contaminants. Increased water temperatures will affect chemical reaction kinetics and, combined with deteriorations in quality, freshwater ecological status. With increased flows there will be changes in stream power and, hence, sediment loads with the potential to alter the morphology of rivers and the transfer of sediments to lakes, thereby impacting freshwater habitats in both lake and stream systems. This paper reviews such impacts through the lens of UK surface water quality. Widely accepted climate change scenarios suggest more frequent droughts in summer, as well as flash-flooding, leading to uncontrolled discharges from urban areas to receiving water courses and estuaries. Invasion by alien species is highly likely, as is migration of species within the UK adapting to changing temperatures and flow regimes. Lower flows, reduced velocities and, hence, higher water residence times in rivers and lakes will enhance the potential for toxic algal blooms and reduce dissolved oxygen levels. Upland streams could experience increased dissolved organic carbon and colour levels, requiring action at water treatment plants to prevent toxic by-products entering public water supplies. Storms that terminate drought periods will flush nutrients from urban and rural areas or generate acid pulses in acidified upland catchments. Policy responses to climate change, such as the growth of bio-fuels or emission controls, will further impact freshwater quality.

Trace element partitioning and soil particle characterisation around mining and smelting areas at Tharsis, Riotinto and Huelva, SW Spain

Trace elements may present an environmental hazard in the vicinity of mining and smelting activities. However, the factors controlling trace element distribution in soils around ancient and modem mining and smelting areas are not always clear. Tharsis, Riotinto and Huelva are located in the Iberian Pyrite Belt in SW Spain. Tharsis and Riotinto mines have been exploited since 2500 B.C., with intensive smelting taking place. Huelva, established in 1970 and using the Flash Furnace Outokumpu process, is currently one of the largest smelter in the world. Pyrite and chalcopyrite ore have been intensively smelted for Cu. However, unusually for smelters and mines of a similar size, the elevated trace element concentrations in soils were found to be restricted to the immediate vicinity of the mines and smelters, being found up to a maximum of 2 kin from the mines and smelters at Tharsis, Riotinto and Huelva. Trace element partitioning (over 2/3 of trace elements found in the residual immobile fraction of soils at Tharsis) and soil particles examination by SEM-EDX showed that trace elements were not adsorbed onto soil particles, but were included within the matrix of large trace element-rich Fe silicate slag particles (i.e. 1 min circle divide at least 1 wt.% As, Cu and Zn, and 2 wt.% Pb). Slag particle large size (I mm 0) was found to control the geographically restricted trace element distribution in soils at Tharsis, Riotinto and Huelva, since large heavy particles could not have been transported long distances. Distribution and partitioning indicated that impacts to the environment as a result of mining and smelting should remain minimal in the region. (c) 2006 Elsevier B.V. All rights reserved.

Cross-hemispheric transport of Central African biomass burning pollutants: Implications for downwind ozone production

Pollutant plumes with enhanced concentrations of trace gases and aerosols were observed over the southern coast of West Africa during August 2006 as part of the AMMA wet season field campaign. Plumes were observed both in the mid and upper troposphere. In this study we examined the origin of these pollutant plumes, and their potential to photochemically produce ozone (O3) downwind over the Atlantic Ocean. Their possible contribution to the Atlantic O3 maximum is also discussed. Runs using the BOLAM mesoscale model including biomass burning carbon monoxide (CO) tracers were used to confirm an origin from central African biomass burning fires. The plumes measured in the mid troposphere (MT) had significantly higher pollutant concentrations over West Africa compared to the upper tropospheric (UT) plume. The mesoscale model reproduces these differences and the two different pathways for the plumes at different altitudes: transport to the north-east of the fire region, moist convective uplift and transport to West Africa for the upper tropospheric plume versus north-west transport over the Gulf of Guinea for the mid-tropospheric plume. Lower concentrations in the upper troposphere are mainly due to enhanced mixing during upward transport. Model simulations suggest that MT and UT plumes are 16 and 14 days old respectively when measured over West Africa. The ratio of tracer concentrations at 600 hPa and 250 hPa was estimated for 14–15 August in the region of the observed plumes and compares well with the same ratio derived from observed carbon dioxide (CO2) enhancements in both plumes. It is estimated that, for the period 1–15 August, the ratio of Biomass Burning (BB) tracer concentration transported in the UT to the ones transported in the MT is 0.6 over West Africa and the equatorial South Atlantic. Runs using a photochemical trajectory model, CiTTyCAT, initialized with the observations, were used to estimate in-situ net photochemical O3 production rates in these plumes during transport downwind of West Africa. The mid-troposphere plume spreads over altitude between 1.5 and 6 km over the Atlantic Ocean. Even though the plume was old, it was still very photochemically active (mean net O3 production rates over 10 days of 2.6 ppbv/day and up to 7 ppbv/day during the first days) above 3 km especially during the first few days of transport westward. It is also shown that the impact of high aerosol loads in the MT plume on photolysis rates serves to delay the peak in modelled O3 concentrations. These results suggest that a significant fraction of enhanced O3 in mid-troposphere over the Atlantic comes from BB sources during the summer monsoon period. According to simulated occurrence of such transport, BB may be the main source for O3 enhancement in the equatorial south Atlantic MT, at least in August 2006. The upper tropospheric plume was also still photochemically active, although mean net O3 production rates were slower (1.3 ppbv/day). The results suggest that, whilst the transport of BB pollutants to the UT is variable (as shown by the mesoscale model simulations), pollution from biomass burning can make an important contribution to additional photochemical production of O3 in addition to other important sources such as nitrogen oxides (NOx) from lightning.

HOx observations over West Africa during AMMA: impact of isoprene and NOx

Aircraft OH and HO2 measurements made over West Africa during the AMMA field campaign in summer 2006 have been investigated using a box model constrained to observations of long-lived species and physical parameters. "Good" agreement was found for HO2 (modelled to observed gradient of 1.23 ± 0.11). However, the model significantly overpredicts OH concentrations. The reasons for this are not clear, but may reflect instrumental instabilities affecting the OH measurements. Within the model, HOx concentrations in West Africa are controlled by relatively simple photochemistry, with production dominated by ozone photolysis and reaction of O(1D) with water vapour, and loss processes dominated by HO2 + HO2 and HO2 + RO2. Isoprene chemistry was found to influence forested regions. In contrast to several recent field studies in very low NOx and high isoprene environments, we do not observe any dependence of model success for HO2 on isoprene and attribute this to efficient recycling of HOx through RO2 + NO reactions under the moderate NOx concentrations (5–300 ppt NO in the boundary layer, median 76 ppt) encountered during AMMA. This suggests that some of the problems with understanding the impact of isoprene on atmospheric composition may be limited to the extreme low range of NOx concentrations.

Purification and electron cryomicroscopy of coronavirus particles

Intact, enveloped coronavirus particles vary widely in size and contour, and are thus refractory to study by traditional structural means such as X-ray crystallography. Electron microscopy (EM) overcomes some problems associated with particle variability and has been an important tool for investigating coronavirus ultrastructure. However, EM sample preparation requires that the specimen be dried onto a carbon support film before imaging, collapsing internal particle structure in the case of coronaviruses. Moreover, conventional EM achieves image contrast by immersing the specimen briefly in heavy-metal-containing stain, which reveals some features while obscuring others. Electron cryomicroscopy (cryo-EM) instead employs a porous support film, to which the specimen is adsorbed and flash-frozen. Specimens preserved in vitreous ice over holes in the support film can then be imaged without additional staining. Cryo-EM, coupled with single-particle image analysis techniques, makes it possible to examine the size, structure and arrangement of coronavirus structural components in fully hydrated, native virions. Two virus purification procedures are described.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

Synthesis of hyperbranched poly(aryl amine)s via a carbene insertion approach

A novel diazirine functionalised aniline derivative, 3-(3-aminophenyl)-3-methyldiazirine 1, was prepared and employed as an AB(2)-type monomer in the synthesis of hyperbranched polymers; thus providing the first instance in which polyamines have been prepared via carbene insertion polymerisation. Photolysis of the monomer 1 in bulk and in solution resulted in the formation of hyperbranched poly(aryl amine)s with degrees of polymerisation (DP) varying from 9 to 26 as determined by gel permeation chromatography (GPC). In solution, an increase in the initial monomer concentration was generally found to result in a decrease in the molecular weight characteristics of the resulting poly(aryl amine) s. Subsequent thermal treatment of the poly(aryl amine) s caused a further increase in the DP values up to a maximum of 31. Nuclear magnetic resonance (NMR) spectroscopic analysis revealed that the increase in molecular weight upon thermal treatment resulted from hydroamination of styrenic species formed in the initial photopolymerisation or activation of diazirine moieties.

Ultrafast excited state dynamics controlling photochemical isomerization of N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium coordinated to a {Re-1(CO)(3)(2,2 '-bipyridine)} chromophore

Two Multifunctional photoactive complexes [Re(Cl)(CO)(3)-(MeDpe(+))(2)](2+) and [Re(MeDpe(+))(CO)(3)(bpy)](2+) (MeDpe(+) = N-methyl-4-[trans-2-(4-pyridyl)ethenyl]pyridinium, bpy = 2,2'-bipyridine) were synthesized. characterized. and their redox and photonic properties were investigated by cyclic voltammetry: ultraviolet-visible-infrared (UV/Vis/IR) spectroelectrochemistry, stationary UV/Vis and resonance Raman spectroscopy; photolysis; picosecond time-resolved absorption spectroscopy in the visible and infrared regions: and time-resolved resonance Raman spectroscopy. The first reduction step of either complex Occurs at about -1.1 V versus Fc/Fc(+) and is localized at MeDpe(+). Reduction alone does not induce a trans -> cis isomerization of MeDpe(+). [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) is photostable, while [Re(MeDpe(+))(CO)(3)(bpy)](2+) and free MeDpe(+) isomerize under near-UV irradiation. The lowest excited state of [Re(Cl)(CO)(3)(MeDPe(+))(2)](2+) has been identified as the Re(Cl)(CO)(3) -> MeDpe(+) (MLCT)-M-3 (MLCT = metal-to-ligand charge transfer), decaying directly to the ground state with lifetimes of approximate to 42 (73%) and approximate to 430ps (27%). Optical excitation of [Re(MeDpe(+))(CO)(3)(bpy)](2+) leads to population of Re(CO)(3) -> MeDpe(+) and Re(CO)(3) -> bpy (MLCT)-M-3 states, from which a MeDpe(+) localized intraligand 3 pi pi* excited state ((IL)-I-3) is populated with lifetimes of approximate to 0.6 and approximate to 10 ps, respectively. The 3IL state undergoes a approximate to 21 ps internal rotation, which eventually produces the cis isomer on a much longer timescale. The different excited-state behavior of the two complexes and the absence of thermodynamically favorable interligand electron transfer in excited [Re(MeDpe(+))(CO)(3)(bpy)](2+) reflect the fine energetic balance between excited states of different orbital origin, which can be tuned by subtle Structural variations. The complex [Re(MeDpe+)(CO)(3)(bpy)](2+) emerges as a prototypical, multifunctional species with complementary redox and photonic behavior.

A temperature-dependent relative-rate study of the OH initiated oxidation of n-butane: the kinetics of the reactions of the 1-and 2-butoxy radicals

The kinetics of the reactions of 1-and 2-butoxy radicals have been studied using a slow-flow photochemical reactor with GC-FID detection of reactants and products. Branching ratios between decomposition, CH3CH(O-.)CH2CH3 CH3CHO + C2H5, reaction (7), and reaction with oxygen, CH3CH(O-.)CH2CH3 + O-2 -> CH3C(O)C2H5 + HO2, reaction (6), for the 2-butoxy radical and between isomerization, CH3CH2CH2CH2O. -> CH2CH2CH2CH2OH, reaction (9), and reaction with oxygen, CH3CH2CH2CH2O. + O-2 -> C3H7CHO + HO2, reaction (8), for the 1-butoxy radical were measured as a function of oxygen concentration at atmospheric pressure over the temperature range 250-318 K. Evidence for the formation of a small fraction of chemically activated alkoxy radicals generated from the photolysis of alkyl nitrite precursors and from the exothermic reaction of 2-butyl peroxy radicals with NO was observed. The temperature dependence of the rate constant ratios for a thermalized system is given by k(7)/k(6) = 5.4 x 1026 exp[(-47.4 +/- 2.8 kJ mol(-1))/RT] molecule cm(-3) and k(9)/k(8) = 1.98 x 10(23) exp[(-22.6 +/- 3.9 kJ mol(-1))/RT] molecule cm(-3). The results agree well with the available experimental literature data at ambient temperature but the temperature dependence of the rate constant ratios is weaker than in current recommendations.

Higher oxides of chlorine: absorption cross-sections of Cl2O6 and Cl2O4, the decomposition of Cl2O6, and the reactions of OClO with O and O-3

The absorption cross-sections of Cl2O6 and Cl2O4 have been obtained using a fast flow reactor with a diode array spectrometer (DAS) detection system. The absorption cross-sections at the wavelengths of maximum absorption (lambda(max)) determined in this study are those of Cl2O6: (1.47 +/- 0.15) x 10(-17) cm(2) molecule(-1), at lambda(max) = 276 nm and T = 298 K; and Cl2O4: (9.0 +/- 2.0) x 10(-19) cm(2) molecule(-1), at lambda(max) = 234 nm and T = 298 K. Errors quoted are two standard deviations together with estimates of the systematic error. The shapes of the absorption spectra were obtained over the wavelength range 200-450 nm for Cl2O6 and 200-350 nm for Cl2O4, and were normalized to the absolute cross-sections obtained at lambda(max) for each oxide, and are presented at 1 nm intervals. These data are discussed in relation to previous measurements. The reaction of O with OCIO has been investigated with the objective of observing transient spectroscopic absorptions. A transient absorption was seen, and the possibility is explored of identifying the species with the elusive sym-ClO3 or ClO4, both of which have been characterized in matrices, but not in the gas-phase. The photolysis of OCIO was also re-examined, with emphasis being placed on the products of reaction. UV absorptions attributable to one of the isomers of the ClO dimer, chloryl chloride (ClClO2) were observed; some Cl2O4 was also found at long photolysis times, when much of the ClClO2 had itself been photolysed. We suggest that reports of Cl2O6 formation in previous studies could be a consequence of a mistaken identification. At low temperatures, the photolysis of OCIO leads to the formation of Cl2O3 as a result of the addition of the ClO primary product to OCIO. ClClO2 also appears to be one product of the reaction between O-3 and OCIO, especially when the reaction occurs under explosive conditions. We studied the kinetics of the non-explosive process using a stopped-flow technique, and suggest a value for the room-temperature rate coefficient of (4.6 +/- 0.9) x 10(-19) cm(3) molecule(-1) s(-1) (limit quoted is 2sigma random errors). The photochemical and thermal decomposition of Cl2O6 is described in this paper. For photolysis at k = 254 nm, the removal of Cl2O6 is not accompanied by the build up of any other strong absorber. The implications of the results are either that the photolysis of Cl2O6 produces Cl-2 directly, or that the initial photofragments are converted rapidly to Cl-2. In the thermal decomposition of Cl2O6, Cl2O4 was shown to be a product of reaction, although not necessarily the major one. The kinetics of decomposition were investigated using the stopped-flow technique. At relatively high [OCIO] present in the system, the decay kinetics obeyed a first-order law, with a limiting first-order rate coefficient of 0.002 s(-1). (C) 2004 Elsevier B.V. All rights reserved.

Time-resolved gas-phase kinetic study of the germylene addition reaction, GeH2 + C2D2

Time-resolved studies of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate coefficients for its bimolecular reaction with C2D2. The reaction was studied in the gas phase, mainly at a total pressure of 1.3 kPa (in SF6 bath gas) at five temperatures in the range 298-558 K. Pressure variation measurements over the range 0.13-13 kPa ( SF6) at 298, 397 and 558 K revealed a small pressure dependence but only at 558 K. After correction for this, the second-order rate coefficients gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.96 +/- 0.05) + ( 6.16 +/- 0.37 kJ mol(-1))/RT ln 10 Comparison with the reaction of GeH2 + C2H2 (studied earlier) showed a similar behaviour with almost identical rate coefficients. The lack of a significant isotope effect is consistent with a rate-determining addition process and is explained by irreversible decomposition of the reaction intermediate to give Ge(P-3) + C2H4. This result contrasts with that for GeH2 + C2H4/C2D4 and those for the analogous silylene reactions. The underlying reasons for this are discussed.

Peroxy radical partitioning during the AMMA radical intercomparison exercise

Peroxy radicals were measured onboard two scientific aircrafts during the AMMA (African Monsoon Multidisciplinary Analysis) campaign in summer 2006. This paper reports results from the flight on 16 August 2006 during which measurements of HO2 by laser induced fluorescence spectroscopy at low pressure (LIF-FAGE) and total peroxy radicals (RO2* = HO2+ΣRO2, R = organic chain) by two similar instruments based on the peroxy radical chemical amplification (PeRCA) technique were subject of a blind intercomparison. The German DLR-Falcon and the British FAAM-BAe-146 flew wing tip to wing tip for about 30 min making concurrent measurements on 2 horizontal level runs at 697 and 485 hPa over the same geographical area in Burkina Faso. A full set of supporting measurements comprising photolysis frequencies, and relevant trace gases like CO, NO, NO2, NOy, O3 and a wider range of VOCs were collected simultaneously. Results are discussed on the basis of the characteristics and limitations of the different instruments used. Generally, no data bias are identified and the RO2* data available agree quite reasonably within the instrumental errors. The [RO2*]/[HO2] ratios, which vary between 1:1 and 3:1, as well as the peroxy radical variability, concur with variations in photolysis rates and in other potential radical precursors. Model results provide additional information about dominant radical formation and loss processes.

The north Atlantic marine oundary layer experiment (NAMBLEX): Overview of the campaign held at Mace Head, Ireland, in summer 2002

The North Atlantic Marine Boundary Layer Experiment (NAMBLEX), involving over 50 scientists from 12 institutions, took place at Mace Head, Ireland (53.32° N, 9.90° W), between 23 July and 4 September 2002. A wide range of state-of-the-art instrumentation enabled detailed measurements of the boundary layer structure and atmospheric composition in the gas and aerosol phase to be made, providing one of the most comprehensive in situ studies of the marine boundary layer to date. This overview paper describes the aims of the NAMBLEX project in the context of previous field campaigns in the Marine Boundary Layer (MBL), the overall layout of the site, a summary of the instrumentation deployed, the temporal coverage of the measurement data, and the numerical models used to interpret the field data. Measurements of some trace species were made for the first time during the campaign, which was characterised by predominantly clean air of marine origin, but more polluted air with higher levels of NOx originating from continental regions was also experienced. This paper provides a summary of the meteorological measurements and Planetary Boundary Layer (PBL) structure measurements, presents time series of some of the longer-lived trace species (O3, CO, H2, DMS, CH4, NMHC, NOx, NOy, PAN) and summarises measurements of other species that are described in more detail in other papers within this special issue, namely oxygenated VOCs, HCHO, peroxides, organo-halogenated species, a range of shorter lived halogen species (I2, OIO, IO, BrO), NO3 radicals, photolysis frequencies, the free radicals OH, HO2 and (HO2+Σ RO2), as well as a summary of the aerosol measurements. NAMBLEX was supported by measurements made in the vicinity of Mace Head using the NERC Dornier-228 aircraft. Using ECMWF wind-fields, calculations were made of the air-mass trajectories arriving at Mace Head during NAMBLEX, and were analysed together with both meteorological and trace-gas measurements. In this paper a chemical climatology for the duration of the campaign is presented to interpret the distribution of air-mass origins and emission sources, and to provide a convenient framework of air-mass classification that is used by other papers in this issue for the interpretation of observed variability in levels of trace gases and aerosols.