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 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 established 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 centred on Delos Island, not only due to 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, 3rd 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 (
Figure 1c).
Delos has an irregular oblong shape, measuring 5 km N–S by an average of 1.30 km E–W, with a total area of 6.85 km
2 (
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) dominate 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 southern 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 centre. 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 of 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 few centuries AD, there was a Christian community there, but after the 7th 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 rising relative sea level (henceforth referred to as rsl) since the period they were initially constructed and used, as first reported by Négris in the early 20th century [
10] and later studied in detail by relevant scholars (e.g., [
11,
12,
13]).
On the one hand, this paper aims to illustrate on very high-resolution digital surface models 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 6300 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 activities: aerial photogrammetry to produce a digital surface model (DSM) of the survey area; the acquisition and processing of GNSS data to estimate the Vertical Land Motion (VLM); and the 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 stereovision. 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. Stereovision, 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°, 24 mm lens equivalent, 20MPixel, ¾ sensor size), owned by GeoSense. Image geotags were marked utilizing the post-processing kinematic (PPK) correction methodology, which resolved errors after all the geotagged located data had 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 was set at 80%, and flight speed at 9 m/s. 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 min flights were conducted.
The relative height of the flights was 80 m with terrain-following 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 was calculated according to 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 m2.
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 were adjusted from geometric to orthometric by applying Geoid model EGM08.
The deliverables from the processing process were as follows: (a) a dense 3D point cloud, colorised, in LAS/LAZ format; (b) a 2.5D GeoTIFF DSM with a final resolution of 0.024 m; and (c) an 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 ensures 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 y
−1 was inferred (
Figure A1,
Appendix A). An analogous rate at −1.8 ± 1.1 mm y
−1 was reported by [
15] (Table S1), calculating GPS 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 were provided by the Intergovernmental Panel on Climate Change’s Sixth Assessment Report (IPCC AR6) [
16,
17,
18] and were accessed through the ‘NASA/IPCC Sea Level Projections Tool’ [
19]. According to this, only projections with at least medium confidence relative to the period of 1995–2014 were 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 resulted from the consideration of several geophysical processes, including 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, without being disrupted by any paroxysmal tectonic event. Our assumption is further supported by the aseismic character of the Cycladic metamorphic central core, with scarce and scattered seismicity [
13,
21,
22,
23]. The ‘aseismically deformed’ area of the northern and central Cyclades is probably the result of the existence of a very closely spaced geometric fracture framework within the metamorphic rocks, preventing strain accumulation. Thus, energy release manifests itself in continuous deformations creeping along the fracture planes [
23] (p. 24).
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 were established.
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 levels, in combination with the former sea level stands, as briefly presented in
Section 2.4.
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 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 A3,
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 an 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 A3,
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 for sea defences in the later periods as the ancient city was abandoned until the late 18th century, when an archaeological dig brought it to light again.
In the early 20th century, the debris (~60,000 m
3) 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 2nd 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 groundwater 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
2) (
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 coming 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 orthomosaics produced and on cross-sections drawn at selected positions (
Figure 6a,b,
Figure 7a and
Figure 8).
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 an 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 an 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 9).
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 8), A-A’ (from the shoreline ENE toward the Sanctuary of Apollo) and B-B’ (from the shoreline toward 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 21st 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. Furthermore, we should 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 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, and why local and collective decisions and interventions on natural resources and extreme phenomena are more necessary than ever in our era [
32].