Estimating overall persistence and long-range transport potential of persistent organic pollutants: a comparison of seven multimedia mass balance models and atmospheric transport models

Two different approaches to modeling the environmental fate of organic chemicals have been developed in recent years. The first approach is applied in multimedia box models, calculating average concentrations in homogeneous boxes which represent the different environmental media, based on intermedia partitioning, transport, and degradation processes. In the second approach, used in atmospheric transport models, the spatially and temporally variable atmospheric dynamics form the basis for calculating the environmental distribution of chemicals, from which also exchange processes to other environmental media are modeled. The main goal of the present study was to investigate if the multimedia mass balance models CliMoChem, SimpleBox, EVn-BETR, G-CIEMS, OECD Tool and the atmospheric transport models MSCE-POP and ADEPT predict the same rankings of the overall persistence (P(ov)) and long-range transport potential (LRTP) of POPs, and to explain differences and similarities between the rankings by the mass distributions and inter-compartment mass flows. The study was performed for a group of 14 reference chemicals. For P(ov), the models yield consistent results, owing to the large influence of phase partitioning parameters and degradation rate constants, which are used similarly by all models. Concerning LRTP, there are larger differences between the models than for P(ov), due to different LRTP calculation methods and spatial model resolutions. Between atmospheric transport models and multimedia fate models, no large differences in mass distributions and inter-compartment flows can be recognized. Deviations in mass flows are mainly caused by the geometrical design of the models.     

Modeling of oak pollen dispersal on the landscape level with a mesoscale atmospheric model

We present the extension and application of the mesoscale atmospheric meteorology model METRAS for dispersion of oak pollen. We incorporated functions for pollen emission, pollen viability and pollen deposition into METRAS and simulated pollen dispersal on a scale of up to 200 km. The basis of the simulations is a real landscape structure that includes topography, land use, and the location and size of oak stands. We simulated the oak pollen dispersion of one single oak stand with an estimated annual pollen production of 1 billion pollen grains/m² forest surface on two exemplary days of the flowering season in 2000. Depending on the meteorological situation of the simulated days, a pollen cloud with about 10 pollen/m³ may extend up to 30 km from the source. Downstream of the oak stand, approximately 1,000 pollen/m² deposited up to a distance of 25 km, and lower amounts of pollen deposited up to 100 km away. These values of pollen concentration and deposition lay within the range of published field studies. Overall, it is shown that mesoscale atmospheric models are applicable to simulate pollen dispersal on the landscape level.     

Evaluation of a Gaussian-modified dispersion model for atmospheric release from the Marcoule nuclear site

This paper presents a mathematical model of local pollutant dispersion designed to compute the concentration field above and around the Marcoule nuclear site. The model is based on integrating the classical turbulent diffusion equation, corrected (prior to integration) by experimental wind tunnel data obtained for a scaled-down model of the site. The computed results are compared with full-scale experimental observations at Marcoule in the case of neutral atmosphere. A comparison with the standard Gaussian model is also made. Finally, a critical analysis of the model is presented. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modelling the Distribution and Quality of Sand and Gravel Resources in 3D: a Case Study in the Thames Basin,UK

Three-dimensional (3D) models are often utilised to assess the presence of sand and gravel deposits. Expanding these models to provide a better indication of the suitability of the deposit as aggregate for use in construction would be advantageous. This, however, leads to statistical challenges. To be effective, models must be able to reflect the interdependencies between different criteria (e.g. depth to deposit, thickness of deposit, ratio of mineral to waste, proportion of ‘fines’) as well as the inherent uncertainty introduced because models are derived from a limited set of boreholes in a study region. Using legacy borehole data collected during a systematic survey of sand and gravel deposits in the UK, we have developed a 3D model for a 2400 km² region close to Reading, southern England. In developing the model, we have reassessed the borehole grading data to reflect modern extraction criteria and explored the most suitable statistical modelling technique. The additive log-ratio transform and the linear model of coregionalization have been applied, techniques that have been previously used to map soil texture classes in two dimensions, to assess the quality of sand and gravel deposits in the area. The application of these statistical techniques leads to a model which can be used to generate thousands of plausible realisations of the deposit which fully reflect the extent of model uncertainty. The approach offers potential to improve regional-scale mineral planning by providing an enhanced understanding of sand and gravel deposits and the extent to which they meet current extraction criteria.

     

Aluminium in UK rivers: a need for integrated research related to kinetic factors, colloidal transport, carbon and habitat

Dissolved aluminium concentrations ([Al]) in the <0.45 μm filtered fraction are described for 54 UK river sites covering rural, acidic/acid sensitive, agricultural and urban typologies, and wide pH range (4 to 11). High [Al] occurred under acidic conditions and for acid runoff neutralised by bicarbonate rich groundwater. Thermodynamic analysis indicates Al hydroxide/hydroxy-silicate oversaturation at circumneutral pH across the rivers, but undersaturation at lower/higher pH. The oversaturation reflects in part the presence of Al bearing colloids as indicated by (1) [Al] being correlated with components associated with both lithogenic (Fe, Ti and lanthanides) colloids and organic carbon, (2) baseflow studies using cross-flow ultrafiltration and (3) comparison of our data with Acid Waters Monitoring Network (AWMN) information on labile and non-labile Al. Tree harvesting and emission reductions of SO(x) in acidic and acid sensitive catchments in mid-Wales led to acidification reversal, lower [Al] and changing [H(+)] - [Al] relationships. The [Al] decline was confined to acidic conditions while [Al] increased during the later part of the monitoring period with a peak around 2002 for moorland and forested systems. Colloidal production across the flow range was indicated late in the record by comparison of our data with information collected by the AWMN for a site in mid-Wales. This production seems interlinked with organic carbon and with dissolved CO(2) changes. In order for further understanding of Al hydrogeochemistry in river systems there is a need to integrate research that moves from equilibrium to kinetic and colloidal consideration including the critical issues of organic and inorganic controls within the context of bioavailability and aquatic stress. The colloidal Al may well be of low environmental concern to fish and other factors such as habitat may well be critical.     

Computational Fluid Dynamics Modelling of the Pollution Dispersion and Comparison with Measurements in a Street Canyon in Helsinki

A measuring campaign was conducted in a street canyon (Runeberg St.) in Helsinki in 2003–2004. The concentrations of NO _x , NO₂, PM₁₀ and PM_2.5 were measured at street level and at roof level at an urban background location. This study utilises the data measured from 1 Jan to 30 April, 2004, when wind speed and direction measurements were also conducted on-site at the roof level. The computational fluid dynamics model ADREA-HF was used to compute the street concentrations, and the results were compared with the measurements. The predictions for the selected cases agreed fairly well (within < 25 % for 15 min average values) with the measured data, except for two cases: a windward flow in case of a low wind speed, and a moderate southerly flow parallel to the street canyon. The main reasons for the differences of predictions and measurements are the negligence of traffic-induced turbulence in the modelling and an under-prediction of ventilation of urban background air from a crossing street. Numerical results are presented for various example cases; these illustrate the formation of the vortices in the canyon in terms of the wind direction and speed and the influence of the characteristics of the flow fields on the concentration distributions.     

Testing the validity of a critical sulfur and nitrogen load model in southern Ontario, Canada, using soil chemistry data from MARYP

The validity of a steady-state massbalance model (Arp et al., 1996; referred to asARP) was tested using physicochemical soil data fromthe Monitoring Acid Rain Youth Program (MARYP). FourARP sites were matched with ten MARYP sites accordingto proximity, bedrock type and subsoil pH to test thevalidity of the ARP model for critical loadexceedances. Soil solution pH, base concentration andAl concentration from MARYP sites, which were wellmatched to ARP sites, validated the modelled criticalload exceedances. Higher exceedance areas wereassociated with more acidic pH and lower base andhigher Al concentrations from matched MARYP sites andvice versa. One ARP site was inappropriately matchedwith MARYP sites and could not be validated using baseand Al concentrations. This study also confirmed thesouthern limit of the zero critical load exceedanceisopleth from the model. However, variability of theother exceedance isopleths was noted due to thelimited number of sites used in the model. Thevalidation of these sites in the ARP model and thezero critical load exceedance isopleth nonethelessallows greater confidence in using this model as amanagement tool for acidic deposition.     

Mathematical models for accurate prediction of atmospheric visibility with particular reference to the seasonal and environmental patterns in Hong Kong

Atmospheric visibility impairment has gained increasing concern as it is associated with the existence of a number of aerosols as well as common air pollutants and produces unfavorable conditions for observation, dispersion, and transportation. This study analyzed the atmospheric visibility data measured in urban and suburban Hong Kong (two selected stations) with respect to time-matched mass concentrations of common air pollutants including nitrogen dioxide (NO(2)), nitrogen monoxide (NO), respirable suspended particulates (PM(10)), sulfur dioxide (SO(2)), carbon monoxide (CO), and meteorological parameters including air temperature, relative humidity, and wind speed. No significant difference in atmospheric visibility was reported between the two measurement locations (p > or = 0.6, t test); and good atmospheric visibility was observed more frequently in summer and autumn than in winter and spring (p < 0.01, t test). It was also found that atmospheric visibility increased with temperature but decreased with the concentrations of SO(2), CO, PM(10), NO, and NO(2). The results showed that atmospheric visibility was season dependent and would have significant correlations with temperature, the mass concentrations of PM(10) and NO(2), and the air pollution index API (correlation coefficients mid R: R mid R: > or = 0.7, p < or = 0.0001, t test). Mathematical expressions catering to the seasonal variations of atmospheric visibility were thus proposed. By comparison, the proposed visibility prediction models were more accurate than some existing regional models. In addition to improving visibility prediction accuracy, this study would be useful for understanding the context of low atmospheric visibility, exploring possible remedial measures, and evaluating the impact of air pollution and atmospheric visibility impairment in this region.     

Development of a personal monitoring method for nitrogen dioxide and sulfur dioxide with Sep-Pak C18 cartridge sampling and ion chromatographic determination

Personal monitoring methods for the determination of hourly integrated concentrations of NO2 and SO2 in ambient air have been developed. Triethanolamine (TEA)-impregnated C18 Sep-Pak cartridges were used to collect NO2 and SO2 simultaneously. After sampling, NO2 and SO2 as their nitrite, nitrate, sulfite and sulfate analogues were stripped from the cartridges with a solution of 5% methanol in distilled, deionized water (DDW) and then determined by ion chromatography. Laboratory tests were conducted to evaluate the sampling rate, collection and recovery efficiencies, breakthrough volumes, absorption capacity, interference and sample stability on the cartridge during storage. NO2 and SO2 detection limits of 0.3 and 0.4 ppb respectively for 1 h samples were obtained. Recoveries for both NO2 and SO2 exceeded 85%.     

Performances and application of a passive sampling method for the simultaneous determination of nitrogen dioxide and sulfur dioxide in ambient air

The performances and applicability of a diffusion tubesampler for the simultaneous measurements of NO₂ andSO₂ in ambient air were evaluated. SO₂ andNO₂are collected by the passive sampler using triethanolamine astrapping agent and are determined as sulphate and nitrite withion chromatography. The detection limit (2.3 g m^-3 ofNO₂ and 4.2 g m^-3 of SO₂ for two weekssampling) is adequate for the determination of concentrationsin urban and industrial areas. Precision of the method as RSDis in mean 5% for NO₂ and 12% for SO₂ at theconcentration levels in urban areas. Calibration of the methodwas performed in the field conditions by comparison between theresponses of sampler and the concentrations measured by thecontinuous monitors. High degree of linearity (correlationcoefficients > 0.8) is found between the passive sampler tubeand the continuous monitor data for both NO₂ and SO₂.To reduce the wind velocity influence on passive sampling ofdiffusion tubes, a protective shelter was tested in this study.The overall uncertainty of one measure for the optimised methodis estimated at 5 g m^-3 for NO₂ and 6g m^-3 for SO₂. Suitability of this passivesampling method for air pollution monitoring in urban areas wasdemonstrated by the results shown in this paper on a campaigncarried out in the French agglomeration.     

Calibration and verification of a dissolved oxygen management model for a highly polluted river with extreme flow variations in Pakistan

A dissolved oxygen (DO) model is calibrated and verified for a highly polluted River Ravi with large flow variations. The model calibration is done under medium flow conditions (431.5 m³/s), whereas the model verification is done using the data collected during low flow conditions (52.6 m³/s). Biokinetic rate coefficients for carbonaceous biochemical oxygen demand (CBOD) and nitrogenous biochemical oxygen demand (NBOD) (i.e, K _cr and K _n ) are determined through the measured CBOD and ammonia river profiles. The calculated values of K _cr and K _n are 0.36 day^?1 and 0.34 day^?1, respectively. The close agreement between the DO model results and the field values shows that the verified model can be used to develop DO management strategies for the River Ravi. The biokinetic coefficients are known to vary with degree of treatment (DOT) and therefore need to be adjusted for a rational water quality management model. The effect of this variation on level of treatment has been evaluated by using the verified model to attain a DO standard of 4 mg/L in the river using the biokinetic rate coefficients as determined during the model calibration and verification process. The required DOT in this case is found to be 96 %, whereas the DOT is 86 % if adjusted biokinetic rate coefficients are used to reflect the effect of wastewater treatment. The cost of wastewater treatment is known to increase exponentially as the removal efficiency increases; therefore, the use of appropriate biokinetic coefficients to manage the water quality in rivers is important.     

渤海湾无机氮、活性磷酸盐的变化对海洋初级生产力(叶绿素)的影响

通过三年海洋监测结果显示了渤海湾无机营养盐的变化对初级生产力的影响。其中无机氮严重超标达66%,无机氮与叶绿素的关系在枯水期、丰水期呈正相关;活性磷酸盐与叶绿素关系显著性不明显(仅枯水期显著),活性磷酸盐含量较低,限制了初级生产力。     

Load estimation and source apportionment of nonpoint source nitrogen and phosphorus based on integrated application of SLURP model, ECM, and RUSLE: a case study in the Jinjiang River, China

The nonpoint source (NPS) pollution is difficult to manage and control due to its complicated generation and formation. Load estimation and source apportionment are an important and necessary process for efficient NPS control. Here, an integrated application of semi-distributed land use-based runoff process (SLURP) model, export coefficients model (ECM), and revise universal soil loss equation (RUSLE) for the load estimation and source apportionment of nitrogen and phosphorus was proposed. The Jinjiang River (China) was chosen for the evaluation of the method proposed here. The chosen watershed was divided into 27 subbasins. After which, the SLURP model was used to calculate land use runoff and to estimate loads of dissolved nitrogen and phosphorus, and ECM was applied to estimate dissolved loads from livestock and rural domestic sewage. Next, the RUSLE was employed for load estimation of adsorbed nitrogen and phosphorus. The results showed that the 12,029.06 t?a^?1 pollution loads of total NPS nitrogen (TN) mainly originated from dissolved nitrogen (96.24 %). The major sources of TN were land use runoff, which accounted for 45.97 % of the total, followed by livestock (32.43 %) and rural domestic sewage (17.83 %). For total NPS phosphorous (TP), its pollution loads were 570.82 t?a^?1 and made up of dissolved and adsorbed phosphorous with 66.29 and 33.71 % respectively. Soil erosion, land use runoff, rural domestic sewage, and livestock were the main sources of phosphorus with contribution ratios of 33.71, 45.73, 14.32, and 6.24 % respectively. Therefore, land use runoff, livestock, and soil erosion were identified as the main pollution sources to influence loads of NPS nitrogen and phosphorus in the Jinjiang River and should be controlled first. The method developed here provided a helpful guideline for conducting NPS pollution management in similar watershed.     

Natural forest regrowth as a proxy variable for agricultural land abandonment in the Swiss mountains: a spatial statistical model based on geophysical and socio-economic variables

In many European mountain regions, natural forest regrowth on abandoned agricultural land and the related consequences for the environment are issues of increasing concern. We developed a spatial statistical model based on multiple geophysical and socio-economic variables to investigate the pattern of natural forest regrowth in the Swiss mountain area between the 1980s and 1990s. Results show that forest regrowth occurred primarily in areas with low temperature sum, intermediate steepness and soil stoniness as well as close to forest edges and relatively close to roads. Model results suggest that regions with weak labor markets are favored in terms of land abandonment and forest regrowth. We could not find an effect of population change on land abandonment and forest regrowth. Therefore, we conclude that decision makers should consider non-linearities in the pattern of forest regrowth and the fact that labor markets have an effect on land abandonment and forest regrowth when designing measures to prevent agricultural land abandonment and natural forest regrowth in the Swiss mountains.     

Zoobenthic variability associated with a flood and drought in the Hawkesbury estuary,New South Wales: Some consequences for environmental monitoring

Changes in the Hawkesbury zoobenthos associated with a major flood and drought are described and comparisons made in the patterns of change among estuarine reaches and among sites within reaches. The consequences for environmental monitoring and management are discussed. Replicate grab samples were taken from four sites in each of the lower, middle and upper reaches. Reductions in the mean number of species per grab (S) following the flood were significant only in the lower reaches and at one site in the middle reaches. Increases in S accompanied the drought in all reaches but intra-reach variation in temporal patterns occurred in both S and in the most abundant species found. Thus, major weather events are associated with temporal changes whose patterns differ on both small and large spatial scales. Consequently, the results from fixed-factor sampling designs, which are widely used, may be unrepresentative of other areas. Unfortunately, the alternative approach of stratified random sampling will probably be both prohibitively expensive and difficult to implement in the complex estuarine benthic habitat. Further, short-term studies will probably be grossly unrepresentative of natural temporal variation. Attempts to reduce expenses by using only one or two abundant species as characterising communities or as indicators of physicochemical conditions may be unreliable because of variation in both time and space in dominant species and the lack of pollution-response knowledge for local species.     

Preliminary results of laboratory toxicity tests with the mayfly, Isonychia bicolor (Ephemeroptera: Isonychiidae) for development as a standard test organism for evaluating streams in the Appalachian coalfields of Virginia and West Virginia

Selecting the most appropriate test species for sediment and water column assays has been a primary goal for ecotoxicologists. Standard test organisms and established test guidelines exist, but the USEPA-recommended species may not be the most sensitive organisms to anthropogenic inputs. This paper describes preliminary results of toxicity tests with the mayfly, Isonychia bicolor (Ephemeroptera). Results suggested that Isonychia were moderately sensitive to NaCl after 96 h with an average LC₅₀ value of 3.10 g NaCl per liter. This value decreased after 7 days of exposure, resulting in a mean LC₅₀ value of 1.73 g NaCl per liter. When exposed to a coal-mine-processed effluent, Isonychia generated LC₅₀ values that ranged from 13% to 39% effluent. I. bicolor were more sensitive to the coal processing effluent than Ceriodaphnia dubia with conductivity lowest observable effects concentration (LOEC) values for mayfly survivorship that ranged from 1,508 to 4,101 μS/cm, while LOEC values for C. dubia reproduction ranged from 2,132 to 4,240 μS/cm.     

Comparison of the use of mussels and semipermeable membrane devices for monitoring and assessment of accumulation of mutagenic pollutants in marine environment in combination with a novel microbiological mutagenicity assay

A novel microbiological mutagenicity assay, based on bioluminescence of a marine bacterium Vibrio harveyi mutant strain, potentially suitable for monitoring and assessment of mutagenic pollution of marine environment, has been described recently. Here, we tested the use of this assay, in combination with either mussels (Mytilus sp.) or semipermeable membrane devices (SPMDs), in assessment of accumulation of mutagens in marine water (samples of Baltic Sea water were tested). Either similar results were obtained in both systems or higher signals in the SPMD-based system were detected, depending on the tested water samples. We conclude that the use of both mussels and SPMDs in combination with the V. harveyi bioluminescence mutagenicity assay is a method suitable for monitoring and assessment of accumulation of mutagenic pollutants in marine environment, but in some cases the SPMD-based system may provide a more sensitive test.