Heterogeneous aquitard properties in sedimentary successions in the Apennine chain: case studies in southern Italy

The aim of this study was to analyse the hydrogeologic role of siliciclastic deposits and marly‐clayey‐calcareous successions within the carbonate Apennine chain (southern Italy). The study was carried out along the northern part of the Matese carbonate massif through (1) the hydraulic characterization of siliciclastic rocks in a test site, by means of Lugeon tests, and (2) the identification of the groundwater flow system discharging at an important spring located within a marly‐clayey‐calcareous succession in a second test site, by means of isotopic investigations. The results showed that the investigated siliciclastic deposits and marly‐clayey‐calcareous successions may allow significant groundwater discharge from carbonate aquifers. Thus, they do not everywhere behave as aquitard, contrary to the previous model. Instead, groundwater flows through the upper part of these successions, where stress release fracturing enhanced rock permeability in the near‐surface bedrock. Thus, these successions may locally be a new groundwater source within the southern Apennine chain. Copyright © 2009 John Wiley & Sons, Ltd.     

The effect of low‐permeability fault zones on groundwater flow in a compartmentalized system. Experimental evidence from a carbonate aquifer (Southern Italy)

The hydrogeological behaviour of fault zones in carbonate aquifers is often neglected in conceptual and numerical models. Furthermore, no information is available regarding the relationships between piezometric levels when significant compartmentalization occurs due to the occurrence of low‐flow fault zones. The aim of this study was to refine the conceptualization of subsurface flow in faulted carbonate aquifers and to analyse relationships between sub‐basins within a compartmentalized aquifer system in Southern Italy. The interactions between compartments that straddle low‐flow faults were investigated over four hydrologic years using a statistical approach to compare (i) the hydraulic heads within two wells located up‐ and down‐gradient of tectonic discontinuities as well as (ii) the rainfall and piezometric levels. The results of this study suggest that a set of barriers exists between the wells, and, therefore, the total head loss observed between the wells (approximately 80 m) should be distributed across several aquitards, with one aquitard exhibiting a relatively high permeability or low degree of integrity. Due to slight differences in permeability, transient conditions in aquitards can occur over relatively short periods, which is in agreement with the results of the statistical data analysis. Consequently, rather than being caused by pure aquitards, aquifer system compartmentalization likely results from slight differences in the permeability between lower‐permeability fault zones and adjacent higher‐permeability protoliths. Copyright © 2014 John Wiley & Sons, Ltd.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of beach slope breaks on nearshore groundwater dynamics

Interactions between fresh groundwater and seawater affect significantly the nearshore pore water flow, which in turn influences the fate of nutrients and contaminants in coastal aquifers prior to discharge to the marine environment. Field investigations and numerical simulations were carried out to examine the groundwater dynamics in the intertidal zone of a carbonate sandy aquifer on the tropical island of Rarotonga, Cook Islands. The study site was featured by distinct cross‐shore slope breaks on the beach surface. Measured pore water salinities revealed different distributions under the influences of different beach profiles, inland heads, and tidal oscillations. Fresh groundwater was found to discharge around a beach slope break located in the middle area of the intertidal zone. The results indicate a strong interplay between the slope break beach morphology and tidal force in controlling the nearshore groundwater flow and solute transport. The fresh groundwater discharge location was largely determined by the beach morphology in combination with the tidal force. The nearshore groundwater flow can be very sensitive to beach slope breaks, which induce local circulation and flow instabilities. As slope breaks are a common feature of beaches around the world, these results have important, general implications for future studies of nutrients transport and transformations in nearshore aquifers and associated fluxes via submarine groundwater discharge.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater Availability as Constrained by Hydrogeology and Environmental Flows

Groundwater pumping from aquifers in hydraulic connection with nearby streams has the potential to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes‐St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations may constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin may be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time, regional and local hydrogeology, streambed conductance, and streamflow depletion limits. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.     

Closed Analytic Elements with Flexible Geometry

A simple and fast treatment of hydrogeologic features with irregularly shaped boundaries in two‐dimensional analytic element groundwater flow models is presented. The star domain shapes of the features are restricted to closed shapes represented as smooth and continuous single‐valued functions of distance from a focus point, . The element can be used to treat a variety of boundary and continuity conditions, including those of irregularly shaped lakes or heterogeneities in hydraulic conductivity. The new element is demonstrated via some simple illustrative test cases and shown to be efficient, accurate, and much simpler to implement than existing solutions for irregular shapes.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Numerical modeling of geothermal groundwater flow in karst aquifer system in eastern Weibei, Shaanxi Province, China

The quantitative assessment of geothermal water resources is important to the exploitation and utilization of geothermal resources. In the geothermal water systems the density of groundwater changes with the temperature, therefore the variations in hydraulic heads and temperatures are very complicated. A three-dimensional density-dependent model coupling the groundwater flow and heat transport is established and used to simulate the geothermal water flow in the karst aquifers in eastern Weibei, Shaanxi Province, China. The multilayered karst aquifer system in the study area is cut by some major faults which control the regional groundwater flow. In order to calibrate and simulate the effect of the major faults, each fault is discretized as a belt of elements with special hydrological parameters in the numerical model. The groundwater dating data are used to be integrated with the groundwater flow pattern and calibrate the model. Simulation results show that the calculated hydraulic heads and temperature fit with the observed data well.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Characterizing Groundwater Flow in Monitoring Wells by Altering Dissolved Oxygen

In this study, the dissolved oxygen (DO) alteration method (Chlebica and Robbins 2013) is used to evaluate the patterns of flow into and vertically within shallow screened monitoring wells. The method entails bubbling air into a well, followed by conducting DO profiles with time. Tests were conducted in six standard 2" (5 cm) polyvinyl chloride shallow screened monitoring wells at four test sites in Storrs, Connecticut. Test sites vary in formation permeability, flow patterns, and nearby geographic features influencing flow. The method provides a means for groundwater flow characterization and potentially improved interpretation of contamination sampling results in the absence of detailed three‐dimensional hydrogeologic information.     

Static and dynamic conceptual model of a complexly fractured crystalline rock aquifer

A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies.     

A proposed structure–groundwater model: application to chalky media of the Paris Basin

The groundwater flow in a fissured chalky environment at the northern border of the Paris Basin depends on several geological and hydrogeological parameters. Although the studied sector of the basin presents a homogeneous rock type, it is affected by a fracture network. In this type of environment, in which the permeability is low, the groundwater flow displays significant disruption, which is localized in the Fruges region (northern France). The interconnection of the discontinuities (network of fault and/or joints) is reliant on the structural control of groundwater flow through increases in the hydraulic connection between the unsaturated and the saturated zone. The methodology developed herein makes use of microstructural and regional analysis of the fracture patterns, and allowed consideration of the piezometric variations of the chalk aquifer during periods of low and high groundwater levels (April and October 2001) and a diagraphic representation of the estimated physical parameters (electrical resistivity). This enabled us to construct a ‘flow structure’ conceptual model in which we identify two types of faults: tight walls and flow paths that control the piezometric heads and the flow rate. Model validation was carried out on a similar sector. Copyright © 2007 John Wiley & Sons, Ltd.     

Induced Temperature Gradients to Examine Groundwater Flowpaths in Open Boreholes

Techniques for characterizing the hydraulic properties and groundwater flow processes of aquifers are essential to design hydrogeologic conceptual models. In this study, rapid time series temperature profiles within open‐groundwater wells in fractured rock were measured using fiber optic distributed temperature sensing (FO‐DTS). To identify zones of active groundwater flow, two continuous electrical heating cables were installed alongside a FO‐DTS cable to heat the column of water within the well and to create a temperature difference between the ambient temperature of the groundwater in the aquifer and that within the well. Additional tests were performed to examine the effects of pumping on hydraulic fracture interconnectivity around the well and to identify zones of increased groundwater flow. High‐ and low‐resolution FO‐DTS cable configurations were examined to test the sensitivities of the technique and compared with downhole video footage and geophysical logging to confirm the zones of active groundwater flow. Two examples are presented to demonstrate the usefulness of this new technique for rapid characterization of fracture zones in open boreholes. The combination of the FO‐DTS and heating cable has excellent scope as a rapid appraisal tool for borehole construction design and improving hydrogeologic conceptual models.     

Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.     

Fine‐scale characteristics of groundwater flow in a peatland

Fine‐scale dynamics of groundwater flow were studied in a 1·5 ha peatland in central New York. Measurements of the hydraulic head throughout a detailed network of piezometer clusters revealed spatial and temporal variability in the direction of groundwater flow at a very fine (within a few metres) scale of analysis. Within the small wetland, there were areas of groundwater recharge, discharge and lateral flow. Such patterns of groundwater flow frequently reversed or changed due to fluctuations of only a few centimetres in hydraulic head. Specific conductance, deuterium signatures and calcium concentrations of groundwater corroborated the groundwater flow patterns determined with hydraulic head measurements and illustrated the influence of source water chemistry and evaporation on different layers in the peat column. The control of peat chemistry by such fine‐scale groundwater flow may have important implications for plant community composition and diversity in groundwater‐fed peatlands. Copyright © 1999 John Wiley & Sons, Ltd.     

Aquifer–peatland connectivity in southern Quebec (Canada)

In areas where peatlands are abundant, they are likely to play a significant role in the hydrological and hydrogeological dynamics of a watershed. Although individual case studies are reported in the literature, there is a large range of aquifer–peatland interactions and there is a need to understand the controls of these interactions. The objectives of this study were (1) to better understand aquifer–peatland connections and how these may be predicted by geology and geomorphic location and (2) to provide a variety of reference sites for glacial geological settings. Slope and depression peatlands were studied in the Abitibi‐Témiscamingue region and in the St. Lawrence Lowlands, two contrasting regions of southern Quebec. A total of 12 transects that span a shallow aquifer–peatland interface were instrumented with piezometers. Field investigations included peatland characterization, monthly water level monitoring, and continuous hydraulic head measurements with pressure transducers. The results indicate that 7 of the 12 transects receive groundwater from the surrounding shallow aquifer. At the peatland margin, four lateral flow patterns were identified and associated with slope peatlands (parallel inflow and divergent flow) and with depression peatlands (convergent flow and parallel outflow). Vertical hydraulic gradients suggest that water flows mainly downwards, i.e. from the peatland to the underlying mineral deposits. Vertical connectivity appears to decrease as the distance from the peatland margin increases. All of these exchanges are important components in the sustainability of peatland hydrogeological functions. The regional comparison of aquifer–peatland flow dynamics performed in this study provides a new set of referenced data for the assessment of aquifer–peatland connectivity. Copyright © 2014 John Wiley & Sons, Ltd.     

Meltwater discharge through the subglacial bed and its land‐forming consequences from numerical experiments in the Polish lowland during the last glaciation

Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd.     

A multi‐method field experiment to determine local groundwater flow in a glacier forefield

We implemented multiple independent field techniques to determine the direction and velocity of groundwater flow at a specific stream reach in a glacier forefield. Time‐lapse experiments were conducted using two electrical resistivity tomography (ERT) lines installed in a cross pattern. A circular array of groundwater tubes was also installed to monitor groundwater flow via discrete salt injections. Both inter‐borehole and ERT results confirmed this stream section as a losing reach and enabled quantification of the flow direction. Both techniques yielded advection velocities varying between 5.7 and 21.8 m/day. Estimates of groundwater flow direction and velocity indicated that groundwater infiltrates from the stream nearby and not from the adjacent lateral moraine. Groundwater age estimated from radon concentration measurements supported this hypothesis. Despite uncertainties inherent to each of the methods deployed, the combination of multiple field techniques allowed drawing consistent conclusions about local groundwater flow. We thus regard our multi‐method approach as a reliable way to characterize the two‐dimensional groundwater flow at sites where more invasive groundwater investigation techniques are difficult to carry out and local heterogeneities can make single measurements unreliable. Copyright © 2014 John Wiley & Sons, Ltd.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.