Victor C. Tsai

Earthquake in a Maze: Compressional Rupture Branching During the 2012 Mw 8.6 Sumatra Earthquake

Earthquake in a Maze The 11 April 2012 magnitude 8.6 earthquake offshore of Sumatra was the largest measured earthquake along a strike-slip boundary that modern seismological instruments have ever recorded. Despite its size and proximity to a large population, there was no subsequent tsunami and there were no reported fatalities. Meng et al. (p. 724, published online 19 July) used teleseismic data from seismological networks in Japan and Europe to image the source of high-frequency radiation generated by the earthquake to understand the mechanics of this unique event. The resultant back projections showed that the earthquake slowly ruptured along a complex series of faults. The deeper-than-usual rupture path and large stress drop are both features that may not be unique to this earthquake, suggesting that regions in a similar tectonic environment may have the potential for more complex—or larger—intraplate earthquakes than might have been expected. The mechanics of the largest strike-slip earthquake ever recorded give clues about how intraplate earthquakes rupture. Seismological observations of the 2012 moment magnitude 8.6 Sumatra earthquake reveal unprecedented complexity of dynamic rupture. The surprisingly large magnitude results from the combination of deep extent, high stress drop, and rupture of multiple faults. Back-projection source imaging indicates that the rupture occurred on distinct planes in an orthogonal conjugate fault system, with relatively slow rupture speed. The east-southeast–west-northwest ruptures add a new dimension to the seismotectonics of the Wharton Basin, which was previously thought to be controlled by north-south strike-slip faulting. The rupture turned twice into the compressive quadrant, against the preferred branching direction predicted by dynamic Coulomb stress calculations. Orthogonal faulting and compressional branching indicate that rupture was controlled by a pressure-insensitive strength of the deep oceanic lithosphere.

Multiple CMT source analysis of the 2004 Sumatra earthquake

While it is agreed that the great Sumatra earthquake of December 26, 2004 was among the largest earthquakes of the past century, there has been disagreement on how large it was, which part of the fault ruptured, and how the rupture took place. We present a centroid-moment-tensor (CMT) analysis of the earthquake in which multiple point sources are used in the inversion to mimic a propagating slip pulse. The final model consists of five point sources, with the southernmost sources accounting for the majority of the moment release. The presumed fault planes of the southern sources strike northwest, while those in the north strike northeast, consistent with the geometry of the subduction trench. Slip on the fault is found to be more oblique in the north than in the south. The inversion with five sources leads to a moment magnitude for the Sumatra earthquake of M_W = 9.3, consistent with estimates from long-period normal-mode amplitudes.

Seasonality and Increasing Frequency of Greenland Glacial Earthquakes

Some glaciers and ice streams periodically lurch forward with sufficient force to generate emissions of elastic waves that are recorded on seismometers worldwide. Such glacial earthquakes on Greenland show a strong seasonality as well as a doubling of their rate of occurrence over the past 5 years. These temporal patterns suggest a link to the hydrological cycle and are indicative of a dynamic glacial response to changing climate conditions.

Ice-front variation and tidewater behavior on Helheim and Kangerdlugssuaq Glaciers, Greenland

We used satellite images to examine the calving behavior of Helheim and Kangerdlugssuaq Glaciers, Greenland, from 2001 to 2006, a period in which they retreated and sped up. These data show that many large iceberg-calving episodes coincided with teleseismically detected glacial earthquakes, suggesting that calving-related processes are the source of the seismicity. For each of several events for which we have observations, the ice front calved back to a large, pre-existing rift. These rifts form where the ice has thinned to near flotation as the ice front retreats down the back side of a bathymetric high, which agrees well with earlier theoretical predictions. In addition to the recent retreat in a period of higher temperatures, analysis of several images shows that Helheim retreated in the 20th Century during a warmer period and then re-advanced during a subsequent cooler period. This apparent sensitivity to warming suggests that higher temperatures may promote an initial retreat off a bathymetric high that is then sustained by tidewater dynamics as the ice front retreats into deeper water. The cycle of frontal advance and retreat in less than a century indicates that tidewater glaciers in Greenland can advance rapidly. Greenlands larger reservoir of inland ice and conditions that favor the formation of ice shelves likely contribute to the rapid rates of advance.

The Yellowstone magmatic system from the mantle plume to the upper crust

Yellowstones missing magmatic link Yellowstone is an extensively studied “supervolcano” that has a large supply of heat coming from a pool of magma near the surface and the mantle below. A link between these two features has long been suspected. Huang et al. imaged the lower crust using seismic tomography (see the Perspective by Shapiro and Koulakov). Their findings provide an estimate of the total amount of molten rock beneath Yellowstone and help to explain the large amount of volcanic gases escaping from the region. Science, this issue p. 773; see also p. 758 The Yellowstone supervolcano has a large magma body between the mantle hot spot and the upper crustal magmatic reservoir. [Also see Perspective by Shapiro and Koulakov] The Yellowstone supervolcano is one of the largest active continental silicic volcanic fields in the world. An understanding of its properties is key to enhancing our knowledge of volcanic mechanisms and corresponding risk. Using a joint local and teleseismic earthquake P-wave seismic inversion, we revealed a basaltic lower-crustal magma body that provides a magmatic link between the Yellowstone mantle plume and the previously imaged upper-crustal magma reservoir. This lower-crustal magma body has a volume of 46,000 cubic kilometers, ~4.5 times that of the upper-crustal magma reservoir, and contains a melt fraction of ~2%. These estimates are critical to understanding the evolution of bimodal basaltic-rhyolitic volcanism, explaining the magnitude of CO2 discharge, and constraining dynamic models of the magmatic system for volcanic hazard assessment.

Understanding the amplitudes of noise correlation measurements

Cross correlation of ambient seismic noise is known to result in time series from which station-station travel-time measurements can be made. Part of the reason that these cross-correlation travel-time measurements are reliable is that there exists a theoretical framework that quantifies how these travel times depend on the features of the ambient noise. However, corresponding theoretical results do not currently exist to describe how the amplitudes of the cross correlation depend on such features. For example, currently it is not possible to take a given distribution of noise sources and calculate the cross correlation amplitudes one would expect from such a distribution. Here, we provide a ray-theoretical framework for calculating cross correlations. This framework differs from previous work in that it explicitly accounts for attenuation as well as the spatial distribution of sources and therefore can address the issue of quantifying amplitudes in noise correlation measurements. After introducing the general framework, we apply it to two specific problems. First, we show that we can quantify the amplitudes of coherency measurements, and find that the decay of coherency with station-station spacing depends crucially on the distribution of noise sources. We suggest that researchers interested in performing attenuation measurements from noise coherency should first determine how the dominant sources of noise are distributed. Second, we show that we can quantify the signal-to-noise ratio of noise correlations more precisely than previous work, and that these signal-to-noise ratios can be estimated for given situations prior to the deployment of seismometers. It is expected that there are applications of the theoretical framework beyond the two specific cases considered, but these applications await future work.

A model for turbulent hydraulic fracture and application to crack propagation at glacier beds

Glaciological observations of under-flooding suggest that fluid-induced hydraulic fracture of an ice sheet from its bed sometimes occurs quickly, possibly driven by turbulently flowing water in a broad sheet flow. Taking the approximation of a fully turbulent flow into an elastic ice medium with small fracture toughness, we derive an approximate expression for the crack-tip speed, opening displacement and pressure profile. We accomplish this by first showing that a Manning-Strickler channel model for resistance to turbulent flow leads to a mathematical structure somewhat similar to that for resistance to laminar flow of a power law viscous fluid. We then adapt the plane-strain asymptotic crack solution of Desroches et al. (1994) and the power law self-similar solution of Adachi and Detournay (2002) for that case to calculate the desired quantities. The speed of crack growth is shown to scale as the overpressure (in excess of ice overburden) to the power 7/6, inversely as ice elastic modulus to the power 2/3, and as the ratio of crack length to wall roughness scale to the power 1/6. We tentatively apply our model by choosing parameter values thought appropriate for a basal crack driven by the rapid drainage of a surface meltwater lake near the margin of the Greenland Ice Sheet. Making various approximations perhaps relevant to this setting, we estimate fluid inflow rate to the basal fracture and vertical and horizontal surface displacements and find order-of-magnitude agreement with observations by Das et al. (2008) associated with lake drainage. Finally, we discuss how these preliminary estimates could be improved.

Joint inversion of Rayleigh wave phase velocity and ellipticity using USArray: Constraining velocity and density structure in the upper crust

Rayleigh wave ellipticity, or H/V ratio, observed on the surface is particularly sensitive to shallow earth structure. In this study, we jointly invert measurements of Rayleigh wave H/V ratio and phase velocity between 24–100 and 8–100 sec period, respectively, for crust and upper mantle structure beneath more than 1000 USArray stations covering the western United States. Upper crustal structure, in particular, is better constrained by the joint inversion compared to inversions based on phase velocities alone. In addition to imaging Vs structure, we show that the joint inversion can be used to constrain Vp/Vs and density in the upper crust. New images of uppermost crustal structure (<3 km depth) are in excellent agreement with known surface features, with pronounced low Vs, low density, and high Vp/Vs anomalies imaged in the locations of several major sedimentary basins including the Williston, Powder River, Green River, Denver, and San Juan basins. These results demonstrate not only the consistency of broadband H/V ratios and phase velocity measurements, but also that their complementary sensitivities have the potential to resolve density and Vp/Vs variations.

A physical model for seismic noise generation from sediment transport in rivers

Measuring sediment flux in rivers remains a significant problem in studies of landscape evolution. Recent studies suggest that observations of seismic noise near rivers can help provide such measurements, but the lack of models linking observed seismic quantities to sediment flux has prevented the method from being used. Here, we develop a forward model to describe the seismic noise induced by the transport of sediment in rivers. The model provides an expression for the power spectral density (PSD) of the Rayleigh waves generated by impulsive impacts from saltating particles which scales linearly with the number of particles of a given size and the square of the linear momentum. After incorporating expressions for the impact velocity and rate of impacts for fluvially transported sediment, we observe that the seismic noise PSD is strongly dependent on the sediment size, such that good constraints on grain size distribution are needed for reliable estimates of sediment flux based on seismic noise observations. The model predictions for the PSD are consistent with recent measurements and, based on these data, a first attempt at inverting seismic noise for the sediment flux is provided.

Analysis of glacial earthquakes

In 2003, Ekstrom et al. reported on the detection of a new class of earthquakes that occur in glaciated regions, with the vast majority being in Greenland. The events have a characteristic radiation pattern and lack the high-frequency content typical of tectonic earthquakes. It was proposed that the events correspond to large and sudden sliding motion of glaciers. Here we present an analysis of all 184 such events detected in Greenland between 1993 and 2005. Fitting the teleseismic long-period surface waves to a landslide model of the source, we obtain improved locations, timing, force amplitudes, and force directions. After relocation, the events cluster into seven regions, all of which correspond to regions of very high ice flow and most of which are named outlet glaciers. These regions are Daugaard Jensen Glacier, Kangerdlugssuaq Glacier, Helheim Glacier, the southeast Greenland glaciers, the northwest Greenland glaciers, Rinks Isbrae, and Jakobshavn Isbrae. Event amplitudes range from 0.1 to 2.0 × 10^(14) kg m. Force directions are consistent with sliding in the direction of glacial flow over a period of about 50 s. Each region has a different temporal distribution of events. All glaciers are more productive in the summer, but have their peak activity in different months. Over the study period, Kangerdlugssuaq has had a constant number of events each year, whereas Jakobshavn had most events in 1998–1999, and the number of events in Helheim and the northwest Greenland glaciers has increased substantially between 1993 and 2005. The size distribution of events in Kangerdlugssuaq is peaked above the detection threshold, suggesting that glacial earthquakes have a characteristic size.