Late Quaternary Marine Terraces and Tectonic Uplift Rates of the Broader Neapolis Area (SE Peloponnese, Greece)

<jats:p>Marine terraces are geomorphic markers largely used to estimate past sea-level positions and surface deformation rates in studies focused on climate and tectonic processes worldwide. This paper aims to investigate the role of tectonic processes in the late Quaternary evolution of the coastal landscape of the broader Neapolis area by assessing long-term vertical deformation rates. To document and estimate coastal uplift, marine terraces are used in conjunction with Optically Stimulated Luminescence (OSL) dating and correlation to late Quaternary eustatic sea-level variations. The study area is located in SE Peloponnese in a tectonically active region. Geodynamic processes in the area are related to the active subduction of the African lithosphere beneath the Eurasian plate. A series of 10 well preserved uplifted marine terraces with inner edges ranging in elevation from 8 ± 2 m to 192 ± 2 m above m.s.l. have been documented, indicating a significant coastal uplift of the study area. Marine terraces have been identified and mapped using topographic maps (at a scale of 1:5000), aerial photographs, and a 2 m resolution Digital Elevation Model (DEM), supported by extensive field observations. OSL dating of selected samples from two of the terraces allowed us to correlate them with late Pleistocene Marine Isotope Stage (MIS) sea-level highstands and to estimate the long-term uplift rate. Based on the findings of the above approach, a long-term uplift rate of 0.36 ± 0.11 mm a−1 over the last 401 ± 10 ka has been suggested for the study area. The spatially uniform uplift of the broader Neapolis area is driven by the active subduction of the African lithosphere beneath the Eurasian plate since the study area is situated very close (~90 km) to the active margin of the Hellenic subduction zone.</jats:p>

section_content: Coastal regions above subduction zones are of great interest for studying the active tectonic processes at convergent plate boundaries, especially if the presence of suitable geomorphic markers allows quantifying the deformation of the upper plate [1] [2] [3] [4]. Among such geomorphic markers are palaeoshoreline indicators (beach rocks, marine notches, marine terraces, and palaeo-strandlines) which indicate the former positions of shorelines and are commonly used to determine rates and patterns of long-term rock uplift and deformation in tectonically active settings [1, [5] [6] [7] [8] [9] [10] [11] [12]. 

Marine terraces are palaeo-platforms of marine origin that formed at the coast as relatively flat and generally smooth surfaces that dip gently seaward [13]. They are most commonly observed in regions that are tectonically uplifting. These landforms are prominent along many coastlines, and they can be useful for measuring rates and patterns of tectonic uplift. A marine terrace typically consists of two distinct surfaces: an abrasion ramp dipping gently seaward and an inland-bounding sea-cliff dipping seaward. The junction between the roughly flat paleo-platform and its internal palaeo-cliff is called the palaeoshoreline angle (also called the inner edge of the platform), and it is here that a wave-cut notch may be preserved [13]. The inner edge represents a marker that closely approximates the elevation of the local sea level at the time of terrace formation, typically during sea-level highstands [1, [14] [15] [16]. Hence this feature is an important palaeohorizontal indicator that permits the assessment of spatial variations in vertical motions of the crust [13]. 

The occurrence of a flight of marine terraces (a sequence of marine terraces with step-like morphology) in a coastal region is the result of the interaction between longterm tectonic uplift and Quaternary cyclic sea-level fluctuations, associated with global climatic changes [1, 6, [17] [18] [19] [20] [21]. Quaternary cyclic sea-level changes are typically represented by eustatic (global) curves derived from the Oxygen isotope ratios [22] [23] [24] [25]. These curves show a cyclic trend characterized by peaks corresponding to distinct marine interglacial highstands, represented by the odd-numbered Marine Isotope Stages (MIS), and marine glacial lowstands, indicated by the even-numbered MIS [26]. Paleoplatforms are thought to form during periods of relative sea-level stability (e. g., when the rates of rock uplift and sealevel change balance each other) and are abandoned following sea-level highstands [1, 5, 6, 8]. 

The study of marine terraces started in the 1970s in several places of the world and has since expanded to almost all the uplifting coastal regions [1]. The hypothesis that the mean uplift rate during the late Pleistocene is high enough to preserve the marine terraces permitted the dating of the terrace sequences and the evaluation of the uplift history of several raised coastal regions worldwide [7, 18, [27] [28] [29] [30]. The progress achieved is owed, to a great extent, to both the establishment of the theory of the Pleistocene eustatic sea-level fluctuations and the development of accurate dating techniques. It is obvious that to accurately evaluate vertical movements of coastal regions affected by Quaternary tectonic uplift, appropriate age constraints for marine terraces are required. 

The age of a marine terrace can be estimated through various dating methods. For very recent marine terraces, 14 C dating has been used to determine the age of emplacement of certain markers. Among the absolute dating methods of marine terraces, U-Th disequilibrium provides the least ambiguous one, but the presence of suitable fossils is required. Other, less precise dating methods also rely on scarce fossil material, e. g., amino acid racemization, d 18 O stratigraphy, and faunal assemblage correlations [20]. When organic remains are generally absent, Optically Stimulated Luminescence (OSL) dating methods can be applied to determine the ages of the marine terrace deposits [8, 21, 28, 31, 32]. Another dating method that has great potential for constraining the age of Pleistocene marine bedrock terraces is surface exposure dating, which is based on the analysis of in situ-produced cosmogenic nuclides such as 10 Be, 21 Ne, 26 Al, or 36 Cl in near-surface rocks [11, 33, 34]. Cosmogenic radionuclides (CRNs) produced in situ in the terrace cover deposits provide a relatively new tool for dating marine terraces [35] because 10 Be and 26 Al accumulation in quartz acts as a clock [36]. On flights of terraces where only one terrace can be dated, steady uplift is often assumed, and the altitudinal spacing of the terraces is used to assign sea-level highstand ages [5] [6] [7] 37]. 

In Greece, marine terraces are observed in various places such as the southern coasts of the Gulf of Corinth (the Corinth Rift) and several coastal regions along the emergent fore-arc of the Hellenic Subduction Zone, including the southwestern Peloponnese and the islands of Cephalonia, Kythera, Crete, Karpathos, and Rhodes [6, 10, 29, 30, [38] [39] [40] [41] [42]. The flight of marine terraces along the southern coast of the Gulf of Corinth, an intracontinental rift in Central Greece, has attracted considerable scientific interest since it is one of the fastest extending regions accessible on land worldwide [6, 11, 29, [42] [43] [44] [45] [46] [47]. Another wellstudied area is the island of Crete, which lies along the emergent fore-arc of the Hellenic Subduction Zone. Cretan paleo-shorelines have been extensively used to determine the spatial and temporal pattern of Quaternary uplift and explain the mechanisms of fore-arc deformation [4, [8] [9] [10] 20, [48] [49] [50]. 

This study aims to provide information about the long-term landscape evolution of the broader area of Neapolis, in southeastern Peloponnese, Greece, in response to tectonic uplift during the late Quaternary. For this purpose, a series of ten uplifted marine terraces, which represent former shorelines, have been identified and mapped in detail through Geographic Information System (GIS) techniques and extensive fieldwork, and selected samples from sediments of two marine terraces have been collected and dated using the OSL technique. Marine terraces have long been recognized in the study area by previous authors [51] [52] [53] [54] but have received little attention so far, and their dating and distribution, both in space and time, has never been successfully attempted before. The study of this outstandingly wellpreserved flight of marine terraces provides important information on the Quaternary uplift history, as well as on the process governing orogenesis and landscape evolution, which remains poorly understood for the study area. Age determination of two of the marine terraces allowed us to correlate their inner edges to past sea-level highstands [25, 26, 55], and in this way to estimate the long-term tectonic uplift rate for the study area over the last 401 ± 10 ka.
section_name: Introduction

section_content: The study area is located along the southeastern coast of the Lakonic Gulf in southeast Peloponnese and lies between cape Koulenti to the northwest, and Agia Marina, to the southeast (Figure 1 ). It also includes Elafonissos, a small island with an area of 19 km 2 separated from the mainland by a 1 km wide and 2. 5 m deep canal of seawater [56]. It is a tectonically active area since it is the easternmost terrestrial part of the southwestern Aegean fore-arc, situated approximately 90 km east-northeast of the Hellenic Subduction Zone where the African tectonic plate is subducting beneath the Eurasian plate at a mean rate of ~36 mm a -1 [57, 58]. Rates of crustal motion offshore the Peloponnese and Crete Island with respect to Africa, estimated from GPS measurements, are up to ~40 mm a -1 [59, 60]. Earthquakes, seismic tomography, and seismic reflection data reveal that the Hellenic Subduction Zone is the largest, fastest, and most seismically active subduction zone in the Mediterranean [61, 62]. Seismic tomography has imaged the subducted slab to depths > 600 km [61, 63], whereas Mesozoic to Cenozoic arc volcanism [64] from the eastern Mediterranean and the Aegean show that the Hellenic Subduction Zone has continued to consume the African plate since at least 100 Ma [65, 66]. Hence, the geodynamic processes in the region are related to the active convergent plate boundary, and the study area provides an example of actively deforming continental crust situated in the fore-arc of a subduction zone (Figure 1a ).
section_name: Study Area

section_content: The geotectonic setting of Greece (modified from Reilinger et al. [58], and Vott et al. [67] ) depicting the location of the study area and (b) hill shaded map of the study area produced from topographic maps at 1:5000 scale of the Hellenic Military Geographical Service. 

The Mesozoic basement of the study area consists of geological formations of the geotectonic units of Gavrovo-Tripolis and Olonos-Pindos and a unit of phyllites-quartzites (Figure 2 ). The Olonos-Pindos unit crops out in a very small part of the study area and The Mesozoic basement of the study area consists of geological formations of the geotectonic units of Gavrovo-Tripolis and Olonos-Pindos and a unit of phyllites-quartzites (Figure 2 ). The Olonos-Pindos unit crops out in a very small part of the study area and consists of tectonized radiolarites and limestones of the Triassic to Eocene age [68, 69]. The geotectonic unit of Tripolis is represented by Upper Triassic to Eocene dolomites, limestones, and dolomitic limestones, whereas its lower parts consist of very low-grade metamorphic formations (mainly clastic and carbonate rocks) of Triassic age (Tyros beds) (Figure 2 ). A relatively extensive part of the broader Neapolis area consists of a series of highly metamorphosed rocks consisting of phyllites, quartzites, as well as calcitic and dolomitic marbles of Permian-Triassic age [68, 69]. The Alpine basement is overlain by extensive outcrops of Pliocene and Pleistocene deposits (Figure 2 ). The occurrence of Pliocene-Pleistocene sedimentary sequences of marine origin up to ~180 m above presentday sea level is the most evident consequence of the widespread Quaternary uplift of the study area due to its proximity to the Hellenic Subduction Zone. Upper Pliocene to lower Pleistocene formations are composed of marine deposits, which are mainly pelites, sandstones, conglomerates, calcarenites, and carbonate rocks with red algae [54, 68, 69]. On top of the phyllites-quartzites unit, a series of Pleistocene formations of fluvio-terrestrial origin (clays, sands, loams, and angular rock fragments) occurs [51, 68, 69]. The Holocene deposits consist of talus cones, scree, and unconsolidated alluvial deposits, eluvial mantle materials, whereas coastal sand dunes are observed along the north, northwest, and south shoreline of Elafonissos Island and along the coast east of Cape Punta [56]. 2 ). A relatively extensive part of the broader Neapolis area consists of a series of highly metamorphosed rocks consisting of phyllites, quartzites, as well as calcitic and dolomitic marbles of Permian-Triassic age [68, 69]. The Alpine basement is overlain by extensive outcrops of Pliocene and Pleistocene deposits (Figure 2 ). The occurrence of Pliocene-Pleistocene sedimentary sequences of marine origin up to ~180 m above presentday sea level is the most evident consequence of the widespread Quaternary uplift of the study area due to its proximity to the Hellenic Subduction Zone. Upper Pliocene to lower Pleistocene formations are composed of marine deposits, which are mainly pelites, sandstones, conglomerates, calcarenites, and carbonate rocks with red algae [54, 68, 69]. On top of the phyllites-quartzites unit, a series of Pleistocene formations of fluvio-terrestrial origin (clays, sands, loams, and angular rock fragments) occurs [51, 68, 69]. The Holocene deposits consist of talus cones, scree, and unconsolidated alluvial deposits, eluvial mantle materials, whereas coastal sand dunes are observed along the north, northwest, and south shoreline of Elafonissos Island and along the coast east of Cape Punta [56]. 

Figure 2. Geological map of the study area modified from Theodoropoulos [51], IGME [68], and Lekkas et al. [69]. Fault kinematics are undifferentiated, and most faults are inactive. 

Previous researches led to the identification of uplifted marine terraces carved into the pre-existing Pliocene-Pleistocene formations and limestone bedrock at several locations in the study area. Theodoropoulos [51] reported 3 relatively flat surfaces at an elevation between 10 and 140 m above mean sea level (m. s. l. ). According to him, the 2 lower (at 10-15 m and 25-30 m, respectively) correspond to marine terraces of the Tyrrhenian age. [51], IGME [68], and Lekkas et al. [69]. Fault kinematics are undifferentiated, and most faults are inactive. 

Previous researches led to the identification of uplifted marine terraces carved into the pre-existing Pliocene-Pleistocene formations and limestone bedrock at several locations in the study area. Theodoropoulos [51] reported 3 relatively flat surfaces at an elevation between 10 and 140 m above mean sea level (m. s. l. ). According to him, the 2 lower (at 10-15 m and 25-30 m, respectively) correspond to marine terraces of the Tyrrhenian age. 

A sequence of 5 uplifted marine terraces and remnants of uplifted marine abrasion platforms, ranging in elevation from 4 m to approximately 140 m, mainly carved into the Pliocene-Pleistocene formations, have been recognized by Kelletat et al. [52] at the broader area of Neapolis, Profitis Ilias, and Agia Marina. According to the same authors, there is a sudden southward inclination and submergence of the older marine surfaces, whereas the altitudes of the terraces on the Elafonissos Island are higher than on the opposite mainland. 

Lekkas et al. [53] noted the presence of 3 relatively continuous marine terraces from 6 to 55 m in elevation, as well as several non-continuous higher marine terrace surfaces, up to 120 m, on the Elafonissos Island. Some of the uplifted marine terraces of the island cut into Pleistocene sedimentary formations, but there are also surfaces carved into the limestone bedrock. 

The general trend of the faults in the study area is mainly northwest-southeast, with some secondary ones having a northeast-southwest direction (Figure 2 ). Most of these faults are impossible to distinguish on the hill-shaded map, indicating that they are inactive.
section_name: Figure 1. (a)

section_name: Materials and Methods

section_content: To correlate the marine terraces of the study area in space and time, the first step was the detailed mapping of the distinct marine terrace surfaces and their inner edges (the palaeoshoreline angle of each palaeo-platform). Mapping was carried out using both Geographic Information System (GIS) and field techniques [40, 70, 71]. A GIS spatial geodatabase including various thematic layers (contour lines of 4 m, and 2 m or even 1 m in gently sloping areas, elevation points, and drainage network) was organized and implemented using Esri ArcGIS software. The main sources of these data were detailed analog topographic maps at 1:5000 scale, obtained from the Hellenic Military Geographical Service. The projection system was the Hellenic Geodetic Reference System 1987 (HGRS '87). GIS procedures produced secondary layers including a raster 2 m cell size DEM of the study area (to generate the DEM, the TOPOGRID algorithm in ArcGIS was used), a hill-shade map, a slope-aspect map, a range map, and a map of morphological discontinuities of slopes. Additionally, land surface parameters such as slope, gradient, aspect, curvature, and roughness were calculated. These derivatives, along with the above thematic maps and visual inspection and interpretation of georeferenced aerial photo mosaic (made from aerial photos taken by the Hellenic Military Geographical Service in 1996) and recent Google Earth Pro images, were particularly useful in identifying, delineating, and mapping the uplifted marine terraces of the study area. In this way, preliminary morphological sketches depicting marine terraces' surfaces and their inner edges (when possible) were produced in the form of printed paper maps at the scale of 1:5000, suitable for the field survey. 

The field survey focused upon the spatial arrangement of the uplifted abrasion platforms and the accuracy of the inner edges' elevations at specific locations. Field control was necessary to correct and validate marine terraces' maps/sketches drawn from the DEM analysis and aerial photos and Google Earth Pro images interpretation. The entire study area was investigated, verifying most of the uplifted marine terraces identified during the previous lab mapping phase. Where possible, a handheld Topcon GMS-2 Differential Global Positioning System (DGPS) was used to mark waypoints for the precise location of the inner edges that represent former shorelines. The accuracy of the DGPS used is considered to be suitable for 1:5000 scale field surveys. During the post-field survey mapping phase, the existing GIS database was integrated with the DGPS data, and a comparison of field and remote sensing mapping data took place to validate the initial map of the marine terraces sequence and proceed to the production of the final marine terraces map. The final map depicts a series of ten well-defined marine terraces grouped according to the elevation ranges of their inner edges. 

To verify the elevation of the different terraces and their inner edges, three topographic profiles, one for each locality of the study area, were drawn across the terraces, perpendicular to the coastline. The profiles are based on the 2 m DEM, while the inner shoreline angle elevations of the different marine terraces along the lines of the topographic profiles were measured in the field using the DGPS equipment.
section_name: Mapping

section_content: To correlate the marine terraces in space and time, two sedimentary block samples, designated as M and N, were collected from the caprock outcrops of two different terraces for OSL dating. The samples were collected at the locations of the topographic profiles 1 and 3, as close as possible to the marine terraces' inner edges. Sediments deposited in shallow marine conditions have a high chance of complete bleaching before burial [72, 73]. The OSL samples were referenced by handheld DGPS, sealed, and then transported to the Nuclear Physics Laboratory, Department of Physics, University of Ioannina, Ioannina, Greece, for luminescence dating. 

The outer ~3 cm from each of the 10-15 cm-sized sedimentary block samples was discarded, under dim red light, and then divided into 2 parts for the equivalent dose (ED) and annual dose (AD) analysis. The ED sub-samples were subjected to wet-sieving to obtain the 100-150 µm fractions, from which purified quartz grains were obtained by standard laboratory procedures [74]. The purity of the extracted quartz grains was checked using the OSL-IR depletion ratio [75]. The chemically purified quartz grains were mounted on 10 mm diameter stainless steel discs by covering a part with 3-4 mm diameter at the center of the disc [76] using acetone suspension. Following sample preparation, the OSL signal from the extracted quartz grains was measured using a Risø TL/OSL-DA-20 reader equipped with blue (470 nm) LEDs and a 7 mm Hoya U-340 optical detection filter [77]. The singlealiquot regenerative-dose (SAR) protocol of Murray and Wintle [78] was used to determine all ED values. All OSL measurements were made using stimulation with blue diodes at 125 • C for 40 s. The OSL signal was calculated using the initial 0. 8 s minus a background from the mean of the last 8 s of the stimulation curve at 90% LED power. For all OSL measurements, a preheat temperature of 240 • C (10 s) and a cut-heat temperature of 160 • C (0 s) were selected after performing a dose recovery-preheat plateau test on a representative sample using different preheat temperatures. Measurements were rejected if recycling ratios exceeded 10% from unity and recuperation exceeded 5% of the natural signal. 

The AD rate for each sample was calculated from radionuclide concentrations, measured by high-resolution gamma spectrometry, using the conversion factors of Liritzis et al. [79], and calculated using the DRAC software (v1. 1; [80] ).
section_name: Optically Stimulated Luminescence (OSL) Dating

section_content: It is generally accepted that flights of marine terraces result from the interplay of the rate of eustatic sea-level change and the rate of vertical tectonic motion [13, [81] [82] [83] [84]. Hence the sequence of uplifted marine terraces observed in the landscape of the study area is assumed to be the geomorphic record of late Quaternary periods of relative sealevel stability such as glacio-eustatic sea-level highstands [1, 5], which happened during odd-numbered MIS [1]. 

Late Pleistocene uplift rates were determined by elevation-age relationships using the following equation [1, 8, 28] :

where U is the uplift rate, E is the present-day elevation of the palaeo-shoreline for a particular terrace, A is the age of the respective terrace, and e corresponds to the elevation of the sea level at the time of terrace formation. 

Since many sea-level curves that show a significant variability of sea-level through time have been derived for different parts of the world [26, 55, 85], there are some approaches where apparent uplift rates have been calculated by neglecting the eustatic component (e = 0, e. g., [3, 86, 87] ). In this way, errors attributed to past eustatic sea-level positions can be discarded by obtaining a relative uplift rate with respect to present-day sea level. In our study, we did not follow this approach since it may introduce a substantial bias [16]. To model the deformation of several different palaeo-horizontal lines, the use of a common reference level is necessary. To refer palaeo-shorelines to the present-day sea level requires corrections. In principle, prominent shorelines that formed during past interglacial periods need corrections because they correspond to sea-level highstands that, in the global record, are likely to be within a few or several meters of the present interglacial. In this paper, we calculate uplift rates by using the eustatic sea-level curve proposed by Waelbroeck et al. [25], which is built from a multitude of different δ18O records from all around the globe. We chose this curve because it is a good approximation of global eustatic sea-level change since it is a composite relative sea-level curve, over the last four climatic cycles, based on long benthic isotopic records retrieved at one North Atlantic and one Equatorial Pacific site. Moreover, it provides a confidence interval, which in turn is useful to calculate uplift rate errors. 

Although Equation ( 1 ) is simple, it is difficult to estimate precise long-term uplift rates because of the errors on the different terms involved in it, especially on the eustatic sea-level curve (e) as well as on the accuracy of the elevation measurements of the present-day inner edge of the terrace (E). Therefore, the calculation of a precise uplift rate value is complicated, and we focused on obtaining uplift rate ranges, including the different uncertainties using a modified version of Equation ( 1 ) [28] (Figure 3 ):

where the delta symbols (∆) represent the estimated variability in the different parameters. We applied equations 1, 2, and 3 for the two OSL-dated marine terraces of the study area and a long-term uplift rate has been estimated. The age of the other undated terraces has been inferred from their correlation with past MIS sea-level highstands. Assuming a temporally constant uplift rate [4-6,8,20,89] and taking into account the duration of each MIS [88] as well as the maximum and minimum position of eustatic sea-level for each MIS highstand [25], we estimated the "expected" (or "predicted") elevation ranges of the inner edge for each MIS marine terrace. The predicted elevations were then compared with the observed elevations of the marine Errors regarding the eustatic sea-level curve are mainly associated with uncertainties in the age of the marine terraces (A) since the palaeoshoreline does not represent a specific age but the duration of the sea-level highstands, as well as with the position of the sea level at that time (e) (Figure 3 ). To estimate the uncertainties due to the age extend of the MIS, we adopted the time boundaries of the MIS highstands from Lisiecki and Raymo [88]. 

The errors of estimation of the marine terrace inner edges' elevation are associated with uncertainties on the accuracy of both the 2 m cell size DEM and the DGPS measurements. Another error in elevation measurements is related to the precision of the definition of the inner edges as the shoreline angle of emerged older marine terraces is commonly concealed when the paleo-cliffs and/or paleo-platforms are affected by scarp diffusion, river incision, or aeolian depositional processes, or a combination thereof [16]. For long-term uplift rate estimates, we took into consideration a measurement error of ± 1 m for the elevation of the inner edge associated with the accuracy of the DEM and the DGPS system used, while an additional error range of ±1 m due to mislocation of the inner edges' elevation was also considered [16, 41]. Therefore, the accuracy of the determination of the palaeo-shoreline angle elevations is quite good as the total average vertical uncertainty in elevation estimates of shoreline angles due to both the vertical mislocation error of the inner edge and the vertical standard mean square error on the DEM and DGPS used, is ±2 m. It should be noted that this error is negligible in estimating the long-term Quaternary uplift rates involving periods of tens to hundreds of thousands of years [88]. We applied Equations ( 1 ), ( 2 ) and ( 3 ) for the two OSL-dated marine terraces of the study area and a long-term uplift rate has been estimated. 

The age of the other undated terraces has been inferred from their correlation with past MIS sea-level highstands. Assuming a temporally constant uplift rate [4] [5] [6] 8, 20, 89] and taking into account the duration of each MIS [88] as well as the maximum and minimum position of eustatic sea-level for each MIS highstand [25], we estimated the "expected" (or "predicted") elevation ranges of the inner edge for each MIS marine terrace. The predicted elevations were then compared with the observed elevations of the marine terraces, and for each terrace, an MIS was assigned. This is a model generally used since the dating of marine terraces is sometimes difficult, and very often, only one-if any-terrace can be reliably dated.
section_name: Rates of Coastal Uplift Estimation and Eustatic Correlations

section_name: Results and Discussion

section_content: The application of mapping techniques led to the identification of a sequence of 10 uplifted marine terraces along the coastline between cape Koulendi and Agia Marina, west of cape Maleas, with inner edges' elevation ranging between 8 ± 2 m and 172 ± 2 m (Figures 4 and 5 ). On Elafonissos Island, the marine terraces flight, consisting of 10 emerged abrasion platforms, can be traced with shoreline angles (inner edges) ranging in elevation from 8 ± 2 m to 192 ± 2 m. A detailed map of the marine terraces (numbered T1-T10 starting from the youngest and lowermost towards the highest and oldest) is presented in Figure 4. 

Marine terraces in the study area are mainly characterized by well preserved wave-cut surfaces carved on marine and lacustrine sediments of the Pleistocene age. Some of them, especially those along the coast between Paleokastro and Sfontili and partially those on the Elafonissos Island, have developed on Mesozoic limestone and dolomites in the form of bedrock terraces, whereas between Neapolis and Paleokastro, they are built on top of the metamorphic rocks of the Mesozoic bedrock. At the area of Neapolis, remnants of dissected marine platforms occur on fluvio-terrestrial deposits of lower Pleistocene age. Between Cape Punda and Neapolis, the lateral morphological continuity of the terraces is interrupted due to the fluvial erosion of the marine-lacustrine Late Pleistocene sedimentary formations caused by the channels of well-developed drainage networks that dissect the palaeo-surfaces. 

T10 starting from the youngest and lowermost towards the highest and oldest) is presented in Figure 4. The lithology of the geological formations and their relative erodibility seem to play an important role in the extension of the marine terraces' surface. Marine terraces developed on Pliocene-Pleistocene formations as well as on metamorphic rocks are more extensive in relation to those carved into limestone and dolomites. The surfaces of some of the bedrock marine terraces eroded into Mesozoic limestones and dolomites have also been affected by karst processes. 

Three topographic profiles were drawn across the terraces and roughly perpendicular to the coastline at three locations of the study area (Figure 6 ). None of the three profiles contains the entire terrace sequence. The elevation values measured for the shoreline angles of the different terraces along the profiles are illustrated in Figure 6. Topographic Profile A was constructed 2. 3 km northwest of Cape Punta (Figure 5a, b ). It is 2. 2 km long and extends from sea level to a maximum elevation of 165 m. It contains seven marine terraces (terrace T7, along with the highest terraces T9 and T10, are not represented in this profile. A sedimentary block sample (N) for OSL dating was collected from the caprock of T5, very close to the inner edge of this terrace (Figures 5e and 6 ). Profile 2, constructed at the western part of the Elafonissos Island, has a length of 740 m and reaches a maximum elevation of 160 m (Figure 6 ). It revealed a staircase of seven marine terraces (T1, T2, T3, T4, T6, T7, and T8) separated by relatively steep palaeo-cliffs. Samples taken from the sedimentary blocks of the two lower terraces along this profile did not give reliable OSL ages. Topographic Profile 3 was constructed 4. 6 km south of Neapolis. It is 1. 21 km long and extends from the coastline to a maximum elevation of 140 m. This profile cut a sequence of six terraces (T1, T2, T3, T4, T6, and T7). A sedimentary block sample (M) was taken from the caprock close to the inner edge of T3 for OSL dating. The elevation of the inner edges of the terraces remains roughly constant in all three profiles (Figure 6 ), indicating that there is no evidence for coast-parallel tilting of the marine terraces in the study area. Three topographic profiles were drawn across the terraces and roughly perpendicular to the coastline at three locations of the study area (Figure 6 ). None of the three profiles contains the entire terrace sequence. The elevation values measured for the shoreline angles of the different terraces along the profiles are illustrated in Figure 6. Topographic Profile A was constructed 2. 3 km northwest of Cape Punta (Figure 5a, b ). It is 2. 2 km long and extends from sea level to a maximum elevation of 165 m. It contains seven marine terraces (terrace of the two lower terraces along this profile did not give reliable OSL ages. Topographic Profile 3 was constructed 4. 6 km south of Neapolis. It is 1. 21 km long and extends from the coastline to a maximum elevation of 140 m. This profile cut a sequence of six terraces (T1, T2, T3, T4, T6, and T7). A sedimentary block sample (M) was taken from the caprock close to the inner edge of T3 for OSL dating. The elevation of the inner edges of the terraces remains roughly constant in all three profiles (Figure 6 ), indicating that there is no evidence for coast-parallel tilting of the marine terraces in the study area. 4. The marine terraces are numbered consecutively, starting with the lowest one. The elevations of the inner edges for the terraces, as derived from the topographic profiles, are provided. The location and the age of the OSL samples, as well as the MIS attributed to each marine terrace, are also given.
section_name: Mapping

section_content: The results of the application of the OSL dating technique to the samples collected from the caprock outcrops of two of the terraces are depicted in Table 1. 

To estimate mean uplift rates, we modeled the elevations of the inner edges of the terraces T3 and T5 with the OSL dates of the sedimentary block samples collected from the caprock of these palaeo-platforms. Our preferred correlation of these two OSL dated marine terraces with the MIS highstands of the sea-level curve proposed by Waelbroeck et al. [25] allowed us to quantify rate ranges of vertical uplift by applying the above-given formulas 1, 2, and 3 (Figure 3, Table 2 ). In the following, this approach is described for the two terraces, T3 and T5. a Water content expressed as a percentage of the mass of dry sediment, calculated using field values. b Concentrations of 238 U, 232 Th, and 40 K were determined from laboratory measurements using high-resolution gamma spectrometry. c Beta, gamma, and cosmic dose rates were combined to calculate the total dose rates after Aitken [90]. Beta and gamma dose rates were calculated using the conversion factors of Liritzis et al. [79]. Beta dose rates have been corrected for the effect of the water content, grain size, and chemical etching. Gamma dose rates have been corrected for the effect of the water content. Cosmic dose rates were calculated according to Prescott and Hutton (1994) [91]. d The number of aliquots accepted/measured. e Equivalent dose values were determined using the central age model (CAM) of Galbraith et al. [92]. f Ages before 2016 AD. 

Table 2. Calculation of uplift rates of the inner edges of terraces T3 and T5 along Profiles 1 and 3 at the locations where the samples for OSL dating were collected. For the calculations of the errors associated with the position of the sea-level at the time of the formation of the terraces, we used the sea-level curve of Waelbroeck et al. 

[25], while for the uncertainties due to the duration of the MIS, we adopted the time boundaries of the MIS highstands from Lisiecki and Raymo [88]. The OSL age of the sample M collected from the T3 outcrop (topographic profile 3) ranges between 130. 1 and 113. 3 ka, strongly suggesting that T3 was formed during the prominent sea-level highstand corresponding to MIS 5e at 124 ± 6 ka ago (Table 2, Figure 7 ). Terrace T3 is relatively well preserved, and its formation during MIS 5e is supported by the observation that marine terraces formed during that stage are generally well preserved in the geomorphologic record of many tectonically active regions [4, 5, 12, 93, 94]. According to Lisiecki and Raymo [88], MIS 5e began possibly 130 ka ago, and the high sea level was sustained until 118 ka ago and then fell rapidly. Taking into consideration the eustatic sea-level curve of Waelbroeck et al. [25], during this highstand, sea-level reached +6 m above present m. s. l. (with a minimum of -7 m below present-day m. s. l. and a maximum of +19 m above present-day m. s. l. ). As the elevations of the inner edge of this terrace at the sampling site are located at 50 ± 2 m above m. s. l. (Profile 3 in Figure 6 ), the maximum and minimum corrected elevations, taking into account the maximum and minimum eustatic sea-level at its formation time, are 30 and 59 m above m. s. l., respectively. Hence, the resulting maximum and minimum uplift rates would be 0. 50 and 0. 23 mm a -1, respectively, which means a mean long-term uplift rate of 0. 37 ± 0. 13 mm a -1 over the last 124 ± 6 ka (Table 2 ).
section_name: Optically Stimulated Luminescence (OSL) Ages and Tectonic Uplift

section_content: The sedimentary block sample N, collected from a caprock outcrop of marine terrace T5, 2. 3 km northwest of Cape Punta (Profile 1 in Figure 6 ), gave an OSL age of 218. 77 ± 13. 3 Ka. This date indicates that T6 formed during the prominent sea-level highstand of MIS 7c (Table 2, Figure 7 ). This sea-level highstand began 220. 5 ka ago and lasted till about 209. 5 ka ago [88]. The elevation of the inner edge of T5 at the location where the OSL sample was collected is 80 ± 2 m above m. s. l. (Profile 1 in Figure 6 ), and sea-level during the correlative isotope stage 7c was -4 m lower than at present (with a minimum of -17 m below and a maximum of +9 m above present-day m. s. l. ). Taking into consideration the uncertainties regarding the inner edge elevation (±2 m), the eustatic sea-level at the time of the formation of T5, and the duration of the MIS 7c, maximum and minimum uplift rates of 0. 47 and 0. 31 mm a -1, respectively, are estimated. Thus, the long-term mean uplift rate for the last 215 ± 5. 5 ka is 0. 39 ± 0. 08 mm a -1 (Table 2 ). This value is in good agreement with the long-term uplift rate of 0. 37 ± 0. 13 mm a -1 over the last 124 ± 6 ka resulting from the correlation of terrace T3 with MIS 5e. It should be noted that to gain more insight on uplift variations through time, different terrace surfaces within the same profiles must be dated and their uplift rates compared. The sedimentary block sample N, collected from a caprock outcrop of marine terrace T5, 2. 3 km northwest of Cape Punta (Profile 1 in Figure 6 ), gave an OSL age of 218. 77 ± 13. 3 Ka. This date indicates that T6 formed during the prominent sea-level highstand of MIS 7c (Table 2, Figure 7 ). This sea-level highstand began 220. 5 ka ago and lasted till about 209. 5 ka ago [88]. The elevation of the inner edge of T5 at the location where the OSL sample was collected is 80 ± 2 m above m. s. l. (Profile 1 in Figure 6 ), and sea-level during the correlative isotope stage 7c was -4 m lower than at present (with a minimum of -17 m below and a maximum of +9 m above present-day m. s. l. ). Taking into consideration the uncertainties regarding the inner edge elevation (±2 m), the eustatic sea-level at the time of the formation of T5, and the duration of the MIS 7c, maximum and minimum uplift rates of 0. 47 and 0. 31 mm a -1, respectively, are estimated. Thus, the long-term mean uplift rate for the last 215 ± 5. 5 ka is 0. 39 ± 0. 08 mm a -1 (Table 2 ). This value is in good agreement with the long-term uplift rate of 0. 37 ± 0. 13 mm a -1 over the last 124 ± 6 ka resulting from the correlation of terrace T3 with MIS 5e. It should be noted that to gain more insight on uplift variations through time, different terrace surfaces within the same profiles must be dated and their uplift rates compared. 

Along the coastline of the southeastern part of the study area, between Profitis Ilias and Agia Marina, uplifted abrasion platforms up to 2 ± 0. 5 m above present-day m. s. l. are observed (Figures 4 and 8 ). In some cases, these uplifted surfaces are backed by a low, steep cliff with a tidal notch at its base, indicating the position of the recent palaeoshoreline (Figure 8a, b ). Although we did not date these recently raised sea-level indicators, it is evident by their presence that the prevailing late Pleistocene tectonic movement remained positive (uplift) during the Holocene. Along the coastline of the southeastern part of the study area, between Profitis Ilias and Agia Marina, uplifted abrasion platforms up to 2 ± 0. 5 m above present-day m. s. l. are observed (Figures 4 and 8 ). In some cases, these uplifted surfaces are backed by a low, steep cliff with a tidal notch at its base, indicating the position of the recent palaeo-shoreline (Figure 8a, b ). Although we did not date these recently raised sea-level indicators, it is evident by their presence that the prevailing late Pleistocene tectonic movement remained positive (uplift) during the Holocene. 

Mean long-term uplift rates estimated in this study are higher than the mean uplift rates estimated for the study area by Kourampas [54], who focused on the study of the late Pleistocene sedimentary facies at the broader region of the Lakonic Peninsula. He suggested a long-term uplift rate ranging from 0. 08 to 0. 33 mm a -1 during the last 1 Ma. It should be noted that the results of our work are in accordance with regional upper Quaternary uplift rates proposed for Crete Island, a tectonically active region very close to the Hellenic Subduction Zone (Figure 1a ). Along the coast of the western part of Crete, as well as along the western-central coasts of the southern part of the island, Pleistocene rock uplift rates range between 0. 5 and 1. 1 mm a -1, whereas for the eastern coast of the southern part of the island a mean uplift rate < 0. 5 mm a -1 is suggested [8, 10, 20]. In this area, the role of the local active faults is quite important and differentiates the observed uplift rates from place to place. Mean long-term uplift rates estimated in this study are higher than the mean uplift rates estimated for the study area by Kourampas [54], who focused on the study of the late Pleistocene sedimentary facies at the broader region of the Lakonic Peninsula. He suggested a long-term uplift rate ranging from 0. 08 to 0. 33 mm a -1 during the last 1 Ma. It should be noted that the results of our work are in accordance with regional upper Quaternary uplift rates proposed for Crete Island, a tectonically active region very close to the Hellenic Subduction Zone (Figure 1a ). Along the coast of the western part of Crete, as well as along the western-central coasts of the southern part of the island, Pleistocene rock uplift rates range between 0. 5 and 1. 1 mm a -1, whereas for the eastern coast of the southern part of the island a mean uplift rate < 0. 5 mm a -1 is suggested [8, 10, 20]. In this area, the role of the local active faults is quite important and differentiates the observed uplift rates from place to place.
section_name: Sample

section_content: To assign an age to the undated marine terraces, we have tried to correlate them with past sea-level highstands observed in the eustatic curve [25] using a mean uplift rate of 0. 36 ± 0. 11 mm a -1, estimated from the dated terraces T3 and T5, based on the assumption that this long-term rate remains temporally constant during the late Pleistocene. This approach is consistent with many similar previous studies that linked marine terraces to past sea-level highstands in tectonically active regions [6, 10, 14, 20, 95, 96]. Assuming a
section_name: Correlation of the Undated Marine Terraces with Past Sea-Level Highstands

section_content: To assign an age to the undated marine terraces, we have tried to correlate them with past sea-level highstands observed in the eustatic curve [25] using a mean uplift rate of 0. 36 ± 0. 11 mm a -1, estimated from the dated terraces T3 and T5, based on the assumption that this long-term rate remains temporally constant during the late Pleistocene. This approach is consistent with many similar previous studies that linked marine terraces to past sea-level highstands in tectonically active regions [6, 10, 14, 20, 95, 96]. Assuming a constant uplift rate, we estimated the predicted (or expected) elevation ranges of the inner edges for each MIS marine terrace (Table 3 ). This approach is justified by the fact that the ages of the two OSL dated samples gave similar uplift rates over the last ~125 and ~220 ka. The predicted elevations were compared with the observed real inner edge elevations of the marine terraces of the study area, and for each terrace, an MIS was assigned (Table 3, Figures 9 and 10 ). 

The lowest and thus youngest marine terrace T1, shows good overall morphological continuity, particularly along the north and west shores of Elafonissos Island, as well as along the opposite coast between Cape Pounta and Asproudia (Figure 4 ). The widest surfaces of this terrace are developed on marine Pleistocene deposits north and northwest of Cape Punta, and at the southeastern part of the Elafonissos Island, T1 is bordered by an inner edge with an elevation of 8 ± 2 m above m. s. l. By applying the above estimated average long-term uplift rate of 0. 36 ± 0. 11 mm a -1 to MIS 5a, the younger MIS 5 highstand at 81 ± 5 ka, and taking into consideration that this highstand was lower than the present mean sea-level by about -19 ± 13 m [25], a predicted elevation of 12 ± 4 m is obtained for the inner edge. The observed present-day elevation of the T1 inner edge (8 ± 2 m above m. s. l. ) falls within the predicted elevation range of the MIS 5a marine terrace palaeoshoreline angle (Table 3, Figure 9 ). Hence we can assign an MIS 5a age to this terrace (Figure 10 ). This terrace probably corresponds to the surface carved into the Pliocene-Pleistocene formations, recognized at the broader area of Neapolis, Profitis Ilias, and Agia Marina by previous researchers [51, 52, 97] at elevations ranging from 4 to 7 m above m. s. l. It was correlated by them with the "Neo-Tyrrhenian" sea-level cycle, which corresponds to MIS 5a; ca. 80 ka [54, 98], which is in agreement with our proposal. Table 3. Name, assigned MIS, age, and observed present-day elevations of the inner edges of the marine terraces in three different locations of the study area. The predicted (or expected) elevation ranges of the inner edges (palaeo-shorelines) for each MIS were calculated under the assumption that the long-term uplift rate (0. 36 ± 0. 11 mm a -1 ) was constant over the last 401 ± 10 ka and taking into account the duration of each MIS (according to Lisiecki and Raymo [88], as well as the maximum and minimum position of eustatic sea-level for each MIS highstand (according to sea-level curve by Waelbroeck et al. [25]. constant uplift rate, we estimated the predicted (or expected) elevation ranges of the inner edges for each MIS marine terrace (Table 3 ). This approach is justified by the fact that the ages of the two OSL dated samples gave similar uplift rates over the last ~125 and ~220 ka. The predicted elevations were compared with the observed real inner edge elevations of the marine terraces of the study area, and for each terrace, an MIS was assigned (Table 3, Figures 9 and 10 ). The lowest and thus youngest marine terrace T1, shows good overall morphological continuity, particularly along the north and west shores of Elafonissos Island, as well as along the opposite coast between Cape Pounta and Asproudia (Figure 4 ). The widest surfaces of this terrace are developed on marine Pleistocene deposits north and northwest of Cape Punta, and at the southeastern part of the Elafonissos Island, T1 is bordered by an inner edge with an elevation of 8 ± 2 m above m. s. l. By applying the above estimated average long-term uplift rate of 0. 36 ± 0. 11 mm a -1 to MIS 5a, the younger MIS 5 highstand at 81 ± 5 ka, and taking into consideration that this highstand was lower than the present mean sea-level by about -19 ± 13 m [25], a predicted elevation of 12 ± 4 m is obtained for the inner edge. The observed present-day elevation of the T1 inner edge (8 ± 2 m above m. s. l. ) falls within the predicted elevation range of the MIS 5a marine terrace palaeoshoreline angle (Table 3, Figure 9 ). Hence we can assign an MIS 5a age to this terrace (Figure 10 ). This terrace probably corresponds to the surface carved into the Pliocene-Pleistocene formations, recognized at the broader area of Neapolis, Profitis Ilias, and Agia Marina by previous researchers [51, 52, 97] at elevations ranging from 4 to 7 m above m. s. l. It was correlated by them with the "Neo-Tyrrhenian" sea-level cycle, which corresponds to MIS 5a; ca. 80 ka [54, 98], which is in agreement with our proposal. 

The next higher marine terrace, T2, is relatively well preserved and shows good morphological continuity, following the present coastline. The widest surface of this terrace is carved into marine Pleistocene deposits between Profitis Ilias and Agia Marina, The next higher marine terrace, T2, is relatively well preserved and shows good morphological continuity, following the present coastline. The widest surface of this terrace is carved into marine Pleistocene deposits between Profitis Ilias and Agia Marina, west of Cape Maleas (Figure 4 ). The current elevation of the shoreline angle of T2 varies between 16 ± 2 and 24 ± 2 m above m. s. l. For this terrace, there are no dates. However, an age corresponding to MIS 5c of 101 ± 5 ka was adopted as the most appropriate. By using the mean uplift rate estimated by the dated terraces and considering that during the MIS 5c sea-level was located at -21 m (with a minimum of -34 m and a maximum of -8 m below present m. s. l. [25]), the elevation of the inner edge of the terrace formed during MIS 5c should be at 15 ± 4 m above m. s. l. which is in agreement with the lower limit of the observed inner edge elevation of T2 (Figure 9 ). This terrace is possibly the surface reported by Kelletat et al. [52] at elevations ranging between 8 and 20 m. They assigned an "Eu-Tyrrhenian" age (which corresponds to MIS 5e; ca 124 ± 6 ka) to this surface. However, we propose the age of MIS 5c (101 ± 5 ka) as appropriate for T2 as we correlated the higher terrace T3 with MIS 5e sea-level highstand based on the OSL date of sample M1. 

Terrace T4 is a discontinuous surface that was not identified at the broader area of Profitis Ilias-Agia Marina during the field survey (Figure 4 ). The absence of T4 palaeo-surfaces from this area can be attributed to its destruction by wave activity during subsequent sea-level highstands, when the landward retreat of sea-cliffs formed younger and lower terraces, or to the morphology of the area at the time of its formation. Where present, the elevation of the inner edges of T4, is observed between 64 ± 2 and 68 ± 2 m above m. s. l. and remains roughly constant along the entire coast of the study area. T4 most likely correlates with the MIS 7a sea-level highstand at 197 ± 6 ka. Assuming the 0. 36 ± 0. 11 mm a -1 long-term uplift rate value as constant and that during MIS 7a sea-level was lower than the present m. s. l. by about -10 ± 13 m [25], an elevation of 61 ± 6 m is expected for the marine terrace of this highstand (Table 3, Figure 9 ) which fits well with the observed elevations of T7 inner edges. 

Marine terrace T6 is generally a discontinuous surface, not recorded at the southeastern part of the study area (Profitis Ilias-Agia Marina) (Figure 4 ). This terrace was presumably destroyed in this area due to surface erosional processes after the time of its formation. Along the western coast of Elafonissos Island and northwest of Cape Punta (Profiles 2 and 3 of Figure 6 ), T6 has the shape of a narrow, relatively morphologically continued bedrock marine terrace eroded into Mesozoic limestones and dolomites. T6 surfaces are bordered by inner edges extending at elevations between 100 ± 2 and 108 ± 2 m above m. s. l. For this terrace, we assign an age of 287 ± 5 ka, which corresponds to the previous prominent highstand MIS 9a (Figure 10 ). Considering the long-term tectonic uplift rate of the area and that this highstand was about 16 ± 13 m below present m. s. l., this implies that the inner edges of the terrace that corresponds to MIS 9a should be located at elevations of about 87 ± 17 m, which is not far from the elevations of the inner edges measured for this terrace (Figure 9 ). The small difference in elevation among the predicted and the measured inner edge elevation may be due to a possible slight change in the uplift rate or due to the long period of erosion of the inner edges of the terrace. 

Marine terrace T7 is represented by remnants of extensive, gently seaward-dipping surfaces carved mainly into Mesozoic bedrock (limestones, dolomites, and metamorphic rocks) at the eastern part of the study area and partly into limestones and dolomites on Elafonissos Island. In the area north of Cape Punta, these platforms are exclusively developed on Pleistocene marine sedimentary formations (Figures 2 and 4 ). T7 palaeoplatform is bounded landwards by inner edges of elevations ranging from 128 ± 2 to 140 ± 2 m above m. s. l. An age of 330 ± 4 ka, which corresponds to the MIS 9e, is suggested for this terrace (Figure 10 ). Assuming that the study area has been uplifted with a constant mean rate of 0. 36 ± 0. 11 mm a -1, and considering that sea level during this period was +5 m above the present-day m. s. l. (with a minimum of -8 m below and a maximum of +14 m above m. s. l. ) [25], the elevation of the inner edge of the terrace of this highstand should be between 103 m and 145 m (Table 3, Figure 9 ). This predicted elevation coincides with the measured inner edge elevations of T7. 

Marine terrace T8 is characterized by extensive, relatively plain surfaces developed on marine Pleistocene deposits, as well as on carbonate and metamorphic rocks of the Mesozoic bedrock (Figures 2 and 4 ). T8 was not identified in the eastern part of the study area. The palaeoshoreline angle of this terrace is well defined and has been recognized at elevations between 152 ± 2 and 156 ± 2 m above m. s. l. The MIS 11c at 407 ± 10 ka was adopted for T8 by applying the previous reasoning (Figure 10 ). Considering the constant mean long-term tectonic uplift rate proposed by Waelbroeck et al. [25] sea-level position at the interglacial period of MIS 11c (+6 m above m. s. l. -between -6 m below and +19 m above m. s. l. ), a correlation with MIS 11c seems reasonable. The inner edge of the marine terraces of this highstand should be at 150 ± 27 m, which is similar to the observed inner edge elevation of T8 (Table 3, Figure 9 ). T9, the second-highest marine terrace, exists in Elafonissos Island and at the opposite coast between Cape Punta and Cape Koulenti (Figure 4 ). It is developed exclusively on the Mesozoic bedrock, and its inner edges were recorded at 172 ± 2 m above m. s. l. On top of Elafonissos Island, the highest and oldest surface (T10) occurs with the inner edge at about 192 ± 2 m above m. s. l. (Figure 4 ). T10 is not present on the opposite continental coasts of the study area. The morphological discontinuity of these 2 highest terraces (T9 and T10), in addition to their heavily eroded inner edges and the lack of highly reliable global sea-level curves before 450 ka, made precise mapping and correlation difficult (Figure 10 ). Probably they were formed during the MIS sea-level highstands before 407 ± 10 ka. 

As already mentioned, previous studies [52, 54] have reported a sudden southward inclination and submergence of the older marine palaeo-platforms, whereas the terraces are supposedly higher on Elafonisos Island than on the opposite Punta-Neapolis area. The results of our approach are not in agreement with this trend. The fact that the two successfully OSL dated marine terraces had different ages and were collected from different locations of the study area, but gave similar rates of long-term tectonic uplift, suggests that the area is being emerged uniformly in both space and time, at least for the past ~400 ka. Another strong evidence for the almost uniform tectonic uplift of the region is that according to the results of the analysis of the detailed DEM and their validation, where possible, in the field, we detected a uniform morphological distribution of the surfaces and inner edges of each marine terrace at similar elevation ranges throughout the study area. The spatially uniform uplift around Neapolis is very different from the well-known examples in Crete and Corinth [4, 6, [8] [9] [10] [11] 20, 29, [42] [43] [44] [45] [46] [47] [48] [49] [50], where terrace elevations are heavily influenced by local active faults on and offshore and vary in elevations over distances of a few kilometers. The DEM of the area suggests there is not much active faulting onshore, but the uniform terrace elevations give the suggestion that, on a broader scale, the uplift is rather a subduction signal than an offshore normal fault signal. Hence, the driving factor of the tectonic uplift is associated with the active subduction of the African lithosphere beneath the Eurasian plate since the area is situated very close (~90 km) to the active margin of the Hellenic subduction zone.
section_name: Correlation of the Undated Marine Terraces with Past Sea-Level Highstands

section_content: The application of GIS analysis, DGPS techniques, and conventional field-mapping methods along the coasts of the broader area of Neapolis in southeastern Peloponnese, enabled us to recognize and map a sequence of 10 tectonically uplifted marine terraces with inner edges ranging in elevations from 8 ± 2 to 192 ± 2 m above m. s. l. The presence of marine deposits of Pliocene-Lower Pleistocene age up to the elevation of ~180 m above m. s. l., along with the presence of the flight of preserved marine terraces, demonstrate that the area was affected by a vigorous tectonic uplift during the late Quaternary. Additionally, the presence of recently raised wave-cut shore platforms up to +2. 5 m above m. s. l. indicates that the positive tectonic movement of the area continues during the Holocene. 

Two sedimentary block samples collected from 2 of the marine terraces (T3 and T5) have been OSL dated and gave burial ages (with associated standard errors) of 121. 7 ± 8. 4 and 218. 8 ± 13. 3, respectively. These OSL ages yielded constraints for correlating T3 and T5 to the prominent sea-level highstand of MIS 5e (124 ± 6 ka) and 7c (215 ± 5. 5 ka) and to estimate a mean long-term tectonic uplift rate of 0. 36 ± 0. 11 mm a -1 over the last 215 ± 5. 5 ka. However, to gain more insight on uplift variations through time, more terrace levels must be dated and their uplift rates compared. 

The assumption that the long-term uplift rate remained temporally constant during the Pleistocene permitted the correlation of six of the undated terraces of the sequence with the past sea-level highstands observed in the global eustatic curve. The lowermost Terrace, T1, can be correlated to sea-level highstand MIS 5a, whereas T1, T2, T4 T6, T7, and T8 were correlated with MIS 5a, 5c, 7a, 9a, 9e, and 11c sea-level highstands, respectively. We were not able to correlate the two higher marine terraces, which probably were formed before 401 ± 1 ka ago. 

The results of this study showed that the broader Neapolis area is being emerged uniformly in both space and time, at least for the past ~400 ka, while the driver of tectonic uplift is the active subduction of the African lithosphere beneath the Eurasian plate.
section_name: Conclusions

section_name: Acknowledgments:

section_content: We would like to thank Gino de Gelder and an anonymous reviewer for their comments and suggestions that significantly improved the final version of the paper.
section_name: Acknowledgments:

section_name: Data Availability Statement:

section_content: The data presented in this study are available on request from the corresponding author.
section_name: Data Availability Statement:

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: