Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 1

Copyright Ó 2010, Paper 14-004; 46114 words, 7 Figures, 0 Animations, 0 Tables. http://EarthInteractions.org

Storminess and Environmental Changes in the Mediterranean Central Area

Nazzareno Diodato and Gianni Bellocchi* Met European Research Observatory, GEWEX-CEOP Network, World Climate Research Programme, Benevento,

Received 11 October 2009; accepted 7 April 2010

ABSTRACT: Earth ecosystems are not static, and they respond to environ- mental changes, particularly climatic and anthropogenic. Precipitation varying in its extremeness, with shifts to greater or lesser intensity of individual storms and/or to change in the length and frequency of wet and dry periods, can ad- versely affect both urban and rural ecosystems. Here, the authors review long- term precipitation records of the central Mediterranean area and employ a Web geographical information system (GIS)-based analytical approach to compare current rainfall impact with historical data on different spatial and temporal scales. Autumn (September–November) was recognized as the most hazardous season that marks the evidence of a changing climate, with a shift toward more intense rainfalls in recent times. In the first decade of the third millennium, areas of peninsular and insular Italy have been especially affected by extreme rains. A focus was put on the island of Sicily, where extraordinary rain events occurred in September 2009, discussed in the context of upcoming trends and climate histories. An improved knowledge and understanding of the scale at which changes on extremes occur is essential for dealing with the forthcoming challenges regarding soil and water conservation practices. The characteristics of changes in natural rainfall, its role on terrestrial ecosystems, and its effect on

* Corresponding author address: G. Bellocchi, Met European Research Observatory, via Monte Pino, 82100 Benevento, Italy. E-mail address: [email protected]

DOI: 10.1175/2010EI306.1

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 2 surface water erosion dynamics are discussed. It is argued that understand- ing these issues are major priorities for future research to promote a better understanding of the Earth interaction with water resources and related hydro- logical issues. KEYWORDS: Climate change; Extreme rainfall; Mediterranean central area

1. Introduction A march of waters; mustering from above, Ne` di rado dal ciel scende dirotto, D’acque il The clouds roll up the tempest, heaped and diluvio: in procelloso nembo; grim Vedi repente le aggruppate nubi, Fosche With angry showers: down falls the height addensarsi, e con opaco velo; of heaven, Tutto il mondo oscurar. Disciolto in acqua, And with a great rain floods the smiling Sembra cadere il ciel; scorre a crops, torrenti, La pioggia, e inonda i seminati, e Theoxen’slabour:nowthedikesfillfast, strugge Le fatiche de’ buoi: And the void river-beds swell thunderously, s’empion le fosse, Gonfiansi i fiumi, e And all the panting firths of Ocean boil. la sonante piena Cacciano al mar (...). Publio Virgilio Marone, Georgiche, Book I, written 29 BCE There is growing evidence at global, regional, and local scales that intra-annual precipitation regimes have already become more extreme (Easterling et al. 2000; Groisman et al. 2005; Knapp et al. 2008). However, the uncertainty of climate information poses challenges for the analysis of observed rain data because the heaviest areas of precipitation may fall between recording stations (Willmott and Legates 1991). For instance, assertions about extreme precipitation changes may be more reliable for regions with dense networks because of the small radius of correlation for many intense rainfall events (Groisman et al. 2005). Few literature sources are available worldwide regarding the extreme precipitation, especially about rainstorm effects on terrestrial ecosystems and water resources (e.g., Clarke and Rendell 2007; Curtis et al. 2007; Zolina et al. 2009). This also poses another question related to the development of dynamic hydrological models, which is hampered by incomplete understanding of spatially varying processes and the lack of adequate datasets to spatially characterize varying rain inputs. According to Wei et al. (Wei et al. 2009), secular 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 rain- fall patterns with marked interannual variability expose landforms to exacerbated, damaging hydrological processes (van Leeuwen and Sammons 2003; Sa´nchez et al. 2004; van Rompaey et al. 2005) and also to spur the emergence of new hazards, such as coastal and urban flooding (Papathoma and Dominey-Howes 2003; Barroca et al. 2006). Mediterranean storms and cyclones tend to be characterized by short life cycles, with average radius ranging from 300 to 500 km (after Lionello et al. 2006), many of which are a combination of both frontal and convective rainstorms. 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 (Mazzarella 1999).

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 3

Figure 1. (a) Economic losses from damaging hydrological events worldwide (from APFM) and (b) expected value of the expected cost as a function of dif- ferent daily rainfall threshold values (Martina et al. 2006; modified).

For the Holocene—the geological epoch started around 10 000 BCE and con- tinues to the present—there is continued debate about the relative impact of anthropogenic activity, but there is also increasing recognition that Mediterranean- type ecosystems should not be regarded as fragile, degraded landscapes but instead as disturbance-adapted systems (Allen 2003). Nevertheless, conservation measures face increasing challenges from contemporary human forcing and rainstorm events. The increase in sealed surfaces resulting from changes in land use together with a decrease in forest cover has increased the frequency and size of storm runoff too, causing flooding, mudflows, and landslides. These events may be grouped in some particularly rainy years or months according to storms climatic variability over interannual to century scales (Garcı´a-Oliva et al. 1995; D’Odorico et al. 2001; Peterson et al. 2002; Cavazos and Rivas 2004). Rainfall variables such as depth, duration, and intensity (and its erosive power) are becoming more and more changeable (Diodato and Bellocchi 2009a) and with a regimen that would still markedly change over time and space in the following years (Richard 2007). The environmental and economic impact of these events can be high in several regions (Kunkel et al. 1999; Alca´ntara-Ayala 2002), especially in agricultural and river-torrential areas (Thornes and Alca´ntara-Ayala 1998; Camarasa Belmonte and Segura Beltra´n 2001; Ramos and Mulligan 2003). The European southernmost regions are especially more sensitive to erratic rainfalls and are currently threat- ened by land degradation leading to a lowering in water resource availability and agricultural productivity (Poesen and Hooke 1997). Over the 1980–90 decade, for example, flash floods caused billions of euros of damage in Europe (Gruntfest and Handmer 2001). This is also evident with more recent data that show (Figure 1a) the extent to which economic losses from worldwide hydrological disasters tend to increase [World Meteorological Organization Advanced Programme on Flood Management (APFM); available online at http://www.apfm.info/index.htm]. In Figure 1b, one can see instead how the shape of expected-value cost of disasters affected by a given cumulated volume of rainfall during a storm event (i.e., 50 mm day21) separates decreasing cost from increasing cost. This is so because dry periods are associated with expected-value costs related to environmental drought

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 4 stress, whereas excessive rains (wet periods) turn into expected-value costs associ- ated with multiple damaging hydrological events. Threshold behavior manifests itself as an indicator of economic losses, and its characteristics are of importance for understanding and extrapolating the dynamics and stability of climate systems (Pitman and Stouffer 2006), the dynamics and resilience of geoecosystems, and the dynamics of fluxes in hydrologic systems (Zehe and Sivapalan 2009). Despite the general recognition that land degradation is a serious and widespread problem in Mediterranean countries, the same has not been quantified for many locations, and its geographical distribution and real extent are not accurately known, because pre- vailing studies only target subregional scales or isolated places (Clarke and Rendell 2007). There is a need to update and estimate the current overall hazard related to rainstorms and to assess in what measure hazardous rainstorms evolved with cli- mate variability. The present work aims at prospecting, for the central Mediter- ranean and for multiple spatial and temporal scales, rain-intensity trends related to recorded extreme rainfalls. To achieve this goal, precipitation data were extracted from the Web with the extension of Internet-based geographic information system (Web GIS), currently available at different spatial and temporal scales and with the support of graphical external functionalities. This review was necessary because, with the advent of weather digital networks, many datasets are indeed archived and publicly open, but often only short records (latest years) are available at subre- gional and local scales. Two of the most prominent reanalysis data sources that may serve the purposes of time–spatial modeling are the National Oceanic and Atmospheric Administration (NOAA)/National Climatic Data Center and National Aeronautics and Space Administration (NASA) Goddard Earth Sciences Data and Information Services Center. Review of online meteorological resources that are basic to prospective studies of climate variability is provided in the next section of this paper, together with a description of the study area. Section 3 assesses tem- poral variability of extreme rainfall observations (from 1801 to 2009) to reveal whether climate change may have implications on the extreme rainfall hazard. The case of Sicily (Italy), the largest island in the Mediterranean Sea and centrally situated in the basin, was examined with more detail (section 3.2) because its lands have been the theatre of storms of exceptional intensities in the autumn of 2009. Section 4 discusses possible implications between extreme rain events and their ecological impacts and concludes this paper.

2. Study area and data 2.1. Description of the study area The study area, hereafter called the Mediterranean central area (MCA), is lo- cated between 308–458N, 88–208E (Figure 2, large square box). The area of interest is approximately centered over the Tyrrhenian Sea, between the western coast of the Italian mainland and the northern coast of Sicily. Westward, the area is de- limited by southeastern France and northeastern Algeria. The Alps mountain chain and the Ligurian Sea are relevant geographic features to the north. The MCA extends eastward over two small basins (the Adriatic Sea and the Aegean Sea) up to the Balkans and includes former Yugoslavia, Albania, and a large part of Greece.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 5

Figure 2. Geographical setting and morphology of the Mediterranean basin. MCA (large square box) and Sicily (small box) are highlighted (base map from the International Water Management Institute; available online at http:// www.iwmi.cgiar.org).

The Strait of Sicily is a bridge between Tunisia and northern Libya (delimiting the southern margin of the MCA) and southern Europe, and it also subdivides the Mediterranean Sea into a western and an eastern basin. Water masses are indeed exchanged through both the Strait of Gibraltar (between Spain and Morocco) and the Strait of Sicily by eastward and westward flows. Water masses originating in the east of the basin flow westward and penetrate the Adriatic Sea and the Ionian Sea; on approaching the Strait of Sicily, part of the water mass is recirculated back within the eastern basin, whereas the remainder continues to enter the western Mediterranean basin (Rohling et al. 2009). Two broad areas of the Mediterranean region were left out of this study for their distinct climate patterns (Lionello et al. 2006). In a western area, including the Iberian Peninsula and Morocco (with the Atlas and Rif Mountains) and forming part of the eastern Atlantic–European region, orographic cyclogenesis is often triggered by the passage of Atlantic cyclones. Instead, areas of intense cyclogenic activity occur in the eastern Mediterranean region (including parts of Greece and Turkey) around Cyprus and the Middle East, to a large extent controlled by a large- scale Atlantic–European region. In the MCA, cyclogenesis is most frequent over the Gulf of Genoa and the rest of Ligurian Sea, but the Aegean Sea is a major center for cyclogenesis as well (e.g., Trigo et al. 1999). The majority of Genoan depressions tracks southeastward down the coast of Italy and then eastward or northeastward across the Aegean Sea. The variegated morphology of the MCA (basins and gulfs, mountainous groups, and peninsulas of various sizes) has important consequences on both sea and atmospheric circulations, which determine a nonuniform distribution of weather types (Lionello et al. 2006) and a large spectrum of associated hydrogeomorphological events (Petrucci and Polemio 2003; Sivakumer 2005). Three principal impact ways can be

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 6 described for precipitation (Diodato 2004): 1) very large-scale and low-intensity continuous rainstorms, 2) large frontal heavy thunderstorms with associated large erosive phenomena, and 3) short-period and localized convective rain showers leadingtoflashfloodsandmudflows.

2.2. Generation of precipitation data The rainfall datasets used in this study were derived from official records of data collected from a variety of national and international agencies for the Mediterra- nean region. The data used are mostly from the Royal Meteorological Institute (KNMI) Web site (available online at http://climexp.knmi.nl), which is a gateway to a number of online global regional databases and supplies the Web GIS KNMI Climate Explorer (van Oldenborgh and Burgers 2005), which is applied in this study to arrange for reanalysis fields and climate change experiments. We have accessed data archived at the NOAA/National Paleoclimatology Program (available online at http://www.ncdc.noaa.gov/paleo), the NOAA/Earth System Research Lab- oratory (ESRL) Physical Sciences Division Web site (Boulder, Colorado; available online at http://www.esrl.noaa.gov/psd), and the NASA Goddard Earth Sciences Data and Information Services Center (available online at http://disc2.nascom. nasa.gov). NOAA provides consistent long-term weather observations and has the most complete and updated global database of online weather records of daily res- olution (grid resolution of about 100 km). Paleoclimatology data with monthly res- olution come from natural sources (e.g., tree rings, ice cores, corals, and ocean and lake sediments) and extend the archives back to hundreds of years (e.g., Pauling et al. 2006) with grid resolutions of about 50 km. Specifically, reanalysis data are provided by the NOAA/ESRL Physical Sciences Division Web site. Since 1998, NASA archives and distributes high-resolution (;25 km) gridded precipitation data from the Tropical Rainfall Measuring Mission (TRMM) platform, providing satellite-based information on the intensity and dis- tribution of the rain type on the storm depth and on the height at which the snow melts into rain (e.g., Huffman et al. 2007). For particular and localized years or events that occurred in Italy, suitable datasets are those derived from regional monitoring centers. In this study, we integrated data detected by the Sicilian Agro- meteorological Information Service (SIAS; available online at http://www.sias. regione.sicilia.it) to refer precipitation data relative to the island of Sicily.

3. Extreme rainfall changes By comparing daily satellite observations with model data during a 20-yr period, Richard Allan and his collaborators at Reading University recently confirmed that heavy precipitation events are more frequent during warm, moist periods and less frequent during cold, dry periods (press release available online at http://www. reading.ac.uk/about/newsandevents/releases/PR16410.asp). To verify if the above hypothesis could be valid for a time long enough to include the period between the last part of the Little Ice Age (period of cooling confined to approximately the sixteenth century to the mid-nineteenth century; Lamb 1977) and the current phase of warming, we have compared some particularly wet decades that have affected the Mediterranean central basin since the early nineteenth century.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 7

Figure 3. (top) Gumbel plot (GEV based) of the autumn precipitation vs return period over MCA in (a) 1801–40 (cold period), (b) 1881–1920 (transitional period), and (c) 1970–2009 (latest and warmer period). The interpolating curve is reported with 95% confidence intervals for the most uncertain return pe- riods (latest interpolation and extrapolated periods). (bottom) The re- spective 98th percentile maps are compared [dataset sources: Pauling et al. 2006; KNMI (available online at http://www.knmi.nl/about_knmi), with updating until 2009].

3.1. Mediterranean central area The work of Casty et al. (Casty et al. 2007) on the European pattern climatology (1766–2000), which has been the basis for the development of this study, indicates no clear trend for precipitation totals in central Mediterranean when aggregated at annual and seasonal scales. However, from reviewing and updating that study for autumn (September–November) precipitation, we turned to a different elaboration of data, capturing the rainfall amounts via generalized extreme value (GEV) dis- tributions as Gumbel plots (e.g., Coles 2001) and comparing three time slices of about four decades each (Figures 3a–c), reflecting a cold period (1801–40), a transitional period (1881–1920), and a warmer period (1970–2009). From these plots, it is evident a gradual exacerbation of the extreme rainfalls over fixed return periods during the last two centuries. This trend was also em- phasized by the 98th percentile–based maps (Figures 3a–c, top). Although the above illustration delineates three distinct climate periods, most of the greater hydrological impacts have occurred in recent years, which is known to be the result of unusual deluges occurred with extraordinary intensity. A more enlarged view of the Mediterranean central area over the decade 2000–09 shows that autumn seasons are in fact prone to intensified precipitation rates with positive anomalies over many zones (Figure 4a). This agrees with the results of Cislaghi et al. (Cislaghi et al. 2005), who found that the average rain rate significantly increased in the 1800–2000 period in Italy, with associated shorter durations of rain

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 8

Figure 4. Rainfall rate anomalies for the decade 2000–09 compared to the period 1951–99 over the MCA (from NCEP reanalysis and NOAA/ESRL): (a) au- tumn and (b) September. For both autumn and September, average amounts were calculated dividing total rainfall by the number of days with measurable rainfall. episodes and with an evident effect on rainfall extremes. This is reflected in the surprising result of the extent to which temperature rise widens precipitation ex- tremes. According to Maheras et al. (Maheras et al. 1999), the widespread distri- bution of high temperatures over the Mediterranean tends to be associated with a negative sea level pressure anomaly in the Mediterranean basin that produces southwesterly flow from the ocean into the western Mediterranean, which may penetrate eastward or allow warmer air from North Africa to influence eastern parts of the basin. Anomalous warm conditions in parts of the Mediterranean and central Europe were also recognized by Luterbacher et al. (Luterbacher et al. 2007) as related to advection of warm air masses from the eastern subtropical Atlantic as well as strong anticyclonic conditions over large parts of the continent. Increased temperatures may result in nearly unchanged (but decreasing) precipitation totals but with a shift toward heavy, intense rainfalls (Alpert et al. 2002). In this context, it is relevant to learn how the past warming affected the changes of precipitation extremes (Klein Tank and Ko¨nnen 2003). The major increase of precipitation rates was observed, especially affecting the month of September over peninsular Italy, where the anomaly was about 11.5 mm day21 with an elliptical core reaching the island of Sicily (Figure 4b).

3.2. The case of Sicily As noted in Figure 4, the rise of rainfall is not limited to Italy, but also includes the island of Sicily. A focus on this island was meant to display the influence exerted by the power of rainfall on the climate regime (Figure 5). The rainfall amount

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 9

Figure 5. (a) Spatial pattern of rainfall anomalies in the months of September of the decade 1998–2008 compared to the respective monthly climatology 1951–99 (arranged from the TRMM earth science data; see also Acker and Leptoukh 2007); (b) rainfall amount that occurred on 16 Sep 2009 (;2-km grid, from Meteosat Second Generation rainfall estimation; available on- line at http://www.cespevi.it/previs2.htm). anomalies in the recent decade (1998–2008), compared to the climatological period 1951–99 (Figure 5a), show that only the central part of Sicily (includes the Etna volcano and the cities of Enna and Caltanissetta) is apparently not affected by rainfall increase. The positive anomalies depicted in Figure 5a are often the result of rainfalls concentrated in very few days because those occurred in September 2009, when the monthly rain amount registered at the Palermo Astronomical Observatory (388N, 138E), 230 mm, exceeded the previous monthly record of 190 mm dating back to 1820 (V. Iuliano 2009, personal communication; http://www.astropa.unipa.it). To give an idea of the magnitude of the extreme events that characterized September 2009, we also illustrate the rainfall occurred on the 16th day of the same month, when a storm depth of about 100 mm was recorded over and around Palermo district and south of the town of Messina (388N, 158E), with a main storm core of 200 mm, northeast of Palermo (Figure 5b). The extraordinary nature of the situation was the result of repeated deluges that afflicted several lands of Sicily in September 2009, combined with those occurred in the months of September of past years. For instance, deluges were recorded on 17 September 2003 at some stations of Syracuse province, reaching about 600 mm of rain, 398 mm of which fell in only 6 h. Phenomena patterns at these localities show that subgrid-scale convection and intensification are dominating the rain- producing mechanisms (Mazzarella 1999; Du¨nkeloh and Jacobeit 2003) and are shared with several rain showers releasing in few hours as much energy as equal to or higher than the annual amount. This is in agreement with the results of Bonaccorso et al. (Bonaccorso et al. 2005), by which subgrid-scale convection and intensification phenomena indicate for Sicily an increasing trend toward shorter rain durations (about 1 h) during the period 1927–2004. In Figure 6a, time evolution is shown of daily precipitation rates averaged over Sicily during the September months from 1948 to 2009. The graph shows changes ramping only during the last decade according to a climate shift. This indicates an

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 10

Figure 6. Temporal pattern of precipitation rate over Sicily: (a) September trend over 1948–2009 and (b) monthly regimen during the decade 2000–09 (gray line), compared to the climatological period 1948–99 [bold line; arranged by NCEP reanalysis and NOAA-ESRL (available online at http://www.cdc. noaa.gov/cgi-bin/PublicData/getpage.pl); Kalnay et al. 1996].

accelerated power of rainfall with severe and important impacts on ecosystems because tilled soils in the beginning of autumn are expanding with increased vulnerability of lands. Based on Figure 6b, it can be argued that the intensifying precipitation affecting the Sicilian lands could not be limited to the month of September only, but could also affect the months of November and December (arrows in Figure 6b). It is also noteworthy to mention the violent storms that struck Sicily between 1 and 2 October 2009. On this occasion, the province of Messina was the most affected, with Fiumedinisi being the rainiest site with 159 mm on 1 October (SIAS data). To interpret the regional patterns of rainfall–climate interaction over Sicily in the context of the Mediterranean area, the evolution of the standardized precipitation index (SPI; McKee et al. 1993) was also examined (Figure 7), as arranged from the (Ispra, Italy; available online at http://www.jrc.it). SPI is based on the probability of recording a given amount of precipitation occurring over a given prior time period (which may vary from 1 to 36 months), and the probabilities are standardized so that a value of zero indicates the median precipitation amount. The index is negative for drought and positive for wet conditions, and becomes more negative or positive as the dry or wet conditions become more severe. For the month of September, the yearly evolution of SPI values calculated at Palermo over the prior 3 months (Figure 7a) indicates mani- festations of climate change resulting from a gradual (although irregular) shift from dry to wet conditions. The year 2009 is remarkable after a previous peak observed in 2000. For September 2009, the map of SPI over the Mediterranean area (Figure 7b) indicates large areas of dryness at northern and central latitudes, whereas signals of wetter climate are apparent at most sites of the southern Mediterranean, including Sicily. This evidence indicates an unusual reversal of precipitation patterns be- tween north and south in response to increased warmth.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 11

Figure 7. SPI: (a) 3-month-based SPI evolution (1973–2009) for September at Palermo and (b) map of SPI in September 2009 over Mediterranean area (arranged from the European Commission Joint Research Centre Institute for Environment and Sustainability; available online at http://desert. jrc.ec.europa.eu).

4. Summary and environmental implications Earth ecosystems are not static and respond to environmental changes on many scales and from a variety of causes, particularly climatic and anthropogenic, which for the Mediterranean area can be among the most severe worldwide (Allen 2003). This paper presents an approach, based on the investigation of past and current extreme rain events, to deepen the consequences of this underappreciated aspect of climate change over the Mediterranean central area. Reviewing online precipitation resources was basic to this research in view of prospective studies of climate vari- ability. Historical datasets available for the Mediterranean area were not always updated until today. The datasets used in this study to explain over 200 years of precipitation range from rain gauges, reanalysis data, and satellite inferences. The National Centers for Environmental Prediction (NCEP) reanalysis represents an updated source providing relatively coarse spatial data (about 18318), not allowing in-depth studies at a subregional level. To overcome such a limitation, we based some analyses on the TRMM platform, a high-resolution source also updated in time. Where possible, some comparisons between these two main sources showed substantially similar results (data not shown) that support using the two sources together in this climatological study. Other sources referenced in our study (e.g., regional monitoring centers) are minor datasets used in support to the two main sources at either early or late assessment. Without limiting our study to investigate important rain extremes recently occurred (e.g., September 2009), attempts were made to place them in the context of the climate history of the MCA. Precipitation is, in fact, a key link to the global water cycle, and a proper understanding of its temporal and spatial character will have broad implications in ongoing climate diagnostics and predictions, global water and energy cycles, analysis and modeling, weather forecasting, freshwater resource management, and land–atmosphere–ocean interface processes (Shepherd and Burian 2003). Our study relies on the importance of concepts such as sequence and position, as well as magnitude and frequency of storm events, stemming out of Wolman and

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 12 Miller (Wolman and Miller 1960) and employed by Richards (Richards 1999). Because climatic events (such as droughts and storms) are sometimes clustered into longer-term groups, their impacts may not vary simply in relation to their size, thus complicating any simple magnitude and frequency relationship. Climatic vari- ability is recognized as one key cause of these clusterings; thus, identification of modes of climatic variability occurring at the decadal scale going back into the Holocene may help testing the effectiveness of different rain events within real event histories and within the broad sweep of environmental and ecological history (Viles and Goudie 2003). If the Little Ice Age is identified as the stormiest period of the last four to five centuries in the Mediterranean, important historical climate fluctuations call for high-time-resolution studies to assess seasonal changes of rain extremeness and their relation to soil conservation (e.g., Diodato et al. 2008). Based on the findings of the present study, focused on September months, alternating centuries of limited and highly hazardous rainfalls spaced out by a transitional century seem to have occurred in the recent past (three centuries) at the Mediterranean central area as a response to the changing temperatures. In a previous study on the environmental implications of erosive rainfall across the Mediterranean (Diodato and Bellocchi 2009a), large erosive rainfalls were observed to occur, especially in summer (continental areas) and autumn (along-coast and near-coast reliefs). It was also shown, from detailed explorations of temporal patterns of erosivity (Diodato and Bellocchi 2009a; Diodato and Bellocchi 2009b), that increasing the number of extreme events in autumn does not cause seasonal rainfall totals to deviate from the historical range of climate variation, but it tends to generate more disproportion between Mediterranean dry and wet periods, which could bring soil loss to higher rates (Diodato and Bellocchi 2008). If this rainfall regime were to continue, it could result in ever-increasing exacerbated erosive hazard affecting Mediterranean countries in an erratic way. The magnitude and even the direction of climate change may differ among regions, so it appears reasonable to focus on a single region in attempting to assess probable effects of climate change when that region has some intrinsic meaning. Based on its central position in the Mediterranean, we have identified the island of Sicily as a special case for this study. Placed at the crossroads for Mediterranean flows, Sicily can be considered representative of a geographic scenario of regional climate change associated with warming, reflected in the Mediterranean cyclo- genicity, and that would have contributed to the exceptional rainfall rates observed in recent times. During September, practically the whole island has experienced remarkable positive anomalies in the first decade of the twenty-first century in comparison to the monthly climatology, with a trend toward increasing daily rates. Also, similar evidence in other months indicates that precipitation patterns in this region are actually on the verge of important changes. This is also indicated by the information conveyed by the evolution of the SPI at individual Sicilian places (e.g., Palermo) and over the Mediterranean region. Following all of this evidence, we argue that the assessment of current and future management systems should not only be based on the average rainfall for a period, but it should also include the hazard of extreme precipitation events, which likely cause accelerated urban and rural land degradation. This must be considered to- gether with the findings that, although rainfall amounts are not always increasing,

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 13 erratic spatial and temporal storm patterns in some seasons or months drive the erosive power of rain to increase its hazard. According to D’Asaro et al. (D’Asaro et al. 2007), investigations are needed, especially in terms of statistical analysis of long pluviometric series, rain intensity, and erosivity occurring at decadal and in- terdecadal scales, which mark substantial climatic fluctuations and changes during last centuries. More generally, identification of links between nonlinear climatic, ecological, and geomorphological systems should lead to search for increased un- derstanding of their mutual interactions and behavior and not to use parts of them in simple cause–effect relationships. Quantitative reconstruction of the impacts of climate change on organisms and ecosystems represents an important and chal- lenging line of enquiry. Improved conceptualization of the impacts of decadal-scale climatic variability can help interpreting the projections of the ecological impacts of future climate change. The identification of climate boundaries (areas with homo- geneous climate) of rain scenarios at regional scales (as the MCA in this study) is important for understanding the climatic features and their potential impact on ecosystems. From all the evidence presented in this paper and in the recent literature, it becomes apparent that variability of events in time will be an important part of the environmental future of the MCA. Rain intensity and erosive power on decadal and interdecadal scales remain to mark substantial climatic fluctuations and changes during the last centuries, although the perception of disasters by public domain is not only determined by objective data but also from the relevance that media give to it (Lastoria et al. 2006). This perception is central to the process of the resource allocation for the risk mitigation; in general, it somehow influences the type and entity of measures taken by authorities and people involved in flood risk man- agement. An increased understanding of the properties of decadal and interdecadal climatic variability also puts us in a better position to understand the possibilities and limitations of extrapolating future climate events and impacts. The results gained from the present study encourage further studies on climatic change by simulation models and statistical approaches to incorporate the characteristics offered by finer time and spatial scales.

Acknowledgments. This work is part of the Hydrological Cycle in the Mediterranean Experiment (HYMEX; available online at http://www.hymex.org) aiming at hydromete- orological modeling for the whole Mediterranean basin.

References Acker, J. G., and G. Leptoukh, 2007: Online analysis enhances use of NASA Earth science data. Eos, Trans. Amer. Geophys. Union, 88, doi:10.1029/2007EO020003. Alca´ntara-Ayala, I., 2002: Geomorphology, natural hazards, vulnerability and prevention on natural disaster in developing countries. Geomorphology, 47, 107–124. Allen, H. D., 2003: Response of past and present Mediterranean ecosystems to environmental change. Prog. Phys. Geogr., 27, 359–377. Alpert, P., and Coauthors, 2002: The paradoxical increase of Mediterranean extreme daily rainfall in spite of decrease in total values. Geophys. Res. Lett., 29, 1536, doi:10.1029/2001GL013554. Barroca, B., P. Bernardara, J. M. Mouchel, and G. Hubert, 2006: Indicators for identification of urban flooding vulnerability. Nat. Hazards Earth Syst. Sci., 6, 553–561.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 14 Bonaccorso, B., A. Cancelliere, and G. Rossi, 2005: Detecting trends of extreme rainfall series in Sicily. Adv. Geosci., 2, 7–11. Camarasa Belmonte, A., and F. Segura Beltra´n, 2001: Flood events in Mediterranean ephemeral streams (ramblas) in Valencia region, Spain. Catena, 45, 229–249. Casty, C., C. C. Raible, T. F. Stocker, H. Wanner, and J. Luterbacher, 2007: A European pattern climatology 1766–2000. Climate Dyn., 29, 791–805. Cavazos, T., and D. Rivas, 2004: Variability of extreme precipitation events in Tijuana, Mexico. Climate Res., 25, 229–243. Cislaghi, M., T. C. De Michele, A. Ghezzi, and R. Rosso, 2005: Statistical assessment of trends and oscillations in rainfall dynamics: Analysis of long daily Italian series. Atmos. Res., 77, 188–202. Clarke, M., and H. Rendell, 2007: Climate, extreme events and land degradation. Climate and Land Degradation, M. V. K. Sivakumar and N. Ndiang’ui, Eds., Springer, 137–152. Coles, S., 2001: An Introduction to Statistical Modeling of Extreme Values. Springer, 224 pp. Curtis, S., A. Salahuddin, R. F. Adler, G. J. Huffman, G. Gu, and Y. Hong, 2007: Precipitation extremes estimated by GPCP and TRMM: ENSO relationships. J. Hydrometeor., 8, 678–689. D’Asaro, F., L. D’Agostino, and V. Bagarello, 2007: Assessing changes in rainfall erosivity in Sicily during the twentieth century. Hydrol. Processes, 21, 2862–2871. Diodato, N., 2004: Local models for rainstorm-induced hazard analysis on Mediterranean river- torrential geomorphological systems. Nat. Hazards Earth Syst. Sci., 4, 389–397. ——, and G. Bellocchi, 2008: Drought stress patterns in Italy using agro-climatic indicators. Climate Res., 36, 53–63. ——, and ——, 2009a: Environmental implications of erosive rainfall across Mediterranean. Envi- ronmental Impact Assessments, G. T. Halley and Y. T. Fridian Eds., Nova Science Publishers, 225–253. ——, and ——, 2009b: Assessing and modelling changes in rainfall erosivity at different climate scales. Earth Surf. Processes Landforms, 34, 969–980. ——, M. Ceccarelli, and G. Bellocchi, 2008: Decadal and century-long changes in the recon- struction of erosive rainfall anomalies in a Mediterranean fluvial basin. Earth Surf. Processes Landforms, 33, 2078–2093. D’Odorico, P., J. Yoo, and T. M. Over, 2001: An assessment of ENSO-induced patterns of rainfall erosivity in the southwestern United States. J. Climate, 14, 4230–4242. Du¨nkeloh, A., and J. Jacobeit, 2003: Circulation dynamics of Mediterranean precipitation vari- ability 1948-1998. Int. J. Climatol., 23, 1843–1866. Easterling, D. R., J. L. Evans, P. Ya Groisman, T. R. Karl, K. E. Kunkel, and P. Ambenje, 2000: Observed variability and trends in extreme climate events: A brief review. Bull. Amer. Meteor. Soc., 81, 417–425. Garcı´a-Oliva, F., J. M. Maass, and L. Galicia, 1995: Rainstorm analysis and rainfall erosivity of a seasonal tropical region with a strong cyclonic influence on the Pacific coast of Mexico. J. Appl. Meteor., 34, 2491–2498. Groisman, P. Ya., R. W. Knight, D. R. Easterling, T. R. Karl, G. C. Hegerl, and V. N. Razuvaev, 2005: Trends in intense precipitation in the climate record. J. Climate, 18, 1326–1350. Gruntfest, E., and J. Handmer, 2001: Coping with Flash Floods. NATO Science Series, Vol. 77, Springer, 322 pp. Huffman, G. J., and Coauthors, 2007: The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydro- meteor., 8, 38–55. Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471. Klein Tank, A. M. G., and G. P. Ko¨nnen, 2003: Trends in indices of daily temperature and pre- cipitation extremes in Europe, 1946–99. J. Climate, 16, 3665–3680. Knapp, A. K., and Coauthors, 2008: Consequences of more extreme precipitation regimes for terrestrial ecosystems. Bioscience, 58, 811–821.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 15 Kunkel, K. E., R. A. Pielke, and S. A. Changnon, 1999: Temporal fluctuations in weather and climate extremes that cause economic and human health impacts: A review. Bull. Amer. Meteor. Soc., 80, 1077–1098. Lamb, H. H., 1977: Climatic History and the Future. Vol. 2, Climate: Present, Past and Future, Methuen, 603 pp. Lastoria, B., M. R. Simonetti, M. Casaioli, S. Mariani, and G. Monacelli, 2006: Socio-economic impacts of major floods in Italy from 1951 to 2003. Adv. Geosci., 7, 223–229. Lionello, P., and Coauthors, 2006: Cyclones in the Mediterranean region: Climatology and effects on the environment. Mediterranean Climate Variability, P. Lionello, P. Malanotte-Rizzoli, and R. Boscolo, Eds., Developments in Earth and Environmental Sciences, Vol. 4, Elsevier, 325–372. Luterbacher, J., M. A. Liniger, A. Menzel, N. Estrella, P. M. Della-Marta, C. Pfister, T. Rutishauser, and E. Xoplaki, 2007: Exceptional European warmth of autumn 2006 and winter 2007: Historical context, the underlying dynamics, and its phenological impacts. Geophys. Res. Lett., 34, L12704, doi:10.1029/2007GL029951. Maheras, P., E. Xoplaki, T. Davies, J. Martı´n-Vide, M. Barriendos, and M. J. Alcoforado, 1999: Warm and cold monthly anomalies across the Mediterranean basin and their relationship with circulation; 1860-1990. Int. J. Climatol., 19, 1697–1715. Martina, M. L. V., E. Todini, and A. Libralon, 2006: A Bayesian decision approach to rainfall thresholds based flood warning. Hydrol. Earth Syst. Sci., 2, 413–426. Mazzarella, A., 1999: Multifractal dynamic rainfall processes in Italy. Theor. Appl. Climatol., 63, 73–78. McKee, T. B., N. J. Doesken, and J. Kleist, 1993: The relationship of drought frequency and duration to time scales. Proc. Eighth Conf. on Applied Climatology, Anaheim, CA, Amer. Meteor. Soc., 5.1. Papathoma, M., and D. Dominey-Howes, 2003: Tsunami vulnerability assessment and its impli- cations for coastal hazard analysis and disaster management planning, Gulf of Corinth, Greece. Nat. Hazards Earth Syst. Sci., 3, 733–747. Pauling, A., J. Luterbacher, C. Casty, and H. Wanner, 2006: Five hundred years of gridded high- resolution precipitation reconstructions over Europe and the connection to large-scale cir- culation. Climate Dyn., 26, 387–405. Peterson, T. C., and Coauthors, 2002: Recent changes in climate extremes in the Caribbean region. J. Geophys. Res., 107, 4601, doi:10.1029/2002JD002251. Petrucci, O., and M. Polemio, 2003: The use of historical data for the characterisation of multiple damaging hydrogeological events. Nat. Hazards Earth Syst. Sci., 3, 17–30. Pitman, A. J., and R. J. Stouffer, 2006: Abrupt change in climate and climate models. Hydrol. Earth Syst. Sci., 10, 903–912. Poesen, J. W. A., and J. M. Hooke, 1997: Erosion, flooding and channel management in Medi- terranean environments of southern Europe. Prog. Phys. Geogr., 21, 157–199. Ramos, M. C., and M. Mulligan, 2003: Impacts of climate variability and extreme events on soil hydrological processes. Geophys. Res., 5, 92–115. Richard, A. K., 2007: Global warming is changing the world. Science, 316, 188–190. Richards, K., 1999: The magnitude-frequency concept in fluvial geomorphology: A component of a degenerating research programme? Z. Geomorphol., 155, 1–18. Rohling, E. J., R. Abu-Zied, C. S. L. Casford, A. Hayes, and B. A. A. Hoogakker, 2009: The Mediterranean Sea: Present and past. Physical Geography of the Mediterranean Basin, J. C. Woodward, Ed., Oxford University Press, 33–68. Sa´nchez, J., E. C. Gallardo, M. Gaertner, A. Arribas, and M. Castro, 2004: Future climate extreme events in the Mediterranean simulated by a regional climate model: A first approach. Global Planet. Change, 44, 163–180. Shepherd, J. M., and S. J. Burian, 2003: Detection of urban-induced rainfall anomalies in a major coastal city. Earth Interactions, 7. [Available online at http://EarthInteractions.org.]

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC Earth Interactions d Volume 14 (2010) d Paper No. 5 d Page 16 Sivakumer, M. V. K., 2005: Impacts of natural disasters in agriculture, rangeland and forestry: An overview. Natural Disasters and Extreme Events in Agriculture, M. V. K. Sivakumer, R. P. Motha, and H. P. Das, Eds., Springer-Verlag, 1–22. Thornes, J. B., and I. Alca´ntara-Ayala, 1998: Modelling mass failure in a Mediterranean mountain environment: Climatic, geological, topographical and erosional controls. Geomorphology, 24, 87–100. Trigo, I. F., T. D. Davies, and G. R. Bigg, 1999: Objective climatology in the Mediterranean region. J. Climate, 12, 1685–1696. van Leeuwen, W. J. D., and G. Sammons, 2003: Seasonal land degradation risk assessment for Arizona. Proc. 30th Int. Symp. on Remote Sensing of Environment, ISPRS, 378–381. [Available online at http://wildfire.arid.arizona.edu/methods.htm.]. van Oldenborgh, G. J., and G. Burgers, 2005: Searching for decadal variations in ENSO precipi- tation teleconnections. Geophys. Res. Lett., 32, L15701, doi:10.1029/2005GL023110. van Rompaey, A., P. Bazzoffi, R. J. A. Jones, and L. Montanarella, 2005: Modeling sediment yields in Italian catchments. Geomorphology, 65, 157–169. Viles, H. A., and A. S. Goudie, 2003: Interannual, decadal and multidecadal scale climatic vari- ability and geomorphology. Earth Sci. Rev., 61, 105–131. Wei, W., L. Chen, and B. Fu, 2009: Effects of rainfall change on water erosion processes in terrestrial ecosystems: A review. Prog. Phys. Geogr., 33, 1–12. Willmott, C. J., and D. R. Legates, 1991: Rising estimates of terrestrial and global precipitation. Climate Res., 1, 179–186. Wolman, M. G., and J. P. Miller, 1960: Magnitude and frequency of forces in geomorphic processes. J. Geol., 68, 54–74. Zehe, E., and M. Sivapalan, 2009: Threshold behaviour in hydrological systems as (human) geo- ecosystems: Manifestations, controls, implications. Hydrol. Earth Syst. Sci., 13, 1273–1297. Zolina, O., C. Simmer, K. Belyaev, A. Kapala, and S. Gulev, 2009: Improving estimates of heavy and extreme precipitation using daily records from European rain gauges. J. Hydrometeor., 10, 701–716.

Earth Interactions is published jointly by the American Meteorological Society, the American Geophysical Union, and the Association of American Geographers. Permission to use figures, tables, and brief excerpts from this journal in scientific and educational works is hereby granted provided that the source is ac- knowledged. Any use of material in this journal that is determined to be ‘‘fair use’’ under Section 107 or that satisfies the conditions specified in Section 108 of the U.S. Copyright Law (17 USC, as revised by P.IL. 94- 553) does not require the publishers’ permission. For permission for any other from of copying, contact one of the copublishing societies.

Unauthenticated | Downloaded 09/25/21 07:57 AM UTC