Frontal wedge variations and controls of submarine landslides in the Negros–Sulu Trench System, Philippines

Frontal wedge characteristics provide clues to the efficiency of the overriding slab for large displacement during megathrust and upper-plate earthquakes, whereas submarine landslides along active margins may trigger or amplify tsunamis. The lack of clear precursors of submarine failures poses difficulty in monitoring and providing real-time alert warning systems. With that, delineating submarine features along active margins, their spatial distribution, and controls provide valuable information in identifying regions susceptible to large submarine landslides and tsunami hazard assessments. In this study, we performed terrain and morphometric analyses on 20 m resolution bathymetry data to map submarine landslides, submarine canyons, and lineaments in the forearc margin of the Negros–Sulu Trench System in the Philippines. Lineaments are distributed mainly along the frontal wedge, where previous seismic surveys revealed that the mapped ridges are morphotectonic expressions of thrusted sediments. The morphological variations of the four frontal wedge segments were attributed to heterogeneous sediment influx, convergence rates, and subduction processes. More than 1,200 submarine landslides and their morphometric parameters were delineated, and exploratory spatial analyses indicate clustering and underlying controls. The tendencies of prolate submarine landslides (high L/W) to significantly cluster along submarine canyons while oblate morphologies (low L/W) along the frontal wedge reflect the different environments and geomorphological conditions to form these contrasting shapes. Ubiquitous small submarine landslides are mainly controlled by submarine canyon systems at relatively shallow depths of <2 km, where high sediment influx from inland sources preconditions instability. Large submarine landslides (>0.5 km3), on the other hand, are significantly most clustered where the Cagayan Ridge seamount collides and subsequently subducts beneath the northernmost frontal wedge. This suggests the dominant role of seamount subduction and related tectonic processes causing slope steepening to mainly induce large submarine landslides. This study unveiled how submarine landslides vary morphologically depending on their spatial, geomorphological, and tectonic controls in the active margin. This new information provides clues in identifying offshore areas susceptible to large submarine landslides that may induce damaging tsunamis in the Negros–Sulu Trench System as well as in other active margins of similar underlying controls.


Introduction
Submarine landslides occur when driving stresses on sediment layers increase or when shearing resistance decreases and are mainly triggered and preconditioned by the following: earthquakes, slope steepening, volcanic eruptions, rapid loading causing overpressure, wave and tidal loading, rapid sediment loading, gas hydrates, fluid migration, and weak stratigraphic layers (Vanneste et al., 2013;Scarselli 2020;Watson et al., 2020). Moreover, submarine landslides can induce tsunamis. This occurs when seafloor vertical displacement is translated into the seawater column and is a function of submarine landslide volume, depth, acceleration, velocity, and fluid density (Harbitz et al., 2014; Yavari-Ramshe and Ataie-Ashtiani 2016; Løvholt et al., 2017;Baba et al., 2019). Complex triggering factors with no clear precursors are the main challenges in creating an effective real-time early-warning system for submarine landslide-driven tsunamis.
In the past, submarine landslides were overlooked to cause damaging tsunamis until the 1998 Papua New Guinea submarine slump. This event caused a 15 m wave height leading to more than 2,200 fatalities that reshaped the perspective on the impacts of submarine landslides (Heinrich et al., 2001;Tappin 2021). Advancements in seafloor mapping and tsunami modeling have also led to revisiting past tsunamis that were once attributed only to earthquakes. Tsunami modeling of the following events was found to either be amplified or mainly induced by an earthquake-triggered submarine landslide: 1908 Mw 7.1 Messina (12 m maximum tsunami, 80,000 fatalities) (Schambach et al., 2020), 1945 (Rastgoftar andSoltanpour 2016), 1946 Ms 7.1 Aleutian (16 m, 167) (Fryer et al., 2004), 1964 Mw 9.2 Alaska (12 m, 131) (Parsons et al., 2014), 1992 Ms 7.5 Flores Island (26 m, 1,713) (Hidayat et al., 1995;Imamura et al., 1995), 1994 Mw 7.1 Mindoro (7 m, 78) (Ramirez et al., 2022), as well as the more recent 2011 Mw 9.0 Tohoku (40 m, 18,490) (Tappin et al., 2014) and 2018 Mw 7.5 Sulawesi (10 m, 2,000) (Takagi et al., 2019;Nakata et al., 2020) tsunami events. The 1945Makran, 1946Aleutian, 1964Alaska, 1992Flores Island, and 2011 Tohoku tsunami events occurred along active margins, where subduction processes predominantly trigger submarine failures. Recent studies linked the distribution of submarine landslides in forearc basins to seismic loading, oversteepening of the frontal wedge, seamount subduction, and submarine canyon development (Mountjoy et al., 2009;Pedley et al., 2010;Watson et al., 2020). The frontal wedge is of particular interest as its configuration implicates the variability of subduction processes along the trench. This includes whether a trench segment undergoes accretion or erosion (e.g., Kao et al., 2000;Pedley et al., 2010;Armada et al., 2020) as well as the relationship between the frontal wedge morphometry and maximum tsunami wave heights for megathrust earthquakes (Qiu and Barbot 2022). The link between the frontal wedge and subsequent submarine landslides was also documented in other active margins including the Nankai Trough (Moore et al., 2019), Hikurangi Trough (Pedley et al., 2010), and Java Trench (Kopp et al., 2006). These studies highlight the role of variability in subduction processes along the trench and the consequent submarine slope failures.
Mapping both the frontal wedge and other submarine features as well as conducting exploratory spatial statistics can be significant initial steps toward providing informed decisions in minimizing the impacts and inferring seafloor areas that are susceptible to tsunamigenic submarine landslides (e.g., Völker, 2010;Kawamura et al., 2012;Parsons et al., 2014;Watson et al., 2020;Gamboa et al., 2022).
The Philippines is an archipelago surrounded by subduction zones with the potential to generate strong to great earthquakes that can trigger submarine landslides. Despite its active seismicity and tectonics, the Negros-Sulu Trench System (NSTS) remains understudied, and its submarine geomorphological features are poorly constrained. In this study, we mapped the frontal wedge variations, submarine canyons, lineaments, and submarine landslides in the active margin of the Sulu Sea using highresolution multibeam bathymetry data. Furthermore, this study aims to infer the controls of submarine landslides through their morphometric parameters and exploratory spatial analyses. For the first time, the morphological characteristics of these submarine features in the NSTS are examined in detail, with a particular focus on the morphological variations and spatial distribution of submarine landslides.

Sulu Sea tectonic setting and historical tsunamis
The NSTS is an active margin where the Sulu Sea subducts beneath the Philippine Mobile Belt (PMB). This marginal basin is enclosed by the Zamboanga Peninsula and the islands of Sulu, Panay, Negros, and Palawan. The Sulu Sea is divided into four major tectonic terranes: northwest Sulu basin (NWS), Cagayan Ridge (CR), southeast Sulu basin (SES), and the accretionary prism (AP) along the NSTS (Rangin, 1989). The Palawan Microcontinental Block (PCB) terrane occupies north of the Sulu Sea. The NWS has thicker crust (>10 km) than the SES (6 km), and between these two sediment-filled sub-basins is the Cagayan Ridge seamount composed of basalt, tuff, and andesite (Rangin and Silver 1990) (Figure 1). GPS surveys revealed that this trench system generally has low convergence rates (<60 mm yr −1 ) (Rangin et al., 1999;Simons et al., 1999). Two prominent models suggest that the Sulu Sea was formed due to back-arc rifting from the subduction of either the proto-South China Sea along the Cagayan Arc (e.g., Bellon & Rangin, 1991;Rangin & Silver, 1991;Schlüter et al., 2001;Liu et al., 2014) or the Celebes Sea along the Sulu and Zamboanga Peninsula Arc (Rangin, 1989;Hall 2002Hall , 2013Hall , 2012Lai et al., 2021).
Since 1589, three historical tsunamigenic earthquakes occurred along the NSTS: 1897 Ms 7.5, 1925 Ms 6.8, and 1948 Ms 8.2 events (Bautista et al., 2012). The 1897 earthquake event had an epicenter near the Zamboanga Peninsula and generated up to a 6 m tsunami that widely affected the coastal areas in the Sulu Sea. The basin-wide distribution and maximum tsunami wave height for a Ms 7.5 earthquake pose questions regarding its mechanism. Whether this event was a tsunami earthquake or triggered by a submarine landslide is still unknown. On the other hand, the 1925 and 1948 earthquake events generated a tsunami of up to 2 m that only affected southeastern Negros and Panay islands, respectively. The 1948 event is distinct among the three as its epicenter is inland. The rupture parameters of these tsunamigenic earthquakes, however, are still poorly constrained.

Delineating submarine landslides and their morphological parameters, submarine canyons, and lineaments
The bathymetry data of the National Mapping and Resource Information Authority (NAMRIA) has a 20 m resolution and a depth accuracy of 1 m. These were collated from the multibeam bathymetric surveys in the Sulu Sea between 1999 and 2018 by the NAMRIA and the Naval Oceanographic Office (NAVOCEANO, United States). Terrain analyses were applied to the bathymetry data including slope, roughness, and hillshade. Slope was derived using the algorithm of Horn (1981). By showing both the slope and depth through transparency adjustments, submarine landslides, canyons, and lineaments were mapped. To delineate submarine canyons, drainage extraction algorithms in ArcGIS were utilized including flow direction and flow accumulation.
The submarine landslides were mapped based on the following key morphological features: head scarp (scar) characterized by a relatively steep slope at the top of the landslide, arcuate or concave morphology of the landslide area, and the presence of mass transport deposits (MTDs) at the toe of the landslide (Watson et al., 2020). Majority of the mapped landslides, however, are based mainly on head scarp and arcuate morphologies due to the rare occurrence of MTDs observed. The apparent rare occurrence of MTDs is attributed to the resolution of the bathymetry and the low preservation of these deposits once disintegrated. Interactive 2D cross-sectional profiles in QGIS also aided in recognizing the concave morphology of submarine landslides, especially at smaller scales. The submarine landslides are further categorized into four levels (1-4) of confidence scale (Watson et al., 2020) based on the presence of key morphological features of submarine landslides (Figure 2).  Rangin, 1989). Arrows indicate the velocity vectors (mm yr −1 ) from continuous GPS campaigns (1994)(1995)(1996) across the Philippines (Rangin et al., 1999;Simons et al., 1999). Red circles indicate the epicenter of the three historical tsunamigenic earthquakes in the Sulu Sea (Bautista et al., 2012). The onshore faults, lineaments (black lines), and catalog (1589-2020) of surface-wave magnitude (Ms) earthquakes (EQ) are from PHIVOLCS. The inset map shows the location of the Sulu Sea (SS) and the other surrounding basins (West Philippine Sea, WPS; Cotabato Sea, CS; Philippine Sea, PS) and trenches (Manila Trench, MT; Negros Trench, NT; Sulu Trench, ST; Cotabato Trench, CT; Philippine Trench, PT; East Luzon Trough, ELT).

Frontiers in Earth Science
frontiersin.org Table 1 shows the delineated morphometric parameters of submarine landslides modified after Gamboa et al. (2022) and Watson et al. (2020). Volume, in particular, was estimated by first creating a pre-landslide digital elevation model (DEM) similar to the methods of Gamboa et al. (2022). The prelandslide DEM is estimated by sampling elevations around the perimeter of the landslide and applying a multi-level b-spline interpolation at a larger 100 x 100 m raster grid. Volume is calculated based on the surface difference between the prelandslide DEM as the top surface and the present bathymetry data as the base. Figure 3 shows a sample of the modeled prelandslide surface and the present scarp morphology for volume calculation.

Statistical analyses of submarine features
Two-dimensional (2D) histograms were applied among the morphometric parameters to infer correlations. Ripley's L function (Ripley 1977) was calculated through the spatstat package in R to determine whether the distribution of submarine landslides has significant clustering across different distance scales.
The L function is a linear transformation of the K function expressed in distance and calculated by: where A = area of the bathymetry coverage, n = total number of mapped submarine landslides, r = distance, and k(i, j) is the weighting factor for edge detection. A neighboring (j) submarine landslide is included in the summation when its distance from the reference submarine landslide (i) is ≤ r. The centroid of each submarine landslide served as an input point. Observed L(r) greater than the expected values indicate significant clustering, whereas lower observed L(r) suggest dispersion. A total of 999 permutations were run to calculate the confidence interval (CI) envelope of the L(r) signifying complete spatial randomness. Kernel density analysis was performed in ArcGIS to estimate the 2D probability density distribution (Silverman 1986) of the delineated submarine landslides, submarine canyons, and lineaments. The distribution of earthquakes (1589-2020) from the Philippine Institute of Volcanology and Seismology (PHIVOLCS) was also included in the kernel estimation. Kernel Frontiers in Earth Science frontiersin.org density can provide insights into the intensity variation, spatial correlation among submarine features, and underlying processes that may be linked to the distribution of submarine landslides. These kernel densities are further evaluated through correlation matrices of the raster bands between submarine landslides and submarine canyons, lineaments, and earthquakes.
To examine statistically significant clusters of relatively largeand small-volume submarine landslides, we identified their spatial distribution using the Getis-Ord (Gi*) statistic (Getis and Ord 1992) that calculates the z-score and equivalent p-values of each submarine landslide: where G * i = z-score, x j = volume, w i,j = weighting factor based on the spatial relationship of i and j neighbors at distance d, x = mean sample volume, and s = volume standard deviation among neighbors. The optimized hot spot analysis in ArcGIS based on Getis-Ord assigns a fixed distance d based on the first peak clustering (=14.3 km). Submarine landslides were then grouped into seven: one that belongs to null groups (1), three that belong to significant clusters of high values at 99% confidence interval (CI) (2), 95% CI (3), 90% CI (4), and three that belong to significant clusters of low values at 99% CI (5), 95% CI (6), and 90% CI (7). These groups are then compared based on select morphometric parameters using boxplots. In addition to volume, another morphometric parameter that was calculated with Gi* is the length-to-width ratio (L/W) to explore significant clusters of prolate (high L/W) and oblate (low L/W) submarine landslides relative to the flow direction. Altogether, these exploratory spatial analyses aim to Distance from submarine canyon The shortest distance of the submarine landslide from a submarine canyon (km)

FIGURE 3
Three-dimensional (3D) perspective of the modeled pre-landslide and the present surface of the submarine landslides. (A) The modeled prelandslide is calculated based on the multi-level b-spline calculation at a coarser 100 x 100 m resolution around the perimeter of the submarine landslide that represents the intact slope. (B) The present surface DEM serves as the base, whereas the modeled pre-landslide is the top surface for estimating the slump volume.
Frontiers in Earth Science frontiersin.org unravel spatial patterns, providing insights into the occurrence and controls of submarine landslides in the NSTS.

Submarine landslides, submarine canyons, lineaments, and frontal wedge variations
A total of 1,214 submarine landslides were mapped ( Figure 4) with 64.7% under category 4, 31.6% category 3, 3.1% category 2, and 0.6% category 1. The multitude of categories 3 and 4 reflects the limitation of the bathymetry resolution in readily discerning head scarp and arcuate morphologies but without traces of MTDs. Large networks of submarine canyons were mapped between Negros Island and Zamboanga Peninsula and between Panay and Negros islands. Submarine gullies develop at the shelf edge of about 180 m below sea level and interconnect with larger submarine canyons at deeper portions as slope abruptly increases (Figures 5D, E). Lineaments are mainly situated along the frontal wedge, occurring as parallel ridges, with the majority following a north-to-northeast trend. In addition, the northwest-trending   Figure 5D.
Interpreted seismic reflection profiles ( Figure 6) by Schlüter et al. (1996) across the trenches revealed that these mapped frontal ridges are seafloor manifestations of accreted trench sediments thrusted by  Figure 5B. The identified oceanic splinter is linked to intense thrusting and folding along this segment.

Frontiers in Earth Science
frontiersin.org 08 reverse faults. Profile A-A′ shows the collision and subsequent subduction of the Cagayan Ridge beneath the northern Negros Trench forearc. In profiles B-B′ and C-C′, the décollement and subduction of the southeast Sulu Sea along the frontal wedge are shown. Notably, profile C-C′ also depicts a crustal splinter that is associated with intense thrusting and folding along this segment.
Four segments were delineated based on frontal wedge variations: (1) The northern Negros Trench (NT1) is northeast trending (15°azimuth) and transects between 10.86°N and 9.76°N with a length of about 126 km. This segment is bounded by the Cagayan Ridge seamount chain and southern Panay Island. Its frontal wedge has a maximum width of about 30 km and steepens subvertically northward.

Morphometry of submarine landslides
Morphometric parameters that are size-dependent, i.e., area, volume, width, height, length, and head scarp length, all follow a power-law distribution with ubiquitous low values ( Figures 7A-F). The distance from the submarine canyon also follows the same distribution. On the other hand, L/W, maximum slope, and mean, top, and maximum depth have skewed binomial distributions that peak around 1.4-1.8 km. Of the 1,214 submarine landslides, 78% have an area and volume of <5 km 2 and 0.5 km 3 , respectively (Figure 7). The mean area and volume are 1.8 km 2 (median = 4.0 km 2 ) and 0.03 km 3 (median = 0.18 km 3 ), respectively. The average maximum slope is 38°, and the majority (23%) of distances from the trenches were <20 km. On the other hand, 77% of submarine landslides are less than 5 km away from submarine canyons. Area, volume, length, width, distance from submarine canyons, head scarp length, and height were scaled logarithmically in the 2D histogram to normalize highly skewed power-law distributions. A base of 10 was arbitrarily chosen for the logarithmic normalization. Figure 8 and Figure 9 show the 2D histogram correlation relative to the volume and L/W, respectively, among morphometric parameters. Area, length, width, headscarp length, and height show positive correlations with volume as these are typically size-dependent. The maximum slope is also positively correlated with volume. Although the plots were highly scattered, an apparent inverse trend is observed between volume and the distance from the trench and L/W. Top and mean depths show an apparent positive trend where volume increases with depth up to around 3 km and then shifts to an inverse pattern toward 5 km depth. For L/W correlations, the rest of the parameters have highly scattered plots, except for the apparent negative trend with respect to area, width, maximum slope, and distance from submarine canyons. In addition, an apparent positive trend is observed between length and L/W.

Spatial distribution and correlation of submarine features
Ripley's L function (Figure 10) shows higher observed L (L obs(r) ) than the expected (L theo(r) ) values, indicating a non-random distribution and clustering at a wide range of distances. This is Other parameters (E-L) have a weak association with volume, but an apparent inverse trend can be observed with (F) distance from the trench and (I) L/W. Volume is observed to increase with (K, L) depth up to around 3 km and then shifts to a negative trend toward 5 km depth.

Frontiers in Earth Science
frontiersin.org exemplified by the inhomogeneous kernel density of submarine landslides in Figure 11A. The highest intensity distribution of submarine landslides is located offshore between southeastern Negros Island and northern Zamboanga Peninsula, as well as between southern Panay and southwestern Negros islands. Submarine canyons also have a similar distribution to submarine landslides ( Figure 11C), whereas lineaments ( Figure 11B) are concentrated in the frontal wedge of the trench. Matrix correlation among kernel densities including earthquakes reveals that submarine canyons have the highest correlation (0.76) with the distribution of submarine landslides. Regions with significant clusters of high and low volume and L/W based on Getis-Ord statistics are shown in Figure 12. Significantly large-volume (>0.5 km 3 ) submarine landslides are clustered offshore of southern Panay Island as well as offshore southeastern Negros, Zamboanga Peninsula, and Sulu islands. On the other hand, those with a significantly small volume (<0.05 km 3 ) are concentrated where the highest probability density of submarine landslides and canyons occurs ( Figures 11A,C, respectively). Significantly high L/W submarine landslides ( Figure 12B) are clustered along areas with high kernel density of submarine canyons ( Figure 11C). Those with significantly low L/W are clustered along the frontal wedge and abyssal plains.

Significant clusters of high and low volume and L/W
Significant clusters of submarine landslides were grouped based on the degree of statistical significance (% CI), while their depth and

Frontiers in Earth Science
frontiersin.org slope parameters were plotted in boxplots ( Figure 13). Significant clusters of large-volume submarine landslides (90%, 95%, and 99% CI) have a wide interquartile range (IQR) of top and mean depths ( Figures 13A, B), whereas significantly small-volume (90% and 95% CI) submarine landslides have a narrow IQR around 1-2 km depths. An opposite trend is observed for those with significantly high and low L/W ratios ( Figures 13E, F). Submarine landslides with high L/W ratios (90%, 95%, 99% CI) have high distribution occurring at shallower depths, whereas those with low L/W ratios (90%, 95%, 99% CI) occur at deeper levels with a wider IQR. In terms of slope, those with significantly high volume have generally steeper mean and maximum slope than those with significantly low volume ( Figures 13C, D). On the other hand, groups in terms of L/W have a similar range of mean slope ( Figure 13G), but the maximum slope decreases with higher L/W ratios ( Figure 13H).

Frontal wedge variations along the NSTS
The well-developed frontal wedge of the three segments (NT1, NT2, ST1) of the NSTS is linked to relatively high sedimentation rates of this active margin demonstrated by dense networks of submarine canyons (i.e., conduits of inland clastic deposits) and low overall convergence rates (54 mm yr −1 , NT; 28 mm yr −1 , ST) based on previous GPS surveys (Rangin et al., 1999;Simons et al., 1999). High sediment influx and slow convergence rates (<60 mm yr −1 ) are conducive conditions in developing accretionary margins and frontal wedge (e.g., Huene and Scholl, 1991;Lallemand et al., 1994;Clift and Vannucchi 2004;Simpson 2010). In addition, morphological variations of the frontal wedge along the NSTS are attributed to the heterogeneity of sediment distribution from nearby islands and variability of convergence rates and forearc deformation. NT1, NT2, and ST1 segments are relatively closer to Panay, Guimaras, Negros, and Zamboanga Peninsula. These segments are presumed to receive higher sediment influx due to larger land areas than the ST2 segment which is closer to the smaller Sulu Group of Islands (Figure 1). Similar findings along the northern segments of both the Manila (Armada et al., 2020) and Japan (Tsuru et al., 2002;Kodaira et al., 2017) trenches also attributed the well-developed accretionary prism to higher sedimentation influx than the southern segments.
Another variation observed on the trench segments is their steepness, most notably on NT1, which is marked by an oversteepened frontal wedge. While well-developed frontal wedges have been correlated to high sediment influx and slow convergence rates, oversteepening is mainly driven by seamount subduction. These bathymetric highs cause intense uplift, deformation, and thrust faulting of the overriding plate inducing steepening and erosion (Yang et al., 2022). Examples of these frontal wedge oversteepening include segments of the Hikurangi Trough (Pedley et al., 2010), Java Trench (Masson et al., 1990), Middle America Trench (Hühnerbach et al., 2005), and Nankai Trough (Bangs et al., 2006;Moore et al., 2019). Similarly, the collision and subsequent subduction of the northeastern portion of the Cagayan Ridge along NT1 presumably induce oversteepening of the frontal wedge as well as its orientation change to the northeast.
The characteristics of the four trench segments based on frontal wedge variations provide essential information not only on the potential impacts of submarine landslides but also for modeling future megathrust tsunami scenarios. Furthermore, the width of the frontal wedge is indicative of its tsunami hazard potential as tsunami run-ups of shallow megathrust ruptures have been correlated to the width of the frontal wedge (Qiu and Barbot 2022). This is further emphasized by historical records where great offshore earthquakes occur mostly along accretionary prisms with thick sediment cover (Bilek and Lay 2018).

Controls of submarine landslide distribution
Exploratory spatial analyses revealed non-random and significant clustering of submarine landslides in which ubiquitous and small-volume landslides are controlled predominantly by submarine canyon systems. These conduits of terrestrial sedimentation influx induce overloading and instability to precondition slope failures. In addition, significant clustering of prolate (high L/W) submarine landslides coincides with the distribution of submarine canyon systems. The resemblance of prolate morphology to submarine canyons implies that these prolate submarine landslides influence the development of submarine canyon networks. Small submarine landslides may eventually interconnect forming incisions for sediment transportation (Baztan et al., 2005;Micallef et al., 2012). This relationship between prolate submarine landslides and submarine canyons is further supported by the apparent inverse trend between L/W and the distance from submarine canyons ( Figure 9F). These

FIGURE 10
Ripley's L function of submarine landslide distribution. The observed L (L obs(r) ) values are significantly higher than the expected (L theo(r) ) across different distances indicating significant clustering and inhomogeneity over various distances. The confidence interval (CI; L hi(r) and L hi(r) ) envelope is calculated based on 999 permutations to compare L(d) at a certain distance where points are distributed randomly.

Frontiers in Earth Science
frontiersin.org clustered submarine landslides with high L/W are prevalent in relatively shallower depths around 1-3 km where submarine canyons typically form. On the other hand, significant clusters of oblate (low L/W) submarine landslides mainly occur along the frontal wedge and in abyssal plains ( Figure 12B). This is due to the subparallel and steep ridges along the frontal wedge that are suitable morphotectonic structures to form oblate submarine landslides. While small landslides are mainly controlled by submarine canyons, large ones are significantly clustered along steep slopes and at varying depths. Most notably, the NT1 frontal wedge is marked by the densest cluster of large submarine landslides. As discussed in Section 5.1, this segment has an oversteepened and wide frontal wedge linked to the collision and subsequent subduction of Cagayan Ridge seamounts. Thus, we associate the seamount subduction and subsequent slope oversteepening as dominant preconditioning factors for the occurrence of large submarine landslides along this segment. Modeling of seamount subduction beneath the accretionary prism, in its initial stage, resulted to oversteepening of the frontal wedge, reactivation of frontal thrust faults, and large submarine landslides (Dominguez et al., 2000;Ruh 2016;Morgan and Bangs 2017). This is supported by other accretionary prisms where seamount subduction underlies large submarine landslides (e.g., Hühnerbach et al., 2005;Pedley et al., 2010). While it is unclear what causes the clustering of large submarine landslides in other portions, the role of nearby fault structures (e.g., SCDL), intense tectonic deformation and fully subducted seamounts, together with sediment overloading could all play complex roles to induce steepening. Nonetheless, the identified regions with clusters of large submarine landslides

FIGURE 11
Kernel density variations (per km 2 ) of (A) submarine landslides, (B) earthquakes, (C) submarine canyons, and (D) lineaments across the NSTS. Earthquake kernel density is distributed mainly along the Negros Trench segments, whereas that of lineaments along the frontal wedge. Matrix band correlations showed that submarine canyons have the highest correlation to submarine landslide distribution of 0.76, followed by earthquakes (0.46) and then lineaments (0.20).

Frontiers in Earth Science
frontiersin.org

FIGURE 12
Getis-Ord (Gi*) analysis of submarine landslide distribution with respect to volume and L/W. (A) Significant clusters of large-volume submarine landslides (red dots) are found in different regions of the active margin, but the densest clusters are in the northern NT segment (NT1), where collision with the Cagayan Ridge occurs. Significant clusters of small-volume submarine landslides (blue dots) are located offshore of Negros Island and Zamboanga Peninsula, corresponding to areas with the densest submarine canyon systems ( Figure 11A). (B) Significant clusters of high L/W (red dots) are also in regions with dense submarine canyon systems. On the other hand, the frontal wedge and abyssal plain host significant clusters of low L/W (blue dots).

FIGURE 13
Boxplots showing the depth and slope of submarine landslides that are grouped based on the Gi* statistic (blue boxplots, based on volume; yellow green, based on L/W) in Figure 13 (red hue, hot spots; blue hue, cold spots). Red dots indicate the outliers, while broken lines inside the interquartile range (IQR) show the mean value. (A-D) Box plots based on the volume show that significant clusters of large submarine landslides have depths with wide IQR (0.1-4 km), while small ones are in a narrow IQR of about 0.1-2 km. Mean slope (C) tends to concentrate between 10°and 30°with slightly lower values for significantly small-volume submarine landslides. (D) Maximum slope showed a positive association with clusters of significantly large-volume submarine landslides, whereas small-volume landslides have a lower maximum slope. (E-H) Clusters of significantly high L/W occur in shallower depths and slightly lower slopes than those of significantly low L/W.
Frontiers in Earth Science frontiersin.org provide significant insights into modeling landslide-driven tsunamis, where the estimated volume is highly correlated with potential tsunami wave height (Sabeti and Heidarzadeh 2022). This study highlights the dominant role of submarine canyon systems due to high sediment influx and subduction-induced oversteepening to precondition slope instability. As the slope becomes unstable, gravity and earthquakes can readily trigger submarine failures along active margins (Masson et al., 2006;Scarselli 2020). The possible presence of gas hydrates may also be involved as they occur along margins with high sedimentation influx (e.g., Gee et al., 2007;Mountjoy et al., 2014;Crutchley et al., 2016). Thus, we recommend dense seismic reflection surveys to further map the underlying structures and their variations across the active margin in detail. In addition, the current resolution of bathymetry data is inadequate to map even smaller (m-scale) submarine landslides. This has been shown by the absence of significant clusters of small submarine landslides at 99% CI. As bathymetry resolution has been a common limitation in mapping submarine landslides (e.g., Gazioǧlu et al., 2005;Watson et al., 2020;Gamboa et al., 2022), higher-resolution multibeam bathymetry surveys will further improve the characterization of submarine geomorphological features.

Conclusion
Four frontal wedge segments were delineated in the NSTS, with varying widths, lengths, and associated morphotectonic features. These variations are linked to the disparity of convergence rates along the trench, subduction processes, and heterogeneous distribution of sediments from inland sources. A total of 1,214 submarine landslides and their morphometric parameters along this active margin were delineated. Prolate submarine landslides (high L/W) tend to form in submarine canyons, whereas oblate morphologies (low L/W) generally occur along the frontal wedge and abyssal plains. These opposing submarine landslide morphologies reflect their different environments and geomorphological conditions. Ubiquitous small submarine landslides are mainly controlled by submarine canyon systems where high sediment influx from inland sources preconditions instability. These small-volume submarine landslides also occur in relatively shallow depths of about 1-2 km and generally have a lower mean and maximum slope. On the other hand, seamount subduction and related tectonic processes that induce oversteepening play dominant role to precondition large submarine landslides. This study revealed the influence of spatial, geomorphological, and tectonic controls to the morphological variations of submarine landslides in the active margin. Furthermore, the identified regions with clustered large submarine landslides and the segments of the NSTS provide essential information in modeling fault and submarine failuredriven tsunamis and identification of highly exposed coastal areas. This study also demonstrates the use of morphological and exploratory spatial analyses to elucidate underlying controlling factors and to evaluate the hazard potential of areas with limited geological and geophysical datasets.

Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions
LN and NR contributed to the conceptualization and structure of the study. LN wrote the first draft of the manuscript and performed the terrain and statistical analyses. Both authors contributed to the manuscript writing and revision and have read and approved the submitted version.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.