The Importance of Geotechnical Evaluation and Shoreline Evolution in Coastal Vulnerability Index Calculations

The aim of this specific study is to present a new weighted Coastal Vulnerability Index (CVIWF), with an emphasis given to the geotechnical evaluation and shoreline evolution rate measured through high-resolution remote sensing, which seem to be the most interfering variables in CVI calculations. As a pilot area for the application of the new CVIWF, the Gulf of Patras in Western Greece was selected, which is suffering erosion problems due to climate change, the sea level rising and human intervention. The new CVIWF, which was applied in this research, includes the following innovations: (1) the use of geotechnical characterization instead of geological–geomorphological characterization, (2) the use of high-resolution remote sensing data for the detection of shoreline evolution rate and (3) the insertion of a specific weighted geotechnical factor in the CVIWF formula. The results from the application of the unweighted CVI show that percentages of 20.13%, 20.47%, 24.56%, 29.39% and 5.45% of the gulf’s shoreline are under the regime of very low, low, moderate, high and very high vulnerability, respectively. On the other hand, the corresponding results from the application of the weighted CVIWF show a percentage of 14.59%, 25.91%, 20.04, 36.48% and 2.98, respectively.


Introduction
The coastal zone can be considered as the geographical area of interaction between marine and land area (land-sea interface), which coexists with human societies, as well as socioeconomic and human activities. The rising sea level and extreme events related to climate change are causing severe threats to coastal areas, affecting both natural and human systems worldwide [1]. One of the most common and diachronic problems in coastal zones is sea-wave erosion, named coastal erosion. In Greece, more than 20% of its total shoreline [2] is threatened by coastal erosion and rising sea levels, ranking the country fourth in terms of coastal vulnerability among the 22 coastal EU member states. In this already complicated scenario, coastal areas are also experiencing relevant pressures resulting from multiple human-induced stressors linked with coastal economic development (e.g., touristic activities and infrastructures along the shoreline) and the connected land use changes (e.g., urbanization) [3]. Due to the continuous development and the increasing population living in the coastal zone (half of EU population lives in regions within 50 km of the coast) [4], the construction of various technical works is required, such as buildings, infrastructures, roads, ports and coastal defense works. However, coastal areas are vulnerable and susceptible to coastal hazards like storm surges, tsunamis, floods, coastal erosion, sea level rise, etc. Climate change, global warming and, as a result, the sea level rise (SLR) can cause accelerated erosion, shoreline retreat in coastal zones and saltwater intrusion into groundwater aquifers [5]. In a very recent study [6], it has been indicated that most coastal countries in Europe (including Greece) are planning for SLR with the horizon the year of 2100. For these reasons, the assessment of vulnerability to the impact of the rising sea levels and coastal hazards is an important aspect of coastal zone management and coastal protection [7].
In many countries and institutions, different risk assessment methodologies have been developed to understand the processes underpinning coastal erosion risks [8][9][10][11][12]. As a result, several tools and methods supporting policy and decision makers in the implementation of recommendations, and directives for coastal zone management have been developed in several countries [13][14][15], with one of the most important aspects being coastal vulnerability. The assessment and calculation of coastal vulnerability, focusing on coastal erosion processes, is accomplished with the following methods: (a) methods based on a specific vulnerability index, (b) methods based on a sequence of indicators, (c) methods based in geographical information system (GIS)-based decision-making systems and (d) methods based on dynamic computational models [3]. The most common used method is the vulnerability index, which is considered as a consistent, easy and fast way to calculate the vulnerability of coasts. This approach combines the "sensitivity" of the coastal system to changes (shoreline differentiation due to rising sea levels) with its natural ability to adapt to changing environmental conditions [16]. The most used and accepted method of calculating coastal vulnerability is the Coastal Vulnerability Index (CVI). The CVI was developed by Gornitz [16,17] and modified by Thieler and Hammer-Klose [18,19] for the evaluation of vulnerability throughout the United States. The most common parameters used in the CVI methodology are (a) geology and geomorphology, (b) coastal slope, (c) shoreline change rate, (d) average tidal range, (e) change in the sea level and (f) average significant wave height.
Many researchers have used vulnerability indexes or similar modified models for the assessment of coastal vulnerability to erosion in coastal zones [20][21][22][23][24][25][26][27]. For example, Diez et al. [21] changed the variable of coastal slope with elevation and inserted local tectonics affections, while Pramanik et al. [28] used seven parameters and different conventional and remotely sensed data. Rao et al. [7] added weighted values in the calculation formulas, and Ruzic et al. [29] improved the CVI for complex geologic conditions by adding a weighted factor in the parameter of geology. De Serio et al. [30] and Zhu et al. [31] performed a multicriteria evaluation method, named the Analytical Hierarchy Process (AHP). in the Adriatic Sea, Italy and in China, respectively. Greco and Martino [32] integrated in their CVI approach the study of time-dependent vulnerability together with information on the morphology and socioeconomic features, showing that vulnerability is not static in time but evolves according to changing patterns in the climate system. Lopez et al. [33] used four different sea level rise values corresponding to the four greenhouse gas emission scenarios given by the IPCC [34] to produce CVI maps on a national scale in Spain. Rangel-Buitrago [35] developed a more complex model using forcing, susceptibility and hazard indexes. A systematic review on various Coastal Vulnerability Indexes and Mapping was given by Kantameni et al. [36], Bukvic et al. [37] and Anfuso et al. [38].
Models of the vulnerability index have also been developed with different approaches that determine the coastal vulnerability using socioeconomic variables (e.g., population density, poverty levels and land use) [39][40][41][42][43], while other models use socioeconomic and physical-geographical parameters, such as the SIVECA index and CERA GIS-based tool [32,[44][45][46][47]. Szlafsztein and Sterr [48] introduced a Composite Vulnerability Index, allowing the combination of eight variables reflecting the natural dimension of vulnerability (e.g., coastline length, coastal protection measures and fluvial drainage) with seven other indicators reflecting the socioeconomic dimensions. Furthermore, Toressan et al. [49] introduced the methodology by integrating physical, geological and socioeconomic indicators in the CVI model and also taking into account the expected SLR and (storm surge level) flooding projections for the 2100 timeframe (i.e., climate-improved CVI).
In Greece, numerous studies have been conducted with the use of CVI or other adopted methods to assess the vulnerability to erosion. Gaki-Papanastasiou et al. [22] and Karymbalis et al. [24,25] assessed the vulnerability to an anticipated future sea-level rise, using the Coastal Vulnerability Index (CVI) and the Coastal Sensitivity Index (CSI), respectively. Doukakis [20] added in the CVI the parameter of subsidence for Western Peloponnese, while Tragaki et al. [50] considered socioeconomic and demographic parameters for the CVI model, for the evaluation of coastal hazard vulnerability in the same area. Alexandrakis and Poulos [51] developed a new index for the evaluation of vulnerability (Beach Vulnerability Index or BVI) mainly considering hydro and sediment dynamic parameters (longshore sediment transport, riverine inputs, storm surge etc.). Finally, Alexandrakis et al. [52] combined the environmental and economic approaches along the geographical space to relate the beach erosion vulnerability with the expected land loss and the relevant values from economic activities in Crete.
Several types of methods based on remote sensing low-resolution data derived from Landsat or Sentinel images have been used in order to estimate the shoreline change rate [53][54][55][56]. Accordingly, two recent studies sited in Northwestern Peloponnese by Apostolopoulos and Nikolakopoulos [57] and Apostolopoulos et al. [58] concluded that the accuracy of the vector data derived from low-resolution data are not suitable either for shoreline forecasting or shoreline mapping. An overall review and meta-analysis of the remote-sensing data, GIS methods, materials and indices used for monitoring the coastline evolution over the last twenty years was given by Apostolopoulos and Nikolakopoulos [59].
Another important variable that proved to be very important in the determination of coastal vulnerability is shoreline evolution, otherwise known as the historical shoreline change rate, with the application of extensive remote-sensing techniques. Generally, there are two main categories of studies regarding remote-sensing monitoring: those that use free-of-charge low-resolution data and those that use high-resolution data. Over time, many studies have tried to process shoreline evolution with low-resolution satellite images, but their results always lack accuracy. Different methods have been used to overpass this drawback regarding the accuracy of the measurements and the results they provide. Pardo Pasqual et al. [53] developed a very complicated methodology in order to ameliorate the Landsat data accuracy using a high-resolution air photo mosaic by co-registering the lowresolution coastline derived from Landsat data to the respective high-resolution shoreline. In another similar study, Almonacid et al. [60] computed the average annual coastline movement of many Landsat images and yielded a mean error of −4.7 m, showing that the shorelines derived from Landsat images were placed seaward in comparison with the respective high-precision shorelines. In a recent study in Northwestern Peloponnese, Apostolopoulos and Nikolakopoulos [57] followed a statistical analysis using Digital Shoreline Analysis System (DSAS) software and concluded that low-resolution data are not suitable either for shoreline forecasting or shoreline mapping, as their accuracy are very low, and they present a divergence that ranges between 6 and 11 m. In a similar study, Apostolopoulos et al. [58] concluded that the accuracy of the vector data derived from Sentinel-2 10-m spatial resolution images ranged between 1 and 5 m using the NDVI and NDWI indices, respectively. On the other hand, the synergy of high-resolution satellite imagery and high-scale aerial photos has become a very accurate combination in order to study the coastal zone changes. In the current study the shoreline rate of change was computed based on high-resolution aerial and satellite images, improving the rate's accuracy. These datasets were processed in a GIS platform, and the statistical analysis was performed through DSAS application by computing the end point rate (EPR). The use of high-resolution data provides an improved accuracy during (CVI) analysis and gives better results in finding the respective erosion or accretion trends along the shoreline.
On the other hand, the combination of historical cartography, aerial photographs and GPS surveys has become a very accurate combination in order to study the coastal zone changes. In a recent study, Antonello Aiello et al. [61] used this combination to assess the changes of the littoral environment between 1870 and 2005 on the coastal area of the Jonian Sea in Italy. In another study, Murray Ford [62] used vertical aerial photographs along with high-resolution satellite imagery such as IKONOS (0.8-m pixel size), QuickBird (0.6-m pixel size), WordView-2 (0.5-m pixel size) and GeoEye-1 (0.5-m pixel size) to estimate the shoreline changes in the Marshall Islands. Additionally, historical changes in the shoreline position along the Bay of Jijel, Algeria were studied by reference [63] with the combination of multi-dated aerial photographs and QuickBird (0.6m pixel size) satellite image. Furthermore, ref. [64] a combined semi-automatic technique for coastal risk assessment and monitoring integrated ground data with Earth Observation (EO) data in Apulia Region (South Italy). Finally, reference [65] used aerial photography, Google Earth Pro and Topographic survey measurements in order to determine changes in the shoreline position in four beaches on Valparaíso Bay in Central Chile.
Taking into consideration the above methodologies, it was decided, for the purposes of the current research, to investigate more thoroughly the interference of geotechnical evaluation and shoreline change rate in the preformed CVI calculations. For this reason, the calculation of coastal vulnerability was based on the original formula of Gornitz [16,17] with the following modifications: (1) the use of geotechnical characterization instead of geological-geomorphological characterization, (2) the use of high-resolution remotesensing data and (3) insertion of a specific geotechnical weighted factor into the CVI formula. All calculations were performed in the GIS; therefore, it finally produced a combined CVI-GIS model. The new CVI WF model was also compared with an unweighted (CVI) model, and for this comparison as a pilot area, the Gulf of Patras in Western Greece was used. The reason for selecting the Gulf of Patras as a pilot area had to do with the detailed data that has been collected from recent coastal monitoring activities in the frame of the TRITON project [66]. In order to proceed with the relevant CVI WF calculations, the pilot area was divided in eight coastal zones, and the geotechnical behavior of the prevailing geological formations was thoroughly investigated. The geotechnical investigation that was performed included borehole drilling, laboratory tests, Standard Penetration Tests (SPT), sediment sampling and engineering-geological mapping. This type of investigation improved the quantification of the geotechnical variable and led to a geotechnical characterization of the pilot area. Moreover, the shoreline rate of change was computed by using high-resolution aerial and satellite images, as well as the Digital Shoreline Analysis System [67] and by computing the end point rates (EPR) along the eight coastal zones of the pilot area. These rates were imported to the (CVI WF ) model and improved its accuracy. The final purpose of this work was mainly to highlight the importance of geotechnical evaluation in Coastal Vulnerability Index calculations, along with shoreline evolution.

Study Area
The Gulf of Patras is in Western Greece, Northwest Peloponnese and extends from the Rio-Antirio Bridge in the eastern part to the Cape of Araxos in the western part ( Figure 1). The total length of the southern shoreline of the gulf is around 45 km, with the maximum depth of water at 130 m. The Gulf of Patras is a Plio-Quaternary graben [68]. In the west it is connected with the Ionian Sea and, in the east, with the Gulf of Corinth. The main rivers discharging into the gulf are the Acheloos River and the smaller Evinos and Peiros rivers. The evolution of the gulf in the Quaternary occurs because of the interactions between the tectonic subsidence (rates at 3-5 mm/year at the central graben), global sea level rise and river sediment supply [69].
One of the most significant hazards that threatens the area of study is coastal erosion, and according to the Hellenic legislative framework, the Gulf of Patras has been characterized as a water body at risk. Due to the above description, it is easily understandable that the assessment of vulnerability to erosion is very important for this particular area, because it will assist local authorities with important information for its future protection against erosion.

Materials and Methods
The main features of the presented (CVI WF ) model are: (1) the use of geotechnical data in order to characterize the coastal formations as geotechnical units, (2) the significant mean wave height with a return period of 10 years, (3) the coastal slope, which has been calculated by using a Digital Elevation Model (DEM) (2X2 m), (4) the shoreline evolution by computing the relevant end point rates (EPR), (5) the average tidal range and (6) the sea level rise (SLR).
The calculations of the CVI were performed with two different methods: (1) with the use of the formula proposed by Gornitz [16] by computing the CVI as the square root of the product of the variables divided by their total number (n) and (2) with the use of the new CVI WF by adding a weighted factor in the geotechnical variable, which seems to be the most representative of the resistance to erosion. The calculation of the new (CVI WF ) model performed with the use of data that was derived from a marine equipment that was installed in the pilot area for the purposes of the TRITON project [66]. The data acquired from the TRITON project was (a) the significant wave height and (b) the tidal range measurements. The data were derived from a tide gauge system and a wave buoy that were installed in the frame of the TRITON project. The significant wave height also was acquired by performing a numerical simulation of wave propagation for wind speed with a return period of 10 years for each wind direction, using the MIKE 21, DHI software, which was also performed in the frame of the TRITON project.

Geotechnical Characterization
The most common variables used for the description of resistance to erosion are coastal geology and geomorphology. During this research, this variable was replaced with the geotechnical characterization parameter. For this purpose, a thorough geotechnical investigation was carried out in the research area, including the drilling of forty (40) boreholes, laboratory and in-situ geotechnical tests (SPT), a numerous coastal sediment sampling and field engineering-geological mapping. An SPT in-situ geotechnical test is widely used for a soil strength characterization. The calculation of SPT can be described as the number of blows required for each 15-cm penetration of the sampler into the soil, until it penetrates 45 cm. Finally, the result of this test is the number of blows required for the penetration of the last 30 (15+15) cm and is indicated as N or N SPT blow numbers [70].
The use of geotechnical data improves the quantification of the coastal resistance to the erosion parameter and is more reliable and representative, since it categorizes the total area according to the geotechnical behavior of the prevailing coastal formations. In Table  1 are presented some of the most representative geotechnical data that were used for the geotechnical characterization of the eight (8) zones of the pilot area. This method led to a new categorization depending on the geotechnical behaviors of coastal formations, which were ranked as "very low", "low", "moderate", "high" and "very high" vulnerability. Medium dense silty sand with gravels (SP-SM) and moderately dense silty sand with few gravels (SP-SM) 4 Kato Achaia Beach

Significant Mean Wave Height
The main variable that affects and forces the sediment to move and transport (longshore and cross-shore sediment transports) and causes erosion is the wave energy. High waves have bigger energy compared with low waves and, as a result, cause erosion with a higher rate. A combination of wave height, sea level rise and tidal range form the main hydrodynamic forces that cause coastal erosion and, especially, the wave height has very a high impact on coastal erosion and calculation of vulnerability. Wave height values for the coastal zone of the Gulf of Patras calculated with a numerical simulation of wave propagation for wind speed with a return period of 10 years for each wind direction using the MIKE 21, DHI software.

Coastal Slope
A shoreline with a high coastal slope is not vulnerable to SLR and shoreline retreat. However, in a flat coast (small coastal slope), a possible inundation has significant consequences, such as shoreline retreat, since a possible SLR will cause a big inundation of the sea inside the land and cover large areas. Ruzic et al. [29] inverted this rating based on the conception that high coastal slopes cause rockfalls and extreme erosion in the shoreline. The parameter of the coastal slope, during this research, was determined through a Digital Elevation Model (DEM) (2X2 m) acquired through the National Greek Cadastre and Mapping Agency and ranked with "low" vulnerability on steep coasts and "high" vulnerability on gentle coasts.

Shoreline Evolution
The shoreline evolution or change rate is an important parameter that indicates the accretion/erosion pattern in an area. The computation of a historical trend and monitoring of shoreline evolution in the pilot area was performed by using remote-sensing data processing for a 10-year period (from 2008 to 2018). High-resolution satellite data and orthomosaics were used to monitor the coastline's changes in the littoral area. In the current study, official datasets of orthomosaics for the years of 2008 and 2016 were acquired through the National Greek Cadastre and Mapping Agency with a spatial resolution of 0.5 m, which are the most accurate datasets for the Greek territory. Moreover, Worldview-2 high-resolution imagery (0.5-m spatial resolution) of the year 2018 was orthorectified and used to extract the shoreline for that year. All remote-sensing data used for the current research are illustrated in Table 2. The uncertainties, limitations and errors related to the shorelines extracted from maps and aerial photograph techniques have occupied the research community [71][72][73]. These refer to the five main routes of uncertainty during the shoreline vectorize procedure when using aerial photographs: seasonal error, tidal fluctuation error, digitizing error, pixel error and rectification error [72].
Seasonal and tidal fluctuation errors could considerate as negligible, since the images acquired during the summer seasons and the tide ranged from 0 to ± 10 cm, according to data derived from the online platform (https://www.worldtides.info/; accessed on: 30 September 2020) and the tide gauge system installed in the Gulf of Patras. As a relative example, we mention the high tide on the 18th of September 2018, the day of Worldview-2 image acquisition, which was estimated at almost 0.00 m. Considering the previous statement, any uncertainty is limited into the last three referred digitizing errors.
According to reference [74], the error associated with the pixel size of high-resolution imagery ranged from ± 0.3 to 1.0 m. As the datasets consist of high-resolution images with pixel sizes up to 0.5 m, we chose to proceed with the visual interpretation method. ERDAS IMAGINE 2014 and ArcMap v10.5 software were used for the processing of shoreline vectorizing, while DSAS application was used for the statistics and the assessment of the rate of change in every zone of the coastal area by computing the end point rates (EPR). The EPR calculates the rate of shoreline change dividing the distance between the oldest and most recent shorelines by the time elapsed [67].
After that, all digitized shorelines were produced in shapefiles and imported into a database in the ArcGIS platform in accordance with the Digital Shoreline Analysis System (DSAS v5.0) standards, such as the acquisition date, identity, shape, length and uncertainty [75]. The Digital Shoreline Analysis System is a tool that cooperates with the ArcMap software package and has been created by the United States Geological Survey. It is a popular tool, which calculates the changes in shoreline movements due to some various statistical rates [67]. Moreover, a smooth baseline was carefully created landward and used for the calculations. Several transects were created perpendicular to the baseline with fiftymeter intervals. The outlier rates were removed, and the remaining rates were statistically processed. With this tool, the vectorized shorelines for each zone were processed and computed as end point rates (EPR) for the period of 2008-2018.
The current study focused on the recent (last twelve years) evolution of the coastline, and for the specific period, we selected the most accurate available datasets. As described in Section 3.4, the three datasets used in the current study have a spatial resolution of 0.5m. The 2008 and 2016 are the official datasets of the Greek Cadastral. They are the most accurate datasets for the Greek territory. The 2018 worldview image also had a 0.5-m spatial resolution.

Average Tidal Range
According to Gornitz et al. [16], macro-tidal environments (big tidal range) are vulnerable to erosion, because they can create stronger currents and cause bigger sediment transportation in comparison with the micro-tidal environments. According to Thieler and Hammer-Klose [19] and Pendleton et al. [5], macro-tidal environments are considered as nonvulnerable due to the fact that during a storm surge in a micro-tidal coastal zone, the sea level is all the time close to a high tide and is more vulnerable to flooding.
The Gulf of Patras has been characterized as a micro-tidal environment, since the tidal change on the coasts is small. During this specific study, the small tide range ranked with low vulnerability, since the calculation of the CVI is against erosion and sea level rise and is not vulnerability against flooding and storm surges. The data and measurements of the tidal range in the Gulf of Patras were taken from a tide gauge system that has been installed in the port of Patras, with the average tidal range being equal to 0.24 m.

Sea Level Rise (SLR)
SLR is one of the consequences of climate change and global warming. According to the IPCC, the predicted sea level rise is expected to reach 1 m until 2100, and due to this, the coastal areas are high hazards. SLR is the annual mean increase or decrease water elevation measured at tide gauge stations, and the sea-level change along the Mediterranean coast is the sum of eustatic, glaciohydro-isostatic and tectonic factors [74].
In the Mediterranean Sea, the rate of sea level rise is 1.8 mm/year, as inferred from the analysis of tidal gauge data for the last century [76]. Galassi and Spada [77] extracted data from tide gauges placed along the Mediterranean shoreline (6 different cities) and observed SLR rates of a minimum of 0.64 to maximum 2.44 mm/year. In Greece, SLR due to eustatism is lower than 2 mm per year [78]. Based on these researches, the value of 1 mm/year was determined and taken as the most representative for the entire coastline of the study area. Same SLR value used in the research of Karymbalis et al. [24] for the assessment of the Coastal Sensitivity Index (CSI) in the Gulf of Corinth, which was located next to the Gulf of Patras.

CVI Calculation and Ranking
The initial step before the calculation of CVI, is to define the range values of each parameter to rank them with a specific vulnerability value; for example, "high or very high vulnerability. The ranges for the vulnerability ranking in Table 3 were derived from Karymbalis et al. [24] and Tragaki et al. [50]. However, the geotechnical characterization and its vulnerability ranking in Table 3 were presented for first time and gave a mostly accurate quantification of the resistance to erosion for this variable. As t mentioned before, for the calculations of the Coastal Vulnerability Index, two different methodologies and equations were used: (a) the general equation (CVI) of Gornitz et al. [16] and (b) a new (CVI WF ) by adding a weighted factor into the geotechnical characterization variable. Comparisons between the two different methodologies (CVI and CVI WF ) were performed. The two equations for the calculation of CVI and CVI WF are presented below: where a: Geotechnical characterization, b: Significant mean height wave, c: Coastal slope, d: Shoreline evolution, e: Average tidal range and f : Sea Level Rise (SLR). The final stage before the production of vulnerability maps is the ranking of the CVI values in the vulnerability classes. In the bibliography, researchers used different classification methods for giving vulnerability ranks in CVI values, such as equal intervals [7,16,30,[79][80][81], natural breaks [22,24,25,50,82], percentiles [5,28,33,[83][84][85] and standard deviation [51,86]. On this specific research, the equal interval classification method was chosen for ranking the CVI and CVI WF values.

Geotechnical Characterization
The results for this variable show that more than half of the coastline of the Gulf of Patras is under a regime of very high vulnerability (52.29%), while high vulnerability covers 34.49% of the area (Figures 3 and 4). The results reveal that 52.29% of the coastline of the Gulf of Patras consists of materials that have very high vulnerability, i.e., materials of small grain sizes (sands with D 50 < 2.50 mm) and poor geotechnical characteristics, while 34.49% is highly vulnerable and consists of materials of larger grain sizes (gravel with D 50 > 2.50 mm). As moderately vulnerable characterizes 2.17% of the total coastline, while low and very low vulnerability characterizes 6.88% and 4.18% of the total coastline, respectively, consisting mainly of sandstone coasts with medium cliffs and rocky coasts with high cliffs. (1) * * * * * 6 (2) where a: Geotechnical characterization, b: Significant mean height wave, c: Coastal slope, d: Shoreline evolution, e: Average tidal range and f: Sea Level Rise (SLR). The final stage before the production of vulnerability maps is the ranking of the CVI values in the vulnerability classes. In the bibliography, researchers used different classification methods for giving vulnerability ranks in CVI values, such as equal intervals [7,16,30,[79][80][81], natural breaks [22,24,25,50,82], percentiles [5,28,33,[83][84][85] and standard deviation [51,86]. On this specific research, the equal interval classification method was chosen for ranking the CVI and CVIWF values.

Geotechnical Characterization
The results for this variable show that more than half of the coastline of the Gulf of Patras is under a regime of very high vulnerability (52.29%), while high vulnerability covers 34.49% of the area (Figures 3 and 4). The results reveal that 52.29% of the coastline of the Gulf of Patras consists of materials that have very high vulnerability, i.e., materials of small grain sizes (sands with D50 < 2.50 mm) and poor geotechnical characteristics, while 34.49% is highly vulnerable and consists of materials of larger grain sizes (gravel with D50 > 2.50 mm). As moderately vulnerable characterizes 2.17% of the total coastline, while low and very low vulnerability characterizes 6.88% and 4.18% of the total coastline, respectively, consisting mainly of sandstone coasts with medium cliffs and rocky coasts with high cliffs.

Significant Mean Wave Height
The survey area was divided into zones, and each of them received a specific value using the average value of the significant wave height occurring 40 m from the shoreline. Generally, the wave height in the Gulf of Patras is low, since it is a narrow gulf, where a high wind speed cannot be developed and generate waves with high heights. The results of the numerical simulation of wave propagation for wind speed with a return period of 10 years for each wind direction in the eight zones revealed that the maximum value of the wave heights was 0.53 m, whereas the minimum was 0.23 m. Similar values (0.3-0.4 m) of wave heights were extracted from the wave buoy, which was installed in the port of Patras. According to Table 3 (vulnerability ranks), values of wave heights lower than 0.3 m and values between 0.3 m and 0.6 m are ranked with a very low and low vulnerability, respectively. Based on this ranking, the 50.20% of the area is characterized with a very low vulnerability, while 49.80% with low vulnerability (Figures 3 and 4). This is logical, since the gulf of Patras is not a wide sea to create higher waves.

Significant Mean Wave Height
The survey area was divided into zones, and each of them received a specific value using the average value of the significant wave height occurring 40 m from the shoreline. Generally, the wave height in the Gulf of Patras is low, since it is a narrow gulf, where a high wind speed cannot be developed and generate waves with high heights. The results of the numerical simulation of wave propagation for wind speed with a return period of 10 years for each wind direction in the eight zones revealed that the maximum value of the wave heights was 0.53 m, whereas the minimum was 0.23 m. Similar values (0.3-0.4 m) of wave heights were extracted from the wave buoy, which was installed in the port of Patras. According to Table 3 (vulnerability ranks), values of wave heights lower than 0.3 m and values between 0.3 m and 0.6 m are ranked with a very low and low vulnerability, respectively. Based on this ranking, the 50.20% of the area is characterized with a very low vulnerability, while 49.80% with low vulnerability (Figures 3 and 4). This is logical, since the gulf of Patras is not a wide sea to create higher waves.

Coastal Slope
Coastal slope values lower than 3% are characterized with very high vulnerability, and areas with values higher than 12% are characterized with very low vulnerability (Table 3). Values between 3% and 6% are characterized with high vulnerability, while values between 6-9% and 9-12% are characterized with moderate and low vulnerability, respectively ( Table 3). The results from the calculation of the coastal slope variable show that most of the coastline (42.57%) is very high in vulnerability (gentle slopes), while 33.19% is described as very low in vulnerability (high slopes and rock masses). The remaining shoreline was calculated at 12.13%, 6.47% and 5.64% and was described as having high, medium and low vulnerabilities, respectively. The spatial variation of the coastal slope parameter along the study area is presented in the next paragraphs (Figures 3 and 4).

Shoreline Evolution
Shoreline evolution was estimated with the use of DSAS by computing a total of 773 transects along a land baseline in the eight zones of Figure 2  In the western of zone 1, accretion with a rate over 1.50 m/y took place, whereas the majority of the transects showed a deposition (green color) with a rate that ranged from 0.50 to 1.00 m/y. Only in a small part of the eastern coast was erosion observed (orange color) with a rate that ranged from −1.00 to −1.50 m/y (Figure 9). Zones 2, 3 and 4 show erosion with rates of 81.67%, 80.00% and 70.18% of the relative transects, corresponding to 2.47 km, 5.85 km and 3.11 Km of littoral length. Erosion procedure dominated in zones 2 and 3 with a rate range from −1.00 to −1.50 m/y and from −0.50 to −1.00 m/y throughout almost the entire coast, while there were some places with severe erosion up to −1.50 m in both coasts (red color transects) ( Figure 9). The similar rates of erosion were observed in zone 4 too, while there was a small area where accretion prevailed with a rate that ranges from 1.00 to 1.50 m/y ( Figure 9). Furthermore, in zones 5, 6, 7 and 8, it seems erosion rates are lower instead. Moreover, in these zones, there are parts where transects record values of EPR between −0.10 and +0.10 m/y (blue color) with significant rates such as 28.05%, 30.59%, 30.77% and 20.54%, which correspond to a distance of 1.21, 1.33, 1.26 and 1.37 km, respectively ( Figure 10). These rates could be revealed stably, as they are very close to deviations of the statistical error. The highest shoreline change (erosion-accretion) according to the EPR rates was recorded as −2.57 m/year adjacent to zone 5 and +5.95 m/year in zone 1. Finally, although the overall trend of the shoreline according to the EPR rates is in agreement with those revealed from each zone, in zone 8, the shoreline trend recorded as accretion according to the mean EPR rates (+0.24 m/y) ( Figure 10).
According to the ranges of vulnerability ranking for this variable (Table 3), the results revealed that the biggest part of the coastline (65.27%) has a moderate vulnerability, while a significant percentage (19.57%) has been assessed with a high vulnerability. The other vulnerability classes were estimated with a small percentage (very low 3.58%, low 9.09% and very high 2.49%) (Figures 3 and 4).

Average Tidal Range and SLR
The value of the tidal range for the entire coastline is equal to 0.24 m and is characterized with a low vulnerability (Figures 3 and 4). The value of the SLR is equal to 1 mm/year and is characterized with a very low vulnerability (Figures 3 and 4). Both the tidal range and SLR variables do not show any spatial variances, since they remain the same along the entire coastline (Figures 3 and 4).

Comparison between the CVI and CVI WF Results
The percentage in each vulnerability class with the CVI application along the coastline of the research area was estimated as tidal range and SLR variables do not show any spatial variances, since they remain the same along the entire coastline (Figures 3 and 4).

Comparison between the CVI and CVIWF Results
The percentage in each vulnerability class with the CVI application along the coastline of the research area was estimated as  In general, with the insertion of a weighted factor in the geotechncial parameter, an increase was observed in the low vulnerability class (CVI Class 2) and the high vulnerability class (CVI Class 4). At the same time, a decrease was observed in the very low vulnerability class (CVI Class 1), the moderate vulnerability class (CVI Class 3) and the very high vulnerability class (CVI Class 5) ( Figure 6). In general, with the insertion of a weighted factor in the geotechncial parameter, an increase was observed in the low vulnerability class (CVI Class 2) and the high vulnerability class (CVI Class 4). At the same time, a decrease was observed in the very low vulnerability class (CVI Class 1), the moderate vulnerability class (CVI Class 3) and the very high vulnerability class (CVI Class 5) ( Figure 6). The results of the spatial analysis from the CVI model application revealed that the areas that have very low vulnerability at erosion are Zones 2 and 5, low vulnerability, Zones 3 and 6, moderate vulnerability, Zones 1 and 4 and, finally, high vulnerability, Zones 7 and 8. On the other hand, the results of the spatial analysis from the CVIWF model application revealed that the areas that have very low vulnerability at erosions are Zone 3; low vulnerability, Zones 2, 5 and 6 and, finally, high vulnerability, Zones 1, 4, 7 and 8.
The results from the application of both models in the entire area of study (Zones 1-8) as a percentage (%) and length in Km for each one of the vulnerability classes are presented in Tables 4 and 5. A column chart of the vulnerability classes in Km for Zones 1-8 of the research area with the CVI and CVIWF applications is also presented in Figure 7. Table 4. CVI and CVIWF results in Zones 1-4 of the research area.  The results of the spatial analysis from the CVI model application revealed that the areas that have very low vulnerability at erosion are Zones 2 and 5, low vulnerability, Zones 3 and 6, moderate vulnerability, Zones 1 and 4 and, finally, high vulnerability, Zones 7 and 8. On the other hand, the results of the spatial analysis from the CVI WF model application revealed that the areas that have very low vulnerability at erosions are Zone 3; low vulnerability, Zones 2, 5 and 6 and, finally, high vulnerability, Zones 1, 4, 7 and 8.

CVI Zone 1(km) Zone 1(%) Zone 2(Km) Zone 2(%) Zone 3(Km) Zone 3(%) Zone 4(km) Zone 4(%)
The results from the application of both models in the entire area of study (Zones 1-8) as a percentage (%) and length in Km for each one of the vulnerability classes are presented in Tables 4 and 5. A column chart of the vulnerability classes in Km for Zones 1-8 of the research area with the CVI and CVI WF applications is also presented in Figure 7.       (Table 6).

Analysis on the Weighted Geotechnical Factor Applied on CVI WF Equation
As it has been mentioned before, a weighted geotechnical factor was added in the CVI WF equation in order to examine the importance of this parameter in the assessment of coastal vulnerability. The comparison between the CVI and CVI WF results (Figures 7 and 8) showed that the areas mostly affected by the weighted factor are (a) the Papa Lagoon-Karnari (Zone 1) and Alikes-Gialos (Zone 4), since they were assessed from a moderate vulnerability with the CVI to a high vulnerability with the CVI WF in their biggest parts, (b) the Karnari-Ioniki Akti (Zone 2) and Peiros-West Kaminia (Zone 5), since they were assessed from a very low vulnerability with the CVI to a low vulnerability with CVI WF in their biggest parts and (c) Ioniki Akti-West Alikes (Zone 3), since it was assessed from a low vulnerability with the CVI to a very low vulnerability with the CVI WF in its biggest part. The remaining Zones 6, 7 and 8 (from West kaminia to Glafkos) remained in their biggest parts with the same vulnerability (low in Zone 6 and high in Zones 7 and 8), regardless of the CVI or CVI WF applications. The spatial distribution of the CVI and CVI WF results alongside the coastline of the gulf of Patras is presented in detail in Figure 8. The geotechnical weighted factor decreased the vulnerability class on some segments of the coastline, which were mapped as cobble and gravel beaches (D50 > 2.50 mm) from high to moderate. On the contrary, an increase of the vulnerability class was observed on some segments of the coastline that were mapped as sandy beaches (D50 < 2.50 mm) from moderate to high. The addition of a weighted geotechnical factor in the CVIWF equation By comparing the results from both equations, it seems that the CVI WF equation is more realistic for this specific area, since it characterizes the eight zones more accurately in terms of their vulnerability at erosion. Especially, it characterizes more accurately the western part of the gulf (Zones 1, 2, 3 and 4), since this part is mostly affected by the SE winds; thus, the coastal erodibility and vulnerability should be higher, and this is evident with the CVI WF results.
The geotechnical weighted factor decreased the vulnerability class on some segments of the coastline, which were mapped as cobble and gravel beaches (D 50 > 2.50 mm) from high to moderate. On the contrary, an increase of the vulnerability class was observed on some segments of the coastline that were mapped as sandy beaches (D 50 < 2.50 mm) from moderate to high. The addition of a weighted geotechnical factor in the CVI WF equation changed the coastal vulnerability to a lower class in gravel and cobble beaches and a higher class in sandy beaches, which is more accurate for this particular area.

Analysis on the Computed EPR Values and the Respective CVI WF Results
In Zone 1, the EPR shows that most of the coastline is under accretion, despite the fact of the existence of a geotechnical unit (fine sand) and a gentle coastal slope of very high vulnerability (CVI rank 5). However, in the eastern part of Zone 1 (Karnari Beach), erosion is observed with a rate ranging from −1.00 to −1.50 m/y, and this agrees mostly with the CVI WF results (Figure 9). In Zone 2, the EPR reveals that most of the coastline is under erosion, with a rate ranging from −1.00 to −1.50 m/y. The CVI and CVI WF results show this area as lowly vulnerable in its western part and moderate-to-highly vulnerable in its eastern part. Those differences between the CVI WF and EPR is due to the steeper coastal slope and the rocky coast existing in the western part of Zone 2 (Figure 9). A lower erosion rate with an EPR −0.40 to −0.80 m/y exists in Zone 3, especially in the western and eastern parts, which decrease at −0.30 to 0.00 m/y in the central part. Both the CVI and CVI WF results reveal high and very high vulnerabilities in the western and eastern parts of the zone, respectively, being in agreement with their respective EPR, which becomes a low vulnerability in the central part of the zone due to the existence of a higher coastal slope with soft rocks/hard soils (CVI rank 1) (Figure 9). In Zone 4, erosion dominates, apart from a very small accretion area in its central part. EPR is very high in the eastern part (−1.00 to −1.50 and > −1.50 m/y), where a dam was constructed upstream; thus, it blocks the sediment movement towards the coastline, and this agrees mostly with the CVI WF results ( Figure 9). Zones 5, 6, 7 and 8 have smaller rates of erosion compared with Zones 1, 2, 3 and 4, with some parts characterized as stable shorelines (−0.10 to 0.10 m/y). Specifically, in Zone 5, its western part has a high rate of EPR (−1.00-−1.50 m/y), while most of the remaining area has smaller (−0.50-−0.00 m/y) and stabler rates (−0.10-0.10 m/y). Both the CVI and CVIWF results revealed a low vulnerability because of the existence of a high coastal slope (CVI rank 1 and 2) with stiff-hard cohesive soils and medium cliffs (CVI rank 3). Zone 6 EPR has generally moderate rates and, in some parts, is stable (−0.10 to 0.10 m/y), while the results from both the CVI and CVIWF reveal a low-to-moderate vulnerability, which, in general, agrees with the EPR values ( Figure 10). The EPR in Zone 7 is similar with Zone 6, with most of its parts having moderate and high rates of erosion and some parts a stable rate of erosion. A comparison between the EPR and CVIWF shows that most of the shoreline is under moderate/high erosion and vulnerability and generally match, apart from a Zones 5, 6, 7 and 8 have smaller rates of erosion compared with Zones 1, 2, 3 and 4, with some parts characterized as stable shorelines (−0.10 to 0.10 m/y). Specifically, in Zone 5, its western part has a high rate of EPR (−1.00-−1.50 m/y), while most of the remaining area has smaller (−0.50-−0.00 m/y) and stabler rates (−0.10-0.10 m/y). Both the CVI and CVI WF results revealed a low vulnerability because of the existence of a high coastal slope (CVI rank 1 and 2) with stiff-hard cohesive soils and medium cliffs (CVI rank 3). Zone 6 EPR has generally moderate rates and, in some parts, is stable (−0.10 to 0.10 m/y), while the results from both the CVI and CVI WF reveal a low-to-moderate vulnerability, which, in general, agrees with the EPR values ( Figure 10). The EPR in Zone 7 is similar with Zone 6, with most of its parts having moderate and high rates of erosion and some parts a stable rate of erosion. A comparison between the EPR and CVI WF shows that most of the shoreline is under moderate/high erosion and vulnerability and generally match, apart from a few segments, which were calculated with stable rates of the EPR (−0.10 to 0.10 m/y) ( Figure 10). In Zone 8, the biggest part of the coastline has low EPR values, but in the eastern part of the zone is observed accretion. The CVI WF values agree only in some parts of this zone with the EPR, since the EPR rates show a low or stable erosional trend, and the vulnerability is calculated as moderate to high ( Figure 10). few segments, which were calculated with stable rates of the EPR (−0.10 to 0.10 m/y) (Figure 10). In Zone 8, the biggest part of the coastline has low EPR values, but in the eastern part of the zone is observed accretion. The CVIWF values agree only in some parts of this zone with the EPR, since the EPR rates show a low or stable erosional trend, and the vulnerability is calculated as moderate to high ( Figure 10). Weighted and unweighted vulnerability models have been tested for their appliance in the study area. For this purpose, the CVI and CVIWF models were compared with the calculated EPR values in each coastal segment. The main conclusion is that the historical evolution rate along the shoreline better suits the CVIWF model. As a result, the second estimation (with CVIWF) seems to be more realistic, since it coordinates better with the EPR and historical evolution rate of the shoreline.
Specifically, in Zone 1, the eastern part (Karnari), the EPR values are in the biggest agreement with the CVIWF rather than the CVI results ( Figure 9). In Zone 2, the EPR and Weighted and unweighted vulnerability models have been tested for their appliance in the study area. For this purpose, the CVI and CVI WF models were compared with the calculated EPR values in each coastal segment. The main conclusion is that the historical evolution rate along the shoreline better suits the CVI WF model. As a result, the second estimation (with CVI WF ) seems to be more realistic, since it coordinates better with the EPR and historical evolution rate of the shoreline.
Specifically, in Zone 1, the eastern part (Karnari), the EPR values are in the biggest agreement with the CVI WF rather than the CVI results ( Figure 9). In Zone 2, the EPR and CVI WF are closer to their results, especially in the central part of the zone (Figure 9). Zones 3 and 4 have a similar behavior between the EPR, CVI and CVI WF results, apart from a few segments where the EPR and CVI WF have better agreement (Figure 9). In Zones 5 and 6, the results of the CVI and CVI WF are similar with the EPR values, except for a few segments where the CVI has a better agreement with the EPR (Figure 10). In Zone 7, the EPR values are in the biggest agreement with the CVI WF rather than CVI (Figure 10), while in Zone 8, the EPR values in some parts are closer with the CVI and, in other parts, are closer to the CVI WF ( Figure 10).

Discussion
Among a variety of other indicators and methodologies, the Coastal Vulnerability Index (CVI) was chosen as a method for estimation and calculation of the vulnerability against sea erosion in the pilot area of the Gulf of Patras in Western Greece, due to its multiple advantages compared with other techniques. Another reason is that it is a very productive, low-cost and efficient tool for coastal hazard mapping and management. Thus, it gives a quick identification of erosion hotspots in order to enhance coastal engineers and local authorities to facilitate suitable management and coastal protection plans, the regulation of coastal zones and adaption strategies against erosion. However, the CVI approach reveals many uncertainties and doubts and, consequently, can bring wrong results in the vulnerability analysis. These uncertainties must be defined, improved and modified before any vulnerability calculations in order to increase the accuracy and reliability of the results. In the present study, to overcome the aforementioned problems, the following steps were conducted: (1) the use of geotechnical data instead of geologicalgeomorphological data, (2) the use of high-resolution remote-sensing data and (3) the insertion of a geotechnical weighted factor in the CVI formula. Finally, for the better understanding and evaluation of erosion analysis and vulnerability detection, as well as the efficiency of the CVI WF on these calculations, comparisons between the CVI and CVI WF results and with the EPR rates were performed.
A first attempt towards the improvement and quantification of resistance to erosion with CVI calculations was performed in this research with the use of the geotechnical parameter instead of the geological-geomorphological parameter. Instead of the classical qualitative methodologies using just the rock type or the geomorphology feature, the use of geotechnical data for the characterization of the coastal formations has proven to be a more reliable and accurate method to assess the resistance of a coastal formation to erosion. Only a few researchers have used geotechnical methods for the estimation of the CVI, such as Ruzic et al. [29], with the application of the Geological Strength Index (GSI) in order to calculate the quality of limestone rock masses in Croatia. Due to that, it is necessary to conduct more research on this variable, since it is maybe the most important variable that affects coastal erosion. The proposed modified CVI WF model is a first attempt to introduce geotechnical aspects in the erosion vulnerability detection, since only a small number of papers have been published on this topic. From this research, it can be stated that the geotechnical characterization for CVI calculations is mainly based on the physical and mechanical characteristics of soils derived from borehole drilling, insitu and laboratory tests, as well as sediment sampling. This type of characterization uses the soil's granulometry, geotechnical designation and the resistance to SPT measurements, which, in total, provide a more accurate measure of the soil's penetration resistance.
Another important variable that proved to be very important in the determination of the coastal vulnerability is the shoreline evolution, otherwise known as the historical shoreline change rate, with the application of extensive remote-sensing techniques. Generally, there are two main categories of studies regarding remote-sensing monitoring: those that they use free-of-charge low-resolution data and those that use high-resolution data. Over time, many studies have tried to process shoreline evolution with low-resolution satellite images, but their results that emerged always had a lack of accuracy. Different methods have been used to overpass this drawback regarding the accuracy of the measurements and the results they provide. Pardo Pasqual et al. [53] developed a very complicated methodology in order to ameliorate the Landsat data accuracy using high-resolution air photo mosaic by co-registering the low-resolution coastline derived from Landsat data to the respective high-resolution shoreline. In another similar study, Almonacid et al. [60] computed the average annual coastline movement of many Landsat images and yield a mean error of −4.7 m, showing that the shorelines derived from Landsat images are placed seaward in comparison to the respective high-precision shorelines. In a recent study in Northwestern Peloponnese, Apostolopoulos and Nikolakopoulos [57] followed a statistical analysis using the Digital Shoreline Analysis System (DSAS) software and concluded that low-resolution data are not suitable either for shoreline forecasting or shoreline mapping, as their accuracy were very low, as they presented a divergence that ranges between 6 and 11 m. In a similar study, Apostolopoulos et al. [58] concluded that the accuracy of the vector data derived from Sentinel-2 10-m spatial resolution images ranged between 1 and 5 m using the NDVI and NDWI indices, respectively. On the other hand, the synergy of high-resolution satellite imagery and high-scale aerial photos has become a very accurate combination in order to study the coastal zone changes. In the current study, the shoreline rate of change was computed based on high-resolution aerial and satellite images, improving the rate's accuracy. These datasets were processed in a GIS platform, and the statistical analysis was performed through DSAS application by computing the EPR. The use of high-resolution data provided an improved accuracy in the CVI analysis and gave better results in finding the respective erosion or accretion trends along the shoreline.
The original formula for the CVI computation did not include weighted factors, and this was mentioned as a big disadvantage of the method, since the parameters do not all have the same importance. However, not a big number of researchers put in their vulnerability analysis weighted variables. The calculation of weighted factors can be done either arbitrarily by the specialist who develops the model or with various semiquantitative methods, such as the Analytical Hierarchical process (AHP). In the current study, the weighted factor was added arbitrarily to the geotechnical variable only, since that specific variable seemed to be one of the most important in the CVI calculations. Initially, the CVI WF was calculated with a squared power (a 2 ) weighted factor of the geotechnical variable, but the results were similar with those of the unweighted factor index. Hence, the cubed weighted factor (a 3 ) was chosen in the CVI WF formula. By comparing the results from both equations CVI and CVI WF (a 3 ), it was concluded that the CVI WF is more realistic, since it characterized the eight zones of the investigation area more accurately in terms of their vulnerability at erosion. This was also analyzed in Section 4.8 (analysis on the computed EPR values and the respective CVI WF results). The same procedure was followed by Rao et al. [7], who added weighted parameters in the CVI formula in the variables of geomorphology, coastal slope and coastal change, while Ruzic et al. [29] used weighted factor only in the geological parameter. New studies have calculated the weighted factors by using more advanced techniques, such as the Analytical Hierarchical process (AHP) [30,31,[87][88][89][90]. In any case, the introduction of weighted factors in the CVI models is now necessary in order to calculate more accurately the coastal vulnerability; thus, the results can be derived by coastal engineers and managers for making the right decisions to mitigate the erosion phenomenon.
In order to compare which model CVI or CVI WF is more accurate and closer to reality, a direct comparison between the two vulnerability indexes (unweighted and weighted) and the EPR values was performed (Figures 9 and 10). Both methods predict the degree of vulnerability to erosion very well, but the second model CVI WF with the geotechnical weighted factor seemed to be more realistic, since it agreed with the rate of shoreline change in most of the coastal segments. This means that the geotechnical parameter affected the coastal vulnerability to a very high degree and especially in areas without high waves and SLR (such as the Gulf of Patras) and, therefore, may be one of the most important variables that controls vulnerability to erosion. Hence, the geotechnical properties and characterization of the coastal zones must be defined for the better and higher accuracy of the predicted vulnerability indexes and must be included in the CVI computations instead of the parameter of geology.
De Serio et al. [30] calculated the CVI with and without weighted factors and concluded that the unweighted CVI underestimated the vulnerability, while the weighted seemed to be more realistic, which agrees with the conclusions of this research. Bagdanaviciute et al. [88] also ran the CVI with weighted factors that led to an increase in the very high vulnerability class provided with this way higher accuracy and more consistent results compared with the unweighted CVI.
Similar studies based on indexes for the estimation of vulnerability to coastal erosion, to SLR, to flooding, etc. have been performed in many countries with reliable results. Toressan et al. [49] assessed the relative vulnerability of the Mediterranean shoreline of Egypt by using a Climate-improved Coastal Vulnerability Index (CCVI) evaluating the potential effects of climate change in the assessment of coastal hazards. Perez and Selvaraj [84] applied the same model in the district of Buenaventura in West Colombia using eight variables, while Hereher et al. [91] used the Coastal Vulnerability Index for the estimation of susceptibility in Oman with the implementation of five physical parameters. This index-based methodology was applied also in Spain, where a comparison study between different CVI approaches was performed in Barcelona [27], and in the Valdelagrana area, Rizzo et al. [92] determined the coastal susceptibility to erosion and flooding processes by means of an index-based method. Audere and Robin [93] presented an assessment of the vulnerability to erosion of sandy coasts of Vendee (Western France) using long-term measurements of the historical evolution of the shoreline between 1950 and 2016 and short-term measurements caused by extreme events. In China, Zhang et al. [94] developed a composite Vulnerability Index in a coastline of 5627 Km using biophysical, sensitivity and adaptivity variables.
The extensive descriptions in the literature of these models prove that vulnerability models with the combined action of an index and Geographical Information Systems (GIS) are among the most widely used methodologies for the assessment of coastal vulnerability. Their use in many areas and their reference by many researchers shows that it is a reliable method, but it also needs many improvements that need to be clarified, some of which were attempted to be clarified in this work.

Conclusions
The aim of the present study was the estimation of Coastal Vulnerability Index using only physical factors in the pilot area of the Gulf of Patras in Western Greece in order to give prominence to the importance of geotechnical evaluation and shoreline evolution in CVI calculations. For this reason, the CVI index modified the replacement of the geological parameter with that of geotechnical characterizations. In addition, the modifications included the introduction of a weighted factor in the original CVI formula in the new geotechnical parameter and the use of remote-sensing data of high resolution, for better and more reliable results in the calculation of the historical shoreline changes with the EPR analysis. The data inserted in the Geographical Information System (GIS) and the CVI calculated for the entire coastline of the pilot area was divided into 50-m coastal segments.
According to the results, with the insertion of a weighted factor in the geotechnical parameter, an increase in CVI Class rates was observed in the low vulnerability class (from 20.47% to 25.91% in Class 2) and the high vulnerability class (from 29.39% to 36.48% in Class 4). At the same time, a decrease was observed in the very low vulnerability class (from 20.13% to 14.59% in Class 1), the moderate vulnerability class (from 24.56% to 20.04% in Class 3) and the very high vulnerability class (from 5.45% to 2.98% in Class 5). By comparing the results from both equations, it seems that the CVI WF equation is more realistic for this specific area, since it characterizes the eight zones of the pilot area more accurately in terms of their vulnerability at erosion. Especially, it characterizes more accurately the western part of the gulf (Zones 1, 2, 3 and 4), since this part is mostly affected by the SE winds; thus, the coastal erodibility and vulnerability were found to be higher in those zones with the performed CVI WF calculations.
Furthermore, EPR calculated rates revealed that most of the coastline in the study area is under erosion, with mean rates that range from −0.27 to −0.90 m/y, but in some parts, higher erosion rates were computed. Zone 1 is the only area that accretion processes dominate against erosion, whereas the erosion process dominates in all the other zones (Zones 2-8). In general, the results of the CVI WF model are identical to the erosion hotspot calculated with the EPR values apart from a few parts where the variable of the coastal slope is high and, as a result, decreases the vulnerability. A comparison between the EPR values and CVI-CVI WF models revealed that the second model of the unweighted method is closer and more realistic for a vulnerability assessment. This fact shows that the variable of geotechnical characterization is important and must be included in vulnerability calculations. However, geotechnical data are not easy to obtain, since they demand experienced personnel and in-situ and lab equipment. In addition, it is difficult to evaluate geotechnical data that have been performed along a very long shoreline for CVI calculations because (a) it is a very complex issue and (b) there are not many papers in the bibliography with the integration of the geotechnical parameter in CVI models. Therefore, on this research, a first attempt towards this direction was made in order to demonstrate the importance of the geotechnical parameter on coastal erosion processes and vulnerability. The next stage of this methodology could be the addition of geotechnical data with higher accuracy in order to develop a holistic approach of the geotechnical aspect in CVI models.
Although the pilot area's coastline is largely composed of materials prone to erosion and with gentle slopes, it was not assessed with a high percentage of high and very high vulnerability on the entire coastline, because the Gulf of Patras is a micro-tidal environment with low wave heights and a small sea level rise rate. However, in a possible change of those hydrodynamic parameters, due to the climate change, the regime of vulnerability will increase dramatically. Due to that, it is significant and necessary to systematically monitor this phenomenon and record in detail the observed values from all the parameters.