To investigate the millimetre-scale surface processes caused by natural
rainfall, an undisturbed sample of badlands soil (1 m long, 0.5 m wide and
0.15 m thick) was carefully extracted. The sample is composed of black marl
soil from a badlands area of the Draix Observatory (SE France). After
extraction, the undisturbed sample was placed at the same slope angle
(45
Small black marl watersheds have strong responses to climate forcing (Malet et al., 2007). This study aims to better understand the erosional behaviour and the micro-morphological evolution of these materials caused by micro-scale mass movements in response to precipitation. The goal of this experiment was to extract an undisturbed sample of black marl and expose this material to a natural rainfall event in order to monitor the sample's surface evolution using a 3-D laser scanner. This monitoring has permitted us to study the micro-topographic surface deformation and erosion processes that may have an impact on infiltration rates during rainfall events (Mitchell and van Genuchten, 1993; Römkens and Prasad, 2006). In addition, these surface changes may also play a role in the triggering of landslides (Galeandro et al., 2014). Contrary to the original intention, the duration of the experiment and the rainfall intensity did not permit investigation of splash erosion, which can be important when rainfall is intense (Selby, 1993), because its effect was below the resolution of our acquisition and because the intensity of the rainfall was rather low.
In the past, several studies of artificial rain simulations have been performed at the Draix Observatory (ORE Draix), which is dedicated to the research of mountain hydrology and erosional processes. Previous studies focused on the measurement of sediment transport (Oostwoud Wijdenes and Ergenzinger, 1998) and runoff as a function of precipitation (Mathys et al., 2005). The badlands ground surface evolution has already been monitored during rainfall simulations by using a pin-type micro-relief meter and photography (Torri et al., 1999; Mathys et al., 2005), but these tools were not able to reach the same precision as a laser scanner. Recent studies have shown the potential for observing and characterising erosional processes at the level of micro-topography using a laser scanner (Schmid et al., 2004; Barneveld et al., 2013).
The soil sample was extracted from the Draix experimental site (ORE Draix, IRSTEA), a badlands area near the city of Digne-les-Bains in the southern French Alps. The experiment was performed in Lausanne (Switzerland), where the precipitation is similar to the Alpine region. For the monitoring of the sample surface, a terrestrial laser scanner (TLS) was used to (1) identify the millimetre-scale processes that control erosion and infiltration by (2) quantifying the swelling of the material in the marls during rainfall events and (3) mapping all possible modifications of the terrain surface.
The badlands near the village of Draix are composed of weathered black marls
of Middle Jurassic age (Callovian–Oxfordian). The black marl formation is more
than 2000 m thick in certain places (Antoine et al., 1995). The study site has
no vegetal cover, and the regolith is usually approximately 40 cm to 1 m
thick (Maquaire et al., 2002; Antoine et al., 1995). The upper approximately
10 cm of the regolith is loose detrital material composed of local clasts
and platelets (Maquaire et al., 2002). This regolith corresponds to a sandy
loam when it has been exposed for long enough to be disaggregated by
weathering (Antoine et al., 1995). When the upper regolith layer is fresh,
it can be considered loamy sand or sand. Below the regolith is a layer of
plate-like unstructured rock. Finally, compacted regolith 10–20 cm thick
lies in contact with the bedrock (Maquaire et al., 2002). The clay size
fraction of the black marls located in this region has been found to be
35
The sample of soil used for this experiment was extracted from a marl
outcrop with a 45
On 31 May 2011, the soil sample was exposed to natural rainfall from 11:01
to 17:47 LT (local time used throughout the text). During the experiment, the soil sample was kept in its metal
extraction casing and was tilted at 45
The soil sample was scanned every 30 min using a ground-based TLS Leica
ScanStation II, which produces point clouds (
Summary of the TLS scans in the TLS campaign experiment on 31 June 2010. This table compiles the information concerning the scan time, the number of points before and after cleaning, the mean spacing between the points and the density of points.
The sample was dry at the beginning of the experiment. The total precipitation during the experiment was 5.5 mm, which corresponds to approximately 1.3 L entering the metal case. Most of the rainfall was absorbed by the sediment, limiting the transport of sediment by runoff. Consequently, the sample did not reach full saturation. We suspect that a small quantity of evaporation may have occurred during the experiment; however, we did not weigh the box before and after the experiment, and this minor influence has been neglected in our study.
The surface evolution was monitored for 6 h and 46 min by 12 successive laser scan acquisitions. The first scan and the last scan (acquired at 11:01 and 17:47, respectively) had a point spacing of approximately 0.001 m with a duration time of 8 min. All the other scans had a 0.002 m point spacing for a duration of 2 min. From 0 to 50 m, 50 % of the laser beam is at a diameter of 4 mm at full width half height (FWHH) (Leica Geosystems AG, 2007). The scans were acquired vertically from left to right, which means that 1 cm in width is scanned in less than 2.4 s for the 2 min scans and less than 9.6 s for the 8 min scans. The scan distance was 2 m. The instrument was not moved throughout the experiment, meaning that the position and the orientation of the scans are identical for all acquisitions. Therefore, no scan alignment was necessary to compare the obtained point clouds.
The point clouds were “manually” cleaned of the points that were not imaging the surface of the sample, i.e. the metal casing sides, the background of the scenery and several artefacts, including points not located on the terrain's surface, such as rain drops. Only the surface of the sample was kept. Every cleaned scan has a very high point density: more than 400 000 points for the first and last acquisitions and a minimum of 110 000 points for the other scans (Table 1).
Each TLS point cloud was first rotated by 45
A 10 cm long downhill section using the points from a 1 cm wide swath of the original TLS data point cloud compared with the 1 mm DEM profile (17:47 scan data). The profile has been processed with an inverse distance interpolation with a power of 1.
A comparison of TLS scans. The colours represent topographic changes between 11:01 and 17:47. The figure also includes a detailed view of the surface. The locations of Figs. 4, 6 and 7 are indicated.
The precipitation lasted from 13:30 until the end of the experiment, with a
30 min hiatus starting at approximately 16:00. The maximum rainfall
intensity (5.2 mm h
The micro-topographic changes that occurred along transect 1–2 are shown in Figs. 3 and 4. Figure 4c shows a depression (in red) formed at the top of a small ridge, and the surface elevation of a depression below this ridge increased (in blue). This change can be interpreted as a downward mass movement at the millimetre scale. The small particles moved downward and filled a desiccation crack. This phenomenon occurred 30 min after the highest intensity rain event, between 17:00 and 17:30 (Fig. 5).
When the rain intensity reached 1 mm h
Profile changes showing mass movement occurring between 17:00
and 17:30.
Another observed process is linked to changes in the soil surface by lateral expansion (Figs. 6, 7). This process occurred continuously through the closing of desiccation cracks, which were present on the surface of the sample.
The final observed process was the stripping of soil particles by the kinetic energy of the raindrops. Although these changes were observed by the authors during the experiment, their magnitude was unfortunately too small to be significantly monitored by the TLS used in this study.
Observed phenomenon relative to the precipitation data. The total rain amount in millimetres, the cumulative precipitation and the rain intensity are presented. The rise is shown in millimetres, and the lateral expansion intensity and mass movement processes are presented in relative scales.
The topographic changes observed along transect 1–2 (Figs. 3 and 4) occurred almost instantaneously, which excludes the rain splashing process. These changes are interpreted as a micro-landslide. Based on the duration of the scans and the size of the transported mass, we assume that it is highly unlikely that a scan had just crossed the region of interest when the movements occurred within it (0.05 m scanned in approximately 50 s). In addition, the two profiles before and after the occurrence of the changes mimic the profile changes usually observed in real landslides. The observed process is likely related to the initiation of miniature debris flows (MDFs) observed by Oostwoud Wijdenes and Ergenzinger (1998); however, MDFs require heavy rainfall.
3-D views of the evolution of lateral expansion, which closes a
crack between 11:01
The observed rise in the surface, which occurred following a 30 min delay after the initiation of rainfall, is certainly linked to the swelling of the material. The soil sample experienced swelling after the first rainfall intensity peak, contraction 30 min after the cessation of rainfall and renewed rising once the rain started again. This measured cyclic behaviour demonstrates that the swelling and contraction of the soil surface is a reversible process. We can assume that this process is linked to moderate rainfall intensities, which allow water to infiltrate the fine-grained material, causing swelling. This process is not necessarily caused by clay minerals because they are present only in small quantities (primarily illite) in the study area (Antoine et al., 1995). The swelling dissipates rapidly for moderate rainfall events because of the diffusion of water when the rainfall stops, which leads to a decrease in the effect of the water.
The crack-closing lateral expansion of the surface is also certainly linked to swelling but is not reversed when the rain stops. Figures 6 and 7 clearly show that the material has expanded and that the material has not been transported because no deposition was observed at the bottom. We have not found any definitive explanation that accounts for the difference between the rising process and the lateral expansion process. However, the difference must be related to gravity, which increases the effect of swelling downward. When the two sides of the crack make contact, the moistened zone doubles its thickness, decreasing the water diffusion in the material. In addition, these processes must be components of creeping (Selby, 1993), likely because the retreat by drying is less effective in the downslope portion, leading to slow progressive downward movements.
The above interpretations are also important for the understanding of the infiltration process. Cracks clearly play an important role in the infiltration rates (Mitchell and van Genuchten, 1993; Römkens and Prasad, 2006) and consequently in the destabilisation of slopes (Stumpf et al., 2013; Galeandro et al., 2014). The processes analysed here play a role in the closure of cracks, as shown in Fig. 7. In the present case study, we demonstrated how micro-scale infiltration can influence the degradation of soil surface by inducing downward mass movements that are not reversible.
We have also shown here the great potential of high-resolution three-dimensional TLS or photogrammetry point clouds for the analysis of the processes that lead to erosion through surface mass movements at the millimetre scale. Investigations of erosional processes using point clouds are increasing in number. These studies use laser scanners for either micro-scale surface imaging (Schmid et al., 2004; Barneveld et al., 2013) or measuring crack apertures (Sanchez et al., 2013). In addition, photogrammetry and structure from motion (SfM) methods are now being developed to analyse soil surfaces (Snapir et al., 2014).
This paper shows that monitoring the changes at the millimetre scale to
examine soil surface changes and erosion is now possible. This development
will aid in designing future experiments to analyse certain processes, such
as swelling, crack closure, micro-landslides, and initiation of MDFs. With
heavier rainfall, sediments will be mobilised and transported across longer
distances, enabling the study of MDFs and the formation of rills. This study
also demonstrates that material and rain intensity must be suitable to
permit the efficient detection of rain splash processes and associated
erosion; specifically, a rainfall intensity of greater than 20 mm h
The authors would like to thank the Observatoire de Recherche en Environnement (ORE) of Draix for letting us sample soil in a protected area. We would also like to thank the UNIBAT service of the University of Lausanne (UNIL) for providing meteorological data from their weather station. We also greatly appreciate the comments and suggestions of our colleagues M.-H. Derron, Benjamin Rudaz and Antonio Abellán. We thank American Journal Experts for the improvements of the English language of this paper. We dedicate this paper to the first author Jacques Bechet, who died in a snow avalanche 3 days after the acceptance of this paper. The content of this paper is an expression of its great ingenuity, curiosity and research spirit he shared with his co-worker Julien Duc. We will ever remember his enthusiasm. Edited by: M. Mikos