Original Research ARTICLE
Complex and Cascading Triggering of Submarine Landslides and Turbidity Currents at Volcanic Islands Revealed From Integration of High-Resolution Onshore and Offshore Surveys
- 1National Oceanography Centre, University of Southampton Waterfront Campus, Southampton, United Kingdom
- 2National Oceanography Centre, School of Ocean and Earth Science, University of Southampton, Southampton, United Kingdom
- 3Department of Geography & Environment, University of Southampton, Southampton, United Kingdom
- 4Department of Geography, Durham University, Durham, United Kingdom
- 5Department of Earth and Planetary Sciences, Macquarie University, Sydney, NSW, Australia
- 6School of Environment, University of Auckland, Auckland, New Zealand
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 large-scale 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 m3/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 highlights the often-ignored importance of outburst floods, 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.
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 km3) 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, deep-water 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 high-flux, 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 low-explosivity volcano (Stromboli: Aeolian Archipelago; Casalbore et al., 2010, 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). Morphologically-similar 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 (Clare et al., 2016; Hughes Clarke, 2016); dense river-water plunging (Casalbore et al., 2011, 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.
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 strombolian- and 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).
Figure 1. (A) Location of Vanuatu Arc and (B) Tanna Island. Location and extent of survey area offshore Tanna Island, Vanuatu (C) shown in relation to onshore features. Onshore geomorphological and structural mapping based on Firth et al. (2014) and recreated under Copyright License 4433040166144. Three possible locations of the offshore extent of the Siwi Ring Fracture are annotated in white. Terrestrial photography from Google: Digital Globe. Aerial drone photograph (D) of Siwi River and Sulfur Bay taken toward the north-east at point D annotated on panel C. Satellite photograph (E) taken in 2001, following outburst flood in the Siwi River in 2000. December 2000 coastline is annotate in yellow.
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.
Figure 2. Overview of seafloor topography and main features in the study area. (A) Colourwash bathymetry overlain on greyscale slope map. Terrestrial data from Google: Digital Globe. (B) 3D rendering (3 × vertical exaggeration) of hillshaded bathymetry (illumination from the north-west) annotated with main geomorphologic features and north-arrow (green).
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 cm3) 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. Aliquot samples (1 g) were dispersed in 30 ml 0.05% sodium hexametaphosphate solution and shaken for 24 h. Dispersed aliquots were analyzed using a Malvern Mastersizer 2000 using laser diffraction of suspended sediment grains (10,000 counts) to measure grain size distributions. Grain size distributions were measured three times per aliquot. Aliquots showed intra-sample variations of <3%. Standard reference materials showed intra-sample variations of up to 3% and accuracy toward reference values of 1.5%. Scanning Electron Microscopy (SEM) was performed using a Hitachi TM-1000 Microscope at the British Ocean Sediment Core Research Facility (BOSCORF) on selected samples to investigate micro-textural properties of the sediments.
Offshore Morphology of the Yenkahe Resurgent Dome, Tanna Island
Analysis of bathymetric data (ground-truthed based on observations made from video footage acquired from the ROV) generally reveals a smooth low-lying seafloor (carbonate platform) or a rougher, textured seafloor (fringing or patch 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.
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).
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).
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).
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 down-slope of an arcuate bight-like feature on the northern flank of the carbonate platform in Sulfur Bay (Figure 3D). The largest of these blocks has an estimated volume of 9.29 × 10−6 km3, but all were significantly smaller than the scar from which they potentially originated. Three larger (<3.82 × 10−4 km3) angular blocky features form localized prominent positive relief that deflect the course of seafloor channels (Figures 3B,E,F, 4, 5). 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 km3. The two debris aprons have estimated volumes of 1.28 × 10−3 km3 and 9.67 × 10−4 km3, based on their planform area and assuming that they have a wedge-shaped cross-sectional geometry [in line with the approach used by McAdoo et al. (2000); Chaytor et al. (2009)].
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 down-slope 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.
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).
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 back-angled 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 bight-like features, tilted blocks, nor the coalescent sinuous channels to the west of Port Resolution so no comment can be made on the seafloor sediments in those areas.
Figure 7. Location of sediment samples presented in this study (A) overlain on drone-acquired photographic survey of Siwi River and multibeam bathymetry offshore. Grain size distributions (B) from Siwi River (gray lines) are similar to those from the proximal part of the sinuous channel with crescentic bedforms (black lines – 1–3), but distinct from samples from linear gullies (yellow lines). Details of SEM and transmitted light microscopy are shown for specific locations (C,D).
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, 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 km3). 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 small-scale slope collapses. Similar piecemeal failures of carbonate-dominated 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 large-magnitude 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 caldera-modifying 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).
Figure 8. Photograph (A) from survey vessel of steep coastal cliffs at location annotated on (C). Photograph taken using ROV (B) showing boulders (~30–50 cm in diameter) below the cliffs. Satellite photography prior to (C,D) and immediately following (E,F) tropical cyclone Pam. (D,F) zoomed in to Siwi River and Sulfur Bay. Map images acquired by NASA and exported from Google: Digital Globe.
Table 1. Review of possible triggers for turbidity currents and submarine landslides at volcanic islands, with specific reference to documented events at Tanna Island (parenthesized numbers are cross-references to Figure 11).
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 ROV-video 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.
Figure 9. Comparison of satellite photography and bathymetry of Resolution Bay acquired in 2017, compared with outline of coastline from mapping by Captain Belcher in 1840 in yellow (Hydrographic Office of the Admiralty, 1843). Filled circles illustrate comparison between bathymetric soundings in 1840 and 2017, which show an average elevation difference of +6.8 m; equating to an annual average rise of 40 mm/year.
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 (Clare et al., 2016; 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 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).
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.
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 m3 of water over only 2 days, triggering catastrophic erosion of more than 1.1 million m3 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/m3 (Wilson et al., 2012). A 73% freshwater (1,000 kg/m3) and 27% sediment mixture equates to a flow density of 1,365–1,446 kg/m3. 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/m3 above seawater (~1,030 kg/m3) 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 m3/s is estimated, with an upper bound (99th percentile) limit of 7,000 m3/s (Figure 11). Turbidity currents (up to ~4 m/s) triggered by plume settling (with a density surfeit of <1 kg/m3 above seawater) have been directly observed to occur frequently offshore from bedload-dominated rivers at a discharge threshold of >250 m3/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 m2, which equates to a bankfull discharge of between 210 m3/s and 630 m3/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 channel offshore Sulfur Bay, which are similar to those associated with jokulhaups with observed peak discharges of ~1,000 m3/s; Russell, 1993).
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 m3/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 m3/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 m3/s, with an upper bound estimate of 7,000 m3/s.
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 × 106 m3; Manville, 2010). Caldera lakes can be larger still, such as Lake Toba in Indonesia (2.4 × 1011 m3) or Lake Taupo in New Zealand (6 × 1010 m3, where an outburst flood had an estimated peak discharge of 17,000–35,000 m3/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), 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 cyclone-triggered turbidity currents (Pope et al., 2017). Multiple powerful cyclones have been documented in Vanuatu in recent years, including tropical cyclones Uma in 1987 A.D., Fran in 1993 A.D., Prema in 1993 A.D., Paula in 2001 A.D. and Ivy in 2004 A.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 sediment-laden 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 m3/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 run-off 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.
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 base-level 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). 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).
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.
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).
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 m3/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.
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.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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 Douglas Kiri of the Vanuatu Environmental Science Society. We also thank the Department of Geology and Mines for support with fieldwork and logistics on Tanna. The staff of the BOSCORF are acknowledged for assistance with SEM analysis. We acknowledge NERC funding grants NE/M007138/1, NE/M017540/1, NE/P009190/1 and NE/P005780/1. MJB is supported by a Royal Society Dorothy Hodgkin Fellowship (RF1504449).
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/
Babonneau, N., Delacourt, C., Cancouët, R., Sisavath, E., Bachèlery, P., Mazuel, A., et al. (2013). Direct sediment transfer from land to deep-sea: insights into shallow multibeam bathymetry at La Réunion Island. Mar. Geol. 346, 47–57. doi: 10.1016/j.margeo.2013.08.006
Brothelande, E., Lénat, J. F., Normier, A., Bacri, C., Peltier, A., Paris, R., et al. (2015). Insights into the evolution of the Yenkahe resurgent dome (Siwi caldera, Tanna Island, Vanuatu) inferred from aerial high-resolution photogrammetry. J. Volcanol. Geother. Res. 322, 212–224. doi: 10.1016/j.jvolgeores.2015.04.006
Brothelande, E., Peltier, A., Got, J. L., Merle, O., Lardy, M., and Garaebiti, E. (2016). Constraints on the source of resurgent doming inferred from analogue and numerical modeling—Implications on the current feeding system of the Yenkahe dome–Yasur volcano complex (Vanuatu). J. Volcanol. Geother. Res. 322, 225–240. doi: 10.1016/j.jvolgeores.2015.11.023
Calder, E. S., Cole, P. D., Dade, W. B., Druitt, T. H., Hoblitt, R. P., Huppert, H. E., et al. (1999). Mobility of pyroclastic flows and surges at the Soufriere Hills Volcano, Montserrat. Geophys. Res. Lett. 26, 537–540. doi: 10.1029/1999GL900051
Caminade, P., Charlie, D., Kanoglu, U., Koshimura, S. I., Matsutomi, H., Moore, A., et al. (2000). Vanuatu earthquake and tsunami cause much damage, few casualties. EOS Trans. Am. Geophys. Union 81, 641–647. doi: 10.1029/EO081i052p00641-02
Carey, S., Morelli, D., Sigurdsson, H., and Bronto, S. (2001). Tsunami deposits from major explosive eruptions: an example from the 1883 eruption of Krakatau. Geology 29, 347–350. doi: 10.1130/0091-7613(2001)029andlt;0347:TDFMEEandgt;2.0.CO;2
Carlino, S., Cubellis, E., Luongo, G., and Obrizzo, F. (2006). On the mechanics of caldera resurgence of Ischia Island (southern Italy). Geol. Soc. Lond. Spec. Public. 269, 181–193. doi: 10.1144/GSL.SP.2006.269.01.12
Casalbore, D., Chiocci, F. L., Mugnozza, G. S., Tommasi, P., and Sposato, A. (2011). Flash-flood hyperpycnal flows generating shallow-water landslides at Fiumara mouths in Western Messina Strait (Italy). Mar. Geophys. Res. 32:257. doi: 10.1007/s11001-011-9128-y
Casalbore, D., Ridente, D., Bosman, A., and Chiocci, F. L. (2017). Depositional and erosional bedforms in Late Pleistocene-Holocene pro-delta deposits of the Gulf of Patti (southern Tyrrhenian margin, Italy). Mar. Geol. 385, 216–227. doi: 10.1016/j.margeo.2017.01.007
Casalbore, D., Romagnoli, C., Bosman, A., and Chiocci, F. L. (2014). Large-scale seafloor waveforms on the flanks of insular volcanoes (Aeolian Archipelago, Italy), with inferences about their origin. Mar. Geol. 355, 318–329. doi: 10.1016/j.margeo.2014.06.007
Casalbore, D., Romagnoli, C., Chiocci, F., and Frezza, V. (2010). Morpho-sedimentary characteristics of the volcaniclastic apron around Stromboli volcano (Italy). Mar. Geol. 269, 132–148. doi: 10.1016/j.margeo.2010.01.004
Chen, J. K., Taylor, F. W., Edwards, R. L., Cheng, H., and Burr, G. S. (1995). Recent emerged reef terraces of the Yenkahe resurgent block, Tanna, Vanuatu: implications for volcanic, landslide and tsunami hazards. J. Geol. 103, 577–590. doi: 10.1086/629777
Chouet, B. A. (1996). “New methods and future trends in seismological volcano monitoring,” in Monitoring and Mitigation of Volcano Hazards, eds R. Scarpa and R. I. Tilling (Berlin; Heidelberg: Springer), 23–97.
Clare, M., Chaytor, J., Dabson, O., Gamboa, D., Georgiopoulou, A., Eady, H., et al. (2018). A consistent global approach for the morphometric characterization of subaqueous landslides. Geol. Soc. Lond. Spec. Public 477, SP477–15. doi: 10.1144/SP477.15
Clare, M. A., Clarke, J. H., Talling, P. J., Cartigny, M. J. B., and Pratomo, D. G. (2016). Preconditioning and triggering of offshore slope failures and turbidity currents revealed by most detailed monitoring yet at a fjord-head delta. Earth Planetary Sci. Lett. 450, 208–220. doi: 10.1016/j.epsl.2016.06.021
Clare, M. A., Vardy, M. E., Cartigny, M. J., Talling, P. J., Himsworth, M. D., Dix, J. K., et al. (2017). Direct monitoring of active geohazards: emerging geophysical tools for deep-water assessments. Near Surface Geophys. 15, 427–444. doi: 10.31223/osf.io/m3gkj
Comte, I., Colin, F., Whalen, J. K., Grünberger, O., and Caliman, J.-P. (2012). “Agricultural practices in oil palm plantations their impact on hydrological changes, nutrient fluxes water quality in indonesia: a review,” in Advances in Agronomy, ed D. L. Sparks (Burlington: Academic Press), 71–124.
Counts, J. W., Jorry, S. J., Leroux, E., Miramontes, E., and Jouet, G. (2018). Sedimentation adjacent to atolls and volcano-cored carbonate platforms in the Mozambique Channel (SW Indian Ocean). Mar. Geol. 404, 41–59. doi: 10.1016/j.margeo.2018.07.003
Coussens, M., Wall-Palmer, D., Talling, P., Watt, S., Cassidy, M., Jutzeler, M., et al. (2016). The relationship between eruptive activity, flank collapse, and sea level at volcanic islands: a long-term (> 1 Ma) record offshore Montserrat, Lesser Antilles. Geochem. Geophys. Geosyst. 17, 2591–2611. doi: 10.1002/2015GC006053
Covault, J. A., Kostic, S., Paull, C. K., Sylvester, Z., and Fildani, A. (2017). Cyclic steps and related supercritical bedforms: building blocks of deep-water depositional systems, western North America. Mar. Geol. 393, 4–20. doi: 10.1016/j.margeo.2016.12.009
Cronin, S. J., Gaylord, D. R., Charley, D., Alloway, B. V., Wallez, S., and Esau, J. W. (2004). Participatory methods of incorporating scientific with traditional knowledge for volcanic hazard management on Ambae Island, Vanuatu. Bull. Volcanol. 66, 652–668. doi: 10.1007/s00445-004-0347-9
Cronin, S. J., Neall, V. E., Lecointre, J. A., and Palmer, A. S. (1997). Changes in Whangaehu River lahar characteristics during the 1995 eruption sequence, Ruapehu volcano, New Zealand. J. Volcanol. Geother. Res. 76, 47–61. doi: 10.1016/S0377-0273(96)00064-9
Cronin, S. J., and Sharp, D. S. (2002). Environmental impacts on health from continuous volcanic activity at Yasur (Tanna) and Ambrym, Vanuatu. Int. J. Environ. Health Res. 12, 109–123. doi: 10.1080/09603120220129274
Delmelle, P., Henley, R. W., and Bernard, A. (2015). “Volcano-related lakes,” in The Encyclopedia of Volcanoes, 2nd Edn., eds H. Sigurdsson, B. Houghton, S. McNutt, H. Rymer and J. Stix (San Diego, CA: Elsevier) 851–864.
Duller, R. A., Mountney, N. P., Russell, A. J., and Cassidy, N. C. (2008). Architectural analysis of a volcaniclastic jökulhlaup deposit, southern Iceland: sedimentary evidence for supercritical flow. Sedimentology 55, 939–964. doi: 10.1111/j.1365-3091.2007.00931.x
Firth, C. W., Cronin, S. J., Turner, S. P., Handley, H. K., Gaildry, C., and Smith, I. (2015). Dynamics and pre-eruptive conditions of catastrophic, ignimbrite-producing eruptions from the Yenkahe Caldera, Vanuatu. J. Volcanol. Geother. Res. 308, 39–60. doi: 10.1016/j.jvolgeores.2015.10.012
Firth, C. W., Handley, H. K., Cronin, S. J., and Turner, S. P. (2014). The eruptive history and chemical stratigraphy of a post-caldera, steady-state volcano: Yasur, Vanuatu. Bull. Volcanol. 76:837. doi: 10.1007/s00445-014-0837-3
Gavey, R., Carter, L., Liu, J. T., Talling, P. J., Hsu, R., Pope, E., et al. (2017). Frequent sediment density flows during 2006 to 2015, triggered by competing seismic and weather events: observations from subsea cable breaks off southern Taiwan. Mar. Geol. 384, 147–158. doi: 10.1016/j.margeo.2016.06.001
Gombiner, J. H., Hemming, S. R., Hendy, I. L., Bryce, J. G., and Blichert-Toft, J. (2016). Isotopic and elemental evidence for Scabland Flood sediments offshore Vancouver Island. Quat. Sci. Rev. 139, 129–137. doi: 10.1016/j.quascirev.2016.02.026
Graham, I. J. A. G., Reyes, I. C., Wright, K. M., Peckett, I. E. M., Smith, R. J., and Arculus (2008) Structure petrology of newly discovered volcanic centers in the northern Kermadec–southern Tofua arc, South Pacific Ocean, J. Geophys. Res. 113:B08SB02. doi: 10.1029/2007JB005453
Gudmundsson, M. T., Thordarson, T., Höskuldsson, A., Larsen, G., Björnsson, H., Prata, F. J., et al. (2012). Ash generation and distribution from the April-May 2010 eruption of Eyjafjallajökull, Iceland. Sci. Rep. 2:572. doi: 10.1038/srep00572
Hage, S., Cartigny, M. J., Clare, M. A., Sumner, E. J., Vendettuoli, D., Hughes Clarke, J. E., et al. (2018). How to recognize crescentic bedforms formed by supercritical turbidity currents in the geologic record: insights from active submarine channels. Geology 46, 563–566. doi: 10.1130/G40095.1
Hizzett, J. L., Hughes Clarke, J. E., Sumner, E. J., Cartigny, M. J. B., Talling, P. J., and Clare, M. A. (2018). Which triggers produce the most erosive, frequent, and longest runout turbidity currents on deltas?. Geophys. Res. Lett. 45, 855–863. doi: 10.1002/2017GL075751
Hodgson, D. M., Bernhardt, A., Clare, M. A., da Silva, A. C., Fosdick, J., Mauz, B., et al. (2018). Grand challenges (and Great Opportunities) in sedimentology, stratigraphy, and diagenesis research. Front. Earth Sci. 6:173. doi: 10.3389/feart.2018.00173
Hoffmann, G., Silver, E., Day, S., Morgan, E., Driscoll, N., and Orange, D. (2008). Sediment waves in the Bismarck volcanic arc, Papua New Guinea. Spec. Papers Geol. Soc. Am. 436:491. doi: 10.1130/2008.2436(05)
Hong, I., Pilarczyk, J. E., Horton, B. P., Fritz, H. M., Kosciuch, T. J., Wallace, D. J., et al. (2018). Sedimentological characteristics of the 2015 tropical cyclone pam overwash sediments from Vanuatu, South Pacific. Mar. Geol. 396, 205–214. doi: 10.1016/j.margeo.2017.05.011
Hunt, J. E., Wynn, R. B., Talling, P. J., and Masson, D. G. (2013). Multistage collapse of eight western Canary Island landslides in the last 1.5 Ma: sedimentological and geochemical evidence from subunits in submarine flow deposits. Geochem. Geophys. Geosyst. 14, 2159–2181. doi: 10.1002/ggge.20138
International Cable Protection Committee (2016). Submarine Cables and BBNJ. Available online at: http://www.un.org/depts/los/biodiversity/prepcom_files/ICC_Submarine_Cables_&_BBNJ_August_2016.pdf
Jazi, S. D., and Wells, M. (2018). Settling-driven convection limits the spatial scale of deposition beneath sediment-laden buoyant flows in lakes and the coastal ocean. EarthArXiv. doi: 10.31223/osf.io/9xymn
Kosciuch, T. J., Pilarczyk, J. E., Hong, I., Fritz, H. M., Horton, B. P., Rarai, A., et al. (2018). Foraminifera reveal a shallow nearshore origin for overwash sediments deposited by Tropical Cyclone Pam in Vanuatu (South Pacific). Mar. Geol. 396, 171–185. doi: 10.1016/j.margeo.2017.06.003
Kroon, F. J., Kuhnert, P. M., Henderson, B. L., Wilkinson, S. N., Kinsey-Henderson, A., Abbott, B., et al. (2012). River loads of suspended solids, nitrogen, phosphorus and herbicides delivered to the Great Barrier Reef lagoon. Mar. Pollut. Bull. 65, 167–181. doi: 10.1016/j.marpolbul.2011.10.018
Kudrass, H. R., Michels, K. H., Wiedicke, M., and Suckow, A. (1998). Cyclones and tides as feeders of a submarine canyon off Bangladesh. Geology 26, 715–718. doi: 10.1130/0091-7613(1998)026andlt;0715:CATAFOandgt;2.3.CO;2
Kuleshov, Y. L., Qui, L., Fawcett, R., and Jones, D. (2008). On tropical cyclone activity in the Southern Hemisphere: trends and the ENSO connection. Geophys. Res. Lett. 35, 1–5. doi: 10.1029/2007GL032983
Lavigne, F., and Thouret, J. C. (2003). Sediment transportation and deposition by rain-triggered lahars at Merapi Volcano, Central Java, Indonesia. Geomorphology 49, 45–69. doi: 10.1016/S0169-555X(02)00160-5
Lavigne, F., Thouret, J. C., Voight, B., Young, K., LaHusen, R., Marso, J., et al. (2000). Instrumental lahar monitoring at Merapi Volcano, Central Java, Indonesia. J. Volcanol. Geother. Res. 100, 457–478. doi: 10.1016/S0377-0273(00)00151-7
Leat, P. T., Day, S. J., Tate, A. J., Martin, T. J., Owen, M. J., and Tappin, D. R. (2013). Volcanic evolution of the South Sandwich volcanic arc, South Atlantic, from multibeam bathymetry. J. Volcanol. Geother. Res. 265, 60–77. doi: 10.1016/j.jvolgeores.2013.08.013
Lee, T. Y., Huang, J. C., Lee, J. Y., Jien, S. H., Zehetner, F., and Kao, S. J. (2015). Magnified sediment export of small mountainous rivers in Taiwan: chain reactions from increased rainfall intensity under global warming. PLoS ONE 10:e0138283. doi: 10.1371/journal.pone.0138283
Lintern, D. G., Hill, P. R., and Stacey, C. (2016). Powerful unconfined turbidity current captured by cabled observatory on the Fraser River delta slope, British Columbia, Canada. Sedimentology 63, 1041–1064. doi: 10.1111/sed.12262
Liu, J. T., Wang, Y. H., Yang, R. J., Hsu, R. T., Kao, S. J., Lin, H. L., et al. (2012). Cyclone-induced hyperpycnal turbidity currents in a submarine canyon. J. Geophys. Res. 117:C4. doi: 10.1029/2011JC007630
Lonergan, L., Jamin, N. H., Jackson, C. A. L., and Johnson, H. D. (2013). U-shaped slope gully systems and sediment waves on the passive margin of Gabon (West Africa). Mar. Geol. 337, 80–97. doi: 10.1016/j.margeo.2013.02.001
Manville, V., White, J. D. L., Houghton, B. F., and Wilson, C. J. N. (1999). Paleohydrology and sedimentology of a post−1.8 ka breakout flood from intracaldera Lake Taupo, North Island, New Zealand. Geol. Soc. Am. Bull. 111, 1435–1447. doi: 10.1130/0016-7606(1999)111andlt;1435:PASOAPandgt;2.3.CO;2
Maria, A., Carey, S., Sigurdsson, H., Kincaid, C., and Helgadóttir, G. (2000). Source and dispersal of jokulhlaup sediments discharged to the sea following the 1996 Vatnajokull eruption. Geol. Soc. Am. Bull. 112, 1507–1521. doi: 10.1130/0016-7606(2000)112andlt;1507:SADOJKandgt;2.0.CO;2
Meheux, K., and Parker, E. (2006). Tourist sector perceptions of natural hazards in Vanuatu and the implications for a small island developing state. Tour. Manage. 27, 69–85. doi: 10.1016/j.tourman.2004.07.009
Merle, O., Brothelande, E., Lénat, J. F., Bachèlery, P., and Garaébiti, E. (2013). A structural outline of the Yenkahe volcanic resurgent dome (Tanna Island, Vanuatu Arc, South Pacific). J. Volcanol. Geother. Res. 268, 64–72. doi: 10.1016/j.jvolgeores.2013.10.009
Milliman, J. D., Snow, J., Jaeger, J., and Nittrouer, C. A. (1996). Catastrophic discharge of fluvial sediment to the ocean: evidence of Jökulhlaups events in the Alsek Sea Valley, southeast Alaska (USA). IAHS Publications Series Proc. Rep. Int. Assoc. Hydrol. Sci. 236, 367–380.
Mitchell, N. C., Masson, D. G., Watts, A. B., Gee, M. J., and Urgeles, R. (2002). The morphology of the submarine flanks of volcanic ocean islands: a comparative study of the Canary and Hawaiian hotspot islands. J. Volcanol. Geother. Res. 115, 83–107. doi: 10.1016/S0377-0273(01)00310-9
Moore, J. G., Clague, D. A., Holcomb, R. T., Lipman, P. W., Normark, W. R., and Torresan, M. E. (1989). Prodigious submarine landslides on the Hawaiian Ridge. J. Geophys. Res. 94, 17465–17484. doi: 10.1029/JB094iB12p17465
Mulder, T., Syvitski, J. P., Migeon, S., Faugeres, J. C., and Savoye, B. (2003). Marine hyperpycnal flows: initiation, behavior and related deposits. A review. Mar. Petrol. Geol. 20, 861–882. doi: 10.1016/j.marpetgeo.2003.01.003
Nomikou, P., Carey, S., Papanikolaou, D., Bell, K. C., Sakellariou, D., Alexandri, M., et al. (2012). Submarine volcanoes of the Kolumbo volcanic zone NE of Santorini Caldera, Greece. Glob Planet Change 90, 135–151. doi: 10.1016/j.gloplacha.2012.01.001
Normandeau, A., Lajeunesse, P., Poiré, A. G., and Francus, P. (2016). Morphological expression of bedforms formed by supercritical sediment density flows on four fjord-lake deltas of the south-eastern Canadian Shield (Eastern Canada). Sedimentology 63, 2106–2129. doi: 10.1111/sed.12298
Parsons, J. D., Bush, J. W., and Syvitski, J. P. (2001). Hyperpycnal plume formation from riverine outflows with small sediment concentrations. Sedimentology 48, 465–478. doi: 10.1046/j.1365-3091.2001.00384.x
Partin, J. W., Quinn, T. M., Shen, C. C., Emile-Geay, J., Taylor, F. W., Maupin, C. R., et al. (2013). Multidecadal rainfall variability in South Pacific Convergence Zone as revealed by stalagmite geochemistry. Geology 41, 1143–1146. doi: 10.1130/H34718.1
Paull, C. K., Talling, P. J., Maier, K. L., Parsons, D., Xu, J., Caress, D. W., et al. (2018). Powerful turbidity currents driven by dense basal layers. Nat. Commun. 9:4114. doi: 10.1038/s41467-018-06254-6
Pope, E. L., Jutzeler, M., Cartigny, M. J., Shreeve, J., Talling, P. J., Wright, I. C., et al. (2018). Origin of spectacular fields of submarine sediment waves around volcanic islands. Earth Planet. Sci. Lett. 493, 12–24. doi: 10.1016/j.epsl.2018.04.020
Pope, E. L., Talling, P. J., Carter, L., Clare, M. A., and Hunt, J. E. (2017). Damaging sediment density flows triggered by tropical cyclones. Earth Planet. Sci. Lett. 458, 161–169. doi: 10.1016/j.epsl.2016.10.046
Puga-Bernabéu, Á., Webster, J. M., Beaman, R. J., and Guilbaud, V. (2013). Variation in canyon morphology on the Great Barrier Reef margin, north-eastern Australia: the influence of slope and barrier reefs. Geomorphology 191, 35–50. doi: 10.1016/j.geomorph.2013.03.001
Quartau, R., Ramalho, R. S., Madeira, J., Santos, R., Rodrigues, A., Roque, C., et al. (2018). Gravitational, erosional and depositional processes on volcanic ocean islands: Insights from the submarine morphology of Madeira Archipelago. Earth Planet. Sci. Lett. 482, 288–299. doi: 10.1016/j.epsl.2017.11.003
Rothwell, R. G., Hoogakker, B., Thomson, J., Croudace, I. W., and Frenz, M. (2006). Turbidite emplacement on the southern Balearic Abyssal Plain (western Mediterranean Sea) during Marine Isotope Stages 1–3: an application of ITRAX XRF scanning of sediment cores to lithostratigraphic analysis. Geol. Soc. Lond. Spec. Public 267, 79–98. doi: 10.1144/GSL.SP.2006.267.01.06
Russell, A. J. (1993). Obstacle marks produced by flow around stranded ice blocks during a glacier outburst flood (jökulhlaup) in west Greenland. Sedimentology 40, 1091–1111. doi: 10.1111/j.1365-3091.1993.tb01381.x
Silver, E., Day, S., Ward, S., Hoffmann, G., Llanes, P., Driscoll, N., et al. (2009). Volcano collapse and tsunami generation in the Bismarck volcanic arc, Papua New Guinea. J. Volcanol. Geother. Res. 186, 210–222. doi: 10.1016/j.jvolgeores.2009.06.013
Stephens, S. A., and Ramsay, D. (2014). Extreme cyclone wave climate in the Southwest Pacific Ocean: influence of the El Niño Southern Oscillation and projected climate change. Glob. Planet Change 123, 13–26. doi: 10.1016/j.gloplacha.2014.10.002
Stow, D. A., Hernández-Molina, F. J., Llave, E., Sayago-Gil, M., Díaz del Río, V., and Branson, A. (2009). Bedform-velocity matrix: the estimation of bottom current velocity from bedform observations. Geology 37, 327–330. doi: 10.1130/G25259A.1
Symons, W. O., Sumner, E. J., Talling, P. J., Cartigny, M. J., and Clare, M. A. (2016). Large-scale sediment waves and scours on the modern seafloor and their implications for the prevalence of supercritical flows. Mar. Geol. 371, 130–148. doi: 10.1016/j.margeo.2015.11.009
Tappin, D. R., Watts, P., McMurtry, G. M., Lafoy, Y., and Matsumoto, T. (2001). The Sissano, Papua New Guinea tsunami of July 1998—offshore evidence on the source mechanism. Mar. Geol. 175, 1–23. doi: 10.1016/S0025-3227(01)00131-1
Terry, J. P., Garimella, S., and Kostaschuk, R. A. (2002). Rates of floodplain accretion in a tropical island river system impacted by cyclones and large floods. Geomorphology 42, 171–182. doi: 10.1016/S0169-555X(01)00084-8
Terry, J. P., Kostaschuk, R. A., and Wotling, G. (2008). Features of tropical cyclone-induced flood peaks on Grande Terre, New Caledonia. Water Environ. J. 22, 177–183. doi: 10.1111/j.1747-6593.2007.00098.x
Terry, J. P., and Lau, A. A. (2018). Magnitudes of nearshore waves generated by tropical cyclone Winston, the strongest landfalling cyclone in South Pacific records. Unprecedented or unremarkable? Sediment. Geol. 364, 276–285. doi: 10.1016/j.sedgeo.2017.10.009
Toomey, M., Donnelly, J., and Woodruff, J. (2013). Reconstructing mid-late Holocene cyclone variability in the Central Pacific using sedimentary records from Tahaa, French Polynesia. Quat. Sci. Rev. 77, 181–189. doi: 10.1016/j.quascirev.2013.07.019
Vanuatu Ministry of Lands and Natural Resources (2014). Available online at: https://web.archive.org/web/20171016152101/https://mol.gov.vu/index.php/en/water/223-the-lake-isiwi-disaster-sulfur-bay-tanna (Accessed November 5, 2018).
Watson, S. J., Whittaker, J. M., Lucieer, V., Coffin, M. F., and Lamarche, G. (2017). Erosional and depositional processes on the submarine flanks of Ontong Java and Nukumanu atolls, western equatorial Pacific Ocean. Mar. Geol. 392, 122–139. doi: 10.1016/j.margeo.2017.08.006
Willems, B. A., Powell, R. D., Cowan, E. A., and Jaeger, J. M. (2011). Glacial outburst flood sediments within Disenchantment Bay, Alaska: implications of recognizing marine jökulhlaup deposits in the stratigraphic record. Mar. Geol. 284, 1–12. doi: 10.1016/j.margeo.2011.03.004
Wilson, T. M., Stewart, C., Sword-Daniels, V., Leonard, G. S., Johnston, D. M., Cole, J. W., et al. (2012). Volcanic ash impacts on critical infrastructure. Phys. Chem. Earth Parts A B C 45, 5–23. doi: 10.1016/j.pce.2011.06.006
Wright, I. C., Worthington, T. J., and Gamble, J. A. (2006). New multibeam mapping and geochemistry of the 30–35 S sector, and overview, of southern Kermadec arc volcanism. J. Volcanol. Geother. Res. 149, 263–296. doi: 10.1016/j.jvolgeores.2005.03.021
Zhang, Y., Liu, Z., Zhao, Y., Colin, C., Zhang, X., Wang, M., et al. (2018). Long-term in situ observations on typhoon-triggered turbidity currents in the deep sea. Geology 46, 675–678. doi: 10.1130/G45178.1
Keywords: cascading hazards, turbidity current, submarine landslide, tropical cyclone, outburst flood, volcanic island, crescentic bedforms, Strombolian volcano
Citation: Clare MA, Le Bas T, Price DM, Hunt JE, Sear D, Cartigny MJB, Vellinga A, Symons W, Firth C and Cronin S (2018) Complex and Cascading Triggering of Submarine Landslides and Turbidity Currents at Volcanic Islands Revealed From Integration of High-Resolution Onshore and Offshore Surveys. Front. Earth Sci. 6:223. doi: 10.3389/feart.2018.00223
Received: 19 September 2018; Accepted: 21 November 2018;
Published: 13 December 2018.
Edited by:Ivar Midtkandal, University of Oslo, Norway
Reviewed by:Gijs Allard Henstra, University of Bergen, Norway
Miquel Poyatos Moré, University of Oslo, Norway
Copyright © 2018 Clare, Le Bas, Price, Hunt, Sear, Cartigny, Vellinga, Symons, Firth and Cronin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael A. Clare, email@example.com