ALBERT

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Infrared spectroscopy in the visible to near‐infrared (vis–NIR) and mid‐infrared (MIR) regions is a well‐established approach for the prediction of soil properties. Different data fusion and training approaches exist, and the optimal procedures are yet undefined and may depend on the heterogeneity present in the set and on the considered scale. The objectives were to test the usefulness of partial least squares regressions (PLSRs) for soil organic carbon (SOC), total carbon (C〈sub〉t〈/sub〉), total nitrogen (N〈sub〉t〈/sub〉) and pH using vis–NIR and MIR spectroscopy for an independent validation after standard calibration (use of a general PLSR model) or using memory‐based learning (MBL) with and without spiking for a national spectral database. Data fusion approaches were simple concatenation of spectra, outer product analysis (OPA) and model averaging. In total, 481 soils from an Austrian forest soil archive were measured in the vis–NIR and MIR regions, and regressions were calculated. Fivefold calibration‐validation approaches were carried out with a region‐related split of spectra to implement independent validations with n ranging from 47 to 99 soils in different folds. MIR predictions were generally superior over vis–NIR predictions. For all properties, optimal predictions were obtained with data fusion, with OPA and spectra concatenation outperforming model averaging. The greatest robustness of performance was found for OPA and MBL with spiking with 〈italic toggle="no"〉R〈/italic〉〈sup〉2〈/sup〉 ≥ 0.77 (N), 0.85 (SOC), 0.86 (pH) and 0.88 (C〈sub〉t〈/sub〉) in the validations of all folds. Overall, the results indicate that the combination of OPA for vis–NIR and MIR spectra with MBL and spiking has a high potential to accurately estimate properties when using large‐scale soil spectral libraries as reference data. However, the reduction of cost‐effectiveness using two spectrometers needs to be weighed against the potential increase in accuracy compared to a single MIR spectroscopy approach.〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉In recent years, many two‐dimensional (2D) hydrodynamic models have been extended to include the direct rainfall method (DRM). This allows their application as a hydrological‐hydrodynamic model for the determination of floodplains in one model system. In previous studies on DRM, the role of catchment hydrological processes (CaHyPro) and its interaction with the calibration process was not investigated in detail. In the present, case‐oriented study, the influence of the spatiotemporal distribution of the processes precipitation and runoff formation in combination with the 2D model HEC‐RAS is investigated. In a further step, a conceptual approach for event‐based interflow is integrated. The study is performed on the basis of a single storm event in a small rural catchment (low mountain range, 38 km〈sup〉2〈/sup〉) in Hesse (Germany). The model results are evaluated against six quality criteria and compared to a simplified baseline model. Finally, the calibrated improved model is contrasted with a calibrated baseline model. The results show the enhancement of the model results due to the integration of the CaHyPro and highlight its interplay with the calibrated model parameters.〈/p〉

Description: Spatiotemporal characterisation of the soil redox status within the capillary fringe (CF) is a challenging task. Air‐filled porosities (ε), oxygen concentration (O〈sub〉2〈/sub〉) and soil redox potential (EH) are interrelated soil variables within active biogeochemical domains such as the CF. We investigated the impact of water table (WT) rise and drainage in an undisturbed topsoil and subsoil sample taken from a Calcaric Gleysol for a period of 46 days. We merged 1D (EH and matric potential) and 2D (O〈sub〉2〈/sub〉) systems to monitor at high spatiotemporal resolution redox dynamics within self‐constructed redoxtron housings and complemented the data set by a 3D pore network characterization using X‐ray microtomography (X‐ray μCT). Depletion of O〈sub〉2〈/sub〉 was faster in the organic matter‐ and clay‐rich aggregated topsoil and the CF extended 〉10 cm above the artificial WT. The homogeneous and less‐aggregated subsoil extended only 4 cm above the WT as indicated by ε–O〈sub〉2〈/sub〉–EH data during saturation. After drainage, 2D O〈sub〉2〈/sub〉 imaging revealed a fast aeration towards the lower depths of the topsoil, which agrees with the connected ε derived by X‐ray μCT (ε〈sub〉CT_conn〈/sub〉) of 14.9% of the total porosity. However, small‐scaled anoxic domains with O〈sub〉2〈/sub〉 saturation 〈5% were apparent even after lowering the WT (down to 0.25 cm〈sup〉2〈/sup〉 in size) for 23 days. These domains remained a nucleus for reducing soil conditions (E〈sub〉H〈/sub〉 〈 −100 mV), which made it challenging to characterise the soil redox status in the CF. In contrast, the subsoil aeration reached O〈sub〉2〈/sub〉 saturation after 8 days for the complete soil volume. Values of ε〈sub〉CT_conn〈/sub〉 around zero in the subsoil highlighted that soil aeration was independent of this parameter suggesting that other variables such as microbial activity must be considered when predicting the soil redox status from ε alone. The use of redoxtrons in combination with localised redox‐measurements and image based pore space analysis resulted in a better 2D/3D characterisation of the pore system and related O〈sub〉2〈/sub〉 transport properties. This allowed us to analyse the distribution and activity of microbiological niches highly associated with the spatiotemporal variable redox dynamics in soil environments. Highlights: The time needed to turn from reducing to oxidising (period where all platinum electrodes feature E〈sub〉H〈/sub〉 〉 300 mV) condition differ for two samples with contrasting soil structure. The subsoil with presumably low O〈sub〉2〈/sub〉 consumption rates aerated considerably faster than the topsoil and exclusively by O〈sub〉2〈/sub〉 diffusion through medium‐ and fine‐sized pores. To derive the soil redox status based upon the triplet ε–O〈sub〉2〈/sub〉–E〈sub〉H〈/sub〉 is challenging at present in heterogeneous soil domains and larger soil volumes than 250 cm〈sup〉3〈/sup〉. Undisturbed soil sampling along with 2D/3D redox measurement systems (e.g., redoxtrons) improve our understanding of redox dynamics within the capillary fringe.

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈sec xmlns:mml="http://www.w3.org/1998/Math/MathML" id="ejss13362-sec-1003" xml:lang="en"〉 〈p xml:lang="en"〉Long‐term experiments (LTEs) have provided data to modellers and agronomists to investigate changes and dynamics of soil organic carbon (SOC) under different cropping systems. As treatment changes have occurred due to agricultural advancements, so too have analytical soil methods. This may lead to method bias over time, which could affect the robust interpretation of data and conclusions drawn. This study aims to quantify differences in SOC due to changes in dry combustion methods over time, using soil samples of a LTE established in 1963 that focuses on mineral and organic fertilizer management in the temperate zone of Northeast Germany. For this purpose, 1059 soil samples, collected between 1976 and 2008, have been analysed twice, once with their historical laboratory method right after sampling, and a second time in 2016 when all samples were analysed using the same elementary analyser. In 9 of 11 soil sampling campaigns, a paired 〈italic toggle="no"〉t〈/italic〉‐test provided evidence for significant differences in the historical SOC values when compared with the re‐analysed concentrations of the same LTE sample. In the sampling years 1988 and 2004, the historical analysis obtained about 0.9 g kg〈sup〉−1〈/sup〉 lower SOC compared with the re‐analysed one. For 1990 and 1998, this difference was about 0.4 g kg〈sup〉−1〈/sup〉. Correction factors, an approach often used to correct for different analytical techniques, could only be applied for 5 of 11 sampling campaigns to account for constant and proportional systematic method error. For this particular LTE, the interpretation of SOC changes due to agronomic management (here fertilization) deviates depending on the analytical method used, which may weaken the explanatory power of the historical data. We demonstrate that analytical method changes over time present one of many challenges in the interpretation of time series data of SOC dynamics. Therefore, LTE site managers need to ensure providing all necessary protocols and data in order to retrace method changes and if necessary recalculate SOC.〈/p〉 〈/sec〉〈sec xmlns:mml="http://www.w3.org/1998/Math/MathML" id="ejss13362-sec-0003" xml:lang="en"〉 〈title〉Highlights〈/title〉 〈p xml:lang="en"〉〈list list-type="bullet" id="ejss13362-list-0001"〉 〈list-item id="ejss13362-li-0001"〉〈p〉A total of 1059 LTE soil samples taken between 1976 and 2008 were re‐analysed for SOC in 2016〈/p〉〈/list-item〉 〈list-item id="ejss13362-li-0002"〉〈p〉Several methodological changes for SOC determination led to significant different SOC concentration in the same sample〈/p〉〈/list-item〉 〈list-item id="ejss13362-li-0003"〉〈p〉Interpretation and time series of LTE soil data suffer from consideration of analytical method changes and poor documentation of the same〈/p〉〈/list-item〉 〈list-item id="ejss13362-li-0004"〉〈p〉Soil archive establishment, thorough method protocols and diligent proficiency testing after soil method changes ameliorate the dilemma〈/p〉〈/list-item〉 〈/list〉〈/p〉 〈/sec〉

Description: Brandenburger Staatsministerium für Wissenschaft, Forschung und Kultur http://dx.doi.org/10.13039/501100004581

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100004937

Description: https://doi.org/10.4228/zalf-acge-b683

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Flood risk assessments require different disciplines to understand and model the underlying components hazard, exposure, and vulnerability. Many methods and data sets have been refined considerably to cover more details of spatial, temporal, or process information. We compile case studies indicating that refined methods and data have a considerable effect on the overall assessment of flood risk. But are these improvements worth the effort? The adequate level of detail is typically unknown and prioritization of improvements in a specific component is hampered by the lack of an overarching view on flood risk. Consequently, creating the dilemma of potentially being too greedy or too wasteful with the resources available for a risk assessment. A “sweet spot” between those two would use methods and data sets that cover all relevant known processes without using resources inefficiently. We provide three key questions as a qualitative guidance toward this “sweet spot.” For quantitative decision support, more overarching case studies in various contexts are needed to reveal the sensitivity of the overall flood risk to individual components. This could also support the anticipation of unforeseen events like the flood event in Germany and Belgium in 2021 and increase the reliability of flood risk assessments.〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: BMBF http://dx.doi.org/10.13039/501100002347

Description: Federal Environment Agency http://dx.doi.org/10.13039/501100010809

Description: http://howas21.gfz-potsdam.de/howas21/

Description: https://www.umwelt.niedersachsen.de/startseite/themen/wasser/hochwasser_amp_kustenschutz/hochwasserrisikomanagement_richtlinie/hochwassergefahren_und_hochwasserrisikokarten/hochwasserkarten-121920.html

Description: https://download.geofabrik.de/europe/germany.html

Description: https://emergency.copernicus.eu/mapping/list-of-components/EMSN024

Description: https://data.jrc.ec.europa.eu/collection/id-0054

Description: https://oasishub.co/dataset/surface-water-flooding-footprinthurricane-harvey-august-2017-jba

Description: https://www.wasser.sachsen.de/hochwassergefahrenkarte-11915.html

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The increasing demand for biomass for food, animal feed, fibre and bioenergy requires optimization of soil productivity, while at the same time, protecting other soil functions such as nutrient cycling and buffering, carbon storage, habitat for biological activity and water filter and storage. Therefore, one of the main challenges for sustainable agriculture is to produce high yields while maintaining all the other soil functions. Mechanistic simulation models are an essential tool to fully understand and predict the complex interactions between physical, biological and chemical processes of soils that generate those functions. We developed a soil model to simulate the impact of various agricultural management options and climate change on soil functions by integrating the relevant processes mechanistically and in a systemic way. As a special feature, we include the dynamics of soil structure induced by tillage and biological activity, which is especially relevant in arable soils. The model operates on a 1D soil profile consisting of a number of discrete layers with dynamic thickness. We demonstrate the model performance by simulating crop growth, root growth, nutrient and water uptake, nitrogen cycling, soil organic matter turnover, microbial activity, water distribution and soil structure dynamics in a long‐term field experiment including different crops and different types and levels of fertilization. The model is able to capture essential features that are measured regularly including crop yield, soil organic carbon, and soil nitrogen. In this way, the plausibility of the implemented processes and their interactions is confirmed. Furthermore, we present the results of explorative simulations comparing scenarios with and without tillage events to analyse the effect of soil structure on soil functions. Since the model is process‐based, we are confident that the model can also be used to predict quantities that have not been measured or to estimate the effect of management measures and climate states not yet been observed. The model thus has the potential to predict the site‐specific impact of management decisions on soil functions, which is of great importance for the development of a sustainable agriculture that is currently also on the agenda of the ‘Green Deal’ at the European level.〈/p〉

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: https://git.ufz.de/bodium/bodium_v1.0

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Gas transport in soils is usually assumed to be purely diffusive, although several studies have shown that non‐diffusive processes can significantly enhance soil gas transport. These processes include barometric air pressure changes, wind‐induced pressure pumping and static air pressure fields generated by wind interacting with obstacles. The associated pressure gradients in the soil can cause advective gas fluxes that are much larger than diffusive fluxes. However, the contributions of the respective transport processes are difficult to separate. We developed a large chamber system to simulate pressure fields and investigate their influence on soil gas transport. The chamber consists of four subspaces in which pressure is regulated by fans that blow air in or out of the chamber. With this setup, we conducted experiments with oscillating and static pressure fields. CO〈sub〉2〈/sub〉 concentrations were measured along two soil profiles beneath the chamber. We found a significant relationship between static lateral pressure gradients and the change in the CO〈sub〉2〈/sub〉 profiles (R〈sup〉2〈/sup〉 = 0.53; 〈italic toggle="no"〉p〈/italic〉‐value 〈2e‐16). Even small pressure gradients between −1 and 1 Pa relative to ambient pressure resulted in an increase or decrease in CO〈sub〉2〈/sub〉 concentrations of 8% on average in the upper soil, indicating advective flow of air in the pore space. Positive pressure gradients resulted in decreasing, negative pressure gradients in increasing CO〈sub〉2〈/sub〉 concentrations. The concentration changes were probably caused by an advective flow field in the soil beneath the chamber generated by the pressure gradients. No effect of oscillating pressure fields was observed in this study. The results indicate that static lateral pressure gradients have a substantial impact on soil gas transport and therefore are an important driver of gas exchange between soil and atmosphere. Lateral pressure gradients in a comparable range can be induced under windy conditions when wind interacts with terrain features. They can also be caused by chambers used for flux measurements at high wind speed or by fans used for head‐space mixing within the chambers, which yields biased flux estimates.〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Temperature and soil moisture are known to control pesticide mineralization. Half‐life times (DT〈sub〉50〈/sub〉) derived from pesticide mineralization curves generally indicate longer residence times at low soil temperature and moisture but do not consider potential changes in the microbial allocation of pesticide‐derived carbon (C). We aimed to determine carbon use efficiency (CUE, formation of new biomass relative to total C uptake) to better understand microbial utilization of pesticide‐derived C under different environmental conditions and to support the conventional description of degradation dynamics based on mineralization. We performed a microcosm experiment at two MCPA (2‐methyl‐4‐chlorophenoxyacetic acid) concentrations (1 and 20 mg kg〈sup〉−1〈/sup〉) and defined 20°C/pF 1.8 as optimal and 10°C/pF 3.5 as limiting environmental conditions. After 4 weeks, 70% of the initially applied MCPA was mineralized under optimal conditions but MCPA mineralization reached less than 25% under limiting conditions. However, under limiting conditions, an increase in CUE was observed, indicating a shift towards anabolic utilization of MCPA‐derived C. In this case, increased C assimilation implied C storage or the formation of precursor compounds to support resistance mechanisms, rather than actual growth since we did not find an increase in the 〈italic toggle="no"〉tfdA〈/italic〉 gene relevant to MCPA degradation. We were able to confirm the assumption that under limiting conditions, C assimilation increases relative to mineralization and that C redistribution, may serve as an explanation for the difference between mineralization and MCPA dissipation‐derived degradation dynamics. In addition, by introducing CUE to the temperature‐ and moisture‐dependent degradation of pesticides, we can capture the underlying microbial constraints and adaptive mechanisms to changing environmental conditions.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Changing environmental conditions alter the MCPA degradation dynamics and the allocation of pesticide‐derived carbon to anabolic or catabolic metabolism.〈boxed-text position="anchor" content-type="graphic" id="ejss13417-blkfxd-0001" xml:lang="en"〉 〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:13510754:media:ejss13417:ejss13417-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: Collaborative Research Center 1253 CAMPOS (DFG)

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: DFG Priority Program 2322 “Soil System”

Description: Ellrichshausen Foundation

Description: Research Training Group “Integrated Hydrosystem modeling”

Description: https://doi.org/10.5281/zenodo.5081655

Description: Charcoal‐rich Technosols on century‐old relict charcoal hearths (RCHs) are the subject of ongoing research regarding potential legacy effects that result from historic charcoal production and subsequent charcoal amendments on forest soil properties and forest ecosystems today. RCHs consist mostly of Auh horizons that are substantially enriched in soil organic carbon (SOC), of which the largest part seems to be of pyrogenic origin (PyC). However, the reported range of SOC and PyC contents in RCH soil also suggests that they are enriched in nonpyrogenic SOC. RCH soils are discussed as potential benchmarks for the long‐term influence of biochar amendment and the post‐wildfire influences on soil properties. In this study, we utilised a large soil sample dataset (n = 1245) from 52 RCH sites in north‐western Connecticut, USA, to quantify SOC contents by total element analysis. The contents of condensed highly aromatic carbon as a proxy for black carbon (BC) were predicted by using a modified benzene polycarboxylated acid (BPCA) marker method in combination with diffuse reflectance infrared Fourier transform (DRIFT) spectroscopy‐based partial least square regression (r2 = 0.89). A high vertical spatial sampling resolution allowed the identification of soil organic matter (SOM) enrichment and translocation processes. The results show an average 75% and 1862% increase in TOC and BPCA‐derived carbon, respectively, for technogenic Auh horizons compared to reference soils. In addition to an increase in aromatic properties, increased carboxylic properties of the RCH SOC suggest self‐humification effects of degrading charcoal and thereby the continuing formation of leachable aromatic carbon compounds, which could have effects on pedogenic processes in buried soils. Indeed, we show BPCA‐derived carbon concentrations in intermediate technogenic Cu horizons and buried top/subsoils that suggest vertical translocation of highly aromatic carbon originating in RCH Auh horizons. Topmost Auh horizons showed a gradual decrease in total organic carbon (TOC) contents with increasing depth, suggesting accumulation of recent, non‐pyrogenic SOM. Lower aliphatic absorptions in RCH soil spectra suggest different SOM turnover dynamics compared to reference soils. Furthermore, studied RCH soils featured additional TOC enrichment, which cannot be fully explained now. Highlights BC to TOC ratio and high resolution vertical SOC distribution in 52 RCH sites were studied. RCH soils non‐BC pool was potentially different to reference soils. RCH soils feature TOC accumulation in the topmost horizon. There is BC translocation into buried soils on RCH sites.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Little research attention has been given to validating clusters obtained from the groundwater geochemistry of the waterworks' capture zone with a prevailing lake‐groundwater exchange. To address this knowledge gap, we proposed a new scheme whereby Gaussian finite mixture modeling (GFMM) and Spike‐and‐Slab Bayesian (SSB) algorithms were utilized to cluster the groundwater geochemistry while quantifying the probability of the resulting cluster membership against each other. We applied GFMM and SSB to 13 geochemical parameters collected during different sampling periods at 13 observation points across the Barnim Highlands plateau located in the northeast of Berlin, Germany; this included 10 observation wells, two lakes, and a gallery of drinking production wells. The cluster analysis of GFMM yielded nine clusters, either with a probability ≥0.8, while the SSB produced three hierarchical clusters with a probability of cluster membership varying from 〈0.2 to 〉0.8. The findings demonstrated that the clustering results of GFMM were in good agreement with the classification as per the principal component analysis and Piper diagram. By superimposing the parameter clustering onto the observation clustering, we could identify discrepancies that exist among the parameters of a certain cluster. This enables the identification of different factors that may control the geochemistry of a certain cluster, although parameters of that cluster share a strong similarity. The GFMM results have shown that from 2002, there has been active groundwater inflow from the lakes towards the capture zone. This means that it is necessary to adopt appropriate measures to reverse the inflow towards the lakes.

Description: Article impact statement: The probability of cluster membership quantified using an algorithm should be validated against another probabilistic‐based classifier.

Description: Federal Ministry of Education and Research http://dx.doi.org/10.13039/501100002347

Description: Copper (Cu) is an essential element for plants and microorganisms and at larger concentrations a toxic pollutant. A number of factors controlling Cu dynamics have been reported, but information on quantitative relationships is scarce. We aimed to (i) quantitatively describe and predict soil Cu concentrations (CuAR) in aqua regia considering site‐specific effects and effects of pH, soil organic carbon (SOC) and cation exchange capacity (CEC), and (ii) study the suitability of mixed‐effects modelling and rule‐based models for the analysis of long‐term soil monitoring data. Thirteen uncontaminated long‐term monitoring soil profiles in southern Germany were analysed. Since there was no measurable trend of increasing CuAR concentrations with time in the respective depth ranges of the sites, data from different sampling dates were combined and horizon‐specific regression analyses including model simplifications were carried out for 10 horizons. Fixed‐ and mixed‐effects models with the site as a random effect were useful for the different horizons and significant contributions (either of main effects or interactions) of SOC, CEC and pH were present for 9, 8 and 7 horizons, respectively. Horizon‐specific rule‐based cubist models described the CuAR data similarly well. Validations of cubist models and mixed‐effects models for the CuAR concentrations in A horizons were successful for the given population after random splitting into calibration and validation samples, but not after independent validations with random splitting according to sites. Overall, site, CEC, SOC and pH provide important information for a description of CuAR concentrations using the different regression approaches. Highlights: Information on quantitative relationships for factors controlling Cu dynamics is scarce. Site, CEC, SOC and pH provide important information for a description of Cu concentrations. Validations of cubist models and mixed‐effects models for A horizons were successful for a closed population of sites.

Description: Bavarian State Ministry of the Environment and Consumer Protection http://dx.doi.org/10.13039/501100010219

Description: Ministry of Agriculture and Environment Mecklenburg‐Western Pomerania

Description: Soil aeration is a critical factor for oxygen‐limited subsoil processes, as transport by diffusion and advection is restricted by the long distance to the free atmosphere. Oxygen transport into the soil matrix is highly dependent on its connectivity to larger pore channels like earthworm and root colonised biopores. Here we hypothesize that the soil matrix around biopores represents different connectivity depending on biopore genesis and actual coloniser. We analysed the soil pore system of undisturbed soil core samples around biopores generated or colonised by roots and earthworms and compared them with the pore system of soil, not in the immediacy of a biopore. Oxygen partial pressure profiles and gas relative diffusion was measured in the rhizosphere and drilosphere from the biopore wall into the bulk soil with microelectrodes. The measurements were linked with structural features such as porosity and connectivity obtained from X‐ray tomography and image analysis. Aeration was enhanced in the soil matrix surrounding biopores in comparison to the bulk soil, shown by higher oxygen concentrations and higher relative diffusion coefficients. Biopores colonised by roots presented more connected lateral pores than earthworm colonised ones, which resulted in enhanced aeration of the rhizosphere compared to the drilosphere. This has influenced biotic processes (microbial turnover/mineralization or root respiration) at biopore interfaces and highlights the importance of microstructural features for soil processes and their dependency on the biopore's coloniser.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Application of farmyard manure (FYM) is common practice to improve physical and chemical properties of arable soil and crop yields. However, studies on effects of FYM application mainly focussed on topsoils, whereas subsoils have rarely been addressed so far. We, therefore, investigated the effects of 36‐year FYM application with different rates of annual organic carbon (OC) addition (0, 469, 938 and 1875 g C m−2 a−1) on OC contents of a Chernozem in 0–30 cm (topsoil) and 35–45 cm (subsoil) depth. We also investigated its effects on soil structure and hydraulic properties in subsoil. X‐ray computed tomography was used to analyse the response of the subsoil macropore system (≥19 μm) and the distribution of particulate organic matter (POM) to different FYM applications, which were related to contents in total OC (TOC) and water‐extractable OC (WEOC). We show that FYM‐C application of 469 g C m−2 a−1 caused increases in TOC and WEOC contents only in the topsoil, whereas rates of ≥938 g C m−2 a−1 were necessary for TOC enrichment also in the subsoil. At this depth, the subdivision of TOC into different OC sources shows that most of the increase was due to fresh POM, likely by the stimulation of root growth and bioturbation. The increase in subsoil TOC went along with increases in macroporosity and macropore connectivity. We neither observed increases in plant‐available water capacity nor in unsaturated hydraulic conductivity. In conclusion, only very high application of FYM over long periods can increase OC content of subsoil at our study site, but this increase is largely based on fresh, easily degradable POM and likely accompanied by high C losses when considering the discrepancy between OC addition rate by FYM and TOC response in soil. Highlights A new image processing procedure to distinguish fresh and decomposed POM. The increase of subsoil C stock based to a large extend on fresh, labile POM. Potential of arable subsoils for long‐term C storage by large FYM application rates is limited. The increase in TOC has no effect on hydraulic properties of the subsoil.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Erosion is a severe threat to the sustainable use of agricultural soils. However, the structural resistance of soil against the disruptive forces steppe soils experience under field conditions has not been investigated. Therefore, 132 topsoils under grass‐ and cropland covering a large range of physico‐chemical soil properties (sand: 2–76%, silt: 18–80%, clay: 6–30%, organic carbon: 7.3–64.2 g kg−1, inorganic carbon: 0.0–8.5 g kg−1, pH: 4.8–9.5, electrical conductivity: 32–946 μS cm−1) from northern Kazakhstan were assessed for their potential erodibility using several tests. An adjusted drop‐shatter method (low energy input of 60 Joule on a 250‐cm3 soil block) was used to estimate the stability of dry soil against weak mechanical forces, such as saltating particles striking the surface causing wind erosion. Three wetting treatments with various conditions and energies (fast wetting, slow wetting, and wet shaking) were applied to simulate different disruptive effects of water. Results indicate that aggregate stability was higher for grassland than cropland soils and declined with decreasing soil organic carbon content. The results of the drop‐shatter test suggested that 29% of the soils under cropland were at risk of wind erosion, but only 6% were at high risk (i.e. erodible fraction 〉60%). In contrast, the fast wetting treatment revealed that 54% of the samples were prone to become “very unstable” and 44% “unstable” during heavy rain or snowmelt events. Even under conditions comparable to light rain events or raindrop impact, 53–59% of the samples were “unstable.” Overall, cropland soils under semi‐arid conditions seem much more susceptible to water than wind erosion. Considering future projections of increasing precipitation in Kazakhstan, we conclude that the risk of water erosion is potentially underestimated and needs to be taken into account when developing sustainable land use strategies. Highlights Organic matter is the important binding agent enhancing aggregation in steppe topsoils. Tillage always declines aggregate stability even without soil organic carbon changes. All croplands soil are prone to wind or water erosion independent of their soil properties. Despite the semi‐arid conditions, erosion risk by water seems higher than by wind.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Stable hydrogen isotope ratios (δ2H values) in structural hydroxyl groups of pedogenic clay minerals are inherited from the surrounding water at the time of their formation. Only non‐exchangeable H preserves the environmental forensic and paleoclimate information (δ2Hn value). To measure δ2Hn values in structural H of clay minerals and soil clay fractions, we adapted a steam equilibration method by accounting for high hygroscopicity. Our δ2Hn values for USGS57 biotite (−95.3 ± SD 0.9‰) and USGS58 muscovite (30.7 ± 1.4‰) differed slightly but significantly from the reported δ2H values (−91.5 ± 2.4‰ and −28.4 ± 1.6‰), because the minerals contained 1.1%–4.4% of exchangeable H. The low SD of replicate measurements (n = 3) confirmed a high precision. The clay separation method including destruction of Fe oxides, carbonates and soil organic matter, and dispersion did not significantly change the δ2Hn values of five different clay minerals. However, we were unable to remove all organic matter from the soil clay fractions resulting in an estimated bias of 1‰ in two samples and 15‰ in the carbon‐richest sample. Our results demonstrate that δ2Hn values of structural H of clay minerals and soil clay fractions can be reliably measured without interference from atmospheric water and the method used to separate the soil clay fraction. Highlights We tested steam equilibration to determine stable isotope ratios of structural H in clay. Gas‐tight capsule sealing in Ar atmosphere was necessary to avoid remoistening. Our steam equilibration method showed a high accuracy and precision. The clay separation method did not change stable isotope ratios of structural H in clay.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Soil fauna drives crucial processes of energy and nutrient cycling in agricultural systems, and influences the quality of crops and pest incidence. Soil tillage is the most influential agricultural manipulation of soil structure, and has a profound influence on soil biology and its provision of ecosystem services. The objective of this study was to quantify through meta‐analyses the effects of reducing tillage intensity on density and diversity of soil micro‐ and mesofaunal communities, and how these effects vary among different pedoclimatic conditions and interact with concurrent management practices. We present the results of a global meta‐analysis of available literature data on the effects of different tillage intensities on taxonomic and functional groups of soil micro‐ and mesofauna. We collected paired observations (conventional vs. reduced forms of tillage/no‐tillage) from 133 studies across 33 countries. Our results show that reduced tillage intensity or no‐tillage increases the total density of springtails (+35%), mites (+23%), and enchytraeids (+37%) compared to more intense tillage methods. The meta‐analyses for different nematode feeding groups, life‐forms of springtails, and taxonomic mite groups showed higher densities under reduced forms of tillage compared to conventional tillage on omnivorous nematodes (+53%), epedaphic (+81%) and hemiedaphic (+84%) springtails, oribatid (+43%) and mesostigmatid (+57%) mites. Furthermore, the effects of reduced forms of tillage on soil micro‐ and mesofauna varied with depth, climate and soil texture, as well as with tillage method, tillage frequency, concurrent fertilisation, and herbicide application. Our findings suggest that reducing tillage intensity can have positive effects on the density of micro‐ and mesofaunal communities in areas subjected to long‐term intensive cultivation practices. Our results will be useful to support decision making on the management of soil faunal communities and will facilitate modelling efforts of soil biology in global agroecosystems. HIGHLIGHTS Global meta‐analysis to estimate the effect of reducing tillage intensity on micro‐ and mesofauna Reduced tillage or no‐tillage has positive effects on springtail, mite and enchytraeid density Effects vary among nematode feeding groups, springtail life forms and mite suborders Effects vary with texture, climate and depth and depend on the tillage method and frequency

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: https://doi.org/10.20387/bonares-eh0f-hj28

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Hydrogeological information about an aquifer is difficult and costly to obtain, yet essential for the efficient management of groundwater resources. Transferring information from sampled sites to a specific site of interest can provide information when site‐specific data is lacking. Central to this approach is the notion of site similarity, which is necessary for determining relevant sites to include in the data transfer process. In this paper, we present a data‐driven method for defining site similarity. We apply this method to selecting groups of similar sites from which to derive prior distributions for the Bayesian estimation of hydraulic conductivity measurements at sites of interest. We conclude that there is now a unique opportunity to combine hydrogeological expertise with data‐driven methods to improve the predictive ability of stochastic hydrogeological models.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉〈italic〉Article impact statement〈/italic〉: This article introduces hierarchical clustering as a method for defining a notion of site similarity; the aim of this method is to improve the derivation of prior distributions in Bayesian methods in hydrogeology.〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://github.com/GeoStat-Bayesian/geostatDB

Description: https://github.com/GeoStat-Bayesian/exPrior

Description: https://github.com/GeoStat-Bayesian/siteSimilarity

Description: Efforts to collaboratively manage the risk of flooding are ultimately based on individuals learning about risks, the decision process, and the effectiveness of decisions made in prior situations. This article argues that much can be learned about a governance setting by explicitly evaluating the relationships through which influential individuals and their immediate contacts receive and send information to one another. We define these individuals as “brokers,” and the networks that emerge from their interactions as “learning spaces.” The aim of this article is to develop strategies to identify and evaluate the properties of a broker's learning space that are indicative of a collaborative flood risk management arrangement. The first part of this article introduces a set of indicators, and presents strategies to employ this list so as to systematically identify brokers, and compare their learning spaces. The second part outlines the lessons from an evaluation that explored cases in two distinct flood risk management settings in Germany. The results show differences in the observed brokers' learning spaces. The contacts and interactions of the broker in Baden‐Württemberg imply a collaborative setting. In contrast, learning space of the broker in North Rhine‐Westphalia lacks the same level of diversity and polycentricity.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: MWK Baden‐Württemberg

Description: Most common machine learning (ML) algorithms usually work well on balanced training sets, that is, datasets in which all classes are approximately represented equally. Otherwise, the accuracy estimates may be unreliable and classes with only a few values are often misclassified or neglected. This is known as a class imbalance problem in machine learning and datasets that do not meet this criterion are referred to as imbalanced data. Most datasets of soil classes are, therefore, imbalanced data. One of our main objectives is to compare eight resampling strategies that have been developed to counteract the imbalanced data problem. We compared the performance of five of the most common ML algorithms with the resampling approaches. The highest increase in prediction accuracy was achieved with SMOTE (the synthetic minority oversampling technique). In comparison to the baseline prediction on the original dataset, we achieved an increase of about 10, 20 and 10% in the overall accuracy, kappa index and F‐score, respectively. Regarding the ML approaches, random forest (RF) showed the best performance with an overall accuracy, kappa index and F‐score of 66, 60 and 57%, respectively. Moreover, the combination of RF and SMOTE improved the accuracy of the individual soil classes, compared to RF trained on the original dataset and allowed better prediction of soil classes with a low number of samples in the corresponding soil profile database, in our case for Chernozems. Our results show that balancing existing soil legacy data using synthetic sampling strategies can significantly improve the prediction accuracy in digital soil mapping (DSM). Highlights Spatial distribution of soil classes in Iran can be predicted using machine learning (ML) algorithms. The synthetic minority oversampling technique overcomes the drawback of imbalanced and highly biased soil legacy data. When combining a random forest model with synthetic sampling strategies the prediction accuracy of the soil model improves significantly. The resulting new soil map of Iran has a much higher spatial resolution compared to existing maps and displays new soil classes that have not yet been mapped in Iran.

Description: Alexander von Humboldt‐Stiftung http://dx.doi.org/10.13039/100005156

Description: German Research Foundation http://dx.doi.org/10.13039/501100001659

Description: Soil and Water Research Institute, Agricultural Research, Education and Extension Organization, Karaj, Iran

Description: The application of biochar to agricultural soils to increase nutrient availability, crop production and carbon sequestration has gained increasing interest but data from field experiments on temperate, marginal soils are still under‐represented. In the current study, biochar, produced from organic residues (digestates) from a biogas plant, was applied with and without digestates at low (3.4 t ha−1) and intermediate (17.1 t ha−1) rates to two acidic and sandy soils in northern Germany that are used for corn (Zea mays L.) production. Soil nutrient availability, crop yields, microbial biomass and carbon dioxide (CO2) emissions from heterotrophic respiration were measured over two consecutive years. The effects of biochar application depended on the intrinsic properties of the two tested soils and the biochar application rates. Although the soils at the fallow site, with initially low nutrient concentrations, showed a significant increase in pH, soil nutrients and crop yield after low biochar application rates, a similar response was found at the cornfield site only after application of substantially larger amounts of biochar. The effect of a single dose of biochar at the beginning of the experiment diminished over time but was still detectable after 2 years. Whereas plant available nutrient concentrations increased after biochar application, the availability of potentially phytotoxic trace elements (Zn, Pb, Cd, Cr) decreased significantly, and although slight increases in microbial biomass carbon and heterotrophic CO2 fluxes were observed after biochar application, they were mostly not significant. The results indicate that the application of relatively small amounts of biochar could have positive effects on plant available nutrients and crop yields of marginal arable soils and may decrease the need for mineral fertilizers while simultaneously increasing the sequestration of soil organic carbon. Highlights A low rate of biochar increased plant available nutrients and crop yield on marginal soils. Biochar application reduced the availability of potentially harmful trace elements. Heterotrophic respiration showed no clear response to biochar application. Biochar application may reduce fertilizer need and increase carbon sequestration on marginal soils.

Description: German Academic Exchange Service http://dx.doi.org/10.13039/501100001655

Description: Institute Strategic Programme grants, “Soils to Nutrition”

Description: Social inequalities lead to flood resilience inequalities across social groups, a topic that requires improved documentation and understanding. The objective of this paper is to attend to these differences by investigating self‐stated flood recovery across genders in Vietnam as a conceptual replication of earlier results from Germany. This study employs a regression‐based analysis of 1,010 respondents divided between a rural coastal and an urban community in Thua Thien‐Hue province. The results highlight an important set of recovery process‐related variables. The set of relevant variables is similar across genders in terms of inclusion and influence, and includes age, social capital, internal and external support after a flood, perceived severity of previous flood impacts, and the perception of stress‐resilience. However, women were affected more heavily by flooding in terms of longer recovery times, which should be accounted for in risk management. Overall, the studied variables perform similarly in Vietnam and Germany. This study, therefore, conceptually replicates previous results suggesting that women display slightly slower recovery levels as well as that psychological variables influence recovery rates more than adverse flood impacts. This provides an indication of the results' potentially robust nature due to the different socio‐environmental contexts in Germany and Vietnam.

Description: Nitrogen (N) fertilization is the major contributor to nitrous oxide (N2O) emissions from agricultural soil, especially in post‐harvest seasons. This study was carried out to investigate whether ryegrass serving as cover crop affects soil N2O emissions and denitrifier community size. A microcosm experiment was conducted with soil planted with perennial ryegrass (Lolium perenne L.) and bare soil, each with four levels of N fertilizer (0, 5, 10 and 20 g N m−2; applied as calcium ammonium nitrate). The closed‐chamber approach was used to measure soil N2O fluxes. Real‐time PCR was used to estimate the biomass of bacteria and fungi and the abundance of genes involved in denitrification in soil. The results showed that the presence of ryegrass decreased the nitrate content in soil. Cumulative N2O emissions of soil with grass were lower than in bare soil at 5 and 10 g N m−2. Fertilization levels did not affect the abundance of soil bacteria and fungi. Soil with grass showed greater abundances of bacteria and fungi, as well as microorganisms carrying narG, napA, nirK, nirS and nosZ clade I genes. It is concluded that ryegrass serving as a cover crop holds the potential to mitigate soil N2O emissions in soils with moderate or high NO3− concentrations. This highlights the importance of cover crops for the reduction of N2O emissions from soil, particularly following N fertilization. Future research should explore the full potential of ryegrass to reduce soil N2O emissions under field conditions as well as in different soils. Highlights This study was to investigate whether ryegrass serving as cover crop affects soil N2O emissions and denitrifier community size; Plant reduced soil N substrates on one side, but their root exudates stimulated denitrification on the other side; N2O emissions were lower in soil with grass than bare soil at medium fertilizer levels, and growing grass stimulated the proliferation of almost all the denitrifying bacteria except nosZ clade II; Ryegrass serving as a cover crop holds the potential to mitigate soil N2O emissions.

Description: China Scholarship Council http://dx.doi.org/10.13039/501100004543

Description: The National Science Project for University of Anhui Province

Description: High‐performance numerical codes are an indispensable tool for hydrogeologists when modeling subsurface flow and transport systems. But as they are written in compiled languages, like C/C++ or Fortran, established software packages are rarely user‐friendly, limiting a wider adoption of such tools. OpenGeoSys (OGS), an open‐source, finite‐element solver for thermo‐hydro‐mechanical–chemical processes in porous and fractured media, is no exception. Graphical user interfaces may increase usability, but do so at a dramatic reduction of flexibility and are difficult or impossible to integrate into a larger workflow. Python offers an optimal trade‐off between these goals by providing a highly flexible, yet comparatively user‐friendly environment for software applications. Hence, we introduce ogs5py, a Python‐API for the OpenGeoSys 5 scientific modeling package. It provides a fully Python‐based representation of an OGS project, a large array of convenience functions for users to interact with OGS and connects OGS to the scientific and computational environment of Python.

Description: German Federal Environmental Foundation http://dx.doi.org/10.13039/100007636

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment–water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier–Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in‐ and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high‐resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions.

Description: Technische Universität Berlin, Germany

Description: German Research Foundation http://dx.doi.org/10.13039/501100001659

Description: Sustainable arable cropping relies on repeated liming. Yet, the associated increase in soil pH can reduce the availability of iron (Fe) to plants. We hypothesized that repeated liming, but not pedogenic processes such as lessivage (i.e., translocation of clay particles), alters the Fe cycle in Luvisol soil, thereby affecting Fe isotope composition in soils and crops. Hence, we analysed Fe concentrations and isotope compositions in soil profiles and winter rye from the long‐term agricultural experimental site in Berlin‐Dahlem, Germany, where a controlled liming trial with three field replicates per treatment has been conducted on Albic Luvisols since 1923. Heterogeneity in subsoil was observed at this site for Fe concentration but not for Fe isotope composition. Lessivage had not affected Fe isotope composition in the soil profiles. The results also showed that almost 100 years of liming lowered the concentration of the HCl‐extractable Fe that was potentially available for plant uptake in the surface soil (0–15 cm) from 1.03 (standard error (SE) 0.03) to 0.94 (SE 0.01) g kg−1. This HCl‐extractable Fe pool contained isotopically lighter Fe (δ56Fe = −0.05 to −0.29‰) than the bulk soil (δ56Fe = −0.08 to 0.08‰). However, its Fe isotope composition was not altered by the long‐term lime application. Liming resulted in relatively lower Fe concentrations in the roots of winter rye. In addition, liming led to a heavier Fe isotope composition of the whole plants compared with those grown in the non‐limed plots (δ56FeWholePlant_ + Lime = −0.12‰, SE 0.03 vs. δ56FeWholePlant_‐Lime = −0.21‰, SE 0.01). This suggests that the elevated soil pH (increased by one unit due to liming) promoted the Fe uptake strategy through complexation of Fe(III) from the rhizosphere, which favoured heavier Fe isotopes. Overall, the present study showed that liming and a related increase in pH did not affect the Fe isotope compositions of the soil, but may influence the Fe isotope composition of plants grown in the soil if they alter their Fe uptake strategy upon the change of Fe availability. Highlights Fe concentrations and stocks, but not Fe isotope compositions, were more heterogeneous in subsoil than in topsoil. Translocation of clay minerals did not result in Fe isotope fractionation in the soil profile of a Luvisol. Liming decreased Fe availability in topsoil, but did not affect its δ56Fe values. Uptake of heavier Fe isotopes by graminaceous crops was more pronounced at elevated pH.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: We present a workflow to estimate geostatistical aquifer parameters from pumping test data using the Python package welltestpy. The procedure of pumping test analysis is exemplified for two data sets from the Horkheimer Insel site and from the Lauswiesen site, Germany. The analysis is based on a semi‐analytical drawdown solution from the upscaling approach Radial Coarse Graining, which enables to infer log‐transmissivity variance and horizontal correlation length, beside mean transmissivity, and storativity, from pumping test data. We estimate these parameters of aquifer heterogeneity from type‐curve analysis and determine their sensitivity. This procedure, implemented in welltestpy, is a template for analyzing any pumping test. It goes beyond the possibilities of standard methods, for example, based on Theis' equation, which are limited to mean transmissivity and storativity. A sensitivity study showed the impact of observation well positions on the parameter estimation quality. The insights of this study help to optimize future test setups for geostatistical aquifer analysis and provides guidance for investigating pumping tests with regard to aquifer statistics using the open‐source software package welltestpy.

Description: Article impact statement: We present a workflow to infer parameters of subsurface heterogeneity from pumping test data exemplified at two sites using welltestpy.

Description: German Federal Environmental Foundation (DBU) http://dx.doi.org/10.13039/100007636

Description: In designed experiments, different sources of variability and an adequate scale of measurement need to be considered, but not all approaches in common usage are equally valid. In order to elucidate the importance of sources of variability and choice of scale, we conducted an experiment where the effects of biochar and slurry applications on soil properties related to soil fertility were studied for different designs: (a) for a field‐scale sampling design with either a model soil (without natural variability) as an internal control or with composited soils, (b) for a design with a focus on amendment variabilities, and (c) for three individual field‐scale designs with true field replication and a combined analysis representative of the population of loess‐derived soils. Three silty loam sites in Germany were sampled and the soil macroaggregates were crushed. For each design, six treatments (0, 0.15 and 0.30 g slurry‐N kg−1 with and without 30 g biochar kg−1) were applied before incubating the units under constant soil moisture conditions for 78 days. CO2 fluxes were monitored and soils were analysed for macroaggregate yields and associated organic carbon (C). Mixed‐effects models were used to describe the effects. For all soil properties, results for the loess sites differed with respect to significant contributions of fixed effects for at least one site, suggesting the need for a general inclusion of different sites. Analysis using a multilevel model allowed generalizations for loess soils to be made and showed that site:slurry:biochar and site:slurry interactions were not negligible for macroaggregate yields. The use of a model soil as an internal control enabled observation of variabilities other than those related to soils or amendments. Experiments incorporating natural variability in soils or amendments resulted in partially different outcomes, indicating the need to include all important sources of variability. Highlights Effects of biochar and slurry applications were studied for different designs and mixed‐effects models were used to describe the effects. Including an internal control allowed observation of, e.g., methodological and analytical variabilities. The results suggested the need for a general inclusion of different sites. Analysis using a multilevel model allowed generalizations for loess soils. The results indicated the need to include all important sources of variability.

Description: Temperate forest soils are often considered as an important sink for atmospheric carbon (C), thereby buffering anthropogenic CO2 emissions. However, the effect of tree species composition on the magnitude of this sink is unclear. We resampled a tree species common garden experiment (six sites) a decade after initial sampling to evaluate whether forest floor (FF) and topsoil organic carbon (Corg) and total nitrogen (Nt) stocks changed in dependence of tree species (Norway spruce—Picea abies L., European beech—Fagus sylvatica L., pedunculate oak—Quercus robur L., sycamore maple—Acer pseudoplatanus L., European ash—Fraxinus excelsior L. and small‐leaved lime—Tilia cordata L.). Two groups of species were identified in terms of Corg and Nt distribution: (1) Spruce with high Corg and Nt stocks in the FF developed as a mor humus layer which tended to have smaller Corg and Nt stocks and a wider Corg:Nt ratio in the mineral topsoil, and (2) the broadleaved species, of which ash and maple distinguished most clearly from spruce by very low Corg and Nt stocks in the FF developed as mull humus layer, had greater Corg and Nt stocks, and narrow Corg:Nt ratios in the mineral topsoil. Over 11 years, FF Corg and Nt stocks increased most under spruce, while small decreases in bulk mineral soil (esp. in 0–15 cm and 0–30 cm depth) Corg and Nt stocks dominated irrespective of species. Observed decadal changes were associated with site‐related and tree species‐mediated soil properties in a way that hinted towards short‐term accumulation and mineralisation dynamics of easily available organic substances. We found no indication for Corg stabilisation. However, results indicated increasing Nt stabilisation with increasing biomass of burrowing earthworms, which were highest under ash, lime and maple and lowest under spruce. Highlights We studied if tree species differences in topsoil Corg and Nt stocks substantiate after a decade. The study is unique in its repeated soil sampling in a multisite common garden experiment. Forest floors increased under spruce, but topsoil stocks decreased irrespective of species. Changes were of short‐term nature. Nitrogen was most stable under arbuscular mycorrhizal species.

Description: Deutsche Forschungsgemeinschaff (DFG)

Description: Coping with the growing impacts of flooding in EU countries, a paradigm shift in flood management can be observed, moving from safety‐based towards risk‐based approaches and holistic perspectives. Flood resilience is a common denominator of most of the approaches. In this article, we present the ‘Flood Resilience Rose’ (FRR), a management tool to promote harmonised action towards flood resilience in European regions and beyond. The FRR is a result of a two‐step process. First, based on scientific concepts as well as analysis of relevant policy documents, we identified three ‘levels of operation’. The first level refers to the EU Floods Directive and an extended multi‐layer safety approach, comprising the four different layers of protection, prevention, preparedness and recovery, and related measures to be taken. This level is not independent but depends both on the institutional (second level) and the wider (third level) context. Second, we used surveys, semi‐structured interviews and group discussions during workshops with experts from Belgium, Denmark, Germany, the Netherlands and the United Kingdom to validate the definitions and the FRR's practical relevance. The presented FRR is thus the result of rigorous theoretical and practical consideration and provides a tool capable to strengthen flood risk management practice.

Description: European Regional Development Fund http://dx.doi.org/10.13039/501100008530

Description: In recent years, German cities were heavily impacted by pluvial flooding and related damage is projected to increase due to climate change and urbanisation. It is important to ask how to improve urban pluvial flood risk management. To understand the current state of property level adaptation, a survey was conducted in four municipalities that had recently been impacted by pluvial flooding. A hybrid framework based on the Protection Motivation Theory (PMT) and the Protection Action Decision Model (PADM) was used to investigate drivers of adaptive behaviour through both descriptive and regression analyses. Descriptive statistics revealed that participants tended to instal more low‐ and medium‐cost measures than high‐cost measures. Regression analyses showed that coping appraisal increased protection motivation, but that the adaptive behaviour also depends on framing factors, particularly homeownership. We further found that, while threat appraisal solely affects protection motivation and responsibility appraisal affects solely maladaptive thinking, coping appraisal affects both. Our results indicate that PMT is a solid starting point to study adaptive behaviours in the context of pluvial flooding, but we need to go beyond that by, for instance, considering factors of the PADM, such as responsibility, ownership, or respondent age, to fully understand this complex decision‐making process.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: 〈span〉〈div〉Abstract〈/div〉At the Antamina deposit, Peru, accurate classification of exoskarns and endoskarns can be problematic when textures are mottled. In this study, we use whole-rock geochemical compositions (62 elements) of 221 samples to differentiate texturally similar endoskarns and exoskarns by comparing their compositions to least altered precursors (wall rocks and intrusive rocks). We present a simple method for discriminating these skarn types using immobile element bivariate plots. The most effective discriminators partition endoskarn and exoskarn into distinct domains defined by the composition of each precursor; these include Al〈sub〉2〈/sub〉O〈sub〉3〈/sub〉 versus heavy rare earth elements and some high field strength elements. Using these geochemical parameters, undifferentiated skarn samples can be more reliably classified as endoskarn or exoskarn. The effectiveness of these element pairs is attributed to their significantly different initial concentrations in wall rocks versus igneous precursors and their immobility during skarn formation. While immobile elements can differentiate the skarns, mobile element gains and losses (quantified using isocon analysis) provide insight on the bulk mineralogical and mass changes that take place during skarn formation.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Tongkuangyu copper deposit in the Zhongtiaoshan region, southern Trans-North China orogen, is hosted by a poorly constrained sequence of Paleoproterozoic volcano-sedimentary (quartz-sericite schist and biotite schist) and granitic rocks that have been metamorphosed to lower greenschist facies and variably deformed. The deposit has previously been proposed to be either a porphyry-type or a sediment-hosted stratiform Cu deposit, and its age of formation has been debated.The quartz-sericite schist is interpreted to be a felsic crystal tuff and consists of angular quartz crystals in a fine-grained sericite-altered matrix. Two quartz-sericite schist samples yielded zircon U-Pb upper concordia intercept ages of 2512 ± 12 (2〈span〉σ〈/span〉, mean square of weighted deviates [MSWD] = 0.19) and 2335 ± 16 Ma (2〈span〉σ〈/span〉, MSWD = 0.80). Biotite schist, which is interleaved locally with the quartz-sericite schist and is interpreted to be a basaltic-andesitic sill, yielded a younger zircon U-Pb upper concordia intercept age of 2191 ± 10 Ma (2〈span〉σ〈/span〉, MSWD = 1.7). Five samples of granodiorite and granodiorite porphyry that intruded the schist sequence yielded similar zircon U-Pb ages, with a weighted mean upper concordia intercept age of 2182 ± 7 Ma (2〈span〉σ〈/span〉, MSWD = 1.3). These results suggest that the volcanic sequence was deposited between ∼2.5 and 2.3 Ga and was intruded by basaltic-andesitic sills and a suite of granodiorite and granodiorite porphyry intrusions at ∼2.19 to 2.18 Ga.Two stages of copper mineralization are interpreted to have formed after pervasive sericite alteration of the felsic volcanic rocks. Stage 1 mineralization includes disseminated and deformed quartz veinlets containing chalcopyrite ± pyrite ± magnetite ± molybdenite associated with biotite ± K-feldspar alteration in granodiorite porphyry and schist. Stage 2 comprises undeformed quartz-chlorite-carbonate veins with bornite ± chalcopyrite ± magnetite associated with local chlorite and silicic alteration. Allanite crystals intergrown with chalcopyrite in the granodiorite porphyry yielded an approximate concordia U-Pb age of 2115 ± 31 Ma (2〈span〉σ〈/span〉, MSWD = 2.3). Two molybdenite samples in a deformed quartz-chalcopyrite-molybdenite vein yielded Re-Os model ages of 2106 ± 9 and 2089 ± 9 Ma (2〈span〉σ〈/span〉), consistent with previously published results.Hydrothermal monazite grains with Cu-Fe sulfide inclusions in the granodiorite porphyry, quartz-sericite schist, and undeformed chlorite-bornite-quartz veins yielded much younger U-Pb upper concordia intercept ages of 1832 ± 16 (2〈span〉σ〈/span〉, MSWD = 0.48), 1810 ± 14 (2〈span〉σ〈/span〉, MSWD = 0.92), and 1809 ± 12 Ma (2〈span〉σ〈/span〉, MSWD = 0.38), respectively. The results are in agreement with four Re-Os model ages for pyrite mineral separates from undeformed quartz-sulfide veins, which yielded a weighted mean age of 1807 ± 4 Ma (2〈span〉σ〈/span〉, 〈span〉n〈/span〉 = 4, MSWD = 0.42). In contrast, hydrothermal rutile crystals in the quartz-sericite schist and biotite schist yielded a range of roughly concordant ages between 2.1 and 1.8 Ga, reflecting isotopic disturbance.We interpret these results to indicate original copper mineralization at ∼2.1 Ga that is significantly later than the granodiorite (∼2.18 Ga) and schists (∼2.5–2.2 Ga), followed by hydrothermal remobilization and metamorphism at ∼1.8 Ga. The metavolcanic and granodiorite porphyry host rocks, alteration styles, and disseminated and veinlet form of the earlier mineralization are strongly reminiscent of porphyry Cu deposits, and ages of ∼2.1 Ga have been reported for one intrusion and three volcanic rock samples from the district. The Tongkuangyu, therefore, represents one of the oldest known porphyry copper deposits. Remobilization of copper occurred at ∼1.8 Ga during the Zhongtiao orogeny.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉This study integrates new geologic mapping and structural analysis with previous work near Walnut Canyon and Telegraph Canyon to address the style and magnitude of shortening and the relationship between contractional structures and porphyry preservation and localization between the Ray and Resolution porphyry copper deposits. Cenozoic extensional structures were superimposed on earlier contractional structures formed during the Laramide orogeny, which dates from ~80 to 50 Ma. This superposition requires that Cenozoic normal faults be restored prior to analysis of Laramide contractional structures and their relationship to nearby porphyry copper deposits. Five distinct sets of normal faults within the study area progressively tilted the region 65° east. The amount of extension was 10.3 km or 276%. Using key constraints such as offset strata, cutoff angles between faults and various units, and Laramide fault geometries, the study area was structurally reconstructed and verified using 2-D kinematic modeling of reverse fault offset and related folding. Total shortening is 7.2 km or 98%. Laramide reverse faults are interpreted as thick-skinned basement-cored uplifts, because they restore to moderate angles, have related fault-propagation folds, and involve significant crystalline basement rock. The Telegraph Canyon reverse fault has at least 5.3 km of offset, and the Walnut Canyon reverse fault has 3.2 km. The preferred estimate of the total vertical uplift for the fault system is 5.2 km but could be several kilometers greater. The restored strike direction of these faults combined with mid-Cenozoic erosion surfaces throughout the region suggests that this fault system may be responsible for the Laramide uplift of the Tortilla Mountains and Black Hills. In addition, most major porphyry centers appear to have been intruded into the footwall of this large uplift, with local examples including Ray and Resolution, suggesting that topography generated from this uplift may have been critical to preservation of these ore systems. Though definitive crosscutting relationships do not exist in the immediate map area, geologic relationships in a broader area suggest that shortening here began after 74 Ma and, in the Ray area, had ended by ~69 Ma and that porphyry formation postdated reverse faulting by as much as 5 m.y. to as little as 〈1 m.y.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Mt. Carlton Au-Ag-Cu deposit, northern Bowen basin, northeastern Australia, is an uncommon example of a sublacustrine hydrothermal system containing economic high-sulfidation epithermal mineralization. The deposit formed in the early Permian and comprises vein- and hydrothermal breccia-hosted Au-Cu mineralization within a massive rhyodacite porphyry (V2 open pit) and stratabound Ag-barite mineralization within volcano-lacustrine sedimentary rocks (A39 open pit). These orebodies are all associated with extensive advanced argillic alteration of the volcanic host rocks. Stable isotope data for disseminated alunite (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 6.3–29.2‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉OSO〈sub〉4〈/sub〉 = –0.1 to 9.8‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉OH〈/sub〉 = –15.3 to –3.4‰; 〈span〉δ〈/span〉D = –102 to –79‰) and pyrite (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = –8.8 to –2.7‰), and void-filling anhydrite (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 17.2–19.2‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉SO〈sub〉4〈/sub〉〈/sub〉 = 1.8–5.7‰), suggest that early advanced argillic alteration formed within a magmatic-hydrothermal system. The ascending magmatic vapor (〈span〉δ〈/span〉〈sup〉34〈/sup〉S〈sub〉ΣS〈/sub〉 ≈ –1.3‰) was absorbed by meteoric water (~50–60% meteoric component), producing an acidic (pH ≈ 1) condensate that formed a silicic → quartz-alunite → quartz-dickite-kaolinite zoned alteration halo with increasing distance from feeder structures. The oxygen and hydrogen isotope compositions of alunite-forming fluids at Mt. Carlton are lighter than those documented at similar deposits elsewhere, probably due to the high paleolatitude (~S60°) of northeastern Australia in the early Permian. Veins of coarse-grained, banded plumose alunite (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 0.4– 7.0‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉SO〈sub〉4〈/sub〉〈/sub〉 = 2.3–6.0‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉OH〈/sub〉 = –10.3 to –2.9‰; 〈span〉δ〈/span〉D = –106 to –93‰) formed within feeder structures during the final stages of advanced argillic alteration. Epithermal mineralization was deposited subsequently, initially as fracture- and fissure-filling, Au-Cu–rich assemblages within feeder structures at depth. As the mineralizing fluids discharged into lakes, they produced syngenetic Ag-barite ore. Isotope data for ore-related sulfides and sulfosalts (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = –15.0 to –3.0‰) and barite (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 22.3–23.8‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉SO〈sub〉4〈/sub〉〈/sub〉 = –0.2 to 1.3‰), and microthermometric data for primary fluid inclusions in barite (Th = 116°– 233°C; 0.0–1.7 wt % NaCl), are consistent with metal deposition at temperatures of ~200 ± 40°C (for Au-Cu mineralization in V2 pit) and ~150 ± 30°C (Ag mineralization in A39 pit) from a low-salinity, sulfur- and metal-rich magmatic-hydrothermal liquid that mixed with vapor-heated meteoric water. The mineralizing fluids initially had a high-sulfidation state, producing enargite-dominated ore with associated silicification of the early-altered wall rock. With time, the fluids evolved to an intermediate-sulfidation state, depositing sphalerite- and tennantite-dominated ore mineral assemblages. Void-filling massive dickite (〈span〉δ〈/span〉〈sup〉18〈/sup〉O = –1.1 to 2.1‰; 〈span〉δ〈/span〉D = –121 to –103‰) with pyrite was deposited from an increasingly diluted magmatic-hydrothermal liquid (≥70% meteoric component) exsolved from a progressively degassed magma. Gypsum (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 11.4–19.2‰; 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉SO〈sub〉4〈/sub〉〈/sub〉 = 0.5–3.4‰) occurs in veins within postmineralization faults and fracture networks, likely derived from early anhydrite that was dissolved by circulating meteoric water during extensional deformation. This process may explain the apparent scarcity of hypogene anhydrite in lithocaps elsewhere. While the Mt. Carlton system is similar to those that form subaerial high-sulfidation epithermal deposits, it also shares several key characteristics with magmatic-hydrothermal systems that form base and precious metal mineralization in shallow-submarine volcanic arc and back-arc settings. The lacustrine paleosurface features documented at Mt. Carlton may be useful as exploration indicators for concealed epithermal mineralization in similar extensional terranes elsewhere.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉New SHRIMP U-Pb data from dioritic to granodioritic synmineral intrusions associated with the Jebel Ohier porphyry copper deposit (mineral inventory, including NI43-101-compliant total inferred and indicated resources, of 593 million tonnes [Mt] at 0.33% Cu and 0.05 ppm Au, for 1.953 Mt of contained Cu and 933,600 oz of Au at 0.15% Cu cutoff) in the Red Sea Hills of northeastern Sudan have bracketed the age of mineralization to ca. 816 to 812 Ma. This age range, as well as constraints from new and existing lithogeochemical data, is consistent with the deposit’s formation from a productive parental magma source during the early stages in the evolution of an intra-Mozambique Ocean island arc. The Jebel Ohier porphyry copper deposit bears many similarities to well-documented Phanerozoic analogues elsewhere in terms of (1) the mapped style and zonation of hydrothermal alteration (i.e., proximal K-silicate–dominated, to sericitic, to distal propylitic alteration), (2) the occurrence of intense Cu-bearing A- and B-type vein stockwork, as well as sulfide-only C-type veins, anhydrite veins, and younger, peripheral D-type veins, and (3) the geochemical fingerprint of the associated porphyry, which is akin to those of ore-related Tertiary porphyries in the Escondida area in northern Chile. The multiphase intrusion hosting the Jebel Ohier porphyry copper deposit has been intruded by several generations of mafic to felsic postmineralization dikes and voluminous plutons, with a peak in magmatic activity coinciding with the suturing of the Gebeit terrane at ca. 724 Ma. In spite of, or perhaps because of, the occurrence of extensive postmineralization magmatism, and regardless of subsequent deformation, regional metamorphism, uplift, and erosion, the deposit has remained remarkably intact. The discovery of a relatively ancient, yet well-preserved porphyry copper deposit in the Neoproterozoic Arabian-Nubian Shield has major implications for the exploration potential of this resource-rich geologic terrain.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The newly discovered Shuangjianzishan Ag-Pb-Zn deposit, with 145 Mt of ore grading 128.5 g/t Ag (locally up to 32,000 g/t) and 2.2 wt % Pb + Zn, is located in the Great Hinggan Range metallogenic belt, northeastern China, and is currently the largest Ag deposit in Asia. The Ag-Pb-Zn orebodies occur as veins and are hosted primarily by a Permian slate. Recent drilling and core logging have identified a partially Mo mineralized granite porphyry intrusion adjacent to the Ag-Pb-Zn mineralized veins. This well-preserved magmatic-hydrothermal system therefore offers an excellent opportunity to evaluate the possible temporal and genetic relationship between Mo-mineralized porphyry intrusions and Ag-Pb-Zn veins. Three primary paragenetic stages of veining have been recognized: (I) early pyrite + quartz ± K-feldspar, (II) main ore sulfide + sulfosalt + quartz + calcite + sericite + chlorite ± epidote, and (III) post-ore quartz. The silver mineralization occurs mainly in the late paragenetic part of Stage II, in which canfieldite (Ag〈sub〉8〈/sub〉SnS〈sub〉6〈/sub〉), argentite (Ag〈sub〉2〈/sub〉S) and freibergite [(Ag, Cu)〈sub〉12〈/sub〉Sb〈sub〉4〈/sub〉S〈sub〉13〈/sub〉] are the dominant Ag-bearing ore minerals. A combination of ore mineral chemical and sulfur isotope geothermometers and physicochemical calculations suggest that the Ag-Pb-Zn mineralization took place at a temperature of 250° to 200°C, a pH of 6.7 to 5.6, and a Δlog〈span〉f〈/span〉o〈sub〉2〈/sub〉 (HM) of –2.4 to –8.7.A conspicuous enrichment of Sn and Se in the ore, which is represented by minerals containing the metal suite Ag-Pb-Zn-(Cu-Sn-Se-Sb), likely reflects a close genetic association between the base metal mineralization and a magma. In situ analyses show that the 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values of the sulfides and Ag-bearing sulfosalts from the Ag-Pb-Zn mineralized veins vary from –4.67 to +2.44‰; the mean value is –2.11 ± 1.49‰ (〈span〉n〈/span〉 = 77). The calculated mean 〈span〉δ〈/span〉〈sup〉34〈/sup〉S〈sub〉H2S〈/sub〉 value of the ore-forming fluid is –1.65 ± 0.83‰, which is indicative of a magmatic sulfur source. In situ Pb isotope analyses of the ore minerals yielded a narrow range of values (〈sup〉206〈/sup〉Pb/〈sup〉204〈/sup〉Pb of 18.243–18.310, 〈sup〉207〈/sup〉Pb/〈sup〉204〈/sup〉Pb of 15.503–15.563 and 〈sup〉208〈/sup〉Pb/〈sup〉204〈/sup〉Pb of 38.053–38.203, 〈span〉n〈/span〉 = 59). Comparisons to corresponding isotopic data for the various rock units in the area and sulfides from nearby ore deposits indicate that there were substantial contributions of Pb and other metals (e.g., Ag and Zn) to the Shuangjianzishan deposit from a Mesozoic granitic source.Diorite-granodiorite dikes and dacite are crosscut by the Ag-Pb-Zn veins, and therefore, predate ore formation. These rock units have zircon U-Pb ages of 250.2 ± 2.0 and 133.9 ± 1.4 Ma, respectively. A concealed, weakly Mo mineralized granite porphyry intrusion proximal to the Ag-Pb-Zn mineralized vein system yielded zircon U-Pb ages of 134.4 ± 1.0 (MSWD = 0.1) and 134.4 ± 1.0 Ma (MSWD = 0.2), for coarse- and fine-grained facies, respectively. These ages are indistinguishable within the uncertainty from the zircon ages for the dacite and a granite intrusion ~2 km north of the mineralized veins, which has a weighted mean zircon U-Pb age of 135.2 ± 1.4 Ma (MSWD = 0.78). Molybdenite from three quartz vein/veinlet samples hosted by slate immediately above the porphyry intrusion yielded Re-Os model ages from 136.3 ± 0.9 to 133.7 ± 1.2 Ma and a weighted mean Re-Os age of 134.9 ± 3.4 Ma. Finally, three pyrite samples separated from the Ag-Pb-Zn mineralized veins have a weighted mean Re-Os model age of 135.0 ± 0.6 Ma. The very similar zircon U-Pb ages for the Mo-mineralized granite porphyry and dacite, and Re-Os ages for molybdenite and pyrite in the Shuangjianzishan ore district indicate that the Mesozoic magmatic-hydrothermal activity was restricted to a relatively short time interval (~136–133 Ma). They also suggest that the weakly Mo mineralized granite porphyry was likely the source of the fluids and metals that produced the Ag-Pb-Zn mineralization.Based on our geological observations and an extensive analytical database, a model is proposed for the genesis of the giant Shuangjianzishan Ag-Pb-Zn deposit in which the ore-forming fluid and its metals (i.e., Ag, Pb, and Zn) were exsolved during crystallization of the final phase of a composite granite porphyry intrusion. This fluid transported metals to the distal parts of the system, where they were deposited in preexisting faults or fractures created by the withdrawal of magma during the waning stages of the magmatic-hydrothermal event. The present study of the Shuangjianzishan Ag-Pb-Zn deposit and those of other magmatic-hydrothermal ore deposits in the region provide compelling evidence that the widespread Mesozoic felsic magmatism and Ag-Pb-Zn mineralization in the southern Great Hinggan Range took place in an intracontinental extensional tectonic setting, which was synchronous with, and spatially associated to, Paleo-Pacific slab rollback and lithospheric delamination and thinning.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Orogenic Au deposits have contributed the majority of Au recovered globally throughout history. However, the mechanism that concentrates Au to extremely high bonanza grades in small domains within these deposits remains enigmatic. The volume of fluid required to provide extreme Au endowments in localized occurrences is not reflected in field observations (e.g., in the extent of quartz veining or hydrothermal alteration). Detailed optical, scanning and transmission electron microscopy, nanoscale secondary ion mass spectrometry, and 3-D neutron tomography have been used to investigate the processes responsible for development of anomalously high grade ore (upward of 3% Au) found in quartz veins at Fosterville gold mine (Victoria, Australia). Distinct textural settings of visible Au include (1) Au concentrated along pressure solution seams associated with wall-rock selvages, (2) as nano- to microscale dusty Au seams parallel to pressure solution seams, and (3) in microscale tension fractures perpendicular to stylolitic seams. The distribution of Au in arsenopyrite and pyrite hosted within pressure solution seams changes as a function of the extent of deformation. Sulfides in highly deformed pressure solution seams exclusively host Au as nano- to micrometer-sized clusters within features associated with corrosion and brittle failure, whereas sulfides in mildly deformed pressure solution seams have Au bound in the crystal structure. It is proposed that Au supersaturation in fluids introduced during seismic periods led to the deposition of abundant Au nanoparticles in quartz-carbonate veins. Subsequent pressure dissolution of vein quartz and carbonate during interseismic intervals allowed for episodic increase in the Au/quartz ratio and permitted liberation and migration of Au nanoparticles, promoting Au grain growth in favorable textural settings. Galvanic corrosion and brittle fracturing of auriferous sulfides during the interseismic period allowed additional remobilization and/or enrichment of sulfide-hosted Au. Repetition of this mechanism over the time scale of deposit formation acted to concentrate Au within the lodes. This Au ore upgrading model, referred to as “aseismic refinement,” provides a new insight for the genesis of ultrarich Au mineralization and, based on textures reported from many Au deposits, may be a globally significant component in the formation of orogenic Au deposits.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉To better understand processes leading to porphyry Mo deposit formation, the metal content, volatile content, and crystallization conditions of melt inclusions from pre- and synmineralization intrusions in six porphyry(-skarn) Mo deposits of northeastern China (Aolunhua, Hashitu, Lanjiagou, Songbei, Wanbaoyuan, and Yangjiazhangzi) were investigated by means of laser ablation-inductively coupled plasma-mass spectrometry and electron microprobe analysis. The ore-forming silicate melts were one to four times more evolved than average granite with 1 to 7 ppm Mo. The ore-related intrusions crystallized predominantly at 760° to 690°C and 3.7 to 1.0 kbar, except for the one at Hashitu, which crystallized at 770° to 740°C and lower pressures (2.0–1.0 kbar). Fertile silicate melts at Hashitu contain up to 0.4 wt % F, 0.03 to 0.09 wt % Cl, 5.0 to 7.0 wt % H〈sub〉2〈/sub〉O, 10 to 24 ppm Cs, and 200 to 500 ppm Rb, whereas those at Yangjiazhangzi and Wanbaoyuan contain less Cs (3–6 ppm and 5–7 ppm, respectively), less Rb (180–220 ppm and 200–240 ppm, respectively), and negligible F (〈0.15 wt %) but have similar Cl (0.03–0.05 wt %) and H〈sub〉2〈/sub〉O (5.3–6.5 wt % and 4.0–5.2 wt %, respectively) contents. Calculated melt viscosities in fertile magmas (log 〈span〉η〈/span〉 = 4.3–6.1 Pa s) are at the lower end of the values reported for felsic melts at the same temperature.Comparison between syn- and premineralization intrusions in individual deposits reveals that the ore-related intrusions were similarly evolved and had similar Mo contents and crystallization conditions as the nonmineralizing intrusions. The only difference is that the premineralization intrusions tend to occur as batholiths. The key to porphyry Mo mineralization lies in the focusing of fluid into and through a small rock volume on the top of the intrusion.For the studied porphyry Mo deposits, the mineralizing magmas are all Mo poor, indicating Mo enrichment is not required to form porphyry Mo deposits. Metal endowments in porphyry Mo deposits have no direct relationship with the composition and crystallization condition of mineralizing melts but are linked with the fluid flux released from the underlying magma chamber through a cupola.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉High-grade hematite mineralization is widely developed in banded iron formations (BIFs) worldwide. However, in the North China craton where Neoarchean-Paleoproterozoic BIFs are abundant, economic high-grade hematite ores are scarce. High-grade hematite ores hosted in the Paleoproterozoic Yuanjiacun BIFs represent the largest occurrence of this type of ore in the North China craton. The orebodies are fault controlled and show sharp contacts with lower greenschist facies metamorphic BIFs. In situ U-Pb geochronology of monazite and xenotime intergrown with microplaty hematite and martite in high-grade ore established two episodes of metamorphic-hydrothermal monazite/xenotime growth after deposition of the BIFs. The earlier episode at ca. 1.94 Ga is interpreted as the timing of lower greenschist-facies metamorphism, and the later episode at 1.41 to 1.34 Ga represents the timing of high-grade hematite mineralization. Petrography and microthermometry of primary fluid inclusion assemblages indicate that the high-grade hematite ore formed from hot (313°–370°C), CO〈sub〉2〈/sub〉-rich, and highly saline (~20 wt % NaCl equiv) hydrothermal fluids. These fluids channeled along faults, which concentrated iron through interaction with the BIFs—a process similar to typical hematite mineralization elsewhere. The deposition of hematite was probably related to tectonic extension in the North China craton related to the breakup of the Columbia/Nuna supercontinent. Our results challenge a previously proposed model ascribing the scarcity of high-grade hematite ores in the North China craton to the lack of prolonged weathering conditions. Rather, we argue that the high-grade ore formed in lower metamorphic-grade BIFs at shallower depths than magnetite mineralization and was largely eroded during later exhumation and uplift of the craton.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Ore assemblages in uranium roll-front deposits are highly variable and heavily dependent on Eh/pH conditions. Sulfur isotopes in pyrite traditionally have been employed to distinguish between biogenic or abiogenic redox pathways as drivers of roll-front propagation. However, the extent of and constraints on bacterial productivity have never been quantified, nor have the chemical conditions imposed by either primary formation mechanism. Moreover, this approach implicitly assumes that deposits form via one process or the other and disregards the possibility that both processes participate simultaneously in generating some orebodies. In this study, we analyzed sulfur isotopes from pyrite coprecipitated with uranium in two Wyoming roll-front deposits: Lost Creek and Willow Creek Mine Unit 10. The results document contrasting isotopic fractionation that correlates with pyrite morphology. Both deposits evolved with both abiogenic and biogenic redox mechanisms as active contributors to ore formation. In the past, bimodal fractionation behavior with pyrite morphology has been attributed to distinct temporal episodes of pyrite formation, driven by either a change in redox mechanism or multiple independent fluid events with unique isotopic signatures. However, neither explanation is appropriate for the isotopic trends identified in this study, where the two pyrite morphologies appear coeval in both deposits. Moreover, the contemporaneous formation of both pyrite morphologies cannot occur under the same conditions by the same precipitation mechanism because of the difference in their free energies of formation. The data suggest a third alternative in which pyrite morphology correlates to its biogenic or abiogenic mode of formation. Given the isotopic composition of pre-ore pyrite, sulfur isotope fractionation trends within the ore zone can be applied to establish prolificacy of bacteria and chemical conditions of the ore-forming solution.In both study sites, framboidal pyrite occurred as the primary by-product of sulfur-reducing bacteria, and the corresponding fractionation pattern constrains the sulfur availability and bacterial productivity. Euhedral to anhedral pyrite precipitated from abiogenic redox, the sulfur fractionation recording Eh/pH gradients during ore evolution. At Lost Creek, framboidal pyrite produced 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values from –50.8 ± 0.5‰ to +142.8 ± 0.3‰, while subhedral pyrite ranged from –68.1 ± 0.4‰ to +33.8 ± 0.3‰ with 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values increasing toward the barren, unaltered contact. Pre-ore pyrite at Lost Creek ranged from –0.8 ± 0.5‰ to +70.6 ± 0.3‰. 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values from biogenically derived pyrite at Lost Creek indicate a closed system with limited sulfate availability and a slow rate of bacterial reduction, implying restricted bacterial activity. Abiogenic fractionation behavior indicates a system driven by an Eh drop under neutral or basic pH conditions, and pyrite distribution across the roll identifies abiogenic pyrite recycling as the dominant redox mechanism at Lost Creek. At Willow Creek Mine Unit 10, framboidal pyrite ranged from –32.5 ± 0.4‰ to +68.2 ± 0.4‰, and subhedral pyrite ranged from –45.1 ± 0.4‰ to +5.4 ± 0.4‰. The subhedral pyrite 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values initially increased into the center of the roll and subsequently decreased again approaching the barren, unaltered contact. Pre-ore pyrite ranged from –48.1 ± 0.4‰ to +15.6 ± 0.5‰. Willow Creek Mine Unit 10 biogenically produced 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values show minimal fractionation from pre-ore pyrite, indicating an open system with abundant sulfate and rapid reduction from prolific bacterial activity. The abiogenic trends indicate an Eh drop and low pH at the barren, altered contact progressively neutralized across the orebody. This correlates to the anticipated Eh/pH gradients in a system dominated by biogenic redox.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Karangahake deposit is the third largest gold producer in the Hauraki goldfield, Coromandel, New Zealand. Production was mainly from the Maria and Welcome/Crown veins that were mined over a strike length of 〉1,300 m and vertical interval of 700 m. The veins are hosted in andesite that is strongly altered over an area of approximately 4.2 × 2.7 km. By using petrography, qualitative and quantitative X-ray diffraction (XRD), and automated mineralogy from a scanning electron microscope (SEM), hydrothermal alteration has been characterized in detail for surface rocks exposed along the Ohinemuri and Waitawheta rivers and along the underground Keillors crosscut, which links the Maria and Welcome/Crown veins. The chemistry of hydrothermally altered rocks was analyzed using portable X-ray fluorescence (pXRF). The automated mineralogy reveals alteration mineral occurrences, abundances, and some textural relationships that are not as apparent from conventional studies.Hydrothermal alteration of the andesite is typically strong to intense, with 100% conversion of primary to secondary minerals, but decreases to moderate to weak, comprising 〈25% secondary minerals toward the west and east, where the margins of alteration are exposed. The andesite is altered by combinations of quartz, chlorite, adularia, albite, illite, pyrite, calcite, mixed-layered illite-smectite, smectite, and rare epidote. Quartz and chlorite are common and present in both strongly and weakly altered rocks. Adularia is widespread and extends laterally up to 500 m from the Maria and Welcome/Crown veins. Illite is broadly coextensive with and locally occurs beyond the adularia halo. The zone of illite alteration is bordered to the east by patchy mixed-layered illite-smectite and outlying smectite alteration. Along the Keillors crosscut, the rocks between the Maria and Welcome/Crown veins are altered to adularia, albite, illite, chlorite, and rare calcite and epidote.Automated mineralogy and quantitative XRD mineral estimates indicate many rocks have been altered by 〉30% adularia (max. 74%) and that the amount of adularia can exceed that of quartz. Adularia and albite occur throughout the Keillors crosscut, with adularia most abundant within 20 to 50 m of the Maria, Mystery, and Welcome/Crown veins, whereas albite is more abundant between these veins.Hydrothermal alteration has resulted in the variable gains in K and Rb forming anomalous halos that surround veins, which are attributed to adularia and illite. By contrast, Ca, Sr, and to a lesser extent Fe are variably lost around veins in broadly overlapping zones of depletion.The key hydrothermal alteration minerals are illite, illite-smectite, smectite, adularia, and albite. The systemscale zonation of illite, illite-smectite, and smectite outlines the broad thermal structure, with illite stability indicating 〉230°C, illite-smectite indicating 130° to 230°C, and smectite indicating 〈130°C. Adularia indicates zones of inferred high permeability and the central upflow of boiling hydrothermal fluids, whereas albite indicates lower inferred permeability and peripheral-style alteration. Adularia represents a broad-scale vector toward veins, and its occurrence and abundance can be determined using field portable XRD or chemically by proxy from the K/Al value.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Chengchao Fe skarn deposit (280 Mt @ 45% Fe; nearly half of the total Fe reserve has an ore grade higher than 53 wt %) is the largest high-grade magnetite skarn deposit in the Middle-Lower Yangtze River metallogenic belt. The magnetite orebodies and skarns occur mainly along the contact zones between granitic (129 ± 2 Ma) rocks and diorite (133 ± 1 Ma) and Triassic strata. The paragenesis of the Chengchao deposit includes the following five stages: Na-K alteration stage (albite + K-feldspar), prograde skarn stage (garnet + pyroxene), retrograde skarn stage (magnetite + phlogopite + amphibole), sulfate-sulfide stage (anhydrite + pyrite ± garnet ± magnetite), and carbonate stage (calcite).Four generations of garnet and magnetite with distinct textures are recognized. The elemental compositions and fluid inclusion data of the three generations of garnet from the prograde skarn stage indicate fluctuating changes in the oxygen fugacity during their precipitation. Some magnetite grains from the retrograde skarn stage are characterized by well-developed zonation with respect to Al, K, Ca, and Si, with concentrations of up to hundreds of parts per million. The higher contents of (V + Ti) and (Al + Mn) in Mt2 relative to Mt1 (0.8 and 2,321 ppm for Mt1, and 182 and 10,315 ppm for Mt2, respectively) may be indicative of increasing fluid temperatures from episodic magmatic events. Fluid inclusion data show that the fluids responsible for the prograde skarn stage have high temperature (〉750°C), high salinity (〉50 wt % NaCl equiv) and high Fe concentration (e.g., magnetite and pyrite daughter crystals within garnet and pyroxene). Oxygen and hydrogen isotope data suggest that magmatic water was dominant during the prograde skarn stage and that increasing percentages of meteoric water were involved during later stages. Multiple episodic incursions of Fe-rich fluids, which were probably related to the emplacement of granitic rocks and diorite, were crucial for controlling the formation of the high-grade Chengchao Fe skarn deposit.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉High- and low-temperature platinum nuggets from two historical localities, Chocó in Colombia and Córrego Bom Sucesso in Brazil, are compared with respect to trace elements. Supergene platinum nuggets from Córrego Bom Sucesso are enriched in Se and Hg, but depleted in siderophile and chalcophile elements, and have fractionated platinum group element patterns, in comparison with magmatic platinum nuggets from Chocó. In particular, Se concentrations over ~100 〈span〉μ〈/span〉g/g Se and S/Se ratios above unity indicate Se recycling in a supergene environment with abundant organic matter.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The transport and deposition of gold from colloidal suspensions in hydrothermal fluids has been a persistent theme in ore deposits research. Studies of active geothermal systems show that a complete model of gold transport must include both dissolved and particulate forms. However, samples of the hydrothermal fluids are commonly spiked with aqua regia after collection in order to put any solids back into solution, thus preventing a quantitative assessment of the particle load. Although attempts have been made to filter the solids, gold nanoparticles (Au NPs) will mostly pass the 0.2-〈span〉μ〈/span〉m filters that are in common use, and a simple technique for analyzing suspended particles in the liquids has been lacking. In this study, we demonstrate how time-resolved acquisition of mass 197 in a conventional inductively coupled plasma-mass spectrometer (ICP-MS) can be used to detect and measure Au NPs in the filtered liquids, with an example of well-characterized fluids from the Reykjanes geothermal field on Iceland. The technique allows for precise monitoring of the solution as it is introduced into the plasma with the capability of identifying individual particles carried in suspension. Results show that Au particles passing the 0.2-〈span〉μ〈/span〉m filters are abundant in the studied samples, and measurements of the individual particles can be used to determine their size. The experiment highlights the potential of emerging ICP-MS techniques, including very fast data acquisition and multielement analysis of single particles in time-of-flight mode, for characterization of NPs in hydrothermal fluids.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉We report the first Ag isotope data for the Paleozoic orogenic Au deposits in the Victorian Goldfields, southeast Australia, a world-class province with a historic production of 2,400 tonnes of Au. Despite their relatively uniform geology—similar host-rock types, age, mineralization style—deposits in Victoria show a wide range in 〈sup〉107〈/sup〉Ag/〈sup〉109〈/sup〉Ag ratios in native Au (〈span〉ε〈/span〉〈sup〉107〈/sup〉Ag −6.6 to +8.3, relative to the NIST SRM 978a Ag standard), comparable to the entire previously known terrestrial range (−9.4 to +5.3). The data show no correlation with mineralization age or host-rock composition, and there is no obvious isotopic link to established “mantle” or “crustal” Ag isotope values, implying that source rock signatures are unlikely to be the main control on Ag isotope variations. Instead, it is suggested that the Ag isotope variation is primarily related to physicochemical processes, particularly Ag isotope fractionation during redox reactions such as conversion of Ag〈sup〉0〈/sup〉 in native Au to Ag〈sup〉+〈/sup〉 in dissolved Ag(HS)〈sub〉2〈/sub〉〈sup〉−〈/sup〉 or sulfide-borne Ag. Repeated Ag〈sup〉0〈/sup〉〈strong〉←→〈/strong〉 Ag〈sup〉+〈/sup〉 reactions along transport pathways and at sites of ore accumulation could generate a wide range in 〈span〉ε〈/span〉〈sup〉107〈/sup〉Ag, and evidence of this range is presented in the data here. Silver isotope fractionation via numerous deposition-dissolution cycles provides a different perspective into large-scale ore genesis that has not previously been recognized for orogenic gold systems; multistaged metal remobilization along fluid transport pathways is standard during their formation. Detailed Ag isotope studies have considerable potential for understanding the history of episodic metal addition and within-deposit redistribution.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Four porphyry Cu-Mo systems were investigated by Re-Os molybdenite geochronology to constrain their timing with respect to the geodynamic and magmatic evolution of the eastern Pontides, Turkey. Molybdenite from the Ispir-Ulutaş deposit yielded an Re-Os age of 131.0 ± 0.7 Ma, which is consistent with Early Cretaceous U-Pb laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) zircon ages of local calc-alkaline intrusions. It demonstrates that porphyry deposits were already formed during Early Cretaceous subduction of the Neotethys along the eastern Pontides, and that they can be correlated with porphyry Cu events in the adjacent Lesser Caucasus. Molybdenite Re-Os ages of 76.0 ± 0.4 and 75.7 ± 0.4 Ma at the Elbeyli prospect and 77.2 ± 1.0 Ma at the Emeksen prospect overlap with U-Pb LA-ICP-MS zircon ages of shoshonitic to high-K calc-alkaline intrusions in the region, which were emplaced during Late Cretaceous Neotethys subduction. A 50.7 ± 0.3 Ma molybdenite Re-Os age at the Güzelyayla deposit confirms porphyry Cu-Mo emplacement coeval with Eocene postcollisional, calc-alkaline adakitic magmatism of the eastern Pontides.An electron microprobe study of molybdenite samples, supplemented by data obtained during Re-Os dating, shows that the Eocene Güzelyayla deposit and the Late Cretaceous Emeksen prospect have the highest Re enrichment. Postcollisional melting of a thickened mafic lower continental crust and melting of a metasomatized lithospheric mantle with little to no interaction with upper crustal rocks may explain the Re enrichment at Güzelyayla and Emeksen, respectively.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The trace element composition of zircon, especially in tandem with U-Pb geochronology, has become a powerful tool for tracing magmatic processes associated with the formation of porphyry copper deposits. However, the use of the redox-sensitive Eu and Ce anomalies as a potential mineral exploration proxy is controversial. This study presents a comprehensive, temporally constrained data set of zircon trace element compositions (n = 645) for three compositionally distinct magmatic series identified in the Meghri-Ordubad pluton, southernmost Lesser Caucasus. The 30 million years of Cenozoic magmatism in the Meghri-Ordubad pluton are associated with several ore-forming pulses leading to the formation of porphyry copper deposits and epithermal-style mineralization. Our zircon geochemical data constrain the thermal and chemical evolution of this complex intrusive suite and allow an evaluation of the usefulness of zircon as a mineral exploration proxy for porphyry copper deposits. Our results combined with Rayleigh fractionation modeling indicate that the trace element composition of zircon (Th/U, Hf, Ti, Yb〈sub〉N〈/sub〉/Dy〈sub〉N〈/sub〉, Eu anomalies) is influenced by the composition and the water concentration of the parental magma, as well as by co-crystallizing titanite and apatite. In contrast, the variations of Ce anomalies remain difficult to explain by magmatic processes and could rather be ascribed to relative fluctuations of the redox conditions. In the Meghri-Ordubad pluton, we do not observe any systematic patterns between the trace element composition in zircons and the different ore-forming pulses. This questions the reliability of using the trace element composition in zircon as an exploration mineral proxy, and it rather emphasizes that a good knowledge of the entire magmatic evolution of a metallogenic province is required.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Muratdere Cu-Mo (Au) porphyry deposit in western Turkey contains elevated levels of rhenium and is hosted within granodioritic intrusions into an ophiolitic mélange sequence in the Anatolian belt. The deposit contains several stages of mineralization: early microfracture-hosted molybdenite and chalcopyrite, followed by a quartz-pyrite-chalcopyrite vein set associated with Cu-Au grade, a quartz-chalcopyrite-pyrite-molybdenite vein set associated with Cu-Mo-Re grade, and a later polymetallic quartz-barite-sphalerite-galena-pyrite vein set. The rhenium in Muratdere is hosted within two generations of molybdenite: early microfracture-hosted molybdenite and later vein-hosted molybdenite. In situ laser ablation-inductively coupled plasma-mass spectrometry analysis of sulfides shows that the later molybdenite has significantly higher concentrations of Re (average 1,124 ppm, 〈span〉σ〈/span〉 = 730 ppm, n = 43) than the early microfracture-hosted molybdenite (average 566 ppm, 〈span〉σ〈/span〉 = 423 ppm, n = 28). Pyrite crystals associated with the Re-rich molybdenite have higher Co and As concentrations than those in other vein sets, with Au associated with As. The microfracture-hosted sulfides have 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values between −2.2‰ and +4.6‰, consistent with a magmatic source. The vein-hosted sulfides associated with the high-Re molybdenite have a 〈span〉δ〈/span〉〈sup〉34〈/sup〉S signature of 5.6‰ to 8.8‰, similar to values found in peridotite lenses in the Anatolian belt. The later enrichment in Re and 〈span〉δ〈/span〉〈sup〉34〈/sup〉S-enriched S may be sourced from the surrounding ophiolitic country rock or may be the result of changing redox conditions during deposit formation.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Kapan mining district in the southernmost Lesser Caucasus is one of the few locations along the central Tethyan metallogenic belt where ore-forming processes were associated with magmatic arc growth during Jurassic Tethys subduction along the Eurasian margin. Three ore deposits of the Kapan district were investigated in this study: Centralni West, Centralni East, and Shahumyan. The ore deposits are hosted by Middle Jurassic andesitic to dacitic volcanic and volcaniclastic rocks of tholeiitic to transitional affinities below a late Oxfordian unconformity, which is covered by calc-alkaline to transitional Late Jurassic-Early Cretaceous volcanic rocks interlayered with sedimentary rocks.The mineralization consists of veins, subsidiary stockwork, and partial matrix replacement of breccia host rocks, with chalcopyrite, pyrite, tennantite-tetrahedrite, sphalerite, and galena as the main ore minerals. Centralni West is a dominantly Cu deposit, and its host rocks are altered to chlorite, carbonate, epidote, and sericite. At Centralni East, Au is associated with Cu, and the Shahumyan deposit is enriched in Pb and Zn as well as precious metals. Both deposits contain high-sulfidation mineral assemblages with enargite and luzonite. Dickite, sericite, and diaspore prevail in altered host rocks in the Centralni East deposit. At the Shahumyan deposit, phyllic to argillic alteration with sericite, quartz, pyrite, and dickite is dominant with polymetallic veins, and advanced argillic alteration with quartz-alunite ± kaolinite and dickite is locally developed.The lead isotope composition of sulfides and alunite (〈sup〉206〈/sup〉Pb/〈sup〉204〈/sup〉Pb = 18.17–18.32, 〈sup〉207〈/sup〉Pb/〈sup〉204〈/sup〉Pb = 15.57–15.61, 〈sup〉208〈/sup〉Pb/〈sup〉204〈/sup〉Pb = 38.17–38.41) indicates a common metal source for the three deposits and suggests that metals were derived from magmatic fluids that were exsolved upon crystallization of Middle Jurassic intrusive rocks or leached from Middle Jurassic country rocks. The 〈span〉δ〈/span〉〈sup〉18〈/sup〉O values of hydrothermal quartz (8.3–16.4‰) and the 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values of sulfides (2.0–6.5‰) reveal a dominantly magmatic source at all three deposits. Combined oxygen, carbon, and strontium isotope compositions of hydrothermal calcite (〈span〉δ〈/span〉〈sup〉18〈/sup〉O = 7.7–15.4‰, 〈span〉δ〈/span〉〈sup〉13〈/sup〉C = −3.4−+0.7‰, 〈sup〉87〈/sup〉Sr/〈sup〉86〈/sup〉Sr = 0.70537–0.70586) support mixing of magmatic-derived fluids with seawater during the last stages of ore formation at Shahumyan and Centralni West.〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar dating of hydrothermal muscovite at Centralni West and of magmatic-hydrothermal alunite at Shahumyan yield, respectively, a robust plateau age of 161.78 ± 0.79 Ma and a disturbed plateau age of 156.14 ± 0.79 Ma. Re-Os dating of pyrite from the Centralni East deposit yields an isochron age of 144.7 ± 4.2 Ma and a weighted average age of the model dates of 146.2 ± 3.4 Ma, which are younger than the age of the immediate host rocks. Two different models are offered, depending on the reliability attributed to the disturbed 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar alunite age and the young Re-Os age. The preferred interpretation is that the Centralni West Cu deposit is a volcanogenic massive sulfide deposit and the Shahumyan and Centralni East deposits are parts of porphyryepithermal systems, with the three deposits being broadly coeval or formed within a short time interval in a nascent magmatic arc setting, before the late Oxfordian. Alternatively, but less likely, the three deposits could represent different mineralization styles successively emplaced during evolution and growth of a magmatic arc during a longer time frame between the Middle and Late Jurassic.〈/span〉

Description: 〈span〉〈span〉Economic Geology〈/span〉 thanks the following individuals for book reviews that they prepared for the 2018 volume of the journal. Their reviews appear on the page indicated by the number that follows each name.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Igarapé Cinzento (GT-46) iron oxide copper-gold deposit is located in the northwestern part of the Cinzento shear zone in Carajás province, Brazil. Copper mineralization in the deposit is hosted by a hydrothermally altered metavolcanosedimentary sequence that consists of interlayered amphibolites, (actinolite)-biotite schists, almandine-biotite schists, and metamorphosed banded iron formations, intruded by granitoids (tonalite to granite) and mafic dikes.The complex evolution of the deposit involved (1) deposition of a volcanosedimentary sequence (2774 ± 19 Ma U-Pb zircon), (2) a tectonometamorphic event (D〈sub〉n〈/sub〉 = upper amphibolite facies), restricted in time between 2.77 and 2.72 Ga, (3) development of the Cinzento shear zone (D〈sub〉n + 1〈/sub〉) at 2.72 Ga, with associated hydrothermal activity and deposit formation, (4) superimposed hydrothermal events (ca. 2.6 and 2.5 Ga), and (5) granitoid emplacement at ca 2.5 Ga in the region.The GT-46 deposit is the result of overprinting hydrothermal stages with at least two distinct mineralization events. Variations in alteration styles are a consequence of changes in the mechanisms of fluid control, first regulated by ductile deformation and later by brittle fracturing. Early hydrothermal alteration is characterized by zones of sodic-calcic and potassic alteration as well as iron metasomatism with associated disseminated mineralization I, controlled by ductile structures (chalcopyrite-bornite-magnetite, 2718 ± 56 Ma, Re-Os molybdenite). Veins and breccias with chalcopyrite-magnetite(-chlorite-calcite-quartz) represent an overprinting episode of mineralization (II) at ca. 2.6 Ga. Late alteration stages (ca. 2.5 Ga) represent fluid remobilization with chloritization zones and crystal growth spatially related to pegmatite intrusions. Re-Os molybdenite ages of 2503 ± 51 and 2449 ± 44 Ma provide evidence for younger episodes of hydrothermal fluid circulation in the deposit.Fluid inclusion study reveals aqueous-saturated, highly saline (32−〉50 wt % NaCl equiv, total homogenization temperature [T〈sub〉ht〈/sub〉] = 203°−〉450°C) and carbonic fluids (CO〈sub〉2〈/sub〉 melting temperature [T〈sub〉mco〈sub〉2〈/sub〉〈/sub〉] = –56.3° to –54.5°C) for early alteration zones and saline (2.1–27.5 wt % NaCl equiv, T〈sub〉ht〈/sub〉 = 153.3°–215.5°C) solutions for late alteration zones. These data provide evidence for several alternative fluids involved in the evolution of the deposit: (1) hypersaline fluids derived from hydrothermal-magmatic brines or evaporative brines/evaporite-derived fluids, (2) saline solutions generated by progressive dilution of primary hypersaline fluids with shallower, colder, and less saline fluids (e.g., meteoric, seawater), and (3) carbonic components sourced from magmatic fluids, devolatilization of carbonate units and, partially, metamorphic fluids.Sulfur is interpreted to have been leached from country rocks in the earlier hydrothermal stages (mineralization I, 〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 0.01–1.99‰), and to have come from mixed granitic and external oxidized sources for later overprinting events (mineralization II, 〈span〉δ〈/span〉〈sup〉34〈/sup〉S = 9.75–11.25‰). These data provide evidence for the recurrence of hydrothermal activity, indicating the existence of different potential sources of metals and sulfur during the evolution of the deposit.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The world-renowned Pliocene El Laco iron deposit in northern Chile is the youngest and the best-preserved Kiruna-type deposit in the world. The genesis of the El Laco Kiruna-type iron deposit—i.e., whether it is magmatic or hydrothermal in origin—is a long-standing controversy. The interstitial Fe-P phase lined by early formed magnetite crystals observed in massive ores at El Laco is morphologically and geochemically consistent with that produced by recent immiscible experiments, confirming that the massive magnetite ores were products of complete solidification of an iron-rich mush comprising early crystallized magnetite and an interstitial immiscible Fe-P melt. The hydrothermal features in geochemistry of massive magnetite are similar to the features of magnetite in altered andesites, implying a superimposed hydrothermal process. The occurrence of melt inclusions with high homogenization temperatures (〉700°C) hosted by the apatite in Cristales Grandes ores indicates that the veined ores formed shortly after the massive ores genetically related to the magmatic system. Some vesicles in the massive ores and the magnetite scoriae, previously interpreted as compelling evidence of volcanic structures, are demonstrated to be the residual pore spaces formed after interstitial phases between magnetite grains were removed by postmagmatic hydrothermal fluids and the accumulated debris eroded from massive ores with minor altered andesites, respectively. Thus, we propose a two-step genetic model in which massive ores resulted from magmatic processes and then were modified by subsequent magmatic-hydrothermal alteration. This model could be applied to the metallogenesis of most Kiruna-type iron deposits.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Reconstructing the thermal evolution of hydrothermal ore deposits mainly relies on their fluid inclusion record, which is limited by favorable trapping conditions, calling for alternative temperature constraints. Muscovite and tourmaline coexist in many hydrothermal ore deposits and in the granites or pegmatites related with them. Whereas in situ analyses of boron isotopes in tourmaline are widely applied to constrain fluid sources and evolution, muscovite has seldom been used in this way, and the potential for isotope exchange thermometry with this mineral pair is unexplored. The different boron coordination in muscovite and tourmaline causes a significant, temperature-dependent isotopic fractionation between them, which has been determined experimentally. We used this relationship to study mineralization conditions and fluid evolution at the Panasqueira W-Sn-Cu deposit in Portugal, where the source and evolution of the mineralizing fluids are still debated. The difference in 〈sup〉11〈/sup〉B/〈sup〉10〈/sup〉B ratios of coexisting muscovite and tourmaline, expressed as Δ〈sup〉11〈/sup〉B〈sub〉mica-tourmaline〈/sub〉 = 〈span〉δ〈/span〉〈sup〉11〈/sup〉B〈sub〉mica –〈/sub〉〈span〉δ〈/span〉〈sup〉11〈/sup〉B〈sub〉tourmaline〈/sub〉, yields median temperatures for vein selvages from 400° to 460°C within a total range of 350° to 600°C, which agrees with published Ti-in-quartz temperatures. Mineral pairs from a late fault zone yield a lower median temperature of about 250°C (range 220°–320°C), which fits with published homogenization temperatures of quartz-hosted fluid inclusions from the veins. Taking these temperatures into account, the calculated fluid composition of the early and late muscovite generations is about 〈span〉δ〈/span〉〈sup〉11〈/sup〉B〈sub〉fluid〈/sub〉 = −6 ± 2‰, which indicates that the recurrent fluid pulses had a uniform composition and a magmatic-hydrothermal origin.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Faults are important structures in the formation of many mineral deposits, often acting as conduits for ore-forming fluids and sometimes providing, or generating, the bounding structures to associated mineralizing sites. Using 3-D analysis and modeling of the Lisheen and Silvermines deposits within the Irish ore field, we investigate the geometry of normal fault systems and their implications on the origin and nature of associated deposits. These Irish-type deposits are carbonate hosted and developed within the hanging walls of normal faults arising from an Early Carboniferous episode of north-south rifting, with relatively limited amounts of later deformation. Structural analysis of high-quality mine datasets indicates that fault segmentation is ubiquitous with left-stepping segments arising from north-south stretching developed above generally ENE-NE-trending fault arrays, which are subparallel to older Caledonian penetrative fabrics and structure within underlying Silurian and Ordovician rocks. Fault segments occur on different scales and have a profound impact on structural evolution, with larger scale segments and intervening relay ramps defining distinct orebodies within deposits and smaller scale segments and relays potentially providing paths for upfault fluid flow. The difference in behavior is attributed to the integrity of associated relay ramps where intact ramps represent orebody-bounding structures, and smaller breached ramps provide enhanced associated hydraulic properties and act as vertical conduits. Hanging-wall deformation along the rheological boundary between host-rock limestones and underlying shales has an important control on the localization of earlier dolomitization and/or brecciation and later mineralization adjacent to this contact, and on the migration pathways for basinal brines and mineralizing fluids.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Active sea floor hydrothermal systems are modern analogues of ancient volcanogenic massive sulfide deposits. Hydrothermal sulfide chimneys are one of the important components in these systems and are formed by rapid mixing between seawater and metal-rich hydrothermal fluids venting onto the sea floor. Previous models of chimney growth have been built up mainly based on studies including optical petrographic observations and scanning electron microscopy (SEM). The present study, conducted on a sample from the PACMANUS hydrothermal field (Manus basin, Papua New Guinea), reports for the first time the electron backscattered diffraction (EBSD) observations on a modern sea floor chimney, coupled with synchrotron X-ray fluorescence (SXRF) and scanning electron microscopy-backscattered electron (SEM-BSE) imagery to reveal fine-scale and primary microstructures of sphalerite, allowing the reconstruction of crystal growth history. The results show that sphalerite clusters are formed via the coalescence of multiple smaller sphalerite globules. Electron backscattered diffraction images also highlight that each globule includes an inner zone with fine-grained particles (〈1 〈span〉μ〈/span〉m) and an outer zone with elongate blade-shaped crystals (length up to 40 〈span〉μ〈/span〉m), in some cases showing branching dendritic habit. Both zones are dominated by sphalerite with minor other sulfides, such as chalcopyrite, pyrite, and wurtzite. The individual globules are interpreted as forming under conditions of supersaturation within high-temperature gradients, for example, during the mixing between high-temperature (e.g., 300°C) hydrothermal fluids and ambient cold seawater. The occurrence of these contrasting inner and outer zones reflects fluctuation between two regimes: bursts of rapid nucleation from supersaturated fluids that occurred almost instantaneously during the initial mixing and the skeletal crystal growth of particular crystal faces from limited nucleation sites where diffusion-limited boundary layers developed around the growing crystals. Those growing globules then coalesced into columnar aggregates. These observations could potentially have important implications for identifying fossil chimneys in ancient ore deposits. Moreover, this study emphasizes the benefits of advanced techniques, such as EBSD and SXRF, in order to characterize sulfides in various hydrothermal chimneys to reveal crystal growth and fluid mixing history and gain further insight into chimney growth processes.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Combined whole-rock lead isotope and high-precision trace element data are reported for shales from the Wolverine deposit and surrounding prospects (e.g., Puck, Fisher, Sable) to elucidate the source of Pb in the shales and their potential utility as an exploration vector. Shales exhibit distinct variations in time-integrated U/Th/Pb evolution with proximity to mineralization. Distal samples resemble the local crustal Pb isotope composition of the northern Cordilleran continental crust (i.e., shale curve). More proximal samples exhibit evidence for hydrothermal activity, but were not significantly mineralized. They have very high U/Th and U/Pb and anomalous present-day 〈sup〉206〈/sup〉Pb/〈sup〉204〈/sup〉Pb (〉20) coincident with elevated Y/Ho, Zr/Hf, and carbonate alteration. These signatures were derived from oxygenated seawater and imprinted on the shales by oxygenated, CO〈sub〉2〈/sub〉-rich hydrothermal fluids that contained uranyl carbonate ions (e.g., UO2(CO3)34−). It is envisioned that this was due to shallow-level (within 50 m of the seafloor), near-vent ingress of oxygenated seawater associated with hydrothermal venting. A third population of samples is most proximal to mineralization. These, too, have trace element signatures (i.e., high Y/Ho, Zr/Hf, CO〈sub〉2〈/sub〉-U enrichment) inherited from seawater, but they do not exhibit anomalous present-day 〈sup〉206〈/sup〉Pb/〈sup〉204〈/sup〉Pb, being more juvenile in 〈sup〉206〈/sup〉Pb/〈sup〉204〈/sup〉Pb, 〈sup〉207〈/sup〉Pb/〈sup〉204〈/sup〉Pb, and 〈sup〉208〈/sup〉Pb/〈sup〉204〈/sup〉Pb. This isotope composition coincides with enrichments in Pb, Zn, Cu, Tl, Sb, W, and Sn (±Cu, Mo). The Pb enrichment associated with mineralization effectively buffered these samples from appreciable Pb isotope evolution, despite overprinting by uranyl carbonate. The juvenile Pb associated with the shales proximal to mineralization argues for deep-penetrative circulation of hydrothermal fluids and leaching of Pb from mafic igneous rocks at depth and/or influence from juvenile magmas at depth (i.e., magmatic fluids). Both of these processes were likely critical for the genesis of mineralization at Wolverine.This study illustrates the utility of quadrupole inductively coupled plasma-mass spectrometry (ICP-MS) for providing sufficiently precise and accurate Pb isotope determinations for district-scale metallogenic evaluations. Furthermore, the coupling of solution ICP-MS Pb isotopes with high-precision trace element geochemistry from the same digest provides a powerful tool for regional exploration studies for hydrothermal mineralization in shale basins.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Carbonate-hosted Au deposits (12 million tonnes at 5.0 g/t Au) are located less than 1.5 km away from Cu-Au skarn deposits in the Fengshan area, Middle-Lower Yangtze River metallogenic belt, eastern China. The skarn deposits occur in the contact zones of late Mesozoic granodiorite porphyry intrusions. Bismuth tellurides and coloradoite (HgTe) occur in Cu-Au skarn and peripheral Au-As ores, respectively, whereas tellurides, including coloradoite, and lorandite (TlAsS〈sub〉2〈/sub〉) are present in carbonate-hosted Au ores outside the marble front. New zircon U-Pb ages indicate that the granodiorite porphyry intrusions spatially related to both Cu-Au skarn and carbonate-hosted Au deposits were emplaced between 149.9 ± 1.3 and 145.9 ± 0.7 Ma, which overlaps with molybdenite Re-Os ages of 149.1 ± 2.1 to 145.4 ± 2.9 Ma from the Cu-Au skarn deposits. Sulfur isotope compositions of sulfide minerals from the Cu-Au skarn deposits (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = –2.5 to 6.4〈span〉‰〈/span〉) are similar to those from carbonate-hosted Au deposits (〈span〉δ〈/span〉〈sup〉34〈/sup〉S = –4.8 to 2.6〈span〉‰〈/span〉), suggesting a similar magmatic source of sulfur. Geologic, mineralogical, and sulfur isotope evidence collectively indicate that the carbonate-hosted Au deposits are the distal product of Cu-Au skarn mineralization rather than Carlin-type Au deposits, as has been previously suggested. This study proposes a genetic model to illustrate possible linkages between oxidized Cu-Au skarn deposits (associated with oxidized magmas) and carbonate-hosted Au deposits, which may have important implications for Au deposits in the more distal parts of Cu-Au skarn systems in the Middle-Lower Yangtze River metallogenic belt and elsewhere.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉High-resolution trace element mapping (2 〈span〉µ〈/span〉m beam) was performed by synchrotron radiation X-ray fluorescence microscopy (SXFM) on Cu- and Zn-rich chimneys from Brothers volcano, Kermadec arc. The maps cover 84- to 136-mm〈sup〉2〈/sup〉 cross sections of the inner chimney walls and document the distribution of Fe, Cu, Zn, As, Se, Sr, Pb ± Ga, Au, Bi, and U. Comparative element maps were generated by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) at a lower resolution (47 〈span〉µ〈/span〉m beam), which additionally measured Co, Ni, Mo, Ag, Cd, In, Sn, Sb, Ba, Au, Tl ± V, Te, and Bi. In two Zn-rich chimneys, Cu distribution varies from Cu-bearing sphalerite with chalcopyrite disease to distinct chalcopyrite-lined channels, implying a progression of chalcopyrite replacement of sphalerite. Conversely, the two Cu-rich chimneys have different styles of massive chalcopyrite lining their conduits. The first displays elongate chalcopyrite grains that radiate into and infill the conduit; these radial grains may be related to decreased fluid flow during sealing of the chimney at both its top and bottom. The second style involves multiple, concentric chalcopyrite laminations (0.25–1 mm) inside the conduit, which progressively narrowed the orifice as they were deposited. Fine (15–40 〈span〉µ〈/span〉m) trace element rings are revealed within and between laminations that exhibit variable contents and distributions of Co, Ni, Zn, As, Se, Mo, Ag, Cd, Sn, Te, Au, Tl, Pb, Bi, and U. The presence of U in the rings specifically indicates seawater ingress into the chimney interior despite it vigorously discharging 274°C fluids at the time of sampling. During these periodic seawater incursions, rapidly changing chemical gradients within the chimney wall induce the instantaneous precipitation of metals from the vent fluid. Thus, the trace element rings are a proxy for the secular evolution of vent fluid compositions. We compared enrichment factors of trace element rings to those of fumarole condensates studied at subaerial arc and rift volcanoes and molten S pooled atop a submarine back-arc volcano. Our enrichment factors show remarkable consistency with the other volcanoes and indicate Au, Te, Bi, Cu, Ag, and Se in Brothers’ chimneys were magmatically derived.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Although magmatic-hydrothermal ore deposits may contain as much as 1 billion tonnes (Gt) of sulfur as hydrothermal sulfides and sulfates, whether the genetically related magmas are particularly sulfur rich remains an open question. In oxidized melts, the presence of magmatic anhydrite would provide some evidence of this, but its preservation potential in the rock record is low, owing to the high solubility of sulfate in low-temperature aqueous fluids. In this study we have examined granitic and porphyry samples from porphyry copper ore deposits using quantitative evaluation of materials by scanning electron microscopy (QEMSCAN) and scanning electron microscopy-energy dispersive spectrometry (SEM-EDS) to identify small and sparse inclusions of anhydrite preserved in a range of igneous minerals including apatite, amphibole, plagioclase, quartz, and titanite. In total, anhydrite inclusions were present in 11 different samples from six sulfur-rich magmatic-hydrothermal ore deposits. Electron microprobe analysis (EMPA) of anhydrite inclusions suggests that primary magmatic anhydrite typically has elevated Sr (〉2,500 ppm) and P (〉500 ppm) compared to hydrothermal anhydrite; however, in many cases this is obscured by postentrapment exchange between the inclusion and host mineral or hydrothermal alteration of host minerals. Where the composition of inclusions is inconclusive, we have used cathodoluminescence imagery of zoning patterns in host minerals to help distinguish between igneous and hydrothermal origins. This study provides evidence that anhydrite saturation may be a common feature of arc magmas, particularly those associated with ore deposits, even though little evidence of magmatic anhydrite is preserved. Widespread anhydrite saturation has significant implications for any consideration of magmatic sulfur budgets, particularly in low-temperature, evolved magmas, which have little capacity to carry sulfur species in the melt.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Irish-type mineralization is commonly attributed to fault-controlled mixing of a seawater-derived, sulfur-rich fluid and basement-derived, metal-rich fluid. However, maar-diatreme volcanoes discovered in close spatial and temporal association with Zn-Pb mineralization at Stonepark in the Limerick basin (southwest Ireland) bring a new dimension to established geologic models and may increase the deposit-scale prospectivity in one of the world’s greatest Zn-Pb districts. Stonepark exhibits many incidences of dolomitic black matrix breccias with associated Zn-Pb mineralization, the latter typically occurring within 150 m of the diatremes. Highly negative 〈span〉δ〈/span〉〈sup〉34〈/sup〉S pyrite values within country rock-dominated black matrix breccias (–12 to –34〈span〉‰〈/span〉) are consistent with sulfide precipitation from bacteriogenic sulfur reduction in seawater-derived brines. However, 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values of Zn-Pb sulfides replacing black matrix breccias (–10 to 1〈span〉‰〈/span〉) reflect multiple sulfur sources. Diatreme emplacement both greatly enhanced country rock fracture permeability and produced conduits that are filled with porous volcaniclastic material and extend down to basement rock types. Our 〈span〉δ〈/span〉〈sup〉34〈/sup〉S data suggest that diatremes provide more efficient fluid pathways for basement-derived fluids. The diatremes introduce another potential sulfur source and facilitate a greater input of metal-rich basement-derived hydrothermal fluid into the system compared to other Irish-type deposits such as Navan and Lisheen, evidenced by Stonepark’s more positive modal 〈span〉δ〈/span〉〈sup〉34〈/sup〉S value of –4〈span〉‰〈/span〉. Irish-type deposits are traditionally thought to form in association with extensional basement faults and are considered unrelated to extensive Carboniferous magmatism. Our results indicate that a direct link exists between diatreme volcanism and Zn-Pb mineralization at Limerick, prompting a reevaluation of the traditional Irish-type ore formation model, in regions where mineralization is spatially associated with volcanic pipes.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Textural and compositional data of magnetite from Igarapé Bahia, Alemao, Sossego, Salobo, and Candelaria iron oxide copper-gold (IOCG) and El Romeral Kiruna-type iron oxide-apatite (IOA) deposits show that some magnetite grains display oscillatory zoning or have been reequilibrated by oxy-exsolution, coupled dissolution and reprecipitation (CDR) reactions, and/or recrystallization. Textures formed via CDR are most widespread in the studied samples. The original oscillatory zoning was likely derived from the crystal growth during fluctuating fluid compositions rather than from variation in temperature and oxygen fugacity. The oxy-exsolution of ilmenite in magnetite is attributed to increasing oxygen fugacity and decreasing temperature with alteration and mineralization, resulting in product magnetite with lower Ti and higher V contents. Recrystallization of some magnetite grains is commonly due to high-temperature annealing that retained primary compositions. Two different types of CDR processes are defined according to textures and chemical compositions of different generations of magnetite. The first generation of magnetite (Mag-1) is an inclusion-rich and trace element-rich core, which was replaced by an inclusion-poor and trace element-poor rim (Mag-2). The third generation of magnetite (Mag-3), inclusion poor but trace element rich, occurs as veins replacing Mag-2 along fractures or grain margins. Type 1 CDR process transforming Mag-1 to Mag-2 is more extensive and is similar to processes reported in skarn deposits, whereas type 2 CDR process is local, transforming Mag-2 to Mag-3. During type 1 CDR process, minor and trace elements Si, K, Ca, Mg, Al, and Mn in magnetite are excluded, and Fe contents increase to various extents, in contrast to type 2 CDR process, which is characterized by increased contents of Si, K, Ca, Mg, Al, and Mn. Type 1 CDR process is possibly induced by the changing fluid composition and/or decreasing temperature during progressive alteration and ore formation, whereas type 2 CDR process can be interpreted as post-ore replacement due to a new pulse of magmatic-hydrothermal fluids.The identification of magnetite core (Mag-1) with igneous origin and rim (Mag-2) with magmatic-hydrothermal origin in the Sossego IOCG and El Romeral IOA deposits supports a fluid changing from magmatic to magmatic-hydrothermal during IOCG and IOA formation and indicates a genetic link between these two deposit types. The large data set here further demonstrates that magnetite is susceptible to textural and compositional reequilibration during high-temperature magmatic and magmatic-hydrothermal processes. Reequilibrated magnetite, particularly that formed by CDR processes, has a chemical composition that can be different from that of primary magnetite. Modified magnetite, therefore, cannot be used to discriminate its primary origin or to interpret its provenance in overburden sediments. Therefore, in situ chemical analysis of magnetite combined with textural characterization is necessary to understand the origin of magnetite in IOCG and IOA deposits.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Tinakula is the first seafloor massive sulfide deposit described in the Jean Charcot troughs and is the first such deposit described in the Solomon Islands—on land or the seabed. The deposit is hosted by mafic (basaltic-andesitic) volcaniclastic rocks within a series of cinder cones along a single eruptive fissure. Extensive mapping and sampling by remotely operated vehicle, together with shallow drilling, provide insights into deposit geology and especially hydrothermal processes operating in the shallow subsurface. On the seafloor, mostly inactive chimneys and mounds cover an area of ~77,000 m〈sup〉2〈/sup〉 and are partially buried by volcaniclastic sand. Mineralization is characterized by abundant barite- and sulfide-rich chimneys that formed by low-temperature (〈250°C) venting over ~5,600 years. Barite-rich samples have high SiO〈sub〉2〈/sub〉, Pb, and Hg contents; the sulfide chimneys are dominated by low-Fe sphalerite and are high in Cd, Ge, Sb, and Ag. Few high-temperature chimneys, including zoned chalcopyrite-sphalerite samples and rare massive chalcopyrite, are rich in As, Mo, In, and Au (up to 9.26 ppm), locally as visible gold. Below the seafloor, the mineralization includes buried intervals of sulfide-rich talus with disseminated sulfides in volcaniclastic rocks consisting mainly of lapillistone with minor tuffaceous beds and autobreccias. The volcaniclastic rocks are intensely altered and variably cemented by anhydrite with crosscutting sulfate (± minor sulfide) veins. Fluid inclusions in anhydrite and sphalerite from the footwall (to 19.3 m below seafloor; m b.s.f.) have trapping temperatures of up to 298°C with salinities close to, but slightly higher than, that of seawater (2.8–4.5 wt % NaCl equiv). These temperatures are 10° to 20°C lower than the minimum temperature of boiling at this depth (1,070–1,204 m below sea level; m b.s.l.), suggesting that the highest-temperature fluids boiled below the seafloor. The alteration is distributed in broadly conformable zones, expressed in order of increasing depth and temperature as (1) montmorillonite/nontronite, (2) nontronite + corrensite, (3) illite/smectite + pyrite, (4) illite/smectite + chamosite, and (5) chamosite + corrensite. Zones of argillic alteration are distinguished from chloritic alteration by large positive mass changes in K〈sub〉2〈/sub〉O (enriched in illite/smectite), MgO (enriched in chamosite and corrensite), and Fe〈sub〉2〈/sub〉O〈sub〉3〈/sub〉 (enriched in pyrite associated with illite/smectite alteration). The 〈span〉δ〈/span〉〈sup〉18〈/sup〉O and 〈span〉δ〈/span〉D values of clay minerals confirm increasing temperature with depth, from 124° to 256°C, and interaction with seawater-dominated hydrothermal fluids at high water/rock ratios. Leaching of the volcanic host rocks and thermochemical reduction of seawater sulfate are the primary sources of sulfur, with 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values of sulfides, from –0.8 to 3.4‰, and those of sulfate minerals close to seawater sulfate, from 19.3 to 22.5‰.The mineralization and alteration at Tinakula are typical of a class of ancient massive sulfide deposits hosted mainly by permeable volcaniclastic rocks with broad, semiconformable alteration zones. Processes by which these deposits form have never been documented in modern seafloor massive sulfide systems, because they mostly develop below the seafloor. Our study shows how hydrothermal fluids can become focused within permeable rocks by progressive, low-temperature fluid circulation, leading to a large area (〉150,000 m〈sup〉2〈/sup〉) of alteration with reduced permeability close to the seafloor. In our model, overpressuring and fracturing of the sulfate- and clay-cemented volcaniclastic rocks produced the pathways for higher-temperature fluids to reach the seafloor, present now as sulfate-sulfide veins within the footwall. In the geologic record, the sulfate (anhydrite) is not preserved, leaving a broad zone of intense alteration with disseminated and stringer sulfides typical of this class of deposits.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The trace element composition of igneous and hydrothermal magnetite from 19 well-studied porphyry Cu ± Au ± Mo, Mo, and W-Mo deposits was measured by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and then classified by partial least squares-discriminant analysis (PLS-DA) to constrain the factors explaining the relationships between the chemical composition of magnetite and the magmatic affinity and porphyry deposit subtypes. Igneous magnetite can be discriminated by relatively high P, Ti, V, Mn, Zr, Nb, Hf, and Ta contents but low Mg, Si, Co, Ni, Ge, Sb, W, and Pb contents, in contrast to hydrothermal magnetite. Compositional differences between igneous and hydrothermal magnetite are mainly controlled by the temperature, oxygen fugacity, cocrystallized sulfides, and element solubility/mobility that significantly affect the partition coefficients between magnetite and melt/fluids. Binary diagrams based on Ti, V, and Cr contents are not enough to discriminate igneous and hydrothermal magnetite in porphyry deposits.Relatively high Si and Al contents discriminate porphyry W-Mo hydrothermal magnetite, probably reflecting the control by high-Si, highly differentiated, granitic intrusions for this deposit type. Relatively high Mg, Mn, Zr, Nb, Sn, and Hf but low Ti and V contents discriminate porphyry Au-Cu hydrothermal magnetite, most likely resulting from a combination of mafic to intermediate intrusion composition, high chlorine in fluids, relatively high oxygen fugacity, and low-temperature conditions. Igneous or hydrothermal magnetite from Cu-Mo, Cu-Au, and Cu-Mo-Au deposits cannot be discriminated from each other, probably due to similar intermediate to felsic intrusion composition, melt/fluid composition, and conditions such as temperature and oxygen fugacity for the formation of these deposits.The magmatic affinity of porphyritic intrusions exerts some control on the chemical composition of igneous and hydrothermal magnetite in porphyry systems. Igneous and hydrothermal magnetite related to alkaline magma is relatively rich in Mg, Mn, Co, Mo, Sn, and high field strength elements (HFSEs), perhaps due to high concentrations of chlorine and fluorine in magma and exsolved fluids, whereas those related to calc-alkaline magma are relatively rich in Ca but depleted in HFSEs, consistent with the high Ca but low HFSE magma composition. Igneous and hydrothermal magnetite related to high-K calc-alkaline magma is relatively rich in Al, Ti, Sc, and Ta, due to a higher temperature of formation or enrichment of these elements in melt/fluids.Partial least squares-discriminant analysis on hydrothermal magnetite compositions from porphyry Cu, iron oxide copper-gold (IOCG), Kiruna-type iron oxide-apatite (IOA), and skarn deposits around the world identify important discriminant elements for these deposit types. Magnetite from porphyry Cu deposits is characterized by relatively high Ti, V, Zn, and Al contents, whereas that from IOCG deposits can be discriminated from other types of magnetite by its relatively high V, Ni, Ti, and Al contents. IOA magnetite is discriminated by higher V, Ti, and Mg but lower Al contents, whereas skarn magnetite can be separated from magnetite from other deposit types by higher Mn, Mg, Ca, and Zn contents. Decreased Ti and V contents in hydrothermal magnetite from porphyry Cu and IOA, to IOCG, and to skarn deposits may be related to decreasing temperature and increasing oxygen fugacity. The relative depletion of Al in IOA magnetite is due to its low magnetite-silicate melt partition coefficient, immobility of Al in fluids, and earlier, higher-temperature magmatic or magmatic-hydrothermal formation of IOA deposits. The relative enrichment of Ni in IOCG magnetite reflects more mafic magmatic composition and less competition with sulfide, whereas elevated Mn, Mg, Ca, and Zn in skarn magnetite results from enrichment of these elements in fluids via more intensive fluid-carbonate rock interaction.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Haigou lode gold deposit (〉40 tons [t] at 3.4 g/t), which is located near the eastern boundary of the Central Asian orogenic belt and the North China craton, is one of the largest gold deposits in northeastern China. Native gold is intergrown with molybdenite and pyrite in auriferous quartz veins hosted by a monzogranite-monzonite stock and locally by Proterozoic gneiss, thereby offering an excellent opportunity to directly date the mineralizing event. Uranium-Pb age determinations for zircon yielded ages for the monzogranite and monzonite of 327.1 ± 1.1 and 329.5 ± 1.0 Ma, respectively. Numerous mafic to felsic dikes, which are crosscut by ore veins (pre-ore), parallel to these veins (possibly synore), or crosscut by them (post-ore), were carefully examined and dated. Their zircon 〈sup〉206〈/sup〉Pb/〈sup〉238〈/sup〉U ages are 318.3 ± 1.0, 310.9 ± 1.1, and 134.9 ± 0.4 Ma, respectively, thereby placing the timing of gold mineralization within the relatively large interval of 318.3 ± 1.0 to 134.9 ± 0.4 Ma. The age of mineralization was determined directly using the Re-Os method applied to molybdenite. A total of 19 molybdenite samples separated from auriferous quartz veins yielded widely differing Re-Os model ages of 467 to 155 Ma, and replicate analyses of individual samples also yielded widely differing ages. Significantly, the wide range is attributable entirely to the results obtained for some coarse-grained molybdenite samples and is interpreted to be due to Re and Os isotope decoupling, the considerable spatial Re heterogeneity, the analytical procedure (e.g., use of small sample aliquots), and the post-ore deformation. Nine of the samples, which are all fine grained, yielded a robust weighted mean model age of 310 ± 3 Ma and an isochron age of 309 ± 8 Ma. Thus, the molybdenite Re-Os ages are identical, within uncertainty, to those of the dikes that are parallel to the ore veins, indicating that these dikes were emplaced contemporaneously with the ore and that they and the Haigou gold mineralization are of late Paleozoic age (ca. 310 Ma). Finally, a sericite sample obtained from an auriferous vein returned a 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar plateau age of 165.3 ± 1.2 Ma, which is much younger than the age of the mineralization constrained by Re-Os age determinations of molybdenite. This indicates that the 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar isotope system was reset by post-ore thermal events.Our new geochronological data provide evidence for late Paleozoic gold mineralization in Haigou, which makes it the oldest known lode gold deposit in the easternmost Central Asian orogenic belt, a finding that has important implications for precious metal mineral exploration in the eastern part of the Solonker-Xar Moron-Changchun-Yanji suture zone between the Central Asian orogenic belt and the North China craton. This study also indicates that accurate and reproducible molybdenite Re-Os ages representing the true timing of ore deposition need an integrated combination of careful petrography, proper sampling procedures, sufficiently large analyzed aliquots, multiple analyses of individual samples, and multiple dating methods.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Lone Tree deposit is located in the northern Battle Mountain mining district, Nevada. Prior to mine closure in 2006, Santa Fe Pacific Gold and Newmont produced 4.2 Moz of gold at an average grade of 2.06 g/t at Lone Tree, primarily from the N-S– to NNW-SSE–striking Wayne zone. The ore is located between the Roberts Mountain and Golconda thrusts in siliciclastic rocks of the Ordovician Valmy Formation and in the Pennsylvanian-Permian Battle Mountain and Edna Mountain Formations, and above the Golconda thrust in siliciclastic and carbonate rocks of the Mississippian to Permian Havallah sequence. Ore is also hosted by rhyolitic dikes that were emplaced at 40.95 ± 0.06 Ma based on zircon U-Pb chemical abrasion-thermal ionization mass spectrometry.The gold is associated with sericitic and argillic alteration of the siliciclastic rocks and dikes and with decarbonatization and Fe carbonate alteration of the carbonate-bearing units, as well as in Fe-As sulfide and finegrained quartz alteration of all rock types. Oxidation affects 30 to 45% of the deposit, penetrating into the stratigraphy along numerous steeply dipping north-south, east-west, and north-northeast–south-southwest structures. Gold is positively correlated with Ag, As, Hg, and Sb. The highest Au grades occur in quartz-sulfide ore hosted in siliciclastic and carbonate sedimentary rocks and rhyolitic intrusions. In this ore style, fine-grained quartz and sericite are intergrown with disseminated sulfide minerals (quartz-sericite-pyrite alteration), constituting cores of weakly mineralized pyrite or marcasite, which are surrounded by fuzzy arsenopyrite rims that contain up to ~2,000 ppm Au. Low gold grades occur in late-stage banded pyrite breccias consisting of a finely zoned Au-poor pyrite matrix surrounding jigsaw-fit clasts of quartz-, illite-, barite-, and adularia-altered siliciclastic rock. The timing of main-stage mineralization is bracketed between the emplacement of the dikes and an adularia 〈sup〉40〈/sup〉Ar/〈sup〉39〈/sup〉Ar age of 40.14 ± 0.74 Ma.Sericite intergrown with arsenopyrite-rimmed pyrite in phenocrysts of the rhyolite dikes gave 〈span〉δ〈/span〉〈sup〉18〈/sup〉O values of 1.6 to 9.5‰ and 〈span〉δ〈/span〉D values of –105 to –145‰. For temperatures of 300 ± 100°C, the calculated fluid isotopic compositions are consistent with felsic magmatic water and minor modifications by mixing with meteoric water and exchange with wall rocks. In the silica-sulfide ore, in situ isotopic laser ablation-multicollector-inductively coupled plasma-mass spectrometry (LA-MC-ICP-MS) analyses of pyrite cores yielded 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values ranging from 3.4 to 7.7‰, with average values of 5.6‰ in the felsic dikes, 4.5‰ in the siliciclastic rocks, and 5.3‰ in the carbonate rocks. These values match conventional pyrite 〈span〉δ〈/span〉〈sup〉34〈/sup〉S data reported for Eocene porphyry systems elsewhere in the district. Nanoscale secondary ion mass spectrometry analyses show that gold and associated trace elements occur in submicron-scale zones within arsenopyrite rims on pyrite. The average 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values of the arsenopyrite rims are 5.3 to 6.5‰ heavier than the pyrite cores, indicating cooling and an increasing H〈sub〉2〈/sub〉S/SO〈sub〉2〈/sub〉 ratio. The highest grades resulted from episodic pulses of a gold-rich fluid that was partly derived from, or exchanged with, the sedimentary host rocks. In situ LA-MC-ICP-MS 〈span〉δ〈/span〉〈sup〉34〈/sup〉S values for the late-stage banded pyrite breccia become progressively lighter from veinlet margin to center, reaching a low of –32‰. These veinlets indicate a shift from main-stage quartz-sericite-pyrite and intermediate argillic alteration to more neutral pH and oxidizing conditions during late-stage mineralization, indicating either increasing interaction between the fluid and sedimentary sulfur sources in the host-rock package or bacterial sulfate reduction and supergene sulfide precipitation.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Nanling region of South China hosts the largest W-Sn metallogenic province in the world, accounting for more than 54% of global tungsten resources as well as important resources of tin and rare metals. An important feature of this province, which is shared by a number of other W-Sn metallogenic provinces, is that W deposits occur separately from Sn and Sn-W deposits, with the latter concentrated in the western part of the region (especially along the deep, NE-trending Chenzhou-Linwu fault) and the W deposits to the east of them. All the deposits are associated with ilmenite series, peraluminous granites. However, the granites associated with the Sn and Sn-W deposits can be distinguished from the W granites by their higher bulk-rock 〈span〉ε〈/span〉Nd values and their higher zircon 〈span〉ε〈/span〉Hf values. Most importantly, the Sn and Sn-W granites are characterized by higher zircon saturation temperatures (800 ± 20°C) than the W granites (650°–750°C). The Sn and Sn-W granites also contain abundant mantle-derived mafic microgranular enclaves, whereas such enclaves are rare in the W granites. A model is proposed in which the protolith to the W granites released W to the melt as a result of the breakdown of muscovite. The temperature of melting, however, was too low for biotite to melt. In the west, particularly along the Chenzhou-Linwu fault (the location of the Sn and Sn-W deposits), higher temperatures enabled the breakdown of both muscovite and biotite and the consequent release of both Sn and W to form Sn and Sn-W granites. This model, which is based on differences in the protolith melting temperature and thus mobilization temperatures for Sn and W, is potentially applicable to any Sn-W metallogenic province in which the Sn and Sn-W deposits are spatially separated from the W deposits.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The use of indium in modern technologies has grown in recent decades, creating a growth in indium demand; thus, there is a need to constrain the spatial and temporal distribution of indium-bearing, granite-related deposits. Toward this end, a conceptual model and exploration vectors for the formation of granite-related indium deposits have been developed. The magmatic-hydrothermal system is modeled by consideration of crystal-melt and vapor-melt equilibria. The model calculates the efficiency of removal of indium from a melt into a volatile phase, which may serve as a component of an ore-forming fluid. The results of the model suggest that as the proportion of ferromagnesian minerals increases in the associated granites, the probability of indium ore formation decreases. Further, for a given modal proportion of ferromagnesian minerals, as the modal proportion of amphibole increases, the probability of indium ore formation decreases. Lastly, for a given modal proportion of biotite, as the magnesium content of the biotite increases (as would result from increasing oxidation of the magmatic system), the probability of indium ore formation decreases. Granites with the highest probability of being associated with indium ore formation will typically be part of A- or S-type igneous systems and will likely be highly fractionated (e.g., A-type topaz granites). I-type granites will generally have a lower potential of being associated with indium-bearing deposits. However, some I-type granites may be associated with indium-bearing deposits if the deposits contain granites (sensu stricto) or other related rocks (e.g., alaskites) that lack amphibole or other ferromagnesian phases.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Molybdenite-bearing porphyry deposits are the predominant supplier of molybdenum to industrialized society and one of the main hosts of Mo in the upper continental crust. The Mo isotope compositions (δ〈sup〉98/95〈/sup〉Mo, normalized to NIST3134 equals 0‰) of molybdenite show considerable variation (–1.62 to +2.27‰), but the factors controlling this variability remain poorly constrained. This information is critical for underpinning genetic models of porphyry deposits, understanding elemental cycling, and utilizing the δ〈sup〉98/95〈/sup〉Mo of marine sediments as a paleoredox proxy. Using the well-characterized Qulong porphyry Cu-Mo deposit (Tibet) as an example, here we discuss how rapid cooling, facilitated by mixing hot magmatic fluid with cold meteoric water, can be a controlling factor on efficient mineralization, and then tackle how fluid evolution regulates molybdenum isotope fractionation. Molybdenites, which preferentially partition isotopically light Mo (Rayleigh fractionation), precipitated from a single fluid will develop a heavier δ〈sup〉98/95〈/sup〉Mo composition over time, and this also creates heterogeneous δ〈sup〉98/95〈/sup〉Mo between molybdenite grains. Whereas a fluid undergoing multiple episodes of intensive boiling will gradually lose its isotopically heavy Mo to the vapor phase, molybdenites crystallizing successively from the residual liquid will then have lighter δ〈sup〉98/95〈/sup〉Mo over time. However, when mineralization efficiency becomes too low, a negligible variation in δ〈sup〉98/95〈/sup〉Mo of molybdenite is observed. Given that the mineralization efficiency (i.e., the amount of Mo crystallized as molybdenite from the fluid) rarely reaches 100% and molybdenite favors isotopically light Mo, the presence of a residual fluid with isotopically heavy Mo is inevitable. This residual fluid may then become trapped in alteration halos; hence, δ〈sup〉98/95〈/sup〉Mo has the potential to aid in locating the mineralization center (e.g., lighter δ〈sup〉98/95〈/sup〉Mo toward the orebody). The residual fluid may also feed surface hydrological systems and eventually impact Mo cycling. Our study highlights that understanding the controls of isotope fractionation is critical to bridge the gap between ore formation and elemental cycling, and that other transition metals (e.g., Cu, Fe, and Zn) may follow similar trajectories.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Zangezur-Ordubad mining district of the southernmost Lesser Caucasus is located in the central segment of the Tethyan metallogenic belt and consists of porphyry Cu-Mo and epithermal Au and base metal systems hosted by the composite Cenozoic Meghri-Ordubad pluton. Ore-hosting structures and magmatic intrusions are predominantly confined to a central N-S–oriented corridor 40 km long and 10 to 12 km wide, located between two regional NNW-oriented right-lateral faults, the Khustup-Giratagh and Salvard-Ordubad faults. The anatomy and kinematics of the main fault network are consistent with dextral strike-slip tectonics controlled by the NNW-oriented Khustup-Giratagh and Salvard-Ordubad faults.Dextral strike-slip tectonics was initiated during the Eocene, concomitantly with final subduction of the Neotethys, and controlled the emplacement of the Agarak, Hanqasar, Aygedzor, and Dastakert porphyry Cu-Mo and Tey-Lichkvaz and Terterasar epithermal Au and base metal deposits. The Eocene structures were repeatedly reactivated during subsequent Neogene evolution in transition to a postsubduction geodynamic setting. Ore-bearing structures at the Oligocene world-class Kadjaran porphyry Cu-Mo deposit were also controlled by dextral strike-slip tectonics, as well as porphyry mineralization and its epithermal overprint hosted by an early Miocene intrusion at Lichk.Eocene to early Miocene dextral strike-slip tectonics took place during NE- to NNE-oriented compression related to Paleogene Eurasia-Arabia convergence and subsequent Neogene postcollision evolution. Paleostress reconstruction indicates major reorganization of tectonic plate kinematics since the early Miocene, resulting in N-S– to NW-oriented compression. Early Miocene epithermal overprint at the Kadjaran porphyry deposit and left-lateral reactivation of faults and mineralized structures are linked to this late Neogene tectonic plate reorganization.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The Au-rich polymetallic massive sulfide orebodies of the Kassandra mining district belong to the intrusion-related carbonate-hosted replacement deposit class. Marble lenses contained within the Stratoni fault zone host the Madem Lakkos and Mavres Petres deposits at the eastern end of the fault system, where paragenetically early skarn and massive sulfide are spatially associated with late Oligocene aplitic and porphyritic dikes. Skarn transitions into predominant massive and banded replacement sulfide bodies, which are overprinted by a younger assemblage of boulangerite-bearing, quartz-rich sulfide and late quartz-rhodochrosite vein breccias. The latter style of mineralization is most abundant at the Piavitsa prospect at the western end of the exposed fault system. The sulfide orebodies at the Olympias deposit are hosted by marble in association with the Kassandra fault, where textural and mineralogical similarities to the sulfide bodies within the Stratoni fault zone suggest a genetic relationship. Estimated trapping temperatures and pressures based on fluid inclusion data indicate that carbonate replacement mineralization took place at depths less than about 5.9 km.Carbon and oxygen isotope patterns in carbonate from the Stratoni fault zone support isotopic exchange principally through fluid–wall-rock interaction, whereas decarbonation and fluid-rock exchange reactions were important at the Olympias deposit. Carbonate minerals associated with skarn and replacement sulfide throughout the district have isotopic compositions that are consistent with formation from a hydrothermal fluid of magmatic origin. Lower homogenization temperatures and salinities in the younger quartz-rich sulfide assemblage and quartz-rhodochrosite vein breccias, together with low 〈span〉δ〈/span〉〈sup〉18〈/sup〉O values of gangue carbonate, suggest dilution of a primary magmatic fluid with meteoric water late in the evolution of the hydrothermal system in both the Olympias area and the Stratoni fault zone. The replacement sulfide orebodies in the district likely inherited their uniform Pb isotope composition from a late Oligocene igneous source and the isotopically heterogeneous metamorphic basement units.Metal distribution patterns at the scale of the Stratoni fault zone show diminishing Cu concentration with decreasing Pb/Zn and Ag/Au ratios from Madem Lakkos to Mavres Petres and the Piavitsa prospect in the west. The sulfide orebodies at the Olympias deposit exhibit elevated Cu values in the east with increasing Pb/Zn and Ag/Au ratios down-plunge to the south-southwest. Metal concentration and ratios support zoning related to temperature and solubility changes with increasing distance from a probable magmatic source. Structural and igneous relationships, together with fluid inclusion microthermometric and carbon-oxygen isotope data and metal distribution patterns, are supportive of a zoned hydrothermal system that exceeded 12 km along the Stratoni fault zone, sourced by an igneous intrusion to the southeast of the Madem Lakkos deposit. The Olympias replacement sulfide orebodies, associated with the Kassandra fault, resulted from a local hydrothermal system that was likely derived from a concealed igneous intrusion to the east of the deposit.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Major Au and Cu deposits in the Western Tethyan magmatic belt formed during two main periods of Cretaceous and Cenozoic magmatism. The Cretaceous deposits are dominantly Cu-Au porphyry, high-sulfidation epithermal, and volcanic massive sulfide deposits, whereas in the Cenozoic Cu is significant only in porphyry systems. However, the Cenozoic contains approximately three times greater total Au endowment (for Au deposits 〉0.5 million ounces), and also has a greater deposit diversity, including porphyry Au-Cu and Au-only deposits, high-, intermediate-, and low-sulfidation epithermal Au systems, and Au-rich carbonate replacement and sediment-hosted styles. The differences in endowment and deposit styles likely reflect regional-scale tectono-magmatic processes as well as local preservation and emplacement levels. The Cu ± Au endowment of the Cretaceous is consistent with typical subduction-related arc environments and generation of calc-alkaline porphyry to high-sulfidation epithermal systems, whereas Au enrichment related to Cenozoic magmatism appears to be related to high-K calc-alkalic to shoshonitic compositions. In many of the Au-rich Cenozoic magmatic belts, there is geochemical evidence for sourcing subcontinental lithospheric mantle that was previously enriched by Cretaceous subduction-related metasomatism. Additional differences in Au endowment may reflect the preservation of shallow-level systems in the Cenozoic, particularly for the Au-rich Miocene porphyry deposits such as Kışladağ and Bierly Vrch and the Apuseni porphyry Au-Cu deposits. However, in both the Cretaceous and Cenozoic, crustal exposure levels vary across the belt and cannot explain all the differences in Cu and Au endowment.A compilation of exploration discovery methods highlights the importance of historic workings in addition to geochemistry and geology as an initial vector, whereas geophysics has had limited involvement in direct discovery, primarily due to its limited application historically. Geologic models for well-understood systems such as porphyry and proximal epithermal systems provide excellent guides for explorers; however, more distal deposits such as Au-rich carbonate replacement deposits and deposits with poorly constrained models such as sedimentary rock-hosted and intermediate-sulfidation deposits are more challenging for exploration.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The newly defined, 1,500-km-long, late Cenozoic Anatolian metallogenic trend of Turkey is the central segment of the western Tethyan metallogenic belt and formed after the closure of the southern Neotethys Ocean. Mineral deposit discoveries along this trend show that the Oligocene to Miocene igneous units are highly prospective for gold-rich porphyry- and epithermal-style mineralization (~27 Moz) but that copper endowment is poor. However, the temporal and spatial constraints on late Cenozoic gold districts and isolated prospects and their tectonic affinity are poorly known, despite recent efforts. We herein provide new U-Pb and Re-Os age data and field observations from Miocene gold prospects and deposits throughout the Anatolian trend, which we interpret together with previously published age data in the region. We define nine new porphyry and epithermal districts: Simav, İzmir, Uşak, Bodrum, Konya, Aksaray, Kayseri, Tunceli, and Ağri. Gold-rich porphyry and epithermal systems peaked at (1) 25 to 17 Ma in eastern Anatolia, (2) 21 to 9 Ma in western Anatolia, and (3) 10 to 3 Ma in central Anatolia. The westward migration of porphyry and epithermal mineralization from eastern to central Anatolia is interpreted to reflect slab break-off propagation and gap opening after the onset of the Arabian continental collision. Conversely, the southwestward migration of the magmatic front and associated mineralization in western Anatolia resulted from the acceleration of the Aegean slab rollback and subsequent lateral tearing (15–8 Ma). Thus, the bulk of gold mineralization formed in response to the slab segmentation and thermal events at 25 and 15 Ma.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The formation of ore deposits in the Lavrion Pb-Zn-Ag district was associated with Miocene detachment that accommodated orogenic collapse and exhumation of high-grade nappes across the ductile-brittle transition. This district consists of (1) low-grade porphyry Mo style, (2) Cu-Fe skarn, (3) high-temperature carbonate replacement Pb-Zn-Ag, and (4) vein and breccia Pb-Zn-Ag mineralization. The vein and breccia mineralization locally contains high-grade silver in base metal sulfides that are cemented by fluorite and carbonate gangue. The rare earth element contents of these gangue minerals, chondrite-normalized patterns, and fluid inclusion studies suggest that they precipitated from a low-temperature hydrothermal fluid. Primary and pseudosecondary fluid inclusions in fluorite and calcite are characterized by a wide range of homogenization temperatures (92°–207°C) and salinities of up to 17.1 wt % NaCl equiv. Secondary fluid inclusions only represent 〈5 vol % of the total fluid trapped. Fluids extracted from inclusions in fluorite have values of 〈span〉δ〈/span〉D = –82.1 to –47.7‰ (Vienna-standard mean ocean water [V-SMOW]) and 〈span〉δ〈/span〉〈sup〉18〈/sup〉O = –10.4 to –5.1‰ (V-SMOW). These data and low ratios of Cl/Br measured by crush-leach analyses for fluids in fluorite (102–315) and calcite (162–188) are compatible with the ore fluid being the result of mixing of meteoric water with evaporated seawater. These data suggest that fluids leading to the deposition of late Pb-Zn-Ag–rich vein- and breccia-style mineralization in Lavrion were related to circulation of mixed evaporated seawater and meteoric fluids that was enhanced by brittle deformation. This contrasts with the fluids of magmatic origin related to the formation of low-grade porphyry Mo, Cu-Fe skarn, and high-temperature carbonate replacement deposits spatially related to the Plaka granodiorite.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The calc-alkalic Mitchell Au-Cu-Ag-Mo porphyry deposit, hosted in intrusive rocks of the Stikine volcanic arc terrane of northwestern British Columbia, is the largest undeveloped gold resource in Canada, with 40.72 Moz of contained gold. It lies within the KSM trend, a 12-km-long linear porphyry array in the Sulphurets district. It is genetically related to Early Jurassic Sulphurets stocks: phase 1 diorite to monzodiorite hosts Cu-Au mineralization in potassic assemblages (stage 1), a phase 2 granodiorite plug cores a molybdenum halo (190.3 ± 0.8 Ma, 191.3 ± 0.7 Ma; Re-Os, molybdenite) that is accompanied by phyllic alteration (stage 2), and a poorly mineralized phase 3 diorite plug temporally overlaps with quartz-pyrophyllite alteration at shallow levels (stage 3). Two deformation events (D〈sub〉1〈/sub〉 and D〈sub〉2〈/sub〉), attributed to orogen-scale mid-Cretaceous transpression, structurally modified Mitchell. D〈sub〉1〈/sub〉 deformation, expressed as steep E-striking pressure solution cleavage (S〈sub〉1〈/sub〉) and related W-plunging folded veins (F〈sub〉1〈/sub〉), is heterogeneously developed as a function of alteration type. D〈sub〉2〈/sub〉 is divided into two progressive events: D〈sub〉2a〈/sub〉, defined by N-plunging, folded veins (F〈sub〉2a〈/sub〉), and D〈sub〉2b〈/sub〉, defined by E-vergent thrust faults, including the Mitchell thrust fault, which offsets the updip continuation of Mitchell (the Snowfield deposit, 192.0 ± 1.0 Ma, 191.1 ± 0.8 Ma; Re-Os, molybdenite) ~1,600 m to the east-southeast. Host structures for the KSM trend may have been long-lived, N-striking basement lineaments that provided transcrustal magma and fluid pathways. East-trending intrusions, hydrothermal veins, alteration and metal distribution at Mitchell are attributed to subsidiary E-striking cross faults. These original anisotropies in turn influenced the geometry of Cretaceous faults and flattening domains within the deposit.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Regional exploration for rare and precious metal enrichment traditionally uses multielement geochemical analysis of the fines (〈150 〈span〉μ〈/span〉m) fraction of stream sediments; however, these data can be misleading in glaciated regions with complex geology. Here, we compare the spatial distributions of data from sediment fines with distributions of heavy mineral concentrates, including detrital gold, in the same region. Gold grains are characterized according to abundance and morphology plus a microchemical signature from the combination of alloy analysis and systematic identification of mineral inclusions revealed in polished sections. These inclusions survive indefinitely and are indicative of hypogene mineralogy, whereas the metal loadings of fines may be affected by weathering or anthropogenic activity. All three approaches are mutually supportive: sediment fines analyses provide a basis for more labor-intensive targeted gold-grain studies, which in turn highlight specific gold-element associations useful for interpretation of geochemical data sets. Spatial distributions of resistate heavy mineral suites constrain the directions and extents of glacial transport, which facilitates more confident interpretations of placer-lode relationships from gold-grain studies.Characterization of 2,160 gold grains from 40 localities in the auriferous region of southeast Ireland provided a clear indication of proximity of gold to source and identified gold derived from different episodes of mineralization. A distinction is apparent between gold in the south of the region (Wexford), likely derived from widespread stratabound Au-As-Fe-S mineralization, and that in the north (Wicklow), where the historical placer mining district of the Goldmines River yielded gold with inclusions exhibiting a distinctive Pb-Bi-As association. The Goldmines River placers formed by the efficient accumulation and preservation of detrital gold derived from several discrete intravalley sources. We recommend that a combination of classic stream-sediment geochemistry, heavy mineral analysis, and gold-grain studies is used more widely to provide additional insights on the location and nature of gold mineralization and regional metallogeny in regions of poor exposure and complex geology.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The giant Pulang porphyry Cu-Au district (446.8 million tonnes at 0.52% Cu and 0.18 g/tonne Au) is in the southern segment of the Yidun arc (Zhongdian arc), part of the Sanjiang Tethyan orogenic belt in southwest China. The district consists of three deposits: South Pulang (~96% of the total ore reserves) and the smaller East and North Pulang deposits. Four intrusive phases host the three Pulang deposits, which are, in order of emplacement, premineralization fine-grained quartz diorite and coarse-grained quartz diorite, intermineralization quartz monzonite, and late-mineralization diorite porphyry. The complex intruded carbonaceous rocks of the Late Triassic Tumugou slates. Zircon U-Pb laser ablation-inductively coupled plasma-mass spectrometry dating shows that intrusive activity occurred at about 216 ± 2 Ma.Hydrothermal alteration of the intrusions at Pulang includes five main types: K-silicate, epidote-chlorite, chlorite-illite, quartz-illite, and clay alteration. K-silicate alteration, subdivided into early K-feldspar alteration and late biotite alteration (dominant), mainly affected the central quartz monzonite and adjacent coarse-grained quartz diorite and fine-grained quartz diorite. Epidote-chlorite alteration, the most widespread alteration in the district, extends from the deposit core outward and has overprinted K-silicate alteration at South and North Pulang. Late chlorite-illite, quartz-illite, and clay alteration have overprinted preexisting K-silicate and epidotechlorite alteration assemblages and are locally developed in all four intrusive phases. Copper and gold are positively correlated and are mainly (90%) associated with epidote-chlorite alteration and, to a lesser degree, with K-silicate and chlorite-illite alteration. Hypogene pyrrhotite is intergrown with chalcopyrite and mainly occurs in chlorite-illite– and quartz-illite–altered, coarse-grained quartz diorite at East Pulang. Molybdenite Re-Os dating shows that mineralization in the district occurred at 216.54 ± 0.87 to 216.13 ± 0.86 Ma.The sequence of intrusion emplacement, alteration and veining, and sulfide associations at the three deposits suggests that South and North Pulang are two separate porphyry Cu-Au deposits, whereas East Pulang is probably a distal part of South Pulang. The dominance of primary magnetite over ilmenite and the assemblage titanite + magnetite + quartz in the causative quartz monzonite, and the abundant hydrothermal anhydrite veins associated with early K-silicate and main-mineralization epidote-chlorite alteration indicate the oxidized nature of the felsic intrusion and resultant early hydrothermal fluids. The pyrrhotite related to late chlorite-illite and quartz-illite alteration suggests local reduction due to interaction with the carbonaceous Tumugou slates. The atypical association of epidote-chlorite alteration and Cu mineralization at Pulang either is due to fluids from another porphyry deposit nearby overprinting epidote-chlorite alteration onto preexisting copper mineralization and K-silicate alteration at Pulang or is the result of collapse of epidote-chlorite–stable fluids into the K-silicate-altered core during waning hydrothermal activity.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Carbonate-hosted hydrothermal deposits typically show narrow visible mineralogical and textural alteration halos, which inhibit exploration targeting. In contrast, hydrothermal modification of the country rock’s stable isotope composition usually extends far beyond the limited visible alteration. Hence, stable isotope studies should be an effective tool to aid exploration for carbonate-hosted deposits. Here we present new insight into the development of a large stable isotope alteration halo based on 910 O and C isotope analyses of carbonate veins and hydrothermally altered limestone hosting the Cinco de Mayo Pb-Zn-Ag (Au, Cu) carbonate replacement deposit (CRD), in Chihuahua, Mexico. Our results demonstrate that stable isotope alteration is consistent with reactive, magmatic fluid flow into unaltered limestone and represents a powerful tool for the characterization of these hydrothermal ore systems. Synmineralization veins are texturally and isotopically distinct from those formed during pre- and postmineralization diagenesis and fluid flow and show distinct gradients along the direction of mineralizing fluid flow: this appears to be a promising exploration vectoring tool. Downhole variations in wall-rock isotope values reveal aquifers and aquicludes and outline the principal hydrothermal flow paths. Furthermore, wall-rock 〈span〉δ〈/span〉〈sup〉18〈/sup〉O〈sub〉VSMOW〈/sub〉 systematically decreases toward mineralization from ~23‰ to 〈17‰ over a distance of ~10 km, providing another vectoring tool. The extent of the stable isotope alteration halo likely reflects the overall fluid volume and areal extent of a fossil hydrothermal system, which may be expected to scale with the mineral endowment. This suggests that constraining the size, shape, and degree of isotopic alteration has direct application to mineral exploration by outlining the system and indicating the potential size of a deposit.〈/span〉