Application of Artificial Neural Networks and Particle Swarm Optimization for the Management of Groundwater Resources

Ground management problems are typically solved by the simulation-optimization approach where complex numerical models are used to simulate the groundwater flow and/or contamination transport. These numerical models take a lot of time to solve the management problems and hence become computationally expensive. In this study, Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) models were developed and coupled for the management of groundwater of Dore river basin in France. The Analytic Element Method (AEM) based flow model was developed and used to generate the dataset for the training and testing of the ANN model. This developed ANN-PSO model was applied to minimize the pumping cost of the wells, including cost of the pipe line. The discharge and location of the pumping wells were taken as the decision variable and the ANN-PSO model was applied to find out the optimal location of the wells. The results of the ANN-PSO model are found similar to the results obtained by AEM-PSO model. The results show that the ANN model can reduce the computational burden significantly as it is able to analyze different scenarios, and the ANN-PSO model is capable of identifying the optimal location of wells efficiently.

Modelling the impact of extensive irrigation on the groundwater resources

Drastic groundwater resource depletion due to excessive extraction for irrigation is a major concern in many parts of India. In this study, an attempt was made to simulate the groundwater scenario of the catchment using ArcSWAT. Due to the restriction on the maximum initial storage, the deep aquifer component in ArcSWAT was found to be insufficient to represent the excessive groundwater depletion scenario. Hence, a separate water balance model was used for simulating the deep aquifer water table. This approach is demonstrated through a case study for the Malaprabha catchment in India. Multi-site rainfall data was used to represent the spatial variation in the catchment climatology. Model parameters were calibrated using observed monthly stream flow data. Groundwater table simulation was validated using the qualitative information available from the field. The stream flow was found to be well simulated in the model. The simulated groundwater table fluctuation is also matching reasonably well with the field observations. From the model simulations, deep aquifer water table fluctuation was found very severe in the semi-arid lower parts of the catchment, with some areas showing around 60m depletion over a period of eight years. Copyright (c) 2012 John Wiley & Sons, Ltd.

Multiobjective fuzzy optimization for sustainable groundwater management using particle swarm optimization and analytic element method

Groundwater management involves conflicting objectives as maximization of discharge contradicts the criteria of minimum pumping cost and minimum piping cost. In addition, available data contains uncertainties such as market fluctuations, variations in water levels of wells and variations of ground water policies. A fuzzy model is to be evolved to tackle the uncertainties, and a multiobjective optimization is to be conducted to simultaneously satisfy the contradicting objectives. Towards this end, a multiobjective fuzzy optimization model is evolved. To get at the upper and lower bounds of the individual objectives, particle Swarm optimization (PSO) is adopted. The analytic element method (AEM) is employed to obtain the operating potentio metric head. In this study, a multiobjective fuzzy optimization model considering three conflicting objectives is developed using PSO and AEM methods for obtaining a sustainable groundwater management policy. The developed model is applied to a case study, and it is demonstrated that the compromise solution satisfies all the objectives with adequate levels of satisfaction. Sensitivity analysis is carried out by varying the parameters, and it is shown that the effect of any such variation is quite significant. Copyright (c) 2015 John Wiley & Sons, Ltd.

Implications of using On-Farm Flood Flow Capture to recharge groundwater and mitigate flood risks along the Kings River, CA

Two large hydrologic issues face the Kings Basin, severe and chronic overdraft of about 0.16M ac-ft annually, and flood risks along the Kings River and the downstream San Joaquin River. Since 1983, these floods have caused over $1B in damage in today’s dollars. Capturing flood flows of sufficient volume could help address these two pressing issues which are relevant to many regions of the Central Valley and will only be exacerbated with climate change. However, the Kings River has high variability associated with flow magnitudes which suggests that standard engineering approaches and acquisition of sufficient acreage through purchase and easements to capture and recharge flood waters would not be cost effective. An alternative approach investigated in this study, termed On-Farm Flood Flow Capture, involved leveraging large areas of private farmland to capture flood flows for both direct and in lieu recharge. This study investigated the technical and logistical feasibility of best management practices (BMPs) associated with On-Farm Flood Flow Capture. The investigation was conducted near Helm, CA, about 20 miles west of Fresno, CA. The experimental design identified a coordinated plan to determine infiltration rates for different soil series and different crops; develop a water budget for water applied throughout the program and estimate direct and in lieu recharge; provide a preliminary assessment of potential water quality impacts; assess logistical issues associated with implementation; and provide an economic summary of the program. At check locations, we measured average infiltration rates of 4.2 in/d for all fields and noted that infiltration rates decreased asymptotically over time to about 2 – 2.5 in/d. Rates did not differ significantly between the different crops and soils tested, but were found to be about an order of magnitude higher in one field. At a 2.5 in/d infiltration rate, 100 acres are required to infiltrate 10 CFS of captured flood flows. Water quality of applied flood flows from the Kings River had concentrations of COC (constituents of concern; i.e. nitrate, electrical conductivity or EC, phosphate, ammonium, total dissolved solids or TDS) one order of magnitude or more lower than for pumped groundwater at Terranova Ranch and similarly for a broader survey of regional groundwater. Applied flood flows flushed the root zone and upper vadose zone of nitrate and salts, leading to much lower EC and nitrate concentrations to a depth of 8 feet when compared to fields in which more limited flood flows were applied or for which drip irrigation with groundwater was the sole water source. In demonstrating this technology on the farm, approximately 3,100 ac-ft was diverted, primarily from April through mid-July, with about 70% towards in lieu and 30% towards direct recharge. Substantial flood flow volumes were applied to alfalfa, wine grapes and pistachio fields. A subset of those fields, primarily wine grapes and pistachios, were used primarily to demonstrate direct recharge. For those fields about 50 – 75% of water applied was calculated going to direct recharge. Data from the check studies suggests more flood flows could have been applied and infiltrated, effectively driving up the amount of water towards direct recharge. Costs to capture flood flows for in lieu and direct recharge for this project were low compared to recharge costs for other nearby systems and in comparison to irrigating with groundwater. Moreover, the potentially high flood capture capacity of this project suggests significant flood avoidance costs savings to downstream communities along the Kings and San Joaquin Rivers. Our analyses for Terranova Ranch suggest that allocating 25% or more flood flow water towards in lieu recharge and the rest toward direct recharge will result in an economically sustainable recharge approach paid through savings from reduced groundwater pumping. Two important issues need further consideration. First, these practices are likely to leach legacy salts and nitrates from the unsaturated zone into groundwater. We develop a conceptual model of EC movement through the unsaturated zone and estimated through mass balance calculations that approximately 10 kilograms per square meter of salts will be flushed into the groundwater through displacing 12 cubic meters per square meter of unsaturated zone pore water. This flux would increase groundwater salinity but an equivalent amount of water added subsequently is predicted as needed to return to current groundwater salinity levels. All subsequent flood flow capture and recharge is expected to further decrease groundwater salinity levels. Second, the project identified important farm-scale logistical issues including irrigator training; developing cropping plans to integrate farming and recharge activities; upgrading conveyance; and quantifying results. Regional logistical issues also exist related to conveyance, integration with agricultural management, economics, required acreage and Operation and Maintenance (O&M).

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.

Borehole Optimisation System (BOS) - a case study assessing options for abstraction of urban groundwater in Nottingham, UK

The recognition that urban groundwater is a potentially valuable resource for potable and industrial uses due to growing pressures on perceived less polluted rural groundwater has led to a requirement to assess the groundwater contamination risk in urban areas from industrial contaminants such as chlorinated solvents. The development of a probabilistic risk based management tool that predicts groundwater quality at potential new urban boreholes is beneficial in determining the best sites for future resource development. The Borehole Optimisation System (BOS) is a custom Geographic Information System (GIs) application that has been developed with the objective of identifying the optimum locations for new abstraction boreholes. BOS can be applied to any aquifer subject to variable contamination risk. The system is described in more detail by Tait et al. [Tait, N.G., Davison, J.J., Whittaker, J.J., Lehame, S.A. Lerner, D.N., 2004a. Borehole Optimisation System (BOS) - a GIs based risk analysis tool for optimising the use of urban groundwater. Environmental Modelling and Software 19, 1111-1124]. This paper applies the BOS model to an urban Permo-Triassic Sandstone aquifer in the city centre of Nottingham, UK. The risk of pollution in potential new boreholes from the industrial chlorinated solvent tetrachloroethene (PCE) was assessed for this region. The risk model was validated against contaminant concentrations from 6 actual field boreholes within the study area. In these studies the model generally underestimated contaminant concentrations. A sensitivity analysis showed that the most responsive model parameters were recharge, effective porosity and contaminant degradation rate. Multiple simulations were undertaken across the study area in order to create surface maps indicating areas of low PCE concentrations, thus indicating the best locations to place new boreholes. Results indicate that northeastern, eastern and central regions have the lowest potential PCE concentrations in abstraction groundwater and therefore are the best sites for locating new boreholes. These locations coincide with aquifer areas that are confined by low permeability Mercia Mudstone deposits. Conversely southern and northwestern areas are unconfined and have shallower depth to groundwater. These areas have the highest potential PCE concentrations. These studies demonstrate the applicability of BOS as a tool for informing decision makers on the development of urban groundwater resources. (c) 2007 Elsevier Ltd. All rights reserved.

Hydraulic Subsurface Measurements and Hydrodynamic Modeling as Indicators for Groundwater Flow Systems in the Rotondo Granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.

Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire)

To assess the efficiency of different agro-environmental strategies used to reduce groundwater pollution by nitrates, transport modelling in soils and groundwater has been carried out on two withdrawal areas in an alluvial plain. In a first time, the agro-environmental model AgriFlux allowed the simulation of water and nitrates fluxes flowing to groundwater. This model was calibrated for each agro-pedological unit of the studied territory. In a second time, the application of the hydrogeological model MODFLOW-MT3D allowed the simulation of nitrate transport in groundwater for the 1980-2004 period. This soil-groundwater coupled modelling has shown that soil nature is the first factor that conditions the vulnerability to nitrates. Thus, nitrate leaching occurs preferentially under sandy soils. Efficiency of different agro-environmental operations for groundwater quality recovery was quantified. The best results are obtained by combination of (1) grassland re-installation on sandy agricultural lots located in near well protection perimeter and (2) fertilization reduction on sandy agricultural lots located in the well alimentation area upstream the near protection perimeter. On other soils, the effect of grassland on groundwater quality improvement is more limited. Nevertheless, the control of nitrate fertilisation remains essential and is justified in both near and far well protection perimeters. Modelling thus allows optimising and priorizing agro-environmental actions in alluvial agricultural zones. [Comte J.-C., Banton O., Kockmann F., Villard A., Creuzot G. (2006), Assessment of groundwater quality recovery strategies using nitrate transport modelling. Application to the Saône alluvial formations (Tournus, Saône-et-Loire), Ingénieries Eau-Agriculture-Territoires, 45, 15-28]

An Integrated Approach to Quantifying Groundwater and Surface Water Contributions of Stream Flow

The management of water resources in Ireland prior to the Water Framework Directive (WFD) has focussed on surface water and groundwater as separate entities. A critical element to the successful implementation of the
WFD is to improve our understanding of the interaction between the two and flow mechanisms by which groundwaters discharge to surface waters. An improved understanding of the contribution of groundwater to surface water is required for the classification of groundwater body status and the determination of groundwater quality thresholds. The results of the study will also have a wider application to many areas of the WFD.
A subcommittee of the WFD Groundwater Working Group (GWWG) has been formed to develop a methodology to estimate the groundwater contribution to Irish Rivers. The group has selected a number of analytical techniques to quantify components of stream flow in an Irish context (Master Recession Curve, Unit Hydrograph, Flood Studies Report methodologies and
hydrogeological analytical modelling). The components of stream flow that can be identified include deep groundwater, intermediate and overland. These analyses have been tested on seven pilot catchments that have a variety of hydrogeological settings and have been used to inform and constrain a mathematical model. The mathematical model used was the NAM (NedbØr-AfstrØmnings-Model) rainfall-runoff model which is a module of DHIs MIKE 11 modelling suite. The results from these pilot catchments have been used to develop a decision model based on catchment descriptors from GIS datasets for the selection of NAM parameters. The datasets used include the mapping of aquifers, vulnerability and subsoils, soils, the Digital Terrain Model, CORINE and lakes. The national coverage of the GIS datasets has allowed the extrapolation of the mathematical model to regional catchments across Ireland.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Coupling ground and airborne geophysical data with upscaling techniques for regional groundwater modeling of heterogeneous aquifers: Case study of a sedimentary aquifer intruded by volcanic dykes in Northern Ireland

In highly heterogeneous aquifer systems, conceptualization of regional groundwater flow models frequently results in the generalization or negligence of aquifer heterogeneities, both of which may result in erroneous model outputs. The calculation of equivalence related to hydrogeological parameters and applied to upscaling provides a means of accounting for measurement scale information but at regional scale. In this study, the Permo-Triassic Lagan Valley strategic aquifer in Northern Ireland is observed to be heterogeneous, if not discontinuous, due to subvertical trending low-permeability Tertiary dolerite dykes. Interpretation of ground and aerial magnetic surveys produces a deterministic solution to dyke locations. By measuring relative permeabilities of both the dykes and the sedimentary host rock, equivalent directional permeabilities, that determine anisotropy calculated as a function of dyke density, are obtained. This provides parameters for larger scale equivalent blocks, which can be directly imported to numerical groundwater flow models. Different conceptual models with different degrees of upscaling are numerically tested and results compared to regional flow observations. Simulation results show that the upscaled permeabilities from geophysical data allow one to properly account for the observed spatial variations of groundwater flow, without requiring artificial distribution of aquifer properties. It is also found that an intermediate degree of upscaling, between accounting for mapped field-scale dykes and accounting for one regional anisotropy value (maximum upscaling) provides results the closest to the observations at the regional scale.

Integrating Multiple Point Statistics with Aerial Geophysical Data to assist Groundwater Flow Models

The process of accounting for heterogeneity has made significant advances in statistical research, primarily in the framework of stochastic analysis and the development of multiple-point statistics (MPS). Among MPS techniques, the direct sampling (DS) method is tested to determine its ability to delineate heterogeneity from aerial magnetics data in a regional sandstone aquifer intruded by low-permeability volcanic dykes in Northern Ireland, UK. The use of two two-dimensional bivariate training images aids in creating spatial probability distributions of heterogeneities of hydrogeological interest, despite relatively ‘noisy’ magnetics data (i.e. including hydrogeologically irrelevant urban noise and regional geologic effects). These distributions are incorporated into a hierarchy system where previously published density function and upscaling methods are applied to derive regional distributions of equivalent hydraulic conductivity tensor K. Several K models, as determined by several stochastic realisations of MPS dyke locations, are computed within groundwater flow models and evaluated by comparing modelled heads with field observations. Results show a significant improvement in model calibration when compared to a simplistic homogeneous and isotropic aquifer model that does not account for the dyke occurrence evidenced by airborne magnetic data. The best model is obtained when normal and reverse polarity dykes are computed separately within MPS simulations and when a probability threshold of 0.7 is applied. The presented stochastic approach also provides improvement when compared to a previously published deterministic anisotropic model based on the unprocessed (i.e. noisy) airborne magnetics. This demonstrates the potential of coupling MPS to airborne geophysical data for regional groundwater modelling.

A field and reactive transport model study of arsenic in a basaltic rock aquifer

The use of geothermal energy as a source for electricity and district heating has increased over recent decades. Dissolved As can be an important constituent of the geothermal fluids brought to the Earth's surface. Here the field application of laboratory measured adsorption coefficients of aqueous As species on basaltic glass surfaces is discussed. The mobility of As species in the basaltic aquifer in the Nesjavellir geothermal system, Iceland was modelled by the one-dimensional (1D) reactive transport model PHREEQC ver. 2, constrained by a long time series of field measurements with the chemical composition of geothermal effluent fluids, pH, Eh and, occasionally, Fe- and As-dissolved species measurements. Di-, tri- and tetrathioarsenic species (As(OH)S22-, AsS₃H^2-, AsS33- and As(SH)4-) were the dominant form of dissolved As in geothermal waters exiting the power plant (2.556μM total As) but converted to some extent to arsenite (H₃AsO₃) and arsenate HAsO42- oxyanions coinciding with rapid oxidation of S2- to S2O32- and finally to SO42- during surface runoff before feeding into a basaltic lava field with a total As concentration of 0.882μM following dilution with other surface waters. A continuous 25-a data set monitoring groundwater chemistry along a cross section of warm springs on the Lake Thingvallavatn shoreline allowed calibration of the 1D model. Furthermore, a series of ground water wells located in the basaltic lava field, provided access along the line of flow of the geothermal effluent waters towards the lake. The conservative ion Cl^- moved through the basaltic lava field (4100m) in less than10a but As was retarded considerably due to surface reactions and has entered a groundwater well 850m down the flow path as arsenate in accordance to the prediction of the 1D model. The 1D model predicted a complete breakthrough of arsenate in the year 2100. In a reduced system arsenite should be retained for about 1ka. © 2011 Elsevier Ltd.

Project, modelling and simulation of air stripping columns for the treatment of groundwater contaminated with volatile organic compounds

Volatile organic compounds are a common source of groundwater contamination that can be easily removed by air stripping in columns with random packing and using a counter-current flow between the phases. This work proposes a new methodology for the column design for any particular type of packing and contaminant avoiding the necessity of a pre-defined diameter used in the classical approach. It also renders unnecessary the employment of the graphical Eckert generalized correlation for pressure drop estimates. The hydraulic features are previously chosen as a project criterion and only afterwards the mass transfer phenomena are incorporated, in opposition to conventional approach. The design procedure was translated into a convenient algorithm using C++ as programming language. A column was built in order to test the models used either in the design or in the simulation of the column performance. The experiments were fulfilled using a solution of chloroform in distilled water. Another model was built to simulate the operational performance of the column, both in steady state and in transient conditions. It consists in a system of two partial non linear differential equations (distributed parameters). Nevertheless, when flows are steady, the system became linear, although there is not an evident solution in analytical terms. In steady state the resulting system of ODE can be solved, allowing for the calculation of the concentration profile in both phases inside the column. In transient state the system of PDE was numerically solved by finite differences, after a previous linearization.

Spectral solution for the air stripping pollutants removal dynamic model with non linear steady state conditions

This work deals with the numerical simulation of air stripping process for the pre-treatment of groundwater used in human consumption. The model established in steady state presents an exponential solution that is used, together with the Tau Method, to get a spectral approach of the solution of the system of partial differential equations associated to the model in transient state.