Friction and slip measured at the bed of an Antarctic ice stream

The slip of glaciers over the underlying bed is the dominant mechanism governing the migration of ice from land into the oceans, contributing to sea-level rise. Yet glacier slip remains poorly understood or constrained by observations. Here we observe both frictional shear-stress and slip at the bed of an ice stream, using 100,000 repetitive stick-slip icequakes from Rutford Ice Stream, Antarctica. Basal shear-stresses and slip-rates vary from 10^4 to 10^7 Pa and 0.2 to 1.5 m day^(-1), respectively. Friction and slip vary temporally over the order of hours and spatially over 10s of meters, caused by corresponding variations in ice-bed interface material and effective-normal-stress. Our findings also suggest that the bed is substantially more complex than currently assumed in ice stream models and that basal effective-normal-stresses may be significantly higher than previously thought. The observations also provide previously unresolved constraint of the basal boundary conditions of ice dynamics models. This is critical for constraining the primary contribution of ice mass loss in Antarctica, and hence the endeavour to reduce uncertainty in sea-level rise projections.

Strike-slip overprinting of initial co-axial shortening within the toe region of a submarine landslide: a case study from the Angoche Basin, offshore Mozambique

Submarine landslides (slides) are some of the most voluminous sediment gravity-flows on Earth and they dominate the stratigraphic record of many sedimentary basins. Their general kinematics and internal structure are relatively well-understood. However, how slides increase in volume and internally deform as they evolve, and how these processes relate, in time and space, to the growth of their basal (shear) zone, are poorly understood. We here use three high-resolution 3D seismic surveys from the Angoche Basin, offshore Mozambique to map strain within a shallowly buried, large, and thus seismically well-imaged slide (c. 530 km3). We document several key kinematic indicators, including broadly NW-trending lateral margins and longitudinal shears bounding and within the slide body, respectively, and broadly NE-trending symmetric pop-up blocks in the slide toe. Approximately 7 km downdip of the slide toe wall, thrusts and related folds also occur within otherwise undeformed slope material, with thrusts detaching downwards onto the downslope continuation of the basal shear zone underlying the slide body. Based on the style, trend, and distribution of these features, and their cross-cutting relationships, we propose an emplacement model involving two distinct phases of deformation: (i) bulk shortening, parallel to the overall SE-directed emplacement direction, with contractional shear strains reaching c. 8%; and (ii) the development of broadly emplacement direction-parallel shear zones that offset the earlier-formed shortening structures. We infer that the contractional strains basinward of the slide body formed due to cryptic basinward propagation of the basal shear zone ahead of and to accommodate updip sliding and shortening associated with, the entire slide mass. Our study demonstrates the value of using 3D seismic reflection data to reveal slide emplacement kinematics, especially the multiphase, non-coaxial nature of deformation, and the dynamics of basal shear zone growth.

Rock Glacier Dynamics by a Thermo-Elastic-Viscoplastic Constitutive Relationship

As a result of mountain permafrost creep, rock glaciers are common features in high-altitude periglacial areas. From a practical point of view, beyond their localization and inventorying, both the monitoring and prediction of their evolution due to climate changes are crucial. One of the effects of climate change is the thickening of the basal shear zone (the portion of the rock glacier where most deformations are localized), eventually leading to the development of unexpected and unprecedented (in terms of location, magnitude, frequency, and timing) instability phenomena. These phenomena bear consequences for the understanding of landscape evolution, natural hazards, and the safe and sustainable operation of high-mountain infrastructures. Most of the studies about active rock glaciers are focused on the analysis of monitoring data, while just a few studies are focused on modeling their behavior to understand their possible further evolution. The active rock glacier response is characterized by a viscous (rate-dependent) behavior, influenced by seasonal temperature oscillations, and characterized by a seasonal transition from slow to fast. In this work, a new thermo-mechanical model based on the delayed plasticity theory and calibrated on experimental results is proposed. The model is employed to evaluate the influence of geometry and forcing (air temperature) on a real rock glacier (Murtèl-Corvatsch rock glacier) creep behavior.

Geographic-information-system-based topographic reconstruction and geomechanical modelling of the Köfels rockslide

The Köfels rockslide in the Ötztal Valley (Tyrol, Austria) represents the largest known extremely rapid landslide in metamorphic rock masses in the Alps. Although many hypotheses for the trigger were discussed in the past, until now no scientifically proven trigger factor has been identified. This study provides new data about the (i) pre-failure and failure topography, (ii) failure volume and porosity of the sliding mass, and (iii) numerical models on initial deformation and failure mechanism, as well as shear strength properties of the basal shear zone obtained by back-calculations. Geographic information system (GIS) methods were used to reconstruct the slope topographies before, during and after the event. Comparing the resulting digital terrain models leads to volume estimates of the failure and deposition masses of 3100 and 4000 million m3, respectively, and a sliding mass porosity of 26 %. For the 2D numerical investigation the distinct element method was applied to study the geomechanical characteristics of the initial failure process (i.e. model runs without a basal shear zone) and to determine the shear strength properties of the reconstructed basal shear zone. Based on numerous model runs by varying the block and joint input parameters, the failure process of the rock slope could be plausibly reconstructed; however, the exact geometry of the rockslide, especially in view of thickness, could not be fully reproduced. Our results suggest that both failure of rock blocks and shearing along dipping joints moderately to the east were responsible for the formation or the rockslide. The progressive failure process may have taken place by fracturing and loosening of the rock mass, advancing from shallow to deep-seated zones, especially by the development of internal shear zones, as well as localized domains of increased block failure. The simulations further highlighted the importance of considering the dominant structural features of the rock mass. Considering back-calculations of the strength properties, i.e. the friction angle of the basal shear zone, the results indicated that under no groundwater flow conditions, an exceptionally low friction angle of 21 to 24∘ or below is required to promote failure, depending on how much internal shearing of the sliding mass is allowed. Model runs considering groundwater flow resulted in approximately 6∘ higher back-calculated critical friction angles ranging from 27 to 30∘. Such low friction angles of the basal failure zone are unexpected from a rock mechanical perspective for this strong rock, and groundwater flow, even if high water pressures are assumed, may not be able to trigger this rockslide. In addition, the rock mass properties needed to induce failure in the model runs if no basal shear zone was implemented are significantly lower than those which would be obtained by classical rock mechanical considerations. Additional conditioning and triggering factors such as the impact of earthquakes acting as precursors for progressive rock mass weakening may have been involved in causing this gigantic rockslide.

A Structural Design Comparison Between Two Reinforced Concrete Regular 6-Level Buildings using Soil-Structure Interaction in Linear Range

Structural engineers commonly design superstructures as fixed at the base and transmit the reactions to the infrastructure in order to design the foundation system and estimate the displacement of the soil while disregarding the change in seismic response that this induces. In this article, the foundation system was transformed into equivalent springs, and the seismic response in the linear range was compared and quantified, obtaining results such as increased periods, increased amounts of steel reinforcement in beams (between 7% and 25%) and columns (between 29% and 39%), an increase in the number of stirrups per linear meter (between 3% and 11% in columns and between 5% and 45% in beams) and drifts (between 1% and 14%), and a decrease in basal shear (up to 20%), which directly affects the design of the structure. This study concludes that the inclusion of the soil-structure interaction is necessary for structural design in the linear range.

Microscale evidence of the role of water during the emplacement of mass flows in glacial environments

Microscale analysis of unlithified glacial soils can provide far greater detail regarding their depositional and deformation histories than can be obtained from macroscale studies alone. This paper presents the results of three detailed case studies which examine the processes occurring during overriding and emplacement of mass flows in glacial environments: (i) laminated soils deposited in a proglacial lake setting at Heinabergsjökull, Iceland; (ii) a channelised, ice-marginal to submarginal mass flow at Whitburn (County Durham), England; and (iii) a mass flow exposed at Carstairs, Central Scotland which was emplaced in a glaciofluvial to glaciolacustrine setting. Microscale evidence from all three sites is combined to develop a conceptual model of the role played by water during mass flow; from the initial soil disruption under and/or in front of an advancing mass flow, to the formation of a basal shear zone facilitating mass flow transport and emplacement, through to the decoupling of the flow from the underlying substrate as a result of the injection of fluidised soil along its base. The development of these detachments during mass flow has the potential to increase the velocity of the flow and its run out distance, increasing the potential impact of these geohazards.Thematic collection: This article is part of the Role of water in destabilizing slopes collection available at: https://www.lyellcollection.org/cc/Role-of-water-in-destabilizing-slopes

Basal Pressure Variations Induced by a Turbulent Flow Over a Wavy Surface

Turbulent flows over wavy surfaces give rise to the formation of ripples, dunes and other natural bedforms. To predict how much sediment these flows transport, research has focused mainly on basal shear stress, which peaks upstream of the highest topography, and has largely ignored the corresponding pressure variations. In this article, we reanalyze old literature data, as well as more recent wind tunnel results, to shed a new light on pressure induced by a turbulent flow on a sinusoidal surface. While the Bernoulli effect increases the velocity above crests and reduces it in troughs, pressure exhibits variations that lag behind the topography. We extract the in-phase and in-quadrature components from streamwise pressure profiles and compare them to hydrodynamic predictions calibrated on shear stress data.

Generalized sliding law applied to the surge dynamics of Shisper Glacier and constrained by timeseries correlation of optical satellite images

Understanding fast ice flow is key to assess the future of glaciers. Fast ice flow is controlled by sliding at the bed, yet that sliding is poorly understood. A growing number of studies show that the relationship between sliding and basal shear stress transitions from an initially rate-strengthening behavior to a rate-independent or rate-weakening behavior. Studies that have tested a glacier sliding law with data remain rare. Surging glaciers, as we show in this study, can be used as a natural laboratory to inform sliding laws because a single glacier shows extreme velocity variations at a sub-annual timescale. The present study has two parts: (1) we introduce a new workflow to produce velocity maps with a high spatio-temporal resolution from remote sensing data combining Sentinel-2 and Landsat 8 and use the results to describe the recent surge of Shisper glacier, and (2) we present a generalized sliding law and provide a first-order assessment of the sliding-law parameters using the remote sensing dataset. The quality and spatio-temporal resolution of the velocity timeseries allow us to identify a gradual amplification of spring speed-up velocities in the two years leading up to the surge that started by the end of 2017. We also find that surface velocity patterns during the surge can be decomposed in three main phases, and each phase appears to be associated with hydraulic changes. Using this dataset, we are able to constrain the sliding law parameter range necessary to encompass the sliding behavior of Shisper glacier, before and during the surge. We document a transition from rate-strengthening to rate-independent or rate-weakening behavior. A range of parameters is probably necessary to describe sliding at a single glacier. The approach used in this study could be applied to many other sites in order to better constrain glacier sliding in various climatic and geographic settings.