| Citation: |
Lan-chi Kang, Xing Jin, Yong-xiang Wei (2013). Modifying the intensity distribution by assessing the reliability. Earthq Sci 26(6): 407-413. DOI: 10.1007/s11589-013-0058-4
|
This article presents an application of a procedure to modify the intensity distribution by assessing the reliability. There are two potential possibilities that may influence the intensity distribution: (1) For the interpolation error, we generate a measured grid across the calculation region. When the point to station spacing is < 5 km, we consider the results precise; however, some points have less precision because these are farther from the corresponding stations. When the spacing is between 5 and 50 km, we consider the results imprecise and define a reliability factor that correlates with the distance. (2) Some records may have errors that result from local site conditions, equipment problems, or some disturbance such as lightning stroke, which will lead to some grid points having an incorrect intensity. We regress the attenuation relation for sites with abnormal intensities and consider the results to be accurate when the standard deviation (STD) is < σ and inaccurate when the STD is > 2σ. We then define a reliability factor to correlate with STD between σ and 2σ, such that the intensity distribution is in accord with both wave propagation theory and the investigation intensity.
For rapid emergency response, such as resource dispatch management, information regarding the spatial distribution of seismic intensity is particularly useful, and a near real-time or quick damage assessment in terms of seismic intensity is greatly in demand. One practical way to fulfill this demand is to correlate seismic intensity with peak ground motions because they can be acquired in near real time via a rapid reporting system (Lind et al. 1996; Yamazaki et al. 1997; David et al. 1999; Kanamori 1993; Ma et al. 2003; Jin et al. 2008; Wei et al. 2010).
On May 12, 2008, the large magnitude Wenchuan earthquake struck the densely populated Sichuan Province of China (31.0°N, 103.4°E). More than 80, 000 people were killed, ~370, 000 people were injured, and economic losses have been estimated to be in excess of 800 billion RMB. The large moment magnitude (MW 7.9) (Hao et al. 2009; Shen Sun et al. 2009; Xu et al. 2009; Zhang et al. 2009) of the event attests to it being a great intercontinental earthquake. The Wenchuan earthquake occurred on the Longmen Shan fault zone, which delineates the eastern margin of the Tibet Plateau. Located on the western circum-Pacific seismic belt, Taiwan is situated in the collision zone between the Philippine Sea plate (PSP) and the Eurasian plate (EP). As a result, seismicity in Taiwan is considerably high, as demonstrated by the catastrophic 1999 Jiji (Chi-Chi) earthquake. For rapid reporting of felt earthquakes within Taiwan, a real-time strong motion network, the Taiwan Rapid Earthquake Information Release System (TREIRS), has been operated by the Central Weather Bureau (CWB) since 1997 (Shin and Teng 2001; Wang et al. 2005). In addition to routine magnitude and hypocenter determinations following a felt earthquake, peak ground motions recorded by the TREIRS system are used to generate seismic intensity.
In this paper, we use the records from the Wenchuan and Jiji (Chi-Chi) earthquakes to simulate the progress of rapid reporting of seismic intensity. After mesh generation in the calculation region, we apply the interpolation method to calculate the peak ground acceleration (PGA) for each grid. In accordance with the intensity scale, we transform the PGA to intensity. However, owing to some problems, following the analysis of the intensity distribution, we present a reliability index to modify the intensity distribution.
In general, PGA is calculated in the reporting system of seismic intensity; however, estimating the intensity from PGA is a complex academic problem. Therefore, in this paper, we analyze the relation between intensity and PGA, which are given on the Chinese seismic intensity scale (GB/T 17742-2008). The relation can be written as
|
|
(1) |
where Am is the PGA, I is the intensity, and a and b are the regressed coefficients.
Figure 1 shows the relation between Am and I. In the areas below I = IV, PGAs corresponding to intensity are not given, and the regression is shown as a dash-dot line. Based on this result, we can directly calculate the intensity from PGA.
In this study, we divide the calculated regions into 16, 416 grids with a grid spacing of 0.02° and 21, 091 grids with a grid spacing of 0.05° for the Jiji (Chi-Chi) and Wenchuan earthquakes, respectively. From the catalog, 379 records distributed acro ss the whole of Taiwan have been selected for the Chi-Chi earthquake (Fig. 2b), and 87 records have been selected from the Wenchuan earthquake, including all stations in Sichuan Province and additional stations in Gansu and Shanxi provinces (Fig. 2a).
Peak ground acceleration and peak ground velocity (PGV) grid points are calculated by an interpolation methodology as shown in Fig. 3 (Wei et al. 2010; Zhang 2008). In this interpolation methodology, the weighted is refined as the inverse square of the distance, as follows:
|
|
(2) |
The methodology can be written as follows:
|
|
(3) |
where Wa, Wb, Wc is the weighted, Aa, Ab, Ac is the PGA of a, b, and c stations, and f (R, Ra), f (R, Rb), f (R, Rc) is the coefficient of attenuation relations.
By transforming PGA or PGV at each grid point into intensity, we can generate an intensity contour (Fig. 4).
Figure 4 depicts the spatial distribution of intensity for the calculated region. In general, the spatial distribution of intensity is closely related to station distribution, the quality of the records, or site conditions. Two notable abnormal regions are observed in the intensity distribution. This may be due to a limited distribution of stations in these regions, so we could interpret that the interpolation precision is related to the distance from the stations to the grid points. Alternatively, local site conditions may influence PGA. Some grid points are controlled by abnormal PGA, for example, Fig. 5 showing three records, where the PGAs are higher contrast to some records with a closer epicentral distance; therefore, the region must be showing abnormal intensities. These two reasons lead to abnormal intensity estimates; therefore, we present a methodology to modify the intensity distribution by assessing the reliability.
The most commonly used attenuation model in China that is recommended for Evaluation of Seismic Safety for Engineering Sites (GB17741-2005) is
|
|
(4) |
where Y is the peak ground motion or response spectrum at a variable period; Ci (i = 1, 2, …7) is the regression coefficient; M is the local magnitude; R is the fault distance, which can be written as D; R0 (M) = C5 exp (C6M) is a saturation factor that has a correlation with the earthquake intensity. In general, R0 (M) is the hypocentral depth and σ is the random variable. The C2M+C3M2 is a constant if the attenuation is response to particular earthquake. In general, the coefficient C7R is very small, which shows that the coefficient C7 has less relation to Y; therefore, we remove C7R, so that Eq. (3) can be written as
|
|
(5) |
Figure 6 and Table 1 shows the regression results and the standard deviation (STD) of PGA. We generally accept the PGA when the STD < σ and assume that stations with STD > 2σ are not valid; thus, the PGA reliability can be refined as
 |
DownLoad:
CSV
|
|
(6) |
where η is residuals between observation and the attenuation relation, σ is the STD, P(Am) is the PGA reliability, and the modifying methodology can be written as:
|
|
(7) |
where Am is the observation PGA, A is the calculated PGA by the attenuation relation.
We use an interpolation method of inverse weighted square distance, presented by Jin and Zhang, where the interpolation precision is related to the interpolation distance. This method is very accurate when the distance from point to the station is < 5 km (Rmin). However, if the distance from point to the station is > 50 km (Rmax), the method is inaccurate and these grid points will be calculated by the attenuation relation, and we define a reliability index with which the distance correlates, when the distance is between 5 and 50 km. The reliability P (R) of interpolation distance can be refined as:
|
|
(8) |
The modifying intensity distribution can be written as:
|
|
(9) |
where Am' is the modified PGA, Am is the unmodified PGA, and A is calculated by the attenuation.
The synthetical reliability can be refined as:
|
|
(10) |
Figure 7 depicts the spatial distributions of the synthetic reliability; we find that the error regions of previous intensity distribution have a very low reliability, so we confirm that this method could be used to modify the errors. The modifying method can be written as:
|
|
(11) |
where Am' is the modified PGA, Am is the unmodified PGA, and A is calculated by the attenuation Ps is the synthetic reliability.
Figure 8 shows the modified intensity distribution, which can be differentiated from Fig. 3 by the absence of regions with abnormal intensities, and the grid points in these regions are calculated by the attenuation relation. This methodology mainly modifies the points with greater separation between station and points, with the aid of which the results accord not only with the theory of the wave propagation, but also with the investigation intensity.
We analyze the intensity distribution of the Wenchuan and Jiji (Chi-Chi) earthquakes. We found that some interpolation errors and abnormal stations influence the intensity distribution. First, we find that some stations with an abnormal PGA have affects on the intensity distribution and introduce the reliability of PGA, refining the PGA reliability as STD, which are calculated by the attenuation relation. We then introduce the interpolation reliability; however, some researchers suggest that the interpolation methodology has only high precision if the distance is less than 5 km, but gives large errors when the distances are greater than 50 km. So we adopt the PGA attenuation relation to modify some error grid points in the intensity distribution. The modified results are in accord with wave propagation theory and also meet the investigation intensity. All of these tasks can be completed in 2 min; thus, this methodology can be applied to the rapidly reporting system of seismic intensity.
This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author (s) and the source are credited.
|
David JW, Quitoriano V, Heaton TH, Kanamori H, Scrivner CW, Worden CB (1999) TriNet "ShakeMaps": Rapid generation of peak ground motion and intensity maps for earthquake in Southern California. Earthq Spectra 15:537-556 doi: 10.1193/1.1586057
|
|
General administration of quality supervision, inspection and quarantine of the People's Republic of China (2008) The Chinese seismic intensity scale (GB/T 17742-2008). China Zhijian Publishing House, Beijing
|
|
Hao KX, Si HJ, Fujiwara H, Ozawa T (2009) Coseismic surface ruptures and crustal deformations of the 2008 Wenchuan earthquake MW 7.9, China. Geophys Res Lett 36:L11303 doi: 10.1029/2009GL037971
|
|
Jin X, Zhang HC, Wei YX (2008) Shaking Intensity Criterion Based on Rapid Release Information by Seismic Monitoring Network and Its Application. Recent Dev World Seismol 10:20-27 (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTotal-GJZT200810006.htm
|
|
Kanamori H (1993) Location earthquake with amplitude application to real-time seismology. Bull Seismol Soc Am 83:264-268
|
|
Lind SG, Douglas SN, Douglas SD, Michael EP (1996) Real-time seismology at UC Berkeley: the rapid earthquake integration project. Bull Seismol Soc Am 86:936-945
|
|
Ma Q, Jin X, Li SY (2003) Real-time algorithm of seismic dynamic response for SDOF system. Earthq Eng Eng Vib 23:61-68 (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTOTAL-dggc200305011.htm
|
|
Shen Sun Z, Zhang J, Wan P, Wang Y, Bürgmann M, Zeng R, Gan YW, Liao H, Wang Q (2009) Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake. Nat Geosci 2(10):718-724 doi: 10.1038/ngeo636
|
|
Shin TC, Teng TL (2001) An overview of the 1999 Chi-Chi Taiwan earthquake. Bull Seismol Soc Am 91:895-913
|
|
Wang WM, Zhao LF, Li J, Yao ZX (2005) Rupture process of the Chichi (Taiwan) earthquake in 1999. Chin J Geophys 48:132-147 (in Chinese with English abstract)
|
|
Wei YX, Li J, Jin X, Zhang HC (2010) Preliminary study of near real-time interpolation calculation of ShakeMap based on seismic monitoring networks information. J Disaster Prev Reduct 26:2-11 (in Chinese with English abstract) http://en.cnki.com.cn/Article_en/CJFDTotal-DDYJ201001002.htm
|
|
Xu X, Wen X, Yu G, Chen G, Klinger Hubbard YJ, Shaw J (2009) Coseismic reverse-and oblique-slip surface faulting generated by the 2008 MW 7.9 Wenchuan earthquake, China. Geology 37:515-518 doi: 10.1130/G25462A.1
|
|
Yamazaki F, Noda S, Meguro K (1997) Developments of monitoring and early damage assessment system in Japan. Bull ERS 30:15-31
|
|
Zhang HC (2008) Research on ShakeMaps and Shaking Intensity Based on Records for Monitoring Network. Master Dissertation, Harbin: Institute of Engineering Mechanics State Seismological Bureau, China (in Chinese)
|
|
Zhang Y, Feng WP, Xu LS, Zhou CH, Chen YT (2009) Spatio-temporal rupture process of the 2008 great Wenchuan earthquake. Sci China Earth Sci 52:145-154 doi: 10.1007/s11430-008-0148-7
|
震球期刊
EQS
Supported by: Beijing Renhe Information Technology Co., Ltd. 
