Reconstruction of a Storm Map and New Approach in the Definition of Categories of the Extreme Rainfall, Northeastern Sicily

After more than 350 mm of rainfall fell in a few hours on 22 November 2011, thousands of landslides and floods were induced in two main zones of Northeastern Sicily. The total rainfall has been reconstructed integrating available rain gauge data with Tropical Rainfall Measuring Mission (TRMM) satellite data from NASA (National Aeronautics and Space Administration); the landslide distribution in the field has confirmed the pattern of rainfall accumulated on 22 November 2011. Precipitation maxima of 1, 3, 6, 12, and 24 h was recognized as the hazardous events, which marks the evidence of a changing climate, with a shift toward more intense rainfalls in recent times. To investigate the sequence of the annual maxima, the historical time series have been transformed in the Standard normal distribution, from the cumulative probability of the GEV (Generalized Extreme Value) distribution. Following a similar definition of the Standard Precipitation Index (SPI), the transformation of the historical data in the standardized values allows the definition of categories of hourly maxima in term of extreme, severe, moderate, or mild. This transformation allows to eliminate the asymmetry of the time series, so that trends and fluctuations have been highlighted by the progressive accumulation of data (Rescaled Adjust Partial Sum). This statistical approach allows the improvement of the interpretability of the hydrological extreme events, and could also be used in other cases.


Introduction
Despite hydroclimate millennia reconstruction not supporting a general unprecedented intensification of the Northern Hemisphere hydrological cycle in the twentieth century associated with both more extreme wet and dry conditions [1], rainstorms and their injuries are becoming more dangerous as population and infrastructure continue to increase and to occupy areas exposed to flood risk [2,3].However, few literature sources are available worldwide regarding extreme precipitation and, especially, about rainstorm effects on regional and subregional terrestrial ecosystems and water resources [4][5][6][7].This also poses another question related to the development of dynamic hydrological models, hampered by incomplete understanding of spatially-varying processes and the lack of adequate datasets to spatially characterize varying rain inputs.According to [8], secular and more records of historical precipitation data are required to deal with long-term studies and cross-site comparisons.This is especially so in the Mediterranean region, where human pressure and erratic Water 2016, 8, 330 2 of 15 rainstorm patterns with marked inter-annual variability expose landforms to exacerbated, damaging hydrological processes [9][10][11][12].This is also so to spur the emergence of new hazards, such as coastal and urban flooding [13,14].Mediterranean rainstorms and cyclones tend to be characterized by short life cycles, with average radii ranging from 300 to 500 km (after [15]), many of which are a combination of both frontal and convective storms.Heavy flooding and storms occurring at Mediterranean sites were found to be characterized by a complex property, known as multifractality, which is the spatial distribution organized into clusters of high-rainfall localized cells embedded within a larger cloud system or clusters of lower intensity [16].
Sicily has been hit recently by intense rainstorms, which caused floods and landslides in many places.Particularly affected by recent and intense storms has been the northwestern sector, where a high and continuous ridge, the Peloritani Mountains, extends between the Ionian and Tyrrhenian seas.During recent years, the most catastrophic event occurred on 1 October 2009 along the Ionian side, with more than 200 mm recorded in a few hours, and hit the village of Giampilieri, where landslides caused 36 deaths [17,18].After 13 November 1855, the event of 2009 was the most catastrophic to have occurred in the Messina province, and follows other catastrophic events that have occurred recently.In particular, on 25 October 2007 an intense rainstorm induced many landslides and floods along the Ionian coast [19,20].
The storm of 22 November 2011 hit a wide area the Tyrrhenian sector of the Peloritani mountains, and caused three deaths in Saponara village (Figure 1).It induces thousands of shallow landslides in two distinct, wide areas (Figure 2), one located in the middle catchment of the Saponara torrent, and the other across both the catchments of Mela and Longano torrents.Huge alluvial phenomena occurred in the urban area of Barcellona Pozzo di Gotto and Milazzo.Recent landslides and flash floods hit the same area on 11 December 2008, 2 November 2010, and 10 October 2015.
Water 2016, 8, 330 2 of 15 rainstorm patterns with marked inter-annual variability expose landforms to exacerbated, damaging hydrological processes [9][10][11][12].This is also so to spur the emergence of new hazards, such as coastal and urban flooding [13,14].Mediterranean rainstorms and cyclones tend to be characterized by short life cycles, with average radii ranging from 300 to 500 km (after [15]), many of which are a combination of both frontal and convective storms.Heavy flooding and storms occurring at Mediterranean sites were found to be characterized by a complex property, known as multifractality, which is the spatial distribution organized into clusters of high-rainfall localized cells embedded within a larger cloud system or clusters of lower intensity [16].
Sicily has been hit recently by intense rainstorms, which caused floods and landslides in many places.Particularly affected by recent and intense storms has been the northwestern sector, where a high and continuous ridge, the Peloritani Mountains, extends between the Ionian and Tyrrhenian seas.During recent years, the most catastrophic event occurred on 1 October 2009 along the Ionian side, with more than 200 mm recorded in a few hours, and hit the village of Giampilieri, where landslides caused 36 deaths [17,18].After 13 November 1855, the event of 2009 was the most catastrophic to have occurred in the Messina province, and follows other catastrophic events that have occurred recently.In particular, on 25 October 2007 an intense rainstorm induced many landslides and floods along the Ionian coast [19,20].
The storm of 22 November 2011 hit a wide area the Tyrrhenian sector of the Peloritani mountains, and caused three deaths in Saponara village (Figure 1).It induces thousands of shallow landslides in two distinct, wide areas (Figure 2   In this study we estimate the rainfall that fell on 22 November 2011 by the rain gauge data available, integrated with satellite data; the time of the landslide occurrences has helped in an attempt to define the main features of the storm.
The historical time series of the annual maxima of rainfall, for different time duration, have been analyzed to estimate the frequency of extreme rainfall events and to highlight trend and fluctuations.

Main Physical and Geomorphological Characteristics of the Area Hit by the Storm
Due to its climate, affected by long, dry periods followed by heavy rainstorms, the Mediterranean region is particularly prone to multiple damaging hydrological events (MDHE; [21]) since it is characterized by high land vulnerability and poor vegetation coverage.From this point of view, Sicily is one the recognized areas particularly sensitive to changes in extreme precipitation events and currently threatened by flash-floods.The climate of Sicily, described in detail by [22], is the typical Mediterranean climate (Köppen's Csa type), with rainy winters and long, dry periods in summer.In Sicily, remarkable differences can be recognized between the southwestern portion, characterized by dry sub-humid to semi-arid conditions, and the northeastern portion, where subhumid/humid conditions prevail.The average annual temperatures range between 11 °C and 20 °C, with summer maxima higher than 40 °C, while rainfalls, mainly concentrated in autumn and winter (Figure 3), show mean annual values ranging from 385 mm up to 1200 mm and more.
The 22 November 2011 storm hit a wide area along the northern sector of the Peloritani Mountains, between Barcellona Pozzo di Gotto and Saponara villages.The main floods occurred in the Longano, Idria, Mela, and Saponara catchments (Figure 2), associated to diffuse landslide phenomena; however, other minor catchments were also affected by flash-floods and landslides, and the entire area of the Milazzo plain was affected by generalized flood phenomena.In this study we estimate the rainfall that fell on 22 November 2011 by the rain gauge data available, integrated with satellite data; the time of the landslide occurrences has helped in an attempt to define the main features of the storm.
The historical time series of the annual maxima of rainfall, for different time duration, have been analyzed to estimate the frequency of extreme rainfall events and to highlight trend and fluctuations.

Main Physical and Geomorphological Characteristics of the Area Hit by the Storm
Due to its climate, affected by long, dry periods followed by heavy rainstorms, the Mediterranean region is particularly prone to multiple damaging hydrological events (MDHE; [21]) since it is characterized by high land vulnerability and poor vegetation coverage.From this point of view, Sicily is one the recognized areas particularly sensitive to changes in extreme precipitation events and currently threatened by flash-floods.The climate of Sicily, described in detail by [22], is the typical Mediterranean climate (Köppen's Csa type), with rainy winters and long, dry periods in summer.In Sicily, remarkable differences can be recognized between the southwestern portion, characterized by dry sub-humid to semi-arid conditions, and the northeastern portion, where sub-humid/humid conditions prevail.The average annual temperatures range between 11 ˝C and 20 ˝C, with summer maxima higher than 40 ˝C, while rainfalls, mainly concentrated in autumn and winter (Figure 3), show mean annual values ranging from 385 mm up to 1200 mm and more.
The 22 November 2011 storm hit a wide area along the northern sector of the Peloritani Mountains, between Barcellona Pozzo di Gotto and Saponara villages.The main floods occurred in the Longano, Idria, Mela, and Saponara catchments (Figure 2), associated to diffuse landslide phenomena; however, other minor catchments were also affected by flash-floods and landslides, and the entire area of the Milazzo plain was affected by generalized flood phenomena.The area hit by the storm extends between the Tyrrhenian sea and water divide of the Peloritani Mountains.The coastal sector is characterized by a wide flat area, constituted by recent (Quaternary) alluvial and marine deposits, with ground elevation up to a few tens of meters above sea level (a.s.l.), and includes the Milazzo plain.Inland, ground elevation rapidly increases and substratum outcrops; here the slope angle can be higher than 45°, and characterizes the steep slopes of the Peloritani mountains.
Following the recent geological map of Italy [23], scale 1:50,000 (Map n. 600 "Barcellona Pozzo di Gotto" and Map n. 601 "Messina and Reggio Calabria"), the substratum belongs to a metamorphic complex of Paleozoic age, and to marine sedimentary complexes of Tertiary age (Flysch di Capo d'Orlando and Argille scagliose), both involved in Tertiary tectonic activity of the Apennine chain.Marine sinorogenetic deposits of Miocene (Formation of San Pier Niceto) and Pliocene-Pleistocene (Formation of Rometta) lie on the older terrains, and outcrop between the coastal and inland sectors.
A complex of normal faults have caused the general uplift of the area during the Pliocene and Pleistocene; marine deposits of the lower Pleistocene actually lie at 560 m a.s.l.(Rometta village), testifying the high uplift rate during the Quaternary.As a consequence, the hydrographic network is cached in the substratum, and determines steep and unstable slopes.
A diffuse soil mantle covers the slopes, constituted by colluvial and eluvial deposits originated from the weathering and erosion of the substratum; this soil mantle constitutes the main material involved in the landslide processes.

Main Features of the 22 November 2011 Storm
The meteorological conditions leading to 22 November 2011 storm were connected to an anomalous barometric gradient occurring between the northern side (Tyrrenian basin) and southwestern side (Ionian sea) of the Peloritani Mountains, inducing high-speed sirocco winds.These winds generally cause intense rainfall on the Ionian side of the Peloritani Mountains, where the highest accumulated rainfalls were expected the day before the storm.However, due to the very high intensity of winds during 22 November 2011, the rainfall accumulated mainly on the other side of the Peloritani water divide, causing unexpected landslides and flash floods along the Tyrrhenian side of the Peloritani Mountains.This storm behavior is known locally as the "Alcantara effect", because the wind from the south rises along the valley between Etna volcano and Peloritani Mountains and crosses the watershed of Peloritani thanks also to strong winds in altitude [24].
The storm few hours later reached the Southern Calabria, inducing further landsides and flash floods; a train derailed near Lamezia Terme, fortunately without further deaths.The area hit by the storm extends between the Tyrrhenian sea and water divide of the Peloritani Mountains.The coastal sector is characterized by a wide flat area, constituted by recent (Quaternary) alluvial and marine deposits, with ground elevation up to a few tens of meters above sea level (a.s.l.), and includes the Milazzo plain.Inland, ground elevation rapidly increases and substratum outcrops; here the slope angle can be higher than 45 ˝, and characterizes the steep slopes of the Peloritani mountains.
Following the recent geological map of Italy [23], scale 1:50,000 (Map n. 600 "Barcellona Pozzo di Gotto" and Map n. 601 "Messina and Reggio Calabria"), the substratum belongs to a metamorphic complex of Paleozoic age, and to marine sedimentary complexes of Tertiary age (Flysch di Capo d'Orlando and Argille scagliose), both involved in Tertiary tectonic activity of the Apennine chain.Marine sinorogenetic deposits of Miocene (Formation of San Pier Niceto) and Pliocene-Pleistocene (Formation of Rometta) lie on the older terrains, and outcrop between the coastal and inland sectors.
A complex of normal faults have caused the general uplift of the area during the Pliocene and Pleistocene; marine deposits of the lower Pleistocene actually lie at 560 m a.s.l.(Rometta village), testifying the high uplift rate during the Quaternary.As a consequence, the hydrographic network is cached in the substratum, and determines steep and unstable slopes.
A diffuse soil mantle covers the slopes, constituted by colluvial and eluvial deposits originated from the weathering and erosion of the substratum; this soil mantle constitutes the main material involved in the landslide processes.

Main Features of the 22 November 2011 Storm
The meteorological conditions leading to 22 November 2011 storm were connected to an anomalous barometric gradient occurring between the northern side (Tyrrenian basin) and southwestern side (Ionian sea) of the Peloritani Mountains, inducing high-speed sirocco winds.These winds generally cause intense rainfall on the Ionian side of the Peloritani Mountains, where the highest accumulated rainfalls were expected the day before the storm.However, due to the very high intensity of winds during 22 November 2011, the rainfall accumulated mainly on the other side of the Peloritani water divide, causing unexpected landslides and flash floods along the Tyrrhenian side of the Peloritani Mountains.This storm behavior is known locally as the "Alcantara effect", because the wind from the south rises along the valley between Etna volcano and Peloritani Mountains and crosses the watershed of Peloritani thanks also to strong winds in altitude [24].The storm few hours later reached the Southern Calabria, inducing further landsides and flash floods; a train derailed near Lamezia Terme, fortunately without further deaths.
Figure 4 shows the development of the atmospheric disturbance from satellite, where, from 6:00, a plume begins to develop, with its vertex on Northeastern Sicily and it spreads to the southern part of the Italian peninsula.The storm has a typical enhanced-V shape from the satellite images.
Water 2016, 8, 330 5 of 15 Figure 4 shows the development of the atmospheric disturbance from satellite, where, from 6:00, a plume begins to develop, with its vertex on Northeastern Sicily and it spreads to the southern part of the Italian peninsula.The storm has a typical enhanced-V shape from the satellite images.Figure 2 shows the rain gauge network available during the storm and Table 1 the rainfall accumulated.The highest rainfall was recorded at Castroreale where 351 mm were recorded; with the exception of Torregrotta, with 188.8 mm, the other rain gauges recorded lower rainfall, especially for those located along the Ionian side of the Peloritani Mountains.In Messina town, only a few millimeters of rainfall were recorded.However, the distribution of accumulated rainfall of Figure 5 cannot explain the diffuse landslide phenomena and flash floods occurred in the Saponara catchment, but only that of Longano, Idria, and Mela catchments (cf. Figure 2).This is due to low rain gauge density distribution, which is unable to record the areal peaks distribution of the rainfall, connected to storm development (Figure 4).Additionally, in the area of Saponara, landslides developed during the afternoon of 22 November 2011, several hours later than those occurring in the Mela and Longano catchments.These characteristics indicate that storm migrated toward the northeast, according to the satellite images.
Figure 6 shows the accumulated hourly rainfall recorded at rain gauges during the 22 November 2011 storm.The Castroreale rain gauge recorded the most intense rainfall; it began to rain since the early hours of 22 November, but between 6:00 and 15:00 the storm caused the major part of the precipitation.In particular, 72.5 mm fell between 7:00 and 8:00 and caused the first flood phenomena in the minor gullies.Towards the northeast, (Torregrotta, San Pier Niceto, and Ziriò) the most intense rainfall occurred one or few hours later.Many eyewitnesses described that, in the area of Saponara, the major precipitation occurred in the afternoon of 22 November 2011, highlighting that mainly in this zone the rain gauge network was insufficient to catch the spatial distribution of the rainfall.Figure 2 shows the rain gauge network available during the storm and Table 1 the rainfall accumulated.The highest rainfall was recorded at Castroreale where 351 mm were recorded; with the exception of Torregrotta, with 188.8 mm, the other rain gauges recorded lower rainfall, especially for those located along the Ionian side of the Peloritani Mountains.In Messina town, only a few millimeters of rainfall were recorded.However, the distribution of accumulated rainfall of Figure 5 cannot explain the diffuse landslide phenomena and flash floods occurred in the Saponara catchment, but only that of Longano, Idria, and Mela catchments (cf. Figure 2).This is due to low rain gauge density distribution, which is unable to record the areal peaks distribution of the rainfall, connected to storm development (Figure 4).Additionally, in the area of Saponara, landslides developed during the afternoon of 22 November 2011, several hours later than those occurring in the Mela and Longano catchments.These characteristics indicate that storm migrated toward the northeast, according to the satellite images.
Figure 6 shows the accumulated hourly rainfall recorded at rain gauges during the 22 November 2011 storm.The Castroreale rain gauge recorded the most intense rainfall; it began to rain since the early hours of 22 November, but between 6:00 and 15:00 the storm caused the major part of the precipitation.In particular, 72.5 mm fell between 7:00 and 8:00 and caused the first flood phenomena in the minor gullies.Towards the northeast, (Torregrotta, San Pier Niceto, and Ziriò) the most intense Water 2016, 8, 330 6 of 15 rainfall occurred one or few hours later.Many eyewitnesses described that, in the area of Saponara, the major precipitation occurred in the afternoon of 22 November 2011, highlighting that mainly in this zone the rain gauge network was insufficient to catch the spatial distribution of the rainfall.

Methods Used
To reconstruct a more likely spatial distribution of rainfall occurring on 22 November 2011, rain gauge density has been integrated with data recorded from the Tropical Rainfall Measuring Mission satellite (TRMM).TRMM is a satellite launched by NASA for research purposes and provides rainfall data as a function of the heat release which manifests contextually to precipitation.Data download is through a portal that can be accessed through the NASA website, with an acquisition interval of 95.2 min.Acquisition points of TRMM provide coverage to the entire Earth's surface, they form a regular grid of dots arranged at intervals of 0.25° of latitude and longitude.
Available rain gauge data series on the annual maxima rainfall (1, 3, 6, 12, and 24 h) have been analyzed to find some statistical characteristics, as the return time and trend in the time series.
The return time has been estimated by the Generalized Extreme Value (GEV) [25] frequency distribution, where the cumulative probability, P, for values below a specific threshold (X ≤ x), is: where ξ is the shape parameter (dimensionless), μ and σ are the position and scale parameters.
The actual values were plotted by their order in the time series: where N is the number of the series data, and i is the position in the ordered set of data.The cumulative probability found by the GEV frequency distribution (Equation ( 1)) has been transformed in the Z value of the Standard normal distribution (Figure 7).This transformation has been carried out for the annual maxima series of rainfall (1, 3, 6, 12, and 24 h), and obtains the symmetric series, with zero mean and standard deviation equal to 1.After the transformation in the normal standardized distribution, the progressive summation of the values highlights trends and fluctuations of the time series; this summation coincides with the Rescaled Adjusted Partial Sum (RAPS; [26]).The RAPS of the time series are defined as follows: where is sample mean, SY is the standard deviation, k is counter limit of the current summation.The plot of the RAPS versus time is the visualization of the trends and fluctuations of Yt [26].Decreasing patterns in the RAPS are the result of periodicity of mostly below-average Yt values,

Methods Used
To reconstruct a more likely spatial distribution of rainfall occurring on 22 November 2011, rain gauge density has been integrated with data recorded from the Tropical Rainfall Measuring Mission satellite (TRMM).TRMM is a satellite launched by NASA for research purposes and provides rainfall data as a function of the heat release which manifests contextually to precipitation.Data download is through a portal that can be accessed through the NASA website, with an acquisition interval of 95.2 min.Acquisition points of TRMM provide coverage to the entire Earth's surface, they form a regular grid of dots arranged at intervals of 0.25 ˝of latitude and longitude.
Available rain gauge data series on the annual maxima rainfall (1, 3, 6, 12, and 24 h) have been analyzed to find some statistical characteristics, as the return time and trend in the time series.
The return time has been estimated by the Generalized Extreme Value (GEV) [25] frequency distribution, where the cumulative probability, P, for values below a specific threshold (X ď x), is: where ξ is the shape parameter (dimensionless), µ and σ are the position and scale parameters.The actual values were plotted by their order in the time series: where N is the number of the series data, and i is the position in the ordered set of data.The cumulative probability found by the GEV frequency distribution (Equation ( 1)) has been transformed in the Z value of the Standard normal distribution (Figure 7).This transformation has been carried out for the annual maxima series of rainfall (1, 3, 6, 12, and 24 h), and obtains the symmetric series, with zero mean and standard deviation equal to 1.After the transformation in the normal standardized distribution, the progressive summation of the values highlights trends and fluctuations of the time series; this summation coincides with the Rescaled Adjusted Partial Sum (RAPS; [26]).The RAPS of the time series are defined as follows: where Y is sample mean, S Y is the standard deviation, k is counter limit of the current summation.The plot of the RAPS versus time is the visualization of the trends and fluctuations of Y t [26].Decreasing patterns in the RAPS are the result of periodicity of mostly below-average Y t values, whereas increasing patterns are the result of periods of mostly above-average Y t values.The RAPS highlights the tendency of values to cluster together; that is, the consecutive values below or above the mean provides a decreasing or increasing trend in the RAPS, respectively.
Water 2016, 8, 330 8 of 15 whereas increasing patterns are the result of periods of mostly above-average Yt values.The RAPS highlights the tendency of values to cluster together; that is, the consecutive values below or above the mean provides a decreasing or increasing trend in the RAPS, respectively.

Reconstruction of the Storm Map
To improve the rainfall spatial distribution on 22 November 2011, rain gauge data have been integrated with data from the TRMM satellite platform.The original values provided by TRMM has been multiplied be a constant, obtained from the ratio between rainfall recorded at several rain gauges and their nearest TRMM measure points; the constant found has a value of 6.On the network formed by rain gauges and TRMM points, a new spatial rainfall distribution has been mapped by kriging interpolation.Figure 8 depicts probably a more realistic spatial distribution of rains from 22 November 2011 with respect to Figure 5; it seems to agree with what was derived from EUMETSAT satellite images of Figure 4, even if it is also based on rainfall data derived indirectly from satellite sensors.The map shows that the storm event of 22 November 2011 is characterized by two rainy centers, elongated in the direction of Southwest-Northeast, one located transversely to Longano, Idria and Mela catchments, and the other in the area of Saponara catchment.The result of interpolation of the data in Figure 8, with the two centers of rain, is in fact the result of the stationing of the plume vertex (mentioned before) for most of the time on these areas.It can be observed that, in the area of Milazzo plain at least 200 mm of rain, would have fallen, according to flash-floods occurring in this area.However, the distribution of rainfall in the Saponara area appears still be coarse (locally, it is also missing a TRMM point, Figure 8), as suggested by the landslide distribution in this zone (Figure 2).

Reconstruction of the Storm Map
To improve the rainfall spatial distribution on 22 November 2011, rain gauge data have been integrated with data from the TRMM satellite platform.The original values provided by TRMM has been multiplied be a constant, obtained from the ratio between rainfall recorded at several rain gauges and their nearest TRMM measure points; the constant found has a value of 6.On the network formed by rain gauges and TRMM points, a new spatial rainfall distribution has been mapped by kriging interpolation.Figure 8 depicts probably a more realistic spatial distribution of rains from 22 November 2011 with respect to Figure 5; it seems to agree with what was derived from EUMETSAT satellite images of Figure 4, even if it is also based on rainfall data derived indirectly from satellite sensors.The map shows that the storm event of 22 November 2011 is characterized by two rainy centers, elongated in the direction of Southwest-Northeast, one located transversely to Longano, Idria and Mela catchments, and the other in the area of Saponara catchment.The result of interpolation of the data in Figure 8, with the two centers of rain, is in fact the result of the stationing of the plume vertex (mentioned before) for most of the time on these areas.It can be observed that, in the area of Milazzo plain at least 200 mm of rain, would have fallen, according to flash-floods occurring in this area.However, the distribution of rainfall in the Saponara area appears still be coarse (locally, it is also missing a TRMM point, Figure 8), as suggested by the landslide distribution in this zone (Figure 2).

Frequency Analyses
Specific to the area hit by the 22 November 2011 storm, the rain gauge at Castroreale is the only one where a long historical time series of maxima exists (since 1930), even if some missing data are present.This rain gauge is located in one of two zones characterized by the maximum amount of the rainfall that fell during the storm, as can be deduced from Figure 8 and the landslide distribution occurring in the field (Figure 2).Thus, the maximum intensity of the storm can be analyzed statistically, and provides details on the recurrence of the intense storms in this specific point.
Even if a considerable accumulated rainfall was also recorded at the Torregrotta rain gauge, it has only operated since 2002.Table 2 shows that the rainfall recorded at the Castroreale rain gauge during the 22 November 2011 was the highest for each hourly duration (from 1 to 24 hourly rainfall).In particular, the rainfall of 6 and 12 h has been well higher than previous historical maxima (Table 2).Figure 9 shows the statistical analyses carried out on annual maxima for 1, 3, 6, 12, and 24 cumulative rainfalls, plotting results as described in [27].If annual maxima of 2011 are excluded from the dataset (Figure 9A), each statistical curve (GEV) fits the historical data well; the return time of the 22 November 2011 would reach values of 450-600 years between three and 12 h of accumulated rainfall.
If the annual maxima of 2011 are included in the dataset (Figure 9B), the statistical curves (GEV) are not able to fit the extreme value of the series for 3, 6, and 12 h cumulative rainfall (Figure 9B); the return time of the 22 November 2011 storm fall to 160-175 years between 3 and 12 h accumulated rainfall.Thus, the presence of annual maxima of 2011 in data set conditions heavily the estimation of the return time.

Frequency Analyses
Specific to the area hit by the 22 November 2011 storm, the rain gauge at Castroreale is the only one where a long historical time series of maxima exists (since 1930), even if some missing data are present.This rain gauge is located in one of two zones characterized by the maximum amount of the rainfall that fell during the storm, as can be deduced from Figure 8 and the landslide distribution occurring in the field (Figure 2).Thus, the maximum intensity of the storm can be analyzed statistically, and provides details on the recurrence of the intense storms in this specific point.
Even if a considerable accumulated rainfall was also recorded at the Torregrotta rain gauge, it has only operated since 2002.Table 2 shows that the rainfall recorded at the Castroreale rain gauge during the 22 November 2011 was the highest for each hourly duration (from 1 to 24 hourly rainfall).In particular, the rainfall of 6 and 12 h has been well higher than previous historical maxima (Table 2).Figure 9 shows the statistical analyses carried out on annual maxima for 1, 3, 6, 12, and 24 cumulative rainfalls, plotting results as described in [27].If annual maxima of 2011 are excluded from the dataset (Figure 9A), each statistical curve (GEV) fits the historical data well; the return time of the 22 November 2011 would reach values of 450-600 years between three and 12 h of accumulated rainfall.
If the annual maxima of 2011 are included in the dataset (Figure 9B), the statistical curves (GEV) are not able to fit the extreme value of the series for 3, 6, and 12 h cumulative rainfall (Figure 9B); the return time of the 22 November 2011 storm fall to 160-175 years between 3 and 12 h accumulated rainfall.Thus, the presence of annual maxima of 2011 in data set conditions heavily the estimation of the return time.
Water 2016, 8, 330 10 of 15 The exceptional rainfall of the 22 November 2011 storm caused a considerable shift from each frequency distribution adopted to fit the historical data, and raises the well-known complication in this type of statistical analysis [28,29].This aspect highlights the outlier characteristic of the 2011 event, and will be considered later.

Transformation in the Standard Normal Distribution
To investigate the sequence of the annual maxima, the historical time series have been transformed in the Standard normal distribution, from the cumulative probability of the GEV distribution (Figure 10).The transformed series of Figure 10 are symmetric with respect to the zero value of the mean and have a standard deviation equal 1.The exceptional rainfall of the 22 November 2011 storm caused a considerable shift from each frequency distribution adopted to fit the historical data, and raises the well-known complication in this type of statistical analysis [28,29].This aspect highlights the outlier characteristic of the 2011 event, and will be considered later.

Transformation in the Standard Normal Distribution
To investigate the sequence of the annual maxima, the historical time series have been transformed in the Standard normal distribution, from the cumulative probability of the GEV distribution (Figure 10).The transformed series of Figure 10 are symmetric with respect to the zero value of the mean and have a standard deviation equal 1.With the exception of 1978, which appears randomly distributed in the time series, the most intense storms are concentrated in recent years (Figure 10), especially considering durations of 1, 3, and 6 h.In particular, large storms occurred in 2008 and 2010, and, locally, induced landslides and floods, but the rainfall accumulated for different durations was lower than the 2011 storm.For durations of 12 and 24 h, the storms of 1943 and 1949 appear intense in the historical series.
In the transformed series of Figure 10, the values that exceed 2 have a probability of occurrence lower than 0.023; as commonly used in hydrology, for example in the definition of drought categories by the SPI (Standard Precipitation Index; [30]), the values that exceed 2 can be considered as extreme events (in drought analyses the negative values of Z are considered).Following a similar definition of the SPI, but considering positive values of Z in this case, the categories of hourly maxima are defined as extreme, severe, moderate, or mild.The transformation in the Standard normal distribution provides dimensionless values, and allows a direct comparison between the time series of Figure 10 (1, 3, 6, 12 and 24 h).Besides, it also allows the comparison of different hydrological parameters, such as spring discharge and rainfall [31].The transformed historical data of Figure 10 have been plotted in Figure 11, fitted by the Standard normal distribution.Values are defined as extremes for Z ≥ 2, severe for 1.5 ≤ Z < 2, moderate for 1 ≤ Z < 1.5, and mild for 0 ≤ Z < 1.The Z = 2.43 limit has been found following the Chauvenet's method, and would allow to identify the outliers in With the exception of 1978, which appears randomly distributed in the time series, the most intense storms are concentrated in recent years (Figure 10), especially considering durations of 1, 3, and 6 h.In particular, large storms occurred in 2008 and 2010, and, locally, induced landslides and floods, but the rainfall accumulated for different durations was lower than the 2011 storm.For durations of 12 and 24 h, the storms of 1943 and 1949 appear intense in the historical series.
In the transformed series of Figure 10, the values that exceed 2 have a probability of occurrence lower than 0.023; as commonly used in hydrology, for example in the definition of drought categories by the SPI (Standard Precipitation Index; [30]), the values that exceed 2 can be considered as extreme events (in drought analyses the negative values of Z are considered).Following a similar definition of the SPI, but considering positive values of Z in this case, the categories of hourly maxima are defined as extreme, severe, moderate, or mild.The transformation in the Standard normal distribution provides dimensionless values, and allows a direct comparison between the time series of Figure 10 The 1943 and 1949 storms were severe only for 12 and 24 h, indicating that the hourly distribution of these storms were characterized by the absence of particular peaks.In contrast, some recent storms (1995 and 2004) were severe only for a one hour duration, and they lost their power for a longer duration.
As the annual maximum hourly rainfall time series have been standardized, their progressive accumulation provides the RAPS of the series (Figure 10).Looking at the RAPS of the annual series of Figure 10, the clustering appears in two main periods, the first occurred between 1965 and 1974 with values below the mean; the second period has occurred since 2007 and it is characterized by values well above the mean.
The latter clustering highlights the general increasing of the occurrence of the intense storms, which seems common in all northern parts of Sicily.This is in agreement with the results of [32], by which sub-grid scale convection and intensification phenomena indicate, for Sicily, an increasing trend towards shorter rain durations (about one hour) during the period 1927-2004.Additionally, the increase in the intense precipitation has been recently highlighted in Western Sicily by [33].
The recent increasing of the intense rainfall, marked by the inconstant mean and variance in time series, highlights the difficulties in the estimation of the frequency of recent storms.

Conclusions
Placed at the crossroads of Mediterranean flows, Sicily could be considered representative of a geographic scenario of regional climate change associated with warming, reflected in the Mediterranean cyclogenesis, and that would have contributed to the exceptional rainfall rates observed in recent times [2,15].The assessment of current and future management systems should focus on the hazard of extreme precipitation events, which likely cause accelerated urban and rural land degradation.This must be considered together with the findings that, although rainfall amounts are not always increasing, erratic spatial and temporal storm patterns in some seasons or months drive the power of rainfall to increase its hazard.According to [34], investigations are needed especially in terms of statistical analysis of long pluviometric series, rain intensity and erosivity occurring at decadal and interdecadal scales, which mark substantial climatic fluctuations and changes during last centuries.However, the increase of extremes of precipitation cannot be assumed as a rule in Southern Italy, and in the Calabria region this increase was not observed [35].
Specific to the investigated area, on the basis of the landslide distribution in the field, and on the numerous eyewitnesses describing the storm of 22 November 2001, it clearly appears that the gauge network is insufficient to define the map of intense storms.Thus, some historical storms could be unrecorded or underestimated by the actual rain gauge network, as well as floods occurring in the ungauged catchment of Mela and Saponara.
Historical data available for the Castroreale rain gauge, the only one useful for a statistic, have shown the intensification of the powerful storms in this area, as occurred in 2008, 2010, 2011, and 2015.This intensification could be connected to new and different atmospheric circulation, and the consequent develop of the "Alcantara effect" in this area; before 2008, on the basis of the historical records, this phenomenon seems to have occurred more rarely in the past.
To estimate the probability of occurrence of intense storms, the transformation of the historical data in the standardized values, passing through their best distribution fit (GEV in our case), has allowed to define several categories of the intensity of storms, similarly to SPI.
Even if the effect of storms in the field (floods and landslides) strongly depends on the previous wet conditions of the soil that is on the antecedent rainfall [36], the power of storms appear useful to be expressed in term of the Z value.
This approach can help in the evaluation of extreme data, as the use of the Standard normal distribution allows a more simple statistical management of data.Trends and fluctuations in time series need specific considerations, because the effect of the recent increase in frequency of intense storms could cause an overestimation of the return time values.
), one located in the middle catchment of the Saponara torrent, and the other across both the catchments of Mela and Longano torrents.Huge alluvial phenomena occurred in the urban area of Barcellona Pozzo di Gotto and Milazzo.Recent landslides and flash floods hit the same area on 11 December 2008, 2 November 2010, and 10 October 2015.

Figure 1 .
Figure 1.(A) The mouth of Longano River on the morning of 23 November 2011, and the bridge collapsed after the flood; (B) tree trunks obstruct the river flow in town center of Barcellona Pozzo di Gotto; (C) debris flows hit Varela village, Mela river catchment; and (D) a debris avalanche struck houses in Scarcelli village, Saponara municipality, and caused three deaths.

Figure 1 .Figure 2 .
Figure 1.(A) The mouth of Longano River on the morning of 23 November 2011, and the bridge collapsed after the flood; (B) tree trunks obstruct the river flow in town center of Barcellona Pozzo di Gotto; (C) debris flows hit Varela village, Mela river catchment; and (D) a debris avalanche struck houses in Scarcelli village, Saponara municipality, and caused three deaths.

Figure 2 .
Figure 2. The northeastern sector of Sicily, characterized by the Peloritani ridge, extended between the Tyrrhenian and Ionian sea; (1) rain gauge network working and rainfall recorded on 22 November 2011; (2) catchment boundary of the Peloritani ridge, Tyrrhenian side; (3) catchment of Longano (A), Idria (B), Mela (C) and Saponara (D) torrent affected by floods on 22 November 2011; and (4) the area hit by landslides on 22 November 2011.Damage occurred mainly in the villages of Barcellona Pozzo di Gotto and Saponara, which are highlighted by the white circles.
Apr May Jun Jul Aug Sep Oct Nov Dec 50°percentile 95°percentile monthly rainfall, mm

Figure 3 .
Figure 3. Seasonal regime of different percentile values for the Messina municipality (mean data, 1965-1994 period).

Figure 4 .
Figure 4. Satellite imagines (EUMETSAT, IR 10.8), every 2 h of the 22 November 2011, between 03:00 and 19:00.The white zones are the colder clouds, and the enhanced-V shape develops on Northeastern Sicily.

Figure 4 .
Figure 4. Satellite imagines (EUMETSAT, IR 10.8), every 2 h of the 22 November 2011, between 03:00 and 19:00.The white zones are the colder clouds, and the enhanced-V shape develops on Northeastern Sicily.

Figure 6 .
Figure 6.Hourly rainfall of 22 November 2011 for several rain gauges.

Figure 6 .
Figure 6.Hourly rainfall of 22 November 2011 for several rain gauges.

Figure 7 .
Figure 7. Transformation from GEV distribution to Standard normal distribution.

Figure 7 .
Figure 7. Transformation from GEV distribution to Standard normal distribution.

Figure 9 .
Figure 9.Return time for different rainfall durations (Castroreale rain gauge).Historical data for 1, 3, 6, 12, and 24 h (white and blue circles) are fitted by the GEV (cumulative) frequency distribution.Triangles show the return times of 22 Novenber 2011 on the GEV curves.(A) Analysis with n.65 data for the 1930-2010 period (without 2011 data); and (B) analysis with n.66 data of the 1930-2011 period (with 2011 data).

Figure 9 .
Figure 9.Return time for different rainfall durations (Castroreale rain gauge).Historical data for 1, 3, 6, 12, and 24 h (white and blue circles) are fitted by the GEV (cumulative) frequency distribution.Triangles show the return times of 22 Novenber 2011 on the GEV curves.(A) Analysis with n.65 data for the 1930-2010 period (without 2011 data); and (B) analysis with n.66 data of the 1930-2011 period (with 2011 data).

Table 1 .
Figure 5. 22 November 2011 rainfall reconstructed by kriging interpolation on available rain gauge ofFigure 5. 22 November 2011 rainfall reconstructed by kriging interpolation on available rain gauge of Table1.

Table 1 .
Rainfall recorded on 22 November 2011 by the rain gauges network.

Table 2 .
Historical maxima recorded at the Castroreale rain gauge, for different time intervals.

Table 2 .
Historical maxima recorded at the Castroreale rain gauge, for different time intervals.

Table 10 .
Annual maximum hourly rainfall, 1, 3, 6, 12, and 24 h time series (blue circles), expressed as Z value of the Standard normal distribution; progressive accumulation (RAPS) is also plotted (black line).