Geophysical and hydrogeological characterisation of the impacts of on-site wastewater treatment discharge to groundwater in a poorly productive bedrock acquifer

Contaminants discharging from on-site wastewater treatment systems (OSWTSs) can impact groundwater quality, threatening human health and surface water ecosystems. Risk of negative impacts becomes elevated in areas of extreme vulnerability with high water tables, where thin unsaturated intervals limit vadose zone attenuation. A combined geophysical/hydrogeological investigation into the effects of an OSWTS, located over a poorly productive aquifer (PPA) with thin subsoil cover, aimed to characterise effluent impacts on groundwater. Groundwater, sampled from piezometers down-gradient of the OSWTS percolation area displayed spatially erratic, yet temporally consistent, contaminant distributions. Electrical resistivity tomography identified an area of gross groundwater contamination close to the percolation area and, when combined with seismic refraction and water quality data, indicated that infiltrating effluent reaching the water table discharged to a deeper more permeable zone of weathered shale resting on more competent bedrock. Subsurface structure, defined by geophysics, indicated that elevated chemical and microbiological contaminant levels encountered in groundwater samples collected from piezometers, down-gradient of sampling points with lower contaminant levels, corresponded to those locations where piezometers were screened close to the weathered shale/competent rock interface; those immediately up-gradient were too shallow to intercept this interval, and thus the more impacted zone of the contaminant plume. Intermittent occurrence of faecal indicator bacteria more than 100 m down gradient of the percolation area suggested relatively short travel times. Study findings highlight the utility of geophysics as part of multidisciplinary investigations for OSWTS contaminant plume characterisation, while also demonstrating the capacity of effluent discharging to PPAs to impact groundwater quality at distance. Comparable geophysical responses observed in similar settings across Ireland suggest the phenomena observed in this study are more widespread than previously suspected.

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.

Using geostatistical Bayesian Updating to integrate airborne radiometrics and soil geochemistry to improve mapping for mineral exploration

Mineral exploration programmes around the world use data from remote sensing, geophysics and direct sampling. On a regional scale, the combination of airborne geophysics and ground-based geochemical sampling can aid geological mapping and economic minerals exploration. The fact that airborne geophysical and traditional soil-sampling data are generated at different spatial resolutions means that they are not immediately comparable due to their different sampling density. Several geostatistical techniques, including indicator cokriging and collocated cokriging, can be used to integrate different types of data into a geostatistical model. With increasing numbers of variables the inference of the cross-covariance model required for cokriging can be demanding in terms of effort and computational time. In this paper a Gaussian-based Bayesian updating approach is applied to integrate airborne radiometric data and ground-sampled geochemical soil data to maximise information generated from the soil survey, to enable more accurate geological interpretation for the exploration and development of natural resources. The Bayesian updating technique decomposes the collocated estimate into a production of two models: prior and likelihood models. The prior model is built from primary information and the likelihood model is built from secondary information. The prior model is then updated with the likelihood model to build the final model. The approach allows multiple secondary variables to be simultaneously integrated into the mapping of the primary variable. The Bayesian updating approach is demonstrated using a case study from Northern Ireland where the history of mineral prospecting for precious and base metals dates from the 18th century. Vein-hosted, strata-bound and volcanogenic occurrences of mineralisation are found. The geostatistical technique was used to improve the resolution of soil geochemistry, collected one sample per 2 km2, by integrating more closely measured airborne geophysical data from the GSNI Tellus Survey, measured over a footprint of 65 x 200 m. The directly measured geochemistry data were considered as primary data in the Bayesian approach and the airborne radiometric data were used as secondary data. The approach produced more detailed updated maps and in particular maximized information on mapped estimates of zinc, copper and lead. Greater delineation of an elongated northwest/southeast trending zone in the updated maps strengthened the potential to investigate stratabound base metal deposits.

Towards groundwater security in coastal East Africa: initiating a regional research network and integrated hydrogeologic, climatic & socio-economic observatories in coastal aquifers of the Comoros Islands, Kenya and Tanzania

African coastal regions are expected to experience the highest rates of population growth in coming decades. Fresh groundwater resources in the coastal zone of East Africa (EA) are highly vulnerable to seawater intrusion. Increasing water demand is leading to unsustainable and ill-planned well drilling and abstraction. Wells supplying domestic, industrial and agricultural needs are or have become, in many areas, too saline for use. Climate change, including weather changes and sea level rise, is expected to exacerbate this problem. The multiplicity of physical, demographic and socio-economic driving factors makes this a very challenging issue for management. At present the state and probable evolution of coastal aquifers in EA are not well documented. The UPGro project 'Towards groundwater security in coastal East Africa' brings together teams from Kenya, Tanzania, Comoros Islands and Europe to address this knowledge gap. An integrative multidisciplinary approach, combining the expertise of hydrogeologists, hydrologists and social scientists, is investigating selected sites along the coastal zone in each country. Hydrogeologic observatories have been established in different geologic and climatic settings representative of the coastal EA region, where focussed research will identify the current status of groundwater and identify future threats based on projected demographic and climate change scenarios. Researchers are also engaging with end users as well as local community and stakeholder groups in each area in order to understanding the issues most affecting the communities and searching sustainable strategies for addressing these.

Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Modelling the response of fresh groundwater to climate and vegetation changes in coral islands

In coral islands, groundwater is a crucial freshwater resource for terrestrial life, including human water supply. Response of the freshwater lens to expected climate changes and subsequent vegetation alterations is quantified for Grande Glorieuse, a low-lying coral island in the Western Indian Ocean. Distributed models of recharge, evapotranspiration and saltwater phytotoxicity are integrated into a variable-density groundwater model to simulate the evolution of groundwater salinity. Model results are assessed against field observations including groundwater and geophysical measurements. Simulations show the major control currently exerted by the vegetation with regards to the lens morphology and the high sensitivity of the lens to climate alterations, impacting both quantity and salinity. Long-term changes in mean sea level and climatic conditions (rainfall and evapotranspiration) are predicted to be responsible for an average increase in salinity approaching 140 % (+8 kg m-3) when combined. In low-lying areas with high vegetation density, these changes top +300 % (+10 kg m-3). However, due to salinity increase and its phytotoxicity, it is shown that a corollary drop in vegetation activity can buffer the alteration of fresh groundwater. This illustrates the importance of accounting for vegetation dynamics to study groundwater in coral islands.

Adsorption Study using Optimised 3D Organised Mesoporous Silica Coated with Fe and Al Oxides for Specific As(III) and As(V) Removal from Contaminated Synthetic Groundwater

This work presents the possibility of optimising 3D Organised Mesoporous Silica (OMS) coated with both iron and aluminium oxides for the optimal removal of As(III) and As(V) from synthetic contaminated water. The materials developed were fully characterised and were tested for removing arsenic in batch experiments. The effect of total Al to Fe oxides coating on the selective removal of As(III) and As(V) was studied. It was shown that 8% metal coating was the optimal configuration for the coated OMS materials in removing arsenic. The effect of arsenic initial concentration and pH, kinetics and diffusion mechanisms was studied, modelled and discussed. It was shown that the advantage of an organised material over an un-structured sorbent was very limited in terms of kinetic and diffusion under the experimental conditions. It was shown that physisorption was the main adsorption process involved in As removal by the coated OMS. Maximum adsorption capacity of 55 mg As(V).g-1 was noticed at pH 5 for material coated with 8% Al oxides while 35 mg As(V).g-1 was removed at pH 4 for equivalent material coated with Fe oxides.

Coupling ground and airborne geophysical data with upscaling techniques for regional groundwater modeling of heterogeneous aquifers: Case study of a sedimentary aquifer intruded by volcanic dykes in Northern Ireland

In highly heterogeneous aquifer systems, conceptualization of regional groundwater flow models frequently results in the generalization or negligence of aquifer heterogeneities, both of which may result in erroneous model outputs. The calculation of equivalence related to hydrogeological parameters and applied to upscaling provides a means of accounting for measurement scale information but at regional scale. In this study, the Permo-Triassic Lagan Valley strategic aquifer in Northern Ireland is observed to be heterogeneous, if not discontinuous, due to subvertical trending low-permeability Tertiary dolerite dykes. Interpretation of ground and aerial magnetic surveys produces a deterministic solution to dyke locations. By measuring relative permeabilities of both the dykes and the sedimentary host rock, equivalent directional permeabilities, that determine anisotropy calculated as a function of dyke density, are obtained. This provides parameters for larger scale equivalent blocks, which can be directly imported to numerical groundwater flow models. Different conceptual models with different degrees of upscaling are numerically tested and results compared to regional flow observations. Simulation results show that the upscaled permeabilities from geophysical data allow one to properly account for the observed spatial variations of groundwater flow, without requiring artificial distribution of aquifer properties. It is also found that an intermediate degree of upscaling, between accounting for mapped field-scale dykes and accounting for one regional anisotropy value (maximum upscaling) provides results the closest to the observations at the regional scale.

Integrating Multiple Point Statistics with Aerial Geophysical Data to assist Groundwater Flow Models

The process of accounting for heterogeneity has made significant advances in statistical research, primarily in the framework of stochastic analysis and the development of multiple-point statistics (MPS). Among MPS techniques, the direct sampling (DS) method is tested to determine its ability to delineate heterogeneity from aerial magnetics data in a regional sandstone aquifer intruded by low-permeability volcanic dykes in Northern Ireland, UK. The use of two two-dimensional bivariate training images aids in creating spatial probability distributions of heterogeneities of hydrogeological interest, despite relatively ‘noisy’ magnetics data (i.e. including hydrogeologically irrelevant urban noise and regional geologic effects). These distributions are incorporated into a hierarchy system where previously published density function and upscaling methods are applied to derive regional distributions of equivalent hydraulic conductivity tensor K. Several K models, as determined by several stochastic realisations of MPS dyke locations, are computed within groundwater flow models and evaluated by comparing modelled heads with field observations. Results show a significant improvement in model calibration when compared to a simplistic homogeneous and isotropic aquifer model that does not account for the dyke occurrence evidenced by airborne magnetic data. The best model is obtained when normal and reverse polarity dykes are computed separately within MPS simulations and when a probability threshold of 0.7 is applied. The presented stochastic approach also provides improvement when compared to a previously published deterministic anisotropic model based on the unprocessed (i.e. noisy) airborne magnetics. This demonstrates the potential of coupling MPS to airborne geophysical data for regional groundwater modelling.

Quantitative Analysis of Fault and Fracture Systems and their Impact on Groundwater Flow in Irish Bedrock Aquifers

Fault and fracture systems are the most important store and pathway for groundwater in Ireland’s bedrock aquifers, either directly as conductive flow structures, or indirectly as the locus for the development of dolomitised limestone and karst. This article presents the preliminary results of a study involving the quantitative analysis of fault and fracture systems in the broad range of Irish bedrock types and a consideration of their impact on groundwater flow. The principal aims of the project are to develop generic conceptual models for different fault/fracture systems in different lithologies and at different depths, and to link them to observed groundwater behaviour. Here we briefly describe the geometrical characteristics of the main post-Devonian fault/fracture systems controlling groundwater flow from field observations at outcrops, quarries and mines. The structures range from Lower Carboniferous normal faults through to Variscan-related faults and veins, with the most recent structures including Tertiary strike-slip faults and ubiquitous uplift-related joint systems. The geometrical characteristics of different fault/fracture systems combined with observations of groundwater behaviour in both quarry and mine localities, can be linked to general flow and transport conceptualisations of Irish fractured bedrock. Most importantly they also provide a basis for relating groundwater flow to particular fault/fracture systems and their expression with depth and within different lithological sequences, as well as their regional variability.