Presence of boulders associated with an extreme wave event in the western Mediterranean (Cape Cope, Murcia, Spain): possible evidence of a tsunami

The region of Murcia, in the southeast of the Iberian Peninsula, experiences moderate tectonic activity, with earthquakes of up to 6.2–6.9 Mw recorded. Even with seismic activity of this scale there is no geological or instrumental record of tsunamis affecting the area. The presence at Cape Cope, Murcia of a ridge of metre-sized imbricated boulders (comprised of Upper Pleistocene sediments) reaching a height of up to 4 m above sea level, indicates that there has been an extreme wave event in the area during the Holocene. Through studying the wave conditions generated during large storms in this area, the boulder ridges appear to have been caused by extreme waves associated with a seismic event, as a tsunami.


Introduction
The processes resulting in the deposition of coarse-clast ridges in coastal settings has been discussed in recent years. Scheffers et al. (2009) in their studies of storm/tsunami deposits in the Caribbean listed the main diagnostic features discriminating between storm-induced and tsunamigenic coarse deposits and concluded that boulder rampart/ridges formed by tsunami include: boulders of over 300 tons; megaclasts diminishing in size landward; a seaward strip of bare rock; imbrication strongly present; the presence of seaward slopes and smooth landward slopes or boulders being mostly 1 3 angular. Some studies (Noormets et al., 2004;Nott, 2004;Scheffers, 2005) have shown that extreme storm waves are not as efficient as tsunamis in the detachment of, and transport of large boulders. Mastronuzzi et al., (2007aMastronuzzi et al., ( , 2007b argued that the use of boulder accumulations as indicators of tsunamis has been a matter of debate as boulder accumulations can also result from storm wave action. Noormets et al. (2004) indicated that tsunamis, as well as large swell waves, are capable of quarrying large boulders, provided that enough initial fracturing is present.
Spike and Baldburgh (2011) reviewed published examples of boulder transport by tsunamis and concluded that various criteria, for example, the size and mass of boulders, as well as the distance of the boulder deposits to the coast and their position above sea-level, were frequently used to discriminate between a tsunami or storm origin (Bryant & Nott, 2001;Mastronuzzi & Sansò, 2000;Scheffers et al., 2009;Whelan & Kelletat, 2005) and that boulders entrained and transported by tsunamis were supposed to have larger sizes and greater mass, and should be transported further inland than storm boulders (Goto et al., 2009;Whelan & Kelletat, 2005).
Scheffers and Kinis (2014) described multiple imbricated boulder deposits (most of them as ramparts or ridges) with storm or tsunami origins. They completed a comprehensive review of existing literature and concluded that all of these imbricated boulder deposits were associated with tsunamis. They also concluded that generally, strong imbrication, even in large boulders (from 10 to > 200 tons in weight), is best developed in coarse tsunami deposits. Etienne et al. (2011) argued that there are no published accounts of extensive boulder ridge formations by tsunami in any case studies of recent events, but some authors have described boulder trains originated by the 2004 Indian Ocean tsunami (Goto et al., 2007;Nandasena et al., 2011b;Paris et al., 2009). Yamanda et al. (2014) described concentrations of boulders (instead of ridges or alignments) in three major clusters, in relation to the 2011 Tohoku-oki tsunami at Miyako City, Japan. Also, Goto et al. (2012) described boulder deposits that originated from the same tsunami. Recently Lau et al. (2018) identified lines of large coral reef boulders in Fiji, associated with a 1953 tsunami. Spike and Bahlburg (2011) found that if the different waves of a tsunami wave train had sufficient energy to entrain clasts, the transport would occur step by step in a landward direction The process of entrainment, transport, and deposition will be repeated two or three times thus building the groups which then give the appearance of a single transport sequence. The proposed multi-step process model has also been reported in other tsunami boulder transport studies (Nandasena et al., 2011b). Cox et al. (2019) reported wave-tank experiments testing whether storm events could generate imbricated boulder ridges and they found that storm waves can produce all of the features of imbricated boulder deposits and concluded that these types of deposits cannot be used as de facto tsunami indicators. Therefore, the discussion about the origin of boulder ridges due to extreme wave events is still open and a range of approaches and interpretations exist in the literature.
De Martini et al. (2021) founded the Mediterranean area as an ideal open laboratory for paleotsunami research for multiple reasons and conclude that the area appears to be ideal also to face the problematic strom vs. tsunami source differentiation that is very critical in the reconstruction of a paleotsunami history at a specific site. The presence along the Mediterranean coast of large boulders was firstly reported in Greece by Pirazzoli et al (1999). Other studies have described boulders accumulations in several sites of the long coastal perimeter of the Mediterranean (i.e.: 2019a; Alvarez-Gómez et al., 2011;Barbano et al., 2010;Kelletat & Schellmann, 2002;Lario et al., 2017;Maouche et al., 2009;Mastronuzzi & Sansò, 2000, 2004Mastronuzzi & Pignatelli, 2006;Mastronuzzi et al., 2007aMastronuzzi et al., , 2007bRoig-Munar et al., 2018a;Scheffers et al., 2008;Scicchitano et al., 2012;Shah-Hosseini et al., 2013;Vött et al., 2008Vött et al., , 2019. Tsunami generation models show that in the Iberian Peninsula, Murcia is the province affected by the greatest tsunamis, mainly generated by large seismic events in Northern Algeria. The Balearic Islands Ibiza and Minorca also recorded the highest maximum elevations, with heights that may pass the four-meters mark locally. (Alvarez-Gómez et al., 2011).
In the Mediterranean Spanish Coast has been described scarce geological records of EWE (Extreme Wave Events) such as tsunamis or storm surges. Even some historical tsunamis have been reported in the region, their impact at the coast was negligible (De Martini et al., 2021). During the past years some boulders deposits, usually with ridge morphology, have been studied in the area. Roig-Munar et al. (2019a) founded large boulders on marine cliffs in more than 20 study sites on Majorca Island Deposits consist of large imbricated boulders of up to 25 t located on platforms that form the rocky coastline of the island, several tens of meters from the edge of the cliff, up to 15 m above sea level. They are mainly located on the eastern and southern coasts of Majorca, where wave height and energy are low compared with those from the N and NE. The study conclude that a tsunami appears to be a more likely mechanism for transportation and deposition of the boulder field than storm waves. Along the rocky coastline of Minorca Island circa 25 sites with large imbricated boulders deposits have been found on platforms located several tens of meters from the edge of the cliff, up to 15 m above the sea level (Gomez-Pujol & Roig-Munar, 2013;Martín-Prieto et al., 2019;Roig-Munar et al., 2018a). Roig-Munar et al. (2018a) analysed 3000 large boulders along the studied 1 3 sites and concluded that they should be dislodged and positioned by the action of tsunami waves, although some of these boulders were also displaced by storm waves. Also, in seven study sites on Ibiza and Formentera Islands large boulder deposits were located several tens of meters from the edge of the cliff, up to 11 m above sea level and several kilometers away from any inland escarpment (Roig-Munar et al., 2019b). Six sites were identified in North Castellon, Valencia Region with boulders ridges with sedimentary characteristics typical of tsunamis flows (Roig-Munar et al., 2018b). In 1522 there was an earthquake M > 6.5 in Almería that affected large areas of the western Mediterranean. The epicentral area is located offshore on the platform of the Gulf of Almería. The earthquake triggered several underwater landslides that could have produced a tsunami (Reicherter & Becker-Heidmann, 2009). Tomassetti et al. (2021) described in Malaga a tsunami deposit in an archaeological context assigning it to an earthquake documented in historical sources in the year 881 AD.
The study of the boulder ridge located at Cope Cape complements the studies carried out so far and confirms the occurrence of tsunamis in eastern Spanish coast.

Geodynamical framework
The structural framework of the eastern Betic Cordilleras (SE Spain) is linked to the occurrence of a large left-lateral transcurrent zone known as the Eastern Betics Shear Zone (EBSZ) (Larouzière et al., 1988;Silva et al., 1993) that evolved since the middle Miocene because of the continuous northward indentation of the Aguilas Arc block (Coppier et al., 1989). The Aguilas Arc is a differentially up-thrusted cortical structure characterized by a large-scale bending driven by two main fault systems-N10°-20° E system (Palomares-Los Arejos faults) and N90°-100° E (El Saladillo-Las Moreras faults) that separates the arc north Africa (IGN). White circles-historical earthquakes. Depth of instrumental earthquakes-red circles (0-30 km); yellow circles (30-60 km); green circles (60-180 km); blue circles (≥ 600 km). Diameter of circles is related to the magnitude of recorded events, with maximum instrumental magnitude never exceeding M7.5 itself from the Guadalentín Tectonic Depression, what in fact is the central segment of the EBSZ (Silva et al., 1993) (Fig. 1). Repeated rotation of the regional stress field from N-S to NNW-SSE since Tortonian times (Coppier et al., 1989;Grievaud, 1989;Montenat et al., 1987;Ott d'Estevou & Montenat, 1985;Silva et al., 1993) promoted the opening of different sedimentary basins within the Aguilas Arc ( Fig. 1), that have evolved differentially depending on their location within the arc and on the direction and kinematics of their inner faults (Bardají, 1999). Cope Basin has been defined as a littoral detachment basin (Bardají, 1999) given its extensional behaviour within the general compressive tectonic framework. N120° E and N60° E faults drove its opening during the early Pliocene and later evolution of this basin (Fig. 1A).
The Eastern Betics Shear Zone is one of the most seismically active areas of the Iberian Peninsula (Fig. 1b), however most of this activity is concentrated in the periphery of the Aguilas Arc, along those fault systems oriented perpendicularly or slightly oblique to the main stress field. The strongest earthquakes, reaching magnitude 5.2 Mw (instrumental) or 6.2 to 6.9 Mw (estimated magnitude for historical earthquakes; Alfaro et al., 2012;Martínez Díaz et al., 2011) are associated to these ENE-WSW to E-W faults (e.g. Lorca-Alhama and Bajo Segura faults) that outline the periphery of the Arc. Other active faults that run in parallel to the main N-S trending stress field, such as the Palomares fault system (N10°-20° E) yield very high slip rates (30 km lateral offset since the early Pleistocene; Bousquet, 1979;Weijermars, 1987) but very low seismic activity (≤ 3Mw) even along its N120° E prolongation (Saladillo fault). In the Aguilas Arc itself, the main master faults show high displacement rates but scant seismic activity (minimal instrumental earthquakes ≤ 4Mw; and only two historic earthquakes reported in Aguilas (AD1596 and AD1882) with maximum MSK intensity of IV-V (Mezcua, 1982). Liquefaction structures and shoreline vertical displacements have been identified in the Cope Basin (Bardají et al., 2015) and related to seismic activity during MIS 5e-Holocene. Both features point to a maximum seismic intensity between VII and IX, after the INQUA ESI-07 macroseismic scale (Environmental Seismic Intensity Scale; Guerreri & Vittori, 2007).
Regarding to submarine faults that could act as a potential tsunamigenic source, no one has been reported close to the studied area but mainly in the Alboran Sea where the 2016 earthquake (Mw 6.4) is the largest ever recorded event. Main active faults capable of generate tsunamis in the southern Iberian coasts are the Al-Idrissi, the Serrata-Carboneras and the Yusuf-Habibas fault zones (Gràcia et al., 2019;Somoza et al., 2021;Vázquez et al., 2021, and references there in).

Sedimentary filling of the Aguilas basin
The studied boulder deposits are located in the Aguilas basin. The earliest sedimentary filling of this basin was strongly dominated by Plio-Pleistocene fossiliferous yellow marine calcarenites, the extension and outcrop of which are controlled by the above mentioned N60° E and N120° E faults (Fig. 2). This unit is unconformably overlain by a well-cemented, sea-level-controlled, off-lapping Pleistocene sedimentary sequence, gently descending to the sea, where coastal gravel sediments alternate with subaerial alluvial fans (Bardají, 1999;Bardají et al., 1986Bardají et al., , 2015Dabrio et al., 1991;García-Tortosa et al., 2004;Zazo et al., 2013). The whole sequence is composed of up to ten marine units, with the three more recent ones presenting a staircase arrangement within the general off-lapping trend. These three younger units have been attributed to MIS7, MIS 5e and MIS 5 c/a based on the faunal content (Strombus bubonius) and specific lithology (oolithic dune-beach associated to the intermediate unit) (Bardají et al., 2009(Bardají et al., , 2015. The younger of these three units (quartzous beach-dune system) is scarcely represented along the basin but the two older ones are widely developed along the coast usually recorded on top of a wave-cut platform carved into the lower yellow calcarenites (Fig. 3).
MIS7 unit (B in Fig. 3) is characterised by fine bedded (decimetre scale) calcarenites alternating with gravel layers, deposited in a foreshore-shoreface environment. MIS5 unit (D in Fig. 3) is a grey conglomerate of coastal origin, characterized by abundant white rounded quartz pebbles, cobbles and some boulders. Occasionally along the shore, reddish terrestrial sediments separate this unit form the older one (Unit C in Fig. 3) (Bardají et al., 2015).
In the studied zone (shaded in grey in Fig. 2), these units are eroded by the present-day wave-cut platform where the extreme wave event (EWE) displaced boulders are located, the provenance and source lithology of which are the MIS7 and MIS5 units (described above) as well as the underlying yellow calcarenites.

Meteorological conditions
Due to its latitude and location in the southeast of the Iberian Peninsula, the location is affected by very diverse atmospheric conditions, associated with the interaction of cold polar and warm tropical air masses. During much of the year it is subjected to the action of aerological shelter exerted by the Azores anticyclone that is responsible for long and intense droughts. However, the characteristics of the weather and, above all, the production of rain, depends more closely on polar mechanisms than on tropical ones (Capel Molina, 1991). Although the territory of Murcia is outside the area of greatest turbulence and cyclogenesis, located to the north of the 40º parallel, it remains under the influence of meridian circulation flows from the east, isolated cold irruptions in height (vortices detached from the jetstream) and thermoconvective processes developed in the western Mediterranean that trigger significant torrential rain events (Conesa & Alonso, 2006).
The presence of storms in the Mediterranean associated with convection phenomena gives rise to a regime of easterly winds (levantes), typically humid on the eastern flank of the region and drying out as they move inland. This situation predominates during spring and summer, extending until autumn. The frequency distribution of the wind directions in the Murcia region varies according to the geographical  (Bardají, 1999) The 'studied zone' is shaded in grey, and white star in red circle mark the location of boulder ridge outcrop studied in detail in this work (Figs. 4 and 5) 1 3 situation and the topographic relief of the area. Although the winds of the first quadrant dominate in a large part of the region, it is in the eastern coastal zone where they are more persistent, especially in summer, with observed frequencies of 26.5% for the E direction, 29.6% for the winds from the NE and 16% for those from the SE. In winter, on the other hand, winds from the north are more frequent: NW 19.2% and N 15.8% (Conesa & Alonso, 2006).
Since there are no oceanographic parameter buoys around Cape Cope, SIMAR (operational wave prediction system) data was used to establish the maritime climate (www. puert osdel estado. es). The SIMAR dataset consists of time series data (in this case covering the period 1958-2021) of wind and wave parameters from numerical modelling. The source is therefore synthetic data and does not come from direct measurements. This study utilised SIMAR point 2,067,089, the nearest to the study area (5.5 km east of Cape Cope). Maximum monthly significant wave heights (Hs) at this point during 1958-2022 period have been calculated as 2.23-4.27 m at this point. At the SIMAR point the significant swell recorded resulted mainly from the E (50% of the time) and from the SSW (28% of the time). For a real maximum wave height that has been recorded from direct measurement, data from the nearest sea buoy, Cabo de Palos (code 1613), located 75 km NE of Cabo Cope, has been utilised. During the period 1985-2012 the measured monthly maximum heights varied between 5.66 m and 14.8 m, and the significant heights (Hs) between 2.94 m and 5.64 m.

Boulder accumulation of Cape Cope
Deposits consisting of heterometric boulders were recorded along the coast to the north of Cape Cope (Fig. 2). The boulders formed a ridge parallel to the coastline, but also isolated boulders were present in the study area. The main boulder ridge reached heights of 4.0 m asl and were composed largely of angular boulders, but also with some sub-rounded and rounded boulders present. Some boulders appeared to be imbricated and some individual boulders had long axes of > 4.15 m, volumes of up to 7.1 m 3 and weights of 17.7 tonnes (Fig. 4).
The boulders under investigation in this study were made up of three different lithologies (Bardají et al., 2015): (1) Yellow fossiliferous calcarenites of medium-coarse grain size and massive character (Upper Pliocene-Lower Pleistocene); (2) Calcarenites of medium-fine grain size alternating with layers of gravels of variable thickness, and layers with frequent bivalve molds (Upper Pleistocene); (3) Highly cemented grey conglomerates in which two facies have been differentiated, a conglomerate of very rounded gravels (1-2 cm) mainly of white quartz, and a heterometric conglomerate with very rounded gravels, pebbles, and boulders of different lithologies (Upper Pleistocene). Three younger units present a staircase arrangement in the general offlapping sequence with the youngest further staircased into the two older units. The youngest was only barely visible in a single outcrop at 0.5 m asl but the two older ones were easy to trace along the basin and have been attributed to MIS 7 and MIS 5 (Bardají et al., 2009(Bardají et al., , 2015. The boulder deposits appeared to originate from joint bounded blocks (JBB) originally forming part of a wave-cut platform carved into the Plio-Pleistocene and the Upper Pleistocene marine terraces.
For the geolocation of the boulder ridges, several flights were made with a DJI Phantom 4 unmanned aerial vehicle (drone), that carried out a photogrammetric flight and subsequent 3D modelling utilising Pix4D software. To obtain a centimetric resolution of the geolocation of the boulders and their dimensions, various control points were obtained with a Trimble Rover R8 GNSS GPS system. To obtain millimetric detail, terrestrial photogrammetry of the boulders was carried out using the commercial software Agisoft Metashape. Three-dimensional restitution and modelling of the boulders and integrating the resulting products of both the photogrammetric flight and the terrestrial photogrammetry was completed using GIS facilitated through the free software QGIS (Fig. 5). In this way, the morphometry, orientation, altitude, and the distance from the coastline of each boulder was determined.
To calculate the densities of the different lithologies identified, a minimum of 2-3 samples of each was taken. In the laboratory, each was submerged in water for at least 24 h to achieve total saturation to accurately simulate the submerged environment. Once saturated each sample was weighed (m i ) and the volume (v i ) of each sample was calculated by immersing them in water and measuring the volume of displacement. All measurements were made twice to minimize errors. Finally, the density (ρi) was calculated using the basic formula ρi = mi/vi. Average density has been estimated at 2.5 gr/cm 3 .
Of the outcrops mentioned in Fig. 2, the central outcrop was selected to measure the boulder parameters (south Rambla Elena, star in Fig. 2). A total of 500 boulders were measured, for each boulder the following parameters were recorded; the major axis (a), the middle axis (b) and the minor axis (c), the distance to the coastline, measured to the high-water mark (d), the height above sea level (h), and the angle of the slope of the surface on which it is located (Table 1). Usually, to calculate volume the formula V = a · b · c is used. Several authors used automatic volume calculation from scanned boulder images and concluded that scan volume variates between 40-60% of the volume calculated by axis estimation (Engel & May, 2012;Lario et al., 2017Lario et al., , 2020. Thus, a correction factor of 60% was applied to the volume calculated by using the V = a · b · c formula. 48 boulders were selected as representative, including the heaviest and the lightest, the furthest from the coast and the closest, as well as the highest and the closest to sea level, plus a set of randomly selected boulders (Table 1).
The average distance of the boulder deposits to the coast was 19.4 m (ranging from 0.2-39.5 m), the average altitude was 2.5 m asl (with a range of 0.5-4.0 m asl), average corrected volume was 1.1 m 3 (ranging from 0.04-7.1 m 3 ) and the average weight was 2796 kg (range being 95.1-17,754 kg). A first approximation demonstrates that, as would be expected, heavier boulders were distributed nearest to the present coast while the lighter boulders were deposited further inland (Fig. 6).
Likewise altitudinal data showed a similar relationship with heavier boulders present at lower heights above sea level than the lighter ones. (Fig. 7).  1 3

Methodology for modelling boulder transport
Several studies have tried to infer the height of waves required to transport boulders at the coast. Nott (2003) presented formulas used to calculate wave heights required to move boulders by both tsunami and storm waves in three different pre-transport situations: submerged boulders; subaerial boulders; and joint-bounded blocks (JBB). To incorporate the fact that a boulder has one side facing the wave, one top surface, and is limited by four sides. Later Pignatielli et al. (2009) proposed some modifications to Nott´s equations. Also, Barbano et al. (2010) modified these equations after testing wave transport formulas on the coast of Sicily. Nandasena et al. (2011a) improved Nott´s equations further and found that the minimum flow velocity, derived from their revised equation, required to initiate the transport of submerged boulders was less than that inferred from Nott's equations. Assuming wave heights from the initiation-of-motion approach is sometimes considered questionable, although it has been used in previous work. Mastronuzzi et al., (2007a, 2007b) used Nott's (2003 equations for hydrodynamic calculations relating to boulder accumulations in south-eastern Salento (Italy). This methodology has been used to study the movement of boulders related to extreme wave events (Spike et al., 2008, in the Caribbean;Maouche et al., 2009, in  Algeria; Switzer & Burston, 2010, in southeast Australian;Costa et al., 2011, in southern Portugal;Spike & Baldburgh, 2011, in Chile;Engel & May, 2012, in the Caribbean;Lario et al., 2017 in SW Spain andLario et al., 2020, in Mexico).
Taking all of these formulas into account, and considering that most of the boulders examined in this study corresponded to a JBB scenario, the equations proposed by Engel and May (2012), that used a modification of Nott's equations proposed by Nandasena et al. (2011a) for the JBB scenario have been applied, where: where Ht = tsunami wave height; Hs = storm wave height; C L = coefficient of lift (0.178); a = A axis of boulder; ρw = density of sea water = 1.02 g/ml; ρb = boulder density (in this study an average density of 2.5 g/cm 3 was used); V = corrected boulder volume; q = boulder area coefficient (0.73); b = B axis of boulder; µ = coefficient of static friction (0.65); θ = bed slope angle of the pre-transport setting.
If the initiation of movement of most of the boulders was in subaerial conditions, the Barbano et al. (2010) modification of Nott´s equations, the following formulas, was applied: where C d the coefficient of drag = 2; C m the coefficient of mass = 1; ü = 1 m/s 2 is the flow acceleration; g is the acceleration due to gravity (9.81 m/s 2 ). Nandasena et al. (2011a) developed a "boulder transport histogram" to represent the range of flow velocities that satisfied the requirements for initial transport of a boulder via different modes of transport: sliding, rolling, saltation or lifting. The boulder transport histogram can be used to predict the possible initial transport modes of a dislodged boulder from the flow velocity.
In the case of a submerged or subaerial boulder the following scenarios exist: Transport initiates with sliding.
Transport initiates with saltation.
Or in the case of a joint-bounded boulder where transport only begins with saltation/lifting: Wave height necessary to move the boulders by storms (Hs) or by tsunamis (Ht) in a subaerial scenario, applying the Barbano (2010) equations In Eqs. (5) to (8), u is the flow velocity, g is gravitational acceleration (9.81 m/s 2 ), µ s is the coefficient of static friction between the boulder and the bed (0.7) and C d is the coefficient of drag (1.95) (Nandasena et al., 2011a(Nandasena et al., , 2014.
A relationship between flow depth and flow velocity described by the Froude number as Fr = u∕ √ gh , where h is the flow depth, was used to estimate the flow departing depth from the minimum flow velocity required to move a boulder (Nandasena et al., 2014). The Froude number for past tsunamis has been determined to be between 0.7 and 2.0; in this study values of 1.0 and 1.5 were used, as proposed by Nandasena et al., (2012Nandasena et al., ( , 2014.

Results
The results obtained by applying formulas (3) and (4) to the boulders of the Cabo Cope site are presented in Table 1. They show that to move the heaviest and largest boulders in subaerial conditions, tsunami wave height would need to exceed 1.75 to 2.45 m (based on mass of the five largest boulders found in different locations in the study area). The height of storm waves capable of moving those boulders would need to be at least 6.90 to 9.85 m (Fig. 8).
The results obtained using Engel and May (2012) formula for JBB scenarios (Eqs. (1) and (2)) indicated that to move the heaviest and largest boulders, tsunami wave heights would need to have exceeded 2.7 to 3.1 m (based on the mass of the five largest boulders found in different locations in the study area). The height of storm waves capable of moving the same boulders would have needed to be at least 10.9 to 12.3 m. (Fig. 9).
When considering the flow velocity, the results obtained by applying Eqs. (5) to (8) to the boulders investigated in this study (considering weight, distance to the coast and altitude) are presented in Table 2 and Fig. 10 which shows the corresponding boulder transport histogram. These figures also display the minimum flow velocity and range of flow velocities that would have been required to initiate boulder transport at each site. In the case of submerged/ subaerial boulders, the average results indicate that a wave celerity of not less than 3.65 m/s would have been required to transport boulders from their pre-transport location. If the average flow velocity was greater than 3.6 m/s, the boulders would have been transported by sliding. If the average flow velocity was > 5.7 m/s, the boulders would have been transported by rolling; if average flow velocity was > 7.8 m/s, boulder transport would have been via saltation.
In the case of the original position of the boulders being a JBB scenario (Fig. 11), the boulder would remain stable until the average flow velocity for transport by saltation/ lifting is higher than 8.0 m/s. Bosnic et al. (2021) modelled the 1755 Lisbon tsunami in the Algarve to assess tsunami onshore flow characteristics. Results from the inverse model demonstrated tsunami onshore average flow velocities ranged from 7.3 up to 9.3 m/s and forward modelling results show tsunami onshore velocities of 7 m/s. These values are similar to the 5 m/s estimation derived from boulder deposits by the 1755 CE Lisbon tsunami in Costa et al. (2011) and also are comparable to the data of Bujan and Cox (2020) where the average onshore velocity for a 3 m offshore tsunami wave was estimated to be around 5 m/s. Consequently, in this study, it appears that the flow velocities calculated as necessary to move the boulders at Cape Cope are similar to the results of previous tsunami studies.   Table 2 presents flow depths necessary to initiate boulder transport using the formulas applied in this paper. Nandasena et al. (2014) indicated the difficulty in establishing the Froude number at each site. Using a Fr = 1.0, the minimum flow depth necessary to allow for: (i) sliding averages at 1.4 m; (ii) rolling has an average of 3.5 m); and (iii) for saltation averages at 6.5 m). Lifting, as is the case of JBB, required flow average of 7.0 m. When using Fr = 1.5 the

Discussion: extreme event origin and age
The presence of imbricated boulder ridges on the coast of Cape Cope is related to an extreme wave event, either a tsunami or a strong storm. Although there are various formulae relating the deep-water significant wave height to the runup of a wave breaking on the shore that consider seabed topography, wave frequency, or wave amplitude, some approximations from direct measurements suggest a relationship between runup and the significant wave measured offshore (Guza & Thorton, 1982). Senechal et al. (2011) calculated runup from empirical data observed in extreme storm conditions and found that the total vertical runup elevation is significantly correlated to H0 with this relationship R = 2.14*tanh (0.4H 0 ). Thus, the maximum offshore wave heights predicted by the SIMAR network in the area, corresponding to 4.27 m (maximum significant wave height at SIMAR point 2,067,089) and 5.64 m and 14.8 m (maximum significant wave height and maximum wave height at Cabo de Palos buoy). Such waves would have a maximum coastal runup of 2.00 to 2.14 m during extreme storms, far away from the necessary storm waves calculated to move the boulders. It is a fact that the SIMAR record corresponds only to 1951-2021 period, and there may be previous episodes of storms larger than those recorded, although they do not seem to have been observed in the geological record. The necessary waves to move the studied boulders during storms (Hs) has been calculated to be between 10.9 to 12.3 m in JBB scenario and between 6.9 to 9.85 m in subaerial setting, therefore, attributing the transport of the boulders to storm events is ruled out utilising the historically recorded extreme sea conditions in this region. Therefore, it is necessary to try to determine whether a possible tsunami could have affected the coast and thus be responsible for this deposit. In relation to tsunami generation, the decision matrix for tsunamis on the Spanish Mediterranean coast indicates that earthquakes from Mw 6.0-6.5 generated less than 40 km from the coast and at depths of less than 100 km have the potential to generate a destructive tsunami on this coast (IOC, 2011). Tsunami generation models also show that the Murcia region is the area most affected by large tsunamis generated mainly in the north of Algeria . Even so, references to tsunamis occurring in the area in historical times are very scarce. Bardají et al. (2015) reported a tsunami event in Cape Cope occurring during the MIS 5e. In North Castellon (Valencia Region, Spain) six sites were identified with boulder ridges present that have sedimentary characteristics typical of tsunami flows (Roig-Munar et al., 2018, 2019a, 2019b. The boulders had an average weight of 1.5 t and were located at an average distance from the coastline of 15.5 m and at 2.3 m asl. Using the same approaches as used in previous studies the results indicated that local storms were not capable of moving the boulders to their present situation and thus a tsunami wave would appear to be the main cause of transport and deposition. Also, the average boulder orientation indicated that this coastline would be affected by tsunami waves generated in north Algeria passing through the Ibiza-Majorca channel, as has been demonstrated by the models of Alvarez-Gomez et al. (2011). In tsunami catalogues (Lario et al., 2011) there are no data available that provided an age of the events.
Data from Estepona (south Malaga) demonstrates the presence of a sedimentary layer rich in marine and continental materials located at 2.0 m asl, that appears to be representative of a tsunami backwash deposit that is very likely to correspond to the tsunami reported in historical documents that occurred in the 881 AD (Tomasseti et al., 2021).
In 1522, an earthquake of M > 6.5 in Almeria affected large parts of the western Mediterranean, causing more than 1000 deaths. The epicentral area was located on the Gulf of Almeria shelf and produced several submarine landslides that could have produced tsunami (Reicherter & Becker-Heidmann, 2009). These authors found sediments associated with this (and previous) events that indicated wave heights of at least 2-3 m, although evidence of this tsunami has not been found north of Cabo de Gata.
With this context of other tsunami activity in this region the EWE explored in this study can be assigned to a tsunami generated in the western Mediterranean and which may have reached at least 4 m asl on the coast of Cape Cope. Regarding the age of the event, it has not been possible to date the boulders studied but given their disposition with respect to the underlying deposits, their age must be Holocene post-maximum transgressive (ca. 6000 yrBP). The historical earthquakes and tsunamis studied in the area show two possible tsunamis (at 881 AD and 1522 AD) that could be responsible for the boulder ridges present in Cape Cope, although the wave heights referred to in these historical events are far from those necessary to move the boulders. Even in some cases have been possible to suggest a Tsunami Intensity value regard with the environmental impact of the tsunami (Lario et al., 2016), the parameters of the Cope Cape deposits does not allow assigning a scale to the event.

Conclusions
The characteristics of the deposits found on the coast of Cape Cope conform to the features of an EWE. The application of different equations for the study of wave heights, flow velocity and depth necessary to mobilize these boulders indicate that the extreme storm conditions recorded in the area could not form these boulder ridges, assigning their formation to a high-energy tsunami-type event. These types of deposits are not unique and their presence in other areas of eastern Spain indicates that the occurrence of tsunamis in the area are not isolated events but have occurred on several occasions at least during the Holocene. Once again, the geological record of these events provides data for the study of the seismic hazards on the Mediterranean coast.