Description: Santorini volcanic field has had 12 major (1–10 km3 or more of magma), and numerous minor, explosive eruptions over the last ~ 200 ka. Deposits from these eruptions (Thera Pyroclastic Formation) are well exposed in caldera-wall successions up to 200 m thick. Each of the major eruptions began with a pumice-fall phase, and most culminated with emplacement of pyroclastic flows. Pyroclastic flows of at least six eruptions deposited proximal lag deposits exposed widely in the caldera wall. The lag deposits include coarse-grained lithic breccias (andesitic to rhyodacitic eruptions) and spatter agglomerates (andesitic eruptions only). Facies associations between lithic breccia, spatter agglomerate, and ignimbrite from the same eruption can be very complex. For some eruptions, lag deposits provide the only evidence for pyroclastic flows, because most of the ignimbrite is buried on the lower flanks of Santorini or under the sea. At least eight eruptions tapped compositionally heterogeneous magma chambers, producing deposits with a range of zoning patterns and compositional gaps. Three eruptions display a silicic–silicic + mafic–silicic zoning not previously reported. Four eruptions vented large volumes of dacitic or rhyodacitic pumice, and may account for 90% or more of all silicic magma discharged from Santorini. The Thera Pyroclastic Formation and coeval lavas record two major mafic-to-silicic cycles of Santorini volcanism. Each cycle commenced with explosive eruptions of andesite or dacite, accompanied by construction of composite shields and stratocones, and culminated in a pair of major dacitic or rhyodacitic eruptions. Sequences of scoria and ash deposits occur between most of the twelve major members and record repeated stratocone or shield construction following a large explosive eruption. Volcanism at Santorini has focussed on a deep NE–SW basement fracture, which has acted as a pathway for magma ascent. At least four major explosive eruptions began at a vent complex on this fracture. Composite volcanoes constructed north of the fracture were dissected by at least three caldera-collapse events associated with the pyroclastic eruptions. Southern Santorini consists of pryoclastic ejecta draped over a pre-volcanic island and a ridge of early- to mid-Pleistocene volcanics. The southern half of the present-day caldera basin is a long-lived, essentially non-volcanic, depression, defined by topographic highs to the south and east, but deepened by subsidence associated with the main northern caldera complex, and is probably not a separate caldera.

Description: The 18,500 yr. b.p. Cape Riva (CR) eruption of Santorini vented several km3 or more of magma, generating four eruption units: a basal Plinian fall deposit (CR-A) and three pyroclastic flow deposits (CR-B to CR-D upwards). CR-B and CR-D are welded ignimbrites; CR-C consists predominantly of up to 25 m thick coarse, lithic-rich co-ignimbrite lag breccias resulting from a climactic phase of the eruption. The initial Plinian phase occurred from a localized vent in N Santorini, and subsequent column collapse resulted in emplacement of CR-B. Towards the end of CR-B, new conduits were activated and pyroclastic flows discharged from multiple vents to generate the lag breccias (CR-C). CR-D probably records a return to a localized vent as the eruption waned. The eruption sampled a zoned magma chamber containing rhyodacite overlying andesite, and leaks of these magmas were manifested as the Skaros-Therasia lavas preceding the CR eruption. Plinian and initial ignimbrite stages occurred while the magma chamber was overpressured; subsequent underpressuring, due to magma discharge, caused fracturing of the chamber roof, caldera collapse, and eruption of pyroclastic flows from multiple vents. Activation and widening of new conduits during collapse resulted in the rapid escalation of discharge rate favoring the formation of lag breccias by: (i) promoting erosion of lithic debris at the surface vent; and (ii) raising surface exit pressures, thereby resulting in a dramatic increase in the grain size of the ejecta.

Description: The Milos, Christiana-Santorini-Kolumbo (CSK) and Kos-Yali-Nisyros (KYN) volcanic complexes of the Aegean Volcanic Arc have repeatedly produced highly explosive eruptions from at least ∼360 ka into historic times and still show recent unrest. We present the marine tephra record from an array of 50, up to 7.4 m long, sediment cores along the arc collected in 2017 during RV Poseidon cruise POS513, which complements earlier work on distal to ultra-distal eastern Mediterranean sediment cores. A unique set of glass-shard trace element (LA-ICPMS) compositions complements our major element (EMP) data on 220 primary ash layers and 40 terrestrial samples to support geochemical fingerprinting for correlations with 19 known tephras from all three volcanic complexes and with the 39 ka Campanian Ignimbrite from the Campi Flegrei, Italy. The correlations include eleven eruptions from CSK (Kameni, Kolumbo 1650, Minoan, Cape Riva, Cape Tripiti, Upper Scoriae 1 and 2, Middle Pumice, Cape Thera, Lower Pumice, Cape Therma 3). We identify a previously unknown widespread tephra from a plinian eruption on Milos (Firiplaka Tephra). Near the KYN we correlate marine tephras with the Kos Plateau Tuff, the Yali 1 and Yali 2 tephras, and the Upper and Lower Pumice on Nisyros. Between these two major tephras, we found two tephras from Nisyros not yet observed on land. The four Nisyros tephras form a systematic trend toward more evolved magma compositions. In the companion paper we use the tephrostratigraphic framework established here to constrain new eruption ages and magnitudes as a contribution to volcanic hazard assessment.

Description: We use the tephrostratigraphic framework along the Aegean Volcanic Arc established in part 1 of this contribution to determine hemipelagic sedimentation rates, calculate new tephra ages, and constrain the minimum magnitudes of (sub)plinian eruptions of the last 200 kyrs. Hemipelagic sedimentation rates range from ∼0.5 cm/kyr up to ∼40 cm/kyr and vary laterally as well as over time. Interpolation between dated tephras yields an eruption age of ∼37 ka for the Firiplaka tephra, showing that explosive volcanism on Milos is ∼24 kyrs younger than previously thought. The four marine Nisyros tephras (N1 to N4) identified in part 1 (including the Upper (N1) and Lower (N4) Pumice) have ages of ∼57 ka, ∼63 ka, ∼69 ka, and ∼76 ka, respectively. Eruption ages for the Yali-1 and Yali-2 tephras are ∼55 ka and ∼34 ka, respectively. The Yali-2 tephra comprises two geochemically and laterally distinct marine facies. The southern facies is identical to the Yali-2 fall deposit on land but the western facies has slightly less evolved glass compositions. Overall, erupted plinian and co-ignimbrite fall tephra volumes range from 〈1 to 56 km3 (excluding possible caldera fillings and ignimbite volumes), and 80% of the eruptions had magnitude 5.5〈M≤7.2 (M=log(m)-7; m = erupted magma mass in kg). Twenty percent of the tephras represent 3.2〈M〈5.5 eruptions. The long-term average tephra magma mass flux through highly explosive eruptions of Santorini is estimated at ∼40 kg/s. The analogous data for the Kos-Yali-Nisyros volcanic complex is less-well constrained but similar to Santorini.

Description: The Christiana-Santorini-Kolumbo volcanic field (CSKVF) in the Aegean Sea is one of the most active volcano-tectonic lineaments in Europe. Santorini has been an iconic site in volcanology and archaeology since the 19th century, and the onshore volcanic products of Santorini are one of the best-studied volcanic sequences worldwide. However, little is known about the chronology of volcanic activity of the adjacent submarine Kolumbo volcano, and even less is known about the Christiana volcanic island. In this study, we exploit a dense array of high-resolution marine seismic reflection profiles to link the marine stratigraphy to onshore volcanic sequences and present the first consistent chronological framework for the CSKVF, enabling a detailed reconstruction of the evolution of the volcanic rift system in time and space. We identify four main phases of volcanic activity, which initiated in the Pliocene with the formation of the Christiana volcano (phase 1). The formation of the current southwest-northeast–trending rift system (phase 2) was associated with the evolution of two distinct volcanic centers, the newly discovered Poseidon center and the early Kolumbo volcano. Phase 3 saw a period of widespread volcanic activity throughout the entire rift. The ongoing phase 4 is confined to the Santorini caldera and Kolumbo volcano. Our study highlights the fundamental tectonic control on magma emplacement and shows that the CSKVF evolved from a volcanic field with local centers that matured only recently to form the vast Santorini edifice.

Description: Santorini caldera has had a long history of plinian eruptions and caldera collapses, separated by 20–40 kyr interplinian periods. We have carried out a study to constrain magma storage/extraction depths beneath the caldera. We analysed H 2 O in 138 olivine-, pyroxene- and plagioclase-hosted melt inclusions from plinian and interplinian products from the last 200 kyr, and CO 2 , S, Cl, F and D in various subsets of these. The dataset includes 64 inclusions in products of the Minoan plinian eruption of the late 17th century BCE. All the melt inclusions were ellipsoidal and isolated, with no textural evidence for volatile leakage. Mafic melt inclusions contain 1–4 wt % H 2 O and up to 1200 ppm CO 2 , 1200 ppm S, 2000 ppm Cl and 400 ppm F; silicic inclusions contain 2–7 wt % H 2 O, up to 150 ppm CO 2 , up to 400 ppm S, 2000–6000 ppm Cl and 600–1000 ppm F. The D values of 27 representative inclusions (–37 to –104) are intermediate between mantle and slab values and rule out significant H 2 O loss by hydrogen diffusion from olivine-hosted inclusions. H 2 O, S and Cl behave compatibly in melt inclusion suites varying from mafic to silicic in composition, showing that entrapment of many melt inclusions took place under volatile-saturated conditions. Most Santorini melts are saturated in a free COHSCl vapour phase at depths of less than ~10 km; the only exceptions are basaltic melts from a single interplinian eruption, which were volatile-undersaturated up to K 2 O contents of ~1 wt %. The rhyolitic melt of the Minoan eruption probably contained a free hypersaline liquid phase. H 2 O + CO 2 saturation pressures were calculated using suitably calibrated solubility models to estimate pre-eruptive magma storage depths. Magmas feeding plinian eruptions were stored at 〉4 km (〉100 MPa) and extracted over depth intervals of several kilometres. Plagioclase phenocrysts in rhyodacitic pumice from the Minoan eruption have cores containing melt inclusions trapped at depths up to 10–12 km (320 MPa), and rims (also orthopyroxene and clinopyroxene) containing inclusions trapped at 4–6 km (100–160 MPa). This records late-stage silicic replenishment of a 〈2 km thick shallow magma chamber, rather than extraction of melts syn-eruptively over the entire depth range. The plagioclase cores were carried from depth in the ascending melt, then overgrown by the rims in the shallow chamber. Exsolution of volatiles during ascent may have caused the replenishment melt to inject as a bubbly plume, causing mixing prior to eruption. This would explain (1) the homogeneity of the Minoan rhyodacitic magma, and (2) extraction of melt inclusions from the entire pressure spectrum during the first eruptive phase. Most silicic magmas feeding eruptions of the interplinian periods were stored in reservoirs at shallow depths (2–3 km) compared with those feeding the plinian eruptions (〉4 km). Melt inclusions from the ad 726 eruption of Kameni Volcano yield a pre-eruptive storage depth of ~4 km, which is similar to that estimated from geodetic data for the inflation source during the 2011–2012 period of caldera unrest; this supports a magmatic origin of the unrest. The level of pre- ad 726 magma storage beneath Kameni was deeper than that of earlier silicic interplinian eruptions, perhaps owing to changes in crustal stress caused by the Minoan eruption. Combined with previously published results, the melt inclusion data provide a time-integrated image of the crustal plumbing system. Mantle-derived basalts are injected into the lower crust, where they fractionate to produce evolved melts in bodies of hot crystal mush. Evolved residual melts separate from their parent mushes in the 8 to 〉15 km depth interval, then ascend rapidly into the upper crust, where they either crystallize or accumulate as bodies of eruptible, crystal-poor magma.

Description: Santorini volcano in the Aegean region (Greece) is characterized by andesitic- to silicic-dominated explosive activity and caldera-forming eruptions, sourced from magmatic reservoirs located at various structural levels beneath the volcano. There is a good understanding of the silica-rich magmatism of the island whereas the andesite-dominated volcanism and the petrogenesis of the parental mafic magmas are still poorly understood. To fill this gap we have performed crystallization experiments on a representative basalt from Santorini with the aim of determining the conditions of differentiation (pressure, temperature, volatile fugacities) and the parental magma relationship with the andesitic eruptive rocks. Experiments were carried out between 975 and 1040°C, in the pressure range 100–400 MPa, f O 2 from QFM to NNO + 3·5 (where QFM is quartz–fayalite–magnetite and NNO is nickel–nickel oxide), with H 2 O melt contents varying from saturation to nominally dry conditions. The results show that basalt phenocrysts within the basalt crystallized at around 1040°C in a magma storage reservoir located at a depth equivalent to 200–400 MPa pressure, with 3–5 wt % dissolved H 2 O, and f O 2 around QFM. Comparison with the xenocryst and phenocryst assemblages of the Upper Scoria 1 andesite shows that andesitic liquids are produced by fractionation of a similar basalt at 1000°C and 400 MPa, following 60–80 wt % crystallization of an ol + cpx + plag + Ti-mag + opx ± pig–ilm assemblage, with melt water contents around 4–6 wt %. At Santorini, the andesitic low-viscosity and water-rich residual liquids produced at these depths segregate from the parent basaltic mush and feed the shallow magma reservoirs, eventually erupting upon mixing with resident magma. Changes in prevailing oxygen fugacity may control the tholeiitic–calc-alkaline character of Santorini magmas, explaining the compositional and mineralogical differences observed between the recent Thyra and old eruptive products from Akrotiri.