After the 2016 Gyeongju earthquake and the 2017 Pohang earthquake struck the Korean peninsula, securing financial stability regarding earthquake risks has become an important issue in South Korea. Many domestic researchers are currently studying potential earthquake risk. However, empirical analyses and statistical approaches are ambiguous in the case of South Korea because no major earthquake has ever occurred on the Korean peninsula since the Korean Meteorological Agency started monitoring earthquakes in 1978. This study focuses on evaluating possible losses due to earthquake risk in Seoul, the capital of South Korea, by using a catastrophe model methodology integrated with GIS (Geographic Information Systems). Building information, such as structure and location, is taken from the building registration database and the replacement cost for buildings is obtained from insurance information. As the seismic design code in the KBC (Korea Building Code) is similar to the seismic design code of the UBC (Uniform Building Code), the damage functions provided by HAZUS-Multi-hazard (HAZUS-MH) are used to assess the damage state of each building in event of an earthquake. A total of 12 earthquake scenarios are evaluated by considering the distribution and characteristics of active fault zones on the Korean peninsula and damages, with total loss amounts are calculated for each of the scenarios. The results of this study show that loss amounts due to potential earthquakes are significantly lower than those of previous studies. The challenge of this study is to implement an earthquake response spectrum and to reflect the actual asset value of buildings in Seoul.
On 15 November 2017, an earthquake of
Modified Mercalli Intensity (MMI) map of the Pohang earthquake
The earthquakes that occurred in Gyeongju and Pohang were caused
by the Yangsan fault zone, classified as an active fault, on the Korean
Peninsula, which has the ability to generate earthquakes with a maximum intensity of
This study uses catastrophe model methodology to predict losses and damage to buildings and their contents from a potential earthquake that could occur in Seoul. This study differs from the previous studies in that it implements the actual building and insurance data and the observed seismic characteristics of South Korea. The detailed information of approximately 630 000 buildings across Seoul is acquired through the building registration database. The replacement cost of each building and its contents is statistically estimated by using an insurance database that is classified by occupancy to simulate the real situation in South Korea.
Predicting loss amount from a potential disaster using a catastrophe model differs from the actuarial approach model. While the actuarial technique estimates the loss based on empirical data, the catastrophe model generates disaster scenarios based on a scientific understanding of disasters and assesses the loss amount from an event scenario. For possible earthquakes in South Korea, it is appropriate to use the catastrophe model for predicting losses because empirical data from earthquakes on the Korean peninsula are too scarce to enable actuarial processing. The definition and procedure of the catastrophe model can marginally differ between researchers or suppliers, but the conventional procedure can be summarized as in Fig. 2.
Procedure of a catastrophe model. Source: modified from Parodi (2014).
As shown in Fig. 2, the catastrophe model has a four step process. (i) An information database of that property that may be exposed to disaster should be constructed. However, the exposure data sets in previous studies are typically available at relatively coarse resolutions because they are accompanied by difficulties related to limited resources or privacy issues, among others (Dell'Acqua et al., 2012; Figueiredo and Martina, 2016). In order to overcome these limitations, this study used the building registration database of South Korea to build an exposure data set. Detailed information of the building must be recorded, which is registered in the building registration database whenever the building is constructed or reconstructed according to the Building Act of South Korea. In this study, the detailed information needed to evaluate the vulnerability of all buildings in Seoul was extracted from the building registration database. The extracted data are classified into 36 structure types and 33 occupancies (the same as the building types of HAZUS-MH) and divided into three seismic codes estimated based on a comprehensive consideration of the construction year, total building area, and occupancy. The details of the classification of the 36 structure types and 33 occupancies are shown in Tables 1 and 2, respectively.
Model building types.
Source: Federal Emergency Management Agency (2013)
Building occupancy classes.
Source: Federal Emergency Management Agency (2013)
(ii) A hazard module for the generation of a physical hazard map from a simulated earthquake event should be developed. For example, the peak ground acceleration can be represented as hazard intensity in an earthquake hazard map. The seismic events are usually generated by stochastic methodologies, such as Monte Carlo simulation. However, it has been less than 40 years since earthquakes began to be monitored on the Korean Peninsula, and there was no large-scale earthquake in Seoul during the monitoring period. In this study, synthetic earthquakes were generated considering the activity of the active faults passing through Seoul, and the seismic-hazard map was prepared by selecting the attenuation relation that most closely resembles the Gyeongju earthquake and the Pohang earthquake from many attenuation relations that have been proposed by many Korean and international researchers.
(iii) A vulnerability module for assessing the damage state of individual properties should be prepared by combining the information of exposed property and hazard intensity. The probabilities of each damage state should be estimated from the spectral displacement of each building due to the seismic impact on the vulnerability module. The spectral displacement is determined by performance point, which is the intersection of the demand curve and the capacity spectrum.
(iv) The financial module is implemented to quantify the damage of individual buildings into a monetary loss to predict a total loss amount. In order to estimate the repair cost of a building due to a seismic impact, it is necessary to ascertain the replacement value or the current assets of the building calculated in cost mode. In this study, the values of building, contents, and inventories of the representative building in each category were estimated by statistical processing of the appraisal data for insuring property.
The detailed information of each building, such as location, structure, size, floor area, construction year, occupancy, and other minor considerations influencing the seismic response, is obtained from a computerized database of building registration records. There are presently about 630 000 buildings within Seoul city as of 2016, according to building registration records. These buildings are classified as residential (76 %), commercial (20.3 %), industrial (0.5 %), and other (3.2 %), which includes government and educational institution buildings. The residential category is dominated by masonry structures of fewer than five stories and concrete structures are dominant in the commercial category. A total of 82 % of the buildings in Seoul were built before 1988 when the seismic building code began being implemented. Table 3 summarizes the statistical characteristics of buildings in Seoul.
A summary of the statistical characteristics of buildings in Seoul.
The replacement costs for each building damaged by the earthquake and their contents are estimated based on statistical processing of 1500 records of asset evaluation data for property insurance and a construction cost table (Korea Appraisal Board, 2016). After processing, the total replacement cost of buildings and contents is estimated to be about USD 900 billion in Seoul, and approximately 72 % of buildings in Seoul were estimated to have a replacement cost of USD 0.1–1 million. The indirect costs and losses attributed to land, intangible assets, and business interruption are not considered in this study.
In the circum-Pacific seismic zone, South Korea is considered safer and
less prone to earthquakes compared to neighboring areas, such as
Japan, China, and Taiwan. However, many domestic researchers insist that
there are two representative active faults in South Korea. One of these, the
Yangsan Fault, caused the Gyeongju earthquake. The second fault, Chugaryeong
Fault, is centrally located on the Korean peninsula (Choi et al., 2012; South Korean Ministry of Public Safety and Security,
2012; Chung et al., 2014). The Chugaryeong Fault crosses the eastern side of
Seoul and is believed to have caused the 2010 earthquake, with an intensity of
Chugaryeong fault zone in the middle of the Korean peninsula. Source: Modified from Chung et al. (2014).
The Chugaryeong Fault has similar activity to Yangsan Fault, which has the
capacity to cause an earthquake of
To construct each hazard map from each earthquake event scenario, it is important to understand the attenuation relationship of ground motions from epicentral distance. The ground motion can be characterized by peak ground acceleration (PGA) and spectral responses based on a response spectrum shape.
A lot of experimental attenuation formulas for estimating PGA have been
developed by means of regression analysis (Atkinson and Boore, 1997; Toro et
al., 1997; Atkinson and Silva, 2000; Lee and Kim, 2002; South Korean Ministry of Public Safety and Security, 2012).
However, in choosing the attenuation formula, a careful approach is needed
since the impact of the formula on the estimated amount of
earthquake loss is very large. The South Korean Ministry of Public Safety and Security (2012) proposed three attenuation formulas for the Korean
Peninsula, which are expressed following Eq. (1). The
attenuation equation (Eq. 1) proposed by the South Korean Ministry of Public Safety and Security (2012) is used
in this study. The attenuation formula of the MPSS requires four coefficients
(or fitting parameters). In this study, the four coefficients in Eq. (1) of
The influence of the seismic attenuation equation on the seismic-hazard map
is great, but the reliability of the attenuation equations presented so
far remains controversial. In this study, we tried to utilize the results of
domestic studies reflecting the seismic characteristics in South Korea, and the
attenuation equation of Eq. (1), which is considered to be the most
conservative formula because its attenuation of the seismic wave is the lowest of
the formulas proposed by the South Korean Ministry of Public Safety and Security (2012), is chosen for building the earthquake
hazard map from the event scenario. The hazard maps, according to the each
scenario, are shown in Fig. 4. PGA in Seoul ranges from 0.06 to 0.7
PGA hazard map of Seoul according to each scenario event.
The severity of the vibratory response of buildings to an earthquake impact depends
on the relationship between the characteristics of ground motion, described
as the response spectrum (which has a different shape according to ground
conditions) and structural characteristics of the building. But since the design
response spectrum currently used in South Korea is based on a high-seismicity
region like California, a lot of domestic researchers insist that the
spectrum is different due to the characteristics of earthquakes on the Korean
peninsula (Kim et al., 1998; Han, 2003; Hwang et al., 2015; Lee and Ju,
2017). In general, the ground conditions of South Korea, including Seoul, are
characterized by shallow bedrock, and the earthquakes that have occurred in South Korea have
a duration of strong motion that is shorter than those in
high-seismicity regions. In the case of the Gyeongju earthquake, strong
motion with a short duration of 0.1–0.2 s was also
observed. As shown in Fig. 5, the shape of the standard response spectrum is
described by four transition points.
Shape of the standard response spectrum.
Factors of the standard response spectrum in this study.
In 1988, when an earthquake occurred in Mexico, the seismic design code in South Korea
was first mandated for buildings with six and more stories or floor areas of
100 000 m
The seismic design codes in the KBC were established based on the Uniform Building Code (UBC), the Applied Technology Council (ATC), and the International Building Code (IBC) (Seoul Metropolitan Government, 2012; Lee, 2015). While the seismic zone in the UBC is divided into six zones that each have a regional factor, which is defined as the peak ground acceleration (PGA). The zone in KBC is divided into two zones: Zone I, which includes the Seoul area, and Zone II. The regional factor of Zone I was 0.11 before 2009 but strengthened to 0.22 after that, and the seismic design code of buildings built in Seoul before 2009 is similar to Zone 2A of the UBC and the code of buildings built after 2009 is similar to Zone 2B of the UBC.
Alteration of seismic design code in South Korea.
As mentioned above, since the seismic design code of South Korea is similar to the UBC
and the ATC code, the damage functions proposed by HAZUS-MH can be applied
to estimate building damage due to seismic impact. The damage function for
each building type in HAZUS-MH includes two types of damage curves: capacity
curve and fragility curve. The capacity curve is used to determine the peak
building response from the capacity spectrum method. This method is a
schematic procedure for comparing the capacity curve obtained by pushover
analysis with the demand spectrum of ground motion on the Acceleration
Displacement Response Spectrum (ADRS). Thus, the response spectrum has to be
converted to the demand spectrum for representing the relationship between
spectral displacement and spectral acceleration. Equation (2), proposed by
HAZUS-MH, can relate spectral acceleration with spectral displacement for a
given period value (Federal Emergency Management Agency, 2013).
The intersection of the capacity curve and the demand spectrum is a performance point, which can evaluate the associated damage state for the structure and compare that damage state for different earthquakes (Fig. 6).
Performance points according to the intersection of the capacity curve and demand spectrum.
The fragility curves estimate the probability of exceeding different damage
states given peak building response, which is represented as spectral displacement or
spectral acceleration at the performance point. The damage state is divided into
four states: slight, moderate, extensive, and complete. Each fragility curve
is expressed as a lognormal function defined by a median value of peak
building response, corresponding to the mean threshold of associated damage
state by a logarithmic standard deviation (
Example of the fragility curve for structural component.
The nonstructural components of the building are divided into drift-sensitive components and acceleration-sensitive components. In general, while architectural components such as interior or exterior walls are more drift-sensitive, the mechanical and electrical components of the building are acceleration-sensitive. Therefore, the function of inter-story drift is used to estimate the damage state of drift-sensitive components and the function of floor acceleration is used to estimate the damage state of acceleration sensitive components or contents in the building.
The capacity curve and the fragility curve in the HAZUS-MH are classified into high-code, moderate-code, low-code and precode buildings as per seismic design codes (Federal Emergency Management Agency, 2013). When all buildings in Seoul are classified by comparing the seismic design code of the HAZUS-MH, it is estimated to be approximately 91.7 % precode, 5.4 % low-code and 2.9 % moderate-code buildings.
Estimated mean loss ratio of building based on occupancy type.
Estimated mean loss ratio of building based on structure type.
Using estimates that include the structural and nonstructural repair costs caused by building damage and the associated loss of building contents and business inventory, provided by HAZUS-MH, the probability of exceeding different damage states for the each component can be converted to the loss ratio of replacement cost for the purpose of evaluation of direct economic loss. Tables 6 and 7 summarize the estimated mean loss ratio of buildings, which includes structural and nonstructural components, and contents, depending on occupancy and structure type. It is a common pattern that the majority of damage is inflicted on low-rise residential buildings made of masonry in South Korea, where earthquakes are characterized by strong short-period components are dominant. However, this common pattern is not clearly shown in the result of this study and there are two main reasons for that result.
The first reason is that although the number of low-rise residential buildings made of masonry is much larger, the total asset value is much lower than high-rise residential buildings made of reinforced concrete, such as apartments, which are generally classified as a luxury residence, while low-rise masonry houses are classified as a low-priced residence in Seoul. The second reason is that the nonstructural elements, such as mechanical and electrical components that are more vulnerable to ground shaking than structural components, are of a higher proportion in the buildings made of concrete and steel than in masonry.
Loss ratio map for each building by the scenario.
Figure 8 is a map that shows the loss ratio of each building in the Gangnam
district, located 3 km from the epicenter, with an intensity of
Aggregated loss amount due to each scenario.
The total loss amount for each scenario can be simulated by replacement
cost, seismic intensity, damage function, and other factors mentioned above.
But since systematic data issues or biases across a portfolio can result in
losses being consistently under- or over-simulated (Lloyd's Market Association, 2017), the results
need to be corrected by comparing empirical data. The Linear Scaling Method
(LSM), which is one of the common methods of correcting systematic errors, can
be used to calibrate pre-simulated loss amount. The LSM reflects the difference
between pre-simulated results and observed results in the simulated results
as shown by Eq. (4).
The existence of active fault zones on the Korean peninsula and the recent earthquakes that affected Gyeongju and Pohang have made experts question whether current overall practices would still be adequate if a similar earthquake were to occur in Seoul. The concentration of major industrial and commercial facilities carries a significant inherent risk of causing catastrophic loss of life and economic damage and poses a significant administrative challenge for disaster management in South Korea. The disaster management is divided into four phases: (1) mitigation, (2) preparedness, (3) response, and (4) recovery. At each phase, which each have particular needs and problems, different strategies and support are required to force social resilience against each natural disaster. It is also important that the activities at each phase generate a virtuous cycle of practices and assist in making each phase stronger.
The development of the insurance industry could be used as a good example for
explaining virtuous cycles in disaster management. The insurance industry as disaster
risk financing commonly plays a major role in securing financial stability for
a smooth recovery from natural disasters. However, it also helps these
activities to perform more effectively during the other phases such as
mitigation, preparedness, and response. The Sichuan earthquake of 2008 is in
stark contrast with the New Zealand earthquake of 2010 in terms of disaster
management efficiency, due to limited insurance coverage. The Sichuan, China,
earthquake of
Most domestic insurers believe that it is impossible to predict loss amount from potential earthquakes and that it is difficult to quantify the earthquake risk in South Korea. This belief of insurers is a major obstacle to the development of earthquake insurance programs. However, as mentioned above, various studies required for the catastrophe model methodology have either been completed or are in progress by various domestic researchers and a lot of databases related to potential earthquake risk in South Korea is being accumulated. Compared to other studies, this study is different in that real insurance information and building registration databases are used to predict the loss amount from a potential earthquake. It not only helps advance the prediction process but also serves the insurer to better understand and estimate the earthquake risk. This study shows that risks due to potential earthquakes in South Korea are significant and that the insurance industry can support more detailed studies for a better understanding of insurance risk and the expanding scope of current insurance practices for earthquake risks. Because of this, insurance companies have an opportunity to further explore currently underused areas of business in property insurance.
Access to the Seoul building and asset data used for the analysis are restricted as they are the property of a private company. Part of the calculation coding used for data analysis can be found in the Supplement.
The supplement related to this article is available online at:
WC conceived and wrote the entire paper. JWP and JK analyzed the results and revised the paper.
The authors declare that they have no conflict of interest.
This work was supported by the National Research Council of Science and Technology (NST) grant by the South Korean government (MSIP) (no. CRC-16-02-KICT).
This paper was edited by Maria Ana Baptista and reviewed by two anonymous referees.