Self-similar stochastic slip distributions on a non-planar fault for tsunami scenarios for megathrust earthquakes

doi:10.21203/rs.3.rs-32273/v1

Download PDF

Research article

Self-similar stochastic slip distributions on a non-planar fault for tsunami scenarios for megathrust earthquakes

https://doi.org/10.21203/rs.3.rs-32273/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 31 Aug, 2020

Read the published version in Progress in Earth and Planetary Science →

You are reading this older preprint version

Read the latest preprint version →

Megathrust earthquakes that occur repeatedly along the plate interface of subduction zones can cause severe damage due to strong ground motion and the destructive tsunamis they can generate. We developed a set of scenario earthquakes to evaluate tsunami hazards and tsunami early warning systems for such devastating earthquakes. Although it is known that the slip distribution on a fault strongly affects the tsunami height distribution in near-field coastal areas, the slip distribution of future earthquakes cannot be exactly predicted. One way to resolve this difficulty is to create a set of scenario earthquakes in which a set of heterogeneous slip distributions on the source fault is stochastically generated based on a given slip probability density function (SPDF). The slip distributions generated in this manner differ from event to event, but their average over a large ensemble of models converges to a predefined SPDF resembling the long-term average of ruptures on the target fault zone. We created a set of SPDF-based scenario earthquakes for an expected future M_w 8.2 Tonankai earthquake in the eastern half of the Nankai trough, off southwest Japan, and computed the ensuing tsunamis. We found that the estimated peak coastal amplitudes among the ensemble of tsunamis along the near-field coast differed by factors of 3 to 9 and the earliest and latest arrivals at each observation site differed by 400 to 700 s. The variations in both peak tsunami amplitude and arrival time at each site were well approximated by a Gaussian distribution. For cases in which the slip distribution is unknown, the average and standard deviation of these scenario datasets can provide first approximations of forecast tsunami height and arrival time and their uncertainties, respectively. At most coastal observation sites, tsunamis modelled similarly but using a uniform slip distribution underpredicted tsunami amplitudes but gave earlier arrival times than those modeled with a heterogeneous slip distribution. Use of these earlier arrival times may be useful for providing conservative early warnings of tsunami arrivals. Therefore, tsunami computations for both heterogeneous and uniform slip distributions are important for tsunami disaster mitigation.

Planetary Science

Nankai trough

Tsunami early warning

Scenario earthquakes

Megathrust earthquakes along the Nankai trough off southwestern Japan, such as the 1944 Tonankai and 1946 Nankai events, have caused severe damage due to strong ground motion and large tsunamis. Devastating earthquakes such as these have occurred at intervals of 100 to 150 years (e.g., Ando 1975; Ishibashi 2004), so we should be prepared for more of them in the future. Probabilistic analyses of earthquake and tsunami hazards in this region have been based on sets of scenario earthquakes of various magnitudes and rupture areas (e.g., Earthquake Research Committee, Headquarters for Earthquake Research Promotion in Japan 2013, 2020). These analyses have also been used to develop tsunami early warning systems to mitigate the disastrous effects of the tsunamis generated (e.g. Yamamoto et al. 2016b).

Tsunami warning systems generally estimate coastal tsunami heights and arrival times on the basis of earthquake source parameters, crustal deformation, and initial tsunami heights estimated by inversion of seismic, geodetic, or offshore ocean-bottom pressure gauge (OBP) records, respectively (Tatehata 1997; Titof et al. 2005; Blewitt et al. 2009; Tsushima et al. 2009, 2011, 2012, 2014). Tsunami early warning systems of this type are operated by the Japan Meteorological Agency (JMA) (Tatehata 1997, Tsushima et al. 2012) and NOAA (Titof et al. 2005).

In recent decades, to allow earlier and more precise issuance of earthquake and tsunami warnings, dense offshore seismic and OBP networks have been deployed on the ocean floor above the source areas of expected large earthquakes; these include DONET along the Nankai trough (Kaneda et al. 2015; Kawaguchi et al. 2015) and S-net along the Japan Trench (Kanazawa et al. 2016; Mochizuki et al. 2016; Uehira et al. 2016). Yamamoto et al. (2016b) developed a method to rapidly estimate initial tsunami heights by using dense OBP networks. The deployment of OBP networks in the near field of seismic sources has accelerated the development of other approaches for tsunami early warnings that need to estimate neither earthquake nor tsunami sources. For example, Maeda et al. (2015) developed a method to use real-time OBP data to estimate coastal tsunami heights based on data assimilation.

Simple methods have also been proposed that estimate coastal tsunami heights by applying amplification factors to offshore tsunami heights recorded at OBP sites (Baba et al. 2014a; Igarashi et al. 2016), or by selecting candidate tsunami sources based on spatial correlations of observed tsunami height distributions at OBP sites with those expected from scenario earthquakes (Yamamoto et al. 2016a; Takahashi et al. 2017, 2018). These approaches require a pre-computed database of tsunami parameters for a set of scenario earthquakes in the target area of various magnitudes, source locations, and source geometries. After creation of such a tsunami database, the computational cost of estimating coastal heights and arrival times is lower. Because of the low cost and simplicity of implementation of this type of system, without the need for near-real time source estimations, they are used by some local governments and suppliers of lifeline services (e.g., Takahashi et al. 2017, 2018; Ishibashi et al. 2018).

The accuracy of tsunami heights estimated using scenario earthquakes depends on the ensemble of scenarios providing adequate coverage of the range of possible future earthquakes. To date, most studies have assumed uniform slip distributions on planar faults in their scenario earthquakes, whereas the spatial distribution of slip is in reality heterogeneous, and this heterogeneity affects coastal tsunami wave height distributions in the near field (Geist 2002; Mueller et al. 2015; Li et al. 2016; Melger et al. 2019). Although it is impossible to predict the slip distributions of future earthquakes, the use of large suites of rapidly produced stochastic slip distributions offers a means to account for some of the inherent uncertainties.

In general, stochastic fault slip distributions can be well approximated by the k-square model, where the wavenumber spectrum decays as k^− 2 beyond the corner wavenumber k_c, which is dependent on earthquake magnitude (Herrero and Bernard 1994; Mai and Beroza 2002; Causse et al. 2010). There are two principal methods used to produce stochastic self-similar slip distributions. One is based on describing the slip spectra in Fourier space, which is then converted to the spatial domain using fast Fourier transforms (Herrero and Bernard, 1994). The alternative, which is the method we used, is the composite source model. It involves the placement of a large number of circular subevents (i.e., patches of slip) whose size-frequency distribution on the fault plane follows a power law (Frankel 1991; Zeng et al. 1994; Ruiz et al. 2011), such that the summation of all the subevents produces a self-similar slip distribution. The placement of these subevents is based on a probability density function (PDF) defined spatially across the fault plane, hereafter referred to as the slip probability density function (SPDF). Individual earthquakes generated using a particular SPDF will differ from event to event, but over a large ensemble of earthquakes, the average slip will converge to the predefined SPDF (Murphy and Herrero 2020). Therefore, the choice of the function used to describe the SPDF depends on a priori knowledge or assumptions of the likely first-order spatial distribution of slip. For example, in situations where there is no knowledge of preferential locations for slip, a uniform PDF could be used (Murphy and Herrero 2020); whereas if a specific earthquake or a particular rupture feature is assumed to occur (e.g., shallow slip amplification) a more specific function can be applied (Murphy et al. 2016; Scala et al. 2020).

Our aim in this study was to create a family of megathrust earthquakes that adequately represents possible scenarios that might occur in the Nankai trough. The problem that remained was how to set the SPDF. Recent studies have indicated that the coseismic slip distributions of megathrust earthquakes are possibly correlated with plate interface coupling (Hashimoto et al. 2009, 2012; Perfettini et al. 2010; Moreno et al. 2010), but in some instances, stress accumulated during an interseismic period is partially released, resulting in partial rupture of a strongly coupled region (Konca et al. 2008). These observations suggest that the probability of slip during a megathrust earthquake may be proportional to the interplate coupling. Based on this hypothesis, for the SPDF we used the slip deficit rate (SDR) obtained along the Nankai trough by Yokota et al. (2016). Using this technique, we generated 200 scenario earthquakes for an anticipated M_w 8.2 Tonankai earthquake, which would rupture the eastern half of the Nankai trough. By computing the ensuing tsunamis, we evaluated the probability distributions of peak tsunami amplitudes and arrival times along the coast facing the source area to evaluate potential inherent errors in early warnings. We also compared these scenarios with tsunamis computed on the basis of a uniform slip distribution on the same fault area. Our aim was to determine the relative importance of tsunami computations based on heterogeneous and uniform slip distributions for effective tsunami early warnings.

Fault geometry of the Nankai trough subduction zone

The geometry we used for the subducting Philippine Sea plate is based on that of the Japan Integrated Velocity Structure Model (JIVSM; Koketsu et al. 2012) with the deep plate interface placed at the boundary between oceanic and continental crust (Figure 1). As a rupture on the plate interface approaches the surface the rupture plane rises from the top of the oceanic crust and cuts through the overlying sedimentary sequence to reach the trough axis at the seafloor (e.g., Moore et al. 2007, 2009). Following the method of Ide et al. (2010), we represented the plate interface by using an adjustable tension, continuous-curvature, surface-gridding algorithm (Smith and Wessel 1990) to define a surface that smoothly connects the top of the oceanic plate and the trough axis. We represented the non-planar fault of the Nankai subduction zone defined above by using a triangular mesh with nodes at 5 km intervals that extended from the seafloor to a depth of 50 km (Figure 1).

Stochastic scenarios for a future Tonankai earthquake

As a test case, we generated a set of scenario earthquakes for a future Tonankai earthquake that ruptures the northeastern half of the Nankai trough in an area similar to those of the 1854 and 1944 events (Ando 1975). We assumed an earthquake of M_w 8.2 in a source area in the northeastern part of the Nankai trough (Earthquake Research Committee, the Headquarters for Earthquake Research Promotion in Japan, 2013) (Figure 1). We assumed the SPDF to be proportional to the SDR obtained by Yokota et al. (2016) and set the probability to 0 beyond the southwestern boundary of the source area shown in Figure 1. We used common rupture area to generate tsunamigenic scenario earthquakes so that we could evaluate the uncertainties of tsunami heights and arrival times in that region without knowledge of the slip distribution on the earthquake fault. We also considered the different effects if the tsunamis were generated by an earthquake with uniform slip on the same rupture area.

To stochastically generate heterogeneous slip distributions on the non-planar fault, we used the method of Herrero and Murphy (2018), but adapted it to facilitate surface ruptures on the seafloor, that is, to allow large fault slip near the trough axis (Murphy and Herrero 2020). Because an M_w 8.2 earthquake would occur in an area smaller than the entire Nankai subduction zone, we needed to subdivide the meshed fault plane into areas of slip (rupture areas) and no slip. This was done by first randomly choosing a single cell from the SPDF generated from the SDR of Yokota et al. (2016). Cells from the mesh were then added to the rupture area one at a time until the rupture area marginally exceeded the empirical area–magnitude scaling relationship of Strasser et al. (2010). Cells for which at least one of its three sides would be common with a cell of the pre-existing rupture area were selected. When the rupture area had been established, families of circular subevents that had a power-law size distribution, each with an Eshelby slip profile (Eshelby 1957; Ruiz et al. 2011), were placed within the rupture area on the fault plane. The distribution of these subevents was based on the original SPDF after it had been re-normalized over the now-defined rupture area.

To generate stochastic slip distributions, we used two parameter sets: = (0.2, 60) and (0.25, 20), where is the maximum radius of subevents relative to fault width and is the total number of subevents used to generate one slip distribution. We generated 100 scenarios for each parameter set and used the resultant 200 scenarios in our analysis.

Figure 2 shows examples of individual slip distributions for the two parameter sets and the average slip distribution for all models. Correlations between the stacks of increasing numbers of slip distributions and the SDR for the two parameter sets (Figure 3) show that the correlations converge when there are more than about 50 slip distributions. This confirms the validity of our assumption that the average of many slip distributions converges to the SDR distribution and that the ensemble of slip distributions is therefore representative of the SDR.

We also generated a uniform slip distribution for the same rupture area, in which we averaged 100 slip distributions randomly generated from = (0.1, 50,000) with a uniform SPDF.

Tsunami computations

We computed the tsunamis generated by each scenario earthquake. Vertical displacements due to dip-slip motion of the triangular fault elements were computed using the dislocation method of Meade (2007) and the total seafloor displacement for each scenario was computed by summing the contributions from all the elements within the rupture area. Tsunamis due to these seafloor displacements were computed using the JAGURS code (Baba et al. 2014b). We solved the linear shallow-water equations in a finite difference scheme with a nesting algorithm comprising three layers of which the finest grid spacing was 2 arc-seconds. We used bathymetry from the M7000 Digital Bathymetric Chart provided by the Japan Hydrographic Association. We assumed a 40 s rise-time for crustal deformation for M_w 8.2 events. Tsunami observation sites were placed at 2.5 km intervals on the 3 m isobath along the coastline facing the source area (Figure 4). We included vertical crustal deformations in the tsunami heights because we focused on the near-field effect of the tsunamis. We simulated the tsunami propagations for 2 h with a time step of 0.1 s to satisfy the stability conditions solving for the finite-difference methods.

We measured the peak near-shore tsunami amplitude (PNTA) and arrival time of each of the tsunamis generated by the 200 scenario earthquakes. We defined the arrival time as the time of arrival of the first tsunami for which the wave-height exceeded 0.2 m above sea level. This threshold is based on the JMA ‘tsunami advisory’, whereby a warning is issued if a tsunami height is expected to exceed 0.2 m above sea level at a target site.

The conditional probability of exceedance of PNTA at each site (Figure 5a) and the ratio of the maximum and minimum values of PNTA at each site (Figure 5b) indicate that the variations in computed tsunami heights are related to the heterogeneous slip distributions. At most sites, PNTA values differed by factors of 3 to 4, depending on the slip distribution, but they differed by as much as a factor of 9 at some sites. After removing the highest and lowest PNTA values to consider the middle 95% of scenarios, these differences were reduced to a factor of about 3 at most sites (Figure 5b), implying that the largest variations occurred at only a few extreme cases.

The distribution of the cumulative probability of exceedance of PNTA at each site (Figure 6a) is well approximated by a Gaussian distribution, which indicates that the average and standard deviation can be used to make first approximation estimates of tsunami height and its error, respectively, for early warning systems when the slip distribution is unknown.

The differences between the earliest and latest tsunami arrival times at each observation site varied by 400 to 700 s depending on the scenario (Figure 7) but, as was the case for PNTA, were reduced (to 300 to 600 s) after removal of the most extreme cases. The distribution of the probability of exceedance at each site is again well approximated by a Gaussian distribution (Figure 6b).

To investigate the magnitude of the effect of slip distribution on PNTA, we compared the scenarios that caused the largest and smallest PNTAs near the Shima Peninsula (location in Figure 4). The slip distributions that generated both the largest and smallest PNTAs are broadly similar (Figures 8a and 8c), with slip concentrated mainly at a large asperity close to the Shima Peninsula, but the patterns of initial seafloor deformation for these two scenarios differ considerably. The areas of seafloor uplift and subsidence (Figures 8b and 8d) are distributed roughly symmetrically and parallel to the coast for the largest tsunami, whereas for the smallest tsunami the uplift area partially wraps around the area of subsidence.

For the largest tsunami, the propagation of the initial wave is relatively simple. The wave packet of greatest amplitude propagated directly towards the Shima Peninsula with some focusing due to bathymetry (Figure 4), which led to the largest tsunami reaching the coast of the Shima Peninsula (Figure 9a, Movie01). For the smallest tsunami, the largest amplitude wave packet propagated toward the Enshu Nada coast (location in Figure 4) east of the Shima Peninsula (Figure 9b, Movie02), because the largest asperity was farther east than that of the largest tsunami (Figure 8c). Furthermore, the later coalescence of tsunami waves with positive and negative initial amplitudes west and east, respectively, of the Shima Peninsula (t = 120 s in Figure 9b) reduced the tsunami amplitude at the peninsula. These results show that the near-field PNTA distribution can be strongly affected by both fault-slip distributions and interference of tsunami waves coming from different directions. Bathymetry can also affect PNTA distributions. This modelling demonstrates how difficult it can be to estimate PNTA distributions for future tsunamigenic earthquakes and shows that for effective tsunami early warning and disaster prevention it is important to consider a wide variety of scenarios.

Comparison of the PNTAs of the scenario with uniform fault slip distribution and those with heterogeneous slip distributions shows that at almost all observation sites the PNTA for the uniform slip distribution is smaller than those of about 50% of the heterogeneous slip scenarios (Figure 5a). At several sites, all the PNTAs for the heterogeneous slip scenarios were greater than that of the uniform slip scenario (e.g., in Figure 5a the uniform slip PNTA curve is at or below the 100% probability of exceedance for the heterogeneous case at stations 110–170). This result is consistent with previous studies that have shown that compared with scenarios with heterogeneous slip distributions, a uniform slip distribution mostly underpredicts the near-field PNTA (Geist 2002; Mueller et al. 2015; Li et al. 2016; Melger et al. 2019). This underprediction is the result of asperities concentrating large slip within limited areas in the source region that in turn excite larger amplitude tsunami waves than would be the case with a uniform slip distribution within the same source area. In contrast, the uniform slip distribution scenario produced relatively small but widely distributed slip in the source area, thus causing smaller initial seafloor displacements over a wider area (Figure 8f). Consequently, the tsunami wavefield generated was also wider with smaller amplitude, resulting in smaller PNTAs (Figure 10, Movie03).

At most observation sites, the tsunami arrival time for the uniform slip scenario was earlier than those for 50% of the heterogeneous slip scenarios, and earlier than all those scenarios at about 20% of sites (Figure 7a). This is because the wider distribution of seafloor uplift for the scenario with uniform slip (Figure 8f) resulted in shorter propagation distances and therefore earlier tsunami arrivals at most sites. Accordingly, the earliest arrival times for both the heterogeneous and uniform slip scenarios should be considered in examinations of tsunami hazards.

Clearly, near-field PNTA and arrival time distributions depend on the scenarios that were used to generate them. We based all of our scenarios on M_w 8.2 events on the same source area in the Tonankai earthquake source region. Coastal PNTA distributions will change if scenarios with different magnitudes and source areas are added. Moreover, we used the SDR distribution of Yokota et al. (2016) to obtain the SPDF we used. Use of different SDRs (e.g., Noda et al. 2018; Nishimura et al. 2018) for the Nankai trough will also affect the source area and slip distributions for the scenario earthquakes. Both Noda et al. (2018) and Nishimura et al. (2018) obtained smaller SDRs beneath the Enshu Nada than did Yokota et al. (2016), the use of which may in turn affect PNTAs along the coast. In the future, these (and perhaps other) SDR distributions, with appropriate weightings, should be incorporated in PNTA and travel time calculations. The geometry of the plate interface, especially the shallow part, will also affect the patterns of crustal deformation and the tsunamis they generate. The slip on splay faults can also affect the generation of tsunamis (Park et al. 2002). Further careful investigations of these factors are necessary to properly evaluate tsunami hazards and prepare appropriate ensembles of scenario earthquakes for use in providing effective tsunami early warnings.

We developed a procedure to stochastically generate slip distributions on a non-planar fault surface to generate scenario earthquakes for potential use in tsunami early warning applications. We generated 200 scenario earthquakes for an anticipated M_w 8.2 Tonankai earthquake in the Nankai trough and found that the estimated near-field tsunami amplitudes and coastal arrival times were strongly dependent on slip distribution. The cumulative probability distributions of the variations of both coastal tsunami amplitude and arrival time among the scenarios were well approximated by Gaussian distributions, thus indicating that if slip distributions are unknown the average tsunami amplitude could be used as a first approximation forecast of tsunami height and the standard deviation as its uncertainty. For a uniform slip distribution on the fault, at most observation sites the estimates of peak tsunami amplitudes were underestimated and the arrival times were earlier than those for which a heterogeneous slip distribution scenario was used. Our results highlight the importance of the use of heterogeneous slip distributions in estimations of tsunami amplitude, but also showed that the use of uniform slip distributions may be useful for calculating conservative arrival times for tsunami early warnings.

The detailed statistical characteristics of tsunami amplitudes and arrival times clearly depend on assumptions about earthquake magnitude, source region, and the probability density function used to generate slip distributions. Therefore, a comprehensive ensemble of scenarios covering a wide range of earthquake magnitudes and source regions are necessary for accurate evaluations of tsunami hazards and effective tsunami early warning systems.

JIVSM: Japan Integrated Velocity Structure Model; JMA: Japan Meteorological Agency; OBP: ocean-bottom pressure gauge; PDF: probability density function; PNTA: peak near-shore tsunami amplitude; SDR: slip deficit rate; SPDF: slip probability density function

Availability of data and material

The code to compute slip distribution on a non-planar fault surface with seafloor rupture is available at https://github.com/s-murfy/k223d. The JAGURS code to compute tsunami wave propagations is available at https://github.com/jagurs-admin/jagurs. We used a M7000 Digital Bathymetric Chart provided by the Japan Hydrographic Association (https://www.jha.or.jp/en/jha/) for our modelling of tsunami wave propagation. We used coastline data from the National Land Numerical Information provided by the Japanese Ministry of Land, Infrastructure, Transport and Tourism (https://nlftp.mlit.go.jp/ksj/index.html).

Competing interests

The authors declare that they have no competing interest.

Funding

This work was supported by the CREST project “Establishing the most advanced disaster reduction management system by fusion of real-time disaster simulation and big data assimilation” of the Japan Science and Technology Agency.

Authors' contributions

MN generated the scenario earthquake, computed the tsunami models, and analyzed the results. SM developed the computation code for the stochastic slip distributions. RA developed the plate geometry model along the Nankai trough. YI and MO designed the research plan. HT proposed the research concept and organized the project. All authors read and approved the final manuscript.

Authors' information

Acknowledgements

We thank Dr. T. Iinuma for sharing the FORTRAN version of the algorithm of Meade (2007). All figures were drawn using Generic Mapping Tools (Wessel and Smith 1998). We thank Dr. T. Baba for careful review and constructive comments that allowed us to improve the manuscript.

Ando M (1975) Source mechanisms and tectonic significance of historical earthquakes along the Nankai trough, Japan. Tectonophysics 27:119–140.

Baba T, Takahashi N, Kaneda Y (2014a) Near-field tsunami amplification factors in the Kii Peninsula, Japan for Dense Oceanfloor Network for Earthquakes and Tsunamis (DONET). Mar Geophys Res 35:319–325. doi:10.1007/s11001-013-9189-1

Baba T, Takahashi N, Kaneda Y, Inazawa Y, Kikkojin M (2014b) Tsunami Inundation Modeling of the 2011 Tohoku earthquake using three-dimensional building data for Sendai, Miyagi Prefecture, Japan. In: Kontar YA, Santiago-Fandiño V, Takahashi T (eds) Tsunami Events and Lessons Learned, Advances in Natural and Technological Hazards Research. Vol. 35. Dordrecht: Springer, pp. 89–98.

Blewitt G, Hammond WC, Kreemer C, Plag HP, Stein S, Okal E (2009) GPS for real-time earthquake source determination and tsunami warning systems. J Geod 83:335–343. doi:10.1007/s00190-008-0262-5

Causse M, Cotton F, Mai PM (2010) Constraining the roughness degree of slip heterogeneity. J Geophys Res 115:B05304. doi:10.1029/2009JB006747

Earthquake Research Committee, the Headquarters for Earthquake Research Promotion in Japan (2013) Evaluations of Long-Term Seismic Activity along the Nankai Trough, Second Ed. available at https://www.jishin.go.jp/main/chousa/13may_nankai/nankai_gaiyou.pdf. Accessed 26 August 2019 (in Japanese).

Earthquake Research Committee, the Headquarters for Earthquake Research Promotion in Japan (2020) Probabilistic Hazard Assessment of Tsunami due to Large Earthquakes along the Nankai Trough. available at http://www.jishin.go.jp/main/chousa/20jan_tsunami/nankai_tsunami.pdf. Accessed 24 January 2020 (in Japanese).

Eshelby J (1957) The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc Roy Soc London Ser A 241:376–396. doi: 10.1098/rspa.1957.0133

Frankel A (1991) High-frequency spectral falloff of earthquakes, fractal dimension of complex rupture, b value, and the scaling of strength on faults. J Geophys Res 96:6291–6302.

Geist EL (2002) Complex earthquake rupture and local tsunamis. J Geophys Res 107:2086. doi:10.1029/2000JB000139

Hashimoto C, Noda A, Sagiya T, Matsu’ura M (2009) Interplate seismogenic zones along the Kuril–Japan trench inferred from GPS data inversion. Nature Geosci 2:141–144. doi:10.1038/NGEO421

Hashimoto C, Noda A, Matsu’ura M (2012) The Mw 9.0 northeast Japan earthquake: total rupture of a basement asperity. Geophys J Int 189:1–5. doi:10.1111/j.1365-246X.2011.05368.x

Herrero A, Bernard P (1994) A kinematic self-similar rupture process for earthquakes. Bull Seism Soc Am 84:1216–1228

Herrero A, Murphy S (2018) Self-similar slip distributions on irregular shaped faults. Geophy J Int 213:2060–2070. doi:10.1093/gji/ggy104

Ide S, Shiomi K, Mochizuki K, Tonegawa T, Kimura G (2010) Split Philippine Sea plate beneath Japan. Geophys Res Lett 37:L21304. doi:10.1029/2010GL044585

Igarashi Y, Hori T, Murata S, Sato K, Baba T, Okada M (2016) Maximum tsunami height prediction using pressure gauge data by a Gaussian process at Owase in the Kii Peninsula, Japan. Mar Geophys Res 37:361–370. doi:10.1007/s11001-016-9286-z

Ishibashi K (2004) Status of historical seismology in Japan. Annals Geophys 47:339–368.

Ishibashi M, Baba T, Takahashi N, Imai K (2018) Social Implementation of Tsunami Prediction System on Wakayama by Using DONET Information. J JSNDS 37:125–142.

Kanazawa T, Uehira K, Mochizuki M, Shinbo T, Fujimoto H, Noguchi S, Kunugi T, Shiomi K, Aoi A, Matsumoto T, Sekiguchi S, Okada Y (2016) S-net project, cabled observation network for earthquakes and tsunamis. Abstract WE2B–3 presented at SubOptic 2016, Suboptic, Dubai, 18–21 Apr.

Kaneda Y, Kawaguchi K, Araki E, Matsumoto H, Nakamura T, Kamiya S, Ariyoshi K, Hori T, Baba T, Takahashi N (2015) Development and application of an advanced ocean floor network system for megathrust earthquakes and tsunamis. In Favali P, Beranzoli L, De Santis A (eds) Seafloor Observatories. Springer, Heidelberg, pp 643–662. doi:10.1007/978-3-642-11374-1_252

Kawaguchi K, Kaneko S, Nishida T, Komine T (2015) Construction of the DONET real-time seafloor observatory for earthquakes and tsunami monitoring, In Favali P, Beranzoli L, De Santis A (eds) Seafloor Observatories. Springer, Heidelberg, pp. 211–228. doi:10.1007/978-3-642-11374-1_10

Koketsu K, Miyake H, Suzuki H (2012) Japan Integrated Velocity Structure Model Version 1. In: Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September.

Konca AO, Avouac JP, Sladen A, Meltzner AJ, Sieh K, Fang P, Li Z, Galetzka J, Genrich J, Chlieh M, Natawidjaja DH, Bock Y, Fielding EJ, Ji C, Helmberger DV (2008) Partial rupture of a locked patch of the Sumatra megathrust during the 2007 earthquake sequence. Nature 456:631–635. doi:10.1038/nature07572

Li L, Switzer AD, Chan CH, Wang Y, Weiss R, Qiu Q (2016). How heterogeneous coseismic slip affects regional probabilistic tsunami hazard assessment: A case study in the South China Sea, J Geophys Res 121:6250–6272. doi:10.1002/2016JB013111

Maeda T, Obara K, Shinohara M, Kanazawa T, Uehira K (2015) Successive estimation of a tsunami wavefield without earthquake source data: A data assimilation approach toward real-time tsunami forecasting. Geophys Res Lett 42:7923–7932. doi:10.1002/2015GL065588

Mai PM, Beroza GC (2002) A spatial random field model to characterize complexity in earthquake slip. J Geophys Res 107:2308. doi:10.1029/2001JB000588

Meade BJ (2007) Algorithms for the calculation of exact displacements, strains, and stresses for triangular dislocation elements in a uniform elastic half space. Comp Geosci 33:1064–1075. doi:10.1016/j.cageo.2006.12.003

Melger DM, Williamson AL, Salazar-Monroy EF (2019) Differences between heterogeneous and homogenous slip in regional tsunami hazards modelling. Geophys J Int 219:553–562. doi: 10.1093/gji/ggz299

Mochizuki M, Kanazawa T, Uehira K, Shimbo T, Shiomi K, Kunugi T, Aoi S, Matsumoto T, Sekiguchi S, Yamamoto-Chikasada N, Takahashi N, Shinohara M, Yamada T (2016). S-net project: Construction of large scale seafloor observatory network for tsunamis and earthquakes in Japan. Abstract NH43B-1840 presented at 2016 AGU Fall Meeting, American Geophysical Union, San Francisco, Calif., 12–16 Dec.

Moore GF, Bangs NL, Taira A, Kuramoto S, Pangborn E, Tobin HJ (2007) Three-Dimensional Splay Fault Geometry and Implications for Tsunami Generation. Science 318:1128–1131. doi:10.1126/science.1147195

Moore GF, Park JO, Bangs NL, Gulick SP, Tobin HJ, Nakamura Y, Sato S, Tsuji T, Yoro T, Tanaka H, Uraki S, Kido Y, Sanada Y, Kuramoto S, Taira A (2009) Structural and seismic stratigraphic framework of the NanTroSEIZE Stage 1 transect. In Kinoshita M, Tobin H, Ashi J, Kimura G, Lallemant S, Screaton EJ, Curewitz D, Masago H, Moe KT, the Expedition 314/315/316 Scientists, Proc. IODP, 314/315/316: Washington, DC (Integrated Ocean Drilling Program Management International, Inc.). doi:10.2204/iodp.proc.314315316.102.2009

Moreno M, Rosenau M, Oncken O (2010) 2010 Maule earthquake slip correlates with

pre-seismic locking of Andean subduction zone. Nature 467:198–202. doi:10.1038/nature09349

Mueller C, Power W, Fraser S, Wang X (2015) Effects of rupture complexity on local tsunami inundation: Implications for probabilistic tsunami hazard assessment by example. J Geophys Res 120:488–502. doi:10.1002/2014JB011301

Murphy S, Scala A, Herrero A, Lorito S, Festa G, Trasatti E, Tonini R, Romano F, Molinari I, Nielsen S (2016) Shallow slip amplification and enhanced tsunami hazard unravelled by dynamic simulations of mega-thrust earthquakes. Sci Rep 6:35007. doi:10.1038/srep35007

Murphy S, Herrero A, (2020) Surface rupture in stochastic slip models. Geophys J Int 221:1081–1089. doi: 10.1093/gji/ggaa055

Nishimura T, Yokota Y, Tadokoro K, Ochi T (2018) Strain partitioning and interplate coupling along the northern margin of the Philippine Sea plate, estimated from Global Navigation Satellite System and Global Positioning System‑Acoustic data. Geosphere 14:535–551. doi:10.1130/GES01529.1

Noda A, Saito T, Fukuyama E (2018) Slip-deficit rate distribution along the Nankai trough, Southwest Japan, with elastic lithosphere and viscoelastic asthenosphere. J Geophys Res 123:8125–8142. doi:10.1029/2018JB015515

Park JO, Tsuru T, Kodaira S, Cummins PR, Kaneda Y (2002) Splay Fault Branching Along the Nankai Subduction Zone. Science 297:1157–1160. doi:10.1126/science.1074111.

Perfettini H, Avouac JP, Tavera H, Kositsky A, Nocquet JM, Bondoux F, Chlieh M, Sladen A, Audin L, Farber DL, Soler P (2010) Seismic and aseismic slip on the Central Peru megathrust. Nature 465:78–81. doi:10.1038/nature09062

Ruiz JA, Baumont D, Bernard P, Berge-Thierry C (2011) Modelling directivity of strong ground motion with a fractal, k−2, kinematic source model, Geophys J Int 186:226–244. doi:10.1111/j.1365-246X.2011.05000.x

Scala A, Lorito S, Romano F, Murphy S, Selva J, Basili R, Babeyko A, Herrero A, Hoechner A, Løvholt F, Maesano FE, Perfetti P, Tiberti MM, Tonini R, Volpe M, Davies G, Festa G, Power W, Piatanesi A, Cirella A (2020) Effect of shallow slip amplification uncertainty on probabilistic tsunami hazard analysis in subduction zones: Use of long-term balanced stochastic slip models. Pure Appl Geophys 177:1497–1520. doi:10.1007/s00024-019-02260-x

Smith WHF, Wessel P (1990) Gridding with continuous curvature splines in tension. Geophysics 55:293–305.

Strasser FO, Arango MC, Bommer JJ (2010) Scaling of the source dimensions of interface and intraslab Subduction-zone earthquakes with moment magnitude. Seismo Res Lett 81:941–950. doi:10.1785/gssrl.81.6.941

Takahashi N, Imai K, Ishibashi M, Sueki K, Obayashi R, Tanabe T, Tamazawa F, Baba T, Kaneda Y (2017) Real-Time Tsunami Prediction System Using DONET. J Disast Res 12:766–774.

Takahashi N, Imai K, Sueki K, Obayashi R, Ishibashi M, Tanabe T, Baba T, Kaneda Y (2018) Real-Time Tsunami Prediction System Based on Seafloor Observatory Data Applied to the Inland Sea, Japan. Mar Tech Soc J 52:120–127.

Tatehata H (1997) The new tsunami warning system of the Japan Meteorological Agency. In Perspectives on Tsunami Hazard Reduction, edited by Hebenstreit G, pp. 175–188, Springer, New York.

Titof VV, González FI, Bernard EN, Eble MC, Mofjeld HO, Newman JC, Venturato AJ (2005) Real-time tsunami forecasting: Challenges and solutions. Nat Hazards 35:41–58.

Tsushima H, Hino R, Fujimoto H, Tanioka Y, Imamura F (2009) Near-field tsunami forecasting from cabled ocean bottom pressure data. J Geophys Res 114:B06309. doi:10.1029/2008JB005988

Tsushima H, Hirata K, Hayashi Y, Tanioka Y, Kimura K, Sakai S, Shinohara M, Kanazawa Y, Hino R, Maeda K (2011) Near-field tsunami forecasting using offshore tsunami data from the 2011 off the Pacific coast of Tohoku Earthquake. Earth Planet Space 63:821–826. doi:10.5047/eps.2011.06.052

Tsushima H, Hino R, Tanioka Y, Imamura F, Fujimoto H (2012) Tsunami waveform inversion incorporating permanent seafloor deformation and its application to tsunami forecasting. J Geophys Res 117:B03311. doi:10.1029/2011JB008877

Tsushima H, Hino R, Ohta Y, Iinuma T, Miura S (2014). tFISH/RAPiD: Rapid improvement of near-field tsunami forecasting based on offshore tsunami data by incorporating onshore GNSS data. Geophys Res Lett 41:3390–3397. doi:10.1002/2014GL059863

Uehira K, Kanazawa T, Mochizuki M, Fujimoto H, Noguchi S, Shinbo T, Shiomi K, Kunugi T, Aoi S, Matsumoto T, Sekiguchi S, Okada Y, Shinohara M, Yamada T (2016) Outline of seafloor observation network for earthquakes and tsunamis along the Japan Trench (S-net), Abstract EGU2016-13832 presented at EGU General Assembly 2016, European Geosciences Union, Vienna, Austria, 17–22 Apr.

Wessel P, Smith WHF (1998) New, Improved Version of Generic Mapping Tools Released. EOS 79:579.

Yamamoto N, Aoi S, Hirata K, Suzuki W, Kunugi T, Nakamura H (2016a) Multi‑index method using offshore ocean‑bottom pressure data for real‑time tsunami forecast. Earth Planet Space 68:128. doi:10.1186/s40623-016-0500-7

Yamamoto N, Hirata K, Aoi S, Suzuki W, Nakamura H, Kunugi T (2016b) Rapid estimation of tsunami source centroid location using a dense offshore observation network. Geophys Res Lett. 43:4263–4269. doi:10.1002/2016GL068169

Yokota Y, Ishikawa T, Watanabe S, Tashiro T, Asada A (2016) Seafloor geodetic constrains on interplate coupling of the Nankai Trough megathrust zone, Nature 534:374–377. doi:10.1038/nature17632.

Zeng Y, Anderson JG, Yu G (1994) A composite source model for computing realistic synthetic strong ground motions. Geophys Res Lett 21:725–728.

supplement15.avi
{"primaryId":"undefined","secondaryId":"PEPS-D-20-00058","acdcId":"undefined","revision":"undefined","timestamp":"2020-06-02T14:28:35","document":"supplements","linkRel":"supplement"}
supplement16.avi
{"primaryId":"undefined","secondaryId":"PEPS-D-20-00058","acdcId":"undefined","revision":"undefined","timestamp":"2020-06-02T14:28:35","document":"supplements","linkRel":"supplement"}
supplement17.avi
{"primaryId":"undefined","secondaryId":"PEPS-D-20-00058","acdcId":"undefined","revision":"undefined","timestamp":"2020-06-02T14:28:35","document":"supplements","linkRel":"supplement"}

Download PDF

Journal Publication

published 31 Aug, 2020

Read the published version in Progress in Earth and Planetary Science →

Editorial decision: Minor revision
05 Jul, 2020
Review #1 received at journal
03 Jul, 2020
Review #2 received at journal
21 Jun, 2020
Reviewer #2 agreed at journal
14 Jun, 2020
Reviewers invited by journal
05 Jun, 2020
Reviewer #1 agreed at journal
05 Jun, 2020
Editor assigned by journal
02 Jun, 2020
Submission checks completed at journal
01 Jun, 2020
Editor invited by journal
01 Jun, 2020
First submitted to journal
28 May, 2020

You are reading this older preprint version

Read the latest preprint version →

Self-similar stochastic slip distributions on a non-planar fault for tsunami scenarios for megathrust earthquakes

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Methods/experimental

Results And Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1