Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal     

A general solution for tide‐induced groundwater fluctuation in an estuarine‐coastal confined and unconfined aquifer system

This paper presents an analytical solution to tide‐induced head fluctuations in a two‐dimensional estuarine‐coastal aquifer system that consists of an unconfined aquifer and a heterogeneous confined aquifer extending under a tidal river with a semipermeable layer between them. This study considers the joint effects of tidal‐river leakage, inland leakage, dimensionless transmissivity between the tidal‐river and inland confined aquifer, and transmissivity anisotropic ratios. The analytical solution for this model is obtained via the separation of variables method. Three existing solutions related to head fluctuation in one‐ or two‐dimensional leaky confined aquifers are considered as special cases in the present solution. This study shows that there is a threshold of tidal‐river confined aquifer length. When the tidal‐river length is greater than the threshold length, the inland head fluctuations remain sensitive to the leakage effect but become insensitive to the tidal‐river width and dimensionless transmissivity. Considering leakage and transmissivity anisotropy, this study also demonstrates that at a location farther from the river–inland boundary, head fluctuations increase with increasing leakage and transmissivity anisotropy; the maximum head fluctuation occurs when leakage and transmissivity anisotropy are both at their maximum values. The combined action of the 3 effects of loading, tidal‐river aquifer leakage, and inland aquifer leakage differs significantly according to various aquifer parameters. The analytical solution in this paper can be applied to demonstrate the behaviours of the head fluctuations of an estuarine‐coastal aquifer system, and the head fluctuations can be clearly described when the tidal and hydrogeological parameters are derived from field measurement data or hypothetical cases.     

Study of the feasibility of an aquifer storage and recovery system in a deep aquifer in Belgium

Abstract

The feasibility of aquifer storage and recovery (ASR) was tested in a deep aquifer near Koksijde, Belgium. To achieve this, oxic drinking water was injected into a deep aquifer (the Tienen Formation) that contains anoxic brackish water. The hydraulic properties of the aquifer were determined using a step-drawdown test. Chemical processes caused by the injection of the water were studied by two push—pull tests. The step-drawdown test was interpreted by means of an inverse numerical model, resulting in a transmissivity of 3.38 m²/d and a well loss coefficient of 0.00038 d²/m⁵. The push—pull tests identified mixing between the injection and pristine waters, and cation exchange, as the major processes determining the quality of the recovered water. Mobilization of DOC, aerobic respiration, denitrification and mobilization of phosphate were also observed.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Determination of the conduit structure in a karst aquifer based on tracer data—Lurbach system,Austria

A structure model was used to analyse solute‐transport parameter estimates based on tracer breakthrough curves. In the model system, groundwater flow is envisioned to be organised in a complex conduit network providing a variety of short circuits with relative small carrying capacities along different erosion levels. The discharge through the fully filled conduits is limited owing to void geometries and turbulent flow; thus, a hierarchic overflow system evolves where conduits are (re‐)activated or dried up depending on the flow condition. Exemplified on the Lurbach–Tanneben karst aquifer, the applicability of the model approach was tested. Information derived from multi‐tracer experiments performed at different volumetric flow rates enabled to develop a structural model of the karst network, under constraint of the geomorphological and hydrological evolution of the site. Depending on the flow rate, groundwater is divided into up to eight flow paths. The spatial hierarchy of flow paths controls the sequence of flow path activation. Conduits of the topmost level are strongly influenced by reversible alteration processes. Sedimentation or blocking causes an overflow of water to the next higher conduit. Flow path specific dissolutional denudation rates were estimated using the temporal development of the partial discharge. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of benzene transport in a two‐dimensional aquifer model

In this study, the fate and transport of aqueous benzene was investigated in a laboratory‐scale homogeneous aquifer by conducting a two‐dimensional plume test. Benzene solution was introduced as a pulse type along the width of the aquifer model through a recharge zone situated at the upper‐left part of the model and followed by a steady state flow. Solution samples were collected at various locations on the front side of the model to capture two‐dimensional plumes at discrete time intervals. The benzene plumes showed a moderate retardation relative to chloride plumes observed from the previous study conducted for the same aquifer model. The retardation factor was obtained from the ratio of travel distances of benzene peaks to chloride peaks from the injection point, computed using a line integral method. Mass recovery of aqueous benzene revealed that there was a significant reduction of benzene mass, indicating the occurrence of volatilization and/or irreversible sorption during transport. Thus, retardation along with volatilization and/or irreversible sorption may be important processes affecting the fate and transport of aqueous benzene in the aquifer model. Copyright © 2005 John Wiley & Sons, Ltd.     

Spatio‐temporal analysis of groundwater recharge and mound dynamics in an unconfined aquifer: a GIS‐based approach

Groundwater recharge and mounding of water‐table is a complex phenomenon involving time‐ and space‐dependent hydrologic processes. The effect of long‐term groundwater mounding in the aquifer depends on soil, aquifer geometry and the area contributing to recharge. In this paper, a GIS‐based spatio‐temporal algorithm has been developed for the groundwater mound dynamics to estimate the potential rise in the water‐table and groundwater volume balance residual in an unconfined aquifer. The recharge and mound dynamics as predicted using the methodology recommended here were compared with those using the Hantush equation, and the differences were quite significant. The significance of the study is to assess the effectiveness of the basin in terms of its hydrologic and hydraulic properties for sustainable management of groundwater recharge. Copyright © 2007 John Wiley & Sons, Ltd.     

The effect of low‐permeability fault zones on groundwater flow in a compartmentalized system. Experimental evidence from a carbonate aquifer (Southern Italy)

The hydrogeological behaviour of fault zones in carbonate aquifers is often neglected in conceptual and numerical models. Furthermore, no information is available regarding the relationships between piezometric levels when significant compartmentalization occurs due to the occurrence of low‐flow fault zones. The aim of this study was to refine the conceptualization of subsurface flow in faulted carbonate aquifers and to analyse relationships between sub‐basins within a compartmentalized aquifer system in Southern Italy. The interactions between compartments that straddle low‐flow faults were investigated over four hydrologic years using a statistical approach to compare (i) the hydraulic heads within two wells located up‐ and down‐gradient of tectonic discontinuities as well as (ii) the rainfall and piezometric levels. The results of this study suggest that a set of barriers exists between the wells, and, therefore, the total head loss observed between the wells (approximately 80 m) should be distributed across several aquitards, with one aquitard exhibiting a relatively high permeability or low degree of integrity. Due to slight differences in permeability, transient conditions in aquitards can occur over relatively short periods, which is in agreement with the results of the statistical data analysis. Consequently, rather than being caused by pure aquitards, aquifer system compartmentalization likely results from slight differences in the permeability between lower‐permeability fault zones and adjacent higher‐permeability protoliths. Copyright © 2014 John Wiley & Sons, Ltd.     

Subspace system identification of support‐excited structures—part I: theory and black‐box system identification

This paper reviews the theoretical principles of subspace system identification as applied to the problem of estimating black‐box state‐space models of support‐excited structures (e.g., structures exposed to earthquakes). The work distinguishes itself from past studies by providing readers with a powerful geometric interpretation of subspace operations that relates directly to theoretical structural dynamics. To validate the performance of subspace system identification, a series of experiments are conducted on a multistory steel frame structure exposed to moderate seismic ground motions; structural response data is used off‐line to estimate black‐box state‐space models. Ground motions and structural response measurements are used by the subspace system identification method to derive a complete input–output state‐space model of the steel frame system. The modal parameters of the structure are extracted from the estimated input–output state‐space model. With the use of only structural response data, output‐only state‐space models of the system are also estimated by subspace system identification. The paper concludes with a comparison study of the modal parameters extracted from the input–output and output‐only state‐space models in order to quantify the uncertainties present in modal parameters extracted from output‐only models. Copyright © 2012 John Wiley & Sons, Ltd.     

Multireservoir system optimization in the Han River basin using multi‐objective genetic algorithms

In this study, NSGA‐II is applied to multireservoir system optimization. Here, a four‐dimensional multireservoir system in the Han River basin was formulated. Two objective functions and three cases having different constraint conditions are used to achieve nondominated solutions. NSGA‐II effectively determines these solutions without being subject to any user‐defined penalty function, as it is applied to a multireservoir system optimization having a number of constraints (here, 246), multi‐objectives, and infeasible initial solutions. Most research by multi‐objective genetic algorithms only reveals a trade‐off in the objective function space present, and thus the decision maker must reanalyse this trade‐off relationship in order to obtain information on the decision variable. Contrastingly, this study suggests a method for identifying the best solutions among the nondominated ones by analysing the relation between objective function values and decision variables. Our conclusions demonstrated that NSGA‐II performs well in multireservoir system optimization having multi‐objectives. Copyright © 2005 John Wiley & Sons, Ltd.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of complex aquifer geometry on groundwater storage in high‐elevation meadows of the Sierra Nevada Mountains,CA

Recent research has indicated that Sierra Nevada meadows are hydrologically more complex than previously considered. Improved understanding of the effects of aquifer parameters and climate change on water resources in and downstream of meadows is critically needed to effectively manage mountain meadows for ecosystem services and watershed contributions. This research investigates the roles of bedrock geometry, saturated hydraulic conductivity, and meadow gradient in affecting groundwater storage dynamics and surface‐water outflows in site‐scale high‐elevation meadows. Under current and projected lower snowpack conditions, we modeled groundwater flow in representative high‐elevation meadows considering 2 conceptual aquifer thickness models: uniform and variable thickness. Spatially, variable aquifer thicknesses interpreted from bedrock depths (0–28 m) were identified from a high‐resolution ground‐penetrating radar survey conducted at Tuolumne Meadows, CA. Our interpreted bedrock surface indicated several buried U‐shaped valleys including a buried ridge that separates 2 U‐shaped valleys. Groundwater flow simulations show that an increase in meadow gradient and hydraulic conductivity led to a decrease in seasonal storage and an increase in surface‐water outflow. However, models with varying bedrock geometries change the magnitude and timing of these processes. Uniform thickness models overestimated storage at the model edges and resulted in higher projected volumes of water being released to streams earlier than previously observed.     

Non‐linear system identification of the versatile‐typed structures by a novel signal processing technique

Non‐linear structural identification problems have raised considerable research efforts since decades, in which the Bouc–Wen model is generally utilized to simulate non‐linear structural constitutive characteristic. Support vector regression (SVR), a promising data processing method, is studied for versatile‐typed structural identification. First, a model selection strategy is utilized to determine the unknown power parameter of the Bouc–Wen model. Meanwhile, optimum SVR parameters are selected automatically, instead of tuning manually. Consequently, the non‐linear structural equation is rewritten in linear form, and is solved by the SVR technique. A five‐floor versatile‐type structure is studied to show the effectiveness of the proposed method, in which both power parameter known and unknown cases are investigated. Copyright © 2007 John Wiley & Sons, Ltd.     

Near‐infrared spectroscopy of a hydroecological indicator: new tool for determining sustainable yield for Floridan aquifer system

Pond‐cypress (Taxodium ascendens Brong.) is a dominant canopy species in depressional wetlands of the south‐eastern Coastal Plain. Unsustainable withdrawals from the karst Floridan aquifer system have caused premature decline and death of pond‐cypress trees, presumably owing to altered hydroperiods (which alter the flow of water and nutrients in trees). There has been no scientifically based means to determine sustainable yield from this regional aquifer system or to detect early stages of physical/ecological damage associated with groundwater mining and aquifer storage and recovery (ASR, which also can alter natural hydroperiods). In this study, the relationship between visual symptoms (indicators) of stress or premature decline, and spectral reflectance was evaluated using dried, milled branch tips collected from natural stands of mature pond‐cypress. Depressional systems evaluated represented four of the six aquifer system subregions where subsurface perturbations from groundwater mining: (i) were presumed not to be occurring (reference wetlands); (ii) may be occurring but are not documented; and (iii) have been confirmed. Sampled trees were assigned to one of three stress classes (1, no/minimal; 2, moderate; 3, severe) based on the visual indicators. Partial least squares–linear discriminant analysis of second derivative spectral transformations in the visible/shortwave near‐infrared (NIR) region (400–1100 nm) and the NIR region (1100–2500 nm) was used to evaluate the samples in assigned classes. Class 1 samples were discriminated from combined class 2 and 3 samples in the NIR region with 100% and 97% accuracy for consecutive winter sample periods (before bud‐break). The percentage of correctly classified samples in this spectral region was lower (85%) for summer samples (full leaf‐out). Second‐derivative models for the NIR region developed from the winter data sets predicted assigned classes for alternate winter's samples with an accuracy of 97% and 100%. High correlation between spectral reflectance of dried, milled branch tips collected from mature pond‐cypress in winter and visual indicators of premature decline suggests in situ pond‐cypress are hydroecological indicators of anthropogenic subsurface hydroperiod perturbations. This approach provides objective means for early detection of unsustainable aquifer yield and adverse impacts from ASR activities in the south‐eastern Coastal Plain. Used in conjunction with hydrological monitoring and modelling, the hydroecological indicators should provide the means with which sustainable yield in the south‐eastern Coastal Plain can be achieved and maintained. Copyright © 2003 John Wiley & Sons, Ltd.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Fluctuation of groundwater in an urban coastal city of India: a GIS‐based approach

Groundwater is the most important and valuable natural resources especially in coastal urban environment where surface water is insufficient to satisfy the water requirement. Puri city is located on the east coast of India where groundwater is the only source available to meet city water supply. As the city is situated on the sandy aquifer, quality of groundwater is deteriorating because of anthropogenic activities, lack of sewerage system, etc. The objective of the study was to assess the groundwater fluctuation during post‐monsoon and summer with respect to hydrogeological conditions, topography, and groundwater consumption pattern of the city. For this assessment and analysis, Geographic Information System (GIS) was used to visualize topography of the area through digital elevation model (DEM) and distribution of groundwater contours spatially and temporally. The probable areas prone to contamination were identified based on aquifer property and depths to water table below ground. Copyright © 2010 John Wiley & Sons, Ltd.     

Design of store‐release covers to minimize deep drainage in the mining and waste‐disposal industries: results from a modelling analyses based on ecophysiological principles

Sustainable long‐term storage of municipal waste and waste rock from mining activities in waste dumps (either above or below the land surface) requires minimization of percolation of rainwater into and then through stored waste material. There has been increasing attention given to the use of store‐release covers (transpirational covers) to achieve this. However, the design of such covers remains problematic because of the unique combinations of weather, vegetation composition, soils and their interactions that determine the efficacy of each design that could be available for the construction of the covers. The aim of the work described here was to use ecophysiological knowledge of soil‐plant‐atmosphere (SPA) interactions through the application of a detailed mechanistic model of the SPA continuum. We examined the relative influence of soil depth, soil texture, leaf area index and rainfall as determinants of rates of evapotranspiration and water budget for several different theoretical cover designs. We show that minimizing deep drainage requires a cover that has the following attributes: (i) a water storage capacity that is large enough to store the volume of water that is received as rainfall in above‐average wet months/seasons; (ii) a root distribution that explores the entire depth of the cover; (iii) a leaf area index that is present all year sufficient to evapotranspire monthly rainfall; and (iv) takes into account the intra‐annual and inter‐annual variability in rainfall and other climatic variables that drive ET. Copyright © 2012 John Wiley & Sons, Ltd.     

Parameter identification of framed structures using an improved finite element model‐updating method—Part I: formulation and verification

In this study, we formulate an improved finite element model‐updating method to address the numerical difficulties associated with ill conditioning and rank deficiency. These complications are frequently encountered model‐updating problems, and occur when the identification of a larger number of physical parameters is attempted than that warranted by the information content of the experimental data. Based on the standard bounded variables least‐squares (BVLS) method, which incorporates the usual upper/lower‐bound constraints, the proposed method (henceforth referred to as BVLSrc) is equipped with novel sensitivity‐based relative constraints. The relative constraints are automatically constructed using the correlation coefficients between the sensitivity vectors of updating parameters. The veracity and effectiveness of BVLSrc is investigated through the simulated, yet realistic, forced‐vibration testing of a simple framed structure using its frequency response function as input data. By comparing the results of BVLSrc with those obtained via (the competing) pure BVLS and regularization methods, we show that BVLSrc and regularization methods yield approximate solutions with similar and sufficiently high accuracy, while pure BVLS method yields physically inadmissible solutions. We further demonstrate that BVLSrc is computationally more efficient, because, unlike regularization methods, it does not require the laborious a priori calculations to determine an optimal penalty parameter, and its results are far less sensitive to the initial estimates of the updating parameters. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of base motion in instrumented steel buildings using output‐only system identification

A procedure to estimate the seismic motion at the base of a building from measured acceleration response at two or more floors is presented. The proposed method is comprised of two steps. In the first step, the dynamic characteristics of the building are inferred by using an output‐only system identification procedure. In the second step, the motion of the base of the building is estimated by using the transfer function of a simplified building model consisting of a shear and flexural continuous beam together with dynamic properties obtained in the first step. The proposed method is validated first with an analytical model subjected to the 1940 El Centro ground motion and then with an instrumented building in California that experienced the 1994 Northridge earthquake, and the ground motions at the base of the building are available. It is shown that the proposed method is capable of providing very good estimates of the motion at the base. The use of the proposed method is finally illustrated on an instrumented building, where the sensor at the base of the building did not function during the 1994 Northridge earthquake. Copyright © 2013 John Wiley & Sons, Ltd.     

Efficacy of neural network and genetic algorithm techniques in simulating spatio‐temporal fluctuations of groundwater

Groundwater modelling has emerged as a powerful tool to develop a sustainable management plan for efficient groundwater utilization and protection of this vital resource. This study deals with the development of five hybrid artificial neural network (ANN) models and their critical assessment for simulating spatio‐temporal fluctuations of groundwater in an alluvial aquifer system. Unlike past studies, in this study, all the relevant input variables having significant influence on groundwater have been considered, and the hybrid ANN technique [ANN‐cum‐Genetic Algorithm (GA)] has been used to simulate groundwater levels at 17 sites over the study area. The parameters of the ANN models were optimized using a GA optimization technique. The predictive ability of the five hybrid ANN models developed for each of the 17 sites was evaluated using six goodness‐of‐fit criteria and graphical indicators, together with adequate uncertainty analyses. The analysis of the results of this study revealed that the multilayer perceptron Levenberg–Marquardt model is the most efficient in predicting monthly groundwater levels at almost all of the 17 sites, while the radial basis function model is the least efficient. The GA technique was found to be superior to the commonly used trial‐and‐error method for determining optimal ANN architecture and internal parameters. Of the goodness‐of‐fit statistics used in this study, only root‐mean‐squared error, r² and Nash–Sutcliffe efficiency were found to be more powerful and useful in assessing the performance of the ANN models. It can be concluded that the hybrid ANN modelling approach can be effectively used for predicting spatio‐temporal fluctuations of groundwater at basin or subbasin scales. Copyright © 2014 John Wiley & Sons, Ltd.