Hydrochemical characteristics and salinization processes of groundwater in the shallow aquifer of Eastern Laizhou Bay,China

Groundwater salinization has become a crucial environmental problem worldwide and is considered the most widespread form of groundwater contamination in the coastal zone. In this study, a hydrochemical investigation was conducted in the eastern coastal shallow aquifer of Laizhou Bay to identify the hydrochemical characteristics and the salinity of groundwater using ionic ratios, deficit or excess of each ions, saturation indices and factor analysis. The results indicate that groundwater in the study area showed wide ranges and high standard deviations for most of hydrochemical parameters and can be classified into two hydrochemical facies, Ca²⁺‐Mg²⁺‐Cl^‐ facies and Na⁺‐Cl^‐ facies. The ionic ratio, deficit or excess of each ions and SI were applied to evaluate hydrochemical processes. The results obtained indicate that the salinization processes in the coastal zones were inverse cation exchange, dissolution of calcite and dolomite, and intensive agricultural practices. Factor analysis shows that three factors were determined (Factor 1: TDS, EC, Cl^‐, Mg²⁺, Na⁺, K⁺, Ca²⁺ and SO₄^2‐; Factor 2: HCO₃^‐ and pH; Factor 3: NO₃^‐ and pH), representing the signature of seawater intrusion in the coastal zone, weathering of water–soil/rock interaction, and nitrate contamination, respectively. Copyright © 2011 John Wiley & Sons, Ltd.     

A quadratic element method for evaluating groundwater flow by the inversion of surface tilt with application to the Tono Area, Japan

By modifying a previous method with constant elements, we developed a quadratic element method for more accurately estimating groundwater flow by the inversion of tilt data. In this method: (1) a region of groundwater flow is divided into quadratic elements in which the change in groundwater volume per unit volume of rock (Δv) and the Skempton coefficient (B) vary in a quadratic manner with the coordinates, (2) the values of Δv are set to zero at the boundaries of the region of groundwater flow and (3) the sum of the squared second derivatives of Δv is adopted as a constraining condition that is weighted and added to the sum of the squared errors in tilt. First, analyses were performed for a flow model to determine the accuracy of this method for estimating groundwater flow and also to clarify the effect of the assumed size of a region of groundwater flow. These analyses showed that the quadratic element method proposed in this study gives a much better estimation of Δv than the constant element method and that a large region of groundwater flow should be assumed, rather than a small region, since the values of Δv at points outside of the actual region of groundwater flow are estimated to be nearly zero when a large region is assumed while these values are greatly overestimated when an excessively small region is assumed. Finally, the quadratic element method was applied to the site of the Mizunami Underground Research Laboratory in the Tono area, Japan. Inverse analyses were performed for tilt data measured by four tiltmeters with a resolution of 10⁻⁹ radians during the excavation of two shafts under the assumption that the rock mass is an isotropic and homogeneous half- space. The results showed that the method proposed in this study reproduced the tilt data very accurately. Thus, the distribution of Δv was estimated without sacrificing the reproducibility of the tilt data. The contour maps of B(1 + ν)Δv (ν: Poisson’s ratio) showed that the heterogeneous flow of groundwater occurred at the site and that groundwater volume decreased mainly in the area surrounded by two faults. The latter result is consistent with the finding obtained by previous investigations that these faults have low permeability in the direction perpendicular to the strike and may act as a flow barrier.     

In-situ neutron flux,Cl production and groundwater evolution in crystalline rocks at Stripa, Sweden

The meteoric input of³⁶Cl due to cosmogenic or nuclear-weapon-produced³⁶Cl cannot contribute significantly to the³⁶Cl present in the saline groundwaters (up to 700 mg l⁻¹ Cl⁻) from the Stripa granite. The extent of in-situ production of³⁶Cl has been estimated on the basis of the neutron fluxes within the granite and its surrounding leptite. The³⁶Cl present in the groundwaters is attributed to either admixture of labelled Cl⁻ from the leptite with Cl⁻ from the granite or to the total derivation of groundwater chlorinity within the leptite followed by radiochemical ingrowth of³⁶Cl during subsequent groundwater residence within the granite. The chloride derived from the leptite may be either matrix chloride or chloride from an external source which has had a long residence time within the leptite. The implications of³⁶Cl in-situ production for the estimation of groundwater residence times and for the geochemical evolution of groundwater chlorinity are discussed.     

In-situ neutron flux,36Cl production and groundwater evolution in crystalline rocks at Stripa,Sweden

The meteoric input of³⁶Cl due to cosmogenic or nuclear-weapon-produced³⁶Cl cannot contribute significantly to the³⁶Cl present in the saline groundwaters (up to 700 mg l⁻¹ Cl⁻) from the Stripa granite. The extent of in-situ production of³⁶Cl has been estimated on the basis of the neutron fluxes within the granite and its surrounding leptite. The³⁶Cl present in the groundwaters is attributed to either admixture of labelled Cl⁻ from the leptite with Cl⁻ from the granite or to the total derivation of groundwater chlorinity within the leptite followed by radiochemical ingrowth of³⁶Cl during subsequent groundwater residence within the granite. The chloride derived from the leptite may be either matrix chloride or chloride from an external source which has had a long residence time within the leptite. The implications of³⁶Cl in-situ production for the estimation of groundwater residence times and for the geochemical evolution of groundwater chlorinity are discussed.     

Hydrogeochemical modeling of groundwater chemical environmental evolution in Hebei Plain

Using the hydrogeochemical modeling method, the groundwater chemical environmental problems of the Hebei Plain which involve increasing of hardness and total dissolved solids in piedmont area and mixing of saline water with fresh water in middle-eastern area are studied. The water-rock interactions and mass transfer along a ground-water flow path and in mixing processes are calculated. Thus the evolution mechanisms of the groundwater chemical environment are brought to light. Project supported by the National Natural Srience Foundation of China.     

Hydrochemical and isotopic characterization of groundwater in the southeastern part of the Plateaux Region,Togo

Hydrochemical, inverse geochemical modelling and isotopic approaches are used to assess the hydrogeochemical evolution of groundwater from the basement aquifers in the southeastern part of the Plateaux Region, Togo. Groundwater originates from present-day rainwater infiltration and is mostly fresh and slightly acidic to neutral. Hydrochemical facies are predominantly mixed cations-HCO₃ associated with Ca/Mg-Cl, Na-HCO₃ and Na-Cl water types in equilibrium with kaolinite and Ca/Mg-smectites. They are related to silicates hydrolysis, anthropogenic contamination, nitrification/denitrification, mixing along flowpaths and dissolution/precipitation of secondary minerals. The pattern of flow paths is in accordance with an increasing trend in total dissolved solids (TDS) toward the potentiometric depression located in the central and southern parts of the aquifer system. Inverse geochemical modelling using the NETPATH-WIN model showed the relative importance of biotite, plagioclase and amphibole weathering and dissolution of secondary carbonate minerals along the flowpath, suggesting that an abundance of minerals is not necessarily the main factor controlling the groundwater chemistry evolution.     

Thermal evolution of the Ryoke metamorphic belt, southwestern Japan: Tectonic and numerical modeling

Abstract The Ryoke metamorphic belt of southwestern Japan is composed of Cretaceous Ryoke granitoids and associated metamorphic rocks of low-pressure facies series. The Ryoke granitoids are divided into sheet-like bodies (e.g. Gamano granodiorite) and stock-like bodies. The Gamano granodiorite intruded concordantly into the high-grade metamorphic rocks without development of a contact metamorphic aureole, and the intrusion ages of the granodiorite are similar to the ages of thermal peak of the low pressure (low-P) metamorphism. It is suggested that the low-P Ryoke metamorphism resulted from the intrusion of the Gamano granodiorite. In this study, a simple 1-D numerical model of conductive heat transfer was used to evaluate the thermal effects of emplacement of the Gamano granodiorite. Calculated temperature-time ( T-t ) paths are characterized by a rapid increase of metamorphic temperature and a relatively short-lived period of high temperature. For example, the T-t path at the 15-km depth is characterized by a rapid average increase in temperature of 1.4 × 10-³°C/year and high temperatures for < ca 0.5 Ma. The calculated peak temperature for each depth is nearly equal to the petrologically estimated value for each correlated metamorphic zone. The results suggest that the magma-intrusion model is one possible thermal model for low-pressure facies series metamorphism.     

On the nature of cosmic dust in sedimentary rocks

In the paper, the chemical composition and origin of cosmic dust, as well as its distribution in the geological past, are discussed. Special attention is given to the problem of the difference in cosmic dust from volcanic dust on the basis of indicators such as cosmogenic helium and iridium.     

煤矿采掘过程中煤岩体电磁辐射特征及应用

运用研制的KBD5矿用本安型电磁辐射监测仪测试了煤矿采掘过程中工作面煤岩体的电磁辐射，分析了电磁信号与采掘工艺及煤岩动力灾害危险性等的关系。研究表明，煤矿采掘过程伴随产生电磁辐射，电磁辐射是煤岩体受到采动影响后应力重新分布或变形破裂趋向新平衡的结果；有煤与瓦斯突出和冲击地压危险时，有明显的电磁异常前兆；采取防治措施后，电磁辐射显著下降。电磁辐射技术在煤矿可以用于预测煤与瓦斯突出、冲击地压等煤岩动力灾害。     

Formation and deformation processes of late Paleogene sedimentary basins in southern central Hokkaido,Japan: Paleomagnetic and numerical modeling approach

Formation and deformation processes of the late Paleogene sedimentary basins related to a strike–slip fault system in southern central Hokkaido are described by a combination of paleomagnetic study and numerical analysis. After correction of the Miocene counter‐clockwise rotation associated with back‐arc opening of the Japan Sea, paleomagnetic declination data obtained from surface outcrops in the Umaoi and Yubari areas show significant easterly deflections. Although complicated differential rotation is anticipated as a result of recent thrust movements, clockwise rotation in the study areas is closely linked with development of the Paleogene Minami‐naganuma Basin as a pull‐apart depression along the north–south fault system. Numerical modeling suggests that 30 km of strike–slip is required to restore the distribution and volume of the Minami‐naganuma Basin. The relative slip rate on the long‐standing fault system is about 10 mm/yr, which corresponds to global‐scale plate motion. It has inevitably caused regional rearrangement of the eastern Eurasian margin. A rotation field simulated by simplified dextral motion using dislocation modeling basically accords with the paleomagnetic data around the pull‐apart basin.     

Hydrochemical characteristics of the Chita River basin

Human-induced factors responsible for the formation of the Chita River hydrochemical regime in the Amur River basin are discussed. An environmental map of the Chita River basin is compiled. The available results of hydrochemical observations over surface water and groundwater discharged from reclaimed lands of Chita region are generalized. The main factors affecting the hydrochemical regime of the Chita River are specified.     

Tectonic evolution of lower crustal rocks in an exposed magmatic arc section in the Hidaka metamorphic belt, Hokkaido, northern Japan

Abstract : The Hidaka metamorphic belt consists of an island-arc assembly of lower to upper crustal rocks formed during early to middle Paleogene time and exhumed during middle Paleogene to Miocene time. The tectonic evolution of the belt is divided into four stages, D₀rs, D₁, D₂rs, and D₃, based on their characteristic deformation, metamorphism, and igneous activity. The premetamorphic and igneous stage (D₀) involves tectonic thickening of an uppermost Cretaceous and earliest Tertiary accretionary complex, including oceanic materials in the lower part of the complex. D₁ is the stage of prograde metamorphism with increasing temperatures at a constant pressure during an early phase, and with a slight decrease of pressure at the peak metamorphic phase, accompanying flattening of metamorphic rocks and intrusions of mafic to intermediate igneous rocks. At the peak, incipient partial melting of pelitic and psammitic gneisses took place in the amphibolite–granulite facies transition zone, the melt and residuals cutting the foliations formed by flattening. In the deep crust, large amounts of S-type tonalite magma formed by crustal anatexis, intruded into the granulite facies gneiss zone and also into the upper levels of the metamorphic sequence during the subsequent stage. During D₁ stage, mafic and intermediate magmas supplied and transported heat to form the arc-type crust and at the same time, the magmatic underplating caused extensional doming of the crust, giving rise to flattening and vertical uplifting of the crustal rocks. D₂ stage is characterized by subhorizontal top-to-the-south displacement and thrusting of lower to upper crustal rocks, forming a basal detachment surface (décollement) and duplex structures associated with intrusions of S-type tonalite. Deformation structures and textures of high-temperature mylonites formed along the décollement, as well as the duplex structures, show that the D₂ stage movement occurred under a N-S trending compressional tectonic regime. The depth of intra-crustal décollement in the Hidaka belt was defined by the effect of multiplication of two factors, the fraction of partial melt which increases downward, and the fluid flux which decreases downward. The crustal décollement, however, might have extended to the crust-mantle boundary and/or to the lithosphere and asthenosphere boundary. The subhorizontal movement was transitional to a dextral-reverse-slip (dextral transpression) movement accompanied by low-temperature mylonitization with retrograde metamorphism, the stage defined as D₃. The crustal rocks from the basal décollement to the upper were tilted eastward on the N–S axis and exhumed during the D₃ stage. During D₂ and D₃ stages, the intrusion of crustal acidic magmas enhanced the crustal deformation and exhumation in the compressional and subsequent transpressional tectonic regime.     

On the theory of depositional remanent magnetization in sedimentary rocks

The effects of various factors such as thermal agitation, coagulation, anisotropy of susceptibility, and shape irregularity on the alignment of magnetic carrier grains during the process of acquisition of depositional remanent magnetization in sedimentary rocks is discussed.     

The groundwater recharge regime of some slightly metamorphosed neoproterozoic sedimentary rocks: an application of natural environmental tracers

Isotope data of precipitation and groundwater in parts of the Voltaian Basin in Northern Ghana were used to explain the groundwater recharge regime in the area. Groundwater recharge is an important parameter in the development of a decision support system for the management and efficient utilization of groundwater resources in the area. It is therefore important to establish the processes and sources of groundwater recharge. δ¹⁸O and δ²H data for local precipitation suggest enrichment relative to the Global Meteoric Water Line (GMWL) and indicate that precipitation takes place at a relative humidity less than 100%. The groundwater data plot on an evaporation line with a slope of 5, suggesting a high degree of evaporative enrichment of the precipitation in the process of vertical infiltration and percolation through the unsaturated zone into the saturated zone. This finding is consistent with the observation of high evapotranspiration rates in the area and ties in with the fact that significant clay fraction in the unsaturated zone limits vertical percolation and thus exposes the percolating rainwater to the effects of high temperatures and low humidities resulting in high evapotranspiration rates. Groundwater recharge estimates from the chloride mass balance, CMB, method suggest recharge in the range of 1.8–32% of the annual average precipitation in the form of rainfall. The highest rates are associated with areas where open wells encourage significant amount of groundwater recharge from precipitation in the area. In the northern parts of the study area, groundwater recharge is lower than 12%. The recharge so computed through the application of the CMB methodology takes on a spatial distribution akin to the converse of the spatial pattern of both δ¹⁸O and δ²H in the area. As such, the locations of the highest recharge are associated with the most depleted values of the two isotopes. This observation is consistent with the assertion that low vertical hydraulic conductivities slow down vertical percolation of precipitation down to the groundwater water. The percolating precipitation water thus gets enriched in the heavier isotopes through high evapotranspiration rates. At the same time, the amount of water that finally reaches the water table is considerably reduced. Copyright © 2013 John Wiley & Sons, Ltd.     

Multivariate statistical analysis of geochemical data as indicative of the hydrogeochemical evolution of groundwater in a sedimentary rock aquifer system

The study of groundwater hydrogeochemistry of the Paleozoic Basses-Laurentides sedimentary rock aquifer system in Québec produced a large geochemical dataset. Groundwater samples were collected at 153 sites over a 1500 km² study area and analyzed for major and minor ions. The large number of data can lead to difficulties in the integration, interpretation and representation of the results. Two multivariate statistical methods, hierarchical cluster analysis (HCA) and principal components analysis (PCA), were applied to a subgroup of the dataset to evaluate their usefulness to classify the groundwater samples, and to identify geochemical processes controlling groundwater geochemistry. This subgroup consisted of 144 samples and 14 parameters (Ca²⁺, Mg²⁺, Na⁺, K⁺, , Cl⁻, , Fe²⁺, Mn²⁺, Br⁻, Sr²⁺, F⁻, Ba²⁺, HS⁻). Seven geochemically distinct clusters, C1–C7, resulted from the HCA. Samples from clusters C3, C4, C6 and C7 are mostly located in preferential recharge areas. The majority of these samples have Ca–Mg–HCO₃ recharge groundwater (C3, C6, C7) and Na–HCO₃ evolved groundwater (C4). Samples from the other three clusters (C1, C2, C5) are characteristic of an aquifer system under confined conditions. The majority of these samples have Na–HCO₃ evolved groundwater (C1, C5) and Na–Cl ancient groundwater that exhibits elevated concentrations in Br⁻ (C2). In addition to recognizing the importance of hydrogeological conditions on groundwater geochemistry, the distribution of clusters also showed the importance of the geological formations on minor and trace elements, such as Fe²⁺, Mn²⁺, Sr²⁺, F⁻ and Ba²⁺. The first five components of the PCA account for 78.3% of the total variance in the dataset. Component 1 is defined by highly positive loadings in Na⁺, Cl⁻ and Br⁻ and is related to groundwater mixing with Champlain Sea water and solute diffusion from the marine clay aquitard. The high positive loadings in Ca²⁺ and Mg²⁺ of component 2 suggest the importance of dissolution of carbonate rocks in this aquifer system. From their characteristic loadings, the first two components are defined as the “salinity” and “hardness” components, respectively. Components 3–5 are related to more local and geological effects. The integration of the HCA and the PCA, with conventional classification of groundwater types, as well as with the hydrogeological and geological contexts, allowed the division of the region into four main geochemical areas, providing an improved regional picture of the aquifer system dynamics and hydrogeochemical evolution of groundwater. The following factors were recognized as influencing the evolution of groundwater identified in every geochemical area: (1) geological characteristics including sedimentary rock type and till mineralogy; (2) hydrogeological characteristics represented by the level of confinement and the hydraulic gradient; and (3) the geological history including the latest glaciation and the Champlain Sea invasion. With its integrated approach, this hydrogeochemical study contributes to the characterization and understanding of complex groundwater flow systems, and provides an example of the long-term geochemical evolution of hydrogeological systems after a major perturbation, in this case seawater invasion.     

辽宁省地震监测站地下水化学类型及成因分析

通过测试辽宁省15个地震监测站泉水或井水的氢、氧同位素组成及主要离子组分含量,讨论了泉水或井水的化学类型及其成因.测得泉水或井水的δD和δ~(18)O值范围分别为-82.5‰--54.4‰和-11.4‰—-7.3‰,表明所测泉水和井水的主要来源为大气降水.研究区低温泉水为低矿化度的Ca-SO_4·HCO_3型或Ca-HCO_3·SO_4型;而较高温度的花岗岩裂隙水中溶解了较多的钠硅酸盐矿物,水化学类型主要为Na-HCO_3·SO_4型;碳酸盐岩及含砾砂岩含水层分别分布于辽宁省西部及中部地区,水温略低于高温花岗岩裂隙水,水化学类型为Na·Ca-HCO_3型.水中F~-含量较高,且F~-含量与温度、pH值、Na~+和HCO_3~-的浓度呈正相关,表明泉水或井水的化学类型是深部富CO_2流体及大气降水与花岗岩反应的结果.     

The relationships between the velocities,attenuations and petrophysical properties of reservoir sedimentary rocks1

We have measured the velocities and attenuations of compressional and shear waves in 29 water-saturated samples of sandstones and shales at a confining pressure of 60 MPa and at frequencies of about 0.85 MHz. The measurements were made using a pulse echo method in which the samples (diameter 5 cm, length 1.5 cm to 2.5 cm) were placed between perspex buffer rods inside a high-pressure cell. The velocity of each seismic wave was determined from the traveltime difference of equivalent phase points (corrected for diffraction effects) of the signals reflected from the top and from the base of each sample. Attenuation was determined in a similar way by comparison of the diffraction corrected amplitudes of the signals. The attenuation data are presented as ‘quality factors’: Q_p and Q_s for compressional and shear waves respectively. The results show that Q_s is strongly correlated with V_s, that Q_p is weakly correlated with V_p, and that Q_p is strongly correlated with Q_s. Q_p is strongly dependent on the volume percentage of the assemblage of intra-pore minerals, whether they are clays or carbonates. It is concluded that the attenuation mechanism is due to the local fluid flow arising from the differential dilation of the solid rock frame and the intra-pore mineral assemblage, which is a result of their very different elastic moduli.