skip to main content

Title: Cosmogenic 3He production rate in ilmenite and the redistribution of spallation 3He in fine-grained minerals
Cosmogenic nuclide surface exposure dating and erosion rate measurements in basaltic landscapes rely primarily on measurement of 3He in olivine or pyroxene. However, geochemical investigations using 3He have been impossible in the substantial fraction of basalts that lack separable olivine or pyroxene crystals, or where such crystals were present, but have been chemically weathered. Fine-textured basalts often contain small grains of ilmenite, a weathering-resistant mineral that is a target for cosmogenic 3He production with good He retention and straightforward mineral separation, but with a poorly constrained production rate. Here we empirically calibrate the cosmogenic 3He production rate in ilmenite by measuring 3He concentrations in basalts with fine-grained (~20 lm cross-section) ilmenite and co-existing pyroxene or olivine from the Columbia River and Snake River Plain basalt provinces in the western United States. The concentration ratio of ilmenite to pyroxene and olivine is 0.78 ± 0.02, yielding an apparent cosmogenic 3He production rate of 93.6 ± 7.7 atom g-1 yr-1 that is 20–30% greater than expected from prior theoretical and empirical estimates for compositionally similar minerals. The production rate discrepancy arises from the high energy with which cosmic ray spallation reactions emit tritium and 3He and the associated long stopping distances that cause more » them to redistribute within a rock. Fine-grained phases with low cosmogenic 3He production rates, like ilmenite, will have anomalously high production rates owing to net implantation of 3He from the surrounding, higher 3He production rate, matrix. Semi-quantitative modeling indicates implantation of spallation 3He increases with decreasing ilmenite grain size, leading to production rates that exceed those in a large grain by ~10% when grain radii are <150 lm. The modeling predicts that for the ilmenite grain size in our samples, implantation causes production rates to be ~20% greater than expected for a large grain, and within uncertainty resolves the discrepancy between our calibrated production rate, theory, and rates from previous work. The redistribution effect is maximized when the host rock and crystals differ substantially in mean atomic number, as they do between whole-rock basalt and ilmenite. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Geochimica et cosmochimica acta
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this

    Quaternary lavas of the Stardalur Caldera, 20 km northeast of Reykjavik, Iceland, create a 27 300 nT magnetic anomaly visible in both ground and aeromagnetic surveys. Here, we provide a comprehensive mineralogical and rock magnetic data set to analyse NRM intensities and Koenigsberger ratios of 57 drill-core samples from the critical zone (CZ) of the anomaly high at depths between 41 and 131 m. This extends previous studies and verifies that the anomaly is due to an unusually high intensity of remanent magnetization carried by magnetite. The NRM of the CZ samples was acquired during the Olduvai subchron in a field of at most today’s strength. NRM intensities range from 20 to 128 A m–1 with a median of 55 A m–1, and an average of 61 A m–1, respectively, approximately 13–15 times higher than in typical Icelandic basalts (AIB) with an NRM intensity of 4 A m–1. Our new data set shows that the magnetite concentration throughout the CZ basalts is at most twofold higher than in AIB lavas. New data on domain state and TRM efficiency prove that these properties account for an additional factor of at most 2.3. Because magnetite is the most abundant remanence carrier in rocks on Earth, and its remanence acquisition is considered to be extremely wellmore »understood, we assert that the remaining discrepancy is a critical enigma in rock magnetism. Results from scanning electron microscopy show that a significant fraction of all CZ magnetite particles have dendritic shapes with grain sizes <1 μm, indicating rapid crystallization. Most large magnetite grains are heavily subdivided by very fine oxidation-exsolution lamellae of ilmenite, and subordinate amount of exsolved spinel as needles, blebs and blades. These common microstructures found throughout the CZ subdivide the initially homogeneous mineral into separate cubicles, here denoted as compartments. The magnetite compartments then have sizes below 1 μm. Hysteresis data, Preisach maps and FORC data consistently confirm that the coercivity distribution is dominated by values above 10 mT, such that multidomain behaviour is of little relevance in the CZ. Between 5 and 20 per cent of the IRM is carried by coercivities above 100 mT, which for magnetite indicates unusually high anisotropy effects in the individual particles. Based on the quantitative analysis of all magnetic contributions to the NRM, we can demonstrate that the average efficiency of NRM acquisition in the CZ Stardalur basalts must be at least a factor 3 higher than in typical basalts. We speculate that this is related to the observed focused compartment size distribution <1 μm, and indicates thermochemical remanence acquisition below the Curie temperature of magnetite. Yet, a detailed physical mechanism for the extreme overefficiency of NRM acquisition remains enigmatic.

    « less
  2. Hydrogen is a rapidly diffusing monovalent cation in nominally anhydrous minerals (NAMs, such as olivine, orthopyroxene, and clinopyroxene), which is potentially re-equilibrated during silicate melt-rock and aqueous fluid-rock interactions in massif and abyssal peridotites. We apply a 3D numerical diffusion modeling technique to provide first-order timescales of complete hydrogen re-equilibration in olivine, clinopyroxene, and orthopyroxene over the temperature range 600-1200°C. Model crystals are 1-3 mm along the c-axis and utilize H+ diffusion coefficients appropriate for Fe-bearing systems. Two sets of models were run with different boundary compositions: 1) “low-H models” are constrained by mineral-melt equilibrium partitioning with a basaltic melt that has 0.75 wt% H2O and 2) “high-H models,” which utilize the upper end of the estimated range of mantle water solubility for each phase. Both sets of models yield re-equilibration timescales that are identical and are fast for all phases at a given temperature. These timescales have strong log-linear trends as a function of temperature (R2 from 0.97 to 0.99) that can be used to calculate expected re-equilibration time at a given temperature and grain size. At the high end of the model temperatures (1000-1200°C), H+ completely re-equilibrates in olivine, orthopyroxene, and clinopyroxene within minutes to hours, consistent withmore »previous studies. These short timescales indicate that xenolith NAM mantle water contents are likely to be overprinted prior to eruption. The models also resolve the decoupled water-trace element relationship in Southwest Indian Ridge peridotites, in which peridotite REE abundances are reproduced by partial melting models whereas the relatively high NAM H2O contents require later re-equilibration with melt. At temperatures of 600-800°C, which correspond to conditions of hydrothermal alteration of pyroxene to amphibole and talc, H+ re-equilibration typically occurs over a range of timescales spanning days to years. These durations are well within existing estimates for the duration of fluid flow in oceanic hydrothermal systems, suggesting that peridotite NAM water contents are susceptible to diffusive overprinting during higher temperature hydrothermal alteration. Thus, diffusion during aqueous fluid-rock interactions may also explain NAM H2O contents that are too high to reflect residues of melting. These relatively short timescales at low temperatures suggest that the origin of water contents measured in peridotite NAMs requires additional constraints on sample petrogenesis, including petrographic and trace element analyses. Our 3D model results also hint that H+ may diffuse appreciably during peridotite serpentinization, but diffusion coefficients at low temperature are unconstrained and additional experimental investigations are needed.« less
  3. Abstract We investigate the shallow plumbing system of the Deccan Traps Large Igneous Province using rock and mineral data from Giant Plagioclase Basalt (GPB) lava flows from around the entire province, but with a focus on the Saurashtra Peninsula, the Malwa Plateau, and the base and top of the Western Ghats (WG) lava pile. GPB lavas in the WG typically occur at the transition between chemically distinct basalt formations. Most GPB samples are evolved basalts, with high Fe and Ti contents, and show major and trace elements and Sr-Nd-Pb isotopic compositions generally similar to those of previously studied Deccan basalts. Major element modeling suggests that high-Fe, evolved melts typical of GPB basalts may derive from less evolved Deccan basalts by low-pressure fractional crystallization in a generally dry magmatic plumbing system. The basalts are strongly porphyritic, with 6–25% of mm- to cm-sized plagioclase megacrysts, frequently occurring as crystal clots, plus relatively rare olivine and clinopyroxene. The plagioclase crystals are mostly labradoritic, but some show bytownitic cores (general range of anorthite mol%: 78–55). A common feature is a strong Fe enrichment at the plagioclase rims, indicating interaction with an Fe-rich melt similar to that represented by the matrix compositions (FeOt up tomore »16–17 wt%). Plagioclase minor and trace elements and Sr isotopic compositions analyzed by laser ablation inductively coupled plasma mass spectrometry show evidence of a hybrid and magma mixing origin. In particular, several plagioclase crystals show variable 87Sr/86Sri, which only partially overlaps with the 87Sr/86Sri of the surrounding matrix. Diffusion modeling suggests residence times of decades to centuries for most plagioclase megacrysts. Notably, some plagioclase crystal clots show textural evidence of deformation as recorded by electron back-scatter diffraction analyses and chemical maps, which suggest that the plagioclase megacrysts were deformed in a crystal-rich environment in the presence of melt. We interpret the plagioclase megacrysts as remnants of a crystal mush originally formed in the shallow plumbing system of the Deccan basalts. In this environment, plagioclase acquired a zoned composition due to the arrival of chemically distinct basaltic magmas. Prior to eruption, a rapidly rising but dense Fe-rich magma was capable of disrupting the shallow level crystal mush, remobilizing part of it and carrying a cargo of buoyant plagioclase megacrysts. Our findings suggest that basaltic magmas from the Deccan Traps, and possibly from LIPs in general, are produced within complex transcrustal magmatic plumbing systems with widespread crystal mushes developed in the shallow crust.« less
  4. International Ocean Discovery Program Expedition 397T sought to address the shortage of drilling time caused by COVID-19 mitigation during Expedition 391 (Walvis Ridge Hotspot) by drilling at two sites omitted from the earlier cruise. A week of coring time was added to a transit of JOIDES Resolution from Cape Town to Lisbon, which would cross Walvis Ridge on its way north. These two sites were located on two of the three seamount trails that emerge from the split in Walvis Ridge morphology into several seamount chains at 2°E. Site U1584 (proposed Site GT-6A) sampled the Gough track on the east, and Site U1585 (proposed Site TT-4A) sampled the Tristan track on the west. Together with Site U1578, drilled on the Center track during Expedition 391, they form a transect across the northern Walvis Ridge Guyot Province. The goal was to core seamount basalts and associated volcanic material for geochemical and isotopic, geochronologic, paleomagnetic, and volcanologic study. Scientifically, one emphasis was to better understand the split in geochemical and isotopic signatures that occurs at the morphologic split. Geochronology would add to the established age progression but also give another dimension to understanding Walvis Ridge seamount formation by giving multiple ages atmore »the same sites. The paleomagnetic study seeks to establish paleolatitudes for Walvis Ridge sites for comparison with those published from hotspot seamount chains in the Pacific, in particular to test whether a component of true polar wander affects hotspot paleolatitude. Hole U1584A cored a 66.4 m thick sedimentary and volcaniclastic section with two lithostratigraphic units. Unit I is a 23 m thick sequence of bioturbated clay and nannofossil chalk with increasing volcaniclastic content downhole. Unit II is a >43 m thick sequence of lapillistone with basalt fragments. Because the seismic section crossing the site shows no evidence as to the depth of the volcaniclastic cover, coring was terminated early. Because there were no other shallow nearby sites with different character on existing seismic lines, the unused operations time from Site U1584 was shifted to the next site. The seismic reflector interpreted as the top of igneous rock at Site U1585 once again resulted from volcaniclastic deposits. Hole U1585A coring began at 144.1 mbsf and penetrated a 273.5 m thick sedimentary and volcaniclastic section atop a 81.2 m thick series of massive basalt flows. The hole was terminated at 498.8 mbsf because allotted operational time expired. The sedimentary section contains four main units. Unit I (144.1–157.02 mbsf) is a bioturbated nannofossil chalk with foraminifera, similar to the shallowest sediments recovered at Site U1584. Unit II (157.02–249.20 mbsf), which is divided into two subunits, is a 92.2 m thick succession of massive and bedded pumice and scoria lapillistone with increased reworking, clast alteration, and tuffaceous chalk intercalations downhole. Unit III (249.20–397.76 mbsf) is 148.6 m thick and consists of a complex succession of pink to greenish gray tuffaceous chalk containing multiple thin, graded ash turbidites and tuffaceous ash layers; intercalated tuffaceous chalk slumps; and several thick coarse lapilli and block-dominated volcaniclastic layers. Befitting the complexity, it is divided into eight subunits (IIIA–IIIH). Three of these subunits (IIIA, IIID, and IIIG) are mainly basalt breccias. Unit IV (397.76–417.60 mbsf) is a volcanic breccia, 19.8 m thick, containing mostly juvenile volcaniclasts. The igneous section, Unit V (417.60–498.80 mbsf) is composed of a small number of massive basaltic lava flows. It is divided into three lithologic units, with Unit 2 represented by a single 3 cm piece of quenched basalt with olivine phenocrysts in a microcrystalline groundmass. This piece may represent a poorly recovered set of pillow lavas. Unit 1 is sparsely to highly olivine-clinopyroxene ± plagioclase phyric massive basalt and is divided into Subunits 1a and 1b based on textural and mineralogical differences, which suggests that they are two different flows. Unit 3 also consists of two massive lava flows with no clear boundary features. Subunit 3a is a 10.3 m thick highly clinopyroxene-plagioclase phyric massive basalt flow with a fine-grained groundmass. Subunit 3b is a featureless massive basalt flow that is moderately to highly clinopyroxene-olivine-plagioclase phyric and >43.7 m thick. Alteration of the lava flows is patchy and moderate to low in grade, with two stages, one at a higher temperature and one at a low temperature, both focused around fractures. The Site U1585 chronologic succession from basalt flows to pelagic sediment indicates volcanic construction and subsidence. Lava eruptions were followed by inundation and shallow-water volcaniclastic sediment deposition, which deepened over time to deepwater conditions. Although the massive flows were probably erupted in a short time and have little variability, volcaniclasts in the sediments may provide geochemical and geochronologic data from a range of time and sources. Chemical analyses indicate that Site U1585 basalt samples are mostly alkalic basalt, with a few trachybasalt flow and clast samples and one basaltic trachyandesite clast. Ti/V ratios lie mostly within the oceanic island basalt (OIB) field but overlap the mid-ocean-ridge basalt (MORB) field. Only a handful of clasts from Site U1584 were analyzed, but geochemical data are similar. Paleomagnetic data from Site U1585 indicate that the sediments and basalt units are strongly magnetic and mostly give coherent inclination data, which indicates that the basaltic section and ~133 m of overlying volcaniclastic sediment is reversely polarized and that this reversal is preserved in a core. Above this, the rest of the sediment section records two normal and two reversed zones. Although there are not enough basalt flows to give a reliable paleolatitude, it may be possible to attain such a result from the sediments.« less
  5. Abstract The Bishop Tuff (BT), erupted from the Long Valley caldera in California, displays two types of geochemical gradients with temperature: one is related to magma mixing, whereas the other is found in the high-SiO2 rhyolite portion of the Bishop Tuff and is characterized by twofold or lower concentration variations in minor and trace elements that are strongly correlated with temperature. It is proposed that the latter zonation, which preceded phenocryst growth, developed as a result of mineral–melt partitioning between interstitial melt and surrounding crystals in a parental mush, from which variable melt fractions were segregated. To test this hypothesis, trends of increasing vs decreasing element concentrations with temperature (as a proxy for melt fraction), obtained from published data on single-clast pumice samples from the high-SiO2 rhyolite portion of the Bishop Tuff, were used to infer their relative degrees of incompatibility vs compatibility between crystals and melt in the parental mush. Relative compatibility values (RCVi) for all elements i, defined as the concentration slope with temperature divided by average concentration, are shown to be linearly correlated with their respective bulk partition coefficients (bulk Di). Mineral–melt partition coefficients from the literature were used to constrain the average stoichiometry of the crystallization/meltingmore »reaction in the parental mush: 32 % quartz + 34 % plagioclase + 31 % K-feldspar + 1·60 % biotite + 0·42 % titanomagnetite + 0·34 % ilmenite + 0·093 % allanite + 0·024 % zircon + 0·025 % apatite = 100 % liquid. The proportions of tectosilicates in the reaction (i.e. location of eutectic) are consistent with depths of melt segregation of ~400–550 MPa and an activity of H2O of ~0·4–0·6. Temperatures of <770–780 °C are constrained by allanite in the reaction. Evidence that a fluid phase was present in the parental mush is seen in the decreasing versus increasing H2O and CO2 contents with temperature in the segregated interstitial melt that formed the high-SiO2 rhyolite portion of the Bishop Tuff. The presence of an excess fluid phase, which strongly partitions CO2 relative to the melt, is required to explain the compatible behavior of CO2, whereas the fluid abundance must have been low to explain the incompatible behavior of H2O. Calculated degassing paths for interstitial melts, which segregated from the parental mush and ascended to shallower depths to grow phenocrysts, match published volatile analyses in quartz-hosted melt inclusions and constrain fluid abundances in the mush to be ≤1 wt%. The source of volatiles in the parental mush, irrespective of whether it formed by crystallization or partial melting, must have been primarily from associated basalts, as granitoid crust is too volatile poor. Approximately twice as much basalt as rhyolite is needed to provide the requisite volatiles. The determination of bulk Di for several elements gives the bulk composition of the parental leucogranitic mush and shows that it is distinct from Mesozoic Sierran arc granitoids, as expected. Collectively, the results from this study provide new constraints for models of the complex, multi-stage processes throughout the Plio-Quaternary, involving both mantle-derived basalt and pre-existing crust, that led to the origin of the parental body to the Bishop Tuff.« less