Items

Advanced search

4D-Var assimilation of MIPAS chemical observations: ozone and nitrogen dioxide analyses

This paper discusses the global analyses of stratospheric ozone (O3) and nitrogen dioxide (NO2) obtained by the Belgian Assimilation System for Chemical Observations from Envisat (BASCOE). Based on a chemistry transport model (CTM) and the 4-dimensional variational (4D-Var) method, BASCOE has assimilated chemical observations of O3, NO2, HNO3, N2O, CH4 and H2O, made between July 2002 and March 2004 by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) onboard the European Space Agency (ESA) Environment Satellite (ENVISAT). This corresponds to the entire period during which MIPAS was operating at its nominal resolution. Our analyses are evaluated against assimilated MIPAS data and independent HALOE (HALogen Occultation Experiment) and POAM-III (Polar Ozone and Aerosol Measurement) satellite data. A good agreement is generally found between the analyses and these datasets, in both cases within the estimated error bars of the observations. The benefit of data assimilation is also evaluated by comparing a BASCOE free model run with MIPAS observations. For O3, the gain from the assimilation is significant during ozone hole conditions, and in the lower stratosphere. Elsewhere, the assimilation does not provide significant improvement. For NO2, the gain from the assimilation is realized through most of the stratosphere. Using the BASCOE analyses, we estimate the differences between MIPAS data and independent data from HALOE and POAM-III, and find results close to those obtained by classical validation methods involving only direct measurement-to-measurement comparisons. Our results extend and reinforce previous MIPAS data validation efforts by taking into account a much larger variety of atmospheric states and measurement conditions. This study discusses possible further developments of the BASCOE data assimilation system; these concern the horizontal resolution, a better filtering of NO2 observations, and the photolysis calculation near the lid of the model. The ozone analyses are part of the PROMOTE project and are publicly available via the BASCOE website (<a href="http://www.bascoe.oma.be/promote/" target="_blank">http://www.bascoe.oma.be/promote/</a>).
Surface and boundary layer exchanges of volatile organic compounds, nitrogen oxides and ozone during the GABRIEL campaign

We present an evaluation of sources, sinks and turbulent transport of nitrogen oxides, ozone and volatile organic compounds (VOC) in the boundary layer over French Guyana and Suriname during the October 2005 GABRIEL campaign by simulating observations with a single-column chemistry and climate model (SCM) along a zonal transect. Simulated concentrations of O3 and NO as well as NO2 photolysis rates over the forest agree well with observations when a small soil-biogenic NO emission flux was applied. This suggests that the photochemical conditions observed during GABRIEL reflect a pristine tropical low-NOx regime. The SCM uses a compensation point approach to simulate nocturnal deposition and daytime emissions of acetone and methanol and produces daytime boundary layer mixing ratios in reasonable agreement with observations. The area average isoprene emission flux, inferred from the observed isoprene mixing ratios and boundary layer height, is about half the flux simulated with commonly applied emission algorithms. The SCM nevertheless simulates too high isoprene mixing ratios, whereas hydroxyl concentrations are strongly underestimated compared to observations, which can at least partly explain the discrepancy. Furthermore, the model substantially overestimates the isoprene oxidation products methlyl vinyl ketone (MVK) and methacrolein (MACR) partly due to a simulated nocturnal increase due to isoprene oxidation. This increase is most prominent in the residual layer whereas in the nocturnal inversion layer we simulate a decrease in MVK and MACR mixing ratios, assuming efficient removal of MVK and MACR. Entrainment of residual layer air masses, which are enhanced in MVK and MACR and other isoprene oxidation products, into the growing boundary layer poses an additional sink for OH which is thus not available for isoprene oxidation. Based on these findings, we suggest pursuing measurements of the tropical residual layer chemistry with a focus on the nocturnal depletion of isoprene and its oxidation products.
Absolute rate constant and O(3P) yield for the O(1D)+N2O reaction in the temperature range 227 K to 719 K

The absolute rate constant for the reaction that is the major source of stratospheric NOx, O(1D)+N2O → products, has been determined in the temperature range 227 K to 719 K, and, in the temperature range 248 K to 600 K, the fraction of the reaction that yields O(3P). Both the rate constants and product yields were determined using a recently-developed chemiluminescence technique for monitoring O(1D) that allows for higher precision determinations for both rate constants, and, particularly, O(3P) yields, than do other methods. We found the rate constant, kR1, to be essentially independent of temperature between 400 K and 227 K, having a value of (1.37±0.11)×10−10 cm3 s−1, and for temperatures greater than 450 K a marked decrease in rate constant was observed, with a rate constant of only (0.94±0.11)×10−10 cm3 s−1 at 719 K. The rate constants determined over the 227 K–400 K range show very low scatter and are significantly greater, by 20% at room temperature and 15% at 227 K, than the current recommended values. The fraction of O(3P) produced in this reaction was determined to be 0.002±0.002 at 250 K rising steadily to 0.010±0.004 at 600 K, thus the channel producing O(3P) can be entirely neglected in atmospheric kinetic modeling calculations. A further result of this study is an expression of the relative quantum yields as a function of temperature for the chemiluminescence reactions (kCL1)C2H + O(1D) → CH(A) + CO and (kCL2)C2H + O(3P) → CH(A) + CO, both followed by CH(A) → CH(X) + hν, as kCL1(T)/kCL2(T)=(32.8T−3050)/(6.29T+398).
Sources of organic carbon in fine particulate matter in northern European urban air

A major fraction of fine particle matter consists of organic carbon (OC) but its origin is still inadequately known. In this study the sources of OC were investigated in the northern European urban environment in Helsinki, Finland. Measurements were carried out over one year and they included both filter (PM1) and online methods. From the filter samples OC, elemental carbon (EC), water-soluble OC (WSOC), levoglucosan and major ions were analyzed. Filter data together with the concentrations of inorganic gases were analyzed by Positive matrix factorization (PMF) in order to find the sources of OC (and WSOC) on an annual as well as on a seasonal basis. In order to study the diurnal variation of sources, OC and EC were measured by a semicontinuous OC/EC analyzer and major ions were determined by a Particle-into-Liquid Sampler coupled to ion chromatographs. According to PMF, OC concentrations were impacted by four sources: biomass combustion, traffic, long-range transport and secondary production. On an annual basis the OC concentration was dominated by secondary organic aerosol (SOA). Its contribution to OC was as high as 64% in summer, which besides anthropogenic sources may also result from the large biogenic volatile organic carbon (VOC) emissions in the boreal region. In winter biomass combustion constituted the largest fraction in OC due to domestic wood combustion for heating purposes. Traffic contributed to OC from 15 to 27%. Regarding the diurnal variation, the contribution from traffic was higher from 08:00 to 18:00 on weekdays than on weekends. The contribution from long-range transport to OC was 24% on average. All four sources also influenced the WSOC concentrations, however, the contribution of SOA was significantly larger for WSOC than OC.
Trans-Pacific dust events observed at Whistler, British Columbia during INTEX-B

The meteorology and physico-chemical characteristics of aerosol associated with two new cases of long range dust transport affecting western Canada during spring 2006 are described. Each event showed enhancements of both sulfate aerosol and crustal material of Asian origin. However, the events were of quite different character and demonstrate the highly variable nature of such events. The April event was a significant dust event with considerable enhancement of fine particle sulfate while the May event was a weaker dust event, also with significant fine particle sulfate enhancement. The latter event was notable in the sense that it was of short duration and was quickly followed by a large increase of organic material likely of regional origin. Comparison of these two events with other documented cases extending back to 1993, suggests that all dust events show coincident enhancements of sulfate and crustal aerosol. However, events vary across a wide continuum based on the magnitude of aerosol enhancements and their sulfate to calcium ratios. At one extreme, events are dominated by highly significant crustal enhancements (e.g. the well-documented 1998 and 2001 "dust" events) while at the other are events with some dust transport, but where sulfate enhancements are of very high magnitude (e.g. the 1993 event at Crater Lake and the 15 May 2006 event at Whistler). Other events represent a "mix". It is likely that this variability is a function of the comparative strengths of the dust and anthropogenic SO2 sources, the transport pathway and in particular the extent to which dust is transported across industrial SO2 sources, and finally, meteorological and chemical processes.
Reconstruction of erythemal UV-doses for two stations in Austria: a comparison between alpine and urban regions

The aim of this study is the reconstruction of past UV-doses for two stations in Austria, Hoher Sonnblick and Vienna, using a physical radiation transfer model. The method uses the modeled UV-dose under clear-sky conditions, cloud modification factors and a correction factor as input variables. To identify the influence of temporal resolution of input data and modification factors, an ensemble of four different modelling approaches was calculated, each with hourly or daily resolution. This is especially important because we found no other study describing the influence of the temporal resolution of input data on model performance. Following the results of the statistical analysis of the evaluation period the model with the highest temporal resolution (HMC) was chosen for the reconstruction of UV-doses. A good agreement between modelled and measured values of erythemally effective UV-doses was found at both stations. In relation to the reference period 1976–1985 an increase in the erythemal UV-dose in Vienna of 11% is visible in the period 1986–1995 and an increase of 17% in the period 1996–2005 can be seen. At Hoher Sonnblick the corresponding increase is 2% and 9%. For the different seasons the strongest increase in erythemal UV-dose has been found for winter and spring season at both stations. Further the influences of total ozone and cloudiness on changes in erythemal UV-doses were analyzed. This analysis showed for both stations, that changes in total ozone had a larger influence on erythemal UV-doses than changes in cloudiness.
Effects of aerosol organics on cloud condensation nucleus (CCN) concentration and first indirect aerosol effect

Aerosol microphysics, chemical composition, and CCN properties were measured on the Department of Energy Gulfstream-1 aircraft during the Marine Stratus/Stratocumulus Experiment (MASE) conducted over the coastal waters between Point Reyes National Seashore and Monterey Bay, California, in July 2005. Aerosols measured during MASE included free tropospheric aerosols, marine boundary layer aerosols, and aerosols with high organic concentration within a thin layer above the cloud. Closure analysis was carried out for all three types of aerosols by comparing the measured CCN concentrations at ~0.2% supersaturation to those predicted based on size distribution and chemical composition using Köhler theory. The effect of aerosol organic species on predicted CCN concentration was examined using a single hygroscopicity parameterization. For aerosols with organics volume fraction up to 70%, such as the marine boundary layer and free troposphere aerosols, CCN concentration and the corresponding first indirect aerosol effect are insensitive to the properties of organics, and can be accurately predicted with a constant hygroscopicity for all organic species. This simplification can facilitate the prediction of indirect aerosol effects using physically-based parameterizations in large scale models. However, for the aerosols within the thin layers above clouds, organics contributed up to 90% of the total aerosol volume, and an accurate knowledge of the overall organic hygroscopicity is required to accurately predict CCN concentrations. Derivations of organic properties in future closure studies, when aerosols are dominated by organic species, would help constrain the descriptions of organics and aerosol-cloud parameterizations in large scale models.
Variations of O3 and CO in summertime at a rural site near Beijing

Large intra-season differences in mixing ratios of CO and O3 were detected at Miyun, a rural site north of Beijing, in summer 2006. Despite an increase in mean daytime mixing ratio of CO from 500 ppbv in June to 700 ppbv in July, mean daytime O3 dropped from 67 ppbv in June to 50 ppbv in July and August. The observed changes in CO and O3 are attributed to the influence of the summer monsoonal circulation that develops over the North China Plain in July. Photochemical production of O3 is reduced as a consequence of increased cloudiness during July and August, as indicated by the strong negative correlation observed between O3 and satellite observations of cloud optical depth, with cloudiness having little effect on CO. The analysis suggests a strategy for emission controls that could be implemented in an economically efficient manner to minimize the frequency of high levels of O3 during summer in Beijing.
Four-dimensional variational data assimilation for inverse modelling of atmospheric methane emissions: method and comparison with synthesis inversion

A four-dimensional variational (4D-Var) data assimilation system for inverse modelling of atmospheric methane emissions is presented. The system is based on the TM5 atmospheric transport model. It can be used for assimilating large volumes of measurements, in particular satellite observations and quasi-continuous in-situ observations, and at the same time it enables the optimization of a large number of model parameters, specifically grid-scale emission rates. Furthermore, the variational method allows to estimate uncertainties in posterior emissions. Here, the system is applied to optimize monthly methane emissions over a 1-year time window on the basis of surface observations from the NOAA-ESRL network. The results are rigorously compared with an analogous inversion by Bergamaschi et al. (2007), which was based on the traditional synthesis approach. The posterior emissions as well as their uncertainties obtained in both inversions show a high degree of consistency. At the same time we illustrate the advantage of 4D-Var in reducing aggregation errors by optimizing emissions at the grid scale of the transport model. The full potential of the assimilation system is exploited in Meirink et al. (2008), who use satellite observations of column-averaged methane mixing ratios to optimize emissions at high spatial resolution, taking advantage of the zooming capability of the TM5 model.
Retrieval of stratospheric aerosol size information from OSIRIS limb scattered sunlight spectra

Recent work has shown that the retrieval of stratospheric aerosol vertical profiles is possible using limb scattered sunlight measurements at optical wavelengths. The aerosol number density profile is retrieved for an assumed particle size distribution and composition. This result can be used to derive the extinction at the measured wavelength. However, large systematic error can result from the uncertainty in the assumed size distribution when the result is used to estimate the extinction at other wavelengths. It is shown in this work that the addition of information obtained from the near infrared limb radiance profile at 1530 nm measured by the imaging module of the OSIRIS instrument yields an indication of the aerosol size distribution profile that can be used to improve the fidelity of the retrievals. A comparison of the estimated extinction profile at 1020 nm with two coincident occultation measurements demonstrates agreement to within approximately 15% from 12 to 27 km altitude.