FINITE ELEMENT MODELLING FOR LAND DISPLACEMENTS DUE TO PUMPING

Equations of equilibrium (force balance) and flow in multidimensions were coupled in this paper to describe land displacements due to pressure decline in aquifers. A Galerkin finite element model based on these equations was developed. The saturated/unsaturated behaviour and the isotropic/anisotropic properties of permeability and elasticity were considered when the model was formulated. This model was verified by comparing its simulation results with those of known analytical solutions for simplified cases. The simulation of displacements due to pressure decline in unsaturated media was also performed. Those results demonstrated that the choice of boundary ranges for an aquifer with infinite domain may significantly affect the estimated horizontal and vertical displacements. To obtain a good estimation of land displacements, the boundary ranges should be carefully chosen. The displacements occurring in unconfined aquifers are caused by the drop of the water table and the change in body force in the dewatering zone. Simulation results also indicated that the change in body force should be considered once an unconfined aquifer has been pumped. Otherwise, the horizontal and vertical displacements in unconfined aquifers would be overestimated and underestimated, respectively. The behaviour of land displacements due to pumping was shown to be affected by changes in the total stresses in aquifers.     

Hydrochemistry and isotope geochemistry as tools for groundwater hydrodynamic investigation in multilayer aquifers: a case study from Lomellina, Po plain, South-Western Lombardy, Italy

A multicriteria approach in studying hydrodynamics of a multilayer aquifer system has been used in the Lomellina region (Northern Italy). It involves the reconstruction of the hydrogeological framework coupled to the definition of the hydrochemical and isotopic features of the aquifers. A shallow phreatic aquifer, reaching depths of about 60–80 m from the surface, and deeper aquifers containing confined groundwater, were distinguished. Groundwater generally shows mineralisation decreasing with depth; dissolved ions depict calcium-bicarbonate hydrochemical facies and stable isotopes define the recharge mechanisms, the origin of groundwater, and the hydraulic confinement of deep aquifers. The phreatic aquifer is fed by local infiltration and by streams and irrigation channels. Tritium and Carbon-14 groundwater dating indicate long residence times (on the order of thousands of years) for confined aquifers. The confined aquifers show essentially passive hydrodynamic conditions and maintain a higher piezometric level than the phreatic aquifer. This inhibits the possibility of recent water penetrating far below the surface. The hydrogeological setting of the Lomellina region displays features which are common to other sectors of the Po plain. As a consequence, the results of this study, although conducted on a restricted area, are highly illustrative of groundwater hydrodynamics in large sedimentary aquifers.     

长江河口地区地下水咸、淡水界面的动态变化研究

本文对长江河口地区1960～1985年间78口钻井中的各层承压水进行了245个水化学样品测试，分析内容包括K+、Na+、Ca²+、Mg²+、NH4+、Al³+、Cl-、HCO₃-、CO₃₂-、SO₄₂-离子含量以及酸碱性、温度等项目。其中第I层有样品24个；II层有63个；III层有53个；IV层有70个；V层有35个。测试时采用平行双样进行精度控制，水质分析的结果，按阴阳离子毫克当量总和的平衡关系进行检查，超出允许误差及时进行复检。
为了研究各承压含水层咸、淡水界面分布的变化，我们将矿化度作为水质指标，并利用人工神经网络技术进行了计算机模拟。结果表明，本区第四纪地层中各含水层水质的时空分布格局主要受第四纪海平面波动控制；IV、V含水层中的微咸水、咸水的分布可能受上层咸水通过越流污染下层淡水所致。对II、III含水层在二十世纪六十年代和八十年代两阶段水质分布的分析发现，含水层内淡水区在上海市区有向东南和西北方向明显扩展的趋势。造成这种现象的主要原因在于二十世纪七十和八十年代，上海市区地下水开采量逐年压缩，并开展了人工回灌（均为淡水），从而阻止了由于二十世纪六十年代市区过量抽取地下水而引起的地面沉降等问题。     

扬-泰-靖地区地下水系统水力联系与硫酸盐污染特征

本文对长江三角洲扬-泰-靖地区第四系松散层地下水中环境同位素(D、18O、34S)的分布特征进行了分析,旨在揭示大气降水、长江水、潜水及承压水之间的水力联系,辨别地下水中硫酸盐的来源及其污染状况。研究结果表明潜水含水层接受大气降水及长江水的补给,硫酸盐主要为农业污染来源或与海源硫酸盐的混合。承压含水层主要接受大气降水的补给,与潜水含水层及长江之间的水力联系较差,硫酸盐来源不同。在研究区顶部和沿江地段的浅层孔隙承压水中,硫酸盐来源于硫化物的氧化;在东部的深层孔隙承压水中,硫酸盐主要来源于硫酸盐岩的溶解或海源硫酸盐的滞留,基本未受到潜水或地表水中硫酸盐的污染。     

Hydrochemical and isotopic behaviour of a Saharan phreatic aquifer suffering severe natural and anthropic constraints (case of Oued-Souf region, Algeria)

An approach combining the use of water dissolved chemical species and isotopic fingerprints has been used to understand the behavior of a phreatic aquifer and to determine the origin of its different water components. This aquifer is located in the large sedimentary basin of the Great Oriental Erg (Algeria) and overlies two deeper aquifers: the Complexe Terminal (CT) and the Continental Intercalaire (CI). Besides the deterioration of its groundwater quality, its water table has risen during the last 20 years. A water budget surplus between 950 and 2500 l s^?1 was estimated. Down-gradient groundwater evolution (south-north) has shown that the mineralisation increases from 1.23 to 5.20 g l^?1 due to evaporite minerals dissolution. Chemical and isotopic data demonstrated that in addition to rainfall there is a contribution from the CT and CI aquifers. The latter are tritium-free and less mineralized than the phreatic aquifer. Their radiocarbon contents are very low (<10 pmC, percent modern Carbon) (Pleistocene recharge) whereas quite the contrary is observed for the superficial aquifer which exhibits fairly high and variable C-14 activities (50–100 pmC), evidence of recent recharge. On the basis of tritium contents, two groundwater groups were identified for the phreatic aquifer.     

含水层中海水入侵的数学模型

提出了考虑过渡带的适合潜水含水层和承压含水层的海水入侵数学模型,导出了相应的方程和定解条件.经实例验证,模型是合理和可靠的.     

Assessing groundwater availability and the response of the groundwater system to intensive exploitation in the North China Plain by analysis of long-term isotopic tracer data

The use of isotope tracers as a tool for assessing aquifer responses to intensive exploitation is demonstrated and used to attain a better understanding of the sustainability of intensively exploited aquifers in the North China Plain. Eleven well sites were selected that have long-term (years 1985–2014) analysis data of isotopic tracers. The stable isotopes δ¹⁸O and δ²H and hydrochemistry were used to understand the hydrodynamic responses of the aquifer system, including unconfined and confined aquifers, to groundwater abstraction. The time series data of ¹⁴C activity were also used to assess groundwater age, thereby contributing to an understanding of groundwater sustainability and aquifer depletion. Enrichment of the heavy oxygen isotope (¹⁸O) and elevated concentrations of chloride, sulfate, and nitrate were found in groundwater abstracted from the unconfined aquifer, which suggests that intensive exploitation might induce the potential for aquifer contamination. The time series data of ¹⁴C activity showed an increase of groundwater age with exploitation of the confined parts of the aquifer system, which indicates that a larger fraction of old water has been exploited over time, and that the groundwater from the deep aquifer has been mined. The current water demand exceeds the sustainable production capabilities of the aquifer system in the North China Plain. Some measures must be taken to ensure major cuts in groundwater withdrawals from the aquifers after a long period of depletion.     

Optimization of water resources management using SWOT analysis: the case of Zakynthos Island, Ionian Sea, Greece

Zakynthos, an island of 408 km² in the Ionian Sea, is completely dependent on its groundwater resources for fulfilling the demands of the water supplies. The use of groundwater resources has become particularly intensive during the last decades because of the intense urbanization, the tourist development and the irrigated land expansion that took place. The main aquifers are developed in limestones (karstic), sandstones of neogene deposits (confined) and alluvial deposits (phreatic). This paper focuses on the assessment of their hydrogeological characteristics and the groundwater quality. For this investigation, groundwater level measurements, drilling data, pumping tests and chemical analyses of groundwater samples were used. The average annual consumption that is abstracted from the aquifer systems, is 4.9 × 10⁶ m³ year⁻¹. The exploitable groundwater reserves were estimated to be 3.3 × 10⁶ m³ year⁻¹. In the last decades, the total abstractions exceed the natural recharge, due to the tourist development; therefore the aquifer systems are not used safely. The results of chemical analyses showed a deterioration of the groundwater quality. According to the analyses the shallow alluvial aquifer and the confined aquifer are polluted by nitrates at concentrations in excess of 25 mg L⁻¹. High sulphate concentrations might be related to the dissolution of gypsum. Seawater intrusion phenomena are recorded in coastal parts of aquifer systems. The increased Cl⁻ concentrations in karstic aquifer indicate signs of overexploitation. Strengths, weaknesses, opportunities and threats (SWOT) analysis was applied in order to evaluate the SWOT of the groundwater resources. Moreover, some recommendations are made to assist the rational management that aim at improving the sustainability of the groundwater resources of Zakynthos Island.     

吉林中部平原区地下水温动态变化影响因素的关联度分析

以吉林中部平原区为例,分析了孔隙潜水、孔隙承压水以及白垩系孔隙裂隙承压水温度的多年变化特征,并采用灰色关联度分析方法,确定了各含水层地下水温度与地下水埋深、气温、降水量之间的关联性,分析了影响其变化的主要因素。结果表明,对于3组含水层,从20世纪80年代初至2001年的近20年间,地下水温度都有不同程度的上升。对于孔隙潜水和孔隙承压水,地下水温度呈上升变化的本质原因是全球气候逐年变暖,而地下水活动和大气降水的增强起到了促进作用,特别是对孔隙承压水温度的上升作用更加明显;对于白垩系孔隙裂隙承压水,其温度呈上升变化则主要是由该含水层被普遍超采所引起的。     

The integrated impacts of natural processes and human activities on groundwater salinization in the coastal aquifers of Beihai,southern China

Salinization in coastal aquifers is usually related to both seawater intrusion and water–rock interaction. The results of chemical and isotopic methods were combined to identify the origin and processes of groundwater salinization in Daguansha area of Beihai, southern China. The concentrations of the major ions that dominate in seawater (Cl^?, Na⁺, Ca²⁺, Mg²⁺ and SO ₄ ^2– ), as well as the isotopic content and ratios (²H, ¹⁸O, ⁸⁷Sr/⁸⁶Sr and ¹³C), suggest that the salinization occurring in the aquifer of the coastal plain is related to seawater and that the prevailing hydrochemical processes are evaporation, mixing, dissolution and ion exchange. For the unconfined aquifer, groundwater salinization has occurred in an area that is significantly influenced by land-based sea farming. The integrated impacts of seawater intrusion from the Beibuwan Gulf and infiltration of seawater from the culture ponds are identified in the shallowest confined aquifer (I) in the middle of the area (site BBW2). Leakage from this polluted confined aquifer causes the salinization of groundwater in the underlying confined aquifer (II). At the coastal monitoring site (BBW3), confined aquifer I and lower confined aquifer II are heavily contaminated by seawater intrusion. The weak connectivity between the upper aquifers, and the seaward movement of freshwater, prevents saltwater from encroaching the deepest confined aquifer (III). A conceptual model is presented. Above all, understanding of the origin and processes of groundwater salinization will provide essential information for the planning and sustainable management of groundwater resources in this region.     

Hydrogeochemical characteristics of groundwater in coastal phreatic aquifers of Alleppey district,Kerala

Hydrogeochemical characteristics of groundwater in phreatic aquifers of Alleppey district were studied. Factor analysis has been applied to the chemical analysis data of 32 water samples collected from dug wells to extract the principal factors corresponding to the sources of variation in the hydrochemistry. 12 hydrochemical parameters were correlated and statistically examined. Varimax rotation was used to define the factor scores and percentage of variance in the hydrogeochemistry. A four-factor model is extracted and explains over 80.394% of the total groundwater quality variation. Factor-1 has high loading values of Electrical Conductivity (EC), Ca⁺⁺, and Cl⁻, and reflects the signature of saline water. Similarly strong correlation exists between F3 score and pH. The correlation coefficient matrix between EC and Na⁺, Cl⁻, SO₄⁻⁻ is significant. The mineralogy of coastal aquifers and the marine aerosol are playing significant role in the hydrogeochemistry of groundwater in the phreatic aquifer system.     

Environmental impact assessment using FORM and groundwater system reliability concept: case study Jining, China

In this paper, fist-order reliability method (FORM) is used to evaluate the impacts of uncertainties posed by traditional deterministic models on the environment in Jining, China. Because of groundwater contamination in shallow aquifer, and an increase in water demand, the new wells target the confined aquifer with constant pumping rate of 5,000 m³/d. Using Theis equation, the groundwater drawdown is analyzed to determine whether the confined aquifer will be contaminated. Although the piezometric level is higher than the phreatic level by 11.0 m, the risk of drawdown is still 19.49% when the pumping rate of 5,000 m³/d is maintained for 2 years. The deterministic model indicates a drawdown of 8.94 m which is lower than the maximum tolerance drawdown of 11.0 m. The sensitivity and uncertainty analysis reveal that the model result is more sensitive to transmissivity than specific yield, while the reliability analysis offers significant information for the decision makers. This approach exposes and minimizes the risk of undesirable consequences such as groundwater contamination.     

化隆—循化盆地不同类型含水层组高氟地下水的分布及形成过程

天然成因的高氟地下水是世界范围内备受关注的环境问题和饮用水安全问题。前人对高氟地下水的形成过程已开展了大量研究,但是对于高原盆地复杂水文地质条件下不同类型含水层组(第四系松散层含水层、基岩裂隙或岩溶含水层以及新生代古近纪以来的碎屑岩含水层)高氟地下水的分布和形成过程尚不明确。本文以化隆—循化盆地为研究区,通过采集、测试研究区内的各类地下水样品,分析研究区内不同类型含水层中地下水的化学特征及同位素特征。结果表明,高氟地下水(1.007.73 mg/L)主要分布在沿黄河的河谷区域和巴燕低山丘陵区域的泉水和潜水中以及深部的承压水中,在垂向上高氟地下水无明显分布规律。接受黄河水入渗补给的河谷潜水中氟离子浓度较低,补给黄河的河谷潜水中氟离子浓度较高。贫钙富钠的弱碱性苏打型水有利于地下水中氟的富集。泉水和潜水中氟主要来源于萤石的溶解,而承压水中氟除了来源于萤石外,还来源于其他含氟矿物。对于潜水和第四系松散层泉水,蒸发浓缩作用促进了地下水中氟的富集。另外,阴离子竞争吸附作用、阳离子交换吸附作用是泉水(第四系松散层泉水和基岩裂隙泉水)和潜水中氟元素富集的主要原因,而承压水中氟离子浓度受竞争吸附作用影响较大,阳离子交换吸附作用影响较小。研究成果可为化隆—循化盆地低氟地下水的勘查和开发提供科学依据。     

Identifying the recharge sources and age of groundwater in the Songnen Plain (Northeast China) using environmental isotopes

The recharge sources and groundwater age in the Songnen Plain, Northeast China, were confirmed using environmental isotopes. The isotopic signatures of the unconfined aquifers in the southeast elevated plain and the north and west piedmont, cluster along local meteoric water lines (LMWLs) with a slope of about 5. The signature of source water was obtained by the intersection of these LMWLs with the regional meteoric water line (RMWL). This finding provides evidence that the recharge water for these areas originate from the Changbai Mountains and the Low and High Hingan Mountains, respectively. Groundwater in the unconfined aquifer in the low plain yields a LMWL with a slope of 4.4; its nitrate concentration indicates the admixture of irrigation return flow. The δ-values of the unconfined aquifer in the east elevated plain plot along the RMWL, reflecting recharge by local precipitation. The mean residence time of groundwater in these aquifers is less than 50?years. However, the ¹⁴C age of the groundwater in the confined Quaternary aquifer ranges from modern to 19,500?years, and in the Tertiary confined aquifer from 3,100 to 24,900?years. Modern groundwater is mainly recharged to the Quaternary confined aquifer on the piedmont by local precipitation and lateral subsurface flow.     

Application of genetic algorithm technique to inverse modeling of tide–aquifer interaction

Modeling of tide–aquifer interaction plays a vital role in the management of coastal aquifer systems. A novel and robust methodology is presented in this paper for estimating aquifer parameters of coastal aquifers from tide–aquifer interaction data using tide–aquifer interaction model and genetic algorithm (GA). Two stand-alone computer programs were developed to optimize hydraulic diffusivities of unconfined and confined coastal aquifers at multiple sites using GA technique and tide–aquifer interaction model and considering two approaches (‘lumped tidal component approach’ and ‘multi-tidal component approach’). Five sets of real-world tide–aquifer interaction data at two sites of an unconfined aquifer and one set of tide–aquifer interaction data at three sites of a confined aquifer were used to demonstrate the efficacy of the methodology. The analysis of the GA-based inverse modeling results indicated that the ‘multi-tidal component approach’ yields more accurate and reliable hydraulic diffusivities for the unconfined aquifer (RMSE = 0.0129–0.0521 m, NSE = 0.70–0.97, and d1 = 0.91–0.99) as well as for the confined aquifer (RMSE = 0.0204–0.0545 m, NSE = 0.95–0.97, and d1 = 0.99) compared with the ‘lumped tidal component approach’. A comparative evaluation of data-size revealed that the short-duration datasets of the unconfined aquifer provide more reliable estimate of hydraulic diffusivity than the long-duration datasets. Further, it was found that the spring and neap tidal data yield unreasonable values of hydraulic diffusivity with considerably high values of RMSE and very low values of r ², NSE, and d1, thereby suggesting that spring and neap tidal data are not suitable for aquifer parameter estimation. Overall, it is concluded that the GA-based tide–aquifer interaction model following ‘multi-tidal component approach’ is the most efficient tool for estimating aquifer parameters of unconfined and confined aquifers from tide–aquifer interaction data. The developed methodology is also applicable to other coastal basins of the world irrespective of hydrogeological settings.     

Hydrogeochemical investigations and groundwater provinces of the Friuli Venezia Giulia Plain aquifers,northeastern Italy

Water resources are a key factor, particularly for the planning of the sustainable regional development of agriculture, as well as for socio-economic development in general. A hydrochemical investigation was conducted in the Friuli Venezia Giulia aquifer systems to identify groundwater evolution, recharge and extent of pollution. Temperature, pH, electric conductivity, total dissolved solids, alkalinity, total hardness, SAR, Ca²⁺, Na⁺, K⁺, Mg²⁺, Cl⁻, SO₄ ²⁻, NO₃ ⁻, HCO₃ ⁻ _, water quality and type, saturation indexes and the environmental stable isotope δ¹⁸O were determined in 149 sampling stations. The pattern of geochemical and oxygen stable isotope variations suggests that the sub-surface groundwater (from phreatic and shallow confined aquifers) is being recharged by modern precipitations and local river infiltrations. Four hydrogeological provinces have been recognised and mapped in the Friuli Venezia Giulia Plain having similar geochemical signatures. These provinces have different degrees of vulnerability to contamination. The deep confined groundwater samples are significantly less impacted by surface activities; and it appears that these important water resources have very low recharge rates and would, therefore, be severely impacted by overabstraction.     

Using geochemical indicators to investigate groundwater mixing and residence time in the aquifer system of Djeffara of Medenine (southeastern Tunisia)

Stable isotopes (??²H, ??¹⁸O and ??¹³C) and radiocarbon (¹⁴C) have been used in conjunction with chemical data to evaluate recharge mechanisms and groundwater residence time, and to identify inter-aquifer mixing in the Djeffara multi-aquifer in semi-arid southeastern Tunisia. The southern part of this basin, the Djeffara of Medenine aquifer system, is comprised of two main aquifers of Triassic and Miocene sandstone. The Triassic aquifer presents two compartments; the first one (west of the Medenine fault system) is unconfined with a well-defined isotope fingerprint; the second compartment is deeper and confined. Multi-tracer results show groundwater of different origins, ages and salinities, and that tectonic features control groundwater flows. Fresh and brackish groundwater from the unconfined part of the Triassic aquifer was mostly recharged during the Holocene. The recharge rates of this aquifer, inferred by ¹⁴C ages, are variable and could reach 3.5?mm/year. Brackish water of the deep confined part of the Triassic aquifer has stable isotope composition and ¹⁴C content that indicates earlier recharge during late Pleistocene cold periods. Brackish to saline water of the Miocene aquifer presents variable isotope composition. Groundwater flowing through the Medenine fault system is mainly feeding the Miocene aquifer rather than the deep confined part of the Triassic aquifer.     

Groundwater flow patterns in the San Luis Valley,Colorado, USA revisited: an evaluation of solute and isotopic data

Groundwater systems in the San Luis Valley, Colorado, USA have been re-evaluated by an analysis of solute and isotopic data. Existing stream, spring, and groundwater samples have been augmented with 154 solute and isotopic samples. Based on geochemical stratification, three groundwater regimes have been identified within 1,200 m of the surface: unconfined, upper active confined, and lower active confined with maximum TDS concentrations of 35,000, 3,500 and 600 mg/L, respectively. The elevated TDS of northern valley unconfined and upper active confined systems result from mineral dissolution, ion exchange and methanogenesis of organic and evaporate lake sediments deposited in an ancient lake, herein designated as Lake Sipapu. Chemical evolutions along flow paths were modeled with NETPATH. Groundwater ages, and δ¹³C, δ²H and δ¹⁸O compositions and distributions, suggest that mountain front recharge is the principle recharge mechanism for the upper and lower confined aquifers with travel times in the northern valley of more than 20,000 and 30,000 ¹⁴C years, respectively. Southern valley confined aquifer travel times are 5,000 ¹⁴C years or less. The unconfined aquifer contains appreciable modern recharge water and the contribution of confined aquifer water to the unconfined aquifer does not exceed 20%.     

承压含水层井流-盖层弯曲效应的解析理论

指出承压含水层盖层的弯曲变形与开采井周围的径向地下水运动存在相互作用, 而这一效应在传统的井流理论中没有被认识到.通过引入弹性薄板理论, 建立了无越流的承压含水层井流-顶板弯曲效应的解析模型, 同时考虑了含水层和水的压缩性, 结果表明Theis井流方程给出的抽水降深偏小.在此基础上推导了有越流承压含水层井流-盖层弯曲效应的偏微分方程, 求出了解析解, 并与传统理论的结果进行了对比, 表明Hantush-Jacob公式计算的降深也是偏小的.在抽水井附近和抽水初期, 传统理论可能导致显著的相对误差.      

Methodology to evaluate the renewal period of carbonate aquifers: a key tool for their management in arid and semiarid regions,with the example of Becerrero aquifer,Spain

A methodological procedure is proposed for determining the renewal period (RP), which expresses the ratio of total storage to recharge of carbonate aquifers, and it was applied to the overexploited moderate-size Becerrero carbonate aquifer (southern Spain). To this end, geological and subsurface data—time domain electromagnetic (TEM) soundings and borehole logs—were integrated to construct a three-dimensional (3D) geological model of the aquifer. The interconnected porosity was estimated by analyzing 73 rock samples. The resulting 3D geometrical model makes it possible to quantify the fractions of the aquifer having a confined or unconfined behaviour. Based on the total storage capacity (179?·?10⁶–514?·?10⁶ m³) and available aquifer recharge estimation (4.8?·?10⁶–6.4?·?10⁶ m³/year), an RP between 37 and 106 years is obtained. In view of the RP, an exploitation rate slightly lower than the average recharge of the system is recommended, so that the piezometric level will be stable but below the discharge head that is produced through the springs in natural conditions. The proposed methodology to obtain an aquifer RP and the management strategies designed accordingly are of broad interest, especially for carbonate aquifers, which are abundant in arid and semiarid regions.