Unsupervised machine learning reveals slab hydration variations from deep earthquake distributions beneath the northwest Pacific

Although transformational faulting in the rim of the metastable olivine wedge is hypothesized as a triggering mechanism of deep-focus earthquakes, there is no direct evidence of such rim. Variations of the b value – slope of the Gutenberg-Richter distribution – have been used to decipher triggering and rupture mechanisms of deep earthquakes. However, detection limits prevent full understanding of these mechanisms. Using the Japan Meteorological Agency catalog, we estimate b values of deep earthquakes in the northwestern Pacific Plate, clustered in four regions with unsupervised machine learning. The b-value analysis of Honshu and Izu deep seismicity reveals a kink at magnitude 3.7–3.8, where the b value abruptly changes from 1.4–1.7 to 0.6–0.7. The anomalously high b values for small earthquakes highlight enhanced transformational faulting, likely catalyzed by deep hydrous defects coinciding with the unstable rim of the metastable olivine wedge, the thickness of which we estimate at ∼1 km