Finite-fault Model Research Articles

Impact of earthquake source complexity and land elevation data resolution on tsunami hazard assessment and fatality estimation

This study investigates the impact of model complexity in source characterization and digital elevation model (DEM) resolution on the accuracy of tsunami hazard assessment and fatality estimation through a case study in Padang, Indonesia. Two types of earthquake source models, i.e. complex and uniform slip models, are adopted by considering three resolutions of DEMs, i.e. 150 m, 50 m, and 10 m. For each of the three grid resolutions, 300 complex source models are generated using new statistical prediction models of earthquake source parameters developed from extensive finite-fault models of past subduction earthquakes, whilst 100 uniform slip models are constructed with variable fault geometry without slip heterogeneity. The results highlight that significant changes to tsunami hazard and fatality estimates are observed with regard to earthquake source complexity and grid resolution. Coarse resolution (i.e. 150 m) leads to inaccurate tsunami hazard prediction and fatality estimation, whilst 50-m and 10-m resolutions produce similar results. However, velocity and momentum flux are sensitive to the grid resolution and hence, at least 10-m grid resolution needs to be implemented when considering flow-based parameters for tsunami hazard and risk assessments. In addition, the results indicate that the tsunami hazard parameters and fatality number are more sensitive to the complexity of earthquake source characterization than the grid resolution. Thus, the uniform models are not recommended for probabilistic tsunami hazard and risk assessments. Finally, the findings confirm that uncertainties of tsunami hazard level and fatality in terms of depth, velocity and momentum flux can be captured and visualized through the complex source modeling approach. From tsunami risk management perspectives, this indeed creates big data, which are useful for making effective and robust decisions.

Overview of Ground-Motion Issues for Cascadia Megathrust Events: Simulation of Ground-Motions and Earthquake Site Response

Ground motions for earthquakes of M7.5 to 9.0 on the Cascadia subduction interface are simulated based on a stochastic finite-fault model, and used to estimate average response spectra for reference firm soil conditions. The simulations are first validated by modeling the wealth of ground-motion data from the the 2011 M9.0 Tohoku earthquake of Japan. Adjustments to the calibrated model are then made to consider average source, attenuation and site parameters for the Cascadia region. This includes an evaluation of the likely variability in stress drop for large interface earthquakes, and an assessment of regional attenuation and site effects. We perform best-estimate simulations for a preferred set of input parameters. Typical results suggest mean values of 5%-damped pseudo-acceleration in the range from about 100 to 200 cm/s2, at frequencies from 1 to 4 Hz, for firm-ground conditions in Vancouver. Uncertainty in most-likely value of the parameter representing stress drop causes variability in simulated response spectra of about ±50%. Uncertainties in the attenuation model produce even larger variability in response spectral amplitudes – a factor of about two at a closest distance to the rupture plane (Rcd) of 100 km, becoming even larger at greater distances. It is thus important to establish the regional attenuation model for ground-motion simulations, and to bound the source properties controlling radiation of ground motion. We calculate theoretical 1D spectral amplification estimates for four selected Fraser River Delta sites to show how the presence of softer sediments in the region may alter the predicted ground motions. The amplification functions are largely consistent with observed spectral amplification at Fraser River delta sites, suggesting amplification by factors of 2.5 to 5 at the peak frequency of the site; we note that deep sites in the delta have a low peak frequency, ~0.3 Hz. This work will aid in seismic hazard assessment and mitigation efforts in the active Cascadia region of southwestern B.C. An important consideration is that the uncertainties are large and present a dominant unknown when assessing seismic risk. We find that variability in the expected motions exceeds a factor of two even on rock-like sites, with uncertainty

Uncertainties in the 2004 Sumatra-Andaman source through nonlinear stochastic inversion of tsunami waves.

Numerical inversions for earthquake source parameters from tsunami wave data usually incorporate subjective elements to stabilize the search. In addition, noisy and possibly insufficient data result in instability and non-uniqueness in most deterministic inversions, which are barely acknowledged. Here, we employ the satellite altimetry data for the 2004 Sumatra–Andaman tsunami event to invert the source parameters. We also include kinematic parameters that improve the description of tsunami generation and propagation, especially near the source. Using a finite fault model that represents the extent of rupture and the geometry of the trench, we perform a new type of nonlinear joint inversion of the slips, rupture velocities and rise times with minimal a priori constraints. Despite persistently good waveform fits, large uncertainties in the joint parameter distribution constitute a remarkable feature of the inversion. These uncertainties suggest that objective inversion strategies should incorporate more sophisticated physical models of seabed deformation in order to significantly improve the performance of early warning systems.

Imaging the 2016 Mw 7.8 Kaikoura, New Zealand, earthquake with teleseismic P waves: A cascading rupture across multiple faults

AbstractThe 13 November 2016 Mw 7.8 Kaikoura, New Zealand, earthquake was investigated using teleseismic P waves. Backprojection of high‐frequency P waves from two regional arrays shows unilateral rupture of at least two southwest‐northeast striking faults with an average rupture speed of 1.4–1.6 km/s and total duration of ~100 s. Guided by these backprojection results, 33 globally distributed low‐frequency P waves were inverted for a finite fault model (FFM) of slip. The FFM showed evidence of several subevents; however, it lacked significant moment release near the epicenter, where a large burst of high‐frequency energy was observed. A local strong‐motion network recorded strong shaking near the epicenter; hence, for this earthquake the distribution of backprojection energy is superior to the FFM as a guide of strong shaking. For future large earthquakes that occur in regions without strong‐motion networks, initial shaking estimates could benefit from backprojection constraints.

Stochastic ground motion simulation of the 2016 Meinong, Taiwan earthquake

We applied a stochastic method for the finite-fault modeling of strong ground motions to the 2016 Meinong, Taiwan earthquake. Newly developed attenuation models in Southern Taiwan with the frequency-dependent Q = 86.4f0.73 and the high-frequency decay factor κ0 were used in the synthetic model. The horizontal-to-vertical spectral ratios (HVSR) were calculated from weak motions and the Meinong mainshock and used for the site amplification correction of the synthetic waveforms produced by the stochastic ground motion simulation. Simulations incorporating the attenuation models and site correction improved the prediction of the S-wave envelope, duration, and peak ground acceleration (PGA). The nonlinear site response during the Meinong mainshock was identified by the degree of nonlinear site response (DNL), which is a summation of HVSR differences between weak motions and the Meinong mainshock as recorded by the Taiwan Strong Motion Instrument Program. The DNL showed a positive correlation with ground motion intensity. The surface site conditions influenced DNL strength. The percentage of PGA reduction calculated in this study can be an indicator of the spatial distribution of the degree of nonlinear soil effects on the Meinong earthquake in the time domain. Areas that had high levels of PGA reduction overlap with areas that had high liquefaction potential. Based on the residual analysis, forward directivity was identified in a 105° range in the northwestward direction. The amplification of forward rupture directivity was three times greater than the backward rupture directivity.Graphical .

Coseismic changes of gravitational potential energy induced by global earthquakes based on spherical‐Earth elastic dislocation theory

AbstractWe compute the coseismic gravitational potential energy Eg change using the spherical‐Earth elastic dislocation theory and either the fault model treated as a point source or the finite fault model. The rate of the accumulative Eg loss produced by historical earthquakes from 1976 to 2016 (about 42,000 events) using the Global Centroid Moment Tensor Solution catalogue is estimated to be on the order of −2.1 × 1020 J/a, or −6.7 TW (1 TW = 1012 W), amounting to 15% in the total terrestrial heat flow. The energy loss is dominated by the thrust faulting, especially the megathrust earthquakes such as the 2004 Sumatra earthquake (Mw 9.0) and the 2011 Tohoku‐Oki earthquake (Mw 9.1). It is notable that the very deep focus events, the 1994 Bolivia earthquake (Mw 8.2) and the 2013 Okhotsk earthquake (Mw 8.3), produced significant overall coseismic Eg gain according to our calculation. The accumulative coseismic Eg is mainly lost in the mantle of the Earth and also lost in the core of the Earth but with a relatively smaller magnitude. By contrast, the crust of the Earth gains gravitational potential energy cumulatively because of the coseismic deformations. We further investigate the tectonic signature in the coseismic crustal Eg changes in some complex tectonic zone, such as Taiwan region and the northeastern margin of the Tibetan Plateau. We found that the coseismic Eg change is consistent with the regional tectonic character.

The finite, kinematic rupture properties of great-sized earthquakes since 1990

Here, I present a database of >160 finite fault models for all earthquakes of M 7.5 and above since 1990, created using a consistent modeling approach. The use of a common approach facilitates easier comparisons between models, and reduces uncertainties that arise when comparing models generated by different authors, data sets and modeling techniques.I use this database to verify published scaling relationships, and for the first time show a clear and intriguing relationship between maximum potency (the product of slip and area) and average potency for a given earthquake. This relationship implies that earthquakes do not reach the potential size given by the tectonic load of a fault (sometimes called “moment deficit,” calculated via a plate rate over time since the last earthquake, multiplied by geodetic fault coupling). Instead, average potency (or slip) scales with but is less than maximum potency (dictated by tectonic loading). Importantly, this relationship facilitates a more accurate assessment of maximum earthquake size for a given fault segment, and thus has implications for long-term hazard assessments. The relationship also suggests earthquake cycles may not completely reset after a large earthquake, and thus repeat rates of such events may appear shorter than is expected from tectonic loading. This in turn may help explain the phenomenon of “earthquake super-cycles” observed in some global subduction zones.

Factors controlling high-frequency radiation from extended ruptures

Small-scale slip heterogeneity or variations in rupture velocity on the fault plane are often invoked to explain the high-frequency radiation from earthquakes. This view has no theoretical basis, which follows, for example, from the representation integral of elasticity, an exact solution for the radiated wave field. The Fourier transform, applied to the integral, shows that the seismic spectrum is fully controlled by that of the source time function, while the distribution of final slip and rupture acceleration/deceleration only contribute to directivity. This inference is corroborated by the precise numerical computation of the full radiated field from the representation integral. We compare calculated radiation from four finite-fault models: (1) uniform slip function with low slip velocity, (2) slip function spatially modulated by a sinusoidal function, (3) slip function spatially modulated by a sinusoidal function with random roughness added, and (4) uniform slip function with high slip velocity. The addition of “asperities,” both regular and irregular, does not cause any systematic increase in the spectral level of high-frequency radiation, except for the creation of maxima due to constructive interference. On the other hand, an increase in the maximum rate of slip on the fault leads to highly amplified high frequencies, in accordance with the prediction on the basis of a simple point-source treatment of the fault. Hence, computations show that the temporal rate of slip, not the spatial heterogeneity on faults, is the predominant factor forming the high-frequency radiation and thus controlling the velocity and acceleration of the resulting ground motions.

Comparisons of Co-Seismic Gravity Changes between GRACE Observations and the Predictions from the Finite-Fault Models for the 2012 Mw = 8.6 Indian Ocean Earthquake Off-Sumatra

Comparisons of Co-Seismic Gravity Changes between GRACE Observations and the Predictions from the Finite-Fault Models for the 2012 Mw = 8.6 Indian Ocean Earthquake Off-Sumatra

REGARD: A new GNSS‐based real‐time finite fault modeling system for GEONET

AbstractThe short‐period seismometer‐based magnitude saturation problem, especially for events with magnitude > 8, can be improved by a real‐time Global Navigation Satellite System (GNSS) positioning technique, which has enabled rapid estimation of a finite fault model for a large earthquake without any saturation. A new real‐time fault modeling system based on the GNSS Earth Observation Network (GEONET) is developed and is under experimental operation in Japan. In this paper, we present the newly developed system REGARD (the Real‐time GEONET Analysis system for Rapid Deformation monitoring), which consists of real‐time GNSS positioning, automatic detection of coseismic displacement by the event, and quasi real‐time finite fault model inversion routines. The performance of the automatic event detection system is tested through experimental real‐time operation based on GEONET data for 2 months. Furthermore, we also test the reliability of the finite fault model inversion routines using real raw GNSS data observed for past large earthquakes: the 2003 Tokachi‐oki earthquake (moment magnitude (Mw) 8.3), the 2011 Tohoku earthquake (Mw 9.0), and the 2011 off‐Ibaraki earthquake (Mw 7.7). A simulated 1707 Hoei‐type Nankai Trough earthquake (Mw 8.7) is also tested. The real‐time experimental operation shows that real‐time GNSS positioning is precise enough to detect all the tested earthquakes, and the inversion results demonstrate that the REGARD can reliably estimate the earthquake size and its extent within 3 min after the origin time. These results suggest that the REGARD system will complement the seismometer‐based magnitude determination system.

Hybrid broadband simulation of strong-motion records from the September 16, 1978, Tabas, Iran, earthquake (M w 7.4)

This paper presents a simulation of three components of near-field ground shaking recorded during the main shock at three stations of the September 16, 1978, Tabas (M w = 7.4), Iran, earthquake, close to the causative fault. A hybrid method composed of a discrete wavenumber method developed by Bouchon (Bouchon in Bull Seismol Soc Am 71:959–971, 1981; Cotton and Coutant in Geophys J Int 128:676–688, 1997) and a stochastic finite-fault modeling based on a dynamic corner frequency proposed by Motazedian and Atkinson (Bull Seismol Soc Am 95:995–1010, 2005), modified by Assatourians and Atkinson (Bull Seismol Soc Am 97:935–1949, 2007), is used for generating the seismograms at low (0.1–1.0 Hz) and high frequencies (1.0–20.0 Hz), respectively. The results are validated by comparing the simulated peak acceleration, peak velocity, peak displacement, Arias intensity, the integral of velocity squared, Fourier spectrum and acceleration response spectrum on a frequency-by-frequency basis, the shape of the normalized integrals of acceleration and velocity squared, and the cross-correlation with the observed time-series data. Each characteristic is compared on a scale from 0 to 10, with 10 being perfect agreement. Also, the results are validated by comparing the simulated ground motions with the modified Mercalli intensity observations reported by reconnaissance teams and showed reasonable agreement. The results of the present study imply that the damage distribution pattern of the 1978 Tabas earthquake can be explained by the source directivity effect.

Strong Aftershock Study Based on Coulomb Stress Triggering—A Case Study on the 2016 Ecuador Mw 7.8 Earthquake

The 2016 Ecuador M 7.8 earthquake ruptured the subduction zone boundary between the Nazca plate and the South America plate. This M 7.8 earthquake may have promoted failure in the surrounding crust, where six M ≥ 6 aftershocks occurred following this mainshock. These crustal ruptures were triggered by the high coulomb stress changes produced by the M 7.8 mainshock. Here, we investigate whether the six M ≥ 6 aftershocks are consistent with the positive coulomb stress region due to the mainshock. To explore the correlation between the mainshock and the aftershocks, we adopt a recently published high-quality finite fault model and focal mechanisms to study the coulomb stress triggers during the M 7.8 earthquake sequence. We compute the coulomb failure stress changes (ΔCFS) on both of the focal mechanism nodal planes. We compare the ΔCFS imparted by the M 7.8 mainshock on the subsequent aftershocks with the epicenter location of each aftershock. In addition, the shear stress, normal stress, and coulomb stress changes in the focal sources of each aftershock are also computed. Coulomb stress changes in the focal source for the six M ≥ 6 aftershocks are in the range of −2.17–7.564 bar. Only one computational result for the M 6.9 aftershock is negative; other results are positive. We found that the vast majority of the six M ≥ 6 aftershocks occurred in positive coulomb stress areas triggered by the M 7.8 mainshock. Our results suggest that the coulomb stress changes contributed to the development of the Ecuador M 7.8 earthquake sequence.

Empirical fatality model for Indonesia based on a Bayesian approach

Abstract An empirical fatality model for Indonesia has been developed by relating the macroseismic intensity to the fatality rate using compiled sub-district level fatality rate data and the numerically simulated ground-shaking intensity for four recent damaging events. The fatality rate data were compiled by collecting population and fatality statistics of the regions affected by the selected events. The ground-shaking intensity was numerically estimated by incorporating a finite-fault model of each event and local site conditions approximated by topographically based site amplifications. The macroseismic intensity distribution of each event was generated using ShakeMap software, combining a selected pair of ground motion prediction equations and ground motion to intensity conversion equations. The developed fatality model is a Bayesian generalized linear model in which the fatality rate is assumed to follow a mixture of Bernoulli and gamma distributions. The model was validated by calculating the fatalities in past events from the EXPO-CAT catalogue and comparing the estimates with the EXPO-CAT fatality records. Although the model can provide an estimate of the range of fatalities for future events, it needs ongoing refinement by the incorporation of additional fatality rate data from past and future events.

Tsunami Hazard Analysis of Future Megathrust Sumatra Earthquakes in Padang, Indonesia Using Stochastic Tsunami Simulation

This study assesses the tsunami hazard potential in Padang, Indonesia probabilistically using a novel stochastic tsunami simulation method. The stochastic tsunami simulation is conducted by generating multiple earthquake source models for a given earthquake scenario, which are used as input to run Monte Carlo tsunami simulation. Multiple earthquake source models for three magnitude scenarios, i.e. Mw 8.5, Mw 8.75, and Mw 9.0, are generated using new scaling relationships of earthquake source parameters developed from an extensive set of 226 finite-fault models. In the stochastic tsunami simulation, the effect of incorporating and neglecting the prediction errors of earthquake source parameters is investigated. In total, 600 source models are generated to assess the uncertainty of tsunami wave characteristics and maximum tsunami wave height profiles along coastal line of Padang. The results highlight the influence of the uncertainty of the scaling relationships on tsunami simulation results and provide a greater range of tsunamigenic scenarios produced from the stochastic tsunami simulation. Additionally, the results show that for the future major earthquakes in the Sunda megathrust, the maximum tsunami wave height in Padang areas can reach 20 m and therefore, significant damage and loss may be anticipated in this region.

Kinematic finite fault and 3D seismic wave propagation of the 24 August, 2016, Mw 6.0 central Italy earthquake

The magnitude Mw 6.0 earthquake of 24th August 2016 caused severe damages and nearly 300 fatalities in the central Italy region. Initial reports revealed an asymmetrical distribution of damage and coseismic effects, suggesting a major role of heterogeneities, both in the rupture history and in the geological structure of the region. Near realtime availability of seismological data afforded a timely determination of a finite fault model (Tinti et al., 2016). Here we test this source model by performing a 3D simulation of seismic wave propagation within a 3D structural model containing the major geological features of the region. Agreement between modeled seismograms and observed seismograms suggests that some complexities in the waveforms, such as high amplification in the region of the Mt. Vettore fault system, can be accounted for by complexities in the fault rupture and 3D structural models. Finally, the consistency of the hypothesis of two distinct events has been analyzed.

The first month of the 2016 central Italy seismic sequence: fast determination of time domain moment tensors and finite fault model analysis of the ML 5.4 aftershock

<p>We present the revised Time Domain Moment Tensor (TDMT) catalogue for earthquakes with M_L larger than 3.6 of the first month of the ongoing Amatrice seismic sequence (August 24th - September 25th). Most of the retrieved focal mechanisms show NNW–SSE striking normal faults in agreement with the main NE-SW extensional deformation of Central Apennines. We also report a preliminary finite fault model analysis performed on the larger aftershock of this period of the sequence (M_w 5.4) and discuss the obtained results in the framework of aftershocks distribution.</p>

Simulations of long‐period ground motions from a large earthquake using finite rupture modeling and the ambient seismic field

AbstractLong‐period ground motions generated by large earthquakes slowly attenuate with distance and can be significantly amplified by local velocity structures even at large distances. We take advantage of the wave propagation information carried by the ambient seismic field to simulate the long‐period ground motions of the 2008 Mw 6.9 Iwate‐Miyagi Nairiku earthquake, which occurred in the Tohoku region, Japan. We extract Green's functions between pairs of stations by regarding seismometers located in the vicinity of the main shock fault plane as virtual sources and others as receivers. We calibrate the amplitude of the extracted Green's functions using records of a Mw 5.0 aftershock with a reverse‐faulting focal mechanism similar to that of the main shock. We use scaling relations between small and large earthquakes to construct several finite fault models that are used together with the extracted Green's functions to simulate the ground motions of the main shock. Among the different models, a uniform finite source model with a rupture velocity of 1.95 km/s combined with the Green's functions extracted using one virtual source station located at the center of the fault plane yields the best fit between the simulated and observed long‐period ground motions. This study demonstrates the potential of the ambient seismic field combined with finite source modeling to assess the seismic hazard related to long‐period ground motions that could be generated by forthcoming large earthquakes.

First result from the GEONET real-time analysis system (REGARD): the case of the 2016 Kumamoto earthquakes

We present the initial results of rapid fault estimations for the 2016 Kumamoto earthquake on April 16 (M j 7.3), and coseismic displacements caused by the two large foreshocks that occurred on April 14 (M j 6.5) and April 15 (M j 6.4) from the GEONET real-time analysis system (REGARD), which is based on a Global Navigation Satellite System (GNSS) kinematic positioning technique. The real-time finite-fault estimate (M w 6.85) was obtained within 1 min and converged to M w 6.96 within 5 min of the origin time of the mainshock (M j 7.3). The finite-fault estimate shows right-lateral strike-slip fault along the Futagawa fault segment, which is consistent with the finite-fault model inferred from post-processed GNSS and InSAR analysis. Furthermore, significant coseismic displacements were observed due to the April 14 and April 15 foreshocks at nearby sites, though these earthquakes were smaller than the pre-assigned system threshold. Our results also demonstrate the potential for the GNSS-based earthquake early warning system for inland earthquakes.

Stochastic ground-motion simulations for the 2016 Kumamoto, Japan, earthquake

On April 15, 2016, Kumamoto, Japan, was struck by a large earthquake sequence, leading to severe casualty and building damage. The stochastic finite-fault method based on a dynamic corner frequency has been applied to perform ground-motion simulations for the 2016 Kumamoto earthquake. There are 53 high-quality KiK-net stations available in the Kyushu region, and we employed records from all stations to determine region-specific source, path and site parameters. The calculated S-wave attenuation for the Kyushu region beneath the volcanic and non-volcanic areas can be expressed in the form of Q s = (85.5 ± 1.5)f 0.68±0.01 and Q s = (120 ± 5)f 0.64±0.05, respectively. The effects of lateral S-wave velocity and attenuation heterogeneities on the ground-motion simulations were investigated. Site amplifications were estimated using the corrected cross-spectral ratios technique. Zero-distance kappa filter was obtained to be the value of 0.0514 ± 0.0055 s, using the spectral decay method. The stress drop of the mainshock based on the USGS slip model was estimated optimally to have a value of 64 bars. Our finite-fault model with optimized parameters was validated through the good agreement of observations and simulations at all stations. The attenuation characteristics of the simulated peak ground accelerations were also successfully captured by the ground-motion prediction equations. Finally, the ground motions at two destructively damaged regions, Kumamoto Castle and Minami Aso village, were simulated. We conclude that the stochastic finite-fault method with well-determined parameters can reproduce the ground-motion characteristics of the 2016 Kumamoto earthquake in both the time and frequency domains. This work is necessary for seismic hazard assessment and mitigation.

On the scale dependence of earthquake stress drop

We discuss the debated issue of scale dependence in earthquake source mechanics with the goal of providing supporting evidence to foster the adoption of a coherent interpretative framework. We examine the heterogeneous distribution of source and constitutive parameters during individual ruptures and their scaling with earthquake size. We discuss evidence that slip, slip-weakening distance and breakdown work scale with seismic moment and are interpreted as scale dependent parameters. We integrate our estimates of earthquake stress drop, computed through a pseudo-dynamic approach, with many others available in the literature for both point sources and finite fault models. We obtain a picture of the earthquake stress drop scaling with seismic moment over an exceptional broad range of earthquake sizes (−8 < MW < 9). Our results confirm that stress drop values are scattered over three order of magnitude and emphasize the lack of corroborating evidence that stress drop scales with seismic moment. We discuss these results in terms of scale invariance of stress drop with source dimension to analyse the interpretation of this outcome in terms of self-similarity. Geophysicists are presently unable to provide physical explanations of dynamic self-similarity relying on deterministic descriptions of micro-scale processes. We conclude that the interpretation of the self-similar behaviour of stress drop scaling is strongly model dependent. We emphasize that it relies on a geometric description of source heterogeneity through the statistical properties of initial stress or fault-surface topography, in which only the latter is constrained by observations.