Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.     

Dehydration and melting during continental collision: Constraints from element and isotope geochemistry of low-T/UHP granitic gneiss in the Dabie orogen

A combined study of petrography, whole-rock major and trace elements as well as Rb?Sr and Sm?Nd isotopes, and mineral oxygen isotopes was carried out for two groups of low-T/UHP granitic gneiss in the Dabie orogen. The results demonstrate that metamorphic dehydration and partial melting occurred during exhumation of deeply subducted continent. Zircon δ¹⁸O values of ? 2.8 to + 4.7‰ for the gneiss are all lower than normal mantle values of 5.3 ± 0.3‰, consistent with ¹⁸O depletion of protolith due to high-T meteoric-hydrothermal alteration at mid-Neoproterozoic. Most samples have extremely low ⁸⁷Sr/⁸⁶Sr ratios at t₁ = 780 Ma, but very high ⁸⁷Sr/⁸⁶Sr ratios at t₂ = 230 Ma. This suggests intensive fluid disturbance due to the hydrothermal alteration of protoliths during Neoproterozoic magma emplacement and the metamorphic dehydration during Triassic continental collision. Rb–Sr isotopes, Th/Ta vs. La/Ta and Th/Hf vs. La/Nb relationships suggest that Group I gneiss experienced lower degrees of hydrothermal alteration, but higher degrees of dehydration, than Group II gneiss. The two groups of gneiss have similar patterns of REE and trace element partition. Group I gneiss displays good correlations between Nb and LREEs but no correlations between Nb and LILEs (Rb, Ba, Pb, Th and U), indicating differential mobilities of LILEs during the dehydration. Thus the correlation between Nb and LREEs is inherited from protolith rather than caused by metamorphic modification. Relative to Group I gneiss, Group II gneiss has stronger negative Eu anomaly, lower contents of Sr and Ba but higher contents of Rb, Th and U. In particular, Nb correlates with LILEs (e.g., Rb, Sr, Ba, Th and U), but not with LREEs (La and Ce). This may indicate decoupling between the dehydration and LILEs transport during continental collision. Furthermore, dehydration melting may have occurred due to breakdown of muscovite during “hot” exhumation. Group II gneiss has extremely low contents of FeO + MgO + TiO₂ (1.04 to 2.08 wt.%), high SiO₂ contents of 75.33 to 78.23 wt%, and high total alkali (Na₂O + K₂O) contents (7.52 to 8.92 wt.%), comparable with compositions predicted from partial melting of felsic rocks by experimental studies. Almost no UHP metamorphic minerals survived; felsic veins of fine-grain minerals occurs locally between coarse-grain minerals, resulting in a kind of metatexite migmatites due to dehydration melting without considerable escape of felsic melts from the host gneiss. In contrast, Group I gneiss only shows metamorphic dehydration. Therefore, the two groups of gneiss show contrasting behaviors of fluid–rock interaction during the continental collision.     

Estimation of groundwater recharge using the soil moisture budget method and the base-flow model

A detailed characterization of the site is crucial to designing an efficient method of managing the risks associated with tailings from abandoned mines. Therefore, samples collected from various depths within tailings in Guryong mine, Korea, were analyzed for their chemical, physical and mineralogical characteristics. All samples of the Guryong tailings had acid-generating potential. However, in the oxidation zone, the net acid generation (NAG) was low (30 kg H₂SO₄ t⁻¹) although the acid neutralization potential (ANP) was less than zero. The ANP values in the unoxidation zone were higher (> −56.0 kg CaCO₃ t⁻¹) than in the other zones. As a result, the amount of alkali ions that are needed to neutralize the acid needs to be considered. In this experiment G3, G4 and G6 drill cores containing fine tailings particles near the unoxidation zone were observed to contain calcite (CaCO₃) with acid-neutralizing capacity. A low pH (2−4) in the oxidation zone of the tailings changed to a neutral pH in the unoxidation zone of the tailings. These results suggest that the acid-neutralizing capacity of the tailings was controlled by particle and mineral composition of tailings.     

Occurrence and evolution of the Xiaotangshan hot spring in Beijing, China

Thermal groundwater occurs in bedrock aquifers consisting of the dolomite of the Wumishan Group of the Jixianin System and the Cambrian carbonate in the Xiaotangshan geothermal field near the northern margin of the North China Plain, China. The hot water in the geothermal field of basin-type discharges partly in the form of the Xiaotangshan hot spring under natural conditions. The hot water has TDS of less than 600 mg/L and is of Na·Ca-HCO₃ type. The geothermal water receives recharge from precipitation in the mountain area with elevation of about 500 m above sea level to the north of the spring. Thermal groundwater flows slowly south and southeast through a deep circulation with a residence time of 224 years estimated with the Ra–Rn method. The Xiaotangshan hot spring dried up in the middle of the 1980s owing to the increasing withdrawal of the hot water in the geothermal field in the past decades. The water level of the geothermal system still falls continually at an annual average rate of about 2 m, although water temperature changes very little, indicating that the recharge of such a geothermal system of basin-type is limited. Over-exploitation has a dramatic impact on the geothermal system, and reduction in exploitation and reinjection are required for the sustainable usage of the hot water.     

Response of groundwater chemistry to water deliveries in the lower reaches of Tarim River,Northwest China

In this paper, we analysed the monitored data from nine groundwater-monitoring transects in the lower reaches of Tarim River during the five times of stream water deliveries to the river transect where the stream flow ceased. The results showed that the groundwater depth in the lower reaches of Tarim River rose from −9.30 m before the conveyances to −8.17 and −6.50 m after the first and second conveyances, −5.81 and −6.00 m after the third and fourth the conveyance, and −4.73 m after the fifth. The horizontal extent of groundwater recharge was gradually enlarged along both sides of the channel of conveyance, i.e., from 250 m in width after the first conveyance to 1,050 m away from the channel after the fourth delivery. With the rising groundwater level, the concentrations of major anions Cl⁻, SO₄²⁻ and cations Ca²⁺, Mg²⁺, Na⁺, as well as total dissolved solids (TDS) in groundwater underwent a significant change. The spatial variations in groundwater chemistry indicated that the groundwater chemistry at the transect near Daxihaizi Reservoir changed earlier than that farther from it. In the same transect, the chemical variations were earlier in the monitoring well close to watercourse than that farther away from the stream. In general, the concentration of the major ions and TDS at each monitoring well increased remarkably when the water delivery started, and decreased with the continued water delivery, and then increased once again at the end of the study period. Hence, the whole study period may be divided into three stages: the initial stage, the intermediate stage and the later stage. According to the three stages of groundwater chemistry reaction to water delivery and the relationships between groundwater chemical properties and groundwater depths, we educe that under the situation of water delivery, the optimum groundwater depth in the lower reaches of the Tarim River should be −5 m.     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

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.     

Evolution of subcontinental lithospheric mantle beneath eastern China: Re–Os isotopic evidence from mantle xenoliths in Paleozoic kimberlites and Mesozoic basalts

The ages of subcontinental lithospheric mantle beneath the North China and South China cratons are less well-constrained than the overlying crust. We report Re–Os isotope systematics of mantle xenoliths entrained in Paleozoic kimberlites and Mesozoic basalts from eastern China. Peridotite xenoliths from the Fuxian and Mengyin Paleozoic diamondiferous kimberlites in the North China Craton give Archean Re depletion ages of 2.6–3.2 Ga and melt depletion ages of 2.9–3.4 Ga. No obvious differences in Re and Os abundances, Os isotopic ratios and model ages are observed between spinel-facies and garnet-facies peridotites from both kimberlite localities. The Re–Os isotopic data, together with the PGE concentrations, demonstrate that beneath the Archean continental crust of the eastern North China Craton, Archean lithospheric mantle of spinel- to diamond-facies existed without apparent compositional stratification during the Paleozoic. The Mesozoic and Cenozoic basalt-borne peridotite and pyroxenite xenoliths, on the other hand, show geochemical features indicating metasomatic enrichment, along with a large range of the Re–Os isotopic model ages from Proterozoic to Phanerozoic. These features indicate that lithospheric transformation or refertilization through melt-peridotite interaction could be the primary mechanism for compositional changes during the Phanerozoic, rather than delamination or thermal-mechanical erosion, despite the potential of these latter processes to play an important role for the loss of garnet-facies mantle. A fresh garnet lherzolite xenolith from the Yangtze Block has a Re depletion age of ∼1.04 Ga, much younger than overlying Archean crustal rocks but the same Re depletion ages as spinel lherzolite xenoliths from adjacent Mesozoic basalts, indicating Neoproterozoic resetting of the Re–Os system in the South China Craton.