The Volcanic Relief within the Kos-Nisyros-Tilos Tectonic Graben at the Eastern Edge of the Aegean Volcanic Arc, Greece and Geohazard Implications

: The active Kos-Nisyros-Tilos volcanic ﬁeld is located in the eastern sector of the Aegean Volcanic Arc resulting from the subduction of the African plate beneath the Aegean plate. The volcanic activity is developed since Middle Pleistocene and it occurs within a tectonic graben with several volcanic outcrops both onshore and offshore. Data obtained from previous offshore geophysical surveys and ROV exploration, combined with geospatial techniques have been used to construct synthetic maps of the broader submarine area. The volcanic relief is analyzed from the base of the volcanic structures offshore to their summits onshore reaching 1373 m of height and their volumes have been computed with 24.26 km 3 for Nisyros Island and a total volume of 54.42 km 3 for the entire volcanic area. The volcanic structures are distinguished in: (1) volcanic cones at the islands of Nisyros (older strato-volcano), Pergousa, Yali and Strongyli, (2) volcanic domes at the islands of Pachia, East Kondeliousa and Nisyros (younger Prophitis Ilias domes), (3) submarine volcanic calderas (Avyssos and Kefalos). Submarine volcanic debris avalanches have been also described south of Nisyros and undulating features at the eastern Kefalos bay. Submarine canyons and channels are developed along the Kos southern margin contrary to the Tilos margin. Ground truth campaigns with submarine vessels and ROVs have veriﬁed the previous analysis in several submarine volcanic sites. The geohazards of the area comprise: (1) seismic hazard, both due to the activation of major marginal faults and minor intra-volcanic faults, (2) volcanic hazard, related to the recent volcanic structures and long term iconic eruptions related to the deep submarine calderas, (3) tsunami hazard, related to the seismic hazard as well as to the numerous unstable submarine slopes with potential of gravity sliding.


Introduction
The Kos-Nisyros volcanic field is located at the eastern edge of the Aegean Volcanic Arc, developed behind the Hellenic arc and trench system in Greece [1][2][3][4][5][6] (Figure 1). The onshore volcanic outcrops on Kos, Nisyros and Yali islands have been studied by several researchers [7][8][9][10][11][12][13][14][15][16][17][18]. More recently, offshore evidence for volcanic activity in the area of Kos-Nisyros has been reported on the basis of volcanic tephra layers occurring within the uppermost sediments obtained from coring [19][20][21][22]. Since the late 1990's offshore systematic research has been carried out in the subsided graben structure area between Kos and Tilos islands resulting in the discovery of several submarine volcanic outcrops [23][24][25][26][27]. This subsided part of the graben structure is made up of marine basins with depths around 500 m-700 m and volcanic islands in between mainly along the central axis of the Modern oceanographic survey techniques involving high-resolution swath bathymetry, side-scan sonar and ROV (remotely operated vehicle) exploration greatly facilitated the recognition and research of submarine morphological features in the volcanic environments of the area using digital data similarly to other insular volcanoes and seamounts [31][32][33][34][35].
In this paper, we describe the volcanic relief within the Kos-Tilos tectonic graben, combining onshore and offshore data, based on detailed bathymetric maps and morphological slope distribution maps. The various volcanic features such as volcanic domes, craters and calderas, are mapped and distinguished within the sedimentary basins where volcanosedimentary formations of more than 700 m thickness have been detected through seismic profiling [24,27,36]. The overall area of the volcanic submarine outcrops has been measured and the overall volume of extruded volcanic rocks has been computed. Additionally, the submarine channel network system has been analyzed using geospatial algorithms. Finally, some images of submarine outcrops which were taken from the submersible «THETIS» or from ROV systems have been used for ground truth. The morphological analysis of the Nisyros volcanic field can be used for understanding the magnitude and the rates of the active geodynamic processes at the eastern edge of the Aegean Volcanic Arc and to make a comparison with other geodynamic regimes.

Geological-Geodynamic Setting
The Aegean volcanic arc is the result of the northeastward subduction of the African plate beneath the active European margin throughout Late Mesozoic-Cenozoic [3]. Since the Late Miocene-Pliocene the actual subduction of the East Mediterranean oceanic basin started with an increasing subduction rate contrary to the previous slow rate of the continental subduction of the External Carbonate Platform of the Hellenides [37]. Thus, during the Plio-Quaternary the Aegean volcanic arc was formed approximately at its present location, comprising several volcanic centers from Methana to the NW, through Milos and Santorini to Nisyros to the east ( Figure 1). The first discovery of the Paphsanias submarine volcano in the Epidaurus Basin [38,39] was followed by the discovery of several other submarine volcanoes especially around the Kos-Nisyros islands [24].
The general tectonic structure of the Kos-Nisyros volcanic field is developed within a neotectonic graben of ENE-WSW direction, bordered by the adjacent tectonic horsts of Kos to the NW and Tilos to the SE [25,40]. Alpine basement rocks crop out on both sides of the Nisyros graben forming the mountainous relief with several hundred m elevation (900 m in Dikeos Mt. of Kos and 600 m in Tilos). The Kos alpine basement comprises several tectonic units with a Palaeozoic metamorphosed unit at the base occurring at the uplifted horst structure of Dikeos Mt., overlain by the non-metamorphosed nappes of the neritic Zia (Tripolis type), the pelagic Prophitis Ilias (Pindos Type) and the Eastern Kos mélange [8,41]. The subsided area of the tectonic graben comprises several volcanic islands and a number of submarine volcanic features within the marine basins. The volcanic activity started in Late Miocene in Kos Island with the intrusion of the "Kos Monzonite" at the southwestern slopes of Dikeos Mt. [8,42] and some small outcrops of Upper Miocene-Pliocene volcanics. A distinct volcanic activity occurred during Middle Pleistocene at about 500 Ka at Kefalos peninsula in Western Kos, related to the Kefalos caldera [43]. However, large volumes of volcanic products were extruded in the form of the "Kos Ignimbrite" which represents the largest eruption in the Eastern Mediterranean during the Late Middle Pleistocene at approximately 0.16 Ma ago [6,14,44]. The explosion center of this eruption is thought to be located at the "Avyssos Caldera" northeast of the Strongyli volcanic Islet [26,36] (see Figures 3 and 14). Contrary to the Kos long volcanic history there is no volcanic activity in Tilos, which lies in a more external position outside the volcanic arc, except for some pumice layers deposited during the Kos and/or Nisyros strong explosive events [45].
The volcanic island of Nisyros features large volumes of volcanic products composed mostly of Late Pleistocene volcanic rocks represented by alternating lava flows, pyroclastic layers and viscous lava domes [6,8,9,11,17,46]. The major distinction in the volcanic history of Nisyros is the first period of the strato-volcano formation which ended by a major eruption (Nikia rhyolites) and caldera formation, followed by the second period of formation of volcanic domes, disrupting the former caldera rim (presently observed at about 300 m of elevation), and forming the highest actual mountain of Prophitis Ilias (698 m). In the submarine area around Nisyros recent volcanic formations have created the volcanic centers of Pergousa, Yali, Strongyli, Pachia and East Kondeliousa [18,24]. However, the largest part of the submarine area is mostly composed of more than 700 m thick marine volcano-sedimentary deposits which overlie the alpine basement of the Nisyros graben [24,36].
The Kos-Tilos graben is subdivided in two sub-grabens/sub-basins separated by the Kondeliousa intermediate horst structure, where the alpine basement crops out. Immediately east of the Kondeliousa Island the submarine East Kondeliousa volcanic domes have been detected and their prolongation to the ENE reaches Pachia Island and Prophitis Ilias domes of Nisyros Island [24,36]. Three basins are developed towards Kos Island to the north and two basins to the south, whereas west of the Kondeliousa Islet a shallow water platform is developed on top of the alpine basement [28].
The vertical displacements between the neotectonic blocks, as estimated from the throws of the marginal faults, are about 1-2 km [25,36]. Some of the major faults display high escarpments reaching tens to hundreds of meters high and relative uplift of about 1 km, such as the marginal faults observed on both sides of the Kondeliousa tectonic horst [25,36]. The volcanic extrusions have penetrated the central part of the tectonic graben forming the present-day Quaternary volcanic structures of Nisyros, Pachia, Pergousa, Yali and Strongyli, which are aligned in the NE-SW direction [24,27]. Small intra-volcanic marine basins with less than 350 m depth have been formed between Nisyros and the other volcanic centers. Nisyros Volcano is dissected by major fault zones forming a radial pattern around the Nisyros caldera [11,13,40,[46][47][48].

Collection of Data and Methodology
The large data set was obtained from three multibeam bathymetric surveys carried out on R/V «AEGAEO» during 2000 in the area between Kos, Nisyros and Tilos islands. During the first mission, the area of Nisyros Island and the surrounding small islets has been mapped using the SEABEAM 1180 (180 kHz) system, which is suitable for seabed mapping in shallow-middle depths (<500 m). The SEABEAM 1180 is mobile, shallow water, compact system integrating transmitter, receiver, interfaces and power stage within a single unit. The other two missions completed the mapping of the whole area of Kos-Nisyros-Tilos using the SEABEAM 2120 (20 kHz) system, which is suitable for seabed mapping in depths >500 m. The SEABEAM 2120 is a relatively new swath system that has been especially designed to accomplish survey requirements exceeding 6000 m of water-depth, providing a satisfactory resolution (up to 1 • × 1 • ), without mounting a very large array. By operating the systems for a total period of 12 days with an average speed of 10 knots, 3500 km 2 were covered from very shallow depths to depths of 2200 m.
Several dives with the submersible «THETIS» supported by the R/V «AEGAEO» were carried out in November 2000. THETIS carries one pilot and one scientist/observer and its diving period can be 3-6 h, with an operational depth of 610 m.
ROV exploration of the submarine volcanism around Nisyros volcano took place in October 2010 using the Exploration Vessel (E/V) Nautilus. The dives of ROVs "Hercules" and "Argus" aimed to explore the submarine volcanic centers and active structures previously detected around Nisyros with the submersible THETIS of HCMR, in 2000. The E/V Nautilus of O.E.T. (Ocean Exploration Trust) is equipped with the ROVs Hercules and Argus, which are state-of-the-art, deep-sea robotic vehicle systems capable of exploring depths up to 4000 m [49]. Both vehicles have high-definition video cameras which provide a clear and high resolution imaging of the seafloor.
The previously described surveys within the Kos-Nisyros-Tilos submarine environment have resulted in very large data sets based on high resolution offshore data. A part of the processing of the bathymetric data so as to produce the Digital Terrain Model (DTM) of the area, has been achieved using NEANIAS's thematic service U1-BAT. The U1 [UW_BAT] "Bathymetry mapping from acoustic data service implementation" of NEANIAS provides an advanced user-friendly, cloud-based version of the popular MB-System software, for processing bathymetric data through Jupyter Notebooks [50]. After the processing of the bathymetric data, a specific workflow was applied in order to identify the morphological features and the general structure of the area. In order to produce valuable and comparable outcomes, preparation and filtering of the original geospatial data were applied, before the thorough spatial analysis. According to the following schematic diagram (Figure 2) a specific workflow was implemented hierarchically using spatial analysis and particularly hydrology analysis toolset within the environment of ArcGIS 10.4 version.
Most of the morphological features of the studied area were identified and digitized semi-automatically or manually based on the post processing data. The first step was the investigation of the multibeam bathymetric data (with resolution of 20 m × 20 m) in order to remove possible sinks and peaks from the grid using the fill algorithm. According to ESRI sinks are often errors due to the resolution of the data or rounding of elevations to the nearest integer value. The sinks should be filled to ensure proper delineation of basins and streams. If the sinks are not filled, a derived drainage network may be discontinuous. The fill tool uses the equivalents of several tools, such as Focal Flow, Flow Direction, Sink, Watershed, and Zonal Fill, to locate and fill the sinks [51]. Thus the tool iterates until all sinks within the specified z limit are filled. Subsequently, after the implementation of the sink tool more than 29,000 sink areas had been identified and removed from the initial multibeam bathymetric dataset. After the correction of the bathymetry, the extraction of the submarine drainage network was performed regarding the Strahler technique with the use of the flow accumulation algorithm [52]. Flow accumulation is a tool to generate a stream network based upon the amount of pixels flowing downslope. In our work the selected threshold pixel value for the derived flow accumulation dataset was set to 10,000 according to [53]. This preferred value was applied in order to filter data due to the inherent noise of the multibeam data. In some areas the derived drainage network was modified manually in order to correct topological errors between the interconnected branches. An orientation analysis of the tectonic lineaments was performed and the results were presented in a rose diagram (Figure 3) produced by a script written in the Python programming language designed for the geometric and topological analysis of fracture networks [54].

Seabed Morphology
The compiled bathymetric map was combined with onshore topographic maps of the islands and the result was a synthetic topographic map of the entire area both onshore and offshore ( Figure 3). The processed bathymetric data were classified and merged with onshore data ranging from +698 m (Nisyros highest elevation point) to more than −1000 m (part of Karpathos Basin, further south of the Kondeliousa horst). The marine area under investigation has an estimated extent of 2564 km 2 based on geospatial measurements according to the Hellenic projected coordinate system EGSA'87.
The subsided area between the Kos-Tilos graben structure comprises five main basins, three basins north of the median Kondeliousa horst/Nisyros volcanic ridge zone and two south of this zone. The Eastern Kos Basin with maximum depth 667 m, the Western Kos Basin with maximum depth 514 m, the Kondeliousa Basin with maximum depth 594 m, the Northern Tilos Basin with more than 500 m depth and the Southern Nisyros Basin with more than 800 m depth, which continues to the south towards the major Karpathos Basin (more than 2500 m depth). Prominent morphological structures within the studied area are the Kondeliousa plateau with average depths around 400 m and the Western Tilos shelf with average depth 300 m. Small basins with depths around 300 m are also observed between Nisyros and the adjacent volcanic islets of Yali, Pergousa and Pachia. The major morphological distinction within the Kos-Tilos graben is outlined by the 400 m isobaths, which separate the three basins to the north from the two basins in the south.
The rose diagram ( Figure 3) presenting the prevailing orientation of the faults, of ENE-WSW direction, parallel to the general orientation of the Kos-Tilos tectonic graben. It is remarkable that there is only a small deviation of 30 • of the orientation with absence of perpendicular trends.

Morphological Slopes
Another important part of the geoprocessing procedure was the conversion of the elevation/depth data in a slope map classified into four categories starting from subhorizontal areas with slope values 0-2 • , low slope areas 2-10 • , medium slope areas 10-20 • and steep slope areas >20 • . After the implementation of the slope algorithm, the next process was to detect the basinal submarine areas and the steep slopes in order to semiautomatically identify and delineate the rest of the submarine geomorphological features around Kos-Nisyros-Tilos area ( Figure 4). This morphological slope map is a good starting point for providing an effective means to identify the presence of active fault zones, various volcanic features (volcanic craters, calderas and domes), major morphological features (basins, plateaus) and gravitational-erosive-depositional features (channels and canyons, debris avalanches, seafloor undulating features).
The sub-horizontal areas define: (1) the five basins previously distinguished on the basis of bathymetry; the Eastern Kos Basin, the Western Kos Basin, the Kondeliousa Basin, the Southern Nisyros Basin and the Northern Tilos Basin, (2) the shallow platform west of Kondeliousa Island and (3) the continental shelf area west of Tilos Island (Figure 4 and Figure 14).
The northern sub-graben between Kos and the median Kondeliousa horst/Nisyros volcanic ridge comprises three of the five basins ( Figure 5).   The steep morphological slopes vary between 25 • north of Eastern Kos Basin due to a major neotectonic fault running parallel to the coast for more than 15 km and 50 • north of Yali-Kos rise (area between Eastern and Western Kos basins). The steep slopes are observed either in the margins of the tectonic graben structure or around the volcanic islands ( Figure 4). The zones of steep slopes occurring along the margins are linear with a general NE-SW orientation. This is observed mainly in the southern margin of central and eastern Kos Island, in both margins of the Kondeliousa median horst/platform and in the western margin of Tilos, which is located at 10 km-15 km distance from the NW Tilos coastline. The exceptions are the northern Tilos margin oriented in the E-W direction and the western part of the Kos margin oriented in the N-S and WNW-ESE direction. On the contrary, the rest steep slope areas are observed around the volcanic islands in the middle of the graben structure with circular forms defining volcanic cones or more elongated forms defining volcanic domes. Two special cases of negative volcanic forms corresponding to volcanic calderas are observed in the Eastern Kos Basin and the Western Kos Basin respectively.  (Table 1). Two main volcanic features are present:

Volcanic Structures and Volcanic Relief
(1) The stratovolcano which culminated with the major eruption following the Nikia rhyolites and the formation of the Nisyros caldera [17] and (2) the Prophitis Ilias rhyodacitic domes (Figure 6), which intruded the western part of the caldera rim. The highest post-caldera volcanic domes reach 698 m of altitude and form the Prophitis Ilias peak delimited between two onshore perpendicular fault zones [11] (Figure 3). The lavas of Prophitis Ilias have flowed offshore in a southwesterly direction extending approximately down to 166 m water depth. Based on geospatial calculations, the base of Nisyros volcano has a mean diameter of 11 km and the total mass volume is estimated to be 24.26 km 3 (including the rhyodacites of Prophitis Ilias). According to the slope map ( Figure 4) the volcanic cone is characterized by more than 20 • slope. In addition, the slopes of the volcanic domes reach a maximum value of 56 • . The NE-SW profile of Nisyros Island ( Figure 6) shows the Prophitis Ilias continuous dome structure from outside the pre-existing Nisyros caldera up to its central area whereas the NW-SE profile shows its position within the pre-existing caldera rim of the strato-volcano. (ii) The Yali volcanic edifice (Figure 7a) is built from 501 m water depth and emerges at 175 m onshore, with an overall volcanic relief of 876 m (Table 1). It has a semi-cone structure forming the western half of the island with pumice layers under exploitation. On the contrary, the eastern half of the island is made of obsidian and rhyolitic lavas. The two islets of Yali are separated by a narrow land corridor where a N-S striking sinistral strike-slip fault has brought in contact two different segments of the overall Yali volcanic edifice [24]. The small Yali-Nisyros submarine basin (275 m of depth) is observed south of Yali with less than 2 • of slope. Its western margin occurs at 155 m of depth, controlled by a N-S fault with more than 120 m of throw (pink line in Figure 7a). This submarine fault continues to the south into the Mandraki fault activated during the 1996-1997 seismic activity [36,55].  (Table 1). It forms a volcanic cone with a small crater of a few hundred m diameter. Immediately northeast of Pergousa islet a submarine ridge structure with elongated shape and E-W axis direction has been characterized as a volcanic dome.     The above volcanic structures and their onshore/offshore volcanic relief are given in the synthetic Table 1.
Additionally, several submarine volcanic intrusions are observed in the form of elongated domes in the NE-SW direction northwest of Kondeliousa Island within the Kondeliousa Basin as previously described.
Besides the above volcanic structures producing positive relief there are also two submarine volcanic calderas creating negative relief, discovered from the integration of swath mapping and air-gun seismic profiling [24,36] (Figures 8 and 9): (i) the Avyssos caldera is located northeast of the base of Strongyli volcanic cone and it was probably formed after the massive eruption of the Kos ignimbrite 0.16 Ma ago [6,14,44]. This assumption has been based on the location of the Avyssos caldera at the central area of the mega-eruption based on the widespread ignimbrite deposits and on the fact that only a thin sedimentary cover was detected on the volcanic basement [24,36]. The caldera is elliptic with 3 km length in the NW-SE direction and 4 km in the NE-SW direction (Figure 8). Its base occurs at around 680 m depth and its caldera rim lies at about 630 m depth, forming a submarine circular cliff of 40 m-50 m topographic difference [25]. The Avyssos caldera rim is characterized by 5-6 • slope, contrasting the 0-2 • of the planar floor ( Figure 4). Some small hills are observed inside the planar base of the caldera with a relative relief of about 50 m-70 m above its base. This internal caldera area is characterized by a thin cover of fine-grained sediments overlying volcanic formations as verified during the dives of the Nautilus cruise NA011, in 2010 [26]. However, no evidence of hydrothermal activity was detected in the Avyssos caldera. The peak of the central volcanic dome lies at 610 m depth. These intra-caldera hills represent later intrusions of volcanic rocks [26] (Figure 8).
(ii) the Kefalos caldera, is observed at the south-western part of Kos Island, to the southeast of the Kefalos Peninsula and its formation occurred during the Middle Pleistocene (0.5 Ma) [43]. The offshore caldera rim is characterized by steep slopes of 20-30 • and the caldera base lies at depths from 220 m to 300 m ( Figure 9). Only the western part of the Kefalos caldera is preserved whereas its eastern continuation has collapsed and the result is an undulating morphology of arcuate terraces reaching 400 m to 450 m water depth, similar to La Fossa Caldera at Vulcano which is partially submerged [56,57]. The undulating feature is characterized by 10-20 • of slope and resembles a radial flow southwards to the basinal area of the Western Kos Basin (see also Figure 14).  A special morphology indicating volcanic debris avalanche deposits is observed on the submarine slopes southeast and southwest of the Nisyros volcanic structure. A number of 364 hills/hummocks were counted through geoprocessing procedure throughout the avalanche field ( Figure 10). They form clusters at varying depths between 196 m to 769 m and extend as far as 20 km away from the Nisyros coast. The shape of these hummocks varies from circular round to elongated and asymmetrical blocks of different orientations, having an average axial length of 350 m (the larger ones). A distinct debris avalanche deposit has been studied in detail on the southeast flank of Nisyros [16,[58][59][60]. It is likely that it originated onshore and covers an area of about 90 km 2 . The overall volume is estimated about 1 km 3 and the average thickness 40-50 m at water depths between 250 m and 410 m. The blocks display elongation parallel to the flow of the debris avalanches from the Nikia rhyolites cropping out at southeast Nisyros Island [16]. The volcanic debris avalanche comprises numerous hills rising up to 60 m above the sea bottom and longitudinal ridges, revealing a horseshoe-shaped structure. Figure 10. Hillshade map of the SE Nisyros submarine area, with 10 m isobaths, depicting a part of the volcanic debris avalanche field. The larger hills/hummocks are shown as closed isobaths in the map (yellow arrows). Slope profiles (in degrees) along the debris avalanches interpret their shape to horizontal distance.

Submarine Drainage System
The streams onshore discharge to the sea, where they continue their route on the seafloor, forming a well-organized submarine drainage network, which ends up into the three submarine basins ( Figure 11).
According to the Strahler' technique four orders of streams are developed in the submarine region. Streams of 1st, 2nd and 3rd class have been formed along the southern continental slopes of Kos Island. Along the slopes of the Kos margin a number of submarine canyons have been formed, probably by fast-flowing turbidity currents carrying terrigenous sediments. The canyons usually correspond to 1st or 2nd order streams/channels, characterized by steep walls and high incision characteristics along steep morphological slopes [61]. They are observed from their head at shallow depth down to the smoothing of the slopes at the margin of the sub-horizontal basinal areas. The canyons usually continue as higher order channels in the basinal area where they may gradually disappear by discharging their transported material in the sea bed. Thus, along the southern slopes of Kos Island, these submarine channels start from 150 m of depth to 550-640 m where the morphological slopes values become very low. The transport and deposition of sediments from the channels form fan-shaped deposits, especially at the Eastern Kos basin. Thus, the drainage network represents the continuation of the onshore one. The dominant drainage pattern at the Kondeliousa Basin is dendritic and follows the NE-SW ridges orientation of the submarine morphology. Streams up to the 4 th order extend along the northern flanks of the Kondeliousa Basin seafloor, forming channels with deposition of unconsolidated sediments in between the parallel submarine volcanic outcrops as detected in the seismic profiles [27,36] The main characteristic of the drainage system along the southern Kos margin is the absence of a major stream, oriented ENE-WSW, receiving the N-S to NNW-SSE channels. This is due to the fact that the transported material is deposited as long as it reaches the basinal areas of the Eastern and Western Kos basins. On the contrary, the minor N-S streams of the northern margin of the Kondeliousa Basin join the higher order streams of ENE-WSW orientation draining to the west, outside the studied area.

Submersible Diving and ROV Exploration
Several submarine volcanic structures have been observed at various depths down to 504 m through submersible diving and ROV exploration. Submersible dives with submersible THETIS of the HCMR (Hellenic Center for Marine Research) were operated during 2000 at the NW flank of the volcanic cone of Strongyli and at the subvertical slopes of the East Kondeliousa volcanic domes [24]. Some volcanic craters were discovered at 240 m of depth in Strongyli and at 430 m of depth in Kondeliousa ( Figure 12). Small volcanic craters were also discovered on the northwest slopes of Strongyli volcanic cone by ROV in 2010 [26,36]. Several dives with submersible THETIS made observations at the northern abrupt slopes of the volcanic cone of Pergousa and Pachia volcanic dome starting from the depth of 380 m and 250 m respectively. Additionally, submarine dives made observations at the western slope of the Nisyros-Yali basin along the N-S striking fault zone where a landslide occurred probably during the seismic activity of July 1996 [55].
The ROVs Hercules and Argus during the "Nautilus" cruise in October 2010, were deployed in five locations (the Avyssos crater, the slopes of Yali, Strongyli, and Kondeliousa Islands, and a suspected debris avalanche deposit southeast of Nisyros) [26]. ROV exploration of the eastern flank of Yali revealed wavy sedimentary structures with dark sediments on the ridges and light sediments in the troughs, as well as linear fractures at various depths. Small craters were discovered on the northwest slopes of Strongyli, aligned with ENE-WSW trending fractures, without any sign of hydrothermal activity. The outcrops of the volcanic rocks on the flanks of Yali and Strongyli volcanoes are covered by heavy biogenic encrustations. Bio-encrustations were observed on the steep cliffs of the volcanic domes of Kondeliousa with large fields of feather corals bearing elasmobranch egg cases on them, as well as patches of biogenic sediments with abundant sponges (Figure 13). The ROV dive on this region confirmed the tectonic boundary between the Mesozoic limestones of Kondeliousa Island to the southwest and the recent volcanic rocks towards the northeast.

Synthesis and Geohazard Implications
The synthesis of the previously presented data is given in the morphotectonic map of the submarine area of the Nisyros graben between Kos and Tilos islands ( Figure 14).
The Kos-Tilos tectonic graben shows ENE-WSW orientation in the eastern part which becomes NE-SW in the western part towards the Karpathos Basin with increasing depth. The recent volcanic activity spreads out along the central area of the graben around Nisyros Island. To the north on Kos Island the volcanic activity is present since late Miocene with more recent eruptions the Kefalos caldera 0.5 Ma ago and the deposition of the Kos ignimbrite 0.161 Ma ago [14,22]. On the contrary, to the south, on Tilos Island [46], there is no volcanic activity, except for the deposition of some ignimbrite pumice outcrops of the Kos ignimbrite mega-eruption [44]. The volcanic structures can be distinguished in: (1) a number of volcanic cones and calderas such as the Nisyros stratovolcano cone and caldera, the Yali, Strongyli and Pergousa cones and the Avyssos caldera. (2) several volcanic domes oriented NE-SW along a distinct zone from the East Kondeliousa, Pachia and Prophitis Ilias domes of rhyodacite composition. Another group of volcanic domes/dikes is observed at the basinal area of the Kondeliousa Basin at depths around 500 m, which are also oriented NE-SW. It is remarkable that submarine debris avalanches made of volcanic rocks are observed around the Nisyros major volcanic structure mainly towards the southwest to southeast slopes at distances up to 15 km-20 km [16,59,60]. The two margins of the graben are very different with numerous channels and canyons dissecting the abrupt slopes of the Kos southern margin ending down to 500-600 m depth in the Eastern and Western Kos basins contrary to the Tilos north-northwest margin where they are absent. Additionally, the Kos margin is very abrupt and is developed very close to the coast line contrary to the Tilos margin which extends much further to the west, forming a sub-horizontal plateau at about 300 m depth. Another special structure is the Kondeliousa horst which towards the west is developed upon the Alpine Mesozoic limestones and forms a sub-horizontal plateau at about 400 m depth. This plateau is interrupted by a fault immediately east of the Kondeliousa islet and at its continuation to the NE it is substituted by the East Kondeliousa, Pachia and Prophitis Ilias volcanic domes.
The volcanic relief emerges from the average depth of the basinal areas of 600 m and reaches 700 m of height on the Nisyros post-caldera domes. Thus, an overall volcanic edifice of about 1300 m is built within the Kos-Nisyros-Tilos graben structure since Middle Pleistocene.
The geological hazards in the Nisyros volcanic area within the Kos-Tilos graben comprise: (1) Seismic hazard resulting from the activation of the faults. This hazard is different in trems of the activation of marginal main faults with a few tens of km length, such as the magnitude 6.6 1933 event on Kos, which affected also Nisyros Island [62] than in the case of intra-volcanic edifice faults of a few km length, such as the magnitude 5.6 1996-1997 events on Mandraki/Nisyros [63], which activated the N-S fault of northern Nisyros and its offshore prolongation towards Yali Island [55]. More recently, in 2017 a 6.3 event occurred at another main fault east of Kos Island between the Greek and Turkish waters with significant losses [64]. In several cases the earthquake epicenters occur offshore along active marginal faults trending ENE-WSW to E-W that border the present marine basins from the islands [27]. This tectonic trend is regional characterizing the eastern part of the Aegean volcanic arc in the back arc area of the actual Hellenic arc and trench system (see also Figure 1) resulting from a regional N-S extension [65]. (2) Volcanic hazard occurs in the active craters of Nisyros, Yali, Strongyli, Kondeliousa and to a lesser degree in Pergousa and Pachia. The last minor eruptions have occurred on Nisyros in 1887 [17]. The recurrence period of such eruptions may be estimated from the voluminous volcanic products during the Late Pleistocene and Holocene in the area during a few thousands of years. Distinction should be made between the Plinian type explosions with deposition of pumice layers both onshore and offshore from the lava extrusions in the form of rhyodacitic domes. However, iconic megaeruption similar to the Kos ignimbrite eruption 161 Ka ago may be expected from the volcanic activity in the submarine Avyssos caldera. This event has more than 100 thousand years cycle and represents a major volcanic threat for the broader area of the eastern Dodecanese islands and the coastal zone of Minor Asia. Submarine volcanic eruptions may also occur along the E-W dikes/domes of the Kondeliousa Basin but their impact is not expected to be important in the onshore areas. (3) Tsunami hazard is expected due to seismic activity of the marginal fault zones that may generate major earthquakes of magnitude 6-7 similar to the 1956 Amorgos event [66][67][68]. Additional tsunami hazard is expected due to submarine sliding, as observed at the broader zones of volcanic debris avalanches especially at the SE and SW slopes of the Nisyros volcanic center and at the undulating slopes of the disrupted Eastern Kefalos caldera. However, tsunami deposits have not been reported up to present in the surrounding coastal areas, although there are several sites with appropriate lagoon type environment, especially in Kos Island.
Since Middle Pleistocene the volcanic activity is taking place within the Kos-Nisyros-Tilos graben structure either on the islands or in the submarine surrounding area. Thus, the base of the volcanic structures occurs at an average depth of 600 m and culminates at several hundred m higher either below sea level or at several hundred m altitude. It is characteristic that on Nisyros Island the first lavas are pillow lavas and the transition to aerial explosions occurs higher in the volcano-stratigraphy of the first stratovolcano period of Nisyros [11,36]. However, this area paleogeographically belonged to the Neogene Aegeis land bridge, developed between continental Greece and Minor Asia [69]. Thus, in Kos there are continental deposits with lacustrine fossiliferous formations of Upper Miocene-Pliocene age [8,41]. The same facies occurs also in the other east Aegean islands both towards the north (e.g., Samos and Chios) and the south (e.g., Rhodes). The first post-Alpine marine outcrops overlying Aegeis land mass occur in southeast Ikaria Island deposited during Early Pliocene [70] and in Kos and Rhodes islands deposited during middle-late Pleistocene [71,72]. Thus, the subsidence observed within the Kos-Nisyros-Tilos tectonic graben is post early Pleistocene and it is accompanied also by submarine volcanism. The volcanic structures comprise submarine basal volcanic formations and onshore upper parts with subaerial explosions, observed in the volcanic islands. The onshore volcanic outcrops of middle-late Pleistocene age observed on Kos, Tilos and surrounding areas are pyroclastic-pumice-tephra products of volcanic activity occurring within the Nisyros volcanic group [15,17,44]. The mid Pleistocene eastern Kefalos volcanic event is located at the present coastal area of western Kos [43]. The Kefalos caldera was combined with its upper part onshore and the middle and lower part offshore. However, only its western part is preserved whereas its eastern part has been disrupted and forms a submarine undulating slope zone towards the Western Kos basinal area. The mixed nature of the Kefalos caldera resembles the multi-stage formation of La Fossa Caldera at Vulcano Island, Italy [54,56,57].
Strongyli volcanic cone is differentiated from all the other volcanic islands because of its geometry and the absence of any pyroclastic/pumice deposit on it. Additionally, its location at the southwestern margin of the Avyssos caldera implies a probable post caldera development, as the contact between the Strongyli base and the SW margin of the Avyssos caldera indicates. Nevertheless, there are human constructions at the small, planar area of its crater with cultivation terraces. The 1.8 km overall subsidence of the Alpine basement in Nisyros Island [29,30] is compensated by the building of the volcanic relief in at least two stages (stratovolcano and post caldera domes). This is observed also in the ENE-WSW oriented zone of volcanic domes observed from the East Kondeliousa submarine domes to the Pachia and Prophitis Ilias domes of Nisyros Island. This volcanic zone is located at the middle of the Kos-Nisyros-Tilos graben structure and follows its geometry and geodynamic characteristics with extension in the NNW-SSE direction.
The main marginal faults with important seismic events are also oriented in ENE-WSW direction and this is also observed in the overall morphological structure of the area. Additionally, the same direction is noticed along the domes of the Kondeliousa Basin. NNW-SSE extension is detected in the GPS rates of the area [65] as well as in the fault plane solutions of large earthquakes [73].
The extended submarine areas of debris avalanches to the SE and SW of Nisyros Island resemble the hummocky terrains formed by the debris avalanches with alignment of the hummocks in the collapse affecting the E flank of Stromboli [74,75]. They represent a potential tsunami hazard triggered by the slope instabilities of the volcanic islands as this is known from several cases worldwide [76,77]. Slope instabilities are observed also along the southern margin of Kos Island and especially at the eastern Kefalos collapsing caldera with undulating zones, where various types of landslides may occur during seismic activity [78,79].
In conclusion, the volcanic hazard of the area is of major concern, although since 1887 there is no activity, contrary to the Santorini intense activity throughout the 19th and first half of 20th century [80][81][82]. An unrest period was observed during 1996-1997 with earthquakes and minor changes of the thermal waters and gas emissions, followed by GPS measurements that indicated uplift and radial opening of Nisyros Island [83]. However, since 1998 this activity ceased and there are no more signs of unrest. In any case, the volcanic hazard is known in the area since pre-historic times and the Nisyros genesis is attributed according to the Greek mythology to the Giant Polyvotis, who has been buried below Nisyros from Poseidon. The name Polyvotis means someone creating a lot of noise, which probably corresponds to the sounds produced during the volcanic eruptions and/or major earthquakes.
Author Contributions: P.N. was involved in the field data collection, participated in results interpretation and paper writing; P.K. was involved in order to implement geospatial techniques and to interpret spatial derived data; S.K. participated by conducting data interpretation and assisted in the paper writing; D.P. and N.K. contributed to the paper editing. All authors have read and agreed to the published version of the manuscript. Acknowledgments: Seafloor mapping was a part of the GEOWARN project (Contract no: IST-1999-12310). The European Commission is acknowledged for their financial contribution to the project. ROV exploration was supported by Institute for Exploration (IFE-USA) and the collaborative project "New Frontiers in the Ocean Exploration". The authors would like to thank the officers and the crew of the R/V AEGAEO and E/V NAUTILUS for their effective contribution to the field work. This work has been [fully/partially] supported by NEANIAS, funded by the European Union's Horizon 2020 research and innovation programme, under grant agreement No 863448. Maps and diagrams throughout this work were created using ArcGIS ® software by Esri. ArcGIS ® and ArcMap™ are the intellectual property of Esri and are used herein under license. Copyright © Esri. All rights reserved.