Review of 3D model development with application to groundwater studies in New Zealand 1996 – 2011

Three dimensional geological modelling techniques have been applied since 1996 with an aim to characterise the lithological and chronological units of New Zealand’s many diverse aquifers. Models of property-scattered data have also been applied to assess physical properties of aquifers and the distribution of groundwater chemistry, including groundwater age, to inform an understanding of groundwater systems. These models, fundamental to understanding groundwater recharge, flow and discharge have found many uses as outlined in this paper.

Groundwater visualisation system (GVS) : a software framework for integrated display and interrogation of conceptual hydrogeological models, data and time-series animation

Management of groundwater systems requires realistic conceptual hydrogeological models as a framework for numerical simulation modelling, but also for system understanding and communicating this to stakeholders and the broader community. To help overcome these challenges we developed GVS (Groundwater Visualisation System), a stand-alone desktop software package that uses interactive 3D visualisation and animation techniques. The goal was a user-friendly groundwater management tool that could support a range of existing real-world and pre-processed data, both surface and subsurface, including geology and various types of temporal hydrological information. GVS allows these data to be integrated into a single conceptual hydrogeological model. In addition, 3D geological models produced externally using other software packages, can readily be imported into GVS models, as can outputs of simulations (e.g. piezometric surfaces) produced by software such as MODFLOW or FEFLOW. Boreholes can be integrated, showing any down-hole data and properties, including screen information, intersected geology, water level data and water chemistry. Animation is used to display spatial and temporal changes, with time-series data such as rainfall, standing water levels and electrical conductivity, displaying dynamic processes. Time and space variations can be presented using a range of contouring and colour mapping techniques, in addition to interactive plots of time-series parameters. Other types of data, for example, demographics and cultural information, can also be readily incorporated. The GVS software can execute on a standard Windows or Linux-based PC with a minimum of 2 GB RAM, and the model output is easy and inexpensive to distribute, by download or via USB/DVD/CD. Example models are described here for three groundwater systems in Queensland, northeastern Australia: two unconfined alluvial groundwater systems with intensive irrigation, the Lockyer Valley and the upper Condamine Valley, and the Surat Basin, a large sedimentary basin of confined artesian aquifers. This latter example required more detail in the hydrostratigraphy, correlation of formations with drillholes and visualisation of simulation piezometric surfaces. Both alluvial system GVS models were developed during drought conditions to support government strategies to implement groundwater management. The Surat Basin model was industry sponsored research, for coal seam gas groundwater management and community information and consultation. The “virtual” groundwater systems in these 3D GVS models can be interactively interrogated by standard functions, plus production of 2D cross-sections, data selection from the 3D scene, rear end database and plot displays. A unique feature is that GVS allows investigation of time-series data across different display modes, both 2D and 3D. GVS has been used successfully as a tool to enhance community/stakeholder understanding and knowledge of groundwater systems and is of value for training and educational purposes. Projects completed confirm that GVS provides a powerful support to management and decision making, and as a tool for interpretation of groundwater system hydrological processes. A highly effective visualisation output is the production of short videos (e.g. 2–5 min) based on sequences of camera ‘fly-throughs’ and screen images. Further work involves developing support for multi-screen displays and touch-screen technologies, distributed rendering, gestural interaction systems. To highlight the visualisation and animation capability of the GVS software, links to related multimedia hosted online sites are included in the references.

Natural free convection in porous media: First field documentation in groundwater

Natural free convection is a process of great importance in disciplines from hydrology to meteorology, oceanography, planetary sciences, and economic geology, and for applications in carbon sequestration and nuclear waste disposal. It has been studied for over a century - but almost exclusively in theoretical and laboratory settings, Despite its importance, conclusive primary evidence of free convection in porous media does not currently exist in a natural field setting. Here, we present recent electrical resistivity measurements from a sabkha aquifer near Abu Dhabi, United Arab Emirates, where large density inversions exist. The geophysical images from this site provide, for the first time, compelling field evidence of fingering associated with natural free convection in groundwater.

Improved modeling of groundwater recharge in agricultural watersheds using a combination of crop model and remote sensing

For improved water management and efficiency of use in agriculture, studies dealing with coupled crop-surface water-groundwater models are needed. Such integrated models of crop and hydrology can provide accurate quantification of spatio-temporal variations of water balance parameters such as soil moisture store, evapotranspiration and recharge in a catchment. Performance of a coupled crop-hydrology model would depend on the availability of a calibrated crop model for various irrigated/rainfed crops and also on an accurate knowledge of soil hydraulic parameters in the catchment at relevant scale. Moreover, such a coupled model should be designed so as to enable the use/assimilation of recent satellite remote sensing products (optical and microwave) in order to model the processes at catchment scales. In this study we present a framework to couple a crop model with a groundwater model for applications to irrigated groundwater agricultural systems. We discuss the calibration of the STICS crop model and present a methodology to estimate the soil hydraulic parameters by inversion of crop model using both ground and satellite based data. Using this methodology we demonstrate the feasibility of estimation of potential recharge due to spatially varying soil/crop matrix.

Isotope Hydrology of the Oldman River basin, southern Alberta, Canada

The partially semi-arid Oldman River basin (OMRB), located in southern Alberta (Canada), has an area of 28 200 km2, is forested in its western headwater part, and is used for agriculture in its eastern part. Hydrometric measurements indicate that flow in the Oldman River has decreased by ~34% between 1913 and 2003, and it is predicted that water withdrawals will increase in the next 20 years. The objective of this study was to determine whether isotope ratio measurements can provide further insight into the water dynamics of the Oldman River and its tributaries. Surface water samples were collected monthly between December 2000 and March 2003. Groundwater samples were taken from 58 wells during one-time sampling trips. Runoff within the OMRB is currently about 70 mm year-1, with a corresponding runoff ratio of 0Ð18. Seasonal flow characteristics are markedly different upstream and downstream of the Oldman River reservoir. Upstream, sharp increases in flow in late spring and early summer are followed by a rapid decrease to base flow levels. Downstream, a prolonged high flow peak is observed due to the storage effect of the Oldman River reservoir. The seasonal variation in the isotopic composition of surface water from upstream sites is small. This suggests that peak runoff is not predominantly generated by melting snow accumulated during the preceding winter, but mainly by relatively well-mixed young groundwater. A significant increase in the d18O and d2H values in the downstream part of the basin was observed. The increase in the isotopic values is partly due to surface water and groundwater influx with progressively higher d18O and d2H values in the eastern part, and partly due to evaporation. Hence, the combination of hydrometric data with isotope measurements yields valuable insights into the water dynamics in the OMRB that may be further refined with more intensive measurement programmes in the future.

Clara Bog: An ombotrophic groundwater dependant terrestrial ecosystem (Submitted)

Comparison of flow duration curves for a weir draining an undrained raised peat with those generated 20 years previously reveal that more recent curves reflect to be flatter with a lower Q95/Q5 ratio. Comparison of the bog topography for the same period revealed that although marginal drainage/peat reclamation had resulted in desiccation of peat around the bog margin and more frequent intense runoff, the central part of the bog had subsided to form an enclosed basin ,resulting in the creation of newly formed lakes that gave the central part of the bog an improved capacity to store, and more slowly discharge, water. Interrogation of groundwater monitoring data revealed a net decline in groundwater levels of up to three metres in the glacial tills underlying the bog associated with deepening and expansion of a marginal drain network which penetrated the base of the peat. Comparing organic carbon levels in peat the central part of the bog over a ten year period revealed an overall increase, with changes being most marked in deeper fen peat layers. These findings suggest that the decline in groundwater levels in the peat substrate resulted in an increase in effective stress in the peat causing greater subsidence in the central part of the bog due to greater overall thickness. Study results highlight how the hydrology of apparently isolated obotrophic raised bog ecosystems may be influenced by groundwater pressures in deeper deposits, and how marginal drains may have the capacity to impact areas at significant distances.

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.

Analysis of groundwater-level response to rainfall and recharge estimates in fractured hard rock aquifers, NW Ireland

Despite fractured hard rock aquifers underlying over 65% of Ireland, knowledge of key processes controlling groundwater recharge in these bedrock systems is inadequately constrained. In this study, we examined 19 groundwater-level hydrographs from two Irish hillslope sites underlain by hard rock aquifers. Water-level time-series in clustered monitoring wells completed at the subsoil, soil/bedrock interface, shallow and deep bedrocks were continuously monitored hourly over two hydrological years. Correlation methods were applied to investigate groundwater-level response to rainfall, as well as its seasonal variations. The results reveal that the direct groundwater recharge to the shallow and deep bedrocks on hillslope is very limited. Water-level variations within these geological units are likely dominated by slow flow rock matrix storage. The rapid responses to rainfall (⩽2 h) with little seasonal variations were observed to the monitoring wells installed at the subsoil and soil/bedrock interface, as well as those in the shallow or deep bedrocks at the base of the hillslope. This suggests that the direct recharge takes place within these units. An automated time-series procedure using the water-table fluctuation method was developed to estimate groundwater recharge from the water-level and rainfall data. Results show the annual recharge rates of 42–197 mm/yr in the subsoil and soil/bedrock interface, which represent 4–19% of the annual rainfall. Statistical analysis of the relationship between the rainfall intensity and water-table rise reveal that the low rainfall intensity group (⩽1 mm/h) has greater impact on the groundwater recharge rate than other groups (>1 mm/h). This study shows that the combination of the time-series analysis and the water-table fluctuation method could be an useful approach to investigate groundwater recharge in fractured hard rock aquifers in Ireland.

Assessment of groundwater contamination in an agricultural peri-urban area (NW Portugal): an integrated approach

The excessive use of pesticides and fertilisers in agriculture has generated a decrease in groundwater and surface water quality in many regions of the EU, constituting a hazard for human health and the environment. Besides, on-site sewage disposal is an important source of groundwater contamination in urban and peri-urban areas. The assessment of groundwater vulnerability to contamination is an important tool to fulfil the demands of EU Directives. The purpose of this study is to assess the groundwater vulnerability to contamination related mainly to agricultural activities in a peri-urban area (Vila do Conde, NW Portugal). The hydrogeological framework is characterised mainly by fissured granitic basement and sedimentary cover. Water samples were collected and analysed for temperature, pH, electrical conductivity, chloride, phosphate, nitrate and nitrite. An evaluation of groundwater vulnerability to contamination was applied (GOD-S, Pesticide DRASTIC-Fm, SINTACS and SI) and the potential nitrate contamination risk was assessed, both on a hydrogeological GIS-based mapping. A principal component analysis was performed to characterised patterns of relationship among groundwater contamination, vulnerability, and the hydrogeological setting assessed. Levels of nitrate above legislation limits were detected in 75 % of the samples analysed. Alluvia units showed the highest nitrate concentrations and also the highest vulnerability and risk. Nitrate contamination is a serious problem affecting groundwater, particularly shallow aquifers, especially due to agriculture activities, livestock and cesspools. GIS-based cartography provided an accurate way to improve knowledge on water circulation models and global functioning of local aquifer systems. Finally, this study highlights the adequacy of an integrated approach, combining hydrogeochemical data, vulnerability assessments and multivariate analysis, to understand groundwater processes in peri-urban areas.

Delineating environmental groundwater vulnerability and protection zones mapping in fractured rock masses

Hard‐rock watersheds commonly exhibit complex geological bedrock and morphological features. Hydromineral resources have relevant economic value for the thermal spas industry. The present study aims to develop a groundwater vulnerability approach in Caldas da Cavaca hydromineral system (Aguiar da Beira, Central Portugal) which has a thermal tradition that dates back to the late 19th century, and contribute to a better understanding of the hydrogeological conceptual site model. In this work different layers were overlaid, generating several thematic maps to arrive at an integrated framework of several key‐sectors in Caldas da Cavaca site. Thus, to accomplish a comprehensive analysis and conceptualization of the site, a multi‐technical approach was used, such as, field and laboratory techniques, where several data was collected, like geotectonics, hydrology and hydrogeology, hydrogeomorphology, hydrogeophysical and hydrogeomechanical zoning aiming the application of the so‐called DISCO method. All these techniques were successfully performed and a groundwater vulnerability to contamination assessment, based on GOD‐S, DRASTIC‐Fm, SINTACS, SI and DISCO indexes methodology, was delineated. Geographical Information Systems (GIS) technology was on the basis to organise and integrate the geodatabases and to produce all the thematic maps. This multi‐technical approach highlights the importance of groundwater vulnerability to contamination mapping as a tool to support hydrogeological conceptualisation, contributing to better decision‐making of water resources management and sustainability.

Coupled atmosphere-biophysics-hydrology models for environmental modeling

The formulation and implementation of LEAF-2, the Land Ecosystem–Atmosphere Feedback model, which comprises the representation of land–surface processes in the Regional Atmospheric Modeling System (RAMS), is described. LEAF-2 is a prognostic model for the temperature and water content of soil, snow cover, vegetation, and canopy air, and includes turbulent and radiative exchanges between these components and with the atmosphere. Subdivision of a RAMS surface grid cell into multiple areas of distinct land-use types is allowed, with each subgrid area, or patch, containing its own LEAF-2 model, and each patch interacts with the overlying atmospheric column with a weight proportional to its fractional area in the grid cell. A description is also given of TOPMODEL, a land hydrology model that represents surface and subsurface downslope lateral transport of groundwater. Details of the incorporation of a modified form of TOPMODEL into LEAF-2 are presented. Sensitivity tests of the coupled system are presented that demonstrate the potential importance of the patch representation and of lateral water transport in idealized model simulations. Independent studies that have applied LEAF-2 and verified its performance against observational data are cited. Linkage of RAMS and TOPMODEL through LEAF-2 creates a modeling system that can be used to explore the coupled atmosphere–biophysical–hydrologic response to altered climate forcing at local watershed and regional basin scales.

Groundwater favorability map using GIS multicriteria data analysis on crystalline terrain, Sao Paulo State, Brazil

This paper presents the groundwater favorability mapping on a fractured terrain in the eastern portion of Sao Paulo State, Brazil. Remote sensing, airborne geophysical data, photogeologic interpretation, geologic and geomorphologic maps and geographic information system (GIS) techniques have been used. The results of cross-tabulation between these maps and well yield data allowed groundwater prospective parameters in a fractured-bedrock aquifer. These prospective parameters are the base for the favorability analysis whose principle is based on the knowledge-driven method. The mutticriteria analysis (weighted linear combination) was carried out to give a groundwater favorabitity map, because the prospective parameters have different weights of importance and different classes of each parameter. The groundwater favorability map was tested by cross-tabulation with new well yield data and spring occurrence. The wells with the highest values of productivity, as well as all the springs occurrence are situated in the excellent and good favorabitity mapped areas. It shows good coherence between the prospective parameters and the well yield and the importance of GIS techniques for definition of target areas for detail study and wells location. (c) 2008 Elsevier B.V. All rights reserved.