Delayed triggering of the 1999 Hector Mine earthquake by viscoelastic stress transfer

Stress changes in the crust due to an earthquake can hasten the failure of neighbouring faults and induce earthquake sequences in some cases. The 1999 Hector Mine earthquake in southern California (magnitude 7.1) occurred only 20 km from, and 7 years after, the 1992 Landers earthquake (magnitude 7.3). This suggests that the Hector Mine earthquake was triggered in some fashion by the earlier event. But uncertainties in the slip distribution and rock friction properties associated with the Landers earthquake have led to widely varying estimates of both the magnitude and sign of the resulting stress change that would be induced at the location of the Hector Mine hypocentre-with estimates varying from -1.4 bar (ref. 6) to +0.5 bar (ref. 7). More importantly, coseismic stress changes alone cannot satisfactorily explain the delay of 7 years between the two events. Here we present the results of a three-dimensional viscoelastic model that simulates stress transfer from the ductile lower crust and upper mantle to the brittle upper crust in the 7 years following the Landers earthquake. Using viscoelastic parameters that can reproduce the observed horizontal surface deformation following the Landers earthquake, our calculations suggest that lower-crustal or upper-mantle flow can lead to postseismic stress increases of up to 1-2 bar at the location of the Hector Mine hypocentre during this time period, contributing to the eventual occurrence of the 1999 Hector Mine earthquake. These results attest to the importance of considering viscoelastic processes in the assessment of seismic hazard.     

Remote triggering of deep earthquakes in the 2002 Tonga sequences

It is well established that an earthquake in the Earth's crust can trigger subsequent earthquakes, but such triggering has not been documented for deeper earthquakes. Models for shallow fault interactions suggest that static (permanent) stress changes can trigger nearby earthquakes, within a few fault lengths from the causative earthquake, whereas dynamic (transient) stresses carried by seismic waves may trigger earthquakes both nearby and at remote distances. Here we present a detailed analysis of the 19 August 2002 Tonga deep earthquake sequences and show evidence for both static and dynamic triggering. Seven minutes after a magnitude 7.6 earthquake occurred at a depth of 598 km, a magnitude 7.7 earthquake (664 km depth) occurred 300 km away, in a previously aseismic region. We found that nearby aftershocks of the first mainshock are preferentially located in regions where static stresses are predicted to have been enhanced by the mainshock. But the second mainshock and other triggered events are located at larger distances where static stress increases should be negligible, thus suggesting dynamic triggering. The origin times of the triggered events do not correspond to arrival times of the main seismic waves from the mainshocks and the dynamically triggered earthquakes frequently occur in aseismic regions below or adjacent to the seismic zone. We propose that these events are triggered by transient effects in regions near criticality, but where earthquakes have difficulty nucleating without external influences.     

Tracking the rupture of the Mw = 9.3 Sumatra earthquake over 1,150 km at teleseismic distance

On 26 December 2004, a moment magnitude Mw = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region. The rapid and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave-height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN), located approximately 9,000 km from the event. Here we present a modification of a standard array-seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.     

Predicting the endpoints of earthquake ruptures

The active fault traces on which earthquakes occur are generally not continuous, and are commonly composed of segments that are separated by discontinuities that appear as steps in map-view. Stress concentrations resulting from slip at such discontinuities may slow or stop rupture propagation and hence play a controlling role in limiting the length of earthquake rupture. Here I examine the mapped surface rupture traces of 22 historical strike-slip earthquakes with rupture lengths ranging between 10 and 420 km. I show that about two-thirds of the endpoints of strike-slip earthquake ruptures are associated with fault steps or the termini of active fault traces, and that there exists a limiting dimension of fault step (3-4 km) above which earthquake ruptures do not propagate and below which rupture propagation ceases only about 40 per cent of the time. The results are of practical importance to seismic hazard analysis where effort is spent attempting to place limits on the probable length of future earthquakes on mapped active faults. Physical insight to the dynamics of the earthquake rupture process is further gained with the observation that the limiting dimension appears to be largely independent of the earthquake rupture length. It follows that the magnitude of stress changes and the volume affected by those stress changes at the driving edge of laterally propagating ruptures are largely similar and invariable during the rupture process regardless of the distance an event has propagated or will propagate.     

Analysing the 1811-1812 New Madrid earthquakes with recent instrumentally recorded aftershocks

Although dynamic stress changes associated with the passage of seismic waves are thought to trigger earthquakes at great distances, more than 60 per cent of all aftershocks appear to be triggered by static stress changes within two rupture lengths of a mainshock. The observed distribution of aftershocks may thus be used to infer details of mainshock rupture geometry. Aftershocks following large mid-continental earthquakes, where background stressing rates are low, are known to persist for centuries, and models based on rate-and-state friction laws provide theoretical support for this inference. Most past studies of the New Madrid earthquake sequence have indeed assumed ongoing microseismicity to be a continuing aftershock sequence. Here we use instrumentally recorded aftershock locations and models of elastic stress change to develop a kinematically consistent rupture scenario for three of the four largest earthquakes of the 1811-1812 New Madrid sequence. Our results suggest that these three events occurred on two contiguous faults, producing lobes of increased stress near fault intersections and end points, in areas where present-day microearthquakes have been hitherto interpreted as evidence of primary mainshock rupture. We infer that the remaining New Madrid mainshock may have occurred more than 200 km north of this region in the Wabash Valley of southern Indiana and Illinois--an area that contains abundant modern microseismicity, and where substantial liquefaction was documented by historic accounts. Our results suggest that future large mid-plate earthquake sequences may extend over a much broader region than previously suspected.     

Evidence for fault weakness and fluid flow within an active low-angle normal fault.

Determining the composition and physical properties of shallow-dipping, active normal faults (dips < 35 degrees with respect to the horizontal) is important for understanding how such faults slip under low resolved shear stress and accommodate significant extension of the crust and lithosphere. Seismic reflection images and earthquake source parameters show that a magnitude 6.2 earthquake occurred at about 5 km depth on or close to a normal fault with a dip of 25-30 degrees located ahead of a propagating spreading centre in the Woodlark basin. Here we present results from a genetic algorithm inversion of seismic reflection data, which shows that the fault at 4-5 km depth contains a 33-m-thick layer with seismic velocities of about 4.3 km s(-1), which we interpret to be composed of serpentinite fault gouge. Isolated zones exhibit velocities as low as approximately 1.7 km s(-1) with high porosities, which we suggest are maintained by high fluid pressures. We propose that hydrothermal fluid flow, possibly driven by a deep magmatic heat source, and high extensional stresses ahead of the ridge tip have created conditions for fault weakness and strain localization on the low-angle normal fault.     

全球大地震破裂空间复杂度特征研究

基于有限断层反演得到的滑移数据,用独立滑移单元的数量表征地震破裂复杂度,据此对大地震破裂进行分组,并研究破裂复杂度与主要震源参数之间的关系,探讨破裂复杂度的全球及区域空间分布特征。结果表明,矩震级很大(Mw≥8.5)的事件,地震破裂复杂度更大;破裂复杂度较高的地震分布在浅层地壳(≤30km)内的概率最大,随着震源深度增加,破裂复杂度对震源深度的敏感性逐渐消失;走滑断层机制占比较大的事件,破裂复杂度较高;破裂复杂度与地震能矩比没有明确的关系;破裂复杂度的空间分布特征与区域地质构造环境相关。破裂复杂度的空间分布特征可以分为3类,第一类是板块之间简单碰撞产生的俯冲带,板块交界处的滑动速率和方向较为一致,这种情形下以比较简单的事件类型为主;第二类是多板块交界处,或者板块交界处的滑动速率和方向在整个区域内存在差异;第三类是大陆内部的强烈挤压地带。与第一类空间分布特征相比,第二类和第三类情形下破裂复杂度相对更高。地震破裂复杂度可以在一定程度上反映区域应力场的复杂性。     

Damage to the shallow Landers fault from the nearby Hector Mine earthquake

Crustal faults have long been identified as sites where localized sliding motion occurs during earthquakes, which allows for the relative motion between adjacent crustal blocks. Although there is a growing awareness that we must understand the evolution of fault systems on many timescales to relate present-day crustal stresses and fault motions to geological structures formed in the past, fault-zone damage and healing have been documented quantitatively in only a few cases. We have been monitoring the healing of damage on the shallow Johnson Valley fault after its rupture in the 1992 magnitude-7.3 Landers earthquake, and here we report that this healing was interrupted in 1999 by the magnitude-7.1 Hector Mine earthquake rupture, which occurred 20-30 km away. The Hector Mine earthquake both strongly shook and permanently strained the Johnson Valley fault, adding damage discernible as a temporary reversal of the healing process. The fault has since resumed the trend of strength recovery that it showed after the Landers earthquake. These observations lead us to speculate that fault damage caused by strong seismic waves may help to explain earthquake clustering and seismicity triggering by shaking, and may be involved in friction reduction during faulting.     

Migration of seismic scatterers associated with the 1993 Parkfield aseismic transient event

The time-varying deformation field within a fault zone, particularly at depths where earthquakes occur, is important for understanding fault behaviour and its relation to earthquake occurrence. But detection of this temporal variation has been extremely difficult, although laboratory studies have long suggested that certain structural changes, such as the properties of crustal fractures, should be seismically detectable. Here we present evidence that such structural changes are indeed observable. In particular, we find a systematic temporal variation in the seismograms of repeat microearthquakes that occurred on the Parkfield segment of the San Andreas fault over the decade 1987-97. Our analysis reveals a change of the order of 10 m in the location of scatterers which plausibly lie within the fault zone at a depth of approximately 3 km. The motion of the scatterers is coincident, in space and time, with the onset of a well documented aseismic transient (deformation event). We speculate that this structural change is the result of a stress-induced redistribution of fluids in fluid-filled fractures caused by the transient event.     

The spatial distribution pattern of landslides triggered by the 20 April 2013 Lushan earthquake of China and its implication to identification of the seismogenic fault

After the 20 April 2013 Lushan MS6.6 earthquake occurred,investigation and identification of the seismogenic fault for this event have become a focused and debatable issue.This work prepared an initial landslide inventory map related to the Lushan earthquake based on field investigations and visual interpretation of high-resolution aerial photographs and provided evidence for solving the issue aforementioned.The analysis of three landslide-density profiles perpendicular to strike direction of the probable seismogenic fault shows that many landslides occurred on the footwall of the Shuangshi–Dachuan fault(SDF),without sudden change of landslide density near the fault.Very few landslides were detected near the Dayi fault(DF)and also no change of landslide density there.While obvious sudden change of landslide density appeared about 1–2 km from the northwest to the western Shangli fault(WSF),and the landslide density on the hanging wall of the fault is obviously higher than that of on the footwall.Therefore,we infer that the seismogenic fault for the Lushan earthquake is neither the SDF nor the DF,rather probably the WSF located between these two faults,which is an evident linear trace on the earth surface.Meanwhile,the coseismic slip did not propagate upward to the ground,implying the Lushan earthquake was spawned by a blind-thrust-fault beneath the WSF.     

A great earthquake doublet and seismic stress transfer cycle in the central Kuril islands

Temporal variations of the frictional resistance on subduction-zone plate boundary faults associated with the stick-slip cycle of large interplate earthquakes are thought to modulate the stress regime and earthquake activity within the subducting oceanic plate. Here we report on two great earthquakes that occurred near the Kuril islands, which shed light on this process and demonstrate the enhanced seismic hazard accompanying triggered faulting. On 15 November 2006, an event of moment magnitude 8.3 ruptured the shallow-dipping plate boundary along which the Pacific plate descends beneath the central Kuril arc. The thrust ruptured a seismic gap that previously had uncertain seismogenic potential, although the earlier occurrence of outer-rise compressional events had suggested the presence of frictional resistance. Within minutes of this large underthrusting event, intraplate extensional earthquakes commenced in the outer rise region seaward of the Kuril trench, and on 13 January 2007, an event of moment magnitude 8.1 ruptured a normal fault extending through the upper portion of the Pacific plate, producing one of the largest recorded shallow extensional earthquakes. This energetic earthquake sequence demonstrates the stress transfer process within the subducting lithosphere, and the distinct rupture characteristics of these great earthquakes illuminate differences in seismogenic properties and seismic hazard of such interplate and intraplate faults.     

Fracture surface energy of the Punchbowl fault, San Andreas system

Fracture energy is a form of latent heat required to create an earthquake rupture surface and is related to parameters governing rupture propagation and processes of slip weakening. Fracture energy has been estimated from seismological and experimental rock deformation data, yet its magnitude, mechanisms of rupture surface formation and processes leading to slip weakening are not well defined. Here we quantify structural observations of the Punchbowl fault, a large-displacement exhumed fault in the San Andreas fault system, and show that the energy required to create the fracture surface area in the fault is about 300 times greater than seismological estimates would predict for a single large earthquake. If fracture energy is attributed entirely to the production of fracture surfaces, then all of the fracture surface area in the Punchbowl fault could have been produced by earthquake displacements totalling <1 km. But this would only account for a small fraction of the total energy budget, and therefore additional processes probably contributed to slip weakening during earthquake rupture.     

Space geodetic evidence for rapid strain rates in the New Madrid seismic zone of central USA

In the winter of 1811-1812, near the town of New Madrid in the central United States and more than 2,000 km from the nearest plate boundary, three earthquakes within three months shook the entire eastern half of the country and liquefied the ground over distances far greater than any historic earthquake in North America. The origin and modern significance of these earthquakes, however, is highly contentious. Geological evidence demonstrates that liquefaction due to strong ground shaking, similar in scale to that generated by the New Madrid earthquakes, has occurred at least three and possibly four times in the past 2,000 years (refs 4-6), consistent with recurrence statistics derived from regional seismicity. Here we show direct evidence for rapid strain rates in the area determined from a continuously operated global positioning system (GPS) network. Rates of strain are of the order of 10(-7) per year, comparable in magnitude to those across active plate boundaries, and are consistent with known active faults within the region. These results have significant implications for the definition of seismic hazard and for processes that drive intraplate seismicity.     

Some thoughts on seismotectonics of major earthquake occurrence zones in China

A major earthquake occurrence zone means a place where M≥6 events have occurred since the Holocene and similar shocks may happen again in the future. The dynamic context of the major earthquake occurrence zones in China is primarily associated with the NNE-directed push of the India plate,next with the westward subduction of the Pacific plate. The Chinese mainland is a grand mosaic structure of many crust blocks bounded by faults and sutures. When it is suffered from boundary stresses,deformation takes place along these faults or sutures while the block interiors remain relatively stable or intact. Since the Quaternary,for example,left slip on the Xianshuihe-Xiaojiang fault zone in southwestern China has produced a number of fault-depression basins in extensional areas during periods Q1 and Q2. In the Q3,the change of stress orientation and enhancement of tectonic movement made faults of varied trends link each other,and continued to be active till present day,producing active fault zones in this region. Usually major earthquakes occur at some special locations on these active fault zones. During these events,in the epicenter areas experience intensive deformation characterized by large-amplitude rise and fall of neighboring sections,generation of horst-graben systems and dammed rivers. The studies on palaeoearthquakes suggest that major shocks of close magnitudes often repeated for several times at a same place. By comparison of the Chi-Chi,Taiwan event in 1999 and Yuza,Yunnan event in 1955,including contours of accelerations and intensities,destruction of buildings,and in contrast to the Xigeda formation in southwestern China,a sandwich model is established to account for the mechanism of deformation caused by major earthquakes. This model consists of three layers,i. e. the two walls of a fault and the ruptured zone intercalated between them. This ruptured zone is just the loci where stress is built up and released,and serves as a channel for seismic waves.     

Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin

On 12 May 2008, the devastating magnitude 7.9 (Wenchuan) earthquake struck the eastern edge of the Tibetan plateau, collapsing buildings and killing thousands in major cities aligned along the western Sichuan basin in China. After such a large-magnitude earthquake, rearrangement of stresses in the crust commonly leads to subsequent damaging earthquakes. The mainshock of the 12 May earthquake ruptured with as much as 9 m of slip along the boundary between the Longmen Shan and Sichuan basin, and demonstrated the complex strike-slip and thrust motion that characterizes the region. The Sichuan basin and surroundings are also crossed by other active strike-slip and thrust faults. Here we present calculations of the coseismic stress changes that resulted from the 12 May event using models of those faults, and show that many indicate significant stress increases. Rapid mapping of such stress changes can help to locate fault sections with relatively higher odds of producing large aftershocks.     

Variation of b value with hypocentral depth in Beijing area： Implications for earthquake nucleation

We relocated 2098 earthquakes that occurred in Beijing area between 1980 and 2000 using a double-difference (DD) earthquake location algorithm, and obtained the high-precision relative locations of 1825 events, b values versus depth were investigated with the relocated hypocenters. The results show that the b values decrease with the increasing hypocentral depth systematically. A dramatic variation in b is observed around the depth of 8 km. It indicates that there are more smaller earthquakes at shallow depth (0-8 km), while more larger earthquakes at greater depth (8--25 km). The physical mechanism behind this phenomenon can be explained by the variations in material heterogeneity and lithostatic stress condition. Large earthquakes are more likely to nucleate at greater depth with more homogenous material and higher lithostatic stress. On the basis of the results, we suggest that future strong earthquakes in Beijing area tend to occur below the depth of 8 km.     

The Lushan M S7.0 earthquake and activity of the southern segment of the Longmenshan fault zone

Following the Lushan M _S7.0 earthquake on 20 April 2013, a topic of much concern is whether events of M _S7 or greater could occur again on the southern segment of the Longmenshan fault zone. In providing evidence to answer this question, this work analyzes the tectonic relationship between the Lushan event and the 2008 Wenchuan earthquake and the rupture history of the southern segment of the Longmenshan fault zone, through field investigations of active tectonics and paleoearthquake research, and our preliminary conclusions are as follows. The activity of the southern segment of the Longmenshan fault zone is much different to that of its central section, and the late Quaternary activity has propagated forward to the basin in the east. The seismogenic structure of the 2008 Wenchuan earthquake is the central-fore-range fault system, whereas that of the 2013 Lushan event is attributed to the fore-range-range-front fault system, rather than the central fault. The southern segment of the Longmenshan fault zone becomes wider towards the south with an increasing number of secondary faults, of which the individual faults exhibit much weaker surface activity. Therefore, this section is not as capable of generating a major earthquake as is the central segment. It is most likely that the 2013 earthquake fills the seismic gap around Lushan on the southern segment of the Longmenshan fault zone.     

近断层长周期地震动前向性效应对峰值影响研究

筛选了断层距不大于15km的具有前向性效应的长周期近断层地震记录,分类汇总建立了统计分析数据样本库,对长周期地震动的竖向分量和水平向分量峰值比进行了分析,研究了震级、脉冲周期、不同分量对峰值加速度比的影响。在此基础上,采用前期建立的回归分析模型,对长周期地震动竖向分量和水平向两分量的峰值加速度进行了回归分析,拟合给出了峰值参数与震级和断层距的衰减关系式,通过拟合曲线与数据样本的峰值实测值对比,说明拟合关系式在一定的范围内可以用于近断层长周期地震峰值的预测。     

利用接收函数反演龙门山断裂带及邻区深部结构

利用接收函数的方法通过接收震中距30°～90°、震级在5.5以上的远震事件反演龙门山断裂带及其邻区的深部结构,探索汶川地震形成原因。结果表明,扬子地台西缘的莫霍面向西侧倾斜缓降;处在龙门山推覆体范围之内的都江堰、汶川一带莫霍面起伏变化不大,在跨过龙门山中央断裂带后开始下降,向北降至黑水县附近后平缓上升。结合2005年10月至2007年4月远震P波波形资料接收函数反演结果:①2条被动源剖面均显示莫霍面在龙门山推覆体中央位置深度约43km的地方出现不同程度的陡降,说明该断裂带是地壳厚度的陡变带,为扬子地台和松潘甘孜地台的构造边界。②莫霍面深度向南陡降至最深约68km处后平缓上升,向北陡降至最深约58km处后平缓上升。表明松潘-甘孜地块东缘地壳厚度呈南深北浅、东深西浅分布。     

The M S7.1 Yushu earthquake surface rupture and large historical earthquakes on the Garzê-Yushu Fault

As revealed by field investigations, the co-seismic surface rupture zone of the 2010 M_S7.1 Yushu earthquake, Qinghai is a char-acteristic sinistral strike-slip feature consisting of three distinct sinistral primary ruptures, with an overall strike of 310°–320° and a total length of 31 km. In addition, an approximately 2-km-long en-echelon tensile fissure zone was found east of Longbao Town; if this site is taken as the north end of the rupture zone, then the rupture had a total length of ~51 km. The surface rupture zone is composed of a series of fissures arranged in an en-echelon or alternating relationship between compressive bulges and tensile fissures, with a measured maximum horizontal displacement of 1.8 m. The surface rupture zone extends along the mapped Garzê-Yushu Fault, which implicates it as the seismogenic fault for this earthquake. Historically, a few earthquakes with a magnitude of about 7 have occurred along the fault, and additionally traces of paleoearthquakes are evident that characterize the short-period recurrence interval of large earthquakes here. Similar to the seismogenic process of the 2008 Wenchuan earthquake, the Yushu earthquake is also due to the stress accumulation and release on the block boundaries resulting from the eastward expansion of Qinghai-Tibet Plateau. However, in contrast with the Wenchuan earthquake, the Yushu earthquake had a sinistral strike-slip mechanism resulting from the uneven eastward extrusion of the Baryan Har and Sichuan-Yunnan fault blocks.