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   
                    
                            
                            Heralds of Future Volcanism: Swarms of Microseismicity Beneath the Submarine Kolumbo Volcano Indicate Opening of Near‐Vertical Fractures Exploited by Ascending Melts
                        
                    
    
            Abstract The Kolumbo submarine volcano in the southern Aegean (Greece) is associated with repeated seismic unrest since at least two decades and the causes of this unrest are poorly understood. We present a ten‐month long microseismicity data set for the period 2006–2007. The majority of earthquakes cluster in a cone‐shaped portion of the crust below Kolumbo. The tip of this cone coincides with a low Vp‐anomaly at 2–4 km depth, which is interpreted as a crustal melt reservoir. Our data set includes several earthquake swarms, of which we analyze the four with the highest events numbers in detail. Together the swarms form a zone of fracturing elongated in the SW‐NE direction, parallel to major regional faults. All four swarms show a general upward migration of hypocenters and the cracking front propagates unusually fast, compared to swarms in other volcanic areas. We conclude that the swarm seismicity is most likely triggered by a combination of pore‐pressure perturbations and the re‐distribution of elastic stresses. Fluid pressure perturbations are induced likely by obstructions in the melt conduits in a rheologically strong layer between 6 and 9 km depth. We conclude that the zone of fractures below Kolumbo is exploited by melts ascending from the mantle and filling the crustal melt reservoir. Together with the recurring seismic unrest, our study suggests that a future eruption is probable and monitoring of the Kolumbo volcanic system is highly advisable. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2023338
- PAR ID:
- 10370510
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 23
- Issue:
- 7
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Deep‐crustal magma plumbing at arc volcanoes controls the volume, frequency, and composition of magma being transported to and stored in the upper crust. However, the mid‐to‐lower crust remains a challenging region to image. We explore the mid‐to‐lower crustal velocity structure beneath the Christiana‐Santorini‐Kolumbo Volcanic Field (CSKVF) to better understand how an established stratovolcano and flanking volcano (Santorini and Kolumbo) are fed through the mid‐to‐lower crust. We use active‐source seismic data to obtain a P‐wave velocity model of the crust below the CSKVF. We invert direct and reflected P phases to cover the entire depth extent of the crust and solve for the Moho interface depth. Our model requires a curved Moho interface representative of crustal thickening via underplating. Results show a highVpanomaly in the lower crust under Santorini and a mid‐crustal lowVpanomaly offset from both Santorini and Kolumbo. We find that accumulation of magma is located under the local extensional basin in the upper mid‐crust (<10 km) but is offset at deeper depths. We find evidence for melt storage at 11–13 km depth feeding volcanism at the Kolumbo volcanic chain. This melt is also a plausible source for the 2025 seismic swarm and dike intrusion. Resolution is limited in the mid‐crust below the Santorini caldera, leaving Santorini's mid‐crustal magma plumbing unconstrained. We think it likely that Santorini and Kolumbo have entirely separate crustal plumbing systems and mantle sources, but allow the possibility of a connection in the mid or lower crust.more » « less
- 
            Santorini volcano has a history of caldera‐forming eruptions, most recently in the Late Bronze Age, at 3.4 kya, and remains volcanically active. The Kolumbo submarine volcano, located 7 km to the northeast of Santorini, erupted in 1650 AD in a deadly phreatomagmatic eruption. Ongoing seismic activity and active hydrothermal venting at Kolumbo indicate this volcano is a significant hazard to the Santorini region. The magma source for Santorini and the Kolumbo edifice are considered separate in the shallow crust, though their deeper magma distribution is not yet constrained. In this study, we improve constraints on the mid‐crustal magma system of Santorini caldera and the nearby Kolumbo volcano using local earthquake tomography. We use 1515 P‐wave and 1435 S‐wave arrival times from 63 local earthquakes with magnitudes from 0.5 to 3.0 that occurred between 5 and 15 km depth together with an existing data set of active source Pg arrivals. The upper crustal magma system beneath Santorini is imaged to at least 6 km depth, and to 12 km depth beneath Kolumbo. We recover a high P‐wave velocity layer (∼6–8 km) under the Kolumbo magma reservoir that we infer is a rheologically strong seismogenic layer. We also recover a mid‐crustal magma body below 8 km depth located to the NE of Santorini and Kolumbo.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
- 
            Abstract The Laguna del Maule (LdM) volcanic field comprises the greatest concentration of postglacial rhyolite in the Andes and includes the products of ~40 km3of explosive and effusive eruptions. Recent observations at LdM by interferometric synthetic aperture radar and global navigation satellite system geodesy have revealed inflation at rates exceeding 20 cm/year since 2007, capturing an ongoing period of growth of a potentially large upper crustal magma reservoir. Moreover, magnetotelluric and gravity studies indicate the presence of fluids and/or partial melt in the upper crust near the center of inflation. Petrologic observations imply repeated, rapid extraction of rhyolitic melt from crystal mush stored at depths of 4–6 km during at least the past 26 ka. We utilize multiple types of surface‐wave observations to constrain the location and geometry of low‐velocity domains beneath LdM. We present a three‐dimensional shear‐wave velocity model that delineates a ~450‐km3shallow magma reservoir ~2 to 8 km below surface with an average melt fraction of ~5%. Interpretation of the seismic tomography in light of existing gravity, magnetotelluric, and geodetic observations supports this model and reveals variations in melt content and a deeper magma system feeding the shallow reservoir in greater detail than any of the geophysical methods alone. Geophysical imaging of the LdM magma system today is consistent with the petrologic inferences of the reservoir structure and growth during the past 20–60 kyr. Taken together with the ongoing unrest, a future rhyolite eruption of at least the scale of those common during the Holocene is a reasonable possibility.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
