Search

Water infiltration into fractures is ubiquitous in crustal rocks. However, little is known about how such a progressive wetting process affects fracture stiffness and seismic wave propagation, which are highly relevant for characterizing fracture systems in situ. We study the acousto‐mechanical behavior of a free‐standing fractured granite subjected to gradual water infiltration with a downward‐moving wetting front over 12 days. We observe significant differences (i.e., by an order of magnitude) in wave amplitudes across the fractured granite compared to an intact granite, with both cases showing a strong correlation between wave amplitudes and wetting front movement. Effects of water infiltration into the fracture and surrounding matrix on seismic attenuation are captured by a numerical model with parameters constrained by experimental data. Back‐calculated fracture stiffness decreases exponentially with the wetting front migration along the fracture. We propose that moisture‐induced matrix expansion around the fracture increases asperity mismatch, leading to reduced fracture stiffness. Plain Language Summary: In the shallow layers of the Earth, hydrological cycles such as snowmelt, fog, dew, and rain have been shown to change the moisture content of crustal rocks, which can alter the elastic properties of natural fractures and affect the propagation of seismic waves. Understanding how seismic waves propagate in the near‐surface environment is crucial for the assessment of earthquake hazards and the characterization of geologic heterogeneities. In this work, we perform well‐controlled laboratory experiments to study the acousto‐mechanical behavior of a single fracture in granitic rock subjected to progressive wetting over 12 days. We report that the fracture stiffness decreases exponentially as the wetting front advances along the fracture. Our research sheds light on an important question in fracture characterization: how elastic waves propagate across a fracture undergoing moisture‐induced expansion. Key Points: A laboratory study establishes a relationship among water imbibition, seismic attenuation, and stiffness evolution in a wetted fractureWave amplitudes across a fracture correlate strongly with the wetting front movement of infiltrated water within the fractureFracture stiffness exponentially decreases with the advance of the wetting front along the fracture [ABSTRACT FROM AUTHOR]

Rock slope failures often result from progressive rock mass damage which accumulates over long timescales. In deglaciating environments, rock slopes are affected by stress perturbations driven by mechanical unloading due to ice downwasting and concurrent changes in thermal and hydraulic boundary conditions. Since in‐situ data are rare, the different processes and their relative contribution to slope damage remain poorly understood. Here, we analyze borehole monitoring data from a rock slope adjacent to the retreating Great Aletsch Glacier (Switzerland) and compare it to englacial water levels, climate data, and decreasing ice levels. Rock slope pore pressures show a seasonal signal controlled by infiltration events as well as effects from the connectivity to the englacial hydrological system. We find that reversible and irreversible strains are driven by: (a) hydromechanical effects caused by englacial pressure fluctuations and infiltration events, (b) stress transfer related to changing mechanical glacial loads from short‐term englacial water level fluctuations and longer term ice downwasting, and (c) thermomechanical effects from annual temperature cycles penetrating the shallow subsurface, which primarily result in reversible deformation. We relate most observed irreversible strain (damage) to mechanical unloading from ice downwasting. Damage is strongest directly at the ice margin and moves through the slope at the pace of glacial retreat and advance. Locations with many retreat/advance cycles are very sensitive to landslide formation. The current climate warming impacts very sensitive valley sectors, which is confirmed by landslide distributions and activity in the study area. Plain Language Summary: The formation of rock slope instabilities is a long‐term process related to continuous rock mass weakening. For rock slopes in glaciated valleys, increased rates of rock mass damage are assumed to occur during ice retreat due to a reduction of the glacier weight and concurrent changes in rock temperatures and slope groundwater pressures. However, because direct field measurements are missing, governing processes and their relative importance are still poorly understood. Here, we analyze borehole monitoring data and compare it to englacial water levels, climate data, and decreasing ice levels. We find that the glacial hydrology influences the dynamics of the groundwater pressures in the adjacent rock slope at various timescales. We record micrometer‐scale reversible (elastic) deformation and irreversible rock mass damage and identify the involved drivers. Several processes contribute to progressive rock mass weakening at different timescales and magnitudes, but most damage is related to mechanical unloading from ice downwasting and is strongest directly at the ice margin. Locations with many retreat/advance cycles are sensitive to increased damage accumulation and landslide formation. The current climate warming impacts very sensitive valley sectors, which is confirmed by landslide distributions and activity in the study area. Key Points: We monitored subsurface pore pressures and micrometer scale strain in a rock slope during strong glacial retreat caused by climate warmingWe identified thermo‐hydro‐mechanical drivers for reversible deformation and irreversible rock mass damage acting at various timescalesSlope damage intensity can be related to landslide spatial and temporal distribution [ABSTRACT FROM AUTHOR]

Cycles of glaciation alter the temperature field in proximal alpine valley flanks, driving rock slope damage through thermomechanical stresses. Here we extend simplified assumptions of glacial debuttressing to quantitatively examine how paraglacial bedrock temperature changes, acting in concert with changing ice loads during Late Pleistocene and Holocene glacial cycles, create damage in adjacent rock slopes and prepare future slope instabilities. When in contact with temperate glacier ice, valley walls maintain near isothermal ~0 °C surface temperatures and are shielded from daily and annual cycles. With retreat, rock walls are rapidly exposed to strongly varying temperature boundary conditions, a transition we term "paraglacial thermal shock." Using detailed, conceptual numerical models based on the Aletsch Glacier in Switzerland, we show that including thermomechanical stresses during simulated glacial cycles creates significantly more rock slope damage than predicted for purely mechanical ice loading and unloading. Glacier advances are especially effective in generating damage as rapid cooling drives contraction of the rock mass reducing joint normal stresses. First time exposure to annual temperature cycles during deglaciation induces a shallow damage front that follows the retreating ice margin, generating damage in a complementary process at shorter time scales. Acting on a reduced strength rock mass, modeled thermomechanical cycles enhance the development of a slope instability with similar attributes as observed in our study area. Our results demonstrate that thermomechanical stresses acting in conjunction with changing ice loads are capable of generating considerable rock slope damage with spatial and temporal patterns controlled by glacier extents. Key Points: Subsurface temperature changes on glacial time scales can generate rock mass damage through thermomechanical coupled stressesGlacier advances are most effective in producing rock slope damage as cooling drives contraction of subglacial bedrockThermal shock during deglaciation creates a shallow damage front as bedrock is first exposed to annual temperature cycles [ABSTRACT FROM AUTHOR]

Within the framework of the Gotthard Base Tunnel Project in the Central Alps, Switzerland, geodetic monitoring networks were installed above the tunnel trajectory in alpine valleys. Natural ground-surface deformation recorded in the years prior to the tunneling excavation was seen to contain an unexpectedly large cyclical component of horizontal strain across the valleys, which was seasonal and appeared to be due to elastic processes. The strain is strongly correlated with snow melt and rainstorm precipitation, suggesting the implied rock-mass deformation is driven by changes in water-table elevation within adjacent mountain slopes. The horizontal strains are of the order of 1-2 · 10, which is close to the design limits that can be accommodated by hydropower arch dams in the study area. This study investigates these processes in detail and describes a new mathematical model (REROD), which is able to accurately reproduce and predict such natural rock-slope displacements. The model implements a transfer-function approach to predict the valley-crossing strains from rainfall and winter snow height data recorded at nearby meteorological stations. It has been used to estimate and remove the natural strain signal from the net recorded deformation so as to resolve the component due to tunneling. [ABSTRACT FROM AUTHOR]

Investigating surface displacements in high alpine environments is often subject to challenges due to the difficult accessibility or harsh climatic conditions. Measurement systems have improved greatly in recent years regarding accuracy, range, or energy consumption. Continuously receiving high-precision, real-time monitoring data from a remote location can still support a better understanding of slope dynamics and risk. We present the design, construction, operation, and performance of a complex surface displacement monitoring system installed in the surroundings of the Great Aletsch Glacier in the Swiss Alps, based on two robotic total stations to continuously measure 3D displacements with high accuracies. In addition, GNSS stations are also considered in order to pass from a local to a geographic reference system, as well as to improve the measurement accuracy. The monitoring network is aimed at studying several types of deformation processes, i.e., (i) gravitationally driven and irreversible rockslide movements around the tongue of the Great Aletsch Glacier, (ii) reversible rock slope deformations caused by annual cycles of groundwater recharge and depletion, and (iii) small irreversible deformations of stable rock slopes resulting from progressive rock damage driven by glacier retreat and cyclic hydraulic and thermal loading. We describe the technical details of the monitoring system, which has been in operation successfully for 6 years, and discuss the system performance in terms of its robustness and accuracy. [ABSTRACT FROM AUTHOR]

Abstract: In 1991, the steep slope above the village of Randa in southern Switzerland failed during two successive rockslide events three weeks apart. The failures, which resulted in the release of approximately 30×106 m3 of material, occurred in fractured gneissic rock whose principal foliation, faults and fracture zones dip into the slope. Above the scarp of the second rockslide, the mountain continues to move at rates of up to 2 cm/year. In an attempt to understand failure mechanisms in fractured crystalline rock masses, we have carried out a comprehensive multidisciplinary study that included detailed geological mapping, borehole logging and imaging, and surface and borehole geophysical surveys. An integrated interpretation of the resultant data has revealed the presence of three prominent fault- and fracture zone sets that define the internal structure of the rock mass. Some faults and fracture zones could be traced up to 85 m into the subsurface using borehole radar reflection methods. The dips of the faults and fracture zones estimated from surface mapping matched those from radar reflection and borehole televiewer surveys. No evidence was found for significant faults or persistent weakness planes that paralleled the 1991 rockslide failure surfaces. Instead, the new observations were consistent with dominant failure mechanisms that involve the progressive development of major shear and brittle fracture zones and the eventual destruction of intact rock bridges immediately prior to catastrophic collapse. This integrated study provides the first reliable model of the 3-D fault structure associated with and contributing to a large rock slope instability in a brittle crystalline rock mass. [Copyright &y& Elsevier]

Crystalline rock has been tested/selected by many countries for developing underground nuclear waste repositories at ~500 m depth to achieve geological isolation of high-level, long-lived radioactive waste. During the assessment period of up to one million years, large earthquakes may occur around the repository and trigger coseismic displacements along secondary fractures, jeopardising the integrity of the buffer-waste canister system. It is, therefore, of great importance to understand the coseismic behaviour of the repository site during large earthquakes. In this paper, we develop a finite element method-based seismo-mechanical model to simulate the response of fractured rocks subject to both in-situ stresses and seismic activities. The model can well capture the fracture displacement behaviour under dynamic loadings involving alternating regimes of contact loss, partial slip and total sliding. We model the earthquake-induced displacement field according to the seismic source theory in combination with a generic source time function. We apply the model to the nuclear waste repository site at Forsmark, Sweden, and analyse the coseismic responses of both single fractures and fracture networks during a potential post-glacial earthquake with a moment magnitude of M w = 5.6. The shear dislocation of a single fracture is strongly dependent on the fracture length and dip angle, while the displacement pattern of a fracture network is dominated by its "backbone" structures. We observe that significant coseismic shear displacement occurs if the fractured rock is close to the hypocentre and located in the dilational quadrant of the primary fault due to reduced shear resistance. However, the earthquake-induced shear displacement decreases drastically with the increased distance to the hypocentre and a distance of ~700 m may be needed for a fracture up to 100 m long to not displace beyond 5 cm. If two repeated earthquakes would occur, a distance of ~1200 m may be necessary. [ABSTRACT FROM AUTHOR]

The Moosfluh slope, located in the vicinity of the current tongue of the Great Aletsch glacier, Switzerland, hosts one of the largest active instabilities of the European Alps. Despite at this location accelerated slope displacement has been detected over the past decade, the rapid evolution observed between September and November 2016, locally reaching tens of meters, was largely unexpected. During this "crisis phase", the rock mass experienced substantial internal deformation mainly composed of toppling, formation of tensile scarps, and basal sliding. The large internal deformation caused also several local failure events in the form of single block falls and/or rock mass collapses of moderate size. Here we focus on the post-crisis, i.e. the period going between spring 2017 until today. The deformation occurring at the Moosfluh slope is still consistent and rock failure events are regularly going on, however, not in a dramatic fashion as observed in fall 2016. We show the results obtained by jointly investigating satellite, airborne, and ground based imagery, exploiting optical and radar sensors, as well as the data recorded from a local seismic network. In particular, we focus on the joint analysis of the information obtained from the optical imagery acquired with high temporal frequencies and the seismic waveforms associated to rock failure events. The combined remote sensing and geophysical dataset presented here is unique, and allows for a detailed interpretation of the rapid spatial and temporal rock mass disintegration in an alpine slope. Our results are of interest to understand mutual influences between predisposing slope instability factors (i.e., mainly due to geological, structural, and geomorphological conditions), and external driving factors (i.e., seismicity and meteo-climatic effects) on the progressive evolution of rock mass, specifically at the Moosfluh slope but also in similar case scenarios. [ABSTRACT FROM AUTHOR]

Fractured rock-masses deform in response to changes in pore pressure. In alpine catchments, groundwater flow dynamics is generally dominated by strong recharge signal from snowmelt infiltration in spring / early summer and in early winter, and punctuated by large rainstorms in snow-free periods. In crystalline rock slopes, infiltration occurs along discontinuities, such as faults and joints. Under constant normal stress joint apertures increase non-linearly during pore pressure increase, increasing the volume of saturated rock masses and leading to surface deformation. The latter can be measured with high-resolution geodetic monitoring systems and remote sensing methods such as DInSAR.In this work, we discuss temporal and spatial variations of reversible ground surface deformation measured in the Aletsch region, Switzerland. The study area is situated in the vicinity of the tongue of the Great Aletsch Glacier, central Swiss Alps, in the gneissic and granitic bedrock of the Aar massif. The valley (1550 to 2300 m.a.s.l.), is roughly oriented parallel to the main Alpine foliation (NE-SW) and to major fault zones. Total annual precipitation is less than 1 m/y, but snow depth frequently reaches 2 m during winter, with strong spatial variability. We use a network of 4 continuous Global Positioning System (cGPS) stations and 90 reflectors linked to 2 total stations (TPS), operating since 2014, to track deformation at millimetric to centimetric scale. We use weather and hydrological data from surrounding stations to estimate the timing of infiltration.At some locations in our study area we observe peak-to-peak magnitude of seasonal deformation over 2 cm, whereas other areas show no significant seasonal deformation. Measured displacement vectors are generally oriented towards the valley center in spring, when deformation rates are the highest (>0.5 mm/d). The start of the high-rate of deformation period also varies spatially. The South-facing slope starts deforming around 20 days earlier than the North-facing slope. The vertical component of motion is generally smaller than the horizontal one, and often below the noise level. Finally, we discuss possible factors (e.g. topography and slope orientation, discontinuities' orientation and density, lithology, groundwater recharge magnitudes and durations) that might influence the spatial and temporal variations of reversible deformation on the slopes adjacent to the tongue of the Great Aletsch Glacier. [ABSTRACT FROM AUTHOR]

The Gotthard Base Tunnel (GBT), constructed between 2000 and 2011, is a 57 km long and up to 2.5 km deep high-speed railway tunnel located in the Swiss Alps. Significant ground surface displacements reaching about 10 cm were observed during and after the tunnel construction. To gain a better understanding of the causal mechanism of such conspicuous ground deformations, we develop a three-dimensional (3D) fully-coupled hydro-mechanical model to simulate GBT-induced groundwater drainage, stress redistribution, rock mass consolidation, and fault zone deformation at the regional scale. First, we construct a geological model with topographical features, lithological units, and natural faults realistically represented. We constrain the material properties of fault zones and rock masses based on available extensive laboratory testing results and site characterisation datasets. We then simulate the tunnelling process over time, with the resulting coupled hydro-mechanical responses of faulted rock masses well captured and ground surface/subsurface displacements quantitatively analysed. The simulation results in general show a good agreement with the field monitoring data of ground surface displacement, subsurface tunnel settlement, and groundwater inflow into the tunnel. Our model indicates that ground surface displacements originate from GBT-induced water drainage and rock mass consolidation in the deep subsurface. Our results also show that the GBT construction could trigger faults to shear via drainage-induced pressure diffusion and poroelastic stressing. The research findings from our work have important implications for many groundwater drainage-related geoengineering activities such as underground excavation in alpine mountains and fluid withdrawal in subsurface reservoirs. [Display omitted] • 3D simulation captures tunnelling-induced groundwater flow and rock mass deformation. • The simulation explains ground subsidence in crystalline rock during deep tunnelling. • Drainage and consolidation in deep subsurface cause ground surface displacements. • Tunnelling-induced pressure diffusion and poroelastic stressing drive fault movement. [ABSTRACT FROM AUTHOR]

In this work we present a methodology based on remote sensing data to map and classify unstable slopes in high alpine areas and apply this methodology to the Himalaya of northwestern Bhutan. We leverage on radar images acquired by the ESA Envisat and the JAXA ALOS-1 satellites between 2007 and 2011. Synthetic aperture radar differential interferometry (DInSAR) is used to identify and quantify recent ground surface displacements associated with potentially unstable slopes. Optical satellite images and a high-resolution digital surface model is used for the geomorphological classification of individual landforms and deposits associated with displacements. The analysis of the geomorphological characteristics reveals that displacements are associated with: rock slides, mountain slope deformation and rock slope deformation (77.5%), soil creep (4%), soil slides (1.9%) and rock glaciers (12.3%); 4.3% of the detected displacements remain of unknown origin. The method applied here for the assessment of activity likelihood uses a very large number of interferograms obtained with images acquired by different satellites. The method, which is new and replicable, takes advantage of the combination of different acquisition geometries and then combines indices of vegetation cover, layover, velocity estimate and temporal sampling associated with each detected active area to assign a weight. The analysis shows that a cluster of unstable slopes exists in the northwest, potentially related to the presence of an important regional lineament and the distribution of permafrost. Comparatively little ongoing activity is detected in the south, but field validation documents substantial landslide activity in the recent geological past. The dataset here generated, containing 700 potentially unstable areas with an associated activity likelihood, is a novel dataset that makes an essential basis of future analysis to investigate rock slope instability distribution and predisposing factors. Moreover, this dataset represents a new important step towards future landslide hazard assessment in the Himalaya of Bhutan. • 693 previously unknown unstable slopes in NW Bhutan were found with DInSAR • Geomorphological classification revealed that almost 78% are rock slope instabilities • The spatial distribution of potentially unstable slopes in the region is clearly non uniform • Activity of unstable slopes is likely influenced by regional tectonic structures • Unstable slopes seem to be less abundant in the region in recent times than in the past [ABSTRACT FROM AUTHOR]

Abstract: The monitoring of slope displacements over time provides the basis for most rockslide early warning systems, yet the prediction of catastrophic failure from these records is highly problematic and tenuous, especially if the underlying kinematics and instability mechanism are poorly understood. An example is the moving slope above the town of Randa in the Swiss Alps. This slope is considered typical of those in crystalline rock that lacks a natural, highly persistent, weakness plane dipping towards the valley that can serve as a through-going detachment surface for kinematic release. In Part I (the companion paper to this), the findings from a comprehensive geological and geophysical investigation to image the internal structures of the unstable rock mass were presented. In this paper we develop a kinematic model that describes the pattern of displacement vectors for the rock mass along these structures, both on surface and at depth. The displacements were estimated from 5 years of data from an extensive monitoring system which included surface geodetic and crackmeter measurements, borehole inclinometer and extensometer measurements (up to 120 m depth), and microseismicity. The results showed that the displacement field is highly heterogeneous. Internal deformation is accommodated by both shear and opening-mode dislocation of faults and fracture zones which dip moderately to steeply into the slope. Microseismicity is most intense near the front of the scarp where surface translation in large part reflects the toppling of blocks accommodated by slip along the steeper dipping fractures. The small displacement rates at the study site of up to 2 cm/year coupled with a modest deviation of up to 15° from vertical of the boreholes posed severe problems for the estimation of horizontal displacements from the inclinometer data. We describe the error analysis of these data in some detail since it is relevant for similar installations elsewhere. Conclusions drawn from this work highlight the importance of integrating various types of data in order to better understand the complex block-kinematic processes whose evolution governs the long-term progressive failure of unstable rock slopes. [Copyright &y& Elsevier]

In February 2017, a total of six in-situ hydraulic shearing (HS) experiments were conducted in the framework of the decameter scale In-situ Stimulation and Circulation (ISC) project in crystalline rocks at the Grimsel Test Site, Switzerland. During these experiments, three ductile (containing very few distinct fractures) and three brittle-ductile (containing multiple fractures) shear zones were the target for high-pressure fluid injections. The injected fluid volume was ~1 m³ per experiment. The experiment volume has an overburden of ~480 m and the maximum principal stress is dipping 40° towards East with a magnitude of 13 MPa. In this contribution, we present the observations and preliminary interpretations of the hydromechanical rock mass response during the different experiments and compare them to the target shear zone type. To monitor the hydromechanical response, six boreholes were drilled into the rock mass. Three of them are dedicated to pressure monitoring in seven intervals covering the target shear zones. In the other three boreholes, a total of sixty longitudinal Fiber-Bragg Grating strain sensors were installed. These sensors have a base length of one meter and are installed across intact rock, fractures or shear zones. During most of the HS experiments, the near-wellbore injectivity was successfully increased to a magnitude that is similar for all experiments. The same was observed for the near-wellbore transmissivity. For both parameters, a stronger increase was observed at the ductile shear zones, where the initial parameter magnitudes were smaller. The jacking pressure was reduced during most experiments and the induced slip dislocation at the injection point was in the range of 0.5 - 1.5 mm. We observed a pressure increase up to a radial distance of ~16 m from the injection well during all experiments. Further away from injection location the pressure monitoring data suggest poro-elastic effects. Pressure fronts that were expected to propagate along the target shear zone were only observed during two experiments. The strain monitoring revealed a zone of transient variable deformation (i.e., mix of compression and extension) close to the injection location and a compressive zone further away from the injection location. The extent of the variable zone seems to depend on the target shear zone type and is larger for the ductile shear zone. During the experiments, a meter scale Poisson effect was directly measured in the rock mass. Additionally, the formation of new fractures was detected.Our data highlight the complexity of the rock mass and fracture response during a decameter scale hydraulic stimulation and provide a basis for the development of future large scale hydromechanical-coupled numerical models of stimulation experiments. [ABSTRACT FROM AUTHOR]