Impact of the Knudsen number and mass-transfer expression on multi-phase kinetic modeling

Three different mass-transfer expressions are employed within the Model of Aerosol, Gas, and Interfacial Chemistry (MAGIC) to study gas-phase molecular chlorine and bromine production from NaCl and NaBr aerosols, respectively. Simulations of chamber experiments are performed in which NaCl aerosols react with gas-phase ozone in the presence of UV light, in order to identify the importance of the Knudsen number and mass-transfer expression in systems with varying contributions from gas-phase, aqueous-phase, and interfacial chemistry. In the case of NaBr aerosols, simulations are performed of both dark and photolytic conditions. A range of Knudsen numbers spanning the continuum, transition and free-molecular regimes is studied. Particle size is varied over three orders of magnitude, and particle concentration is changed to keep either (a) total aerosol volume or (b) total aerosol surface area constant. When total aerosol volume is constant, the total amount of surface area available for interfacial reaction increases linearly with Knudsen number. Consequently peak gas-phase Cl₂ and Br₂ concentrations increase by two orders of magnitude from the continuum regime to the free-molecular regime. When total aerosol surface area is constant, total aerosol volume is inversely proportional to Knudsen number, with lesser volume being available at higher Knudsen numbers. Consequently Cl^? depletion in the kinetic regime leads to most gas-phase Cl₂ being produced in the transition regime. Gas-phase Br₂ concentration trends are determined by aqueous-phase reaction mechanisms, leading to a monotonic decrease in production with Knudsen number. At all Knudsen numbers, more gas-phase bromine is produced in the photolytic case than in the dark case, the difference being significant in the transition regime. Results of this study suggest that halogen production is insensitive to the mass-transfer expression used in the simulations.     

Characteristics of springtime profiles and sources of ozone in the low troposphere over northern Taiwan

To quantify the possible sources of the high ambient ozone concentration in the low troposphere over Taiwan, ozone sounding data from a two-year intensive field measurement program conducted in April and early May of 2004 and 2005 in northern Taiwan has been examined. We found that the vertical ozone distributions and occurrence of enhanced ozone in the lower troposphere (below 6 km) mainly resulted from (1)Type NE: the long-range transport of ozone controlled by the prevailing northeasterly winds below 2 km, (2)Type LO: the local photochemical ozone production process, and (3)Type SW: the strong southwest/westerly winds aloft (2–6 km). In the boundary layer (BL), where Asian continental outflow prevails, the average profile for type NE is characterized by a peak ozone concentration of nearly 65 ppb at about 1500 m altitude. For type LO, high ozone concentration with an average ozone concentration greater than 80 ppb was also found in the BL in the case of stagnant atmospheric and sunny weather conditions dominated. For type SW, significant ozone enhancement with average ozone concentration of 70–85 ppb was found at around 4 km altitude. It is about 10 ppb greater than that of the types NE and LO at the same troposphere layer owing to the contribution of the biomass burning over Indochina. Due to Taiwan's unique geographic location, the complex interaction of these ozone features in the BL and aloft, especially features associated with northeasterly and south/southwesterly winds, have resulted in complex characteristics of ozone distributions in the lower troposphere over northern Taiwan.     

Aerosol speciation and mass prediction from toluene oxidation under high NOx conditions

A kinetically based gas-particle partitioning box model is used to highlight the importance of parameter representation in the prediction of secondary organic aerosol (SOA) formation following the photo-oxidation of toluene. The model is initialized using experimental data from York University's indoor smog chamber and provides a prediction of the total aerosol yield and speciation. A series of model sensitivity experiments were performed to study the aerosol speciation and mass prediction under high NO_x conditions (VOC/NO_x = 0.2). Sensitivity experiments indicate vapour pressure estimation to be a large area of weakness in predicting aerosol mass, creating an average total error range of 70 μg m^?3 (range of 5–145 μg m^?3), using two different estimation methods. Aerosol speciation proved relatively insensitive to changes in vapour pressure. One species, 3-methyl-6-nitro-catechol, dominated the aerosol phase regardless of the vapour pressure parameterization used and comprised 73–88% of the aerosol by mass. The dominance is associated with the large concentration of 3-methyl-6-nitro-catechol in the gas-phase. The high NO_x initial conditions of this study suggests that the predominance of 3-methyl-6-nitro-catechol likely results from the cresol-forming branch in the Master Chemical Mechanism taking a significant role in secondary organic aerosol formation under high NO_x conditions. Further research into the yields and speciation leading to this reaction product is recommended.     

Global scale emission and distribution of sea-spray aerosol: Sea-salt and organic enrichment

The chemical composition of marine aerosols as a function of their size is an important parameter for the evaluation of their impact on the global climate system. In this work we model fine particle organic matter emitted by sea spray processes and its influence on the aerosol chemical properties at the global scale using the off-line global Chemistry-Transport Model TM5. TM5 is coupled to a microphysical aerosol dynamics model providing size resolved information on particle masses and numbers. The mass of the emitted sea spray particles is partitioned between water insoluble organic matter (WIOM) and sea salt components in the accumulation mode using a function that relates the emitted organic fraction to the surface ocean chlorophyll-a concentrations. The global emission in the sub-micron size range of organic matter by sea spray process is 8.2 Tg yr^?1, compared to 24 Tg fine yr^?1 sea-salt emissions. When the marine sources are included, the concentrations of modelled primary particulate organic matter (POM) increase mainly over the oceans. The model predictions of WIOM and sea salt are evaluated against measurements carried out at Mace Head (Northern Hemisphere) and Amsterdam Island (Southern Hemisphere), showing that in clean marine conditions WIOM marine emissions contribute significantly to POM values.     

Comparison of land-use regression models between Great Britain and the Netherlands

Land-use regression models have increasingly been applied for air pollution mapping at typically the city level. Though models generally predict spatial variability well, the structure of models differs widely between studies. The observed differences in the models may be due to artefacts of data and methodology or underlying differences in source or dispersion characteristics. If the former, more standardised methods using common data sets could be beneficial. We compared land-use regression models for NO₂ and PM_10, developed with a consistent protocol in Great Britain (GB) and the Netherlands (NL).Models were constructed on the basis of 2001 annual mean concentrations from the national air quality networks. Predictor variables used for modelling related to traffic, population, land use and topography. Four sets of models were developed for each country. First, predictor variables derived from data sets common to both countries were used in a pooled analysis, including an indicator for country and interaction terms between country and the identified predictor variables. Second, the common data sets were used to develop individual baseline models for each country. Third, the country-specific baseline models were applied after calibration in the other country to explore transferability. The fourth model was developed using the best possible predictor variables for each country.A common model for GB and NL explained NO₂ concentrations well (adjusted R² 0.64), with no significant differences in intercept and slopes between the two countries. The country-specific model developed on common variables for NL but not GB improved the prediction.The performance of models based upon common data was only slightly worse than models optimised with local data. Models transferred to the other country performed substantially worse than the country-specific models. In conclusion, care is needed both in transferring models across different study areas, and in developing large inter-regional LUR models.     

Establishing ozone flux–response relationships for winter wheat: Analysis of uncertainties based on data for UK and Polish genotypes

The work outlined in this paper had three objectives. The first was to explore the effects of ozone pollution on grain yield and quality of commercially-grown winter wheat cultivars. The second was to derive a stomatal ozone flux model for winter wheat and compare with those already developed for spring wheat. The third was to evaluate exposure- versus flux–response approaches from a risk assessment perspective, and explore the implications of genetic variation in modelled ozone flux.Fifteen winter wheat cultivars were grown in open-top chambers where they were exposed to four levels of ozone. During fumigation, stomatal conductance measurements were made over the lifespan of the flag leaf across a range of environmental conditions. Although significant intra-specific variation in ‘ozone sensitivity’ (in terms of impacts on yield) was identified, yield was inversely related (R² = 0.63, P < 0.001) to the accumulated hourly averaged ozone exposure above 40 ppb during daylight hours (AOT40) across the dataset. The adverse effect of ozone on yield was principally due to a decline in seed weight. Algorithms defining the influence of environmental variables on stomatal uptake were subtly different from those currently in use, based on data for spring wheat, to map ozone impacts on pan-European cereal yield. Considerable intra-specific variation in phenological effects was identified. This meant that an ‘average behaviour’ had to be derived which reduced the predictive capability of the derived stomatal flux model (R² = 0.49, P < 0.001, 15 cultivars included). Indeed, given the intra-specific variability encountered, the flux model that was derived from the full dataset was no better in predicting O₃ impacts on wheat yield than was the AOT40 index. The study highlights the need to use ozone risk assessment tools appropriate to specific vegetation types when modelling and mapping ozone impacts at the regional level.     

Characterization of particulate matter concentrations during controlled indoor activities

Indoor sources have been identified as a major contributor to the increase of particle concentration in indoor environments. The work presented here is a study of the characteristics of particulate matter number size distribution and mass concentration under controlled indoor activities in a laboratory room. The objective is to characterize particulate matter concentrations indoors resulted under the influence of specific sources. Measurements were performed in an empty laboratory (period September–October 2006) using a GRIMM SMPS+C system (particle size range between 11.1 and 1083.3 nm), a DustTrak Aerosol Monitor (TSI) and a P-Trak Ultrafine Particle Counter (TSI). The studied indoor activities included candle burning, hot plate heating, water boiling, onion frying, vacuuming, hair drying, hair spraying, smoking and burning of incense stick. The AMANpsd computer algorithm was used to evaluate the modal structure of measured particle number size distribution data. Furthermore, the change of the particle number size distribution shape under the influence of different emission sources was studied versus time. Finally the particle emission rates were computed. High particle number concentrations were observed during smoking, onion frying, candle burning and incense stick burning. The highest particle mass concentrations were measured during smoking and hair spraying. The shift of the particle size distribution to larger diameters suggests the presence of strong coagulation effect during candle burning, incense stick burning, smoking and onion frying. The size distribution was mainly bimodal during onion frying and candle burning, whereas the size distribution remained unimodal during incense stick burning and smoking experiments.     

Improved land cover and emission factors for modeling biogenic volatile organic compounds emissions from Hong Kong

This paper describes a study of local biogenic volatile organic compounds (BVOC) emissions from the Hong Kong Special Administrative Region (HKSAR). An improved land cover and emission factor database was developed to estimate Hong Kong emissions using MEGAN, a BVOC emission model developed by Guenther et al. (2006). Field surveys of plant species composition and laboratory measurements of emission factors were combined with other data to improve existing land cover and emission factor data. The BVOC emissions from Hong Kong were calculated for 12 consecutive years from 1995 to 2006. For the year 2006, the total annual BVOC emissions were determined to be 12,400 metric tons or 9.82 × 10⁹ g C (BVOC carbon). Isoprene emission accounts for 72%, monoterpene emissions account for 8%, and other VOCs emissions account for the remaining 20%. As expected, seasonal variation results in a higher emission in the summer and a lower emission in the winter, with emission predominantly in day time. A high emission of isoprene occurs for regions, such as Lowest Forest-NT North, dominated by broadleaf trees. The spatial variation of total BVOC is similar to the isoprene spatial variation due to its high contribution. The year to year variability in emissions due to weather was small over the twelve-year period (?1.4%, 2006 to 1995 trendline), but an increasing trend in the annual variation due to an increase in forest land cover can be observed (+7%, 2006 to 1995 trendline). The results of this study demonstrate the importance of accurate land cover inputs for biogenic emission models and indicate that land cover change should be considered for these models.     

Using a reduced Common Representative Intermediates (CRIv2-R5) mechanism to simulate tropospheric ozone in a 3-D Lagrangian chemistry transport model

A reduced chemical scheme (CRIv2-R5) which describes ozone formation from the tropospheric degradation of methane and 22 emitted non-methane hydrocarbons and oxygenated volatile organic compounds has been applied in a global-3D chemistry transport model (STOCHEM). The scheme, which contains 220 species in 609 reactions, has been used to simulate ozone and its precursors for the meteorological year of 1998 and the results have been compared with those from STOCHEM runs with its original chemistry. Compared with the original chemistry scheme, the degradation of a larger number of more reactive VOCs in the CRI scheme results in the formation (and their consequent transportation) of more NO_x active reservoirs thus leading to formation of more ozone away from land-based sources. Conversely, the more reactive VOCs also lead to greater removal of OH in continental areas and greater formation of OH in marine environments. STOCHEM run with the CRI scheme simulates more ozone (by up to 10 ppb), which results in better agreement with observed vertical ozone profiles. The CRI scheme transforms the globally and annually integrated ozone budget for the considered year in STOCHEM from a net loss of ?55 Tg yr^?1 to a net gain of +50 Tg yr^?1.     

Kinetic study of the OH reaction with some hydrochloroethers under simulated atmospheric conditions

Using the relative rate technique, rate constants for the gas-phase reactions of hydroxyl radicals with 2-chloroethyl methyl ether (k₁), 2-chloroethyl ethyl ether (k₂) and bis(2-chloroethyl) ether (k₃) have been measured. Experiments were carried out at (298 ± 2) K and atmospheric pressure using synthetic air as bath gas. Using n-pentane and n-heptane as reference compounds, the following rate constants were derived: k₁ = (5.2 ± 1.2) × 10^?12, k₂ = (8.3 ± 1.9) × 10^?12 and k₃ = (7.6 ± 1.9) × 10^?12, in units of cm³ molecule^?1 s^?1. This is the first experimental determination of k₂ and k₃ under atmospheric pressure. The rate constants obtained are compared with previous literature data and the observed trends in the relative rates of reaction of hydroxyl radicals with the ethers studied are discussed. The atmospheric implications of the results are considered in terms of lifetimes and fates of the hydrochloroethers studied.