Publications

2021
Koutsovitis P, Soukis K, Voudouris P, Lozios S, Ntaflos T, Stouraiti C, Koukouzas N. {The Late Cretaceous magmatic arc of the south Aegean: Geodynamic implications from petrological and geochemical studies of granitoids from Anafi island (Cyclades–Greece)}. International Geology Review [Internet]. 2021;00:1–24. WebsiteAbstract
In the Aegean region (Cyclades–Greece), a large part of the Island of Anafi consists of exhumed high-grade metamorphic sequences (amphibolites, serpentinites and metasediments) that have been intruded by Late Cretaceous intermediate and felsic granitoids. These correspond to I-type arc-related rocks, displaying calc-alkaline geochemical affinities. Variability in their petrographic features, mineral composition and geochemistry is assigned to differentiation processes that mostly involved plagioclase and/or K-feldspar, and to a lesser extent amphibole and biotite. Mineral chemistry and geochemical results suggest that the Anafi granitoids are highly comparable with the Late Cretaceous granitoid intrusive rocks from East Crete and Donousa island. The amphibole and zircon saturation thermometry yields relatively moderate temperature crystallization conditions ($\sim$790°C) for the intermediate granitoids and lower for the felsic granitoids ($\sim$630°C). Geobarometric calculations based upon the chemistry of magnesiohornblende, as well as the normative (Qz), (Ab) and (Or) contents clearly point to shallow intrusion conditions ($\sim$2.0–6.5 kbar), corresponding to an estimated depth of crystallization of $\sim$12 ± 4 km. The thrust sheets that overly the Palaeogene flysch in Anafi (metasediments with serpentinized peridotites, amphibolites and metabasites), constitute a subducted and metamorphosed oceanic sequence. These metamorphic units likely represent a part of the Pindos–Cycladic Blueschist Unit domain that was subducted at an earlier pre-Campanian stage. In the hydrated mantle wedge, incorporation of shallow level granitoids within metamorphic units was likely facilitated via corner flow intrusion mechanisms. Following the intrusion, the granitoid rocks were exhumed in a syn-convergent setting, along with the metamorphic thrust sheets, by continuous underplating of more external units, thus, escaping penetrative ductile deformation.
Ross J, Voudouris P, Melfos V, Vaxevanopoulos M, Soukis K, Merigot K. {The Lavrion silver district: Reassessing its ancient mining history}. Geoarchaeology. 2021:1–26.Abstract
Lavrion was an important silver mining district in the ancient world, exploited almost continuously from the fourth millennium BC. Its mining history is central to understanding the availability of silver, lead and copper in Greece and the Aegean through the Bronze age. This history is also relevant to the phenomenal rise of Athens in the 5th century BC. Our reassessment of the mineralisation and ancient mining history at Lavrion provides a clearer understanding of its capacity to generate wealth and the scale of its impact on Athens. Integration of new geochemical and geological data with existing information has produced outcomes that help address these two issues. Foremost is the contrasting character of mineralisation at the upper 'first' and lower 'third' contacts and their marked differences in silver content, spatial distribution and supergene alteration. Discovery of the concealed, bonanza mineralisation of the third contact at Kamariza, possibly early in the 5th century BC, followed about 3000 years of mining, which appears to have been largely restricted to discontinuous, low-grade mineralisation at the first contact. Undoubtedly, this later discovery had enduring impacts on Lavrion and Athens, and it most likely funded the trireme fleet, which brought victory over the Persians at Salamis in 480 BC.
Voudouris P, Melfos V, Mavrogonatos C, Photiades A, Moraiti E, Rieck B, Kolitsch U, Tarantola A, Scheffer C, Morin D, et al. {The Lavrion Mines : A Unique Site of Geological and}. Minerals. 2021;11:76.Abstract
The Lavrion area corresponds to the western part of the Attic-Cycladic metamorphic belt, in the back-arc region of the active Hellenic subduction zone. Between the Eocene and the Miocene, metamorphic rocks (mainly marbles and schists) underwent several stages of metamorphism and deformation due to collision and collapse of the Cycladic belt. Exhumation during the Miocene was accommodated by the movement of a large-scale detachment fault system, which also enhanced emplacement of magmatic rocks, leading to the formation of the famous Lavrion silver deposits. The area around the mines shows the stacking of nappes, with ore deposition mainly localized within the marbles, at marble-schist contacts, below, within, or above the detachment. The Lavrion deposit comprises five genetically-related but different styles of mineralization, a feature never observed in another ore deposit elsewhere, containing the highest number of different elements of any known mining district. The local geology, tectonic, and magmatic activity were fundamental factors in determining how and when the mineralization formed. Other key factors, such as the rise and the fall of sea level, which resulted from climate change over the last million years, were also of major importance for the subsequent surface oxidation at Lavrion that created an unmatched diversity of secondary minerals. As a result, the Lavrion deposit contains 638 minerals of which Lavrion is type-locality for 23 of them, which is nearly 12% of all known species. Apart from being famous for its silver exploitation, this mining district contains more minerals than any other district on Earth.
2020
Coleman MJ, Schneider DA, Grasemann B, Soukis K, Lozios S, Hollinetz MS. {Lateral Termination of a Cycladic-Style Detachment System (Hymittos, Greece)}. Tectonics. 2020;39.Abstract
The bedrock of Hymittos, Attic peninsula, Greece, exposes a pair of low-angle crustal-scale ductile-then-brittle detachment faults. The uppermost detachment fault separates sub-greenschist facies phyllite and marble of a Pelagonian Zone hanging wall, from greenschist facies metasedimentary schist, calc-schist, and marble correlated to the Cycladic Blueschist Unit. A second, structurally lower detachment fault subdivides the metamorphic rocks of the Cycladic blueschist unit footwall into middle and lower units. There is a marked step in metamorphic grade between the sub-greenschist facies uppermost package, and the middle-to-upper greenschist facies middle and lower packages. A suite of new white mica 40Ar/39Ar and zircon (U-Th)/He dates indicates accommodation of deformation along these faults occurred from the late Oligocene to the late Miocene with both faults active during the middle Miocene. The structures have clear top-S/SSW kinematics determined from flanking folds, sigmoids, shear bands, stair-stepping of strain shadows on porphyroclasts, and SCC' fabrics. The ductile-to-brittle deformation of the structures, morphology of the massif, and the increase in metamorphic grade suggest these low-angle structures are part of a major, crustal-scale extensional complex, located at the northwest end of the West Cycladic Detachment System, that accommodated Miocene bivergent exhumation of Attic-Cycladic metamorphic core complexes in the central Aegean. Taken together, the above data suggest that multiple coeval detachment branches may form in areas with high strain gradients to accommodate the mechanically necessary termination of Cycladic-style detachment systems.
Koutsovitis P, Magganas A, Ntaflos T, Koukouzas N, Rassios AE, Soukis K. {Petrogenetic constraints on the origin and formation of the Hellenic Triassic rift-related lavas}. Lithos. 2020;368-369.Abstract
Triassic rift-related volcanic rocks outcrop over all mainland Greece, comprising of trachybasalts and basaltic trachyandesites. Relatively immobile to the effects of alteration processes major and trace element abundances classify the volcanics into OIB and E-MORB lavas. They have mainly been distinguished based upon their: i) LREE contents, ii) silica-saturation index (S·I.), iii) Zr/Nb and Nb/Y ratio values; iv) Th, U, and Ta contents v) geotectonic discrimination diagrams. Their geochemistry indicates that most rocks were affected by moderate to extensive differentiation processes, mostly expressed by clinopyroxene fractionation. Some of the OIB and E-MORB volcanics are considered as primitive undersaturated, displaying (low SiO2, Zr/Nb and S.I. values, enhanced CaO/Al2O3 ratios). Calculated average mantle potential temperatures are comparable (1410 °C OIB; 1370 °C E-MORB), with melt fractions estimated at 3–5% for primary OIB magmas and 6–8% for primary E-MORB magmas. An asthenospheric origin is inferred for the OIB lavas, with melting in the garnet stability field (75–95 km; 2.5–3.0 GPa), whereas E-MORB parent magmas were formed with melting in the garnet/spinel (transitional) stability field (55–70 km; 1.8–2.2 GPa). The Hellenic Triassic rift-related lavas were most likely generated and erupted after lithospheric attenuation and extension, followed by subsequent asthenospheric upwelling of the mantle. The high calculated partial melting degrees and the observed thick and voluminous lava formations account for fast-spreading of the Tethys ridge during the Triassic. Temperature results suggest that the Hellenic Triassic rift-related magmas were generated from mantle at ambient temperature, render a mantle plume-based scenario improbable.
Christina Stouraiti, Vassiliki Angelatou SPKS, and Eliopoulos D. Effect of Mineralogy on the Beneficiation of REE from Heavy Mineral Sands: The Case of Nea Peramos, Kavala, Northern Greece. Minerals. 2020;10:387.
2019
Coleman M, Dubosq R, Schneider DA, Grasemann B, Soukis K. {Along-strike consistency of an extensional detachment system, West Cyclades, Greece}. Terra Nova. 2019;31:220–233.Abstract
The Lavrion district of the Attica Peninsula, Greece, exposes the West Cycladic Detachment System (WCDS), a low-angle crustal-scale extensional fault system separating tectonostratigraphic units of the Cycladic Blueschist Unit. New multiple single-grain fusion and step-heated white mica 40Ar/39Ar ages integrated with existing (U–Th)/He ages and independent paleothermometry resolves a syn- to post-orogenic deformation history. A structurally higher unit records Oligocene greenschist facies deformation that evolved into brittle conditions by the middle Miocene, and shares a similar history to Serifos at the southern end of the detachment system. The structurally lowest unit remained ductile until the late Miocene, preserving pervasive post-orogenic structures, similar to along-strike structures at the centre, deepest part of the fault. The similarities of structural styles and the timing of deformation across > 150 km of strike length of the detachment system indicates tens of kilometers of offset and extraordinary potential for correlating observations along Cycladic-style detachment systems.
Flansburg ME, Stockli DF, Poulaki EM, Soukis K. {Tectono-magmatic and Stratigraphic Evolution of the Cycladic Basement, Ios Island, Greece}. Tectonics. 2019;38:2291–2316.Abstract
The Cycladic Basement (CB) and the overlying Cycladic Blueschist Unit (CBU) are part of the Paleogene Cycladic subduction complex exposed in Miocene metamorphic core complexes in the distended back-arc of the retreating Hellenic subduction zone of the southern Aegean. While the Cenozoic tectono-metamorphic evolutions of the CB and the CBU have been the foci of numerous studies, this study presents new laser ablation inductively coupled plasma mass spectrometry bedrock and detrital zircon (DZ) U-Pb ages that place robust constraints on the presubduction tectonic, magmatic, and paleogeographic evolution of the CB. Zircon U-Pb ages of crystalline CB are $\sim$306-330 Ma, demonstrating local plutonism associated with regional voluminous, protracted Carboniferous magmatism related to Paleo-Tethys subduction. The plutons intruded the CB metasedimentary host-rock sequence, characterized by distinct Gondwanan DZ U-Pb provenance, Neoproterozoic to early Paleozoic maximum depositional ages, and synmagmatic, contact metamorphic zircon rims ($\sim$300-330 Ma). DZ U-Pb dating revealed postmagmatic Permian metasedimentary rocks ($\sim$270-295 Ma) that unconformably overlie the CB and have unimodal DZ spectra that indicate exhumation of the CB prior to Permian deposition within extensional basins, as well as mark the onset of CBU deposition prior to formation of the Pindos rift domain. These U-Pb results clarify the late Paleozoic-early Mesozoic evolution of the CB as a peri-Gondwanan terrane composed of Neoproterozoic and early Paleozoic metasedimentary rocks, intruded by voluminous Carboniferous arc magmatism, and exhumed in the Permian, prior to Triassic rifting and CBU deposition. Additionally, these data provide a chronostratigraphic framework and illuminate subduction-related juxtaposition within the CB metasedimentary sequence.
Coleman M, Dubosq R, Schneider DA, Grasemann B, Soukis K. {Along-strike consistency of an extensional detachment system, West Cyclades, Greece}. Terra Nova. 2019;31.Abstract
© 2019 John Wiley & Sons Ltd The Lavrion district of the Attica Peninsula, Greece, exposes the West Cycladic Detachment System (WCDS), a low-angle crustal-scale extensional fault system separating tectonostratigraphic units of the Cycladic Blueschist Unit. New multiple single-grain fusion and step-heated white mica 40 Ar/ 39 Ar ages integrated with existing (U–Th)/He ages and independent paleothermometry resolves a syn- to post-orogenic deformation history. A structurally higher unit records Oligocene greenschist facies deformation that evolved into brittle conditions by the middle Miocene, and shares a similar history to Serifos at the southern end of the detachment system. The structurally lowest unit remained ductile until the late Miocene, preserving pervasive post-orogenic structures, similar to along-strike structures at the centre, deepest part of the fault. The similarities of structural styles and the timing of deformation across > 150 km of strike length of the detachment system indicates tens of kilometers of offset and extraordinary potential for correlating observations along Cycladic-style detachment systems.
Stouraiti C, Soukis K, Voudouris P, Mavrogonatos C, Lozios S, Lekkas S, Beard A, Strauss H, Palles D, Baziotis I, et al. {Silver-rich sulfide mineralization in the northwestern termination of the Western Cycladic Detachment System, at Agios Ioannis Kynigos, Hymittos Mt. (Attica, Greece): A mineralogical, geochemical and stable isotope study}. Ore Geology Reviews. 2019;111.Abstract
© 2019 Elsevier B.V. Identification of silver-rich sulfide mineralization in the abandoned mines at northern Hymittos Mt (central Attica, Greece) district has made this region important due to its structural relation to the world-class Lavrion mining district (25 km SE of Hymittos Mt.). At Agios Ioannis Kynigos mine veins of oxidized Fe-Pb-Zn-Ag ore are hosted in mylonitic marbles that were predominantly affected by brittle-ductile deformation associated with the West Cycladic Detachment System. The major sulfides are galena and arsenian pyrite which develop as open-space filling (fractures and vugs), disseminations and veinlets. Galena contains abundant inclusions of Pb-As-Sb-Cu-Ag sulfosalts (tetrahedrite group minerals, bournonite, jalpaite). An unusual silver-rich galena ore i.e from 192 to 1500 ppm Ag, is recorded from the present study, with Ag hosted mainly in Ag-tetrahedrite, acanthite and jalpaite. Compositionally pyrite shows As enrichment and the highest As concentration i.e. up to 3.9 wt% is found in the core of the oscillatory zoned pyrites. Values of $δ$34S of galena are isotopically light (−11.2 to −12.2‰), which suggests a sedimentary component to the ore fluid, that is likely derived from leaching of the meta-sedimentary (calc-mica schists) wall rocks. Carbon and oxygen isotope compositions of ore-stage carbonates show values i.e $δ$13CV-PDB (‰) = −8.2 to −11.7‰ and $δ$18O = 24.3 to 25.9‰ and for the non-mineralized carbonates $δ$13CV-PDB (‰) = +0.33 to 2.52‰ and $δ$18O = 20.8 to 26.1‰. The C and O isotopic composition of the carbonate wall rock distal to the orebody is typical of marine carbonates, whereas hydrothermal carbonates (brown calcite and dolomite) are dominated by an external, light C source. The lowest $δ$13CV-PDB values are recorded from cerussite and correspond to the low temperature supergene alteration from meteoric waters. The homogenized values of $δ$18O in the mineralized carbonates indicate interaction of fluid with host carbonates under a low fluid/rock ratio. The Pb-Zn-Ag mineralization at Agios Ioannis Kynigos bears many similarities with the low-temperature carbonate replacement-style mineralization in Kamariza, Lavrion area. Ore deposition in the northern Hymittos detachment fault is focused in the structurally upper parts of the detachment fault where circulation of meteoric fluids is dominant in the hydrothermal system.
Flansburg ME, Stockli DF, Poulaki EM, Soukis K. {Tectono-magmatic and Stratigraphic Evolution of the Cycladic Basement, Ios Island, Greece}. Tectonics. 2019;38.Abstract
©2019. American Geophysical Union. All Rights Reserved. The Cycladic Basement (CB) and the overlying Cycladic Blueschist Unit (CBU) are part of the Paleogene Cycladic subduction complex exposed in Miocene metamorphic core complexes in the distended back-arc of the retreating Hellenic subduction zone of the southern Aegean. While the Cenozoic tectono-metamorphic evolutions of the CB and the CBU have been the foci of numerous studies, this study presents new laser ablation inductively coupled plasma mass spectrometry bedrock and detrital zircon (DZ) U-Pb ages that place robust constraints on the presubduction tectonic, magmatic, and paleogeographic evolution of the CB. Zircon U-Pb ages of crystalline CB are $\sim$306-330 Ma, demonstrating local plutonism associated with regional voluminous, protracted Carboniferous magmatism related to Paleo-Tethys subduction. The plutons intruded the CB metasedimentary host-rock sequence, characterized by distinct Gondwanan DZ U-Pb provenance, Neoproterozoic to early Paleozoic maximum depositional ages, and synmagmatic, contact metamorphic zircon rims ($\sim$300-330 Ma). DZ U-Pb dating revealed postmagmatic Permian metasedimentary rocks ($\sim$270-295 Ma) that unconformably overlie the CB and have unimodal DZ spectra that indicate exhumation of the CB prior to Permian deposition within extensional basins, as well as mark the onset of CBU deposition prior to formation of the Pindos rift domain. These U-Pb results clarify the late Paleozoic-early Mesozoic evolution of the CB as a peri-Gondwanan terrane composed of Neoproterozoic and early Paleozoic metasedimentary rocks, intruded by voluminous Carboniferous arc magmatism, and exhumed in the Permian, prior to Triassic rifting and CBU deposition. Additionally, these data provide a chronostratigraphic framework and illuminate subduction-related juxtaposition within the CB metasedimentary sequence.
Poulaki EM, Stockli DF, Flansburg ME, Soukis K. {Zircon U-Pb Chronostratigraphy and Provenance of the Cycladic Blueschist Unit and the Nature of the Contact With the Cycladic Basement on Sikinos and Ios Islands, Greece}. Tectonics. 2019;38.Abstract
©2019. American Geophysical Union. All Rights Reserved. Sikinos and Ios Islands, located in the Southern Cyclades, represent part of a Cenozoic metamorphic core complex system that exposes subduction-related metamorphic rocks in the highly extended back-arc region of the Hellenic subduction zone. These exhumed HP-LT metamorphic units are composed of Mesozoic metasedimentary rocks of the Cycladic Blueschist Unit (CBU) and the Paleozoic Cycladic Basement (CB). The magmatic and stratigraphic evolution of these units, as well as the nature of the contact between the CBU and CB, have remained poorly understood. We used zircon U-Pb dating to determine crystallization ages of the CB on Sikinos and the maximum deposition ages and detrital provenance of the metasedimentary units to reconstruct the Mesozoic to early Cenozoic stratigraphic and tectonic evolution of the CBU on both islands. The results reveal that the CB in Sikinos is composed of Cambrian-Silurian metasedimentary rocks intruded by Carboniferous granites and is overlain by metasedimentary rocks of the CBU with depositional ages spanning from Permo-Triassic to Late Cretaceous. The provenance data from the CBU records a long-lived tectonic evolution from Paleo-Tethys subduction and rifting, to passive margin formation, and to subduction of the Neo-Tethyan Pindos basin. The continuous stratigraphic record and provenance evolution from the CB into the CBU imply a para-autochthonous relationship. On NE Sikinos and Ios, stratigraphic constraints suggest older-over-younger relationships along cryptic-thrusts, supporting premetamorphic or synmetamorphic structural repetition of the CBU by imbrication, likely during subduction underplating.
Voudouris P, Mavrogonatos C, Melfos V, Spry PG, Magganas A, Alfieris D, Soukis K, Tarantola A, Periferakis A, Kołodziejczyk J, et al. {The geology and mineralogy of the Stypsi porphyry Cu-Mo-Au-Re prospect, Lesvos Island, Aegean Sea, Greece}. Ore Geology Reviews [Internet]. 2019;112:103023. WebsiteAbstract
The Stypsi Cu-Mo-Au-Re prospect, Lesvos island, is a shallow porphyry-epithermal system hosted within a middle Miocene microgranite porphyry, which was emplaced along NNE-, NW- and NE-trending structures, within trachyandesites to trachydacites and felsic pyroclastics of the Stypsi caldera. The mineralization comprises three stages: Stage I is characterized by magnetite-actinolite ± quartz in transitional and banded quartz veins mostly developed within the microgranite and the surrounding lavas related to calcic-potassic and propylitic alteration of the host rocks. In the banded veins, quartz is botryoidal, suggesting crystallization from a gel. Abundant vapor-rich inclusions in the bands produce a dark grey to black color. Within the veins, magnetite, chalcopyrite, bornite and native gold were followed in time by pyrite, hematite, sphalerite and galena. Stage I mineralization was synchronous and also postdates quartz formation in the veins, since it crosscuts and/or fills vugs in the center of the veins. Sulfides are associated with various combinations of K-feldspar, actinolite, epidote, chlorite, and calcite gangue minerals. Pyrite-molybdenite-chalcopyrite (Stage II) and late intermediate sulfidation epithermal veins (Stage III) overprint earlier mineralization and alteration, and are associated with sericite-calcite and sericite-kaolinite alteration of the porphyry system and spatially associated lavas, respectively. Molybdenite is widespread in the first two mineralization stages, with Re content ranging from about 0.3 wt% in Stage I to 1.96 wt% in Stage II. A barren silicic and advanced argillic (alunite-kaolinite) lithocap, is exposed on top of the porphyry-style mineralization and is crosscut by non-mineralized high-sulfidation epithermal chalcedony-barite veins. Bulk ore analyses of surface samples from the Stypsi prospect yielded values of up to 276 ppm Mo, 978 ppm Cu, up to 0.5 g/t Au, up to 3 g/t Ag, and up to 70 ppb Pd and Re. Fluid inclusion data indicate that the Stage I transitional and banded quartz veins were deposited at 420 °C–530 °C and at pressures up to 450 bars, from boiling hydrothermal fluids. The fluid in the veins consists of a brine (40–61 wt% NaCl equiv) that coexists with a lower salinity (6–14 wt% NaCl equiv) liquid phase and a low-density vapor-rich fluid. Fluid inclusions in quartz of the Stage III intermediate-sulfidation epithermal veins are characterized by relatively low homogenization temperatures (231 °C–288 °C) and salinities (up to 1.9 wt% NaCl equiv), which was the result of subsequent dilution of the moderately saline fluids by circulating meteoric water. This study verifies earlier works suggesting that Au-enriched felsic magmas are able to crystallize Re-rich molybdenite, and that Re may also be redistributed and enriched in later stages during the deposition of porphyry-style mineralization. The Stypsi prospect ressembles in many respects (e.g., Au grades, Cu/Mo ratios, the Re content of molybdenite and the presence of ore-grade calcic-potassic alteration), other porphyry Cu-Mo-Re-Au systems hosted in calc-alkaline rocks in northeastern Aegean, Greece.
S} {L, K} {S, V} {A. {A guide to lithostratigraphy and structural inventory of the Cycladic Blueschist Unit rocks of NE Attica (Greece)}. Journal of the VIRTUAL EXPLORER [Internet]. 2019;50. Website
2018
Schneider DA, Grasemann B, Lion A, Soukis K, Draganits E. {Geodynamic significance of the Santorini Detachment System (Cyclades, Greece)}. Terra Nova. 2018;30:414–422.
Voudouris P, Mavrogonatos C, Rieck B, Kolitsch U, Spry PG, Scheffer C, Tarantola A, Vanderhaeghe O, Galanos E, Melfos V, et al. {The gersdorffite-bismuthinite-native gold association and the Skarn-Porphyry mineralization in the Kamariza Mining district, Lavrion, Greece}. Minerals. 2018;8.Abstract
Vein-type Pb-Ni-Bi-Au-Ag mineralization at the Clemence deposit in the Kamariza and “km3” in the Lavrion area, was synchronous with the intrusion of a Miocene granodiorite body and related felsic and mafic dikes and sills within marbles and schists in the footwall of (and within) the Western Cycladic detachment system. In the Serpieri deposit (Kamariza area), a porphyry-style pyrrhotite-arsenopyrite mineralized microgranitic dike is genetically related to a garnet-wollastonite bearing skarn characterized by a similar base metal and Ni (up to 219 ppm) enrichment. The Ni–Bi–Au association in the Clemence deposit consists of initial deposition of pyrite and arsenopyrite followed by an intergrowth of native gold-bismuthinite and oscillatory zoned gersdorffite. The zoning is related to variable As, Ni, and Fe contents, indicating fluctuations of arsenic and sulfur fugacity in the hydrothermal fluid. A late evolution towards higher sulfur fugacity in the mineralization is evident by the deposition of chalcopyrite, tennantite, enargite, and galena rimming gersdorffite. At the “km3” locality, Ni sulfides and sulfarsenides, vaesite, millerite, ullmannite, and polydymite, are enclosed in gersdorffite and/or galena. The gersdorffite is homogenous and contains less Fe (up to 2 wt.%) than that from the Clemence deposit (up to 9 wt.%). Bulk ore analyses of the Clemence ore reveal Au and Ag grades both exceeding 100 g/t, Pb and Zn > 1 wt.%, Ni up to 9700 ppm, Co up to 118 ppm, Sn > 100 ppm, and Bi > 2000 ppm. The “km3” mineralization is enriched in Mo (up to 36 ppm), Ni (>1 wt.%), and Co (up to 1290 ppm). Our data further support a magmatic contribution to the ore-forming fluids, although remobilization and leaching of metals from previous mineralization and/or host rocks, through the late involvement of non-magmatic fluid in the ore system, cannot be excluded.
Skilodimou HD, Bathrellos GD, Koskeridou E, Soukis K, Rozos D. {Physical and anthropogenic factors related to landslide activity in the northern Peloponnese, Greece}. Land. 2018;7.Abstract
The geological, geomorphic conditions of a mountainous environment along with precipitation and human activities influence landslide occurrences. In many cases, their relation to landslide events is not well defined. The scope of the present study is to identify the influence of physical and anthropogenic factors in landslide activity. The study area is a mountainous part of the northern Peloponnesus in southern Greece. The existing landslides, lithology, slope angle, rainfall, two types of road network (highway-provincial roads and rural roads) along with land use of the study area are taken into consideration. Each physical and anthropogenic factor is further divided into sub-categories. Statistical analysis of landslide frequency and density, as well as frequency and density ratios, are applied and combined with a geographic information system (GIS) to evaluate the collected data and determine the relationship between physical and anthropogenic factors and landslide activity. The results prove that Plio-Pleistocene fine-grained sediments and flysch, relatively steep slopes (15°-30°) and a rise in the amount of rainfall increase landslide frequency and density. Additionally, Plio-Pleistocene fine-grained sediments and flysch, as well as schist chert formations, moderate (5°-15°) and relatively steep slopes (15°-30°), along with the amount of rainfall of > 700 mm are strongly associated with landslide occurrences. The frequency and magnitude of landslides increase in close proximity to roads. Their maximum values are observed within the 50 m buffer zone. This corresponds to a 100 m wide zone along with any type of road corridors, increasing landslide occurrences. In addition, a buffer zone of 75 m or 150 m wide zone along highway and provincial roads, as well as a buffer zone of 100 m or 200 m wide zones along rural roads, are strongly correlated with landslide events. The extensive cultivated land of the study area is strongly related to landslide activity. By contrast, urban areas are poorly related to landslides, because most of them are located in the northern coastal part of the study area where landslides are limited. The results provide information on physical and anthropogenic factors characterizing landslide events in the study area. The applied methodology rapidly estimates areas prone to landslides and it may be utilized for landslide hazard assessment mapping as well as for new and existing land use planning projects.
Bathrellos GD, Skilodimou HD, Soukis K, Koskeridou E. {Temporal and spatial analysis of flood occurrences in the drainage basin of Pinios River (Thessaly, Central Greece)}. Land. 2018;7:1–18.Abstract
Historic data and old topographic maps include information on historical floods and paleo-floods. This paper aims at identifying the flood hazard by using historic data in the drainage basin of Pinios (Peneus) River, in Thessaly, central Greece. For this purpose, a catalogue of historical flood events that occurred between 1979 and 2010 and old topographic maps of 1881 were used. Moreover, geomorphic parameters such as elevation, slope, aspect and slope curvature were taken into account. The data were combined with the Geographical Information System to analyze the temporal and spatial distribution of flood events. The results show that a total number of 146 flood events were recorded in the study area. The number of flood events reaches its maximum value in the year 1994, while October contains the most flood events. The flood occurrences increased during the period 1990-2010. The flooded area reaches its maximum value in the year 1987, and November is the month with the most records. The type of damages with the most records is for rural land use. Regarding the class of damages, no human casualties were recorded during the studied period. The annual and monthly distribution of the very high category reaches the maximum values, respectively, in the year 2005 and in June. The analysis of the spatial distribution of the floods proves that most of the occurrences are recorded in the southern part of the study area. There is a certain amount of clustering of flood events in the areas of former marshes and lakes along with the lowest and flattest parts of the study area. These areas are located in the central, southern, south-eastern and coastal part of the study area and create favorable conditions for flooding. The proposed method estimates the localization of sites prone to flood, and it may be used for flood hazard assessment mapping and for flood risk management.
2017
Seman S, Stockli DF, Soukis K. {The provenance and internal structure of the Cycladic Blueschist Unit revealed by detrital zircon geochronology, Western Cyclades, Greece}. Tectonics. 2017;36:1407–1429.Abstract
This study presents detrital zircon U-Pb analyses of 23 samples of the Cycladic Blueschist Unit (CBU) from Kea, Kythnos, and Serifos islands, as well as the Lavrion Peninsula of SE Attica. The maximum depositional ages (MDA) and age distributions of detrital zircon U-Pb dates are used to correlate metasediments between the islands considered herein and infer their provenance. Two distinct detrital zircon U-Pb age distributions are found in CBU metasediments: “Proterozoic,” comprised of >40% Neoproterozoic zircons with Triassic-Early Jurassic maximum depositional ages and “Paleozoic,” containing >30% Paleozoic zircons and yielding Late Jurassic-Cretaceous MDAs. Proterozoic affinity metasediments are rift margin deposits derived from the northern Gondwanan margin. Paleozoic metasediments are flysch sediments most probably sourced from the Internal Hellenides. This metamorphosed flysch forms a distinct marker horizon found in a similar structural position in Lavrion, Kythnos, and Serifos. Based on lithologic correlation, sediment provenance, and MDA estimates, the CBU of Kythnos is correlative to the Lavrion Schists of Attica. On the islands of Serifos and Kythnos and within the Lavrion Schists only young-on-old relationships exist between rocks based on MDA estimates.
Malandri C, Soukis K, Maffione M, Özkaptan M, Vassilakis E, Lozios S, van Hinsbergen DJJ. {Vertical-axis rotations accommodated along the Mid-Cycladic lineament on Paros Island in the extensional heart of the Aegean orocline (Greece)}. Lithosphere. 2017;9:78–99.Abstract
The Aegean-west Anatolian orocline formed due to Neogene opposite rotations of its western and eastern limbs during opening of the Aegean back-arc basin. Stretching lineations in exhumed metamorphic complexes in this basin mimic the regional vertical-axis rotation patterns and suggest that the oppositely rotating domains are sharply bounded along a Mid-Cycladic lineament, the tectonic nature of which is enigmatic. Some have proposed this lineament to be an extensional fault accommodating orogen-parallel extension, while others have considered it to be a transform fault. The island of Paros hosts the only exposure of the E- to NE-trending lineations characterizing the NW Cyclades and the N-trending lineations of the SE Cyclades. Here, we show new paleomagnetic results from isotropic, ca. 16 Ma granitoids that intruded both domains and demonstrate that the trend difference resulted from post-16 Ma $\sim$90° clockwise and 10° counterclockwise rotation of the NW and SE blocks, respectively. We interpret the semiductile to brittle, low-angle, SE-dipping Elitas shear zone that accommodated this rotation difference to reflect the Mid-Cycladic lineament. We conclude a two-stage exhumation history for Paros that is consistent with regional Aegean reconstructions. Between ca. 23 and 16 Ma, the metamorphic rocks of Paros were exhumed from amphibolite-facies to greenschist-facies conditions along a top-to-the-N detachment. The Elitas shear zone then started to exhume the northwestern, clockwiserotating domain from below the southeastern, counterclockwise rotating domain since 16 Ma. From this, we infer that the Mid-Cycladic lineament is an extensional shear zone, consistent with geometric predictions that Aegean oroclinal bending was accommodated by orogennormal and orogen-parallel extension.
2013
Soukis K, Stockli DF. {Structural and thermochronometric evidence for multi-stage exhumation of southern Syros, Cycladic islands, Greece}. Tectonophysics. 2013;595-596:148–164.Abstract
The Attic-Cycladic metamorphic-core complex province was formed during Miocene to recent regional back-arc extension in response to slab-rollback and gravitational collapse. The development of several metamorphic-core complexes along extensional low-angle detachment faults resulted in the exhumation of lower-plate rocks and the tectonic denudation of the upper plate that is only sparsely exposed throughout the Aegean Sea. One of these areas is in southern Syros where the tectonic configuration includes three tectonometamorphic units separated by two detachment faults. The lower-plate Cycladic Blueschist Unit underwent late Cretaceous-Eocene HP/LT metamorphism and subsequent greenschist-facies overprint during Oligo-Miocene times. The tectonically overlying greenschist-facies Upper Unit is a tectonic sliver between the Cycladic Blueschist Unit and the structurally highest unit, the Vari Unit, which is dominated by quartzofeldspathic gneisses that are affected by four distinct deformation events (D1-D4). The late Cretaceous upper greenschist-facies D1 episode is related to late Cretaceous-Eocene subduction while D2-D4 fabrics were formed in response to late Cenozoic extension and exhumation. Field relations and structural data suggest that the late Miocene brittle Vari Detachment juxtaposes both Vari Unit and Upper Unit against the Cycladic Blueschist Unit. Mylonitic rocks observed at the base of Vari Unit along the contact with Upper Unit and zircon and apatite (U-Th)/He thermochronometry corroborate this scenario with lower-plate Cycladic Blueschist Unit rocks recording rapid cooling at $\sim$. 8-10. Ma. In contrast, cooling ages from the Vari Unit cluster at $\sim$. 13-15. Ma, indicating rapid exhumation coeval with activity along the Livadi Detachment on the neighboring island of Tinos. These data document new evidence that Vari Unit was affected by a major middle Miocene phase of exhumation and an older semi-brittle episode, both unrelated to late Miocene detachment faulting, revealing that crustal extension occurred in multiple stages and in a protracted episodic fashion similar to other parts of the Cyclades. © 2012 Elsevier B.V.
Bargnesi EA, Stockli DF, Mancktelow N, Soukis K. {Miocene core complex development and coeval supradetachment basin evolution of Paros, Greece, insights from (U-Th)/He thermochronometry}. Tectonophysics. 2013;595-596:165–182.Abstract
The Aegean region of Greece hosts a series of crustal-scale extensional detachment systems that have accommodated the southward retreating Hellenic subduction zone. Extension has overprinted and dissected the Alpine nappe pile and locally exhumed Cordilleran-type metamorphic core complexes. On the island of Paros, a low-angle extensional detachment fault separates metamorphic footwall rocks from an unmetamorphosed sedimentary succession of the hanging wall. Basement orthogneisses were extensionally sheared in the footwall of the detachment until after 16. Ma (zircon U-Pb age of a slightly deformed granite), but pervasive ductile deformation had ceased by 7. Ma (zircon U-Pb age of an undeformed rhyolite dike that intrudes gneisses). Apatite and zircon (U-Th)/He ages from the gneisses confirm a period of cooling at rates >. 100. °C/Ma from 16 to 7. Ma. In the upper-plate, the basal sedimentary unit yields reset detrital apatite (U-Th)/He (DAHe) ages from 17 to 7. Ma and detrital zircon (U-Th)/He (DZHe) ages ranging from 270 to 18. Ma. DAHe ages from the stratigraphically higher fanglomerate units are reset to 10-7. Ma. The DZHe data have a primary thermal signature of 12-7. Ma, but preserve ages up to 113. Ma. The uppermost conglomerates exhibit completely reset DAHe ages of 15-9. Ma and reset DZHe ages from 10 to 8. Ma, with DZHe ages up to 104. Ma. Reset DAHe ages indicate late exposure of the footwall and constrain the depositional age of most sedimentary rocks on Paros to be from 14 to 7. Ma. Unreset DZHe ages preserve thermal signatures from the major Mesozoic-Tertiary tectonic events in the Aegean Region: [1] Cretaceous Pelagonian-type metamorphism; [2] Eocene peak HP metamorphism; and [3] Miocene Barrovian overprinting. Preservation of these signatures indicates long-term upper-plate recycling prior to syn-extensional deposition. The Paros supradetachment basin represents a classic inverted unroofing sequence deposited during progressive core complex exhumation in the Middle to Late Miocene. © 2012 Elsevier B.V.
2006
Poscolieri M, Lagios E, Gregori GP, Paparo G, Sakkas VA, Parcharidis I, Marson I, Soukis K, Vassilakis E, Angelucci F, et al. {Crustal stress and seismic activity in the Ionian archipelago as inferred by satellite- and ground-based observations, Kefallinìa, Greece}. Geological Society Special Publication. 2006;261:63–78.Abstract
Different observational techniques are compared in order to investigate possible correlations in seismic activity. The study site is the island of Kefallin_gnia (Greece), where measurements available included (1) DInSAR, DGPS, and DEM data, (2) soil exhalation measured by monitoring Radon (Rn) well content, and (3) acoustic emissions (AE) at high and low frequency (point-like records with high temporal resolution). AE records provide: (1) relative time variation of the applied stress intensity and (2) the state of fatigue of stressed rock volumes, the AE source. Our results indicate that the large spatial scale (poor time resolution) may be considered quite satisfactory, whereas fractal analysis of the AE time series displayed some discrepancies when compared to analogous investigations in the Italian Peninsula. Therefore, some refinement is needed in order to reach more precise interpretations of the relevant information available with this kind of data. However, both sets of observations appear in agreement with each other, although more exhaustive investigations would require a suitable array of point-like AE and Rn (or other) measuring sites, as well as longer data series. The latter are particulary helpful for detailed interpretations of the different occurrences within tectonically complex settings where crustal stress crises are marked by various types of geological phenomena. © The Geological Society of London 2006.