the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 22 Jun 2026
The 2022 Rainfall-Triggered Landslide-Tsunami Disaster in Pilar, Abuyog, Leyte, Philippines and its implications for coastal hazard assessment
Abstract. On 12 April 2022, a rainfall-triggered landslide in Barangay Pilar, Abuyog, Leyte, Philippines generated destructive tsunamis upon reaching the coast, resulting in one of the deadliest landslide disasters in the country’s history. The event occurred during prolonged and intense rainfall associated with the interaction of a low pressure area, Severe Tropical Storm Megi (Agaton), and Typhoon Malakas (Basyang). The cascading rainfall–induced landslide–tsunami hazards caused 54 fatalities, injured 49 individuals, and left 33 persons missing and presumed dead, with tsunami inundation responsible for a substantial proportion of the losses. We reconstructed the event chronology and failure mechanisms using an integrated dataset comprising satellite imagery, unmanned aerial vehicle (UAV) surveys, field mapping, sedimentological analyses, and three-dimensional terrain modelling, eyewitness accounts, social media videos, and incident reports. Results indicate that two landslides with a total volume of around 5.5 x 105 m3 were mobilized as debris flows. The landslides were composed of volcanic materials that originated from a fault-controlled ridge south of Pilar that reached the coast of Leyte Gulf. The initial and larger debris flow generated tsunami waves with runup heights of about 12 m above sea level, devastating residential areas along the coast and overwhelming people fleeing the landslide. Analysis of landslide–tsunami coupling indicates that wave generation was governed primarily by high-momentum impact conditions, with debris flow velocities of up to 25 m/s and short source-to-shore distances limiting energy dissipation prior to coastal entry. These conditions produced strongly supercritical flow and demonstrate that relatively moderate-volume landslides can generate significant near-field tsunamis when velocity and proximity are favorable. This study is the first well-documented subaerial landslide-generated tsunami in the Philippines and provides a data-rich example of a rainfall-induced landslide–tsunami cascade in a tropical island setting. The Pilar disaster highlights the extreme hazard posed by rainfall-triggered landslides in steep coastal volcanic terrains and underscores the need to explicitly incorporate landslide–tsunami cascade scenarios into hazard assessment, early warning systems, and land-use planning in similar coastal regions.
Competing interests: Richard Ybanez is a PhD student of NHESS Executive Editor Dr. Bruce D. Malamud in the University of the Philippines National Institute of Geological Sciences. Dr. Alfredo Mahar Francisco A. Lagmay is the main supervisor of Mr. Ybañez in the same research. No other competing interest is known to the authors.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.- Preprint
(7469 KB) - Metadata XML
-
Supplement
(34051 KB) - BibTeX
- EndNote
Status: open (until 03 Aug 2026)
-
RC1: 'Comment on egusphere-2026-2543', David Tappin, 17 Jul 2026
reply
-
AC1: 'Reply on RC1', Likha G. Minimo, 18 Jul 2026
reply
We summarized our response to RC1 in this table.
RC 1 AC It’s an interesting paper with a great deal of important information, but in its present form it’s hard to see the wood for the trees.
Thank you for your comments! We’ll work on the details to make the presentation clearer. Paper strengths are the mapping of the landslide area and the images presented and the eyewitness interviews. This is a great data set and the basis of the paper. However, the interpretations are poorly shown in the key figures (e.g. 3), which are too small and with too low resolution to see the features described in the text – e.g. Lines 200-201 “…the Pilar landslide originated from a steep, fault-cut volcanic ridge 590m SSW of Barangay Pilar (Fig. 3A)”. The ridge is not identified on Fig 3A, and the resolution is too low to distinguish the features. Lines 202-203 “These features were more defined in the post-landslide DSM (Fig. 3B). Line 204 “Field mapping identified prominent landslide scarps and arcuate failure surfaces upslope of the deposit fan (Fig. 3B)”. I disagree; these features cannot be discerned on Fig 3B.
We can enlarge Fig.3A to maximize the width of the page (with margins) for the drone orthophoto, and provide an inset pre-landslide DEM that highlights the geomorphological features of the study area. We’ll label these features in the map. There are three debris flows, DF1a, DF1b and DF2, but there no figures showing these, so I found it hard to understand the distribution of the different debris flows – they need be shown on a map overview? DF1a and DF1b are described in the text but not distinguished on the figures – just DF1. The models in Figure 4 are much too small to distinguish between them. DF1a and DF1b were distinguished by some eyewitnesses; others described DF1 as one flow event that continued as the T1 was triggered as DF1 initially hit the water. In the field, we could not distinguish between DF1a and DF1b; these deposits are also mostly buried under DF2 (supported by Fig.7). We labelled both DF1 and DF2 in Fig.4, but we’ll add these extents in the Fig. 3 revision. There are no figures showing the results of the geological models from Leapfrog Geo (Lines 172-173); especially the changes in heights down the valley from erosion at the top to deposition at the bottom initiation, entrainment/transport and deposition areas need to be shown on a map and in figures. There is no validation presented for the modelling – Lines 177-178 – “Model outputs were cross-checked against field observations and mapped deposit extents to ensure consistency”.
Leapfrog Geo models are shown in Fig. 3B and Fig.5, we’ll add that reference to the figures in the text to better guide the readers. The elevation differences are shown in the Supplementary File 2 as mentioned in line 170-171. We’ll add the thickness measurements from the fieldwork that were used to create the 3D deposit model (currently Fig. 3B). We'll also add the photos of the outcrops from which these were measured as a Supplementary file. The different landslides are composed of different sediments DF1a and DF1b – grey and DF2 brown – why, how do these relate to their source volcanic sediments and the location of these? We clarified these in the description of the initiation zone (5.4.1, lines 269-271). Both the eastern and western side of the zone has the sequence of “light gray source rock that becomes reddish (oxidized) and weathered into brown soil towards higher elevation. Near the top of the soil mantle is a distinct boulder-rich layer upon which vegetation was anchored, as seen both on the slumping units and the scarp.”
We illustrated 3 different models that could fit the varying eyewitness descriptions and the composition of the deposits (lines 274-287). In all these models, DF1 could have more materials from the gray source rock unit (Fig.4A), while the DF2 deposit and the slumping unit could have more from the overlying sediments and the reddish soil mantle.
Regarding the three alternative models how do these relate to the fault cut volcanic ridge - - Line 201 – “…..the Pilar landslide originated from a steep, fault-cut volcanic ridge 590m SSW of Barangay Pilar (Fig. 3A)”. This is not identified on any figure (as I note above) so it is hard to understand the interpretation of the different three models, (which are difficult to discern on Figure 5). In the interpretation of the three models, A failed on the western crown and B and C on the eastern crown – how do these different models relate to the fault cut volcanic ridge? All these models have the source rock fractures near the western edge of the failed volume of DF1 or DF1a. Both DF1 and DF2 passed through the fault identified at S6 (shown in Supplemental File 3). We’ll improve the resolution of Fig. 3 to present this better. Different terms used – the landslides are complicated but the presentation in the paper is not helped by the different terms used. On Figure 3 there are DF1 and DF2, then on Figure 4 there is added the ‘entrained volume which looks like DF2 and a slumping unit? Entrained and slumping units were only indicated in Figure 5. Figures 3 and 4 show post-landslide imagery and interpretation of the deposits in the study area, while Figure 5 was based on the pre-landslide topography, delineating the possible source areas of the slope failures, according to the eyewitness accounts and field observations. There are sections on the rainfall, local geological and geomorphological context, then the landslide event chronology, then the landslide morphology – a bit back to front? Key to the interpretations is video evidence, but this is not shown. The sequence of results presentation was according to the sequence of each method applied:
• Documentary/reports review led to the rainfall data interpretation
• Remote Sensing led to the local geological and geomorphological context analysis
• Interviews and video review led to Landslide-Tsunami event chronology
• Fieldwork, laboratory analysis and 3D modelling led to Landslide and tsunami descriptionsA link to the publicly available video of the second debris flow event, as uploaded online by the owner (Bureau of Fire Protection, 2022) was provided in the reference (https://www.facebook.com/watch/?v=2119315464893710). The rest of the online copies were only posted by secondary sources (vloggers who got the videos from eyewitnesses). We opted to show the annotated frames because the videos were taken by witnesses moving away from the landslide tsunami.
Citation: https://doi.org/10.5194/egusphere-2026-2543-AC1
-
AC1: 'Reply on RC1', Likha G. Minimo, 18 Jul 2026
reply
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 22 | 39 | 1 | 62 | 6 | 0 | 1 |
- HTML: 22
- PDF: 39
- XML: 1
- Total: 62
- Supplement: 6
- BibTeX: 0
- EndNote: 1
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Overview
The paper addresses the cause of the April 12th, 2022, rainfall-triggered landslide and tsunami disaster in Pilar, Leyte, Philippines, addressing the event chronology, scientific insights, and recommendations for hazard assessment and disaster risk reduction. It is based on detailed fieldwork, eye-witness interviews and video footage, data analysis, and comparison with global events.
The landslide was generated by prolonged, intense rainfall (from a low-pressure area, Severe Tropical Storm Megi, and Typhoon Malakas) and there were 54 fatalities, 49 injuries, and 33 missing persons presumed dead, most from the tsunamis. The event was reconstructed from satellite imagery, UAV surveys, field mapping, sedimentological analyses, 3D terrain modelling, eyewitness accounts, social media videos, and incident reports. The fieldwork included interviews with survivors and responders.
There were two major landslides (total volume ~550,000 m³) composed of volcanic materials originating from a fault-controlled ridge and reaching the Leyte Gulf coast as debris flows. The first and larger debris flow generated tsunami waves with runup heights up to 12 meters above sea level, devastating coastal areas and overwhelming evacuees.
The tsunamis caused great destruction, with houses and infrastructure swept away, and the coastal landscape transformed with the shoreline advancing by up to 95 meters. The majority of Barangay Pilar was destroyed, and the community was permanently relocated. The loss of life was due to short warning times and high landslide flow velocities. The event is the first well-documented subaerial landslide-generated tsunami in the Philippines providing important information on the tsunami hazard in the country.
The disaster exposes the limitations of treating landslides and tsunamis as independent hazards; cascading and compound hazards can amplify disaster impacts. It demonstrates how rainfall-driven landslides in steep, volcanic coastal terrains can rapidly cascade into high-impact tsunamis, with devastating consequences for communities. The event provides lessons for improving hazard assessment, early warning, and disaster risk reduction in similar environments.
Comments
It’s an interesting paper with a great deal of important information, but in its present form it’s hard to see the wood for the trees.
Paper strengths are the mapping of the landslide area and the images presented and the eyewitness interviews. This is a great data set and the basis of the paper. However, the interpretations are poorly shown in the key figures (e.g. 3), which are too small and with too low resolution to see the features described in the text – e.g. Lines 200-201 “…the Pilar landslide originated from a steep, fault-cut volcanic ridge 590m SSW of Barangay Pilar (Fig. 3A)”. The ridge is not identified on Fig 3A, and the resolution is too low to distinguish the features. Lines 202-203 “These features were more defined in the post-landslide DSM (Fig. 3B). Line 204 “Field mapping identified prominent landslide scarps and arcuate failure surfaces upslope of the deposit fan (Fig. 3B)”. I disagree; these features cannot be discerned on Fig 3B.
There are three debris flows, DF1a, DF1b and DF2, but there no figures showing these, so I found it hard to understand the distribution of the different debris flows – they need be shown on a map overview? DF1a and DF1b are described in the text but not distinguished on the figures – just DF1. The models in Figure 4 are much too small to distinguish between them.
There are no figures showing the results of the geological models from Leapfrog Geo (Lines 172-173); especially the changes in heights down the valley from erosion at the top to deposition at the bottom initiation, entrainment/transport and deposition areas need to be shown on a map and in figures. There is no validation presented for the modelling – Lines 177-178 – “Model outputs were cross-checked against field observations and mapped deposit extents to ensure consistency”.
The different landslides are composed of different sediments DF1a and DF1b – grey and DF2 brown – why, how do these relate to their source volcanic sediments and the location of these?
Regarding the three alternative models how do these relate to the fault cut volcanic ridge - - Line 201 – “…..the Pilar landslide originated from a steep, fault-cut volcanic ridge 590m SSW of Barangay Pilar (Fig. 3A)”. This is not identified on any figure (as I note above) so it is hard to understand the interpretation of the different three models, (which are difficult to discern on Figure 5). In the interpretation of the three models, A failed on the western crown and B and C on the eastern crown – how do these different models relate to the fault cut volcanic ridge?
Different terms used – the landslides are complicated but the presentation in the paper is not helped by the different terms used. On Figure 3 there are DF1 and DF2, then on Figure 4 there is added the ‘entrained volume which looks like DF2 and a slumping unit?
There are sections on the rainfall, local geological and geomorphological context, then the landslide event chronology, then the landslide morphology – a bit back to front? Key to the interpretations is video evidence, but this is not shown.
Conclusions, Recommendations
Great data, great figures, but presentation is poor and results presented confusing. Three landslides identified but no figures showing their location and relationships. There seems to be a lot of video evidence, that would be good to see. Three alternative models are presented, but do not seem to relate to the location of their origin on the fault-cut volcanic ridge? No figures to support the relationship(s) between the modelling and the field observations and mapped deposit extents. The order of presentation - chronology before landslide description is back to front. Geographical relationships between the landslides not shown, so impossible to understand. Can the different sedimentary composition of the landslides, DF1a and DF1b – grey and DF2 brown – be related to different volcanic source sediments and the location of these?
Fundamental to the paper is a description and interpretation of the landslides, and this is missing. More figures on the landslides, their distribution and thickness, and their relationships are required and a fuller presentation of the results of geological models together with justification of these results from the landslide field mapping. The overall results from the research need improved explanation and justification. Figure 3 is critical in introducing the context of the interpretations but the images here are poorly annotated and too low resolution to identify their description in the text – so higher resolution images and descriptions required here. Is it possible to include some of the video evidence?
Overall, this promises to be a significant paper, but major revision is required to answer my comments and recommendations before publication.