Characterizing flow pathways in a sandstone aquifer: Tectonic vs sedimentary heterogeneities

Sandstone aquifers are commonly assumed to represent porous media characterized by a permeable matrix. However, such aquifers may be heavy fractured when rock properties and timing of deformation favour brittle failure and crack opening. In many aquifer types, fractures associated with faults, bedding planes and stratabound joints represent preferential pathways for fluids and contaminants. In this paper, well test and outcrop-scale studies reveal how strongly lithified siliciclastic rocks may be entirely dominated by fracture flow at shallow depths (≤ 180 m), similar to limestone and crystalline aquifers. However, sedimentary heterogeneities can primarily control fluid flow where fracture apertures are reduced by overburden pressures or mineral infills at greater depths. The Triassic St Bees Sandstone Formation (UK) of the East Irish Sea Basin represents an optimum example for study of the influence of both sedimentary and tectonic aquifer heterogeneities in a strongly lithified sandstone aquifer-type. This fluvial sedimentary succession accumulated in rapidly subsiding basins, which typically favours preservation of complete depositional cycles including fine grained layers (mudstone and silty sandstone) interbedded in sandstone fluvial channels. Additionally, vertical joints in the St Bees Sandstone Formation form a pervasive stratabound system whereby joints terminate at bedding discontinuities. Additionally, normal faults are present through the succession showing particular development of open-fractures. Here, the shallow aquifer (depth ≤ 180 m) was characterized using hydro-geophysics. Fluid temperature, conductivity and flow-velocity logs record inflows and outflows from normal faults, as well as from pervasive bed-parallel fractures. Quantitative flow logging analyses in boreholes that cut fault planes indicates that zones of fault-related open fractures characterize ~ 50% of water flow. The remaining flow component is dominated by bed-parallel fractures. However, such sub-horizontal fissures become the principal flow conduits in wells that penetrate the exterior parts of fault damage zones, as well as in non-faulted areas. The findings of this study have been compared with those of an earlier investigation of the deeper St Bees Sandstone aquifer (180 to 400 m subsurface depth) undertaken as part of an investigation for a proposed nuclear waste repository. The deeper aquifer is characterized by significantly lower transmissivities. High overburden pressure and the presence of mineral infillings, have reduced the relative impact of tectonic heterogeneities on transmissivity here, thereby allowing matrix flow in the deeper part of the aquifer. The St Bees Sandstone aquifer contrasts the hydraulic behaviour of low-mechanically resistant sandstone rock-types. In fact, the UK Triassic Sandstone of the Cheshire Basin is low-mechanically resistant and flow is supported both by matrix and fracture. Additionally, faults in such weak-rocks are dominated by granulation seams representing flow-barriers which strongly compartmentalize the UK Triassic Sandstone in the Cheshire Basin

Forensic Excavation of Rock Masses: A Technique to Investigate Discontinuity Persistence

True persistence of rock discontinuities (areas with insignificant tensile strength) is an important factor controlling the engineering behaviour of fractured rock masses, but is extremely difficult to quantify using current geological survey methodologies, even where there is good rock exposure. Trace length as measured in the field or using remote measurement devices is actually only broadly indicative of persistence for rock engineering practice and numerical modelling. Visible traces of discontinuities are treated as if they were open fractures within rock mass classifications, despite many such traces being non-persistent and actually retaining considerable strength. The common assumption of 100% persistence, based on trace length, is generally extremely conservative in terms of strength and stiffness, but not always so and may lead to a wrong prediction of failure mechanism or of excavatability. Assuming full persistence would give hopelessly incorrect predictions of hydraulic conductivity. A new technique termed forensic excavation of rock masses is introduced, as a procedure for directly investigating discontinuity persistence. This technique involves non-explosive excavation of rock masses by injecting an expansive chemical splitter along incipient discontinuities. On expansion, the splitter causes the incipient traces to open as true joints. Experiments are described in which near-planar rock discontinuities, through siltstone and sandstone, were opened up by injecting the splitter into holes drilled along the lines of visible traces of the discontinuities in the laboratory and in the field. Once exposed the surfaces were examined to investigate the pre-existing persistence characteristics of the incipient discontinuities. One conclusion from this study is that visible trace length of a discontinuity can be a poor indicator of true persistence (defined for a fracture area with negligible tensile strength). An observation from this series of experiments was that freshly failed surfaces through pre-existing rock bridges were relatively rough compared to sections of pre-existing weaker areas of geologically developed (though still incipient) discontinuities. Fractographic features such as hackle and rib marks were typical of the freshly broken rock bridges, whereas opened-up areas of incipient discontinuity were smoother. Schmidt hammer rebound values were generally higher for the rock bridge areas, probably reflecting their lower degree of chemical and physical weathering

The geometry of fluvial channel bodies: Empirical characterization and implications for object-based models of the subsurface

The distribution of channel deposits in fluvial reservoirs is commonly modeled with object-based techniques, constrained on quantities describing the geometries of channel bodies. To ensure plausible simulations, it is common to define inputs to these models by referring to geologic analogs. Given their ability to reproduce complex geometries and to draw upon the analog experience, object-based models are considered inherently realistic. Yet this perceived realism has not hitherto been tested by assessing the outputs of these techniques against sedimentary architectures in the stratigraphic record. This work presents a synthesis of data on the geometry of channel bodies, derived from a sedimentologic database, with the following aims: (1) to provide tools for constraining stochastic models of fluvial reservoirs in data-poor situations, and (2) to test the intrinsic realism of object-based modeling algorithms by comparing characteristics of the modeled architectures against analogs. An empirical characterization of the geometry of fluvial channel bodies is undertaken that describes distributions in (and relationships among) channel-body thickness, cross-stream width, and planform wavelength and amplitude. Object-based models are then built running simulations conditioned on six alternative, analog-informed parameter sets, using four algorithms according to nine different approaches. Closeness of match between analogs and models is then determined on a statistical basis. Results indicate which modeling approaches return architectures that more closely resemble the organization of fluvial depositional systems known from nature and in what respect. None of the tested algorithms fully reproduce characteristics seen in natural systems, demonstrating the need for subsurface modeling methods to better incorporate geologic knowledge

Laboratory study of spectral induced polarization responses of magnetite - Fe2+ redox reactions in porous media

Spectral Induced Polarization (SIP) phase anomalies in field surveys at contaminated sites have previously been shown to correlate with the occurrence of chemically reducing conditions and/or semiconductive minerals, but the reasons for this are not fully understood. We report a systematic laboratory investigation of the role of the semiconductive mineral magnetite and its interaction with redoxactive versus redox-inactive ions in producing such phase anomalies. The SIP responses of quartz sand with 5% magnetite in solutions containing redox-inactive Ca2+, and Ni2+, versus redox-active Fe2+ were measured across the pH ranges corresponding to adsorption of these metals to magnetite. With redox inactive ions Ca2+ and Ni2+, SIP phase response showed no changes across the pH range 4 to10, corresponding to their adsorption, showing ~30 mrad anomalies peaking at ~59 to 74 Hz. These large phase anomalies are probably caused by polarization of the magnetite-solution interfaces. With the redox-active ion Fe2+, frequency of peak phase res onse decreased progressively from ~46 to ~3 Hz as effluent pH increased from 4 to 7, corresponding to progressive adsorption of Fe2+ to the magnetite surface. The latter frequency (3 Hz) corresponds approximately with those of phase anomalies detected in field surveys reported elsewhere. We conclude that pH sensitivity arises from redox reactions between Fe2+ and magnetite surfaces, with transfer of electrical charge through the bulk mineral, as reported in other laboratory investigations. Our results confirm that SIP measurements are sensitive to redox reactions involving charge transfers between adsorbed ions and semiconductive minerals. Phase anomalies seen in field surveys of groundwater contamination and biostimulation may therefore be indicative of iron-reducing conditions, when semiconductive iron minerals such as magnetite are present

Characterisation of subglacial water using a constrained transdimensional Bayesian transient electromagnetic inversion

Subglacial water modulates glacier-bed friction and therefore is of fundamental importance when characterising the dynamics of ice masses. The state of subglacial pore water, whether liquid or frozen, is associated with differences in electrical resistivity that span several orders of magnitude; hence, liquid water can be inferred from electrical resistivity depth profiles. Such profiles can be obtained from inversions of transient (time-domain) electromagnetic (TEM) soundings, but these are often non-unique. Here, we adapt an existing Bayesian transdimensional algorithm (Multimodal Layered Transdimensional Inversion – MuLTI) to the inversion of TEM data using independent depth constraints to provide statistical properties and uncertainty analysis of the resistivity profile with depth. The method was applied to ground-based TEM data acquired on the terminus of the Norwegian glacier, Midtdalsbreen, with depth constraints provided by co-located ground-penetrating radar data. Our inversion shows that the glacier bed is directly underlain by material of resistivity 102 Ωm ± 1000 %, with thickness 5–40 m, in turn underlain by a highly conductive basement (100 Ωm ± 15 %). High-resistivity material, 5×104 Ωm ± 25 %, exists at the front of the glacier. All uncertainties are defined by the interquartile range of the posterior resistivity distribution. Combining these resistivity profiles with those from co-located seismic shear-wave velocity inversions to further reduce ambiguity in the hydrogeological interpretation of the subsurface, we propose a new 3-D interpretation in which the Midtdalsbreen subglacial material is partitioned into partially frozen sediment, frozen sediment/permafrost and weathered/fractured bedrock with saline water

Numerical simulation of spring hydrograph recession curves for evaluating behavior of the East Yorkshire chalk aquifer

The Cretaceous Chalk aquifer is the most important in the UK for the provision of water to public supply and agriculture. The Chalk has both matrix and fracture porosity and is thus best considered as a dual porosity aquifer system. Although the matrix porosity is large, typically around 0.35 in the study area of East Yorkshire, UK (ESI, 2010), pore diameters are typically very small, and the water contained in them is virtually immobile. The high permeability fracture network is responsible for the ability of water to drain; spatial variations in fracture network properties mean conventional approaches to aquifer characterization such as borehole pumping tests are of limited utility. Hence this study attempts to better understand the flow system and characterise aquifer properties from the recession response seen at springs during the spring/summer period when recharge is minimal. This approach has the advantage that spring hydrographs represent the sum of the response from entire catchments. This paper reports numerical modeling for simulating aquifer and spring responses during hydrological recession. Firstly, available geological and hydrogeological information for the study area was used to develop hydrogeological conceptual models. Four different numerical models have been constructed representing four possible scenarios that could represent the aquifer in the selected area. These are: single reservoir aquifer, double reservoir aquifer, single reservoir aquifer with highly permeable vertical zone intersecting the spring location and single reservoir aquifer containing tunnel shaped highly permeable zone at the spring elevation respectively. The sensitivity of spring recession response to various external and internal parameter values was investigated, to understand relations between spring recession, hydrological inputs (recharge) and aquifer structure. Spring hydrographs from the real aquifer were compared with the hydrographs generated from models, in order to estimate aquifer properties. The work aims to identify the utility of spring hydrographs in eliciting aquifer permeability structure, as well as identifying the conceptual scenario which best represents the Chalk Aquifer in East Yorkshire, UK

Groundwater flow velocities in a fractured carbonate aquifer-type: Implications for contaminant transport

Contaminants that are highly soluble in groundwater are rapidly transported via fractures in mechanically resistant sedimentary rock aquifers. Hence, a rigorous methodology is needed to estimate groundwater flow velocities in such fractured aquifers. Here, we propose an approach using borehole hydraulic testing to compute flow velocities in an un-faulted area of a fractured carbonate aquifer by applying the cubic law to a parallel plate model. The Cadeby Formation (Yorkshire, NE England) – a Permian dolostone aquifer present beneath the University of Leeds Farm - is the fractured aquifer selected for this hydraulic experiment. The bedding plane fractures of this dolostone aquifer, which are sub-horizontal, sub-parallel and laterally persistent, largely dominate the flow at shallow (<~40 mBGL) depths. These flowing bedding plane discontinuities are separated by a rock matrix which is relatively impermeable (Kwell-test/Kcore-plug~104) as is common in fractured carbonate aquifers. In the workflow reported here, the number of flowing fractures - mainly bedding plane fractures - intersecting three open monitoring wells are found from temperature/fluid conductivity and acoustic/optical televiewer logging. Following well installation, average fracture hydraulic apertures for screened intervals are found from analysis of slug tests. For the case study aquifer, this workflow predicts hydraulic apertures ranging from 0.10 up to 0.54 mm. However, groundwater flow velocities range within two order of magnitude from 13 up to 242 m/day. Notably, fracture apertures and flow velocities rapidly reduce with increasing depth below the water table; the upper ~10 m shows relatively high values of hydraulic conductivity (0.30–2.85 m/day) and corresponding flow velocity (33–242 m/day). Permeability development around the water table in carbonate aquifer-types is common, and arises where high pCO2 recharge water from the soil zone causes calcite/dolomite dissolution. Hence, agricultural contaminants entering the aquifer with recharge water are laterally transported rapidly within this upper part. Computation of groundwater flow velocities allows determination of the Reynolds number. Values of up ~1, indicating the lower limit of the transition from laminar to turbulent flow, are found at the studied site, which is situated away from major fault traces. Hence, turbulent flow is likely to arise in proximity to tectonic structures, such as normal faults, which localize flow and enhance karstification. The occurrence of turbulent flow in correspondence of such tectonic structures should be represented in regional groundwater flow simulations

Tensile strength of incipient rock discontinuities

This paper addresses the issue of the tensile strength of incipient discontinuities in rock and presents preliminary results from a series of laboratory studies. In most rock masses rock discontinuities, as veins or incipient fractures, often retain some tensile strength that may approach that of the parent rock. This fact is of high importance to rock mass strength but is generally ignored, neglected or underestimated. Samples of incipient rock discontinuities including joints, bedding and mineral veins have been tested in direct tension. It has been confirmed that incipient, visible and discrete discontinuities, that might be recorded as ‘joints’ in a rock mass characterisation programme can indeed have high tensile strength, approaching that of parent rock. Others are of course far weaker. The factors contributing to tensile strength have been examined. It is concluded that the degree of incipiency of rock discontinuities needs to be differentiated in the process of rock mass classification and engineering design and this can best be done with reference to the tensile strength relative to that of the parent rock

Socioeconomic variation in the incidence of childhood coeliac disease in the UK.

BACKGROUND: Serological studies indicate that evidence of coeliac disease (CD) exists in about 1% of all children, but we lack estimates of current diagnostic patterns among children and how they vary by socioeconomic group. METHODS: We identified all children aged 0-18 years between 1993 and 2012 who were registered with general practices across the UK that contribute to a large population-based general practice database. The incidence of CD was evaluated in each quintile of the Townsend index of deprivation and stratified by age, sex, country and calendar year. RESULTS: Among 2,063,421 children, we identified 1247 CD diagnoses, corresponding to an overall CD incidence of 11.9 per 100,000 person-years, which was similar across the UK countries and higher in girls than in boys. We found a gradient of CD diagnosis across socioeconomic groups, with the rate of diagnosis being 80% higher in children from the least-deprived areas than in those from the most-deprived areas (incident rate ratio 1.80, 95% CI 1.45 to 2.22). This pattern held for both boys and girls and across all ages. Across all four countries of the UK, we found similar associations between CD and socioeconomic status. While CD incidence up to age 2 remained stable over the study period, diagnoses at older ages have almost tripled over the past 20 years. CONCLUSIONS: Children living in less socioeconomically deprived areas in the UK are more likely to be diagnosed with CD. Increased implementation of diagnostic guidelines could result in better case identification in more-deprived areas

A molecular insight into algal-oomycete warfare : cDNA analysis of Ectocarpus siliculosus infected with the basal oomycete Eurychasma dicksonii

Peer reviewedPublisher PD