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 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.  more » « less
Award ID(s):
1529528
NSF-PAR ID:
10120685
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Geochimica et cosmochimica acta
Volume:
265
ISSN:
0016-7037
Page Range / eLocation ID:
19-31
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. SUMMARY

    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 well 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.

     
    more » « less
  2. Abstract

    Caleta el Cobre (CeC) 022 is a Martian meteorite of the nakhlite group, showing an unbrecciated cumulate texture, composed mainly of clinopyroxene and olivine. Augite shows irregular core zoning, euhedral rims, and thin overgrowths enriched in Fe relative to the core. Low‐Ca pyroxene is found adjacent to olivine. Phenocrysts of Fe‐Ti oxides are titanomagnetite with exsolutions of ilmenite/ulvöspinel. Intercumulus material consists of both coarse plagioclase and fine‐grained mesostasis, comprising K‐feldspars, pyroxene, apatite, ilmenite, Fe‐Ti oxides, and silica. CeC 022 shows a high proportion of Martian aqueous alteration products (iddingsite) in olivine (45.1 vol% of olivine) and mesostasis. This meteorite is the youngest nakhlite with a distinct Sm/Nd crystallization age of 1.215 ± 0.067 Ga. Its ejection age of 11.8 ± 1.8 Ma is similar to other nakhlites. CeC 022 reveals contrasted cooling rates with similarities with faster cooled nakhlites, such as Northwest Africa (NWA) 817, NWA 5790, or Miller Range 03346 nakhlites: augite irregular cores, Fe‐rich overgrowths, fine‐grained K‐feldspars, quenched oxides, and high rare earth element content. CeC 022 also shares similarities with slower cooled nakhlites, including Nakhla and NWA 10153: pyroxene modal abundance, pyroxenes crystal size distribution, average pyroxene size, phenocryst mineral compositions, unzoned olivine, and abundant coarse plagioclase. Moreover, CeC 022 is the most magnetic nakhlite and represents an analog source lithology for the strong magnetization of the Martian crust. With its particular features, CeC 022 must originate from a previously unsampled sill or flow in the same volcanic system as the other nakhlites, increasing Martian sample diversity and our knowledge of nakhlites.

     
    more » « less
  3. 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 with 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. 
    more » « less
  4. 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 to 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. 
    more » « less
  5. SUMMARY The occurrence of plate tectonics on Earth is rooted in the physics of lithospheric ductile weakening and shear-localization. The pervasiveness of mylonites at lithospheric shear zones is a key piece of evidence that localization correlates with reduction in mineral grain size. Most lithospheric mylonites are polymineralic and the interaction between mineral phases, such as olivine and pyroxene, especially through Zener pinning, impedes normal grain growth while possibly enhancing grain damage, both of which facilitate grain size reduction and weakening, as evident in lab experiments and field observations. The efficacy of pinning, however, relies on the mineral phases being mixed and dispersed at the grain scale, where well-mixed states lead to greater mylonitization. To model grain mixing between different phases at the continuum scale, we previously developed a theory treating grain-scale processes as diffusion between phases, but driven by imposed compressive stresses acting on the boundary between phases. Here we present a new model for shearing rock that combines our theory for diffusive grain mixing, 2-D non-Newtonian flow and two-phase grain damage. The model geometry is designed specifically for comparison to torsional shear-deformation experiments. Deformation is either forced by constant velocity or constant stress boundary conditions. As the layer is deformed, mixing zones between different mineralogical units undergo enhanced grain size reduction and weakening, especially at high strains. For constant velocity boundary experiments, stress drops towards an initial piezometric plateau by a strain of around 4; this is also typical of monophase experiments for which this initial plateau is the final steady state stress. However, polyphase experiments can undergo a second large stress drop at strains of 10–20, and which is associated with enhanced phase mixing and resultant grain size reduction and weakening. Model calculations for polyphase media with grain mixing and damage capture the experimental behaviour when damage to the interface between phases is moderately slower or less efficient than damage to the grain boundaries. Other factors such as distribution and bulk fraction of the secondary phase, as well as grain-mixing diffusivity also influence the timing of the second stress drop. For constant stress boundary conditions, the strain rate increases during weakening and localization. For a monophase medium, there is theoretically one increase in strain rate to a piezometric steady state. But for the polyphase model, the strain rate undergoes a second abrupt increase, the timing for which is again controlled by interface damage and grain mixing. The evolution of heterogeneity through mixing and deformation, and that of grain size distributions also compare well to experimental observations. In total, the comparison of theory to deformation experiments provides a framework for guiding future experiments, scaling microstructural physics to geodynamic applications and demonstrates the importance of grain mixing and damage for the formation of plate tectonic boundaries. 
    more » « less