Quantitative basin characterisation to refine debris-flow triggering criteria and processes: an example from the Italian Western Alps

This paper investigates rapid channelized debris flow related to rainfalls in small alpine basins. Its goal is to evaluate and correlate different geological and technical aspects with predisposing and triggering factors that can control these phenomena. The study area is the upper part of the Susa Valley where 12 small basins were selected. For each of them, lithological, geomorphological, climatic and technical information were mapped and analysed. Debris-flow triggering conditions, flow and depositional processes were related to physical characteristics of the basin that can be easily measured and quantified. At least three different groups of basins were found: G1) basins with one event each 4–6 years, characterised by massive or blocky calcareous rocks, G2) basins with more than one event per year that show an abundance of layered or sheared fine-grained rocks and G3) basins with recurrence levels exceeding 10 years, activated only by heavy and prolonged rainfalls, marked by massive or blocky coarse-grained igneous rocks. Furthermore, important morphometric differences were found. These considerations are useful in terms of hazard zonation and risk mitigation.     

Elastic properties and stability of coexisting 3T and 2M 1 phengite polytypes

The elastic properties of coexisting natural 3T and 2M ₁ phengite samples (Cima Pal, Sesia Zone; Val Savenca; Western Alps, Italy) with similar chemical compositions have been studied by room temperature–high pressure powder diffraction, using synchrotron radiation on the ID9A beam-line at ESRF (Grenoble, France). The P–V curves have been modelled by the Birch–Murnaghan model; a third-order expansion fitted to the experimental data yields for 3T and 2M ₁ K ₀=60.4(±0.7) GPa, K′=5.79(±0.11) at V ₀=703.8851 ų, and K ₀=57.3(±1.0) GPa, K′=6.97(±0.24) at V ₀=938.8815 ų, respectively. The relative stability of 3T vs. 2M ₁ has been explored as a function of pressure and temperature in terms of configuration and deformation contributions to the Gibbs energy, using the elastic properties determined here and other thermodynamic parameters from earlier investigations. The results presented agree with the hypothesis of stability of the 3T polytype in the high pressure regime.     

A semi-empirical approach to model pressure dependence of elastic moduli in granular media accounting for variations of coordination-number and Poisson-ratio

The effective medium theory based on the Hertz–Mindlin contact law is the most popular theory to relate dynamic elastic moduli (or elastic velocities) and confining pressure in dry granular media. However, many experimental results proved that the effective medium theory predicts pressure trends lower than experimental ones and over-predicts the shear modulus. To mitigate these mispredictions, several evolutions of the effective medium theory have been presented in the literature. Among these, the model named modified grain contact theory is an empirical approach in which three parametric curves are included in the effective medium theory model. Fitting the parameters of these curves permits to adjust the pressure trends of the Poisson ratio and the bulk modulus. In this paper, we present two variations of the modified grain contact theory model. First, we propose a minor modification in the fitting function for the porosity dependence of the calibration parameters that accounts for non-linearity in the vicinity of the critical porosity. Second, we propose a major modification that reduces the three-step modified grain contact theory model to a two-step model, by skipping the calibration parameter–porosity fit in the model and directly modelling the calibration parameter–pressure relation. In addition to an increased simplicity (the fitting parameters are reduced from 10 to 6), avoiding the porosity fit permits us to apply the model to laboratory data that are not provided with accurate porosity measurements. For this second model, we also estimate the uncertainty of the fitting parameters and the elastic velocities. We tested this model on dry core measurements from literature and we verified that it returns elastic velocity trends as good as the original modified grain contact theory model with a reduced number of fitting parameters. Possible developments of the new model to add predictive power are also discussed.     

Jitter radiation images, spectra and light curves from a relativistic spherical blastwave

We consider radiation emitted by the jitter mechanism in a Blandford–McKee self-similar blastwave. We assume the magnetic field configuration throughout the whole blastwave meets the condition for the emission of jitter radiation and we compute the ensuing images, light curves and spectra. The calculations are performed for both a uniform and a wind environment. We compare our jitter results to synchrotron results. We show that jitter radiation produces slightly different spectra than synchrotron, in particular between the self-absorption and the peak frequency, where the jitter spectrum is flat, while the synchrotron spectrum grows as  ν^1/3  . The spectral difference is reflected in the early decay slope of the light curves. We conclude that jitter and synchrotron afterglows can be distinguished from each other with good quality observations. However, it is unlikely that the difference can explain the peculiar behaviour of several recent observations, such as flat X-ray slopes and uncorrelated optical and X-ray behaviour.     

Predictive Power Evaluation of a Physically Based Model for Shallow Landslides in the Area of Oltrepò Pavese,Northern Italy

The use of real-time landslide early warning systems is attracting the attention of the scientific community, since it allows to assess “where” and “when” a shallow rainfall-induced landslide might occur by coupling rainfall amounts, hydrological models and slope-stability analysis. The paper deals with the main results of a back analysis, which refers to the application of a physically based stability model [Shallow Landslides Instability Prediction (SLIP)] on regional scale. The analysis concerns the occurrence of some recent rainfall-induced shallow landslides in the municipal territory of Broni, in the area of Oltrepò Pavese (Northern Italy). The study area is a hilly region 2.4 km² wide, where more than 40 % of the territory has slopes steeper than 15° and altitudes are between 90 and 250 m a.s.l. As regards the geologic setting, clayey-silty shallow colluvial deposits, with a maximum thickness of about 3 m, overlap a bedrock made of clayey shales, calcareous flysch and marls. The SLIP model is based on the limit equilibrium method applied to an infinite slope and on the Mohr–Coulomb strength criterion for the soil. By assuming that the main hydro-geotechnical process that leads to failure is the saturation of parts of the soil, the model allows to take into account the condition of partial saturation of the soil. The safety factor (F _S ) of a slope is also function of previous rainfalls. After the implementation of the model at territory scale, the input data have been introduced through a geographic information systems platform. In the current paper we mainly intend to evaluate the performance of SLIP at catchment scale, by comparison to (1) observed landslide events and (2) another well-established physically based model (TRIGRS). Further, we want to assess the suitability of the model as early warning tool. The results produced by the model are analyzed both in terms of safety factor maps, corresponding to some particular rainfall events, and in terms of the time-varying percentage of unstable areas over a 2-year span period. The paper shows the comparison between observed landslide localizations and model predictions. A quantitative comparison between the SLIP model and TRIGRS is presented, only for the most important event that occurred during the analyzed period. Overall, the results of the stability analyses based on observed rainfalls show the capability of the SLIP model to predict, in real-time and on a wide area, the occurrence of the analyzed phenomena.     

Acoustic and electromagnetic properties of soils saturated with salt water and NAPL

Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.     

Estimation of gas-hydrate concentration and free-gas saturation at the Norwegian-Svalbard continental margin

We estimate the concentration of gas hydrate and free gas at an area located to the north of the Knipovich Ridge (western Svalbard margin). The method is based on P-wave velocities computed by reflection tomography applied to multicomponent ocean-bottom seismometer data. The tomographic velocity field is fitted to theoretical velocities obtained from a poro-elastic model based on a Biot-type approach (the interaction between the rock frame, gas hydrate and fluid is modelled from first physical principles). We obtain average hydrate concentrations of 7% and maximum free-gas saturations of 0.4% and 9%, depending on the saturation model.     

Turbulent Intensities and Velocity Spectra for Bare and Forested Gentle Hills: Flume Experiments

To investigate how velocity variances and spectra are modified by the simultaneous action of topography and canopy, two flume experiments were carried out on a train of gentle cosine hills differing in surface cover. The first experiment was conducted above a bare surface while the second experiment was conducted within and above a densely arrayed rod canopy. The velocity variances and spectra from these two experiments were compared in the middle, inner, and near-surface layers. In the middle layer, and for the canopy surface, longitudinal and vertical velocity variances () were in phase with the hill-induced spatial mean velocity perturbation (Δu) around the so-called background state (taken here as the longitudinal mean at a given height) as predicted by rapid distortion theory (RDT). However, for the bare surface case, and remained out of phase with Δu by about L/2, where L is the hill half-length. In the canopy layer, wake production was a significant source of turbulent energy for , and its action was to re-align velocity variances with Δu in those layers, a mechanism completely absent for the bare surface case. Such a lower ‘boundary condition’ resulted in longitudinal variations of to be nearly in phase with Δu above the canopy surface. In the inner and middle layers, the spectral distortions by the hill remained significant for the background state of the bare surface case but not for the canopy surface case. In particular, in the inner and middle layers of the bare surface case, the effective exponents derived from the locally measured power spectra diverged from their expected  − 5/3 value for inertial subrange scales. These departures spatially correlated with the hill surface. However, for the canopy surface case, the spectral exponents were near  − 5/3 above the canopy though the minor differences from  − 5/3 were also correlated with the hill surface. Inside the canopy, wake production and energy short-circuiting resulted in significant departures from  − 5/3. These departures from  − 5/3 also appeared correlated with the hill surface through the wake production contribution and its alignment with Δu. Moreover, scales commensurate with Von Karman street vorticies well described wake production scales inside the canopy, confirming the important role of the mean flow in producing wakes. The spectra inside the canopy on the lee side of the hill, where a negative mean flow delineated a recirculation zone, suggested that the wake production scales there were ‘broader’ when compared to their counterpart outside the recirculation zone. Inside the recirculation zone, there was significantly more energy at higher frequencies when compared to regions outside the recirculation zone.