How to Constrain Multi‐Objective Calibrations of the SWAT Model Using Water Balance Components

Accurate discharge simulation is one of the most common objectives of hydrological modeling studies. However, a good simulation of discharge is not necessarily the result of a realistic simulation of hydrological processes within the catchment. We propose an evaluation framework that considers both discharge and water balance components as evaluation criteria for calibration of the Soil and Water Assessment Tool (SWAT). In this study, we integrated average annual values of surface runoff, groundwater flow, and evapotranspiration in the model evaluation procedure to constrain the selection of good model runs for the Little River Experimental Watershed in Georgia, United States. For evaluating water balance and discharge dynamics, the Nash‐Sutcliffe efficiency (NSE) and percent bias (PBIAS) were used. In addition, the ratio of root mean square error and standard deviation of measured data (RSR) was calculated for individual segments of the flow duration curve to identify the best model runs in terms of discharge magnitude. Our results indicate that good statistics for discharge do not guarantee realistic simulations of individual water balance components. Therefore, we recommend constraining the ranges of water balance components to achieve a more realistic simulation of the entire hydrological system, even if tradeoffs between good statistics for discharge simulations and reasonable amounts of the water balance components are unavoidable. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Does the Temporal Resolution of Precipitation Input Influence the Simulated Hydrological Components Employing the SWAT Model?

This study aimed to evaluate the influence of sub‐daily precipitation time steps on model performance and hydrological components by applying the Green and Ampt infiltration method using the Soil and Water Assessment Tool (SWAT). Precipitation was measured at a resolution of 0.1 mm and aggregated to 5‐, 15‐, 30‐, and 60‐min time steps. Daily discharge data over a 10‐year period were used to calibrate and validate the model. Following a global sensitivity analysis, relevant parameters were optimized through an automatic calibration procedure using SWAT‐CUP for each time step. Daily performance statistics were almost equal among all four time steps (NSE ≈ 0.47). Discharge mainly consisted of groundwater flow (55%) and tile flow (42%), in reasonable proportions for the investigated catchment. In conclusion, model outputs were almost identical, showing simulations responded nearly independently of the chosen precipitation time step. This held true for (1) the selection of sensitive parameters, (2) performance statistics, (3) the shape of the hydrographs, and (4) flow components. However, a scenario analysis revealed that the precipitation time step becomes important when saturated hydraulic conductivities are low and curve numbers are high. The study suggests that there is no need in using precipitation time steps <1 h for lowland catchments dominated by soils with a low surface runoff potential if daily flow values are being considered. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

White clover clones as a cost-effective indicator of phytotoxic ozone: 10 years of experience from central Italy

Data collected at one site in central Italy using the NC-S/NC-R clover (Trifolium repens) biotype system during 1997-2007 were analysed in order to assess: (a) its performance under Mediterranean conditions; (b) variations of ozone damage linked with meteorological conditions; (c) if critical level approach is a good predictor of ozone risk on vegetation. NC-S dry biomasses were systematically lower than those of NC-R, the mean ratio being 0.7. Relevant relationship between ozone visible injury and cumulated values of AOT40 were also reported. Temperature and number of rainy days were the most important factors associated with ozone presence and, as a consequence, with leaf injury index. Photosynthetic gas exchange properties indicate that NC-S has higher values of stomatal conductance.     

A decision framework for possible remediation of contaminated sediments in the River Kymijoki, Finland

Background, aim, and scope  The paper describes the spatial contamination of the River Kymijoki, South-Eastern Finland, and the coastal region of the Gulf of Finland with PCDD/Fs and mercury. The findings of ecotoxicologial and human health studies are also reported, including environmental and human risk assessments. Sediments from the River Kymijoki, draining into the Gulf of Finland, have been heavily polluted by the pulp and paper industry and by chemical industries. A wood preservative, known as Ky-5, was manufactured in the upper reaches of the river between 1940 and 1984 causing severe pollution of river sediments with polychlorinated dibenzo-p-dioxins (PCDD) and dibenzofurans (PCDF). Moreover, the sediments have been polluted with mercury (Hg) from chlor-alkali production and the use of Hg as a slimicide in pulp and paper manufacturing. Materials and methods  An extensive sediment survey was conducted as well as sediment transport modeling, toxicity screening of sediment invertebrates, and a survey of contaminant bioaccumulation in invertebrates and fish. Studies on human exposure to PCDD/Fs and the possible effects on hypermineralization of teeth as well as an epidemiological study to reveal increased cancer risk were also conducted. An assessment of the ecological and human health risks with a null hypothesis (no remediation) was undertaken. Results  The sediment survey revealed severe contamination of river and coastal sediments with PCDD/Fs and Hg. The total volume of contaminated sediments was estimated to reach 5 × 10⁶ m³ and hot spots with extremely high concentrations (max 292,000 ng g⁻¹ or 1,060 ng I-TEQ g⁻¹ d.w.) were located immediately downstream from the pollution source (approximately 90,000 m³). Sediment contamination was accompanied by changes in benthic assemblages, but direct effects were masked by many factors. The fish showed only slightly elevated PCDD/F levels in muscle, but orders of magnitude higher in the liver compared with reference freshwater sites and the Baltic Sea. The concentrations in human fat did not reveal high human exposure in the Kymijoki area in general and was lower than in sea fishermen. The relative risk for total cancer among farmers was marginally higher (RR = 1.13) among those living close to the river, compared with farmers living further away, and the possibility of increased cancer risk cannot be ruled out. A conservative risk assessment revealed that the present probability of exceeding the WHO upper exposure limit of 4 pg WHO-TEQ kg⁻¹ d⁻¹ for PCDD/Fs and DL-PCBs was 6%. The probability of exceeding the WHO limit value of 0.23 μg kg⁻¹ d⁻¹ for methyl mercury was estimated to be notably higher at 62%. Based on these studies and the estimated risks connected with different remediation techniques a general remediation plan with cost benefit analysis was generated for several sub-regions in the river. Dredging, on-site treatment, and a close disposal of the most contaminated sediments (90,000 m³) was suggested as the first phase of the remediation. The decision regarding the start of remediation will be made during autumn 2008. Conclusions  The sediments in the River Kymijoki are heavily polluted with PCDD/Fs and mercury from earlier chlorophenol, chlor-alkali, and pulp and paper manufacturing. A continuous transport of contaminants is taking place to the Gulf of Finland in the Baltic Sea. The highly increased PCDD/F and Hg levels in river sediments pose an ecotoxicological risk to benthic fauna, to fish-eating predators and probably to human health. The risks posed by mercury exceed those from PCDD/Fs and need to be evaluated for (former) chlor-alkali sites and other mercury releasing industries as one basis for remediation decision making. Recommendations and perspectives  The studies form the basis of a risk management strategy and a plan for possible remediation of contaminated sediments currently under consideration in the Southeast Finland Regional Environment Centre. It is recommended that a detailed restoration plan for the most seriously contaminated areas should be undertaken. Based on current knowledge, the restoration of the whole river is not feasible, considering the current risk caused by the contaminated sediment in the river and the costs of an extensive restoration project. The experiences gained in the present case should be utilized in the evaluation of PCDD/F- and mercury-contaminated sites in other countries. The case demonstrates that the historic reservoirs are of contemporary relevance and should be addressed, e.g., in the national implementation plans of the Stockholm Convention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Occurrence, source diagnosis, and biological effect assessment of DDT and its metabolites in various environmental compartments of the Pearl River Delta, South China: A review

This review summarizes and analyzes available data in the literature (mostly after 2000) on the occurrence of dichlorodiphenyltrichloroethane (DDT) and its main metabolites, dichlorodiphenyldichloroethane (DDD) and chlorodiphenyldichloroethylene (DDE), in the environment of the Pearl River Delta (PRD), South China. Generally, the concentration levels of the sum of DDT, DDD, and DDE, designated as DDTs thereafter, have not significantly declined since 1983. However, the composition of DDTs residues has changed with time. DDTs in soil, freshwater sediment and freshwater fish species were mainly residues from chronological use. There is evidence that new point sources, such as dicofol and anti-fouling paint, may have contributed DDTs to various environmental compartments. Risk assessment against existing criteria indicated that the levels of DDTs in water and some fish species may pose adverse effects to humans or wildlife, and those in sediment/soil may also cause negative impacts to the eco-environment of the PRD.     

Model evaluation of the phytoextraction potential of heavy metal hyperaccumulators and non-hyperaccumulators

Evaluation of the remediation ability of zinc/cadmium in hyper- and non-hyperaccumulator plant species through greenhouse studies is limited. To bridge the gap between greenhouse studies and field applications for phytoextraction, we used published data to examine the partitioning of heavy metals between plants and soil (defined as the bioconcentration factor). We compared the remediation ability of the Zn/Cd hyperaccumulators Thlaspi caerulescens and Arabidopsis halleri and the non-hyperaccumulators Nicotiana tabacum and Brassica juncea using a hierarchical linear model (HLM). A recursive algorithm was then used to evaluate how many harvest cycles were required to clean a contaminated site to meet Taiwan Environmental Protection Agency regulations. Despite the high bioconcentration factor of both hyperaccumulators, metal removal was still limited because of the plants' small biomass. Simulation with N. tabacum and the Cadmium model suggests further study and development of plants with high biomass and improved phytoextraction potential for use in environmental cleanup.     

Long-range atmospheric transport of three toxaphene congeners across Europe. Modeling by chained single-box FATEMOD program

Background, aims, and scope  Since toxaphene (polychlorocamphene, polychloropinene, or strobane) mixtures were applied for massive insecticide use in the 1960s to replace the use of DDT, some of their congeners have been found at high latitudes far away from the usage areas. Especially polychlorinated bornanes have demonstrated dominating congeners transported by air up to the Arctic areas. Environmental fate modeling has been applied to monitor this phenomenon using parallel zones of atmosphere around the globe as interconnected environments. These zones, shown in many meteorological maps, however, may not be the best way to configure atmospheric transport in air trajectories. The latter could also be covered by connecting a chain of simple model boxes. We aim to study this alternative approach by modeling the trajectory chain using catchment boxes of our FATEMOD model. Polychlorobornanes analyzed in biota of the Barents Sea offered one case to study this modeling alternative, while toxaphene has been and partly still is used massively at southern East Europe and around rivers flowing to the Aral Sea. Materials and methods  Pure model substances of three polychlorobornanes (toxaphene congeners P26, P50, and P62) were synthesized, their environmentally important thermal properties measured by differential scanning calorimetry, as evaluated from literature data, and their temperature dependences estimated by the QSPR programs VPLEST, WATSOLU, and TDLKOW. The evaluated property parameters were used to model their atmospheric long-range transport from toxaphene heavy usage areas in Ukraine and Aral/SyrDarja/AmuDarja region areas, through East Europe and Northern Norway (Finnmarken) to the Barents Sea. The time period used for the emission model was June 1997. Usual weather conditions in June were applied in the model, which was constructed by chaining FATEMOD model boxes of the catchment’s areas along assumed maximal air flow trajectories. Analysis of the three chlorobornanes in toxaphene mixtures function as a basis for the estimates of emission levels caused by its usage. High estimate (A) was taken from contents in a Western product chlorocamphene and low estimate (B) from mean contents in Russian polychloroterpene products to achieve modeled water concentrations. Bioaccumulation to analyzed lipid of aquatic biota at the target region was estimated by using statistical calculation for persistent organic pollutants in literature. Results  The results from model runs A and B (high and low emission estimate) for levels in sea biota were compared to analysis results of samples taken in August 1997 at Barents Sea. The model results (ng g⁻¹ lw): 4–95 in lipid of planktovores and 7–150 in lipid of piscivores, were in fair agreement with the analysis results from August 1997: 21–31 in Themisto libellula (chatka), 26–42 in Boreocadus saida (Polar cod), and 5–27 in Gadus morhua (cod) liver. Discussion  The modeling results indicate that the application of chained simple multimedia catchment boxes on predicted trajectory is a useful method for estimation of volatile airborne persistent chemical exposures to biota in remote areas. For hazard assessment of these pollutants, their properties, especially temperature dependences, must be estimated by a reasonable accuracy. That can be achieved by using measurements in laboratory with pure model compounds and estimation of properties by thermodynamic QSPR methods. The property parameters can be validated by comparing their values at an environmental temperature range with measured or QSPR-estimated values derived by independent methods. The chained box method used for long-range air transport modeling can be more suitable than global parallel zones modeling used earlier, provided that the main airflow trajectories and properties of transported pollutants are predictable enough. Conclusions  Long-range air transport modeling of persistent, especially photo-resistant organic compounds using a chain of joint simple boxes of catchment’s environments is a feasible method to predict concentrations of pollutants at the target area. This is justified from model results compared with analytical measurements in Barents Sea biota in August 1997: three of six modeled values were high and the other three low compared to the analysis results. The order of magnitude level was similar in both modeled (planktovore and piscivore) and observed (chatka and polar cod) values of lipid samples. The obtained results were too limited to firm validation but are sufficient to justify feasibility of the method, which prompts one to perform more studies on this modeling system. Recommendations and perspectives  For assessment of the risk of environmental damages, chemical fate determination is an essential tool for chemical control, e.g., for EU following the REACH rules. The present conclusion of applicability of the chained single-box multimedia modeling can be validated by further studies using analyses of emissions and target biota in various other cases. To achieve useful results, fate models built with databases having automatic steps for most calculations and outputs accessible to all chemical control professionals are essential. Our FATEMOD program catchments at environments and compound properties listed in the database represent a feasible tool for local, regional, and, according our present test results, for global exposure predictions. As an extended use of model, emission estimates can be achieved by reversed modeling from analysis results of samples corresponding to the target area. This article is dedicated to the memory of Professor Alexander B Terentiev (who passed away in November 2006), our true friend. With his Institute of Organo-Element Compounds, Russian Academy of Science, Moscow, he was an important main organizer of the six joint Finnish–Russian seminars (every third year since 1989) on the field (‘Chemistry and Ecology of Organo-Element Compounds’). He prompted us especially to search properties and environmental fates for various polyhalogen compounds. We remember him for his friendly character and great sense of humor.     

The effect of alkyl chain length on the degradation of alkylimidazolium- and pyridinium-type ionic liquids in a Fenton-like system

Background, aim, and scope  Ionic liquids are regarded as essentially “green” chemicals because of their insignificant vapor pressure and, hence, are a good alternative to the emissions of toxic conventional volatile solvents. Not only because of their attractive industrial applications, but also due to their very high stability, ionic liquids could soon become persistent contaminants of technological wastewaters and, moreover, break through into natural waters following classical treatment systems. The removal of harmful organic pollutants has forced the development of new methodologies known as advanced oxidation processes (AOPs). Among them, the Fenton and Fenton-like reactions are usually modified by the use of a higher hydrogen peroxide concentration and through different catalysts. The aim of this study was to assess the effect of hydrogen peroxide concentration on degradation rates in a Fenton-like system of alkylimidazolium ionic liquids with alkyl chains of varying length and 3-methyl-N-butylpyridinium chloride. Materials and methods  The ionic liquids were oxidized in dilute aqueous solution in the presence of two different concentrations of hydrogen peroxide. All reactions were performed in the dark to prevent photoreduction of Fe(III). The concentrations of ionic liquids during the process were monitored with high-performance liquid chromatography. Preliminary degradation pathways were studied with the aid of ¹H NMR. Results  Degradation of ionic liquids in this system was quite effective. Increasing the H₂O₂ concentration from 100 to 400 mM improved ionic liquid degradation from 57–84% to 87–100% after 60 min reaction time. Resistance to degradation was weaker, the shorter the alkyl chain. Discussion  The compound omimCl was more resistant to oxidation then other compounds, which suggests that the oxidation rates of imidazolium ionic liquids by OH· are structure-dependent and are correlated with the n-alkyl chain length substituted at the N-1-position. The level of degradation was dependent on the type of head group. Replacing the imidazolium head group with pyridinium increased resistance to degradation. Nonetheless, lengthening the alkyl chain from four to eight carbons lowered the rate of ionic liquid degradation to a greater extent than changing the head group from imidazolium to pyridinium. 1H-NMR spectra show, in the first stage of degradation, that it is likely that radical attack is nonspecific, with any one of the carbon atoms in the ring and the n-alkyl chain being susceptible to attack. Conclusions  The proposed method has proven to be an efficient and reliable method for the degradation of imidazolium ionic liquids by a Fenton-like reagent deteriorated with lengthening n-alkyl substituents and by replacing the imidazolium head group with pyridinium. The enhanced resistance of 1-butyl-3-methylpyridinium chloride when the resistance of imidazolium ionic liquids decreases with increasing H₂O₂ concentration is probably indicative of a change in the degradation mechanism in a vigorous Fenton-like system. H-NMR spectra showed, in the first stage of degradation, that radical attack is nonspecific, with any one of the carbon atoms in the ring and the n-alkyl chain being susceptible to attack. Recommendations and perspectives  Since ionic liquids are now one of the most promising alternative chemicals of the future, the degradation and waste management studies should be integrated into a general development research of these chemicals. In the case of imidazolium and pyridinium ionic liquids that are known to be resistant to bio- or thermal degradation, studies in the field of AOPs should assist the future structural design as well as tailor the technological process of these chemicals     

Use of a food web model to evaluate the factors responsible for high PCB fish concentrations in Lake Ellasjøen, a high Arctic Lake

Background, aim, and scope

Lake Ellasjøen, located in the Norwegian high arctic, contains the highest concentrations of polychlorinated biphenyls (PCBs) ever recorded in fish and sediment from high arctic lakes, and concentrations are more than 10 times greater than in nearby Lake Øyangen. These elevated concentrations in Ellasjøen have been previously attributed, in part, to contaminant loadings from seabirds that use Ellasjøen, but not Øyangen, as a resting area. However, other factors, such as food web structure, organism growth rate, weight, lipid content, lake morphology, and nutrient inputs from the seabird guano, also differ between the two systems. The aim of this study is to evaluate the relative influence of these factors as explanatory variables for the higher PCB fish concentrations in Ellasjøen compared with Øyangen, using both a food web model and empirical data.

Methods

The model is based on previously developed models but parameterized for Lakes Ellasjøen and Øyangen using measured data wherever possible. The model was applied to five representative PCB congeners (PCB 105, 118, 138, 153, and 180) using measured sediment and water concentrations as input data and evaluated with previously collected food web data.

Results

Modeled concentrations are within a factor of two of measured concentrations in 60% and 40% of the cases in Lakes Ellasjøen and Øyangen, respectively, and within a factor of 10 in 100% of the cases in both lakes. In many cases, this is comparable to the variability associated with the data as well as the efficacy of the predictions of other food web model applications.

Discussion

We next used the model to quantify the relative importance of five major differences between Ellasjøen and Øyangen by replacing variables representing each of these factors in the Ellasjøen model with those from Øyangen, in separate simulations. The model predicts that the elevated PCB concentrations in Ellasjøen water and sediment account for 49%–58% of differences in modeled fish PCB concentrations between lakes. These elevated sediment and, to a lesser extent, water concentrations in Ellasjøen are due to PCB loadings from seabird guano. However, sediment–water fugacity ratios of PCBs are consistently greater in Ellasjøen compared with Øyangen, which suggests that internal lake processes also contribute to differences in sediment and water concentrations. We hypothesize that the nutrients associated with guano influence sediment–water fugacity ratios of PCBs by increasing the stock of pelagic algae. As both these algae and the guano settle, their organic carbon content is degraded faster than PCBs, which causes an extra magnification step in Ellasjøen before these detrital particles are consumed by benthic organisms, which are in turn consumed by fish. The model predicts that the remaining ～50% of the differences in PCB concentrations observed between the fish of these lakes are due to other subtle differences in their food web structures.

Conclusions

In conclusion, based on the results of a food web model, we found that the most dominant factors influencing the higher PCB fish concentrations in Lake Ellasjøen compared with Øyangen are the higher sediment and water concentrations in Ellasjøen, caused by seabird guano. Together, sediment and water are predicted to account for 49%–58% of differences in fish concentrations between lakes. Although seabird guano provides a source of nutrients to the lake, in addition to contaminants, empirical data and indirect model results suggest that nutrients are not leading to decreased bioaccumulation, in contrast to what has been observed in temperate, pelagic food webs.

Recommendations and perspectives

The results of this study emphasize the importance of considering even small differences in food web structure when comparing bioaccumulation in two lakes; although the food web structures of Ellasjøen and Øyangen differ only slightly, the model predicts that these differences account for most of the remaining ～50% of the differences in PCB fish concentrations between the two lakes. This study further demonstrates the utility of food web models as we were able to predict and tease apart the influence of various factors responsible for the elevated concentrations in the fish from Lake Ellasjøen, which would have been difficult using the field data alone.