The volcanic center of Santorini Island is the most active volcano of the southern Aegean volcanic arc. Α dense seismic array consisting of fourteen portable broadband seismological stations has been ...deployed in order to monitor and study the seismo-volcanic activity at the broader area of the Santorini volcanic center between March 2003 and September 2003. Additional recordings from a neighbouring larger scale temporary network (CYCNET) were also used for the relocation of more than 240 earthquakes recorded by both arrays. A double-difference relocation technique was used, in order to obtain optimal focal parameters for the best-constrained earthquakes.
The results indicate that the seismic activity of the Santorini volcanic center is strongly associated with the tectonic regime of the broader Southern Aegean Sea area as well as with the volcanic processes. The main cluster of the epicenters is located at the Coloumbo Reef, a submarine volcano of the volcanic system of Santorini Islands. A smaller cluster of events is located near the Anydros Islet, aligned in a NE–SW direction, running almost along the main tectonic feature of the area under study, the Santorini–Amorgos Fault Zone. In contrast, the main Santorini Island caldera is characterized by the almost complete absence of seismicity. This contrast is in very good agreement with recent volcanological and marine studies, with the Coloumbo volcanic center showing an intense high-temperature hydrothermal activity, in comparison to the corresponding low-level activity of the Santorini caldera.
The high-resolution hypocentral relocations present a clear view of the volcanic submarine structure at the Coloumbo Reef, showing that the main seismic activity is located within a very narrow vertical column, mainly at depths between 6 and 9 km. The focal mechanisms of the best-located events show that the cluster at the Coloumbo Reef is associated with the “Kameni–Coloumbo Fracture Zone”, which corresponds to the western termination of the major ENE–WSW Santorini–Amorgos Fault Zone. Stress–tensor inversion of the available fault plane solutions from Coloumbo Reef, as well as existing neotectonic fault information from NE Santorini (Coloumbo peninsula), suggests that the NE Santorini–Coloumbo faults belong to a single rupture system, with a ~
30° rotation of the local stress field with respect to the NNW–SSE regional extension field of the southern Aegean Sea. The observed change of the fault plane solutions shows that local conditions at the Coloumbo submarine volcano area control the observed faulting pattern.
In pancreaticobiliary maljunction (PBM), reflux of pancreatic juice and bile produces various pathological conditions in the biliary tract and pancreas. Clinical features according to the ...classification of PBM by confluence between the distal bile duct and the main pancreatic duct proposed in 2015 were evaluated in children.
Clinical features and complicating diseases according to the PBM classification were evaluated in 168 adult PBM patients. Patency of Santorini duct and associated biliary carcinomas were evaluated in 123 patients.
Similar to children, there were significant differences in age (P < 0.01) and type of common bile duct (P < 0.01) between the groups of the classification. Unlike in children, there was no significant difference in the incidence of abdominal pain and hyperamylasemia. There were 87 associated biliary carcinomas (79 gallbladder carcinomas and eight cholangiocarcinomas). PBM patients with a cudgel-type Santorini duct, which is greater than 2 mm in diameter, did not develop biliary carcinomas, compared to 61.1% of those with other types of Santorini duct (P < 0.01).
Clinical features according to the PBM classification in adults were different from those in children. Although biliary carcinomas were frequently seen in adult PBM patients, none of those with a cudgel-type Santorini duct developed biliary carcinoma.
Arc volcanoes are underlain by complex systems of molten‐rock reservoirs ranging from melt‐poor mush zones to melt‐rich magma chambers. Petrological and satellite data indicate that eruptible magma ...chambers form in the topmost few kilometres of the crust. However, very few chambers have ever been definitively located, suggesting that most are too short‐lived or too small to be imaged, which has direct implications for hazard assessment and modeling of magma differentiation. Here we use a high‐resolution technology based on inverting full seismic waveforms to image a small, high‐melt‐fraction magma chamber that was not detected with standard seismic tomography. The melt reservoir extends from ∼2 to at least 4 km below sea level (b.s.l.) at Kolumbo—a submarine volcano near Santorini, Greece. The chamber coincides with the termination point of the recent earthquake swarms and may be a missing link between a deeper melt reservoir and the high‐temperature hydrothermal system venting at the crater floor. The chamber poses a serious hazard as it could produce a highly explosive, tsunamigenic eruption in the near future. Our results suggest that similar reservoirs (relatively small but high‐melt‐fraction) may have gone undetected at other active volcanoes, challenging the existing eruption forecasts and reactive‐flow models of magma differentiation.
Plain Language Summary
Arc volcanoes, which mark the curved boundaries between converging tectonic plates, host the most explosive events on Earth. The associated hazard depends on how much mobile magma is currently present shallow beneath a volcano. Standard tomographic methods used so far have relatively low resolution and give a blurred picture of only the largest molten‐rock bodies. In particular, they struggle to distinguish between mobile magma and melt spread between tightly packed mineral grains. This study, a first in volcanology, combines a next‐generation tomographic method with extraordinarily dense seafloor recordings of controlled marine sound sources. This state‐of‐the‐art experiment at Kolumbo volcano, offshore of Santorini allowed us to detect a body of mobile magma which has been growing at an average rate of 4 × 106 m3 per year since the last eruption in 1650 CE. This rate is large enough to counteract the effect of cooling and crystallization. Our results show that Kolumbo poses a serious threat and deserves a real‐time monitoring facility. Despite the excellent data coverage, the small magma body was missed by standard tomography. This suggests that applying next‐generation imaging methods to already‐well‐studied volcanoes may lead to similar discoveries. We envision that small‐volume, high‐melt‐fraction reservoirs may be more widespread than previously thought.
Key Points
A shallow, very strong negative Vp anomaly imaged under the explosive, submarine Kolumbo volcano, Greece, using full‐waveform inversion
The high‐fidelity image and petrologic data indicate the anomaly is a small (∼0.6‐km wide, ∼2‐km deep), magma chamber with ∼42% of melt
The chamber was missed by travel‐time tomography indicating similar reservoirs may have gone undetected at other volcanoes
Santorini is the most active volcano of the Southern Aegean Volcanic Arc, with volcanic material that has accumulated during the last 600+ Kyrs on top of the pre-volcanic Santorini island. The ...geometry of the pre-volcanic basement not only provides constraints on the volcanic history, but is also fundamental for seismic hazard assessment, as it partly controls site-effects on strong ground motions. We investigate the geometry of the metamorphic (Cycladic) basement of Santorini using information from both passive (noise) and active source data. We performed a large number (∼280) of single-station ambient noise measurements and collected additional data from 42 previous studies to compute the Horizontal-to-Vertical Spectral Ratio (HVSR) curves. The majority of the recovered HVSR curves show prominent spectral ratio peaks (locally larger than 7–8) indicating a large impedance contrast between pyroclastic volcanic strata and bedrock metamorphic formations. We also performed active (MASW) and passive (Noise Array) surface wave measurements at 4 selected sites to determine local dispersion curves. To find the 1D Vs structure from the surface down to the metamorphic bedrock we performed a joint Monte Carlo inversion of the dispersion curves and the fundamental frequency (f0HVSR) values from the HVSR curves. These models, as well as deep borehole information and HVSR data along the caldera cliff, where the thickness and stratigraphy of the volcanic formations above the bedrock can be directly observed, allowed us to determine a linear relationship between the thickness of the pyroclastic formations and their fundamental resonance period for central-southern Santorini. The resulting geometry of pre-volcanic island shows deep basins around the pre-Alpine bedrock outcrop (Profitis Ilias) and generally agrees with recent results from larger-scale tomographic and gravity data but with some significant local differences. The basins are 200+ meters thick in the Kamari and Perissa area (southeastern Santorini) and >400 m thick in the central (Fira-Imerovigli) area, while there is a thick (∼70-130 m) pyroclastic layer on top of a metamorphic basement “plateau” in the Megalochori area.
•Prominent HVSR peaks due to the high impedance contrast of pyroclastics and bedrock.•Several VS profiles derived from joint inversion of noise arrays and active seismics.•Soil resonant period correlates with bedrock depth from Vs profiles and geology.•The observed correlation allowed to map the pre-volcanic metamorphic basement geometry.•The metamorphic bedrock model provides control on the volcanic history of Santorini.
Understanding the complex dynamics of volcanic systems demands a multidimensional approach that combines geophysics, geology, and geodetics. In this study, we examine observed spatiotemporal gravity ...changes within the Santorini volcanic complex from 1975 to 2014. The historical data indicates that gravity has been increasing continuously since at least 1966 until our latest measurements in 2014, albeit with a decreasing rate of increase over time. Utilizing gravity inversion of various gravity datasets and evidence from other studies, we explore different scenarios to shed light on the underlying processes. Our preferred interpretation involves both a magmatic episode and continuous evolution of the shallow structure. We find that the 2011-12 unrest period resulted from the intrusion of ∼3.3x1011 kg of basaltic magma at 3 km depth near the previously identified Mogi source. We attribute the continuous gravity increase beneath Nea Kameni to a density increases at about 1350 m depth. We infer these are a result of hydrothermal fluctuations, degassing, and/or vesicle collapse within the stored magma.
1mGal = 10-5 m/s2 (SI)
Pre-eruptive processes and their timescales are critical information for risk management at explosive volcanoes, and Santorini caldera (Greece) provides an excellent context in which to approach this ...subject. We ask two questions. First, are pre-eruptive processes the same for small and big eruptions? To investigate, we performed a multi-mineral diffusion timescale study of a small explosive eruption of Kameni Volcano and compared the results with those published for larger caldera-forming eruptions at Santorini. The Kameni dacite resembles products of larger eruptions in being crystal-poor, containing plagioclase with antecrystic cores and autocrystic rims, bearing orthopyroxene with sector zoning and phantom skeletal morphologies, and showing evidence for mixing of different silicic magmas prior to eruption. Diffusion timescales from Mg-Fe profiles in orthopyroxene and clinopyroxene phenocrysts are <1–23 years, and Mg diffusion modelling in plagioclase gives <10 years. Our physical model for the Kameni eruption is similar to those proposed for larger eruptions, where silicic melt produced in gabbroic to dioritic lower to middle crustal mush bodies is transferred (along with entrained mafic magma) to an upper crustal reservoir. Crystals grow in the hydrous silicic melts due to decompression, cooling, and magma mixing during ascent and injection into upper crust. We propose that large eruptions are preceded by similar processes as small ones, but on a larger scale. Our second question: do diffusion timescales relate to eruptive volume or position in a caldera cycle? For this, we obtained orthopyroxene Mg-Fe diffusion timescales for three additional eruptions, growing our orthopyroxene timescale database to seven eruptions of different sizes and cycle timings. No clear relationship exists between diffusion timescale and volume; however, timescales are systematically shorter (<0.01–10 years) early in a cycle and longer (1–5,000 years) late in a cycle. Thermal maturation and H 2 O-flushing of the crustal magma reservoir through the caldera cycle could explain this, as the reservoir would change from a rigid to more mushy state as the cycle progresses. This would change the mechanical response to melt input and allow accumulation of progressively larger melt layers in the upper crust, resulting in increasing crystal residence times.
The Santorini-Amorgos zone is located in the central part of the Hellenic volcanic arc and is hosting eight large faults as well as Kolumbo and Santorini volcanic centers. The largest earthquake ...(Mw ~ 7.1) in the southern Aegean during the 20th century also occurred in this area on 9 July 1956. A total of 1868 crustal events were recorded by temporary networks during September 2002 to July 2004 and October 2005 to March 2007, and also by the permanent network from 2011 to 2019. We relocated 1455 of these events by using HypoDD and revealed clusters of earthquakes beneath Kolumbo, Anydros graben, and Santorini-Amorgos ridge. Only the faults in the SW of Anydros, SE of Ios, and along the south coast of Amorgos were delineated by the relocated events. Nearly vertical clusters were observed beneath the island of Anydros, south of Amorgos, and in NE end of Amorgos fault, indicating possible pathways of upward migrating fluids. The seismogenic layer thickness calculated based on the depth distribution of the relocated events was 12.5 km. We combined this thickness with geometrical properties of the faults to calculate the expected moment magnitude of future earthquakes, resulting in a range of 6.3–7.2. In an effort to map the distribution of fluids, the Vp/Vs ratio distribution was estimated by utilizing the event-station travel time data along with crack density, fluid saturation, and Poisson's ratio. The petrophysical parameters observed in the northern part of the Santorini caldera suggest the existence of melt, while those observed in Anydros and in the NE of Amorgos fault support the suggestion of upward migrating fluids in these areas.
•A total of 1455 crustal events in Santorini-Amorgos zone have been relocated.•The seismogenic layer along the zone is found to be 12.5 km thick.•Expected moment magnitude of future earthquakes is in the range of 6.3 to 7.2.•High Vp/Vs ratios in northern part of Santorini caldera indicate the presence of melt.•Upward migrating fluids exist at areas with vertical earthquake clusters.
The destructive tsunami on 22 December 2018 due to the flank collapse of the Anak Krakatau volcano was a bitter reminder of large tsunami risks and of the shortcomings of the existing tsunami warning ...systems for atypical sources (tsunamis generated by non-seismic and complex sources). In the Mediterranean, several tsunamis were generated by landslides associated with volcanic systems in the past.The volcanic unrest experienced in 2011–2012 on the Santorini volcanic island in the Southern Aegean Sea pointed out the need to identify and quantify tsunami hazard and risk due to possible flank instability which may be triggered as a result of volcanic unrest or nearby seismotectonic activities. Inspired from this need, in this study we examined three possible landslide scenarios in Santorini Island with tsunamigenic potential. The results show that the scenarios considered in our study are able to generate significant local tsunamis impacting Santorini and the nearby islands, as well as producing significant impact along the coasts of the Southern Aegean Sea. While maximum tsunami amplitudes/arrival time ranges are 1.2 m/30-90 min for locations in the Greek-Turkish coasts in the far field, they are in the order of ≈60 m/1-2 min for some locations at the Santorini Island. The extreme tsunami amplitudes and short arrival times for locations inside the Santorini Island is a major challenge in terms of tsunami hazard warning and mitigation. As an effort to address this challenge, a discussion on the requirements for local tsunami warning system addressing atypical sources in the context of multi-hazard disaster risk reduction is also provided.
The explosive eruption at Santorini in the Aegean Sea during the second millennium BCE was the largest Holocene volcanic upheaval in the Eastern Mediterranean region. The eruption was disastrous for ...the Minoan settlements at Santorini, but the effect on human society in the neighbouring islands and regions is still clouded in uncertainty. Tsunami generation was suggested, but comparatively little evidence was found. The lack of firm tsunami traces is particularly puzzling in Crete with its coastal settlements of the Late Minoan IA period, during which the Santorini eruption occurred. Here, we report the discovery of extensive geoarchaeological tsunami deposits at Palaikastro in north-eastern Crete. These deposits are characterized by a mixture of geological materials, including volcanic Santorini ash, and archaeological settlement debris. Various tsunami signatures were identified: (1) erosional contact with the underlying strata, (2) volcanic ash intraclasts in the lower part of the deposit, (3) reworked building stone material in the lower part of the deposit, (4) individual marine shells, (5) marine micro-fauna, (6) imbrication of rounded beach pebbles, settlement debris, ceramic sherds and even bones, (7) multi-modal chaotic composition. Late Minoan human settlement activities at Palaikastro provided architectural and stratigraphic frameworks in space and time that recorded and preserved tsunami evidence as geoarchaeological deposits. Such stratigraphic resolution and preservation may not occur in the natural landscape. Volcanic ash transported by wind from Santorini south-east to Crete preceded the tsunami. Geological, archaeological and radiocarbon dating criteria all converge, indicating that the tsunami deposits are coeval with the Minoan Santorini eruption. Field evidence suggests that tsunami waves at Palaikastro were at least 9
m high. Inverse tsunami modeling was attempted, based on these newly discovered tsunamigenic deposits. The initial wave in the generation region at Santorini that best fits the stratigraphic data is a wave with +35 to −15
m initial amplitude and a crest length of about 15
km.
The South Aegean Volcanic Arc overlies a slowly subducting, cool slab of oceanic-to-transitional crust, and hosts the hazardous Christiana–Santorini–Kolumbo volcanic field. In order to investigate ...the primitive melts feeding the volcanic field, we present major and trace element analyses of 130 olivine-hosted melt inclusions from Santorini, integrated with previously published H
2
O and CO
2
data. Following post-entrapment corrections, we identify four endmember primitive melt types preserved in Fo ≥ 80 olivines, ranging from low-K island-arc basalts with La/Yb ~ 1.5 and 1.5–3.0 wt% H
2
O to andesites with La/Yb ~ 6–10 and 3.0–3.5 wt% H
2
O. They are consistent with melting at 1.3 to 2.3 GPa and 1350–1440 °C of variably depleted peridotitic mantle fluxed by slab-derived melts and fluids. The chemical signatures of sediment melts dominate, while those of fluids derived from the ocean crust are low compared to global datasets. This is consistent with thick sediment accumulations observed in the Hellenic trench, and with low calculated fluid fluxes from the downgoing slab. The low H
2
O contents estimated for the primary melts (0.8–1.8 wt%) may imply a component of decompression melting beneath the arc. Coupled with a well-constrained chronostratigraphic context, the melt inclusion archive provides a time series of mantle-derived input into the silicic crustal magmatic system over the last 530 ka. Primitive melts with La/Yb ≤ 5 have been erupted encased in olivines over the last 530 ky, without any evident time variation. Melt inclusions with La/Yb > 5 have, on the other hand, been restricted to two periods: (1) prior to the onset of major explosive volcanism at ~ 360 ka, and (2) the products of the 3.6 ka Late-Bronze-Age eruption and the 22-to-3.6 ka inter-Plinian period immediately preceding it. The observations may be explained by time-varying differential extraction of melts from deep storage zones in the mantle or lower crust, related to lithospheric rifting and caldera collapse events. Temporal variations in the supplies of slab-derived melts and fluids may also play a role.