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 (≤180m), 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≤180m) 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 indicate 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 400m 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.

Geothermal reservoirs are highly anisotropic and heterogeneous, and thus require a variety of structural geology, geomechanical, remote sensing, geophysical and hydraulic techniques to inform Discrete Fracture Network flow models. Following the Paris Agreement on reduction of carbon emissions, such reservoirs have received more attention and new techniques that support Discrete Fracture Network models were developed. A comprehensive review is therefore needed to merge innovative and traditional technical approaches into a coherent framework to enhance the extraction of geothermal energy from the deep subsurface. Traditionally, statistics extracted from structural scanlines and unmanned aerial vehicle surveys on analogues represent optimum ways to constrain the length of joints, bedding planes, and faults, thereby generating a model of the network of fractures. Combining borehole images with seismic attributes has also proven to be an excellent approach that supports the stochastic generation of Discrete Fracture Network models by detecting the orientation, density, and dominant trends of the fractures in the reservoirs. However, to move forward to flow modelling, computation of transmissivities from pumping tests, and the determination of hydraulically active fractures allow the computation of the hydraulic aperture in permeable sedimentary rocks. The latter parameter is fundamental to simulating flow in a network of discrete fractures. The mechanical aperture can also be estimated based on the characterization of geomechanical parameters (Poisson’s ratio, and Young’s modulus) in Hot Dry Rocks of igneous-metamorphic origin. Compared with previous review studies, this paper will be the first to describe all the geological and hydro-geophysical techniques that inform Discrete Fracture Network development in geothermal frameworks. We therefore envisage that this paper represents a useful and holistic guide for future projects on preparing DFN models.

Fractured aquifers are heterogeneous due to the variable frequency, orientation, and intersections of rock discontinuities. A ~100-m-thick Silurian dolostone sequence provides a bedrock aquifer supplying the city of Guelph, Canada. Here, fracture network characteristics and associated influences on hydraulic head were examined using several data types obtained from 24 cored holes in a study that is novel for the quantity and quality of data. High (50–90°) angle joint orientations, heights, and terminations relative to bedding features were determined from acoustic televiewer logs and outcrop scanlines. These data were compared to high-resolution hydraulic head profiles showing head loss over depth-discrete intervals identifying zones with lower vertical hydraulic conductivity. This study reveals that the marl-rich Vinemount Member, traditionally considered the principal aquitard, corresponds to head loss in only 62% of the 24 boreholes. The vertical position of head loss varies across the 90-km² study area and occurs in any of the lithostratigraphic units of the Lockport Group. Within this sedimentary sequence, aquitards are laterally discontinuous or “patchy” at variable depths and relate to: (1) the frequency of the high-angle joints; (2) shorter joint height; and (3) the type of joint terminations. The head loss occurs in thin (2–2.5 m) intervals where the frequency of the high-angle joints is low. Where a large proportion of small joints cross-cut marl bedding planes, head loss is negligible, suggesting that the vertical hydraulic conductivity is not reduced. Overall, these findings are potentially applicable to assessing aquitard and cap rock integrity in carbonate sedimentary sequences worldwide.

Sandstones are characterized by different hydraulic behaviours and need to be modelled in various ways to represent groundwater flow and contaminant transport. This review shows how sandstone aquifers within the UK Triassic Sherwood Sandstone Group can be represented using three modelling approaches: the Conduit Network, Discrete Fracture Network and Equivalent Porous Medium. The Sherwood Sandstone aquifer is dominated by matrix flow in the Eastern England Shelf, Worcester, Needwood and Staffordshire basins. Here, the aquifers are modelled as Equivalent Porous Media at different spatial scales. Fractures represent the principal flow pathways in the Cheshire Basin. In this basin, Discrete Fracture Network models that account for diffusivity in the matrix can be used where the domain scale is small. The Sherwood Sandstone aquifer across northwestern England shows evidence of intense groundwater alteration and high flow velocities in solutionally enlarged fractures. Turbulent flowing pipe-elements can be inserted in the modelling domain represented by Equivalent Porous Medium at specific sites. The review shows how the Sherwood Sandstone aquifer as well as other siliciclastic deposits across the world need to be represented using a range of modelling approaches, as they behave as matrix or fracture flow aquifers, or in specific cases show a karst-like behaviour. Thematic collection: This article is part of the Hydrogeology of Sandstone collection available at: https://www.lyellcollection.org/cc/hydrogeology-of-sandstone

Fluvio-aeolian sedimentary successions host groundwater aquifers at shallow depths (<~0.15 km), which overlie geothermal and shale-gas reservoirs, and nuclear waste repositories at intermediate depths (~0.15–2.0 km). Additionally, such deposits represent petroleum reservoirs at greater depths (~2.0–4.0 km). The need to improve conceptual understanding of the hydraulic behaviour of fluvial-aeolian sandstone successions over a large depth interval (~0–4 km) is important for socio-economic reasons. Thus, the hydraulic properties of the Triassic Sherwood Sandstone aquifer in the UK have been reviewed and compared to similar fluvio-aeolian successions. The ratio between well-scale and core-plug-scale permeability (K_well-test/K_core-plug) acts as a proxy for the relative importance of fracture versus intergranular flow. This ratio (which typically varies from ~2 to 100) indicates significant contribution of fractures to flow at relatively shallow depths (<~0.15 km). Here, permeability development is controlled by dissolution of calcite-dolomite in correspondence of fractures. The observed ratio (K_well-test/K_core-plug) decreases with depth, approaching unity, indicating that intergranular flow dominates at ~1 km depth. At depths ≥ ~1 km, dissolution of carbonate cement by rock alteration due to groundwater flow is absent and fractures are closed. Aeolian and fluvial deposits behave differently in proximity to normal faults in the Sherwood Sandstone aquifer. Deformation bands in aeolian dune deposits strongly compartmentalize this aquifer. The hydro-structural properties of fluvio-aeolian deposits are also controlled by mineralogy in fault zones. A relative abundance of quartz vs. feldspar and clays in aeolian sandstones favours development of low-permeability deformation bands.