Teh-Ru Alex Song

Low-velocity zone atop the 410-km seismic discontinuity in the northwestern United States

The seismic discontinuity at 410 km depth in the Earths mantle is generally attributed to the phase transition of (Mg,Fe)2SiO4 (refs 1, 2) from the olivine to wadsleyite structure. Variation in the depth of this discontinuity is often taken as a proxy for mantle temperature owing to its response to thermal perturbations. For example, a cold anomaly would elevate the 410-km discontinuity, because of its positive Clapeyron slope, whereas a warm anomaly would depress the discontinuity. But trade-offs between seismic wave-speed heterogeneity and discontinuity topography often inhibit detailed analysis of these discontinuities, and structure often appears very complicated. Here we simultaneously model seismic refracted waves and scattered waves from the 410-km discontinuity in the western United States to constrain structure in the region. We find a low-velocity zone, with a shear-wave velocity drop of 5%, on top of the 410-km discontinuity beneath the northwestern United States, extending from southwestern Oregon to the northern Basin and Range province. This low-velocity zone has a thickness that varies from 20 to 90 km with rapid lateral variations. Its spatial extent coincides with both an anomalous composition of overlying volcanism and seismic ‘receiver-function’ observations observed above the region. We interpret the low-velocity zone as a compositional anomaly, possibly due to a dense partial-melt layer, which may be linked to prior subduction of the Farallon plate and back-arc extension. The existence of such a layer could be indicative of high water content in the Earths transition zone.

Subducting slab ultra-slow velocity layer coincident with silent earthquakes in southern Mexico.

Seismic mapping suggests that silent earthquakes may be related to an ultralow velocity layer on top of a subducting slab. Hot Silent Quakes Subduction zones tend to produce the largest and potentially most destructive earthquakes. Recent observations show that some deformation in several subduction zones seems to be occurring through small or “silent” quakes. The origin of these silent quakes, and their effect on the seismic hazard, is uncertain. Song et al. (p. 502) use a specific seismic signal to map out thin regions with low seismic velocities on the subduction zone beneath southern Mexico. The regions seem to occur at depths below the seismogenic zone where temperatures are higher. These high temperatures and the silent quakes may reflect the release and episodic trapping of fluids from metamorphic reactions. Great earthquakes have repeatedly occurred on the plate interface in a few shallow-dipping subduction zones where the subducting and overriding plates are strongly locked. Silent earthquakes (or slow slip events) were recently discovered at the down-dip extension of the locked zone and interact with the earthquake cycle. Here, we show that locally observed converted SP arrivals and teleseismic underside reflections that sample the top of the subducting plate in southern Mexico reveal that the ultra-slow velocity layer (USL) varies spatially (3 to 5 kilometers, with an S-wave velocity of ~2.0 to 2.7 kilometers per second). Most slow slip patches coincide with the presence of the USL, and they are bounded by the absence of the USL. The extent of the USL delineates the zone of transitional frictional behavior.

Rupture kinematics of the 2005 Mw 8.6 Nias-Simeulue earthquake from the joint inversion of seismic and geodetic data

The 2005 Mw 8.6 Nias-Simeulue earthquake was caused by rupture of a portion of the Sunda megathrust offshore northern Sumatra. This event occurred within an array of continuous Global Positioning System (GPS) stations and produced measurable vertical displacement of the fringing coral reefs above the fault rupture. Thus, this earthquake provides a unique opportunity to assess the source character- istics of a megathrust event from the joint analysis of seismic data and near-field static co-seismic displacements. Based on the excitation of the normal mode data and geodetic data we put relatively tight constraints on the seismic moment and the fault dip, where the dip is determined to be 8� to 10� with corresponding moments of 1.24 10 22 to 1.00 10 22 N m, respectively. The geodetic constraints on slip distribution help to eliminate the trade-off between rupture velocity and slip kine- matics. Source models obtained from the inversion of various combinations of the teleseismic body waves and geodetic data are evaluated by comparing predicted and observed long-period seismic waveforms (100-500 sec). Our results indicate a rela- tively slow average rupture velocity of 1.5 to 2.5 km/sec and long average rise time of up to 20 sec. The earthquake nucleated between two separate slip patches, one beneath Nias and the other beneath Simeulue Island. The gap between the two patches and the hypocentral location appears to be coincident with a local geological disrup- tion of the forearc. Coseismic slip clearly tapers to zero before it reaches the trench probably because the rupture propagation was inhibited when it reached the accre- tionary prism. Using the models from joint inversions, we estimate the peak ground velocity on Nias Island to be about 30 cm/sec, an order of magnitude slower than for thrust events in continental areas. This study emphasizes the importance of util- izing multiple datasets in imaging seismic ruptures.

Evidence for fault lubrication during the 1999 Chi-Chi, Taiwan, earthquake (Mw7.6)

The ground motion data of the 1999 Chi-Chi, Taiwan, earthquake exhibit a striking difference in frequency content between the north and south portions of the rupture zone. In the north, the ground motion is dominated by large low-frequency displacements with relatively small high-frequency accelerations. The pattern is opposite in the south, with smaller displacements and larger accelerations. We analyze the fault dynamics in light of a fault lubrication mechanism using near-field seismograms and a detailed rupture model. The fault zone contains viscous material (e.g., gouge), in which pressure increases following the Reynolds lubrication equation. When the displacement exceeds a threshold, lubrication pressure becomes high enough to widen the gap, thereby reducing the area of asperity contact. With less asperity contact, the fault slips more smoothly, suppressing high-frequency radiation.

Spatial slip distribution of the September 20, 1999, Chi‐Chi, Taiwan, Earthquake (MW7.6) —Inverted from teleseismic data

The teleseismic waveforms of the MW7.6 September 20, 1999 Chi-Chi earthquake were examined to obtain the quick information on the fault rupture process. The deconvolution results show the fault ruptured from south to the north, and revealed the west movement of the hanging wall to the footwall on the eastern dipping plane. The spatial slip distribution shows that the earthquake was mainly composed by a large asperity with a dimension of about 45km × 15km. The maximum slip was about 8 m located at about 45 km to the north of the epicenter. The slip distributions obtained in this study have a good agreement with the observed surface breaks. The comparable static and the dynamic stress drops of about 11 Mpa in the large slip region indicate that the melting/fluid pressurization might have taken place in this earthquake. It reduces the dynamic friction and results in the large slip observed.

Slip distribution and tectonic implication of the 1999 Chi-Chi, Taiwan, earthquake

We report on the fault complexity of the large (M_w = 7.6) Chi‐Chi earthquake obtained by inverting densely and well‐distributed static measurements consisting of 119 GPS and 23 doubly integrated strong motion records. We show that the slip of the Chi-Chi earthquake was concentrated on the surface of a ”wedge shaped” block. The inferred geometric complexity explains the difference between the strike of the fault plane determined by long period seismic data and surface break observations. When combined with other geophysical and geological observations, the result provides a unique snapshot of tectonic deformation taking place in the form of very large (>10m) displacements of a massive wedge‐shaped crustal block which may relate to the changeover from over‐thrusting to subducting motion between the Philippine Sea and the Eurasian plates.

Two‐stage composite megathrust rupture of the 2015 Mw8.4 Illapel, Chile, earthquake identified by spectral‐element inversion of teleseismic waves

The Mw8.4 Illapel earthquake occurred on 16 September was the largest global event in 2015. This earthquake was not unexpected because the hypocenter was located in a seismic gap of the Peru-Chile subduction zone. However, the source model derived from 3-D spectral-element inversion of teleseismic waves reveals a distinct two-stage rupture process with completely different slip characteristics as a composite megathrust event. The two stages were temporally separated. Rupture in the first stage, with a moment magnitude of Mw8.32, built up energetically from the deeper locked zone and propagated in the updip direction toward the trench. Subsequently, the rupture of the second stage, with a magnitude of Mw8.08, mainly occurred in the shallow subduction zone with atypical repeating slip behavior. The unique spatial-temporal rupture evolution presented in this source model is key to further in-depth studies of earthquake physics and source dynamics in subduction systems.

Low Velocity Zone Atop the Transition Zone in the Western US From S Waveform Triplication

Song et al. [2004] modeled regional S wave triplications in the northwestern US and found a low velocity zone atop the 410 seismic discontinuity. Strong azimuthal variation in waveforms associated with paths sampling the western edge of this structure are observed on the TriNet array for several events. Here, we model this data with a new 3D simulation technique which combines 2D finite-difference with Kirchhoff diffraction operators to include responses off the great circle. To reconcile such sharp changes in waveforms requires a sharp western edge less than 100 km across a boundary with a change of 3-5% in velocity. Though the geometry of the LVZ is not unique due to limited data analyzed, the sharp edge of the LVZ is robustly constrained with available array data. Such a LVZ is consistent with the existence of water in Earths transition zone, at least locally.

Source Characteristics of the 2016 Meinong (ML 6.6), Taiwan, Earthquake, Revealed from Dense Seismic Arrays: Double Sources and Pulse‐like Velocity Ground MotionSource Characteristics of the 2016 Meinong, Taiwan, Earthquake

The 5 February 2016, Meinong, Taiwan, earthquake brought extensive damage to nearby cities with significant pulse‐like velocity ground motions. In addition to the spatial slip distribution determination using filtered strong‐motion data, we show that, with the advantage of the densely distributed seismic network as a seismic array, we can project the earthquake sources (asperities) directly using nearly unfiltered data, which is crucial to the understanding of the generation of the pulse‐like velocity ground motions. We recognize that the moderate but damaging M_L 6.6 Meinong earthquake was a composite of an M_w 5.5 foreshock and an M_w 6.18 mainshock with a 1.8–5.0 s time delay. The foreshock occurred at the hypocenter reported by the official agency, followed by the mainshock with a centroid located at 12.3 km to the north‐northwest of the hypocenter and at a depth of 15 km. This foreshock–mainshock composition is not distinguishable in the finite‐fault inversion because it filtered the seismic data to low frequencies. Our results show that the pulse‐like velocity ground motions are mainly attributed to the source of mainshock with its directivity and site effects, resulting in the disastrous damages in the city of Tainan. Although finite‐fault inversion using filtered seismic data for spatial slip distribution on the fault has been a classic procedure in understanding earthquake rupture processes, using a dense seismic network as a seismic array for unfiltered records helps us delineate the earthquake sources directly and provide more delicate information for future understanding of earthquake source complexity.