1. Introduction
In recent decades, the coastal zone has been continuously sustaining increased tourist activities, but it is also directly impacted by climate-change factors (i.e., continuous sea-level rise, increased frequency and magnitude of storms) as well as non-climate change related factors (i.e., land use, population increase, water management, subsidence) [
1]. Climate change, particularly associated with the ongoing sea-level rise, is intensifying coastal erosion and inundation, resulting in significant environmental changes, habitat loss (directly or indirectly), as well as significant financial loss through its effects on local populations and tourism in the area. The consequences of natural phenomena, such as tsunamis, hurricanes, and storm surges, are greatly magnified due to increasing population densities at the coast. On the other hand, increasing human pressure in the coastal zone is reducing natural buffer zones and causing the so called “coastal squeeze” where valuable habitats, coastal processes, coastal ecosystem services, etc. are all squeezed into narrow zones. The observed climate change and the resultant sea level rise, the increase of storm surges and other extreme events that afflict the coastal zone and, also, the rapidly increasing population of the coastal zone will surely have a significant impact not only on the resilience of the coastal systems, but also on the strategies of adaptive services and goods [
2,
3,
4].
In the case of the Mediterranean region, approximately 70% of the countries’ population lives in the coastal areas [
5]. The Mediterranean coastal zone hosts a notable number of visitors, with increasing demand for access to infrastructure and accommodation. This concentrated urbanization in combination with industrial development has led to significant pressures on the coastal environment, including habitat loss, degradation and fragmentation and hence loss of biodiversity. In the absence of proper management, these pressures threaten the values that have led to such visitation numbers and development.
Sandy coastal ecosystems, lying in the transition zone between marine and sedimentary terrestrial environments, are subject to continuous changes in their morphology and vegetation [
6,
7], while showing great biodiversity, in terms of plant species and communities [
8]. These fragile ecosystems hosting beaches, wetlands, and dune habitats are threatened by coastal erosion and human impact [
9]. Coastal dunes constitute an unstable, but dynamic geomorphological unit of the coastal zone. Sand dune systems hold the dynamic role in the coastal zone, as they constitute a natural barrier protecting the inland from eroding wave action, while they also act as water filters, enhancing the quality of coastal water. Sand dune habitats provide niches for plants and animals, with many rare and often endangered species. Nowadays, the loss and degradation of the coastal dunes is mainly owed to human activities [
10], which typically differentiate the natural processes on the coastal zone and therefore modify the system dynamics. These activities may be categorized into six sectors [
10,
11]: (a) housing and recreation, (b) industrial and commercial use, (c) waste disposal, (d) agriculture, (e) mining, and (f) military activities.
At the same time, coastal wetlands are very important and valued features due to their high biodiversity and ecosystem services, such as preventing floods, cleansing waters, and recharging groundwater aquifers. Wetlands also host a wide variety of flora and fauna, providing habitat for rare and endangered species. Many of the Mediterranean coastal wetlands are threatened by a number of factors, such as modifications of their natural functions, over-fishing and hunting, excessive pumping, and partial drainage [
12]. Hence, these coastal environments are particularly vulnerable and worthy of attention and protection for their specific ecosystem functions.
In this context, our study aims to assess the vulnerability and pressures on the western coasts of Naxos, considering both natural and anthropogenic factors. This is achieved through a series of methodologies that include (a) the study of aerial photographs from 1960 until today, (b) systematic seasonal aerial monitoring by drone since 2015, (c) coastal zone mapping and topographic sections through the use of DGPS, and (d) mapping of coastal dunes.
2. Study Area
Naxos Island is located in the central part of the Aegean Sea, 200 km SE from Athens (
Figure 1). The island has a mountainous relief, and the midland is dominated by a mountain chain from the north to the south across the island. The western part of Naxos is characterized by low relief, with sandy beaches along an extensive coastal zone and, consistently, concentrates the majority of the communities in the vicinity of our study area. Typically, the western coasts are composed of a sandy beach, bordered by low lying sand dunes, lagoons and an alluvial plain. Dune-fields are well-developed along the western coast, with an almost continuous extent from Chora (Naxos city) to Pyrgaki, 20 km further south. These dunes were already mapped both during the realization of the geological map of Naxos [
13] and in recent research focusing on coastal evolution [
14]. According to past field investigations [
15], six main dune fields are found (
Figure 1), commonly separated transversal by granodiorite promontories, jutting out at the coastal area.
The wind, which plays a key role in the formation and evolution of the dune fields, is characterized by a predominance (>60%) of northerly winds (Etesian Meltemi) [
16]. The southerly winds represent about 25%, while the westerly and easterly winds occur only occasionally. The trend is similar when considering the storm winds (>7 on Beaufort scale). According to 10 years of data (1995–2004), on an annual basis, the most common propagating ripples originate from the 0° and 225° areas, with an incidence of 38% and 25% respectively (
Figure 2a). The average value of the significant wave height is 0.73 m, with typical values between 0–0.5 m, while the maximum value of the significant wave height is 5.6 m with direction N–NW (area 352°) (
Figure 2b). The highest significant wave heights originate from the N and SW directions (area 345–10° and 225°) [
17]. With respect to fetch length, it extends to 40 km for the northern direction, (extending up to Mykonos Island), 55 km for the southeastern direction (extending up to Sikinos and Folegandros islands), and to 30 km for the southern direction, limited by Ios island.
Late Holocene Evolution of the Western Coasts of Naxos
On the western coasts of Naxos, several landforms testify to the late Holocene evolution of the area and the relative sea level changes. A series of slabs of submerged beachrocks have been identified at five different sites at western Naxos, reaching depths of −6.3 m, providing evidence on the morphological configuration of the coastline and the coastal evolution in the last few millennia [
18].
At Agios Georgios, three boreholes have been accomplished, providing insights into the paleogeography of the bay for the last 6000 years [
19]. Around 6000 BP, the bay of St. Georgios was protected by an elongated beachrocks reef, which extends almost parallel to the modern coastline. For the period between 6144 BP until 232 BP, St. Georgios was an active lagoon. Sedimentological and palaeontological data have shown that the embayment was changing from a coastal environment to an environment frequently alternating from shallow marine with freshwater input to a brackish mesohaline one [
20]. Parts of the beachrocks reef, during that period were periodically submerged and sea water would enter into the bay [
19].
A few kilometers to the south, the Mikri Vigla area was also protected by a beachrock bench. Analysis of sediments and micro-faunal content from three boreholes revealed that between 3800 ± 50 BP to 1625 ± 45 BP, the site was an active lagoon [
19]. The deepest beachrock partially protected the lagoon 2000 years ago, while the shallower beachrock was partially protecting the area up to 1000 AD [
21]. At least during the Bronze Age (3300–1200 BC), the site of Vigla may have been used as a place where boats could anchor, even though bound for other destinations [
19].
Fossil palaeoshorelines in the form of tidal notches have been identified at Naxos island, as well as other Cycladic islands. At least seven former shorelines were found at depths between 280 ± 20 and 30 ± 5 cm below the present sea-level [
22]. The vertical sequence of several submerged notches indicates rapid subsidence events, probably of seismic origin that have taken place after 3300 BP, suggesting a subsidence rate of about 0.8–0.9 mm per year [
22].
The presence of the submerged beachrocks [
18,
23], the coastal corings from western Naxos [
19,
20] along with the multiple submerged notches testify to relative sea level rise during the late Holocene. Similar analysis from the neighboring island of Paros [
24] and comparisons with sea-level predictions from two GIA models suggest an average rate of tectonic subsidence ~1.0 ± 0.4 mm/yr since 5500 cal BP [
24]. This subsidence trend in the central Cyclades is most likely owed to gradual long-term subsidence [
24,
25,
26] in combination with seismic events [
22], due to the dominance of an extensional structural pattern [
27].
3. Methodology
For this work, we used a series of methodologies for the environmental and geomorphological study of the western coastal zone and the shallow submarine area of Naxos, which included: (a) photointerpretation of aerial photographs from 1960 until today. The aerial photographs from 1960 and 1988 were obtained from the Hellenic Military Geographical Service, were photointerpreted through stereoscopic observation and were used for the documentation of coastal features; (b) systematic seasonal aerial monitoring by drone, since 2015, (c) mapping of the coastal zone and topographic sections using DGPS (Differential Global Positioning System). Specifically, a DGPS-GNSS Spectra SP60 was used for this purpose, which offers an accuracy of 5 cm; (d) dune mapping, which was also accomplished using DGPS.
Dune-fields are well-developed along Naxos western coast, with an almost continuous extent from Chora (Naxos city) to Pyrgaki, 20 km further south. These dunes were already mapped both during the realization of the geological map of Naxos [
13], the geomorphological map [
28] and in recent research focusing on the coastal evolution [
19,
28]. In this study, we performed detailed topographic mapping of the coastal dunes, using a DGPS-GNSS Spectra SP60 with accuracy of 5 cm. We performed cross sections on the dune-fields defining the back-shore boundary, the boundaries of the sand dunes, the vegetation cover, the coastline, and the corresponding sea level. The cross sections were performed on the dune-fields of Agios Prokopios, Plaka, Mikri Vigla, Glyfada, and Pyrgaki (
Figure 3). The sand dunes are named according to their proximity to the coastline. Thus, the closest to the coastline series of dune fields were coded as T1 and the more remote are coded ascendingly (T2, T3, etc.).
4. Results
4.1. Agios Georgios Tombolo
The development of Agios Georgios tombolo was studied using aerial photographs from 1960 until today. From the study of the aerial photographs of 1960, 1988 and 2015 the tombolo shape has changed. In 1960 the small island of Manto was connected to the rocky islet on the west. Beachrocks are also present on the north shore of Manto. In 1960 the islet’s distance from the shore was 22.90 m and covered an area of ~3230.2 m
2. During 1988, its morphology has changed and the distance from the shore was decreased to 16.60 m and its covering area was reduced to ~3000.6 m
2 (
Figure 4a). Additionally, it can be observed that a narrow land that was once connecting Manto and the rocky islet was submerged. During fieldwork in 1997–2000 [
28], the rocky islet was still submerged along with three beachrock slabs. Today, the tombolo is formed between the Manto islet and the main coast of the island and it is even possible to have access by foot (
Figure 4b). The area of Manto has been reduced between 1960 and 1988 by 53.47% [
28]. Nowadays, the area reaches 25,861.7 m
2.
Human driven changes in the area are also noted in the aerial photograph of 1988. A coastal road has been built near the shore as well as and the construction of Naxos airport. Today, many human activities affect the area, however they seem to have no impact on the tombolo.
4.2. Agios Prokopios
The coastal zone of Agios Prokopios has a length of 1100 m and a width ranging between 15 m and 124 m, while the coast has a gentle dip of 8°. The coastal area is characterized by the presence of sand dunes and a lagoon. The dunes are extensive and cover an area equal to 2/3 of the beach. This area, during the summer months, is used as a parking space. Behind the dunes there is a lagoonal system, parallel to the coastline (
Figure 5a). During summer, when the water table is lower, the lagoon dries up and it is separated into three lagoons. The presence of sand dunes also functions as a natural barrier to the water communication of the three seasonal lagoons with the sea. The coastal zone of Agios Prokopios is partially delimited by the arbitrary and irregular road network that has been developed on the dunes, which limits the width of the land beach to 15 m and by a sequence of low dunes with vegetation.
The dune ridge is located only in the northern half of Agios Prokopios beach, in contrast to the southern part, which is replaced by dense urbanization. The presence of anthropogenic pressure at the specific area is strong despite the attempts of preservation, which seem insufficient. The composition of the sand dunes is characterized as sand, which is provided by the eroded local granodiorite.
The performed cross section is located almost at the middle of the three-lagoon system (
Figure 5b). The T1 dune has a width of 40 m with a maximum elevation of 2.48 m. This dune is characterized by vegetation on its landward part which is protected by the sea and sea spray. The T2 dune has a width of 43 m with a maximum elevation of 2.3 m while it is fully covered by vegetation. Between the T2 and T1 dune, there is a manmade unconformity serving as a road for light weight vehicles. At the back border of the T2, there is a coastal road which separates the dune field. The last T3 dune has width of 20 m and elevation of 0.55 m. The three-lagoon system is located behind the dunes with a width of 102 m and the deepest part reaches 9 cm below sea level, which is located 240 m from the shoreline and 12 m from the back border of the lagoon.
4.3. Plaka
The coastal zone of Plaka area is characterized by well-developed dunes and seasonal streams appear at the southern part (
Figure 6a–c). The shape of the sand dune systems is rectilinear and changes to an arc shape on the southern part, characterized by complexity and two series (T1, T2). Their length is located for about 3 km with a width of 15 m. The presence of beachrocks in front of the beach leads to the deposition of the finest parts of sediments, as the beachrocks prevent the transport of coarser sediments by the wave. Therefore, the wind transports the finer grains more easily and forms, in a shorter time, higher sand dunes without vegetation. The fine-grain material of which the sand dunes are composed leads to the erosion of their lower layers. It is noted that the area of Plaka shows absence of any lagoonal features.
Due to the extended Plaka dune-fields and their morphological variety, we created four cross sections: Plaka 1 (
Figure 6a,d), Plaka 2 (
Figure 6a,e), Plaka 3 (
Figure 6b,f), and Plaka 4 (
Figure 6c,g).
On the cross section, Plaka 1 is located at the north part of Plaka beach. The T1 dune has a width of 9 m with a maximum elevation of 3.9 m. Following the T2 dune has a width of 15 m with a maximum elevation of 6.8 m. Both dune series are stabilized by vegetation. The human impact is observable on the shore as there are organized public sun loungers. Furthermore, a coastal road is cutting through the T2 dune, which is possibly to be extended even further to its landward part.
The cross section Plaka 2, is southern to Plaka 1 in a distance of 45 m. This cross section reveals the complexity of the dune formation as the two series are not easily distinguished. It is believed that human impact may have done alternations on the dune morphology. The dune has a width of 24 m forming two ridges with a maximum elevation of 2.2 m above sea level and fully covered by vegetation. The dune is cut off by the aforementioned coastal road and continues with the T2 which has a width of 32 m forming multi-ridges with the highest one at 4.9 m. The T2 dune is fully covered by vegetation. At the back boundary of T2 dune, many human alternations can be observed as ongoing contractions.
The Plaka 3 cross section is located southern to Plaka 2 at a distance of 400 m. Dune T1 has a width of 7 m and a maxim elevation 2 m, while T2 has a width of 30 m with the highest elevation at 5.8 m. A parking and the coastal road are noted in the landward part of T2. Furthermore, intense human activity is present, such as houses, stores and hotels.
Finally, Plaka 4 cross section is located 1 km southern to Plaka 3. On this cross section, only a T1 dune can be observed with width of 48 m. The dune is fully covered by vegetation and its back boarders the coastal road.
4.4. Mikri Vigla
The coast of Mikri Vigla has a total area of 2 km
2 and features two natural coastal lagoons clustered on a 0.2 km stretch of the coast, extensive transitional sand dune systems and a temporary pond formed between the lagoons (
Figure 7a). Also, the planting of trees by humans in the sand dunes reflects the desire to protect the dunes or houses near the dunes, especially at the eastern part of the study area, where human constructions were built directly on the dunes. The dunes’ shape is semicircular, with a total length of 678 m and a width 20 m.
Three cross sections were performed at Mikri Vigla area (
Figure 7b–d). The first cross section, Mikri Vigla 1 (
Figure 7b), was accomplished on Parthena beach, north of Mikri Vigla village. This cross section indicates two series (T1/T2) of dunes. The T1 dune has a width of 6 m and a maximum elevation at 3.7 m. T2 dune follows with a width of 24 m forming two ridges with maximum elevation at 6.8 m. Both dune series are covered by vegetation.
At Mikri Vigla beach, the characteristic cross section, Mikri Vigla 2 (
Figure 7c), indicates the two-series dune system protecting the lagoonal system. T1 dune has a width of 55 m with a maximum elevation at 2.7 m. T2 dune has a width of 25 m but with a maximum elevation at 0.88 m. Both dunes are fully covered by vegetation. The lagoon extends behind T2 dune for 170 m and has the same elevation as the sea level.
The third cross section, Mikri Vigla 3 (
Figure 7d), is located 500 m south of Mikri Vigla 2. At this cross section, three series of dunes are noted. T1 is characterized by a width of 10 m and maximum elevation at 2.6 m, while T2, in close proximity to T1, has a width of 6 m and 1.9 m elevation. At the landward part of T2, a coastal road can be observed. The T3 dune has a width of 10 m and an elevation at 3.4 m. All dunes are covered by vegetation.
4.5. Glyfada
The coastal area of Glyfada is characterized by extensive sand dunes, a seasonal lagoon while and the beach slope is about 12° (
Figure 8a). The total area covered by dunes is 370 m2 and features two natural coastal lagoons clustered on a 0.13 km stretch east of Glyfada beach. There is a communication channel at the southern part of the lagoon where the dune-fields end. The lagoons are isolated from the sea by extensive sand dune systems, at a distance of about 9.40 m from the coastline. Their shape is defined as a circular arc with a total length of 2 km and width of 20 m until the beach.
The typical geomorphological features of Glyfada coastal area are interpreted on a cross section (
Figure 8b). The Glyfada dune has 62 m width forming multiple ridges of 4.75 m maximum elevation. Behind the dune, the lagoon extends for 77 m with an elevation of 5 cm below sea level.
4.6. Pyrgaki
The site of Pyrgaki has a total area of 0.12 km
2 and the sediment composition of the beach consists of sand mixed with pebbles (
Figure 9a). The presence of pebbles can be explained by the presence of an old stream, which seems to be no longer active. It is noted that the dunes cover them.
The Pyrgaki cross section indicates an extended two series dune field (
Figure 9b). The T1 dune has 20 m width and an elevation of 4.9 m. The T2 dune has a width of more than 80 m and a maximum elevation of 10.1 m. This dune is disturbed by a road. All dunes are covered by vegetation while the T2 is fully covered by a cinder field, which creates a unique scenery.
5. Discussion
A number of coastal and marine landforms testify to the late Holocene and present-day evolution of the area and the relative sea level changes, providing evidence concerning the morphological configuration of the coastline and the coastal evolution. The coastal zone of the western coasts of Naxos is typically composed of a sandy beach, bordered by low lying sand dunes, lagoons and an alluvial plain, while the submarine area is characterized by the presence of a series of beachrock slabs [
18,
19,
20,
28] (
Table 1). From North to South, dune fields are developed, at an almost continuous area of 20 km at the sites of Agios Prokopios, Plaka, Mikri Vigla, Kastraki, Glyfada, and Pyrgaki, defining the back-shore boundary, relative to the coastline. At the same time, natural coastal wetlands (lagoons) extend behind the dunes, inland, below sea level. Sand dunes are one of the most vulnerable and threatened ecosystems of the island. The sand dunes have average dimensions 3–10 m in height and 30–150 m in width, showing high slopes on the windward and mainly on the leeward side, while the maximum value of the significant coastal wave height is 5.6 m with N –NW direction. Based on the distribution of the dune fields, it seems that the relief, the morphological characteristics of the area, along with the sand source supply have played the major role on their location. On the other hand, the influence of the wind can be recognized in the constant NNW-SSE to N-S orientation of the dune fields. The natural coastal lagoon system extends behind the sand dunes, for 77–170 m, with the deepest part reaching 9 cm below sea level.
These systems are becoming increasingly vulnerable, due to sea level rise, considering both natural (waves/storms/extreme weather events) and anthropogenic pressures, leads to their erosion and degradation. Anthropogenic pressures include uncontrolled parking, vegetation uprooting, heavy touristic traffic, etc. The study of past aerial photographs reveals that the main pressures on the coastal zone derive from the recent human interventions and the development of touristic and residential constructions. Specifically, at Agios Prokopios area, the development of sand dunes has been affected by human activity, mainly at the southern part of the coast as they have been replaced by dense urbanization, exposing the area to extreme natural phenomena. At Plaka and Pyrgaki, anthropogenic pressures located at the back of sand dunes (e.g., parking area, coastal roads) reduce their width and make them prone to erosion, in combination also with the reduction of clastic feeding from the seasonal streams, making these conditions irreversible. Mikri Vigla is an example of what happens when human removes and degrades sand dunes, as after human constructions built directly on sand dunes, they were forced to take protection measures by planting trees in the sand dunes, reflecting the desire to protect the dunes and constructions near the dunes, especially from the southeastern natural events. The lagoon of Agios Georgios and the seasonal lagoons of Glyfada, Mikri Vigla and Agios Prokopios are characterized as significant habitat areas both in terms of biodiversity and population as well as scientific importance. Human intervention in these areas (construction of an airport and creation of cultural fields at Agios Georgios, and a parking area and touristic activities at Glyfada, Mikri Vigla, and Agios Prokopios) led to their degradation, causing significant problems in the biodiversity of the areas and consequently the environment.
The location and morphology of the dune systems make them important variables in the coastal zone and specifically in the advance of seawater in a possible sea level rise, because of global climate change. Sand dunes act as a natural barrier that protects the inland from strong natural phenomena, preventing the penetration of water bodies. In case of sea level rise only by 0.5 m, according to the conservative scenario of the IPCC (2007, 2013, and 2019) [
29,
30,
31] for the next 100 years, it is estimated that the land to be covered by the sea will be limited in area, due to the high altitudes and slopes of the sand dunes, which contribute to the area as a measure to prevent coastal flooding [
32]. Moreover, sand dunes act as filters for water, enhancing and maintaining coastal water quality, while provide niches for highly specialized plants and animals, including numerous rare species and some endangered species. At the same time, coastal wetlands are very important and valued features due to their high biodiversity and ecosystem services.
The combined tectonic subsidence in the study area (submerged beachrocks/notches) along with the ongoing global climate change particularly associated with the sea-level rise, poses a threat to a variety of coastal landforms as well as on several touristic activities that are hosted on the western coasts of Naxos. In the absence of sand dunes, marine flooding will follow not only to the beach, but also to the lagoons and many acres of rural and residential land, resulting in significant environmental changes and financial loss through its effects on local populations and tourism in the area. In addition, sand dunes can also protect from coastal hazards, such as tsunamis, as they can partially absorb tsunami wave energy [
33], and in fact the dunes fields of Glyfada offer better protection in comparison to those of Agios Prokopios, as the former have been less impacted by anthropogenic pressures [
33].
However, sand dunes of the specific areas are subject to significant environmental pressures, direct and indirect, mainly from human activity, exposing significant parts of coastal areas to the risk of flooding by sea water. In particular, the effects of uncontrolled tourism development (e.g., irregular road network and passages on the sand dunes, sand extractions and clearings), in combination with the reduced supply of material for the development of the sand dunes, results in their leveling and consequently in their destruction. Long-term loss of sand from the beach system will lead to coastal retreat, with the outgoing amount of sand being greater than the incoming. At the same time, the presence of the underwater beachrocks along the shallow submarine area of the western coasts, which extend from many kilometers parallel to the present-day coastline, act as natural breakwaters, protecting the coastal zone from erosion [
34]. Their presence and preservation are essential for the protection of the coastal zone. The potential loss of sand dunes is a threat to the wider areas, especially the low-lying coastal areas, from marine floods, resulting in the social, cultural, and economic degradation of the areas.
6. Conclusions
The coastal areas of western Naxos are characterized by several landforms such as wetlands that border with low-lying sand dunes. The coastal zone is typically composed of a sandy beach, bordered by low lying sand dunes, lagoons, and an alluvial plain. These systems are increasingly vulnerable due to sea level rises, considering both natural (waves/storms/extreme weather events) and anthropogenic pressures, which lead to their erosion and degradation. The main pressures on the western coasts of Naxos are directly related to the effects of uncontrolled tourism development (e.g., irregular road network and passages on the sand dunes, sand extractions, and clearings). The coastal wetlands of the study area are also burdened with human interventions, impacting their biodiversity. Regarding the sand dunes, these pressures in combination with the reduced supply of material, results in their leveling and consequently in their destruction. Therefore, the protection of sand dunes, their restoration where they have been eroded, and the definition of a zone of control and restriction of human activities must be a key priority. The existence of sand dunes and a sufficient coastal width (determined by coastal geomorphological-oceanographic study) represent a necessary condition for limiting the inflow of seawater.