A mineral and cumulate perspective to magma differentiation at Nisyros volcano, Aegean arc

Lavas and pyroclastic products of Nisyros volcano (Aegean arc, Greece) host a wide variety of phenocryst and cumulate assemblages that offer a unique window into the earliest stages of magma differentiation. This study presents a detailed petrographic study of lavas, enclaves and cumulates spanning the entire volcanic history of Nisyros to elucidate at which levels in the crust magmas stall and differentiate. We present a new division for the volcanic products into two suites based on field occurrence and petrographic features: a low-porphyricity andesite and a high-porphyricity (rhyo)dacite (HPRD) suite. Cumulate fragments are exclusively found in the HPRD suite and are predominantly derived from upper crustal reservoirs where they crystallised under hydrous conditions from melts that underwent prior differentiation. Rarer cumulate fragments range from (amphibole-)wehrlites to plagioclase-hornblendites and these appear to be derived from the lower crust (0.5–0.8 GPa). The suppressed stability of plagioclase and early saturation of amphibole in these cumulates are indicative of high-pressure crystallisation from primitive hydrous melts (≥ 3 wt% H2O). Clinopyroxene in these cumulates has Al2O3 contents up to 9 wt% due to the absence of crystallising plagioclase, and is subsequently consumed in a peritectic reaction to form primitive, Al-rich amphibole (Mg# > 73, 12–15 wt% Al2O3). The composition of these peritectic amphiboles is distinct from trace element-enriched interstitial amphibole in shallower cumulates. Phenocryst compositions and assemblages in both suites differ markedly from the cumulates. Phenocrysts, therefore, reflect shallow crystallisation and do not record magma differentiation in the deep arc crust.


Introduction
Arc magmas typically represent mixtures of one or more distinct volatile, melt and mineral components that are not necessarily in equilibrium with each other.A wide variety of macro-and microscopic disequilibrium textures is particularly common in more evolved arc rocks (andesites, dacites and rhyodacites) and includes, amongst others, resorbed or sieve textured phenocrysts, strong compositional zoning, replacement reactions and co-occurrence of olivine and quartz (e.g., Eichelberger 1978;Streck 2008).In addition to textural observations, in situ geochemical analyses of minerals and melt inclusions provide abundant evidence for open-system differentiation of arc rocks (e.g., Davidson et al. 2005;Ginibre et al. 2007).Mineral and melt heterogeneity in arc rocks implies that any bulk rock analysis will provide an average composition of multiple geochemically distinct components and inherently obscures a major part of the variation and complexity of a sample (e.g., Reubi and Blundy 2009).In contrast, the composition and textural relationships of mineral phases in arc magmas provide a wealth of information that can be used to investigate petrogenetic processes (e.g., Davidson et al. 2005).
There exists, however, a common mismatch between observed phenocrysts in extrusive products and cumulate assemblages that are inferred from whole rock geochemical trends (e.g., Arculus and Wills 1980;Annen et al. 2006).For example, the important role of amphibole in the evolution of arc magmas, as deduced from lava trace element compositions, is at odds with the common absence of amphibole phenocrysts (e.g., Davidson et al. 2007).It is increasingly recognised that arc magmas can undergo significant modification at the base of the arc crust without leaving a clear trace in their crystal record (e.g., Annen et al. 2006), which is supported by extensive evidence for cumulate-melt reactions in exposed arc root complexes (e.g., Jagoutz et al. 2007;Dessimoz et al. 2012;Bouilhol et al. 2015).Hence, there has been renewed interest in the use of cumulate suites to constrain differentiation processes at different levels in the arc crust (e.g., Tollan et al. 2012;Smith 2014;Stamper et al. 2014;Cooper et al. 2016).In particular, arc cumulates have the potential to record the first stages of magma differentiation that are rarely preserved in the erupted lavas.
We have undertaken a detailed petrographic study of Nisyros volcano, the easternmost volcanic centre of the Aegean arc (Greece), to further the understanding of its magmatic plumbing system and the differentiation of arc magmas in general.Nisyros is particularly well suited for such a study as it has erupted basaltic andesites to rhyodacites with a wide variety of mineral contents and textures over a < 160 kyr period (e.g., Wyers and Barton 1989;Francalanci et al. 1995;Vanderkluysen et al. 2005a).Moreover, abundant mafic inclusions are present in the Nisyros volcanic products, ranging from quenched enclaves to cumulate fragments.Cumulates vary from microcumulates, consisting of 5-50 crystals within a single thin section, to coarse plagioclase-hornblendite cumulate blocks.Earlier petrographic studies of Nisyros have provided representative mineral chemistry data (Di Paola 1974;Wyers and Barton 1989;Seymour and Vlassopoulos 1992) or were restricted to selected volcanic units (Spandler et al. 2012;Braschi et al. 2014;Zouzias and St Seymour 2014).We expand on these studies and aim to provide a more comprehensive overview of the subaerial volcanic history of Nisyros.In addition, we provide the first detailed petrographic study of (micro)cumulate fragments from Nisyros and evaluate the insights they can provide into the structure of its plumbing system.Of particular interest are the hornblende-gabbro and hornblende-wehrlite cumulates that are potentially derived from the lower arc crust and hence offer a unique opportunity to study cumulate-melt reactions that drive the differentiation of calc-alkaline arc magmas.

Geological setting
Nisyros is the easternmost volcanic centre of the Aegean arc, which has formed in response to the northward subduction of the African plate underneath the Turkish-Aegean microplate.Volcanism in the eastern section of the arc started with the emplacement of dacitic to rhyolitic domes on Kos from ca. 2.7-0.5 Ma (Pe-Piper and Moulton 2008) and culminated in the large explosive eruption that deposited the Kos Plateau Tuff at 161 ± 1 ka (Smith et al. 1996;Bachmann et al. 2007).The entirely volcanic island of Nisyros is situated ca.20 km to the south of Kos and has a diameter of ca. 8 km with a central caldera with a diameter of ca. 4 km (Fig. 1).The subaerial volcanic products of Nisyros are inferred to be younger than 161 ka as the Kos Plateau Tuff was found immediately overlying the non-volcanic basement in drill cores in Nisyros' caldera (Volentik et al. 2005a) and references therein).A short summary of the eruptive history of Nisyros is given here, following the detailed account provided by Vanderkluysen et al. (2005b) and using the classification and abbreviations (in capital letters, e.g., EMB, LF7) of the volcanic units as presented by Volentik et al. (2005b).
The oldest volcanic rocks on Nisyros are the Holaki pillow lavas and hyaloclastites (unit HO), which possibly predate the KPT eruption (Keller 1971).Subaerial activity started with the construction of a small stratocone composed of basaltic andesitic to dacitic lava flows (LF1-LF6) and pyroclastic deposits (e.g., MS) that are exposed in the caldera wall.The presence of intercalated lacustrine sediment deposits suggests that small caldera collapses were followed by periods of volcanic quiescence.The subsequent phase of activity was characterised by the eruption of lava flows and explosive events (LF7, PA, LF8, XO) from multiple different mono-and polygenetic eccentric vents, including the Lies tuff cone on the eastern side of the island.Following these paroxysmal events, the style of volcanism on Nisyros underwent a gradual change from effusive and mildly explosive activity to the emplacement of (rhyo)dacitic domes and lava flows and larger Plinian eruptions.These porphyritic (rhyo)dacites are characterised by the conspicuous presence of mafic enclaves with a large range in phenocryst compositions and disequilibrium textures, which are indicative of magma mingling (Braschi et al. 2012(Braschi et al. , 2014;;Zouzias and St Seymour 2014).Lava flow 7 (LF7), underlying the Lies tuff cone, is the oldest studied enclave-bearing porphyritic dacite, followed by the Emborio domes (EMB) and the calderaforming Lower and Upper Pumice eruptions (LP and UP, respectively; Limburg and Varekamp 1991;Vanderkluysen Page 3 of 23 95 et al. 2005b;Longchamp et al. 2011).The latter are separated by the emplacement of a ca. 1 km 3 rhyodacitic, highly viscous block lava flow that covers the southeast of Nisyros: the Nikia flow (NLF; Volentik et al. 2005a).After the eruption of the Upper Pumice, the newly formed caldera was partly filled by the emplacement of six dacitic domes (the Profitis Ilias domes; PFI) that comprise up to 20 vol% of basaltic-andesitic enclaves.

Analytical methods
The volcanic deposits that are analysed in this study were sampled on the basis of the geological map of Volentik et al. (2005b) and the volcanic units defined therein, with the aim to cover most of the volcanic history.Whole rock major element compositions of the samples were determined by X-ray fluorescence spectroscopy (XRF) on fused glass beads.LOI was determined by igniting sample powders for 1 h at 1000 °C.Fused glass beads were made by mixing ignited sample material with Li 2 B 4 O 7 /LiBO 2 (1:4 dilution) and measured on a Panalytical MagixPro XRF instrument at Vrije Universiteit Amsterdam.Results are expressed on a volatile-free basis with Fe as total ferrous iron (FeO*).Mineral compositions were measured on carbon-coated polished sections by wavelength dispersive electron microprobe (EMP) using a JEOL JXA-8800 M at the Vrije Universiteit Amsterdam and a JEOL JXA-8530F at the University of Bristol.Operating conditions were similar for the two instruments.The EMP was operated at an acceleration voltage of 15-20 kV and a beam intensity of 20-25 nA.A focussed beam was used for the analysis of olivine, pyroxene and amphibole while plagioclase was analysed with a spot size of 10 µm.All analyses were matrix-corrected using the ZAF (VU) or φ(rZ) method (Bristol) and calibrated against natural and synthetic mineral standards.Accuracy of the calibration was tested by regular analysis of mineral standards and calibration was accepted in case the standard values were reproduced at least within 1% (RSD) for each element.Mineral compositions were checked for stoichiometry by normalising to four atoms per formula unit (apfu) O for olivine, six apfu O for pyroxene and eight apfu O for plagioclase.Amphibole was classified on the basis of 13 apfu O (Leake et al. 1997).Analyses were excluded if the sum of cations was < 99.5 or > 100.5% than the stoichiometric value (< 1% of the analyses).
Mineral trace element concentrations of selected samples were determined by laser-ablation inductively coupled plasma mass spectrometry (LA-ICPMS) at the Westfälische Wilhems-Universität Münster.Ablation was achieved with a 193-nm excimer laser-ablation system (Analyte G2, Photon Machines) connected to a ThermoFisher Element 2 mass spectrometer.The repetition rate of the laser was set to 5-10 Hz with an energy density of ~ 4 J/cm 2 and a spot size of 35 µm.Optimisation of the system was aimed at PA comprises pyroclastic deposits with rapid switching between layers with HPRD and LPA affinity and is thus transitional between the two suites maximum sensitivity for 139 La and 232 Th on the NIST 612 reference material while keeping oxide production below 0.1% ( 232 Th 16 O/ 232 Th).A single measurement lasted 60 s including 20 s background measurements.GLITTER software (Griffin et al. 2008) was used to calculate concentrations for the measured elements with NIST 612 reference material as an external standard and internal standardisation to 43 Ca.Laser-ablation spots were put on top of EMP spots as much as possible and hence Ca contents obtained by EMP could be used for normalisation.Plagioclase is characterised by significant heterogeneity in Ca content on a scale much finer than the LA spot size and these variations would introduce a significant uncertainty in the internal normalisation.Based on the stoichiometry of feldspar, we used measured Si/Ca to derive an integrated Ca and anorthite content for each LA spot and used this Ca content to obtain corrected trace element concentrations (see supplementary material).Glass reference material BCR-2G was used to monitor precision and accuracy, and reproducibility was better than 10% (2 RSD) for elements of interest (see supplementary material).

Petrography
Two suites of volcanic rocks can be identified on Nisyros that are easily distinguished on the basis of petrography and field relationships.The main distinguishing feature is the clear evidence for magma mingling, such as the abundance of magmatic enclaves and disequilibrium textures, in the first suite while the second suite is characterised by an on average lower crystal content and lack of macroscopic evidence for mingling.The two suites have been somewhat haphazardly classified as the "pre-" and "post-caldera" magmatic series by Buettner et al. (2005) and Zellmer and Turner (2007).This terminology is, however, not very exact as there have been multiple caldera-forming events on Nisyros.In addition, we have identified multiple units that predate the last caldera collapse as part the magma-mingling suite, including lava flow 7 (LF7), the Emborio domes (EMB) and the Nikia flow (NLF).Hence, we prefer to classify the two suites based on their petrography and field occurrence in a lowporphyricity andesite (LPA) suite and high-porphyricity (rhyo)dacite (HPRD) suite (Fig. 1).The LPA suite comprises most of the oldest volcanic units emplaced during the early cone-building stage, including the lower part of the succession of lava flows that is exposed in the caldera wall, and is dominated by lava flows and mildly explosive pyroclastic and scoria deposits.After the episode of violent eccentric activity following the initial cone building stage, the enclave-bearing HPRD suite became the dominant magmatic series.HPRD units occur as domes (EMB, PFI), pyroclastic deposits related to explosive eruptions (LP and UP) and viscous block lava flows (LF7, NLF).

LPA suite
The low-porphyricity andesite suite comprises basalticandesites, andesites and rare dacites.The majority of the LPA units form a group of petrographically homogeneous andesites.These vesicular (0-70 vol% vesicles), low-porphyricity andesites (0-10 vol% phenocrysts) have a glassy to fine-grained, plagioclase-dominated microlite groundmass with phenocrysts of plagioclase, olivine, clinopyroxene, Fe-Ti-oxides and rare orthopyroxene (Fig. 2).Plagioclase is generally euhedral but displays a wide variety of zoning patterns.Rim compositions vary largely between and within samples in the range An 57-81 and samples often contain both normally and reversely zoned plagioclase crystals.Crystals with complex zoning are typically more subhedral and contain sieve textured cores up to An 88 .Olivine is the most common phenocryst phase after plagioclase, is subto euhedral and displays normal zoning from homogeneous core compositions (Fo 69-74 ) to variable, more Fe-rich rims.Sub-to euhedral clinopyroxene crystals are in general normally zoned from Mg# 85 cores to Mg# 69-76 rims, but can also show reverse zoning to Mg# 80 rims.
Dacites of the LPA suite are restricted to only two samples (LF1 and a dyke therein) and contain an equilibrium mineral assemblage of plagioclase, amphibole, Fe-Ti-oxides and orthopyroxene (Fig. 2c).Apart from opacite reaction rims related to decompression, amphibole phenocrysts are euhedral and in apparent equilibrium with the matrix.Plagioclase is generally euhedral and shows normal zoning with a limited range of An 40-45 , but a few crystals have resorbed, sieve textured cores.Orthopyroxene is sub-to euhedral and homogeneous at Mg# 60-62, and Fe-Ti-oxides are more abundant than in the andesites (~ 2 and ~ 1 vol%, respectively).Unit ALF represents a low-porphyricity rhyodacite that is hence treated here as part of the LPA suite.Basaltic andesites in the LPA suite are different from the andesites in that they generally have a higher crystal content (10-40 vol%;Francalanci et al. 1995).Clinopyroxene is the dominant Fe-Mg phase and occurs as euhedral, up to 5 mm large normally zoned crystals in LF4 and the Holaki pillow lavas (HO).Core domains of these large crystals are the most magnesian on Nisyros (Mg# 91).Clinopyroxene in HO commonly has resorbed, sieve textured cores with normal zonation from Mg# 80 to 75 with euhedral rims that range from Mg# 85 to 80, as also described by Spandler et al. (2012).The basaltic andesites either contain olivine, zoned from Fo 81-83 cores to Fo 70 rims in the Aghia Basilei lava flows (BLF), or orthopyroxene with Mg# 78-80 in LF4 and HO, but the two minerals have not been observed together in one thin section.Fe-Ti-oxides are very rare in the basaltic andesites.

HPRD suite
The high-porphyricity dacites and rhyodacites are distinguished from the LPA suite by their higher crystal content (20-40 vol%) and more felsic composition (> 64 wt% SiO 2 ).The most striking macroscopic feature of the HPRD suite is the abundance of enclaves and cumulate fragments in these units.Dacites cannot be easily distinguished from rhyodacites on the basis of petrography as they have a common texture of coarse plagioclase (up to 1 cm), orthopyroxene, Fe-Ti-oxides and subordinate clinopyroxene in a plagioclase-dominated microlite groundmass (Fig. 2).Plagioclase is the most abundant phenocryst phase (10-35 vol%) and can be roughly grouped into euhedral equilibrium plagioclase and sub-to anhedral disequilibrium plagioclase.Braschi et al. (2014) recognise five distinct plagioclase types in the Profitis Ilias domes (PFI) on the basis of crystal habit, zoning patterns and presence of sieve textured cores and/ or dusty rims.These authors relate the wide variety in plagioclase textures to magma mingling and the transfer of plagioclase crystals from the host rock to the enclaves and vice versa.Compared to the other HPRD units, plagioclase crystals in the PFI domes are more anhedral and disequilibrium textures are significantly more common.Euhedral, oscillatory zoned plagioclase is dominant in LF7, the Emborio domes (EMB) and the Nikia flow (NLF).Typical plagioclase rim compositions in the HPRD dacites (LF7, EMB) are An 40-52 , compared to ca.An 30 in the rhyodacites (NLF, PFI); core compositions and zoning patterns are variable, with compositions ranging up to An 90 .Orthopyroxene crystals are generally euhedral, commonly associated with Fe-Tioxides and homogeneous at Mg# 58-63.Clinopyroxene is subordinate to orthopyroxene in the HPRD suite and occurs as an-to euhedral isolated grains with Mg# between 67 and 73, sometimes with a reaction rim.Amphibole has not been found as phenocrysts and occurs exclusively in disaggregated fragments of enclaves.

Enclaves and cumulates in the HPRD suite
Mafic inclusions in the HPRD suite display a large textural and mineralogical variation and can be divided into nearly holocrystalline plutonic xenoliths and quenched enclaves.The latter commonly have chilled margins, contain up to 15 vol% vesicles and display a fine-grained texture of acicular, occasionally spherulitic, amphibole (up to 1 cm in length) and plagioclase with groundmass Fe-Ti-oxides and diktytaxitic voids (Fig. 2).This texture indicates that the enclaves formed in response to quench crystallisation of a hot, mafic melt injected into a cooler reservoir (e.g., Bacon 1986).Enclaves are generally rounded in EMB and LF7 but can be more irregular in PFI and range in size from < 1 mm to 50 cm (Fig. 2e).Of the studied units, the largest enclaves and highest abundance are found in the PFI domes (see Braschi et al. 2012 for a detailed description).Crystals in the enclaves that appear to predate the quench crystallisation are referred to as phenocrysts.The PFI enclaves contain less than 5 vol% phenocrysts that predominantly comprise dusty, sieve textured plagioclase and minor olivine and clinopyroxene.Enclaves in LF7, however, contain reasonably abundant clinopyroxene, olivine and calcic plagioclase phenocrysts (in total 20-30 vol%) set in a groundmass of acicular amphibole, plagioclase and orthopyroxene (Fig. 3c).Plagioclase phenocrysts are strongly zoned from An 87-91 cores to An 58-63 rims.Olivine also displays normal zoning from Fo 80 cores to Fo 67 rims.Orthopyroxene occurs both as acicular crystals in vesicles as well as in complex, vermicular intergrowths with plagioclase.
On the other side of the mafic inclusion spectrum are plutonic xenoliths.Textural and chemical evidence indicates that minerals in these xenoliths formed through fractional crystallisation and that their bulk composition forms a subtractive assemblage that is unlike any of the whole rock compositions found on Nisyros.Hence, these plutonic xenoliths are referred to as cumulates following the classification of Irvine (1982), which is applicable irrespective of the presence of intercumulus crystals or melt.Cumulates vary in size from > 50 to < 0.5 mm and occur both as discrete nodules found in pyroclastic deposits (Lies tuff cone, Lower Pumice) and as microcumulates.The latter are defined as clusters of 5-50 crystals in HPRD samples that are only visible in thin section.Plagioclase, clinopyroxene and olivine are the dominant mineral phases, together with subordinate amphibole, Fe-Ti-oxides and trapped interstitial melt.Microcumulates typically show strong disequilibrium textures such as normal zoning and orthopyroxene mantles around olivine when hosted in a (rhyo)dacitic melt (Fig. 3d).Crystal-rich (~ 70 vol%) plagioclase-hornblendite cumulates are restricted to the Nikia flow (NLF).These cumulate blocks are rounded and characterised by a general absence of chilled margins, but sometimes display a vesicle-rich halo around the enclave in the rhyodacitic host rock.Euhedral, coarse (up to 1 cm) amphibole constitutes 25-30 vol% of the minerals and lath-shaped plagioclase, clinopyroxene (0-5 vol%), Fe-Ti-oxides, interstitial glass and ca.10% vesicles fill up the space between the coarse amphibole grains.Amphibole is zoned but lacks an opacite rim and is found to be growing at the expense of clinopyroxene (Fig. 3f).The different types of cumulates (Type 1-3) are described in detail in the discussion.

Whole rock composition
Major element data for the studied samples overlap with previously reported data for Nisyros (Di Paola 1974;Wyers and Barton 1989;Seymour and Vlassopoulos 1992;Francalanci et al. 1995;Buettner et al. 2005;Vanderkluysen et al. 2005a;Zellmer and Turner 2007;Braschi et al. 2012;Spandler et al. 2012).Nisyros lavas range from medium-K basaltic andesites to high-K rhyodacites (Fig. 4).The samples show typical fractional crystallisation trends with decreasing Al 2 O 3 , CaO, MgO and Mg# (molar Mg/[Mg + Fe]) and increasing Na 2 O and K 2 O with SiO 2 content.Basaltic andesites in the LPA suite have Mg# up to 68 while major element systematics rule out significant accumulation of Mg-rich phases such as clinopyroxene or olivine in these samples.For this reason, the high Mg# is compatible with the LPA basaltic andesites being primitive melts despite SiO 2 contents of 53.5-57.5 wt%, and hence they can be classified as primitive andesites (following the definition of Kelemen et al. 2003).Mafic enclaves hosted in the HPRD suite largely overlap in composition with LPA basaltic andesites, although the former generally have lower Mg# (50-60) and SiO 2 contents, and Al 2 O 3 contents up to 20.5 wt% are found in enclaves in the Profitis Ilias domes (PFI; Braschi et al. 2012).The lack of significant trace element and isotopic variations between the most primitive LPA and HPRD samples suggests that there is no pronounced difference in primary melts feeding the two suites (Klaver et al. 2016).
Cumulates in the HPRD suite have lower SiO 2 content than any of the extrusive products.The bulk composition of the plagioclase-amphibolite heteradcumulates in the Nikia flow is relatively similar to HPRD enclaves due to the presence of intercumulus rhyolitic melt but is nevertheless distinct in, for instance, higher TiO 2 contents.A plagioclase-clinopyroxene-olivine and an amphibole-wehrlite orthocumulate are characterised by higher Mg# and lower SiO 2 contents than the extrusive products, but have highly variable Al 2 O 3 contents depending on the presence of plagioclase.

Mineral compositions
An overview of the variation in plagioclase, pyroxene and olivine compositions in the Nisyros volcanic rocks and cumulates is shown in Fig. 5 and amphibole major element variations are presented in Fig. 6.The samples have been grouped according to their petrographic features in a lowporphyricity andesite (LPA) and high-porphyricity rhyodacites (HPRD) suite and further subdivided on the basis of bulk rock composition.Detailed major element variation diagrams for all minerals are provided in the online supplementary material.

Clinopyroxene
Nisyros clinopyroxenes are diopsides to augites that range in composition from Mg# 91-65 and are characterised by a large variation in Al 2 O 3 , TiO 2 , CaO and Cr contents.In general, clinopyroxene with the highest Mg# and up to 6000 ppm Cr is restricted to the LPA basaltic andesites; clinopyroxene in HPRD suite cumulates has lower Mg# (89-75) and Cr contents (< 3000 ppm).The LPA andesites and HPRD (rhyo)dacites host clinopyroxene with Mg# generally < 80 and Cr contents below detection limit.Clinopyroxene TiO 2 contents vary between the two suites: LPA andesite clinopyroxene contains 0.5-1.5 wt% TiO 2 while the TiO 2 content is < 0.5 wt% in the HPRD suite.Aluminium concentrations are variable: mostly < 3.5 wt% Al 2 O 3 in the LPA suite, but up to 9 wt% in cumulate crystals.

Amphibole
Amphibole displays a limited variation in SiO 2 (40-45 wt%) and Al 2 O 3 (8-15 wt%) concentrations, predominantly at the low-SiO 2 /high-Al 2 O 3 end of the range displayed by amphibole in arc volcanoes (e.g., Ridolfi et al. 2010).A bimodal amphibole population with low-temperature magnesio-hornblende and high-temperature tschermakitic pargasite, which Fig. 4 Major element variation of Nisyros eruptive products.Cumulates for which whole rock data are available, either by direct bulk analysis or by point counting in the case of homogeneous, unzoned cumulate phases, are included.The Nikia flow (NLF) Type 3 cumulates represent bulk analyses that are affected by the infiltration of rhyolitic melt and do, therefore, not reflect the true bulk composition of the subtractive mineral assemblage.Grey squares are literature data, see text for references.Data are reported on a volatile-free basis with all Fe as FeO*.In a, the low-K to high-K classification is after Le Maitre et al. (1989).In b, the grey bar represents the range of primary melts in equilibrium with mantle olivine (Fo 90 ;Roeder and Emslie 1970) is characteristic for many intermediate arc lavas (e.g., Mt.Hood, Unzen, Soufrière Hills; see Kent 2014 for a compilation), has not been observed (Fig. 6).Instead, amphibole compositions are restricted to high Al/Si magnesio-hastingsites and tschermakitic pargasites.Amphibole in the HPRD suite enclaves overlaps with the high-temperature population commonly found in arcs, while Al 2 O 3 is lower and TiO 2 and Na 2 O contents are higher in LPA dacite amphibole (Fig. 6).In addition, acicular amphibole in the Emborio dome enclave is compositionally similar to the coarse amphibole in the LPA dacites, indicating that the chemical composition of the amphiboles does not correlate with their morphology.

Plagioclase
Plagioclase varies from An 24 to An 91 in the Nisyros sample suite.Orthoclase contents (Or 0.1-5 ) are higher in LPA suite (basaltic) andesites than in HPRD enclaves and host lavas in the range An 50-80 but overlap at An >80 .The highest An% is found in cumulates hosted in the HPRD suite and in LPA basaltic andesite phenocrysts.In the cumulates, plagioclase can display zoning and resorption from calcic (An >80 ) cores to more sodic rims.Reverse zoning is rare, in contrast to plagioclase phenocrysts in the LPA suite that often preserve complex and oscillatory zoning patterns (Fig. 2d).In the LPA suite, crystals with normal and reverse zoning to a common rim composition frequently occur within the same thin section.HPRD suite host lavas and enclaves often contain similar, often bimodal plagioclase populations.

Olivine and orthopyroxene
Olivine in the Nisyros sample suite ranges from Fo 86 to Fo 60 .The most forsteritic olivine (Fo 86 ) is found in a Two separate groups of orthopyroxene compositions can be distinguished that constitute a bimodal distribution in Mg#.Orthopyroxene is most abundant in the HPRD hosts and LPA dacites where it displays a narrow range in Mg# from 60 to 65.Although orthopyroxene Mg# in the HPRD and LPA dacites overlaps, the two groups are separated by higher Al 2 O 3 and TiO 2 in the LPA dacite orthopyroxenes, which is similar to the clinopyroxene systematics.Orthopyroxene ranges from Mg# 78 to 83 in the LPA basaltic andesites; orthopyroxene in ALF and cumulate fragments has Mg# intermediate between the two groups.The Al 2 O 3 content of Nisyros orthopyroxene ranges from 1.0 to 1.5 wt% in the LPA basaltic andesites to 0.3-1.0wt% in the HPRD suite and LPA dacites; TiO 2 contents are generally lower than 0.4 wt%.

Distinction between the LPA and HPRD suites
We propose a new division of the Nisyros volcanic deposits based on field relationships and petrographic features.The main distinction between the two suites is the pronounced macro-and microscopic evidence for disequilibrium and magma mixing and mingling in the HPRD suite, and lack thereof in the LPA suite.For example, both suites include dacites, but the rare LPA dacites contain < 15 vol% phenocrysts set in a glassy groundmass while HPRD dacites have a porphyritic texture with 20-40 vol% crystals up to 1 cm in size and a coarse microlite matrix.Phenocryst phases also differ between the two as illustrated by the lack of primary amphibole and presence of clinopyroxene in the HPRD dacites.Moreover, the HPRD dacites contain mafic enclaves and a wide variety of microcumulates and xenocrysts.Distinguishing the two suites is thus straightforward in the field and on the hand specimen to thin section scale.
Porphyritic andesites and dacites are common eruptive products of arc volcanoes and often show textural and geochemical evidence for (extensive) hybridisation (e.g., Humphreys et al. 2009;Reubi and Blundy 2009;Cooper and Kent 2014;Kent 2014;Blum-Oeste and Wörner 2016).Such hybrid andesites are believed to tap shallow, crystalrich mush systems that have stalled in the upper crust due to decompression-induced degassing and rapid crystallisation (e.g., Blundy and Cashman 2005;Annen et al. 2006;Bachmann and Huber 2016).Due to their high viscosity, these mush systems spend the majority of the time close to the solidus in an uneruptible state, but can be rapidly remobilised and erupted in response to recharge by a volatilerich mafic magma (e.g., Cooper and Kent 2014;Parmigiani et al. 2014;Bergantz et al. 2015).Field relationships and textures of the HPRD (rhyo)dacites are very similar to these hybrid, reactivated mushes.High crystal contents, a coarse microlite groundmass and plagioclase-dominated mineralogy are diagnostic features.Moreover, the HPRD lavas include abundant mafic enclaves with a texture of acicular amphibole and plagioclase that attests to quench crystallisation upon injection in the cooler (rhyo)dacitic host (Bacon 1986;Vanderkluysen et al. 2005a;Braschi et al. 2014) and hence supports mafic recharge as a likely eruption trigger (Bachmann et al. 2012;Braschi et al. 2012Braschi et al. , 2014)).Magma mingling and disaggregation of the enclaves leads to crystal transfer from the enclaves to the (rhyo)dacitic host and vice versa and can thus account for the complex zoning patterns recorded by HPRD crystals as well as similar, often bimodal plagioclase distributions in the enclaves and hosts (Humphreys et al. 2009;Braschi et al. 2014;Bergantz et al. 2015).Given their textural and compositional similarity, it is likely that all HPRD units originated from a shallow crystal-rich mush system that was reactivated by mafic recharge and that this process was not restricted solely to the youngest PFI domes on Nisyros.
The crystal-poor, enclave-free LPA andesites and dacites form a striking contrast with the HPRD suite.
Andesites can be generated through extensive hybridisation of felsic and mafic melts in evolved mush systems (e.g., Reubi and Blundy 2009;Kent 2014) but this process is not tenable for the LPA andesites.The main argument against a hybrid origin for the LPA andesites and dacites is the lack of geochemical evidence for mixing (e.g., Wyers and Barton 1989;Vanderkluysen et al. 2005a).The LPA andesites and dacites do not fall on mixing lines between HPRD rhyodacites and primitive samples as, for instance, in a K 2 O versus SiO 2 diagram (Fig. 4).In addition, a hybrid origin is excluded by the low crystal content and lack of macro-and microscopic evidence for magma mixing between the required mafic and felsic components.Disequilibrium textures related to mixing between primitive and evolved melts, such as seen in the HPRD suite, are not present in the LPA andesites and dacites.There is no evidence for a felsic mixing component as no crystals with a HPRD (rhyo)dacite signature are found in the LPA crystal record.Plagioclase and clinopyroxene in the LPA andesites have An% and Mg# intermediate between the LPA basaltic andesites and HPRD rhyodacites, but minor and trace element contents in these phases are distinct.This is best illustrated by the TiO 2 content of clinopyroxene, which is higher in LPA andesites than in both the LPA basaltic andesites and HPRD suite and thus inconsistent with crystallisation from a hybrid andesitic melt (Fig. 7).In addition, clinopyroxene in LPA andesites has higher Al 2 O 3 and lower La/Sm (see supplementary material) compared to the HPRD suite.

Equilibrium mineral assemblages in the HPRD suite
The HPRD suite displays a wide variation in crystal components that range from phenocrysts that are in apparent equilibrium with the melt to antecrysts, xenocrysts, microcumulates and crystals derived from disaggregated enclaves, which are in textural and compositional disequilibrium.As these disequilibrium assemblages carry important information on the structure of Nisyros' plumbing system, we endeavour to identify mineral assemblages that appear to be in inter-mineral equilibrium, even though these might not be in equilibrium with the melt in which they are hosted.Cumulate fragments and microcumulates in the HPRD suite can be divided into three groups that are present in varying abundance, which are discussed below.An overview of these mineral assemblages is given in Table 1.

HPRD phenocryst assemblage
The HPRD (rhyo)dacites have a common equilibrium phenocryst assemblage of plagioclase (An 20-60 ), orthopyroxene (Mg# 60-64), ilmenite and clinopyroxene that often occur in clusters (Fig. 2).In contrast, amphibole phenocrysts are conspicuously absent despite being a common phase in similar (rhyo)dacitic arc magmas worldwide.Such magmas often display a bimodal distribution of amphibole compositions with a low-temperature population in equilibrium with the felsic host and a high-temperature population derived from the mafic recharge magma (Fig. 6; e.g., Coombs et al. 2013;Kent 2014;Kiss et al. 2014).Rare amphibole crystals in the HPRD (rhyo)dacites are resorbed and compositionally distinct from low-T amphibole but overlap with hightemperature compositions found in the enclaves, and thus appear to be derived from the mafic enclaves (Wyers and Barton 1989).The lack of primary amphibole likely results from a combination of higher magmatic temperatures and lower H 2 O contents compared to other hybrid arc magmas that place the HPRD (rhyo)dacites outside the stability field of amphibole (e.g., Bachmann et al. 2012).

Type 1 cumulates
The most common Type 1 cumulates invariably contain Anrich plagioclase (An 80-90 ) and clinopyroxene (Mg# 73-86), with or without olivine, amphibole and Fe-Ti-oxides.Their main distinguishing feature is the co-occurrence of calcic plagioclase (~ An 90 ) and relatively evolved (Mg# ≤ 85) clinopyroxene and olivine.The presence of interstitial amphibole and bimodal plagioclase compositions suggest that several of the Type 1 cumulates have reacted with percolating melts.Olivine and clinopyroxene Mg# are generally lower in the cumulates with interstitial amphibole.Type 1 cumulates occur as microcumulates in the HPRD enclaves and lavas and as fragments in pyroclastic deposits of the Lies tuff cone (Fig. 3).When entrained in a (rhyo)dacitic melt, the microcumulates are in disequilibrium as indicated by pronounced normal zoning of olivine, clinopyroxene and plagioclase and the presence of orthopyroxene mantles around olivine towards contacts with the host melt (Fig. 3d).
The preservation of Fe-Mg zoning profiles in olivine indicates that cooling was sufficiently rapid to inhibit complete re-equilibration.Diffusion modelling of Mg-Fe zoning in olivine can be used to constrain timescales of magma mixing (e.g., Costa and Dungan 2005;Turner and Costa 2007), but this was not undertaken due the scarcity of suitable olivine grains.Nevertheless, a timescale of ≪ 1 kyr between mixing and eruption can safely be assumed based on the pronounced zoning in olivine in these micro-enclaves.Similar zoned olivines from the Nea Kameni dacites, Santorini, suggest that less than 1 year has elapsed between mingling of the olivine crystals into the (rhyo)dacitic host and eruption (Martin et al. 2008).

Type 2 cumulates
Rare Type 2 cumulates comprise plagioclase-free or -poor assemblages that contain more primitive (higher Mg#) olivine than Type 1 cumulates.Type 2 microcumulates typically occur as clusters of anhedral crystals (< 10 grains) and are predominantly found in the PFI domes and enclaves.The exception is a larger hornblende-wehrlite cumulate fragment (NCX-06; Fig. 3e).This fragment consists of roughly equal proportions of olivine and clinopyroxene and ~ 10 vol% amphibole.It has a non-equilibrated texture characterised by euhedral olivine and poikilitic clinopyroxene.Amphibole typically fills the intercumulus space and the cusps between adjacent grains, but does not appear to be forming at the expense of clinopyroxene.No replacement of clinopyroxene with amphibole along cracks or cleavage planes is observed, in contrast to the Type 3 cumulates described below.A single isolated plagioclase crystal (An 86 ) is present, as are rare plagioclase inclusions in amphibole, but the sample is essentially plagioclase-free (< 0.5 vol%).Apart from the general scarcity of plagioclase, Type 2 cumulates are distinguished by high Al 2 O 3 contents in clinopyroxene (up to 7 wt% in individual crystals; Fig. 8).Type 1 and 2 cumulates are clearly distinct, but their occurrence is not mutually exclusive and they can occur in the same unit (e.g., in PFI and LF7).

Type 3 cumulates
The Type 2 cumulates transition into Type 3 plagioclase + amphibole ± clinopyroxene cumulates that have been found exclusively in the Nikia flow (NLF).Amphibole in these cumulates occurs in various forms.It replaces clinopyroxene along cracks and cleavage planes but also mantles clinopyroxene with less clear textural evidence for reaction-replacement (Fig. 9).Amphibole, clinopyroxene, magnetite and plagioclase cores are in apparent equilibrium with each other in the centre of crystal clusters.Euhedral plagioclase and amphibole crystals extend outwards and are in disequilibrium with the rhyolitic host melt, as indicated by resorption of amphibole and strong zoning of plagioclase from An 80 cores to An <55 rims (Fig. 9).The preservation of zoning profiles indicates that infiltration of the rhyolitic melt occurred shortly before eruption.Similar clinopyroxene-amphibole replacement textures have been observed in arc root complexes (e.g., Dessimoz et al. 2012;Bouilhol et al. 2015) and cumulate nodules in arc lavas (Smith 2014;Cooper et al. 2016).High Mg# (67-76; with all Fe as Fe 2+ ), low trace element contents and Eu/Eu* > 0.9 (Fig. 6) suggest that hornblende in the Type 3 cumulates crystallised from relatively primitive melts, in contrast to interstitial amphibole in the Type 1 cumulates.Clinopyroxene in the Type 3 cumulates has the highest Al 2 O 3 content found on Nisyros (up to 9 wt%).

Constraining the origin of the Nisyros cumulates
Phenocrysts in the LPA and HPRD suites broadly overlap with cumulate crystals in terms of An% (plagioclase), Fo% (olivine) and Mg# (pyroxenes and amphibole; Fig. 5).When considered in more detail, however, there are striking differences between phenocryst and cumulate assemblages in terms of paragenesis and mineral compositions (Fig. 8), as is often found in arc suites (e.g., Arculus and Wills 1980;Stamper et al. 2014).Key characteristics of the cumulate suite include the co-occurrence of high-An plagioclase with low-Fo olivine in Type 1 cumulates, the pronounced increase in clinopyroxene Al 2 O 3 content with decreasing Mg# in Type 2 and 3 cumulates, and overall higher clinopyroxene CaO contents compared to phenocrysts in the LPA suite.These differences likely reflect variations in depth, temperature, water content and oxygen fugacity at which magmas were stored and differentiated in the crust.Constraining these parameters is crucial to understand the magmatic plumbing system of Nisyros.Although the identification of equilibrium mineral assemblages is an important first step, determining the intensive parameters and in particular crystallisation pressures of cumulate suites and arc lavas is notoriously difficult (e.g., Stamper et al. 2014).The singleamphibole thermobarometer of Ridolfi et al. (2010) is appropriate to the Nisyros samples, but we refrain from using this barometer in line with concerns raised by Erdmann et al. (2014) and Putirka (2016).To put constraints on the depth of magma storage and differentiation, we therefore follow the approach of Blundy and Cashman (2008) and Stamper et al. (2014) of correlating observed equilibrium mineral assemblages with those produced in crystallisation experiments.An overview of experimental studies performed on appropriate starting compositions at conditions relevant to subduction zone magmatism is given in Table 2, and Fig. 10 shows a phase diagram for hydrous primary arc magmas compiled from these studies.

Hydrous high-pressure crystallisation
The Type 2 and 3 cumulates form a textural and compositional continuum from (hornblende-)wehrlites to plagioclase-hornblendites that we interpret as the result of progressive fractional crystallisation of a primitive, hydrous melt.Olivine is the first phase to crystallise followed by clinopyroxene, as indicated by the inclusion of euhedral olivine in larger clinopyroxene crystals (Fig. 3e).Amphibole is the next phase to appear, although the rare (< 0.5 vol%) presence of plagioclase in hornblende-wehrlite NCX-06 suggests it can be accompanied by minor plagioclase crystallisation.The early stabilisation of amphibole is reflected by its high Mg# (all Fe as Fe 2+ ), which ranges from 73 in plagioclasehornblendites up to 83 in hornblende-wehrlite NCX-06.Such magnesian amphibole is not commonly found in arc phenocrysts or cumulates (e.g., in the well-studied Lesser Antilles arc cumulate suites; Stamper et al. 2014;Cooper et al. 2016), but has been observed in similar amphibolebearing ultramafic lithologies in exposed arc root complexes (e.g., Dessimoz et al. 2012;Bouilhol et al. 2015).Textural evidence indicates that amphibole forms through a peritectic reaction consuming clinopyroxene.This is particularly evident in the Type 3 hornblende-gabbro AAN-007 where clinopyroxene is replaced by amphibole along rims, cracks and cleavage planes (Fig. 9a), similar to textures observed by Smith (2014).Plagioclase saturation is delayed and it first appears as a major crystallising phase in the hornblendegabbros where An 77-80 plagioclase is in apparent equilibrium with relatively evolved, aluminous clinopyroxene (Mg# 75-80, up to 9 wt% Al 2 O 3 ) and amphibole (Mg# 73-76, up to 15 wt% Al 2 O 3 ).The suppression of plagioclase crystallisation is recorded by a marked increase in clinopyroxene Al 2 O 3 content from ~ 3 wt% at Mg# 88 to 9 wt% Al 2 O 3 at Mg# 75 (Fig. 8).Elevated Al 2 O 3 contents reflect an increase in Al 2 O 3 /SiO 2 of the residual melt and increased Ca-Tschermak's substitution in clinopyroxene in the absence of plagioclase fractionation (e.g., Müntener et al. 2001;Villiger et al. 2007).Upon plagioclase saturation, a decrease in clinopyroxene Al 2 O 3 content would be expected.Clinopyroxene crystallisation, however, stops before the appearance of plagioclase and hence the plagioclase-in point is not recorded in clinopyroxene compositions.Amphibole crystallising together with plagioclase in Type 3 cumulates records near-constant Al 2 O 3 contents over a range of Mg# from 76 to < 70 (Fig. 9).The two key characteristics of the Type 2 and 3 cumulates are thus the early crystallisation of amphibole and the delayed appearance of plagioclase.The stability of amphibole and plagioclase in arc magmas is mainly a function of melt H 2 O content and pressure.Experimental studies addressing the differentiation of hydrous arc magmas at pressures ranging from 0.2 to 1.3 GPa indicate that the stability of plagioclase is suppressed with increasing pressure and water content (Fig. 10; e.g., Sisson and Grove 1993;Müntener et al. 2001;Müntener and Ulmer 2006;Pichavant and Table 2 Overview of equilibrium mineral assemblages in the Nisyros samples compared with experimental studies at appropriate conditions that reproduce these assemblages  et al. (2015) to include data from additional experimental studies (Müntener et al. 2001;Pichavant and Macdonald 2007;Blatter et al. 2013;Nandedkar et al. 2014;Andújar et al. 2015).The Nisyros Type 2 and 3 cumulates evolve from (hornblende-)wehrlites to hornblendegabbros and plagioclase-hornblendites, consistent with crystallisation at 0.5-0.8GPa (blue arrow).See text for further discussion Page 17 of 23 95 Macdonald 2007; Nandedkar et al. 2014;Melekhova et al. 2015).In addition, these studies conclude that high pressure and water content favour the crystallisation of amphibole at higher temperature (earlier in the crystallisation sequence) and potentially before plagioclase.Specifically, crystallisation experiments of primitive arc magmas at a pressure of 0.7 GPa and initial water content ≥ 3 wt% (Blatter et al. 2013;Nandedkar et al. 2014;Melekhova et al. 2015) reproduce the observed range of Type 2 and 3 cumulate lithologies (Fig. 10) as well as diagnostic compositional aspects such as the high Al 2 O 3 content of clinopyroxene and high Mg# and Al VI content of amphibole (Fig. 8).These studies find (hornblende-)wehrlite crystallisation at 1150-1050 °C while hornblende-gabbros and plagioclase-hornblendites progressively form from 1050 to 950 °C (Fig. 10).Although the total Al 2 O 3 content of amphibole is mainly dependent on melt composition and temperature, Al VI content has a significant pressure dependency (e.g., Holloway and Burnham 1972;Putirka 2016).The combination of high Al VI and Mg# in amphibole is diagnostic of crystallisation at lower crustal pressures (0.7 GPa).Crystallisation experiments undertaken at lower pressure (0.2-0.4 GPa; Pichavant et al. 2002;Pichavant and Macdonald 2007;Andújar et al. 2015) invariable yield amphibole saturation at lower temperatures and after plagioclase, producing amphibole with lower Mg# (< 70) and Al VI contents < 0.3 apfu (Fig. 8).
On the basis of the convincing overlap with experimental assemblages, we conclude that the Nisyros Type 2 and 3 cumulates constitute a high-pressure (0.5-0.8 GPa) hydrous crystallisation trend.Such a trend is predicted by experimental studies (e.g., Nandedkar et al. 2014;Melekhova et al. 2015) and reported for exposed arc root complexes (e.g., Dessimoz et al. 2012;Bouilhol et al. 2015), but has not been previously found in arc cumulate suites.A pressure of 0.7 GPa equates to a depth of 25-30 km, roughly the crustal thickness underneath Nisyros (e.g., Tirel et al. 2004).As such, the Type 2 and 3 cumulates provide a rare glimpse into a major magma differentiation level at the base of the crust that is not preserved in the erupted phenocrysts.

Water-saturated shallow crystallisation
The Type 1 cumulates are characterised by the association of calcic plagioclase (> An 85 ) with relatively evolved clinopyroxene (Mg# 73-86, < 1000 ppm Cr) and olivine (Fo 70-82 , < 800 ppm Ni; Fig. 8).A similar assemblage is found in cumulates from the Lesser Antilles Arc where it is generally interpreted as the result of a two-stage crystallisation history (e.g., Arculus and Wills 1980;Tollan et al. 2012;Stamper et al. 2014;Cooper et al. 2016).The calcic composition of plagioclase reflects the suppression of plagioclase stability under hydrous conditions (e.g., Sisson and Grove 1993;Müntener et al. 2001;Pichavant and Macdonald 2007;Melekhova et al. 2015), whereas the abundance of plagioclase and relatively evolved composition of olivine and clinopyroxene in these cumulates suggests that they have crystallised from already differentiated, lower Mg# mafic melts.In line with the interpretation of the Antilles cumulates, we therefore conclude that the Type 1 cumulates are the likely crystallisation products of hydrous, mildly evolved basaltic andesitic melts that underwent prior differentiation.It is probable that the mildly evolved character of these melts results from fractionation of the Type 2 and 3 cumulates in the lower arc crust.Upon ascent to shallow levels in the crust, the mildly evolved melts reach water saturation, thereby promoting the crystallisation of calcic plagioclase (e.g., Arculus and Wills 1980).This interpretation is supported by crystallisation experiments at 0.2 GPa using a mildly evolved (Mg# 60) starting compositions under H 2 O-saturated conditions (Sisson and Grove 1993) that reproduce the co-occurrence of high-An plagioclase and evolved olivine an clinopyroxene (grey plusses in Fig. 8).In addition, it agrees with the suggestion of Braschi et al. (2012Braschi et al. ( , 2014) ) that mafic melts pond at the base of a (rhyo) dacitic mush reservoir in the upper crust, leading to crystallisation of Type 1 cumulate assemblages that can subsequently be entrained in HPRD host melts or enclaves.Textural relationships indicate that amphibole is a late crystallising phase in Type 1 cumulates and typically associated with the development of lower An rims (~ An 60 ) in response to interaction with (host) melts, which is discussed in more detail below.

Shallow level crystallisation of LPA suite phenocrysts
There is a clear dichotomy between the early and late appearance of plagioclase in the LPA and HPRD suite, respectively.In contrast to the delayed plagioclase saturation seen in the Type 2 and 3 cumulates, textural and compositional evidence indicates that plagioclase is an early crystallising phase in the LPA suite.Plagioclase (up to An 90 ) occurs as phenocrysts together with primitive clinopyroxene (Mg# > 85), which has Eu/Eu* < 1 in all but the most primitive crystals (see supplementary material) and low Al 2 O 3 contents.Clinopyroxene in LPA basaltic andesites shows increasing Al 2 O 3 content down to Mg# ~ 86.At this point, an inflection in the Al 2 O 3 versus Mg# trend marks the appearance of plagioclase as a crystallising phase (Fig. 8).Low-Al clinopyroxene is reproduced by fractional crystallisation experiments of high-Mg# andesitic melts at upper crustal pressures (0.2 GPa; Grove et al. 2003), whereas experiments at higher pressures (0.4-0.7 GPa) invariably produce more aluminous clinopyroxene (Fig. 8; Müntener et al. 2001;Pichavant and Macdonald 2007;Melekhova et al. 2015).
In contrast to the Type 1 cumulates, LPA suite phenocryst assemblages are characterised by lower plagioclase An% at a given olivine composition and lower CaO contents in clinopyroxene (Fig. 8).These characteristics likely result from a combination of lower pressure, and possibly higher temperature and lower water contents, compared to the cumulate assemblages in the HPRD suite.A comparison with experimentally produced phenocryst assemblages suggests that the LPA phenocrysts crystallised from hydrous melts (2-4 wt% H 2 O) in the upper crust (0.2-0.4 GPa).Anhydrous crystallisation experiments produce more Na-rich plagioclase (Fig. 8; Almeev et al. 2013) which is at odds with the observed plagioclase compositions and the appearance of amphibole in the LPA dacites.High whole rock and clinopyroxene phenocryst Mg# (up to 68 and 91, respectively) suggest that the LPA basaltic andesites represent relatively primitive melts.The absence of systematic trace element or isotopic differences compared to mafic enclaves in the HPRD suite (Klaver et al. 2016) effectively rules out that the LPA suite derives from compositionally distinct or less hydrous primary melts than those from which the Type 2 and 3 cumulates crystallised in the lower crust.Hence, we propose that the dichotomy in mineral assemblages and compositions between the LPA and HPRD suites predominantly reflects depth of crystallisation.Whereas the HPRD suite records crystallisation at the base of the crust, the LPA suite appears to have differentiated from relatively primitive (high Mg#) melts that were injected into mid-to upper crustal reservoirs without pronounced modification at the base of the crust.

A deep or shallow amphibole sponge?
Whole rock data of global arc lavas suggest that amphibole is a common, but often cryptic, phase in the differentiation of arc magmas (e.g., Davidson et al. 2007).Although studies of arc root complexes (e.g., Dessimoz et al. 2012;Bouilhol et al. 2015) and arc cumulates (e.g., Smith 2014;Cooper et al. 2016) have established that amphibole forms through a peritectic reaction consuming clinopyroxene, it is still unclear at which level in the volcanic plumbing system this "amphibole sponge" exerts the strongest control on melt compositions.Nisyros amphibole can be roughly divided into three groups based on its composition: (1) amphibole phenocrysts in the LPA dacites, (2) late-stage, interstitial amphibole in Type 1 cumulates and (3) peritectic amphibole in Type 2 and 3 cumulates.Phenocrysts in the LPA dacites match amphibole crystallising from andesitic to dacitic liquids at 0.2 GPa and > 6 wt% H 2 O in the melt (Andújar et al. 2016) and are in equilibrium with their host melt.Therefore, they do not represent reaction-replacement amphibole and are not considered further here.
Textural evidence suggests that amphibole has formed through the reaction of clinopyroxene with primitive melts in Type 2 and 3 cumulates, which is supported by high Mg#, low trace element contents and a lack of Eu-anomalies in these amphiboles (Fig. 6).Amphibole in these cumulates displays a progressive crystallisation trend with slightly increasing tetrahedral Al (Al IV ) and TiO 2 contents and decreasing Mg# from Type 2 to Type 3 cumulates.A similar trend is present in core-to-rim traverses of euhedral amphibole in hornblende-gabbro 11).Entrainment of this cumulate nodule in a rhyolitic melt caused the Fig. 11 Diagram of Al IV versus Mg# for amphibole in the cumulates.Amphibole in the Type 2 hornblende-wehrlite (NCX-06) displays a well-defined linear trend along a Mg-Tschermak's substitution vector that results from subsolidus re-equilibration during cooling (Ernst and Liu 1998;Dessimoz et al. 2012).The Type 3 cumulates display a small core-to-rim increase in amphibole Al IV with decreasing Mg# that is consistent with progressive crystallisation.Reaction with a felsic melt leads to resorption, decreasing Al IV and a small increase in Mg#.See text for further discussion Page 19 of 23 95 pronounced disequilibrium in the amphibole crystals, as indicated by a marked decrease in Mg#, TiO 2 , Na 2 O and Al IV contents in the resorbed rims (Fig. 9).Interestingly, Mg# increases again in the outermost rims of the resorbed crystals.The resorbed rims of AAN-007 have similar low Al IV content and Mg# as interstitial amphibole in Type 1 cumulates (Fig. 11).The latter has an enriched trace element composition, including negative Eu-anomalies and high La/ Yb and Ba contents (Fig. 6), as also found by Smith (2014) and Cooper et al. (2016) for interstitial amphibole in similar arc cumulate nodules.Such amphibole compositions have not been reproduced in crystallisation experiments and hence suggest growth in response to interaction with an infiltrating melt rather than equilibrium crystallisation.A common process of reaction with an evolved, hydrous melt is, therefore, evoked for the crystallisation of the low Al IV amphiboles in Type 1 cumulates and the resorbed rims of AAN-007 (Smith 2014).In addition, the low Al VI contents of interstitial amphibole (Fig. 8) corroborate our conclusion that the Type 1 cumulates crystallised at shallow crustal levels.
The peritectic reaction in which amphibole forms at the expense of clinopyroxene is instrumental in driving derivative melts to andesitic-rhyolitic calc-alkaline compositions over a relatively short crystallisation interval (e.g., Holloway and Burnham 1972;Cawthorn and O'Hara 1976;Davidson et al. 2007).We find evidence for the crystallisation of amphibole at different crustal levels and under different conditions.Although the infiltration of a felsic melt in Type 1 cumulates at shallow levels can lead to the crystallisation of interstitial amphibole, we propose that this process is volumetrically insignificant and cannot account for the important role of amphibole in generating calc-alkaline differentiation trends.The Type 2 and 3 cumulates, on the other hand, record early stabilisation of amphibole in equilibrium with mafic melts and experimental studies demonstrate that amphibole will continue to crystallise over a large temperature interval (e.g., Nandedkar et al. 2014).A detailed assessment of the role of amphibole in generating the (rhyo)dacites on Nisyros is outside the scope of this work.We do note, however, the striking resemblance in lithology between the Type 3 cumulates and the hornblende-gabbroic assemblage found to be in equilibrium with Mount St. Helens dacites (Blatter et al. 2017).This study demonstrates that porphyritic dacites, compositionally and texturally similar to the Nisyros HPRD suite, are multiply saturated in plagioclase, pyroxene, amphibole and magnetite at 0.7 GPa and form along a peritectic reaction involving amphibole and clinopyroxene.Hence, we tentatively propose that the lower crustal Type 3 cumulates can be complementary to the Nisyros HPRD suite magmas.

Synthesis: structure of Nisyros' plumbing system
Our interpretation of the structure of Nisyros' volcanic plumbing system based on the petrology and composition of cumulates and phenocrysts is shown schematically in Fig. 12.This model is in general agreement with ideas proposed by Wyers and Barton (1989) and Seymour and Vlassopoulos (1992), but provides further refinement and, for the first time, documents the importance of differentiation at the base of the arc crust.We interpret the Type 2 cumulates with Al-rich clinopyroxene as the most primitive cumulate assemblage found in the HPRD suite.Derivative melts formed by the fractionation of Type 2 cumulates will have lower Mg# and higher H 2 O contents than primitive melts and will, upon ascent and water saturation, crystallise the common Type 1 cumulates.In the lower crust, Type 2 cumulates will mature to Type 3 hornblende-gabbros and plagioclase-hornblendites along a peritectic reaction in which clinopyroxene is consumed and amphibole and dacitic melt are generated (e.g., Blatter et al. 2017).The latter is a prime candidate to feed a shallow (rhyo)dacitic mush system that, upon reactivation by mafic melts, can be erupted as porphyritic, enclave-bearing domes and lava flows of the HPRD suite (this study; Braschi et al. 2012Braschi et al. , 2014;;Zouzias and St Seymour 2014).As such, we infer that differentiation of the HPRD suite is predominantly controlled by cumulate-melt reaction processes in the lower crust, but that most of their textural features are obtained through shallow crystallisation and mixing processes.
The LPA suite, on the other hand, does not show evidence for protracted differentiation in the lower crust.Textural and chemical evidence indicates that plagioclase is an early crystallising phase and in equilibrium with very primitive (Mg# up to 91) clinopyroxene.Thus, the LPA (basaltic) andesites appear to have differentiated from primitive melts predominantly at shallow crustal levels.Given the lack of hybridisation between the LPA and HPRD suites, we argue that the two suites evolved in isolated parts of the volcanic plumbing system.Although there is a gradual change from LPA-dominated to HPRD-dominated magmatism with time, LPA andesites and basaltic andesites have also been emplaced during the caldera-forming stage (Fig. 1) without any associated major changes in HPRD suite petrography.In addition, the eccentric activity stage is characterised by rapid and repeated switching between eruption of HPRD dacites and LPA basaltic andesites as seen in unit PA (Lies tuff cone; Fig. 1).These lines of evidence contradict the assertion that the LPA and HPRD suites can be accommodated in a single shallow reservoir through gradual differentiation with time, as for instance proposed by Francalanci et al. (1995).

Implications for island arc magmatism
One of the implications of this study is that the HPRD (rhyo)dacites could be generated through a peritectic reaction at the base of the Nisyros arc crust as recorded in the Type 3 cumulates.In contrast, they acquire their petrographic features through crystallisation and mixing in a mush system in the mid-to upper crust.Hence, Nisyros is in good agreement with the notion of Annen et al. (2006) that significant decoupling between the geochemical and petrographic characteristics of arc magmas is probably the norm rather than the exception in volcanic rocks in arcs built on (thinned) continental crust.Whereas differentiation in the lower arc crust can exert a dominant control on the geochemical composition of the magmas (e.g., Hildreth and Moorbath 1988;Davidson et al. 2007), textural and petrographic aspects are predominantly obtained during shallow storage (e.g., Eichelberger 1978;Annen et al. 2006;Bachmann and Huber 2016).This provides an elegant solution to the common misfit between observed phenocryst phases and cumulate assemblages in arc lavas and the fractionating mineral assemblages required to reproduce whole rock major element variations in arc suites (e.g., Annen et al. 2006;Melekhova et al. 2015).Type 1 cumulates, the products of shallow crystallisation of evolved melts at water saturation, dominate in arc cumulate suites (e.g., Arculus and Wills 1980;Tollan et al. 2012;Cooper et al. 2016).Together with phenocryst assemblages, such cumulates rarely reflect the early stages of magmatic differentiation in the lower arc crust (e.g., Arculus and Wills 1980;Tollan et al. 2012).Moreover, many arc suites have a geochemical signature of phases that are stable at pressures pertinent to the lower crust but are scarcely present as phenocryst phases, such as amphibole (e.g., Davidson et al. 2007) and garnet (e.g., Macpherson et al. 2006;Alonso-Perez et al. 2009).The Nisyros Type 2 and 3 cumulates, however, provide a rare insight into differentiation of primitive, hydrous melts Fig. 12 Schematic cross-section illustrating the relationship between the different cumulate and phenocryst suites found on Nisyros.The HPRD suite preserves evidence for differentiation at the base of the arc crust where the reaction of wehrlitic (Type 2) cumulates with hydrous melts stabilises amphibole to form the Type 3 cumulates.Type 1 cumulates, however, crystallise from already differentiated melts at shallower levels in the plumbing system, likely at the base of a shallow mush system where the HPRD host lavas acquire their textural features.In contrast, the LPA suite does not show evidence for protracted modification at the base of the arc crust and differentiated from high-Mg# melts at shallow levels.There is no evidence for interaction between the two suites in their crystal records Page 21 of 23 95 at lower crustal conditions in arcs.The replacement of clinopyroxene by amphibole in these cumulates provides direct evidence for the fundamental importance of amphibole fractionation in generating evolved melts (e.g., Holloway and Burnham 1972;Cawthorn and O'Hara 1976), consistent with textural evidence from arc root complexes (e.g., Dessimoz et al. 2012;Bouilhol et al. 2015).Whether amphibole is actually present as a phenocryst phase in these evolved rocks depends on the P-T-pH 2 O conditions in the shallow storage region and does not reflect differentiation in a deep crustal hot zone.

Conclusions
This study reports petrographic descriptions and mineral major and trace element data of lavas, pyroclastic rocks, enclaves and cumulates from Nisyros volcano, Greece.We propose a new division of the Nisyros volcanic rocks into a low-porphyricity andesite (LPA) and high-porphyricity (rhyo)dacite (HPRD) suite, which have evolved along distinct petrogenetic pathways and show no evidence for hybridisation.The preservation of primitive, lower crustal cumulate fragments in the HPRD suite provide a glimpse of the earliest stages of magma differentiation in arcs.The main conclusions that can be drawn from this study are: 1. Type 2 and 3 cumulates define a continuum from (hornblende-)wehrlites to plagioclase-hornblendites that form a striking resemblance to the products of fractionation experiments of hydrous primitive melts at lower crustal (0.7 GPa) pressures and to exposed arc root complexes.2. Clinopyroxene Al 2 O 3 contents can be used in tandem with textural relationships to infer the timing of plagioclase saturation: clinopyroxene phenocrysts in the LPA suite have low Al 2 O 3 contents, consistent with early plagioclase saturation at low pressure, whereas clinopyroxene in the Type 2 and 3 cumulates shows Al 2 O 3 contents up to 9 wt% due to the suppression of plagioclase crystallisation.3. Reaction of the wehrlitic Type 2 cumulates with a primitive, hydrous melt leads to the formation of amphibole at the expense of clinopyroxene to form Type 3 cumulates, a key reaction in producing (rhyo)dacitic derivative melts.The major and trace element composition of amphibole in these cumulates confirms crystallisation at high pressure (~ 0.7 GPa), in contrast to late-stage interstitial amphibole in Type 1 cumulates.4. The most common cumulates in the HPRD suite are plagioclase-clinopyroxene-olivine Type 1 cumulates.Anorthite-rich plagioclase in combination with low-Mg# clinopyroxene and olivine attest to crystallisation from an already mildly evolved melt at upper crustal pressures and water saturated conditions.Hence, the Type 1 cumulates are not complementary to the (rhyo)dacites in which they are hosted, but rather an accumulation of phenocrysts from a mildly evolved melt in a shallow reservoir.

Fig. 2
Fig. 2 Petrographic features of the volcanic rocks of Nisyros.Photomicrographs are in plain-polarised light unless otherwise stated; a LPA suite Lava Flow 8 (unit LF8) and associated scoria deposits (XO) that were deposited during the eccentric activity stage at Aghia Irini; backpack (~ 60 cm) for scale; b typical texture of a low-porphyricity andesite (sample AAN-029-LF6), in cross-polarised light (XPL): very fine microlite to glassy groundmass and low crystal content with sub-to euhedral plagioclase, olivine and clinopyroxene (not in this photo) phenocrysts.Note the faint flow banding and alignment of plagioclase phenocrysts; c glomeroporphyritic aggregate of euhedral, unzoned hornblende, plagioclase, orthopyroxene and Fe-Ti-oxides in a very fine grained groundmass: the typical texture of LPA dacites (AAN-031-LF1); d complexly zoned plagioclase crystal in an LPA andesite; spots show An% (AAN-015-LF3);

Fig. 5
Fig.5Overview of the composition (An% for plagioclase, Mg# for olivine, amphibole, clino-and orthopyroxene) of major phases in the magmatic products of Nisyros.The samples have been grouped based on petrography and field occurrence in a lowporphyricity andesite (LPA) and high-porphyricity (rhyo)dacite (HPRD) suite and are sorted by increasing SiO 2 content downward for the LPA and HPRD suite lavas and enclaves, and median clinopyroxene Mg# for the cumulates.Abbreviations for the volcanic units afterVolentik et al. (2005b)

Fig. 6
Fig. 6 Amphibole major and trace element composition.The samples have been grouped based on petrography and field occurrence in a low-porphyricity andesite (LPA) and high-porphyricity (rhyo)dacite (HPRD) suite.A typical example of a bimodal amphibole population in intermediate arc rocks (Kiss et al. 2014) is shown for comparison in a and b.Nisyros amphibole lacks such a bimodal distribution

Fig. 7
Fig. 7 Clinopyroxene TiO 2 content versus Mg# showing the difference between the LPA basaltic andesites (LPA-BA) and andesites (LPA-A) and the HPRD suite.Andesites in the LPA suite have higher TiO 2 contents than the other groups and are not consistent with crystallisation from a hybrid melt.See supplementary material for an extended diagram showing individual samples

Fig. 8
Fig. 8 Composition of Nisyros lava and cumulate mineral assemblages, compared with published results of crystallisation experiments: a olivine Fo% versus plagioclase An%; b amphibole Mg# versus plagioclase An%; c clinopyroxene Al 2 O 3 content versus Mg#; d clinopyroxene CaO content versus Mg#; e amphibole Al VI versus Mg#.Symbols for the Nisyros samples depict the median composition; error bars indicate twice the interquartile range.In a and b, the "< 0.5% plagioclase" field at the right depicts (near-)plagioclase-free

Table 1
Overview of equilibrium mineral assemblages that can be recognised in the LPA and HPRD suites Mineral abbreviations: Pl plagioclase, Opx orthopyroxene, Cpx clinopyroxene, Ox Fe-Ti-oxides, Hbl hornblende, Ol olivine, n/p not present