Numerical analysis and evaluation of groundwater recession in a flood detention basin

This study analyzes the groundwater environment in the Yangzhuang flood detention basin in Henan Province, China. A numerical model of groundwater flow is established based on the hydrogeological conditions in the basin and changes of groundwater level in the flood detention basin under flood detention and recession conditions. The results show that during flood diversion and storage, the groundwater level in the basin rises, mainly in four flood detention zones, with a maximum rise of 1.0 m. After the floodwater recedes, the groundwater level slowly drops in the detention basin along with flood discharge, finally returning to its original level fifty days later. This study indicates that groundwater recession in the flood detention basin is a slow process, where the rise of groundwater level may cause environmental problems such as soil swamping.     

Numerical simulation of groundwater flow in regional rock aquifers, southwestern Quebec, Canada

The St. Lawrence Lowlands platform, Quebec, Canada, is a densely-populated area, heavily dependent on groundwater resources. In 1999, the Geological Survey of Canada initiated a large-scale hydrogeological assessment study over a 1,500 km² region northwest of Montreal. The objectives were to define the regional groundwater flow, and to give quantitative estimates of the groundwater dynamic parameters and of the available groundwater resources. The applied approach consisted of defining the hydrogeologic framework, hydraulic properties of the aquifer units, and groundwater dynamic components. Lower Paleozoic sedimentary rocks represent regional aquifer units. Coarse Quaternary fluvio-glacial sediments locally overlay the rock sequence and constitute an interface aquifer unit. Fine marine sediments confine most of the regional aquifers. Collected GIS based information was synthesized in a finite element numerical model. The regional saturated steady-state flow was calibrated under current stress conditions assuming an equivalent porous medium approach. Water budget calculations show that the total groundwater flow in regional aquifers amounts to 97.7 Mm³/y. Infiltration from precipitation provides 86.6% of the groundwater supply, while 9.6% comes from subsurface inflow and the remaining 3.8% is induced recharge from surface waters. Discharge from regional aquifers occurs through flow to streams (76.9%), groundwater withdrawal (18.4%), and underground outflow (4.7%).

---

Resumen La plataforma de Tierras Bajas San Lorenzo, Quebec, Canadá, es un área densamente poblada que depende fuertemente de recursos de agua subterránea. En 1999 el Servicio Geológico de Canadá inició un estudio de evaluación hidrogeológica a gran escala sobre un área de 1,500 km² en la región noroeste de Montreal. Los objetivos fueron definir el flujo regional de agua subterránea y aportar estimados cuantitativos de los parámetros dinámicos de agua subterráne y de los recursos disponibles de agua subterránea. El enfoque aplicado consistió en definir el marco hidrogeológico, propiedades hidráulicas de las unidades acuíferas, y los componentes dinámicos de agua subterránea. Rocas sedimentarias del Paleozoico Inferior representan unidades regionales de acuíferos. Sedimentos fluvio-glaciares Cuaternarios gruesos sobreyacen localmente la secuencia rocosa y constituyen una unidad acuífera interfacial. Sedimentos marinos finos confinan la mayoría de acuíferos regionales. Información colectada basada en SIG se sintetizó en un modelo numérico de elemento finito. El flujo regional saturado en régimen permanente se calibró bajo condiciones de stress asumiendo un enfoque de medio poroso equivalente. Los cálculos de balance hídrico muestran que el flujo total de agua subterránea en acuíferos regional alcanza 97.7 Mm³/año. Infiltración a partir de lluvia aporta el 86.6% del abastecimiento al agua subterránea, mientras que el 9.6% proviene de entradas subsuperficiales y el restante 3.8% consiste de recarga inducida a partir de aguas superficiales. La descarga proveniente de acuíferos regionales ocurre a través de flujo a ríos (76.9%), utilización de agua subterránea (18.4%), y salida subterránea (4.7%).

---

Résumé La plateforme Lowllands du Saint-Laurent, Québec, Canada, est une aire densément peuplée, dépendant grandement des ressources en eau souterraine. En 1999, le Service Géologique du Canada a initié une étude hydrogéologique à grande-échelle sur 1500 km² au Nord-Ouest de Montréal. Le sobjectifs ont été de définir la dynamique de lécoulement régional, et de donner des estimations quantitatives des paramètres dynamiques des ressources disponibles en eau souterraine. Lapproche utilisée consista à définir le cadre hydrogéologique de travail, les propriétés hydrauliques des unités aquifères, et les composantes dynamiques des eaux souterraines. Les roches sédimentaires du Paléozoïque Inférieur représentent les unités aquifères régionales. Les sédiments marins fins confinent la plus grande partie des aquifères régionaux. Les informations de base, collectées dans un SIG, ont été synthétisées dans un modèle numérique aux éléments finis. Lécoulement permanent régional, en zone saturée, a été calibré en conditions de stress en assumant une approche de milieu poreu équivalent. Les calculs du bilan hydrologique montrent que lécoulement total des eaux souterraines dans les aquifères régionaux atteind 97,7 Mm³/an. Linfiltration à partir des précipitations apporte 86.6% de leau souterraine exploitée; sachant que 9,6% proviennent découlement de subsurfaces, et que les 3,8% restants proviennent de la recharge via les eaux de surface. Le débit pompé des aquifères régionaux apparaît à travers lécoulement des cours deau, le rabattement des eaux souterraines (18.4%), et lécoulement ascendant (4.7%).

     

Radon in groundwater as a tracer to assess flow velocities: two test cases from Israel

 Two test cases from Israel are presented herein employing the decay rate of radon along the flow path to assess groundwater flow velocities. Groundwater flow reaching the fault zone emerges in several places along the rift fault zone as thermal springs because of deep water confinement. The high water temperature of the surface is indicative of high vertical flow velocities, which maintains the original high temperatures. Knowing the Rn content at a source point and at a given down-gradient, and assuming no Rn addition from the water itself or along the flow path, one can calculate the flow velocity based on the Rn half-life time. The decay of Rn in western Galilee was found to be ∼570–150 pCi/l, and in the Dead Sea area from 5000–2000 pCi/l along a respective flow path of 1000 and 200 m, Based on the above, the calculated flow velocities were compared with those obtained from pumping tests in the study area. The method is applicable, because of the short Rn half-life, to cases of high Rn contents, short distances and high flow velocities. Received: 18 January 2000 · Accepted: 21 March 2000     

Estimating groundwater storage changes in the Mississippi River basin (USA) using GRACE

Based on satellite observations of Earth’s time variable gravity field from the Gravity Recovery and Climate Experiment (GRACE), it is possible to derive variations in terrestrial water storage, which includes groundwater, soil moisture, and snow. Given auxiliary information on the latter two, one can estimate groundwater storage variations. GRACE may be the only hope for groundwater depletion assessments in data-poor regions of the world. In this study, soil moisture and snow were simulated by the Global Land Data Assimilation System (GLDAS) and used to isolate groundwater storage anomalies from GRACE water storage data for the Mississippi River basin and its four major sub-basins. Results were evaluated using water level records from 58 wells set in the unconfined aquifers of the basin. Uncertainty in the technique was also assessed. The GRACE-GLDAS estimates compared favorably with the well based time series for the Mississippi River basin and the two sub-basins that are larger than 900,000 km². The technique performed poorly for the two sub-basins that have areas of approximately 500,000 km². Continuing enhancement of the GRACE processing methods is likely to improve the skill of the technique in the future, while also increasing the temporal resolution.     

格尔木河流域地下水数值模拟

格尔木河流域是西北干旱区水资源开发利用的典型地区,为满足可持续开发利用的需要,针对《青海省格尔木河流域水资源综合利用规划》提出的水资源开发利用需求量,分析其是否对下游盐湖生态环境产生破坏影响。根据流域概况建立格尔木河流域水资源数学模拟模型,对平均水文序列变化趋势及各种典型年的水资源状况进行预测,探讨在不影响下游盐湖生态环境的前提下,格尔木河流域水资源合理开发利用的方案。     

基于MODFLOW-CFP的贵州大井流域岩溶地下水数值模拟

贵州大井流域岩溶分布广泛,岩溶水是当地人民生产和生活的主要来源。由于对岩溶水资源的不合理开发利用,水资源短缺现象经常发生。大井流域水文地质条件复杂,管道-多孔介质双重介质特征明显。文章采用MODFLOW-CFP耦合模型对大井流域展开数值模拟,进而掌握大井流域地下水运动规律、准确评价岩溶水资源,促进其合理开发利用。结果表明:大井流域管道与多孔介质交换量为6 719.1 m³·a^-1,主要集中在上游和中游;总补给水量为10 977.3×10⁴m³·a^-1,补给模数为133.495m³·km^-2·a^-1,其中降雨汇入量和降雨入渗量占总补给量的81.35%,而总排泄量为10 813.47×10⁴m³·a^-1,主要在地下河出口排泄。     

Numerical simulation of groundwater flow in the Peshawar intermontane basin, northwest Himalayas

The hypothesis that abnormal fluid pressure is generated in basins under tectonic compression is tested. The study site, between the Main Karakoram Thrust (MKT), Main Mantle Thrust (MMT), Main Central Thrust (MCT), Main Boundary Thrust (MBT) and Salt Range Thrust (SRT) in northwest Pakistan, is experiencing a tectonic compression of 90 MPa. The Peshawar basin is a broad, oval depression comprising a thick sequence of lacustrine, deltaic and fluvial sediments overlain by loess and alluvial deposits dated at 2.8–0.6 Ma. It is surrounded by the Precambrian and Tertiary intrusive and metamorphic rocks on the north and sedimentary rocks of Paleogene and Neogene to the south. The basin was divided into four hydrostratigraphic units for numerical simulations using the three-dimensional finite-element model FEMWATER within groundwater modeling system (GMS) ver. 5.1. Simulated pressure head data have been compared with the field measurements of hydraulic heads. Transient simulations indicate that topography alone is not sufficient to induce the pressure heads observed in the field, generating consistently positive residuals 0.98–2.90 m over the topography-driven flow. The residuals disappeared after inclusion of the elastic properties of the four hydrostratigraphic units in the model, suggesting the effect of tectonic compression.     

Tracing stable isotope values from meteoric water to groundwater in the southwestern part of the Chad basin

The southwestern Chad basin is a semi-arid region with annual rainfall that is generally less than 500 mm and over 2,000 mm of evapotranspiration. Surface water in rivers is seasonal, and therefore groundwater is the perennial source of water supply for domestic and other purposes. Stable isotope has been measured for rainwater, surface water and groundwater samples in this region. The stable isotope data have been used to understand the inter-relationships between the rainwater, surface water, shallow and deep groundwater of this region. This is being used in a qualitative sense to demonstrate present day recharge to the groundwater. Stable isotope in rainwater for the region has an average value of –4‰ δ¹⁸O and –20‰ δ²H. Surface water samples from rivers and Lake Chad fall on the evaporation line of this average value. The Upper Zone aquifer water samples show stable isotope signal with a wide range of values indicating the complex character of the aquifer Zone with three distinguishable units. The wide range of values is attributable to waters from individual unit and/or mixture of waters of different units. The Middle and Lower aquifers Zones’ waters show similar stable isotopes values, probably indicating similarity in timing and/or mechanism of recharge. These are palaeowaters probably recharged under a climate that is different from today. The Upper Zone aquifer is presently being recharged as some of its waters show stable isotope compositions similar to those of average rainfall waters of the region.     

双Packer系统在鄂尔多斯盆地地下水勘查中的应用

双Packer系统在鄂尔多斯盆地地下水勘查中的应用,为在该盆地开展深层地下水水化学同位素的研究,提供了一个有效的工具。了解该系统的工作原理、应用条件和在鄂尔多斯盆地地下水勘查中的具体应用,对在我国水文界普及该技术具有重要意义。     

Regional hydrostratigraphy and groundwater flow modeling in the arsenic-affected areas of the western Bengal basin, West Bengal, India

The first documented interpretation of the regional-scale hydrostratigraphy and groundwater flow is presented for a ～21,000-km² area of the arsenic-affected districts of West Bengal [Murshidabad, Nadia, North 24 Parganas and South 24 Parganas (including Calcutta)], India. A hydrostratigraphic model demonstrates the presence of a continuous, semi-confined sand aquifer underlain by a thick clay aquitard. The aquifer thickens toward the east and south. In the south, discontinuous clay layers locally divide the near-surface aquifer into several deeper, laterally connected, confined aquifers. Eight 22-layer model scenarios of regional groundwater flow were developed based on the observed topography, seasonal conditions, and inferred hydrostratigraphy. The models suggest the existence of seasonally variable, regional, north–south flow across the basin prior to the onset of extensive pumping in the 1970s. Pumping has severely distorted the flow pattern, inducing high vertical hydraulic gradients across wide cones of depression. Pumping has also increased total recharge (including irrigational return flow), inflow from rivers, and sea water intrusion. Consequently, downward flow of arsenic contaminated shallow groundwater appears to have resulted in contamination of previously safe aquifers by a combination of mechanical mixing and changes in chemical equilibrium.     

基于环境同位素分析吐鲁番盆地地下水流系统

吐鲁番盆地是中国重要的煤炭基地之一,研究盆地地下水流系统对于区域水资源评价具有重要意义。通过在野外开展取样工作,首先根据取样点井深的不同,将所有样点分为4组,结合当地水文地质条件、同位素值的高程效应以及取样点位置,将同位素值以及位置相近的点归为同一个水流系统。经分析,在火焰山以北存在一个相同型式的地下水流动系统,而在火焰山以南随着山区坎尔其水库的修建,水利工程设施的不断完善和农业灌溉制度的改革,存在多级地下水流动系统,其形成条件主要受自然和人类活动影响。其中同位素和TDS值都较高的T1点曾是艾丁湖湖面的边缘,当时湖水水位浅,同位素蒸发效应明显。同时,火焰山地下基岩存在缺口,存在北盆地向南盆地的地下水补给。     

Numerical simulation of groundwater flow for a coastal plain in Japan: data collection and model calibration

Using a three-dimensional finite element model, this study characterizes groundwater flow in a costal plain of the Seto Inland Sea, Japan. The model characterization involved taking field data describing the aquifer system and translating this information into input variables that the model code uses to solve governing equations of flow. Geological geometry and the number of aquifers have been analyzed based on a large amount of geological, hydrogeological and topographical data. Results of study demonstrate a high correlation between the ground surface elevation and the groundwater level in the shallow coastal aquifer. For calibrating the numerical groundwater model, the groundwater flow was simulated in steady state. In addition, the groundwater level and trend in the transient state has also been elucidated. The numerical result provides excellent visual representations of groundwater flow, presenting resource managers and decision makers with a clear understanding of the nature of the types of groundwater flow pathways. Results build a base for further analysis under different future scenarios.     

Deuterium composition and flow path analysis as additional calibration targets to calibrate groundwater flow simulation in a coastal wetlands system

A conceptual model of the Lake Warden coastal wetlands system, Western Australia, was developed using hydraulic, chemical and stable isotopic data, and formed the basis for a groundwater flow model using the finite element numerical code (FEFLOW). The system to be modeled is complex. The surface water and groundwater within the wetlands system show varying salinity and isotopic composition over short distances and time frames. As a first step, the flow model was calibrated to observed groundwater levels measured since 2001 for both steady state and transient stresses. Particle tracking analysis was conducted to test the source areas of water discharging to the lakes within the wetlands system. The analysis was able to delineate the connectivity between the lakes in the wetland and the flow path. Enrichment of isotopic concentration is evident along a NE–SW transect and the data set provides a means for calibrating a detailed transport model. The study incorporates the varying deuterium composition of the water bodies directly into a transport model and a good match between observed and simulated temporal variations along the transect indicates that the model closely simulated the dynamics of water exchange between the lakes and groundwater within the system.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

西北内陆盆地冻结-冻融期的地下水补给与损耗

天山北麓平原区,每年冻结过程的时间远大于融化过程。在昌吉地下水均衡试验场,选择该平原的3种代表性土壤(粉质轻粘土、细砂、砂砾石),利用地中渗透仪观测冻结冻融期的地下水补给与损耗。根据不同岩性的地下水蒸发极限深度,设计了不同的地下水埋深。在冻结期,砂砾石、细砂岩性分布地区有利于浅埋型地下水资源的保护。在冻融期,无论何种地下水埋深水平,3种岩性的地下水都获得冻融水的补给,但细颗粒岩性的补给量相对较少且存在滞后效应,相反,粗颗粒岩性更有利于地下水资源的形成。冻结期一个月的地下水最大损耗量不超过253mm(砂砾石),而冻融期一个月的地下水最大补给量高达1133mm(细砂),冻融期是年内地下水的重要补给时期。     

Characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems in the North China coal basin

A conceptual framework for characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems is presented, based on the features of the site conditions and analyses of the water bursting and mine inundation events of the North China coal basin. Comprehensive analyses of the hydrogeological conditions for establishment of the three-dimensional groundwater flow systems of the North China coal basin revealed different vertical hydraulic connection paths or channels. These connections include karst collapse columns, fault or fracture zones, buried weathering zones, and fracture networks in aquitards, within the multilayered groundwater flow systems. Also examined, was the effect of the primary features of the flow system and those different connection paths on water bursting and discharge into the underground coal mines. It was demonstrated that appropriately identifying and adequately understanding the three-dimensional multilayered groundwater flow systems and those various vertical hydraulic connection mechanisms are critical for appropriately preventing water bursting hazards and predicting water discharge into underground mines. A valuable guideline is presented on characterization of water bursting and discharge into underground mines with multilayered groundwater flow systems.

---

Resumen Un cadre de travail conceptuel pour la caractérisation du jaillissement et de la vidange des eaux dans les mines à plusieurs niveaux d'écoulement est présenté, basé sur les conditions et l'analyse sur site du jaillissement des eaux et des inondations dans les bassins charbonneux du Nord de la Chine. L'analyse compréhensive des conditions hydrogéologiques des systèmes d'écoulement tridimensionnel de ces bassins, révèle des trajets ou des chenaux qui sont autant de connections hydrauliques verticales. Ces connections comprennent des effondrement karstiques en colonne, des zones de fractures et de failles, des zones d'altération, des réseaux de fractures dans les aquitards, dans le système hydrogéologique multicouche. L'influence des caractéristiques primaires du système d'écoulement et de ces différentes connections sur le jaillissement et la vidange dans les mines, est également examiné. Il a été démontré que l'identification appropriée et la compréhension adéquate du système d'écoulement souterrain multicouche et tridimensionnel, et des mécanismes variés de connection hydraulique verticale est critiquée pour prévenir de manière appropriée les risques de jaillissement et de vidange dans les mines. Des guidelines valorisables sont présentées, pour la caractérisation du jaillissement et de la vidange dans de tels systèmes.

---

Résumé Se presenta un marco conceptual para la de caracterización de la despresurización y descarga de agua a minas subterráneas con sistemas de acuíferos de múltiples niveles, el cual se basa en las condiciones del sitio, análisis de la despresurización de agua e inundaciones de minas en la cuenca de carbón del norte de la China.Un análisis comprehensivo de las condiciones hidrogeológicas para el establecimiento de los sistemas de flujo de agua de tri-dimensionales en la cuenca de carbón del norte de la China ha revelado diferentes canales hidráulicos verticales. Estas conexiones incluyen columnas de karst colapsadas, fallas o zonas de fracturas, zonas de regalita hundida y redes de fracturas en zonas de poca permeabilidad dentro de los sistemas de flujo de múltiples niveles.También se investigó el efecto de los rasgos primarios del sistema de flujo y los diferentes canales conectivos sobre la despresurización y descarga de agua a minas subterráneas. Se demostró que el entendimiento e identificación adecuada de los sistemas de acuíferos de múltiples niveles en tres dimensiones es crítico para prevenir apropiadamente la despresurización y descarga de agua a minas subterráneas. Se presentan lineamientos valiosos en relación con la despresurización y descarga de agua a minas subterráneas en sistemas de flujo de múltiple nivel.

     

基于Matlab小型潜水盆地地下水流动系统模拟研究

托特等人基于小型潜水盆地研究逐渐发展起来的地下水流动系统理论,是水文地质学发展中的一个重大突破.为了方便教学以及进一步完善小型潜水盆地模型中相关问题,尝试用Matlab的PDE工具箱解决这些问题.结果表明,运用Matlab的PDE Toolbox求解小型潜水盆地模型便捷高效;结果还表明不同流动系统水流相向汇流处及盆地下部两个角落均会形成水动力圈闭带.     

河西内陆河流域冰川融水近期变化

基于第一次中国冰川编目数据库,以国家气象台站的月降水与月气温资料为驱动数据,利用度日模型重建了河西内陆河流域各水系1961～2006年冰川物质平衡和融水径流序列,并与短期观测资料的对比验证表明,模型具有较高的可信度.河西内陆河流域1961～2006年平均冰川物质平衡为-49.5 mm/a,从东段石羊河流域到西段疏勒河、...