5:15 PM - 6:30 PM
[MTT27-P04] Joint estimation of parameters and initial conditions of numerical models of fault slip using data assimilation methods
Keywords: data assimilation, fault slip, friction parameters
In this study, we employ a two-dimensional planar fault in a three-dimensional elastic half-space. The evolution of slip rate and shear stress on the fault is governed by the force balance equation and rate- and state-dependent friction law. We assume that time-series of displacements at stations on the ground surface are observable. We employ the ensemble Kalman filter and smoother (EnKF) and adjoint method to simultaneously estimate cumulative slip, slip rate, and frictional state variable for each subdivided fault patch and friction parameters.
We conduct synthetic tests to assess the validity of the proposed methods. In these synthetic tests, we focus on afterslip, which is transient, decelerating, aseismic slip triggered by stress changes due to a large earthquake. We impose an earthquake using a circular crack with constant stress drop and calculate coseismic shear stress changes on the fault. Slip rate and frictional state variables immediately after the earthquake are determined based on the calculated coseismic shear stress changes. Using these slip rate and frictional state variables as the initial condition, we employ the rate-state friction model to simulate the evolution of afterslip assuming that the fault has a velocity-strengthening frictional property. We then compute displacement time series at stations on the ground surface using the simulated slip history. Finally, we add normally distributed observation errors to generate simulated data. This simulated data set is used in the synthetic tests.
For the case of the EnKF, we find that the prediction step of the EnKF becomes numerically unstable if some of the initial ensemble members have values significantly different from the true values. These numerical instabilities occur when, for example, the initial ensemble members are randomly generated. In contrast, if the initial ensemble members have similar values to the true values, such numerical instabilities do not occur and consequently the EnKF can successfully be implemented. These results indicate that it is important to develop a method to generate appropriate initial ensemble to successfully implement the EnKF for our model. We will also show the results obtained from the adjoint method in the presentation.