Polycyclic aromatic hydrocarbons in Norwegian forest soils: impact of long range atmospheric transport

Levels of nine selected polycyclic aromatic hydrocarbons, PAHs, in surface soils from areas in southern and central Norway are presented. Levels in central Norway are generally low, while the southern Norway soils are about ten-fold higher with respect to 4 and 5 ring PAHs. Comparison with air quality data indicates long-range atmospheric transport to be the major source of the excess 4 and 5 ring PAHs in the south. Analyses of peat cores from ombrotrophic bogs support this assumption, and these provide a potentially useful approach for temporal studies of atmospheric PAH deposition. Analytical data for naphthalene in soils depend very much on the sampling and storage of the samples before analysis.     

Effects of broadleaf woodland cover on streamwater chemistry and risk assessments of streamwater acidification in acid-sensitive catchments in the UK

Streamwater was sampled at high flows from 14 catchments with different (0-78%) percentages of broadleaf woodland cover in acid-sensitive areas in the UK to investigate whether woodland cover affects streamwater acidification. Significant positive correlations were found between broadleaf woodland cover and streamwater NO3 and Al concentrations. Streamwater NO3 concentrations exceeded non-marine SO4 in three catchments with broadleaf woodland cover>or=50% indicating that NO3 was the principal excess acidifying ion in the catchments dominated by woodland. Comparison of calculated streamwater critical loads with acid deposition totals showed that 11 of the study catchments were not subject to acidification by acidic deposition. Critical loads were exceeded in three catchments, two of which were due to high NO3 concentrations in drainage from areas with large proportions of broadleaved woodland. The results suggest that the current risk assessment methodology should protect acid-sensitive catchments from potential acidification associated with broadleaf woodland expansion.     

Development and application of a three-dimensional aerosol chemical transport model,PMCAMx

A three-dimensional chemical transport model (PMCAMx) is used to simulate PM mass and composition in the eastern United States for a July 2001 pollution episode. The performance of the model in this region is evaluated, taking advantage of the highly time and size-resolved PM and gas-phase data collected during the Pittsburgh Air Quality Study (PAQS). PMCAMx uses the framework of CAMx and detailed aerosol modules to simulate inorganic aerosol growth, aqueous-phase chemistry, secondary organic aerosol formation, nucleation, and coagulation. The model predictions are compared to hourly measurements of PM_2.5 mass and composition at Pittsburgh, as well as to measurements from the AIRS and IMPROVE networks. The performance of the model for the major PM_2.5 components (sulfate, ammonium, and organic carbon) is encouraging (fractional errors are in general smaller than 50%). Additional improvements are possible if the rainfall measurements are used instead of the meteorological model predictions. The modest errors in ammonium predictions and the lack of bias for the total (gas and particulate) ammonium suggest that the improved ammonia inventory used is reasonable. The significant errors in aerosol nitrate predictions are mainly due to difficulties in simulating the nighttime formation of nitric acid. The concentrations of elemental carbon (EC) in the urban areas are significantly overpredicted. This is a problem related to both the emission inventory but also the different EC measurement methods that have been used in the two measurement networks (AIRS and IMPROVE) and the actual development of the inventory. While the ability of the model to reproduce OC levels is encouraging, additional work is necessary to confirm that that this is due to the right reasons and not offsetting errors in the primary emissions and the secondary formation. The model performance against the semi-continuous measurements in Pittsburgh appears to be quite similar to its performance against daily average measurements in a wide range of stations across the Eastern US. This suggests that the skill of the model to reproduce the diurnal variability of PM_2.5 and its major components is as good as its ability to reproduce the daily average values and also the significant value of high temporal resolution measurements for model evaluation.     

Environmental effects of transport: a model of optimal pricing and investment for the UK

This paper models the effects of imposing efficient pricing on interurban passenger travel in the United Kingdom and passenger travel in London. Both the internal and external costs of the main modes are estimated. The resulting efficient prices are used to predict the changes in demand for each mode, together with implied tax revenues, operators' aggregate surplus or deficit, and the implied investment needs. Efficient pricing and taxation of externalities does not lead to substantial shifts to modes of transport with lower external costs and thus does not substantially reduce the impact of transport on the environment. In the long run, demand changes are driven mainly by the growth of income. The assumption made about the appropriate scale economies present is an important determinant of the degree of change and overall effects. In the London model, sufficient revenue is generated to cover the deficits incurred by public transport operators.     

Groundwater flow, multicomponent transport and biogeochemistry: development and application of a coupled process model

A research tool for modeling the reactive flow and transport of groundwater contaminants in multiple dimensions is presented. Arbitrarily complex coupled kinetic–equilibrium heterogeneous reaction networks, automatic code generation, transfer-function based solutions, parameter estimation, high-resolution methods for advection, and robust solvers for the mixed kinetic–equilibrium chemistry are some of the features of reactive flow and transport (RAFT) that make it a versatile research tool in the modeling of a wide variety of laboratory and field experiments. The treatment of reactions is quite general so that RAFT can be used to model biological, adsorption/desorption, complexation, and mineral dissolution/precipitation reactions among others. The integrated framework involving automated code generation and parameter estimation allows for the development, characterization, and evaluation of mechanistic process models. The model is described and used to solve a problem in competitive adsorption that illustrates some of these features. The model is also used to study the development of an in situ Fe(II)-zone by encouraging the growth of an iron-reducing bacterium with lactate as the electron donor. Such redox barriers are effective in sequestering groundwater contaminants such as chromate and TCE.     

Occurrence of polybrominated diphenyl ethers in soil from the central Loess Plateau, China: role of regional range atmospheric transport

Very few studies were conducted in highland and depositional areas in studying the transport and behavior of polybrominated diphenyl ethers (PBDEs). In this study, surface soils were collected from Huan County to investigate the level, profile, and potential influence of PBDEs via regional range atmospheric transport in the central part of the Loess Plateau (CLP) of China, one of the most extensive areas of loess deposition in the world. PBDEs were ubiquitous and log-normally distributed in soils from the CLP with mean concentrations of 0.91 and 0.54 ng g⁻¹ for ΣPBDEs (sum of PBDE congeners except for BDE-209) and BDE-209, respectively. BDE-209 was predominated congener (43.5%), followed by BDE-47 (15.7%), 99 (10.7%), and 153 (7.5%). Further principal component analysis on congener profiles showed that PBDEs in the CLP originated from similar source(s). Additionally, significant differences in the ratios of BDE-47 to 99 and BDE-153 to 154 were found between soil samples and commercial products, indicating that they have undergone fractionation during the process of regional range atmospheric transport. The deposition of PBDEs in the CLP could be influenced by the sources from surrounding regions. For example, Xi’an may have potential influence to the CLP based on geographical analysis and concentrations comparison of PBDEs in gaseous. Therefore, more studies are needed to clarify the atmospheric transport and fate of PBDEs in this region.     

Development and application of a regional-scale atmospheric mercury model based on WRF/Chem: a Mediterranean area investigation

The emission, transport, deposition and eventual fate of mercury (Hg) in the Mediterranean area has been studied using a modified version of the Weather Research and Forecasting model coupled with Chemistry (WRF/Chem). This model version has been developed specifically with the aim to simulate the atmospheric processes determining atmospheric Hg emissions, concentrations and deposition online at high spatial resolution. For this purpose, the gas phase chemistry of Hg and a parametrised representation of atmospheric Hg aqueous chemistry have been added to the regional acid deposition model version 2 chemical mechanism in WRF/Chem. Anthropogenic mercury emissions from the Arctic Monitoring and Assessment Programme included in the emissions preprocessor, mercury evasion from the sea surface and Hg released from biomass burning have also been included. Dry and wet deposition processes for Hg have been implemented. The model has been tested for the whole of 2009 using measurements of total gaseous mercury from the European Monitoring and Evaluation Programme monitoring network. Speciated measurement data of atmospheric elemental Hg, gaseous oxidised Hg and Hg associated with particulate matter, from a Mediterranean oceanographic campaign (June 2009), has permitted the model’s ability to simulate the atmospheric redox chemistry of Hg to be assessed. The model results highlight the importance of both the boundary conditions employed and the accuracy of the mercury speciation in the emission database. The model has permitted the reevaluation of the deposition to, and the emission from, the Mediterranean Sea. In light of the well-known high concentrations of methylmercury in a number of Mediterranean fish species, this information is important in establishing the mass balance of Hg for the Mediterranean Sea. The model results support the idea that the Mediterranean Sea is a net source of Hg to the atmosphere and suggest that the net flux is ≈30 Mg year^?1 of elemental Hg.     

Micropollutants in urban stormwater: occurrence,concentrations, and atmospheric contributions for a wide range of contaminants in three French catchments

This study aimed at: (a) providing information on the occurrence and concentration ranges in urban stormwater for a wide array of pollutants (n?=?77); (b) assessing whether despite the differences between various catchments (land use, climatic conditions, etc.), the trends in terms of contamination level are similar; and (c) analyzing the contribution of total atmospheric fallout (TAF) with respect to sources endogenous to this contamination. The studied contaminants include conventional stormwater contaminants (polycyclic aromatic hydrocarbons (PAHs), Zn, Cu, Pb, etc.), in addition to poorly or undocumented pollutants such as nonylphenol and octylphenol ethoxylates (NPnEO and OPnEO), bisphenol A (BPA), polybrominated diphenyl ethers (PBDEs), a wide variety of pesticides, and various metals of relevance (As, Ti, Sr, V). Sampling and analysis were performed using homogeneous methods on three urban catchments with different land use patterns located in three distinct French towns. For many of these pollutants, the results do not allow highlighting a significant difference in stormwater quality at the scale of the three urban catchments considered. Significant differences were, however, observed for several metals (As, Cr, Cu, Ni, Sr and Zn), PAHs, and PBDEs, though this assessment would need to be confirmed by further experiments. The pollutant distributions between dissolved and particulate phases were found to be similar across the three experimental sites, thus suggesting no site dependence. Lastly, the contributions of TAF to stormwater contamination for micropollutants were quite low. This finding held true not only for PAHs, as previously demonstrated in the literature, but also for a broader range of molecules such as BPA, NPnEO, OPnEO, and PBDEs, whose high local production is correlated with the leaching of urban surfaces, buildings, and vehicles.     

The sensitivity of PM2.5 source-receptor relationships to atmospheric chemistry and transport in a three-dimensional air quality model

Air quality model simulations constitute an effective approach to developing source-receptor relationships (so-called transfer coefficients in the risk analysis framework) because a significant fraction of particulate matter (particularly PM2.5) is secondary (i.e., formed in the atmosphere) and, therefore, depends on the atmospheric chemistry of the airshed. In this study, we have used a comprehensive three-dimensional air quality model for PM2.5 (SAQM-AERO) to compare three approaches to generating episodic transfer coefficients for several source regions in the Los Angeles Basin. First, transfer coefficients were developed by conducting PM2.5 SAQM-AERO simulations with reduced emissions of one of four precursors (i.e., primary PM, sulfur dioxide (SO2), oxides of nitrogen (NOx), and volatile organic compounds) from each source region. Next, we calculated transfer coefficients using two other methods: (1) a simplified chemistry for PM2.5 formation, and (2) simplifying assumptions on transport using information limited to basin-wide emission reductions. Transfer coefficients obtained with the simplified chemistry were similar to those obtained with the comprehensive model for VOC emission changes but differed for NOx and SOz emission changes. The differences were due to the parameterization of the rates of secondary PM formation in the simplified chemistry. In 90% of the cases, transfer coefficients estimated using only basin-wide information were within a factor of two of those obtained with the explicit source-receptor simulations conducted with the comprehensive model. The best agreement was obtained for VOC emission changes; poor agreement was obtained for primary PM2.5.     

Development and application of a Lagrangian model to determine the origins of ozone episodes in the UK

This paper describes the further development and application of the Edinburgh–Lancaster Model for Ozone (ELMO). We replace straight-line back-trajectories with trajectories and associated meteorology supplied by the US National Oceanic and Atmospheric Administration Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) service to allow more realistic modelling of specific UK ozone episodes. We call this ELMO-2. Model performance is rigorously tested against observed ozone concentrations for two episodes recorded across 14 rural UK monitoring stations during the spring and summer of 1995. For both episodes, the afternoon concentrations (usually coinciding with the daily maxima) are captured well by the model and the diurnal ozone cycle is reproduced, although the amplitude in concentrations is generally smaller than the observed. The summer episode is investigated further through indicator species analysis and source attribution, and found to be mainly VOC-limited. European emissions account for the majority of ozone production. We demonstrate how improved modelling leads to better understanding of regional and local ozone production across the UK under episodic conditions.     

Sulfur processing in the marine atmospheric boundary layer: A review and critical assessment of modeling uncertainties

Sulfur is an extremely motile and vital element in the Earth's biogeochemical environment, one whose active redox chemistry maintains small reservoirs in the atmosphere and biosphere yet large fluxes through both. Essential for life, intimately linked to the climate state, and an important component of air quality, sulfur and its transport and processing in the atmosphere have been the subject of active research for several decades. This review article describes the current state of our understanding of the atmospheric sulfur cycle, focusing on the marine atmospheric boundary layer, with the aim of identifying the largest roots of uncertainty that most inhibit accurate simulation of sulfur cycling in the atmosphere. An overview of the emissions by phytoplankton and shipping, dispersion and entrainment in the marine boundary layer, and chemical processing by aerosols, clouds, and dry deposition is presented. Analysis of 20 contemporary modeling studies suggests that the greatest ambiguity in global sulfur cycling derives from (in descending order) wet deposition of aerosol sulfate, dry deposition of sulfur dioxide to the Earth's surface, and the heterogeneous oxidation of SO₂ in aerosols and clouds.     

A model of oscillatory transport in granular soils, with application to barometric pumping and earth tides

A simple algebraic model is proposed to estimate the transport of a volatile or soluble chemical caused by oscillatory flow of fluid in a porous medium. The model is applied to the barometric pumping of vapors in the vadose zone, and to the transport of dissolved species by earth tides in an aquifer. In the model, the fluid moves sinusoidally with time in the porosity of the soil. The chemical concentration in the mobile fluid is considered to equilibrate with the concentration in the surrounding matrix according to a characteristic time governed by diffusion, sorption, or other rate processes. The model provides a closed form solution, to which barometric pressure data are applied in an example of pore gas motion in the vadose zone. The model predicts that the additional diffusivity due barometric pumping in an unfractured vadose zone would be comparable to the diffusivity in stagnant pore gas if the equilibration time is 1 day or longer. Water motion due to the M2 lunar tide is examined as an example of oscillatory transport in an aquifer. It is shown that the tidal motion of the water in an aquifer might significantly increase the vertical diffusivity of dissolved species when compared to diffusion in an absolutely stagnant aquifer, but the hydrodynamic dispersivity due to tidal motion or gravitational flow would probably exceed the diffusivity due to oscillatory advection.     

Local and distant residence times of contaminants in multi-compartment models. Part II: application to assessing environmental mobility and long-range atmospheric transport

In Part I, the concepts of inherent, local and distant residence times (DRTs) were reviewed as metrics of the extent to which chemical discharges or emissions in one region or box are transported to distant regions. In this second part, the concepts are applied to geographically relevant systems to illustrate their applicability to the assessment of chemicals for long-range transport potential (LRTP). It is shown that the relative ranking of chemicals as characterized by the DRT method is similar to that of the characteristic travel distance concept. A DRT source-receptor matrix is developed that can express the chemical-specific potential of source regions to contaminate a specific receptor region of concern such as the Arctic. The matrix can be modified to identify for a specific source region the likely destinations of emissions as well as to assess the relative vulnerability of regions in the global environment to contaminants of concern.     

Dissolved and colloidal transport of Cd, Pb, and Zn in a silt loam soil affected by atmospheric industrial deposition.

As a result of processing of metal ores, trace metals have contaminated large areas of northern France. Metal migration from the soil to groundwater presents an environmental risk that depends on the physico-chemical properties of each contaminated soil. Soil water samples were obtained over the course of 1 year with zero-tension lysimeters from an acidic, loamy, metal contaminated soil. The average trace metal concentrations in the soil water were high (e.g. for Zn 11 mg l-1 under the surface horizon), but they varied during the sampling period. Zn concentrations were not correlated with pH or total organic carbon in the solutions but were correlated with Cd concentrations. On average, 95% of the Zn and Cd but only 50% of Pb was present in a dissolved form. Analytical transmission electron microscopy was used to identify the Zn or Pb carriers. Colloids containing Pb and Zn were biocolloids, whereas colloids containing only Zn were smectites.     

Illuminating reactive microbial transport in saturated porous media: demonstration of a visualization method and conceptual transport model

We demonstrate a method to study reactive microbial transport in saturated translucent porous media using the bacteria Pseudomonas fluorescens 5RL genetically engineered to carry a plasmid with bioluminescence genes inducible by salicylate. Induced bacteria were injected into a cryolite grain filled chamber saturated with a sterile non-growth-promoting (phosphorus limited) chemical mixture containing salicylate as an aromatic hydrocarbon analogue. The amount of light produced by the bacteria serves as an estimator of the relative efficiency of aerobic biodegradation since bioluminescence is dependent on both salicylate and oxygen but only consumes oxygen. Bioluminescence was captured with a digital camera and analyzed to study the evolving spatial pattern of the bulk oxygen consuming reactions. As fluid flow transported the bacteria through the chamber, bioluminescence was observed to initially increase until an oxygen depletion zone developed behind the advective front. Bacterial transport was modeled with the advection dispersion equation and oxygen concentration was modeled assuming bacterial consumption via Monod kinetics with consideration of additional effects of rate-limited mass transfer from residual gas bubbles. Consistent with previous measurements, bioluminescence was considered proportional to oxygen consumed. Using the observed bioluminescence, model parameters were fit that were consistent with literature values and produced results in good agreement with the experimental data. These findings demonstrate potential for using this method to investigate the complex spatial and temporal dynamics of reactive microbial transport in saturated porous media.     

Uncertainty analysis using a fugacity-based multimedia mass-balance model: Application of the updated EQC model to decamethylcyclopentasiloxane (D5)

The EQuilibrium Criterion (EQC) model developed and published in 1996 was recently revised to include improved treatment of input partitioning and reactivity data, temperature dependence and an easier sensitivity and uncertainty analysis. This New EQC model was used to evaluate the multimedia, fugacity-based fate of decamethylcyclopentasiloxane (D5; CAS No. 541-02-6) in the environment over a temperature range of 1–25 °C. In addition, Monte Carlo uncertainty analysis was used to quantitatively determine the influence of temperature and input partitioning and reactivity data on the behavior of D5 under various emission scenarios. Results indicated that emission mode was the most influential factor determining the fate and distribution of D5 in the model environment. When emitted to air and soil, D5 partitioned to and remained in the air compartment where rates of removal from degradation and advection processes were relatively rapid. In contrast, D5 emitted to water resulted in a substantial mass fraction of D5 being accumulated in the sediment compartment, where rates of removal from degradation and advection processes were slow. The mass distributions and fate of D5 in the model environment were strongly influenced by multiple input parameters, including temperature, the mode of emission (especially emission rate to water), K_OC and half-life in air. As temperature decreased from 25 °C to 1 °C, K_OC and half-life in air became increasingly more influential such that the mass distribution of D5 increased in air and decreased in sediment, resulting in decreased overall persistence.     

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.     

The arrival and discharge of conjugated estrogens from a range of different sewage treatment plants in the UK

The occurrence of free and conjugated estrogens was examined in a survey of eleven sewage treatment plants (STPs) and their discharge water in the United Kingdom using grab sampling. The STPs included trickling filter with and without tertiary treatment, and activated sludge with tertiary treatment. For three activated sludge plants both influent and effluent samples were compared. For a further 8 STPs only the effluent was examined. The estrone-3-sulphate, estradiol-3-sulphate and estriol-3-sulphate concentrations (up to 20 ng L⁻¹) were typically 5-fold that of the respective free estrogen concentration in the effluents. This represents a substantial additional ‘potential’ estrogen load arriving in the receiving waters. Estrone-3-glucuronide was found at 9 ng L⁻¹, estradiol-3-glucuronide at 7 ng L⁻¹, and estriol-3-glucuronide at 32 ng L⁻¹ in sewage influent. Except on one occasion, no glucuronide conjugates could be found in the effluent. The results suggest in most cases glucuronide conjugates will be completely transformed in sewage treatment whilst sulphate conjugates will only be partially removed.