Physics-Informed Deep Learning for Estimating the Spatial Distribution of Frictional Parameters in Slow Slip Regions
Abstract
Slow slip events (SSEs) have been observed in many subduction zones and are understood to result from frictional unstable slip on the plate interface. The diversity of their characteristics and the fact that interplate slip can also be seismic suggest that frictional properties are heterogeneous. We are however lacking methods to determine spatial variations of frictional properties. In this paper, we employ a Physics‐Informed Neural Network (PINN) to achieve this goal using a synthetic model inspired by the long‐term SSEs observed in the Bungo channel. PINN is a deep learning technique that can be used to solve the differential equations representing the physics of the problem and determine the model parameters from observations. We start with an idealized case where it is assumed that fault slip is directly observed. We next move to a more realistic case where the observations consist of synthetic surface displacement velocity data measured by virtual GNSS stations. We find that the geometry and friction properties of the velocity weakening region, where the slip instability develops, are well estimated, especially if surface displacement velocity above the velocity weakening region is observed. Our PINN‐based method can be seen as an inversion technique with the regularization constraint that fault slip obeys a particular friction law. This approach remediates the issue that standard regularization techniques are based on non‐physical constraints. Our results show that the PINN‐based method is a promising approach for estimating the spatial distribution of friction parameters from GNSS observations.
Copyright and License
© 2025. The Author(s). This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.
Acknowledgement
This work was supported by the MEXT Project for Seismology toward Research Innovation with Data of Earthquake (STAR-E) [Grant JPJ010217] and by the JSPS KAKENHI [Grants 23K03552, 24K02951, 24H01019, 23H00466, 21H05203, and 21K03694]. The travel of Rikuto Fukushima to the California Institute of Technology was supported by the Summer Undergraduate Research Fellowship (SURF) program of the California Institute of Technology.
Data Availability
Python code developed in this study is available via https://doi.org/10.5281/zenodo.14977672 (Fukushima, 2025).
Supplemental Material
Files
Name | Size | Download all |
---|---|---|
md5:8422d46684769b6f802ef24b16c47cba
|
4.0 MB | Preview Download |
md5:6ef8ec62e7fc7a736acaed7d4ccb15f3
|
906.3 kB | Preview Download |
Additional details
- Ministry of Education, Culture, Sports, Science and Technology
- JPJ010217
- Japan Society for the Promotion of Science
- 23K03552
- Japan Society for the Promotion of Science
- 24K02951
- Japan Society for the Promotion of Science
- 24H01019
- Japan Society for the Promotion of Science
- 23H00466
- Japan Society for the Promotion of Science
- 21H05203
- Japan Society for the Promotion of Science
- 21K03694
- California Institute of Technology
- Summer Undergraduate Research Fellowship (SURF) -
- Accepted
-
2025-04-29
- Available
-
2025-05-09Version of record online
- Available
-
2025-05-09Issue online
- Caltech groups
- Center for Geomechanics and Mitigation of Geohazards (GMG), Seismological Laboratory, Division of Geological and Planetary Sciences (GPS)
- Publication Status
- Published