Next Article in Journal
Variations of Drinking Water Quality Influenced by Seasons and Household Interventions: A Case Study from Rural Maharashtra, India
Previous Article in Journal
Biotrickling Filtration of Air Contaminated with 1-Butanol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Upscaling Plot-Scale Soil Erosion under Mediterranean Climate Variability

1
Met European Research Observatory, HyMex Network, 82100 Benevento, Italy
2
Department of Sciences and Technologies, University of Sannio, 82100 Benevento, Italy
3
French National Institute for Agricultural Research, Grassland Ecosystem Research Unit, 63000 Clermont-Ferrand, France
*
Author to whom correspondence should be addressed.
Submission received: 24 June 2017 / Revised: 6 August 2017 / Accepted: 13 August 2017 / Published: 21 August 2017

Abstract

:
Soil erosion is an issue in the Mediterranean slopes. Erosion plots are useful to quantify erosion rates, but data are difficult to scale up to a slope level. Moreover, short observational frameworks are generally established, making it difficult to represent multi-year fluctuations. This paper deals with the potential of parsimonious modelling to upscale plot erosion (~23 m2) at Monte Pino Met European Research Observatory (South Italy) from 2001 to 2006. Under the assumption that the slope is fractal and contains plots, monthly gross soil erosion was modeled by lumping together the erosivity factor (runoff component), Normalized Difference Vegetation Index (vegetation cover factor), and the spatial scale dependence (slope length factor). This model was applied to reconstruct monthly gross soil erosion rates for the period of 1986–2006, for which hydrological inputs were available with sufficient detail. Pronounced interannual variations, with two distinct patterns, were observed: increasing rates of erosion were visible in 1995–2006 (peaking in November 1997, 50 Mg·ha−1·month−1), while in previous years only a few peaks slightly exceeded the average of the whole period (1 Mg·ha−1·month−1). Hydrological conditions indicate that important erosional processes have been triggered during low-frequency, short rainfall events occurring in spring–summer (e.g., May 2001, June 2003), or during longer, less intense events occurring in autumn–winter (e.g., November 1997) seasons. It is likely that increased precipitation amounts associated with more frequent convective storms created conditions for higher energy events triggering erosion. For the recent warm period, investigations at a higher than monthly resolution are required to better assess the seasonal changes of erosion rates and their relationship with soil conservation.

1. Introduction

Soil erosion is a form of land degradation caused by complex and not always well known interactions among rainfall impact and variability, overland flow, land-use and land-cover, soil properties, topography, and conservation management [1]. In agricultural areas, tillage impact is noticeable on soil’s physical, chemical, and biological properties, which includes the effect on the soil environment in the form of runoff and soil erosion [2,3]. Extreme precipitation events can lead to considerable soil loss, with implications for downstream flood risk, and shallow landslides [4]. Soil erosion is a great concern all over the world, and a primary problem in many Mediterranean-climate areas of Southern Europe, which are characterized by typical within-year seasonality in precipitation distribution, with heavy rainfalls in fall and spring, and a relatively large year-to-year climate variability [5,6]. This feature affects, in turn, the local hydrological conditions with related instabilities in the vegetation ecosystems when severe perturbations occur [7]. In Mediterranean environments, soil erosion is mainly driven by raindrop impact and overland flow, and induces a loss of cropland available for food production, and is a major cause of various forms of land degradation [8,9,10,11]. Such meteorological and hydrological events exacerbate soil erosion phenomena, especially in those locations where the soil parent material is formed by a weakly consolidated substrate, and observed changes in local climate conditions reveal an increase in rainfall intensity [12,13,14]. A rapid evolution towards the conditions of severe soil erosion degradation is ongoing in the Mediterranean belt [15,16,17], including Southern Italy [18,19,20].
Soil erosion assessment and mitigation requires large datasets to understand the fluctuations in time and space, and is particularly challenging under climate change-driven shifts in precipitation [21]. For the Euro-Mediterranean region, only long-term averages of sediment yield have been predicted by the PESERA (Pan-European Soil Erosion Risk Assessment) project (see Figure 1; [22]). Available studies about the recurrence of dry periods and extreme rainfall events, which affect the variations in monthly soil erosion rates [23], should be extended over longer time intervals, and carried out at different temporal scales (e.g., inter-annual and inter-monthly scales) in order to be able to highlight not only the differences between the spring–summer and autumn precipitation, but also the seasonal variability in rainfall frequency and intensity. Month-time intervals seem an appropriate timescale to identify both short-and-long-term changes in soil erosion and its effects on ecosystem evolution. However, the complex interaction that exists between factors governing soil erosion processes (e.g., soil physical and hydraulic properties, local terrain attributes, and land-use) makes soil erosion predictions very challenging [24]. Whereas the available literature and databases are sources of good climatic data at regional and sub-regional scales for the last century, they do not provide enough information on space-time scale, location, and timing of geomorphologic changes (except in the general sense of rivers or coastline modifications [25]). Soil losses were quantified in various Mediterranean landscapes by using erosion plots [26,27,28,29,30,31]. The results of these studies have provided important reference outcomes that have had a fruitful repercussion in a range of related disciplines, such as ecology, geomorphology, and hydrology [32,33]. However, monitoring of soil erosion phenomena is often restricted to a few years of observations only, implying that the measurements are not necessarily representative of middle-term fluctuations occurring in ecosystems, especially under a Mediterranean climate [34,35].
We investigated the possibility of reconstructing middle-term monthly gross soil erosion on the basis of discrete plot measurements and subsequent upscaling of plot-scale modeling results to a slope-scale. This is a key issue for soil loss predictions also considering that soil erosion, as well as other environmental phenomena, is often monitored only across restricted areas. For this purpose, an empirical model is developed and calibrated at a plot-scale, and then evaluated to assess its effectiveness of transfer information to the slope scale though a GIS-integrated procedure. The model is used to reconstruct monthly gross soil erosion rates for the period of 1986–2006, and investigates possible signals of monthly erosion change at the Monte Pino plot near the city of Benevento in Southern Italy (Figure 2). Climate extremes are also considered. Under the conditions of climate change, a divergence of precipitation patterns is predicted. Some regions of Mediterranean Central Area (MCA) will experience more annual precipitation and a more regular seasonal distribution of rainfall events, while other areas will experience stationarity in annual amounts of precipitation and a higher rainfall intensity.

2. The Monte Pino Experimental Site (MPES)

The Monte Pino experimental site (MPES) is located in the Calore River Basin, in the Samnium territory (14°42′30′′ long. East and 41°06′30′′ lat. North, 184 m above the sea level; Benevento, southern Italy; Figure 2a,b). It is placed between the Benevento valley (northeast) and Montemauro Mountain (southwest; dotted square in Figure 2c) along the eastern slope of the Monte Pino (hence the name of the experimental site). The MPES comprises of an experimental plot of 23 m2 and a meteorological station. The meteorological station is a part of the Met European Research Observatory (MetEROBS) of the HyMex network [37], and is located, together with the erosion plot, in the extreme east mountainside (little white square in Figure 2c) that is also part of the regional park of the Taburno and Camposauro Mountains. The HyMex project develops and facilitates hydrological monitoring in European areas. This site offers a unique opportunity for modeling erosion responses to climatic variability in an area that has not been subject to drastic land-use change over the last century.
The images of Figure 3 illustrate the eastern slope of MPES as a heterogeneously structured entity, having a common climate (Figure 3a) with complex pedo-topographic and geological features, and a variable land-use (Figure 3b) controlling their hydrologic dynamics. The area has a typical Mediterranean climate, which allows growing high-value crops such as grapes and olives to produce high-quality wine and olive oil. In particular, mean annual rainfall is 950 mm (±150 mm Std.Dev.), and mean annual temperature is 14 °C (±0.5 °C Std.Dev.), both computed using the available datasets from 1986 to 2006. The MP site has a surface area of about 10 km2 and a mean elevation of 275 m above sea level (Figure 4a). The mean slope steepness is 16%, with hillsides mainly oriented towards the northwest, sloping to around 30% (Figure 4b). A plot of about 23 m2 was established along the hillslope, in an area characterized by a slope gradient of 32% and cultivated with olives trees (Figure 4c).

3. Data and Model

3.1. Data

We used climatic data retrieved from 1986 to 2006 from the meteorological station, installed within the Monte Pino experimental site, and soil erosion rates collected from 2001 to 2006 at the 23 m2 experimental plot. Rainfall intensity was measured using a tipping bucket rain gauge every 24 h, whereas the soil erosion rate was measured after each rainfall event. The tipping bucket was installed on a pole that was 1 m from the ground. The plot was tilled once each year (usually in spring) to reproduce the MPES land-use. This also allows for reducing the risk for the plot to exhaust the available material for soil detachment over long-term studies. Runoff was drained into a sedimentation drainage installed in the ground in the lower part of the plot. Eroded soil (kg·m−2) was collected once per month. Measured soil erosion rates are reported in Table 1. Both sediment and rainfall data were aggregated on a monthly basis. For instance, soil loss measured in the plot, during the years 2001 and 2002 was of 28 kg·m−2 and 0.55 kg·m−2, respectively. From 2001–2006 the median erosion values were equal to 2.52 kg·m−2. This example indicates the importance of sparse disastrously rainy days or months in the determination of medium-term erosion, and the need to explore the entire distribution of the erosion rates by months over a long period. Model upscaling was completed into a GIS environment using data derived from a 20 m2 square-cell digital elevation model (DEM) and a color orthophoto. The DEM was used to derive topographic attributes (Figure 4b), whereas the orthophoto was used to obtain values for the Normalized Difference Vegetation Index (−1 ≤ NDVI ≤ 1), a metric of vegetation greenness (usually derived from satellite-driven information), with higher values indicating highly vegetated areas.

3.2. Model

The Universal Soil Loss Equation (USLE) and its modified versions, such as RUSLE or MUSLE, are definitely the most known and used regression-based soil erosion models [38,39,40], and also provide the basis of process-based hillslope models, such as the USDA Water Erosion Prediction Project (WEPP; [41] and Eurosem [42]). Process-based models are data-demanding, and their applications at relatively large spatial scales can become a drawback developed a model [43] that is aimed at predicting monthly soil erosion caused by rainfall and overland flow across a plot-cell of 25 m × 25 m and is cast as follows:
E    =    k Q 2 S 1.67 e 0.07 v
where E is monthly soil erosion (mm·month−1), k is the erodibility index, Q is overland flow (mm·month−1), S is slope steepness (m·m−1), and v is vegetation cover (%).
In our study, a scale transfer tool is required to expand Thornes’s model [43] from the plot scale (23 m2) to the slope scale (an area of about 10 km2). Since the spatial averaging is basically equivalent to the spatial integration [44,45], a viable possible upscaling technique is by averaging the plot-scale model over the entire mountainside area [46]. This approach is employed in our study and based on the assumption that the plot is a fractal-element of the slope. Accordingly, Thornes’s model [43] is revised to generate an adapted RUSLE-Thornes model, which replaces runoff, Q, with RUSLE erosivity-factor, EI, and vegetation cover, v, with the Normalized Difference Vegetation Index, NDVIm. Therefore, the time-scale dependence of erodibility was expanded to include spatial scale-dependence of slope length factor (L).
These attributes are lumped together into the following semi-empirical, nonlinear relation:
E ( ART ) m    =    [ k m ( E I m ) η L ( γ d o m s 1 D ) ] e ( ω 100 N D V I m )
where E(ART)m is the gross soil erosion rate in kg·m−2·month−1. In Equation (2), the roughness erosive-resistance (the exponential term in the right-hand side), described as a function of vegetation cover, is balanced by the interaction between the potential transport capacity (the term in square brackets), described as a function of soil, rainfall, and topography. For the sake of simplification, Equation (2) does not explicitly account for rainfall-runoff processes (with the so called rainfall-runoff erosivity factor, R, calculated as a summation of event-based energy-intensity values, EI30), under consideration that in the USLE some relationship does exist between runoff and EI30, and between the peak rate of runoff rate and the maximum rainfall intensity in 30 min (I30) [47]. The variables between square brackets are explained as follows:
k m    =    K W ( d w e t + C ) + Z
where km is the erodibility index, dwet is the number of rainy days in a month with rain depth ≥1 mm·d−1, K, W, C, and Z are empirical parameters;
E I m    =    0 . 1174 ( p r d r 0.53 h r 1.18 )
where EIm is the rainfall erosivity-factor (MJ·mm·ha−1·h−1·month−1) that provides the forces applied to the soil and causing water erosion (indirect RUSLE procedure; details in [48]). In Equation (4), pr is the monthly precipitation amount (mm), dr is the daily maximum rain in each month (mm), and hr is the hourly maximum rain in each month (mm);
( γ d o m s 1 D )
is the variogram equation to calculate the fractal dimension of topography [49], which can be converted to the equality;
| Z p Z q d o m s |    =    γ d o m s 1 D
where Zp and Zq are the elevations at the locations p and q, respectively. doms is the distance between p and q (that can be viewed as the original modelling scale), γ is a coefficient, and D is the fractal dimension. The following equation:
| Z p Z q d o m s |
represents the surface slope, so it can be assumed that the slope value, S (in Equation (1)), is associated to its corresponding scale–grid size doms [49]. As the USLE was originally developed for agricultural fields, its application to landscape scale erosion modeling is often inappropriate. Developed models based on the unit stream power theory [50,51] include the influence of terrain shape and are more suitable for complex topographic conditions, permitting to accommodate issues of upscaling with soil erosion models. USLE and RUSLE, which consider soil erosion along the flow line (such as in a plot field), have disadvantages because they neglect the influence of flow convergence/divergence, and their application is limited only to areas experiencing net erosion. The general form of the sediment transport equation can be used to describe the effects of terrain and scaling on soil erosion [52]. In particular, the following term in Equation (2):
L ( γ d o m s 1 D )
can be viewed as the RUSLE-slope length factor LS (m). According to Mitasova et al. [51], the unit stream power varies with the runoff contributing area (A, in m2·m−1), and the slope angle β (degrees), as:
L ( γ d o m s 1 D )    =    ( m + 1 ) ( A a o ) m ( sin β b o ) n
where ( γ d o m s 1 D ) is the variogram equation to calculate the fractal dimension of topography in the upscaling equation, ao = 22.13 is the RUSLE unit plot length (m), bo = 0.0896, and m is the slope length exponent equal to (β/(1 + β)) being related to the ratio β of rill erosion (caused by overland flow) to interrill erosion (principally caused by raindrop impact; [51]). Parameters m and n are set equal to 0.5 and 1.3, respectively [51]
The upslope contributing area is calculated from the Digital Elevation Model having a cell size of 20 m, and using the SINMAP tool of ArcView GIS [53,54]. This tool employs the D-∞ technique proposed by Tarboton [55], that computes the contributing area more accurately on divergent hillslopes and uses a multiple-flow direction algorithm allowing for (i) continuous flow angles, and (ii) flow partitioning between one or two neighbor pixels. For the plot-scale based model:
  • the contributing area is equal to the plot area; and,
  • the term e(−ω∙100∙NDVIm) is the Thornes vegetation exponential function, where NDVIm is the Normalized Difference Vegetation Index derived from the red near-infrared reflectance ratio (NDVI = (NIRRED)/(NIR + RED)) (where NIR and RED are the amounts of near-infrared and red light, respectively, reflected by the vegetation and captured by the sensor of the satellite).
NDVI values were estimated indirectly according to [56]. NDVI takes on values of 0 for bare soil, and 0.8 for maximum vegetation greenness. Parameter ω is a function of the ratio β defined above [43]. Exponent η of Equation (2) depends on the type of sediments transported by the erosion process. In similar problems, Carson and Kirkby [57] showed that various suspended sediment and the surface erosion equations correspond to η = 1.0, 1.5, 1.6, and 2.0.

3.3. Model Calibration

The model for predicting monthly gross soil erosion rates was run at MetEROBS over the years from 1986 to 2006. The data gathered over the overlapping period between the modeled time series and the plot measurement (i.e., from 2001 to 2006) have been used to calibrate the model parameters. Model calibration was performed using the measurements collected in the 23 m2 experimental plot of the site. Instead, model validation was carried out by comparing the modeled data with the datasets collected by the Italian Hydrographic and Geographic Service within the Calore River Basin (having a drainage area of about 3000 km2), during the years of 1959 and 1961 [58]. We used only these two years for model validation because of the absence of longer overlapping time-series.
The values of the empirical parameters were determined at the calibration step using mean-monthly measurements through a least-square minimization technique. An iterative calibration process was employed to identify the terms of Equation (2). First, the set of parameters was determined for km Equation (3), fitting the erosion values and keeping constant the parameters of Equation (2). Next, the parameters of Equation (2) were calibrated against the variables measured at the plot scale. The process was reiterated up to reach a converging solution. The parameterization obtained (see Table 2) roughly matched the soil erosion data of both calibration and validation datasets (Figure 5a,b). A set of performance statistics—Percent Relative Root Mean Square Error (RRMSE%, from 0, best, to positive infinite [59], Modelling Efficiency (EF), from negative infinite to 1, best [60], Coefficient of Residual Mass (CRM), from negative infinite, over-estimation, to positive infinite, under-estimation [61])—computed on the basis of the calibration dataset indicated a good agreement between modeled and measured data (RRMSE = 12%; EF = 0.81; CRM = 0.15).

4. Results and Discussion

Equation (2) is used to reconstruct the historical time variations of the monthly gross soil erosion rates along the eastern slope of the Monte Pino site. The bar graph in Figure 6 displays the amount of particles mobilized within the Monte Pino site, whereas the blue lines show the erosivity density as computed by the ratio of simulated erosivity, versus the observed rainfall depth from January 1986 to December 2006. In the graph showing the gross soil erosion time series, one can easily detect two main periods. The first period spans from 1986 to 1994, when the monthly gross erosion values take on values near the median and attain the maximum values close to the average. The second period spans from 1995 to 2006, and is characterized by extreme values of monthly gross erosion abruptly increased with the highest peak located in November 1997. Although such increments match the general trend of variation of the erosivity density for the same period, the November 1997 peak does not correspond to a peak of the erosivity density. This can be interpreted by the fact that during this precipitation event the larger amount of soil particles are mobilized by overland flow.
Gross soil erosion time series is depicted in Figure 6, and characterized by a non-Normal probability distribution with a relatively great standard deviation value (not shown). The difference between the mean and median values of this series indicates also an erratic distribution of erosion-events that, in particular months, might also exceed the tolerable soil loss long-term annual-threshold (about 7 Mg·ha−1; [58]). This was the case of the events that occurred in January, August, and November 1997, May 2000, and 2001, June and September 2003.
It is important to note that our model runs at a monthly timescale and is not able to predict soil erosion rates caused by daily storm events. For instance, more than 70% of sediment produced in November 1997 derived from an erosion process that occurred on 13 November, when the rainfall depth recorded at MetEROBS was 111 mm, of which about 86 mm fell in only three hours. Due to this event, remarkable soil erosion, floods, and also landslides occurred in the entire Benevento Province. During this event, in the Serretelle River Basin (including the MPES) the upper parts of high order channels were extensively incised (laterally and vertically) in response to heavy rainfall (Figure 7a). For this area, the analysis of historical precipitation data, referring to three-hour rainfall (Gumbel method, [62]), shows that a rain event of 86 mm has a return period of 20 years. On 13 November 1997 the rainfall erosivity was 933 MJ·mm·ha−1·h−1, a value comparable with the annual average value recorded for this area. This gives an idea of the high erosive-and overland-flow potential of this storm.
Most of the intense rainstorms with a low spatial extent, including the event of May 2001, occurred between the end of spring and summer. In May 2001, a heavy downpour hit the agricultural lands of the MPES, inducing extreme topographical changes connected with triggered erosion-processes (Figure 7b). In particular, precipitation characterizing this meteorological event amounted to 65 mm in 35 min, with a maximum intensity in 30 min higher than 100 mm·h−1. Field observations revealed important changes in the landscape caused by this event, like deepening of the rill network and development of new rills even in low and flat places. Since extreme rainfall events (e.g., greater than 100 mm·d−1) are not frequent at MetEROBS, the soil material forming the hillslope surface is easily transported along the slope also as consequence of moderate events (rainfall intensity roughly less than 20 mm·h−1), producing changes and exporting the material out of agricultural fields.
Figure 7 illustrates two situations in which remarkable land degradation phenomena at MPES have been caused by extreme rainfall events. On November 1997 and May 2001, which can be considered as representative months with respect to modeling erosion processes under Mediterranean climate conditions, the estimated sediment was close to 50 Mg·ha−1 and 30 Mg·ha−1, respectively. Such values are much more above the long-term monthly average of gross soil erosion of 4 Mg·ha−1 and 2 Mg·ha−1. Based on both field observations and modeling outcomes (i.e., see Figure 7, left) we put forward the hypothesis that during the event occurring on November 1997, soil particles moving within MPES were successively transferred to the floodplain via drainage channels, and finally into Serretelle river. Field observations, on the other hand, indicated that during the extreme event that occurred on 14 May 2001, the soil material mobilized along the Monte Pino slope was prevalently trapped within concave subsections of the area. Our interpretation is that the hydrologic condition of the upper pedogenized layer (i.e., water content), forming the ground surface, influence the persistence and duration of runoff (together with morphology and land-use/cover of the slope). In other words, during the November 1997 event runoff duration was longer than the time of concentration (tc) for this part of the basin, and thus mobilized sediment that reached the river. Conversely, during the May 2001 event, runoff was concentrated within the MPES area, and its duration was shorter than the time of concentration. Thus, most of the mobilized sediment did not reach the river outlet.
On these bases, erosional processes in the MPES can be triggered by meteorological events of magnitude around the average (for the site), but also during low-frequency, high-magnitude short events. This is the case of the extreme meteorological events occurring in May 2001 and June 2003. Field observations indicated that during high-magnitude phenomena, most of the soil lost along the hillslope was deposited before reaching the main stream outlet. Therefore, soil erosion within the MPES was controlled by spring-summer high-intensity rainfall events, and autumn and winter prolonged rainfall inducing persistent runoff.
The seasonal variability in erosivity and gross erosion is very important for evaluating the seasonal hazard of erosive rainstorm events. Figure 8 highlights the seasonal variability that affects the ratio of monthly to annual erosivity and gross cumulative erosion. The relative proportion between the monthly erosivity and erosion ratios remains similar, except for April, August, September, November and December. Especially the months of September and November present a higher disproportionality connected to the ability of heavy rainfall events to trigger erosive processes, and to soil condition and soil use in these months.

Seasonal Timing of the Intensifying Storminess

A more enlarged view of the Mediterranean Central Area (MCA) over the recent decades shows that autumn seasons are prone to intensified precipitation rates, with positive anomalies over many zones of central and southern Italy. Increased temperatures may result, paradoxically, in a generally stationary trend of total annual rainfall, in spite of the fact that daily rainfall increases [63]. In this context, it is relevant to learn how past warming has affected the changes of extreme precipitation [64,65].
To identify which monthly or seasonally rainfall could have led to a major increase in the annual storm erosivity in Southern Italy, we compared the rainfall intensity of the more recent period 1991–2014, with that one of the previous period 1961–1990. As shown in Figure 9, the most positive anomalies of rainfall intensity are observed mainly during September and December.
A positive core was extended from central Italy toward southern Italy in September, which might have triggered a high-velocity of soil loss, and heavy rain in the months in which the higher frequency of soil tillage was detected in the study area.
The observed increasing trend in the rain rates corresponds to the climate change prognosis that predicts reductions in average summer and autumn precipitation, combined with an increase in high-intensity rainy events for many parts of European agricultural lands [67,68]. Nearing et al. [69] point out that just because the sensitivity values for runoff and erosion were generally greater for rainfall changes as opposed to cover changes, this does not imply that future changes in rainfall will dominate over changes in land-use. Predictions by climate change experts suggest that the possibility of rainfall changes in the order of a few percent values in rainfall depths and intensities. In this context, with global and local climate change, an increasing number of ecosystems exist that are experiencing meteorological events of novel magnitudes, with timing and durations that are out of synchrony with the stress thresholds of organismic communities [70,71].

5. Conclusions

We have combined the RUSLE (Revised Universal Soil Loss Equation) and Thornes’ models to simulate gross soil erosion rates at the Met European Research Observatory of Monte Pino experimental site (Southern Italy), using plot-scale measurements and weather observations. The proposed method, resulting from a GIS-based upscaling procedure shows the fractal pattern that mirrors at a larger scale (sloped side surface) the soil loss data collected from plot-scale measurements. Model upscaling through scale-related parameters is thus a key to represent slope dynamics with plot-based, semi-empirical modeling. This approach proved effective to detect magnitude and timing of disproportional and nonlinear variations of soil erosion over a decadal time-scale in the Calore River Basin. Results from our modeling study indicate a clear increase in the peak values of gross erosion after 1995, with the highest peak value occurring in November 1997. Other peaks occurred in May 2001, and June 2003. Field observations revealed that it was mostly during long and moderately intense rainfall events occurring in winter that the sediment mobilized along the sloped side surfaces reached the main river, contributing to the sediment yield. These results thus suggest that, in the study area, the control exerted on soil erosion by spring and summer rainfall events of short duration and high intensity might be a less important part of the hydrological control on sediment yield. As revealed by some discrepancies between basin-wide estimates and measurements, hidden non-linear processes occur in down-and-up-scaling ranges [72], not easy to capture through parsimonious modeling. Presentation of these results is intended to provide an initial understanding for more detailed numerical modeling. For instance, the capacity to predict event soil losses could be improved by considering runoff explicitly [47,73]. This implies reporting discharge information or runoff data from the experimental plot, which were not available in this study.

Author Contributions

Nazzareno Diodato conceived and designed the study, collected data, did basic data analysis, and drafted the paper, figures, and tables; Luigi Guerriero, Gianni Bellocchi supervised the work and discussed the results. All authors actively participated in revising the paper, contributed to the discussion, and gave approval of the final version. Nazzareno Diodato is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fu, B.J.; Zhao, W.; Chen, L.-D.; Lü, Y.-H. Eco-hydrological effects of landscape pattern change. Landsc. Ecol. Eng. 2005, 1, 25–32. [Google Scholar] [CrossRef]
  2. Busari, M.A.; Kukal, S.S.; Kaur, A.; Bhatt, R.; Dulazy, A.A. Conservation tillage impacts on soil, crop and environment. Int. Soil Water Conserv. Res. 2015, 3, 119–129. [Google Scholar] [CrossRef]
  3. Mikha, M.M.; Benjamin, J.G.; Vigil, M.F.; Poss, D.J. Manure and tillage use in remediation of eroded land and impacts on soil chemical properties. PLoS ONE 2017, 12, e0175533. [Google Scholar] [CrossRef] [PubMed]
  4. Diodato, N.; Soriano, M.; Bellocchi, G.; Fiorillo, F.; Cevasco, A.; Revellino, P.; Guadagno, F.M. Historical evolution of slope instability in the Calore River Basin, Southern Italy. Geomorphology 2017, 282, 74–84. [Google Scholar] [CrossRef]
  5. Naylor, L.A.; Spencer, T.; Lane, S.N.; Darby, S.E.; Magilligan, F.J.; Macklin, M.G.; Möller, I. Stormy geomorphology: Geomorphic contributions in an age of climate extremes. Earth Surf. Proc. Landf. 2017, 42, 166–190. [Google Scholar] [CrossRef]
  6. Diodato, N.; Bellocchi, G. Introduction. In Storminess and Environmental Change; Diodato, N., Bellocchi, G., Eds.; Springer: Dordrecht, The Netherlands, 2014; pp. 1–14. [Google Scholar]
  7. Borrelli, P.; Märker, M.; Panagos, P.; Schütt, B. Modeling soil erosion and river sediment yield for an intermountain drainage basin of the Central Apennines, Italy. Catena 2014, 114, 45–58. [Google Scholar] [CrossRef]
  8. Thornes, J. Land degradation. In The Physical Geography of the Mediterranean; Woodward, J., Ed.; Oxford University Press: Oxford, UK, 2009; pp. 563–581. [Google Scholar]
  9. Poesen, J.W.A.; Hooke, M. Erosion, flooding and channel management in Mediterranean environments of southern Europe. Prog. Phys. Geog. 1997, 21, 157–199. [Google Scholar] [CrossRef]
  10. Boukheir, R.; Abdallah, C. Conceptualization of GIS field prediction regional soil-erosion Mediterranean models, case study Lebanon. Geophys. Res. Abstr. 2006, 8, 01289. [Google Scholar]
  11. Sundquist, B. Topsoil Loss and Degradation—Causes, Effects and Implications: A Global Perspective. Available online: http://www.civilizationsfuture.com/bsundquist/se0.html (accessed on 16 August 2017).
  12. Borrelli, P.; Diodato, N.; Panagos, P. Rainfall erosivity in Italy: A national scale spatio-temporal assessment. Int. J. Digit. Earth 2016, 9, 835–850. [Google Scholar] [CrossRef]
  13. Nearing, M.A. Potential changes in rainfall erosivity in the US with climate change during the 21st century. J. Soil Water Conserv. 2001, 56, 220–232. [Google Scholar]
  14. Zhang, X.C.; Nearing, M.A. Impact of climate change on soil erosion, runoff, and wheat productivity in central Oklahoma. Catena 2005, 61, 185–195. [Google Scholar] [CrossRef]
  15. Diodato, N.; Bellocchi, G. Enhanced propagation of rainfall kinetic energy in the UK. Theor. Appl. Climatol. 2017, 129, 1335–1340. [Google Scholar] [CrossRef]
  16. Martínez-Murillo, J.F.; Nadal-Romero, E.; Regüés, D.; Poesen, J. Soil erosion and hydrology of the western Mediterranean badlands throughout rainfall simulation experiments: A review. Catena 2013, 106, 101–112. [Google Scholar] [CrossRef]
  17. Morugán-Coronado, A.; Garciá-Orenes, F.; Cerdiá, A. Effect of land management on soil properties in flood irrigated citrus orchards in Eastern Spain. Soil Discuss. 2015, 2, 1–27. [Google Scholar] [CrossRef]
  18. Aiello, A.; Adamo, M.; Canora, F. Remote sensing and GIS to assess soil erosion with RUSLE3D and USPED at river basin scale in southern Italy. Catena 2015, 131, 174–185. [Google Scholar] [CrossRef]
  19. Conforti, M.; Buttafuoco, G. Interplay between land use changes and soil erosion in a small mountainous catchment in southern Italy. Rendiconti Online Societa Geologica Italiana 2016, 38, 25–28. [Google Scholar] [CrossRef]
  20. Rodolfi, G.; Agnesi, V.; Augelli, P.P.C.; Conoscenti, C.; Del Monte, M.; Fredi, P.; Gentili, B.; Lupia Palmieri, E.; Märker, M.; Materazzi, M.; et al. Soil erosion by water in Mediterranean environment: Italian assessment network of test areas and catchments (EROMED). Geophys. Res. Abstr. 2005, 7, 07527. [Google Scholar]
  21. FAO (Food and Agriculture Organization). Climate Change Adaptation and Mitigation in the Food and Agriculture Sector; Technical Background Document from the Expert Consultation; FAO: Rome, Italy, 2008. [Google Scholar]
  22. Mantel, S.; van Lynden, G.J.; Huting, J. Deliverable 10: Scenario Analysis. Pan-European Soil Erosion Risk Assessment; European Commission, Joint Research Center: Ispra, Italy, 2017. [Google Scholar]
  23. Synder, K.A.; Tartowski, S.L. Multi-scale temporal variation in water availability: Implications for vegetation dynamics in arid and semi-arid ecosystems. J. Arid Environ. 2006, 65, 219–234. [Google Scholar] [CrossRef]
  24. Le Bissonnais, Y.; Montier, C.; Jamagne, M.; Daroussin, J.; King, D. Mapping erosion risk for cultivated soil in France. Catena 2002, 46, 207–220. [Google Scholar] [CrossRef]
  25. Brunsden, D. Back A’long: A Millennial Geomorphology. In Geomorphological Processes and Landscape Change; Higgitt, D.L., Lee, E.M., Eds.; Blackwell Publisher: Oxford, UK, 2001; pp. 27–60. [Google Scholar]
  26. Wainwright, J.; Parsons, A.J.; Abrahams, A.D. Plot-scale studies of vegetation, overland flow and erosion interactions: Case studies from Arizona and New Mexico. Hydrol. Process. 2000, 14, 2921–2943. [Google Scholar] [CrossRef]
  27. Chmelová, R.; Šarapatka, B. Soil erosion by water: Contemporary research methods and their use. Geographica 2002, 37, 23–30. [Google Scholar]
  28. Ollesch, G.; Vacca, A. Influence of time on measurement results of erosion plot studies. Soil Tillage Res. 2002, 67, 23–39. [Google Scholar] [CrossRef]
  29. Cerdan, O.; Le Bissonnais, Y.; Govers, G.; Lecomte, V.; van Oost, K.; Couturier, A.; King, C.; Dubreuil, N. Scale effect on runoff from experimental plots to catchments in agricultural areas in Normandy. J. Hydrol. 2004, 299, 4–14. [Google Scholar] [CrossRef]
  30. Assouline, S.; Ben-Hur, M. Effects of rainfall intensity and slope gradient on the dynamics of interrill erosion during soil surface sealing. Catena 2006, 66, 611–620. [Google Scholar] [CrossRef]
  31. Boix-Fayos, C.; Martinez-Mena, M.; Arnau-Rosalèn, E.; Calvo-Cases, A.; Castillo, V.; Albaladejo, J. Measuring soil erosion by field plots: Understanding the sources of variation. Earth Sci. Rev. 2006, 78, 267–285. [Google Scholar] [CrossRef]
  32. Parsons, A.J.; Brazier, R.E.; Wainwright, J.; Powell, D.M. Scale relationships in hillslope runoff and erosion. Earth Surf. Proc. Landf. 2006, 31, 1384–1393. [Google Scholar] [CrossRef]
  33. Di Stefano, C.; Ferro, V.; Pampalone, V. Applying the USLE family of models at the Sparacia (South Italy) experimental site. Land Degrad. Dev. 2017, 28, 994–1004. [Google Scholar] [CrossRef]
  34. Basso, F.; Pisante, M.; Basso, B. Soil erosion and land degradation. In Mediterranean Desertification; Geeson, N.A., Brandt, C.J., Thornes, J.B., Eds.; Wiley: Chichester, NJ, USA, 2002; pp. 347–359. [Google Scholar]
  35. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses from Cropland East of the Rocky Mountains: Guide for Selection of Practices for Soil and Water Conservation; Agricultural Research Service, US Department of Agriculture in Cooperation with Purdue Agricultural Experiment Station: Washington, DC, USA, 1965; p. 47.
  36. Pan European Soil Erosion Risk Assessment—PESERA. Available online: http://esdac.jrc.ec.europa.eu/public_path/Pesera.pdf (accessed on 16 August 2017).
  37. HyMeX Database. Available online: http://mistrals.sedoo.fr/HyMeX (accessed on 16 August 2017).
  38. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses—A Guide to Conservation Planning; US Department of Agriculture: Washington, DC, USA, 1978; pp. 1–69.
  39. Williams, J.R. Sediment-yield prediction with universal equation using runoff energy factor. In Present and Prospective Technology for Predicting Sediment Yield and Sources; US Department of Agriculture, Agriculture Research Service: Washington, DC, USA, 1975; pp. 244–252. [Google Scholar]
  40. Renard, K.G.; Foster, G.R.; Weesies, G.A.; McCool, D.K.; Yoder, D.C. Predicting Soil Erosion by Water: A Guide to Conservation Planning with the Revised Universal Soil Loss Equation (RUSLE); Agriculture Handbook No. 703; US Department of Agriculture: Washington, DC, USA, 1997; p. 384.
  41. Laflen, J.M.; Lane, L.J.; Foster, G.R. WEPP: A new generation of erosion prediction technology. J. Soil Water Conserv. 1991, 46, 34–38. [Google Scholar]
  42. Morgan, R.P.C.; Quinton, J.N.; Smith, R.E.; Govers, G.; Posen, J.W.A.; Auerswald, K.; Chisci, G.; Torri, D.; Styczen, M.E. The European soil erosion model (EUROSEM): A dynamic approach for predicting sediment transport from fields and small catchments. Earth Surf. Proc. Landf. 1998, 23, 527–544. [Google Scholar] [CrossRef]
  43. Thornes, J.B. The ecology of erosion. Geography 1985, 70, 222–234. [Google Scholar]
  44. Chen, Z.-Q.; Govindaraju, R.S.; Kavvas, M.L. Spatial averaging of unsaturated flow equations under infiltration conditions over areally heterogeneous fields: 1. Development of models. Water Resour. Res. 1994, 30, 523–533. [Google Scholar] [CrossRef]
  45. Mulligan, M.; Wainwright, J. Modelling and model building. In Environmental Modelling; Wainwright, J., Mulligan, M., Eds.; John Wiley and Sons, Ltd.: Chichester, NJ, USA, 2013; pp. 7–73. [Google Scholar]
  46. Kavvas, M.L. On the coarse-graining of hydrologic processes with increasing scales. J. Hydrol. 1999, 217, 191–202. [Google Scholar] [CrossRef]
  47. Kinnell, P.I.A. Event soil loss, runoff and the universal soil loss equation family of models: A review. J. Hydrol. 2010, 385, 384–397. [Google Scholar] [CrossRef]
  48. Diodato, N. Predicting RUSLE (Revised Universal Soil Loss Equation) monthly erosivity index from readily available rainfall data in Mediterranean area. Environmentalist 2005, 25, 63–70. [Google Scholar] [CrossRef]
  49. Zhang, X.; Drake, N.A.; Wainwright, J. Scaling issues in environmental modelling. In Environmental Modelling; Mulligan, M., Wainwright, J., Eds.; Wiley & Sons, Ltd.: Chichester, NJ, USA, 2004; pp. 7–73. [Google Scholar]
  50. Mitasova, H.; Iverson, L. Erosion and sedimentation potential analysis for Hunter lake. In An Environmental Assessment of the Hunter Lake Project Area; Brigham, W.U., Brigham, A.R., Eds.; Illinois Natural History Survey: Champaign, IL, USA, 1992. [Google Scholar]
  51. Mitasova, H.; Hofierka, J.; Zlocha, M.; Iverson, L.R. Modelling topographic potential for erosion and deposition using GIS. Int. J. Geogr. Inf. Sci. 1996, 10, 629–641. [Google Scholar] [CrossRef]
  52. Moore, I.D.; Wilson, J.P. Length-slope factors for the Revised Universal Soil Loss Equation: Simplified method of estimation. J. Soil Water Conserv. 1992, 47, 423–428. [Google Scholar]
  53. Pack, R.T.; Tarboton, D.G.; Goodwin, C.N. Terrain Stability Mapping with SINMAP; Technical Description and Users Guide for Version 1.00, Report Number 4114-0; Terratech Consulting Ltd.: Salmon Arm, BC, Canada, 1998. [Google Scholar]
  54. Pack, R.T.; Tarboton, D.G.; Goodwin, C.N. Assessing terrain stability in a GIS using SINMAP. In Proceedings of the 15th Annual GIS Conference, GIS 2001, Vancouver, BC, Canada, 19–22 February 2001. [Google Scholar]
  55. Tarboton, D.G. A new method for the determination of flow directions and contributing areas in grid digital elevation models. Water Resour. Res. 1997, 33, 309–319. [Google Scholar] [CrossRef]
  56. Diodato, N.; Bellocchi, G. Modelling NDVI responses to climate variability in Mediterranean terrestrial ecosystem. Environ. Monit. Assess. 2008, 44, 147–159. [Google Scholar] [CrossRef] [PubMed]
  57. Carson, M.A.; Kirkby, M.J. Hillslope Form and Process; Cambridge University Press: New York, NY, USA, 1972; p. 475. [Google Scholar]
  58. Diodato, N. Modelling net erosion responses to enviroclimatic changes recorded upon multisecular timescales. Geomorphology 2006, 80, 164–177. [Google Scholar] [CrossRef]
  59. Jørgensen, S.E.; Kamp-Nielsen, L.; Christensen, T.; Windolf-Nielsen, J.; Westergaard, B. Validation of a prognosis based upon a eutrophication model. Ecol. Modell. 1986, 35, 165–182. [Google Scholar] [CrossRef]
  60. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  61. Loague, K.; Green, R.E. Statistical and graphical methods for evaluating solute transport models: Overview and application. J. Contam. Hydrol. 1991, 7, 51–73. [Google Scholar] [CrossRef]
  62. Gumbel, E. Statistics of Extremes; Columbia University Press: New York, NY, USA, 1958; p. 400. [Google Scholar]
  63. Alpert, P.; Ben-Gai, T.; Baharad, A.; Manes, A. The paradoxical increase of Mediterranean extreme daily rainfall in spite of decrease in total value. Geophys. Res. Lett. 2002, 29, 1536–1539. [Google Scholar] [CrossRef]
  64. Klein Tank, A.M.G.; Peterson, T.C.; Quadir, D.A.; Dorji, S.; Zou, X.; Tang, H.; Santhosh, K.; Joshi, U.R.; Jaswal, A.K.; Kolli, R.K.; et al. Changes in daily temperature and precipitation extremes in central and south Asia. J. Geophys. Res. 2006, 111, D16105. [Google Scholar] [CrossRef]
  65. Trenberth, K.E. Changes in precipitation with climate change. Clim. Res. 2011, 47, 123–138. [Google Scholar] [CrossRef]
  66. Earth System Research Laboratory—Physical Science Division. Available online: http://www.cdc.noaa.gov (accessed on 16 August 2017).
  67. Christensen, J.H.; Christensen, O.B. Climate modelling: Severe summertime flooding in Europe. Nature 2003, 421, 805–806. [Google Scholar] [CrossRef] [PubMed]
  68. Meusburger, K.; Steel, A.; Panagos, P.; Montanarella, L.; Alewell, C. Spatial and temporal variability of rainfall erosivity factor for Switzerland. Hydrol. Earth Syst. Sci. 2012, 16, 167–177. [Google Scholar] [CrossRef] [Green Version]
  69. Nearing, M.A.; Jetten, V.; Baffaut, C.; Cerdan, O.; Couturier, A.; Hernandez, M.; Bissonnais, Y.L.; Nichols, M.H.; Nunes, J.P.; Renschler, C.S.; et al. Modeling response of soil erosion and runoff to changes in precipitation and cover. Catena 2005, 61, 131–154. [Google Scholar] [CrossRef]
  70. Reusch, T.B.H.; Ehlers, A.; Hämmerli, A.; Worm, B. Ecosystem recovery after climatic extremes enhanced by genotypic diversity. Proc. Natl. Acad. Sci. USA 2005, 102, 2826–2831. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Jentsch, A.; Beierkuhnlein, C. Research frontiers in climate change: Effects of extreme meteorological events on ecosystems. C. R. Geosci. 2008, 340, 621–628. [Google Scholar] [CrossRef]
  72. Lane, L.J.; Hernandez, M.; Nichols, M. Processes controlling sediment yield from watersheds as functions of spatial scale. Environ. Modell. Softw. 1997, 12, 355–369. [Google Scholar] [CrossRef]
  73. Kinnel, P.I.A.; Risse, M. USLE-M: Empirical modeling rainfall erosion through runoff and sediment concentration. Soil Sci. Soc. Am. J. 1998, 62, 1667–1672. [Google Scholar] [CrossRef]
Figure 1. Soil erosion in Mediterranean area from PESERA (Pan-European Soil Erosion Risk Assessment) project of the European Commission—Joint Research Center, Soil Bureau (arranged from the PESERA Map [36]), with the location of study area (white square).
Figure 1. Soil erosion in Mediterranean area from PESERA (Pan-European Soil Erosion Risk Assessment) project of the European Commission—Joint Research Center, Soil Bureau (arranged from the PESERA Map [36]), with the location of study area (white square).
Environments 04 00058 g001
Figure 2. Geographic map of Benevento Province (a,b) and topography of the southwest Benevento District (c), with relative MPES site (square dotted) and Monte Pino Met European Research Observatory (MetEROBS).
Figure 2. Geographic map of Benevento Province (a,b) and topography of the southwest Benevento District (c), with relative MPES site (square dotted) and Monte Pino Met European Research Observatory (MetEROBS).
Environments 04 00058 g002
Figure 3. Spatial pattern of the annual average precipitation around the Monte Pino experimental site (little square; (a)) and orthophoto of the hillslope large area in spring season (b).
Figure 3. Spatial pattern of the annual average precipitation around the Monte Pino experimental site (little square; (a)) and orthophoto of the hillslope large area in spring season (b).
Environments 04 00058 g003
Figure 4. Views of (a) hillshade-map and (b) slope-map at 50 m grid-cell resolution of the Monte Pino experimental site (MPES) with drainage network and MetEROBS (a,b), and experimental plot placed on the eastern slope of the MetEROBS (c).
Figure 4. Views of (a) hillshade-map and (b) slope-map at 50 m grid-cell resolution of the Monte Pino experimental site (MPES) with drainage network and MetEROBS (a,b), and experimental plot placed on the eastern slope of the MetEROBS (c).
Environments 04 00058 g004
Figure 5. Scatterplots of measured and estimated monthly soil loss values at the plot scale of Monte Pino on the calibration data set (a), and between estimated and measured sediment at the scale of Calore River Basin, on the validation data set (b).
Figure 5. Scatterplots of measured and estimated monthly soil loss values at the plot scale of Monte Pino on the calibration data set (a), and between estimated and measured sediment at the scale of Calore River Basin, on the validation data set (b).
Environments 04 00058 g005
Figure 6. Reconstructed monthly gross soil erosion rate (bar) and erosivity density (>1 MJ·ha−1·h−1) averaged on the running of three months (line) from Monte Pino hillslopes large area during the period January 1986–December 2006. The two horizontal black lines depict the gross erosion mean and median values, respectively.
Figure 6. Reconstructed monthly gross soil erosion rate (bar) and erosivity density (>1 MJ·ha−1·h−1) averaged on the running of three months (line) from Monte Pino hillslopes large area during the period January 1986–December 2006. The two horizontal black lines depict the gross erosion mean and median values, respectively.
Environments 04 00058 g006
Figure 7. View of the riverbank erosion along the Serretelle river after the intensive rainstorm occurred on 13 November 1997 (a), and gully erosion in tilled vineyard surrounding Monte Pino site after the downpour of 14 May 2001 (b).
Figure 7. View of the riverbank erosion along the Serretelle river after the intensive rainstorm occurred on 13 November 1997 (a), and gully erosion in tilled vineyard surrounding Monte Pino site after the downpour of 14 May 2001 (b).
Environments 04 00058 g007
Figure 8. Ratio between median monthly erosivity (EI30mon; blue bars) and net erosion (NEmon; brown bars) and the respective annual cumulative values for Monte Pino hillslopes area.
Figure 8. Ratio between median monthly erosivity (EI30mon; blue bars) and net erosion (NEmon; brown bars) and the respective annual cumulative values for Monte Pino hillslopes area.
Environments 04 00058 g008
Figure 9. Monthly rainfall intensity anomalies in 1991–2014, compared with the period 1961–1990 in a Mediterranean Central Area on September (a) and November (b). The little square indicates the study area (Arranged from National Centers for Environmental Prediction reanalysis, NOAA-ESRL [66]).
Figure 9. Monthly rainfall intensity anomalies in 1991–2014, compared with the period 1961–1990 in a Mediterranean Central Area on September (a) and November (b). The little square indicates the study area (Arranged from National Centers for Environmental Prediction reanalysis, NOAA-ESRL [66]).
Environments 04 00058 g009
Table 1. Summary of rainfall and erosion rates for each rainstorm event between May 2001 and November 2006 recorded at Monte Pino—MetEROBS.
Table 1. Summary of rainfall and erosion rates for each rainstorm event between May 2001 and November 2006 recorded at Monte Pino—MetEROBS.
YearMonthDayP (mm·day−1)EI30 (MJ·mm·ha−1·h−1)Erosion (kg·m−2)Events
2001May2367.0230020.000Downpours with flash-floods and surface landslides
2426.62001.740Thunderstorm
August2217.01770.668Thunderstorm with isolated downpours
September1514.8330.056Thunderstorms and showers
1710.0500.160Intense rainfall
November1350.02610.800Intense rainfall with wind gust
2818.0400.223Rainfall
2002April441.0470.044Intense rainfall
May2013.0880.018Thunderstorms and showers
278.6470.031Short thunderstorm
June217.2650.013Thunderstorms and showers
September2266.04940.578Continuous thunderstorms and showers with wind
2326.61130.036Thunderstorms and showers (afternoon)
October1046.41710.033Showers at morning and thunderstorm at evening
2003April626.6200.016Continuous and moderate rain, somewhat with showers
June360.015002.990Afternoon deluge with quiet rain, after stormy
July717.0800.044Moderate rain with isolated thunderstorms in the afternoon
3127.21720.333Heavy thunderstorm at evening
October529.61880.261Heavy thunderstorm followed by moderate rain at nighttime
2927.81360.036Thunderstorms in the evening
2005August2134.0860.267Afternoon thunderstorm
September114.0800.220Afternoon thunderstorm
18–2032.7600.600Showers
2006July2518.21860.187Evening thunderstorm
August610.8380.049Afternoon thunderstorm
820.01480.190Afternoon thunderstorm
1338.23241.320Afternoon thunderstorm
September1534.82180.222Showers
2639.62770.266Thunderstorm and showers
November2232.5450.022Showers and isolated thunderstorms
P: Precipitation; EI30: Rainfall erosivity.
Table 2. Calibrated parameters for the study area.
Table 2. Calibrated parameters for the study area.
ParameterValues
η1.66
K0.001
W0.16
C5
Z1
ω0.05

Share and Cite

MDPI and ACS Style

Diodato, N.; Guerriero, L.; Bellocchi, G. Modeling and Upscaling Plot-Scale Soil Erosion under Mediterranean Climate Variability. Environments 2017, 4, 58. https://0-doi-org.brum.beds.ac.uk/10.3390/environments4030058

AMA Style

Diodato N, Guerriero L, Bellocchi G. Modeling and Upscaling Plot-Scale Soil Erosion under Mediterranean Climate Variability. Environments. 2017; 4(3):58. https://0-doi-org.brum.beds.ac.uk/10.3390/environments4030058

Chicago/Turabian Style

Diodato, Nazzareno, Luigi Guerriero, and Gianni Bellocchi. 2017. "Modeling and Upscaling Plot-Scale Soil Erosion under Mediterranean Climate Variability" Environments 4, no. 3: 58. https://0-doi-org.brum.beds.ac.uk/10.3390/environments4030058

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

Article Metrics

Back to TopTop