Numerical Study of Tide-Induced Airflow and Salt–Fresh Water Dynamics in Unconfined Aquifers

This study employed a coupled water-air two-phase flow and salt water transport model to analyze the behaviors of generated airflow in unsaturated zones and the fluctuations of salinity at the salt–fresh water interface in a two-layered unconfined aquifer with a sloping beach surface subjected to tidal oscillations. The simulation results show that as the new dynamic steady state including effects of tidal fluctuations is reached through multiple tidal cycles, the dispersion zone in the lower salt water wedge is broadened because fresh water/salt water therein flows continuously landward or seaward during tidal cycles. The upper salt–fresh water interface exhibits more vulnerable to the tidal fluctuations, and the variation of salinity therein is periodic, which is irrelevant to the hydraulic head but is influenced by the direction and velocity of surrounding water-flow. With the tidal level fluctuating, airflow is mainly concentrated in the lower permeable layer due to the restraint of the upper semi-permeable layer, and the time-lag between the pore-air pressure and the tidal level increases with distance from the coastline. The effect of airflow in unsaturated zones can be transmitted downward, causing both the magnitude of salinity and its amplitude in the upper salt–fresh water interface to be smaller for the case with airflow than without airflow due to the resistance of airflow to water-flow. Sensitivity analysis reveal that distributions of airflow in unsaturated zones are affected by the permeability of the upper/lower layer and the van Genuchten parameter of the lower layer, not by the van Genuchten parameter of the upper layer, whereas the salinity fluctuations in the salt–fresh water interface are affected only by soil parameters of the lower layer.     

Modeling Coastal Salinity in Quasi 2D and 3D Using a DUALEM‐421 and Inversion Software

Rising sea levels, owing to climate change, are a threat to fresh water coastal aquifers. This is because saline intrusions are caused by increases and intensification of medium‐large scale influences including sea level rise, wave climate, tidal cycles, and shifts in beach morphology. Methods are therefore required to understand the dynamics of these interactions. While traditional borehole and galvanic contact resistivity (GCR) techniques have been successful they are time‐consuming. Alternatively, frequency‐domain electromagnetic (FEM) induction is potentially useful as physical contact with the ground is not required. A DUALEM‐421 and EM4Soil inversion software package are used to develop a quasi two‐ (2D) and quasi three‐dimensional (3D) electromagnetic conductivity images (EMCI) across Long Reef Beach located north of Sydney Harbour, New South Wales, Australia. The quasi 2D models discern: the dry sand (<10 mS/m) associated with the incipient dune; sand with fresh water (10 to 20 mS/m); mixing of fresh and saline water (20 to 500 mS/m), and; saline sand of varying moisture (more than 500 mS/m). The quasi 3D EMCIs generated for low and high tides suggest that daily tidal cycles do not have a significant effect on local groundwater salinity. Instead, the saline intrusion is most likely influenced by medium‐large scale drivers including local wave climate and morphology along this wave‐dominated beach. Further research is required to elucidate the influence of spring‐neap tidal cycles, contrasting beach morphological states and sea level rise.     

鲁14井数字化水位潮汐响应特性初析

以井—含水层系统潮汐理论为基础,利用Baytap-G潮汐分析程序计算鲁14井数字化水位的潮汐响应特征参数,提取井水位潮汐因子和相位差,分析远场大震前后井水位潮汐响应特征参数的变化情况,同时对固体潮和气压影响量进行分离,以期深入了解井水位微动态变化特征,为该井水位观测资料的动态评价和效能评估提供参考。结果表明：鲁14井潮汐地下水流类型以径向流为主,潮汐因子和相位差变化相对稳定,但受远场大震影响显著,其含水层渗透性随潮汐参数的增大（减小）而增大（减小）;地震波作用可使含水层渗透性发生变化,但并非导致含水层渗透性变化的主要因素。     

Rapid recharge of fresh water to the halite‐hosted brine aquifer of Salar de Atacama,Chile

We document, analyse, and interpret direct and rapid infiltration of precipitation to the southern margin of the Salar de Atacama halite‐hosted brine aquifer during two intense precipitation events in 2012–2013. We present physical, geochemical, and stable and radioactive isotope data to detail this influx of water. The two events differ distinctly in the mechanisms of recharge. The 2012 event did not produce direct precipitation onto the salar surface, while the 2013 event did. Both events are recorded by abrupt changes in head in observation wells along the halite aquifer margin. Spatially distributed water levels rose by over 30 cm during the larger 2013 event consistent with remotely sensed observations of surface water extent. The lithium concentration and stable isotopic composition of water indicate dilution of brine and dissolution of salt with fresh water. Tritium measurements of precipitation, surface water, and groundwater all indicate modern influx of water to the halite aquifer along the southern margin. We extend these observations by examining the response of the halite aquifer as a whole to precipitation events during the period of 2000–2010. This study suggests that local recharge to the aquifer during sporadic precipitation onto the halite nucleus is an important component of the modern water budget in this hyper‐arid environment. The rapid dissolution and salinization along the southern margin of the salar halite nucleus are aided by such precipitation events contributing a modern fresh water component to the water budget of the economically valuable lithium‐rich brine. Copyright © 2016 John Wiley & Sons, Ltd.     

Forest fire effects on groundwater in a coastal aquifer (Ravenna,Italy)

We examined the fire‐induced changes in groundwater recharge rate. This aspect is particularly important in the case of large forested areas growing over a coastal aquifer affected by saltwater intrusion. In the Ravenna coastal area (Italy), pine forests grow on coastal dune belts, overlying a sandy unconfined aquifer, which is strongly affected by marine ingression. Three groundwater profiles across the forest and perpendicular to the coastline were monitored for groundwater level, physical, and chemical parameters. The aims were to define groundwater quality, recharge rate, freshwater volume, and highlight change, which occurred after a forest fire with reference to pre‐fire conditions. Analytical solutions based on Darcy Law and the Dupuit Equation were applied to calculate unconfined flow and compare recharge rates among the profiles. The estimated recharge rates increased in the partially and completely burnt areas (219 and 511 mm year^?1, respectively) compared with the pristine pine forest area (73 mm year^?1). Although pre‐fire conditions were similar in all monitored profiles, a post‐fire decrease in salinity was observed across the burnt forest, along with an increase in infiltration and freshwater lens thickness. This was attributed to decrease canopy interception and evapotranspiration caused by vegetation absence after the fire. This research provided an example of positive forest fire feedback on the quantity and quality of fresh groundwater resources in a lowland coastal aquifer affected by saltwater intrusion, with limited availability of freshwater resources. The fire provided an opportunity to evaluate a new forest management approach and consider the restoration and promotion of native dune herbaceous vegetation.     

地下水潮汐现象的物理机制和统一数学方程

用新的分层承压含水层模式 ,不但考虑含水层的力学压缩性质 ,而且考虑含水层的渗流特性 ,并结合扰动信息源的频率特性 ,分别研究扰动源地球固体潮、大气潮和地表负荷潮对承压井水位和流量的影响机理 ,给出相应的偏微分方程。从方程的解释或数值解讨论扰动源与承压井含水层的力学压缩参数、渗流特性参数及与频率特性频数的关系 ,进而给出承压井水位和流量对地球固体潮、大气潮和地表负荷潮汐响应的统一数学方程及其潮汐响应函数 ,并揭示了上述几类潮汐扰动信息源对承压井水位和流量影响的物理机理     

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.     

Submarine Ground Water Discharge Driven by Tidal Pumping in a Heterogeneous Aquifer

Submarine ground water discharge (SGD) is now recognized as an important water pathway between land and sea. It is difficult to quantitatively predict SGD owing to its significant spatial and temporal variability. This study focuses on quantitative estimation of SGD caused by tidally induced sea water recirculation and a terrestrial hydraulic gradient. A two-dimensional hydrogeological model was developed to simulate SGD from a coastal unconfined aquifer in the northeastern Gulf of Mexico, where previous SGD studies were performed. A density-variable numerical code, SEAWAT2000, was applied to simulate SGD. To accurately predict discharge, various influencing factors such as heterogeneity in conductivity, uncertain boundary conditions, and tidal pumping were systematically assessed. The tidally influenced sea water recirculation zone and the fresh water–salt water mixing zone under various tidal patterns, tidal ranges, and water table heights were also investigated. The model was calibrated and validated from long-term, intensive measurements at the study site. The percentage of fresh SGD relative to total SGD ranged from 4% to 50% under normal conditions. Based on simulations of two field measurements in summer and spring, respectively, the fresh water ratios were 9% and 15%, respectively. These results support the hypothesis that the SGD induced by tidally driven sea water recirculation is much larger than terrestrial fresh ground water discharge at this site. The estimates of total and fresh SGD are at the low and high ends, respectively, of the estimation ranges obtained from geochemical tracers (e.g., ²²²Rn).     

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.     

云南曲靖井水位潮汐动态特征分析

以井-含水层系统潮汐理论为基础, 结合云南曲靖井水位观测资料实际情况, 利用Baytap-G提供的潮汐分析方法, 分别计算了曲靖井水位M2、 S2、 O1、 S1K1和M3共5个潮汐分波的振幅比、 相位和计算误差, 分析了井水位观测资料不同潮汐波振幅比和相位的动态变化特征, 认为能反映含水层参数的井水位潮汐振幅比和相位变化与附近的抽水和远距离大震引起的同震响应有关, 具体表现为, 2001年11月和2003年8月的振幅比和相位下降变化可能与抽水有关, 而几次大震后的井水位潮汐振幅比和相位小幅度上升和缓慢恢复过程, 则可能与中远场大震的地震波动力引起的含水层渗透系数增加有关．      

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.     

Mapping Coastal Aquifers by Joint Inversion of DC and TEM Soundings-Three Case Histories

Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

A New Operational Paradigm for Small‐Scale ASR in Saline Aquifers

A new operational paradigm is presented for small‐scale aquifer storage and recovery systems (ASR) in saline aquifers. Regular ASR is often not feasible for small‐scale storage in saline aquifers because fresh water floats to the top of the aquifer where it is unrecoverable. In the new paradigm, fresh water storage is combined with salt water extraction from below the fresh water cone. The salt water extraction counteracts the buoyancy due to the density difference between fresh water and salt water, thus preventing the fresh water from floating up. The proposed approach is applied to assess the feasibility of ASR for the seasonal storage of fresh water produced by desalination plants in tourist resorts along the Egyptian Red Sea coast. In these situations, the continuous extraction of salt water can be used for desalination purposes. An analytical Dupuit solution is presented for the steady flow of salt water toward a well with a volume of fresh water floating on top of the cone of depression. The required salt water discharge for the storage of a given volume of fresh water can be computed with the analytical solution. Numerical modeling is applied to determine how the stored fresh water can be recovered. Three recovery approaches are examined. Fresh water recovery rates on the order of 70% are achievable when salt water is extracted in high volumes, subsurface impermeable barriers are constructed at a distance from the well, or several fresh water recovery drains are used. The effect of ambient flow and interruptions of salt water pumping on the recovery efficiency are reported.     

Locating the zone of saline intrusion in a coastal karst aquifer using springflow data

Coastal fresh water aquifers are an increasingly desirable resource. In a karstic aquifer, sea water intrusion occurs as a salt water wedge, like in porous media. However, preferential flow conduits may alter the spatial and temporal distribution of the salt water. This is typically the case when the outlet of the aquifer is a brackish spring. This paper shows that salinity and flow rate variations at a spring, where salinity is inversely proportional to discharge, can help to understand the hydrodynamic functioning of the aquifer and to locate the fresh water-sea water mixing zone deep inside the aquifer. The volume of water-filled conduit between the sea water intrusion zone and the spring outlet is calculated by the integral over time of the flow rate during the time lag between the flow rate increase and the salinity decrease as measured at the spring. In the example of the spring at Almyros of Heraklio (Crete, Greece), this time lag is variable, depending on the discharge, but the volume of water-filled conduit appears to be constant, which shows that the processes of salt water intrusion and mixing in the conduit are constant throughout the year. The distance between the spring and the zone where sea water enters the conduit is estimated and provides an indication of the position where only fresh water is present in the conduit.     

Rapid nutrient load reduction during infiltration of managed aquifer recharge in an agricultural groundwater basin: Pajaro Valley,California

Artificial recharge of groundwater is an increasingly important method for augmenting groundwater supply and can have a positive or negative influence on the quality of water resources. We instrumented a managed aquifer recharge (MAR) pond in central coastal California to assess how patterns of infiltration and recharge affect the load of nitrate delivered to the underlying aquifer. The concentration of nitrate in infiltrating water consistently decreased during passage through the first metre of subsurface soils. Enrichment of ¹⁸O and ¹⁵ N in the residual nitrate in infiltrating water proceeded in a ratio of 1:2, indicating that denitrification plays a significant role in the quantitative reduction of nutrients exported during infiltration through shallow soils. The extent and rate of nitrate removal was spatially and temporally variable across the bottom of the recharge pond, with 30% to 60% of the nitrate load being removed over the first 6 weeks of managed aquifer recharge operation. During the period of highest N loading to the system, when the average infiltration rate was > 1 m/day, the recharge pond achieved a load reduction efficiency of 7 kg NO₃^?‐N/day/ha, which compares favourably to nitrate load reductions achieved by treatment wetlands. Groundwater mounding and water composition below the recharge pond suggest that recharge and subsequent lateral transport occur heterogeneously in the underlying aquifer. Nitrate concentrations in the aquifer following infiltration were lowered primarily by dilution, with little evidence for additional denitrification occurring in the aquifer in comparison to high rates documented during shallow infiltration. Copyright © 2011 John Wiley & Sons, Ltd.     

华蓥山地区4口井水位潮汐动态特征研究

本文选择沿华蓥山断裂带分布的荣昌等4口观测井,利用Baytap-G潮汐分析方法,计算各井水位和气压及理论固体潮的潮汐振幅谱,比较其潮汐频谱差异,通过对主要潮汐分波振幅的回归计算定量分析各井水位受气压潮和固体潮影响的大小。基于对井水位正常动态的认识,选择各井水位潮汐的主要分波,对井水位长时序数据进行分析计算,提取水位潮汐响应特征参数(振幅比和相位差),进而探讨特征参数动态变化特征。最后对井水位受气压潮和固体潮影响的差异原因进行了初步探讨。结果表明,荣昌井水位主要受气压作用的影响,北碚、大足、南溪三口井水位受固体潮-气压潮综合作用的影响,而荣昌井水位只受气压潮影响可能与该井所处含水层裂隙发育且该井未下设止水套管有关;荣昌井P_1S_1K_1波和南溪井M_2波振幅比和相位差在几次大震后没有明显变化,说明地震波没有使井孔与含水层之间的水流交换发生显著变化,而北碚井和大足井M_2波振幅比和相位差分别在汶川和芦山地震时发生变化,反映了地震波的疏通影响。     

Electromagnetic‐Induction Logging to Monitor Changing Chloride Concentrations

Water from the San Joaquin Delta, having chloride concentrations up to 3590 mg/L, has intruded fresh water aquifers underlying Stockton, California. Changes in chloride concentrations at depth within these aquifers were evaluated using sequential electromagnetic (EM) induction logs collected during 2004 through 2007 at seven multiple‐well sites as deep as 268 m. Sequential EM logging is useful for identifying changes in groundwater quality through polyvinyl chloride‐cased wells in intervals not screened by wells. These unscreened intervals represent more than 90% of the aquifer at the sites studied. Sequential EM logging suggested degrading groundwater quality in numerous thin intervals, typically between 1 and 7 m in thickness, especially in the northern part of the study area. Some of these intervals were unscreened by wells, and would not have been identified by traditional groundwater sample collection. Sequential logging also identified intervals with improving water quality—possibly due to groundwater management practices that have limited pumping and promoted artificial recharge. EM resistivity was correlated with chloride concentrations in sampled wells and in water from core material. Natural gamma log data were used to account for the effect of aquifer lithology on EM resistivity. Results of this study show that a sequential EM logging is useful for identifying and monitoring the movement of high‐chloride water, having lower salinities and chloride concentrations than sea water, in aquifer intervals not screened by wells, and that increases in chloride in water from wells in the area are consistent with high‐chloride water originating from the San Joaquin Delta rather than from the underlying saline aquifer.     

Chemical water types and their distribution in space and time in the Amsterdam dune-water catchment area with artificial recharge

The geology, hydrology and types of groundwater exploitation of the Amsterdam dune-water catchment area with artificial recharge of water from the river Rhine are briefly introduced. A classification of natural and artificial watertypes in the catchment is presented together with their genesis and distribution in space and time. The natural, vertical, hydrochemical zonation in fresh, brackish and salt groundwaters has been disturbed by decades of overdraft and the subsequent start of artificial recharge since 1957. Various mixed-water types resulted. A detailed case study of a part of the catchment area (around the Barnaart-Schuster Canal) covers a period of about 50 years. The decisive influence of geological inhomogeneities on the patterns of saltwater encroachment during the period of overdraft and on the patterns of infiltration of river water into the deep aquifer after the start of artificial recharge are demonstrated. The shifts in the subsurface positions of the various water types due to changes in the groundwater management of the catchment area are elucidated. They give an example of displacement and dispersion phenomena under field conditions.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 1. Small catchment water balance model

A long-term water balance model has been developed to predict the hydrological effects of land-use change (especially forest clearing) in small experimental catchments in the south-west of Western Australia. This small catchment model has been used as the building block for the development of a large catchment-scale model, and has also formed the basis for a coupled water and salt balance model, developed to predict the changes in stream salinity resulting from land-use and climate change. The application of the coupled salt and water balance model to predict stream salinities in two small experimental catchments, and the application of the large catchment-scale model to predict changes in water yield in a medium-sized catchment that is being mined for bauxite, are presented in Parts 2 and 3, respectively, of this series of papers. The small catchment model has been designed as a simple, robust, conceptually based model of the basic daily water balance fluxes in forested catchments. The responses of the catchment to rainfall and pan evaporation are conceptualized in terms of three interdependent subsurface stores A, B and F. Store A depicts a near-stream perched aquifer system; B represents a deeper, permanent groundwater system; and F is an intermediate, unsaturated infiltration store. The responses of these stores are characterized by a set of constitutive relations which involves a number of conceptual parameters. These parameters are estimated by calibration by comparing observed and predicted runoff. The model has performed very well in simulations carried out on Salmon and Wights, two small experimental catchments in the Collie River basin in south-west Western Australia. The results from the application of the model to these small catchments are presented in this paper.