Tobias Diehl

Tomography from 26 years of seismicity revealing that the spatial extent of the Yellowstone crustal magma reservoir extends well beyond the Yellowstone caldera

The Yellowstone volcanic field has experienced three of Earths most explosive volcanic eruptions in the last 2.1 Ma. The most recent eruption occurred 0.64 Ma forming the 60 km long Yellowstone caldera. We have compiled earthquake data from the Yellowstone Seismic Network from 1984 to 2011 and tomographically imaged the three-dimensional P wave velocity (Vp) structure of the Yellowstone volcanic system. The resulting model reveals a large, low Vp body, interpreted to be the crustal magma reservoir that has fueled Yellowstones youthful volcanism. Our imaged magma body is 90 km long, 5–17 km deep, and 2.5 times larger than previously imaged. The magma body extends ~15 km NE of the caldera and correlates with the location of the largest negative gravity anomaly, a −80 mGal gravity low. This new seismic image provides important constraints on the dynamics of the Yellowstone magma system and its potential for future volcanic eruptions and earthquakes.

Automatic S-Wave Picker for Local Earthquake Tomography

High-resolution seismic tomography at local and regional scales requires large and consistent sets of arrival-time data. Algorithms combining accurate pick- ing with an automated quality classification can be used for repicking waveforms and compiling large arrival-time data sets suitable for tomographic inversion. S-wave velocities represent a key parameter for petrological interpretation, improved hypocenter determination, as well as for seismic hazard models. In our approach, we combine three commonly used phase detection and picking methods in a robust S-wave picking procedure. Information from the different techniques provides an in situ estimate of timing uncertainty and of the reliability of the automatic phase identification. Automatic picks are compared against manually picked reference picks of selected earthquakes in the Alpine region. The average accuracy of automatic picks and their classification is comparable with the reference picks, although a higher number of picks is downgraded to lower quality classes by the automatic picker. In the production-mode, we apply the picker to a data set of 552 earthquakes in the Alps recorded at epicentral distances ≤150 km. The resulting data set includes about 2500 S phases with an upper error bound of 0.27 sec. Online Material: Details on the proposed automatic S-wave picking algorithm.

Splay faults imaged by fluid-driven aftershocks of the 2004 Mw 9.2 Sumatra-Andaman earthquake

High-precision teleseismic double-difference locations and focal mechanisms of aftershocks of the A.D. 2004 M w 9.2 Sumatra-Andaman earthquake illuminate an active imbricate fault system in the accretionary prism offshore northern Sumatra. They reveal repeated failure of a shallow northeast-dipping thrust fault above the megathrust, which we interpret to be the reactivation of a splay fault that rises from the megathrust at ∼55 km depth and cuts through the overriding Sunda plate. The projected intersection of the splay fault with the seafloor correlates with a recently active thrust fault seen in postseismic bathymetry data west of the Aceh basin. A spatiotemporal analysis of the streaky aftershock distribution indicates that ascending fluids released from the subducting oceanic crust along inherited seafloor fabric may control brittle fracture in the overriding plate. We speculate that if the splay fault was active coseismically, it may have led to amplified vertical uplift of the forearc ridge and contributed to generating the cataclysmic near-field tsunami that struck the northwest Sumatra coast following the 2004 rupture.

The effects of data quality in local earthquake tomography: Application to the Alpine region

Despite the increase in quality and number of seismic stations in many parts of the world, accurate timing of individual arrival times remains crucial for many tomographic applications. To achieve a data set of high quality, arrival times need to be picked with high accuracy, including a proper assessment of the uncertainty of timing and phase identification, and a high level of consistency. We have investigated the effects of data quantity and quality on the solution quality in local earthquake tomography. We have compared tomographic results obtained with synthetic and real data of two very different data sets. The first data set consisted of a large set of arrival times of low precision and unknown accuracy taken from the International Seismological Centre (ISC) Bulletin for the greater Alpine region. The second high-quality data set for the same region was seven times smaller and was obtained by automated quality-weighted repicking. During a first series of inversions, synthetic data resembling the two data sets were inverted with the same amount of Gaussian distributed noise added. Subsequently, during a second series of inversions, the noise level was increased successively for ISC data to study the effect of larger Gaussian distributed error on the solution quality. Finally, the real data for both data sets were inverted. These investigations showed that, for Gaussian distributed error, a smaller data set of high quality could achieve a similar or better solution quality than a data set seven times larger but about four times lower in quality. Our results further suggest that the quality of the ISC Bulletin is degraded significantly by inconsistencies, strongly limiting the use of this large data set for local earthquake tomography studies.