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.     

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.     

Assessment of frozen-ground conditions for engineering geology along the Qinghai–Tibet highway and railway, China

The Qinghai–Tibet Highway and Railway (the Corridor) across the Qinghai–Tibet Plateau traverses 670 km of permafrost and seasonally frozen-ground in the interior of the Plateau, which is sensitive to climatic and anthropogenic environmental changes. The frozen-ground conditions for engineering geology along the Corridor is complicated by the variability in the near-surface lithology, and the mosaic presence of warm permafrost and talik in a periglacial environment. Differential settlement is the major frost-effect problem encountered over permafrost areas. The traditional classification of frozen ground based on the areal distribution of permafrost is too generalized for engineering purposes and a more refined classification is necessary for engineering design and construction. A proposed classification of 51 zones, sub-zones, and sections of frozen ground has been widely adopted for the design and construction of foundations in the portion of the Corridor studied. The mean annual ground temperature (MAGT), near-surface soil types and moisture content, and active faults and topography are most commonly the primary controlling factors in this classification. However, other factors, such as local microreliefs, drainage conditions, and snow and vegetation covers also exert important influences on the features of frozen ground. About 60% of the total length of the Corridor studied possesses reasonably good frozen-ground conditions, which do not need special mitigative measures for frost hazards. However, other sections, such as warm and ice-rich or -saturated permafrost, particularly in the sections in wetlands, ground improvement measures such as elevated land bridges and passive or proactive cooling techniques need to be applied to ensure the long-term stability of thermally unstable, thick permafrost subsoils, and/or refill with non-frost-susceptible soils. Due to the long-history of the construction and management of the Corridor by various government departments, adverse impacts of construction and operation on the permafrost environment have been resulted. It is recommended that an integrated, executable plan for the routing of major construction projects within this transportation corridor be established and long-term monitoring networks installed for evaluating and mitigating the impact from anthropogenic and climatic changes in frozen-ground conditions.     

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.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

Sulfide-driven arsenic mobilization from arsenopyrite and black shale pyrite

We examined the hypothesis that sulfide drives arsenic mobilization from pyritic black shale by a sulfide-arsenide exchange and oxidation reaction in which sulfide replaces arsenic in arsenopyrite forming pyrite, and arsenide (As−1) is concurrently oxidized to soluble arsenite (As+3). This hypothesis was tested in a series of sulfide-arsenide exchange experiments with arsenopyrite (FeAsS), homogenized black shale from the Newark Basin (Lockatong formation), and pyrite isolated from Newark Basin black shale incubated under oxic (21% O₂), hypoxic (2% O₂, 98% N₂), and anoxic (5% H₂, 95% N₂) conditions. The oxidation state of arsenic in Newark Basin black shale pyrite was determined using X-ray absorption-near edge structure spectroscopy (XANES). Incubation results show that sulfide (1 mM initial concentration) increases arsenic mobilization to the dissolved phase from all three solids under oxic and hypoxic, but not anoxic conditions. Indeed under oxic and hypoxic conditions, the presence of sulfide resulted in the mobilization in 48 h of 13-16 times more arsenic from arsenopyrite and 6-11 times more arsenic from isolated black shale pyrite than in sulfide-free controls. XANES results show that arsenic in Newark Basin black shale pyrite has the same oxidation state as that in FeAsS (−1) and thus extend the sulfide-arsenide exchange mechanism of arsenic mobilization to sedimentary rock, black shale pyrite. Biologically active incubations of whole black shale and its resident microorganisms under sulfate reducing conditions resulted in sevenfold higher mobilization of soluble arsenic than sterile controls. Taken together, our results indicate that sulfide-driven arsenic mobilization would be most important under conditions of redox disequilibrium, such as when sulfate-reducing bacteria release sulfide into oxic groundwater, and that microbial sulfide production is expected to enhance arsenic mobilization in sedimentary rock aquifers with major pyrite-bearing, black shale formations.     

Origin of Fe-Ti Oxide Ores in Mafic Intrusions: Evidence from the Panzhihua Intrusion, SW China

Economic concentrations of Fe–Ti oxides occur as massive,conformable lenses or layers in the lower part of the Panzhihuaintrusion, Emeishan Large Igneous Province, SW China. Mineralchemistry, textures and QUILF equilibria indicate that oxidesin rocks of the intrusion were subjected to extensive subsolidusre-equilibration and exsolution. The primary oxide, reconstructedfrom compositions of titanomagnetite in the ores and associatedintergrowths, is an aluminous titanomagnetite (Usp₄₀) with 40wt % FeO, 34 wt % Fe₂O₃, 16·5 wt % TiO₂, 5·3 wt% Al₂O₃, 3·5 wt % MgO and 0·5 wt % MnO. This compositionis similar to the bulk composition of the oxide ore, as inferredfrom whole-rock data. This similarity strongly suggests thatthe ores formed from accumulation of titanomagnetite crystals,not from immiscible oxide melt as proposed in earlier studies.The occurrence of oxide ores in the lower parts of the Panzhihuaintrusion is best explained by settling and sorting of densetitanomagnetite in the ferrogabbroic parental magma. This magmamust have crystallized Fe–Ti oxides relatively early andabundantly, and is likely to have been enriched in Fe and Tibut poor in SiO₂. These features are consistent with fractionationof mantle-derived melts under relatively high pressures (10kbar), followed by emplacement of the residual magma at 5 kbar.This study provides definitive field and geochemical evidencethat Fe–Ti oxide ores can form by accumulation in ferrogabbro.We suggest that many other massive Fe–Ti oxide depositsmay have formed in a similar fashion and that high concentrationsof phosphorus or carbon, or periodic fluctuation of fO₂ in themagma, are of secondary importance in ore formation. KEY WORDS: ELIP; Fe–Ti oxide ore; layered intrusion; Panzhihua; QUILF     

The state of land subsidence and prediction approaches due to groundwater withdrawal in China

This article gives a general introduction to land subsidence with the prediction approaches due to withdrawal of groundwater in three subsided/subsiding regions in China: the deltaic plain of Yangtse River (YRDP), North China Plain (NCP), and Fenwei Plain (FP). On YRDP, Shanghai is the typical subsided/subsiding city; on NCP Tianjin is the typical subsided/subsiding city, and on FP Taiyuan is the typical subsided/subsiding city. The subsided area with subsidence over 200 mm on YRDP is about 10,000 km² and the maximum subsided value reached 2.9 m at Shanghai; on NCP the subsided area reached 60,000 km² with the maximum subsidence of 3.9 m at Tianjing; on FP the subsided area is relatively smaller than that on the other two plains and is about 1,135 km² with maximum subsidence of 3.7 m at Taiyuan city. In order to protect the civil and industrial facilities, it is necessary to predict the future development of land subsidence based on present state. Many researchers proposed several approaches to predict the land subsidence due to groundwater withdrawal according to different geological conditions and groundwater withdrawal practice. This article classifies these approaches into five categories: (i) statistical methods; (ii) 1D numerical method; (iii) quasi-3D seepage model; (iv) 3D seepage model; (v) fully coupled 3D model. In China, the former four categories are presently employed in the prediction practice and their merits and demerits are discussed. According to the prediction practice, 3D seepage model is the best method presently.