Past and Future Impacts of the Relative Sea Level Rise on the Seafront of Ancient Delos (Cyclades, Greece) and Flooding Scenarios by 2150

1. Introduction

Rising sea levels due to global warming constitute one of the most crucial impacts of climate change, one which will dramatically affect the coastal lowlands, causing future flooding of their social, economic, cultural, and touristic infrastructures. The very insularity of its territory makes Greece one of the Mediterranean countries that will be strongly burdened by the future sea level rise. Although some studies have presented the Greek coastal areas that will, in general, be subject to flooding up to 2100 and 2150 due to sea level rise [1,2], only one has focused on the Cycladic archipelago [3], where a large number of islands are to be found (Figure 1a). This paper is centered on Delos Island, not only for its geographic position in the middle of the Cyclades in the Aegean Sea (Figure 1b), but also, indeed mainly, because of its universal archaeological value, as it was considered the most sacred of all islands (Callimachus, To Delos, 3^rd c. BC) in ancient Greek culture, and has been inscribed on the Unesco World Heritage List since 1990 [4].

Delos is situated between Mykonos and Rheneia Islands, separated from the latter by the Delos Strait, a narrow channel <1 km-wide, and two rocky islets: Mikros (Little) Rheumatiaris to the north and Megalos (Great) Rheumatiaris to the south (Fig. 1c).

Delos has an irregular oblong shape, measuring 5 km N-S by average 1.30 km E-W, and a total area of 6.85 km² (Figure 1c). It has 14.20 km of coastline, some narrow gravelly beaches, and mostly rocky steep shores with the exception of the coastal lowland on the western side, where the ancient city is located. Narrow peninsulas have formed on the northern and southern edges of the island, as have small bays on the NE (Gourna Bay), NW (Skardanas Bay), and SW (Fourni Bay) sides. Delos is mostly characterised by a rocky, hilly landscape, disrupted by small valleys. Mount Cynthus (+112 m) dominates the middle of the island, and Gamila Hill (+53 m) and Kato Vardia (+82 m) the northern and southern parts, respectively. The Inopos River originated from a basin southwest of Mount Cynthus, initially flowed in a northerly direction and then turned westward to end in the lowland where the ancient city was situated. Around 400 BC, a dam was constructed at the river bend to create a public reservoir for water collection [6].

The northern part of Delos Island consists of Miocene granitoid rocks, psammitic schist, marble, and amphibolite, which are part of the Miocene plutonic assemblage of the Cyclades [5,7]. The main rock in the central part of the island is megacrystic granodiorite, while in the south part medium-grained granodiorite and granite with abundant enclaves are the prevailing rock types. These are separated by three main shear zones: Skardanas to the north, dipping steeply to the south, Cynthus in the central part, and Fourni to the south, dipping northward at lower angles. The shear zones are hundreds of metres wide and consist of thin screens of metasedimentary rocks and parallel sheets of strained plutonic rocks striking ENE [5] (Figure 1d).

The island was first inhabited in the 3rd millennium BC on the summit of Mount Cynthus, while in the late 15th century BC the Mycenaeans settled in the valley on the western shore of the island (Figure 1c). The cult of Apollo on Delos was established by about the 9th century BC and reached a peak in the Archaic and Classical periods, when Delos became a famous religious center. In 478 BC, the Delian League of Greek cities was formed on Delos. By the end of the 3rd century BC, magnificent buildings and many sanctuaries of most of the Greek Gods and many foreign deities were constructed on Delos. In late Hellenistic times (after 166 BC), the Romans proclaimed Delos a tax-free port, which resulted in an increase in population, urban growth toward the sea, and extensive construction and commercial activity, rendering Delos the biggest hub of transit trade in the Aegean, serving all the eastern Mediterranean. Delos was sacked twice, first by Mithridates in 88 BC and a little later in 69 BC by the pirates of Athenodorus. The port became unsafe and the city fell into decline and was gradually abandoned. In the first centuries AD, there was a Christian community there, but after the 7^th century AD Delos was completely deserted [8,9]. The coastal structures (sea defences and buildings) all along the seafront of the ancient city were gradually submerged due to the relative sea level (henceforth referred to as rsl) rise since the period they were initially constructed and used, as first reported by Négris in the early 20^th century [10] and later studied in detail by relevant scholars [e.g. 11-13].

On the one hand, this paper aims at illustrating on a very high-resolution digital surface model those areas on the seafront of ancient Delos at risk of flooding, evaluating both the –mainly anthropogenic– sea level rise under different climatic projections by 2150, and the ongoing land subsidence. On the other hand, through correlation of the determined rsl changes during the last 6,300 years with ancient sea defences on the seafront of ancient Delos, a verification of future flooding scenarios is attempted, while also demonstrating the efforts of ancient societies to adapt to natural phenomena by making their structures resilient, and consequently conveying a message of effective sustainability when facing the challenges of modern climate change.

2. Materials and Methods

Rsl rise projections and flooding scenarios constitute a demanding multidisciplinary project that necessarily follow specific steps [e.g. 14]. To achieve this goal for Delos Island, which despite its local scale is certainly of particular value, we conducted the following: aerial photogrammetry to produce the Digital Surface Model (DSM) of the survey area; acquisition and processing of GNSS data for estimating the Vertical Land Motion (VLM); incorporation of geospatial and geodetic data into the updated regional sea level projections up to 2150, as provided by the Intergovernmental Panel on Climate Change's Sixth Assessment Report (IPCC AR6). At the same time, previous results on the former sea level changes were reassessed and the ancient sea defences were interpreted anew. The methodology followed for each specific step is described below.

2.1. Geospatial Surveying

Photogrammetry is based on the principles of triangulation and stereo-vision. Triangulation is the process of determining the position of a point by measuring angles to it from known points at either end of a fixed baseline. Stereo-vision, on the other hand, is the process of creating a 3D representation of an object by combining multiple 2D images taken from different angles. Today, drones equipped with cameras able to tag the images taken with GNSS coordinates have opened up new opportunities in several scientific fields, such as surveying, mapping, engineering, archaeology, etc.

Drone digital photogrammetry was applied to produce very high-resolution digital surface models (DSM) for the survey area. The aerial surveys were conducted on 10 November 2023, using a DJI Mavic 3E with global shutter camera (FOV 84ο, 24mm lens equivalent, 20MPixel, ¾ sensor size), owned by GeoSense. Image geotags have been marked utilizing post-processing kinematic (PPK) correction methodology, which resolves errors after all the geotagged located data has been collected, as RTK correction was not available.

All surveys were conducted under extremely calm wind conditions (<1 bf), in order to achieve the best shoreline reconstruction and avoid any distortion in the coastline position created by sea waves and water refraction.

Side overlap was set at 70%, frontal overlap at 80%, and flight speed at 9 m/sec. Furthermore, and despite the fact that the drone used is capable of collecting images at an even higher speed, we decided to switch to safer speeds, thus minimizing motion blur and providing crisper and more detailed images. Two successive 38-minute flights were conducted.

The relative height of the flights was 80 m with terrain follow adjustments. At this flight level, the resulting ground sample distance (GSD) was 2.2 cm/pixel, providing a very detailed and crisp orthomosaic. The geometrical accuracy of the missions is calculated by the following formulas: Horizontal accuracy = ± (GSD × 2), Vertical accuracy = ± (GSD × 3). A total of 1013 images were collected, with the surveyed area amounting to 630,000 m².

The photogrammetry software used for the image processing was PIX4Dmatic. WGS84/UTM zone 34N was selected as a generic coordinate system compatible with any other third party source of data (e.g. Google Earth/Map, OpenStreetMap), but with the flexibility to transform to local coordinate systems such as GGRS87 and HTRS08. The heights have been adjusted from geometric to orthometric by applying Geoid model EGM08.

The deliverables from the processing process were as follows: (a) Dense 3D Point cloud – colorized, in LAS/LAZ format, (b) 2.5D GeoTIFF DSM with a final resolution of 0.024 m, and (c) Orthomosaic (Figure 2).

In the absence of the official coastline position on a scale matching the needs of this study, and to precisely reconstruct the contemporary coastline position so as to apply the rsl rise projections, we carried out topogeodetic measurements with tidal correction at the time of the survey. This enabled us to define the zero-elevation position and therefore adjust the contour lines generated from DSM.

The above methodological approach assures the accuracy of geospatial data, thus achieving high reliability in applying the flooding scenarios to estimate the areas that will be potentially inundated by 2150 due to rsl rise.

2.2. Geodetic Data (VLM Rate)

The current rate of VLM was estimated using the closest GNSS station operated by METRICA on the neighbouring island of Mykonos, located in Chora of Mykonos. The station details (code, coordinates) and time interval for the data provided by METRICA are given in Table 1. Data covering the period were computed for the purposes of this study by GAMIT software at the INGV (Rome, Italy) - Geodetic Analysis Data Centre following the methodology described in [14, p.6 and relevant references therein]. According to this, from the vertical velocity for the time interval 2018-2024 (Table 1), a subsidence rate of -1.736 ± 0.752 mm yr-1 is inferred (Figure A1, Appendix A). An analogous rate at -1.8 ± 1.1 mm yr-1 was reported by [15, Table S1], calculating data for the time interval 20016-2020 from the MYKN station.

2.3. Rsl Rise Projections and Flooding Scenarios by 2150

The updated regional sea level projections up to 2150 are provided by the Intergovernmental Panel on Climate Change's Sixth Assessment Report (IPCC AR6) [16,17,18] and have been accessed through the ‘NASA/IPCC Sea Level Projections Tool’ [19]. According to this, only projections with at least medium confidence relative to the period 1995–2014 are considered for five future Shared Socioeconomic Pathway (SSP) scenarios (SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, SSP5-8.5) and five different future Global Mean Surface Temperatures (1.5°C, 2°C, 3°C, 4°C and 5°C). The IPCC AR6 sea level projections result from the consideration of several geophysical processes, indicatively: both the Antarctic and Greenland Ice Sheets, glaciers, land water storage, ocean thermal expansion, and long-term VLM inferred by the longest time series of sea level. More importantly, the IPCC AR6 clearly communicates the projection uncertainty and ambiguity of future sea level assessment [20].

Into the regional sea level projections (Mykonos and Delos Islands, Aegean Sea), we inserted the estimated current VLM rate, as computed by the time series ground displacement data obtained from the continuously operating GNSS station on Mykonos (see section 2.2.), on the proviso that the subsidence trend continues linearly until 2150 at the same rate as today’s, without being disrupted by any paroxysmal tectonic event. Our assumption is further supported by the aseismic character of the Cycladic metamorphic central core [21], with scarce and scattered seismicity [22].

Finally, the rsl rise projections were calculated for all SSP scenarios and subsequently the areas that will potentially be at risk of flooding by 2150.

2.4. Past rsl Changes in Delos Island

The systematic recording, mapping and measurement of the contemporary depth of submerged geomorphological features related to former sea levels and their dating based on archaeological sea level indicators and radiocarbon ages, throughout the coasts of the island group of Mykonos, Delos, and Rheneia, were used by [13] for the determination of five former sea levels during the Late Holocene and consequently the periods of relative sea level stability, which are consistent with the inferred former sea level stands for the northern and central Cyclades [13,23]. The rsl curve for the island group of Mykonos, Delos and Rheneia constituted the outcome of this geoarchaeological project (Figure 3). In summary, the sea level was at -4.80 ± 0.10 m between 4300 BC and 14C cal. age 1972 BC; at -3.70 ± 0.20 m between 1972 BC and 130 BC; at -2.40 ± 0.25 m between 130 BC and 14C cal. age AD 1272; at -1.55 ± 0.25 m between AD 1272 and 1650 and, finally, at 0.80 ± 0.10 m between AD 1650 and AD 1891 [13]. The sea level stands refer to depths below mean sea level (henceforth referred to as bmsl).

2.5. Mapping of the Submerged Ancient Remains and Interpretation of Ancient Adaptation Measures

The orthomosaic generated from the geospatial survey was used to demarcate the boundaries of the near-shore zone where remains of ancient structures are now submerged (Figure 2). Furthermore, previously published depths at selected points all along this zone and the functional elevations of ancient harbourworks [12,13] were taken into consideration in order to produce a detailed plan of the seafront of the ancient city of Delos and interpret the sea defences that ancient people made in order to protect the seafront of the ancient city against rising sea level, in combination with the former sea level stands, as briefly presented in section 2.4.

3. Results

Rsl Changes on the Seafront of Ancient Delos and Sea Defences during the Classical and Hellenistic Periods (5th to 1st Century BC)

The mapping of the near-shore zone, according to the methodology presented in the above section 2.5., enabled us to produce a detailed plan of the seafront of the ancient city of Delos and distinguish the geomorphological sea level indicators (beachrock generations) and submerged ancient remains (harbourworks, seawalls, protective rockfills, building remains) that are presented in Figure 4 and described below.

Along the seafront of the ancient city of Delos, north of the Main Harbour (or the so-called Sacred Harbour) and up to the southern cliff of the small Skardanas Bay, important harbourworks were undertaken, possibly during Classical Antiquity (4th to 3rd century BC) (Figure 1c,4). These consist of a rubble mound construction, an extended rockfill, stretching for 430 m parallel to the shoreline and at its southernmost end curving as it enters the Main Harbour. It has a maximum width of 115 m, a maximum thickness of 4 m, and an estimated volume of 100,000 m³. It consists of granite boulders ranging in size from 0.60 x 0.30 m to 1.80 x 2 m, resting on the rocky granite bedrock. It seems to have been built to protect the coast from waves and erosion while at the same time gaining vital space.

On top of the rubble mound harbourwork, four distinct beachrock generations have been formed, today submerged at -3.70 ± 0.20 m, -2.40 ± 0.25 m, -1.55 ± 0.25 m, -0.80 ± 0.10 m bmsl (Figure 4), the seaward end of each defining the shoreline position during the periods when each was cemented [12,13,24].

In the southern part of the harbourwork, the deepest beachrock generation (-3.70 ± 0.20 m bmsl) has incorporated a wall made of rough stones bound together with mortar; this appears to have served as a seawall to protect the busy coast of the time from coastal flooding and erosion. In front of it, at a distance of 3 m toward the open sea, a rockfill had been placed parallel to the coastline to reduce the wave action and protect the ancient coast. The rockfill is 80 m long and 40 m wide with its axis slightly curved and its surface at a depth of -3.60 m bmsl. Two almost parallel rockfills have also been placed in the central part of the harbourwork: the deepest, 100 m long and 11 m wide, and, less that one meter away, the shallowest, 50 m long and 5 m wide. The tops of the rockfills are at a depth of 3.90 m bmsl and 3.30 m bmsl, respectively, with the deepest coastline passing in between (Figure 4).

Another line of coastal protection is observed at a higher elevation toward the coast, along a shoreline, defined by the beachrock generation at 2.40 ± 0.25 m bmsl (Figure 4). It is a seawall integrated into the beachrock, 40 m long and 1.30 m wide, made of rough stones bound together with mortar. The seawall is 0.60 m high with its top surface at a depth of 1 m bmsl. A total of 16 blocks, 3 x 2 m in size, have been placed in front of the seawall at a distance of 4 m to 7 m seaward, clearly for the purpose of reducing wave energy. A 6.30 m-wide road, now submerged and visible for 30 m, seems to have ended there (Figure 4). To the south of the seawall, the protective works consist of large boulders, measuring 3.50 x 2 m, aligned for a length of 47 m in contact each other, followed by a wall made of stone blocks 0.50 x 1 m in size that stretches for a length of 53 m. To the north of the seawall, 18 protective boulders, measuring 2.50 x 2 m, had been placed along the then shoreline. Further north, although a large number of submerged Hellenistic buildings has been recorded [25,26], extending to a distance of approximately 10 m off the then shoreline, there seems to have been no provision made for sea defences/coastal protection. At the northernmost end, however, a submerged paved area approximating 900 m², demarcated by the then shoreline and resembling a dock, appears to have been protected by 25 boulders, measuring 2.50 x2 m, randomly placed about 10 m in front of it facing seaward (Figure 4). All the submerged building remains that have been integrated into the beachrock generation at 2.40 ± 0.25 m bmsl are dated back to the Hellenistic period of Delos [25,26]. The two younger beachrock generations at 1.55 ± 0.25 m and 0.80 ± 0.10 m bmsl (Figure 4) have incorporated, along with the coarse beach deposits, ruins and architectural members of the Hellenistic buildings, thus postdating their formation to subsequent periods of rsl stability [13].

3.2. Rsl Rise Projections and Flooding Scenarios by 2150

The relative sea level rise projections that combine IPCC AR6 data for all the SSP scenarios (SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, SSP5-8.5) and the inferred VLM rate for the years 2030, 2050, 2070, 2100, and 2150 were computed (as explained in sections 2.2. and 2.3.), and are presented in Table A1 (Appendix A) and Figure 5. The projected sea levels for each year and SSP scenario were applied to the DSM generated from geospatial surveying (as detailed in section 2.1.) and the areas potentially at risk of flooding were calculated and are also presented in Table A1 (Appendix A). The last data set in Table A1 (labelled: observation) represents the rsl rise projections by 2150 deduced from the rsl curve for Delos (Figure 3, section 2.4.), on the assumption that, despite the accelerating increase in global temperature, the sea level will be rising linearly during the next decades at the same rate as that of the period between the Second Industrial Revolution (late 19th - early 20th century) and 2018, for which observational data substantiate a rsl rise of 0.80 ± 0.10 m, including the actual magnitude and direction of VLMs whenever they occurred during that period.

A recent analysis of the tide gauge data from the station operated by the Hellenic Navy Hydrographic Service (HNHS) in the port of Syros (central Cyclades, Aegean Sea), performed for a time series ranging from 1969 to 2019 with data completeness reaching 86%, despite revealing a negative linear trend (slope: -7.5 mm yr^-1) for the time span 1969 – 1995, confirm a positive trend for the time span 1995 – 2019 at a rate of 3.1 mm yr^-1 /year [27]. The decreasing trend is mainly attributed to cooling up until the North Atlantic Oscillation (NAO) index shift that induced the Eastern Mediterranean Transient (EMT) phenomenon and an increase in salinity in the Aegean. The increasing trend observed afterwards can most probably be linked with the increasing sea surface temperatures [27].

4. Discussion

The continuous shoreline retreat caused by the gradual sea level rise during the late Holocene has drastically changed the coastal landscape through the ages, thus affecting human coastal activity and compelling the ancient societies to adopt protective measures against rising sea level, waves and erosion [12].

The sea level during the 5th century BC was 3.70 ± 0.20 m lower than at present [13]. At that time, sea defence works were constructed along the shoreline. During this period of rsl stability, when the deepest beachrock generation that is now submerged at -3.70 ± 0.20 m bmsl was forming, the seawall and a few ruins of that period were integrated into it. The cove formed by the then shoreline in front of the Agora of the Competaliasts, where there was a coastal swamp at the mouth of the Inopos River (Figure 4 and Figure A2, Appendix A), allowed ships to approach, disembark visitors heading to the sanctuaries, and unload goods [13,28]. The suggestion that the so-called Sacred Harbour was the Main Harbour of ancient Delos should be reassessed, as the major part of the small bay was occupied by the ancient protective harbourwork (rockfill), leaving no space for ships to approach. With a sea level at -3.70 ± 0.20 m bmsl, the coast of the Sacred (Main) Harbour would have been 200 m wider than the current coastal zone, while the cove currently being created there would have been land [13].

The recent borehole findings [29], which show that the sedimentary floor of the Sacred (Main) Harbour dated back to the end of the Hellenistic period (212-92 BC), was found at a depth of -2.25 m bmsl, that is to say about 1.50 m higher than the sea level of the Classical and Early Hellenistic period (-3.70 ± 0.20 m bmsl), and at approximately the same level of the Late Hellenistic period onward (-2.40 ± 0.25 m), with the rsl change between the two stands having occurred at some time between 180 and 130 BC [13] (Figure 3 and Figure 4).

Evidence of a rsl change by 1.30 m that shifted the sea level to -2.40 ± 0.25 m bmsl before the end of the Hellenistic period (at some time between 180 BC and 130 BC) [13] justifies the new sea defences constructed along the coast of that time (Figure 4 and Figure A2, Appendix A). This is further substantiated by epigraphic evidence, reporting that the Athenian Theophrastus, superintendent of Delos in 126/125 BC, constructed the homonymous Agora and surrounded the harbour with backfilling (“κατασκευάσαντα τὴν ἀγορὰν καὶ τὰ χώματα περιβαλόντα τῶι λιμένι”, ID 1645, as referred to in [6]). This sea level seems to have remained relatively stable even after the successive devastations of Delos in 88 BC and 69 BC. The Hellenistic ruins integrated into the beachrock generation today submerged at 2.40 ± 0.25 m bmsl prove that this beachrock was cemented after the destruction of Delos. The once busy coastal zone of the thriving ancient city is now deserted, left to the mercy of marine transgression and erosion. The sea level continued to rise gradually with three intermediate short periods of rsl stability between the 13th and 20th centuries AD before it reached its current stand (Figure 3). There was no longer any need of sea defences in the later periods as the ancient city was abandoned until the late 18th century, when the archaeological dig brought it to light again.

In the early 20th century, the debris (~60,000 m³) from the extensive archaeological excavations carried out by the French School were dumped into the Sacred Harbour and a jetty, 280 m long and max 28 m wide, was constructed there at an angle of approx. 45° to the coast, which serves the small ships coming from the old port of Mykonos (Figure 4). This human intervention dramatically altered the coastal configuration, shaping two new coves to the north and south, and causing silting in the south cove at least.

The seafloor of the near-shore zone on the western coast of Delos, where the remains of the ancient harbourworks and sea defences lie in shallow water (depths up to 4 m), slopes seaward at low gradients, 3° - 4° north of the Sacred Harbour and up to 2° to the south. The lack of fine beach deposits and sandy beaches along this coastal zone, where coarse deposits, rubble and debris prevail, mostly fosters the shoreline retreat and inland inundation to a greater extent. The Sanctuary of Apollo, temples, public spaces (agora), porticos, and other important public buildings, commercial shops, workshops, houses, etc.—generally speaking, the entire city and the port—were extensively developed along the low-lying coastal valley from the 2^nd century BC. Their ruins are located on the main route visited by the thousands of tourists who flock to this sacred island every year.

The island, due to its metamorphic geological structure, lacks natural water sources. In the fragmented granitic bedrock, a shallow ground water horizon has formed, feeding the wells that have been there since antiquity for the provision of fresh water. To collect rainwater, a dam with a reservoir was constructed uphill on Mount Cynthus across the Inopos River. Irrigation channels would direct the water flow to the low-lying valley where the city and the harbour were located. Rainwater gathering at lower elevations created the famous Sacred Lake, also favouring the growth of aquatic plants [6,8]. As a result of the rsl rise, the ground water horizon is elevated, thus shaping extended wetlands, now permanently or temporarily inundated by the brackish water, over an area estimated at approx. 1.8 ha (18,000 m²) (Figure 2). The ongoing rsl rise poses a threat to the coastal lowland, inducing extended inundation inland and creation of larger wetlands in the coming decades. To mitigate further inundation and flood damage to the archaeological site, a protection dyke was constructed by the competent Hellenic Ministry of Culture, mainly along the coast of the Sacred (Main) Harbour (Figure 2), which also serves as a dirt road connecting the modern jetty with the path used to visit the archaeological site. Additionally, some low-lying areas, mainly in front of the portico of Philip V, have been filled with rubble.

In summary, the seafront of the ancient city of Delos is subject to the effects of the rising sea level, which have caused significant coastline retreat, exposure to the northerly winds and waves, erosion, further submersion of the ancient coastal structures, and inland inundation, and is expected to be further affected by the rsl rise over the next decades. The rsl projections are, sadly, premonitory, and large areas will be at potential risk of flooding, while the modern harbour infrastructure will become non-functional. For a better understanding of the future impacts, the potential flooding areas for IPCC AR6 scenarios SSP1-2.6, SSP3-7.0, and SSP5-8.5 corrected for the VLM rate for the years 2050, 2100, and 2150 are illustrated on the orthomosaic produced and on cross-sections drawn at selected positions (Figure 6,Figure 7,Figure 9).

Indicatively, with the favourable SSP1-2.6 climatic scenario corrected for VLM, the predicted rsl rise is 0.24 m by 2050, 0.57 m by 2100, and 0.87 m by 2150 (with respect to 2020). A total area of 1.17 ha by 2050, 3.13 ha by 2100, and 4.15 ha by 2150 will have been inundated due to the rsl rise. The coastline will retreat by 190 m inland until 2150. Accordingly, the SSP3-7.0 climatic scenario provides a rsl rise of 0.26 m by 2050, 0.78 m by 2100, and 1.33 m by 2150, while the areas at risk of flooding are estimated at 1.77 ha by 2050, 3.84 ha by 2100, and 5.21 ha by 2150. The coastline will shift 210 m inland by 2150. Finally, the most unfavourable, and certainly undesired, climatic scenario SSP5-8.5 gives a rsl rise of 0.28 m by 2050, 0.87 m by 2100, and 1.48 m by 2150, while the areas at risk of flooding are estimated at 1.85 ha by 2050, 4.15 ha by 2100, and 5.45 by 2150. The coastline will shift 220 m inland by 2150 (Figure 6a,b and Figure 7a). In all cases, at least 50% of the potentially flooded areas are spread across the low-lying plain where the ruins of the ancient city are located.

What is of particular importance for validating the rsl projections and climatic scenarios are the results yielded by applying the actual rate of rsl rise that occurred in the last century, which represents the sum of the glacio-hydro-isostatic, eustatic, and tectonic components (Figure 3, Table A1) and is based on field observations and interpretation of all available geomorphological and archaeological sea level indicators [13]. The projected rsl rise is 0.23 m by 2050, 0.58 m by 2100, and 0.93 m by 2150, while the areas that will be at risk of flooding are estimated at 1.64 ha by 2050, 3.17 ha by 2100 and 4.34 ha by 2150 (Figure 7b). Comparison with the corresponding values deduced from the SSP scenarios corrected for VLM indicates that, even if none of these scenarios were to happen, the sea level would rise and its effects would impact on almost as many areas as in the moderate SSP2-4.5 scenario (Figure 8).

Regarding the modern jetty, which is the main harbour infrastructure of Delos with a maximum elevation of 1.35 m, it seems that by 2100 it will be non-functional under the scenarios SSP2-4.5, SSP3-7.0, and SSP5-8.5, albeit less affected under the more favourable scenarios.

The two cross-sections (Figure 9), A-A’ (from the shoreline ENE toward the Sanctuary of Apollo) and B-B’ (from the shoreline at the commercial district SSW landward) clearly show that, in the worst-case climatic scenario SSP5-8.5 corrected for VLM and without taking any protection/adaptation measures, both areas will be seriously flooded due to rsl rise even earlier than the end of the 21^st century, with the archaeological remains being harshly affected.

Analysis of all the data shows, beyond doubt, that the planning and implementation of protection measures for the seafront of the ancient city of Delos and adaptation solutions to the effects of the future rsl rise are now imperative. Even more so should we cogitate the consequences of storm surge events, which, although not considered in this study, could temporarily increase the reported values of rsl rise and the areas of potential flooding. Moreover, the impacts of storms and tsunamis will be intensified in sea level rise conditions, causing permanent and temporary flooding that will dramatically affect the coastal area. Additionally, increasing wave energy and wave run-up onto the western coast of Delos due to extreme storm surge events will accelerate coastal erosion and retreat [14].

Coastal areas worldwide and regionally are, undeniably, vulnerable to risk of flooding due to rising sea levels accelerated by global warming, with serious potential effects on society and the environment. Recent surveys carried out on a national scale in Greece pointed out that the awareness of Greek citizens is very low in terms of the rsl rise and consequently its potential threat is not recognised [30,31]. The results reveal a gap in the dissemination of information and knowledge among citizens, state and scientists, mainly due to the lack of centrally-organised and locally-implemented large-scale information and awareness-raising campaigns. It is important for scientists and journalists to be properly trained in the comprehensive communication of the complex aspects of the climate crisis. Scientific results concerning the potential risks of the rsl rise should be seriously considered by stakeholders and policy-makers in planning future infrastructure and imposing protection measures, resilience strategies, and adaptation solutions, along with mitigation state-policies. Digitised data and visualisation tools (apps, maps, videos etc.), exclusively based on reliable research data, will surely enhance our ability to comprehend straightforwardly why and how the climate is changing, why rising sea levels are a potential risk, why actions for a sustainable planet can no longer be delayed, why local and collective decisions and interventions on natural resources and extreme phenomena are more than ever necessary in our era [32].

5. Conclusions

Knowledge of the rsl history of a specific area is essential for a deeper understanding of, on the one hand, the regional geological and seismotectonic setting, and the extreme climatic events, and on the other hand, the interaction between humans and their coastal environment. It also contributes to the identification, interpretation and protection of the underwater traces of human existence, ensuring the preservation of cultural heritage and ecotouristic promotion of the natural seascape. Moreover, it is crucial in predicting the regional sea level rise and verifying climate modelling and flooding scenarios.

In this study, special emphasis is given to the methodology and the use of new scientific tools and techniques, which in combination with climatic models and sea level studies, amplify our perception of the extent of the past and future marine transgression and coastline retreat induced by the former and ongoing rsl rise all along the seafront of the ancient city of Delos. This approach is crucial for the monitoring and future management of the flooding risk of one of the most prominent monumental sites of our cultural heritage, one which has been seriously affected and will be further impacted by rising sea levels. Beyond methodological and scientific aspects, the results of this study will be decisive in the design of early response systems to coastal retreat, future coastal and harbour infrastructure, and civil protection and adaptation measures, mitigating the catastrophic effects of floods and preserving both our cultural heritage and ecosystems.