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.     

An automatic approach for rill network extraction to measure rill erosion by terrestrial and low‐cost unmanned aerial vehicle photogrammetry

For an erosion event (October 2016) occurred at the Sparacia experimental area (Southern Italy), both terrestrial and low‐altitude aerial surveys were carried out by consumer grade camera and quadcopter (low‐cost unmanned aerial vehicle [UAV]) to measure rill erosion on two plots with steepness of 22% and 26%. Applying the structure from motion (SfM) technique, the three‐dimensional digital terrain models (3D‐DTMs) and the quasi three‐dimensional models (2.5D‐digital elevation model [DEM]) were obtained by the two surveys. Furthermore, 3D‐DTM and DEM were built using the available aerial photographs (166) and adding 40 terrestrial photographs. For the first time, the convergence index was applied to high‐resolution rill data for extracting the rill network, and a subsequent separation into contributing and non‐contributing rills was carried out. The comparison among the three surveys (terrestrial, UAV, and UAV + terrestrial) was developed using two morphometric parameters of the rill network (drainage density and drainage frequency). Moreover, using as reference the weight of sediment stored on the tanks located downstream of the plots, the reliability of soil loss measurement by 3D models was tested. For both contributing and non‐contributing rills, the morphometric parameters were higher for the terrestrial than for UAV and UAV + terrestrial surveys. For both plots, SfM always provided reliable soil loss measurements, which were affected by errors ranging from ?8% to 13%. Although the applied technique used a low‐cost UAV and a consumer grade camera, the obtained results demonstrated that a reliable estimate of rill erosion can be obtained in an area of interest.     

Multiparametric Approach in Investigating Volcano-Hydrothermal Systems: the Case Study of Vulcano (Aeolian Islands, Italy)

Seismic activity, ground deformation, and soil and fumarole temperatures acquired during 2004–2007 at Vulcano (Aeolian Islands) are analysed and the time relations among the different time series are discussed. Changes in temperature of fumarolic gases took place during four “anomalous” periods (November 2004–March 2005; October 2005–February 2006; August–October 2006; July–December 2007) at the same time as an increasing number of volcano-seismic events. In particular, the temperatures at high temperature vents and at steam heated soil ranged in time from 180 to 440°C and from 20 to 90°C, respectively. The maximum daily number of volcano-seismic events was 57, reached during the second anomalous period. This seismicity, characterised by focal depth generally lower than 1?km below sea level (b.s.l.) and composed of different kinds of events associated to both resonance and shear failure processes, is related to the shallow dynamics of the hydrothermal system. During the analysed period, very few volcano-tectonic earthquakes took place and tilt recordings showed no sharp or important changes. In light of such observations, the increases in both temperature and volcano-seismic events number were associated to increases in the release of gas from a deep and stable magma body, without magma intrusions within the shallow hydrothermal system. Indeed, a greater release of gas from depth leads to increased fluid circulation, that can promote increases in volcano-seismic events number by both fracturing processes and resonance and vibration in cracks and conduits. The different trends observed in the measured geochemical and geophysical series during the anomalous periods can be due to either time changes in the medium permeability or a changing speed of gas release from a deep magma body. Finally, all the observed variations, together with the changing temporal distribution of the different seismic event kinds, suggest that the hydrothermal system at Vulcano can be considered unsteady and dynamic.     

Shape factors of cylindrical piezometers in uniform soil

This article presents the analytical solution of Laplace equation for the steady flow around open-ended cylindrical piezometers located in an infinite, isotropic, and incompressible saturated soil. Shape factors are obtained for piezometers with varying length-to-diameter ratios. Comparisons are made with published factors obtained by means of approximate analytical solutions, numerical approaches, and measurements in electric analog models. It is shown that some expressions that are currently used in practice are inadequate and should be abandoned.     

Structural setting of the Bay of Naples (Italy) seismic reflection data: implications for Campanian volcanism

This paper focuses on the recent tectonic evolution of the Bay of Naples with the aim of exploring the connection between local tectonics and volcanism. We reprocessed the seismic reflection dataset acquired in the area in the late 1973. The new processing was highly successful in obtaining a decisive strong reduction of random noise, removal of coherent noise and reduction of spatial aliasing. Classical interpretative schemes and complex attributes of seismic traces were used to reconstruct fault kinematics and reflector patterns. The results show that the faults affecting the Bay of Naples exhibit prevailing NE structural strikes, with the exception of the Pozzuoli Caldera where NW patterns are also common. Many faults are subvertical and show seismic evidence of volcanic activity along them. A main alignment of conjugate NE–SW faults, named here as “Magnaghi–Sebeto line”, intersects several submarine volcanic banks and separates the bay into two sectors, characterized by important geological, geophysical and petrochemical differences. The structural configuration of the bay may reflect the occurrence of either oblique extension or a transfer zone of the NW–SE fault system, along which, in the Campanian–Lucanian Apennine chain, great vertical displacements occur.