Understanding the Aquifier : 3D Visualisation and Groundwater Time-Series Animation in Lockyer Valley

The Lockyer Valley, southeast Queensland, hosts intensive irrigated agriculture using groundwater from over 5000 alluvial bores. A current project is considering introduction of PRW (purified recycled water) to augment groundwater supplies. To assess this, a valley-wide MODFLOW simulation model is being developed plus a new unsaturated zone flow model. To underpin these models and provide a realistic understanding of the aquifer framework a 3D visualisation model has been developed using Groundwater Visualisation System (GVS) software produced at QUT.

Geotermometria das águas hipertermais do Sistema Aquífero Guarani no estado de São Paulo

Several geothermometers have been employed to determine the expected temperature in the groundwater reservoirs, most part of them based on the chemical equilibrium between thermal groundwaters and minerals which constitute the aquifer framework. Based on this principle, groundwater from Guarani Aquifer System, which presents temperatures upon to 38°C, classified as hyperthermal by the Brazilian Mineral Water Code, had their chemical composition evaluated with the main purpose to determine the expected temperatures in the reservoir based on the application of several geothermometers. These groundwaters can be classified as Na-HCO3 type evolving to Na-HCO3-Cl-SO4 type due the increase in sulphate and chloride concentrations. Measured groundwater temperatures show a good correlation with the depth of the aquifer, indicating an average geothermal gradient of 27.7°C/km. The lowest geothermal gradients, bellow 25°C/km, are observed in the central portion of the study area; meanwhile the highest gradients are observed in NE and W portions of the study area. Quartz equilibrium based geothermometers return temperature over the measured temperatures, meanwhile chalcedony geothermometer returns temperatures near the measured in the wells, except for the wells located in the area which presents the lowest geothermal gradients. The geothermometers based on the equilibrium between Na/K have returned temperatures diverse from the measured at the field, meanwhile the geothermometer based on the equilibrium between Na/K/Ca has presented good match with groundwater which present temperatures under 50°C, due the thermodynamic equilibrium between these waters and feldspars. Geothermometers show evidences about mixtures between groundwaters from GAS and underlying units, especially in wells located in the SW portion of the study area.

Hidrogeologia e hidroquímica do Sistema Aqüífero Guarani na porção ocidental da bacia sedimentar do Paraná

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

An analytical framework for quantifying aquifer response time scales associated with transient boundary conditions

A major challenge in studying coupled groundwater and surface-water interactions arises from the considerable difference in the response time scales of groundwater and surface-water systems affected by external forcings. Although coupled models representing the interaction of groundwater and surface-water systems have been studied for over a century, most have focused on groundwater quantity or quality issues rather than response time. In this study, we present an analytical framework, based on the concept of mean action time (MAT), to estimate the time scale required for groundwater systems to respond to changes in surface-water conditions. MAT can be used to estimate the transient response time scale by analyzing the governing mathematical model. This framework does not require any form of transient solution (either numerical or analytical) to the governing equation, yet it provides a closed form mathematical relationship for the response time as a function of the aquifer geometry, boundary conditions, and flow parameters. Our analysis indicates that aquifer systems have three fundamental time scales: (i) a time scale that depends on the intrinsic properties of the aquifer; (ii) a time scale that depends on the intrinsic properties of the boundary condition, and; (iii) a time scale that depends on the properties of the entire system. We discuss two practical scenarios where MAT estimates provide useful insights and we test the MAT predictions using new laboratory-scale experimental data sets.

Relation between sedimentary framework and hydrogeology in the Guarani Aquifer System in Sao Paulo state, Brazil

This paper presents the results of a new investigation of the Guarani Aquifer System (SAG) in Sao Paulo state. New data were acquired about sedimentary framework, flow pattern, and hydrogeochemistry. The flow direction in the north of the state is towards the southwest and not towards the west as expected previously. This is linked to the absence of SAG outcrop in the northeast of Sao Paulo state. Both the underlying Piramboia Formation and the overlying Botucatu Formation possess high porosity (18.9% and 19.5%, respectively), which was not modified significantly by diagenetic changes. Investigation of sediments confirmed a zone of chalcedony cement close to the SAG outcrop and a zone of calcite cement in the deep confined zone. The main events in the SAG post-sedimentary history were: (1) adhesion of ferrugineous coatings on grains, (2) infiltration of clays in eodiagenetic stage, (3) regeneration of coatings with formation of smectites, (4) authigenic overgrowth of quartz and K-feldspar in advanced eodiagenetic stage, (5) bitumen cementation of Piramboia Formation in mesodiagenetic stage, (6) cementation by calcite in mesodiagenetic and telodiagenetic stages in Piramboia Formation, (7) formation of secondary porosity by dissolution of unstable minerals after appearance of hydraulic gradient and penetration of the meteoric water caused by the uplift of the Serra do Mar coastal range in the Late Cretaceous, (8) authigenesis of kaolinite and amorphous silica in unconfined zone of the SAG and cation exchange coupled with the dissolution of calcite at the transition between unconfined and confined zone, and (9) authigenesis of analcime in the confined SAG zone. The last two processes are still under operation. The deep zone of the SAG comprises an alkaline pH, Na-HCO(3) groundwater type with old water and enriched delta(13)C values (<-3.9), which evolved from a neutral pH, Ca-HCO(3) groundwater type with young water and depleted delta(13)C values (>-18.8) close to the SAG outcrop. This is consistent with a conceptual geochemical model of the SAG, suggesting dissolution of calcite driven by cation exchange, which occurs at a relatively narrow front recently moving downgradient at much slower rate compared to groundwater flow. More depleted values of delta(18)O in the deep confined zone close to the Parana River compared to values of relative recent recharged water indicate recharge occur during a period of cold climate. The SAG is a ""storage-dominated"" type of aquifer which has to be managed properly to avoid its overexploitation. (C) 2011 Elsevier Ltd. All rights reserved.

Three-dimensional geological modelling and multivariate statistical analysis of water chemistry data to analyse and visualise aquifer structure and groundwater composition in the Wairau Plain, Marlborough District, New Zealand

Concerns regarding groundwater contamination with nitrate and the long-term sustainability of groundwater resources have prompted the development of a multi-layered three dimensional (3D) geological model to characterise the aquifer geometry of the Wairau Plain, Marlborough District, New Zealand. The 3D geological model which consists of eight litho-stratigraphic units has been subsequently used to synthesise hydrogeological and hydrogeochemical data for different aquifers in an approach that aims to demonstrate how integration of water chemistry data within the physical framework of a 3D geological model can help to better understand and conceptualise groundwater systems in complex geological settings. Multivariate statistical techniques(e.g. Principal Component Analysis and Hierarchical Cluster Analysis) were applied to groundwater chemistry data to identify hydrochemical facies which are characteristic of distinct evolutionary pathways and a common hydrologic history of groundwaters. Principal Component Analysis on hydrochemical data demonstrated that natural water-rock interactions, redox potential and human agricultural impact are the key controls of groundwater quality in the Wairau Plain. Hierarchical Cluster Analysis revealed distinct hydrochemical water quality groups in the Wairau Plain groundwater system. Visualisation of the results of the multivariate statistical analyses and distribution of groundwater nitrate concentrations in the context of aquifer lithology highlighted the link between groundwater chemistry and the lithology of host aquifers. The methodology followed in this study can be applied in a variety of hydrogeological settings to synthesise geological, hydrogeological and hydrochemical data and present them in a format readily understood by a wide range of stakeholders. This enables a more efficient communication of the results of scientific studies to the wider community.

How long does it take for aquifer recharge or aquifer discharge processes to reach steady state?

Groundwater flow models are usually characterized as being either transient flow models or steady state flow models. Given that steady state groundwater flow conditions arise as a long time asymptotic limit of a particular transient response, it is natural for us to seek a finite estimate of the amount of time required for a particular transient flow problem to effectively reach steady state. Here, we introduce the concept of mean action time (MAT) to address a fundamental question: How long does it take for a groundwater recharge process or discharge processes to effectively reach steady state? This concept relies on identifying a cumulative distribution function, $F(t;x)$, which varies from $F(0;x)=0$ to $F(t;x) \to \infty$ as $t\to \infty$, thereby providing us with a measurement of the progress of the system towards steady state. The MAT corresponds to the mean of the associated probability density function $f(t;x) = \dfrac{dF}{dt}$, and we demonstrate that this framework provides useful analytical insight by explicitly showing how the MAT depends on the parameters in the model and the geometry of the problem. Additional theoretical results relating to the variance of $f(t;x)$, known as the variance of action time (VAT), are also presented. To test our theoretical predictions we include measurements from a laboratory–scale experiment describing flow through a homogeneous porous medium. The laboratory data confirms that the theoretical MAT predictions are in good agreement with measurements from the physical model.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

Modelling solute and microorganism transport through a peri-alpine gravel aquifer by considering the aquifer as a dual porosity fractured medium.

Approach:
In-situ passive gradient comparative artificial tracer testing, undertaken using solutes (Uranine and Iodide), Bacteria (E.coli and P.putida) and bacteriophage (H40/1), permitted comparison of the mobility of different sized microorganisms relative to solutes in the sand and gravel aquifer underlying Dornach, Germany.
Tracer breakthrough curves reveal that even though uranine initially arrived at observation wells at the same time as microbiological tracers, maximum relative concentrations were sometimes less than those of microbiological tracers, while solute breakthrough curves proved more disperse.
Monitoring uranine breakthrough with depth suggested tracers arrived in observation wells in discrete 0.5m-1m thick intervals, over the aquifer’s 12m saturated thickness. Nearby exposures of aquifer material suggested that the aquifer consisted of sandy gravels enveloping sequences of open framework (OW) gravel up to 1m thick. Detailed examination of OW units revealed that they contained lenses of silty sand up to 1m long x 30cm thick., while granulometric data suggested that the gravel was two to three orders of magnitude more permeable than the enveloping sandy gravel.
Solute and microorganism tracer responses could not be simulated using conventional advective-dispersive equation solutions employing the same velocity and dispersion terms. By contrast solute tracer responses, modelled using a dual porosity approach for fractured media (DP-1D) corresponded well to observed field data. Simulating microorganism responses using the same transport terms, but no dual porosity term, generated good model fits and explained the higher relative concentration of the bacteria, compared to the non-reactive solute, even with first order removal to account for lower RR. Geologically, model results indicate that the silty units within open framework gravels are accessible to solute tracers, but not to microorganisms.
Importance:
Results highlight the benefits of geological observations developing appropriate conceptual models of solute and micro organism transport and in developing suitable numerical approaches to quantifying microorganism mobility at scales appropriate for the development of groundwater supply (wellhead) protection zones.

Using solute/ micro organism comparative tracer testing to investigate hydrogeological heterogeneity in a fluvioglacial sand and gravel aquifer.

Groundwater drawn from fluvioglacial sand and gravel aquifers form the principal source of drinking water in many part of central Western Europe. High population densities and widespread organic agriculture in these same areas constitute hazards that may impact the microbiological quality of many potable supplies. Tracer testing comparing two similarly sized bacteria (E.coli and P. putida) and the smaller bacteriophage (H40/1) with the response of non-reactive solute tracer (uranine) at the decametre scale revealed that all tracers broke through up to 100 times more quickly than anticipated using conventional rules of thumb. All microbiological tracer responses were less disperse than the solute, although bacterial peak relative concentrations consistently exceeded those of the solute tracer at one sampling location reflecting exclusion processes influencing micro biological tracer migration. Relative recoveries of H40/1 and E.coli proved consistent at both monitoring wells, while responses of H40/1 and P.putida differed. Examination of exposures of the upper reaches of the aquifer in nearby sand and gravel quarries revealed the aquifer to consist of laterally extensive layers of open framework (OW) gravel enveloped in finer grained gravelly sand. Granulometric analysis of these deposits suggested that the OW gravel was up to two orders of magnitude more permeable than the surrounding deposits giving rise to the preferential flow paths. By contrast fine grained lenses of silty sand within the OW gravels are suspected to play an important role in the exclusion processes that permit solutes to access them but exclude larger micro organisms.

Characterizing Aquifer Heterogeneity Using Bacterial and Bacteriophage Tracers

Gravel aquifers act as important potable water sources in central western Europe yet they are subject to numerous contamination pressures. Compositional and textural heterogeneity makes protection zone delineation around groundwater supplies in these units challenging; artificial tracer testing aids characterization. This paper reappraises previous tracer test results in light of new geological and microbiological data. Comparative passive gradient testing, using a fluorescent solute (Uranine), virus (H40/1 bacteriophage), and comparably sized bacterial tracers Escherichia coli and Pseudomonas putida, was used to investigate a calcareous gravel aquifer’s ability to remove microbiological contaminants at a test site near Munich, Germany. Test results revealed E. coli relative recoveries could exceed those of H40/1 at monitoring wells 10 m and 20 m from an injection well by almost four times; P. putida recoveries varied by a factor of up to three between wells. Application of filtration theory suggested greater attenuation of H40/1 relative to similarly charged E. coli occurred due to differences in microorganism size, while estimated collision efficiencies appeared comparable. By contrast, more positively charged P. putida experienced greater attenuation at one monitoring point, while lower attenuation rates at the second location indicated the influence of geochemical heterogeneity. Test findings proved consistent with observations from nearby fresh outcrops that suggested thin open framework gravel beds dominated mass transport in the aquifer, while discrete intervals containing stained clasts reflect localized geochemical heterogeneity. Study results highlight the utility of reconciling outcrop observations with artificial tracer test responses, using microbiological tracers with well-defined properties, to characterize aquifer heterogeneity.

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Groundwater geochemical evolution in the northern portion of the Guarani Aquifer System (Brazil) and its relationship to diagenetic features

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stochastic simulation of time-series models combined with geostatistics to predict water-table scenarios in a Guarani Aquifer System outcrop area, Brazil

Stochastic methods based on time-series modeling combined with geostatistics can be useful tools to describe the variability of water-table levels in time and space and to account for uncertainty. Monitoring water-level networks can give information about the dynamic of the aquifer domain in both dimensions. Time-series modeling is an elegant way to treat monitoring data without the complexity of physical mechanistic models. Time-series model predictions can be interpolated spatially, with the spatial differences in water-table dynamics determined by the spatial variation in the system properties and the temporal variation driven by the dynamics of the inputs into the system. An integration of stochastic methods is presented, based on time-series modeling and geostatistics as a framework to predict water levels for decision making in groundwater management and land-use planning. The methodology is applied in a case study in a Guarani Aquifer System (GAS) outcrop area located in the southeastern part of Brazil. Communication of results in a clear and understandable form, via simulated scenarios, is discussed as an alternative, when translating scientific knowledge into applications of stochastic hydrogeology in large aquifers with limited monitoring network coverage like the GAS.

Balancing groundwater conservation and rural livelihoods under water and climate uncertainties: An integrated hydro-economic modeling framework

In arid countries worldwide, social conflicts between irrigation-based human development and the conservation of aquatic ecosystems are widespread and attract many public debates. This research focuses on the analysis of water and agricultural policies aimed at conserving groundwater resources and maintaining rurallivelihoods in a basin in Spain's central arid region. Intensive groundwater mining for irrigation has caused overexploitation of the basin's large aquifer, the degradation of reputed wetlands and has given rise to notable social conflicts over the years. With the aim of tackling the multifaceted socio-ecological interactions of complex water systems, the methodology used in this study consists in a novel integration into a common platform of an economic optimization model and a hydrology model WEAP (Water Evaluation And Planning system). This robust tool is used to analyze the spatial and temporal effects of different water and agricultural policies under different climate scenarios. It permits the prediction of different climate and policy outcomes across farm types (water stress impacts and adaptation), at basin's level (aquifer recovery), and along the policies’ implementation horizon (short and long run). Results show that the region's current quota-based water policies may contribute to reduce water consumption in the farms but will not be able to recover the aquifer and will inflict income losses to the rural communities. This situation would worsen in case of drought. Economies of scale and technology are evidenced as larger farms with cropping diversification and those equipped with modern irrigation will better adapt to water stress conditions. However, the long-term sustainability of the aquifer and the maintenance of rurallivelihoods will be attained only if additional policy measures are put in place such as the control of illegal abstractions and the establishing of a water bank. Within the policy domain, the research contributes to the new sustainable development strategy of the EU by concluding that, in water-scarce regions, effective integration of water and agricultural policies is essential for achieving the water protection objectives of the EU policies. Therefore, the design and enforcement of well-balanced region-specific polices is a major task faced by policy makers for achieving successful water management that will ensure nature protection and human development at tolerable social costs. From a methodological perspective, this research initiative contributes to better address hydrological questions as well as economic and social issues in complex water and human systems. Its integrated vision provides a valuable illustration to inform water policy and management decisions within contexts of water-related conflicts worldwide.