Abstract

The 2016 Meinong Earthquake hit southern Taiwan and many shallow foundation structures were damaged due to soil liquefaction. In response, the government initiated an investigation project to construct liquefaction potential maps for metropolitans in Taiwan. These maps were used for the preliminary safety assessment of infrastructures or buildings. However, the constructed liquefaction potential map used the pseudo-probabilistic approach, which has inconsistent return period. To solve the inconsistency, the probabilistic liquefaction hazard analysis (PLHA) was introduced. However, due to its complicated calculation procedure, PLHA is not easy and convenient for engineers to use without a specialized program, such as in Taiwan. Therefore, PLHA is not a popular liquefaction evaluation procedure in practice. This study presents a simple PLHA program, HAZ45PL Module, customized for Taiwan. Sites in Tainan City and Yuanlin City are evaluated using the HAZ45PL Module to obtain the hazard curve and to construct the liquefaction probability map. The liquefaction probability map provides probabilities of different liquefaction potential levels for engineers or owners to assess the performance of an infrastructure or to design a mitigation plan.

1. Introduction

The 2016 Meinong Earthquake (Mw = 6.5) struck south Taiwan on February 6, 2016. The strong ground motion caused huge loss of lives (117 people dead and 504 people injured) [1]. Many shallow foundation structure in Tainan area were damaged due to soil liquefaction which thus gained considerable public attention. In response, the government initiated a project to evaluate and construct the liquefaction potential map for the urban areas of Taiwan. This map was constructed based on the liquefaction potential index (LPI) [2] which is calculated using the factor of safety against liquefaction (FSL) determined from a simplified liquefaction evaluation procedure. Engineers and owners use the liquefaction potential map to perform a preliminary safety evaluation on infrastructures or buildings and to determine the need of a mitigation plan.

The procedures for the LPI calculation in Taiwan are (1) select the design level (e.g., 475-year return period) peak ground acceleration (PGA); (2) pair PGA with a single earthquake magnitude (Mw) obtained from the deaggregation process of the probabilistic seismic hazard analysis (PSHA) at the same return period; (3) decide a representative water table level; (4) calculate FSL using a simplified liquefaction evaluation procedure with the selected PGA, Mw, and water table level; and (5) estimate LPI according to the depth and FSL. The estimated LPI or other related hazards are assumed to correspond to the same return period of PGA. This approach is named the “pseudo-probabilistic” approach [35].

In the pseudo-probabilistic approach, LPI at a site is estimated from a selected pair of PGA and Mw , and Mw is assumed to have the same return period with the selected PGA. However, from the PSHA concept, a site could be subjected to numerous pairs of PGA and Mw. Each pair has its own rate of occurrence. Therefore, LPI estimated from the pseudo-probabilistic approach and the selected PGA will not have a consistent return period as expected [2, 68]. As a solution to this inconsistency issue, a probabilistic approach, the performance-based earthquake engineering (PBEE) or the probabilistic liquefaction hazard analysis (PLHA), was introduced to incorporate the occurrence of the ground motion (PGA and Mw pairs) obtained from PSHA into the calculation of LPI or other liquefaction hazards [911].

In PLHA, engineers no longer need to select one representative pair of PGA and Mw as analysis inputs. All scenarios (PGA and Mw pairs) and their occurrence frequencies are considered in PSHA by accounting for the uncertainties of earthquake occurrence, ground-motion attenuation, and site effects. The calculation of PSHA is complicated, and thus PLHA is not easy and convenient for engineers to use without a specialized program. Recently, several simplified performance-based methods [8, 9, 1214] were established for engineers to perform PLHA in an approximated way. Available seismic hazard curves and deaggregation results provided by the United States Geological Survey (USGS) are used to obtain the occurrence frequency of PGA and Mw pairs in these simplified methods.

In Taiwan, without tools for PLHA, the current liquefaction analysis approach is still based on the pseudo-probabilistic approach. In addition, the simplified performance-based method mentioned above cannot be utilized because the seismic hazard curve and the deaggregation result are not yet available. This study presents a simple PLHA program, HAZ45PL Module, customized for engineers in Taiwan to obtain the LPI hazard curve. Moreover, a new application of PLHA results, the probability of exceedance map of LPI, is demonstrated. The probability information can help engineers, owners, or government agencies to carry out a performance assessment on infrastructures or prioritize a mitigation plan. In addition, because of the importance of the water table level on LPI, the HAZ45PL Module includes the uncertainty of the seasonal water table level in the LPI calculation.

2. Methodology of HAZ45PL Module

The HAZ45PL Module consists of HAZ45PL and R language code [15]. HAZ45PL is an executable file calculating LPI of each soil layer based on a PGA and Mw pair obtained from earthquake scenarios in PSHA. HAZ45PL is modified from HAZ45, an open source FORTRAN code developed by Dr. Norm Abrahamson for PSHA [16]. The steps of HAZ45PL Module to construct the LPI hazard curve are explained as follows: (1) generate an earthquake scenario following procedures in HAZ45; (2) provide Mw and the occurrence frequency of the earthquake scenario; (3) calculate PGA using reasonable ground-motion prediction equations (GMPE) for the earthquake scenario; (4) calculate FSL of each soil layer using PGA and Mw from Steps (2) and (3); (5) estimate LPI according to FSL from Step (4); (6) export LPI and the occurrence frequency of each soil layer; (7) repeat Steps (1)–(6) for earthquake scenarios from all seismic sources; (8) use R language code as a post-process platform to process LPI and the occurrence frequency of the corresponding PGA and Mw pairs and generate the LPI hazard curve.

The input data files of the HAZ45PL Module include the seismic source and the soil property data files. Input data of seismic sources are obtained from Central Geological Survey [17], Cheng et al. [18], and Cheng et al. [19]. In Taiwan, seismogenic sources can be categorized as follows: (1) shallow regional source, (2) deep regional source, (3) crustal active fault source, (4) subduction interface source, and (5) subduction intraslab source. Because the contribution of the deep regional source is much less comparing to other seismic sources [20], in order to shorten the analysis time, the deep regional source is not included in the current HAZ45PL Module. The needed input data of soil properties are unit weight, SPT N, fines content, plastic index, and clay content.

GMPEs used in the current version of the HAZ45PL Module are built in equations for regional, crustal, and subduction seismic sources. Two crustal GMPEs, Lin11 [21] and CY14 [22], are selected for regional and crustal sources. Lin11 was built using a GMPE model for the tectonic setting of Taiwan and the mostly used GMPE in Taiwan. Lin11 includes the site response effect and the hanging wall and footwall effect on GMPE. CY14 is a GMPE from NGA-West2 selected for considering the ground-motion prediction model uncertainty (epistemic uncertainty) in the HAZ45PL Module. CY14 was built using a GMPE model for the tectonic setting similar to Taiwan and also includes the site condition effect and the hanging wall and footwall effect on GMPE. LL08 [23] and BCH16 [24] are two GMPEs selected for subduction sources. LL08 was built using data of subduction earthquakes in Taiwan and is also the mostly used GMPE in Taiwan. BCH16 is for worldwide subduction earthquakes and selected for considering the ground-motion prediction model uncertainty (epistemic uncertainty) in the HAZ45PL Module. In the current version of HAZ45PL Module, in order to simplify the analysis process, the variation of GMPE on PGA is estimated using the three-point approximation method [25].

Five frequently used simplified liquefaction evaluation procedures for FSL calculation in Taiwan are included in the HAZ45PL Module as follows: National Center for Earthquake Engineering Research (NCEER) procedure [26], Hyperbolic Function (HBF) procedure [27], Japan Road Association (JRA) procedure [28], Tokimatsu and Yoshimi (T-Y) procedure [29], and Architecture Institute of Japan (AIJ) procedure [30]. FSL estimated from the above procedures are used to obtain LPI following equations from Iwasaki et al. [2]:where i is the layer number (total number of soil layers is n), (Fi = 0 for ), Wi = 10–0.5 Zi is the weighting factor for ith layer and Zi is the depth (in meters) of ith layer, and Hi is thickness (in meters) of ith layer. The liquefaction potential is categorized in four levels: (1) LPI = 0, the liquefaction potential is very low; (2) 0 < LPI ≤ 5, the liquefaction potential is low; (3) 5 < LPI ≤ 15, the liquefaction potential is high; and (4) 15 < LPI, the liquefaction potential is very high. The uncertainty of the evaluation procedure of FSL on LPI is considered via the logic tree approach.

Apart from ground motion and soil properties, the variation of water table level is another important factor affecting the liquefaction potential. In the HAZ45PL Module, the variation of the water table level is included via the logic tree approach. The logic tree used in the HAZ45PL Module is shown in Figure 1. Weights and variations of GMPE are fixed while those of the water table level and the liquefaction potential calculation procedure can be adjusted by users.

3. Applications of HAZ45PL Module

Applications on a single location using PLHA results are demonstrated in Tainan City, and applications of PLHA results on liquefaction potential maps are demonstrated in Changhua County.

3.1. Tainan City

First, the HAZ45PL Module is applied to two sites (Annan and Xinshi Sites) in Tainan City, where a detailed investigation with 18 boreholes was performed after the 2016 Meinong Earthquake. In this study, six boreholes (BH-A1 to BH-A3 of Annan Site and BH-S1 to BH-S3 of Xinshi Site) are used for analyses. Figure 2 and Table 1 show the site locations and their information, respectively. Soil layers of Annan and Xinshi Sites are mostly Silty Sand (SM) layers or Silt (ML) layers with few Clay (CL) layers, as shown in Figure 3. Shallow SM or ML layers (<10 m) have relatively low SPT N values (<10) and water table levels are high (between 0.6 and 0.9 m below the ground surface). These make Annan and Xinshi Sites susceptible to liquefaction during earthquakes. Detailed investigation information refers to NCEER [31].

Hazard curves of LPI for boreholes in Annan and Xinshi Sites are shown in Figure 4 and the 475-year return period LPIs are listed in Table 2. Results of the pseudo-probabilistic approach using PGA and Mw pairs in Table 1 are also marked for comparisons. The water table used here is the water table during the site investigation, as listed in Figure 3. Effects of water table on LPI are discussed later. All listed liquefaction potential calculation procedures in Figure 1 are used and equal weights are assigned to each procedure in this analysis for simplicity. LPI hazard curves at the same site from different boreholes (BH-A1 to BH-A3 of Annan Site or BH-S1 to BH-S3 Xinshi Site) show different levels of variations because of the spatial variation and randomness of soil properties. Table 2 shows that LPIs at 475-year return period (475 yr RP) differ from those of the pseudo-probabilistic approach with 475 yr RP PGA. In addition, return periods of LPIs of the pseudo-probabilistic approach are quite different between boreholes. The inconsistencies of LPIs and return periods are caused by the incapability of the pseudo-probabilistic approach to include the occurrence frequency of the PGA and Mw pairs. PGA and Mw pairs and their rate contributions to PGA at 475 yr RP at two sites are deaggregated and shown in Figure 5. Clearly, each PGA and Mw pair contributes differently to the total hazard and cannot be represented by only one PGA and Mw pair, such as used in the pseudo-probabilistic approach. The only way to obtain the correct LPI at 475 yr RP is to include the occurrence frequency of all PGA and Mw pairs in the LPI estimation.

Next, BH-S1 is used to demonstrate the effects of the water table on LPI. As a continuous random variable, the variation of the water table also can be approximated using the three-point approximation method [25]. Weights are assigned to different water table levels (shown in Figure 1) as follows: (1) High Water Table: average water table (sets at the ground surface) of the high water table period with a weighting of 0.185. (2) Yearly Average: long term water table (sets at the water table in Figure 3) with a weighting of 0.630. (3) Low Water Table: average water table (sets at 2.0 m below the ground surface) of the low water table period with a weighting of 0.185. Figure 6 and Table 3 show LPI hazard curves and the pseudo-probabilistic approach results of BH-S1 at different water table levels. It is obvious that the water table level has a significant impact on LPI and the return period. In the current practice, LPI is estimated using the pseudo-probabilistic approach with the high water table because engineers presume that the LPI at this water table level is relatively conservative. However, the conservativeness cannot be reasonably quantified using the pseudo-probabilistic approach. In PLHA, the water table fluctuation can be considered via a probabilistic approach and the LPI hazard curve can reasonably reflect effects of the water table fluctuation.

Other than the inconsistency in return period, another disadvantage of the pseudo-probabilistic approach is the provision of one LPI value (Tables 2 or 3) for engineers or owners to evaluate the need and the priority of mitigation schemes for infrastructures. For example, BH-S1 and BH-A1 have LPI = 20.6 and 36.4 from the pseudo-probabilistic approach. When LPI > 15, the liquefaction potential is very high and the mitigation plan on soil liquefaction is needed [2]. Therefore, the mitigation is necessary for infrastructures in both sites. On the contrary, PLHA can provide not only LPI but also the return period which can be converted to the probability of exceedance for ≥LPI in certain exposure period (e.g., lifetime of infrastructure) using the Poisson process. Table 4 lists the probabilities of LPI > 15 for BH-S1 and BH-A1 during different exposure periods as a reference. Engineers or owners can easily and rationally make the final design decision according to the probability information.

3.2. Changhua County

The HAZ45PL Module is applied to Yuanlin City, Changhua County (Figure 7), to demonstrate the construction of the liquefaction potential contour maps (Figures 8 and 9) using the Kriging method. In Figure 8, liquefaction potential levels follow the guidelines in Iwasaki et al. [2]. Comparing Figures 8(a) and 8(b), the pseudo-probabilistic approach obviously underestimates the liquefaction potential. In Figure 8(b), the east side of Yuanlin City is at a very high liquefaction potential. For future infrastructure or mitigation projects in this area, engineers and owners must optimize the mitigation plan and prioritize the site. As mentioned earlier, PLHA can provide the probability of exceedance in different exposure periods. Figure 9 shows probabilities of exceedance having LPI ≥5 (high liquefaction potential) and LPI ≥ 15 (very high liquefaction potential) in 50 years. Inside the very high liquefaction potential level area (the east side of Yuanlin City), probabilities in Figure 9(b) are divided into three groups 20%–40%, 40%–60%, and 60%–80%. Comparing to use only Figure 8(b) as a reference, engineers or owners can make reasonable choices on the construction site or the mitigation plan for infrastructures or other buildings based on the probability information shown in Figure 9(b). The probability information also provides a valuable reference for the government to make an appropriate emergency response and evacuation plan.

4. Conclusions

After the 2016 Meinong Earthquake, in response of the public concern on the soil liquefaction hazard, the central government initiated a project to evaluate and construct the liquefaction potential map. The common procedure used for constructing these maps is called the pseudo-probabilistic approach, which uses a pair of PGA and Mw at the selected return period with a “representative” water table level to estimate LPI for the map. The use of the pseudo-probabilistic approach has two weakness points: (1) the inconsistent return period of the estimated LPI and (2) the difficulty to select a representative water table level. PLHA was introduced to incorporate the occurrence of the ground motion (PGA and Mw pairs) into LPI calculation to solve the inconsistency issue. In addition, the difficulty of water table level selection can also be solved in PLHA via the logic tree or the probability density function. However, given the lack of an appropriate tool, PLHA is not popular in Taiwan. In this study, a simple PLHA program, HAZ45PL Module, is developed for engineers in Taiwan to obtain the LPI in a probabilistic way.

The HAZ45PL Module is applied to two selected sites to demonstrate the use of PLHA results, the hazard curve of LPI and the liquefaction probability map. With the LPI hazard curve and the Poisson process, PLHA can provide the probability of exceedance information to construct liquefaction probability maps, which enable engineers or owners to easily assess the performance of an infrastructure or to design a mitigation plan. The government can make an emergency response and evacuation plan according to the liquefaction probability map too. The liquefaction probability map was not introduced or used in previous PLHA studies.

In the current HAZ45PL Module, to keep the calculation simple and easy, the variability of the liquefaction resistance of soil is not included and the variations of PGA from GMPE are estimated using the three-point approximation. The ongoing modification of the HAZ45PL Module addresses these issues. Moreover, using the liquefaction potential only can indirectly assess the infrastructure or facility safety. In the next version of the HAZ45PL Module, the probability of liquefaction-related hazards (e.g., settlement of structure and uplift of underground structure) can be incorporated to provide a direct assessment.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by Central Geological Survey, MOEA, Taiwan, and Sinotech Engineering Consultants, Inc. The authors would like to thank staff of Central Geological Survey for providing necessary assistance and valuable expertise.