Examination of a subgrid-scale parameterization for the transport of pollutants in a nonprecipitating cumulus cloud ensemble

The mass flux based subgrid-scale parameterization technique of Gidel (1983, J. geophys. Res.88, 6587–6599) is re-examined for use in Eulerian long-range transport models. Specifically, the parameterization scheme is incorporated into the STEM-II Eulerian transport/transformation/removal model and the model is used to investigate pollutant transport in a nonprecipitating cumulus cloud ensemble. The effects of entrainment, detrainment, evaporation and the transport by subsidence, updrafts and turbulent diffusion are included in the analysis. Presented simulation results indicate that the parameterization is able to treat the rapid vertical transport by cloud updrafts, enables the calculation of reaction rates based on subgrid-scale concentrations, and is readily adopted by Eulerian models.     

Fully coupled “online” chemistry within the WRF model

A fully coupled “online” Weather Research and Forecasting/Chemistry (WRF/Chem) model has been developed. The air quality component of the model is fully consistent with the meteorological component; both components use the same transport scheme (mass and scalar preserving), the same grid (horizontal and vertical components), and the same physics schemes for subgrid-scale transport. The components also use the same timestep, hence no temporal interpolation is needed. The chemistry package consists of dry deposition (“flux-resistance” method), biogenic emission as in [Simpson et al., 1995. Journal of Geophysical Research 100D, 22875–22890; Guenther et al., 1994. Atmospheric Environment 28, 1197–1210], the chemical mechanism from RADM2, a complex photolysis scheme (Madronich scheme coupled with hydrometeors), and a state of the art aerosol module (MADE/SORGAM aerosol parameterization).The WRF/Chem model is statistically evaluated and compared to MM5/Chem and to detailed photochemical data collected during the summer 2002 NEAQS field study. It is shown that the WRF/Chem model is statistically better skilled in forecasting O₃ than MM5/Chem, with no appreciable differences between models in terms of bias with the observations. Furthermore, the WRF/Chem model consistently exhibits better skill at forecasting the O₃ precursors CO and NO_y at all of the surface sites. However, the WRF/Chem model biases of these precursors and of other gas-phase species are persistently higher than for MM5/Chem, and are most often biased high compared to observations. Finally, we show that the impact of other basic model assumptions on these same statistics can be much larger than the differences caused by model differences. An example showing the sensitivity of various statistical measures with respect to the treatment of biogenic volatile organic compounds emissions illustrates this impact.     

Chemical transport models

Background, aim, and scope  Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. Materials and methods  The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. Results  This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. Discussion  In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO₂ by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. Conclusions  Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15–22%). Recommendations and perspectives  To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO₂) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.     

Measurements and analyses of nitrogen oxides and ozone in the yard and on the roof of a street-canyon in Suzhou

The concentrations of air pollutants such as nitrogen oxides and ozone characterised by very fast chemical reactions can significantly vary within urban street-canyon due to the short distances between sources and receptor. With the primary objective to analyse this issue, NO, NO₂, NO_x, O₃, BTX, and wind flow field were continuously measured for 1 week at two heights (a street-level yard and a 25-m-high rooftop) in an urban canyon in Suzhou (China). The yard ozone concentrations were found to be up to six times lower than on the roof. Different frequency distributions (FD), dynamical and chemical processes of the pollutant variations from yard to roof are discussed to explain the findings. The predominant factors for the dissimilar pollutant vertical diffusion at the two measurement locations were associated to dissimilar fluid-dynamic and heterogeneous removal effects that likely induced dissimilar ozone chemical processes relative to NO_x and BTX precursors.     

The contribution of in-cloud oxidation of SO2 to wet scavenging of sulfur in convective clouds

The sensitivity of in-cloud oxidation of SO₂ in corrective clouds to a number of chemical and physical parameters is examined. The parameterization of precipitation growth processes is based on the work of Scott (1978) and Hegg (1983). A chemical model predicts gas and aqueous phase distributions of soluble gases and in-cloud uncatalyzed oxidation of SO₂ by O₃ and H₂O₂. Sulfate aerosol and SO₂, CO₂, NH₃, H₂O₂ and O₃ gases and their aqueous phase dissociation products are treated.The results indicate that in-cloud conversion is an important removal mechanism for SO₂ and accounts for a significant fraction of the precipitation sulfate. However, except at low SO₂ concentrations, the precipitation sulfate concentration is insensitive to the initial SO₂ concentration; the sulfate concentration is most sensitive to the initial H₂O₂ and NH₃ concentrations. At low SO₂ concentrations, the precipitation sulfate concentration is determined primarily by the initial sulfate aerosol concentration. The feedback between sulfate production and pH is important in limiting SO₂ oxidation by O₃. If gas phase H₂O₂ of order 1 ppb is the major source of aqueous phase H₂O₂ for S(IV) oxidation, it is likely that the oxidation reaction is oxidant limited. The sulfate concentration is a decreasing function of the precipitation rate. At low rainfall rates (< 1 mm h⁻¹), ice phase growth decreases the sulfate concentration. However, the results are insensitive to an ice phase origin at moderate and high rainfall rates.     

A computer modelling study of the chemistry occurring during cloud formation over hills

A model has been developed which describes the chemical processes occurring during the formation of hill cloud. This model has been applied to the study of four cases which should be typical of the chemical conditions which can exist during the formation of such cloud over hills in the N of England.The results show that a wide range of chemical behaviour is possible in cloud droplets, depending on the history of the airmass. The composition of the water-soluble aerosol is of importance in determining the chemical composition of cloud droplets, with most of the sulphate and large proportions of the other ions being derived from this source. Gas-phase chemistry is also important, not only because it dictates the initial concentrations of oxidants such as H₂O₂ and O₃ and of acidic gases such as HNO₃ at the onset of condensation, but also because of the chemical processes which occur concurrently with those in clouds. These are important sources of radicals, as well as of nitric acid, which continue to dissolve in solution. Gasphase HCl, HNO₃ and NH₃ are all of primary importance in describing the acidity of cloud droplets. The various oxidation mechanisms for SO₂ can all be important, with their relative contributions varying from case to case. H₂O₂ does not invariably dominate this process, especially in winter when its concentration may be low.The chemistry of hill clouds is far from trivial and oversimplification can result in a misunderstanding of the behaviour which occurs. There is a need, rather, to consider a wide range of reactions if the nature of the overall processes and the impact of manmade emissions on the environment via this route are to be understood.     

An assessment of the impact of climate change on air quality at two UK sites

Possible effects of climate change on air quality are studied for two urban sites in the UK, London and Glasgow. Hourly meteorological data were obtained from climate simulations for two periods representing the current climate and a plausible late 21st century climate. Of the meteorological quantities relevant to air quality, significant changes were found in temperature, specific humidity, wind speed, wind direction, cloud cover, solar radiation, surface sensible heat flux and precipitation. Using these data, dispersion estimates were made for a variety of single sources and some significant changes in environmental impact were found in the future climate. In addition, estimates for future background concentrations of NO_x, NO₂, ozone and PM₁₀ upwind of London and Glasgow were made using the meteorological data in a statistical model. These showed falls in NO_x and increases in ozone for London, while a fall in NO₂ was the largest percentage change for Glasgow. Other changes were small. With these background estimates, annual-average concentrations of NO_x, NO₂, ozone and PM₁₀ were estimated within the two urban areas. For London, results averaged over a number of sites showed a fall in NO_x and a rise in ozone, but only small changes in NO₂ and PM₁₀. For Glasgow, the changes in all four chemical species were small. Large-scale background ozone values from a global chemical transport model are also presented. These show a decrease in background ozone due to climate change. To assess the net impact of both large scale and local processes will require models which treat all relevant scales.     

Evaluation of a reactive plume model with power plant plume data—application to the sensitivity analysis of sulfate and nitrate formation

A reactive plume model that treats secondary aerosol formation is used to investigate the major physical and chemical processes that affect the rate of sulfate and nitrate aerosol formation in power plant plumes. The reactive plume model is evaluated with experimental data collected in three power plant plumes, and model performance is shown to be quite satisfactory. One of these case studies is used to perform singleparameter and multi-parameter analyses of the sensitivity of sulfate and nitrate aerosol concentrations to various meteorological, air quality and chemical kinetic parameters. The results suggest that sulfate aerosol concentrations are most sensitive to relative humidity and temperature at high relative humidity, whereas nitrate aerosol concentrations are most sensitive to temperature, particularly at low relative humidity. The importance of the NO_x/reactive hydrocarbon chemistry to sulfate and nitrate aerosol formation is examined.     

Effects of chemical oxidation on sorption and desorption of PAHs in typical Chinese soils

In situ chemical oxidation is a commonly applied soil and groundwater remediation technology, but can have significant effects on soil properties, which in turn might affect fate and transport of organic contaminants. In this study, it was found that oxidation treatment resulted mainly in breakdown of soil organic matter (SOM) components. Sorption of naphthalene and phenanthrene to the original soils and the KMnO₄-treated soils was linear, indicating that hydrophobic partitioning to SOM was the predominant mechanism for sorption. Desorption from the original and treated soils was highly resistant, and was well modeled with a biphasic desorption model. Desorption of residual naphthalene after treating naphthalene-contaminated soils with different doses of KMnO₄ also followed the biphasic desorption model very well. It appears that neither changes of soil properties caused by chemical oxidation nor direct chemical oxidation of contaminated soils had a noticeable effect on the nature of PAH-SOM interactions.     

A study of the ozone formation by ensemble back trajectory-process analysis using the Eta–CMAQ forecast model over the northeastern U.S. during the 2004 ICARTT period

The integrated process rates (IPRs) estimated by the Eta–CMAQ model at grid cells along the trajectory of the air mass transport path were analyzed to quantitatively investigate the relative importance of physical and chemical processes for O₃ formation and evolution over the northeastern U.S. during the 2004 International Consortium for Atmospheric Research on Transport and Transformation (ICARTT) period. The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model is used to determine the back trajectory of air masses reaching the northeast by linking a downwind receptor to upwind source areas. The process analysis is applied to a high O₃ episode occurring on July 22, 2004 at three selected sites in the northeastern U.S. The process analysis at the location of the site shows that during the daytime, the O₃ concentrations in the surface layer are mainly enhanced by the vertical diffusion of O₃-rich air from aloft, followed by horizontal advection (HADV) and chemical production (CHEM), whereas dry deposition (DDEP) and vertical advection (ZADV) mainly deplete O₃ concentrations at the sites of Valley Central (VC), NY and Castle Spring (CS), NH. By integrating the effects of each process over the depth of the daytime planetary boundary layer (PBL), it was found that at the VC site, CHEM and HADV contributed about 53% and 41%, respectively, to O₃ levels within the PBL. This confirms the significance of regional transport of O₃ from the industrialized areas into the Northeast. On the other hand, the process analysis results for O₃ formation in moving air masses indicate that on July 22, large chemical production of O₃ along the transport path over the polluted urban regions leads to significant increase in O₃ in the air mass reaching the VC site, whereas the low chemical production of O₃ along the transport path over the low emission regions leads to the low O₃ concentration at the site of Belleayre Mountain (BM), NY. The dramatic buildup of O₃ concentration from 50 ppb to 102 ppb in the air masses before reaching the VC site after 12:00 EST on 7/22 indicates the significant impact of pollution from the northeastern urban corridor at this site. On the basis of the results at the CS site, it was found that high NO_x emissions along the transport path led to large chemical production of O₃ in the air mass reaching the CS site on July 22. In contrast, the low chemical production of O₃ associated with low emission (relatively clean conditions) along the transport path over the northern portions of the domain is responsible for the low O₃ concentration at the CS site on July 26.     

Comprehensive sensitivity analysis of an SO2 deposition model for three measurement sites: Consequences for SO2 deposition fluxes

A parameterization for the dry deposition of SO₂ is evaluated by means of four statistical sensitivity studies using a sampling-based approach to sensitivity analysis. Measurements over a heathland, grassland and coniferous forest are used to identify those meteorological and vegetation-specific factors which drive the deposition component of the model (model-based studies) and which are most important in the determination of model-measurement discrepancies. Sensitivity indices and scatter plots for monthly average fluxes indicate that the dominant factors are for most months related to non-stomatal components of the deposition flux.     

Evaluation of dynamic subgrid-scale models in large-eddy simulations of neutral turbulent flow over a two-dimensional sinusoidal hill

Large-eddy simulation (LES) is used to simulate neutral turbulent boundary-layer flow over a rough two-dimensional sinusoidal hill. Three different subgrid-scale (SGS) models are tested: (a) the standard Smagorinsky model with a wall-matching function, (b) the Lagrangian dynamic model, and (c) the recently developed scale-dependent Lagrangian dynamic model [Stoll, R., Porté-Agel, F., 2006. Dynamic subgrid-scale models for momentum and scalar fluxes in large-eddy simulation of neutrally stratified atmospheric boundary layers over heterogeneous terrain. Water Resources Research 42, W01409. doi:10.1029/2005WR003989]. The simulation results obtained with the different models are compared with turbulence statistics obtained from experiments conducted in the meteorological wind tunnel of the AES (Atmospheric Environment Service, Canada) [Gong, W., Taylor, P.A., Dörnbrack, A., 1996. Turbulent boundary-layer flow over fixed aerodynamically rough two-dimensional sinusoidal waves. Journal of Fluid Mechanics 312, 1–37]. We find that the scale-dependent dynamic model is able to account, without any tuning, for the local changes in the eddy-viscosity model coefficient. It can also capture the scale dependence of the coefficient associated with regions of the flow with strong mean shear and flow anisotropy. As a result, the scale-dependent dynamic model yields results that are more realistic than the ones obtained with the scale-invariant Lagrangian dynamic model.     

Transport and transformation of sulfur compounds over East Asia during the TRACE-P and ACE-Asia campaigns

On the basis of the recently estimated emission inventory for East Asia with a resolution of 1×1°, the transport and chemical transformation of sulfur compounds over East Asia during the period of 22 February through 4 May 2001 was investigated by using the Models-3 Community Multi-scale Air Quality (CMAQ) modeling system with meteorological fields calculated by the regional atmospheric modeling system (RAMS). For evaluating the model performance simulated concentrations of sulfur dioxide (SO₂) and aerosol sulfate (SO₄²⁻) were compared with the observations on the ground level at four remote sites in Japan and on board aircraft and vessel during the transport and chemical evolution over the Pacific and Asian Pacific regional aerosol characterization experiment field campaigns, and it was found that the model reproduces many of the important features in the observations, including horizontal and vertical gradients. The SO₂ and SO₄²⁻ concentrations show pronounced variations in time and space, with SO₂ and SO₄²⁻ behaving differently due to the interplay of chemical conversion, removal and transport processes. Analysis of model results shows that emission was the dominant term in regulating the SO₂ spatial distribution, while conversion of SO₂ to SO₄²⁻ in the gas phase and the aqueous phase and wet removal were the primary factors that controlled SO₄²⁻ amounts. The gas phase and the aqueous phase have the same importance in oxidizing SO₂, and about 42% sulfur compounds (25% in SO₂) emitted in the model domain was transported out, while about 57% (35% by wet removal processes) was deposited in the domain during the study period.     

Effects of dry deposition on the concentration-distributions of atmospheric pollutants within land- and sea-breeze circulations

In a land- and sea-breeze situation, effects of dry deposition on the dynamics of the concentrations of chemically reacting air pollutants are investigated using a transport/transformation/removal model with diurnally varying deposition velocities modeled in terms of the aerodynamic, surface, and residual resistances. The results show that the diurnally varying flows and eddy diffusivities, which are characteristic of the landand sea-breeze system, transfer the effects of dry deposition on the concentrations quickly to the upper layer over the land and sea surfaces. The dry deposition effect on one species can be transmitted to others through the network of chemical reactions, e.g. inclusion of dry deposition into the simulation resulted in the increase of hydrocarbon concentrations. It is also predicted that the dry deposition processes could remove a considerable part of emitted NO_x, and SO₂ from the local circulations, e.g. for 2 days about 40% of the emitted NO_x was removed by the dry deposition of NO, NO₂, HNO₃ and PAN and in the case of SO₂, 25 % by that of SO₂ and SO₄²⁻.     

Mercury cycling and species mass balances in four North American lakes

A mass balance model for mercury based on the fugacity concept is applied to Lake Superior, Lake Michigan, Onondaga Lake and Little Rock Lake to evaluate model performance, analyze cycling of three mercury species groups (elemental, divalent and methyl mercury), and identify important processes that determine the source-to-concentration relationship of the three mercury species groups in these lakes. This model application to four disparate ecosystems is an extension of previous applications of fugacity-based models describing mercury cycling. The model performs satisfactorily following site-specific parameterization, and provides an estimate of minimum rates of species interconversion that compare well with literature. Volatilization and sediment burial are the main processes removing mercury from the lakes, and uncertainty analyses indicate that air-water exchange of elemental mercury and water-sediment exchange of divalent mercury attached to particles are influential in governing mercury concentrations in water. Any new model application or field campaign to quantify mercury cycling in a lake should consider these processes as important.     

Numerical studies of acidification processes within and below clouds with a flow-through chemical reactor model

A flow-through chemical reactor model has been exercised to assess the importance of various oxidation reactions and cloud processes on wet removal and redistribution of atmospheric pollutants and to investigate the effect of in-cloud acidification on precipitation chemistry at the surface. Preliminary results indicate that in-cloud acidification accounts for more than 60% of the wet deposited acids derived from acidification of initial SO₂, that 42–57% of water-soluble, non-reactive NH₃ and HNO₃ are removed by wet deposition. The pseudo-first-order conversion rate of SO₂ to SO₄²⁻ ranges from 3 to 25% h ⁻¹ depending on initial and boundary conditions.Sensitivity studies have been carried out to test the importance of time evolution of clouds on partitioning of pollutants in the atmosphere and to investigate the variability of precipitation chemistry due to changes in rate constants. The distributions of NH₃ and HNO₃ are found to be dependent largely on the cloud microphysical parameters, while the distributions of H₂O₂ and SO₂ depend largely on initial conditions of both species. Individual physical and chemical mechanisms can determine the overall rate of sulfate wet deposition at different stages of cloud evolution.     

A statistical approach for estimating uncertainty in dispersion modeling: An example of application in southwestern USA

A method based on a statistical approach of estimating uncertainty in simulating the transport and dispersion of atmospheric pollutants is developed using observations and modeling results from a tracer experiment in the complex terrain of the southwestern USA. The method takes into account the compensating nature of the error components by representing all terms, except dispersion error and variance of stochastic processes. Dispersion error and the variance of the stochastic error are estimated using the maximum likelihood estimation technique applied to the equation for the fractional error. Mesoscale Model 5 (MM5) and a Lagrangian random particle dispersion model with three optional turbulence parameterizations were used as a test bed for method application. Modeled concentrations compared well with the measurements (correlation coefficients on the order of 0.8). The effects of changing two structural components (the turbulence parameterization and the model grid vertical resolution) on the magnitude of the dispersion error also were examined. The expected normalized dispersion error appears to be quite large (up to a factor of three) among model runs with various turbulence schemes. Tests with increased vertical resolution of the atmospheric model (MM5) improved most of the dispersion model statistical performance measures, but to a lesser extent compared to selection of a turbulence parameterization. Method results confirm that structural components of the dispersion model, namely turbulence parameterizations, have the most influence on the expected dispersion error.     

Monitoring of sources and atmospheric processes controlling air quality in an urban Mediterranean environment

Different monitoring parameters (PM mass concentrations, number–size distribution, black carbon, gaseous pollutants, and chemical composition, among others) are currently used in air quality studies. Urban aerosols are the result of several sources and atmospheric processes, which suggests that a single monitoring technique is insufficient to quantitatively evaluate all of them.This study assesses the suitability of a number of monitoring techniques (PM mass concentrations, number and size distribution of ultra-fine particles, levels of gaseous pollutants, and a complete chemical characterization of PM₁₀ and PM_2.5) by examining the response of those techniques to the different emission sources and/or atmospheric processes affecting an urban Mediterranean area (Barcelona, NE Spain).The results of this work reveal that the PM mass, the number concentration and the chemical composition give different, but complementary, information. Whereas the mineral matter, a key atmospheric aerosol component across the Mediterranean, is not properly quantitatively assessed by measuring sub-micrometric particles, the monitoring of the number concentration is indispensable to interpret the origin of specific aerosol episodes. Furthermore, the chemical composition yields very relevant information to deduce the causes of specific pollution episodes.The number concentration of ultra-fine particles in urban areas is strongly dependent upon vehicle exhaust emissions, which may cause adverse health impacts. Moreover, urban Mediterranean environments are favourable to produce nucleation-mode particles (<20 nm) with photochemical origin. In those cases, these particles are expected to be of high solubility and consequently their toxicity may differ from that of traffic-generated ultra-fine particles. Thus, the use of a single monitoring parameter to evaluate the health effects seems to be not enough.     

Analysis of a Farquhar-von Caemmerer-Berry leaf-level photosynthetic rate model for Populus tremuloides in the context of modeling and measurement limitations

The balance of mechanistic detail with mathematical simplicity contributes to the broad use of the Farquhar, von Caemmerer and Berry (FvCB) photosynthetic rate model. Here the FvCB model was coupled with a stomatal conductance model to form an [A,g_s] model, and parameterized for mature Populus tremuloides leaves under varying CO₂ and temperature levels. Data were selected to be within typical forest light, CO₂ and temperature ranges, reducing artifacts associated with data collected at extreme values. The error between model-predicted photosynthetic rate (A) and A data was measured in three ways and found to be up to three times greater for each of two independent data sets than for a base-line evaluation using parameterization data. The evaluation methods used here apply to comparisons of model validation results among data sets varying in number and distribution of data, as well as to performance comparisons of [A,g_s] models differing in internal-process components.