Next Article in Journal
Autonomous Underwater Vehicles: Localization, Navigation, and Communication for Collaborative Missions
Next Article in Special Issue
Prediction of Structural Type for City-Scale Seismic Damage Simulation Based on Machine Learning
Previous Article in Journal
Experimental and Simulation Study of Micro-Channel Backplane Heat Pipe Air Conditioning System in Data Center
Previous Article in Special Issue
Study on Seismic Performance of a Mold Transformer through Shaking Table Tests
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Numerical Model for Simulating Ground Motions for the Korean Peninsula

Department of Architectural Engineering, Hanyang University, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(4), 1254; https://doi.org/10.3390/app10041254
Submission received: 7 January 2020 / Revised: 29 January 2020 / Accepted: 7 February 2020 / Published: 13 February 2020

Abstract

:

Featured Application

Developing a numerical model for simulating ground motions for low-to-moderate seismic regions such as the Korean Peninsula.

Abstract

Ground motions are used as input for the response history analyses of a structure. However, the number of ground motions recorded at a site located in low-to-moderate seismic regions such as the Korean Peninsula is few. In this case, artificial ground motions need to be used, which can reflect the characteristics of source mechanism, travel path, and site geology. On 15 November, 2017, the Pohang earthquake, with a magnitude of 5.4 and a focal depth of 9 km, occurred near the city of Pohang. This earthquake caused the most significant economic loss among the earthquakes that occurred in the Korean Peninsula. During the Pohang earthquake, valuable ground motions were recorded at stations distributed in the Korean Peninsula. In this study, a ground motion model is proposed based on ground motions recorded during the 2017 Pohang earthquake. The accuracy of the proposed model is verified by comparing measured and simulated ground motions at 111 recording stations.

1. Introduction

The Korean Peninsula is located in stable continental regions about 400 km from the boundaries of four plates: the Philippine Sea, Pacific, North American, and Eurasian Plates [1,2]. The seismic activity in this region is lower than that in regions located near plate boundaries such as Japan, Indonesia, or the Western United States. However, Korean historical documents reported that big earthquakes with a magnitude of 6 or larger occurred in the peninsula from 2 A.D. to 1904 A.D. [3]. The Pohang earthquake, with a magnitude of 5.4 and focal depth of 9 km, occurred near the city of Pohang, located in the Korean South-East province. This earthquake caused the largest casualties and economic loss among the earthquakes that have occurred in Korea.
To protect structures from earthquakes, it is necessary to conduct the seismic performance evaluation of structures with reliable input ground motions, and to retrofit the structures based on the results of the seismic performance evaluation [4,5,6]. However, it is difficult to collect available ground motions recorded from mid- to large-size earthquakes that have occurred in sites located in low-to-moderate seismic regions. In this case, artificial ground motions can be simulated and used in seismic performance evaluation. Simulated ground motions should retain the characteristics of the local source mechanism, travel path, and geology of a site.
Previous studies [7,8,9,10,11,12] have developed and improved the ground motion simulation model to generate artificial ground motions by considering key components of seismological characteristics. Ground motions can be simulated using the stochastic point-source model in the frequency domain and the shaping window model in the time domain.
In this study, a ground motion simulation model is proposed. The ground motion parameters are determined based on the ground motions recorded from the mainshock of the 2017 Pohang earthquake.
To verify the accuracy of the proposed model, the key components of measured and simulated ground motions are compared, such as peak ground acceleration (PGA) and 5% damped pseudo spectral acceleration (PSA).
In order to develop the proper ground motion simulation model for the Korean Peninsula, this article is organized into three sections. In Section 2, the collected ground motions in three orthogonal directions (East–West, North–South, and Vertical) are presented, which were recorded from 111 stations in the inland Korean Peninsula during past earthquakes. In Section 3, a stochastic point-source model and shape window model are proposed to simulate the ground motions recorded during past earthquakes that occurred in the Korean Peninsula. In Section 4, the proposed model is verified by comparing the generated artificial ground motions and ground motions recorded during the 2017 Pohang earthquake event.

2. Ground Motions Collected for Developing the Ground Motion Simulation Model

To develop a numerical model for simulating ground motions in the Korean Peninsula, the ground motions recorded at 111 seismic stations during the 2017 Pohang earthquake were collected from the National Earthquake Comprehensive Information System (NECIS) of Korea Meteorological Administration (KMA). Figure 1 shows the distribution of 111 stations out of 157 stations in Korea, which provided ground motion records with a high signal-to-noise ratio (SNR). To exclude ground motions with signal distortion due to noise, we only collected recordings with a signal-to-noise ratio (SNR) greater than 2.0. Each station provided three components (North–South, East–West, and Vertical) of ground motions. Thus, the total number of recorded ground motions is 333 (=3 × 111).
All collected ground motions were recorded with a sampling rate of 0.01 s (100 samples per second), and the NECIS provided these records after applying a high cut anti-aliasing filter to remove the noise with a frequency greater than 40 Hz. In this study, low-frequency noise was also removed by using the 0.1 Hz (corresponding period 10 s) low-cut filter, and baseline correction was done using a technique proposed by Papazafeiropoulos and Plevris [13]. Table 1 summarizes the information on the 2017 Pohang earthquake.
Figure 2 shows the East–West (E–W) direction component of ground acceleration ( u ¨ g ), velocity ( u ˙ g ), and displacement ( u g ) at station BUS2. Figure 2a–c show raw record data for ground acceleration, velocity, and displacement. Figure 2d–f show the records after applying a low-pass filter, whereas Figure 2g–i show ground motions after applying base line correction. Figure 2j–l show ground motions after applying a low-cut filter, as well as base correction. It can be observed that ground motions are properly adjusted by using both a low-cut filter and base correction.
The site amplification effect was also removed from the ground motion records because a ground motion model has been developed for hard rock conditions without considering site effects. In general, site effects are considered in numerical models by applying a site amplification factor to the ground motions generated for hard rock site conditions.
In this study, the site amplification effect was removed from the recorded ground motions by using a horizontal-to-vertical spectral ratio (HVSR) technique, which has been widely used for site amplification factor calculation [14,15,16,17].
The site amplification factor, A M P ( f ) , can be calculated from Equation (1).
A M P ( f ) = P S A E W ( f ) × P S A N S ( f ) P S A V ( f )
where P S A E W ( f ) , P S A N S ( f ) , and P S A V ( f ) are the 5% damped pseudo spectral acceleration (PSA) from the East–West, North–South, and Vertical components of ground motions, respectively. Figure 3 shows the calculated site amplification factor function estimated from the ground motions recorded at stations BUS2, CEA, and CEJA during the 2017 Pohang earthquake. This function was applied to the Fourier spectrum of a ground motion at each frequency. Figure 4 shows the ground acceleration time histories and Fourier spectra before and after applying Equation (1).

3. Ground Motion Simulation Model

The numerical model to simulate ground motions in the Korean Peninsula is proposed based on the stochastic point-source model in the frequency domain and the shaping window model in the time domain. The point source model is an effective tool for generating ground motions for a site with limited seismological information.
The Fourier amplitude spectrum (FAS) of a ground motion can be determined using the stochastic point-source model [8,9,10,11,12], which was developed based on the omega-square model [18,19]. This model consists of three main seismological characteristics, namely source, path, and site effects [12].
The shaping window model consists of a ground motion envelope and duration time. Saragoni and Hart [7] proposed the envelope of ground motion records, which was adopted in this study. The duration of a ground motion affects the incidence of the collapse of a structure [20].

3.1. Stochastic Point-Source Model Estimation

In the stochastic point-source model, the Fourier amplitude spectrum [ A ( f ) ] of a ground motion is calculated using Equation (2).
A ( f ) = S o u r c e ( M 0 , f , f c ) × P a t h ( R H , f ) × S i t e ( κ 0 , f )
where S o u r c e ( M 0 , f , f c ) is the earthquake source effect function, P a t h ( R H , f ) is the path effect function, S i t e ( κ 0 , f ) is the site effect function, M 0 is the seismic moment, R H is the hypocentral distance, f is the frequency, and f c is the corner frequency. A detailed explanation of individual functions is summarized in Table 2.
The component parameters κ o , M o , and f c to calculate A ( f ) were estimated from the ground motions recorded at individual stations during the 2017 Pohang earthquake, whereas the values of other component parameters were adopted from previous studies, as listed in Table 2.
To effectively determine the values for κ o , M o , and f c , the A ( f ) of the recorded ground motions were smoothed using a technique proposed by Konno and Ohmachi [21]. Figure 5 shows the original and smoothed A ( f ) of ground motions recorded at three sample stations (BUS2, CEA, and CEJA) during the 2017 Pohang earthquake.
In this study, the site attenuation coefficient, κ 0 , was determined first. Since the contribution of κ 0 in Equation (2) is determined without interaction with path effect and site amplification, the contribution of path effect and site amplification were removed from Equation (2). The Fourier amplitude spectrum without path effect and site amplification [ A ( f ) ] was calculated using Equation (3).
A ( f ) = A ( f ) P a t h ( R H , f ) × A M P ( f ) .
Anderson and Hough [25] estimated the site attenuation coefficient ( κ 0 ) using a residual slope of A ( f ) between 10 Hz and 40 Hz, which was also used in this study. Figure 6 shows A ( f ) and A ( f ) . The slope ( κ 0 ) of A ( f ) is also plotted in Figure 6. The values of κ 0 estimated for individual stations are summarized in Table A1 of Appendix A.
Previous studies [9,26,27,28] estimated M o and f c by using the omega-square model proposed by Brune [18,19] based on the FAS of the recorded ground motions. To determine M o and f c without interaction with site attenuation ( κ o ), the Fourier amplitude spectrum without the effect of κ o [ A ( f ) ] was calculated using Equation (4). Figure 7 shows A ( f ) and A ( f ) .
A ( f ) = A ( f ) exp ( π κ 0 f )
Corner frequency, f c , can be estimated using Equations (5)–(8), which were proposed by Andrews [29] and Jo and Baag [27].
f c = ( J 2 π 3 Ω 0 2 )
Ω 0 = 2 ( K 3 J ) 1 4
J = 2 3 ( Ω 0 ω 1 ) 2 f 1 + 2 f 1 f 2 | ω A ( ω ) | 2 d f + 2 | ω 2 A ( ω 2 ) | 2 f 2
K = 2 | ω 1 A ( ω 1 ) | 2 f 1 + 2 f 1 f 2 | A ( ω ) | 2 d f + 2 3 | A ( ω 2 ) | 2 f 2
Seismic moment, M 0 , can be also calculated from Equation (9), which was proposed by Joshi et al. [28].
M 0 = 4 π ρ β S 3 × Ω 0 × R r e f R θ ϕ × F × V
Figure 8 shows the values of f c and M 0 for 111 stations distributed within the Korean Peninsula. Table A2 in Appendix A lists the values of f c and M 0 estimated for 111 stations. In this study, the median values ( f ^ c , M ^ o ) of f c and M 0 were calculated as 0.58 Hz and 8.39 × 10 24 dyne-cm, respectively, which were used to calculate A ( f ) (Equations (2)–(4)). Figure 9 shows the measured and calculated A ( f ) for three sample stations. It can be observed that the calculated A ( f ) matches the measured A ( f ) , by which the accuracy of the proposed procedure is verified.

3.2. Shaping Window Model Estimation

In order to generate a ground acceleration in the time domain ( u ¨ g ( t ) ) from A ( f ) , a proper shaping window function should be applied [12]. A shaping window model consists of a ground motion envelope shape and duration time. Boore [12] simulated ground motion recordings using an envelope shape proposed by Saragoni and Hart [7] and a duration time proposed by Atkinson and Boore [11]. However, the shaping window model used by Boore [12] may not properly reflect the envelope shape of the ground motions recorded in the Korean Peninsula. In this study, a shaping window model ( W ( t ) ) is proposed based on the ground motions recorded during the 2017 Pohang earthquake (Equation (10)).
ln ( W ( t , T D ) ) = c 0 + c 1 ln ( t T D ) + c 2 t T D
where t is the time, and T D is the duration time. The duration time is calculated from the S-wave arrival time to the time corresponding to the 95% normalized arias intensity energy ( E ( t ) ). E ( t ) is estimated using Equation (11).
E ( t ) = 0 t u ¨ g 2 ( t ) d t E T
where E T is the total energy of u ¨ g ( t ) . Figure 10a–c show the estimated duration time of a ground motion recorded at the BUS2 station, E ( t ) , and estimated shaping window, respectively. Figure 11 shows the estimated shaping window and duration of ground motions recorded at 111 stations during the 2017 Pohang earthquake. The coefficients c 0 , c 1 , and c 2 in Equation (10) are estimated to match the median envelope shape denoted with a thick solid line in Figure 7a.
Regression analyses were conducted, by which the values of c 0 , c 1 , and c 2 were determined to be 1.6546, 0.6227, and −3.2663, respectively. The envelope shape obtained from the 2017 Pohang earthquake is significantly different from that calculated using the equations proposed by Boor (2003). The envelope shape constructed from the equation proposed by Boor (2003) significantly underestimates the measured envelope shape of ground motions recorded during the 2017 Pohang earthquake.
Figure 11b shows the estimated duration time ( T D ) at each station. Previous studies [10,11,30] reported that duration time is mainly affected by corner frequency ( f c ) and source-to-site distance ( R H ). In this study, the following equation is proposed for T D , according to R H and f c . Figure 11b shows the duration time calculated using Equation (12). This figure also shows that T D calculated from Equation (12) distinctively differs from that calculated using the equations proposed by Atkinson and Boore [11]. The equation of T D developed for the northeastern United States could underestimate T D obtained from the 2017 Pohang earthquake. This observation reveals the importance of why numerical models should be developed for low-to-moderate seismic regions, such as the Korean Peninsula.
T D = 1 f c + { 3.256 ( R H 10   km ) 0.247 + 0.350 R h y p o ( 10   km < R H 50   km ) 19.522 0.045 R h y p o ( 50   km < R H 100   km ) 9.005 + 0.060 R h y p o ( R H > 100   km )

4. Ground Motion Simulation for the 2017 Pohang Earthquake

Ground motions were generated using the proposed numerical model, consisting of the stochastic point-source model and shaping window model. The procedure is briefly summarized as follows:
(1)
White noise is first generated in time domain for a duration time of the ground motion (Equation (12)).
(2)
The noise is then windowed using the shaping window model (Equation (10)).
(3)
The windowed noise is transformed into the frequency domain.
(4)
The A ( f ) of the noise is normalized by the square-root of the mean square of A ( f ) in all frequencies.
(5)
The normalized A ( f ) is then multiplied by the stochastic simulation model (Equation (2)).
(6)
The resulting A ( f ) is transformed back to the time domain, which is a simulated ground motion.
Figure 12a,b show two horizontal components’ ground motions at station BUS2 during the 2017 Pohang earthquake, whereas Figure 12c presents one ground motion simulated using the proposed procedure. Figure 12d presents the geomean of A ( f ) for recorded ground motions. In this figure, A ( f ) calculated using Equation (2) is also included, which shows that the calculated A ( f ) matches that obtained from recorded ground motions. Figure 12e shows the geomean of 5% damped-pseudo spectral acceleration P S A ( T ) and the median P S A ( T ) values of 1000 simulated ground motions. The difference in P S A ( T ) between recorded and simulated ground motions is generally small. Similar observations were made for the ground motions for stations CEA and CEJA.
To verify the accuracy of the proposed numerical model, ground motions were simulated at 111 stations for the 2017 Pohang earthquake and their P S A ( T ) values were calculated. Residuals induced by the difference in simulated and recorded P S A ( T ) ( P S A ( T ) s i m P S A ( T ) m e a s u r e d ) were calculated using Equation (13), with a period ( T ) range between 0.01 s and 10 s.
Residual ( T ) = log [ P S A ( T ) s i m ] log [ P S A ( T ) m e a s u r e d ]
Residuals were calculated for all ground motions of the 111 stations. Figure 13 shows the residuals according to hypocentral distance. As shown in Figure 13a, the mean value of residuals is near zero, which indicates that the ground motions in the Korean Peninsula can be adequately simulated using the proposed numerical model with the source, path effect, and site amplification functions.
If the path effect function proposed by Boor (2003) is used in the numerical model, the mean residual deviates from zero, as shown in Figure 13b. The degree of deviation varies according to the period. In the case of using the source function and site amplification function proposed by others when simulating ground motions with the proposed numerical model, similar observations were made, which are not included in the paper.

5. Conclusions

In this study, a numerical model was proposed to simulate ground motions in the Korean Peninsula. The model consisted of a stochastic point-source model and a shaping window model developed based on ground motions recorded at 111 stations during the 2017 Pohang earthquake. Conclusions obtained from this study are as follows:
(1)
Source, path, and site effect functions were developed for the stochastic point-source model, which reflected the seismological characteristics of the Korean Peninsula.
(2)
To generate ground motions in the time domain that represents ground motions recorded in the Korean Peninsula, an envelope shape and a duration time function were proposed based on the ground motions recorded at 111 stations.
(3)
In order to verify the accuracy of the proposed numerical model, residuals measuring the difference in P S A ( T ) between recorded and simulated ground motions were calculated for 111 stations. It was observed that ground motions in the Korean Peninsula were simulated accurately using the proposed numerical model, which included proper source, path, and site effect functions.
(4)
The results of this study reveal the potential of the proposed numerical model to simulate input ground motions for low-to-moderate seismic regions, such as the Korean Peninsula.

Author Contributions

Methodology, S.W.H.; investigation, H.W.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Korea Agency for Infrastructure Technology Advancement, grant number 20CTAP-C152179-02.

Acknowledgments

This research was supported by grants from the Korea Agency for Infrastructure Technology (20CTAP-C152179-02).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. κ 0 values for 111 stations.
Table A1. κ 0 values for 111 stations.
No.Station κ 0 No.Station κ 0 No.Station κ 0 No.Station κ 0
1ADO20.009229GOCB0.020357JEU2−0.003585SCHA0.0161
2ADOA0.028430GSG−0.011258JINA0.019786SEHB0.0259
3BON0.018231GUM0.044759JMJ0.014987SEO20.0059
4BOSB0.015032GUS0.001760JNPA0.022188SES20.0132
5BSA0.020833GUWB0.023561JUR0.012089SHHB0.0149
6BURB0.025434GWJ0.028962KAWA0.022790SKC20.0102
7BUS0.010735GWYB0.014263KCH20.031291SMKB0.0267
8BUYB0.009336HAC0.022664KMSA0.023492SUCA0.0149
9CEA0.008337HAD0.035265KOJ20.024193SWO0.0255
10CEJA0.024238HALB0.025366KWJ20.011294TBA20.0352
11CHC20.016139HANB0.015067MAS20.008795TEJ20.0370
12CHJ20.011140HCNA0.024568MGY20.017796TOHA0.0183
13CHO0.025341HES0.016769MIYA0.028797ULJ20.0118
14CHR0.025442HWCA0.019270MOP0.021898USN20.0275
15CHYB0.021943HWCB0.025171MUS20.026199WJU20.0472
16CIGB0.016044ICN0.013172NAJ0.0186100YAPA0.0208
17CPR20.007745IJA20.019773NAWB0.0122101YAY0.0227
18CSO0.027746IJAA0.023574NOW−0.0065102YAYA0.0242
19CWO20.009647IKSA0.014075OKCB0.0109103YCH0.0179
20DAG20.019748IMSA0.023576OKEB0.0269104YEG0.0051
21DAU0.006149IMWB0.020577PHA20.0130105YEYB0.0190
22DGY20.005550INCA0.012578PORA0.0216106YINB0.0344
23DUSB0.012151JAHA0.028879PTK0.0113107YNCB0.0214
24EMSB0.027952JASA0.016280PUAA0.0182108YOA0.0140
25EURB0.013053JECB0.012281PYC0.0237109YOCB0.0058
26EUSB0.021254JEJB0.016282PYCA0.0199110YODB0.0222
27GAPB0.028855JEO20.002183SACA0.0155111YOJB0.0170
28GIC0.011056JES0.034084SAJ0.0227Median0.0192
Table A2. Calculated corner frequency and seismic moment at all stations used in this study.
Table A2. Calculated corner frequency and seismic moment at all stations used in this study.
No. f c , Hz M 0 , Dyne-cmNo. f c , Hz M 0 , Dyne-cmNo. f c , Hz M 0 , Dyne-cmNo. f c , Hz M 0 , Dyne-cm
10.59 1.03 × 10 25 290.70 6.37 × 10 24 570.44 6.82 × 10 24 850.76 2.53 × 10 25
20.46 5.42 × 10 25 300.22 1.92 × 10 25 580.50 2.21 × 10 25 860.48 1.46 × 10 25
30.84 8.99 × 10 24 311.42 6.82 × 10 24 590.36 1.27 × 10 25 870.25 1.36 × 10 25
40.55 5.94 × 10 24 320.66 4.21 × 10 24 600.64 2.53 × 10 25 880.63 8.39 × 10 24
51.05 4.21 × 10 24 330.49 8.99 × 10 24 610.63 1.03 × 10 25 890.39 8.99 × 10 24
60.43 1.27 × 10 25 341.46 5.94 × 10 24 620.55 8.99 × 10 24 900.35 1.27 × 10 25
70.56 4.83 × 10 24 350.58 5.94 × 10 24 631.37 4.21 × 10 24 911.17 3.93 × 10 24
80.56 5.17 × 10 24 360.44 3.12 × 10 25 640.63 5.94 × 10 24 920.53 6.37 × 10 24
90.32 1.11 × 10 25 371.43 8.39 × 10 24 650.74 5.54 × 10 24 930.89 7.83 × 10 24
100.43 4.11 × 10 25 380.72 4.83 × 10 24 660.78 4.21 × 10 24 941.46 9.64 × 10 24
110.48 4.83 × 10 24 390.49 7.31 × 10 24 671.21 4.51 × 10 24 951.29 6.82 × 10 24
120.50 9.64 × 10 24 400.80 1.79 × 10 25 680.36 8.39 × 10 24 960.59 2.72 × 10 25
131.05 4.21 × 10 24 410.84 8.99 × 10 24 690.80 4.83 × 10 24 970.28 1.46 × 10 25
141.46 4.51 × 10 24 420.58 2.72 × 10 25 700.94 5.17 × 10 24 980.97 4.83 × 10 24
150.32 1.46 × 10 25 430.68 7.83 × 10 24 711.31 8.39 × 10 24 991.37 1.56 × 10 25
160.44 5.94 × 10 24 440.45 1.03 × 10 25 721.23 5.54 × 10 24 1000.48 9.64 × 10 24
170.58 5.17 × 10 24 450.66 1.03 × 10 25 730.35 6.37 × 10 24 1010.87 8.99 × 10 24
180.74 7.83 × 10 24 460.43 6.22 × 10 25 740.32 1.27 × 10 25 1020.44 4.72 × 10 25
190.43 1.03 × 10 25 470.46 2.53 × 10 25 750.40 7.31 × 10 24 1030.80 8.39 × 10 24
200.78 5.54 × 10 24 480.66 1.11 × 10 25 760.56 7.83 × 10 24 1041.34 5.17 × 10 24
210.63 7.83 × 10 24 490.50 5.54 × 10 24 770.50 3.42 × 10 24 1050.30 1.79 × 10 25
220.23 1.67 × 10 25 500.32 3.84 × 10 25 780.63 5.94 × 10 24 1060.48 9.64 × 10 24
230.53 4.51 × 10 24 510.92 5.94 × 10 24 790.52 8.39 × 10 24 1070.50 9.64 × 10 24
240.84 6.82 × 10 24 520.61 5.94 × 10 24 800.78 1.36 × 10 25 1081.14 5.54 × 10 24
250.61 4.83 × 10 24 530.40 9.64 × 10 24 811.27 8.39 × 10 24 1090.25 7.31 × 10 24
260.45 1.11 × 10 25 540.17 1.46 × 10 25 820.61 4.72 × 10 25 1100.28 9.64 × 10 24
270.52 7.83 × 10 24 550.43 4.83 × 10 24 830.48 5.94 × 10 24 1110.39 9.64 × 10 24
281.17 5.54 × 10 24 561.14 1.11 × 10 25 841.05 7.83 × 10 24 Median0.58 8.39 × 10 24

References

  1. Johnston, A.C. Seismic Moment Assessment of Earthquakes in Stable Continental Regions-I. Instrumental Seismicity. Geophys. J. Int. 1996, 124, 381–414. [Google Scholar] [CrossRef]
  2. Schulte, S.M.; Mooney, W.D. An Updated Global Earthquake Catalogue for Stable Continental Regions: Reassessing the Correlation with Ancient Rifts. Geophys. J. Int. 2005, 161, 707–721. [Google Scholar] [CrossRef] [Green Version]
  3. Lee, K.; Yang, W.S. Historical Seismicity of Korea. Bull. Seismol. Soc. Am. 2006, 96, 846–855. [Google Scholar] [CrossRef]
  4. Wen, Y.K.; Collins, K.R.; Han, S.W.; Elwood, K.J. Dual-level designs of buildings under seismic loads. Struct. Saf. 1996, 18, 195–224. [Google Scholar] [CrossRef]
  5. Lee, L.H.; Lee, H.H.; Han, S.W. Method of selecting design earthquake ground motions for tall buildings. Struct. Des. Tall. Spec. 2000, 9, 201–213. [Google Scholar] [CrossRef]
  6. Han, S.W.; Seok, S.W. Efficient Procedure for Selecting and Scaling Ground Motions for Response History Analysis. J. Struct. Eng. 2014, 140, 06013004. [Google Scholar] [CrossRef]
  7. Saragoni, G.; Hart, G. Simulation of Artificial Earthquakes. Earthq. Eng. Struct. Dyn. 1974, 2, 249–267. [Google Scholar] [CrossRef]
  8. Hanks, T.C.; McGuire, R.K. The Character of High-frequency Strong Ground Motion. Bull. Seismol. Soc. Am. 1981, 71, 2071–2095. [Google Scholar]
  9. Boore, D.M. Stochastic Simulation of High-frequency Ground Motions Based on Seismological Models of the Radiated Spectra. Bull. Seismol. Soc. Am. 1983, 73, 1865–1894. [Google Scholar]
  10. Boore, D.M.; Atkinson, G.M. Stochastic Prediction of Ground Motion and Spectral Response Parameters at Hard-rock Sites in Eastern North America. Bull. Seismol. Soc. Am. 1987, 77, 440–467. [Google Scholar]
  11. Atkinson, G.M.; Boore, D.M. Ground-motion Relations for Eastern North America. Bull. Seismol. Soc. Am. 1995, 85, 17–30. [Google Scholar]
  12. Boore, D.M. Simulation of Ground Motion Using the Stochastic Method. Pure. Appl. Geophys. 2003, 160, 635–676. [Google Scholar] [CrossRef] [Green Version]
  13. Papazafeiropoulos, G.; Plevris, V. OpenSeismoMatlab: A New Open-source Software for Strong Ground Motion Data Processing. Heliyon 2018, 4, e00784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Nakamura, Y. A Method for Dynamic Characteristics Estimation of Subsurface Using Microtremor on the Ground Surface. Q. Rep. Railw. Tech. Res. Inst. 1989, 30, 25–33. [Google Scholar]
  15. Atkinson, G.M.; Cassidy, J.F. Integrated Use of Seismograph and Strong-Motion Data to Determine Soil Amplification: Response of the Fraser River Delta to the Duvall and Georgia Strait Earthquakes. Bull. Seismol. Soc. Am. 2000, 90, 1028–1040. [Google Scholar] [CrossRef]
  16. Zhao, J.X.; Irikura, K.; Zhang, J.; Yoshimitsu, F.; Somerville, P.G.; Asano, A.; Ohno, Y.; Oouchi, T.; Takahashi, T.; Ogawa, H. An Empirical Site-Classification Method for Strong-Motion Stations in Japan Using H/V Response Spectral Ratio. Bull. Seismol. Soc. Am. 2006, 96, 914–925. [Google Scholar] [CrossRef]
  17. Zandieh, A.; Pezeshk, S. Investigation of Geometrical Spreading and Quality Factor Functions in the New Madrid Seismic Zone. Bull. Seismol. Soc. Am. 2010, 100, 2185–2195. [Google Scholar] [CrossRef]
  18. Brune, J. Tectonic Stress and the Spectra of Seismic Shear Waves. J. Geophys. Res. 1970, 75, 4997–5009. [Google Scholar] [CrossRef] [Green Version]
  19. Brune, J. Correction: Tectonic Stress and the Spectra of Seismic Shear Waves. J. Geophys. Res. 1971, 76, 5002. [Google Scholar]
  20. Chandramohan, R.; Baker, J.W.; Deierlein, G.G. Quantifying the Influence of Ground Motion Duration on Structural Collapse Capacity Using Spectrally Equivalent Records. Earthq. Spectra 2016, 32, 927–950. [Google Scholar] [CrossRef] [Green Version]
  21. Konno, K.; Ohmachi, T. Ground-Motion Characteristics Estimated from Spectral Ratio between Horizontal and Vertical Components of Microtremor. Bull. Seismol. Soc. Am. 1998, 88, 228–241. [Google Scholar]
  22. Boore, D.M.; Boatwright, J. Average Body-wave Radiation Coefficients. Bull. Seismol. Soc. Am. 1984, 74, 1615–1621. [Google Scholar]
  23. Kim, S.K. A Study on the Crustal Structure of the Korean Peninsula. J. Geol. Soc. Korea 1995, 8, 393–403. [Google Scholar]
  24. Jee, H.W.; Han, S.W. Equations of Path Effects to Simulate Ground Motions in Korean Peninsula Using Point-source Model. In Proceedings of the The 2019 World Congress on Advances in Structural Engineering and Mechanics (ASEM19), Jeju, Korea, 17–21 September 2019. [Google Scholar]
  25. Anderson, J.G.; Hough, S.E. A Model for the Shape of the Fourier Amplitude Spectrum of Acceleration at High Frequencies. Bull. Seismol. Soc. Am. 1984, 74, 1969–1993. [Google Scholar]
  26. Atkinson, G.M.; Boore, D.M. Evaluation of Models for Earthquake Source Spectra in Eastern North America. Bull. Seismol. Soc. Am. 1998, 88, 917–934. [Google Scholar]
  27. Jo, N.D.; Baag, C.E. Stochastic Prediction of Strong Ground Motions in Southern Korea. J. Earthq. Eng. Soc. Korea 2001, 5, 17–26. [Google Scholar]
  28. Joshi, A.; Tomer, M.; Lal, S.; Chopra, S.; Singh, S.; Prajapati, S.; Sharma, M.L.; Sandeep. Estimation of the Source Parameters of the Nepal Earthquake from Strong Motion Data. Nat. Hazards. 2016, 83, 867–883. [Google Scholar] [CrossRef]
  29. Andrews, D.J. Objective Determination of Source Parameters and Similarity of Earthquakes of Different Size. Earthq. Source Mech. 1986, 37, 259–267. [Google Scholar]
  30. Boatwright, J.; Choy, G.L. Acceleration source spectra anticipated for large earthquakes in northeastern North America. Bull. Seismol. Soc. Am. 1992, 82, 660–682. [Google Scholar]
Figure 1. The 2017 Pohang earthquake and seismic recording stations with ray path lines.
Figure 1. The 2017 Pohang earthquake and seismic recording stations with ray path lines.
Applsci 10 01254 g001
Figure 2. Original and adjusted ground motions at station BUS2: (a) Original data-acceleration, (b) Original data-velocity, (c) Original data-displacement, (d) Low cut filter-acceleration, (e) Low cut filter-velocity, (f) Low cut filter-displacement, (g) Baseline-acceleration, (h) Baseline-velocity, (i) Baseline-displacement, (j) Signal processing-acceleration, (k) Signal processing-velocity, (l) Signal processing-displacement.
Figure 2. Original and adjusted ground motions at station BUS2: (a) Original data-acceleration, (b) Original data-velocity, (c) Original data-displacement, (d) Low cut filter-acceleration, (e) Low cut filter-velocity, (f) Low cut filter-displacement, (g) Baseline-acceleration, (h) Baseline-velocity, (i) Baseline-displacement, (j) Signal processing-acceleration, (k) Signal processing-velocity, (l) Signal processing-displacement.
Applsci 10 01254 g002
Figure 3. Site amplification factor function from the 2017 Pohang earthquake at sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Figure 3. Site amplification factor function from the 2017 Pohang earthquake at sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Applsci 10 01254 g003
Figure 4. Ground motion recorded at station BUS2 with and without site effects: (a) Ground motion acceleration with site effects, (b) Ground motion acceleration without site effects, (c) Fourier amplitude spectrum with site effects, (d) Fourier amplitude spectrum without site effects.
Figure 4. Ground motion recorded at station BUS2 with and without site effects: (a) Ground motion acceleration with site effects, (b) Ground motion acceleration without site effects, (c) Fourier amplitude spectrum with site effects, (d) Fourier amplitude spectrum without site effects.
Applsci 10 01254 g004
Figure 5. Original and smoothed Fourier amplitude spectrum: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Figure 5. Original and smoothed Fourier amplitude spectrum: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Applsci 10 01254 g005
Figure 6. A ( f ) , A ( f ) , and κ 0 values of three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Figure 6. A ( f ) , A ( f ) , and κ 0 values of three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Applsci 10 01254 g006
Figure 7. A ( f ) versus A ( f ) for three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Figure 7. A ( f ) versus A ( f ) for three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Applsci 10 01254 g007
Figure 8. Distributions of the values of f c and M 0 , according to hypocentral distance: (a) Measured corner frequencies at 111 stations, (b) Measured seismic moments at 111 stations.
Figure 8. Distributions of the values of f c and M 0 , according to hypocentral distance: (a) Measured corner frequencies at 111 stations, (b) Measured seismic moments at 111 stations.
Applsci 10 01254 g008
Figure 9. Measured and calculated Fourier amplitude spectrum (FAS) at three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Figure 9. Measured and calculated Fourier amplitude spectrum (FAS) at three sample stations: (a) BUS2 station, (b) CEA station, (c) CEJA station.
Applsci 10 01254 g009
Figure 10. Estimated shaping window and duration at the BUS2 station for the 2017 Pohang earthquake: (a) Ground motion acceleration, (b) Normalized arias intensity, (c) Shaping window.
Figure 10. Estimated shaping window and duration at the BUS2 station for the 2017 Pohang earthquake: (a) Ground motion acceleration, (b) Normalized arias intensity, (c) Shaping window.
Applsci 10 01254 g010
Figure 11. Estimated shaping window and duration for the 2017 Pohang earthquake: (a) Shaping window, (b) Duration with considering corner frequency.
Figure 11. Estimated shaping window and duration for the 2017 Pohang earthquake: (a) Shaping window, (b) Duration with considering corner frequency.
Applsci 10 01254 g011
Figure 12. Recorded and simulated ground motion at sample stations: (a) BUS2 station-EW direction, (b) BUS2 station-NS direction, (c) BUS2 station-simulated, (d) BUS2 station-FAS, (e) BUS2 station-PSA, (f) CEA station-EW direction, (g) CEA station-NS direction, (h) CEA station-simulated, (i) CEA station-FAS, (j) CEA station-PSA, (k) CEJA station-EW direction, (l) CEJA station-NS direction, (m) CEJA station-simulated, (n) CEJA station-FAS, (o) CEJA station-PSA.
Figure 12. Recorded and simulated ground motion at sample stations: (a) BUS2 station-EW direction, (b) BUS2 station-NS direction, (c) BUS2 station-simulated, (d) BUS2 station-FAS, (e) BUS2 station-PSA, (f) CEA station-EW direction, (g) CEA station-NS direction, (h) CEA station-simulated, (i) CEA station-FAS, (j) CEA station-PSA, (k) CEJA station-EW direction, (l) CEJA station-NS direction, (m) CEJA station-simulated, (n) CEJA station-FAS, (o) CEJA station-PSA.
Applsci 10 01254 g012
Figure 13. Estimated residuals, mean residual, and slope residual at various periods for the 2017 Pohang earthquake: (a) Path effect function used in this study, (b) Path effect function by Boore (2003).
Figure 13. Estimated residuals, mean residual, and slope residual at various periods for the 2017 Pohang earthquake: (a) Path effect function used in this study, (b) Path effect function by Boore (2003).
Applsci 10 01254 g013
Table 1. Information on the 2017 Pohang earthquake (National Earthquake Comprehensive Information System (NECIS)).
Table 1. Information on the 2017 Pohang earthquake (National Earthquake Comprehensive Information System (NECIS)).
Event NameLocal Date-TimeLongitude (East)Latitude (North)Focal Depth (km)Magnitude M L
the 2017 Pohang earthquake14:29 15 November, 2017129.3736.1195.4
Table 2. Information about the functions and component parameters of the stochastic point-source model.
Table 2. Information about the functions and component parameters of the stochastic point-source model.
FunctionsParameters
Source effect function S o u r c e ( M 0 , f , f c ) = M o 1 + ( f / f c ) 2 × R θ ϕ F V 4 π ρ β S 3 × ( 2 π f ) p × 1 R r e f R θ ϕ ( = 0.63 ) : S-wave averaged radiation pattern coefficient [22]
F ( = 2 ) : free surface effect [9]
V ( = 1 / 2 ) : partition coefficient of a vector into the horizontal component [9]
ρ ( = 2.7   g / cm 3 ) : near source soil density [23]
β S ( = 3.36   km / s ) : near source shear wave velocity [23]
p : ground motion type coefficient (0, 1, and 2 for displacement, velocity, and acceleration, respectively) [9]
R r e f ( = 1   km ) : reference source-to-site distance for seismic source
Path effect function P a t h ( R H , f ) = G ( R H ) × exp ( π f R H / Q ( f ) β s ) G ( R H ) ( = { R H 1.3 ( R H 70   km ) 70 1.3 ( R H / 70 ) 0.3 ( 70   km < R H 100   km ) 70 1.3 ( 100 / 70 ) 0.3 ( R H / 100 ) 0.5 ( R H > 100   km ) )
: geometrical attenuation function [24]
Q ( f ) ( = 348 f 0.48 ) : quality factor of the anelastic attenuation function [24]
Site effect function S i t e ( κ 0 , f ) = A M P ( f ) × exp ( π κ 0 f ) A M P ( f ) : site amplification factor function at each station
κ 0 : site attenuation coefficient of the site attenuation function [25]

Share and Cite

MDPI and ACS Style

Han, S.W.; Jee, H.W. A Numerical Model for Simulating Ground Motions for the Korean Peninsula. Appl. Sci. 2020, 10, 1254. https://doi.org/10.3390/app10041254

AMA Style

Han SW, Jee HW. A Numerical Model for Simulating Ground Motions for the Korean Peninsula. Applied Sciences. 2020; 10(4):1254. https://doi.org/10.3390/app10041254

Chicago/Turabian Style

Han, Sang Whan, and Hyun Woo Jee. 2020. "A Numerical Model for Simulating Ground Motions for the Korean Peninsula" Applied Sciences 10, no. 4: 1254. https://doi.org/10.3390/app10041254

APA Style

Han, S. W., & Jee, H. W. (2020). A Numerical Model for Simulating Ground Motions for the Korean Peninsula. Applied Sciences, 10(4), 1254. https://doi.org/10.3390/app10041254

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop