Assessing hydrogeologic controls on dynamic groundwater storage using long‐term instrumental records of water table levels

This study analyzes a long‐term regional compilation of water table response to climate variability based on 124 long‐term groundwater wells distributed across New England, USA, screened in a variety of geologic materials. The New England region of the USA is located in a humid‐temperature climate underlain by low‐storage‐fractured metamorphic and crystalline bedrock dissected by north–south trending valleys filled with glacial and post‐glacial valley fill sediments. Uplands are covered by thin glacial till that comprises more than 60% of the total area. Annual and multi‐annual responses of the water table to climate variability are assessed to understand how local hydraulic properties and hydrogeologic setting (located in recharge/discharge region) of the aquifer influence the hydrologic sensitivity of the aquifer system to climate variability. This study documents that upland aquifer systems dominated by thin deposits of surface till comprise ~70% of the active and dynamic storage of the region. Total aquifer storage changes of +5 to ?7 km³ occur over the region during the study interval. The storage response is dominated by thin and low permeability surficial till aquifer that fills and drains on a multi‐annual basis and serves as the main mechanism to deliver water to valley fill aquifers and underlying bedrock aquifers. Whereas the till aquifer system is traditionally neglected as an important storage reservoir, this study highlights the importance of a process‐based understanding of how different landscape hydrogeologic units contribute to the overall hydrologic response of a region.     

Some relationships between lithology,basin form and hydrology: a case study from the Thames basin,UK

The role of lithology in influencing basin form and function is explored empirically by investigating correlations between a range of catchment variables, where the spatial unit of analysis is not surface catchments but lithologically coherent groundwater units. Using the Thames basin, UK, as a case study, nine groundwater units have been identified. Values for 11 hydrological and geomorphological variables, including rainfall, drainage density, Baseflow Index, aquifer porosity, storage coefficient and log‐hydraulic conductivity, aquifer and drainage elevation, river incision, and hypsometric integral have been estimated for each of the groundwater units in the basin, and Pearson correlation coefficients calculated for all pairs of variables. Seven of the correlation coefficients are found to be significant at a confidence level of > 99%. Negative correlations between drainage density and log aquifer hydraulic conductivity, and between drainage density and river incision, and positive correlations between log‐hydraulic conductivity and river incision, log‐hydraulic conductivity and Baseflow Index, and between Baseflow Index and river incision are inferred to have consistent causal explanations. For example, incision of rivers into aquifers leads to relative increases in hydraulic gradients in the vicinity of rivers which, in turn, promotes the development of secondary porosity increasing both aquifer hydraulic conductivity and, hence, Baseflow Index. The implication of this interpretation is that the geomorphological evolution of basins is intimately linked to the evolution of hydraulic conductivity of the underlying aquifers. This is consistent with, and supports the notion of a coupled complexly evolving surface water‐groundwater system. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of Groundwater Hydraulic Gradient on Bank Storage Metrics

The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.     

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.     

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.     

Modeling the spatial distribution of landslide‐prone colluvium and shallow groundwater on hillslopes of Seattle,WA

Landslides in partially saturated colluvium on Seattle, WA, hillslopes have resulted in property damage and human casualties. We developed statistical models of colluvium and shallow‐groundwater distributions to aid landslide hazard assessments. The models were developed using a geographic information system, digital geologic maps, digital topography, subsurface exploration results, the groundwater flow modeling software VS2DI and regression analyses. Input to the colluvium model includes slope, distance to a hillslope–crest escarpment, and escarpment slope and height. We developed different statistical relations for thickness of colluvium on four landforms. Groundwater model input includes colluvium basal slope and distance from the Fraser aquifer. This distance was used to estimate hydraulic conductivity based on the assumption that addition of finer‐grained material from down‐section would result in lower conductivity. Colluvial groundwater is perched so we estimated its saturated thickness. We used VS2DI to establish relations between saturated thickness and the hydraulic conductivity and basal slope of the colluvium. We developed different statistical relations for three groundwater flow regimes. All model results were validated using observational data that were excluded from calibration. Eighty percent of colluvium thickness predictions were within 25% of observed values and 88% of saturated thickness predictions were within 20% of observed values. The models are based on conditions common to many areas, so our method can provide accurate results for similar regions; relations in our statistical models require calibration for new regions. Our results suggest that Seattle landslides occur in native deposits and colluvium, ultimately in response to surface‐water erosion of hillslope toes. Regional groundwater conditions do not appear to strongly affect the general distribution of Seattle landslides; historical landslides were equally dispersed within and outside of the area potentially affected by regional groundwater conditions. Published in 2007 by John Wiley & Sons, Ltd.     

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.     

Steady state flow with hydraulic conductivity change around large diameter wells

Darcian flow law in aquifers assumes that the aquifer hydraulic conductivity is constant and the groundwater movement is due only to the piezometric level changes through hydraulic gradient. In practice, after the well development the aquifer just around the well has comparatively larger hydraulic conductivity and gradient. Patchy aquifer solutions in the literature consider sudden hydraulic conductivity changes with distance for the steady state flow. The change of transmissivity is demonstrated by the application of slope‐matching procedure to actual field data. It is the main purpose of this paper to derive simple analytical expressions for aquifer parameter evaluations with steadily decreasing hydraulic conductivity around the well. Spatial nonlinear hydraulic conductivity changes around a large‐diameter well within the depression cone of a confined aquifer are considered as exponentially decreasing functions of the radial distance. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Sea‐level rise impact on fresh groundwater lenses in two‐layer small islands

The fresh groundwater lenses (FGLs) of small islands can be highly vulnerable to climate change impacts, including sea‐level rise (SLR). Many real cases of atoll or sandy islands involve two‐layer hydrogeological conceptualizations. In this paper, the influential factors that affect FGLs in two‐layer small islands subject to SLR are investigated. An analytical solution describing FGLs in circular islands, composed of two geological layers, is developed for the simplified case of steady‐state and sharp‐interface conditions. An application of the developed model is demonstrated to estimate the FGL thickness of some real‐world islands by comparison with existing FGL thickness data. Furthermore, numerical modelling is applied to extend the analysis to consider dispersion effects and to confirm comparable results for both cases. Sensitivity analyses are used to assess the importance of land‐surface inundation caused by SLR, relative to other parameters (i.e. thickness of aquifer layers, hydraulic conductivity, recharge rate and land‐surface slope) that influence the FGL. Dimensionless parameters are used to generalize the findings. The results demonstrate that land‐surface inundation has a considerable impact on a FGL influenced by SLR, as expected, although the FGL volume is more sensitive to recharge, aquifer thickness and hydraulic conductivity than SLR impacts, considering typical parameter ranges. The methodology presented in this study provides water resource managers with a rapid‐assessment tool for evaluating the likely impacts of SLR and accompanying LSI on FGLs.     

Investigating Groundwater‐Lake Interactions by Hydraulic Heads and a Water Balance

Discharge of groundwater into lakes (lacustrine groundwater discharge, LGD) can play a major role in water balances of lakes. Unfortunately, studies often neglect this input path because of methodological difficulties in its determination. Direct measurements of LGD are labor‐consuming and prone to error. The present study uses both spatially variable hydraulic‐head data and meteorological data to estimate groundwater input by LGD and lake water output through infiltration. The study sites are two shallow, groundwater‐fed lakes without any surface inflows or outflows. Horizontally interpolated groundwater heads were combined with lake water levels to obtain vertical hydraulic gradients between the aquifer and the lake, which are separated by a thick layer of lake bed sediment which has a much lower hydraulic conductivity than the underlying aquifer. By fitting the hydraulic gradient to the results of a simple mass balance and considering the process of clogging, we were able to estimate the hydraulic conductivity of the lake bed sediments. We calculated groundwater inputs by LGD and lake water outputs by infiltration on an annual basis. Although our method requires several assumptions, the results are reasonable and provide useful information about the exchange between the aquifer and the lake, which can, for example, be used for the calculation of nutrient mass balances.     

Stage‐dependent transient storage of phosphorus in alluvial floodplains

Models for contaminant transport in streams commonly idealize transient storage as a well mixed but immobile system. These transient storage models capture rapid (near‐stream) hyporheic storage and transport, but do not account for large‐scale, stage‐dependent interaction with the alluvial aquifer. The objective of this research was to document transient storage of phosphorus (P) in coarse gravel alluvium potentially influenced by large‐scale, stage‐dependent preferential flow pathways (PFPs). Long‐term monitoring was performed at floodplain sites adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping which was correlated to hydraulic conductivity data, observation wells were installed both in higher hydraulic conductivity and lower hydraulic conductivity subsoils. Water levels in the wells were monitored over time, and water samples were obtained from the observation wells and the stream to document P concentrations at multiple times during high flow events. Contour plots indicating direction of flow were developed using water table elevation data. Contour plots of total P concentrations showed the alluvial aquifer acting as a transient storage zone, with P‐laden stream water heterogeneously entering the aquifer during the passage of a storm pulse, and subsequently re‐entering the stream during baseflow conditions. Some groundwater in the alluvial floodplains had total P concentrations that mirrored the streams' total P concentrations. A detailed analysis of P forms indicated that particulate P (i.e. P attached to particulates greater than 0·45 µm) was a significant portion of the P transport. This research suggests the need for more controlled studies on stage‐dependent transient storage in alluvial systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Seasonal variation of infiltration rates through pond bed in a managed aquifer recharge system in St-André, Belgium

In Belgium, IWVA uses managed aquifer recharge (MAR) to recharge the aquifer with treated wastewater generated from the communities to sustain the potable water supply on the Belgian coast. This MAR facility is faced with a challenge of reduced infiltration rates during the winter season when pond water temperatures near 4°C. This study involves the identification of the predominant factor influencing the rate of infiltration through the pond bed. Several factors, including pumping rates, natural recharge, tidal influences of the North Sea and pond-water temperature, were identified as potential causes for variation of the recharge rate. Correlation statistics and linear regression analysis were used to determine the sensitivity of the infiltration rate to the aforementioned factors. Two groundwater flow models were developed in visual MODFLOW to simulate the water movement under the pond bed and to obtain the differences in flux to track the effects of variation of hydraulic conductivity during the two seasons. A 32% reduction in vertical hydraulic gradient in the top portion of the aquifer was observed in winter, causing the recharge rates to fluctuate. Results showed that water temperature caused a 30% increase in hydraulic conductivity in summer as compared with winter and has the maximum impact on infiltration rate. Cyclic variations in water viscosity, occurring because of seasonal temperature changes, influence the saturated hydraulic conductivity of the pond bed. Results from the models confirm the impact on infiltration rate by temperature-influenced hydraulic conductivity.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

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.     

The role of antecedent groundwater heads in controlling transient aquifer storage and flood peak attenuation in karst watersheds

Transient storage of floodwaters in aquifers is known to attenuate peak flows in rivers and drive subsurface dissolution. Transient aquifer storage could be enhanced in watersheds overlying karst aquifers where caves facilitate surface and groundwater exchange. Few studies, however, have examined controls on, or magnitudes of, transient aquifer storage or flood peak attenuation in karstic watersheds. Here we evaluate flood peak attenuation with multiple linear regression analyses of 10 years of river and groundwater data from the Suwannee River, which flows over the karstic upper Floridan aquifer in north-central Florida and experiences frequent flooding. Regressions show antecedent river stage exerts the dominant control on magnitudes of transient aquifer storage, with recharge and time to peak having secondary controls. Specifically, low antecedent stages result in larger magnitudes of transient aquifer storage and thus greater flood attenuation than conditions of elevated antecedent stage. These findings suggest subsurface weathering, including cave formation and enlargement, caused by transient aquifer storage could occur on a more frequent basis in aquifers where groundwater table elevation is lowered due to anthropogenic or climatic influences. Our work also shows that measures of groundwater table elevation prior to an event could be used to improve predictive flood models. © 2018 John Wiley & Sons, Ltd.     

Integrated water resource management in a major canal command in eastern India

The present rice‐dominated cropping system in the Hirakud canal command (eastern India) is under severe threat due to imbalance between irrigation water supply and demand. The canal water supply, which is the only source of irrigation, only meets 54% of the demand at 90% probability of exceedance (PE). In order to mitigate the irrigation water deficit from canal water, groundwater is considered as a supplemental source. Quasi‐three‐dimensional groundwater flow simulation modelling was, therefore, carried out by using Visual MODFLOW to detect the change in hydraulic head due to transient pumping stresses. The simulation model was calibrated and validated satisfactorily. Sensitivity analysis of the model parameters shows that groundwater recharge is most sensitive followed by aquifer hydraulic conductivity at almost all the sites of the command area, whereas the model is comparatively less sensitive to specific storage and specific yield. Enhanced pumping scenarios showed that groundwater extraction can be increased up to 50 times of the existing pumping without causing any adverse effect to the aquifer but the aquifer does not permit to exploit water in order to fulfill the irrigation water demand even at 10% PE. Hence, it is imperative to develop an optimal land and water resources management plan of the command area. Copyright © 2011 John Wiley & Sons, Ltd.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Compartmentalisation and groundwater–surface water interactions in a prospective shale gas basin: Assessment using variance analysis and multivariate statistics on water quality data

An environmental concern with hydraulic fracturing for shale gas is the risk of groundwater and surface water contamination. Assessing this risk partly involves the identification and understanding of groundwater–surface water interactions because potentially contaminating fluids could move from one water body to the other along hydraulic pathways. In this study, we use water quality data from a prospective shale gas basin to determine: if surface water sampling could identify groundwater compartmentalisation by low-permeability faults; and if surface waters interact with groundwater in underlying bedrock formations, thereby indicating hydraulic pathways. Variance analysis showed that bedrock geology was a significant factor influencing surface water quality, indicating regional-scale groundwater–surface water interactions despite the presence of an overlying region-wide layer of superficial deposits averaging 30–40 m thickness. We propose that surface waters interact with a weathered bedrock layer through the complex distribution of glaciofluvial sands and gravels. Principal component analysis showed that surface water compositions were constrained within groundwater end-member compositions. Surface water quality data showed no relationship with groundwater compartmentalisation known to be caused by a major basin fault. Therefore, there was no chemical evidence to suggest that deeper groundwater in this particular area of the prospective basin was reaching the surface in response to compartmentalisation. Consequently, in this case compartmentalisation does not appear to increase the risk of fracking-related contaminants reaching surface waters, although this may differ under different hydrogeological scenarios.