Characterising the hydrological regime of an ungauged temporary river system: a case study

Temporary streams are characterised by specific hydrological regimes, which influence ecosystem processes, groundwater and surface water interactions, sediment regime, nutrient delivery, water quality and ecological status. This paper presents a methodology to characterise and classify the regime of a temporary river in Southern Italy based on hydrological indicators (HIs) computed with long-term daily flow records. By using a principal component analysis (PCA), a set of non-redundant indices were identified describing the main characteristics of the hydrological regime in the study area. The indicators identified were the annual maximum 30- and 90-day mean (DH4 and DH5), the number of zero flow days (DL6), flow permanence (MF) and the 6-month seasonal predictability of dry periods (SD6). A methodology was also tested to estimate selected HIs in ungauged river reaches. Watershed characteristics such as catchment area, gauging station elevation, mean watershed slope, mean annual rainfall, land use, soil hydraulic conductivity and available water content were derived for each site. Selected indicators were then linked to the catchment characteristics using a regression analysis. Finally, MF and SD6 were used to classify the river reaches on the basis of their degree of intermittency. The methodology presented in this paper constitutes a useful tool for ecologists and water resource managers in the Water Framework Directive implementation process, which requires a characterisation of the hydrological regime and a ‘river type’ classification for all water bodies.     

Integrated methodology for assessing the HCH groundwater pollution at the multi-source contaminated mega-site Bitterfeld/Wolfen

A large-scale groundwater contamination characterises the Pleistocene groundwater system of the former industrial and abandoned mining region Bitterfeld/Wolfen, Eastern Germany. For more than a century, local chemical production and extensive lignite mining caused a complex contaminant release from local production areas and related dump sites. Today, organic pollutants (mainly organochlorines) are present in all compartments of the environment at high concentration levels. An integrated methodology for characterising the current situation of pollution as well as the future fate development of hazardous substances is highly required to decide on further management and remediation strategies. Data analyses have been performed on regional groundwater monitoring data from about 10 years, containing approximately 3,500 samples, and up to 180 individual organic parameters from almost 250 observation wells. Run-off measurements as well as water samples were taken biweekly from local creeks during a period of 18 months. A kriging interpolation procedure was applied on groundwater analytics to generate continuous distribution patterns of the nodal contaminant samples. High-resolution geological 3-D modelling serves as a database for a regional 3-D groundwater flow model. Simulation results support the future fate assessment of contaminants. A first conceptual model of the contamination has been developed to characterise the contamination in regional surface waters and groundwater. A reliable explanation of the variant hexachlorocyclohexane (HCH) occurrence within the two local aquifer systems has been derived from the regionalised distribution patterns. Simulation results from groundwater flow modelling provide a better understanding of the future pollutant migration paths and support the overall site characterisation. The presented case study indicates that an integrated assessment of large-scale groundwater contaminations often needs more data than only from local groundwater monitoring. The developed methodology is appropriate to assess POP-contaminated mega-sites including, e.g. HCH deposits. Although HCH isomers are relevant groundwater pollutants at this site, further organochlorine pollutants are present at considerably higher levels. The study demonstrates that an effective evaluation of the current situation of contamination as well as of the related future fate development requires detailed information of the entire observed system.     

Groundwater vulnerability and risk mapping using GIS, modeling and a fuzzy logic tool

A groundwater vulnerability and risk mapping assessment, based on a source-pathway-receptor approach, is presented for an urban coastal aquifer in northeastern Brazil. A modified version of the DRASTIC methodology was used to map the intrinsic and specific groundwater vulnerability of a 292 km(2) study area. A fuzzy hierarchy methodology was adopted to evaluate the potential contaminant source index, including diffuse and point sources. Numerical modeling was performed for delineation of well capture zones, using MODFLOW and MODPATH. The integration of these elements provided the mechanism to assess groundwater pollution risks and identify areas that must be prioritized in terms of groundwater monitoring and restriction on use. A groundwater quality index based on nitrate and chloride concentrations was calculated, which had a positive correlation with the specific vulnerability index.     

Application of Unknown Groundwater Pollution Source Release History Estimation Methodology to Distributed Sources Incorporating Surface-Groundwater Interactions

Sources of contamination of groundwater are often difficult to characterize. However, it is essential for effective remediation of polluted groundwater resources. This study demonstrates an application of the linked simulation-optimization based methodology to estimate the release history from spatially distributed sources of pollution at an illustrative abandoned mine-site. In linked simulation-optimization approaches a numerical groundwater flow and transport simulation model is linked to the optimization model. In this study, topographic and geologic characteristics of the abandoned mine-site were simulated using a three-dimensional (3D) numerical groundwater flow model. Transport of contaminant in the groundwater was simulated using a 3D transient advective-dispersive contaminant transport model. Adsorption or chemical reaction of the contaminant was not considered in the contaminant transport model. Adaptive simulated annealing (ASA) was employed for solving the optimization problem. An optimization algorithm generates the candidate solutions corresponding to various unknown groundwater source characteristics. The candidate solutions are used as input in the numerical groundwater transport simulation model to generate the concentration of pollutant in the study area. This information is used to calculate the objective function value, which is utilized by the optimization algorithm to improve the candidate solution. This process continues until an optimal solution is obtained. Optimal solutions obtained in this study show that the linked simulation-optimization based methodology is potentially applicable for the characterization of spatially distributed pollutant sources, typically present at abandoned mine-sites.     

Development and application of a methodology for determining background groundwater quality at the Savannah River Site

A statistical methodology formulated for defining background or baseline levels of constituents of concern in groundwater is presented. The methodology was developed for the case where prior delineation of unimpacted areas is not possible because of site history and a large set of groundwater monitoring measurements exists. Consideration was given to spatial and temporal trends, outliers, and final segregation of wells into impacted or unimpacted categories to develop probability distributions and summary statistics for each constituent evaluated. The formulated approaches were applied to groundwater monitoring data for the U.S. Department of Energy Savannah River Site facility, and results for four representative constituents (aluminum, arsenic, mercury, and tritium) are discussed.     

Analysis of water flow paths: methodology and example calculations for a potential geological repository in Sweden

Safety assessment related to the siting of a geological repository for spent nuclear fuel deep in the bedrock requires identification of potential flow paths and the associated travel times for radionuclides originating at repository depth. Using the Laxemar candidate site in Sweden as a case study, this paper describes modeling methodology, data integration, and the resulting water flow models, focusing on the Quaternary deposits and the upper 150 m of the bedrock. Example simulations identify flow paths to groundwater discharge areas and flow paths in the surface system. The majority of the simulated groundwater flow paths end up in the main surface waters and along the coastline, even though the particles used to trace the flow paths are introduced with a uniform spatial distribution at a relatively shallow depth. The calculated groundwater travel time, determining the time available for decay and retention of radionuclides, is on average longer to the coastal bays than to other biosphere objects at the site. Further, it is demonstrated how GIS-based modeling can be used to limit the number of surface flow paths that need to be characterized for safety assessment. Based on the results, the paper discusses an approach for coupling the present models to a model for groundwater flow in the deep bedrock.     

Modelling assessment of regional groundwater contamination due to historic smelter emissions of heavy metals

Historic emissions from ore smelters typically cause regional soil contamination. We developed a modelling approach to assess the impact of such contamination on groundwater and surface water load, coupling unsaturated zone leaching modelling with 3D groundwater transport modelling. Both historic and predictive modelling were performed, using a mass balance approach for three different catchments in the vicinity of three smelters. The catchments differ in their hydrology and geochemistry. The historic modelling results indicate that leaching to groundwater is spatially very heterogeneous due to variation in soil characteristics, in particular soil pH. In the saturated zone, cadmium is becoming strongly retarded due to strong sorption at neutral pH, even though the reactivity of the sandy sediments is low. A comparison between two datasets (from 1990 to 2002) on shallow groundwater and modelled concentrations provided a useful verification on the level of statistics of "homogeneous areas" (areas with comparable land use, soil type and geohydrological situation) instead of comparison at individual locations. While at individual locations observations and the model varies up to two orders of magnitude, for homogeneous areas, medians and ranges of measured concentrations and the model results are similar. A sensitivity analysis on metal input loads, groundwater composition and sediment geochemistry reveals that the best available information scenario based on the median value of input parameters for the model predicts the range in observed concentrations very well. However, the model results are sensitive to the sediment contents of the reactive components (organic matter, clay minerals and iron oxides). Uncertainty in metal input loads and groundwater chemistry are of lesser importance. Predictive modelling reveals a remarkable difference in geochemical and hydrological controls on subsurface metal transport at catchment-scale. Whether the surface water load will peak within a few decades or continue to increase until after 2050 depends on the dominant land use functions in the areas, their hydrology and geochemical build-up.     

Biomonitoring for metal contamination near two Superfund sites in Woburn, Massachusetts, using phytochelatins

Characterizing the spatial extent of groundwater metal contamination traditionally requires installing sampling wells, an expensive and time-consuming process in urban areas. Moreover, extrapolating biotic effects from metal concentrations alone is problematic, making ecological risk assessment difficult. Our study is the first to examine the use of phytochelatin measurements in tree leaves for delimiting biological metal stress in shallow, metal-contaminated groundwater systems. Three tree species (Rhamnus frangula, Acer platanoides, and Betula populifolia) growing above the shallow groundwater aquifer of the Aberjona River watershed in Woburn, Massachusetts, display a pattern of phytochelatin production consistent with known sources of metal contamination and groundwater flow direction near the Industri-Plex Superfund site. Results also suggest the existence of a second area of contaminated groundwater and elevated metal stress near the Wells G&H Superfund site downstream, in agreement with a recent EPA ecological risk assessment. Possible contamination pathways at this site are discussed.     

The ChimERA project: coupling mechanistic exposure and effect models into an integrated platform for ecological risk assessment

Current techniques for the ecological risk assessment of chemical substances are often criticised for their lack of environmental realism, ecological relevance and methodological accuracy. ChimERA is a 3-year project (2013–2016), funded by Cefic’s Long Range Initiative (LRI) that aims to address some of these concerns by developing and testing mechanistic fate and effect models, and coupling of these models into one integrated platform for risk assessment. This paper discusses the backdrop against which this project was initiated and lists its objectives and planned methodology.     

SALSA: a simulation tool to assess ecological sustainability of agricultural production

In order to assess the ecological sustainability of agricultural production systems, there is a need for effective tools. We describe an environmental systems analysis tool called SALSA (Systems Ana/ysis for Sustainable Agriculture). It consists of substance/material flow models in which the simulation results are interpreted with life-cycle assessment methodology. The application of SALSA is demonstrated in a case study in which three different ways of producing pigs are compared with respect to energy input and the environmental impacts of global warming, eutrophication, and acidification. The scenario that combined a low-protein diet without soy meal with an improved manure-management technique with low nitrogen losses was the best for all impact categories studied. The strength of the SALSA models was their capacity to capture consequences of management options that had an influence on several processes on a farm, which enabled the type of complex studies we describe.     

Influence of soil hydrological pathways on stream aluminium chemistry at Llyn Brianne, mid-wales

Storm runoff in afforested catchments at Llyn Brianne is acidic and Al-bearing. At baseflows, stream water is well-buffered with low Al levels. This paper presents the results of a study into how hydrological pathways account for these variations in stream-water chemistry. The investigation was carried out in the LI1 catchment; a 0.4-ha subcatchment covered by stagnohumic gley soils was monitored between October 1988 and September 1989. An instrumented hill-slope was established to identify the hydrological pathways that control the hydrochemistry of storm runoff draining from the subcatchment. Perched watertables developed in the surface horizons of the soil during storm episodes and produced lateral flow above the impeding subsoil. This near-surface flow path was responsible for generating acid, Al-rich storm runoff. Some water drained vertically through the soil profile into the underlying slope drift; seepage from groundwater in the drift sustained baseflows. Buffering reactions in the groundwater zone reduced the acidity and Al levels of baseflows. These hydrochemical characteristics are likely to be representative of other areas of stagnohumic gley soils, which cover 19% of the LI1 catchment: these soils may therefore provide a substantial source of acid, Al-bearing storm runoff in LI1 and similar afforested catchments.     

Breakthrough Curves Characterization and Identification of an Unknown Pollution Source in Groundwater System Using an Artificial Neural Network (ANN)

Contamination of groundwater constrains its uses and poses a serious threat to the environment. Once groundwater is contaminated, the cleanup may be difficult and expensive. Identification of unknown pollution sources is the first step toward adopting any remediation strategy. The proposed methodology exploits the capability of a universal function approximation by a feed-forward multilayer artificial neural network (ANN) to identify the sources in terms of its location, magnitudes, and duration of activity. The back-propagation algorithm is utilized for training the ANN to identify the source characteristics based on simulated concentration data at specified observation locations in the aquifer. Uniform random generation and the Latin hypercube sampling method of random generation are used to generate temporal varying source fluxes. These source fluxes are used in groundwater flow and the transport simulation model to generate necessary data for the ANN model-building processes. Breakthrough curves obtained for the specified pollution scenario are characterized by different methods. The characterized breakthrough curves parameters serve as inputs to ANN model. Unknown pollution source characteristics are outputs for ANN model. Experimentation is also performed with different number of training and testing patterns. In addition, the effects of measurement errors in concentration measurements values are used to show the robustness of ANN based methodology for source identification in case of erroneous data.     

From Site Data to Safety Assessment: Analysis of Present and Future Hydrological Conditions at a Coastal Site in Sweden

This paper presents an analysis of present and future hydrological conditions at the Forsmark site in Sweden, which has been proposed as the site for a geological repository for spent nuclear fuel. Forsmark is a coastal site that changes in response to shoreline displacement. In the considered time frame (until year 10 000 ad), the hydrological system will be affected by landscape succession associated with shoreline displacement and changes in vegetation, regolith stratigraphy, and climate. Based on extensive site investigations and modeling of present hydrological conditions, the effects of different processes on future site hydrology are quantified. As expected, shoreline displacement has a strong effect on local hydrology (e.g., groundwater flow) in areas that change from sea to land. The comparison between present and future land areas emphasizes the importance of climate variables relative to other factors for main hydrological features such as water balances.     

Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

Mobility of heavy metals from tailings to stream waters in a mining activity contaminated site

In this paper the results of a recent characterization of Rio Piscinas (SW of Sardinia, Italy) hydrological basin are reported. In such area (about 50 km2), previous mining activities caused a serious heavy metal contamination of surface waters, groundwater, soils and biota. Acid mine drainage phenomena were observed in the area. The main sources of contamination are the tailings stored in mine tunnels and abandoned along fluvial banks. A methodological approach was adopted in order to identify relations between tailings and water contamination. Representative samples of tailings and stream sediments samples were collected. XRD analyses were performed for mineralogical characterization, while acid digestion was carried out for determining metal contents. Batch sequential leaching tests were performed in order to assess metal mobility. Also groundwater and stream water were sampled in specific locations and suitably characterized. All information collected allowed the understanding of the effect of tailings on water contamination, thus contributing to the qualitative prediction of pollution evolution on the basis of metal mobility. Finally, a potential remediation strategy of stream water is proposed.     

Hydrological change--climate change impact simulations for Sweden

Climate change resulting from the enhanced greenhouse effect is expected to give rise to changes in hydrological systems. This hydrological change, as with the change in climate variables, will vary regionally around the globe. Impact studies at local and regional scales are needed to assess how different regions will be affected. This study focuses on assessment of hydrological impacts of climate change over a wide range of Swedish basins. Different methods of transferring the signal of climate change from climate models to hydrological models were used. Several hydrological model simulations using regional climate model scenarios from Swedish Regional Climate Modelling Programme (SWECLIM) are presented. A principal conclusion is that subregional impacts to river flow vary considerably according to whether a basin is in northern or southern Sweden. Furthermore, projected hydrological change is just as dependent on the choice of the global climate model used for regional climate model boundary conditions as the choice of anthropogenic emissions scenario.     

Dynamic factor analysis of groundwater quality trends in an agricultural area adjacent to Everglades National Park

The extensive eastern boundary of Everglades National Park (ENP) in south Florida (USA) is subject to one of the most expensive and ambitious environmental restoration projects in history. Understanding and predicting the water quality interactions between the shallow aquifer and surface water is a key component in meeting current environmental regulations and fine-tuning ENP wetland restoration while still maintaining flood protection for the adjacent developed areas. Dynamic factor analysis (DFA), a recent technique for the study of multivariate non-stationary time-series, was applied to study fluctuations in groundwater quality in the area. More than two years of hydrological and water quality time series (rainfall; water table depth; and soil, ground and surface water concentrations of N-NO3-, N-NH4+, P-PO4(3-), Total P, F-and Cl-) from a small agricultural watershed adjacent to the ENP were selected for the study. The unexplained variability required for determining the concentration of each chemical in the 16 wells was greatly reduced by including in the analysis some of the observed time series as explanatory variables (rainfall, water table depth, and soil and canal water chemical concentration). DFA results showed that groundwater concentration of three of the agrochemical species studied (N-NO3-, P-PO4(3-)and Total P) were affected by the same explanatory variables (water table depth, enriched topsoil, and occurrence of a leaching rainfall event, in order of decreasing relative importance). This indicates that leaching by rainfall is the main mechanism explaining concentration peaks in groundwater. In the case of N-NH4+, in addition to leaching, groundwater concentration is governed by lateral exchange with canals. F-and Cl- are mainly affected by periods of dilution by rainfall recharge, and by exchange with the canals. The unstructured nature of the common trends found suggests that these are related to the complex spatially and temporally varying land use patterns in the watershed. The results indicate that peak concentrations of agrochemicals in groundwater could be reduced by improving fertilization practices (by splitting and modifying timing of applications) and by operating the regional canal system to maintain the water table low, especially during the rainy periods.     

A screening method for preliminary assessment of risk to groundwater from land-applied chemicals

A simple mathematical model for initial screening is presented that can aid in evaluating the relative risk to groundwater from applying nonpolar synthetic organic chemicals to soil. The basic premise is that the magnitude of the quotient of the chemical concentration of the water entering the aquifer and the maximum allowable concentration (as established by EPA or Health Departments) represents the health risk of a chemical. The chemical concentration of the soil water is estimated based on conservative, simplifying assumptions and requires only readily available data such as: basic soil properties (organic matter and saturated hydraulic conductivity), organic chemical properties (octanol-water partition coefficient and degradation rate) and environmental factors (recharge rate and depth to groundwater).The methodology was applied to assess the relative risk of organic chemicals in municipal sewage sludge and pesticides applied to agricultural land. The results are realistic.     

The comparison of rapid bioassays for the assessment of urban groundwater quality

Groundwater is a complex mixture of chemicals that is naturally variable. Current legislation in the UK requires that groundwater quality and the degree of contamination are assessed using chemical methods. Such methods do not consider the synergistic or antagonistic interactions that may affect the bioavailability and toxicity of pollutants in the environment. Bioassays are a method for assessing the toxic impact of whole groundwater samples on the environment. Three rapid bioassays, Eclox, Microtox and ToxAlert, and a Daphnia magna 48-h immobilisation test were used to assess groundwater quality from sites with a wide range of historical uses. Eclox responses indicated that the test was very sensitive to changes in groundwater chemistry; 77% of the results had a percentage inhibition greater than 90%. ToxAlert, although suitable for monitoring changes in water quality under laboratory conditions, produced highly variable results due to fluctuations in temperature and the chemical composition of the samples. Microtox produced replicable results that correlated with those from D. magna tests.     

深圳市河流沉积物重金属污染特征及评价

分析了深圳市的深圳河、布吉河、龙岗河和茅洲河等4条河流沉积物中重金属的含量和富集状况,并对重金属污染的潜在生态危害进行了评价。结果表明：（1）1996-2003年,深圳河、布吉河、龙岗河和茅洲河沉积物中重金属的含量都有上升的趋势;布吉河和龙岗河沉积物中的重金属含量较高,特别是Cu、Zn和Cr的含量比其他两条河流高出较多,深圳河居中,茅洲河较小;这4条河流都以Cd、Cu、Cr和Zn的富集为主。（2）在研究时段内,这4条河流沉积物重金属污染的潜在生态危害指数（RI）都有上升的趋势;其中深圳河,布吉河和龙岗河的RI在多数年份都达到强或很强的程度,Cd、Hg、Cu为主要贡献元素;4条河流中,茅洲河受重金属污染最小,其RI为轻微或中等,只有Cd和Hg的生态危害系数为中等或强。针对深圳市河流沉积物中重金属的累积机理,建议通过改善区域水文条件、控制污染物的排放等措施加以控制和治理。