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.

Bacterially mediated iron cycling and associated biogeochemical processes in a subtropical shallow coastal aquifer: implications for groundwater quality

Bacterially mediated iron redox cycling exerts a strong influence on groundwater geochemistry, but few studies have investigated iron biogeochemical processes in coastal alluvial aquifers from a microbiological viewpoint. The shallow alluvial aquifer located adjacent to Poona estuary on the subtropical Southeast Queensland coast represents a redox-stratified system where iron biogeochemical cycling potentially affects water quality. Using a 300 m transect of monitoring wells perpendicular to the estuary, we examined groundwater physico-chemical conditions and the occurrence of cultivable bacterial populations involved in iron (and manganese, sulfur) redox reactions in this aquifer. Results showed slightly acidic and near-neutral pH, suboxic conditions and an abundance of dissolved iron consisting primarily of iron(II) in the majority of wells. The highest level of dissolved iron(III) was found in a well proximal to the estuary most likely a result of iron curtain effects due to tidal intrusion. A number of cultivable, (an)aerobic bacterial populations capable of diverse carbon, iron, or sulfur metabolism coexisted in groundwater redox transition zones. Our findings indicated aerobic, heterotrophic respiration and bacterially mediated iron/sulfur redox reactions were integral to carbon cycling in the aquifer. High abundances of dissolved iron and cultivable iron and sulfur bacterial populations in estuary-adjacent aquifers have implications for iron transport to marine waters. This study demonstrated bacterially mediated iron redox cycling and associated biogeochemical processes in subtropical coastal groundwaters using culture-based methods.

Water–rock interactions : an investigation of the relationships between mineralogy and groundwater composition and flow in a subtropical basalt aquifer

A holistic study of the composition of the basalt groundwaters of the Atherton Tablelands region in Queensland, Australia was undertaken to elucidate possible mechanisms for the evolution of these very low salinity, silica- and bicarbonate-rich groundwaters. It is proposed that aluminosilicate mineral weathering is the major contributing process to the overall composition of the basalt groundwaters. The groundwaters approach equilibrium with respect to the primary minerals with increasing pH and are mostly in equilibrium with the major secondary minerals (kaolinite and smectite), and other secondary phases such as goethite, hematite, and gibbsite, which are common accessory minerals in the Atherton basalts. The mineralogy of the basalt rocks, which has been examined using X-ray diffraction and whole rock geochemistry methods, supports the proposed model for the hydrogeochemical evolution of these groundwaters: precipitation + CO 2 (atmospheric + soil) + pyroxene + feldspars + olivine yields H 4SiO 4, HCO 3 -, Mg 2+, Na +, Ca 2+ + kaolinite and smectite clays + amorphous or crystalline silica + accessory minerals (hematite, goethite, gibbsite, carbonates, zeolites, and pyrite). The variations in the mineralogical content of these basalts also provide insights into the controls on groundwater storage and movement in this aquifer system. The fresh and weathered vesicular basalts are considered to be important in terms of zones of groundwater occurrence, while the fractures in the massive basalt are important pathways for groundwater movement.

Use of hierarchical cluster analysis to assess the representativeness of a baseline groundwater quality monitoring network: comparison of New Zealand’s national and regional groundwater monitoring programs

Baseline monitoring of groundwater quality aims to characterize the ambient condition of the resource and identify spatial or temporal trends. Sites comprising any baseline monitoring network must be selected to provide a representative perspective of groundwater quality across the aquifer(s) of interest. Hierarchical cluster analysis (HCA) has been used as a means of assessing the representativeness of a groundwater quality monitoring network, using example datasets from New Zealand. HCA allows New Zealand's national and regional monitoring networks to be compared in terms of the number of water-quality categories identified in each network, the hydrochemistry at the centroids of these water-quality categories, the proportions of monitoring sites assigned to each water-quality category, and the range of concentrations for each analyte within each water-quality category. Through the HCA approach, the National Groundwater Monitoring Programme (117 sites) is shown to provide a highly representative perspective of groundwater quality across New Zealand, relative to the amalgamated regional monitoring networks operated by 15 different regional authorities (680 sites have sufficient data for inclusion in HCA). This methodology can be applied to evaluate the representativeness of any subset of monitoring sites taken from a larger network.

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.

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.

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.