Next Article in Journal
Meteorological Patterns Linked to Landslide Triggering in Asturias (NW Spain): A Preliminary Analysis
Next Article in Special Issue
Assessment of Educational Methods for Improving Children’s Awareness of Tsunamis and Other Natural Disasters: Focusing on Changes in Awareness and Regional Characteristics in Japan
Previous Article in Journal
Earthquake Magnitude and Shaking Intensity Dependent Fragility Functions for Rapid Risk Assessment of Buildings
Previous Article in Special Issue
Development of a Tsunami Inundation Analysis Model for Urban Areas Using a Porous Body Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Assessment of Epistemic Uncertainties in Tsunami Hazard Effects on Building Risk Assessments

1
College of Science and Engineering, Kanto Gakuin University, Yokohama 236-8501, Japan
2
International Research Institute of Disaster Science, Tohoku University, Sendai 980-0845, Japan
*
Author to whom correspondence should be addressed.
Geosciences 2018, 8(1), 17; https://doi.org/10.3390/geosciences8010017
Submission received: 30 September 2017 / Revised: 20 December 2017 / Accepted: 5 January 2018 / Published: 10 January 2018
(This article belongs to the Special Issue Interdisciplinary Geosciences Perspectives of Tsunami)

Abstract

:
Based on the definition of tsunami risk, we quantitatively evaluated the annual expected tsunami loss ratio (Tsunami Risk Index) and clarified the quantitative effects of epistemic uncertainties in tsunami hazard assessments on the tsunami risk of buildings by combining probabilistic information regarding tsunami inundation depths at target points and tsunami fragility assessments of buildings. For the risk assessment, we targeted buildings with four different structures (reinforced concrete, steel, brick, wood) located in three different areas (Soma, Sendai, Kesennuma). In conclusion, we demonstrated that the expected tsunami risk could vary by approximately two orders of magnitude when considering tsunami hazard uncertainties between the 95th percentile and the 5th percentile. In addition, we quantitatively clarified the fact that we cannot properly understand the tsunami risk by evaluating the tsunami fragility alone. For example, the analysis results indicate that the tsunami risk of a wood building located in Kesennuma is lower than that of a reinforced concrete building located in either Soma or Sendai.

1. Introduction

In general, the risks of natural disasters can be understood by multiplying the magnitude of the hazard and the evaluated vulnerability [1]. The magnitude of the hazard must include the severity of the hazard and its occurrence probability. Therefore, it is necessary to carry out probabilistic natural hazard assessments in order to properly evaluate such risks. Based on this concept, tsunami risks can be understood by multiplying the probabilistic tsunami hazard assessment and the tsunami fragility assessment, which represents the vulnerability of the tsunami. In other words, neither the tsunami hazard nor the tsunami fragility alone can capture the tsunami risk for an appropriate evaluation. Numerous studies on probabilistic tsunami hazard assessment methods and tsunami fragility evaluation methods have been conducted individually but few studies have quantitatively assessed tsunami risks by multiplying them together. Moreover, probability estimations of tsunami inundation depths at inland locations constitute necessary input information for the quantitative evaluation of tsunami risks; however, not many studies have conducted probabilistic tsunami inundation assessments (e.g., [2,3,4,5,6,7,8]). González et al. [2] applied a probabilistic tsunami hazard assessment methodology to Seaside, Oregon, that combines tsunami inundation modeling with probabilistic concepts. Goda et al. [4] developed stochastic random-field slip models for the 2011 Tohoku earthquake and conducted a tsunami inundation simulation using those developed models, after which they concluded that the inundation heights at a coastal location are major sources of uncertainties in the prediction of tsunami risks. Park and Cox [5] demonstrated an approach for assessing probabilistic near-field tsunami hazards within the Cascadia Subduction Zone using the logic tree method.
The logic tree method has been widely utilized to conduct probabilistic tsunami hazard analysis (e.g., [5,6,7,8,9,10]). In this method, the uncertainties in tsunami hazard assessments are classified as either epistemic uncertainties or aleatory uncertainties. Historically, Cornell [11] categorized the types of uncertainties that should be considered during an evaluation of strong earthquake motion into epistemic uncertainties and aleatory uncertainties. It has since become commonplace to evaluate the epistemic and aleatory uncertainties in probabilistic tsunami hazard assessments. As shown in Figure 1, the epistemic uncertainty can be captured using a logic tree, while the aleatory uncertainty can be evaluated through a probability distribution function of the tsunami height.
In this study, we first demonstrate a method for evaluating the probabilistic tsunami inundation depths at inland locations by employing a probabilistic tsunami hazard assessment model using the logic tree method. The demonstrated probabilistic tsunami inundation depths are combined with information regarding the tsunami fragility of buildings and is therefore capable of quantitatively evaluating the tsunami risk of a specific building located at a specific location. If the tsunami risk can be evaluated quantitatively, it is possible to objectively compare the risks among buildings with different structures in different areas. Since the method proposed in this study is relatively simple, its calculation cost is smaller than those of the methods employed within previous studies.
After the 2011 Great East Japan Tsunami, many improvements of tsunami countermeasures have been made in both structural measures and nonstructural measures but large future challenges remain such as probabilistic risk prediction and an advanced simulation technique and system for real-time hazard etc. [12]. The method proposed in this study will be used to produce tsunami risk index using objective judgment for various situations in disaster prevention scenarios. Consequently, we can also examine how the evaluation of the tsunami risk index varies depending on the structures of buildings and other regional differences. The method is intended to quantitatively clarify the contribution from the epistemic uncertainty in the tsunami hazard evaluation to the quantitative amount of risk.

2. Methodology for the Tsunami Risk Assessment and Quantification of Hazard Uncertainty

Figure 2 shows the workflow for the quantitative assessment of the risk posed by a tsunami of the direct damage of a building located at a specific location in consideration of tsunami hazard uncertainties. First, the tsunami wave height ( h ) simulated from each fault source exceeding a specific level ( h t h ) is determined as follows:
q i j k ( h t h ) = h t h p i j k ( h ) d h
where i is the number of assumed fault sources in one earthquake region, j is the number of logic tree branches, k is the number of assumed earthquake regions and p i j k ( h ) is the aleatory uncertainty, which is defined as a lognormal normal distribution. Then, the annual exceedance frequency of the tsunami wave height in one earthquake region ( Q j k ( h t h ) ) is calculated as follows using the annual exceedance frequency of each fault source μ i :
Q j k ( h t h ) = i q i j k ( h t h ) × μ i
If we assume a Poissonian arrival time distribution, the annual exceedance probability of the tsunami wave height in one earthquake region ( f j k ( h t h ) ) is given using the following exponential function:
f j k ( h t h ) = 1 e x p ( Q j k ( h t h ) )
By considering the annual excess probability of the tsunami wave height created for each logic tree branch and the weight of each logic tree branch ( w j ) , the annual exceedance probability of the tsunami wave height for each percentile (perc) in each earthquake region ( F k ,   p e r c ( h t h ) ) can be obtained as follows:
F k ,   p e r c ( h t h ) = f j k ( h t h ,   w j )
After calculating these curves for each earthquake region, we can obtain the annual exceedance probability of the tsunami wave height for each percentile (perc) for the entire earthquake region ( P p e r c ( h t h ) ) through the following formula (Figure 2a):
P p e r c ( h t h ) = 1 k ( 1 F k ,   p e r c ( h t h ) )
We convert this annual exceedance probability of the tsunami wave height into the annual exceedance probability of the tsunami inundation depth using the results from the tsunami numerical simulation (Figure 2b). Meanwhile, assuming that the tsunami-induced damage level of each building is n, the damage probability for each damage level ( d n ( h ) ) can be written as follows (Figure 2c):
d n ( h ) = Φ ( ln ( h ) λ ζ ) n
where Φ ( · ) n represents the standardized normal distribution function for the damage levels and λ and ζ represent the mean and standard deviation of ln ( h ) , respectively. Then, by using the damage probability coefficient for each damage level ( a n ) , the damage probability function ( D ( h ) ) can be written as follows:
D ( h ) = l = 1 6 a l × ( d l ( h ) d l 1 ( h ) )
In this study, we consider six damage levels: minor damage, moderate damage, major damage, complete damage, collapsed damage and washed away. Consequently, by using the annual exceedance probability of the inundation depth (Equation (5)) and the damage probability (Equation (7)), the tsunami risk curve for each percentile hazard can be written as follows (Figure 2d):
R p e r c ( P p e r c ,   D )
We should note that we have not performed a tsunami risk estimation that includes variability in the building response and uncertainties in the loss estimates because the main objective of this study is to quantify the epistemic uncertainties in tsunami hazards affecting the building risk. Future research will include other uncertainties in addition to those of hazards within the risk assessment. Finally, we can quantitatively evaluate the tsunami risk by integrating Equation 8 based on the definition of the tsunami risk. We can consider this index representative of the annual expected tsunami loss ratio (%/year), that is, a Tsunami Risk Index (TRI) (Figure 2e):
T s u n a m i   R i s k   I n d e x   p e r c = 0 1 R p e r c ( P p e r c ,   D ) d D

3. Application to the Tohoku Area

3.1. Assessment Targets

For the targets of the risk assessment, we considered wooden, brick, steel and reinforced concrete buildings located in the city of Soma in Fukushima Prefecture and the cities of Sendai and Kesennuma in Miyagi Prefecture. Figure 3 and Table 1 show the information for the locations of the buildings in each city. Kesennuma is in the ria coast, Sendai is inside of the Sendai bay and Soma is directly facing the Pacific Ocean. We selected these three regions to capture the difference of the tsunami risk due to different geography. The elevation and distance from coastline of the three cities also indicated in Table 1.

3.2. Probabilistic Tsunami Wave Height

3.2.1. Construction of the Logic Trees

In this section, we aim to evaluate the probabilistic coastal tsunami height in the Tohoku region adjacent to the Pacific Ocean. First, we need to select an appropriate earthquake-generating fault to produce a tsunami. We select ten regions and one region of interlocking earthquakes from among the occurrence areas of trench-type earthquakes along the Japan Trench used in the probabilistic earthquake prediction map released by the NIED [14], all of which are shown in Figure 4 and Table 2 as earthquake faults that could generate tsunamis. These selections exclude both earthquakes with moment magnitudes (Mw) reaching 7.4 or less in consideration of a variation of ±0.1 Mw and earthquakes for which the source fault is unlikely to be predicted beforehand. To evaluate the epistemic uncertainties for these eleven fault regions, we use the logic tree method proposed in Annaka’s study [10]. Figure 5 shows the logic tree constructed for these regions. We establish five branches within the logic tree: The Mw range of the earthquake, the asperity position of the earthquake fault, the average occurrence interval (return period) of the earthquake, the standard deviation of the lognormal distribution followed by the error of the tsunami wave height and the truncation range of the lognormal distribution. Except for the asperity position of the earthquake fault, the other four branches follow the branches shown in Annaka’s study [10]. The total number of branches in the logic tree constructed using this approach is 3384 branches. The outline of setting for each branch is as follows.
The Mw range of each earthquake is varied by ±0.1. This Mw variation is accomplished by changing the average slip amount along the entire fault. For the asperity position of the fault, three types of branches with asperities located at the center of the fault and near both ends of the fault are established when the fault length is 150 km or more. Since only Tohoku-type earthquakes have long lengths (approximately 500 km), we set five branches by adding two cases where an additional asperity is located between the three asperities. The method utilized to establish the asperity positions along the fault is detailed in Fukutani et al. [15]. Regarding the occurrence probability of the earthquake, we construct three types of branches that take into consideration the confidence interval of the occurrence probability determined by the probabilistic seismic motion prediction evaluation published by the NIED [14]. Table 3 shows the model name for the generation interval of the earthquake, α value of the BPT distribution, average return period, sample period, earthquake generation time within the period used to determine the average occurrence interval and lower and upper limits of the confidence interval for each earthquake fault, which are shown in the Headquarters for Earthquake Research Promotion [16]. See Appendix A for the method used to establish the confidence intervals. Although the probability of the occurrence of an earthquake with a relatively small magnitude can be evaluated with high accuracy using the general Gutenberg-Richter rule, it is known that the probability of occurrence of a relatively high-magnitude earthquake that causes a large tsunami cannot be perfectly evaluated using the Gutenberg-Richter rule. Based upon this background, in this study, we note that data of the return period collected as the result of a detailed examination of the historical earthquake record are used. Regarding the standard deviation of the lognormal distribution (i.e. aleatory uncertainty) followed by the error of the tsunami wave height, we use the modeling error in the tsunami numerical simulation results and observation records, that is, the geometrical standard deviation κ of Aida [17], for the past eleven historical earthquakes represented by the following expression based on Annaka’s study [10]:
l o g β = i l o g β i
l o g κ = 1 n i ( l o g β i ) 2 ( l o g β ) 2
where n is the number of observation points, i is the observation point, βi = (Ri/Hi), Ri is the observed tsunami height at the i-th point and Hi is the simulated value at the i-th point. The values of κ are evaluated from the past eleven historical earthquakes. The minimum value is the result for the 1707 Hoei earthquake, where κ = 1.35 (σ = log κ = 0.300) and the maximum value is the result for the 1946 Nankai earthquake, where κ = 1.60 (σ = log κ = 0.470). Finally, two types of branches with ±2.3σ and ±10σ are established as truncation values at both ends of the lognormal distribution.
The numbers attached to the branches of the logic trees are the weights of each branch and the sum of all of the weights is set to 1.0. The weights of the branches containing the Mw range and asperity positions along the fault are set by equally dividing their weights. Regarding the branch consisting of the occurrence probability of the earthquakes, we use a weight of 0.50 for the central branch and 0.25 for the branches at both ends to consider the confidence interval. The weight value adopted in Annaka’s study [10] is also adopted for the weights of the branches for the standard deviation of the lognormal distribution and the truncation range.

3.2.2. Tsunami Numerical Simulation

For each branch comprising either the Mw range or asperity positions in the logic trees constructed in the previous section, the fault parameters for the tsunami numerical simulation are determined. The fault parameters for the reference magnitudes of each fault are shown in Table 4. These data, which are related to the position of the fault, were published by the NIED [14]. To provide heterogeneous asperities along the fault planes, we generated 10 km mesh points covering the Japan Trench area shown in Figure 6 while assuming the presence of small faults with lengths and widths of 10 km and setting a slip amount to each fault. As described above, the method used to establish the asperity along the 3.11 Tohoku earthquake-type fault is described in Fukutani et al. [15]. The depth of each small fault was set considering its length, width and dip estimated from the depth along the entire fault. The strike, dip and rake of each small fault were set to the same values as those along the entire fault. The calculation conditions of the tsunami numerical simulation are shown in Table 5 and the terrain data used for the simulation are shown in Figure 7. The terrain data were generated from 30 s gridded depth data (J-TOPO 30, Japan Hydrographic Association) near the islands of Japan and a mesh of 500 m bathymetry data (J-EGG 500, Japan Oceanographic Data Center). We obtained the initial displaced water height using the formula of Okada [18] from these fault parameters. Using the calculated initially displaced water heights as input values, time integration was performed in each mesh using the continuous equation (Equation (12)) and the equations of motion (Equations (13) and (14)) based on a nonlinear longwave equation using the TUNAMI model (Tohoku University Numerical Analysis Modeling for Inundation) [19]:
η t + M x + N y = 0
M t + x ( M 2 D ) + y ( M N D ) + g D η x + g n 2 M D 7 / 3 M 2 + N 2 = 0
N t + x ( M N D ) + y ( N 2 D ) + g D η y + g n 2 N D 7 / 3 M 2 + N 2 = 0
where η is the water level, D is the total water depth, g is the gravitational acceleration, n is Manning’s roughness coefficient and M and N are the flow fluxes in the x and y directions. Taking into consideration the number of branches comprising the Mw range and asperity positions along the earthquake faults, there are 258 cases in which the tsunami numerical simulation can be performed.

3.2.3. Tsunami Hazard Curves at the Offshore Points

For each branch of the constructed logic tree, the maximum wave height in each mesh is determined according to the results of the tsunami numerical simulation. If we assume the simulated maximum wave height is the median value μ and if we use the lognormal standard deviation σ in the logic tree, we can obtain a probability density function of the tsunami wave height represented by the following (Equation (15)):
f ( x ) = 1 2 π σ x e x p { ( l o g x μ ) 2 2 σ 2 } ,   0 x
Next, by converting the probability density function into the exceedance probability distribution under an ergodic assumption, it is possible to obtain the tsunami hazard curve expressed by the relationship between the tsunami wave height and the annual exceedance probability for each branch of the logic tree. The ergodic assumption is a statistical assumption that the spatial variation is equal to the temporal variation. By evaluating the annual exceedance probability distribution within an earthquake area and drawing curves along percentile paths in consideration of the weights of the logic tree branches, it is possible to estimate the tsunami hazard curve (i.e. a fractile curve) for that earthquake area. Finally, we integrate each tsunami hazard curve for each earthquake area. The evaluation results for 10 m water depths off Soma, Sendai and Kesennuma are shown in Figure 8.

3.3. Tsunami Inundation Assessment

3.3.1. Probabilistic Tsunami Hazard Map

The results shown in Table 6 for Soma, Table 7 for Sendai and Table 8 for Kesennuma represent the values converted from the tsunami wave heights at water depths of 10 m generated for each earthquake to the annual exceedance probability and the return period using each tsunami hazard curve. It should be noted that we used the average values of the tsunami hazard curves in this section. Using these tables, we can consider the possibility of specifying the return period of an earthquake-generating fault by focusing on the tsunami height in a coastal area. Then, if a tsunami run-up simulation is carried out using the parameters of the fault specifying the return period, a tsunami inundation area and a tsunami inundation depth for every return period on land are obtained by running a tsunami simulation. The tsunami run-up simulations were conducted using the nonlinear longwave equations under the calculation conditions shown in Table 9 (Soma) and Table 10 (Sendai and Kesennuma). The calculation regions for each area are shown in Figure 9 (Soma) and Figure 10 (Sendai and Kesennuma). The tsunami inundation areas for each return period simulated using the above method are shown in Figure 11 and the tsunami inundation heights and tsunami inundation depths at each target point within the risk assessment are shown in Table 6 through Table 8 (g), (h). This method of calculating the inundation area for each return period is advantageous because it is possible to easily perform the inundation calculations by appropriately changing the information regarding the initial tide levels or artificial structures (e.g., dikes and buildings).

3.3.2. Tsunami Hazard Curves at Inland Points

Based on the annual exceedance probability of the tsunami wave height and the numerical simulation results of the inundation depths, we evaluate the annual exceedance probability curves of inland inundation depths.
First, in each case of the tsunami numerical simulation, we assumed that the probability density of the tsunami wave height was equal to the probability density of the inundation depth. Then, after normalizing the probability density of the inundation depth so that the maximum value of the annual exceedance probability of the inundation depth represents the annual exceedance probability of the tsunami wave height in the case where the inundation depth is zero, we calculated the annual exceedance probability of the inundation depth from the corresponding probability density data (see Figure 12). From Table 6, Table 7 and Table 8, the tsunami wave height (in the case where the inundation depth at each evaluation point is zero) and the annual exceedance probability are respectively 2.01 m and 0.007986 (average value) for Soma, 4.30 m and 0.006729 (average value) for Sendai and 3.63 m and 0.004062 (average value) for Kesennuma. In this study, although the procedure is not entirely probabilistic because the tsunami numerical simulations are performed only discretely with a limited number, we constructed the annual exceedance probability curve of the inundation depth by regressing each plot with a straight line. The average values in addition to the 5th percentile, 50th percentile and 95th percentile curves created using the same method are also shown in Figure 12.

3.4. Fragility Assessment

To evaluate the fragilities of buildings with respect to tsunamis, we used the fragility curve studied by Suppasri et al. [20], who performed regression analyses using damage data from buildings that suffered from tsunami inundation during the 3.11 Tohoku earthquake. Their study proposed various fragility curves for different building structures (i.e. reinforced concrete, steel, brick and wood) according to six different damage levels: minor damage, moderate damage, major damage, complete damage, collapsed damage and washed away (Figure 13). Although there are numerous studies about tsunami fragility functions [21], we selected the before-mentioned curves in this time. That is, why we also need to investigate uncertainties of tsunami risk due to the difference of tsunami fragility functions in the future study.

3.5. Risk Assessment and Quantitative Effects of the Hazard Assessment

To quantify the tsunami risk with regard to the direct damage of a building located at a risk assessment point, we eliminated the tsunami inundation depth axis from Figure 12 and Figure 13 and obtained the tsunami risk curve represented by the relationship between the damage probability of a building and the annual exceedance probability. The tsunami risk curve for each type of building structure is shown in Figure 14. In calculating the damage probability of a building, the six different types of destruction were united using each damage probability. We set the damage probabilities of minor damage to 0.1, moderate damage to 0.3, major damage to 0.5, complete damage to 0.8, collapsed damage to 0.9 and washed away to 1.0. As indicated in the abovementioned methodology, we note that we have not performed a tsunami risk estimation that includes variability in the building response and uncertainties in the loss estimates because the main objective of this study is to quantify epistemic uncertainties in the tsunami hazards affecting the building risk.
Generally, risk is expressed by the product of the degree of loss with its occurrence probability. Applying this idea, if we consider that the tsunami risk for direct damage to a building is the product of “the probability of destruction of a building by a tsunami” and its “generation probability”, we can estimate the tsunami risk using the product of the horizontal axis and the vertical axis of the risk curve. Therefore, we can quantify the tsunami risk by integrating the derived risk curve, calculating the area under the curve and estimating the annual expected tsunami loss ratio (i.e. the Tsunami Risk Index, or TRI). Figure 15 shows the calculation results obtained by integrating the area under the regressed exponential function for (a) the average value of the tsunami hazard and for each percentile value ((b) the 5th percentile value, (c) the 50th percentile value and (d) the 95th percentile value). Figure 16 shows the calculation results for each area.

3.6. Discussion

Based on the results shown in Figure 15, we can quantitatively understand the tsunami risk imposed on different types of buildings, that is, a reinforced concrete building located in Kesennuma is the safest (the lowest risk) while a wooden building located in Soma is the insecure (the highest risk) except for the case when average hazard values are used.
Focusing on the magnitude, we can see that the results using the 5th percentile hazard values and the results using the 50th percentile hazard values differ by approximately one order of magnitude and that the results using the 50th percentile values and the results using the 95th percentile values also differ by approximately one order of magnitude. A difference of approximately two orders of magnitude is observed between the 5th percentile values and the 95th percentile values. These results clearly show that tsunami risks have substantial hazard uncertainties.
Overall, the risk of a tsunami tends to increase successively from Soma to Sendai and then to Kesennuma because the impact of the tsunami hazard is large in that same order. In addition, a wood building located in Kesennuma has a lower tsunami risk than a reinforced concrete building located in either Soma or Sendai and a steel building located in Soma is at a higher tsunami risk than a wood building located in Kesennuma. We understand that we cannot properly comprehend the extent of the tsunami risk by evaluating the tsunami fragility only.
From the results shown in Figure 16, we can compare the results from evaluating the Tsunami Risk Index As mentioned previously, the tsunami risk changes dramatically depending on the percentile hazard that is considered. Among them, the difference between the tsunami risk using the 50th percentile hazard and the 5th percentile hazard is much greater than the difference between the tsunami risk using the 95th hazard and the 50th percentile hazard. As a matter of course, these results show that it is necessary to use hazard information from the same percentile when comparing the tsunami risks among different regions.
In this way, we can stochastically interpret and compare the evaluation results of the Tsunami Risk Index among several target regions.

4. Conclusions

In this study, through the process of quantifying the uncertainty in tsunami hazard effecting on building risk assessment, we proposed two new evaluation methods that are essential to implement tsunami risk assessment.
We first proposed a method that can be used to comprehensively address the uncertainties (epistemic uncertainty and aleatory uncertainty) in tsunami hazard assessments by probabilistically evaluating the tsunami inundation area and tsunami inundation depth using tsunami hazard curves. In this method, we can estimate the tsunami inundation area with each return period by performing tsunami numerical simulations using fault parameters after specifying the return period of the fault from the coastal tsunami hazard curve. This calculation method is advantageous because it is possible to easily perform the calculations by appropriately changing the information regarding the initial tide levels or artificial structures. In addition, we visualized the uncertainty in the tsunami hazard assessment in an easy-to-understand manner by indicating several tsunami inundation areas corresponding to the return periods of several earthquakes rather than using conventional tsunami hazard maps generated from single earthquakes.
Next, we proposed a method to evaluate the annual expected tsunami loss ratio (i.e. the tsunami risk index, or TRI) targeting buildings located in Soma, Sendai and Kesennuma by combining probabilistic tsunami inundation information and building fragility information. We used the fragility curves of different construction building derived from using the damaged data of buildings during the 3.11 Tohoku earthquake, which were regressed by log-normal distribution. In addition, then, we derived various tsunami risk curves, which represent the relationship between the damage probability of the building and the annual exceedance probability. Based on the definition of the risk, we considered the tsunami risk for direct damage to buildings as the product of the damage probability of a building destroyed by the tsunami and its generation probability. Therefore, we can quantify the tsunami risk by integrating the derived risk curve, that is, calculating the area under the curve and estimating the annual expected tsunami loss ratio (Tsunami Risk Index: TRI).
Focusing on the magnitude of the calculated TRI, we clearly showed that the results using the 5th percentile hazard values and the results using the 95th percentile hazard values differ by approximately two orders of magnitude. Furthermore, for example, based on the results using the average values, A wood building located in Kesennuma is at a lower tsunami risk than a reinforced concrete building located in Soma and Sendai. We cannot appropriately understand the extent of the tsunami risk if the information regarding tsunami hazards and tsunami fragility evaluations are individually captured. Through this study, we quantitatively showed for that we can understand the tsunami risk by combining hazard and fragility information. The proposed method in this study can be applied to various regions regardless of the area and therefore, the Tsunami Risk Index can be an effective index for ranking priority investments in disaster prevention endeavors by comparing the magnitudes of the tsunami risks of several targets located across several regions.

Acknowledgments

This research was supported by Specific Project Research from the International Research Institute of Disaster Science (IRIDeS) at Tohoku University. This research was also supported by funding from Tokio Marine & Nichido Fire Insurance Co., Ltd. through the IRIDeS. The authors also greatly appreciate the questions and valuable comments of the anonymous reviewers that helped improve the manuscript.

Author Contributions

Yo Fukutani conceived and designed the experiments, analyzed the data and wrote the paper; Anawat Suppasri and Fumihiko Imamura advised the experimental design.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

If we assume a Poisson process, the confidence interval of the earthquake occurrence interval is established based on the study of Weichert [22], as shown in Table A1. That is, if the number of records of the earthquake occurrence is N (N in the table), μU and μL are determined based on the information in Table A1. If the sample period is T, the confidence interval is as follows:
T μ U T μ L
Meanwhile, the confidence interval of the earthquake occurrence interval assuming an updating process using a BPT distribution is as follows:
exp ( α n )   ~   exp ( + α n )
where the variation coefficient of the earthquake occurrence interval is α and the number of records of the earthquake occurrence interval is n.
Table A1. Lower and upper ±1.0 standard deviation confidence intervals for a Poisson variable [22].
Table A1. Lower and upper ±1.0 standard deviation confidence intervals for a Poisson variable [22].
μ U N μ L
1.8400
3.3010.173
4.6420.708
5.9231.37
7.1642.09
8.3852.84
9.5863.62
10.874.42
12.085.23
13.196.06
14.3106.89

References

  1. United Nations Department of Humanitarian Affairs. Glossary: Internationally Agreed Glossary of Basic Terms Related to Disaster Management, Geneva, Switzerland. Available online: https://reliefweb.int/sites/reliefweb.int/files/resources/004DFD3E15B69A67C1256C4C006225C2-dha-glossary-1992.pdf (accessed on 28 September 2017).
  2. González, F.I.; Geist, E.L.; Jaffe, B.; Kânoğlu, U.; Mofjeld, H.; Synolakis, C.E.; Titov, V.V.; Arcas, D.; Bellomo, D.; Carlton, D.; et al. Probabilistic tsunami hazard assessment at seaside, Oregon, for near-and far-field seismic sources. J. Geophys. Res. Oceans 2009, 114, C11023. [Google Scholar] [CrossRef]
  3. Løvholt, F.; Pedersen, G.; Bazin, S.; Kuhn, D.; Bredesen, R.E.; Harbitz, C. Stochastic analysis of tsunami runup due to heterogeneous coseismic slip and dispersion. J. Geophys. Res. 2012, 117, C03047. [Google Scholar] [CrossRef]
  4. Goda, K.; Mai, P.M.; Yasuda, T.; Mori, N. Sensitivity of tsunami wave profiles and inundation simulations to earthquake slip and fault geometry for the 2011 Tohoku earthquake. Earth Planets Space 2014, 66, 105. [Google Scholar] [CrossRef] [Green Version]
  5. Park, H.; Cox, D.T. Probabilistic assessment of near-field tsunami hazards: Inundation depth, velocity, momentum flux, arrival time, and duration applied to Seaside, Oregon. Coast. Eng. 2016, 117, 79–96. [Google Scholar] [CrossRef]
  6. Thio, H.K.; Somerville, P.G.; Polet, J. Probabilistic Tsunami Hazard in California; College of Engineering, University of California: Los Angeles, CA, USA, 2010. [Google Scholar]
  7. Horspool, N.; Pranantyo, I.; Griffin, J.; Latief, H.; Natawidjaja, D.H.; Kongko, W.; Cipta, A.; Bustaman, B.; Anugrah, S.D.; Thio, H.K. A probabilistic tsunami hazard assessment for Indonesia. Nat. Hazards Earth Syst. Sci. 2014, 14, 3105–3122. [Google Scholar] [CrossRef]
  8. Davies, G.; Griffin, J.; Løvholt, F.; Glimsdal, S.; Harbitz, C.; Thio, H.K.; Lorito, S.; Basili, R.; Selva, J.; Geist, E.; et al. A global probabilistic tsunami hazard assessment from earthquake sources. Geol. Soc. Lond. Spec. Publ. 2017, 456, SP456-5. [Google Scholar] [CrossRef]
  9. Geist, E.L.; Parsons, T. Probabilistic analysis of tsunami hazards. Nat. Hazards 2006, 37, 277–314. [Google Scholar] [CrossRef]
  10. Annaka, T.; Satake, K.; Sakakiyama, T.; Yanagisawa, K.; Shuto, N. Logic-tree Approach for Probabilistic Tsunami Hazard Analysis and its Applications to the Japanese Coasts. Pure Appl. Geophys. 2007, 164, 577–592. [Google Scholar] [CrossRef]
  11. Cornell, C.A. Engineering Seismic Risk Analysis. Bull. Seism. Soc. Am. 1968, 58, 1583–1606. [Google Scholar]
  12. Suppasri, A.; Latcharote, P.; Bricker, J.D.; Leelawat, N.; Hayashi, A.; Yamashita, K.; Makinoshima, F.; Roeber, V.; Imamura, F. Improvement of Tsunami Countermeasures Based on Lessons from The 2011 Great East Japan Earthquake and Tsunami—Situation After Five Years. Coast. Eng. J. 2016, 58, 1640011. [Google Scholar] [CrossRef]
  13. Japan Society of Civil Engineers. A Method for Probabilistic Tsunami Hazard Analysis. The Tsunami Evaluation Subcommittee, The Nuclear Civil Engineering Committee, Japan Society of Civil Engineers: Tokyo, Japan. Available online: http://committees.jsce.or.jp/ceofnp/system/files/PTHA20111209_0.pdf (accessed on 28 September 2017). (In Japanese).
  14. National Research Institute for Earth Science and Disaster Resilience (NIED). Japan Seismic Hazard Information Station (JSHIS 2013). Available online: http://www.j-shis.bosai.go.jp/map/ (accessed on 28 September 2017).
  15. Fukutani, Y.; Suppasri, A.; Imamura, F. Stochastic analysis and uncertainty assessment of tsunami wave height using a random source parameter model that targets a Tohoku-type earthquake fault. Stoch. Environ. Res. Risk. Assess. 2015, 29, 1763–1779. [Google Scholar] [CrossRef]
  16. Headquarters for Earthquake Research Promotion, Probabilistic Seismic Hazard Maps for Japan. Available online: http://www.jishin.go.jp/evaluation/seismic_hazard_map/shm_report/ (accessed on 28 September 2017).
  17. Aida, I. Reliability of a tsunami source model derived from fault parameters. J. Phys. Earth 1978, 26, 57–73. [Google Scholar] [CrossRef]
  18. Okada, Y. Surface deformation due to shear and tensile faults in a half-space. Bull. Seism. Soc. Am. 1985, 75, 1135–1154. [Google Scholar]
  19. IUGG/IOC TIME Project. Numerical Method of Tsunami Simulation with the Leap-Frog Scheme; Intergovernmental Oceanographic Commission Manuals and Guides 35; UNESCO: Paris, France, 1997. [Google Scholar]
  20. Suppasri, A.; Mas, E.; Charvet, I.; Gunasekera, R.; Imai, K.; Fukutani, Y.; Abe, Y.; Imamura, F. Building damage characteristics based on surveyed data and fragility curves of the 2011 Great East Japan tsunami. Nat. Hazards 2013, 66, 319–341. [Google Scholar] [CrossRef]
  21. Charvet, I.; Macabuag, J.; Rossetto, T. Estimating Tsunami-Induced Building Damage through Fragility Functions: Critical Review and Research Needs. Front. Built Environ. 2017, 3, 36. [Google Scholar] [CrossRef]
  22. Weichert, D.H. Estimation of the earthquake recurrence parameters for unequal observation periods for different magnitudes. Bull. Seism. Soc. Am. 1980, 70, 1337–1356. [Google Scholar]
Figure 1. A diagram of the evaluation methods for the (a) epistemic uncertainty using a logic tree [10] and (b) aleatory uncertainty using a probabilistic density function [13].
Figure 1. A diagram of the evaluation methods for the (a) epistemic uncertainty using a logic tree [10] and (b) aleatory uncertainty using a probabilistic density function [13].
Geosciences 08 00017 g001
Figure 2. Flow of the tsunami risk quantification. (a) The probabilistic tsunami wave height is converted to (b) the probabilistic tsunami inundation information, and (d,e) the risk assessment is carried out by combining the hazard information with (c) fragility assessment.
Figure 2. Flow of the tsunami risk quantification. (a) The probabilistic tsunami wave height is converted to (b) the probabilistic tsunami inundation information, and (d,e) the risk assessment is carried out by combining the hazard information with (c) fragility assessment.
Geosciences 08 00017 g002
Figure 3. Locations of the target buildings.
Figure 3. Locations of the target buildings.
Geosciences 08 00017 g003
Figure 4. The selected earthquakes that could generate a tsunami based on the occurrence region of the trench-type earthquake along the Japan Trench, which is used in the probabilistic seismic motion prediction map [14].
Figure 4. The selected earthquakes that could generate a tsunami based on the occurrence region of the trench-type earthquake along the Japan Trench, which is used in the probabilistic seismic motion prediction map [14].
Geosciences 08 00017 g004
Figure 5. Logic trees constructed for the eleven fault regions. The numbers appended onto the branches of the logic trees are the weights of the branches.
Figure 5. Logic trees constructed for the eleven fault regions. The numbers appended onto the branches of the logic trees are the weights of the branches.
Geosciences 08 00017 g005aGeosciences 08 00017 g005b
Figure 6. Mesh consisting of 10 km points encompassing the Japan Trench area (black dots).
Figure 6. Mesh consisting of 10 km points encompassing the Japan Trench area (black dots).
Geosciences 08 00017 g006
Figure 7. Calculation area for the tsunami numerical simulation (450 m mesh).
Figure 7. Calculation area for the tsunami numerical simulation (450 m mesh).
Geosciences 08 00017 g007
Figure 8. Relationships between the annual exceedance probability and tsunami wave height at (a) a water depth of 10 m (37.85127° N, 140.96812° E) off the port of Soma; (b) a water depth of 10 m (37.85127° N, 140.96812° E) off the coast of Sendai Plain; and (c) a water depth of 10 m (37.85127° N, 140.96812° E) off the coast of the city of Kesennuma.
Figure 8. Relationships between the annual exceedance probability and tsunami wave height at (a) a water depth of 10 m (37.85127° N, 140.96812° E) off the port of Soma; (b) a water depth of 10 m (37.85127° N, 140.96812° E) off the coast of Sendai Plain; and (c) a water depth of 10 m (37.85127° N, 140.96812° E) off the coast of the city of Kesennuma.
Geosciences 08 00017 g008
Figure 9. Five regions of terrain data used for the tsunami numerical run-up simulations in the Soma area (Domain 1, Domain 2, Domain 3, Domain 4 and Domain 5). The numbers in parentheses are the east-west and north-south mesh dimensions.
Figure 9. Five regions of terrain data used for the tsunami numerical run-up simulations in the Soma area (Domain 1, Domain 2, Domain 3, Domain 4 and Domain 5). The numbers in parentheses are the east-west and north-south mesh dimensions.
Geosciences 08 00017 g009
Figure 10. Five regions of terrain data used for the tsunami numerical run-up simulations in the Sendai and Kesennuma areas (Domain 1, Domain 2, Domain 3, Domain 4 and Domain 5). The numbers in parentheses are the east-west and north-south mesh dimensions.
Figure 10. Five regions of terrain data used for the tsunami numerical run-up simulations in the Sendai and Kesennuma areas (Domain 1, Domain 2, Domain 3, Domain 4 and Domain 5). The numbers in parentheses are the east-west and north-south mesh dimensions.
Geosciences 08 00017 g010
Figure 11. Results of tsunami inundation assessments for (a) Soma; (b) Sendai and (c) Kesennuma at return periods of approximately 200, 700 and 1500 years.
Figure 11. Results of tsunami inundation assessments for (a) Soma; (b) Sendai and (c) Kesennuma at return periods of approximately 200, 700 and 1500 years.
Geosciences 08 00017 g011
Figure 12. Relationships between the annual exceedance probability and tsunami inundation height for (a) Soma; (b) Sendai and (c) Kesennuma.
Figure 12. Relationships between the annual exceedance probability and tsunami inundation height for (a) Soma; (b) Sendai and (c) Kesennuma.
Geosciences 08 00017 g012aGeosciences 08 00017 g012b
Figure 13. Tsunami fragility curves for different building structures ((a) Reinforced concrete; (b) Steel; (c) Brick and (d) Wood) and different damage levels, which were created by regressing the damage data from the 3.11 Tohoku earthquake [20].
Figure 13. Tsunami fragility curves for different building structures ((a) Reinforced concrete; (b) Steel; (c) Brick and (d) Wood) and different damage levels, which were created by regressing the damage data from the 3.11 Tohoku earthquake [20].
Geosciences 08 00017 g013
Figure 14. Tsunami risk curves for different types of building construction with different percentile of tsunami hazard for three target points.
Figure 14. Tsunami risk curves for different types of building construction with different percentile of tsunami hazard for three target points.
Geosciences 08 00017 g014
Figure 15. Expected value of the tsunami risk using each tsunami hazard value ((a) average value; (b) 5th percentile value; (c) 50th percentile value and (d) 95th percentile value) indicated according to the structure (upper horizontal axis) and the region (lower horizontal axis).
Figure 15. Expected value of the tsunami risk using each tsunami hazard value ((a) average value; (b) 5th percentile value; (c) 50th percentile value and (d) 95th percentile value) indicated according to the structure (upper horizontal axis) and the region (lower horizontal axis).
Geosciences 08 00017 g015
Figure 16. Expected value of the tsunami risk at each risk assessment point ((a) Soma; (b) Sendai; and (c) Kesennuma) according to the structure (upper horizontal axis) and the percentile value of the tsunami hazard (lower horizontal axis).
Figure 16. Expected value of the tsunami risk at each risk assessment point ((a) Soma; (b) Sendai; and (c) Kesennuma) according to the structure (upper horizontal axis) and the percentile value of the tsunami hazard (lower horizontal axis).
Geosciences 08 00017 g016
Table 1. Location information for the target buildings in each city.
Table 1. Location information for the target buildings in each city.
Target LocationLatitude (°)Longitude (°)Elevation (m)Distance from Coastline (m)
Soma37.84829140.950523.7390
Sendai38.23349140.984791.11040
Kesennuma38.87447141.587751.490
Table 2. Abbreviation explanations.
Table 2. Abbreviation explanations.
AbbreviationEarthquake Name
JTN1-1Large interplate earthquakes in Northern Sanriku-Oki (repeating earthquakes)
JTN1-2Large interplate earthquakes in Northern Sanriku-Oki (other than repeating earthquakes)
JTN2Miyagi-ken-Oki earthquake (repeating earthquakes)
JTN3-1Earthquakes close to the offshore trenches in Southern Sanriku-Oki (repeating earthquakes)
JTN3-2Earthquakes close to the offshore trenches in Southern Sanriku-Oki (other than repeating earthquakes)
JTN2 + JTN3Miyagi-ken-Oki, earthquakes close to the offshore trenches in Southern Sanriku-Oki consolidated-type-earthquake
TOHOKUGreat East Japan Earthquake (2011 Tohoku-type earthquake)
JTTLarge interplate earthquakes close to the offshore trenches in the Sanriku-Oki to Boso-Oki regions (tsunami earthquakes)
JTNRLarge intraplate earthquakes close to the offshore trenches in the Sanriku-Oki to Boso-Oki regions (normal fault-type)
JTS1Interplate earthquakes in Fukushima-ken-Oki
IBRKInterplate earthquakes in Ibaraki-ken-Oki (other than repeating earthquakes)
Table 3. Model names for the generation intervals of earthquakes, α values of the BPT distribution, average return periods, sample periods, earthquake generation times within the periods used to determine the average occurrence intervals and lower and upper limits of the confidence intervals for the eleven earthquakes, which are shown in the Headquarters for Earthquake Research Promotion [16].
Table 3. Model names for the generation intervals of earthquakes, α values of the BPT distribution, average return periods, sample periods, earthquake generation times within the periods used to determine the average occurrence intervals and lower and upper limits of the confidence intervals for the eleven earthquakes, which are shown in the Headquarters for Earthquake Research Promotion [16].
Earthquake Name (Abbreviation)Model for Generation Interval of Earthquakeα Value of BPT DistributionAverage Return Period (Year)Sample Period (Year)Earthquake Generation TimeLower Limit of Confidence Interval for the Return PeriodUpper Limit of Confidence Interval for the Return Period
JTN1-1BPT0.0897412493101
0.1897412489106
0.2897412484112
JTN2Poisson process-3811041553
JTN3-1BPT0.121092203102117
0.22109220396124
0.32109220391131
JTN2 + JTN3Poisson process-2182181661260
JTTPoisson process-103400456191
JTNRPoisson process-57557511743324
JTN1-2Poisson process-1412791021
JTN3-2Poisson process-4212732193
JTS1Poisson process-206412289582
IBRKPoisson process-2612751545
TOHOKUBPT0.1460024004559644
0.2460024004532676
0.3460024004506711
Table 4. Fault parameters for the reference magnitude of each assumed fault.
Table 4. Fault parameters for the reference magnitude of each assumed fault.
Earthquake Name (Abbreviation)Moment Magnitude (Mw)Number of Earthquake FaultsEarthquake Fault ParameterEarthquake Moment Mo (Nm)Shear Modulus μ (N/m2)Average Slip (m)
Longitude (°)Lattitude (°)Depth (km)Length (km)Width (km)Strike (°)Rake (°)Dip (°)
JTN1-18.2 1 143.09641.6039.017010015690202.51 × 10213.50 × 10103.30
JTN27.6 1 142.38838.45430.8606019490222.72 × 10207.00 × 10101.10
JTN3-17.9 1 143.20338.7416.85016519190109.43 × 10205.00 × 10103.30
JTN2+JTN38.1 1 The fault parameters were assumed as JTN2, JTN3-1 Consolidated Type Earthquake1.78 × 10215.00 × 10102.20
TOHOKU9.0 1 144.13939.94.150020019390133.98 × 1022Shallow part: 3.60 × 1010, Deep part: 5.20 × 1010Shallow part: 18.1, Deep part: 6.0
JTT8.01144.72941.0876.5200501929061.26 × 10213.50 × 10103.60
2144.42841.1878.1200501929091.26 × 10213.50 × 10103.60
3144.27839.3516.8200501909071.26 × 10213.50 × 10103.60
4143.92839.40110.3200501909081.26 × 10213.50 × 10103.60
5143.91537.467.0200502119071.26 × 10213.50 × 10103.60
6143.637.57210.3200502109091.26 × 10213.50 × 10103.60
7142.4436.0587.0200501899081.26 × 10213.50 × 10103.60
8142.2136.03110.82005019390101.26 × 10213.50 × 10103.60
JTNR8.31144.70641.0880.0200100192270453.55 × 10217.00 × 10102.50
2144.22639.3530.0200100190270453.55 × 10217.00 × 10102.50
3143.85637.4580.0200100211270453.55 × 10217.00 × 10102.50
4142.40736.0580.0200100189270453.55 × 10217.00 × 10102.50
5145.48640.9880.0200100192270453.55 × 10217.00 × 10102.50
6145.00639.2530.0200100190270453.55 × 10217.00 × 10102.50
7144.73637.3580.0200100211270453.55 × 10217.00 × 10102.50
8143.18635.9580.0200100189270453.55 × 10217.00 × 10102.50
JTN1-27.51143.71741.4317.5606020690132.24 × 10205.00 × 10101.50
2142.96641.72634.4606020790172.24 × 10205.00 × 10101.50
3143.40541.09919.3606018590142.24 × 10205.00 × 10101.50
4142.65841.30435.4606018490192.24 × 10205.00 × 10101.50
5143.31340.65420.1606018490142.24 × 10205.00 × 10101.50
6142.60540.75435.4606018490192.24 × 10205.00 × 10101.50
7143.22240.19421.3606018690142.24 × 10205.00 × 10101.50
8142.55640.234.9606018590202.24 × 10205.00 × 10101.50
JTN3-27.41143.26738.72716505019290101.58 × 10205.00 × 10101.30
2143.12138.33216.5505019190101.58 × 10205.00 × 10101.30
3142.99737.92216.9505019390101.58 × 10205.00 × 10101.30
JTS17.41142.80437.29817.3505020490101.58 × 10205.00 × 10101.30
2142.37237.44523.4505020190151.58 × 10205.00 × 10101.30
3141.94837.56434.6505020490231.58 × 10205.00 × 10101.30
4142.55636.87715.7505021290131.58 × 10205.00 × 10101.30
5142.1237.03624.6505020590161.58 × 10205.00 × 10101.30
6141.69737.14935.3505020790221.58 × 10205.00 × 10101.30
IBRK7.51142.236.51516555520790152.24 × 10205.00 × 10101.40
2141.46736.7734.7555520590212.24 × 10205.00 × 10101.40
3141.88136.14917.5555520190152.24 × 10205.00 × 10101.40
4141.22236.36736.1555519490192.24 × 10205.00 × 10101.40
5141.71435.9918.4555519790172.24 × 10205.00 × 10101.40
6141.29936.12731.6555519490202.24 × 10205.00 × 10101.40
Table 5. Calculation conditions for the tsunami numerical simulation.
Table 5. Calculation conditions for the tsunami numerical simulation.
ItemCalculation Condition
Governing equation2D non-linear shallow water equation (Tohoku University TUNAMI model) [19]
Numerical integration methodStaggered leap-frog differential method
Initial conditionOkada equation [18]
Boundary conditionOpen boundary
Coordination systemSpherical coordinate system
Tidal settingT.P. +0.0 m
Mesh size450 m
Time step0.9 s
Calculation time3 h
Table 6. Calculation results for Soma. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
Table 6. Calculation results for Soma. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
(a)(b)(c)(d)(e)(f)(g)(h)
Earthquake SourceMoment Magnitude
(Mw)
Position of AsperityTsunami Height (m)
(10 m Water Depth Point)
Annual Exceedance ProbabilityReturn Period
(year)
Tsunami Inundation Height (m)
(Risk Assessment Point)
Tsunami Inundation Depth (m)
(Risk Assessment Point)
TOHOKU9.1Center8.580.000294340513.639.96
9.1Between south and center7.840.000423236310.336.66
9.1South7.840.000425235512.548.87
9.0Center7.100.000579172811.547.87
9.1Between north and center6.780.000655152610.566.89
9.0South6.580.000702142511.137.46
9.1North6.510.000720138810.687.01
9.0Between south and center6.480.000725137910.416.74
9.0Between north and center5.490.00099210089.385.71
8.9Center5.440.0010099919.275.60
9.0North5.420.0010159869.355.68
8.9South5.280.0010529509.655.98
8.9Between south and center4.930.0011718548.544.87
8.9Between north and center4.400.0013817248.084.41
8.9North4.080.0015386507.894.22
JTNR8.4North2.500.0043162324.300.63
8.4Center2.470.0044642243.970.30
8.4South2.370.0050141993.840.17
JTN2 + JTN38.2South2.290.0054591834.510.84
8.2Center2.260.0056671764.360.69
JTT8.1Center2.250.0057381744.150.48
8.1North2.230.0058851704.100.43
8.1Center2.230.0059611683.930.26
JTNR8.4South2.230.0059611684.280.61
8.4North2.200.0061951614.881.21
8.4Center2.190.0062751594.350.68
JTT8.1North2.140.0066121513.780.11
8.0Center2.010.0079861250.000.00
8.1South1.960.0084411180.000.00
8.0Center1.960.0084411180.000.00
Table 7. Calculation results for Sendai. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
Table 7. Calculation results for Sendai. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
(a)(b)(c)(d)(e)(f)(g)(h)
Earthquake SourceMoment Magnitude
(Mw)
Position of AsperityTsunami Height (m)
(10 m Water Depth Point)
Annual Exceedance ProbabilityReturn Period
(year)
Tsunami Inundation Height (m)
(Risk Assessment Point)
Tsunami Inundation Depth (m)
(Risk Assessment Point)
TOHOKU9.1Center9.890.00022045366.845.74
9.0Center8.500.00045921775.974.87
9.1Between north and center7.950.00060316575.063.96
9.1South7.580.00071214044.663.56
8.9Center6.820.0010239784.273.17
9.0Between north and center6.780.0010419604.263.16
9.1Between south and center6.640.0011238903.842.74
9.0South6.520.0011978363.772.67
9.1North6.380.0013397473.852.75
8.9Between north and center5.620.0023924182.741.64
9.0Between south and center5.560.0024934012.151.05
9.0North5.520.0025923862.951.85
8.9South5.300.0030473282.251.15
8.9North4.490.0057331741.760.66
8.9Between south and center4.300.0067291490.000.00
JTN2 + JTN38.2South4.160.0074651340.000.00
8.2Center4.040.0082851210.000.00
8.2North3.790.010163980.000.00
JTN3-18.0South3.760.010331970.000.00
8.0Center3.610.011687860.000.00
JTN2 + JTN38.1Center3.410.013891720.000.00
JTN3-18.0North3.360.014356700.000.00
JTN2 + JTN38.1South3.310.014959670.000.00
JTT8.1South3.160.017056590.000.00
JTN2 + JTN38.1North3.130.017479570.000.00
JTT8.1Center3.010.019274520.000.00
JTN3-17.9Center2.930.020564490.000.00
JTNR8.4South2.900.021068470.000.00
8.4North2.870.021410470.000.00
8.4Center2.850.021757460.000.00
Table 8. Calculation results for Kesennuma. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
Table 8. Calculation results for Kesennuma. (a) Earthquake source name; (b) moment magnitude of the earthquake; (c) position of the asperity; (d) tsunami wave height at a water depth of 10 m; (e) annual exceedance probability estimated using the tsunami hazard curve; (f) return period calculated from the annual exceedance probability; (g) tsunami inundation height at the risk assessment point simulated using a nonlinear longwave equation with the fault parameters.
(a)(b)(c)(d)(e)(f)(g)(h)
Earthquake SourceMoment Magnitude
(Mw)
Position of AsperityTsunami Height (m)
(10 m Water Depth Point)
Annual Exceedance ProbabilityReturn Period
(year)
Tsunami Inundation Height (m)
(Risk Assessment Point)
Tsunami Inundation Depth (m)
(Risk Assessment Point)
TOHOKU9.1 Between north and center19.44 0.00013275742.91 1.51
9.1 Center17.01 0.0002546342.97 1.57
9.0 Between north and center15.66 0.0003345313.18 1.78
9.0 Center14.21 0.0004444403.16 1.76
9.1 Between south and center12.41 0.0006164162.73 1.33
9.1 North12.24 0.0006274056.08 4.68
8.9 Between north and center10.50 0.0008533273.26 1.86
8.9 Center10.48 0.0008553263.08 1.68
9.0 North10.11 0.0008963185.13 3.73
9.0 Between south and center9.30 0.0009923092.78 1.38
9.1 South9.13 0.0010193094.35 2.95
8.9 North7.82 0.0011913094.34 2.94
9.0 South6.85 0.0013203073.71 2.31
8.9 Between south and center6.44 0.0013833052.46 1.06
8.9 South4.92 0.0019511943.14 1.74
JTT8.1 South3.63 0.004062780.00 0.00
8.1 South3.61 0.004093770.00 0.00
8.1 Center3.56 0.004288730.00 0.00
8.1 North3.46 0.004608680.00 0.00
8.1 Center3.37 0.005007610.00 0.00
JTNR8.4 South3.27 0.005458550.00 0.00
8.4 North3.26 0.005506550.00 0.00
JTT8.1 North3.19 0.005808510.00 0.00
JTNR8.4 Center3.18 0.005860500.00 0.00
JTN2 + JTN38.2 North3.05 0.006667440.00 0.00
JTNR8.4 South2.86 0.008003350.00 0.00
8.4 Center2.82 0.008243340.00 0.00
JTT8.0 South2.82 0.008325340.00 0.00
8.0 Center2.78 0.008662320.00 0.00
8.0 South2.78 0.008662320.00 0.00
Table 9. Calculation conditions for the nonlinear longwave equation used in the tsunami numerical simulation for Soma.
Table 9. Calculation conditions for the nonlinear longwave equation used in the tsunami numerical simulation for Soma.
ItemCalculation Condition
Governing equation2D non-linear shallow water equation (Tohoku University TUNAMI model) [19]
Numerical integration methodStaggered leap-frog differential method
Initial conditionOkada equation [18]
Boundary conditionRun-up boundary
Coordinate systemPlane rectangular coordinate system IX
Tidal settingT.P. +0.0 m
Mesh size810 m, 270 m, 90 m, 30 m, 10 m
Time step0.9 s, 0.3 s, 0.1 s, 0.03 s, 0.01 s
Calculation time3 h
Table 10. Calculation conditions for the nonlinear longwave equation used in the tsunami numerical simulations for Sendai and Kesennuma.
Table 10. Calculation conditions for the nonlinear longwave equation used in the tsunami numerical simulations for Sendai and Kesennuma.
ItemCalculation Condition
Governing equation2D non-linear shallow water equation (Tohoku University TUNAMI model) [19]
Numerical integration methodStaggered leap-frog differential method
Initial conditionOkada equation [18]
Boundary conditionRun-up boundary
Coordinate systemPlane rectangular coordinate system X
Tidal settingT.P. +0.0 m
Mesh size1215 m, 405 m, 135 m, 45 m, 15 m
Time step0.9 s, 0.3 s, 0.1 s, 0.03 s, 0.01 s
Calculation time3 h

Share and Cite

MDPI and ACS Style

Fukutani, Y.; Suppasri, A.; Imamura, F. Quantitative Assessment of Epistemic Uncertainties in Tsunami Hazard Effects on Building Risk Assessments. Geosciences 2018, 8, 17. https://doi.org/10.3390/geosciences8010017

AMA Style

Fukutani Y, Suppasri A, Imamura F. Quantitative Assessment of Epistemic Uncertainties in Tsunami Hazard Effects on Building Risk Assessments. Geosciences. 2018; 8(1):17. https://doi.org/10.3390/geosciences8010017

Chicago/Turabian Style

Fukutani, Yo, Anawat Suppasri, and Fumihiko Imamura. 2018. "Quantitative Assessment of Epistemic Uncertainties in Tsunami Hazard Effects on Building Risk Assessments" Geosciences 8, no. 1: 17. https://doi.org/10.3390/geosciences8010017

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop