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.     

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.     

Magma-maintained rift segmentation at continental rupture in the 2005 Afar dyking episode

Seafloor spreading centres show a regular along-axis segmentation thought to be produced by a segmented magma supply in the passively upwelling mantle. On the other hand, continental rifts are segmented by large offset normal faults, and many lack magmatism. It is unclear how, when and where the ubiquitous segmented melt zones are emplaced during the continental rupture process. Between 14 September and 4 October 2005, 163 earthquakes (magnitudes greater than 3.9) and a volcanic eruption occurred within the approximately 60-km-long Dabbahu magmatic segment of the Afar rift, a nascent seafloor spreading centre in stretched continental lithosphere. Here we present a three-dimensional deformation field for the Dabbahu rifting episode derived from satellite radar data, which shows that the entire segment ruptured, making it the largest to have occurred on land in the era of satellite geodesy. Simple elastic modelling shows that the magmatic segment opened by up to 8 m, yet seismic rupture can account for only 8 per cent of the observed deformation. Magma was injected along a dyke between depths of 2 and 9 km, corresponding to a total intrusion volume of approximately 2.5 km3. Much of the magma appears to have originated from shallow chambers beneath Dabbahu and Gabho volcanoes at the northern end of the segment, where an explosive fissural eruption occurred on 26 September 2005. Although comparable in magnitude to the ten year (1975-84) Krafla events in Iceland, seismic data suggest that most of the Dabbahu dyke intrusion occurred in less than a week. Thus, magma intrusion via dyking, rather than segmented normal faulting, maintains and probably initiated the along-axis segmentation along this sector of the Nubia-Arabia plate boundary.     

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.     

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.     

Low-frequency earthquakes in Shikoku, Japan, and their relationship to episodic tremor and slip

Non-volcanic seismic tremor was discovered in the Nankai trough subduction zone in southwest Japan and subsequently identified in the Cascadia subduction zone. In both locations, tremor is observed to coincide temporally with large, slow slip events on the plate interface downdip of the seismogenic zone. The relationship between tremor and aseismic slip remains uncertain, however, largely owing to difficulty in constraining the source depth of tremor. In southwest Japan, a high quality borehole seismic network allows identification of coherent S-wave (and sometimes P-wave) arrivals within the tremor, whose sources are classified as low-frequency earthquakes. As low-frequency earthquakes comprise at least a portion of tremor, understanding their mechanism is critical to understanding tremor as a whole. Here, we provide strong evidence that these earthquakes occur on the plate interface, coincident with the inferred zone of slow slip. The locations and characteristics of these events suggest that they are generated by shear slip during otherwise aseismic transients, rather than by fluid flow. High pore-fluid pressure in the immediate vicinity, as implied by our estimates of seismic P- and S-wave speeds, may act to promote this transient mode of failure. Low-frequency earthquakes could potentially contribute to seismic hazard forecasting by providing a new means to monitor slow slip at depth.     

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.     

中国陆域近10年地震时空分布统计特征

对中国陆域近10年4级以上地震震中的空间分布和季节分布特征进行了统计分析，从1994年1月到2005年2月间的4级以上地震，共计966次。近10年的地震空间分布与中国陆域历史地震分布基本吻合。中国陆域地震活动与印度板块、太平洋板块及欧亚板块间的相互作用有关，与前者相关的地震位于青藏高原及其外围地区；与后者密切的地震活动集中在中国陆域南海，尤其是台湾省。同时地震活动多围绕断块的周围分布。地震震中较集中地分布在活动断裂带附近、断裂带密集分布地带以及大构造区域的边界地带。地震年际活动的季节统计分析表明，地震高峰期集中出现在春季和秋季，各年地震活动波动性最大，地震距平曲线与地球自转速度变化曲线协同，地球自转速度的季节性变化是影响地震活动的基本因素之一。     

Discovery of a magma chamber and faults beneath a Mid-Atlantic Ridge hydrothermal field

Crust at slow-spreading ridges is formed by a combination of magmatic and tectonic processes, with magmatic accretion possibly involving short-lived crustal magma chambers. The reflections of seismic waves from crustal magma chambers have been observed beneath intermediate and fast-spreading centres, but it has been difficult to image such magma chambers beneath slow-spreading centres, owing to rough seafloor topography and associated seafloor scattering. In the absence of any images of magma chambers or of subsurface near-axis faults, it has been difficult to characterize the interplay of magmatic and tectonic processes in crustal accretion and hydrothermal circulation at slow-spreading ridges. Here we report the presence of a crustal magma chamber beneath the slow-spreading Lucky Strike segment of the Mid-Atlantic Ridge. The reflection from the top of the magma chamber, centred beneath the Lucky Strike volcano and hydrothermal field, is approximately 3 km beneath the sea floor, 3-4 km wide and extends up to 7 km along-axis. We suggest that this magma chamber provides the heat for the active hydrothermal vent field above it. We also observe axial valley bounding faults that seem to penetrate down to the magma chamber depth as well as a set of inward-dipping faults cutting through the volcanic edifice, suggesting continuous interactions between tectonic and magmatic processes.     

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.     

Boundary-layer mantle flow under the Dead Sea transform fault inferred from seismic anisotropy

Lithospheric-scale transform faults play an important role in the dynamics of global plate motion. Near-surface deformation fields for such faults are relatively well documented by satellite geodesy, strain measurements and earthquake source studies, and deeper crustal structure has been imaged by seismic profiling. Relatively little is known, however, about deformation taking place in the subcrustal lithosphere--that is, the width and depth of the region associated with the deformation, the transition between deformed and undeformed lithosphere and the interaction between lithospheric and asthenospheric mantle flow at the plate boundary. Here we present evidence for a narrow, approximately 20-km-wide, subcrustal anisotropic zone of fault-parallel mineral alignment beneath the Dead Sea transform, obtained from an inversion of shear-wave splitting observations along a dense receiver profile. The geometry of this zone and the contrast between distinct anisotropic domains suggest subhorizontal mantle flow within a vertical boundary layer that extends through the entire lithosphere and accommodates the transform motion between the African and Arabian plates within this relatively narrow zone.     

A sharp lithosphere-asthenosphere boundary imaged beneath eastern North America

Plate tectonic theory hinges on the concept of a relatively rigid lithosphere moving over a weaker asthenosphere, yet the nature of the lithosphere-asthenosphere boundary remains poorly understood. The gradient in seismic velocity that occurs at this boundary is central to constraining the physical and chemical properties that create differences in mechanical strength between the two layers. For example, if the lithosphere is simply a thermal boundary layer that is more rigid owing to colder temperatures, mantle flow models indicate that the velocity gradient at its base would occur over tens of kilometres. In contrast, if the asthenosphere is weak owing to volatile enrichment or the presence of partial melt, the lithosphere-asthenosphere boundary could occur over a much smaller depth range. Here we use converted seismic phases in eastern North America to image a very sharp seismic velocity gradient at the base of the lithosphere-a 3-11 per cent drop in shear-wave velocity over a depth range of 11 km or less at 90-110 km depth. Such a strong, sharp boundary cannot be reconciled with a purely thermal gradient, but could be explained by an asthenosphere that contains a few per cent partial melt or that is enriched in volatiles relative to the lithosphere.     

Focal depth research of earthquakes in Mainland China: Implication for tectonics

Focal depth data of earthquakes in Mainland China are processed and analyzed in this paper, as well as the relationship between the focal depths and large-scale tectonic structures. As a basic parameter for earthquakes, focal depth is used to investigate deep environment of seismogenic regions, tectonic backgrounds for concentration and release of seismic energy, the inner crustal deformation and its mechanic features. Depth data of 31282 ML≥2.0 events with 1st class and 2nd class precision in Mainland China from Jan. 1, 1970 to May 31, 2000 are used to get spatial features of earthquakes distributed with depth and to provide average depth for each grid area throughout China. Researches show that the average depth (D-) for all the earthquakes used in this paper is (16±7) km, and (13±6) km and (18±8) km for the events in eastern China and western China, respectively. The area with the deepest focal depth is located in southwest Xinjiang region, near the western and southwestern ends of the Tarim Basin. The focal depth related to large-scale tectonic structures, for instance, = (33±12), (21±10), (14±7), (11±5) and (10±4) km in Tibet plateau block, Xinjiang block, North China, Northeastern China and South China, respectively. The earthquakes are deeper at the bounders of the integrated tectonic blocks, including the southwestern and northern brims of the Tarim Basin, southern brim of the Zhunge'r Basin and that of the Alashan block, as well as the eastern and western sides of the Edos block and the western brim of the Sichuan Basin. The earthquakes at the newly ruptured belts are relatively shallower, for instance, at the southwestern Yunnan seismic belt and the Zhangjiakou-Bohai seismic belt. The mechanic behavior, deformation and features for the crust and mantle structures are also discussed.     

Evidence for seismogenic fracture of silicic magma

It has long been assumed that seismogenic faulting is confined to cool, brittle rocks, with a temperature upper limit of approximately 600 degrees C (ref. 1). This thinking underpins our understanding of volcanic earthquakes, which are assumed to occur in cold rocks surrounding moving magma. However, the recent discovery of abundant brittle-ductile fault textures in silicic lavas has led to the counter-intuitive hypothesis that seismic events may be triggered by fracture and faulting within the erupting magma itself. This hypothesis is supported by recent observations of growing lava domes, where microearthquake swarms have coincided with the emplacement of gouge-covered lava spines, leading to models of seismogenic stick-slip along shallow shear zones in the magma. But can fracturing or faulting in high-temperature, eruptible magma really generate measurable seismic events? Here we deform high-temperature silica-rich magmas under simulated volcanic conditions in order to test the hypothesis that high-temperature magma fracture is seismogenic. The acoustic emissions recorded during experiments show that seismogenic rupture may occur in both crystal-rich and crystal-free silicic magmas at eruptive temperatures, extending the range of known conditions for seismogenic faulting.     

含有粘土矿物断层泥地震与非地震标志

实验变形的含有粘土矿物断层泥的显微结构观察提供了判别自然界断层滑动习性的标志. 含有粘土矿物断层泥的摩擦实验显示了2种模式, 一种是稳态滑动,另一种是粘滑,即有可能发生大地震. 此两种实验产物的显微结构记录了断层的滑动性质,均匀变形的断层泥是断层非地震( 稳滑) 的结果, 局部强烈变形是断层地震( 粘滑) 的产物. 均匀变形使粘土矿物产生优选方位( S 叶理) , 在一定程度上它反映了剪应变的大小, 局部变形形成吕德剪切和 Y 剪切,云南小湾 F7 断层中具有类似的结构特征,计算结果表明 F7 的稳滑量约 43%~ 52%.     

Pliocene eclogite exhumation at plate tectonic rates in eastern Papua New Guinea

As lithospheric plates are subducted, rocks are metamorphosed under high-pressure and ultrahigh-pressure conditions to produce eclogites and eclogite facies metamorphic rocks. Because chemical equilibrium is rarely fully achieved, eclogites may preserve in their distinctive mineral assemblages and textures a record of the pressures, temperatures and deformation the rock was subjected to during subduction and subsequent exhumation. Radioactive parent-daughter isotopic variations within minerals reveal the timing of these events. Here we present in situ zircon U/Pb ion microprobe data that dates the timing of eclogite facies metamorphism in eastern Papua New Guinea at 4.3 +/- 0.4 Myr ago, making this the youngest documented eclogite exposed at the Earth's surface. Eclogite exhumation from depths of approximately 75 km was extremely rapid and occurred at plate tectonic rates (cm yr(-1)). The eclogite was exhumed within a portion of the obliquely convergent Australian-Pacific plate boundary zone, in an extending region located west of the Woodlark basin sea floor spreading centre. Such rapid exhumation (> 1 cm yr(-1)) of high-pressure and, we infer, ultrahigh-pressure rocks is facilitated by extension within transient plate boundary zones associated with rapid oblique plate convergence.     

Seismic evidence for convection-driven motion of the North American plate

Since the discovery of plate tectonics, the relative importance of driving forces of plate motion has been debated. Resolution of this issue has been hindered by uncertainties in estimates of basal traction, which controls the coupling between lithospheric plates and underlying mantle convection. Hotspot tracks preserve records of past plate motion and provide markers with which the relative motion between a plate's surface and underlying mantle regions may be examined. Here we show that the 115-140-Myr surface expression of the Great Meteor hotspot track in eastern North America is misaligned with respect to its location at 200 km depth, as inferred from plate-reconstruction models and seismic tomographic studies. The misalignment increases with age and is consistent with westward displacement of the base of the plate relative to its surface, at an average rate of 3.8 +/- 1.8 mm yr(-1). Here age-constrained 'piercing points' have enabled direct estimation of relative motion between the surface and underside of a plate. The relative displacement of the base is approximately parallel to seismic fast axes and calculated mantle flow, suggesting that asthenospheric flow may be deforming the lithospheric keel and exerting a driving force on this part of the North American plate.     

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.     

The role of fluids in lower-crustal earthquakes near continental rifts

The occurrence of earthquakes in the lower crust near continental rifts has long been puzzling, as the lower crust is generally thought to be too hot for brittle failure to occur. Such anomalous events have usually been explained in terms of the lower crust being cooler than normal. But if the lower crust is indeed cold enough to produce earthquakes, then the uppermost mantle beneath it should also be cold enough, and yet uppermost mantle earthquakes are not observed. Numerous lower-crustal earthquakes occur near the southwestern termination of the Taupo Volcanic Zone (TVZ), an active continental rift in New Zealand. Here we present three-dimensional tomographic imaging of seismic velocities and seismic attenuation in this region using data from a dense seismograph deployment. We find that crustal earthquakes accurately relocated with our three-dimensional seismic velocity model form a continuous band along the rift, deepening from mostly less than 10 km in the central TVZ to depths of 30-40 km in the lower crust, 30 km southwest of the termination of the volcanic zone. These earthquakes often occur in swarms, suggesting fluid movement in critically loaded fault zones. Seismic velocities within the band are also consistent with the presence of fluids, and the deepening seismicity parallels the boundary between high seismic attenuation (interpreted as partial melt) within the central TVZ and low seismic attenuation in the crust to the southwest. This linking of upper and lower-crustal seismicity and crustal structure allows us to propose a common explanation for all the seismicity, involving the weakening of faults on the periphery of an otherwise dry, mafic crust by hot fluids, including those exsolved from underlying melt. Such fluids may generally be an important driver of lower-crustal seismicity near continental rifts.