Showing papers in "Atmospheric Chemistry and Physics in 2003"

CO 2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport

Abstract: . Based on about 20 years of NOAA/CMDL's atmospheric CO2 concentration data and a global atmospheric tracer transport model, we estimate interannual variations and spatial patterns of surface CO2 fluxes in the period 01/1982-12/2000, by using a time-dependent Bayesian inversion technique. To increase the reliability of the estimated temporal features, particular care is exerted towards the selection of data records that are homogeneous in time. Fluxes are estimated on a grid-scale resolution (~8o latitude x 10o longitude), constrained by a-priori spatial correlations, and then integrated over different sets of regions. The transport model is driven by interannually varying re-analyzed meteorological fields. We make consistent use of unsmoothed measurements. In agreement with previous studies, land fluxes are estimated to be the main driver of interannual variations in the global CO2 fluxes, with the pace predominantly being set by the El Nino/La Nina contrast. An exception is a 2-3 year period of increased sink of atmospheric carbon after Mt. Pinatubo's volcanic eruption in 1991. The largest differences in fluxes between El Nino and La Nina are found in the tropical land regions, the main share being due to the Amazon basin. The flux variations for the Post-Pinatubo period, the 1997/1998 El Nino, and the 1999 La Nina events are exploited to investigate relations between CO2 fluxes and climate forcing. A rough comparison points to anomalies in precipitation as a prominent climate factor for short-term variability of tropical land fluxes, both through their role on NPP and through promoting fire in case of droughts. Some large flux anomalies seem to be directly related to large biomass burning events recorded by satellite observation. Global ocean carbon uptake shows a trend similar to the one expected if ocean uptake scales proportional to the anthropogenic atmospheric CO2 perturbation. In contrast to temporal variations, the longterm spatial flux distribution can be inferred with lesser robustness only. The tentative pattern estimated by the present inversion exhibits a northern hemisphere land sink on the order of 0.4 PgC/yr (for 01/1996-12/1999, non-fossil fuel carbon only) that is mainly confined to North America. Southern hemisphere land regions are carbon neutral, while the tropical land regions are taking up carbon (e.g., at a rate of 0.8 PgC/yr during 01/1996-12/1999). Ocean fluxes show larger uptake in the Northern mid to high latitudes than in the Southern mid latitude regions, in contrast to the estimates by Takahashi et al. (1999) based on in-situ measurements. On a regional basis, results that differ the most from previous estimates are large carbon uptake of 1 to 1.5 PgC/yr by the Southern temperate Pacific ocean region, weak outgassing from the Southern ocean, and a carbon source from eastern Europe.

A revised parameterization for gaseous dry deposition in air-quality models

Abstract: A parameterization scheme for calculating gaseous dry deposition velocities in air-quality models is revised based on recent study results on non-stomatal uptake of O 3 and SO 2 over 5 different vegetation types. Non-stomatal resistance, which includes in-canopy aerodynamic, soil and cuticle resistances, for SO 2 and O 3 is parameterized as a function of friction velocity, relative humidity, leaf area index, and canopy wetness. Non-stomatal resistance for other chemical species is scaled to those of SO 2 and O 3 based on their chemical and physical characteristics. Stomatal resistance is calculated using a two-big-leaf stomatal resistance sub-model for all gaseous species of interest. The improvements in the present model compared to its earlier version include a newly developed non-stomatal resistance formulation, a realistic treatment of cuticle and ground resistance in winter, and the handling of seasonally-dependent input parameters. Model evaluation shows that the revised parameterization can provide more realistic deposition velocities for both O 3 and SO 2 , especially for wet canopies. Example model output shows that the parameterization provides reasonable estimates of dry deposition velocities for different gaseous species, land types and diurnal and seasonal variations. Maximum deposition velocities from model output are close to reported measurement values for different land types. The current parameterization can be easily adopted into different air-quality models that require inclusion of dry deposition processes.

Multi axis differential optical absorption spectroscopy (MAX-DOAS)

Abstract: . Multi Axis Differential Optical Absorption Spectroscopy (MAX-DOAS) in the atmosphere is a novel measurement technique that represents a significant advance on the well-established zenith scattered sunlight DOAS instruments which are mainly sensitive to stratospheric absorbers. MAX-DOAS utilizes scattered sunlight received from multiple viewing directions. The spatial distribution of various trace gases close to the instrument can be derived by combining several viewing directions. Ground based MAX-DOAS is highly sensitive to absorbers in the lowest few kilometres of the atmosphere and vertical profile information can be retrieved by combining the measurements with Radiative Transfer Model (RTM) calculations. The potential of the technique for a wide variety of studies of tropospheric trace species and its (few) limitations are discussed. A Monte Carlo RTM is applied to calculate Airmass Factors (AMF) for the various viewing geometries of MAX-DOAS. Airmass Factors can be used to quantify the light path length within the absorber layers. The airmass factor dependencies on the viewing direction and the influence of several parameters (trace gas profile, ground albedo, aerosol profile and type, solar zenith and azimuth angles) are investigated. In addition we give a brief description of the instrumental MAX-DOAS systems realised and deployed so far. The results of the RTM studies are compared to several examples of recent MAX-DOAS field experiments and an outlook for future possible applications is given.

A new feedback mechanism linking forests, aerosols, and climate

Abstract: . The possible connections between the carbon balance of ecosystems and aerosol-cloud-climate interactions play a significant role in climate change studies. Carbon dioxide is a greenhouse gas, whereas the net effect of atmospheric aerosols is to cool the climate. Here, we investigated the connection between forest-atmosphere carbon exchange and aerosol dynamics in the continental boundary layer by means of multiannual data sets of particle formation and growth rates, of CO2 fluxes, and of monoterpene concentrations in a Scots pine forest in southern Finland. We suggest a new, interesting link and a potentially important feedback among forest ecosystem functioning, aerosols, and climate: Considering that globally increasing temperatures and CO2 fertilization are likely to lead to increased photosynthesis and forest growth, an increase in forest biomass would increase emissions of non-methane biogenic volatile organic compounds and thereby enhance organic aerosol production. This feedback mechanism couples the climate effect of CO2 with that of aerosols in a novel way.

Weekly cycle of NO 2 by GOME measurements: a signature of anthropogenic sources

Abstract: . Nitrogen oxides (NO+NO 2 =NO x and reservoir species) are important trace gases in the troposphere with impact on human health, atmospheric chemistry and climate. Besides natural sources (lightning, soil emissions) and biomass burning, fossil fuel combustion is estimated to be responsible for about 50% of the total production of NO x . Since human activity in industrialized countries largely follows a seven-day cycle, fossil fuel combustion is expected to be reduced during weekends. This "weekend effect" is well known from local, ground based measurements, but has never been analysed on a global scale before. The Global Ozone Monitoring Experiment (GOME) on board the ESA-satellite ERS-2 allows measurements of NO 2 column densities. By estimating and subtracting the stratospheric column, and considering radiative transfer, vertical column densities (VCD) of tropospheric NO 2 can be determined (e.g. Leue et al., 2001). We demonstrate the statistical analysis of weekly cycles of tropospheric NO 2 VCDs for different regions of the world. In the cycles of the industrialized regions and cities in the US, Europe and Japan a clear Sunday minimum of tropospheric NO 2 VCD can be seen. Sunday NO 2 VCDs are about 25-50% lower than working day levels. Metropolitan areas with other religious and cultural backgrounds (Jerusalem, Mecca) show different weekly patterns corresponding to different days of rest. In China, no weekly pattern can be found. The presence of a weekly cycle in the measured tropospheric NO 2 VCD may help to identify the different anthropogenic source categories. Furthermore, we estimated the lifetime of tropospheric NO 2 by analysing the mean weekly cycle exemplarily over Germany, obtaining a value of about 6 h in summer and 18-24 h in winter.

Global modelling of secondary organic aerosol in the troposphere: a sensitivity analysis

Abstract: . A global 3-dimensional chemistry/transport model able to describe O3, NOx, Volatile Organic Compounds (VOC), sulphur and NH3 chemistry has been extended to simulate the temporal and spatial distribution of primary and secondary carbonaceous aerosols in the troposphere focusing on Secondary Organic Aerosol (SOA) formation. A number of global simulations have been performed to determine a possible range of annual global SOA production and investigate uncertainties associated with the model results. The studied uncertainties in the SOA budget have been evaluated to be in decreasing importance: the potentially irreversible sticking of the semi-volatile compounds on aerosols, the enthalpy of vaporization of these compounds, the partitioning of SOA on non-carbonaceous aerosols, the conversion of aerosols from hydrophobic to hydrophilic, the emissions of primary carbonaceous aerosols, the chemical fate of the first generation products and finally the activity coefficient of the condensable species. The large uncertainties associated with the emissions of VOC and the adopted simplification of chemistry have not been investigated in this study. Although not all sources of uncertainties have been investigated, according to our calculations, the above factors within the experimental range of variations could result to an overall uncertainty of about a factor of 20 in the global SOA budget. The global annual SOA production from biogenic VOC might range from 2.5 to 44.5 Tg of organic matter per year, whereas that from anthropogenic VOC ranges from 0.05 to 2.62 Tg of organic matter per year. These estimates can be considered as a lower limit, since partitioning on coarse particles like nitrate, dust or sea-salt, together with the partitioning and the dissociation of the semi-volatile products in aerosol water has been neglected. Comparison of model results to observations, where available, shows a better agreement for the upper budget estimates than for the lower ones.

Kinetics of the gas-phase reactions of OH radicals with alkanes and cycloalkanes

Abstract: . The available database concerning rate constants for gas-phase reactions of the hydroxyl (OH) radical with alkanes through early 2003 is presented over the entire temperature range for which measurements have been made (~180-2000 K). Measurements made using relative rate methods are re-evaluated using recent rate data for the reference compound (generally recommendations from this review). In general, whenever more than one study has been carried out over an overlapping temperature range, recommended rate constants or temperature-dependent rate expressions are presented. The recommended 298 K rate constants, temperature-dependent parameters, and temperature ranges over which these recommendations are applicable are listed in Table 1.

Urban aerosol number size distributions

Abstract: . Aerosol number size distributions have been measured since 5 May 1997 in Helsinki, Finland. The presented aerosol data represents size distributions within the particle diameter size range 8-400nm during the period from May 1997 to March 2003. The daily, monthly and annual patterns of the aerosol particle number concentrations were investigated. The temporal variation of the particle number concentration showed close correlations with traffic activities. The highest total number concentrations were observed during workdays; especially on Fridays, and the lowest concentrations occurred during weekends; especially Sundays. Seasonally, the highest total number concentrations were observed during winter and spring and lower concentrations were observed during June and July. More than 80% of the number size distributions had three modes: nucleation mode ( 30nm), Aitken mode (20-100nm) and accumulation mode ( }$">90nm). Less than 20% of the number size distributions had either two modes or consisted of more than three modes. Two different measurement sites were used; in the first (Siltavuori, 5.5.1997-5.3.2001), the arithmetic means of the particle number concentrations were 7000cm , 6500cm , and 1000cm respectively for nucleation, Aitken, and accumulation modes. In the second site (Kumpula, 6.3.2001-28.2.2003) they were 5500cm , 4000cm , and 1000cm . The total number concentration in nucleation and Aitken modes were usually significantly higher during workdays than during weekends. The temporal variations in the accumulation mode were less pronounced. The lower concentrations at Kumpula were mainly due to building construction and also the slight overall decreasing trend during these years. During the site changing a period of simultaneous measurements over two weeks were performed showing nice correlation at both sites.

Averaging kernels for DOAS total-column satellite retrievals

Abstract: . The Differential Optical Absorption Spectroscopy (DOAS) method is used extensively to retrieve total column amounts of trace gases based on UV-visible measurements of satellite spectrometers, such as ERS-2 GOME. In practice the sensitivity of the instrument to the tracer density is strongly height dependent, especially in the troposphere. The resulting tracer profile dependence may introduce large systematic errors in the retrieved columns that are difficult to quantify without proper additional information, as provided by the averaging kernel (AK). In this paper we discuss the DOAS retrieval method in the context of the general retrieval theory as developed by Rodgers. An expression is derived for the DOAS AK for optically thin absorbers. It is shown that the comparison with 3D chemistry-transport models and independent profile measurements, based on averaging kernels, is no longer influenced by errors resulting from a priori profile assumptions. The availability of averaging kernel information as part of the total column retrieval product is important for the interpretation of the observations, and for applications like chemical data assimilation and detailed satellite validation studies.

Hygroscopic properties of water-soluble matter and humic-like organics in atmospheric fine aerosol

Abstract: . Ambient continental-rural fine aerosol (K-puszta, Hungary, PM1.5) was sampled on quartz fibre filters in winter and summer 2001. Water-soluble matter (WSM) was extracted in MilliQ-water, and, in a second step, solid phase extraction was used to isolate the less hydrophilic fraction (ISOM) of the water-soluble organic matter (WSOM) from remaining inorganic salts and "most" hydrophilic organic matter (MHOM). This approach allowed ISOM, which constitutes the major fraction of WSOM, to be isolated from ambient aerosols and investigated in pure form. Hygroscopic properties of both WSM and ISOM extracts as well as of aquatic reference fulvic and humic acids were investigated using a Hygroscopicity Tandem Differential Mobility Analyser (H-TDMA). ISOM deliquesced between 30% and 60% relative humidity (RH), and hygroscopic growth factors at 90% RH ranged from 1.08 to 1.17. The hygroscopicity of ISOM is comparable to secondary organic aerosols obtained in smog chamber experiments, but lower than the hygroscopicity of highly soluble organic acids. The hygroscopic behaviour of investigated fulvic and humic acids had similarities to ISOM, but hygroscopic growth factors were slightly smaller and deliquescence was observed at higher RH (75-85% and 85-95% RH for fulvic acid and humic acid, respectively). These differences probably originate from larger average molecular mass and lower solubility of fulvic and humic acids. Inorganic composition data, measured ISOM hygroscopicity, and a presumed value for the hygroscopicity of the small remaining MHOM fraction were used to predict hygroscopic growth of WSM extracts. Good agreement between model prediction and measured water uptake was observed with differences (by volume) ranging from +1% to -18%. While deliquescence properties of WSM extracts were mainly determined by the inorganic salts (42-53 wt % of WSM), the WSOM accounted for a significant fraction of particulate water. At 90% RH, according to model predictions and measurements, about 80-62% of particulate water in the samples are associated with inorganic salts and about 20-38% with WSOM. The relative contributions of both distinguished WSOM fractions, ISOM and MHOM, remains uncertain since MHOM was not available in isolated form, but the results suggest that the less abundant MHOM is also important due to its presumably larger hygroscopicity.

Inorganic bromine in the marine boundary layer: a critical review

Abstract: . The cycling of inorganic bromine in the marine boundary layer (mbl) has received increased attention in recent years. Bromide, a constituent of sea water, is injected into the atmosphere in association with sea-salt aerosol by breaking waves on the ocean surface. Measurements reveal that supermicrometer sea-salt aerosol is substantially depleted in bromine (often exceeding 50%) relative to conservative tracers, whereas marine submicrometer aerosol is often enriched in bromine. Model calculations, laboratory studies, and field observations strongly suggest that the supermicrometer depletions reflect the chemical transformation of particulate bromide to reactive inorganic gases that influence the processing of ozone and other important constituents of marine air. Mechanisms for the submicrometer enrichments are not well understood. Currently available techniques cannot reliably quantify many Br containing compounds at ambient concentrations and, consequently, our understanding of inorganic Br cycling over the oceans and its global significance are uncertain. To provide a more coherent framework for future research, we have reviewed measurements in marine aerosol, the gas phase, and in rain. We also summarize sources and sinks, as well as model and laboratory studies of chemical transformations. The focus is on inorganic bromine over the open oceans outside the polar regions. The generation of sea-salt aerosol at the ocean surface is the major tropospheric source producing about 6.2 Tg/a of bromide. The transport of Br from continents (as mineral aerosol, and as products from biomass-burning and fossil-fuel combustion) can be of local importance. Transport of degradation products of long-lived Br containing compounds from the stratosphere and other sources contribute lesser amounts. Available evidence suggests that, following aerosol acidification, sea-salt bromide reacts to form Br2 and BrCl that volatilize to the gas phase and photolyze in daylight to produce atomic Br and Cl. Subsequent transformations can destroy tropospheric ozone, oxidize dimethylsulfide (DMS) and hydrocarbons in the gas phase and S(IV) in aerosol solutions, and thereby potentially influence climate. The diurnal cycle of gas-phase Br and the corresponding particulate Br deficits are correlated. Higher values of Br in the gas phase during daytime are consistent with expectations based on photochemistry. We expect that the importance of inorganic Br cycling will vary in the future as a function of both increasing acidification of the atmosphere (through anthropogenic emissions) and climate changes. The latter affects bromine cycling via meteorological factors including global wind fields (and the associated production of sea-salt aerosol), temperature, and relative humidity.

Interaction of aerosol particles composed of protein and saltswith water vapor: hygroscopic growth and microstructural rearrangement

Abstract: . The interaction of aerosol particles composed of the protein bovine serum albumin (BSA) and the inorganic salts sodium chloride and ammonium nitrate with water vapor has been investigated by hygroscopicity tandem differential mobility analyzer (H-TDMA) experiments complemented by transmission electron microscopy (TEM) and Kohler theory calculations (100-300nm particle size range, 298K, 960hPa). BSA was chosen as a well-defined model substance for proteins and other macromolecular compounds, which constitute a large fraction of the water-soluble organic component of air particulate matter. Pure BSA particles exhibited deliquescence and efflorescence transitions at 35% relative humidity ( ) and a hygroscopic diameter increase by up to 10% at 95% in good agreement with model calculations based on a simple parameterisation of the osmotic coefficient. Pure NaCl particles were converted from near-cubic to near-spherical shape upon interaction with water vapor at relative humidities below the deliquescence threshold (partial surface dissolution and recrystallisation), and the diameters of pure NH4NO3 particles decreased by up to 10% due to chemical decomposition and evaporation. Mixed NaCl-BSA and NH4NO3-BSA particles interacting with water vapor exhibited mobility equivalent diameter reductions of up to 20%, depending on particle generation, conditioning, size, and chemical composition (BSA dry mass fraction 10-90%). These observations can be explained by formation of porous agglomerates (envelope void fractions up to 50%) due to ion-protein interactions and electric charge effects on the one hand, and by compaction of the agglomerate structure due to capillary condensation effects on the other. The size of NH4NO3-BSA particles was apparently also influenced by volatilisation of NH4NO3, but not as much as for pure salt particles, i.e. the protein inhibited the decomposition of NH4NO3 or the evaporation of the decomposition products NH3 and HNO3. The efflorescence threshold of NaCl-BSA particles decreased with increasing BSA dry mass fraction, i.e. the protein inhibited the formation of salt crystals and enhanced the stability of supersaturated solution droplets. The H-TDMA and TEM results indicate that the protein was enriched at the surface of the mixed particles and formed an envelope, which inhibits the access of water vapor to the particle core and leads to kinetic limitations of hygroscopic growth, phase transitions, and microstructural rearrangement processes. The Kohler theory calculations performed with different types of models demonstrate that the hygroscopic growth of particles composed of inorganic salts and proteins can be efficiently described with a simple volume additivity approach, provided that the correct dry solute mass equivalent diameter and composition are known. A parameterisation for the osmotic coefficient of macromolecular substances has been derived from an osmotic pressure virial equation. For its application only the density and molar mass of the substance have to be known or estimated, and it is fully compatible with traditional volume additivity models for salt mixtures.

Model study of multiphase DMS oxidation with a focus on halogens

Abstract: . We studied the oxidation of dimethylsulfide (DMS) in the marine boundary layer (MBL) with a one-dimensional numerical model and focused on the influence of halogens. Our model runs show that there is still significant uncertainty about the end products of the DMS addition pathway, which is especially caused by uncertainty in the product yield of the reaction of the intermediate product methyl sulfinic acid (MSIA) with OH. BrO strongly increases the importance of the addition branch in the oxidation of DMS even when present at mixing ratios smaller than 0.5pmol mol-1. The inclusion of halogen chemistry leads to higher DMS oxidation rates and smaller DMS to SO2 conversion efficiencies. The DMS to SO2 conversion efficiency is also drastically reduced under cloudy conditions. In cloud-free model runs between 5 and 15% of the oxidized DMS reacts further to particulate sulfur, in cloudy runs this fraction is almost 100%. Sulfate production by HOClaq and HOBraq is important in cloud droplets even for small Br- deficits and related small gas phase halogen concentrations. In general, more particulate sulfur is formed when halogen chemistry is included. A possible enrichment of HCO3- in fresh sea salt aerosol would increase pH values enough to make the reaction of S(IV)* (=SO2,aq+HSO3-+SO32-) with O3 dominant for sulfate production. It leads to a shift from methyl sulfonic acid (MSA) to non-sea salt sulfate (nss-SO42-) production but increases the total nss-SO42- only somewhat because almost all available sulfur is already oxidized to particulate sulfur in the base scenario. We discuss how realistic this is for the MBL. We found the reaction MSAaq+OH to contribute about 10% to the production of nss-SO42- in clouds. It is unimportant for cloud-free model runs. Overall we find that the presence of halogens leads to processes that decrease the albedo of stratiform clouds in the MBL.

The North Atlantic Oscillation controls air pollution transport to the Arctic

Abstract: . This paper studies the interannual variability of pollution pathways from northern hemisphere (NH) continents into the Arctic. Using a 15-year model simulation of the dispersion of passive tracers representative of anthropogenic emissions from NH continents, we show that the North Atlantic Oscillation (NAO) exerts a strong control on the pollution transport into the Arctic, particularly in winter and spring. For tracer lifetimes of 5 (30) days, surface concentrations in the Arctic winter are enhanced by about 70% (30%) during high phases of the NAO (in the following referred to as NAO+) compared to its low phases (NAO-). This is mainly due to great differences in the pathways of European pollution during NAO+ and NAO- phases, respectively, but reinforced by North American pollution, which is also enhanced in the Arctic during NAO+ phases. In contrast, Asian pollution in the Arctic does not significantly depend on the NAO phase. The model results are confirmed using remotely-sensed NO2 vertical atmospheric columns obtained from seven years of satellite measurements, which show enhanced northward NO2 transport and reduced NO2 outflow into the North Atlantic from Central Europe during NAO+ phases. Surface measurements of carbon monoxide (CO) and black carbon at high-latitude stations further corroborate the overall picture of enhanced Arctic pollution levels during NAO+ phases

Global distribution of tropospheric ozone from satellite measurements using the empirically corrected tropospheric ozone residual technique: Identification of the regional aspects of air pollution

Abstract: . Using coincident observations of total ozone from the Total Ozone Mapping Spectrometer (TOMS) and stratospheric ozone profiles from the Solar Backscattered Ultraviolet (SBUV) instruments, detailed maps of tropospheric ozone have been derived on a daily basis over a time period spanning more than two decades. The resultant climatological seasonal depictions of the tropospheric ozone residual (TOR) show much more detail than an earlier analysis that had used coincident TOMS and Stratospheric Aerosol and Gas Experiment (SAGE) ozone profiles, although there are many similarities between the TOMS/SAGE TOR and the TOMS/SBUV TOR climatologies. In particular, both TOR seasonal depictions show large enhancements in the southern tropics and subtropics in austral spring and at northern temperate latitudes during the summer. The much greater detail in this new data set clearly defines the regional aspect of tropospheric ozone pollution in northeastern India, eastern United States, eastern China, and west and southern Africa. Being able to define monthly climatologies for each year of the data record provides enough temporal resolution to illustrate significant interannual variability in some of these regions.

Effects of the physical state of tropospheric ammonium-sulfate-nitrate particles on global aerosol direct radiative forcing

Abstract: . The effect of aqueous versus crystalline sulfate-nitrate-ammonium tropospheric particles on global aerosol direct radiative forcing is assessed. A global three-dimensional chemical transport model predicts sulfate, nitrate, and ammonium aerosol mass. An aerosol thermodynamics model is called twice, once for the upper side (US) and once for lower side (LS) of the hysteresis loop of particle phase. On the LS, the sulfate mass budget is 40% solid ammonium sulfate, 12% letovicite, 11% ammonium bisulfate, and 37% aqueous. The LS nitrate mass budget is 26% solid ammonium nitrate, 7% aqueous, and 67% gas-phase nitric acid release due to increased volatility upon crystallization. The LS ammonium budget is 45% solid ammonium sulfate, 10% letovicite, 6% ammonium bisulfate, 4% ammonium nitrate, 7% ammonia release due to increased volatility, and 28% aqueous. LS aerosol water mass partitions as 22% effloresced to the gas-phase and 78% remaining as aerosol mass. The predicted US/LS global fields of aerosol mass are employed in a Mie scattering model to generate global US/LS aerosol optical properties, including scattering efficiency, single scattering albedo, and asymmetry parameter. Global annual average LS optical depth and mass scattering efficiency are, respectively, 0.023 and 10.7 m2 (g SO4-2)-1, which compare to US values of 0.030 and 13.9 m2 (g SO4-2)-1. Radiative transport is computed, first for a base case having no aerosol and then for the two global fields corresponding to the US and LS of the hysteresis loop. Regional, global, seasonal, and annual averages of top-of-the-atmosphere aerosol radiative forcing on the LS and US (FL and FU, respectively, in W m-2) are calculated. Including both anthropogenic and natural emissions, we obtain global annual averages of FL=-0.750, FU=-0.930, and D FU,L=24% for full sky calculations without clouds and FL=-0.485, FU=-0.605, and D FU,L=25% when clouds are included. Regionally, D FU,L=48% over the USA, 55% over Europe, and 34% over East Asia. Seasonally, D FU,L varies from 18% in DJF to 75% in SON over the USA. The global annual average contribution from anthropogenic aerosol is FL=-0.314 and FU=-0.404, which yield normalized direct radiative forcings (G) of GL=-205 W (g SO4-2)-1 and GU=-264 W (g SO4-2)-1.

An improved infrared carbon monoxide analyser for routine measurements aboard commercial Airbus aircraft: technical validation and first scientific results of the MOZAIC III programme

Abstract: The European-funded MOZAIC programme (Measurements of ozone and water vapour by Airbus in-service aircraft) has been operational since 1994 aboard 5 commercial Airbus A340. It has gathered ozone and water vapour data between the ground and an altitude of 12 km from more than 20 000 long-range flights. A new infrared carbon monoxide analyser has been developed for installation on the MOZAIC equipped aircraft. Improvements in the basic characteristics of a commercial CO analysers have achieved performance suitable for routine aircraft measurements : ±5 ppbv, ±5% precision for a 30 s response time. The first year of operation on board 4 aircraft with more than 900 flights has proven the reliability and the usefulness of this CO analyser. The first scientific results are presented here, including UTLS exchange events and pollution within the boundary layer.

Inverse modeling of CO 2 sources and sinks using satellite data: A synthetic inter-comparison of measurement techniques and their performance as a function of space and time

Abstract: . Currently two polar orbiting satellite instruments measure CO2 concentrations in the Earth's atmosphere, while other missions are planned for the coming years. In the future such instruments might become powerful tools for monitoring changes in the atmospheric CO2 abundance and to improve our quantitative understanding of the leading processes controlling this. At the moment, however, we are still in an exploratory phase where first experiences are collected and promising new space-based measurement concepts are investigated. This study assesses the potential of some of these concepts to improve CO2 source and sink estimates obtained from inverse modelling. For this purpose the performance of existing and planned satellite instruments is quantified by synthetic simulations of their ability to reduce the uncertainty of the current source and sink estimates in comparison with the existing ground-based network of sampling sites. Our high resolution inversion of sources and sinks (at 8°x10°) allows us to investigate the variation of instrument performance in space and time and at various temporal and spatial scales. The results of our synthetic tests clearly indicate that the satellite performance increases with increasing sensitivity of the instrument to CO2 near the Earth's surface, favoring the near infra-red technique. Thermal infrared instruments, on the contrary, reach a better global coverage, because the performance in the near infrared is reduced over the oceans owing to a low surface albedo. Near infra-red sounders can compensate for this by measuring in sun-glint, which will allow accurate measurements over the oceans, at the cost, however, of a lower measurement density. Overall, the sun-glint pointing near infrared instrument is the most promising concept of those tested. We show that the ability of satellite instruments to resolve fluxes at smaller temporal and spatial scales is also related to surface sensitivity. All the satellite instruments performed relatively well over the continents resulting mainly from the larger prior flux uncertainties over land than over the oceans. In addition, the surface networks are rather sparse over land increasing the additional benefit of satellite measurements there. Globally, challenging satellite instrument precisions are needed to compete with the current surface network (about 1ppm for weekly and 8°x10° averaged SCIAMACHY columns). Regionally, however, these requirements relax considerably, increasing to 5ppm for SCIAMACHY over tropical continents. This points not only to an interesting research area using SCIAMACHY data, but also to the fact that satellite requirements should not be quantified by only a single number. The applicability of our synthetic results to real satellite instruments is limited by rather crude representations of instrument and data retrieval related uncertainties. This should receive high priority in future work.

Sensitivities in global scale modeling of isoprene

Abstract: . A sensitivity study of the treatment of isoprene and related parameters in 3D atmospheric models was conducted using the global model of tropospheric chemistry MATCH-MPIC. A total of twelve sensitivity scenarios which can be grouped into four thematic categories were performed. These four categories consist of simulations with different chemical mechanisms, different assumptions concerning the deposition characteristics of intermediate products, assumptions concerning the nitrates from the oxidation of isoprene and variations of the source strengths. The largest differences in ozone compared to the reference simulation occured when a different isoprene oxidation scheme was used (up to 30-60% or about 10 nmol/mol). The largest differences in the abundance of peroxyacetylnitrate (PAN) were found when the isoprene emission strength was reduced by 50% and in tests with increased or decreased efficiency of the deposition of intermediates. The deposition assumptions were also found to have a significant effect on the upper tropospheric HOx production. Different implicit assumptions about the loss of intermediate products were identified as a major reason for the deviations among the tested isoprene oxidation schemes. The total tropospheric burden of O3 calculated in the sensitivity runs is increased compared to the background methane chemistry by 26±9 Tg( O3) from 273 to an average from the sensitivity runs of 299 Tg(O3). % revised Thus, there is a spread of ± 35% of the overall effect of isoprene in the model among the tested scenarios. This range of uncertainty and the much larger local deviations found in the test runs suggest that the treatment of isoprene in global models can only be seen as a first order estimate at present, and points towards specific processes in need of focused future work.

Organic acids as cloud condensation nuclei: Laboratory studies of highly soluble and insoluble species

Abstract: . The ability of sub-micron-sized organic acid particles to act as cloud condensation nuclei (CCN) has been examined at room temperature using a newly constructed continuous-flow, thermal-gradient diffusion chamber (TGDC). The organic acids studied were: oxalic, malonic, glutaric, oleic and stearic. The CCN properties of the highly soluble acids - oxalic, malonic and glutaric - match very closely Kohler theory predictions which assume full dissolution of the dry particle and a surface tension of the growing droplet equal to that of water. In particular, for supersaturations between 0.3 and 0.6, agreement between the dry particle diameter which gives 50% activation and that calculated from Kohler theory is to within 3nm on average. In the course of the experiments, considerable instability of glutaric acid particles was observed as a function of time and there is evidence that they fragment to some degree to smaller particles. Stearic acid and oleic acid, which are both highly insoluble in water, did not activate at supersaturations of 0.6% with dry diameters up to 140nm. Finally, to validate the performance of the TGDC, we present results for the activation of ammonium sulfate particles that demonstrate good agreement with Kohler theory if solution non-ideality is considered. Our findings support earlier studies in the literature that showed highly soluble organics to be CCN active but insoluble species to be largely inactive.

The roles of dynamical variability and aerosols in cirrus cloud formation

Abstract: . The probability of occurrence of ice crystal number densities in young cirrus clouds is examined based on airborne measurements. The observations have been carried out at midlatitudes in both hemispheres at equivalent latitudes (52-55°N/S) during the same season (local autumn in 2000). The in situ measurements considered in the present study include temperatures, vertical velocities, and total ice crystal concentrations, the latter determined with high precision and accuracy using a counterflow virtual impactor. Most young cirrus clouds typically contain high number densities (1-10 cm-3) of small (diameter m m) ice crystals. This mode dominates the probability distributions and is shown to be caused by rapid cooling rates associated with updraft speeds in the range 10-100 cm s-1. A second mode containing larger crystals extends from ~1 cm-3 to low concentrations close to the detection threshold (~3 x 10-4 cm-3) and could be associated with lower updraft speeds. Results of a statistical analysis provide compelling evidence that the dynamical variability of vertical air motions on the mesoscale is the key factor determining the observed probability distributions of pristine ice crystal concentrations in cirrus. Other factors considered are changes of temperature as well as size, number, and ice nucleation thresholds of the freezing aerosol particles. The variability in vertical velocities is caused by atmospheric gravity waves leading to small-scale temperature fluctuations. Inasmuch as gravity waves are widespread, mesoscale variability in vertical velocities can be viewed as a universal feature of young cirrus clouds. Large-scale models that do not account for this subgrid-scale variability yield erroneous predictions of the variability of basic cirrus cloud properties. Climate change may bring about changes in the global distribution of updraft speeds, mean air temperatures, and aerosol properties. As shown in this work, these changes could significantly modify the probability distribution of cirrus ice crystal concentrations. This study emphasizes the key role of vertical velocities and mesoscale variability in vertical velocities in controlling cirrus properties. The results suggest that, in any effort to ascribe cause to trends of cirrus cloud properties, a careful evaluation of dynamical changes in cloud formation should be done before conclusions regarding the role of other anthropogenic factors, such as changes in aerosol composition, are made.

Quantification of topographic venting of boundary layer air to the free troposphere

Abstract: . Net vertical air mass export by thermally driven flows from the atmospheric boundary layer (ABL) to the free troposphere (FT) above deep Alpine valleys was investigated. The vertical export of pollutants above mountainous terrain is presently poorly represented in global chemistry transport models (GCTMs) and needs to be quantified. Air mass budgets were calculated using aircraft observations obtained in deep Alpine valleys. The results show that on average 3 times the valley air mass is exported vertically per day under fair weather conditions. During daytime the type of valleys investigated in this study can act as an efficient "air pump" that transports pollutants upward. The slope wind system within the valley plays an important role in redistributing pollutants. Nitrogen oxide emissions in mountainous regions are efficiently injected into the FT. This could enhance their ozone (O3) production efficiency and thus influences tropospheric pollution budgets. Once lifted to the FT above the Alps pollutants are transported horizontally by the synoptic flow and are subject to European pollution export. Forward trajectory studies show that under fair weather conditions two major pathways for air masses above the Alps dominate. Air masses moving north are mixed throughout the whole tropospheric column and further transported eastward towards Asia. Air masses moving south descend within the subtropical high pressure system above the Mediterranean.

Size-segregated aerosol mass closure and chemical composition in Monte Cimone (I) during MINATROC

Abstract: . Physical and chemical characterizations of the atmospheric aerosol were carried out at Mt. Cimone (Italy) during the 4 June-4 July 2000 period. Particle size distributions in the size range 6nm-10µm were measured with a differential mobility analyzer (DMA) and an optical particle counter (OPC). Size-segregated aerosol was sampled using a 6-stage low pressure impactor. Aerosol samples were submitted to gravimetric and chemical analyses. Ionic, carbonaceous and refractory components of the aerosol were quantified. We compared the sub- and superµm aerosol mass concentrations determined by gravimetric measurements (mGM), chemical analyses (mmCA), and by converting particle size distribution to aerosol mass concentrations (mmSD). Mean random uncertainties associated with the determination of mmGM, mmCA, and mmSD were assessed. The three estimates of the sub-µm aerosol mass concentration agreed, which shows that within experimental uncertainty, the sub-µm aerosol was composed of the quantified components. The three estimates of the super-µm aerosol mass concentration did not agree, which indicates that random uncertainties and/or possible systematic errors in aerosol sampling, sizing or analyses were not adequately accounted for. Aerosol chemical composition in air masses from different origins showed differences, which were significant in regard to experimental uncertainties. During the Saharan dust advection period, coarse dust and fine anthropogenic particles were externally mixed. No anthropogenic sulfate could be found in the super-µm dust particles. In contrast, nitrate was shifted towards the aerosol super-µm fraction in presence of desert dust.

Freezing thresholds and cirrus cloud formation mechanisms inferred from in situ measurements of relative humidity

Abstract: . Factors controlling the microphysical link between distributions of relative humidity above ice saturation in the upper troposphere and lowermost stratosphere and cirrus clouds are examined with the help of microphysical trajectory simulations. Our findings are related to results from aircraft measurements and global model studies. We suggest that the relative humidities at which ice crystals form in the atmosphere can be inferred from in situ measurements of water vapor and temperature close to, but outside of, cirrus clouds. The comparison with concomitant measurements performed inside cirrus clouds provides a clue to freezing mechanisms active in cirrus. The analysis of field data taken at northern and southern midlatitudes in fall 2000 reveals distinct differences in cirrus cloud freezing thresholds. Homogeneous freezing is found to be the most likely mechanism by which cirrus form at southern hemisphere midlatitudes. The results provide evidence for the existence of heterogeneous freezing in cirrus in parts of the polluted northern hemisphere, but do not suggest that cirrus clouds in this region form exclusively on heterogeneous ice nuclei, thereby emphasizing the crucial importance of homogeneous freezing. The key features of distributions of upper tropospheric relative humidity simulated by a global climate model are shown to be in general agreement with both, microphysical simulations and field observations, delineating a feasible method to include and validate ice supersaturation in other large-scale atmospheric models, in particular chemistry-transport and weather forecast models.

One year boundary layer aerosol size distribution data from five nordic background stations

Abstract: Size distribution measurements performed at five different stations have been investigated during a one-year period between 01 June 2000 and 31 May 2001 with focus on diurnal, seasonal and geographical differences of size distribution properties. The stations involved cover a large geographical area ranging from the Finnish Lapland (67o N) down to southern Sweden (56o N) in the order Varrio, Pallas, Hyytiala, Aspvreten and Vavihill. The shape of the size distribution is typically bimodal during winter with a larger fraction of accumulation mode particles compared to the other seasons. Highest Aitken mode concentration is found during summer and spring during the year of study. The maximum of nucleation events occur during the spring months at all stations. Nucleation events occur during other months as well, although not as frequently. Large differences were found between different categories of stations. Northerly located stations such as Pallas and Varrio presented well-separated Aitken and accumulation modes, while the two modes often overlap significantly at the two southernmost stations Vavihill and Aspvreten. A method to cluster trajectories was used to analyse the impact of long-range transport on the observed aerosol properties. Clusters of trajectories arriving from the continent were clearly associated with size distributions shifted towards the accumulation mode. This feature was more pronounced the further south the station was located. Marine- or Arctic-type clusters were associated with large variability in the nuclei size ranges. A quasi-lagrangian approach was used to investigate transport related changes in the aerosol properties. Typically, an increase in especially Aitken mode concentrations was observed when advection from the north occurs, i.e. allowing more continental influence on the aerosol when comparing the different measurement sites. When trajectory clusters arrive to the stations from SW, a gradual decrease in number concentration is experienced in all modes as latitude of measurement site increases.

Laboratory studies of the homogeneous nucleation of iodine oxides

Abstract: . Laboratory experimental results of iodine oxide nucleation are presented. Nucleation was induced following UV photolysis of CF3I or CH2I2 in the presence of excess ozone. Measurements were performed in a 70 L Teflon reactor with new particles detected using an Ultrafine Condensation Particle Counter, UCPC. The experimental results are interpreted using a coupled chemical - aerosol model to derive model parameters assuming single component homogeneous nucleation of OIO. The aerosol model results have been applied in an atmospheric box-model to interpret the possible implications of iodine oxide nucleation in the marine boundary layer. The model calculations demonstrate that IO and OIO concentrations reported in recent field measurements using long path absorption (Allan et al., 2000, 2001) are not sufficient to account for significant aerosol production either in the coastal or open ocean marine boundary layer using the mechanism presented. We demonstrate that inhomogeneous sources of iodine oxides, i.e. "hot" spots with elevated iodine species emissions, could account for the aerosol production bursts observed in the coastal region near Mace Head, Ireland.

Halogen cycling and aerosol pH in the Hawaiian marine boundary layer

Abstract: . Halogen species (HCl* (primarily HCl), Cl* (including Cl2 and HOCl), BrO, total gaseous inorganic Br and size-resolved particulate Cl- and Br -) and related chemical and physical parameters were measured in surface air at Oahu, Hawaii during September 1999. Aerosol pH as a function of particle size was inferred from phase partitioning and thermodynamic properties of HCl. Mixing ratios of halogen compounds and aerosol pHs were simulated with a new version of the photochemical box model MOCCA that considers multiple aerosol size bins. Inferred aerosol pHs ranged from 4.5 to 5.4 (median 5.1, n=22) for super- m m (primarily sea-salt) size fractions and 2.6 to 5.3 (median 4.6) for sub- m m (primarily sulphate) fractions. Inferred daytime pHs tended to be slightly lower than those at night, although daytime median values did not differ statistically from nighttime medians. Simulated pHs for most sea-salt size bins were within the range of inferred values. However, simulated pHs for the largest size fraction in the model were somewhat higher (oscillating around 5.9) due to the rapid turnover rates and relatively larger infusions of sea-salt alkalinity associated with fresh aerosols. Measured mixing ratios of HCl* ranged from Mixing ratios of total gaseous inorganic Br varied from Increasing the sea-salt mixing ratio in MOCCA by ~25% (within observed range) led to a decrease in O3 of ~16%. The chemistry leading to this decrease is complex and is tied to NOx removal by heterogeneous reactions of BrNO3 and ClNO3. The sink of O3 due to the catalytic Cl-ClO and Br-BrO cycles was estimated at -1.0 to -1.5 nmol mol-1 day-1, a range similar to that due to O3 photolysis in the MOCCA simulations.

Climatological aspects of aerosol optical properties in Northern Greece

Abstract: . Measurements of aerosol optical properties (aerosol optical depth, scattering and backscattering coefficients) have been conducted at two ground-based sites in Northern Greece, Ouranoupolis (40° 23' N, 23° 57' E, 170 m a.s.l.) and Thessaloniki (40° 38' N, 22° 57' E, 80 m a.s.l.), between 1999 and 2002. The frequency distributions of the observed parameters have revealed the presence of individual modes of high and low values, indicating the influence from different sources. At both sites, the mean aerosol optical depth at 500 nm was 0.23. Values increase considerably during summer when they remain persistently between 0.3 and 0.5, going up to 0.7-0.8 during specific cases. The mean value of 65±40 Mm-1 of the particle scattering coefficient at 550 nm reflects the impact of continental pollution in the regional boundary layer. Trajectory analysis has shown that higher values of aerosol optical depth and the scattering coefficient are found in the east sector (former Soviet Union countries, eastern Balkan countries), whereas cleaner conditions are found for the NW direction. The influence of Sahara dust events is clearly reflected in the Angstrom exponents. About 45-60% of the observed diurnal variation of the optical properties was attributed to the growth of aerosols with humidity, while the rest of the variability is in phase with the evolution of the sea-breeze cell. The contribution of local pollution is estimated to contribute 35±10% to the average aerosol optical depth at the Thessaloniki site during summer. Finally, the aerosol scale height (aerosol optical depth divided by scattering coefficient) was found to be related to the height of the boundary layer with values between 0.5-1 km during winter and up to 2.5-3 km during summer.

Particle size distribution and particle mass measurements at urban,near-city and rural level in the Copenhagen area and Southern Sweden

Abstract: . Particle size distribution (size-range 3-900nm) and PM10 was measured simultaneously at an urban background station in Copenhagen, a near-city background and a rural location during a period in September-November 2002. The study investigates the contribution from urban versus regional sources of particle number and mass concentration. The total particle number (ToN) and NOx are well correlated at the urban and near-city level and show a distinct diurnal variation, indicating the common traffic source. The average ToN at the three stations differs by a factor of 3. The observed concentrations are 2500#cm , 4500#cm and 7700#cm at rural, near-city and urban level, respectively. PM10 and total particle volume (ToV) are well correlated between the three different stations and show similar concentration levels, in average within 30% relative difference, indicating a common source from long-range transport that dominates the concentrations at all locations. Measures to reduce the local urban emissions of NOx and ToN are likely to affect both the street level and urban background concentrations, while for PM10 and ToV only measurable effects at the street level are probable. Taking into account the supposed stronger health effects of ultrafine particles reduction measures should address particle number emissions. The traffic source contributes strongest in the 10-200nm particle size range. The maximum of the size distribution shifts from about 20-30nm at kerbside to 50-60nm at rural level. Particle formation events were observed in the 3-20nm size range at rural location in the afternoon hours, mainly under conditions with low concentrations of pre-existing aerosol particles. The maximum in the size distribution of the "traffic contribution" seems to be shifted to about 28nm in the urban location compared to 22nm at kerbside. Assuming NOx as an inert tracer on urban scale allows to estimate that ToN at urban level is reduced by 15-30% compared to kerbside. Particle removal processes, e.g. deposition and coagulation, which are most efficient for smallest particle sizes ( 20nm) and condensational growth are likely mechanisms for the loss of particle number and the shift in particle size.

DRIFTS and Knudsen cell study of the heterogeneous reactivity of SO 2 and NO 2 on mineral dust

Abstract: . The heterogeneous oxidation of SO2 by NO2 on mineral dust was studied using Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS) and a Knudsen cell. This made it possible to characterise, kinetically, both the formation of sulfate and nitrate as surface products and the gas phase loss of the reactive species. The gas phase loss rate was determined to be first order in both SO2 and NO2. From the DRIFTS experiment the uptake coefficient, g , for the formation of sulfate was determined to be of the order of 10-10 using the BET area as the reactive surface area. No significant formation of sulfate was seen in the absence of NO2. The Knudsen cell study gave uptake coefficients of the order of 10-6 and 10-7 for SO2 and NO2 respectively. There was no significant difference in uptake when SO2 or NO2 were introduced individually compared to experiments in which SO2 and NO2 were present at the same time.