Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

定水头注水引起的含水层水平运动和应变

基于含水层固体颗粒与孔隙水不可压缩的假设 ,本文导出了单井注水情况下泰斯承压含水层水平运动速度与水头之间的基本关系式。然后利用注水井壁处的应力、应变边界条件 ,进一步导出了单井定水头注水引起的泰斯承压含水层水平运动速度、位移和应变解析表达式。该水平位移与应变由两部分组成 :一部分为由注水压力本身引起的经典弹性力学解项 ,它仅随半径而变化 ,与注水时间无关 ;另一部分为由地下水头变化引起的水动力学位移和应变解项。其中 ,含水层水动力学水平位移随时间加长呈指数增长特征 ,水动力学径向应变则表现为近井处拉张、远井处挤压的分区特征 ,且近井拉张区随时间加长逐渐向外扩展。单井注水含水层水动力学水平位移、应变解的导出 ,完善和发展了单孔内压经典弹性平面力学问题解     

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.     

A comparative study of four vulnerability mapping methods in a detritic aquifer under mediterranean climatic conditions

This paper presents the results of a comparative study relating to the application of four vulnerability mapping methods, GOD, AVI, DRASTIC and SINTACS, in a pilot detritic aquifer situated in NW Morocco, known as the Martil–Alila aquifer. The principal objective of this work is to determine the most suitable such methods for this aquifer type within a Mediterranean context, and to show the effect of the rainfall variations that are characteristic of the Mediterranean climate on the degree of vulnerability. The methods applied distinguish five classes of vulnerability, these being irregularly divided up in space, with the division varying according to the method in question. The vulnerability maps obtained by the different methods strongly suggest that the eastern half of the aquifer is more vulnerable to contamination than the western half, for all hydrological situations. The effect of climatic conditions on the degree of vulnerability is well represented by the DRASTIC, according to which the aquifer is moderately to strongly vulnerable during humid hydrological years and weakly to moderately vulnerable during dry ones. For the other methods, this climatic effect is limited to the area occupied by the two predominant classes (“High” and “Low” for GOD and “High” and “Moderate” for SINTACS) while it is null for AVI. In conclusion, DRASTIC appears the most suitable for mapping the vulnerability to contamination of Mediterranean coastal detritic aquifers such as the Martil–Alila aquifer.     

Occurrence and evolution of the Xiaotangshan hot spring in Beijing, China

Thermal groundwater occurs in bedrock aquifers consisting of the dolomite of the Wumishan Group of the Jixianin System and the Cambrian carbonate in the Xiaotangshan geothermal field near the northern margin of the North China Plain, China. The hot water in the geothermal field of basin-type discharges partly in the form of the Xiaotangshan hot spring under natural conditions. The hot water has TDS of less than 600 mg/L and is of Na·Ca-HCO₃ type. The geothermal water receives recharge from precipitation in the mountain area with elevation of about 500 m above sea level to the north of the spring. Thermal groundwater flows slowly south and southeast through a deep circulation with a residence time of 224 years estimated with the Ra–Rn method. The Xiaotangshan hot spring dried up in the middle of the 1980s owing to the increasing withdrawal of the hot water in the geothermal field in the past decades. The water level of the geothermal system still falls continually at an annual average rate of about 2 m, although water temperature changes very little, indicating that the recharge of such a geothermal system of basin-type is limited. Over-exploitation has a dramatic impact on the geothermal system, and reduction in exploitation and reinjection are required for the sustainable usage of the hot water.     

Deciphering potential groundwater zone in hard rock through the application of GIS

Remote Sensing and Geographic Information System has become one of the leading tools in the field of hydrogeological science, which helps in assessing, monitoring and conserving groundwater resources. It allows manipulation and analysis of individual layer of spatial data. It is used for analysing and modelling the interrelationship between the layers. This paper mainly deals with the integrated approach of Remote Sensing and geographical information system (GIS) to delineate groundwater potential zones in hard rock terrain. The remotely sensed data at the scale of 1:50,000 and topographical information from available maps, have been used for the preparation of ground water prospective map by integrating geology, geomorphology, slope, drainage-density and lineaments map of the study area. Further, the data on yield of aquifer, as observed from existing bore wells in the area, has been used to validate the groundwater potential map. The final result depicts the favourable prospective zones in the study area and can be helpful in better planning and management of groundwater resources especially in hard rock terrains.     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.     

Groundwater recharge of carbonate aquifers of the Silesian-Cracow Triassic (southern Poland) under human impact

A Triassic carbonate unit has been intensively drained by zinc and lead ore mines and numerous borehole fields since the nineteenth century. Its groundwater recharge has increased due to: pumping of water from boreholes, mining activity, and urbanization. An approach to determine the amounts of the recharge at a variety of spatial scales is presented in the paper. Different methods were used to identify and quantify recharge components on a regional and local scale: mathematical modelling was performed for four aquifers included in an aquifer system, an analytical estimation based on the assumption that an average recharge is equal to the average discharge of the hydrogeological system—for six man-made drainage centres, and the method of water level fluctuation (WLF) was applied in one observation borehole. Results of modelling have been supplemented by observation of environmental tracers (δ¹⁸O, δ²H, ³H), noble gases temperatures, and ⁴He_exc in groundwater. The regional aquifer’s current recharge according to estimations performed by means of modelling varies from 39 to 101 mm/year on average. Depending on the aquifer site the average precipitation ranges from 779 to 864 mm/year. In the confined part of the aquifer average recharge ranges from 26 to 61 mm/year. Within outcrops average recharge varies from 96 to 370 mm/year. Current recharge estimated by the analytical method for man-made drainage centres varies from 158 up to 440 mm/year. High values are caused by different recharge sources like precipitation, induced leakage from shallow aquifers, and water losses from streams, water mains and sewer systems. Pumping of water, mining and municipal activities constitute additional factors accounting for the intensified recharge.     

渤海湾盆地新生界生油岩系底界面温度分布

依据渤海湾盆地2000余口测温井的地温梯度数据、地层岩性描述、分层数据以及钻井资料,计算了该盆地各生油岩系底界面的温度。统计结果表明:渤海湾盆地沙河街组大部分凹陷区地层底界面温度介于90℃至150℃,目前仍具有大量生油的温度条件,而在隆起或一些凸起地区,该地层组段温度普遍小于90℃,未能达到生油的温度指标。东营组和孔店组地层也仍具有一定的生油温度条件。研究还表明:地层温度与地层界面埋深密切相关,温度随界面埋深的增大而升高,沉积厚度大的凹陷区地层界面温度大于沉积厚度小的凸起区或斜坡地带,说明地层界面埋深是决定地温高低的主要控制因素,而地温梯度对地层界面温度的影响相对较小。