Complex and Cascading Triggering of Submarine Landslides and Turbidity Currents at Volcanic Islands Revealed From Integration of High-Resolution Onshore and Offshore Surveys

highlights the often-ignored importance of outburst ﬂoods, non-linear responses to land-use and climatic changes, and the complex interactions between a range of coastal and tectonic processes that may overshadow volcanic regimes.

Submerged flanks of volcanic islands are prone to hazards including submarine landslides that may trigger damaging tsunamis and sediment-laden seafloor flows (called "turbidity currents"). These hazards can break seafloor infrastructure which is critical for global communications and energy transmission. Small Island Developing States are particularly vulnerable to these hazards due to their remote and isolated nature, small size, high population densities, and weak economies. Despite their vulnerability, few detailed offshore surveys exist for such islands, resulting in a geohazard "blindspot," particularly in the South Pacific. Understanding how these hazards are triggered is important; however, pin-pointing specific triggers is challenging as most studies have been unable to link continuously between onshore and offshore environments, and focus primarily on largescale eruptions with sudden production of massive volumes of sediment. We address these issues by integrating the first detailed (2 × 2 m) bathymetry data acquired from Tanna Island, Vanuatu with a combination of terrestrial remote sensing data, onshore and offshore sediment sampling, and documented historical events. Mount Yasur on Tanna has experienced low-magnitude Strombolian activity for at least the last 600 years. We find clear evidence for submarine landslides and turbidity currents, yet none of the identified triggers are related to major volcanic eruptions, in contrast to conclusions from several previous studies. Instead we find that cascades of non-volcanic events (including outburst floods with discharges of >1,000 m 3 /s, and tropical cyclones), that may be separated by decades, are more important for preconditioning and triggering of landslides and turbidity currents in oversupplied sedimentary regimes such as at Tanna. We conclude with a general model for how submarine landslides and turbidity currents are triggered at volcanic and other heavily eroding mountainous islands. Our model INTRODUCTION Active volcanic islands can create a variety of subaerial hazards including explosive eruptions that disrupt air transport (e.g., Gudmundsson et al., 2012), emission of gases harmful to health (e.g., Horwell and Baxter, 2006), fast-moving pyroclastic flows and lahars (e.g., Cronin et al., 1997;Calder et al., 1999), and ash falls that destroy agriculture and pollute water supplies (e.g., Wilson et al., 2012). Given their high relief in the surrounding deep ocean (up to 7 km above the surrounding seafloor), the subaerial extents of volcanic islands are typically dwarfed by their submerged flanks (Watt et al., 2014). These submarine slopes are often affected by dynamic sediment transport processes that can also pose a major hazard (Watt et al., 2014). Subsea flank collapses can be prodigious in scale (>100 s of km 3 ) and trigger damaging tsunamis (Moore et al., 1989;Keating and McGuire, 2000;Carey et al., 2001;Tappin et al., 2001;Coussens et al., 2016). Strategically important seafloor infrastructure, such as the subsea telecommunications cable network that transmits more than 95% of all digital data traffic worldwide, is vulnerable to submarine landslides or powerful sediment avalanches (called "turbidity currents") that occur offshore from volcanic islands (Carter et al., 2014;Pope et al., 2017). Small Island Developing States are disproportionately vulnerable to both subaerial and submarine hazards; largely due to their remote and isolated nature, small size, high population densities at or near sea-level and weak economies (Briguglio, 1995;Pelling and Uitto, 2001;Cronin et al., 2004;Terry and Goff, 2012;Hodgson et al., 2018). Understanding the threats posed to seafloor cables is particularly important for these islands, as telecommunication links underpin many critical areas for development, including access to regional markets, overseas bank transactions and booking for tourism (ICPC, 2016). Despite their vulnerability, remarkably few detailed offshore surveys exist for Small Island Developing States in the South Pacific; largely due to geographic and economic constraints (CClare et al., 2018). Therefore, the South Pacific has been identified as a "hazard blind-spot" with respect to submarine landslides and associated tsunamis (Terry and Goff, 2012;Goff and Terry, 2016). Furthermore, the link between onshore and offshore sediment transport processes, and hence, identification of the triggers for offshore hazards, is often unclear as integrated subaerial and submarine surveys are limited to relatively few volcanic islands worldwide (e.g., Casalbore et al., 2010;Babonneau et al., 2013).
The advent of modern multibeam bathymetry has enabled detailed imaging of the seafloor, thus providing insights into submarine landslides at volcanic islands (e.g., Mitchell et al., 2002). In addition to evidence of past slope failures, deepwater seafloor surveys (>25 m cell size) of active volcanoes have revealed that previously undocumented, crescentic bedforms are common from shallow to deep-water modern marine volcaniclastic systems (Wright et al., 2006;Hoffmann et al., 2008;Silver et al., 2009;Gardner, 2010;Leat et al., 2013;Pope et al., 2018). These bedforms may be diagnostic of turbidity currents triggered by major, episodic volcanic events, including: (i) sector or flank collapse; (ii) powerful Vulcanian eruptions; or (iii) sustained Plinian eruptions that can produce highflux, sustained pyroclastic density currents (Pope et al., 2018). These past studies focussed on large and powerful scenarios, but what should we expect if and when volcanoes erupt at lower rates or have long-term low but steady outputs? These cases are arguably most common, with large explosive volcanoes undergoing centuries or millennia of quiescence between events and many less-explosive volcanoes having regular small eruptions (e.g., throughout Vanuatu; Tonga, Solomon Islands, Papua New Guinea and New Zealand). Only one existing bathymetric study is known for a regularly erupting lowexplosivity volcano (Stromboli: Aeolian Archipelago; Casalbore et al., 2010Casalbore et al., , 2014, which also reveals similar crescentic bedforms. Thus, it is plausible that crescentic bedforms offshore from volcanic islands may signify sudden catastrophic collapses that originate after long-term preconditioning and by a range of multiple possible triggering mechanisms-directly related to volcanism or otherwise. With smaller individual eruptions known at such sites, crescentic bedforms this bedforms thus may not constitute evidence for major volcanic events. Crescentic bedforms offshore from some volcanic islands have been tentatively linked to turbidity currents triggered by non-volcanic processes, with documented examples including ephemeral sediment-laden river floods that plunge directly into the sea (Babonneau et al., 2013;Quartau et al., 2018). Morphologicallysimilar crescentic bedforms have also been described in many other subaqueous, non-volcanic settings worldwide and have been related to a wide range of triggers, including: subaqueous delta collapses Hughes Clarke, 2016); dense river-water plunging (Casalbore et al., 2011(Casalbore et al., , 2017; settling of sediment from river plumes (Hizzett et al., 2018); wave and storm resuspension (Xu et al., 2004;Normandeau et al., 2016); and glacial outburst floods (Duller et al., 2008). Thus, there is a large degree of ambiguity in linking bedform morphology at volcanic islands with triggering mechanisms. One constraint to our understanding has been the challenge of acquiring detailed multibeam data in water depths of <100 m; hence, few studies have acquired data that extends shallow enough to link the subaerial volcano to offshore and of appropriate resolution to image bedforms (Casalbore et al., 2010;Quartau et al., 2018). There is a pressing need to acquire detailed data offshore from volcanic islands to better understand the nature and triggers of offshore hazards, and the link between terrestrial and marine environments; in particular at Small Island Developing States and offshore from long-lived volcanoes. Here, we present the first detailed survey (2 × 2 m cell size) offshore from Yasur volcano on Tanna Island, Vanuatu in the South Pacific. Typical volcanic activity at Yasur involves low-magnitude Strombolian eruptions (Nairn et al., 1988;Firth et al., 2014), making it an ideal location for this study.

AIMS
We integrate our offshore data with existing and new onshore data to address the following specific aims. First, what is the offshore morphology of a continuously active and rapidly uplifting volcanic island, and what processes caused that morphology? We investigate whether arcuate-bight like features can be linked to slope failure, as suggested previously by Goff and Terry (2016), and whether offshore sediment transport pathways can be identified, such as the trains of crescentic bedforms observed on other volcanic islands. Second, we ask whether submarine landslides and crescentic bedforms found offshore from volcanic islands are always directly linked to major eruptive volcanic activity or flank collapses? We identify possible volcanic and non-volcanic triggers for submarine landslides and turbidity currents offshore Tanna Island based on documented historical events, and through integration of onshore and offshore analysis. Third, we ask how important is the role of cascades of events, which may be separated by decades, compared to instantaneous triggers? Finally, we outline a general model for the preconditioning and triggering of submarine landslides and turbidity currents at volcanic islands worldwide, based on insights from this and other studies.

Geological and Physiographic Setting for Tanna Island, Vanuatu
Tanna is one of 83 islands making up the 1,200 km-long Vanuatu volcanic arc in the south-west Pacific (Brothelande et al., 2016; Figure 1). Vanuatu has been affected by a wide range of natural hazards in recent and historical times, including earthquakes, tropical cyclones, and tsunamis (Meheux and Parker, 2006). Given its vulnerability to these hazards, Vanuatu has been described as the most disaster-prone country in the South Pacific (Meheux and Parker, 2006). Tanna was formed approximately 2.5 Ma by successive episodes of volcanism and reef growth (Carney and Macfarlane, 1979). Volcanism is currently focussed on Yasur; one of the most active volcanoes in the archipelago. Yasur is a scoria cone, which formed from repeated strombolianand vulcanian-style eruptions that occur every few minutes (Nairn et al., 1988;Merle et al., 2013). These are fed by a steady-state magma reservoir, which has been providing basaltic trachy-andesitic magma to drive eruptions at Yasur for at least the last 600 years (Nairn et al., 1988;Merle et al., 2013;Firth et al., 2014). Previous phases of volcanic activity in this vicinity were more dramatic, however. The Siwi ring fracture (Figure 1), defines the previous limit of a compound caldera, which collapsed during at least two major ignimbrite-forming eruptions at approximately 43 and 3-8 ka (Firth et al., 2015). In its lower-most course, the Siwi River drains along the northern edge of the Siwi ring fracture, until it reaches Sulfur Bay where it meets a back barrier-type beach (Figure 1). Until recently, the Siwi River fed into Lake Isiwi, which was dammed by a lava flow that was emplaced prior to 1800 A.D. (Figure 1; Firth et al., 2015). Heavy rainfall in 2000 A.D. led to the breaching of the dam, triggering a major outburst flood that cut a new channel and flowed into Sulfur Bay (Kanas et al., 2000).
The present-day eruptive activity of Mount Yasur shows continuous low-level explosivity; however, shallow magma intrusion drives significant post-caldera uplift on Tanna which may contribute to a range of potential geohazards (Merle et al., 2013). The Yenkahe Resurgent Dome is among the fastest resurgent calderas worldwide (Merle et al., 2013), with uplift rates of 156 mm/year calculated over the last 1,000 years from dating of uplifted coral terraces (Chen et al., 1995; Figure 1). Two strong earthquakes in 1878 A.D. caused up to 12 m of co-seismic uplift at the coast by Port Resolution (Nairn et al., 1988;Merle et al., 2013). Photogrammetric surveys provide possible evidence for several subaerial collapse scars, each with estimated volumes of a few million cubic meters (Brothelande et al., 2015). Some of these potential headscars abut the coastline between Sulfur Bay and Port Resolution, and form steep, often-overhanging cliffs cut into weathered basaltic sands (Brothelande et al., 2015; Figure 1C). Recommendations were made by Brothelande et al. (2015) to perform bathymetric surveys offshore from these features to understand whether such features, and their run-out, extend offshore. Our study focuses specifically on this area (Figures 1, 2) to understand the links between onshore and offshore sediment transport at a dynamic volcanic island through integration with previous land-based studies.

Data
A multibeam survey was performed by EGS Survey on behalf of the UK Hydrographic Office in March 2017. The survey covers an area of ∼6.5 × 3.2 km, and extends from the coastline to 292 m water depth (Figure 2A). Multibeam bathymetry data were acquired using a Kongsberg EM2040 system (200-400 kHz range) and processed into 2 × 2 m bins; hence features smaller than 2 m across cannot be resolved. An onshore photogrammetry survey of the distal part of the Siwi River and the beach at Sulfur Bay (∼900 × 900 m composed from two flights) was performed in October 2017 using a DJI Phantom 4 unmanned aerial vehicle. Pix4Dcapture was used to predefine a flight plan at 100 m altitude. Agisoft Photoscan was used to create an orthomosaic with a pixel size of <4 cm. Offshore sediment sampling was performed using a two-disc grabber-cup (10 cm 3 ) mounted on a small portable Deep Trekker DTG2 Remotely Operated Vehicle (ROV) equipped with an additional high resolution camera (GoPro HERO4 silver) and deployed from the MV Escape (a 12.9 m catamaran) in October 2017. Offshore sediment samples were targeted within linear gullies (seven locations) and a submarine channel (three locations). Onshore sediment samples were hand-excavated from five locations in the Siwi River during the same survey in October 2017. Grain size analysis followed the procedures in Rothwell et al. (2006). Sediment was sieved at 2 mm to remove rare over-sized particles then three aliquots of each sub-sample were taken for measuring grain size. FIGURE 3 | Overview map (A) annotated with panels that illustrate evidence for localized slope instability. Line of sight is illustrated by yellow arrow for the following panels. Slope rendering (B, where black is steepest slope) and seafloor roughness (C, where red is roughest and blue is smoothest) maps show an elongated, tilted block that may be a failed slab of carbonate platform that has subsequently been eroded around. Arcuate bight-like structures may reflect past episodes of slope instability, and feature smaller blocks that appear to have slumped (D). High slope angles (B) and roughness (C) on the flanks of the elongated, tilted block may indicate smaller-scale incipient slope failures, which are annotated in 3D renderings (E,F).
reefs) in shallow (<60 m) waters to the north of Sulfur Bay and offshore from Port Resolution (Figure 2). Unlike its expression onshore, the edge of the Siwi Ring Fracture (i.e., the collapsed caldera margin) is difficult to trace offshore. Three different types of geomorphic character indicative of sediment transport are found locally obscuring the carbonate platform and extending into deeper water. These include: (i) arcuate bight-like features and tilted or displaced blocks ( Figure 3); (ii) linear gullies, which are either isolated or coalescent in form ( Figure 2); and (iii) trains of crescentic bedforms within sinuous channels, locally associated with scours (Figures 4, 5). These three geomorphic characters now form the observational basis of this paper.

Arcuate Bight-Like Features and Tilted Blocks
Steep-flanked (up to 60 • ) arcuate bight-like features were identified locally cutting back into the carbonate platform ( Figure 3D). At least five tilted blocks occur immediately downslope of an arcuate bight-like feature on the northern flank of Frontiers in Earth Science | www.frontiersin.org 6 December 2018 | Volume 6 | Article 223 FIGURE 4 | Overview map (A) and 3D rendering (B) annotated with dashed line along the axis of sinuous channel with cresentic bedforms that originates in Sulfur Bay. Along-channel profile (black and gray lines in C) shows how bedforms generally increase in wavelength and amplitude with increasing water depth. Average slope gradients are annotated on the black profile, while local slope gradient that highlights bedforms and knickpoints is shown as a red line in (C). Localized obstacle and scour features are shown in (D).
the carbonate platform in Sulfur Bay ( Figure 3D). The largest of these blocks has an estimated volume of 9.29 × 10 −6 km 3 , but all were significantly smaller than the scar from which they potentially originated. Three larger (<3.82 × 10 −4 km 3 ) angular blocky features form localized prominent positive relief that deflect the course of seafloor channels (Figures 3B,E,F, 4, 5).
Frontiers in Earth Science | www.frontiersin.org 7 December 2018 | Volume 6 | Article 223 FIGURE 5 | Overview map (A) and 3D rendering (B) annotated with dashed lines along the axis of sinuous coalescent channels with cresentic bedforms that originate ∼2.5 km to the east of Sulfur Bay. Along-channel profiles (black, gray, blue, orange and green lines in C) show how bedforms generally increase in wavelength and amplitude to 1,500 m water depth, but then decrease in response to the constriction of the channel. Average slope gradients are annotated on the black profile, while local slope gradient that highlights bedforms is shown as a red line in (C).
Frontiers in Earth Science | www.frontiersin.org 8 December 2018 | Volume 6 | Article 223 FIGURE 6 | Isolated linear gullies are much shorter in length than coalescent gullies, but attain similar slope angles, which are higher on average than those for channels with crescentic bedforms (A). Local slope measurements along all mapped gullies and channels (B) shows that channels with crescentic bedforms only occur at much lower slope angles (mean of 3 degrees).
These blocky features have a low-angle, tilted upper surface (5-10 • ), which is otherwise similar in seafloor expression to the surrounding flat-lying carbonate platform. Their flanks are often steep (up to 60 • ) and are cut by arcuate incisions, downslope of which apron-like accumulations of roughly textured seafloor debris are found ( Figure 3E). The largest of these blocky features has an estimated volume of 2.29 × 10 −3 km 3 . The two debris aprons have estimated volumes of 1.28 × 10 −3 km 3 and 9.67

Linear Gullies That Lack Bedforms
We observe two types of gully morphology: (i) isolated and (ii) coalescent forms. Isolated linear gullies initiate in water depths of 20-30 m on the steepest slopes in the survey area (Figure 2), with a mean slope of 9 • but can locally reach up to 30-40 • ( Figure 6A). Such slopes are immediately downslope of areas with a limited extent of fringing coral (extending no more than 130-200 m seaward from the high water mark; Figure 2A) and with an abundance of boulders (observed from ROV dives). These gullies are most abundant on the eastern flank of Sulfur Bay; the flanks of the Yenkahe dome which is undergoing most rapid uplift (Figure 1; Chen et al., 1995). Linear gullies are up to 500 m in length, and terminate as slope angles reduce; typically where they intersect sinuous channels. Linear gullies maintain a near-continuous width along their straight course, which ranges from 20 to 60 m. Bedforms are absent from linear gullies. Coalescent gullies initiate at water depths of ∼50 m, to the north-east of Port Resolution, in an area of less dramatic uplift outside of the Yenkahe dome (Figures 1,  2). Similarly to isolated linear gullies, they initiate on slopes of up to 30 • , with a mean gradient of 10 • ( Figure 6A). Unlike isolated linear gullies, coalescent forms become adjoined downslope from their initiation points, in an amphitheater-shaped depression (Figure 2). Another difference is that they initiate >500 m offshore from the high water mark, downslope of a more extensive patch of coral reef (Figure 2A). Bedforms are also absent from these features.

Sinuous Channels With Crescentic Bedforms
One of the most extensive bathymetric features in the survey area is a sinuous channel that initiates as a series of small (8-20 m wavelength, 1-2 m amplitude) bedforms in 30 m water depth, immediately offshore from Siwi River at Sulfur Bay (Figure 4), and extends to the north-east beyond the limits of the survey area (Figure 4). Unlike linear gullies, this channel forms on much lower angle slopes (mean of 3 • ; Figure 6A). The channel contains abundant crescentic bedforms, which generally increase in wavelength and amplitude with increasing water depth, where the channel broadens out (to 200 m) on lower angle slopes (Figure 4). The bedforms show a backstepping asymmetry, featuring steep lee (down-stream) faces and lower angle backangled stoss (up-stream) faces ( Figure 4C). The channel is also punctuated by steeper and deeper scours with gradients of up to 20-30 degrees on their lee (down-stream) face, and obstacle and comet structures (Stow et al., 2009) oriented parallel with the axis of the channel. A series of channels also initiates in water depths of between 40 and 50 m, to the west of Port Resolution. These commence individually as approximately 10 m-wide channels, until they coalesce at approximately 125 m water depth to form one channel that broadens to approximately 130 m ( Figure 5). This combined channel then adjoins the single broad sinuous channel and extends beyond the limits of the survey area ( Figure 5). In common with the broad sinuous channel (Figure 4), these channels feature an abundance of similar back-stepping crescentic bedforms ( Figure 5C). The bedforms generally increase in size with increasing water depth; however, they locally attain lower amplitudes and wavelengths where the channel is constricted or steepened by seafloor relief. Channel orientation is strongly controlled by features that present prominent seafloor relief, such as tilted blocks.

Composition of Seafloor Sediments
Grain size analysis from crescentic bedforms in the submarine channel at Sulfur Bay reveals a very similar distribution to samples from the Siwi River, with clear bimodality at the most proximal location, becoming progressively finer offshore (Figure 7). The grain size within linear gullies is distinctly different to both samples from Siwi River and crescentic bedforms offshore Sulfur Bay, showing a finer and broader grain size distribution. Transmitted light and scanning electron microscopy show that sediment is dominantly comprised of basaltic lithics with a small component of volcanic glass. Small amounts of carbonate and coralline debris are incorporated further offshore. Samples were not taken from the arcuate bightlike features, tilted blocks, nor the coalescent sinuous channels to Frontiers in Earth Science | www.frontiersin.org 9 December 2018 | Volume 6 | Article 223 the west of Port Resolution so no comment can be made on the seafloor sediments in those areas.

DISCUSSION
We first discuss the origin of the bathymetric features observed offshore Tanna and whether they relate to caldera collapse (along ring fractures due to catastrophic volcanic eruption) or offshore sediment transport processes. Second, based on the evidence for slope instability offshore Tanna, we discuss whether similar features on volcanic islands elsewhere in the world represent tsunamigenic collapse of carbonate platforms. Third, we discuss the potential triggers for slope instability and turbidity currents on volcanic islands, initially focusing on the range of plausible triggering events at Tanna. We conclude by proposing a general model for their preconditioning and triggering at volcanic islands worldwide, invoking a complex interplay of both volcanic and non-volcanic processes.

Challenges in Delineating the Offshore Extent of Caldera Margins
Caldera collapses at many volcanic islands have a distinct outer margin; the extent of which can be continuously mapped from onshore to offshore (e.g., Pantelleria, Italy; Nisyros and Santorini, Greece; Deception Island, Antarctica; Rabaul, Papua New Guinea; Aira, Japan -Walker, 1984;Nomikou et al., 2012). Such clarity is not apparent offshore Tanna, however. The submerged outline of the Siwi Ring Fracture is difficult to define (Figure 1). The overall morphological complexity of the caldera at Tanna is probably due to its formation by at least two major caldera-collapse episodes, with other modification possible during smaller intervening eruptions (Firth et al., 2014(Firth et al., , 2015. In addition, the northern margin of the Siwi Ring Fracture is possibly erased by the rapid uplift of the Yenkahe block, with higher uplift rates in the NE possibly "popping-out" the northern caldera margin (Brothelande et al., 2015). Resurgent calderas (pushed piston-like back up along or within caldera fractures) are common in submarine volcanic settings where magma rises into them following collapse (e.g., submarine Tonga arc; Graham et al., 2008), where crust is thin. They are also possible in other areas, e.g., Ischia Island off Naples is mapped as a fully resurgent caldera, where past caldera fill has been uplifted to form a steep island (Carlino et al., 2006). Furthermore, it is likely that the seafloor was strongly affected by a combination of sediment deposition, transport, and slope failures in the period since (or during) caldera-formation. These processes have sculpted and reworked both the caldera margin and carbonate platform. That a feature which is so distinct onshore, can be almost entirely reworked or masked by offshore sediment transport processes, has implications for the recognition and interpretation of partially or entirely submerged caldera collapses in areas of active seafloor sediment transport processes.

Are Arcuate Bight-Like Features Formed by Slope Failures and if So, Were They
Single-Event or Multi-Phase in Nature?
Based on coarse resolution (>100 m) regional bathymetry, Terry and Goff (2013) identified arcuate bight-like features incised into carbonate platforms on a number of volcanic islands and atolls in the South Pacific, and proposed a submarine slope failure for their origin. They further suggested that such events may be very large in volume, and could trigger significant tsunamis if failure occurs in one displacement event. Indeed, many studies of volcanic islands have revealed prodigious-volume landslides of their submerged flanks (e.g., Moore et al., 1989;Masson, 1996;Coussens et al., 2016). The high resolution bathymetry offshore Tanna reveals arcuate bight-like features cut into the carbonate platform. Much of the large-scale (>km) "scalloping" of the carbonate platform is attributed here to caldera collapse, rather than slope failure. There is, however, compelling evidence of smaller-scale submarine slope failure within arcuate bights with perimeter lengths of 100-1,000s of meters. These slope failures are superimposed on the post-caldera collapse relief ( Figure 3D). We interpret the arcuate bights to the north and east of Sulfur Bay as the up-slope limit of collapse events. The tilted blocks found downslope (partially-translated and/or rotated debris) are much smaller than the scars from which they originated (Figure 3). Thus, it is likely that submarine slope failures offshore from Tanna occurred progressively, as multiple phases of small volume collapses and partially-rotated blocks (<2.9 × 10 −3 km 3 ). The heterogeneous nature of the mixed carbonate platform and patch reefs into which these bights are incised presumably results in localized zones of weaker material that may fail preferentially due to erosion, undercutting, or from external cyclic loading (e.g., earthquake, storm waves; Keating and McGuire, 2000). These smaller bights are interpreted to arise from a combination of retrogression and lateral unloading (when adjacent areas of seafloor are removed) during multiple phases of relatively smallscale slope collapses. Similar piecemeal failures of carbonatedominated shelf breaks and slopes are common in the Bahama Banks, Great Barrier Reef and a number of volcanic-cored atolls in the South Indian Ocean (Puga-Bernabéu et al., 2013;Jo et al., 2015;Watson et al., 2017;Counts et al., 2018), suggesting that this situation may be similar for many other atolls and volcanic islands flanked by carbonate platform-reefs. The landslide-origin hypothesis for bight-like features (Terry and Goff, 2013) is generally supported by our findings; however, we suggest that landslide-related bights may form progressively in multiple stages rather than as one event. Multi-stage slope failures typically relate to a much lower tsunami hazard than one-off en-masse collapses, due to the smaller event volumes involved and the time elapsed between displacements (Hunt et al., 2013). Furthermore, without high resolution data, it may be challenging to attribute arcuate bight-like morphology to slope failure rather than caldera collapse. These complexities thus underline the importance of acquiring high resolution multi-beam bathymetry and the value of future efforts to map the offshore regions of Small Island Developing States for local and regional hazard assessments, particularly in the South Pacific.
What Processes Are Responsible for Creating Gullies and Submarine Channels With Crescentic Bedforms?
The linear gullies and seafloor channels observed offshore Tanna (Figures 2, 4, 5) are morphologically similar to those observed in many settings worldwide where density currents transport sediment to deeper waters (e.g., Micallef and Mountjoy, 2011;Babonneau et al., 2013;Lonergan et al., 2013;Symons et al., 2016;Casalbore et al., 2017;Covault et al., 2017). In particular, the scale, morphology and grain-size of the crescentic bedforms within the sinuous channels are very similar (i.e., meters in amplitude, tens of meters in wavelength, fine to coarse sand) to those where repeat seafloor surveys and direct flow monitoring have demonstrated the occurrence of density-stratified turbidity currents that undergo a series of hydraulic jumps (Hughes Clarke, 2016;Normandeau et al., 2016;Hage et al., 2018;Paull et al., 2018). In such sites, flows switch between super-and subcritical flow regimes that drive the up-stream migration of crescentic bedforms (Hughes Clarke, 2016;Normandeau et al., 2016;Hage et al., 2018). So why do linear gullies without bedforms occur, as well as channels with crescentic bedforms? Slope gradient appears to exert a strong control occur, as linear gullies have a significantly higher gradient (mean of 9-10 • ; Figure 6) than channels containing crescentic bedforms (mean of 3 • ). This is in line with observations by Micallef and Mountjoy (2011) who identified a minimum slope threshold (5 • ) for the formation of linear gullies, arguing that a critical bed shear stress can only be attained on such steep slopes. Quartau et al. (2018) only found linear gullies on the volcanic islands of the Madeira archipelago at slopes of >15 • . On the submarine flanks of Stromboli volcano, Casalbore et al. (2010) only observed crescentic bedforms on slopes of <5 • . We suggest therefore that both gullies and crescentic bedforms offshore Tanna were created by turbidity currents and that slope angle dictates the nature of the flow-seafloor interaction and thus the resultant morphology (Kostic, 2011;Zhong et al., 2015). But what processes were responsible for triggering these flows? Were volcanic events solely responsible? We now explore these questions.

Do Gullies and Crescentic Bedforms on Volcanic Islands Only Result From Flows Triggered by Major Volcanic Events?
Trains of crescentic bedforms occur on the submarine flanks of many volcanic islands globally, including Stromboli, Reunion Island, the Canary Islands, and islands in the Bismark, West Mariana, Kermadec, and South Sandwich island arcs (see database in Symons et al., 2016;Pope et al., 2018 and references therein). At most of these sites, it has been inferred that these seafloor features result from major volcanic events: either largemagnitude explosive eruptions, or large flank/sector collapses (Pope et al., 2018). This is highly unlikely to be the case for the features observed offshore Tanna. The most recent Plinian eruption on Tanna occurred ∼3-8 ka, modifying the Siwi Caldera and emplacing widespread ignimbrite deposits (Firth et al., 2015). Onshore, these deposits radiate out from the Siwi Ring Fracture, but are absent within the caldera (Firth et al., 2015). Crescentic bedforms and linear gullies are found within the inferred offshore caldera margin, suggesting that they must post-date the calderamodifying eruption. More recent volcanic activity from Yasur has involved continuous, low magnitude Strombolian and Vulcanian eruptions over at least the last 600 years (Nairn et al., 1988;Chen et al., 1995). This style of activity produces high rates of sediment input into the surrounding areas with ash fall and also contributes to a large, devegetated or sparsely vegetated area downwind of the volcano. Eruptions with major flow events powerful enough to scour gullies onshore Tanna have not been recorded in the recent geological record. Sea cliffs between Sulfur Bay and Port Resolution ( Figure 8A) form the north and eastern margins of the rapidly uplifting Yenkahe Dome. Based on dating of uplifted coral terraces on this block (Figure 1) by Chen et al. (1995), these cliffs have formed over the last 1-2 ka. Sea cliffs to the north of the Siwi River incise into deposits from Plinian eruptions dated at 3-8 and ∼43 ka (Firth et al., 2015), hence must also significantly pre-date the features observed at seafloor. We can therefore rule out these Plinian eruptions as a trigger and conclude that bedforms on the flanks of volcanic islands do not necessarily relate to major volcanic events. If this is the case, then what are the other plausible triggers? We now explore potential mechanisms, initially considering those that are indirectly related to volcanic activity, and then those that are unrelated ( Table 1).

Low Magnitude Volcanic Activity and Related Preconditioning Effects
While the effect of major eruptive volcanic activity may not necessarily be directly responsible, the cumulative or antecedent conditions resulting from past or ongoing low magnitude volcanic activity may play a key role in preconditioning slopes to failure or setting up a successive chain of events that may trigger turbidity currents at volcanic islands. For instance, the accumulation of relatively weak and laminated volcanic deposits, and the effects of shallow hydrothermal circulation, are plausible contributing factors to promote slope instability (Brothelande et al., 2015). Dynamic topographic changes also play a potentially important role. The uplift rates calculated for the Yenkahe Dome, abutting the area between Sulfur Bay and Port Resolution are among the highest for any resurgent dome worldwide (156 mm/year), with two earthquake events in 1878 A.D. leading to up to 10 m co-seismic uplift at Port Resolution (Chen et al., 1995). Elevation differences between a hydrographic survey performed in 1840 A.D. (Hydrographic Office of the Admiralty, 1843) and our 2017 A.D. survey, indicate the seafloor rose by between 2.6 and 11.3 m in Resolution Bay (an area located east of the even more rapidly uplifting Yenkahe Dome). This equates to an average rise of 40 mm /year; however, most of the elevation change was likely due to the 1878 A.D. earthquakes as evidenced by eye witness accounts. The Rev. Lawrie wrote in 1898: "On the island of Tanna there was a great earthquake on the 10th January 1878, which caused a surge of the water at Port Resolution to rise forty feet, and to sweep everything before it, destroying all the canoes of the natives. Two minutes after the earthquake a rise of the land took place on the while west side of the harbour, to the extent of about twenty feet. This narrowed considerably the effective anchorage of the harbour, and a lost anchor came into view where a ship had ridden safely some years previously. About a month afterwards another earthquake caused a further elevation, so that rocks which were formerly covered with seven or eight fathoms of water are now above high-water mark" (Lawrie, 1898). These uplift events caused the subaerial exposure of parts of Resolution Bay, forming the present day Lake Eweya (Figure 9). The coupling of geotechnically weak volcanic deposits and rapid uplift has been used to explain the presence of multiple onshore slope failures (Brothelande et al., 2016) and presumably explains the existence of steep, often-overhanging cliffs between Sulfur Bay and Port Resolution ( Figure 8A). Underwater ROVvideo footage reveals accumulations of boulders and other debris below these cliffs on the carbonate platform ( Figure 8B), downslope of which a series of isolated linear gullies is found. Thus, cliff collapses may transition to sediment-laden density flows as they disaggregate and mix with seawater at the edge of the carbonate platform to create linear gullies.

The Role of Outburst Floods-An Under-appreciated Hazard at Volcanic Islands?
Crescentic bedforms similar to those observed offshore Tanna have been identified at many non-volcanic locations where rivers directly feed submarine canyons or channels (Symons et al., 2016 and references therein). At river-fed locations, it is hypothesized that turbidity currents initiate from a number of possible mechanisms, during, or shortly following periods of elevated river discharge. First, if sediment-laden river water is dense enough it may directly plunge upon entering the sea, initiating a hyperpycnal flow (Mulder et al., 2003). Second, settling from a buoyant sediment-laden river plume may settle more diffusively via a process known as convective fingering, periodically initiating turbidity currents (Parsons et al., 2001;  Hizzett et al., 2018;Jazi and Wells, 2018). Third, sediment delivered by a river flood rapidly accumulates at the river mouth and periodically becomes unstable, thus triggering a delayed slope failure that initiates a turbidity current Hughes Clarke, 2016). Submarine channels with crescentic bedforms occur offshore from river outflows on volcanic islands in La Reunion and the Madeira Archipelago, and have been tentatively linked to ephemeral periods of flash flooding that may trigger hyperpycnal flow (Babonneau et al., 2013;Quartau et al., 2018). This may indicate that river floods are a potential trigger; however, Pope et al. (2018) explicitly ruled out river floods as a potential explanation for the formation of crescentic bedforms on volcanic islands of the Kermadec Arc, on the basis of a small hydrologic system that prevents  Last known pyroclastic-forming eruptions occurred at 3-8 ka and 34 ka (Firth et al., 2015) so unlikely to be responsible for features observed on present day seafloor.
Flank/dome/sector collapse [4] Turbidity current evolves from slope failure or lava dome collapse due to eruption or growth of magma chamber.
Last major sector collapse was 34 ka so unlikely responsible for features on present day seafloor. No major collapses have occurred in recent times, but local slope failure may have been triggered due to fast uplift rates (156 mm/year, averaged over 2 ka) for the Yenkahe resurgent dome (Chen et al., 1995;Firth et al., 2015).

INDIRECTLY RELATED TO VOLCANIC ACTIVITY
Crater lake or lava-dammed lake outburst flood [5] Sudden discharge of heavily sediment-laden outwash [A] plunges to directly trigger hyperpycnal flow, [B] indirectly triggers turbidity current due to settling out of sediment from a sediment-laden plume, or [C] leads to rapid accumulation of sediment at shelf break setting up delayed slope failure.
Heavy rainfall in 2000 triggered outburst flood from Lake Isiwi releasing 4.1 million m 3 of water and erosion of 1.1 million m 3 of sediment that cut a new channel and flowed out to sea at Sulfur Bay (Kanas et al., 2000;Firth et al., 2014). Estimated peak discharge of >1,000 m 3 /s ( Figure 11).

Lahar
[6] Heavy rainfall mobilizes volcanic sediments and washes them offshore and triggers turbidity currents (in the same manner as 5A-C above).
No recorded evidence, but may occur due to tropical cyclone [8].
Earthquake/ground movement [7] Ground shaking and/or uplift triggers slope or cliff collapse of weathered subaerial volcanic sediments or at the steep fringes of carbonate platform.
Strong ground movement (from witnesses) and co-seismic uplift of up to 12 m due to two earthquakes in 1878 (Chen et al.,1995).

Tsunami or storm surge [8] Loading by waves triggers [A] slope instability or cliff collapse or [B] resuspension of seafloor sediments.
7.1-7.5 M w earthquake triggered a tsunami on November 26th 1999 with 6.6 m run-up height on Tanna. Tsunamis also recorded in 1875, 1961, 1965(Caminade et al., 2000.  .
Onshore failures/cliff collapse [10] Onshore slope failures or collapse of coastal cliffs enters the sea triggered by climatic, erosion, or other non-volcanic processes (Terry and Goff, 2013).
Two historic landslides (1919 and 1975) occurred on the slopes of Yasur's cone but did not enter the sea (Carney and Macfarlane, 1979;Merle et al., 2013). Steep overhanging sea cliffs cut into weathered basalt may represent the headscars of subaerial failures (Brothelande et al., 2015) prone to incipient failure.
large-scale fluvial output to the ocean. At Tanna, we observe a submarine channel with crescentic bedforms offshore from the Siwi River, so are the two systems linked here? Grain size analysis points to a connection, given the similarities between samples hand-excavated from the river and those acquired from the submarine channel using an ROV (Figure 7). However, analysis of satellite photography since 2001, and our new aerial photography, does not indicate a river plume (with the exception of the aftermath of Tropical Cyclone Pam as discussed in the following section), and the river discharge is generally very low (based on visual observations and the presence of a back barrier at the river outflow). Therefore, it is unlikely that the background river discharge is capable of directly triggering turbidity currents (Figures 8C,D).
A series of events that culminated in May 2000 A.D. offers a likely mechanism for the submarine morphology (Figure 10). In 2000 A.D., above average rainfall triggered an outburst flood from Lake Isiwi, which was previously impounded by a tephra barrier on top of a lava flow (Kanas et al., 2000;Vanuatu Ministry of Lands and Natural Resources, 2014). The lake had no regular surface water outlet. The effect of heavy rainfall was exacerbated by the loss of storage capacity, caused by lake-wide deposition of over 1 m (average of 2.3 m) thickness of sediments that were eroded from the flanks of Mount Yasur and the upper reaches of the catchment during tropical cyclone Uma in 1987 (Kanas et al., 2000). The Vanuatu Ministry of Lands and Natural Resources provide a summary of observations from islanders: "Water began overflowing the corner of the lake closest to Sulfur Bay (northeast). On morning of May 2, Provincial Officials and the Local Police warned the residents of Sulfur Bay to move to higher ground to avoid flooding. Once the lake began overflowing its banks it quickly eroded through the soft volcanic deposits situated downgradient of the Lake. At approximately 6:00 pm on May 2 a huge volume of water began flowing from the Lake toward Sulphur Bay. The sound of the rushing water could be heard for miles. Villagers alerted by the roar of the flood waters, ran for higher ground. Luckily no lives were lost. Approximately 10 houses and a Nakamal were destroyed. An unknown number of livestock were reported to have been swept into the sea" (Vanuatu Ministry of Lands and Natural Resources, 2014). The outburst flood ran approximately along the course of the Siwi River, cutting a new channel of up to 40 m depth, until it reached the ocean at Sulfur Bay; washing several cattle and ten houses out to sea (Kanas et al., 2000). The village at Sulfur Bay was severely damaged by the outflow and was buried with up to 1.5 m of sediment (Kanas et al., 2000). Satellite data reveal that the back-barrier was breached, but reformed within at least a year ( Figure 1E).
As the outburst flood discharged to the ocean at the outflow of the Siwi River, it is a highly plausible candidate for creating the crescentic bedforms in the submarine channel offshore of Sulfur Bay. So what were the likely flow conditions? The flood released 4.1 million m 3 of water over only 2 days, triggering catastrophic erosion of more than 1.1 million m 3 of sediment (Kanas et al., 2000;Firth et al., 2014). Based on these values, the flow averaged sediment concentration may have been up to 27% by volume (i.e., hyperconcentrated flow). This concentration is not unreasonable in light of estimates for glacial outbursts ("jokulhaups"; Russell, 1993;Duller et al., 2008), and direct measurements of subaerial debris flows (up to 60%- Weirich, 1989) and remobilised tephra lahars (up to 62% -Cronin et al., 1997;Lavigne and Thouret, 2003) triggered by heavy rainfall. Sediments from the flood, visited by Shane Cronin only months after the event were sand-dominated levees alongside the river, consistent with hyperconcentrated flow. No debris-flow deposits were found. Microscope and SEM analysis of river sediments indicate dominantly basaltic lithics with some volcanic glasses that have a density of 2,350-2,650 kg/m 3 (Wilson et al., 2012). A 73% freshwater (1,000 kg/m 3 ) and 27% sediment mixture equates to a flow density of 1,365-1,446 kg/m 3 . Thus, should the flow have maintained this concentration when it entered the ocean, the density of the flow would have far exceeded the 40 kg/m 3 above seawater (∼1,030 kg/m 3 ) required for hyperpycnal flow (Mulder et al., 2003). As a result the flow could have plunged directly to trigger a turbidity current; as observed from subaerial debris flows and lahars which transform into turbidity currents that last for many hours (Weirich, 1989;Mulder et al., 2003). Lahars can maintain hyperconcentrated conditions for over 40 km (Cronin et al., 1997), thus it is possible that the flow maintained this density to the coastline. We do not have a record of oceanographic conditions during the time of the flood, hence, it is unclear as to whether waves would have inhibited and/or dispersed the plunging of sediment-laden water.
Even if the flow had deposited much of the suspended sediment prior to reaching the ocean at Sulfur Bay, or it was partially redistributed as a homopycnal plume, it is likely that this outburst flood could still have triggered a turbidity current at lower concentrations, particularly given its discharge. Analysis of a global collation of outburst floods (that also includes jokulhaups, artificial dam, and moraine bursts) identified a power-law relationship between volume of water released and peak discharge (Manville, 2010 ; Figure 11). On the basis of the volume released from Lake Isiwi, a peak discharge of 1,000 m 3 /s is estimated, with an upper bound (99th percentile) limit of 7,000 m 3 /s (Figure 11). Turbidity currents (up to ∼4 m/s) triggered by plume settling (with a density surfeit of <1 kg/m 3 above seawater) have been directly observed to occur frequently offshore from bedload-dominated rivers at a discharge threshold of >250 m 3 /s and form similar bedforms (Bornhold et al., 1994;Clare et al., 2016;Hughes Clarke, 2016). Therefore, at its peak, the estimated discharge value for the outburst flood on Tanna is more than that required for triggering turbidity currents. The discharge of the turbidity current itself is more challenging to estimate, however. The cross-sectional area of the submarine channels proximal to the river mouth is ∼70 m 2 , which equates to a bankfull discharge of between 210 m 3 /s and 630 m 3 /s assuming velocities of turbidity currents based on measurements at locations with similar-scale bedforms (3 m/s at Squamish Delta -Hughes Clarke, 2016; 9 m/s at Fraser Delta - Lintern et al., 2016). These estimates may be supported by the localized presence of comet and tail scoured features within the sinuous FIGURE 10 | Sequential chain of cascading effects (from 1 to 15) that may have led to the triggering of turbidity currents that formed crescentic bedforms offshore Sulfur Bay, Tanna Island. channel offshore Sulfur Bay, which are similar to those associated with jokulhaups with observed peak discharges of ∼1,000 m 3 /s; Russell, 1993).
Evidence for outburst floods is increasingly being identified on volcanic islands, where craters, calderas, or past lava flows trap water without a surface outlet (Manville, 2010;Delmelle et al., 2015). These floods are only exceeded in discharge volume by the breaching of glacial impoundments, which are the largest known terrestrial floods on Earth (Manville, 2010). Intracaldera lakes have been identified from more than 100 Holocene volcanoes, with similar water storage volumes to that of the impounded Lake Isiwi (1-10 × 10 6 m 3 ; Manville, 2010). Caldera lakes can be larger still, such as Lake Toba in Indonesia (2.4 × 10 11 m 3 ) or Lake Taupo in New Zealand (6 × 10 10 m 3 , where an outburst flood had an estimated peak discharge of 17,000-35,000 m 3 /s in 232 AD; Manville et al., 1999), and result in far greater discharges than estimated for the outburst flood in 2000 A.D. on Tanna (Figure 11). While several studies have focused on the marine records of jokulhaups (e.g., Milliman et al., 1996;Maria et al., 2000;Willems et al., 2011;Gombiner et al., 2016), FIGURE 11 | Plot of peak discharge in relation to volume of water released from outburst floods based on global data from Manville (2010;reproduced under Copyright License 4433040862707. Black dashed lined is power-law trend for all outburst floods, and light gray parallel fill is the 99th percentile range. Lettered annotations refer to peak river discharges for: (A) Gaoping River following Typhoon Morakot in 2009 [28,000 m 3 /s; (Carter et al., 2014)] when hyperpycnal flow triggered a turbidity current; (B) highest recorded river discharge in New Caledonia on Grand Terre due to Cyclone Anne in 1988 [4,583 m 3 /s; (Terry et al., 2008)]; and (C and dark gray fill) the range of discharge on the Squamish River, British Columbia when turbidity currents are known to occur Clare et al. (2016). The estimated discharge for the 2,000 Lake Isiwi outburst flood is ∼1,000 m 3 /s, with an upper bound estimate of 7,000 m 3 /s. to our knowledge, none have studied the offshore effects of non-glacial outburst floods at volcanic islands. Given the high discharges involved, we suggest that such outburst floods may be an under-appreciated hazard and a potentially important mechanism for initiating long run-out turbidity currents at many volcanic islands.
Triggers Unrelated to Volcanic Activity: Land Use, Extreme Weather Events and the Role of Climate Change While volcanic processes may often be indirectly responsible, a number of non-volcanic events are also capable of providing the sediment discharges and preconditioning required for submarine landslides and turbidity currents to occur. Changes in land cover resulting from human activities in coastal tropical catchments substantially increase suspended sediment loads to the coastal zone (e.g., by 5.5 times, Kroon et al., 2012) and may dramatically increase the likelihood of hillslope failures and other terrestrial landslides (Froude and Petley, 2018). Historically, plantation growth and clearance by the arrival of humans on Pacific islands has increased sediment delivery from river systems as vegetation cover was disturbed by burning and cropping practices. Similarly, changes from forest cover to plantations or agriculture increase storm runoff (Comte et al., 2012). Persistent volcanism exacerbates this, attested by tens of square kilometers of de-vegetated areas in areas affected by volcanic ash and acid rains (Cronin and Sharp, 2002) in downwind areas of volcanoes such as Yasur and Ambrym in Vanuatu.
Tropical cyclones are an important type of non-volcanic event that can enhance preconditioning or directly trigger submarine landslides or turbidity currents due to: (i) storm wave-induced resuspension of shelf sediments; (ii) cyclic loading of unstable slope sediments; (iii) undercutting of coastal cliffs; or (iv) extreme rainfall triggering sediment-laden river floods and surface water run-off that discharge to the ocean (Kudrass et al., 1998;Liu et al., 2012;Carter et al., 2014;Pope et al., 2017). Recent analysis of a global database of telecommunications cable breaks revealed that the Pacific Ocean is a hotspot for tropical cyclonetriggered turbidity currents . Multiple powerful cyclones have been documented in Vanuatu in recent years, including tropical cyclones Uma in 1987A.D., Fran in 1993A.D., Prema in 1993A.D., Paula in 2001A.D. and Ivy in 2004A.D. (Kosciuch et al., 2018. Most recently, tropical cyclone Pam (13th March 2015) made landfall on Tanna Island, traveling at up to 270 km/h with up to 5.3 m-high storm surges, resulting in up to $449M USD in damages Kosciuch et al., 2018. The direct impact of storm waves by events such as Pam is a further plausible explanation for downslope submarine sediment transport in the zone between Sulfur Bay and Port Resolution. Retrogressive or undercutting erosion of the steep coastal cliffs on Tanna by both storm waves and surface water run-off could result in cliff collapse and seaward transport of sediment; perhaps explaining the downslope location of linear gullies and coalescent sinuous channels. The power of such events is demonstrated by tropical cyclones on Fiji that were capable of eroding and transporting carbonate boulders (weighing up to 61 tons; Terry and Lau, 2018). Enhanced turbidity and the presence of a sedimentladen plume was visible around Sulfur Bay and Port Resolution from satellite photography in the days following tropical cyclone Pam, and the beach at Sulfur Bay was eroded landward by tens of meters (Figures 8E,F). Enhanced river outflow also caused breaching of the barrier at the mouth of the Siwi River. River discharges following tropical cyclones can be orders of magnitude higher than background conditions, such as Cyclone Anne in 1988 A.D. which triggered a peak river discharge of more than 4,500 m 3 /s on Grand Terre in New Caledonia (Figure 11; Terry et al., 2008). Sediment loads during tropical cyclone floods have been recorded at 200-500 g/l in Fiji (Terry et al., 2002); far exceeding normal concentrations. Thus, it is likely that tropical cyclone Pam may also have contributed to, or triggered a turbidity current that formed or modified bedforms in the submarine channel initiating in Sulfur Bay (Figure 10). The increasing frequency of El-Nino-Southern Oscillation (ENSO) cycles due to climate change appear to be modifying the intensity of tropical cyclones, their migration tracks, and slowing the rate of their passage, which will result in increased surface water runoff and river discharge (Emanuel, 2005;Kossin et al., 2014;Lee et al., 2015;Mei and Xie, 2016;Chand et al., 2017;Gavey et al., 2017;Pope et al., 2017). Thus, we may expect such events to be a more likely preconditioning and/or triggering mechanism for submarine landslides and turbidity currents offshore from volcanic islands, at least in the Pacific.
FIGURE 12 | General model of processes that contribute to preconditioning or drive triggering of submarine slope failures and turbidity currents at volcanic islands, related directly (A), indirectly (B) and unrelated to volcanic processes. Unlike in Figure 10, the numbering here is not sequential and simply refers to isolated processes.

A General Model for Triggering Submarine Landslides and Turbidity Currents at Volcanic Islands
We found that a single triggering mechanism is often unlikely for submarine landslides and turbidity currents offshore from volcanic islands, and instead that a combination of preconditioning and triggering processes is responsible (Figure 10). On Tanna for example, the series of cascading events that commenced with the closure of Lake Isiwi by a lava flow (pre-1800 A.D.), was compounded by the rising of lake baselevel due to sediment in-wash during tropical cyclone Uma in 1987 A.D., and culminated in the flushing of sediment to the ocean following the outburst flood triggered by elevated rainfall in 2000 A.D (Kanas et al., 2000). The outburst flood contributed to construction of the beach at Sulfur Bay, which was then eroded in 2015 A.D. during tropical cyclone Pam (Figure 10).
FIGURE 13 | Illustration of the potential complexity of interacting processes at volcanic islands that may precondition and trigger submarine landslides and turbidity currents. Plus signs and arrowed lines indicate how an increase in a variable may make a subsequent process more likely. Figure shows scenarios where volcanic factors may dominate (dark gray) and where climatic or anthropogenic factors may be more important (light gray). An animated version of this figure is available in the online material. Yellow circle refers to ultimately triggered event (landslide or turbidity current).
While separated by years to decades in time, these events each served to sequentially modify baseline conditions, setting up a cascade of hazards (Gill and Malamud, 2016). Similarly, the two earthquakes in 1878 A.D. that co-seismically uplifted sea cliffs by up to 12 m (Nairn et al., 1988;Merle et al., 2013), made them steeper, and more prone to wave erosion during severe storms and tropical cyclones in the following decades. Cascading or compound effects of volcanic, climatic and anthropogenic factors should therefore not be overlooked for the triggering of slope failures and turbidity currents offshore from volcanic islands. Land cover and climate changes, in particular, are relatively slow processes that change the background state of the land surface and runoff regime, and may be punctuated by extreme events such as cyclones, earthquakes and eruptions. The result is likely to be a non-linear response over time for given individual or multiple drivers for increased sediment delivery to the coast. Positive-phase Interdecadal Pacific Oscillation (IPO) countries such as Vanuatu, Fiji, and Samoa lie within the South Pacific Convergence Zone under normal conditions (Partin et al., 2013). During El-Nino, however, IPO-positive phase regions experience markedly increased tropical cyclone activity (Kuleshov et al., 2008;Toomey et al., 2013;Stephens and and Ramsay, 2014). If periods of increased land use change (e.g., deforestation), or volcanic eruptions occur coincident with future enhanced ENSO and tropical cyclone intensity [as is predicted for Vanuatu and other Pacific SIDs (Partin et al., 2013;Stephens and and Ramsay, 2014)], we posit that the compounded increase in sediment loads from rivers discharging to the coastal zone will create hotspots for turbidity current generation.
We now conclude with a general model of processes that may contribute to preconditioning and instantaneous triggering of submarine landslides and turbidity currents at volcanic islands (Figure 12), that includes: (A) Processes that are directly related to volcanic activity that are mostly attributed to major eruptive or collapse events, e.g., Plinian eruptions (Manville et al., 1999;Pope et al., 2018); (B) Processes that are indirectly-related to volcanic activity that mostly relate to the preconditioning effects of past volcanism, or progressive ongoing low-magnitude events that may typify quiescent volcanoes and steady-state low-explosivity centers; and (C) Processes that are unrelated to volcanic activity, which include oceanographic and extreme weather events that can affect any type of volcanic island, but are most pronounced in tropical oceans. Each of these processes may play a contributing role in instantaneous triggering, or may continue to precondition the system to enhance the likelihood of offshore sediment transport; hence understanding those interplays is important. The role of cascading hazards may be much more important than that attributed to instantaneous events; particularly for volcanoes under constant low-explosivity conditions where climatic, oceanographic and anthropogenic processes may dominate. We highlight the potentially complex interrelationships between different processes in Figure 13 (animated examples of feedback loops are shown in online Video S1). Because of these compound and/or cascading relationships, attempting to identify one specific triggering mechanism for submarine landslides or turbidity currents is challenging, and may be impossible in many cases. Therefore, determining links between triggers and offshore sediment transport requires careful integration of onshore and marine datasets, and may require direct monitoring of changes in onshore environmental baselines as well as offshore sediment transport processes. Such monitoring is challenging, but new technologies now enable measurement of both the environmental conditions and seafloor processes, thus opening up new opportunities to better understand these complex links to improve offshore hazard assessments in Small Island Developing States (Chouet, 1996;Lavigne et al., 2000;Clare et al., 2017;Zhang et al., 2018).

CONCLUSIONS
We presented the first detailed (2 × 2 m) bathymetric data acquired offshore from Tanna Island, Vanuatu and identified evidence for submarine slope failure and seafloor turbidity currents. These data, coupled with sediment sampling, help to address important knowledge gaps concerning seafloor hazards at Small Island Developing States in the South Pacific, and more generally on the flanks of Strombolian volcanoes, both of which are under-represented in the literature. We found that arcuate bight-like features, incised into the carbonate and reef platform, can be linked to slope collapses that occurred in multiple phases, and thus pose a lower tsunami hazard than if they occurred as one-off, larger failures. Integration of onshore and offshore surveys, with documented historical events, enabled identification of a number of potential triggers for slope failures and turbidity currents offshore Tanna. None of these triggers are related to major volcanic eruptions or collapses, in contrast to conclusions from several previous studies. One highly plausible triggering event was an outburst flood with an estimated discharge of >1,000 m 3 /s. We suggest that outburst floods from crater lakes, caldera lakes and lava flow-impounded features may be under-recognized triggers at many other volcanic islands. Non-volcanic processes, such as tropical cyclones, were also identified as a plausible trigger for triggering slope collapses and turbidity currents, due to storm loading and elevated river discharge to the sea. Tropical cyclones may become more important triggers at islands such as Tanna, due to global warming-induced changes to the El-Nino Southern Oscillation. Finally, we presented a general model for the triggering of submarine landslides and turbidity currents at volcanic islands, underlining the often-ignored importance of non-volcanic processes, and the complex interactions between a range of processes that may precondition the system. We propose that compounding effects, and cascading chains of events, may be more important than instantaneous triggers in many volcanic islands; particularly those in quiescent or Strombolian regimes.

AUTHOR CONTRIBUTIONS
MC led the research, wrote the manuscript and created the figures. All authors contributed to discussions to form the basis of this paper, and provided feedback on the manuscript and figures. TL and DP performed offshore and onshore fieldwork. JH performed grain size, SEM and grain size analysis. DS assisted with analysis of outburst floods and links to river hydrology. MC, AV, and WS assisted with morphodynamic analysis. CF and SC assisted with analysis of the Lake Isiwi outburst flood and interpretation of volcanic aspects of the landscape.

ACKNOWLEDGMENTS
Fieldwork and analysis was supported by the Commonwealth Marine Economies Program which aims to enable safe and sustainable marine economies across Commonwealth Small Island Developing States. We thank EGS Survey Ltd, Justin Jenkin and the crew of the MV Escape, Camillia Garae of the Vanuatu Department of Geology and Mines, and

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.

2018.00223/full#supplementary-material
Video S1 | Conceptual interplay of volcanic and non-volcanic processes on the preconditioning and triggering of offshore landslides and turbidity currents illustrated using a simple process interaction model. Model created using code at https://ncase.me/loopy/