On the 7 August 2019, a 195 km2 raft of andesitic pumice was produced at 200 m below sea level at an unnamed submarine volcano in the Tonga Islands (Southwest Pacific Ocean). Drifting chiefly ...westward, the raft reached the Fiji Islands on the 19 September. Yachts that crossed the raft as early as 2 days post‐eruption provided an outstanding data set of raft characteristics and pristine samples. Further, exceptional tracking of raft dispersal by satellite images allows us to contrast virtual particle tracking methods with ocean model currents to explore the relative influence of surface currents, wind, and wave action on pumice flotsam dispersal over up to 2 years. Attenuation of ocean waves by large and compact pumice rafts appears to reduce the effect of Stokes drift. The coupling of real‐time satellite observations with oceanographic Lagrangian simulations allows near‐real time forecasting for global maritime hazard mitigation.
Plain Language Summary
Although 70% of volcanism occurs underwater, submarine eruptions are seldom witnessed owing to their remoteness and absence of technologies enabling their detection. On the 7 August 2019, a submarine volcanic eruption 200 m below sea level in the Tonga Islands formed a 195 km2 pumice raft. The raft was first encountered and sampled by yacht crews, and is visible on satellite images. Pumice rafts are dispersed by ocean currents, wind, and waves, and can cross entire oceans or get stranded on coasts. After 7 weeks of dispersal, the pumice raft reached the Fiji Islands, and a fraction continued its westward route toward Vanuatu and eastern Australia. Excellent imaging of the raft by satellites permitted reconstruction of daily raft dispersal for the first 8 weeks. Here we show that drift calculations, including components of ocean current, wind, and wave action, can usefully forecast raft dispersal. We tested and tuned these drift calculations by comparing the simulated drift with daily satellite images of the raft. Further, a combination of satellite images and drift calculations based on oceanographic models were used for maritime hazard mitigation in near‐real time.
Key Points
A >30 million m3 pumice raft formed at an unnamed volcano in Tonga midday on the 7 August 2019 and was dispersed by ocean currents and wind
Forecasts of pumice raft dispersal, based on particle drift with model currents, winds, and waves, are evaluated with daily satellite images
Stokes drift does not improve short‐term predictions on large and compact pumice rafts because these rafts attenuate ocean waves
•High resolution seismic data of volcanic island bedforms is presented.•Bedforms related to eruption-fed flows migrated up-slope.•Diagnostic criteria distinguishing eruption-fed and landslide ...bedforms are proposed.•Subaerial and submarine bedform contrasts relate to eruption-fed flow dynamics.
Understanding how large eruptions and landslides are recorded by seafloor morphology and deposits on volcanic island flanks is important for reconstruction of volcanic island history and geohazard assessment. Spectacular fields of bedforms have been recognised recently on submerged flanks of volcanic islands at multiple locations worldwide. These fields of bedforms can extend over 50 km, and individual bedforms can be 3 km in length and 150 m in height. The origin of these bedform fields, however, is poorly understood. Here, we show that bedforms result from eruption-fed supercritical density flows (turbidity currents) in some locations, but most likely rotational landslides at other locations. General criteria are provided for distinguishing between submarine bedforms formed by eruptions and landslides, and emphasise a need for high resolution seismic datasets to prevent ambiguity. Bedforms associated with rotational landslides have a narrower source, with a distinct headscarp, they are more laterally confined, and internal bedform structure does not suggest upslope migration of each bedform. Eruption-fed density currents produce wide fields of bedforms, which extend radially from the caldera. Internal layers imaged by detailed seismic data show that these bedforms migrated up-slope, indicating that the flows that produced them were Froude supercritical. Due to the low density contrast between interstitial fluid and sediment, the extent and dimensions of submarine eruption-fed bedforms is much greater than those produced by pyroclastic density currents on land.
Pumice rafts are floating mobile accumulations of low-density pumice clasts generated by silicic volcanic eruptions. Pumice in rafts can drift for years, become waterlogged and sink, or become ...stranded on shorelines. Here we show that the pumice raft formed by the impressive, deep submarine eruption of the Havre caldera volcano (Southwest Pacific) in July 2012 can be mapped by satellite imagery augmented by sailing crew observations. Far from coastal interference, the eruption produced a single >400 km(2) raft in 1 day, thus initiating a gigantic, high-precision, natural experiment relevant to both modern and prehistoric oceanic surface dispersal dynamics. Observed raft dispersal can be accurately reproduced by simulating drift and dispersal patterns using currents from an eddy-resolving ocean model hindcast. For future eruptions that produce potentially hazardous pumice rafts, our technique allows real-time forecasts of dispersal routes, in addition to inference of ash/pumice deposit distribution in the deep ocean.
A long-standing conceptual model for deep submarine eruptions is that high hydrostatic pressure hinders degassing and acceleration, and suppresses magma fragmentation. The 2012 submarine rhyolite ...eruption of Havre volcano in the Kermadec arc provided constraints on critical parameters to quantitatively test these concepts. This eruption produced a >1 km3 raft of floating pumice and a 0.1 km3 field of giant (>1 m) pumice clasts distributed down-current from the vent. We address the mechanism of creating these clasts using a model for magma ascent in a conduit. We use water ingestion experiments to address why some clasts float and others sink. We show that at the eruption depth of 900 m, the melt retained enough dissolved water, and hence had a low enough viscosity, that strain-rates were too low to cause brittle fragmentation in the conduit, despite mass discharge rates similar to Plinian eruptions on land. There was still, however, enough exsolved vapor at the vent depth to make the magma buoyant relative to seawater. Buoyant magma was thus extruded into the ocean where it rose, quenched, and fragmented to produce clasts up to several meters in diameter. We show that these large clasts would have floated to the sea surface within minutes, where air could enter pore space, and the fate of clasts is then controlled by the ability to trap gas within their pore space. We show that clasts from the raft retain enough gas to remain afloat whereas fragments from giant pumice collected from the seafloor ingest more water and sink. The pumice raft and the giant pumice seafloor deposit were thus produced during a clast-generating effusive submarine eruption, where fragmentation occurred above the vent, and the subsequent fate of clasts was controlled by their ability to ingest water.
•Havre magma entered the ocean before fragmenting.•Clasts were produced by quenching buoyant magma in the ocean.•Buoyant >1 m diameter pumice blocks floated to the ocean surface.•Clasts with enough isolated porosity and trapped gas floated in a raft while the rest sank.
Transition from effusive to explosive volcanism is common during subaerial eruptions, and here we demonstrate that this behavior is also possible underwater. The pyroclastic facies produced ...underwater are distinctive and can be used to distinguish subaqueous from subaerial eruptions and depositional settings. The Pliocene Dogashima Formation (Izu Peninsula, Japan) is a pumice-rich succession originally deposited in an open-marine, below wave-base setting (Jutzeler et al., 2014a). A thick, clast-supported, gray andesite breccia composed of very coarse, dense andesite clasts with quenched margins was sourced from disintegration of an active lava dome. Overall, the gray andesite breccia is gradationally to sharply overlain by thick, graded, clast-supported white pumice breccia chiefly composed of angular pumice clasts and free broken crystals. Regional setting and distinctive facies show that this succession was produced by a fully underwater, magmatic volatile-driven, pumice-forming explosive eruption. The gradational contact between the two breccias, compositional similarities, rare mingled clasts, and fluidal textures in the gray andesite clasts suggest that the explosive eruption destroyed a hot lava dome and generated an eruption-fed, high-concentration density current. In most places, the coarsest hot lava dome fragments were deposited first, followed by the lower density white pumice clasts. The low amount of fine (<2mm) components, well-developed hydraulically controlled grading and sorting, clast angularity, and very coarse dome-derived clasts, some including well-defined quenched margins and common fluidal textures, distinguish the products of subaqueous effusive-to-explosive eruptions from their subaerial counterparts.
•We record subaqueous effusive-to-explosive behavior in the Dogashima Fm, Japan.•Products of subaqueous effusive-to-explosive eruptions are characterized.•Products from subaqueous and subaerial effusive-to-explosive eruptions are distinct.
The International Ocean Discovery Program Expedition 350 drilled between two Izu rear‐arc seamount chains at Site U1437 and recovered the first complete succession of rear‐arc rocks. The drilling ...reached 1806.5 m below seafloor. In situ hyaloclastites, which had erupted before the rear‐arc seamounts came into existence at this site, were recovered in the deepest part of the hole (∼15–16 Ma). Here it is found that the composition of the oldest rocks recovered does not have rear‐arc seamount chain geochemical signatures, but instead shows affinities with volcanic front or some of the extensional zone basalts between the present volcanic front and the rear‐arc seamount chains. It is suggested that following the opening of the Shikoku back‐arc Basin, Site U1437 was a volcanic front or a rifting zone just behind the volcanic front, and was followed at ~ 9 Ma by the start of rear‐arc seamount chains volcanism. This geochemical change records variations in the subduction components with time, which might have followed eastward moving of hot fingers in the mantle wedge and deepening of the subducting slab below Site U1437 after the cessation of Shikoku back‐arc Basin opening.
The IODP Expedition 350 drilled between two Izu rear‐arc seamount chains at Site U1437, which reached 1806.5 m below seafloor and recovered in situ hyaloclastites of ∼15–16 Ma. The composition of the oldest rocks shows affinities with volcanic front. This geochemical change records variations in the subduction components with time, which might have followed eastward moving of hot fingers in the mantle wedge and deepening of the subducting slab below Site U1437 after the cessation of Shikoku Basin opening.
Most of Earth’s volcanic eruptions occur underwater, and these submarine eruptions can significantly impact large-scale Earth systems (e.g., enhancing local primary production by phytoplankton). ...However, detecting submarine eruptions is challenging due to their remote locations, short eruption durations, lack of sea surface signature (if eruptions do not breach the surface), and the transient nature of the surface manifestations of an eruption (e.g., floating pumice clasts, hydrothermal fluids). We can utilize global satellite imagery of 10–30 m resolution (e.g., Landsat 8, Sentinel-2) to detect new eruptions; however, the large data volumes make it challenging to systematically analyze satellite imagery globally. In this study, we address these challenges by developing a new semi-automated analysis framework to detect submarine eruptions through supervised classification of satellite images on Google Earth Engine. We train our algorithm using images from rafts produced by the August 2019 eruption of Volcano F in the Tofua Arc and present a case study using our methodology on satellite imagery from the Rabaul caldera region in Papua New Guinea. We potentially find a large number of new unreported pumice rafts (in ∼16% of images from 2017–present). After analysis of the spatial pattern of raft sightings and ancillary geophysical and visual observations, we interpret that these rafts are not the result of a new eruption. Instead, we posit that the observed rafts represent remobilization of pumice clasts from previous historical eruptions. This novel process of raft remobilization may be common at near-shore/partially submarine caldera systems (e.g., Rabaul, Krakatau) and may have significant implications for new submarine eruption detection and volcanic stratigraphy.
Silicic effusive eruptions in deep submarine environments have not yet been directly observed and very few modern submarine silicic lavas and domes have been described. The eruption of Havre caldera ...volcano in the Kermadec arc in 2012 provided an outstanding database for research on deep submarine silicic effusive eruptions because it produced 15 rhyolite (70–72 wt.% SiO2) lavas and domes with a total volume of ∼0.21 km3 from 14 separate seafloor vents. Moreover, in 2015, the seafloor products were observed, mapped and sampled in exceptional detail (1-m resolution) using AUV Sentry and ROV Jason2 deployed from R/V Roger Revelle. Vent positions are strongly aligned, defining NW-SE and E-W trends along the southwestern and southern Havre caldera margin, respectively. The alignment of the vents suggests magma ascent along dykes which probably occupy faults related to the caldera margin. Four vents part way up the steeply sloping southwestern caldera wall at 1,200–1,300 m below sea level (bsl) and one on the caldera rim (1,060 m bsl) produced elongate lavas. On the steep caldera wall, the lavas consist of narrow tongues that have triangular cross-section shapes. Two of the narrow-tongue segments are connected to wide lobes on the flat caldera floor at ∼1,500 m bsl. The lavas are characterized by arcuate surface ridges oriented perpendicular to the propagation direction. Eight domes were erupted onto relatively flat sea floor from vents at ∼1,000 m bsl along the southern and southwestern caldera rim. They are characterized by steep margins and gently convex-up upper surfaces. With one exception, the domes have narrow spines and deep clefts above the inferred vent positions. One dome has a relatively smooth upper surface. The lavas and domes all consist of combinations of coherent rhyolite and monomictic rhyolite breccia. Despite eruption from deep-water vents (most >900 m bsl), the Havre 2012 rhyolite lavas and domes are very similar to subaerial rhyolite lavas and domes in terms of dimensions, volumes, aspect ratio, textures and morphology. They show that lava morphology was strongly controlled by the pre-existing seafloor topography: domes and wide lobes formed where the rhyolite was emplaced onto flat sea floor, whereas narrow tongues formed where the rhyolite was emplaced on the steep slopes of the caldera wall.
Hole U1395B, drilled southeast of Montserrat during Integrated Ocean Drilling Program Expedition 340, provides a long (>1 Ma) and detailed record of eruptive and mass‐wasting events (>130 discrete ...events). This record can be used to explore the temporal evolution in volcanic activity and landslides at an arc volcano. Analysis of tephra fall and volcaniclastic turbidite deposits in the drill cores reveals three heightened periods of volcanic activity on the island of Montserrat (∼930 to ∼900 ka, ∼810 to ∼760 ka, and ∼190 to ∼120 ka) that coincide with periods of increased volcano instability and mass‐wasting. The youngest of these periods marks the peak in activity at the Soufrière Hills volcano. The largest flank collapse of this volcano (∼130 ka) occurred toward the end of this period, and two younger landslides also occurred during a period of relatively elevated volcanism. These three landslides represent the only large (>0.3 km3) flank collapses of the Soufrière Hills edifice, and their timing also coincides with periods of rapid sea level rise (>5 m/ka). Available age data from other island arc volcanoes suggest a general correlation between the timing of large landslides and periods of rapid sea level rise, but this is not observed for volcanoes in intraplate ocean settings. We thus infer that rapid sea level rise may modulate the timing of collapse at island arc volcanoes, but not in larger ocean‐island settings.
Key Points
Heightened volcanic activity on Montserrat at 120–190, 760–810, and 900–930 ka
Large landslides coincide with rapid sea level rise at island arc volcanoes