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.

The first cycle (ab amplitude) of teleseismic short-period P waves from underground nuclear explosions at Pahute Mesa (NTS) show a systematic azimuthal amplitude pattern that can possibly be explained by tectonic release. The amplitudes vary by a factor of three, with diminished amplitudes being recorded at azimuths around N25°E. This azimuthal pattern has a strong sin(2φ) component and is observed, to varying degrees, for 25 Pahute Mesa events, but not for events at other sites within the NTS. Events that are known to have large tectonic release have more pronounced sin(2φ) amplitude variations. A synthesis of long-period body and surface wave investigations of tectonic release for Pahute Mesa events shows that, in general, the nonisotropic radiation is equivalent to nearly vertical, right-lateral strike-slip faulting trending from N20°W to due north. Long-period P waves at upper mantle distances demonstrate that there is a significant high-frequency component to the tectonic release. Using the long-period constraints on orientation, moment, and frequency content of the tectonic release, the expected short-period P wave effects are predicted. For models in which the downgoing P wave from the explosion triggers tectonic release within a few kilometers below the shot point, a factor of 2.5 amplitude variation with azimuth is predicted for the short-period ab amplitudes, with the lowest amplitudes expected near N25°E. Rather subtle azimuthal variations in the waveforms are expected, particularly for downward propagating ruptures, which is consistent with the absence of strong variations in the data. The occurrence of the azimuthal pattern, albeit with varying strength, for all of the Pahute Mesa events suggests a tectonic release model in which the shatterzone surrounding the explosion cavity is extended preferentially downward by driving a distributed network of faults and joints underlying the Mesa several kilometers beneath the surface. In this model, all events could have a component of tectonic release which would reflect the regional stress regime, although there may be slight spatial and temporal variations in the tectonic release contribution. Some events may trigger slip on larger throughgoing faults as well. While it is shown that tectonic release can affect teleseismic short-period signals significantly, and may contribute to the Pahute Mesa amplitude pattern, other possible explanations are considered.

The Santa Barbara earthquake of 13 August 1978, provides an opportunity to perform a broadband investigation of body waves for a well-recorded, moderate size (M_L = 5.1) event. The long- and short-period teleseismic body waves are modeled in the time domain to construct a source time function which is consistent in the period range of 1 to 20 sec. The long-period records indicate an overall duration of 6 sec while the short-period records reveal the fine-scale character of the slip history consisting of two sharp pulses separated by about 1 sec. The source mechanism determined from this analysis is a moderately dipping (30°NE) thrust with significant left-lateral slip. The moment was determined to be 1.1 × 10^(25) dyne-cm. The earthquake was also reasonably well recorded on accelerographs in the near-field. The modeling of the strong motion displacements was a two step procedure: (1) the displacements were modeled alone, and (2) in an attempt to achieve consistency between the local and far-field time functions, the qualitative features of the teleseismic short-period time function were used to predict the displacements. If the two sources in the short-period time function are allowed to have different mechanisms, the displacements can be modeled quite well. This suggests that the overall faulting process was rough, and the multiple source character suggested at high frequencies is due to high-stress drop asperities. The two sources are modeled as asperities separated by 1.5 km; the first source has a mechanism consistent with the teleseismic solution while the second source is more steeply dipping. The total moment determined from the strong motion data is 3.5 × 10^(24) dyne-cm or one-third the long-period moment. This is consistent with other recent studies which suggest that the high-frequency strong ground motion is controlled by the distribution of asperities even though the sum of their moments may be small compared to the overall moment. This study also shows the importance of teleseismic short periods in predicting the local displacements.

Source models are determined for the three underground nuclear explosions at the Amchitka test site using seismic observations in the period range 0.5 to 20.0 sec. Empirical yield-scaling relations are inferred from the source models and compared with the predictions of the Haskell, Mueller-Murphy, and finite difference numerical models. Several recent studies of high-frequency, near-field signals and teleseismic short-period P waves for LONGSHOT, MILROW, and CANNIKIN constrain the source functions at periods of 0.5 to 2.0 sec. Teleseismic pS and Rayleigh wave observations are used to constrain the source functions at longer periods. Using a modified Haskell source time function representation given by ψ(t) = ψ_∞ {1-e^(-kt)[1 + Kt + (Kt)^2/2-B(Kt)^3]}, the data are best-fit if the corner frequency parameter, K, scales as predicted by the Mueller-Murphy model, and if the amount of overshoot in the reduced displacement potential, which is proportional to B, decreases with increasing yield (depth of burial). The latter behavior is opposite to that predicted by the Mueller-Murphy model and follows from the observation that the long-period level of the explosion potential, ψ_∝, increases with yield, W, by ψ_∝ ∝ W^(0.90), or with yield and depth by ψ_∝ ∝ W/h^(1/3). This long-period and overshoot scaling is consistent with that found for some numerical models, and allowing for the depth dependence of the Rayleigh wave excitation, results in the observed M_S versus log(W) slope of ∼1. The decrease in overshoot with increasing depth of burial may be the result of the increase in shear strength with increasing overburden pressure. If yield or depth dependence of the source potential overshoot proves to be a general phenomenon, a possibility supported by a preliminary investigation of Pahute Mesa observations, accurate yield estimation will require broadband seismic data. The source function representation adopted is shown to provide an excellent fit to the rise time of very near-in velocity recordings to the rise time with frequencies of 10 Hz and higher.