Decay of aftershock density with distance indicates triggering by dynamic stress

The majority of earthquakes are aftershocks, yet aftershock physics is not well understood. Many studies suggest that static stress changes trigger aftershocks, but recent work suggests that shaking (dynamic stresses) may also play a role. Here we measure the decay of aftershocks as a function of distance from magnitude 2-6 mainshocks in order to clarify the aftershock triggering process. We find that for short times after the mainshock, when low background seismicity rates allow for good aftershock detection, the decay is well fitted by a single inverse power law over distances of 0.2-50 km. The consistency of the trend indicates that the same triggering mechanism is working over the entire range. As static stress changes at the more distant aftershocks are negligible, this suggests that dynamic stresses may be triggering all of these aftershocks. We infer that the observed aftershock density is consistent with the probability of triggering aftershocks being nearly proportional to seismic wave amplitude. The data are not fitted well by models that combine static stress change with the evolution of frictionally locked faults.     

Interaction between adjacent left-lateral strike-slip faults and thrust faults: the 1976 Songpan earthquake sequence

Based on the published focal mechanisms we have built the fault model of the main shocks of the 1976 Songpan earthquake sequence and calculated the coseismic Coulomb stress changes in the region. The results show that most of the aftershocks had occurred in the region where the Coulomb stresses had been increased, indicating a triggering relationship between the main shocks and the aftershocks. We also show that the first main shock （Ms = 7.2）, which is a left-lateral slip event, had increased the Coulomb stresses by 5×10^5 Pa at the second main shock （a thrust event with Ms = 6.7）. Therefore, we conclude that the first main shock had triggered the second main shock. The third main shock is also a left-lateral event, however, the triggering relationship between the third main shock and the previous two events is less obvious. General model calculations show that there is a good triggering relationship between adjacent left-lateral slip fault and thrust fault, but triggering between parallel slip faults is rather weak.     

Foreshock sequences and short-term earthquake predictability on East Pacific Rise transform faults

East Pacific Rise transform faults are characterized by high slip rates (more than ten centimetres a year), predominantly aseismic slip and maximum earthquake magnitudes of about 6.5. Using recordings from a hydroacoustic array deployed by the National Oceanic and Atmospheric Administration, we show here that East Pacific Rise transform faults also have a low number of aftershocks and high foreshock rates compared to continental strike-slip faults. The high ratio of foreshocks to aftershocks implies that such transform-fault seismicity cannot be explained by seismic triggering models in which there is no fundamental distinction between foreshocks, mainshocks and aftershocks. The foreshock sequences on East Pacific Rise transform faults can be used to predict (retrospectively) earthquakes of magnitude 5.4 or greater, in narrow spatial and temporal windows and with a high probability gain. The predictability of such transform earthquakes is consistent with a model in which slow slip transients trigger earthquakes, enrich their low-frequency radiation and accommodate much of the aseismic plate motion.     

日本Mw9.0级地震的库仑应力变化及对余震的触发

采用已公布的1个前震和3个主震震源可变滑动模型,计算了日本Mw7.3级前震、M_w9.0级主震产生的同震库仑应力变化,并采用哈佛Global CMT和日本F-net余震目录,在不同有效摩擦系数取值下计算了主震对余震的触发效应,分析了应力空间分布与余震活动空间分布的相关性。研究表明,日本地震前震对主震有触发作用,前震在主震处产生的正应力变化为0.404 bar,剪应力变化0.282 bar,同震库仑应力变化为0.517 bar,应力变化超过了应力触发阈值0.1 bar,达5倍之多,属于典型的地震应力触发,而非诱发事件。主震对后续余震有触发效应,最小触发率为57%,最大达到75%。余震多分布在主震产生的同震库仑应力红色加载区。同时,分析了可变滑动模型、有效摩擦系数、余震目录、震级及节面选取等对计算结果的影响。可变滑动模型、有效摩擦系数变化对计算结果影响不大,余震目录、震级及节面选取对结果有一定影响,尤其是当余震目录足够多、震级足够大时,计算结果更加可靠。为了深入讨论节面选取对计算结果的影响,对余震两个节面上库仑应力、剪切应力及正应力变化差值进行了频次统计。分析显示节面选取对正应力有一定影响,而对库仑应力及剪应力影响不大。     

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.     

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.     

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.     

Post-earthquake ground movements correlated to pore-pressure transients

Large earthquakes alter the stress in the surrounding crust, leading to triggered earthquakes and aftershocks. A number of time-dependent processes, including afterslip, pore-fluid flow and viscous relaxation of the lower crust and upper mantle, further modify the stress and pore pressure near the fault, and hence the tendency for triggered earthquakes. It has proved difficult, however, to distinguish between these processes on the basis of direct field observations, despite considerable effort. Here we present a unique combination of measurements consisting of satellite radar interferograms and water-level changes in geothermal wells following two magnitude-6.5 earthquakes in the south Iceland seismic zone. The deformation recorded in the interferograms cannot be explained by either afterslip or visco-elastic relaxation, but is consistent with rebound of a porous elastic material in the first 1-2 months following the earthquakes. This interpretation is confirmed by direct measurements which show rapid (1-2-month) recovery of the earthquake-induced water-level changes. In contrast, the duration of the aftershock sequence is projected to be approximately 3.5 years, suggesting that pore-fluid flow does not control aftershock duration. But because the surface strains are dominated by pore-pressure changes in the shallow crust, we cannot rule out a longer pore-pressure transient at the depth of the aftershocks. The aftershock duration is consistent with models of seismicity rate variations based on rate- and state-dependent friction laws.     

Effects of acoustic waves on stick-slip in granular media and implications for earthquakes

It remains unknown how the small strains induced by seismic waves can trigger earthquakes at large distances, in some cases thousands of kilometres from the triggering earthquake, with failure often occurring long after the waves have passed. Earthquake nucleation is usually observed to take place at depths of 10-20 km, and so static overburden should be large enough to inhibit triggering by seismic-wave stress perturbations. To understand the physics of dynamic triggering better, as well as the influence of dynamic stressing on earthquake recurrence, we have conducted laboratory studies of stick-slip in granular media with and without applied acoustic vibration. Glass beads were used to simulate granular fault zone material, sheared under constant normal stress, and subject to transient or continuous perturbation by acoustic waves. Here we show that small-magnitude failure events, corresponding to triggered aftershocks, occur when applied sound-wave amplitudes exceed several microstrain. These events are frequently delayed or occur as part of a cascade of small events. Vibrations also cause large slip events to be disrupted in time relative to those without wave perturbation. The effects are observed for many large-event cycles after vibrations cease, indicating a strain memory in the granular material. Dynamic stressing of tectonic faults may play a similar role in determining the complexity of earthquake recurrence.     

Great Himalayan earthquakes and the Tibetan plateau

It has been assumed that Himalayan earthquakes are driven by the release of compressional strain accumulating close to the Greater Himalaya. However, elastic models of the Indo-Asian collision using recently imaged subsurface interface geometries suggest that a substantial fraction of the southernmost 500 kilometres of the Tibetan plateau participates in driving great ruptures. We show here that this Tibetan reservoir of elastic strain energy is drained in proportion to Himalayan rupture length, and that the consequent growth of slip and magnitude with rupture area, when compared to data from recent earthquakes, can be used to infer a approximately 500-year renewal time for these events. The elastic models also illuminate two puzzling features of plate boundary seismicity: how great earthquakes can re-rupture regions that have already ruptured in recent smaller earthquakes and how mega-earthquakes with greater than 20 metres slip may occur at millennia-long intervals, driven by residual strain following many centuries of smaller earthquakes.     

Triggering of earthquake aftershocks by dynamic stresses

It is thought that small 'static' stress changes due to permanent fault displacement can alter the likelihood of, or trigger, earthquakes on nearby faults. Many studies of triggering in the near-field, particularly of aftershocks, rely on these static changes as the triggering agent and consider them only in terms of equivalent changes in the applied load on the fault. Here we report a comparison of the aftershock pattern of the moment magnitude Mw = 7.3 Landers earthquake, not only with static stress changes but also with transient, oscillatory stress changes transmitted as seismic waves (that is, 'dynamic' stresses). Dynamic stresses do not permanently change the applied load and thus can trigger earthquakes only by altering the mechanical state or properties of the fault zone. These dynamically weakened faults may fail after the seismic waves have passed by, and might even cause earthquakes that would not otherwise have occurred. We find similar asymmetries in the aftershock and dynamic stress patterns, the latter being due to rupture propagation, whereas the static stress changes lack this asymmetry. Previous studies have shown that dynamic stresses can promote failure at remote distances, but here we show that they can also do so nearby.     

Lushan M S7.0 earthquake: A special earthquake occurs on curved fault

Relocation result shows that the aftershocks of the Lushan M _S7.0 earthquake spatially distribute in a shape like “half bowl”, indicating that the rupture structure of the mainshock is a highly curved surface. Kinematic analysis reveals that a laterally varied dislocation pattern occurs on this curved fault even though a single relative horizontal movement controls slip on this fault. Reverse slip prevails on curved fault. However, significant normal slip is predicted near the edge of north flank. Moreover, the north flank features left-lateral slip while the south flank contrarily features right-lateral slip. The relative scope of aftershock distribution implies inadequate breaking of the curved fault during the mainshock, calling for the attention to potential earthquake risk on the neighboring portions of the coseismic rupture due to significant increase of the coseismic Coulomb stress. Coseismic stress modeling also reveals that it is unnecessary for the stress on ruptured part to be unloaded following the earthquakes on the curved fault. The coseismic stress loading on ruptured elements unveils the specialty of faulting for the Lushan earthquake and we conclude that this specialty is due to the highly curved fault geometry.     

Preseismic velocity changes observed from active source monitoring at the Parkfield SAFOD drill site

Measuring stress changes within seismically active fault zones has been a long-sought goal of seismology. One approach is to exploit the stress dependence of seismic wave velocity, and we have investigated this in an active source cross-well experiment at the San Andreas Fault Observatory at Depth (SAFOD) drill site. Here we show that stress changes are indeed measurable using this technique. Over a two-month period, we observed an excellent anti-correlation between changes in the time required for a shear wave to travel through the rock along a fixed pathway (a few microseconds) and variations in barometric pressure. We also observed two large excursions in the travel-time data that are coincident with two earthquakes that are among those predicted to produce the largest coseismic stress changes at SAFOD. The two excursions started approximately 10 and 2 hours before the events, respectively, suggesting that they may be related to pre-rupture stress induced changes in crack properties, as observed in early laboratory studies.     

地震发生与日月运行之关联

 基于第23、24太阳活动周发生的全球M≥7.8级以上大震数据、汶川余震和智利余震数据,发现：地震的发生与日月位置存在着关联性;地震的发生与太阳风磁场到达地球的优势聚集方向存在着关联性。说明地震的发生与日月运行有关联。所用地震事件有两类：一类为余震——2008-05-12 M8.0汶川地震M≥5.0强余震与2010-02-27 M8.8智利地震M≥5.5强余震;另一类为1997—2010年全球M≥7.8大地震。虽然参与分析的地震事件数不算很多,但从两类各自不同的代表性以及所使用地震震级与地域的跨度,仍可望其结论具有普适性。文中提供了众多地震发生与日月运行关联的证据,其随机发生概率多为10^-4—10^-6,小概率为不发生事件,但竟然发生了,表明它们之间必定有其非随机的统一成因机制控制。这是本文作者对地震发生的宇宙环境——“天外来客”特别给予重视的立论所在。     

论构造应力不是油气初次运移的主要动力

油气初次运移是一个尚未完全解决的学术难题。石油地质学界认为,构造应力是油气初次运移的主要动力之一,笔者对该观点进行了质疑。根据有效应力理论和弹性力学理论,分析了构造应力对孔隙压力的影响,建立了孔隙压力增量的计算模型以及孔隙压力释放驱动油气运移尺度的计算模型。结果表明,相同构造应力条件下,烃源岩的孔隙压力增量主要受杨氏模量影响,杨氏模量越大,孔隙压力增量越小;构造应力导致的孔隙压力增量一般不会超过60 MPa;孔隙压力释放只需要排出少量地层流体即可,驱动油气运移的距离很短,相对于烃源岩厚度几乎可以忽略不计;构造应力对油气初次运移的影响十分微弱,不可能是初次运移的主要动力。     

Decay of aftershock density with distance does not indicate triggering by dynamic stress

Resolving whether static or dynamic stress triggers most aftershocks and subsequent mainshocks is essential to understand earthquake interaction and to forecast seismic hazard. Felzer and Brodsky examined the distance distribution of earthquakes occurring in the first five minutes after 2?≤?M?相似文献   

四川芦山Ms 7.0级地震导致周边断层的应力变化

 2013年4月20日,四川雅安市芦山县发生M_s 7.0级强烈地震,造成重大的人员伤亡和经济损失。地震造成的区域库仑应力变化、对周围断层的影响及后续地震带发展趋势是应该关注的问题。利用USGS震源机制解,根据地震静态触发原理,基于弹性位错理论和分层地壳模型,计算得出芦山地震引起同震库仑应力变化从断层的1.0MPa量级减小到200km外的0.1kPa,研究认为在地震之后大部分区域的应力得到释放,鲜水河断裂道孚—康定段和玉龙希断裂南段危险性增加。     

Short-lived orogenic cycles and the eclogitization of cold crust by spasmodic hot fluids

Collision tectonics and the associated transformation of continental crust to high-pressure rocks (eclogites) are generally well-understood processes, but important contradictions remain between tectonothermal models and petrological-isotopic data obtained from such rocks. Here we use 40Ar-39Ar data coupled with a thermal model to constrain the time-integrated duration of an orogenic cycle (the burial and exhumation of a particular segment of the crust) to be less than 13 Myr. We also determine the total duration of associated metamorphic events to be approximately 20 kyr, and of individual heat pulses experienced by the rocks to be as short as 10 years. Such short timescales are indicative of rapid tectonic processes associated with catastrophic deformation events (earthquakes). Such events triggered transient heat advection by hot fluid along deformation (shear) zones, which cut relatively cool and dry subducted crust. In contrast to current thermal models that assume thermal equilibrium and invoke high ambient temperatures in the thickened crust, our non-steady-state cold-crust model satisfactorily explains several otherwise contradictory geological observations.     

Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake

Most large earthquakes occur along an oceanic trench, where an oceanic plate subducts beneath a continental plate. Massive earthquakes with a moment magnitude, M(w), of nine have been known to occur in only a few areas, including Chile, Alaska, Kamchatka and Sumatra. No historical records exist of a M(w) = 9 earthquake along the Japan trench, where the Pacific plate subducts beneath the Okhotsk plate, with the possible exception of the ad 869 Jogan earthquake, the magnitude of which has not been well constrained. However, the strain accumulation rate estimated there from recent geodetic observations is much higher than the average strain rate released in previous interplate earthquakes. This finding raises the question of how such areas release the accumulated strain. A megathrust earthquake with M(w) = 9.0 (hereafter referred to as the Tohoku-Oki earthquake) occurred on 11 March 2011, rupturing the plate boundary off the Pacific coast of northeastern Japan. Here we report the distributions of the coseismic slip and postseismic slip as determined from ground displacement detected using a network based on the Global Positioning System. The coseismic slip area extends approximately 400?km along the Japan trench, matching the area of the pre-seismic locked zone. The afterslip has begun to overlap the coseismic slip area and extends into the surrounding region. In particular, the afterslip area reached a depth of approximately 100?km, with M(w) = 8.3, on 25 March 2011. Because the Tohoku-Oki earthquake released the strain accumulated for several hundred years, the paradox of the strain budget imbalance may be partly resolved. This earthquake reminds us of the potential for M(w)?≈?9 earthquakes to occur along other trench systems, even if no past evidence of such events exists. Therefore, it is imperative that strain accumulation be monitored using a space geodetic technique to assess earthquake potential.     

Nonlinear dynamics, granular media and dynamic earthquake triggering

The 1992 magnitude 7.3 Landers earthquake triggered an exceptional number of additional earthquakes within California and as far north as Yellowstone and Montana. Since this observation, other large earthquakes have been shown to induce dynamic triggering at remote distances--for example, after the 1999 magnitude 7.1 Hector Mine and the 2002 magnitude 7.9 Denali earthquakes--and in the near-field as aftershocks. The physical origin of dynamic triggering, however, remains one of the least understood aspects of earthquake nucleation. The dynamic strain amplitudes from a large earthquake are exceedingly small once the waves have propagated more than several fault radii. For example, a strain wave amplitude of 10(-6) and wavelength 1 m corresponds to a displacement amplitude of about 10(-7) m. Here we show that the dynamic, elastic-nonlinear behaviour of fault gouge perturbed by a seismic wave may trigger earthquakes, even with such small strains. We base our hypothesis on recent laboratory dynamic experiments conducted in granular media, a fault gouge surrogate. From these we infer that, if the fault is weak, seismic waves cause the fault core modulus to decrease abruptly and weaken further. If the fault is already near failure, this process could therefore induce fault slip.