Main content

Home

Menu

Loading wiki pages...

View
Wiki Version:
Numerical data of results shown in *Source time functions of earthquakes based on a stochastic differential equation*: https://doi.org/10.1038/s41598-022-07873-2 See the above paper for details of case A and case B. For each cases, the Bessel processes ($X_t^{(1)}$ and $X_t^{(2)}$) and their convolution (STF) are given in ascii format. $X_t^{(1)}$ has shorter duration than $X_t^{(2)}$. $n$-th STF (`STF/${n}.txt`) is a convolution of $n$-th $X_t^{(1)}$ (`Xt1/${n}.txt`) and $n$-th $X_t^{(2)}$ (`Xt2/${n}.txt`), where $n=1, ..., 1000$. STFs are not normalized by their duration or total moment, so the normalization is required to obtain plots in the paper. The Bessel process as a 1-D array is obtained by the following Julia-1.6.1 function using DifferentialEquations.jl package, function BesselProcess(Tmin,Tmax) Xₜ = zeros(Tmax) (d, σ, u₀, dt) = (0, 1, 1e-3, 1e-6) while true f(u,p,t) = ifelse(u≤0, 0, 0.5*(d-1)/u) g(u,p,t) = ifelse(u≤0, 0, σ) prob = SDEProblem(f,g,u₀,(0, Tmax*dt)) sol = solve(prob,SRIW1(),dt=dt,adaptive=false) Xₜ = max.(0, sol.u) if Xₜ[Tmin-1] > 0 && Xₜ[Tmax] == 0 break end end return Xₜ end where `(Tmin, Tmax) = (1000, 2000)` and `(Tmin, Tmax) = (200, 2000)` for case A and B, respectively, and `d=0` is equivalent to $b=1$ for GR law (see eq.(5) in the paper). Note that $d$ must be close to $0$, otherwise, `break` may not work because the duration $T$ does not satisfy $T_\textrm{min} < T < T_\textrm{max}$ probabilistically (see eq.(3) in the paper).