A new version of the USLE‐MM for predicting bare plot soil loss at the Sparacia (South Italy) experimental site

Improving empirical prediction of plot soil erosion at the event temporal scale has both scientific and practical importance. In this investigation, 492 runoff and soil loss data from plots of different lengths, λ (11 ≤ λ ≤ 44 m), and steepness, s (14.9 ≤ s ≤ 26.0%), established at the Sparacia experimental station, in Sicily, South Italy, were used to derive a new version of Universal Soil Loss Equation (USLE)‐MM model, by only assuming a value of one for the topographic length, L, and steepness, S, factors for λ = 22 m and s = 9%, respectively. An erosivity index equal to (Q_REI₃₀)^b1, Q_R and EI₃₀ being the runoff coefficient and the event rainfall erosivity index, respectively, with b₁ > 1 was found to be an appropriate choice for the Sparacia area. The specifically developed functions for L and S did not differ appreciably from other, more widely accepted relationships (maximum differences by a factor of 1.22 for L and 1.09 for S). The new version of the USLE‐MM performed particularly well for highly erosive events, because predicted soil loss differed by not more than a factor of 1.19 from the measured soil loss for measured values of more than 100 Mg ha^?1. The choice of the relationships to predict topographic effects on plot soil loss should not represent a point of particular concern in the application of the USLE‐MM in other environments. However, tests of the empirical approach should be carried out in other experimental areas in an attempt to develop analytical tools, usable at the event temporal scale, reasonably simple and of wide validity. Copyright © 2015 John Wiley & Sons, Ltd.     

Measuring rill erosion at plot scale by a drone‐based technology

The traditional direct method (i.e. metric ruler and rillmeter) of monitoring rill erosion at plot scale is time consuming and invasive because it modifies the surface of the rilled area. Measuring rill features using a drone‐based technology is considered a non‐invasive method allowing a fast field relief. In the experimental Sparacia area a survey by a quadricopter Microdrones md4‐200 was carried out, and this relief allowed the generation of a Digital Elevation Model (DEM), with a mesh size of 1 cm and a resolution elevation equal to 2 mm, for three plots (L, G and C) affected by rill erosion. At first for the experimental L plot, which is 44 m long, the rill features were surveyed by a ‘manual’ method which was carried out by drawing on the PC screen the rill paths obtained by a visual orthophoto interpretation. This manual method was not applicable for the plots in which rills of limited depth occurred and were not detectable. Then, for both L plot and the other experimental plots having a length ranging from 22 to 44 m, an ‘automatic’ extraction method of rills from DEM was applied. Using an appropriate calculation routine, a vector coverage of transects orthogonal to the main flow direction (i.e. the maximum slope steepness path) was generated. The intersection of each plot DEM with the transect coverage allowed to obtain both the cross sections and the main rill morphological features. For the L plot the comparison between the rill morphological features obtained by the two different methods (manual, automatic) was carried out. Finally, the length–volume relationship and a dimensionless relationship proposed in literature were tested for all studied experimental plots. Copyright © 2015 John Wiley & Sons, Ltd.     

The volcanic and geothermally active Campi Flegrei caldera: an integrated multidisciplinary image of its buried structure

The Campi Flegrei caldera in southern Italy is one of the greatest geohazard areas on Earth. Evidence of an active magmatic and geothermal system is provided by ongoing ground uplift, with volcano-tectonic and long-period (LP) seismicity, the persistent degassing of ~1500 tonnes of CO₂ per day, the presence of hot fumaroles at temperatures of 90–150 °C, brine-rich aquifers (with total dissolved solids up to 33 g l^?1) and high thermal gradients in the crust (with temperatures reaching 420 °C at 3,050 m b.s.l.). Since the 1940s, more than 100 exploratory boreholes have been drilled in the area to depths of 80–3,100 m by the Azienda Geologica Italiana Petroli (AGIP) and the Società Anonima Forze Endogene Napoletane (SAFEN). To date, however, no systematic reanalysis of the drilling data has been carried out, and the buried volcanic structure has not been updated using the most recent scientific results and previous findings. By integrating unpublished data from the AGIP and SAFEN reports with published information from geological, volcanological, petrological, petrophysical and geophysical studies, this paper presents an improved picture of the Campi Flegrei caldera that will be useful for volcanic hazard assessment and mitigation in the Naples area and for future research planning. The results suggest that intra-caldera activity has been influenced by how the magmatic system at depths greater than about 4 km has determined the transfer of magma, volatiles, and heat to the overlying geothermal system and, ultimately, to the surface. In particular, intriguing is that the most volcanically active central-eastern sector of the caldera, which is subject to intense bradyseismic ground movement and gas emission, coincides with a structurally delimited subsurface rock volume characterized by an uprising of the 100 °C isotherm, a deep water supply to the shallower aquifer, the early disappearance of secondary calcite, LP seismicity and high seismic S-wave attenuation. In this area, we also document evidence of repeated injection at depths of c. 1.5–3.0 km of isolated and small-volume batches of magma, where occurred their crystallization and degassing. Shallow intrusions and degassing of magma are thus identified as two of the key processes that drive unrest in Campi Flegrei.     

Precipitable water evaluations from infrared sun‐photometric measurements analyzed using the atmospheric hygrometry technique

Measurements of direct solar irradiance were taken employing 4 different sun‐photometers at near infrared wavelengths, suitable for use in atmospheric hygrometry. This technique utilising a set of spectral ratios, in and out of selected water vapour absorption bands, was applied to the measurements to obtain accurate evaluations of precipitable water. For all the hygrometric ratios given by the 4 sun‐photometers used at the 3 stations of Sagres, Monchique and Mt. Foia, during the CLEARCOLUMN experiment, we determined the calibration curves by correcting them for the Rayleigh scattering effects and, then, plotting the natural logarithms of such corrected ratios versus the square root of the water vapour mass present along the atmospheric slant path. The regression lines drawn for the various scatter diagrams were estimated to give evaluations of precipitable water with an uncertainty of less than 5%, 3% and 10% at the 3 stations, respectively. The calibration curves of the sun‐photometer located at the Sagres station were determined using the precipitable water evaluations obtained from the local radiosounding measurements taken on 5 clear‐sky days. Those of the sun‐photometers used at the Monchique and Mt. Foia stations were instead determined through intercomparison between subsets of measurements simultaneously taken with various instruments at Sagres and Mt. Foia. Using these calibration curves, we examined all the field measurements determining the time‐patterns of precipitable water at the 3 stations. During the period from 16 June to 25 July 1997, precipitable water was found to vary between 1.1 and 3.7 g cm⁻² at the Sagres station (with an accuracy within ±13%), between 1.0 and 2.8 g cm⁻² at Monchique (±11%) and between 0.8 and 3.0 g cm⁻² at the top of Mt. Foia (±26%).     

Laboratory Studies on Silicates Relevant for the Physics of TNOs

Silicates are one of the principal components present in Solar System objects.Silicates evolve in space modifying their physical properties according to theastronomical environments they go through. To characterise the nature of TNOsin the framework of the formation and evolution of the Solar System, experimentson structural transitions of silicates have been performed in the laboratoryto simulate some of the processing suffered by the dust. The infrared spectralproperties of possible silicate candidates thought to be present in TNOs have beenstudied. The results of thermal annealing of amorphous silicates and amorphisationof crystalline forsterite (pure-Mg olivine) by ion irradiation are presented. Theobservable properties of TNOs surfaces are inferred.     

Comparing theoretically supported rainfall‐runoff erosivity factors at the Sparacia (South Italy) experimental site

Interpreting rainfall‐runoff erosivity by a process‐oriented scheme allows to conjugate the physical approach to soil loss estimate with the empirical one. Including the effect of runoff in the model permits to distinguish between detachment and transport in the soil erosion process. In this paper, at first, a general definition of the rainfall‐runoff erosivity factor REF_e including the power of both event runoff coefficient Q_R and event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) is proposed. The REF_e factor is applicable to all USLE‐based models (USLE, Modified USLE [USLE‐M] and Modified USLE‐M [USLE‐MM]) and it allows to distinguish between purely empirical models (e.g., Modified USLE‐M [USLE‐MM]) and those supported by applying theoretical dimensional analysis and self‐similarity to Wischmeier and Smith scheme. This last model category includes USLE, USLE‐M, and a new model, named USLE‐M based (USLE‐MB), that uses a rainfall‐runoff erosivity factor in which a power of runoff coefficient multiplies EI₃₀. Using the database of Sparacia experimental site, the USLE‐MB is parameterized and a comparison with soil loss data is carried out. The developed analysis shows that USLE‐MB (characterized by a Nash–Sutcliffe Efficiency Index NSEI equal to 0.73 and a root mean square error RMSE = 11.7 Mg ha^?1) has very similar soil loss estimate performances as compared with the USLE‐M (NSEI = 0.72 and RMSE = 12.0 Mg ha^?1). However, the USLE‐MB yields a maximum discrepancy factor between predicted and measured soil loss values (176) that is much lower than that of USLE‐M (291). In conclusion, the USLE‐MB should be preferred in the context of theoretically supported USLE type models.     

Testing the Universal Soil Loss Equation‐MB equation in plots in Central and South Italy

Planning soil conservation strategies requires predictive techniques at event scale because a large percentage of soil loss over a long‐time period is due to relatively few large storms. Considering runoff is expected to improve soil loss predictions and allows relation of the process‐oriented approach with the empirical one, furthermore, the effects of detachment and transport on soil erosion processes can be distinguished by a runoff component. In this paper, the empirical model USLE‐MB (USLE‐M based), including a rainfall‐runoff erosivity factor in which the event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) multiplies the runoff coefficient Q_R raised to an exponent b₁ > 1 is tested by the measurements carried out for the Masse (10 plots) and Sparacia (22 plots) experimental stations in Italy. For the Masse experimental station, an exponent b₁ > 1 was also estimated by tests carried out by a nozzle‐type rainfall simulator. For each experimental site in fallow conditions, the effect of the sample size of the plot soil loss measurements on the estimate of the b₁ coefficient was also studied by the extraction of a fixed number N of randomly obtained pairs of the normalized soil loss and runoff coefficient. The analysis showed that the variability of b₁ with N is low and that 350 pairs are sufficient to obtain a stable estimate of b₁. A total of 1,262 soil loss data were used to parameterize the model both locally and considering the two sites simultaneously. The b₁ exponent varied between the two sites (1.298–1.520), but using a common exponent (1.386) was possible. Using a common b₁ exponent for the two experimental areas increases the practical interest for the model and allows the estimation of a baseline component of the soil erodibility factor, which is representative of the at‐site soil intrinsic and quasi‐static properties. Development of a single USLE‐MB model appears possible, and sampling other sites is advisable to develop a single USLE‐MB model for general use.     

A comprehensive analysis of Universal Soil Loss Equation-based models at the Sparacia experimental area

Improving Universal Soil Loss Equation (USLE)-based models has large interest because simple and reliable analytical tools are necessary in the perspective of a sustainable land management. At first, in this paper, a general definition of the event rainfall- runoff erosivity factor for the USLE-based models, REF_e = (Q_R)^b1(EI₃₀)^b2, in which Q_R is the event runoff coefficient, EI₃₀ is the single-storm erosion index, and b₁ and b₂ are coefficients, was introduced. The rainfall-runoff erosivity factors of the USLE (b₁ = 0 and b₂ = 1), USLE-M (b₁ = b₂ = 1), USLE-MB (b₁ ≠ 1 and b₂ = 1), USLE-MR (b₁ = 1 and b₂ ≠ 1), USLE-MM (b₁ = b₂ ≠ 1), and USLE-M2 (b₁ ≠ b₂ ≠ 1) can be defined using REF_e. Then the different expressions of REF_e were simultaneously tested against a data set of normalized bare plot soil losses, A_eN, collected at the Sparacia (south Italy) site. As expected, the poorest A_eN predictions were obtained with the USLE. The observed tendency of this model to overestimate small A_eN values and underestimate high A_eN values was reduced by introducing in the soil loss prediction model both Q_R and an exponent for the erosivity term. The fitting to the data was poor with the USLE-MR as compared with the USLE-MB and the USLE-MM. Estimating two distinct exponents (USLE-M2) instead of a single exponent (USLE-MB, USLE-MR, and USLE-MM) did not appreciably improve soil loss prediction. The USLE-MB and the USLE-MM were recognized to be the best performing models among the possible alternatives, and they performed similarly with reference to both the complete data set and different sub-data sets, only including small, intermediate, and severe erosion events. In conclusion, including the runoff coefficient in the soil loss prediction model is important to improve the quality of the predictions, but a great importance has to be paid to the mathematical structure of the model.     

Comparing flow resistance law for fixed and mobile bed rills

Rills caused by run‐off concentration on erodible hillslopes have very irregular profiles and cross‐section shapes. Rill erosion directly depends on the hydraulics of flow in the rills, which may differ greatly from hydraulics of flow in larger and regular channels. In this paper, a recently theoretically deduced rill flow resistance equation, based on a power–velocity profile, was tested experimentally on plots of varying slopes (ranging from 9% to 26%) in which mobile and fixed bed rills were incised. Initially, measurements of flow velocity, water depth, cross‐section area, wetted perimeter, and bed slope, carried out in 320 reaches of mobile bed rills and in 165 reaches of fixed rills, were used for calibrating the theoretical flow resistance equation. Then the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number was separately calibrated for the mobile bed rills and for the fixed ones. The measurements carried out in both conditions (fixed and mobile bed rills) confirmed that the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach. For mobile bed rills, the data were supportive of the slope independence hypothesis of velocity, due to the feedback mechanism, stated by Govers. The feedback mechanism was able to produce quasicritical flow conditions. For fixed bed rills, obtained by fixing the rill channel, by a glue, at the end of the experimental run with a mobile bed rill, the slope independence of the flow velocity measurements was also detected. Therefore, an experimental run carried out by a rill bed fixed after modelling flow action is useful to detect the feedback mechanism. Finally, the analysis showed that, for the investigated conditions, the effect of sediment transport on the flow resistance law can be considered negligible respect to the grain roughness effect.