Tourmaline boron and strontium isotope systematics reveal magmatic fluid pulses and external fluid influx in a giant iron oxide-apatite (IOA) deposit

Zhi-kun Su; Xinfu Zhao; Li-ping Zeng; Kui-dong Zhao; Albert H. Hofstra

doi:10.1016/j.gca.2019.06.006

Home/Archive/USGS

USGS · 70223691

Tourmaline boron and strontium isotope systematics reveal magmatic fluid pulses and external fluid influx in a giant iron oxide-apatite (IOA) deposit

Zhi-kun Su·Xinfu Zhao·Li-ping Zeng·Kui-dong Zhao·Albert H. Hofstra

View original source ↗Export metadata ↓

Abstract

Tourmaline is a common boron-bearing mineral in hydrothermal system and has been widely used as a mineral probe to reconstruct geological processes because of its broad range in composition and resistance to metasomatic alteration. The origin of Kiruna-type iron oxide-apatite (IOA) deposits, commonly linked to andesitic subvolcanic or volcanic rocks, is highly controversial. Constraints on the evolution of these mineralizing systems are needed to advance understanding of the ore-forming process. In this study, we apply in situ elemental and combined B-Sr isotopic analyses of tourmaline to elucidate the nature and evolution of the subsurface hydrothermal system associated with IOA mineralization in the giant Taocun deposit, eastern China. Taocun is hosted at the top of a diorite intrusion and exhibits three stages of hydrothermal alteration that contain tourmaline: pre-ore Na alteration (Tur I), syn-ore magnetite formation and associated Ca-Fe alteration (Tur II), and post-ore Ca-Mg alteration with sulfide veins (Tur III). Compositional data for each stage of tourmaline plot along the “oxy-dravite”–povondraite join, which is indicative of precipitation from relatively oxidizing fluids. Ranges of Sr-isotopic compositions in Tur I (0.7065–0.7078) and Tur II (0.7068–0.7076) are identical to those of the igneous host rocks, indicating precipitation from magmatic-hydrothermal fluids. The range of B-isotopic compositions in Tur I (δ 11 B values of −6.3‰ to −1.2‰) is also consistent with a magmatic source. Higher δ 11 B values (−2.4‰ to 5.4‰) obtained from Tur II are mainly ascribed to Rayleigh fractionation in the magmatic-hydrothermal system as tourmaline precipitated. Post-ore Tur III has a wide range of mostly lower B-isotopic compositions (−8.5‰ to 0.8‰) that record another pulse of magmatic fluid input. This interpretation is supported by the enrichment of Na, Li, Be, W, Sn, V, and Ti in Tur III, relative to Tur I and II. However, the higher Sr-isotope composition (0.7076–0.7086) of Tur III and available O-isotope composition (−7‰ to 3.5‰) of fluids of this stage record the infiltration of meteoric ground water from adjacent sedimentary country rocks. The results suggest that the Taocun IOA deposit formed in a magmatic-hydrothermal system characterized by two (or more) pulses of magmatic fluid discharge from subvolcanic diorite intrusions, followed by the influx of external ground water as the system waned. This study highlights the utility of tourmaline as a robust geochemical and isotopic monitor of ore-forming processes in such systems.

Explore related subjects

Geochimica et Cosmochimica Acta Article

Source-reported bounding extent: 30.097613277217132° to 31.704803074739214° latitude; 117.48779296875° to 118.58642578124999° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Zhi-kun Su, Xinfu Zhao, Li-ping Zeng, Kui-dong Zhao, Albert H. Hofstra. 2019. Tourmaline boron and strontium isotope systematics reveal magmatic fluid pulses and external fluid influx in a giant iron oxide-apatite (IOA) deposit. https://doi.org/10.1016/j.gca.2019.06.006

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Organic and isotopic indicators for sorting of sedimentary organic matter along a marginal submarine canyon

Submarine canyons are incised features of many continental margins that can have significant influence on the hydrodynamic distribution of sediments and organic matter eroded and deposited from the continents. Baltimore Canyon, on the mid-Atlantic margin of the United States, contains a complex set of sedimentary processes that simultaneously create unique benthic habitats and control the deposition of organic matter. Along the canyon axis, loci of net erosion, net deposition, and intense winnowing each host diverse faunal assemblages and varying mixtures of sedimentary organic matter derived both from production in the overlying water column and from mobilized sediments. Bioavailable components of this deposited organic matter sustain benthic communities, while recalcitrant components can contribute to long-term carbon burial in the deep sea. However, commonly employed bulk geochemical analyses provide little information about the relative bioavailability or depositional history of sedimentary organic matter. Here we employ a range of organic and isotopic analyses to explore in more detail how canyon-specific sediment dynamics determine the sorting of organic matter from shelf to open ocean. In combination with bulk geochemical characteristics, we subjected surface sediments from water depths of ∼200–1200 m in Baltimore Canyon to a sequential extraction procedure, isolating nonpolar and polar lipid classes, an acid-soluble fraction, and an acid-insoluble fraction. Each class was analyzed for carbon and nitrogen quantities and stable isotope ratios, and radiocarbon content where possible, along with compound-specific carbon and nitrogen isotope analysis of individual amino acids in the acid-hydrolysed fraction. We find different organic matter sources and depositional history recorded in the properties of younger, bioavailable organic matter components (polar lipids, amino acids) in comparison to the older, more recalcitrant components (nonpolar lipids, acid-insoluble fraction). These differences in source and bioavailability of organic matter vary along the canyon, correlating with grain size and erosion/deposition dynamics, and may help shape the benthic faunal assemblages. Additionally, our results suggest that determining the relative concentrations of acid-soluble and acid-insoluble organic matter may provide an easily accessible method to improve our understanding of the nutritional quality of sediments for benthic fauna than more commonly used bulk carbon or nitrogen concentrations.

Delaware, New Jersey↗

USGS2025

Bayesian calibration of the 40K decay scheme with implications for 40K-based geochronology

The K/Ar and 40 Ar/ 39 Ar geochronometers are based on the naturally occurring radionuclide 40 K. Their precision and accuracy are limited by uncertainties on the 40 K decay constants and, in the case of the 40 Ar/ 39 Ar geochronometer, the isotopic composition of neutron fluence monitors. To address these limitations, we introduce a Bayesian calibration of the 40 K decay scheme. We formulate robust priors for all model parameters including partial 40 K decay constants, 238 U and 235 U decay constants, and age offset parameters to account for phenomena that can perturb apparent U-Pb and 40 Ar/ 39 Ar ages. We then harness a set of complementary 40 Ar/ 39 Ar, 238 U/ 206 Pb, and 235 U/ 207 Pb data from well- characterized geological samples with ages from 1.919 ka to 2000 Ma to derive Bayesian estimates of the 40 K decay constants. Posterior values for the partial 40 K decay constants are λ β - "> λ β - = (4.9252 ± "> ± 0.0054) × "> × 10 −10 yr −1 , λ β + "> λ β + = (5.6658 ± "> ± 0.1543) × "> × 10 −15 yr −1 , λ EC ∗ "> λ EC0 = (5.7404 ± "> ± 0.0053) × "> × 10 −11 yr −1 , and λ EC 0 "> λ EC0 = (4.9060 ± "> ± 0.2942) × "> × 10 −13 yr −1 (uncertainties reported at the 68 % (1 σ "> σ ) credible interval). These combine to a total 40 K decay constant λ tot "> λ tot = (5.5042 ± "> ± 0.0054) × "> × 10 −10 yr −1 . Model estimates of the 238 U and 235 U decay constants are statistically indistinguishable from those reported by Jaffey et al. (1971) . Posterior values of the 40 K decay constants and the 40 Ar*/ 40 K isotopic composition of Fish Canyon sanidine (FCs) define a K/Ar FCs age of 28.183 ± "> ± 0.017 Ma (1 σ "> σ ). Significantly, Bayesian calibrated 40 Ar/ 39 Ar ages align with astronomically tuned ages throughout the Cenozoic and with 238 U/ 206 Pb and 235 U/ 207 Pb ages in the Mesozoic, Paleozoic, and Proterozoic, as well as having comparable precision to the 238 U/ 206 Pb method. Thus, Bayesian calibration of the 40 K decay scheme and the K/Ar age of FCs reconciles the 40 Ar/ 39 Ar, U-Pb, and astronomical chronometers.

Geochimica et Cosmochimica Acta↗

USGS2024

Characterizing sulfur redox state and geochemical implications in deep-time using mineral chemistry network analysis

Sulfur (S) is a central element in global biogeochemical cycling and Earth’s redox evolution. Minerals that contain S are an important record of local environmental conditions at the time of their formation based on chemical speciation and redox. However, the oxidation state of S for hundreds of different S-containing minerals and thousands of S-containing mineral localities is unknown, largely sulfides and sulfosalts, and the redox state alone does not fully capture mineral chemistry diversity, thus limiting understanding of S redox evolution. Here, we use mineral chemistry network analysis and the weighted Mineral Element Electronegativity Coefficient of Variation (wMEE CV ) metric to investigate the element interactions and localities of S-containing minerals from the Mineral Evolution Database (MED) to infer the redox state of S in minerals where the redox state is unknown (S U ). Louvain community detection of the S mineral chemistry redox network reveals that there are three main network communities that are separated by redox state. The S 6+ community includes minerals that contain the S 6+ redox state and a small number of S 4+ and S 2+ minerals, the S 2− community includes S 2− -containing minerals, and the S U community includes minerals in which the redox state of S is unknown. The wMEE CV values of the S U community closely overlap with the wMEE CV values of the S 2− community, and do not overlap with the wMEE CV values of the S 6+ community, indicating the S U community minerals contain predominately reduced S. Assuming that S U community minerals contain reduced S, as supported by their network chemical associations and wMEE CV values, then reduced S-containing minerals make up approximately 81 % of S-containing mineral localities in the S mineral chemistry network, even though the majority of all mineral localities (S-containing and non-S-containing) are oxygen (O)-containing minerals. Additionally, reduced S-containing minerals make up the majority (∼75 %) of all non-O containing mineral localities in the MED, representing the importance of reduced S as an electron source and substrate in the evolution of microbial metabolic networks. The range wMEE CV values of S 6+ community minerals expands through time due primarily to formation of chemically diverse sulfate minerals, coinciding with crustal oxidation from the late Proterozoic to Phanerozoic and the expansion of the marine sulfate reservoir. The intersection of shared constituent elements among reduced and oxidized S in the mineral chemistry network represents redox convergence of weathered S in the geosphere that was crucial in the formation of natural resource deposits and the evolution of biogeochemical cycles.

Geochimica et Cosmochimica Acta↗