Predictive geochemical mapping using environmental correlation

The distribution of chemical elements at and near the Earth's surface, the so-called critical zone, is complex and reflects the geochemistry and mineralogy of the original substrate modified by environmental factors that include physical, chemical and biological processes over time.Geochemical data typically is illustrated in the form of plan view maps or vertical cross-sections, where the composition of regolith, soil, bedrock or any other material is represented. These are primarily point observations that frequently are interpolated to produce rasters of element distributions. Here we propose the application of environmental or covariate regression modelling to predict and better understand the controls on major and trace element geochemistry within the regolith. Available environmental covariate datasets (raster or vector) representing factors influencing regolith or soil composition are intersected with the geochemical point data in a spatial statistical correlation model to develop a system of multiple linear correlations. The spatial resolution of the environmental covariates, which typically is much finer (e.g. ∼90 m pixel) than that of geochemical surveys (e.g. 1 sample per 10-10,000 km²), carries over to the predictions. Therefore the derived predictive models of element concentrations take the form of continuous geochemical landscape representations that are potentially much more informative than geostatistical interpolations.Environmental correlation is applied to the Sir Samuel 1:250,000 scale map sheet in Western Australia to produce distribution models of individual elements describing the geochemical composition of the regolith and exposed bedrock. As an example we model the distribution of two elements – chromium and sodium. We show that the environmental correlation approach generates high resolution predictive maps that are statistically more accurate and effective than ordinary kriging and inverse distance weighting interpolation methods. Furthermore, insights can be gained into the landscape processes controlling element concentration, distribution and mobility from analysis of the covariates used in the model. This modelling approach can be extended to groups of elements (indices), element ratios, isotopes or mineralogy over a range of scales and in a variety of environments.     

Laboratory Simulation of the Formation Process of Surface Geochemical Anomalies Applied to Hydrocarbon Exploration

The formation mechanisms and processes of geochemical anomalies used as proxies in surface geochemistry exploration (SGE) have not been well understood. Previous studies cannot realize 3D measurement of microseeping hydrocarbons from reservoirs to the surface, which made it difficult to understand the features and pathways of deep hydrocarbon microseepages. Understanding the processes of hydrocarbon microseepages will contribute to the acceptance and effectiveness of surface geochemistry. Based on a simplified geological model of hydrocarbon microseepages, including hydrocarbon reservoir, direct caprock, overlying strata and Quaternary sediments, this work established a 3D experimental system to simulate the mechanisms and processes of deep hydrocarbon microseepes extending to the surface. The dispersive halos of microseeping hydrocarbons in the subsurface were adequately described by using this 3D experimental system. Results indicate that different migration patterns of hydrocarbons above the point gas source within the simulated caprock and overlying strata can be reflected by the ratio of i-butane to n-butane (i-C4/n-C4), which follow diffusion and infiltration (buoyancy) mechanisms. This is not the case for vertical measurement lines far from the point gas source. A vertical gas flow in the form of a plume was found during hydrocarbon microseepage. For sampling methods, the high-density grid sampling is favorable for delineating prospecting targets. Hydrocarbon infiltration or buoyancy flow occurs in the zones of infiltration clusters, coupling with a diffusion mechanism at the top of the water table and forming surface geochemical anomalies. These results are significant in understanding hydrocarbon microseepage and interpreting SGE data.     

华南大规模低温成矿的主要科学问题

大面积低温成矿是全球独特的重要成矿事件。华南地区扬子地块西南部面积约50万平方千米的广大范围,低温矿床广泛发育,包括卡林型金矿床、MVT型Pb-Zn矿床和脉型锑、汞、砷等矿床,构成华南低温成矿域。该成矿域由川滇黔接壤区的Pb-Zn、右江盆地Au-Sb-As-Hg和湘中盆地Sb-Au等三个矿集区组成。其中的矿床主要赋存于沉积岩中,受断裂构造控制,形成于100~250℃的低温条件下。以往的研究取得了重要进展,但成矿时代、成矿动力学背景、成矿物质基础和成矿过程等关键科学问题一直悬而未决。近年的研究表明,这些矿床可能形成于200~230Ma(印支期)和130~160Ma(燕山期)两个时期,除川滇黔Pb-Zn矿集区由盆地卤水演化成矿外,其余两个矿集区的形成可能均受印支期和燕山期的深部岩浆活动驱动,尽管成矿流体中的水主要是大气成因地下水。本专题报道了近年来这些方面的最新研究进展,包括20篇文章,主要涉及这些低温矿床的地质地球化学特征、成矿时代、矿床成因和成矿动力学背景等。     

青岛冷水团的消亡机理研究

本文基于多年月平均水温资料,分析了青岛冷水团的长消过程,并利用气候态月平均大气数据和数值模拟结果,探讨了青岛冷水团的消亡机理。结果表明,青岛冷水团3月出现,4月成型,5月最盛,6月减弱,7月消失;南黄海6-7月间偏南风的增强和温跃层以下反气旋涡的减弱是青岛冷水团消亡的动力机制,而海面净热通量的下传和水平热量的输入则是青岛冷水团消亡的热力机制。     

Foundations of petroleum prediction based on the activation of geochemical processes by direct current

Petroleum prediction in local objects and the applicability of geologic environment parameters in this regard are considered. The proposed guidelines for petroleum potential prediction are based on the analysis of the distribution of typomorphic epigenetic multivalent elements loosely fixed in rocks before and after electric current treatment. The advantage of geoelectrochemical methods providing data on the chemical composition of above-productive deposits is demonstrated.     

Origin of natural sulfur-metal chimney in the Tangyin hydrothermal field,Okinawa Trough: constraints from rare earth element and sulfur isotopic compositions

For the first time, we present the rare earth element (REE) and sulfur isotopic composition of hydrothermal precipitates recovered from the Tangyin hydrothermal field (THF), Okinawa Trough at a water depth of 1206 m. The natural sulfur samples exhibit the lowest ΣREE concentrations (ΣREE= 0.65×10^–6–4.580×10^–6) followed by metal sulfides (ΣREE=1.71×10^–6–11.63×10^–6). By contrast, the natural sulfur-sediment samples have maximum ΣREE concentrations (ΣREE=11.54×10^–6–33.06×10^–6), significantly lower than those of the volcanic and sediment samples. Nevertheless, the δEu, δCe, (La/Yb)_N, La/Sm, (Gd/Yb)_N and normalized patterns of the natural sulfur and metal sulfide show the most similarity to the sediment. Most hydrothermal precipitate samples are characterized by enrichments of LREE (LREE/HREE=10.09–24.53) and slightly negative Eu anomalies or no anomaly (δEu=0.48–0.99), which are different from the hydrothermal fluid from sediment-free mid-oceanic ridges and back-arc basins, but identical to the sulfides from the Jade hydrothermal field. The lower temperature and more oxidizing conditions produced by the mixing between seawater and hydrothermal fluids further attenuate the leaching ability of hydrothermal fluid, inducing lower REE concentrations for natural sulfur compared with metal sulfide; meanwhile, the negative Eu anomaly is also weakened or almost absent. The sulfur isotopic compositions of the natural sulfur (δ³⁴S=3.20‰–5.01‰, mean 4.23‰) and metal sulfide samples (δ³⁴S=0.82‰–0.89‰, mean 0.85‰) reveal that the sulfur of the chimney is sourced from magmatic degassing.     

黄河口清水沟河道的冲淤过程与模拟

为揭示黄河口清水沟河道长时段的冲淤演变规律并建立其冲淤计算方法,分析了清水沟1976—2015年的时空冲淤演变过程,采用河床演变的滞后响应模型,考虑河口来水来沙及河道延伸与蚀退的影响,建立了清水沟累计冲淤量的计算方法。结果表明:1976—1980年改道初期清水沟改道点上游先冲后淤,改道点下游淤滩塑槽,淤积量随着下游河道展宽而增加,1980年后改道点上、下游河道冲淤过程趋于一致;受水沙条件等因素影响,1980—1986年清水沟主槽冲刷展宽,之后主槽淤积萎缩;1996年清八改汊和2002年小浪底水库"调水调沙"原型试验以来,河道转淤为冲,2002年后河道冲刷速率随时间指数衰减;河床演变的滞后响应模型可计算清水沟长时段的冲淤过程,该方法可为预测未来清水沟冲淤演变趋势提供科学参考。     

四川盆地云和降水观测科学试验设计与实践

2014—2017年四川地区开展了大范围云系观测科学试验。观测对象以盆地层状云系和积层混合云系为主,积状云(对流云)为辅。本文围绕试验目标、区域、观测要素、观测布局设计、观测方案设计、设备技术参数和典型个例等7个方面进行介绍,根据不同类型云系和降水变化特征,设计有针对性的观测方案,获得了不同类型云系和降水的多尺度连续性观测数据,为四川地区开展云和降水关系研究提供详实的综合型外场观测数据。层状降水云典型个例云高可达8km,云强核心部位的云雷达反射率可达28dBz,径向速度可达-6m·s~(-1),在0℃附近,反射率和退偏振因子LDR上有一条明显的亮带,表现为极大值,液态水主要集中于4.5km以下,随着雨强增大,液水含量增加,降水滴谱分布较窄,随着雨强减小,雨滴谱和速度谱变窄,但小粒径数浓度增加,说明对层状云雨强起主导作用的是雨滴直径,而不是数浓度;无降水层状云典型个例云层厚度为3.2km,云顶约为4km,云雷达反射率不超过0dBz,径向速度不超过5m·s~(-1),层状云内整层水汽含量和液水含量较为稳定,云中主要为液态粒子且粒径偏小、小粒径数浓度较高。     

一次拖曳型飑线过程雨滴谱演变特征研究

利用Thies激光雨滴谱仪观测资料和CINRAD/SA多普勒雷达观测资料,分析了2017年7月18日一次典型中纬度拖曳型飑线过程不同发展阶段雨滴谱和积分参数的演变特征,主要结果为：1）成熟飑线回波包括对流带、过渡区和拖曳层状云区三部分,对流带前侧不断有对流带生成并合并到主对流带中,使得对流带的前沿具有强的反射率因子,并且有多个雨强大值中心。2）垂直穿过飑线对流带,雨强增加阶段有较少的小粒子（直径小于1 mm）和特大粒子（直径大于5 mm）,以及较低的雨滴浓度和反射率因子,而雨强减弱阶段有较多的小粒子和特大粒子,以及较高的雨滴浓度和反射率因子;飑线加强阶段,雨滴谱有较大的峰值直径（0.44 mm）、较多的大（直径大于3 mm）和特大粒子,而飑线减弱阶段,雨滴谱有较小的峰值直径（0.19 mm）、较少的大粒子。3）对流带、过渡区和层状云降水雨滴谱的Gamma谱三参数N0、μ、λ随雨强增大有明显的分层特征,相同雨强时,对流云和过渡区降水的三参数比层状云降水的数值大;而飑线不同发展阶段、不同降水类型的λ-μ关系具有一致性,二次多项式可以很好地拟合λ-μ关系。4）归一化雨滴谱参数NW和D0的分布可以用来区分对流云和层状云降水,并给出了新的分离线方程;另外,飑线在发展和减弱阶段的雨滴谱特征有明显差异,表明飑线演变过程降水形成的微物理机制发生变化,前期冷云过程有重要影响,而后期暖云过程起主导作用。     

线性化物理过程对GRAPES 4DVAR同化的影响

线性化物理过程能够改善四维变分同化中极小化收敛的稳定性和增加极小化过程中对大气物理过程和动力更加精确的描述,它是四维变分同化中非常重要的一部分。通过在GRAPES全球模式中研究线性化物理过程,尤其是两个湿线性化物理过程,改善切线性模式预报精度,来提高GRAPES全球四维变分同化的分析和预报效果。线性化物理过程的开发首先需要简化原非线性化物理过程中的强非线性项,然后对线性化物理过程进行规约化,以抑制切线性扰动的异常增长。目前GRAEPS全球模式中的线性化物理过程主要包括次网格尺度地形参数化、垂直扩散、积云深对流和大尺度凝结。线性化物理过程预报精度的检验方法是通过选择合适大小的初始扰动（同化分析增量）,来比较非线性模式和切线性模式中的扰动演化的纬向平均误差。然后以绝热版本的切线性模式为基础,通过冬、夏两个个例试验来分别检验4个线性化物理过程的12 h预报效果。试验结果表明,通过添加次网格地形参数化和垂直扩散两个干线性化物理过程方案,可以有效抑制住绝热版本切线性模式低层扰动的异常增长,大幅度改善切线性模式预报效果。通过添加积云深对流和大尺度凝结两个湿线性化物理过程,可以在热带区域和中、高纬度地区提高切线性模式中湿变量和温度变量的近似精度,提高GRAPES全球四维变分同化的分析和预报效果。