Yellowstone Hotspot Research Articles

Plume-modified mantle flow in the northern Basin and Range and southern Cascadia back-arc region since ca. 12 Ma

AbstractBimodal volcanism and rhyolite migration along the High Lava Plains in central Oregon (United States) lie above a broader feature defined by low seismic velocity in the upper mantle that emanates from the Yellowstone hotspot (northwest United States) and extends westward across the northern Basin and Range. It was emplaced by a westward current, driven in part by rapid buoyancy-driven flow across the east-west cratonic boundary of North America. Geothermometry studies and geochemical considerations suggest that the low-velocity feature may be composed of moderately hot, low-density mantle derived from the Yellowstone plume but diluted by thermomechanical erosion and entrainment of colder mantle lithosphere. Finger-like conduits of plume-modified mantle beneath Quaternary eruption sites delineate flow-line channels that have developed across the broader mantle structure since 2 Ma. These channels have allowed low-density mantle to accumulate against the Cascades arc, thus providing a heated mantle source for mafic magmatism in the Newberry (Oregon) and Medicine Lake (California) volcanic fields.

40Ar/39Ar and paleomagnetic constraints on the age and areal extent of the Picabo volcanic field: Implications for the Yellowstone hotspot

AbstractThe Picabo volcanic field is one of the key silicic volcanic fields in the time-transgressive track of the Yellowstone hotspot. The Picabo volcanic field is also one of the most poorly defined volcanic fields along the track of the Yellowstone hotspot. Determining the age and areal extent of the Picabo volcanic field ignimbrites is one of the primary objectives of this study. In our effort to correlate ignimbrites within the Picabo volcanic field as well as identify those from the neighboring Twin Falls and Heise volcanic fields, we present new petrographic, 40Ar/39Ar, and paleomagnetic data. With these data, we correlated several ignimbrites within the Picabo volcanic field. In some cases, we correlate units previously thought to be in the Picabo volcanic field to older volcanic fields. This includes the Picabo Tuff, which we suggest originates from the Twin Falls volcanic field rather from its namesake volcanic field. The first and best documented major silicic eruption of the volcanic field, the Arbon Valley Tuff, is also the largest ignimbrite in the Picabo volcanic field. There is disagreement as to whether the Arbon Valley Tuff is the result of a single ignimbrite eruption or multiple eruptions. We previously have suggested that the Arbon Valley Tuff is the result of two eruptions, one at 10.41 ± 0.01 Ma and the other at 10.22 ± 0.01 Ma (Anders et al., 2014). Those combining radiometric dates into a single eruption age report ages of 10.2 Ma, 10.27 ± 0.01 Ma, 10.34 ± 0.03 Ma, and 10.44 ± 0.27 Ma. We also suggest the final eruption of the Picabo volcanic field was the tuff of American Falls dated at 7.58 ± 0.02 Ma. Estimates of the location of Picabo volcanic field have been used to mark a major change in the migration rate of the Yellowstone–Snake River Plain silicic volcanic system. Based on our new data, we found only minor changes of the boundaries of the Picabo volcanic field from previous studies. Using the age of the Arbon Valley Tuff (10.41 Ma), we calculated an extension-corrected migration rate of 2.27 ± 0.2 cm/yr between the position of the Picabo volcanic field and that of the Yellowstone volcanic field over the past ∼10 m.y. This estimate is close to the extension corrected 2.38 ± 0.21 cm/yr value based on the migration of the hotspot deformation field. These rates are consistent with independent estimates of North American plate velocity over the past 10 m.y. and therefore consistent with a fixed reference frame for the Yellowstone hotspot. These results stand in contrast with several recent models for the evolution of the Yellowstone–Snake River Plain volcanic system.We also discovered a new ignimbrite from the Heise volcanic field, the 4.37 ± 0.08 Ma tuff of Birch Creek Sinks, in core from the U.S. Geological Survey (USGS) borehole 2-2A, which now represents the youngest outflow ignimbrite of the Heise volcanic field. Although recently, several intracaldera ignimbrites younger than 4 Ma have been identified in the volcanic field, the age range of outflow ignimbrites from the Heise volcanic field is now extended from 6.66 Ma to at least 4.37 Ma.

Elemental abundance patterns and Sr-, Nd- and Hf-isotope systematics for the Yellowstone hotspot and Columbia River flood basalts: Bearing on petrogenesis

The Yellowstone hotspot track in the northwestern US has a 17-myr history that began with eruption of the voluminous Columbia River basalts (CRB), and produced an eastward track with the bimodal rhyolite and basaltic eruptive centers of the Snake River Plain (SRP) as the North American Plate drifted westwards. We have generated the first Hf isotopic data on basalts in this large igneous province along with bulk-rock trace-element abundance patterns and Sr-Nd isotopic ratios to assess the hypothesized kinship with co-located rhyolites throughout the hotspot track. The isotopic compositions for all three elements span much of the compositional range observed in mantle-derived basalts worldwide. Varying between a high of +11.3 in parts of the CRB and a low of −8.3 in some of the late SRP basalts, the new εHf values are geographically controlled, similar to the shift first documented with Sr and Nd isotopes across the “0.706 Sr line,” which marks the boundary between thin- accreted young ocean lithosphere and arc terranes in the west to thick, old cratonized crust in the east. This observation shows that there is significant lithospheric control on the isotopic compositions of both basalts and rhyolites associated with the Yellowstone hotspot track.Modeling of assimilation and fractional crystallization (AFC) processes in seven basaltic lava flows filling the 2-Ma Big Bend Ridge Caldera at Mesa Falls, Idaho, near the eastern end of the hotspot track suggests that the crustal contribution to the Nd and Hf concentration budgets increases systematically with time and was significant. A comparison of new and previously published bulk-rock major-oxide concentrations (e.g., the alkalis Na2O and K2O) with our previously published data for H2O, Cl, F and S in olivine-hosted melt inclusions suggests that the volatiles are sub-crustal, possibly from the lithospheric mantle, and degas extensively during magma ascent and crystallization.

Understanding the isotopic and chemical evolution of Yellowstone hot spot magmatism using magmatic-thermomechanical modeling

Large-scale melting and reworking of the crust and the upper mantle occurs in areas of mantle plume-crust interaction, generating voluminous and bimodal basaltic and silicic magmatism. Continental hot spot track magmas preserve geochemical and temporal records of these processes. Here, we apply and further develop the large scale (crust and the upper mantle) I2VIS magmatic-thermomechanical modeling code to investigate the potential origins of chemical and isotopic trends in the Yellowstone hot spot track, perhaps the best studied and constrained subcontinental mantle plume system on Earth. We model the propagation of melt via dikes in solid crust and via percolation in partially molten crust by using Lagrangian markers which also track the chemical and isotopic compositions of the melts that eventually erupt at the surface. Confirming the results of earlier geophysical and geochemical studies, we show that the eruptive activity at Yellowstone is best explained in terms of the formation of a ~15 km-thick mid-crustal mafic sill complex with its top at a depth of ~8 km that develops over a period of ~3 Myr. This sill complex releases rhyolitic fractionates and creates additional rhyolites by melting the surrounding crust, driving the voluminous rhyolitic volcanism which characterizes the Yellowstone system.We recognize three key trends in the evolution of a continental hot spot volcanic center such as Yellowstone: (1) frequent, small eruptions results in in larger erupted rhyolite volumes compared to rare and large eruptions, (2) erupted magmas tend to be initially comprised of a large degree of crustal melt, and the relative contribution of fractionates of basalts from the mantle increases with time, and (3), the initial depth of the crust which melts to produce eruptible rhyolites becomes shallower with time. The first of these trends is simply a function of the continuous supply of new basalt to the crust from the mantle. The latter trend, which is responsible for low-δ18O rhyolites, is produced by a combination of progressive melting of shallower crust as the system becomes hotter, repeated caldera collapses advecting shallow crust to a depth where it can melt, and the overplating and burial of shallow crust by repeated intrusions of basalt. The mid- and lower crust is not a significant source of erupted rhyolitic melts, as both the crustal melting and basalt fractionation takes place in the sill complex which forms near the bottom of the upper crust. The resulting modeled isotopic evolution of erupted magmas are a good match with the actual stable and radiogenic isotopic record of Yellowstone hot spot track volcanism preserved in zircon phenocrysts, and further replicates the results of recent geophysical imaging campaigns, giving us confidence that the assumptions and results of these thermomechanical models are broadly correct.

Using Fast Marching Eikonal Solver to Compute 3‐D Pds Traveltime for Deep Receiver‐Function Imaging

AbstractIn common‐conversion‐point stacking of receiver functions, most current studies compute 3‐D relative Pds traveltime corrections by integrating traveltime anomalies along 1‐D raypaths. This ray tracing approach is generally time‐consuming and less accurate when prominent velocity anomalies exist and effects of 3‐D raypaths become significant. In this study we introduce a new scheme that utilizes a fast‐marching method eikonal solver to improve both the efficiency and accuracy of 3‐D Pds traveltime computation. We first employ a 1‐D raytracing method and the iasp91 model to calibrate the accuracy of the new scheme and optimize the parameters of the numerical solver. We then apply the new scheme to compute a massive number of Pds traveltimes using two sets of 3‐D synthetic models, one set with a high‐velocity slab and another set with a low‐velocity plume, and compare these 3‐D traveltimes with those computed with the raypath integrating approach. We find 2.7% and 11.8% raypaths in the two slab models and 7.8% and 12.0% raypaths in the two plume models show a 3‐D traveltime difference of ≥0.5 s. We apply the proposed scheme to a subset of transportable array receiver functions that sample the transition zone beneath the Yellowstone hotspot and find that a common‐conversion‐point stack using 3‐D Pds traveltimes computed by the eikonal solver method has the best focused P660s. Finally, we illustrate that computational times can be reduced by 1 to 2 orders of magnitude with the new scheme to compute 3‐D Pds traveltimes of 20,000–200,000 receiver functions.

Caldera Life-Cycles of the Yellowstone Hotspot Track: Death and Rebirth of the Heise Caldera

Caldera Life-Cycles of the Yellowstone Hotspot Track: Death and Rebirth of the Heise Caldera

Anomalous mantle transition zone beneath the Yellowstone hotspot track

The origin of the Yellowstone and Snake River Plain volcanism has been strongly debated. The mantle plume model successfully explains the age-progressive volcanic track, but a deep plume structure has been absent in seismic imaging. Here I apply diffractional tomography to receiver functions recorded at USArray stations to map high-resolution topography of mantle transition-zone discontinuities. The images reveal a trail of anomalies that closely follow the surface hotspot track and correlate well with a seismic wave-speed gap in the subducting Farallon slab. This observation contradicts the plume model, which requires anomalies in the mid mantle to be confined in a narrow region directly beneath the present-day Yellowstone caldera. I propose an alternative interpretation of the Yellowstone volcanism. About 16 million years ago, a section of young slab that had broken off from a subducted spreading centre in the mantle first penetrated the 660 km discontinuity beneath Oregon and Idaho, and pulled down older stagnant slab. Slab tearing occurred along pre-existing fracture zones and propagated northeastward. This reversed-polarity subduction generated passive upwellings from the lower mantle, which ascended through a water-rich mantle transition zone to produce melting and age-progressive volcanism. The mantle transition zone in the western United States is perturbed along a path that mirrors the line of the Yellowstone hotspot track at the surface, according to analysis of tomographic data.

Lower-mantle plume beneath the Yellowstone hotspot revealed by core waves

The Yellowstone hotspot, located in North America, is an intraplate source of magmatism the cause of which is hotly debated. Some argue that a deep mantle plume sourced at the base of the mantle supplies the heat beneath Yellowstone, whereas others claim shallower subduction or lithospheric-related processes can explain the anomalous magmatism. Here we present a shear wave tomography model for the deep mantle beneath the western United States that was made using the travel times of core waves recorded by the dense USArray seismic network. The model reveals a single narrow, cylindrically shaped slow anomaly, approximately 350 km in diameter that we interpret as a whole-mantle plume. The anomaly is tilted to the northeast and extends from the core–mantle boundary to the surficial position of the Yellowstone hotspot. The structure gradually decreases in strength from the deepest mantle towards the surface and if it is purely a thermal anomaly this implies an initial excess temperature of 650 to 850 °C. Our results strongly support a deep origin for the Yellowstone hotspot, and also provide evidence for the existence of thin thermal mantle plumes that are currently beyond the resolution of global tomography models. The Yellowstone hotspot could be fed by a thin, thermal mantle plume that extends from the core–mantle boundary to the surface position of the hotspot, according to analyses of seismic data.

Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology, oxygen isotope, and hafnium isotope analyses of zircon

We present new high-precision CA-ID-TIMS and in situ U–Pb ages together with Hf and O isotopic analyses (analyses performed all on the same grains) from four tuffs from the 15−10 Ma Bruneau–Jarbidge center of the Snake River Plain and from three rhyolitic units from the Kimberly borehole in the neighboring 10−6 Ma Twin Falls volcanic center. We find significant intrasample diversity in zircon ages (ranges of up to 3 Myr) and in δ18O (ranges of up to 6‰) and eHf (ranges of up to 24 e units) values. Zircon rims are also more homogeneous than the associated cores, and we show that zircon rim growth occurs faster than the resolution of in situ dating techniques. CA-ID-TIMS dating of a subset of zircon grains from the Twin Falls samples reveals complex crystallization histories spanning 104–106 years prior to some eruptions, suggesting that magma genesis was characterized by the cyclic remelting of buried volcanic rocks and intrusions associated with previous magmatic episodes. Age-dependent trends in zircon isotopic compositions show that rhyolite production in the Yellowstone hotspot track is driven by the mixing of mantle-derived melts (normal δ18O and eHf) and a combination of Precambrian basement rock (normal δ18O and eHf down to − 60) and shallow Mesozoic and Cenozoic age rocks, some of which are hydrothermally altered (to low δ18O values) by earlier stages of Snake River Plain magmatism. These crustal melts hybridize with juvenile basalts and rhyolites to produce the erupted rhyolites. We also observe that the Precambrian basement rock is only an important component in the erupted magmas in the first eruption at each caldera center, suggesting that the accumulation of new intrusions quickly builds an upper crustal intrusive body which is isolated from the Precambrian basement and evolves towards more isotopically juvenile and lower-δ18O compositions over time.

The leading wisps of Yellowstone: Post–ca. 5 Ma extension-related magmatism in the upper Wind River Basin, Wyoming (USA), associated with the Yellowstone hotspot tectonic parabola

The leading wisps of Yellowstone: Post–ca. 5 Ma extension-related magmatism in the upper Wind River Basin, Wyoming (USA), associated with the Yellowstone hotspot tectonic parabola

Spatial autocorrelation of Neogene-Quaternary lava along the Snake River Plain, Idaho, USA

The sequence of eruption, spatial pattern, and spatio-temporal relationships among the Neogene-Quaternary rhyolitic and basaltic lava along the Snake River Plain (SRP) in Idaho are analyzed applying the spatial methods of global and local Moran’s I, standard deviational ellipse, and Ripley’s K-function. The results of the analyses by the Moran’s I and K-function methods indicate a higher spatial autocorrelation, hence clustering, of rhyolitic lava compared to the more dispersed basaltic lava in each center of eruption along the SRP. The clustered nature of rhyolitic lava around each caldera either reflects the original spread and large thickness of the rhyolitic lava, or the absence of younger cover strata or lava like the distribution of rhyolite in the present caldera at the Yellowstone National Park. The standard deviational ellipses (SDEs) of the lavas indicate that younger basaltic lava that erupted from newer calderas overlapped older rhyolitic and basaltic lava as the position of the Yellowstone hotspot progressively migrated to the northeast along the SRP. The less eccentric SDEs of rhyolitic lava in each caldera probably reflect the original caldera-scale spread of viscous felsic lava, compared to the more eccentric and larger SDEs of basaltic lava which represent basalt’s wider and more directed spread due to its higher fluidity and ability to flow longer distances along the trend of the SRP. The alignment of the long axes of the lava SDEs with the trend of the Eastern SRP and the trend of systematic spatial overlap of older lava by successively younger basaltic lava corroborate the previously reported migration of the centers of eruption along the ESRP as the Yellowstone hotspot migrated to the northeast.

Discovery of Paleogene Deposits of the Central High Plains Aquifer In the Western Great Plains, U.S.A.

Abstract: Recent drilling in southwestern Kansas has produced intact cores to investigate the subsurface lithostratigraphy, sedimentary provenance, and chronostratigraphy of strata constituting the central High Plains Aquifer (HPA). In large portions of the HPA, groundwater withdrawals greatly exceed rates of recharge, leading to dramatic water-level declines and growing concerns for long-term sustainability. Two cores, HP1A (98 m) and CMC1 (50 m), are the first and deepest intact cores of the HPA ever collected. The cores show decameter-scale intercalations between suspended-load fluvial deposits composed of fine-grained sands with pedogenically modified overbank deposits, and very coarse-grained sands and gravels consistent with high-energy, bedload-dominated fluvial systems. Six intervals in HP1A and one from CMC1 were analyzed for detrital-zircon U-Pb ages by LA-ICP-MS. The HP1A samples show maximum depositional ages (MDAs) ranging from ∼ 27.5 to 36.4 Ma with depth. A MDA of ∼ 30 Ma was measured in the CMC1 core. These Paleogene zircons likely originated from explosive volcanism associated with the mid-Cenozoic ignimbrite flare-up (44–18 Ma), which blanketed much of western North America with high-volume air-fall tuffs. We propose that a large evaporite-dissolution basin in southwestern Kansas provided the accommodation space to preserve a record of Paleogene strata. The lack of younger middle–late Miocene zircons from cores in southwestern Kansas is striking given that such grains, likely derived from the Snake River Plain–Yellowstone hotspot volcanic provinces (16.1–0.6 Ma), are readily identified in the Ogallala Formation in north-central Kansas and Nebraska. The MDAs suggest Eocene to Oligocene age deposits that are time-equivalent to the Arikaree and White River groups in Nebraska are also present in Kansas.

Geology and 40Ar/39Ar geochronology of the middle Miocene McDermitt volcanic field, Oregon and Nevada: Silicic volcanism associated with propagating flood basalt dikes at initiation of the Yellowstone hotspot

Geology and 40Ar/39Ar geochronology of the middle Miocene McDermitt volcanic field, Oregon and Nevada: Silicic volcanism associated with propagating flood basalt dikes at initiation of the Yellowstone hotspot

Dual stable isotopes of CH4 from Yellowstone hot-springs suggest hydrothermal processes involving magmatic CO2

Volcanism and post-magmatism contribute significant annual methane (CH4) fluxes to the atmosphere (on par with other natural sources such as forest fire and wild animal emissions) and have been implicated in past climate-change events. The Yellowstone hot spot is one of the largest volcanic systems on Earth and is known to emit CH4 (as well as carbon dioxide (CO2) and other gases), but the ultimate sources of this CH4 flux have not been elucidated. Here we use dual stable isotope analysis (δ2H and δ13C) of CH4 sampled from ten high-temperature geothermal pools in Yellowstone National Park along with other isotopic and gas analyses to evaluate potential sources of methane. The average δ13C and δ2H values of CH4 emitted from hot springs (26.7 (±2.4) and −236.9 (±12.0) ‰, respectively) are inconsistent with microbial methanogenesis but do not allow distinction between thermogenic and abiotic sources. Correlation between δ13CCH4 and δ13C of dissolved inorganic C (DIC) is consistent with DIC as the parent C source for the observed CH4, or with equilibration of CH4 and DIC. Methane formation temperatures estimated by isotopic geothermometry based on δ13CCH4 and δ13CCO2 ranged from ~250–350°C, which is just below previous temperature estimates for the hydrothermal reservoir. Further, the δ2HH2O of the thermal springs and the measured δ2HCH4 values are consistent with equilibration between the source water and the CH4 at the formation temperatures. Though the ultimate origin of the CH4 could be attributed to either abiotic of themorgenic processes with subsequent isotopic equilibration, the C1/C2+ composition of the gases is more consistent with abiotic origins for most of the samples. Thus, our data support the hypothesis that subsurface rock-water interactions are responsible for at least a significant fraction of the CH4 flux from the Yellowstone National Park volcanic system.

The role of mantle‐derived magmas in the isotopic evolution of Yellowstone's magmatic system

AbstractInjection of mantle‐derived magmas into the Earth's crust provides the heat necessary to develop and maintain large silicic magmatic systems. However, the role of mantle‐derived magmas in controlling the compositional evolution of large silicic systems remains poorly understood. Here we examine the role of mantle‐derived magmas in the postcaldera magmatic system at Yellowstone Plateau, the youngest magmatism associated with the Yellowstone hotspot. Using microbeam techniques, we characterize the age and Hf isotope composition of single zircon crystals hosted in rhyolites from the most recent eruptive episode at Yellowstone Plateau, which produced the Central Plateau Member rhyolites. We place these zircon data into context by comparing them to new solution Hf isotope data for the Central Plateau Member glasses, Yellowstone basalts, and potential local crustal sources. Zircons in the Central Plateau Member rhyolites record a wide range of Hf isotope compositions relative to their host melts and extend from values similar to previously erupted Yellowstone rhyolites to values similar to Yellowstone basalts. Most zircons (∼90%) are in isotopic equilibrium with their host melt, but a significant proportion show εHf values higher than their host melt, thus providing the direct evidence that silicic derivatives of mantle‐derived basalts have recharged Yellowstone's magmatic system. Mixing models confirm that the isotopic characteristics of the youngest Yellowstone rhyolites can be explained by recharge of Yellowstone's magma reservoir with silicic derivatives of underplating, mantle‐derived basalts (∼5–10% material added by mass). This process helps drive the long‐term isotopic evolution of Yellowstone's magmatic system.

A crustal-scale integrated geophysical and tectonic study of the Snake River Plain region, northwestern U.S.A.

ABSTRACTThe Pacific Northwest region of North America is a site of very complex tectonomagmatic activity. This activity is due to subduction of the Pacific plate, the associated Cascade chain of volcanoes, micro-plate interactions, and mantle plume activity to the east of the plate margin that produced the Yellowstone hotspot track along the Eastern Snake River Plain (ESRP). A number of recent geophysical and geological studies have produced new results that have drawn attention to the complex tectonic setting of the region east of the Cascade Range, and its tectonic evolution is the subject of considerable scientific interest and debate. Numerous seismic studies have specifically focused on the crustal and upper mantle structure of the ESRP and Yellowstone area. However, crustal-scale studies of the Western Snake River Plain (WSRP) are limited. We undertook an integrated analysis of new and existing geophysical data and geologic constraints to study the crustal structure of the WSRP and generated two-dimensional crustal models across it. We observed both differences and similarities in the structural and tectonic evolution of the eastern and western arms of the SRP based on our integrated analysis. From a broader perspective based on recent geological and geophysical studies in the surrounding region, the intersection of the two arms of the SRP emerges as a major element of a complex tectonic intersection that includes the High Lava Plains of eastern Oregon, the Northern Nevada rift, a southwestern extension of the ESRP into northern Nevada, as well as, faulting and volcanism extending north-westward to connect with the Columbia River basalt plateau region. Thus, the goal of this study is to advance our understanding of the tectonomagmatic evolution of the region and to encourage further studies in the region.

Oceanic mafic magmatism in the Siletz terrane, NW North America: Fragments of an Eocene oceanic plateau?

The Siletz terrane, a predominantly mafic accreted oceanic terrane, is located in the Cascadia forearc region of Oregon, Washington and Vancouver Island. The terrane represents a late Palaeocene–Eocene large igneous province that consists of pillow lavas, massive flows and intrusive sheets. Previously it has been proposed that the Siletz terrane represents either an accreted oceanic plateau, hotspot island chain, backarc basin, island arc, or a sequence of slab window volcanics. A province-wide geochemical reassessment of the terrane, including new high precision Sr-Pb-Nd-Hf isotope data, has been used to assess the validity of the proposed tectonomagmatic models for the Siletz terrane. The trace element data show little evidence of crustal contamination, or an arc signature, and the samples have rare earth element (REE) patterns that are flat to light REE enriched. These features are similar to other oceanic plateaus such as the Ontong Java and the Caribbean. Initial isotope ratios range from 206Pb/204 Pb: 18.751 to 19.668, 207Pb/204Pb: 15.507 to 15.661, 208Pb/204Pb: 38.294 to 39.2128, 176Hf/177Hf: 0.28300 to 0.28316 (εHf: 9.0 to 14.5), 143Nd/144Nd: 0.51282 to 0.51299 (εNd: 5.0 to 8.1) and 87Sr/86Sr: 0.70302 to 0.70380. These data are consistent with a mantle source of the Siletz terrane that appears to have been heterogeneous and slightly enriched. The enriched signature has characteristics of both EM2 and HIMU components and this, combined with a calculated mantle potential temperature well above ambient mantle, indicates derivation of the Siletz magmatism from a mantle plume, possibly the Yellowstone Hotspot. We therefore conclude that the Siletz terrane represents an accreted oceanic plateau.

Geochronological and isotopic records of crustal storage and assimilation in the Wolverine Creek–Conant Creek system, Heise eruptive centre, Snake River Plain

Understanding the processes of differentiation of the Yellowstone–Snake River Plain (YSRP) rhyolites is typically impeded by the apparent lack of erupted intermediate compositions as well as the complex nature of their shallow interaction with the surrounding crust responsible for their typically low O isotopic ratios. A pair of normal-δ18O rhyolitic eruptions from the Heise eruptive centre in eastern Idaho, the Wolverine Creek Tuff and the Conant Creek Tuff, represent unique magmatic products of the Yellowstone hotspot preserving abundant vestiges of the intermediate differentiation steps leading to rhyolite generation. We address both shallow and deep processes of magma generation and storage in the two units by combining high-precision ID–TIMS U–Pb zircon geochronology, trace element, O and Hf isotopic studies of zircon, and Sr isotopic analyses of individual high-Mg# pyroxenes inherited from lower- to mid-crustal differentiation stages. The zircon geochronology confirms the derivation of both tuffs from the same rhyolitic magma reservoir erupted at 5.5941 ± 0.0097 Ma, preceded by at least 92 ± 14 ky of continuous or intermittent zircon saturation approximating the length of pre-eruptive magma accumulation in the upper crust. Some low-Mg# pyroxenes enclosing zircons predate the eruption by at least 45 ± 27 ky, illustrating the co-crystallisation of major and accessory phases in the near-liquidus rhyolitic melts of the YSRP over a significant period of time. Coeval zircon crystals are isotopically heterogeneous (two populations at eHf ~−5 and −13), requiring the assembly of isotopically distinct melt pockets directly prior to, or during, the eruption. The primitive Mg# 60–90 pyroxenes are out of isotopic equilibrium with the host rhyolitic melt (87Sr/86Sri = 0.70889), covering a range of 87Sr/86Sri = 0.70705–0.70883 corresponding to ratios typical of the most radiogenic YSRP basalts to the least radiogenic YSRP rhyolites. Together with the low eHf in zircon, the Sr isotopic ratios illustrate limited assimilation dominated by radiogenic Archean crustal source materials incorporated into variably evolved YSRP melts as they progress towards rhyolitic compositions by assimilation–fractional crystallisation.

Destruction of the Wyoming craton: Seismic evidence and geodynamic processes

Abstract Cratons are old and strong continental cores where the lithosphere is thick and remains largely undeformed for 2–3 b.y. Unlike typical cratons, the Wyoming craton underwent pervasive deformation ca. 80–55 Ma during the Laramide orogeny in the west-central United States, and has been subsequently encroached upon by the Yellowstone hotspot since 2.0 Ma. However, the mechanism for the deformation and the craton-hotspot interaction are not well understood. We present here a three-dimensional shear wave velocity model beneath the Wyoming craton constrained from Rayleigh wave data, which reveal new details about the cratonic lithosphere. The average lithosphere thickness beneath the craton is ∼150 km, significantly thinner than a normal cratonic root (>200 km). Continuous low velocities are observed beneath the Yellowstone hotspot and the Cheyenne belt. A low-velocity column is also present in the central-eastern craton at depths of 115–250 km. These low velocities can be explained by hot temperature and partial melting, implying mantle upwelling. A high-velocity anomaly with a dripping shape in central Wyoming extends to 200–250 km depth, indicating mantle downwelling and lithosphere erosion. Our model provides the first seismic evidence for complex small-scale mantle convection beneath the Wyoming craton. The convection probably developed during the subduction of the Farallon plate and has been reinforced by the Yellowstone hotspot. We propose that the combination of flat-slab subduction, small-scale convection, and hotspot activity can lead to massive destruction of a cratonic lithosphere.

Oxygen isotope evolution of the Lake Owyhee volcanic field, Oregon, and implications for the low-δ18O magmatism of the Snake River Plain–Yellowstone hotspot and other low-δ18O large igneous provinces

The Snake River Plain–Yellowstone (SRP-Y) hotspot track represents the largest known low-δ18O igneous province; however, debate persists regarding the timing and distribution of meteoric hydrothermal alteration and subsequent melting/assimilation relative to hotspot magmatism. To further constrain alteration relations for SRP-Y low-δ18O magmatism, we present in situ δ18O and U–Pb analyses of zircon, and laser fluorination δ18O analyses of phenocrysts, from the Lake Owyhee volcanic field (LOVF) of east-central Oregon. U–Pb data place LOVF magmatism between 16.3 and 15.4 Ma, and contain no evidence for xenocrystic zircon. LOVF δ18O(Zrc) values demonstrate (1) both low-δ18O and high-δ18O caldera-forming and pre-/post-caldera magmas, (2) relative increases in δ18O between low-δ18O caldera-forming and post-caldera units, and (3) low-δ18O magmatism associated with extension of the Oregon-Idaho Graben. The new data, along with new compilations of (1) in situ zircon δ18O data for the SRP-Y, and (2) regional δ18O(WR) and δ18O(magma) patterns, further constrain the thermal and structural associations for hydrothermal alteration in the SRP-Y. Models for low-δ18O magmatism must be compatible with (1) δ18O(magma) trends within individual SRP-Y eruptive centers, (2) along axis trends in δ18O(magma), and (3) the high concentration of low-δ18O magmas relative to the surrounding regions. When considered with the structural and thermal evolution of the SRP-Y, these constraints support low-δ18O magma genesis originating from syn-hotspot meteoric hydrothermal alteration, driven by hotspot-derived thermal fluxes superimposed on extensional tectonics. This model is not restricted to continental hotspot settings and may apply to several other low-δ18O igneous provinces with similar thermal and structural associations.