The Philippine archipelago is tectonically complex and seismically hazardous, yet few seismic hazard assessments have provided national coverage. This article presents an updated probabilistic seismic hazard analysis for the nation. Active shallow crustal seismicity is modeled by faults and gridded point sources accounting for spatially variable occurrence rates. Subduction interfaces are modeled with faults of complex geometry. Intraslab seismicity is modeled by ruptures confined to the slab volume. Source geometries and earthquake rates are derived from seismicity catalogs, geophysical data sets, and historic-to-paleoseismic constraints on fault slip rates. The ground motion characterization includes models designed for global use, with partial constraint by residual analysis. Shallow crustal faulting near metropolitan Manila, Davao, and Cebu dominates shaking hazard. In a few places, peak ground acceleration with 10% probability of exceedance in 50 years on rock reaches 1.0g. The results of this study may have utility for defining the design base shear in the National Structural Code of the Philippines.
AbrahamsonNAGregorNAddoK (2016) BC Hydro ground motion prediction equations for subduction earthquakes. Earthquake Spectra32(1): 23–44.
2.
AbrahamsonNASilvaWJKamaiR (2014) Summary of the ASK14 ground motion relation for active crustal regions. Earthquake Spectra30(3): 1025–1055.
3.
AkkarSBommerJJ (2010) Empirical equations for the prediction of PGA, PGV, and spectral accelerations in Europe, the Mediterranean region, and the Middle East. Seismological Research Letters81(2): 195–206.
4.
AllenTIRyuHBautistaBBautistaMLNaragISevillaWIMelosantosMLPPapionaKBonitaJ (2014) Enhancing Risk Analysis Capacities for Flood, Tropical Cyclone Severe Wind and Earthquake for the Greater Metro Manila Area (Component 5). Canberra, ACT, Australia: Geoscience Australia.
5.
AmbraseysNNDouglasJSarmaSKSmitPM (2005) Equations for the estimation of strong ground motions from shallow crustal earthquakes using data from Europe and the Middle East: Horizontal peak ground acceleration and spectral acceleration. Bulletin of Earthquake Engineering3(1): 1–53.
6.
American Society of Civil Engineers/Structural Engineering Institute (ASCE/SEI) (2013) Minimum Design Loads and Associated Criteria for Buildings and Other Structures (ASCE/SEI 7-16). Reston, VA: ASCE.
7.
AtkinsonGMAdamsJ (2013) Ground motion prediction equations for application to the 2015 Canadian national seismic hazard maps. Canadian Journal of Civil Engineering40: 988–998.
8.
AtkinsonGMBooreDM (2003) Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions. Bulletin of the Seismological Society of America93(4): 1703–1729.
9.
AurelioMA (2000) Shear partitioning in the Philippines: Constraints from Philippine fault and global positioning system data. Island Arc9(4): 584–597.
10.
BacolcolTC (2003) Etude geodesique de la faille Philippine dans les Visayas. Unpublished Doctoral Dissertation, Universite Pierre et Marie Curie, Paris.
11.
BacolcolTCBarrierEDuquesnoyTAguilarAJorgioRdelaCruzRLasalaM (2005) GPS constraints on Philippine fault slip rate in Masbate Island, central Philippines. Journal of the Geological Society of the Philippines60: 1–7.
12.
BacolcolTCSaplaGPelicanoALuisAJrJorgioRMarfitoBParasMRagadioZBalotoAAbsinO (2017) Interseismic velocities and GPS motions associated with 6 July 2017 M6.5 Philippine fault (Leyte Segment) earthquake. PHIVOLCS open-file report no. 18–02, 10 pp. Quezon City, Philippines: Philippine Institute of Volcanology and Seismology.
13.
BarrierEHuchonPAurelioM (1991) Philippine fault: A key for Philippine kinematics. Geology19: 32–35.
14.
BautistaMLPOikeK (2000) Estimation of the magnitudes and epicenters of Philippine historical earthquakes. Tectonophysics317: 137–169.
15.
BesanaGMAndoM (2005) The central Philippine Fault Zone: Location of great earthquakes, slow events and creep activity. Earth Planets and Space57(10): 987–994.
16.
BindiDCottonFKothaSRBosseCStromeyerDGrünthalG (2017) Application-driven ground motion prediction equation for seismic hazard assessments in non-cratonic moderate-seismicity areas. Journal of Seismology21(5): 1201–1218.
17.
BooreDMAtkinsonGM (2008) Ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5%-damped PSA at spectral periods between 0.01 s and 10.0 s. Earthquake Spectra24(1): 99–138.
18.
BooreDMStewartJPSeyhanEAtkinsonGM (2014) NGA-West2 equations for predicting PGA, PGV, and 5% damped PSA for shallow crustal earthquakes. Earthquake Spectra30(3): 1057–1085.
19.
CampbellKWBozorgniaY (2014) NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra. Earthquake Spectra30(3): 1087–1115.
20.
ChiouBSJYoungsRR (2008) An NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra24(1): 173–215.
21.
ChiouBSJYoungsRR (2014) Update of the Chiou and Youngs NGA model for the average horizontal component of peak ground motion and response spectra. Earthquake Spectra30(3): 1117–1153.
22.
DaligdigJA (1997) Recent faulting and paleoseismicity along the Philippine fault zone, north central Luzon, Philippines. Unpublished Doctoral Dissertation, Kyoto University, Kyoto, Japan, 73 pp.
23.
DziewonskiAMChouTAWoodhouseJH (1981) Determination of earthquake source parameters from waveform data for studies of global and regional seismicity. Journal of Geophysical Research: Solid Earth86: 2825–2852.
24.
EkströmGNettlesMDziewonskiAM (2012) The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes. Physics of the Earth and Planetary Interiors200–201: 1–9.
25.
FarrTGRosenPACaroECrippenRDurenRHensleySKobrickMPallerMRodriguezERothLSealDShafferSShimadaJUmlandJWernerMOskinMBurbankDAlsdorfD (2007) The shuttle radar topography mission. Reviews of Geophysics45(2): RG2004.
26.
FieldEHArrowsmithRJBiasiGPBirdPDawsonTEFelzerKRJacksonDDJohnsonKMJordanTHMaddenCMichaelAJMilnerKRPageMTParsonsTPowersPMShawBEThatcherWRWeldonRJZengY (2014) Uniform California earthquake rupture forecast, version 3 (UCERF3)—the time-independent model. Bulletin of the Seismological Society of America104: 1122–1180.
27.
FieldEHDawsonTEFelzerKRFrankelADGuptaVJordanTHParsonsTPetersenMDSteinRSWeldonRJWillsCJ (2009) Uniform California earthquake rupture forecast, version 2 (UCERF 2). Bulletin of the Seismological Society of America99: 2053–2107.
28.
FitchTJ (1972) Plate convergence, transcurrent faults, and internal deformation adjacent to Southeast Asia and the western Pacific. Journal of Geophysical Research77(23): 4432–4460.
29.
FletcherJMTeranOJRockwellTKOskinMEHudnutKWMuellerKJSpelzRMAkcizSOMasanaEFanerosGFieldingEJLeprinceSMorelanAEStockJLynchDKElliottAJGoldPLiu-ZengJGonzález-OrtegaAHinojosa-CoronaAGonzález-GarcíaJ (2014) Assembly of a large earthquake from a complex fault system: Surface rupture kinematics of the 4 April 2010 El Mayor–Cucapah (Mexico) Mw 7.2 earthquake. Geosphere10(4): 797–827.
30.
FrankelA (1995) Mapping seismic hazard in the central and eastern United States. Seismological Research Letters66(4): 8–21.
31.
FukushimaYTanakaT (1990) A new attenuation relation for peak horizontal acceleration of strong earthquake ground motion in Japan. Bulletin of the Seismological Society of America80(4): 757–783.
32.
GalganaGHamburgerMMcCaffreyRCorpuzEChenQ (2007) Analysis of crustal deformation in Luzon, Philippines using geodetic observations and earthquake focal mechanisms. Tectonophysics432(1–4): 63–87.
33.
GarciaCLValenzuelaRArnoldEPMacalinagTGAmbubuyogGFLanceNTCordetaJDDoniegoAGDabiACBalceGRSuFS (1985) Series on Seismology: Philippines. In: Southeast Asia Association of Seismology and Earthquake Engineering, vol. 4, pp. 743–792. US Geological Survey-USAID Singapore.
34.
GardnerJKKnopoffL (1974) Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian?Bulletin of the Seismological Society of America64(5): 1363–1367.
35.
General Bathymetric Chart of the Oceans (GEBCO) (2008) The GEBCO_08 Grid. International Hydrographic Organization and the Intergovernmental Oceanographic Commission of UNESCO, p. 19. Available at: http://www.gebco.net/
36.
GutenbergBRichterCF (1944) Frequency of earthquakes in California. Bulletin of the Seismological Society of America34: 185–188.
37.
HamlingIJHreinsdóttirSClarkKElliottJLiangCFieldingELitchfieldNVillamorPWallaceLWrightTJD’AnastasioEBannisterSBurbidgeDDenysPGentlePHowarthJMuellerCPalmerNPearsonCPowerWBarnesPBarrellDJVan DissenRLangridgeRLittleTNicolAPettingaJRowlandJStirlingM (2017) Complex multifault rupture during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand. Science356(6334): eaam7194.
38.
HayesGPMooreGLPortnerDEHearneMFlammeHFurtneyMSmoczykGM (2018) Slab2, a comprehensive subduction zone geometry model. Science362(6410): 58–61.
39.
HsuYJYuSBLovelessJPBacolcolTSolidumRLuisAJrPelicanoAWoessnerJ (2016) Interseismic deformation and moment deficit along the Manila subduction zone and the Philippine Fault system. Journal of Geophysical Research: Solid Earth121(10): 7639–7665.
40.
HsuYJYuSBSongTRABacolcolT (2012) Plate coupling along the Manila subduction zone between Taiwan and northern Luzon. Journal of Asian Earth Sciences51: 98–108.
41.
International Code Council (ICC) (2018) International Residential Code. Country Club Hills, IL: ICC.
42.
KannoTNaritaAMorikawaNFujiwaraHFukushimaY (2006) A new attenuation relation for strong ground motion in Japan based on recorded data. Bulletin of the Seismological Society of America96(3): 879–897.
43.
KaverinaANLanderAVProzorovAG (1996) Global creepex distribution and its relation to earthquake-source geometry and tectonic origin. Geophysical Journal International125(1): 249–265.
44.
LaskeGMastersGTMaZPasyanosM (2013) Update on CRUST1.0—A 1-Degree Global Model of Earth’s Crust, vol. 15 (Geophysical research abstracts). Vienna: EGU General Assembly.
45.
LeonardM (2010) Earthquake fault scaling: Self-consistent relating of rupture length, width, average displacement, and moment release. Bulletin of the Seismological Society of America100(5A): 1971–1988.
46.
MolasGLYamazakiF (1994) Seismic macrozonation of the Philippines based on seismic hazard analysis. Journal of Structural Mechanics and Earthquake Engineering11: 33s–43s.
47.
National Structural Code of the Philippines (NSCP) (2015) NSCP C101-15, vol. 1. 7th ed.Manila, Philippines: Association of Structural Engineers of the Philippines (ASEP).
48.
PaganiMMMarcelliniA (2007) Seismic-hazard disaggregation: A fully probabilistic methodology. Bulletin of the Seismological Society of America97(5): 1688–1701.
49.
PaganiMMMonelliDWeatherillGDanciuLCrowleyHSilvaVHenshawPButlerLNastasiMPanzeriLSimionatoMViganoD (2014) OpenQuake engine: An open hazard (and risk) software for the Global Earthquake Model. Seismological Research Letters85(3): 692–702.
50.
PapionaKLKinugasaY (2008) Trenching surveys along the Masbate segment of the Philippine fault zone. Journal of the Geological Society of the Philippines64: 19–53.
51.
PasyanosMEMastersGTLaskeGMaZ (2014) LITHO1.0: An updated crust and lithospheric model of the Earth. Journal of Geophysical Research: Solid Earth119(3): 2153–2173.
52.
PeñarubiaHC (2017) Seismic coefficient exceedance in the recent destructive earthquake events in the Philippines. In: Proceedings of the 18th Association of Structural Engineers of the Philippines (ASEP) international convention, Quezon City, Philippines, 25–27 May.
53.
PeñarubiaHCGrutasRNDeocampoJB (2017) Probabilistic seismic hazard analysis of the Philippines. In: Proceedings of the 18th Association of Structural Engineers of the Philippines (ASEP) international convention, Quezon City, Philippines, 25–27 May.
54.
PerezJSTsutsumiH (2011) Paleoseismic studies along the Philippine fault zone, eastern Mindanao, Philippines. In: Proceedings of the symposium on Japan Geoscience Union meeting, Chiba, Japan, 22–27 May.
55.
PerezJSTsutsumiH (2017) Tectonic geomorphology and paleoseismology of the Surigao segment of the Philippine fault in northeastern Mindanao Island, Philippines. Tectonophysics699: 244–257.
56.
PerezJSTsutsumiHCahuloganMTCabanlitDPAbiganiaMITNakataT (2015) Fault distribution, segmentation and earthquake generation potential of the Philippine fault in eastern Mindanao, Philippines. Journal of Disaster Research10(1): 74–82.
57.
Philippine Institute of Volcanology and Seismology (PHIVOLCS) (2017) The Philippine Earthquake Model: A Probabilistic Seismic Hazard Assessment of the Philippines and of Metro Manila. Available at: https://www.phivolcs.dost.gov.ph/images/PEM.pdf
58.
PubellierMBaderAGRanginCDeffontainesBQuebralR (1999) Upper plate deformation induced by subduction of a volcanic arc: The Snellius Plateau (Molucca Sea, Indonesia and Mindanao, Philippines). Tectonophysics304(4): 345–368.
59.
RanginC (2016) Rigid and non-rigid micro-plates: Philippines and Myanmar-Andaman case studies. Comptes Rendus Geoscience348(1): 33–41.
60.
RanginCLe PichonXMazzottiSPubellierMChamot-RookeNAurelioMWalpersdorfAQuebralR (1999) Plate convergence measured by GPS across the Sundaland/Philippine sea plate deformed boundary: The Philippines and eastern Indonesia. Geophysical Journal International139(2): 296–316.
61.
RongYPaganiMMMagistraleHWeatherillG (2017) Modeling seismic hazard by integrating historical earthquake, fault, and strain rate data. In: Proceedings of the 16th world conference on earthquake engineering, Santiago, 9–13 January, paper no. 448.
62.
SadighKChangCYEganJAMakdisiFYoungsRR (1997) Attenuation relationships for shallow crustal earthquakes based on California strong motion data. Seismological Research Letters68(1): 180–189.
63.
StiphoutTVan ZhuangJMarsanD (2012) Theme V—models and techniques for analysing seismicity. Technical report, Community Online Resource for Statistical Seismicity Analysis, p. 22. Available at: http://www.corssa.org
64.
SuSS (1988) Seismic hazard analysis for the Philippines. Natural Hazards1: 27–44.
65.
ThenhausPCHansonSLAlgermissenSTBautistaBCBautistaMLPunongbayanBJPunongbayanRS (1994) Estimates of the Regional Ground Motion Hazard of the Philippines. National Disaster Mitigation of the Philippines, pp. 45–60.
66.
ThingbaijamKKSMaiPMGodaK (2017) New empirical earthquake source-scaling laws. Bulletin of the Seismological Society of America107(5): 2225–2246.
67.
TorregosaRFSugitoMNojimaN (2002) Assessment of seismic hazard and microzoning in the Philippines. Journal of Structural Mechanics and Earthquake Engineering19(2): 81s–98s.
68.
TsutsumiHPerezJS (2013) Large-scale digital mapping of the Philippine fault zone based on aerial photograph interpretation. Active Faults Research39: 29–37.
69.
TsutsumiHDaligdigJAGotoHTungolNMKondoHNakataTOkunoMSugitoN (2006) Timing of surface-rupturing earthquakes on the Philippine fault zone in Central Luzon Island, Philippines. EOS Transactions American Geophysical Union87: 52.
70.
TsutsumiHPerezJSMarjesJUPapionaKLRamosNT (2015) Coseismic displacement and recurrence interval of the 1973 Ragay Gulf earthquake, southern Luzon, Philippines. Journal of Disaster Research10(1): 83–90.
71.
TurcotteDLSchubertG (2014) Geodynamics. Cambridge: Cambridge University Press, 636 pp.
72.
UhrhammerRA (1985) Seismicity record, M> 2.5, for the central. Clement F Shearer, p. 199
73.
WeatherillGA (2014) Openquake Ground Motion Toolkit—User Guide (Technical report). Pavia: Global Earthquake Model (GEM) Foundation.
74.
WeatherillGAPaganiMMGarciaJ (2016) Exploring earthquake databases for the creation of magnitude-homogeneous catalogues: Tools for application on a regional and global scale. Geophysical Journal International206(3): 1652–1676.
75.
WeichertDH (1980) Estimation of the earthquake recurrence parameters for unequal observation periods for different magnitudes. Bulletin of the Seismological Society of America70(4): 1337–1346.
76.
WellsDLCoppersmithKJ (1994) New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America84(4): 974–1002.
77.
WesselPSmithWHFScharrooRLuisJWobbeF (2013) Generic mapping tools: Improved version released. Eos Transactions American Geophysical Union94: 409–410.
78.
YoungsRRCoppersmithKJ (1985) Implications of fault slip rates and earthquake recurrence models to probabilistic seismic hazard estimates. Bulletin of the Seismological Society of America75: 939–964.
79.
YoungsRRChiouSJSilvaWJHumphreyJR (1997) Strong ground motion attenuation relationships for subduction zone earthquakes. Seismological Research Letters68(1): 58–73.
80.
YuS-BKuoL-CPunongbayanRSRamosEG (1999) GPS observation of crustal deformation in the Taiwan-Luzon region. Geophysical Research Letters26(7): 923–926.
81.
YumulGPDimalantaCBTamayoRA (2005) Indenter-tectonics in the Philippines: Example from the Palawan microcontinental block—Philippine mobile belt collision. Resource Geology55(3): 189–198.
82.
YumulGPJrDimalantaCBTamayoRAJrMauryRC (2003) Collision, subduction and accretion events in the Philippines: A synthesis. Island Arc12(2): 77–91.
83.
ZhangQGuoFZhaoLWuY (2017) Geodynamics of divergent double subduction: 3-D numerical modeling of a Cenozoic example in the Molucca Sea region, Indonesia. Journal of Geophysical Research: Solid Earth122(5): 3977–3998.
84.
ZhaoJXZhouSZhouJZhaoCZhangHZhangYGaoPLanXRhoadesDFukushimaYSomervillePGIrikuraK (2016) Ground-motion prediction equations for shallow crustal and upper-mantle earthquakes in Japan using site class and simple geometric attenuation functions. Bulletin of the Seismological Society of America106(4): 1552–1569.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
