Simulating the fine and coarse inorganic particulate matter concentrations in a polluted megacity

A three dimensional chemical transport model (PMCAMx) is applied to the Mexico City Metropolitan Area (MCMA) in order to simulate the chemical composition and mass of the major PM₁ (fine) and PM_1–10 (coarse) inorganic components and determine the effect of mineral dust on their formation. The aerosol thermodynamic model ISORROPIA-II is used to explicitly simulate the effect of Ca, Mg, and K from dust on semi-volatile partitioning and water uptake. The hybrid approach is applied to simulate the inorganic components, assuming that the smallest particles are in thermodynamic equilibrium, while describing the mass transfer to and from the larger ones. The official MCMA 2004 emissions inventory with improved dust and NaCl emissions is used. The comparison between the model predictions and measurements during a week of April of 2003 at Centro Nacional de Investigacion y Capacitacion Ambiental (CENICA) “Supersite” shows that the model reproduces reasonably well the fine mode composition and its diurnal variation. Sulfate predicted levels are relatively uniform in the area (approximately 3 μg m^?3), while ammonium nitrate peaks in Mexico City (approximately 7 μg m^?3) and its concentration rapidly decreases due to dilution and evaporation away from the urban area. In areas of high dust concentrations, the associated alkalinity is predicted to increase the concentration of nitrate, chloride and ammonium in the coarse mode by up to 2 μg m^?3 (a factor of 10), 0.4 μg m^?3, and 0.6 μg m^?3 (75%), respectively. The predicted ammonium nitrate levels inside Mexico City for this period are sensitive to the physical state (solid versus liquid) of the particles during periods with RH less than 50%.     

A new condensed toluene mechanism for Carbon Bond: CB05-TU

Toluene is ubiquitous in urban atmospheres and is a precursor to tropospheric ozone and aerosol (smog). An important characteristic of toluene chemistry is the tendency of some degradation products (e.g., cresols and methyl-catechols) to form organic nitro and nitrate compounds that sequester NOx (NO and NO₂) from active participation in smog formation. Explaining the NOx sinks in toluene degradation has made mechanism development more difficult for toluene than for many other organic compounds. Another challenge for toluene is explaining sources of radicals early in the degradation process. This paper describes the development of a new condensed toluene mechanism consisting of 26 reactions, and evaluates the performance of CB05 with this new toluene scheme (Toluene Update, TU) against 38 chamber experiments at 7 different environmental chambers, and provides recommendations for future developments. CB05 with the current toluene mechanism (CB05-Base) under-predicted the maximum O₃ and O₃ production rate for many of these toluene–NOx chamber experiments, especially under low-NOx conditions ([NOx]_t=0 < 100 ppb). CB05 with the new toluene mechanism (CB05-TU) includes changes to the yields and reactions of cresols and ring-opening products, and showed better performance than CB05-Base in predicting the maximum O₃, O₃ formation rate, NOx removal rate and cresol concentration. Additional environmental chamber simulations with xylene–NOx experiments showed that the TU mechanism updates tended to improve mechanism performance for xylene.     

Prymnesium parvum: The Texas Experience1

Southard, Gregory M., Loraine T. Fries, and Aaron Barkoh, 2010. Prymnesium parvum: The Texas Experience. Journal of the American Water Resources Association (JAWRA) 46(1):14-23. DOI: 10.1111/j.1752-1688.2009.00387.x Abstract: Golden alga Prymnesium parvum was first identified in Texas during a fish kill investigation on the Pecos River in 1985. Since then golden alga kills occurred sporadically in a variety of waters in the western part of the state until 2001 when the alga became endemic in the Brazos, Canadian, Colorado, Red, and Rio Grande river systems, including the water supplies of two public fish hatcheries, the Possum Kingdom and Dundee state fish hatcheries. The increasing area adversely affected by the alga and frequent massive fish kills heightened public and political awareness and concerns regarding the ecological and economic impacts of P. parvum blooms. The Texas Parks and Wildlife Department (TPWD), the wildlife conservation agency of the state, responded to these concerns with a program to assess the ecological and economic impacts and to develop management options. To date 33 water bodies have been affected and losses are conservatively estimated at 34 million fish valued at US$13 million. Several sport fisheries, including smallmouth bass Micropterus dolomieu, striped bass Morone saxatilis, channel catfish Ictalurus punctatus, and blue catfish Ictalurus furcatus, have been severely affected. Additionally, 26 imperiled fish species occur in the affected water basins and some have been adversely affected. Economic losses associated with reduced fishing and other water-based recreational activities appear considerable. The combined economic losses to three counties (Palo Pinto, Stephens, and Young) surrounding Possum Kingdom reservoir for 2001 and 2003 were estimated at US$2.8 million and US$1.1 million, respectively. This paper describes how the TPWD responded to public and political concerns relative to the emergence of golden alga, its harmful effects to fisheries, and its historic and current statewide distribution.     

Reclaimed Mineland Curve Number Response to Temporal Distribution of Rainfall1

Warner, Richard C., Carmen T. Agouridis, Page T. Vingralek, and Alex W. Fogle, 2010. Reclaimed Mineland Curve Number Response to Temporal Distribution of Rainfall. Journal of the American Water Resources Association (JAWRA) 46(4): 724-732. DOI: 10.1111/j.1752-1688.2010.00444.x Abstract: The curve number (CN) method is a common technique to estimate runoff volume, and it is widely used in coal mining operations such as those in the Appalachian region of Kentucky. However, very little CN data are available for watersheds disturbed by surface mining and then reclaimed using traditional techniques. Furthermore, as the CN method does not readily account for variations in infiltration rates due to varying rainfall distributions, the selection of a single CN value to encompass all temporal rainfall distributions could lead engineers to substantially under- or over-size water detention structures used in mining operations or other land uses such as development. Using rainfall and runoff data from a surface coal mine located in the Cumberland Plateau of eastern Kentucky, CNs were computed for conventionally reclaimed lands. The effects of temporal rainfall distributions on CNs was also examined by classifying storms as intense, steady, multi-interval intense, or multi-interval steady. Results indicate that CNs for such reclaimed lands ranged from 62 to 94 with a mean value of 85. Temporal rainfall distributions were also shown to significantly affect CN values with intense storms having significantly higher CNs than multi-interval storms. These results indicate that a period of recovery is present between rainfall bursts of a multi-interval storm that allows depressional storage and infiltration rates to rebound.     

Uncertainty in health risks due to anthropogenic primary fine particulate matter from different source types in Finland

The emission-exposure and exposure-response (toxicity) relationships are different for different emission source categories of anthropogenic primary fine particulate matter (PM_2.5). These variations have a potentially crucial importance in the integrated assessment, when determining cost-effective abatement strategies. We studied the importance of these variations by conducting a sensitivity analysis for an integrated assessment model. The model was developed to estimate the adverse health effects to the Finnish population attributable to primary PM_2.5 emissions from the whole of Europe. The primary PM_2.5 emissions in the whole of Europe and in more detail in Finland were evaluated using the inventory of the European Monitoring and Evaluation Programme (EMEP) and the Finnish Regional Emission Scenario model (FRES), respectively. The emission-exposure relationships for different primary PM_2.5 emission source categories in Finland have been previously evaluated and these values incorporated as intake fractions into the integrated assessment model. The primary PM_2.5 exposure-response functions and toxicity differences for the pollution originating from different source categories were estimated in an expert elicitation study performed by six European experts on air pollution health effects. The primary PM_2.5 emissions from Finnish and other European sources were estimated for the population of Finland in 2000 to be responsible for 209 (mean, 95% confidence interval 6–739) and 357 (mean, 95% CI 8–1482) premature deaths, respectively. The inclusion of emission-exposure and toxicity variation into the model increased the predicted relative importance of traffic related primary PM_2.5 emissions and correspondingly, decreased the predicted relative importance of other emission source categories. We conclude that the variations of emission-exposure relationship and toxicity between various source categories had significant impacts for the assessment on premature deaths caused by primary PM_2.5.     

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.     

Influence of boundary conditions on CMAQ simulations over the Iberian Peninsula

The influence of chemical boundary conditions (BC) on the response of the Community Multiscale Air Quality (CMAQ) model over the Iberian Peninsula was investigated in this study. Three strategies to supply boundary conditions in the context of the Integrated Assessment Modelling System for the Iberian Peninsula (SIMCA) were tested. Alternative methods consist in providing BC from (1) fixed, time-independent, concentration profiles, (2) concentrations predicted in a CMAQ mother domain (48 km, 1 h resolution) and (3) concentration values from the GEOS-Chem chemical-transport global model (2 × 2.5°, 3 h resolution). High resolution (3 km) simulated concentrations of the main pollutants (NO₂, NO, SO₂, O₃, PM₁₀ and PM_2.5) were compared through a comprehensive statistical analysis including observational data from 165 monitoring stations all over the Iberian Peninsula. It was found that model sensitivity to BC for nitrogen and sulphur oxides was limited, being restricted to the vicinity of model boundaries. However, significant domain-wide differences were found when modelling ozone and PM depending on the BC provided to run the tests. Although model performance was affected by spatial and seasonal factors, the results indicate that model-derived, dynamic BC improved CMAQ predictions when compared to those based on static concentrations prescribed in the boundaries. Aggregated statistics suggest that the GEOS-Chem produced the best results for O₃ and PM_2.5 while NO₂ and PM₁₀ were slightly better predicted under the CMAQ nesting approach. Besides the statistical evaluation some other relevant issues in the context of Integrated Assessment Modelling (IAM) are discussed to gain a better insight into the suitability of the methods analyzed and limitations of downscaling methods. Despite being useful to get a better understanding of the role of BC in SIMCA, this study contributed to highlight model deficiencies and therefore to point out future research needs to improve IAM activities in the Iberian Peninsula.     

Sea salt concentrations across the European continent

The oceans are a major source for particles that play an important role in many atmospheric processes. In Europe sea salt may contribute significantly to particulate matter concentrations. We have compiled sodium concentration data as a tracer for sea salt for 89 sites in Europe to provide more insight in the distribution of sea salt across Europe. The annual average sea salt concentrations above land were estimated to range between 0.3 and almost 13 μg m^?3. Maximum concentrations are found at the Irish coast. At coastal sites along the Atlantic and North Sea coast concentrations tend to be around 5 μg m^?3. More inland locations up to about 300 km away from the coast tend to show concentrations between 2 and 5 μg m^?3, whereas sites further away from the coast are characterized by lower concentrations. An analysis of the representativity of the data with respect to a long term average showed that the long average is associated with a standard deviation of around 15%. The compilation of observations provides an improved overview of sea salt concentrations in Europe as well as an improved basis for model validation. Verification of the results of the LOTOS-EUROS model learned that the model represents well the spatial variability of the observed sea salt concentrations very well. However, the absolute concentrations are significantly overestimated due to large uncertainties in the emission and dry deposition parameterizations. Using the high explained variability in the gradients across Europe, the bias-corrected modelled distribution serves as a best estimate of the sea salt distribution across Europe for 2005.     

Evaluation of aerosol sources at European high altitude background sites with trajectory statistical methods

This study has investigated the influence of synoptic weather patterns and long-range transport episodes on the concentrations of several compounds related to different aerosol sources (EC, OC, SO₄^2?, Ca²⁺, Na⁺, K⁺, ²¹⁰Pb, levoglucosan and dicarboxylic acids) registered in PM10 or PM2.5 aerosol samples collected at three remote background sites in central Europe. Air mass back-trajectories arriving at these sites have been analysed by statistical methods. Firstly, air mass back-trajectories have been grouped into clusters. Each cluster corresponds to specific meteorological scenarios, which were extracted and discussed. Finally, redistributed concentration fields have been computed to identify the main potential source regions of the different key aerosol components. A marked seasonal pattern is observed in the occurrence of the different clusters, with fast westerly and northerly Atlantic flows during winter and weak circulation flows in summer. Spring and fall were characterised by advection of moderate flows from northeastern and eastern Europe. Significant inter-cluster differences were observed for concentrations of receptor aerosol components, with the highest concentrations of EC, OC, SO₄^2?, K⁺ and ²¹⁰Pb associated with local and mesoscale aerosol sources located over central Europe related to enhanced photochemical processes. Emissions produced by fossil fuel and biomass burning processes from the Baltic countries, Byelorussia, western regions of Russia and Kazakhstan in spring and fall also contribute to elevated levels of EC, OC, SO₄^2?, K⁺ and ²¹⁰Pb. In the summer period long-range transport episodes of mineral dust from North-African deserts were also frequently detected, which caused elevated concentrations of coarse Ca²⁺ at sites. The baseline aerosol concentrations in central Europe at the high altitude background sites were registered in winter, with the exception of coarse Na⁺. While the relatively high concentrations of Na⁺ can be explained by sea salt advected from the Atlantic, the low levels of other aerosol components are caused by efficient aerosol scavenging associated to advections of Atlantic air masses, as well as lower emissions of these species over the Atlantic compared to those over the European continent and very limited vertical air mass exchange over the continent.     

Polycyclic aromatic hydrocarbons in air on small spatial and temporal scales – II. Mass size distributions and gas-particle partitioning

Polycyclic aromatic hydrocarbons (PAHs) were measured together with inorganic air pollutants at two urban sites and one rural background site in the Banja Luka area, Bosnia and Hercegovina, during 72 h in July 2008 using a high time resolution (5 samples per day) with the aim to study gas-particle partitioning, aerosol mass size distributions and to explore the potential of a higher time resolution (4 h-sampling).In the particulate phase the mass median diameters of the PAHs were found almost exclusively in the accumulation mode (0.1–1.0 μm of size). These were larger for semivolatile PAHs than for non-volatile PAHs. Gas-particle partitioning of semivolatile PAHs was strongly influenced by temperature. The results suggest that the Junge–Pankow model is inadequate to explain the inter-species variation and another process must be significant for phase partitioning which is less temperature sensitive than adsorption. Care should be taken when interpreting slopes m of plots of the type log K_p = m log p_L⁰ + b based on 24 h means, as these are found sensitive to the time averaging, i.e. tend to be higher than when based on 12 h-mean samples.