The DPSIR Approach for Coastal Risk Assessment under Climate Change at Regional Scale: The Case of Apulian Coast (Italy)

In the present paper, the coastal risk is assessed, at a regional scale, to produce a risk map that will help to prioritize policies and economic resources in coastal management and planning activities. The DPSIR method is here used taking into account processes and cause/effect relationship between Drivers and Pressures, which induce an alteration of actual State and, hence, Impacts on the environment, society and economy. The study area is located in South Italy (Apulia region), where the Risk Index is calculated and mapped for all municipalities facing the Adriatic and Ionian Seas. Both coastal Vulnerability and Exposure Indexes are firstly calculated according to the procedure suggested in the EUrosion project (EU model) as the product of specific indicators describing the state of coastal zones, their natural characteristics and both natural and anthropic pressures. Based on both EU model results and knowledge of states and pressures of the study area, a new modified model is then proposed (Mod.E.M.) and final risk maps compared. The comparison shows that new Vulnerability and Exposure indexes better describe the ongoing coastal processes and pressures and allow us to identify hot-spot sites where more detailed analyses could be further focused on. The Mod.E.M. has more than just a local significance since the case study includes coastal areas with so different characteristics, that it can be easily applied to other coastal regions. Moreover, to take into account climate change effects, Risk Index is evaluated under Representative Concentration Pathways RCP4.5 and RCP8.5, mainly affecting the sea level rise and the storm surge level.


Introduction
Over the last decades, coastal risk assessment has become the first step in all strategic planning and management activities. Coastal areas are particularly exposed to risk drivers and are subject to a gradual weakening of the environmental asset, affecting also economic development and coastal community lifestyle (e.g., [1,2]). The growing urbanization and the massive exploitation are also threatening these areas, by reducing the resilience of territories after hazardous extreme events, exacerbated by climate change [3]. Moreover, the wide extension of coastal areas that must be protected requires economic investments that are often unsustainable for public administrations [4].
Hence, to support decision-makers in managing and prioritizing the interventions, the basic principles of the Integrated Coastal Zone Management (ICZM) have to be followed through the definition of a coastal Risk Index.
The analysis of critical environmental or safety situations is not only strictly related to the recognition of the physical processes driving erosion and flooding, but it depends primarily on characteristics of both natural and human environments and their interactions [5]. To create effective sustainable development policies, a correct balance between human, environmental, infrastructural, social and economic factors must be ensured, according to a cause-conditions-effect sequence, enabling a multidisciplinary and integrated vision of the overall environmental processes.
In the past years, several research projects aimed to assess coastal risk have been carried out and numerous approaches and GIS-based tools have been proposed, namely the DIVA model (DINAS-COAST project) [6] and the RVA method [7]. More recently, the THESEUS Decision Support System (DSS) [8], starting from an interdisciplinary approach, provided a coastal system status assessment exploiting an open-access approach. Moreover, the RISC-KIT project [9] proposed a two-step methodological approach in which a large scale Coastal Risk Assessment Framework (CRAF) is applied to detect coastal hot-spots where further in-depth higher resolution analyses are required.
Most of the risk assessment approaches are based on the Drivers-Pressures-State-Impact-Response (DPSIR) conceptual framework, that enables us to understand and explain the issues of actual or predicted situations and their potential evolution across different scales [10]. Indeed, through multivariate analysis, the DPSIR model can be used to examine mechanisms underlying environmental issues and interaction between society and environment [11]. The method is built on the selection of specific indicators, grouped by algorithms referring to a pyramid scheme of themes and objectives, aimed at framing the phenomenology under study.
A concise and reliable representation of processes underway can be only achieved if appropriate environmental indicators of the economic, social and environmental history of the study area and the connection patterns with potential hazards are known. Accordingly, over the years the DPSIR model has undergone significant transformations, mainly focused on deriving different definitions of States and Impacts, depending on the main aim of risk evaluation, available historical data and expert knowledge of processes and cause-effect relationships. Up to 2016, 25 versions of the DPSIR framework have been counted (e.g., [12][13][14][15][16]).
A significant research project funded by the European Union (EU) dealing with the risk assessment for the European coast based on DPSIR was introduced within the EUrosion project (www.eurosion. org). It follows a simplified DPSIR approach to map and classify coastal risk due to erosion and flooding at a regional scale, suggesting the use of specific parameters and indicators which are able to represent the state or a trend of coastal conditions.
In the present work, the DPSIR framework is used for assessing and mapping erosion and flooding risk in Apulia littoral (South Italy). Based on the analysis of both actual and past state of coastal territories, the risk profile is drawn for the study area to explain the relationship between causes and effects of the hazards affecting coasts. Consistently with EUrosion project guidelines, specific indicators are firstly used to quantify the Coastal Vulnerability Index (CVI) and the Coastal Exposure Index (CEI) and map the normalized Risk Index at a regional scale for all municipalities. Both CVI and CEI are then modified according to the risk profile. A Modified EUrosion Model (Mod.E.M.) is here proposed, involving different indicators and, for some indicators, integrating different variables. The main aim of the research is to develop a simple tool available for public administrations, useful to identify coastal hot spots at the regional scale and properly manage defense strategies priorities.
The paper is structured as follows. Section 2 recalls the models applied in the present work for quantifying the risk, according to the indicators and databases adopted. Section 3 briefly describes the study area, focusing the attention on natural and anthropogenic processes mainly responsible for the negative effects induced by flooding and erosion. According to the DPSIR framework, risk profile is defined for Drivers, Pressures, States and Impacts, specifically recognized in the Apulian coast. In Section 4 Vulnerability, Exposure and Risk maps derived from EUrosion and Mod.E.M. are compared and discussed. Conclusions on the adopted methodology for risk assessment at the regional scale and the risk level of Apulia coasts are finally drawn in Section 5.

EUrosion Model
The methodology proposed in the EUrosion project is one of the most popular models for coastal risk assessment widely used to map coastal hazard in the European littorals thanks to its simplicity and effectiveness [17,18].
Risk Index (RI) (Equation (1)) is estimated as a product of coast predisposition to erosion and/or flood, termed as Coastal Vulnerability Index (CV I) and the exposure of the surrounding environment termed as Coastal Exposure Index (CEI). Since risk levels depend on the interaction of the two indicators, regions with low exposure and high vulnerability are considered low-risk. It has to be highlighted that in this paper, according to the latest terminology suggested by [2], the terms "Coastal Sensitivity" and "Coastal Vulnerability", originally used in EUrosion project, have been replaced, respectively, by "Coastal Vulnerability" and "Coastal Exposure".
In Table 1 both Vulnerability and Exposure indicators used by the EUrosion approach are reported, together with the relative score ranges. CVI represents morphological, geological and sediment characteristics of the study area, as well as wave energy approaching the coast, and CEI represents the level and class of human activities, infrastructures and socio-economic development. CVI and CEI values are calculated by the simple addition of the scores assigned to the selected indicators, while various index-based approaches have been proposed to assess coastal vulnerability by estimating the geometric mean of some selected indicators [19,20]. Over time, the Gornitz's CVI, probably the most popular index [21] has been greatly modified including multi-hazard and exposure indicators [7,22]. Pressure factors associated with Sea Level Rise (SLR) and storm surge combined with wave setup during particularly severe marine weather events (termed as Highest Water Level (HWL)) are the main cause of coastal flooding. Sea level data can be extracted from EUrosion project data-sets that cover the period 1985-2001. HWL analysis can be successfully integrated by using sea level data recorded from national and local tidal stations that provide continuous and reliable datasets. Coastal evolution trend (TEV) globally represents an index of beach dynamics. It considers both positive (accretion) and negative (erosion) shoreline movements, deemed as pressure factors that increase coastal risk. Coastal strip geology (GEC) allows estimating coastal potential exposure to erosion reflecting potential instability, through the analysis of geomorphology, by defining erodable (sandy and mixed sediment shores) and non-erodable (rock and cliffs) coast. Another pressure factor which, according to EUrosion, is the expression of coast instability and, hence, could contribute to increased risk level is the percentage of "engineered frontage" of the coast, due to the presence of protection structures (ODC). The last indicator considered in CVI evaluation identifies the coastal municipal area (A RICE ), corresponding to the percentage of municipal areas in the Radium of Influence of Coastal Erosion (RICE). The RICE area includes all areas potentially subject to erosion and flooding over the next 100 years, located within 500 m from the coastline and that can be extended to areas lying under +5 m above sea level (a.s.l). Impact factors are quantified by estimating the Coastal Exposure Index in the RICE, representative of coastal natural and anthropogenic value. Urban and industrial areas in the RICE (U RICE ) and their growth within 10 km from the shoreline (U 10km ) have to be considered, as well as resident population (P RICE ) and high ecological value areas (E RICE ).

Modified EUrosion Model
A critical analysis of the previous indicators has led to the development of the Modified EUrosion Model (Mod.E.M.) to better fit the local detailed Path of Effects (PoE). The path analysis allows us to describe the connection among the involved variables and are helpful in developing effective management plans to reduce the coastal risk [23]. The need for a risk map according to different sea-level rise scenarios (Flood Directive 2007/60/EC) has been taken into account in the proposed changes to include climate change effects.
Since risk assessment at regional scale requires a detailed analysis of local phenomena and site-specific metrics [24], the analysis is based on a preliminary study of causes leading erosion and flooding and the related effects on coastal zones, by means of the construction of the local PoE.
Hence, Coastal Vulnerability and Exposure Indexes are amended. The new indicators considered in the Mod.E.M. are presented in Equations (2) and (3), where: (i) starred indexes have the same meaning as in the EUrosion model but refer to different databases; (ii) new indicators are highlighted in bold; (iii) neglected indicators are strikethrough and (iv) normal font relates to the unchanged indicators.
Climate-change-related flooding and erosion impacts are here considered by means of sea-level rise, not only due to stationary contributions of tides, waves and storm surges, but also considering extreme sea levels, thus affecting both SLR and HWL indicators. During the 21st century, according to the fifth Assessment Report (AR5) of Intergovernmental Panel on Climate Change [2], due to the climate changes affecting waves, storm surges and tides, the global mean sea level rise with medium confidence is expected to exceed the rate of the past 30 years with a rate of 8 ÷ 16 mm/yr, depending on the different long-term Representative Concentration Pathway (RCP) adopted scenario (0. 40  To take into account even different scenarios, sea level rise projections (termed as RSLs inLISCoAsT database) are considered [25], under RCP4.5 and RCP8.5 scenarios, corresponding to a likely global mean temperature increase of 2.0-3.6°C and 3.2-5.4°C in 2081-2100 above the 1850-1900 levels. The LIScoast database also reports 100-year Extreme Storm Surge Levels (termed as ESSLs) estimated by means of the combined effect of tides and water level fluctuations due to waves and storm surges [26].
As reported in Equation (2), the shoreline evolution trend expressed by TEV indicator in EUrosion model is replaced by the new Eroded Shoreline indicator (ESRL) which only concerns with landward shoreline movements. It is evaluated by the weighted average of Past Shoreline Erosion (PSE) and Recent Shoreline Erosion (RSE) as follows: where RSE and PSE are estimated from the eroded shoreline ratio assessed in the last 15 years and the previous decades, respectively. The ranking scores of RSE and PSE are reported in Table 2.
The analysis of PSE together with the RSE, helps to understand the long-term shoreline evolution trend and the most recent evolution. Analyzing only RSE could lead to a wrong interpretation of the process affecting the coast. Low RSE values indicate that negligible shoreline movements occur in the analyzed period. However, such an assumption could lead to misunderstanding the ongoing processes, since it could reveal a stability condition of the beach or, on the contrary, the full inhibition of the shoreline movement if the beach is completely eroded. Accordingly, the PSE helps us to understand the actual beach state, correlating it with its historical background.
Particular attention has to be paid in coastal erosion trend assessment, considering a temporal extent long enough to capture long-shore sediment dynamics neglecting extreme events and local dynamics. A threshold value for shoreline movement should also be applied to minimize errors resulting from operator digitization and exclude seasonal and/or tidal cycles from statistics.
Coastal areas flooding can be exacerbated when high river discharges and severe storm surges occur simultaneously since a storm can drive both high surges and river flows [27][28][29]. For this reason, to consider the further pressure induced by fluvial inundation on the coastal zone, the Hydraulic Hazard indicator (HH RICE ) is considered in the Mod.E.M. The new indicator can easily be estimated from the hydraulic hazard maps following the EC Flood Directive, calculated as the percentage of flooding-prone areas within the RICE.
The coastal protection structures indicator (ODC), considered in the EUrosion model, is neglected in the Mod.E.M. for regional risk assessment since a coastal structure could have different responses on the local scale, by increasing or reducing the sensitivity of the area against flooding and erosion. Such a morphodynamic response in presence of a coastal protection structure can be only quantified at the local scale, by means of numerical predictive calibrated models (e.g., [30][31][32]), able to reproduce hydrodynamics and morphodynamics induced by the presence of protection structures or harbors. EUrosion model considers also an indicator of the decrease in sediment supply from rivers, that has not been included in the discussion due to the major problems in collecting reliable riverine sediment transport data. P * RICE indicator is quantified considering population data extracted from national census data-sets. In the present work, both tourists (seasonal population) and residents are considered since coastal areas host a large number of visitors that, in some municipalities, are more numerous than the residents.
The ecological value indicator E RICE can be adequately estimated exploiting CORINE Land Cover data integrated with Natura 2000 and Regional Parks and Reserves maps. In many cases Natura 2000 areas include High Nature Value Farmland (HNVF) [33], hence, involving in the risk assessment also the exposure of agricultural areas. Table 2 reports the new indicators considered in the present work and their relative weights. For the remaining indicators, the same scores and weights of EUrosion project are assumed.

Description and Analysis of the Study Area
The Apulia region is located in South Italy (East Mediterranean coast) and faces the Adriatic Sea in the North and East sides and the Ionian Sea in the South and West sides ( Figure 1). The regional coastline is about 975 km long, with a large diversity of elevation and geological and morphological features ( Figure 1). Sandy beaches and rocky coasts each represent about a third of the total, while high cliffs are 21% of the total coastline. Coastal cliffs are often characterized by caves, which worsen the effects of wave action, undermining cliffs stability and causing their collapse [34]. Shores of mixed sediments (including pebbly beaches) represent the 12% of the total shoreline. Wave climate analysis is based on data collected offshore Taranto and Tremiti Islands by the Apulia Region Meteomarine Network (http://93.51.158.171/web/simop/home) and offshore Monopoli by the Italian Wave Network (IWN) as reported in Figure 1. The distribution of the percentage frequency of wave direction and height in the selected locations ( Figure 2) shows that the wave climate is moderate, characterized by low-medium energy sea states [35][36][37]. While extreme wave storm events are not very frequent, Apulian coasts are subject to considerable phenomena of coastal erosion and flooding, which lead to negative consequences on both natural and anthropic environmental as well as human life. In the early twentieth century, the realization of coastal railway lines and road system and the massive reclamation of wetlands started the initial transformation of coastline with a relevant erosion of most sandy beaches. Then, in the second half, the growth of tourism increased the demand for land in coastal areas to be used for tourist facilities. Nowadays, the Apulia region enjoys a growing reputation both on the national and international tourism markets. It is identified as a high nature value destination and sea and protected natural areas constitute the main environmental resources attracting thousands of tourists every year. In this regard, the highest priority of regional policy is to enhance natural, cultural and landscape regional resources to increase its attractiveness, without deteriorating natural environment. It follows that a further increase in tourist flows must necessarily be associated with forward-looking management of coastal territories, based on protection and recovery strategies aimed at preserving landscape and environment, and on the improvement of infrastructures and tourism services.
The regional hydrographic network is rather poor and mainly endorheic. Highly permeable soils inhibit the sediment movement through surface runoff and sediment trapping by dams and river regulation works contributes to extensive losses of wetlands and deltas, reducing watershed sediment input to the coastal sand budget [38]. Moreover, low-lying or depressed area, which constitutes a significant portion of regional coastal territories, are often subjected to flooding during storm events even not particularly intense [29,39].
Coastline evolution has been also strongly influenced by the gradual disappearance of dunes belts, which act as a sediment reservoir and a protection system for the back dune areas [40]. All these pressures led to a progressive deterioration of the physical environment with negative effects on the economy [41].
The above conditions forced local administrations to carry out an extensive campaign of beach monitoring [42][43][44][45], protection and restoration interventions [46], testing also alternative strategies [47], since the coastal protection strategies often prove to be not adequate to counteract the erosive action of waves, the reduction of sediments supply and flooding events.
From an administrative point of view, the regional coastline falls within the planning jurisdiction of 68 municipalities, some with constrained access, and this is a measure of management and planning fragmentation. Hence, following the principles outlined in ICZM, in 2007 Apulia region drafted and implemented the Regional Coastal Plan (PRC) [48], which still constitutes a useful regulatory and technical operational tool intended to standardize coastal areas protection and maintenance. The plan aimed at defining the status of the Apulian coast as well as to identify the most critical and sensitive (from an environmental and historical point of view) coastal zones [49]. Significant efforts were focussed on collecting large geospatial datasets that are regularly updated (http://www.sit.puglia.it/).

Risk Profile
According to the above characterization of Apulian coasts, processes are summarized through the Path of Effects (PoE), defined as a cause-effects chain, where processes are categorized in Drivers, Pressures and the corresponding States and Impacts (e.g., [23]). Figure 3 shows the path of risk flow chart developed for the presented case of study. Since differences in terminology and definitions in the DPSIR model could lead to a misunderstanding and ambiguities in the interpretation of results (e.g., [12,16,50]), a declaration of the adopted terms is required. In the present work, a risk profile is drawn by starting from the identification of the main drivers of erosion and flooding hazards typically considered in coastal risk assessment. Two main groups of drivers are defined (considered relevant) for the hazard characterization. The first group (anthropogenic) involves basic human needs-related drivers occurring with social and economic development. Human activities and policies aimed at satisfying such needs lead human-induced pressures on the coastal environment. The second group includes physical drivers exerted by natural phenomena, combining sea storms and meteorological events (e.g., precipitation) with climate change. The latter, even if actually derived from long-term effects of human existence, are here recognized as exogenous drivers, as they increase the magnitude and frequency of coastal storms, extreme discharges and considered as the main cause of sea level rise. The relative pressures are strictly related to the processes responsible for coastal erosion and flooding, resulting from both anthropogenic and physical drivers. Urbanization, growing population and subsidence are the negative pressures induced by social and economic development of coastal zones and their heavy exploitation. Extreme surge levels, sea level rise and river floods are here considered as the pressures induced by physical drivers, exacerbated by the effects of climate changes. Consequently, changes in coastal environment states result in impacts on human health, ecosystems and economy. It has to be stressed that in the Figure 3 quantitative/qualitative changes of States are reported. Negative consequences are associated with alteration of sediment budget balance, losses of dune belts, reduction of the emerged beach extension, land cover transformation and losses of wetlands and deltas. Impacts, in turn, affect drivers and exacerbate the process. Hence, decisive responses are required to prevent negative impacts and preclude additional pressures. In such a way, potentially damaging events can be proactively managed by means of contrasting adaptation measures.

CVI indicators
In the following analyses, the extent of the RICE area has a minimum width of 500 m extended to coastal areas below 5 m a.s.l. to expand risk analysis to low-lying areas.
The analysis of the A RICE indicator shows that several coastal zones, mostly located in PU2 and PU7 (Figure 1), are prone to coastal flooding due to the natural conformation of the most inland territories, characterized by very low depths above the mean sea level. Some of those low-lying zones are flooded even during ordinary sea storms, giving some areas a typical lagoon appearance with considerable damage to both productive and touristic sites.
Compared to other coastal European countries, in the Apulia region, both SLR and HWL indexes are constant and relatively low (Table 3). According to [2] offshore the Apulian coast, compared to the period 1986-2005, the mean sea level is expected to rise for 2081-2100 about 35 cm (RPC2.6). HW L 100 , obtained by combining the projections of astronomical tides and the 100-year meteorological tide (storm surge and wave setup), has been extracted from LISCoast database [25]. Present-day HW L 100 values show a good agreement with other available datasets (e.g., EUrosion [11]) and values estimated carrying out the analysis of extreme sea levels and tidal residuals recorded by the Apulia Region Meteomarine Network and the Italian National Tide Gauge Network [39,51]. Nevertheless, due to the limitations of the spatial and temporal scale in the ocean modeling [26], local surge and wave prediction are required at smaller scale risk assessments.
According to to [26], the mean sea level is expected to increase around 45 cm on average for Apulia by 2100 under RCP4.5 and about 70 cm under RCP8.5 (Table 3). While HW L 100 shows a very slight increase by 2100 since, as observed by [26], rising total extreme sea levels are mainly due to sea-level rise because of the bigger effect of climate change on RSLR than on extreme meteorologic events. The HH RICE indicator is estimated exploiting the official flooding areas mapped in the Hydro-Geologic Safety Plan produced by the Apulian basin authority [52]. In the Apulian coastal region, high-risk river flooding areas can be observed in PU2 and PU7.
Passive and active remote sensing techniques are an essential element in the early identification of critical areas [40,49,53,54], since extensive in situ campaigns are too cost-and time-expensive in wide areas. Aerial images taken in 1992, 2005 and 2017 are used for the aim. A polygon overlay analysis is used to determine areas where shoreline deviation is greater than 10 m and shoreline evolution trend (TEV) is estimated. The TEV analysis shows significant erosive phenomena in the Apulian coastline: the evolution analysis in the period 2005-2017 (Table 4) reveals that about 34% of sandy coast (127 km) suffered from erosion processes. Further analyses carried out by exploiting orthophotos taken in 2008, 2010 and 2011 highlight that sandy coasts suffered a sharp and sudden deterioration between 2008 and 2010, leading to the erosion of about 31% of the total coastline.
The coastal erosion of the Apulian coastline is certainly not a recent phenomenon, since 130 km of sandy coastline receded over the period 1960-1992 ( Table 4). The last analysis is carried out superimposing the IGMI cartographic maps (scale 1:25,000) dated back to 1960 with the 1992 coastline and determining only areas where shoreline deviations were greater than 30 m to take in account geolocation and large-scale map accuracy. Whereas in the period between 1992 and 2005 erosive phenomena affected only the 16% of sandy coast (63 km). The erosive processes are expanding not only to previously eroded areas but new critical areas are emerging in the entire regional coastline. Moreover, it has to be highlighted that in many cases the trend evolution has been limited by the presence of coastal protection structures that prevent shoreline recession and irreversible coastline changes. Information on both coastal geomorphology (GEC) and coastal protection structures (ODC) is achieved by further detailing and updating up to 2017 the data used for the Apulian Regional Coastal Plan [48]. Coastal areas with a high percentage of engineered coastline are located in the Gulf of Manfredonia (PU2), where almost 55% of the coastal protection works were placed to counteract and mitigate severe erosion processes affecting the area. However, except the above-mentioned area and those municipalities that host important port activities, the regional coastline evolution seems to be little affected by coastal protection structures.

CEI Indicators
Impacts are then quantified by estimating the CEI representative of coastal natural and anthropogenic assets, derived from CORINE Land Cover (CLC) soil maps 1990, 2000, 2012 and 2018. Both land use analysis (U RICE ) and land-use change over 1990-2018 in a 10 km wide coastal area (U 10km ) show a significant and growing urbanization in the Apulian coastline with a consequent increased demographic concentration (P RICE ). The evaluation of U 10km over two different periods (1990-2018 and 2006-2018) shows that for most of the examined municipalities, the increasing trend has been reduced in the second decade, especially for those areas facing the Ionian Sea. On the other hand, an increase of urban/industrial areas involves some of the most important urban centers and towns with a strong tourist vocation. Over the last decades, due to the expansion of industrial/commercial areas, together with high population density, regional coastal areas have suffered a major loss of natural areas and habitats that serve as buffer zones between the sea and the hinterland, resulting in a coastal squeeze. P RICE indicator is estimated exploiting the data collected by the Italian National Statistical Institute (ISTAT) for years 1991, 2001 and 2011. Among the population data, the touristic population is also considered, since in the last years an increase of seasonal population occurred. E RICE indicator is derived from CORINE Land Cover used datasets integrated with Natura 2000 and Regional Parks and Reserves maps. In the Apulia region, the network of protected areas, established under the EU Birds and Habitats Directives, occupies 21% of the regional area. It should be also emphasized that a large percentage of Natura 2000 areas in Apulia region includes High Nature Value Farmland (HNVF) [33] allowing to include in the risk assessment also the exposure of agricultural areas. In Table 5 the databases used in the present work are reported.  CLC (1990CLC ( , 2000CLC ( , 2012CLC ( , 2018; Natura 2000 Network, Regional Parks and reserves U 10km CLC (1990CLC ( , 2000CLC ( , 2012CLC ( , 2018

CVI and CEI Indicators
The quantification of the CVI and CEI indicators is performed for all coastal regional municipalities. Both vulnerability and exposure indexes are represented through a color map (Figure 4) that easily allows identifying causes and impacts on the regional coastline at the municipality scale. With respect to Figure 4, it can be seen that the highest scores in GEC and ESRL indicators are present in PU1, PU2, PU3 and PU7, where the greater concentration of population (U RICE and P RICE can be seen in PU2 and PU3. The A RICE indicator reaches the high or medium score in almost all the municipalities, while the highest scores in HH RICE can be found in PU2 and PU7.
CVI and CEI under RCP4.5 and RCP 8.5 conditions are not reported in Figure 4, since only a variation in SW L score occurs. Climate change, according to [26], seems to affect only sea-level rise rather than storm surge levels.

Risk Levels
Risk levels have been mapped at a regional scale with levels normalized with respect the maximum value theoretically derived from the model, according to the EUrosion procedure ( Figure 5a). Accordingly, 5 classes of risk are identified: (i) very low, 0-5 (ii) low, 5-15 (iii) medium, 15-30 (iv) high, 30-50 (v) very high, 50-100.
Since the presence of coastal protection structures and harbors should not be considered a priori as a pressure variable, according to the discussion presented in Section 2, the indicator ODC has been neglected in the EUrosion model. Figure 5a does not highlight significant risk situations, reaching medium and low values, if compared to the risk levels observed in the most stressed coastal areas in Europe (e.g., [11]) and in the Northern Adriatic Sea in Italy (e.g., [7,17,55]). However, about 27% of the regional coastline presents a medium risk level. Two coastal hot spots are highlighted in the Gulf of Manfredonia (PU2), where the predominant risk factors are essentially due to the coastal type, characterized by long sandy coastlines and low inland heights, the erosion rates and the presence of high ecological value areas. Other critical areas can be detected in the regional coastline, especially the major cities and the most known touristic areas where coastal exposure has larger values. In Figure 5b the normalized risk indicator (nRI) derived from the Mod.E.M. is mapped for the regional coastline, using the present-day values. The comparison between the two models (Figure 5a,b) shows that both the new indicators and data-sets considered in the present work affect RI calculation. With reference to Mod.E.M. results, an increase of risk level can be observed, mainly due to the HH RICE indicator and seasonal tourist population considered in P RICE , since they increase both CVI and CEI. Almost 25% of the municipalities are affected by river flooding in the coastal area and these events cannot be overlooked since heavy rainfalls are increasingly driven by changing climate [2]. This aspect can be better appreciated by looking at Figure 6, where land cover, shoreline erosion and river flooding areas are mapped in a region within PU2. The highlighted area is chosen because it contains several important aspects: (i) it has been heavily affected by erosion since the second half of the 20th century and the shoreline retreat is currently ongoing, (ii) a large part of that area is covered by the largest saltworks in Europe and it is classified as a protected natural area and (iii) wide flood-prone areas have been identified by [52]. In such a context the EUrosion procedures highly underestimate the coastal vulnerability and exposure.  The inclusion of the seasonal population leads to an increase of the P RICE for 30% of the coastline, highlighting the high tourist pressures on the exposed population and coastal facilities. On the contrary, the new indicator ESRL which takes into account both past and recent shoreline erosion are reduced with respect to TEV for about 7% of the shoreline since beach accretion is not considered, whereas a small increase in ESRL has been estimated for about 7% of the shoreline.
The validation of risk assessment models remains a challenge in multi-hazard environments such as coastal areas [18,55], therefore estimated risk levels have been compared to publicly available information (e.g., Civil Protection bulletins and local news) and in-situ survey campaigns after severe weather events [29,40,56]. Obtained results demonstrate that the proposed Mod.E.M is able to reproduce the coastal situation more closely to reality, introducing new pressures that cannot be neglected in the study area, considered as a benchmark in the present work, and deleting indicators that cannot be included without a detailed local analysis.
The changes made on the indicators stem from a case study that covers so wide coastal regions, including different main Physiographic Units, that the use of new indicators can be easily extended in different territorial contexts. On the other hand, the suggested modifications don't increase the number of indicators since a good index for risk assessment should be easily applicable and avoid overlapping among the various indicators [7,22,57]. Figure 7 reports the risk map under the RCP4.5 scenario. Under RCP4.5 (Figure 7), the LISCoast projection datasets relative to sea level rise induce an expected increase of CVI overall coastal municipalities, since SLR assume higher values, reaching the maximum score, while HWL score remains low. Low-risk shoreline ratio decreases from 60% to 40% passing from present-day values to the RCP4.5 scenario, whereas the medium risk shoreline ratio increases from 26% to more than 38% (Figure 8). In the worst scenarios, about 14% of the regional shoreline becomes a "high risk" area. No changes between RCP4.5 and RCP8.5 scenario are observed since the expected sea-level rise by 2100 reaches the maximum score already in the RCP4.5 scenario and climate change does not affect HWL scores [26].

Conclusions
The DPSIR model is used here for defining a risk map at the municipal scale to identify coastal areas more exposed to erosion and flooding risk. In this paper, a Modified Eurosion Model (Mod.E.M.) is proposed to map coastal risk at a municipal scale and identify coastal hot-spots. It can be used for rating coastal regions in terms of exposure to coastal erosion and flooding and it is useful for the management and planning of coastal areas.
The Mod.E.M. is here used to map the Risk Index for Apulia coastal areas (South Italy). The new indicators have been proposed after a critical discussion of the original indicators based on the local PoE analysis. Compared to the original EUrosion model, it has to be highlighted the importance to include indicators that also consider past erosive processes, river flood hazards and seasonal population. Since the case study includes coastal areas with different characteristics, the improved model can be easily applied to other coastal regions and it has more than just a local significance. Moreover, climate change scenarios show an increase in the number of coastal hot-spots in the region requiring further analysis. The Mod.E.M. has the advantage of exploiting public data without the need for numerical modeling. Once the coastal hot spots are identified, more in-depth studies are needed to understand the occurring phenomena, by means of more detailed data (also not of public domain), in-situ surveys and monitoring programs, supported by suitable numerical as well as physical modelling.