Ulrich Ofterdinger

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

Comparison of methods used to calculate typical threshold values for potentially toxic elements in soil

The environmental quality of land can be assessed by calculating relevant threshold values, which differentiate between concentrations of elements resulting from geogenic and diffuse anthropogenic sources and concentrations generated by point sources of elements. A simple process allowing the calculation of these typical threshold values (TTVs) was applied across a region of highly complex geology (Northern Ireland) to six elements of interest; arsenic, chromium, copper, lead, nickel and vanadium. Three methods for identifying domains (areas where a readily identifiable factor can be shown to control the concentration of an element) were used: k-means cluster analysis, boxplots and empirical cumulative distribution functions (ECDF). The ECDF method was most efficient at determining areas of both elevated and reduced concentrations and was used to identify domains in this investigation. Two statistical methods for calculating normal background concentrations (NBCs) and upper limits of geochemical baseline variation (ULBLs), currently used in conjunction with legislative regimes in the UK and Finland respectively, were applied within each domain. The NBC methodology was constructed to run within a specific legislative framework, and its use on this soil geochemical data set was influenced by the presence of skewed distributions and outliers. In contrast, the ULBL methodology was found to calculate more appropriate TTVs that were generally more conservative than the NBCs. TTVs indicate what a “typical” concentration of an element would be within a defined geographical area and should be considered alongside the risk that each of the elements pose in these areas to determine potential risk to receptors.

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.

Correlation analysis as a tool to investigate the bioaccessibility of nickel, vanadium and zinc in Northern Ireland soils

Correlation analyses were conducted on nickel (Ni), vanadium (V) and zinc (Zn) oral bioaccessible fractions (BAFs) and selected geochemistry parameters to identify specific controls exerted over trace element bioaccessibility. BAFs were determined by previous research using the unified BARGE method. Total trace element concentrations and soil geochemical parameters were analysed as part of the Geological Survey of Northern Ireland Tellus Project. Correlation analysis included Ni, V and Zn BAFs against their total concentrations, pH, estimated soil organic carbon (SOC) and a further eight element oxides. BAF data were divided into three separate generic bedrock classifications of basalt, lithic arenite and mudstone prior to analysis, resulting in an increase in average correlation coefficients between BAFs and geochemical parameters. Sulphur trioxide and SOC, spatially correlated with upland peat soils, exhibited significant positive correlations with all BAFs in gastric and gastro-intestinal digestion phases, with such effects being strongest in the lithic arenite bedrock group. Significant negative relationships with bioaccessible Ni, V and Zn and their associated total concentrations were observed for the basalt group. Major element oxides were associated with reduced oral trace element bioaccessibility, with Al2O3 resulting in the highest number of significant negative correlations followed by Fe2O3. spatial mapping showed that metal oxides were present at reduced levels in peat soils. The findings illustrate how specific geology and soil geochemistry exert controls over trace element bioaccessibility, with soil chemical factors having a stronger influence on BAF results than relative geogenic abundance. In general, higher Ni, V and Zn bioaccessibility is expected in peat soil types.

The effects of lead sources on oral bioaccessibility in soil and implications for contaminated land risk management.

Lead (Pb) is a non-threshold toxin capable of inducing toxic effects at any blood level but availability of soil screening criteria for assessing potential health risks is limited. The oral bioaccessibility of Pb in 163 soil samples was attributed to sources through solubility estimation and domain identification. Samples were extracted following the Unified BARGE Method. Urban, mineralisation, peat and granite domains accounted for elevated Pb concentrations compared to rural samples. High Pb solubility explained moderate-high gastric (G) bioaccessible fractions throughout the study area. Higher maximum G concentrations were measured in urban (97.6 mg kg(-1)) and mineralisation (199.8 mg kg(-1)) domains. Higher average G concentrations occurred in mineralisation (36.4 mg kg(-1)) and granite (36.0 mg kg(-1)) domains. Findings suggest diffuse anthropogenic and widespread geogenic contamination could be capable of presenting health risks, having implications for land management decisions in jurisdictions where guidance advises these forms of pollution should not be regarded as contaminated land.

Analysis and numerical modelling of large-scale controls on aquifer structure and hydrogeological properties in the African basement (Benin, West Africa)

Abstract The metamorphic basement units of the Upper Ouémé watershed in Benin have been investigated to identify the structural controls on aquifer properties, groundwater flow and water balance at large scale. Spatial analysis of borehole and hydrogeophysical data suggests that large-scale weathering profiles, aquifer transmissivity and storage properties are better correlated to a palaeo-weathering surface. Multi-model analysis, combined with assessment of nine transient numerical groundwater models against observations, suggests the best conceptualizations are those where hydraulic conductivity and specific yield are distributed within a weathered zone determined through interpolation of weathered zone thickness. When compared to previous studies, the general groundwater balance of simulated models suggests the groundwater system contributes, on average, 49.8 m3 s−1 to the river flow (mostly during the rainy season). The same volumetric flow would be lost to groundwater evapo-transpiration and deep/lateral drainage of the catchment. Borehole abstraction (about 7.5 m3 s−1) represents only 6% of the average groundwater recharge and 1% of the average rainfall. This suggests that despite relatively low borehole productivity, the basement aquifer system still has an important unused potential for rural to mid-scale water supply and that, at present, the main external drivers for groundwater resource sustainability are changes in climate and land use.

Catchment-scale heterogeneity of flow and storage properties in a weathered/fractured hard rock aquifer from resistivity and magnetic resonance surveys: implications for groundwater flow paths and the distribution of residence times

Abstract Groundwater pathways and residence times are controlled by aquifer flow and storage properties, which, in weathered/fractured hard rock aquifers, are characterized by high spatial heterogeneity. Building on earlier work in a metamorphic aquifer in NW Ireland, new clay mineralogy and analyses of geophysical data provided high spatial resolution constraints on the variations in aquifer properties. Groundwater storage values derived from magnetic resonance sounding and electrical resistivity tomography were found to largely vary laterally and with depth, by orders of magnitude. The subsequent implementation of hillslope, two-dimensional numerical groundwater models showed that incorporating heterogeneity from geophysical data in model parametrization led to the best fit to observations compared with a reference model based on borehole data alone. Model simulations further revealed that (1) strong spatial heterogeneity produces deeper, longer groundwater flow paths and higher age mixing, in agreement with the mixed sub-modern/modern ages (mostly <50 years) provided by independent tritium data, and (2) areas with extensive weathering/fracturing are correlated with seepage zones of older groundwater resulting from changes in the flow directions and are likely to act as drainage structures for younger groundwater on a catchment or regional scale. Implications for groundwater resilience to climate extremes and surface pollution are discussed together with recommendations for further research.

The influence of bedrock hydrogeology on catchment-scale nitrate fate and transport in fractured aquifers.

Characterising catchment scale biogeochemical processes controlling nitrate fate in groundwater constitutes a fundamental consideration when applying programmes of measures to reduce risks posed by diffuse agricultural pollutants to water quality. Combining hydrochemical analyses with nitrate isotopic data and physical hydrogeological measurements permitted characterisation of biogeochemical processes influencing nitrogen fate and transport in the groundwater in two fractured bedrock aquifers with contrasting hydrogeology but comparable nutrient loads. Hydrochemical and isotopic analyses of groundwater samples collected from moderately fractured, diffusely karstified limestone indicated nitrification controlled dissolved nitrogen fate and delivery to aquatic receptors. By contrast nitrate concentrations in groundwater were considerably lower in a low transmissivity highly lithified sandstone and pyrite-bearing shale unit with patchy subsoil cover. Geophysical and hydrochemical investigations showed shallower intervals contained hydraulically active fractures where denitrification was reflected through lower nitrogen levels and an isotopic enrichment ratio of 1.7 between δ(15)N and δ(18)O. Study findings highlight the influence of bedrock hydrogeological conditions on aqueous nitrogen mobility. Investigation results demonstrate that bedrock conditions need to be considered when implementing catchment management plans to reduce the impact of agricultural practices on the quality of groundwater and baseflow in receiving rivers. Nitrate isotopic signatures in the groundwater of a freely draining catchment underlain by a karstified aquifer and a poorly draining aquifer with a low transmissivity aquifer.

Identifying Subsurface Features in the Lagan and Enler Valleys in Northern Ireland Using Airborne Tellus Geophysics Data

The City of Belfast in the North of Ireland is underlain by the Triassic Sherwood sandstone aquifer. Increasing abstractions from the aquifer over recent years mean that quantitative assessment of the resource is becoming increasingly important. In addition to this physical pressure, the qualitative status of the groundwater body may be affected by landfilling activity in the urban fringe area surrounding Belfast City. New datasets from the Tellus survey give enough information to allow the identification of subsurface features such as individual dykes, dyke swarms, sills, or combinations of such which may be a factor influencing groundwater flow regimes. A GIS of the compiled data was developed, and using total magnetic intensity mapping and variations on this, it was possible to identify and delineate subsurface anomalies. The igneous intrusions in the study area trend almost perpendicularly to the Lagan Valley and the flow of the River Lagan, and preliminary interpretations have determined that they are subvertical. Previous investigations have shown that flow in this region is mainly from the sides of the Valley towards the course of the River Lagan. To this end, topography driven groundwater flow may be enhanced by preferential dyke-parallel flow.