Double Seismic Zone Research Articles

Anisotropic thermal conductivity of antigorite along slab subduction impacts seismicity of intermediate-depth earthquakes

Double seismic zones (DSZs) are a feature of some subducting slabs, where intermediate-depth earthquakes (~70–300 km) align along two separate planes. The upper seismic plane is generally attributed to dehydration embrittlement, whereas mechanisms forming the lower seismic plane are still debated. Thermal conductivity of slab minerals is expected to control the temperature evolution of subducting slabs, and therefore their seismicity. However, effects of the potential anisotropic thermal conductivity of layered serpentine minerals with crystal preferred orientation on slab’s thermal evolution remain poorly understood. Here we measure the lattice thermal conductivity of antigorite, a hydrous serpentine mineral, along its crystallographic b- and c-axis at relevant high pressure-temperature conditions of subduction. We find that antigorite’s thermal conductivity along the c-axis is ~3–4 folds smaller than the b-axis. Our numerical models further reveal that when the low-thermal-conductivity c-axis is aligned normal to the slab dip, antigorite’s strongly anisotropic thermal conductivity enables heating at the top portion of the slab, facilitating dehydration embrittlement that causes the seismicity in the upper plane of DSZs. Potentially, the antigorite’s thermal insulating effect also hinders the dissipation of frictional heat inside shear zones, promoting thermal runaway along serpentinized faults that could trigger intermediate-depth earthquakes.

Evidence From Intermediate‐Depth Earthquakes of Slab‐Derived Fluids Beneath the Taupō Volcanic Zone

AbstractAhi Tupua, the central section of the Taupō Volcanic Zone in the central North Island of New Zealand, encompasses Taupō and Ōkataina calderas and has been the most frequently active and productive silicic magma system worldwide during the Quaternary. The entire Taupō Volcanic Zone is underlain by the Hikurangi Subduction Zone where a Large Igneous Province, the Hikurangi Plateau, is being subducted, but Ahi Tupua exhibits much higher rates of magma output than either the northern or the southern sections. Trace element signatures of Ahi Tupua eruptive products suggest that both decompression melting and flux melting play a role in this rifted arc setting, with fluid flux signatures being more prevalent in regions of active caldera volcanism. Intermediate‐depth (50–300 km) earthquakes provide a means of studying faulting and fluid flow processes within and around the subducting slab beneath Ahi Tupua. Using data from the national seismic network and a 13‐station temporary network (“ECLIPSE”), we repick and relocate 397 intermediate‐depth earthquakes of magnitudes M2.5+ that occurred in a 29‐month interval, and compute focal mechanisms for a subset of 47 earthquakes. We observe some seismicity that may be within the mantle wedge but most of the relocated earthquakes occur within the crust and uppermost mantle of the subducted slab and exhibit a weak transition from predominantly normal‐faulting in the slab crust to predominantly reverse‐faulting in the slab mantle. No double seismic zone is observed beneath Ahi Tupua. These observations are consistent with dehydration of the slab and flow of slab‐derived fluids along existing faults into the mantle wedge, that drives flux melting and accounts for the distinctive geochemical signatures and voluminous output of Ahi Tupua calderas.

Impact of chlorite dehydration on intermediate-depth earthquakes in subducting slabs

Intermediate-depth earthquakes are common in the double seismic structures of many subduction zones under high pressures (~1–4 GPa). Serpentine dehydration exhibits well-established links with double seismic zone earthquakes. Additionally, dehydration of several hydrous minerals including lawsonite and chlorite underlying the upper and lower layers, respectively, may be responsible for intermediate-depth earthquakes. Here, we present experimental evidence suggesting that chlorite dehydration can trigger intermediate-depth earthquakes at the lower plane (~700 °C). We conducted deformation experiments on chlorite peridotite under high-pressure (0.5–2.5 GPa) and high-temperature (500–750 °C) using a modified Griggs apparatus. Experiments revealed the presence of faults in samples that had undergone partial chlorite dehydration with the presence of the dehydration product Ca-amphibole along these faults. Our findings confirm, together with correlation studies between seismicity and mineral stability, that a part of intermediate-depth seismicity in the lower plane of double seismic zones can be attributed to chlorite dehydration.

Scaling of earthquake waiting time distributions in northern Chile

SUMMARY In this study, we examine the scaling of earthquake waiting time distributions in Northern Chile. For this purpose, we analysed 7-yr high spatial resolution and low completeness magnitude IPOC seismic catalogue and 43-yr USGS catalogue. A unified earthquake moment–space–time scaling relation is empirically evaluated by calculating waiting times for different ranges of magnitude and epicentral area linear dimension and estimating the associated scaling coefficients, β analogue to b-value and γ, the correlation fractal dimension. We find a scaling function that can be characterized with three distinct regions, regions whose behaviour depends on whether seismicity is in the coastal area or at intermediate depth. Moreover, high resolution localizations from the IPOC catalogue allows us to further observe differences: the lower plane seismicity of the double seismic zone behaves just like that at intermediate depth, while the crustal upper plate and the upper plane from the double seismic zone seismicity behaves more like interface seismicity. Thus, the earthquake waiting time distribution primarily depends on whether seismicity is located in a crust-on-crust or in a crust–mantle system: having, respectively, high/low correlated behaviour in the short scale region, non-exponential/exponential decay in the transition middle region and, in all cases, long-term clustering with a slower than exponential decay in the long scale.

Spatiotemporal Variations of Intermediate‐Depth Earthquakes Before and After 2011 Tohoku Earthquake Revealed by a Template Matching Catalog

AbstractWe investigate spatiotemporal changes of intermediate‐depth earthquakes in the double seismic zone beneath Central and Northeastern Japan before and after the 2011 magnitude 9 Tohoku earthquake. We build a template‐matching catalog 1 year before and 1 year after the Tohoku earthquake using Hi‐net recordings. The new catalog has a six‐fold increase in earthquakes compared to the Japan Meteorological Agency catalog. Our results show no significant change in the intermediate‐depth earthquake rate prior to the Tohoku earthquake, but a clear increase in both planes following the Tohoku earthquake. The regions with increased intermediate‐depth earthquake activity and the post‐seismic slips following the Tohoku earthquake are spatially separate and complementary with each other. Aftershock productivity of intermediate‐depth earthquakes increased in both planes following the Tohoku earthquake. Overall, aftershock productivity of the upper plane is higher than the lower plane, likely indicating that stress environments and physical mechanisms of intermediate‐depth earthquakes in the two planes are distinct.

Anomalous intraslab structure revealed by the analysis of aftershocks of the Mw 6.7 Coquimbo-La Serena earthquake of 20 January 2019

Seismograms of about 12,000 aftershocks of the Mw6.7 Coquimbo-La Serena earthquake of 20 January 2019 that were recorded during six months by an array of 33 short-period seismic stations are used to generate wavespeed images of the subduction wedge beneath the Chilean forearc between 28.0°-30.6°S. This part of the margin is near the northern terminus of the 2015 Mw8.4 Illapel earthquake and the southern terminus of the 1922 Mw8.5 Copiapo earthquake. The distribution of hypocenters and wavespeeds at the slab interface and within the subduction wedge, are like those determined by Comte et al. (2019) in the Illapel aftershock region, and we infer that the same processes of subduction erosion and re-accretion inferred by them are active here. Moreover, anomalous wavespeeds are found in the immediate vicinity of a dense cluster of aftershocks located in a well-defined double seismic zone, within the subducting Nazca plate. Vp/Vs ratios in the upper seismic zone are uniformly high (∼1.9), those in most of the lower seismic zone are average to low, and those in the intraslab region of the aftershock cluster are exceptionally low (∼1.6). These values suggest dehydration of the slab interior and hydration of the mantle near the slab interface. While the origin of this feature is enigmatic, its limited spatial extent suggests that it most likely is due to an inherited structure, such a seamount or petit-spot volcano, in the Nazca plate, rather than a feature that developed in situ. The high rate of activity coupled with the low wavespeeds suggests a region of local weakness that could eventually evolve into a tear or hole in the slab, and its location at the bounds of two megathrust earthquakes suggests that it plays a significant role in delimiting the rupture of such events.

Tectonic deformation at the outer rise of subduction zones

SUMMARY Fluids associated with subducting slabs play a crucial role in regulating the dynamics of water discharge, subsequent arc magmatism and intermediate-depth earthquakes in subduction zones. The incoming slab mantle hydration is primarily determined by deep normal faulting due to plate bending at the trench. However, the controlling factors on the outer rise faulting pattern, and the correlation between the inherited outer rise deformation and the intermediate-depth earthquakes, remain to be understood. Here we present high-resolution viscoelasto-plastic numerical models of free subduction for slab bending-related faulting prior to subduction. Our model results show that plastic weakening and friction coefficient of the slab mantle exhibit a significant impact on fault pattern, while plate age and elasticity have a minimal bearing for mature slabs. The brittle bending faults result in large positive pressure gradients in the vertical direction, facilitating seawater infiltrating into the subducting slabs, which corroborates previous numerical models. The faults reaching 15–30 km beneath the Moho coincide with the width of the double seismic zone in subduction zones. We anticipate that water pumped into the slab mantle along the faults, with decreasing water content along the depth, can explain the relatively sporadic lower plane earthquakes.

Stress Drops of Intermediate‐Depth and Deep Earthquakes in the Tonga Slab

AbstractMultiple physical mechanisms have been proposed to explain the cause of intermediate‐depth and deep earthquakes, but they are still under debate. Source parameters such as stress drop, have the potential to provide insight into their physical mechanisms. We develop a modified spectral decomposition method to analyze 1‐year seismic data from temporary land‐based and ocean bottom seismographs in a complex subduction zone. By applying this method to investigate 1,083 intermediate‐depth and deep earthquakes in the Tonga slab, we successfully resolve the source spectra and stress drops of 743 MW 2.6–6.0 earthquakes. Although the absolute stress drops are subject to the choices of source model parameters, the relative stress drops are more reliably resolved. The median stress drop of Tonga earthquakes does not change with respect to magnitude but decreases with depth by 2–3 times in two separate depth ranges of 70–250 and 400–600 km, corresponding to intermediate‐depth and deep earthquakes, respectively. The median stress drops show spatial variations, with two high‐stress‐drop (five times higher than the ambient value) regions, coinciding with strong local deformation where the Tonga slab bends or tears. In the Tonga double seismic zone at 120–300 km depths, the median stress drop appears smaller in the lower plane than in the upper plane, suggesting a slower rupture velocity or a higher fluid content in the lower‐plane region. Our results suggest that intermediate‐depth and deep earthquakes in the Tonga slab generally follow the earthquake self‐similar model and favor the fluid‐related embrittlement hypothesis for both groups of earthquakes.

Global Constraints on Intermediate‐Depth Intraslab Stresses From Slab Geometries and Mechanisms of Double Seismic Zone Earthquakes

AbstractDouble seismic zones (DSZs), parallel planes of intermediate‐depth earthquakes inside oceanic slabs, have been observed in a number of subduction zones and may be a ubiquitous feature of downgoing oceanic plates. Focal mechanism observations from DSZ earthquakes sample the intraslab stress field at two distinct depth levels within the downgoing lithosphere. A pattern of downdip compressive over downdip extensive events was early on interpreted to indicate an unbending‐dominated intraslab stress field. In the present study, we show that the intraslab stress field in the depth range of DSZs is much more variable than previously thought. Compiling DSZ locations and mechanisms from literature, we observe that the “classical” pattern of compressive over extensive events is only observed at about half of the DSZ locations around the globe. The occurrence of extensional mechanisms across both planes accounts for most other regions. To obtain an independent estimate of the bending state of slabs at intermediate depths, we compute (un)bending estimates from slab geometries taken from the slab2 compilation of slab surface depths. We find no clear global prevalence of slab unbending at intermediate depths, and the occurrence of DSZ seismicity does not appear to be limited to regions of slab (un)bending. Focal mechanism observations are frequently inconsistent with (un)bending estimates from slab geometries, which may imply that bending stresses are not always prevalent, and that other stress types such as in‐plane tension due to slab pull or shallow compression due to friction along the plate interface may also play an important role.

Double Seismic Zones along the Eastern Aleutian-Alaska Subduction Zone Revealed by a High-Precision Earthquake Relocation Catalog

Abstract The Eastern Aleutian-Alaska Subduction Zone (EAASZ) manifests significant along-strike variations in structure and geometry. The limited spatial resolution in intermediate-depth earthquake locations precludes investigation of small-scale variations in seismic characteristics. In this study, we use an existing 3D seismic velocity model and waveform cross-correlation data to relocate the earthquakes in 2016 near the EAASZ. Our improved absolute and relative earthquake locations reveal complex spatial characteristics of double seismic zones (DSZs). There are significant variations in location, depth, layer separation, and length of the DSZs along the EAASZ. We also observe nonuniform layer separations along the slope of the subducting slab that may imply either rheological or crustal thickness variations. In addition, our results suggest a triple seismic zone (TSZ) beneath Kenai. The interplay among different factors, including dehydration of metamorphic facies, intraslab stress, preexisting structures, and abrupt changes in slab geometry, may explain the observed variations in seismogenesis of the DSZs and TSZs. The comparison of our relocated seismicity with the thermal model for the slab beneath Cook Inlet shows that the intermediate-depth earthquakes occur between 500°C and 900°C isotherms. The 2016 Mw 7.1 Iniskin earthquake and its aftershocks are located at ∼800°C–900°C. The intricate small-scale variations in different characteristics of the DSZs and intermediate-depth seismicity and their correlations with major geometrical and physical controls can provide insight into what governs the seismogenesis of subduction-induced earthquakes.

The Moho, slab and tomography of the East Japan forearc derived from seafloor S-net data

We investigate the 3-D P and S wave velocity structure of the East Japan forearc using a large number of high-quality arrival-time data of local earthquakes recorded by a permanent seafloor seismic network (S-net) and the land-based Kiban seismic network. Geometries of the subducting slab boundary and the Conrad and Moho discontinuities under the forearc are refined by referring to results of previous studies including active-source seismic soundings and seismicity located by the S-net. The obtained P and S wave velocity images are quite similar to each other, suggesting that they reflect robust features of the 3-D seismic structure. The subducting Pacific slab exhibits high velocities, whereas significant low-velocity (low-V) anomalies are revealed in the overlying Okhotsk plate beneath the forearc, which may reflect accretionary materials containing abundant fluids. The subducting oceanic crust atop the Pacific slab is revealed clearly as a thin low-V layer extending from the trench axis down to ~100 km depth under the volcanic front. Our results of earthquake relocation using the 3-D velocity model indicate that the double seismic zone in the subducting Pacific slab occurs from the trench axis down to ~180 km depth. The occurrence and focal mechanisms of the double seismic zone are caused by slab deformation during its subduction, such as slab bending or unbending, as well as other processes including hydration, dehydration, compositional variations, and phase changes in the subducting slab.

Metamorphism-facilitated faulting in deforming orthopyroxene: Implications for global intermediate-depth seismicity

Significance The exothermic metamorphic reaction in orthopyroxene (Opx), a major component of oceanic lithospheric mantle, is shown to trigger brittle failure in laboratory deformation experiments under conditions where garnet exsolution takes place. The reaction product is an extremely fine-grained material, forming narrow reaction zones that are mechanically weak, thereby facilitating macroscopic faulting. Oceanic subduction zones are characterized by two separate bands of seismicity, known as the double seismic zone. The upper band of seismicity, located in the oceanic crust, is well explained by dehydration-induced mechanical instability. Our newly discovered metamorphism-induced mechanical instability provides an alternative physical mechanism for earthquakes in the lower band of seismicity (located in the oceanic lithospheric mantle), with no requirement of hydration/dehydration processes.

Impact of bending-related faulting and oceanic-plate topography on slab hydration and intermediate-depth seismicity

Abstract It is commonly assumed that intermediate-depth seismicity is in some way linked to dehydration reactions inside subducting oceanic plates. Although there is growing evidence that the hydration state of an oceanic plate is controlled by its structure and degree of faulting, we do not have a quantitative understanding of this relationship. Double seismic zones offer the possibility of investigating changes in oceanic-plate hydration not only along strike but also with depth beneath the slab surface. To quantify the impact of oceanic-plate structure and faulting on slab hydration and intermediate-depth seismicity, with a focus on the genesis of double seismic zones, we correlate high-resolution earthquake catalogs and seafloor maps of ship-based bathymetry for the northern Chilean and Japan Trench subduction zones. The correlations show only a weak influence of oceanic-plate structure and faulting on seismicity on the upper plane of the double seismic zone, which may imply that hydration is limited by slow reaction kinetics at low temperatures 5–7 km below the seafloor and by the finite amount of exposed wall rock in the outer-rise region. These factors seem to limit hydration even if abundant water is available. Seismicity in the lower plane is, in contrast, substantially enhanced where deformation of the oceanic plate is high and distributed across intersecting faults. This likely leads to an increase in the volume of damaged wall rock around the faults, thereby promoting the circulation of water to mantle depths where serpentinization is faster due to elevated temperatures. Increased lower-plane seismicity around subducting oceanic features such as seamounts or fracture zones may also be caused by enhanced faulting around these features. Our results provide a possible explanation for the globally observed presence of rather homogeneous upper-plane seismicity in double seismic zones as well as for the commonly patchy and inhomogeneous distribution of lower-plane seismicity.

Complex Structure in the Nootka Fault Zone Revealed by Double‐Difference Tomography and a New Earthquake Catalog

AbstractWe employ an automatic earthquake detection algorithm to seismic waveforms recorded between the years 2000 and 2020 in southwest British Columbia. 32,121 events which possess at least three paired P‐ and S‐wave arrival times are located, compared to 21,538 seismic events in the existing Geologic Survey of Canada catalog. This augmented catalog is employed for double‐difference seismic tomography across Vancouver Island, with particular focus on the Nootka Fault zone (NFZ). The NFZ is a transform boundary that separates the Juan de Fuca and Explorer plates in a zone of distributed left‐lateral strike‐slip faulting. Tomographic results indicate that a double seismic zone exists within the NFZ that parallels Vp/Vs structure typically observed in subduction zones. Specifically, a dipping high Vp/Vs layer is underlain by reduced Vp/Vs material. Structural complexities are revealed by three Mw > 6 events and their aftershocks. The 2004 Mw 6.3 and 2011 Mw 6.4 events are interpreted to reside southeast of the NFZ within the Juan de Fuca plate at depths as shallow as 20 km, and the 2014 Mw 6.6 event is interpreted to reside within oceanic lithosphere of the NFZ at depths >35 km. Vp/Vs structure further indicates that underthrust oceanic lithosphere of the Explorer plate extends to 20 km depth beneath Brooks Peninsula. Our results support the hypothesis that the NFZ represents a structurally independent slice of oceanic lithosphere that exhibits (a) a more northerly trajectory than typically depicted, and (b) increased plate curvature when compared to either the Juan de Fuca or Explorer plates.

The thermal structure of subduction zones predicted by plate cooling models with variable thermal properties

SUMMARY Previous modelling studies have investigated the effects of experimentally constrained thermal properties (i.e. thermal conductivity, specific heat and the thermal expansion coefficient) on the thermal structure of subduction zones. However, these studies have not carefully considered whether the assumed thermal structure of the slab before subduction is consistent with geophysical observations. This study investigates the effects of thermal properties on the thermal structure of the Tohoku subduction zone, northeast Japan, by using the slab temperature at the trench determined from plate cooling models. Three types of thermal properties were tested: constant, temperature-dependent and temperature- and lithology-dependent types. For each case, the parameters for the plate cooling models were inferred based on the observed surface heat flow and seafloor depth using Bayes’ theorem. It was found that the predicted temperature and location of phase boundaries in the slab, which are possibly related to intermediate-depth earthquakes, are similar for the three cases. This suggests that, in the Tohoku subduction zone, constant thermal properties can be used in modelling to examine phenomena related to slab dehydration. The depth uncertainties for isotherms in the oceanic plate and slab increase with temperature, and are about ±10 and ±20 km for the 600 and 1200 °C isotherms, respectively. When this uncertainty is considered, the location of the serpentinite-out boundary matches that of the lower plane of double seismic zone, suggesting that dehydration may be important in triggering intermediate-depth seismicity. However, the large uncertainty makes it difficult to discuss in detail the origins of intraplate earthquakes, the lithosphere–asthenosphere boundary, and the lower boundary of the slab in terms of temperature.

TECTONIC STRESS FIELD AT INTERMEDIATE DEPTHS OF THE SOUTHERN FLANK OF THE KURIL-KAMCHATKA SEISMIC ZONE

The study has been done on the orientation of the principal axes of tectonic stress field at intermediate depths of earthquakes of the southern flank of the Kuril-Kamchatka subduction system separately for the upper and lower double seismic focal zones. Use has been made of the NIED and GlobalCMT catalogue data. The computation-based results are presented on schemes of the stressed state of the investigated areas and in tables. In the context of the southern Kuril Islands, evidence has been provided for predominance of the maximum compressive stresses along the slab plane in the upper layer and minimal compressive stresses (deviation extension) in the lower layer. However, the principal axes of maximum and minimum compression are displaced in direction relative to the slab dip: by 30–40° counterclockwise for the compression axis in the upper layer, which coincides with the direction of plate movement, and clockwise for the extension axis in the lower layer. This might be caused by the right-lateral strike-slip component of the Pacific Plate subduction. Unlike the general trend, the orientation of the principal axes of the stress field beneath the central Hokkaido-related segments in the upper layer is almost identical to that in the lower layer. There have also been found the segments exposed to shear stress, with the most extensive located opposite the northern Kunashir Island and beneath the southern Hokkaido Island. The results obtained for major large groups of clusters show good accordance with those published by other authors. The discrepancies relate primarily to small groups of isolated clusters showing local stress field heterogeneities.

Nested regional-global seismic tomography and precise earthquake relocation along the Hikurangi subduction zone, New Zealand

SUMMARY Seismic and geodetic examinations of the Hikurangi subduction zone (HSZ) indicate a remarkably diverse and complex system. Here, we investigate the 3-D P-wave velocity structure of the HSZ by applying an iterative, nested regional-global tomographic algorithm. The new model reveals enhanced details of seismic variations along the HSZ. We also relocate over 57 000 earthquakes using this newly developed 3-D model and then further improve the relative locations for 75 per cent of the seismicity using waveform cross-correlation. Double seismic zone characteristics, including occurrence, depth distribution and thickness change along the strike of the HSZ. An aseismic but fast Vp zone separates the upper and lower planes of seismicity in the southern and northern North Island. The upper plane of seismicity correlates with low Vp zones below the slab interface, indicating fluid-rich channels formed on top and/or within a dehydrated crust. A broad low Vp zone is resolved in the lower part of the subducting slab that could indicate hydrous mineral breakdown in the slab mantle. In the northern North Island and southern North Island, the lower plane of seismicity mostly correlates with the top of these low Vp zones. The comparison between the thermal model and the lower plane of seismicity in the northern North Island supports dehydration in the lower part of the slab. The mantle wedge of the Taupo volcanic zone (TVZ) is characterized by a low velocity zone underlying the volcanic front (fluid-driven partial melting), a fast velocity anomaly in the forearc mantle (a stagnant cold nose) and an underlying low velocity zone within the slab (fluids from dehydration). These arc-related anomalies are the strongest beneath the central TVZ with known extensive volcanism. The shallow seismicity (<40 km depth) correlates with geological terranes in the overlying plate. The aseismic impermeable terranes, such as the Rakaia terrane, may affect the fluid transport at the plate interface and seismicity in the overlying plate, which is consistent with previous studies. The deep slow slip events (25–60 km depths) mapped in the Kaimanawa, Manawatu and Kapiti regions coincide with low Vp anomalies. These new insights on the structure along the HSZ highlight the change in the locus of seismicity and dehydration at depth that is governed by significant variations in spatial and probably temporal attributes of subduction zone processes.

Along-strike variations in intermediate-depth seismicity and arc magmatism along the Alaska Peninsula

The Alaska Peninsula section of the Aleutian-Alaska subduction zone exhibits significant along-strike variations in plate coupling, earthquakes, and arc magmatism. Here we use regional and teleseismic data from 1990 to 2018 to investigate intermediate-depth intra-plate seismicity beneath the Alaska Peninsula. The distribution of these events shows five distinct segments from southwest to northeast: Unimak, Pavlof, Shumagin, Semidi, and Kodiak. We identify two double seismic zones (DSZs) in the Pavlof and Semidi segments but not in others. The magnitude-frequency relationship of these earthquakes (b-value) also varies along strike, suggesting more fluids in the southwestern segments and closer to the slab top interface. We also evaluate trace elements of arc lavas in the Alaska Peninsula. Th/La ratios indicate a significant amount of sediments subducted in the Semidi segment that includes Aniakchak and Veniaminof, and the sediment signature peaks beneath the Black Peak volcano. In contrast, high Pb/Th but low Th/La ratios from southwestern volcanoes suggest high fluids released from the subducted slab. In sum, our observations suggest that the slab is more hydrated and thus undergoing more intensive dehydration in the southwest than in the northeast, whereas more sediments are subducted in the northeast than in the southwest. These along-arc changes are generally in agreement with previous studies of plate coupling and the crustal structure of the incoming plate. Moreover, the complexities beyond the monotonic changes imply that the intermediate-depth intra-plate seismicity and sub-arc melting are influenced by not only the pre-existing fabrics in the incoming plate, but also the subducted sediments and plate stress state.

Seismic imaging of the double seismic zone in the subducting slab in Northern Chile

Double seismic zones are commonly observed in the subducting slabs in a global scale, serving as ideal examples for studying the seismogenetic mechanism of the intermediate-depth earthquakes. In this study, we relocate earthquakes and determine seismic velocity models using the double-difference seismic tomography method in the northern Chile subduction zone where a double seismic zone exists. The results suggest that the double seismic zone in northern Chile is located at about 50–140 ​km depth, with an interval of approximately 20 ​km between the two zones. The upper seismic zone is characterized by relatively low Vp (~7.8–8.0 ​km/s), low Vs (~4.4–4.5 ​km/s) and high Vp/Vs (~1.85) above the depth of ~90 ​km, while the region below ~90 ​km is distinguished by relatively high Vp (~8.2 ​km/s), high Vs (~4.8 ​km/s) and slightly high Vp/Vs (~1.75), which may be related to a series of dehydration reactions of hydrous minerals in the subducted oceanic crust. In comparison, the lower seismic zone is featured by the anomaly of low Vp/Vs (~1.7), although some local areas may consist of relatively high Vp/Vs values (~1.8), possibly due to the dehydration reaction of serpentine. Based on the Vp, Vs, Vp/Vs anomalies combined with previous petrological experiments and thermodynamic models, it can be derived that intermediate-depth earthquakes are mainly related to the dehydration of various hydrous minerals in the subducting slab. The dehydration process of hydrous minerals releases water into the subducting slab and subsequently leads to the increase of pore fluid pressure and the decrease of effective normal stress, thus causing the occurrence of brittle failure and intermediate-depth earthquakes in subduction zones. The imaging results of the northern Chile subduction zone further indicate that the existence of the double seismic zone is related to the dehydration process of different hydrous minerals.

Hydration of the crust and upper mantle of the Hikurangi Plateau as it subducts at the southern Hikurangi margin

Controlled-source seismic studies at most subduction zones show that bending of the subducting plate results in reduced seismic wave-speeds in the crust and upper mantle near the trench. Similar studies also have found unusually high P-wave-speeds (Vp) in the upper mantle under oceanic plateaus. Onshore-offshore seismic profiling at the southern Hikurangi margin, where the ≈120 Ma old oceanic Hikurangi Plateau is subducting, indicates that a fast (Vp≈8.7±0.2 km/s) upper mantle layer lies beneath a ≈25 km thick mantle layer with more regular wave-speeds (Vp≈8.0±0.2 km/s) under the Hikurangi trough. This is consistent with previous findings of upper mantle Vp≈8.7-9.0 km/s in the margin-parallel direction under the North Island (≈100 km northwest of the deformation front) at depths ≈8-10 km below the Moho. Our profiles are margin-perpendicular, thus we show that the upper mantle lid of the subducting Pacific Plate is characterized by unusually high P-wave-speeds along all azimuths. We find an area of lowered Vp in the ≈12±1 km thick Hikurangi Plateau crust beneath the trough. This drop in Vp is ≈10%, and a similar drop in Vp is deduced to depths of 25±2 km into the upper mantle. We interpret that the increase in thickness of the regular mantle beneath the trough results from the formation of a low-velocity zone in the faster upper mantle layer; this zone formed from serpentinisation by hydration through bending-induced normal faults and/or due to crack porosity introduced by thermal cracking, further enhanced by bending-related faulting. Thus the “regular mantle” (Vp≈8 km/s) is not in fact regular, but rather the high-speed mantle has mechanically bent, fractured, and altered. The absolute depth of fast mantle Vp under the Hikurangi trough is around 50 km. The onset of the lower band of seismicity of the double seismic zone and high upper mantle Vp under the North Island is observed at similar depths. This is consistent with the hypothesis that the lower band of earthquakes in a double seismic zone is due to antigorite dehydration processes, a hydrous mineral formed in the low velocity zone in the upper mantle beneath the trough. Our study on the Hikurangi margin is different, as the subducting plate here contains a ≈120 Ma old oceanic plateau with a ≈12 km thick crust, but the results are similar to other subduction margins where regular oceanic crust is subducting.