Testing slope effect on flow resistance equation for mobile bed rills

In this paper, a recently theoretically deduced rill flow resistance equation, based on a power‐velocity profile, is tested experimentally on plots of varying slopes in which mobile bed rills are incised. Initially, measurements of flow velocity, water depth, cross‐sectional area, wetted perimeter and bed slope conducted in 106 reaches of rills incised on an experimental plot having a slope of 14% were used to calibrate the flow resistance equation. Then, the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number, which was calibrated using the 106 rill reach data, was tested using measurements carried out in plots having slopes of 22% and 9%. The measurements carried out in the latter slope conditions confirmed that (a) the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach, and (b) the data were supportive of the slope independence hypothesis of rill velocity stated by Govers.     

Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

Alternative analysis of transient infiltration experiment to estimate soil water repellency

The repellency index (RI) defined as the adjusted ratio between soil‐ethanol, S_e, and soil‐water, S_w, sorptivities estimated from minidisk infiltrometer experiments has been used instead of the widely used water drop penetration time and molarity of ethanol drop tests to assess soil water repellency. However, sorptivity calculated by the usual early‐time infiltration equation may be overestimated as the effects of gravity and lateral capillary are neglected. With the aim to establish the best applicative procedure to assess RI, different approaches to estimate S_e and S_w were compared that make use of both the early‐time infiltration equation (namely, the 1 min, S1, and the short‐time linearization approaches), and the two‐term axisymmetric infiltration equation, valid for early to intermediate times (namely, the cumulative linearization and differentiated linearization approaches). The dataset included 85 minidisk infiltrometer tests conducted in three sites in Italy and Spain under different vegetation habitats (forest of Pinus pinaster and Pinus halepensis, burned pine forest, and annual grasses), soil horizons (organic and mineral), postfire treatments, and initial soil water contents. The S1 approach was inapplicable in 42% of experiments as water infiltration did not start in the first minute. The short‐time linearization approach yielded a systematic overestimation of S_e and S_w that resulted in an overestimation of RI by a factor of 1.57 and 1.23 as compared with the cumulative linearization and differentiated linearization approaches. A new repellency index, RI_s, was proposed as the ratio between the slopes of the linearized data for the wettable and hydrophobic stages obtained by a single water infiltration test. For the experimental conditions considered, RI_s was significantly correlated with RI and WDPT. Compared with RI, RI_s includes information on both soil sorptivity and hydraulic conductivity and, therefore, it can be considered more physically linked to the hydrological processes affected by soil water repellency.     

δ18O and chemical composition of Libyan Desert Glass,country rocks,and sands: New considerations on target material

Oxygen isotope and chemical measurements were carried out on 25 samples of Libyan Desert Glass (LDG), 21 samples of sandstone, and 3 of sand from the same area. The δ¹⁸O of LDG samples range from 9.0‰ to 11.9‰ (Vienna Standard Mean Ocean Water [VSMOW]); some correlations between isotope data and typological features of the LDG samples are pointed out. The initial δ¹⁸O of a bulk parent material may be slightly increased by fusion due to the loss of isotopically light pore water with no isotope exchange with oxygen containing minerals. Accordingly, the δ¹⁸O of the bulk parent material of LDG may have been about 9.0 ± 1‰ (VSMOW). The measured bulk sandstone and sand samples have δ¹⁸O values ranging from 12.6‰ to 19.5‰ and are consequently ruled out as parent materials, matching the results of previous studies. However, separated quartz fractions have δ¹⁸O values compatible with the LDG values suggesting that the modern surface sand inherited quartz from the target material. This hypothesis fits previous findings of lechatelierite and baddeleyite in these materials. As the age of the parent material reported in previous studies is Pan‐African, we measured the δ¹⁸O values of bulk rock and quartz from intrusives of Pan‐African age and the results obtained were compatible with the LDG values. The main element abundances (Fe, Mg, Ca, K, Na) in our LDG samples conform to previous estimates; Fe, Mg, and K tend to be higher in heterogeneous samples with dark layers. The hypothesis of a low‐altitude airburst involving silica‐rich surface materials deriving from weathered intrusives of Pan‐African age, partially melted and blown over a huge surface by supersonic winds matches the results obtained.     

Petrological, geochemical and isotopic characteristics of the lithospheric mantle beneath Sardinia (Italy) as indicated by ultramafic xenoliths enclosed in alkaline lavas

Mantle xenoliths hosted in Miocene-Quaternary mafic alkaline volcanic rocks from Sardinia have been investigated with electron microprobe, laser ablation microprobe-inductively coupled plasma-mass spectrometry and thermal ionization mass spectrometry techniques. The xenoliths are anhydrous clinopyroxene-poor lherzolites and harzburgites, plus very rare websterites and olivine-websterites. Glassy pods having thin subhedral to euhedral microlites of olivine, clinopyroxene and spinel have been found in harzburgites and websterites. Clinopyroxene shows trace element variability, with values of (La/Yb)_N ranging from sub-chondritic (0.01) to supra-chondritic (8.6). The Sr–Nd isotopic ratios of the clinopyroxenes fall mostly in the field of the European lithospheric mantle xenoliths (⁸⁷Sr/⁸⁶Sr from 0.70385 to 0.70568 and ¹⁴³Nd/¹⁴⁴Nd ranging from 0.512557 to 0.512953). The geochemical characteristics of the Sardinian xenoliths testify to the variable degrees of earlier partial melt extraction, followed by metasomatic modification by alkaline melts or fluids. Websterites are considered to represent small lenses or veins of cumulitic (i.e. magmatic) origin within the mantle peridotite.     

Evolutionary polynomial regression to alert rainfall-triggered landslide reactivation

The derivation of an alert model for landslide risk management is a paramount problem for those sites which are affected by complex landslides involving strategic infrastructures as well as towns. This is a quite common scenario all over the world and then it is a primary problem for the management of geomorphological risk. Along the Adriatic Coast of south Italy, Petacciato landslide is peculiar, since it showed 11 reactivations between 1924 and 2009. It is a deep-seated landslide, and the history of its reactivations shows that even if generally related to quite abundant rainfall periods, there is no clear correlation between rainfall events and reactivations. For this reason, here, an analysis based on a data-driven evolutionary modeling technique is attempted, in order to identify an alert model based on cumulative rainfall heights. Modeling results are quite interesting and encouraging, since they are able to provide landslide forecasting whereas no false positive are ever returned. This work shows the results of this attempt as well as an analysis of the input to the modeling approach, in order to identify which are those cumulative rainfall heights which are physically sound with respect to the particular landslide.     

Inferring groundwater system dynamics from hydrological time-series data

Abstract

The problem of identifying and reproducing the hydrological behaviour of groundwater systems can often be set in terms of ordinary differential equations relating the inputs and outputs of their physical components under simplifying assumptions. Conceptual linear and nonlinear models described as ordinary differential equations are widely used in hydrology and can be found in several studies. Groundwater systems can be described conceptually as an interlinked reservoir model structured as a series of nonlinear tanks, so that the groundwater table can be schematized as the water level in one of the interconnected tanks. In this work, we propose a methodology for inferring the dynamics of a groundwater system response to rainfall, based on recorded time series data. The use of evolutionary techniques to infer differential equations from data in order to obtain their intrinsic phenomenological dynamics has been investigated recently by a few authors and is referred to as evolutionary modelling. A strategy named Evolutionary Polynomial Regression (EPR) has been applied to a real hydrogeological system, the shallow unconfined aquifer of Brindisi, southern Italy, for which 528 recorded monthly data over a 44-year period are available. The EPR returns a set of non-dominated models, as ordinary differential equations, reproducing the system dynamics. The choice of the representative model can be made both on the basis of its performance against a test data set and based on its incorporation of terms that actually entail physical meaning with respect to the conceptualization of the system.

Citation Doglioni, A., Mancarella, D., Simeone, V. & Giustolisi, O. (2010) Inferring groundwater system dynamics from hydrological time-series data. Hydrol. Sci. J. 55(4), 593–608.     

Flooding scenarios due to land subsidence and sea‐level rise: a case study for Lipari Island (Italy)

Archaeological and instrumental data indicate that the southern sector of the volcanic island of Lipari has been subsiding for the last 2100 years due to isostatic and tectonic factors, at variable rates of up to ~11 mm a^?1. Based on this data, a detailed marine flooding scenario for 2100 AD is provided for the bay of Marina Lunga in the eastern part of the island from (1) an ultra‐high‐resolution Digital Terrain and Marine Model (DTMM) generated from multibeam bathymetry (MB) and Unmanned Aerial Vehicles (UAV), (2) the rate of land subsidence from Global Positioning System (GPS) data and (3) the regional sea‐level projections of the International Panel on Climate Change (IPCC). When land subsidence is considered, the upper bound of sea‐level rise is estimated at 1.36 m and 1.60 m for RCP4.5 and RCP8.5 climate change scenarios, respectively. Here, we show the expected impact of marine flooding at Lipari for the next 85 years and discuss the hazard implications for the population living along the shore.