Tectonomagmatic Constraints on the Sources of Eastern Mediterranean K-rich Lavas

Potassium-enriched mafic lavas, and in particular basalts of the high-K series, are important end-members of subduction zone volcanism. Two petrogenetic models can explain the generation of mafic K-rich lavas, involving either melting of ancient, enriched lithospheric mantle sources (single-stage model) or melting triggered by recent refertilization by subduction-related components derived from subducted sediments or oceanic crust (multi-stage model). These two models are tested for post-collisional high-K rocks of Eocene-Oligocene age from the Rhodopes, SE Bulgaria, based on new major element, trace element and Sr-Nd-Hf-Pb isotope data. The single-stage model is evaluated by Sr-Nd isotope modelling assuming the presence of ancient enriched lithospheric domains whereas the multi-stage model is assessed by comparing the compositions of the Bulgarian lavas with those of lavas from Santorini. Santorini lavas are considered to sample the current trace element and isotope inventory of the long-lived Aegean subduction zone system. This northward facing system has been active most probably since late Jurassic or early Cretaceous times and was potentially involved in refertilizing the mantle sources of the Bulgarian lavas. In addition to the data for Bulgarian lavas, we present new major element, trace element and Sr-Nd-Hf-Pb isotope data for Santorini. These are evaluated together with previously published data to infer the mode of source enrichment in the Aegean realm. The Bulgarian lavas exhibit a broad range of compositions from medium-K to high-K and shoshonitic with radiogenic Sr-87/Sr-86 (0.706-0.709), and unradiogenic epsilon Nd (-5.7 to -1.9) and epsilon Hf (-3 to +3) isotope signatures. The trace element budget was apparently well buffered against shallow-level crustal assimilation as documented by the trace element and Sr-Nd isotope systematics. Modelling of the Sr-Nd isotope compositions of the Bulgarian lavas argues for recent (Mesozoic to Cenozoic) source enrichment. Therefore single-stage models involving melting of ancient (>1 Ga) lithospheric mantle can be confidently ruled out, in agreement with tectonic models for the region. The enriched isotope signatures, together with a pronounced enrichment of incompatible elements, rather indicate mantle refertilization by a subduction component similar to continent-derived sediments subducted at the Hellenic Trench at present. The Bulgarian lavas record a predominant influx of subduction fluid components in comparison with the Santorini lavas. Collectively, the data presented for the Bulgarian lavas are clearly in favour of a multi-stage model.