Impact of Sea-Level Rise on Sea Water Intrusion in Coastal Aquifers

Despite its purported importance, previous studies of the influence of sea-level rise on coastal aquifers have focused on specific sites, and a generalized systematic analysis of the general case of the sea water intrusion response to sea-level rise has not been reported. In this study, a simple conceptual framework is used to provide a first-order assessment of sea water intrusion changes in coastal unconfined aquifers in response to sea-level rise. Two conceptual models are tested: (1) flux-controlled systems, in which ground water discharge to the sea is persistent despite changes in sea level, and (2) head-controlled systems, whereby ground water abstractions or surface features maintain the head condition in the aquifer despite sea-level changes. The conceptualization assumes steady-state conditions, a sharp interface sea water-fresh water transition zone, homogeneous and isotropic aquifer properties, and constant recharge. In the case of constant flux conditions, the upper limit for sea water intrusion due to sea-level rise (up to 1.5 m is tested) is no greater than 50 m for typical values of recharge, hydraulic conductivity, and aquifer depth. This is in striking contrast to the constant head cases, in which the magnitude of salt water toe migration is on the order of hundreds of meters to several kilometers for the same sea-level rise. This study has highlighted the importance of inland boundary conditions on the sea-level rise impact. It identifies combinations of hydrogeologic parameters that control whether large or small salt water toe migration will occur for any given change in a hydrogeologic variable.     

Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.     

Interaction of Aquifer and River‐Canal Network near Well Field

The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.     

Modeling the Transient Response of Saline Intrusion to Rising Sea‐Levels

Sea levels are expected to rise as a result of global temperature increases, one implication of which is the potential exacerbation of sea water intrusion into coastal aquifers. Given that approximately 70% of the world's population resides in coastal regions, it is imperative to understand the interaction between fresh groundwater and sea water intrusion in order to best manage available resources. For this study, controlled investigation has been carried out concerning the temporal variation in sea water intrusion as a result of rising sea levels. A series of fixed inland head two‐dimensional sea water intrusion models were developed with SEAWAT in order to assess the impact of rising sea levels on the transient migration of saline intrusion in coastal aquifers under a range of hydrogeological property conditions. A wide range of responses were observed for typical hydrogeological parameter values. Systems with a high ratio of hydraulic conductivity to recharge and high effective porosity lagged behind the equilibrium sea water toe positions during sea‐level rise, often by many hundreds of meters, and frequently taking several centuries to equilibrate following a cease in sea‐level rise. Systems with a low ratio of hydraulic conductivity to recharge and low effective porosity did not develop such a large degree of disequilibrium and generally stabilized within decades following a cease in sea‐level rise. This study provides qualitative initial estimates for the expected rate of intrusion and predicted degree of disequilibrium generated by sea‐level rise for a range of hydrogeological parameter values.     

Effect of a Freshwater Layer Overlying Salt Water on the Regime of Replacement in Confined Coastal Aquifers

A mathematical model of fresh groundwater flow from a semi-infinite confined aquifer into a sea (pool, trench, and the like) filled with salt water and having a freshwater layer above its horizon. To study this problem a mixed boundary-value problem of the theory of analytical functions is formulated and solved with the use of the Polubarinova-Kochina method. The obtained exact analytical relationships and numerical calculations are used to perform a detailed hydrodynamic analysis of the effect of the freshwater layer and other physical parameters of the model on the character and the extent of displacement.     

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.     

忻州XZ－2号井水位动态气压效应的研究

对忻州ＸＺ－２井气压系数的确定、排除气压干扰的方法进行了探讨，并收到了较好的处理效果。研究表明，该井孔排除气压干扰后的正常水位动态为一条在稳定水位背景上叠加出现的水位固体潮变化曲线。为识别水位异常，找到了判断依据。     

Use of Groundwater Level Fluctuations near an Operating Water Supply Well to Estimate Aquifer Transmissivity

We developed a method to estimate aquifer transmissivity from the hydraulic-head data associated with the normal cyclic operation of a water supply well thus avoiding the need for interrupting the water supply associated with a traditional aquifer test. The method is based on an analytical solution that relates the aquifer's transmissivity to the standard deviation of the hydraulic-head fluctuations in one or more observation wells that are due to the periodic pumping of the production well. We analyzed the resulting analytical solution and demonstrated that when the observation wells are located near the pumping well, the solution has a simple, Dupuit like form. Numerical analysis demonstrates that the analytical solution can also be used for a quasi-periodic pumping of the supply well. Simulation of cyclic pumping in a statistically heterogeneous medium confirms that the method is suitable for analyzing the transmissivity of weakly or moderately heterogeneous aquifers. If only one observation well is available, and the shift in the phase of hydraulic-head oscillations between the pumping well and the observation well is not identifiable. Prior knowledge of aquifer's hydraulic diffusivity is required to obtain the value of the aquifer transmissivity.     

Maximizing Net Extraction Using an Injection‐Extraction Well Pair in a Coastal Aquifer

In this study, we examine the maximum net extraction rate from the novel arrangement of an injection‐extraction well pair in a coastal aquifer, where fresh groundwater is reinjected through the injection well located between the interface toe and extraction well. Complex potential theory is employed to derive a new analytical solution for the maximum net extraction rate and corresponding stagnation‐point locations and recirculation ratio, assuming steady‐state, sharp‐interface conditions. The injection‐extraction well‐pair system outperforms a traditional single extraction well in terms of net extraction rate for a broad range of well placement and pumping rates, which is up to 50% higher for an aquifer with a thickness of 20 m, hydraulic conductivity of 10 m/d, and fresh water influx of 0.24 m²/d. Sensitivity analyses show that for a given fresh water discharge from an inland aquifer, a larger maximum net extraction is expected in cases with a smaller hydraulic conductivity or a smaller aquifer thickness, notwithstanding physical limits to drawdown at the pumping well that are not considered here. For an extraction well with a fixed location, the optimal net extraction rate linearly increases with the distance between the injection well and the sea, and the corresponding injection rate and recirculation ratio also increase. The analytical analysis in this study provides initial guidance for the design of well‐pair systems in coastal aquifers, and is therefore an extension beyond previous applications of analytical solutions of coastal pumping that apply only to extraction or injection wells.     

昌黎井水位、水氡趋势动态特征分析

分析了昌黎井2004年洗井前后水位、水氡的趋势动态特征,认为洗井是造成水位和水氡两者趋势变化不同步的主要原因之一.经过观测环境调查,判定2009年以来水位趋势下降并打破年变规律与附近温泉井抽水有关.     

地下水超采区井水位非前兆异常分析

根据地下水过量开采的现状,结合每一个观测井所处的地质环境以及井区的水文地质条件、观测系统等因素,对超采地下水造成井水位非前兆异常特征的普遍性与个性进行研究。结果表明,这类非前兆异常与含水层状态密切相关,变化幅度大,空间上不集中,时间上不丛集,一般为单个变化,个性特征较明显。     

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.     

Modeling and Mapping High Water Table for a Coastal Region in Florida using Lidar DEM Data

Predicting and mapping high water table elevation in coastal landscapes is critical for both science application projects like inundation risk analysis and engineering projects like pond design and maintenance. Previous studies of water table mapping focused on the application of geostatistical methods, which cannot predict values beyond an observation spatial domain or generate an ideal pattern for regions with sparse measurements. In this study, we evaluated the multiple linear regression (MLR) and support vector machine (SVM) techniques for high water table prediction and mapping using fine spatial resolution lidar-derived Digital Elevation Model (DEM) data, and designed an application protocol of these two techniques for high water table mapping in a coastal landscape where groundwater, tide, and surface water are related. Testing results showed that SVM largely improved the high water table prediction with a mean absolute error (MAE) of 1.22 feet and root mean square error (RMSE) of 2.22 feet compared to the application of the ordinary Kriging method which could not generate a reasonable water table. MLR was also promising with a MAE of around 2 feet and RMSE of around 3 feet. The study suggests that both MLR and SVM are valuable alternatives to estimate high water table elevation in Florida. Fine resolution lidar DEMs are beneficial for high water table prediction and mapping.     

地下水超采区承压井水位变化机理探讨

通过研究河北平原区地下水超采区第四系第三承压含水层水位的变化机理,依据上覆地层的荷载效应、非稳定流抽水中弱透水层的释水和越流过程,阐明地下水超采区内地下水位大幅下降的内在规律,为科学合理地提取地震前兆信息提供科学依据。