Probabilistic seismic hazards assessments of Sabah ...

5 downloads 0 Views 3MB Size Report
Sep 27, 2017 - Programme by the government separates Sabah into 6 strategic development areas with each area. 44 focusing on a particular aspect of the ...
Journal of Geophysics and Engineering

ACCEPTED MANUSCRIPT

Probabilistic seismic hazards assessments of Sabah, east Malaysia: accounting for local earthquake activity near Ranau To cite this article before publication: Amin Esmail Khalil et al 2017 J. Geophys. Eng. in press https://doi.org/10.1088/1742-2140/aa8d51

Manuscript version: Accepted Manuscript Accepted Manuscript is “the version of the article accepted for publication including all changes made as a result of the peer review process, and which may also include the addition to the article by IOP Publishing of a header, an article ID, a cover sheet and/or an ‘Accepted Manuscript’ watermark, but excluding any other editing, typesetting or other changes made by IOP Publishing and/or its licensors” This Accepted Manuscript is © 2017 Sinopec Geophysical Research Institute.

During the embargo period (the 12 month period from the publication of the Version of Record of this article), the Accepted Manuscript is fully protected by copyright and cannot be reused or reposted elsewhere. As the Version of Record of this article is going to be / has been published on a subscription basis, this Accepted Manuscript is available for reuse under a CC BY-NC-ND 3.0 licence after the 12 month embargo period. After the embargo period, everyone is permitted to use copy and redistribute this article for non-commercial purposes only, provided that they adhere to all the terms of the licence https://creativecommons.org/licences/by-nc-nd/3.0 Although reasonable endeavours have been taken to obtain all necessary permissions from third parties to include their copyrighted content within this article, their full citation and copyright line may not be present in this Accepted Manuscript version. Before using any content from this article, please refer to the Version of Record on IOPscience once published for full citation and copyright details, as permissions will likely be required. All third party content is fully copyright protected, unless specifically stated otherwise in the figure caption in the Version of Record. View the article online for updates and enhancements.

This content was downloaded from IP address 202.170.60.244 on 27/09/2017 at 06:53

Page 1 of 29

pt

1

Probabilistic Seismic Hazards Assessments of Sabah, east Malaysia:

3

Accounting for local earthquake activity near Ranau

4

By

5

Amin E. Khalil1,2, Ismail A. Abir1, Hesham E. Abdel Hafiez3, Hanteh Ginsos3 and Sohail Khan4

6

1

7

2

8

3

us cri

2

School of Physics, Univesiti Sains Malaysia, Pulau Pinang, Malaysia, 11800.

Geology Department, Faculty of Science, Helwan University, Helwan, Cairo, Egypt, 11795 National Research Institute of Astronomy and Geophysics, Helwan, Cairo, Egypt, 11421

9

Correspondent Author: Amin E. Khalil email: [email protected]

11

Abstract

an

10

Sabah state in eastern Malaysia, unlike most of the other Malaysian states, is

13

characterized by a common seismological activity; generally a moderate magnitude

14

earthquake is experienced at a roughly 20 years’ interval originating mainly from two

15

major sources, either local source (e.g. Ranau and Lahad Dato) or regional source (e.g.

16

Kalimantan and South Philippines subductions). The seismicity map of Sabah shows

17

the presence of two zones of distinctive seismicity, where these zones are near Ranau

18

(near Kota Kinabalo) and Lahad Dato to the southwest of Sabah. The seismicity record

19

of Ranau starts in 1991, according to the international seismicity bulletins (e.g. USGS

20

and International Seismological Center ISC), and this short record is not sufficient for

21

seismic source characterization. Fortunately, active Quaternary fault systems are

22

delineated in the area. Henceforth, seismicity of the area is thus determined as line

23

sources referring to these faults. Two main fault systems are believed to be the source

24

of such activities; namely, the Mensaban fault zone (MFZ) and the Crocker fault zone

25

(CFZ) in addition to some other faults in their vicinity. Seismic hazard assessments

26

became a very important and needed study for the extensive developing projects in this

27

state especially with the presence of earthquake activities. Probabilistic seismic hazard

28

assessments are adopted for the present work since it can provide the probability of

29

various ground motion levels during expected future large earthquakes. The output

30

results are presented in terms of spectral acceleration curves and uniform hazard

ce

pte

dM

12

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

curves for periods of 500, 1000 and 2500 years. Since this is the first time that a

32

complete hazard study has been done for the area, the output will be a base and

33

standard for any future strategic plans in the area.

35

us cri

34

pt

31

Keywords: Sabah seismicity, probabilistic seismic hazards, Ranau.

36

Introduction

38

For regions with documented seismic activities, seismic hazard assessments are an important

39

task for the preparedness and mitigation of their effects. For a developing country like Malaysia,

40

such importance is even more vital to ensure the sustainability of the developing projects both

41

for industries and urbanization purposes. Sabah is currently witnessing a considerable number of

42

development projects due to Malaysia’s government national agenda to transform the country

43

into a high-income nation by the year 2020 (SEDIA, 2015). The Economic Transformation

44

Programme by the government separates Sabah into 6 strategic development areas with each area

45

focusing on a particular aspect of the economy, such as oil & gas, tourism (northwest Sabah) and

46

agriculture (SEDIA, 2015). Furthermore, due to the expected increase in the demand for energy,

47

a coal power plant was proposed to be built in Lahad Datu, Sabah by a local company (Koh and

48

Lim, 2010). This and many more development projects motivated the present study as to

49

determine those expected ground motion levels during future large earthquakes.

50

The probabilistic seismic hazards assessments approach (PSHA) is adopted in the

51

present study, which gives the probability of ground motion levels for various magnitude

52

earthquakes.

53

researchers in various localities all over the globe (e.g. Sokolov et al. (2004); Faenza et

54

al. (2007); Gupta (2007); Primer (2008); El-Hussain et al. (2012); Liu et al. (2013);

55

Ordaz et al. (2014); Zahran et al. (2015); Ezzelarab et al. (2016); Liu et al. (2016)). The

56

seismic hazard is calculated based on the procedure that was introduced by Cornell

pte

dM

an

37

Henceforth, PSHA is receiving an increasing interest by many

ce

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 29

57

(1968) and modified by McGuire (1978) and Bender and Perkins (1987).

The

58

probabilistic seismic hazard assessment is calculated using variable data (geological,

Page 3 of 29

seismological and structural) to construct a model of earthquake ground motion at the

60

site of interest. This technique follows four steps; firstly, a complete analysis of the

61

historical and recent seismicity in the area was done, followed by determining and

62

identifying the seismotectonic source model for the area. Secondly, recurrence

63

parameters for seismicity, including the expected maximum earthquake, were estimated

64

for each determined seismic source. Finally, a predictive ground motion is used to

65

describe the expected ground motion as a function of different magnitudes and

66

distances.

67

The region of interest, i.e. Sabah state, is the most seismically active region in Malaysia.

68

Several moderate earthquakes with magnitudes around 6 on the Richter scale shocked

69

the state in the past. These earthquakes caused moderate damage to properties. For

70

Ranau area, damage due to earthquakes is more dangerous. The only life casualties

71

reported are caused by the June 5th, 2015 earthquake with a moment magnitude

72

Mw=6.0. Another event in Ranau area, which took place in 1991, also damaged some

73

buildings. The events near Lahad Dato, on the other hand, are reported to produce

74

moderate damage with no reported life casualties. It is evident that the events near

75

Ranau are more dangerous compared to those near Lahad Dato. The reason may be

76

due to its proximity to Kota Kinabalo which comprises the dense population. Moreover,

77

the only life casualties reported affected the climbers at mount Kinabalo. According to

78

Tjia (2007) and the Meteorological Department of Malaysia (GMDM, 2006), the

79

earthquakes near Ranau are triggered due to the movement of Quaternary active faults

80

in the area. Kundasang-Ranau is located near the intersection of regional fault zones.

81

These faults are namely Crocker fault zone (CFZ) and Mensaban fault zone (MFZ). The

82

widespread mass movements of the Kundasang area can baaine partially attributed to

83

such active movements and intersection. Despite such remarkable fault movements, the

84

seismicity record at the area is short. Previous work conducted in the area (e.g. Leyu et

85

al. (1985); Che Abas (2001); Harith et al. (2014)), used the deterministic seismic

86

hazards approach while considering only the regional sources. Deterministic seismic

87

hazard approach depends on the so-called characteristic earthquake which can be

88

defined as the maximum expected earthquake affecting the site. Hence, the effects of

89

other earthquakes are not considered, which may mask important effects. Furthermore,

ce

pte

dM

an

us cri

pt

59

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

their work didn’t consider the local seismicity near Ranau which may have dangerous

91

effects on neighboring sites like Kota Kinabalo.

92

Seismotectonic Setting

93

In general, the region of northern Borneo has a moderate rate of earthquakes

94

influenced by the local tectonics with the biggest reported earthquake occurred in 1923

95

with a magnitude of 6.9 ((Rangin et al. 1999); (Engdahl and Villasenor 2002); (Simons

96

et al. 2007); (Mark Petersen et al. 2007)). The local tectonics of Sabah is controlled

97

mainly by the rifting episodes of the South China Sea generating NE-SW structures and

98

the opening of the Sulu Sea basin generating NW-SE structures (e.g. (Tan and Lamy

99

1990); (Tongkul 1993)). Moreover, Sabah has experienced episodes of compression

100

since Cretaceous, which can be observed from the existence of folds and thrust faults,

101

however, an episode of extension occurred during the late early Miocene which

102

culminates in NE-SW trending normal faults (Tongkul 1993).

103

Such a complex deformation history controlled mainly by the tectonic evolution of the

104

surrounding deep basins which are the south China Sea basin to the north, the Sulu

105

Sea basin to the northeast and the Celebes Sea basin to the southeast (Hall et al,

106

2008). Consequently, three major episodes of compression along the northwestern

107

boundary of Sabah occurred during the Late Eocene, Early Miocene and Late Miocene

108

periods (Balaguru and Hall, 2008). This is mainly due to the subduction of the south

109

China Sea beneath Borneo resulting in thrust faults and folds. The presence of the

110

Crocker-Rajang mountain belt is a clear indication of the compressional environment

111

along the northwest of Sabah ((Franke et al. 2008)). Moreover, the last tectonic event

112

that is probably still ongoing is the formations of sinistral transpressional faulting in the

113

southeast of Sabah, which is believed to be the result of the propagation of deformation

114

from Sulawesi (Balaguru & Hall, 2008). This ongoing transpressional deformation is

115

observed to affect the local seismicity of Sabah.

116

Furthermore, the seismotectonic studies conducted by a number of researchers (e.g.

ce

pte

dM

an

us cri

pt

90

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 29

117

(Leyu et al. 1985); (Che Abas 2001); (Wah 2011)) shows that the state of Sabah exhibit

118

clear rate of crustal deformation. Sabah is believed to owe its moderate seismicity to the

Page 5 of 29

active Mensaban and Lobou-Lobou fault zones which have brought about earthquakes

120

that caused light damages to infrastructure

121

According to (Tjia 2007), the Crocker zone is extending to about 170 Km in length

122

comprising the Quaternary structural valleys of Tenom, Keningau, and Tambunan.

123

These valleys are trending N-S to N 20o E. The sense of motion for these family of

124

faults is Normal dip-slip with Sinistral Strike slip component. Surface deformations of

125

tarred roads clearly claim that these faults are still active. An important question arises

126

here as for whether these activities are either seismic or Aseismic (creep). Seismic

127

movements are transient and sudden which characterizes areas of brittle rocks. Creep,

128

on the other hand, ductile rocks regions. In ductile zones, the accumulated strains are

129

dissipated in the form of slow movements. The recorded earthquake activities there

130

(Figure 1) claims these faults may have both seismic movements and creep.

131

In addition to Ranau, two other local activity can be observed from the local seismicity

132

map of Sabah (figure 2).

133

Semporna in the southeast. The seismicity at Semporna may be related to the

134

Kalimantan subduction zone and thus will be treated as such.

135

Regional seismicity is encountered also at the subduction zones either in the south

136

Philippines Sea or in the east Indonesia at Kalimantan. In addition to these, other

137

seismic activities are also present but they are excluded because of the large distances

138

separating them from the investigated area.

139

The tectonics of south Philippines Sea has been studied by several researchers (e.g.

140

(Acharya 1979); (Hori 2006); (Galgana et al. 2007); (Ramos and Tsutsumi 2010); (Lin

141

and Lo 2013)). The distance between the south Philippine Sea seismicity and the state

142

of Sabah is around 600 Km indicating that only large earthquakes of magnitudes greater

143

than 7 can produce some damage specifically on the eastern parts (Figure 3). The

144

tectonic activities of the south Philippine Sea are concentrated along the Manilla trench,

145

the Philippine trench, and the Philippine fault system. Acharya (1979) indicated that the

an

us cri

pt

119

pte

dM

These two zones are theSundakan in the northeast and

ce

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

146

west Philippine basin has moderate to high seismicity with magnitudes up to 7.9. From

147

historical data, we cannot identify damage in Sabah related to the activity in the

148

southwestern Philippine region.

us cri an

149

dM

Figure 1: Faults and earthquake epicenters around Ranau northwest Sabah.

ce

pte

150

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 29

pt

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

151 152

Figure 2: Local seismicity of Sabah.

Page 7 of 29

Seismicity from south Indonesia is higher due to the subduction of the Indian and

154

Australian plates beneath the Sunda plate (e.g. (Mark Petersen et al. 2007); (Katili

155

1989)). High seismicity is observed along the Sunda subduction, where shortening and

156

strain are accommodated by thrust faults along the collision boundary. However, such

157

activities are so far from the study region that its effect is neglected. Henceforth, south

158

Indonesia sources will not be considered the present seismic hazards assessments.

159

Probabilistic Seismic Hazard

160

The technique was first introduced to the community by (Cornell 1968), who defined the

161

technique as the plot of the relationship between ground motion levels (e.g. intensity,

162

peak ground velocity, and peak ground acceleration) with their average return periods.

163

The concept has been developed over the years by many researchers (e.g. McGuire

164

and Shedlock 1981; Kijko and A. 1992; Kijko and Graham 1999; McGuire R 2001;

165

Mäntyniemi, Tsapanos, and Kijko 2004; McGuire 2008; Gupta 2013; Wong 2013;

166

Pappin et al. 2015).

167

According to the abovementioned definition, earthquake source zones in the

168

neighborhood of the inverstigated site must be defined. The definition involves the

169

determination of the spatial extention, the frequency of earthquakes with a certain size

170

(e.g magnitudes or seismic moments) and their expected return periods. The spatial

171

extention is delineated via numerous information such as the tectonic regeme, focal

172

mechanism, distribution of earthquake foci, Ietc. Within each earthquake source, the

173

seismic activity is assumed to be homogeneous. The frequency of earthquakes and

174

their return period, on the other hand, is derived from the magnitude frequency relation

175

of (Gutenberg and Richter 1944):

us cri

an

dM

pte

ce

176

pt

153

log =  −

(1)

177

where n is the number of earthquakes, M is their magnitude range, ‘a’ is a measure of

178

seismicity (rarthquake rate) and ‘b’ is constant related to the tectonic setting. The later

179

has a value ranging generally between 0.8 and 1.2. However, since the original relation

180

doesn’t exhibit a maximum earthquake magnitude, a doubly truncated equation was

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

proposed by (Cornell and Vanmarcke 1969). From the frequency-magnitude relation,

182

the return period of each magnitude level can be determined. Such parameter is

183

important for the determination of its probability.

184

After

185

activity

186

must

be

calculated.

187

equations

(GMPE).

188

relating the ground motion at the site to the magnitude and distance of

189

the

190

mathematically in the form (simplified after (Akkar and Bommer 2007):

determination

rates,

the

levels

the

of

their

This

is

ground

done

These

source.

earthquzke

source

zones

and

their

motion

using

the

equations

Such

at

the

ground

are

site

of

motion

prediction

empirically

relation

can

be

interest

generated

represented

an

earthquake

of

us cri

the

pt

181

log  =

+  +   +   +   log +   +   2 Where A is ground motion parameter, M is the size of the earthquake (i.e magnitude or

192

moment), R is the distance between the earthquake source and the site, S is the soil

193

factor and F is the fault type factor. The parameters bi, where i=1 to 7, is the model

194

parameters to be determined by data fitting. GMPE can also be estimated at certain

195

spectral frequencies of the ground motion to determine the spectra of ground motion at

196

the site. The final stage after the calculations of the ground motion levels from different

197

sources and their respective return periods is to build the probability model that

198

determines the overall probability of exceedance of each level.

199

Earthquake sources:

200

Local

201

Malaysian

202

main

203

Semporna.

204

Kota Kinabalo, the capital city of the state. The Ranau activity is most

205

likely

pte

dM

191

seismicity

206 207

of

territories.

ce

locations.

an

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 29

The

Sabah

Moderate

These

first

intraplate

is

is

distinctive

from

earthquakes

locations

are

the

dangerous

more

earthquake

activity

are

near

that

other

experienced Ranau,

as

it’s

resulted

along local quaternary faults (Tjia (2007) and Wah (2011)).

parts

of at

the three

Sandakan, located

from

close

and to

movements

us cri an dM

208

Figure 3: Regional seismicity with possible impacts on Sabah.

210

Seismological data of Ranau shows that the earliest earthquake in this region started by

211

the early 1990’s. From the temporal point of view, three clusters can be identified. The

212

first took place on May 26th, 1991 represented by four events with magnitudes ranging

213

from 4.6 to 5.4. Another event took place in 1995 with magnitude mb of 4.1. Afterwards,

214

the area was quite until 2015 when a large earthquake with magnitude Mw=6.0 shocked

215

the region. The lack of earthquake data before 1991 is odd and may be interpreted as

216

being related to the distribution of seismological stations in the region. Thus, local

217

seismicity prior 1991 maybe either not detected or misallocated. Accordingly, the

218

earthquake data is clearly not sufficient for determining the basic parameters for such

219

seismic source characterization. Alternatively, Wells and Coppersmith (1994) proposed

220

an alternative source description based on the data driven from active faults.

ce

pte

209

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 9 of 29

221

Fortunately, the quaternary active faults in the area were investigated by Tjia (2007).

222

From Figure 1 above, the distribution of recorded seismicity and quaternary faults in the

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

region are presented. The seismicity of this area will be related to the lengths of the

224

active faults there.

225

A similar situation prevails in Sandakan region. Seismicity there is scarce and too few

226

for the fitting of the frequency magnitude relation. Henceforth, it is not possible to

227

determine the recurrence of seismic activity there. Unlike Ranau, no information about

228

active faulting is available. This makes the delineation of required seismic source

229

parameters an uneasy task. Hence, the activity there is omitted for the present study.

230

Activity at Semporna region, on the other hand, is considered as belonging to

231

Kalimantan subduction zone of east Indonesia. Semporna is closely located near

232

Kalimatan with seismicity covering the entire region without identifiable borders.

233

Henceforth, it will be dealt with as part of regional seismicity rather than local one.

234

From the seismotectonic discussion in the previous section, it is identified that we have

235

two potential regional sources that may cause the considerable shaking level to Sabah.

236

These sources are the South Philippines and Kalimantan subduction zone. The

237

identification of seismic sources adopted in the present work was taken from the work of

238

(Irsyam et al. 2010). The needed parameters for the seismic hazards calculations is

239

shown in figure 4. Sources considered are those at distances less than 600 km that are

240

expected to have considerable effects on Sabah.

241

Attenuation model

242

There are several GMPE that can be can be categorized in terms of their tectonic

243

environment (i.e. subduction zone and shallow crustal earthquakes). Henceforth. the

244

selection of an appropriate attenuation relationship is one of the critical factors in PSHA.

245

As the tectonic environments near Sabah comprise both abovementioned types,

246

different forms of GMPE will be used. For subduction zone earthquakes, several

247

attenuation relationships are derived such as (Crouse 1991), (Youngs et al. 1997),

248

(Atkinson and Boore 1997), (Petersen et al. 2004), whereas for shallow crustal

249

earthquakes, attenuation relationships developed by (Campbell 1997); (Campbell

250 251

ce

pte

dM

an

us cri

pt

223

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 29

2003), (Toro, Abrahamson, and Schneider 1997) are frequently used.

us cri an

252

Figure 4: Earthquake source definitions in the region around the area under

254

investigation (after Irsyam et. al., 2010)

255

Previous works assumed that the damage of the earthquake in Sabah Borneo was

256

triggered solely by the large earthquakes from Sulawesi, Indonesia, and Southern

257

Philippines which are all subduction zones. Hence, previous studies (e.g. Harith et al.

258

(2014) and Adnan and Marto (2008)) have used the GMPE model for subduction zone.

259

The later used Campbell (2003) for calculation of distant earthquake while (Sadigh et al.

260

1997) and (Boore, Joyner, and Fumal 1997) for short distance earthquake. Sabah is

261

located on the stable continent of Sunda Plate and there’s no GMPE which have been

262

derived specifically for this area. As an alternative approach, attenuation models had

263

been developed from analogous regions which were considered to possess similar

264

seismo-tectonics and geological conditions. Therefore, in this study, a GMPE which

265

derived for the stable continent for another region was applied. For subduction zone, on

266

the other hand, such as the Philippines and Sulawesi megathrust Young et al. (1997) is

267

applied because it possesses a similar unit and covers the same range of spectral

268

response period with GMPE for the stable continent as mentioned above.

ce

pte

dM

253

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 11 of 29

269

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

Result Interpretation

271

Probabilistic seismic hazard assessments were conducted for six cities belonging to the

272

state of Sabah. For each site, the results are represented in terms of a uniform hazard

273

curves for the return periods of 500, 1000 and 2500 years. Moreover, the exceedance

274

probability curves of expected intensity (acceleration in cm/s2) are also presented for

275

the periods of 0.01, 1.0 and 2.0 seconds (figures 5 thru 10). Finally, the spectral

276

intensity (acceleration) is also mapped for the next 500, 1000 and 2500 years at the

277

periods of 0.5s, 1s, 2s, and 3s (figures 11, 12 and 13.)

278

The results showed that the expected ground motion levels at both Ranau and Lahad

279

Dato are the highest, with the later showing the highest expected ground motion levels

280

in the return periods considered in the present work. The reason for this is the proximity

281

of both sites to active earthquake sources. Compared to previous studies (e.g (Adnan

282

and Marto 2008)) shown in figure 14 here, the results for Ranau deduced from the

283

present study show higher level of expected ground shaking. This is a direct impact of

284

considering local seismic sources for the present work.

285

The results for Ranau (figure 5), shows generally lower annual probability of

286

occurrences (panels a, b and c) compared to their respective in both Tawau and Lahad

287

Dato (figures 7 and 9). This is a consequence of lower rate of local seismicity at Ranau

288

compared to that of the Kalimantan subduction zone. The uniform hazard curve (panel d

289

in the figures), on the other hand, illustrates that the spectral acceleration near Ranau is

290

the second highest after Lahad Dato.

291

Other cities considered in the present study, show the relatively lower level of expected

292

ground motions. However, the magnitude of ground motion shows a general decrease

293

from southeast to the northwest until Ranau and Kota Kinabalo where local seismicity

294

produces anomalous ground motion trend.

us cri

an

dM

pte

ce

295

pt

270

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 29

Ranau 1.0 Sec

1.00E+01 1.00E-01 1.00E-03 1.00E-05 1.00E-07 1.00E+00

2.50E+03

1.00E-08 1.00E-11

c. 296

1.00E+02

302

hazard curve is shown at d)

pte

301

1.00E+04

d.

Fig. 5. Spectral acceleration for Ranau at periods of a) 0.01s, b) 1s and c) 2s. Uniform

ce

300

1.00E+02

an

1.00E-05

dM

Annual Prob. of Exceedance

1.00E-02

Spectral Acceleration (cm/sec2)

299

1.00E-10

Spectral Acceleration (cm/s2)

b.

1.00E-14 1.00E+00

298

1.00E-07

1.00E+00

(cm/s2)

Ranau 2 S

297

1.00E-04

1.00E-13

5.00E+01

Spectral Acceleration

a.

1.00E-01

us cri

Annual propability of exceedance

Annual propability of exceedance

Ranau (0.01 sec)

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 13 of 29

1.00E-04 1.00E-06 1.00E+02

Spectral

1.00E+04

Acceleration (cm/s2)

b.

a.

1.00E+00 1.00E-03 1.00E-06 1.00E-09 1.00E-12 1.00E-15 1.00E+00 1.00E+01 1.00E+02 1.00E+03

c. 303 304 305

Uniform hazard curve is shown at d)

pte

309

ce

308

1.00E-07 1.00E-10

1.00E-13 1.00E+00

1.00E+02

1.00E+04

Spectral Acceleration (cm/s2)

Fig. 6. Spectral acceleration for Kota Kinabalo at periods of a) 0.01s, b) 1s and c) 2s.

306 307

1.00E-04

d.

dM

Spectral acceleration (cm/s2)

1.00E-01

an

Annual Prob. of Exceedance

Kota Kinabalo 2 S

pt

1.00E-02

us cri

1.00E+00

1.00E-08 1.00E+00

Page 14 of 29

Kota Kinabalo 1 S Annual Prob. of Exceedance

Kota Kinabalo 0.01 S

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Annual prob. of exceedane

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

a

Annual Prob. of Exceedance

1.00E+00 1.00E-02 1.00E-04 1.00E-06 1.00E-08 1.00E+02

1.00E+04

b

Spectral

1.00E+02

310

1.00E-09 1.00E-11

1.00E-13 1.00E+00

1.00E+02

1.00E+04

d

Fig. 7. Spectral acceleration for Tawau at periods of a) 0.01s, b) 1s and c) 2s. Uniform

pte

hazard curve is shown at d)

ce

314

1.00E-07

1.00E+04

Acceleration (cm/s2)

c

313

1.00E-05

an

1.00E+01 1.00E-01 1.00E-03 1.00E-05 1.00E-07 1.00E-09 1.00E-11 1.00E-13 1.00E-15 1.00E+00

dM

Annual Prob. of Exceedance

Tawau 2 S

312

1.00E-03

Spectral Acceleration (cm/s2)

1.00E-10 1.00E+00

Spectral Acceleration (cm/s2)

311

1.00E-01

us cri

Annual Prob. of Exceedance

Tawau 1 S

Tawau 0.01 S

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 15 of 29

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

Sandakan 0.01 S

pt

Sandakan 1 S

1.00E-02

1.00E-05

1.00E-08

1.00E-11 1.00E+00

1.00E+02

b

1.00E-03 1.00E-06 1.00E-09 1.00E-12 1.00E-15 1.00E-01

1.00E+01

1.00E-12

1.00E+02

1.00E+04

1.00E+03

Spectral Acceleration (cm/s2)

c 315

d

pte

Fig. 8. Spectral acceleration for Sandakan at periods of a) 0.01s, b) 1s and c) 2s. Uniform hazard curve is shown at d)

ce

319

1.00E-09

an

1.00E+00

dM

Annual Prob. of Exceedance

Sandakan 2 S

318

1.00E-06

Spectral Acceleration (cm/s2)

a

317

1.00E-03

1.00E-15 1.00E+00

1.00E+04

Spectral Acceleration (cm/s2)

316

1.00E+00

us cri

Annual Prob. of Exceedance

Annual Prob. of Exceedance

1.00E+01

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 29

Lahad Dato 1 S

Lahad Dato 0.01 S 1.00E+00 1.00E-01 1.00E-02 1.00E-03 1.00E-04 1.00E-05 1.00E-06 1.00E+00 1.00E+01 1.00E+02 1.00E+03 1.00E+04

b

1.00E-01

1.00E-09

1.00E+02

1.00E+04

1.00E-04 1.00E-07 1.00E-10 1.00E-13 1.00E-16

1.00E+00

dM

Annual Prob. of Exceedance

1.00E-07

an

Lahad Dato 2 S

1.00E+02

1.00E+04

Spectral Acceleration (cm/s2)

pte

c

d

Fig. 9. Spectral acceleration for Lahad Dato at periods of a) 0.01s, b) 1s and c) 2s. Uniform hazard curve is shown at d)

ce

322

1.00E-05

Spectral Acceleration (cm/s2)

a

321

1.00E-03

1.00E-11 1.00E+00

Spectral Acceleration (cm/s2)

320

1.00E-01

us cri

Annual Prob.of Exceedance

Annual prob. of Exceedance

1.00E+01

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 17 of 29

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

1.00E+00 1.00E-01 1.00E-03 1.00E-05 1.00E-07 1.00E-09 1.00E-11 1.00E+00

1.00E+02

1.00E+04

1.00E-02 1.00E-04 1.00E-06

us cri

Annual Prob. of Exceedance

Annual Prob. of Exceedance

pt

Keningau 1 S

Keningau 0.01 S

1.00E-08 1.00E-10 1.00E-12 1.00E-14 1.00E-16

1.00E+00

Spectral Acceleration (cm/s2)

a

1.00E+02

1.00E+04

Spectral Acceleration (cm/s2)

b

an

1.00E-01 1.00E-03 1.00E-05 1.00E-07 1.00E-09 1.00E-11 1.00E-13

dM

Annual Prob. of Exceedance

Keningau 2 S

1.00E-15 1.00E+00 1.00E+01 1.00E+02 1.00E+03 Spectral Acceleration (cm/s2)

c

325 326

d

Fig. 10. Spectral acceleration for Keningau at periods of a) 0.01s, b) 1s and c) 2s. Uniform hazard curve is shown at d)

ce

324

pte

323

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 29

us cri an dM 329 330 331

pte

328

Fig. 11 Spectral ground motion at Sabah for 2500 years return Period at the periods a) 0.5 s, b) 1 s, c) 2 s and d) 3 s.

ce

327

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 19 of 29

332

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

335 336 337 338

Fig. 12. Spectral ground motion at Sabah for 2500 years return Period at the periods a) 0.5 s, b) 1 s, c) 2 s and d) 3 s.

ce

334

pte

dM

an

us cri

pt

333

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 29

339

us cri an 342 343

pte

dM 341

ce

340

Fig 13. Spectral ground motion at Sabah for 2500 years return Period at the periods a) 0.5s, b) 1s, c) 2s. and d) 3s.

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 21 of 29

345 346

us cri an

344

Figure 14: Peak ground acceleration (PGA) maps for 500 years and 2500 years (after (Adnan and Marto 2008))

dM

347

Conclusion

349

Probabilistic seismic hazard assessments are applied for six cities in Sabah state, east Malaysia.

350

The study is fueled by two main reasons. First, Sabah shows extraordinary earthquake activities

351

compared to other regions of Malaysia. Second, the occurrence of a moderate earthquake with

352

magnitude Mw=6.0 in June 2015 which produce some damages to the building and caused a

353

death of some tourists near the Kinabalo mountain. Previous seismic hazard assessments were

354

mostly deterministic and ignored the local earthquake sources near Ranau. The present study,

355

on the other hand, is designed to overcome these limitations of the previous ones.

356

The earthquake records of local seismicity near Ranau is extremely short with the earliest activities

357

dated in 1991. Since then, few events occurred until the largest one took place in June 2015. Such

358

seismic data is insufficient to characterize the source parameters near Ranau. Fortunately, geologic

359

studies revealed that earthquake activities are related to Quaternary active faults. The lengths of the

360

delineated faults are adopted using (Donald L. Wells and Coppersmith 1994) model to characterize the

ce

pte

348

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 29

pt

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

361 362

local seismic activity there.

Page 23 of 29

Furthermore, adopting the PSHA technique was not carried out in Sabah before this work. The

364

advantage is to define the effects of different ground motion levels rather than the maximum

365

earthquake alone. This will change the measures to design earthquake resistant structure considerably.

366

Moreover, the results obtained showed that the ground motion levels estimated for return periods of

367

500, 1000 and 2500 years, are largest near Lahad Dato and decrease in the northwest direction until

368

Ranau. Around Ranau, the ground motion levels show an anomalous increase due to the local activities

369

there. This ground motion anomaly was omitted in the previous studies. Henceforth, the ground motion

370

levels obtained here is more realistic than previous ones. The new assessments maybe useful for both

371

preparedness and mitigation of earthquakes in the state of Sabah.

372

References

373

Acharya, H. K. 1979. “Seismicity of the Southern Philippine Sea.” Marine Geology 29(1–

376 377

us cri

an

375

4):25–32.

Adnan, A. and Aminaton, M.. 2008. “Development of Seismic Hazard Maps of East Malaysia.” Advances in Earthquake Engineering Applications 1–17.

dM

374

pt

363

Akkar, S. and Bommer, J. 2007. “Empirical Prediction Equations for Peak Ground Velocity Derived from Strong-Motion Records from Europe and the Middle East.”

379

Bulletin of the Seismological Society of America 97(2):511–30. Retrieved August

380

11, 2017

381

(https://www.researchgate.net/profile/Sinan_Akkar/publication/228372953_Empiric

382

al_Prediction_Equations_for_Peak_Ground_Velocity_Derived_from_Strong-

383

Motion_Records_from_Europe_and_the_Middle_East/links/0046351ea73d483ad50

384

00000.pdf).

386 387

Atkinson, Gail M. and David M. Boore. 1997. “Some Comparisons Between Recent Ground-Motion Relations.” Seismological Research Letters 68(1):24–40. Retrieved

ce

385

pte

378

(http://srl.geoscienceworld.org/content/68/1/24).

388

Balaguru, Allagu and Robert Hall. 2008. “Tectonic Evolution and Sedimentation of

389

Sabah, North Borneo, Malaysia.” American Association of Petroleum "

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

390

30084(1990):1–7. Retrieved

391

(http://www.searchanddiscovery.com/documents/2009/30084balaguru/images/bala

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

393

guru).

pt

392

Boore, D. M., W. B. Joyner, and T. E. Fumal. 1997. “Equations for Estimating Horizontal Response Spectra and Peak Acceleration from Western North American

395

Earthquakes: A Summary of Recent Work.” Seismological Research Letters

396

68(1):128–53.

397

us cri

394

Campbell, Kenneth W. 1997. “Empirical Near-Source Attenuation Relationships for Horizontal and Vertical Components of Peak Ground Acceleration, Peak Ground

399

Velocity, and Pseudo-Absolute Acceleration Response Spectra.” Seismological

400

Research Letters 68(1):154–79. Retrieved

401

(http://srl.geoscienceworld.org/content/68/1/154).

402

an

398

Campbell, Kenneth W. 2003. “Prediction of Strong Ground Motion Using the Hybrid Empirical Method and Its Use in the Development of Ground-Motion (Attenuation)

404

Relations in Eastern North America.” Bulletin of the Seismological Society of

405

America 93(3):1012–33. Retrieved April 25, 2017

406

(http://www.ce.memphis.edu/7137/PDFs/attenuations/USGS2008/Campbell2003.p

407

df).

dM

403

Che Abas, M. R. 2001. “Earthquake Monitoring in Malaysia.” Seismic Risk Seminar.

409

Cornell, A. 1968. “Engineering Seismic Risk Analysis.” BSSA 58(5):1583–1606.

410

Cornell, C. A. and E. H. Vanmarcke. 1969. “The Major Influences on Seismic Risk.” Pp.

411

69–83 in Proc. of the 4th World Conference on Earthquake Engineering. Retrieved

412

August 11, 2017 (http://www.iitk.ac.in/nicee/wcee/article/4_vol1_A1-69.pdf).

413

Crouse, C. B. 1991. “Ground‐Motion Attenuation Equations for Earthquakes on the

415

Cascadia Subduction Zone.” Earthquake Spectra 7(2):201–36. Retrieved

ce

414

pte

408

(http://dx.doi.org/10.1193/1.1585626).

416 417 418

EL-HUSSAIN, I., DEIF, A., AL-JABRI, K., TOKSOZ, N., EL-HADY, S., AL-HASHMI, S., AL-TOUBI, K., AL-SHIJBI, Y., AL-SAIFI, M. & KULELI, S. 2012. Probabilistic seismic hazard maps for the sultanate of Oman. Natural Hazards, 64, 173-210.

419 420

Engdahl, E. R. and A. Villasenor. 2002. “Global Seismicity: 1900-1999.” Retrieved

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 29

Page 25 of 29

(https://earthquake.usgs.gov/data/centennial/centennial.pdf).

pt

421

EZZELARAB, M., SHOKRY, M. M. F., MOHAMED, A. M. E., HELAL, A. M. A., MOHAMED, A. A. & EL-HADIDY,

423

M. S. 2016. Evaluation of seismic hazard at the northwestern part of Egypt. Journal of African Earth

424

Sciences, 113, 114-125.

us cri

422

425 426 427 428 429

FAENZA, L., HAINZL, S., SCHERBAUM, F. & BEAUVAL, C. 2007. Statistical analysis of time-dependent earthquake occurrence and its impact on hazard in the low seismicity region Lower Rhine Embayment. Geophysical Journal International, 171, 797-806.

430

Franke, Dieter, Udo Barckhausen, Ingo Heyde, Mark Tingay, and Nordin Ramli. 2008. “Seismic Images of a Collision Zone Offshore NW Sabah / Borneo.” Marine and

432

Petroleum Geology 25:606–24. Retrieved April 24, 2017

433

(https://www.researchgate.net/profile/Mark_Tingay/publication/236011871_Seismic

434

_images_of_a_collision_zone_offshore_NW_SabahBorneo/links/0046353a187d49f

435

f33000000.pdf).

dM

436

an

431

Galgana, Gerald, Michael Hamburger, Robert McCaffrey, Ernesto Corpuz, and Qizhi

437

Chen. 2007. “Analysis of Crustal Deformation in Luzon, Philippines Using Geodetic

438

Observations and Earthquake Focal Mechanisms.” Tectonophysics 432(1–4):63–

439

87.

441 442

Gupta, I. D. 2013. “Source-to-Site Distance Distributions for Area Type of Seismic Sources Used in PSHA Applications.” Natural Hazards 66(2):485–99.

pte

440

Gutenberg, B. and C. F. Richter. 1944. “Frequency of Earthquakes in California.” Bulletin of the Seismological Society of America 34:185–188. Retrieved August 10,

444

2017 (http://authors.library.caltech.edu/47734/1/185.full.pdf).

445 446 447

ce

443

Hall, Robert, Marco W. A. van Hattum, and Wim Spakman. 2008. “Impact of India-Asia Collision on SE Asia: The Record in Borneo.” Tectonophysics 451(1–4):366–89.

HARITH, N. S. H., ADNAN, A. & SHOUSHTARI, A. V. 2014. Estimation of peak ground acceleration map of

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

448

Sabah Based on deterministic seismic hazards analysis. National Seminar on Civil Engineering

449

Research (SEPKA). Universiti Teknologi Malaysia (UTM), Johor Bahru, Skudai, Malaysia

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

Hori, Sadaki. 2006. “Seismic Activity Associated with the Subducting Motion of the

pt

450 451

Philippine Sea Plate beneath the Kanto District, Japan.” Tectonophysics 417(1–

452

2):85–100. Retrieved (http://dx.doi.org/10.1016/j.tecto.2005.08.027).

Irsyam, M. et al. 2010. “Development of Seismic Hazard and Risk Maps for New

us cri

453 454

Seismic Building and Infrastructure Codes in Indonesia.” Pp. 1–8 in The 6th Civil

455

Engineering Conference in Asia Region: Embracing the Future through

456

Sustainability.

458 459

Katili, John A. 1989. “Review of Past and Present Geotectonic Concepts of Eastern Indonesia.” Netherlands Journal of Sea Research 24(2–3):103–29.

Kijko, Andrzej and Sellovel M. A. 1992. “Estimation of Earthquake Hazard Parameters

an

457

460

from Incomplete Data Files. Part II. Incorporation of Magnitude Heterogeneity.”

461

Bulletin of the seismological soceity of America 82(1):120–34. Kijko, A. and Graham, G. 1999. “‘“Parametric-Historic”’ Procedure for Probabilistic

dM

462 463

Seismic Hazard Analysis Part II: Assessment of Seismic Hazard at Specified Site.”

464

Pure and applied Geophysics 154:1–22.

Koh, S. L. and Lim, Y. S. 2010. Meeting energy demand in a developping economy

466

without damaging the environment- A case study in Sabah, Malaysia, from

467

technical, environmental and economic perspective. Energy Policy, 38 (8), 4719-

468

4728.

469

pte

465

Leyu, Chong-Hua et al. 1985. Series on Seismology, Volume III -MALAYSIA. edited by

470

E. P. Arnold. SOUTHEAST ASIA ASSOCIATION OF SEISMOLOGYAND

471

EARTHQUAKF ENGINEERINC. LIU, J., WANG, Z., XIE, F. & LV, Y. 2013. Seismic hazard assessment for greater North China from historical intensity observations. Engineering Geology, 164, 117-130.

474 475

LIU, J., XIE, F. & LV, Y. 2016. Seismic hazard assessments for the Ordos Block and its periphery in China.

476

ce

472 473

Soil Dynamics and Earthquake Engineering, 84, 70-82.

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 29

477 478

Lin, Jing Yi and Chung Liang Lo. 2013. “Earthquake-Induced Crustal Gravitational Potential Energy Change in the Philippine Area.” Journal of Asian Earth Sciences

479

66(November 2011):215–23. Retrieved

480

(http://dx.doi.org/10.1016/j.jseaes.2013.01.009).

481

Mäntyniemi, Päivi, Theodoros M. Tsapanos, and Andrzej Kijko. 2004. “An Estimate of Probabilistic Seismic Hazard for Five Cities in Greece by Using the Parametric-

483

Historic Procedure.” Engineering Geology 72(3–4):217–31.

484

us cri

482

Mark Petersen et al. 2007. “Documentation for the Southeast Asia Seismic Hazard

485

Maps.” 67. Retrieved (http://pubs.usgs.gov/circ/1383f/Circ1383-F.pdf).

486

McGuire, R. K. and K. M. Shedlock. 1981. “STATISTICAL UNCERTAINTIES IN

SEISMIC HAZARD EVALUATIONS IN THE UNITED STATES.” BSSA 71(4):1287–

488

1308.

490 491 492 493

McGuire, R. K. 2008. “Probabilistic Seismic Hazard Analysis: Early History.” Earthquake Engineering & Structural Dynamics 37(3):329–38.

McGuire R. 2001. “Deterministic vs. Probabilistic Earthquake Hazards and Risks.” Soil

dM

489

an

487

Dynamics and Earthquake Engineering 21(5):377–84. ORDAZ, M. G., CARDONA, O.-D., SALGADO-GÁLVEZ, M. A., BERNAL-GRANADOS, G. A., SINGH, S. K. &

494

ZULOAGA-ROMERO, D. 2014. Probabilistic seismic hazard assessment at global level. International

495

Journal of Disaster Risk Reduction, 10, 419-427.

Pappin, J. W. et al. 2015. “A Rigorous Probabilistic Seismic Hazard Model for Southeast

pte

496 497

China: A Case Study of Hong Kong.” Bulletin of Earthquake Engineering

498

13(12):3597–3623.

499

Petersen, Mark D. et al. 2004. “Probabilistic Seismic Hazard Analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula.” Tectonophysics 390(1–

501

4):141–58. Retrieved

502 503 504

ce

500

(http://linkinghub.elsevier.com/retrieve/pii/S0040195104002689).

Ramos, Noelynna T. and Hiroyuki Tsutsumi. 2010. “Evidence of Large Prehistoric Offshore Earthquakes Deduced from Uplifted Holocene Marine Terraces in

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

pt

Page 27 of 29

505

Pangasinan Province, Luzon Island, Philippines.” Tectonophysics 495(3–4):145–

506

58. Retrieved (http://dx.doi.org/10.1016/j.tecto.2010.08.007).

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1

508

Rangin, C. et al. 1999. “Cenozoic Pull-Apart Basins in Central Myanmar: The Trace of

pt

507

the Path of India along the Western Margin of Sundaland.” Terra Abstr 4.

Sadigh, K., C. Y. Chang, J. A. Egan, F. Makdisi, and R. R. Youngs. 1997. “Attenuation

510

Relationships for Shallow Crustal Earthquakes Based on California Strong Motion

511

Data.” Seismological Research Letters 68(1):180–89. Retrieved April 25, 2017

512

(http://courses.ce.metu.edu.tr/ce5603/wp-content/uploads/sites/25/2015/03/Sadigh-

513

et.al_.-1997.pdf).

515 516

SEDIA (Sabah Economic Development and Investment Authority) Annual report (2015): http://www.sedia.com.my/SEDIA_AR2015.pdf

Simons, W. J. F. et al. 2007. “A Decade of GPS in Southeast Asia: Resolving

an

514

us cri

509

517

Sundaland Motion and Boundaries.” Journal of Geophysical Research: Solid Earth

518

112(B6):n/a--n/a. Retrieved (http://dx.doi.org/10.1029/2005JB003868). SOKOLOV, V., BONJER, K.-P. & WENZEL, F. 2004. Accounting for site effect in probabilistic assessment of seismic hazard for Romania and Bucharest: a case of deep seismicity in Vrancea zone. Soil Dynamics and Earthquake Engineering, 24, 929-947.

524

Tan, Denis N. K. and J. M. Lamy. 1990. “Tectonic Evolution of the NW Sabah

dM

519 520 521 522 523

Continental Margin since the Late Eocene.” Geol. Soc. Malaysia, Bulletin 27:241–

526

60.

pte

525

527

Tjia, H. D. 2007. “Kundasang (Sabah) at the Intersection of Regional Fault Zones of

528

Quaternary Age.” Geological Society of Malaysia, Bulletin, 53 (Jun) 59–66.

530 531 532 533

Tongkul, F. 1993. “Tectonic Control on the Development of the Neogene Basins in Sabah, East Malaysia.” Bulletin of the Geological Society of Malaysia 33:95–103.

ce

529

Retrieved April 24, 2017 (http://www.gsm.org.my/products/702001-101020PDF.pdf).

Toro, Gabriel R., Norman A. Abrahamson, and John F. Schneider. 1997. “Model of

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 29

534

Strong Ground Motions from Earthquakes in Central and Eastern North America:

535

Best Estimates and Uncertainties.” Seismological Research Letters 68(1):41–57.

Page 29 of 29

537

Wah, Alexander Yan Sze. 2011. “Geological Assessments of the Earthquake Sources and Hazard in Malaysia.” Seminar Teknikal Gempabumi.

us cri

538 539

pt

Retrieved (http://srl.geoscienceworld.org/content/68/1/41).

536

Wells, Donald L. and Kevin J. Coppersmith. 1994. “New Empirical Relationships among

540

Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface

541

Displacement.” Bulletin of the Seismological Society of America 84(4):974–1002.

542

Wong, Ivan G. 2013. “How Big, How Bad, How Often: Are Extreme Events Accounted for in Modern Seismic Hazard Analyses?” Natural Hazards.

543 544

Youngs, R. R., S. .. J. Chiou, W. Silva, and J. R. Humphrey. 1997. “R. R. Youngs, S. -J. Chiou, W. J. Silva and J. R. Humphrey, 1997, Strong Ground Motion Attenuation

546

Relationships for Subduction Zone Earthquakes, , Vol. 68, No. 1, 58-73.”

547

Seismological Research Letters 68(1):58–73.

an

545

549

dM

548

.

550 551

554 555 556 557

ce

553

pte

552

Ac

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

AUTHOR SUBMITTED MANUSCRIPT - JGE-101511.R1