Santorini volcano in the South Aegean Volcanic Arc has a detailed history of ongoing volcanic and seismic activity, making it a prime location for studying magma storage and transport at arc volcanoes. The shallow magmatic system (<5 km depth) is well constrained by geophysical studies, but the deeper crustal structure is not. Located 15 km NE of Santorini, the Kolumbo seamount is also an active edifice, with consistently more seismicity and hydrothermal venting than Santorini. Geochemical studies indicate that Santorini and Kolumbo are fed by separate mantle and crustal magma sources, but prior seismic studies suggest otherwise (Dimitriadis et al, 2010; McVey et al, 2020). This study addresses the nature of lower-crustal magma structure beneath arc volcanoes and whether evolved volcanoes and nearby vents are connected through their plumbing. Tomographic inversion of P-wave Moho reflection (PmP) and turning P-wave (Pg) traveltimes is used to create 3-D models of Moho depth and P-wave velocity (Vp) down to depths of ~25 km. The PROTEUS experiment provides an exceptionally dense and large aperture traveltime dataset from an amphibious array of ~150 seismometers and ~14,000 active marine sources. The data are ~33,000 manually picked PmP arrivals and ~256,000 Pg arrivals from existing studies. Results show a low Vp anomaly extending from the Moho to the surface. This anomaly starts at the base of the crust under the NW Santorini caldera and extends up to the east. It is most pronounced at 10-15 km depth, where it is offset from both Santorini and Kolumbo. Limited resolution prevents imaging of a connection between this mid-crustal anomaly and the known shallow magma storage region under the Santorini caldera. A high-velocity core beneath Santorini is not found, a feature interpreted at other volcanoes as a cooled intrusive complex. Because no additional low Vp anomalies are found in the lower crust, we infer that a common mantle source and mid-crustal plumbing system is actively feeding both Santorini and Kolumbo. The spatial offset and elongated nature of magma storage implies a complex relationship between evolving magmatic structures and tectonics.
more »
« less
Melt distribution and Vp/Vs for the Mid Crust at the Santorini-Kolumbo Volcanic System
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
- Award ID(s):
- 2023338
- PAR ID:
- 10567273
- Publisher / Repository:
- 2023 Fall Meeting, American Geophysical Union
- Date Published:
- Format(s):
- Medium: X
- Location:
- https://agu.confex.com/agu/fm23/meetingapp.cgi/Paper/1316502
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract Despite multidisciplinary evidence for crustal magma accumulation below Santorini volcano, Greece, the structure and melt content of the shallow magmatic system remain poorly constrained. We use three-dimensional (3-D) velocity models from tomographic inversions of active-source seismic P-wave travel times to identify a pronounced low-velocity anomaly (–21%) from 2.8 km to 5 km depth localized below the northern caldera basin. This anomaly is consistent with depth estimates of pre-eruptive storage and a recent inflation episode, supporting the interpretation of a shallow magma body that causes seismic attenuation and ray bending. A suite of synthetic tests shows that the geometry is well recovered while a range of melt contents (4%–13% to fully molten) are allowable. A thin mush region (2%–7% to 3%–10% melt) extends from the main magma body toward the northeast, observed as low velocities confined by tectono-magmatic lineaments. This anomaly terminates northwest of Kolumbo; little to no melt underlies the seamount from 3 to 5 km depth. These structural constraints suggest that crustal extension and edifice loads control the geometry of magma accumulation and emphasize that the shallow crust remains conducive to melt storage shortly after a caldera-forming eruption.more » « less
-
Abstract With substantial postglacial rhyolite eruptions and ongoing rapid uplift, the Laguna del Maule volcanic field in the southern Andes provides an exceptional opportunity to study the dynamics of an active silicic magmatic system. Using 4,093Parrivals from 137 distant earthquakes recorded by 44 local stations over∼2.25 years, we conduct teleseismic tomography to image the crustal structure down to 40 km below the volcanic field. A prominent low‐velocity body at depths between∼0 and 12 km below sea level (b.s.l.), characterized by a volume of∼500 km3and a peak anomaly of−400 m/s (∼9%), overlaps the location of the upper‐crustal magma reservoir detected in recent gravity and surface wave tomography studies. Its estimated averagePwave velocity of∼4.6 km/s corresponds to an average melt fraction of about 14% and a melt volume of∼70 km3. Petrologic observations are also consistent with generation and storage of rhyolitic melts at depths corresponding to the anomalous zone. Moreover, the tomographic results support a lower crust zone of MASH (melting, assimilation, storage, and homogenization) from a depth of∼25 km to the base of the model, which likely reflects a deep crustal source of magma that contributes to and incubates the shallow silicic reservoir.more » « less
-
null (Ed.)Caldera footprints and erupted magma volumes provide a unique constraint on vertical dimensions of upper crustal magma reservoirs that feed explosive silicic eruptions. Here we define a Vertical Separation (VS) ratio in which we compare the geometric vertical extent with the range of depths indicated petrologically by melt inclusion water and CO2 saturation pressures for fifteen caldera-forming eruptions spanning ∼10^0 km3 to ∼10^3 km3 in volume. We supplement melt inclusion saturation pressures with rhyolite-MELTS barometry and plagioclase-melt hygrometry to generate a petrologic image of magma reservoir architecture. We find that pre-eruptive upper crustal magma reservoirs range from contiguous bodies (where petrologic and geometric estimates match closely) to vertically dispersed structures. Vertically dispersed pre-eruptive reservoirs are more common among intermediate-volume eruptions than among the smallest and largest caldera-forming eruptions. We infer that the architecture of magma reservoirs tracks the thermomechanical evolution of large volcanic systems.more » « less
-
The principal aim at Site U1589 (proposed Site CSK-01A) was to reconstruct the evolution of the Anhydros Basin, including its history of subsidence, as well as to document the presence of volcanic event layers in the basin sediments and draw conclusions regarding the links between volcanism and crustal tectonics. The site is located about 10 km southwest of Amorgos Island at 484 meters below sea level (mbsl) (Figure F1). The drill site targeted the volcano-sedimentary fill of the Anhydros Basin. We received permission from the International Ocean Discovery Program (IODP) Environmental Protection and Safety Panel to touch the Alpine basement using an advanced piston corer/extended core barrel/rotary core barrel (APC/XCB/RCB) drilling strategy. The site involved three holes (U1589A–U1589C) and terminated in basement limestone at 612.4 meters below seafloor (mbsf) (all depths below seafloor are given using the core depth below seafloor, Method A [CSF-A] scale, except in Operations, where the drilling depth below seafloor [DSF] scale is used). Core recovery was good in Holes U1589A (78%) and U1589B (87%) and poor in Hole U1589C (24%). Site U1589 was chosen to sample the eruptive histories of both Santorini and the Kolumbo chain and was expected to yield volcaniclastics from many Kolumbo eruptions and the major Santorini eruptions. Many deposits from smaller Santorini eruptions were not expected at this distance from the volcano, in part due to flow blocking by the Kolumbo volcanic chain. Santorini has been active since 0.65 Ma, with many large explosive eruptions since about 0.36 Ma (Druitt et al., 2016). Kolumbo Volcano was known from seismic profiles to have had at least five eruptions (Hübscher et al., 2015), the last of which was in 1650 Common Era (CE) and killed 70 people on Santorini (Fuller et al., 2018). Seismic profiles provided constraints on the relative ages of the Kolumbo cones (Preine et al., 2022c) but not on the absolute ages. The site offered a near-continuous time series of volcanism in the area since rift inception. The site was also chosen to develop a core-log-seismic integration stratigraphy and compare it with the recently published seismic stratigraphy for the basin (Preine et al., 2022a) and the paleotectonic reconstruction of the region (Nomikou et al., 2016, 2018) to promote a holistic view of the Anhydros Basin evolution (Figure F2). The site transects all six seismic packages of the Anhydros rift basin, as well as the onlap surfaces between them (Nomikou et al., 2016, 2018; Preine et al., 2022a). The anticipated lithologies were undisturbed hemipelagic muds, volcaniclastics, turbidites, and finally continental basement rocks. A gravity core recovered 7 km to the east indicated that the uppermost sediments on site would consist of hemipelagic muds and volcaniclastic layers, as well as sapropels (Kutterolf et al., 2021). The Anhydros Basin is crossed by many seismic profiles obtained in campaigns between 2006 and 2019, many of them multichannel (Hübscher et al., 2015; Nomikou et al., 2016, 2018), and its southwestern part is included within the area of the 2015 PROTEUS seismic tomography experiment, during which subbottom profiling, gravity, and magnetic data were also recorded (Hooft et al., 2017). The basin bathymetry had been studied in several marine campaigns, and fault distributions and throws had been mapped (Nomikou et al., 2016; Hooft et al., 2017). Previously published analyses of the seismic data suggested the following possible interpretations (from the bottom up; Preine et al., 2022b, 2022c): Units U1 and U2: sediment packages predating Santorini and Kolumbo volcanism; Unit U3: sediments and the products of the early Kolumbo volcanism and some of the Kolumbo cones; Unit U4: sediments associated with a major rift pulse; and Units U5 and U6: sediments and the products of Santorini activity, some of the Kolumbo cones, and the later eruptions of Kolumbo including the 1650 CE eruption. Units U3–U6 were believed to be of Pleistocene age, and Units U1 and U2 were believed to be of possible Pliocene age. The site enabled us to test these interpretations by using the cores to reconstruct a near-complete volcanic stratigraphy consistent with both onshore and offshore constraints and pinned by chronological markers from biostratigraphy, magnetostratigraphy, and sapropel records. Benthic foraminifera from fine-grained sediments provided estimates of paleowater depths and, via integration with seismic profiles and chronologic data, of time-integrated basin subsidence rates. Coring at Site U1589 in the Anhydros Basin addressed scientific Objectives 1–4 and 6 of the Expedition 398 Scientific Prospectus (Druitt et al., 2022). It was complemented by Site U1592 in the Anafi Basin because each basin taps a different sediment distributary branch of the Christiana-Santorini-Kolumbo volcanic system.more » « less