ShakeMap Ground-Motion Interpolation Using the Multivariate Normal Distribution

Abstract:

A key update in the next major release of ShakeMap is the application of the multivariate normal distribution (MVN) to ground motion interpolation. The MVN provides a complete solution to the problem of ground-motion interpolation, fully accounting for the correlation among the instrumental observations both spatially and across intensity measure types (IMTs), and is therefore an important improvement over the current ShakeMap interpolation scheme. Practical application of the MVN in a ShakeMap setting, however, is challenging for several reasons: 1) the large size of typical ShakeMap output grid makes computation of the complete MVN impractical on commonly available computers; 2) a large number of IMTs in the input may also inflate the required matrices to impractical dimensions; 3) the lack of currently available empirical cross-correlation functions between some ShakeMap IMTs [e.g., between peak ground velocity (PGV) and Modified Mercalli Intensity (MMI)] prevents a straightforward implementation of the MVN in some cases; 4) the computation of the event-specific bias is complicated by the presence of multiple IMTs in the input and heteroskedastic inter-event residuals in some ground motion models. Here, we describe a number of adaptations and simplifying assumptions made to overcome these, and other, obstacles. For instance, we discuss a piecewise application of the MVN that produces output grids of amplitude and variance, without the necessity of computing the full covariance matrix of the output grid. We also show that the down-weighting of uncertain data as described by Worden et al. (BSSA, in press) allows for the inclusion of data converted from one IMT to another, while accounting for the uncertainty inherent in the conversion. This approach enables, for example, the inclusion of PGV data in a map of MMI despite the lack of an empirical cross-correlation function, by first applying an appropriate conversion of PGV to MMI.

Slidecast:

https://vimeo.com/276918642

Remote Dynamic Triggering in Southeast Africa

Abstract:

We present results of the first systematic search for evidence of dynamic triggering of microearthquakes in southeast Africa. We analyzed seismic waveforms from 53 broadband and short-period stations in SE Africa, to identify increases in the rate of microearthquakes and tremor coincident with the passage of the Mw8.6 2012 Indian Ocean earthquake teleseismic waves (Neves et al., 2018). We computed β-values to evaluate the statistical significance of the observed increase in seismic activity, and statistically significant dynamic triggering at 7 stations. Two of these stations are located in the northeast coast of Madagascar. The other 5 stations are located in the Kaapvaal Craton, southern Africa, around one of the most active clusters of seismicity in the region that has been associated with a combination of natural and anthropogenic factors. We found no evidence of dynamically triggered microearthquakes in stations located near the structures of the East African Rift system (EARS). We also did not find any evidence of deep tectonic tremor in the entire region. Although hydrothermal activity has been identified close to the stations that recorded dynamic triggering, it also exists near the EARS structures where no triggering was observed. Our results suggest that factors other than solely tectonic regime and geothermalism are needed to explain the mechanisms underlying earthquake triggering. Further investigation revealed a higher vertical to horizontal peak-ground velocity ratio in the stations where triggering occurred, suggesting a local distortion of the wavefield, possibly due to topography, structure or site effects, that may facilitate triggering. We are in the process of analyzing additional mainshocks located at the African Continent to better constrain the factors that favor dynamic triggering in the region and study the effects of dynamic stresses in the intermediate-distant ranges. Updated results will be presented in the meeting.

Slidecast:

https://vimeo.com/276918506

Strong Motion Simulation for the 2004 Southwest-Off Kii Peninsula Earthquake with Pseudo Point-Source Model

Abstract:

We have been testing the applicability of a new simplified source model called the “pseudo point-source model,” which was aimed to calculate strong ground motions (Nozu 2012). The model is simpler, and involves much less model parameters than the conventional characterized source model, which itself is a simplified expression of an actual earthquake source. In the model, the spatio-temporal distribution of slip within a subevent is not modeled. Instead, the source spectrum associated with the rupture of a subevent is modeled and it is assumed to follow the omega-square model. By multiplying the source spectrum with the path effect and the site amplification factor, the Fourier amplitude at a target site can be obtained. Then, combining it with Fourier phase characteristics of a smaller event, the time history of strong ground motions from the subevent can be calculated. Finally, by summing up contributions from the subevents, strong ground motions from the entire rupture can be obtained. In this study, the model was applied to the main shock of the 2004 southwest-off Kii Peninsula earthquake (Mw7.4), which occurred on September 5, 2004, 23:57 (JST) near the Nankai Trough, within the subducting Philippine Sea plate (e.g., Park and Mori 2005; Bai et al. 2007). Because of its location, strong motion records obtained during the earthquakes are suitable for the validation of strong motion simulation techniques. The pseudo point-source model for the earthquake was developed by appropriately locating point sources on the fault plane, based on the results of the waveform inversion. Then the model was used to simulate strong ground motions in a wide region. It was found that, when the model is used with an appropriate path model, it can accurately simulate strong ground motions at remote stations, even when the stations are located on deep sedimentary basins such as the Osaka, Nobi and Kanto plains.

Slidecast:

https://vimeo.com/276922356

Kinematic Source Models for Earthquake Simulations with Fault-Zone Plasticity

Abstract:

Fault slip and surface deformation patterns are essential factors in seismic hazard assessment. However, slip inversions reveal a widely observed shallow slip deficit (SSD) which has not yet been clearly explained. One possible cause of the SSD is distributed plastic deformation in the fault damage zone near the surface. Roten et al. (2017) performed 3D dynamic rupture simulations of the 1992 M7.3 Landers earthquake in a medium governed by Drucker-Prager plasticity. The study showed that while linear simulations tend to underpredict SSD and off-fault deformation (OFD), nonlinear simulations with moderately fractured rock mass can properly reproduce results that are consistent with the 30-60% SSD and around 46% OFD reported in geodetic observations. Analysis of the Roten et al. (2017) results show that discrepancies between linear and nonlinear simulations are only significant in the top hundreds of meters of the surface-rupturing fault. Although inelastic responses in the fault damage zone provide more realistic representations of earthquake physics, it can be computationally expensive or numerically unfeasible (e.g., in adjoint methods) to include nonlinear effects in ground motion simulations. One solution proposed here is to use an equivalent kinematic source (EKS) model that mimics the fault-zone plastic effects. This method generates source-time-functions by modifying the slip rate time histories based on comparisons of linear and non-linear dynamic rupture models, which are then used as input to kinematic simulations. The EKS model is shown to capture the SSD and OFD patterns observed in dynamic simulations with fault-zone plasticity. Further verification of the method is needed before the anticipated use in practical applications such as the Southern California Earthquake Center (SCEC) CyberShake and Broadband platforms.

Slidecast:

https://vimeo.com/276922362

A Computational Model of the Seismic Response in the Valley of Mexico after the Damaging m7.1 Earthquake of September 19, 2017

Abstract:

Built-up on top of ancient lake deposits, Mexico City experiences some of the largest seismic site effects in the world. The M7.1 intermediate-depth earthquake of September 19, 2017 (S19) collapsed 43 one-to-ten story buildings in the city close to the western edge of the lake-bed sediments, on top of the geotechnically-known transition zone. In this work we explore the physical reasons explaining such a damaging pattern and the long-lasting strong motion records well documented from past events by means of new observations and high performance computational modeling. Besides the extreme amplification of seismic waves, duration of intense ground motion in the lake-bed lasts more than three times those recorded in hard-rock a few kilometers away. Different mechanisms contribute to the long lasting motions, such as the regional dispersion and multiple-scattering of the incoming wavefield all the way from the source. However, recent beamforming observations at hard-rock suggest that duration of the incoming field is significantly shorter than the strong shaking in the lake-bed. We show that despite the highly dissipative shallow deposits, seismic energy can propagate long distances in the deep structure of the valley, promoting also a large elongation of motion. Our 3D simulations reveal that the seismic response of the basin is dominated by surface-waves overtones, and that this mechanism increases the duration of ground motion up to 280% and 500% of the incoming wavefield duration at 0.5 and 0.3 Hz, respectively. Furthermore, our results indicate that the damage pattern of the S19 earthquake is most likely due to the propagation of the fundamental mode in the transition zone of the basin generated as a consequence of a localized basin edge effect. Some conclusions contradicts what has been previously stated from observational and modeling investigations, where the basin itself has been discarded as a preponderant factor promoting long and devastating shaking in Mexico City. Reference: Cruz-Atienza, V. M., J. Tago, J. D. Sanabria-Gómez, E. Chaljub, V. Etienne, J. Virieux and L. Quintanar. Long Duration of Ground Motion in the Paradigmatic Valley of Mexico. Nature – Scientific Reports, 6, 38807; doi:10.1038/srep38807, 2016.

Slidecast:

https://vimeo.com/276921577

Real-Time Earthquake Impact Mapping in Indonesia: The Case of the 2016 Aceh Earthquake (m6.5)

Abstract:

Indonesia is one of the most earthquake disaster-prone countries in the world due to its tectonic activity and high population exposure. To assist with the earthquake response process, Indonesia’s Agency of Meteorology, Climatology and Geophysics (BMKG) generates near-real-time maps of ground motion and shaking intensity (i.e. ShakeMaps) following significant earthquakes. The original ShakeMaps are based on empirical predictions using basic hypocentre and magnitude information, and are generally generated within 5 minutes of the earthquake’s origin. Such maps are updated automatically by incorporating recorded ground motion data from Indonesia’s national Strong-Motion Network. This network, operated and maintained by BMKG, currently includes 278 strong-motion stations, and in 2017 has recorded more than 5000 records from 620 events. The ShakeMap products are subsequently used by Indonesian National Board for Disaster Management (BNPB) to produce near-real time earthquake impact maps in terms of the population exposed to different levels of ground shaking as well as expected number of fatalities. In this paper, we describe the implementation of real-time earthquake impact alerting systems within BMKG and BNPB. We also demonstrate the use of these systems in response to the recent M­W 6.5 Kota-Banda Aceh earthquake that occurred on the 7th of December 2016 at a depth of 15 km. The earthquake was reported to have caused 104 fatalities and ~8000 number of displaced people. The integrated ShakeMap system automatically generated shaking estimates calibrated by BMKG’s strong-motion network within 10 minutes of the event’s origin time. The BMKG Shakemaps are automatically uploaded to the InaSAFE Realtime platform, managed by BNPB, to estimate the number of people exposed. It has been observed that recorded ground motions as well as number of casualties are generally consistent with theoretical models.

Slidecast:

https://vimeo.com/276918887