Tectonothermal events in the Olympic IOCG Province constrained by apatite and REE-phosphate geochronology

The Olympic iron oxide–copper–gold province in South Australia contains numerous deposits and prospects, including the Olympic Dam Cu–U–Au–Ag deposit and the Acropolis prospect. The Acropolis prospect comprises massive, coarse-grained magnetite–apatite veins partly replaced by a hematite-stable assemblage. The apatite grains in the veins contain zones with abundant inclusions of other minerals (including monazite and xenotime) and low trace-element concentrations relative to the inclusion-free zones. The inclusion-rich apatite zones are interpreted to be formed from the recrystallisation of the inclusion-free apatite and remobilisation of U, Th and rare earth element (REE) from apatite into monazite and xenotime. Apatite, monazite and xenotime are all established U–Th–Pb geochronometers and offer the potential to constrain the alteration history of the Acropolis prospect. The LA-ICPMS U–Pb age of inclusion-free apatite is within error of the age of the host volcanic units (ca 1.59 Ga). Inclusion-rich apatite yields both near-concordant analyses that are within error of the inclusion-free apatite as well as highly disturbed (discordant) analyses. The most concordant analyses of monazite (Th–Pb) inclusions and xenotime (U–Pb) inclusions and rim grains indicate an alteration event occurred at ca 1.37 Ga and possibly also at ca 500 Ma. The disparity in age of the inclusion-rich apatite and the REE-phosphate inclusions (and rim grains) is suggested to be owing to the apatite being initially recrystallised at ca 1.59 Ga and modified again by a later event that also formed (or coarsened) most of the inclusions. Partial resetting of the majority of the monazite inclusions as well as the presence of significant amounts of common Pb has complicated the interpretation of the monazite results. In contrast, xenotime is a more robust geochronometer in this setting. The ages of the two post-1.59 Ga events that appear to have affected the Acropolis prospect do not correspond to any events known to have occurred in the Gawler Craton. The earlier (ca 1.37 Ga) age instead corresponds best with metamorphic–magmatic–hydrothermal activity in Laurentia, consistent with the proximity of Laurentia and the Gawler Craton inferred from palaeogeographic reconstructions. The later (ca 500 Ma) event corresponds to the Delamerian Orogeny and has been shown by prior studies to have also affected the Olympic Dam deposit.     

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.     

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.     

模式内部变率引起的1.5℃和2℃升温阈值出现时间模拟的不确定性研究

模式内部变率是模拟结果不确定性的重要来源,然而它对于1.5℃和2℃升温阈值出现时间不确定性的影响尚不清楚。因此,基于耦合模式比较计划第五阶段(CMIP5)的多模式数据研究了模式内部变率对1.5℃和2℃升温阈值出现时间不确定性的影响以及对未来排放情景的敏感性。结果表明,模式内部变率对升温阈值出现时间模拟的影响与外强迫的影响相当,单个模式内部不同成员达到全球平均1.5℃或2℃增温的年份相差2~12年;其影响具有明显的空间差异,影响极大值出现在欧亚大陆以北洋面、白令海峡周围区域、北美东北部及其与格陵兰岛之间的海域、南半球高纬地区等;低排放情景下模式内部变率的影响大于高排放情景。     

基于数值预报模式的乡镇温度预报方法研究

利用2015年8月至2017年7月长兴岛站和交流岛站日最高气温、日最低气温实况资料,对ECMWF细网格2 m温度预报值和日本FSFE02(24 h地面形势场预报)、FSFE03(36 h地面形势场预报)进行了检验。结果表明:根据历史回归统计检验,ECMWF细网格模式24 h的2 m最高气温、最低气温预报效果显著,通过了0. 05信度显著性检验。对各月做相关分析,相关性均较好。利用前一日ECMWF细网格2 m温度预报值与长兴岛站实况差值,根据统计的ECMWF细网格2 m温度预报订正值,做出长兴岛站未来24 h的气温预报。交流岛站温度预报是在长兴岛站温度预报的基础上订正做出,经统计分析,交流岛站和长兴岛站的气温差值与地面形势场和风场有较好的对应关系,根据不同类型的地面形势场和风场订正值,做出交流岛站的温度预报。应用Matlab计算机语言的开发功能,提取ECMWF细网格2 m温度预报的最高、最低气温值,并录入当日长兴岛站和交流岛站最高、最低气温实况值,自动预报各站未来24 h最高气温、最低气温。创建可视化预报工作界面,实现乡镇温度预报自动化。     

基于WACCM+DART的临近空间SABER和MLS臭氧观测同化试验研究

本研究在WACCM+DART（Whole Atmosphere Community Climate Model,Data Assimilation Research Test-Bed）临近空间资料同化预报系统中加入SABER（Sounding of the Atmosphere using Broadband Emission Radiometry）和MLS（Microwave Limb Sounder）臭氧观测同化接口,并以2016年2月一次平流层爆发性增温（SSW）过程为模拟个例进行了SABER和MLS臭氧观测同化试验,得出以下结论:同化SABER和MLS臭氧体积浓度观测得出的WACCM+DART臭氧分析场能够较真实反映SSW期间北极上空平流层臭氧廓线随时间的演变特征,且与ERA5（Fifth Generation of ECMWF Reanalyses）再分析资料描述的臭氧变化特征具有很好的一致性;基于SABER和MLS臭氧观测的WACCM臭氧6 h预报检验表明同化臭氧观测对臭氧分析和预报误差的改善效果主要体现在南半球高纬平流层和北半球中高纬平流层中上层-中间层底部;基于ERA5再分析资料的WACCM+DART分析场检验表明同化SABER和MLS臭氧体积浓度资料可在提高北半球高纬地区上平流层-中间层底部臭氧场分析质量的同时减小该地区上平流层-中间层底部温度场和中间层底部纬向风场的分析误差;基于MLS臭氧资料的臭氧中期预报检验表明相对控制试验同化SABER和MLS臭氧体积浓度资料能更好改善0~5 d下平流层和中间层底部臭氧的预报效果。     

Latitudinal and Scan-dependent Biases of Microwave Humidity Sounder Measurements and Their Dependences on Cloud Ice Water Path

The relationship between differences in microwave humidity sounder(MHS)–channel biases which represent measured brightness temperatures and model-simulated brightness temperatures, and cloud ice water path(IWP) as well as the influence of the cloud liquid water path(LWP) on the relationship is examined. Seven years(2011–17) of NOAA-18 MHS-derived measured brightness temperatures and IWP/LWP data generated by the NOAA Comprehensive Large Array-data Stewardship System Microwave Surface and Precipitation Products System are used. The Community Radiative Transfer Model, version2.2.4, is used to simulate model-simulated brightness temperatures using European Center for Medium-Range Weather Forecasts reanalysis data as background fields. Scan-angle deviations of the MHS window channel biases range from-1.7 K to1.0 K. The relationships between channels 2, 4, and 5 biases and scan angle are symmetrical about the nadir. The latitudedependent deviations of MHS window channel biases are positive and range from 0–7 K. For MHS non-window channels,the latitudinal deviations between measured brightness temperatures and model-simulated brightness temperatures are larger when the detection height is higher. No systematic warm or cold deviations are found in the global spatial distribution of difference between measured brightness temperatures and model-simulated brightness temperatures over oceans after removing scan-angle and latitudinal deviations. The corrected biases of five different MHS channels decrease differently with respect to the increase in IWP. This decrease is stronger when LWP values are higher.