A critical zone observatory dedicated to suspended sediment transport: The meso-scale Galabre catchment (southern French Alps)

The 20 km² Galabre catchment belongs to the French network of critical zone observatories (OZCAR; Gaillardet et al., Vadose Zone Journal, 2018, 17(1), 1–24). It is representative of the sedimentary lithology and meteorological forcing found in Mediterranean and mountainous areas. Due to the presence of highly erodible and sloping badlands on various lithologies, the site was instrumented in 2007 to understand the dynamics of suspended sediments (SS) in such areas. Two meteorological stations including measurements of air temperature, wind speed and direction, air moisture, rainfall intensity, raindrop size and velocity distribution were installed both in the upper and lower part of the catchment. At the catchment outlet, a gauging station records the water level, temperature and turbidity (10 min time-step). Stream water samples are collected automatically to estimate SS concentration-turbidity relationships, allowing quantification of SS fluxes with known uncertainty. The sediment samples are further characterized by measuring their particle size distributions and by applying a low-cost sediment fingerprinting approach using spectrocolorimetric tracers. Thus, the contributions of badlands located on different lithologies to total SS flux are quantified at a high temporal resolution, providing the opportunity to better analyse the links between meteorological forcing variability and watershed hydrosedimentary response. The set of measurements was extended to the dissolved phase in 2017. Both stream water electrical conductivity and major ion concentrations are measured each week and every 3 h during storm events. This extension of measurements to the dissolved phase will allow progress in understanding both the origin of the water during the events and the partitioning between particulate and dissolved fluxes of solutes in the critical zone. All data sets are available at https://doi.osug.fr/public/DRAIXBLEONE_GAL/index.html .     

Assessment of intermediate fine sediment storage in a braided river reach (southern French Prealps)

This paper presents a field investigation on river channel storage of fine sediments in an unglaciated braided river, the Bès River, located in a mountainous region in the southern French Prealps. Braided rivers transport a very large quantity of bedload and suspended sediment load because they are generally located in the vicinity of highly erosive hillslopes. Consequently, these rivers play an important role because they supply and control the sediment load of the entire downstream fluvial network. Field measurements and aerial photograph analyses were considered together to evaluate the variability of fine sediment quantity stored in a 2·5‐km‐long river reach. This study found very large quantities of fine sediment stored in this reach: 1100 t per unit depth (1 dm). Given that this reach accounts for 17% of the braided channel surface area of the river basin, the quantities of fine sediment stored in the river network were found to be approximately 80% of the mean annual suspended sediment yields (SSYs) (66 200 t year^?1), comparable to the SSYs at the flood event scale: from 1000 t to 12 000 t depending on the flood event magnitude. These results could explain the clockwise hysteretic relationships between suspended sediment concentrations and discharges for 80% of floods. This pattern is associated with the rapid availability of the fine sediments stored in the river channel. This study shows the need to focus on not only the mechanisms of fine sediment production from hillslope erosion but also the spatiotemporal dynamics of fine sediment transfer in braided rivers. Copyright © 2010 John Wiley & Sons, Ltd.     

Magnitude–frequency relationships of debris flows and debris avalanches in relation to slope relief

Probability of occurrence, hazard intensity and encounter probability are key parameters in the quantitative risk analysis (QRA) of landslides. All are strongly dependent on magnitude of the landslides. As a result, magnitude–frequency analysis should be a part of QRA. Deriving representative magnitude–frequency relationships for debris avalanches and debris flows, however, is difficult. One key problem is illustrated with the example of a unique database from the coastal region of British Columbia, Canada, which was compiled entirely from detailed ground investigations. The magnitude of debris avalanches and debris flows is not an independent statistical quantity, but a function of the scale of a given slope, as characterized by the slope length. Thus, attempting to derive probability and magnitude for a given location or sub-region from a regionally-derived magnitude–frequency curve may lead to incorrect predictions. The same problem is pertinent to the application of the same approach to any type of landslide in which the largest combined dimension of the source volume (including entrainment) is of the same order as the length of the slope. It is recommended that greater emphasis be placed on site-specific geological observations, at the expense of generalized statistics.     

Test of three methods to detect the overbank flow from water level time‐series analysis

In this study, we propose to identify morphological and hydraulic characteristics related to overbank flows in the water level time‐series available at many gauged stations. The results obtained at 13 river stations (the catchment sizes vary from 10 to 1700 km² with contrasted geology, land use and rainfalls regime) show that overflow mechanisms at the river‐reach scale can be systematically identified in the water level frequency distribution estimated with the peaks‐over‐threshold (POT) method. A first level (Lts1) was in the range of the incipient flooding onto the floodplain. Even if the definition of this level is variable in terms of flooded area at the reach scale, this method can be useful in providing a first estimation of the bankfull level for many gauged stations, without complex and costly field surveys. A second level (Lts2) was systematically detected on average 38 cm above the topographic flat floodplain elevation. The Lts2 inflection in the water‐level frequency distribution is assumed to reflect a composite effect of catchment hydrology and local hydraulics and channel geometry, without possibility to make a clear distinction between both processes at the moment. The local or reach scale effect would probably play an essential role in the frequency distribution as flood attenuation at catchment scale may explain the inflection Lts2 at only three sites. In light of the knowledge acquired in laboratory flumes with compound channels, most of the time Lts2 level would correspond to the level of resumption of flow both in the main channel and the floodplain. Once this method is validated in various physiographic contexts, it should apply to many hydrometric stations for both synchronic (e.g. regional analysis) and diachronic analysis (e.g. evolution over time of the bankfull discharge) to evaluate anthropic impacts on river morphology. Copyright © 2010 John Wiley & Sons, Ltd.     

Morphological methods and dynamic modelling in landslide hazard assessment of the Campania Apennine carbonate slope

Landslide risk of the Campanian carbonate slopes covered by pyroclastic deposits is mainly connected with the occurrence of high-velocity debris avalanches and debris flows. Analyses show that flows initiate as small translational slides in the pyroclastics. The failure process is controlled by the interaction of both natural and human-induced factors. Geomorphological settings play a decisive role in locating the source failures. Therefore, the crucial aspects in landslide hazard and risk assessment are: (a) recognise the geomorphological control factors, (b) determine parameters defining landslide intensity (velocity, volume, depth of deposit) and (c) predict landslide runout distance. An approach combining geomorphology and numerical analysis has been adopted in the work reported here. Potential future landslide intensity scenarios are simulated predicting the runout behaviour of potential instabilities by using a dynamic model previously calibrated by back-analysing observed events of similar scale and type. The selected area is a sector of the Avella Mountains having the same geomorphological environment as the 1998 Sarno landslides (Campania, Southern Italy).     

Velocity and runout simulation of destructive debris flows and debris avalanches in pyroclastic deposits,Campania region,Italy

The Campanian Apennines are characterized by the presence of monocline ridges, mainly formed by limestone. During the periods of volcanic activity of the Somma-Vesuvius and Phlegrean Fields, the ridges were mantled with pyroclastic materials in varying thickness. The pyroclastics have been involved in destructive landslides both in historical time and in the recent past (1997, 1998, 1999). The landslides occur following intense and prolonged rainfalls. In some cases, landslides extended up to 4 km into the surrounding lowlands and reached towns, causing severe destruction and over 200 deaths. Generally, the landslides begin as small debris slides that develop into large, shallow debris avalanches or debris flows involving pyroclastic horizons and colluvial soils (0.5–2 m thick) on steep and vegetated slopes, often at the heads of gullies. During motion, the landslide materials eroded vegetation and soils from the slope, so that the moving material volume tended to increase. Then, proceeding towards and beyond the base of the slopes, the phenomena evolved into hyperconcentrated streamflow due to dilution by incorporating water. The results of motion analyses are described. An empirical rheological relationship was used including two principal terms that depend on the total normal stress and on the flow velocity. On this basis, the model has simulated the velocity and duration of debris avalanches and the distribution of the deposits. The selected areas were those of Sarno/Quindici and Cervinara, where a large amount of data is available both on the material properties and geomorphological setting. It was found that the majority of the cases at the two sites can be simulated successfully with only one specific pair of rheologic parameters. This provides the possibility for first-order predictions to be made of the motion of future landslides. Such predictions will be a valuable tool for outlining potential hazard areas and designing remedial measures.     

Rockfall rebound: comparison of detailed field experiments and alternative modelling approaches

The accuracy of rockfall trajectory simulations mainly rests on the calculation of the rebound of fragments following their impact on the slope. This paper is dedicated to the comparative analysis of two rebound modelling approaches currently used in rockfall simulation using field experiments of single rebounds. The two approaches consist in either modelling the rock as a single material point (lumped mass approach) or in explicitly accounting for the fragment shape (rigid body approach). A lumped mass model accounting for the coupling between translational and rotational velocities and introducing a slope perturbation angle was used. A rigid body approach modelling the rocks as rigid locally deformable (in the vicinity of the contact surface) assemblies of spheres was chosen. The comparative analysis of the rebound models shows that both of them are efficient with only a few parameters. The main limitation of each approach are the calibration of the value of the slope perturbation (‘roughness’) angle, for the lumped mass approach, and the estimation of the rock length and height from field geological and historical analyses, for the rigid body approach. Finally, both rebound models require being improved in a pragmatic manner to better predict the rotational velocities distribution. Copyright © 2012 John Wiley & Sons, Ltd.