Search

The Tonga subduction zone in the south-west Pacific is the fastest convergent plate boundary in the world with the most active mantle seismicity. This zone shows unique tectonic features including Samoan volcanic lineament of plume-driven origin near the northern rim of the Tonga subducting slab. The proximity of the Samoa hotspot to the slab is enigmatic and invokes debates on interactions between the Samoa plume and the Tonga subduction. Based on long-term observations of intermediate and deep-focus Tonga earthquakes reported in the Global Centroid Moment Tensor (CMT) catalog, we provide novel detailed imaging of this region. Accurate traveltime residua of the P- and S-waves recorded at two nearby seismic stations of the Global Seismographic Network are inverted for the P- and S-wave velocities and their ratio and reveal their pronounced lateral variations. In particular, they differ for the southern and northern parts of the Tonga subduction region. While no distinct anomalies are detected in the southern Tonga segment, striking low-velocity anomalies associated with a high Vp/Vs ratio are observed in the northern Tonga segment close to the Samoa plume. These anomalies spread through the whole upper mantle down to depths of ~ 600 km. Together with the fast extension of the northern back-arc Lau Basin, slab deformation and geochemical enrichment in the northern Tonga region, they trace deep-seated magmatic processes and evidence an interaction of the Tonga subduction with the Samoa plume.

The relative contributions of scattering and viscoelastic attenuation to the apparent attenuation of seismic body waves are estimated from synthetic and observed S waves multiply reflected from Earth’s surface and the core-mantle boundary. The synthetic seismograms include the effects of viscoelasticity and scattering from small-scale heterogeneity predicted from both global tomography and from thermodynamic models of mantle heterogeneity that have been verified from amplitude coherence measurements of body waves observed at dense arrays. Assuming thermodynamic models provide an estimate of the maximum plausible power of heterogeneity measured by elastic velocity and density fluctuations, we predict a maximum scattering contribution of 43 % to the total measured attenuation of mantle S waves having a dominant frequency of 0.05 Hz. The contributions of scattering in the upper and lower mantle to the total apparent attenuation are estimated to be roughly equal. The relative strength of the coda surrounding observed ScSn waves from deep focus earthquakes is not consistent with a mantle having zero intrinsic attenuation.

Mantle seismic anisotropy beneath the Amur Plate is studied using the data of ScS waves reflected from the outer core of the Earth from local deep earthquakes in the area of five stations located in Primorye and Priamurye regions. The results of measuring the parameters of the split ScS waves near the stations show the dominance of the polarization azimuths of the fast wave along the eastern directions and are consistent with azimuthal anisotropy and the direction of motion of the Pacific Plate (300°) and Amur Plate (~120°) depending on the epicenter–station direction. In Primorye, in the vicinity of the TEY station, the polarization azimuths of the fast ScS wave dominate in the interval of NE–E directions orthogonally to the lines of the mantle flow along the complex 3D surface of the subsiding Pacific Plate. It is revealed that the delay time of ScS waves increases to 2 s and 3.4 s in the upper mantle and in the mantle transition zone, respectively, as the depth of events increases. The highest degree of anisotropy manifests itself in the upper mantle. In the case of the vertical propagation of waves under conditions of horizontal mantle flow, the difference in the arrival times of the ScS waves is the lowest. The anisotropy in the upper zone of the transition mantle can be related to the wadsleyite texture with the polarization of the fast ScS wave parallel or orthogonally to the motion of the stagnant plate. The anisotropy in the upper part of the lower mantle is associated with the texture of perovskite and periclase with the orientation of the symmetry axis and the polarization of the ScS waves in parallel to the plate subsidence and the subduction direction.

We used recent seismic data and advanced techniques to investigate 3D fault geometry over the transition from the partially coupled to the fully coupled plate interface inboard of the Nankai Trough off the Kii Peninsula, Japan. We found that a gently dipping plate boundary decollement with a thick underthrust layer extends beneath the entire Kumano forearc basin. The 1 April 2016 Off-Mie earthquake (Mw6.0) and its aftershocks occurred, where the plate boundary decollement steps down close to the oceanic crust surface. This location also lies beneath the trenchward edge of an older accretionary prism (∼14 Ma) developed along the coast of the Kii peninsula. The strike of the 2016 rupture plane was similar to that of a formerly active splay fault system in the accretionary prism. Thus, the fault planes of the 2016 earthquake and its aftershocks were influenced by the geometry of the plate interface as well as splay faulting. The 2016 earthquake occurred within the rupture area of large interplate earthquakes such as the 1944 Tonankai earthquake (Mw8.1), although the 2016 rupture area was much smaller than that of the 1944 event. Whereas the hypocenter of the 2016 earthquake was around the underplating sequence beneath the younger accretionary prism (∼6 Ma), the 1944 great earthquake hypocenter was close to oceanic crust surface beneath the older accretionary prism. The variation of fault geometry and lithology may influence the degree of coupling along the plate interface, and such coupling variation could hinder slip propagation toward the deeper plate interface in the 2016 event.

This paper presents the results of reconstruction of the modern tectonic stress field of the northern part of the Kuril–Okhotsk region prior to the May 24, 2013 deep-focus earthquake. This earthquake is the strongest deep-focus earthquake not only in the Okhotsk region but also in the world over the entire period of seismic observations. The tectonic stress field is reconstructed using the method of cataclastic analysis (MCA) of the earthquake source mechanism data. New data on the depth variations of regional tectonic stress field are obtained.

Robust relations correlating 10 different magnitudes of intermediate‐depth and deep‐focus earthquakes to moment magnitude are proposed, in order to be efficiently incorporated into the compilation process of homogeneous (with respect to magnitude) earthquake catalogs. By using global data available from International Seismological Centre (ISC), National Earthquake Information Center (NEIC), Comprehensive Nuclear‐Test‐Ban Treaty Organization’s International Data Centre (IDC), Institute of Physics of the Earth in Moscow, Russia, and China Earthquake Networks Center in Beijing, the performance of several widely used magnitude scales, such as body wave ( m b , m B ) and surface wave ( M s ), is examined with respect to the moment magnitude scale ( M w ). Similarly, appropriate M w ‐calibrated relations are also provided for regional magnitude scales such as the M JMA magnitude calculated by the Japan Meteorological Agency. The analysis also involves the integration of focal depth as an additional variable to some of the above magnitude‐conversion relations. This depth effect proved to be important for ISC/NEIC’s body‐wave ( m b IN) and surface‐wave ( M s IN) magnitudes, leading to significant corrections for the estimated magnitudes. More specifically, a major change in the magnitude residual variation around the depth of 230 km was identified for the case of m b IN. Furthermore, the obtained results provide important observations on the behavior of certain magnitude scales. A typical case is the m b scale reported by IDC, which shows a systematic and large bias, with respect to the published M w data values.

The 2014 M w 7.9 Rat Islands intermediate-depth (∼107 km) earthquake is studied using seismic observations at teleseismic distances. The nodal plane with a steep-dip angle (87°) and a strike of 308° inferred from the Global CMT solution is selected as the causative fault plane because of the consistency with the distribution of the aftershocks relocated using the Joint Hypocenter Determination method. Its slip history is constrained using broadband P and SH waveforms and long-period (3–6 mHz) seismic waves in vertical and transverse components. The rupture process associated with this plane shows a rupture propagating at ∼2 km/s. The rupture essentially breaks one large slip patch (peak slip amplitude of 3.5 m) located 30 km west of the hypocenter. The size of this slip patch extends from 80 to 140 km depth, which would be the largest depth extent of earthquake slip ever documented. We show that this earthquake breaks the entire cold core of the subducting slab defined by the depth of the 600–650 °C isotherm. In addition, the lateral extent of the rupture appears to be bounded by the heterogeneities of subducting slab. Our results further indicate that, along this segment of Aleutian Arc, the intermediate-depth earthquakes might contribute to the westward transporting motion of Rat block. Lastly, the slip distribution of this earthquake as well as that of two other well studied M w > 7.5 intermediate-depth earthquakes suggests that the coldest portion of the subducting slab, outlined by the double Wadati–Benioff zone, is seismic during large intermediate-depth earthquakes. Such a result is important for estimating the seismic hazard of future large intermediate-depth earthquakes.

The existence of earthquakes within continental lithospheric mantle remains a highly controversial topic. Here, we present a detailed set of seismological analyses confirming the occurrence of a mantle earthquake beneath the Wind River Range of central Wyoming. Combining regional waveform inversion with the analysis of the delay and relative amplitudes of teleseismically-observed depth phases, we demonstrate that the 2013 Wind River earthquake – a M W 4.7 highly-oblique thrust-faulting event – occurred at 75 ± 8 km , well beneath the base of the crust. The magnitude, mechanism, and location of this earthquake suggest that it represents simple brittle failure at relatively high temperatures within the mantle lithosphere, as a result of tectonic, rather than magmatic, processes.

We apply a multiple source inversion method to systematically study the source processes of 25 large deep-focus (depth >400 km) earthquakes with Mw > 7.0 from 1994 to 2012, based on waveform modeling of P, pP, SH and sSH wave data. The earthquakes are classified into three categories based on spatial distributions and focal mechanisms of the inferred sub-events: 1) category one, with non-planar distribution and variable focal mechanisms of sub-events, represented by the 1994 Mw 8.2 Bolivia earthquake and the 2013 Mw 8.3 Okhotsk earthquake; 2) category two, with planar distribution but focal mechanisms inconsistent with the plane, including eighteen earthquakes; and 3) category three, with planar distribution and focal mechanisms consistent with the plane, including six earthquakes. We discuss possible physical mechanisms for earthquakes in each category in the context of plane rupture, transformational faulting and shear thermal instability. We suggest that the inferred source processes of large deep-focus earthquakes can be best interpreted by cascading failure of shear thermal instabilities in pre-existing weak zones, with the perturbation of stress generated by a shear instability triggering another and focal mechanisms of the sub-events controlled by orientations of the pre-existing weak zones. The proposed mechanism can also explain the observed great variability of focal mechanisms, the presence of large values of CLVD (Compensated Linear Vector Dipole) and the super-shear rupture of deep-focus earthquakes in the previous studies. In addition, our studies suggest existence of relationships of seismic moment ∼ (source duration) 3 and moment ∼ (source dimension) 3 in large deep-focus earthquakes.

SUMMARY Deep focus earthquakes within the underthrust Indian lower crust beneath the Himalaya occur in very specific regions and have distinct source characteristics. The study of the source mechanisms of these earthquakes provides valuable constraints on the kinematics of deformation of the underthrust Indian Plate, and its influence on the active deformation of the overlying Himalayan wedge. One of the most significant regions of these deep focus earthquakes is beneath the Sikkim and Bhutan Himalaya. We study the source characteristics of the 2011 September 18 (Mw 6.9) deep focus Sikkim main shock and its major aftershocks using global, regional and local waveform data. We determined the focal mechanism of the main shock using moment tensor inversion of global P and SH waveforms, and ascertained the earthquake fault plane using rupture directivity from regional P-wave spectra. The main shock originated at 53 ± 4 km depth and ruptured at least 20 km thickness of the underthrust Indian lower crust. Faulting occurred on a near vertical dextral strike-slip fault oriented NW-SE (strike 127 ◦ ,d ip 81 ◦ and rake 167 ◦ ), oblique to the local strike of the Himalayan arc. The rupture initiated from the SE end of the fault and propagated to the northwest. The main shock was followed by 20 small-to-moderate aftershocks (mb > 3.0), which we relocated using phase arrival times. We computed the focal mechanisms of the larger ones (mb ≥ 3.5) using local waveform inversion. We find that all aftershocks originated SE of the main shock, between depths of 12 and 50 km, and have dominantly strike-slip mechanisms. Our results, combined with the source mechanisms of earthquakes from previous studies, reveals that the entire underthrust Indian crust is seismogenic and deforms by dextral strikeslip motion on oblique structures beneath the Sikkim and Bhutan Himalaya. These active oblique structures with transverse motion possibly mark the western boundary of the clockwise rotating ‘microplates’ in northeast India observed from GPS geodesy.

We described the 2014 January 17 earthquake (M w = 4.3) occurred on the Siberian platform in the area of sublongitudinal part of the Angara river in the zone of possible influence of two large reservoirs—Ust-Ilimsk and Boguchan. This is the first event of such magnitude recorded in this previously aseismic area during the whole period of instrumental observations. A seismic moment, a moment magnitude, a hypocentral depth, and a focal mechanism of the event were calculated on the basis of surface wave amplitude spectra. Analysis of the geological and geophysical data showed that the earthquake origin is connected with high-velocity gradient zone located at the border of the Late Proterozoic cover and the Precambrian basement of the Siberian platform. Some evidences for a natural character of the earthquake were considered.

It has been hypothesized that much of the crustal deformation attributed to the Laramide orogeny of the southwest North American Cordillera was caused by dynamic effects induced by the flat subduction of a large oceanic plateau that was embedded within the Farallon plate. Previous studies have identified within the North American mantle a seismic velocity anomaly that plausibly represents the remnants of the subducted plateau. Coupled plate kinematic and dynamic modeling of the anomaly, as well as surface geological findings identify this anomaly as the subducted conjugate to the Shatsky Rise. Here, we find clear evidence for a northeastward dipping (35° dip), thick (up to 400 km thick) slab-like seismic velocity anomaly within the top of the lower mantle below the central United States. Using a deep focus earthquake below Spain, we find that the observed seismic waveforms recorded with the dense USArray display multi-pathing indicative of sharp surfaces. Plate tectonic reconstructions in which the anomaly is migrated back in time suggest strong coupling of the plateau-thickened slab segment to the overriding North America Plate as it was subducted. In combination with the reconstructions, we interpret the structure as arising from eastward dipping Farallon subduction at the western margin of North America during the Cretaceous, in contrast with some recent interpretations. It appears that the plateau area of the slab has been further thickened, which might undergo a combination of pure shear bulk shortening during flat-slab subduction and/or by a folding instability during penetration into the lower mantle.

The parameters of split shear ScS and S waves from the strongest (Mw = 8.4) Okhotsk earthquake and its strong aftershocks that occurred in the mantle’s transitional zone in the backarc part of the Kamchatka subduction zone are measured. The measurement results show domination of the fast ScS wave’s east azimuths (110°–149°), which are orthogonal to the trench orientation, and the time difference between the arrivals of the ScS waves (δt) in the range of 0.9–1.6 s. The fast S wave from the Okhotsk earthquake is east-directed (89°), and the time difference between its arrivals is up to 2.5 s. As for aftershocks, the azimuths of the fast S waves are oriented in parallel to the trench, and δt are estimated at 1.2–1.3 s. The orientations of the azimuths for the ScS and S waves are consistent with the model of the transverse-isotropic symmetry of the medium with the symmetry axis tilted along the plate dip and/or along the trench strike.

From 2011 to 2012, we conducted marine seismic surveys using a multichannel seismic reflection (MCS) system and ocean bottom seismometers (OBSs) in and around the hypocentral region of the 1983 Nihonkai–Chubu earthquake. The survey covered the area from the Japan coast of the Japan Sea to the Yamato Basin and the Japan Basin. We investigated the relationship between the Nihonkai–Chubu earthquake and the crustal structure obtained from seismic reflection imaging and the P-wave velocity structure near the hypocentral region. The focus in the present paper is on seismic data that were acquired using both the MCS system and OBSs. Based on the results, the crustal structure of the Japan Sea around the hypocentral region of the 1983 Nihonkai–Chubu earthquake is divided into three types: island arc crust, thick oceanic crust, and oceanic crust. The Nihonkai–Chubu earthquake occurred near the boundary between the thick oceanic crust and the island arc crust. There is a strong possibility that the east-dipping reflector detected in the crust was the source fault of the Nihonkai–Chubu earthquake. In addition, the location of anticlines formed by the east-dipping reverse fault occurring in the boundary of the crustal structure associated with the source fault of the Nihonkai–Chubu earthquake was confirmed in the range of latitude from 39 ° 20 ′ to 41 ° 15 ′ N , although they were not found south of this segment. There is also a concentrated area of reverse faults and folds approximately 100 km west of the source region of the Nihonkai–Chubu earthquake. This area is located near the boundary between the thick oceanic crust and the oceanic crust.

On 24 May 2013, the largest ever-recorded deep earthquake occurred beneath Sea of Okhotsk. A multiple point source inversion procedure is applied to constrain source process of this earthquake, based on waveform modeling of both direct P and SH waves and near-surface reflected pP and sSH waves. Our results indicate that the earthquake consists of six major subevents separated in space and time, encompassing a horizontal dimension of 64 ± 4 km along ~ N160°E and a downward depth extension of 35 ± 4 km. The geographic distribution and focal mechanisms of the inferred subevents and foreshock/aftershock locations do not fit into plane rupture. We suggest that the earthquake can be best explained by a cascading failure of shear instability within preexisting weak zones in the region, with the perturbation of stress generated by a shear instability triggering another.

The physical mechanism of deep earthquakes (depth >300 km) remains enigmatic, partly because their rupture dimensions are difficult to estimate due to their low aftershock productivity and absence of geodetic or surface rupture observations. The two largest deep earthquakes, the recent Great 2013 Sea of Okhotsk earthquake (M 8.3, depth 607 km) and the Great 1994 Bolivia earthquake (M 8.3, depth 637 km), together provide a unique opportunity to compare their rupture patterns in detail. Here we extend a travel-time sub-event location method to perform full teleseismic P-waveform inversion. This new method allows us to explain the observed broadband records with a set of sub-events whose model parameters are robustly constrained without smoothing. We find that while the Okhotsk event is mostly unilateral, rupturing 90 km along strike with a velocity over 4 km/s, the Bolivia earthquake ruptured about half this distance at a slow velocity (about 1.5 km/s) and displayed a major change in rupture direction. We explain the observed differences between the two earthquakes as resulting from two fundamentally different faulting mechanisms in slabs with different thermal states. Phase transformational faulting is inferred to occur inside the metastable olivine wedge within cold slab cores whereas shear melting occurs inside warm slabs once triggered.

A sequence of 98 teleseismically recorded earthquakes occurred off the east coast of Kamchatka at depths between 10-90 km around latitude 52.5°N and longitude 160°E on May 16–23, 2013. The swarm occurred along the northern limit of the rupture area of the 1952 Mw 9.0 great Kamchatka earthquake, the fifth largest earthquake in the history of seismic observations. On May 24, 2013 the strongest deep earthquake ever recorded of Mw 8.3 occurred beneath the Sea of Okhotsk at a depth of 610 km in the Pacific slab of the Kamchatka subduction zone, becoming the northernmost deep earthquake in the region. The deep Mw 8.3 earthquake occurred down-dip of the shallow swarm in a transition zone between the southern deep and northern shallow segments of the Pacific slab. Several deep aftershocks followed, covering a large, laterally elongated part of the slab. We suppose that the two described earthquake sequences, the May 16–23 shallow earthquake swarm and the May 24–28 deep mainshock-aftershock series, represent a single tectonic event in the Pacific slab having distinct properties at different depth levels. A low-angle underthrusting of the shallow part of the slab recorded by the shallow earthquake swarm activated the deep part; this process induced the deep mainshock-aftershock series only three days after the swarm. The domain of the subducting slab activated by the May 2013 earthquake occurrence was extraordinarily large both down-dip and along-strike.

S U M M A R Y We explore the use of on-land GPS observations to detect deformation due to tsunami propagation near source regions of large interplate earthquakes. Here, we focus on the Mw 9 Tohoku-oki earthquake, which occurred around 14:46 (JST) on 2011 March 11. We consider GPS data in the time span 14:54–15:22 (JST) along the Sanriku coast, where the tsunami had the largest amplitude. The displacement data shows the signatures of large aftershocks as well as post-seismic fault slip (afterslip). These effects are particularly evident in the east component. From the horizontal displacement vectors, we construct a simple fault model for the early phase of the afterslip. Mean slip velocity of the early afterslip reaches 0.1 mm s−1. By compiling the early afterslip velocity of recent interplate earthquakes around that region, we find its increasing trend with the main shock magnitude. This scaling relation may reflect higher stressing rates at edges of larger main shock faults. Separately, we forward calculate land deformation due to tsunami height changes based on a tsunami simulation. Tsunami-induced deformation is only evident in the vertical direction at coastal GPS stations. The predicted subsidence amounts at some coastal stations can account for a large portion of the residuals between the observation and the modelled deformation due to the fault slip.

This paper discusses results from ongoing research on long-term earthquake prediction for the Kuril–Kamchatka island arc based on the concepts of seismic gaps and the seismic cycle. We developed a forecast for the next 5 years (April 2016 through March 2021) for all segments of the earthquake-generating zone along the Kuril–Kamchatka arc. The 20 segments of the arc were analyzed to develop forecasts of the appropriate phases of the seismic cycle, a normalized parameter of the rate of small earthquakes (A10), the magnitudes of moderate earthquakes that are expected with probabilities of 0.8, 0.5, and 0.15, the maximum expected magnitudes, and the probabilities of great (M ≥ 7.7) earthquakes. We discuss the seismic process in the Kuril–Kamchatka earthquake-generating zone before and after the deep-focus May 24, 2013 M 8.3 earthquake in the Sea of Okhotsk. The results corroborate the high seismic hazard in the area of Petropavlovsk-Kamchatskii and the urgent need for continuing with and expanding the ongoing work of seismic retrofitting and seismic safety enhancement. We quote practical results from applications of the method during 30 years.

We present results for joint kinematic inversion of high-rate GPS, strong motion and InSAR data for the 2016 M7.6 Melinka earthquake. We show that the source is a compact 35 s long rupture with 5.0 ± 0.15 m of peak slip. We find the Melinka earthquake occurs inside the slip region for the 1960 M9.5 Valdivia earthquake and within the highly locked portion of the megathrust inferred from inter-seismic velocity analysis. We show that there is very modest post-seismic deformation at a nearby GPS site QLLN and argue that this indicates the Melinka earthquake ruptures within the intermediate portion of the megathrust and is an isolated asperity surrounded by locked velocity weakening material. Further we find that the peak slip observed during this earthquake is larger than what has been accumulated in the intervening 57 years since the 1960 rupture and conclude that, at least in the area of the Melinka earthquake, this is indicates that the 1960 Valdivia event might not have used all the slip deficit available on the megathrust.

[1] Mass rearrangement within the crust and lithospheric mantle and ocean water redistribution caused by great earthquakes are made visible by their coseismic gravity signature, nowadays detectable by the Gravity Recovery And Climate Experiment (GRACE) space mission. Here we present a novel procedure for estimating the principal seismic source parameters (hypocenter and moment tensor) that relies solely on space gravity data from GRACE, applied for the first time to the 2011 Tohoku earthquake. It yields a seismic source model that is consistent with a thrust earthquake and geological information of the subduction zone. It closely resembles the Global Centroid-Moment-Tensor Project solution based on the teleseismic wave inversion, although the moment magnitude is higher (9.13 ± 0.11 compared to 9.08) and the hypocenter is further offshore by about 40 km, within the oceanic plate. This procedure will become an important tool in seismology as it complements Centroid-Moment-Tensor analysis by exploiting the new gravity data from GRACE.

The 2016 Ecuador M _W 7.8 earthquake ruptured the subduction zone boundary between the Nazca and South American plates. Joint modeling of seismic and tsunami observations indicates an ∼120 km long rupture area beneath the coastline north of the 1998 M_W 7.2 rupture. The slip distribution reveals two discrete asperities near the hypocenter and around the equator. Their locations and the patchy pattern are consistent with the prior interseismic geodetic strain, which showed highly locked patches also beneath the coastline. Aftershocks cluster along two streaks, one aligned nearly parallel to the plate convergence direction up-dip of the main slip patches, and the other on a trench-perpendicular lineation south of the 1958 rupture zone. Comparisons of seismic waveforms and magnitudes show that the 2016 event and 1942 earthquakes have similar surface wave magnitude M_S 7.5), overlapping rupture areas, and similar main pulses of moment rate. The same area ruptured as the southernmost portion of the larger earthquake of 1906 (M_W 8.6, M_S 8.6). The seismic behavior reflects persistent heterogeneous frictional properties of the Colombia–Ecuador megathrust.

Using about 5200 handpicked P, S, pP, sP, PcP, and ScP phases from digitally recorded seismograms, together with International Seismological Centre reported phases, we obtain improved hypocentral locations for ∼1790 earthquakes, ∼983 of them having 90% confidence limits 50 km, in the period 1962 to September 1996, along the Indonesian subduction zone. We use an improved joint hypocenter determination method in which the solution is more stable than in the original one. We show that for deep earthquakes, use of the core-reflected phases is as good as use of depth phases, and since such core phases often have larger amplitudes and sharper onsets than the depth phases, their arrival times can be read very accurately. The positions of the relocated hypocenters show that (1) a portion of the Indonesian arc between ∼ 110° E and 123° E longitude, and deeper than ∼500 km, is dipping southward at an angle of ∼75°, that is, in a direction opposite to the upper part of the north dipping slab, and (2) east of about 108° E, the seismic zone is wider near 670 km than near the 500 km depth. The first suggests southward lateral flow in the mantle, relative to the plate motion vector here. From the contortion of the seismic zone along the eastern portion of the Indonesian arc, we can estimate the average lateral shear strain rate in the 300 to 670 km depth range to be of the order of 10-16S-1, over the last 10-20 Myr. Copyright 1999 by the American Geophysical Union.

SUMMARY Earthquake source parameter determination is of great importance for hazard assessment, as well as for a variety of scientific studies concerning regional stress and strain release and volcano-tectonic interaction. This is especially true for poorly instrumented, densely populated regions such as encountered in Africa, where even the distribution of seismicity remains poorly documented. In this paper, we combine data from satellite radar interferometry (InSAR) and teleseismic waveforms to determine the source parameters of the Mw 5.9 earthquake that occurred on 2008 February 3 near the cities of Bukavu (DR Congo) and Cyangugu (Rwanda). This was the second largest earthquake ever to be recorded in the Kivu basin, a section of the western branch of the East African Rift (EAR). This earthquake is of particular interest due to its shallow depth and proximity to active volcanoes and Lake Kivu, which contains high concentrations of dissolved carbon dioxide and methane. The shallow depth and possible similarity with dyking events recognized in other parts of EAR suggested the potential association of the earthquake with a magmatic intrusion, emphasizing the necessity of accurate source parameter determination. In general, we find that estimates of fault plane geometry, depth and scalar moment are highly consistent between teleseismic and InSAR studies. Centroid-moment-tensor (CMT) solutions locate the earthquake near the southern part of Lake Kivu, while InSAR studies place it under the lake itself. CMT solutions characterize the event as a nearly pure double-couple, normal faulting earthquake occurring on a fault plane striking 350° and dipping 52° east, with a rake of –101°. This is consistent with locally mapped faults, as well as InSAR data, which place the earthquake on a fault striking 355° and dipping 55° east, with a rake of –98°. The depth of the earthquake was constrained by a joint analysis of teleseismic P and SH waves and the CMT data set, showing that the earthquake occurred in the shallow crust, at approximately 8 km depth. Inversions of ENVISAT (Environment Satellite) and ALOS (Advanced Land Observation Satellite) data place the earthquake at 9 km. A comparison of the scalar moment (9.43 ± 0.06 × 1017 Nm from seismology and 8.99 ± 0.010 × 1017 Nm from the joint InSAR solution) shows good agreement between the two data sets. Such an agreement is in contrast to the large discrepancies observed (up to an order of magnitude) in other places along the EAR where similar earthquake sequences are associated with magmatic intrusion. From this, we infer that the rupture was brittle and occurred with little aseismic deformation as might be expected from magma injection. Our results provide insights into the style of rifting occurring in the South Kivu Volcanic Province and hence will aid future studies on seismic risk in the context of Lake Kivu and underline the importance of systematic monitoring of the EAR area.

The 1963 great Kurile earthquake was an underthrust earthquake occurred in the Kurile–Kamchatka subduction zone. The slip distribution of the 1963 earthquake was estimated using 21 tsunami waveforms recorded at tide gauges along the Pacific and Okhotsk Sea coasts. The extended rupture area was divided into 24 subfaults, and the slip on each subfault was determined by the tsunami waveform inversion. The result shows that the largest slip amount of 2.8 m was found at the shallow part and intermediate depth of the rupture area. Large slip amounts were found at the shallow part of the rupture area. The total seismic moment was estimated to be 3.9 × 1021 Nm (Mw 8.3). The 2006 Kurile earthquake occurred right next to the location of the 1963 earthquake, and no seismic gap exists between the source areas of the 1963 and 2006 earthquakes.

SUMMARY We study a large M w = 7.6 earthquake that occurred on 2007 November 14 in the Northern Chile seismic gap near the city of Tocopilla. Using a variety of seismic data we show that this earthquake ruptured only the lower part of the interplate seismic zone and generated a series of plate interface aftershocks. Two large aftershocks on 2007 November 15 ruptured the interplate zone oceanwards of the Mejillones Peninsula, a major geographical feature in the Antofagasta region. On 2007 December 16, a large M w = 6.8 aftershock, that occurred near the southern bottom of the fault plane of the main event, is shown to be a slab-push earthquake located inside the subducted Nazca Plate and triggered by along slab compression. Aftershocks of this event demonstrate that it occurred on an almost vertical fault. The Tocopilla earthquake took place just after the installation of a new seismological network by Chilean, German and French researchers. The accelerometric data combined with far field seismic data provide a quite complete and consistent view of the rupture process. The earthquake broke a long (130 km) and narrow (about 30–50 km) zone of the plate interface just above the transition zone. Using a non-linear kinematic inversion method, we determined that rupture occurred on two well-defined patches of roughly elliptical shape. We discuss the consequences of this event for models of gap filling earthquakes in Chile proposed in the 1970s.

SUMMARY The 2005 June 13, an intraplate, intermediate depth earthquake (M w = 7.8, 98 km depth) occurred in northern Chile. Previous studies show that this earthquake is a slab-pull event with down dip extensional source mechanism. However, the physical origin and the role in seismic cycle of this kind of event are still under debate. We present a seismological source study from strong motion data, based on the simplified slip patch kinematic parametrization. We find two distinct sources. The main source is almost centred around the hypocentre and triggers the secondary one, situated deeper in the slab. The weaknesses and uncertainties of the kinematic models are partly resolved and discussed by constructing dynamic spontaneous rupture models with two radiative asperities and one embedding zone of non-radiative slip process, a required feature to fulfil the seismic moment constraint. Two kinds of models are compatible with observations: a crack mode model (combined with fast rupture velocity and small zone of non-radiative slip) and a pulse mode model (combined with slower rupture velocity and large zone of non-radiative slip). The fast rupture model is preferred for its physical simplicity. Replacing our result in the geometry and thermal contexts of the subduction, we conclude that the rupture has probably broken the whole lithosphere by encountering little resistance and without being very sensitive to the presence of the double seismic zone. This suggests that the rupture happens along a weak plane, on which large slab instabilities can take place.

SUMMARY investigation of the SH wave seismograms indicates that rupture occurred on a shallowly dipping fault, propagating in all directions from the hypocentre. The average rupture velocity in the downdip and updip directions are ∼2.0-3 and ∼1.0 km s −1 , respectively. The maximum slip is found to be ∼12 m and is located downdip of the earthquake's hypocentre. Taken over the fault as a whole, the average stress drop is ∼0.55 MPa. Used in conjunction with relocated aftershocks, we find that two barriers constrain the faulting extent of slip, creating a rupture area of ∼165 km × 140 km. These barriers coincide with submarine topographic features, are domain boundaries for the seismic gaps along the arc, and have high concentrations of large aftershocks near them. The high slip regions align with the orientation of the subducting plate, and resemble in size and shape the many seamounts about to subduct. The 2003 Tokachi-Oki earthquake was the largest earthquake to occur globally in that year and reruptured the portion of the Kuril- Japan Trench that last broke in the great Tokachi-Oki earthquake of 1952 (Hamada & Suzuki 2004). The abundance of local data in this region coupled with the large magnitude of the event has led to a range of models for the rupture being proposed. These solutions are summarised in Table 1. All of these solutions have the total moment for the earthquake significantly below that obtained from the long-period seismological data by the global CMT (GCMT) solution, with some solutions having less than one third of the GCMT moment. This study aims to address this anomaly, by first investigating the long-period moment tensor solution and then by finding a rupture model which is consistent with both broad-band teleseismic data and the long-period solution.

Relationships between relative plate velocity and seismicity differ by plate-boundary class. We test the null hypothesis of linearity of earthquake rates with velocity in each of seven classes. A linear relationship is expected if earthquake rate is proportional toseismicmoment rate, whichis proportional torelative platevelocity. To reducebiasbyaftershocksandswarms,weestimateindependenceprobabilitiesofearth- quakes and use them as weights. We assign shallow earthquakes to boundary steps and classes,thensortboundarystepswithineachclassbyvelocityandplotcumulativeearth- quakesagainstcumulativemodelmomentrate.Weusetwomeasuresofnonlinearityand 10 4 stochastic simulations to assess significance. In subduction zones the relationship betweenseismicityandvelocityisnonlinearwith99.9%confidence.Slowersubduction at

The measurements of the parameters of split shear (S) waves from local deep-focus earthquakes recorded in 2005–2007 by a network of 12 seismic stations in Southern Sakhalin are presented. The results revealed the heterogeneous distribution of the anisotropic properties beneath Southern Sakhalin. The azimuths of the fast S-wave polarization beneath the stations in the central part of the peninsula are oriented along the NNW and NNE-NE directions normal to and along the Kuril Trench. Beneath the stations located along the western and eastern coasts, the azimuths of the fast S-wave polarization change their direction from NNW in the northern area to E-SE in the southern area. The highest anisotropy degree (up to 0.9–1.5%) is recorded beneath the central part of Southern Sakhalin. The maximum values of the discrepancy in the arrival time of the split S-waves are observed when the azimuth of the fast S-wave is oriented along the NNE beneath the active fault zones. The analysis of the variations of the S-wave lag time shows their weak depth dependence. The highest anisotropy is assumed in the upper layers of the medium (down to a depth of about 250 km). The results obtained for the dominating wave frequency of 1–5 Hz represent mainly the medium-scale anisotropy of the top of the studied region.

New Zealand is one of the more seismically active countries in the world, with more than 15,000 earthquakes located each year. Routine moment tensor analysis of regional seismic data for earthquakes with moment magnitude M w>∼3.5–4.0 has recently been implemented in New Zealand. Nearly 330 regional moment tensor (RMT) solutions have been calculated for earthquakes in the New Zealand region dating back to late 2003. This complements local magnitude ( M L), the primary magnitude calculated by GeoNet in New Zealand. The RMT catalog, along with 155 Global Centroid Moment Tensor (Global CMT) catalog solutions, is used to compare M w with M L for New Zealand earthquakes. In addition to M w and M L, there are more than 330 teleseismic body-wave magnitude ( m b) values available from the United States Geological Survey and International Seismological Center catalogs for events that also have an M w. These are used to examine the relationship between M w, M L, and m b for New Zealand earthquakes. There is a clear distinction in the relationship of M L to M w for shallow focus (≤33 km depth) and deep focus (>33 km depth) earthquakes. Shallow focus earthquakes show M L to be fairly consistent with M w, particularly for events with M w>∼4.5 and with ![Graphic][1] . Deep focus earthquakes have M L consistently larger than M w (more than a full magnitude unit for some events) with ![Graphic][2] . M w and m b are in fairly good agreement regardless of the depth, whereas m b estimates are consistently smaller than M L for deep events. This suggests that M L is overestimated for deep focus events and that the large M Ls are not the result of physical properties of the source. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif

SUMMARY We investigate the 2007 Pisco, Peru, earthquake (M8.0) using seismological and geodetic data. Analysis of teleseismic records indicates that the earthquake is characterized by a steadily increasing moment rate, with the maximum occurring between about 55 and 70 s and the associated slip occurring about 40 km south of the epicentre. We combine InSAR data from L-band (ALOS) and C-band (Envisat) satellites to produce a 3-D image of the ground displacements caused by the earthquake. Distributed slip inversions of the interferograms identify a single large patch of slip in the location of the maximum moment rate identified by seismology. Despite a clear spatial coincidence between the location of the coseismic slip and the anomalously high topography of the Paracas Peninsula, we find no evidence that the earthquake was directly responsible for motion on any upper plate faults.

A compiled gravity anomaly map of the Western Himalayan Syntaxis is analysed to understand the tectonics of the region around the epicentre of Kashmir earthquake of October 8, 2005 (Mw = 7.6). Isostatic gravity anomalies and effective elastic thickness (EET) of lithosphere are assessed from coherence analysis between Bouguer anomaly and topography. The isostatic residual gravity high and gravity low correspond to the two main seismic zones in this region, viz. Indus–Kohistan Seismic Zone (IKSZ) and Hindu Kush Seismic Zones (HKSZ), respectively, suggesting a connection between siesmicity and gravity anomalies. The gravity high originates from the high-density thrusted rocks along the syntaxial bend of the Main Boundary Thrust and coincides with the region of the crustal thrust earthquakes, including the Kashmir earthquake of 2005. The gravity low of HKSZ coincides with the region of intermediate–deep-focus earthquakes, where crustal rocks are underthrusting with a higher speed to create low density cold mantle. Comparable EET (∼55 km) to the focal depth of crustal earthquakes suggests that whole crust is seismogenic and brittle. An integrated lithospheric model along a profile provides the crustal structure of the boundary zones with crustal thickness of about 60 km under the Karakoram–Pamir regions and suggests continental subduction from either sides (Indian and Eurasian) leading to a complex compressional environment for large earthquakes.

We present an acoustic emission (AE) monitoring technique to study high-pressure (P > 1 GPa) microseismicity in multi-anvil rock deformation experiments. The application of this technique is aimed at studying fault mechanisms of deep-focus earthquakes that occur during subduction at depths up to 650 km. AE monitoring in multi-anvil experiments is challenging because source locations need to be resolved to a submillimetre scale due to the small size of the experimental assembly. AEs were collected using an 8-receiver array, located on the back truncations of the tungsten carbide anvils. Each receiver consists of a 150-1000 kHz bandwidth PZT transducer assembly. Data were recorded and processed using a high-speed AMSY-5 acquisition system from Vallen-Systems, allowing waveform collection at a 10 MHz sampling rate for each event signal. 3-D hypocentre locations in the assembly are calculated using standard seismological algorithms. The technique was used to monitor fault development in 3 mm long x 1.5 mm diameter olivine cores during axisymmetric compression and extension. The faults were generated during cold compression to similar to 2 GPa confining pressure. Subsequent AEs at 2-6 GPa and 900 degrees C were found to locate near these pre-existing faults and exhibit high pressure stick-slip behaviour.

The Kuril-Kamchatka seismofocal zone was thought to be a single plate approximately 90 km wide and dipping to a depth of 700 km at an angle of 40°–45°. This concept reflects primarily the physical differences (elastic wave velocities, density, temperature, etc.) between the seismofocal zone and the mantle hosting it. Detailed investigations show that the seismofocal zone proper is also heterogeneous with earthquake hypocenters variably concentrated and clustered within this zone, where both seismogenic and aseismic strata, as well as subvertical zones, can be identified. The latter are reflected in the structure and faults of the Earth’s crust and upper mantle.

SUMMARY We present results from forward modelling to study the feasibility of using S-to-P converted phases to image the seismic discontinuity structure of the crust and upper mantle. We show that a significant level of P-wave energy arriving before the direct S-wave arrival can interfere with the S-to-P converted phases of interest and may result in Sp receiver function phases that do not represent true earth structure. The source of this P-wave energy is attributable to a number of phases, including those that have undergone multiple reflections off the Earth’s surface. For deep focus earthquakes (300‐600 km deep), a significant amount of P-wave energy is observed from pPPP, pPPPP and sPPPP phases, and arrives within the same time window as predicted for S-to-P converted phases from the direct S phase arrival. Furthermore, for earthquakes at all depths, interfering P-wave energy arrives within the same time window as predicted for S-to-P converted phases from the SKS phase arrival, limiting the usefulness of SKSp receiver functions for upper mantle imaging. To isolate true Sp receiver function phases from contamination due to other P-wave phases, we find it necessary to stack receiver functions from a range of epicentral distances and depths in order to aid the suppression of noise and other unwanted phases. We provide constraints on the noise levels to be expected as a function of epicentral distance and earthquake depth. We find that the lowest noise levels are achievable by restricting epicentral distance to less than 75 degrees and the depth of earthquakes used to less than 300 km.

Data of the intermediate depth (the Geioyo and the Shizuoka) earthquakes in Japan recorded in a dense network is analysed and compared with various available attenuation relations. The approach of Midorikawa (Midorikawa S., 1993, Technophysics 218, 287–295) based on the empirical Green’s function technique of Irikura (Irikura, K., 1986, Proceedings of the 7th Japan Earthquake Engineering Symposium, pp. 151–156.) has been used to model the rupture responsible for these earthquakes and peak ground acceleration are simulated at selected observation points. The method presented in this work includes the transmission effect in a multiple layer crustal model for a finite fault earthquake source model. Sharp attenuation rate is seen for such intermediate depth earthquakes which is difficult to explain through conventional attenuation relations. Detail study of the methodology and comparison of results shows that the transmission factor plays an important role for the sharp attenuation rate for intermediate-to deep-focus earthquakes.

Between October 12, 1989 and December 31, 2002, the Alaska Volcano Observatory (AVO) located 162 deep long-period (DLP) events beneath 11 volcanic centers in the Aleutian arc. These events generally occur at mid- to lower-crustal depths (10–45 km) and are characterized by emergent phases, extended codas, and a strong spectral peak between 1.0 and 3.0 Hz. Observed wave velocities and particle motions indicate that the dominant phases are P- and S-waves. DLP epicenters often extend over broad areas (5–20 km) surrounding the active volcanoes. The average reduced displacement of Aleutian DLPs is 26.5 cm2 and the largest event has a reduced displacement of 589 cm2 (or ML 2.5). Aleutian DLP events occur both as solitary events and as sequences of events with several occurring over a period of 1–30 min. Within the sequences, individual DLPs are often separated by lower-amplitude volcanic tremor with a similar spectral character. Occasionally, volcano-tectonic earthquakes that locate at similar depths are contained within the DLP sequences. At most, Aleutian volcanoes DLPs appear to loosely surround the main volcanic vent and occur as part of background seismicity. A likely explanation is that they reflect a relatively steady-state process of magma ascent over broad areas in the lower and middle portions of the crust. At Mount Spurr, DLP seismicity was initiated by the 1992 eruptions and then slowly declined until 1997. At Shishaldin Volcano, a short-lived increase in DLP seismicity occurred about 10 months prior to the April 19, 1999 eruption. These observations suggest a link between eruptive activity and magma flux in the mid- to lower-crust and uppermost mantle.

We installed five broadband seismometers around Mt. Usu, Hokkaido, Japan, just before the first surface eruption on March 31, 2000. By using these broadband data with short-period seismograms recorded by Japan Meteorological Agency and Institute of Seismology and Volcanology of Hokkaido University, we located 590 earthquakes associated with the 2000 eruption of Mt. Usu from March 31, 12:58 to April 5, 2000, including 41 low-frequency earthquakes. Low-frequency earthquakes have clear predominant frequencies of about 1.0–1.5 Hz. While the seismicity of tectonic earthquakes was still active on April 5, there were almost no low-frequency events after April 2. All the low-frequency earthquakes occurred within a vertical and narrow zone near the craters of eruption, and they are shallower than 4 km. In contrast, the active area of tectonic earthquakes spreads out in a large area, particularly deep in the south of the craters. From these temporal and spatial patterns, together with GPS and gravity measurements that support the termination of magma activities before April 2, we suggest that low-frequency earthquakes were caused by magma or hydrothermal activities beneath the volcano while tectonic ones were caused by the regional tectonic stress affected by the intrusion of magma body near the craters. Related to long-term magma activities of this eruption, four deep crustal low-frequency events were identified on October 17, 1998 just in the south of Mt. Usu in the depth range of 20–30 km. As proposed for other active volcanoes, these deep crustal events may represent a part of deep magma activities beneath Mt. Usu.

Quantitative estimates of the crack density (ϵ), saturation rate (ξ) and porosity (ψ) parameters from seismic velocities (Vp, Vs) and Poisson’s ratio (σ) in the 2001 Bhuj earthquake area in western India indicate that the 2001 Bhuj earthquake hypocenter is associated with high-ϵ, high-ξ and high-ψ in the depth range of about 23–28 km, extending 15–30 km laterally. These anomalies may be due to a fluid-filled, fractured rock matrix, which might have contributed to trigger the 2001 Bhuj earthquake in the intraplate, stable continental region of the Indian Peninsula. This feature is similar to that of the 1995 Kobe earthquake [D. Zhao et al., Science 274 (1996) 1891–1894]. High-ψ areas are generally consistent with high-ϵ areas, but high-ξ areas have a wider distribution, indicating that micro-cracks exist in more localized areas within the crust, and that permeation of fluids in the hypocenter zone might have occurred extensively through the intergranular and fractured pores due to hidden intersecting fault geometry. Here we suggest the possibility that earthquake occurrence is closely related to in situ material heterogeneities, rather than stress conditions alone.

Several destructive earthquakes have occurred in Tien-Shan region at the beginning of 20th century. However, the detailed seismological characteristics, especially source parameters of those earthquakes are still poorly investigated. The Chon-Kemin earthquake is the strongest instrumentally recorded earthquake in the Tien-Shan region. This earthquake has produced an approximately 200 km long system of surface ruptures along Kemin-Chilik fault zone and killed about similar to 400 people. Several studies presented the different information on the earthquake epicentre location and magnitude, and two different focal mechanisms were also published. The reason for the limited knowledge of the source parameters for the Chon-Kemin earthquake is the complexity of old analogue records processing, digitization and analysis. In this study the data from 23 seismic stations worldwide were collected and digitized. The earthquake epicentre was relocated to 42.996NA degrees and 77.367EA degrees, the hypocentre depth is estimated between 10 and 20 km. The magnitude was recalculated to m(B) 8.05, M-s 7.94 and M-w 8.02. The focal mechanism, determined from amplitude ratios comparison of the observed and synthetic seismograms, was: str = 264A degrees, dip = 52A degrees, rake = 98A degrees. The apparent source time duration was between similar to 45 and similar to 70 s, the maximum slip occurred 25 s after the beginning of the rupture. Two subevents were clearly detected from the waveforms with the scalar moment ratio between them of about 1/3, the third subevent was also detected with less certainty. Taking into account surface rupture information, the fault geometry model with three patches was proposed. Based on scaling relations we conclude that the total rupture length was between similar to 260 and 300 km and a maximum rupture width could reach similar to 70 km.

It has been proposed that intermediate- and deep-focus earthquakes occur on preexisting planes of weakness that were created at shallow depth. The purpose of this study is to test this hypothesis. We examined fault plane solutions for 360 events that occurred in the Tonga subduction zone from January 1976 to February 1999. For events in the slab, each fault plane solution is rotated by the slab dip angle so that the slab is in its original horizontal position. We then compare fault plane solutions for events in the slab with those for events at the outer rise. We find that an asymmetric fault system, corresponding to that found for outer rise events, persists down to about 450 km depth. This suggests that earthquakes down to this depth are caused by the reactivation of preexisting faults created in the oceanic plate before subduction. A similar pattern, though less well constrained, is found in the Kurile subduction zone down to 200 km depth. For earthquakes deeper than 450 km the fault plane solutions are much more scattered. In some places the fault plane pattern suggests the creation of new fault planes along orientations of maximum shear with zero internal friction. The pattern of fault plane solutions in particular localized volumes cannot be easily explained by the creation of new faults in the direction of maximum shear under a coherent regional tectonic stress field. Thus, either a highly heterogeneous stress field exists or preexisting planes of weakness also account for some deep earthquakes.

We investigate the nature of laterally varying P wave speeds (Vp) and their relationship to a group of unusual deep earthquakes beneath the backarc of Tonga. This subhorizontal swath of “outboard” earthquakes occurs within the transition zone of the mantle but is several hundred kilometers away from the inclined Wadati-Benioff zone. Surprisingly, the seismogenic region is not associated with fast Vp. It follows that the effect of cold temperature, as implied by the presence of large earthquakes, must have been counteracted by other factors such as compositional or mineralogical effects. Our results are based on a sequence of seismic profiles in a fan-shot geometry. For a given azimuthal sector, we seek the simplest model of radially varying Vp that explains observed broadband waveforms and travel times. The lateral extent of anomalous Vp is then constrained by comparing the results from a wide range of azimuths. This approach fully utilizes information from triplicate P arrivals that are sensitive to seismic wave speeds at the depths of interest. Our results identify prominent anomalies of fast Vp (∼3% over average Earth models) both in the transition zone and in the lower mantle beneath the Tonga backarc. However, the region of fast Vp in the transition zone is offset by ∼300 km to the south of that in the lower mantle, and the region of offset coincides with the outboard earthquakes. A straightforward interpretation is that a large flux of subducted Pacific lithosphere lingers in the transition zone before descending into the lower mantle. Within this broad remnant of subducted lithosphere, the outboard earthquakes occur in a region where impounding of subducted material is most significant, a condition favorable for the accumulation of met astable olivine or volatiles. Such materials reduce Vp and are candidates for triggering earthquakes, thus accounting for the absence of fast Vp in the source zone of deep earthquakes.

—We analyze the seismic moment-frequency relation in various depth ranges and for different seismic regions, using Flinn-Engdahl's regionalization of global seismicity. Three earthquake lists of centroid-moment tensor data have been used the Harvard catalog, the USGS catalog, and the Huang et al. (1997) catalog of deep earthquakes. The results confirm the universality of the β-values and the maximum moment for shallow earthquakes in continental regions, as well as at and near continental boundaries. Moreover, we show that although fluctuations in earthquake size distribution increase with depth, the β-values for earthquakes in the depth range of 0–500 km exhibit no statistically significant regional variations. The regional variations are significant only for deep events near the 660 km boundary. For declustered shallow earthquake catalogs and deeper events, we show that the worldwide β-values have the same value of 0.60 ± 0.02. This finding suggests that the β-value is a universal constant. We investigate the statistical correlations between the numbers of seismic events in different depth ranges and the correlation of the tectonic deformation rate and seismic activity (the number of earthquakes above a certain threshold level per year). The high level of these correlations suggests that seismic activity indicates tectonic deformation rate in subduction zones. Combined with the universality of the β-value, this finding implies little if any variation in maximum earthquake seismic moment among various subduction zones. If we assume that earthquakes of maximum size are similar in different depth ranges and the seismic efficiency coefficient, χ, is close to 100% for shallow seismicity, then we can estimate χ for deeper earthquakes for intermediate earthquakes χ≈ 5%, and χ≈ 1% for deep events. These results may lead to new theoretical understanding of the earthquake process and better estimates of seismic hazard.

We analyze the coda of teleseismic P waves to deterministically map seismic scatterers in a 1000 km on-a-side cube beneath the Izu-Bonin trench by migrating and stacking waveforms. This method was applied to several hundred short-period vertical-component western United States seismometer recordings of 17 deep-focus Izu-Bonin earthquakes. Except for isolated arrivals, the midmantle appears devoid of sharp discontinuities. S-to-P conversions at the 660-km discontinuity generated the largest signals; the 410-km discontinuity occasionally generated signals. Horizontal discontinuities poorly explain other signals; however, S-to-P conversions generated at a nearly north-south trending, steeply dipping discontinuity at 30°N, 145°E, and 1000 km beneath and parallel to the Izu-Bonin trench explain most additional signals. Observed coherent S-to-P conversions of 2-s period waves limit the width of a gradient zone to