The Santorini arc volcano in the Hellenic subduction zone has a history of caldera-forming Plinian eruptions, most recently in the Late Bronze Age 3.4 kya, and it remains volcanically active. To inform volcanic hazard assessments, it is crucial to understand where melt is distributed. The PROTEUS experiment in 2015 recorded >14,000 controlled marine sound sources on 165 land and seafloor seismic stations. Tomographic inversion of this data revealed low P-wave velocities in the upper 4 kilometers beneath the caldera and nearby Kolumbo seamount interpreted as the magma system (McVey et al., 2020; Chrapkewiecz et al., 2022). However, structure of the magma system was only determined in the upper (<4-6km) crust and melt content is only weakly constrained. In this study we improve constraints on the deeper magma system and subsurface melt content with a tomographic P and S wave velocity structure. To do so, we add to the inverse problem arrival times from ~1500 local earthquakes with magnitudes from 0.5 to 3.0 that occurred between 5 and 20 km depth. The events were recorded on 142 3-component ocean bottom and island seismic stations that span the seafloor ~60 km west and east of the island and the nearby islands. Results beneath Santorini and Kolumbo suggest that the upper crustal magma reservoirs extend deeper than previously found, and we image a high Vp layer (~5-8 km) under the magma reservoir at Kolumbo. We identify this layer as strong, cooled, intruded magma and correlate it to the location of earthquakes, within which, swarms of rapidly upward propagating seismicity support prior inferences of melt conduits traversing a rheologically strong layer (Schmid et al, 2022). We give values for melt content of the upper crustal reservoirs using a scaled Vp/Vs model. Since the number of arrivals, apriori assigned uncertainty, and differences in ray geometry can result in P and S waves with different resolving power, we use measured resolution to scale the Vs perturbations and create a more realistic Vp/Vs model. The addition of earthquake arrivals allows us to map the magma reservoirs beneath the Santorini-Kolumbo magma system to 8 km depth and identify regions of elevated melt content.
more »
« less
Deep long-period earthquakes generated by second boiling beneath Mauna Kea volcano
Deep long-period earthquakes (DLPs) are an enigmatic type of volcanic seismicity that sometimes precedes eruptions but mostly occurs at quiescent volcanoes. These earthquakes are depleted in high-frequency content and typically occur near the base of the crust. We observed a near-periodic, long-lived sequence of more than one million DLPs in the past 19 years beneath the dormant postshield Mauna Kea volcano in Hawaiʻi. We argue that this DLP sequence was caused by repeated pressurization of volatiles exsolved through crystallization of cooling magma stalled beneath the crust. This “second boiling” of magma is a well-known process but has not previously been linked to DLP activity. Our observations suggest that, rather than portending eruptions, global DLP activity may more commonly be indicative of stagnant, cooling magma.
more »
« less
- Award ID(s):
- 1848302
- PAR ID:
- 10183747
- Date Published:
- Journal Name:
- Science
- Volume:
- 368
- Issue:
- 6492
- ISSN:
- 0036-8075
- Page Range / eLocation ID:
- 775 to 779
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The architecture of magma plumbing systems plays a fundamental role in volcano eruption and evolution. However, the precise configuration of crustal magma reservoirs and conduits responsible for supplying eruptions are difficult to explore across most active volcanic systems. Consequently, our understanding of their correlation with eruption dynamics is limited. Axial Seamount is an active submarine volcano located along the Juan de Fuca Ridge, with known eruptions in 1998, 2011, and 2015. Here we present high-resolution images of P-wave velocity, attenuation, and estimates of temperature and partial melt beneath the summit of Axial Seamount, derived from multi-parameter full waveform inversion of a 2D multi-channel seismic line. Multiple magma reservoirs, including a newly discovered western magma reservoir, are identified in the upper crust, with the maximum melt fraction of ~15–32% in the upper main magma reservoir (MMR) and lower fractions of 10% to 26% in other satellite reservoirs. In addition, a feeding conduit below the MMR with a melt fraction of ~4–11% and a low-velocity throat beneath the eastern caldera wall connecting the MMR roof with eruptive fissures are imaged. These findings delineate an asymmetric shallow plumbing system beneath Axial Seamount, providing insights into the magma pathways that fed recent eruptions.more » « less
-
Abstract Classical mechanisms of volcanic eruptions mostly involve pressure buildup and magma ascent towards the surface1. Such processes produce geophysical and geochemical signals that may be detected and interpreted as eruption precursors1–3. On 22 May 2021, Mount Nyiragongo (Democratic Republic of the Congo), an open-vent volcano with a persistent lava lake perched within its summit crater, shook up this interpretation by producing an approximately six-hour-long flank eruption without apparent precursors, followed—rather than preceded—by lateral magma motion into the crust. Here we show that this reversed sequence was most likely initiated by a rupture of the edifice, producing deadly lava flows and triggering a voluminous 25-km-long dyke intrusion. The dyke propagated southwards at very shallow depth (less than 500 m) underneath the cities of Goma (Democratic Republic of the Congo) and Gisenyi (Rwanda), as well as Lake Kivu. This volcanic crisis raises new questions about the mechanisms controlling such eruptions and the possibility of facing substantially more hazardous events, such as effusions within densely urbanized areas, phreato-magmatism or a limnic eruption from the gas-rich Lake Kivu. It also more generally highlights the challenges faced with open-vent volcanoes for monitoring, early detection and risk management when a significant volume of magma is stored close to the surface.more » « less
-
Abstract Seismicity during explosive volcanic eruptions remains challenging to observe through the eruptive noise, leaving first‐order questions unanswered. How do earthquake rates change as eruptions progress, and what is their relationship to the opening and closing of the eruptive vent? To address these questions for the Okmok Volcano 2008 explosive eruption, Volcano Explosivity Index 4, we utilized modern detection methods to enhance the existing earthquake catalog. Our enhanced catalog detected significantly more earthquakes than traditional methods. We located, relocated, determined magnitudes and classified all events within this catalog. Our analysis reveals distinct behaviors for long‐period (LP) and volcano‐tectonic (VT) earthquakes, providing insights into the opening and closing cycle. LP earthquakes occur as bursts beneath the eruptive vent and do not coincide in time with the plumes, indicating their relationship to an eruptive process that occurs at a high pressurization state, that is, partially closed conduit. In contrast, VT earthquakes maintain a steadier rate over a broader region, do not track the caldera deflation and have a largerb‐value during the eruption than before or after. The closing sequence is marked by a burst of LPs followed by small VTs south of the volcano. The opening sequence differs as only VTs extend to depth and migrate within minutes of the eruption onset. Our high‐resolution catalog offers valuable insights, demonstrating that volcanic conduits can transition between partially closed (clogged) and open (cracked) states during an eruption. Utilizing modern earthquake processing techniques enables clearer understanding of eruptions and holds promise for studying other volcanic events.more » « less
-
null (Ed.)Abstract Increased resolution of data constraining topography and crustal structures provides new quantitative ways to assess province-scale surface-subsurface connections beneath volcanoes. We used a database of mapped vents to extract edifices with known epoch ages from digital elevation models (DEMs) in the Cascades arc (western North America), deriving volumes that likely represent ∼50% of total Quaternary eruptive output. Edifice volumes and spatial vent density correlate with diverse geophysical data that fingerprint magmatic influence in the upper crust. Variations in subsurface structures consistent with volcanism are common beneath Quaternary vents throughout the arc, but they are more strongly associated with younger vents. Geophysical magmatic signatures increase in the central and southern Cascade Range (Cascades), where eruptive output is largest and vents are closely spaced. Vents and correlated crustal structures, as well as temporal transitions in the degree of spatially localized versus distributed eruptions, define centers with lateral extents of ∼100 km throughout the arc, suggesting a time-evolving spatial focusing of magma ascent.more » « less