Simulation of groundwater flow in a crystalline rock aquifer system in Southern Ghana – An evaluation of the effects of increased groundwater abstraction on the aquifers using a transient groundwater flow model

Monitored groundwater level data, well logs, and aquifer data as well as the relevant surface hydrological data were used to conceptualise the hydrogeological system of the Densu Basin in Southern Ghana. The objective was to numerically derive the hydraulic conductivity field for better characterization of the aquifer system and for simulating the effects of increasing groundwater abstraction on the aquifer system in the basin. The hydraulic conductivity field has been generated in this study through model calibration. This study finds that hydraulic conductivity ranges between a low of 2 m/d in the middle sections of the basin and about 40 m/d in the south. Clear differences in the underlying geology have been indicated in the distribution of aquifer hydraulic conductivities. This is in consonance with the general assertion that the hydrogeological properties of the aquifers in the crystalline basement terrains are controlled by the degree of fracturing and/or weathering of the country rock. The transient model suggest aquifer specific storage values to range between 6.0 × 10^?5 m^?1 and 2.1 × 10^?4 m^?1 which are within acceptable range of values normally quoted for similar lithologies in the literature. There is an apparent subtle decrease in groundwater recharge from about 13% of the annual precipitation in 2005 to about 10.3% of the precipitation in 2008. The transient model was used to simulate responses of the system to annual increment of groundwater abstraction by 20% at the 2008 recharge rates for the period 2009 – 2024. The results suggest that the system will not be able to sustain this level of abstraction as it would lead to a basin wide drawdown in the hydraulic head by 4 m by the end of the prediction period. It further suggests a safe annual increment in groundwater abstraction by 5% under business as usual recharge conditions. Identification and protection of groundwater recharge areas in the basin are recommended in order to safeguard the integrity of the resource under the scenario of increased abstraction for commercial activities in the basin. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of hydrogeological setting in two Canadian peatlands investigated through 2D steady-state groundwater flow modelling

This study investigated how hydrogeological setting influences aquifer–peatland connections in slope and basin peatlands. Steady-state groundwater flow was simulated using Modflow on 2D transects for an esker slope peatland and for a basin peatland in southern Quebec (Canada). Simulations investigated how hydraulic heads and groundwater flow exported toward runoff from the peatland can be influenced by recharge, hydraulic properties, and heterogeneity. The slope peatland model was strongly dominated by horizontal flow from the esker. This suggests that slope peatlands are dependent on the hydrogeological conditions of the adjacent aquifer reservoir, but are resilient to hydrological changes. The basin peatland produced groundwater outflow to the surface aquifer. Lateral and vertical peat heterogeneity due to peat decomposition or compaction were identified as having a significant influence on fluxes. These results suggest that basin peatlands are more dependent on recharge conditions, and could be more susceptible to land use and climate changes.     

Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems

The hydrological influence of fault zones in tectonic areas is usually difficult to depict from field data. Numerical simulation allows representation of such flow systems and an estimation of flow lines and rates. This paper reports on simulations of the groundwater flow in a range‐and‐basin area affected by a regional fault zone, which may drain or recharge an overlaying alluvial aquifer. Different hydraulic conductivity values for the range rocks, the fault‐zone, and the sedimentary infill of the basin are considered, as well as different fault‐zone widths and boundary conditions. Results show that upward and downward fluxes develop in the upper part of the fault zone, controlled by the action of the alluvial aquifer, influencing the recharge of the sedimentary basin. This paper shows the hydrological efficiency of fault zones as preferential flow; it also analyses the constraints that determine groundwater recharge to the surrounding basins. These results contribute to the understanding of hydrogeological dynamics in tectonic areas. Copyright © 2008 John Wiley & Sons, Ltd.     

Groundwater discharge mechanism in semi‐arid regions and the role of evapotranspiration

The present study examined groundwater recharge/discharge mechanisms in the regional Central Sudan Rift Basins (CSRB). Aquifers in CSRB constitute poorly sorted silisiclastics of sand, clay and gravels deposited in closed hydrologic systems of the Cretaceous–Pleistocene fluviolacustrine environments. CSRB are bounded to the north by the highlands of the Central African Shear Zone (CAZS) that represents the surface and groundwater divides. Sporadic recharge in the peripheries of the basins along the CASZ occurs subsequent to decadal and centennial storm events. Inflow from the Nile into the aquifers represents an additional source of recharge. Thus, groundwater resources cannot be labelled fossil nor can they be readily recharged. Closed hydrologic troughs located adjacent to the influent Nile system mark areas of main groundwater discharge characterized by lower hydraulic heads. This study has examined mechanisms that derive the discharge of the groundwater in these closed basins and concluded that only evapotranspirative discharge can provide a plausible explanation. Groundwater abstraction is mainly through deep‐rooted trees and effective evaporation. The increase of TDS along the flow indicates local recharge at the peripheries of basins and shows the influence of evaporation and rock/water interaction. The decline in groundwater level along a flow path was calculated using Darcy's law to estimate average recharge and evapotranspirative discharge, which are equal under natural equilibrium and make the only fluxes in CSRB. Steady‐state 2D flow modelling has demonstrated that an average recharge of 4–8 mm yr^?1 and evapotranspirative discharge of 1–22 mm yr^?1 will maintain natural equilibrium in CSRB. Sporadic storms provide recharge in the highlands to preserve the current hydraulic gradient and maintain aquifer dynamics. Simulated recharge from the Nile totals about 17·5 mm yr^?1 and is therefore a significant contributor to the water balance. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of Groundwater Quality Trends Resulting from Anthropogenic Changes in Southeast Florida

The effects of surface water flow system changes caused by constructing water‐conservation areas and canals in southeast Florida on groundwater quality under the Atlantic Coastal Ridge was investigated with numerical modeling. Water quality data were used to delineate a zone of groundwater with low total dissolved solids (TDS) within the Biscayne aquifer under the ridge. The delineated zone has the following characteristics. Its location generally coincides with an area where the Biscayne aquifer has high transmissivities, corresponds to a high recharge area of the ridge, and underlies a part of the groundwater mound formed under the ridge prior to completion of the canals. This low TDS groundwater appears to be the result of pre‐development conditions rather than seepage from the canals constructed after the 1950s. Numerical simulation results indicate that the time for low TDS groundwater under the ridge to reach equilibrium with high TDS surface water in the water‐conservation areas and Everglades National Park are approximately 70 and 60 years, respectively. The high TDS groundwater would be restricted to the water‐conservation areas and the park due to its slow eastward movement caused by small hydraulic gradients in Rocky Glades and its mixing with the low TDS groundwater under the high‐recharge area of the ridge. The flow or physical boundary conditions such as high recharge rates or low hydraulic conductivity layers may affect how the spatial distribution of groundwater quality in an aquifer will change when a groundwater flow system reaches equilibrium with an associated surface water flow system.     

Recharge of aquifers by flood events in an arid region

The process of aquifer recharge by flood events in an arid region was investigated as applied to the Hazeva Formation (the Karkom graben, the Wadi Paran watershed, Israel). The hydrological model was established as a complex system, with due regard for groundwater and transmission losses of surface runoff. It was based on a previously outlined hydrogeological model of the Karkom graben and a model of transmission losses in arid watercourses under conditions of data deficiency. Proceeding from calculation of groundwater balances, the contribution of surface runoff as a decisive balance component was confirmed. The main characteristics of aquifer regimes, such as changes in storage volume and groundwater level, as well as lateral flow, were all found to be dependent upon the net extraction rate, i.e. pumpage discounting replenishment by flood events. Analysis and physical interpretation of model parameters enabled assessment of the influence of groundwater extraction on aquifer recharge. This became apparent as increasing absorption capacity and recharge availability of the aquifer as a result of the groundwater abstraction. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical modelling of the impacts of water abstraction for hydraulic fracturing on groundwater–surface water interaction: a case study from northwestern Alberta,Canada

ABSTRACT

The temporal dynamics of groundwater–surface water interaction under the impacts of various water abstraction scenarios are presented for hydraulic fracturing in a shale gas and oil play area (23 984.9 km²), Alberta, Canada, using the MIKE-SHE and MIKE-11 models. Water-use data for hydraulic fracturing were obtained for 433 wells drilled in the study area in 2013 and 2014. Modelling results indicate that water abstraction for hydraulic fracturing has very small (<0.35%) negative impacts on mean monthly and annual river and groundwater levels and stream and groundwater flows in the study area, and small (1–4.17%) negative impacts on environmental flows near the water abstraction location during low-flow periods. The impacts on environmental flow depend on the amount of water abstraction and the daily flow over time at a specific river cross-section. The results also indicate a very small (<0.35%) positive impact on mean monthly and annual groundwater contributions to streamflow because of the large study area. The results provide useful information for planning long-term seasonal and annual water abstractions from the river and groundwater for hydraulic fracturing in a large study area.     

Groundwater flow systems in the arid crystalline province of Southern Sinai

ABSTRACT

A hydrogeological research project was conducted in the arid crystalline province of the Southern Sinai. To gain a better understanding of regional groundwater flow patterns, Landsat images as well as air photos were used. Three main flow systems were investigated: (1) The flow system occurring along longitudinal regional crustal fractures with the direction parallel to the Gulf of Elat Rift; this system captures the flow in the main wadi beds draining eastwards. (2) The flow system found in shallow fractured rocks, capturing the flow in small tributaries in the upper part of the drainage basin. (3) The system of dikes governing the flow regime in wadi beds draining westwards. The preliminary estimate of the recharge into the fractured crystalline aquifers and the interconnected gravel bed aquifers is in the order of magnitude of 10–15% of the total precipitation.     

A Model-Averaging Method for Assessing Groundwater Conceptual Model Uncertainty

This study evaluates alternative groundwater models with different recharge and geologic components at the northern Yucca Flat area of the Death Valley Regional Flow System (DVRFS), USA. Recharge over the DVRFS has been estimated using five methods, and five geological interpretations are available at the northern Yucca Flat area. Combining the recharge and geological components together with additional modeling components that represent other hydrogeological conditions yields a total of 25 groundwater flow models. As all the models are plausible given available data and information, evaluating model uncertainty becomes inevitable. On the other hand, hydraulic parameters (e.g., hydraulic conductivity) are uncertain in each model, giving rise to parametric uncertainty. Propagation of the uncertainty in the models and model parameters through groundwater modeling causes predictive uncertainty in model predictions (e.g., hydraulic head and flow). Parametric uncertainty within each model is assessed using Monte Carlo simulation, and model uncertainty is evaluated using the model averaging method. Two model-averaging techniques (on the basis of information criteria and GLUE) are discussed. This study shows that contribution of model uncertainty to predictive uncertainty is significantly larger than that of parametric uncertainty. For the recharge and geological components, uncertainty in the geological interpretations has more significant effect on model predictions than uncertainty in the recharge estimates. In addition, weighted residuals vary more for the different geological models than for different recharge models. Most of the calibrated observations are not important for discriminating between the alternative models, because their weighted residuals vary only slightly from one model to another.     

Characterization of the groundwater flow regime and hydrochemistry of groundwater from the Buem formation,Eastern Ghana

This study demonstrates the application of multivariate statistical methods in definition of groundwater recharge and discharge areas in a sedimentary basin in Ghana. Q‐mode hierarchical cluster analysis (HCA) was applied to 57 hydrochemical data from the Buem formation in the northern part of the Volta Region in Ghana. R‐mode HCA and R‐mode factor analysis were then applied to the same dataset to reveal the processes controlling the hydrochemistry of groundwater from this hydrogeological formation. Results of both the Q‐ and R‐mode analyses were backed by graphical methods. The analyses revealed two major water types, differentiated by salinity levels into four spatial groundwater associations. The characteristics of the four groundwater types are discussed. The recharge areas are characterized by Ca? HCO₃ low salinity waters which evolve through rock–water interactions to Na? HCO₃ high salinity waters in the discharge areas. This study finds that the hydrochemistry of groundwater from this formation is mainly controlled by the weathering of minerals, principally silicates in the aquifer matrix. The effects of the chemistry of recharging precipitation are higher in the recharge areas, while mineral weathering tends to be severe close to the discharge areas in the groundwater flow regime. All the four spatial groundwater associations have low sodium content, but salinity levels increase towards the discharge areas, such that some of wells in the discharge areas may not be acceptable for irrigation on grounds of high salinities which might affect the osmotic potentials of plants. Copyright © 2011 John Wiley & Sons, Ltd.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Multivariate approaches optimize locations of groundwater pumping facilities for different hydrogeological scales

A number of optimization approaches regarding the design location of groundwater pumping facilities in heterogeneous porous media have elicited little discussion. However, the location of groundwater pumping facilities is an important factor because it affects water resource usage. This study applies two optimization approaches to estimate the best recharge zone and suitable locations of the pumping facilities in southwestern Taiwan for different hydrogeological scales. First, for the regional scale, this study employs numerical modelling, MODFLOW‐96, to simulate groundwater direction and the optimal recharge zone in the study area. Based on the model's calibration and verification results, this study preliminarily utilizes the simulated spatial direction of groundwater and compares the safe yield for each well group in order to determine the best recharge zone. Additionally, for the local scale, the micro‐hydrogeological characteristics are considered before determining the design locations of the pumping facilities. According to drawdown record data from six observation wells derived from pumping tests at the best recharge area, this study further utilizes the modified artificial neural network approach to improve the accuracy of the estimation parameters as well as to analyse the direction and anisotropy of the hydraulic conductivities of an equivalent homogeneous aquifer. The results suggested that the best locations for the pumping facilities are along the more permeable major direction. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of Hydraulic Fracturing on Overlying Aquifers in the Presence of Leaky Abandoned Wells

The association between hydrocarbon‐rich reservoirs and organic‐rich source rocks means unconventional oil and gas plays usually occur in mature sedimentary basins—where large‐scale conventional development has already taken place. Abandoned wells in proximity to hydraulic fracturing could be affected by increased fluid pressures and corresponding newly generated fractures that directly connect (frac hit) to an abandoned well or to existing fractures intersecting an abandoned well. If contaminants migrate to a pathway hydraulically connected to an abandoned well, upward leakage may occur. Potential effects of hydraulic fracturing on upward flow through a particular type of leaky abandoned well—abandoned oil and gas wells converted into water wells were investigated using numerical modeling. Several factors that affect flow to leaky wells were considered including proximity of a leaky well to hydraulic fracturing, flowback, production, and leaky well abandonment methods. The numerical model used historical records and available industry data for the Eagle Ford Shale play in south Texas. Numerical simulations indicate that upward contaminant migration could occur through leaky converted wells if certain spatial and hydraulic conditions exist. Upward flow through leaky converted wells increased with proximity to hydraulic fracturing, but decreased when flowback and production occurred. Volumetric flow rates ranged between 0 and 0.086 m³/d for hydraulic‐fracturing scenarios. Potential groundwater impacts should be paired with plausible transport mechanisms, and upward flow through leaky abandoned wells could be unrelated to hydraulic fracturing. The results also underscore the need to evaluate historical activities.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

Performance Assessments of a Novel Well Design for Reducing Exposure to Bedrock‐Derived Arsenic

Arsenic in groundwater is a serious problem in New England, particularly for domestic well owners drawing water from bedrock aquifers. The overlying glacial aquifer generally has waters with low arsenic concentrations but is less used because of frequent loss of well water during dry periods and the vulnerability to surface‐sourced bacterial contamination. An alternative, novel design for shallow wells in glacial aquifers is intended to draw water primarily from unconsolidated glacial deposits, while being resistant to drought conditions and surface contamination. Its use could greatly reduce exposure to arsenic through drinking water for domestic use. Hypothetical numerical models were used to investigate the potential hydraulic performance of the new well design in reducing arsenic exposure. The aquifer system was divided into two parts, an upper section representing the glacial sediments and a lower section representing the bedrock. The location of the well, recharge conditions, and hydraulic properties were systematically varied in a series of simulations and the potential for arsenic contamination was quantified by analyzing groundwater flow paths to the well. The greatest risk of arsenic contamination occurred when the hydraulic conductivity of the bedrock aquifer was high, or where there was upward flow from the bedrock aquifer because of the position of the well in the flow system.     

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Aquifer interactions with a polluted mountain river of Nicaragua

The interactions between a stream and nearby shallow aquifers were investigated in a mountain basin being polluted by mercury released during mining in central Nicaragua. Hourly data series of water levels and temperatures were analysed using cross‐correlation. Resistivity imaging was used to map the subsurface and to complement the hydrological data interpretation. The results show the complex hydrogeological conditions that characterize the region, with weathering and fractured rock as main contributors to groundwater transport. The resistivity images suggest the presence of two vertical dykes perpendicular to the stream, and zones rich in clay. The data series indicate a rapid response from the aquifers to recharge events, followed by immediate discharge on a yearly basis. Furthermore, alternating periods of stream infiltration and aquifer discharge were identified. This work demonstrates that surface water pollution is a threat to groundwater quality in the area. Copyright © 2007 John Wiley & Sons, Ltd.     

Radionuclide Transport Behavior in a Generic Geological Radioactive Waste Repository

We performed numerical simulations of groundwater flow and radionuclide transport to study the influence of several factors, including the ambient hydraulic gradient, groundwater pressure anomalies, and the properties of the excavation damaged zone (EDZ), on the prevailing transport mechanism (i.e., advection or molecular diffusion) in a generic nuclear waste repository within a clay‐rich geological formation. By comparing simulation results, we show that the EDZ plays a major role as a preferential flowpath for radionuclide transport. When the EDZ is not taken into account, transport is dominated by molecular diffusion in almost the totality of the simulated domain, and transport velocity is about 40% slower. Modeling results also show that a reduction in hydraulic gradient leads to a greater predominance of diffusive transport, slowing down radionuclide transport by about 30% with respect to a scenario assuming a unit gradient. In addition, inward flow caused by negative pressure anomalies in the clay‐rich formation further reduces transport velocity, enhancing the ability of the geological barrier to contain the radioactive waste. On the other hand, local high gradients associated with positive pressure anomalies can speed up radionuclide transport with respect to steady‐state flow systems having the same regional hydraulic gradients. Transport behavior was also found to be sensitive to both geometrical and hydrogeological parameters of the EDZ. Results from this work can provide useful knowledge toward correctly assessing the post‐closure safety of a geological disposal system.     

Hydrogeological investigations in a low permeability claystone formation: the Mont Terri Rock Laboratory

This paper addresses the achievements in understanding the hydrogeological conditions in the low permeability claystone formation of the Opalinus Clay at the Mont Terri Rock Laboratory. The synthesis work consisted of (i) an assessment of clay-specific artifacts which may affect the interpretation of hydraulic tests, (ii) a survey of hydraulic rock properties such as hydraulic conductivity and storage coefficient, and (iii) an assessment of the governing flow laws and definition of hydrogeological units and flow boundaries on the site scale. Drawing on the broad hydrogeological data base, confidence is gained in the hydraulic barrier function of the Opalinus Clay, which effectively prevents groundwater flow between the over- and underlying aquifer systems.