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.     

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

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

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

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.     

湖北徐家山锑矿床流体包裹体特征及其意义

文章利用激光拉曼光谱和显微测温学方法,对湖北徐家山锑矿床成矿期的石英、重晶石和方解石中的流体包裹体进行了研究。研究表明,这些矿物中的流体包裹体主要有纯液体包裹体和液体包裹体(气相 液相)2类,其液体包裹体的气相成分为H2O±CO2±N2;石英、重晶石和方解石的均一温度分别为134~258℃、154~259℃和145~230℃,主要集中于150~200℃;流体的盐度w(NaCleq)和密度分别集中于3%~6%和0.90~0.96g/cm3。流体包裹体资料揭示出该矿床为典型的中低温热液锑矿床,其成矿流体为中低温、低盐度、中等密度热液。结合H、O、Sr、Pb同位素等研究结果,进一步推断该成矿热液主要是经深部循环演化的大气降水。     

新疆阿尔泰蒙库铁矿床的成矿流体及成矿作用

蒙库大型铁矿床赋存于上志留统—下泥盆统康布铁堡组变质火山-沉积岩系中,容矿岩石为石榴子石矽卡岩、变粒岩、浅粒岩和大理岩。矿体总体顺层分布,空间上与矽卡岩密切相关。研究表明,矽卡岩期石榴子石以发育玻璃质熔融包裹体、流体熔融包裹体和流体包裹体为特征,晚期矽卡岩阶段矿物中发育液相包裹体,变质期矿物中主要发育液相包裹体和含子矿物包裹体。矽卡岩期熔融包裹体的均一温度为1100℃,早期矽卡岩阶段流体包裹体均一温度变化于193～499℃,在450℃、350℃和230℃出现峰值。中期矽卡岩阶段均一温度变化于236～550℃,峰值为350℃。区域变质期均一温度介于132～513℃,在350℃、230℃和190℃出现峰值。流体包裹体的盐度w(NaCleq)介于1.23%～60.31%,流体密度变化于0.60～1.16g/cm3。石榴子石、石英和方解石的δ18OSMOW变化于0.2‰～8.4‰,δ18OH2O介于-5.1‰～5.33‰,δD为-127‰～-81‰,表明矽卡岩期成矿流体主要是岩浆水,混合少量大气降水;变质期流体主要为大气降水,为混合变质水。方解石δ13CPDB变化于-6.1‰～-2.3‰,表明流体中碳来自深部或地幔。成矿时代为早泥盆世早期(略晚于404～400Ma),成矿作用与矽卡岩的退化变质作用有关。     

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.     

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.