Earthquake model successfully forecasts location of August 3 magnitude 6.8 event

The earthquake that struck offshore eastern Mindanao in the southern Philippines on August 3, 2024 resulted from renewed movement of an extinct fracture zone subducting into Earth’s mantle. That movement along the subduction zone was responsible for a magnitude 7.4 event that occurred eight months earlier.
 

By Mario Aurelio and Sandra Donna Catugas, Structural Geology and Tectonics Laboratory, University of Philippines National Institute of Geological Sciences, Philippines
 

Citation: Aurelio, M. and Catugas, S. D., 2024, Earthquake model successfully forecasts location of August 3 magnitude 6.8 event, Temblor, http://doi.org/10.32858/temblor.349
 

On August 3, 2024 at 6:23 a.m. local time, a magnitude 6.8 earthquake struck offshore, east of central Mindanao in the southern Philippines. The Philippine Institute of Volcanology and Seismology (PHIVOLCS) located the event about 48 kilometers east-northeast of the town of Lingig (population approximately 35,000) in the province of Surigao del Sur.

The quake ruptured at a depth of 25 kilometers (16 miles). Within four days, this mainshock was succeeded by at least six aftershocks with magnitudes between 6.2 and 4.8, all at shallower depths. Fortunately, no major infrastructure damage or casualties were reported.
 

 Figure 1a. Seismotectonics of the magnitude 6.8 earthquake of August 3, 2024 in the province of Surigao Del Sur, Mindanao (southern Philippines) in relation to the magnitude 7.4 earthquake of December 3, 2023 in the same province. The result of Coulomb Stress Transfer (CST) modeling is shown. In this model, the fault that generated the magnitude 7.4 earthquake of December 3, 2023 serves as source fault (thrust mechanism) and the fault that generated the August 3, 2024 earthquake serves as receiver fault (thrust fault mechanism). Focal mechanism solutions of the August 2024 mainshock and six succeeding aftershocks with magnitudes ranging from magnitude 4.8 to 6.2, and other post-mainshock events of lower magnitudes (small circles) are shown. Location of cross section in Figure 1b is shown (A-A’). Data: Jarvis et al., 2008 for SRTM topography; Weatherall et al., 2020 for bathymetry; Toda et al., 2011 for CST modeling; PHIVOLCS for earthquake data. GMT (Wessel and Smith, 1995) was used to generate the map. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0

Figure 1a. Seismotectonics of the magnitude 6.8 earthquake of August 3, 2024 in the province of Surigao Del Sur, Mindanao (southern Philippines) in relation to the magnitude 7.4 earthquake of December 3, 2023 in the same province. The result of Coulomb Stress Transfer (CST) modeling is shown. In this model, the fault that generated the magnitude 7.4 earthquake of December 3, 2023 serves as source fault (thrust mechanism) and the fault that generated the August 3, 2024 earthquake serves as receiver fault (thrust fault mechanism). Focal mechanism solutions of the August 2024 mainshock and six succeeding aftershocks with magnitudes ranging from magnitude 4.8 to 6.2, and other post-mainshock events of lower magnitudes (small circles) are shown. Location of cross section in Figure 1b is shown (A-A’). Data: Jarvis et al., 2008 for SRTM topography; Weatherall et al., 2020 for bathymetry; Toda et al., 2011 for CST modeling; PHIVOLCS for earthquake data. GMT (Wessel and Smith, 1995) was used to generate the map. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0

 

Seismotectonic setting

The left-lateral Philippine Fault, which stretches more than 1,200 kilometers in distance, has several southern segments that traverse Eastern Mindanao (Aurelio, 1992; Quebral, 1994) (Figure 1a). Offshore and to the east, the Philippine Trench marks the westward subduction of the Philippine Sea Plate underneath the Philippine Mobile Belt (Cardwell et al., 1980; Aurelio, 2000a; Aurelio and Peña, 2010). The Philippine Fault–Philippine Trench pair operates under a shear partitioning mechanism, where the oblique convergence between the Philippine Sea Plate and the Philippine Mobile Belt (Rangin et al., 1999) can be divided into westward frontal subduction along the Philippine Trench and left-lateral slip along the Philippine Fault (Aurelio, 2000b). Both tectonic features are responsible for generating numerous earthquakes in the region, many of which have been devastating.
 

Figure 1b. Cross section showing the vertical distribution of earthquakes, including the magnitude 7.4 event of December 3, 2024, and the August 3, 2024 mainshock (magnitude 6.8) along with six succeeding aftershocks (with magnitudes between 4.8 and 6.2). Focal mechanism solution plots are shown for the main events, whereas small pale dots indicate historical earthquakes. Approximate trace of the subducting Philippine Sea Plate is shown by a dashed line. Data sources same as in Figure 1a. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0
Figure 1b. Cross section showing the vertical distribution of earthquakes, including the magnitude 7.4 event of December 3, 2024, and the August 3, 2024 mainshock (magnitude 6.8) along with six succeeding aftershocks (with magnitudes between 4.8 and 6.2). Focal mechanism solution plots are shown for the main events, whereas small pale dots indicate historical earthquakes. Approximate trace of the subducting Philippine Sea Plate is shown by a dashed line. Data sources same as in Figure 1a. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0

 

The subducting Philippine Sea Plate enters the Philippine Trench at a slight westerly angle (between 15 and 30 degrees). The plate begins to steepen between depths of 10 and 20 kilometers, eventually reaching its steepest angle of about 45 degrees at around 50 kilometers depth. The slab can be detected to a depth of more than 150 kilometers in some places (Cardwell et al., 1980; Aurelio, 2000) (Figure 1b). The 25-kilometer depth of the magnitude 6.8 earthquake places this event at the end of the gently dipping section of the subducting slab at almost the same depth as the magnitude 7.4 earthquake of December 3, 2023 (Aurelio et al., 2023).
 

Triggered earthquakes

We explored the relationship between the December 2023 magnitude 7.4 earthquake and the August 2024 magnitude 6.8 event using a Coulomb Stress Transfer (CST) model. All CST models require a source fault and a receiver fault. In our model, the fault that generated the magnitude 7.4 event served as the source fault. The fault that generated the magnitude 6.8 earthquake served as the receiver fault. Red areas represent zones of increased stress, whereas blue areas represent zones of decreased stress (Figures 1a and 1b). The focus of the August 3 event is located in a red zone, suggesting that this event was triggered by the magnitude 7.4 event.

However, in the same CST model (Figures 1a and 1b), none of the aftershocks of the magnitude 6.8 earthquake plot in the stress increase (red) zones. This may suggest that the magnitude 7.4 earthquake of December 3, 2023 did not influence the generation and distribution of the August 2024 aftershocks.

In search of a more coherent explanation of the distribution of the August 2024 aftershocks, we created another CST model. In this case, the fault that generated the magnitude 6.8 event served as the source fault. The fault that generated the largest aftershock (magnitude 6.2) served as the receiver fault (Figures 2a and 2b).
 

Figure 2a. Seismotectonics of the magnitude 6.8 earthquake of August 3, 2024, in the province of Surigao Del Sur, Mindanao (southern Philippines) and its aftershocks. Here, we show the result of CST modeling in which the fault that generated the mainshock serves as the source fault (thrust mechanism) and the fault that generated the strongest aftershock (magnitude 6.2) serves as the receiver fault (thrust mechanism). We also show the focal mechanism solutions of the mainshock and six succeeding aftershocks with magnitudes ranging from magnitude 4.8 to 6.2, and other post-mainshock events of lower magnitudes (small circles). Location of cross section in Figure 2b is shown (B-B’). Data: Jarvis et al., 2008 for SRTM topography; Weatherall et al., 2020 for bathymetry; Toda et al., 2011 for CST modeling; PHIVOLCS for earthquake data. GMT (Wessel and Smith, 1995) was used to generate the map. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0
Figure 2a. Seismotectonics of the magnitude 6.8 earthquake of August 3, 2024, in the province of Surigao Del Sur, Mindanao (southern Philippines) and its aftershocks. Here, we show the result of CST modeling in which the fault that generated the mainshock serves as the source fault (thrust mechanism) and the fault that generated the strongest aftershock (magnitude 6.2) serves as the receiver fault (thrust mechanism). We also show the focal mechanism solutions of the mainshock and six succeeding aftershocks with magnitudes ranging from magnitude 4.8 to 6.2, and other post-mainshock events of lower magnitudes (small circles). Location of cross section in Figure 2b is shown (B-B’). Data: Jarvis et al., 2008 for SRTM topography; Weatherall et al., 2020 for bathymetry; Toda et al., 2011 for CST modeling; PHIVOLCS for earthquake data. GMT (Wessel and Smith, 1995) was used to generate the map. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0

 

In this model, we observe that all aftershocks have focal depths that are shallower than the mainshock; all aftershocks plot within the stress increase (red) zones. Five of the aftershocks, including the largest (magnitude 6.2) plot east and above the mainshock, whereas only one plots to the west. This strongly suggests that the aftershocks were generated by the release of energy mostly along the upper, gently dipping segments of the subduction thrust, whereas one (magnitude 5.4) aftershock may have been associated with an east-dipping backthrust on the overriding plate.
 

Figure 2b. Cross section showing the vertical distribution of the August 2024 earthquakes, including the mainshock (magnitude 6.8 of August 3, 2024) and six aftershocks (with magnitudes between 4.8 and 6.2) up to five days after the mainshock. Focal mechanism solution plots are shown for the main events; small pale dots indicate smaller-magnitude aftershocks. Approximate trace of the subducting Philippine Sea Plate is shown by a dashed line. Here, seismicity doesn’t clearly show that the slab subducts to a depth of 150 kilometers, as it does elsewhere. Data sources same as in Figure 2a. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0
Figure 2b. Cross section showing the vertical distribution of the August 2024 earthquakes, including the mainshock (magnitude 6.8 of August 3, 2024) and six aftershocks (with magnitudes between 4.8 and 6.2) up to five days after the mainshock. Focal mechanism solution plots are shown for the main events; small pale dots indicate smaller-magnitude aftershocks. Approximate trace of the subducting Philippine Sea Plate is shown by a dashed line. Here, seismicity doesn’t clearly show that the slab subducts to a depth of 150 kilometers, as it does elsewhere. Data sources same as in Figure 2a. Credit: Aurelio and Catugas, 2024, CC BY-NC-ND 4.0

 

Predictive power of CST modeling

The CST models presented above demonstrate the predictive power of this modeling technique, especially when properly taken within the context of complex tectonic settings such as those found in eastern Mindanao. Based on our modeling, the recent magnitude 6.8 mainshock occurred in a region in which stress increased as a result of the December 2023 magnitude 7.4 event. The December 2023 event likely resulted from subduction of an extinct fracture zone (Aurelio et al., 2023). We posit that the recent mainshock resulted from renewed movement along the subduction zone where the extinct fracture zone continues to be consumed at the trench.

In our earlier work, we wrote, “the population in eastern Mindanao should be reminded to be alert against any strong earthquakes in the future that may be triggered by this magnitude 7.4 event.” Indeed, the “future” event occurred eight months later (Aurelio et al., 2023).

At this point, it is important to reiterate that similarly complicated tectonic settings around the Philippines are capable of generating similarly large magnitude earthquakes. Special attention should be given to areas with large populations, such as in metro Manila and environs (approximate population of 25 million); the offshore Manila Trench is another location where an extinct spreading ridge is currently being subducted (Pautot et al., 1986; Aurelio, 2000a).
 

Science editor: Dr. Alka Tripathy-Lang, Ph.D.
Reviewer: Hector Gonzalez-Huizar, Ph.D.
 

References

Aurelio, M.A. (1992). Tectonique du segment central de la faille Philippine: etude structurale, cinématique et evolution géodynamique. Thèse de doctorat de l’Université de Paris 6, Université Pierre et Marie Curie, Académie de Paris, T 92-22, France, 500 p.

Aurelio, M.A. (2000a). Tectonics of the Philippines revisited. J. Geol. Soc. Phil., 55, 1/2, 119-185.

Aurelio, M.A. (2000b). Shear partitioning in an island arc setting: constraints from
Philippine Fault and recent GPS Data. Island-Arc, Thematic Issue, Tokyo, 9, 585-598.

Aurelio, M.A. and Peña, R.E., Editors (2010). Geology of the Philippines, 2nd Edition – Vol. 1: Tectonics and Stratigraphy. 532 pages. Published by the Mines and Geosciences Bureau, Department of Environment and Natural Resources. Quezon City, Philippines.

Aurelio, M.A. and Almeda R.L. (1999). Active deformation and stress state in and around the Philippines: present-day crustal motion from GEODYSSEA. In: Prog. and Abs. GPS 99 – The International Symposium on GPS. Tsukuba, Japan, 18-22 October 1999.

Aurelio, M., Catugas, S.D., Anicete, A. (2023), Eastern Mindanao earthquake sequence strikes on subducting, extinct fracture zone, Temblor,
http://doi.org/10.32858/temblor.331

Cardwell, R.K., Isacks, B.L., and Karig, D.E. (1980). The spatial distribution of earthquakes, focal mechanism solutions and subducted lithosphere in the Philippine and northeastern Indonesian Islands. In: Hayes, D.E. (ed.) The Tectonic and Geologic Evolution of Southeast Asian Seas and Islands, Part 1. Am. Geophys. Union Monograph, 23, 1-35.

Jarvis, A., H.I. Reuter, A. Nelson, E. Guevara (2008). Hole-filled SRTM for the globe Version 4, available from the CGIAR-CSI SRTM 90m Database (http://srtm.csi.cgiar.org).

Pautot, G., Rangin, C., Briais, A., Tapponier, P., Beuzart, P., Lericolais, G., Mathieus, X., Wu, J., Han, S., Li, H., Lu, Y. and Zhao, J. (1986). Spreading direction in the Central South China Sea. Nature, 321, 6066, 150-154.

Philippine Institute of Volcanology and Seismology Earthquake Data (2024). Available at: https://swift1.phivolcs.dost.gov.ph/swiftdb/

Quebral, Ramon D. (1994). Tectonique du segment meridional de la faille philippine, Mindanao Oriental, Philippines: passage d’une zone de collision à une zone de décrochement. These de doctorat, Université Pierre et Marie Curie, Paris, France.

Rangin, C., Le Pichon, X., Mazzotti, S., Pubellier, M., Chamot-Rooke, N., Aurelio, M., Walpersdorf, A. and Quebral, C. (1999). Plate convergences measured by GPS across the Sundaland/Philippine Sea Plate deformed boundary: the Philippines and eastern Indonesia. Geophys. J. Int., 139, 296-316.

Toda, S., Stein, R.S., Sevilgen, V. and Lin, J. (2011). Coulomb 3.3 Graphic-rich deformation and stress-change software for earthquake, tectonic, and volcano research and teaching—user guide: U.S. Geological Survey Open-File Report 2011–1060, 63 p., available at https://pubs.usgs.gov/of/2011/1060/

Weatherall, P., Tozer, B., Arndt, J.E., Bazhenova, E., Bringensparr, C., Castro, C.F., Dorschel, B., Ferrini, V., Hehemann, L., Jakobsson, M., Johnson, P., Ketter, T., Mackay, K., Martin, T.V., Mayer, L.A., McMichael-Phillips, J., Mohammad, R., Nitsche, F.O., Sandwell, D.T., Snaith, H., Viquerat, S. (2020). The GEBCO_2020 Grid – a continuous terrain model of the global oceans and land. British Oceanographic Data Centre, National Oceanography Centre, NERC, UK. doi:10.5285/a29c5465-b138-234d-e053-6c86abc040b9

Wessel, P. and Smith, W.H.F., (1995). New version of the Generic Mapping Tools released. EOS Trans. Am. Geophys. Union 76, 329.
 

Copyright

Text © 2024 Temblor. CC BY-NC-ND 4.0

We publish our work — articles and maps made by Temblor — under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license.

For more information, please see our Republishing Guidelines or reach out to news@temblor.net with any questions.