Search

This review first pays tribute to the famous Croatian oceanographer, Mira Zore-Armanda, and her seminal work on the Adriatic water masses in 1963, and emphasises the importance of the densest Mediterranean water mass: North Adriatic Dense Water (NAddW). This water mass is generated through substantial wintertime surface cooling and evaporation over the wide northern Adriatic and is known to (1) influence the Adriatic-Ionian thermohaline circulation, (2) bring oxygen and carbon to the deep Adriatic layers and, (3) more generally, have a substantial impact on the physics and biogeochemistry of the whole Adriatic. Second, the NAddW physics, from preconditioning, through generation and spreading, to accumulation in Adriatic depressions, is reviewed. Then, the temporal evolution of the NAddW properties influenced and connected to (1) basin-wide interannual and decadal variability and (2) trends towards warmer and saltier source characteristic due to ongoing climate change, is discussed. The importance of long-term observations and atmosphere-ocean modelling in event, decadal and climate studies is then presented. Finally, a review of the identified gaps and perspectives for future research is concluding this article. [ABSTRACT FROM AUTHOR]

Four chronologies of the bivalve species Glycymeris pilosa have been constructed along a 300 km gradient of the eastern coastal Adriatic Sea, all of which span the common period of 1982–2015. The chronologies are compared to local and remote environmental drivers suspected to influence the biology of the system, including air and seawater temperature, precipitation and freshwater discharge. The Adriatic-Ionian Bimodal Oscillating System (BiOS), a key oceanographic feature quantified by satellite-derived absolute dynamic topography, is also compared to the chronologies. The chronologies at the two southern sites are more strongly influenced by local river discharge, while the two northern chronologies are more strongly influenced by BiOS. These results highlight the broadscale importance of BiOS to the Adriatic system as well as the heterogeneity of nearshore environmental and drivers of growth. These G. pilosa chronologies provide unique multidecadal, continuous, biological time series to better understand the ecology and fine-scale variability of the Adriatic with potential for other shallow, semi-enclosed seas. [ABSTRACT FROM AUTHOR]

This paper aims to quantify the productivity of research concerning the Adriatic Sea, with a focus on oceanography and meteorology. Productivity and impact were measured by analysing articles and citations from the Thomson Reuters Web of Science database, spanning the period 1994-2008. The most productive country was Italy but the highest number of citations was achieved by articles from Germany (all Adriatic publications) and Spain ("Oceanography" and "Meteorology and atmospheric sciences", only). By contrast, the second-most productive country, Croatia, had the lowest citation rate. Collaborations between Adriatic researchers were driven not only by the geographical position of a country (e.g., Italy vs. Croatia), but also by investment rates in Adriatic research (e.g., Italy vs. USA and Croatia vs. USA). Such collaborations substantially improved the impact of the research, especially from transitional countries such as Croatia. [ABSTRACT FROM AUTHOR]

This paper analyzes potential recruitment effect of endangered bivalve Pinna nobilis in a small marine protected area-Little Lake (Adriatic Sea). A box model has been applied for the estimation of residence time of the lake, forced by tides only at the open end of the lake system. The residence time of Little Lake was estimated to be 54 days. In our approach we assumed that P. nobilis larvae may be found in the surface waters only during the night, thus restricting their outgoing from the lake to about 7-8 hours per day in the June-August period, as the larvae need some time to reach surface waters. Therefore, the equivalent residence time of the larvae in Little Lake is estimated to be about 160-185 days. Knowing that the bivalve larva phase lasts between 3 and 5 weeks, one may estimate that only 11 to 22% of the larvae, respectively, leave Little Lake. From a conservation standpoint these low values imply a special level of protection within Mljet National Park In addition, future climatic changes will affect the recruitment of P. nobilis, not only due to sea level rise but also with temperature and other changes which may influence its reproduction phase and larvae transport. [ABSTRACT FROM AUTHOR]

In this study, the Adriatic Sea and Coast (AdriSC) kilometre-scale atmosphere–ocean climate model covering the Adriatic Sea and northern Ionian Sea is presented. The AdriSC ocean results of a 31-year-long (i.e. 1987–2017) climate simulation, derived with the Regional Ocean Modeling System (ROMS) 3 km and 1 km models, are evaluated with respect to a comprehensive collection of remote sensing and in situ observational data. In general, it is found that the AdriSC model is capable of reproducing the observed sea surface properties, daily temperatures and salinities, and the hourly ocean currents with good accuracy. In particular, the AdriSC ROMS 3 km model demonstrates skill in reproducing the main variabilities of the sea surface height and the sea surface temperature, despite a persistent negative bias within the Adriatic Sea. Furthermore, the AdriSC ROMS 1 km model is found to be more capable of reproducing the observed thermohaline and dynamical properties than the AdriSC ROMS 3 km model. For the temperature and salinity, better results are obtained in the deeper parts than in the shallow shelf and coastal parts, particularly for the surface layer of the Adriatic Sea. The AdriSC ROMS 1 km model is also found to perform well in reproducing the seasonal thermohaline properties of the water masses over the entire Adriatic–Ionian domain. The evaluation of the modelled ocean currents revealed better results at locations along the eastern coast and especially the northeastern shelf than in the middle eastern coastal area and the deepest part of the Adriatic Sea. Finally, the AdriSC climate component is found to be a more suitable modelling framework to study the dense water formation and long-term thermohaline circulation of the Adriatic–Ionian basin than the available Mediterranean regional climate models. [ABSTRACT FROM AUTHOR]

In this study, the Adriatic Sea and Coast (AdriSC) kilometre-scale atmosphere-ocean climate model covering the Adriatic and northern Ionian Seas is presented. The AdriSC ocean results of a 31-year long (i.e. 1987-2017) climate simulation, derived with the Regional Ocean Modeling System (ROMS) 3-km and 1-km models, are evaluated with respect to a comprehensive collection of remote-sensing and in situ observational data. In general, it is found that the AdriSC model is capable to reproduce the observed sea-surface properties, daily temperatures and salinities and the hourly ocean currents with good accuracy. In particular, the AdriSC ROMS 3-km model demonstrates skill in reproducing the main variabilities of the sea-surface height as well as the sea-surface temperature, despite a persistent negative bias within the Adriatic Sea. Furthermore, the AdriSC ROMS 1-km model is found to be more capable to reproduce the observed thermohaline and dynamical properties than the AdriSC ROMS 3-km model. For the temperature and salinity, better results are obtained in the deeper parts than in the shallow shelf and coastal parts, particularly for the surface layer of the Adriatic Sea. The AdriSC ROMS 1-km model is also found to perform well in reproducing the seasonal thermohaline properties of the water masses over the entire Adriatic-Ionian domain. The evaluation of the modelled ocean currents revealed better results at locations along the eastern coast and especially the north-eastern shelf than in the middle-eastern coastal area and the deepest part of the Adriatic Sea. Finally, the AdriSC climate component is found to be a more suitable modelling framework to study the dense water formation and long-term thermohaline circulation of the Adriatic-Ionian basin than the available Mediterranean regional climate models. [ABSTRACT FROM AUTHOR]

A climate synoptic meteotsunami index has been constructed using synoptic variables for the Balearic Islands. The index allows for the very first assessment of atmospherically driven intense sea level oscillations at the tsunami timescale (<2 hr) in future climates. The index has been computed using outputs from evaluation, historical, and three scenario MED‐11_CNRM Med‐CORDEX regional atmospheric climate runs. The reliability of the index has been verified against reanalysis simulations and on documented meteotsunami events. No significant changes in the index are projected under RCP2.6 and RCP4.5 scenarios, while the yearly number of days with meteotsunamis is expected to increase by 34% under the RCP8.5 scenario by the year 2100. This increase will dominantly occur during the summer season (May–August), being contemporaneous with maximum index values. The presented results are relevant for assessment of sea level extremes worldwide, since high‐frequency sea level oscillations may contribute up to 40% of the total range. Plain Language Summary: The study aims to develop the proxy methodology to be used for assessment of meteotsunamis—atmospherically driven waves in a tsunami frequency band—in the future climate, as the present climate models are far from reproduction of processes on a minute and a kilometer scale. The proxy estimates are based on connecting synoptic parameters and patterns to the high‐frequency sea level observations. The methodology is successfully verified for Ciutadella harbor, a known hot spot for meteotsunamis, on state‐of‐the‐art reanalysis data. It seems that meteotsunamis are likely not to increase in all but business‐as‐usual (RCP8.5) climate scenarios, where the increase is matching the season when meteotsunami reach their maximum (spring–summer). The methodology might be useful for assessing meteotsunami hazard in future climates, while being applied at places where meteotsunamis are driven dominantly by a single atmospheric process. Global and regional analyses indicate that these might encompass the majority of midlatitudal regions. Key Points: A proxy‐based index has been developed for assessing meteotsunamigenic potential in future climateNo significant changes in the index have been found for the RCP2.6 and RCP4.5 scenariosThe number of days with meteotsunamis is expected to increase by 34% by the end of the century for the RCP8.5 scenario [ABSTRACT FROM AUTHOR]

Present investigations of sea level extremes are based on hourly data measured at coastal tide gauges. The use of hourly data restricts existing global and regional analyses to periods larger than 2 h. However, a number of processes occur at minute timescales, of which the most ruinous are tsunamis. Meteotsunamis, hazardous nonseismic waves that occur at tsunami timescales over limited regions, may also locally dominate sea level extremes. Here, we show that nonseismic sea level oscillations at tsunami timescales (<2 h) may substantially contribute to global sea level extremes, up to 50% in low-tidal basins. The intensity of these oscillations is zonally correlated with mid-tropospheric winds at the 99% significance level, with the variance doubling from the tropics and subtropics to the mid-latitudes. Specific atmospheric patterns are found during strong events at selected locations in the World Ocean, indicating a globally predominant generation mechanism. Our analysis suggests that these oscillations should be considered in sea level hazard assessment studies. Establishing a strong correlation between nonseismic sea level oscillations at tsunami timescales and atmospheric synoptic patterns would allow for forecasting of nonseismic sea level oscillations for operational use, as well as hindcasting and projection of their effects under past, present and future climates. [ABSTRACT FROM AUTHOR]

A series of tsunami-like waves of non-seismic origin struck several southern European countries during the period of 23 to 27 June 2014. The event caused considerable damage from Spain to Ukraine. Here, we show that these waves were long-period ocean oscillations known as meteorological tsunamis which are generated by intense small-scale air pressure disturbances. An unique atmospheric synoptic pattern was tracked propagating eastward over the Mediterranean and the Black seas in synchrony with onset times of observed tsunami waves. This pattern favoured generation and propagation of atmospheric gravity waves that induced pronounced tsunami-like waves through the Proudman resonance mechanism. This is the first documented case of a chain of destructive meteorological tsunamis occurring over a distance of thousands of kilometres. Our findings further demonstrate that these events represent potentially dangerous regional phenomena and should be included in tsunami warning systems. [ABSTRACT FROM AUTHOR]

Worldwide tsunamis driven by atmospheric waves—or planetary meteotsunami waves—are extremely rare events. They mostly occur during supervolcano explosions or asteroid impacts capable to generate atmospheric acoustic-gravity waves including the Lamb waves that can circle the globe multiple times. Recently, such ocean waves have been globally recorded after the Hunga Tonga–Hunga Ha'apai volcano eruption on 15 January 2022, but did not pose any serious danger to the coastal communities. However, this study highlights that the mostly ignored destructive potential of planetary meteotsunami waves can be compared to the well-studied tsunami hazards. In practice, several process-oriented numerical experiments are designed to force a global ocean model with the realistic atmospheric response to the Hunga Tonga–Hunga Ha'apai event rescaled in speed and amplitude. These simulations demonstrate that the meteotsunami surges can be higher than 1 m (and up to 10 m) along more than 7% of the world coastlines. Planetary meteotsunami waves thus have the potential to cause serious coastal damages and even human casualties during volcanic explosions or asteroid impacts either releasing intense acoustic energy or producing internal atmospheric gravity waves triggering the deep-ocean Proudman resonance at a speed of ∼212 m s−1. Based on records of catastrophic events in Earth's history, both scenarios are found to be realistic, and consequently, the global meteotsunami hazards should now be properly assessed to prepare for the next big volcanic eruption or asteroid impact even occurring inland. [ABSTRACT FROM AUTHOR]

There is growing evidence that many large coastal boulder deposits found on the exposed rocky ocean shores were deposited by extreme storm waves rather than by catastrophic tsunamis, as previously thought. In addition, before the first discovery in the northern Adriatic a few years ago, such deposits were not expected in relatively shallow semi-enclosed inland basins. Here we report on a large coastal coarse-clast deposit on the central Adriatic island of Mana, which also contains numerous large storm boulders that weigh up to several tons. Large solitary boulders are also located outside of the deposit, closer to the partly submerged sea cliff and the wave impact. The erosion of the cliff top and displacement of the carbonate bedrock fragments began when the extreme waves inundated the lowermost part of the cliff edge, probably during the late Holocene sea-level rise. The UAS photogrammetry-based fragmentation analysis of the storm deposit and the calculated fractal dimension value indicate that the material was fragmented by multiple high-energy events. A comparison of the available photographs indicates that displacements of the most exposed solitary boulders probably occurred during Vaia, the last extreme storm that hit the Adriatic on the 29th of October 2018. However, the modeled maximum wave height south of Mana during the peak of the storm would be insufficient to move these boulders. Yet local geomorphology probably further influenced the increase in wave height that, in combination with specific geological features, caused displacements of the boulders. There is a shorter fetch affecting Mana Island with respect to the northern Adriatic boulder field in southern Istria. Thus, such an active local erosion of the generally stable eastern Adriatic karstic coast depends on the extreme storms that have a weaker impact in the central than in the northern Adriatic. [ABSTRACT FROM AUTHOR]

The paper documents the use of a Digital Elevation Model (DEM) method to reconstruct paleochannels and the paleo-coastline during the Last Glacial Maximum (LGM) in the Central Eastern Adriatic area. We focused on the paleo-coastline and paleo-channels of the Neretva and the Cetina rivers, which were estimated from the 15" bathymetry available for the Adriatic Sea. While being aware of the limitations of the method and the resolution of the bathymetry grid, we successfully reproduced the paleo-channels of both rivers. Results for the Cetina River indicate the presence of depressions that were filled with water along its flow. The configurations of existing seabeds in the vicinity of the Cetina and Neretva River mouths indicate morphologies, similar to river mouths. The vertical profiles suggest that during the LGM the sea level was about 115 m lower than today. The total length of the Neretva riverbed was longer for about 136 km than today, and the Cetina River was approximately 154 km longer. [ABSTRACT FROM AUTHOR]

Large-amplitude seiches (up to 40 cm) were recorded at a tide-gauge inside Bakar. Bay from 19-21 August 2006. They were analysed using available tide-gauge and meteorological (air pressure and wind, vertical sounding,) data. Sea level and air pressure spectra showed simultaneous increases of energy in the periods in which the seiches of both Bakar Bay (about 20 mm) and Rijeka Bay (ca. 2 h) occurred. The propagation speed and direction as well as the horizontal span of the four pronounced pressure disturbances, favourable for coupling with Rijeka Bay waters via the Proudman resonance, were determined. It is argued that only short atmospheric pressure disturbances may generate significant barotropic sea waves in a small basin such as Rijeka Bay; once generated they are amplified when hitting Bakar Bay due to the harbour resonance. Furthermore, examination of the synoptic situation and upper-air soundings revealed the presence of dynamically unstable layer at a height of around 2000 m which was also the height of temperature inversion that developed as a consequence of an advection of warm African air from the south-east. During the afternoon hours of 19 August 2006 wind speed in this layer was close to the estimated speed of propagation of the air pressure disturbances which, together with dynamical instability, made conditions favourable for long distance propagation and trapping of atmospheric gravity waves. [ABSTRACT FROM AUTHOR]

Sea-level observations provide information on a variety of processes occurring over different temporal and spatial scales that may contribute to coastal flooding and hazards. However, global research on sea-level extremes is restricted to hourly datasets, which prevent the quantification and analyses of processes occurring at timescales between a few minutes and a few hours. These shorter-period processes, like seiches, meteotsunamis, infragravity and coastal waves, may even dominate in low tidal basins. Therefore, a new global 1 min sea-level dataset – MISELA (Minute Sea-Level Analysis) – has been developed, encompassing quality-checked records of nonseismic sea-level oscillations at tsunami timescales (T<2 h) obtained from 331 tide-gauge sites (10.14284/456, Zemunik et al., 2021b). This paper describes data quality control procedures applied to the MISELA dataset, world and regional coverage of tide-gauge sites, and lengths of time series. The dataset is appropriate for global, regional or local research of atmospherically induced high-frequency sea-level oscillations, which should be included in the overall sea-level extremes assessments. [ABSTRACT FROM AUTHOR]

In this evaluation study, the coupled atmosphere–ocean Adriatic Sea and Coast (AdriSC) climate model, which was implemented to carry out 31-year evaluation and climate projection simulations in the Adriatic and northern Ionian seas, is briefly presented. The kilometre-scale AdriSC atmospheric results, derived with the Weather Research and Forecasting (WRF) 3 km model for the 1987–2017 period, are then thoroughly compared to a comprehensive publicly and freely available observational dataset. The evaluation shows that overall, except for the summer surface temperatures, which are systematically underestimated, the AdriSC WRF 3 km model has a far better capacity to reproduce surface climate variables (and particularly the rain) than the WRF regional climate models at 0.11 ∘ resolution. In addition, several spurious data have been found in both gridded products and in situ measurements, which thus should be used with care in the Adriatic region for climate studies at local and regional scales. Long-term simulations with the AdriSC climate model, which couples the WRF 3 km model with a 1 km ocean model, might thus be a new avenue to substantially improve the reproduction, at the climate scale, of the Adriatic Sea dynamics driving the Eastern Mediterranean thermohaline circulation. As such it may also provide new standards for climate studies of orographically developed coastal regions in general. [ABSTRACT FROM AUTHOR]

Sea-level observations provide information on a variety of processes occurring over different temporal and spatial scales that may contribute to coastal flooding and hazards. However, global research of sea-level extremes is restricted to hourly datasets, which prevent quantification and analyses of processes occurring at timescales between a few minutes and a few hours. These shorter period processes, like seiches, meteotsunamis, infragravity and coastal waves, may even dominate in low-tidal basins. Therefore, a new global 1-minute sea-level dataset - MISELA (Minute Sea-Level Analysis) - has been developed, encompassing quality-checked records of nonseismic sea-level oscillations at tsunami timescales (T<2h) obtained from 331 tide-gauge sites (https://doi.org/10.14284/456, Zemunik et al., 2021b). This paper describes data quality-control procedures applied to the MISELA dataset, world and regional coverage of tide-gauge sites and lengths of time-series. The dataset is appropriate for global, regional or local research of atmospherically-induced high-frequency sea-level oscillations, which should be included in the overall sea-level extremes assessments. [ABSTRACT FROM AUTHOR]

In this evaluation study, the coupled atmosphere-ocean Adriatic Sea and Coast (AdriSC) climate model, which was implemented to carry out 31-year long evaluation and climate projection simulations in the Adriatic and northern Ionian seas, is briefly presented. The kilometre-scale AdriSC atmospheric results, derived with the Weather Research and Forecasting (WRF) 3-km model for the 1987-2017 period, are then thoroughly compared to a comprehensive publicly and freely available observational dataset. The evaluation shows that overall, except for the summer surface temperatures which are systematically underestimated, the AdriSC WRF 3-km model has a far better capacity to reproduce the surface climate variables (and particularly the rain) than the WRF regional climate models at 0.11° of resolution. In addition, several spurious data have been found in both gridded products and in situ measurements which thus should be used with care in the Adriatic region for climate studies at local and regional scales. Long-term simulations with the AdriSC climate model, which couples the WRF 3-km model with a 1-km ocean model, might thus be a new avenue to substantially improve the reproduction, at the climate scale, of the Adriatic Sea dynamics driving the Eastern Mediterranean thermohaline circulation. As such it may also provide new standards for climate studies of orographically-developed coastal regions in general. [ABSTRACT FROM AUTHOR]

This study quantifies the hazard associated with extreme sea levels due to eastern Adriatic meteotsunamis—long waves generated by traveling atmospheric disturbances—and assesses the sensitivity of the ocean response to the disturbances responsible for those events. In this spirit, a surrogate model of meteotsunami maximum elevation based on generalized polynomial chaos expansion (gPCE) methods, is implemented. The approach relies on the definition of a synthetic pressure disturbance—depending on six different stochastic parameters known to be important for meteotsunami generation, which is used as forcing to produce series of meteotsunami simulations defined with sparse grid methods (up to 10 689 used in this study). The surrogate model and the sensitivity study are then obtained with a pseudo-spectral approximation (PSA) method based on the chosen meteotsunami simulations. This study mainly presents the developed methodology and discusses the feasibility of implementing such gPCE-based surrogate models to assess the hazard and to study the sensitivity of meteorologically driven extreme sea levels. [ABSTRACT FROM AUTHOR]

The paper reports an unusual response of the microzooplankton community to oceanographic conditions observed during the winter of 2015 in the open southern Adriatic. Record-breaking nauplii abundance of 13,734 ind. m-3 was sampled for the open southern Adriatic by 50-μm net sampling. This could be explained by (i) warmer-than-usual surface and intermediate ocean temperatures, (ii) higher precipitation that freshened and widened the surface layer, pushing saline Levantine Intermediate Water below a 400-m depth, and (iii) strong wind episodes that transported nutrients from the coastal zone to the open ocean and induced limited vertical mixing. Neritic tintinnids and viable photoautotrophs well below the photic zone support the last possibility. Average seasonal maximum of zooplankton abundance was shifted from spring and early summer to the late winter months. Our results document large and fast variations in production conditions rarely found to occur in oligotrophic waters such as the southern Adriatic Sea. [ABSTRACT FROM AUTHOR]

Climate predictions of air temperature in coastal regions represent a great challenge due to the complex interactions among the atmosphere, sea, and land. With approximately 1,200 islands, the Adriatic is a region with a strong land‐sea contrast, land‐atmosphere feedback, and intense air‐sea interaction. Because the Mediterranean has been regarded as a "hot spot" for climate change, regional climate models can be used to provide insight into a more realistic representation of small‐scale weather and climate structure and variability. This advantage is due to the better representation of complex topography, developed coastlines, and land‐sea contrasts, which are important for investigating air temperature trends. The use of regional climate models together with high‐resolution reanalyses and observations in assessments of climate variability and climate change is highly valuable for understanding climate processes at regional scales. The present study focused on air temperature and its trends calculated from measurements and simulated by eight regional climate models from the EURO‐CORDEX database; these data were represented by UERRA reanalyses and E‐OBS gridded data. In the evaluation period (1989–2008), the models' RMSEs were fairly small, in the range 0.5–1.5°C, compared to the historical period (1961–2005), with RMSEs greater than 1.75°C. However, the models showed small absolute trend differences (up to 0.12°C·decade−1 for the historical period). The ensemble means in both periods showed an accuracy improvement of 15–20% compared to the individual models. The models exhibited more success in terms of representing the main statistics and variability of the air temperature structure than in reproducing the temperature trends over 45 years, especially in the northern Adriatic, where there is complex coastal topography and significant seasonal variability in the wind regime. The reanalyses well represented temperature structure but showed less success in explaining the temperature trends than the results from the measured data. [ABSTRACT FROM AUTHOR]

This study investigates a coastal boulder deposit that was recently identified along the Northern Adriatic coast (Premantura Promontory, Istria, Croatia). Accumulations of large boulders have not previously been reported in the North Adriatic Sea, which is a semi-enclosed basin elongated in the SE-NW direction. In particular, we devoted our attention to the mechanisms that may be responsible for the detachment and transport of these large limestone fragments from the emergent part of the coast and from the sea bed towards inland areas. We adopted a multidisciplinary approach including geological and geomorphological surveys, climate analysis, hydrodynamic modelling together with the use of an Unmanned Aerial Vehicle (UAV). Moreover, 14C AMS datings have been made in order to possibly constrain the age of the studied boulders.The deposit is composed of 950 clasts and lies on a low-lying limestone promontory where the topography, together with the bedding planes and dense joint pattern constitute the predisposing factors for boulder size and detachment.We focused our attention on an isolated boulder characterized by fresh marine carbonate encrustations. Its arrival was reconstructed by means of historical satellite images, wave data and hydrodynamic models and was ascribed to a severe storm which occurred during early 2014. Concerning the remaining boulders, thanks to 14C AMS datings, their deposition can be explained by multiple past extreme wave events.We explained the boulders detachment mechanism as due to frost and salt weathering during severe bora conditions, when air temperatures are below 0°C and waves and associated sea spray fill limestone cracks with water. The freezing and subsequent thawing caused further ruptures and eventually detachment of limestone fragments. Then, the boulders are carried onshore during major sirocco events when much longer and higher waves are generated. We estimate that sirocco generated wave heights can exceed ~15 m, which is enough to transport most of observed boulders to their present day location, according to hydrodynamic equations results. We hypothesize that the boulder detachment mechanism was more effective during the past, when the coastal slope was intact and the boulders were progressively removed and thrown onto the upper storm berm that represents the maximum inundation limit reached by the boulders. Through the centuries the coast has been eroded through the removal of rocks and has been shaped by truncations of the gently dipping limestone beds that acted as ramps for boulder movements along channels delimited by scarps. The mechanism might have also been favoured over the centuries by slow sea level rise inundating new portions of the rocky coast, bringing seawater into contact with limestone beds that were previously stable, in a subaerial environment, causing their resizing and rupture. [ABSTRACT FROM AUTHOR]