skip to main content


Title: K/U of the MORB Source and Silicate Earth
Abstract

Potassium (K) informs on the radiogenic heat production, atmospheric composition, and volatile element depletion of the Earth and other planetary systems. Constraints on the abundance of K in the Earth, Moon, and other rocky bodies have historically hinged on K/U values measured in planetary materials, particularly comparisons of the continental crust and mid‐ocean ridge basalts (MORBs), for developing compositional models of the bulk silicate Earth (BSE). However, a consensus on the most representative K/U value for global MORB remains elusive despite numerous studies. Here, we statistically analyze a critical compilation of MORB data to determine the K/U value of the MORB source. Covariations in the log‐normal abundances of K and U establish that K is 3–7 times less incompatible than U during melting and/or crystallization processes, enabling inverse modeling to infer the K/U of the MORB source region. These comprehensive data have a mean K/U for global MORB = 13,900 ± 200 (2σm;n = 4,646), and define a MORB source region with a K/U between 14,000 and 15,500, depending on the modeled melting regime. However, this range represents strictly a lower limit due to the undefined role of fractional crystallization in these samples and challenges preserving the signatures of depleted components in the MORB mantle source. This MORB source model, when combined with recent metadata analyses of ocean island basalt (OIB) and continental crust, suggests that the BSE has a K/U value >12,100 and contains >260 × 10−6 kg/kg K, resulting in a global production of3.5 TW of radiogenic heat today and 1.5 × 1017 kg of40Ar over the lifetime of the planet.

 
more » « less
Award ID(s):
1650365
NSF-PAR ID:
10374886
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
125
Issue:
12
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We report the Earth's rate of radiogenic heat production and (anti)neutrino luminosity from geologically relevant short‐lived radionuclides (SLR) and long‐lived radionuclides (LLR) using decay constants from the geological community, updated nuclear physics parameters, and calculations of theβspectra. We track the time evolution of the radiogenic power and luminosity of the Earth over the last 4.57 billion years, assuming an absolute abundance for the refractory elements in the silicate Earth and key volatile/refractory element ratios (e.g., Fe/Al, K/U, and Rb/Sr) to set the abundance levels for the moderately volatile elements. The relevant decays for the present‐day heat production in the Earth (19.9 ± 3.0 TW) are from40K,87Rb,147Sm,232Th,235U, and238U. Given element concentrations in kg‐element/kg‐rock and densityρin kg/m3, a simplified equation to calculate the present‐day heat production in a rock isurn:x-wiley:ggge:media:ggge22244:ggge22244-math-0001

    The radiogenic heating rate of Earth‐like material at solar system formation was some 103to 104times greater than present‐day values, largely due to decay of26Al in the silicate fraction, which was the dominant radiogenic heat source for the first10 Ma. Assuming instantaneous Earth formation, the upper bound on radiogenic energy supplied by the most powerful short‐lived radionuclide26Al (t1/2= 0.7 Ma) is 5.5×1031 J, which is comparable (within a factor of a few) to the planet's gravitational binding energy.

     
    more » « less
  2. Abstract

    Debate continues on the amount and distribution of radioactive heat producing elements (i.e., U, Th, and K) in the Earth, with estimates for mantle heat production varying by an order of magnitude. Constraints on the bulk‐silicate Earth's (BSE) radiogenic power also places constraints on overall BSE composition. Geoneutrino detection is a direct measure of the Earth's decay rate of Th and U. The geoneutrino signal has contributions from the local (40%) and global (35%) continental lithosphere and the underlying inaccessible mantle (25%). Geophysical models are combined with geochemical data sets to predict the geoneutrino signal at current and future geoneutrino detectors. We propagated uncertainties, both chemical and physical, through Monte Carlo methods. Estimated total signal uncertainties are on the order of20%, proportionally with geophysical and geochemical inputs contributing30% and70%, respectively. We find that estimated signals, calculated using CRUST2.0, CRUST1.0, and LITHO1.0, are within physical uncertainty of each other, suggesting that the choice of underlying geophysical model will not change results significantly, but will shift the central value by up to15%. Similarly, we see no significant difference between calculated layer abundances and bulk crustal heat production when using these geophysical models. The bulk crustal heat production is calculated as 7  2 TW, which includes an increase of 1 TW in uncertainty relative to previous studies. Combination of our predicted lithospheric signal with measured signals yield an estimated BSE heat production of 21.5  10.4 TW. Future improvements, including uncertainty attribution and near‐field modeling, are discussed.

     
    more » « less
  3. Fresh samples of basalts were collected by dredging from the Nanyue intraplate seamount in the Southwest sub-basin of the South China Sea (SCS). These are alkali basalts displaying right-sloping, chondrite-normalized rare earth element (REE) profiles. The investigated basalts are characterized by low Os content (60.37–85.13 ppt) and radiogenic 187Os/188Os ratios (~0.19 to 0.21). Furthermore, 40Ar/39Ar dating of the Nanyue basalts showed they formed during the Tortonian (~8.3 Ma) and, thus, are products of (Late Cenozoic) post-spreading volcanism. The Sr–Nd–Pb–Hf isotopic compositions of the Nanyue basalts indicate that their parental melts were derived from an upper mantle reservoir possessing the so-called Dupal isotopic anomaly. Semiquantitative isotopic modeling demonstrates that the isotopic compositions of the Nanyue basalts can be reproduced by mixing three components: the average Pacific midocean ridge basalt (MORB), the lower continental crust (LCC), and the average Hainan ocean island basalt (OIB). Our preferred hypothesis for the genesis of the Nanyue basalts is that their parental magmas were produced from an originally depleted mantle (DM) source that was much affected by the activity of the Hainan plume. Initially, the Hainan diapir caused a thermal perturbation in the upper mantle under the present-day Southwest sub-basin of the SCS that led to erosion of the overlying LCC. Eventually, the resultant suboceanic lithospheric mantle (SOLM) interacted with OIB-type components derived from the nearby Hainan plume. Collectively, these processes contributed crustal- and plume-type components to the upper mantle underlying the Southwest sub-basin of the SCS. This implies that the Dupal isotopic signature in the upper mantle beneath the SCS was an artifact of in situ geological processes rather than a feature inherited from a Southern Hemispheric, upper mantle source. 
    more » « less
  4. Abstract

    The southern US and northern Mexican Cordillera experienced crustal melting during the Laramide orogeny (c. 80–40 Ma). The metamorphic sources of melt are not exposed at the surface; however, anatectic granites are present throughout the region, providing an opportunity to investigate the metamorphic processes associated with this orogeny. A detailed geochemical and petrochronological analysis of the Pan Tak Granite from the Coyote Mountains core complex in southern Arizona suggests that prograde metamorphism, melting, and melt crystallization occurred here from 62 to 42 Ma. Ti-in-zircon temperatures (TTi-zr) correlate with changes in zircon rare earth elements (REE) concentrations, and indicate prograde heating, mineral breakdown, and melt generation took place from 62 to 53 Ma. TTi-zr increases from ~650 to 850 °C during this interval. A prominent gap in zircon ages is observed from 53 to 51 Ma and is interpreted to reflect the timing of peak metamorphism and melting, which caused zircon dissolution. The age gap is an inflection point in several geochemical-temporal trends that suggest crystallization and cooling dominated afterward, from 51 to 42 Ma. Supporting this interpretation is an increase in zircon U/Th and Hf, a decrease in TTi-zr, increasing zircon (Dy/Yb)n, and textural evidence for coupled dissolution–reprecipitation processes that resulted in zircon (re)crystallization. In addition, whole rock REE, large ion lithophile elements, and major elements suggest that the Pan Tak Granite experienced advanced fractional crystallization during this time. High-silica, muscovite± garnet leucogranite dikes that crosscut two-mica granite represent more evolved residual melt compositions. The Pan Tak Granite was formed by fluid-deficient melting and biotite dehydration melting of meta-igneous protoliths, including Jurassic arc rocks and the Proterozoic Oracle Granite. The most likely causes of melting are interpreted to be a combination of (1) radiogenic heating and relaxation of isotherms associated with crustal thickening under a plateau environment, (2) heat and fluid transfer related to the Laramide continental arc, and (3) shear and viscous heating related to the deformation of the deep lithosphere. The characteristics and petrologic processes that created the Pan Tak Granite are strikingly similar to intrusive suites in the Himalayan leucogranite belt and further support the association between the North American Cordilleran anatectic belt and a major orogenic and thermal event during the Laramide orogeny.

     
    more » « less
  5. Abstract

    The distribution of the short‐lived radionuclide26Al in the early solar system remains a major topic of investigation in planetary science. Thousands of analyses are now available but grossite‐bearing Ca‐, Al‐rich inclusions (CAIs) are underrepresented in the database. Recently found grossite‐bearing inclusions inCO3 chondrites provide an opportunity to address this matter. We determined the oxygen and magnesium isotopic compositions of individual phases of 10 grossite‐bearingCAIs in the Dominion Range (DOM) 08006 (CO3.0) andDOM08004 (CO3.1) chondrites. All minerals inDOM08006CAIs as well as hibonite, spinel, and pyroxene inDOM08004 are uniformly16O‐rich (Δ17O = −25 to −20‰) but grossite and melilite inDOM08004CAIs are not; Δ17O of grossite and melilite range from ~ −11 to ~0‰ and from ~ −23 up to ~0‰, respectively. Even within this small suite, in the two chondrites a bimodal distribution of the inferred initial26Al/27Al ratios (26Al/27Al)0is seen, with four having (26Al/27Al)0≤1.1 × 10−5and six having (26Al/27Al)0≥3.7 × 10−5. Five of the26Al‐richCAIs have (26Al/27Al)0within error of 4.5 × 10−5; these values can probably be considered indistinguishable from the “canonical” value of 5.2 × 10−5given the uncertainty in the relative sensitivity factor for grossite measured by secondary ion mass spectrometry. We infer that the26Al‐poorCAIs probably formed before the radionuclide was fully mixed into the solar nebula. All minerals in theDOM08006CAIs, as well as spinel, hibonite, and Al‐diopside in theDOM08004CAIs retained their initial oxygen isotopic compositions, indicating homogeneity of oxygen isotopic compositions in the nebular region where theCOgrossite‐bearingCAIs originated. Oxygen isotopic heterogeneity inCAIs fromDOM08004 resulted from exchange between the initially16O‐rich (Δ17O ~−24‰) melilite and grossite and16O‐poor (Δ17O ~0‰) fluid during hydrothermal alteration on theCOchondrite parent body; hibonite, spinel, and Al‐diopside avoided oxygen isotopic exchange during the alteration. Grossite and melilite that underwent oxygen isotopic exchange avoided redistribution of radiogenic26Mg and preserved undisturbed internal Al‐Mg isochrons. The Δ17O of the fluid can be inferred from O‐isotopic compositions of aqueously formed fayalite and magnetite that precipitated from the fluid on theCOparent asteroid. This and previous studies suggest that O‐isotope exchange during fluid–rock interaction affected mostCAIs in CO ≥3.1 chondrites.

     
    more » « less