Mati City Site Response Atlas
Mati City Site Response Atlas
SITE RESPONSE
ATLAS
2021
Coastal urban areas in Mati City, Davao Oriental. (Photo taken on November 2020)
Cagayan River in Cauayan City, Isabela
Published by the Department of Science and Technology - Philippine Institute of Volcanology and Seismology
(DOST-PHIVOLCS)
ISBN 978-621-96559-3-4
Any part of this document may be used and reproduced as presented provided proper acknowledgement is made.
https://www.phivolcs.dost.gov.ph
Front Cover Photo: Aerial view of the baywalk in Mati City, Davao Oriental. (Photo taken on November 2020)
Congratulations to DOST-PHIVOLCS and DOST-PCIEERD for the successful implementation and funding of the project titled:
"Specific Earthquake Ground Motion Levels to Help Increase the Seismic Resiliency of Government Infrastructures, Residential,
and Medium-to-High Rise Buildings in Pangasinan, Tarlac, Metro Iloilo-Guimaras, Cauayan City, Butuan City, and Mati City",
also known as the Specific Earthquake Project. I would like to commend their grit and passion towards the development and
publication of the Mati City Response Atlas.
This Atlas is another addition to the many disaster mitigation and risk reduction tools produced by DOST-PHIVOLCS. The Atlas will
help local goverment units of concerned study areas to verify areas that are vulnerable to various seismic hazards, and to advise
stakeholders on the parameters to consider in their building designs. This will also provide significant help in the land use planning
and zoning efforts in Mati City, by adjusting the existing development plans with the observed ground motion levels within an
area. Considering the local ground conditions, this Atlas aims to identify areas suitable and safe for infrastructural development.
May the information contained in this atlas benefit all end users in responding to earthquake-related challenges. Let the knowledge
start with us, and let us aim to create more disaster-resilient communities.
I would like to thank DOST and DOST-PCIEERD for their continuous support in the completion of the Specific Earthquake Project.
With their help, we were able to develop and proudly publish the Mati City Site Response Atlas. Finally, I would like to acknowledge
and congratulate our very own, Specific Earthquake Project Team for their hard work and dedication towards the fulfillment and
publication of this Atlas.
This Atlas is an effective means to enlighten everyone of us, and help us prepare for possible future earthquakes by presenting
to us the relevant information essential for increasing the seismic resiliency of structures in Mati City. The Atlas showcases the
following hazard maps generated for this city namely the: Peak Ground Acceleration (PGA) Map, Vs30 Model Map, Short-Period,
and Long-Period Microzonation Map. Both the Vs30 Model Map and PGA Map are widely used by engineers and geologists
to identify the soil profile type on site, and to calculate the design base shear as part of the minimum design requirements in
building structures, respectively. On the other hand, the Short-Period and Long-Period Microzonation Maps can be used as guide
for building designs to avoid soil-structure resonance, which commonly leads to structural damage and collapse. Moreover, with
this new data, it would be possible to identify existing structures that may be in need of further retrofitting schemes due to
potential building resonance. This technology will be a valuable supplemental data for urban planning, zoning, and disaster risk
reduction and management endeavors.
With this Atlas, we are confident that we have surely strengthened not only our structural resilience but also the intellectual
resilience of Filipinos, taking one step at a time for our country. Mabuhay!
The City of Mati would like to congratulate the Department of Science and Technology - Philippine Institute of Volcanology and
Seismology (DOST - PHIVOLCS) for the achievement of the DOST-PCIEERD funded project titled "Specific Earthquake Ground Motion
Levels to Help Increase the Seismic Resiliency of Government Infrastructures, Residential, and Medium-to-High Rise Buildings in
Pangasinan, Tarlac, Metro Iloilo-Guimaras, Cauayan City, Butuan City, and Mati City". We are greatly honored to be one of the
target smart cities of this project, and we are grateful to be provided with this significant data at hand.
We believe that the success of this project will substantially benefit us by providing the necessary information to our citizens
about ground motion levels in the City of Mati. The Mati City Site Response Atlas will be very useful to aspects like urban planning
and zoning, city engineering, planning and development, as well as to our disaster risk reduction and management endeavors.
This will ultimately enhance the resiliency of future infrastructures and development, which in turn will also strengthen the will
of the people of Mati to ground shaking-related adversities with the knowledge and wisdom this atlas brings.
On behalf of the City of Mati, I would like to extend my grattitude to DOST-PHIVOLCS and DOST-PCIEERD helping us in our goal to
achieve a more disasater-resilient city.
I would like to commend the scientists, engineers, and staff of DOST-PHIVOLCS and DOST-PCIEERD for their effort in the development
of the Site Response Atlases of six (6) selected cities/provinces. This is a way forward in the field of earthquake engineering
practice in the Philippines, and can therefore help increase the seismic resiliency and support efficient engineering design of
public infrastructures.
The next generation of engineering codes such as the DPWH Bridge Seismic Design Specifications and the upcoming revision of
the ASEP-NSCP for buildings, adopted by the Bureau of Design, will certainly benefit by adopting input parameters and datasets
from these Atlases. The shift of various codes worldwide towards full adoption of seismic soil parameter maps such as Vs30 map
in the design of buildings and bridges as embodied by the Site Response Atlas, has also laid out the groundwork for a national
scale mapping.
In line with this, let me express this Bureau's full support in further development of Site Response Atlases.
We would like to acknowledge and congratulate the commited scientists and engineers of DOST-PHIVOLCS that worked for the
completion and publication of the six (6) Site Response Atlases of the Specfic Earthquake Project, and the DOST-PCIEERD for the
funding assistance and guidance throughout the implementation of this scientific project.
With the publication of this Atlas, engineers and urban planners of the local government units of cities and municipalities involved
in the atlas can incorporate the results as part of their minimum design requirements and seismic regulations, during the planning
and development phases of infrastructures and urban designs. Further, the information gathered and the methodologies used in
the completion of site response atlases can be utilized as inputs in the amendment of the National Structural Code of the Philippines
(NSCP). With the initiative of this project, we can continue to pave the way to further develop and refine its methodologies of
estimating the ground motion response levels within an area, which is very important, timely, and relevant for a seismically-active
country in its developing stage. If we continue these efforts, we can further improve its current potential in increasing the seismic
resiliency of not only our present urban situation, but also in shaping our future's urban development.
May the accompanied wisdom embedded in this Atlas be used in its utmost potential by the local governments units, stakeholders,
private sectors, and other end users to help us in this campaign of continuously building sustainable communities at the same
time spread and the disaster-resilient mind set and outlook.
Definition of Terms stiff mountain range, and its proximity to earthquake generators, Mati City may exhibit
vulnerability to ground shaking and the effects of site response.
• Active Fault – a fault that has moved within the last 10,000 years. It shows
evidence or has documented history of its recent movement. One of the notable examples of this effect was demonstrated last 1985 in Mexico
• Long-Period – period of ground motion greater than 1 second (band under 1 City when an earthquake with a moment magnitude of 8.0, and an epicenter located
Hertz; low frequency). more than 300 kilometers away rattled the city. As a huge portion of the city was
sitting on top of ancient lake sediments, ground shaking was amplified. This resulted
• Microtremors – one of the passive/natural fields of the earth. They are to huge property damage and casualties (Cassaro & Romero, 1986). Ground shaking
ubiquitous, weak, low-amplitude background vibrations or ambient noises amplification occurs when the period of ground motion matches the natural period
generated from natural phenomena and anthropogenic sources that could be or resonance of the structures on site called soil-structure resonance. The given
recorded on the surface of the earth (Okada, 2003). structures will experience the largest oscillations possible and are more likely to sustain
• Peak Ground Acceleration (PGA) – the maximum ground acceleration that structural damage and collapse due to prolonged and amplified ground shaking.
occurred during an earthquake ground shaking at a given location.
The DOST-PCIEERD funded project titled “Specific Earthquake Ground Motion Levels to
• Seismic shear-wave – the body wave that is commonly perceived as a horizontal Help Increase the Seismic Resiliency of Government Infrastructures, Residential, and
vibration perpendicular to the direction of the wave propagation during an Medium-to-High Rise Buildings in Pangasinan, Tarlac, Metro Iloilo-Guimaras, Cauayan
earthquake. City, Butuan City, and Mati City”, also known as the “Specific Earthquake Project”
• Short-Period – period of ground motion less than 1 second (band over 1 Hertz; of DOST-PHIVOLCS aims to determine site response using the three (3) geophysical
high frequency). methods and processing techniques to generate the three (3) different seismic
microzonation maps. Important details of the survey methodologies are presented in
• Site Response – the influence of local site conditions to alter earthquake Figure 1 and Table 1.
ground motion on a particular site.
• Soil-Structure Resonance – accentuated amplitude of vibration or oscillation of
the structure during an earthquake excitation occurring when its predominant
period approximately matches with the predominant period of the ground.
• Trench – the deepest portion on the sea floor. It is a manifestation of subduction
zone where a tectonic plate moved or is pushed under another tectonic plate.
• Vs30 – the average seismic shear-wave velocity of the upper 30 meters of soil
layer used for characterizing site conditions.
Overview
Mati City is located on the southeastern part of Mindanao Island and the capital of the Figure 1. Specific Earthquake Project Conceptual Flowchart.
Province of Davao Oriental in Region XI (Davao Region). According to the Philippine
Statistics Authority (2021), Mati City has a population of about 147,547 in 2020, making
it the 4th most populous city in the whole region (excluding the independent, highly Table 1. Specific Earthquake Project Method Descriptions.
urbanized Davao City). The growing population in the city results to rapid urbanization Data
Frequency Sampling
and increase in the demand for infrastructural developments. Method Instrument
Response
Acquisition
Frequency
Target Ground Property
Time
Mati City is considered one of the seismically active areas in the country due to its Refraction
DAQLink III
4.5 Hertz 500
Average shear-wave
proximity to several earthquake generators that could induce short-period and long- Microtremor
Seismograph;
(vertical 15 minutes samples
velocity of the upper 30
12-Channel R.T. Clark meters of soil layer (Vs30)
period ground motions. The nearest potential earthquake source is the Mati Segment (ReMi) Survey
Vertical Geophones
geophones) per second
in meter per second (m/s)
of the Philippine Fault traversing the city (pp. 4) (Perez et. al., 2015; Tsutsumi & Perez, Three- OYO McSEIS-MT
100 Amplification at higher
2013). Another is an unnamed offshore fault about 30 kilometers southeast of the city Component Neo Seismometer 0.1 to 200 30 to 50
samples frequency (short-period)
Microtremor (3-channel <Vertical, Hertz minutes
proper. The Philippine Trench is also an active source that can be found more than a Survey Horizontal>)
per second in second (s)
hundred kilometers in the offshore east of the city (pp. 3). Microtremor Seven (7) OYO McSEIS- 100 Amplification at lower
0.1 to 200 40 to 60
Array MT Neo Seismometers samples frequency (long-period) in
Aside from the presence of earthquake generators, site response influences the Measurements (1-channel <Vertical>)
Hertz minutes
per second second (s)
amplitude and duration of ground shaking regardless of its distance to the epicenter
of the earthquake, and dictates the amplification at specific frequencies during an
earthquake. This phenomenon is primarily anchored to the subsurface geology of All surveys were performed under fair weather conditions starting from the 13th of
the area. It is usually manifested by geomorphological features or differences in November to 2nd of December 2020, from 8 AM to 5 PM. Survey sites were distributed
topography (e.g., sedimentary basin versus mountain range/basement). Stiffness widely across the city to ensure good spatial resolution. Gain selection depends
and thickness of layers also play a significant role in the amplification of the ground. on the local site conditions during the survey, since all three (3) survey methods
Decrease in the stiffness and increase in the thickness of a layer (e.g., soil layer) result measure microtremors. The recordings were conducted at quiet environments where
to a decrease in the frequency. Thus, high impedance contrast across the boundary microtremors as a stationary process dominate in the data (e.g., avoiding areas with
of a hard rock and an overlying thick, soft soil layer result to lower frequencies being heavy traffic or industrial/mechanical noise and vibrations from construction sites,
amplified in the soil layer with respect to the hard rock. This is commonly exemplified etc.).
in a sedimentary basin, colloquially known as a bowl of hard rock filled with soft soils
– an amplifying geologic structure. Vs30 Model Map
Mati City is located in a coastal area where most of the city proper and urban area Ground response varies per area, depending on the physio-chemical properties of the
are sitting on the coastal plains fringing the southern tip of the Eastern Mindanao rock or soil layers influencing its stiffness. Thickness, on the other hand, also varies
Ridge (pp. 5). These coastal areas are flat tract of land form by accumulation of from place to place but is commonly dependent on the geology and topography of the
sediments supplied by rivers coming from the nearby mountain range. Several area. Soils or sediments in sedimentary basin or alluvial valleys are commonly thicker
creeks traverse the city and flow towards the bay. Two bays are found in Mati City compared to soils in mountain ranges and other elevated areas. Because of Mati City’s
separated by the Guanguan Peninsula, namely the Pujada Bay and Mayo Bay. Mati geologic setting and topography, there is wide range of soil thickness and variable
City’s downtown urban area, being part of a coastal depositional environment, where stiffness. Thus, obtaining the average shear-wave velocity of the top 30 meters of soil
soft soils are relatively thick, has the tendency to alter or amplify ground motion at known as “Vs30” is necessary, which is one of the main parameters influencing ground
specific frequencies. Thus, because of the city’s topography and subsurface geology motion and amplification.
brought about by its geographical position in a coastal area at the southern tip of a
1
Vs30 values were obtained using the Refraction Microtremor (ReMi) Survey as shown change in lithology have also corresponded to a decrease in short-period values and
in the flowchart in Figure 1 and in Table 1. Across various topographies in Mati City, increase in the Vs30 values. Based on these observations, the geology and topography
27 sites were surveyed using the ReMi method. The Vs30 values were interpolated, and generally affect the period of ground motion as well as the corresponding stiffness or
the resulting map was integrated with the Vs30 map produced from topographic slope Vs30 value.
data using Interferometric Synthetic Aperture Radar (IfSAR) (Wald & Allen, 2007), to
finally generate the Vs30 Model Map of Mati City. Vs30 ranges were also categorized as
site classes in the legend of the map in pp. 7 using Table 2. Long-Period Microzonation Map
Since Mati City is situated in a generally erosional landscape at the foot of the mountains,
Table 2. Vs30 site classification for seismic site response as defined by the National the effects of long-period ground shaking cannot be significantly estimated, compared
Earthquake Hazards Reduction Program (NEHRP, 1994). to a depositional one (e.g. sedimentary basins). Long-period zones delineated using
the method are generally limited to areas that does not significantly vary in elevation.
Site Class Range of Vs30 (m/s) General Description Long-period values were obtained using the Microtremor Array Measurements
performed at 5 sites across the narrow and low-lying coastal plains (occasionally
A >1500 Hard rock
B 760 – 1500 Rock transitions to gently sloping piedmont areas) in Mati City, as shown in the flowchart
C 360 – 760 Very dense soil and soft rock in Figure 1 and in Table 1. The long-period values (in seconds) were interpolated to
D 180 – 360 Stiff soil produce the Long-Period Microzonation Map of Mati City (pp. 9). Based on the results,
E <180 Soil profile long-period zones in Mati City are varied and are divided into three zones that are
roughly localized in the built-up, urban areas of the city. The highest estimated long-
period values (1.51 to 2.00 seconds) in this study coincide with the highest short-
Based on the results, most of the areas in Mati City are dominated by the Site Classes period ranges and lowest Vs30 values (Site Class E). This area is observed to be an
C and D. Generally, Vs30 values increase with elevation. This means that the lowest alluvial valley (with some piedmont area included) and is one of the locations in the
Vs30 values of 180 to 240 m/s pertaining to Site Class D (with small portions of <180 city’s coastal urban area with extensive Quaternary Alluvium (Qal) deposits (pp.5).
m/s, Site Class E) dominate the flat topographies pertaining to the coastal bay areas It is also in this area where the Mati Segment of the Philippine Fault can be traced.
at the foot of the Eastern Mindanao Ridge composed of young and less-consolidated Long-period values of 1.00 to 1.50 seconds roughly coincide and delineate the middle
sediments known as Quaternary Alluvium (pp. 5). These areas also include alluvial ranges of the short-period and Vs30 values. Lastly, the lowest estimated long-period
valleys, river floodplains, and the lowlands of the Guanguan Peninsula. On the other values in this study that roughly covers the flat Qal area of the Guanguan Peninsula
hand, areas classified as Site Class C (Vs30 values of 490 to 760 m/s) correspond to the coincide with one of the lowest short-period values (0.09 to 0.15 seconds) with also
relatively older, rigid, and resistant mountainous terrains of the Eastern Mindanao lower Vs30 values (180 to 240 m/s). This suggests the existence of soft, but shallow or
Ridge (northern Mati City) composed of magmatic arc rock suites and volcaniclastic thin soil layers beneath this tract of lowlands in the peninsula. The setting of this area
sequences. This highest Vs30 ranges estimated in the study also include the ophiolitic does not significantly promote thick deposition of soft sediments since it is generally
rocks of the Pujada Peninsula (southern Mati City), and the uplands at the southern tip an eroding feature. The varying or shifting conditions (e.g. waves and tidal effects) and
of the Guanguan Peninsula. These predominantly Site Class C areas may also include the absence of rivers/creeks in this area that supply sediments may have also affected
areas with Vs30 values classified under Site Class B (>760 m/s) which were not detailed the deposition of soils in this area.
in the study. Lowest Vs30 values in this study with less than 180 m/s (Site Class E) are
found in some localized coastal areas immediately facing Pujada Bay. Another site
with Vs30 value of this range is also localized around the area where the mouth of the Application and Limitation
Mayo River meets Mayo Bay. Design engineers can use the Vs30 model map to identify the soil or rock type in an area
and its stiffness based on the corresponding Vs30 range. Zoning the predominant period
Short-Period Microzonation Map of ground motion is also essential in developing earthquake-resilient infrastructures,
especially for residential and low-rise buildings in the city. In addition, the long-period
Short-period values were obtained using the Three-Component Microtremor Survey microzonation maps are used to verify areas suitable or critical for the development
performed at 70 sites across various topographies in Mati City (Figure 1 and Table 1). of medium-to-high rise buildings. The period maps can be generally used by
The short-period values (in seconds) were interpolated to produce the Short-Period planners or engineers in making sure that the natural period of infrastructures to be
Microzonation Map of Mati City (pp. 8). The resulting map correlates with the Vs30 developed does not coincide with the observed period of ground motion. Generally,
model map. Higher ranges of short-period values dominate areas in the city with the these microzonation maps can be used as a supplementary data for seismic load
lowest relief composed of Quaternary Alluvium located mostly in the coastal areas of design considerations for future structures, as well as in urban planning and zoning.
the city facing Pujada Bay and Mayo Bay, whereas the lower ranges are found in the Moreover, concerned Local Government Units (LGUs) could adapt and incorporate the
mountainous terrains in the northern parts of the city associated with the Eastern results of this atlas into their respective Comprehensive Land Use Plans (CLUPs) and
Mindanao Ridge, the hilly lands at the southern tip of the Guanguan Peninsula, and Comprehensive Development Plans (CDPs), or as a resolution or policy, that in turn
the entire stretch of the mountainous Pujada Peninsula in the southern parts of the would also benefit the private sectors, researchers, stakeholders, and other end-users
city (pp. 5). It is worthy to note the similarities of both maps in some coastal areas to build a more earthquake-resilient and disaster-prepared city.
and alluvial valleys within the city proper directly facing Pujada Bay. The highest short-
period values in this study of more than 0.70 seconds roughly coincide with the lowest The maps included in this atlas are semi-detailed and are produced from regional-
Vs30 values (<180 m/s) pertaining to Site Class E. This suggests that thicker and softer scale assessments, wherein data boundaries are approximate and gradational. Hence,
soils or sediments underlie these coastal areas. There is also a notable similarity in the detailed and site-specific evaluation is recommended for construction of major
coastal areas facing Mayo Bay, wherein lower Vs30 values dominated by Site Class D structures and lifelines, such as borehole drilling and auxiliary geophysical surveys. The
(180 to 240 m/s) also coincide with higher short-period values of 0.36 to 0.4 seconds. results reflected in the maps may be used as preliminary data and may be augmented
The piedmont areas of the city where flatlands generally start to increase slope and by detailed study or through other conventional survey methods.
References
• Cassaro, M. A., and Romero, E. M. (1986). The Mexico earthquakes 1985: factors involved and lesson learned. Proceedings of the International Conference, Mexico
City, September 1986. American Society of Civil Engineers
• MGB, Mines and Geosciences Bureau [Cartographer]. (2021). Geologic Map of Mati City - Compilation of Quadrangle Maps from 1992 to 1995 [Map]. 1:120,000.
Version 1. Quezon City, Philippines: DOST-PHIVOLCS.
• NEHRP, National Earthquake Hazards Reduction Program (1994). Recommended Provisions for Seismic Regulations for New Buildings: Part I – Provisions, FEMA 222A,
Federal Agency Management Agency, Washington, D.C.
• Okada, H. (2003). The Microtremor Survey Method, Geophysical Monograph Series. Society of Exploration Geophysicists.
• Perez, J. S., Tsutsumi, H., Cahulogan, M. T., Cabanlit, D. P., Abigania, M. I., and Nakata, T. (2015). Fault Distribution, Segmentation and Earthquake Generation
Potential of the Philippine Fault in Eastern Mindanao, Philippines. Journal of Disaster Research. 10(1), 74-82
• PHIVOLCS, Philippine Institute of Volcanology and Seismology [Cartographer]. (2018). Active Faults Map of the Province of Davao Oriental [Map]. 1:250,000. Version
1. Quezon City, Philippines: DOST-PHIVOLCS.
• PHIVOLCS, Philippine Institute of Volcanology and Seismology [Cartographer]. (2021). Distribution of Active Faults and Trenches in the Philippines [Map]. Quezon City,
Philippines: DOST-PHIVOLCS.
• PHIVOLCS, Philippine Institute of Volcanology and Seismology [Cartographer]. (2021). Peak Ground Acceleration Map of Mati City [Map]. 1:120,000. Version 1.
Quezon City, Philippines: DOST-PHIVOLCS.
• Philippine Statistics Authority (2021). Highlights of the Region XI (Davao Region) Population 2020 Census of Population and Housing (2020 CPH). Accessed on:
September 2021. Retrieved from: https://psa.gov.ph/population-and-housing/node/165015
• Tsutsumi, H., and Perez, J. S. (2013). Large-scale active fault map of the Philippine fault based on aerial photograph interpretation. Active Fault Research, 39, 29-37
• Wald, D. J. and Allen, T. I. (2007). Topographic Slope as a Proxy for Seismic Site Conditions and Amplification. Bulletin of the Seismological Society of America, 97(5),
1379-1395. doi:10.1785/0120060267
2
3
4
5
6
7
8
9
Survey Methods
B8-Center
10
Several Destructive Earthquakes in Mindanao, Philippines
DOST-PHIVOLCS DOST-PHIVOLCS
A reinforced three-storey commercial building in A two-storey house in Brgy. Poblacion, Magsaysay, Davao del Sur
Brgy. Poblacion, Padada, Davao del Sur had a pancake collapse. damaged due to ground shaking.
DOST-PHIVOLCS DOST-PHIVOLCS
Eva’s Building in Kidapawan City, Cotabato damaged Daig Elementary School in Tulunan, Cotabato collapsed
during the 31 October 2019 M6.5 earthquake event. during the 29 October 2019 M6.6 earthquake event.
DOST-PHIVOLCS DOST-PHIVOLCS
Collapsed Anao-aon Bridge in San Francisco, Surigao del Norte. Damaged San Nicolas High School in San Francisco, Surigao del Norte.
MATI CITY
SITE
RESPONSE
ATLAS
Department of Science and Technology
PHILIPPINE INSTITUTE OF VOLCANOLOGY AND SEISMOLOGY
PHIVOLCS Building, C.P. Garcia Avenue, U.P. Campus
Diliman, Quezon City 1101
https://www.phivolcs.dost.gov.ph