Short communication: Runout of rock avalanches limited by basal friction but controlled by fragmentation
- 1GFZ German Research Centre for Geosciences, Helmholtz Centre Potsdam, Telegrafenberg, 14473 Potsdam, Germany
- 2Njord centre, Department of Geosciences, University of Oslo, P.O. Box 1048, 0316 Oslo, Norway
- 3Van Hall Larenstein University of Applied Sciences, Larensteinselaan 26a, 6882 Velp, the Netherlands
Correspondence: Matthias Rosenau (firstname.lastname@example.org)
Rock avalanches produce exceptionally long run-outs that correlate with their rock volume. This relationship has been attributed to the size-dependent dynamic lowering of the effective basal friction. However, it has also been observed that run-outs of rock avalanches with similar volumes can span several orders of magnitude, suggesting additional controlling factors. Here, we analyse analogue models of rock avalanches, with the experiments designed to test the role of dynamic fragmentation. We show that for a fixed low basal friction, the run-out of experimental rock avalanches varies over 2 orders of magnitude and is determined by their degree of fragmentation, while the basal friction acts only as an upper limit on run-out. We interpret the run-out's dependence on fragmentation as being controlled by the competition between mobility enhancing spreading and energy-consuming fragmentation limited by basal friction. We formalize this competition into a scaling law based on energy conservation, which shows that the variation in the degree of fragmentation can contribute to the large variation in run-out of rock avalanches seen in nature.
With volumes larger than 109 m3 and speeds reported at over 150 km/h (Campbell, 1989) and possibly up to 100 m/s (Legros, 2002), the destructive power of rock avalanches is unprecedented. They are exceptional hazards produced when very large rockslides disintegrate during transport (Hungr et al., 2013). The travel distance of the deposit front, or run-out, is an important measure for hazard assessment (Vaunat and Leroueil, 2002) and is generally found to be more than 10 times longer than the fall height (Hsü, 1975). This suggests low effective basal friction μeff, which is usually attributed to various dynamic weakening processes (e.g. Kent, 1966; Shreve, 1968; Hsü, 1975; Melosh, 1979; Campbell, 1989; Pudasaini and Miller, 2013; Legros, 2002; Lucas et al., 2014; Wang et al., 2017) or additional basal erosion processes (e.g. Hungr and Evans, 2004; Pudasaini and Fischer, 2020).
Field observations of the displacement of rock avalanches are typically given by the ratio of vertical (H) and horizontal (L) distance from the deposit's front to the top of the main scarp. The resulting ratio
is known as the Heim's ratio (Heim, 1882, as cited in Hsü, 1975) and serves as a proxy for μeff, where it is called the “apparent” coefficient of friction (Manzella and Labiouse, 2012). One of the best established but perhaps least understood observations of rock avalanches is the dependence of the Heim's ratio on volume: rockslides below a size of approximately 106 m3 all have a relatively constant Heim's ratio of ∼0.4–0.7, but for larger rockslides Heim's ratio decreases with volume, reaching values <0.1 for volumes larger than 109 m3 (Pudasaini and Miller, 2013; Lucas et al., 2014). This suggests a scale-dependent mechanism of decreasing apparent friction with volume that becomes dominant at large volumes (Davies and McSaveney, 1999). Analytical modelling and numerical simulation involving lubrication mechanisms by Pudasaini and Miller (2013) and Lucas et al. (2014) provided mechanical explanations for this observation. Importantly, however, even within a narrow range of volumes, run-outs are seen to span orders of magnitude, suggesting additional controlling factors on run-out that are insensitive to size. Runout variability also lies in part in the fact that the run-out is defined by the front of the deposits, and it therefore contains the combined effect of both translation and spreading of the rock mass. The additional travel distance caused by spreading can have a profound effect on the run-out (Staron and Lajeunesse, 2009), especially if the effective basal friction is low.
Recently, the process of dynamic fragmentation has received increased attention from the research community, and much progress has been made in our understanding of its role in the dynamics of rock avalanches (Locat et al., 2006; Imre et al., 2010; Bowman et al., 2012; Pudasaini and Miller, 2013; De Blasio and Crosta, 2015; Haug et al., 2016; Zhao et al., 2017, 2018; Lin et al., 2020; Gao et al., 2020; Knapp and Krautblatter, 2020). Firstly, one may expect that the finer the material is, the more flow-like the behaviour will be, increasing its mobility and allowing the rock mass to spread more easily (Locat et al., 2006; Wang et al., 2017; Zhao et al., 2018). Secondly, models of fragmenting rockslides suggest that dynamic fragmentation actively increases the spreading (Bowman et al., 2012; De Blasio and Crosta, 2015; Lin et al., 2020). However, fragmentation has also been shown to consume energy (Haug et al., 2016; Zhao et al., 2017; Lin et al., 2020), potentially at a cost to the run-out length. Clearly, understanding the integrated effect of fragmentation on the run-out dynamics of rock avalanches requires more analysis.
To study the effects of friction and fragmentation on rock avalanche dynamics, here we analyse analogue models of dynamically fragmenting rockslides. We assume that there exists some mechanism (or a set of mechanisms) that causes a low effective coefficient of basal friction, which we set to 0.15–0.20 in our models. To isolate the scale-independent effect of fragmentation we keep both the volume and friction within a narrow range in our models compared to nature. Note that this approach explicitly excludes dynamic weakening mechanisms that are suspected in natural prototypes. Specifically, our models do not include fluids, and frictional heating is insignificant such that lubrication mechanisms (e.g. Pudasaini and Miller, 2013; Lucas et al., 2014) do not play a role. Granular pressurization (e.g. Imre et al., 2010) is also not considered significant in our experiments because of the low energy involved. Other potentially important mechanisms like bedrock erosion (e.g. Hungr and Evans, 2004; Pudasaini and Fischer, 2020) are excluded here for simplicity. The experimental design therefore means that the observed variation in Heim's ratio is due to fragmentation and dry friction. We describe the dependence observed between the run-out and the degree of fragmentation in the form of a scaling law. Finally, we compare our experimental results to a set of natural data and discuss their relevance to natural systems. All data underlying this study and additional relevant data are published as open access in Haug et al. (2021).
In the experiments originally documented in Haug et al. (2016), a block of height h and length l0 (width =l0) of rock analogue material is gravitationally accelerated down a plate held at an angle of 45∘ to the horizontal (Fig. 1). After 1 m of travel, the sample impacts a horizontal plate causing it to fragment. Once the sample fragments have slid onto the horizontal plane, they spread and decelerate due to the internal and basal frictional interaction before finally coming to rest. We use silicate glass as our substrate, upon which the basal friction coefficient is ca. 0.15–0.20 (Haug et al., 2016) – similar to the lowermost values found in natural prototypes (Pudasaini and Miller, 2013; Lucas et al., 2014). The analogue rock material is a cemented fluvial quartz sand. The sand is cemented by mixing it with water and gypsum or potato starch and is left to set for 2 d (for gypsum cement) or heated for 15 min in a 900 W microwave (for potato starch cement). The cohesion of the material can be controlled by the type and amount of cement added to the mixture, allowing control of the strength of the material over several orders of magnitude. The internal friction coefficient relevant for fragmenting intact material is 0.7 and reduces to 0.6 when fragments interact (see Haug et al., 2014, 2016, for details on the experimental setup).
The three main observables from the experiments are (i) the degree of fragmentation (mc), (ii) the Heim's ratio (), and (iii) the normalized deposit length (). We characterize the degree of fragmentation through the total mass of the sample divided by the mass of the largest fragment (). We choose this rather simple parameter, which has been validated and benchmarked against breakage parameters used by previous studies in Haug et al. (2016), as a tradeoff between capturing the process accurately in models and the accessibility of the equivalent information in nature. To define Lspread, we consider the mass-weighted average position of the most proximal and distal 5 % of the total mass. This defines a rim that is a more robust run-out estimate than using single fragment positions, as done by Haug et al. (2016), and is at the same time accessible both experimentally and empirically. We normalize Lspread by fall height H in order to have a parameter describing the conversion of potential energy into spreading that is equivalent to Heim's ratio.
The experimental data analysed here come from two series of experiments with varying degrees of fragmentation: (i) one series of experiments where the thickness-to-length ratio () of the samples has been varied between 0.033 and 0.49 (corresponding to a 1 order of magnitude range in volume), while keeping the cohesion (C) constant at 14 kPa, and (ii) one series of experiments where the cohesion of the material is varied between 4 and 350 kPa, while keeping the thickness to length ratio constant at 0.13. In both series of experiments, the fall height (H) is kept constant at 0.71 m. Interested readers are referred to Haug et al. (2016) for details on the effect of cohesion and geometry on the degree of fragmentation. Additionally, two new experiments were performed to study the moment of fragmentation at high temporal resolution. For these experiments, the fragmentation of two samples with different cohesion but equal geometry (C=4 and 40 kPa, ) is considered. These two experiments have a fall height of 0.35 m, and data are captured by a camera with a frame rate of 500 Hz (see Haug et al., 2021, for movies of these experiments). Combining these sets of data from various experiments allows for covering a wide enough parameter space for the analysis in this study.
3.1 Experimental observations and interpretation
Figure 2 presents snapshots from two representative experiments, one with an intermediate-strong sample and one with a low-strength sample, illustrating the process of fragmentation. The stronger sample (Fig. 2a) is observed to fragment less than the weaker one (Fig. 2b). Thereafter, fragments of the stronger sample spread with limited interaction, while the fragments from the weaker sample collide and/or slide next to each other and deposition starts relatively early. We infer at the first order that while mobility generally increases with fragmentation, a higher amount of internal deformation is experienced, along with increased fragmentation and increased deposition.
To quantitatively analyse the experiments, we focus on the correlation between run-out and fragmentation and neglect all other parameters. This is justified by the collapse of experimental and natural data when plotting Heim's ratio against fragmentation in Fig. 3a. Qualitatively, Heim's ratio decreases rapidly for low to intermediate degrees of fragmentation, reaching a minimum at mc≈5 of about 0.2, and increases again slightly for higher degrees of fragmentation. A similar relation is observed between the length of the deposits (Fig. 3b), which increases with fragmentation until mc≈5 and slightly decreases beyond.
The rapidly decreasing Heim's ratio for mc<5 observed in Fig. 3a is likely linked to the increased spreading with fragmentation seen in Fig. 3b. A similar result was also obtained by previous analogue experiments (Bowman et al., 2012; Haug et al., 2016) and numerical models (De Blasio and Crosta, 2015; Zhao et al., 2017). However, here we show that the Heim's ratio is not simply decreasing with the degree of fragmentation but that it instead displays an optimum for mc≈5. Importantly, the lowest apparent basal friction, equivalent to the lowest Heim's ratio, is close to the implemented basal friction (i.e. friction coefficient of 0.15–0.20 between samples and glass). Therefore, all processes operating in our models (e.g. fragmentation, internal friction between fragments, deposition) tend to consume energy and thereby reduce run-out from its optimum (Haug et al., 2016). Considering the increased internal deformation observed with the degree of fragmentation (Fig. 2), the reduction of run-out and spreading for mc>5 appears to be the result of the increased energy dissipation through internal friction within the rock mass and an increase in basal friction as the sliding surface becomes rougher due to syn-sliding deposition (e.g. Pudasaini and Fischer, 2020). A loss of mass and therefore momentum due to deposition may additionally result in deceleration and reduced run-out as a function of mc (e.g. Pudasaini and Fischer, 2020). Consequently, the minimum of the Heim's ratio observed in Fig. 3a appears as the result of a competition between the spreading enhancing mobility and the energy-consuming fragmentation process.
3.2 A scaling law for run-out
The interplay between fragmentation and friction in a dry environment can be formalized into a scaling law by considering the conservation of energy. Generally, the conservation of energy of a sliding mass M requires that
where g is the gravitational acceleration, H is the vertical fall height, Lp is the entire travel path of the slide, and W is the sum of any other energy dissipating terms. Here, we have assumed a Coulomb friction coefficient μ at the base.
For the geometry of our experimental setup (see Fig. A1), as well as (roughly) for the set of selected rock avalanches, the Lp can be expressed in terms of the horizontal run-out L as
where Ls is the length and θ the angle of the slope, and l0 is the initial length of the slide. It is assumed that the additional travel length due to spreading is equal to half the deposit length (Lspread). Since l0 is expected to be very small compared to the other terms, it is neglected in further analysis. Inserting Eq. (3) into (2) and solving for L gives
where it is emphasized that both Lspread and W are expected to be functions of the basal friction, μ, internal friction, ϕ, and the degree of fragmentation, mc, as well as a possible non-linear dependence between Lspread and W. Rearranging Eq. (4) yields the Heim's ratio in the form of
A direct determination of the two last terms in Eq. (5) is difficult. However, based on the shape of the function of both the Heim's ratio and the Lspread plotted in Fig. 3, it appears that it can be reasonably described by an exponential function of mc:
Additionally, the experimental work by Haug et al. (2016) suggests that dissipative energy loss through fragmentation increases less for higher degrees of fragmentation and therefore can be described with a logarithmic function of mc:
Using these approximations, Heim's ratio can be expressed as
where α, β, and γ are constants to be empirically determined.
This equation describes the competition between spreading (proportional to ) and the increasing energy dissipation (proportional to log (mc)) with mc and its relation to friction. A best fit of this function to the natural and experimental data is presented in Fig. 3a (blue line), where , β=1.5, and γ=0.11. A fit constrained only the experimental data yields very similar results (, β=2.0, γ=0.11; see Fig. B1). This suggests spreading dominates run-out for low degrees of fragmentation (i.e. mc<5) but has little effect at high degrees of fragmentation as the exponential term approaches zero. At high degrees of fragmentation, the energy dissipation related to fragmentation therefore becomes increasingly relevant in controlling run-out. At mc≈5, i.e. when about 80 %–85 % of the volume is fragmented, a state of optimal mobility is reached with a Heim's ratio limited by the basal friction coefficient suggesting that energy is consumed mainly by basal friction, which then is the limiting factor for run-out.
3.3 Application to a natural data set
We compare our experimental results in Fig. 3 with data from nine rock avalanches reported by Locat et al. (2006) that show no clear volume dependence of run-out. This feature makes this data set ideal for testing whether a scale-independent process is operating besides dynamic basal weakening. However, not all the rock avalanches reported in Locat et al. (2006) are comparable to our experimental setup concerning material properties and geometries (Fig. 1). Based on slope geometry, the Queen Elizabeth slide is discarded because of its run-up on the opposite valley wall. Also discarded is the Charmonétier slide because of the sudden free-fall stage at the end of its transport. Additionally, the Arvel slide was observed to bulldoze soft material in front of it, and as such complexities are not considered in our models, this one is also neglected. Note that in all three discarded cases, the late-stage processes tend to increase the expected Heim's ratio, and they consistently plot above the trend of the other data in Fig. 3b.
Figure 3 displays remarkably similar trends between the experimental and the selected natural data that all follow the proposed scaling law. The data points from Jonas Creek (north) and Claps de Luc are observed to extend the trend from the experiments to higher Heim's ratios for low degrees of fragmentation, while the La Madelaine slide is observed to extend the trend of the experimental results of Heim's ratio to higher degrees of fragmentation (Fig. 3a). Its low spreading value (Fig. 3b) suggests that the reduction of spreading indicated by the experiments for mc>5 continues for even higher degrees of fragmentation. The agreement between these slide deposit lengths and the extrapolation of the experimental trend through Eq. (8) (Fig. 3b) supports the validity of our proposed scaling law. The Heim's ratios of the neglected slides are all, as expected, higher than the selected data set for their respective degrees of fragmentation, illustrating the importance of topography (e.g. opposite valley wall) and processes such as bulldozing.
The similarity seen between experimental and natural data suggests some universality concerning the empirical constants. Moreover, the similarity suggests that the rock avalanches considered here all have a close to constant effective basal friction of about 0.15–0.20. This implies that over the range of 2 orders of magnitude (from 2×106 to 90×106 m3) represented by this data set, the effective coefficient of friction of rock avalanches could be considered independent of volume. Consequently, our results suggest that the variation seen in Heim's ratio for these rock avalanches is not (only) caused by scale-dependent basal friction but instead by differing degrees of fragmentation. This shows that fragmentation plays a governing role in the run-out of rock avalanches and should be included in hazard assessments.
We studied the dynamics of fragmenting rock avalanches experimentally to unravel the control of basal friction versus fragmentation on run-out behaviour. We find that fragmentation causes both spreading and frictional interaction – competing processes that control the avalanche dynamics. Based on energy arguments we derive a scaling law with empirical constants that quantifies the relative importance of spreading and frictional interaction as a function of fragmentation. The scaling law approaches an extreme for which run-out is maximized and limited only by basal friction, which itself might be volume-dependent as suggested by earlier studies. The scaling law is validated against a natural data set verifying its applicability.
ØTH designed and run the experiments, derived the scaling law, and wrote the first draft of the manuscript. MRo and MRu assisted in the experiments. MRo, KL, and OO were involved in study design. All authors contributed to discussion and writing.
The authors declare that they have no conflict of interest.
The authors would like to thank to Frank Neumann and Thomas Ziegenhagen for construction and technical assistance. We thank Kirsten Elger and GFZ Data Services for publishing the data. We thank the editor Michael Krautblatter and the two anonymous reviewers, who provided very helpful comments that improved the manuscript. We thank Jon Bedford for providing a thorough language check of the manuscript.
The research has been funded by the Helmholtz Graduate Research School GEOSIM and
the German Ministry for Education and Research (BMBF, FKZ03G0809A) and the Deutsche
Forschungsgemeinschaft (DFG) through grant CRC 1114 “Scaling Cascades in Complex
Systems” (no. 235221301) project B01.
The article processing charges for this open-access publication were covered by the Helmholtz Centre Potsdam – GFZ German Research Centre for Geosciences.
This paper was edited by Michael Krautblatter and reviewed by two anonymous referees.
Bowman, E. T., Take, W. A., Rait, K. L., and Hann, C.: Physical models of rock avalanche spreading behaviour with dynamic fragmentation, Can. Geotech. J., 49, 460–476, https://doi.org/10.1139/t2012-007, 2012. a, b, c
Gao, G., Meguid, M. A., Chouinard, L. E., and Zhan, W.: Dynamic disintegration processes accompanying transport of an earthquake-induced landslide, Landslides, 18, 1612–5118, https://doi.org/10.1007/s10346-020-01508-1, 2020. a
Haug, Ø. T., Rosenau, M., Leever, K., and Oncken, O.: Modelling Fragmentation in Rock Avalanches, Landslide Science for a Safer Geo-Environment, edited by: Sassa, K., Canuti, P., and Yin, Y., in: Landslide Science for a Safer Geoenvironment, Springer, Cham., 2, https://doi.org/10.1007/978-3-319-05050-8_16, 2014. a
Haug, Ø. T., Rosenau, M., Leever, K., and Oncken, O.: On the Energy Budgets of Fragmenting Rockfalls and Rockslides: Insights from Experiments, J. Geophys. Res.-Earth, 121, 1310–1327, https://doi.org/10.1002/2014JF003406, 2016. a, b, c, d, e, f, g, h, i, j, k, l
Haug, Ø. T., Rosenau, M., Rudolf, M., Leever, K., and Oncken, O.: Laboratory model data from experiments on fragmenting analogue rock avalanches, https://doi.org/10.5880/GFZ.4.1.2020.004, 2021. a, b, c, d, e, f
Heim, A.: Der Bergsturz von Elm, Zeitschrift der Deutschen Geologischen Gesellschaft, 34, 74–115, Schweizerbart Science Publishers, Stuttgart, Germany, 1882. a
Imre, B., Laue, J., and Springman, S. M.: Fractal fragmentation of rocks within sturzstroms: insight derived from physical experiments within the ETH geotechnical drum centrifuge, Granul. Matter, 12, 267–285, https://doi.org/10.1007/s10035-009-0163-1, 2010. a, b
Kent, P. E.: The transport mechanism in catastrophic rock falls, J. Geol., 74, 79–83, 1966. a
Lin, Q., Cheng, Q., Li, K., Xie, Y., and Wang, Y.: Contributions of Rock Mass Structure to the Emplacement of Fragmenting Rockfalls and Rockslides: Insights From Laboratory Experiments, J. Geophys. Res.-Earth, 125, e2019JB019296, https://doi.org/10.1029/2019JB019296, 2020. a, b, c
Lucas, A., Mangeney, A., and Ampuero, J. P.: Frictional velocity-weakening in landslides on Earth and on other planetary bodies, Nat. Commun., 5, 3417, https://doi.org/10.1038/ncomms4417, 2014. a, b, c, d, e
Manzella, I. and Labiouse, V.: Empirical and analytical analyses of laboratory granular flows to investigate rock avalanche propagation, 10, 23–36, Landslides, https://doi.org/10.1007/s10346-011-0313-5, 2012. a
Melosh, H. J.: Acoustic Fluidization : A New Geologic Process?, J. Geophys. Res., 84, 7513–7520, 1979. a
Pudasaini, S. P. and Fischer, J.-T.: A mechanical erosion model for two-phase mass flows, Int. J. Multiphas. Flow., 132, 103416, https://doi.org/10.1016/j.ijmultiphaseflow.2020.103416, 2020. a, b, c, d
Shreve, R. L.: Leakage and fluidization in air-layer lubricated avalanches, Geol. Soc. Am. Bull., 79, 653–658, https://doi.org/10.1130/0016-7606(1968)79[653:LAFIAL]2.0.CO;2, 1968. a
Wang, Y. F., Dong, J. J., and Cheng, Q. G.: Velocity-dependent frictional weakening of large rock avalanche basal facies: Implications for rock avalanche hypermobility?, J. Geophys. Res.-Sol., 122, 1648–1676, https://doi.org/10.1002/2016JB013624, 2017. a, b
Zhao, T., Crosta, G. B., Utili, S., and De Blasio, F. V.: Investigation of rock fragmentation during rockfalls and rock avalanches via 3D DEM analyses, J. Geophys. Res.-Earth, 122, 678–695, https://doi.org/10.1002/2016JF004060, 2017. a, b, c
Zhao, T., Crosta, G. B., Dattola, G., and Utili, S.: Dynamic Fragmentation of Jointed Rock Blocks During Rockslide-Avalanches: Insights From Discrete Element Analyses, J. Geophys. Res.-Sol., 123, 3250–3269, https://doi.org/10.1002/2017JB015210, 2018. a, b