Search

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three‐dimensional simulations using the CVM‐S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM‐H do not match the timing of the initial arrivals as well as CVM‐S4.26.M01, and the strong late arrivals in the CVM‐H simulation travel with an apparent velocity slower than observed. A higher‐quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off‐azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long‐period ground motions in basins.

On 12 August 2000 two explosions damaged the Russian submarine, the Kursk. The largest event was well recorded at seismic networks in northern Europe, which we then modeled. We developed a hybrid method based on generalized ray theory that treats an explosive source embedded in a fluid and recorded along continental paths. Matching record sections of observations with synthetics, we obtain an estimate of explosive size of slightly over 4 t. Several earth models determined previously, K8 and a Baltic model, were used to assess accuracy. These results are in general agreement with other investigators using more empirical methods. Knowing the conventional missile yield and the explosion size allows for an estimate of approximately five missiles exploded in the second larger explosion onboard the Kursk.

The fine-scale seismic structure of the central Mexico subduction zone is studied using moderate-sized (M4-6) intraslab earthquakes. Regional waveforms from the Mapping the Rivera Subduction Zone (MARS) seismic array are complicated and contain detailed information about the subduction zone structure, including evidence of lateral heterogeneity. This waveform information is used to model the structure of the subducted plates, particularly along the transition from flat to normal subduction, where recent studies have shown evidence for possible slab tearing along the eastern projection of the Orozco Fracture Zone (OFZ). The lateral extent of a thin ultra-slow velocity layer (USL) imaged atop the Cocos slab in recent studies along the Meso America Subduction Experiment array is examined here using MARS waveforms. We find an edge to this USL which is coincident with the western boundary of the projected OFZ region. Forward modeling of the 2D structure of the subducted Rivera and Cocos plates using a finite difference algorithm provides constraints on the velocity and geometry of each slab's seismic structure in this region and confirms the location of the USL edge. We propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, in agreement with observations of variable Cocos plate motion on either side of the OFZ. This tearing event may be a young analogy to the 10 Ma Rivera-Cocos plate boundary, and may be related to the slab rollback in central Mexico.

Seismic tomography images of the upper mantle structures beneath the Pacific Northwestern United States display a maze of high-velocity anomalies, many of which produce distorted waveforms evident in the USArray observations indicative of the Juan de Fuca (JdF) slab. The inferred location of the slab agrees quite well with existing contour lines defining the slab's upper interface. Synthetic waveforms generated from a recent tomography image fit teleseismic travel times quite well and also some of the waveform distortions. Regional earthquake data, however, require substantial changes to the tomographic velocities. By modeling regional waveforms of the 2008 Nevada earthquake, we find that the uppermost mantle of the 1D reference model AK135, the reference velocity model used for most tomographic studies, is too fast for the western United States. Here, we replace AK135 with mT7, a modification of an older Basin-and-Range model T7. We present two hybrid velocity structures satisfying the waveform data based on modified tomographic images and conventional slab wisdom. We derive P and SH velocity structures down to 660 km along two cross sections through the JdF slab. Our results indicate that the JdF slab is subducted to a depth of 250 km beneath the Seattle region, and terminates at a shallower depth beneath Portland region of Oregon to the south. The slab is about 60 km thick and has a P velocity increase of 5% with respect to mT7. In order to fit waveform complexities of teleseismic Gulf of Mexico and South American events, a slab-like high-velocity anomaly with velocity increases of 3% for P and 7% for SH is inferred just above the 660 discontinuity beneath Nevada.

Upper mantle seismic velocity structures in both vertical and horizontal directions are key to understanding the structure and mechanics of tectonic plates. Recent deployment of the USArray Transportable Array (TA) in the Midwestern United States provides an extraordinary regional earthquake data set to investigate such velocity structure beneath the stable North American craton. In this paper, we choose an M_w5.1 Canadian earthquake in the Quebec area, which is recorded by about 400 TA stations, to examine the P wave structures between the depths of 150 km to 800 km. Three smaller Midwestern earthquakes at closer distance to the TA are used to investigate vertical and horizontal variations in P velocity between depths of 40 km to 150 km. We use a grid-search approach to find the best 1-D model, starting with the previously developed S25 regional model. The results support the existence of an 8° discontinuity in P arrivals caused by a negative velocity gradient in the lithosphere between depths of 40 km to 120 km followed by a small (∼1%) jump and then a positive gradient down to 165 km. The P velocity then decreases by 2% from 165 km to 200 km, and we define this zone as the regional lithosphere-asthenosphere boundary (LAB). Beneath northern profiles, waves reflected from the 410 discontinuity (410) are delayed by up to 1 s relative to those turning just below the 410, which we explain by an anomaly just above the discontinuity with P velocity reduced by ∼3%. The 660 discontinuity (660) appears to be composed of two smaller velocity steps with a separation of 16 km. The inferred low-velocity anomaly above 410 may indicate high water concentrations in the transition zone, and the complexity of the 660 may be related to Farallon slab segments that have yet to sink into the deep mantle.

We determined the location, size, mechanism, and the frequency content of the first 4.0 s of the 2011 Tohoku–Oki earthquake. Since the beginning of this earthquake is very small, we develop a comparative approach against a near-by reference earthquake, the master event. We first determined the water depth near the master event using the differential timing between the water phase pwP reflected from the air–water interface and the depth phase pP reflected from the water–crust interface. Then we located the master event using the well-known ocean bathymetry in the area. After calibrating teleseismic arrays (Δ=30° to 90°) at short periods for timing and amplitude with respect to the master event, we were able to determine the initiation of themain event. It began as a small (M_w=4.9) thrust event located at 38.19°N, 142.68°E at a depth of 21 km, and, a few seconds later, evolved into a slower extremely large slip event up-dip.

Theoretical studies on ambient seismic noise (ASN) predict that complete Green's function between seismic stations can be retrieved from cross correlation. However, only fundamental mode surface waves emerge in most studies involving real data. Here we show that Moho-reflected body wave (SmS) and its multiples can be identified with ASN for station pairs near their critical distances in the short period band (1–5 s). We also show that an uneven distribution of noise sources, such as mining activity and wind–topography interaction, can cause surface wave precursors, which mask weaker body wave phases.

We analyze the rupture process of the 1999 M_w 7.1 Duzce earthquake using seismological, remote sensing, and geodetic data. Ground deformation measured from the subpixel cross correlation of Satellite Pour l'Observation de la Terre (SPOT) images reveals a 55 km long fault trace and smooth surface-slip distribution peaking at 3.5–4 m. The westernmost segment overlaps for over 10 km with ruptures from the M_w 7.4 Izmit earthquake. The 15 km long easternmost segment, which cuts across mountainous topography, had not been reported previously. We determine a well-constrained source model using a four-segment fault geometry using constraints on surface fault slip and inverting Global Positioning System and Interferometric Synthetic Aperture Radar data along with strong-motion records. Our results show that some variability of the rupture velocity and an eastward supershear velocity are required to fit the strong-motion data. The rise time, up to 6 sec, correlates with cumulative slip, suggesting a sliding velocity of about 1 m/sec. The source model predicts teleseismic waveforms well, although early by 2 sec. This time shift is probably due to the weak beginning of the earthquake that is not observable at teleseismic distances. Strong-motion records are relatively well predicted from a source model derived from the teleseismic data using the fault geometry derived from the satellite images. This study demonstrates the benefit of using accurate fault geometries to determine finite-fault source models.

A recent prospective on the Yellowstone Caldera discounts its explosive potential based on inferences from tomographic studies which suggests a high degree of crystallization of the underlying magma body. In this study, we show that many of the first teleseismic P-wave arrivals observed at seismic stations on the edge of the caldera did not travel through the magma body but have taken longer but faster paths around the edge. After applying a number of waveform modeling tools, we obtain much lower seismic velocities than previous studies, 2.3 km/sec (V_p) and 1.1 km/sec (V_s). We estimate the physical state of the magma body by assuming a fluid-saturated porous material consisting of granite and a mixture of rhyolite melt and water and CO_2 at a temperature of 800°C and pressure at 5 km (0.1 GPa). We found that this relatively shallow magma body has a volume of over 4,300 km^3 and is about 32% melt saturated with about 8% water plus CO_2 by volume.

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.

We analyze the M_w 7.6 Kashmir earthquake of October 8, 2005, using sub-pixel correlation of ASTER images to measure ground deformation, and modeling seismic waveforms. The surface rupture is continuous over a distance of 75 km and cuts across the Hazara syntaxis reactivating the Tanda and the Muzaffarabad faults. North of Muzaffarabad the surface rupture coincides approximately with the MBT, on the southwestern flank of the syntaxis, although the two faults have opposite dip angles. The rupture terminates abruptly at the hairpin turn of the MBT showing a strong structural control. The fault offset is 4 m on average and peaks to 7 m northwest of Muzaffarabad. The rupture lasted about 25 s and propagated updip and bi-laterally ~2 km/s, with a rise time of 2–5 s. The shallowness and compactness of the rupture, both in time and space, provide an explanation for the intensity of destructions. This kind of analysis could be achieved as soon as a post-earthquake image is available, and would provide key information for early assessment of damages. The study sheds some light on seismic hazard in the Himalaya, and raises concern regarding the possibility of a repetition of the 1555 event which presumably ruptured the Himalayan front south of the Kashmir basin and may have reached a magnitude M_w > 8.

Recent (ScS-S) results from probing the deep structure beneath southern Africa display strong delays of up to 10 s at distances beyond 90°. Such delays could be explained by long-period tomographic models containing smooth (weak) features with the addition of rough (strong) D″ structure (3–9% drops in shear velocities). However, these structures cannot explain the (SKS-S) differentials sampling the same region. To explain the (SKS-S) and (ScS-S) data sets simultaneously requires instead a large-scale ridge-like structure with a relatively uniform 3% reduction in shear velocity. The structure is about 1000 km wide and extends at least 1200 km above D″. It is orientated roughly NW-SE and leans toward the east at latitudes from 15° to about 30°. It proves difficult to explain such sharp features with thermal effects alone and, thus, the importance of high-resolution waveform modeling to establish their existence. To derive the above results, we developed a special algorithm by matching simulated synthetics to observed broadband waveforms. This is achieved by computing the various arrivals separately using generalized ray theory for a reference model and allowing the arrivals to shift in relative times to match the data. Tomographic models can then be constructed, or existing tomographic models can be altered, to match these data, and new 2-D synthetics can be constructed as well to better fit the waveform data. These updated synthetics can again be decomposed and reassembled, and the process can be repeated. This algorithm is applied to a combination of analog and digital data along a corridor from South America, producing the high-resolution 2-D model described above.