Validating tomographic model with broad-band waveform modelling: an example from the LA RISTRA transect in the southwestern United States

Traveltime tomographic models of the LA RISTRA transect produce excellent waveform fits if we amplify the damped images. We observe systematic waveform distortions across the western edge of the Great Plains from South American events, starting about 300 km east of the centre of the Rio Grande Rift. The amplitude decreases by more than 50 per cent within array stations spanning less than 200 km while the pulse width increases by more than a factor of 2. This feature is not observed for the data arriving from the northwest. While the S-wave tomographic image shows a fast slab-like feature dipping to the southeast beneath the western edge of the Great Plains, synthetics generated from this model do not reproduce the waveform characteristics. However, once we modify the tomographic image by amplifying the velocity contrast between the slab and adjoining mantle by a factor of 2–3, the synthetics produce observed amplitude decay and pulse broadening. In addition to the traveltime delay, amplitude variation due to wave phenomena such as slab diffraction, focusing and defocusing provide much tighter constraints on the geometry of the fast anomaly and its amplitude and sharpness as demonstrated by a forward sensitivity test and snapshots of the seismic wavefield. Our preferred model locates the slab 200 km east of the Rio Grande Rift dipping 70°–75° to the southeast, extending to a depth near 600 km with a thickness of 120 km and a velocity of about 4 per cent fast. In short, adding waveform and amplitude components to regional tomographic studies can help validate and establish structural geometry, sharpness and velocity contrast

Large Trench-Parallel Gravity Variations Predict Seismogenic Behavior in Subduction Zones

We demonstrate that great earthquakes occur predominantly in regions with a strongly negative trench-parallel gravity anomaly (TPGA), whereas regions with strongly positive TPGA are relatively aseismic. These observations suggest that, over time scales up to at least 1 million years, spatial variations of seismogenic behavior within a given subduction zone are stationary and linked to the geological structure of the fore-arc. The correlations we observe are consistent with a model in which spatial variations in frictional properties on the plate interface control trench-parellel variations in fore-arc topography, gravity, and seismogenic behavior

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

The seismic discontinuity at 410 km depth in the Earth's mantle is generally attributed to the phase transition of (Mg,Fe)_2SiO_4 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 Earth's transition zone

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

Validation of the rupture properties of the 2001 Kunlun, China (M_s=8.1), earthquake from seismological and geological observations

We determine the finite-fault slip distribution of the 2001 Kunlun earthquake (M_s = 8.1) by inverting teleseismic waveforms, as constrained by geological and remote sensing field observations. The spatial slip distribution along the 400-km-long fault was divided into five segments in accordance with geological observations. Forward modelling of regional surface waves was performed to estimate the variation of the speed of rupture propagation during faulting. For our modelling, the regional 1-D velocity structure was carefully constructed for each of six regional seismic stations using three events with magnitudes of 5.1–5.4 distributed along the ruptured portion of the Kunlun fault. Our result shows that the average rupture velocity is about 3.6 km s^−1, consistent with teleseismic long period wave modelling. The initial rupture was almost purely strike-slip with a rupture velocity of 1.9 km s^−1, increasing to 3.5 km s^−1 in the second fault segment, and reaching a rupture velocity of about 6 km s^−1 in the third segment and the fourth segment, where the maximum surface offset, with a broad fault zone, was observed. The rupture velocity decelerated to a value of 3.3 km s^−1 in the fifth and final segment. Coseismic slip on the fault was concentrated between the surface and a depth of about 10 km. We infer that significant variations in rupture velocity and the observed fault segmentation are indicative of variations in strength along the interface of the Kunlun fault, as well as variations in fault geometry

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

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

Observation of temporal variations in seismic anisotropy within an active fault‐zone revealed from the Taiwan Chelungpu‐fault Drilling Project Borehole Seismic Array

Temporal fault-zone observations are important to better understand the evolution of fault structure and stress configuration. However, long-term monitoring in the fault-zone is rare after a large earthquake. Here, we use seismic data in the fault-zone at 1-km depth from the Taiwan Chelungpu-fault Drilling Project to study long-term anisotropy after the 1999 Mw7.6 Chi-Chi earthquake. The direct S-wave splitting measurements resolve the overall weak anisotropy in the shallow crust. In order to resolve fault damage zone anisotropy, we perform coda cross-correlation technique for 794 microearthquakes between 2007 and 2013. We estimate the temporal change in background shear-wave velocity, fast shear-wave polarization direction (FSP), and strength of anisotropy (Aani) in the fault damage zone. We show the average FSP direction is N93°E with a significant Aani of about 12%, likely due to the pervasive vertical microcracks created after the earthquake. Temporal variations of anisotropy exhibit seasonal variation with periodicity every 9 to 12 months that correlates with rainfall events. Furthermore, long-term anisotropy shows a gradual rotation of FSP direction of about 15° during the first 4 years of observation. At the same time, the strength of anisotropy reduced from 17 to 10 % and shear-wave velocity increased, suggesting the fault healed after the earthquake. This study reports in-situ evidence for two key observations: (1) long-term, fault-zone healing after a major earthquake, and (2) modulation of 1-km deep fault-zone properties by surficial hydrologic processes. These observations may provide constraints on the response of the fault damage zone in the interseismic period

Rupture Kinematics of the 2005 M_w 8.6 Nias–Simeulue Earthquake from the Joint Inversion of Seismic and Geodetic Data

The 2005 M_w 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 characteristics 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 x 10^(22) to 1.00 x 10^(22) N m, respectively. The geodetic constraints on slip distribution help to eliminate the trade-off between rupture velocity and slip kinematics. 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 relatively 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 disruption 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 accretionary 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 utilizing multiple datasets in imaging seismic ruptures

Evidence for non-self-similarity of microearthquakes recorded at a Taiwan borehole seismometer array

We investigate the relationship between seismic moment M_0 and source duration t_w of microearthquakes by using high-quality seismic data recorded with a vertical borehole array installed in central Taiwan. We apply a waveform cross-correlation method to the three-component records and identify several event clusters with high waveform similarity, with event magnitudes ranging from 0.3 to 2.0. Three clusters—Clusters A, B and C—contain 11, 8 and 6 events with similar waveforms, respectively. To determine how M_0 scales with t_w, we remove path effects by using a path-averaged Q. The results indicate a nearly constant t_w for events within each cluster, regardless of M_0, with mean values of t_w being 0.058, 0.056 and 0.034 s for Clusters A, B and C, respectively. Constant t_w, independent of M_0, violates the commonly used scaling relation t_w ∝ M^(1/3)_0. This constant duration may arise either because all events in a cluster are hosted on the same isolated seismogenic patch, or because the events are driven by external factors of constant duration, such as fluid injections into the fault zone. It may also be related to the earthquake nucleation size

Coseismic Slip and Afterslip of the Great M_w 9.15 Sumatra–Andaman Earthquake of 2004

We determine coseismic and the first-month postseismic deformation associated with the Sumatra–Andaman earthquake of 26 December 2004 from near- field Global Positioning System (GPS) surveys in northwestern Sumatra and along the Nicobar-Andaman islands, continuous and campaign GPS measurements from Thailand and Malaysia, and in situ and remotely sensed observations of the vertical motion of coral reefs. The coseismic model shows that the Sunda subduction megathrust ruptured over a distance of about 1500 km and a width of less than 150 km, releasing a total moment of 6.7–7.0 x 10^(22) N m, equivalent to a magnitude M_w 9.15. The latitudinal distribution of released moment in our model has three distinct peaks at about 4° N, 7° N, and 9° N, which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated T waves. Our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very-long-period surface waves. The tsunami predicted from this model fits relatively well the altimetric measurements made by the JASON and TOPEX satellites. Neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave. The near-field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500-sec-long seismic rupture. The postseismic geodetic moment of about 2.4 x 10^(22) N m (M_w 8.8) is equal to about 30 ± 5% of the coseismic moment release. Evolution of postseismic deformation is consistent with rate-strengthening frictional afterslip