Analysis of Soil Erosion Induced by Heavy Rainfall: A Case Study from the NE Abruzzo Hills Area in Central Italy

Soil erosion induced by heavy rainfall deeply affects landscape changes and human activities. It depends on rainfall distribution (e.g., intensity, duration, cumulative per event) and is controlled by the interactions between lithology, orography, hydrography, land use, and vegetation. The Abruzzo piedmont coastal hilly area has been affected by several heavy rainfall events in the last decades. In this work, we investigated three ~1-day heavy rainfall (>35 mm/h and 100–220 mm/day) events in 2007, 2011, and 2012 that occurred in the clayey hilly coastal NE Abruzzo area, analyzing cumulative rainfall, intensity, and duration while mapping triggered geomorphological effects (soil erosion and accumulation) and evaluating average erosion. The analysis provides contributions to a soil erosion assessment of clayey landscapes that characterizes the Adriatic hilly area, with an estimation of rainfall-triggering thresholds for heavy soil erosion and a comparison of erosion in single events with rates known in the Mediterranean area. The triggering threshold for heavy soil erosion shows an expected value of ~100–110 mm. The estimated average soil erosion is from moderate to high (0.08–3.08 cm in ~1-day heavy rainfall events) and shows a good correlation with cumulative rainfall and a poor correlation with peak rainfall intensity. This work outlines the strong impact of soil erosion on the landscape changes in the Abruzzo and Adriatic hilly areas.


Introduction
Rainfall and heavy rainfall events are among the most important factors inducing soil erosion and landslides, particularly in Mediterranean countries [1][2][3][4]. Their impact on slope instabilities is mainly related to the interaction between lithology, orography, hydrography, land use, and vegetation [5][6][7][8][9][10]. The relationship between rainfall intensity/duration and the geomorphological effects on the landscape, especially in terms of heavy soil erosion, is an open issue [1][2][3]5]. This is particularly true in the framework of changing climate/land use and increasing flash flooding induced by extremely heavy rainfall events [8,[11][12][13]. These events have a strong impact on landscape changes and deeply affect human activities, inducing large losses and even casualties [14][15][16][17]. For these reasons, they are studied with different approaches and data in Mediterranean and arid environments (e.g., rain gauge data, satellite rainfall data, time series statistical analysis, remote sensing, climate modeling, empirical

Regional Setting
The NE Abruzzo region is located in the central eastern part of the Italian peninsula along the piedmont area of the Central Apennines and the hilly coastal area (Figure 1). It includes the lower part of the main SW-NE to W-E fluvial valleys (the Tronto, Vibrata, Salinello, Vomano, and Saline rivers), and the small tributary catchment of the main rivers and those incising the coastal slopes.
The elevation ranges from 0 to ~600 m along the main drainage divides between the main rivers, and the terrain gradient (derived from a 10-m-cell digital elevation model, DEM) ranges from 0° along the river and coastal plains to >40°, with local vertical scarps along the slopes. There are four lithologic bedrock complex, or groups of rock units, in the area, each comprising different marine sedimentary rock types varying in strength from hard to weak and soft rocks: (1) sandstone and claystone (in the inner hilly areas); (2) clay (in the lower slopes of the coastal hilly areas); (3) sand and weak sandstone (in the upper slopes of the hilly coastal areas); and (4) conglomerate and consistent gravel (at the top of the main coastal hills). In the recent official geological map of Italy [44,45], these complexes are pertaining to Neogene turbiditic arenaceous and pelitic rocks, and to Plio-Pleistocene marine clay-sand-sandstone-conglomerate rocks. The bedrock units are covered with Middle Pleistocene-Holocene clastic continental deposits (slope, landslides, fluvial deposits) at the bases of the slopes and in the main valleys ( Figure 2). The hillslopes and coastal slopes are mantled by clay-and sand-rich eluvial and colluvial cover, for which thickness ranges from < 1 m to >10 m depending on morphology and lithology (thicker in clay bedrock and within minor valleys). Finally, the superficial soil cover reflects the lithological types, with thickness ranging from a few centimeters to >1 m. The piedmont and coastal hilly areas, as well as the river and coastal plains, are mostly characterized by luvisoils, vertisoils, and cambisoils ( Figure 3 [46]) developed on clayey and sandy parent materials [47], which are typical of terrains affected by high erosion rates [48]. The structural setting is defined in the inner part by thrust reliefs, which are affected by regional NNW-SSE Pliocene thrusts and minor high-angle normal faults and, in the outer part, by a large slightly NE-dipping homocline [44,45,[49][50][51][52][53].
The piedmont coastal area features a hilly landscape carved by cataclinal valleys (SW-NE) on arenaceous-pelitic thrusted and faulted successions and on a gently NE-dipping homocline on clay, sand, and conglomerate deposits. This landscape is incised by the main river valleys (featuring up to 4-km-wide fluvial plains) and by minor catchments flowing toward the main valleys and the coastal There are four lithologic bedrock complex, or groups of rock units, in the area, each comprising different marine sedimentary rock types varying in strength from hard to weak and soft rocks: (1) sandstone and claystone (in the inner hilly areas); (2) clay (in the lower slopes of the coastal hilly areas); (3) sand and weak sandstone (in the upper slopes of the hilly coastal areas); and (4) conglomerate and consistent gravel (at the top of the main coastal hills). In the recent official geological map of Italy [44,45], these complexes are pertaining to Neogene turbiditic arenaceous and pelitic rocks, and to Plio-Pleistocene marine clay-sand-sandstone-conglomerate rocks. The bedrock units are covered with Middle Pleistocene-Holocene clastic continental deposits (slope, landslides, fluvial deposits) at the bases of the slopes and in the main valleys ( Figure 2). The hillslopes and coastal slopes are mantled by clay-and sand-rich eluvial and colluvial cover, for which thickness ranges from < 1 m to >10 m depending on morphology and lithology (thicker in clay bedrock and within minor valleys). Finally, the superficial soil cover reflects the lithological types, with thickness ranging from a few centimeters to >1 m. The piedmont and coastal hilly areas, as well as the river and coastal plains, are mostly characterized by luvisoils, vertisoils, and cambisoils ( Figure 3 [46]) developed on clayey and sandy parent materials [47], which are typical of terrains affected by high erosion rates [48]. The structural setting is defined in the inner part by thrust reliefs, which are affected by regional NNW-SSE Pliocene thrusts and minor high-angle normal faults and, in the outer part, by a large slightly NE-dipping homocline [44,45,[49][50][51][52][53].
The piedmont coastal area features a hilly landscape carved by cataclinal valleys (SW-NE) on arenaceous-pelitic thrusted and faulted successions and on a gently NE-dipping homocline on clay, sand, and conglomerate deposits. This landscape is incised by the main river valleys (featuring up to 4-km-wide fluvial plains) and by minor catchments flowing toward the main valleys and the coastal plain (featuring radial, trellis, and angular drainage patterns [49]). The coastal area is characterized by wide coastal slopes and a coastal plain up to >1 km wide [49,54]. The geomorphological processes are mainly fluvial, gravity-induced, and mass-wasting. These processes are frequently activated by the heavy rainfall events that affect the region. Fluvial processes affect the main rivers, alternating between channel incisions and flooding. The slope processes due to running water mostly affect the clayey and arenaceous-pelitic hills of piedmont and the coastal areas, generating badlands and minor landforms such as rills, gullies, and mudflows [41]. Mass-wasting processes have induced the formation of a huge number of landslides and mass movements, mostly affecting the hilly piedmont area and the chain area and, locally, the coastal area [54][55][56][57]. Moreover, the clayey hill areas are largely affected by both short-term and long-term heavy soil erosion processes [58]. The long-term average erosion rate is documented to be~0.1-14 mm/year, and the short-term erosion rate is 1-90 mm/year in this area and in the central Italy clayey hills [29,[39][40][41][42].
The hilly piedmont coastal area is characterized by a maritime Mediterranean climate [59]. The average annual precipitation is 600-800 mm/y, with occasional heavy rainfall (>100 mm/day and 30-40 mm/h) [33,34,[60][61][62]. In the last decades, this area was affected by flash flood events induced by heavy rainfall ranging from 60 to 80 mm in a few hours to >200 mm in one day (e.g.  plain (featuring radial, trellis, and angular drainage patterns [49]). The coastal area is characterized by wide coastal slopes and a coastal plain up to >1 km wide [49,54]. The geomorphological processes are mainly fluvial, gravity-induced, and mass-wasting. These processes are frequently activated by the heavy rainfall events that affect the region. Fluvial processes affect the main rivers, alternating between channel incisions and flooding. The slope processes due to running water mostly affect the clayey and arenaceous-pelitic hills of piedmont and the coastal areas, generating badlands and minor landforms such as rills, gullies, and mudflows [41]. Mass-wasting processes have induced the formation of a huge number of landslides and mass movements, mostly affecting the hilly piedmont area and the chain area and, locally, the coastal area [54][55][56][57]. Moreover, the clayey hill areas are largely affected by both short-term and long-term heavy soil erosion processes [58]. The long-term average erosion rate is documented to be ~0.1-14 mm/year, and the short-term erosion rate is 1-90 mm/year in this area and in the central Italy clayey hills [29,[39][40][41][42].
The hilly piedmont coastal area is characterized by a maritime Mediterranean climate [59]. The average annual precipitation is 600-800 mm/y, with occasional heavy rainfall (>100 mm/d and 30-40 mm/h) [33,34,[60][61][62]. In the last decades, this area was affected by flash flood events induced by heavy rainfall ranging from 60 to 80 mm in a few hours to >200 mm in one day (e.g.   [44,45,63]). The red boxes indicate locations of the study areas (grid in UTM WGS84 coordinates).

Figure 2.
Geologic scheme of the Abruzzo region (modified from [44,45,63]). The red boxes indicate locations of the study areas (grid in UTM WGS84 coordinates). In terms of land use, the hilly coastal area of NE Abruzzo is characterized by arable lands, and the number of high-quality plantations (mainly olive groves and vineyards) has increased in the last decades; minor urban areas and sparsely urbanized rural areas are also present. The fluvial plains show alternating cropland and industrial commercial areas. Finally, the coastal plain is mostly seamlessly urbanized with villages, tourist facilities, roads, railroad, ports, etc.
Water 2018, 10, x FOR PEER REVIEW 5 of 28 In terms of land use, the hilly coastal area of NE Abruzzo is characterized by arable lands, and the number of high-quality plantations (mainly olive groves and vineyards) has increased in the last decades; minor urban areas and sparsely urbanized rural areas are also present. The fluvial plains show alternating cropland and industrial commercial areas. Finally, the coastal plain is mostly seamlessly urbanized with villages, tourist facilities, roads, railroad, ports, etc.

Rainfall Event of 6-7 October 2007
The 2007 heavy rainfall event affected, for a short time (14-16 h), a local area in the northern Abruzzo (hilly and coastal Tortoreto area; Figure 2), including tributary catchments between the Salinello and Vibrata rivers and coastal catchment (total area of 43.7 km 2 ). River and coastal plains (slope gradient > 5°) cover 35% of the area, while 65% is hills (slope gradient 5-30°, with local scarps). The event occurred after two dry months and affected a very poorly vegetated landscape. The agricultural areas (arable land, vineyards, and olive groves; Figure 4a) were largely plowed on erodible clay bedrock with sands and conglomerates on the top of the hill, both largely covered with clayey eluvial and colluvial cover and landslides (Figure 4b).

Rainfall Event of 6-7 October 2007
The 2007 heavy rainfall event affected, for a short time (14-16 h), a local area in the northern Abruzzo (hilly and coastal Tortoreto area; Figure 2), including tributary catchments between the Salinello and Vibrata rivers and coastal catchment (total area of 43.7 km 2 ). River and coastal plains (slope gradient > 5 • ) cover 35% of the area, while 65% is hills (slope gradient 5-30 • , with local scarps). The event occurred after two dry months and affected a very poorly vegetated landscape. The agricultural areas (arable land, vineyards, and olive groves; Figure 4a) were largely plowed on erodible clay bedrock with sands and conglomerates on the top of the hill, both largely covered with clayey eluvial and colluvial cover and landslides (Figure 4b).

Rainfall Event of 1-2 March 2011
The 2011 heavy rainfall event affected the entire northeastern Abruzzo piedmont hilly area (hilly and coastal Teramo area between the Vomano, Tordino, Salinello, and Vibrata rivers; Figure 2) for a moderately short time (22-26 h), after~10 days of moderate rainfall (total 50 mm). After this event, two areas were investigated in detail because they suffered the heaviest rainfall and induced effects: the Salinello River valley and the Pineto-Atri hilly area ( Figure 2). The Salinello area includes the lower part of the Salinello River valley and the surrounding hillslopes (including part of the Tortoreto hilly area affected by the 2007 event). Over 50% of the area pertains to the river and the coastal plains; the remaining area incorporates the hilly relief up to~300 m a.s.l. with a slope gradient of 5-30 • (total area 39.7 km 2 ; Figure 5). With respect to the 2007 event, which partially affected the same area (lower Salinello hillslopes), the 2011 event occurred on a moderately vegetated landscape with agricultural areas (arable land, vineyards, and olive groves; Figure 5a) at an initial crop growth stage and during grass development. This event occurred again on clay-sand-conglomerate bedrock blanketed by a clayey eluvial and colluvial cover and landslides (Figure 5b).
The Pineto-Atri area consists of three coastal catchments, the intervening coastal slopes, and an up to <1 km wide coastal plain (total area 60.6 km 2 ; Figure 6). The elevation ranges from sea level tõ 450 m a.s.l. The coastal and river plain areas make up~15% of the area, while the hilly area (slope gradient 5-35 • with local scarps) makes up the remaining 85%. Similar to the others, the Pineto area is characterized by arable land and high-quality plantations in the hilly areas, but also large bare clay areas within widespread badlands affecting the south-facing clay slopes. The 2011 heavy rainfall event affected the entire northeastern Abruzzo piedmont hilly area (hilly and coastal Teramo area between the Vomano, Tordino, Salinello, and Vibrata rivers; Figure 2) for a moderately short time (22-26 h), after ~10 days of moderate rainfall (total 50 mm). After this event, two areas were investigated in detail because they suffered the heaviest rainfall and induced effects: the Salinello River valley and the Pineto-Atri hilly area ( Figure 2). The Salinello area includes the lower part of the Salinello River valley and the surrounding hillslopes (including part of the Tortoreto hilly area affected by the 2007 event). Over 50% of the area pertains to the river and the coastal plains; the remaining area incorporates the hilly relief up to ~300 m a.s.l. with a slope gradient of 5-30° (total area 39.7 km 2 ; Figure 5). With respect to the 2007 event, which partially affected the same area (lower Salinello hillslopes), the 2011 event occurred on a moderately vegetated landscape with agricultural areas (arable land, vineyards, and olive groves; Figure 5a) at an initial crop growth stage and during grass development. This event occurred again on clay-sand-conglomerate bedrock blanketed by a clayey eluvial and colluvial cover and landslides (Figure 5b).
The Pineto-Atri area consists of three coastal catchments, the intervening coastal slopes, and an up to <1 km wide coastal plain (total area 60.6 km 2 ; Figure 6). The elevation ranges from sea level to ~450 m a.s.l. The coastal and river plain areas make up ~15% of the area, while the hilly area (slope gradient 5-35° with local scarps) makes up the remaining 85%. Similar to the others, the Pineto area is characterized by arable land and high-quality plantations in the hilly areas, but also large bare clay areas within widespread badlands affecting the south-facing clay slopes. The coastal area is, again, seamlessly urbanized ( Figure 6a). The hillslopes are on clay bedrock; in the hilltops, sand, sandstone and conglomerate are present. As for the other areas, the slopes are covered with eluvial-colluvial covers and landslides (Figure 6b). The coastal area is, again, seamlessly urbanized ( Figure 6a). The hillslopes are on clay bedrock; in the hilltops, sand, sandstone and conglomerate are present. As for the other areas, the slopes are covered with eluvial-colluvial covers and landslides (Figure 6b).   During September 2012, after a >1-month dry time, the hilly coastal area of the northeastern Apennines (between the Pescara and Vibrata rivers; Figure 2), was affected by two rainfall events a short time from one another: a moderate one on 5-6 September and a heavy one on 13-14 September. After both of these events, the Pineto-Atri area was investigated in detail (the same three minor coastal catchments analyzed in 2011; Figure 6). The 2012 event occurred on poorly vegetated landscape with agricultural areas (arable land, vineyards, and olive groves) that were mainly plowed.

Methods
This work deals with the rainfall-induced soil erosion issue with a geomorphological approach based on direct field observation and remote air-photo analysis. It is based on the comparison of rain gauge data and geomorphological effects on the piedmont and coastal area of the NE Abruzzo region during three heavy rainfall events. We used a rainfall dataset from a network of 18 meteorological pluviometric stations (blue dot in Figure 2) provided by the Functional Center and Hydrographic Office of the Abruzzo Region. In the study area, 5-15 min sampled rain gauges (for at least 6 days around the main events, accuracy 0.2-1 mm) were analyzed, as well as daily data (1-2 months before the event). For each event, the data (from single rain gauges or from networks of rain gauges) enabled the analysis and comparison of: (1) hourly rainfall intensity; (2) event cumulative rainfall; (3) event duration; (4) daily rainfall; and (5) pre-event rainfall [2,43,65]. In this work, we do not consider the uncertainty in the rainfall duration triggering the geomorphological instabilities [65] because all the events are well defined in terms of time, and, focusing on the overall soil erosion, we consider all the effects induced by the entire events as an assumption.
The geomorphological effects of the heavy rainfall events were investigated through different types of available data: (1)  . This allowed us to be confident that the mapped effects were triggered by the events [66]. The field survey investigated the lithological features of bedrock and superficial deposits cover and the field data were compared to the results of borehole investigations and technical reports. The field geomorphological investigations focused on the type and distribution of the geomorphological effects (i.e., sheet-rill areas, gully areas, major gullies, channel incision, mud flows, flooding areas, crevasse splays, bank failures). For all of the landforms, the length, width, and surface area were mapped in the field (~5 m plan resolution by means of GPS measurements and mapping on 1:5000 topographic maps) and digitized in GIS software (ArcMap ® 10.1, ESRI, Redlands, CA, USA). The incision of the erosional landform and the thickness of the depositional landforms were directly measured in the field (~1-5 cm precision) on~55% of the landforms and inferred from remote observation on the remaining~45%. For the 2007 event, the remote observation from aerial photos (0.5-m resolution) allowed us to map, with GIS software, the surface distribution of the landforms and partly interpret their incision and thickness. Some landforms (from four to ten) for each type were observed and measured in the field for ground-truthing. Through these combined field and remote investigations, a dataset including (1) 1192 features for 2007 (~40 points of direct measurements), (2) 299 features for 2011 (~180 points of direct measurements), and (3) 614 features for 2012 (~330 points of direct measurements), was defined including soil erosion and flooding features induced by the three events. This allowed us to define or estimate the average erosion depth or accumulation thickness for all the landforms. Considering the surface area (for areal landforms), length and width (for linear features), and the average erosion depth or accumulation thickness of the landforms, the eroded volume, as well as the sedimented volume, were estimated for each landform for the overall investigated areas and for each event as follows: (1) areal landform erosion = ∑ i=1−n avg.erosion depth i × land f orm area i Water 2018, 10, 1314 10 of 29 The total eroded volume was calculated from (1) and (2), and the total accumulated sediment from (3).
The percentage and areal distribution of these effects were also analyzed for the different events and compared to the land use and crop conditions in the agricultural areas [42]. Averaging the eroded volumes in the investigated areas (average erosion = total eroded volume/hilly surface area) provided the average erosion depth for each event. In the calculation, only the hilly areas (with slope > 5 • ) were considered, excluding valley floors and fluvial and coastal plains.
We exploited soil erosion values and rainfall data in each event (heuristic expert method according to Brunetti et al. [66]) in order to outline a correlation and the expected threshold for the trigger of heavy soil erosion features in the investigated area. Finally, the erosion rates/rainfall ratio were compared to short-term and long-term erosion rates known for the Mediterranean area in order to outline the impact of heavy soil erosion on the clayey and sandy Apennine piedmont and coastal hilly landscape.

Results
The results of the rainfall distribution analysis and post-event investigations are presented for the 2007, 2011, and 2012 events. For each event, only the most significant rainfall graphs are presented as well as the mapping, features, and distribution of the geomorphological effects.

Rainfall Amount and Duration
The rainfall event occurred from about 20:00 to 21:00 on October 6 to about 11:00 to 12:00 on 7 October 2007, for a total duration of some 14-16 h. The cumulative rainfall was from 60-80 mm in the coastal area to >200 mm in the hilly area (205 mm, Nereto station; Figure 7a). The intensity was from moderate to high, with values ranging from 10 mm/h in the coastal area to 40 mm/h in the hilly area (Nereto station; Figure 7a). The daily rainfall was comparable to the cumulative for the event (~100 mm up to a maximum of 210 mm, Tortoreto and Nereto stations, Figures 7b and 8). Along the coast, the recorded daily precipitation was around 60-80 mm ( Figure 7b). Moreover, this event occurred after a fairly dry period, with very low September rainfall (<20 mm; Figure 8).
In summary, the Tortoreto 2007 event had high intensity (10-40 mm/h, up to >200 mm/day) and high cumulative rainfall (up to 210 mm) and occurring after >1 month of very low rainfall. The total eroded volume was calculated from (1) and (2), and the total accumulated sediment from (3).
The percentage and areal distribution of these effects were also analyzed for the different events and compared to the land use and crop conditions in the agricultural areas [42]. Averaging the eroded volumes in the investigated areas (average erosion = total eroded volume/hilly surface area) provided the average erosion depth for each event. In the calculation, only the hilly areas (with slope > 5°) were considered, excluding valley floors and fluvial and coastal plains.
We exploited soil erosion values and rainfall data in each event (heuristic expert method according to Brunetti et al. [66]) in order to outline a correlation and the expected threshold for the trigger of heavy soil erosion features in the investigated area. Finally, the erosion rates/rainfall ratio were compared to short-term and long-term erosion rates known for the Mediterranean area in order to outline the impact of heavy soil erosion on the clayey and sandy Apennine piedmont and coastal hilly landscape.

Results
The results of the rainfall distribution analysis and post-event investigations are presented for the 2007, 2011, and 2012 events. For each event, only the most significant rainfall graphs are presented as well as the mapping, features, and distribution of the geomorphological effects.

Rainfall Amount and Duration
The rainfall event occurred from about 20:00 to 21:00 on October 6 to about 11:00 to 12:00 on 7 October 2007, for a total duration of some 14-16 h. The cumulative rainfall was from 60-80 mm in the coastal area to >200 mm in the hilly area (205 mm, Nereto station; Figure 7a). The intensity was from moderate to high, with values ranging from 10 mm/h in the coastal area to 40 mm/h in the hilly area (Nereto station; Figure 7a). The daily rainfall was comparable to the cumulative for the event (~100 mm up to a maximum of 210 mm, Tortoreto and Nereto stations, Figures 7b and 8). Along the coast, the recorded daily precipitation was around 60-80 mm ( Figure 7b). Moreover, this event occurred after a fairly dry period, with very low September rainfall (<20 mm; Figure 8).
In summary, the Tortoreto 2007 event had high intensity (10-40 mm/h, up to >200 mm/day) and high cumulative rainfall (up to 210 mm) and occurring after >1 month of very low rainfall.

Geomorphological Effects
This event induced heavy soil erosion processes on the slopes (sheet, rill, and gully erosion), rapid mud flows at the bases of slopes and minor drainage basins, and flooding within the main river and coastal plains, mostly at the outlets of minor tributary catchments (Figures 9 and 10). On the downstream side, along the urbanized area and the coastal plain, large volumes of sediment were deposited on roads, at the stream outlet, and along paleovalleys (covered with urban areas). This induced serious problems at the circulation and also affected the safety of the population.
Hilly catchments and slopes were affected by gully erosion (4.36 km 2 , 10.0% of the investigated area; Figures 9 and 10) and sheet-rill erosion (4.84 km 2 , 11.1% of the investigated area; Figures 9 and 10). The low gradient slopes and ridges with vineyards and olive groves, or, in general, those not plowed, were mostly affected by moderate sheet-rill erosion. Rill features were observed as 1-5 cm deep and with 1-10 m spacing in these areas. The average erosion in the sheet-rill erosion areas was calculated to be as much as ~1 cm. The total calculated eroded volume was ~48,400 m 3 . The high gradient slopes, particularly where plowed in a downslope direction, were incised by gullies ( Figure  9). Gullies were measured mostly from 20 cm up to ~100 cm deep and 40-200 cm wide, usually spaced 4-8 m apart. The average erosion in the gully erosion area was calculated to reach ~14 cm. The total eroded volume was calculated as ~610,900 m 2 . Along the main channels of the tributary catchments, channel incision occurred (channels up to 2 m deep and 3-4 m wide). The area affected was 0.23 km 2 and the average depth was ~100 cm. The total eroded volume was calculated as ~228,400 m 3 .
Mud flows were the result of heavy soil erosion on the slopes, inducing the mobilization of huge sediment volumes of clay-silt and minor volumes of sand from the eluvial and colluvial cover. Extensive mud and water flooding affected both coastal plains, coming more from the small catchments than from the main rivers ( Figure 10). In the lower part of the slopes and at the outlets of the catchments, mud flows were very common (total area of ~0.37 km 2 , ~0.9% of the investigated area), showing mud accumulated thickness ranging from 50 to 200 cm. The total estimated sediment volume was 250,900 m 3 . Mud-rich flooding areas occurred at the bases of the slopes, along the channels of tributary catchments, and at the outlets of catchments to the fluvial or coastal plains (~3.97 km 2 , ~9.1% of the investigated area); accumulated mud was 0-10 cm thick, for a total estimated volume of 153,000 m 3 .

Geomorphological Effects
This event induced heavy soil erosion processes on the slopes (sheet, rill, and gully erosion), rapid mud flows at the bases of slopes and minor drainage basins, and flooding within the main river and coastal plains, mostly at the outlets of minor tributary catchments (Figures 9 and 10). On the downstream side, along the urbanized area and the coastal plain, large volumes of sediment were deposited on roads, at the stream outlet, and along paleovalleys (covered with urban areas). This induced serious problems at the circulation and also affected the safety of the population.
Hilly catchments and slopes were affected by gully erosion (4.36 km 2 , 10.0% of the investigated area; Figures 9 and 10) and sheet-rill erosion (4.84 km 2 , 11.1% of the investigated area; Figures 9 and 10). The low gradient slopes and ridges with vineyards and olive groves, or, in general, those not plowed, were mostly affected by moderate sheet-rill erosion. Rill features were observed as 1-5 cm deep and with 1-10 m spacing in these areas. The average erosion in the sheet-rill erosion areas was calculated to be as much as~1 cm. The total calculated eroded volume was~48,400 m 3 . The high gradient slopes, particularly where plowed in a downslope direction, were incised by gullies (Figure 9). Gullies were measured mostly from 20 cm up to~100 cm deep and 40-200 cm wide, usually spaced 4-8 m apart. The average erosion in the gully erosion area was calculated to reach~14 cm. The total eroded volume was calculated as~610,900 m 2 . Along the main channels of the tributary catchments, channel incision occurred (channels up to 2 m deep and 3-4 m wide). The area affected was 0.23 km 2 and the average depth was~100 cm. The total eroded volume was calculated as~228,400 m 3 .
Mud flows were the result of heavy soil erosion on the slopes, inducing the mobilization of huge sediment volumes of clay-silt and minor volumes of sand from the eluvial and colluvial cover. Extensive mud and water flooding affected both coastal plains, coming more from the small catchments than from the main rivers ( Figure 10). In the lower part of the slopes and at the outlets of the catchments, mud flows were very common (total area of~0.37 km 2 ,~0.9% of the investigated area), showing mud accumulated thickness ranging from 50 to 200 cm. The total estimated sediment volume was 250,900 m 3 . Mud-rich flooding areas occurred at the bases of the slopes, along the channels of tributary catchments, and at the outlets of catchments to the fluvial or coastal plains (~3.97 km 2 ,~9.1% of the investigated area); accumulated mud was 0-10 cm thick, for a total estimated volume of 153,000 m 3 . Water 2018, 10, x FOR PEER REVIEW 12 of 28

Rainfall Amount and Duration
The rainfall event occurred from about 02:00 on 1 March to 02:00 on 2 March 2011, for a total duration of some 22-26 h. The cumulative rainfall was 100-130 mm at most of the stations (e.g., Pineto station; Figure 11a), up to a maximum of 211 mm in the hilly area (Nereto station; Figure 11b). The rainfall intensity was around 15-20 mm/h and up to 35 mm/h (e.g., Pineto and Nereto stations; Figure 11). The daily rainfall was~80-120 mm/day in the entire hilly area, up to a maximum of 180 mm (Nereto station; Figure 12a). Moreover, this event occurred after a moderately humid winter period (~50 mm in a 10-day time span before the event, Figure 12b).
In summary, the 2011 event affected a regional area (Pineto 2011 and Salinello 2011) for a moderately short duration (22-26 h) with high intensity (15-35 mm/h, up to >180 mm/day) and high cumulative rainfall (up to 211 mm) and after~10 days of moderate antecedent rainfall.
Water 2018, 10, x FOR PEER REVIEW 13 of 28

Rainfall Amount and Duration
The rainfall event occurred from about 02:00 on 1 March to 02:00 on 2 March 2011, for a total duration of some 22-26 h. The cumulative rainfall was 100-130 mm at most of the stations (e.g., Pineto station; Figure 11a), up to a maximum of 211 mm in the hilly area (Nereto station; Figure 11b). The rainfall intensity was around 15-20 mm/h and up to 35 mm/h (e.g., Pineto and Nereto stations; Figure  11). The daily rainfall was ~80-120 mm/d in the entire hilly area, up to a maximum of 180 mm (Nereto station; Figure 12a). Moreover, this event occurred after a moderately humid winter period (~50 mm in a 10-day time span before the event, Figure 12b).
In summary, the 2011 event affected a regional area (Pineto 2011 and Salinello 2011) for a moderately short duration (22-26 h) with high intensity (15-35 mm/h, up to >180 mm/d) and high cumulative rainfall (up to 211 mm) and after ~10 days of moderate antecedent rainfall.

Geomorphological Effects
The geomorphological analysis performed after the event outlined soil erosion due to sheet and rill erosion, gully erosion, and sedimentation as mud flows and flooding (Figure 13), similar to the 2007 event but with a different distribution and minor total extent. In the hilly catchments and slopes of the Pineto-Atri area (Figures 13a,b, and 14a), soil erosion occurred as gully areas, with gullies 20-50 cm deep and spaced 3-20 m (over an area of 1.13 km 2 , 1.9% of the investigated area, Figure 14); the total eroded volume was 133,000 m 3 . Major gullies up to 1 m deep and 2 m wide (total gully length

Rainfall Amount and Duration
The rainfall event occurred from about 02:00 on 1 March to 02:00 on 2 March 2011, for a total duration of some 22-26 h. The cumulative rainfall was 100-130 mm at most of the stations (e.g., Pineto station; Figure 11a), up to a maximum of 211 mm in the hilly area (Nereto station; Figure 11b). The rainfall intensity was around 15-20 mm/h and up to 35 mm/h (e.g., Pineto and Nereto stations; Figure  11). The daily rainfall was ~80-120 mm/d in the entire hilly area, up to a maximum of 180 mm (Nereto station; Figure 12a). Moreover, this event occurred after a moderately humid winter period (~50 mm in a 10-day time span before the event, Figure 12b).
In summary, the 2011 event affected a regional area (Pineto 2011 and Salinello 2011) for a moderately short duration (22-26 h) with high intensity (15-35 mm/h, up to >180 mm/d) and high cumulative rainfall (up to 211 mm) and after ~10 days of moderate antecedent rainfall.

Geomorphological Effects
The geomorphological analysis performed after the event outlined soil erosion due to sheet and rill erosion, gully erosion, and sedimentation as mud flows and flooding (Figure 13), similar to the 2007 event but with a different distribution and minor total extent. In the hilly catchments and slopes of the Pineto-Atri area (Figures 13a,b, and 14a), soil erosion occurred as gully areas, with gullies 20-50 cm deep and spaced 3-20 m (over an area of 1.13 km 2 , 1.9% of the investigated area, Figure 14); the total eroded volume was 133,000 m 3 . Major gullies up to 1 m deep and 2 m wide (total gully length

Geomorphological Effects
The geomorphological analysis performed after the event outlined soil erosion due to sheet and rill erosion, gully erosion, and sedimentation as mud flows and flooding (Figure 13), similar to the 2007 event but with a different distribution and minor total extent. In the hilly catchments and slopes of the Pineto-Atri area (Figure 13a,b and Figure 14a), soil erosion occurred as gully areas, with gullies 20-50 cm deep and spaced 3-20 m (over an area of 1.13 km 2 , 1.9% of the investigated area, Figure 14); the total eroded volume was 133,000 m 3 . Major gullies up to 1 m deep and 2 m wide (total gully length 10.2 km; Table 1a, Figure 14a) occurred along the channels of tributary catchments and slope undulations, in some cases enlarging natural or manmade notches or plowing incisions (Figure 13b); the total estimated volume was 8900 m 3 . At the bases of the slopes within the main valleys and along the coastal slope rapid mud flows occurred from 20 cm to 150 cm thick (total area of 0.26 km 2 , 0.4% of the investigated area; Table 1a, Figure 14a); the estimated accumulated volume was 60,400 m 3 . Along the main rivers and at the outlet to the coastal plain, crevasse splays occurred, with a sediment thickness of 10-50 cm (total area 0.16 km 2 , 0.3% of the area; estimated volume 20,000 m 3 ), as well as flooding areas, with mud 0-5 cm thick (total area 2.01 km 2 , 3.3% of the area; estimated volume 26,100 m 3 ).  In the Salinello River valley (Figures 13c,d,e and 14b), limited soil erosion occurred on the moderately vegetated slopes, while most of the effects were related to flooding along the main river. Heavy mud and water flooding were the prevailing effects of this event, affecting the coastal plains as well as the river plain almost seamlessly. Extensive overbank flooding along the rivers induced the formation of wide and long crevasse splays on the floodplain (Figure 13c-e). Soil erosion occurred In the Salinello River valley (Figure 13c,d,e and Figure 14b), limited soil erosion occurred on the moderately vegetated slopes, while most of the effects were related to flooding along the main river. Heavy mud and water flooding were the prevailing effects of this event, affecting the coastal plains as well as the river plain almost seamlessly. Extensive overbank flooding along the rivers induced the formation of wide and long crevasse splays on the floodplain (Figure 13c-e). Soil erosion occurred in terms of areas affected by sheet-rill erosion areas (Table 1b, Figure 14b), with an average erosion of 5-10 cm (0.10 km 2 , 0.2% of the area; estimated volume 5700 m 3 ), with gullies 20-100 cm deep and 0.5-2.0 m wide (for a total length of~4.5 km; estimated eroded volume 4400 m 3 ). Along the main channel of the Salinello River bank, failure occurred (for a total length of~2.0 km; Figure 13c,d), inducing severe damage to valley roads. At the bases of the slopes, mud flows were formed (Table 1b, Figure 14b) over a 0.15 km 2 area, as thick as 50-150 cm (total estimated volume 65,500 m 3 ). Flooding areas occurred (Table 1b, Figure 14b)

Rainfall Amount and Duration
This event is the result of two rainfall periods separated by a 7-day dry period and has to be considered as two separated events (according to Peruccacci et al. and Brunetti et al. [7,66]). The first rainfall event occurred from about 17:00 on 5 September to 11:00 on 6 September 2012, for a total duration of some~18 h, with peaks in the evening and the morning; the second event occurred from about 18:00 on 13 September to 18:00 on 14 September 2012, for a total duration of~24 h. The cumulative rainfall was~30-110 mm for the first event and~80-190 mm for the second event (Figure 15a,b). The recorded rainfall intensity was 10 to >60 mm/h for the first event and~15-45 mm/h for the second event (Figure 15a,b). The daily rainfall was~60-110 mm/day for the first event and up to a maximum of 190 mm for the second event (Figure 15c,d). Moreover, these events occurred after more than a month of completely dry conditions (Figure 15c).
In summary, the Pineto 2012 event was a double event, 18 and 24 h long, and affected and coastal hilly area within a 1-week interval, with high intensity (15-65 mm/h, up to >190 mm/day) and combined cumulative rainfall up to 280 mm, after a >1-month completely dry period.

Geomorphological Effects
The geomorphological analysis performed after both the first and second event outlined a large distribution of soil erosion features combined with mud flows, flooding areas, and minor crevasse splays ( Figure 16).
On 13-14 September (Figure 17b, Table 2b), further heavy soil erosion occurred, mostly due to gully areas with an average 20-40 cm incision (over 1.18 km 2 , 2.0% of the area; estimated volume 236,500 m 3 ) and major gullies 20-70 cm deep and 1-2 m wide (total length 48.8 km; estimated volume 23,300 m 3 ; Figure 16 d-f). Sheet-rill erosion areas were observed with a 2-5 cm incision (over 0.92 km 2 , 1.5% of the area; estimated volume 27,600 m 3 ). In the lower part of the slopes and at the junctions of minor tributaries to main streams, mud flows occurred up to 150 cm thick (over 0.18 km 2 , 0.3% of the area; estimated volume 88,000 m 3 ). Along the main fluvial plain and the coastal plain, flooding areas occurred with 0-3 cm of accumulated mud (over 0.36 km 2 , 0.59% of the area; estimated volume ~5000 m 3 ), as well as some local crevasse splay 10-30 cm thick (over 0.03 km 2 , 0.1% of the area; estimated volume 5900 m 3 ).

Geomorphological Effects
The geomorphological analysis performed after both the first and second event outlined a large distribution of soil erosion features combined with mud flows, flooding areas, and minor crevasse splays ( Figure 16).
On 13-14 September (Figure 17b, Table 2b), further heavy soil erosion occurred, mostly due to gully areas with an average 20-40 cm incision (over 1.18 km 2 , 2.0% of the area; estimated volume 236,500 m 3 ) and major gullies 20-70 cm deep and 1-2 m wide (total length 48.8 km; estimated volume 23,300 m 3 ; Figure 16 d-f). Sheet-rill erosion areas were observed with a 2-5 cm incision (over 0.92 km 2 , 1.5% of the area; estimated volume 27,600 m 3 ). In the lower part of the slopes and at the junctions of minor tributaries to main streams, mud flows occurred up to 150 cm thick (over 0.18 km 2 , 0.3% of the area; estimated volume 88,000 m 3 ). Along the main fluvial plain and the coastal plain, flooding areas occurred with 0-3 cm of accumulated mud (over 0.36 km 2 , 0.59% of the area; estimated volumẽ 5000 m 3 ), as well as some local crevasse splay 10-30 cm thick (over 0.03 km 2 , 0.1% of the area; estimated volume 5900 m 3 ).

Discussion
A direct investigation (field mapping and air-photo) analysis of the geomorphological effects induced by three heavy rainfall events in the NE Abruzzo hilly coastal area provided the mapping of >2000 soil erosion and accumulation features. The amount of soil erosion was estimated along the slopes and in the hilly areas for the different types of landforms in all of the investigated events ( Table 3). The distribution of mud flows, flooding areas, and crevasse splay outlined the estimated volume of sediment accumulated at the bases of slopes, at the outlets of minor catchments, and on the river and coastal plains (Table 4).
For the 2007 event, the estimated soil erosion volume on the hillslopes and catchments was about 887,700 m 3 (mostly due to gully erosion; Table 3). The average overall erosion calculated in the hilly area (65% of the investigated area) was~3.08 cm in this single event. At the base of the slopes, in the valley floors, and in the coastal and river plains, mud flows and floods accumulated an estimated 403,800 m 3 volume of sediment. This means that more than half of the eroded soil (55%) was removed from the hilly system to the sea. For the 2011 event in the Pineto-Atri area, the estimated amount of soil erosion was about 141,900 m 3 (mostly due to gully erosion; Table 3). The average overall erosion calculated in the slopes of the hilly area (~85% of the investigated area) was~0.29 cm. At the bases of slopes, in the valley floors, and in the coastal and river plains, mud flows and floods accumulated an estimated 106,500 m 3 volume of sediment. Again, 25% of the eroded soil was removed from the hilly system to the sea. In the Salinello River area, the estimated eroded volume was very low at 12,100 m 3 (overall erosion estimated in the hilly area <0.01 cm), while the sedimented volume was as great as 196,000 m 3 (mostly due to mud flows and crevasse splay; Table 3). This confirms that this case is in a different geomorphological framework (main river valley instead of coastal hilly area) and the main geomorphological effects were related to flooding in the main river (mostly coming from upstream and the main valley) and not to soil erosion in the hilly slopes. For the 2012 double event, the estimated amount of soil erosion (Table 3) was about 36,500 m 3 on 5-6 September (Pineto 2012-1) and 289,900 m 3 on 13-14 September (Pineto 2012-2). The average overall erosion in the slopes of the hilly Pineto-Atri area (~85% of the investigated area) was~0.08 (5-6 September) and~0.60 cm (13)(14). At the bases of slopes, in the valley floors, and in the coastal and river plains, mud flows and floods accumulated an estimated amount of about 37,500 m 3 of sediment (5-6 September, roughly corresponding to the eroded volume) and 98,900 m 3 of sediment (13)(14) September, about one-third of the eroded volume). Again, this means that 66% of the eroded soil was removed from the hilly system to the sea. Moreover, for all the events, considering that large part of the sediment was accumulated in urban/industrial areas and removed after the events, the actual soil loss was even higher.  The events were characterized by variable heaviness (i.e., cumulate rain and intensity) and different areas, but in comparable lithological-geomorphological settings (coastal clay-sand-conglomerate hills covered with eluvial-colluvial and landslide deposits), except for the Salinello River valley (major river valley and hillslopes). The main rainfall parameters of the events in the different areas (i.e., cumulative rainfall, which in this case roughly corresponds to daily rainfall, rainfall peak intensity, and rainfall duration) were compared to the soil erosion estimations (eroded volume, erosion averaged over the affected hilly slope areas, sediment volume, sediment loss) ( Table 4). Despite encompassing a limited number of investigated events, rainfall data and soil erosion estimations were analyzed ( Figure 18) in order to outline the controlling factors and possibly a triggering threshold for the heavy soil erosion features in this area. Considering the distribution of rainfall values, the strong influence on soil erosion of cumulative rainfall and the poor influence of maximum rainfall intensity are outlined in this case. More specifically, eroded volume vs., cumulative rainfall (Figure 18a), which roughly corresponds to daily rainfall for these events as well as average erosion vs. cumulative rainfall (Figure 18b), shows good correlation (the latter not depending on the size of the areas). Otherwise, average erosion vs. rainfall intensity (Figure 18c) shows a scattered distribution and poor control of the peak intensity in soil erosion for these events. Only the 2011 Salinello event is completely incomparable (2a in Figure 18), as expected, being in a different geomorphological framework (main river valley vs., coastal hilly area) and possibly due to the condition of slopes in terms of vegetation cover. The event occurred at the beginning of March, with the cultivated areas already covered by incipient vegetation, which prevented the slopes from undergoing heavy soil erosion. This is partially true also for the Pineto 2011 area (in terms of landform distribution), in which sheet-rill features did not form due to the vegetation cover of the slope, and only gully areas and major gullies incised the slopes.
Several studies outlined that the main parameters controlling the soil erosion are related to the interaction among rainfall parameters and lithology, orography, hydrography, land-use and vegetation [5][6][7][8][9][10]. In this case, considering the relatively homogeneous lithologic-geomorphological settings, we mainly outlined the role of land-use and vegetation cover (as mentioned above) and the control of rainfall parameters (i.e., cumulative rainfall) on the soil erosion.
etation cover. The event occurred at the beginning of March, with the cultivated areas already covered by incipient vegetation, which prevented the slopes from undergoing heavy soil erosion. This is partially true also for the Pineto 2011 area (in terms of landform distribution), in which sheet-rill features did not form due to the vegetation cover of the slope, and only gully areas and major gullies incised the slopes.
Several studies outlined that the main parameters controlling the soil erosion are related to the interaction among rainfall parameters and lithology, orography, hydrography, land-use and vegetation [5][6][7][8][9][10]. In this case, considering the relatively homogeneous lithologic-geomorphological settings, we mainly outlined the role of land-use and vegetation cover (as mentioned above) and the control of rainfall parameters (i.e., cumulative rainfall) on the soil erosion. Figure 18. Comparison of soil erosion values (eroded volume, average erosion, sediment loss) and rainfall distribution (cumulative rainfall, max and mean rainfall intensity). The event numbers refer to Table 3.
We compared the average single-event erosion estimated in the investigated events (0.8-30.8 mm) with the short-and long-term erosion rates known in the literature (0.05-90 mm/year [67][68][69][70][71][72][73][74][75]; Table 5). The values resulting from this work are comparable to the highest short-term values known for badlands (10-30 mm/year) and landslides (30-90 mm/year) from direct measurements at local scale and in small catchments [29,39,41] and about one order of magnitude higher than the values known at basin scale (0.05-8 mm/year) obtained from indirect estimations [37,69,75] and the longterm erosion rates (0.06-15 mm/year) derived from trapped sediment calculations and morphotectonic reconstructions [41,[67][68][69][70] (Table 5). Even if these values are related to different types of processes (episodic vs., continue processes), taking into account that they affect comparable lithologic, geomorphological, and climatic environments (or in some cases, the same area, 10 in Table 5), this comparison highlights that the impact of soil erosion on the landscape is largely underestimated, particularly in the long-term analysis, in respect to what is outlined by the recent heavy rainfall events. Figure 18. Comparison of soil erosion values (eroded volume, average erosion, sediment loss) and rainfall distribution (cumulative rainfall, max and mean rainfall intensity). The event numbers refer to Table 3.
We compared the average single-event erosion estimated in the investigated events (0.8-30.8 mm) with the short-and long-term erosion rates known in the literature (0.05-90 mm/year [67][68][69][70][71][72][73][74][75]; Table 5). The values resulting from this work are comparable to the highest short-term values known for badlands (10-30 mm/year) and landslides (30-90 mm/year) from direct measurements at local scale and in small catchments [29,39,41] and about one order of magnitude higher than the values known at basin scale (0.05-8 mm/year) obtained from indirect estimations [37,69,75] and the long-term erosion rates (0.06-15 mm/year) derived from trapped sediment calculations and morphotectonic reconstructions [41,[67][68][69][70] (Table 5). Even if these values are related to different types of processes (episodic vs., continue processes), taking into account that they affect comparable lithologic, geomorphological, and climatic environments (or in some cases, the same area, 10 in Table 5), this comparison highlights that the impact of soil erosion on the landscape is largely underestimated, particularly in the long-term analysis, in respect to what is outlined by the recent heavy rainfall events. Concerning the issue of the rainfall threshold for geomorphological effects of heavy rainfall events, which is widely debated in the literature ( [43,70,71] and references therein), this work provides a local contribution that could be representative of a larger area. The rainfall thresholds for soil erosion are usually defined for short-term intensities (e.g., from 5.6 to 20 mm/h [70,71]). In this work, even after a few investigated cases, the average erosion vs. cumulative rainfall graph (Figure 18b) suggests a threshold for the triggering of heavy soil erosion features around 100-110 mm of cumulative rainfall in~1-day events, according to the intersection of the correlation line with the x-axis (excluding the 2011 Salinello event, 2a; Figure 18). The hourly intensity threshold is less defined in this case (scattered distribution in Figure 18b) at roughly~30-40 mm/h. These values are more meaningful for agricultural areas (e.g., arable lands, vineyards, olive groves) during the plowed stage (in autumn and winter) than during the unplowed vegetative stage (spring and summer), when soil erosion distribution is less intense. These results can also be representative of larger areas, considering that the geologic and geomorphological features of the investigated areas are typical of large parts of the Adriatic clay hills.
When compared with other direct and indirect methods applied to soil erosion [36][37][38][39][40][41][42], the direct field survey approach appears to be suitable for intermediate-scale (several km 2 ) investigations and is useful for assessing the impact of soil erosion to hilly and agricultural landscapes. Other direct observations provide greater advantages at the local scale (~1 km 2 and less) and in small basin investigations, while indirect estimations are more suitable for large-scale (hundreds of km 2 ) soil erosion assessment of large drainage basins and entire physiographic domains. Moreover, the used approach is more suitable for single or multiple event investigations as in this work, and for applied issues, while other direct and indirect methods provide average long-term and short-term assessments that, however, might underestimate the assessment of soil erosion and its impact on the landscape.

Conclusions
In this work, we investigated three~1-day heavy rainfall events (combined or not with antecedent rainfall) in 2007, 2011, and 2012 in different areas along the clayey hilly-coastal NE Abruzzo area, largely characterized by agricultural areas and urban tourist areas (affected by strong urbanization along coast and hilly area). The events were investigated in terms of cumulative rainfall, intensity and duration, triggered geomorphological effects (soil erosion and accumulation), and average erosion. The investigation offers several contributions to widely debated issues: (1) soil erosion assessment and impact on clayey landscapes typical of the entire Adriatic coastal hilly area; (2) comparison of estimated erosion with short-and long-term erosion rates in the Mediterranean area; and (3) rainfall triggering threshold for heavy soil erosion.
The soil erosion assessment shows average erosion from 0.08 cm to 3.08 cm in a single~1-day heavy rainfall event ( Figure 18). The erosion values show a good comparison (even if over a few events) with cumulative rainfall and not with peak rainfall intensity. The correlation is significant for the hilly coastal areas (Tortoreto, Pineto) but not for the Salinello area, incorporating the lower part of a main river valley and processes that are more connected to the main river than to the hillslopes. More specifically, the geomorphological features of the 2011 event, which occurred in March with partly vegetated cropland, outline the positive contribution (reduced erosion values) of vegetated landscape for preventing slopes from heavy soil erosion with respect to other events occurring in a heavily plowed landscape.
The average erosion values obtained for single events are comparable to the annual erosion rates obtained from previous studies for badlands (10-30 mm/y) and up to about one order of magnitude greater than the basin scale values (0.05-8 mm/y) and the long-term erosion rates (0.06-15 mm/y). Since these events occurred several times in the last decades and heavy rainfall is one of the main issues in future scenarios in terms of climate change, it is possible to outline the underestimation of the impact of soil erosion on the Adriatic hilly landscapes. This is fundamental in areas largely dedicated to agriculture, including high-quality plantations (e.g., olive groves, vineyards), and urban tourist areas, which are also suffering a strong loss of free soil due to land-use change and urbanization in recent years.
In terms of a heavy soil erosion triggering thresholds for clayey hilly landscapes, the correlation of cumulative rainfall and average erosion shows an expected threshold value as high as~100-110 mm/day ( Figure 18). This value, even if it results from a few local cases, could be representative of larger areas of the Adriatic clay hills with a similar morphoclimatic setting and it is comparable with the results of previous works, although mostly defining hourly rainfall intensity thresholds.
Finally, comparing the field mapping approach used in this work with other direct and indirect approaches for soil erosion assessment, it appears to be suitable mostly for intermediate-scale investigations and for single-multiple events analysis and can be used to verify and calibrate the evaluation of erosion rates with other methods.
In further development of this research, we expect to analyze new events, which occasionally affect the Adriatic clayey and sandy hills in order to verify the results. Moreover, we expect to compare the direct investigations with erosion modeling, and to use this approach to outline future scenarios of soil erosion impact in the framework of climate changes in terms of the debated issue of increasing heavy rainfall events.