Land Damage Mapping and Liquefaction Potential Analysis of Soils from the Epicentral Region of 2017 Pohang Mw 5.4 Earthquake, South Korea

Studies on earthquake-induced liquefaction and identification of source unit for causing liquefaction have been a major concern in sustainable land use development especially in low to moderate seismic areas. During the 2017 Mw 5.4 Pohang earthquake, widespread liquefaction was reported around the Heunghae basin, which was the first ever reported case of liquefaction in the modern seismic history of Korea. The epicentral area is one of the major industrial hubs along the SE Korean Peninsula with no detailed liquefaction hazard map. The purpose of this study was to determine the land damage classification on the basis of surface manifestation of liquefaction features and carry out detailed liquefaction potential analysis to delineate the depth of liquefiable soil. This will eventually support developing a liquefaction hazard zonation map and sustainable development of infrastructure to minimize earthquake damages. In this present study, the southern part of the Heunghae basin, which has more field evidences of liquefaction than the northern part, was taken for detailed liquefaction analysis. From the detailed analysis, it was observed that the soils from 1.5 to 15 m depth with the probability of liquefaction varying from 2 to 20 are prone to liquefaction. On the basis of land damage pattern, the epicentral area falls in orange to red zone, which means the necessity of further detailed liquefaction analysis. This study urges more detailed liquefaction zonation should be carried out for the epicentral area and liquefaction hazard should be included in the multi-hazard map in the future for the sustainable land use planning.


Introduction
During earthquake shaking or other rapid loadings, mostly saturated sandy/silty soil loses its strength and stiffness and behaves as a liquid which is known as liquefaction [1][2][3]. The repeated/cyclic loading under undrained conditions leads to a decrease in effective stress and a rise in pore water pressure of the soil [4,5]. Along with earthquake magnitude, several other factors affect the liquefaction susceptibility of the soil, such as the age of the soil, sedimentation process, the depth of the water table, density, burial depth, ground slope, grain size distribution, etc. [6,7]. The basic understanding of the relationship between cyclic loading during an earthquake and soil liquefaction has been largely derived from laboratory studies, such as cyclic triaxial tests, shake table tests, resonant column tests, centrifuge modeling, etc. [8][9][10][11][12][13][14].
of surface manifestation of liquefaction features to make the land damage map. Since there are no specific criteria available for land damage pattern categorization for the Korean Peninsula, we have adopted the methodology applied to the Canterbury basin. Moreover, we have considered a 10 km line along the southern part of the basin, where more localized liquefaction features were observed during the earthquake.

Geological Settings
The southern part of the Korean Peninsula is located on the Eurasian Plate and its tectonic activities are controlled by the ongoing subduction of the adjoining Pacific and Philippine Sea Plates. Moreover, the collision between the Indian Plate with Eurasian Plate controls the tectonism of the Korean Peninsula [47][48][49][50][51]. The Cenozoic deformation along the southeastern part of the Korean Peninsula is accommodated by the Yangsan-Ulsan Fault system. The Oligocene rifting of the Korean Peninsula in NE-SW direction has resulted in the development of several normal faults and extensional basins [52]. The major basin is known as Pohang basin, which is subdivided into several sub-basins. The Heunghae basin is one of the sub basins which was affected by the 2017 Pohang earthquake (Mw 5.4) ( Figure 1). There have been several reports of Quaternary faulting with surface ruptures along the Yangsan-Ulsan Fault System [51,53]. The distribution of Quaternary faults and mechanisms involved are consistent with the maximum principal stress direction (E-W; ENE-WSW) during the Quaternary period [54][55][56]. The recorded seismic data along this area suggest that this area has experienced four moderate magnitude earthquakes (Mw 5-6) during the last 14 years, including the 2016 Gyeongju earthquake (Mw 5.5) and the 2017 Pohang earthquake (Mw 5.4) [47][48][49].  Pohang earthquake and its focal mechanism (modified from [47]).

2017 Pohang Earthquake and Observed Liquefaction Features
The Korean Meteorological Administration (KMA), United States Geological Survey (USGS), Korea Institute of Geosciences and Mineral Resources (KIGAM) reported the main shock (Mw 5.4) and four major aftershocks (˃Mw 3.5). The focal depth of the main shock was estimated to be 4.5 km with the depth of aftershock sequences varying from 3 to 7 km [47,49,57]. The focal mechanism solution of the earthquake indicates that the main shock was an obliqueslip fault (reverse and dextral strike-slip) with NE strike and dip of approximately 60°-70°. Post-earthquake InSAR data analysis suggests that the area underwent 15 cm of coseismic slip during the main shock [49]. The earthquake caused damages to more than 2,500 houses, 227 school buildings, 11 bridges, and many roads around the epicentral area. The estimated total cost of loss is about 52 million USD.
The post-earthquake reconnaissance surveys carried out by several authors reported more than 600 cases of liquefaction along with several examples of lateral spreading [47][48][49][50]. Field photos showing evidence of liquefaction and lateral spreading are shown in Figure 2.  In some places, the liquefaction features were even more than 2 m in diameter ( Figure 2). The shaking was so intense that there was a report of liquefaction of gravelly soil along the dried river beds and rice farms [47]. Due to the presence of soft ground and shallow focal depth of the earthquake, there were several cases of ground settlements of up to 0.39 m. Some residential buildings were tilted by an angle of 1.6 • with differential settlement of 0.15 m [47,50]. The grain size analysis was done as per the American Standard for Testing Materials (ASTM) D6913 [58] for the ejected soil from the sand boil, which is classified as mostly silty sand, silt, and sandy ground ( Figure 2) and falls under liquefiable soil category [47,59,60]. Figure 3 depicts some of the typical damaged buildings during the earthquake. Considering the observed liquefaction, ground cracks, and lateral spreading around the epicentral area, we have proposed a local scale land damage zone as per the guideline for geotechnical engineering practice in New Zealand [19,61]. Table 1 shows the land damage categories in the epicentral area of the Pohang earthquake and Figure 4 shows the land damage distribution with locations of surface manifestations of liquefaction during the earthquake. On the basis of the land damage distribution, it was observed that most of the damaged lands are categorized as moderate to major land damage and are confined between the Gokgang and Chogok Rivers. In some places, the liquefaction features were even more than 2 m in diameter ( Figure 2). The shaking was so intense that there was a report of liquefaction of gravelly soil along the dried river beds and rice farms [47]. Due to the presence of soft ground and shallow focal depth of the earthquake, there were several cases of ground settlements of up to 0.39 m. Some residential buildings were tilted by an angle of 1.6° with differential settlement of 0.15 m [47,50]. The grain size analysis was done as per the American Standard for Testing Materials (ASTM) D6913 [58] for the ejected soil from the sand boil, which is classified as mostly silty sand, silt, and sandy ground ( Figure  2) and falls under liquefiable soil category [47,59,60]. Figure 3 depicts some of the typical damaged buildings during the earthquake. Considering the observed liquefaction, ground cracks, and lateral spreading around the epicentral area, we have proposed a local scale land damage zone as per the guideline for geotechnical engineering practice in New Zealand [19,61]. Table 1 shows the land damage categories in the epicentral area of the Pohang earthquake and Figure 4 shows the land damage distribution with locations of surface manifestations of liquefaction during the earthquake. On the basis of the land damage distribution, it was observed that most of the damaged lands are categorized as moderate to major land damage and are confined between the Gokgang and Chogok Rivers. The land damage pattern suggests if a similar or higher magnitude earthquake occurs in the near future, a similar scale of land will undergo liquefaction. This kind of repeated liquefaction was observed during the Christchurch earthquake and Emilia earthquake. The recent paleo-liquefaction studies from the epicentral area of the Pohang earthquake [48] also suggested that the area have experienced liquefaction in the past. The site which shows liquefaction evidence during this earthquake has undergone similar liquefactions within 730 to 130 years BP [48]. Therefore, detailed The land damage pattern suggests if a similar or higher magnitude earthquake occurs in the near future, a similar scale of land will undergo liquefaction. This kind of repeated liquefaction was observed during the Christchurch earthquake and Emilia earthquake. The recent paleo-liquefaction studies from the epicentral area of the Pohang earthquake [48] also suggested that the area have experienced liquefaction in the past. The site which shows liquefaction evidence during this earthquake has undergone similar liquefactions within 730 to 130 years BP [48]. Therefore, detailed liquefaction analysis, specifically the determination of liquefiable soil depth, will help in reducing liquefaction-related hazard for future construction projects on liquefiable soil around the epicenter.

Ground Motion Data
The 2017 November Pohang earthquake data were taken for the liquefaction analysis. The ground motion parameters recorded at the closest seismic station PHA2 were taken for the analysis. The acceleration time history data for the Pohang earthquake used for the analysis are given in Figure 5. The acceleration values recorded at PHA2 seismic station (0.268 g) were taken for the liquefaction analysis. A detailed liquefaction analysis was carried out using the geotechnical borehole data ( Figure 4) collected from the Integrated Data Base Center of National Geotechnical Information, Korea [62].

Site Investigation and Liquefaction Assessment
From the spatial distribution of the liquefaction features in the epicentral area, it can be inferred that the liquefaction features are mostly concentrated on the southern part of the Heunghae basin [47]. In order to analyze the subsurface soil properties of liquefied sites and determine the thickness of the liquefaction-prone soil around the Heunghae basin, we collected the geotechnical investigation reports from the Integrated Data Base Center of National Geotechnical Information, Korea [62]. The integrated database compiled all the soil classification data and Standard Penetration Test (SPT-N) data sets from the completed or ongoing road, railways, and residential construction projects. The standard penetration test is one of the most widely used in-situ tests carried out for major geotechnical projects due to its simple nature.
In this method, a 65 cm long thick-walled (5 cm ) sampler was driven into the ground at the bottom of the borehole with a sliding hammer of 63.5 kg, falling from a distance of 76 cm [63,64]. The number of blows needed for each of the 15 cm sampling tube penetration is recorded. The sum of the number of blows required for the second and third 30 cm of penetration is known as standard penetration resistance (SPT-N) values. The SPT-N blow count can indicate the density of the soil [63,64]. Typical geotechnical data used for the liquefaction analysis is shown in Table 2. Figure 6 shows the grain size distribution graph for the ejected soil from the sand boils.

Site Investigation and Liquefaction Assessment
From the spatial distribution of the liquefaction features in the epicentral area, it can be inferred that the liquefaction features are mostly concentrated on the southern part of the Heunghae basin [47]. In order to analyze the subsurface soil properties of liquefied sites and determine the thickness of the liquefaction-prone soil around the Heunghae basin, we collected the geotechnical investigation reports from the Integrated Data Base Center of National Geotechnical Information, Korea [62]. The integrated database compiled all the soil classification data and Standard Penetration Test (SPT-N) data sets from the completed or ongoing road, railways, and residential construction projects. The standard penetration test is one of the most widely used in-situ tests carried out for major geotechnical projects due to its simple nature.
In this method, a 65 cm long thick-walled (5 cm) sampler was driven into the ground at the bottom of the borehole with a sliding hammer of 63.5 kg, falling from a distance of 76 cm [63,64]. The number of blows needed for each of the 15 cm sampling tube penetration is recorded. The sum of the number of blows required for the second and third 30 cm of penetration is known as standard penetration resistance (SPT-N) values. The SPT-N blow count can indicate the density of the soil [63,64]. Typical geotechnical data used for the liquefaction analysis is shown in Table 2. Figure 6 shows the grain size distribution graph for the ejected soil from the sand boils.  Considering the clustered liquefaction features and closest available geotechnical datasets, we took a 10.5 km transect to determine the liquefaction potential in terms of thefactor of safety against liquefaction (FSL) of the soil. The variation of the depth of liquefiable soil, estimation of post liquefaction settlement and liquefaction potential index (IL) along the southern part of the Heunghae basin experienced more severe liquefaction than the northern part of the basin. Since there was no detailed analysis of liquefaction potential of soils from the epicentral area, this study will help in determining the liquefiable soil depth and the possible post liquefaction settlement in the area, which is one of the most important factors in seismic-resistant building design.

Result of Selected SPT in the Heunghae Basin
Geotechnical investigation results from 41 boreholes covering the lateral extent of the Heunghae basin were collected. Around the Heunghae basin, the geotechnical investigations were available for varying depths depending upon the bedrock depth. For a better estimation of liquefaction potential and lateral variation of liquefaction-prone soil along the basin, we considered the borehole data with maximum depth from the Integrated Data Base Center of National Geotechnical Information, Korea (https://www.geoinfo.or.kr), which varies from ground surface to 80 m (Figure 4). The typical bore Considering the clustered liquefaction features and closest available geotechnical datasets, we took a 10.5 km transect to determine the liquefaction potential in terms of thefactor of safety against liquefaction (FSL) of the soil. The variation of the depth of liquefiable soil, estimation of post liquefaction settlement and liquefaction potential index (I L ) along the southern part of the Heunghae basin experienced more severe liquefaction than the northern part of the basin. Since there was no detailed analysis of liquefaction potential of soils from the epicentral area, this study will help in determining the liquefiable soil depth and the possible post liquefaction settlement in the area, which is one of the most important factors in seismic-resistant building design.

Result of Selected SPT in the Heunghae Basin
Geotechnical investigation results from 41 boreholes covering the lateral extent of the Heunghae basin were collected. Around the Heunghae basin, the geotechnical investigations were available for varying depths depending upon the bedrock depth. For a better estimation of liquefaction potential and lateral variation of liquefaction-prone soil along the basin, we considered the borehole data with maximum depth from the Integrated Data Base Center of National Geotechnical Information, Korea (https://www.geoinfo.or.kr), which varies from ground surface to 80 m (Figure 4). The typical bore log and variation of SPT-N with depth for the borehole BB24, BB26, and BB 80 are presented in Figure 7. The geotechnical data suggest that the soils are mostly CH, ML, MH, SM, SP, SW, and GW as per the unified soil classification system (USCS). The SPT-N values for the boreholes taken for the analysis vary from 1 to 50, which indicates the soils are in a range of very loose to very dense conditions. The groundwater table varies from ground level to 15 m in the southern part of the basin.

SPT Based Liquefaction Potential Assessment
There are several methodologies for the estimation of liquefaction potential of the soil using in situ field testing parameters like SPT-N, Cone Penetration Test (CPT), shear wave velocity at 30 m (Vs30) as we all several laboratory testing methods such as cyclic triaxial test, resonant column test, shake table test, etc. [65][66][67]. The lab testing methods are complex and lead to difficulties associated with sample disturbance during sampling and testing. The empirical stress based SPT-N approach has gained a wide acceptance in geotechnical engineering communities [68]. Cyclic stress ratio (CSR)

SPT Based Liquefaction Potential Assessment
There are several methodologies for the estimation of liquefaction potential of the soil using in situ field testing parameters like SPT-N, Cone Penetration Test (CPT), shear wave velocity at 30 m (Vs 30 ) as we all several laboratory testing methods such as cyclic triaxial test, resonant column test, shake table test, etc. [65][66][67]. The lab testing methods are complex and lead to difficulties associated with sample disturbance during sampling and testing. The empirical stress based SPT-N approach has gained a wide acceptance in geotechnical engineering communities [68]. Cyclic stress ratio (CSR) is defined as the seismic demand of the soil, whereas the capacity of the soil to resist liquefaction is defined as cyclic resistance ratio (CRR) [6,[69][70][71]. These two parameters are used to determine the factor of safety against liquefaction (FSL).

Estimation of Cyclic Stress Ratio (CSR)
For the estimation of cyclic stress ratio (CSR) of soil, a simplified equation proposed by [69,70] has been used, which is given by Equation (1). The equation proposed by [68] was simplified by [70] to estimate the CSR values.
where a max is the peak horizontal acceleration on the ground surface, g is the acceleration due to gravity, σ v and σ' v is the total and effective stress and r d is the stress reduction factor of the soil. The stress reduction factor was calculated using an equation proposed by [71] and given in Equation (2).
where z is the depth in meters below the ground surface.

Estimation of Cyclic Resistance Ratio
Cyclic Resistance Ratio (CRR) has been used by several researchers for the liquefaction potential estimation. Several researchers have proposed equations for the estimation of CRR using SPT-N data [69,70,72,73]. In this present analysis, the equation proposed by [69] has been used. In the [70] method, the cyclic resistance ratio (CRR) needs fineness content correction and the (N 1 ) 60cs values are the value of (N 1 ) 60 to an equivalent clean sand value. The equation used for CRR estimation is given in Equation (3).
Since the earthquake magnitude is different from M w = 7.5, it is compulsory to modify the CRR value according to the magnitude of the earthquake used so that the CRR values will be more realistic for the earthquake used for the present liquefaction potential analysis. For the same, the magnitude was scaled down using the following equation

Correction for SPT-N Blow Count
SPT-N samplings were carried out by the Korean Industrial Standard (KSF) 2307 (test method specified for penetration test by Korean Industrial Standards) [63], which is similar to ASTM D1586-11 [63,64]. The measured SPT-N values were corrected using the NCEER report 1997 and 2001 [70]. The measured N values were corrected for overburden correction (C N ), hammer energy (C E ), Borehole diameter (C B ), rod length (C R ), liner correction (C S ) and fines content. The corrected values of ((N 1 ) 60cs ) was calculated using Equation (5).

Factor of Safety against Liquefaction (FSL)
The liquefaction potential defined as a factor of safety against liquefaction (FSL) and can be quantified as the ratio between cyclic resistance ratio (CRR) vs. cyclic stress ratio (CSR). The equation used for the estimation of factor of safety against liquefaction is given as Equation (7).
If the FSL value is <1, the soil is prone to liquefaction and if the FSL is >1, the soil is safe against liquefaction.

Assessment of Liquefaction Potential Index
The liquefaction potential index (I L ) quantifies the severity of liquefaction and predicts whether the liquefaction features will reach the ground surface or not [38,[73][74][75]. The liquefaction potential index for a specific site can be calculated by integrating the factor of safety (FSL) of a soil column with depth [38,76]. A weighting function has been added to give more weight to the shallow layer, close to the ground surface. In the present analysis, we estimate the liquefaction potential index (I L ) using Equation (8) given by [38] for all the boreholes taken into consideration for liquefaction analysis along the southern boundary of the Heunghae basin.
where z is the depth of the midpoint of the soil column taken for analysis, D z is the differential increment of depth. The weighting factor, w(z) and the severity factor, F(z) were calculated as: For the soil profile, I L was calculated using Equation (9) which is given below.
Where H i is the thickness of the individual soil layer. N is the number of the total soil layer, FSL is the factor of the safety against liquefaction, w i is the weighting factor (=10-0.5 Hz), and z i is the depth of i th layer (m) [76,77].

Post Liquefaction Settlement Estimation
For the level ground, the post liquefaction settlement can be estimated from the volumetric strain induced due to the rise in excess pore water pressure [42,[77][78][79][80]. Most of these methods are analytical and based on laboratory studies. In this present study, the methodology proposed by Ishihara and Yoshimine [80] was used for the post liquefaction settlement estimation.

Liquefaction Analysis
In total, 41 borehole pieces of data ( Figure 4) were taken from the epicentral area within 3 km radius from the epicenter. The selection of the borehole was designed in such a way that it covered the southern part of the Heunghae basin, which experienced more sand boils than the northern part of the basin. The soil ejected to the ground surface during the liquefaction was mostly sand and silty sand (Figures 2f and 6). The NOVOLIQ [81] software was used for the liquefaction analysis, which was designed for liquefaction potential analysis of layered soil. The input parameters for the liquefaction analysis were soil type, SPT N value, unit weight, percentage of fines, earthquake magnitude, groundwater table depth and distance from the fault, etc.

Results and Discussion
In total, 41 borehole pieces of data ( Figure 4) were taken for the present analysis. A cross-section covering 10 km along the E-W direction of the Heunghae basin was selected on the basis of surface manifestation of sand boils during the Pohang earthquake. Since there were no borehole data available passing thought the center of the basin, the southern part of the basin was taken for the detailed liquefaction analysis as a representative section of the basin. A recent study [47] reported that the southern part of the Heunghae basin had more reported cases of liquefaction than the northern part. An empirical method introduced by the National Center for Earthquake Engineering Research (NCEER) [5] was used for detailed liquefaction analysis to delineate the variation in depth of liquefaction prone soil around the epicentral area.
From the analysis, the silty sand and sandy soil having lower SPT-N values (between 1-20) show a factor of safety against liquefaction lower than 1. Along the 10.5 km transect, the central part of the Heunghae basin ( Figure 8) shows evidences of liquefaction at shallow depth. The silty soil, silty sand, and sandy soil from 1.5 to 15 m show FSL lower than 1 and are prone to liquefaction. The lateral variation of liquefiable soil layer for the southern part of the epicenter well matches with the present analysis ( Figure 8). The variation of FSL with depth for some of the boreholes with evidence of liquefaction is shown in Figure 9.
The liquefaction potential index (I L ) for each borehole was assessed using Equation (8). The estimated I L varies from 2 to 19 for the boreholes showing FSL lower than 1 and zero for boreholes having FSL higher than 1. The estimated I L for the boreholes are given in Table 3. In terms of liquefaction susceptibility, the central part of the basin is more susceptible than the east and west side of the southern part of the basin. This agrees well with the field evidence of liquefaction during the sand boil of the earthquake. The post liquefaction settlement estimated for each borehole varies from 6 to 50 cm and was mostly caused by the soil layers from 1.5 to 15 m, showing lower FSL. The post liquefaction settlement for each borehole is given in Table 3.   Understanding the spatial correlation between the liquefaction potential of soil, the probability of liquefaction, post-liquefaction settlement and its aerial extent with exposure to the population, buildings and infrastructures around the epicentral area will help for sustainable land use management in the future. Since this is the first time ever liquefaction phenomenon observed in the Korean Peninsula, detailed reconnaissance survey, land damage analysis, delineation of liquefiable subsurface soil layers will help in addressing the liquefaction risk around the epicentral area or areas with similar geological settings. On the basis of land damage distribution, most of the damaged lands are categorized as moderate to major land damage and confined between the Gokgang and Chogok Rivers. This may be due to the presence of alluvial soil at shallow depth. The land damage distribution map will help for the quantitative hazard assessment of the areas, especially for heavy construction. This information will also help for further detailed liquefaction hazard mapping with more geotechnical data and lead to liquefaction sustainable land use planning.

Conclusions
The paper presents the liquefaction analysis, probability of liquefaction and post liquefaction settlement for the Heunghae basin, which has experienced extensive liquefaction during the 2017 Pohang earthquake. A detailed post-liquefaction survey was carried out immediately after the earthquake. The result of liquefaction survey shows that the southern part has more sand boil than the northern part. To determine the depth of liquefiable soil along the southern part of the basin, SPT-N values were collected from the Integrated Data Base Center of National Geotechnical Information, Korea. From the present study, it can be inferred that: Understanding the spatial correlation between the liquefaction potential of soil, the probability of liquefaction, post-liquefaction settlement and its aerial extent with exposure to the population, buildings and infrastructures around the epicentral area will help for sustainable land use management in the future. Since this is the first time ever liquefaction phenomenon observed in the Korean Peninsula, detailed reconnaissance survey, land damage analysis, delineation of liquefiable subsurface soil layers will help in addressing the liquefaction risk around the epicentral area or areas with similar geological settings. On the basis of land damage distribution, most of the damaged lands are categorized as moderate to major land damage and confined between the Gokgang and Chogok Rivers. This may be due to the presence of alluvial soil at shallow depth. The land damage distribution map will help for the quantitative hazard assessment of the areas, especially for heavy construction. This information will also help for further detailed liquefaction hazard mapping with more geotechnical data and lead to liquefaction sustainable land use planning.

Conclusions
The paper presents the liquefaction analysis, probability of liquefaction and post liquefaction settlement for the Heunghae basin, which has experienced extensive liquefaction during the 2017 Pohang earthquake. A detailed post-liquefaction survey was carried out immediately after the earthquake. The result of liquefaction survey shows that the southern part has more sand boil than the northern part. To determine the depth of liquefiable soil along the southern part of the basin, SPT-N values were collected from the Integrated Data Base Center of National Geotechnical Information, Korea. From the present study, it can be inferred that: 1.
The liquefaction triggered during the M w 5.4 Pohang earthquake was localized and observed more in the southern part than the northern part of the basin. This may be due to the basin geometry, presence of the Heunghae fault at the center of the basin and the locality lying between the Gokgang River and the Chogok River with a shallow water table.

2.
On the basis of land damage distribution, most of the damaged lands are categorized as moderate to major land damage and confined between the Gokgang and Chogok Rivers.

3.
The grain size analysis of the ejected soil suggests that the soils are silty sand, sand with gravel, and sandy soil, which falls in the range of most liquefied soils reported in the literature. 4.
A 10.5 km long profile was chosen to determine the liquefaction potential, LPI, post-liquefaction settlement along the southern part of the basin, using SPT-N data and the Pohang earthquake ground motion data recorded at PHN2 station. The analysis indicates that the silty soil, silty sand, and sandy soil from 1.5 to 15 m deep are prone to liquefaction.

5.
Liquefaction susceptibility of the soil based on SPT-N value showed that the soils from 1.5 to 15 m possess FSL values from 0.2 to 1.0. The liquefaction potential index (I L ) estimated for the boreholes suggests that the southern part of the Heunghae basin has a low to high liquefaction potential index (I L ) values ranging from 2 to 19, which complements with the field evidence of sand boils during the earthquake. 6.
From the post liquefaction settlement, the soils having FSL values lower than 1 show maximum post-liquefaction settlement ranging from 6 to 50 cm.

7.
This study suggests that the soils of the Heunghae basin are prone to liquefaction for earthquakes higher than M w 5.4 and more comprehensive geotechnical investigations are required for better estimation of liquefiable soil depth along with the preparation of liquefaction hazard zonation maps. 8.
Similar liquefaction potential analysis must be carried out for other basins situated in close proximity to the Yangsan and Ulsan faults, including the Pohang basin, a major industrial center along the SE Korean Peninsula. This study will support the national disaster management authorities in developing a sustainable guideline for liquefaction hazard mitigation.