Non-volcanic tremor driven by large transient shear stresses

Non-impulsive seismic radiation or 'tremor' has long been observed at volcanoes and more recently around subduction zones. Although the number of observations of non-volcanic tremor is steadily increasing, the causative mechanism remains unclear. Some have attributed non-volcanic tremor to the movement of fluids, while its coincidence with geodetically observed slow-slip events at regular intervals has led others to consider slip on the plate interface as its cause. Low-frequency earthquakes in Japan, which are believed to make up at least part of non-volcanic tremor, have focal mechanisms and locations that are consistent with tremor being generated by shear slip on the subduction interface. In Cascadia, however, tremor locations appear to be more distributed in depth than in Japan, making them harder to reconcile with a plate interface shear-slip model. Here we identify bursts of tremor that radiated from the Cascadia subduction zone near Vancouver Island, Canada, during the strongest shaking from the moment magnitude M(w) = 7.8, 2002 Denali, Alaska, earthquake. Tremor occurs when the Love wave displacements are to the southwest (the direction of plate convergence of the overriding plate), implying that the Love waves trigger the tremor. We show that these displacements correspond to shear stresses of approximately 40 kPa on the plate interface, which suggests that the effective stress on the plate interface is very low. These observations indicate that tremor and possibly slow slip can be instantaneously induced by shear stress increases on the subduction interface-effectively a frictional failure response to the driving stress.     

Non-volcanic tremor and low-frequency earthquake swarms

Non-volcanic tremor is a weak, extended duration seismic signal observed episodically on some major faults, often in conjunction with slow slip events. Such tremor may hold the key to understanding fundamental processes at the deep roots of faults, and could signal times of accelerated slip and hence increased seismic hazard. The mechanism underlying the generation of tremor and its relationship to aseismic slip are, however, as yet unresolved. Here we demonstrate that tremor beneath Shikoku, Japan, can be explained as a swarm of small, low-frequency earthquakes, each of which occurs as shear faulting on the subduction-zone plate interface. This suggests that tremor and slow slip are different manifestations of a single process.     

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.     

A scaling law for slow earthquakes

Recently, a series of unusual earthquake phenomena have been discovered, including deep episodic tremor, low-frequency earthquakes, very-low-frequency earthquakes, slow slip events and silent earthquakes. Each of these has been demonstrated to arise from shear slip, just as do regular earthquakes, but with longer characteristic durations and radiating much less seismic energy. Here we show that these slow events follow a simple, unified scaling relationship that clearly differentiates their behaviour from that of regular earthquakes. We find that their seismic moment is proportional to the characteristic duration and their moment rate function is constant, with a spectral high-frequency decay of f(-1). This scaling and spectral behaviour demonstrates that they can be thought of as different manifestations of the same phenomena and that they comprise a new earthquake category. The observed scale dependence of rupture velocity for these events can be explained by either a constant low-stress drop model or a diffusional constant-slip model. This new scaling law unifies a diverse class of slow seismic events and may lead to a better understanding of the plate subduction process and large earthquake generation.     

Earthquakes triggered by silent slip events on Kīlauea volcano, Hawaii

Slow-slip events, or 'silent earthquakes', have recently been discovered in a number of subduction zones including the Nankai trough in Japan, Cascadia, and Guerrero in Mexico, but the depths of these events have been difficult to determine from surface deformation measurements. Although it is assumed that these silent earthquakes are located along the plate megathrust, this has not been proved. Slow slip in some subduction zones is associated with non-volcanic tremor, but tremor is difficult to locate and may be distributed over a broad depth range. Except for some events on the San Andreas fault, slow-slip events have not yet been associated with high-frequency earthquakes, which are easily located. Here we report on swarms of high-frequency earthquakes that accompany otherwise silent slips on Kīlauea volcano, Hawaii. For the most energetic event, in January 2005, the slow slip began before the increase in seismicity. The temporal evolution of earthquakes is well explained by increased stressing caused by slow slip, implying that the earthquakes are triggered. The earthquakes, located at depths of 7-8 km, constrain the slow slip to be at comparable depths, because they must fall in zones of positive Coulomb stress change. Triggered earthquakes accompanying slow-slip events elsewhere might go undetected if background seismicity rates are low. Detection of such events would help constrain the depth of slow slip, and could lead to a method for quantifying the increased hazard during slow-slip events, because triggered events have the potential to grow into destructive earthquakes.     

Reflection signature of seismic and aseismic slip on the northern Cascadia subduction interface

At the northern Cascadia margin, the Juan de Fuca plate is underthrusting North America at about 45 mm x yr(-1) (ref. 1), resulting in the potential for destructive great earthquakes. The downdip extent of coupling between the two plates is difficult to determine because the most recent such earthquake (thought to have been in 1700) occurred before instrumental recording. Thermal and deformation studies indicate that, off southern Vancouver Island, the interplate interface is presently fully locked for a distance of approximately 60 km downdip from the deformation front. Great thrust earthquakes on this section of the interface (with magnitudes of up to 9) have been estimated to occur at an average interval of about 590 yr (ref. 3). Further downdip there is a transition from fully locked behaviour to aseismic sliding (where high temperatures allow ductile deformation), with the deep aseismic zone exhibiting slow-slip thrust events. Here we show that there is a change in the reflection character on seismic images from a thin sharp reflection where the subduction thrust is inferred to be locked, to a broad reflection band at greater depth where aseismic slip is thought to be occurring. This change in reflection character may provide a new technique to map the landward extent of rupture in great earthquakes and improve the characterization of seismic hazards in subduction zones.     

Seismic reflection imaging of two megathrust shear zones in the northern Cascadia subduction zone

At convergent continental margins, the relative motion between the subducting oceanic plate and the overriding continent is usually accommodated by movement along a single, thin interface known as a megathrust. Great thrust earthquakes occur on the shallow part of this interface where the two plates are locked together. Earthquakes of lower magnitude occur within the underlying oceanic plate, and have been linked to geochemical dehydration reactions caused by the plate's descent. Here I present deep seismic reflection data from the northern Cascadia subduction zone that show that the inter-plate boundary is up to 16 km thick and comprises two megathrust shear zones that bound a >5-km-thick, approximately 110-km-wide region of imbricated crustal rocks. Earthquakes within the subducting plate occur predominantly in two geographic bands where the dip of the plate is inferred to increase as it is forced around the edges of the imbricated inter-plate boundary zone. This implies that seismicity in the subducting slab is controlled primarily by deformation in the upper part of the plate. Slip on the shallower megathrust shear zone, which may occur by aseismic slow slip, will transport crustal rocks into the upper mantle above the subducting oceanic plate and may, in part, provide an explanation for the unusually low seismic wave speeds that are observed there.     

A wide depth distribution of seismic tremors along the northern Cascadia margin

The Cascadia subduction zone is thought to be capable of generating major earthquakes with moment magnitude as large as M(w) = 9 at an interval of several hundred years. The seismogenic portion of the plate interface is mostly offshore and is currently locked, as inferred from geodetic data. However, episodic surface displacements-in the direction opposite to the long-term deformation motions caused by relative plate convergence across a locked interface-are observed about every 14 months with an unusual tremor-like seismic signature. Here we show that these tremors are distributed over a depth range exceeding 40 km within a limited horizontal band. Many occurred within or close to the strong seismic reflectors above the plate interface where local earthquakes are absent, suggesting that the seismogenic process for tremors is fluid-related. The observed depth range implies that tremors could be associated with the variation of stress field induced by a transient slip along the deeper portion of the Cascadia interface or, alternatively, that episodic slip is more diffuse than originally suggested.     

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.     

汤加-克马德克俯冲带的地质构造与地震火山特征

2022年1月15日西南太平洋的洪阿哈阿帕伊岛海底火山发生了爆炸式的剧烈喷发，吸引了全球的关注。洪阿哈阿帕伊岛海底火山位于汤加-克马德克俯冲带，综合前期研究结果，对汤加-克马德克俯冲带的地质构造特征、地震和火山分布进行初步分析，发现：（1）从汤加-克马德克俯冲带弧前向海方向直到俯冲的太平洋板块，构造上主要表现为大规模正断层。（2）路易斯维尔海山链的俯冲将汤加-克马德克俯冲带分为北部的汤加俯冲带和南部的克马德克俯冲带，沿汤加俯冲带板块汇聚率为67～84 mm/a，沿克马德克俯冲带板块汇聚率为41～58 mm/a，板块俯冲速度的差异造成汤加俯冲带和克马德克俯冲带目前俯冲深度的不同。（3）在路易斯维尔海山链以北，太平洋板块上覆沉积物厚度不足0.4 km，而在南侧达到1 km左右，由于俯冲板块上覆沉积物厚度的差异而造成北部的汤加俯冲带和南部的克马德克俯冲带孕育地震能力的差异。这些认识对研究该俯冲带的火山喷发机制、大地震成因机理及其灾害风险具有重要意义。     

西太平洋俯冲带岩石圈变形研究

全球海洋岩石圈的最大弯曲与大地震发生在俯冲带。当弯曲应力超过岩石承受范围，就会产生正断层和地震，海水沿着正断层进入上地幔并发生蛇纹石化，引发浅源地震并可能造成灾难性海啸。选取西太平洋最具代表性的日本、伊豆-小笠原、马里亚纳和雅浦俯冲带以及汤加-克马德克俯冲带，归纳近些年的地球物理观测及地球动力学模拟的结果，对比分析了不同俯冲带挠曲正断层的分布特征，并探讨了俯冲板块变形与地震之间的相关性,以揭示俯冲板片弯曲变形及相应的正断层与潜在板块水化特征。     

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.     

Geological record of fluid flow and seismogenesis along an erosive subducting plate boundary

Tectonic erosion of the overriding plate by the downgoing slab is believed to occur at half the Earth's subduction zones. In situ investigation of the geological processes at active erosive margins is extremely difficult owing to the deep marine environment and the net loss of forearc crust to deeper levels in the subduction zone. Until now, a fossil erosive subduction channel-the shear zone marking the plate boundary-has not been recognized in the field, so that seismic observations have provided the only information on plate boundary processes at erosive margins. Here we show that a fossil erosive margin is preserved in the Northern Apennines of Italy. It formed during the Tertiary transition from oceanic subduction to continental collision, and was preserved by the late deactivation and fossilization of the plate boundary. The outcropping erosive subduction channel is approximately 500 m thick. It is representative of the first 5 km of depth, with its deeper portions reaching approximately 150 degrees C. The fossil zone records several surprises. Two décollements were simultaneously active at the top and base of the subduction channel. Both deeper basal erosion and near-surface frontal erosion occurred. At shallow depths extension was a key deformation component within this erosive convergent plate boundary, and slip occurred without an observable fluid pressure cycle. At depths greater than about 3 km a fluid cycle is clearly shown by the development of veins and the alternation of fast (co-seismic) and slow (inter-seismic) slip. In the deepest portions of the outcropping subduction channel, extension is finally overprinted by compressional structures. In modern subduction zones the onset of seismic activity is believed to occur at approximately 150 degrees C, but in the fossil channel the onset occurred at cooler palaeo-temperatures.     

“311”东北日本M9.0大地震可预测性特征探索

 大地震不是任何地方都能发生的,大地震也不是任何时间都能发生的,大地震有其特定的时空结构。前期工作表明,条环交会、差异活动以及深震与“天外来客”事件等是制约发震时空窗口的主要因素。2011年3月11日发生在日本海沟俯冲带的M9.0大地震,其震前地震活动全面、清晰地展现出了其可预测性特征：(1) 1925年开始出现的板间地震长期平静区;(2) 2003年开始出现的位于平静区中部,起始于深震的海沟垂向地震活动条带;(3) 2009年3月(海沟垂向地震活动条带形成后)开始出现的广义前震发震地方时的非随机显著聚集特征;(4) 2011年3月9日(主震前2天)发生在平静区内的前震。主震发生在平静区边缘,海沟垂向地震活动条带与日本海沟板间地震带的交会处的地震活动差异性较大的地方,其断层面与海沟垂向地震活动条带内另外两个板间地震(2005-08-16地震和前震)的断层面共面。主震发生于广义前震发震地方时非随机聚集时段的中位时辰。日本大地震的可预测性特征逐一印证了我们的前期研究。     

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.     

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.     

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.     

Subcontinental-scale crustal velocity changes along the Pacific-North America plate boundary

Transient tectonic deformation has long been noted within approximately 100 km of plate boundary fault zones and within active volcanic regions, but it is unknown whether transient motions also occur at larger scales within plates. Relatively localized transients are known to occur as both seismic and episodic aseismic events, and are generally ascribed to motions of magma bodies, aseismic creep on faults, or elastic or viscoelastic effects associated with earthquakes. However, triggering phenomena and systematic patterns of seismic strain release at subcontinental (approximately 1,000 km) scale along diffuse plate boundaries have long suggested that energy transfer occurs at larger scale. Such transfer appears to occur by the interaction of stresses induced by surface wave propagation and magma or groundwater in the crust, or from large-scale stress diffusion within the oceanic mantle in the decades following clusters of great earthquakes. Here we report geodetic evidence for a coherent, subcontinental-scale change in tectonic velocity along a diffuse approximately 1,000-km-wide deformation zone. Our observations are derived from continuous GPS (Global Positioning System) data collected over the past decade across the Basin and Range province, which absorbs approximately 25 per cent of Pacific-North America relative plate motion. The observed changes in site velocity define a sharp boundary near the centre of the province oriented roughly parallel to the north-northwest relative plate motion vector. We show that sites to the west of this boundary slowed relative to sites east of it by approximately 1 mm yr(-1) starting in late 1999.     

Measuring the onset of locking in the Peru-Chile trench with GPS and acoustic measurements

The subduction zone off the west coast of South America marks the convergence of the oceanic Nazca plate and the continental South America plate. Nazca-South America convergence over the past 23 million years has created the 6-km-deep Peru-Chile trench, 150 km offshore. High pressure between the plates creates a locked zone, leading to deformation of the overriding plate. The surface area of this locked zone is thought to control the magnitude of co-seismic release and is limited by pressure, temperature, sediment type and fluid content. Here we present seafloor deformation data from the submerged South America plate obtained from a combination of Global Positioning System (GPS) receivers and acoustic transponders. We estimate that the measured horizontal surface motion perpendicular to the trench is consistent with a model having no slip along the thrust fault between 2 and 40 km depth. A tsunami in 1996, 200 km north of our site, was interpreted as being the result of an anomalously shallow interplate earthquake. Seismic coupling at shallow depths, such as we observe, may explain why co-seismic events in the Peruvian subduction zone create large tsunamis.     

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.