| Popis: |
Physics-based simulations of earthquake ground motion are useful to complement recorded ground motions. However, the computational expense of performing numerical simulations hinders their applicability to tasks that require real-time solutions or ensembles of solutions for different earthquake sources. To enable rapid physics-based solutions, we present a reduced-order modeling (ROM) approach based on interpolated proper orthogonal decomposition (POD) to predict peak ground velocities (PGVs). As a demonstrator, we consider PGVs from regional 3D wave propagation simulations at the location of the 2008 Mw 5.4 Chino Hills earthquake using double-couple sources with varying depth and focal mechanisms. These simulations resolve frequencies $\leq$ 1.0 Hz and include topography, viscoelastic attenuation, and S-wave speeds $\geq$ 500 m/s. We evaluate the accuracy of the interpolated POD ROM as a function of the approximation method. Comparing the radial basis function (RBF), multilayer perceptron neural network, random forest, and $k$-nearest neighbor, we find that the RBF interpolation gives the lowest error ($\approx$ 0.1 cm/s) when tested against an independent dataset. We also find that evaluating the ROM is $10^7-10^8$ times faster than the wave propagation simulations. We use the ROM to generate PGV maps for one million different focal mechanisms, in which we identify potentially damaging ground motions and quantify correlations between focal mechanism, depth, and accuracy of the predicted PGV. Our results demonstrate that the ROM can rapidly and accurately approximate the PGV from wave propagation simulations with variable source properties, topography, and complex subsurface structure. |
