Calculating fugacity

Exploring the distribution and fate of bisphenol A in an aquatic microcosm combined with a multimedia model

Bisphenol A (BPA), a well-known endocrine-disrupting chemical, has garnered significant attention in environmental science and policy. BPA can enter the aquatic environment through different routes, posing potential risks even at a low concentration. In this study, a four-compartment system [water, sediment, biota (zebrafish), and submerged aquatic vegetation (Vallisneria natans)] of a point source continuous discharge microcosm was established to investigate the distribution and fate of BPA in an aquatic microcosm. The fugacity model generated predicted values were highly consistent with those of the experiments. The distribution of BPA in the model indicates that sediment was the dominant sink. The residence time of reaction and advection was 5.8 and 75.2 d, respectively, which showed that BPA was mainly removed from the aquatic microcosm through the reaction in biota (58 %). Sensitivity analysis revealed that emission data were the most influential parameters for the model output. Transfer processes between the water and biota phases had a closer relationship. This study provides technical support for pollution source management and risk assessment for BPA.

Concentration of 226Ra in soil, water, and sediment of active and abandoned quarries in southern Vietnam and environmental risk assessment: experimental and modeling studies

Open-pit mining disturbs the earth's surface, impacts geological characteristics, and releases many pollutants including heavy metals, radionuclides, and poisonous gases into the environment. This study investigated the difference between the activity levels of 226Ra radionuclide in the abandoned quarry region (region A) and the active quarry region (region B). In the surface soil, the mean values of activity concentrations were 40 Bq kg-1 and 55 Bq kg-1 in region A and region B, respectively. The statistical analysis shows that the obtained values of 226Ra concentrations in Region B are higher and more dispersed than those in Region A. For four study quarry lakes, the ranges of activity concentrations in water and sediment were (0-19.4 mBq L-1) and (25-71 Bq kg-1), respectively. The levels of 226Ra in lake water are within the recommended value of 0.74 Bq L-1 of the United States Environmental Protection Agency (USEPA) for radiation safety. This water can be considered for drinking, agricultural irrigation, and fish farming. The Revised Universal Soil Loss Equation (RUSLE) estimated the soil loss due to soil erosion in the abandoned quarries. It was predicted that the soil erosion degrees in the study area were from very low class to moderately high class according to soil loss classification. The Quantitative Water, Air, Sediment Interaction (QWASI) was validated to predict 226Ra activity concentrations in the water and sediment of quarry lakes. Based on Student's t-test, the predicted values of 226Ra activity concentrations for study quarry lakes agreed with the measured values (confidence interval 95%).

Mercury Dynamics in the Sea of Azov: Insights from a Mass Balance Model

The Sea of Azov, an inland shelf sea bounding Ukraine and Russia, experiences the effects of ongoing and legacy pollution. One of the main contaminants of concern is the heavy metal mercury (Hg), which is emitted from the regional coal industry, former Hg refineries, and the historic use of mercury-containing pesticides. The aquatic biome acts both as a major sink and source in this cycle, thus meriting an examination of its environmental fate. This study collated existing Hg data for the SoA and the adjacent region to estimate current Hg influxes and cycling in the ecosystem. The mercury-specific model "Hg Environmental Ratios Multimedia Ecosystem Sources" (HERMES), originally developed for Canadian freshwater lakes, was used to estimate anthropogenic emissions to the sea and regional atmospheric Hg concentrations. The computed water and sediment concentrations (6.8 ng/L and 55.7 ng/g dw, respectively) approximate the reported literature values. The ongoing military conflict will increase environmental pollution in the region, thus further intensifying the existing (legacy) anthropogenic pressures. The results of this study provide a first insight into the environmental Hg cycle of the Sea of Azov ecosystem and underline the need for further emission control and remediation efforts to safeguard environmental quality.

Evaluating the environmental impacts of smart vineyards through the Life Cycle Assessment

This study aimed to assess the environmental effectiveness of vineyards utilising on-site weather stations integrated with a decision support system (DSS), and to identify the critical hotspots in smart farms that have already obtained integrated or organic certification. For this purpose, Life Cycle Assessment (LCA) methodology was applied. The research comprised three smart farms employing on-site weather stations and a traditional farm without advanced technologies, which served as a benchmark. The analysis revealed variations in environmental footprints driven by differences in farm management practices and soil characteristics. The results highlighted that smart farms, in compliance with integrated or organic certifications, focus on reducing inputs such as agrochemicals or water consumption. However, these reductions could shift the environmental burden to other impacts, such as those related to machinery use, which remained the most critical aspect across all vineyards considered. In some smart farms, critical issues involve other aspects, such as irrigation and fertilisation. The lack of awareness about the potential environmental impacts of the adopted technical options could make smart farms more impactful than traditional farms. Interestingly, this study found that solely implementing advanced technologies could fall short of achieving ecological objectives. This study emphasises the significance of utilising LCA as a valuable tool to support farmers in making informed decisions while adopting technological strategies to achieve environmentally sustainable goals.

Comprehensive evaluation of antibiotics pollution the Yangtze River basin, China: Emission, multimedia fate and risk assessment

Antibiotics have attracted global attention because of their potential ecological and health risks. The emission, multimedia fate and risk of 18 selected antibiotics in the entire Yangtze River basin were evaluated by using a level Ⅳ fugacity model. High antibiotic emissions were found in the middle and lower reaches of the Yangtze River basin. The total antibiotic emissions in the Yangtze River basin exceeded 1600 tons per year between 2013 and 2021. The spatial distribution of antibiotics concentration was the upper Yangtze River > middle Yangtze River > lower Yangtze River, which is positively correlated with animal husbandry size in the basin. Temperature and precipitation increases may decrease the antibiotic concentrations in the environment. Transfer fluxes showed that source emission inputs, advection processes, and degradation fluxes contributed more to the total input and output. High ecological risks in the water environment were found in 2018, 2019, 2020, and 2021. The comprehensive health risk assessment through drinking water and fish consumption routes showed that a small part of the Yangtze River basin is at medium risk, and children have a relatively high degree of health risk. This study provides a scientific basis for the pollution control of antibiotics at the basin scale.

A level IV fugacity-based multimedia model based on steady-state particle/gas partitioning theory and its application to study the spatio-temporal trends of PBDEs in atmosphere of northeast China

Particle/gas (P/G) partitioning can significantly affect the environmental behavior of atmospheric pollutants. In this study, we established a large-scale level IV fugacity-based multimedia model (the S-L4MF Model) based on the steady-state P/G partitioning theory. The spatial and temporal trends with the atmospheric contamination of polybrominated diphenyl ethers (PBDEs) in northeastern China under various climate conditions were simulated by the model. There is a reasonable agreement between the simulated and measured gaseous and particulate concentrations of 3 selected PBDE congeners (BDE-47, -99 and -209). For BDE-47, -99 and -209, 91.9 %, 94.8 % and 86.2 % of data points in the evaluation of the spatial trend, whereas 97.4 %, 98.2 % and 91.6 % of data points in the evaluation of the temporal trend, exhibit discrepancies between the modeled and measured data within 1 order of magnitude. The S-L4MF Model performed better than the other model with the same configuration but an equilibrium-state P/G partitioning assumption. The sensitivity and uncertainty analysis indicated that the air temperature and hexadecane-air partition coefficient were the dominant influencing factors on atmospheric concentrations. In addition, the model was successfully applied to study the inter-annual and seasonal variations of gaseous and particulate concentrations of the three PBDEs during 1971-2020 in Harbin, a northeastern Chinese city. Finally, we illustrated the potential to use the model to understand P/G partitioning behavior and the effects of snow and ice on atmospheric concentrations. In summary, the S-L4MF Model provided a powerful and effective tool for studying the environmental behavior of atmospheric organic pollutants, especially in cold regions.

Understanding Interface Exchanges for Assessing Environmental Sorption of Additives from Microplastics: Current Knowledge and Perspectives

Although the impacts of plastic pollution have long been recognized, the presence, pervasiveness, and ecotoxicological consequences of microplastic-i.e., plastic particles < 5 mm-contamination have only been explored over the last decade. Far less focus has been attributed to the role of these materials and, particularly, microplastics, as vectors for a multitude of chemicals, including those (un)intentionally added to plastic products, but also organic pollutants already present in the environment. Owing to the ubiquitous presence of microplastics in all environmental matrices and to the diverse nature of their chemical and physical characteristics, thoroughly understanding the mechanistic uptake/release of these compounds is inherently complex, but necessary in order to better assess the potential impacts of both microplastics and associated chemicals on the environment. Herein, we delve into the known processes and factors affecting these mechanisms. We center the discussion on microplastics and discuss some of the most prominent ecological implications of the sorption of this multitude of chemicals. Moreover, the key limitations of the currently available literature are described and a prospective outlook for the future research on the topic is presented.

Simulation and risk assessment of typical antibiotics in the multi-media environment of the Yangtze River Estuary under tidal effect

Frequent human activities in estuary areas lead to the release of a large number of antibiotics, which poses a great threat to human health. However, there are very limited studies about the influence of the special natural phenomena on the occurrence and migration of antibiotics in the environment. In this study, we simulated the migration and transformation of six typical antibiotics, including oxytetracycline (OTC), tetracycline (TC), norfloxacin (NOR), ofloxacin (OFX), erythromycin (ETM), and amoxicillin (AMOX), in the environmental media from 2011 to 2019 in the Yangtze River Estuary, by using the level III multi-media fugacity model combined with the factor of tides. The simulation results showed that the most antibiotics mainly existed in soil and sediment while erythromycin were found mainly in water. The concentrations of antibiotics in air, freshwater, seawater, groundwater, sediment, and soil were 10^-23-10^-25, 0.1-12 ng/L, 0.02-7 ng/L, 0.02-16 ng/L, 0.1-13 ng/g, and 0.1-15 ng/g respectively. Sensitivity analysis showed that the degradation rate (K_m) and the soil-to-water runoff coefficient (K_l) were important model parameters, indicating that hydrodynamic conditions had a significant impact on the migration of antibiotics in various environmental phases in estuarine areas. Tide can enhance the exchange between water bodies and cause the transformation of the antibiotics from freshwater to seawater and groundwater, which improved the accuracy of the model, especially the seawater and soil phase. Risk assessments showed that amoxicillin, erythromycin, ofloxacin, and norfloxacin posed a threat to the estuarine environment, but the current source of drinking water did not affect human health. Our findings suggested that, when one would like to exam the occurrence and migration of antibiotics in environment, more consideration should be given to the natural phenomena, in addition to human activities and the nature of antibiotics.

A new and integrated approach to evaluate the environmental and ecotoxicological impact of herbicide mixtures: A case study in maize

An index-based approach for a comprehensive evaluation of the potential risk for active substances and their mixtures to impact the environment was developed. Some of the indices considered already exist (PRISW-1, Priority Index), while others were created ex novo from indicators available on open-source platforms (PESTi, ECOi, AGROi). These indices maybe used for an evaluation before use of pesticides by farmers and advisers. The present approach was initially validated for herbicides in maize crops, but it can readily be applied to other PPPs and crops. PESTi index underline the physical and chemical characteristics as a whole, not considering the impact of other factors such as application rate or period of application. Hence, this index may underestimate the risk associated to a certain chemical. AGROi has a precautionary approach. The risk associated to a specific mixture derives from a combination of intrinsic characteristics of the chemicals, agronomic impacts, regulation restrictions and potential hazard to water compartment. The ECOi index is focused on the ecotoxicological impact against non-target organisms. The helpfulness of this index stands in its ability to easily discriminate the ecotoxicological impact of chemicals using indicators commonly available in literature and without making complex calculations. PRISW-1 Index discriminate active substances according to their risk against three representative non-target organisms. However, due to the intrinsic characteristics of each pesticide, a high PRISW-1 value could not always mean an easy movement of the chemical via runoff waters. The information deriving from Priority index may certainly help public authorities to select chemicals to be detected in water monitoring campaigns. The application of these indices may represent a valid decision tool for public stakeholders in defining agricultural measures to reduce the externalities of pest control.

Source, fate, transport and modelling of selected emerging contaminants in the aquatic environment: Current status and future perspectives

The occurrence of emerging contaminants (ECs), such as pharmaceuticals and personal care products (PPCPs), perfluoroalkyl and polyfluoroalkyl substances (PFASs) and endocrine-disrupting chemicals (EDCs) in aquatic environments represent a major threat to water resources due to their potential risks to the ecosystem and humans even at trace levels. Mathematical modelling can be a useful tool as a comprehensive approach to study their fate and transport in natural waters. However, modelling studies of the occurrence, fate and transport of ECs in aquatic environments have generally received far less attention than the more widespread field and laboratory studies. In this study, we reviewed the current status of modelling ECs based on selected representative ECs, including their sources, fate and various mechanisms as well as their interactions with the surrounding environments in aquatic ecosystems, and explore future development and perspectives in this area. Most importantly, the principles, mathematical derivations, ongoing development and applications of various ECs models in different geographical regions are critically reviewed and discussed. The recommendations for improving data quality, monitoring planning, model development and applications were also suggested. The outcomes of this review can lay down a future framework in developing a comprehensive ECs modelling approach to help researchers and policymakers effectively manage water resources impacted by rising levels of ECs.

Moving persistence assessments into the 21st century: A role for weight-of-evidence and overall persistence

Assessing the persistence of chemicals in the environment is a key element in existing regulatory frameworks to protect human health and ecosystems. Persistence in the environment depends on many fate processes, including abiotic and biotic transformations and physical partitioning, which depend on substances' physicochemical properties and environmental conditions. A main challenge in persistence assessment is that existing frameworks rely on simplistic and reductionist evaluation schemes that may lead substances to be falsely assessed as persistent or the other way around-to be falsely assessed as nonpersistent. Those evaluation schemes typically assess persistence against degradation half-lives determined in single-compartment simulation tests or against degradation levels measured in stringent screening tests. Most of the available test methods, however, do not apply to all types of substances, especially substances that are poorly soluble, complex in composition, highly sorptive, or volatile. In addition, the currently applied half-life criteria are derived mainly from a few legacy persistent organic pollutants, which do not represent the large diversity of substances entering the environment. Persistence assessment would undoubtedly benefit from the development of more flexible and holistic evaluation schemes including new concepts and methods. A weight-of-evidence (WoE) approach incorporating multiple influencing factors is needed to account for chemical fate and transformation in the whole environment so as to assess overall persistence. The present paper's aim is to begin to develop an integrated assessment framework that combines multimedia approaches to organize and interpret data using a clear WoE approach to allow for a more consistent, transparent, and thorough assessment of persistence. Integr Environ Assess Manag 2022;18:868-887. © 2021 ExxonMobil Biomedical Sciences, Inc. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Forests effects on the environmental fates of organic pollutants in a tropical watershed

Poly-brominated diphenyl ethers (PBDEs), dioxins, furans and current-use pesticides (CUPs) are common organic pollutants that have received global scrutiny due to their association with adverse environmental and health impacts. However, there is limited previous work assessing their environmental fates in the context of tropical multi-use watersheds. The aim of this study was to estimate the effect of forests on long-term environmental fate for some PBDEs, dioxins, furans and CUPs for a case study watershed of a tropical developing country (the Rio Cobre River drainage basin, Jamaica). Specifically, a dynamic, 16-compartment environmental multimedia model, RioShed, was developed and applied to calculate compartmental concentrations, as well as some long-term environmental fate metrics. Results indicate that the presence of tropical forests, especially those that are evergreen, reduced the atmospheric concentrations, atmospheric long-range transport potential, and the overall persistence of the study pollutants. Reductions in atmospheric concentrations by tropical forests were most enhanced for the more polar CUPs. Forest parameters that notably influenced soil concentration and/or overall persistence included the canopy drip parameter, the leaf area index and the wax erosion rate. The results of this research are expected to inform land-use and environmental management of the study area and similar tropical regions.

Particle/gas partitioning for semi-volatile organic compounds (SVOCs) in Level III multimedia fugacity models: Gaseous emissions

Atmospheric transport is a global-scale process that moves semi-volatile organic compounds (SVOCs) rapidly from source regions to remote locations, where these chemicals have never been produced or used. Particle/gas (P/G) partitioning of SVOCs during atmospheric transport governs wet and dry deposition, and thereby controls the efficiency and scope of long-range atmospheric transport and fate for these sorts of compounds. Previous work has shown that the assumption of steady state between particulate and gaseous phases in the atmosphere leads to model results that more closely match observations especially for compounds that strongly favor the particulate phase. Here, the practical application of steady-state P/G partitioning in the atmosphere in multimedia fugacity models is presented in greater detail. A method is developed whereby the fugacity of a chemical in the particle-phase is set equal to that in the gaseous phase (a pseudo equilibrium) but still maintains steady state of the chemical between air and aerosols in the atmosphere. This procedure greatly simplifies the application of multimedia fugacity models. Using this approach, a condition of steady state between air and aerosols is developed and applied in a Level III six-compartment six-fugacity model, which becomes a much simpler Level III six-compartment four-fugacity model. This newly-developed model is then applied to data observed during a monitoring program.

Particle/gas partitioning for semi-volatile organic compounds (SVOCs) in level III multimedia fugacity models: Both gaseous and particulate emissions

Multimedia fugacity models have long been used to address the fate of toxic organic chemical emissions by providing a quantitative account of the sources, transport processes, and sinks. Recently, we have examined three level-III fugacity models (E4F (equilibrium six-compartment four-fugacity), S6F (steady-state six-compartment six-fugacity) and S4F (steady-state six-compartment four-fugacity) Models), in the context of their performance set against real-world data, and their practicality of application. Here, we discuss how the balance between gaseous and aerosol phases of emissions assumed for initial conditions affects the different model outcomes. Our results show that the S6F Model predictions closely match those of the S4F Model when chemical emissions are entirely in the gas-phase. As the particulate proportion of the emission increases, the S6F Model predictions diverge from those of the S4F Model and approach those of the E4F Model. Once the particulate portion reaches 100%, the S6F and E4F Models produce identical results: an internally inconsistent system where chemicals are not in a steady state between air and aerosols, and mass balance for both air and aerosols is not achieved. Thus, in terms of practicality, internal consistency, chemical mass balance and agreement with observations, the S4F Model is clearly the best choice.

Regional variation in mercury bioaccumulation among NW Atlantic Golden (Lopholatilus chamaeleonticeps) and Blueline (Caulolatilus microps) Tilefish

Mercury (Hg) concentrations in fishes from the NW Atlantic Ocean pose concern due to the importance of this region to U.S. fisheries harvest. In this study, total Hg (THg) concentrations and nitrogen stable isotope (δ¹⁵N) values were quantified in muscle tissues sampled from Golden (Lopholatilus chamaeleonticeps) and Blueline (Caulolatilus microps) Tilefish collected during a fishery-independent survey conducted in the NW Atlantic to compare bioaccumulation patterns between these species. Total Hg concentrations averaged (±SD) 0.4 ± 0.4 μg/g dry weight (d.w.) for L. chamaeleonticeps and 1.1 ± 0.7 μg/g d.w. for C. microps with <2% of all sampled fish, those >70 cm fork length, exceeding the most restrictive USEPA regulatory guidelines for human consumption (THg > 0.46 μg/g w.w.), when converted to wet weight concentrations. The THg concentrations reported here for individuals from the NW Atlantic stock are comparable to those reported for similarly sized individuals collected from the SW Atlantic stock but notably lower than those reported for Gulf of Mexico L. chamaeleonticeps, indicating different Hg exposure and assimilation kinetics for fish from the NW Atlantic, and highlights the broad geographic variability of Hg bioaccumulation among Tilefish stocks. Caulolatilus microps had higher δ¹⁵N values relative to L. chamaeleonticeps and a pattern of decreasing THg concentrations was also present from south to north across the study range. It is concluded that this trophic difference and spatial pattern in Tilefish THg concentrations emphasizes the habitat and resource partitioning mechanisms described for these sympatric species that permits their coexistence in the continental shelf environment. Importantly, regional variability in THg concentrations accentuate the possible roles of fine-scale biotic and abiotic processes that can act to regulate Hg bioaccumulation among individuals and species.

Treatment of particle/gas partitioning using level III fugacity models in a six-compartment system

In this paper, two level III fugacity models are developed and applied using an environmental system containing six compartments, including air, aerosols, soil, water, suspended particulate matters (SPMs), and sediments, as a "unit world". The first model, assumes equilibrium between air and aerosols and between water and SPMs. These assumptions lead to a four-fugacity model. The second model removes these two assumptions leading to a six-fugacity model. The two models, compared using four PBDE congeners, BDE-28, -99, -153, and -209, with a steady flux of gaseous congeners entering the air, lead to the following conclusions. 1. When the octanol-air partition coefficient (K_OA) is less than 10^11.4, the two models produce similar results; when K_OA > 10^11.4, and especially when K_OA > 10^12.5, the model results diverge significantly. 2. Chemicals are in an imposed equilibrium in the four-fugacity model, but in a steady state and not necessary an equilibrium in the six-fugacity model, between air and aerosols. 3. The results from the six-fugacity model indicate an internally consistent system with chemicals in steady state in all six compartments, whereas the four-fugacity model presents an internally inconsistent system where chemicals are in equilibrium but not a steady state between air and aerosols. 4. Chemicals are mass balanced in air and aerosols predicted by the six-fugacity model but not by the four-fugacity model. If the mass balance in air and aerosols is achieved in the four-fugacity model, the condition of equilibrium between air and aerosols will be no longer valid.

Environmental Partitioning, Spatial Distribution, and Transport of Atmospheric Mercury (Hg) Originating from a Site of Former Chlor-Alkali Plant

Mercury (Hg) is one of the trace toxic and bioaccumulative global pollutants, and due to its long atmospheric lifetime, it presents a significant global challenge. The present study (1) utilizes total gaseous mercury (TGM) measurements made around a former Hg-cell chlor-alkali plant (CAP) located in Pavlodar, Kazakhstan, and predicts the spatial distribution of Hg over its premises and the nearby city. It then (2) estimates the environmental repartition of Hg deposited by the CAP using three fugacity models of varying complexity: Level I, QWASI, and HERMES. Finally, it (3) predicts long-range Hg transport via forward trajectory-based cluster analysis. The atmospheric Hg levels measured in Pavlodar and around Lake Balkyldak were elevated: in the range of 1–37 ng/m3 with an urban background level at 4.9 ng/m3. Specifically, concentrations up to 37 ng/m3 close to Lake Balkyldak and up to 22 ng/m3 nearby the city’s industrial zone (where the CAP was located) had been observed. Interpolation maps created using kriging also suggest these locations as the primary sources of atmospheric Hg in the city. The Level I fugacity model indicated that almost all of Hg is expected to end up in the atmosphere. The modeling results obtained using more complex QWASI and HERMES models showed that some significant quantity of Hg would still be associated with the sediments of Lake Balkyldak (a large wastewater discharge pond nearby the CAP). The forward trajectory-based cluster analysis method revealed the long-range atmospheric transportation routes and local, regional, and global impact zones. Furthermore, a source-receptor relationship using air transportation pathways to identify “areas of impact” was addressed. During both heating and non-heating seasons, the frequency-based analysis identified the distribution of Hg reaching the territories of Mongolia, northwest China, southwest Kazakhstan. The Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT-4) model forward trajectory analysis has confirmed similar patterns during heating and non-heating seasons, except with shorter impact distances during the non-heating period. Even though the CAP was closed more than 30 years ago and those past remediation efforts cleaned up the site, the residual Hg pollution seems significant and should be further investigated in different environmental media.

Ecological Risk Assessment of Oil Spills in Ice-Covered Waters: A Surface Slick Model Coupled with a Food-Web Bioaccumulation Model

The limited knowledge on oil-ice interactions and on the ecological outcomes of oil spills in the Arctic represent sources of uncertainties for shipping and oil and gas activities in polar regions. The present work aims at the definition of the ecological risk posed by oil spills in the Arctic by the integration of an improved surface slick model to a fugacity-based food-web bioaccumulation model for icy waters. The model's outcomes are the representation of transport and weathering processes and the concentrations of a toxic component of oil, namely naphthalene, in the environmental media as a function of ice conditions. Given those concentrations, the associated ecological risk is defined in terms of the bioconcentration factor (BCF). Overall, the model predicted low bioaccumulation and biomagnification potential for naphthalene to a hypothetical Arctic food web, regardless of the ice concentration. Integr Environ Assess Manag 2020;16:729-744. © 2020 SETAC.

ChemFate: A fate and transport modeling framework for evaluating radically different chemicals under comparable conditions

With the ever-increasing development on novel chemicals and materials, with radically different properties and environmental behavior, it is challenging to compare their environmental behavior under similar conditions. For example, pesticides may be non-ionizable organics, ionizable organics, metal ion-based, or nanomaterials. These very different chemicals behave very differently. To date, no single modeling platform can handle all adequately, in a dynamic framework that accounts for actual variability in meteorology, rates of release to the environment, and the actual characteristics of the region of interest. Here we present ChemFate, a framework to address this challenge, by incorporating four different fate and transport models, each to address these four classes of chemicals (non-ionizable organics, ionizable organics, metal ion-based, or nanomaterials). We build upon established models, but have incorporated a number of additional processes. After demonstrating that the individual models comparable favorably with observed data and the previous models, under similar conditions, we conducted a case study with four radically different fungicides, used in the Central Valley, California. We found that although the concentrations of the non-ionizable and ionizable organic fungicide spike after application and temporarily may exceed toxicity thresholds for Daphnia Magna, they do not accumulate over time, while the metal ion and metallic nanoparticle result in increasing accumulation of Cu²⁺, eventually exceeding the toxicity threshold during runoff events. This case study demonstrates the value of a framework that allows the comparison of different classes of chemicals under the same conditions.

A review of multimedia transport and fate models for chemicals: Principles, features and applicability

The frequent use of chemicals has caused ecosystems and humans to be threatened due to their discharge into the environment. Multimedia environmental fate models could provide a comprehensive picture of transport behaviour and fate for organic chemicals in multiple environmental media. They have been designed and widely used for chemical risk assessment, chemical ranking and management support, and determination of chemical bioaccumulation. This study reviewed the principles, features and applicability of recent commonly used multimedia fate models from peer-reviewed literature. Fugacity-based and concentration-based models are now widely adopted for use in chemical fate evaluation, while they are more appropriate for volatile and semi-volatile chemicals. Or the fugacity-based models can use aquivalence equilibrium criterion to cations, anions and involatile chemicals. The MAMI and SESAMe models based on activity approach are applicable to neutral and ionizable molecules. However, interactions of ionic species with other water solutes are not taken into account in these models. Additionally, they could not directionally simulate how chemicals transported form one grid to another. Future attention should be focused on the reliability of transfer behaviour and fate of ionizable chemicals, as integrating the advantages of these two kinds of models into a reconstructed one may be a better choice. In a word, environmental multimedia models have been beneficial tools for chemical control and management, risk and effect estimation, and decision supporting.

Protein and lipid growth rates regulate bioaccumulation of PCBs and Hg in Bighead Carp (Hypophthalmichthys nobilis) and Silver Carp (Hypophthalmichthys molitrix) from the Three Gorges Reservoir, China

This study evaluated the effect of growth of different tissue compartments on the bioaccumulation of mercury (Hg) and polychlorinated biphenyls (PCBs) in Silver Carp (Hypophthalmichthys molitrix) and Bighead Carp (Hypophthalmichthys nobilis) from the Three Gorges Reservoir (TGR), China. A non-steady state bioenergetics/toxicokinetic model was developed to simulate PCB and Hg concentrations in these two species and compared with field data. Simulations using constant whole body growth rate and constant tissue to whole body weight ratios were contrasted against simulations adopting age specific whole body and tissue/age specific growth rates for their goodness of fit to field data. The simulations using age/tissue specific growth rates demonstrated better fit to field data for PCBs compared to the constant growth rate models (22% improved R²), while both models explained similar variation in Hg concentration data. Both species demonstrated higher growth rates of lipids (on a daily basis) relative to whole body and protein contributing to higher growth dilution of PCBs compared to Hg. Although stable isotope data indicated some degree of diet and/or habitat shift, simulations assuming a constant diet concentration explained between 36 and 40% of the variation in fish concentrations for both contaminants and fish species. This study demonstrates that differences in the bioaccumulation rate of PCBs and Hg by Asian carp can be partially explained by differences in the growth rates of key tissue storage compartments associated with each contaminant. These differences in chemical-specific growth dilution subsequently contribute to differences in chemical retention and bioaccumulation patterns of Hg and PCBs by fish.

Seasonal soil/snow-air exchange of semivolatile organic pollutants at a coastal arctic site (Tromsø, 69°N)

Soils are a major reservoir of semivolatile organic pollutants such as polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons (PAHs), and exert a control on their atmospheric occurrence. We present here an assessment of the atmospheric occurrence and seasonality of soil/snow-air partitioning and exchange of PCBs, PAHs, hexachlorobenzene (HCB), and hexachlorocyclohexanes (HCHs) in the arctic city Tromsø, northern Norway. The fugacities of the organic pollutants in soils and snow were determined using a soil fugacity sampler by equilibrating the air concentrations with those in the surface soil/snow. The concentrations in soils did not show a significant seasonality. Conversely, the ambient air concentrations and the soil (or snow) fugacity showed a clear seasonality for PCBs, HCH, HCB and some PAHs, related to temperature. Fugacities in soil/snow were correlated with those in the ambient gas phase, suggesting a close seasonal air-soil/snow coupling. Generally, there was a net deposition or close to equilibrium conditions during the winter, which contrasts with the net volatilization observed during the warmer periods. The chemicals with lower octanol-air partition coefficients showed a larger tendency for being volatilized and thus remobilized from this coastal arctic environment. Conversely, the more hydrophobic compounds were close to air-soil/snow equilibrium or showed a net deposition.

Petroleum hydrocarbon mixture toxicity and a trait-based approach to soil invertebrate species for site-specific risk assessments

Although petroleum hydrocarbons released to the environment typically occur as mixtures, petroleum hydrocarbon remediation guidelines often reflect individual substance toxicity. It is well documented that groups of aliphatic petroleum hydrocarbons act via the same mechanism of action (nonpolar narcosis) and, theoretically, concentration addition mixture toxicity principles apply. To assess this theory, 10 standardized acute and chronic soil invertebrate toxicity tests on a range of organisms (Eisenia fetida, Lumbricus terrestris, Enchytraeus crypticus, Folsomia candida, Oppia nitens, and Hypoaspis aculeifer) were conducted with a refined petroleum hydrocarbon binary mixture. Reference models for concentration addition and independent action were applied to the mixture toxicity data with consideration of synergism, antagonism, and dose level toxicity. Both concentration addition and independent action, without further interactions, provided the best fit with observed response to the mixture. Individual fraction effective concentration values were predicted from optimized, fitted reference models. Concentration addition provided a better estimate than independent action of individual fraction effective concentrations based on comparison with available literature and species trends observed in toxic responses to the mixture. Interspecies differences in standardized laboratory soil invertebrate species responses to petroleum hydrocarbon-contaminated soil was reflected in unique traits. Diets that included soil, large body size, permeable cuticle, low lipid content, lack of ability to molt, and no maternal transfer were traits linked to a sensitive survival response to petroleum hydrocarbon-contaminated soil in laboratory tests. Traits linked to sensitive reproduction response in organisms tested were long life span and small clutch size. By deriving single-fraction toxicity endpoints considerate of mixtures, we can reduce the resources and time required to conduct site-specific risk assessments for the protection of a soil organism's exposure pathway. Environ Toxicol Chem 2018;37:2222-2234. © 2018 SETAC.

Persistent organic pollutants in the Atlantic and southern oceans and oceanic atmosphere

Persistent organic pollutants (POPs) continue to cycle through the atmosphere and hydrosphere despite banned or severely restricted usages. Global scale analyses of POPs are challenging, but knowledge of the current distribution of these compounds is needed to understand the movement and long-term consequences of their global use. In the current study, air and seawater samples were collected Oct. 2007-Jan. 2008 aboard the Icebreaker Oden en route from Göteborg, Sweden to McMurdo Station, Antarctica. Both air and surface seawater samples consistently contained α-hexachlorocyclohexane (α-HCH), γ-HCH, hexachlorobenzene (HCB), α-Endosulfan, and polychlorinated biphenyls (PCBs). Sample concentrations for most POPs in air were higher in the northern hemisphere with the exception of HCB, which had high gas phase concentrations in the northern and southern latitudes and low concentrations near the equator. South Atlantic and Southern Ocean seawater had a high ratio of α-HCH to γ-HCH, indicating persisting levels from technical grade sources. The Atlantic and Southern Ocean continue to be net sinks for atmospheric α-, γ-HCH, and Endosulfan despite declining usage.

Watershed-scale modeling on the fate and transport of polycyclic aromatic hydrocarbons (PAHs)

PAHs are potentially carcinogenic substances that are persistent in the environment. Increasing concentrations of PAHs were observed due to rapid urbanization, thus; monitoring PAHs concentrations is necessary. However, it is expensive to conduct intensive monitoring activities of a large number of PAHs. This study addressed this issue by developing a multimedia model coupled with a hydrological model (i.e., Soil and Water Assessment Tool (SWAT)) for Taehwa River (TR) watershed in Ulsan, the industrial capital of South Korea. The hydrologic module of the SWAT was calibrated, and further used to simulate the fate and transport of PAHs in soil and waterbody. The model demonstrated that the temporal or seasonal variation of PAHs in soil and waterbody can be well reproduced. Meanwhile, the spatial distribution of PAHs showed that urban areas in TR watershed have the highest PAH loadings compared to rural areas. Sensitivity analyses of the PAH soil and PAH water parameters were also able to determine the critical processes in TR watershed: degradation, deposition, volatilization, and wash off mechanism. We hope that this model will be able to aid the stakeholders in: regulating PAH concentrations emitted by various sources; and also apply the model to other Persistent Organic Pollutants (POPs).

On the formulation of environmental fugacity models and their numerical solutions

Multimedia models based on chemical fugacity, solved numerically, play an important role in investigating and quantifying the environmental fate of chemicals such as persistent organic pollutants. These models have been used extensively in studying the local and global distribution of chemicals in the environment. The present study describes potential sources of error that may arise from the formulation and numerical solution of environmental fugacity models. The authors derive a general fugacity equation for the rate of change of mass in an arbitrary volume (e.g., an environmental phase). Deriving this general equation makes clear several assumptions that are often not articulated but can be important for successfully applying multimedia fugacity models. It shows that the homogeneity of fugacity and fugacity capacity in a volume (the homogeneity assumption) is fundamental to formulating discretized fugacity models. It also shows that when using the fugacity rather than mass as the state-variable, correction terms may be necessary to accommodate environmental factors such as varying phase temperatures and volume. Neglecting these can lead to conservation errors. The authors illustrate the manifestation of these errors using heuristic multimedia fugacity models. The authors also show that there are easily avoided errors that can arise in mass state-variable models if variables are not updated appropriately in the numerical integration scheme. Environ Toxicol Chem 2016;35:2182-2191. © 2016 SETAC.

The environmental fate of polybrominated diphenyl ethers (PBDEs) in western Taiwan and coastal waters: evaluation with a fugacity-based model

The environmental fate of polybrominated diphenyl ethers (PBDEs), a group of flame retardants that are considered to be persistent organic pollutants (POPs), around the Zhuoshui River and Changhua County regions of Taiwan was assessed. An investigation into emissions, partitioning, and fate of selected PBDEs was conducted based on the equilibrium constant (EQC) fugacity model developed at Trent University, Canada. Emissions for congeners PBDE 47, PBDE 99, and PBDE 209 to air (4.9-92 × 10(-3) kg/h), soil (0.91-17.4 × 10(-3) kg/h), and water (0.21-4.04 × 10(-3) kg/h), were estimated by modifying previous models on PBDE emission rates by considering both industrial and domestic rates. It was found that fugacity modeling can give a reasonable estimation of the behavior, partitioning, and concentrations of PBDE congeners in and around Taiwan. Results indicate that PBDE congeners have a high affinity for partitioning into sediments then soils. As congener number decreases, the PBDEs then partition more readily into air. As the degree of bromination increases, congeners more readily partition to sediments. Sediments may then act as a long-term source of PBDEs which can be released back into the water column due to resuspension during storm events.

What determines PCB concentrations in soils in rural and urban areas? Insights from a multi-media fate model for Switzerland as a case study

Polychlorinated biphenyls (PCBs) are banned worldwide under the Stockholm Convention on Persistent Organic Pollutants. However, PCBs are still emitted in appreciable amounts from remaining primary sources in urban areas or landfills and are ubiquitous environmental contaminants, inter alia in soil and air. Concentrations of PCBs have been measured in various media by numerous studies worldwide. However, monitoring data do not always provide quantitative information about transport processes between different media, deposition fluxes to ground, or distribution of PCBs between environmental compartments. Also future trends in environmental contamination by PCBs cannot be predicted from monitoring data, but such information is highly relevant for decision-makers. Here, we present a new regionally resolved dynamic multimedia mass balance model for Switzerland to investigate the origin of PCBs in air and to investigate their long-term fate and mass balance in the environment. The model was validated with existing field data for PCBs. We find that advective inflow of PCBs from outside Switzerland into the atmospheric boundary layer is responsible for 80% of PCBs present in air in Switzerland, whereas Swiss emissions cause the remaining 20%. Furthermore, we show that the atmospheric deposition of the higher-chlorinated PCBs is dominated by particle-bound deposition, whereas the deposition of the lower-chlorinated PCBs is a combination of particle-bound and gaseous deposition. The volume fraction of particles in air is in both cases an important factor driving the deposition of PCBs to ground and, thus, contributing to the higher concentrations of PCBs generally observed in populated and polluted areas. Regional emissions influence the deposition fluxes only to a limited extent. We also find that secondary emissions from environmental reservoirs do not exceed primary emissions for all PCB congeners until at least 2036. Finally, we use our model to evaluate the effect of chemical regulation on future environmental contamination by PCBs.

Aquatic Exposure Predictions of Insecticide Field Concentrations Using a Multimedia Mass-Balance Model

Highly complex process-driven mechanistic fate and transport models and multimedia mass balance models can be used for the exposure prediction of pesticides in different environmental compartments. Generally, both types of models differ in spatial and temporal resolution. Process-driven mechanistic fate models are very complex, and calculations are time-intensive. This type of model is currently used within the European regulatory pesticide registration (FOCUS). Multimedia mass-balance models require fewer input parameters to calculate concentration ranges and the partitioning between different environmental media. In this study, we used the fugacity-based small-region model (SRM) to calculate predicted environmental concentrations (PEC) for 466 cases of insecticide field concentrations measured in European surface waters. We were able to show that the PECs of the multimedia model are more protective in comparison to FOCUS. In addition, our results show that the multimedia model results have a higher predictive power to simulate varying field concentrations at a higher level of field relevance. The adaptation of the model scenario to actual field conditions suggests that the performance of the SRM increases when worst-case conditions are replaced by real field data. Therefore, this study shows that a less complex modeling approach than that used in the regulatory risk assessment exhibits a higher level of protectiveness and predictiveness and that there is a need to develop and evaluate new ecologically relevant scenarios in the context of pesticide exposure modeling.

SimpleBox 4.0: Improving the model while keeping it simple…

Chemical behavior in the environment is often modeled with multimedia fate models. SimpleBox is one often-used multimedia fate model, firstly developed in 1986. Since then, two updated versions were published. Based on recent scientific developments and experience with SimpleBox 3.0, a new version of SimpleBox was developed and is made public here: SimpleBox 4.0. In this new model, eight major changes were implemented: removal of the local scale and vegetation compartments, addition of lake compartments and deep ocean compartments (including the thermohaline circulation), implementation of intermittent rain instead of drizzle and of depth dependent soil concentrations, adjustment of the partitioning behavior for organic acids and bases as well as of the value for enthalpy of vaporization. In this paper, the effects of the model changes in SimpleBox 4.0 on the predicted steady-state concentrations of chemical substances were explored for different substance groups (neutral organic substances, acids, bases, metals) in a standard emission scenario. In general, the largest differences between the predicted concentrations in the new and the old model are caused by the implementation of layered ocean compartments. Undesirable high model complexity caused by vegetation compartments and a local scale were removed to enlarge the simplicity and user friendliness of the model.

Use of a vial equilibration technique to measure the change in fugacity capacity of avian food and feces samples for 1,2,3,4-tetrachlorobenzene

A vial equilibration technique was used to estimate the fugacity capacities of food and feces samples for 1,2,3,4-tetrachlorobenzene (TCB). The method was calibrated using different volumes of n-octanol and by comparing the measured and predicted fugacity capacity (Zoct) of n-octanol for TCB. The vial equilibration technique showed linearity with increasing amounts of n-octanol added to the vial. However, the measured Zoct was on average 8.5 times lower than the literature estimate and interpreted to be influenced by co-solvent effects. The ratio of fugacity capacities of food/feces was 2.9 and was consistent with the ratio estimated using Zt calculation methods (4.3) which considers partitioning capacity of both lipids and non-lipid organic matter. These results provide experimental support to the use of lipid equivalent approaches as opposed to lipid normalization when estimating the partition capacity of biological samples containing low lipid contents.

Field determination and QSPR prediction of equilibrium-status soil/vegetation partition coefficient of PCDD/Fs

Characterizing pseudo equilibrium-status soil/vegetation partition coefficient KSV, the quotient of respective concentrations in soil and vegetation of a certain substance at remote background areas, is essential in ecological risk assessment, however few previous attempts have been made for field determination and developing validated and reproducible structure-based estimates. In this study, KSV was calculated based on measurements of seventeen 2,3,7,8-substituted PCDD/F congeners in soil and moss (Dicranum angustum), and rouzi grass (Thylacospermum caespitosum) of two background sites, Ny-Ålesund of the Arctic and Zhangmu-Nyalam region of the Tibet Plateau, respectively. By both fugacity modeling and stepwise regression of field data, the air-water partition coefficient (KAW) and aqueous solubility (SW) were identified as the influential physicochemical properties. Furthermore, validated quantitative structure-property relationship (QSPR) model was developed to extrapolate the KSV prediction to all 210 PCDD/F congeners. Molecular polarizability, molecular size and molecular energy demonstrated leading effects on KSV.

Bioamplification as a bioaccumulation mechanism for persistent organic pollutants (POPs) in wildlife

Persistent organic pollutant bioaccumulation models have generally been formulated to predict bioconcentration and biomagnification. A third bioaccumulation process that can mediate chemical fugacity in an organism is bioamplification.Bioamplification occurs when an organism loses body weight and the chemical partitioning capacity occurs at a rate that is faster than the chemical can be eliminated.Although bioamplification has not been widely recognized as a bioaccumulation process, the potential consequences of this process are significant. Bioamplification causes an increase in chemical fugacity in the animal's tissues and results in there distribution of contaminants from inert storage sites to more toxicologically sensitive tissues. By reviewing laboratory and field studies, we have shown in this paper that bioamplification occurs across taxonomic groups that include, invertebrates,amphibians, fishes, birds, and mammals. Two case studies are presented, and constitute multi-life stage non-steady state bioaccumulation models calibrated for yellow perch and herring gulls. These case studies were used to demonstrate that bioamplification is predicted to occur under realistic scenarios of animal growth and seasonal weight loss. Bioamplification greatly enhances POP concentrations and chemical fugacities during critical physiological and behavioral events in an animal's life history, e.g., embryo development, juvenile stages, metamorphosis, reproduction, migration, overwintering, hibernation, and disease. Consequently,understanding the dynamics of bioamplification, and how different life history scenario scan alter tissue residues, may be helpful and important in assessing wildlife hazards and risks.

Polychlorinated biphenyl concentrations and activity of sea lamprey Petromyzon marinus vary by sex

We determined the polychlorinated biphenyl (PCB) concentrations of 40 male and 40 female adult sea lampreys Petromyzon marinus captured in the Cheboygan River, a tributary to Lake Huron, during May 2011. In addition, we performed a laboratory experiment using passive integrated transponder tags to determine whether male adult sea lampreys were more active than female adult sea lampreys. Sex had a significant effect on PCB concentration, and PCB concentration at a given level of sea lamprey condition was approximately 25 % greater in males than in females. Adjusting for the difference in condition between the sexes, males averaged a 17 % greater PCB concentration compared with females. Results from the laboratory experiment indicated that males were significantly more active than females. The observed sex difference in PCB concentrations was not due to female sea lampreys releasing eggs at spawning because the sea lamprey is semelparous, and we caught the sea lampreys before spawning. Rather, we attributed the sex difference in PCB concentrations to a greater rate of energy expenditure in males compared with females. We proposed that this greater rate of energy expenditure was likely due to greater activity. Our laboratory experiment results supported this hypothesis. A greater resting metabolic rate may also have contributed to a greater rate of energy expenditure. Our findings should eventually be applicable toward improving control of sea lamprey, a pest responsible for considerable damage to fisheries in lakes where it is not native.

Multimedia modeling of the fate of triclosan and triclocarban in the Dongjiang River Basin, South China and comparison with field data

Triclosan (TCS) and triclocarban (TCC) are two active ingredients widely used in many home and personal care products. Multimedia fate of TCS and TCC in the Dongjiang River basin, South China were addressed by the developed level III fugacity model based on their usage. Under the assumption of steady state, the concentrations in air, water, soil, sediment, suspended particulate matter (SPM) and fish as well as transfer flux across the interface between the compartments were simulated. The measured concentrations for the two compounds in water, SPM, and sediment from field monitoring campaigns were then compared to validate the model. The results showed that the model predicted reasonably accurate concentrations and the differences between the measured and modeled concentrations were all less than 0.7 log units. TCS and TCC had a tendency to distribute into the sediment phase, which accounted for more than 66.3% and 90.3% of the total masses, respectively. Wastewater discharge was the main source for the occurrence of the two compounds in the aquatic environment, while degradation was the primary process for the loss in the study area, followed by the advection export. Sensitivity analysis showed that the most influential parameters for the fate of the target chemicals were source term, degradation rates and adsorption coefficients. Monte Carlo simulation could well describe the modeling uncertainty and variability.

Evaluation of the gas stripping technique for calculation of Henry's law constants using the initial slope method for 1,2,4,5-tetrachlorobenzene, pentachlorobenzene, and hexachlorobenzene

Henry's law constant (HLC) is an important factor used in environmental risk assessment and fate and transport models to describe mass transfer of chemical between water and air. HLCs and structure-property relationships were assessed for 1,2,4,5-tetrachlorobenzene (TeCB), pentachlorobenzene (PeCB), and hexachlorobenzene (HCB). HLCs were determined using the volatilization rate (kv) of sparged chemical at 25 °C. Despite the assumption that kv should be constant throughout the stripping duration, results indicated that kv decreased over time according to three separate slope regions. Results of ANCOVA indicate that kv is statistically different in the third slope region, which leads to the conclusion that use of the entire stripping data set would lead to biased HLCs. This decrease in kv may be attributed to desorption from sparger surfaces, which has not been considered widely in the literature. Statistical analysis was possible because of the robustness of the current experimental procedure which included numerous replications (15 total spargers) and extensive data points available to discern key slope changes. HLCs determined using the gas stripping technique were 57, 33, and 30 Pa m(3) mol(-1) for 1,2,4,5-TeCB, PeCB, and HCB, respectively. In comparison to literature values, current TeCB and HCB HLCs were within wide reference ranges spanning approximately an order of magnitude for each chemical. PeCB HLC of the current study was two times lower than the lowest reference data.

Simulation of the fate and seasonal variations of α-hexachlorocyclohexane in Lake Chaohu using a dynamic fugacity model

Fate and seasonal variations of α-hexachlorocyclohexane (α-HCH) were simulated using a dynamic fugacity model in Lake Chaohu, China. Sensitivity analyses were performed to identify influential parameters and Monte Carlo simulation was conducted to assess model uncertainty. The calculated and measured values of the model were in good agreement except for suspended solids, which might be due to disregarding the plankton in water. The major source of α-HCH was an input from atmospheric advection, while the major environmental outputs were atmospheric advection and sediment degradation. The net annual input and output of α-HCH were approximately 0.294 t and 0.412 t, respectively. Sediment was an important sink for α-HCH. Seasonal patterns in various media were successfully modeled and factors leading to this seasonality were discussed. Sensitivity analysis found that parameters of source and degradation were more important than the other parameters. The sediment was influenced more by various parameters than air and water were. Temperature variation had a greater impact on the dynamics of the model output than other dynamic parameters. Uncertainty analysis showed that the model uncertainty was relatively low but significantly increased in the second half of the simulation period due to the increase in the gas-water diffusion flux variability.

Guidance on the selection of efficient computational methods for multimedia fate models

Dynamic multimedia fate models (MFMs) have to deal with the temporal and spatial variation of physical-chemical properties, environmental scenarios, and chemical emissions. In such complex simulation tools, an analytical solution is not practically feasible, and even a numerical approach requires a suitable choice of the method in order to obtain satisfying speed and reliability, particularly when certain combinations of modeling scenarios and chemical properties occur. In this paper, considering some examples of a wide range of realistic chemical and scenario properties, some sources of stiffness in MFM equations are pinpointed. Next, a comparison of the performances of several numerical schemes (chosen as representatives of three wide classes) is performed. The accuracy and the computational effort required by each method is evaluated, illustrating the general effectiveness of automatically adapted timesteps in numerical algorithms and the pros and cons of implicit timestepping. The results show that automatic error control methods can significantly improve the quality of the computed solutions and most often lead to relevant savings in computing time. Additionally, explicit and implicit methods are compared, indicating that an implicit method of medium order (around 5) is the best choice as a general purpose MFM computing engine.

A field studies and modeling approach to develop organochlorine pesticide and PCB total maximum daily load calculations: case study for Echo Park Lake, Los Angeles, CA

Echo Park Lake is a small lake in Los Angeles, CA listed on the USA Clean Water Act Section 303(d) list of impaired water bodies for elevated levels of organochlorine pesticides (OCPs) and polychlorinated biphenyls (PCBs) in fish tissue. A lake water and sediment sampling program was completed to support the development of total maximum daily loads (TMDL) to address the lake impairment. The field data indicated quantifiable levels of OCPs and PCBs in the sediments, but lake water data were all below detection levels. The field sediment data obtained may explain the contaminant levels in fish tissue using appropriate sediment-water partitioning coefficients and bioaccumulation factors. A partition-equilibrium fugacity model of the whole lake system was used to interpret the field data and indicated that half of the total mass of the pollutants in the system are in the sediments and the other half is in soil; therefore, soil erosion could be a significant pollutant transport mode into the lake. Modeling also indicated that developing and quantifying the TMDL depends significantly on the analytical detection level for the pollutants in field samples and on the choice of octanol-water partitioning coefficient and bioaccumulation factors for the model.

Differential kinetics and temperature dependence of abiotic and biotic processes controlling the environmental fate of TNT in simulated marine systems

This work seeks to understand how the balance of abiotic and biotic kinetic processes in sediments control the residual concentration of TNT in marine systems after release from ocean-dumped source. Kinetics of TNT disappearance were followed using marine sediments at different temperatures and under both biotic and presumably abiotic conditions (through sodium azide addition). Sediments exhibiting the highest rate of TNT disappearance under biotic conditions also exhibited the highest sorption affinity for TNT under abiotic conditions. Significant temperature dependence in the abiotic processes was observed in the diffusion coefficient of TNT and not sediment sorption affinity. At higher temperature, kinetics of biotic processes outpaced abiotic processes, but at low temperature, kinetics of abiotic processes were much more significant. We concluded that the differential influence of temperature on the kinetics of abiotic and biotic processes could provide distinguishing predictions for the potential residual concentration of TNT contamination in marine-sediment systems.

A novel environmental indicator for monitoring of pesticides

The identification of an indicator for monitoring pesticides is a very effective analytical approach because it allows one to schedule and simplify the analytical routine. In this study, a new indicator has been designed, which has to be able to define a scale of priorities in the pesticides monitoring. A starting equation was developed involving the escaping tendency of a given substance from a phase (based on the Mackay model I level). The reliability of the indicator was tested by application to a model system, consisting of a defined and limited area, choosing water as matrix over a period of 6 years. A group of marker compounds was also defined to implement the predictive efficacy of the indicator. The results obtained by modeling were compared to those obtained by experimentation of the same model system. The indicator was subsequently and appropriately modified creating a new equation, including a kinetic factor, which considers the environmental degradation processes. The effect of the rectified indicator was consistent with the sales data list of compounds, when applied to the markers. The indicator developed in this study, tested as a model on specific area-phase-period (Province of Siena, water phase, 2000-2006), is applicable to any other area-phase-period, adjusting the partition value of the Mackay model for the case under study.

Novel control and steady-state correction method for standard 28-day bioaccumulation tests using Nereis virens

Evaluation of dredged material for aquatic placement requires assessment of bioaccumulation potentials for benthic organisms using standardized laboratory bioaccumulation tests. Critical to the interpretation of these data is the assessment of steady state for bioaccumulated residues needed to generate biota sediment accumulation factors (BSAFs) and to address control correction of day 0 contaminant residues measured in bioassay organisms. This study applied a novel performance reference compound approach with a pulse-chase experimental design to investigate elimination of a series of isotopically labeled polychlorinated biphenyl ((13)C-PCBs) in the polychaete worm Nereis virens while simultaneously evaluating native PCB bioaccumulation from field-collected sediments. Results demonstrated that all (13)C-PCBs, with the exception of (13)C-PCB209 (> 80%), were eliminated by more than 90% after 28 d. The three sediment types yielded similar (13)C-PCB whole-body elimination rate constants (k(tot)) producing the following predictive equation: log k(tot)  =  - 0.09 × log K(OW)  - 0.45. The rapid loss of (13)C-PCBs from worms over the bioassay period indicated that control correction, by subtracting day 0 residues, would result in underestimates of bioavailable sediment residues. Significant uptake of native PCBs was observed only in the most contaminated sediment and proceeded according to kinetic model predictions with steady-state BSAFs ranging from 1 to 3 and peaking for congeners of log K(OW) between 6.2 and 6.5. The performance reference compound approach can provide novel information about chemical toxicokinetics and also serve as a quality check for the physiological performance of the bioassay organism during standardized bioaccumulation testing.

TICKET-UWM: a coupled kinetic, equilibrium, and transport screening model for metals in lakes

The tableau input coupled kinetic equilibrium transport-unit world model (TICKET-UWM) has been developed as a screening model for assessing potential environmental risks associated with the release of metals into lakes. The model is based on a fully implicit, one-step solution algorithm that allows for simultaneous consideration of dissolved and particulate phase transport; metal complexation to organic matter and inorganic ligands; precipitation of metal hydroxides, carbonates, and sulfides; competitive interactions of metals and major cations with biotic ligands; a simplified description of biogeochemical cycling of organic carbon and sulfur; and dissolution kinetics for metal powders, massives, and other solid forms. Application of TICKET-UWM to a generalized lake in the Sudbury area of the Canadian Shield is presented to demonstrate the overall cycling of metals in lakes and the nonlinear effects of chemical speciation on metal responses. In addition, the model is used to calculate critical loads for metals, with acute toxicity of Daphnia magna as the final endpoint. Model results show that the critical loads for Cu, Ni, Pb, and Zn varied from 2.5 to 39.0 g metal/m(2) -year and were found to be one or more orders of magnitude higher than comparable loads for pesticides (lindane, 4,4'-DDT) and several polyaromatic hydrocarbon (PAH) compounds. In sensitivity calculations, critical metal-loading rates were found to vary significantly as a function of the hydraulic detention time, water hardness, and metal dissolution kinetic rates.

Bioaccumulation of organic contaminants in humans: a multimedia perspective and the importance of biotransformation

Bioaccumulation is an important component of the exposure hazard assessment and risk assessment of organic chemicals. Screening criteria for chemical hazard used in national and international regulations are based on the paradigm that partitioning properties are the primary chemical determinants of bioaccumulation. We use a holistic multimedia perspective to evaluate the partitioning property paradigm with respect to assessing human bioaccumulation. Multimedia bioaccumulation factors (mmBAFs) for humans were modeled for hypothetical chemicals with a wide range of physical-chemical properties. Varying partitioning properties over 12 orders of magnitude (a plausible range for nonionizing organics) resulted in only modest changes in mmBAFs (a factor of ∼ 10) for all but very volatile or hydrophilic chemicals. In contrast, varying biotransformation rate constants over 6 orders of magnitude resulted in substantial differences in mmBAFs (greater than a factor of 10(9)). Our model results are supported by empirical observations of well characterized pollutants, which demonstrate that chemicals with similar partitioning properties can have very different bioaccumulation behavior. Susceptibility to biotransformation clearly determines bioaccumulation in humans for many chemicals. We conclude that a holistic multimedia perspective for bioaccumulation assessment is necessary to develop regulations, criteria, and policies that are protective of human health and the environment.

Fate estimation of polycyclic aromatic hydrocarbons in soils in a rapid urbanization region, Shenzhen of China

Our previous study indicated that the current level of polycyclic aromatic hydrocarbons (PAHs) in Shenzhen soil is in the low-end of world soil PAH pollution. In this study, the fate of PAHs in the soil of Shenzhen was investigated. The mass inventories of Σ(27)PAHs and Σ(15)PAHs (defined as the sum of the 27 or 15 PAH compounds sought) in topsoil of Shenzhen were ∼204 and ∼152 metric tons, respectively. Fate estimation of Σ(15)PAHs shows that air-soil gaseous exchange is the primary environmental process with ∼10,076 kg/year diffusing from soil to air. Rain washing (∼1131 kg/year from air to soil) is the most important input pathway followed by wet (∼17 kg/year) and dry deposition (∼8 kg/year) to soils in Shenzhen. The transport of Σ(15)PAHs by soil erosion is a crucial loss process for soil PAHs in Shenzhen (1918 kg/year for water runoff and 657 kg/year for solid runoff from soil). Moreover, degradation is not ignorable at present (95 kg/year). Comparison of inventory and residue (defined as Σ(15)PAHs left in topsoils after all environmental loss processes) suggested that input and loss of high molecular weight PAHs for Shenzhen's soil reached apparent equilibrium. Soil PAH pollution in Shenzhen will stay in a quasi-steady state for a long period and the natural environmental processes can not significantly reduce the pollution.