Crustal P-Wave Speed Structure Under Mount St. Helens From Local Earthquake Tomography

We used local earthquake data to model the P-wave speed structure of Mount St. Helens with the aim of improving our understanding of the active magmatic system. Our study used new data recorded by a dense array of 19 broadband seismographs that were deployed during the current eruption together with permanent network data recorded since the May 18, 1980 eruption. Most earthquakes around Mount St. Helens during the last 25 years were clustered in a narrow vertical column beneath the volcano from the surface to a depth of about 10 km. Earthquakes also occurred in a well-defined zone extending to the NNW from the volcano known as the St. Helens Seismic Zone (SHZ). During the current eruption, earthquakes have been confined to within 3 km of the surface beneath the crater floor. These earthquakes apparently radiate little shear-wave energy and the shear arrivals are usually contaminated by surface waves. Thus, we focused on developing an improved P- wave speed model. We used two data sources: (1) the short-period, vertical-component Pacific Northwest Seismograph Network and (2) new data recorded on a temporary array between June 2005 and February 2006. We first solved for a minimum one-dimensional model, incorporating the Moho depth found during an earlier wide-aperture refraction study. The three-dimensional model was solved simultaneously with hypocenter locations using the computer code SIMULPS14, extended for full three-dimensional ray shooting. We modified the code to force raypaths to remain below the ground surface. We began with large grid spacing and progressed to smaller grid spacing where the earthquakes and stations were denser. In this way we achieve a 40 km by 40 km regional model as well as a 10 km by 10 km fine-scale model directly beneath Mount St. Helens. The large-scale model is consistent with mapped geology and other geophysical data in the vicinity of Mount St. Helens. For example, there is a zone of relatively low velocities (-2% to -5% lower than background model) from 3 to at least 10 km depth extending NNW from the volcano parallel to the SHZ. The low-wave- speed zone coincides with a linear magnetic low, the western edge of a magnetotelluric conductive anomaly, and a localized gravity low. The coincidence of the volcano and these anomalies indicates this preexisting zone of weakness may control the location of Mount St. Helens, as has been suggested by previous investigators. Prominent high-wave-speed anomalies (+3% to +6% relative to background) on either side of this zone are due to plutons, which are also imaged with other geophysical data. Fine-scale modeling of the upper crust directly beneath Mount St. Helens reveals subtle structures not seen in the larger-scale model. The key structure is a cylindrical volume with speeds almost 10% slower than the background model extending from 6 to at least 10 km depth. The vertical, cylindrical volume of earthquakes, which reaches from the surface to more than 10 km depth, splits around this low-wave-speed volume creating an aseismic zone coincident with the low P-wave speeds. We interpret this volume as a melt-rich reservoir surrounded by hot rock.