The Effects of Urban and Economic Development on Coastal Zone Management

: The land transformation process in the last decades produced the urbanization growth in ﬂat and coastal areas all over the world. The combination of natural phenomena and human pressure is likely one of the main factors that enhance coastal dynamics. These factors lead to an increase in coastal risk (considered as the product of hazard, exposure, and vulnerability) also in view of future climate change scenarios. Although each of these factors has been intensively studied separately, a comprehensive analysis of the mutual relationship of these elements is an open task. Therefore, this work aims to assess the possible mutual interaction of land transformation and coastal management zones, studying the possible impact on local coastal communities. The idea is to merge the techniques coming from urban planning with data and methodology coming from the coastal engineering within the frame of a holistic approach. The main idea is to relate urban and land changes to coastal management. Then, the study aims to identify if stakeholders’ pressure motivated the deployment of rigid structures instead of shoreline variations related to energetic and sedimentary balances. The inﬂuence of coastal protection (described by an indicator taking into account the linear density of the rigid protection extension alongshore) measures have been considered. Finally, the economical performances in the period just after the land transformation have been studied by means of a dimensionless index based on the concept of the local unit. The method has been applied to the case study of the Abruzzo Region. Results reveal a large urbanization growth and a generalized over-protection of the coast with a resulting spatial alternation of shoreline retreat and advance. On the other hand, the analysis of the local economy reveals that the presence of tourism activities does not necessarily lead to an increase in the local economy. Indeed, the most important improvement in local economies in the region can be associated with coastal stretches characterized by a high environmental value.


Introduction
In the last 50 years, there was an increase in anthropic pressure on coastal areas e.g., [1][2][3][4]. As underlined by Mentaschi et al. [5], coastal environments host about 40% of the world's population as well as infrastructures, tourism, and social activities. Coastal municipalities in Italy are 645, covering a total area of about 43,000 km 2 , equivalent to 14% of the entire national territory with a population of about 17,200,000 inhabitants. As confirmed by Romano et al. [6], since the 1950's, the peninsula had seen a strong land transformation process. The increase of urbanization in coastal areas highlights a surprising match in per-capita urbanized areas among many Italian regions, in line with the western European Union average (360 m 2 /inhabitant), although they started off from very different situations in the 50s. As a consequence, this process triggered the threat of extreme form of a dispersion in the key areas of national landscapes and ecosystems (i.e., protected The phenomenon of beach dynamics is originally natural and contributed over the centuries to shape the coastal landscape. However, it is also true that the migration of humans towards the marine areas and strong modification of internal areas (i.e., with the deployment of artificial reservoirs and the protection of watersheds from erosion) induced unavoidable interference in the coastal environments e.g., [11]. Therefore, at present, it is commonly accepted that the phenomenon of coastal dynamics is affected by both natural and human factors. Problems are amplified due to the lack of sediments that often remain trapped in artificial reservoirs [12] and harbors, and in general, due to all the phenomena inducing a decrease in river sediments e.g., [13]. Moreover, large volumes of sediments are lost for the progressive dismantling of dunes because of anthropic pressures e.g., [14,15], as well as for vertical land movements (i.e., subsidence) due to natural phenomena (i.e., tectonics), anthropic activities (e.g., groundwater table drop) and sediment compaction that determines a decrease of sediment supply e.g., [16,17] The damages suffered by dune systems could reach a threshold that makes the consequences irreversible e.g., [18] reducing their fundamental role in coastal dynamics e.g., [1,19].
This situation must be considered also in view of the future possible scenarios of climate change. Indeed, one of the major climate change effects on the marine area is the rise of sea level (SLR). Combining the effect of the SLR with the increase of the frequency of strong storms, wind, and storm surge e.g., [20][21][22][23][24][25], in the near future there will be the increase of flooding e.g., [26] impacting on the hazard (and then on risk) of the urbanized or protected areas on the coast e.g., [21,27]. The Intergovernmental Panel on Climate Change (IPCC) stressed e.g., [21] that, even without the SLR, presently, there is evidence that confirms inter-annual changes in inundations e.g., [28]. The resulting interaction between land use planning and managing is a recent open question in the literature e.g., [29].
The controlling of coastal dynamics is often performed by "soft" and "hard" techniques e.g., [15,30,31]. One of the most employed soft solutions is beach nourishment, devoted to restoring the sediment budget re-introducing sand on the beach e.g., [32][33][34][35] (i.e., paying particular attention to geomorphological and sedimentological characterization e.g., [36][37][38]). Hard solutions are rigid structures parallel (emerged or submerged breakwaters e.g., [39][40][41]) or perpendicular (groins e.g., [42]) to the shoreline aiming at reducing the energy of the incoming waves in the first case and decreasing (totally or partially) the longshore transport [43] intercepting the sediment transport in the latter case e.g., [44]. For a more detailed description of these coastal processes, readers may refer to e.g., [43,45,46]. Despite the difference in aims and impacts, experience has shown that there are no correct or incorrect solutions, but often the best results have been reached by combining soft and hard solutions e.g., [47].
Although the large-scale sediment budget is not influenced by hard structures, they are effective in restoring beaches at their deployment sites. Nevertheless, boundary effects can influence the morphodynamics of downdrift areas e.g., [30]. Therefore, these effects must be taken into account with care at the design and monitoring stages.
On the other hand, nourishment has a beneficial impact on the beach uses but is a temporary solution that needs maintenance. There is no general rule to solve the problem, but a possible way to mitigate the effects comes from coastal management that implies a general analysis of the possible cause-and-effect relationship between urbanization, land-use, and coastal dynamics with attention paid to sustainability.
In this context, this paper aims to analyze how urbanization, and in general the urban development (used as a proxy of the economic growth), can influence coastal management and, therefore, the possible resulting executions of coastal interventions. In particular, attention is paid to how the economy has affected the coastal management policies, rather than the possible effects of coastal dynamics on the economy. The case study of the Abruzzo Region is chosen to highlight the applicability of the proposed method to a real case.
In particular, considering that the possible response of a natural system is a long-term process, the idea of the work is to refer to two different time spans. The evolution of urbanization e.g., [6] has been evaluated in the period 1950-2000 (i.e., the forcing factor) while the shoreline variation is referred to in the period 2001-2018 (i.e., the response to coastal management policies). The urbanization changes in the span of 2001-2018 are not considered, as it is assumed that the effects of coastal management are delayed with respect to the coastal management measures.
Moreover, the shoreline changes are analyzed taking into account the hard defense system built along the coast. A final analysis is devoted to understanding how the development of urbanization has impacted the economy of the municipality of the studied Region.
The paper is structured as follows. Section 2 describes the methodology for urban analysis and morphological evaluation. Section 3 illustrates the features of the case study and available data, Section 4 discusses the obtained results. Concluding remarks close the paper.

Methodology
As previously underlined, the great variation of urban cover in Italy and other countries were in the period ranging from the 50s to the end of the twentieth century. Starting from this evidence, the main idea of the work is to investigate the possible link between urban development (and a consequent change in land use) and coastal dynamics. The possible effects of human activities on the environment are often referred to a long-term scale, then the effects can be observed shifted in time. Then, urban variations and shoreline evolution are refer to different time spans. The first are evaluated considering the difference of the urbanization values in 1950 and in 2000 (hereinafter referred to as ∆ URB [m 2 ]), while the latter is referred to the shoreline extracted for the years 2001 and 2018. The state of urbanization of the Italian coast on a municipal basis and its evolution were determined through the use of data deriving from the study by Romano et al. [7]. This dataset resulted from elaborating data of the Italian National Institute for Environmental Protection and Research (ISPRA) [48]. Following the suggestions given by Syvitski et al. [12], the evaluation of the morphodynamics and urban information have been evaluated using the concept of Radius of Influence of Coastal Erosion (RICE) [10,47]. The influence areas (hereinafter referred to as IA) can be defined considering an area characterized by an extension up to 500 m landward [10,47]. The extension toward the sea is 2500 m, i.e., up to the typical width of morphodynamics active beach zone extending up to the depth of closure e.g., [49,50].
Considering only the sandy coasts (i.e., by neglecting the phenomena related to cliff erosion and coastal dynamics) the shoreline variation has been analyzed by using a diachronic analysis by considering the local total variation (m) of the shoreline positions between two different epochs (∆Y i,i+1 ) for each IA (see Figure 2a). The evaluation of ∆Y i,i+1 requires (for each IA) the definition of transversal direction to the shoreline and the geographical location of each i-th position. In more details, the mean location of the shoreline of each IA for each epoch has been evaluated. Consequently, the transversal direction is the orthogonal line to the mean value of the mean positions (see Figure 2b) between the considered epochs. It is then possible to evaluate the mean displacement of the shoreline in the given j-th IA as: where x is the coordinate in the longitudinal direction (locally parallel to the shoreline) and Y the transversal direction (locally perpendicular to the shoreline). Going deeper into details (see Figure 2b), within each area of influence, the direction along which shoreline displacements are evaluated must be defined (n). Then, for each IA, the direction (n j , being j the j-th area of influence) has been defined as the spatial mean of the local normal directions to the shoreline at epoch i at all the points belonging to the given IA. The local normal directions (n k ), being k the points belonging to the IA) have been estimated by employing a finite-difference first-order scheme with a centered stencil. Figure 2b shows a qualitative sketch representing the methodology described above. As far as the shoreline data are concerned, the data come from available datasets. Through these data, using whatever procedure, the shoreline position in two different epochs must be extracted (see Section 3.3 for details about the extraction of the shoreline for the case study).
Historically, at least at the national scale in Italy, the use of breakwaters is the most used practice to defend coastal areas against erosion e.g., [31]. More recently, groins are deployed along with breakwaters in order to avoid loss of sediment in particular after beach nourishment. Therefore the situation is a prevalent alternation of groins and breakwaters (emerged and submerged). In order to assess the influence of rigid protection (that locally influenced the nearshore hydrodynamic inducing variation in the whole management area) on the changes in the shoreline position, a "Protection Index" (hereinafter referred to as PI) is proposed herein. More in details, although a correct interpretation of the influence of the rigid structures on shoreline variation should be carried out using at least a one-line model e.g., [51][52][53], considering the aim of the work (not focused on the hydro-and morphodynamic details) the index is intended to relate marine structures to shoreline evolution from a phenomenological point of view.
The PI is a dimensionless measure of the length of the rigid protection alongshore (L structures ) against the total extension of the shoreline (L shoreline ) for each IA: This index is lower than one when the shoreline length in the IA is lower than the linear extension of the rigid protection and one when the two values are equal (i.e., the total extension of the shoreline in the IA is protected by a longitudinal breakwater).
Changes in the coastal area (both concerning shoreline variation and urban growth) suggest an investigation into the possible effect (direct or indirect) on the local economy. The analysis has been carried out using the Local Unit variable, as suggested by the European Union [54]. The Local Unit (hereinafter referred to as LU) is defined by the European Law as: "an enterprise or part thereof (e.g., a workshop, factory, warehouse, office, mine or depot) situated in a geographically identified place. At or from this place economic activity is carried out for which-save for certain exceptions-one or more persons work (even if only part-time) for one and the same enterprise" [54]. The specific variables of LUs are the localization, the typology (classified based on the Classification of Economic Activity, ATECO 2007 code) and the number of employees. Only the ATECO codes directly related to as "coastal resource" have been considered herein. This choice is motivated by the aim to consider the "coastal resource" as the main driver of the local scale economy. Therefore only the categories with a direct correlation with touristic activity have been considered i.e., hotels, accommodation services, house rental, car leasing, tour operators, etc.
These data are often available only at the municipality scale. Therefore, to adapt this information to the influence area scale, also in this case an Index of Local Units (hereinafter referred to as ILU) has been proposed: where ∆% is the percentage variation of the LU in a given time span and for each municipality. It has to be stressed that ILU values are assumed constant over the areas of the considered municipality. Table 1 resumes all the proposed indexes. Table 1. List of the proposed index.

PI
Dimensionless measure of the weight of the rigid protection extension both alongshore and cross-shore (L structures ) against the total extension of the shoreline (L shoreline ) for each IA

ILU
Percentage variation of the LU in a given time span and for each municipality. It is constant over the areas in which the municipality is divided

Site Description
The case study is a stretch of the Adriatic coast on the Abruzzo Region, in Central Italy. The coastal area includes 19 municipalities, with a total coastline length of about 130 km. From a morphological point of view, it is heterogeneous with features characterized by a low and sandy coast and others mainly high and rocky, which are not subjected to erosion processes.
The Adriatic coast, like the Tyrrhenian one, was the protagonist of compulsive urbanization that has increased about four times, assuming a current value of over 11%. A rate of increase that was almost similar also along the central-northern Tyrrhenian coast. Soil transformation rates were very pronounced: in the Adriatic coastal municipalities an urban conversion of over 0.04 km 2 /day can be estimated, but, considering the development of the linear front of the coast, an urbanization speed of almost 10 km/year from the postwar period to today can be estimated [55]. Starting from the north, the coastal area is characterized by a preponderance of sandy beaches with the presence of a protected marine area ("Torre del Cerrano"). There are productive and touristic activities that have encouraged an increase in human activity over the years. The presence of the Giulianova harbor (see Figure 3) and the A24 highway make this area one of the most developed of the entire region. At the center, there is another section characterized by low and sandy beaches. This area contains the metropolitan city of Pescara, the main economic hub of the entire Region. This central area ends with the channel port of the Pescara River. The harbor is connected to the highway and the near Abruzzo airport (the only one in the Region). Moving to the south, the Ortona harbor identifies a non-homogeneous area from a morphological point of view. South from Pescara, there is an area with low and sandy beaches, on the contrary near Ortona harbor the coast is mainly composed of high and rocky coasts. On the contrary to the northern and central parts of the Region, in this section, there is a massive presence of agricultural activities not inland from the coast. The remaining part of the coast is mainly rocky (with some long sandy beaches exceptions) and ends in correspondence of the Trigno's river mouth. Here, the regional protected area named "Riserva naturale Punta Aderci" stretches from the attractive beach of Punta Penna, near the Vasto harbor, up to the Sinello's river mouth (at the southern boundary of the municipality of Casalbordino).
The total length of the shoreline prone to be eroded is about 98 km, the 79.3% of the total extension of the shoreline. To counteract the erosion phenomena, during the last 60 years, a series of "hard" and "soft" interventions have been made (see Figure 4). Actually, the total length of the rigid structures (seawalls, groins, and breakwaters, both emerged and submerged) is about 60 km, 49% of the total coastal length (considering only the sandy beaches, about 60%).

Management Unit Definition
Looking at the sediment transport, due to the presence of the four ports, and considering the local economies, the Abruzzo Region could be split into four "management units" (hereinafter referred to as MU) depicted in Figure 5.
On the other hand, looking at the population variation it is possible to divide the Region into two different sectors: a northern and a southern part bounded by the metropolitan city of Pescara. This is mainly due to the different morphological features of the coast. The northern area is prone to be intensively urbanized due to the flat areas near the shoreline. As previously underlined the southern area, instead (south to the metropolitan area of Pescara), is characterized by an alternation of rocky and sandy beaches influencing urbanization and land use. It is possible to find these features also in the economic framework.
The area north of Pescara harbor contains infrastructures for touristic activities, companies, and is well connected with the northern area of the Italian peninsula. The situation is different in the south of the Region where the coast's features do not favor the exploitation of coastal touristic activities. Nevertheless, the morphological configuration of this area has retained the environment, by promoting the institution of areas characterized by high environmental value. Indeed, the Abruzzo Region has nine terrestrial protected areas, two in the northern part and seven in the southern area, confirming the role of urbanization in the progressive change of nature and environment. The areas recognized by the European Commission as a site of community importance (ZPSs) in Abruzzo are 6, all in the southern Area [56][57][58].

Data
The data concerning the urbanization evolution can be referred to the ISPRA database. However, the raw data have been analyzed by using the method proposed by Romano et al. [7]. The original study of [7] analyzed the land transformation for the Italian peninsula. Here, a subset of the analyzed data concerning only the Abruzzo Region data for the period 1950-2000 (consistently with the proposed parameter ∆ URB ) have been used. For a more detailed description, readers may refer to [7].
Following the main aim of this study (see Section 2), the shoreline has been extracted considering the time span 2001-2018 using rectified orthophotos (2001) and satellite images (2018). Both have been re-sampled with a dimension of 0.25 m for each pixel. Therefore, the resulting pixel dimension is 0.25 × 0.25 m. Based on the pixel dimensions and re-sampling the algorithm of the shoreline extraction a total displacement of 5 m (i.e., ±0.25 m/year) has been used to assess the stability of coastal stretches (i.e., neither in advance nor in retreat). The reason for the selection of these two epochs lies in the aim of the proposed indicators to catch the long-term evolution of the shoreline that is assumed to be delayed with respect to the coastal management measures. The automatic procedure was implemented based on the pixel-by-pixel analysis within the RGB space. More in detail, each image has been indexed by using the RGB information (i.e., lower than a given threshold identified with a sensitivity analysis) in order to discriminate the water and the land e.g., [59]. The interface has been then identified. Then, the raw results of the automatic procedure were manually checked, and if needed, amended to meet the real configuration of the shoreline as visualized in the images (either orthophoto or satellite image). Thereafter, the shoreline was further corrected by using the measured tidal level (as extracted from the time series collected at Ortona using the mareographic station belonging to National Tidal Gauges Network) and by estimating a mean emerged beach slope (as estimated based on the LIDAR survey performed by National Authorities). More specifically, the LIDAR data has been used to estimate the slope of the emerged beach. Then the slope has been utilized to correct the shoreline location by taking into account the synchronous tidal level collected in the Ortona station.
The same procedure has been implemented in order to extract all the protection structures insisting on the coast.
Finally, the economic data from the Italian National Institute of Statistics (ISTAT) have been used. The available period was 2011-2017 only at the municipal scale. From all the database, as previously described, only the ATECO 2007 codes referred to activities related to the coastal economy have been used.

Results and Discussion
The coastal stretch of the Abruzzo Region, following the methodology detailed in Section 2, waswas divided into 780 IAs (with an average alongshore extension of 200 m), with 628 IAs referring to the sandy coast. Results reveal that starting from 1950, there was a mean increase in the population of +128% while urbanization reaches +718%. In the northern part, the population variations were only positive, and the municipality of Montesilvano has the highest value of increment equal to +635% while the total variation of urbanization is equal to +1800%. The southern area, instead, was characterized by small variation (either positive or negative), except for Vasto and San Salvo. The reason lies in the recent development of coastal tourism. Indeed, in recent years, the Abruzzo Region has invested in the tourism promotion of its territory through the Strategic Tourism Plan (2017-2019), particularly in the Teatina coast (where Vasto and San Salvo belong). Besides, the mentioned municipalities are those with greater availability of sandy beaches, more suitable for seaside tourism. Looking inside the change rate in the urbanization of the coastal and maritime areas, instead, the situation is growing everywhere. Speaking about the land take speed in the Abruzzo Region it is possible to say that before 2000, the average land take speed was equal to 17,500 m 2 /day. As for the general behavior of the Italian Peninsula, the Abruzzo Region has its maximum increase in urbanization between 1960 and 1980, started to decrease from 1980 and 1990 and has stabilized between 1990 and 2010. At present, the rate is almost unchanged and future inferential projection until 2050 does not show significant variations [60]. Figure 5 shows the variation of urbanization in the whole coastal area of the Region while Figure 6 describes more in detail the Urban change for each management unit. As previously underlined, in order to investigate the long-term effect of the urbanization growth, the shoreline variation between the epochs 2001 and 2018 has been investigated. The variations have been evaluated from north to south using Equation (1).
The obtained results (see Figure 7) reveal the impossibility to define a general trend describing the evolution of the shoreline in all the considered management units. There is a clear alternation of retreats and advances of the shoreline against the always positive trend of urban growth.
In general, sediments tend to remain trapped in the updrift area near each port delimiting each MU. This aspect is expected considering the main direction of the sediment transport in the Abruzzo Region (waves mainly come from the north).
In the northern part, the highest retreat of the shoreline is observed, in particular in the first management unit where the extension and width of the retreats are evident and localized far from the Giulianova harbor (IAs 89-90).
The second MU, between Giulianova and Pescara harbors is characterized by a similar pattern of the first MU, a shoreline advance near the ports. The third management unit shows, instead, a general trend to shoreline growth, except for some local hot-spots.  In order to assess the effect of the coastal defense systems on the shoreline evolution, the Protection Index has been evaluated for the considered case study. Typically for the whole Adriatic coast, the stretch of the Adriatic coast on the Abruzzo Region has a large fraction of the coast protected by groins and breakwaters (both emerged and submerged). Figure 8 shows the protection index. It has to be noticed that the protection structures have been built from the early 60s up until today, therefore it is not possible in this case to indicate a fixed time period as for the variation in urbanization and shoreline.
However, the inspection of Figure 7 reveals that, aside from the first management unit, where only the northern part is not protected by the structures, in other cases the total amount of protection covers almost all their extension. It is also quite clear that in most cases, due to an increase of urbanization (therefore to an increase of exposure of the area) corresponds to a growth in coastal defense. This evidence confirms the direct and close relationship between urban growth, expansion and the following protection of coastal zones. The causal relationship is mutual. Urbanization influenced the natural coastal evolution, i.e., by affecting the coastal dunes system. On the other hand, the presence of building, infrastructures, and economic activities has led to an increase of beach width demand. Shoreline retreat threatened private houses, streets, or economic activities increasing the risk. Therefore, adopted solutions were to plan beach nourishment and rigid defense even if other possible solutions could have been ecological, and biological remediation aimed to increase beach width by promoting natural morphodynamics. However, comparing the results shown in Figures 7 and 8 it is possible to observe that although the PI is in a lot of cases equal or close to one, there is an alternation of retreats and advances of the shoreline. Furthermore, it is difficult to give an unambiguous interpretation of the role of the rigid structure against coastal dynamics. As expected, an increase of protection in the same area, reveals a general tendency in advancement in the shoreline position. Nevertheless, despite a large amount of defended areas, the shoreline changes in response to urban growth are always characterized by an alternation of positive and negative values. This is likely related to two main factors. On the one hand, the coastal defense works were firstly carried out at river mouths: the trigger effect is related to the decrease of sediment transport from rivers. This occurrence is mainly due to the sediment that remains trapped in a series of artificial reservoirs (e.g., Campotosto and Bomba Lakes), and a series of small dams for hydroelectric energy generation insisting on the Region.
On the other hand, boundary effects of coastal works have enhanced the erosion process on the downdrift areas.
Trying to move deeper inside the relationship between protection, shoreline variations, and urbanization, it is possible to qualitatively analyze some evidence in the MUs.
The first MU is characterized by protection only in the northern sector (25 IAs, see Figure 8) and almost all the breakwaters have been built in the period 1989-1997. Looking at the shoreline variation it is clear that, in the first 25 IAs, the structures achieved their goal. Indeed, despite some IAs where a small retreat happened (in the downdrift area of a small protected area for small fishing boats), the shoreline shows positive variations. In that areas, a series of groins have been placed in order to counteract the retreat. However, the downdrift stretch (up to the IA number 40) shows an evident retreat of the shoreline and no protection has been placed in this case. Nevertheless, it has to be stressed that it is difficult to directly correlate the presence of the structures and the shoreline retreat in the downdrift area. Indeed, the shoreline retreat can be related to the lack of sediments rather to the boundary effects of updrift coastal structures. In this sector, the urbanization variation is always positive, but there are some areas with zero variation, probably induced by the lack of available beaches for touristic enjoyment. After this sector (IAs 25-40, first MU), the effects of the harbor in retraining sediment are evident. Indeed, no structures are needed and the shoreline grows everywhere. The urbanization variations are always positive, but some areas near the Salinello's river mouth characterized by the presence of a great number of green areas.
In the second MU there is the presence of two protected areas "Riserva Naturale Del Borsacchio" and "Torre del Cerrano". In both cases, as it is possible to see in the Figure 8 there are no protection structures. In the other stretch of the MU, instead, the PI is always positive. The great amount of the structures were deployed in the decade 1988-1998 and only some new structures have been realized in the last twenty years. As previously underlined it is difficult to assert what is the real benefit of rigid protection. In general, the alternation of advances and retreats combined with the small amount of the advance in shoreline position, does not allow an easy interpretation of the actual situation. It is interesting to note that the area near the Pescara harbor (physical limit of this MU, realized in 1989) is characterized by high values of ∆ URB . The effect of the harbor in trapping the sediment is evident in the positive ∆Y values. The area has been protected since 1989 by a series of breakwaters.
Moving to the south, the third MU shows a non-uniform growth in urbanization. The section with wider development is the nearest to the metropolitan city of Pescara while the remaining part of the MU has been left quite intact. The shoreline variations show a general non-negative behavior, with a general advance. This is probably because the third MU is almost completely protected by rigid structures (realized in the decade 1988-1998). The low urbanization values combined with the protection strategy probably lead to this situation. It is worth noting the newest protection structures, realized in 2007, are located downdrift of the Pescara harbor to counteract the negative values of ∆Y, preserving the economic activities and the enjoyment of the beach.
Finally, the MU 4 is in general characterized by low values of ∆ URB . The ∆Y is characterized by an alternation of advances and retreats. As is possible to observe in Figure 8 the area has a great number of IAs with high PI values. The non-protected areas are related to the regional protected area of "Riserva naturale di Punta Aderci" and to some areas with no sandy beaches. In this case, a great number of structures are referred to 1988-1998. The new constructions (i.e., groins, breakwaters) from the period 2007-2019 are in general localized works, complete an existing protection system and account for 25% of the total amount of rigid structures of the entire coastal stretch. Figure 9 shows a general overview of the temporal evolution of the PI in the entire coastal stretch while Figure 10 shows some details of selected areas for each MU in order to highlight the analysis that has been carried out for taking a census of the rigid protections.  However, a possible way to quantify the effect of the rigid protection in the study area could be the analysis of the sediment budget (i.e., m 3 /year) for each MU. The inspection of the available regional documentation, for the coastal stretch of the Abruzzo Region, the annual value for the closure depth (d c ) is about 6.85 m, while the extension (∆ x ) of the IA is 500 m. Therefore, the yearly sediment volume can be evaluated as ∆Y × ∆ x × d c . Table 2 shows the obtained results. As it is possible to see (e.g., Table 2), the general sediment budget is positive for each MU. It has to be stressed that this does not mean the total coastal stretch of each MU is locally advancing, but that the total yearly budget is positive. However, there are areas of negative volume rates compensated by others with higher values, giving final positive results. In the first MU, about 54.5% of IAs show shoreline retreat, while in the other MUs, the values are about 33.4%, 13.5%, and 33.4% respectively.
The presence of human activities, and the exploitation of coastal areas for economic and touristic purposes, suggested a further investigation on the economy of a coastal area. Considering (as previously underlined) only the activities in the coastal area, there are just more than 16,000 employees (situation at 2017), equivalent to 4.4% of the total residents of the coastal area.
As it is possible to see from Figure 11, analyzing the ranking of the Abruzzo Region in the national framework, appears that it is in the middle place considering the rank for people working in the selected business sectors. The evaluation of the ILU, reveals a generally positive trend. As shown in Figure 12, the first MU is characterized by a general increase in ILU index, except for the municipality of Tortoreto. It has no rigid structures (i.e., PI = 0) but it is one of the municipalities with the highest growth of urbanization. In the second MU, there are two municipalities showing negative ILU values, Pineto and Silvi Marina. In this case, the ∆Y shows an alternation of retreats and advances, with positive values of urbanization growth but, not always evident as for the previous case. The third MU, instead, shows values as low to be considered equal to zero. The southern part (the fourth MU), on the contrary, shows higher values of the considered index indicating a generally positive trend in the growth of the local economy in the coastal area. It has to be stressed that this MU is the only one that shows significant changes in ILU, with (in the case of Fossacesia and San Vito) percentage variation reaching an increase of about 15-20%. Analyzing the local economy, it appears that there is no general trend. However, an increase in urbanization seems to suggest a generalized benefit for local activities. On the other hand, a coastal area prone to shoreline retreat does not always show a decrease in local economy. Furthermore, the attempt to increase the beach width with coastal protection measures does not lead necessarily to an increase in the local economy.
In the Abruzzo Region, the greater performances are linked to the "Costa dei Trabocchi", the most attractive area from a touristic point of view.

Concluding Remarks and Ongoing Research
This paper aims to analyze and quantify the relationship between urban development, coastal dynamics and the possible effect that the variation of these variables can induce in the local economies. The analysis was performed for the coastal area of the Abruzzo Region (Adriatic Sea, central Italy) by using different types of datasets dividing the coast into four management units considering the harbors of the region as physical boundaries for each unit. At a first glance, results reveal that a direct causal relationship between urban development and erosion phenomena is difficult to find and there is no clear evidence that an increase in urbanization induces a shoreline retreat. This occurrence is due to the fact that, during the years, coastal areas have been protected by using breakwaters and groins in order to fix the local position of the shoreline modifying the hydrodynamics of the area and, therefore, the morphodynamics. Then a protection index, able to weigh the total amount of rigid protection structures against the shoreline length was proposed. Results have shown a general over-protection of the coast using different types of engineering solutions. In order to counteract the progressive increase of beach width demand almost every coastal municipality resorts to hard structures. However, despite this over-protection, the effect on the coastal dynamics and, in particular on the increase of beach width, are less evident than expected.
In general, the majority of the structures induced a localized benefit, not always producing appreciable changes in the shoreline position. Results reveal that the total protection of the MU lead to a general benefit for the beach's width, although generating, in some cases, a series of confined cells (i.e., combining groins and breakwaters) with resulting problems in water quality. However, a clear and generalized conclusion on the real benefit of rigid protection is difficult and it can be assessed only after a comparison with the "do nothing option". Due to the practical difficulties in modeling a "do nothing option" scenario (out of the scope of this paper), the conclusion on the performances of the structure can be only argued. In the case of a predominant wave climate and, therefore a uni-modal direction of sediment transport, groins can induce localized benefit in the updrift area, paying attention to the effect on the downdrift sector. Breakwaters seem to be a good attenuation of wave action, useful in particular in the case of protecting beach nourishment avoiding or mitigating the loss of sediments. However, the intensive use of alongshore protection involves benefit, often in the face of sustainable maintenance of natural habitats. Moreover, an analysis of the local economies shows that a generalized protection of the coastal area does not produce a definite benefit to the local economy (i.e., the third MU has high values of PI while ILU can be considered equal to zero). On the contrary, areas with total or partial absence of human presence (i.e., the "Costa dei Trabocchi" area) show significant changes in ILU and, therefore, in the local economies. These results suggest the potential of sustainable tourism in the future development of coastal areas.
The ongoing research is focused on the extension of this work to a national scale in order to have a general picture of the various effect produced by different regional policies in the coastal stretches. The aim will be to understand how different urbanization plans and different engineering solutions act on the coastal system also in view of future climate change scenarios.