John Turnbull Thomson,Fellow of the Royal Geographical Society: An early New Zealand geographer

John Turnbull Thomson is more known in New Zealand for his work as Chief Surveyor of Otago Province (1856–1876) and first Surveyor General of New Zealand (1877–1879). He was also a Fellow of the Royal Geographical Society of London from 1848 until his death in 1884. Thomson was an important and early self‐declared ‘does’ of geographer in New Zealand, prior to the establishment of university geography in 1937. Thomson's contribution included survey and mapping but he also gave public lectures on geography and contributed to debates about glaciation.     

Depositional setting,provenance, and tectonic-volcanic setting of Eocene–Recent deep-sea sediments of the oceanic Izu–Bonin forearc,northwest Pacific (IODP Expedition 352)

New biostratigraphical, geochemical, and magnetic evidence is synthesized with IODP Expedition 352 shipboard results to understand the sedimentary and tectono-magmatic development of the Izu–Bonin outer forearc region. The oceanic basement of the Izu–Bonin forearc was created by supra-subduction zone seafloor spreading during early Eocene (c. 50–51 Ma). Seafloor spreading created an irregular seafloor topography on which talus locally accumulated. Oxide-rich sediments accumulated above the igneous basement by mixing of hydrothermal and pelagic sediment. Basaltic volcanism was followed by a hiatus of up to 15 million years as a result of topographic isolation or sediment bypassing. Variably tuffaceous deep-sea sediments were deposited during Oligocene to early Miocene and from mid-Miocene to Pleistocene. The sediments ponded into extensional fault-controlled basins, whereas condensed sediments accumulated on a local basement high. Oligocene nannofossil ooze accumulated together with felsic tuff that was mainly derived from the nearby Izu–Bonin arc. Accumulation of radiolarian-bearing mud, silty clay, and hydrogenous metal oxides beneath the carbonate compensation depth (CCD) characterized the early Miocene, followed by middle Miocene–Pleistocene increased carbonate preservation, deepened CCD and tephra input from both the oceanic Izu–Bonin arc and the continental margin Honshu arc. The Izu–Bonin forearc basement formed in a near-equatorial setting, with late Mesozoic arc remnants to the west. Subduction-initiation magmatism is likely to have taken place near a pre-existing continent–oceanic crust boundary. The Izu–Bonin arc migrated northward and clockwise to collide with Honshu by early Miocene, strongly influencing regional sedimentation.     

Middle Jurassic–Early Cretaceous tectonic evolution of the Bayanhushuo area,southern Great Xing’an Range,NE China: constraints from zircon U–Pb geochronological and geochemical data of volcanic and subvolcanic rocks

This study presents new zircon U–Pb geochronology, geochemistry, and zircon Hf isotopic data of volcanic and subvolcanic rocks that crop out in the Bayanhushuo area of the southern Great Xing’an Range (GXR) of NE China. These data provide insights into the tectonic evolution of this area during the late Mesozoic and constrain the evolution of the Mongol–Okhotsk Ocean. Combining these new ages with previously published data suggests that the late Mesozoic volcanism occurred in two distinct episodes: Early–Middle Jurassic (176–173 Ma) and Late Jurassic–Early Cretaceous (151–138 Ma). The Early–Middle Jurassic dacite porphyry belongs to high-K calc-alkaline series, showing the features of I-type igneous rock. This unit has zircon ε_Hf(t) values from +4.06 to +11.62 that yield two-stage model ages (T_DM2) from 959 to 481 Ma. The geochemistry of the dacite porphyry is indicative of formation in a volcanic arc tectonic setting, and it is derived from a primary magma generated by the partial melting of juvenile mafic crustal material. The Late Jurassic–Early Cretaceous volcanic rocks belong to high-K calc-alkaline or shoshonite series and have A₂-type affinities. These volcanics have ε_Hf(t) and T_DM2 values from +5.00 to +8.93 and from 879 to 627 Ma, respectively. The geochemistry of these Late Jurassic–Early Cretaceous volcanic rocks is indicative of formation in a post-collisional extensional environment, and they formed from primary magmas generated by the partial melting of juvenile mafic lower crust. The discovery of late Mesozoic volcanic and subvolcanic rocks within the southern GXR indicates that this region was in volcanic arc and extensional tectonic settings during the Early–Middle Jurassic and the Late Jurassic–Early Cretaceous, respectively. This indicates that the Mongol–Okhotsk oceanic plate was undergoing subduction during the Early–Middle Jurassic, and this ocean adjacent to the GXR may have closed by the Late Middle Jurassic–Early Late Jurassic.     

Early Carboniferous ophiolite in central Qiangtang,northern Tibet: record of an oceanic back-arc system in the Palaeo-Tethys Ocean

Zircon U–Pb dating of two samples of metagabbro from the Riwanchaka ophiolite yielded early Carboniferous ages of 354.4 ± 2.3 Ma and 356.7 ± 1.9 Ma. Their positive zircon εHf(t) values (+7.9 to +9.9) indicate that these rocks were derived from a relatively depleted mantle. The metagabbros can be considered as two types: R1 and R2. Both types are tholeiitic, with depletion of high-field-strength elements (HFSE) and enrichment of large-ion lithophile elements (LILE) similar to those of typical back-arc basin basalts (BABB), such as Mariana BABB and East Scotia Ridge BABB. Geochemical and isotopic characteristics indicate that the R1 metagabbro originated from a back-arc basin spreading ridge with addition of slab-derived fluids, whereas the R2 metagabbro was derived from a back-arc basin mantle source, with involvement of melts and fluids from subducted ocean crust. The Riwanchaka ophiolite exhibits both mid-ocean ridge basalts- and arc-like geochemical affinities, consistent with coeval ophiolites from central Qiangtang. Observations indicate that the Qiangtang ophiolites developed during the Late Devonian–early Carboniferous (D₃–C₁) in a back-arc spreading ridge above an intra-oceanic subduction zone. Based on our data and previous studies, we propose that an oceanic back-arc basin system existed in the Longmuco–Shuanghu–Lancang Palaeo-Tethys Ocean during the D₃–C₁ period.     

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.     

青岛连三岛地区“荆山岩群”高压变质岩的变质时代和变质演化特征

青岛连三岛地区原划为古元古界荆山岩群中出露各类片岩、片麻岩。本文在前人研究的基础上,对连三岛地区出露的含榴黑云母斜长片麻岩、含榴云母片岩和含榴黑云母钾长片麻岩,进行详细的岩相学和矿物化学研究,确定其变质温压条件及其P-T演化轨迹,并采用LA-ICP-MS锆石U-Pb定年,获得了3个样品的原岩时代和变质时代,为全面深入认识其变质属性提供了进一步的重要依据。根据岩相学和矿物化学成分可以识别出两期矿物组合:第一期(峰期变质阶段)为Grt1+Kfs+Aln+Ph+Qtz;第二期(退变质阶段)为Grt2+Pl+Ep+Bt+Qtz;依据多硅白云母Si压力计、锆石Ti温度计以及GB-GBPQ矿物温压计,确定其峰期变质和退变质的温压条件分别为T=600～817℃、P=2.4～2.6GPa和T=431～456℃、P=0.48～0.82GPa。结合白云母部分熔融现象,上述两个变质阶段构成了一个折返早期升温降压,后穿过多硅白云母熔融反应线,最后降温降压的顺时针型演化的P-T轨迹。CL图像显示3个样品的锆石均具有典型的岩浆核-变质边结构;结合LA-ICP-MS锆石原位微区U-Pb定年和微量元素分析,3个样品分别获得了769～756Ma、～223Ma和213～216Ma三组年龄,分别与苏鲁造山带其他高压-超高压变质岩的新元古代原岩时代(700～800Ma)、峰期变质时代(240～225Ma)和角闪岩相退变质时代(215～205Ma)一致。对样品含榴黑云母斜长片麻岩(17LSD-1)和含榴云母片岩(18LSD-2)进行锆石Hf同位素分析,获得样品17LSD-1的ε_(Hf)(t)=-23.2～2.8、t_(DM2)~C(Hf)=1712～2845Ma,表明其原岩主要来源于古元古代陆壳重熔;样品18LSD-2的ε_(Hf)(t)=-13.9～8.6、t_(DM2)~C(Hf)=1113～2358Ma,表明样品形成时壳源物质成分占主导地位,同时部分幔源或新生地壳物质的加入,导致少部分ε_(Hf)(t)偏正值。Hf同位素结果表明连三岛变质岩原岩的形成与扬子板块新太古代-早古元古代的陆壳重熔有关。对比前人的相关数据,无论是原岩时代、变质年龄还是变质演化特征,本文研究的连三岛地区片岩/片麻岩与苏鲁造山带的部分变质岩均具有相似的原岩属性和变质属性,因此推断其应归属为苏鲁超高压变质带的一部分,是三叠纪扬子板块向华北板块俯冲碰撞引发的高压-超高压变质事件的产物,不应再作为岩石地层单元划归为"荆山岩群"。