Search

We have determined a source rupture model for the 1992 Landers earthquake (MW 7.2) compatible with multiple data sets, spanning a frequency range from zero to 0.5 Hz. Geodetic survey displacements, near-field and regional strong motions, broadband teleseismic waveforms, and surface offset measurements have been used explicitly to constrain both the spatial and temporal slip variations along the model fault surface. Our fault parameterization involves a variable-slip, multiple-segment, finite-fault model which treats the diverse data sets in a self-consistent manner, allowing them to be inverted both independently and in unison. The high-quality data available for the Landers earthquake provide an unprecedented opportunity for direct comparison of rupture models determined from independent data sets that sample both a wide frequency range and a diverse spatial station orientation with respect to the earthquake slip and radiation pattern. In all models, consistent features include the following: (1) similar overall dislocation patterns and amplitudes with seismic moments of 7 to 8 × 1026 dyne-cm (seismic potency of 2.3 to 2.7 km3); (2) very heterogeneous, unilateral strike slip distributed over a fault length of 65 km and over a width of at least 15 km, though slip is limited to shallower regions in some areas; (3) a total rupture duration of 24 sec and an average rupture velocity of 2.7 km/sec; and (4) substantial variations of slip with depth relative to measured surface offsets. The extended rupture length and duration of the Landers earthquake also allowed imaging of the propagating rupture front with better resolution than for those of prior shorter-duration, strike-slip events. Our imaging allows visualization of the rupture evolution, including local differences in slip durations and variations in rupture velocity. Rupture velocity decreases markedly at shallow depths, as well as near regions of slip transfer from one fault segment to the next, as rupture propagates northwestward along the multiply segmented fault length. The rupture front slows as it reaches the northern limit of the Johnson Valley/Landers faults where slip is transferred to the southern Homestead Valley fault; an abrupt acceleration is apparent following the transfer. This process is repeated, and is more pronounced, as slip is again passed from the northern Homestead Valley fault to the Emerson fault. Although the largest surface offsets were observed at the northern end of the rupture, our modeling indicates that substantial rupture was also relatively shallow (less than 10 km) in this region.

The U.S. Geological Survey (USGS) “Did You Feel It?” (DYFI) system is an automated approach for rapidly collecting macroseismic intensity data from Internet users’ shaking and damage reports and generating intensity maps immediately following earthquakes; it has been operating for over a decade (1999-2011). DYFI-based intensity maps made rapidly available through the DYFI system fundamentally depart from more traditional maps made available in the past. The maps are made more quickly, provide more complete coverage and higher resolution, provide for citizen input and interaction, and allow data collection at rates and quantities never before considered. These aspects of Internet data collection, in turn, allow for data analyses, graphics, and ways to communicate with the public, opportunities not possible with traditional data-collection approaches. Yet web-based contributions also pose considerable challenges, as discussed herein. After a decade of operational experience with the DYFI system and users, we document refinements to the processing and algorithmic procedures since DYFI was first conceived. We also describe a number of automatic post-processing tools, operations, applications, and research directions, all of which utilize the extensive DYFI intensity datasets now gathered in near-real time. DYFI can be found online at the website http://earthquake.usgs.gov/dyfi/.

All quality teleseismic recordings of the great 1906 San Francisco earthquake archived in the 1908 Carnegie Report by the State Earthquake Investigation Commission were scanned and digitized. First order results were obtained by comparing complexity and amplitudes of teleseismic waveforms from the 1906 earthquake with well calibrated, similarly located, more recent earthquakes (1979 Coyote Lake, 1984 Morgan Hill, and 1989 Loma Prieta earthquakes) at nearly co-located modern stations. Peak amplitude ratios for calibration events indicated that a localized moment release of about 1 to 1.5 × 1027 dyne-cm was responsible for producing the peak the teleseismic body wave arrivals. At longer periods (50 to 80 sec), we found spectral amplitude ratios of the surface waves require a total moment release between 4 and 6 × 1027 dyne-cm for the 1906 earthquake, comparable to previous geodetic and surface wave estimates (Thatcher, 1975). We then made a more detailed source analysis using Morgan Hill S body waves as empirical Green's Functions in a finite fault subevent summation. The Morgan Hill earthquake was deemed most appropriate for this purpose as its mechanism is that of the 1906 earthquake in the central portion of the rupture. From forward and inverse empirical summations of Morgan Hill Green's functions, we obtained a good fit to the best quality teleseismic waveforms with a relatively simple source model having two regions of localized strong radiation separated spatially by about 110 km. Assuming the 1906 epicenter determined by Bolt (1968), this corresponds with a large asperity (on the order of the Loma Prieta earthquake) in the Golden Gate/San Francisco region and one about three times larger located northwest along strike between Point Reyes and Fort Ross. This model implies that much of the 1906 rupture zone may have occurred with relatively little 10 to 20 sec radiation. Consideration of the amplitude and frequency content of the 1906 teleseismic data allowed us to estimate the scale length of the largest asperity to be less than about 40 km. With rough constraints on the largest asperity (size and magnitude) we produced a suite of estimated synthetic ground velocities assuming a slip distribution similar to that of the Loma Prieta earthquake but with three times as much slip. For purposes of comparison with the recent, abundant Loma Prieta strong motion data set, we “moved” the largest 1906 asperity into Loma Prieta region. Peak ground velocity amplitudes are substantially greater than those recorded during the Loma Prieta earthquake, and are comparable to those predicted by the attenuation relationship of Joyner and Boore (1988) for a magnitude MW = 7.7 earthquake.

A little more than 50 years ago, on 27 March 1964, the Great Alaska earthquake and tsunami struck. At moment magnitude 9.2, this earthquake is notable as the largest in U.S. written history and as the second-largest ever recorded by instruments worldwide. But what resonates today are its impacts on the understanding of plate tectonics, tsunami generation, and earthquake history as well as on the development of national programs to reduce risk from earthquakes and tsunamis.

The advent of the U.S. Geological Survey (USGS) Prompt Assessment of Global Earthquakes for Response (PAGER) system, in conjunction with several recent advances and trends in related data sources and research efforts, bring to light new opportunities within the overlapping realms of earthquake hazard, earthquake engineering, and earthquake epidemiological studies. While casualty modelling has admittedly often suffered from the lack of epidemiological rigour on the part of earth scientists and engineers, comparable laxity is also evident in some analyses of related hazard complexities on the part of social scientists. These limitations have often been due to insufficient oversight or interaction, or more commonly, insufficient data availability. Thanks to improved data sets, modelling approaches, and collaborations, there are now fewer obstacles to performing comprehensive casualty estimation, though formidable challenges remain. Under the auspices of the PAGER system, a global set of ShakeMaps has been produced for all significant earthquakes in the past 34 years (1973–2007). These event-specific ShakeMaps, constrained by any available data, are then combined with new global population data sets to develop systematic hazard and loss analyses. These and other important advancements, as well as their limitations, and their potential for contributing to casualty modelling are discussed. Example studies and applications are presented.

Since the launch of the USGS’s Prompt Assessment of Global Earthquakes for Response (PAGER) system in fall of 2007, the time needed for the U.S. Geological Survey (USGS) to determine and comprehend the scope of any major earthquake disaster anywhere in the world has been dramatically reduced to less than 30 min. PAGER alerts consist of estimated shaking hazard from the ShakeMap system, estimates of population exposure at various shaking intensities, and a list of the most severely shaken cities in the epicentral area. These estimates help government, scientific, and relief agencies to guide their responses in the immediate aftermath of a significant earthquake. To account for wide variability and uncertainty associated with inventory, structural vulnerability and casualty data, PAGER employs three different global earthquake fatality/loss computation models. This article describes the development of the models and demonstrates the loss estimation capability for earthquakes that have occurred since 2007. The empirical model relies on country-specific earthquake loss data from past earthquakes and makes use of calibrated casualty rates for future prediction. The semi-empirical and analytical models are engineering-based and rely on complex datasets including building inventories, time-dependent population distributions within different occupancies, the vulnerability of regional building stocks, and casualty rates given structural collapse.

We thank Art Frankel for his continuing interest in unraveling the details of the complicated source process of the 1987 Superstition Hills earthquake. His insightful comments and subsequent modeling on our part has helped to resolve some of the differences between the rupture model presented in Wald et al. (1990), hereafter referred to as WEA, and that published by Frankel and Wennerberg (1989), hereafter F&W.

Following a major earthquake, one of Caltrans' most critical tasks is to assess the condition of all potentially impacted bridges and roadway corridors in the state highway system. Timely response is important to ensure public safety, guide emergency vehicle traffic, and re-establish critical lifeline routes. In the past, bridge inspection teams had difficulties setting priorities immediately following an earthquake, because they lacked precise information about where the worst shaking and, most likely, the greatest damage had occurred. In 2005 Caltrans initiated a contract with the United States Geological Survey (USGS) to develop and implement a Caltrans-specific version of their ShakeCast system, a post-event software analysis tool that was anticipated to change the way Caltrans responds following a major earthquake. ShakeCast is a web-based application that automatically retrieves measured earthquake shaking data, analyzes it against individual bridge performance characteristics, and generates inspection prioritization emails and other web-based products for emergency responders within minutes of an event. Benefits have already been realized over the course of the test deployment phase. The ShakeCast system has proved to be a valuable tool for post-earthquake response for Caltrans during real events as well as for scenario planning exercises.

[1] Ongoing developments in earthquake source inversions incorporate nonplanar fault geometries as inputs to the inversion process, improving previous approaches that relied solely on planar fault surfaces. This evolution motivates advancing the existing framework for constraining fault geometry, particularly in subduction zones where plate boundary surfaces that host highly hazardous earthquakes are clearly nonplanar. Here, we improve upon the existing framework for the constraint of the seismic rupture plane of subduction interfaces by incorporating active seismic and seafloor sediment thickness data with existing independent data sets and inverting for the most probable nonplanar subduction geometry. Constraining the rupture interface a priori with independent geological and seismological information reduces the uncertainty in the derived earthquake source inversion parameters over models that rely on simpler assumptions, such as the moment tensor inferred fault plane. Examples are shown for a number of wellconstrained global locations. We expand the coverage of previous analyses to a more uniform global data set and show that even in areas of sparse data this approach is able to accurately constrain the approximate subduction geometry, particularly when aided with the addition of data from local active seismic surveys. In addition, we show an example of the integration of many two-dimensional profiles into a threedimensional surface for the Sunda subduction zone and introduce the development of a new global threedimensional subduction interface model: Slab1.0. Components: 8122 words, 8 figures, 1 table.

The U. S. Geological Survey (USGS) ShakeCast system is a freely-available, post-earthquake situational awareness application that automatically retrieves earthquake shaking data from ShakeMap, compares intensity measures against users' facilities, sends notifications of potential damage to responsible parties, and generates facility damage assessment maps and other Web-based products for emergency managers and responders. The ShakeMap/ShakeCast combined system allows a consistent approach for the evaluation of facility performance using the ShakeMap methodology by combining observations with ground motion predictions and covers major geographic regions and earthquakes, real-time, historical, and scenario, worldwide. We present common approaches to assessing facility impact estimates by combining measured or estimated ground-motion parameters with vulnerability functions. Vulnerability functions can be pre-computed, use one of these ground-motion parameters as input, and produce a multi-state discrete output of damage likelihood (or, inspection priorities). Two custom modules for evaluating Arias Intensity and Newmark landslide displacements are presented as examples of using the ShakeCast application with more advanced, customized damage functions.

A pair of significant earthquakes occurred on conjugate faults in the western Imperial Valley involving the through-going Superstition Hills fault and the Elmore Ranch cross fault. The first event was located on the Elmore Ranch fault, Ms = 6.2, and the larger event on the Superstition Hills fault, Ms = 6.6. The latter event is seen as a doublet teleseismically with the amplitudes in the ratio of 1:2 and delayed by about 8 sec. This 8-sec delay is also seen in about a dozen strong-motion records. These strong-motion records are used in a constrained least-squares inversion scheme to determine the distribution of slip on a 2-D fault. Upon closer examination, the first of the doublets was found to be itself complex requiring two episodes of slip. Thus, the rupture model was allowed to have three separate subevents, treated as separate ruptures, with independent locations and start times. The best fits were obtained when all three events initiated at the northwestern end of the fault near the intersection of the cross-fault. Their respective delays are 2.1 and 8.6 sec relative to the first subevent, and their moments are 0.4, 0.9, and 3.5 × 1025 dyne-cm, which is about half of that seen teleseismically. This slip distribution suggests multi-rupturing of a single asperity with stress drops of 60, 200, and 15 bars, respectively. The first two subevents were confined to a small area around the epicenter while the third propagated 18 km southwestward, compatible with the teleseismic and afterslip observations.

The 26 December 2004 Sumatra-Andaman earthquake initiated slowly, with small slip and a slow rupture speed for the first 40 to 60 seconds. Then the rupture expanded at a speed of about 2.5 kilometers per second toward the north northwest, extending 1200 to 1300 kilometers along the Andaman trough. Peak displacements reached ∼15 meters along a 600-kilometer segment of the plate boundary offshore of northwestern Sumatra and the southern Nicobar islands. Slip was less in the northern 400 to 500 kilometers of the aftershock zone, and at least some slip in that region may have occurred on a time scale beyond the seismic band.

We investigate the rupture process of the 1999 Chi-Chi, Taiwan, earthquake using extensive near-source observations, including three-component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three-plane fault geometry derived from our previous inversion using only static data [ Ji et al., 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al., 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10^20 N m and that most of the slip occurred in a triangular-shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s^(−1), and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi-Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

We present a new procedure for the determination of rupture complexity from a joint inversion of static and seismic data. Our fault parameterization involves multiple fault segments, variable local slip, rake angle, rise time, and rupture velocity. To separate the spatial and temporal slip history, we introduce a wavelet transform that proves effective at studying the time and frequency characteristics of the seismic waveforms. Both data and synthetic seismograms are transformed into wavelets, which are then separated into several groups based on their frequency content. For each group, we use error functions to compare the wavelet amplitude variation with time between data and synthetic seismograms. The function can be an L1 + L2 norm or a correlative function based on the amplitude and scale of wavelet functions. The objective function is defined as the weighted sum of these functions. Subsequently, we developed a finite-fault inversion routine in the wavelet domain. A simulated annealing algorithm is used to determine the finite-fault model that minimizes the objective function described in terms of wavelet coefficients. With this approach, we can simultaneously invert for the slip amplitude, slip direction, rise time, and rupture velocity efficiently. Extensive experiments conducted on synthetic data are used to assess the ability to recover rupture slip details. We, also explore slip-model stability for different choices of layered Earth models assuming the geometry encountered in the 1999 Hector Mine, California, earthquake.

We present a rupture model of the Hector Mine earthquake (M 7.1), determined from the joint inversion of strong-motion records, P and SH teleseismic body waves, Global Positioning System (GPS) displacement vectors, and measured surface offset. We solve for variable local slip, rake angle, rise time, and rupture velocity of a finite-fault model involving multiple segments. The inversion methodology developed in a companion article (Ji et al., 2002) combines a wavelet transform approach with a nonlinear (simulated annealing) algorithm. The final model is checked by forward simulating the Interferometric Synthetic Aperture Radar (InSar) data. Our estimation to the seismic moment is 6.28 × 10^(19) N m, which is distributed along three segments from north to south, releasing 37%, 41%, and 22% of the total moment, respectively. The average slip is 1.5 m, with peak amplitudes as high as 7 m. The fault rupture has an average rise time of 3.5 sec and a relatively slow average rupture velocity (1.9 km/sec) resulting in a 14-sec rupture propagation history. Our approach permits large variation in rupture velocity and rise time, and indicates that rise time appears to be roughly proportional to slip and shorter rise times are associated with the initiation of asperity rupture. We also find evidence for nearly simultaneous rupture of the two northern branches.

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.

We introduce a technique for using broadband seismograms recorded from earthquakes at local and regional distances to refine basin structure. For the region outside the basin, we assume a one-dimensional (1D) crustal model and analytical techniques (GRT) to propagate the energy from sources to the basin edge where the motions are then interfaced with a (2D) finite-difference algorithm (Wen and Helmberger, 1996). We parameterize the basin section by isovelocity layers with linear dipping segments between control points. The control point depths are allowed to vary to improve the modeling of waveform data of stations inside the basin. The comparison between data and synthetics is qualified by a fitness function defined by two factors; the timing shift required for best alignment and the correlation coefficient. The procedure was applied to a strong-motion waveform profile across the extended Los Angeles Basin produced by the 1992 Landers, California earthquake to refine the velocity structure using sensitivity testing and forward modeling. Only the correlation coefficient and amplitude were used because absolute timing was unknown. The procedure was extended to a direct waveform inversion by employing a conjugate gradient approach, which uses numerical derivatives. Numerical tests using the new inversion process with synthetic data demonstrate that it is possible to recover a detailed basin structure, if a sufficient amount of high-quality data exists.

The San Andreas fault is the longest fault in California and one of the longest strike-slip faults anywhere in the world, yet we know little about many aspects of its behavior before, during, and after large earthquakes. We conducted a study to locate and to estimate magnitudes for the largest foreshocks and aftershocks of the 1857 M 7.9 Fort Tejon earthquake on the central and southern segments of the fault. We began by searching archived first-hand accounts from 1857 through 1862, by grouping felt reports temporally, and by assigning modified Mercalli intensities to each site. We then used a modified form of the grid-search algorithm of Bakun and Wentworth, derived from empirical analysis of modern earthquakes, to find the location and magnitude most consistent with the assigned intensities for each of the largest events. The result confirms a conclusion of Sieh that at least two foreshocks (“dawn” and “sunrise”) located on or near the Parkfield segment of the San Andreas fault preceded the mainshock. We estimate their magnitudes to be M ≈ 6.1 and M ≈ 5.6, respectively. The aftershock rate was below average but within one standard deviation of the number of aftershocks expected based on statistics of modern southern California mainshock-aftershock sequences. The aftershocks included two significant events during the first eight days of the sequence, with magnitudes M ≈ 6.25 and M ≈ 6.7, near the southern half of the rupture; later aftershocks included a M ≈ 6 event near San Bernardino in December 1858 and a M ≈ 6.3 event near the Parkfield segment in April 1860. From earthquake logs at Fort Tejon, we conclude that the aftershock sequence lasted a minimum of 3.75 years.

Occurrence of large earthquakes close to cities in California is inevitable. The resulting ground shaking will subject buildings in the near-source region to large, rapid displacement pulses which are not represented in design codes. The simulated Mw7.0 earthquake on a blind-thrust fault used in this study produces peak ground displacement and velocity of 200 cm and 180 cm/sec, respectively. Over an area of several hundred square kilometers in the near-source region, flexible frame and base-isolated buildings would experience severe nonlinear behavior including the possibility of collapse at some locations. The susceptibility of welded connections to fracture significantly increases the collapse potential of steel-frame buildings under strong ground motions of the type resulting from the Mw7.0 simulation. Because collapse of a building depends on many factors which are poorly understood, the results presented here regarding collapse should be interpreted carefully.

High-rise flexible-frame buildings are commonly considered to be resistant to shaking from the largest earthquakes. In addition, base isolation has become increasingly popular for critical buildings that should still function after an earthquake. How will these two types of buildings perform if a large earthquake occurs beneath a metropolitan area? To answer this question, we simulated the near-source ground motions of a M(w) 7.0 thrust earthquake and then mathematically modeled the response of a 20-story steel-frame building and a 3-story base-isolated building. The synthesized ground motions were characterized by large displacement pulses (up to 2 meters) and large ground velocities. These ground motions caused large deformation and possible collapse of the frame building, and they required exceptional measures in the design of the base-isolated building if it was to remain functional.

A preliminary rupture model of the 1994 Northridge, California earthquake, determined from strong motion waveform inversion and analysis, is presented. The fault rupture plane chosen is based on the distributions of aftershocks and teleseismic surface-wave and bodywave point-source solutions. The fault strikes 122°, dips 42°, and has a slip vector of 109°. The average slip is about 1.2 meters over the rupture area and the peak slip reaches nearly 4 meters. Our estimate of the seismic moment is 1.2 ± 0.2 x 10^(26) dyne-cm. The area of rupture is small relative to the aftershock dimensions and is approximately 14 km along strike (west-northwest from the hypocenter) and nearly 20 km in the updip direction. There is little indication of slip shallower than about 7 km. The up-dip, near-source strong-motion velocity waveforms show two distinct, large S-wave arrivals 2-3 sec apart (as do the teleseismic P waves), indicating separate source subevents. An along strike (west-northwest) subevent separation of about 8 km is most consistent with the observation that the two main arrivals are separated more in time to the south and southeast (about 4.5 sec at Stone Canyon Reservoir and Sherman Oaks, for example), than at northern azimuths. The interpretation of secondary arrivals observed at more distant stations to the south and southeast (e.g., Santa Monica) is more tenuous, since several of the aftershocks recorded there indicate later arrivals as well. However, a secondary source contribution is expected based on our model of the closer stations. After placing these constraints on the general nature of the rupture, we predict the characteristics of the long-period (1-10 sec) ground velocities over a grid of stations covering the near-source region. This exercise provides a basis for separating the effects of source radiation (dominated by radiation pattern and directivity) from the complex waveform modifications due to wave-propagation and site response.

The Advanced National Seismic System (ANSS),an initiative begun in 1998 to integrate, expand, and modernize seismic monitoring nationwide, is providing improvements in earthquake monitoring and in the development, production, and delivery of earthquake data, information, and science. The ANSS initiative is an element of the Earthquake Hazard Program (EHP) of the U.S.Geological Survey (USGS). The USGS National Earthquake Information Center (NEIC),an integral part of ANSS, has incorporated into its routine operations many of the developments in data processing and data products supported by the ANSS. The four principal elements of ANSS monitoring operations are, in sequence: data collection, data processing, generation of data products, and product dissemination. This article describes the current status and new developments in each of these elements, and how they have affected NEIC capabilities.

Short period and broadband teleseismic waveform data and three-component strong motion records were analyzed to obtain the source parameters of the 1991 Sierra Madre earthquake. Close-in strong motion velocity records (analyzed from 5 s to 5 Hz) show two distinct pulses about 0.35 s apart, requiring some rupture complexity. The near-field shear wave displacement pulse from this event has a relatively short duration (about 1 s) for the magnitude of the event, requiring a particularly high average stress drop. To further constrain the rupture process, the data were used in a finite fault source inversion to determine the temporal and spatial distribution of slip. We chose a fault plane orientation striking S 62°W and dipping 50° toward the northwest as required by the distribution of aftershocks, the first motion mechanism and the teleseismic bodywave point source inversion. In addition to the aftershock locations, depth constraints are provided by teleseismic short period and broadband recordings which require a centroid depth of 10–11 km. Our inverse modeling results indicate that both the teleseismic and strong motion data sets can be fit with a compact rupture area, about 12 km2, southwest and up-dip from the hypocenter. The average slip is approximately 50 cm, and the maximum slip is 120 cm. The seismic moment obtained from either of the separate data sets or both sets combined is about 2.8 ± 0.3 × 1024 dyn cm and the potency is 0.01 km3. Using the area of significant slip estimated from the finite fault inversion, the resulting stress drop is on the order of 150–200 bars.

The section of the Southwest Indian Ocean Ridge west of the Prince Edward Fracture zone has a large ridge axis offset and a complicated ridge-transform morphology. We have determined the source mechanisms of transform earthquakes along this portion of the ridge from an inversion of long-period P and SH waveforms. The seismicity is characterized by anomalous faulting mechanisms, source complexity and an unexpectedly large seismic moment release. Several earthquakes with dip-slip components of faulting have been recognized on the central section of the Andrew Bain and 32° E transforms suggesting geometrical complexity along the transform. This region has experienced a Mw = 8.0 transform earthquake in 1942, yet we observe a seismic slip rate during the last 20 years that is still comparable to the predicted spreading rate (1.6 cm/yr). The calculated slip rate over a period of 60 years is three times greater than the expected rate of spreading.

Although the fact that sediments can amplify earthquake ground motion was recognized at least 100 years ago (Milne, 1898), there has been a lingering uncertainty as to whether the degree of amplification varies with the level of input motion. This issue remains as one of the most important questions with respect to understanding and predicting earthquake ground motion. In accordance with the conservation of energy, seismicwave amplitudes generally increase in sediments due to lower densities and and/or lower seismic velocities. In addition, resonance effects can occur where abrupt impedance contrasts exist. If sediments were perfectly elastic, their response would be independent of incident-wave amplitudes. As with any real material, however, sediments begin to yield at some level of strain, and this violation of Hooke's law will give rise to a nonlinear response. The engineering community has long believed that sediment nonlinearity is significant. This perspective was based almost entirely on...