Delineation of groundwater protection zones by the backward particle tracking method: theoretical background and GIS-based stochastic analysis

The backward particle tracking method, an effective and powerful tool that can be used to delineate groundwater protection zones, is presented. The theoretical background and insights on the applicability of this method are provided. Moreover, the present work enriches the backward particle tracking method with an uncertainty analysis concerning the porosity values, applying a Monte Carlo (MC) approach, coupled with the use of geographical information systems (GIS). As an application example, a wellfield in the Komotini area, Greece, is investigated. The present study may serve as a potential guideline for wellfield delineation, particularly in areas like Greece where lack of data related to the hydrogeological system is often a problem.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.     

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.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

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.     

内蒙阿右旗铁板井镍矿矿床特征及成因探讨

铁板井绦矿床赋存于辉橄岩体中,辉橄岩体具有全岩矿化的特点,岩体及镍矿体的形成经历了三种不同的成矿过程:首先是岩浆深源液态重力分异作用,而后发生岩浆深源熔离-贯入(成岩、成矿)作用,形成岩体和镍矿体,最后在岩浆期后热液作用的叠加下完成整个成矿过程.该矿床是典型的主要由岩浆熔离作用形成的岩浆矿床.     

引水工程秦巴段隧洞地应力模拟及工程地质问题

利用Ansys有限元软件对三峡引水工程秦巴段线路不同深度、不同截面形态隧洞围岩的应力重分布情况进行模拟计算,得到圆形隧洞、城门形隧洞和马蹄形隧洞围岩的应力数值和等值线图.利用库仑一纳维尔强度准则,对花岗岩区及灰岩区隧洞围岩的剪切破坏进行分析;利用三种岩爆应力判别指标,对岩爆进行初步预测;进而对不同深度、不同截面形态的隧洞围岩的稳定性进行分析;最后对引水工程隧洞设计中截面形态的选择给出了初步建议.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

陆面过程中稳定水同位素逐日变化的数学模拟——引入稳定同位素效应的CLM实例分析

将稳定同位素效应引入CLM(Community Land Model),并对巴西马瑙斯站在平衡年的稳定水同位素的逐日变化进行模拟和分析.结果表明: 降水、水汽和地表径流中δ18O存在明显的季节变化,并与相应的水量存在显著的负相关关系,但凝结物中δ18O与地面凝结量存在显著的正相关关系,蒸发水汽中δ18O与蒸发量之间无显著的相关关系.受土壤贮水削峰功能的影响,表层土壤和根区水中δ18O的季节变化全无.植被层蒸发水汽中稳定同位素的丰度与大气的干湿程度存在密切联系: 当降水量少时,大气干燥,植被层的蒸发较少,植被蒸发中δ18O较高;当降水量较大时,空气湿润,植被层的蒸发量较大,蒸发中δ18O则较低.植被蒸腾中δ18O的变化与源区水体中δ18O的变化保持一致,尤其是与根区水中的δ18O.由于地下径流直接源自根区水的补充,因此,地下径流中δ18O等于根区水中的δ18O.模拟结果还显示,降水MWL (大气水线)的梯度项和常数项均比全球平均MWL略偏小.尽管主要来自降水的贡献,但地表径流和植被层水体的MWLs与降水MWL存在较大的差异,这一方面与两类水体在蒸发过程中的稳定同位素的富集作用有关,另一方面与CLM模拟的水量有关.大气水汽线与降水的MWL的梯度值相近,说明大气水汽与降水近似处于稳定同位素平衡状态.另外,模拟的地面的凝结线与植被层的凝结线均与全球大气水线相近,且具有非常高相关程度,说明CLM的模拟是合理的.