When a major earthquake strikes, the resulting devastation can be compounded or even exceeded by the subsequent cascade of triggered seismicity. As the Nepalese recover from the 25 April 2015 shock, knowledge of what comes next is essential. We calculate the redistribution of crustal stresses and implied earthquake probabilities for different periods from daily to 30 years into the future. An initial forecast was completed before a M=7.3 earthquake struck on May 12, 2015 that enables a preliminary assessment; post-forecast seismicity has so far occurred within a zone of 5-fold probability gain. Evaluation of the forecast performance, using 2 months of seismic data, reveals that stress-based approaches present improved skill in higher magnitude triggered seismicity. Our results suggest that considering the total stress field, rather than only the co-seismic one, improves the spatial performance of the model based on the estimation of a wide range of potential triggered faults following a mainshock.
2
IntroductionThe Himalayas rise by absorbing 18-22 mm/yr of Indian-Eurasian plate convergence (Ader et al., 2012) (Fig. 1). On the 25 th of April, 2015, a M=7.8 earthquake ruptured the low-angle (10˚) fault contact between the two plates, leading to more than 8000 confirmed fatalities, 18,000 people injured, and a million families affected (Nepal Red Cross Society 1 ). Eight UNESCO World Heritage sites were damaged or destroyed. The 25 April shock struck the eastern edge of a 500-km-wide gap between historic earthquakes along the Himalayan front (Fig. 1). The potential of triggered earthquakes beneath highly populated basins in the central Himalaya could compound this catastrophe akin to tragedies in Turkey , China , and New Zealand . Recent studies, based on retrospective experiments, investigate the predictive power of short-term earthquake forecasts within different distance ranges and periods (Segou et al., 2013;Strader and Jackson, 2015), but prospective forecasts still face challenges such as the quality of real-time data, the availability of credible historical/modern earthquake catalogs ,and very short research windows.In this study, we calculate the expected redistribution of stress in the Himalayan crust and develop a method, especially applicable to frontier regions, to determine the probability of triggered earthquakes in both space and time. We make prospective earthquake forecasts using four methods incorporating physics-based and statistical approaches for varying time horizons, and formally evaluate the shortest-term calculations.The results reveal the efficiency of physics-based forecasts in estimating short-term earthquake probabilities and their critical contribution in modeling off-fault triggered events, highlighted by the M=5.6 Xegar (3 hours after the mainshock) and the M=7.3 Kodari (May 12) earthquakes.
Development of Short-Term Earthquake ForecastsShort-term earthquake forecasts can use empirical statistics to anticipate cascades of triggered events, or simulations of stress redistribution fol...