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.     

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.     

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.     

苏北盘石山、练山地幔捕虏体的PGE地球化学

通过锍镍火试金预富集法,分析了位于郯庐断裂带东侧的盘石山、练山地幔橄榄岩包体中铂族元素(PGE)和Au含量.不同于部分熔融残留成因地幔橄榄岩中通常所观察到的负斜率型或平坦型的分布模式,这两地的地幔橄榄岩以Pt、Pd、Ru相对富集,Ir、Rh相对亏损的"燕子型"分布模式为特征.Pt、Pd等不相容元素富集说明上地幔除经历过早期的部分熔融外,还经历了后期富Pt、Pd的高熔/岩比的熔(流)体的层析分离交代作用影响.盘石山地幔橄榄岩的PGE总量比练山高,Os的含量也比原始地幔值高;而练山地幔橄榄岩的Os含量比原始地幔值低,说明交代作用带走了练山地幔橄榄岩中的Os,却没有很大改变盘石山地幔橄榄岩中的Os含量,这可能与交代熔(流)体含硫量饱和程度有关.Rh的负异常可能与部分熔融过程中熔体较低的fo2有关.     

黔东震旦系—下寒武统黑色岩系稀土元素地球化学特征

：通过对黔东丹寨南皋剖面和三都渣拉沟剖面的稀土元素分析发现：两条剖面的梅树村期与筇竹寺期 分界处大量稀土元素明显富集,黑色页岩中存在明显的Ｃｅ负异常和Ｅｕ、Ｙ 正异常,表明为缺氧和热水沉积的产物。 这套黑色岩系总体沉积于缺氧和具热水注入的环境中,但发生缺氧和热水注入的时间和强度在不同地方会有不同 的表现。     

邹平王家庄铜矿床成矿地球化学及成因探讨

王家庄铜矿床的矿化脉石英中流体包裹体均一温度介于116 ~ 566 ℃之间,均值为 289 ℃;盐度介于7.2%～63.2% NaCleq,均值为21.1% NaCleq。流体的气相成分主要为H2O和CO2。在均一温度为240 ~ 440 ℃区间内,出现了富气相的两相水溶液包裹体、富液相的两相水溶液包裹体和含子晶的三相水溶液包裹体共存现象,以及加温后富气相包裹体均一到气相和同期富液相包裹体均一到液相的特征,表明成矿流体曾发生过沸腾作用;其中第一次发生于360 ~ 400 ℃,主要形成高温、高盐度含子晶的三相水溶液包裹体和高温、中盐度富液相的两相水溶液包裹体及高温、低盐度富气相的两相水溶液包裹体;第二次发生于240 ~ 320 ℃,主要形成高—中温、高盐度的含子晶的三相水溶液包裹体和高—中温、中盐度富液相的两相水溶液包裹体及高—中温、低盐度富气相的两相水溶液包裹体;之后主要形成富液相的两相水溶液包裹体,具有中低温和中盐度的特征。矿化脉石英中的δ18OH2O介于-1.14‰ ~ 1.79‰之间,均值为0.94‰;δD介于-63.70‰ ~ -56.50‰之间,均值为-59.8‰;说明王家庄铜矿床的成矿流体主要来源于岩浆,后期混入大气降水。矿石矿物的δ34S变化于-8.80‰ ~ -2.80‰之间,均值为-6.33‰。矿石矿物的n(206Pb)/n(204Pb)介于18.1684 ~ 18.3637之间,均值为18.2892;n(207Pb)/n(204Pb)介于155546 ~ 156342之间,均值为155777;n(208Pb)/n(204Pb)介于38.1286 ~ 38.4840之间,均值为38.2780。说明矿石具有贫重硫和富放射性成因铅的特征,硫、铅主要来源于深部,后期可能受到地壳物质或大气降水的混染。