PINN-based short-term forecasting of fault slip evolution during the 2010 slow slip event in the Bungo Channel, Japan

TL;DR

This paper introduces a physics-informed neural network framework for short-term fault slip forecasting during slow slip events, effectively integrating geodetic data and fault mechanics to improve stability and accuracy.

Contribution

It presents a novel PINN-based data assimilation method that incorporates fault heterogeneity for stable, accurate slip evolution predictions during slow slip events.

Findings

Successfully forecasts slip evolution using limited data

Heterogeneous friction models outperform homogeneous ones

Identifies fault regions critical for slip nucleation and propagation

Abstract

Monitoring and forecasting fault slip evolution are fundamental for understanding earthquake cycles and assessing future seismic hazards. This study proposes a physics-based data assimilation framework that integrates geodetic observations with fault mechanics introducing spatial heterogeneity in frictional properties, with a particular focus on short-term fault slip forecasting. The proposed method employs physics-informed neural networks (PINNs) to calculate fault slip evolutions and to optimize the spatial distribution of frictional properties and is applied to the 2010 slow slip event beneath the Bungo Channel, southwest Japan, by changing the data period to be assimilated. When only the initial phase of slip acceleration is assimilated, a velocity-weakening frictional region is inferred beneath southwest Shikoku, corresponding to the initial nucleation are of the slow slip event. Out results demonstrate that the PINN-based data assimilation framework successfully forecasts slow transient slip even when only slip acceleration data are assimilated, whereas forecasts based on frictionally homogeneous models result in unstable fast slip. This difference can be interpreted as a consequence of introducing frictional heterogeneity, which allows both the characteristic size of the slipping region and the critical nucleation size to be variable, leading to stable slip evolution consistent with observations. When longer observation periods are assimilated, a velocity-strengthening region emerges around the slip-weakening patch, progressively restricting the direction of slip propagation. This velocity-strengthening region is interpreted as a mechanical constraint imposed by fault physics, linking the slip regions required to reproduce the observed geodetic time series. The results highlight the capability of PINN-based data assimilation incorporating geodetic observations and fault mechanics.