Improvement of Earthquake Risk Awareness and Seismic Literacy of Korean Citizens through Earthquake Vulnerability Map from the 2017 Pohang Earthquake, South Korea

: Earthquake activities in and around the Korean Peninsula are relatively low in number and intensity compared with neighboring countries such as Japan and China. However, recent seismic activity caused great alarm and concern among citizens and government authorities, and uncovered the level of preparedness toward earthquake disasters. A survey has been conducted on 1256 participants to investigate the seismic literacy of Korean citizens, including seismic knowledge, awareness and management using a questionnaire of citizen earthquake literacy (CEL). The results declared that the citizens had low awareness and literacy, which means that they are not properly prepared for earthquake hazards. To develop an earthquake risk reduction plan and program efﬁciently and effectively, not only must it appropriately characterize the target audience, but also indicate high potential earthquake zones and potential earthquake damage. Therefore, this study mapped and analyzed the seismic vulnerability in southeast Korea using LogitBoost, logistic model tree (LMT), and logistic regression (LR) machine learning algorithms based on a building damage inventory map. The damaged buildings’ locations were generated after the 2017 Pohang earthquake using the damage proxy map (DPM) method from the Sentinel-1 synthetic aperture radar (SAR) data. DPMs detected coherence loss, which indicates damaged buildings in urban areas in the Pohang earthquake and shows a good correlation with the Korea Meteorological Administration (KMA) report with modiﬁed Mercalli intensity (MMI) scale values of more than VII (seven). The damage locations were randomly divided into two datasets: 50% for training the vulnerability models and 50% for validating the models in terms of accuracy and reliability. Fifteen seismic-related factors were used to construct a model of each algorithm. Model validation based on the area under the receiver operating curve (AUC) was used to determine model accuracy. The AUC values of seismic vulnerability maps using the LogitBoost, LMT, and LR algorithms were 0.769, 0.851, and 0.749, respectively. We suggest that earthquake preparedness efforts should focus on reconstruction, retroﬁtting, renovation, and seismic education in areas with high seismic vulnerability in South Korea. The results of this study are expected to be beneﬁcial for engineers and policymakers aiming at developing disaster risk reduction plans, policies, and programs due to future seismic activity in South Korea. the citizens had low literacy, which means that they are not properly prepared for earthquake hazards. To develop an earthquake risk reduction plan and program efﬁciently and effectively, not only must one appropriately characterize the target audience, but also indicate high potential earthquake zones and potential earthquake damage. Therefore, this study mapped and analyzed the seismic vulnerability in southeast Korea using LogitBoost, logistic model tree (LMT), and logistic regression (LR) machine learning algorithms based on a building damage inventory map. The building damage locations were generated after the 2017 Pohang earthquake using the damage proxy map (DPM) method from the Sentinel-1 synthetic aperture radar (SAR) data. The DPMs manage to detect coherence loss, which indicates damaged buildings in residential and commercial areas due to the Pohang earthquake and show a good correlation with the Korea Meteorological Administration (KMA) report with modiﬁed Mercalli intensity (MMI) scale values of more than VII (seven). The damage locations were randomly divided into two datasets: 50% for training the vulnerability models and 50% for validating the models in terms of accuracy and reliability. Fifteen seismic-related factors were used to construct a model for each algorithm. Model validation based on the area under the receiver operating curve (AUC) was used to determine model accuracy. The AUC values of seismic vulnerability maps using the LogitBoost, LMT, and LR algorithms were 0.769, 0.851, and 0.749, respectively. We suggest that earthquake preparedness efforts should focus on reconstruction, retroﬁtting, renovation, and seismic education in areas with high seismic vulnerability in South Korea. The results of this study are expected to be beneﬁcial for engineers and


Introduction
Earthquake activities in and around the Korean Peninsula are relatively low in number and intensity compared with neighboring countries such as Japan and China, because it is located within the Eurasian intracontinental region [1]. However, seismographs often However, a small number of people in the earthquake area (close distance) were aware of the extreme fear and danger of earthquakes compared to the majority of people according to the survey results. Figure 2 shows that people who live near the epicenter of the Pohang earthquake have higher awareness of earthquakes because they were directly affected by the damage from the Pohang earthquake. Meanwhile, people who live far from the epicenter, such as people living in Seoul, have a lower awareness level of earthquakes because they only felt a slight shock from the Pohang earthquake, which was not fatal. Therefore, we need to inform many people of the dangers caused by such an earthquake, and when earthquakes occur in another area in the future, we need to be able to recognize the earthquake and respond to it. To reduce the effects of earthquake disasters and develop an earthquake risk reduction plan and program efficiently and effectively, not only is performing sustainable preparation, such as seismic literacy of citizens, necessary, but so is indicating high potential I think the place where I live is safe from earthquakes However, a small number of people in the earthquake area (close distance) were aware of the extreme fear and danger of earthquakes compared to the majority of people according to the survey results. Figure 2 shows that people who live near the epicenter of the Pohang earthquake have higher awareness of earthquakes because they were directly affected by the damage from the Pohang earthquake. Meanwhile, people who live far from the epicenter, such as people living in Seoul, have a lower awareness level of earthquakes because they only felt a slight shock from the Pohang earthquake, which was not fatal. Therefore, we need to inform many people of the dangers caused by such an earthquake, and when earthquakes occur in another area in the future, we need to be able to recognize the earthquake and respond to it. However, a small number of people in the earthquake area (close distance) were aware of the extreme fear and danger of earthquakes compared to the majority of people according to the survey results. Figure 2 shows that people who live near the epicenter of the Pohang earthquake have higher awareness of earthquakes because they were directly affected by the damage from the Pohang earthquake. Meanwhile, people who live far from the epicenter, such as people living in Seoul, have a lower awareness level of earthquakes because they only felt a slight shock from the Pohang earthquake, which was not fatal. Therefore, we need to inform many people of the dangers caused by such an earthquake, and when earthquakes occur in another area in the future, we need to be able to recognize the earthquake and respond to it. To reduce the effects of earthquake disasters and develop an earthquake risk reduction plan and program efficiently and effectively, not only is performing sustainable preparation, such as seismic literacy of citizens, necessary, but so is indicating high potential  To reduce the effects of earthquake disasters and develop an earthquake risk reduction plan and program efficiently and effectively, not only is performing sustainable preparation, such as seismic literacy of citizens, necessary, but so is indicating high potential earthquake zones and potential earthquake damage by producing seismic vulnerability maps. Seismic vulnerability assessment involves the comprehensive evaluation of factors that affect risks Remote Sens. 2021, 13,1365 4 of 25 associated with earthquakes within predefined areas. Urban areas are at higher risk of seismic disasters than outlying areas due to their higher building and infrastructure density and larger population. Therefore, in assessing seismic vulnerability, it is essential to select suitable influential factors and methods for the area of interest. Several methodologies have been applied for seismic vulnerability assessment and mapping during the past few decades [7,[12][13][14][15].
Seismic vulnerability assessment studies commonly analyze case studies using a combination of multicriteria decision making (MCDM) and geographic information system (GIS) approaches [12,13,16,17]. Among these, the analytical hierarchy process (AHP) is one of the most widely known MCDM methodologies; it stratifies and quantifies the importance of each applied influential factor to determine its relative importance, and assesses vulnerability by applying weights to all factors [12,16,18]. However, this method can be subjective because the opinion of the researcher can affect the weight assignment process; therefore, it is somewhat unsuitable for objective assessment. To address this problem, recent studies have applied hybrid models that combine various methodologies [14,19].
Many recent studies related to seismic vulnerability assessment and mapping have been conducted using machine learning techniques [12,16,[18][19][20][21]. For example, Han et al. (2019) [20] used a logistic regression (LR) model and applied the support vector machine (SVM) methodology to four kernel models (linear, polynomial, radial basis function, and sigmoid) to derive a suitable model for seismic vulnerability assessment; this study was notable in that the results of several seismic vulnerability models were compared analytically; such analyses are rarely conducted in this field, despite the broad application of machine learning techniques in recent years.
Providing training data plays an important role in the accuracy of the vulnerability map. Here, we used the damage proxy map (DPM) method to extract a building damage map for training and testing datasets. The DPM method is part of an ongoing collaborative effort between the Jet Propulsion Laboratory (JPL) and the California Institute of Technology, called the Advanced Rapid Imaging and Analysis (ARIA) project. The DPM method, using synthetic aperture radar (SAR) satellite data, has been shown to be useful for damage mapping following an earthquake and other natural disaster events, including the 2015 MW 7.8 Gorkha, Nepal earthquake using COSMO-SkyMed and ALOS-2 satellites [20], the 2019 typhoon Hagibis, Japan using Sentinel-1 satellites [21], the 2019 MW 7.1 Ridgecrest earthquake in California using Sentinel-1 satellites [22], and the 2014 eruption of Kelud volcano (Indonesia) using COSMO SkyMED satellites [23].
This study aims to improve the risk awareness and seismic literacy of Korean citizens through an earthquake vulnerability map of all buildings in southeast Korea. To produce the earthquake vulnerability map, we generated a damage proxy map (DPM) after the 2017 Pohang earthquake from the Sentinel-1 SAR dataset as a dependent variable, then applied machine learning to construct models using 15 seismic-related factors as independent variables. Model performances were verified using a receiver operating characteristic (ROC) curve. Finally, dangerous and safe areas were identified in southeast Korea by creating maps based on the model with the highest accuracy for each methodology, and the result were assessed. The results of this study should improve citizen earthquake risk awareness and seismic literacy, especially in high seismic vulnerability areas, and facilitate the construction of seismic vulnerability models that will be useful to reduce future losses due to earthquakes in South Korea.

Study Area
The study area covers three metropolitan cities; Busan, Daegu, and Ulsan, and two provinces; Geyongsangbuk and Gyeongsangnam, surrounding the epicenter of the Pohang earthquake. For simplicity, we refer to it as southeast Korea. The blue line in Figure 3 shows the border of the study area. In total, southeast Korea is home to 12,961,687 people and has an area of 32,285 km 2 [24]. Within the total area, urban areas account for 5.49%, followed by agriculture at 21.88%, forestry (66.71%), and other areas (5.92%). The proportions of males and females in these areas are 50.21% and 49.79%, respectively [25]. quake, an ML 5.8 earthquake shocked Gyeongju at 11:32:55 UTC on 12 September 2016. The epicenter was 8.7 km from the south of the city and 15 km beneath the surface [27]. The Gyeongju earthquake was just 40 km from the site of the Pohang earthquake. The Gyeongju earthquake was accompanied by 600 aftershocks, including an ML 5.1 foreshock that occurred near the mainshock at 10:44:32 UTC [28]. In 2018, 115 earthquakes with magnitudes of more than 2 occurred in the Korean Peninsula; among these, 36 earthquakes (31.13%) occurred in southeast Korea [29]. In 2019, 957 earthquakes with a magnitude of less than 2.0 occurred in the Korean Peninsula; among these, 294 earthquakes (30.72%) occurred in southeast Korea [30]. Among the 88 earthquakes of magnitude 2.0 or higher, 23 (26.17%) occurred in the same area.  The study area was chosen considering the high seismic activity in these areas. A series of aftershocks of the Pohang earthquake was observed to have occurred with a magnitude of 3.0 or more until 31 May 2018. A magnitude of 4.3 occurred near the epicenter, 2 h after the mainshock, and a magnitude 4.6 earthquake occurred 4 km to the southwest at around 20:03:04 UTC on 11 February 2018 [26]. One year prior to the Pohang earthquake, an M L 5.8 earthquake shocked Gyeongju at 11:32:55 UTC on 12 September 2016. The epicenter was 8.7 km from the south of the city and 15 km beneath the surface [27]. The Gyeongju earthquake was just 40 km from the site of the Pohang earthquake. The Gyeongju earthquake was accompanied by 600 aftershocks, including an M L 5.1 foreshock that occurred near the mainshock at 10:44:32 UTC [28]. In 2018, 115 earthquakes with magnitudes of more than 2 occurred in the Korean Peninsula; among these, 36 earthquakes (31.13%) occurred in southeast Korea [29]. In 2019, 957 earthquakes with a magnitude of less than 2.0 occurred in the Korean Peninsula; among these, 294 earthquakes (30.72%) occurred in southeast Korea [30]. Among the 88 earthquakes of magnitude 2.0 or higher, 23 (26.17%) occurred in the same area.
The Korean Peninsula lies at the eastern margin of the Eurasian Plate. About 30-15 million years ago, north-northeast (NNE)-striking strike-slip faults and NNE-to NE-striking normal faults settled predominantly in southeastern Korea and adjacent offshore areas when the East Sea opened in the early to middle Tertiary as a back-arc basin; smaller-scale coetaneous basins also formed, including the Pohang Basin. Although this region is about 400-500 km in length, its seismicity is affected by complex interactions of the Indo-Australian and Eurasian plates, as well as by subduction of Philippine Sea plates beneath the Japan and Ryuku trenches [31]. Several faults are distributed within the study area, including Dongrae, Moryang, Miryang, Ulsan, Wangsan, and Yangsan [32]. Due to seismic history and geographic characteristics, the probability of earthquake occurrence in southeast Korea is considered relatively high, and secondary damage in the event of an earthquake with a medium or higher magnitude constitutes an unusually high risk. Therefore, sustainable preparation and management planning for such events is required.

SAR Datasets
A building damage inventory map for producing the seismic vulnerability map in southeast Korea was generated using Sentinel-1 synthetic aperture radar (SAR) C-band data (5.5 cm wavelength) provided by the European Space Agency (ESA). Pohang is located in between two frames (470 and 475), path 61 of the Sentinel-1 imagery; therefore, we obtained six Sentinel-1 single look complex (SLC) images with vertical transmission and vertical return (VV) polarization for the Pohang earthquake, with four scenes prior to the event on 23 October 2017, and 4 November 2017, and two scenes after the event on 16 November 2017. Furthermore, we needed to merge the scenes before processing. The images were co-registered, with the 4 November 2017 scenes as a reference.

Damage Proxy Map (DPM)
DPMs generated from the comparison of pre-and co-event SAR images can help identify damage caused by earthquakes using remote sensing imagery [20]. The method relies on the reduction in the coherence of the radar echoes between satellite-based SAR images taken before and after the earthquake to identify anomalous changes in ground surface properties. Coherence measures the change in radar backscatter from the ground, a proxy for the ground-surface property changes. A low coherence implies a large change to the ground surface that reflected the SAR radiation [33]. Changes can be caused by damage to the ground itself or damage to structures.
The process started with image co-registration, with the 4 November 2017 scenes as a reference. This co-registration process was done with sub-pixel accuracy to match scenes to one another. We used the complex pixel value, c, of the pre-processed SLCs for the change detection analysis, where damage is inferred from loss of coherence or decorrelation between SAR images [20]. We computed the pre-and co-event interferometric coherences, γ (Equation (1)), from a pair of SLCs before the event and another pair spanning across the event, respectively [21].
where c 1 and c 2 are complex pixel values of two co-registered SAR images and * denotes the complex conjugate. The resulting coherence ranges from 0 (incoherent) to 1 (coherent). The coherence is equal to 1 if the observation is identical in the two images because of the stable object-like buildings in the scene. The pre-event coherence represented change unrelated to the event and was assumed to be the background value. Then, we obtained a coherence difference (COD) by subtracting γ pre-event from γ co-event . Therefore, the process could generate a COD ranging from −1 to 1. A negative COD (or coherence gain) usually indicates surface changes occurring between the pre-event scenes and is associated with changes not related to the event. Coherence gain could happen in an agricultural area when the fields are full of crops. Then, when harvested during the period time of the pre-event interferometric pair, the area has low coherence. After harvesting, leaving an empty field, the area has a greater coherence in the co-event interferometric pair (γ co-event > γ pre-event ), so a negative COD is obtained. A positive COD (or coherence loss) indicates surface changes between the co-event scenes spanning the events, such as major damage to a building significantly that increases the interferometric phase variance, causing a decrease in coherence. Hence, the loss of coherence is most effective for detecting damage in built-up areas caused by earthquake. However, COD is generally less effective and less reliable in vegetated areas where coherence changes may be random. Therefore, this study focused on a DPM in urban, built-up areas as changes can be detected easily in SAR imagery. A greater loss in coherence generally correlates with greater severity of the change, such as a fully collapsed building, causing more significant coherence loss than partial collapse [21]. The threshold for significant coherence loss can be chosen by comparing observed coherence changes with reported damage and areas in which it is known that no damage occurred. Yun et al. (2015) compared a DPM with a National Geospatial-Intelligence Agency (NGA) analysis and the United Nations Operational Satellite Application Programme (UNOSAT) damage assessment map. Tay et al. (2020) used high-resolution aerial imagery from the Geospatial Information Authority of Japan (GSI). Here, coherence loss thresholds for DPMs were chosen by considering the damaged area from the Korea Meteorological Administration (KMA) report of the Pohang earthquake [34].

Selection of Seismic-Related Factors
We selected factors affecting seismic vulnerability based on the result of previous studies [7,14,15,17]. The factors affecting the seismic vulnerability were prepared based on five main indicators, which were geotechnical, physical, structural, social, and capacity; we selected a total of 15 factors corresponding to these categories. Geotechnical factors included slope and altitude; physical factors included peak ground acceleration (PGA), epicenter distance, and fault distance; structural factors included land use, construction materials, building density, and building height; social factors included elderly population (≥65 years), child population (<15 years), and population density; and capacity factors included distances from hospitals, fire stations, and police stations. The factors were organized into raster-based spatial databases (30 m spatial resolution) and were reclassified using the quantile method to identify and analyze the effect of each class. The data used in this present study are shown in Figure 4.
in coherence. Hence, the loss of coherence is most effective for detecting damage in builtup areas caused by earthquake. However, COD is generally less effective and less reliable in vegetated areas where coherence changes may be random. Therefore, this study focused on a DPM in urban, built-up areas as changes can be detected easily in SAR imagery. A greater loss in coherence generally correlates with greater severity of the change, such as a fully collapsed building, causing more significant coherence loss than partial collapse [21].
The threshold for significant coherence loss can be chosen by comparing observed coherence changes with reported damage and areas in which it is known that no damage occurred. Yun et al. (2015) compared a DPM with a National Geospatial-Intelligence Agency (NGA) analysis and the United Nations Operational Satellite Application Programme (UNOSAT) damage assessment map. Tay et al. (2020) used high-resolution aerial imagery from the Geospatial Information Authority of Japan (GSI). Here, coherence loss thresholds for DPMs were chosen by considering the damaged area from the Korea Meteorological Administration (KMA) report of the Pohang earthquake [34].

Selection of Seismic-Related Factors
We selected factors affecting seismic vulnerability based on the result of previous studies [7,14,15,17]. The factors affecting the seismic vulnerability were prepared based on five main indicators, which were geotechnical, physical, structural, social, and capacity; we selected a total of 15 factors corresponding to these categories. Geotechnical factors included slope and altitude; physical factors included peak ground acceleration (PGA), epicenter distance, and fault distance; structural factors included land use, construction materials, building density, and building height; social factors included elderly population (≥65 years), child population (<15 years), and population density; and capacity factors included distances from hospitals, fire stations, and police stations. The factors were organized into raster-based spatial databases (30 m spatial resolution) and were reclassified using the quantile method to identify and analyze the effect of each class. The data used in this present study are shown in Figure 4.  Slope and elevation were extracted from a digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) using basic terrain analysis tools. Slope and elevation are factors affecting the vulnerability of urban environments to earthquakes [14,35]. Degradation in terrain with steep topography, especially at the top of hills and peaks, is greatly enhanced. According to construction standards, a slope of 5 to 9% is suitable for urbanization [16]. The highs and lows of elevation in each area are highly correlated with landslide susceptibility in each area [36]. Therefore, because of the amount of erosion and its relation to human activity, the higher the altitude of an area, the greater the seismic vulnerability. Figure 4a,b shows the slope and elevation maps of southeast Korea, respectively.
Peak ground acceleration (PGA) (Figure 4c) is the degree to which the ground shakes at the Earth's surface and is related to the amount of fault activity [18]. In this study, raw data from the Korea Institute of Geoscience and Mineral Resources (KIGAM) were converted to acceleration data and interpolated throughout southeast Korea [27,37]. The epicenter, or location where an earthquake occurs, is an important factor related to earthquake occurrence; the level of damage is different depending on the ground condition or the structure of the fault plane, on which the greatest damage often occurs at the epicenter. Therefore, we used distance data from earthquake epicenters ( Figure 4d) with a magnitude over 4 from 1978 to 2020, including the 2016 Gyeongju earthquake and the 2017 Pohang earthquake. The epicenter locations were acquired from the United States Geological Survey (USGS). The faults are forms of tectonic factors whose presence or absence can be examined in relation to the seismic hazard of different areas. Fault distance ( Figure 4e) plays a key role in vulnerability to earthquake hazards, as proximity to the fault causes high seismic risk and damage, and distance from it will reduce the risk and consequently provide higher resilience [16].
Anti-seismic design in South Korea was introduced in 1988, and it is mandatory only for buildings that are three stories or higher [38]. As of November 2016, 29.9% of residential buildings and 23.7% of non-residential buildings in Seoul were designed to be antiseismic. As there is no guarantee that future earthquakes will not exceed the magnitude of the Gyeongju earthquake, most buildings in South Korea are considered highly vulnerable. To assess their vulnerability, we identified four structural factors of seismic vulnerability: land use, construction materials, building density, and building height, depicted Slope and elevation were extracted from a digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) using basic terrain analysis tools. Slope and elevation are factors affecting the vulnerability of urban environments to earthquakes [14,35]. Degradation in terrain with steep topography, especially at the top of hills and peaks, is greatly enhanced. According to construction standards, a slope of 5 to 9% is suitable for urbanization [16]. The highs and lows of elevation in each area are highly correlated with landslide susceptibility in each area [36]. Therefore, because of the amount of erosion and its relation to human activity, the higher the altitude of an area, the greater the seismic vulnerability. Figure 4a,b shows the slope and elevation maps of southeast Korea, respectively.
Peak ground acceleration (PGA) (Figure 4c) is the degree to which the ground shakes at the Earth's surface and is related to the amount of fault activity [18]. In this study, raw data from the Korea Institute of Geoscience and Mineral Resources (KIGAM) were converted to acceleration data and interpolated throughout southeast Korea [27,37]. The epicenter, or location where an earthquake occurs, is an important factor related to earthquake occurrence; the level of damage is different depending on the ground condition or the structure of the fault plane, on which the greatest damage often occurs at the epicenter. Therefore, we used distance data from earthquake epicenters (Figure 4d) with a magnitude over 4 from 1978 to 2020, including the 2016 Gyeongju earthquake and the 2017 Pohang earthquake. The epicenter locations were acquired from the United States Geological Survey (USGS). The faults are forms of tectonic factors whose presence or absence can be examined in relation to the seismic hazard of different areas. Fault distance (Figure 4e) plays a key role in vulnerability to earthquake hazards, as proximity to the fault causes high seismic risk and damage, and distance from it will reduce the risk and consequently provide higher resilience [16].
Anti-seismic design in South Korea was introduced in 1988, and it is mandatory only for buildings that are three stories or higher [38]. As of November 2016, 29.9% of residential buildings and 23.7% of non-residential buildings in Seoul were designed to be anti-seismic. As there is no guarantee that future earthquakes will not exceed the magnitude of the Gyeongju earthquake, most buildings in South Korea are considered highly vulnerable. To assess their vulnerability, we identified four structural factors of seismic vulnerability: land use, construction materials, building density, and building height, depicted in Figure 4f-i, respectively [7,17]. The greater the number of floors of a building, despite its quality, the greater the vulnerability. The number of floors in the building, if not in accordance with safety principles, will definitely increase the damage [39]. Even if the height is treated with due diligence and calculations, it is difficult for the evacuation of buildings, and to search for and rescue people. Obviously, structures of high strength and standard materials have good earthquake safety [40]. Proper deployment of land uses on the basis of urban planning principles, such as proper accessibility, proper distance from the biological hotspots, safety, comfort, and utility can substantially reduce the amount of vulnerability, injury, and economic damage [41].
Increasing population growth, population density, and poor distribution of services and infrastructure pose risks to society [42]. In recent earthquakes around the world, it can be said that most of the damage is to humans and with the increase in population, it is predicted that in the future, the mortality rate will be higher. The earthquake hazard coefficient in urban centers is also more complex and riskier due to urbanization without planning and development [43]. In events such as earthquakes, everyone in the community is vulnerable, but older people and children are the most vulnerable groups in a community and more attention is needed to minimize pain and injury [44]. Children do not tolerate disruption well and older people are psychologically fragile because of their disrupted life rhythms. The elderly population, child population, and population density in southeast Korea are presented in Figure 4j-l, respectively.
We identified the locations of social infrastructure facilities that offer aid in the event of an earthquake, and of hazardous facilities that have the potential to cause huge damage. The degree of accessibility following a disaster was analyzed by considering the physical distances to three factors, including social infrastructure facilities (hospital, police station, and fire station). Distance from a hospital ( Figure 4m) and access to health services such as hospitals play a key role in controlling post-emergency complications and providing earthquake rescue and hospitalization services. Proper and quick access to medical facilities will increase earthquake resilience [45]. Distance from a fire station (Figure 4n) and access to police stations (Figure 4o) through the communication networks will speed up rescue operations and service the injured. As such, the greater the distance from fire stations and police stations, the greater the vulnerability [46].

Machine Learning
To map seismic vulnerability using a machine learning algorithm, several steps must be performed. First, the spatial relationships between the damaged buildings from the DPM and related factors (geotechnical, physical, structural, social, and capacity indicators) are calculated using the frequency ratio (FR) method. In this research, we analyzed the spatial relationship between damaged buildings' locations (1623 cells) and 15 factors related to seismic vulnerability, based on the FR value of each factor. When the ratio is greater than 1, this denotes that the class in each factor has a closer relationship with seismic vulnerability [47]. In the case that each factor has a less close relationship with seismic vulnerability, the ratio is less than 1. The FR for each factor was calculated by Equation (2) [48]. FR = % of class of related factor % of total area To apply the machine learning algorithm, we used the 1623 cells of damaged buildings generated from the DPM method. Among these cells, 50% (812) were used as a training dataset and 50% (811) were used as a test dataset. We extracted the same number of cells corresponding to undamaged buildings. All cells were randomly sampled and generating models and the accuracy of each model was done based on training (1624) and test datasets (1622). Several seismic-related maps, including geological maps, were produced at a 30 m resolution. Then, all data were classified as categorical or continuous. Continuous variables include the slope, elevation, PGA, distance from epicenters, distance from faults, building density, building height, child population, elderly population, population density, distance from hospitals, distance from police stations, and distance from fire stations. Categorical variables include construction materials and land use.
Model validation was carried out through ROC curve analysis of the testing dataset (50%). Receiver operating characteristic (ROC) curve analysis, as an index of model performance, is commonly used to assess predictive accuracy [49]. To quantitatively determine the accuracy of the model verification, the area under the curve (AUC) of the ROC curve is calculated for the total area and correct predictive accuracy is obtained. AUC values range between 0.5 and 1; higher values indicate more reliable algorithm performance. The workflow of the seismic vulnerability mapping carried out in this study is provided in Figure 5. faults, building density, building height, child population, elderly population, population density, distance from hospitals, distance from police stations, and distance from fire stations. Categorical variables include construction materials and land use.
Model validation was carried out through ROC curve analysis of the testing dataset (50%). Receiver operating characteristic (ROC) curve analysis, as an index of model performance, is commonly used to assess predictive accuracy [49]. To quantitatively determine the accuracy of the model verification, the area under the curve (AUC) of the ROC curve is calculated for the total area and correct predictive accuracy is obtained. AUC values range between 0.5 and 1; higher values indicate more reliable algorithm performance. The workflow of the seismic vulnerability mapping carried out in this study is provided in Figure 5.

LogitBoost
LogitBoost is a boosting algorithm developed by Friedman et al. [50] to reduce bias and variance. The LogitBoost algorithm was modified from AdaBoost, which was the commonly used boosting method for handling noisy data that executes additive logistic regression with least-square fits for individual classes. LogitBoost reduces training errors and enhances classification accuracy by using additive logistic regression for classification with a base-learning regression scheme and an ability to perform multiclass classification. The damaged building inventory map was divided into two classes: damaged buildings and undamaged buildings, using Equation (3): where D is the number of building damage-dependent factors and βi is the coefficient of the i-th component within input vector x. Probabilities were constructed using the linear logistic regression method with Equation (4):

LogitBoost
LogitBoost is a boosting algorithm developed by Friedman et al. [50] to reduce bias and variance. The LogitBoost algorithm was modified from AdaBoost, which was the commonly used boosting method for handling noisy data that executes additive logistic regression with least-square fits for individual classes. LogitBoost reduces training errors and enhances classification accuracy by using additive logistic regression for classification with a base-learning regression scheme and an ability to perform multiclass classification. The damaged building inventory map was divided into two classes: damaged buildings and undamaged buildings, using Equation (3): where D is the number of building damage-dependent factors and β i is the coefficient of the i-th component within input vector x. Probabilities were constructed using the linear logistic regression method with Equation (4): where C is the number of classes and the least-square fit Lc(x) is resolved such that ∑ C = 1 L C C (x) = 0 to set up the least number of instances per node of the logistic model trees.

Logistic Model Tree (LMT)
The logistic model tree combines the C4.5 algorithm [51] and logistic regression (LR) functions. The information gain ratio (IGR) technique is applied to split the tree into nodes and leaves, and the LogitBoost algorithm [52] is used to fit the logistic regression functions at a tree node. The C4.5 algorithm uses the entropy technique for feature selection because it is the fastest method for providing reliable classification accuracy [53]. The over-fitting problem, which is an important challenge in LMT modeling, is overcome using the CART algorithm, which prunes the tree for modeling the training dataset [54]. The IGR can be formulated using Equation (5): where gain (A) is the information after attribute A is selected as a test for classification of the training samples and split info (A) is the information generated when x training samples are categorized into n subsets [51]. In the next step, the LogitBoost algorithm performs additive logistic regression with least-squares fit for each class Ci (damaged or undamaged building) according to Equation (6) [55]: where L c (x) is the least-squares fit and CF and α i are, respectively, the number of seismicrelated factors and the coefficient of the i-th element of vector x. The a posteriori probabilities in the leaves of the LMT are calculated using the linear logistic regression model with Equation (7) [52]: where c is the number of building damage classes and Lc(x), the least-squares fit, is transformed in such a way that ∑ c c =1 L c (x) = 0.

Logistic Regression (LR)
The logistic regression (LR) model, developed by McFadden (1973) [56], is a multivariate regression analysis model that describes the relationship between a bivariate dependent parameter and several independent parameters [57] through the estimation of an optimal model. The addition of a link function suitable for a general linear regression model allows the parameter type to be continuous, discrete, or mixed, thus obviating the requirement for a normal distribution [58,59]. Some studies have shown that the LR model is more accurate than other types of models constructed for the same purpose [60][61][62]. The LR model based on a general linear model can be derived from Equations (8) and (9): P = e y 1 + e y (9) where y is the linear logistic model, b 0 is the y-intercept, b n is the logistic coefficient of each factor, n is the number of factors controlling a seismic event, x is the earthquake conditioning factor, and P is the probability of damage (ranging from 0 to 1) in the event of an earthquake [60].

Building Damage Inventory Map
The DPMs were generated from the Sentinel-1 dataset and geocoded to the Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM), 1 arcsecond. Figure 6 reveals the map view of Sentinel-1 DPMs over Pohang. The assessment technique is most sensitive to the destruction of the built environment. Pixels are set to be relatively transparent where corresponding to areas where decorrelation did not significantly change during the time spanning the earthquake, suggesting little to no destruction. Increased opacity of the radar image pixels reflects increasing ground and building change or potential damage. The color range from yellow to red indicates an increasingly significant coherence change in the area covered by the pixel. Each pixel in the DPMs was registered to the SRTM DEM and had a corresponding dimension of about 30 m.
parent where corresponding to areas where decorrelation did not significantly change during the time spanning the earthquake, suggesting little to no destruction. Increased opacity of the radar image pixels reflects increasing ground and building change or potential damage. The color range from yellow to red indicates an increasingly significant coherence change in the area covered by the pixel. Each pixel in the DPMs was registered to the SRTM DEM and had a corresponding dimension of about 30 m.
Here, we compared the distribution of coherence loss areas after the earthquakes throughout Gyeongju and Pohang, as seen in Figure 6a. According to the land cover map derived from the Korea Institute of Geoscience and Mineral Resources (KIGAM), Figure  6a reveals the DPM after the Gyeongju earthquake over residential, commercial, agricultural, and vegetated areas. Figure 6b,c show the widespread COD in Seonggeon-dong and Gwangmyeong-dong, respectively, corresponding to collapsed houses (Figure 6d [63]),. These areas consist of residential and commercial areas. Figure 6a shows DPM that indicates COD over residential, commercial, industrial, agricultural, and vegetated areas affected by the Pohang earthquake. Some of the spatially large and strong coherence loss after the Pohang earthquake corresponds to Handong Global University (Figure 6e), associated with collapsed walls (Figure 6f), while Figure 6g shows Songdo-dong district, consisting of residential areas. The DPM of the Pohang earthquake was then used as a building damage inventory map to produce a seismic vulnerability map using machine learning.  Here, we compared the distribution of coherence loss areas after the earthquakes throughout Gyeongju and Pohang, as seen in Figure 6a. According to the land cover map derived from the Korea Institute of Geoscience and Mineral Resources (KIGAM), Figure 6a reveals the DPM after the Gyeongju earthquake over residential, commercial, agricultural, and vegetated areas. Figure 6b,c show the widespread COD in Seonggeon-dong and Gwangmyeong-dong, respectively, corresponding to collapsed houses (Figure 6d [63]). These areas consist of residential and commercial areas. Figure 6a shows DPM that indicates COD over residential, commercial, industrial, agricultural, and vegetated areas affected by the Pohang earthquake. Some of the spatially large and strong coherence loss after the Pohang earthquake corresponds to Handong Global University (Figure 6e), associated with collapsed walls (Figure 6f), while Figure 6g shows Songdo-dong district, consisting of residential areas. The DPM of the Pohang earthquake was then used as a building damage inventory map to produce a seismic vulnerability map using machine learning.

Relationship between Damaged Buildings and Related Factors
FR values can provide information on the relationship between seismic vulnerability and related factors (geotechnical, physical, structural, social, and capacity). A ratio greater than 1 denotes that the class in the related factor has more impact on seismic vulnerability. The FR values calculated in this study are shown in Table 1 (4.18). Areas in these classes are predicted to experience the highest degree of damage due to earthquakes.

Seismic Vulnerability Map
Seismic vulnerability maps were made using the training dataset compiled using the building damage inventory map from the DPM of the Pohang earthquake and applying machine learning algorithms, as discussed above. A combination of 15 seismic-related factors served as the dependent variables, and can mainly be classified as geotechnical, physical, structural, social, and capacity indicators. LogitBoost (Figure 7a), LMT (Figure 7b), and LR (Figure 7b) machine learning algorithms were applied to produce the seismic vulnerability maps. Each pixel in the study area was assigned a specific building damage value using the natural breaks method [7]. The seismic vulnerability maps were classified as safe, low to moderate, high, and very high vulnerability classes.  Figure 8 shows the distribution of pixels in each seismic vulnerability map generated by LogitBoost, LMT, and LR models. In the LogitBoost model, 29.53% were classified as safe, 17.59% as low risk, 13.18% as moderate risk, 22.95% as high risk, and 16.74% as very high risk. Gyeongju and Pohang were found to be the most vulnerable to earthquake damage. Among big cities in the study area, Busan, Daegu, and Ulsan were classified as safe and low risk. For the LMT model, 19.74% were classified as safe, 18.54% as low risk, 8.25% as moderate risk, 18.27% as high risk, 35.20% as very high risk. Gyeongju and Pohang  Figure 8 shows the distribution of pixels in each seismic vulnerability map generated by LogitBoost, LMT, and LR models. In the LogitBoost model, 29.53% were classified as safe, 17.59% as low risk, 13.18% as moderate risk, 22.95% as high risk, and 16.74% as very high risk. Gyeongju and Pohang were found to be the most vulnerable to earthquake damage. Among big cities in the study area, Busan, Daegu, and Ulsan were classified as safe and low risk. For the LMT model, 19.74% were classified as safe, 18.54% as low risk, 8.25% as moderate risk, 18.27% as high risk, 35.20% as very high risk. Gyeongju and Pohang were found to be the most vulnerable to earthquake damage. Daegu was classified as a safe and low-risk city, while Busan and Ulsan were high and very high risk. In the LR model, 19.91% were classified as safe, 20.08% as low risk, 20.71% as moderate risk, 19.46% as high risk, and 19.85% as very high risk. The distribution of pixels in the low-and high-risk classes in Figure 8 shows similar results for each algorithm, although LMT shows larger numbers of pixels in very high-risk areas. The most vulnerable areas were Gyeongju and Pohang, whereas low and moderate risk areas were Busan and Ulsan. Seismic vulnerability classes were evenly distributed in Daegu.

Model Validation
A validation step was conducted to assess the reliability of the seismic vulnerability map from each algorithm. ROC curve analysis is a standard way of validating the probability models used to generate seismic vulnerability maps, according to the area under the curve (AUC) [7,14,15]. Higher values indicate more accurate and reliable models. If the AUC, which ranges from 0 to 1, is lower than 0.5, the model is considered unacceptably inaccurate [49]. The accuracy of the seismic vulnerability maps generated using the three algorithms was then evaluated based on ROC curve analysis of the testing dataset (50% of all data). As seen in Figure 9, the AUC values were 0.769, 0.851, and 0.749 for LogitBoost, LMT, and logistic regression, respectively. Thus, the LMT model generated the best seismic vulnerability map in this study. The results indicate that the algorithms are useful to map seismic vulnerability in the southeastern Korean Peninsula. Since all of the AUC values were higher than 0.5, the seismic vulnerability maps produced by all algorithms used in this study are acceptable for predicting vulnerable buildings in southeast Korea [49].

Model Validation
A validation step was conducted to assess the reliability of the seismic vulnerability map from each algorithm. ROC curve analysis is a standard way of validating the probability models used to generate seismic vulnerability maps, according to the area under the curve (AUC) [7,14,15]. Higher values indicate more accurate and reliable models. If the AUC, which ranges from 0 to 1, is lower than 0.5, the model is considered unacceptably inaccurate [49]. The accuracy of the seismic vulnerability maps generated using the three algorithms was then evaluated based on ROC curve analysis of the testing dataset (50% of all data). As seen in Figure 9, the AUC values were 0.769, 0.851, and 0.749 for LogitBoost, LMT, and logistic regression, respectively. Thus, the LMT model generated the best seismic vulnerability map in this study. The results indicate that the algorithms are useful to map seismic vulnerability in the southeastern Korean Peninsula. Since all of the AUC values were higher than 0.5, the seismic vulnerability maps produced by all algorithms used in this study are acceptable for predicting vulnerable buildings in southeast Korea [49].
Remote Sens. 2021, 13, x FOR PEER REVIEW Figure 9. Receiver operating characteristic (ROC) curves associated with seismic vulnerabil maps generated by LogitBoost, LMT, and LR algorithms.

Building Damage Inventory Map
The building damage inventory map was generated through the DPM metho the Pohang earthquake using Sentinel-1 SAR data. We found a good correlation b the DPM result and the released map from the KMA report about the Pohang earth in our qualitative validation. Therefore, the building damage inventory map is r The KMA map was derived from field survey and damage survey data from loc ernments. The map shows the distribution of the Pohang earthquake's magnitud the modified Mercalli intensity (MMI) scale [34]. The results were similar to the report that these areas suffered MMI VII to VIII. KMA defines MMI VII to VIII as d to major structural parts such as pillars, walls, and roofs, even in well-designed an built buildings. KMA map also shows some areas that suffered MMI V to VI; ho the DPMs did not show any CODs in these areas. The scales showed that the dama inside buildings, such as minor cracks in walls and damage caused by dropping or tiles; therefore, SAR cannot detect a significant coherence change.
Here, the total damaged areas from DPMs were calculated by multiplying th number of pixels by the pixel cell size, which for the Gyeongju earthquake yielded of 1.09 km 2 and the Pohang earthquake yielded an area of 1.32 km 2 . The KMA rep total of 69 damaged buildings by the Gyeongju earthquake, and 504 damaged bu by the Pohang earthquake (MMI VII and VIII). The Pohang earthquake caused mor age than the Gyeongju earthquake, some of which was due to the depth of the epi The shallower the epicenter of an earthquake, the more damage it causes. The ep depth of the Pohang earthquake was 7 km, while the epicenter depth of the Gy earthquake was 15 km. Additionally, the surface deformation in urban areas wit buildings caused by the Pohang earthquake affected the occurrence of damage in P

Building Damage Inventory Map
The building damage inventory map was generated through the DPM method after the Pohang earthquake using Sentinel-1 SAR data. We found a good correlation between the DPM result and the released map from the KMA report about the Pohang earthquake in our qualitative validation. Therefore, the building damage inventory map is reliable. The KMA map was derived from field survey and damage survey data from local governments. The map shows the distribution of the Pohang earthquake's magnitude using the modified Mercalli intensity (MMI) scale [34]. The results were similar to the KMA report that these areas suffered MMI VII to VIII. KMA defines MMI VII to VIII as damage to major structural parts such as pillars, walls, and roofs, even in well-designed and well-built buildings. KMA map also shows some areas that suffered MMI V to VI; however, the DPMs did not show any CODs in these areas. The scales showed that the damage was inside buildings, such as minor cracks in walls and damage caused by dropping objects or tiles; therefore, SAR cannot detect a significant coherence change.
Here, the total damaged areas from DPMs were calculated by multiplying the total number of pixels by the pixel cell size, which for the Gyeongju earthquake yielded an area of 1.09 km 2 and the Pohang earthquake yielded an area of 1.32 km 2 . The KMA reported a total of 69 damaged buildings by the Gyeongju earthquake, and 504 damaged buildings by the Pohang earthquake (MMI VII and VIII). The Pohang earthquake caused more damage than the Gyeongju earthquake, some of which was due to the depth of the epicenter. The shallower the epicenter of an earthquake, the more damage it causes. The epicenter depth of the Pohang earthquake was 7 km, while the epicenter depth of the Gyeongju earthquake was 15 km. Additionally, the surface deformation in urban areas with high buildings caused by the Pohang earthquake affected the occurrence of damage in Pohang.

Seismic Vulnerability Map
Estimating the seismic vulnerability of an area is vital for environmental management and land use planning, among other applications [19]. Although many methods and techniques have been developed to assess earthquake hazards around the world to date, the goals of all these studies are to reduce the economic losses and resulting losses. The method used to create seismic vulnerability maps affects the quality of the mapping. Machine learning techniques are effective. In particular, the method used to generate training and testing data is important. Accurate damage building inventory maps can be obtained using the DPM method; we combined DPM and GIS spatial data to produce an accurate seismic vulnerability map.
The seismic vulnerability maps of all three algorithms used in this study revealed that Pohang and Gyeongju are the most vulnerable areas to earthquake damage. We analyzed seismic-related factors by comparing general patterns of damaged buildings with factor maps. Pohang and Gyeongju are the cities where the two largest recent earthquakes happened. Several factors in these areas had higher values, such as the value of PGA and the distance from the epicenter. Therefore, Gyeongju and Pohang were areas with a high risk of earthquake vulnerability. Buildings within 0-36.90 km of an epicenter corresponded to damaged buildings. This result confirms that most buildings close to epicenters were damaged. The results revealed that areas with a high risk had high population and building density, including Gyeongju and Pohang. Similar to Gyeongju and Pohang, Busan, Daegu, and Ulsan are cities with a high density of buildings and populations. These cities consist of areas of low to moderate risk of seismic vulnerability. One of the main causes of a high risk of seismic damage in countryside areas is wood as a building construction material, while buildings constructed with steel reinforced with concrete have a lower level of seismic vulnerability. Therefore, attention should be paid to reconstruction, retrofitting, or renovation of the buildings in these areas.
Seismic vulnerability maps were validated based on ROC curves and AUC values to assess the accuracy of the maps using testing data, which comprised 50% of the total dataset. The AUC data showed that the LMT algorithm had the highest accuracy of 85.10%, which was 8.2% higher than the LogitBoost algorithm and 10.2% higher than the logistic regression algorithm. Therefore, the LMT algorithm is better to produce seismic vulnerability maps than LogitBoost and logistic regression algorithms.
A survey was conducted with 1256 Korean citizens using a questionnaire of citizen earthquake literacy (CEL) during spring 2020. The survey was based on three dimensions, including citizen knowledge, awareness, and management. We developed 15 questions associated with three dimensions. Table A1 in Appendix A presents the questionnaire that consists of five questions for each dimension, and all questions are configured to be answered on a 5-point Likert scale (strongly agree, agree, neutral, disagree, and strongly disagree). Principal component analysis as an exploratory factor analysis approach was used to examine if the items corresponded to each other and explore the construct validity of the instrument. Cronbach's alpha when an item was deleted was also calculated to see whether the items in the construct reliably measured the same latent variable. Cronbach's alpha values for each dimension were 0.847, 0.822, and 0.849, respectively. Furthermore, Cronbach's alpha when an item was deleted in this study was found to be between 0.77 and 0.84. This result indicated that the survey data used in this study were strongly reliable and consistent [64].
To characterize the profile of participants with higher (or lower) levels of seismic literacy, difference in means analyses using a t-test, ANOVA test, and Tukey's honestly significant difference (HSD) post hoc test were carried out. Table 2 shows the average values associated with earthquake literacy, including seismic knowledge, awareness, and management, broken down by sociodemographic characteristics of the participants. The t-test result shows that male citizens have higher earthquake literacy but not significantly more than female citizens in all three aspects. The ANOVA test was conducted to determine whether there were differences in the subcategories of seismic literacy depending on age, risk awareness, and final educational background. The results show that participants 20 years of age and below, and 60 years of age and above, declared a higher level of seismic literacy than participants between 30 and 50 years of age. Participants in their 30s declared the lowest level of seismic awareness and participants in their 40s declared the lowest level Remote Sens. 2021, 13, 1365 20 of 25 of seismic knowledge and management among the different age groups. However, there was no significant difference in earthquake literacy in all aspects among all ranges of ages (under 20s to over 60s). Only participants who were aware that where they live is absolutely safe from earthquakes had a statistically significantly higher level of seismic management behavior. They also declared a higher level of seismic knowledge. Participants who were aware that where they live is not safe at all from earthquakes declared a higher level of seismic awareness but declared the lowest level of seismic knowledge and management. This result is an important warning sign for local and regulatory authorities to raise the citizen seismic knowledge and management. In the latter case, participants in graduate school and graduate school graduates declared a higher level of seismic literacy in all three aspects. Participants with an education level of less than high school declared the lowest level of seismic knowledge and high school graduates declared the lowest level of seismic awareness and management. Further analysis was conducted to find significant differences through Tukey's HSD post hoc test and showed that citizens who were aware that where they live is absolutely safe from earthquakes had a statistically significantly higher level of seismic management behavior. Participants with a highly educated background had the highest level of seismic knowledge.
The participants of 60 years old and above declared the highest level of seismic literacy. Some authors posit that this could be explained because adults in this stage of life have acquired greater experience and care responsibilities (either for others or their assets), which may give rise to increased interest in involving themselves in preparedness measures [65]. Participants' educational background also influenced their seismic literacy, especially in seismic knowledge, in line with previous studies [66,67]. Groups with less seismic literacy should be a target of intervention in order to raise earthquake risk awareness and motivate them to adopt preparedness actions.
Nevertheless, in general, the participants declared relatively low earthquake literacy, including seismic knowledge, awareness, and preparedness. This means that they are not properly prepared for earthquake hazards. For the areas that were classified as highly vulnerable to earthquakes, the focus should be on building reconstruction and retrofitting [16], whereas for those in moderate-risk areas, the focus should be on building renovation to reduce their seismic vulnerability. Areas with high building and population densities, such as Daegu, Busan, and Ulsan, should develop education programs to improve earthquake risk awareness and seismic literacy. The government could increase earthquake literacy, starting with a conscious approach to Korean citizens, especially in high vulnerability areas.
Finally, the institution responsible for developing local disaster risk reduction plans and programs should appropriately characterize their target audiences and areas if they expect to obtain more effective and efficient results. We expect that the results reported in this study will be useful input to achieve this. However, this study has certain limitations. First, the questionnaire could provide more information of the participants profile such as marital status, children in the house, annual household income, and household type and could be analyzed to provide specific information about participants' preparedness. Future studies could investigate the interaction between these variables to find specific patterns of seismic literacy.

Conclusions
We conducted a survey of 1256 participants to investigate the seismic literacy of Korean citizens, including seismic knowledge, awareness, and management, using a questionnaire of citizen earthquake literacy (CEL) following the 2017 Pohang earthquake. The results declared that the citizens had low literacy, which means that they are not properly prepared for earthquake hazards. To develop an earthquake risk reduction plan and program efficiently and effectively, not only must one appropriately characterize the target audience, but also indicate high potential earthquake zones and potential earthquake damage. Therefore, this study mapped and analyzed the seismic vulnerability in southeast Korea using LogitBoost, logistic model tree (LMT), and logistic regression (LR) machine learning algorithms based on a building damage inventory map. The building damage locations were generated after the 2017 Pohang earthquake using the damage proxy map (DPM) method from the Sentinel-1 synthetic aperture radar (SAR) data. The DPMs manage to detect coherence loss, which indicates damaged buildings in residential and commercial areas due to the Pohang earthquake and show a good correlation with the Korea Meteorological Administration (KMA) report with modified Mercalli intensity (MMI) scale values of more than VII (seven). The damage locations were randomly divided into two datasets: 50% for training the vulnerability models and 50% for validating the models in terms of accuracy and reliability. Fifteen seismic-related factors were used to construct a model for each algorithm. Model validation based on the area under the receiver operating curve (AUC) was used to determine model accuracy. The AUC values of seismic vulnerability maps using the LogitBoost, LMT, and LR algorithms were 0.769, 0.851, and 0.749, respectively. We suggest that earthquake preparedness efforts should focus on reconstruction, retrofitting, renovation, and seismic education in areas with high seismic vulnerability in South Korea. The results of this study are expected to be beneficial for engineers and policymakers aiming at developing disaster risk reduction plans, policies, and programs due to future seismic activity in South Korea.

Conflicts of Interest:
The authors declare no conflict of interest.