Damaging Hydrogeological Events: a Procedure for the Assessment of Severity Levels and an Application to Calabria (southern Italy)

A damaging hydrogeological event (DHE) is characterized by two components: a rainfall event and a subsequent damage event, which is the result of floods and landslides triggered by rainfall. The characteristics of both events depend on climatic, geomorphological and anthropogenic factors. In this paper, a methodology to classify the severity of DHEs is presented. A chart which considers indicators of both the damage (Dscore) and the daily rainfall (Rscore) values recorded in the study area is proposed. According to the chart, the events are classified into four types: ordinary events, with low Dscore and Rscore values; extraordinary events, with high Rscore values but low Dscore values; catastrophic events, characterized by non-exceptional rainfall (low Rscore) and severe damage (high Dscore); major catastrophic events, obtained by both high Dscore and Rscore values. Using this approach, the 2013 DHE that occurred in Calabria (Italy) was classified as an ordinary event, when compared to the previous ones, even though the widespread diffusion of damage data induced the perception of high severity damage. The rainfall that triggered this event confirms the negative trend of heavy daily precipitation detected in Calabria, and the damage can be ascribed more to sub-daily than daily rainfall affecting urbanized flood-prone areas.


Introduction
Damaging hydrogeological events (DHEs) are defined as the occurrence of phenomena triggered by rainfall, such as landslides and floods, which cause damage to people and property [1][2][3].These events can have catastrophic effects [4] or produce more or less severe damage, according to both the socioeconomic context of the affected areas [5,6] and human behaviors.
In some countries, such as Italy, there has been an increase in both the recurrence and the severity of damage caused by rainfall over the last few years.Conventionally, during the past century, DHEs were identified following the year of occurrence (e.g., the flood of 1966); nowadays, the event is identified by means of both the year and the month of occurrence, given that a region can be affected by more than one event per year.Moreover, due to the influence of the mass media, people are inclined to think that current events are triggered by increasingly exceptional rainfall.For example, in Italy, the media have introduced the term bombe d'acqua (literally: water bombs), to label heavy rainfall, and they often tend to amplify the qualitative description of the downpours, even though neither quantitative rainfall data nor a comparison with past events are considered.Nevertheless, it is biased to describe an event as exceptional without any comparison with past ones.One may be reminded of what Rudyard Kipling wrote in [7]: In August was the Jackal born; The Rain fell in September; "Now such a fearful flood as this," Says he, "I can't remember!"Despite the misinformation and the average low awareness on these topics, scientific literature has widely analyzed both rainfall and damage characterizing DHEs.
For example, several studies on rainfall variability have been undertaken using various approaches and obtaining different results.Research carried out in the USA [8,9], Australia [10], New Zealand [11], South Africa [12], the United Kingdom [13], and the Mediterranean Basin [14] has highlighted an increase in both daily rainfall intensity and frequency of extreme rainfall events.Some authors have detected a decrease in total annual rainfall [15][16][17][18][19], although an increase in extreme rainfall has also been evidenced [14].By contrast, other authors [20][21][22] have detected no significant trends in either the number of extreme rainfall events or in rainfall values.In Calabria (southern Italy), some authors (e.g., [23]), considering an observation period of 100 years, have shown a negative trend in the heavy precipitation events, and an opposite trend in the low precipitation events.Generally, daily precipitation is studied using approaches based on indices, percentiles, thresholds, and extreme values [21,[24][25][26].Two approaches are available when focusing on changes in extreme rainfall.The first considers changes in the daily rainfall distribution tail, namely the magnitude of a fixed percentile or the percentage of seasonal precipitation falling on days with rainfall greater than a fixed percentile [27].The second one is based on the evaluation of changes in the number of days over a threshold associated with different rainfall categories [28,29].Several authors (e.g., [30]) have proposed an extreme value approach based on the Gumbel probabilistic model [31], which can be useful to assess the values of the daily rainfall return period.
In order to determine whether damage is currently more frequent/widespread than in the past, one should take into account the effects that intense urbanization, especially in developed countries, has on risk levels, as more and more elements become vulnerable to floods and landslides.There is evidence today that social change and economic development are the main factors responsible for increasing losses from natural disasters [32][33][34].Risk modifications are quite common on floodplains, especially where people built and/or raise levees to reduce the frequency of flooding.In fact, due to the reduction of flooding frequency, people feel safer, and more intense economic development occurs close to rivers.Human settlements close to rivers become vulnerable to high-consequence and low-probability events [35,36].Similarly, even landslide risk has increased as a result of urban expansion.In fact, especially in regions with a rugged morphology, the demand of new spaces for urbanization can lead to building houses and infrastructures along unstable slopes, where landslides can occur during/after heavy rainfall [37].Fortunately, social development increases disaster prevention, which partly reduces the number of casualties and losses: both the total death toll and the number of casualties per event tend to decrease [38].Nevertheless, the increasing wealth causes relatively greater losses in high-income nations [39].Higher incomes increase private demand for safety; higher incomes enable individuals (and countries) to respond to risk by employing additional and costly precautionary measures [40].However, in those countries where concentration of properties is larger than the installation of disaster countermeasures, the relationship between income and vulnerability could even be inverted, especially for disasters whose exposure is related to individual choices, such as floods and landslides.In countries with a GDP per capita below $4,500-$5,500, the death toll resulting from natural disasters appears to increase with income, but it starts to decrease once the GDP of these countries increases above that GDP threshold [41].Moreover, the incessant flow of virtual real-time information and the media coverage of an event tend to bias our perception of the intensity of that event and to overestimate the amount of damage caused by recent rather than older events.For example, whereas in the past century it would take 24 h for newspapers-the most frequently used damage data sources [42][43][44][45]-to report an event, technological progress has allowed an extremely rapid relay of information from the site of a disaster.Furthermore, in the past, only major events received media coverage and were therefore reported.Data concerning less severe events, instead, were primarily available for the most important towns alone, where local reporters were headquartered.By contrast, today, thanks to the internet and social media, every individual can report an event by posting personal photos, videos, and comments-regardless of the intensity of the event and damage caused.
This paper presents a methodology to assess the severity levels of DHEs, based on both rainfall and damage intensities of the past events that occurred in the study area.Particularly, an application of the methodology to a region of Calabria, Italy is presented, and the assessment of the severity of the DHE that occurred on 19 November 2013 is performed.

Materials and Methods
Each DHE is characterized by two components: a rainfall event (RE) and a damage event (DE) resulting from rainfall-triggered phenomena.To assess the relative exceptionality of a DHE, these two components must be investigated, and some indicators of the events' severity must be introduced in order to compare current events to those already occurred in the area during the past.Currently, a chart obtained by comparing damage/rainfall indicators has been created, which allows easy assessment of the severity of current events by comparing them to those that formerly affected the study area.The main limitation of this approach is that the chart is representative of a specific area, based on its climatic, geological, geomorphological, and anthropogenic characteristics.Nevertheless, this methodology can be applied to every area for which the kind of data described in the following sections are available.

Damage Data Gathering and Analysis
Information about the damage caused by DHEs is not directly available and it must be specifically collected, as described in previous papers [1,5].Data sources change according to both the country and the epoch in which the events occurred.For older cases, typical sources are historical archives, newspapers and technical reports.When dealing with printed documents, data gathering requires a careful analysis of huge amounts of documents, and the long-lasting transcription of the outcomes.On the contrary, for the current cases, data are immediately available on the internet and this significantly reduces data acquisition time.Online newspapers and websites dealing with meteorology are among the easiest to consult.Technical agencies publish reports on their websites, and social networks publicize a large amount of pictures and comments about damage events.
Generally, data concern the direct damage, which includes physical effects such as the destruction and the modifications that reduce the functionality of structures, in addition to damage to people (death/injury), buildings and their contents, and vehicles, and all subsequent clean-up and disposal costs.On the contrary, indirect and intangible damages are usually unavailable.Specifically, for past events, the economic figures of direct damage are unavailable or incomplete.By contrast, for current events, media publish loss estimates directly after the event, even if these figures are unreliable due to: (a) assessments being overstated in the hope of mobilizing further emergency aid; (b) the initial evaluation no longer being valid because the loss figures develop as the event evolves [46].
In the proposed methodology, following [47], data mined from the different sources are transcribed in a damage database, in which each record contains the text describing the damage affecting a municipality on a certain date, and this information is then split in a series of fields (event identification, time, damaging phenomena, and damaged elements).These data are then elaborated, aiming at producing four outputs: the first three (a, b, and c) are damage indicators, while the last (d) leads to the identification of the rainfall event that triggered damage: a) The index of damaged area (IDA) is an assessment of the percentage of regional area damaged.It is the sum of the area of damaged municipalities divided by the area of the region.It is a proxy of the extent of the effects, because exact delimitations of affected areas are not available; b) The damage index (DI) is a relative assessment of direct damage.The damageable elements are sorted into nine types (roads; railways; houses; public buildings; services; productive activities; tourist and sport resorts; hydraulic works; people), and their values ranges from 1 to 10, on an arbitrary scale, while the levels of loss is set as 1 = high, 0.5 = medium and 0.25 = low [47].
The sum of the products of damaged elements by the respective levels of loss is the DI value.
For people, the value has been set as 100, and the levels of loss have been set as L1 = more than 10 victims; L2 = 5-10 victims; L3 = less than 5 victims; c) The number of victims (NoV) is used as a further indicator of damage severity; d) The duration of damage event (dDE) is the number of days between the first and the last day during which, throughout the region, some damage related to rainfall occurred.
The values of IDA and DI have been grouped into four classes, individuated by dividing in four equal intervals the difference between the maximum and the minimum observed values.Concerning the number of victims, the four classes have been set as follows: (1) 0-10 victims; (2) 11-20 victims; (3) 21-50; (4) >50 victims.

Rainfall Data Gathering and Analysis
To outline the duration of the triggering RE, an appropriate spatial/temporal resolution is required.From a spatial point of view, with the aim of having a complete description of the event, all the regional rain gauges that were working at the time of each event must be analyzed.Concerning the temporal resolution, rainy days (days with precipitation greater than or equal to 1 mm) occurred during the damage event and in a buffer before and after it (i.e., ten days before and ten after) must be investigated.Thus, for each event, a matrix having a number of columns corresponding to the number of gauges, and a number of rows made by the damage event duration plus the buffering can be obtained.For each gauge, the local rainfall event is defined as the continuous sequence of rainy days either included in, or intersecting, the DE.Then, at a regional scale, the duration of rainfall event (DRE) is the envelope of all the local rainfall events recorded by the available gauges.Actually, it is the number of days between the early beginning and the latest end of rainy days among all the affected gauges [48].As a measure of the magnitude of rainfall events, in the proposed methodology, the return period (T) of the daily rainfall recorded during the event has been evaluated.To assess T using the Gumbel probabilistic model [31], the series of annual maxima of daily rainfall recorded during the hydrological years for a period long enough to allow statistical analyses (no less than 30 years), is required.The Gumbel distribution, also known as the Extreme Value Type I distribution, has the following cumulative distribution function: where z = (x − μ)/σ; μ is the location parameter; and σ is the scale parameter (σ > 0).In extreme-value studies, the probability of exceedance of a certain value is usually expressed in terms of the return period T that is defined as follows:

Rainfall and Damage Analysis
To summarize the severity of a DHE, two indices are proposed: the damage score (Dscore) and the rainfall score (Rscore).
For each DHE, Dscore is the sum of damage indicators obtained from the damage data analysis: score where IDA is the index of damaged area; DI is the damage index and NoV is the number of victims.
Rscore is based on the return period of the daily rainfall, sorted in a series of classes.It considers the percentage of gauges falling within each return period class, using a weight factor that accounts for the importance of the higher values of the return period.
( ) ( ) In Equation ( 4), n is the number classes of return periods; i is a return period class; P(T)i is the percentage of rainfall gauges falling within each i return period class, and c0 is a graphic scale factor.
The different return period classes have been chosen according to the institution responsible for monitoring the hydrogeological risk in the Calabria region, the Basin Authority, which has used several return period classes in its Hydrogeological Structure Plan (Piano di Assetto Idrogeologico or PAI), as required by Italian law.
By plotting these two indices in a chart, four sectors corresponding to four severity levels can be recognized.The name of the sectors are similar to those introduced in [49].Sector A groups less severe cases, which are identified as ordinary events, having low Dscore and Rscore values.Sector B includes extraordinary events, characterized by very high rainfall severity and moderately low damage levels.Sector C includes catastrophic events, showing severe damage caused by not-so-extraordinary rainfall.Finally, Sector D includes major catastrophic events, showing high damage and triggering rainfall values.The realization and the use of this chart by using the data of an Italian case study is explained in the following sections.

The Study Area: Calabria Region
Calabria (southern Italy) has an area of 15,080 km 2 ; its mean altitude is 597 m and the maximum is 2266 m.Although Calabria does not have many high summits, it is one of the most mountainous regions in Italy (Figure 1): 42% of the land is mountainous (elevation greater than 500 m a.s.l.), 49% hilly (elevation between 50 and 500 m a.s.l.) and only 9% is flat (elevation lower than 50 m a.s.l.).
The climate is Mediterranean: during the summer, precipitation is less frequent, while in the winter, rainfall and snow occur.Warm currents coming from Africa affect the east side, causing high temperatures with short and heavy precipitation.By contrast, western air currents affect the west side, producing milder temperatures and frequent rainfall.The inland zones show colder winters, with snow and fresher summers [50,51].Numerous DHEs have caused damages and created victims for several centuries in the region [52][53][54][55][56][57], especially along the east side.

The 2013 Damaging Hydrogeological Event
On 18 November 2013, a low-pressure area on the Balearic Islands channeled hot-humid air from the southwest that merged with both colder high-flying air and water vapor raising from the sea, thus originating the Cleopatra cyclone.Sirocco gusts (Sirocco, also called Leveche in other parts of the Mediterranean area, is a Mediterranean wind that comes from the Sahara and reaches hurricane speeds in North Africa and Southern Europe.) with speed up to 95 km/h hit Sardinia (Italy), where daily rainfall exceeded 400 mm (60% of the annual regional average).As a result, flash floods and landslides damaged 61 municipalities and killed 18 people.
Sirocco gusts also affected the east side of Calabria, with peak velocities ranging between 20 and 100 km/h and averages between 14 and 26 km/h, causing sea waves that reached 9 m in height.Then, between the night of 18th and the morning of 19th, Cleopatra moved towards Calabria.Daily rainfall was not extraordinary: only a gauge on the northeast sector recorded a return period greater than 200 years (Figure 2).Rainfall started in the early morning and lasted until 11:00 a.m.: 1) From midnight to 3:00 a.m.: the rainfall occurred over the middle and southeast side, with a maximum intensity of 26 mm/1 h (Figure 3A).2) From 3:00 to 6:00 a.m.: rainfall affected a larger area, involving the western sector, and its intensity increased up to 70 mm/1 h (Figure 3B). 3) From 6:00 to 9:00 a.m.: the event extended over the north and mid-east sector, with intensities around 50 mm/1 h, and greater than 150 mm/3 h (Figure 3C).4) From 9:00 to 11:00 a.m.: precipitation progressively decreased all over the region (Figure 3D).
The simultaneous occurrence of rainfall and sea storms had a cascading effect: waves obstructed the outlet of floods toward the sea, pushing water upstream.Floods suddenly affected east-side rivers: in 51% of the reported cases, waters overstepped both natural and artificial embankments.Numerous cars were submerged, and some people abandoned their houses just before the latter were flooded.Despite the flood's level exceeding one meter (Figure 4), there were no casualties; two people were slightly injured and some others suffered panic attacks due to flooding.According to the available data, no dangerous behaviors were reported.It is likely that the television broadcasting of the dramatic event that affected Sardinia the day before induced people to behave more cautiously than usual [5].Newspapers (Il Quotidiano della Calabria and La Gazzetta del Sud), news agencies, news websites, weather forecast websites and the social network pages of people living in the affected areas diffused damage data.The event damaged 12.6% of the regional surface, affecting 49 municipalities (12% of the total), located in the provinces of Catanzaro (57%), Crotone (31%) and Vibo Valentia (12%).In 59% of the cases, even landslides caused some damage, and in 12%, wind also produced economic losses.In 90% of the affected municipalities, the road network was the most damaged element (Figure 4), followed by houses (33%); in two municipalities, 166 people were evacuated as a precautionary measure.Traffic, electricity and water supply were interrupted or difficult.Catanzaro, the administrative center of Calabria, was the most affected place: following the aqueduct breakdown, water service was interrupted for eight days, and water distribution for about 70,000 people had to be organized using tank trucks.Shops and houses were flooded; the university, the hospital, and the courthouse were flooded as well, and their activities had to be restricted for some days.Floods caused the loss of cultivated products and garbage dragged by rivers damaged fishing.

The Historical Series of Damaging Hydrogeological Events in Calabria
From the historical database of the Damaging Hydrogeological Events in Calabria (Rainfall measurements started in Calabria in 1920s), twelve of the most severe events which occurred in the region since the 1920s (Table 1) have been selected.
For each event, both the daily rainfall (Figure 5) and the return period of daily rainfall (Figure 6) were mapped, while, for the 2013 event, the same features are shown in Figure 2.For all these DHEs, both Dscore and Rscore (using the same T classes as in Figure 6) were assessed, and the chart of the severity of Calabrian DHEs was created (Figure 7).Table 1.Severe DHEs that have occurred in Calabria since the 1920s, chronologically sorted.N: identification number; DE Start-end: start-end of damage event; IDA (%): index of damaged area; DI: damage index; NoV: number of victims; dDE (days): duration of damage event; RE Start-end: start-end of rainfall event; dRE (days): duration of rainfall event; AdRE (days): average duration of rainfall event; mDR (mm): maximum daily rainfall of the event; ADR (mm): average daily rainfall of the event.(modified from: [58]).Figure 5. Maps of the daily rainfall recorded during the events.Numbers from 1 to 12 are the same as in Table 1.
Figure 6.Calabria shadow relief and maximum return period (T) of daily rainfall recorded during each event at the gauges (circles colored according to the legend).Numbers from 1 to 12 are the same as in Table 1.
In particular, in order to evaluate the Dscore, the values of IDA and DI have been classified in four classes.As regards the IDA, respectively a minimum and a maximum value of 2.83 and 27.84 have been observed.The difference between these values has been divided in four equal intervals and the following four classes have been detected:  Ordinary events represent 25% of the cases and show low values of NoV, which, when combined with the value of both IDA and DI, supplies Dscore values less than six.In regard to the Rscore values, these events show both long and short duration but the return periods are less than 200 years, with the exception of a gauge in the event No.8.In these events, rainfall mostly affect the central sector and the southeastern area of the region, from the Sila Plateau to the Aspromonte.
Among the extraordinary events, only a single concentrated summer event has been identified, characterized by very high return period rainfall (>200 years) which strongly damaged (DI = 84) a small sector of Calabria (IDA = 2.83%).In fact, rainfall mainly affected the city of Vibo Valentia and, partially, the regional plain and hilly sectors located between the Sila Plateau and the Serre Chain.
Catastrophic events represent 42% of the total, and high damage and relatively low Rscore (<6) characterize them.In this class, cases affecting either small (No.3) and large (No.6) regional sectors cause numerous victims and DI values usually over 100.Rainfall generally affects almost the entire region, particularly the east side in the central plain and hilly sectors between Sila Plateau and the Ionian Sea, and the area from the Aspromonte and Serre Chain to the Ionian Sea.
The Major catastrophic events present high values of both Rscore and Dscore.Nevertheless, the events of this group are heterogeneous, because they reach similar Dscore values due to the weight of the different damage indicators.Event No.1, for instance, affected a small sector with high DI and high NoV; No.5 affected a medium-sized sector, but caused the highest number of victims; No.9 caused a relatively low NoV and damage to a large percentage of the regional area.In general, rainfall mainly affects the east regional side, from Aspromonte to Serre Chain and the hilly areas of the Sila Plateau.
The 2013 event falls within the field of ordinary events.It presents both low Rscore and Dscore values because no casualties were recorded and the event affected a relatively restricted regional sector, with rainfall lasting less than one day.This case confirms the negative trend of heavy daily precipitation detected in Calabria [23].Despite the damage affecting the east side of the region, daily rainfall was not extraordinary: the highest T value was 20 years.On the contrary, high values of T were assessed for rainfall lasting 1, 3 and 6 h (Figure 8), equal to 50, more than 200 and more than 100 years, respectively.The impressive flood that affected Catanzaro (Figure 4) in part depends on the intense urbanization of those areas near the rivers.This process started in the 1950s, when population growth caused an increase in the demand for housing.Two maps of the town depicting both 1940 and 1990 configurations show the differences in the urbanized sectors near the rivers.
The S. Maria and Catanzaro Lido neighborhoods show the stronger modifications, with an urbanized area that increased almost tenfold in 50 years (Figure 9).During this period, weirs were built along the rivers, leading to a greater concentration of vulnerable elements in flood-prone areas.

Conclusions
The paper analyzed damaging hydrogeological events, defined as rainfall events able to trigger floods and landslides, thereby causing damage events.The characteristics of both triggering rainfall events and the consequent damage events strictly depend on climatic, geological, geomorphological and anthropogenic factors, which vary from region to region.
A methodology to classify severity of DHEs, based on both rainfall and damage indicators, was presented.A chart expressing the overall severity of the events was created using the values obtained from the indicators of damage and rainfall severity during past events.The scale of the chart is set depending on the values of Dscore and Rscore recorded in the study area.The chart is divided into four sectors: ordinary events, with the lowest values of both Dscore and Rscore; extraordinary events, with high values of Rscore but low values of Dscore; catastrophic events, during which even after non-exceptional rainfall (Rscore low) severe damage occurs (Dscore high); and major catastrophic events, characterized by the highest values of both Dscore and Rscore.
A practical application carried out on Calabria (Italy) highlighted the characteristics of regional events.According to this approach, the event that hit the region in 2013 has been classified as an ordinary event, when compared to the previous ones.The rainfall that triggered this event confirms the negative trend of heavy daily precipitation detected in Calabria, and the damage can be ascribed more to sub-daily than to daily rainfall affecting urbanized flood-prone areas.
This chart has been built using past events.However, the aim is not the analysis of past events.Rather, the objective is to produce a chart that will enable an objective classification of the severity of future events which media accounts and first or second-hand personal accounts fail to provide.
The restriction of the proposed approach is that the chart is representative of a specific area.However, using the proposed approach, when the climatic, geological, geomorphological and anthropogenic characteristics of another study area are available, the characteristic chart of that area can be built.

Figure 1 .
Figure 1.(A): digital elevation model of Calabria and rainfall gauges (circles); (B): regional distribution of the average annual rainfall (P).Small box: localization of the Calabria region (in black) in the Mediterranean Basin with the identification of the Sardinia region (in green) and of the Balearic Islands (in red).

Figure 2 .
Figure 2. 19 November 2013: (A) daily rainfall map; (B) return period of daily rainfall (T) at the gauges (circles with different colors according to the legend), and damaged municipalities (yellow areas).

Figure 4 .
Figure 4. Catanzaro town on 19 November 2014.Cars partially submerged by the floodwaters that overflowed from local rivers.

Figure 7 .
Figure 7. Chart of the severity of Calabrian DHEs.For the explanation of damage score (Dscore) and rainfall score (Rscore) see Section 2.3.Each circle represents one of the event listed in Table 1.(A) Ordinary events; (B) extraordinary events; (C) catastrophic events; (D) major catastrophic events.

Figure 8 .
Figure 8.Return period (T) of maximum hourly rainfalls with duration of (A) 1; (B) 3; and (C) 6 h registered on 19 November 2013 overlapping the map of the damaged municipalities.Different colors indicate different ranges of T, according to the legend.

Figure 9 .
Figure 9. Evolution of urbanized areas (striped areas) in the Catanzaro, from (A) 1940 to (B) 1990 (IGM maps 1940 and 1990).Black and blue lines are roads and rivers, respectively.