skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Upper Crustal Structure and Magmatism in Southwest Washington: V p , V s , and V p / V s Results From the iMUSH Active‐Source Seismic Experiment

Structural details of the crust play an important role in controlling the distribution of volcanic activity in arc systems. In southwest Washington, several different regional structures associated with accretion and magmatism have been invoked to explain the broad distribution of Cascade volcanism in this region. In order to image these regional structures in the upper crust, Pg and Sg travel times from the imaging Magma Under St. Helens (iMUSH) active‐source seismic experiment are inverted forVp,Vs, andVp/Vsmodels in the region surrounding Mount St. Helens. Several features of these models provide new insights into the regional structure of the upper crust. A large section of the Southern Washington Cascades Conductor is imaged as a lowVp/Vsanomaly that is inferred to represent a broad sedimentary/metasedimentary sequence that composes the upper crust in this region. The accreted terrane Siletzia is imaged west of Mount St. Helens as north/south trending highVpandVp/Vsbodies. TheVp/Vsmodel shows relatively highVp/Vsregions near Mount St. Helens and the Indian Heaven Volcanic Field, which could be related to the presence of magmatic fluids. Separating these two volcanic regions below 6‐km depth is a northeast trending series of highVpandVsbodies. These bodies have the same orientation as several volcanic/magmatic features at the surface, including Mount St. Helens and Mount Rainier, and it is argued that these high‐velocity features are a regional‐scale group of intrusive bodies associated with a crustal weak zone that focuses magma ascent.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Page Range / eLocation ID:
p. 7067-7080
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Explosions and earthquakes are effectively discriminated byP/Samplitude ratios for moderate magnitude events (M ≥ 4) observed at regional to teleseismic distances (≥200 km). It is less clear ifP/Sratios are effective explosion discriminants for lower magnitudes observed at shorter distances. We report new tests ofP/Sdiscrimination using a dense seismic array in a continental volcanic arc setting near Mount St. Helens, with 23 single‐fired borehole explosions (ML0.9–2.3) and 406 earthquakes (ML1–3.3). The array provides up to 95 three‐component broadband seismographs, and most source‐receiver distances are <120 km. Additional insight is provided by ~3,000 vertical component geophone recordings of each explosion. Potential controls on local distanceP/Sratios are investigated, including frequency range, distance, magnitude, source depth, number of seismographs, and site effects. A frequency band of about 10–18 Hz performs better than lower or narrower bands because explosion‐inducedSwave amplitudes diminish relative toPfor higher frequencies. Source depth and magnitude exhibited weak influences onP/Sratios. Site responses for earthquakes and explosions are correlated with each other and with shallow crustalVpandVsfrom traveltime tomography. Overall, the results indicate high potential for local distanceP/Sexplosion discrimination in a continental volcanic arc setting, with ≥98% true positives and ≤6.3% false positives when using the array median from ≥16 stations. Performance is reduced for smaller arrays, especially those with ≤4 stations, thereby emphasizing the importance of array data for discrimination of low magnitude explosions.

    more » « less
  2. PhD Dissertation Abstract: The imposing andesite stratovolcano is the characteristic expression of subduction zone magmatism, posing hazards to coastal populations and bearing insight into deep Earth processes. On a map of a typical volcanic arc, one can easily distinguish the approximately linear alignment and regular spacing of these major edifices that stand out from a diffuse distribution of mafic volcanoes (e.g. the Quaternary Cascades; Hildreth, 2007). The andesitic composite volcanoes have a reputation for being complex, open systems: crystal zoning “stratigraphies,” diverse crystal cargoes including antecrysts or xenocrysts, quenched magmatic inclusions, and variations in isotopic signatures are among the many lines of evidence that these systems involve a variety of igneous processes and melt sources. To investigate the development and evolution of such transcrustal magma factories, I have conducted a detailed temporal, spatial, and geochemical characterization of a long-lived arc volcanic center in the southern Washington Cascades, the Goat Rocks volcanic complex. Results from ⁴⁰Ar/³⁹Ar and U/Pb geochronology constrain the lifespan of the Goat Rocks volcanic complex from ~3.1 Ma to ~100 ka. During this time, four major composite volcanoes were built (as well as several smaller volcanoes). From oldest to youngest, these are Tieton Peak, Bear Creek Mountain, Lake Creek volcano, and Old Snowy Mountain. Four volcanic stages are defined based on the lifespans of these centers and distinct compositional changes that occur from one to the next: Tieton Peak stage (3.1-2.6 Ma), Bear Creek Mountain stage (1.6-1.1 Ma), Lake Creek stage (1.1 Ma to 456 ka), and Old Snowy Mountain stage (440 ka to 115 ka). Two lava flow remnants also have ages in the interim between Tieton Peak stage and Bear Creek Mountain stage (2.3 Ma and 2.1 Ma), and their sources are not yet identified. The ages of the Bear Creek Mountain and Lake Creek stages in fact overlap, and the gap between Lake Creek stage and Old Snowy Mountain stage is only on the order of 10⁴ years. Based on supporting compositional evidence, the Bear Creek Mountain, Lake Creek, and Old Snowy Mountain stage volcanoes are considered to be the migrating surface expressions of a continuous magmatic system that was active over at least ~1.5 million years. It remains uncertain whether the gaps between the Tieton Peak stage, scattered early Pleistocene andesites, and Bear Creek Mountain stage are due to incomplete exposure/sampling or real quiescent periods earlier in the development of the Goat Rocks volcanic complex. Throughout the construction of the andesitic complex, mafic volcanoes were active on its periphery. These include the Miriam Creek volcano (~3.6-3.1 Ma), Devils Washbasin volcano (3.0-2.7 Ma), Hogback Mountain (1.1 Ma – 891 ka), Lakeview Mountain (194 ka), and Walupt Lake volcano (65 ka). Two basalt and basaltic andesite units (Qob₁ and Qob₂, 1.4 and 1.3 Ma; Hammond, 2017) also erupted from the Goat Rocks area, likely an older incarnation of Hogback Mountain. The suite of mafic magmas erupted in this region are all calcalkaline basalt (or basaltic andesite; CAB), but two compositional groups emerge from the trace element and isotopic data. Group 1 is LILE and LREE-enriched, with higher ⁸⁷Sr/⁸⁶Sr isotopes, and includes compositions from Devils Washbasin, Lower Hogback Mountain, and Lakeview Mountain. Group 2 is less enriched in LILE and LREE and lower in ⁸⁷Sr/⁸⁶Sr, and includes the compositions of Miriam Creek, Qob1, Upper Hogback Mountain, Walupt Lake, and Coleman Weedpatch. The two groups are recurrent through time and with no geographic distinction; in fact, both types were tapped by the Hogback Mountain volcano. Together both of these groups, alongside CABs from Mount Adams and various basalts from Mount St. Helens, form a compositional array between the basalts of the High Cascades and the intraplate-type basalts (IPB) of Mount Adams and Simcoe volcanic field. These results lead to three conclusions. 1) Variably subduction-modified mantle is distributed across the region, perhaps either as stratified layers or a web-like network of fluid pathways amongst less metasomatized mantle. 2) Transitional compositions between the IPBs and typical “High Cascades” CAB/HAOT signature suggest a broader influence of the mantle domain that feeds IPBs—if asthenospheric mantle through a slab window, as suggested by Mullen et al. (2017), then perhaps it bleeds in smaller quantities over a broader area. This compositional trend solidifies the interpretation of the southern Washington Cascades as a unique and coherent “segment” of the arc (the Washington segment of Pitcher and Kent, 2019). 3) The recurrence of variable mafic magma types through time, and with no geographic boundaries, indicates that the compositional evolution of the Goat Rocks volcanic complex was not likely driven by a change in mafic input. Indeed, the Sr, Nd, Hf, and Pb isotope ratios of the intermediate to felsic suite are closely aligned with the local basalts and suggest a limited role of crustal assimilation. Importantly, several mineral thermometers (zircon, ilmenite-magnetite pairs, and amphibole) align in recording higher crystallization temperatures in Bear Creek Mountain to early Lake Creek time, a cooling trend through the Lake Creek stage, and a more diverse range of temperatures in the transition to Old Snowy Mountain stage. Thus, it is suggested that the compositional evolution at Goat Rocks represents a thermal cycle of waxing and waning magmatic flux: where the period of Bear Creek Mountain to early Lake Creek volcanism was the climactic phase of a vertically extensive magma homogenization factory, then the system waned and cooled, ultimately losing its ability to filter, homogenize, and enrich magmas. 
    more » « less
  3. Abstract

    The lower crust has long been recognized as a place where significant chemical evolution of magma occurs. Understanding the conditions that lead to long‐term storage in or rapid transport through these lower crustal systems has remained a challenge. In this study, active‐source seismic data are processed to image prominent reflectors in the lower crust and upper mantle near Mount St. Helens. Results show strong and weak near‐Moho reflectivity located near previously inferred regions of magma storage and high‐density bodies underlying volcanic centers, respectively. To explain these observations, we hypothesize that magmas intersecting the high‐density lower crustal bodies experience enhanced vertical buoyancy forces/transport due to the increased density contrast between melt and host rock. The reduced reflectivity near high‐density bodies is therefore a result of both decreased impedance contrasts between crust and mantle host rocks and diminished lower crust magma accumulation, which ultimately causes focusing of volcanism above these features.

    more » « less
  4. Abstract

    Volcanic arcs consist of many distinct vents that are ultimately fueled by the common melting processes in the subduction zone mantle wedge. Seismic imaging of crustal‐scale magmatic systems can provide insight into how melt is organized in the deep crust and eventually focused beneath distinct vents as it ascends and evolves. Here, we investigate the crustal‐scale structure beneath a section of the Cascades arc spanning four major stratovolcanoes: Mt. Hood, Mt. St. Helens (MSH), Mt. Adams (MA), and Mt. Rainier, based on ambient noise data from 234 seismographs. Simultaneous inversion of Rayleigh and Love wave dispersion constrains the isotropic shear velocity (Vs) and identifies radially anisotropic structures. IsotropicVsshows two sub‐parallel low‐Vszones (∼3.45–3.55 km/s) at ∼15–30 km depth with one connecting Mt. Rainier to MA, and another connecting MSH to Mt. Hood, which are interpreted as deep crustal magma reservoirs containing up to ∼2.5%–6% melt, assuming near‐equilibrium melt geometry. Negative radial anisotropy, from vertical fractures like dikes, is prevalent in this part of the Cascadia, but is interrupted by positive radial anisotropy, from subhorizontal features like sills, extending vertically beneath MA and Mt. Rainier at ∼10–30 km depth and weaker and west‐dipping positive anisotropy beneath MSH. The positive anisotropy regions are adjacent to rather than co‐located with the isotropic low‐Vsanomalies. Ascending melt that stalled and mostly crystallized in sills with possible compositional differences from the country rock may explain the near‐averageVsand positive radial anisotropy adjacent to the active deep crustal magma reservoirs.

    more » « less
  5. Abstract

    At extensional volcanic arcs, faulting often acts to localize magmatism. Santorini is located on the extended continental crust of the Aegean microplate and is one of the most active volcanoes of the Hellenic arc, but the relationship between tectonism and magmatism remains poorly constrained. As part of the Plumbing Reservoirs Of The Earth Under Santorini experiment, seismic data were acquired across the Santorini caldera and the surrounding region using a dense amphibious array of >14,300 marine sound sources and 156 short‐period seismometers, covering an area 120 km by 45 km. Here aPwave velocity model of the shallow, upper‐crustal structure (<3‐km depth), obtained using travel time tomography, is used to delineate fault zones, sedimentary basins, and tectono‐magmatic lineaments. Our interpretation of tectonic boundaries and regional faults are consistent with prior geophysical studies, including the location of basin margins and E‐W oriented basement faults within the Christiana Basin west of Santorini. Reduced seismic velocities within the basement east of Santorini, near the Anydros and Anafi Basins, are coincident with a region of extensive NE‐SW faulting and active seismicity. The structural differences between the eastern and western sides of Santorini are in agreement with previously proposed models of regional tectonic evolution. Additionally, we find that regional magmatism has been localized in NE‐SW trending basin‐like structures that connect the Christiana, Santorini, and Kolumbo volcanic centers. At Santorini itself, we find that magmatism has been localized along NE‐SW trending lineaments that are subparallel to dikes, active faults, and regional volcanic chains. These results show strong interaction between magmatism and active deformation.

    more » « less