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.

Understanding Connectivity between Groundwater Chemistry Data and Geological Stratigraphy via 3D Sub-surface Visualization and Analysis

This paper describes the 3D Water Chemistry Atlas - an open source, Web-based system that enables the three-dimensional (3D) sub-surface visualization of ground water monitoring data, overlaid on the local geological model. Following a review of existing technologies, the system adopts Cesium (an open source Web-based 3D mapping and visualization interface) together with a PostGreSQL/PostGIS database, for the technical architecture. In addition a range of the search, filtering, browse and analysis tools were developed that enable users to interactively explore the groundwater monitoring data and interpret it spatially and temporally relative to the local geological formations and aquifers via the Cesium interface. The result is an integrated 3D visualization system that enables environmental managers and regulators to assess groundwater conditions, identify inconsistencies in the data, manage impacts and risks and make more informed decisions about activities such as coal seam gas extraction, waste water extraction and re-use.

Impact of pit-toilet leachate on groundwater chemistry and role of vadose zone in removal of nitrate and E-coli pollutants in Kolar District, Karnataka, India

Assessment of chemistry of groundwater infiltrated by pit-toilet leachate and contaminant removal by vadose zone form the focus of this study. The study area is Mulbagal Town in Karnataka State, India. Groundwater level measurements and estimation of unsaturated permeability indicated that the leachate recharged the groundwater inside the town at the rate of 1 m/day. The average nitrate concentration of groundwater inside the town (148 mg/L) was three times larger than the permissible limit (45 mg/L), while the average nitrate concentration of groundwater outside the town (30 mg/L) was below the permissible limit. The groundwater inside the town exhibited E. coli contamination, while groundwater outside the town was free of pathogen contamination. Infiltration of alkalis (Na+, K+) and strong acids (Cl-, SO4 (2-)) caused the mixed Ca-Mg-Cl type (60 %) and Na-Cl type (28 %) facies to predominate groundwater inside the town, while, Ca-HCO3 (35 %), mixed Ca-Mg-Cl type (35 %) and mixed Ca-Na-HCO3 type (28 %) facies predominated groundwater outside/periphery of town. Reductions in E. coli and nitrate concentrations with vadose zone thickness indicated its participation in contaminant removal. A 4-m thickness of unsaturated sand + soft, disintegrated weathered rock deposit facilitates the removal of 1 log of E. coli pathogen. The anoxic conditions prevailing in the deeper layers of the vadose zone (> 19 m thickness) favor denitrification resulting in lower nitrate concentrations (28-96 mg/L) in deeper water tables (located at depths of -29 to -39 m).

Influence of anthropogenic contamination on groundwater chemistry in Mulbagal town, Kolar District, India

Geochemical studies are performed to examine the impact of leachate infiltration from on-site sewage disposal systems on the groundwater chemistry in Mulbagal town, Kolar District, Karnataka State, India. The leachate infiltration imposed nitrate concentrations ranging from 4 mg/L to 388 mg/L in the groundwater samples; it was observed that 79% of the samples exhibited nitrate concentrations in excess of drinking water permissible limit (45 mg/L). The average (of 43 measurements) E. coli levels in the groundwater samples corresponded to 189 MPN/100 mL and 55% of the samples tested exhibit pathogen contamination. Results also showed that the groundwater in the study area is characterized by acidic pH, large calcium + magnesium ion and Na/Cl ratios of < unity causing majority of the ground water samples to classify as Ca-Mg-Cl type and Na-Cl type. Saturation index (SI) computation using Visual MINTEQ program showed that the groundwater samples are under-saturated with respect to calcite. The theoretical SI values (of calcite) however suggested that the groundwater samples ought to be over-saturated with calcite. Under-saturation of the groundwater samples with calcite is attributed to increased dissolution of the mineral in the acidic environment of the groundwater.

Characterization of groundwater chemistry under the influence of lithologic and anthropogenic factors along a climatic gradient in Upper Cauvery basin, South India

Hydrogeological and climatic effect on chemical behavior of groundwater along a climatic gradient is studied along a river basin. `Semi-arid' (500-800 mm of mean annual rainfall), `sub-humid' (800-1,200 mm/year) and `humid' (1,200-1,500 mm/year) are the climatic zones chosen along the granito-gneissic plains of Kabini basin in South India for the present analysis. Data on groundwater chemistry is initially checked for its quality using NICB ratio (<+/- 5 %), EC versus TZ+ (similar to 0.85 correlation), EC versus TDS and EC versus TH analysis. Groundwater in the three climatic zones is `hard' to `very hard' in terms of Ca-Mg hardness. Polluted wells are identified (> 40 % of pollution) and eliminated for the characterization. Piper's diagram with mean concentrations indicates the evolution of CaNaHCO3 (semi-arid) from CaHCO3 (humid zone) along the climatic gradient. Carbonates dominate other anions and strong acids exceeded weak acids in the region. Mule Hole SEW, an experimental watershed in sub-humid zone, is characterized initially using hydrogeochemistry and is observed to be a replica of entire sub-humid zone (with 25 wells). Extension of the studies for the entire basin (120 wells) showed a chemical gradient along the climatic gradient with sub-humid zone bridging semi-arid and humid zones. Ca/Na molar ratio varies by more than 100 times from semi-arid to humid zones. Semi-arid zone is more silicaceous than sub-humid while humid zone is more carbonaceous (Ca/Cl similar to 14). Along the climatic gradient, groundwater is undersaturated (humid), saturated (sub-humid) and slightly supersaturated (semi-arid) with calcite and dolomite. Concentration-depth profiles are in support of the geological stratification i.e., not approximate to 18 m of saprolite and similar to 25 m of fracture rock with parent gneiss beneath. All the wells are classified into four groups based on groundwater fluctuations and further into `deep' and `shallow' based on the depth to groundwater. Higher the fluctuations, larger is its impact on groundwater chemistry. Actual seasonal patterns are identified using `recharge-discharge' concept based on rainfall intensity instead of traditional monsoon-non-monsoon concept. Non-pumped wells have low Na/Cl and Ca/Cl ratios in recharge period than in discharge period (Dilution). Few other wells, which are subjected to pumping, still exhibit dilution chemistry though water level fluctuations are high due to annual recharge. Other wells which do not receive sufficient rainfall and are constantly pumped showed high concentrations in recharge period rather than in discharge period (Anti-dilution). In summary, recharge-discharge concept demarcates the pumped wells from natural deep wells thus, characterizing the basin.

Multivariate relationships between groundwater chemistry and toxicity in an urban aquifer

Multivariate statistical methods were used to investigate file Causes of toxicity and controls on groundwater chemistry from 274 boreholes in an Urban area (London) of the United Kingdom. The groundwater was alkaline to neutral, and chemistry was dominated by calcium, sodium, and Sulfate. Contaminants included fuels, solvents, and organic compounds derived from landfill material. The presence of organic material in the aquifer caused decreases in dissolved oxygen, sulfate and nitrate concentrations. and increases in ferrous iron and ammoniacal nitrogen concentrations. Pearson correlations between toxicity results and the concentration of individual analytes indicated that concentrations of ammoinacal nitrogen, dissolved oxygen, ferrous iron, and hydrocarbons were important where present. However, principal component and regression analysis suggested no significant correlation between toxicity and chemistry over the whole area. Multidimensional Scaling was used to investigate differences in sites caused by historical use, landfill gas status, or position within the sample area. Significant differences were observed between sites with different historical land use and those with different gas status. Examination of the principal component matrix revealed that these differences are related to changes in the importance of reduced chemical species.

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.

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.

Influence of anthropogenic contamination on fluoride concentration in groundwater: A study of Mulbagal town, Kolar district, Karnataka

Groundwater contamination is a serious concern in India. Major geogenic contaminants include fluoride, arsenic and iron, while common anthropogenic contaminants include nitrate, metals, organics and microbial contamination. Besides, known point and diffuse sources, groundwater c ontamination from inf iltration of pit to ilet leachate is an emerging concern. The study area of this paper is Kolar district in Karnataka that is hot spot of fluoride contamination. The absence of fluoride contamination in Mulbagal town and the alterations in groundwater chemistry from infiltration of pit toilet leachate motivated the author to examine the possible linkages between anthropogenic contamination and fluoride concentration in groundwater of Mulbagal town. Analysis of the groundwater chemistry revealed that the groundwater in Mulbagal town is under saturated with respect to calcite that suppresses the dissolution of fluorite and the fluoride concentration in the groundwater. The slightly acidic pH of the groundwater is considered responsible to facilitate calcite dissolution under saturation.

Is there a relationship between soil and groundwater toxicity?

Part IIA of the Environmental Protection Act 1990 requires environmental regulators to assess the risk of contaminants leaching from soils into groundwater (DETR, 1999). This newly introduced legislation assumes a link between soil and groundwater chemistry, in which rainwater leaches contaminants from soil into the saturated zone. As the toxicity of both groundwater and overlying soils is dependent upon the chemicals present, their partitioning and their bioavailability, similar patterns of soil, leachates and groundwater toxicity should be observed at contaminated sites. Soil and groundwater samples were collected from different contaminated land sites in an urban area, and used to determine relationships between soil chemistry and toxicity, mobility of contaminants, and groundwater chemistry and toxicity. Soils were leached using water to mimic rainfall, and both the soils and leachates tested using bioassays. Soil bioassays were carried out using Eisenia fetida, whilst groundwater and leachates were tested using the Microtox(TM) test system and Daphnia magna 48 h acute tests. Analysis of the bioassay responses demonstrated that a number of the samples were toxic to test organisms, however, there were no significant statistical relationships between soil, groundwater and leachate toxicity. Nor were there significant correlations between soil, leachates and groundwater chemistry.

The mechanism of 234U 238U activity ratio enhancement in karstic limestone groundwater

Experimental etch/leach of Carboniferous Limestone gravels on a laboratory time-scale has demonstrated that 234U 238U activity ratios (AR's) greater than the radiochemical equilibrium value may be generated on short time-scales. The molar U/Ca and Mg/Ca ratios show that both U and Mg are leached preferentially relative to Ca whereas the molar U/Mg ratio is only slightly greater than that of the rock matrix. The generation of enhanced AR's is attributed to a two-stage process in which the limestone surface is dissolved by zero-order etch and silicate minerals so released are subjected to first-order chemical leach of U and Mg. The implications of these results for the production of enhanced AR's in Carboniferous Limestone groundwater are discussed. It is suggested that chemical leaching or exchange of U between groundwater and its particulate load or at the aquifer fluid-solid interface is an important mechanism controlling AR changes as groundwater migrates beyond a redox boundary. AR's for dissolved U in groundwater are more probably related to chemical equilibria than to groundwater age. © 1993.

Temporal and Spatial Variations of Ions, Isotopes and Agricultural Contaminants in Surface Waters and Groundwater of Nebraska's Rainwater Basin Wetland Region

The wetlands of south-central Nebraska’s Rainwater Basin region are considered of international importance as a habitat for millions of migratory birds, but are being endangered by agricultural practices. The Rainwater Basin extends across 17 counties and covers 4,000 square miles. The purpose of this study was to assemble baseline chemical data for several representative wetlands across the Rainwater Basin region, and determine the use of these chemical data for investigating groundwater recharge. Eight representative wetlands were chosen across the Rainwater Basin to monitor surface and groundwater chemistry. At each site, a shallow well and deep well were installed and sampled once in the summer of 2009 and again in the spring of 2010. Wetland surface water was sampled monthly from April, 2009 to May, 2010. Waters were analyzed for major ions, nutrients, pesticides and oxygen-18 and deuterium isotopes at the University of Nebraska Water Sciences Laboratory. Geochemical analysis of surface waters presents a range of temporal and spatial variations. Wetlands had variable water volumes, isotopic compositions, ion chemistries and agricultural contaminant levels throughout the year and, except for a few trends, theses variations cannot be predicted with certainty year-to-year or wetland-to-wetland. Isotopic compositions showed evaporation was a contributor to water loss, and thus, did impact water chemistry. Surface water nitrate concentrations ranged from <0.10 to 4.04 mg/L. The nitrate levels are much higher in the groundwater, ranging from <0.10 to 18.4 mg/L, and are of concern because they are found above the maximum contaminant level (MCL) of 10 mg/L. Atrazine concentrations in surface waters ranged from <0.05 to 10.3 ppb. Groundwater atrazine concentrations ranged from <0.05 to 0.28 ppb. The high atrazine concentrations in surface waters are of concern as they are above the MCL of 3 ppb, and the highest levels occur during the spring bird migration. Most sampled groundwaters had detectable tritium indicating a mix of modern (<5 to 10 years old) and submodern (older than 1950s) recharge. The groundwater also had differences in chemical and isotope composition, and in some cases, increased nitrate concentrations, between the two sampling periods. Modern groundwater tritium ages and changes in groundwater chemical and isotopic compositions may indicate connections with surface waters in the Rainwater Basin.

Origin of the groundwater salinization in La Aldea coastal aquifer (Gran Canaria, Canary Islands)

[EN] Groundwater chemistry in La Aldea aquifer (Gran Canaria, Canary Islands) shows high contents of chloride and nitrate ions. The salinization process has been modelled using the geochemical data, taking into account the results of a previous flow model. The results allow to identify the salinity of the recharge from the rainfall under aridity conditions and the irrigation returns like the main causes of the groundwater salinization.

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.

Host Rock Classification (HRC) system for nuclear waste disposal in crystalline bedrock

A new rock mass classification scheme, the Host Rock Classification system (HRC-system) has been developed for evaluating the suitability of volumes of rock mass for the disposal of high-level nuclear waste in Precambrian crystalline bedrock. To support the development of the system, the requirements of host rock to be used for disposal have been studied in detail and the significance of the various rock mass properties have been examined. The HRC-system considers both the long-term safety of the repository and the constructability in the rock mass. The system is specific to the KBS-3V disposal concept and can be used only at sites that have been evaluated to be suitable at the site scale. By using the HRC-system, it is possible to identify potentially suitable volumes within the site at several different scales (repository, tunnel and canister scales). The selection of the classification parameters to be included in the HRC-system is based on an extensive study on the rock mass properties and their various influences on the long-term safety, the constructability and the layout and location of the repository. The parameters proposed for the classification at the repository scale include fracture zones, strength/stress ratio, hydraulic conductivity and the Groundwater Chemistry Index. The parameters proposed for the classification at the tunnel scale include hydraulic conductivity, Q´ and fracture zones and the parameters proposed for the classification at the canister scale include hydraulic conductivity, Q´, fracture zones, fracture width (aperture + filling) and fracture trace length. The parameter values will be used to determine the suitability classes for the volumes of rock to be classified. The HRC-system includes four suitability classes at the repository and tunnel scales and three suitability classes at the canister scale and the classification process is linked to several important decisions regarding the location and acceptability of many components of the repository at all three scales. The HRC-system is, thereby, one possible design tool that aids in locating the different repository components into volumes of host rock that are more suitable than others and that are considered to fulfil the fundamental requirements set for the repository host rock. The generic HRC-system, which is the main result of this work, is also adjusted to the site-specific properties of the Olkiluoto site in Finland and the classification procedure is demonstrated by a test classification using data from Olkiluoto. Keywords: host rock, classification, HRC-system, nuclear waste disposal, long-term safety, constructability, KBS-3V, crystalline bedrock, Olkiluoto