Physical mechanisms of oceanic mantle earthquakes: Comparison of natural and experimental events

Different earthquake nucleation conditions revealed by stress drop and b-value mapping in the northern Chilean subduction zone

Stress drop is an earthquake property indicative for the characteristic relation of slip to fault dimension. It is furthermore affected by fault strength, fault topography, the presence of fluids, rupture size, slip, and velocity. In this article, the stress drop image of an entire subduction zone, namely for the seismically highly active northernmost part of Chile, is combined with mapped b-values and their corresponding magnitude distribution in order to better constrain the conditions under which earthquakes of different provenances may nucleate. The underlying recent earthquake catalog contains over 180,000 events, covering 15 years of seismicity, from which more than 50,000 stress drop estimates were computed. Their spatial average segments the subduction zone into different parts, i.e., average stress drop between seismotectonic areas is different, although this difference is small compared to the natural scatter of stress drop values. By considering stress drop variations, b-value map, magnitude distribution, and thermal models, candidate earthquake nucleation mechanisms are identified which can explain the observed distributions. This is done for two exemplary regions: (1) The plate interface, where principally lower stress drop events are found, while at the same time a high spatial heterogeneity of stress drop values is observed. This indicates relatively smooth or lubricated rupture surfaces, and locally it suggests the existence of alternating regions controlled by strong asperities, weaker material, or creep. (2) The highly active intermediate depth (ID) seismicity region, where the variation of stress drop and b-value point to a gradual change of nucleation mechanism from dehydration embrittlement at the top of the ID cloud, over dehydration driven stress transfer in its central part, to thermal runaway shear mechanisms at its bottom. In both cases, the combination of stress drop and b-value distribution helps to better understand the origin and the differences of the observed seismicity.

A switch from horizontal compression to vertical extension in the Vrancea slab explained by the volume reduction of serpentine dehydration

The Vrancea slab, Romania, is a subducted remnant of the Tethyan lithosphere characterized by a significant intermediate-depth seismicity (60-170 km). A recent study showed a correlation between this seismicity and major dehydration reactions, involving serpentine minerals up to 130 km depth, and high-pressure hydrated talc deeper. Here we investigate the potential link between the triggering mechanisms and the retrieved focal mechanisms of 940 earthquakes, which allows interpreting the depth distribution of the stress field. We observe a switch from horizontal compression to vertical extension between 100 and 130 km depth, where the Clapeyron slope of serpentine dehydration is negative. The negative volume change within dehydrating serpentinized faults, expected mostly sub-horizontal in the verticalized slab, could well explain the vertical extension recorded by the intermediate-depth seismicity. This apparent slab pull is accompanied with a rotation of the main compressive stress, which could favour slab detachments in active subduction zones.

Dehydration-induced earthquakes identified in a subducted oceanic slab beneath Vrancea, Romania

Vrancea, Eastern Romania, presents a significant intermediate-depth seismicity, between 60 and 170 km depth, i.e. pressures from 2 to 6.5 GPa. A debate has been lasting for decades regarding the nature of the seismic volume, which could correspond to the remnant of a subducted slab of Tethyan lithosphere or a delamination of the Carpathians lithosphere. Here we compile the entire seismicity dataset (≈ 10,000 events with 2 ≤ Mw ≤ 7.9) beneath Vrancea for P > 0.55 GPa (> 20 km) since 1940 and estimate the pressure and temperature associated with each hypocenter. We infer the pressure and temperature, respectively, from a depth-pressure conversion and from the most recent tomography-based thermal model. Pressure–temperature diagrams show to what extent these hypocentral conditions match the thermodynamic stability limits for minerals typical of the uppermost mantle, oceanic crust and lower continental crust. The stability limits of lawsonite, chloritoid, serpentine and talc minerals show particularly good correlations. Overall, the destabilization of both mantle and crustal minerals could participate in explaining the observed seismicity, but mantle minerals appear more likely with more convincing correlations. Most hypocentral conditions match relatively well antigorite dehydration between 2 and 4.5 GPa; at higher pressures, the dehydration of the 10-Å phase provides the best fit. We demonstrate that the Vrancea intermediate-depth seismicity is evidence of the current dehydration of an oceanic slab beneath Romania. Our results are consistent with a recent rollback of a W-dipping oceanic slab, whose current location is explained by limited delamination of the continental Moesian lithosphere between the Tethyan suture zone and Vrancea.