Sea Level Rise Scenario for 2100 A.D. in the Heritage Site of Pyrgi (Santa Severa, Italy)

: Sea level rise is one of the main risk factors for the preservation of cultural heritage sites located along the coasts of the Mediterranean basin. Coastal retreat, erosion, and storm surges are posing serious threats to archaeological and historical structures built along the coastal zones of this region. In order to assess the coastal changes by the end of 2100 under the expected sea level rise of about 1 m, we need a detailed determination of the current coastline position based on high resolution Digital Surface Models (DSM). This paper focuses on the use of very high ‐ resolution Unmanned Aerial Vehicles (UAV) imagery for the generation of ultra ‐ high ‐ resolution mapping of the coastal archaeological area of Pyrgi, Italy, which is located near Rome. The processing of the UAV imagery resulted in the generation of a DSM and an orthophoto with an accuracy of 1.94 cm/pixel. The integration of topographic data with two sea level rise projections in the Intergovernmental Panel on Climate Change (IPCC) AR5 2.6 and 8.5 climatic scenarios for this area of the Mediterranean are used to map sea level rise scenarios for 2050 and 2100. The effects of the Vertical Land Motion (VLM) as estimated from two nearby continuous Global Navigation Satellite System (GNSS) stations located as close as possible to the coastline are included in the analysis. Relative sea level rise projections provide values at 0.30 ± 0.15 cm by 2050 and 0.56 ± 0.22 cm by 2100 for the IPCC AR5 8.5 scenarios and at 0.13 ± 0.05 cm by 2050 and 0.17 ± 0.22 cm by 2100, for the IPCC Fifth Assessment Report (AR5) 2.6 scenario. These values of rise correspond to a potential beach loss between 12.6% and 23.5% in 2100 for Representative Concentration Pathway (RCP) 2.6 and 8.5 scenarios, respectively, while, during the highest tides, the beach will be provisionally reduced by up to 46.4%. In higher sea level positions and storm surge conditions, the expected maximum wave run up for return time of 1 and 100 years is at 3.37 m and 5.76 m, respectively, which is capable to exceed the local dune system. With these sea level rise scenarios, Pyrgi with its nearby Etruscan temples and the medieval castle of Santa Severa will be exposed to high risk of marine flooding, especially during storm surges. Our scenarios show that suitable adaptation and protection strategies are required.


Introduction
Observational and instrumental data collected worldwide since the last two-three centuries show that the global sea level is continuously rising with an accelerated trend in recent years, which coincides with the rise in global temperatures [1]. Global mean sea level is expected to rise by about 75 to 200 cm by 2100 in the worst scenarios [2][3][4][5], i.e., the most serious effects of climate change that might occur in future decades. These values will be even larger in subsiding coasts of the Mediterranean, entailing widespread environmental changes, coastal retreat, marine flooding, and loss of land, which will be disadvantages for human activities. The sea level rise will amplify the impacts exerted by a multitude of hazards (i.e., storm surges, flooding, coastal erosion, and tsunamis) on the infrastructure and building integrity, people safety, economic assets, and cultural heritage.
Therefore, it is important to mitigate these risks by providing multi-temporal scenarios of expected inland extension of marine flooding as a consequence of the sea level rise for a cognizant coastal management [6][7][8].
This is particularly true for the Mediterranean region where ancient civilizations were born and developed along its coasts [9,10]. A large number of heritage sites are located at the waterfront or very close to the sea level and are exposed to marine flooding under the effects of ongoing climate change. A large part of these sites, which are dated back to the Greek, Roman, and Medieval ages, are exposed at increasing risks to coastal hazards that are related to a sea-level rise [10].
Aerial photogrammetric surveys performed by small Unmanned Aerial Vehicles (UAVs) can provide accurate topography at low costs and, in short time, for small areas with respect to conventional aerial surveys [11]. Results, when analyzed in combination with sea level rise projections and vertical land movements (VLM), can support the realization of the expected sea-level rise and storm surge scenarios for future decades [12].
In this study, we show an effective application for the coastal archaeological area of Pyrgi, Italy, which is near Rome (Figure 1). High resolution maps of expected flooded areas and coastal positions for 2100, even in storm surge conditions, are reported in this paper. Our approach is to apply a multidisciplinary methodology previously tuned in the savemedcoasts project (www.savemedcoasts.eu), which includes topography, geodesy, sea level data, and climatic projections to estimate realistic sea level rise scenarios for targeted coastal areas. Our approach provides a useful analytical tool to identify the best adaptation and defense strategies against the sea level rise impact, and to protect heritage sites. The results help decision-makers in the selection of the best practical actions aimed at preserving the archaeological and historical sites located in coastal areas that are subjected to sea level-related risks. The proposed methodology can be exported in other areas of the Mediterranean region and beyond its borders.

Geo-Archaeological Setting
The heritage site of Pyrgi is located along the coasts of Northern Latium, between the villages of Santa Severa and Cerveteri, which is about 50 km north of Rome (Italy) (Figure 1). The area includes the Castle of Santa Severa that, with Pyrgi, is one of the most important heritage sites of the Tyrrhenian coast. The area has been settled since the V-IV millennium B.C. [13] and continuously developed during the Neolithic Age, during the Bronze Age (II millennium B.C.), and during the Iron Age (IX-VIII century B.C.), thanks to its good environmental conditions. In the Etruscan phase (VII-IV century B.C.), Pyrgi was the port of the ancient Etruscan city of Caere and played an important role in the maritime commerce being frequented by Greeks and Phoenicians ships. The area includes a sanctuary and the temples of Eileithyia-Leukothea and Apollo, Cavatha, Suri, and the Etruscan Uni analogue to the Phoenician Astarte [13].
After the Romanization of this area (III century B.C.), Pyrgi became a maritime colony and a tall fortress surrounded by a polygonal wall was built on part of the Etruscan settlement. During the Roman imperial age, the city of Pyrgi continued to be frequented until the 5th-6th century A.D. when a Byzantine castrum was built on its remains with the early Christian church of Santa Severa inside.
Later, the medieval and renaissance village became a large farm located in a strategic position between the main harbors of Rome and Civitavecchia [13][14][15][16]. In terms of its geological setting, the coast of Pyrgi belongs to the roman co-magmatic province [17] that underwent major volcanic activity during the Plio-Pleistocene era ( Figure 2). The surface geology is characterized by a bedrock belonging to the allochthonous Outer Ligurids [18], represented by the Tolfa Formation, spanning from the Late Cretaceous to Palaeogene [19], defined as the Pietraforte formation and "comprehensive succession" [20]. The latter consists of a series of marl limestone and grey marl beds that outcrop between Rome and Civitavecchia underlying the Pietraforte unit. Biogenic sandstones of the Early-Middle Pleistocene Panchina Formation partly overlies the Pietraforte [21]. .
The Holocene deposits are represented by travertines, slope debris, alluvial, weathering deposits, gravels, and sandy beaches [22]. The neotectonic of this sector of the Tyrrhenian coast of Italy is marked by the elevation of the MIS 5.5 marine terraces that show stability and a weak uplift of the inland sector for the last 124 years [23,24], which is related to magmatic injections under the Vulsini and Sabatini volcanic complexes [24]. Reference [25] underlines that the long-term uplift may not be an appropriate description for all the past two millennia since some weak subsidence may have occurred at Pyrgi and along the nearby coasts.

Methods
We applied a multidisciplinary approach using coastal topography, geodesy, and climatic-driven estimates of the sea-level rise to provide maps of flooding scenarios for the year 2100 A.D. for the coast of Pyrgi. Our study consists of three main steps: 1) the realization of UAV surveys to obtain an ultra-high-resolution orthophoto and a DSM model of the coastal area to map the current and the projected coastline positions including sea level data in the analysis, 2) the estimation of the current vertical land movements from the analysis of geodetic data from the nearest GPS stations, and, 3) by combining these data with the regional IPCC-AR5 projections (RCP-2.6 and RCP-8.5 scenarios), the calculation of the upper bounds of the expected sea levels for the targeted epochs of 2050 and 2100 A.D. and the corresponding expected inland extent of the marine flooding and shoreline positions were calculated. Lastly, storm surge scenarios were implemented for return times of 1 and 100 years, for sea level rise conditions.

Digital Terrain Model Reconstruction
To realize the ultra-high-resolution Digital Surface Model (DSM), an aerial photogrammetric survey was performed using a radio-controlled multi-rotor Da Jiang Innovation (DJI) Phantom 4pro UAV system, equipped with a high resolution lightweight digital camera, to capture a set of aerial images of the investigated area (Table 1). The UAV was controlled by an autopilot system using waypoints previously planned by PIX4D ® capture IOs App as a Ground Control Station system. To optimize the photogrammetric spatial resolution and coverage of the surveyed area, a constant altitude of 70 m was maintained during the flight and 306 partly overlapping (70% of longitudinal and 70% lateral overlapping between subsequent photos) aerial digital photos were acquired during three successful flights of 13, 11 and 15 minutes of duration each. To scale the aerial images, we used a dual-frequency geodetic Global Positioning System Real Time Kinematic (GPS/RTK) receiver STONEX S900A ® to measure the coordinates of a set of reference Ground Control Points (GCPs) falling in the investigated areas. GCPs positions were estimated in real time by the RTK technique with about 1 to 2 cm of accuracy, with respect to the reference GPS station TOLF. The latter is part of the GNSS network of the Istituto Nazionale di Geofisica e Vulcanologia (INGV) [26] and also pertains to the Leica SmartNet ItalPoS network (https://hxgnsmartnet.com/en-gb/) for real-time positioning services. We used the Agisoft PhotoScan ® software package (http://www.agisoft.com) based on the Structure-from-Motion photogrammetry technique [27] to process the acquired georeferenced images. The analysis included: 1) camera alignment with image position and orientation, 2) generation of a dense points cloud, and 3) generation of an orthophoto covering a land surface of 0.226 km 2 with a ground resolution of 1.94 cm/pixel as well as the DSM creation in the WGS84-UTM32 coordinate system.
The obtained ortho-rectified images (orthophotos) and digital elevation models were also managed by Global Mapper ® and ESRI ArcGis ® software to create cross sections, slope maps, surfaces, and coastline positions as well as calculate the dimension of the potential flooded areas. The extracted Digital Surface Model has a resolution of 15.5 cm/pixel and a point density of 41.5 points/m 2 As targets for GCPs measured by GPS/RTK, we used a set of i) natural markers belonging to fixed structures (e.g., the center of manholes, wall and sidewalk corners, small structures, or large stones), and ii) mobile targets for the time of the flight, such as thin metal crosses of 60 × 60 cm of size, deployed in the investigated area. All these GCPs were chosen as areas recognizable on the images during data analysis.
In total, 22 Ground Control Points (GCPs) and Check Points (CPs) were used to geo-reference the orthophoto with the ground control point toolset of AGISOFT photo scan ( Figure 3). The 3D coordinates of these points have been estimated with a mean RMS of 0.6 and 0.9 cm for the horizontal and vertical components, respectively, and were used to evaluate the vertical accuracy of our final DSM. Figure 3 shows the orthophoto while Figure 4 shows the DSM.

Tidal Correction and Coastline Position
The coast of Pyrgi, similar to most of the coasts of the Mediterranean Sea, is characterized by a micro-tidal environment and tides are generally in the range of ±30 cm. Only in the Gulf of Gabes (Tunisia) and the North Adriatic Sea (Italy), tides may reach amplitudes up to about 2 m [27][28][29][30]. We used the tidal data collected by the Italian tidal network managed by the Italian Institute for Environmental Protection and Research (ISPRA) (data are freely available at www.mareografico.it) at the sea level station of Civitavecchia (located at LAT 42°05'38.25'', LON 11°47'22.73'), which is placed near Pyrgi (Figure 1), to estimate the tide level (TL) and the local mean sea level (LMSL) at the time of the UAV surveys (26 April, 2019 at 07:30 or 05:30 UTC). We considered the complete time series to account for a long-term linear trend, representative of the mean sea level during UAV surveys, with respect to which the TL is defined. This tide gauge station shows a valid recording period of about eight years (2011-2019) with a sea level trend of 0.25 ± 0.1 mm/a, which is calculated from a linear fit on the monthly data ( Figure 5a). From the data analysis, a mean tide amplitude of ~35 cm and a maximum tidal range up to ~60 cm has been estimated ( Figure 5b). To define a reference level for elevation data, the mean sea level was computed propagating the linear trend from the time of surveys, assuming the mean sea level as a reference value for the year 2018. A mean sea level of 4.8 ± 11.8 cm above the topographic benchmark was estimated from the hourly tidal data. Since the UAV surveys have been performed during a high tide of 4 cm, as shown by the tidal recordings (www.mareografico.it), then the reference sea level at the time of surveys corresponds to a position of only 0.8 cm above the mean sea level for 2018.
This value is negligible compared to the DSM accuracy, so we did not consider it for the analysis of sea level rise scenarios. Lastly, we used the local mean sea level calculated at the tide station of Civitavecchia as elevation data, given the small distance between this location and Pyrgi. The obtained value was used to define the position of the coastline during the surveys and the one expected from the sea level rise scenarios to 2100 A.D. We preferred to adopt this local vertical datum instead of the value of the geoid elevation for the Italian region provided by the International Service for the Geoid-ISG-ITG2009 (http://www.isgeoid.polimi.it/Geoid/Europe/Italy/ITG2009_g.html) [31] since it is an independent and more accurate elevation datum. The International Service for the Geoid (ISG) estimates a geoid height in the Italian Geoid (ITG) 2009 for the coastline of Pyrgi at 48.319 m. This elevation corresponds to the contour line equal to zero in the Italian height reference frame. The DSM height reference frame is 0.42 m, as estimated by GPS/RTK data at a GCP located at the sea level along the coastline. Its elevation was corrected for the tidal range at the time of the surveys. The reference mean sea level estimated by the tidal analysis provided a local mean sea level with an uncertainty of ±11.8 cm.

Vertical Land Motion (VLM) at Pyrgi
The current rate of VLM at Pyrgi was estimated by the analysis of the available GPS data collected at the nearest GNSS stations of TOLF, belonging the INGV Rete Integrata Nazionale GPS (RING) network (DOI:10.13127/RING) and MAR8, belonging to the Topcon-NetGeo network (http://www.netgeo.it). These stations, which are located at about 6.5 km and 0.3 km of distance from the study site, respectively (Figure 1 GPS data analysis has been carried out following the procedures already described in Reference [32] and updated in Reference [33] by adopting a three-step procedure using the GAMIT/GLOBK V10.7 [34] and QOCA software. This is part of a continental-scale GPS solution including >3000 stations [34]. The daily positions of TOLF and MAR8 have been estimated in the GPS realization of the ITRF2008 frame [35], i.e., the IGb08 reference. The position time series have been analyzed in order to estimate and correct offsets due to station equipment changes while, simultaneously, estimating annual and semi-annual periodic signals and a linear velocity term, whereas velocity uncertainties have been estimated adopting a power law + white noise stochastic model, as in Reference [36]. The results show that both sites are relatively stable in the IGb08 reference frame with a vertical velocity of −0.061 ± 0.135 mm/year for TOLF and −0.456 ± 0.344 mm/year for MAR8. We remark that uncertainties associated on the vertical velocities are about ±0.5 mm/year and are barely significant in view of unresolved questions about the GPS reference frame stability and additional factors [30]. In addition to GPS data, the tectonic stability of this region is also inferred from the low level of seismicity deduced from historical data [37] and instrumental recordings of earthquakes (www.ingv.it), which do not report the occurrence of significant events for the last 3000 years BP.
Lastly, assuming that the area will continue in the near future to have the same tectonic trend shown in the past, it is reasonable to neglect the contribution of VLM in the sea level rise projections and flooding scenarios for 2100 A.D.

Relative Sea-Level Rise Projections and Flooding Scenarios for 2050 and 2100 A.D.
To estimate the sea-level rise for 2050 and 2100 A.D. at Pyrgi, we referred to the regional IPCC AR5 sea-level projections discussed in the Fifth Assessment Report of the IPCC-AR5 [3, www.ipcc.ch] (data available from the Integrated Climate data Center-ICDC of the University of Hamburg, http://icdc.cen.uni hamburg.de/1/daten/ocean/ar5-slr.html). These data consist of the sea-level ensemble mean values and upper/lower 90% confidence bounds of the sea level on a global grid (spatial resolution 1°×1°), obtained by adding the contributions of several geophysical sources driving long-term sea-level changes: 1) the thermosteric/dynamic contribution (from 21 CMIP5 coupled atmosphere-ocean general circulation models AOGCMs), 2) the surface mass balance and dynamic ice sheet contributions from Greenland and Antarctica, 3) the glacier and land water storage contributions, 4) the Glacial Isostatic Adjustment (GIA), and 5) the inverse barometer effect [1]. Projections, which are based on two different Representative Concentration Pathways RCP 2.6 and RCP 8.5 while providing the least and most amounts of future sea level rise, respectively, were used. The IPCC regional sea-level rate at the grid point closest to the location of the tide gauge station (Civitavecchia) was considered. By accounting for VLM from GPS data, very high-resolution DSM and regional IPCC sea level projections at the grid point corresponding to the investigated area, the first marine flooding scenarios for Pyrgi for 2050 and 2100 A.D. have been realized. To include the VLM effect in sea-level projections, we substituted the modelled GIA contribution to the IPCC rates with the GPS vertical velocities, which includes both GIA and tectonic components. Uncertainties for the sea-level estimations were calculated by combining lower and upper sea level bounds from IPCC projection and errors from GPS measurements. In any case, given the tectonic stability of the area, the VLM have a null contribution in the analysis. The relative sea-level rise in RCP2. 6    The projected coastline positions for 2100 A.D., corresponding to the RCP2.6 and RCP8.5 scenarios, are obtained from the DSM for the sea levels listed in Table 2 and shown in Figure 8. The computed and the represented scenarios correspond to the local mean sea level (estimated with the uncertainty of ±11.8 cm), and are obtained by neglecting the periodical contribution due to diurnal and semidiurnal tides. To account for it, in order to estimate the maximum inundation scenarios, the time series of daily hydrometric sea levels at Civitavecchia tidal station from January 2010 to December 2018 were included in the analysis. The average half amplitude of the daily tides was estimated as high as about 30 cm (Figure 5b). Consequently, in the RCP 8.5 scenario for 2100, if we take into account both the sea-level rise (56 cm) and the mean daily tides (30 cm), we can infer a maximum water level of 86 cm. Since the highest sea levels may reach 45 cm for a few hours a year during temporary meteorological conditions when a decrease in atmospheric pressure and wind push water from the sea onto the land (Figure 5b), the maximum expected water level may reach up to 101 cm above the current level. For the RCP 2.6 scenario for 2100, considering a sea-level rise of 30 cm and the same mean daily tides (30 cm), a maximum sea-level rise of 60 cm can be inferred. During the rare maximum highest sea levels (45 cm) (Figure 5b), the maximum expected water level would reach 75 cm above the 2019 level.
The sea level rise scenarios for 2100 A.D. depicts the impact of the marine flooding for the coast surrounding the castle of Santa Severa, the temples of Pyrgi, and nearby beaches. The cross sections shown in Figure 9 traced along a direction that is normal to the coast and across the archaeological area of Pyrgi, which highlights the presence of a soft dune system running parallel to the shore. This system is well developed south of the castle for the whole bay of Pyrgi, while, in the north side, it is buried by the modern buildings of the Santa Severa village.
Since the Etruscan time (2450 years BP), the coastline has retreated at about 80 m, at a mean rate of 3.3 cm/yr [38]. The temple of Pyrgi was originally placed at more than 120-125 m from the present-day shore and well protected from the sea by the dune system. Today, it is placed at about 40-45 m from the coastline and is undergoing water intrusion and frequent flooding, especially during storm surges (Figures 8 and 9). The proximity of the sea and the accelerated erosion of the dune system, which is highly exposed to waves especially during storm surges, are causing a continuous coastal erosion and beach retreat. Based on our sea level rise projections, a beach retreat of about 10 m is expected for 2100 A.D. (see cross sections B, C, and D in Figures 8 and 9).
This rapid retreat will accelerate the erosion process and the likely dismantling of the soft dune system with the consequent direct exposure of the temples to the sea. At the same time, the beach located along the north-west side of the castle, which is lacking a dune system and characterized by a gentle slope, is expected to retreat up to about 25 m (see cross section A), exposing the foot of the castle to severe erosion and flooding.
Lastly, Table 3 shows the expected flooded areas in 2050 and 2100 A.D. for the RCP 8.5 and RCP 2.6 climatic scenarios, together with the corresponding percentage of beach loss, which may reach up to 46% during the highest tides. Table 3. IPCC scenarios, relative sea level rise projections for 2100, expected land loss (m 2 ), and percentage of land loss with respect to the current surface for the given Relative Sea Level Rise (RSLR) projections. Given the vertical tectonic stability of the area, the VLM were considered in the analysis.

Storm Surge Scenarios
To evaluate the storm surge scenario for Pyrgi, a local climate wave assessment was performed using the Weibull distribution and the Equivalent Triangular Storm model originally proposed by Boccotti [39]. The return values of the significant wave heights and related peak wave periods (Table  4) have been evaluated on the parameters of the Weibull distribution [40] using data collected at the wave buoy of Ponza (wave data are freely available at http://dati.isprambiente.it/), which is located in the central Tyrrhenian sea in front of the coasts of Latium (at coordinates Latitude 40.867°N and Longitude 12.950°E). Table 4. Return values RT (years) of the significant wave height HS (m) and the related peak wave period TP (s). The climate wave data for ordinary and extreme storm wave conditions in terms of the return period (RT) for RT = 1 year and RT = 100 years in a sea level rise condition for 2100 in the RCP8.5 climatic scenario reported in Tables 2 and 3 were used to estimate wave setup and wave run-up during the event for both return times ( Figure 10).

RT (Years) HS (m) TP (s)
In the present analysis, wind set up was neglected and the precautionary assumption that the direction of wave travel does not produce refraction was considered. The expression derived by Weggel [41] and the relation of Komar and Gaughan estimated breaker depth index and breaker height index [42]. The breaker type was correlated with the surf similarity parameter, which allows the run-up estimation by the use of the predictive equations [43,44]. Based on our topographic surveys, a mean beach slope of 3.7% for the investigated coast, was assumed for the wave run-up assessment.  The results of the wave run-up analysis in ordinary and extreme storm conditions (Table 5 and Figure 11), were estimated at 3.37 m and 5.76 m for return time of 1 and 100 years, respectively.

Discussion and Conclusions
In this study, we assessed the marine flooding scenarios for 2050 and 2100 for the relevant heritage sites of Pyrgi on the basis of (a) a high-resolution DSM, (b) rates of VLM from GPS data, and (c) tidal data from the nearby tide gauge and the RCP 2.6 and 8.5 climatic scenarios released by the IPCC. The results highlight how the flooding might impact both the beach and archaeological sites with a potential local sea-level rise for the IPCC 8.5 scenario of about 56 cm for 2100 A.D. During the highest sea levels estimated by the tidal data at the nearby Civitavecchia station, the sea-level could potentially reach more than 86 cm and up to 101 cm. In addition, 2450 years ago, the shoreline was extended seaward of about 80 m more than today while, during the ancient Roman period, the sea level was at about 1.2 m below the current position, as estimated by the fish tank data [34]. Hence, a continuous sea level rise is occurring since the last millennia with an acceleration that started about 100±53 years ago, at the beginning of the Industrial Era [25]. Due to the potential significant impacts on both the coast and the heritage sites of the Temple of Pyrgi and the Santa Severa Castle with a beach retreat up to about 25 m, the expected scenario reported in this study can support adaptation plans at different time scales, which are in agreement with the Protocol on Integrated Coastal Zone Management (ICZM, https://eur-lex.europa.eu/eli/prot/2009/89/oj) in the Mediterranean.
Our results detail previous studies for the Italian [7][8][9]12,[45][46][47] and the Mediterranean [10,48] regions and can contribute to raise awareness of policymakers and heritage managers toward the coastal hazard by highlighting the need to adapt actions to protect Pyrgi from marine flooding and erosion under the current conditions due to the expected sea level rise and storm surge scenarios.
We remark that the unceasing coastal erosion since the Etruscan time [38] and the retreat of the soft cliffs characterizing the Pyrgi coastline mostly occurs during high energy marine events. During storm surges, the waves approaching the coast from the Northwest, West, and Southwest marine sectors are particularly dangerous for Pyrgi. The waves approaching the coast with fetches up to about 400 km are weakly slowed down by the seafloor morphology, which results in an increase of wave energy in sea level rise conditions [49] and leads to a maximum wave run-up of 5.76 m for a return time of 100 years. The impact of extreme events may be significantly amplified by the effects of climate change. In this scenario, the intensities of the strongest future storms will exceed the strength of any in the past in the Mediterranean Sea and the oceans [50][51][52][53]. The effects of the sea-level rise in the projected scenarios include severe coastal erosion and potential enhanced damages to the heritage site of Pyrgi, which, as many other coastal areas of the Mediterranean, deserves rapid actions for their future preservation.