Territorial Vulnerability Assessment Supporting Risk Managing Coastal Areas Due to Tsunami Impact

Portugal’s coastline extends 1187 km. It is characterized by social, economic and physical conditions that differentiate it from the rest of the territory, including population density, location of infrastructure and support of tourism activities. Therefore, it has a significant exposure if a tsunami occurs. Six coastal study sites with varying characteristics were selected for evaluation in this paper, including two core beach-use areas, two residential areas and two industrial areas. These sites are located in the municipalities of Figueira da Foz, Setúbal and Vila do Bispo. The analysis began with the calculation of the potential tsunami inundation area for each site using the 1755 Lisbon tsunami. Next, a methodology distinguished by its multidimensional character was applied to assess local vulnerability to tsunamis. This methodology assesses vulnerabilities associated with morphological, structural, social and tax factors. These four vulnerability components were combined to obtain a Composite Vulnerability Index (CVI), which enabled us to identify the most vulnerable areas and to determine the distinguishing characteristics of each area.


Contextualization
The 1 November 1755 earthquake (magnitude Mw 8.7) and tsunami was the most catastrophic natural disaster to occur in Portugal's history; this disaster caused approximately 12,000 deaths and resulted in severe damage to many coastal areas, particularly the city of Lisbon [1].Historical records identify a total of 17 tsunamis that have taken place in Portugal from 60 BC to the present day [2,3].
In recent decades, Portugal's coastal areas have been subject to social, economic, land-use and environmental transformations, similar to coastal areas worldwide [4][5][6].The concentration of population in urban and coastal areas in recent decades has put great socio-economic and environmental pressures on these zones.Worldwide, coastal areas play a key role in supporting human survival; about half of the world's population depends directly on these areas for livelihoods related to fisheries, tourism, trade, transport and energy [6], and 40% of the world population lives less than 100 km from the coast [7].
If we narrow our analysis and only consider areas designated Low Elevation Coastal Zone (LECZ), which is defined as a contiguous area along the coast less than 10 meters above sea level [8,9], we see that this area represents 2% of the world's land area but supports 10% of the global population and 13% of the world's urban population [9].
In Portugal, the coastal zone is 1187 km long and is characterized by social, economic, physical and environmental conditions that distinguish it from the rest of the territory.In the last five decades, urbanization, especially along the coastline between Viana do Castelo, Setúbal and the Algarve coast, has put strong human pressure on this zone, which is represented by the fact that 75% of the country's population and 85% of its GDP are now concentrated in this area [10].These data are supported by the European Environmental Agency (EEA), which notes that between 1990 and 2005 Portugal had one of the highest percentages of construction per 1 km of coastline of any European country [11].These facts have decisively increased the vulnerability and exposure of the coastal zone to natural hazards in general, and tsunamis in particular, as revealed by the damage caused by the 6 January 2014 Hercules storm, which caused more than a million in damage and about eight injuries [12].
The heterogeneity that characterizes the Portuguese coast is reflected in the existence of numerous areas with differentiating features.Industrial areas, a succession of predominantly residential areas, as well as beach-use areas for vacationers and environmentally protected areas are all found along the coast.This diversity requires the adoption of a multidimensional methodology to analyze and evaluate the tsunami vulnerability of different regions.This study analyzed six sites based on the following distinctive features: two beach-use areas, two residential neighborhoods and two industrial areas.
The beach use areas are characterized by strong coastal tourism.They are distinguished from other areas by their high transient populations and low number of long-term residents, a high number of secondary vacation homes, and a strong seasonal flux in the human population that puts amount of strain on the environment and landscape.The two residential neighborhoods were chosen based on their strong anthropic pressures corresponding to population clusters of regional (Cova Gala) and national (Setúbal) importance.These areas distinguished by the variety of services and functions that are provided within their territories.The two industrial areas selected for this study are also characterized by the presence of residential, retail and other services in addition to industry.
The present study will focus on different areas within three Portuguese coastal municipalities, based on historical accounts of the 1755 Lisbon tsunami and territorial characteristics (Figure 1).Even though it was triggered by an earthquake, the methodological approach only explores the tsunami parameters and their relation with coastal attributes for the potential impact assessment.In this study, the impacts from earthquake shaking are not considered, but this would likely increase the damage and number of injuries in the areas studied.This could also potentially cause changes in the land morphology as well as weaken building structures.The analysis began with a worst-case scenario based on a deterministic approach, focused on a scenario similar to the 1 November 1755 tsunami event.Table 1 presents the tsunami parameters for the three study municipalities.Furthermore, according to the historical accounts, Vila do Bispo municipality was the first municipality hit by the 1755 tsunami, less than 20 min after the earthquake [1].In addition, at Vila do Bispo the tsunami run-up was extreme, being the highest values with at the cliffs.As we move away from the source (eastern and northern Portuguese coastal zones), the tsunami travel times increase and tsunami run-up decreases, as confirmed by the tsunami numerical model [2,13].->2000 (1)   Vila do Bispo 10-144 16-17 0 Note: (1) damage/fatalities due to earthquake, fire and tsunami.
Based on the tsunami numerical model carried out by [2,13] at the regional scale, the study uses the inundated areas calculated by [1] at Figueira da Foz, and the outputs model upgrades at Setubal and Vila do Bispo.
A multidimensional methodology was applied to all six study sites to assess their vulnerability to tsunamis, expressing the potential for casualty, destruction, damage, disruption or other form of loss.The selection of this methodology was justified by the fact that the six areas have distinctive morphological structural, social and tax features.Previous studies evaluating the vulnerability to tsunamis show that the distribution of vulnerability over a territory does not occur in a uniform manner and is instead dependent on a number of factors [14][15][16].Thus, this paper applies a multidimensional methodology that integrates morphological, structural, social and fiscal vulnerability components into the analysis.A vulnerability index was calculated for each component.Finally, a Composite Vulnerability Index (CVI) was calculated to create a single index that accounted for the four vulnerabilities previously considered.
The objectives of this work include: 1 Analysis and assessment of the characteristics of beach-use areas, residential areas and industrial areas. 2 Assessment of the characteristics of the buildings potentially affected by tsunami.With this approach, we intend to assess the characteristics of the territory subject to potential tsunami inundation impact.We seek to evaluate both structural elements potentially exposed as the social and economic dimension of the affected study areas.This knowledge will allow the resource scale for prevention and mitigation, as well increase the warning, alert, rescue and emergency effectiveness.

Study Areas
The study area was defined by overlaying the worst-case scenario inundation results (where the maximum water level is related to the mean sea level) from the numerical tsunami model obtained by [2], with the statistical subsections from Censos 2011 [17].Three distinctive land-uses were selected for analysis: beach-use areas (Cabedelo and Martinhal), residential neighborhood (Cova Gala and city of Setúbal) and industrial (Outão and industrial zone of Setúbal).Each land-use was represented by two sites from Figueira da Foz, Setúbal and Vila do Bispo municipalities.

Beach-Use Areas
Cabedelo (Figueira da Foz Municipality) Cabedelo is a seaside area located in central Portugal in the parish of São Pedro (Figure 1b).The parish is located in the northern Cova Gala population cluster (Figure 2) of Figueira da Foz municipality.The municipality consists of 18 parishes distributed over an area of 379 km 2 .The municipality is home to 62,125 inhabitants and has a population density of 163.3 hab/km 2 , a total of 24,893 buildings and 43,300 lodgings [17].The parish of S. Pedro has an area of approximately 7 km 2 with a resident population of 2910 people and a population density of 415.1 hab/km 2 [17].However, the number of residents and non-residents increases considerably during the summer season [18].
Cabedelo is distinguished by a plurality of deployed services.In particular, it has a camping site, coffee shops, restaurants and a surf school, and it is home to the fishing port of Figueira da Foz, shipyards and a group of companies related to the maritime sector.

Martinhal (Vila do Bispo Municipality)
The Martinhal beach (Figure 3) is located in Sagres parish (Figure 1d), which is part of Vila do Bispo municipality in the most southwestern zone of Portugal.Martinhal is characterized by seaside tourism dominated by single-family residences that cater to tourists or serve as second homes.
Vila do Bispo municipality contains five parishes and is 179 km 2 .According to the 2011 Census, the municipality has 5381 inhabitants, with a population density of 29.2 hab/km 2 [17].The municipality has a total of 5184 buildings and 5979 lodgings [17].Sagres parish has 1909 inhabitants, is 34.28 km 2 , and has a population density of 55.7 hab/km 2 [16].Martinhal stands out as a tourist retreat and is almost exclusively devoted to tourism.The Martinhal Beach Resort & Hotel is almost exclusively responsible for all of the buildings within the analyzed area.This area is also characterized by a marked seasonal population.

Residential Neighborhoods
Cova Gala (Figueira da Foz Municipality) Cova Gala is a population cluster located on the left bank of the Mondego River, in the São Pedro parish (the same parish where Cabedelo is located) of Figueira da Foz municipality (Figure 2).In recent decades, the parish of São Pedro, especially the area of Cova Gala, has undergone considerable territorial and socio-economic changes.Once characterized by its intrinsic link to fishing, the importance of this sector has been reduced over the last four decades.The implementation of various industries in the parish and neighboring areas, as well as the opening of the Hospital District of Figueira da Foz, the implementation of a set of services and the consequent development of the tertiary sector and the development of tourism, have all led to profound changes in the area and an increase in the resident population.These facts led to the elevation of São Pedro-Cova Gala from an agglomeration to a town in 2009.

City of Setúbal (Setúbal Municipality)
The other residential neighborhood analyzed in this study is the riverside urban area of the Sado estuary in the city of Setúbal (Figure 4).The municipality of Setúbal is located in the metropolitan area of Lisbon, is 172 km 2 [17] and is divided into eight parishes.The municipality has a total of 121,185 inhabitants and a population density of 698.9 hab/km 2 .In national terms the city of Setúbal is a medium-sized city and is one of the primary national economic centers.It contains commercial and industrial sectors and is an example of the diversity of infrastructure and distribution of facilities along the coast [12].
In terms of use and occupied land, the study area contains residential features, as well as areas allocated to services, commerce and activities related to industry and fishing.According to [19], the area can be divided into four sectors based on deployed equipment and infrastructure: (1) Near the riverfront, which is dominated by support infrastructure for fishing and commercial activities; (2) From the riverfront to Av. Luisa Todi, which is dominated by equipment and infrastructure related to the port of Setúbal, civil protection, services, catering and some residential buildings; (3) Av.Luisa Todi, which is dominated by green spaces; and (4) The historic center of Setúbal, which is characterized by residential buildings and spaces for food services and commerce.

Industrial Areas Outão (Setúbal Municipality)
Outão is located in Arrábida Natural Park in the parish of Nossa Senhora da Anunciada in the southwestern section of Setúbal municipality.The focus of this area is its cement industry, which clearly delineates the region from an environmental, landscape and socio-economic point of view (Figure 5).The area is characterized by industry and services, including the Setúbal orthopaedic hospital, but and also the camping of Outão.It has a distinct non-resident population, and there is significant daily flow into and out of the area.

Industrial Zone of Setúbal (Setúbal Municipality)
The other industrial area analyzed (Figure 6), which includes the riverfronts of São Sebastião and Sado parishes and extends from the end of the port to the beginning of the peninsula of Mitrena.It has an area of 1.25 km 2 .Diverse deployed infrastructure, ranging from a power station and a wide variety of chemical industries to support infrastructure and various port activities characterize this area.This area as well as neighboring areas has a high density of road infrastructure.The resident population is unique among the studies sites in that it is composed of a set of informal settlements and structurally degraded buildings in the surrounding areas of the old industrial complex.

Data and Methods
A new multidimensional approach has been developed to tsunami risk management by combining four components of vulnerability analysis: morphological, structural, social and tax components of vulnerability.These factors were selected based on previously studies: morphologic component [20,21]; structural component [22]; social component [23][24][25][26] and allocation component based on property taxes [27].The analysis of the different vulnerability occurred at the municipal level (social and tax component) and local level (morphological and structural component).At the municipal level, the unit of analysis is the parish, and at the local level the unit of analysis is based on the inundated area as identified from the numerical tsunami model.
The analysis will culminate with the calculation of a Composite Vulnerability Index (CVI), which combines the four vulnerability components discussed above.All vulnerability indices (the four individual components and the final composite index) vary from very low to very high in accordance with the standard deviation and the following categories: "very low", <1 SD; "low", ≥−1 to −0.5 SD; "moderate", ≥−0.5 to +0.5 SD; "high", ≥0.5 to 1 SD; "very high", ≥1 SD [23].characterized based on the explanatory variables included therein.Each of the explanatory variables can positively (+) or negatively (−) influence the FAC, as appropriate.The remaining six components were not considered because they contributed residual values to the total variance.

Taxable Property Assessment (TPv)
A set of four allocation coefficients (Figure 8) was considered for the analysis of the tax factor and served as a basis for the calculation of the Municipal Property Tax.The coefficients were derived from various activities implemented throughout the territory, such as commercial, industrial, housing and services.These data were collected from the Tax and Customs Authority.The incorporation of these new elements reinforced the analysis by providing information about variations within the analyzed territory.The different allocation coefficients for activities enabled the classification and differentiation of each territory based on integrated activities.
The product of the coefficients mentioned above was used to calculate the Taxable Property Vulnerability (TPv) index.After the calculation was performed, the values were adjusted to a range of 0 to 1 using an algorithm.This resulted in five vulnerability classes based on the standard deviation classes mentioned above.

Morphological Assessment (Mv)
Five parameters were considered in the calculation of morphological vulnerability (Figure 9).A matrix based on fieldwork complemented with satellite imagery analysis and the use of GIS software was employed to collect and analyze the attributes of each parameter.The matrix gives values ranging between +1 and −1 to each attribute.Positive values represent an increase in vulnerability and negative values represent a decrease in vulnerability.However, not all parameters have the same influence on morphological vulnerability.Therefore, we used the M-Macbeth ® (Hainaut, Belgium) multi-criteria analysis software [33,34] to assign different weights to each parameter (Figure 9).Morphological vulnerability was calculated using the formula presented in Figure 9, where "pn" is the value assigned to each parameter and "wn" is the weight assigned to each parameter.After the calculation was performed, and the algorithm in Figure 1 was used to adjust the final values to a range between 0 and 1, where "x" represents the direct value of the vulnerability, "y" is the adjusted vulnerability value in the range 0 to 1, "a" is the maximum value of the vulnerability of the sample and "b" is the minimum value of the sample.

Structural Building Assessment (Bv)
A structural component analysis was applied to buildings potentially affected by the numerically modeled tsunami inundation [2].We considered a total of 13 parameters (Figure 10) that characterize the structural, architectural, occupational and functional aspects of the buildings.The identification of these buildings and the collection of their various intrinsic characteristics and those of the surrounding areas occurred through fieldwork.A matrix was used to aggregate all parameters representing the structural aspects of vulnerability.Each parameter was given a value between +1 and -1.
For this assessment, we introduced a new set of parameters (bold in the Figure 10) related to the occupation, functionality and structure of each building [31,35] that have not been addressed in previous works, including those by [21,22].The introduction of these new parameters strengthened our analysis and enabled the introduction of new data that had not been considered in previous methodologies.These new parameters distinguish buildings according to their function, the seasonality of use, capacity, and rate and type of occupation.In terms of function, buildings were divided into six categories: exclusively residential, commercial, services/industry, outbuilding/warehouse, mixed type (e.g., residential and commercial building), and critical buildings (e.g., police, fire station, and school).
Each of the 13 parameters contributed differently to structural vulnerability.M-Macbeth ® multi-criteria analysis software [33,34] was used to assign different weights to each parameter.
After weighting each parameter and assigning values to each of its attributes, structural vulnerability was calculated (Bv) using the formula presented in Figure 10, where "ppn" is the value assigned to each parameter and "wn" the weight assigned to each parameter.After the calculation, the adjustment algorithm was applied to convert the data to values between 0 and 1.

Composite Vulnerability Index (CVI)
After calculating the morphological (Mv), structural (Bv), social (Sv) and taxable property vulnerabilities (TPv), a Composite Vulnerability Index (CVI) was calculated.This index is the sum of the four component indices and is presented in Figure 11.
The CVI varies between 0 (minimal vulnerability) and 4 (maximum vulnerability), with vulnerability ranging between very low to very high according to the five standard deviation levels mentioned above.

Social Vulnerability
Social vulnerability was calculated at the parish and municipality levels for the different areas of study: Figueira da Foz (Cabedelo and Cova Gala), Setúbal (Outão, city of Setúbal and industrial zone) and Vila do Bispo (Martinhal).
Figure 12 is a spatial representation of the social vulnerability of each municipality.In Figueira da Foz (Figure 12a), vulnerability varies between low (0.310) and very high (0.730) (Table 3).Cova Gala and the Cabedelo, which are located on the left bank of the river Mondego, in the parish of São Pedro exhibit very high social vulnerability.
In the district of Setúbal, social vulnerability varies between very low (0.000) and high (0.670) (Table 3).Of the areas analyzed within the municipality (Figure 12b), Outão, which is located in the moderate social vulnerability.Overall, the lowest vulnerability values are observed in Setúbal municipality, where values for the two lowest vulnerability classes (very low and low) ranged from 0.000 to 0.400.The highest vulnerability values were observed in the municipality of Vila do Bispo, where values for the two highest vulnerability classes (high and very high) ranged from 0.580 to 0.770 (Table 3).

Taxable Property Vulnerability
Taxable property vulnerability (Figure 13 and Table 4) varies from very low to very high in Figueira da Foz (Figure 13a), with marked differences across the municipality.Cova Gala and Cabedelo were characterized by very low vulnerability in the south of analyzed areas and moderate vulnerabilities in the remaining areas.These results can be explained by the fact that this is a peri-urban area with strong natural components, with the exception of Cova Gala, which is a consolidated urban area, and the fishing port, which concentrates trade activities for the territory.
In Setubal (Figure 13b), taxable property vulnerability also varies between very low and very, high with a clear distinction between urban and natural/rural areas.The highest vulnerability values are concentrated in the center of the municipality, at the heart of Setúbal's urban area.In terms of the localities analyzed, Outão exhibited moderate vulnerability.This can be explained by the heterogeneity of the area.Despite the existence of infrastructure linked to industry and services, the area also has a strong natural component and is located in the middle of the Arrábida Natural Park.The city of Setúbal, which is a consolidated urban area, is completely characterized by very high vulnerability.The riverfront, which is primarily geared toward industrial uses, is characterized by varying taxable property vulnerability, with the port of Setubal classified as highly vulnerable and a general decrease in vulnerability to the east, where the vulnerability is very low.Finally, the vulnerability of taxable property in Vila do Bispo municipality (Figure 13c) ranges from very low to very high, with a clear distinction between urban and rural areas.Martinhal is classified as very highly vulnerable.Although Martinhal is surrounded by rural/natural features, its taxable property vulnerability value is high because it is a center for seaside tourism, and therefore its location coefficients are given higher values.
Table 4 allows us to assess the numerical variation of different vulnerability classes for the municipalities in question.Overall, it can be concluded that lower levels of vulnerability are observed in the municipality of Vila do Bispo, where the lower class values ranged between 0.000 and 0.630.Higher taxable property vulnerability values for Vila do Bispo varied between 0.717 and 1, highlighting the heterogeneity of vulnerability across the municipality with respect to tax characteristics

Morphological Vulnerability
Figure 14 shows the morphological vulnerability in the area under study, with vulnerabilities ranging between very low and very high.Cabedelo (Figure 14a) and Cova Gala (Figure 14b), located in Figueira da Foz, represent 6.1% of the total flooded area of the municipality (Table 5).The highest morphological vulnerability values are concentrated in the coastal zone, especially on the beach of Cabedelo, in the fishing port and on the beach of Cova Gala.It is apparent that morphological vulnerability is strongly correlated with distance to the coastline, as decreases in vulnerability are associated with increasing distances from the coastline.However, within this factor, slope and morphology of the study area also play a decisive role in the final value of vulnerability.With regard to areas located in Setúbal, the urban area of Setubal accounts for 2.7% of the total flooded area within the municipality (Table 5), and its morphological vulnerability varies between very low and very high (Figure 14c).Lower values are located in the boundary areas in the northern part of the study area, especially at the edge just beyond the north of the marina.The highest values of morphological vulnerability are located along the riverfront of the Sado estuary to approximately 200 m inland.Beyond this band, the region is dominated by moderate values of vulnerability, which affect the vast majority of the study area.The final values of morphological vulnerability are clearly influenced by the following parameters: use and type of soil, slope and distance to the coastline.
In Outão (Figure 14d), vulnerability varies between very low and very high; this area contains 0.5% of the total area affected by flooding within the municipality (Table 5).Vulnerability of this area varies greatly, and the distance to the shoreline parameter has less weight compared to the areas discussed above.The highest vulnerability values are located in the Outão camping location and the areas surrounding the Setúbal orthopedic hospital.The final values of morphological vulnerability are clearly influenced by land-use, land cover and slope.
Regarding the industrial area (Figure 14e), morphological vulnerability varies between very low and very high, with the greatest vulnerability values located along the edge of the Sado estuary and the lowest located inland.In this area, the final values of morphological vulnerability are clearly influenced by the distance to the shoreline and use and land cover.This industrial zone accounts for a total of 4.8% of the total flood-affected area within the municipality.
Morphological vulnerability in Martinhal (Figure 14f) ranges from low to very high and is greatly influenced by distance to the coastline.Figure 14f clearly shows that vulnerability in Martinhal decreases with increasing distance from the shoreline.The highest values are located in the area of Martinhal beach, in the southern zone of the Martinhal beach resort hotel and on the sandy beach of Rebolinhos.Martinhal represent 12.3% of the total flooded area of the municipality (Table 5).
Table 5 shows the total percent of flooded area per morphological vulnerability class.Higher values (high and very high class) are observed in Cabedelo (80%), Cova Gala (80%) and Martinhal (57%).
In contrast, lower levels of vulnerability (very low and low classes) are located in Outão, where the two lowest classes represent 22% of the flooded area.Also noteworthy is the fact that in the city and the industrial zone of Setubal, the vast majority of the flooded area is characterized by only moderate morphological vulnerability (78% and 62%, respectively).If we analyze the data in terms of affected area, we see that the high and very high vulnerability classes in the Cabedelo zone are some of the largest affected areas (56 ha).With regard to the very low and low vulnerability classes, represented by the industrial area, 5 ha are affected.In terms of moderate vulnerability, 55 ha of urban Setubal are affected.

Structural Vulnerability
Figure 15 shows the structural vulnerability, with vulnerability varying between very low and very high.By overlaying the 1755 Lisbon tsunami inundation model with current buildings footprints, we identified a total of 22 buildings in Cabedelo that would potentially be affected by a tsunami (Figure 15a).The structural vulnerability of these buildings varied from low to very high.Table 6 shows that most of the buildings exhibited moderate to high vulnerability.Most buildings of concern were in the outbuilding/warehouse, services/industrial and commercial categories.With regard to residential buildings, only one was potentially affected by the worst-case scenario tsunami.Despite a strong industrial and commercial component, Cabedelo is also an important seaside tourism area characterized by a high influx of people on weekends and during the summer.The camping area has a maximum capacity of 1000 people and is open from January to November.This zone is also much appreciated by surfers because of the natural beach and sea conditions and is a popular choice for their activities.In Cova Gala, a total of 15 buildings (Table 6) were potentially affected by a worst-case scenario tsunami (Figure 15b), with vulnerabilities ranging between very low and high.An analysis by vulnerability class shows that most buildings at this site have low vulnerability.The majority of buildings (13) are exclusively residential, one has a commercial function and the other is an outbuilding.Of the 13 buildings that are exclusively residential, eight exhibit low vulnerability and only two are highly vulnerable.
In the urban area of Setubal, a total of 1259 buildings (Table 6) were identified as potentially affected by the worst-case scenario tsunami (Figure 15c), with vulnerabilities ranging between very low and very high.Of the total, 18% of buildings have a very low structural vulnerability, 16% low, 41% moderate, 13% high and 12% very high vulnerability.Fifty-five percent of the buildings are mixed use, 29% of the buildings are exclusively residential and 22% have a secondary function in addition to a residential function.5% of potentially affected buildings relate to critical infrastructure and are of strategic importance for emergency response and normal societal functioning.Some of these buildings are shown in Figure 15c and include: (1) a fire station with moderate vulnerability; (2) a kindergarten with high vulnerability; (3) a courthouse with low vulnerability; (4) a port captaincy office with high vulnerability; (5) a customs house with very high vulnerability; (6) a police station with high vulnerability; and (7) a port administration building with moderate vulnerability (the methodological approach only considered the port buildings for the structural vulnerability analysis, because the on-water port structures are presented in the morphological vulnerability, namely through the land use occupation).
In Outão (Figure 15d), a total of 32 buildings were identified as potentially affected by the worst-case scenario tsunami (Table 6), with vulnerabilities varying between very low and very high.Most of the buildings exhibited a moderate vulnerability.Typologically, the analysis of potentially affected buildings shows that the vast majority of buildings (22 buildings) are associated with existing cement factory and are industrial in nature.Eight buildings are located in the camping area, including a reception building, a coffee shop as well as a number of support buildings.There are only two exclusively residential buildings in the potentially flooded area.
Regarding the industrial area (Figure 15e), a total of 325 potentially affected buildings were identified, with vulnerabilities varying between very low and very high (Table 6).The analysis of buildings by vulnerability class shows that 11.7% have very low structural vulnerability, 8.6% low, 26.4% moderate, 18% high and 35% very high structural vulnerability.In typological terms, the buildings affected include outbuilding/warehouses and support buildings in the industrial and commercial zones, which represent 31% of all potentially affected buildings (102 buildings); exclusively residential buildings (86 buildings, 26% of total) and buildings related to services/industry (65 buildings, 20% of total).Additional affected buildings include critical infrastructure such as power stations, industrial warehouses and support buildings for the Setúbal port's activity.
Finally, at Martinhal (Figure 15f) a total of seven potentially affected buildings were identified (Table 6), with vulnerability levels varying between very low and very high.Three buildings are exclusively residential, three have a commercial function, one is an outbuilding and one is a wastewater pumping station.Although the number of buildings potentially affected may seem low, Martinhal is a natural area with a strong beach tourism presence and is characterized by fluctuating population flows that increase on weekends and during the bathing season, when the entire beach has the potential to be completely inundated.

Discussions and Conclusions
This study analyzed areas with different characteristics (beach-use, residential and industrial areas) to understand the different dynamics associated with these land-use characteristics.After calculating the morphological, structural, social and taxable property vulnerabilities, the Composite Vulnerability Index (CVI) that results from the combination of these components was calculated, and is presented in the Figure 16.
The results highlight different behaviors and dynamics between study regions.Table 7 shows the percent of affected area by CVI class for all three areas categories analyzed.Beach-use areas clearly have a higher CVI when compared to the industrial and residential areas.Seventy-four percent of affected beach-use areas is classified with a high or very high CVI, with a clear emphasis on morphological component of vulnerability associated with distance to the coastline and slope.This is evidenced by the existence of the highest CVI values near the coast and in areas of lower slopes, as well as decreases in vulnerability with increasing distance from the coastline (Figure 16a,f).It is also analyzed; areas with buildings largely had CVI values between moderate and very high, while the rest of area had CVI values from very low to low (Figure 16d,e).This clearly indicates that structural elements are an important component differentiating vulnerability.Taxable property, structural vulnerability and the parameter use and land cover, belonging to the morphological component of vulnerability, clearly delineate the composite vulnerability of industrial areas.
Most residential areas (73%) had a CVI ranging from very low and moderate.In these areas, the CVI value is strongly influenced by structural vulnerability, with the highest values coinciding with building sites (Figure 16b,c).In addiotion, noteworthy is the influence of the social and tax elements on the final value of CVI in residential areas.
Clearly, the intrinsic characteristics of the region and the types of activities carried out there influence the composite vulnerability of the territory.In the analyzed areas, morphology decisively influences the vulnerability of seaside areas, whereas residential areas are strongly affected by structural, tax and social vulnerability.On the other hand, the vulnerability of industrial areas is largely governed by the tax and structural vulnerability of the region.
With regard to buildings potentially affected by a tsunami, the calculation of the structural vulnerability and the analysis of typology strengthened the overall results, allowing us to identify and highlight differences between the study areas (Table 8).For example, 51% of potentially affected buildings located in the industrial areas exhibit high and very high vulnerabilities.In beach-use areas this value is 45%, and in residential areas it is 25%.The high value observed for industrial areas can be explained by the dominant building types (Service/Industry and Outbuilding/warehouse) because many of these buildings are prefabricated (including buildings to support industrial activities) and many also have significant structural limitations.Conversely, the highest percentage of buildings classified with very low and low vulnerabilities are located in residential areas (34%) followed by industrial areas (20%) and seaside areas (18%).
A typological analysis of potentially affected buildings allowed us to identify functional differences in terms of the building of the analyzed areas.In the beach-use areas (69%) and industrial areas (55%), Services/Industry and outbuilding/warehouse buildings dominate the landscape.In residential areas, mixed-use buildings dominate, with 37% of the 55% of buildings designated as mixed-used categorized as not exclusively residential.In terms of exclusively residential buildings, the highest percent of affected buildings are located in residential areas (30%).In industrial areas, exclusively residential buildings account for 25% of all buildings and most of these buildings are part of informal settlements built near old industrial complexes (e.g., Quinta da Parvoíce and Vila Maria).As a result, these buildings have serious structural constraints that contribute to the high structural vulnerability values of the area.In beach-use areas, this value is 14%.Finally, potentially affected critical infrastructure was also identified and analyzed.In residential areas, critical infrastructure accounts for 4% of all potentially affected buildings.For example, in the urban area of Setúbal, a set of strategic buildings associated with city management and crisis response (e.g., police, firefighters, schools, custom house, and captaincy) was identified.In industrial areas, 4% of potentially affected buildings were related to critical infrastructure, primarily power transformer stations.potentially be hit and destroyed, and the port activities and estuarine agriculture would be affected.A potential tsunami inundation impact in this small-sized city would cause issues in societal function, namely related with access to the hospital, which is located in inundation area proximity.
The results also allow the adoption of risk prevention and mitigation measures related to the impact of tsunamis.The vulnerability outputs enable the selection of planning strategies related to the land-use and occupation, defining areas with several occupation restrictions or removal from affected areas any infrastructure critical for society to function.Still in the planning framework processes, the results can support specific regulation for building, according to the vulnerability zoning.The social vulnerability component of data supports a warning and alert system adjusted for different risk groups and establishes forms of communication tailored to the sectorial activities that take place in the areas directly affected by the flood or the surrounding areas [41][42][43].
With the knowledge of vulnerability resulting from urban occupation, type of occupation of the buildings and the characteristics of users, civil protection resources and emergency procedures, namely through the plan management, the warning capabilities and alerts, evacuation support and rescue ranking, as mention by different authors, can be adjusted and scaled [38,44,45].
The introduction of this composite methodology allows for a more precise idea about a tsunami's local impact and the major societal disruptions enable an integrated view of municipal, sites and buildings data.These made possible a management and planning process by identifying areas of potentially affected critical infrastructure and urban fabric, allowing timely intervention in these areas to reduce vulnerabilities or plan for their relocation outside of the danger area.This local approach is more recognized by land use planning or emergency municipal managers [46,47], in result of the distinctive display in the territory of the impacts resulting from a tsunami inundation.Based on these results, it is determined that authorities in the Setúbal study area should carry out of preliminary training exercises and consider relocating the firehouse outside of the impact area.
This multidimensional methodology can be applied to different geographical contexts (seaside, residential and industrial areas) and enables the identification of differences within and among regions.Finally, the results may also assist critical discussion for scientists and managers about methodologically complementary of territorial vulnerability, for example earthquakes and tsunamis on coastal areas.

Figure 1 .
Figure 1.Location of the two study areas: (a) framework for Portugal; (b) municipality of Figueira da Foz; (c) municipality of Setúbal; and (d) municipality of Vila do Bispo.

3
Calculation and cartographic representation of the vulnerability indices for the four analysis factors (morphological, structural, social and tax characteristics).4 Calculation and cartographic representation of a Composite Vulnerability Index (CVI) using different classification levels.5 A comparison of the tsunami vulnerability classification for the different studied areas and their relation to territorial features.

Figure 6 .
Figure 6.Industrial zone of Setúbal study area: (a) Satellite image from study area; (b) Tsunami inundation from numerical model results for the industrial zone of Setúbal [2].

Figure 8 .
Figure 8. Graphical process of taxable property vulnerability calculation.

Figure 11 .
Figure 11.Graphical process of Composite Vulnerability Index calculation.

Table 5 .
Percentage of area by classes of vulnerability.

Table 6 .
Percentage and number of buildings by classes of vulnerability.