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.     

Distribution of Cd,Pb, Zn and Cu and their chemical speciations in soils from a peri-smelter area in northeast China

An exploratory study on soil contamination of heavy metals was carried out surrounding Huludao zinc smelter in Liaoning province, China. The distribution of total heavy metals and their chemical speciations were investigated. The correlations between heavy metal speciations and soil pH values in corresponding sites were also analyzed. In general, Cd, Zn, Pb, Cu and As presented a significant contamination in the area near the smelter, comparied with Environmental Quality Standards for Soils in China. The geoaccumulation index showed the degree of contamination: Cd > Zn > Pb > Cu > As. There was no obvious pollution of Cr and Ni in the studied area. The speciation analysis showed that the dominant fraction of Cd and Zn was the acid soluble fraction, and the second was the residual fraction. Pb was mostly associated with the residual fraction, which constituted more than 50% of total concentration in all samples. Cu in residual fraction accounted for a high percentage (40–80%) of total concentration, and the proportion of Cu in the oxidizable fraction is higher than that of other metals. The distribution pattern of Pb and Zn was obviously affected by soil pH. It seemed that Pb and Zn content in acid solution fraction increased with increasing soil pH values, while Cd content in acid soluble fraction accounted for more proportion in neutral and alkaline groups than acidic one. The fraction distribution patterns of Cu in three pH groups were very similar and independent of soil pH values. And the residual fraction of Cu took a predominant part (50%) of the total content.     

Glacier change and glacier runoff variation in the Tuotuo River basin,the source region of Yangtze River in western China

Glaciers in the Tuotuo River basin, western China, have been monitored in recent decades by applying topographical maps and high-resolution satellite images. Results indicate that most of glaciers in the Tuotuo River basin have retreated in the period from 1968/1971 to 2001/2002, and their shrinkage area is 3.2% of the total area in the late 1960s. To assess the influence of glacier runoff on river runoff, a modified degree–day model including potential clear-sky direct solar radiation has been applied to the glaciated regions of the river basin over the period 1961–2004. It was found that glacier runoff has increased in the last 44 years, especially in the 1990s when a two-thirds increase in river runoff was derived from the increase in glacier runoff caused by loss of ice mass in the entire Tuotuo River basin.     

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.     

Experimental study on the deterioration and natural remediation of the Ariake Sea tidal mud caused by the sea laver treatment acid practice and the upward seepage of pore water liquid

This study presents a laboratory study of the following two aspects: (1) the influence of sea laver treatment acid on the geoenvironmental properties of Ariake Sea tidal mud, and (2) the natural remediation effect on the sea laver treatment acid contaminated Ariake Sea tidal mud caused by the upward seepage of pore water liquid in the mud. Firstly, the mechanisms of the transport of sea laver treatment acid in the Ariake Sea tidal mud and the generation mechanisms of the upward seepage flow in the Ariake Sea tidal mud are discussed. Secondly, a series of one-dimensional laboratory infiltration tests were carried out to investigate the deterioration of the Ariake Sea tidal mud caused by the sea laver acid treatment practice. Test results reveal that the acid treatment practice caused considerable change in the geochemical properties of the mud in terms of increase in sulfide content and decrease in pH value. After the treatment by the sea laver treatment acid, the sulfide content of the mud even exceeded the safe limit value for the benthos, which represents undesirable living condition for benthos. Thirdly, series of laboratory fresh seawater infiltration tests for the deteriorated Iida site mud were conducted to illustrate this natural remediation efficiency. It is found that with the infiltration of the fresh seawater, the sulfide content of the Iida site mud was considerably reduced and pH value increased to an acceptable range for benthos living in the tidal flat mud. With the increase in the infiltration time and the hydraulic gradient, the remediation efficiency could be increased.     

Triassic Nb-enriched basalts, magnesian andesites, and adakites of the Qiangtang terrane (Central Tibet): evidence for metasomatism by slab-derived melts in the mantle wedge

New chronological, geochemical, and isotopic data are reported for Triassic (219–236 Ma) adakite-magnesian andesite-Nb-enriched basaltic rock associations from the Tuotuohe area, central Qiangtang terrane. The adakites and magnesian andesites are characterized by high Sr/Y (25–45), La/Yb (14–42) and Na₂O/K₂O (12–49) ratios, high Al₂O₃ (15.34–18.28 wt%) and moderate to high Sr concentrations (220–498 ppm) and ε_ND (t) (+0.86 to +1.21) values. Low enrichments of Th, Rb relative to Nb, and subequal normalized Nb and La contents, and enrichments of light rare earth elements combine to distinguish a group of Nb-enriched basaltic rocks (NEBs). They have positive ε_ND (t) (+2.57 to +5.16) values. Positive correlations between Th, La and Nb and an absence of negative Nb anomalies on mantle normalized plots indicate the NEBs are products of a mantle source metasomatized by a slab melt rather than by hydrous fluids. A continuous compositional variation between adakites and magnesian andesites confirms slab melt interaction with mantle peridotite. The spatial association of the NEBs with adakites and magnesian andesites define an “adakitic metasomatic volcanic series” recognized in many demonstrably subduction-related environments (e.g., Mindanao arc, Philippines; Kamchatka arc, Russia; and southern Baja California arc, Mexico). The age of the Touhuohe suite, and its correlation with Triassic NEB to the north indicates that volcanism derived from subduction-modified mantle was abundant prior to 220 Ma in the central Qiangtang terrane.     

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.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

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