Vertical Hydraulic Conductivity Research Articles

A semi-analytical solution for transient confined-unconfined flow toward the partially penetrating well

In this paper, a semi-analytical solution has been proposed for the transient confined-unconfined flow toward a partially penetrating well in a confined aquifer. It can be used to characterize the dynamic development of drawdown and the transient unconfined region during the transient confined-unconfined conversion. The proposed solution takes into account the delayed response recharge from the unsaturated region by the Neuman's upper surface recharge boundary condition, which is derived by the Laplace and Fourier transformations. The reliability of the semi-analytical solution has been proved by comparisons with numerical solutions from Visual MODFLOW Flex. In addition, a new model for the pumping flow in a pure unconfined aquifer is also developed by degradation analysis of the proposed solution, which has been used to fit well with the measured data from the pumping test in Cape Cod aquifer. The results indicate that both the use of partially penetrating pumping well and aquifer anisotropy have negative effects on the drawdown delayed responses and the development of transient unconfined region, which are disappeared at times >105Srb/Kr. The delayed response is positively related to the gravity storage and vertical hydraulic conductivity.

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

AbstractLateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

A Large-Sized Permeameter for Studying Suffusion Characteristics of Anisotropic Soils

Abstract Seepage-induced suffusion involves the migration of fine particles within a soil matrix. Seepage flow is affected by the soil permeability anisotropy of anisotropic soil fabric; however, suffusion anisotropy is unclear because of the limited function of existing permeameters. In recent studies, the effect of seepage direction has been investigated under only low hydraulic gradients because the control of seepage direction relies merely on gravity. In this study, a new, large-sized permeameter is developed with which suffusion tests can be conducted along horizontal or vertical seepage directions under high hydraulic gradients. Correspondingly, the permeameter can accommodate a specimen of 540 × 500 × 470 or 540 × 540 × 440 mm3 (length × width × height). The seepage direction is switched by changing the boundary conditions of the specimen with detachable perforated plates that allow pressurized water originating from different inlets to flow along horizontal or vertical directions. Two repeated pairs of tests were performed on a gap-graded clayey gravel to investigate the suffusion anisotropy of saturated clayey gravel. The results show that the maximum relative deviations of measurements for initial hydraulic conductivity, initiation, and failure hydraulic gradients are less than 3.5 %, demonstrating satisfactory reliability. The ratio of the initial horizontal hydraulic conductivity to vertical hydraulic conductivity for the test soil is 13.87, indicating a significantly anisotropic fabric induced by compaction. The ratios of horizontal initiation and failure hydraulic gradients to vertical initiation and failure hydraulic gradients are 0.52 and 0.59, respectively. This implies that suffusion anisotropy should not be neglected for evaluating the internal instability of anisotropic soils.

Unsteady flow modeling of low-velocity non-Darcian flow to a partially penetrating well in a leaky aquifer system

Previous studies on the flow dynamics of leaky aquifer systems have primarily relied on Darcy's law, assuming fully penetrating wells in confined aquifers. However, when dealing with aquitards dominated by clay, the flow behavior often deviates from Darcy's law. Moreover, in the majority of cases, this entails partially penetrating well pumping. This paper introduces a mathematical model for non-Darcian flow in a confined and partially penetrating well system within a leaky aquifer system. The model is based on the low-velocity non-Darcian flow equation, incorporating the threshold pressure gradient. In the confined aquifer, the flow exhibits two-dimensional Darcian flow behavior, while in the aquitard, the flow is characterized by one-dimensional non-Darcian flow in the vertical direction. In the shallow aquifer, the flow is one-dimensional Darcian flow in the radial direction. The mathematical model is solved using the finite difference method, and the obtained results are compared with calculations based on traditional Darcian flow theory. The findings indicate that the confined groundwater head difference, as predicted by the Darcian and non-Darcian flow theories, diminishes radially as the distance from the pumping well increases and vertically as one moves away from the top of the confined aquifer. As pumping progresses, the confined groundwater head difference initially rises and subsequently stabilizes gradually. Moreover, the magnitude of the confined groundwater head difference is more pronounced when the threshold pressure gradient, vertical hydraulic conductivity of the aquitard, pumping rate, hydraulic conductivity and specific yield of shallow aquifer are larger, or when the hydraulic conductivity and specific storage of the confined aquifer are smaller. The effect of well screen length on the confined groundwater head difference is minimal.

Numerical modeling of vadose zone electrical resistivity to evaluate its hydraulic parameters

Studying and determining the physical properties and hydraulic parameters of vadose zone sediments is an important key to evaluate the infiltration rate into them and assessing the extent to which aquifer sediments benefit from rainwater harvesting in arid and semi-arid areas. Due to the lack of sufficient data on the characteristics of this zone depths, a numerical modeling was used to simulate the electrical resistivity of these sediments by applying the electrical resistivity method, because it is the most affected by the physical properties of dry and wet sediments. This study was applied as a proposal for application in northwestern KSA to calculate the vertical hydraulic conductivity and transmissivity for the vadose zone. This was implemented by assuming a three-layer model using COMSOL Multiphysics model with different electrical resistivity values depending on some in situ electrical resistivity measurements for shallow depths. Hence, the infiltration rate of sediments in this area can be predicted with depth and its effect on aquifer recharge. The focus was on calculating the vertical hydraulic parameters of the most widespread surface sediments with depth and comparing the results of calculating these parameters for some sediments laboratory-wise to ensure their accuracy. Then, their infiltration rate was inferred separately with depth, predicting their ability to aquifer recharge and make the most of rainwater harvesting. Finally, this study can be considered as a preliminary study to determine the expected forward model of electrical resistivity and hydraulic parameters values for the vadose zone sediments with depth along the area and in any other areas, and then apply them accurately in situ to estimate the extent of its usefulness in rainwater harvesting, especially aquifer recharge.

Influence of landscape position and climatic seasonality on soil water and gas conductivity properties in agricultural soils

Agricultural landscape management and climate seasonality can influence soil structure, hydraulic conductivity, and air permeability within the context of soil water and soil gas mobility. To investigate this, in situ and laboratory-based data were collected from three agricultural landscape positions within a watershed in eastern Ontario, Canada during a growing season. Macropore classification, water infiltration tests, and air permeability measurements were conducted in situ and standard soil characterizations were carried out on soil samples. Hydraulic conductivity of the soil matrix, based on grain size data, indicated that the highest values were consistently measured in the B horizon at each landscape setting. Macropores were found to be more abundant within uncultivated drainage ditch bank soils, compared to the adjacent cropped fields. Macropores in the ditch bank soils were exclusively consisted of circular biopores, while both circular and linear macropores were observed in the cultivated field soils. Air permeability, vertical hydraulic conductivity, and horizontal hydraulic conductivity were also greater in the uncultivated soils, relative to the cultivated soils. Field saturated hydraulic conductivity measurements offered evidence of anisotropy, likely due to the vertical nature of the macropore features. Macropore disposition and extent varied over the growing season, especially in the cultivated field soils where tillage and field trafficking are physically disruptive. Seasonality of macropore development will influence temporal changes in advection-based mass exchange of gas and water in the vadose zone. Modeling of mass exchange in agricultural soils should consider time variability in macroporosity to more realistically characterize infiltration and soil gas emissions.

Kaolinite Deposition Dynamics and Streambed Clogging During Bedform Migration Under Losing and Gaining Flow Conditions

AbstractClogging of streambeds due to clay deposition influences the stream‐subsurface exchange flux and thus directly modulates hyporheic ecological and biogeochemical processes. Clogging of sandy streambeds has previously been studied under losing and gaining flows and during streambed movement, but not when these two flow conditions coincided. We conducted flume experiments to quantify the combined effect of moving bedforms and losing or gaining flows on kaolinite deposition and streambed clogging. The experiments were conducted by adding pulses of kaolinite in a flume packed with sand under a stream water velocity of 25 cm/s. We measured the deposition rates, dynamics of hyporheic exchange flux (HEF) and vertical hydraulic conductivity (Kv), and the vertical distribution of kaolinite at the end of the experiments under two losing and two gaining flows (Darcy velocity of 10 and 20 cm/day). Kaolinite deposition led to clogging and reduction in Kv and HEF under all flow conditions. Deposition occurred faster under losing flow conditions than under gaining flow conditions. However, the changes in Kv and HEF were similar under losing and gaining flow conditions for similar kaolinite concentrations in the bed. Our results indicate that the deposition patterns of kaolinite were more influenced by bedform movement than by losing or gaining flow conditions, which is markedly different from the behavior observed under losing and gaining conditions for stationary bedforms. This implies that bedform morphodynamics control local‐scale clogging of sandy streambeds and should be accounted for when studying the hydrology of catchments at larger scales.

An Application of Inverse Problem and Universal Solutions for Pumping Wells in Unconfined Aquifers

As far as we know, universal solutions (or type-curves) for scenarios of flow through anisotropic unconfined aquifers due to pumping wells cannot be found in the literature. On the contrary, those theoretical solutions in hydrogeological manuals are commonly based on Dupuit solutions for isotropic soils or simplifying other characteristics of the chosen medium. In this study, the application of the discriminated nondimensionalization technique allowed for the inclusion of vertical and radial hydraulic conductivities in the data set, with which the monomials ruling unknown variables of the problem, pumping flow and seepage surface in their dimensionless form are obtained. One of the main findings of this research is depicting these relationships as type-curves from a large number of precise numerical simulations based on the Network Simulation Method. The other main finding is an easy-to-apply methodology to estimate vertical and radial hydraulic conductivities employing these type-curves. This methodology can be considered as an inverse problem. In addition, an example of the problem is presented, in which the influence that measure deviations may have on the estimated values of the hydraulic conductivities in anisotropic soils is also studied and discussed.

Numerical investigation of coupled density-driven flow and hydrochemical processes in coastal aquifers influenced by storm surges

Storm surges can result in significant damage to coastal aquifers. Studies have investigated groundwater salinization caused by storm surge inundation. However, the phenomenon of evaporation-driven mineral precipitation in saltwater lakes has not received much attention in the literature. To address this gap, this study simulated the impact of saltwater lakes formed after a storm surge on coastal aquifers by coupling the PHREEQC and SEAWAT models. Six scenarios were considered, with varying lake depth, area, evaporation rate and vertical hydraulic conductivity of lakebed. The results indicate that the total dissolved solids (TDS) and water level of saltwater lake are subject to constant changes due to evaporation. Scenario 2, with the smallest lake area exhibited the least susceptibility to evaporation, with a maximum TDS concentration in the lake of ∼86 g/L, and only dolomite and calcite reached saturation. However, infiltration of saline water into the aquifer was higher than for the other scenarios. A comparing of scenarios with different vertical hydraulic conductivities of lakebed (scenarios 5 and 6 have conductivities that are 0.1 × and 0.01 × that of scenario 1, respectively) shows that a smaller vertical hydraulic conductivity leads to less saline infiltration into the aquifer, but a greater amount of salt precipitates (precipitation ratios of 1.03%, 1.92%, and 5.73% for scenario 1, 5 and 6, respectively). The results of scenario 4, which has double the evaporation rate of the other scenarios, indicate that evaporation leads to significantly greater salt precipitation and a corresponding decrease in the amount of salt that enters the aquifer. The results of secondary salinization indicate that salt precipitated on the surface may infiltrate the aquifer through rainfall dissolution, causing groundwater to become saline again. These findings emphasize the importance of accounting for the impact of evaporation on saltwater lake dynamics and coastal groundwater salinization in future studies.

Improvement of in situ hydraulic conductivity tests for riverbeds under the influence of infiltration intake: A case study of the Oława River, SW Poland

The vertical hydraulic conductivity (KV) of riverbeds plays a pivotal role in controlling water exchanges between the surface water and groundwater. The recent papers have focused on the spatial and temporal variability of KV (K measured in the “z” direction), while point measurement schemes seem to require more attention. The most popular measurement technique is the falling head method, in which a PVC pipe is placed in the riverbed. The most commonly used formulas to determine hydraulic conductivity are those based on the Darcy's law and its modifications. This article presents how the falling head method should be carried out and what errors need to be considered in conditions of the riverbank intake. The simultaneous measurements and analyses of the falling head vs. the riverbed head can produce satisfying outcomes and reduce errors. A full utility MS Excel spreadsheet based on an automated formula with a solver was developed.

Parameter identification and range restriction through sensitivity analysis for a high-temperature heat injection test

In order to compensate for the variable mismatch between heat demand and heat production from renewable sources or waste heat, high-temperature aquifer thermal energy storage (HT-ATES) is a promising option. A reliable prediction of the energetic performance as well as thermal and hydraulic impacts of a HT-ATES requires a suitable model parameterization regarding the subsurface properties. In order to identify the subsurface parameters on which investigation efforts should be focused, we carried out an extensive sensitivity analysis of the thermal and hydraulic parameters for a high-temperature heat injection test (HIT) using numerical modeling of the governing coupled thermo-hydraulic processes. The heat injection test was carried out in a quaternary shallow aquifer using injection temperatures of about 75 °C over 5 days, accompanied by an extensive temperature monitoring. The sensitivity analysis is conducted for parameter ranges based on literature values, based on site investigation at the HIT site and based on a model calibrated to the measured temperature distribution following the heat injection. Comparing the parameter ranges thus obtained in this three-step approach allows to identify those parameters, for which model prediction uncertainty decreased most, which are also the parameters, that strongly affect the thermal behavior. The highest sensitivity is found for vertical and horizontal hydraulic conductivity as well as for groundwater flow velocity, indicating that investigation efforts for HT-ATES projects should focus on these parameters. Heat capacity and thermal conductivity have a smaller impact on the temperature distribution. Our work thus yields a consistent approach to identifying the parameters which can be best restricted by field investigations and subsequent model calibration. Focusing on these during field investigations thus enable improved model predictions of both HT-ATES operation and induced impacts.

An intrinsic vulnerability approach to assess an overburden alluvial aquifer exposure to sinkhole-prone area; results from a Central Iran case study

DRASTIC is a model used to reliably assess the pollution potential of aquifers at global location. Despite various criticisms of this model or its optimization, sometimes the territorial or geological conditions still must rely on this model approach based on its simplicity. The existence of sinkholes in the alluvial aquifer with an underlying karst aquifer is a clear example of the need to change or optimize the DRASTIC model to improve accuracy of prediction. In this study, the DRASTIC model was changed by adding two factors, the distance to a sinkhole and the sinkhole catchment factors. The modified DRASTIC model was then used to evaluate the pollution potential of Abarkouh aquifer in Yazd province, Central Iran. The intrinsic vulnerability index of Abarkouh aquifer was calculated by the weighted sum method in GIS. The resulting numerical values ranged between 62 and 176 for full study area. Without adding the two layers related to Sin-DRASTIC, the generic DRASTIC Index map showed summed values ranging between 54 and 130. The final Sin-DRASTIC map was divided it into 4 orders very low (<106; 722 km2, 78%), low (107-133; 186 km2, 20%), low to moderate (133-160; 15.3 km2, 1.7%), and moderate vulnerability (161-176; 1.7 km2, 0.3%) vulnerability. Both of the created maps showed high Area Under the Curve values at 76% and 78% for the generic DRASTIC and Sin-DRASTIC respectively. However, the area of the aquifer containing the sinkholes showed low and very low zones of contamination vulnerability on the generic DRASTIC map. The sensitivity analysis showed that the aquifer vulnerability to pollution is mainly controlled by the vadose zone impacts and the sinkhole catchment factor. Therefore, it is concluded that for geologic settings containing sinkholes that are open or contain higher vertical hydraulic conductivity sediment compared to the general matrix of the unconfined aquifer, the Sin-DRASTIC model method yields a superior prediction.

Improved Darcian streambed measurements to quantify flux and mass discharge of volatile organic compounds from a contaminated aquifer to an urban stream

Quantifying VOC transport from contaminated groundwater to streams is challenging and important for understanding off-site migration of VOCs, cross-media contamination (groundwater to surface water and eventually air), and potential impacts on downstream ecosystems and human populations. A streambed point sampling approach was used to quantify fluxes of water and 14 VOCs from groundwater to an urban stream in North Carolina, USA, during summer (June 2015) and winter (January 2016). The approach is unique in coupling measurements of vertical hydraulic conductivity, vertical hydraulic head gradient, and groundwater VOC concentration at each individual sampling point, reducing or eliminating some potential concerns with other Darcian methods for quantifying VOC inputs to streams.Most results were consistent with discharge of two main VOC plumes on opposite sides of the stream. Plume 1 from the west side was dominated by cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC) at mean concentrations of 19 and 11 μg L−1, respectively. Plume 2 from the east side was dominated by benzene (mean concentration 56 μg L−1). Plume 2 was not previously known, and the improved sampling approach allowed VOC discharge from both plumes to be quantified simultaneously.For 13 of the 14 detected VOCs, the mean VOC flux from groundwater to the stream (fVOC) was higher in January 2016 than in June 2015, mainly because groundwater flux was higher in January. The only exception was cDCE, the most abundant VOC in Plume 1, which had mean fVOC values of 9.8 and 9.5 mg m−2 d−1 in June 2015 and January 2016, respectively. Benzene was the most abundant VOC in Plume 2 and had mean fVOC values of 11 and 37 mg m−2 d−1 in June 2015 and January 2016, respectively. High groundwater flux drove almost all the occurrences of high VOC flux. For a given VOC, the flow-weighted mean concentration (with each VOC concentration weighted by the upward groundwater flux at the VOC sampling point) was generally larger than the unweighted mean concentration. Thus, flow-weighting of concentrations gave a more accurate indication of the average VOC concentration in net groundwater discharge to the stream.An estimate of total VOC mass discharge from groundwater to the study reach of the stream, 3.6 kg of VOC per year, was based on the fVOC results and streambed area in the reach. The bulk of this discharge was due to benzene, cDCE, and VC, with individual mass discharges of 2.1, 0.83, and 0.40 kg yr−1, respectively. Estimates of maximum potential VOC degradation in the streambed suggest that the 3.6 kg yr−1 estimate of mass discharge was not sensitive to potential degradation of VOCs in the streambed sediments above the groundwater sampling depth.

Temporary clogging effects induced by a sustainable anti-icing hydrogel on the hydraulic conductivity and inertia coefficient of open-graded asphalt pavements

Open-graded asphalt pavements require special winter maintenance procedures to ensure the effectiveness of anti-icing and deicing measures. A recently developed solution for this type of pavement is represented by an environmental friendly thickened bio-based salt hydrogel, which consists of a thermo-sensitive sodium chloride (NaCl) brine admixed with a gellant agent (seaweed fibers) having the ability to form a gel-like structure when hot-sprayed on a cold surface. This thin salt layer results in a long-lasting residual efficacy of the winter maintenance product and in a reduction of salt consumption, but affects the drainage capacity of the pavement. These transient clogging effects were evaluated on laboratory-made specimens through the measurement of their vertical hydraulic conductivity variations before and after the hydrogel application, using a 1-D air permeameter and a constant head water permeameter. The analyses, based on Darcy–Forchheimer model and Darcy model for lower Reynold numbers, revealed a significant short-term reduction of the hydraulic conductivity, that, anyway, tends to be almost totally restored once the critical weather event is over (after about 6 h).

Spatial variability of the vertical saturated hydraulic conductivity of sediments around typical thermokarst lakes

The formation and development of thermokarst lakes is a key factor causing permafrost degradation, in which the vertical hydraulic conductivity (Kv) of lakebed sediments plays an important role as it affects the water and heat transfer between thermokarst lakes and permafrost. However, the spatial variability of Kv of lakebed sediments in the Qinghai–Tibet Plateau (QTP) remains unclear. Hence, in this study, the Kv values of sediments in two thermokarst lakebeds were measured in situ using the standpipe method, and the grain size of the sediments was analyzed. The statistical results showed that BLH-A and B (lakes A and B in Beiluhe) had moderate and strong spatial variabilities, respectively. Lilliefors test results showed a log-normal distribution of the Kv measurements of the sediment samples. From the results of the geostatistical method, it was concluded that the Gaussian model best fitted the variograms of the Kv values. The d50 (grain size for which 50 % of grains are finer) had the strongest correlation with Kv, and the relationship between Kv and the grain size of d50 was a quadratic function. The estimated Kv value (0.03 m/d) was close to the measurements (0.005–0.58 m/d), and the specific discharge of BLH-B ranged from − 0.0038 to 0.00056 m/d. BLH-A and BLH-B would gradually expand under the influence of future climate due to the special geology (low Kv). This basic research on thermokarst lake regions can help understand the properties of thermokarst lakebed sediments and provide parameters for the numerical simulations of the hydrological cycle in thermokarst lake regions.

Groundwater flow numerical model to evaluate the water mass balance and flow patterns in Groundwater Circulation Wells (GCW) with varying aquifer parameters

Groundwater Circulation Wells (GCW) can be an effective in-situ remediation option allowing high mass recovery of contaminants in cases where contamination hotspots are located in saturated soil having low hydraulic conductivity. Traditional treatment options such as Pump&amp;Treat, Air Sparging (AS)/Soil Vapor Extraction (SVE) and Multi Phase Extraction (MPE) typically require long operation times and significant costs for long-term plume management. GCWs induce meaningful changes in the groundwater flow introducing vertical flows both downward and upward, generating a “circulation cell”, which facilitates contaminant desorption from the soil. This study aims to understand the effects of a GCW on an aquifer in terms of both groundwater flow directions and water balance. A groundwater numerical model was built using MODFLOW-2005 to simulate the effect of the hydraulic parameters of the aquifer on the hydraulic circulation pattern of the GCW. The use of particle tracking simulated by MODPATH 7 showed the circulation cells and the impact on groundwater directions induced by different configurations of hydraulic parameters. The water flowing into the cell comes from both the injection well and the surrounding aquifer and the model shows how the hydraulic parameters of the aquifer, in particular the horizontal and vertical hydraulic conductivity, have a paramount influence in determining the shape and dimension of the circulation cell. A water mass balance analysis was carried out. It allowed to predict the groundwater flows exchanges between the GCW system and the surrounding aquifer, and to verify the sensitivity of the water budget to specific aquifer parameters. The results of this study are useful for further understanding the hydraulics of a GCW remediation system in order to support the design and to predict its performance.

Seasonal change of groundwater response to Earth tides

• Tidal response of groundwater may change from being confined-like to unconfined-like. • The tidal change may be attributed to seasonal vertical conductivity changes in aquitard. • Saturation changes in fractures could result in the vertical conductivity changes. Groundwater is an important natural resource and its security depends largely on the confinement of the aquifer. Here we use ten-years of continuous groundwater data from a well in Southwest China to show that seasonal precipitation may change the aquifer's response to Earth tides from being confined-like during dry seasons to being unconfined-like during rainy seasons. Comparing the response of water level in this well with the theoretical prediction from a leaky aquifer model, we show that the observed changes of aquifer confinement reflect a seasonal increase of the vertical hydraulic conductivity from being low in the dry season when the water level is low to being high in the rainy season when the water level is high. Such seasonal change of aquifer confinement has not been reported before and may impact our understanding on the subsurface transport, the security of groundwater resources, and the safety of toxic waste repositories.

Influence of inter-aquifer leakage on well-injection capacity: Theory and aquifer-scale mapping for artificial recharge

Aquifer storage and recovery (ASR) is an important water resources management technique that involves the injection of a large volume of water underground. For the successful implementation of an ASR project, a target aquifer should have a sufficient injection capacity, which is the maximum volume of water that can be safely injected. In nature, no aquitard is perfectly impermeable, and inter-aquifer leakage may have a major impact on injection capacity. Despite the importance of determining the injection capacity for ASR planning, there is no quantitative methodology that estimates the injection capacity of leaky aquifers. In this study, we first develop a solution for injection capacity with inter-aquifer leakage based on the Hantush - Jacob solution, and conduct a comprehensive sensitivity analysis to elucidate the influence of inter-aquifer leakage on injection capacity. From the sensitivity analysis, we show that inter-aquifer leakage can impact injection capacity by more than one order of magnitude, depending on the hydrogeological and operational parameters. We then develop a practical mapping methodology that estimates the injection capacity of leaky aquifers. We demonstrate the proposed methodology by applying it to a potential ASR site in Minnesota, USA, where ASR is considered as a solution to alleviate groundwater contamination by PFAS chemicals. The case study results reveal significant spatial variability in injection capacity over the study area and show an average increase in the injection capacity of about 26% compared to that in the nonleaky scenario. We also analyze the uncertainty in the estimated injection capacity due to the variability of aquitard properties and show that the variability of aquitard vertical hydraulic conductivity leads to a larger uncertainty in the estimated injection capacity than does the variability of aquitard thickness. This study elucidates the effects of inter-aquifer leakage on injection capacity and provides a practical methodology for injection capacity mapping.

Parametric study of fluidization conditions of bed sediments in the hyporheic zone

AbstractThe phenomenon of fluidization of the bed sediment of the hyporheic zone involves the production of dynamic suspension of solid particles in a stream of water flowing in the direction opposite to gravity. Fluidization may be permanent; then, it becomes a spring or a headwater zone, where groundwater discharge rate is almost unchanged within an entire hydrological year. But it can also occur periodically, depending on the changing hydrodynamic conditions. In such cases, the fluidization of the bed sediment of the hyporheic zone may disturb its functioning by changing the conditions of surface water and groundwater exchange. The disturbances affect hydrodynamic and hydrochemical processes connected with water exchange, but first of all, the living conditions of hyporheic zone flora and fauna. Therefore, the phenomenon of fluidization is significant for broadly understood ecology of the hyporheic zone. The paper presents the analysis of conditions of fluidization occurrence depending on the main hydrodynamic parameters that may affect this process. The authors analysed the influence of the difference in water level between the drained aquifer and the river, the hydraulic conductivity and anisotropy of the aquifer, the vertical hydraulic conductivity in the hyporheic zone, thickness of the hyporheic zone, and the width of the river on the possibility of fluidization of the bed sediment of the hyporheic zone. The range of variability of particular parameters corresponds to the actual variability of those parameters in natural conditions. It was proved that the fluidization of the bed sediment of the hyporheic zone mostly depends on: water level difference, hydraulic conductivity of the aquifer, and value of vertical hydraulic conductivity in the hyporheic zone. The other parameters have little influence on this process.

A low-cost model for slug tests in a confined aquifer with skin-zone effect

While the effect of small-scale skin has been extensively investigated, numerical models incorporating the effect require expensive computational cost due to fine spatial discretization of the skin and have not been widely available. This study develops a new model for slug tests at a partially-penetrating well in a two-zone confined aquifer of skin and formation zones. A new basic formulation defines transient flow in the skin with an additional skin factor reflecting the skin storage effect. The formulation assumes horizontal flow in the skin with excluding the vertical skin hydraulic conductivity. The analytical solution and finite element solution of the model are developed. The analytical expression relating the skin factor to the skin specific storage is derived. Results suggest the formulation is a useful alternative to conventional flow equation for the skin in achieving coarse spatial discretization and low computational cost. The present solution agrees well with the field head data obtained from Butler (2019). To conclude, this study presents theoretical advance and efficient simulation in modeling of well hydraulics.