Evaluation and comparison of multimedia mass balance models of chemical fate: application of EUSES and ChemCAN to 68 chemicals in Japan

The European Union System for Evaluation of Substances (EUSES) and the ChemCAN chemical fate model are applied to describe the fate of 68 chemicals on two spatial scales in Japan. Emission information on the chemicals has been obtained from Japan's Pollutant Release and Transfer Registry and available monitoring data gathered from government reports. Environmental concentrations calculated by the two models for the four primary environmental media of air, water, soil and sediment agree within a factor of 3 for over 70% of the data, and within a factor of 10 for over 87% of the data. Reasons for certain large discrepancies are discussed. Concentrations calculated by the models are generally consistent with the lower range of concentrations that are observed in the environment. Agreement between modeled and observed concentrations is considerably improved by including an estimate of the advective input of chemicals in air from outside Japan. The agreement between the EUSES and ChemCAN models suggests that results of individual chemical assessments are not likely to be significantly affected by the choice of chemical fate model. Primary sources of discrepancy between modeled and observed concentrations are believed to be uncertainties in emission rates, degradation half-lives, and the lack of data on advective inflow of contaminants in air.     

Operator-splitting procedures for reactive transport and comparison of mass balance errors

Operator-splitting (OS) techniques are very attractive for numerical modelling of reactive transport, but they induce some errors. Considering reactive mass transport with reversible and irreversible reactions governed by a first-order rate law, we develop analytical solutions of the mass balance for the following operator-splitting schemes: standard sequential non-iterative (SNI), Strang-splitting SNI, standard sequential iterative (SI), extrapolating SI, and symmetric SI approaches. From these analytical solutions, the operator-splitting methods are compared with respect to mass balance errors and convergence rates independently of the techniques used for solving each operator. Dimensionless times, NOS, are defined. They control mass balance errors and convergence rates. The following order in terms of decreasing efficiency is proposed: symmetric SI, Strang-splitting SNI, standard SNI, extrapolating SI and standard SI schemes. The symmetric SI scheme does not induce any operator-splitting errors, the Strang-splitting SNI appears to be O(N2OS) accurate, and the other schemes are first-order accurate.     

Very high correlations between fresh weight and lipid-adjusted PCB-153 serum concentrations: irrespective of fasting status, age, body mass index, gender, or exposure distributions

Background and aim

Persistent organochlorine pollutants (POPs) have in epidemiological studies been associated with health hazards. The exposure window might be crucial and it is, accordingly, important to analyze exposure concentrations in samples collected at the window of interest. This makes the use of stored blood samples extremely interesting. However, small blood volumes from stored samples in biobanks means that the researchers in epidemiological studies most often have to select a limited number of biomarkers representing the POP exposure and abstain from analyses of lipid concentrations. This raises the question of the necessity of lipid adjustment. The aim of the present study was to evaluate the correlation between fresh weight and lipid adjusted concentrations.

Methods

We included data from studies were we have used the PCB congener 2,2′,4,4′,5,5′-hexachlorobiphenyl (PCB-153) as a biomarker for POP exposure and where fresh weight as well as lipid adjusted PCB-153 concentrations are known (n = 1628). The variety of the cohort members in the different studies gave us the possibility to better understand how well these measures correlate.

Results

When all individuals were included in the analyses a correlation coefficient of 0.95 between fresh weight and lipid adjusted serum concentrations was obtained. Restricting the analyses to specific subgroups (based on gender, fasting status, age and BMI) gave correlation coefficients between 0.85 and 0.97.

Conclusions

This study does strongly support the use of samples from biobanks in epidemiological studies for analyzing PCBs, at least regarding PCB-153, even though it is not possible to lipid adjust the results.     

Strategies for including vegetation compartments in multimedia models

The incentives for including vegetation compartments in multimedia Level I, II and III fugacity calculations are discussed and equations and parameters for undertaking the calculations suggested. Model outputs with and without vegetation compartments are compared for 12 non-ionic organic chemicals with a wide variety of physical-chemical properties. Inclusion of vegetation compartments is shown to have a significant effect on two classes of chemicals: (1) those that are taken up by atmospheric deposition and (2) those that are taken up by transpiration through the plant roots. It is suggested that uptake from the atmosphere is important for chemicals with logK(OA) greater than 6 and a logK(AW) of greater than -6. Plant uptake by transpiration is important for chemicals with logK(OW) less than 2.5 and a logK(AW) of less than -1. At logK(OA) > 9 atmospheric uptake is dominated by particle-bound deposition and the importance of partitioning to vegetation is largely dependent on the relative magnitude of the particle deposition velocities to soil and vegetation. These property ranges can be used to determine if a chemical will significantly partition to vegetation. If the chemical falls outside the property ranges of the two classes it will probably be unnecessary to include vegetation in models for assessing environmental fate. The amount of chemical predicted to partition to vegetation compartments in the model is shown to be highly sensitive to certain model assumptions. Further experimental research is recommended to obtain more reliable equations describing equilibrium partitioning and uptake/depuration kinetics.     

Systematic comparison of ILWAS, MAGIC, and ETD watershed acidification models: 3. Mass balance budgets for acid neutralizing capacity

Three watershed acidification models-ILWAS, MAGIC, and ETD-were quantitatively compared to determine model structural differences by using a combination of input mapping and ANC mass balance budgets. Input mapping is a set of rules and algorithms to ensure that consistent input values were simultaneously derived for all three models. ANC budget analysis under current SO4(2-) deposition and a 70% reduction in SO4(2-) deposition allows examination of the relative importance of biogeochemical processes in affecting predictions of ANC or predicted changes in ANC. Model inputs were based on two dissimilar watersheds having characteristics typical of watersheds in the northeastern US. After mapping inputs, the three models predicted values of outflow ANC fluxes that were similar among the models for each watershed and deposition scenario. Within each watershed, the changes in outflow ANC fluxes between the scenarios were similar for the three models. Terrestrial weathering was the major source of ANC for all three models for both watersheds and deposition scenarios. The contributions of other processes to the ANC of the two watersheds were, under certain conditions, model-specific. Cation exchange was responsible for changes in ANC when deposition decreased for the three models. Other processes responsible for changes in ANC between scenarios were SO4(2-) sorption (for MAGIC) and in-lake weathering (for ETD). The processes responsible for the change in ANC from a change in deposition (cation exchange, SO4(2-) sorption, and in-lake weathering) were different from the processes contributing to the absolute ANC for a given deposition scenario (terrestrial weathering). The budget analysis complements an earlier Monte Carlo analysis that showed that the three models are structurally different and that predictions viewed on a relative scale are more similar than absolute scale predictions.     

Development of continental scale multimedia contaminant fate models: integrating GIS

The incentives and approaches for modelling chemical fate at a continental scale are discussed and reviewed. It is suggested that a multi-media model consisting of some 20-30 regions, each of which contains typically seven environmental compartments represents a reasonable compromise between the issues of the need for detailed resolution, avoidance of excessive data demands and inherent complexity and transparency. Strategies adopted in compiling the Berkley-Trent (BETR) model for North America are discussed and used to illustrate the issues of selecting appropriate number and nature of segments, treatment of air and water flows and the acquisition of environmental data. It is suggested that GIS software can play a valuable role in gathering and processing such data and in the display and interpretation of the results of the model assessment. The BETR model will be a useful tool for describing the nature of persistence and long-range transport of chemicals of concern in the North American environment.     

Mercury in the Mediterranean. Part 2: processes and mass balance

Mass balance of contaminants can provide useful information on the processes that influence their concentrations in various environmental compartments. The most important sources, sinks and the equilibrium or non-equilibrium state of the contaminant in individual environmental compartments can also be identified. Using the latest mercury speciation data, the results of numerical models and the results of recent studies on mercury transport and transformation processes in the marine environment, we have re-evaluated the total mercury (HgT) mass balance in the Mediterranean Sea. New calculations have been performed employing three distinct marine layers: the surface layer, the thermocline and the deep sea. New transport mechanisms, deep water formation and density-driven sinking and upwelling, were included in the mass balance calculations. The most recent data have even enabled the calculation of an approximate methylmercury (MeHg) mass balance. HgT is well balanced in the entire Mediterranean, and the discrepancies between inputs and outputs in individual layers do not exceed 20 %. The MeHg balance shows larger discrepancies between gains and losses due to measurement uncertainties and gaps in our knowledge of Hg species transformation processes. Nonetheless, the main sources and sinks of HgT (deposition and evasion) and MeHg (fluxes from sediment, outflow through the Gibraltar Strait) are in accordance with previous studies on mercury in the Mediterranean Basin. Mercury in the Mediterranean fish harvest is the second largest MeHg sink; about 300 kg of this toxic substance is consumed annually with sea food.     

Estimating half-lives of pesticides in/on vegetation for use in multimedia fate and exposure models

Degradation half-lives in/on vegetation are needed in environmental risk assessment of pesticides, but these data are often not available for most active ingredients. To address this, we first correlated experimental soil degradation half-life data of 41 pesticides obtained from the reviewed literature with the corresponding experimental half-lives on plant surface. Degradation half-lives in soil were found to be four times slower compared with half-lives on plant surfaces. In a second step, we explored measured plant surface half-lives directly with those in vegetation. The results were validated by comparing computed values with results obtained from an experimental set-up. The uptake and dissipation of alpha-cypermethrin (insecticide) and bromopropylate (acaricide) was studied by detecting pesticide residues in whole and peeled tomato fruits using gas chromatography. Half-lives within vegetation were found to be four times faster compared with plant surface half-lives. Using this experimental based approach, it is concluded that the estimation of degradation half-lives of pesticides in/on vegetation to be used as input data in environmental mass balance models can be directly correlated from the more abundant ready experimental degradation half-life data for soil.     

应用化学质量平衡模型解析烟台市污染源的排放贡献率

根据烟台地区污染排放特点,应用化学质量平衡模型（CMB8）进行污染源解析分析,得出每个污染源对于受体地区颗粒物体样本的排放贡献率。最终的模型结果显示,建筑尘和居民烯煤锅炉排放,海洋尘和冶炼厂排放对于大气环境污染也有一定的贡献,以上分析可以为城市污染源治理提供一定的依据。     

Critical loads for Finnish lakes: comparison of three steady-state models

Three simple steady-state water-chemistry models are used to calculate critical loads of sulfur for lakes in Finland. Because of the high concentrations of organic matter in Finnish lakes, the influence of organic anions on the calculation of critical loads has been given special attention. The first two methods are well known ion-balance methods which have been used in many previous lake-acidification studies. The third method, developed for this study, includes the numerical solution of equilibrium equations for organic anions, inorganic carbon species and inorganic monomeric aluminum. The original pH and aluminum concentration of the lakes are estimated with this model, and a method to estimate the original acid neutralizing capacity (ANC) by simulating a Gran-titration is also tested on the lake data. Uncertainty in the predictions is estimated by varying the most critical model parameters.     

Calculation of environmental concentration and comparison of output for existing chemicals using regional multimedia modeling

The environmental fate of 40 existing chemicals is discussed using the EUSES multimedia distribution and risk assessment model with site-specific parameter setting in an urban area of Japan including a highly industrial region. There has been a strong need to assess the environmental fate of a huge number of existing chemicals. Data on the emission amounts of chemicals are essential for such prediction, and PRTR surveys may yield this data. The study delivered the following results: (1) Volatile compounds with large amounts of emission showed higher predicted concentrations in air, and the concentrations of several compounds agreed well with averaged monitoring data within an order of magnitude. (2) A close relationship was found between the concentration of water and that of sediment, suggesting that the fate of chemicals in sediment essentially depended on the water environment. (3) A group of volatile solvents had high mass distribution ratio to air. Some compounds having high solubility in water were also included in that group due to the high ratio of air emission. Highly hydrophobic compounds with logK(OW) larger than 6.0 showed a high distribution ratio to soil and sediment. (4) Volatile compounds were mostly taken through air. The exposure through fish is a dominant pathway for highly hydrophobic compounds. (5) Exposure ratio could be gathered from physicochemical properties. The exposure from fish intake was roughly estimated by logK(OW), whereas exposure from air and water intake was difficult to estimate simply by vapor pressure and solubility in water, respectively.     

Occurrence of pharmaceuticals in a municipal wastewater treatment plant: mass balance and removal processes

Occurrence and removal efficiencies of fifteen pharmaceuticals were investigated in a conventional municipal wastewater treatment plant in Michigan. Concentrations of these pharmaceuticals were determined in both wastewater and sludge phases by a high-performance liquid chromatograph coupled to a tandem mass spectrometer. Detailed mass balance analysis was conducted during the whole treatment process to evaluate the contributing processes for pharmaceutical removal. Among the pharmaceuticals studied, demeclocycline, sulfamerazine, erythromycin and tylosin were not detected in the wastewater treatment plant influent. Other target pharmaceuticals detected in wastewater were also found in the corresponding sludge phase. The removal efficiencies of chlortetracycline, tetracycline, sulfamerazine, acetaminophen and caffeine were >99%, while doxycycline, oxytetracycline, sulfadiazine and lincomycin exhibited relatively lower removal efficiencies (e.g., <50%). For sulfamethoxazole, the removal efficiency was approximately 90%. Carbamazepine manifested a net increase of mass, i.e. 41% more than the input from the influent. Based on the mass balance analysis, biotransformation is believed to be the predominant process responsible for the removal of pharmaceuticals (22% to 99%), whereas contribution of sorption to sludge was relatively insignificant (7%) for the investigated pharmaceuticals.     

Prediction of overall persistence and long-range transport potential with multimedia fate models: robustness and sensitivity of results

The hazard indicators persistence (P) and long-range transport potential (LRTP) are used in chemicals assessment to characterize chemicals with regard to the temporal and spatial extent of their environmental exposure. They are often calculated based on the results of multimedia fate models. The environmental and substance-specific input parameters of such models are subject to a range of methodological uncertainties and also influenced by natural variability. We employed probabilistic uncertainty analysis to quantify variance in P and LRTP predictions for chemicals with different partitioning and transport behavior. Variance found in the results is so large that it prevents a clear distinction between chemicals. Additionally, only small improvements are observed when evaluating the results relative to a benchmark chemical. This can be explained by the dominance of substance-specific parameters and the only small direct influence of environmental parameters on P and LRTP as model outcomes. The findings underline the importance of learning how environmental conditions cause variability in substance behavior for improved substance ranking and classification.     

Electrolytic transformation of ordinance related compounds (ORCs) in groundwater: laboratory mass balance studies

Electrolytic reactive barriers (e(-) barriers) consist of closely spaced permeable electrodes installed across a groundwater contaminant plume in a permeable reactive barrier format. Application of sufficient potential to the electrodes results in sequential oxidation and reduction of the target contaminant. The objective of this study was to quantify the mass distribution of compounds produced during sequential electrolytic oxidation and reduction of ordinance related compounds (ORCs) in a laboratory analog to an e(-) barrier. In this study, a series of column tests were conducted using RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) and TNT (2,4,6-trinitrotoluene) as representative ORCs. The experimental setup consisted of a plexiglass column packed with quartz-feldspar sand to simulate aquifer conditions. A single set of porous electrodes consisting of expanded titanium-mixed metal oxide mesh was placed at the midpoint of the sand column as a one-dimensional analog to an e(-) barrier. Constant current of 20mA (variable voltage) was applied to the electrode set. Initial studies involved quantification of reaction products using unlabeled RDX and TNT. Approximately 70% of the influent concentration was transformed, in one pass, through sequential oxidation-reduction for both contaminants. Following the unlabeled studies, (14)C labeled RDX and TNT were introduced to determine the mass balance. An activity balance of up to 96% was achieved for both (14)C-RDX and (14)C-TNT. For both contaminants, approximately 21% of the influent activity was mineralized to (14)CO(2). The proportion of the initial activity in the dissolved fraction was different for the two test contaminants. Approximately 30% of the initial (14)C-RDX was recovered as unreacted in the dissolved phase. The balance of the (14)C-RDX was recovered as non-volatile, non-nitroso transformation products. None of the (14)C-RDX was sorbed to the column sand packing. For (14)C-TNT approximately 51% of the initial activity was recovered in the dissolved phase, the majority was unreacted TNT. The balance of the (14)C-TNT was either sorbed to the sand packing (approximately 24%) or dissolved/mineralized as unidentified ring cleavage products ( approximately 4%).     

Mercury mass balance study in Wujiangdu and Dongfeng Reservoirs, Guizhou, China

From October 2003 to September 2004, we conducted a detailed study on the mass balance of total mercury (THg) and methylmercury (MeHg) of Dongfeng (DF) and Wujiangdu (WJD) reservoirs, which were constructed in 1992 and 1979, respectively. Both reservoirs were net sinks for THg on an annual scale, absorbing 3319.5 g km⁻² for DF Reservoir, and 489.2 g km⁻² for WJD Reservoirs, respectively. However, both reservoirs were net sources of MeHg to the downstream ecosystems. DF Reservoir provided a source of 32.9 g MeHg km⁻² yr⁻¹, yielding 10.3% of the amount of MeHg that entered the reservoir, and WJD Reservoir provided 140.9 g MeHg km⁻² yr⁻¹, yielding 82.5% of MeHg inputs. Our results implied that water residence time is an important variable affecting Hg methylation rate in the reservoirs. Our study shows that building a series of reservoirs in line along a river changes the riverine system into a natural Hg methylation factory which markedly increases the %MeHg in the downstream reservoirs; in effect magnifying the MeHg buildup problem in reservoirs.     

A comparison of nonlinear regression and neural network models for ground-level ozone forecasting

A hybrid nonlinear regression (NLR) model and a neural network (NN) model, each designed to forecast next-day maximum 1-hr average ground-level O3 concentrations in Louisville, KY, were compared for two O3 seasons--1998 and 1999. The model predictions were compared for the forecast mode, using forecasted meteorological data as input, and for the hindcast mode, using observed meteorological data as input. The two models performed nearly the same in the forecast mode. For the two seasons combined, the mean absolute forecast error was 12.5 ppb for the NLR model and 12.3 ppb for the NN model. The detection rate of 120 ppb threshold exceedances was 42% for each model in the forecast mode. In the hindcast mode, the NLR model performed marginally better than the NN model. The mean absolute hindcast error was 11.1 ppb for the NLR model and 12.9 ppb for the NN model. The hindcast detection rate was 92% for the NLR model and 75% for the NN model.     

Completely autotrophic nitrogen-removal over nitrite in lab-scale constructed wetlands: evidence from a mass balance study

A mass-balance study was carried out to investigate the transformation of nitrogenous pollutants in vertical flow wetlands. Landfill leachate containing low BOD, but a high concentration of ammonia, was treated in four wetland columns under predominately aerobic conditions. Influent total nitrogen in the leachate consisted mainly of ammonia with less than 1% nitrate and nitrite, and negligible organic nitrogen. There was a substantial loss of total nitrogen (52%) in one column, whereas other columns exhibited zero to minor losses (<12%). Net nitrogen loss under study conditions was unexpected. Correlations between pH, nitrite and nitrate concentrations indicated the removal of nitrogen under study conditions did not follow the conventional, simplistic, chemistry of autotrophic nitrification. Through mass-balance analysis, it was found that CANON (Completely Autotrophic Nitrogen-removal Over Nitrite) was responsible for the transformation of nitrogen into gaseous form, thereby causing the loss of nitrogen mass. The results show that CANON can be native to aerobic engineered wetland systems treating wastewater that contains high ammonia and low BOD.     

Spatio-temporal behaviour and mass balance of fluorine in forest soils near an aluminium smelting plant: short- and long-term aspects

With the help of a large number of monitoring sites, the behaviour of F in forested ecosystems of a formerly polluted area (M?hlin and Rheinfelden, Switzerland) could be studied over 30 years. An aluminium smelting plant originated the pollution of this area in the past: however, after the installation of a filtration plant in 1958 the F emissions were reduced and since 1991 almost absent the primary production of aluminium was stopped. The present-day area with elevated F contents (water-soluble F with >20 mg/kg) is restricted to a radius of about 1000 m from the plant. In 1969 this area had a radius of about 3 and 5 km. Between 1969 and 2000 a significant decrease in the soluble F content in the soil was observed together with a substantial decrease of F in the vegetation. The net losses of water-soluble F in the soil were in the range of 35 up to more than 70% of the original concentration and the F losses in the vegetation between 60 up to more than 80%. After the reduction of high F deposition rates the accumulated SOM was decomposed within the observation period 1969-1993. The combined decrease in F and humus led to chain reactions with losses of major elements and a dealumination of clay minerals, i.e. removal of interlayered Al of 2:1 minerals and consequent formation of smectites.     

Occurrence of UV filter compounds from sunscreens in surface waters: regional mass balance in two Swiss lakes

Consumer care products often contain UV filters, organic compounds which absorb ultraviolet light. These compounds may enter surface waters directly (when released from the skin during swimming and bathing) or indirectly via wastewater treatment plants (when released during showering or washed from textiles). Predicted and measured UV filter concentrations were compared in a regional mass balance study for two Swiss lakes: Lake Zurich, a typical midland lake which is also an important drinking water resource, and Hüttnersee, a small bathing lake. Both lakes are extensively used for recreational activities and considerable direct input of UV filters is thus expected. This input was estimated from the number of visitors at swimming areas around the lakes and a survey of the usage of sunscreen products among these visitors. Possible additional indirect input via wastewater treatment plants was not considered in this study. The quantitatively most important UV filters, as indicated by the survey data, ethylhexyl methoxycinnamate, octocrylene, 4-methylbenzylidene camphor, butyl methoxydibenzoylmethane, and benzophenone-3, all lipophilic compounds, were selected for analysis by gas chromatography-mass spectrometry. Concentrations of individual UV filters in water from Lake Zurich were low, ranging from <2 ng l(-1) (detection limit) to 29 ng l(-1), and somewhat higher at Hüttnersee, ranging from <2 to 125 ng l(-1), with the highest concentrations found in summer, consistent with direct inputs to the lakes during this time. The concentrations were clearly lower than predicted from input estimates based on the surveys. This may be in part due to (i) an overestimation of these inputs (e.g. less than the 50% wash-off of UV filters assumed to occur during swimming), and (ii) some removal of these compounds from the lakes by degradation and/or sorption/sedimentation. UV filters were also detected in semipermeable membrane devices (SPMDs) deployed at Lake Zurich and Greifensee, another midland lake, at concentrations of 80-950 ng SPMD(-1), confirming the presence of the compounds in surface waters and indicating a certain potential for bioaccumulation. SPMD-derived water concentrations were in the range of 1-10 ng l(-1) and thus corresponded well with those determined in water directly. No UV filters were detected above blank levels in SPMDs deployed at a remote mountain lake used for background measurements.     

Sorption of metolachlor and atrazine in fly ash amended soils: comparison of optimized isotherm models

Adsorption of metolachlor and atrazine was studied in the fly ash (Inderprastha and Badarpur)- amended Inceptisol and Alfisol soils using batch method. Results indicated that sorption of both the herbicides in soil+fly ash mixtures was highly nonlinear and sorption decreased with a higher herbicide concentration in the solution. Also, nonlinearity increased with an increase in the level of fly ash amendment from 0-5%. Three two-parameter monolayer isotherms viz. Langmuir, Temkin, Jovanovic and one imperical Freundlich models were used to fit the experimental data. Data analysis and comparison revealed that the Temkin and the Freundlich isotherms were best-suited to explain the sorption results and the observed and the calculated adsorption coefficient values showed less variability. The study suggested that sorption mechanism of metolachlor and atrazine involved the physical association at the sorbate surface and the nonlinearity in the sorption at higher pesticide or fly ash concentration was due to a decrease in the heat of adsorption and higher binding energy.