Study of shallow groundwater quality evolution under saline intrusion with environmental isotopes and geochemistry

Evolution of the shallow groundwater quality under saline intrusion in porous aquifer system has been studied with environmental isotopes and geochemistry in the Laizhou Bay area, China. Two campaigns of water sampling from various sources were carried out in spring and winter for environmental isotopic and chemical analyses. The origin of groundwater salinity from intrusion of both modern seawater and deep brine water was identified by analysing the correlations between ¹⁸O, D, T, Cl⁻, SO₄²⁻ and electrical conductivity. The results indicate that the brine is originated from evaporating and concentrating of intruded seawater and its δD and δ¹⁸O are different from modern seawater but similar to those of mixture of seawater with fresh groundwater. It is hard to distinguish the salinity origin in this area by the δD–δ¹⁸O relationship alone. The relations between δ¹⁸O and conductivity, Cl⁻ and SO₄²⁻ have been used to identify the salinity origin due to the distinct difference in salinity between the brine and seawater, conjunctively with use of T. A threshold of T = 12 TU was adopted to identify the origin of saline groundwater.     

Conceptual model of a coastal aquifer system in northern Greece and assessment of saline vulnerability due to seawater intrusion conditions

This paper refers to the development of a conceptual model for the management of a coastal aquifer in northern Greece. The research presents the interpretation and analysis of the quantitative (groundwater level recordings and design of piezometric maps) regime and the formation of the upcone within the area of investigation. Additionally it provides the elaboration of the results of chemical analyses of groundwater samples (physicochemical parameters, major chemical constituents and heavy metals and trace elements) of the area which were taken in three successive irrigation periods (July–August 2003, July–August 2004 and July 2005), in order to identify areas of aquifer vulnerability. The study identifies the areas where ion exchange phenomena occur, as well as the parts of the aquifer where the qualitative degradation of the aquifer system is enhanced. The paper, finally, assesses the lack of any scientific groundwater resources management of the area by the local water authorities, as well as the current practices of the existing pumping conditions scheme as applied by groundwater users.An erratum to this article can be found at     

Dynamic process analysis for the initiation and movement of the Donghekou landslide-debris flow triggered by the Wenchuan earthquake

The Donghekou landslide-debris flow was a remarkable geological disaster triggered by the Wenchuan earthquake in 2008. The dynamic process of a rapid landslide-debris flow is very complicated and can be divided into two aspects: the slope dynamic response of the earthquake and the mass movement and accumulation process. A numerical method combined with a finite difference method (FDM) and discrete element method (DEM) for simulation of landslide-debris flow under seismic loading is presented. The FDM and DEM are coupled through the critical sliding surface, initiation time and velocity. The dynamic response of the slope is simulated by the finite difference method, and critical sliding surface is determined using the earthquake response spectrum method. The landslide initiation time and the velocity are determined by time–history analysis. The mass movement and accumulation process is simulated using the discrete element method. Simulation results demonstrate that the maximum amplification coefficient of dynamic acceleration for the Donghekou slope is approximately 3.909, the initiation time of landslide is approximately 6.0 s, and the average initial velocity of the sliding mass is approximately 0.85 m/s. The failure of the slope is the result of elevation-orientated amplification effect and the sliding mass triggered with a small initial velocity. The numerical simulated result of the maximum sliding velocity is approximately 66.35 m/s, and the mass is disintegrated rapidly because of collision and free fall. The landslide velocity decreases when the flowing mass reaches a lower slope angle and gradually comes to a stop, and the total travel distance is approximately 2400 m.     

Numerical modelling of deformation and fluid flow in the Shuikoushan district, Hunan Province, South China

The Shuikoushan district, in Hunan Province, South China, contains major Pb–Zn–Au–Ag polymetallic mineralisation. Two groups of numerical models have been constructed to study the interactions between deformation and fluid flow in the district during the Yanshanian compression event (180 to 90 Ma). The first group includes district-scale conceptual models of coupled deformation and fluid flow during folding. The models show that fluid flow patterns are controlled by deformation within the fold system inferred for the district. During regional shortening and folding, fluids are generally focused towards the fold hinge/core areas along higher permeability layers (in particular Permian limestone units), in preference to flowing across the low permeability seal units (Permian and Jurassic terrestrial sequences). The efficiency of this fluid focusing can only be significantly increased if these folded seal units are allowed to undergo permeability increase as a result of tensile failure. The modelling results show that permeability enhancement localises mostly at fold hinges, dominantly within the silicified zone on the top of the Permian limestone unit. This process results in increased flow velocities and facilitates fluid focusing towards fold hinge/core locations at this silicified rock horizon. The second group includes deposit scale models for the Kangjiawan deposit, which is one of the two major deposits in the Shuikoushan district. The models show patterns of tensile failure, permeability creation, fluid focusing and mixing, and fracture development along a selected exploration cross section through the deposit. These results are consistent with the observed brecciation and mineralisation features. Regions of maximum brecciation in the district are associated with: (1) a combination of fold hinge and fault intersection locations (structural); and (2) the silicified zone and Permian limestone unit (lithological). Such brecciation zones are associated with extensive fluid focusing and mixing, and therefore represent the most favourable locations for mineralisation in the district. On the basis of this work, ideas for future research work and mineral exploration in the district are proposed.     

径流驱动土壤分离过程的影响因素及机制研究进展

土壤侵蚀包括土壤分离、泥沙输移和泥沙沉积3个子过程,研究这些过程发生、发展的水力、地形、土壤及地表特性等临界条件及各过程间相互影响、相互制约的机制,是建立土壤侵蚀过程模型继而准确预报土壤侵蚀的基础。为加强对径流驱动土壤分离过程的认识,推动土壤侵蚀过程与机理的研究,从径流水动力学特性、土壤特性、近地表特性、土壤分离过程的时空变异特性等方面对径流驱动的土壤分离过程影响因素及机制进行了总结和梳理,在此基础上,从土壤分离过程研究的时空尺度、土壤分离过程的动力学机制、细沟发育及其形态特征、网络结构对土壤分离过程的影响机制、土壤分离与泥沙输移的耦合机制、土壤分离测定方法与标准、土壤侵蚀阻力参数获取与预报等方面,展望了该领域有待深化的问题及今后的发展方向。     

袋装砂井爆夯处理软土地基的数值模拟方法及现场试验验证

袋装砂井爆夯法处理软土地基是利用炸药在设置有排水通道的软土中爆炸产生冲击和振动而使土体加固的方法。针对该法进行理论研究,提出了一种将袋装砂井爆夯处理软土地基的三维问题转化为二维平面应变问题的数值模拟方法：袋装砂井转化为等价砂墙;利用等效冲量原理,炮孔爆炸压力则转化为等效压力墙。数值模拟中考虑了土体骨架变形与孔隙水非达西渗流的耦合。对数值模拟的现场试验验证分析表明,沉降数值分析的结果与铁路宁启线软基处理现场测试结果具有很好的可比性。所提出的数值分析方法可模拟袋装砂井爆夯处理软土地基的超静孔隙水压产生和消散以及土体沉降变形的动态过程。     

玄武岩全岩元素含量快速、准确分析新技术:双铱带高温炉与LA-ICP-MS联用法

全岩样品主元素和微量元素的准确测试是地球科学研究的重要手段和依据.传统溶液进样ICP-MS元素测试方法制样和测试所需周期较长(两周以上).相对溶液进样法而言,激光剥蚀系统与ICP-MS联用测试全岩粉末熔融玻璃的方法可以大大地提高测试效率.在前人研究基础上,针对前人熔融制样中Pb、Zn元素的严重丢失问题,研究了Pb、Zn...     

High‐pressure metamorphic evolution of eclogite and associated metapelite from the Chuacús complex (Guatemala Suture Zone): Constraints from phase equilibria modelling coupled with Lu‐Hf and U‐Pb geochronology

As is common in suture zones, widespread high‐pressure rocks in the Caribbean region reached eclogite facies conditions close to ultrahigh‐pressure metamorphism. Besides eclogite lenses, abundant metapelitic rocks in the Chuacús complex (Guatemala Suture Zone) also preserve evidence for high‐pressure metamorphism. A comprehensive petrological and geochronological study was undertaken to constrain the tectonometamorphic evolution of eclogite and associated metapelite from this area in central Guatemala. The integration of field and petrological data allows the reconstruction of a previously unknown segment of the prograde P–T path and shows that these contrasting rock types share a common high‐pressure evolution. An early stage of high‐pressure/low‐temperature metamorphism at 18–20 kbar and 530–580°C is indicated by garnet core compositions as well as the nature and composition of mineral inclusions in garnet, including kyanite–jadeite–paragonite in an eclogite, and chloritoid–paragonite–rutile in a pelitic schist. Peak high‐pressure conditions are constrained at 23–25 kbar and 620–690°C by combining mineral assemblages, isopleth thermobarometry and Zr‐in‐rutile thermometry. A garnet/whole‐rock Lu‐Hf date of 101.8 ± 3.1 Ma in the kyanite‐bearing eclogite indicates the timing of final garnet growth at eclogite facies conditions, while a Lu‐Hf date of 95.5 ± 2.1 Ma in the pelitic schist reflects the average age of garnet growth spanning from an early eclogite facies evolution to a final amphibolite facies stage. Concordant U‐Pb LA‐ICP‐MS zircon data from the pelitic schist, in contrast, yield a mean age of 74.0 ± 0.5 Ma, which is equivalent to a U‐Pb monazite lower‐intercept age of 73.6 ± 2.0 Ma in the same sample, and comparable within errors with a less precise U‐Pb lower‐intercept age of 80 ± 13 Ma obtained in post‐eclogitic titanite from the kyanite‐bearing eclogite. These U‐Pb metamorphic ages are interpreted as dating an amphibolite facies overprint. Protolith U‐Pb zircon ages of 167.1 ± 4.2 Ma and 424.6 ± 5.0 Ma from two eclogite samples reveal that mafic precursors in the Chuacús complex originated in multiple tectonotemporal settings from the Silurian to Jurassic. The integration of petrological and geochronological data suggests that subduction of the continental margin of the North American plate (Chuacús complex) beneath the Greater Antilles arc occurred during an Albian‐Cenomanian pre‐collisional stage, and that a subsequent Campanian collisional stage is probably responsible of the amphibolite facies overprint and late syncollisional exhumation.