3D finite-difference simulations of strong ground motions in the Yanhuai area, Beijing, China during the 1720 Shacheng earthquake (Ms 7.0) using a stochastic finite-fault model

Three-dimensional finite-difference (3D-FD) simulations of earthquake wave propagations in the Yanhuai area were performed for the 1720 Shacheng earthquake (Ms 7.0) using a stochastic finite-fault model, running on a parallel supercomputer Hitachi-SR8000 (http://www.lrz-muenchen.de). A stochastic finite-fault model was implemented into the 3D-FD program. The basic idea of the stochastic finite-fault model is that the fault plane can be subdivided into several subfaults (or elements, sources). Radiation from a large earthquake is the sum of contributions from all subfaults with proper time delays, each of which acts as a small independent double-couple point source. Heterogeneity of the fault rupture process was modeled by randomizing the location of initial rupture (hypocenter), slip vectors (slip, rake), and rise-times of subfaults in this study. A 3D velocity model of the Yanhuai area was constructed based on studies that analyzed available geological and geophysical information. A grid increment of 75 m in three directions was used in the 3D-FD simulation, which made it possible to capture the short period information with a resolution as low as 0.5 s in sediment regions. The uncertainties of simulated results caused by the stochastic finite-fault model were studied with a homogeneous 3D model. We found that the effects of the randomness of source on simulated ground motions are only limited in near-fault-region including the surface exposure of the fault and its vicinities, which occupies about 5% of the whole study area. This article presents an integrated approach for simulating the strong ground motions for engineering purpose using the 3D-FD method. Such simulations would be useful for hazard mapping, land using planning, insurance rate assessment, particularly in planning, preparedness, and coordinating emergency response, which must be based on realistic situations induced by concrete (historic or scenario) earthquakes.