Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

《干旱区地理》的引文统计分析

本文使用文献计量学中的引文分析方法,以《干旱区地理》为调查样本进行了统计分析。《干旱地理》１９８５－１９９４年共版４３期,刊登论文５５８篇,其中附有引文的论文４３６篇,引文２７６６条。     

处理水环境中微量铬的菌种筛选研究

从云南各地采集微生物菌种进行筛选、驯化,找到一株对水环境中的微量铬具有良好吸附去除能力的菌种,并研究了碳源、氮源、培养液pH值和培养温度等因素对其去除铬的影响.试验结果表明,去除水环境中微量铬的适宜菌种的培养条件为以可溶性淀粉为碳源、（NH4）2SO4为氮源、培养温度25℃、培养液pH 8,培养液中碳：氮=10：1,利用该条件培养的菌种对水环境中微量铬具有较好的去除能力.     

岩溶区集中式饮用水源地水质评价与保护——以毕节市倒天河水库和利民水库为例

采用SPSS软件，对2003年01月至2004年12月期间，倒天河水库和利民水库逐月水质监测数据进行数理统计分析，结果表明毕节市饮用水水质污染有较明显的季节变化和空间差异；且不同污染物对水体污染的贡献率各不相同，总磷和粪大肠菌群所占份额最大，方差累积贡献率达77％。在此基础上，利用加权综合污染指数对水体水质污染现状进行综合评价，结果显示倒天河水库和利民水库，水质污染总体较轻．饮用水水质符合《地表水环境质量标准（GB3838—2002）》中Ⅲ类水质标准；最后提出加大水环境保护的建议。     

Study on the kinetics of iron oxide leaching by oxalic acid

The presence of iron oxides in clay or silica raw materials is detrimental to the manufacturing of high quality ceramics. Although iron has been traditionally removed by physical mineral processing, acid washing has been tested as it is more effective, especially for extremely low iron (of less than 0.1% w/w). However, inorganic acids such as sulphuric or hydrochloric acids easily contaminate the clay products with SO₄²⁻ and Cl⁻, and therefore should be avoided as much as possible. On the other hand, if oxalic acid is used, any acid left behind will be destroyed during the firing of the ceramic products. The characteristics of dissolution of iron oxides were therefore investigated in this study.The dissolution of iron oxides in oxalic acid was found to be very slow at temperatures within the range 25–60 °C, but its rate increases rapidly above 90 °C. The dissolution rate also increases with increasing oxalate concentration at the constant pH values set within the optimum range of pH2.5–3.0. At this optimum pH, the dissolution of fine pure hematite (Fe₂O₃) (105–140 μm) follows a diffusion-controlled shrinking core model. The rate expression expressed as 1 − (2 / 3)x − (1 − x)^2 / 3 where x is a fraction of iron dissolution was found to be proportional to [oxalate]^1.5.The addition of magnetite to the leach liquor at 10% w/w hematite was found to enhance the dissolution rate dramatically. Such addition of magnetite allows coarser hematite in the range 0.5–1.4 mm to be leached at a reasonable rate.     

CQK-50型车载化钻机的研制

在分析研究小型地质车载钻机的形式、功能与特点的基础上，提出了钻机车载化的新思路。根据地质调查工作现状，结合国内现有小型地质取样钻机，研制了CQK-50型全液压多功能车载取样钻机，并在野外物化探遥感异常验证钻探中取得了良好效果。     

遥感技术在煤炭地质中的应用现状及前景

介绍了遥感技术在煤炭资源调查和评价、煤矿区地质灾害调查和评估、生态环境及污染调查评价和监测以及煤矿三维可视化仿真规划系统建设方面的应用研究现状及所取得的主要成果，并结合遥感技术的发展方向，展望了其在煤炭地质中的应用前景。     

Comparison of magnetic parameters with vehicular Br levels in urban roadside soils

Chemical methods are generally chosen to monitor soil pollution but magnetic measurements proved to yield additional information at low cost and less time consumption. In this investigation, the novel use of rapid and non-destructive magnetic measurements to characterize Br levels in Xuzhou (China) urban roadside soils was reported. X-ray fluorescence spectrometry (XRF) was used to quantify Br in the soil samples. Data from 21 roadside soil samples confirm Br contamination, with a mean level of 4.36 mg kg⁻¹ and a range of 2.4–8.7 mg kg⁻¹. These values are higher than that of unpolluted soils in Xuzhou that averages 1.1 mg kg⁻¹. Hierarchical clustering analysis shows Br in Xuzhou roadside soils is mainly from road traffic. Clear correlations between Br levels and simple magnetic parameters [mass specific susceptibility (χ _LF), susceptibility of anhysteretic remanent magnetization (χ _ARM), saturation isothermal remanence (SIRM)] are observed. The present study shows that these three magnetic parameters can be used as a proxy for Br levels in Xuzhou urban roadside soils.     

青藏高原东缘龙门山晚新生代剥蚀厚度与弹性挠曲模拟

龙门山是青藏高原东缘边界山脉,具有青藏高原地貌、龙门山高山地貌和山前冲积平原三个一级地貌单元。利用数字高程模式图像和裂变径迹年代测定方法研究和计算龙门山晚新生代剥蚀厚度与剥蚀速率,结果表明:3.6 Ma以来龙门山的剥蚀厚度介于1.91-2.16 km之间,剥蚀速率介于0.53-0.60 mm/a之间。在此基础上,开展了该地区岩石圈的弹性挠曲模拟,结果表明龙门山的隆升机制具有以构造缩短隆升和剥蚀卸载隆升相叠合的特点。3.6 Ma之前,龙门山的隆升与逆冲推覆构造负载有关,以构造缩短驱动的构造隆升为特色;3.6 Ma之后,龙门山的隆升与剥蚀卸载驱动的抬升有关,并以剥蚀卸载隆升为特色,进而提出了龙门山晚新生代以来的隆升机制以剥蚀成山作用为主的认识。     

甘肃花岗岩类成矿作用研究与找矿方向

甘肃省花岗岩分布广泛，与其有关的矿产资源也十分丰富。笔者将甘肃省划分为北山、敦煌一阿拉善、祁连山及西秦岭4个构造-花岗岩-矿化区；进一步在北山、祁连山及西秦岭3个区划分出20个构造-花岗岩-矿化带；并就加强花岗岩成矿作用的研究和找矿工作提出了具体的工作思路与建议。