Geology topics

P. Dawson

Publications and source records attributed to P. Dawson.

11 recordsLinked to original sources

USGS2009

The SCEC/USGS dynamic earthquake rupture code verification exercise

Numerical simulations of earthquake rupture dynamics are now common, yet it has been difficult to test the validity of these simulations because there have been few field observations and no analytic solutions with which to compare the results. This paper describes the Southern California Earthquake Center/U.S. Geological Survey (SCEC/USGS) Dynamic Earthquake Rupture Code Verification Exercise, where codes that simulate spontaneous rupture dynamics in three dimensions are evaluated and the results produced by these codes are compared using Web-based tools. This is the first time that a broad and rigorous examination of numerous spontaneous rupture codes has been performed—a significant advance in this science. The automated process developed to attain this achievement provides for a future where testing of codes is easily accomplished. Scientists who use computer simulations to understand earthquakes utilize a range of techniques. Most of these assume that earthquakes are caused by slip at depth on faults in the Earth, but hereafter the strategies vary. Among the methods used in earthquake mechanics studies are kinematic approaches and dynamic approaches. The kinematic approach uses a computer code that prescribes the spatial and temporal evolution of slip on the causative fault (or faults). These types of simulations are very helpful, especially since they can be used in seismic data inversions to relate the ground motions recorded in the field to slip on the fault(s) at depth. However, these kinematic solutions generally provide no insight into the physics driving the fault slip or information about why the involved fault(s) slipped that much (or that little). In other words, these kinematic solutions may lack information about the physical dynamics of earthquake rupture that will be most helpful in forecasting future events. To help address this issue, some researchers use computer codes to numerically simulate earthquakes and construct dynamic, spontaneous rupture (hereafter called “spontaneous rupture”) solutions. For these types of numerical simulations, rather than prescribing the slip function at each location on the fault(s), just the friction constitutive properties and initial stress conditions are prescribed. The subsequent stresses and fault slip spontaneously evolve over time as part of the elasto-dynamic solution. Therefore, spontaneous rupture computer simulations of earthquakes allow us to include everything that we know, or think that we know, about earthquake dynamics and to test these ideas against earthquake observations.

Seismological Research Letters↗

USGS2006

Real-time monitoring and massive inversion of source parameters of very long period seismic signals: An application to Stromboli Volcano, Italy

We present a comprehensive processing tool for the real‐time analysis of the source mechanism of very long period (VLP) seismic data based on waveform inversions performed in the frequency domain for a point source. A search for the source providing the best‐fitting solution is conducted over a three‐dimensional grid of assumed source locations, in which the Green's functions associated with each point source are calculated by finite differences using the reciprocal relation between source and receiver. Tests performed on 62 nodes of a Linux cluster indicate that the waveform inversion and search for the best‐fitting signal over 100,000 point sources require roughly 30 s of processing time for a 2‐min‐long record. The procedure is applied to post‐processing of a data archive and to continuous automatic inversion of real‐time data at Stromboli, providing insights into different modes of degassing at this volcano.

Stromboli volcano↗

USGS2004

Application of near real-time radial semblance to locate the shallow magmatic conduit at Kilauea Volcano, Hawaii

Radial Semblance is applied to broadband seismic network data to provide source locations of Very-Long-Period (VLP) seismic energy in near real time. With an efficient algorithm and adequate network coverage, accurate source locations of VLP energy are derived to quickly locate the shallow magmatic conduit system at Kilauea Volcano, Hawaii. During a restart in magma flow following a brief pause in the current eruption, the shallow magmatic conduit is pressurized, resulting in elastic radiation from various parts of the conduit system. A steeply dipping distribution of VLP hypocenters outlines a region extending from sea level to about 550 m elevation below and just east of the Halemaumau Pit Crater. The distinct hypocenters suggest the shallow plumbing system beneath Halemaumau consists of a complex plexus of sills and dikes. An unconstrained location for a section of the conduit is also observed beneath the region between Kilauea Caldera and Kilauea Iki Crater.

Hawaii↗

USGS2003

Shallow-velocity models at the Kilauea Volcano, Hawaii, determined from array analyses of tremor wavefields

The properties of the surface wavefield at Kilauea Volcano are analysed using data from small-aperture arrays of short-period seismometers deployed in and around the Kilauea caldera. Tremor recordings were obtained during two Japan-US cooperative experiments conducted in 1996 and 1997. The seismometers were deployed in three semi-circular arrays with apertures of 300, 300 and 400 m, and a linear array with length of 1680 m. Data are analysed using a spatio-temporal correlation technique well suited for the study of the stationary stochastic wavefields of Rayleigh and Love waves associated with volcanic activity and scattering sources distributed in and around the summit caldera. Spatial autocorrelation coefficients are obtained as a function of frequency and are inverted for the dispersion characteristics of Rayleigh and Love waves using a grid search that seeks phase velocities for which the L-2 norm between data and forward modelling operators is minimized. Within the caldera, the phase velocities of Rayleigh waves range from 1400 to 1800 m s -1 at 1 Hz down to 300-400 m s -1 at 10 Hz, and the phase velocities of Love waves range from 2600 to 400 m s -1 within the same frequency band. Outside the caldera, Rayleigh wave velocities range from 1800 to 1600 m s -1 at 1 Hz down to 260-360 m s -1 at 10 Hz, and Love wave velocities range from 600 to 150 m s -1 within the same frequency band. The dispersion curves are inverted for velocity structure beneath each array, assuming these dispersions represent the fundamental modes of Rayleigh and Love waves. The velocity structures observed at different array sites are consistent with results from a recent 3-D traveltime tomography of the caldera region, and point to a marked velocity discontinuity associated with the southern caldera boundary.

Geophysical Journal International↗

USGS2002

Source mechanism of very-long-period signals accompanying dome growth activity at Merapi volcano, Indonesia

Very-long-period (VLP) pulses with period of 6–7s, displaying similar waveforms, were identified in 1998 from broadband seismographs around the summit crater. These pulses accompanied most of multiphase (MP) earthquakes, a type of long-period event locally defined at Merapi Volcano. Source mechanisms for several VLP pulses were examined by applying moment tensor inversion to the waveform data. Solutions were consistent with a crack striking ∼70° and dipping ∼50° SW, 100m under the active dome, suggest pressurized gas transport involving accumulation and sudden release of 10–60 m 3 of gas in the crack over a 6s interval.

Merapi volcano↗

USGS2002

Mapping the sources of the seismic wave field at Kilauea volcano, Hawaii, using data recorded on multiple seismic Antennas

Seismic antennas constitute a powerful tool for the analysis of complex wave fields. Well-designed antennas can identify and separate components of a complex wave field based on their distinct propagation properties. The combination of several antennas provides the basis for a more complete understanding of volcanic wave fields, including an estimate of the location of each individual wave-field component identified simultaneously by at least two antennas. We used frequency-slowness analyses of data from three antennas to identify and locate the different components contributing to the wave fields recorded at Kilauea volcano, Hawaii, in February 1997. The wave-field components identified are (1) a sustained background volcanic tremor in the form of body waves generated in a shallow hydrothermal system located below the northeastern edge of the Halemaumau pit crater; (2) surface waves generated along the path between this hydrothermal source and the antennas; (3) back-scattered surface wave energy from a shallow reflector located near the southeastern rim of Kilauea caldera; (4) evidence for diffracted wave components originating at the southeastern edge of Halemaumau; and (5) body waves reflecting the activation of a deeper tremor source between 02 hr 00 min and 16 hr 00 min Hawaii Standard Time on 11 February.

Hawaii↗

USGS2001

Spatial extent of a hydrothermal system at Kilauea Volcano, Hawaii, determined from array analyses of shallow long-period seismicity 1. Method

We present a probabilistic method to locate the source of seismic events using seismic antennas. The method is based on a comparison of the event azimuths and slownesses derived from frequency-slowness analyses of array data, with a slowness vector model. Several slowness vector models are considered including both homogeneous and horizontally layered half-spaces and also a more complex medium representing the actual topography and three-dimensional velocity structure of the region under study. In this latter model the slowness vector is obtained from frequency-slowness analyses of synthetic signals. These signals are generated using the finite difference method and include the effects of topography and velocity structure to reproduce as closely as possible the behavior of the observed wave fields. A comparison of these results with those obtained with a homogeneous half-space demonstrates the importance of structural and topographic effects, which, if ignored, lead to a bias in the source location. We use synthetic seismograms to test the accuracy and stability of the method and to investigate the effect of our choice of probability distributions. We conclude that this location method can provide the source position of shallow events within a complex volcanic structure such as Kilauea Volcano with an error of ±200 m.

Hawai'i↗

USGS2001

Wavefield properties of a shallow long-period event and tremor at Kilauea Volcano, Hawaii

The wavefields of tremor and a long-period (LP) event associated with the ongoing eruptive activity at Kilauea Volcano, Hawaii, are investigated using a combination of dense small-aperture (300 m) and sparse large-aperture (5 km) arrays deployed in the vicinity of the summit caldera. Measurements of azimuth and slowness for tremor recorded on the small-aperture array indicate a bimodal nature of the observed wavefield. At frequencies below 2 Hz, the wavefield is dominated by body waves impinging the array with steep incidence. These arrivals are attributed to the oceanic microseismic noise. In the 2-6 Hz band, the wavefield is dominated by waves propagating from sources located at shallow depths (<1 km) beneath the eastern edge of the Halemaumau pit crater. The hypocenter of the LP event, determined from frequency-slowness analyses combined with phase picks, appears to be located close to the source of tremor but at a shallower depth (<0.1 km). The wavefields of tremor and LP event are characterized by a complex composition of body and surface waves, whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region. Analyses of the directional properties of the wavefield in the 2-6 Hz band point to the directions of main scattering sources, which are consistent with pronounced velocity contrasts imaged in a high-resolution three-dimensional velocity model of the caldera region. The frequency and Q of the dominant peak observed in the spectra of the LP event may be explained as the dominant oscillation mode of a crack with scale length 20-100 m and aperture of a few centimeters filled with bubbly water. The mechanism driving the shallow tremor appears to be consistent with a sustained excitation originating in the oscillations of a bubbly cloud resulting from vesiculation and degassing in the magma. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

USGS1999

Broadband measurements of the sources of explosions at Stromboli Volcano, Italy

During September–October 1997, 21 three-component broadband seismometers were deployed on Stromboli Volcano at radial distances of 0.3–2.2 km from the active crater to investigate the source mechanisms of Strombolian explosions. In the 2–50 s band, the very-long period (VLP) signals associated with explosions are consistent with two stationary sources repeatedly activated in time. VLP particle motions are essentially linear and analyses of semblance and particle motions are consistent with a source centroid offset 300 m beneath and 300 m northwest of the active vents. Similar VLP waveforms are observed at all 21 stations, indicating that the seismograms are source-dominated. The VLP ground displacement response to each explosion may be qualitatively interpreted as: (1) pressurization of the conduit associated with the ascent of a slug of gas; (2) depressurization of the conduit in response to mass withdrawal during the eruption; and (3) repressurization of the conduit associated with the replenishment of the source with fluid.

Geophysical Research Letters↗

USGS1998

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii

We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 μrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5–10 s followed by an exponential recovery lasting 1–3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located 1 km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half-space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 10 14 N m and 3000 m 3 , respectively. Our results also point to the existence of a single force component with amplitude of 10 9 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5-hour-long activation of the pulsating source is roughly 500,000 m 3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Shallow velocity structure of Stromboli Volcano, Italy, derived from small-aperture array measurements of Strombolian tremor

The properties of the tremor wave field at Stromboli are analyzed using data from small-aperture arrays of short-period seismometers deployed on the north flank of the volcano. The seismometers are configured in two semi-circular arrays with radii of 60 and 150 m and a linear array with length of 600 m. The data are analyzed using a spatiotemporal correlation technique specifically designed for the study of the stationary stochastic wave field of Rayleigh and Love waves generated by volcanic activity and by scattering sources distributed within the island. The correlation coefficients derived as a function of frequency for the three components of motion clearly define the dispersion characteristics for both Rayleigh and Love waves. Love and Rayleigh waves contribute 70% and 30%, respectively, of the surface-wave power. The phase velocities of Rayleigh waves range from 1000 m/sec at 2 Hz to 350 m/sec at 9 Hz, and those for Love waves range from 800 to 400 m/sec over the same frequency band. These velocities are similar to those measured near Puu Oo on the east rift of Kilauea Volcano, Hawaii, although the dispersion characteristics of Rayleigh waves at Stromboli show a stronger dependence on frequency. Such low velocities are consistent with values expected for densely cracked solidified basalt. The dispersion curves are inverted for a velocity model beneath the arrays, assuming those dispersions represent the fundamental modes of Rayleigh and Love waves.

Stromboli volcano↗