Application of Geochemistry and VNIR Spectroscopy in Mapping Heavy Metal Pollution of Stream Sediments in the Takab Mining Area, NW of Iran

This study considered the possibility of using visible and near infrared (VNIR) spectral absorption feature parameters (SAFPs) in predicting the concentration and mapping the distribution of heavy metals in sediments of the Takab area. In total, 60 sediment samples were collected along main streams draining from the mining districts and tailing sites, in order to measure the concentration of As, Co, V, Cu, Cr, Ni, Hg, Ti, Pb and Zn and the reflectance spectra (350–2500 nm). The quantitative relationship between SAFPs (Depth500nm, R610/500nm, R1344/778nm, Area500nm, Depth2200nm, Area2200nm, Asym2200nm) and geochemical data were assessed using stepwise multiple linear regression (SMLR) and enter multiple linear regression (EMLR) methods. The results showed a strong negative correlation between Ni and Cr with Area2200nm, a significant positive correlation between As and Asym2200nm, Ni and Co with Depth2200nm, as well as Co, V and total values with Depth500nm. The EMLR method eventuated in a significant prediction result for Ni, Cr, Co and As concentrations based on spectral parameters, whereas the prediction for Zn, V and total value was relatively weak. The spatial distribution pattern of geochemical data showed that mining activities, along with the natural weathering of base metal occurrences and rock units, has caused high concentrations of heavy metals in sediments of the Sarough River tributaries.     

Development of groundwater favourability map using GIS-based driven data mining models: an approach for effective groundwater resource management

The development of groundwater favourability map is an effective tool for the sustainability management of groundwater resources in typical agricultural regions, such as southern Perak Province, Malaysia. Assessing the potentiality and pollution vulnerability of groundwater is a fundamental phase of favourability mapping. A geographic information system (GIS)-based Boolean operator of a spatial analyst module was applied to combine a groundwater potentiality map (GPM) model and a groundwater vulnerability to pollution index (GVPI) map, thereby establishing the favourable zones for drinking water exploration in the investigated area. The area GPM model was evaluated by applying a GIS-based Dempster–Shafer–evidential belief function model. In the evaluation, six geoelectrically determined groundwater potential conditioning factors (i.e. overburden resistivity, overburden thickness, aquifer resistivity, aquifer thickness, aquifer transmissivity and hydraulic conductivity) were synthesized by employing the probability-based algorithms of the model. The generated thematic maps of the seven hydrogeological parameters of the DRASTIC model were considered as pollution potential conditioning factors and were analysed with the developed ordered weighted average–DRASTIC index model algorithms to construct the GVPI map. Approximately 88.8 and 85.71% prediction accuracies for the Groundwater Potentiality and GVPI maps were established using the reacting operating characteristic curve method and water quality status–vulnerability zone relationship scheme, respectively. Finally, the area groundwater favourability map (GFM) model was produced by applying a GIS-based Boolean operator on the Groundwater Potentiality and GVPI maps. The GFM model reveals three distinct zones: ‘not suitable’, ‘less suitable’ and ‘very suitable’ zones. The area analysis of the GFM model indicates that more than 50% of the study area is covered by the ‘very suitable’ zones. Results produce a suitability map that can be used by local authorities for the exploitation and management of drinking water in the area. The study findings can also be applied as a tool to help increase public awareness of groundwater issues in developing countries.     

Contrasting water–rock interaction behaviors of antimony and arsenic in contaminated rivers around an antimony mine,Hunan Province,China

Although antimony (Sb) and arsenic (As) exhibit similar geochemical behavior and toxicity in the environment, growing evidence suggests that their water–rock interaction behavior in contaminated rivers is quite different. Twenty-nine river water samples were collected between September and November 2018 from contaminated rivers around an antimony mine in Hunan Province, China. The concentrations of As and Sb were inversely proportional to the water flow distance. The rates and magnitudes of Sb decrease were more prominent than those of As. Silicate mineral dissolution from rocks such as silicified limestone increased the As and Sb concentration of in-mine-district (IMD) water. Dissolution of carbonate minerals, ion exchange, and competitive adsorption were the major water–rock interactions, resulting in rapidly decreasing As and Sb concentration in IMD direct impacted water and IMD indirect impacted water. The behaviors of As and Sb during water–rock interaction were dissimilar for areas dominated by carbonate and silicate minerals.     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

矿石体重回归分析在新疆乌宗布拉克矿区的应用研究

在资源储量估算的过程中,矿石小体重的准确与否将直接影响到资源储量估算的客观程度。对矿石小体重的准确预测将需要通过大量的统计分析工作,传统的方法是利用矿石小体重的算术平均值进行估算矿床资源储量,并没有考虑到多矿种对矿石小体重的影响程度。本文利用Excel软件"数据分析"的"回归"功能模块,对实验室所测定的矿石小体重值与其对应品位进行二元线性回归分析,快速、准确地构建矿石体重与其品位的数学模型,从而为资源储量估算提供了更客观、更科学的矿石体重模型。     

Integration of zircon and apatite U–Pb geochronology and geochemical mapping of the Wude basalts(Emeishan large igneous province):A tool for a better understanding of the tectonothermal and geodynamic evolution of the Emeishan LIP

Radiogenic isotopic dating and Lu–Hf isotopic composition using laser ablation-inductively coupled plasma-mass spectrometry(LA-ICP-MS)of the Wude basalt in Yunnan province from the Emeishan large igneous province(ELIP)yielded timing of formation and post-eruption tectonothermal event.Holistic lithogeochemistry and elements mapping of basaltic rocks were further reevaluated to provide insights into crustal contamination and formation of the ELIP.A zircon U–Pb age of 251.3±2.0 Ma of the Wude basalt recorded the youngest volcanic eruption event and was consistent with the age span of 251-263 Ma for the emplacement of the ELIP.Such zircons hadε_Hf(t)values ranging from7.3 to+2.2,identical to those of magmatic zircons from the intrusive rocks of the ELIP,suggesting that crust-mantle interaction occurred during magmatic emplacement,or crust-mantle mixing existed in the deep source region prior to deep melting.The apatite U–Pb age at 53.6±3.4 Ma recorded an early Eocene magmatic superimposition of a regional tectonothermal event,corresponding to the Indian–Eurasian plate collision.Negative Nb,Ta,Ti and P anomalies of the Emeishan basalt may reflect crustal contamination.The uneven Nb/La and Th/Ta values distribution throughout the ELIP supported a mantle plume model origin.Therefore,the ELIP was formed as a result of a mantle plume which was later superimposed by a regional tectonothermal event attributed to the Indian–Eurasian plate collision during early Eocene.     

Geochronological and geochemical evidence for a Late Ordovician to Silurian arc-back-arc system in the northern Great Xing'an Range,NE China

The early Paleozoic tectonic evolution of the Xing'an-Mongolian Orogenic Belt is dominated by two oceanic basins on the northwestern and southeastern sides of the Xing'an Block,i.e.,the Xinlin-Xiguitu Ocean and the Nenjiang Ocean.However,the early development of the Nenjiang Ocean remains unclear.Here,we present zircon U-Pb geochronology and whole-rock elemental and Sr-Nd isotopic data on the gabbros in the Xinglong area together with andesitic tuffs and basalts in the Duobaoshan area.LA-ICP-MS zircon U-Pb dating of gabbros and andesitic tuffs yielded crystallization ages of 443-436 Ma and 452-451 Ma,respectively.The Early Silurian Xinglong gabbros show calc-alkaline and E-MORB affinities but they are enriched in LILEs,and depleted in HFSEs,with relatively low U/Th ratios of 0.18-0.36 andεNd(t)values of-1.6 to+0.5.These geochemical features suggest that the gabbros might originate from a mantle wedge modified by pelagic sediment-derived melts,consistent with a back-arc basin setting.By contrast,the andesitic tuffs are characterized by high MgO(>5 wt.%),Cr(138-200 ppm),and Ni(65-110 ppm)contents,and can be termed as high-Mg andesites.Their low Sr/Y ratios of 15.98-17.15 and U/Th values of 0.24-0.25 and moderate(La/Sm)_n values of 3.07-3.26 are similar to those from the Setouchi Volcanic Belt(SW Japan),and are thought to be derived from partial melting of subducted sediments,and subsequent melt-mantle interaction.The Duobaoshan basalts have high Nb(8.44-10.30 ppm)and TiO2 contents(1.17-1.60 wt.%),typical of Nb-enriched basalts.They are slightly younger than regional adakitic rocks and have positiveεNd(t)values of+5.2 to+5.7 and are interpreted to be generated by partial melting of a depleted mantle source metasomatized by earlier adakitic melts.Synthesized with coeval arc-related igneous rocks from the southeastern Xing'an Block,we propose that the Duobaoshan high-Mg andesitic tuffs and Nbenriched basalts are parts of the Late Ordovician and Silurian Sonid Zuoqi-Duobaoshan arc belt,and they were formed by the northwestern subduction of the Nenjiang Ocean.Such a subduction beneath the integrated Xing'an-Erguna Block also gave rise to the East Ujimqin-Xinglong igneous belt in a continental back-arc basin setting.Our new data support an early Paleozoic arc-back-arc model in the northern Great Xing'an Range.     

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

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

GRAPES全球奇异向量方法改进及试验分析

基于总能量模的奇异向量扰动常用于构造集合预报的初始条件。以建立GRAPES（Global and Regional Assimilation PrEdiction System）全球集合预报系统为目的,基于前期研发的GRAPES全球模式奇异向量方法,在GRAPES全球切线性模式和伴随模式2.0版的框架下,开展了引入线性化边界层方案来改善奇异向量结构,并提高奇异向量计算效率的研究。通过连续试验,从奇异向量的扰动能量结构、扰动能量谱及扰动空间分布等方面,综合分析改进GRAPES全球奇异向量的结构及演变特征。试验结果表明,改进后的GRAPES奇异向量方法有效抑制了之前扰动能量在近地面层不合理的快速增长,同时,奇异向量最优扰动的结构更客观地体现了中高纬度区域大气初始条件中的斜压不稳定扰动及其演变,如在初始时刻奇异向量扰动能量主要位于对流层中层,并呈现出随高度向西倾斜的大气斜压特征;经过线性化演变,扰动能量向较大水平尺度转移,并在垂直结构上表现出向对流层高层上传及向对流层低层下传的特征等。针对GRAPES奇异向量迭代求解中伴随模式计算耗时为主的情况,改进伴随模式中广义共轭余差方案的调用方式,并采用大内存存储法来提高其计算效率,进而将奇异向量总计算时间缩短了25%。总之,改进后的GRAPES奇异向量方法,可应用于构建面向业务应用的GRAPES全球集合预报系统。      

基于镜像图的LRC和CRC偏差结合的人脸识别

为了提高人脸识别率及更好地显示人脸特征,本文提出了一种基于镜像图的LRC和CRC偏差结合的人脸识别方法.该方法首先生成一种镜像人脸,再通过融合原始人脸和镜像人脸形成新的混合训练样本,最后利用LRC和CRC偏差结合进行人脸识别.新方法增加了训练样本的数目,克服了由于光照和姿态等外部因素带来的影响.实验结果表明,镜像图与LRC和CRC偏差结合的人脸识别方法提高了人脸识别的准确性.