Earthquake Cycle Modelling of Multi-segmented Faults: Dynamic Rupture and Ground Motion Simulation of the 1992 M_w 7.3 Landers Earthquake

We perform earthquake cycle simulations with the goal of studying the characteristics of source scaling relations and strong ground motions in multi-segmented fault ruptures. The 1992 M_w 7.3 Landers earthquake is chosen as a target earthquake to validate our methodology. The model includes the fault geometry for the three-segmented Landers rupture from the SCEC community fault model, extended at both ends to a total length of 200 km, and limited to a depth to 15 km. We assume the faults are governed by rate-and-state (RS) friction, with a heterogeneous, correlated spatial distribution of characteristic weakening distance Dc. Multiple earthquake cycles on this non-planar fault system are modeled with a quasi-dynamic solver based on the boundary element method, substantially accelerated by implementing a hierarchical-matrix method. The resulting seismic ruptures are recomputed using a fully-dynamic solver based on the spectral element method, with the same RS friction law. The simulated earthquakes nucleate on different sections of the fault, and include events similar to the M_w 7.3 Landers earthquake. We obtain slip velocity functions, rupture times and magnitudes that can be compared to seismological observations. The simulated ground motions are validated by comparison of simulated and recorded response spectra

Earthquake Cycle Modelling of Multi-segmented Faults: Dynamic Rupture and Ground Motion Simulation of the 1992 M_w 7.3 Landers Earthquake

We perform earthquake cycle simulations with the goal of studying the characteristics of source scaling relations and strong ground motions in multi-segmented fault ruptures. The 1992 M_w 7.3 Landers earthquake is chosen as a target earthquake to validate our methodology. The model includes the fault geometry for the three-segmented Landers rupture from the SCEC community fault model, extended at both ends to a total length of 200 km, and limited to a depth to 15 km. We assume the faults are governed by rate-and-state (RS) friction, with a heterogeneous, correlated spatial distribution of characteristic weakening distance Dc. Multiple earthquake cycles on this non-planar fault system are modeled with a quasi-dynamic solver based on the boundary element method, substantially accelerated by implementing a hierarchical-matrix method. The resulting seismic ruptures are recomputed using a fully-dynamic solver based on the spectral element method, with the same RS friction law. The simulated earthquakes nucleate on different sections of the fault, and include events similar to the M_w 7.3 Landers earthquake. We obtain slip velocity functions, rupture times and magnitudes that can be compared to seismological observations. The simulated ground motions are validated by comparison of simulated and recorded response spectra

Validation of the Osaka basin velocity model by waveform simulation using sources inverted with reciprocal Green’s functions

Abstract To advance the methodology for validating velocity models by waveform comparison, we estimated source parameters for small magnitude earthquakes that can be approximated by point sources. Instead of using published source models, we used the reciprocity method to calculate 3D Green’s functions using the target velocity structure itself, and then inverted the earthquake sources. This method greatly reduces the calculations required compared to a full inversion of the source mechanism, depth, and source duration (risetime), making it possible to reproduce input ground motions in the target basin. Here, we validated the Japan Integrated Velocity Structure model of the Osaka basin using five earthquakes around the Osaka basin; these earthquakes allowed us to investigate the impact of the incident waves’ propagation directions. We first estimated source parameters using records at control bedrock sites surrounding the Osaka basin to reproduce the input ground motions. Then, we conducted a 3D finite-difference simulation for sites within the basin. By mapping the distributions of misfit values for individual waveforms, peak ground velocities, and response spectra, we were able to identify areas in the basin that require additional tuning of the velocity model. Graphical Abstrac

Characterized source model of the M7.3 2016 Kumamoto earthquake by the 3D reciprocity GFs inversion with special reference to the velocity pulse at KMMH16

Abstract The 2016 Kumamoto earthquakes caused severe damage centering on the Mashiki residential area. The velocity waveforms at station KMMH16 in Mashiki, during the M7.3 mainshock, show large pulses. We found that severe damage in Mashiki may be the result of the strong westward velocity pulse. The question raised is how the near-fault ground motions with strong velocity pulse at KMMH16 were generated during the mainshock. We focus on the characterized source model with Strong Motion Generation Areas (SMGA). Empirical Green’s function (EGF) method is widely used for source modeling in this case. However, in case that the target site is located just near the fault in nodal plane of source mechanism (like KMMH16), mechanism of the EGF event should perfectly fit mechanism of the mainshock, which is a rare case. Therefore, instead of using EGFs, we used theoretical 3D Green’s functions. Our approach is a nonlinear source inversion. This method requires calculation of waveforms and comparison with observations for many source models. To accelerate these calculations, we use pre-calculated GFs by the reciprocity method in the JIVSM velocity structure model. By comparison with aftershock records, we validated this structure for periods as short as 1.5 s. Target sites are limited to sites close to the fault: KMM005, KMM006, KMMH14, and KMMH16. First, we look for an initial SMGA source model by the grid search method applied to relatively long-period (> 3 s) waveforms and coarse grid of source parameters. Then, we tune that source model to fit observed short-period waveforms with the simplex search method. Necessary physical constraints for the range of the source parameters are applied here. The important point in our inversion scheme is to describe the Kostrov-like slip velocity functions inside each SMGAs by using two triangles. The resulting source model agrees well with other inversion results. We found that the observed westward pulse at KMMH16 is the result of the constructive interference of two pulses from SMGA1 and SMGA2, located in Hinagu fault and southwestern segment of Futagawa fault. Graphical Abstrac

Earthquake Cycle Modelling of Multi-segmented Faults: Dynamic Rupture and Ground Motion Simulation of the 1992 Mw 7.3 Landers Earthquake

International audienceWe perform earthquake cycle simulations with the goal of studying the characteristics of source scaling relations and strong ground motions in multi-segmented fault ruptures. The 1992 Mw 7.3 Landers earthquake is chosen as a target earthquake to validate our methodology. The model includes the fault geometry for the three-segmented Landers rupture from the SCEC community fault model, extended at both ends to a total length of 200 km, and limited to a depth to 15 km. We assume the faults are governed by rate-and-state (RS) friction, with a heterogeneous, correlated spatial distribution of characteristic weakening distance Dc. Multiple earthquake cycles on this non-planar fault system are modeled with a quasi-dynamic solver based on the boundary element method, substantially accelerated by implementing a hierarchical-matrix method. The resulting seismic ruptures are recomputed using a fully-dynamic solver based on the spectral element method, with the same RS friction law. The simulated earthquakes nucleate on different sections of the fault, and include events similar to the Mw 7.3 Landers earthquake. We obtain slip velocity functions, rupture times and magnitudes that can be compared to seismological observations. The simulated ground motions are validated by comparison of simulated and recorded response spectra

Multicycle Simulation of Strike-Slip Earthquake Rupture for Use in Near-Source Ground-Motion Simulations

International audienceABSTRACT Realistic dynamic rupture modeling validated by observed earthquakes is necessary for estimating parameters that are poorly resolved by seismic source inversion, such as stress drop, rupture velocity, and slip rate function. Source inversions using forward dynamic modeling are increasingly used to obtain earthquake rupture models. In this study, to generate a large number of physically self-consistent rupture models, rupture process of which is consistent with the spatiotemporal heterogeneity of stress produced by previous earthquakes on the same fault, we use multicycle simulations under the rate and state (RS) friction law. We adopt a one-way coupling from multicycle simulations to dynamic rupture simulations; the quasidynamic solver QDYN is used to nucleate the seismic events and the spectral element dynamic solver SPECFEM3D to resolve their rupture process. To simulate realistic seismicity, with a wide range of magnitudes and irregular recurrence, several realizations of 2D-correlated heterogeneous random distributions of characteristic weakening distance (Dc) in RS friction are tested. Other important parameters are the normal stress, which controls the stress drop and rupture velocity during an earthquake, and the maximum value of Dc, which controls rupture velocity but not stress drop. We perform a parametric study on a vertical planar fault and generate a set of a hundred spontaneous rupture models in a wide magnitude range (Mw 5.5–7.4). We validate the rupture models by comparison of source scaling, ground motion (GM), and surface slip properties to observations. We compare the source-scaling relations between rupture area, average slip, and seismic moment of the modeled events with empirical ones derived from source inversions. Near-fault GMs are computed from the source models. Their peak ground velocities and peak ground accelerations agree well with the ground-motion prediction equation values. We also obtain good agreement of the surface fault displacements with observed values

Preliminary determination of the interdependence among strong motion amplitude, earthquake magnitude and hypocentral distance for the Himalayan region

Consiglio Nazionale delle Ricerche - Biblioteca Centrale - P.le Aldo Moro, 7 Rome / CNR - Consiglio Nazionale delle RichercheSIGLEITItal