Next Article in Journal
An Assessment of Hydropeaking Metrics of a Large-Sized Hydropower Plant Operating in a Lowland River, Lithuania
Next Article in Special Issue
Seven Good Reasons for Integrating Terrestrial and Marine Spatial Datasets in Changing Environments
Previous Article in Journal
CFD Modelling of Particle-Driven Gravity Currents in Reservoirs
Previous Article in Special Issue
Sound Velocity in a Thin Shallowly Submerged Terrestrial-Marine Quaternary Succession (Northern Adriatic Sea)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coastal Vulnerability Assessment along the North-Eastern Sector of Gozo Island (Malta, Mediterranean Sea)

1
REgional Models and geo-Hydrological Impacts (REMHI Division), Centro Euro-Mediterraneo sui Cambiamenti Climatici (CMCC), 73100 Lecce, Italy
2
Department of Chemical and Geological Sciences, University of Modena and Reggio Emilia, 41125 Modena, Italy
3
Research and Planning Section, Marine and Storm Water Unit, Public Works Department, FRN 1700 Floriana, Malta
4
Euro-Mediterranean Centre on Insular Coastal Dynamics (ICoD), University of Malta, MSD 2080 Msida, Malta
5
Institute of Earth Systems, University of Malta, MSD 2080 Msida, Malta
*
Author to whom correspondence should be addressed.
Water 2020, 12(5), 1405; https://doi.org/10.3390/w12051405
Submission received: 2 April 2020 / Revised: 11 May 2020 / Accepted: 12 May 2020 / Published: 15 May 2020
(This article belongs to the Special Issue Landscapes and Landforms of Terrestrial and Marine Areas)

Abstract

:
The coastal landscape of the Maltese Islands is the result of long-term evolution, influenced by tectonics, geomorphological processes, and sea level oscillations. Due to their geological setting, the islands are particularly prone to marine-related and gravity-induced processes, exacerbated by climate change. This study aligns different concepts into a relatively concise and expedient methodology for overall coastal vulnerability assessment, taking the NE sector of Gozo Island as a test case. Geomorphological investigation, integrated with analysis of marine geophysical data, enabled characterization of coastal dynamics, identifying this stretch of coast as being potentially hazardous. The study area features a high economic value derived from tourist and mining activities and natural protected areas, that altogether not only make coastal vulnerability a major concern but also the task of assessing it complex. Before introducing the methodology proposed for overall vulnerability assessment, an in-depth revision of the vulnerability concept is provided. The evaluation was carried out by using a set of key indicators related to local land use, anthropic and natural assets, economic activities, and social issues. Results show that the most critical areas are located east of Marsalforn including Ramla Bay, an important tourist attraction hosting the largest sandy beach in Gozo. The method combines physical exposure and social vulnerability into an overall index. It proves to be cost effective in data management and processing and is suitable for the identification and assessment of overall vulnerability of coastal areas to consequences of climate- and marine-related processes, such as coastal erosion, landslides and sea level rise.

1. Introduction

Coastal zones are host to very dynamic and complex environmental systems, subject to the direct and indirect influence of a number of factors that have contributed to their evolution over time. The present-day landscape of coastal areas is the combined result of interactions between natural agents that include physical factors inherent in the system and external climatic and marine forces [1,2]. Human activity also plays an important role in shaping coastal dynamics, often exerting additional pressures that may dominate over natural processes [3]. During the past decades, the rate of human occupation along littoral areas has risen significantly [4] and currently, in Europe, about 86 million people are estimated to live within 10 km from the coastline [5]. Furthermore, approximately one-third of the Mediterranean population is concentrated along its coastal regions and around 120 million inhabitants are concentrated in coastal hydrological basins located in the southern region of the Mediterranean Sea [6]. Coastal areas are the transitional zone between the aquatic and the terrestrial ecosystems and they have an environmental intrinsic value on account of their high level of biological diversity, which supports the provision of several ecosystem services essential for human well-being [7,8]. In view of these considerations, the sustainable conservation of coastal areas is a worldwide issue and coastal vulnerability evaluation and risk assessment are of paramount importance for integrated coastal management [9].
Coastal hazards, including marine-related and gravity-induced processes such as landslides, coastal erosion, storm water runoff and coastal flooding, have different impacts on coastal values, due to differing local geomorphological and anthropogenic settings. The impact of these processes is expected to increase with climate change. The assessment reports of the Intergovernmental Panel on Climate Change [10,11,12,13] emphatically forecast increase in both the frequency and the severity of extreme weather/climate events, combined with sea level rise, which altogether will undoubtedly impact coastal areas more adversely.
The Sendai Framework for Disaster Risk Reduction 2015–2030 [14], whose implementation is overseen by the United Nations Office for Disaster Risk Reduction (UNDRR), has recognized risk assessment as an important preliminary action on the basis of its relevance internationally. This points towards the prioritisation of the role of “understanding disaster risk” and “enhancing disaster preparedness” and promotes the use of “foresight” and “scenarios” for improving the level of preparedness for existing, emerging and new types of risk. At the same time, the identification of suitable climate adaptation strategies, which also account for future scenarios, is essential for increasing the resilience of coastal areas and ensuring the long-term conservation of coastal natural services and anthropic activities [15].
At the European level, risk assessment and mapping have been included in a number of legislative instruments, such as the Floods Directive [16], which requires member states to provide maps and indications for risk management, prescribing specific requirements for climate change impact evaluation.
As highlighted in the Special Report on “managing the risks of extreme events and disasters to advance climate change adaptation” published by the Intergovernmental Panel on Climate Change (IPCC) [17], vulnerability has been a pivotal concept for disaster risk reduction since the 1970s. Initial studies in this field, carried out by Baird et al. [18], O’Keefe et al. [19], Lewis [20], and Hewitt [21], introduced the concept of disaster risk as the combination of the probability of occurrence of a hazardous event with its negative consequences. Over time, the concept of vulnerability has been interpreted in different ways [22,23,24,25] reflecting how the approaches and conceptual models for its assessment differ among the scientific research communities [23,25].
This paper proposes and applies a methodological approach that integrates physical exposure and social vulnerability into an “overall vulnerability index”, to identify areas that can be negatively affected by climate- and marine-related processes, such as coastal erosion, landslides and sea level rise.

2. Conceptual Framework

Given the spectrum of possible conceptual interpretations of “vulnerability” and other related terms such as “exposure”, it is deemed essential that, as of primacy, and prior to proposing the methodology adopted in this research for the overall assessment, the range of most-commonly used concepts of vulnerability and exposure coined in the contexts of both climate change and disaster risk management are first introduced here, before clarifying which concept and definition are adopted by this study. The disaster risk community identifies the notion of vulnerability as one of the elements required for risk assessment [26]. In this scientific context, “risk” is defined by exposure to a hazard, which is a potentially damaging physical event or phenomenon, and vulnerability, which denotes the relationship between the severity of the hazard and the degree of damage caused to the exposed element [27]. Specifically, according to the United Nations International Strategy for Disaster Reduction [26], vulnerability is defined in a qualitative way as “the characteristics and circumstances of a community, system or asset that makes it susceptible to the damaging effects of a hazard”. In the same document, exposure is also defined qualitatively as “people, property, systems, or other elements present in hazard zones that are thereby subject to potential losses”.
In earlier definitions, vulnerability was considered as the degree of loss of an element at risk after the occurrence of a natural process [28]. On one hand, the above-mentioned and more recent (but not the latest) definitions consider the notion of vulnerability as denoting a pre-existing condition related to the characteristics of the elements at risk and gives less emphasis to the process [29]. On the other hand, in the Climate Change Adaptation (CCA) context, more focus is placed on the concept of vulnerability defined as a function of exposure, sensitivity and adaptive capacity [10,30]. Therefore, in the risk context, the definition of vulnerability is quite similar to that describing the sensitivity of the system’s components in the climate approach [10], while in the climate change community it is defined in a similar way to the concept of risk used in the Disaster Risk Reduction (DRR) context [31].
Only in more recent years have the United Nation Institutions (UNDRR, IPCC) had a key role in converging on a common and shared definition of the vulnerability concept in CCA and DRR fields. In fact, the IPCC integrates the different conceptualizations of vulnerability and provides its upgraded and clearer definition in the glossary of the Fifth Assessment Report (AR-5, [11]), in which vulnerability expresses “the propensity or predisposition to be adversely affected. Such predisposition constitutes an internal characteristic of the affected elements (or societies) and includes the concepts to cope with, resist, and recover and the lack of capacity to cope and adapt to the adverse effects of a physical event”. Indeed, in this latter AR, the definition of risk (as a whole) is expressed as being the result of the interaction between vulnerability, exposure and hazard, whereby the distinction between the contribution of vulnerability and exposure to risk has been made clearer, which is a view also shared with the DRR community. This paper adopts a method of combining “physical exposure” information and “social vulnerability” data, which is in line with this evolution and convergence of concepts of vulnerability and exposure, across the above discourses.
In parallel to the vulnerability concept development, several actions taken in the last couple of decades have also focused on the development of methods and tools for supporting decision-makers in the reduction of coastal hazards’ impacts. Particularly noteworthy is the index-based approach, which was introduced by Gornitz at the beginning of the 1990s [32,33]. This is mainly based on indirect analysis supported by photointerpretation and topographic maps and is still the most commonly used method for coastal vulnerability assessment at regional levels [34,35,36,37,38,39,40].
At the European level, a number of transnational research projects have been funded with the aim of improving scientific knowledge and increasing the awareness of decision-makers on coastal vulnerability, the potential damages of marine processes and the effectiveness of coastal adaptation strategies (such as EUROSION, MICORE, RISC-KIT, ANYWHERE, OPERANDUM). In particular, the main outputs of the RISC-KIT Project are represented by a set of tools that allow identification of the most vulnerable areas (hot-spots) by means of a coastal index [41,42] and then selection of suitable measures for increasing coastal adaptation and therefore favouring risk reduction [43].
The evolution and emergence of concepts of vulnerability and exposure as being different features is warranted in certain situations, to provide a better understanding of “why” and “how” societal assets and values may be under threat or at risk under different scenarios of exposure. This can provide guidance to planners and policy makers in identifying suitable adaptation measures and actions, but it also makes risk assessment far more complex.
Some factors that contribute to the vulnerability of society can be expressed in terms of the dependence of societal wellbeing on certain exposed physical assets or other values that are at risk. Examples of this are infrastructure and utilities, which underpin societal wellbeing to different degrees, depending on their nature. This component of vulnerability is distinct from vulnerability associated with the inherent characteristics of society. A community relying on physical elements for its wellbeing can be deemed as being vulnerable to the extent (or level) of its dependence on those specific physical elements located in a hazard-prone area. This creates different levels of “physical exposure”, in direct relation to their level of importance and to the degree to which societal wellbeing is dependent on and would be affected by losses to them. Thus, “physical vulnerability” is considered as one part of “overall vulnerability”, while “social vulnerability” (attributed to inherent societal factors) is considered as another component of “overall vulnerability”.
In this context, the research presented here aims to evaluate the overall vulnerability, and to refine the analysis and understanding of two distinct components, physical vulnerability (representing exposure from dependence on physical elements) and social vulnerability to a given set of external climate- and marine-related processes. Mirroring physical exposure as physical vulnerability and combining it with social vulnerability into an overall vulnerability index, an expedient and cost-effective method for the identification (and assessment) of the overall vulnerability of coastal areas at risk is provided.
The selected coastal study area, located along the NE sector of the Island of Gozo (Maltese archipelago), is one for which considerable applicative research has already been undertaken in order to showcase the high geological and geomorphological significance of its coastal landscapes [44]. Although the area is known to be particularly susceptible to several coastal hazards (erosion, storm surges, floods, sea level rise, landslides), it still lacks a detailed and refined appraisal of its natural and anthropic coastal elements, based on an analysis of assets and values that could potentially be at risk. Approaching the assessment of vulnerability of this coastal area from the perspective of the two components of social vulnerability and physical vulnerability provides a sharper analytical insight into overall vulnerability. Furthermore, the analytical methods and tools used in this study provide a methodological template for overall vulnerability assessment in the form of an Overall Vulnerability Index (OVI) that can be computed and used, with relative ease at different scales.
The scientific literature review carried out as part of this study reveals a gap in the area of assessment of vulnerability. Studies carried out to date have focused on the overall evaluation of hazard and susceptibility along the Maltese coastal sectors [45], rather than on the methodological assessment of vulnerability per se. Integrated approaches have been applied along the NW coast of Malta for landslide hazard assessment [46,47,48,49]. Landslide susceptibility assessment assisted by Persistent Scatterers Interferometry (PSI) was carried along the same stretch of the NW coast of Malta [49,50]. Based on the assumption that the identification and analytical mapping of the exposed elements represents a key step for the evaluation of the coastal risk, this study tries to bridge this gap by proposing a relatively simple yet reliable, cost-effective and easily replicable procedure for the assessment of coastal vulnerability. Nonetheless, while not being overly complex, the proposed approach provides an analytical discernment of different components of overall vulnerability that, when mapped geographically, provide sufficient guidance as to which areas are most vulnerable and on what account.
Following on from the methodological approaches proposed in previous studies [42,51,52,53,54], this study formulates and applies a research method for overall vulnerability assessment, which is based on the following steps: (i) identification of the main exposed elements (natural and anthropic) located in the investigated area; (ii) definition of their relative exposure level, in economic and ecological terms; (iii) assessment of the social vulnerability of the population living in the investigated area; (iv) calculation of the overall vulnerability by means of a combined index. As highlighted in [55] the use of an index as an evaluation tool requires the definition, weighting and aggregation of a number of indicators, which are defined as variables that are “an operational representation of a characteristic or quality of a system” [56,57].
In detail, the method applied here foresees a set of indicators for the evaluation of the exposure level of local land uses (considering both the presence of economic activities as well as natural protected areas) and anthropic assets (such as transport networks and main utilities). Furthermore, the social context that characterizes the investigated area is also taken into account, in order to include social vulnerability in terms of population capacity to respond to and cope with a hazardous event in the overall evaluation of vulnerability. The combination of physical and social vulnerability provides the overall index, which expresses the level of overall vulnerability.
The use of distinct physical and social vulnerability indicators, and the definition of an Overall Vulnerability Index (OVI), has the main advantage of summarizing complex issues making them more easily interpretable, facilitating decision making and communication among stakeholders [58]. Meanwhile, the distinction between the two components at analysis stage provides a more in-depth understanding of specific anthropic elements that contribute to different risk levels.

3. Study Area

3.1. Geological and Geomorphological Setting

The Maltese archipelago is located in the central Mediterranean Sea and comprises the main islands of Malta, Gozo and Comino (Figure 1). The coastal landscape (Figure 2) is the result of long-term evolution under the influence of tectonic activity, geomorphological processes and sea level oscillations. Due to their geological and geomorphological setting, these islands are particularly prone to different marine-related and gravity-induced processes such as landslides, coastal erosion, storm water runoff and sea level rise, enhanced by ongoing climate change. Multidisciplinary research and integrated investigations have been carried out to better understand the evolution of the geomorphological processes within the archipelago in the wider dynamic scenario of ongoing climate change as a way towards the reduction of risks associated with these processes [47,59,60,61].
Previous studies of the Island of Gozo have focused on general aspects of coastal features [60,63,64,65,66]. Only a few papers deal with the specific geomorphological aspects of the island [67,68,69], some of them referring to its rich geoheritage [70,71,72,73]. A detailed geomorphological map of the investigated stretch of coast (NE Gozo), based on geomorphological and geological field surveys integrated with the analysis of marine geophysical data, has been published recently [44]. The geomorphological map is accompanied by two other maps that show land use and the distribution of coastal geomorphotypes (Figure 3). Such documents constituted the basis for carrying out the vulnerability analysis presented below.
From a geomorphological perspective, the investigated coastal stretch is characterized by limestone plateaus bounded by steep structural scarps which are progressively reshaped by gravitational and/or degradation processes. Clayey slopes, located at the foot of the limestone plateau, accommodate terraced fields of actively used or abandoned agricultural land. Numerous blocks of rock are strewn over the clayey terrain (a unique landscape known as rdum in Maltese) that slopes more gently away from the plateau edge.
The investigated coastline is characterized by the alternation of inlets and promontories. The accumulation of sand and mixed grainsize deposits results in the formation of pocket beaches where this corresponds with bays and coves. The large sandy beach of Ramla il-Ħamra Bay (‘red sandy beach’) is a particularly noteworthy example. It is also partly bounded by dunes on the landward side and is protected as a Special Area of Conservation (SAC) within the Natura 2000 network. The main inhabited centre of the investigated stretch of coast is Marsalforn, with its homonym bay, the latter being intensely urbanised and anthropized, comprising a predominantly built-up coastline.
The landscape of the investigated stretch of coast, and in general of the entire Maltese archipelago, is largely controlled by the different erodibility of the exposed lithostratigraphic units constituted by marine limestones and marls of the late Oligocene-Miocene [74,75]. The outcropping geological formations comprise (from the oldest to the youngest): Lower Coralline Limestone Formation (late Oligocene, Chattian), Globigerina Limestone Formation (late Oligocene—middle Miocene, late Chattian—Langhian), Blue Clay Formation (middle-late Miocene, Serravallian—Tortonian) and Upper Coralline Limestone Formation (late Miocene, Tortonian—early Messinian).
The Lower Coralline Limestone Formation consists of pale grey, hard, shallow marine biomicrites and biosparites [76], and outcrops in a restricted coastal stretch forming subvertical cliffs. Few Lower Coralline plunging cliffs (sensu [77]) are found in the investigated stretch of coast. They host roof notches with an asymmetric shape, recognized by Furlani et al. [78]. Lower Coralline Limestone is more commonly found in sloping coast formations, as typified in the stretch between Dahlet Qorrot Bay and Ras il-Qala.
The Lower Coralline Limestone Formation underlies the Globigerina Limestone Formation, the latter being younger and more erodible with respect to the former. The Globigerina Limestone Formation consists of flattened areas along the coast, and includes yellowish, fine-grained, planktonic foraminiferal limestones. This formation features prominently as shore platforms along the investigated stretch of coast, where it outcrops above sea level and where the overlying softer blue clay layer has been eroded away, with examples at the Rdum il-Kbir promontory and on the eastern side of Marsalforn Bay. The Blue Clay Formation, overlying the Globigerina Limestone Formation, consists of grey, soft marls, clays and silty sands, forming gentle slopes. Sea cliffs shaped in Blue Clay can be also found within the study area. Finally, the youngest unit outcropping in the study area is the Upper Coralline Limestone Formation. This layer forms the plateaus at the top and is frequently weathered into steep cliffs and well-developed karst topography. It is also the source of the blocks of rock strewn onto the clayey slopes that slope away from the plateau edge more gently. A schematic stratigraphic coloumn is shown in Figure 1.
The NE coast of Gozo is particularly prone to gravity-induced processes such as rock spreads, block slides, rock falls and earth flows/slides. This is mainly due to the tectonic and geological features of the area which is characterized by a dense network of joints and fractures [79] and by the superposition of terrains with different geomechanical behaviour (cf. [80]). In fact, the brittle limestone plateaus overlying clayey terrains enhance the fracturing of the plateaus and the development of lateral spreading locally evolving into block sliding [44,46,81,82,83,84,85].
Moreover, clayey slopes are more prone to shallow earth flows and earth slides [86], while limestone cliffs are affected by rock falls, which has caused scarp retrogression over time.
Geomorphological investigations, relevant to the submerged area along the investigated stretch of coast [44], have revealed that block slides and earth flow/slide runout continue locally below the sea level, reaching ca. −20 m of depth. This is also in agreement with evidence from other submerged areas along the Maltese Islands [44,87]. Evidence from landslide dating in the NW coast of Malta infer that these deposits were emplaced in a subaerial environment during sea level lowstands and subsequently became submerged during the post-glacial marine transgression [61].

3.2. Social, Economic and Tourist Setting

Studies that include geoheritage assessment and geosite inventory highlight that the Maltese archipelago is considered as an attractive geotourist destination due to the strong interaction between natural and cultural aspects (cf. [71,88]). Data from the World Travel and Tourism Council (WTTC) for the year 2017, retrieved from Selmi et al. [88], show that 27.1% of Malta’s GDP and 28.3% of total Maltese employment (corresponding to 55,000 jobs), were accounted for by activities directly related to and induced by travel and tourism [89,90]. Taking the type of tourism into consideration, data referring to 2017 show that the majority of tourists (almost 85%) visited Malta for holidays, while a very low percentage visited the country for business and other purposes (8% and 7%, respectively). More recent data show that the number of tourists is constantly increasing (standing at 2.6 million in 2018). The WTTC forecast that the travel and tourism contribution to national GDP will rise to 34.6% by 2027.
The tourism on the Island of Gozo is highly dependent on tourism activity of the main island, Malta [91], given that all of Gozo’s tourist traffic necessarily passes through mainland Malta. Tourism is a significant source of income and employment and it is one of the primary contributors to the Gozo economy. Gozo’s tourism also relies heavily on domestic tourists (with 400,000 domestic tourists or visitors per year coming from Malta) and on one-day trips of international tourists [92]. In 2018, almost 100,000 guests including resident and non-resident spent an average of three to four nights in one of the accommodation facilities of Gozo [93].
The Gozo Island attracts many tourists, especially during the summer period, for its environmental, cultural and geological heritage. Moreover, senior and middle-aged foreign residents more commonly chose the Island of Gozo as a place for their retirement, mainly for its relative peacefulness and quiet, combined with its array of scenic features (Gozo is colloquially known as “the place where time stood still”).
The coastal sector investigated by this study includes four administrative districts: Żebbuġ (7.6 km2), Xagħra (7.6 km2), Nadur (7.2 km2) and Qala (5.9 km2), which are characterized by extensive urban development. The study area is host to two renowned tourist destinations, Marsalforn Bay and Ramla Bay, both of high importance to Gozo tourism. Marsalforn Bay is one of the main inhabited centres located on the Gozo coast. The availability of a significant number of accommodation facilities, shops, restaurants and several diving centres contribute to a dense tourist population and to lucrative businesses in Marsalforn, especially in summer.
Total guests in the Gozo and Comino region increased by 12.6% to 97,781 in 2017, while total nights spent went up by 11.1% to 347,943 when compared to the previous year [94]. Within a continuous upward trend in the last five years or so, the Gozo and Comino region recorded a strong growth in terms of inbound tourist arrivals in 2017 [94]. It is worth noting that among the top five localities where inbound tourists to Gozo stayed longest [94], two are within the study area, with important tourist destinations such as Marsalform and Ramla Bay.
Ramla Bay is the largest sandy beach in Gozo. It is characterized by golden-reddish sand, which lends it its name in Maltese (Ramla il-Ħamra) and makes this beach peculiar and unique to the Maltese Islands. The area around the beach is also of archaeological and historical interest, hosting Roman remains lying beneath the sand, a submerged seawall, fugass and battery defensive structures constructed by the Knights of St. John of Jerusalem in the mid-18th century, and the famous Calypso Cave looking over the western side of the beach.
Furthermore, the study area as a whole includes areas of high natural and ecological importance, hosting two wide natural protected zones, Għajn Barrani (located west of Marsalforn, including Ramla Bay) and Il-Qortin tal-Magun u l-Qortin il-Kbir (located close to Dahlet Qorrot Bay), both included in the Natura 2000 network as SAC. Finally, it should be noted that the high economic value of the study area is also derived from specific land uses in certain areas, such as quarrying, which occupy a surface of 0.13 km2.

4. Materials and Methods

The method proposed and applied in this research for the assessment of coastal vulnerability is in line with the most recent index-based approaches generally used for this kind of analysis. In detail, the method relies on the outcomes of other research in the field [41,95] and is based on a new approach for the evaluation of indicators referring to: (i) the potentially exposed categories of assets (defined as “physical indicators”); and (ii) a number of parameters related to the social context (defined as “social indicators”).
The indicators are interpreted here as operational representations (cf. [56,57]) of the physical and social characteristics of the area. Each physical indicator comprises different categories of assets, to which a score representing an increasing level of exposure was attributed. In order to tailor the approach to the local settings of the study area, the exposure level assigned to each category of assets was defined based on expert judgment. It is a widely shared opinion that the scoring and aggregation of indicators into indices may have a large impact on the resulting rankings and, consequently, on decision-making [55].
The spatial overlay of the physical indicators provides an estimate of the physical vulnerability level, which ranges from very low to very high vulnerability. Meanwhile, the various social indicators refer to vulnerability of the socio-economic aspects. These are also classified into five levels ranging from very low to very high, to provide a measure of the social vulnerability level for the investigated area. The overlay of the two sets of aggregated physical and social indicators enables an overall zonation that shows grades of the combined physical-socioeconomic vulnerability, defined as “overall vulnerability”, thus identifying which are the most vulnerable stretches of the investigated coastal area.
Specifically, the method applied here comprises the following four main steps:
  • Definition of the landward limit of the coastal area to be investigated;
  • Classification of physical and social indicators into levels;
  • Data overlay and computation of an Overall Vulnerability Index;
  • Overall coastal vulnerability zonation and representation on a map.
The analyses are supported by GIS tools, which allow identification of the exposed coastal assets (i.e., natural and semi-natural environments, buildings, infrastructure, and agriculture), their combination with social data, and the calculation and mapping out of the overall vulnerability levels.

4.1. Definition of the Landward Limit of the Coastal Area to be Investigated

The area investigated was delimited according to the definition of a RICE area (Radius of Influence of Coastal Erosion and Flooding) proposed in the framework of the EUROSION project [95] for the identification of the coastal areas potentially impacted by marine-related process. The limit was set at a maximum distance of 500 m inland, or reaching a maximum of 5 m a.s.l. from the coastline. Furthermore, a buffer area of 100 m inland from the edge of the scarps formed in the Upper Coralline Limestone—which are extensively affected by landslide processes—was also added to the area defined in the above manner, in order to specifically account for landslide hazard. Information about landslide distribution provided by the geomorphological map of the study area (cf. [44]) was used specifically for this purpose, so as to identify the areas prone to slope instability.
The resulting landward limit was simplified and shown as a dashed red line in Figure 3, in order to provide a more linear indication of the investigated coastal sector.

4.2. Classification of Physical and Social Indicators

This step comprises the identification of data required for the evaluation and classification into five different levels, from very low to very high, for each of the two sets of indicators (physical and social indicators).
The evaluation of the anthropic and natural assets potentially exposed to coastal hazards is based on the elaboration of the following indicators: land use, transport network, and utilities. A GIS layer has to be created for each indicator in order to identify polygons representing the spatial unit to which the related exposure level is assigned, indicating its specific numerical value, ranging from 1 (very low exposure) to 5 (very high exposure). These layers are first converted to a raster format (5 × 5 m) and then overlaid to estimate the physical vulnerability level, by considering the maximum exposure value assigned to each cell.
At the detailed level, the land use information was collected from the land use map available in Prampolini et al. [44], in which land use categories were classified as: (i) natural and semi-natural areas, (ii) agricultural areas, and (iii) artificial surfaces. This land use map was supplemented with additional detailed information in order to include strategic elements, such as civil protection posts, police stations, emergency posts, fire corps posts, port authorities’ posts, hospitals, and schools. The identification of these elements was supported by photointerpretation of the most accurate maps available for the study area and verified in the field during the investigation.
An exposure level (ranging from 1, very low exposure, to 5, very high exposure) was assigned to each land use category based on expert judgement. The highest exposure level was attributed to the areas occupied by settlements, constructions, and human activities while the lowest level was assigned to abandoned agricultural zones. Details concerning the land use indicator and the exposure levels related to each category are shown in Table 1.
The transport network information was retrieved from Open Street Map downloaded from [96] as a shapefile. In this case, the roads were classified as follows: primary road, secondary road, tertiary road, inhabited street, residential, services, footway, path, track, and steps; the greater exposure level was attributed to roads of national or international importance. For each polyline element, a buffer distance was built (up to 20 m of total width for the primary road) in order to take into account a proper area of pertinence. Details concerning the exposure levels defined for the transport indicator are shown in Table 1.
The electricity network was taken into account for the utilities. The spatial distribution of these elements is available on Malta Inspire Geoportal [97]. In this case, a buffer zone up to 10 m width (10 m for substations, feeder pillars, and overhead lines, and 5 m for street lighting) was considered for each element. Details concerning the exposure levels associated with each electricity element are indicated in Table 1.
The social indicators allow characterization of the districts located in the investigated coastal area by evaluating, directly and indirectly, the social characteristics of the population living in the zone and therefore prone to be affected by coastal hazards. Generally, the social vulnerability information is obtained from census data. In this study, social data was downloaded from the Inspire Geoportal [97] and integrated with information available in the special report recently published by the Malta National Statistics Office (NSO) [98]. The social indicators are gauged from the economic budget allocated by the Government for supporting the population. The social vulnerability analysis carried out here is based on the assumption that higher allocation of social schemes corresponds to higher vulnerability of the population living in the investigated districts.
Specifically, the following indicators were taken into account: health care, disability, old age, children, and unemployment. These aspects represent some of the higher-risk groups of society, which social services and budget allocations aim to protect. The NSO provides a classification of these indicators in five classes that were here converted in a numerical value ranging from 1 to 5. The available data is provided by the NSO at district level and thus a numerical value was assigned to each district as a whole. The social vulnerability was then calculated as geometric mean of the values attributed to each indicator and then classified into five levels (from 1, for very low vulnerability, to 5, for very high vulnerability) by means of an equal interval classification method. A detailed description of the social indicators used in this study is reported in Appendix A.
The physical and social indicators identified are weighted equally, as in most of the indicator-based studies [99], meaning that each individual indicator has the same influence on the final calculation of the overall vulnerability.

4.3. Data Overlay and Computation of an Overall Vulnerability Index

Once all the required data related to the physical and social vulnerability are collected and expressed in five levels, they are overlaid by means of a specific GIS tool. The overall vulnerability calculation is therefore the result of the aggregation of the physical vulnerability and social vulnerability data, the combination of which provides the definition of the Overall Vulnerability Index (OVI) as follows:
Overall Vulnerability Index = (Physical vulnerability × Social vulnerability) ^ 0.5
Finally, the Overall Vulnerability Index is classified into five levels of increasing vulnerability (from very low to very high vulnerability) by means of the equal interval classification method.

5. Results

The results of this study are represented by: (i) GIS-based maps that show the spatial distribution of the natural and anthropic exposed elements and their classification in terms of physical vulnerability; (ii) social vulnerability level estimated for each investigated coastal district; (iii) overall vulnerability map of the NE part of the Gozo Island; (iv) areal extent and relative percentage of each vulnerability level. The land use classification (Figure 4) enables the evaluation of the surface occupied by different categories of natural and anthropic elements (cf. Table 1 in Section 4.2).
In detail, 3.6 km2 are occupied by agricultural areas, 1.2 km2 by abandoned agricultural areas, 1 km2 by terraced agricultural field, 1.2 km2 by natural protected areas, 0.8 km2 by bare rock, 0.7 km2 by natural and vegetated areas, and, 0.5 km2 by residential areas. Finally, quarrying and land occupied mainly by agriculture with significant areas of natural vegetation occupy 0.1 km2. The remaining surface area (corresponding to 602,268.2 m2) is occupied by the following categories: strategic elements (including a Police station), historical/archaeological sites (The Tower of Ta Sopu, west of Daħlet Qorrot Bay, and Saint Anthony’s Battery at Ras il-Qala), green urban areas, beaches, dunes, sand, abandoned terraced fields, entertainment zones, and a cemetery. These values are indicated as percentages in Figure 5.
The interpretation and spatial representation of the overlaid physical indicators (cf. Section 4.2) (cf. Table 1) is representative of the different levels of physical vulnerability, given that the physical exposure levels are correlated directly to the level of value of the assets at risk.
This enables the zonation and calculation of what are, in effect, the extent of the areas with different levels of physical vulnerability, expressed in square kilometres and in percentages of the total surface investigated. The spatial interpretation of physical vulnerability is mapped in Figure 6 while the corresponding numerical values are reported in Table 2.
The numerical values assigned to each of the social vulnerability indicators for the investigated coastal districts are summarised in Table 3.
The aggregation of the social vulnerability indicators (cf. Table 3) enables the evaluation of the social vulnerability level for each of the investigated districts (Figure 7). Specifically, two districts (Nadur and Xagħra) are characterized by medium vulnerability and two (Qala and Żebbuġ) by high vulnerability.
Finally, the overlay of the physical vulnerability levels with the social vulnerability levels enables a computation (and mapping) of the Overall Vulnerability Index for different polygons in the study area (cf. Section 4.3). It provides the basis for evaluating the variation in the overall vulnerability levels across the investigated coastal sector and enables this to be represented spatially (Figure 8).
In sum, 1.9 km2 are occupied by areas showing a low overall vulnerability level, 5.6 km2 are occupied by areas showing a medium overall vulnerability level, 0.7 km2 are occupied by areas with a high overall vulnerability level and, finally, areas characterized by a very high overall vulnerability level cover 1 km2 (Table 4).

6. Discussion

The proposed index-based method allows zoning of the investigated coastal stretch into different levels of vulnerability to climate- and marine-related processes. Results show that the study area is divided in four zones only (from low to very high vulnerability), since there are no areas with a very low level of vulnerability.
The coastal sectors located east of Marsalforn Bay, which include Ramla Bay and the area on the western side of Daħlet Qorrot Bay, show the highest overall vulnerability level. This is due to the combination of very high physical vulnerability levels pertaining to the presence of two Natura 2000 sites (cf. Section 3.1 and Section 5) and high social vulnerability levels for the Xagħra and Nadur districts. The high social vulnerability of Xagħra and Nadur is explained by the fact that each of the latter districts accounts for more than 4000 inhabitants, which is twice the number of inhabitants of the two other districts considered, Żebbuġ and Qala, for which a medium social vulnerability is assigned.
The coastal sector surrounding Marsalforn Bay is mainly characterized by a medium level of overall vulnerability as a result of the combination of a medium physical vulnerability level, explainable by the presence of residential areas, and, a medium social vulnerability level obtained for the Żebbuġ district. Furthermore, the innermost sector of the study area shows a prevailing medium overall vulnerability level, resulting from the combination of a low to very low physical exposure level, owing to the presence of cultivated and abandoned agricultural fields respectively, and a high social vulnerability level resulting for the districts of Xagħra and Nadur, as discussed above. The eastern part of the investigated sector, with bare outcrops of rock, hosts areas characterized by a low overall vulnerability.
The index-based method here proposed can be considered a suitable approach for the identification of coastal areas that are most vulnerable to consequences of different climate- and marine-related processes. Since the method adopted here combines vulnerability and exposure, the results represent an overall vulnerability assessment. It relies on the evaluation of the two main components that are used to define the overall vulnerability: physical and the social vulnerability, thus accounting for both the physical assets potentially exposed to climate and marine processes and the social aspects of the local population. This approach is in line with the convergence of different terms by the IPCC in its Fifth Assessment Report (AR-5, [11]).
The application of the proposed OVI method has allowed identification of a number of operational advantages. First of all, as the analysis is based on a sequence of steps, the method is relatively simple and easy to apply. It has been found to be cost-effective on account of the possibility of managing and processing the acquired data in a GIS environment with relative ease. Furthermore, it does not require intense field work, since it relies on indirect analysis, supported by existing geological-geomorphological information that is either readily available in the literature or easily collected. It can therefore be applied to wider coastal stretches, as in the case of other index-based approaches [40,100,101,102,103]. It should be underlined that the applied method can take advantage of data sets and information concerning the natural and anthropic assets which are generally freely accessible and downloadable, e.g., from national geoportals available in most of the countries (at least across Europe), and from open databases, such as OpenStreetMap. Therefore, the method may have a wide range of applicability at different scales. Large scale analyses are based on detailed information about exposed/vulnerable elements and population (such as the land use and vulnerability maps proposed in this study). Meanwhile, small scale vulnerability maps concern wide areas, taking into account only the most prominent and spatially persistent exposed/vulnerable elements and regional data about the exposed/vulnerable population (as in the case of EUROSION Project), which accounted for all the European countries and provided a comprehensive European-level data repository at scale 1:100,000 [95].
The OVI method provides a sound approach suitable for the identification (and assessment) of vulnerable areas and sectors, even as an expedient and cost-effective scoping stage to identify which area may be analysed more specifically and at greater expense. It provides a useful tool for an overall time-efficient and cost-effective approach to focus limited research resources onto where they are needed most.
The presented research is a pilot-study and a first-ever combined overall vulnerability assessment carried out for any area in the Maltese archipelago, and it is promises to provide new contributions at different levels and in different ways. At the localised level, the proposed method, and the results obtained from it, are promising as they reveal the potential applicability of this method for other (and potentially wider) coastal areas of the Maltese archipelago. Moreover, even at the local level, the overall vulnerability assessment could be easily updated by using any newly available or more detailed data regarding the social vulnerability of the investigated districts, as and when this becomes available. The relatively simple combined approach makes it more possible to keep the OVI up to date, even at the localised level. By way of example, it is relevant that during 2016, Malta’s social protection outlay rose by €70.3 million in comparison to 2015 and that Old Age and Sickness/Health Care witnessed the biggest increases in social outlay [98]. The NSO [94] published statistical data on the basis of six different regions, from which localities can be studied specifically as socio-economic parameters change. This means that the social vulnerability level of the investigated district could be expected to rise in the next decades, influencing the overall vulnerability assessment. At the national level, the method could be applied to provide a first-round assessment of overall coastal vulnerability, to identify the most vulnerable stretches of coastal areas and values, and then to follow this up with a more detailed qualitative and quantitative analysis and comparison of risk assessments for specific areas and sectors, for different hazard types.
Running parallel to other research at the European level, further studies could include the spatial susceptibility analysis, including each of the coastal hazardous processes already identified here as affecting the investigated stretch of coast (e.g., erosion, sea level rise, landslides). As already done for other coastal zones in Europe [39,40,54], the investigated area should be classified into zones with different susceptibility, accounting for their proneness to be potentially affected by the specific impacts from extreme events and related processes pertaining to the particular hazard types (e.g., erosion, sea level rise, landslides). The susceptibility to different climate- and marine-related hazards could then be coupled with the vulnerability data derived by this study to perform a complete risk analysis (as done for example in [25,53]).
In this context, it is worth noting that the Maltese archipelago is situated centrally in the Mediterranean Sea, which has been classified as one of the regions most sensitive to climate change and, therefore, it is considered as a hot-spot area [104]. Climate change projections for the Mediterranean region [105] show that the area is experiencing an increasing temperature with consequent change in spatial and temporal distribution of weather/climate extremes [106]. More intense events, with alternating and more severe drought and precipitation, are expected to trigger and exacerbate erosive and mass movement processes [11,13], affecting the spatial distribution of the susceptibility to these types of events and hazards.
A direct consequence of global temperature increase is the rise in sea level, the direct impacts of which on coastal systems become manifest with larger rates of erosion, increase in flooding events, wetland loss, and saltwater intrusion [107]. Studies regarding the reconstruction of sea level changes in the Mediterranean Sea during the last 2000 years have shown that, in tectonically stable Mediterranean areas, the sea rose about 1.1 m [108,109] while for the last two decade the estimated rise was of about 3 cm/decade [110]. However, differential vertical land movements, including uplift and subsidence, characterize the Mediterranean coasts and, for this reason, the trends of sea level rise in the Mediterranean Sea have a large spatial variability [111]. Taking into account the role of terrestrial ice melt, steric effects and glacial isostatic adjustment, the future total Mediterranean averaged sea level rise has been estimated to be between 9.8 and 25.6 cm by 2040–2050, depending on the Representative Concentration Pathways scenario [112]. Sea level projections, obtained by coupling modelled eustatic trends with local ground movements, are in general used for supporting the assessment of the future coastal risk to sea level rise, which is aimed at reducing its impacts on natural coastal environments and human economic activities [3].
Based on the assumption that the investigated coastal area could potentially be affected by all the above mentioned climate-related hazards, further research activities should focus on risk analysis and mapping, as already done with reference to hazards related to sea level rise in several parts of Europe [113,114] and worldwide ([115,116,117,118,119] and reference therein). In this context, the results obtained in this study can prepare the ground for a comprehensive risk assessment, combining the identification of the exposed assets and the evaluation of vulnerability levels with a quantitative assessment of the hazardous processes affecting the area.
Finally, it is worth noting that the study is in line with the methodological approach proposed by the European Environment Agency (EEA) for the identification and implementation of climate change adaptation strategies (European Climate Adaptation Platform Climate-ADAPT). Specifically, in the EEA approach, six steps of analysis are indicated: (i) preparing the ground for adaptation; (ii) assessing risk and vulnerability; (iii) identifying adaptation options; (iv) assessing adaptation options; (v) implementation; (vi) monitoring and evaluation. The vulnerability analysis shown here represents a useful tool for addressing Step 1 and Step 2, while further activities should be planned with the aim of developing and tailoring the most suitable adaptation actions for the protection of the natural ecosystem and the maintenance of the anthropic activities.

7. Conclusions

This research represents a first attempt at the evaluation of coastal overall vulnerability in the Island of Gozo (Maltese archipelago) to the main climate- and marine-related processes, such as coastal erosion, landslides and sea level rise, to which the island is particularly prone. The analysis method developed and tested in this study is based on a conceptual interpretation of overall vulnerability, defined as the combination of two distinct components, namely, physical vulnerability and social vulnerability.
The study departs from a discussion of different notions of risk, including commonly used notions of exposure and vulnerability, which are defined qualitatively in the most authoritative sources (cf. Section 2). Given the identified, often conflicting understanding of the concept of vulnerability by various sectors of researchers and policy makers, it was determined that an in-depth consideration of the vulnerability concept was required prior to determining the vulnerability assessment methodology.
The study proposes a method for the assimilation, analysis and computation of overall vulnerability, and for its graphical spatial representation. The method manages to converge different qualitative notions of physical exposure and social vulnerability, and also makes it possible to derive a spatial quantitative distribution of the Overall Vulnerability Index. The latter can be represented in a map showing different coastal vulnerability levels that are easily readable spatially, making it simple to communicate comprehensibly to decision makers.
The Overall Vulnerability Index is calculated by means of an index-based method that is proposed along the lines of approaches developed in the framework of previous research projects and adapted specifically for the context of the study area and to the available information.
The study area was identified as possessing features of high economic value derived from tourist and mining activities and natural protected areas, altogether making coastal vulnerability a major concern. The final results of the analysis reveal that most of the investigated area (61.3%) is characterized by a medium level of overall vulnerability. A very high overall vulnerability level (11%) was obtained for the areas located east of Marsalforn Bay and close to Daħlet Qorrot Bay, including Ramla Bay, mainly owing to the presence of two sites protected as Special Areas of Conservation within the Natura 2000 network (Għajn Barrani and Il-Qortin tal-Magun u l-Qortin il-Kbir sites). A high vulnerability (7.3%) resulted for the main roads, while 20.4% of the total surface is characterized by low vulnerability areas mainly corresponding to abandoned agricultural fields and bare rocks outcrops.
The type of analysis presented here can be easily replicated for other coastal regions because the data sets required for the proposed method are almost always freely available or relatively easy to compile. Furthermore, the wider application of this method to the entire coastal area of the Gozo Island as well as to the coastal zones of the Maltese Islands has the potential to contribute to an overall risk characterization for the entire Maltese archipelago and it provides a useful tool for the identification of the most exposed and vulnerable zones (hotspot areas) that require action for their protection as a matter of priority. Starting from readily available data, which can be processed and mapped in a GIS environment with relative ease, the method proposed and tested here can provide policy makers, as well as coastal management agencies, with a graphical representation of easily comprehensible and useful data to support decision-making processes at both operational (short-term) and strategic (medium-long-term) level, enabling them to devise adaptation and structural protection measures in a more specific and effective manner.
Moreover, the method is cost-effective and time-efficient on account of ease of data processing, and it can also be a precursor for more detailed qualitative and quantitative analysis of risk (from different hazard types) at different scales for the areas or sectors considered to be at greater risk.
Finally, beyond the site-specific results, the method represents an important contribution toward more comprehensive risk assessment, also in terms of its potential transferability and replicability to different hazard types, including the local effects of climate change from extreme weather/climate events and sea level rise that need to be taken into consideration in a timely and effective manner.

Author Contributions

Conceptualization: A.R., V.V., G.B., A.S.M., M.S.; methodology: A.R., V.V., G.B., A.S.M., M.S.; formal analysis: V.V. and A.R.; investigation: M.S., V.V., A.S.M.; data curation: A.R. and V.V.; writing original draft preparation: A.R. and V.V.; writing review and editing: A.R., V.V., G.B., A.S.M., M.S.; visualization: V.V. and A.R.; supervision: M.S.; project administration: M.S. and A.S.M.; funding acquisition: M.S. and A.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

The study has been carried out in the frame of the Project “Coastal risk assessment and mapping” funded by the EUR-OPA Major Hazards Agreement of the Council of Europe (2020–2021). Grant Number: GA/2020/06 n° 654503. Scientific responsible: Anton Micallef (ICoD) and Mauro Soldati (Unimore).

Acknowledgments

The authors are grateful to Chris Gauci (Research and Planning Section, Marine and Storm Water Unit, Public Works Department, Floriana, Malta) for information provided about the social vulnerability aspects. The authors would also like to thank Chiara Bordin for helping with data processing in GIS environment.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The evaluation of the social vulnerability for the coastal districts located in the investigated area was based on the calculation of a synthetic index that takes into account a number of indicators that represent social protection schemes that provide, where possible, protection against a single risk or need. For a comprehensive description of the social protection gross expenditure, readers can refer to the social protection report for the period 2012–2016, published in 2019 and directly downloadable at: https://nso.gov.mt/en/publicatons/Publications_by_Unit/Documents/A2_Public_Finance/Social%20Protection%202016.pdf
The social indicators taken into account in this study are:
  • The Health Care function, which consists of all those benefits paid to persons during temporary periods of unemployment due to sickness or injury, and health care provided in the framework of social protection;
  • The Disability function, which mainly covers cash benefits paid to persons who are below the retirement age and unable to work because of a mental or physical disability;
  • The Old Age function, which covers all interventions against the risks linked to retirement and ageing. These include pensions given to a person once they retire from the labour market, lodging in specialized retirement homes and any services provided to persons unable to independently care for themselves.
  • The Family/Children function, which includes cash benefits provided to households with children, various childcare services available to families and other social services provided with the specific intention to assist families with children.
  • The Unemployment function, which represents benefits paid to either compensate an individual for the loss of his/her gainful employment or to cover the income of persons who retire from employment prior to the statutory age.

References

  1. Davidson–Arnott, R. An Introduction to Coastal Processes and Geomorphology; Cambridge University Press: Cambridge, UK, 2010; p. 458. [Google Scholar]
  2. Masselink, G.; Gehrels, R. Coastal Environments and Global Change; John Wiley & Sons: West Sussex, UK, 2014; p. 448. [Google Scholar]
  3. Nicholls, R.J.; Wong, P.P.; Burkett, V.R.; Codignotto, J.O.; Hay, J.E.; McLean, R.F.; Ragoonaden, S.; Woodroffe, C.D. Coastal systems and lowlying areas. Climate Change 2007. Impacts, Adaptation and Vulnerability. In Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Parry, M.L., Canziani, O.F., Palutikof, J.P., van der Linden, P.J., Hanson, C.E., Eds.; Cambridge University Press: Cambridge, UK, 2007; pp. 315–356. [Google Scholar]
  4. World Resources Institute. Decision Making in a Changing Climate; United Nations Development Programme World Bank; World Resources Institute: Washington, DC, USA, 2010. [Google Scholar]
  5. European Environment Agency. The Changing Faces of Europe’s Coastal Areas; EEA Report; European Environment Agency: Copenhagen, Denmark, 2006. [Google Scholar]
  6. European Environment Agency. Mediterranean Sea Region Briefing—The European Environment—State and Outlook; EEA Report; European Environment Agency: Copenhagen, Denmark, 2015. [Google Scholar]
  7. Reid, W.V.; Mooney, H.A.; Cropper, A.; Capistrano, D.; Carpenter, S.R.; Chopra, K.; Dasgupta, P.; Dietz, T.; Duraiappah, A.K.; Hassan, R.; et al. Ecosystems and Human Well–Being: Synthesis; Millennium Ecosystem Assessment: Washington, DC, USA, 2005. [Google Scholar]
  8. Maes, J.; Teller, A.; Erhard, M.; Liquete, C.; Braat, L.; Berry, P.; Egoh, B.; Puydarrieux, P.; Fiorina, C.; Santos, F.; et al. Mapping and Assessment of Ecosystems and Their Services. An Analytical Framework for Ecosystem Assessments under Action 5 of the EU Biodiversity Strategy to 2020; Publications Office of the European Union: Luxembourg, Luxembourg, 2013; pp. 1–58. [Google Scholar]
  9. Gallina, V.; Torresan, S.; Critto, A.; Sperotto, A.; Glade, T.; Marcomini, A. A review of multi–risk methodologies for natural hazards: Consequences and challenges for a climate change impact assessment. J. Environ. 2016, 168, 123–132. [Google Scholar] [CrossRef] [PubMed]
  10. Intergovernmental Panel on Climate Change. Climate Change 2007: Mitigation. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Metz, B., Davidson, O.R., Bosch, P.R., Dave, R., Meyer, L.A., Eds.; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  11. Intergovernmental Panel on Climate Change. Climate Change 2014: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Part A: Global Aspects; Field, C.B., Barros, V.R., Dokken, D.J., Mach, K.J., Mastrandrea, M.D., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  12. Hoegh–Guldberg, O.; Jacob, D.; Taylor, M.; Bindi, M.; Brown, S.; Camilloni, I.; Diedhiou, A.; Djalante, R.; Ebi, K.L.; Engelbrecht, F.; et al. 2018: Impacts of 1.5 °C Global Warming on Natural and Human Systems. In Global Warming of 1.5 °C; Masson–Delmotte, V., Zhai, P., Pörtner, H.-O., Roberts, D., Skea, J., Shukla, P.R., Pirani, A., Moufouma–Okia, W., Péan, C., Pidcock, R., et al., Eds.; An IPCC Special Report on the Impacts of Global Warming of 1.5 °C above Pre–Industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty; Intergovernmental Panel on Climate Change: Geneva, Switzerland, 2018. [Google Scholar]
  13. Intergovernmental Panel on Climate Change. Summary for Policymakers. In Climate Change and Land: An IPCC Special Report on Climate Change, Desertification, Land Degradation, Sustainable Land Management, Food Security, and Greenhouse Gas Fluxes in Terrestrial Ecosystems; Shukla, P.R., Skea, J., Buendia, E.C., Masson-Delmotte, V., Pörtner, H.-O., Roberts, D.C., Zhai, P., Slade, R., Connors, S., Diemen, R., et al., Eds.; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  14. UNISDR (United Nations International Strategy for Disaster Reduction). Sendai Framework for Disaster Risk Reduction 2015–2030; The United Nations Office for Disaster Risk Reduction: Geneva, Switzerland, 2015; Available online: https://www.preventionweb.net/files/43291_sendaiframeworkfordrren.pdf (accessed on 26 April 2020).
  15. European Commission. An EU Strategy on Adaptation to Climate Change; The European Commission: Brussels, Belgium, 2013; Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=COM:2013:0216:FIN:EN:PDF (accessed on 26 April 2020).
  16. Directive 2007/60/EC of the European Parliament and of the Council of 23 October 2007 on the Assessment and Management of Flood Risks. Available online: https://eur–lex.europa.eu/legal–content/EN/TXT/PDF/?uri=CELEX:32007L0060&from=EN (accessed on 26 April 2020).
  17. Field, C.B.; Barros, V.; Stocker, T.F.; Qin, D.; Dokken, D.J.; Ebi, K.L.; Mastrandrea, M.D.; Mach, K.J.; Plattner, G.-K.; Allen, S.K.; et al. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation; A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change (IPCC); Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  18. Baird, A.; O’Keefe, P.; Westgate, K.; Wisner, B. Towards an Explanation of and Reduction of Disaster Proneness; Occasional Paper Number 11; Disaster Research Unit, University of Bradford: Bradford, UK, 1975. [Google Scholar]
  19. O’Keefe, P.; Westgate, K.; Wisner, B. Taking the naturalness out of natural disasters. Nature 1976, 260, 566–577. [Google Scholar] [CrossRef] [Green Version]
  20. Lewis, J. The vulnerable state: An alternative view. In Disaster Assistance: Appraisal, Reform and New Approaches; Stephens, L., Green, S.J., Eds.; New York University Press: New York, NY, USA, 1976; pp. 104–129. [Google Scholar]
  21. Hewitt, K. Interpretations of Calamity; Allen & Unwin: London, UK, 1983. [Google Scholar]
  22. O’Brien, K.; Eriksen, S.; Schjolen, A.; Nygaard, L. What’s in a Word? Conflicting Interpretations of Vulnerability in Climate Change Research; CICERO Working Paper 2004:04; CICERO, Oslo University: Oslo, Norway, 2004. [Google Scholar]
  23. Romieu, E.; Welle, T.; Schneiderbauer, S.; Pelling, M.; Vinchon, C. Vulnerability assessment within climate change and natural hazard contexts: Revealing gaps and synergies through coastal applications. Sustain. Sci. 2010, 5, 159–170. [Google Scholar] [CrossRef]
  24. Jurgilevich, A.; Räsänen, A.; Groundstroem, F.; Juhola, S. A systematic review of dynamics in climate risk and vulnerability assessments. Environ. Res. Lett. 2017, 12, 013002. [Google Scholar] [CrossRef]
  25. Mysiak, J.; Torresan, S.; Bosello, F.; Mistry, M.; Amadio, M.; Marzi, S.; Furlan, E.; Sperotto, A. Climate risk index for Italy. Philos. Trans. R. Soc. 2018, 376, 20170305. [Google Scholar] [CrossRef] [PubMed]
  26. UNISDR. UNISDR Terminology on Disaster Risk Reduction. 2009. Available online: https://www.unisdr.org/files/7817_UNISDRTerminologyEnglish.pdf (accessed on 26 April 2020).
  27. Cavallin, A.; Marchetti, M.; Panizza, M.; Soldati, M. The role of geomorphology in the environmental impact assessment. Geomorphology 1994, 9, 143–153. [Google Scholar] [CrossRef]
  28. UNDRO. Disaster Prevention and Mitigation—Compendium of Current Knowledge; United Nations: New York, NY, USA, 1984. [Google Scholar]
  29. Papathoma–Köhle, M.; Gems, B.; Sturm, M.; Fuchs, S. Matrices, curves and indicators: A review of approaches to assess physical vulnerability to debris flows. Earth Sci. Rev. 2017, 171, 272–288. [Google Scholar] [CrossRef]
  30. Intergovernmental Panel on Climate Change. Climate Change 1995. Impacts, Adaptations and Mitigation of Climate Change: Scientific-Technical Analyses; Watson, R.T., Zinyowera, M.C., Moss, H.J.D., Eds.; Second Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University: Cambridge, UK, 1995. [Google Scholar]
  31. Costa, L.; Kropp, J.P. Linking components of vulnerability in theoretic frameworks and case studies. Sustain. Sci. 2013, 8, 1–9. [Google Scholar] [CrossRef]
  32. Gornitz, V. Vulnerability of the East Coast, USA to future sea level rise. J. Coast. Res. 1990, 9, 201–237. [Google Scholar]
  33. Gornitz, V.M.; White, T.W.; Cushman, R.M. Vulnerability of the US to future sea level rise, Coastal Zone 91. In Proceedings of the 7th Symposium on Coastal and Ocean Management, Long Beach, CA, USA, 8–12 July 1991; American Society of Civil Engineers: New York, NY, USA, 1991; pp. 1345–1359. [Google Scholar]
  34. Ojeda–Zújar, J.; Álvarez–Francosi, J.I.; Martín–Cajaraville, D.; Fraile–Jurado, P. El uso de las TIG para el cálculo del índice de Vulnerabilidad costera (CVI) ante una potencial subida del nivel del mar en la costa andaluza (España). GeoFocus 2009, 9, 83–100. [Google Scholar]
  35. Özyurt, G.; Ergin, A. Application of sea level rise vulnerability assessment model to selected coastal areas of Turkey. J. Coast. Res. 2009, 56, 248–251. [Google Scholar]
  36. Özyurt, G.; Ergin, A. Improving coastal vulnerability assessments to sea–level rise: A new indicator–based methodology for decision makers. J. Coast. Res. 2010, 26, 265–273. [Google Scholar] [CrossRef]
  37. McLaughlin, S.; Cooper, J.A.G. A multi-scale coastal vulnerability index: A tool for coastal managers? Environ. Hazard 2010, 9, 233–248. [Google Scholar] [CrossRef]
  38. Di Paola, G.; Iglesias, J.; Rodríguez, G.; Benassai, G.; Aucelli, P.P.C.; Pappone, G. Estimating coastal vulnerability in a meso-tidal beach by means of quantitative and semi–quantitative methodologies. J. Coast. Res. 2011, 61, 303–308. [Google Scholar] [CrossRef]
  39. Santos, M.; Del Río, L.; Benavente, J. GIS–based approach to the assessment of coastal vulnerability to storms. Case study in the Bay of Cádiz (Andalusia, Spain). J. Coast. Res. 2013, 65, 826–831. [Google Scholar] [CrossRef]
  40. Armaroli, C.; Duo, E. Validation of the Coastal storm Risk Assessment Framework along the Emilia–Romagna coast. Coast. Eng. 2018, 134, 159–167. [Google Scholar] [CrossRef] [Green Version]
  41. Van Dongeren, A.; Ciavola, P.; Martinez, G.; Viavattene, C.; Bogaard, T.; Ferreira, O.; McCall, R. Introduction to RISC–KIT: Resilience–increasing strategies for coasts. Coast. Eng. 2018, 134, 2–9. [Google Scholar] [CrossRef] [Green Version]
  42. Viavattene, C.; Jiménez, J.A.; Ferreira, O.; Priest, S.; Owen, D.; McCall, R. Selecting coastal hotspots to storm impacts at the regional scale: A Coastal Risk Assessment Framework. Coast. Eng. 2018, 134, 33–47. [Google Scholar] [CrossRef]
  43. Stelljes, N.; Martinez, G.; McGlade, K. Introduction to the RISC–KIT web based management guide for DRR in European coastal zones. Coast. Eng. 2018, 134, 73–80. [Google Scholar] [CrossRef]
  44. Prampolini, M.; Gauci, C.; Micallef, A.S.; Selmi, L.; Vandelli, V.; Soldati, M. Geomorphology of the north–eastern coast of Gozo (Malta, Mediterranean Sea). J. Maps 2018, 14, 402–410. [Google Scholar] [CrossRef]
  45. Micallef, S.; Micallef, A.; Galdies, C. Application of the Coastal Hazard Wheel to assess erosion on the Maltese coast. Ocean Coast. Manag. 2018, 156, 209–222. [Google Scholar] [CrossRef] [Green Version]
  46. Mantovani, M.; Devoto, S.; Forte, E.; Mocnik, A.; Pasuto, A.; Piacentini, D.; Soldati, M. A multidisciplinary approach for rock spreading and block sliding investigation in the north–western coast of Malta. Landslides 2013, 10, 611–622. [Google Scholar] [CrossRef]
  47. Mantovani, M.; Devoto, S.; Piacentini, D.; Prampolini, M.; Soldati, M.; Pasuto, A. Advanced SAR interferometric analysis to support geomorphological interpretation of slow–moving coastal landslides (Malta Mediterranean Sea). Remote Sens. 2016, 8, 443. [Google Scholar] [CrossRef] [Green Version]
  48. Soldati, M.; Devoto, S.; Foglini, F.; Forte, E.; Mantovani, M.; Pasuto, A.; Piacentini, D.; Prampolini, M. An integrated approach for landslide hazard assessment on the NW coast of Malta. In Proceedings of the International Conference: Georisks in the Mediterranean and Their Mitigation, Valletta, Malta, 20–21 July 2015; Galea, P., Borg, R.P., Farrugia, D., Agius, M.R., D’Amico, S., Torpiano, A., Bonello, M., Eds.; Gutemberg Press Ltd.: Tarxien, Malta, 2015; pp. 160–167. [Google Scholar]
  49. Piacentini, D.; Devoto, S.; Mantovani, M.; Pasuto, A.; Prampolini, M.; Soldati, M. Landslide susceptibility modeling assisted by Persistent Scatterers Interferometry (PSI): An example from the northwestern coast of Malta. Nat. Hazards 2015, 78, 681–697. [Google Scholar] [CrossRef] [Green Version]
  50. Mantovani, M.; Piacentini, D.; Devoto, S.; Prampolini, M.; Pasuto, A.; Soldati, M. Landslide susceptibility analysis exploiting Persistent Scatterers data in the northern coast of Malta. In Proceedings of the International Conference Analysis and Management of Changing Risks for Natural Hazards, Padua, Italy, 18–19 November 2014; pp. 1–7. [Google Scholar]
  51. Viavattene, C.; Jimenez, J.A.; Owen, D.; Priest, S.; Parker, D.; Micou, A.P.; Ly, S. Coastal Risk Assessment Framework Guidance Document. Deliverable No: D.2.3—Coastal Risk Assessment Framework Tool, Risc-Kit Project (G.A. No. 603458). 2015. Available online: http://www.risckit.eu/np4/file/23/RISC_KIT_D2.3_CRAF_Guidance.pdf (accessed on 15 October 2018).
  52. Ferreira, O.; Viavattene, C.; Jiménez, J.; Bole, A.; Plomaritis, T.; Costas, S.; Smets, S. CRAF Phase 1, A framework to identify coastal hotspots to storm impacts. Risk Evaluation & Assessment. In Proceedings of the FLOODrisk 2016—3rd European Conference on Flood Risk Management, Lyon, France, 17–21 October 2016. [Google Scholar] [CrossRef] [Green Version]
  53. Aucelli, P.P.; Di Paola, G.; Rizzo, A.; Rosskopf, C.M. Present day and future scenarios of coastal erosion and flooding processes along the Italian Adriatic coast: The case of Molise region. Envrion. Earth Sci. 2018, 77, 371. [Google Scholar] [CrossRef]
  54. Ballesteros, C.; Jiménez, J.A.; Viavattene, C. A multi–component flood risk assessment in the Maresme coast (NW Mediterranean). Natl. Hazards 2018, 90, 265–292. [Google Scholar] [CrossRef] [Green Version]
  55. Papathoma–Köhle, M.; Cristofari, G.; Wenk, M.; Fuchs, S. The importance of indicator weights for vulnerability indices and implications for decision making in disaster management. Int. J. Disaster Risk. Reduct. 2019, 36, 101103. [Google Scholar] [CrossRef]
  56. Birkmann, J. Indicators and criteria for measuring vulnerability: Theoretical bases and requirements. In Measuring Vulnerability to Natural Hazards: Towards Disaster Resilient Societies; Birkmann, J., Ed.; UNU Press: Tokyo, Japan, 2006. [Google Scholar]
  57. Fuchs, S.; Frazier, T.; Siebeneck, L. Vulnerability and Resilience to Natural Hazards; Fuchs, S., Thaler, T., Eds.; Cambridge University Press: Cambridge, UK, 2018. [Google Scholar]
  58. Nardo, M.; Saisana, M.; Saltelli, A.; Tarantola, S. Tools for Composite Indicators Building; European Commission: Brussels, Belgium, 2005. [Google Scholar]
  59. Soldati, M.; Maquaire, O.; Zezere, J.L.; Piacentini, D.; Lissak, C. Coastline at risk: Methods for multi–hazard assessment. J. Coast. Res. 2011, 61, 335–339. [Google Scholar] [CrossRef]
  60. Foglini, F.; Prampolini, M.; Micallef, A.; Angeletti, L.; Vandelli, V.; Deidun, A.; Soldati, M.; Taviani, M. Late Quaternary Coastal Landscape Morphology and Evolution of the Maltese Islands (Mediterranean Sea) Reconstructed from High–Resolution Seafloor Data. In Geology and Archaeology: Submerged Landscapes of the Continental Shelf; Har, J., Bailey, G., Lüth, L., Eds.; Geological Society, Special Publication: London, UK, 2016; Volume 411, pp. 77–95. [Google Scholar]
  61. Soldati, M.; Barrows, T.T.; Prampolini, M.; Fifield, K.L. Cosmogenic exposure dating constraints for coastal landslide evolution on the Island of Malta (Mediterranean Sea). J. Coast. Conserv. 2018, 22, 831–844. [Google Scholar] [CrossRef] [Green Version]
  62. Oil Exploration Directorate. Geological Map of the Maltese Islands; Office of the Prime Minister: Valletta, Malta, 1993.
  63. Magri, O. A geological and geomorphological review of the Maltese Islands with special reference to the coastal zone. Territoris 2006, 6, 7–26. [Google Scholar]
  64. Micallef, A.; Foglini, F.; Le Bas, T.; Angeletti, L.; Maselli, V.; Pasuto, A.; Taviani, M. The submerged palaeolandscape of the Maltese Islands: Morphology evolution and relation to Quaternary environmental change. Mar. Geol. 2013, 335, 129–147. [Google Scholar] [CrossRef]
  65. Paskoff, R.; Sanlaville, P. Observations geomorphologiques sur le cotes de l’Archipel Maltaise. Z. Geomorphol. 1978, 22, 310–328. [Google Scholar]
  66. Said, G.; Schembri, J. Malta. In Encyclopaedia of the World’s Coastal Landforms; Bird, E., Ed.; Springer: Dordrecht, The Netherlands, 2010; pp. 751–759. [Google Scholar]
  67. Galve, J.P.; Tonelli, C.; Gutiérrez, F.; Lugli, S.; Vescogni, A.; Soldati, M. New insights into the genesis of the Miocene collapse structures of the island of Gozo (Malta, central Mediterranean Sea). J. Geol. Soc. 2015, 172, 336–348. [Google Scholar] [CrossRef] [Green Version]
  68. Mottershead, D.; Bray, M.; Soar, P.; Farres, P.J. Extreme wave events in the central Mediterranean: Geomorphic evidence of tsunami on the Maltese Islands. Z. Geomorphol. 2014, 58, 385–411. [Google Scholar] [CrossRef] [Green Version]
  69. Soldati, M.; Tonelli, C.; Galve, J.P. Geomorphological evolution of palaeosinkhole features in the Maltese archipelago (Mediterranean Sea). Geogr. Fis. Din. Quat. 2013, 36, 189–198. [Google Scholar] [CrossRef]
  70. Coratza, P.; Galve, J.P.; Soldati, M.; Tonelli, C. Recognition and assessment of sinkholes as geosites: Lessons from the Island of Gozo (Malta). Quaest. Geogr. 2012, 31, 22–35. [Google Scholar] [CrossRef] [Green Version]
  71. Coratza, P.; Gauci, R.; Schembri, J.A.; Soldati, M.; Tonelli, C. Bridging Natural and Cultural Values of Sites with Outstanding Scenery: Evidence from Gozo, Maltese Islands. Geoheritage 2016, 8, 91–103. [Google Scholar] [CrossRef]
  72. Cappadonia, C.; Coratza, P.; Agnesi, V.; Soldati, M. Malta and Sicily Joined by Geoheritage Enhancement and Geotourism within the Framework of Land Management and Development. Geosci. J. 2018, 8, 253. [Google Scholar] [CrossRef] [Green Version]
  73. Satariano, B.; Gauci, R. Landform loss and its effect on health and well–being: The collapse of the Azure Window (Gozo) and the resultant reactions of the media and the Maltese community. In Landscapes and Landforms of the Maltese Islands; Gauci, R., Schembri, J.A., Eds.; Springer: Cham, Switzerland, 2019; Volume 14, pp. 289–303. [Google Scholar]
  74. Baldassini, N.; Di Stefano, A. Stratigraphic features of the Maltese Archipelago: A synthesis. Nat. Hazards 2017, 86, 203–231. [Google Scholar] [CrossRef]
  75. Gauci, R.; Schembri, J.A. (Eds.) Landscapes and Landforms of the Maltese Islands; Springer: Cham, Switzerland, 2019. [Google Scholar]
  76. Pedley, M.; Clarke, M.H. Limestone Isles in a Crystal Sea: The Geology of the Maltese Islands; Publishers Enterprises Group: San Gwann, Malta, 2002. [Google Scholar]
  77. Biolchi, S.; Furlani, S.; Devoto, S.; Gauci, R.; Castaldini, D.; Soldati, M. Geomorphological identification, classification and spatial distribution of coastal landforms of Malta (Mediterranean Sea). J. Maps 2016, 12, 87–99. [Google Scholar] [CrossRef] [Green Version]
  78. Furlani, S.; Antonioli, F.; Gambin, T.; Gauci, R.; Ninfo, A.; Zavagno, E.; Micallef, A.; Cucchi, F. Marine notches in the Maltese islands (central Mediterranean Sea). Quat. Int. 2017, 439, 158–168. [Google Scholar] [CrossRef] [Green Version]
  79. Alexander, D. A review of the physical geography of Malta and its significance for tectonic geomorphology. Quat. Sci. Rev. 1988, 7, 41–53. [Google Scholar] [CrossRef]
  80. Soldati, M.; Devoto, S.; Prampolini, M.; Pasuto, A. The spectacular landslide–controlled landscape of the northwestern coast of Malta. In Landscapes and Landforms of the Maltese Islands; Gauci, R., Schembri, J.A., Eds.; Springer: Cham, Switzerland, 2019; Volume 14, pp. 167–178. [Google Scholar]
  81. Devoto, S.; Biolchi, S.; Bruschi, V.M.; Furlani, S.; Mantovani, M.; Piacentini, D.; Soldati, M. Geomorphological map of the NW coast of the Island of Malta (Mediterranean Sea). J. Maps 2012, 8, 33–40. [Google Scholar] [CrossRef]
  82. Devoto, S.; Biolchi, S.; Bruschi, V.M.; González, D.A.; Mantovani, M.; Pasuto, A.; Soldati, M. Landslides Along the North–West Coast of the Island of Malta. In Landslide Science and Practice: Landslide Inventory and Susceptibility and Hazard Zoning; Margottini, C., Canuti, P., Sassa, K., Eds.; Springer: Berlin, Germany, 2013; Volume 1, pp. 57–63. [Google Scholar]
  83. Devoto, S.; Forte, E.; Mantovani, M.; Mocnik, A.; Pasuto, A.; Piacentini, D.; Soldati, M. Integrated Monitoring of Lateral Spreading Phenomena Along the North–West Coast of the Island of Malta. In Landslide Science and Practice: Early Warning, Instrumentation and Monitoring; Margottini, C., Canuti, P., Sassa, K., Eds.; Springer: Berlin, Germany, 2013; Volume 2, pp. 235–241. [Google Scholar]
  84. Magri, O.; Mantovani, M.; Pasuto, A.; Soldati, M. Geomorphological investigation and monitoring of lateral spreading along the north–west coast of Malta. Geogr. Fis. Din. Quat. 2008, 31, 171–180. [Google Scholar]
  85. Pasuto, A.; Soldati, M. Lateral Spreading. In Treatise on Geomorphology; Shroder, J.F., Ed.; Academic Press: San Diego, CA, USA, 2013; Volume 7, pp. 239–248. [Google Scholar]
  86. Dykes, A.P. Mass movements and conservation management in Malta. J. Environ. Manag. 2002, 66, 77–89. [Google Scholar] [CrossRef]
  87. Prampolini, M.; Foglini, F.; Biolchi, S.; Devoto, S.; Angelini, S.; Soldati, M. Geomorphological mapping of terrestrial and marine areas, northern Malta and comino (Central Mediterranean sea). J. Maps 2017, 13, 457–469. [Google Scholar] [CrossRef]
  88. Selmi, L.; Coratza, P.; Gauci, R.; Soldati, M. Geoheritage as a Tool for Environmental Management: A Case Study in Northern Malta (Central Mediterranean Sea). Resources 2019, 8, 168. [Google Scholar] [CrossRef] [Green Version]
  89. Central Bank of Malta. The Evolution of Malta’s Tourism Product over Recent Years. Available online: https://www.centralbankmalta.org/file.aspx?f=72256 (accessed on 17 October 2019).
  90. World Travel and Tourism Council. Travel & Tourism Economic Impact 2018 Malta; World Travel and Tourism Council: London, UK, 2018. [Google Scholar]
  91. Chaperon, S.; Bramwell, B. Dependency and agency in peripheral tourism development. Ann. Tour. Res. 2013, 40, 132–154. [Google Scholar] [CrossRef]
  92. Ebejer, J.; Mangion, M.J.; Bingül, M.B.; Kwiatkowska, D. Rural Landscape and Tourism: A Proposed Policy for Sustainable Tourism in Gozo. In Proceedings of the 7th LE: NOTRE Landscape Forum 2018, Gozo, Malta, 20–24 March 2018; Available online: https://www.um.edu.mt/library/oar/bitstream/123456789/40147/1/Article%20on%20Gozo%20tourism%20for%20LeNOTRE%20Landscape%20FINAL%20DRAFT%20for%20OAR.pdf (accessed on 22 November 2019).
  93. Malta Tourism Authority. Tourism in Malta—Facts and Figures. 2017. Available online: https://www.mta.com.mt/en/file.aspx?f=32328 (accessed on 16 October 2019).
  94. National Statistics Office. Regional Statistics—Malta, 2019th ed.; National Statistics Office: Valletta, Malta, 2019.
  95. Salman, A.; Lombardo, S.; Doody, P. Living with Coastal Erosion in Europe: Sediment and Space for Sustainability; Technical Report; EUCC: Warnemünde, Germany, 2004. [Google Scholar]
  96. Geofabrik. Available online: http://download.geofabrik.de/europe/malta.htm (accessed on 3 July 2019).
  97. Malta Inspire Geoportal. Available online: https://msdi.data.gov.mt/geoportal.html (accessed on 3 July 2019).
  98. National Statistics Office. Social Protection—Reference Years 2012–2016. Malta, 2019. Available online: https://nso.gov.mt/en/publicatons/Publications_by_Unit/Documents/A2_Public_Finance/Social%20Protection%202016.pdf (accessed on 3 July 2019).
  99. Beccari, B. A comparative analysis of disaster risk, vulnerability and resilience composite indicators. PLoS Curr. 2016, 8. [Google Scholar] [CrossRef]
  100. Del Río, L.; Gracia, F.J. Erosion risk assessment of active coastal cliffs in temperate environments. Geomorphology 2009, 112, 82–95. [Google Scholar] [CrossRef] [Green Version]
  101. Di Paola, G.; Aucelli, P.P.C.; Benassai, G.; Iglesias, J.; Rodríguez, G.; Rosskopf, C.M. The assessment of the coastal vulnerability and exposure degree of Gran Canaria Island (Spain) with a focus on the coastal risk of Las Canteras Beach in Las Palmas de Gran Canaria. J. Coast. Conserv. 2018, 22, 1001–1014. [Google Scholar] [CrossRef]
  102. Mattei, G.; Rizzo, A.; Anfuso, G.; Aucelli, P.P.C.; Gracia, F.J. A tool for evaluating the archaeological heritage vulnerability to coastal processes: The case study of Naples Gulf (southern Italy). Ocean Coast. Manag. 2019, 179, 104876. [Google Scholar] [CrossRef]
  103. Rizzo, A.; Aucelli, P.P.C.; Gracia, F.J.; Anfuso, G. A novelty coastal susceptibility assessment method: Application to Valdelagrana area (SW Spain). J. Coast. Conserv. 2018, 22, 973–987. [Google Scholar] [CrossRef]
  104. Giorgi, F. Climate change hot-spots. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  105. Giorgi, F.; Lionello, P. Climate change projections for the Mediterranean region. Glob. Planet. Chang. 2008, 63, 90–104. [Google Scholar] [CrossRef]
  106. Hov, Ø.; Cubasch, U.; Fischer, E.; Höppe, P.; Iversen, T.; Gunnar Kvamstø, N.; Zbigniew, W.K.; Rezacova, D.; Rios, D.; Santos, F.D.; et al. Extreme Weather Events in Europe: Preparing for Climate Change Adaptation. Norwegian Meteorological Institute: Oslo, Norway, 2013.
  107. Nicholls, R.J. Impacts of and Responses to Sea-Level Rise. Understanding Sea-Level Rise and Variability; Church, J.A., Woodworth, P.L., Aarup, T., Wilson, W.W., Eds.; Wiley–Blackwell: Hoboken, NJ, USA, 2010; pp. 17–51. [Google Scholar]
  108. Lambeck, K.; Antonioli, F.; Anzidei, M.; Ferranti, L.; Leoni, G.; Scicchitano, G.; Silenzi, S. Sea level change along the Italian coast during the Holocene and projections for the future. Quat. Int. 2011, 232, 250–257. [Google Scholar] [CrossRef]
  109. Aucelli, P.P.C.; Cinque, A.; Mattei, G.; Pappone, G. Historical sea level changes and effects on the coasts of Sorrento Peninsula (Gulf of Naples): New constrains from recent geoarchaeological investigations. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2016, 463, 112–125. [Google Scholar] [CrossRef]
  110. Tsimplis, M.N.; Calafat, F.M.; Marcos, M.; Jordá, G.; Gomis, D.; Fenoglio-Marc, L.; Struglia, S.; Josey, S.; Chambers, D.P. The effect of the NAO on sea level and on mass changes in the Mediterranean Sea. J. Geophys. Res. Oceans 2013, 118, 944–952. [Google Scholar] [CrossRef] [Green Version]
  111. Anzidei, M.; Lambeck, K.; Antonioli, F.; Furlani, S.; Mastronuzzi, G.; Serpelloni, E.; Vannucci, G. Coastal structure, sea–level changes and vertical motion of the land in the Mediterranean. Geol. Soc. 2014, 388, 453–479. [Google Scholar] [CrossRef]
  112. Galassi, G.; Spada, G. Sea-level rise in the Mediterranean Sea by 2050: Roles of terrestrial ice melt, steric effects and glacial isostatic adjustment. Glob. Planet. Chang. 2014, 123, 55–66. [Google Scholar] [CrossRef]
  113. Antonioli, F.; Anzidei, M.; Amorosi, A.; Presti, V.L.; Mastronuzzi, G.; Deiana, G.; Marsico, A. Sea–level rise and potential drowning of the Italian coastal plains: Flooding risk scenarios for 2100. Quat. Sci. Rev. 2017, 158, 29–43. [Google Scholar] [CrossRef] [Green Version]
  114. Antonioli, F.; Defalco, G.; Moretti, L.; Anzidei, M.; Bonaldo, D.; Carniel, S.; Leoni, G.; Furlani, S.; Presti, V.I.; Mastronuzzi, G.; et al. Relative sea level rise and potential flooding risk for 2100 on 15 coastal plains of the Mediterranean Sea. Geophys. Res. Abstr. 2019, 21, 5274. [Google Scholar]
  115. Aucelli, P.P.C.; Di Paola, G.; Incontri, P.; Rizzo, A.; Vilardo, G.; Benassai, G.; Buonocuore, B.; Pappone, G. Coastal inundation risk assessment due to subsidence and sea level rise in a Mediterranean alluvial plain (Volturno coastal plain–southern Italy). Estuar. Coast. Shelf Sci. 2017, 198, 597–609. [Google Scholar] [CrossRef]
  116. Di Paola, G.; Alberico, I.; Aucelli, P.P.C.; Matano, F.; Rizzo, A.; Vilardo, G. Coastal subsidence detected by Synthetic Aperture Radar interferometry and its effects coupled with future sea-level rise: The case of the Sele Plain (Southern Italy). J. Flood Risk Manag. 2018, 11, 191–206. [Google Scholar] [CrossRef] [Green Version]
  117. Nicholls, R.J.; Cazenave, A. Sea–level rise and its impact on coastal zones. Science 2010, 328, 1517–1520. [Google Scholar] [CrossRef] [PubMed]
  118. Revell, D.L.; Battalio, R.; Spear, B.; Ruggiero, P.; Vandever, J. A methodology for predicting future coastal hazards due to SLR on the California Coast. Clim. Chang. 2011, 109, 251–276. [Google Scholar] [CrossRef]
  119. Kulp, S.A.; Strauss, B.H. New elevation data triple estimates of global vulnerability to sea–level rise and coastal flooding. Nat. Commun. 2019, 10, 1–12. [Google Scholar]
Figure 1. Geological setting of the Maltese archipelago (cf. [62]) and location of the study area, bounded by the NE coastline of Gozo and the red dashed line parallel to it on the inland side (after [44], modified).
Figure 1. Geological setting of the Maltese archipelago (cf. [62]) and location of the study area, bounded by the NE coastline of Gozo and the red dashed line parallel to it on the inland side (after [44], modified).
Water 12 01405 g001
Figure 2. Main geomorphic features of the study area: (a) block slides (west of Daħlet Qorrot Bay); (b) rock fall at the bottom of a limestone plateau and earth flow/slide affecting the underlying clayey terrain (between San Blas Bay and Daħlet Qorrot Bay); (c) plunging cliff between Daħlet Qorrot Bay and Ras il-Qala; (d) sloping coast between Daħlet Qorrot Bay and Ras il-Qala; (e) shore platform east of Marsalforn Bay; (f) Ramla Bay pocket beach; (g) cliff shaped in Blue Clay east of Marsalforn Bay; (h) built-up coast of Marsalforn Bay. After [44], modified.
Figure 2. Main geomorphic features of the study area: (a) block slides (west of Daħlet Qorrot Bay); (b) rock fall at the bottom of a limestone plateau and earth flow/slide affecting the underlying clayey terrain (between San Blas Bay and Daħlet Qorrot Bay); (c) plunging cliff between Daħlet Qorrot Bay and Ras il-Qala; (d) sloping coast between Daħlet Qorrot Bay and Ras il-Qala; (e) shore platform east of Marsalforn Bay; (f) Ramla Bay pocket beach; (g) cliff shaped in Blue Clay east of Marsalforn Bay; (h) built-up coast of Marsalforn Bay. After [44], modified.
Water 12 01405 g002
Figure 3. Distribution of coastal morphotypes in the study area (after [44], modified).
Figure 3. Distribution of coastal morphotypes in the study area (after [44], modified).
Water 12 01405 g003
Figure 4. Detailed land use map. The red dashed line identifies the inland boundary of the study area. The districts’ boundaries are also indicated.
Figure 4. Detailed land use map. The red dashed line identifies the inland boundary of the study area. The districts’ boundaries are also indicated.
Water 12 01405 g004
Figure 5. Percentage of the area occupied by each land use category according to the classification shown in Figure 4.
Figure 5. Percentage of the area occupied by each land use category according to the classification shown in Figure 4.
Water 12 01405 g005
Figure 6. Physical vulnerability map resulting from data overlay of the physical indicators (land use, transport and utilities). The red dashed line identifies the inland boundary of the study area.
Figure 6. Physical vulnerability map resulting from data overlay of the physical indicators (land use, transport and utilities). The red dashed line identifies the inland boundary of the study area.
Water 12 01405 g006
Figure 7. Social vulnerability map of the four administrative districts located in NE Gozo. The red dashed line indicates the inland boundary of the study area.
Figure 7. Social vulnerability map of the four administrative districts located in NE Gozo. The red dashed line indicates the inland boundary of the study area.
Water 12 01405 g007
Figure 8. Overall vulnerability map resulting from the spatial aggregation of the physical vulnerability levels and the social vulnerability levels over the area. The red dashed line indicates the inland boundary of the study area.
Figure 8. Overall vulnerability map resulting from the spatial aggregation of the physical vulnerability levels and the social vulnerability levels over the area. The red dashed line indicates the inland boundary of the study area.
Water 12 01405 g008
Table 1. Physical indicators and expert-based exposure level (ranging from value 1 to value 5) assigned to each category of assets.
Table 1. Physical indicators and expert-based exposure level (ranging from value 1 to value 5) assigned to each category of assets.
Physical Indicators
(Value)
Very Low/Null Exposure
(1)
Low Exposure
(2)
Medium Exposure
(3)
High Exposure
(4)
Very High Exposure
(5)
Land UseAbandoned agricultural area.
Bare rock
Agricultural area;
Green urban areas;
Natural and vegetated area,
Terraced agricultural field; Land principally occupied by agriculture with significant areas of natural vegetation
Residential area; Dump; Quarry; CemeteryBeaches; Dunes; SandHistorical and archaeological site; Natural protected area (SCI); Strategic elements; Entertainment (commerce, finance, business, recreational, leisure, and sport)
TransportAbsence of transport network or highly degraded transport networkFootway; Path; Track; StepsTertiary road; Living street; Residential; ServicesSecondary roadPrimary road
UtilitiesAbsence of utilitiesMainly local and small utilitiesStreet lighting11kV overhead lineSubstations; Feeder pillar
Table 2. Areal distribution of each physical vulnerability level resulting from the overlay and aggregation of the physical indicators (land use, transport network and utilities).
Table 2. Areal distribution of each physical vulnerability level resulting from the overlay and aggregation of the physical indicators (land use, transport network and utilities).
Physical VulnerabilitySurface (km2)Surface (%)
Very low2.021.7
Low5.457.8
Medium0.77.4
High0.0050.1
Very high1.213.0
Table 3. Social indicators evaluated for the four coastal districts of the study area. The method for attributing the numerical value assigned to each indicator is indicated in the methodological paragraph (cf. Section 4.3).
Table 3. Social indicators evaluated for the four coastal districts of the study area. The method for attributing the numerical value assigned to each indicator is indicated in the methodological paragraph (cf. Section 4.3).
Social IndicatorsŻebbuġ DistrictXagħra DistrictNadur DistrictQala District
Health care indicator3233
Disability indicator3545
Old age indicator2333
Family/Children indicator3333
Unemployment indicator5222
Population (number of inhabitants, 2016)2043402940011885
Table 4. Areal distribution of overall vulnerability levels that are the result of the combination and overlay of spatial distributions of the physical and social vulnerability factors.
Table 4. Areal distribution of overall vulnerability levels that are the result of the combination and overlay of spatial distributions of the physical and social vulnerability factors.
Overall VulnerabilitySurface (km2)Surface (%)
Very low--
Low1.920.4
Medium5.661.3
High0.77.3
Very high111

Share and Cite

MDPI and ACS Style

Rizzo, A.; Vandelli, V.; Buhagiar, G.; Micallef, A.S.; Soldati, M. Coastal Vulnerability Assessment along the North-Eastern Sector of Gozo Island (Malta, Mediterranean Sea). Water 2020, 12, 1405. https://doi.org/10.3390/w12051405

AMA Style

Rizzo A, Vandelli V, Buhagiar G, Micallef AS, Soldati M. Coastal Vulnerability Assessment along the North-Eastern Sector of Gozo Island (Malta, Mediterranean Sea). Water. 2020; 12(5):1405. https://doi.org/10.3390/w12051405

Chicago/Turabian Style

Rizzo, Angela, Vittoria Vandelli, George Buhagiar, Anton S. Micallef, and Mauro Soldati. 2020. "Coastal Vulnerability Assessment along the North-Eastern Sector of Gozo Island (Malta, Mediterranean Sea)" Water 12, no. 5: 1405. https://doi.org/10.3390/w12051405

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop