Next Article in Journal
Construction of Ecological Network Based on Multi-Scale Conversion and Nesting
Previous Article in Journal
Detection of Spatial Shift in Flood Regime of the Kabul River Basin in Pakistan, Causes, Challenges, and Opportunities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrating Geomorphological Data, Geochronology and Archaeological Evidence for Coastal Landscape Reconstruction, the Case of Ammoudara Beach, Crete

by
George Alexandrakis
*,
Stelios Petrakis
and
Nikolaos A. Kampanis
Coastal and Marine Research Lab, Institute of Applied & Computational Mathematics, Foundation for Research & Technology-Hellas, 70013 Heraklion, Greece
*
Author to whom correspondence should be addressed.
Water 2021, 13(9), 1269; https://doi.org/10.3390/w13091269
Submission received: 12 April 2021 / Revised: 27 April 2021 / Accepted: 27 April 2021 / Published: 30 April 2021

Abstract

:
Understanding the processes that govern the transformation of the landscape through time is essential for exploring the evolution of a coastal area. Coastal landscapes are dynamic sites, with their evolution strongly linked with waves and sea level variations. Geomorphological features in the coastal area, such as beachrock formations and dune fields, can function as indicators of the coastal landscape evolution through time. However, our knowledge of the chronological framework of coastal deposits in the Aegean coasts is limited. Optically Stimulated Luminescence dating techniques are deemed to be very promising in direct dating of the coastal sediments, especially when they are linked with archaeological evidence. The dating of the sediments from different sediment core depths, determined by the method of luminosity, allowed us to calculate the rate of sediment deposition over time. More recent coastal evolution and stability were examined from 1945 to 2020 with the use of aerial photographs and satellite images. This paper presents the 6000 ka evolution of a coastal landscape based on geomorphological, archaeological, and radio-chronological data. Based on the results, the early stages of the Ammoudara beach dune field appears to have been formed ~9.0–9.6 ka BP, while the OSL ages from 6 m depth represented the timing of its stabilization (OSL ages ~5–6 ka). This indicates that the dune field appears to already have been formed long before the Bronze Age (5–10 ka BP). It became stabilized with only localized episodes of dune reactivation occurring. In contrast, while high coastal erosion rates were calculated for modern times.

1. Introduction

Understanding and modelling geomorphological and depositional processes, resulting in the transformation of the landscape from a diachronic perspective, are essential for exploring the cultural biography of a coastal area. Coastal landscapes are dynamic landscapes, with their evolution strongly linked with waves and sea level variations. There are several challenges in the study of the cultural landscape of coastal areas. In terms of environmental change, these challenges are principally related to the dynamic character of coastal areas. Wave and wind action result in processes such as coastal erosion and the formation of coastal dunes. Coastal erosion often leads to an irreparable loss of information used in the identification of maritime structures and activities. The Mediterranean, as a major example of a strictly coastal related inhabited area since the ancient times, is a hotspot of both opportunity and concern. Consequently, these changes in the coastline can generate new social relations as newly formed areas can be suitable for agriculture or already inhabited ones can be abandoned due to flood by seawater [1].

1.1. Tectonics of Crete

The island of Crete is located along the transition zone between the African and Eurasian plates in the fore arc of the Hellenic Subduction Zone. Differential tectonic behavior is observed in the continental area of Crete. Compression prevails in western Crete, while tectonic extension dominates in eastern Crete [2,3]. By the late Upper Miocene to early Pliocene, the conditions of deposition changed from open sea to coastal or subaerial environment [2,3]. The island uplifted gradually during the Quaternary, as evidenced by marine Pleistocene terraces and coastal deposits formed during the corresponding interglacial stages [4,5,6,7,8,9,10]. The fragmentation and uplift of the island continued during the Pleistocene [6]. During the Upper Holocene, subsidence affected the entire island of Crete and tectonic activity produced intense faulting [11,12,13,14]. Three subsequent phases of tectonic activity formed initially E–W-trending faults that affected Miocene sediments. Next, N–S-striking faults cut the E–W faults which had previously formed, while the younger phase included NE–SW- and NW–SE-dipping faults [15,16]. Arc-parallel (E–W) and arc-normal (N–S) extension faults are related to the Holocene tectonic activity and are possibly connected with large historical earthquakes and recent moderate seismicity (M = 6.0–6.5) [17,18,19,20]. An evolution scenario for the coast of central and eastern Crete [21] indicated that the eastern part of Crete has subsided gradually toward NE to a depth greater than 4.0 m in at least three subsiding phases over the last 4000 years [14].
Moreover, the eruption of Santorini in the late Bronze Age resulted in tsunami waves, based on the evidence of tsunami deposits along various coastal areas of Greece [22,23]. Numerical models have shown that tsunami wave amplitudes were significant, with inland inundation of up to 200 m along the coast of Crete [24]. More recently, evidence from the AD 365 earthquake, the largest known seismic event in the Eastern Mediterranean, has supported the idea that the respective tsunami not only reached the NW part of Crete but also severely affected the inhabited coastal areas in Heraklion [25], where vertical movements in the area reached −1 m [26]. The north coast of the island appears to be rather stable after the AD 365 earthquake event [27] (Figure 1).

1.2. Indicators of Sea Level Changes on Crete

Geomorphological features such as tidal notches, beachrock formations, and archeological evidence along the coastline of Crete are evidence of relative sea level changes over time (during the last 5000 years). In central and eastern Crete, the relative sea level change evidence from tidal notches and beachrock formations revealed five distinct sea level stands at −6.55 ± 0.55 m between 4200 ± 90 BP and 3930 ± 90 BP, −3.95 ± 0.35 m between 3900–3700 BP, −2.70 ± 0.15 m between 3600–3450 BP, −1.25 ± 0.05 m, and −0.55 ± 0.05 m. The lowest sea level stand can be identified by the oldest dated tidal notch of western Crete. Historical sources report a relative sea level rise of 0.70 m during the AD 1604 paroxysmal event. Over the last 400 years, the relative sea level rose by 0.55 m [14].

1.3. Beachrocks and Aeolian Deposits

Beachrock formations originate from the rapid cohesiveness of coastal loose sediments (sands, gravels) in the inland beach area. Their formation takes place either superficially and/or under a thin layer of sediments in the subtidal zone and includes sediments of various diameters and of classical and/or biogenic origin [28,29,30]. The dimensions of the beachrocks range from a few meters to a few kilometers in length, with a width from less than 1 meter to 300–400 m and thickness from 0.3 to 3 m. These formations present flat layers that slope toward the sea, either with the general slope of the beach [31] or with a greater slope than that of the rest of the beach [32,33,34]. Appearances have been reported on many coasts of the world [35,36,37] and in Greece [38,39]. Because they form in the tidal zone, beachrock formations are considered indicators of old coastlines and are often used as a fully accepted and valid method by the scientific community in research to detect sea level changes in earlier periods [36,40,41]. Beachrock formations in Crete are composed of sand and gravels (quartz, limestone lithoclasts, metamorphic, and volcanic) with generally less than 10% of bioclasts. The intertidal cement is cryptocrystalline to microcrystalline, mainly Mg-calcite and in some cases pure calcite [14,42,43]. The determination of the age and position of the respective past coastline and the corresponding sea level by the depth of occurrence of a beachrock formation and its age, using stratigraphic or radio-chronological methods, is now widely used [44,45,46,47,48]. On the other hand, aeolianite is any rock formed by the lithification of sediment deposited by aeolian processes. Coastal aeolianites are usually coastal formations consisting of carbonate sediment of shallow marine biogenic origin, formed into coastal dunes by the wind and subsequently lithified [49]. Quaternary aeolianites are typically coastal deposits composed of fine-to-medium-grained, well-sorted sand. The character of the sand grains significantly depends on the local environmental setting, although the dominant constituents of the aeolianites are quartz and feldspar grains and marine carbonate particles. Both the detrital mineral grains and the carbonate particles are often derived from local sources. Several studies have used dating methods to establish the time framework of aeolianites and beachrock formation and to assess their connection with sea level changes [50,51,52,53]. Aeolianites are mainly considered to be deposited around the world at sea level high stands [54,55], but several recent studies have indicated the cold and low sea level scenario as a genetic mechanism for the formation and accumulation aeolian deposits [56,57,58,59,60,61,62,63]. Differences in the timing of deposition and morphology of aeolian dunes during major sea level changes (high stands, low stands) among regions appear to be dependent on the interactions caused by minor-secondary sea level changes, local climatic conditions, sediment availability, and differences in the inland-near-shore topography [64]. However, our knowledge on the chronological framework of coastal deposits on the coasts of the Aegean is limited. Early studies have found sedimentary deposits on raised marine terraces along the coasts of Crete, which reveals the existence of characteristic Pleistocene marine fauna such as Cantharus viverratus, (Kien), Conus testudinarius (Mart), Polinices lacteus (Guild), Patella safiana, (Lmk), and Strobus bubonius (Lmk) [65,66,67,68,69].

1.4. Coastal Erosion

Coastlines in Greece are generally in a state of constant erosion due to anthropogenic interventions and climate change. Studies on the area [70,71] have estimated that about 25% of the Aegean coast consists mainly of low-lying coasts and deltas in a state of erosion. The northern coastline of Crete experiences extensive erosion, the highest severity in Greece, with erosion affecting 65.8% of its coastal zones. This large percentage of erosion is mainly due to the presence of extensive coastal zones with generally low beach slope, which are thus more vulnerable to erosion and the high frequency of large waves created by north and northeast winds in the Aegean [72].

1.5. Archaeological Setting

The coastal area of the Gazi-Ammoudara region is 8 km northwest of the palatial center of Knossos and not far from the Psiloritis mountain region. The region belonged to the administrative Minoan center of Tylissos, which is among the most poorly archaeologically investigated areas considering its proximity to some of the most important Minoan sites of Crete. This gap in archaeological research provides bias in regional and inter-regional archaeological interpretations and affects our understanding of the political, social, and economic networks in north-central Crete. Consequently, for heritage management, this gap is even more pronounced, since the lack of data and limited knowledge about existing archaeological sites and features implies difficulties in implementing effective strategies for their preservation. The archaeological importance of the region has been documented by early scholars [73,74]. Several archaeological remains in the area were destroyed by natural processes and modern human activity (agriculture, grazing, and construction). Evidence from rescue excavations, surface remains, and chance finds dating from the Neolithic to the Ottoman period have been found in many parts of the area, suggesting intensive patterns of habitation and other human activity. However, the only excavation in the coastal zone was carried out in the late 1970s. The excavation revealed a building, which was only partially excavated (Figure 2). It seems that it was a domestic unit, which most likely belonged to an extensive Neopalatial (MM III-LM IA) settlement that was part of the wider social contact network in the northern coastal zone of Minoan Crete. Pottery from the excavation is dated predominately to MM III-LM IA, though it includes Pre- and Postpalatial sherds [75,76].

1.6. Optically Stimulated Luminescence (OSL) Dating

Luminescence measurements are most frequently carried out on the quartz–rich fraction of aeolian dune deposits [77,78]. Environments rich in coastal and wind sediments [77] such as loess formations and dunes [79] have been shown to be suitable materials for Optically Stimulated Luminescence dating as they have been adequately exposed to sunlight before their deposition. However, in other types of sedimentary environments where there is doubt about their adequate exposure to sunlight, an investigation of the “zeroing” of the original luminance signal is required. Relatively recent applications of luminescence techniques in a variety of sedimentary environments [80], such as alluvial and colluvial sediments, glacial and peripheral deposits [81], riverbed sediments [82], sediments associated with cracks and seismicity [83,84], and terrace materials [85], have proven to be quite reliable, as the obtained ages were consistent with the chronological data found in existing reports. Optically Stimulated Luminescence (OSL) dating [86] has become a promising technique for dating fluvial, deltaic, and coastal sediments [27,87,88,89,90]. Thus, the Luminosity method is preferred as the most appropriate for investigating the location of the coastline [89,91,92,93,94,95]. The ages of the sediments from different depths, as they are determined by the method of luminosity, allow us to calculate the rate of sediment deposition over time and, therefore, the advance or erosion of an area such as the area of Ammoudara.
This paper presents the evolution of a coastal environment based on geomorphological, archaeological, and radio-chronological data that enabled the reconstruction of the coastal area of Ammoudara beach in northern Crete during the Bronze Age.

2. Study Area

The coastal area of Ammoudara (Figure 2) is located on the northern coast of Crete Island, approximately 2 km to the west of the city of Heraklion, Crete, Greece. The area is an alluvial plain and is characterized by low-lying regions, some parts of which have elevations either close to or even below the mean sea level. The entire coastal area has been shaped by significant geomorphological and tectonic processes as well as human activity [96]. The total coastline length is approximately 8 km, with a W to E trend, being slightly curved and having a subaerial width of up to 60 m. The beach is backed by a narrow dune field. where the dunes do not exceed 3 m in height. Most of the dune fields have been destroyed by anthropogenic activities. Three small rivers, the Gazanos (catchment area 192 km2), Xiropotamos (c.a. 35 km2), and Giofiros (c.a. 279 km2), discharge along the shoreline. At the west end of the beach, the Almiros stream drains the homonymous karstic spring (Figure 2). The Ammoudara coastal zone is microtidal with an astronomical tidal range of less than 10 cm [97], while the tide gauge installed in the harbor of Heraklion has recorded sea level fluctuations of more than 1 m over the period 1984–2020. These fluctuations incorporate the effects of astronomical tides and meteorological forcing. The Ammoudara beach zone is exposed primarily to NW, N, and NE wind-induced waves, with the NW direction being most frequent on an annual basis (23.6%), followed by the NW (12.4%) and the NE (3.7%). Wave heights are usually smaller than 2 m (88%), with 36% of them being less than 0.5 m and only 2% being over 4 m on an annual basis [98]. Likewise, 77.3% of the offshore wave periods are less than 5 s, while periods larger than 11 s have an annual frequency of only 0.74%. The maximum wave conditions are induced by NW and N winds during January and March, with a mean annual occurrence frequency of 0.01% and 0.03%, respectively. The maximum wave characteristics of the NW wind-induced waves are 2.3 m/4.6 s, while those of the prevailing N waves are 6 m/11 s [99].
A shallow submerged reef characterizes the shallow marine Ammoudara beach. The reef is ~4 km-long and has a mean width of around 35 m, while its rather planar crest lies, on average, at 0.6 m below the sea surface. Water depths exceed 3 m at its seaward side, while they are less than 2 m at its landward side. Previous studies in the area [99] have concluded that the reef represents a submerged beach rock formation, initially formed concurrently with the beach rocks along the western part of the present beach, based on the morphological and sedimentological characteristics of the reef. The reef consists of two distinct sedimentary layers: The upper layer, with a thickness of ~0.3, which consists of coarse-grained beach material having beackrock morphological characteristics (e.g., runnels, potholes, and grooves), and the lower layer, with a thickness greater than 0.5 m, which consists of fine-grained sandy material with cross bedding. The coastal zone of Ammoudara Beach can be subdivided into two sectors based on the presence or absence of the reef: The western sector, where the reef is absent, and the eastern sector, which is characterized by the continuous presence of the reef at less than 40 m from the shoreline.

3. Materials and Methods

Geomorphological and sedimentological mapping of the area was conducted with the use of satellite images, drone imagery, and field surveys. For the geomorphological mapping, topographic diagrams from the Hellenic Army Geographical Service and Sentinel satellite images were used. The geomorphological surveys included elevation and slope measurements taken with a RTK GPS. Additionally, depth soundings were acquired using a single beam echo-sounder reaching 20 m depth within the water. Such data were synthesized to create current bathymetry and beach elevations. Furthermore, surficial sediment samples were collected along the beach and analyzed according to Folk’s procedure [100]. To determine the depositional environment of the coastal sediments, samples from four 10 m-deep research cores were utilized (Figure 2) at approximately 60 m from the present-day shoreline. From the resulting stratigraphic sections, the different sedimentary horizons were studied macroscopically, the cores were subsampled, and the organisms found within them were identified to determine the deposition environment. For OLS dating, lightproof plastic tubes were used to avoid light contamination.

3.1. Optically Stimulated Luminescence (OSL)

Dating and radioactivity measurements were implemented in 12 core samples using OSL and Equivalent Doses calculations (for locations and depth of the samples, see Figure 2; for stratigraphic characteristics, section). For the radioactivity estimations, the top layer (2–4 cm) of the core samples was used, and 4–8 cm of the inner part of the core was collected, dried at 50 °C, and passed through sieves sequentially decreasing mesh from 150–125 μm to 100–63 μm. Grain sizes from 100–150 μm were chemically treated for quartz purification. Parts of the samples used for radioactivity assessment were dried at 105 °C for 24 h and passed through a 500 μm mesh sieve. About 60–70 g of this sieved material was sealed in plastic containers and kept for radon equilibrium for 1 month. Afterward, the containers were measured for natural radioactivity. For chemical pretreatment, small portions of 100–150 μm grain size of the samples were placed in separate plastic tubes. All samples were treated successively with HCl (8% w/w) to remove carbonates, H2O2 (30% w/w organic residuals and HF (40% w/w), and feldspars. The outer parts of the quartz grains were etched off the samples. Several water washes were performed between chemical treatments.

3.2. OSL Measurements and Equivalent Dose (EDs) Calculations

Purified quartz grains (aliquots) were placed on small stainless disks for luminescence measurements. Typically, 18 disks in total were used for each sample. The SAR protocol was followed [78,101], and Paleodoses were measured after IR and blue-OSL signal recording. The IR protocol was used to check for feldspar impurities. Successively increasing doses were induced on disks to obtain “dose response curves,” and the same protocol was used for these doses as described for natural ones. To eliminate sensitivity differences of the aliquots, the same dose was induced, referred as the test dose, and the ratios of dose to test dose signals were used to obtain the “dose-response curves.” Recycling and recuperation tests, as well as dose recovery tests, were performed for quality assurance. All measurements were performed in the Archaeometry Center of the University of Ioannina. Approximately 100 g of dried material was ground to powder and sieved through a 500 μm mesh sieve at 105 °C until reaching constant weight. About 60–70 g of each sample was placed in a plastic container of standard geometry and hermetically closed with a silicone sealant. The samples could reach the equilibrium of radon daughters (for about 1 month) and were measured for natural radioactivity. The age of a sample is calculated in terms of the equivalent dose ratio or Dose Rate (Gy) [102]. The equivalent dose is the total accumulated radioactivity to which the precipitate is exposed during the time it is deposited and, in the laboratory, and is called the Paleodose (De) [102]. The unit of measurement of absorbed radioactivity in the International System of Units (S.I) is Gray (Gy) ≥ 1 Gy = 1 Joule/kg. The dose rate is the rate at which natural radioactivity is absorbed by the crystalline granules per year. The age was calculated by the ratio equivalent dose (De) per Dose Rate [102].

3.3. Coastline Evolution

Aerial photographs and satellite images were used to study the stability of the coastline and to measure its displacement over time. Based on these data, the coastlines were digitized with the utmost accuracy through ArcGIS 10.1 software, through which conclusions were drawn for the shoreline displacement over time. This was followed by processing with Digital Shoreline Analysis System DSAS v.5 [103] to quantify the changes. Three different types of remote sensing data—historical analogue panchromatic aerial photographs, orthophotos, and natural-color satellite images—were used for the identification of the beach areas alterations during the last 74 years. The acquisition of the aerial photographs took place during 1945, 1968, and 1998 from Hellenic Military Geographical Service (HMGS); the orthophotos were acquired during 2005 and 2010 from Ktimatologio S.A; and the satellite images during 2009, 2012, 2013, 2016, 2017, and 2019 (Table 1). The aerial photos were scanned with a high-resolution scanner (1200 dpi), and all image data wee georeferenced in ArcGIS 10.1, with the 2010 Ktimatologio S.A. aerial photo as reference dataset, using approximately 40 control points along the beach and 20 landmarks. For all spatial data, the WGS 84/UTM zone 35N Reference System was used. The beach area and shorelines of the georeferenced images were manually digitized with the maximum accuracy, determined by each image’s pixel size. The average accuracy of all images was about 0.5 m. Transects perpendicular to the shoreline were created through DSAS with 50 m intervals. The statistic parameter extracted by the DSAS and presented in this work is the Net Shoreline Movement (NSM), which is the total distance between the earliest and most recent shorelines for each transect, in meters. Table 1 shows the date, type, scale, and resolutions of each image used in the analysis.

3.4. Terrain Visualisation

The contemporary landscape was created by draping the 1968 georectified aerial photography onto the DEM using 1968 topographic diagrams from the Hellenic Military Geographical Service–HMGS with 1:5000 scale. The 1968 topographic diagrams and the 1968 aerial photograph were used, since they are the only combined data of the area before the first constructions in the coastal zone (Figure 3). After this, the correlation of sea level changes, beach erosion rates, sedimentological rates, and the geomorphological characteristics were considered for the possible reconstruction of the area, with the elevation ranging from −10 m to +10 m.

4. Results

4.1. Textural Characteristics

The subaerial part of Ammoudara Beach consists mostly of sand (S) and slightly gravely sand ((g)S). The inner nearshore zone (between the shoreline and the reef) is generally characterized as slightly gravely sand ((g)S), while seaward of the reef and in the nearshore zone of the western sector the seabed is predominantly covered by sands (S). At the vicinity of the mouth of the rivers Almiros, Gazanos, and Xiropotamos and toward the eastern end of Ammoudara beach (close to the channelized Giofiros River), the sediments become coarser, varying from gravely sand (gS) to sandy gravel (sG). Relatively coarser material (i.e., sandy gravel (sG) and gravel (G)) is also observed along the landward side of the reef. In the western sector of Ammoudara Beach, the sediments become progressively finer from W to E. The same trend exists in the central and eastern parts of the beach, with the exceptions of the small areas with coarser sediments (Figure 4)

4.2. Coastal Evolution

For the study of the stability of the coastline from 1945 until today, aerial photographs and satellite images were used. From the examination of the aerial photographs for the years 1945, 196,8 and 1998, a clear retreat of the coast can be observed due to the constant and permanent erosion induced by the waves. The retreat of the coastline was estimated to be 45 m on average. In the diagrams of all checkpoints, a displacement rate of ± 1.5 m was calculated, which was due to the errors of georeferencing and digitization of the coastline of aerial photographs and satellite images. Small shoreline displacements (< ± 5 m) were considered within the limits of the expected transient changes in the terrain of the coastal zone due to ephemerous changes in hydrodynamic conditions and displacements due to transitions of the beach from cross section to storm cross section and, therefore, were not rated as significant. The results of the comparisons are presented in Figure 5 and Table 2.
From the 1945 aerial photograph (Figure 5), it is seen that the dune field in 1945 was more extensive, in respect to current state, with the wider coastal area having more accommodation space and availability of sediments. This indicates that during past periods, the area also had an extensive wide terrestrial coastal zone and an abundance of fine-grained sediment. The 1945 aerial photograph shows the boundary of the stabilized field of dunes (by sparse vegetation). Furthermore, the boundary of the fully developed field of dunes with denser vegetation is clearly visible. The analysis of the 1945 aerial photograph shows that, at that time, the area of the transect 37 and Drill G3 was not affected by the waves. Instead, the affected area was located at the marsh of Almiros karstic spring delta, which is in line with the reports of early residents of the area.

4.3. Stratigraphy

To determine the stratigraphy of the area, four research cores were drilled. The geographical characteristics of these drills are presented in Table 3 while their locations are shown in Figure 5.
The stratigraphy of the first drill (G1) indicates that the soil layer extends from the surface down to 3 m depth. It consisted of red sand soils and modern anthropogenic materials (e.g., building material, plastics). The next layer (3–6 m) was about 3 m-thick and consisted of sand with a small presence of clay (Mz: In the next 1.5 m (6–7.5 m), there were sandy deposits (Mz: 164.84 μm) in the upper part until 7.5 m. At a depth of 7.5 m, a layer of sandy gravel (Mz: 270.20 μm) appeared which extended, at least, down to 10 m. At a depth of 3.5 m, 5.8–6 m, and 8–8.3 m, blind samples with codes G1CAD1, G1CAD2, and G1CAD3 were collected in a plastic tube for optical luminosity analysis.
In the second drill (G2), the soil horizon was found at the 0–2 m drilling depth, which consisted of red sands and anthropogenic materials. Then, at the depths of 2–4 m, the sediments were sandy (Mz: 156.62 μm). At a depth of 3.8 m and for 20 cm (3.8 –4 m), a layer of slightly coarse-grained material was found in the presence porous igneous volcanic rock that was identified as Pumice. After that, and at least until the depth of 10 m, the sediments were gravelly sandy deposits (Mz: 266.51 μm). Five blind samples were taken at depths of 4.1–4.4 m, 4.9–5.2 m, 6.2–6.4 m, 7.2–7.4 m, and 9.5–10 m, with sample codes G2CBD1, G2CBD2, G2CBD3, G2CBD4, and G2CBD5. In the same area in 2016, an initial research drill (G2′) was performed for the collection of soil samples. The drilling data are in line with the latest data. The soil horizon, which consisted of red sands and carried anthropogenic materials, was found at the drilling depth of 0–1.8 m. In the next about 20 cm (1.8–2.0 m), the material was sand with a small presence of clay (Mz: 109.73 μm). In the next 1.5 m (2.0–3.5 m), there wee sandy deposites (Mz: 132.87 μm). At 3.5 m depth and for 20 cm (3.5–3.7 m), a layer with slightly coarse material and the presence of pumice was found. From 3.7 m depth to 4.5 m, the material was gravelly sand (Mz: 261.42 μm), while at 4.5 m depth, the material was aquifer, and a small horizon with clay sand at appeared the depth of 4.5–4.7 m. From this depth until the end of the drilling (4.7–8 m), the material was slightly gravelly sand (Mz: 259.72 μm). At a depth of 6.5 m (4.2 m below sea level), a remnant of pottery was found, which, after a conversation with the archaeologist who has the license to study the area, we discovered a possibility that it is of Minoan origin, as it has similarities with those found on the adjacent plot and at the same altitude [104].
In the third drilling (G3), the stratigraphy was slightly different due to the significant presence of organic materials. The soil horizon, which consisted of grey sand soils and anthropogenic materials, was visible at a drilling depth 0–0.4 m. In the next 1.6 m (approximately 0.4–2.0 m), the material was sandy (Mz: 124.57 μm), with a small presence of roots. In the next 4.2 m (2.0–6.m), there were grey sandy (Mz: 119.36 μm) and gravelly sandy (Mz: 221.31 μm) deposits, similar to coastal land dunes, in the vicinity of a river system. From 6 m to the end of the drilling depth of 10 m, there were grey sands (Mz: 216.12 μm) with the presence of cobbles (0.5–1 cm). At depths of 3.3 m, 6–6.2 m, and 9.8–10 m, blind samples with codes G3CARD1, G3CARD2, and G3CARD3 were sampled in plastic tubes for optical luminosity analyses. G3CARD1 and G3CARD2 consisted of slightly gravelly sand, while G1CAD3 was grey sandy gravel. The stratigraphy of all drills is shown in Figure 6.
In all three locations, each core showed that the deposits in the study area had the typical characteristics of subaerial coastal land and aeolian deposits in the upper sections, and no evidence of marine deposits was found. However, the coarser material in the lower section indicated that other coastal processes may affected the area. The stratigraphy of Ammoudara beach indicates that the dune system in the area is very old, since sand deposits were found in all cores with ages reaching 11.5 ± 3.15 ka. Sedimentological finds between 1.8 m and 8 m depth consisted of sandy and gravely sandy deposits, which indicate typical coastal deposits. Few different layers were recognized, e.g., in G2 and G2′, a layer between 3.5 and 3.7 m was found with slightly coarse material and the presence of pumice.

4.3.1. OLS Results

The equivalent doses (De) for the samples were estimated by interpolation of the natural dose in the dose response curve produced following the above protocols. The average values, with the corresponding standard deviation, are given in Table 4.
Quality assurance parameters of the OSL measurements of the samples were also calculated. The recycle ratio was well below 10%, and the recovery (signal of 0 Gy-induced dose compared to the natural dose) was less than 7.0% in all east borehole samples, 5.5% in central drill samples, and 3% in western drill samples. Finally, the dose recovery ratio (the estimated dose from a new dose response curve to a known given dose, i.e., ~ 15 Gy) was also less than 10%. All these values ensure the quality of the estimated equivalent doses for the samples. Recycling and dose recovery values are presented in Table 5.

4.3.2. Radioactivity Results

Measurements of the activity of isotopes of 238U, 235U, the 232Th series, and 40K and 87RB in Bq/kg and the calculated dose rate in Gy/ka are presented in Table 6.
Dose rates were corrected for the moisture content of the raw samples. Using the equivalent doses given in Table 6 (averages ± 1sd), the calculated ages and the range (min–max) are given in Table 7.
Distributions of the calculated ages for the samples are given in Figure 7, while individual ED values are given in Table 8.
Considering the values of the total of the disks used for the OSL measurements, the averages and the associated 1sd are shown in the 3rd, 7th, 11th, and 15th rows of Table 8. The methodology used to calculate the mean age and relative standard deviation considers the total number of tablets (i.e., pellet samples, placed in suitable tablets for passage through the machine), which is N = 18. Then, the methodology applies appropriate exclusion criteria and rejects some noncontributing tablets, namely N = 6, and from the rest N = 12, and gives values that are rounded to multiples of 50 years. Considering the values of the total number of tablets used for the OSL measurements, the averages and the relative sd were obtained.
Excluding the mean calculations, the top three (G1CAD1, G2CBD1, G3CARD1) and bottom three ages (G1CAD3 G1CAD4 G2CBD5), as well as the calculated averages, were still in the 1sd range. Samples G1CAD1, G1CAD2, G1CAD3, and G1CAD4 presented relative errors of 33.9%, 15.3%, 12.3%, and 10.5%, respectively, which were lower than 53.0%, 42.7%, 27.7%, and 31.5% of the initial values (Table 9). For samples G2CBD1, G2CBD2, G2CBD3, G2CBD4, and G2CBD5, the calculated averages were still in the 1sd range and showed relative errors of 20.0%, 14.1%, 9.6%, 14.4%, and 14.0%, respectively, which were lower than 43.6%, 25.4%, 19.0%, 32.4%, and 25.9% of the initial values (Table 9). For samples G3CARD1, G3CARD2, and G3CARD3, the calculated averages were still in the 1sd range and showed relative errors of 10.1%, 11.5%, and 9.9%, which were lower than the 16.8%, 19.6%, and 16.1% of the initial values (Table 9). Finally, typical OSL curves of the samples are given in Figure 8.

5. Discussion

In the Eastern part, core samples consistently produced slightly younger ages, mainly around 6ka. In the central part, a clearer stratigraphy is present. Ages started from around 5 ka ago at 4 m depth and 6.5ka ago at 6.5 m depth, where archaeological remains dating back to the Bronze age were found. Note that the Minoan civilization was a Bronze Age Aegean civilization on the island of Crete, and other Aegean Islands were flourishing from c. 3000 BC to c. 1450 BC. At the depth of 7.2 m, the age of the sediments appeared to be around 9ka, while at the depth of 9.5 m, the age was around 10 ka. The older age cluster of sediment ages fell between 10.5 and 12.0 ka and at a depth of 7.5–8.5 m in the western area of the beach. This period encompassed the period of the initial growth of the dune field while it was affected by the Almiros karstic spring delta. Based on those findings, the sedimentological rate for the formation of the dune fields in each area was 0.84 m/1000 y in the east part, 1.04 m / 1000 y in the central part, and 1.4 m / 1000 y in the western part, which indicates that a faster sedimentological process that may have been related to the availability of river sediments (Figure 9). To examine the evolution of the coastal landscape of the area, this work aimed to determine the rate and extent of beach evolution and how it might affect the discoverability of archaeological sites. If the sediment transfer and coastal erosion were rapid enough, they may not have been easily discernible using OSL dating as the progression of ages would fall within the ages’ standard errors, especially at specific locations. However, it appears that any dune evolution in the area, represented by the wider range of ages for the samples from 4 m depth, as well as the younger age cluster from shallower depth, did not provide any indications of age inconsistencies. At all drill sites, OSL ages indicate a correlation between depth and ages. Notably, the 4.1 m depth sample in Drill G2 is consistent with the depth of the Minoan site excavation, and the pumice layer indicates a volcanic event [105].
Based on the above findings, a potential DEM of the area was generated by inputting the tectonic movements and the sea level rise rates for the last 6000 years and combining these with archaeological data from the same area. The proposed DEM can be used as the future basis for a more detailed reconstruction to better understand the landscape evolution of the area [106]. However, there are several methodological limitations that need to be addressed. The topographic diagram of 1960, although of relatively large scale (1:5000), requires careful interpretation of boundaries and features. The map is also a record of property and agricultural land. This raises the important question of whether the mapping should be constrained by the actual features shown on a particular map or whether other sources of documentary evidence and physical factors (soil type, topography) should be used to interpret the past and create a representative rather than an actual landscape. Based on the topographic and bathymetric field measurements, the current elevations and bathymetry were created. The Ammoudara beach zone of 6ka BP appears wider by 80 m in the reconstructed DEM, and the dune field is more developed (Figure 10). The area behind the dune system had a higher altitude, while the area where the Minoan building was found at a distance from the dune system. During the last decades, there has been a noticeable loss of beach sediments, as indicated by the degradation of the dunes and the overall decline of the coastline.
Erosion is also associated with the increased distance between the reef and the coastline in the central and eastern parts of the coastal zone and the presence of coastal rocks in its western part. The loss of large amounts of sediment is intricately linked to the circulation of water on land and the relative transport of sediment, as well as human intervention and the relative rise in sea level. The former includes reductions in land sediment flows because of infrastructure development, coastal urbanization, and arbitrary disturbance of the dunes, while the latter may be related to climate change and sea level rise of >3 mm/year since 1980 [107], as well as a possible natural reduction of the coastal zone through the condensation of river deposits.
This gradual loss of marine sediments in recent decades has caused changes in the topography of the sea (depth, reef development) to which the hydrodynamic conditions have been dynamically adapted, continuing the removal of marine sediments. This loss of sediment is generally associated with sediment transport to the east, while the loss of sediment in the offshore zone is favored by reef openings in front of river estuaries. These findings agree with the conceptual model of Alexandrakis et al. [100]. In this model, it is reported that as the sea level rose, newer dunes eroded, and stable wind reserves came to the surface. Then, as the sea level reached its current height, coastal erosion removed the sediment from the beach, and finally, the beachrock separated from the beach, sank, and formed the current reef. In this model, there were multiple lines of dunes, with the older being quite stable and hosting an alluvial area behind it. As the sea level rose, younger dunes were eroded, and stable aeolian deposits were brought to the surface. Concurrently, a beachrock formation was developing along the coastline. In the areas of Giofiros and Almiros, the low topographic relief and muddy sediments act as indicators of a potential lagoon system behind the dune field.
The proposed DEM and shoreline position of the Ammoudara beach 6 ka ago is in agreement with findings from archeological excavations in the central Northern Crete. The sea level rise of 1.20 m between the Minoan Protopalatial and Neopalatial phase could be associated with a large seismic event that most likely occurred around 1700 BC [105]. The change of 1.45 m between the sea level stand of the Neopalatial phase (−2.70 ± 0.15 m) and the sea level stand with an upper age dated to the fourth century BC (−1.25 ± 0.25 m) would be either contemporary with the second destruction of the Minoan centers (~1450 BC) or in one or more paroxysmal events between 1450 BC and 400 BC. After the Roman conquest of Crete and during the subsequent centuries, 35 earthquakes from moderate-to-large magnitude (MO 6.5 to 8.2) have been reported in this region [108,109,110,111,112,113]. Particularly, the earthquake of 365 (MO 8.2) split the island and caused crustal uplift in western Crete. The tsunami triggered by this event struck coasts of most parts of the Mediterranean basin [25,26,114]. Sea level indicators depict the submerged ancient coasts of Crete. The presence of submerged tidal notches and beachrocks observed in the nearby coasts provide evidence for the existence of old shores [115]. The proximity of installations such as breakwaters, jetties, ship sheds, and fish tanks provide sufficient sea level indicators. Measurements of the emerged part with respect to the average sea level, the estimated time of construction and abandonment, and the evaluation of both the height and functional depth to the mean sea level play important roles in understanding mean sea level changes in correlation to archaeological remains [116]. In Sicily, the sluice gate of a sunken fish tank helped archaeologists deduce that the functional height of the Roman quarry’s floor (nowadays sunken) would be estimated at around 0.60 m above sea level [117]. The reported rock-cut ship sheds and fish tanks from Crete are in the vicinity of coastal quarries. The fish tanks at Chersonissos, 35 km away from the study area, could be correlated with the beachrock quarry. The floors of three fish tanks at Chersonissos, arranged in a row, now lie between 1.95m and 2.20m below modern sea level [69,118]. To the north at Agia Pelagia, 10km north from the study area, hollow channels and some holes in the shape of tubes have been spotted, which sometimes protruded up to 0,15 m from the surface of the sandstone. Those were interpreted as remnants of previous dunes that were covered with plants [119]. The same applies to the coastal quarry of Malia 40 km from Ammoudara [120,121]. Some of these holes could have been opened by the Minoans quarrymen for the insertion of poles and wooden beams [121]. Moreover, 10 km east of Amnisos, archeological findings of coastal infrastructures indicate a −2m sea level at 3500 BP [122] (for locations, see Figure 11).
Similar approaches of proposing conceptual geomorphological evolution models were proposed for the coastal region of the Sylt Island, a barrier island in the southern North Sea of North China that presented the sediment material sources during glacial and interglacial periods [123]. Moreover, in a case study in the Itapicuru coastal plain, north coast of Bahia State, paleoenvironmental information and OSL dating were used to define the possible geomorphological evolution scenario during the Upper-Middle Holocene of marine terraces and floodplain, where it indicated a slow and gradual regressive phase that was responsible for the highest marine terrace’s formation [124]. Additionally, Naples coast since the Middle Holocene was reconstructed using geomorphological analysis based on historical sources, ancient pictures and maps, and high-resolution digital terrain model (DTM) combined with iconographic research. The resconstruction led researchers to propose three new paleoenvironmental scenarios of the study area during the Holocene, which were derived from the evaluation of the relative sea level changes and vertical ground movements of volcano-tectonic origin [125].
When applying OSL dating in archaeology, one challenge is how to properly interpret the ages and their contexts. OSL could be used to date deposits that are related to otherwise defined cultural layers [126,127,128] or deposits that are directly associated with cultural material, where it essentially dates human disturbance of that sediment. Consequently, to assess the reliability of an OSL age, a detailed understanding of the archaeological context is especially critical to be able to evaluate whether enough of that disturbance had occurred to fully reset the luminescence signal. The archaeological site in the area near drills G2 and G2′ is a well-dated site which contains numerous archaeological features and ceramic artifacts [86]. Moreover, in the research drill G2′, a small ceramic was found at the 6 m depth (~4 m from current sea level) [104]. The OSL samples, which came from drill G2 above and below the dated layer, were dated to 1.4 ± 0.1 ka at both 0.3 and 0.6 m depth. The OSL results date the disturbed sand deposit. It is not necessarily surprising that OSL ages are slightly older. The assumption that the samples were disturbed by human activity, combined with the potential lack of complete sunlight exposure of the sand, make it likely that some residual signal may have been retained in the grains, causing some overestimation of the OSL age [129]. Also, it is possible that there was no more extensive human presence there and that the ~2 ka OSL ages may indicate some natural sand remobilization that occurred at that time. Whereas localized dune formation cannot be discounted, it is not possible to fully interpret these dates without further archaeological testing, as well as further OSL dating at the same depth in other parts of the Ammoudara dune field. This illustrates how important a broader archaeological and sedimentary context is for the interpretation of OSL data. Often, sparse artifact scatters are incorporated into sand deposits without clear indication of whether they are in their original context or perhaps within a buried lag surface, making the interpretation of OSL dates and their relationship with the original occupation time extremely difficult. Even if that is the case, careful examination of OSL data can provide insight into the conceptual development of a coastal landscape [95,130,131,132].

6. Summary and Conclusions

The puzzling scarcity of archaeological sites on the Ammoudara beach prompted an investigation into the chronology of the development of this dune field. The collected soil samples were dated using Optically Stimulated Luminescence to refine the chronology of the area and to assess the variability of ages with depth. The use of drilling in the dune field enabled the collection of samples from a 10 m-long stratigraphic column and the capturing of the early stages of the dune formation (OSL ages ~9.0–9.6 ka), while the OSL ages from 6 m depth represented the timing of its stabilization (OSL ages ~5–6 ka). The dune field appears to have formed long before the Bronze Age and became stabilized with only localized episodes of dune reactivation occurring thereafter. Based on the sediment texture, these localized episodes occurred mainly near river mouths. The shore-parallel reef in the Ammoudara coastal zone is a submerged beachrock that was formed during the Upper Holocene [99,133]. Its current position, about 0.5 m below sea level, indicates a relative sea level rise of more than 0.5 m, whereas the latter has caused a coastline retreat of approximately 50 m. During the last decade, the erosion rates have decreased and the beach zone of Ammoudara seems to be more stable, reaching a morphodynamics equilibrium with the nearshore hydrodynamic conditions. Luminescence dating of the aeolian deposits in the coastal area of Ammoudara yielded ages ranging from the Bronze Age to before the Bronze Age, suggesting that the coastal dunes in the area accumulated during the past 10,000 years. The study also revealed that, at a depth of 6 m below the ground surface, sediments dated back as old as the Bronze Age. Whatever the extent of the more recent dune reactivation in the study area, it was probably localized and not substantial enough to cause significant sand accumulation that would place archaeological deposits outside of the reach of standard proportionate methods. Since the geologic setting does not appear to preclude the archaeologists’ access to potential resources, the sparse distribution pattern of human occupation on the Ammoudara Beach might be better explained by exploring cultural habitation choices of past populations. Overall, the landscape reconstruction suggests that, geomorphologically, the area was rather different in the Bronze Age. Based on the archaeological and geomorphological findings, alluvial plain sediments seem to have been appropriate for agriculture since they were higher than sea level, while the location of shoreline was indicated approximately 200 m away from its current position. This can be also supported by the agricultural activities that are visible in the 1945 aerial photograph. Future work will need to address the depth to which OSL ages are adversely impacted by turbation due to wave and river processes.

Author Contributions

Conceptualization, G.A., S.P and N.A.K.; methodology, G.A. and S.P.; software, G.A. and S.P.; validation, G.A., S.P and N.A.K.; formal analysis, G.A. and S.P.; investigation, G.A. and S.P.; resources, G.A. and S.P.; data curation, G.A. and S.P.; writing—original draft preparation, G.A.; writing—review and editing, G.A., S.P and N.A.K.; visualization, G.A. and S.P.; supervision, G.A. and N.A.K. project administration, N.A.K.; funding acquisition, N.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors would like to thank Antonios Parasiris, from Institute of Applied & Computational Mathematics, Foundation for Research & Technology-Hellas, for his effort in English language editing, and the anonymous reviewers which by their comments and suggestions improved our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Walsh, K. (Ed.) Sea-Level Change and Coastal Settlement—Human Engagements with Littoral Environments. In The Archaeology of Mediterranean Landscapes: Human-Environment Interaction from the Neolithic to the Roman Period; Cambridge University Press: Cambridge, UK, 2013; pp. 30–67. [Google Scholar]
  2. Meulenkamp, J.; Jonkers, A.; Spaak, P.P.; VI Colloquium. Late Miocene to early Pliocene development of Crete. In Proceedings of the VI Colloquium Geology of the Aegean Region, Athens, Greece; 1977; pp. 137–149. [Google Scholar]
  3. Meulenkamp, J.E.; van der Zwaan, G.J.; van Wamel, W.A. On late miocene to recent vertical motions in the Cretan segment of the Hellenic arc. Tectonophysics 1994, 234, 53–72. [Google Scholar] [CrossRef]
  4. Angelier, J.; Lybéris, N.; Le Pichon, X.; Barrier, E.; Huchon, P. The tectonic development of the hellenic arc and the sea of crete: A synthesis. Tectonophysics 1982, 86, 159–196. [Google Scholar] [CrossRef]
  5. Peters, J.M. Neogene and Quaternary Vertical Tectonics in the South Hellenic Arc and Their Effect on Concurrent Sedimentation Processes; Gua: Amsterdam, The Netherlands, 1985; p. 247. [Google Scholar]
  6. Mourtzas, N. Tectonic Movement of the Coasts of Eastern Crete during the Quaternary. Ph.D. Thesis, Technical University of Athens, Athens, Greece, 1990; 480p. (In Greek). [Google Scholar]
  7. Kelletat, D. Perspectives in coastal geomorphology of western Crete, Greece. Z. Geomorphol. Suppl. 1996, 102, 1–19. [Google Scholar]
  8. Peterek, A.; Schwarze, J. Architecture and Late Pliocene to recent evolution of outer-arc basins of the Hellenic subduction zone (south-central Crete, Greece). J. Geodyn. 2004, 38, 19–55. [Google Scholar] [CrossRef]
  9. Wegmann, K.W. Tectonic Geomorphology above Mediterranean Subduction Zones: Northeastern Apennines of Italy and Crete, Greece; Lehigh University: Bethlehem, PA, USA, 2008. [Google Scholar]
  10. Tiberti, M.M.; Basili, R.; Vannoli, P. Ups and downs in western Crete (Hellenic subduction zone). Sci. Rep. 2014, 4, 5677. [Google Scholar] [CrossRef] [Green Version]
  11. Boekschoten, G.J. Some geological observations on the coasts of Crete. Geol. Mijnb. 1963, 42, 241–247. [Google Scholar]
  12. Hafemann, D. Die Niveauveränderungen an den Küsten Kretas seit dem Altertum; Akad. d. Wiss. ud Literatur; Steiner in Komm.: Wiesbaden, Germany, 1966. [Google Scholar]
  13. Kelletat, D. Geomorphologische Studien an den Kuesten Kretas. Beitraege zur Regionalen Kuestenmorphologie des Mittelmeerraumes; Vandenhoeck & Ruprecht: Göttingen, German, 1979. [Google Scholar]
  14. Mourtzas, N.; Kolaiti, E.; Anzidei, M. Vertical land movements and sea level changes along the coast of Crete (Greece) since Late Holocene. Quat. Int. 2016, 401, 43–70. [Google Scholar] [CrossRef]
  15. Le Pichon, X.; Angelier, J. The Hellenic arc and trench system: A key to the neotectonic evolution of the eastern Mediterranean area. Tectonophysics 1979, 60, 1–42. [Google Scholar] [CrossRef]
  16. Ten Veen, J.H.; Meijer, P.T. Late Miocene to Recent tectonic evolution of Crete (Greece): Geological observations and model analysis. Tectonophysics 1998, 298, 191–208. [Google Scholar] [CrossRef]
  17. Fytrolakis, N. The Geological Structure of Crete. Problems, Observations and Conclusions. Habilitation Thesis, National Technical University of Athens, Athens, Greece, 1980; 143p. [Google Scholar]
  18. Caputo, R.; Catalano, S.; Monaco, C.; Romagnoli, G.; Tortorici, G.; Tortorici, L. Active faulting on the island of Crete (Greece). Geophys. J. Int. 2010, 183, 111–126. [Google Scholar] [CrossRef]
  19. Caputo, R.; Monaco, C.; Tortorici, L. Multiseismic cycle deformation rates from Holocene normal fault scarps on Crete (Greece). Terra Nova 2006, 18, 181–190. [Google Scholar] [CrossRef]
  20. Mouslopoulou, V.; Moraetis, D.; Benedetti, L.; Guillou, V.; Bellier, O.; Hristopulos, D. Normal faulting in the forearc of the Hellenic subduction margin: Paleoearthquake history and kinematics of the Spili Fault, Crete, Greece. J. Struct. Geol. 2014, 66, 298–308. [Google Scholar] [CrossRef]
  21. Flemming, N.C.; Pirazzoli, P.A. Archéologie des côtes de la Crète. Histoire Archéologie Dossiers 1981, 50, 66–81. [Google Scholar]
  22. McCoy, F.W.; Heiken, G. Tsunami generated by the Late Bronze age eruption of Thera (Santorini), Greece. Pure Appl. Geophys. 2000, 157, 1227–1256. [Google Scholar] [CrossRef]
  23. Sigurdsson, H.; Carey, S.; Alexandri, M.; Vougioukalakis, G.; Croff, K.; Roman, C.; Sakellariou, D.; Anagnostou, C.; Rousakis, G.; Ioakim, C. Marine investigations of Greece’s Santorini volcanic field. Trans. Am. Geophys. Union 2006, 87, 337–342. [Google Scholar] [CrossRef] [Green Version]
  24. Bruins, H.J.; MacGillivray, J.A.; Synolakis, C.E.; Benjamini, C.; Keller, J.; Kisch, H.J.; Klügel, A.; Van Der Plicht, J. Geoarchaeological tsunami deposits at Palaikastro (Crete) and the Late Minoan IA eruption of Santorini. J. Archaeol. Sci. 2008, 35, 191–212. [Google Scholar] [CrossRef] [Green Version]
  25. Flouri, E.T.; Kalligeris, N.; Alexandrakis, G.; Kampanis, N.A.; Synolakis, C.E. Application of a finite difference computational model to the simulation of earthquake generated tsunamis. Appl. Numer. Math. 2013, 67, 111–125. [Google Scholar] [CrossRef]
  26. Stiros, S.C. The 8.5+ magnitude, AD365 earthquake in Crete: Coastal uplift, topography changes, archaeological and historical signature. Quat. Int. 2010, 216, 54–63. [Google Scholar] [CrossRef]
  27. Ott, R.F.; Gallen, S.F.; Wegmann, K.W.; Biswas, R.H.; Herman, F.; Willett, S.D. Pleistocene terrace formation, Quaternary rock uplift rates and geodynamics of the Hellenic Subduction Zone revealed from dating of paleoshorelines on Crete, Greece. Earth Planet. Sci. Lett. 2019, 525, 115757. [Google Scholar] [CrossRef]
  28. Cooper, J.A.G. Beachrock formation in low latitudes: Implications for coastal evolutionary models. Mar. Geol. 1991, 98, 145–154. [Google Scholar] [CrossRef]
  29. Neumeier, U. Experimental modelling of beachrock cementation under microbial influence. Sediment. Geol. 1999, 126, 35–46. [Google Scholar] [CrossRef]
  30. Kneale, D.; Viles, H.A. Beach cement: Incipient CaCO3-cemented beachrock development in the upper intertidal zone, North Uist, Scotland. Sediment. Geol. 2000, 132, 165–170. [Google Scholar] [CrossRef]
  31. Moore, C.H. Intertidal carbonate cementation Grand Cayman, West Indies. J. Sediment. Petrol. 1973, 43, 591–602. [Google Scholar]
  32. Beier, J.A. Diagenesis of Quaternary Bahamian beachrock: Petrographic and isotopic evidence. J. Sediment. Petrol. 1985, 55, 755–761. [Google Scholar] [CrossRef]
  33. Chowdhury, S.Q.; Fazlul Haq, A.T.M.; Hasan, K. Beachrock in St. Martin’s Island, Bangladesh: Implication of Sea Level Changes on Beachrock cementation. Mar. Geod. 1997, 20, 89–104. [Google Scholar] [CrossRef]
  34. Plomaritis, T. Morphology and Geochemistry of the Beachrocks of Sifnos (Greece). Master’s Thesis, University of Southampton, Southampton, UK, 1999. [Google Scholar]
  35. Holail, H.; Rashed, M. Stable isotopic composition of carbonate-cemented recent beachrock along the Mediterranean and the Red Sea coasts of Egypt. Mar. Geol. 1992, 106, 141–148. [Google Scholar] [CrossRef]
  36. Tatumi, S.H.; Kowata, E.A.; Gozzi, G.; Kassab, L.R.P.; Suguio, K.; Barreto, A.M.F.; Bezerra, F.H.R. Optical dating results of beachrock, eolic dunes and sediments applied to sea-level changes study. J. Lumin. 2003, 102–103, 562–565. [Google Scholar] [CrossRef]
  37. Calvet, F.; Cabrera, M.C.; Carracedo, J.C.; Mangas, J.; Perez-Torrado, F.J.; Recio, C.; Trave, A. Beachrocks from the island of La Palma (Canary Islands, Spain). Mar. Geol. 2003, 197, 75–93. [Google Scholar] [CrossRef]
  38. Font, Y.; Calvet, F. ‘Beachrocks’ Holocenos de la Isla de La Reunion, Oceano Indico. Cuad. Geol. Iber. 1997, 22, 87–102. [Google Scholar]
  39. Poulos, S.E.; Ghionis, G.; Verykiou, E.; Roussakis, G.; Sakellariou, D.; Karditsa, A.; Alexandrakis, G.; Petrakis, S.; Sifnioti, D.; Panagiotopoulos, I.P.; et al. Hydrodynamic, neotectonic and climatic control of the evolution of a barrier beach in the microtidal environment of the NE Ionian Sea (eastern Mediterranean). Geo Mar. Lett. 2015, 35, 37–52. [Google Scholar] [CrossRef]
  40. Vousdoukas, M.I.; Velegrakis, A.F.; Plomaritis, T.A. Beachrock occurrence, characteristics, formation mechanisms and impacts. Earth Sci. Rev. 2007, 85, 23–46. [Google Scholar] [CrossRef]
  41. Hopley, D. Beachrock as a sea-level indicator. In Sea-Level Research: A Manual for the Collection and Evaluation of Data; van der Plassche, O., Ed.; Geo Books: Norwich, UK, 1986; pp. 157–173. [Google Scholar]
  42. Dermitzakis, M.; Theodoropoulos, D. Study of beach-rocks in the Aegean Sea. Observations on occurrences in SE Crete, Rhodes and Metopi. Ann. Geol. Pays Hell. 1975, 26, 275–305. [Google Scholar]
  43. Neumeier, U. Le Role de L’Activite Mocrobienne dans la Cementation Precoce des Beachrocks (Sediments Intertidaux). Ph.D. Thesis, University of Geneva, Geneva, Switzerland, 1998. (In French with English Abstract). [Google Scholar]
  44. Bezerra, F.H.R.; Barreto, A.M.F.; Suguio, K. Holocene sea-level history on the Rio Grande do Norte State coast, Brazil. Mar. Geol. 2003, 196, 73–89. [Google Scholar] [CrossRef]
  45. Falkenroth, M.; Schneider, B.; Hoffmann, G. Beachrock as sea-level indicator—A case study at the coastline of Oman (Indian Ocean). Quat. Sci. Rev. 2019, 206, 81–98. [Google Scholar] [CrossRef]
  46. Rovere, A.; Raymo, M.E.; Vacchi, M.; Lorscheid, T.; Stocchi, P.; Gómez-Pujol, L.; Harris, D.L.; Casella, E.; O’Leary, M.J.; Hearty, P.J. The analysis of Last Interglacial (MIS 5e) relative sea-level indicators: Reconstructing sea-level in a warmer world. Earth Sci. Rev. 2016, 159, 404–427. [Google Scholar] [CrossRef] [Green Version]
  47. Karkani, A.; Evelpidou, N.; Vacchi, M.; Morhange, C.; Tsukamoto, S.; Frechen, M.; Maroukian, H. Tracking shoreline evolution in central Cyclades (Greece) using beachrocks. Mar. Geol. 2017, 388, 25–37. [Google Scholar] [CrossRef]
  48. Desruelles, S.; Fouache, É.; Ciner, A.; Dalongeville, R.; Pavlopoulos, K.; Kosun, E.; Coquinot, Y.; Potdevin, J.-L. Beachrocks and sea level changes since Middle Holocene: Comparison between the insular group of Mykonos–Delos–Rhenia (Cyclades, Greece) and the southern coast of Turkey. Glob. Planet. Chang. 2009, 66, 19–33. [Google Scholar] [CrossRef]
  49. Fairbridge, R.W.; Johnson, D.L. Eolianite. In The Encyclopedia of Sedimentology; Fairbridge, R.W., Bourgeois, J., Eds.; Dowden Hutchinson Ross: Stroudsburg, PA, USA, 1978; pp. 279–282. [Google Scholar]
  50. Bateman, M.D.; Holmes, P.J.; Carr, A.S.; Horton, B.P.; Jaiswal, M.K. Aeolianite and barrier dune construction spanning the last two glacial–interglacial cycles from the southern Cape coast, South Africa. Quat. Sci. Rev. 2004, 23, 1681–1698. [Google Scholar] [CrossRef] [Green Version]
  51. Jacob, T.; Wahr, J.; Pfeffer, W.T.; Swenson, S. Recent contributions of glaciers and ice caps to sea level rise. Nature 2012, 482, 514–518. [Google Scholar] [CrossRef]
  52. Mauz, B.; Hijma, M.P.; Amorosi, A.; Porat, N.; Galili, E.; Bloemendal, J. Aeolian beach ridges and their significance for climate and sea level: Concept and insight from the Levant coast (East Mediterranean). Earth Sci. Rev. 2013, 121, 31–54. [Google Scholar] [CrossRef]
  53. Mauz, B.; Vacchi, M.; Green, A.; Hoffmann, G.; Cooper, A. Beachrock: A tool for reconstructing relative sea level in the far-field. Mar. Geol. 2015, 362, 1–16. [Google Scholar] [CrossRef]
  54. Hearty, P.J.; Pascal, K. The Stratigraphy and Surficial Geology of New Providence and Surrounding Islands, Bahamas. J. Coast. Res. 1997, 13, 798–812. [Google Scholar]
  55. Murray-Wallace, C.V.; Bourman, R.P.; Prescott, J.R.; Williams, F.; Price, D.M.; Belperio, A.P. Aminostratigraphy and thermoluminescence dating of coastal aeolianites and the later Quaternary history of a failed delta: The River Murray mouth region, South Australia. Quat. Geochronol. 2010, 5, 28–49. [Google Scholar] [CrossRef]
  56. Pye, K.; Tsoar, H. Aeolian Sand and Sand Dunes; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  57. Pye, K.; Tsoar, H. The Formation of Sand Seas and Dune Fields. In Aeolian Sand and Sand Dunes; Springer: Berlin/Heidelberg, Germany, 2009; pp. 141–173. [Google Scholar]
  58. Andreucci, S.; Clemmensen, L.B.; Murray, A.S.; Pascucci, V. Middle to late Pleistocene coastal deposits of Alghero, northwest Sardinia (Italy): Chronology and evolution. Quat. Int. 2010, 222, 3–16. [Google Scholar] [CrossRef]
  59. Frechen, M.; Neber, A.; Tsatskin, A.; Boenigk, W.; Ronen, A. Chronology of Pleistocene sedimentary cycles in the Carmel Coastal Plain of Israel. Quat. Int. 2004, 121, 41–52. [Google Scholar] [CrossRef]
  60. Rodríguez-Vidal, J.; Finlayson, G.; Finlayson, C.; Negro, J.J.; Cáceres, L.M.; Fa, D.A.; Carrión, J.S. Undrowning a lost world—The Marine Isotope Stage 3 landscape of Gibraltar. Geomorphology 2013, 203, 105–114. [Google Scholar] [CrossRef]
  61. Fornós, J.J.; Clemmensen, L.B.; Gómez-Pujol, L.; Murray, A.S. Late Pleistocene carbonate aeolianites on Mallorca, Western Mediterranean: A luminescence chronology. Quat. Sci. Rev. 2009, 28, 2697–2709. [Google Scholar] [CrossRef]
  62. Bardají, T.; Goy, J.L.; Zazo, C.; Hillaire-Marcel, C.; Dabrio, C.J.; Cabero, A.; Ghaleb, B.; Silva, P.G.; Lario, J. Sea level and climate changes during OIS 5e in the Western Mediterranean. Geomorphology 2009, 104, 22–37. [Google Scholar] [CrossRef] [Green Version]
  63. Del Valle Villalonga, L.; Timar-Gabor, A.; Fornós, J.J. Geomorphological Processes and Environmental Interpretation at Espalmador islet (Western Mediterranean). J. Mar. Sci. Eng. 2019, 7, 144. [Google Scholar] [CrossRef] [Green Version]
  64. Vacher, H.; Hearty, P.; Rowe, M.; Curran, H.; White, B. Stratigraphy of Bermuda: Nomenclature, concepts, and status of multiple systems of classification. Spec. Pap. Geol. Soc. Am. 1995, 300, 271–294. [Google Scholar]
  65. Boekschoten, G.J. Beachrock at Limani Chersonisos, Crete. Geol. Mijnb. 1962, 41, 3–7. [Google Scholar]
  66. Gaki-Papanastassiou, K.; Karymbalis, E.; Papanastassiou, D.; Maroukian, H. Quaternary marine terraces as indicators of neotectonic activity of the Ierapetra normal fault SE Crete (Greece). Geomorphology 2009, 104, 38–46. [Google Scholar] [CrossRef]
  67. Dermitzakis, M.; Michail, C.; Mpasiakos, G.; Tripolitsioutou, F. Contricution to the absolute dating of beachrock by the means of the Thermoluminescence technique. In Proceedings of the 4th National Symposium on Oceanography and Fisheries, Rhodes Island, Greece, 26–29 April 1993; pp. 259–267, (In Greek with English Abstract). [Google Scholar]
  68. Kelletat, D. Zonality of Modern Coastal Processes and Sea-Level Indicators. Palaeogeogr. Palaeoclimatol. Palaeoecol. 1988, 68, 219–230. [Google Scholar] [CrossRef]
  69. Mourtzas, N.D. Fish tanks of eastern Crete (Greece) as indicators of the Roman sea level. J. Archaeol. Sci. 2012, 39, 2392–2408. [Google Scholar] [CrossRef]
  70. Alexandrakis, G.; Karditsa, A.; Poulos, S.; Ghionis, G.; Kampanis, N. An assessment of the vulnerability to erosion of the coastal zone due to a potential rise of sea level: The case of the Hellenic Aegean coast. In Εncyclopedia of Life Support Systems (EOLSS). Developed under the Auspices of the UNESCO; Sydow, A., Ed.; Eolss Publishers: Oxford, UK, 2010; Volume 1. [Google Scholar]
  71. Alexandrakis, G.; Ghionis, G.; Poulos, S.E.; Kampanis, N.A. Greece. In Coastal Erosion and Protection in Europe; Pranzini, E., Williams, A., Eds.; Routledge: London, UK, 2013; p. 355. [Google Scholar]
  72. Alexandrakis, G.; Poulos, S.Ε. An holistic approach to beach erosion vulnerability assessment. Sci. Rep. 2014, 4, 6078. [Google Scholar] [CrossRef] [Green Version]
  73. Hatzidakis, I. Minoikoi tafoi stin Kriti (in Greek). Archaiologikon Delt. 1918, 4, 45–87. [Google Scholar]
  74. Evans, A.J. The Palace of Minos at Knossos; Macmillan: London, UK, 1928; Volume 2, 232p. [Google Scholar]
  75. Hadzi-Vallianos, D. Prefecture of Heraklion, Gazi—Ammoudara. Archaiologikon Delt. 1979, 34, 384–385. [Google Scholar]
  76. Athanasaki, K. A coastal Neopalatial settlement in Ammoudara Heraklion: A preliminary approach. In Archaeological Work in Crete, Proceedings of the Third Meeting; University of Crete: Rethymnon, Greece, 2015; pp. 271–280. [Google Scholar]
  77. Murray, A.S.; Clemmensen, L.B. Luminescence dating of Holocene aeolian sand movement, Thy, Denmark. Quat. Sci. Rev. 2001, 20, 751–754. [Google Scholar] [CrossRef]
  78. Banerjee, D.; Murray, A.S.; Botter-Jensen, L.; Lang, A. Equivalent dose estimation using a single aliquot of polymineral fine grains. Radiat. Meas. 2001, 33, 73–94. [Google Scholar] [CrossRef]
  79. Roberts, R.G.; Galbraith, R.; Yoshida, H.; Laslett, G.; Olley, J.M. Distinguishing dose populations in sediment mixtures: A test of single-grain optical dating procedures using mixtures of laboratory-dosed quartz. Radiat. Meas. 2000, 32, 459–465. [Google Scholar] [CrossRef]
  80. Kitis, G.; Polymeris, G.; Pagonis, V.; Tsirliganis, N. Anomalous fading of OSL signals originating from very deep traps in Durango apatite. Radiat. Meas. 2013, 49, 73–81. [Google Scholar] [CrossRef]
  81. Duller, G.A.T. Recent developments in luminescence dating of Quaternary sediments. Prog. Phys. Geogr. Earth Environ. 1996, 20, 127–145. [Google Scholar] [CrossRef]
  82. Jain, C.; Singhal, D.; Sharma, M. Adsorption of zinc on bed sediment of River Hindon: Adsorption models and kinetics. J. Hazard. Mater. 2004, 114, 231–239. [Google Scholar] [CrossRef] [PubMed]
  83. Porat, N.; Amit, R.; Enzel, Y.; Zilberman, E.; Avni, Y.; Ginat, H.; Gluck, D. Abandonment ages of alluvial landforms in the hyperarid Negev determined by luminescence dating. J. Arid Environ. 2010, 74, 861–869. [Google Scholar] [CrossRef]
  84. Bowman, D.; Korjenkov, A.; Porat, N. Late-Pleistocene seismites from Lake Issyk-Kul, the tien shan range, Kyrghyzstan. Sediment. Geol. 2004, 163, 211–228. [Google Scholar] [CrossRef]
  85. Choi, J.H.; Murray, A.S.; Cheong, C.S.; Hong, D.G.; Chang, H.W. The resolution of stratigraphic inconsistency in the luminescence ages of marine terrace sediments from Korea. Quat. Sci. Rev. 2003, 22, 1201–1206. [Google Scholar] [CrossRef]
  86. Huntley, D.J.; Godfrey-Smith, D.I.; Thewalt, M.L.W. Optical dating of sediments. Nature 1985, 313, 105–107. [Google Scholar] [CrossRef]
  87. Jacob, J.; Gama, C.; Salgado, R.; Liu, J.T.; Silva, A. Shadowing Effects on Beach Morphodynamics during Storm Events on Tróia-Sines Embayed Coast, Southwest Portugal. J. Coast. Res. 2009, 73–77. [Google Scholar] [CrossRef]
  88. Madsen, A.T.; Murray, A.S. Optically stimulated luminescence dating of young sediments: A review. Geomorphology 2009, 109, 3–16. [Google Scholar] [CrossRef]
  89. Tsakalos, E.; Athanassas, C.; Tsipas, P.; Triantaphyllou, M.; Geraga, M.; Papatheodorou, G.; Filippaki, E.; Christodoulakis, J.; Kazantzaki, M. Luminescence geochronology and paleoenvironmental implications of coastal deposits of southeast Cyprus. Archaeol. Anthropol. Sci. 2018, 10, 41–60. [Google Scholar] [CrossRef]
  90. Nian, X.M.; Zhang, W.G.; Qiu, F.Y.; Qin, J.T.; Wang, Z.H.; Sun, Q.L.; Chen, J.; Chen, Z.Y.; Liu, N.K. Luminescence characteristics of quartz from Holocene delta deposits of the Yangtze River and their provenance implications. Quat. Geochronol. 2019, 49, 131–137. [Google Scholar] [CrossRef]
  91. Roberts, R.; Bird, M.; Olley, J.; Galbraith, R.; Lawson, E.; Laslett, G.; Yoshida, H.; Jones, R.; Fullagar, R.; Jacobsen, G.; et al. Optical and radiocarbon dating at Jinmium rock shelter in northern Australia. Nature 1998, 393, 358–362. [Google Scholar] [CrossRef]
  92. Bristow, C.S.; Duller, G.A.T.; Lancaster, N. Age and dynamics of linear dunes in the Namib Desert. Geology 2007, 35, 555–558. [Google Scholar] [CrossRef]
  93. Costas, S.; Ferreira, O.; Plomaritis, T.A.; Leorri, E. Coastal barrier stratigraphy for Holocene high-resolution sea-level reconstruction. Sci. Rep. 2016, 6, 38726. [Google Scholar] [CrossRef] [Green Version]
  94. Tamura, T.; Oliver, T.S.N.; Cunningham, A.C.; Woodroffe, C.D. Recurrence of Extreme Coastal Erosion in SE Australia beyond Historical Timescales Inferred from Beach Ridge Morphostratigraphy. Geophys. Res. Lett. 2019, 46, 4705–4714. [Google Scholar] [CrossRef] [Green Version]
  95. Wang, F.; Zhang, W.G.; Nian, X.M.; Ge, C.; Zhao, X.Q.; Cheng, Q.Z.; Chen, J.; Hutchinson, S.M. Refining the late-Holocene coastline and delta development of the northern Yangtze River delta: Combining historical archives and OSL dating. Holocene 2019, 29, 1439–1449. [Google Scholar] [CrossRef]
  96. Fassoulas, C. The tectonic development of a Neogene basin at the leading edge of the active European margin: The Heraklion basin, Crete, Greece. J. Geodyn. 2001, 31, 49–70. [Google Scholar] [CrossRef]
  97. Tsimplis, M.N. Tidal oscillations in the Aegean and Ionian Seas. Estuar. Coast. Shelf Sci. 1994, 39, 201–208. [Google Scholar] [CrossRef]
  98. Soukissian, T.; Hatzinaki, M.; Korres, G.; Papadopoulos, A.; Kallos, G.; Anadranistakis, E. Wind and Wave Atlas of the Hellenic Seas; Hellenic Centre for Marine Research Publ.: Athens, Greece, 2007. [Google Scholar]
  99. Alexandrakis, G.; Ghionis, G.; Poulos, S. The Effect of Beach Rock Formation on the Morphological Evolution of a Beach. The Case Study of an Eastern Mediterranean Beach: Ammoudara, Greece. J. Coast. Res. 2013, 69 (Suppl. 1), 47–59. [Google Scholar] [CrossRef]
  100. Folk, R.L. Petrology of the Sedimentary Rocks; Hemphill Publishing Company: Austin, TX, USA, 1980; 182p. [Google Scholar]
  101. Murray, A.S.; Wintle, A.G. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiat. Meas. 2000, 32, 57–73. [Google Scholar] [CrossRef]
  102. Aitken, M.J. Thermoluminescence Dating; Academic Press: London, UK, 1985. [Google Scholar]
  103. Himmelstoss, E.A.; Henderson, R.E.; Kratzmann, M.G.; Farris, A.S. Digital Shoreline Analysis System (DSAS) Version 5.0 User Guide; Open-File Report 2018–1179; U.S. Geological Survey: Reston, VA, USA, 2018. [Google Scholar]
  104. Alexandrakis, G.; Athanasaki, K.; Kampanis, N. Towards a GIS reconstruction of the Ammoudara coastal landscape, integrating archaeological evidence and geomorphological data. In Proceedings of the 2nd Conference on Computer Applications and Quantitative Methods in Archaeology Greek Chapter, Athens, Greece, 22–23 December 2016; pp. 45–50. [Google Scholar]
  105. Pareschi, M.T.; Favalli, M.; Boschi, E. Impact of the Minoan tsunami of Santorini: Simulated scenarios in the eastern Mediterranean. Geophys. Res. Lett. 2006, 33, L18607. [Google Scholar] [CrossRef]
  106. Ashmore, W. Social archaeologies of landscape. In A Companion to Social Archaeology; Blackwell Publishing Ltd.: Hoboken, NJ, USA, 2004; pp. 255–271. [Google Scholar]
  107. Solomon, S. Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel On Climate Change. Cambridge University Press: Cambridge, UK, 2007; 987p. [Google Scholar]
  108. Papazachos, V. The Earthquakes of Greece; Editions Ziti: Thessaloniki, Greece, 1997. [Google Scholar]
  109. Di Vita, A. Archaeologists and earthquakes: The case of 365 A.D. Ann. Geophys. 1995, 38, 5–6. [Google Scholar]
  110. Di Vita, A. Earthquakes and civil life at Gortyn (Crete) in the period between Justinian and Constant II (6–7th century AD). Archaeoseismology 1996, 7, 45–50. [Google Scholar]
  111. Spyropoulos, P.J.A.D. Chronicle of the Earthquakes of Greece from Antiquity until Today; Ziti Publications: Thessaloniki, Greece, 1997. [Google Scholar]
  112. Ambraseys, N. Earthquakes in the Mediterranean and Middle East: A Multidisciplinary Study of Seismicity up to 1900; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  113. Papadopoulos, G. A Seismic History of Crete; Ocelotos Publications: Athens, Greece, 2011; Volume 65. [Google Scholar]
  114. Guidoboni, E.; Comastri, A. The large earthquake of 8 August 1303 in Crete: Seismic scenario and tsunami in the Mediterranean area. J. Seismol. 1997, 1, 55–72. [Google Scholar] [CrossRef]
  115. Pirazzoli, P.; Evelpidou, N. Tidal notches: A sea-level indicator of uncertain archival trustworthiness. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2013, 369, 377–384. [Google Scholar] [CrossRef]
  116. Auriemma, R.; Solinas, E. Archaeological remains as sea level change markers: A review. Quat. Int. 2009, 206, 134–146. [Google Scholar] [CrossRef]
  117. Scicchitano, G.; Lo Presti, V.; Spampinato, C.; Gasparo Morticelli, M.; Antonioli, F.; Auriemma, R.; Ferranti, L.; Monaco, C. Millstones as indicators of relative sea-level changes in northern Sicily and southern Calabria coastlines, Italy. Quat. Int. 2011, 232, 92–104. [Google Scholar] [CrossRef] [Green Version]
  118. Leatham, J.; Hood, S. Sub-Marine Exploration in Crete, 1955. Annu. Br. Sch. Athens 1959, 53, 263–280. [Google Scholar] [CrossRef]
  119. Epting, M. Geologie der östlichen Talea Ori, Kreta. Diploma Thesis, The University of Freiburg, Freiburg im Breisgau, Germany, 1969. [Google Scholar]
  120. Shaw, J.W. Minoan Architecture: Materials and Techniques; Bottega d’Erasmo Università di Catania, Centro di Archeologia Cretese: Catania, Italy, 2009. [Google Scholar]
  121. Guest-Papamanoli, A. Les carrières de grès dunaire à Malia. Bull. Corresp. Hell. 1989, 113, 113–122. [Google Scholar] [CrossRef]
  122. Mouslopoulou, V.; Begg, J.; Nicol, A.; Oncken, O.; Prior, C. Formation of Late Quaternary paleoshorelines in Crete, Eastern Mediterranean. Earth Planet. Sci. Lett. 2015, 431, 294–307. [Google Scholar] [CrossRef]
  123. Tian, S.; Sun, J.; Lü, L.; Cao, M.; Zhang, Z.; Lü, T. Optically stimulated luminescence dating of late Quaternary loess deposits in the coastal region of North China: Provenance and paleoclimatic implications. Quat. Sci. Rev. 2019, 218, 160–177. [Google Scholar] [CrossRef]
  124. Lima, K.C.; Perez Filho, A. Evolução Geomorfológica de Planícies Costeiras no Holoceno. Mercator (Fortaleza) 2020, 19, e19010. [Google Scholar] [CrossRef]
  125. Mattei, G.; Aucelli, P.P.C.; Caporizzo, C.; Rizzo, A.; Pappone, G. New Geomorphological and Historical Elements on Morpho-Evolutive Trends and Relative Sea-Level Changes of Naples Coast in the Last 6000 Years. Water 2020, 12, 2651. [Google Scholar] [CrossRef]
  126. Bueno, L.; Feathers, J.; Blasis, P. The formation process of a paleoindian open-air site in Central Brazil: Integrating lithic analysis, radiocarbon and luminescence dating. J. Archaeol. Sci. 2013, 40, 190–203. [Google Scholar] [CrossRef]
  127. Sommerville, A.A.; Hansom, J.D.; Housley, R.A.; Sanderson, D.C.W. Optically stimulated luminescence (OSL) dating of coastal aeolian sand accumulation in Sanday, Orkney Islands, Scotland. Holocene 2016, 17, 627–637. [Google Scholar] [CrossRef]
  128. Sommerville, A.A.; Hansom, J.; Sanderson, D.C.W.; Housley, R. Optically stimulated luminescence dating of large storm events in Northern Scotland. Quat. Sci. Rev. 2003, 22, 1085–1092. [Google Scholar] [CrossRef]
  129. Bateman, M.D.; Boulter, C.H.; Carr, A.; Frederick, C.; Peter, D.; Wilder, M. Preserving the palaeoenvironmental record in Drylands: Bioturbation and its significance for luminescence-derived chronologies. Sediment. Geol. 2007, 195, 5–19. [Google Scholar] [CrossRef] [Green Version]
  130. Schokker, J.; Cleveringa, P.; Murray, A. Palaeoenvironmental reconstruction and OSL dating of terrestrial Eemian deposits in the southeastern Netherlands. J. Quat. Sci. 2004, 19, 193–202. [Google Scholar] [CrossRef]
  131. Dörschner, N.; Reimann, T.; Wenske, D.; Lüthgens, C.; Tsukamoto, S.; Frechen, M.; Böse, M. Reconstruction of the Holocene coastal development at Fulong Beach in north-eastern Taiwan using optically stimulated luminescence (OSL) dating. Quat. Int. 2012, 263, 3–13. [Google Scholar] [CrossRef]
  132. Shtienberg, G.; Dix, J.K.; Roskin, J.; Waldmann, N.; Bookman, R.; Bialik, O.M.; Porat, N.; Taha, N.; Sivan, D. New perspectives on coastal landscape reconstruction during the Late Quaternary: A test case from central Israel. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2017, 468, 503–519. [Google Scholar] [CrossRef]
  133. Dalrymple, R.A.; Kirby, J.T.; Hwang, P.A. Wave diffraction due to areas of energy dissipation. J. Waterw. Port Coast. Ocean Eng. 1984, 110, 67–79. [Google Scholar] [CrossRef]
Figure 1. Tectonic uplift as a result of the AD 365 earthquake (modified from [27]).
Figure 1. Tectonic uplift as a result of the AD 365 earthquake (modified from [27]).
Water 13 01269 g001
Figure 2. Location of area and hydrological basins in the study area. Known Minoan sites and research drill positions (a). Hydrological basins, digital image, and elevation provided by a 0.5 m DEM from Ktimatologio S.A. (b). Excavated Minoan site and location of G2 and G2′ research drill positions (c).
Figure 2. Location of area and hydrological basins in the study area. Known Minoan sites and research drill positions (a). Hydrological basins, digital image, and elevation provided by a 0.5 m DEM from Ktimatologio S.A. (b). Excavated Minoan site and location of G2 and G2′ research drill positions (c).
Water 13 01269 g002
Figure 3. The 1968 aerial photograph and the elevations of the beach area based on the 1968 HMGS 1:5000 topographic diagram.
Figure 3. The 1968 aerial photograph and the elevations of the beach area based on the 1968 HMGS 1:5000 topographic diagram.
Water 13 01269 g003
Figure 4. Textural characteristics (sand (S), slightly gravely sand ((g)S), gravely sand (gS), sandy gravel (sG), gravel (G)) and locations of the shoreline beachrock and reef of the Ammoudara beach.
Figure 4. Textural characteristics (sand (S), slightly gravely sand ((g)S), gravely sand (gS), sandy gravel (sG), gravel (G)) and locations of the shoreline beachrock and reef of the Ammoudara beach.
Water 13 01269 g004
Figure 5. Coastal retreat in meters during the period 1945–2020. Location of the research cores and the selected transects of Table 2. Expansion of the dune filed for the years 1945 and 2020 (above image). Net shoreline movements in meters per transect based on the DSAS assessments.
Figure 5. Coastal retreat in meters during the period 1945–2020. Location of the research cores and the selected transects of Table 2. Expansion of the dune filed for the years 1945 and 2020 (above image). Net shoreline movements in meters per transect based on the DSAS assessments.
Water 13 01269 g005
Figure 6. Stratigraphy of the drilling cores depth of samples and correlations between the various layers of sediments.
Figure 6. Stratigraphy of the drilling cores depth of samples and correlations between the various layers of sediments.
Water 13 01269 g006
Figure 7. Distribution of the calculated ages (in ka) for the number of samples.
Figure 7. Distribution of the calculated ages (in ka) for the number of samples.
Water 13 01269 g007
Figure 8. Typical OSL dose response curve curves. The y-axis shows the ratios obtained for each of these artificial doses (Lx/Tx) and for the natural sample (Ln/Tn).
Figure 8. Typical OSL dose response curve curves. The y-axis shows the ratios obtained for each of these artificial doses (Lx/Tx) and for the natural sample (Ln/Tn).
Water 13 01269 g008
Figure 9. OSL ages with 1 σ errors vs. depth for samples.
Figure 9. OSL ages with 1 σ errors vs. depth for samples.
Water 13 01269 g009
Figure 10. Current elevations, current position of the reef and the coastline, and extent of the dune fields for 1945 and 2020 (a). Estimated elevations and estimation of the extent of the dune field 6 ka (b). Representative profiles of current and estimated elevation, and estimated shoreline position 6 ka (c).
Figure 10. Current elevations, current position of the reef and the coastline, and extent of the dune fields for 1945 and 2020 (a). Estimated elevations and estimation of the extent of the dune field 6 ka (b). Representative profiles of current and estimated elevation, and estimated shoreline position 6 ka (c).
Water 13 01269 g010
Figure 11. Locations of archeological sea level indicators in the proximity of Ammoudara beach.
Figure 11. Locations of archeological sea level indicators in the proximity of Ammoudara beach.
Water 13 01269 g011
Table 1. Date, type, scale, resolution, and source of the images used for coastal evolution assessment.
Table 1. Date, type, scale, resolution, and source of the images used for coastal evolution assessment.
DateTypeScalePixel SizeSource
17/09/1945Aerial photo1:42,0003 mHMGS
1968Aerial photo1:15,0001 mHMGS
23/07/1998Aerial photo1:10,0001 mHMGS
2005Orthophoto1:15,0000.6 mKtimatologio S.A.
14/06/2009Satellite image-0.5 mGoogle Earth (maxar technologies)
20/06/2010Aerial photo1:30,0001.5 mHMGS
2011Orthophoto1:15,0000.6 mKtimatologio S.A.
06/11/2012Satellite image-0.2 mGoogle Earth (maxar technologies)
29/04/2013Satellite image-0.5 mGoogle Earth
27/06/2016Satellite image-0.2 mGoogle Earth (CNES/Airbus)
28/06/2017Satellite image-0.3 mGoogle Earth
24/03/2019Satellite image-0.2 mGoogle Earth (CNES/Airbus)
Table 2. Changes in the position of the broader beach coastline in meters, on selected transects (for location see Figure 5), during the period 1945–2020.
Table 2. Changes in the position of the broader beach coastline in meters, on selected transects (for location see Figure 5), during the period 1945–2020.
Transects1945–19981998–20022002–20092009–20152015–20201945–2020
200.0−1.93.4−3.3−4−3.4
37−17.8−11.03.614.7−12.4−23.1
50−33.0−15.815.2−7.3−0.6−41.4
75−40.9−8.70.6−0.3−4.7−53.9
90−49.2−10.23.3−2.93.3−55.6
110−53.7−7.22.6−3.511.2−50.5
125−42.4−7.4−2.5−3.99.8−31.45
Table 3. Topographical characteristics of the research drills.
Table 3. Topographical characteristics of the research drills.
NoLatLonAltitude Distance from the CoastlineDepth
G135°20′15.62″ N25°6′7.79″ E2.45 m67.5 m10 m
G235°20′14.73″ N25°4′53.74″ E2.33 m55 m10 m
G2′35°20′14.40″ N25°04′53.57″ E2.37 m51 m8 m
G335°20′22.75″ N25°3′33.88″ E2,62 m76.6 m10 m
Table 4. Equivalent doses estimated with OSL dose response curve lines.
Table 4. Equivalent doses estimated with OSL dose response curve lines.
G1G2G3
Sample CodeG1CA D1G1CAD2G1CAD3G1CA D4G2CBD1G2CBD2G2CBD3G2CBD4G2CB D5G3CARD1G3CARD2G3CARD3
Equivalent Dose (Gy)Mean ± 1sd (N = 18)3.4 ± 1.89.2 ± 3.99.9 ± 2.710.4 ± 3.33.7 ± 1.66.5 ± 1.75.6 ± 1.17.0 ± 2.36.9 ± 1.85.2 ± 1.78.0 ± 3.96.3 ± 1.0
Table 5. Quality assurance parameters of the OSL measurements.
Table 5. Quality assurance parameters of the OSL measurements.
Sample CodeRecycling RatioRecuperation (%)Dose Recovery Ratio
G1G1CAD11.06 ± 0.093.89 ± 3.290.93 ± 0.17
G1CAD21.04 ± 0.102.18 ± 1.580.96 ± 0.13
G1CAD31.02 ± 0.256.97 ± 6.020.91 ± 0.21
G1CAD40.98 ± 0.132.15 ± 2.440.98 ± 0.19
G2G2CBD11.08 ± 0.095.24 ± 5.880.92 ± 0.15
G2CBD21.11 ± 0.111.54 ± 3.350.94 ± 0.24
G2CBD31.08 ± 0.103.08 ± 2.720.97 ± 0.17
G2CBD41.01 ± 0.140.67 ± 2.970.91 ± 0.13
G2CBD51.01 ± 0.142.88 ± 5.410.95 ± 0.19
G3G3CARD11.06 ± 0.052.76 ± 2.270.95 ± 0.14
G3CARD21.02 ± 0.102.12 ± 2.000.93 ± 0.11
G3CARD31.04 ± 0.052.92 ± 1.680.98 ± 0.19
Table 6. Radioactivity measurements and calculated dose rates.
Table 6. Radioactivity measurements and calculated dose rates.
Isotope238U235U232TH40K87RBTotal
G1G1CAD1Radioactivity (Bq/kg)17.5 ± 11.52 ± 0.5911.9 ± 0.9185.9 ± 7.9
Dose rate (Gy/ka)0.337 ± 0.0220.012 ± 0.0050.175 ± 0.0130.494 ± 0.0210.0191.037 ± 0.033
G1CAD2Radioactivity (Bq/kg)18.4 ± 1.10.87 ± 0.6315.7 ± 1232.1 ± 8.5
Dose rate (Gy/ka)0.341 ± 0.0220.007 ± 0.0050.25 ± 0.0160.668 ± 0.0250.0191.285 ± 0.037
G1CAD3Radioactivity (Bq/kg)14 ± 0.80.91 ± 0.449.3 ± 0.7141.8 ± 5.9
Dose rate (Gy/ka)0.274 ± 0.0180.008 ± 0.0040.148 ± 0.0110.41 ± 0.0170.0190.860 ± 0.027
G1CAD4Radioactivity (Bq/kg)14.6 ± 0.90.66 ± 0.339.4 ± 0.7133.8 ± 6.9
Dose rate (Gy/ka)0.283 ± 0.0190.006 ± 0.0030.156 ± 0.0120.402 ± 0.0210.0190.866 ± 0.031
G2G2CBD1Radioactivity (Bq/kg)12.8 ± 0.70.97 ± 0.366.8 ± 0.5104.9 ± 4.5
Dose rate (Gy/ka)0.244 ± 0.0150.008 ± 0.0030.106 ± 0.0080.297 ± 0.0130.0190.675 ± 0.021
G2CBD2Radioactivity (Bq/kg)13.8 ± 0.81.17 ± 0.517.8 ± 0.6117.9 ± 6.6
Dose rate (Gy/ka)0.272 ± 0.0180.01 ± 0.0040.127 ± 0.010.344 ± 0.0190.0190.771 ± 0.028
G2CBD3Radioactivity (Bq/kg)9.9 ± 0.70.37 ± 0.298.9 ± 0.6168.1 ± 11.2
Dose rate (Gy/ka)0.176 ± 0.0130.003 ± 0.0020.14 ± 0.010.477 ± 0.0320.0190.815 ± 0.036
G2CBD4Radioactivity (Bq/kg)11.3 ± 0.70.36 ± 0.38 ± 0.6146.9 ± 6.6
Dose rate (Gy/ka)0.212 ±0.0150.003 ± 0.0020.125 ± 0.010.413 ± 0.0190.0190.772 ± 0.026
G2CBD5Radioactivity (Bq/kg)11 ± 0.70.41 ± 0.298.1 ± 0.6148.6 ± 16.1
Dose rate (Gy/ka)0.189 ± 0.0140.003 ± 0.0020.117 ± 0.0090.389 ± 0.0420.0190.717 ± 0.045
G3G3CARD1Radioactivity (Bq/kg)16.6 ± 0.91.07 ± 0.568 ± 0.6136.8 ± 6.7
Dose rate (Gy/ka)0.243 ± 0.0160.008 ± 0.0040.114 ± 0.0090.355 ± 0.0170.0190.739 ± 0.026
G3CARD2Radioactivity (Bq/kg)16.7 ± 0.90.84 ± 0.4211.3 ± 0.8190.7 ± 7.2
Dose rate (Gy/ka)0.301 ± 0.0190.007 ± 0.0040.181 ± 0.0120.55 ± 0.0210.0191.058 ± 0.031
G3CARD3Radioactivity (Bq/kg)13.2 ± 0.80.68 ± 0.469.3 ± 0.7155.7 ± 6.4
Dose rate (Gy/ka)0.255 ± 0.0170.006 ± 0.0040.151 ± 0.0110.457 ± 0.0190.0190.887 ± 0.028
Table 7. Equivalent doses, dose rates, and calculated ages in y before present.
Table 7. Equivalent doses, dose rates, and calculated ages in y before present.
Sample CodeEquivalent Dose (GY)Dose Rates (GY/KA)Ages (Y)Min–Max (Y)Moisture (%)
G1G1CAD13.4 ± 1.81.037 ± 0.0333300 ± 17501500–740021.280
G1CAD29.2 ± 3.91.285 ± 0.0377150 ± 30503600–17,55013.120
G1CAD39.9 ± 2.70.860 ± 0.02711,500 ± 31507800–21,00012.523
G1CAD410.4 ± 3.30.866 ± 0.03112,050 ± 38007900–25,3008.893
G2G2CBD13.7 ± 1.60.675 ± 0.0215500 ± 24002400–12,90014.591
G2CBD25.6 ± 1.10.815 ± 0.0366850 ± 13004450–950014.519
G2CBD36.5 ± 1.70.771 ± 0.0288450 ± 21505300–13,05011.732
G2CBD47.0 ± 2.30.772 ± 0.0269100 ± 29506150–16,45015.477
G2CBD56.9 ± 1.80.717 ± 0.0459650 ± 25005250–13,85022.991
G3G3CARD15.2 ± 1.70.739 ± 0.0266550 ± 11004850–880023.923
G3CARD28.0 ± 3.91.058 ± 0.0316900 ± 13503800–940012.839
G3CARD36.3 ± 1.00.887 ± 0.0287150 ± 11505600–950011.221
Table 8. Individual OSL ages (y).
Table 8. Individual OSL ages (y).
OSLG1G2G3
Age (Y)N CAD1CAD2CAD3CAD4CBD1CBD2CBD3CBD4CBD5CARD1CARD2CARD3
1151136047781789223775313444761705246485438085586
2155148628010866635055671507067025632541652515638
3179648868506927138426728581670097083545756816080
4207250148775934540006873593970547547567561536153
5216756719312998041146878625870648322575761766205
621925738991310,64443407102635571238659577663096217
72195583010,21911,08544957107636276188762598263616333
82211588810,34411,10145277461643378788960609864066630
92373631610,96911,25045627901644379489083623765067023
102709664711,31911,33646378065644682829265669468557273
113069694311,44011,367493384936816896910,100671070027456
123397723311,46311,979498886127071900910,264676571237487
133689726012,29311,998574988317082911311,037689180807775
144144785412,71812,655667010,0167286928411,380719482477825
155191837613,54313,949718010,232841011,09212,131778983197965
165202883214,36314,390804811,706847811,09512,743779283928266
176410978614,87715,033820311,899887816,06613,644828184489137
18742517,55820,97725,28212,90413,062949216,42713,835879994019516
Table 9. Calculated age averages with the associated standard deviation considering the total number of disks (N = 18) and excluding some of the values (N = 12). Values are rounded to multiples of 50 y.
Table 9. Calculated age averages with the associated standard deviation considering the total number of disks (N = 18) and excluding some of the values (N = 12). Values are rounded to multiples of 50 y.
Sample Code 1812
G1 G1CAD1Average ± 1sd (y)3300 ± 17502950 ± 1000
Relative error (%)53.033.9
G1CAD2Average ± 1sd (y)7150 ± 30506550 ± 1000
Relative error (%)42.715.3
G1CAD3Average ± 1sd (y)11,500 ± 315011,050 ± 1400
Relative error (%)27.412.7
G1CAD4Average ± 1sd (y)12,050 ± 380011,400 ± 1200
Relative error (%)31.510.5
G2 G2CBD1Average ± 1sd (y)5500 ± 24005000 ± 1000
Relative error (%)43.620.0
G2CBD2Average ± 1sd (y)6850 ± 13006750 ± 650
Relative error (%)19.09.6
G2CBD3Average ± 1sd (y)8450 ± 21508150 ± 1150
Relative error (%)25.414.1
G2CBD4Average ± 1sd (y)9100 ± 29508350 ± 1200
Relative error (%)32.414.4
G2CBD5Average ± 1sd (y)9650 ± 25009650 ± 1350
Relative error (%)25.914.0
G3G3CARD1Average ± 1sd (y)6550 ± 11006450 ± 650
Relative error (%)16.810.1
G3CARD2Average ± 1sd (y)6900 ± 13506950 ± 800
Relative error (%)19.611.5
G3CARD3Average ± 1sd (y)7150 ± 11507050 ± 700
Relative error (%)16.19.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alexandrakis, G.; Petrakis, S.; Kampanis, N.A. Integrating Geomorphological Data, Geochronology and Archaeological Evidence for Coastal Landscape Reconstruction, the Case of Ammoudara Beach, Crete. Water 2021, 13, 1269. https://doi.org/10.3390/w13091269

AMA Style

Alexandrakis G, Petrakis S, Kampanis NA. Integrating Geomorphological Data, Geochronology and Archaeological Evidence for Coastal Landscape Reconstruction, the Case of Ammoudara Beach, Crete. Water. 2021; 13(9):1269. https://doi.org/10.3390/w13091269

Chicago/Turabian Style

Alexandrakis, George, Stelios Petrakis, and Nikolaos A. Kampanis. 2021. "Integrating Geomorphological Data, Geochronology and Archaeological Evidence for Coastal Landscape Reconstruction, the Case of Ammoudara Beach, Crete" Water 13, no. 9: 1269. https://doi.org/10.3390/w13091269

APA Style

Alexandrakis, G., Petrakis, S., & Kampanis, N. A. (2021). Integrating Geomorphological Data, Geochronology and Archaeological Evidence for Coastal Landscape Reconstruction, the Case of Ammoudara Beach, Crete. Water, 13(9), 1269. https://doi.org/10.3390/w13091269

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop