Understanding hydrochemical processes using 3D models and hydrochemical facies in Wairau Plain groundwater, Marlborough

Three dimensional models and groundwater quality are combined to better understand and conceptualise groundwater systems in complex geological settings in the Wairau Plain, Marlborough. Hydrochemical facies, which are characteristic of distinct evolutionary pathways and a common hydrologic history of groundwaters, are identified within geological formations to assess natural water-rock interactions, redox potential and human agricultural impact on groundwater quality in the Wairau Plain.

Assessing the representativeness of a groundwater quality monitoring network

The monitoring sites comprising a state of the environment (SOE) network must be carefully selected to ensure that they will be representative of the broader resource. Hierarchical cluster analysis (HCA) is a data-driven technique that can potentially be employed to assess the representativeness of a SOE monitoring network. The objective of this paper is to explore the use of HCA as an approach for assessing the representativeness of the New Zealand National Groundwater Monitoring Programme (NGMP), which is comprised of 110 monitoring sites across the country.

Assessment of interactions between alluvial, Great Artesian Basin (GAB) and volcanic aquifers using 3D visualisation and environmental tracers, Lockyer Valley, southeast Queensland, Australia

A detailed 3D lithological model framework was developed using GOCAD software to understand interactions between alluvial, volcanic and GAB aquifers and the spatial and temporal distribution of groundwater recharge to the alluvium of the Lockyer Valley. Groundwater chemistry, isotope data (H20-δ2H and δ18O , 87Sr/86Sr, 3H and 14C) and groundwater level time-series data from approximately 550 observation wells were integrated into the catchment-wide 3D model to assess the recharge processes involved. This approach enabled the identification of zones where recharge to the alluvium primarily occurs from stream water during episodic flood events. Importantly, the study also demonstrates that in some sections of the alluvium recharge is also from storm rainfall and seepage discharge from the underlying GAB aquifers. These other sources of recharge are indicated by (a) the absence of a response of groundwater levels to flooding in some areas, (b) old radiocarbon ages, and (c) distinct bedrock water chemistry and δ2H and δ18O signatures in alluvial groundwater at these locations. Integration of isotopes, water chemistry and time-series displays of groundwater levels before and after the 2010/2011 flood into the 3D model suggest that the spatial variations in the alluvial groundwater response are mostly controlled by valley morphology and lithological (i.e. permeability) variations within the alluvium. Examination of the groundwater level variations in the 3D model also enabled quantification of the volumetric change of groundwater stored in the unconfined alluvial aquifer prior to and post-flood events.

Chemical variability of groundwater samples collected from a Coal Seam Gas exploration well, Maramarua, New Zealand

A pilot study has produced 31 groundwater samples from a coal seam gas (CSG) exploration well located in Maramarua, New Zealand. This paper describes sources of CSG water chemistry variations, and makes sampling and analytical recommendations to minimize these variations. The hydrochemical character of these samples is studied using factor analysis, geochemical modelling, and a sparging experiment. Factor analysis unveils carbon dioxide (CO2) degassing as the principal cause of sample variation (about 33%). Geochemical modelling corroborates these results and identifies minor precipitation of carbonate minerals with degassing. The sparging experiment confirms the effect of CO2 degassing by showing a steady rise in pH while maintaining constant alkalinity. Factor analysis correlates variations in the major ion composition (about 17%) to changes in the pumping regime and to aquifer chemistry variations due to cation exchange reactions with argillaceous minerals. An effective CSG water sampling program can be put into practice by measuring pH at the well head and alkalinity at the laboratory; these data can later be used to calculate the carbonate speciation at the time the sample was collected. In addition, TDS variations can be reduced considerably if a correct drying temperature of 180°C is consistently implemented.

Groundwater overuse and farm-level technical inefficiency : evidence from Sri Lanka

Extraction of groundwater for onion and other cash crop production has been increasing rapidly during the last two decades in the dry zone areas of Sri Lanka. As a result of overuse, the quantity of available groundwater is gradually declining, while water quality is deteriorating. The deteriorating water quality has a negative impact on agricultural production, especially for crops (such as onions) that are sensitive to increases in salinity levels. This issue is examined with respect to onion production in Sri Lanka. A stochastic frontier production function (SFPF) is used, in which technical efficiency and the determinants of inefficiencies are estimated simultaneously. The results show that farmers are overusing groundwater in their onion cultivation, which has resulted in decreasing yields. Factors contributing to inefficiency in production are also identified. The results have important policy implications.

Linking three-dimensional geological modelling and multivariate statistical analysis to define groundwater chemistry baseline and identify inter-aquifer connectivity within the Clarence-Moreton Basin, Southeast Queensland, Australia

The Clarence-Moreton Basin (CMB) covers approximately 26000 km2 and is the only sub-basin of the Great Artesian Basin (GAB) in which there is flow to both the south-west and the east, although flow to the south-west is predominant. In many parts of the basin, including catchments of the Bremer, Logan and upper Condamine Rivers in southeast Queensland, the Walloon Coal Measures are under exploration for Coal Seam Gas (CSG). In order to assess spatial variations in groundwater flow and hydrochemistry at a basin-wide scale, a 3D hydrogeological model of the Queensland section of the CMB has been developed using GoCAD modelling software. Prior to any large-scale CSG extraction, it is essential to understand the existing hydrochemical character of the different aquifers and to establish any potential linkage. To effectively use the large amount of water chemistry data existing for assessment of hydrochemical evolution within the different lithostratigraphic units, multivariate statistical techniques were employed.

Image-based flow visualisation (IBFV) to enhance interpretation of complex flow patterns within a shallow tidal barrier estuary

We applied a texture-based flow visualisation technique to a numerical hydrodynamic model of the Pumicestone Passage in southeast Queensland, Australia. The quality of the visualisations using our flow visualisation tool, are compared with animations generated using more traditional drogue release plot and velocity contour and vector techniques. The texture-based method is found to be far more effective in visualising advective flow within the model domain. In some instances, it also makes it easier for the researcher to identify specific hydrodynamic features within the complex flow regimes of this shallow tidal barrier estuary as compared with the direct and geometric based methods.

How much geological detail is needed for gasifier success?

The incipient Underground Coal Gasification (UCG) industry in Queensland, Australia, undertook three trial projects in two Mesozoic basins of southeast Queensland. The experiences of these three operations provide useful retrospective insight into gasifier productivity. This paper identifies key output measures of gasifier ‘success’ including output gas composition, presence of contaminants in groundwater and consistency of chamber operation. Likewise, a review of the geological and hydrogeological understanding of each site prior to gasifier commissioning was undertaken. Productivity parameters from gasification were then correlated against the level of baseline geological/hydrogeological understanding for each site. The aim of the study was to identify the optimum scope of geological and hydrogeological understanding required at the site assessment phase to ensure safe, maximum gasifier output during production phase. This approach allows identification of poor or unexpected performance that is attributable to pre-existing uncertainty. A historical review of gasifier conditions inferred from the three trial projects is presented. Hence from the Queensland experiences it is possible to identify what aspects of baseline geological understanding should be clearly understood at the site selection phase in order to limit anomalous gasifier performance and undesirable deviations, and maximise production output.

Multivariate statistical analysis of hydrochemical data to assess alluvial aquifer–stream connectivity during drought and flood : Cressbrook Creek, southeast Queensland, Australia

A catchment-scale multivariate statistical analysis of hydrochemistry enabled assessment of interactions between alluvial groundwater and Cressbrook Creek, an intermittent drainage system in southeast Queensland, Australia. Hierarchical cluster analyses and principal component analysis were applied to time-series data to evaluate the hydrochemical evolution of groundwater during periods of extreme drought and severe flooding. A simple three-dimensional geological model was developed to conceptualise the catchment morphology and the stratigraphic framework of the alluvium. The alluvium forms a two-layer system with a basal coarse-grained layer overlain by a clay-rich low-permeability unit. In the upper and middle catchment, alluvial groundwater is chemically similar to streamwater, particularly near the creek (reflected by high HCO3/Cl and K/Na ratios and low salinities), indicating a high degree of connectivity. In the lower catchment, groundwater is more saline with lower HCO3/Cl and K/Na ratios, notably during dry periods. Groundwater salinity substantially decreased following severe flooding in 2011, notably in the lower catchment, confirming that flooding is an important mechanism for both recharge and maintaining groundwater quality. The integrated approach used in this study enabled effective interpretation of hydrological processes and can be applied to a variety of hydrological settings to synthesise and evaluate large hydrochemical datasets.

The origin of a large (>3km) maar volcano by coalescence of multiple shallow craters: Lake Purrumbete maar, southeastern Australia

Lake Purrumbete maar is located in the intraplate, monogenetic Newer Volcanics Province in southeastern Australia. The extremely large crater of 3000. m in diameter formed on an intersection of two fault lines and comprises at least three coalesced vents. The evolution of these vents is controlled by the interaction of the tectonic setting and the properties of both hard and soft rock aquifers. Lithics in the maar deposits originate from country rock formations less than 300. m deep, indicating that the large size of the crater cannot only be the result of the downwards migration of the explosion foci in a single vent. Vertical crater walls and primary inward dipping beds evidence that the original size of the crater has been largely preserved. Detailed mapping of the facies distributions, the direction of transport of base surges and pyroclastic flows, and the distribution of ballistic block fields, form the basis for the reconstruction of the complex eruption history,which is characterised by alternations of the eruption style between relatively dry and wet phreatomagmatic conditions, and migration of the vent location along tectonic structures. Three temporally separated eruption phases are recognised, each starting at the same crater located directly at the intersection of two local fault lines. Activity then moved quickly to different locations. A significant volcanic hiatus between two of the three phases shows that the magmatic system was reactivated. The enlargement of especially the main crater by both lateral and vertical growth led to the interception of the individual craters and the formation of the large circular crater. Lake Purrumbete maar is an excellent example of how complicated the evolution of large, seemingly simple, circular maar volcanoes can be, and raises the question if these systems are actually monogenetic.

Estimating specific yield and transmissivity with magnetic resonance sounding in an unconfined sandstone aquifer (Niger)

The unconfined aquifer of the Continental Terminal in Niger was investigated by magnetic resonance sounding (MRS) and by 14 pumping tests in order to improve calibration of MRS outputs at field scale. The reliability of the standard relationship used for estimating aquifer transmissivity by MRS was checked; it was found that the parametric factor can be estimated with an uncertainty a parts per thousand currency sign150% by a single point of calibration. The MRS water content (theta (MRS)) was shown to be positively correlated with the specific yield (Sy), and theta (MRS) always displayed higher values than Sy. A conceptual model was subsequently developed, based on estimated changes of the total porosity, Sy, and the specific retention Sr as a function of the median grain size. The resulting relationship between theta (MRS) and Sy showed a reasonably good fit with the experimental dataset, considering the inherent heterogeneity of the aquifer matrix (residual error is similar to 60%). Interpreted in terms of aquifer parameters, MRS data suggest a log-normal distribution of the permeability and a one-sided Gaussian distribution of Sy. These results demonstrate the efficiency of the MRS method for fast and low-cost prospection of hydraulic parameters for large unconfined aquifers.

Editor's message: the sunk cost fallacy of deep drilling

At the beginning of 2008, I visited a watershed, located in Karkinatam village in the state of Karnataka, South India, where crops are intensively irrigated using groundwater. The water table had been depleted from a depth of 5 to 50 m in a large part of the area. Presently, 42% of a total of 158 water wells in the watershed are dry. Speaking with the farmers, I have been amazed to learn that they were drilling down to 500 m to tap water. This case is, of course, not isolated.

Optimizing a piezometric network in the estimation of the groundwater budget: a case study from a crystalline-rock watershed in southern India

An estimate of the groundwater budget at the catchment scale is extremely important for the sustainable management of available water resources. Water resources are generally subjected to over-exploitation for agricultural and domestic purposes in agrarian economies like India. The double water-table fluctuation method is a reliable method for calculating the water budget in semi-arid crystalline rock areas. Extensive measurements of water levels from a dense network before and after the monsoon rainfall were made in a 53 km(2)atershed in southern India and various components of the water balance were then calculated. Later, water level data underwent geostatistical analyses to determine the priority and/or redundancy of each measurement point using a cross-validation method. An optimal network evolved from these analyses. The network was then used in re-calculation of the water-balance components. It was established that such an optimized network provides far fewer measurement points without considerably changing the conclusions regarding groundwater budget. This exercise is helpful in reducing the time and expenditure involved in exhaustive piezometric surveys and also in determining the water budget for large watersheds (watersheds greater than 50 km(2)).

Analytical solutions for transient temperature distribution in a geothermal reservoir due to cold water injection

An analytical solution to describe the transient temperature distribution in a geothermal reservoir in response to injection of cold water is presented. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. Results show that heat loss to the confining rock layers plays a vital role in slowing down the cooling of the reservoir. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters. The effects of injection rate and thermal conductivity have been found to be profound on the results.