USGS Science⌕ Search

Geology topics

P.B. Dawson

Publications and source records attributed to P.B. Dawson.

15 recordsLinked to original sources

Imaging the Juan de Fuca plate beneath southern Oregon using teleseismic P wave residuals

The Juan de Fuca plate subducts under the western margin of the North American plate in southern British Columbia, Washington, Oregon, and northern California. Benioff zone seismicity delineates shallow parts of the plate in Washington and northern California, whereas in Oregon there is a lack of seismicity. This study images the Juan de Fuca plate in southern Oregon using seismic tomography. We inverted P wave travel time residuals from a 366-km-long seismic array operated in southern Oregon in 1982. The southeast striking array extended from the Coast ranges to the Modoc Plateau and crossed the High Cascades at Crater Lake, Oregon. We imaged three features under the array: one high-velocity zone and two low-velocity zones. The high-velocity zone is 3–4% faster than the surrounding upper mantle. It dips steeply at 65° to the east beneath the Cascade Range and extends down to at least 200 km. We propose that this high-velocity feature is subducted Juan de Fuca plate. We also imaged two low-velocity zones, both of which are 3–4% slower than the surrounding earth structure. The southeastern low-velocity zone may be caused by partially molten crust underlying the Crater Lake volcano region.

Journal of Geophysical Research↗

Teleseismic tomography of the compressional wave velocity structure beneath the Long Valley region, California

In 1982 and 1984 the U.S. Geological Survey used several seismic networks, totaling over 90 stations, to record teleseismic P waves and measure travel time residuals in an area centered on the Long Valley caldera. We inverted the travel time residuals to obtain a three-dimensional image of the velocity structure with resolution of 5–6 km to depths of 70 km beneath the array. Direct inversion of these data indicates that the 2- to 4-km-thick low-velocity caldera fill contaminates the signal from any midcrustal velocity anomalies beneath the caldera. Thus two methods were used to strip the effects of the upper crust from the travel time residuals: (1) ray tracing through upper crustal velocity models provided by seismic refraction experiments and gravity surveys, and (2) an iterative stripping scheme using the inversion itself. The methods produce essentially identical results and adequately remove the effects of the shallowest crustal structures, including the caldera fill and hydrothermal alteration effects. The resulting “stripped” models show two well-resolved midcrustal low-velocity bodies in the Long Valley region. The first body is centered between 7 and 20 km depth beneath the resurgent dome of the Long Valley caldera and has a volume of 150–600 km 3 . The second, with a similar volume, is centered between 10 and 20 km depth beneath the Mono Craters, about 10 km north of Long Valley. Velocity contrasts in both of these bodies are about 6–10%, and the features are interpreted as silicic magma chambers. This experiment does not preclude the presence of additional pockets of magma smaller than 5 km across in the upper crust, particularly beneath the resurgent dome of the caldera (which would be removed with the stripping methods). It is likely that the midcrust anomaly beneath the resurgent dome is a remnant of the caldera forming magma chamber. Upper mantle velocities are lower than average beneath Mono Craters and higher than average beneath Long Valley. We tentatively interpret these anomalies as mantle partial melt beneath the Mono Craters and refractory upper mantle residuum beneath Long Valley. The high eruptive rate of the Mono Craters and these upper mantle structures suggest that the focus of volcanism is shifting north from Long Valley to the Mono Craters.

Journal of Geophysical Research Solid Earth↗

Determining the seismic source mechanism and location for an explosive eruption with limited observational data: Augustine Volcano, Alaska

Waveform inversions of the very-long-period components of the seismic wavefield produced by an explosive eruption that occurred on 11 January, 2006 at Augustine Volcano, Alaska constrain the seismic source location to near sea level beneath the summit of the volcano. The calculated moment tensors indicate the presence of a volumetric source mechanism. Systematic reconstruction of the source mechanism shows the source consists of a sill intersected by either a sub-vertical east-west trending dike or a sub-vertical pipe and a weak single force. The trend of the dike may be controlled by the east-west trending Augustine-Seldovia arch. The data from the network of broadband sensors is limited to fourteen seismic traces, and synthetic modeling confirms the ability of the network to recover the source mechanism. The synthetic modeling also provides a guide to the expected capability of a broadband network to resolve very-long-period source mechanisms, particularly when confronted with limited observational data. Copyright 2011 by the American Geophysical Union.

Geophysical Research Letters↗

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma-filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Journal of Geophysical Research B: Solid Earth↗

Source process of a long-period event at Kilauea volcano, Hawaii

We analyse a long-period (LP) event observed by a dense seismic network temporarily operated at Kilauea volcano, Hawaii, in 1996. We systematically perform spectral analyses, waveform inversions and forward modeling of the LP event to quantify its source process. Spectral analyses identify two dominant spectral frequencies at 0.6 and 1.3 Hz with associated Q values in the range 10-20. Results from waveform inversions assuming six moment-tensor and three single-force components point to the resonance of a horizontal crack located at a depth of approximately 150 m near the northeastern rim of the Halemaumau pit crater. Waveform simulations based on a fluid-filled crack model suggest that the observed frequencies and Q values can be explained by a crack filled with a hydrothermal fluid in the form of either bubbly water or steam. The shallow hydrothermal crack located directly above the magma conduit may have been heated by volcanic gases leaking from the conduit. The enhanced flux of heat raised the overall pressure of the hydrothermal fluid in the crack and induced a rapid discharge of fluid from the crack, which triggered the acoustic vibrations of the resonator generating the LP waveform. The present study provides further support to the idea that LP events originate in the resonance of a crack. 2005 RAS.

Hawaii↗

Tomographic inversion of P-wave velocity and Q structures beneath the Kirishima volcanic complex, Southern Japan, based on finite difference calculations of complex traveltimes

We estimate the P-wave velocity and attenuation structures beneath the Kirishima volcanic complex, southern Japan, by inverting the complex traveltimes (arrival times and pulse widths) of waveform data obtained during an active seismic experiment conducted in 1994. In this experiment, six 200-250 kg shots were recorded at 163 temporary seismic stations deployed on the volcanic complex. We use first-arrival times for the shots, which were hand-measured interactively. The waveform data are Fourier transformed into the frequency domain and analysed using a new method based on autoregressive modelling of complex decaying oscillations in the frequency domain to determine pulse widths for the first-arrival phases. A non-linear inversion method is used to invert 893 first-arrival times and 325 pulse widths to estimate the velocity and attenuation structures of the volcanic complex. Wavefronts for the inversion are calculated with a finite difference method based on the Eikonal equation, which is well suited to estimating the complex traveltimes for the structures of the Kirishima volcano complex, where large structural heterogeneities are expected. The attenuation structure is derived using ray paths derived from the velocity structure. We obtain 3-D velocity and attenuation structures down to 1.5 and 0.5 km below sea level, respectively. High-velocity pipe-like structures with correspondingly low attenuation are found under the summit craters. These pipe-like structures are interpreted as remnant conduits of solidified magma. No evidence of a shallow magma chamber is visible in the tomographic images.

Geophysical Journal International↗

Three-dimensional velocity structure of the Kilauea caldera, Hawaii

High-resolution velocity models (0.5 km resolution) of the Kilauea caldera region are obtained by the tomographic inversion of both P - and S -wave arrival times. Data are from the permanent Hawaiian Volcano Observatory (HVO) seismic network, a broadband seismic network, and a temporary array of stations centered on the southern boundary of the caldera. A low-velocity P -wave anomaly is imaged centered on the southeastern edge of the caldera, with a velocity contrast of about 10% and a volume of 27 km³. The V P /V S model mimics the spatial extent of the P-wave anomaly, but is partitioned into two discrete anomalous volumes centered on the southern boundary of the caldera and on the upper east rift of the volcano. The corresponding Poisson's ratio in these zones is high (ν=0.25–0.32) which is consistent with a densely-cracked, hot volume which may contain partial melt. The large-scale features of the models are consistent with results obtained from an earlier, larger-scale (2 km resolution) tomographic image of Kilauea Volcano based on HVO network data.

Geophysical Research Letters↗

Three-dimensional P and S wave velocity structure of Redoubt Volcano, Alaska

The three‐dimensional P and S wave structure of Redoubt Volcano, Alaska, and the underlying crust to depths of 7–8 km is determined from 6219 P wave and 4008 S wave first‐arrival times recorded by a 30‐station seismograph network deployed on and around the volcano. First‐arrival times are calculated using a finite‐difference technique, which allows for flexible parameterization of the slowness model and easy inclusion of topography and source‐receiver geometry. The three‐dimensional P wave velocity structure and hypocenters are determined simultaneously, while the three‐dimensional S wave velocity model is determined using the relocated seismicity and an initial S wave velocity model derived from the P wave velocity model assuming an average Vp / Vs ratio of 1.78. Convergence is steady with approximately 73% and 52% reduction in P and S wave arrival time RMS, respectively, after 10 iterations. The most prominent feature observed in the three‐dimensional velocity models derived for both P and S waves is a relative low‐velocity, near‐vertical, pipelike structure approximately 1 km in diameter that extends from 1 to 6 km beneath sea level. This feature aligns axially with the bulk of seismicity and is interpreted as a highly fractured and altered zone encompassing a magma conduit. The velocity structure beneath the north flank of the volcano between depths of 1 and 6 km is characterized by large lateral velocity variations. High velocities within this region are interpreted as remnant dikes and sills and low velocities as regions along which magma migrates. No large low‐velocity body suggestive of a magma chamber is resolved in the upper 7–8 km of the crust.

Journal of Geophysical Research B: Solid Earth↗