Abstract
The prediction of earthquake ground motion is the first priority in the seismic design of a building. This investigation is aimed at proposing regression models for average peak horizontal ground acceleration of earthquakes recorded by seismometers installed at a station in Chiang Mai, Thailand. The majority of earthquakes measured in Chiang Mai occur in seven areas; namely, the regions around Sumatra, Nicobar Island, the Andaman Sea, Myanmar, Laos, China’s western region, and China’s southern region. The earthquakes’ epicenters range from about 10 to 2600 km away from Chiang Mai. The proposed model used 73 earthquakes recorded from 2006 to 2012 and was subdivided according to the magnitudes of the events and earthquake source zones. It was found that the average peak horizontal ground acceleration by regression models was attenuated by the distance from the epicenters. The results of the regression model were in agreement with the records of seven recent earthquakes obtained from Chiang Mai’s seismic station. The regression model has been used in the design of buildings.
Keywords
- prediction
- earthquake
- ground motion
- regression model
- average peak horizontal ground acceleration
1. Introduction
Southeast Asia was shaped as a result of the interaction between the Indo-Australian, Eurasian, Philippine, and West Pacific tectonic plates. Thailand located inside the Eurasian plate, is delimited by the Andaman thrust within the west, the Sunda Arc within the south, and also the Philippine trench within the east [1]. As a result of recent awareness of earthquake hazards, the Thai Meteorological Department has established quite 20 seismometer stations Thailand. Figure 1 shows samples of ground motion records at a Chiang Mai station (in northern Thailand) from an earthquake event in the regions around Chiang Rai (M 5.9). The distance from the epicenter of the earthquake is about 128 km. One will observe a small peak ground acceleration (PGA) (≈2.29 × 10−3 g) in the vertical component, and a peak horizontal ground acceleration in the horizontal components (E-W component ≈ 1.63 × 10−3 g and N-S component ≈ 1.19 × 10−3 g). The researcher calculated average peak horizontal ground acceleration (PHAavg) ≈ 1.04 × 10−3 g.

Figure 1.
Earthquake sampling E-W components of the Chiang Mai station (UTC: 5/6/2014 00:50:16, Lat = 19.73 N, Lon = 99.69E, Mw = 5. 9, distance = 128 km, and epicenter: Chiang Rai) (Thai Meteorological Department).
The attenuation relationship for PHAavg for the Chiang Mai station in Thailand has been proposed in this investigation. The attenuation relationship for PHAavg was determined based on the multiple linear regression (MLR) models by using 132 components of strong motion data from 66 earthquake events. By using the MLR models, the researcher constructed site-specific attenuation relation based on a previously proposed attenuation relation [2, 3]. The specific site in this investigation is a Chiang Mai seismic station (CMMT) station located on the rock site.
2. Characteristics of ground motion in Thailand
Because Thailand is not a country prone to earthquakes, most of the earthquakes felt in Thailand occur outside the country, in places such as Sumatra, Nicobar Island, the Andaman Sea, Myanmar, Laos, and in the west and south areas of China. However, inside the country, Thailand has active faults, which can cause small-to-moderate earthquakes in areas such as Chiang Mai.
Pairojn and Wasinrat [4] presented the observed earthquakes recorded by the Chiang Mai station in northern Thailand from 07/10/2006 to 11/11/2012, as shown in Figure 2.

Figure 2.
The location of (a) 41 seismic stations in Thailand, and (b) 73 earthquake events used in this investigation (Thai Meteorological Department).
3. Site classification
In this study, the researcher chose to investigate data from shear wave velocity (Vs) recorded at the Chiang Mai seismic station (CMMT) in northern Thailand, located on a rocky site (NEHRP site class). The CMMT station used Trillium 120 velocity sensors and TSA100S accelerometers, as seen in Figure 3.

Figure 3.
Installation of velocity sensor and accelerometer at CMMT seismic station (Thai Meteorological Department) [5].
4. Ground motion prediction
Recently, several ground-motion prediction equations (GMPEs) have been developed. Douglas [6] summarizes all empirical ground-motion prediction equations (GMPEs) used to estimate earthquake peak ground acceleration (PGA) and elastic response spectral ordinates published between 1964 and 2017. Most empirical ground motion attenuation relations are derived from numerical analyses.
For Thailand, Idriss’s model [2, 3] was selected because it has an appropriate attenuation relation. Idriss suggested that the attenuation relation for motions in western North America [2, 3]. The available data included 572 individual horizontal components in rock sites that were used to derive this attenuation model. This study focuses on M ≤ 6.0 and M > 6.0 using local ML and surface wave magnitude MS scales, respectively. The range of applicability is 1–100 km for distance and 4.6–7.4 for M [2, 3]. The peak ground acceleration (Y) at rock sites was derived as following equation:
where Y is in g (m/s2), a = 0.2, for M ≤ 6 𝛼0 = −0.150, 𝛼1 = 2.261, 𝛼2 = −0.083, 𝛽0 = 0, 𝛽1 = 1.602, 𝛽2 = −0.142 and 𝜎 = 1.39–0.14 M and for M > 6 𝛼0 = −0.050, 𝛼1 = 3.477, 𝛼2 = −0.284, 𝛽0 = 0, 𝛽1 = 2.475, 𝛽2 = −0.286 and for M < 7.25 𝜎 = 1.39–0.14 M, and for M ≥ 7.25 𝜎 = 0.38. (F = 0 for Strike slip, F = 0.5 for Oblique, F = 1 for Reverse).
This chapter uses the method of ground-motion prediction under regression analysis [5].
4.1. Concept of multiple linear regression models
Regression analysis is a conceptually simple method for investigating functional relationships among variables. The relationship is explained by an equation or a model combining the response (Y) and predictor variables (X1, X2, …, Xp) [7]. This relationship can be derived from the regression model.
Where
while an example of a non-linear function is:
Note that the term linear here does not describe the relationship between Y and X. This model is linear because in each case the parameters enter linearly, although the relationship between Y and X is non-linear. This can be transformed as follows:
where in the equation we have X1 = InX. The variables here are transformed [7].
A regression model containing only one predictor variable is called a simple linear regression. The model containing more than one predictor variable is called a multiple linear regression (MLR). The multiple linear regression model is:
where
The value of
4.2. The multiple linear regression for predicting earthquake ground motion
Pairojn and Wasinrat [4] proposed the MLR model to predict average peak horizontal ground acceleration, PHAavg (
where M represents the magnitude and R represents the distance. The natural log (ln) of Y and R are used to transform the linear relationship between Y and R, and lnY and M. The value of
No. | Date | Time UTC | Epicenter | Magnitude, M (Mw) | Distance, R (km) | PHAavg, Y (g) |
---|---|---|---|---|---|---|
1 | 5/5/2014 | 12:06:19 | Chiang Rai | 5.1 | 121 | 0.0002990 |
2 | 5/5/2014 | 12:20:57 | Chiang Rai | 5.2 | 140 | 0.0003755 |
3 | 5/5/2014 | 21:17:05 | Chiang Rai | 5.1 | 119 | 0.0005555 |
4 | 5/5/2014 | 23:04:55 | Chiang Rai | 5.2 | 121 | 0.0003380 |
5 | 6/5/2014 | 00:50:16 | Chiang Rai | 5.9 | 128 | 0.0014140 |
6 | 6/5/2014 | 0:58:19 | Chiang Rai | 5.6 | 116 | 0.0010405 |
7 | 12/5/2014 | 11:05:29 | Chiang Rai | 5.0 | 137 | 0.0002955 |
Table 1.
Seven new earthquake events.

Figure 4.
The comparison of MLR and Idriss’s models at CMMT seismic station [4].
5. The practicality of ground motion prediction
The MLR model was used on seven new earthquake events that occurred from 05/05/2014 to 05/12/2014. The data are shown in Table 1.
The predicted PHA values from Eq. (8) and Eq. (1) are presented in Table 2 (Figures 5, 6, 7).
Model | Predicted PHA values | ||||||
---|---|---|---|---|---|---|---|
No. 1 | No. 2 | No. 3 | No. 4 | No. 5 | No. 6 | No. 7 | |
MLR | 0.000410 | 0.000346 | 0.000426 | 0.000488 | 0.001446 | 0.001080 | 0.000257 |
Idriss | 0.002751 | 0.002751 | 0.002751 | 0.002751 | 0.002751 | 0.002751 | 0.002751 |
Table 2.
The predicted peak horizontal ground acceleration from MLR model and Idriss’s model.

Figure 5.
The predicted peak horizontal ground acceleration from MLR model.

Figure 6.
The predicted peak horizontal ground acceleration from Idriss’s model.

Figure 7.
A comparison of MLR and Idriss’s model for new earthquake events.
The root mean square error (RMSE) is used for comparing the predicted PHA values from MLR model to the observed PHA values. The RMSE measures the error between a predicted value and an observed value, defined as:
where
6. Conclusions
The MLR model is presented for developing a previous attenuation relationship based on observations. The MLR model is suitable for measuring the attenuation relationship for Thailand because Thailand has few instances of motion data, with most peak ground acceleration being measured at less than 0.1 g. It is expected that this method can be applied to observation sites throughout the country. The MLR model has been used for probabilistic hazard analysis, risk analysis, building design analysis, and many other fields, that is, and construction of nuclear power plants, dams, bridges, and high-rise buildings.
Acknowledgments
The authors would like to acknowledge the Chandrakasem Rajabhat University and Thai Meteorological Department, which collectively supported this project.