Tide-induced airflow in a two-layered coastal land with atmospheric pressure fluctuations

Tide-induced airflow is commonly seen in coastal lands and affects ground stability especially with a less permeable pavement on the ground surface. A tide-induced airflow model in a two-layered unsaturated zone consisting of a highly permeable layer underneath a less permeable layer was established by Li and Jiao [Li HL, JJ Jiao. One-dimensional airflow in unsaturated zone induced by periodic water table fluctuation. Water Resour Res 2005;41:W04007. doi:10.1029/2004WR003916] to describe the one-dimensional airflow with constant atmospheric pressure at the ground surface. In this study, we expand the Li and Jiao model by considering the realistic atmospheric pressure fluctuations and the initial condition. A new transient solution to the airflow model is developed for an initial boundary value problem (IBVP). The transient solution can be used not only to calculate the subsurface air pressure at a future time with a known initial condition, but also to evaluate the asymptotic air pressure variations when time becomes long. The amplitude ratio and phase lag of the subsurface air pressure relative to the tide-induced hydraulic head variations inside the unconfined aquifer below the unsaturated zone are investigated. The results reveal that effect on the subsurface pressure due to changes of atmospheric pressure amplitude depends on the configurations of air resistance in the less permeable layer and the air-filled porosity difference in the two layers. The introduction of atmospheric pressure fluctuations into the airflow model leads to insignificant influence on water table level. A field application of the new solution at Hong Kong International Airport in Hong Kong, China is demonstrated. It indicates that the new transient solution can be conveniently used to evaluate the subsurface air pressure with discrete atmospheric pressure data at the ground surface.     

Modelling coastal ground‐ and surface‐water interactions using an integrated approach

Beach water table fluctuations have an impact on the transport of beach sediments and the exchange of solute and mass between coastal aquifer and nearby water bodies. Details are given of the refinement of a dynamically integrated ground‐ and surface‐water model, and its application to study ground‐ and surface‐water interactions in coastal regions. The depth‐integrated shallow‐water equations are used to represent the surface‐water flow, and the extended Darcy's equation is used to represent the groundwater flow, with a hydrostatic pressure distribution being assumed to apply for both these two types of flows. At the intertidal region, the model has two layers, with the surface‐water layer being located on the top of the groundwater layer. The governing equations for these two types of flows are discretized in a similar manner and they are combined to give one set of linear algebraic equations that can be solved efficiently. The model is used to predict water level distributions across sloping beaches, where the water table in the aquifer may or may not decouple from the free water surface. Five cases are used to test the model for simulating beach water table fluctuations induced by tides, with the model predictions being compared with existing analytical solutions and laboratory and field data published in the literature. The numerical model results show that the integrated model is capable of simulating the combined ground‐ and surface‐water flows in coastal areas. Detailed analysis is undertaken to investigate the capability of the model. Copyright © 2009 John Wiley & Sons, Ltd.     

A semi‐analytical solution for groundwater responses to stream‐stage variations and tidal fluctuations in a coastal aquifer

A two‐dimensional semi‐analytical solution to analyse stream–aquifer interactions in a coastal aquifer where groundwater level responds to tidal effects is presented. The conceptual model considered is a two‐dimensional subsurface system with stream and coastline boundaries at right angles. The dimensional and non‐dimensional boundary value problems were solved for water level in the aquifer by successive application of Laplace and Fourier transform techniques, and the results were obtained by numerical inversion of the transformed solution. The solution was then verified by reducing the solutions to one‐dimensional known problems and comparing the results with those from previous studies. Hypothetical examples were used to examine the characteristics of water‐level variations due to the variations in stream stage and the fluctuations in tide level. Sensitivity analysis indicated that streambed leakance has no influence over the amplitude of groundwater fluctuations, but that the effect of stream stage increases with increasing leakance. Little difference was observed in the water level for different aquifer penetration ratios with narrow stream width. Increases in streambed leakance caused increases in the effect of aquifer penetration by the stream on the water level. An increased specific yield value resulted in decreased amplitude of water fluctuations and mean water level, and showed that water‐level variations due to stream and tidal boundaries are sensitive to specific yield. Copyright © 2006 John Wiley & Sons, Ltd.     

A general solution for tide‐induced groundwater fluctuation in an estuarine‐coastal confined and unconfined aquifer system

This paper presents an analytical solution to tide‐induced head fluctuations in a two‐dimensional estuarine‐coastal aquifer system that consists of an unconfined aquifer and a heterogeneous confined aquifer extending under a tidal river with a semipermeable layer between them. This study considers the joint effects of tidal‐river leakage, inland leakage, dimensionless transmissivity between the tidal‐river and inland confined aquifer, and transmissivity anisotropic ratios. The analytical solution for this model is obtained via the separation of variables method. Three existing solutions related to head fluctuation in one‐ or two‐dimensional leaky confined aquifers are considered as special cases in the present solution. This study shows that there is a threshold of tidal‐river confined aquifer length. When the tidal‐river length is greater than the threshold length, the inland head fluctuations remain sensitive to the leakage effect but become insensitive to the tidal‐river width and dimensionless transmissivity. Considering leakage and transmissivity anisotropy, this study also demonstrates that at a location farther from the river–inland boundary, head fluctuations increase with increasing leakage and transmissivity anisotropy; the maximum head fluctuation occurs when leakage and transmissivity anisotropy are both at their maximum values. The combined action of the 3 effects of loading, tidal‐river aquifer leakage, and inland aquifer leakage differs significantly according to various aquifer parameters. The analytical solution in this paper can be applied to demonstrate the behaviours of the head fluctuations of an estuarine‐coastal aquifer system, and the head fluctuations can be clearly described when the tidal and hydrogeological parameters are derived from field measurement data or hypothetical cases.     

Air and water flows induced by pumping tests in unconfined aquifers with low‐permeability zones

Air flows from the atmosphere into an unconfined aquifer when the water table falls during pumping tests. Pumping test results in unconfined aquifers may be significantly affected by low‐permeability zones (LPZs) near the initial water table position, because they restrict the downward movement of air. A transient, three‐dimensional air–water two‐phase flow model is employed to investigate numerically the effects of local heterogeneity on pumping test results in unconfined aquifers. Two cases of local heterogeneities are considered herein: a LPZ around the pumping well and on one side of the pumping well. Results show that the drawdown with the LPZ is significantly greater than that of the homogeneous aquifer. The differences in drawdown are the most significant at intermediate times and gradually diminish at later times. The LPZ significantly reduces air flow from the atmosphere to the aquifer. The pore air velocity in the LPZ is very low. The air pressure at the observation point under the LPZ when air begins to enter is significantly lower than the air pressure of the homogeneous aquifer at the same point. After that, the air pressure increases quickly and then increases slowly. The time for the air pressure to reach the atmospheric pressure is significantly longer. Copyright © 2013 John Wiley & Sons, Ltd.     

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi‐steady‐state rise in the mean water‐table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed and this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water‐table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright © 2001 John Wiley & Sons, Ltd.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Detecting Drawdowns Masked by Environmental Stresses with Water‐Level Models

Detecting and quantifying small drawdown at observation wells distant from the pumping well greatly expands the characterized aquifer volume. However, this detection is often obscured by water level fluctuations such as barometric and tidal effects. A reliable analytical approach for distinguishing drawdown from nonpumping water‐level fluctuations is presented and tested here. Drawdown is distinguished by analytically simulating all pumping and nonpumping water‐level stresses simultaneously during the period of record. Pumping signals are generated with Theis models, where the pumping schedule is translated into water‐level change with the Theis solution. This approach closely matched drawdowns simulated with a complex three‐dimensional, hypothetical model and reasonably estimated drawdowns from an aquifer test conducted in a complex hydrogeologic system. Pumping‐induced changes generated with a numerical model and analytical Theis model agreed (RMS as low as 0.007 m) in cases where pumping signals traveled more than 1 km across confining units and fault structures. Maximum drawdowns of about 0.05 m were analytically estimated from field investigations where environmental fluctuations approached 0.2 m during the analysis period.     

Analytical model for fully three‐dimensional radial dispersion in a finite‐thickness aquifer

Analytical solutions for contaminant transport in a non‐uniform flow filed are very difficult and relatively rare in subsurface hydrology. The difficulty is because of the fact that velocity vector in the non‐uniform flow field is space‐dependent rather than constant. In this study, an analytical model is presented for describing the three‐dimensional contaminant transport from an area source in a radial flow field which is a simplest case of the non‐uniform flow. The development of the analytical model is achieved by coupling the power series technique, the Laplace transform and the two finite Fourier cosine transform. The developed analytical model is examined by comparing with the Laplace transform finite difference (LTFD) solution. Excellent agreements between the developed analytical model and the numerical model certificate the accuracy of the developed model. The developed model can evaluate solution for Peclet number up to 100. Moreover, the mathematical behaviours of the developed solution are also studied. More specifically, a hypothetical convergent flow tracer test is considered as an illustrative example to demonstrate the three‐dimensional concentration distribution in a radial flow field. The developed model can serve as benchmark to check the more comprehensive three‐dimensional numerical solutions describing non‐uniform flow contaminant transport. Copyright © 2009 John Wiley & Sons, Ltd.     

Underdamped slug tests with unsaturated‐saturated flows by considering effects of wellbore skins

An analytical solution is presented for the slug tests conducted in a partially penetrating well in an unconfined aquifer affected from above by an unsaturated zone. The solution considers the effects of wellbore skin and oscillatory responses on underdamped slug tests. The flow in the saturated zone is described by a two‐dimensional, axially symmetric governing equation, and the flow in the unsaturated zone above the water table by a linearized one‐dimensional Richards' equation. The unsaturated medium properties are represented by the exponential constitutive relationships. A Laplace domain solution is derived using the Laplace and finite Fourier transform and the solution in the real‐time domain is evaluated using the numerical inverse Laplace transform method. The solution derived in this study is more general and reduces to the most commonly used solutions for slug tests in their specified conditions. It is found that the unsaturated flow has a significant impact on the slug test conducted in an unconfined aquifer. The impact of unsaturated flow on such a slug test is enhanced with a larger anisotropy ratio, a shorter well screen length, a shorter distance between the well screen and the water table, or a larger well screen radius. The impact of unsaturated flow on slug tests decreases as the degree of penetration (the length of well screen) increases. For a fixed well screen length, the impact of unsaturated flow on slug tests decreases as the distance between the centre of screen and the water table increases. A large dimensionless well screen radius (>0.01) leads to significant effects of unsaturated flow on slug tests. The unsaturated flow reduces the oscillatory responses to underdamped slug tests. The unsaturated zone has significant impact on slug test under high‐permeability wellbore skin.     

An analytical study on three‐dimensional versus two‐dimensional water table‐induced flow patterns in a Tóthian basin

Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.     

Quasi-Saturated Layer: Implications for Estimating Recharge and Groundwater Modeling

This study presents an extension of the concept of “quasi-saturation” to a quasi-saturated layer, defined as the uppermost dynamic portion of the saturated zone subject to water table fluctuations. Entrapped air here may cause substantial reductions in the hydraulic conductivity (K) and fillable pore water. Air entrapment is caused by a rising water table, usually as a result of groundwater recharge. The most significant effects of entrapped air are recharge overestimation based on methods that use specific yield (S_y), such as the water table fluctuation method (WTF), and reductions in K values. These effects impact estimation of fluid flow velocities and contaminant migration rates in groundwater. In order to quantify actual groundwater recharge rates and the effects of entrapped air, numerical simulations with the FEFLOW (Version 7.0) groundwater flow model were carried out using a quasi-saturated layer for a pilot area in Rio Claro, Brazil. The calculated recharge rate represented 16% of the average precipitation over an 8-year period, approximately half of estimates using the WTF method. Air entrapment amounted to a fillable porosity of 0.07, significant lower that the value of 0.17 obtained experimentally for S_y. Numerical results showed that the entrapped air volume in the quasi-saturated layer can be very significant (0.58 of the air fraction) and hence can significantly affect estimates of groundwater recharge and groundwater flow rates near the water table.     

Discussion about the validity of sharp‐interface models to deal with seawater intrusion in coastal aquifers

Analytical models have been exhaustively used to study simple seawater intrusion problems and the sustainable management of groundwater resources in coastal aquifers because of its simplicity, easy implementation, and low computational cost. Most of these models are based on the sharp‐interface approximation and the Ghyben–Herzberg relation, and their governing equations are expressed in terms of a single potential theory to calculate critical pumping rates in a coastal pumping scenario. The Ghyben–Herzberg approach neglects mixing of fresh water and seawater and implicitly assumes that salt water remains static. Therefore, the results of the analytical solutions may be inaccurate and unacceptable for some real‐complex case studies. This paper provides insight into the validity of sharp‐interface models to deal with seawater intrusion in coastal aquifers, i.e. when they can be applied to obtain accurate enough results. For that purpose, this work compares sharp‐interface solutions, based on the Ghyben–Herzberg approach, with numerical three‐dimensional variable‐density flow simulations for a set of heterogeneous groundwater flow and mass transport parameters, and different scenarios of spatially distributed recharge values and spatial wells placement. The numerical experiment has been carried out in a 3D unconfined synthetic aquifer using the finite difference numerical code SEAWAT for solving the coupled partial differential equations of flow and density‐dependent transport. This paper finds under which situations the sharp‐interface solution gives good predictions in terms of seawater penetration, transition zone width and critical pumping rates. Additionally, the simulation runs indicate to which parameters and scenarios the results are more sensitive. Copyright © 2013 John Wiley & Sons, Ltd.     

Tide‐induced groundwater head fluctuation in a coastal aquifer system with a submarine outcrop covered by a thin silt layer

We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 John Wiley & Sons, Ltd.     

Terrestrial‐originated submarine groundwater discharge through deep multilayered aquifer systems beneath the seafloor

A large quantity of submarine groundwater discharge (SGD) of about 1000 m³ day^?1 m^?1 of the 600‐km‐long shoreline of South Atlantic Bight has been estimated by Moore (Global Biogeochemical Cycles, 2010b, 24, GB4005, doi: 10.1029/2009GB003747 ). However, there is great uncertainty in estimating the percentage of net, land‐originated groundwater recharge of SGD. Moreover, most previous studies considered the homogeneous case for the coastal superficial aquifers. Here, we investigated the terrestrial‐originated SGD through a multilayered submarine aquifer system, which comprises two confined aquifers and two semi‐permeable layers. The inland recharge includes a constant part representing the annual average and a periodical part representing its seasonal variation. An analytical solution was derived and used to analyse the distributions of the terrestrial‐originated SGD from the multilayered aquifers along the Winyah Bay transect, South Atlantic Bight. It is found that the width of the zone of SGD from the upper aquifer ranges from ~0.8 to ~8.0 km depending on the leakance of the seabed semi‐permeable layer. A head of the upper aquifer at a coastline 1.0 m higher than the mean sea level will cause a SGD of 1.82– 18.3 m³ day^?1 m^?1 from that aquifer as the seabed semi‐permeable layer's leakance varies from 0.001 to 0.1 day^?1, providing considerable possibility for considerable land‐originated SGD. Seasonal terrestrial‐originated SGD variations predicted by the analytical model provide consistent explanation of the seasonal variation of ²²⁶Ra observed by Moore (Journal of Geophysics, 2007, 112, C10013, doi: 10.1029/2007JC004199 ). The contribution of the lower aquifer to SGD is only 1.2–12% of that of the upper aquifer. Copyright © 2014 John Wiley & Sons, Ltd.     

Two‐dimensional flow response to tidal fluctuation in a heterogeneous aquifer‐aquitard system

In alluvial coastal aquifers, finer sediments are preferentially deposited along the downstream direction, so the hydraulic conductivity is generally heterogeneous and changes with distance from the coastline. To investigate the influence of aquifer heterogeneity on seawater‐groundwater interaction, a new two‐dimensional model characterising groundwater flow in an aquifer‐aquitard system was developed assuming that the hydraulic conductivity of the aquifer linearly increases with the distance from the coastline along the inland direction. A closed‐form analytical solution was derived using the separation‐of‐variables method. Comparing the new solution with the numerical solution by comsol Multiphysics (Sweden) based on the finite‐element method, one can see that the new solution agreed with the numerical solution very well except at the early time. We found that both aquitard leakance and the heterogeneity factor (b) could result in the propagation bias. The propagation bias represents the inconsistency between the theoretical calculation and the observed strong attenuation and small time lag between the head and tide fluctuations. The attenuation decreased with perpendicular distance from the coastline (x‐axis), whereas the time lag increased with distance along the x‐axis. The relationship between the time lag and the distance along the x‐axis seemed to be linear when b was 0.001 m^?1, whereas it obeyed a power function when b was greater than 0.01 m^?1. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of Salt Water Upconing Using a Steady‐State Solution for Partial Completion of a Pumped Well

A new steady‐state analytical solution to the two‐dimensional radial‐flow equation was developed for drawdown (head) conditions in an aquifer with constant transmissivity, no‐flow conditions at the top and bottom, constant head conditions at a known radial distance, and a partially completed pumping well. The solution was evaluated for accuracy by comparison to numerical simulations using MODFLOW. The solution was then used to estimate the rise of the salt water‐fresh water interface (upconing) that occurs under a pumping well, and to calculate the critical pumping rate at which the interface becomes unstable, allowing salt water to enter the pumping well. The analysis of salt water‐fresh water interface rise assumed no significant effect on upconing by recharge; this assumption was tested and supported using results from a new steady‐state analytical solution developed for recharge under two‐dimensional radial‐flow conditions. The upconing analysis results were evaluated for accuracy by comparison to those from numerical simulations using SEAWAT for salt water‐fresh water interface positions under mild pumping conditions. The results from the equation were also compared with those of a published numerical sharp‐interface model applied to a case on Cape Cod, Massachusetts. This comparison indicates that estimating the interface rise and maximum allowable pumping rate using the analytical method will likely be less conservative than the maximum allowable pumping rate and maximum stable interface rise from a numerical sharp‐interface model.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed‐form approximations are presented for calculating the steady‐state groundwater age distribution in two‐dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi‐infinite and of finite thickness. The approximations are developed on the basis of existing one‐dimensional analytical solutions for travel‐time calculation in fresh water lenses and approximate streamline formulations. The two‐dimensional age distributions based on the closed‐form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed‐form solutions perform well over a range of conditions found in natural systems.     

Changes in Entrapped Gas Content and Hydraulic Conductivity with Pressure

Water table fluctuations continuously introduce entrapped air bubbles into the otherwise saturated capillary fringe and groundwater zone, which reduces the effective (quasi‐saturated) hydraulic conductivity, K_quasi, thus impacting groundwater flow, aquifer recharge and solute and contaminant transport. These entrapped gases will be susceptible to compression or expansion with changes in water pressure, as would be expected with water table (and barometric pressure) fluctuations. Here we undertake laboratory experiments using sand‐packed columns to quantify the effect of water table changes of up to 250 cm on the entrapped gas content and the quasi‐saturated hydraulic conductivity, and discuss our ability to account for these mechanisms in ground water models. Initial entrapped air contents ranged between 0.080 and 0.158, with a corresponding K_quasi ranging between 2 and 6 times lower compared to the K_s value. The application of 250 cm of water pressure caused an 18% to 26% reduction in the entrapped air content, resulting in an increase in K_quasi by 1.16 to 1.57 times compared to its initial (0 cm water pressure) value. The change in entrapped air content measured at pressure step intervals of 50 cm, was essentially linear, and could be modeled according to the ideal gas law. Meanwhile, the changes in K_quasi with compression–expansion of the bubbles because of pressure changes could be adequately captured with several current hydraulic conductivity models.