降雨入渗补给地下水数值模拟研究

针对在地下水资源评价中入渗补给量采用经验估算存在较大误差的问题,以神东矿区为例,在广泛调查矿区包气带岩性结构并结合野外取样、室内参数测定的基础上,采用数值模拟的方法建立包气带水分运移数值模拟模型,定量模拟矿区内不同地段降雨入渗强度,探讨影响降雨入渗强度的主要因素,计算得出研究区降雨入渗系数大致在0.18~0.27,分析认为影响降雨入渗强度的因素有降雨量、潜水埋深、包气带岩性等。其中在研究区广泛分布的风积沙对地下水资源起到了一定的保护作用。     

Predicting lake water quality responses to load reduction: a three-dimensional modeling approach for total maximum daily load

Water quality restoration efforts often suffer the risk of ineffectiveness and failure due to lack of quantitative decision supports. During the past two decades, the restoration of one of China’s most heavily polluted lakes, Lake Dianchi, has experienced costly decision ineffectiveness with no detectable water quality improvement. The governments are planning to invest tremendous amount of funds in the next 5 years to continue the lake restoration process; however, without a quantitative understanding between the load reduction and the response in lake water quality, it is highly possible that these planned efforts would suffer the similar ineffectiveness as before. To provide scientifically sound decision support for guiding future load reduction efforts in Lake Dianchi Watershed, a sophisticated quantitative cause-and-effect response system was developed using a three-dimensional modeling approach. It incorporates the complex three dimensional hydrodynamics, fate and transport of nutrients, as well as nutrient-algae interactions into one holistic framework. The model results show that the model performs well in reproducing the observed spatial pattern and temporal trends in water quality. The model was then applied to three total maximum daily load scenarios and two refined restoration scheme scenarios to quantify phytoplankton responses to various external load reduction intensities. The results show that the algal bloom in Lake Dianchi responds to load reduction in a complex and nonlinear way, therefore, it is necessary to apply the developed system for future load reduction and lake restoration schemes for more informed decision making and effective management.     

A method for mineral prospectivity mapping integrating C4.5 decision tree,weights-of-evidence and m-branch smoothing techniques: a case study in the eastern Kunlun Mountains,China

In this study, a novel method that integrates C4.5 decision tree, weights-of-evidence and m-branch smoothing techniques was proposed for mineral prospectivity mapping. First, a weights-of-evidence model was used to rank the importance of each evidential map and determine the optimal buffer distance. Second, a classification technique that uses a C4.5 decision tree in data mining was used to construct a decision tree classifier for the grid dataset. Finally, an m-branch smoothing technique was used as a predictor, which transformed the decision tree into a probability evaluation tree. The method makes no conditional independence assumption and can be applied for class imbalanced datasets like those collected during mineral exploration for prospectivity mapping of an area. The traits of comprehensibility, accuracy and efficiency were derived from the C4.5 decision tree. In addition, a case study for iron prospectivity mapping was performed in the eastern Kunlun Mountains, China. Sixty-two Skarn iron deposits and eight evidential maps related to iron mineralization were studied. On the final map, areas of low, moderate and high potential for iron deposit occurrence covered areas of 71,491, 14,298, and 9,532 km², respectively. For the goodness-of-fit test, 91.94 % of the total 62 iron deposits were within a high-potential area, 8.06 % were within a moderate-potential area and 0 % were within a low-potential area. For ten-fold cross-validation, 82.26 % were within a high-potential area, 14.52 % were within a moderate-potential area and 3.22 % were within a low-potential area. To evaluate the predictive accuracy, Receiver Operating Characteristic (ROC) curves and Area Under ROC Curve (AUC) were employed. The accuracy of the goodness-of-fit test reached 97.07 %, and the accuracy of the ten-fold cross-validation was 95.10 %. The majority of the iron deposits were within high-potential and moderate-potential areas, which covered a small proportion of the study area.     

Defining physico-chemical variables, chlorophyll-a and Secchi depth reference conditions in northeast eco-region lakes, China

Establishing reference conditions for regional lakes is necessary to assess human impact on aquatic ecosystems, protect water quality and biotic integrity. The northeast eco-region lakes often present a marked seasonal variability in hydrological, biological and geochemical processes, which could affect physico-chemical variables, chlorophyll-a and Secchi depth reference conditions. Reference conditions were calculated using two different methods: lake population distribution method and trisection method. General physico-chemical reference conditions (pH, conductivity, BOD and nutrients) and chlorophyll-a, and Secchi depth reference conditions were established using the two approaches combined with the analysis of seasonal variability for each parameter for shallow and deep lake types in the northeast eco-region lakes. The results indicated that depth was an important factor affecting physico-chemical variables, chlorophyll-a and Secchi depth reference values. Moreover, the seasonal differences in parameters were also assessed using a non-parametric repeated measure ANOVA (Friedman test) in the shallow and deep lakes. Some parameters showed a seasonal variability for shallow and deep lakes. The least restrictive reference values of the four seasons were proposed to simplify the water quality assessment for regional lake management goals. These results indicated that the 50th percentile (median) value from the best one-third of the nutrient-concentration data distributions is more appropriate for the northeast eco-region lakes.     

Controlling factors of the δ 11B-pH proxy and its research direction

Significant boron isotope fractionation occurs in nature (?70 ‰ to +75 ‰) due to the high geochemical reactivity of boron and the large relative mass difference between ¹⁰B and ¹¹B. Since the 1990s, reconstruction of ancient seawater pH using the isotopic composition of boron in bio-carbonates (δ ¹¹B_carb), and then calculation of the past pCO₂ have become important issues for the international isotope geochemistry community, and are called the δ ¹¹B-pH proxy. Although many achievements have been made by this proxy, various aspects of boron systematics require rigorous evaluation. Based on the previous researches, mechanism of boron isotope fractionation, variation of boron isotope (δ ¹¹B) in nature (especially in bio-carbonates) and controlling factors of the δ ¹¹B-pH proxy, such as the dissociation constant of B(OH)₃ in seawater (pK_a), the δ ¹¹B of seawater (δ ¹¹B_SW), the boron isotopic fractionation factor between B(OH) ₄ ^? and B(OH)₃ (α _4–3), and the incorporated species of boron into bio-carbonates, are reviewed in detail and the research directions of this proxy are proposed. Generally, the controversy about pK_a, δ ¹¹B_sw, and α _4–3 is relatively less, but whether boron incorporated into bio-carbonates only in the form of B(OH) ₄ ^? remains doubtful. In the future, it is required that the physicochemical processes that control boron incorporation into carbonates be rigorously characterized and that the related chemical and isotopic fractionation be quantified. It is also necessary and important to establish a “best-fit empirically equation” between δ ¹¹B_carb and pH of seawater based on the precipitation experiments of inorganic or culture experiments of corals or foraminifera. In addition, extended application of the δ ¹¹B-pH proxy to the earlier part of the Phanerozoic relying on the Brachiopods is worthy of studying. Like other geochemical indicators, there are limiting factors of δ ¹¹B; however, it remains a very powerful tool in the reconstruction of past seawater pH at present.     

Antimony adsorption on kaolinite in the presence of competitive anions

Antimony (Sb) emissions to the environment are increasing, and there is a dearth of knowledge regarding Sb fate and behavior in natural systems. In natural systems, the presence of competitive anions may compete with Sb for adsorption sites on mineral surfaces, hence increasing its potential bioavailability. Accordingly, the adsorption of Sb(III) on kaolinite was investigated in the presence of competitive anions. Kinetic studies suggest that adsorption reaction of Sb(III) on kaolinite is rapid initially and becoming slow after 12 h both in binary Sb(III)–kaolinite system and in ternary Sb(III)-competitive anion–kaolinite system. The presence of PO₄ ^3? has a much stronger and more obvious promotive effect on the adsorption of Sb(III) on kaolinite compared with the other two anions. The adsorption data of Sb(III) on kaolinite in the absence and presence of competitive anions at three temperatures were successfully modeled using Langmuir (r ² > 0.95) and Freundlich (r ² > 0.95) isotherms. Accompanied the adsorption of Sb(III) on kaolinite, significant oxidation of Sb(III) to Sb(V) had occurred under the experimental conditions used in this study. The presence of kaolinite which has a larger specific surface area could increase the contact area between the adsorbed Sb(III) and oxygen in the bulk solution, which promoted the oxidation rate of Sb(III) to Sb(V).     

Decoupling of the East Asian summer monsoon and Indian summer monsoon between 20 and 17 ka

Marine Oxygen Isotope Stage (MIS) 2, with its profound environmental and climatic changes from before the last glacial maximum (LGM) to the last deglaciation, is an ideal period for understanding the evolution of the East Asian summer monsoon (EASM) and Indian summer monsoon (ISM), two Asian monsoon sub-systems. With 875 stable oxygen isotope ratios and 43 ²³⁰Th dates from stalagmites in Sanxing Cave, southwestern China, we construct and interpret a new, replicated, Asian summer monsoon (ASM) record covering 30.9–9.7 ka with decadal resolution. δ¹⁸O records from this site and other reported Chinese caves display similar long-term orbitally dominated trends and synchronous millennial-scale strong and weak monsoonal events associated with climate changes in high northern latitudes. Interestingly, Sanxing δ¹⁸O and Arabian Sea records show a weakening ISM from 22 to 17 ka, while the Hulu and Qingtian records from East and Central China express a 3-ka intensifying EASM from 20 to 17 ka. This decoupling between EASM and ISM may be due to different sensitivities of the two ASM sub-systems in response to internal feedback mechanisms associated with the complex geographical or land-ocean configurations.     

Changes in the critical erosion velocity for sediment colonized by biofilm

In flowing water the incipient motion of sediment can be affected by the presence of microbial biofilm growth. This article documents a series of flume experiments using non‐uniform sediments, in which sediment entrainment was investigated for cases where the sediment was immersed in deionized water, so that no biofilm developed, and for cases where a bio‐sediment was cultivated by placing the sediment in a mixture of natural water and nutrient solution. Differences in entrainment and the velocity at incipient motion were measured over an eight week period, as biofilm grew. It was found that the incipient motion phenomena were quite distinct between the two kinds of sediment. Sediment with biofilm was more stable and, over time, incipient velocity increased to a threshold level, before declining. Biofilm development is clearly an important control on the stability of sediments, especially in eutrophic water bodies. Two incipient velocity formulas were derived for sliding and rolling conditions. Film water theory was utilized to describe the cohesive force between sediment particles and the adhesive force generated by biofilm was introduced into the formula derivation; the time variation characteristics of biofilm strength and the features of the substrate were also taken into consideration. Such analyses can help to predict sediment transport changes due to biofilm presence in nutrient‐rich water bodies.     

塔中地区奥陶系喜马拉雅山期S干气充注的证据

塔中地区奥陶系天然气成因多样;Ⅰ号坡折带中东部奥陶系天然气以高干燥系数、 甲烷同位素值重为特征;与塔深1井寒武系原油裂解气接近;应主要来自寒武系原油裂解气成因。寒武系贫H₂S、 高成熟原油裂解气在喜马拉雅山期时;气侵奥陶系油气藏;得到了以下主要证据的支持: 1）天然气甲烷δ¹³C值大多比Chung et al.（1988）天然气模式甲烷δ¹³C值计算值高3‰以上;2）干燥系数与甲烷δ¹³C值大体上具有正相关关系;3）天然气干燥系数与H₂S含量大体上具有负相关关系。这些特征表明;存在贫H₂S、 相对富¹³C甲烷为主的干气与富H₂S、 相对贫¹³C甲烷的湿气混合作用。奥陶系中H₂S-δ³⁴S 值为14‰～20‰;远低于中深1井寒武系原地热化学硫酸盐还原作用（TSR）成因的H₂S（33‰）;支持了奥陶系中H₂S并不是来源于寒武系古油气藏。于是提出;来自寒武系贫H₂S的干气在喜马拉雅山期对良里塔格组和鹰山组油气藏发生了气洗;油气藏的气/油比值增大、 导致了原油蜡含量增高、 甲烷δ¹³C值发生正偏移。     

中国地热研究的进展与展望(1995～2014)

地热是主要地球物理场之一,地热能是地球的本土能源。李四光在20世纪60年代开创了我国地热科学。到90年代,学科体系基本建立。在过去20年里,我国地热研究得到了进一步深化和拓展。本文从大地热流、岩石圈热结构、地热系统、油气盆地地热、矿山地热、天然气水合物以及气候变化等方面回顾地热研究代表性的创新进展,并对深层地热、海洋地热、环境地热等研究方向作了展望。本文认为,过去20年我国地热研究成果丰硕,国际影响力得到提高,未来发展势头强劲。在经历了由浅入深,从今到古的成长之后,地热研究还将不断拓展领域,为我国地球科学,特别是能源与环境安全做出更大贡献。