Slip zone and energetics of a large earthquake from the Taiwan Chelungpu-fault Drilling Project

Determining the seismic fracture energy during an earthquake and understanding the associated creation and development of a fault zone requires a combination of both seismological and geological field data. The actual thickness of the zone that slips during the rupture of a large earthquake is not known and is a key seismological parameter in understanding energy dissipation, rupture processes and seismic efficiency. The 1999 magnitude-7.7 earthquake in Chi-Chi, Taiwan, produced large slip (8 to 10 metres) at or near the surface, which is accessible to borehole drilling and provides a rare opportunity to sample a fault that had large slip in a recent earthquake. Here we present the retrieved cores from the Taiwan Chelungpu-fault Drilling Project and identify the main slip zone associated with the Chi-Chi earthquake. The surface fracture energy estimated from grain sizes in the gouge zone of the fault sample was directly compared to the seismic fracture energy determined from near-field seismic data. From the comparison, the contribution of gouge surface energy to the earthquake breakdown work is quantified to be 6 per cent.     

Earthquake slip weakening and asperities explained by thermal pressurization

An earthquake occurs when a fault weakens during the early portion of its slip at a faster rate than the release of tectonic stress driving the fault motion. This slip weakening occurs over a critical distance, D(c). Understanding the controls on D(c) in nature is severely limited, however, because the physical mechanism of weakening is unconstrained. Conventional friction experiments, typically conducted at slow slip rates and small displacements, have obtained D(c) values that are orders of magnitude lower than values estimated from modelling seismological data for natural earthquakes. Here we present data on fluid transport properties of slip zone rocks and on the slip zone width in the centre of the Median Tectonic Line fault zone, Japan. We show that the discrepancy between laboratory and seismological results can be resolved if thermal pressurization of the pore fluid is the slip-weakening mechanism. Our analysis indicates that a planar fault segment with an impermeable and narrow slip zone will become very unstable during slip and is likely to be the site of a seismic asperity.     

Seismology: speed and size of the Sumatra earthquake

Our seismological results reveal that Indonesia's devastating Sumatra-Andaman earthquake on 26 December 2004 was 2.5 times larger than initial reports suggested--second only to the 1960 Chilean earthquake in recorded magnitude. They indicate that it slowly released its energy by slip along a 1,200-km fault, generating a long rupture that contributed to the subsequent tsunami. Now that the entire rupture zone has slipped, the strain accumulated from the subduction of the Indian plate beneath the Burma microplate has been released, and there is no immediate danger of a similar tsunami being generated on this part of the plate boundary, although large earthquakes on segments to the south still present a threat.     

Fault lubrication during earthquakes

The determination of rock friction at seismic slip rates (about 1?m?s(-1)) is of paramount importance in earthquake mechanics, as fault friction controls the stress drop, the mechanical work and the frictional heat generated during slip. Given the difficulty in determining friction by seismological methods, elucidating constraints are derived from experimental studies. Here we review a large set of published and unpublished experiments (～300) performed in rotary shear apparatus at slip rates of 0.1-2.6?m?s(-1). The experiments indicate a significant decrease in friction (of up to one order of magnitude), which we term fault lubrication, both for cohesive (silicate-built, quartz-built and carbonate-built) rocks and non-cohesive rocks (clay-rich, anhydrite, gypsum and dolomite gouges) typical of crustal seismogenic sources. The available mechanical work and the associated temperature rise in the slipping zone trigger a number of physicochemical processes (gelification, decarbonation and dehydration reactions, melting and so on) whose products are responsible for fault lubrication. The similarity between (1) experimental and natural fault products and (2) mechanical work measures resulting from these laboratory experiments and seismological estimates suggests that it is reasonable to extrapolate experimental data to conditions typical of earthquake nucleation depths (7-15?km). It seems that faults are lubricated during earthquakes, irrespective of the fault rock composition and of the specific weakening mechanism involved.     

联合slipBERI与D-InSAR技术的青海省玛多县地震形变场提取与模拟

针对2021年5月22日青海玛多县Ms7.4地震震区形变信息、形变特征以及滑动断裂特性的提取与模拟等问题,获取了玛多县地震区的Sentinel-1A影像,采用双轨差分干涉法,并优化各项参数,提取出同震形变场,利用slipBERI(slip from Bayesian Regularized Inversion)方法对断层的几何参数及形变场进行反演和模拟。结果表明：玛多地震同震形变场的形状近似于一个椭圆,断层整体呈西北-东南走向,其上部为沉降区,下部为隆升区,最大LOS(line of sight,视线向)形变分别为0.65m和0.81m。地震形变场的运动主要以东西方向的水平运动为主,并伴有明显的左旋走滑,断层上下方相对视线向运动可达1.50m,表明此次地震的地表破裂有明显的错位移动。通过分析形变信息和地表破裂特征,可以判断该破裂带位于巴颜喀拉块体,为昆仑山口-江口断裂,反演结果与观测结果相符,这表明观测结果较可靠。     

Particle size and energetics of gouge from earthquake rupture zones

Grain size reduction and gouge formation are found to be ubiquitous in brittle faults at all scales, and most slip along mature faults is observed to have been localized within gouge zones. This fine-grain gouge is thought to control earthquake instability, and thus understanding its properties is central to an understanding of the earthquake process. Here we show that gouge from the San Andreas fault, California, with approximately 160 km slip, and the rupture zone of a recent earthquake in a South African mine with only approximately 0.4 m slip, display similar characteristics, in that ultrafine grains approach the nanometre scale, gouge surface areas approach 80 m2 g(-1), and grain size distribution is non-fractal. These observations challenge the common perception that gouge texture is fractal and that gouge surface energy is a negligible contributor to the earthquake energy budget. We propose that the observed fine-grain gouge is not related to quasi-static cumulative slip, but is instead formed by dynamic rock pulverization during the propagation of a single earthquake.     

The deterministic nature of earthquake rupture

Understanding the earthquake rupture process is central to our understanding of fault systems and earthquake hazards. Multiple hypotheses concerning the nature of fault rupture have been proposed but no unifying theory has emerged. The conceptual hypothesis most commonly cited is the cascade model for fault rupture. In the cascade model, slip initiates on a small fault patch and continues to rupture further across a fault plane as long as the conditions are favourable. Two fundamental implications of this domino-like theory are that small earthquakes begin in the same manner as large earthquakes and that the rupture process is not deterministic--that is, the size of the earthquake cannot be determined until the cessation of rupture. Here we show that the frequency content of radiated seismic energy within the first few seconds of rupture scales with the final magnitude of the event. We infer that the magnitude of an earthquake can therefore be estimated before the rupture is complete. This finding implies that the rupture process is to some degree deterministic and has implications for the physics of the rupture process.     

Insight into the 2004 Sumatra-Andaman earthquake from GPS measurements in southeast Asia

Data collected at approximately 60 Global Positioning System (GPS) sites in southeast Asia show the crustal deformation caused by the 26 December 2004 Sumatra-Andaman earthquake at an unprecedented large scale. Small but significant co-seismic jumps are clearly detected more than 3,000 km from the earthquake epicentre. The nearest sites, still more than 400 km away, show displacements of 10 cm or more. Here we show that the rupture plane for this earthquake must have been at least 1,000 km long and that non-homogeneous slip is required to fit the large displacement gradients revealed by the GPS measurements. Our kinematic analysis of the GPS recordings indicates that the centroid of released deformation is located at least 200 km north of the seismological epicentre. It also provides evidence that the rupture propagated northward sufficiently fast for stations in northern Thailand to have reached their final positions less than 10 min after the earthquake, hence ruling out the hypothesis of a silent slow aseismic rupture.     

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.     

Earthquakes as beacons of stress change

Aftershocks occurring on faults in the far-field of a large earthquake rupture can generally be accounted for by changes in static stress on these faults caused by the rupture. This implies that faults interact, and that the timing of an earthquake can be affected by previous nearby ruptures. Here we explore the potential of small earthquakes to act as 'beacons' for the mechanical state of the crust. We investigate the static-stress changes resulting from the 1992 Landers earthquake in southern California which occurred in an area of high seismic activity stemming from many faults. We first gauge the response of the regional seismicity to the Landers event with a new technique, and then apply the same method to the inverse problem of determining the slip distribution on the main rupture from the seismicity. Assuming justifiable parameters, we derive credible matches to slip profiles obtained directly from the Landers mainshock. Our results provide a way to monitor mechanical conditions in the upper crust, and to investigate processes leading to fault failure.     

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

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

形成新生断裂的临界角度及其与强烈地震发生的关系

选取具有薄弱面或破裂面的细粒花岗岩块状露头标本,在实验室模拟高温高压环境条件进行测试,研究岩石产生新破裂面的受力临界角度。在温度100℃、围压100 MPa和温度200℃、围压200 MPa的环境条件下,按照与主应力夹角的不同度数变化分别进行实验。实验结果表明,当夹角α处于0°~30°时,岩石沿原破裂面稳滑;夹角α处于30°~60°时,以粘滑为主;夹角α处于60°~90°时,以新破裂面的产生为主;可能产生新破裂面的临界角度范围为50°~65°,当夹角α大于65°时,新破裂产生。对岩石出现的稳滑-粘滑-新破裂现象,通过分析切面所受应力(σ)与应变率(ε)的关系进行了解释。根据莫尔圆理论认为,复杂断裂背景下,岩石的受力角度对于新破裂面的产生将有着决定性作用,而作用力状态一定时,新破裂产生的临界角度60°是可能导致强烈地震发生的危险值。     

'Melt welt' mechanism of extreme weakening of gabbro at seismic slip rates

Laboratory studies of frictional properties of rocks at slip velocities approaching the seismic range (～0.1-1?m?s(-1)), and at moderate normal stresses (1-10?MPa), have revealed a complex evolution of the dynamic shear strength, with at least two phases of weakening separated by strengthening at the onset of wholesale melting. The second post-melting weakening phase is governed by viscous properties of the melt layer and is reasonably well understood. The initial phase of extreme weakening, however, remains a subject of much debate. Here we show that the initial weakening of gabbro is associated with the formation of hotspots and macroscopic streaks of melt ('melt welts'), which partially unload the rest of the slip interface. Melt welts begin to form when the average rate of frictional heating exceeds 0.1-0.4?MW?m(-2), while the average temperature of the shear zone is well below the solidus (250-450?°C). Similar heterogeneities in stress and temperature are likely to occur on natural fault surfaces during rapid slip, and to be important for earthquake rupture dynamics.     

Earthquake rupture dynamics frozen in exhumed ancient faults

Most of our knowledge about co-seismic rupture propagation is derived from inversion and interpretation of strong-ground-motion seismograms, laboratory experiments on rock and rock-analogue material, or inferred from theoretical and numerical elastodynamic models. However, additional information on dynamic rupture processes can be provided by direct observation of faults exhumed at the Earth's surface. Pseudotachylytes (solidified friction-induced melts) are the most certain fault-rock indicator of seismicity on ancient faults. Here we show how the asymmetry in distribution and the orientation of pseudotachylyte-filled secondary fractures around an exhumed fault can be used to reconstruct the earthquake rupture directivity, rupture velocity and fracture energy, by comparison with the theoretical dynamic stress field computed around propagating fractures. In particular, the studied natural network of pseudotachylytes is consistent with a dominant propagation direction during repeated seismic events and subsonic rupture propagation close to the Rayleigh wave velocity.     

Earthquake slip on oceanic transform faults

Oceanic transform faults are one of the main types of plate boundary, but the manner in which they slip remains poorly understood. Early studies suggested that relatively slow earthquake rupture might be common; moreover, it has been reported that very slow slip precedes some oceanic transform earthquakes, including the 1994 Romanche earthquake. The presence of such detectable precursors would have obvious implications for earthquake prediction. Here we model broadband seismograms of body waves to obtain well-resolved depths and rupture mechanisms for 14 earthquakes on the Romanche and Chain transform faults in the equatorial Atlantic Ocean. We found that earthquakes on the longer Romanche transform are systematically deeper than those on the neighbouring Chain transform. These depths indicate that the maximum depth of brittle failure is at a temperature of approximately 600 degrees C in oceanic lithosphere. We find that the body waves from the Romanche 1994 earthquake can be well modelled with relatively deep slip on a single fault, and we use the mechanism and depth of this earthquake to recalculate its source spectrum. The previously reported slow precursor can be explained as an artefact of uncertainties in the assumed model parameters.     

Plateau 'pop-up' in the great 1897 Assam earthquake

The great Assam earthquake of 12 June 1897 reduced to rubble all masonry buildings within a region of northeastern India roughly the size of England, and was felt over an area exceeding that of the great 1755 Lisbon earthquake. Hitherto it was believed that rupture occurred on a north-dipping Himalayan thrust fault propagating south of Bhutan. But here we show that the northern edge of the Shillong plateau rose violently by at least 11 m during the Assam earthquake, and that this was due to the rupture of a buried reverse fault approximately 110 km in length and dipping steeply away from the Himalaya. The stress drop implied by the rupture geometry and the prodigious fault slip of 18 +/- 7 m explains epicentral accelerations observed to exceed 1g vertically and surface velocities exceeding 3 m s-1 (ref. 1). This quantitative observation of active deformation of a 'pop-up' structure confirms that faults bounding such structures can penetrate the whole crust. Plateau uplift in the past 2-5 million years has caused the Indian plate to contract locally by 4 +/- 2 mm yr-1, reducing seismic risk in Bhutan but increasing the risk in northern Bangladesh.     

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.     

Resonant slow fault slip in subduction zones forced by climatic load stress

Global Positioning System (GPS) measurements at subduction plate boundaries often record fault movements similar to earthquakes but much slower, occurring over timescales of approximately 1 week to approximately 1 year. These 'slow slip events' have been observed in Japan, Cascadia, Mexico, Alaska and New Zealand. The phenomenon is poorly understood, but several observations hint at the processes underlying slow slip. Although slip itself is silent, seismic instruments often record coincident low-amplitude tremor in a narrow (1-5 cycles per second) frequency range. Also, modelling of GPS data and estimates of tremor location indicate that slip focuses near the transition from unstable ('stick-slip') to stable friction at the deep limit of the earthquake-producing seismogenic zone. Perhaps most intriguingly, slow slip is periodic at several locations, with recurrence varying from 6 to 18 months depending on which subduction zone (or even segment) is examined. Here I show that such periodic slow fault slip may be a resonant response to climate-driven stress perturbations. Fault slip resonance helps to explain why slip events are periodic, why periods differ from place to place, and why slip focuses near the base of the seismogenic zone. Resonant slip should initiate within the rupture zone of future great earthquakes, suggesting that slow slip may illuminate fault properties that control earthquake slip.     

Slip model for the 2011 M_w 9.0 Sendai （Japan） earthquake and its M_w 7.9 aftershock derived from GPS data

Using GPS-measured coseismic and post-seismic displacements for the 8 h following the M w 9.0 Sendai earthquake of March 11, 2011, coseismic and post-seismic fault slip models were developed based on a layered crustal model. The geodetic moment magnitude of the main shock was measured as approximately M w 8.98. The slip exhibits clear reverse characteristics, with a maximum near the hypocenter, and a magnitude of about 23.3 m. Some strike-slip behavior may occur on the two sides of the peak rupture zone. Almost 90% of the seismic moments released by the main shock occurred at depths less than 40 km. The energy released by the fault slip in the 8 h following the main shock is approximately equal to an earthquake of M w 8.13. With a maximum of ~1.5 m, the post-seismic slip was concentrated in the southwestern part of the coseismic rupture fault, which agrees well with the location and behavior of the M w 7.9 aftershock. This implies that the post-seismic deformation in the 8 h after the main shock was mainly induced by the M w 7.9 aftershock. In addition, a post-seismic slip of 0.2-0.4 m was observed at the down-dip extension of the coseismic rupture, which may have been caused by the effect of after-slip during this period.     

Rupture process of the Wenchuan <Emphasis Type="Italic">M</Emphasis>8.0 earthquake of Sichuan China: the segmentation feature

Moment tensor solution, rupture process and rupture characteristics of the great Wenchuan M8.0 earthquake are studied by using 39 long-period P and SH waveforms with evenly azimuth coverage of stations. Our results reveal that the Wenchuan M8.0 event consisted of 5 sub-events of Mw≥7.3 occurring succesively in time and space. Rupture started with a Mw7.3 introductory strike-slip faulting in the first 12 s, then within 12?40 s, two sub-events with Mw7.6 and Mw7.4 occurred within 80 km northeast from the init...