Alternative explanations for rising dissolved organic carbon export from organic soils

Since 1988, there has been, on average, a 91% increase in dissolved organic carbon (DOC) concentrations of UK lakes and streams in the Acid Waters Monitoring Network (AWMN). Similar DOC increases have been observed in surface waters across much of Europe and North America. Much of the debate about the causes of rising DOC has, as in other studies relating to the carbon cycle, focused on factors related to climate change. Data from our peat-core experiments support an influence of climate on DOC, notably an increase in production with temperature under aerobic, and to a lesser extent anaerobic, conditions. However, we argue that climatic factors may not be the dominant drivers of DOC change. DOC solubility is suppressed by high soil water acidity and ionic strength, both of which have decreased as a result of declining sulphur deposition since the 1980s, augmented during the 1990s in the United Kingdom by a cyclical decline in sea-salt deposition. Our observational and experimental data demonstrate a clear, inverse and quantitatively important link between DOC and sulphate concentrations in soil solution. Statistical analysis of 11 AWMN lakes suggests that rising temperature, declining sulphur deposition and changing sea-salt loading can account for the majority of the observed DOC trend. This combination of evidence points to the changing chemical composition of atmospheric deposition, particularly the substantial reduction in anthropogenic sulphur emissions during the last 20 years, as a key cause of rising DOC. The implications of rising DOC export for the carbon cycle will be very different if linked primarily to decreasing acid deposition, rather than to changes in climate, suggesting that these systems may be recovering rather than destabilising.

Buffering of recovery from acidification by organic acids

In the United Kingdom, as in other regions of Europe and North America, recent decreases in surface water sulphate concentrations, due to reduced sulphur emissions, have coincided with marked increases in dissolved organic carbon (DOC) concentrations. Since many of the compounds comprising DOC are acidic, the resulting increases in organic acidity may have the potential to offset the benefits of a decrease in mineral (sulphate) acidity. To test this, we used a triprotic model of organic acid dissociation to estimate the proportional organic acid buffering of reduced mineral acidity as measured in the 22 lakes and streams monitored by the UK Acid Waters Monitoring Network. For an average non-marine sulphate decrease of 30 μeq l− 1 over 15 years from 1988–2003, we estimate that around 28% was counterbalanced by rising strong organic acids, 20% by rising alkalinity (partly attributable to an increase in weak organic acids), 11% by falling inorganic aluminium and 41% by falling non-marine base cations. The situation is complicated by a concurrent decrease in marine ion concentrations, and the impact this may have had on both DOC and acidity, but results clearly demonstrate that organic acid increases have substantially limited the amount of recovery from acidification (in terms of rising alkalinity and falling aluminium) that have resulted from reducing sulphur emissions. The consistency and magnitude of sulphate and organic acid changes are consistent with a causal link between the two, possibly due to the effects of changing acidity, ionic strength and aluminium concentrations on organic matter solubility. If this is the case, then organic acids can be considered effective but partial buffers to acidity change in organic soils, and this mechanism needs to be considered in assessing and modelling recovery from acidification, and in defining realistic reference conditions. However, large spatial variations in the relative magnitude of organic acid and sulphate changes, notably for low-deposition sites in northwestern areas where organic acid increases apparently exceed non-marine sulphate decreases, suggest that additional factors, such as changes in sea-salt deposition and climatic factors, may be required to explain the full magnitude of DOC increases in UK surface waters.

Night-time tropospheric chemistry of the unsaturated alcohols (Z)-pent-2-en-1-ol and pent-1-en-3-ol: Kinetic studies of reactions of NO3 and N2O5 with stress-induced plant emissions

The night-time tropospheric chemistry of two stress-induced volatile organic compounds (VOCs), (Z)-pent-2-en-1-ol and pent-1-en-3-ol, has been studied at room temperature. Rate coefficients for reactions of the nitrate radical (NO3) with these pentenols were measured using the discharge-flow technique. Because of the relatively low volatility of these compounds, we employed off-axis continuous-wave cavity-enhanced absorption spectroscopy for detection of NO3 in order to be able to work in pseudo first-order conditions with the pentenols in large excess over NO3. The rate coefficients were determined to be (1.53 +/- 0.23) x 10(-13) and (1.39 +/- 0.19) x 10(-14) cm(3) molecule(-1) s(-1) for reactions of NO3 with (Z)-pent-2-en-1-ol and pent-1-en-3-ol. An attempt to study the kinetics of these reactions with a relative-rate technique, using N2O5 as source of NO3 resulted in significantly higher apparent rate coefficients. Performing relative-rate experiments in known excesses of NO2 allowed us to determine the rate coefficients for the N2O5 reactions to be (5.0 +/- 2.8) x 10(-19) cm(3) molecule(-1) s(-1) for (Z)-pent-2-en-1-ol, and (9.1 +/- 5.8) x 10(-19) cm(3) molecule(-1) s(-1) for pent-1-en-3-ol. We show that these relatively slow reactions can indeed interfere with rate determinations in conventional relative-rate experiments.

A very rare self-assembled ribbon of fused cyclic water pentamers encapsulated in a Cu-II based metal-organic framework

A metal organic framework of Cu-II, tartarate (tar) and 2,2'-bipyridyl (2,2'-bipy)], {[Cu(tar)(2,2'-bipy)]center dot 5H(2)O}(n)} (1) has been synthesized at the mild ambient condition and characterized by single crystal X-ray crystallography. In the compound, the Cu(2,2'-bipy) entities are bridged by tartarate ions which are coordinated to Cu-II by both hydroxyl and monodentate carboxylate oxygen to form a one-dimensional chain. The non-coordinated water molecules form ID water chains by edge-sharing cyclic water pentamers along with dangling water dimers. It shows reversible water expulsion upon heating. The water chains join the ID coordination polymeric chains to a 31) network through hydrogen-bond interactions.

How will organic carbon stocks in mineral soils evolve under future climate? Global projections using RothC for a range of climate change scenarios

We use a soil carbon (C) model (RothC), driven by a range of climate models for a range of climate scenarios to examine the impacts of future climate on global soil organic carbon (SOC) stocks. The results suggest an overall global increase in SOC stocks by 2100 under all scenarios, but with a different extent of increase among the climate model and emissions scenarios. The impacts of projected land use changes are also simulated, but have relatively minor impacts at the global scale. Whether soils gain or lose SOC depends upon the balance between C inputs and decomposition. Changes in net primary production (NPP) change C inputs to the soil, whilst decomposition usually increases under warmer temperatures, but can also be slowed by decreased soil moisture. Underlying the global trend of increasing SOC under future climate is a complex pattern of regional SOC change. SOC losses are projected to occur in northern latitudes where higher SOC decomposition rates due to higher temperatures are not balanced by increased NPP, whereas in tropical regions, NPP increases override losses due to higher SOC decomposition. The spatial heterogeneity in the response of SOC to changing climate shows how delicately balanced the competing gain and loss processes are, with subtle changes in temperature, moisture, soil type and land use, interacting to determine whether SOC increases or decreases in the future. Our results suggest that we should stop looking for a single answer regarding whether SOC stocks will increase or decrease under future climate, since there is no single answer. Instead, we should focus on improving our prediction of the factors that determine the size and direction of change, and the land management practices that can be implemented to protect and enhance SOC stocks.

Isoprene emissions and climate

Biogenic volatile organic compounds (BVOCs) play an important role in atmospheric chemistry and the carbon cycle. Isoprene is quantitatively the most important of the non-methane BVOCs (NMBVOCs), with an annual emission of about 400–600 TgC; about 90% of this is emitted by terrestrial plants. Incorporating a mechanistic treatment of isoprene emissions within land-surface schemes has recently become a focus for the modelling community, the aim being to quantify the potential magnitude of associated climate feedbacks. However, these efforts are hampered by major uncertainties about why plants emit isoprene and the relative importance of different environmental controls on isoprene emission. The availability and reliability of observations of isoprene fluxes from different types of vegetation is limited, and this also imposes constraints on model development. Nevertheless, progress is being made towards the development of mechanistic models of isoprene emission which, in conjunction with atmospheric chemistry models, will ultimately allow improved quantification of the feedbacks between the terrestrial biosphere and climate under past and future climate states.

Simulated climate change impact on summer dissolved organic carbon release from peat and surface vegetation: implications for drinking water treatment

Uncertainty regarding changes in dissolved organic carbon (DOC) quantity and quality has created interest in managing peatlands for their ecosystem services such as drinking water provision. The evidence base for such interventions is, however, sometimes contradictory. We performed a laboratory climate manipulation using a factorial design on two dominant peatland vegetation types (Calluna vulgaris and Sphagnum Spp.) and a peat soil collected from a drinking water catchment in Exmoor National Park, UK. Temperature and rainfall were set to represent baseline and future conditions under the UKCP09 2080s high emissions scenario for July and August. DOC leachate then underwent standard water treatment of coagulation/flocculation before chlorination. C. vulgaris leached more DOC than Sphagnum Spp. (7.17 versus 3.00 mg g−1) with higher specific ultraviolet (SUVA) values and a greater sensitivity to climate, leaching more DOC under simulated future conditions. The peat soil leached less DOC (0.37 mg g−1) than the vegetation and was less sensitive to climate. Differences in coagulation removal efficiency between the DOC sources appears to be driven by relative solubilisation of protein-like DOC, observed through the fluorescence peak C/T. Post-coagulation only differences between vegetation types were detected for the regulated disinfection by-products (DBPs), suggesting climate change influence at this scale can be removed via coagulation. Our results suggest current biodiversity restoration programmes to encourage Sphagnum Spp. will result in lower DOC concentrations and SUVA values, particularly with warmer and drier summers.

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

Radiative forcing and climate metrics for ozone precursor emissions: the impact of multi-model averaging

Multi-model ensembles are frequently used to assess understanding of the response of ozone and methane lifetime to changes in emissions of ozone precursors such as NOx, VOCs (volatile organic compounds) and CO. When these ozone changes are used to calculate radiative forcing (RF) (and climate metrics such as the global warming potential (GWP) and global temperature-change potential (GTP)) there is a methodological choice, determined partly by the available computing resources, as to whether the mean ozone (and methane) concentration changes are input to the radiation code, or whether each model's ozone and methane changes are used as input, with the average RF computed from the individual model RFs. We use data from the Task Force on Hemispheric Transport of Air Pollution source–receptor global chemical transport model ensemble to assess the impact of this choice for emission changes in four regions (East Asia, Europe, North America and South Asia). We conclude that using the multi-model mean ozone and methane responses is accurate for calculating the mean RF, with differences up to 0.6% for CO, 0.7% for VOCs and 2% for NOx. Differences of up to 60% for NOx 7% for VOCs and 3% for CO are introduced into the 20 year GWP. The differences for the 20 year GTP are smaller than for the GWP for NOx, and similar for the other species. However, estimates of the standard deviation calculated from the ensemble-mean input fields (where the standard deviation at each point on the model grid is added to or subtracted from the mean field) are almost always substantially larger in RF, GWP and GTP metrics than the true standard deviation, and can be larger than the model range for short-lived ozone RF, and for the 20 and 100 year GWP and 100 year GTP. The order of averaging has most impact on the metrics for NOx, as the net values for these quantities is the residual of the sum of terms of opposing signs. For example, the standard deviation for the 20 year GWP is 2–3 times larger using the ensemble-mean fields than using the individual models to calculate the RF. The source of this effect is largely due to the construction of the input ozone fields, which overestimate the true ensemble spread. Hence, while the average of multi-model fields are normally appropriate for calculating mean RF, GWP and GTP, they are not a reliable method for calculating the uncertainty in these fields, and in general overestimate the uncertainty.

Climate responses to anthropogenic emissions of short-lived climate pollutants

Policies to control air quality focus on mitigating emissions of aerosols and their precursors, and other short-lived climate pollutants (SLCPs). On a local scale, these policies will have beneficial impacts on health and crop yields, by reducing particulate matter (PM) and surface ozone concentrations; however, the climate impacts of reducing emissions of SLCPs are less straightforward to predict. In this paper we consider a set of idealised, extreme mitigation strategies, in which the total anthropogenic emissions of individual SLCP emissions species are removed. This provides an upper bound on the potential climate impacts of such air quality strategies. We focus on evaluating the climate responses to changes in anthropogenic emissions of aerosol precursor species: black carbon (BC), organic carbon (OC) and sulphur dioxide (SO2). We perform climate integrations with four fully coupled atmosphere-ocean global climate models (AOGCMs), and examine the effects on global and regional climate of removing the total land-based anthropogenic emissions of each of the three aerosol precursor species. We find that the SO2 emissions reductions lead to the strongest response, with all three models showing an increase in surface temperature focussed in the northern hemisphere high latitudes, and a corresponding increase in global mean precipitation and run-off. Changes in precipitation and run-off patterns are driven mostly by a northward shift in the ITCZ, consistent with the hemispherically asymmetric warming pattern driven by the emissions changes. The BC and OC emissions reductions give a much weaker forcing signal, and there is some disagreement between models in the sign of the climate responses to these perturbations. These differences between models are due largely to natural variability in sea-ice extent, circulation patterns and cloud changes. This large natural variability component to the signal when the ocean circulation and sea-ice are free-running means that the BC and OC mitigation measures do not necessarily lead to a discernible climate response.

Effects of forage source and extruded linseed supplementation on methane emissions from growing dairy cattle of differing body weights

Changes in diet carbohydrate amount and type (i.e., starch vs. fiber) and dietary oil supplements can affect ruminant methane emissions. Our objectives were to measure methane emissions, whole-tract digestibility, and energy and nitrogen utilization from growing dairy cattle at 2 body weight (BW) ranges, fed diets containing either high maize silage (MS) or high grass silage (GS), without or with supplemental oil from extruded linseed (ELS). Four Holstein-Friesian heifers aged 13 mo (BW range from start to finish of 382 to 526 kg) were used in experiment 1, whereas 4 lighter heifers aged 12 mo (BW range from start to finish of 292 to 419 kg) were used in experiment 2. Diets were fed as total mixed rations with forage dry matter (DM) containing high MS or high GS and concentrates in proportions (forage:concentrate, DM basis) of either 75:25 (experiment 1) or 60:40 (experiment 2), respectively. Diets were supplemented without or with ELS (Lintec[AU1: Add manufacturer name and location.]; 260 g of oil/ kg of DM) at 6% of ration DM. Each experiment was a 4 × 4 Latin square design with 33-d periods, with measurements during d 29 to 33 while animals were housed in respiration chambers. Heifers fed MS at a heavier BW (experiment 1) emitted 20% less methane per unit of DM intake (yield) compared with GS (21.4 vs. 26.6, respectively). However, when repeated with heifers of a lower BW (experiment 2), methane yield did not differ between the 2 diets (26.6 g/kg of DM intake). Differences in heifer BW had no overall effect on methane emissions, except when expressed as grams per kilogram of digestible organic matter (OMD) intake (32.4 vs. 36.6, heavy vs. light heifers). Heavier heifers fed MS in experiment 1 had a greater DM intake (9.4 kg/d) and lower OMD (755 g/kg), but no difference in N utilization (31% of N intake) compared with heifers fed GS (7.9 kg/d and 799 g/kg, respectively). Tissue energy retention was nearly double for heifers fed MS compared with GS in experiment 1 (15 vs. 8% of energy intake, respectively). Heifers fed MS in experiment 2 had similar DM intake (7.2 kg/d) and retention of energy (5% of intake energy) and N (28% of N intake), compared with GS-fed heifers, but OMD was lower (741 vs. 765 g/kg, respectively). No effect of ELS was noted on any of the variables measured, irrespective of animal BW, and this was likely due to the relatively low amount of supplemental oil provided. Differences in heifer BW did not markedly influence dietary effects on methane emissions. Differences in methane yield were attributable to differences in dietary starch and fiber composition associated with forage type and source.

Equations to predict methane emissions from cows fed at maintenance energy level in pasture-based systems

Ruminant production is a vital part of food industry but it raises environmental concerns, partly due to the associated methane outputs. Efficient methane mitigation and estimation of emissions from ruminants requires accurate prediction tools. Equations recommended by international organizations or scientific studies have been developed with animals fed conserved forages and concentrates and may be used with caution for grazing cattle. The aim of the current study was to develop prediction equations with animals fed fresh grass in order to be more suitable to pasture-based systems and for animals at lower feeding levels. A study with 25 nonpregnant nonlactating cows fed solely fresh-cut grass at maintenance energy level was performed over two consecutive grazing seasons. Grass of broad feeding quality, due to contrasting harvest dates, maturity, fertilisation and grass varieties, from eight swards was offered. Cows were offered the experimental diets for at least 2 weeks before housed in calorimetric chambers over 3 consecutive days with feed intake measurements and total urine and faeces collections performed daily. Methane emissions were measured over the last 2 days. Prediction models were developed from 100 3-day averaged records. Internal validation of these equations, and those recommended in literature, was performed. The existing in greenhouse gas inventories models under-estimated methane emissions from animals fed fresh-cut grass at maintenance while the new models, using the same predictors, improved prediction accuracy. Error in methane outputs prediction was decreased when grass nutrient, metabolisable energy and digestible organic matter concentrations were added as predictors to equations already containing dry matter or energy intakes, possibly because they explain feed digestibility and the type of energy-supplying nutrients more efficiently. Predictions based on readily available farm-level data, such as liveweight and grass nutrient concentrations were also generated and performed satisfactorily. New models may be recommended for predictions of methane emissions from grazing cattle at maintenance or low feeding levels.

Direct observation by time-resolved infrared spectroscopy of the bright and the dark excited states of the [Ru(phen)2(dppz)]2+ light-switch compound in solution and when bound to DNA

The [Ru(phen)2(dppz)]2+ complex (1) is non-emissive in water but is highly luminescent in organic solvents or when bound to DNA, making it a useful probe for DNA binding. To date, a complete mechanistic explanation for this “light-switch” effect is still lacking. With this in mind we have undertaken an ultrafast time resolved infrared (TRIR) study of 1 and directly observe marker bands between 1280–1450 cm-1, which characterise both the emissive “bright” and the non-emissive “dark” excited states of the complex, in CD3CN and D2O respectively. These characteristic spectral features are present in the [Ru(dppz)3]2+ solvent light-switch complex but absent in [Ru(phen)3]2+, which is luminescent in both solvents. DFT calculations show that the vibrational modes responsible for these characteristic bands are predominantly localised on the dppz ligand. Moreover, they reveal that certain vibrational modes of the “dark” excited state couple with vibrational modes of two coordinating water molecules, and through these to the bulk solvent, thus providing a new insight into the mechanism of the light-switch effect. We also demonstrate that the marker bands for the “bright” state are observed for both L- and D enantiomers of 1 when bound to DNA and that photo-excitation of the complex induces perturbation of the guanine and cytosine carbonyl bands. This perturbation is shown to be stronger for the L enantiomer, demonstrating the different binding site properties of the two enantiomers and the ability of this technique to determine the identity and nature of the binding site of such intercalators.

Vehicular particulate matter emissions in road tunnels in Sao Paulo, Brazil

In the metropolitan area of Sao Paulo, Brazil, ozone and particulate matter ( PM) are the air pollutants that pose the greatest threat to air quality, since the PM and the ozone precursors ( nitrogen oxides and volatile organic compounds) are the main source of air pollution from vehicular emissions. Vehicular emissions can be measured inside road tunnels, and those measurements can provide information about emission factors of in-use vehicles. Emission factors are used to estimate vehicular emissions and are described as the amount of species emitted per vehicle distance driven or per volume of fuel consumed. This study presents emission factor data for fine particles, coarse particles, inhalable particulate matter and black carbon, as well as size distribution data for inhalable particulate matter, as measured in March and May of 2004, respectively, in the Janio Quadros and Maria Maluf road tunnels, both located in Sao Paulo. The Janio Quadros tunnel carries mainly light-duty vehicles, whereas the Maria Maluf tunnel carries light-duty and heavy-duty vehicles. In the Janio Quadros tunnel, the estimated light-duty vehicle emission factors for the trace elements copper and bromine were 261 and 220 mu g km(-1), respectively, and 16, 197, 127 and 92 mg km(-1), respectively, for black carbon, inhalable particulate matter, coarse particles and fine particles. The mean contribution of heavy-duty vehicles to the emissions of black carbon, inhalable particulate matter, coarse particles and fine particles was, respectively 29, 4, 6 and 6 times higher than that of light-duty vehicles. The inhalable particulate matter emission factor for heavy-duty vehicles was 1.2 times higher than that found during dynamometer testing. In general, the particle emissions in Sao Paulo tunnels are higher than those found in other cities of the world.

Ozone formation potentials of volatile organic compounds and ozone sensitivity to their emission in the megacity of Sao Paulo, Brazil

In the present study, a three-dimensional Eulerian photochemical model was employed to estimate the impact that organic compounds have on tropospheric ozone formation in the Metropolitan Area of Sao Paulo (MASP). In the year 2000, base case simulations were conducted in two periods: August 22-24 and March 13-15. Based on the pollutant concentrations calculated by the model, the correlation coefficient relative to observations for ozone ranged from 0.91 to 0.93 in both periods. In the simulations employed to evaluate the ozone potential of individual VOCs, as well as the sensitivity of ozone to the VOC/NO(x) emission ratio, the variation in anthropogenic emissions was estimated at 15% (according to tests performed previously variations of 15% were stable). Although there were significant differences between the two periods, ozone concentrations were found to be much more sensitive to VOCs than to NO(x) in both periods and throughout the study domain. In addition, considering their individual rates of emission from vehicles, the species/classes that were most important for ozone formation were as follows: aromatics with a kOH>2x 10(4) ppm(-1) min(-1); olefins with a kOH 7 x 10(4) ppm(-1) min(-1); olefins with a kOH 7 x 10(4) ppm(-1) min(-1); ethene; and formaldehyde, which are the principal species related to the production, transport, storage and combustion of fossil fuels.