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.     

渔洞水库地区地震活动性和构造应力场分析

统计了渔洞水库的水位变化规律,分析了水库蓄水前后库区及其附近地区的地震活动性和构造应力场特征。结果表明,水库蓄水前后研究区地震活动的图像、强度和频次都没有明显变化;水库蓄水和水位变化对研究区构造应力场没有明显的影响。     

印尼8.7、8.5级远震对云南地区小震的激发和地下流体异常响应初探

2004年12月26日和2005年3月29日,印度尼西亚苏门答腊西北海域相继发生8.7和8.5级地震。通过对这两次地震与云南地震活动关系和地下流体记录异常变化的分析,认为:印尼两次8级地震的发生可能会使云南提前进入新一轮强震活跃期;远震对云南地区小震的激发和导致地下流体记录出现异常变化的因素虽然都与大地震面波的强弱、地质构造环境等有关,但机理不同。     

Hydrodynamic performance of a T-shaped floating breakwater

The comprehensive utilization of floating breakwaters, specially acting as a supporting structure for offshore marine renewable energy explorations, has received more and more attention recently. Based on linear water-wave theory, the hydrodynamic performance of a T-shaped floating breakwater is semi-analytically investigated through the matched eigenfunction expansion method (MEEM). Auxiliary functions, to speed up the convergence and improve the accuracy in the numerical computations, are introduced to represent the singular behavior of fluid field near the lower salient corners of the structure. The effects of the height and installation position of the vertical screen on the reflection and transmission coefficients, dynamic response and wave forces are examined. It is found that the presence of the screen shifts the resonance frequency of RAO for both surge and pitch modes to the low-frequency area, while has no effect on heave mode. The identical added masses, damping and transmission coefficients can be obtained in the cases where the screen holds the same distance away from the longitudinal central axis of the upper box-type structure. Moreover, a relatively small pitch response can be achieved in a wide wave–frequency range, when the breakwater is Γ-shaped.     

An efficient and robust approach to predict ship self-propulsion coefficients

Achieving a reliable and accurate numerical prediction of the self-propulsion performance of a ship is still an open problem that poses some relevant issues. Several CFD methods, ranging from boundary element methods (BEM) to higher-fidelity viscous Reynolds averaged Navier–Stokes (RANS) based solvers, can be used to accurately analyze the separate problems, i.e. the open water propeller and the hull calm water resistance. However, when the fully-coupled self-propulsion problem is considered, i.e. the hull advancing at uniform speed propelled by its own propulsion system, several complexities rise up. Typical flow simplifications adopted to speed-up the simulations of the single analysis (hull and propeller separately) lose their validity requiring a more complex solver to tackle the fully-coupled problem. The complexity rises up further when considering a maneuver condition. This aspect increases the computational burden and, consequently, the required time which becomes prohibitive in a preliminary ship design stage.The majority of the simplified methods proposed in literature to include propeller effects, without directly solve the propeller flow, in a high-fidelity viscous solver are not able to provide all the commonly required self-propulsion coefficients. In this work, a new method to enrich the results from a body force based approach is proposed and investigated, with the aim to reduce as much as possible the computational burden without losing any useful result. This procedure is tested for validation on the KCS hull form in self-propulsion and maneuver conditions.     

CFD modelling of a small–scale fixed multi–chamber OWC device

Wave Energy Converters (WECs) have excellent potential as a source of renewable energy that is yet to be commercially realised. Recent attention has focused on the installation of Oscillating Water Column (OWC) devices as a part of harbor walls to provide advantages of cost–sharing structures and proximity of power generation facilities to existing infrastructure. In this paper, an incompressible three–dimensional CFD model is constructed to simulate a fixed Multi–Chamber OWC (MC–OWC) device. The CFD model is validated; the simulation results are found to be in good agreement with experimental results obtained from a scale physical model tested in a wave tank. The validated CFD model is then used for a benchmark study of 96 numerical tests. These investigate the effects of the PTO damping caused by the power take–off (PTO) system on device performance. The performance is assessed for a range of regular wave heights and periods. The results demonstrate that a PTO system with an intermediate damping can be used for all chambers in the MC–OWC device for most wave period ranges, except for the long wave periods. These require a higher PTO damping. An increased incident wave height reduces the device capture width ratio, but there is a noticeable improvement for long wave periods.