Pre-eruptive conditions of the Hideaway Park topaz rhyolite: Insights into metal source and evolution of magma parental to the Henderson porphyry molybdenum deposit, Colorado

Celestine N. Mercer; Albert H. Hofstra; Todor I. Todorov; Julie Roberge; Alain Burgisser; David T. Adams; Michael A. Cosca

doi:10.1093/petrology/egv010

Home/Archive/USGS

USGS · 70162170

Pre-eruptive conditions of the Hideaway Park topaz rhyolite: Insights into metal source and evolution of magma parental to the Henderson porphyry molybdenum deposit, Colorado

Celestine N. Mercer·Albert H. Hofstra·Todor I. Todorov·Julie Roberge·Alain Burgisser·David T. Adams·Michael A. Cosca

View original source ↗Export metadata ↓

Abstract

The Hideaway Park tuff is the only preserved extrusive volcanic unit related to the Red Mountain intrusive complex, which produced the world-class Henderson porphyry Mo deposit. Located within the Colorado Mineral Belt, USA, Henderson is the second largest Climax-type Mo deposit in the world, and is therefore an excellent location to investigate magmatic processes leading to Climax-type Mo mineralization. We combine an extensive dataset of major element, volatile, and trace element abundances in quartz-hosted melt inclusions and pumice matrix glass with major element geochemistry from phenocrysts to reconstruct the pre-eruptive conditions and the source and evolution of metals within the magma. Melt inclusions are slightly peraluminous topaz rhyolitic in composition and are volatile-charged (≤6 wt % H 2 O, ≤600 ppm CO 2 , &sim;0·3–1·0 wt % F, &sim;2300–3500 ppm Cl) and metal-rich (&sim;7–24 ppm Mo, &sim;4–14 ppm W, &sim;21–52 ppm Pb, &sim;28–2700 ppm Zn, <0·1–29 ppm Cu, &sim;0·3–1·8 ppm Bi, &sim;40–760 ppb Ag, &sim;690–1400 ppm Mn). Melt inclusion and pumice matrix glass chemistry reveal that the Hideaway Park magma evolved by large degrees of fractional crystallization (≤60–70%) during quartz crystallization and melt inclusion entrapment at pressures of ≤300 MPa (≤8 km depth), with little to no crystallization upon shallow ascent and eruption. Filter pressing, crystal settling, magma recharge and mixing of less evolved rhyolite melt, and volatile exsolution were important processes during magma evolution; the low estimated viscosities (&sim;10 5 –10 10  Pa s) of these H 2 O- and F-rich melts probably enhanced these processes. A noteworthy discrepancy between the metal contents in the pumice matrix glass and in the melt inclusions suggests that after quartz crystallization ceased upon shallow magma ascent and eruption, the Hideaway Park magma exsolved an aqueous fluid into which Mo, Bi, Ag, Zn, Mn, Cs, and Y strongly partitioned. Given that the Henderson deposit contains anomalous abundances of not only Mo, but also W, Pb, Zn, Cu, Bi, Ag, and Mn, we suggest that these metals were sourced from similar fluids exsolved from unerupted portions of the same magmatic system. Trace element ratios imply that Mo was sourced deep, from either the lower crust or metasomatized mantle. The origin of sulfur remains unresolved; however, given the extremely low S solubility of rhyolite melts in the shallow crust we favor the possibility that another source of S might supplement or account for that present in the ore deposit, probably the comagmatic, mantle-derived lamprophyres that occur in minor quantities with the voluminous topaz rhyolites in the area. To account for the 437 Mt of MoS 2 (&sim;1·0 × 10 6  t Mo) present in the Henderson ore deposit, a volume of &sim;45 km 3 of Hideaway Park rhyolite magma would have been necessary to supply the Mo (a cylindrical pluton measuring 3·1 km × 6·0 km) along with sparging of &sim;6·8 × 10 5  t of S from &sim;0·05 km 3 of lamprophyre magma. Based on a weighted mean 40 Ar/ 39 Ar age of 27·58 ± 0·24 Ma, similar melt geochemistry, and characteristically F-rich biotite phenocrysts, we conclude that the Hideaway Park tuff was cogenetic with the intrusions at Red Mountain that formed the Henderson deposit.

Explore related subjects

Journal of Petrology Article

Source-reported bounding extent: 39.71986348549764° to 39.930800820752765° latitude; -105.94940185546875° to -105.72555541992188° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Celestine N. Mercer, Albert H. Hofstra, Todor I. Todorov, Julie Roberge, Alain Burgisser, David T. Adams, Michael A. Cosca. 2015-05-12. Pre-eruptive conditions of the Hideaway Park topaz rhyolite: Insights into metal source and evolution of magma parental to the Henderson porphyry molybdenum deposit, Colorado. https://doi.org/10.1093/petrology%2Fegv010

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

KEEP EXPLORING

Related USGS reports

USGS2026

Petrological, textural, and volatile constraints on the architecture of Mauna Loa’s magmatic system over 170 years (1852–2022)

Mauna Loa’s reawakening in 2022 highlighted the limited availability of petrological constraints characterizing the geometry and dynamics of pre-eruptive storage at this volcano. In particular, there has been relatively little prior work constraining pre-eruptive volatile budgets, mineral zoning and textures. To address this, we examine olivine- and orthopyroxene-hosted melt inclusions (MIs), their mineral hosts, and matrix glasses from seven Mauna Loa eruptions spanning 170 years (1852–2022). We present the first extensive MI dataset from Mauna Loa with direct Raman measurements of vapor bubble CO 2 . Olivine- and orthopyroxene-hosted MIs contain 0.3–0.5 wt % H 2 O and 100–600 μg/g CO 2 . MI saturation pressures reveal that the dominant reservoir is located at ~2–5 km depth, consistent with fluid inclusion barometry and geophysical constraints. This depth has remained consistent over the past 170 years, although some eruptions appear to have tapped a narrower range of depths than others. Even the high-Mg# orthopyroxene phenocryst-hosted MIs record saturation pressures of <150 MPa (<6 km depth), which is inconsistent with a deep sub-Moho origin proposed by previous studies that rely on thermodynamic modeling. Our Raman measurements also provide an opportunity to assess the accuracy of different methods to reconstruct the amount of CO 2 in vapor bubbles, with applications to MI datasets published without direct Raman measurements. In general, models that partition CO 2 between the melt and bubble using the bubble volume prior to pre-eruptive quenching (e.g. Melt Inclusion Modification Corrections, MIMiC) systematically underestimate CO 2 contents by a median of ~100–200 μg/g (depending on the model and K D choice) and thus underestimate storage depths by ~0.5–1.5 km. While these methods are useful for first-order estimates of bubble CO 2 in published datasets, new studies should measure bubble CO 2 by Raman spectroscopy or perform MI rehomogenization. The zoning and textures of the crystals that host the MIs provide additional constraints on the plumbing system to supplement MI barometry. Orthopyroxene crystals exhibit varied zoning patterns, including oscillatory zoning that records repeated fluctuations in magmatic conditions, and sector zoning indicating rapid, kinetically controlled crystal growth. While faster Fe-Mg diffusion in olivine results in a simpler zoning record, complexly zoned olivines with higher-Fo ‘shoulders’ from 1852, 1868, and 1950 also record the injection of more primitive melts at shallow levels. Electron backscatter diffraction reveals that touching crystals in orthopyroxene glomerocrysts exhibit only very weakly systematic crystallographic orientation relationships, perhaps due to rapid aggregation in a dynamic system. Overall, our data are consistent with a large, shallow (~2–5 km) pre-eruptive reservoir repeatedly recharged by primitive magma injections that homogenized prior to eruption due to their small volume relative to the overall size of this reservoir. This chemical complexity is effectively recorded in mineral and MI compositions (orthopyroxene zoning in particular) but is obscured in glass and whole-rock analyses.

Hawaii↗

USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗

USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗