Andrea Zanini

Simultaneous identification of the pollutant release history and the source location in groundwater by means of a geostatistical approach

This paper describes an innovative procedure that is able to simultaneously identify the release history and the source location of a pollutant injection in a groundwater aquifer (simultaneous release function and source location identification, SRSI). The methodology follows a geostatistical approach: it develops starting from a data set and a reliable numerical flow and transport model of the aquifer. Observations can be concentration data detected at a given time in multiple locations or a time series of concentration measurements collected at multiple locations. The methodology requires a preliminary delineation of a probably source area and results in the identification of both the sub-area where the pollutant injection has most likely originated, and in the contaminant release history. Some weak hypotheses have to be defined about the statistical structure of the unknown release function such as the probability density function and correlation structure. Three case studies are discussed concerning two-dimensional, confined aquifers with strongly non-uniform flow fields. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix and the recent step input function procedure has been successfully applied.

Laboratory sandbox validation of pollutant source location methods

Inverse methods can be used to recover the pollutant source location from concentration data. In this paper, the relative effectiveness of two proposed methods, simultaneous release function and source location identification (SRSI) and backward probability model based on adjoint state method (BPM-ASM) are evaluated using real data collected by using experimental equipment. The device is a sandbox that reproduces an unconfined aquifer in which all the variables are controlled. A numerical model was calibrated using experimental observations. The SRSI is a stochastic procedure which finds the source location and the release history by means of a Bayesian geostatistical approach (GA). The BPM-ASM provides the backward probability location of the pollutant detected at a monitoring point by means of a reverse transport simulation. The results show that both methods perform well. While the simultaneous release function and SRSI method requires a preliminary delineation of a probable source area and some weak hypotheses about the statistical structure of the unknown release function, the backward probability model requires some hypothesis about the contaminant release time. A case study was performed using two observation points only, and despite the scarcity of data, both methodologies were able to accurately reconstruct the true source location. The GA has the advantage to recover the release history function too, whilst the backward probability model works well with fewer data. If there are many observations, both methodologies may be computationally heavy. A transfer function approach has been adopted for the numerical definition of the sensitivity matrix in the SRSI method. The reliability of the experimental equipment was tested in previous laboratory works, conducted under several different conditions.

Evaluation of dispersivity coefficients by means of a laboratory image analysis

This paper describes the application of an innovative procedure that allows the estimation of longitudinal and transverse dispersivities in an experimental plume devised in a laboratory sandbox. The phenomenon of transport in porous media is studied using sodium fluorescein as tracer. The fluorescent excitation was achieved by using blue light and the concentration data were obtained through the processing of side wall images collected with a high resolution color digital camera. After a calibration process, the relationship between the luminosity of the emitted fluorescence and the fluorescein concentration was determined at each point of the sandbox. The relationships were used to describe the evolution of the transport process quantitatively throughout the entire domain. Some check tests were performed in order to verify the reliability of the experimental device. Numerical flow and transport models of the sandbox were developed and calibrated comparing computed and observed flow rates and breakthrough curves. The estimation of the dispersivity coefficients was carried out by analyzing the concentration field deduced from the images collected during the experiments; the dispersivity coefficients were evaluated in the domain zones where the tracer affected the porous medium under the hypothesis that the transport phenomenon is described by advection-dispersion equation (ADE) and by computing the differential components of the concentration by means of a numerical leap-frog scheme. The values determined agree with the ones referred in literature for similar media and with the coefficients obtained by calibrating the numerical model. Very interesting considerations have been made from the analysis of the performance of the methodology at different locations in the flow domain and phases of the plume evolution.

A conceptual hydrogeological model of ophiolitic aquifers (serpentinised peridotite): The test example of Mt. Prinzera (Northern Italy)

&NA; The main aim of this study is the experimental analysis of the hydrogeological behaviour of the Mt. Prinzera ultramafic massif in the northern Apennines, Italy. The analysed multidisciplinary database has been acquired through (a) geologic and structural survey; (b) geomorphologic survey; (c) hydrogeological monitoring; (d) physico‐chemical analyses; and (e) isotopic analyses. The ultramafic medium is made of several lithological units, tectonically overlapped. Between them, a low‐permeability, discontinuous unit has been identified. This unit behaves as an aquitard and causes a perched groundwater to temporary flow within the upper medium, close to the surface. This perched groundwater flows out along several structurally controlled depressions, and then several high‐altitude temporary springs can be observed during recharge, together with several perennial basal (i.e., low altitude) springs, caused by the compartmentalisation of the system because of high‐angle tectonic discontinuities.

Contaminant release history identification in 2-D heterogeneous aquifers through a minimum relative entropy approach

The minimum relative entropy (MRE) method has been applied in a wide variety of fields since it was first introduced. Woodbury and Ulrych (Water Resour Res 29(8): 2847–2860, 1993, Water Resour Res 32(9): 2671–2681, 1996) adopted and improved this method to solve linear inverse problems in aquifers. In this work, the MRE method was improved to detect the source release history in 2-D aquifer characterized by a non-uniform flow-field. The approach was tested on two cases: a 2-D homogeneous conductivity field and a heterogeneous one (the hydraulic conductivity presents three orders of magnitude in terms of variability). In the latter case the transfer function cannot be described with an analytical formulation, thus, the transfer functions were estimated by means of a numerical procedure. In order to analyze the method performance in different conditions, two datasets have been used: observations collected at the same time at 20 different monitoring points, and observations collected at 2 monitoring points at several times. The observed data have been processed with and without a random error and the Boxcar and Gaussian probability distribution functions were considered as a priori information. The agreement between the true and the estimated data has been evaluated through the calculation of the normalized Root Mean Square error. The approach was able to recover the release history even in the most difficult case.

A Python Script to Compute Isochrones for MODFLOW

MODFLOW constitutes today the most popular modeling tool in the study of water flow in aquifers and in modeling aquifers. To simplify the interface to MODFLOW various GUI have been developed for the creation of model definition files and for the visualization and interpretation of results. Recently Bakker et al. (2016) developed the FloPy interface to MODFLOW that allows to import and use the produced simulation data using Python. This allows to construct model input files, run the models, read and plot simulations results through Python scripts. In this note, we present a Python program (that uses FloPy) interface that allows us to generate time-related capture zones (isochrones) for confined 2D steady-state groundwater flow in unbounded domains, with one or more wells. As an application, we show a validation of the approach and the results of four basic test cases: a homogenous aquifer with one well, a heterogeneous aquifer with one well, an aquifer with four wells located both longitudinal and perpendicular to the flow direction.

A multidisciplinary procedure to evaluate and optimize the efficacy of hydraulic barriers in contaminated sites: a case study in Northern Italy

In the last decades, hydraulic barriers have been activated in a large number of polluted sites with the aim of preventing groundwater pollution outside the contaminated area. From a regulatory point of view, there is the need of evaluating the efficacy of these barriers. For this reason, the goal of the present study is to apply a coupled experimental modelling approach aimed at evaluating the efficacy of the barrier and providing management strategies. In particular, a case study in Italy is investigated. The study case is of main interest because of its complexity due to a heterogeneous aquifer and the presence of surface water that interacts with the below aquifer. The study has been carried out through the experimental characterization of the aquifer system (coupling the classic stratigraphic techniques with the results of radiocarbon dating, as well as through pumping and injection tests) and its hydrogeological behaviour (by means of hydraulic- and the stream-head measurements, as well as some isotopic investigations), and the implementation of a numerical model (through MODFLOW 2005). The results show the effectiveness of the coupled experimental modelling approach to analyse and simulate the hydrodynamics within the test aquifer system, as well as to evaluate the efficacy of the hydraulic barrier. Based on the results of the numerical model, easy solutions were designed to manage the barrier.

Hydrodynamics in evaporate-bearing fine-grained successions investigated through an interdisciplinary approach : A test study in southern Italy-hydrogeological behaviour of heterogeneous low-permeability media

Messinian evaporates are widely distributed in the Mediterranean Sea as outcropping sediments in small marginal basins and in marine cores. Progressive filling of subbasins led to the formation of complex aquifer systems in different regions where hypersaline and fresh water coexist and interact in different manner. It also generates a significant diversification of groundwater hydrochemical signature and different microbial communities. In the case study, the hydrogeology and hydrochemistry of the whole system are influenced by good hydraulic connection between the shallower pyroclastic horizon and the underlying evaporate-bearing fine-grained Messinian succession. This is demonstrated by the merge of hydrogeological, chemical, isotopic, and microbiological data. No mixing with deep ascending waters has been observed. As shown by geophysical, hydraulic, and microbiological investigations, the hydraulic heterogeneity of the Messinian bedrock, mainly due to karstified evaporitic interstrata/lenses, causes the hydraulic head to significantly vary with depth. Somewhere, the head increases with the depth’s increase and artesian flow conditions are locally observed. Moreover, the metagenomic investigations demonstrated the existence of a poor hydraulic connection within the evaporate-bearing fine-grained succession at metric and decametric scales, therefore leading to a patchwork of geochemical (and microbiological) subenvironments.