ABSTRACT
Rapid seismic intensity maps for damaging earthquakes enable the swift implementation of earthquake disaster mitigation action, issuance of accurate tsunami warnings, and prevention of associated secondary disasters. However, many countries lack dense local seismic observation networks, making it infeasible to obtain accurate seismic intensity maps of earthquakes within a few hours, particularly for earthquakes that have considerable source extents. In this study, we developed a new algorithm for rapidly obtaining seismic intensity maps of damaging earthquakes. With our model, source energy radiation is acquired using backprojection, and then the locations and relative amplitudes of the fault geometry and subevents are determined. Peak ground accelerations and peak ground velocities (PGVs) are subsequently calculated based on ground‐motion prediction equations and the distribution of the estimated subevents. PGVs are then further site‐corrected using the database (Wald and Allen, 2007; Heath et al., 2020). The algorithm was applied to the 2008 7.9 Wenchuan and 2010 6.9 Yushu earthquakes, and the resulting seismic intensity maps were highly similar to those generated by field surveys. The algorithm is simple and straightforward to use, and local real‐time instrument observations are not required. Calculations can be performed automatically, and reliable seismic intensity maps can be issued within 30 min following damaging earthquakes. The model’s application may assist greatly with rescue and recovery efforts, and enable tsunami hazards to be evaluated immediately following earthquakes, particularly in regions lacking dense observation networks.