-
A simple inverse model is proposed to deduce hemisphere-integrated COS flux based on published time series of total column COS. The global atmosphere is divided into two boxes representing the Northern and Southern Hemispheres, and the total column COS data from several stations are used to calculate hemispheric COS loadings. The integrated flux within each hemisphere is calculated as a linear combination of a steady-state solution and time-varying perturbation. The nature of the time-varying perturbation is deduced using two different approaches: an analytic solution based on a cosine function that was fitted to the original total column COS measurement time series and a Simplex optimization with no underlying assumption about the functional form of the total column time series. The results suggest that there is a steady-state COS flux from the Northern to the Southern Hemisphere. There is a seasonal variation superimposed on this flux that in the Southern Hemisphere has a maximum rate of COS input into the atmosphere around January and a maximum rate of COS removal from the atmosphere around August--September. In the Northern Hemisphere, the maximum rate of COS input into the atmosphere is around May--June, and the maximum rate of COS removal is either August or January, depending on which station in the Northern Hemisphere is considered. The results of the inverse model are compared with the outcome of a forward approach on the temporal and spatial variation of the dominant global sources and sinks published earlier. In general, the deduced hemisphere-integrated flux estimates showed good agreement with the database estimates, though it remains uncertain whether COS removal from the atmosphere in the Northern Hemisphere is dominated by plant and soil uptake in the boreal summer or by oceanic uptake in boreal winter.
-
The relationship between nucleation events and spectral solar irradiance was analysed using two years of data collected at the Station for Measuring Forest Ecosystem-Atmosphere Relations (SMEAR II) in Hyytiälä, Finland. We analysed the data in two different ways. In the first step we calculated ten nanometer average values from the irradiance measurements between 280 and 580 nm and explored if any special wavelengths groups showed higher values on event days compared to a spectral reference curve for all the days for 2 years or to reference curves for every month. The results indicated that short wavelength irradiance between 300 and 340 nm is higher on event days in winter (February and March) compared to the monthly reference graph but quantitative much smaller than in spring or summer. By building the ratio between the average values of different event classes and the yearly reference graph we obtained peaks between 1.17 and 1.6 in the short wavelength range (300--340 nm). In the next step we included number concentrations of particles between 3 and 10 nm and calculated correlation coefficients between the different wavelengths groups and the particles. The results were quite similar to those obtained previously; the highest correlation coefficients were reached for the spectral irradiance groups 3--5 (300--330 nm) with average values for the single event classes around 0.6 and a nearly linear decrease towards higher wavelengths groups by 30%. Both analyses indicate quite clearly that short wavelength irradiance between 300 and 330 or 340 nm is the most important solar spectral radiation for the formation of newly formed aerosols. In the end we introduce a photochemical mechanism as one possible pathway how short wavelength irradiance can influence the formation of SOA by calculating the production rate of excited oxygen. This mechanism shows in which way short wavelength irradiance can influence the formation of new particles even though the absolute values are one to two magnitudes smaller compared to irradiance between 400 and 500 nm.
-
Using historical observations and model simulations, we investigate ozone trends prior to the mid-1970s onset of halogen-induced ozone depletion. Though measurements are quite limited, an analysis based on multiple, independent data sets (direct and indirect) provides better constraints than any individual set of observations. We find that three data sets support an apparent long-term stratospheric ozone trend of -7.2 ± 2.3 DU during 1957-1975, which modeling attributes primarily to water vapor increases. The results suggest that 20th century stratospheric ozone depletion may have been roughly 50% more than is generally supposed. Similarly, three data sets support tropospheric ozone increases over polluted Northern Hemisphere continental regions of 8.2 ± 2.1 DU during this period, which are mutually consistent with the stratospheric trends. As with paleoclimate data, which is also based on indirect proxies and/or limited spatial coverage, these results must be interpreted with caution. However, they provide the most thorough estimates presently available of ozone changes prior to the coincident onset of satellite data and halogen dominated ozone changes. If these apparent trends were real, the radiative forcing by stratospheric ozone since the 1950s would then have been -0.15 ± 0.05 W/m<sup>2</sup>, and -0.2 W/m<sup>2</sup> since the preindustrial. For tropospheric ozone, it would have been 0.38 ± 0.10 W/m<sup>2</sup> since the late 1950s. Combined with even a very conservative estimate of tropospheric ozone forcing prior to that time, this would be larger than current estimates since 1850 which are derived from models that are even less well constrained. These calculations demonstrate the importance of gaining a better understanding of historical ozone changes.
-
Biomass burning has long been recognised as an important source of trace gases and aerosols in the atmosphere. The burning of vegetation has a repeating seasonal pattern, but the intensity of burning and the exact localisation of fires vary considerably from year to year. Recent studies have demonstrated the high interannual variability of the emissions that are associated with biomass burning. In this paper I present a methodology using active fire counts from the Along-Track Scanning Radiometer (ATSR) sensor on board the ERS-2 satellite to estimate the seasonal and interannual variability of global biomass burning emissions in the time period 1996--2000. From the ATSR data, I compute relative scaling factors of burning intensity for each month, which are then applied to a standard inventory for carbon monoxide emissions from biomass burning. The new, time-resolved inventory is evaluated using the few existing multi-year burned area observations on continental scales.
-
A tracer model, DREAM (the Danish Rimpuff and Eulerian Accidental release Model), has been developed for modelling transport, dispersion and deposition (wet and dry) of radioactive material from accidental releases, as the Chernobyl accident. The model is a combination of a Lagrangian model, that includes the near source dispersion, and an Eulerian model describing the long-range transport. The performance of the transport model has previously been tested within the European Tracer Experiment, ETEX, which included transport and dispersion of an inert, non-depositing tracer from a controlled release. The focus of this paper is the model performance with respect to the total deposition of <sup>137</sup>Cs, <sup>134</sup>Cs and <sup>131</sup>I from the Chernobyl accident, using different relatively simple and comprehensive parameterizations for dry- and wet deposition. The performance, compared to measurements, of using different combinations of two different wet deposition parameterizations and three different parameterizations of dry deposition has been evaluated, using different statistical tests. The best model performance, compared to measurements, is obtained when parameterizing the total deposition combined of a simple method for dry deposition and a<br> subgrid-scale averaging scheme for wet deposition based on relative humidities. The same major conclusion is obtained for all the three different radioactive isotopes and using two different deposition measurement databases. Large differences are seen in the results obtained by using the two different parameterizations of wet deposition based on precipitation rates and relative humidities, respectively. The parameterization based on subgrid-scale averaging is, in all cases, performing better than the parameterization based on precipitation rates. This indicates that the in-cloud scavenging process is more important than the below cloud scavenging process for the submicron particles and that the precipitation rates are relatively uncertain in the meteorological model compared to the relative humidity. Relatively small differences are, however, seen in the statistical tests between the three different parameterizations of dry deposition.
-
The time resolved chemical composition of aerosol particles, formed by the oxidation of alpha-pinene has been investigated by liquid chromatography/mass spectrometry (LC-MS) using negative and positive ionisation methods (ESI(-) and APCI(+)). The experiments were performed at the EUPHORE facility in Valencia (Spain) under various experimental conditions, including dark ozone reactions, photosmog experiments with low NO<sub>x</sub> mixing ratios and reaction with OH radicals in the absence of NO<sub>x</sub> (H<sub>2</sub>O<sub>2</sub>-photolysis). Particles were sampled on PTFE f ilters at different stages of the reaction and extracted with methanol. The predominant products from alpha-pinene in the particulate phase are cis-pinic acid, cis-pinonic acid and hydroxy-pinonic acid isomers. Another major compound with molecular weight 172 was detected, possibly a hydroxy-carboxylic acid. These major compounds account for 50% to 80% of the identified aerosol products, depending on the time of sampling and type of experiment. In addition, more than 20 different products have been detected and structures have been tentatively assigned based on their molecular weight and responses to the different ionisation modes. The different experiments performed showed that the aerosol formation is mainly caused by the ozonolysis reaction. The highest aerosol yields were observed in the dark ozone experiments, for which also the highest ratios of mass of identified products to the formed aerosol mass were found (30% to 50%, assuming a density of 1 g cm<sup>-3</sup>).
-
In recent years a number of chemistry-climate models have been developed with an emphasis on the stratosphere. Such models cover a wide range of time scales of integration and vary considerably in complexity. The results of specific diagnostics are here analysed to examine the differences amongst individual models and observations, to assess the consistency of model predictions, with a particular focus on polar ozone. For example, many models indicate a significant cold bias in high latitudes, the "cold pole problem", particularly in the southern hemisphere during winter and spring. This is related to wave propagation from the troposphere which can be improved by improving model horizontal resolution and with the use of non-orographic gravity wave drag. As a result of the widely differing modelled polar temperatures, different amounts of polar stratospheric clouds are simulated which in turn result in varying ozone values in the models.<br> <br> The results are also compared to determine the possible future behaviour of ozone, with an emphasis on the polar regions and mid-latitudes. All models predict eventual ozone recovery, but give a range of results concerning its timing and extent. Differences in the simulation of gravity waves and planetary waves as well as model resolution are likely major sources of uncertainty for this issue. In the Antarctic, the ozone hole has probably reached almost its deepest although the vertical and horizontal extent of depletion may increase slightly further over the next few years. According to the model results, Antarctic ozone recovery could begin any year within the range 2001 to 2008.<br> <br> The limited number of models which have been integrated sufficiently far indicate that full recovery of ozone to 1980 levels may not occur in the Antarctic until about the year 2050. For the Arctic, most models indicate that small ozone losses may continue for a few more years and that recovery could begin any year within the range 2004 to 2019. The start of ozone recovery in the Arctic is therefore expected to appear later than in the Antarctic.<br> <br> Further, interannual variability will tend to mask the signal for longer than in the Antarctic, delaying still further the date at which ozone recovery may be said to have started. Because of this inherent variability of the system, the decadal evolution of Arctic ozone will not necessarily be a direct response to external forcing.
-
The millimeter and sub-millimeter waves have been attracting a lot of attention recently in the cloud remote sensing community. This is largely because of their potential use in measuring cirrus cloud parameters with airborne or space-borne radiometers. In this study, we examine the possibility of using polarization measurements in this frequency range to get information on the microphysical properties of cirrus clouds. By using a simple radiative transfer model, we calculated the brightness temperature differences at the vertical and horizontal polarization channels for the following seven frequencies: 90, 157, 220, 340, 463, 683, and 874 GHz. The ice crystals in cirrus clouds are modeled with nearly spherical particles, circular cylinder, and circular plate, as well as with mixtures of these types. We found that the polarization difference signal shows a unique "resonance' feature with the change of ice particle characteristic size: it has a strong response only in a certain range of ice particle size, beyond that range it approaches zero. The size range where this resonance happens depends to a large extent on particle shape and aspect ratio, but to a much less extent on particle orientation. This resonance feature appears even when ice clouds are composed of a mixture of ice crystals in different shapes, although the magnitude and the position of the resonance peak may change, depending on how the mixture is made. Oriented particles generally show larger polarization difference than randomly oriented ones, and plates have larger polarization difference than cylinders. However, the state of particle orientation has a significantly stronger effect on the polarization difference than the particle shape (cylinder or plate). This makes it difficult to distinguish particle shapes using millimeter and sub-millimeter radiometric measurements, if there is no information available on particle orientations. However, if the state of particle shape mixture can be predetermined by other approaches, polarization measurements can help to determine ice particle characteristic size and orientation. This information, in turn, will benefit our retrieval of the ice water path of cirrus clouds.
-
We present a condensed-mass advection based model (MADVEC) designed to simulate the condensation/evaporation of liquid polar stratospheric cloud (PSC) particles. A (Eulerian-in-radius) discretization scheme is used, making the model suitable for use in global or mesoscale chemistry and transport models (CTMs). The mass advection equations are solved using an adaption of the weighted average flux (WAF) scheme. We validate the numerical scheme using an analytical solution for multicomponent aerosols. The physics of the model are tested using a test case designed by Meilinger et al. (1995). The results from this test corroborate the composition gradients across the size distribution under rapid cooling conditions that were reported in earlier studies.
-
The global sulphur cycle has been simulated using a general circulation model with a focus on the source and oxidation of atmospheric dimethylsulphide (DMS). The sensitivity of atmospheric DMS to the oceanic DMS climatology, the parameterisation of the sea-air transfer and to the oxidant fields have been studied. The importance of additional oxidation pathways (by O<sub>3</sub> in the gas- and aqueous-phases and by BrO in the gas phase) not incorporated in global models has also been evaluated. While three different climatologies of the oceanic DMS concentration produce rather similar global DMS fluxes to the atmosphere at 24-27 Tg S yr <sup>-1</sup>, there are large differences in the spatial and seasonal distribution. The relative contributions of OH and NO<sub>3</sub> radicals to DMS oxidation depends critically on which oxidant fields are prescribed in the model. Oxidation by O<sub>3</sub> appears to be significant at high latitudes in both hemispheres. Oxidation by BrO could be significant even for BrO concentrations at sub-pptv levels in the marine boundary layer. The impact of such refinements on the DMS chemistry onto the indirect radiative forcing by anthropogenic sulphate aerosols is also discussed.
