THE PROJECT FOR STUDY ON IMPROVEMENT OF BRIDGES THROUGH DISASTER MITIGATING MEASURES FOR LARGE SCALE EARTHQUAKES IN THE REPUBLIC OF THE PHILIPPINES

Transcription

1 THE REPUBLIC F THE PHILIPPINES DEPARTMENT F PUBLIC WRKS AND HIGHWAYS DPWH) THE PRJECT FR STUDY N IMPRVEMENT F BRIDGES THRUGH DISASTER MITIGATING MEASURES FR LARGE SCALE EARTHQUAKES IN THE REPUBLIC F THE PHILIPPINES FINAL REPRT APPENDIX 2-A GENERALIZED ACCELERATIN RESPNSE SPECTRA DEVELPMENT BY PRBABILISTIC SEISMIC HAZARD ANALISIS PSHA) DECEMBER 2013 JAPAN INTERNATINAL CPERATIN AGENCY JICA) CTI ENGINEERING INTERNATINAL C., LTD CHDAI C., LTD. NIPPN KEI C., LTD.

2

3 APPENDIX 2-A GENERALIZED ACCELERATIN RESPNSE SPECTRA DEVELPMENT BY PRBABILISTIC SEISMIC HAZARD ANALYSIS PSHA)

4

5 Seismic acceleration spectral maps for the Philippines to be incorporated into the draft BSDS for the following return periods: 1) 1,000 years; 2) 500 years; 3) 100 years; and 4) 50 years are generated. At each return period, spectral maps are developed for three key spectral acceleration parameters: 1) peak ground acceleration PGA); 2) spectral acceleration at sec; and spectral acceleration at 1. sec. The maps are developed for AASHT site class B equivalent to Vs30 = 760 m/s). The flow procedure is shown in Figure 2A-1. Active faults as presently identified by Phivolcs Philippine Institute of Volcanology and Seismology) are shown plotted in Figure 2A-2. Also shown plotted are instrumentally recorded earthquake events from 1907 to 2012 with magnitude greater than 4 and focal depth of less than 100 kms. which are compiled consisting of about 26,000+ events) from Phivolcs and ISC International Seismological Centre) websites into an earthquake catalog. The magnitude scale is homogenized in a common magnitude scale moment magnitude scale in this study for the reasons that moment magnitude does not suffer from saturation during large earthquakes; and is now the most commonly adopted in most ground motion estimation models that are presently being proposed. Declustering algorithm based on Gardner and Knopoff 1974) is applied to retain only independent main shocks into events as shown plotted in Figure 2A-3), removing aftershocks and foreshocks. Completeness analysis based on the method of Stepp 1972) is applied to the catalog to remove possible biases towards bigger events in subsequent regression analysis for temporal characterization of earthquake occurrences for each defined seismic source model since it is known that lower magnitude earthquake events had been under-reported in the early part of the instrumental era; and become less so with progressively improved instruments. Seismic source modeling consisting of fault models and background seismicity models are shown in Figure 2A-3. Background seismicity modeling is used to model seismic occurrences into areal zones where the observed seismicity exhibits a more or less diffused pattern that cannot be clearly identified with a specific fault. This may include earthquake occurrences in the future that could be attributable to blind thrusts or faults with no previous ground surface fault manifestations. Each earthquake event in the declustered set is identified to be associated with one of the fault models or background seismicity seismogenic areal zones. Bigger events are preferably made to be associated with the fault models. For source-to-site distance uncertainty modeling, earthquakes in this study are assumed to be uniformly distributed within a particular source zone i.e., earthquakes are considered equally likely to occur at any location within a source). Rupture may occur with equal likelihood anywhere in the fault plane in the fault zone and anywhere in the seismogenic areal zone. The spatial source-tosite distance) uncertainty can be described by a probability density function PR) which may be approximated by a normalized frequency distribution histogram. For fault models characterizing crustal earthquakes, maximum potential earthquake size capable to be produced within the source is computed using the empirical method of Wells and Coppersmith 1994). n the other hand, the method of Papazachos et al 2004) is used to compute maximum potential earthquake size for sources due to trenches. For seismogenic areal zones modeling back- 1

6 ground seismicity, the highest recorded or documented magnitude plus is used. List of historically documented earthquakes from 1589 to 1895 is based on the study by Bautista and ike 2000). Two types of earthquake recurrence models are used in this study: bounded Gutenberg-Richter recurrence model and characteristic earthquake recurrence model. The more commonly used bounded Gutenberg-Richter recurrence model in most PSHA implementation is expressed as: λ m = ν exp[ βm m 0)] exp[ βm max m 0 )] for m 0 m m max 1 exp[ βm max m 0 )] where m 0 is the lowest magnitude considered to be of engineering significance and m max is the maximum magnitude based on seismological and geological considerations as discussed earlier. Characteristic earthquake recurrence model using data based on paleoseimological observation is preferred but limited in use in this study due to the scarcity of data) due to the short history of instrumental recording in the world relative to geological period over which earthquakes recurred. Probabilistic seismic hazard analysis PSHA) provides a framework in which uncertainties in the size, location, and rate of recurrence of earthquakes and in the variation of ground motion characteristics with earthquake size and location can be identified, quantified, and combined in a rational manner Thenhaus and Campbell, 2003). The probability that an observed ground motion parameter X spectral acceleration, in this study) will be greater than or equal to the value x in the next t years the exposure period) given the annual exceedance rate λ [X x] is computed as: P [X x] =1 exp tλ[x x]) λ [X x] mmax v i sources i m 0 R M P [X x M,R] f M m) f R M r m) dr dm where λ [X x] the annual frequency that ground motion at a site exceeds the chosen level X = x; v i m 0 m max P [X x M,R] f M m) f R M r m) the annual rate of occurrence of earthquakes on seismic source i having magnitudes between m 0 and m max ; the minimum magnitude of engineering significance taken to be 5.0 in this study); the maximum magnitude assumed to occur on the source; the conditional probability that the chosen ground motion level is exceeded for a given magnitude M and distance R; probability density function of earthquake magnitude; probability density function of distance from the earthquake source to the site of interest. 2

7 Ground motion estimation models used in this study are based on Boore-Atkinson NGA 2007) applied to crustal earthquake sources, Young et al model 1997) applied to subduction sources; and Zhao et al 2006) applied to both crustal and subduction sources. The iterative analysis is carried out for a grid interval of 10 kms covering the whole Philippines for a total of 16,471 points. Interpolation and smoothing of the contours are made using the nearestneighbor algorithm. Results showing contours interposed with the source models and declustered earthquake events are shown in Figures 2A-4 to 2A-15. Contour maps of seismic acceleration values which constitute four corresponding to return periods of 1000 years, 500 years, 100 years, and 50 years) sets of maps at 3 key periods 0. sec, sec, and 1. sec) are finally generated: year return period a) peak ground acceleration PGA) Fig. 2A-16 and 17 regional maps in Figs. 2A-17 to 2A-33 b) spectral acceleration at sec S a at s) Fig. 2A-34 and 17 regional maps in Figs. 2A-35 to 2A-51 c) spectral acceleration at 1. sec S a at 1.s) Fig. 2A-52 and 17 regional maps in Figs. 2A-53 to 2A year return period a) peak ground acceleration PGA) Fig. 2A-70 and 17 regional maps in Figs. 2A-71 to 2A-87 b) spectral acceleration at sec S a at s) Fig. 2A-88 and 17 regional maps in Figs. 2A-89 to 2A-105 c) spectral acceleration at 1. sec S a at 1.s) Fig. 2A-106 and 17 regional maps in Figs. 2A-107 to 2A year return period a) peak ground acceleration PGA) Fig. 2A-124 and 17 regional maps in Figs. 2A-125 to 2A-141 b) spectral acceleration at sec S a at s) Fig. 2A-142 and 17 regional maps in Figs. 2A-143 to 2A-159 c) spectral acceleration at 1. sec S a at 1.s) Fig. 2A-160 and 17 regional maps in Figs. 2A-161 to 2A year return period a) peak ground acceleration PGA) Fig. 2A-178 and 17 regional maps in Figs. 2A-179 to 2A-195 b) spectral acceleration at sec S a at s) Fig. 2A-196 and 17 regional maps in Figs. 2A-197 to 2A-213 c) spectral acceleration at 1. sec S a at 1.s) Fig. 2A-214 and 17 regional maps in Figs. 2A-215 to 2A-231 3

8 The spectral acceleration contour maps by regions have been requested by DPWH for easier use of the bridge designers. References: D. M. Boore and G. M. Atkinson, Boore-Atkinson NGA ground motion relations for the geometric mean horizontal component of peak and spectral ground motion parameters, Technical Report PEER 2007/01, Pacific Earthquake Engineering Research Center, Maria Leonila P. Bautista and Kazuo ike, Estimation of the magnitudes and epicenters of Philippine historical earthquakes, Tectonophysics, 31712): , Maria Leonila P. Bautista and Bartolome C. Bautista, The Philippine historical earthquake catalog: its development, current state and future directions, Annals of Geophysics, 472/3), April/June J. K. Gardner and L. Knopoff, Is the sequence of earthquakes in Southern California, with aftershocks removed, poissonian? Bulletin of the Seismological Society of America, 645), Thomas H. Heaton, Fumiko Tajima, and Ann Wildenstein Mori, Estimating ground motions using recorded accelerograms, Surveys in Geophysics, 8:2583, A. R. Nelson, S. F. Personius, R. E. Rimando, R. S. Punongbayan, N. Tuñgol, H. Mirabueno, and A. Rasdas. Multiple large earthquakes in the past 1500 years on a fault in Metropolitan Manila, the Philippines. Bulletin of Seismological Society of America, 90:7385, B. C. Papazachos, E. M. Scordilis, D. G. Panagiotopoulos, C. B. Papazachos, and G. F. Karakaisis, Global relations between seismic fault parameters and moment magnitude of earthquakes, Bulletin of the Geological Society of Greece, XXXVI, L. Reiter. Earthquake Hazard Analysis: Issues and Insights. Columbia University Press, Carl J. Stepp, Analysis of completeness of the earthquake sample in the Puget Sound area and its effects on statistical estimates of earthquake hazard, In Proceedings of the International Conference on Microzonation for Safer Construction Research and Application, Seattle, ct 30 to Nov 3, 1972, vol. 2, P. C. Thenhaus and K. W. Campbell. Seismic Hazard Analysis, in Earthquake Engineering Handbook, edited by W.-F. Chen and C. Scawthorn, CRC Press, D. L. Wells and K. J. Coppersmith, New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement, Bulletin of the Seismological Society of America, 844): , August

9 R. R. Youngs, S.-J. Chiou, W. J. Silva, and J. R. Humphrey, Strong ground motion attenuation relationships for subduction zone earthquakes, Seismological Research Letters, 681), January/February J. X. Zhao, J. Zhang, A. Asano, Y. hno, T. ouchi, T. Takahashi, H. gawa, K. Irikura, H. K. Thio, P. G. Somerville, Y. Fukushima, and Y. Fukushima, Attenuation relations of strong ground motion in Japan using site classification based on predominant period, Bulletin of the Seismological Society of America, 963):898913, June

10 Figure 2A-1 Procedure of PSHA study for spectral mapping of PGA, S a at s and 1. s at base rock equivalent to AASHT site class B V s30 = 760 m/s) corresponding to return periods of 50, 100, 500, and 1000 years 6

11 Km PHILIPPINES LEGEND Active Faults - Trace Certain Approximate ffshore Projection Trace Approximate Trench Trough Depth Magnitude Earthquakes ) m m m Figure 2A-2 Seismological and Tectonic Setting of the Philippines instrumentally recorded earthquake events ; M w > 4, d<100 kms) 7

12 PHILIPPINES Seismic Source Modeling LEGEND Fault Models Background Seismicity Models Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-3 Seismic source modeling fault models and background seismicity models) for this PSHA study of the Philippines 8

13 PHILIPPINES PGA g) for 1000-yr Return Period LEGEND Fault Models Background Seismicity Models PGA g) for 1000-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-4 Contoured result of PGA for 1,000-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 9

14 PHILIPPINES SA g) at sec for 1000-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at sec for 1000-yr RP 0 0 0,70 0, , , Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-5 Contoured result of S a at sec. for 1,000-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 10

15 PHILIPPINES SA g) at 1 sec for 1000-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at 1 sec for 1000-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km 5 Figure 2A-6 Contoured result of S a at 1. sec. for 1,000-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 11

16 PHILIPPINES PGA g) for 500-yr Return Period LEGEND Fault Models Background Seismicity Models 5 PGA g) for 500-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-7 Contoured result of PGA for 500-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 12

17 PHILIPPINES SA g) at sec for 500-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at sec for 500-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km 1 Figure 2A-8 Contoured result of S a at sec. for 500-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 13

18 PHILIPPINES SA g) at 1 sec for 500-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at 1 sec for 500-yr RP ,35 0 0,45 Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-9 Contoured result of S a at 1. sec. for 500-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 14

19 PHILIPPINES PGA g) for 100-yr Return Period 0.12 LEGEND 0.04 Fault Models Background Seismicity Models 0.08 PGA g) for 100-yr RP Magnitude Declustered Earthquakes ) Depth m m m Km Figure 2A-10 Contoured result of PGA for 100-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 15

20 5 PHILIPPINES SA g) at sec for 100-yr Return Period LEGEND 0.1 Fault Models Background Seismicity Models SA g) at sec for 100-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-11 Contoured result of S a at sec. for 100-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 16

21 PHILIPPINES SA g) at 1 sec for 100-yr Return Period 0.1 LEGEND 0.06 Fault Models Background Seismicity Models SA g) at 1 sec for 100-yr RP Magnitude Declustered Earthquakes ) Depth m m m Km Figure 2A-12 Contoured result of S a at 1. sec. for 100-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 17

22 PHILIPPINES PGA g) for 50-yr Return Period LEGEND 0.02 Fault Models Background Seismicity Models PGA g) for 50-yr RP Magnitude Declustered Earthquakes ) Depth m m m Km Figure 2A-13 Contoured result of PGA for 50-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 18

23 PHILIPPINES SA g) at sec for 50-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at sec for 50-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km Figure 2A-14 Contoured result of S a at sec. for 50-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 19

24 PHILIPPINES SA g) at 1 sec for 50-yr Return Period LEGEND Fault Models Background Seismicity Models SA g) at 1 sec for 50-yr RP Magnitude Declustered Earthquakes ) Depth 0-30 m m m Km 0.05 Figure 2A-15 Contoured result of S a at 1. sec. for 50-year return period superimposed on seismic source models and declustered earthquake plots used in analysis 20

25 PHILIPPINES PGAg) for 1000-yr Return Period Km Figure 2A-16 Contour map of peak ground acceleration for 1,000-year return period 21

26 REGIN I PGAg) for 1000-yr Return Period 5 Ilocos Norte 5 Ilocos Sur 5 5 La Union Pangasinan Figure 2A-17 Region I map of peak ground acceleration for 1,000-year return period 22

27 REGIN II 5 PGAg) for 1000-yr Return Period 5 Cagayan 5 Isabela 5 Nueva Vizcaya Quirino Figure 2A-18 Region II map of peak ground acceleration for 1,000-year return period 23

28 CAR 5 PGAg) for 1000-yr Return Period 5 Apayao Abra Kalinga Mountain Province 5 5 Ifugao Benguet Figure 2A-19 CAR map of peak ground acceleration for 1,000-year return period 24

29 5 5 5 Aurora Nueva Ecija Tarlac Zambales Pampanga Bulacan Bataan Figure 2A-20 Region III map of peak ground acceleration for 1,000-year return period REGIN III PGAg) for 1000-yr Return Period 25

30 NCR PGAg) for 1000-yr Return Period Kalookan City Valenzuela Navotas Malabon Quezon City Marikina Manila San Juan Mandaluyong Pasig City Pasay City Makati City Pateros Taguig Parañaque Las Piñas Muntinlupa Km 5 Figure 2A-21 NCR map of peak ground acceleration for 1,000-year return period 26

31 5 5 Quezon Rizal 5 Cavite Laguna Batangas Quezon 5 Figure 2A-22 Region IV-A map of peak ground acceleration for 1,000-year return period REGIN IVA PGAg) for 1000-yr Return Period 27

32 5 5 5 Camarines Norte 5 Catanduanes 5 Camarines Sur Albay 5 Sorsogon Masbate Figure 2A-23 Region V map of peak ground acceleration for 1,000-year return period REGIN V PGAg) for 1000-yr Return Period 28

33 REGIN IVB PGAg) for 1000-yr Return Period 5 5 Marinduque riental Mindoro ccidental Mindoro Romblon 5 Palawan Km Figure 2A-24 Region IV-B map of peak ground acceleration for 1,000-year return period 29

34 REGIN VI PGAg) for 1000-yr Return Period 5 5 Aklan Capiz Antique Iloilo Guimaras Negros ccidental Figure 2A-25 Region VI map of peak ground acceleration for 1,000-year return period 30

35 REGIN VII PGAg) for 1000-yr Return Period Cebu Bohol Negros riental 5 5 Siquijor Figure 2A-26 Region VII map of peak ground acceleration for 1,000-year return period 31

36 REGIN VIII PGAg) for 1000-yr Return Period Northern Samar 5 5 Samar Biliran 5 Eastern Samar 5 Leyte Southern Leyte 5 Figure 2A-27 Region VIII map of peak ground acceleration for 1,000-year return period 32

37 5 REGIN IX PGAg) for 1000-yr Return Period 5 5 Zamboanga del Norte Zamboanga Sibugay Zamboanga del Sur Figure 2A-28 Region IX map of peak ground acceleration for 1,000-year return period 33

38 REGIN X PGAg) for 1000-yr Return Period Camiguin 5 Misamis ccidental Lanao del Norte Misamis riental Bukidnon 5 5 Figure 2A-29 Region X map of peak ground acceleration for 1,000-year return period 34

39 5 5 REGIN XI PGAg) for 1000-yr Return Period Davao del Norte Compostela Valley Davao riental Davao del Sur 5 5 Figure 2A-30 Region XI map of peak ground acceleration for 1,000-year return period 35

40 5 REGIN XIII PGAg) for 1000-yr Return Period 5 Dinagat Islands Surigao del Norte Agusan del Norte 5 Surigao del Sur Agusan del Sur Figure 2A-31 Region XIII map of peak ground acceleration for 1,000-year return period 36

41 ARMM PGAg) for 1000-yr Return Period Lanao del Sur Shariff Kabunsuan Maguindanao Figure 2A-32 ARMM map of peak ground acceleration for 1,000-year return period 37

42 REGIN XII PGAg) for 1000-yr Return Period North Cotabato Sultan Kudarat South Cotabato Sarangani Figure 2A-33 Region XII map of peak ground acceleration for 1,000-year return period 38

43 PHILIPPINES Spectral Acceleration g) at sec 1000-yr Return Period , Km Figure 2A-34 Contour map of spectral acceleration at sec. for 1,000-year return period 39

44 REGIN I 0.85 SA g) at.2sec for 1000-yr Return Period Ilocos Norte Ilocos Sur La Union 1.05 Pangasinan 1.1 Figure 2A-35 Region I map of spectral acceleration at sec. for 1,000-year return period 40

45 REGIN II SA g) at.2sec for 1000-yr Return Period Cagayan Isabela 1.4 Nueva Vizcaya Quirino Figure 2A-36 Region II map of spectral acceleration at sec. for 1,000-year return period 41

46 CAR 0.75 SA g) at.2sec for 1000-yr Return Period Apayao Abra Kalinga 5 Mountain Province Ifugao Benguet Figure 2A-37 CAR map of spectral acceleration at sec. for 1,000-year return period 42

47 Aurora Nueva Ecija Tarlac Zambales Pampanga Bulacan Bataan REGIN III SAg) at.2sec for 1000-yr Return Period Figure 2A-38 Region III map of spectral acceleration at sec. for 1,000-year return period 43

48 NCR SA g) at.2sec for 1000-yr Return Period Kalookan City 1.45 Valenzuela 1.4 Navotas Malabon Quezon City Marikina Manila San Juan Mandaluyong Pasig City 1.35 Pasay City Makati City Pateros Taguig Parañaque Las Piñas 1.25 Muntinlupa Km Figure 2A-39 NCR map of spectral acceleration at sec. for 1,000-year return period 44

49 1 REGIN IVA SAg) at.2sec for 1000-yr Return Period 0.85 Quezon 1.05 Rizal Cavite Laguna Batangas Quezon Figure 2A-40 Region IV-A map of spectral acceleration at sec. for 1,000-year return period 45

50 Camarines Norte Catanduanes Camarines Sur Albay Sorsogon 0.7 Masbate REGIN V SAg) at.2sec for 1000-yr Return Period 1.15 Figure 2A-41 Region V map of spectral acceleration at sec. for 1,000-year return period 46

51 1 REGIN IVB SA g) at.2sec for 1000-yr Return Period Marinduque riental Mindoro ccidental Mindoro Romblon Palawan Km 1.5 Figure 2A-42 Region IV-B map of spectral acceleration at sec. for 1,000-year return period 47

52 1 0.7 REGIN VI SA g) at.2sec for 1000-yr Return Period Aklan Capiz Antique Iloilo 0.75 Guimaras Negros ccidental Figure 2A-43 Region VI map of spectral acceleration at sec. for 1,000-year return period 48

53 1.3 REGIN VII SA g) at.2sec for 1000-yr Return Period Cebu Bohol 1.5 Negros riental Siquijor Figure 2A-44 Region VII map of spectral acceleration at sec. for 1,000-year return period 49

54 REGIN VIII SAg) at.2sec for 1000-yr Return Period Northern Samar Eastern Samar Samar Biliran Leyte Southern Leyte Figure 2A-45 Region VIII map of spectral acceleration at sec. for 1,000-year return period 50

55 REGIN IX SAg) at.2sec for 1000-yr Return Period Zamboanga del Norte Zamboanga Sibugay Zamboanga del Sur Figure 2A-46 Region IX map of spectral acceleration at sec. for 1,000-year return period 51

56 1 REGIN X SAg) at.2sec for 1000-yr Return Period Camiguin Misamis ccidental 0.7 Lanao del Norte Misamis riental Bukidnon Figure 2A-47 Region X map of spectral acceleration at sec. for 1,000-year return period 52

57 REGIN XI SA g) at.2sec for 1000-yr Return Period Davao del Norte Compostela Valley 1.1 Davao riental Davao del Sur Km 1.5 Figure 2A-48 Region XI map of spectral acceleration at sec. for 1,000-year return period 53

58 REGIN XIII SAg) at.2sec for 1000-yr Return Period 0.8 Dinagat Islands Surigao del Norte Agusan del Norte Surigao del Sur 0.85 Agusan del Sur Figure 2A-49 Region XIII map of spectral acceleration at sec. for 1,000-year return period 54

59 ARMM SAg) at.2sec for 1000-yr Return Period 0.7 Lanao del Sur Shariff Kabunsuan Maguindanao Figure 2A-50 ARMM map of spectral acceleration at sec. for 1,000-year return period 55

60 REGIN XII SA g) at.2sec for 1000-yr Return Period North Cotabato Sultan Kudarat South Cotabato Sarangani Figure 2A-51 Region XII map of spectral acceleration at sec. for 1,000-year return period 56

61 PHILIPPINES Spectral Acceleration g) at 1 sec 1000-yr Return Period Km Figure 2A-52 Contour map of spectral acceleration at 1. sec. for 1,000-year return period 57

62 5 5 REGIN I SA g) at 1sec for 1000-yr Return Period 5 Ilocos Norte Ilocos Sur 5 5 La Union Pangasinan Figure 2A-53 Region I map of spectral acceleration at 1. sec. for 1,000-year return period 58

63 REGIN II SA g) at 1sec for 1000-yr Return Period 5 5 Cagayan 5 5 Isabela Nueva Vizcaya Quirino 5 5 Figure 2A-54 Region II map of spectral acceleration at 1. sec. for 1,000-year return period 59

64 5 CAR SA g) at 1sec for 1000-yr Return Period 5 Apayao Abra 5 Kalinga Mountain Province Ifugao 5 Benguet 5 Figure 2A-55 CAR map of spectral acceleration at 1. sec. for 1,000-year return period 60

65 Aurora Nueva Ecija Tarlac Zambales 5 Pampanga Bulacan Bataan Figure 2A-56 Region III map of spectral acceleration at 1. sec. for 1,000-year return period REGIN III SAg) at 1sec for 1000-yr Return Period 61

66 NCR SA g) at 1sec for 1000-yr Return Period Kalookan City 5 Valenzuela Navotas Malabon Quezon City Marikina Manila San Juan Mandaluyong Pasig City Pasay City Makati City Pateros Taguig Parañaque Las Piñas 5 Muntinlupa Km Figure 2A-57 NCR map of spectral acceleration at 1. sec. for 1,000-year return period 62

67 Quezon Rizal Cavite Laguna 5 Batangas Quezon 5 5 REGIN IVA SAg) at 1sec for 1000-yr Return Period 5 Figure 2A-58 Region IV-A map of spectral acceleration at 1. sec. for 1,000-year return period 63

68 Camarines Norte Catanduanes Camarines Sur 5 Albay Sorsogon Masbate 5 Figure 2A-59 Region V map of spectral acceleration at 1. sec. for 1,000-year return period REGIN V SAg) at 1sec for 1000-yr Return Period 64

69 5 REGIN IVB SA g) at 1sec for 1000-yr Return Period Marinduque 5 riental Mindoro ccidental Mindoro 5 Romblon 5 Palawan Km Figure 2A-60 Region IV-B map of spectral acceleration at 1. sec. for 1,000-year return period 65

70 REGIN VI SA g) at 1sec for 1000-yr Return Period 5 Aklan Capiz Antique 5 Iloilo Guimaras Negros ccidental Figure 2A-61 Region VI map of spectral acceleration at 1. sec. for 1,000-year return period 66

71 REGIN VII SA g) at 1sec for 1000-yr Return Period Cebu Bohol 5 5 Negros riental 5 Siquijor Figure 2A-62 Region VII map of spectral acceleration at 1. sec. for 1,000-year return period 67

72 REGIN VIII SAg) at 1sec for 1000-yr Return Period Northern Samar 5 Samar 5 5 Biliran 5 Eastern Samar Leyte 5 Southern Leyte Figure 2A-63 Region VIII map of spectral acceleration at 1. sec. for 1,000-year return period 68

73 5 5 5 REGIN IX SAg) at 1sec for 1000-yr Return Period Zamboanga del Norte 5 5 Zamboanga Sibugay Zamboanga del Sur 5 5 Figure 2A-64 Region IX map of spectral acceleration at 1. sec. for 1,000-year return period 69

74 REGIN X SAg) at 1sec for 1000-yr Return Period Camiguin 5 Misamis ccidental Lanao del Norte Misamis riental Bukidnon 5 Figure 2A-65 Region X map of spectral acceleration at 1. sec. for 1,000-year return period 70

75 5 REGIN XI SA g) at 1sec for 1000-yr Return Period Davao del Norte Compostela Valley Davao riental Davao del Sur Km Figure 2A-66 Region XI map of spectral acceleration at 1. sec. for 1,000-year return period 71

76 5 5 REGIN XIII SAg) at 1sec for 1000-yr Return Period Dinagat Islands Surigao del Norte 5 Agusan del Norte Surigao del Sur Agusan del Sur Figure 2A-67 Region XIII map of spectral acceleration at 1. sec. for 1,000-year return period 72

77 ARMM SAg) at 1sec for 1000-yr Return Period Lanao del Sur Shariff Kabunsuan 5 Maguindanao 5 5 Figure 2A-68 ARMM map of spectral acceleration at 1. sec. for 1,000-year return period 73

78 REGIN XII SA g) at 1sec for 1000-yr Return Period 5 North Cotabato Sultan Kudarat 5 South Cotabato Sarangani 5 Figure 2A-69 Region XII map of spectral acceleration at 1. sec. for 1,000-year return period 74

79 PHILIPPINES PGA g) for 500-yr Return Period , Km 5 Figure 2A-70 Contour map of peak ground acceleration for 500-year return period 75

80 REGIN I PGAg) for 500-yr Return Period Ilocos Norte 5 Ilocos Sur 5 La Union Pangasinan 5 Figure 2A-71 Region I map of peak ground acceleration for 500-year return period 76

81 5 REGIN II PGAg) for 500-yr Return Period 5 5 Cagayan Isabela Nueva Vizcaya Quirino Figure 2A-72 Region II map of peak ground acceleration for 500-year return period 77

82 5 CAR PGAg) for 500-yr Return Period 5 5 Apayao Abra Kalinga 5 Mountain Province Ifugao Benguet Figure 2A-73 CAR map of peak ground acceleration for 500-year return period 78

83 Aurora Nueva Ecija Tarlac Zambales Pampanga Bulacan Bataan Figure 2A-74 Region III map of peak ground acceleration for 500-year return period REGIN III PGAg) for 500-yr Return Period 79

84 NCR PGAg) for 500-yr Return Period 5 Kalookan City Valenzuela Navotas Malabon Quezon City Marikina Manila San Juan Mandaluyong Pasig City 5 5 Pasay City Makati City Pateros Taguig Parañaque Las Piñas Muntinlupa Km 5 Figure 2A-75 NCR map of peak ground acceleration for 500-year return period 80

85 5 5 5 Quezon Rizal 5 Cavite Laguna 5 Batangas Quezon 5 REGIN IVA PGAg) for 500-yr Return Period Figure 2A-76 Region IV-A map of peak ground acceleration for 500-year return period 81

86 5 5 Camarines Norte 5 Catanduanes 5 Camarines Sur Albay Sorsogon 5 Masbate 5 5 Figure 2A-77 Region V map of peak ground acceleration for 500-year return period REGIN V PGAg) for 500-yr Return Period 5 82

87 5 REGIN IVB PGAg) for 500-yr Return Period riental Mindoro ccidental Mindoro Marinduque 5 Romblon 5 Palawan Km 5 Figure 2A-78 Region IV-B map of peak ground acceleration for 500-year return period 83

88 5 REGIN VI PGAg) for 500-yr Return Period 5 Aklan Capiz Antique Iloilo Guimaras Negros ccidental 5 5 Figure 2A-79 Region VI map of peak ground acceleration for 500-year return period 84

89 REGIN VII PGAg) for 500-yr Return Period 5 Cebu Bohol Negros riental 5 Siquijor 5 5 Figure 2A-80 Region VII map of peak ground acceleration for 500-year return period 85

90 5 5 REGIN VIII PGAg) for 500-yr Return Period Northern Samar 5 Samar Biliran Eastern Samar 5 Leyte 5 5 Southern Leyte Figure 2A-81 Region VIII map of peak ground acceleration for 500-year return period 86

91 5 5 REGIN IX PGAg) for 500-yr Return Period 5 Zamboanga del Norte Zamboanga Sibugay Zamboanga del Sur 5 Figure 2A-82 Region IX map of peak ground acceleration for 500-year return period 87

92 REGIN X PGAg) for 500-yr Return Period 5 Misamis ccidental Lanao del Norte Camiguin Misamis riental Bukidnon 5 5 Figure 2A-83 Region X map of peak ground acceleration for 500-year return period 88

93 5 5 REGIN XI PGAg) for 500-yr Return Period Davao del Norte Compostela Valley 5 Davao riental Davao del Sur 5 5 Figure 2A-84 Region XI map of peak ground acceleration for 500-year return period 89

94 REGIN XIII PGAg) for 500-yr Return Period Dinagat Islands Surigao del Norte Agusan del Norte Surigao del Sur 5 Agusan del Sur 5 5 Figure 2A-85 Region XIII map of peak ground acceleration for 500-year return period 90

95 ARMM PGAg) for 500-yr Return Period 5 Lanao del Sur Shariff Kabunsuan Maguindanao 5 5 Figure 2A-86 ARMM map of peak ground acceleration for 500-year return period 91

96 REGIN XII PGAg) for 500-yr Return Period North Cotabato Sultan Kudarat South Cotabato 5 Sarangani Figure 2A-87 Region XII map of peak ground acceleration for 500-year return period 92