Regional groundwater flow modelling of Upper Chaj Doab of Indus Basin, Pakistan using finite element model (Feflow) and geoinformatics

A 3-D finite element model (Feflow) has been used for regional groundwater flow modelling of Upper Chaj Doab in Indus Basin, Pakistan. The thematic layers of soils, landuse, hydrology, infrastructure and climate were developed using Geographic Information System (GIS). The numerical groundwater flow model is developed to configure the groundwater equipotential surface, hydraulic head gradient and estimation of the groundwater budget of the aquifer. Integration of GIS with groundwater modelling and satellite remote sensing capabilities has provided an efficient way of analysing and monitoring groundwater status and its associated land conditions. The Arcview GIS software is used as additive tool to develop supportive data for numerical groundwater modelling, integration and presentation of image processing and modelling results.
The groundwater behaviour of the regional model shows a gradual decline in watertable from year 1999 onward. The persistent dry condition and high withdrawal rates play an influential role in lowering down the groundwater levels. Different scenarios were developed to study the impact of extreme climatic conditions (drought/flood) and variable groundwater abstraction on the regional groundwater system. The results of the study provide useful information regarding the behaviour of aquifer in order to organize management schemes on local and regional basis to monitor future groundwater development in the area.     

Groundwater vulnerability and risk mapping for the Basaltic aquifer of the Azraq basin of Jordan using GIS, Remote sensing and DRASTIC

Water consumption in Jordan already exceeds renewable freshwater resources by more than 20% and, after the year 2005, freshwater resources are likely to be fully utilised. Over 50% of supply derives from groundwater and this paper focuses on a small part of the northern Badia region of Jordan that is underlain by the Azraq groundwater basin where it has been estimated that annual abstraction stands at over 100% of the projected safe yield. While water supply is a crucial issue, there is also evidence to suggest that the quality of groundwater supplies is also under threat as a result of salinisation and an increase in the use of agrochemicals. Focusing on this area, this paper attempts to produce groundwater vulnerability and risk maps. These maps are designed to show areas of greatest potential for groundwater contamination on the basis of hydro-geological conditions and human impacts. All of the major geological and hydro-geological factors that affect and control groundwater movement into, through, and out of the study area were incorporated into the DRASTIC model. Parameters included; depth to groundwater, recharge, aquifer media, soil media, topography, and impact of the vadose zone. The hydraulic conductivity of the aquifer was not included in calculating the final DRASTIC index for potential contamination due to a lack of sufficient quantitative data. A Geographical Information System (GIS) was used to create a groundwater vulnerability map by overlaying the available hydro-geological data. The resulting vulnerability map was then integrated with a land use map as an additional parameter in the DRASTIC model to assess the potential risk of groundwater to pollution in the study area. The final DRASTIC model was tested using hydrochemical data from the aquifer. Around 84% of the study area was classified as being at moderate risk while the re mainder was classified as low risk. While the analysis of groundwater chemistry was not conclusive, it was encouraging to find that no well with high nitrate levels was found in the areas classified as being of low risk suggesting that the DRASTIC model for this area provided a conservative estimate of low risk areas. It is recognised that the approach adopted to produce the DRASTIC index was limited by the availability of data. However, in areas with limited secondary data, this index provides important objective information that could be used to inform local decision making.     

Comparative study of specific groundwater vulnerability of a karst aquifer in central Florida

The Floridan aquifer system (FAS) is known to be one of the most productive aquifer systems in the USA. With the FAS being a karst aquifer, it presents unique challenges to land use planners because of inherent vulnerabilities to contamination through direct connections between the aquifer and the surface. In this study a new Geographic Information Systems (GIS) -based index, the Karst Aquifer Vulnerability Index (KAVI), incorporates geologic layers used in intrinsic groundwater vulnerability models (GVMs) plus an epikarst layer specific to karst, with land use coverages to create a specific groundwater vulnerability model. The KAVI model was compared to another specific vulnerability model, the Susceptibility Index (SI). Tabulation of the percentage areas of vulnerability classes reveals major differences between the two models with SI suggesting greater vulnerability for the study area than KAVI. Validation of these two models found that KAVI vulnerability levels best reproduced spatially varying concentrations of nitrate in the aquifer. Sensitivity analysis, the application of a variation index and measuring the effective weights for each parameter included in KAVI confirmed the importance of closed depressions but also aquifer hydraulic conductivity. The inclusion of land use was justified; however, effective weight analysis determined its assigned weight was too high as used in the initial calculation of KAVI.     

毛乌素沙地海流兔河流域水文地质特征

海流兔河流域是毛乌素沙地的代表性流域,开展以流域为单元的研究,将为毛乌素沙地的流域水资源配置提供科学依据.通过系统分析海流兔河流域的含水层展布和地下水分布特征,创新性地将含水层厚度和地下水埋深的分区进行小波图像融合后,得到了新的水文地质参数分区.结果 表明:海流兔河流域的第四系及白垩系为该区域的巨厚含水层系统,其中白垩...     

Development and Implementation of a Bayesian-based Aquifer Vulnerability Assessment in Florida

The Florida Aquifer Vulnerability Assessment (FAVA) was designed to provide a tool for environmental, regulatory, resource management, and planning professionals to facilitate protection of groundwater resources from surface sources of contamination. The FAVA project implements weights-of-evidence (WofE), a data-driven, Bayesian-probabilistic model to generate a series of maps reflecting relative aquifer vulnerability of Florida’s principal aquifer systems. The vulnerability assessment process, from project design to map implementation is described herein in reference to the Floridan aquifer system (FAS). The WofE model calculates weighted relationships between hydrogeologic data layers that influence aquifer vulnerability and ambient groundwater parameters in wells that reflect relative degrees of vulnerability. Statewide model input data layers (evidential themes) include soil hydraulic conductivity, density of karst features, thickness of aquifer confinement, and hydraulic head difference between the FAS and the watertable. Wells with median dissolved nitrogen concentrations exceeding statistically established thresholds serve as training points in the WofE model. The resulting vulnerability map (response theme) reflects classified posterior probabilities based on spatial relationships between the evidential themes and training points. The response theme is subjected to extensive sensitivity and validation testing. Among the model validation techniques is calculation of a response theme based on a different water-quality indicator of relative recharge or vulnerability: dissolved oxygen. Successful implementation of the FAVA maps was facilitated by the overall project design, which included a needs assessment and iterative technical advisory committee input and review. Ongoing programs to protect Florida’s springsheds have led to development of larger-scale WofE-based vulnerability assessments. Additional applications of the maps include land-use planning amendments and prioritization of land purchases to protect groundwater resources.     

Evaluation of Recharge and Groundwater Dynamics of a Shallow Alluvial Aquifer in Central Ganga Basin,Kanpur (India)

Evaluation of recharge and groundwater dynamics of an aquifer is an important step for finding a proper groundwater management scenario. This has been performed on the basis of statistical Kendall Tau test to find a relationship between groundwater levels and hydro-meteorological parameters (e.g., precipitation, temperature, evaporation). Recharge to the aquifer was estimated for identification of critical areas/locations based on the analytical Soil and Water Assessment Tool. Moreover, spatiotemporal variability of groundwater levels has been quantified using space–time variogram. The overall characterization method has been applied to the shallow alluvial aquifer of Kanpur city in India. The analysis was performed using groundwater level data from 56 monitoring piezometer locations in Kanpur from March 2006 to June 2011. Groundwater level shows relatively higher correlation with temperature. Performance of the geostatistical model was evaluated by comparing with the observed values of groundwater level from January 2011 to June 2011 for two scenarios: “with limited spatiotemporal data” and “without spatiotemporal data.” It is evident that spatiotemporal prediction of groundwater level can be performed even for the unmonitored/missing data. This analysis demonstrates the potential applicability of the method for a general aquifer system.     

Characterization of Mechanisms and Processes Controlling Groundwater Recharge and its Quality in Drought-Prone Region of Central India (Buldhana,Maharashtra) Using Isotope Hydrochemical and End-Member Mixing Modeling

A thorough study on understanding of groundwater recharge sources and mechanisms was attempted by integrating the hydrogeological, geochemical and isotopic information along with groundwater dating and end-member mixing analysis (EMMA). This study was necessitated due to prolonged dryness and unavailability of freshwater in semi arid Deccan trap regions of Central India. In addition, groundwater resources are not characterized well in terms of their geochemical nature and recharge sources. The hydrogeochemical inferences suggest that aquifer I consists of recently recharged water dominated by Ca–Mg–HCO₃ facies, while groundwater in aquifer II shows water–rock interaction and ion exchange processes. Presence of agricultural contaminant, nitrate, in both aquifers infers limited hydraulic interconnection, which is supported by unconfined to semi-confined nature of aquifers. Groundwater in both aquifers is unsaturated with respect to carbonate and sulfate minerals indicating lesser water–rock interaction and shorter residence time. This inference is corroborated by tritium age of groundwater (aquifer I: 0.7–2 years old and aquifer II: 2–4.2 years old). Stable water isotopes (δ²H, δ¹⁸O) suggest that groundwater is a mixture of rainwater and evaporated water (surface water and irrigation return flow). EMMA analysis indicates three groundwater recharge sources with irrigation return flow being the dominant source compared to others (rainwater and surface waters). A conceptual model depicting groundwater chemistry, recharge and dynamics is prepared based on the inferences.

     

Geospatial Information and Environmental Isotopes for Hydrogeological Evaluation: Ras Alam El Rum,Northwestern Coast of Egypt

The present study demonstrates the potential of the geospatial information as well as isotopes and chemical analyses to elucidate the hydrogeological setting of the study area. The main objective of this research is to identify factors affecting the occurrence and quality of groundwater in the study area. The Ras Alam El Rum is the pilot area chosen in the northwestern coast of Egypt to evaluate the carbonate aquifer with emphasis on the Middle Miocene aquifer, which has been explored recently. The average of annual precipitation and surface runoff amounts were estimated to indicate possibility of groundwater recharge, where the water-bearing rocks are mainly fractured. The isotope data indicate that the groundwater was derived from rainwater, the recharge source is local precipitation, and no indication of seawater intrusion was recorded. The salinity of the groundwater is attributed to interaction between the water and aquifer rocks represented by limestone, dolomitic limestone, and marl. The mineral indices obtained from the geochemical model reveal that the groundwater is over-saturated with aragonite, calcite, dolomite, strontianite, and barite. The hydrogeological evaluation of the study area led to recommendations to determine of the suitable areas for exploration of the Middle Miocene aquifer and to protect groundwater from deterioration.     

From data gathering to resources assessment: A holistic view of petroleum geology

To integrate geological data to obtain an interpretation of the geology and natural resources of an area, we need a methodology that provides a holistic approach, addressing all basic issues in geological relations as well as uncertainties that arise with the evolution of basic geological knowledge of an area.In spite of the major role played by geoinformatics—the application of mathematics, statistics, and computer science to solve geological problems—we do not yet have a properly designed method for organizing geological data, including raw data, conceptual models, modeling results, and geological integration. Such a rapid method should provide for updating of existing interpretations based on new data or new theories.A systematic view of data integration and interpretation is important in petroleum exploration and petroleum engineering. We point out the more striking tools already available but of restricted use and some of the possible solutions for known problems that still require research and development.     

Groundwater Assessment at Manukan Island,Sabah: Multidisplinary Approaches

A rapid increase in the number of tourists has placed a heavy demand for freshwater on Manukan Island, a small island located offcoast of Kota Kinabalu, Sabah. Hydrochemical and numerical approaches have been applied in this study to evaluate the groundwater quality of unconfined aquifer lying in Manukan Island, Sabah, East Malaysia. This is vital to enhance better understanding about groundwater management. Hydrochemical analysis output indicated NaCl water type in sampling locations. Seawater intrusion is marked by its relatively high Na⁺, Mg²⁺, Cl⁻ and SO₄ ²⁻ concentrations. Hydrochemical analysis output clearly showed the influence of seawater in groundwater of Manukan Island. The numerical model output proved the influence of seawater in groundwater of Manukan Island by indicating the upconing process at the beneath of the pumping well. Current status of seawater intrusion in Manukan Island is about 14.6% of freshwater and seawater mixing ratio in low lying area of Manukan Island as simulated by SEAWAT-2000 model output. Numerical model SEAWAT-2000 output showed clearly that the upconing process is the possible route of seawater to influence the fresh groundwater aquifer chemistry in Manukan Island. The results have enhanced the current understanding of seawater intrusion in the study area. Future studies will focus on using numerical models to simulate and suggest suitable groundwater management plans in Manukan Island.     

Application of Magnetotelluric (MT) Resistivity to Imaging of Regional Three-Dimensional Geologic Structures and Groundwater Systems

Widespread definition of a groundwater system in three dimensions is necessary for the management and maintenance of groundwater resources. A magnetotelluric (MT) survey can be an effective geophysical prospecting method for imaging regional geological structures by measuring both shallow and deep resistivity. To demonstrate the capability of an MT survey to characterize a groundwater system, the Kumamoto area of central Kyushu in southwestern Japan was selected as a case study site because of its rich groundwater resources. Three-dimensional (3D) MT resistivity structure to a depth of 5?km was modeled by 1D inversion analysis of raw MT data and 3D interpolation of the resultant resistivity column data by the optimization principle method. Consequently, both deep and shallow aquifers were detected. A high-resistivity zone appears at depths between 500 and 2,000?m between the Futagawa?CHinagu faults and the Usuki?CYatsushiro tectonic line, which supports the existence of an aquiclude under the aquifer. The most important characteristic inferred from the 3D resistivity model is that the deep groundwater system below a depth of 1,000?m has two main flow paths. One path is likely to be through porous rocks because the low resistivity zone is regarded as tuff with sand and gravel, and the other flow path is interpreted to be through fractured zones along the Hinagu faults. Hence, the path and direction of the groundwater flows are probably controlled by geologic structures and the configuration of the active faults. These findings support the effectiveness of the MT method for investigating groundwater systems.     

Hydrogeologic Behavior of an Alluvial Aquifer,Salta Province,Argentina: Simulations of Hydraulic Conductivity Field,Groundwater Flow,and Chloride Migration

This work was designed to analyze the hydrogeological behavior of an alluvial aquifer in the River Mojotoro basin site in the Province of Salta, Argentina. The study area presents coarse-grained sediments with high infiltration capacity. The hydraulic conductivity field is affected by the physical heterogeneity of the medium and a geostatistical method, kriging, was used to construct this field from local measurements. The groundwater flow velocity field is computed using a hybrid mixed finite-element method and small variations in the flow directions were determined. Given the existence of a point source of chloride in the region under study, the spatial distribution of the generated plume in the saturated zone is analyzed. The advection-diffusion equation describing the transport of chloride in vertical and horizontal cross-sectional planes of the aquifer is discretized employing a method of characteristics combined with a mixed finite-element procedure, which provides an efficient method to eliminate spurious numerical oscillations and handle the convective term in the equation. The migration of chloride shows patterns, which are affected by the litologic heterogeneities of the aquifer and the regional groundwater flow.     

双重介质模型中地下水运动的环境示踪

为评价博茨瓦纳帕拉匹村非均质裂隙含水层的地下水资源，采用双重介质法及增加阻水水平屏障来模拟含水层，用排水和蒸发双重边界来模拟季节性河流，用溶质运移模型来再现示踪剂氟氯烃化合物（CFC）在地下水中的分布，用环境同位素氢和氧分析地下水的来源，用氯离子质量平衡法估计人渗量。对双重介质的不同有效空隙度的确定则采用了地下水龄示踪剂CFC的实测值分析得出。这些为地下水模型的真实、准确提供了几种独立、有效的分析手段。     

地表水和地下水耦合模型研究进展

地表水和地下水的转换及其过程耦合模拟,是水资源开发利用和科学评价的基础。本文归纳了国内外16 种地表水和地下水耦合模型,并根据耦合模型的特点和耦合方式,划分了三类耦合类型,即：紧密耦合型、半松散耦 合型和松散耦合型。对其中比较典型的“四水”转化模型、SWATMOD、MIKE- SHE、MODBRANCH 模型进行了对比分 析,剖析了不同耦合模型的优缺点,为地表水和地下水耦合模拟研究提供了很好参考。     

Application of Artificial Neural Networks to Complex Groundwater Management Problems

As water quantity and quality problems become increasingly severe, accurate prediction and effective management of scarcer water resources will become critical. In this paper, the successful application of artificial neural network (ANN) technology is described for three types of groundwater prediction and management problems. In the first example, an ANN was trained with simulation data from a physically based numerical model to predict head (groundwater elevation) at locations of interest under variable pumping and climate conditions. The ANN achieved a high degree of predictive accuracy, and its derived state-transition equations were embedded into a multiobjective optimization formulation and solved to generate a trade-off curve depicting water supply in relation to contamination risk. In the second and third examples, ANNs were developed with real-world hydrologic and climate data for different hydrogeologic environments. For the second problem, an ANN was developed using data collected for a 5-year, 8-month period to predict heads in a multilayered surficial and limestone aquifer system under variable pumping, state, and climate conditions. Using weekly stress periods, the ANN substantially outperformed a well-calibrated numerical flow model for the 71-day validation period, and provided insights into the effects of climate and pumping on water levels. For the third problem, an ANN was developed with data collected automatically over a 6-week period to predict hourly heads in 11 high-capacity public supply wells tapping a semiconfined bedrock aquifer and subject to large well-interference effects. Using hourly stress periods, the ANN accurately predicted heads for 24-hour periods in all public supply wells. These test cases demonstrate that the ANN technology can solve a variety of complex groundwater management problems and overcome many of the problems and limitations associated with traditional physically based flow models.     

Conceptualizing Anthropogenic Change in Fluvial Systems: Drainage Development on the Upper Embarras River,Illinois

While the geomorphic structure of many landscapes is materially molded, shaped and transformed by human practices, it is at the same time experienced and perceived through the mediation of cultural discourse and representation. The objective of this article is twofold: (1) to establish a conceptual framework for integrated analysis of human and geomorphic dimensions of landscape change; and (2) to illustrate the use of this framework in the assessment of a physical system severely affected by human agency. This perspective offers one way in which human agency can be synthesized and integrated into aspects of physical geography— specifically fluvial geomorphology— from both a conceptual and a pragmatic point of view. The upper Embarras River in east central Illinois, used as a case study, has been and continues to be physically affected by the cultural practice of agricultural drainage over the past century and a half.     

Agent-based modeling within a cyberinfrastructure environment: a service-oriented computing approach

Agent-based models (ABM) allow for the bottom-up simulation of dynamics in complex adaptive spatial systems through the explicit representation of pattern–process interactions. This bottom-up simulation, however, has been identified as both data- and computing-intensive. While cyberinfrastrucutre provides such support for intensive computation, the appropriate management and use of cyberinfrastructure (CI)-enabled computing resources for ABM raise a challenging and intriguing issue. To gain insight into this issue, in this article we present a service-oriented simulation framework that supports spatially explicit agent-based modeling within a CI environment. This framework is designed at three levels: intermodel, intrasimulation, and individual. Functionalities at these levels are encapsulated into services, each of which is an assembly of new or existing services. Services at the intermodel and intrasimulation levels are suitable for generic ABM; individual-level services are designed specifically for modeling intelligent agents. The service-oriented simulation framework enables the integration of domain-specific functionalities for ABM and allows access to high-performance and distributed computing resources to perform simulation tasks that are often computationally intensive. We used a case study to investigate the utility of the framework in enabling agent-based modeling within a CI environment. We conducted experiments using supercomputing resources on the TeraGrid – a key element of the US CI. It is indicated that the service-oriented framework facilitates the leverage of CI-enabled resources for computationally intensive agent-based modeling.     

Modeling Short-Term Spatial and Temporal Variability of Groundwater Level Using Geostatistics and GIS

Continuous depletion of groundwater levels from deliberate and uncontrolled exploitation of groundwater resources lead to the severe problems in arid and semi-arid hard-rock regions of the world. Geostatistics and geographic information system (GIS) have been proved as successful tools for efficient planning and management of the groundwater resources. The present study demonstrated applicability of geostatistics and GIS to understand spatial and temporal behavior of groundwater levels in a semi-arid hard-rock aquifer of Western India. Monthly groundwater levels of 50 sites in the study area for 36-month period (May 2006 to June 2009; excluding 3 months) were analyzed to find spatial autocorrelation and variances in the groundwater levels. Experimental variogram of the observed groundwater levels was computed at 750-m lag distance interval and the four most-widely used geostatistical models were fitted to the experimental variogram. The best-fit geostatistical model was selected by using two goodness-of-fit criteria, i.e., root mean square error (RMSE) and correlation coefficient (r). Then spatial maps of the groundwater levels were prepared through kriging technique by means of the best-fit geostatistical model. Results of two spatial statistics (Geary’s C and Moran’s I) indicated a strong positive autocorrelation in the groundwater levels within 3-km lag distance. It is emphasized that the spatial statistics are promising tools for geostatistical modeling, which help choose appropriate values of model parameters. Nugget-sill ratio (<0.25) revealed that the groundwater levels have strong spatial dependence in the area. The statistical indicators (RMSE and r) suggested that any of the three geostatistical models, i.e., spherical, circular, and exponential, can be selected as the best-fit model for reliable and accurate spatial interpolation. However, exponential model is used as the best-fit model in the present study. Selection of the exponential model as the best-fit was further supported by very high values of coefficient of determination (r ² ranging from 0.927 to 0.994). Spatial distribution maps of groundwater levels indicated that the groundwater levels are strongly affected by surface topography and the presence of surface water bodies in the study area. Temporal pattern of the groundwater levels is mainly controlled by the rainy-season recharge and amount of groundwater extraction. Furthermore, it was found that the kriging technique is helpful in identifying critical locations over the study area where water saving and groundwater augmentation techniques need to be implemented to protect depleting groundwater resources.     

黑河下游地下水波动带地下水时空分布模拟研究-FEFLLOW模型应用

从黑河下游地下水波动带水文循环出发,确定影响地下水变化的重要补给排泄项,根据水文地质结构对研究区水力传导系数和给水度进行确定,并获得研究区参数分区.通过运行FEFL-LOW模型,得到任意时段研究区地下水水头的空间分布.利用观测数据对模拟结果验证并进行参数调整.结果表明:(1)15个观测点平均误差约为0.59 m,最小平均误差为0.18 m,最大平均误差1.09 m;(2)地下水的空间变化为以东西河为中心向两侧加深,研究区南部岩石山地和北部东西居延海之间地下水位较深,其它区域地下水位在1～4 m之间;(3)从模拟初始(1990年),研究区每年水量都处于一种负均衡状态,2000年之后略有改观,2003年出现了正均衡.但不同区域地下水年际变化的趋势不同,研究区内上下游的水位有上升趋势,而中游绿洲区水位以下降趋势为主.     

区域地下水系统过程辨识

本文从水文地质的基本原理出发,从几个不同的侧面导出了地下水系统状态转移方程的一般形式。以山西省太原市3个地下水源地为例,用长系列观测资料对地下水系统进行了过程辨识。