skip to main content


Title: Puhahonu: Earth’s biggest and hottest shield volcano
New bathymetric and gravity mapping, refined volume calculations and petrologic analyses show that the Hawaiian volcano P¯uh¯ahonu is the largest and hottest shield volcano on Earth. This ∼12.5-14.1 Ma volcano in the northwest Hawaiian Ridge (NWHR) is twice the size of Mauna Loa volcano (148 ±29 vs. 74.0 ×103km3), which was assumed to be not only the largest Hawaiian volcano but also the largest known shield volcano. We considered four testable mechanisms to increase magma production, including 1) thinner lithosphere, 2) slower propagation rate, 3) more fertile source, and 4) hotter mantle. The first three of these have been ruled out. The lithosphere was old (∼88 Myrs) when P¯uh¯ahonu was formed, and thus, too thick and cold to allow for greater extents of partial melting. The propagation rate was relatively fast when it erupted (87 km/Myr), so this is another unlikely reason. Source fertility was Kea-like and no more fertile than for other much smaller NWHR volcanoes. A hotter mantle remains the best mechanism to produce the large magma volumes and is consistent with the high forsteritic olivine phenocryst compositions (up to 91.8%) and the calculated high percent of melting (24%). Thus, the gargantuan size of P¯uh¯ahonu reflects its high melting temperature, the highest reported for any Cenozoic basalt. A solitary wave within the Hawaiian plume is the probable cause of P¯uh¯ahonu’s higher melting temperature and the resulting increased volume flux given the absence of a more fertile source for P¯uh¯ahonu basalts, as found for many basalts from the Hawaiian Islands.  more » « less
Award ID(s):
1834758
NSF-PAR ID:
10283332
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Chemical geology
Volume:
542
ISSN:
0168-9622
Page Range / eLocation ID:
116296
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Oceanic island basalts are targeted for geochemical study because they provide a direct window into mantle composition and a wealth of information on the dynamics and timescales associated with Earth mixing. Previous studies mainly focused on the shield volcanic stage of oceanic islands and the more fusible, enriched mantle components that are easily distinguished in those basalts. Mantle depleted compositions are typically more difficult to resolve unless large amounts of this material participated in mantle melting (e.g., mid-ocean ridges), or unique processes allow for their compositions to be erupted undiluted, such as very small degrees of melting of a source with minimal fusible enriched components (e.g., rejuvenated basalts) or as xenoliths (e.g., abyssal peridotites). Mantle depleted components, defined here as material with low time-integrated Rb/Sr (low 87Sr/86Sr) and high time-integrated Sm/Nd and Lu/Hf ratios (high 143Nd/144Nd and 176Hf/177Hf) relative to primitive mantle, derive from a potentially very large volume reservoir (up to 80% of the mantle), and therefore need adequate characterization in order estimate the composition of the Earth and mantle-derived melts. This review focuses on mantle depleted compositions in oceanic island basalts using the Hawaiian-Emperor chain as a case study. The Hawaiian-Emperor chain is the ∼6000 km long geological record of the deeply sourced Hawaiian mantle plume, active for>81 Myr. Hawaiian volcanism evolves through four volcanic stages as a volcano traverses the Hawaiian plume: alkalic preshield, tholeiitic shield (80–90% volcano volume), alkalic postshield (∼1%), and silica undersaturated rejuvenated (< 0.1%). We report Pb-Sr-Nd-Hf isotope compositions and trace element concentrations of three rejuvenated Northwest Hawaiian Ridge basalts and compare them to an exhaustive compiled dataset of basalts from the Hawaiian Islands to the Emperor Seamounts. The Northwest Hawaiian Ridge (NWHR) includes 51 volcanoes spanning ∼42 m.y. between the bend in the Hawaiian-Emperor chain and the Hawaiian Islands where there is no high-precision isotopic data published on the rejuvenated-stage over ∼47% of the chain. NWHR and Hawaiian Island rejuvenated basalts are geochemically similar, indicating a consistent source for rejuvenated volcanism over ∼12.5 million years. In contrast, shield-stage basalts from the oldest Emperor Seamounts are more depleted in isotopic composition (i.e., higher 176Hf/177Hf, and 143Nd/144Nd with lower 87Sr/86Sr and 208Pb*/206Pb*) and trace element concentrations (i.e., much lower concentrations of highly incompatible elements) than all other depleted Hawaiian basalts younger than the bend, including NWHR rejuvenated basalts. The strongly depleted source for the oldest Emperor Seamounts (> 70 Ma) was likely related to interaction with the Kula-Pacific-Izanagi mid-ocean ridge spreading system active near the Hawaiian plume in the Late Cretaceous. In contrast, the incompatible trace element ratios of NWHR rejuvenated basalts require a distinct source in the Hawaiian mantle plume that was imprinted by ancient (> 1 Ga) partial melting, likely ancient recycled oceanic lithosphere. This review of the geochemistry of Hawaiian depleted components documents the need for the sampling of multiple distinctive depleted compositions, each preferentially melted during specific periods of Hawaiian plume activity. This suggests that the composition of depleted components can evolve during the lifetime of the mantle plume, as observed for enriched components in the Hawaiian mantle plume. Changes in the composition of depleted components are dominantly controlled by the upper mantle tectonic configurations at the time of eruption (i.e., proximity to a mid-ocean ridge), as this effect overwhelms the signal imparted by potentially sampling different lower mantle components through time. 
    more » « less
  2. Abstract

    The break-up of supercontinents is often temporally and spatially associated with large outpourings of basaltic magmas in the form of large igneous provinces (LIPs) and seaward dipping reflectors (SDRs). A widespread view is that the upwelling of hot mantle plumes drives both continental break-up and generation of associated LIPs. This is supported by petrologic estimates of the temperature from olivine-melt thermometers applied to basaltic magmas. These thermometers must be applied to a primary mantle-derived magma, requiring the selection of an appropriate primitive magma and an assumption of how much olivine is to be back-added to correct for fractional crystallization. We evaluated the effects of these assumptions on formation temperatures by compiling and analyzing a database of North Atlantic igneous province (NAIP) and Central Atlantic magmatic province (CAMP) lavas and olivines. Ni and FeOT systematics suggest that many picrite magmas have undergone olivine addition and are not true liquids, requiring careful selection of primitive magmas. The maximum amount of back-added olivine was determined by constraining mantle peridotite melt fractions for a range of possible mantle potential temperatures and continental lithosphere thicknesses. Using an empirical relationship between melting degree and forsterite (Fo) content, we show that the possible maximum olivine forsterite content in equilibrium with NAIP magmas is 90.9, which is lower than the maximum olivine forsterite content observed in the NAIP olivine population. We infer primary magmas that lead to mantle potential temperatures of 1420°C for the NAIP and 1330°C for CAMP. Using a similar approach for consistency, we estimate a mantle potential temperature of 1350°C for mid-ocean ridge basalts (MORB). Our results suggest that LIPs associated with continental break-up are not significantly hotter than MORB, which suggests that continental break-up may not be driven by deep-seated thermal plumes. Instead, we suggest that such voluminous magmatism might be related to preferential melting of fertile components within the lithosphere triggered by far-field extensional stresses.

     
    more » « less
  3. Abstract

    The timing and rate of decompression melting of a compositionally heterogeneous mantle during continental rifting are assessed with a new one‐dimensional geodynamic code, MELT1D. MELT1D computes pressure and temperature in the extending lithosphere and rising asthenosphere and calculates the resulting melt fraction and eruption rates for different lithologies. A series of models simulate syn‐rift melt production from (a) dry and wet depleted lherzolite similar to the mantle that underlies most mid‐ocean ridges, (b) dry and wet relatively fertile ultramafic compositions representing plume or primitive mantle material, (c) pyroxenite representing recycled ultramafic oceanic crust or magmatic metasomes, and (d) basalt representing recycled mafic crust or metasomes. The models predict sequential melting of the different compositions that is broadly consistent with basalt eruption histories in many Phanerozoic rifts. Results show a progressive transition in magma sources as the lithosphere thins, beginning with melting of wet mantle and compositionally fertile mafic components near the lithosphere‐asthenosphere boundary during the earliest stages of extension. This transitions to magmatism dominated by melting of relatively fertile ultramafic components (pyrolitic and pyroxenitic compositions) as extension progresses, and finally to melting of ambient lherzolite asthenosphere as lithosphere thinning approaches breakup. Mantle composition, pre‐rift lithosphere thickness, and mantle temperature exert the greatest controls on the timing and volumes of magmas produced from each lithology. In general, a cool or thick lithosphere has a greater capacity to sequester fertile lithologies than thin or warm lithosphere, and thus has a greater capacity to produce early syn‐rift magmas without requiring a hot mantle plume.

     
    more » « less
  4. Abstract

    Many lines of evidence from high P–T experiments, thermodynamic models, and natural observations suggest that slab-derived aqueous fluids, which flux mantle wedges contain variable amounts of dissolved carbon. However, constraints on the effects of H2O–CO2 fluids on mantle melting, particularly at mantle wedge P–T conditions, are limited. Here, we present new piston cylinder experiments on fertile and depleted peridotite compositions with 3.5 wt.% H2O and XCO2 [= molar CO2 / (CO2 + H2O)] of 0.04–0.17. Experiments were performed at 2–3 GPa and 1350°C to assess how temperature, peridotite fertility, and XCO2 of slab-derived fluid affects partial melting in mantle wedges. All experiments produce olivine + orthopyroxene +7 to 41 wt.% partial melt. Our new data, along with previous lower temperature data, show that as mantle wedge temperature increases, primary melts become richer in SiO2, FeO*, and MgO and poorer CaO, Al2O3, and alkalis when influenced by H2O–CO2 fluids. At constant P–T and bulk H2O content, the extent of melting in the mantle wedge is largely controlled by peridotite fertility and XCO2 of slab-fluid. High XCO2 depleted compositions generate ~7 wt.% melt, whereas, at identical P–T, low XCO2 fertile compositions generate ~30 to 40 wt.% melt. Additionally, peridotite fertility and XCO2 have significant effects on peridotite partial melt compositions. At a constant P–T–XCO2, fertile peridotites generate melts richer in CaO and Al2O3 and poorer in SiO2, MgO + FeO, and alkalis. Similar to previous experimental studies, at a constant P–T fertility condition, as XCO2 increases, SiO2 and CaO of melts systematically decrease and increase, respectively. Such distinctive effects of oxidized form of dissolved carbon on peridotite partial melt compositions are not observed if the carbon-bearing fluid is reduced, such as CH4-bearing. Considering the large effect of XCO2 on melt SiO2 and CaO concentrations and the relatively oxidized nature of arc magmas, we compare the SiO2/CaO of our experimental melts and melts from previous peridotite + H2O ± CO2 studies to the SiO2/CaO systematics of primitive arc basalts and ultra-calcic, silica-undersaturated arc melt inclusions. From this comparison, we demonstrate that across most P–T–fertility conditions predicted for mantle wedges, partial melts from bulk compositions with XCO2 ≥ 0.11 have lower SiO2/CaO than all primitive arc melts found globally, even when correcting for olivine fractionation, whereas partial melts from bulk compositions with XCO2 = 0.04 overlap the lower end of the SiO2/CaO field defined by natural data. These results suggest that the upper XCO2 limit of slab-fluids influencing primary arc magma formation is 0.04 < XCO2 < 0.11, and this upper limit is likely to apply globally. Lastly, we show that the anomalous SiO2/CaO and CaO/Al2O3 signatures observed in ultra-calcic arc melt inclusions can be reproduced by partial melting of either CO2-bearing hydrous fertile and depleted peridotites with 0 < XCO2 < 0.11 at 2–3 GPa, or from nominally CO2-free hydrous fertile peridotites at P > 3 GPa.

     
    more » « less
  5. SUMMARY

    Along age-progressive hotspot volcano chains, the emplacement rate of igneous material varies through time. Time-series analysis of changing emplacement rates at a range of hotspots finds that these rates vary regularly at periods of a few to several tens of millions of years, indicative of changing melt production within underlying mantle plumes. Many hotspots exhibit at least one period between ∼2 and 10 Myr, consistent with several proposed mechanisms for changing near-surface plume flux, and thus melting rate, such as small-scale convection, solitary waves and instability formation in tilted plume conduits. Here, we focus on quantifying instability growth within plumes tilted by overlying plate motion. Previous studies using fluids with constant or temperature-dependent viscosity suggest that such instabilities should not form under mantle conditions. To test this assertion, we use a modified version of the finite element code ASPECT to simulate 400 Myr of evolution of a whole-depth mantle plume rising through the transition zone and spreading beneath a moving plate. In a 2-D spherical shell geometry, ASPECT solves the conservation equations for a compressible mantle with a thermodynamically consistent treatment of phase changes in the mantle transition zone and subject to either a temperature- and depth-dependent linear rheology or a temperature-, depth- and strain-rate dependent non-linear rheology. Additionally, we examine plume evolution in a mantle subject to a range of Clapeyron slopes for the 410 km (1–4 MPa K–1) phase transitions. Results suggest that plume conduits tilted by >67° become unstable and develop instabilities that lead to initial pulses in the transition zone followed by repeated plume pulsing in the uppermost mantle. In these cases, pulse size and frequency depend strongly on the viscosity ratio between the plume and ambient upper mantle. Based upon our results and comparison with other studies, we find that the range of statistically significant periods of plume pulsing in our models (∼2–7 Myr), the predicted increase in melt flux due to each pulse (3.8–26 × 10−5 km3 km−1 yr−1), and the time estimated for a plume to tilt beyond 67° in the upper mantle (10–50 Myr) are consistent with observations at numerous hotspot tracks across the globe. We suggest that pulsing due to destabilization of tilted plume conduits may be one of several mechanisms responsible for modulating the melting rate of mantle plumes as they spread beneath the moving lithosphere.

     
    more » « less