Coastal Boulders on the SE Coasts of Cyprus as Evidence of Palaeo-Tsunami Events

Cyprus has a long history of tsunami events, as noted by archaeological and geological records. At Cape Greco (southeastern Cyprus) large boulders have been noted, however, no detailed geomorphological research has taken place so far and the related high energy event was undated until now. Our research aims to record in detail and interpret these large boulders deposits. The boulders, located between ≈3 and 4.5 m a.m.s.l., are fragments of an upper Pleistocene aeolianite, which is overlaying unconformly a lower Pleistocene calcarenite. Dimensions and spatial distribution of 272 small, medium, and large boulders were documented, while their precise distance from the coastline was recorded by field mapping and remote sensing, using Differential GPS (DGPS), drone, and Geographic Information Systems (GIS) technics. Field data were subsequently combined with hydrodynamic equations, in order to determine the extreme event(s) that caused their transport inland, and radiocarbon dating was accomplished on three samples of Vermetus sp. to determine the chronological context. Our findings appear to broadly correlate with the 1303 AD tsunami, which has displaced at least part of the studied boulders, and one other undocumented event at AD 1512-1824. The large number of boulders and sizes in our study area further indicate that their dislocation is most likely owed to multiple events from various sources.


Introduction
Large boulder accumulations along the coasts have been widely used in the literature to reconstruct and determine the high energy event that displaced them. Such boulders or mega-blocks are often found isolated or in groups bearing evidence of their transportation from their original location. Characteristics include their placement in altitudes of few meters above mean sea level (a.m.s.l.), individually, or in clusters, with the presence of features that certify their transfer from the intertidal or supratidal zone. Boulder deposits have been linked with storm events [1][2][3][4] as well as tsunamis [5][6][7][8][9][10][11].
The interpretation of the high energy events, capable of dislocating large boulders, remains controversial [12] and have been interpreted as both "storm" and "tsunami" deposits and they often The interpretation of the high energy events, capable of dislocating large boulders, remains controversial [12] and have been interpreted as both "storm" and "tsunami" deposits and they often display similar characteristics [13]. The investigation of modern tsunami deposits can provide key features to better distinguish such events [14], such as the 2004 Indian ocean tsunami [7,15], or the 2011 Japan tsunami [16][17][18]. Landward imbrication has been noted by various researchers from recent tsunami deposits [17,19], although this feature may be also found in storm deposits. According to Etienne and Paris [2] elongated boulders are characterized by their imbrication or long axis being tangent to the tsunami wave train, and thus allow to determine the direction of the wave. Tsunami boulders further do not display a particular trend in size distribution along the flow direction [18,20].
Different hydrodynamic equations [18,[21][22][23][24] have been developed to decipher the origin of the deposits through the calculation of wave heights and velocities necessary to transport them. The micro-tidal regime and limited fetch in the Mediterranean is ideal for the application of such equations, as maximum possible storm wave heights are smaller in comparison to areas where significant storm wave heights may occur [20]. When the calculated storm wave heights differ greatly from the wave-storm regime of the study area, then a tsunamigenic origin is deduced [6].
Research on palaeo-high energy events is fundamental to understand the coastal responses of tsunami or storm impact, for estimating the impact of past events for coastal hazard mitigation. Cyprus has a large history of seismic wave activity as described by archaeological and geological records. There is abundance of historical records in literature, which are linked to the development of tsunamis in the Eastern Mediterranean [25,26]. The presence of various geomorphological evidence indicative of tsunamis has been noted by several researchers [27,28] in various areas in Cyprus (e.g., Kormakitis peninsula, Gialoussa coast, Cape Greco, Ayia Napa, Pafos airport, Kissonergas, Lara and Akamas peninsula) ( Figure 1). The most interesting indicators are boulder deposits, which may be found in various locations in Cyprus, primarily on the west and southeast coasts ( Figure 1) [27,28]. In this paper, we focus on the boulders found on the coastal zone of Cape Greco. Their presence has been previously noted by Whelan and Kelletat [28], estimating a weight of 10-20 tons.
The aim of our paper is to demonstrate that the boulders on the coastal zone of Cape Greco were transported and deposited by tsunamis and to provide a chronological framework, in order to correlate them with known earthquakes and tsunamis in the wider area.  [27][28][29]). At the southeastern part is Cape Greco, the study area. Base map sources: Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community.  [27][28][29]). At the southeastern part is Cape Greco, the study area. Base map sources: Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community.

Geological Setting
The island of Cyprus is located in the south east part of the Mediterranean Sea. The morphology of the southern coast of Cyprus is complex, being partly a result of Miocene-Quaternary tectonic uplift reaching 5 m/ka [30]. A relatively simple straight coastline profile is observed southwest of Larnaca Bay. This simple morphology changes gradually into a rocky, tectonically uplifted shoreline with several bays and pools towards the east. The presence of a low-lying rocky shoreline continues through Cape Greco until the shores of Ammochostos Bay [31].
Our study area, Cape Greco, lies at the south-eastern part of the island and forms a narrow peninsula between Ayia Napa and Paralimni area ( Figure 1). It mainly consists of Pliocene marls of the Nicosia Formation and Miocene Terra reefal limestone sandshales, part of the Pakhna Formation [32]. These formations are overlain by calcarenites and beachrocks, which are abundant in Pleistocene terrace deposits [31]. The studied boulders are located on an almost horizontal rocky platform of Pleistocene calcarenites, at a mean altitude of ≈3 m (Figures 2 and 3).  (c) (d) Figure 3. Characteristic boulder deposits from the study area at Cape Greco. (a), (b) Two of the largest boulders in the study area; the boulder of (a) is also visible in the background of (c) and (d).

Seismicity
The easternmost part of the Mediterranean Sea is seismotectonically dominated by the Cyprean arc and the left-lateral strike-slip Levantine rift ( Figure 4). In the Cyprean arc shallow-and intermediate-depth earthquakes occur mostly in the submarine environment and, therefore, highmagnitude shocks are expected to generate strong tsunamis [25].
The spatial distribution of tsunami events in the Cyprus-Levantine Sea region is not random and relates to the tectonic setting of the wider region. According to Fokaefs and Papadopoulos [25], they are concentrated along two zones: (a) Zone of south Cyprus, where earthquake activity follows the south part of Cyprus along the Cyprean arc and the seismogenic sources under the sea favors the generation of tsunamis by co-seismic fault dislocations and (b) Levantine Sea from Gaza northward zone in which causative earthquakes are mainly associated with the left-lateral strike-slip fault system of the Levantine rift, which implies that the seismogenic sources are not located in a submarine environment but on land and, therefore, the tsunami generation mechanism remains unexplained. Characteristic boulder deposits from the study area at Cape Greco. (a, b) Two of the largest boulders in the study area; the boulder of (a) is also visible in the background of (c) and (d).

Seismicity
The easternmost part of the Mediterranean Sea is seismotectonically dominated by the Cyprean arc and the left-lateral strike-slip Levantine rift ( Figure 4). In the Cyprean arc shallow-and intermediate-depth earthquakes occur mostly in the submarine environment and, therefore, high-magnitude shocks are expected to generate strong tsunamis [25].
The spatial distribution of tsunami events in the Cyprus-Levantine Sea region is not random and relates to the tectonic setting of the wider region. According to Fokaefs and Papadopoulos [25], they are concentrated along two zones: (a) Zone of south Cyprus, where earthquake activity follows the south part of Cyprus along the Cyprean arc and the seismogenic sources under the sea favors the generation of tsunamis by co-seismic fault dislocations and (b) Levantine Sea from Gaza northward zone in which causative earthquakes are mainly associated with the left-lateral strike-slip fault system of the Levantine rift, which implies that the seismogenic sources are not located in a submarine environment but on land and, therefore, the tsunami generation mechanism remains unexplained.

Significant Wave Height
Wave heights in Cyprus waters are generally lower than in the large ocean basins because of generally weaker storms. Published data for the period 1961-1980 from ship observations near Cape Greco suggests average significant wave heights between 0.25 and 0.75 m [44]. According to the same dataset, extreme wave heights with a modelled 100-years return period reach 4.7-5 m, with SSE to SSW direction [44]. Buoy data for the period between July 2005 to February 2008 at an open ocean sea

Significant Wave Height
Wave heights in Cyprus waters are generally lower than in the large ocean basins because of generally weaker storms. Published data for the period 1961-1980 from ship observations near Cape Greco suggests average significant wave heights between 0.25 and 0.75 m [44]. According to the same dataset, extreme wave heights with a modelled 100-years return period reach 4.7-5 m, with SSE to SSW direction [44]. Buoy data for the period between July 2005 to February 2008 at an open ocean sea point west of the Eratosthenes Seamount suggest significant wave heights less than 1.5 m, with a frequency percentage of ≈94% [45]. The highest recorded value is 4.75, with a frequency percentage of 0.0167% for the same time period.

Fieldwork
A geomorphological survey was conducted along the coastal zone of Cape Greco (Figures 2 and 3). Detailed mapping of the boulders was realized in the field using a hand-held GPS and DGPS-GNSS system for their precise location. Topographic profiles were performed perpendicular to the coast, using DGPS-GNSS ( Figure 5). The measurement of all boulder dimensions was realized using a measuring tape, i.e., a-b-c axis (length-width-height). The a-axis direction was also recorded, their elevation from sea level (m), and their distance from the coastline. Aerial photography was accomplished using a DJI mavic pro drone. An orthomosaic was developed using Agisoft PhotoScan by combining 840 aerial photographs taken from the drone. point west of the Eratosthenes Seamount suggest significant wave heights less than 1.5 m, with a frequency percentage of ≈94% [45]. The highest recorded value is 4.75, with a frequency percentage of 0.0167% for the same time period.

Fieldwork
A geomorphological survey was conducted along the coastal zone of Cape Greco (Figures 2 and  3). Detailed mapping of the boulders was realized in the field using a hand-held GPS and DGPS-GNSS system for their precise location. Topographic profiles were performed perpendicular to the coast, using DGPS-GNSS ( Figure 5). The measurement of all boulder dimensions was realized using a measuring tape, i.e., a-b-c axis (length-width-height). The a-axis direction was also recorded, their elevation from sea level (m), and their distance from the coastline. Aerial photography was accomplished using a DJI mavic pro drone. Αn orthomosaic was developed using Agisoft PhotoScan by combining 840 aerial photographs taken from the drone.

Radiocarbon Dating
A thorough field investigation was accomplished in order to locate fossils (for 14 C dating) that may reflect when the boulders were transported. Marine fossils that may reflect when the boulders were transported are very rare in the study area, and the investigated boulders are composed of calcarenite. Nevertheless, three samples of marine gastropods (Vermetus sp.) were collected from the studied boulders for radiocarbon dating. These marine gastropods suggest an offshore to inland transportation during the dislocation of the boulders from their in-situ position and, therefore, they may constrain the age of the high energy events.
Radiocarbon dating was carried out at CEntro di DAtazione e Diagnostica (CEDAD), in Lecce, Italy, on three samples of Vermetus sp. Radiocarbon ages were calibrated through the online software Calib 7.10 [46], with the Marine13 curve. The samples were corrected for the local marine reservoir effect using a mean ΔR value of 53±43 for the Eastern Mediterranean [47].

Estimation of Boulder Weight and Volume
For the determination of the boulders' mass (M), volume (V), and density (ρb), we used the water displacement method, in the laboratory, in two samples collected from the calcarenite and the

Radiocarbon Dating
A thorough field investigation was accomplished in order to locate fossils (for 14 C dating) that may reflect when the boulders were transported. Marine fossils that may reflect when the boulders were transported are very rare in the study area, and the investigated boulders are composed of calcarenite. Nevertheless, three samples of marine gastropods (Vermetus sp.) were collected from the studied boulders for radiocarbon dating. These marine gastropods suggest an offshore to inland transportation during the dislocation of the boulders from their in-situ position and, therefore, they may constrain the age of the high energy events.
Radiocarbon dating was carried out at CEntro di DAtazione e Diagnostica (CEDAD), in Lecce, Italy, on three samples of Vermetus sp. Radiocarbon ages were calibrated through the online software Calib 7.10 [46], with the Marine13 curve. The samples were corrected for the local marine reservoir effect using a mean ∆R value of 53±43 for the Eastern Mediterranean [47].

Estimation of Boulder Weight and Volume
For the determination of the boulders' mass (M), volume (V), and density (ρb), we used the water displacement method, in the laboratory, in two samples collected from the calcarenite and the subaerially-eroded red horizon, respectively. The data from the two samples were combined with the measurements of the three major axes to determine the weight of the boulders. For the calculation of the boulders volume, the maximum length of a-, b-, and c-axes, i.e., length, width, and height, respectively, was used: V = a × b × c.

Hydrodynamic Equations
A numerical approach was used to quantify the wave height (storm or tsunami) necessary for the detachment and transportation of the boulders inland. The hydrodynamic equations for storm and tsunamis and the boulder transportation scenarios were initially derived by Nott [21]. The developed hydrodynamic equations are covering a part of all likely scenarios for boulder transportation. The scenarios relate to the possible original position of the boulder. Thus, the main scenarios are focusing on a possible submerged boulder prior to transportation called the "submerged boulder (SMB) scenario". A subaerial boulder scenario (SAB) is considered for cases when the boulders are moved by waves on shore. Finally, a joint bounded block scenario (JBB) was taken into account in order to include a possible movement of boulders that were dislocated from the rocky coast or for boulders that resulted from a rock fall on a sea-cliff. For our case study we focused on the hydrodynamic approach by Pignatelli et al. [22] and Nandasena et al. [48] for the joint-bounded block scenario (JBB) and submerged boulder scenario (SMB). For comparison on the resulting wave height values, other hydrodynamic approaches were considered [23,24,49].
For our calculations, a mean coefficient of drag (C d ) of 1.2 was considered, according to Helley [50] and Benner et al. [23] through a graph based on the shape factor (SF) of the boulders. Furthermore, since the boulders are composed of calcareous material, the coefficient of mass (C m ) was given a value of 2. Additionally, we used a coefficient of lift (C L ) of 0.178 according to Noormets et al. [51], coefficient of friction of the boulders (µ s ) of 0.65, the instantaneous flow acceleration of the boulder as 1 m s −2 [23], and the sea water density of 1.02 g cm −3 .

Field Investigation and Boulders Distribution
During field work, a total number of 272 boulders were identified and studied ( Figure 6). The measured density from the two main lithologies in the study area was calculated at 1.94 × 10 3 kg m −3 for the calcarenite and 1.97 × 10 3 kg m −3 for subaerially-eroded red horizon. Based on the boulders volume calculations, 154 have sizes less than 2 m 3 , 81 boulders up to 6 m 3 , 24 boulders up to 10 m 3 , and 13 very large boulders over 10 m 3 with a remarkable note of the largest up to 41.68 m 3 ( Figure 6).
During field work, a total number of 272 boulders were identified and studied ( Figure 6). The measured density from the two main lithologies in the study area was calculated at 1.94 × 10 3 kg m −3 for the calcarenite and 1.97 × 10 3 kg m −3 for subaerially-eroded red horizon. Based on the boulders volume calculations, 154 have sizes less than 2 m 3 , 81 boulders up to 6 m 3 , 24 boulders up to 10 m 3 , and 13 very large boulders over 10 m 3 with a remarkable note of the largest up to 41.68 m 3 ( Figure 6).  The vast majority of the boulders are located at an altitude between 3 and 4.6 m a.m.s.l., with the exception of 10 boulders lying at the erosional platform at about 1 m elevation (Figures 7 and 8). Their distance from the coastline varies between 5 and 100 m, but most are arranged at a distance of 70 m from the shoreline (Figures 7 and 8). Two main clusters of boulder distribution may be recognized, one at the western part and the second at the eastern part of the study area, while in the central part only few boulders were found sporadically. On the eastern part, a notable accumulation of boulders was in imbricated clusters ( Figure 9). The vast majority of the boulders are located at an altitude between 3 and 4.6 m a.m.s.l., with the exception of 10 boulders lying at the erosional platform at about 1 m elevation (Figures 7 and 8). Their distance from the coastline varies between 5 and 100 m, but most are arranged at a distance of 70 m from the shoreline (Figures 7 and 8). Two main clusters of boulder distribution may be recognized, one at the western part and the second at the eastern part of the study area, while in the central part only few boulders were found sporadically. On the eastern part, a notable accumulation of boulders was in imbricated clusters ( Figure 9).

Wave Height Calculation
Taking into account the likely boulder transportation scenarios (see Section 3.4) and based on field observations, a scenario of joint-bounded blocks (JBB) was considered for boulders located at areas 1 and 2b, while for the boulders of the western part (area 3) a scenario of submerged boulder (SMB) was additionally taken in account. Considering the scenario and using the hydrodynamic equations by Pignatelli et al. [22] and Nandasena et al. [48], the wave heights for storm and tsunami capable of transporting boulders, are presented in Table 2. Table 2 further includes calculations derived from Barbano et al. [49], Benner et al. [23], and Engel and May [24], for comparison.
According to hydrodynamic equations of Pignatelli et al. [22] for a JBB scenario, the displacement of the largest boulder could have taken place by a minimum storm wave (Hs) of ≈19.      Table S1.

Five Largest Values of Storm Wave Height for JBB Scenario SAB Scenario Boulder
No Axes (m) Volume (m 3 )

Radiocarbon Results
Amongst the three dated shell samples, sample LTL19211A was not considered as it provided an age younger than 1950 AD. The calibrated ages of the two other samples (LTL19209A, LTL19210A) are, respectively, 1512-1824 AD and 1309-1496 AD (Table 3).

Wave Height Calculation
Taking into account the likely boulder transportation scenarios (see Section 3.4) and based on field observations, a scenario of joint-bounded blocks (JBB) was considered for boulders located at areas 1 and 2b, while for the boulders of the western part (area 3) a scenario of submerged boulder (SMB) was additionally taken in account. Considering the scenario and using the hydrodynamic equations by Pignatelli et al. [22] and Nandasena et al. [48], the wave heights for storm and tsunami capable of transporting boulders, are presented in Table 2. Table 2 further includes calculations derived from Barbano et al. [49], Benner et al. [23], and Engel and May [24], for comparison. Table 2. Selected boulder measurements and storm wave height calculations based on different hydrodynamic equations [22][23][24]48,49]. All boulder measurements and calculations are provided in Table S1.   According to hydrodynamic equations of Pignatelli et al. [22] for a JBB scenario, the displacement of the largest boulder could have taken place by a minimum storm wave (H s ) of ≈19.57 m or a tsunami wave (H t ) of ≈4.89 m, while according to Nandasena et al. [48] a minimum storm of ≈15.24 m or a tsunami wave of ≈3.81 m would be necessary. In case of the boulders that were submerged at the t0 of the wave impact, a boulder of a-axis 230 cm, b-axis 210 cm, and c-axis 80 cm, would need a minimum storm wave (H s ) of ≈3.6 m or a tsunami wave (H t ) of ≈0.90 m. In agreement with Piscitelli et al. [3], the equations by Benner et al. [23] and Engel and May [24] provide very high values.

Radiocarbon Results
Amongst the three dated shell samples, sample LTL19211A was not considered as it provided an age younger than 1950 AD. The calibrated ages of the two other samples (LTL19209A, LTL19210A) are, respectively, 1512-1824 AD and 1309-1496 AD (Table 3).

Geomorphological Evidence and Boulder Origin
The geomorphological context of the study area suggests a complex pre-transport location for the studied boulders. Erosional and weathering processes are active in the area, predominantly owed to sea waves and sea spray. On the lower part of the rocky platform, near sea level, a well-developed wave cut notch exists and it appears that the origin of some boulders of area 1 is owed to the collapse and detachment of the notch's roof, due to the force of waves. On the same site (area 1), between the calcarenite and the aeolianite formations (Figure 10a,b), a subaerially-eroded red horizon was identified, which appears to have acted as a detachment zone. This is particularly characteristic on boulder 189, which is evident that has been turned upside down, and its original location is clearly visible (Figure 8b), supporting a joint-bounded scenario. Some large boulders have irregular shapes and sharp edges, as well as micro-karstic solution on their surface, indicating that they most likely derived from the edge of the rocky platform. Marine bio-constructions are scarce in the studied boulders, which further support a supratidal origin.

Geomorphological Evidence and Boulder Origin
The geomorphological context of the study area suggests a complex pre-transport location for the studied boulders. Erosional and weathering processes are active in the area, predominantly owed to sea waves and sea spray. On the lower part of the rocky platform, near sea level, a well-developed wave cut notch exists and it appears that the origin of some boulders of area 1 is owed to the collapse and detachment of the notch's roof, due to the force of waves. On the same site (area 1), between the calcarenite and the aeolianite formations (Figure 10a,b), a subaerially-eroded red horizon was identified, which appears to have acted as a detachment zone. This is particularly characteristic on boulder 189, which is evident that has been turned upside down, and its original location is clearly visible (Figure 8b), supporting a joint-bounded scenario. Some large boulders have irregular shapes and sharp edges, as well as micro-karstic solution on their surface, indicating that they most likely derived from the edge of the rocky platform. Marine bio-constructions are scarce in the studied boulders, which further support a supratidal origin. A particular cluster of very large boulders and an accumulation of smaller ones lies in a heavily eroded area (area 2a), the only part of the study area at a low elevation between 0.33 and 1.16 m (Figures 5 and 7). At this site, we consider that the origin of the boulders is associated with the collapse of the landward part of the rocky platform (Figure 11), and do not owe their current location to a high energy event. Therefore, these 12 boulders were not considered in the application of hydrodynamic equations. A particular cluster of very large boulders and an accumulation of smaller ones lies in a heavily eroded area (area 2a), the only part of the study area at a low elevation between 0.33 and 1.16 m (Figures 5 and 7). At this site, we consider that the origin of the boulders is associated with the collapse of the landward part of the rocky platform (Figure 11), and do not owe their current location to a high energy event. Therefore, these 12 boulders were not considered in the application of hydrodynamic equations.
(a) (b) Figure 10. (a) The subaerially eroded red horizon that has acted as a detachment zone is indicated by a yellow arrow; (b) detailed view of this detachment zone.
A particular cluster of very large boulders and an accumulation of smaller ones lies in a heavily eroded area (area 2a), the only part of the study area at a low elevation between 0.33 and 1.16 m (Figures 5 and 7). At this site, we consider that the origin of the boulders is associated with the collapse of the landward part of the rocky platform (Figure 11), and do not owe their current location to a high energy event. Therefore, these 12 boulders were not considered in the application of hydrodynamic equations. Figure 11. The boulders owe their dislocation to the collapse of the landward part of the rocky platform, which is visible in the back. At the western part of the study area (area 3) (Figure 12), the dislocated boulders may have a mixed origin, some dislocated from the submarine environment, supporting a submerged scenario while others may have been part of a roof notch collapsed into the sea that were dislocated on the rocky platform during the high energy event(s). At the western part of the study area (area 3) (Figure 12), the dislocated boulders may have a mixed origin, some dislocated from the submarine environment, supporting a submerged scenario while others may have been part of a roof notch collapsed into the sea that were dislocated on the rocky platform during the high energy event(s).

Storm or Tsunami Origin of Deposits?
It has been noted by several researchers, that the computed volumes may be overestimated, because the shape of the boulders is not a perfect Cuboid [20,24,52]. According to Hoffmeister et al. [53] and Hoffmann et al. [8] boulder volume can be exaggerated even by 30-50%, leading to an

Storm or Tsunami Origin of Deposits?
It has been noted by several researchers, that the computed volumes may be overestimated, because the shape of the boulders is not a perfect Cuboid [20,24,52]. According to Hoffmeister et al. [53] and Hoffmann et al. [8] boulder volume can be exaggerated even by 30-50%, leading to an overestimation by 28-83% on the calculated wave heights. Shah-Hosseini et al. [54] reduced the calculated volume by 25%, in order to reduce potential overestimation to boulders with irregular/prismatic shapes. For the calculation of the boulders volume we considered the potential overestimations based on the boulders shape and irregularities and considered 75% of the total calculated volume.
Two main scenarios were considered, joint-bounded block and submerged boulder, based on the geomorphological evidence of the study area. For the JBB scenario we used the equations of Pignatelli et al. [22] and Nandasena et al. [48], while for the SMB scenario the equations of Barbano et al. [49] and Nandasena et al. [48]. Barbano et al. [49] equations were based on Pignatelli et al. [22] and provided equations for the SMB and SAB scenarios.
According to the results of the JBB scenario for the eastern part of the study area (areas 1 and 2b), a vast majority of the boulders (≈60%) would require a storm wave higher than 5 m to trigger their transportation. In fact, wave height calculations indicate storm waves ranging from 5 m to values higher than 15 m. Such values are significantly higher than those suggested from records of the significant wave height of our study area for the periods 1961-1980 [44], and July 2005-February 2008 [45], where the highest values reaches 4.7-5 m, with a modelled return period of 100 years [44]. On the other hand, in a case of a tsunami, it would require a 4.4 m height in order to flip a 41 m 3 rhombohedric boulder. Thus, the possibility of a tsunami impact is highly possible.
Area 1 hosts many lines of evidence of a tsunami impact. The high-density boulder clusters are oriented perpendicularly (a-axis) to the coast while most are near the detachment zone (see above). Approximately 60% of the boulders of this area would need storm waves greater than 5 m for their dislocation, i.e., values between 5 and 11.2 m, values that far exceed the possible storm wave height of the area.
Area 2b further supports the possibility of a tsunami event. According to the results of JBB scenario, only a tsunami could have dislocated the majority of boulders from their original place and lift them on the coastal slope. In case of a storm, 57.1% of the boulders would require wave height between 5.7 and 17.93 m, for their dislocation, a value incompatible with the available wave records.
On area 3, according to the SMB scenario, 7% of the boulders would require storm wave heights of 5.3-7.40 m to be dislocated to their present location. Taking into account the available data of the significant wave heights of the Levantine Sea [44,45,55], there could be a possibility, although rare for the Mediterranean wave regime, for a storm to move submerged boulders and roll them up to the rocky coast slope. In the case of the JBB scenario 60% of the boulders would require a storm wave height greater than 5m, with values reaching 15 m.

Correlation with Known Tsunami Events
A number of earthquake and tsunami catalogs exist in the literature for the Mediterranean region [25,26,36,40]. According to Fokaefs and Papadopoulos [25], considering tsunamigenic sources within the Levantine sea, the strongest tsunamis to have stricken the coasts of the Levantine Sea are those of 551 AD, 749 AD, 1068 AD, 1202 AD, 1546 AD, and 1759 AD, while on Cyprus they distinguish in 1222 AD a moderate-strong tsunami ( Figure 13). Another known well documented earthquake and tsunami is that of 1303 AD, which was generated near the east coast of Crete and affected a large area of the eastern Mediterranean [56].
According to our radiocarbon results, at least two high energy events appear to have impacted the study area. The age from sample LTL19210A appears to broadly correlate with the 1303 AD tsunami that has displaced at least part of the studied boulders. The 1303 AD event is amongst the largest earthquakes in the Mediterranean (M~8.0), with its epicenter near the eastern coasts of Crete Island [26]. The earthquake brought about severe damage in many locations, such as Crete, Rhodes Island, Palestine, Damascus, Cairo, and was felt even in Constantinople and, possibly, in Venice [36]. The generated tsunami was amongst the largest in the Mediterranean, and struck many coastal areas, such as Heraklion in Crete where " . . . The sea rushed upon the city with such an impetus that many buildings collapsed and many people were drowned . . . " [57], Rhodes, Cyprus, Alexandria in Egypt, Palestine, Lebanon, Antalya, and even reached the Adriatic sea [36,43,58]. Tsunami simulations of the historical tsunami of 1303 AD performed by Yolsal et al. [43] Figure 13. Location of the earthquakes cited in the text (black stars) and sites affected by the 1303 earthquake and tsunami (in red) (from [36,43,58]).

Conclusions
In this paper we analyzed 272 boulders located at Cape Greco, at the southeastern coast of Cyprus. These boulders are located on a rocky platform of Pleistocene calcarenite, at elevations between 0.33 and 4.5 m above sea level. Their size, distribution, and geomorphic characteristics suggests that at least some of the studied boulders are due to a tsunami event. The application of hydrodynamic equations to calculate the minimum storm wave height necessary to dislocate these boulders further supports that some of the boulders were transported by a tsunami event. Radiocarbon dating suggests that at least two high energy events have impacted the study area, one associated with the 1303 AD earthquake and tsunami and the second unrelated to any well-known associated tsunami, based on historical sources. Given the large number and variety of sizes in the studied boulders of Cape Greco, it is highly likely that they owed their dislocation to multiple events from various sources.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1: Boulder measurements and wave height calculations for all the studied boulders.  Acknowledgments: The authors would like to thank Mariolopoulos-Kanaginis Foundation for funding and supporting the activities of this research. We also thank three anonymous reviewers whose comments improved an earlier version of this manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. Figure 13. Location of the earthquakes cited in the text (black stars) and sites affected by the 1303 earthquake and tsunami (in red) (from [36,43,58]).
Another event has been dated at 1512-1824 AD, which does not appear to correlate with a well-known earthquake and associated tsunami. Considering the large number of boulders and sizes in our study area, it is highly likely that they owed their dislocation to multiple events from various sources.

Conclusions
In this paper we analyzed 272 boulders located at Cape Greco, at the southeastern coast of Cyprus. These boulders are located on a rocky platform of Pleistocene calcarenite, at elevations between 0.33 and 4.5 m above sea level. Their size, distribution, and geomorphic characteristics suggests that at least some of the studied boulders are due to a tsunami event. The application of hydrodynamic equations to calculate the minimum storm wave height necessary to dislocate these boulders further supports that some of the boulders were transported by a tsunami event. Radiocarbon dating suggests that at least two high energy events have impacted the study area, one associated with the 1303 AD earthquake and tsunami and the second unrelated to any well-known associated tsunami, based on historical sources. Given the large number and variety of sizes in the studied boulders of Cape Greco, it is highly likely that they owed their dislocation to multiple events from various sources.
Funding: This research was funded by MARIOLOPOULOS-KANAGINIS FOUNDATION.