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.     

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.     

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.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

河南洛宁段河金矿流体包裹体研究和矿床成因

河南省洛宁县段河石英脉型金矿主要包括石寨沟和岭东两个矿区,分别由3～4条含金石英脉构成。矿化过程从早到晚包括石英-黄铁矿、石英-多金属硫化物和石英-碳酸盐等3个阶段．其中中阶段金矿化最强,次为早阶段。各阶段石英中流体包裹体以气液两相包裹体为主．次为纯液体包裹体。激光拉曼测试表明,气液两相包裹体的液相为H2O,气相主要为Ho和CO2混合、纯H2O,次为纯CO2;纯液体包裹体为纯H2O。石寨沟矿区包裹体均一温度从早到晚依次为240．9～315．9℃．188．7～304．5℃,137．3～259．3℃：流体盐度变化依次为（6．74～12．85）wt%NaCl．eq,（2．41～8．68）wt％NaCl．eq,（2．24-7．86）wt%NaCl．eq。岭东矿区均一温度从早到晚依次为303．7-343．1℃,251．8-325．4℃,305．7～355．0℃：流体盐度变化依次为（5．11～11．70）wt％NaCl．eq,（2．74-10．11）wt％NaCl．eq,（0．53-6．74）wt％NaCl．eq。两矿区主成矿期流体均为中温、低盐度,早阶段流体为改造热液和变质热液的混合体,含一定量CO2,且流体CO2含量和盐度从早到晚逐渐降低。石寨沟矿区包裹体均一温度逐渐降低,而岭东矿区包裹体均一温度先降后升,加之岭东矿区各阶段成矿温度均高于石寨沟矿区．表明成矿流体系统主要受岩浆热驱动,岭东矿区更靠近岩体,且在晚阶段又有脉动性的岩浆加热．段河金矿区南部存在隐伏岩体。     

无外加流体、350℃和差异应力条件下硫化物再活化实验研究

红透山块状硫化物矿石主要成分为黄铁矿、磁黄铁矿、黄铜矿、闪锌矿和石英、角闪石、黑云母等脉石矿物。将此矿石烘干后作为试料置于岩石三轴应力试验机,在13h内将轴压、围压和温度分步升至1276MPa、414MPa和350℃,然后在空气中自然冷却至室温。实验产物中黄铁矿、石英、长石等以脆性碎裂为主,而磁黄铁矿、黄铜矿、闪锌矿和云母等以塑性变形为主。再活化的黄铜矿、磁黄铁矿和少量闪锌矿呈脉和网脉穿插黄铁矿碎斑。实验结果表明,即使没有外来流体加入．构造动力作用所导致的流体包裹体破坏所释出的流体,就足以使成矿物质发生活化转移和再活化。     

Validation of a finite difference method for the simulation of vortex-induced vibrations on a circular cylinder

The Karman Vortex Street generated by a circular cylinder is investigated by the numerical solution of the compressible Navier–Stokes equations in the incompressible Mach number range (Mach<0.3) using the Beam and Warming implicit scheme. The agreement with the fully incompressible projection method (Chorin, 1968) is fairly good while convergence time is very much better. The investigation suggests that the compressible Navier–Stokes equations may be used as an efficient alternative to study incompressible flows as well. Mach numbers just below 0.3 are enough to simulate incompressible flow behavior and at the same time do not cause numerical ill-conditioning in the solution. In addition, some relevant features of the vortices generated and carried by the wake of the cylinder could be fairly well captured.