Spatial variability in compartmental fate modelling

A new approach is presented which is designed to address the spatial heterogeneity of the environment in compartmental mass balance models of chemical fate in the environment. It rests on the assumption of chemical equilibration within one phase despite prevailing environmental heterogeneity. Composite D- and Z-values are derived from sub-unit specific environmental parameters and are used to solve mass balance equations which can be adopted essentially unchanged from existing compartmental fugacity models. With the resulting common fugacity value for each compartment, sub-unit specific concentrations and process rates can be calculated. The approach is illustrated using the QWASI lake model to calculate the fate of hexachlorobenzene in a hypothetical lake sub-divided in four distinct sub-units. The approach allows the subdivision of each compartment in a large number of sub-units with distinct environmental characteristics without substantially increasing model complexity. This is a necessary condition for linking fugacity models to geographical information systems.     

Implications of geographic variability on Comparative Toxicity Potentials of Cu, Ni and Zn in freshwaters of Canadian ecoregions

Current methods of estimating potential environmental impacts of metals in hazard and Life Cycle Impact Assessment (LCIA) do not consider differences in chemistry and landscape properties between geographic sites. Here, we developed and applied a model for regional aquatic impact characterization of metals using an updated method for estimating environmental fate factor (FF), bioavailability factor (BF) and aquatic ecotoxicity factor (EF). We applied the model to analyze differences in Comparative Toxicity Potentials (CTPs) of Cu, Ni and Zn for 24 Canadian ecoregions. The combined impacts of regional variability in ambient chemistry (in particular DOC, pH and hardness) and landscape properties (water residence time) can change the CTPs of these metals for freshwater by up to three orders of magnitude and change the relative ranking of metal hazard between ecoregions. Variation among Canadian freshwater chemistries and landscape characteristics influence the FFs within two orders of magnitude, BFs within two orders of magnitude for Ni and Zn and four orders of magnitude for Cu, and EFs within one order of magnitude. Sensitivity of metal FFs to environmental parameters alone spans three orders of magnitude when a constant water chemistry was used for all ecoregions. These results indicate that application of regionalised metal CTPs can have a significant influence in the analysis of ecotoxicological impacts in the life cycle assessment of products and processes.     

A multimedia environmental model of chemical distribution: fate, transport, and uncertainty analysis

This paper presented a framework for analysis of chemical concentration in the environment and evaluation of variance propagation within the model. This framework was illustrated through a case study of selected organic compounds of benzo[alpha]pyrene (BAP) and hexachlorobenzene (HCB) in the Great Lakes region. A multimedia environmental fate model was applied to perform stochastic simulations of chemical concentrations in various media. Both uncertainty in chemical properties and variability in hydrometeorological parameters were included in the Monte Carlo simulation, resulting in a distribution of concentrations in each medium. Parameters of compartmental dimensions, densities, emissions, and background concentrations were assumed to be constant in this study. The predicted concentrations in air, surface water and sediment were compared to reported data for validation purpose. Based on rank correlations, a sensitivity analysis was conducted to determine the influence of individual input parameters on the output variance for concentration in each environmental medium and for the basin-wide total mass inventory. Results of model validation indicated that the model predictions were in reasonable agreement with spatial distribution patterns, among the five lake basins, of reported data in the literature. For the chemical and environmental parameters given in this study, parameters associated to air-ground partitioning (such as moisture in surface soil, vapor pressure, and deposition velocity) and chemical distribution in soil solid (such as organic carbon partition coefficient and organic carbon content in root-zone soil) were targeted to reduce the uncertainty in basin-wide mass inventory. This results of sensitivity analysis in this study also indicated that the model sensitivity to an input parameter might be affected by the magnitudes of input parameters defined by the parameter settings in the simulation scenario. Therefore, uncertainty and sensitivity analyses for environmental fate models was suggested to be conducted after the model output was validated based on an appropriate input parameter settings.     

Fugacity superposition: a new approach to dynamic multimedia fate modeling

The fugacities, concentrations, or inventories of pollutants in environmental compartments as determined by multimedia environmental fate models of the Mackay type can be superimposed on each other. This is true for both steady-state (level III) and dynamic (level IV) models. Any problem in multimedia fate models with linear, time-invariant transfer and transformation coefficients can be solved through a superposition of a set of n independent solutions to a set of coupled, homogeneous first-order differential equations, where n is the number of compartments in the model. For initial condition problems in dynamic models, the initial inventories can be separated, e.g. by a compartment. The solution is obtained by adding the single-compartment solutions. For time-varying emissions, a convolution integral is used to superimpose solutions. The advantage of this approach is that the differential equations have to be solved only once. No numeric integration is required. Alternatively, the dynamic model can be simplified to algebraic equations using the Laplace transform. For time-varying emissions, the Laplace transform of the model equations is simply multiplied with the Laplace transform of the emission profile. It is also shown that the time-integrated inventories of the initial conditions problems are the same as the inventories in the steady-state problem. This implies that important properties of pollutants such as potential dose, persistence, and characteristic travel distance can be derived from the steady state.     

Modelling the fate of persistent organic pollutants in Europe: parameterisation of a gridded distribution model

A regionally segmented multimedia fate model for the European continent is described together with an illustrative steady-state case study examining the fate of gamma-HCH (lindane) based on 1998 emission data. The study builds on the regionally segmented BETR North America model structure and describes the regional segmentation and parameterisation for Europe. The European continent is described by a 5 degrees x5 degrees grid, leading to 50 regions together with four perimetric boxes representing regions buffering the European environment. Each zone comprises seven compartments including; upper and lower atmosphere, soil, vegetation, fresh water and sediment and coastal water. Inter-regions flows of air and water are described, exploiting information originating from GIS databases and other georeferenced data. The model is primarily designed to describe the fate of Persistent Organic Pollutants (POPs) within the European environment by examining chemical partitioning and degradation in each region, and inter-region transport either under steady-state conditions or fully dynamically. A test case scenario is presented which examines the fate of estimated spatially resolved atmospheric emissions of lindane throughout Europe within the lower atmosphere and surface soil compartments. In accordance with the predominant wind direction in Europe, the model predicts high concentrations close to the major sources as well as towards Central and Northeast regions. Elevated soil concentrations in Scandinavian soils provide further evidence of the potential of increased scavenging by forests and subsequent accumulation by organic-rich terrestrial surfaces. Initial model predictions have revealed a factor of 5-10 underestimation of lindane concentrations in the atmosphere. This is explained by an underestimation of source strength and/or an underestimation of European background levels. The model presented can further be used to predict deposition fluxes and chemical inventories, and it can also be adapted to provide characteristic travel distances and overall environmental persistence, which can be compared with other long-range transport prediction methods.     

Environmental Persistence of Chemicals

- DOI: http://dx.doi.org/10.1065/espr2006.01.008 Background The hazard criterion of persistence as it applies to chemicals in the environment is reviewed and discussed. This quantity can not be measured directly in the environment, thus it must be estimated using models that synthesise information on chemical half-lives and partitioning properties, the nature of the environment and how the chemical is released into the environment. Main Features It is suggested that the preferred criterion is the average residence time of the chemical in the environment, i.e. conceptually the sum of the life-times of all molecules (attributable only to losses by degrading reactions) divided by the number of molecules. If all chemical fate processes are first order, this persistence is independent of the quantity of chemical introduced and whether introduction is steady- or unsteady-state in nature. It is shown that in a multimedia environment persistence is affected not only by degradation kinetics, but also by mode-of-entry and partitioning. For screening level purposes a Level II equilibrium model may be adequate but a Level III model is generally preferable for estimating the average persistence. If a distribution of persistences is required a dynamic Level IV model must be used. Discussion The implications for regulating chemicals on the basis of persistence are discussed. Conclusion It is concluded that the preferred strategy is to use Level II, III, and IV models and that the use of only degradation kinetics or media-specific half-lives can be misleading and uneconomical.     

Molecular orbital studies on brominated diphenyl ethers. Part I--conformational properties

Polybrominated diphenyl ethers (PBDEs) are widely used as additive flame retardants and quantities in the environment are on the rise. Because they are structurally related to polychlorinated biphenyls and also to thyroid hormones, there is serious concern that PBDEs may pose a danger to human health. Knowledge of their conformational properties is key to assessing their environmental fate and risk. The conformational properties of PBDEs were investigated by quantum chemical methods including semiempirical self-consistent field molecular orbital (SCF-MO), ab initio SCF-MO and density functional theory (DFT). Conformational analyses of model congeners 2,2',4,6'-tetrabromodiphenyl ether and 2,3,4,4',5,6-hexabromodiphenyl ether, based on energy maps calculated by semiempirical AM1 method, may indicate that all PBDE congeners except those with the tetra-ortho-bromination are conformationally flexible (or soft) due to low energy barriers for interconversion of stable conformers. The results of the conformational analyses are in conformity with recently published X-ray crystallographic data. For comparison with the results of the semiempirical method, higher level ab initio and DFT models were applied as well. The optimized geometries all lie well inside low energy regions on the maps and thus also ascertain the semiempirical calculations. According to computed geometric parameters and net atomic charges, the model B3LYP/3-21G* seemed to give better results than B3LYP/6-31G* and HF/6-31G*.     

BETR North America: a regionally segmented multimedia contaminant fate model for North America

We present the Berkeley-Trent North American contaminant fate model (BETR North America), a regionally segmented multimedia contaminant fate model based on the fugacity concept. The model is built on a framework that links contaminant fate models of individual regions, and is generally applicable to large, spatially heterogeneous areas. The North American environment is modeled as 24 ecological regions, within each region contaminant fate is described using a 7 compartment multimedia fugacity model including a vertically segmented atmosphere, freshwater, freshwater sediment, soil, coastal water and vegetation compartments. Inter-regional transport of contaminants in the atmosphere, freshwater and coastal water is described using a database of hydrological and meteorological data compiled with Geographical Information Systems (GIS) techniques. Steady-state and dynamic solutions to the 168 mass balance equations that make up the linked model for North America are discussed, and an illustrative case study of toxaphene transport from the southern United States to the Great Lakes Basin is presented. Regionally segmented models such as BETR North America can provide a critical link between evaluative models of long-range transport potential and contaminant concentrations observed in remote regions. The continent-scale mass balance calculated by the model provides a sound basis for evaluating long-range transport potential of organic pollutants, and formulation of continent-scale management and regulatory strategies for chemicals.     

Modeling habitat suitability for moose in coastal northern Sweden: empirical vs process-oriented approaches

Habitat models may provide viable tools for co-management of large ungulates and forest resources, yet their applicability has not been comprehensively evaluated in managed forest. We examined 2 inherently different approaches to model the relative winter habitat suitability for moose (Alces alces) in the coastal area of northern Sweden. An empirical approach based on GPS positions of 15 female moose was used to scrutinize the assumptions and functional mechanisms of a process-oriented, conceptual approach, based on published material on the species' preferences for habitat components related to food and cover. For both model approaches habitat was described using estimates of forest-stand characteristics based on satellite imagery. The empirical model also included variables relating to topographic properties of the landscape as well as distances to landscape features. The output from both models was a habitat suitability index (HSI) score, enabling the models to be compared with each other. The models showed different results, highlighting the need to include the spatially explicit distribution of environmental variables in future conceptual, process-oriented models.     

BETR-World: a geographically explicit model of chemical fate: application to transport of alpha-HCH to the Arctic

The Berkeley-Trent (BETR)-World model, a 25 compartment, geographically explicit fugacity-based model is described and applied to evaluate the transport of chemicals from temperate source regions to receptor regions (such as the Arctic). The model was parameterized using GIS and an array of digital data on weather, oceans, freshwater, vegetation and geo-political boundaries. This version of the BETR model framework includes modification of atmospheric degradation rates by seasonally variable hydroxyl radical concentrations and temperature. Degradation rates in all other compartments vary with seasonally changing temperature. Deposition to the deep ocean has been included as a loss mechanism. A case study was undertaken for alpha-HCH. Dynamic emission scenarios were estimated for each of the 25 regions. Predicted environmental concentrations showed good agreement with measured values for the northern regions in air, and fresh and oceanic water and with the results from a previous model of global chemical fate. Potential for long-range transport and deposition to the Arctic region was assessed using a Transfer Efficiency combined with estimated emissions. European regions and the Orient including China have a high potential to contribute alpha-HCH contamination in the Arctic due to high rates of emission in these regions despite low Transfer Efficiencies. Sensitivity analyses reveal that the performance and reliability of the model is strongly influenced by parameters controlling degradation rates.     

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.     

Mass balance modelling of contaminants in river basins: application of the flexible matrix approach

It is useful to have available a variety of catchment-scale water quality models that range in complexity, spatial resolution and data requirements. In a previous paper [Warren, C., Mackay, D., Whelan, M., Fox, K., 2005. Mass balance modelling of contaminants in river basins: a flexible matrix approach. Chemosphere 61, 1458-1467] a series of simple to intermediately complex mass balance models was presented which can be used for tiered exposure assessments in river basins. The connectivity of the segments is expressed using a matrix that permits flexibility in application, enabling the model to be re-segmented and applied to different catchments as required. In this paper, the intermediate models, QWASI matrix-rate constant (QMX-R) and QWASI matrix-fugacity (QMX-F) are used to estimate concentrations of linear alkylbenzene sulfonates (LAS) in the rivers Aire and Calder, UK, and of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in the Fraser River basin, Canada. The results compare satisfactorily with monitoring data, suggesting that these QWASI-based models for exposure and risk assessment may be applicable under data-limited conditions. The use of QWASI-based models for regulatory purposes in an evaluative river system is also discussed with reference to assessments of para-dichlorobenzene (pDCB), trichloroethylene (TCE), bis(2-ethylhexyl) phthalate (DEHP) and toluene. It is shown that multi-media QWASI model predictions can be usefully depicted graphically on chemical space diagrams and used to highlight regions in which advection, partitioning to sediments and volatilization may be important determinants of chemical fate in river systems.     

Effects of input uncertainty and variability on the modelled environmental fate of organic pollutants under global climate change scenarios

Global climate change (GCC) is expected to influence the fate, exposure and risks of organic pollutants to wildlife and humans. Multimedia chemical fate models have been previously applied to estimate how GCC affects pollutant concentrations in the environment and biota, but previous studies have not addressed how uncertainty and variability of model inputs affect model predictions. Here, we assess the influence of climate variability and chemical property uncertainty on future projections of environmental fate of six polychlorinated biphenyl congeners under different GCC scenarios using a spreadsheet version of the ChemCAN model and the Crystal Ball® software. Regardless of emission mode, results demonstrate: (i) uncertainty in degradation half-lives dominates the variance of modelled absolute levels of PCB congeners under GCC scenarios; (ii) when the ratios of predictions under GCC to predictions under present day climate are modelled, climate variability dominates the variance of modelled ratios; and (iii) the ratios also indicate a maximum of about a factor of 2 change in the long-term average environmental concentrations due to GCC that is forecasted between present conditions and the period between 2080 and 2099. We conclude that chemical property uncertainty does not preclude assessing relative changes in a GCC scenario compared to a present-day scenario if variance in model outputs due to chemical properties and degradation half-lives can be assumed to cancel out in the two scenarios.     

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.     

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.     

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.