Inhomogeneous distribution of deep slow earthquake activity along the strike of the subducting Philippine Sea Plate

The spatial distribution of deep slow earthquake activity along the strike of the subducting Philippine Sea Plate in southwest Japan is investigated. These events usually occur simultaneously between the megathrust seismogenic zone and the deeper free-slip zone on the plate interface at depths of about 30 km. Deep low-frequency tremors are weak prolonged vibrations with dominant frequencies of 1.5–5 Hz, whereas low-frequency earthquakes correspond to isolated pulses included within the tremors. Deep very-low-frequency earthquakes have long-period (20 s) seismic signals, and short-term slow-slip events are crustal deformations lasting for several days. Slow earthquake activity is not spatially homogeneous but is separated into segments some of which are bounded by gaps in activity. The spatial distribution of each phase of slow earthquake activity is usually coincident, although there are some inconsistencies. Very-low-frequency earthquakes occur mainly at edges of segments. Low-frequency earthquakes corresponding to tremors of relatively large amplitude are concentrated at spots where tremors are densely distributed within segments. The separation of segments by gaps suggests large differences in stick-slip and stable sliding caused by frictional properties of the plate interface. Within each segment, variations in the spatial distribution of slow earthquakes reflected inhomogeneities corresponding to the characteristic scales of events.     

Non-volcanic deep low-frequency tremors accompanying slow slips in the southwest Japan subduction zone

Non-volcanic deep low-frequency tremors in southwest Japan exhibit a strong temporal and spatial correlation with slow slip detected by the dense seismic network. The tremor signal is characterized by a low-frequency vibration with a predominant frequency of 0.5–5 Hz without distinct P- or S-wave onset. The tremors are located using the coherent pattern of envelopes over many stations, and are estimated to occur near the transition zone on the plate boundary on the forearc side along the strike of the descending Philippine Sea plate. The belt-like distribution of tremors consists of many clusters. In western Shikoku, the major tremor activity has a recurrence interval of approximately six months, with each episode lasting over a week. The tremor source area migrates during each episode along the strike of the subducting plate with a migration velocity of about 10 km/day. Slow slip events occur contemporaneously with this tremor activity, with a coincident estimated source area that also migrates during each episode. The coupling of tremor and slow slip in western Shikoku is very similar to the episodic tremor and slip phenomenon reported for the Cascadia margin in northwest North America. The duration and recurrence interval of these episodes varies between tremor clusters even on the same subduction zone, attributable to regional difference in the frictional properties of the plate interface.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Plate subduction,and generation of earthquakes and magmas in Japan as inferred from seismic observations: An overview

A dense nationwide seismic network recently constructed in Japan has been yielding large volumes of high-quality data that have made it possible to investigate the seismic structure in the Japanese subduction zone with unprecedented resolution. In this article, recent studies on the subduction of the Philippine Sea and Pacific plates beneath the Japanese Islands and the mechanism of earthquake and magma generation associated with plate subduction are reviewed. Seismic tomographic studies have shown that the Philippine Sea plate subducting beneath southwest Japan is continuous throughout the entire region, from Kanto to Kyushu, without disruption or splitting even beneath the Izu Peninsula as suggested in the past. The contact of the Philippine Sea plate with the Pacific plate subducting below has been found to cause anomalously deep interplate and intraslab earthquake activity in Kanto. Detailed waveform inversion studies have revealed that the asperity model is applicable to interplate earthquakes. Analyses of dense seismic and GPS network data have confirmed the existence of episodic slow slip accompanied in many instances by low-frequency tremors/earthquakes on the plate interface, which are inferred to play an important role in stress loading at asperities. High-resolution studies of the spatial variation of intraslab seismicity and the seismic velocity structure of the slab crust strongly support the dehydration embrittlement hypothesis for the generation of intraslab earthquakes. Seismic tomography studies have shown that water released by dehydration of the slab and secondary convection in the mantle wedge, mechanically induced by slab subduction, are responsible for magma generation in the Japanese islands. Water of slab origin is also inferred to be responsible for large anelastic local deformation of the arc crust leading to inland crustal earthquakes that return the arc crust to a state of spatially uniform deformation.     

Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, Japan

We found a characteristic space–time pattern of the tidal triggering effect on earthquake occurrence in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, central Japan, where a large interplate earthquake may be impending. We measured the correlation between the Earth tide and earthquake occurrence using microearthquakes that took place in the Philippine Sea plate for about two decades. For each event, we assigned the tidal phase angle at the origin time by theoretically calculating the tidal shear stress on the fault plane. Based on the distribution of the tidal phase angles, we statistically tested whether they concentrate near some particular angle or not by using Schuster's test. In this test, the result is evaluated by p-value, which represents the significance level to reject the null hypothesis that earthquakes occur randomly irrespective of the tidal phase angle. As a result of analysis, no correlation was found for the data set including all the earthquakes. However, we found a systematic pattern in the temporal variation of the tidal effect; the p-value significantly decreased preceding the occurrence of M ≥ 4.5 earthquakes, and it recovered a high level afterwards. We note that those M ≥ 4.5 earthquakes were considerably larger than the normal background seismicity in the study area. The frequency distribution of tidal phase angles in the pre-event period exhibited a peak at the phase angle where the tidal shear stress is at its maximum to accelerate the fault slip. This indicates that the observed small p-value is a physical consequence of the tidal effect. We also found a distinctive feature in the spatial distribution of p-values. The small p-values appeared just beneath the strongly coupled portion of the plate interface, as inferred from the seismicity rate change in the past few years.     

Influence of interaction between small asperities on various types of slow earthquakes in a 3-D simulation for a subduction plate boundary

Recently, the occurrence of slow earthquakes such as low-frequency earthquakes and very low-frequency earthquakes have been recognized at depths of about 30 km in southwest Japan and Cascadia. These slow earthquakes occur sometimes in isolation and sometimes break into chain-reaction, producing tremor that migrates at a speed of about 5–15 km/day and suggesting a strong interaction among nearby small asperities. In this study, we formulate a 3-D subduction plate boundary model with two types of small asperities chained along the trench at the depth of 30 km. Our simulation succeeds in representing various types of slow earthquakes including low-frequency earthquakes and rapid slip velocity in the same asperity, and indicates that interaction between asperities may cause the very low-frequency earthquakes. Our simulation also shows chain reaction along trench with propagation speed that can be made consistent with observations by adjusting model parameters, which suggests that the interactions also explain the observed migration of slow earthquakes.     

Recurrence of great earthquakes at subduction zones

Statistics of the recurrence times of great earthquakes at the Pacific subduction margins are made. The mean return period of great earthquakes is different from zone to zone, ranging from 27 to 117 years. The standard deviation of the return period proves to be very small, several years say, in some cases. The probabilities of a great earthquake recurring in each zone are estimated on the basis of Weibull distribution analysis.The mean return periods thus estimated are combined with the relative plate velocities at respective zones as obtained in the plate tectonics in order to estimate the ultimate displacement to rupture at the interface of the continental plate and the downgoing oceanic plate. It is presumed that great earthquakes at subduction zones occur as a result of a rebound of the continental plate at the time of rupture. The ultimate displacement thus estimated ranges from 2 to 8 m, and seems somewhat larger than that estimated on the basis of seismic observations, although the value of ultimate displacement seems to harmonize roughly with estimates based on geodetic observations on land. However, the ultimate displacement at the Aleutian—Alaska zone as estimated here seems much smaller than that estimated from actual observations.The ultimate strains, which are deduced from the displacements obtained on the assumption that the logarithmic extent of the deformed area is proportional to earthquake magnitude, are then calculated, and compared with those estimated for large inland earthquakes as revealed by repetition of geodetic surveys. The mean ultimate strain is estimated as 4.3 · 10⁻⁵ for subduction-zone earthquakes while that for inland earthquakes has been estimated as 4.7 · 10⁻⁵. As the agreement between both the ultimate strains is fairly good, it is tentatively concluded that the strength of the plate interface under the sea bottom is more or less the same as that in the crust on land.     

Repeating earthquake activities associated with the Philippine Sea plate subduction in the Kanto district, central Japan: A new plate configuration revealed by interplate aseismic slips

We detect repeating earthquakes associated with the Philippine Sea plate subduction to reveal the plate configuration. In the Kanto district, we find 140 repeating earthquake groups with 428 events by waveform similarity analysis. Most repeating earthquakes in the eastern part of the Kanto district occur with a regular time interval. They have thrust-type focal mechanisms and are distributed near the upper surface of the Philippine Sea plate. These observations indicate that the repeating earthquakes there occur as a repetition of ruptures on the isolated patches distributed on the plate boundary owing to the concentration of stress caused by aseismic slips in the surrounding areas. This shows that the distributions of repeating earthquakes suggest the aseismic slips in the surrounding areas of small patches. We determine spatial distributions of repeating earthquakes in the eastern part of the Kanto district and find that they correspond to the upper boundary of the Philippine Sea plate, that is, the upper boundary of the oceanic crust layer of the Philippine Sea plate. The plate geometry around Choshi is newly constrained by repeating earthquake data and a rather flat geometry in the eastern part of the Kanto district is revealed. The obtained geometry suggests uplift of the Philippine Sea plate due to the collision with the Pacific plate beneath Choshi.Repeating earthquakes in the western part of the Kanto district have extremely shorter recurrence times, and their focal mechanisms are not of the thrust types. These repeating earthquakes are classified as “burst type” activity and likely to occur on the preexistent fault planes which are distributed around the “collision zone” between the Philippine Sea plate and the inland plate. The variation among the repeating earthquake activities in the Kanto district indicates that regular repetition of repeating earthquakes is possible only on the plate boundary with a smooth and simple geometry.     

Space-time seismicity and earthquake swarms: Certain observations along the slow-spreading mid-Indian Ocean ridges

The slow spreading mid-Indian Ocean ridge system containing the Carlsberg, Central and Southwest Indian ridges is seismically very active. In the present study, a detailed analysis has been carried out of the data of earthquake sources along different ridge segments in order to investigate the spatial and temporal clustering patterns and to evaluate crustal processes related to the swarm occurrences along these ridges. The spatial and temporal clustering pattern of the recent earthquakes (1980–1990) pertaining to nine major spreading segments and eight fracture zones suggests that the events cluster in greater proportion along the spreading segments than along the fracture zones. We performed a systematic search of earthquake catalogue during the period 1964–1990 by examining the spatio-temporal hypocentral clusters in order to identify the swarm occurrences along these ridges. The search included eighteen prominent sequences, of which, thirteen were earthquake swarms. Except two, all other swarms were found to be occurring mainly on the spreading segments. The maximum magnitude observed in these swarms is m_b = 5.4 and have many events predominantly showing normal faulting mechanisms. The spatial disposition and temporal activity of the events in swarms is much similar to the foreshock-mainshock-aftershock sequences observed along the spreading rift valley zones. These characteristics help us to support that swarms along the slow spreading mid-Indian Ocean ridges are the result of extensional tectonic activity, leading to the development of the median valley topography, a mechanism similar to that proposed by Bergman and Solomon (1990) for the Mid-Atlantic Ridge.     

阿尔金断裂带东段地区深浅部构造综合分析

阿尔金断裂带东段地区的地质构造特征及其动力学机制一直是地学工作者关注的焦点。近年来小震资料越来越多应用到活动断裂空间展布、深浅构造分析及动力学机制研究领域。本文应用双差定位法获得研究区域2008～2017年间6013次地震事件的精确定位数据,通过多条小震深度剖面清晰刻画出断裂系统的空间展布形态。综合石油地震剖面、人工地震宽角反射/折射剖面、人工地震深反射剖面,充分利用小震精确定位信息以及浅表活动构造研究成果,建立研究区断裂系统的深浅部构造模型。研究区莫霍面由北往南逐渐加深,存在三处断错,呈阶梯状展布,地壳内存在一条厚约10km的低速层,在该层以上为地震多发区,断裂系统总体呈"Y"字型,上部为一系列叠瓦状逆冲断裂,造成祁连山的隆升,向下并入一条主干断层。最后探讨了青藏高原东北缘地区构造运动的动力学机制,亚洲板块俯冲至祁连山前,上地壳以逆冲推覆构造模式造成上地壳增厚现象,而中下地壳主要为亚洲岩石圈地幔下插,上地幔的拖曳作用下发生流动引起地壳增厚,上下地壳整体增厚。     

京津冀地区历史地震事件时空特征研究*

受环太平洋地震带影响,华北平原地区地震频发,尤其是处于中国首都经济圈的京津冀地区的地震事件备受关注。通过对历史文献资料及地震台网记录中的地震事件统计、分析,重建该地区地震事件历史并获取其潜在的空间分布特征及时间规律,对未来地震事件的早期预警具有重要参考意义。分析结果表明,公元前231年至公元2018年期间京津冀地区发生的1044起地震事件中,以有感地震和中强地震为主,小地震、强烈地震以及大地震发生频次较低。地震记录完整性分析结果表明,除小地震外,其他等级地震记录自公元1400年以来基本完整。在空间分布上,京津冀地区历史地震呈“T”字形分布,沿1条北西—南东走向地震带和1条北东—南西走向地震带分布。在时间上,京津冀地区地震事件呈现出阶段性的变化,在公元1480—1680年间以及1950年以来2个时间段内较为活跃,发生频率较高,频谱分析结果进一步表明地震记录存在45年的复发周期。在月际尺度上,地震事件同样存在季节性差异且多发于夏秋季节,同时地震密集区域在年内呈现出自西向东迁移的现象。最后,根据历史地震事件发生的时间规律,在未来一段时间内京津冀地区仍将处于地震活跃期,存在发生强震的风险。     

Influence of interaction between small asperities on various types of slow earthquakes in a 3-D simulation for a subduction plate boundary

Recently, the occurrence of slow earthquakes such as low-frequency earthquakes and very low-frequency earthquakes have been recognized at depths of about 30 km in southwest Japan and Cascadia. These slow earthquakes occur sometimes in isolation and sometimes break into chain-reaction, producing tremor that migrates at a speed of about 5–15 km/day and suggesting a strong interaction among nearby small asperities. In this study, we formulate a 3-D subduction plate boundary model with two types of small asperities chained along the trench at the depth of 30 km. Our simulation succeeds in representing various types of slow earthquakes including low-frequency earthquakes and rapid slip velocity in the same asperity, and indicates that interaction between asperities may cause the very low-frequency earthquakes. Our simulation also shows chain reaction along trench with propagation speed that can be made consistent with observations by adjusting model parameters, which suggests that the interactions also explain the observed migration of slow earthquakes.     

Modeling a tsunami from the Nicoya,Costa Rica,seismic gap and its potential impact in Puntarenas

Although subduction zones around the world are known to be the source of earthquakes and/or tsunamis, not all segments of these plate boundaries generate destructive earthquakes and catastrophic tsunamis. Costa Rica, in Central America, has subduction zones on both the Pacific and the Caribbean coasts and, even though large earthquakes (Mw = 7.4–7.8) occur in these convergent margins, they do not produce destructive tsunamis. The reason for this is that the seismogenic zones of the segments of the subduction zones that produce large earthquakes in Costa Rica are located beneath land (Nicoya peninsula, Osa peninsula and south of Limón) and not off shore as in most subduction zones around the world. To illustrate this particularity of Costa Rican subduction zones, we show in this work the case for the largest rupture area in Costa Rica (under the Nicoya peninsula), capable of producing Mw ～ 7.8 earthquakes, but the tsunamis it triggers are small and present little potential for damage even to the largest port city in Costa Rica.The Nicoya seismic gap, in NW Costa Rica, has passed its ～50-year interseismic period and therefore a large earthquake will have to occur there in the near future. The last large earthquake, in 1950 generated a tsunami which slightly affected the southwest coast of the Nicoya Peninsula. We present here a simulation to study the possible consequences that a tsunami generated by the next Nicoya earthquake could have for the city of Puntarenas. Puntarenas has a population of approximately eleven thousand people and is located on a 7.5 km long sand bar with a maximum height of 2 m above the mean sea level. This condition makes Puntarenas vulnerable to tsunamis.     

Mechanism of the 2011 Tohoku-oki earthquake (Mw 9.0) and tsunami: Insight from seismic tomography

To clarify the generating mechanism of the 2011 Tohoku-oki earthquake (Mw 9.0) and the induced tsunami, we determined high-resolution tomographic images of the Northeast Japan forearc. Significant lateral variations of seismic velocity are visible in the megathrust zone, and most large interplate thrust earthquakes are found to occur in high-velocity (high-V) areas. These high-V zones may represent high-strength asperities at the plate interface where the subducting Pacific plate and the overriding Okhotsk plate are coupled strongly. A shallow high-V zone with large coseismic slip near the Japan Trench may account for the mainshock asperity of the 2011 Tohoku-oki earthquake. Because it is an isolated asperity surrounded by low-velocity patches, most stress on it was released in a short time and the plate interface became decoupled after the Mw 9.0 earthquake. Thus the overriding Okhotsk plate there was shot out toward the Japan Trench and caused the huge tsunami.     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

Seismic slip rates,potential subsurface rupture areas and seismic potential of the Vienna Basin Transfer Fault

The Vienna Basin Transfer Fault (VBTF) is a slow active fault with moderate seismicity (I _max~8–9, M _max~5.7) passing through the most vulnerable regions of Austria and Slovakia. We use different data to constrain the seismic potential of the VBTF including slip values computed from the seismic energy release during the 20th century, geological data on fault segmentation and a depth-extrapolated 3-D model of a generalized fault surface, which is used to define potential rupture zones. The seismic slip of the VBTF as a whole is in the range of 0.22–0.31 mm/year for a seismogenic fault thickness of 8 km. Seismic slip rates for individual segments vary from 0.00 to 0.77 mm/year. Comparing these data to geologically and GPS-derived slip velocities (>1 mm/year) proofs that the fault yields a significant seismic slip deficit. Segments of the fault with high seismic slip contrast from segments with no slip representing locked segments. Fault surfaces of segments within the seismogenic zone (4–14 km depth) vary from 55 to 400 km². Empirical scaling relations show that these segments are sufficiently large to explain both, earthquakes observed in the last centuries, and the 4th century Carnuntum earthquake, for which archeo-seismological data suggest a magnitude of M ≥ 6. Based on the combination of all data (incomplete earthquake catalog, seismic slip deficits, locked segments, potential rupture areas, indications of strong pre-catalog earthquakes) we argue, that the maximum credible earthquake for the VBTF is in the range M _max = 6.0–6.8, significantly larger than the magnitude of the strongest recorded events (M = 5.7).     

The 2003 earthquake swarm in Anjalankoski, south-eastern Finland

During May 2003 a swarm of 16 earthquakes (M_L = 0.6–2.1) occurred at Anjalankoski, south-eastern Finland. The activity lasted for three weeks, but additional two events were observed at the same location in October 2004. A comparison of the waveforms indicated that the source mechanisms and the hypocentres of the events were nearly identical.A relative earthquake location method was applied to better define the geometry of the cluster and to identify the fault plane associated with the earthquakes. The relocated earthquakes aligned along an ENE–WSW trending zone, with a lateral extent of about 1.0 km by 0.8 km. The relative location and the waveform-modelling of depth sensitive surface wave (Rg) and S-to-P converted body wave (sP) phases indicated that the events were unusually shallow, most likely occurring within the first 2 km of the surface. The revised historical earthquake data confirm that shallow swarm-type seismicity is characteristic to the area.The focal mechanism obtained as a composite solution of the five strongest events corresponds to dip-slip motion along a nearly vertical fault plane (strike = 250°, dip = 80°, rake = 90°). The dip and strike of this nodal plane as well as the relocated hypocentres coincide with an internal intrusion boundary of the Vyborg rapakivi batholith.The events occur under a compressive local stress field, which is explained by large gravitational potential energy differences and ridge-push forces. Pore-pressure changes caused by intrusion of ground water and/or radon gas into the fracture zones are suggested to govern the swarm-type earthquake activity.     

Earthquake generation in Cyprus revealed by the evolving stress field

Strong earthquake occurrence (M ≥ 6.0) onshore and offshore the Cyprus Island constitutes significant seismic hazard because they occur close to populated areas. Seismicity is weak south of the Island along the Cyprean Arc and strong events are aligned along the Paphos transform fault and Larnaka thrust fault zone that were already known and the Lemessos thrust fault zone that defined in the present study. By combining the past history of strong (M ≥ 6.0) events and the long-term tectonic loading on these major fault zones, the evolution of the stress field from 1896 until the present is derived. Although uncertainties exist in the location, magnitude and fault geometries of the early earthquakes included in our stress evolutionary model, the resulting stress field provides an explanation of later earthquake triggering. It was evidenced that the locations of all the strong events were preceded by a static stress change that encouraged failure. The current state of the evolved stress field may provide evidence for the future seismic hazard. Areas of positive static stress changes were identified in the southwestern offshore area that can be considered as possible sites of future seismic activity.     

Flexure of the Indian plate and intraplate earthquakes

The flexural bulge in central India resulting from India's collision with Tibet has a wavelength of approximately 670 km. It is manifest topographically and in the free-air gravity anomaly and the geoid. Calculations of the stress distribution within a flexed Indian plate reveal spatial variations throughout the depth of the plate and also a function of distance from the Himalaya. The wavelength (and therefore local gradient) of stress variation is a function of the effective elastic thickness of the plate, estimates of which have been proposed to lie in the range 40–120 km. The imposition of this stress field on the northward moving Indian plate appears fundamental to explaining the current distribution of intraplate earthquakes and their mechanisms. The current study highlights an outer trough south of the flexural bulge in central India where surface stresses are double the contiguous compressional stresses to the north and south. The Bhuj, Latur and Koyna earthquakes and numerous other recent reverse faulting events occurred in this compressional setting. The N/S spatial gradient of stress exceeds 2 bars/km near the flexural bulge. The overall flexural stress distribution provides a physical basis for earthquake hazard mapping and suggests that areas of central India where no historic earthquakes are recorded may yet be the locus of future damaging events.     

Seismic rupture during the 1960 Great Chile and the 2010 Maule earthquakes limited by a giant Pleistocene submarine slope failure

Determining factors that limit coseismic rupture is important to evaluate the hazard of powerful subduction zone earthquakes such as the 2011 Tohoku‐Oki event (Mw = 9.0). In 1960 (Mw = 9.5) and 2010 (Mw = 8.8), Chile was hit by such powerful earthquakes, the boundary of which was the site of a giant submarine slope failure with chaotic debris subducted to seismogenic zone depth. Here, a continuous décollement is absent, whereas away from the slope failure, a continuous décollement is seismically imaged. We infer that underthrusting of inhomogeneous slide deposits prevents the development of a décollement, and thus the formation of a thin continuous slip zone necessary for earthquake rupture propagation. Thus, coseismic rupture during the 1960 and 2010 earthquakes seems to be limited by underthrusted upper plate mass‐wasting deposits. More generally, our results suggest that upper plate dynamics and resulting surface processes can play a key role for determining rupture size of subduction zone earthquakes.