Engineering Applications Using Probabilistic Aftershock Hazard Analyses : Aftershock Hazard Map and Load Combination of Aftershocks and Tsunamis

After the Tohoku earthquake in 2011, we observed that aftershocks tended to occur in a wide region after such a large earthquake. These aftershocks resulted in secondary damage or delayed rescue and recovery activities. In addition, it has been reported that there are regions where the intensity of the vibrations owing to the aftershocks was much stronger than those associated with the main shock. Therefore, it is necessary to consider the seismic risk associated with aftershocks. We used the data regarding aftershocks that was obtained from the Tohoku earthquake and various other historically large earthquakes. We investigated the spatial and temporal distribution of the aftershocks using the Gutenberg–Richter law and the modified Omori law. Subsequently, we previously proposed a probabilistic aftershock occurrence model that is expected to be useful to develop plans for recovery activities after future large earthquakes. In this study, the probabilistic aftershock hazard analysis is used to create aftershock hazard maps. We propose a hazard map focusing on the probability of aftershocks on the scale of the main shock for use with a recovery activity plan. Following the lessons learned from the 2011 Tohoku earthquake, we focus on the simultaneous occurrence of tsunamis and aftershocks just after a great subduction earthquake. The probabilistic aftershock hazard analysis is used to derive load combination equations of the load and resistance factor design. This design is intended to simultaneously consider tsunamis and aftershocks for tsunami-resistant designs of tsunami evacuation buildings.


Introduction
An earthquake with a moment magnitude of 9.0 struck Tohoku by the coast near the Pacific Ocean on 11 March 2011.Subsequently, thousands of aftershocks occurred and the distribution of their hypocenters ranged from the Iwate Sea to the Ibaraki Sea.Therefore, it is clear that it is necessary for emergency workers to consider the occurrence of aftershocks so as to be prepared to face such a scenario in the future.
The importance of including recovery activities in contingency plans has increasingly been recognized.The reports of damages due to the aftershocks caused by the 2011 Tohoku earthquake included collapses of buildings that were damaged by the main shock, delays in recovery activities owing to water cutoffs or no-passing zones, and delays in rescue activities.However, the Central Disaster Management Council of Japan [1] does not consider the effect of aftershocks.Probabilistic methods to evaluate the hazards of aftershocks have been deemed to be particularly effective in improving decision-making and planning due to various uncertainties that are associated with earthquakes in terms of the location, scale, and frequency of aftershocks [2,3].
Multiple studies of aftershock hazard analyses were conducted after the Niigata Chuetsu earthquake in 2004 [4][5][6][7].Standard analysis methods were proposed by the Headquarters for Earthquake Research Promotion to design probabilistic aftershock occurrence models [8].However, aftershocks with a magnitude greater than 7.0 after the occurrence of a main shock with a magnitude of 9.0 were not reviewed.Previously, aftershock risk analyses were conducted immediately after the occurrence of a main shock [9].However, it is difficult to evaluate the hazards of an aftershock before the main shock due to various uncertainties.For possible great earthquakes, we must make decisions based on such uncertainties, and it is important to quantify the various uncertainties.Therefore, a probabilistic model that takes into account the various uncertainties is necessary to review the safety measures before the occurrence of an earthquake.We previously proposed a probabilistic aftershock hazard analysis method for trench-type mega-earthquakes [10].
In this paper, engineering applications of the proposed approach for probabilistic aftershock hazard analysis are shown for demonstration purposes.One application is to use aftershock hazard maps to plan recovery activities.Another application is to derive load combination equations of the load and resistance factor design (LRFD) considering the simultaneous occurrence of tsunamis and aftershocks for the tsunami-resistant design of tsunami evacuation buildings.

Aftershock Definition
The classification of aftershocks [11] used in our study is illustrated in Figure 1.In general, many earthquakes subsequently occur near the hypocenter of a large earthquake after its occurrence.If the subsequent earthquakes are smaller in magnitude than the initial earthquake, the initial one is classified as the main shock and the subsequent ones are classified as aftershocks.Conversely, if a subsequent earthquake is larger in magnitude than the initial earthquake, the initial earthquake is classified as a foreshock and the largest subsequent earthquake is classified as the main shock.Typically, subduction earthquakes stronger than a Richter magnitude of 7.0 are assumed to be main shocks because the majority of such earthquakes are historically observed to be main shocks.
The area where the aftershocks occur is the aftershock region, which contains the fault plane as well as the expanded fault plane.The distribution of aftershocks resulting from the Tohoku earthquake as per the Japan Meteorological Agency [12] is depicted in Figure 2. Aftershocks that occur at a large distance from the fault plane, i.e., induced earthquakes (e.g., on 12 March 2011, or 15 March 2011), are not included in this study.
Geosciences 2018, 8, 1 2 of 22 methods to evaluate the hazards of aftershocks have been deemed to be particularly effective in improving decision-making and planning due to various uncertainties that are associated with earthquakes in terms of the location, scale, and frequency of aftershocks [2,3].Multiple studies of aftershock hazard analyses were conducted after the Niigata Chuetsu earthquake in 2004 [4][5][6][7].Standard analysis methods were proposed by the Headquarters for Earthquake Research Promotion to design probabilistic aftershock occurrence models [8].However, aftershocks with a magnitude greater than 7.0 after the occurrence of a main shock with a magnitude of 9.0 were not reviewed.Previously, aftershock risk analyses were conducted immediately after the occurrence of a main shock [9].However, it is difficult to evaluate the hazards of an aftershock before the main shock due to various uncertainties.For possible great earthquakes, we must make decisions based on such uncertainties, and it is important to quantify the various uncertainties.Therefore, a probabilistic model that takes into account the various uncertainties is necessary to review the safety measures before the occurrence of an earthquake.We previously proposed a probabilistic aftershock hazard analysis method for trench-type mega-earthquakes [10].
In this paper, engineering applications of the proposed approach for probabilistic aftershock hazard analysis are shown for demonstration purposes.One application is to use aftershock hazard maps to plan recovery activities.Another application is to derive load combination equations of the load and resistance factor design (LRFD) considering the simultaneous occurrence of tsunamis and aftershocks for the tsunami-resistant design of tsunami evacuation buildings.

Aftershock Definition
The classification of aftershocks [11] used in our study is illustrated in Figure 1.In general, many earthquakes subsequently occur near the hypocenter of a large earthquake after its occurrence.If the subsequent earthquakes are smaller in magnitude than the initial earthquake, the initial one is classified as the main shock and the subsequent ones are classified as aftershocks.Conversely, if a subsequent earthquake is larger in magnitude than the initial earthquake, the initial earthquake is classified as a foreshock and the largest subsequent earthquake is classified as the main shock.Typically, subduction earthquakes stronger than a Richter magnitude of 7.0 are assumed to be main shocks because the majority of such earthquakes are historically observed to be main shocks.
The area where the aftershocks occur is the aftershock region, which contains the fault plane as well as the expanded fault plane.The distribution of aftershocks resulting from the Tohoku earthquake as per the Japan Meteorological Agency [12] is depicted in Figure 2. Aftershocks that occur at a large distance from the fault plane, i.e., induced earthquakes (e.g., on 12 March 2011, or 15 March 2011), are not included in this study.In this study, we assume that the relationship between the magnitude and frequency of aftershocks can be expressed using the Gutenberg-Richter (GR) law [13].Further, we observe that the relationship between the time elapsed after the main shock and the occurrence of the aftershocks agree with the modified Omori's (MO) law [14][15][16].Consequently, Reasenberg and Jones proposed a model that can describe the rate of occurrence of aftershocks with a magnitude greater than or equal to M at a time t immediately after a main shock of magnitude Mm, and this can be expressed as follows [17]: ( , , ) = 10 where a is proportional to the occurrence of earthquakes, b describes the magnitude distribution of the earthquakes, p is the decay rate of the aftershocks, and c is the time-offset parameter, which is assumed to be a constant of 0.1 in this study.The total number of aftershocks is calculated by integrating Equation ( 1): where T1 is the elapsed time from the main shock.In Equation (2), the relationship between a and the magnitude of the main shock Mm is not specified.However, it is constant regardless of the elapsed time from the main shock.Therefore, the change in the value of b owing to the elapsed time is not considered.

Reasenberg and Jones Model
In this study, we assume that the relationship between the magnitude and frequency of aftershocks can be expressed using the Gutenberg-Richter (GR) law [13].Further, we observe that the relationship between the time elapsed after the main shock and the occurrence of the aftershocks agree with the modified Omori's (MO) law [14][15][16].Consequently, Reasenberg and Jones proposed a model that can describe the rate of occurrence of aftershocks with a magnitude greater than or equal to M at a time t immediately after a main shock of magnitude M m , and this can be expressed as follows [17]: where a is proportional to the occurrence of earthquakes, b describes the magnitude distribution of the earthquakes, p is the decay rate of the aftershocks, and c is the time-offset parameter, which is assumed to be a constant of 0.1 in this study.The total number of aftershocks is calculated by integrating Equation (1): where T 1 is the elapsed time from the main shock.In Equation (2), the relationship between a and the magnitude of the main shock M m is not specified.However, it is constant regardless of the elapsed time from the main shock.Therefore, the change in the value of b owing to the elapsed time is not considered.

Proposed Probabilistic Aftershock Occurrence Model
We previously proposed Equation (3) [10], which contains a parameter a as a function of the magnitude of the main shock and a parameter b as a function of the elapsed time from the main shock T 1 .In addition, the minimum magnitude of the aftershocks (M min ) is set to 4.0 [10].
Further, when the probability distribution of the earthquake scale is designed using the GR law, the probability that the aftershock magnitude M becomes m i at T 1 is where ∆m is assumed to be 0.1.In Equation (3), it is assumed that the maximum aftershock magnitude (M max ) is equal to the main shock magnitude (M m ).However, we propose the following Equation ( 5) according to previous studies [18,19].Hereafter, M m is represented in terms of M max using Equation ( 5): (5)

Model Parameters
Historically large earthquakes including the Tohoku earthquake of 2011, the Tokachi-Oki earthquake of 2003, and other such major earthquakes were compiled to construct an aftershock occurrence model for large earthquakes that may occur in the future.We define aftershocks as earthquakes that occur within 90 days [20] in the area ruptured by the main shock with a focal depth shallower than 100 km and a magnitude (M J ) larger than or equal to 4.0.The hypocenter data (time, location, and magnitude) obtained from the Japan Meteorological Agency were used [21,22].The spatial distribution of the aftershocks is illustrated in Figure 3.We previously proposed Equation (3) [10], which contains a parameter a as a function of the magnitude of the main shock and a parameter b as a function of the elapsed time from the main shock T1.In addition, the minimum magnitude of the aftershocks (Mmin) is set to 4.0 [10].
Further, when the probability distribution of the earthquake scale is designed using the GR law, the probability that the aftershock magnitude M becomes mi at T1 is where Δm is assumed to be 0.1.In Equation (3), it is assumed that the maximum aftershock magnitude (Mmax) is equal to the main shock magnitude (Mm).However, we propose the following Equation ( 5) according to previous studies [18,19].Hereafter, Mm is represented in terms of Mmax using Equation ( 5):

Model Parameters
Historically large earthquakes including the Tohoku earthquake of 2011, the Tokachi-Oki earthquake of 2003, and other such major earthquakes were compiled to construct an aftershock occurrence model for large earthquakes that may occur in the future.We define aftershocks as earthquakes that occur within 90 days [20] in the area ruptured by the main shock with a focal depth shallower than 100 km and a magnitude (MJ) larger than or equal to 4.0.The hypocenter data (time, location, and magnitude) obtained from the Japan Meteorological Agency were used [21,22].The spatial distribution of the aftershocks is illustrated in Figure 3.

Target Aftershocks
For the 2011 Tohoku earthquake, we used the fault model of the Geospatial Information Authority of Japan [23] to estimate the expanded fault plane, which acts as the aftershock area.For the 2003 Tokachi-Oki earthquake, the Yagi fault model [24] was used to estimate the expanded fault

Target Aftershocks
For the 2011 Tohoku earthquake, we used the fault model of the Geospatial Information Authority of Japan [23] to estimate the expanded fault plane, which acts as the aftershock area.For the 2003 Tokachi-Oki earthquake, the Yagi fault model [24] was used to estimate the expanded fault plane.The aftershock areas of other major earthquakes are assumed to be circular regions centered on the hypocenter of the main shock.The aftershock area A (km 2 ) is calculated as follows [20]:

Relationships between Parameters
The parameters of the GR and MO models obtained from the previous aftershock data are depicted in Table 1, which includes the parameters used in previous studies [25,26].Figure 4 illustrates the relationship between M m and each declared parameter.The number of aftershocks (N) increases with increasing M m .However, we cannot confirm the relationships between the other parameters.We use regression to propose the following equations (average and standard deviation) as the aftershock occurrence model parameters.
log 10 A a f ter = 0.778M m − 1.60 ( 7) p = 1.05, σ p = 0.17    In addition, Figure 5a depicts the GR model at various elapsed times (T 1 ) of the Tohoku earthquake.The b-value (the slope of the GR equation) increases with increasing elapsed time.Figure 5b illustrates the relationship between elapsed time (T 1 ) and b(T 1 ) for major earthquakes in the past.As depicted in Figure 5b, b(T 1 ) is given by the following equation: where N is the number of aftershocks that occurred in the aftershock region within a period of 90 days having a magnitude greater than 4.0 (the results of Guo and Ogata [25] and Utsu [26] were excluded due to their having different conditions than the ones required), b(90) is the b-value of the GR equation 90 days after the occurrence of the main shock, p is the decay rate of the MO law, and D 1 is the magnitude difference between the main shock and the largest aftershock.

Aftershock Hazard Analysis Equation
Seismic hazard analyses are used to assess the relationship between the expected earthquake ground motion at a certain period of time and the probability or frequency of exceedance.In this study, we concentrate on aftershocks due to large subduction earthquakes.The average number of aftershock occurrences for days with a ground motion intensity (Ysurf) greater than y is [24] where ( , , ) is the average number of aftershocks having magnitudes greater than 4.0 within the aftershock region, > , is the conditional probability of a ground motion intensity ( ) that is greater than y when the aftershock magnitude is and the distance is , ( | ) is the probability function of the aftershock magnitude in the aftershock region at an elapsed time , and is the probability function of the distance at an aftershock magnitude .We assume a uniform probability distribution in the aftershock region.In addition, we use the following attenuation relation [27] to calculate the peak ground velocity (PGV) as the aftershock hazard, which can more stably estimate the structural damage than the peak ground acceleration (PGA): where is the PGV (cm/s) on the engineering bedrock, M is the aftershock magnitude, D is the focal depth (km), X is the shortest source distance (km), and ε is the deviation between the observed and predicted values with a logarithm standard deviation of 0.23.X is the distance from a spherical surface of radius 1/2L around the hypocenter and is calculated using Equation (15) in accordance with the study by Utsu [28].In addition, we set X to 3 km when the distance is shorter than 1/2L.
Site effects are calculated using the following equations:

Aftershock Hazard Analysis Equation
Seismic hazard analyses are used to assess the relationship between the expected earthquake ground motion at a certain period of time and the probability or frequency of exceedance.In this study, we concentrate on aftershocks due to large subduction earthquakes.The average number of aftershock occurrences for T 1 days with a ground motion intensity (Y surf ) greater than y is [24] where N(T 1 , M min , M m ) is the average number of aftershocks having magnitudes greater than 4.0 within the aftershock region, P Y sur f > y m i , r j is the conditional probability of a ground motion intensity (Y sur f ) that is greater than y when the aftershock magnitude is m i and the distance is r j , P(m i |T 1 ) is the probability function of the aftershock magnitude in the aftershock region at an elapsed time T 1 , and P r j m i is the probability function of the distance at an aftershock magnitude m i .We assume a uniform probability distribution in the aftershock region.In addition, we use the following attenuation relation [27] to calculate the peak ground velocity (PGV) as the aftershock hazard, which can more stably estimate the structural damage than the peak ground acceleration (PGA): where Y bed is the PGV (cm/s) on the engineering bedrock, M is the aftershock magnitude, D is the focal depth (km), X is the shortest source distance (km), and ε is the deviation between the observed and predicted values with a logarithm standard deviation of 0.23.X is the distance from a spherical surface of radius 1/2L around the hypocenter and is calculated using Equation ( 15) in accordance with the study by Utsu [28].In addition, we set X to 3 km when the distance is shorter than 1/2L.log 10 L = 0.5M − 1.85 (15) Site effects are calculated using the following equations: where Y sur f is the PGV (cm/s) on the surface, R is the site amplification, and AVS30 (m/s) is the average S-wave velocity down to 30 m.The aftershock region is assumed to be the expanded fault plane, which is an extension of the fault plane in accordance with Equation (7).We divide the expanded fault plane into a 10 km × 10 km mesh, as shown in Figure 6.The aftershock hazard at each site is calculated by integrating the results of each mesh.where is the PGV (cm/s) on the surface, R is the site amplification, and AVS30 (m/s) is the average S-wave velocity down to 30 m.
The aftershock region is assumed to be the expanded fault plane, which is an extension of the fault plane in accordance with Equation (7).We divide the expanded fault plane into a 10 km × 10 km mesh, as shown in Figure 6.The aftershock hazard at each site is calculated by integrating the results of each mesh.

Verification of the Proposed Method
We calculated the PGV hazard curve for the aftershocks of the Tohoku earthquake to validate the proposed model using the proposed aftershock hazard analysis method.The site location and aftershock region of the model are exhibited in Figure 7 and the site-specific parameters are listed in Table 2.

Verification of the Proposed Method
We calculated the PGV hazard curve for the aftershocks of the Tohoku earthquake to validate the proposed model using the proposed aftershock hazard analysis method.The site location and aftershock region of the model are exhibited in Figure 7 and the site-specific parameters are listed in Table 2.
where is the PGV (cm/s) on the surface, R is the site amplification, and AVS30 (m/s) is the average S-wave velocity down to 30 m.
The aftershock region is assumed to be the expanded fault plane, which is an extension of the fault plane in accordance with Equation (7).We divide the expanded fault plane into a 10 km × 10 km mesh, as shown in Figure 6.The aftershock hazard at each site is calculated by integrating the results of each mesh.

Verification of the Proposed Method
We calculated the PGV hazard curve for the aftershocks of the Tohoku earthquake to validate the proposed model using the proposed aftershock hazard analysis method.The site location and aftershock region of the model are exhibited in Figure 7 and the site-specific parameters are listed in Table 2.We compared the actual observations to the results of the aftershock hazard analysis to validate the proposed model.The aftershock hazard analysis for a period of 90 days after the major shock is calculated by considering the average as well as the average ± one standard deviation of each parameter, and the actual aftershock observation records at each site obtained within 90 days after the main shock are depicted in Figure 8.To compare the aftershock hazard and the main shock PGV, the PGV observations for the main shock are depicted by the dashed line, the points where the aftershock hazard curves and main shock PGVs intersect are checked, and the frequency of exceedance is confirmed.For example, in MYG010, the observations and aftershock hazard analysis are in good agreement.The PGV of the main shock is 56 cm/s.Therefore, the maximum expected aftershock occurrence number within 90 days of the main shock is approximately 0.1 times the average parameter values and approximately 0.5 times the average + one standard deviation.At other sites (MYG004, FKS010, and IBR013), the observations and the aftershock hazard analysis match well at relatively small PGVs but not so well with increasing PGVs.The PGVs of the main shock are 110, 62, and 69 cm/s, whereas the maximum expected aftershocks within 90 days of the main shock are approximately 0.07, 0.4, and 0.4 times the results of the average + one standard deviation, respectively.Therefore, the aftershock hazard analysis matches the observations well using the proposed parameters.However, there are deviations between the expected and observed values.Therefore, we need to evaluate the hazards considering the variations in the parameters.We compared the actual observations to the results of the aftershock hazard analysis to validate the proposed model.The aftershock hazard analysis for a period of 90 days after the major shock is calculated by considering the average as well as the average ± one standard deviation of each parameter, and the actual aftershock observation records at each site obtained within 90 days after the main shock are depicted in Figure 8.To compare the aftershock hazard and the main shock PGV, the PGV observations for the main shock are depicted by the dashed line, the points where the aftershock hazard curves and main shock PGVs intersect are checked, and the frequency of exceedance is confirmed.For example, in MYG010, the observations and aftershock hazard analysis are in good agreement.The PGV of the main shock is 56 cm/s.Therefore, the maximum expected aftershock occurrence number within 90 days of the main shock is approximately 0.1 times the average parameter values and approximately 0.5 times the average + one standard deviation.At other sites (MYG004, FKS010, and IBR013), the observations and the aftershock hazard analysis match well at relatively small PGVs but not so well with increasing PGVs.The PGVs of the main shock are 110, 62, and 69 cm/s, whereas the maximum expected aftershocks within 90 days of the main shock are approximately 0.07, 0.4, and 0.4 times the results of the average + one standard deviation, respectively.Therefore, the aftershock hazard analysis matches the observations well using the proposed parameters.However, there are deviations between the expected and observed values.Therefore, we need to evaluate the hazards considering the variations in the parameters.

Sensitivity Analysis
As previously mentioned, the parameters used to calculate the aftershock hazards (log 10 N(90, 4.0, 9.0), b(90), D 1 , and p) have prediction errors.We performed a sensitivity analysis of these parameters to evaluate the effect of the variations on the accuracy of the analysis.We consider the MYG010 site, using the proposed parameters of Equations ( 7)-( 12), and calculate the aftershock hazards for the average and the average ± one standard deviation of each parameter for 90 days after the earthquake (Figure 9).
The aftershock hazard increases with increasing N-values as well as with decreasing b(90), D 1 , and p-values.The contribution of each parameter to the aftershock hazard varies.

Sensitivity Analysis
As previously mentioned, the parameters used to calculate the aftershock hazards ( (90, 4.0, 9.0), b(90), , and p) have prediction errors.We performed a sensitivity analysis of these parameters to evaluate the effect of the variations on the accuracy of the analysis.We consider the MYG010 site, using the proposed parameters of Equations ( 7)-( 12), and calculate the aftershock hazards for the average and the average ± one standard deviation of each parameter for 90 days after the earthquake (Figure 9).

Short-Term Versus Mid-Term
We compared the short-term (0-3 days) and mid-term (3-30 days) aftershock hazards at MYG010 (Figure 10) to evaluate the differences in the aftershock hazards owing to the elapsed time from the main shock.The mid-term aftershock hazard is calculated by the following equation: We can observe from the mean μ in Figure 10a that the PGVs of the aftershocks that are expected to occur more than once are approximately 18 cm/s in the short-term and approximately 2 cm/s in the mid-term.Figure 10b shows that the maximum expected aftershock occurrence number for μ + σ is 0.6 times in the short-term and nearly 0 times in the mid-term.Therefore, the mid-term hazard is

Short-Term Versus Mid-Term
We compared the short-term (0-3 days) and mid-term (3-30 days) aftershock hazards at MYG010 (Figure 10) to evaluate the differences in the aftershock hazards owing to the elapsed time from the main shock.The mid-term aftershock hazard is calculated by the following equation: We can observe from the mean µ in Figure 10a that the PGVs of the aftershocks that are expected to occur more than once are approximately 18 cm/s in the short-term and approximately 2 cm/s in the mid-term.Figure 10b shows that the maximum expected aftershock occurrence number for µ + σ is 0.6 times in the short-term and nearly 0 times in the mid-term.Therefore, the mid-term hazard is smaller than the short-term hazard; this could be used to estimate the starting time for recovery activities.

Aftershock Hazard Map for Recovery Activities
The strong aftershocks that occur after a large subduction earthquake will adversely affect recovery activities [29][30][31].For example, the gas, water, and electricity supply were disconnected over a wide area owing to an aftershock with a magnitude of 7.2 on 07 April, 2011, nearly a month after the Tohoku earthquake.Further, railway facilities and highways were damaged by the strong aftershocks, affecting the train schedules and resulting in a closure of the highways.In addition, it is crucial to decide whether or not to discontinue operations at nuclear facilities.If large aftershocks occur during rescue and lifeline recovery activities after the main shock, there is a further possibility of secondary damages.As explained above, we need to make decisions as per prescribed priorities and implementation plans regarding rescue activities, lifeline recovery activities, and evacuation plans quickly owing to the subsequent aftershocks following a great earthquake.In such a scenario, a framework that can reasonably predict the aftershock activity after a great earthquake is necessary to make such decisions in a timely manner.

The 1923 Kanto Earthquake
We assume the aftershock hazards for a large earthquake that may occur in the future in the Tokyo metropolitan area [32].Subsequently, we propose a hazard map using the results of the aftershock hazard analysis for the recovery activities.
Initially, we evaluated the aftershock hazards by applying the proposed model to the Kanto earthquake.The aftershock region is assumed to be the expanded fault plane and the main shock magnitude is in the fault plane [33] (Equation ( 6)).The main shock magnitude (Mm) is 7.9, the assumed site is KNG005 on K-NET, and the AVS30 of KNG005 is 443 m/s.The site and expanded fault plane are depicted in Figure 11.
Figure 12 illustrates the short-term aftershock hazards and the main shock PGV.The predicted main shock PGV is 48 cm/s, and the maximum aftershock frequency of three days is approximately 0.1 for μ + σ and approximately 0.02 for μ.

Aftershock Hazard Map for Recovery Activities
The strong aftershocks that occur after a large subduction earthquake will adversely affect recovery activities [29][30][31].For example, the gas, water, and electricity supply were disconnected over a wide area owing to an aftershock with a magnitude of 7.2 on 7 April 2011, nearly a month after the Tohoku earthquake.Further, railway facilities and highways were damaged by the strong aftershocks, affecting the train schedules and resulting in a closure of the highways.In addition, it is crucial to decide whether or not to discontinue operations at nuclear facilities.If large aftershocks occur during rescue and lifeline recovery activities after the main shock, there is a further possibility of secondary damages.As explained above, we need to make decisions as per prescribed priorities and implementation plans regarding rescue activities, lifeline recovery activities, and evacuation plans quickly owing to the subsequent aftershocks following a great earthquake.In such a scenario, a framework that can reasonably predict the aftershock activity after a great earthquake is necessary to make such decisions in a timely manner.

The 1923 Kanto Earthquake
We assume the aftershock hazards for a large earthquake that may occur in the future in the Tokyo metropolitan area [32].Subsequently, we propose a hazard map using the results of the aftershock hazard analysis for the recovery activities.
Initially, we evaluated the aftershock hazards by applying the proposed model to the Kanto earthquake.The aftershock region is assumed to be the expanded fault plane and the main shock magnitude is in the fault plane [33] (Equation ( 6)).The main shock magnitude (M m ) is 7.9, the assumed site is KNG005 on K-NET, and the AVS30 of KNG005 is 443 m/s.The site and expanded fault plane are depicted in Figure 11.

Aftershock Hazard Map
We propose a hazard map to efficiently use the aftershock hazard analysis for recovery activities.We focused on the short-term aftershock hazards and calculated the expected aftershock frequency at 113 sites (Figure 13).Figure 13a illustrates the results for μ + σ, whereas Figure 13b illustrates the results for μ. Figure 13 illustrates that the aftershock frequencies of the main shock PGV are considerably larger at sites that are closer to the southwestern area where the depth of the expanded fault plane is shallow.The maximum value is 0.52 for μ + σ and 0.20 for μ.The aftershock hazards increase near the fault plane.However, the increase in the PGV of the main shock is more significant, whereas the aftershock frequencies of the main shock PGV are negligible.In addition, the aftershock hazards, main shock PGV, and the aftershock frequencies are considerably smaller at sites located on the northeast side with the deep expanded fault plane.Conversely, the aftershock hazards are significant at sites located near the southwest side with the shallow expanded fault plane because the main shock PGV is small due to the small distance from the fault plane.As a result, the aftershock frequencies of the main shock PGV in this area are large.
The hazard map can be used to plan evacuation and infrastructure recovery after a large earthquake.

Aftershock Hazard Map
We propose a hazard map to efficiently use the aftershock hazard analysis for recovery activities.We focused on the short-term aftershock hazards and calculated the expected aftershock frequency at 113 sites (Figure 13).Figure 13a illustrates the results for μ + σ, whereas Figure 13b illustrates the results for μ. Figure 13 illustrates that the aftershock frequencies of the main shock PGV are considerably larger at sites that are closer to the southwestern area where the depth of the expanded fault plane is shallow.The maximum value is 0.52 for μ + σ and 0.20 for μ.The aftershock hazards increase near the fault plane.However, the increase in the PGV of the main shock is more significant, whereas the aftershock frequencies of the main shock PGV are negligible.In addition, the aftershock hazards, main shock PGV, and the aftershock frequencies are considerably smaller at sites located on the northeast side with the deep expanded fault plane.Conversely, the aftershock hazards are significant at sites located near the southwest side with the shallow expanded fault plane because the main shock PGV is small due to the small distance from the fault plane.As a result, the aftershock frequencies of the main shock PGV in this area are large.
The hazard map can be used to plan evacuation and infrastructure recovery after a large earthquake.

Aftershock Hazard Map
We propose a hazard map to efficiently use the aftershock hazard analysis for recovery activities.We focused on the short-term aftershock hazards and calculated the expected aftershock frequency at 113 sites (Figure 13).Figure 13a illustrates the results for µ + σ, whereas Figure 13b illustrates the results for µ. Figure 13 illustrates that the aftershock frequencies of the main shock PGV are considerably larger at sites that are closer to the southwestern area where the depth of the expanded fault plane is shallow.The maximum value is 0.52 for µ + σ and 0.20 for µ.The aftershock hazards increase near the fault plane.However, the increase in the PGV of the main shock is more significant, whereas the aftershock frequencies of the main shock PGV are negligible.In addition, the aftershock hazards, main shock PGV, and the aftershock frequencies are considerably smaller at sites located on the northeast side with the deep expanded fault plane.Conversely, the aftershock hazards are significant at sites located near the southwest side with the shallow expanded fault plane because the main shock PGV is small due to the small distance from the fault plane.As a result, the aftershock frequencies of the main shock PGV in this area are large.
The hazard map can be used to plan evacuation and infrastructure recovery after a large earthquake.

Use of the Aftershock Hazard Map
In this study, we propose a hazard map of the aftershock frequencies beyond the main shock PGV immediately following a large earthquake.However, if the aftershock frequencies following the main shock PGV are sufficiently small after a certain elapsed time after the main shock, it is possible to use hazard maps based on various indices.An example would be the standard seismic intensity, which is adopted by specific applications and can be applied to various decision-making situations, including rescue activities, evacuations, lifeline recoveries, nuclear facility operations, and traffic stops by setting the evaluation period to suit the intended use.

Load Combination of Aftershock and Tsunami for Tsunami-Resistant Designs
Immediately after the main shock of the 2011 Tohoku earthquake, many structures along the coast of the Tohoku region collapsed due to the simultaneous effect of the tsunami and the strong ground motions of the aftershocks.Even though the primary cause of the collapses of the structures was considered to be the tsunami load, the influence of the strong ground motions of the main shock and aftershocks cannot be ignored.In particular, critical structures, such as tsunami refuge buildings, seawalls, and nuclear power plants, which are located on the seashore, must be designed to withstand tsunamis and strong ground motions.It is necessary to quantitatively evaluate and consider the effects of load combinations of aftershock and tsunami.Even if the seismic design and tsunamiresistant design are appropriate, they may not account for the simultaneous occurrence of both aftershocks and tsunamis.In tsunami-resistant designs, we need to expect tsunamis to arrive at any time after the main shock ranging from a few minutes to a few hours based on the experience of the Tohoku earthquake.Therefore, both aftershocks and tsunamis should be considered simultaneously to design critical facilities located at the seashore.A practical assessment of the combined load of aftershocks and tsunamis is needed if both aftershocks and tsunamis occur simultaneously.

Load Combination Framework
The concept of the simultaneous occurrence of aftershocks and tsunamis [34] is depicted in Figure 14.Initially, the maximum tsunami arrival time at the evaluation site is expressed as the elapsed time ta and the tsunami completion time is expressed as te from the time of occurrence of the main shock.Therefore, the tsunami duration is represented as te − ta.According to the analysis of the Tohoku earthquake, ta was short and te − ta was a few minutes.Conversely, the load effects of the aftershocks occur immediately after the main shock and simultaneous effects of the aftershocks and the tsunami on structures occur within te − ta.Therefore, we propose the use of load and resistance combination factors in tsunami-resistant designs.The maximum height of the tsunami will arrive at the specified site at time ta after the main shock, and the maximum value is continuous within the

Use of the Aftershock Hazard Map
In this study, we propose a hazard map of the aftershock frequencies beyond the main shock PGV immediately following a large earthquake.However, if the aftershock frequencies following the main shock PGV are sufficiently small after a certain elapsed time after the main shock, it is possible to use hazard maps based on various indices.An example would be the standard seismic intensity, which is adopted by specific applications and can be applied to various decision-making situations, including rescue activities, evacuations, lifeline recoveries, nuclear facility operations, and traffic stops by setting the evaluation period to suit the intended use.

Load Combination of Aftershock and Tsunami for Tsunami-Resistant Designs
Immediately after the main shock of the 2011 Tohoku earthquake, many structures along the coast of the Tohoku region collapsed due to the simultaneous effect of the tsunami and the strong ground motions of the aftershocks.Even though the primary cause of the collapses of the structures was considered to be the tsunami load, the influence of the strong ground motions of the main shock and aftershocks cannot be ignored.In particular, critical structures, such as tsunami refuge buildings, seawalls, and nuclear power plants, which are located on the seashore, must be designed to withstand tsunamis and strong ground motions.It is necessary to quantitatively evaluate and consider the effects of load combinations of aftershock and tsunami.Even if the seismic design and tsunami-resistant design are appropriate, they may not account for the simultaneous occurrence of both aftershocks and tsunamis.In tsunami-resistant designs, we need to expect tsunamis to arrive at any time after the main shock ranging from a few minutes to a few hours based on the experience of the Tohoku earthquake.Therefore, both aftershocks and tsunamis should be considered simultaneously to design critical facilities located at the seashore.A practical assessment of the combined load of aftershocks and tsunamis is needed if both aftershocks and tsunamis occur simultaneously.

Load Combination Framework
The concept of the simultaneous occurrence of aftershocks and tsunamis [34] is depicted in Figure 14.Initially, the maximum tsunami arrival time at the evaluation site is expressed as the elapsed time t a and the tsunami completion time is expressed as t e from the time of occurrence of the main shock.Therefore, the tsunami duration is represented as t e − t a .According to the analysis of the Tohoku earthquake, t a was short and t e − t a was a few minutes.Conversely, the load effects of the aftershocks occur immediately after the main shock and simultaneous effects of the aftershocks and the tsunami on structures occur within t e − t a .Therefore, we propose the use of load and resistance combination factors in tsunami-resistant designs.The maximum height of the tsunami will arrive at the specified site at time t a after the main shock, and the maximum value is continuous within the tsunami duration and will be combined with the maximum value of the aftershock loads during the tsunami duration.tsunami duration and will be combined with the maximum value of the aftershock loads during the tsunami duration.

Analytical Model and Conditions
To evaluate the load combination factors of aftershocks and tsunamis in the Tohoku region, we assumed the main shocks to be the 2011 Tohoku earthquake (Mm = 9.0), the 1933 Sanriku earthquake (Mm = 8.1), and the 1896 Sanriku earthquake (Mm = 8.5).The evaluation sites are the K-NET locations of AOM012, IWT001, IWT002, IWT003, IWT005, and IWT007.The aftershock area is the expanded fault plane, which is enlarged due to the relationship between the main shock magnitude and the aftershock area [35].The evaluation sites and aftershock area are depicted in Figure 15.

Analytical Model and Conditions
To evaluate the load combination factors of aftershocks and tsunamis in the Tohoku region, we assumed the main shocks to be the 2011 Tohoku earthquake (M m = 9.0), the 1933 Sanriku earthquake (M m = 8.1), and the 1896 Sanriku earthquake (M m = 8.5).The evaluation sites are the K-NET locations of AOM012, IWT001, IWT002, IWT003, IWT005, and IWT007.The aftershock area is the expanded fault plane, which is enlarged due to the relationship between the main shock magnitude and the aftershock area [35].The evaluation sites and aftershock area are depicted in Figure 15.tsunami duration and will be combined with the maximum value of the aftershock loads during the tsunami duration.

Analytical Model and Conditions
To evaluate the load combination factors of aftershocks and tsunamis in the Tohoku region, we assumed the main shocks to be the 2011 Tohoku earthquake (Mm = 9.0), the 1933 Sanriku earthquake (Mm = 8.1), and the 1896 Sanriku earthquake (Mm = 8.5).The evaluation sites are the K-NET locations of AOM012, IWT001, IWT002, IWT003, IWT005, and IWT007.The aftershock area is the expanded fault plane, which is enlarged due to the relationship between the main shock magnitude and the aftershock area [35].The evaluation sites and aftershock area are depicted in Figure 15.The target structure is a ten-story reinforced concrete structure, with a width of 40 m, a depth of 20 m, and a height of 40 m.Based on the experience of the Tohoku earthquake [36], t a is set to 0 min, 30 min, and 60 min after the main shock and the tsunami duration is set to 30 min.The inundation height is assumed to remain constant during the entire duration of the tsunami.The limiting state is that in which the horizontal load of the aftershocks and the tsunami exceed the lateral resistance of the structure.
The design equation for the limiting state function [37] is where Q An is the conditional representative value of the aftershock load during the tsunami duration under the specified occurrence of the main shock, which we assume is a mean value based on the aftershock hazard analysis; Q Tn is the conditional representative value of the tsunami load under the specified occurrence of the main shock, which we assume is the mean value estimated based on the Abe law [38]; Q Mn is the conditional representative value of the main shock load under the specified occurrence of the main shock, which we assume is the median value estimated from the attenuation relation [27]; R n is the nominal value of the lateral resistance; γ A is the load factor for the aftershocks; γ T is the load factor for the tsunami; γ M is the load factor for the aftershocks; ϕ is the resistance factor; and r = Q An /Q Mn is the ratio of the nominal value of the load effect owing to aftershocks during the tsunami duration to the nominal value of the load effect owing to the main shock.
• Aftershock load First, we design the load effect of the aftershocks.We focus on the horizontal base shear force of n-story structures.According to the Recommendations for Loads on Buildings [39], the seismic base shear force Q A is approximately 2 times that of the ground seismic response: where w k is the weight of the floor k (kN), g is the gravitational acceleration [m/s 2 ], n is the number of stories of the target building, and PGA is the peak ground acceleration owing to the aftershocks [m/s 2 ], which is a random variable determined using Equation (13).Note that PGA is used instead of PGV.Therefore, we use the following attenuation relations (Equation ( 21)) for the calculation of PGA as an aftershock hazard and the site amplification R is assumed to be 1.4.
We use σ to model the error of the earthquakes.The mean (µ) and µ + σ of the aftershock PGA at each evaluation site are summarized in Table 3.In addition, the PGA distribution type is assumed to be lognormal.

• Tsunami load
Next, according to the tsunami-resistant design [40], the tsunami load Q T is where B is the width of the target structure [m], γ is the specific weight of water (kN/m 3 ), z 1 is the minimum height of the tsunami pressure-receiving surface (0 ≤ z 1 ≤ z 2 ) (m), z 2 is the maximum height of the tsunami pressure-receiving surface (z 1 ≤ z 2 ≤ 3h) (m), and h is the design tsunami inundation depth (m).The mean value of the design tsunami inundation depth is calculated by assuming a lognormal distribution for h and a standard deviation equal to 0.42 [38,41]: where h is the mean value of the design tsunami inundation depth [m], M m is the main shock magnitude, ∆ is the shortest tsunami propagation distance from the hypocenter to the observation point [km], and C is a correction term that is 0 toward the Pacific Ocean and 0.2 toward the Japan Sea.
• Load-resistance capacity The mean load-resistance capacity, R n is estimated using a target reliability index.Its coefficient of variation (COV) is assumed to be 0.2, and its distribution is assumed to be lognormal.

Analytical Method
The advanced first-order second moment (AFOSM) [42] is used as the method for the analysis.In this study, β T is the conditional target reliability index under the specified main shock occurrence.Therefore, it is necessary to determine β T by considering the occurrence frequency of the main shock.According to the Headquarters for Earthquake Research Promotion [43], the return period of the Tohoku earthquake is 600 years and that of the Sanriku earthquake is 97 years.Therefore, we set the conditional target reliability index β T to take values from 0.0 to 1.0 and calculated the design value of the aftershock and tsunami load combination.In addition, the load combination factor was estimated from the ratio of the representative and design values.

Analytical Results
Figure 16 illustrates the aftershock and tsunami load combination factor at IWT007, whereas Figure 17 depicts the results at all the evaluation sites.
According to Figure 16, regardless of the main shock, the tsunami load factor increases with increasing target reliability index.However, the aftershock load factor remains nearly constant, and the resistance factor decreases.The aftershock load factor tends to decrease with increasing elapsed time from the main shock.The difference between 30-60 min and 60-90 min is smaller than that between 0-30 min and 30-60 min because the occurrence rate of the aftershocks decreases with increasing elapsed time.However, the effect of the elapsed time on the tsunami load factor and the resistance factor is negligible.
Figure 16a,c,e are compared to Figure 16b,d,f to confirm the effect of variations in the aftershock hazard analysis.The variation is significant in the aftershock load factor, which was approximately three times the variation.Simultaneously, we compared the 2011 Tohoku earthquake (M m = 9.0), the 1896 Sanriku earthquake (M m = 8.5), and the 1933 Sanriku earthquake (M m = 8.1) to confirm the effect of the magnitude of the main shock.It is difficult to make a direct comparison because the shortest distance from an evaluation site to the fault plane is different at different sites.The aftershock load factor and tsunami load factor of the 1933 Sanriku earthquake were affected negligibly compared to the factors for the 2011 Tohoku and 1896 Sanriku earthquakes.In addition, the aftershock load factor decreases and the tsunami load factor increases with increasing elapsed time after the main shock.The variation in the tsunami load factor for the 1933 Sanriku earthquake is remarkable.The main shock magnitude of the 1933 Sanriku earthquake is the smallest of the target earthquakes.Therefore, the load effect of the tsunami is considerably smaller than the load effect of the aftershocks.

Load Combination Factor for Tsunami-resistant Designs
We consider the very rare condition that large subduction earthquakes (Tohoku and Sanriku earthquakes) occur.Therefore, the conditional target reliability index is set to 0 and 0.5.The LRFD format for the tsunami-resistant design of the conditional target reliability index is shown in Table 4. Here, we summarize the results at all the evaluation sites.The relationship between the fault distance and the load combination factor for the target reliability index is illustrated in Figure 17.The aftershock load factor and the tsunami load factor for each target reliability index are set high and the resistance factor is set low to evaluate the most conservative load combination factor, including the effects of elapsed time.Figure 16 indicates that the load combination factors are nearly constant regardless of the fault distance.However, the aftershock load factor for the 1933 Sanriku earthquake is the smallest of the three earthquakes.This is because the aftershock occurrence rate is related to the main shock magnitude.Conversely, the resistance factor is nearly constant regardless of the main shock magnitude.Moreover, the aftershock load factor increases when the aftershock hazard increases owing to the variation.However, there is little change in the tsunami load factor owing to the variation in the aftershock hazard for the 2011 Tohoku earthquake and the tsunami load factor decreases with increasing aftershock hazard for the 1896 and 1933 Sanriku earthquakes.Presumably, this is because the tsunami to aftershock load ratios in the 1896 and 1933 Sanriku earthquakes are smaller than that of the 2011 Tohoku earthquake.

Load Combination Factor for Tsunami-resistant Designs
We consider the very rare condition that large subduction earthquakes (Tohoku and Sanriku earthquakes) occur.Therefore, the conditional target reliability index is set to 0 and 0.5.The LRFD format for the tsunami-resistant design of the conditional target reliability index is shown in Table 4.

Conclusions
In this study, we introduced a probabilistic aftershock hazard analysis that is expected to be useful to develop plans for recovery activities after any large earthquake in the future.The probabilistic aftershock hazard analysis depends on the elapsed time after the main shock using model parameters based on compiled statistics of historical aftershock data.The proposed model was validated by comparing the results with observed aftershock ground motions during the 2011 Tohoku earthquake.In addition, we conducted a sensitivity analysis and confirmed the effects of variations in each parameter.Then, two engineering applications of the proposed approach for probabilistic aftershock hazard analysis were shown for demonstration purposes.It was proposed that aftershock hazard maps be used for the effective planning of the recovery activity for future large earthquakes that may occur in Japan.It was also proposed that aftershock hazard analysis be used for the development of load combination equations of LRFD considering the simultaneous occurrence of tsunamis and aftershocks, which is considered to be useful for the design of tsunami evacuation buildings.
In the future, we will consider the effects of aftershocks for safety evaluations of nuclear facilities to contribute to the advancement of probabilistic seismic risk analyses of nuclear facilities.

Figure 2 .
Figure 2. Aftershock distribution due to the 2011 Tohoku earthquake.

Figure 2 .
Figure 2. Aftershock distribution due to the 2011 Tohoku earthquake.

Figure 4 .
Figure 4. Relationship between Mm and (a) the aftershock area; (b) the number of aftershocks (N); (c) the b-value after 90 days; (d) the p-value; and (e) the D1-value.Figure 4. Relationship between M m and (a) the aftershock area; (b) the number of aftershocks (N); (c) the b-value after 90 days; (d) the p-value; and (e) the D 1 -value.

Figure 4 .
Figure 4. Relationship between Mm and (a) the aftershock area; (b) the number of aftershocks (N); (c) the b-value after 90 days; (d) the p-value; and (e) the D1-value.Figure 4. Relationship between M m and (a) the aftershock area; (b) the number of aftershocks (N); (c) the b-value after 90 days; (d) the p-value; and (e) the D 1 -value.

Figure 5 .
Figure 5. (a) Gutenberg-Richter (GR) model at different elapsed times (T 1 ) for the 2011 Tohoku earthquake and (b) b(T 1 ) versus T 1 for past major earthquakes.

Figure 6 .
Figure 6.Setting of the aftershock region.

Figure 6 .
Figure 6.Setting of the aftershock region.

Figure 6 .
Figure 6.Setting of the aftershock region.

Figure 12 22 Figure 11 .
Figure 12 illustrates the short-term aftershock hazards and the main shock PGV.The predicted main shock PGV is 48 cm/s, and the maximum aftershock frequency of three days is approximately 0.1 for µ + σ and approximately 0.02 for µ.Geosciences 2018, 8, 1 12 of 22

Figure 11 .
Figure 11.Site and aftershock region of the 1923 Kanto earthquake.

Figure 11 .
Figure 11.Site and aftershock region of the 1923 Kanto earthquake.

Figure 13 .
Figure 13.Hazard map of the short-term aftershock frequency of the main shock PGV: (a) μ + σ and (b) μ.

Figure 13 .
Figure 13.Hazard map of the short-term aftershock frequency of the main shock PGV: (a) µ + σ and (b) µ.

Figure 14 .
Figure 14.Concept of aftershock and tsunami load combination.

Figure 15 .
Figure 15.Evaluation sites and the aftershock area.The target structure is a ten-story reinforced concrete structure, with a width of 40 m, a depth of 20 m, and a height of 40 m.Based on the experience of the Tohoku earthquake [36], ta is set to 0 min, 30 min, and 60 min after the main shock and the tsunami duration is set to 30 min.The inundation

Figure 14 .
Figure 14.Concept of aftershock and tsunami load combination.

Figure 14 .
Figure 14.Concept of aftershock and tsunami load combination.

Figure 15 .
Figure 15.Evaluation sites and the aftershock area.The target structure is a ten-story reinforced concrete structure, with a width of 40 m, a depth of 20 m, and a height of 40 m.Based on the experience of the Tohoku earthquake [36], ta is set to 0 min, 30 min, and 60 min after the main shock and the tsunami duration is set to 30 min.The inundation

Figure 15 .
Figure 15.Evaluation sites and the aftershock area.
the variation in the aftershock hazard for the 2011 Tohoku earthquake and the tsunami load factor decreases with increasing aftershock hazard for the 1896 and 1933 Sanriku earthquakes.Presumably, this is because the tsunami to aftershock load ratios in the 1896 and 1933 Sanriku earthquakes are smaller than that of the 2011 Tohoku earthquake.

Table 1 .
Aftershock occurrence parameters for past earthquakes.

Table 3 .
Summary of the input values in the reliability analysis.

Table 4 .
Summary of the load and resistance factor design factor.