# Creep enhancement and sliding in a temperate, hard-bedded alpine glacier

Juan-Pedro Roldán-Blasco, Adrien Gilbert, Luc Piard, Florent Gimbert, Christian Vincent, Olivier Gagliardini, Anuar Togaibekov, Andrea Walpersdorf, Nathan Maier

► To cite this version:

Juan-Pedro Roldán-Blasco, Adrien Gilbert, Luc Piard, Florent Gimbert, Christian Vincent, et al. Creep enhancement and sliding in a temperate, hard-bedded alpine glacier. The Cryosphere, 2025, 19 (1), pp.267-282. 10.5194/tc-19-267-2025. hal-04945328

HAL Id: hal-04945328

https://hal.science/hal-04945328v1

Submitted on 31 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Creep enhancement and sliding in a temperate, hard-bedded alpine glacier

Juan-Pedro Roldán-Blasco<fnref n="1" />, Adrien Gilbert<fnref n="1" />, Luc Piard<fnref n="1" />, Florent Gimbert<fnref n="1" />, Christian Vincent<fnref n="1" />, Olivier Gagliardini<fnref n="1" />, Anuar Togaibekov<fnref n="1" />,<fnref n="2" />, Andrea Walpersdorf<fnref n="2" />, and Nathan Maier<fnref n="1" />

**Correspondence:** Adrien Gilbert (adrien.gilbert@univ-grenoble-alpes.fr)

Received: 29 May 2024 – Discussion started: 1 July 2024

Revised: 30 October 2024 – Accepted: 18 November 2024 – Published: 22 January 2025

**Abstract.** Glacier internal deformation is usually described by Glen's flow law using two material parameters: the creep factor (A) and the flow law exponent (n). However, the values of these parameters and their spatial and temporal variability are rather uncertain due to the difficulty in quantifying internal strain and stress fields at natural scales. In this study, we combine 1-year-long continuous measurements of borehole inclinometry and surface velocity with three-dimensional full-Stokes ice flow modeling to infer ice rheologies and sliding velocities for the ablation zone of the Argentière Glacier, a temperate glacier in the French Alps. We demonstrate that the observed deformation rate profile has limited sensitivity to the flow law exponent (n) and instead mainly reflects an increase in the creep factor (A) with depth, with A departing from its surface value by up to a factor of 2.5 below 160 m depth. We interpret this creep factor enhancement as an effect of increasing interstitial water content with depth (from 0% to 1.3%), which results in an average value of A = 148 MPa⁻³ a⁻<fnref n="1" />. We further observe that internal ice deformation exhibits seasonal variability similar to that concerning surface velocity, indicating that the local basal sliding velocity exhibits no significant seasonal variation. We suggest that these changes in deformation rate are due to variations in the stress field, driven by contrasting changes in subglacial hydrology conditions between the sides and center of the glacier. Our study provides further evidence that borehole inclinometry, combined with full-Stokes flow modeling, allows for the constraining of both ice rheology and basal friction at scales that cannot be inferred from surface velocity measurements alone.

## 1 Introduction

Glacier dynamics depend on both internal deformation and basal sliding. Given the scarcity of direct observations of ice rheology and basal sliding speed at natural scales, ice flow models commonly use inverse methods to estimate the material parameters that best fit surface velocities (e.g., Arthern and Gudmundsson, 2010; Fürst et al., 2015; Mosbeux et al., 2016). However, the problem is largely undetermined due to limited knowledge of ice thicknesses, and it does not allow for the identification of model weaknesses as model errors are compensated for by material parameter adjustments. Independent and accurate estimates of ice material parameters, as well as subglacial basal-friction changes, are key to better representing glacier dynamics in models, particularly the fraction of surface velocities due to basal sliding versus that due to internal ice deformation.

Ice deformation is commonly assumed to follow Glen's flow law (Glen, 1955), with the creep parameter (A) being primarily dependent on ice temperature (Barnes et al., 1971; Weertman, 1983). This dependency has been studied extensively in both the laboratory and the field (Cuffey and Paterson, 2010) and can be reasonably well accounted for in models given a temperature profile. Typical values for ice viscosity and exponents for temperate glaciers and ice caps correspond to A ≈ 75 MPa⁻³ a⁻<fnref n="1" /> and n = 3, respectively (Cuffey and Paterson, 2010). The applicability of these values for describing glacier dynamics at natural scales, however, is still quite uncertain, mainly due to other controls coming into play, such as ice texture (orientation and microstructure) (Barnes et al., 1971; Goldsby and Kohlst-

<footnote n="1">¹IGE, Univ. Grenoble Alpes, CNRS, INRAE, IRD, Grenoble INP, 38000 Grenoble, France</footnote>

<footnote n="2">²ISTerre, Univ. Grenoble Alpes, CNRS, IRD, UGE, 38000 Grenoble, France</footnote>
edt, 2001; Chauve et al., 2024), impurities (Jones and Glen, 1969), and water content (Lliboutry and Duval, 1985; Duval, 1977). These controls are challenging to evaluate as they may vary quite extensively in time and space depending on stress and deformation conditions (Chauve et al., 2024; Rathmann and Lilien, 2022). This is especially true under temperate-ice conditions, which are hard to study in the laboratory, where water content is not only difficult to control but also hard to characterize in the field, requiring the drilling of ice cores or indirect geophysical measurements, which are often difficult to interpret. Geophysical methods include ground-penetrating radar observations (Young et al., 2021; Ogier et al., 2023), seismic tomography (Benjumea et al., 2003; Endres et al., 2009), and surface nuclear magnetic resonance (Legchenko et al., 2014). However, the first method is mainly sensitive to integrated water content, and its spatial variations are poorly retrieved. Analysis of seismic and radar wave velocities provides relevant information on relative changes in water content but remains more uncertain in absolute terms (Murray et al., 2007). We note that the recent use of fiber optics in boreholes provides accurate measurements of variations in seismic velocity with depth and might be a promising method for the in situ characterization of temperate-ice properties, although until now, it has mostly been applied to cold-ice settings (Fichtner et al., 2023; Booth et al., 2020). As temperate-ice deformation is expected to primarily control ice dynamics, including on ice sheets such as the Greenland Ice Sheet, where most of the deformation is concentrated in its basal temperate layer (Law et al., 2023), it is crucial to study it specifically.

Basal sliding on hard beds is also known to be a function of ice deformation as it is enhanced near the bed, thus involving the same sources of uncertainty as presented above, with the added complexity that subglacial hydrology also affects it. Water pockets, commonly referred to as cavities, can form on the lee side of bedrock bumps, reducing the apparent bed roughness, facilitating creep, and thus enhancing basal sliding (Lliboutry, 1959, 1968). Subglacial channels can also form under sufficient turbulent-induced melt, causing the opposite effect of reducing basal sliding by lowering the basal water pressure (Röthlisberger, 1972; Schoof, 2010). Ultimately, the resultant effect of subglacial hydrology on overall glacier dynamics depends on the type of subglacial hydrology system at play at any given time and location. Of fundamental importance is the ability to evaluate the spatiotemporal evolution of basal sliding to improve our representation of the evolution of subglacial hydrology and its link to overall glacier dynamics in glacier and ice sheet models.

Borehole inclinometry is a unique means for investigating ice deformation and basal sliding simultaneously and can be achieved either by measuring changes in borehole orientation through repeated surveys (e.g., Perutz, 1949; Shreve and Sharp, 1970; Raymond, 1971; Hooke, 1973; Hooke and Hanson, 1986; Hooke et al., 1992; Harper et al., 2001; Mar-

shall et al., 2002; Chandler et al., 2008) or through continuous englacial tiltmeter recordings (e.g., Gudmundsson et al., 1999; Lüthi et al., 2002; Willis et al., 2003; Amundson et al., 2006; Ryser et al., 2014; Keller and Blatter, 2012; Doyle et al., 2018; Lee et al., 2019; Maier et al., 2019, 2021). These observations provide deformation rate profiles that allow for in situ estimates of ice rheology and basal velocity, as well as their respective temporal variations, with basal velocity determined by integrating the deformation rate profile with depth and subtracting it from the surface velocity (Hooke et al., 1992; Maier et al., 2021). However, a major challenge in evaluating ice rheology based on these observations is the retrieval of stress fields against which deformation rates can be compared. Since stress variations with depth along the borehole cannot be measured, they must be estimated independently using a modeling approach, which, depending on model assumptions or prescribed boundary conditions, can introduce large uncertainties into the derived creep factor or flow exponent, especially with respect to valley glaciers (e.g., Harper et al., 2001; Chandler et al., 2008), where the highly three-dimensional geometry necessitates an evaluation of the full stress tensor. Even on large ice caps, strong spatial variations in the measured deformation rate between different boreholes have been shown to reflect a complex stress field influenced by spatial variability in bed friction (Ryser et al., 2014). This, together with variations in ice temperature, complicates the accurate assessment of ice rheology. This stress field complexity can also influence the temporal variations observed in deformation rates, due to either flow over changing bed topography (Maier et al., 2019) or hydrologically driven temporal changes in basal-drag patterns (Hooke et al., 1992; Willis et al., 2003).

In this paper, we infer internal rheological parameters and reconstruct basal velocity over time by combining continuous borehole inclinometry observations with full-Stokes three-dimensional modeling of the stress field throughout an entire melting season. The study focuses on the Argentière Glacier (temperate ice), located in the French Alps, which has been intensively monitored for several decades (e.g., with regard to mass balance, surface velocity and topography, sliding velocity, and bedrock geometry) and provides a unique and well-constrained environment for a natural-scale study of ice deformation. We first describe the study site and the measurement methods. We then analyze our observations in terms of material parameters using the flow model and, finally, provide the observed time series of both deformation and basal velocities. With this methodology, we identify a depth dependency of ice viscosity, which we attribute to changes in interstitial water content, as well as temporal changes in ice deformation, which we argue are due to subglacial-hydrology-driven changes in basal-friction conditions. These novel observations contribute to a better understanding of the complex interplay between basal sliding and internal ice deformation while providing constraints on ice rheology in a natural setting.
## 2 Field site and instrumentation

### 2.1 The Argentière Glacier

The Argentière Glacier is a temperate glacier located in the Mont Blanc range of the French Alps (45°10′N, 6°10′E). The glacier rests on hard bedrock (Vivian and Bocquet, 1973) and extends for 9 km within an altitude range of 1600 to 3400 m, separated by an icefall at 2300 m. Glacier dynamics have been continuously monitored since the 1970s, particularly with regard to their basal sliding velocities, thanks to direct access to the glacier bed ~ 700 m downstream of the drilling site (Vincent and Moreau, 2016; Gimbert et al., 2021a; Gilbert et al., 2022). Measurements are conducted using a cavitometer installed in a natural ice cavity, which records sliding velocities at a 30 min resolution. Surface velocities exhibit seasonal variability, with the lowest velocities observed in fall (October–November) and sustained high velocities observed during summer (June–August) (Vincent et al., 2022). The glacier also benefits from subglacial-runoff monitoring conducted by the power company Electricité d'Emosson SA, along with continuous records of yearly surface mass balance, topography, and velocity from the GLACIOCLIM monitoring program (https://glacioclim.osug.fr/, last access: 15 January 2025). Several campaigns of ground-penetrating radar measurements have also provided a good understanding of basal topography (Rabatel et al., 2018; Gimbert et al., 2021b), further extended by seismic investigations (Sergeant et al., 2020).

### 2.2 Borehole deformation instrumentation

Drilling operations took place between 12 and 14 September 2019. The drilling sites are located in the central part of the ablation area at an elevation of 2380 m, where the ice thickness along the center flow line is about 230–250 m and the bedrock forms an overdeepening, where most of the boreholes are located (Fig. 1a and Table 1). We used a custom-built hot water drill operating at 70 °C to drill boreholes 10 cm in diameter at an average speed of 60 m h⁻¹. The positions of the final completed and instrumented boreholes (BH1, BH2, BH3, BH4, and BH5) are indicated by the blue dots in Fig. 1a.

The deformation rate sensors, which we refer to as tiltmeters, consist of a high-end triaxial gravity sensor (SCL3300, Murata) and a triaxial magnetic sensor (LSM303, STMicroelectronics). The gravity sensors are used to determine the tilt (θ), the angle with respect to the vertical, with an estimated accuracy of 0.01°. Lab calibrations have shown that tilt readings above 45° become increasingly unreliable. The magnetic sensors do not provide reliable absolute measures of orientation relative to north because they are very sensitive to parasitic magnetic fields. The tiltmeters are grouped in chains of 20 and are more densely concentrated toward the bottom of the glacier (Fig. 1b). For each

<table><caption>**Table 1.** Summary of the borehole instrumentation and depths after installation. The instrumented depth corresponds to the depth of the deepest sensor. The bedrock depth was previously estimated using ground-penetrating radar measurements (Rabatel et al., 2018; Sergeant et al., 2020).</caption><thead><tr><th>Borehole</th><th>Tiltmeters (no.)</th><th>Borehole depth (m)</th><th>Instrumented depth (m)</th><th>Bedrock depth (m)</th></tr></thead><tbody><tr><td>BH1</td><td>18</td><td>208</td><td>190</td><td>253 ± 10</td></tr><tr><td>BH2</td><td>19</td><td>238</td><td>234</td><td>237 ± 20</td></tr><tr><td>BH3</td><td>17</td><td>216</td><td>174</td><td>235 ± 20</td></tr><tr><td>BH4</td><td>19</td><td>237</td><td>211</td><td>234 ± 20</td></tr><tr><td>BH5</td><td>17</td><td>194</td><td>190</td><td>234 ± 10</td></tr></tbody></table>

borehole (i), we refer to each tiltmeter (j) using the format BHi#j, starting with a j value of 1 for the deepest tiltmeter. All sensors acquired data every 30 min, but high-frequency noise was removed by smoothing the data using an exponential filter with a 1 d time window.

We estimate the initial shape of the boreholes using tilt and azimuth data averaged over the first month of the record (15 September to 15 October 2019). The estimated shape is used to correct for borehole and sensor depth errors resulting from the boreholes not being perfectly vertical. Although the actual shapes are three-dimensional, uncertainties in the azimuth measurements lead us to restrict our analysis of borehole shape to the estimated horizontal distance between each inclinometer and a vertical line starting at the surface (Fig. 1b). In Fig. 1b, we also show the final shape attained by the BH2 borehole before the chain of tiltmeters broke in October 2020.

Tiltmeter array performance varied across boreholes. Sensor arrays at BH1 and BH5 stopped working after a few days and provided no useful data. At BH2, BH3, and BH4, all sensors recorded data for more than a year until late October 2020. However, both BH3 and BH4 exhibit a tortuous shape in their deeper sections, suggesting problems during drilling that subsequently affected the quality of the measurements. We provide a summary of borehole and sensor-corrected depths in Table 1. Since BH2 reaches the bed and provides good data quality, which we attribute to the nearly vertical shape of the borehole, most of the results presented in this paper are based solely on BH2.

The tilt time series for each sensor in the three boreholes – BH2, BH3, and BH4 – can be found in Sect. S1 in the Supplement. Tilt data recorded during the first few months are affected by insufficient mechanical coupling between the tiltmeters and the ice, resulting in unstable signals until steady tilt occurs, which we interpret as the point at which the tiltmeters become well coupled to the ice (Fig. 1c–f). The duration of the transition between unstable and steady tilt change varies from sensor to sensor, ranging from about 2 months for deeper sensors, where borehole closure from creep is faster, to about 6 months for shallower sensors, where borehole clo-
<figure><img src="image_2.png" /><figcaption>**Figure 1.** (a) Map of the ablation area of the Argentière Glacier (EPSG:27572 projection), showing ice thickness (black contours) and instrument locations. (b) The estimated initial shape, drilled depth (black triangles), and instrumented depth (with each star representing a tiltmeter) of the five boreholes. The dashed orange line represents the estimated shape of BH2 at the end of the tilt record. (c–f) Tilt ($\theta$) recorded at four example inclinometers at BH2 (as labeled in panel b). The dashed vertical lines indicate 15 February 2020, the starting date of our analysis, when all sensors at BH2 had stabilized. GNSS: Global Navigation Satellite System.</figcaption></figure>

sure from creep takes longer. At BH2, most tiltmeters became fully coupled by 15 February 2020 – that is, 5 months after installation (see Sect. S1). We thus start our analysis from this point. Certain sensors reached a minimum tilt after installation, either shortly after installation (e.g., BH2#8) or long after (e.g., BH2#1 and BH2#14, with the minimum occurring in September 2019), indicating that they initially tilted against the flow.

## 2.3 Surface motion instrumentation

Five Global Navigation Satellite System (GNSS) stations were deployed in the ablation zone of the Argentière Glacier in February 2019, with an additional seven installed in Febru-

ary 2020. We employ multi-frequency Leica GR25 receivers and Leica AS10 antennas, which continuously record GNSS signals at a 1 Hz sampling interval. This GNSS receiver network (green dots in Fig. 1) covers the borehole sites, with the ARG1 station located ~ 30 m from BH2. See Togaibekov et al. (2024) for more details.

# 3 Methods

## 3.1 GNSS processing

GNSS signals are processed using a static approach with a double-difference processing technique and linear, ionosphere-free combination-phase observables (Bock et al.,
1986), incorporated into the software packages GAMIT and GLOBK (Herring et al., 2018). Daily GNSS position time series are then converted into horizontal velocities, which are used in this study. We empirically determined a velocity error of 0.9 m a⁻¹ using a stationary GNSS station located on bedrock approximately 500 m away from the survey network, which is exposed to a similar multi-path scattering environment as the stations on the glacier.

## 3.2 Calculating internal-deformation rates from tiltmeter observations

We use a three-dimensional reference system with x as the main along-flow direction and z as the upward vertical, with its origin at the surface (see Fig. 1a and b). Velocities in the x, y, and z directions are denoted by u, v, and w, respectively. Assuming that the temporal evolution of the tilt (θ) occurs entirely in the along-flow direction, the changes in θ are controlled by the horizontal shear strain (du/dz), as well as the compressive strain (du/dx) and extensive strain (dw/dz) (Keller and Blatter, 2012). To detect the potential effect of compressive/extensive strain, we fit the tilt curves recorded at each BH2 inclinometer using the analytical model from Keller and Blatter (2012) (see Sect. S2). We find that the best data fit can always be obtained by neglecting du/dx and dw/dz, except with respect to the two deepest sensors (BH2#1 and BH2#2), where a better fit is obtained with nonzero compressive strain. For these deepest sensors, neglecting compressive strain leads to an overestimation of du/dz by about 30% at BH2#1 and 20% at BH2#2. The presence of significant compressive strain near the glacier bed is likely related to a local effect of bed roughness, as shown by the local evaluation of strain over a rough bed presented in Sect. 5.1. Hereafter, we refer to this layer as the basal boundary layer. Within this layer, the derived du/dz values are expected to be strongly biased and are thus ignored in the ice rheology interpretation.

With the hypothesis that du/dz dominates the flow gradient outside of the boundary layer, the internal-deformation rate (du/dz) from the temporal evolution of the tilt (θ) is computed as follows (Lüthi et al., 2002; Ryser et al., 2014; Doyle et al., 2018; Maier et al., 2019):

$$ \frac{du}{dz} = \frac{1}{dt} \frac{dx}{dz} \approx \frac{1}{\Delta t} \Delta \tan(\theta), \qquad (1) $$

where Δt is a given time period and Δtan(θ) is the change in the tangent of the tilt during that time period. In our particular implementation, we calculate the least-squares linear approximation of tan(θ) as a function of time for each time period (Δt), set to 1 d, meaning the regression coefficient directly yields du/dz. The estimated maximum uncertainty for the daily averaged deformation rate is, for most of the sensors, lower than 0.1 a⁻¹ (see Sect. S3).

## 3.3 Computation of time series of surface, internal, and basal velocity

As it is not always possible to directly interpolate the surface velocities obtained at borehole locations from the GNSS network due to data gaps, we construct a continuous time series of the surface velocities recorded at the GNSS station close to the boreholes by filling gaps using a linear model traditionally employed to derive continuous surface mass balances from sparse data on alpine glaciers (Lliboutry, 1974; Vincent et al., 2017). In this linear model, we assume similar temporal variability in surface velocity across the stations, meaning that the surface velocity at each GNSS station (i) can be expressed as

$$ u_{si}(t) = \alpha_i + \beta(t), \qquad (2) $$

where αi is the average surface velocity at the station (i) over the measurement period and β(t) represents the temporal variability, which is assumed to be identical for all stations and satisfies Σβ(t) = 0. We first solve the system of equations by finding the values of αi and β(t) that best approximate the observations while satisfying Σβ(t) = 0. We then compute the residuals between the modeled and observed velocities, along with their standard deviation (sres) values, and classify as outliers all observations with residuals greater than 3sres. We then solve the system again, excluding these outliers, to obtain the final values of αi and β(t). The residuals between the observed and reconstructed velocities follow a normal distribution centered around 0, which validates the initial assumption (Sect. S4). The surface velocity time series uS = αARG1 + β(t) (sres = 3.2 m a⁻¹) is thus used to fill the gaps at the ARG1 station, which is the closest station to BH2, as seen in Fig. 1.

The internal velocities are computed by integrating the deformation rate over depth as follows:

$$ u_d(z, t) = \int_{z_{\text{bed}}}^{z} \frac{du}{dz}(z, t)dz, \qquad (3) $$

where zbed is the bedrock elevation (expressed in meters). Finally, the basal velocity, uB(t), is computed as the difference between the reconstructed surface velocity, uS(t), and the integrated deformation rate over the entire ice thickness, uD(t). The time series are computed from the daily averaged deformation rate at a daily resolution.

## 3.4 Quantifying the parameters of Glen's flow law

In addition to deformation measurements, the study of ice rheology requires knowledge of the stress field within the ice, which is evaluated in this study using numerical modeling. We use the three-dimensional, full-Stokes finite-element model Elmer/Ice (Gagliardini et al., 2013) to solve the equation of conservation of momentum for a given glacier geometry and ice rheology. The glacier geometry is prescribed
based on bedrock topography and surface topography measurements derived from Pleiades satellite imagery obtained on 25 August 2019 (Beraud et al., 2022). The ice rheology is described by Glen's flow law:

$$ \dot{\epsilon}_{ij} = A \tau_e^{n-1} \tau_{ij}, \quad (4) $$

where $\dot{\epsilon}_{ij}$ and $\tau_{ij}$ are the components of the strain rate tensor (a$^{-1}$) and the deviatoric stress tensor (MPa), respectively; A is the creep factor (MPa$^{-n}$ a$^{-1}$); $\tau_e = \sqrt{\frac{1}{2}\tau_{ij}\tau_{ij}}$ is the effective stress (MPa); and n is an exponent. We assume a stress-free upper surface boundary condition and a basal boundary condition as defined by the Weertman friction law (Weertman, 1957),

$$ A_s \tau_b^m = u_b, \quad (5) $$

where $\tau_b$ is the basal shear stress (MPa), m is an exponent setting the nonlinearity of the law and is taken to be 3 (Gilbert et al., 2023), $A_s$ is the sliding coefficient at the bed (m a$^{-1}$ MPa$^{-m}$), and $u_b$ is the sliding velocity (m a$^{-1}$). We use values of $A_s$ for m = 3 as inferred by Gilbert et al. (2023) from surface velocity inversion measurements. To avoid stress anomalies caused by uncertainties in the ice thickness, we relax the surface topography for 1 year using surface mass balance forcing from Gilbert et al. (2023) before extracting the deformation rate and stress tensor from the model.

We run several simulations to test the sensitivity of the deformation rate profile to different values of A and n. We run a set of simulations with n = 3, 4, or 5 and a constant, uniform creep factor (A). The value of A for each n value is chosen such that the numerically calculated total deformation velocity at the location of BH2 matches the observed value. To minimize the potential effect of temporal changes in basal friction on the stress field (Hooke et al., 1992; Willis et al., 2003), we use an averaged deformation rate profile that was recorded during the last month of the time series (1 to 31 October 2020), when basal friction was least affected by subglacial hydrology and the inclinometers were well coupled to the ice. We run another set of simulations using n = 3 and a depth-dependent creep factor, A = A(z), ensuring that the computed du/dz profile matches the observed one. We infer A(z) using Glen's flow law (Eq. 4), based on the observed mean du/dz value at BH2 and the stress tensor computed from the Elmer/Ice simulation, such that

$$ A(z) = \frac{1}{2} \frac{du}{dz} \tau_{E, \text{num}}^{-2} \tau_{xz, \text{num}}^{-1} \quad (6) $$

We then approximate A(z) using a piecewise linear function, $A_{fit}(z)$. Since changing the creep factor slightly modifies the overall stress balance, we repeatedly run the numerical model, updating the $A_{fit}(z)$ value inferred from the numerical solution of the previous iteration during each iteration, until the modeled stress field converges. The depth-dependent

<figure><img src="image_3.png" /><figcaption>Figure 2. Average measured deformation rate profiles at BH2, BH3, and BH4, along with the monthly minima and maxima for BH2. The continuous lines show the average measured deformation rate profile for each borehole for the period between 15 February and 15 October 2020, and the shaded region represents the range between the monthly averaged minimum and maximum deformation rate values (shown only for BH2). Each symbol represents a tilt-meter.</figcaption></figure>

creep factor, $A_{fit}(z)$, determined at BH2 is applied uniformly over the entire domain by normalizing it with depth and applying this normalization everywhere.

# 4 Results

## 4.1 Observed mean deformation rate profiles

In Fig. 2, we show the average deformation rate profiles computed from tilt measurements (see Sect. 3.2 (Methods)) at BH2, BH3, and BH4 between 15 February and 15 October 2020. The orange-shaded region represents the range associated with the monthly averaged deformation rate (du/dz) for BH2. The retrieved profiles show a similar depth-increasing deformation rate, with only a few outliers, associated with the high deformation rate at ~ -160 m at BH3 and the low deformation rate at ~ -200 m at BH4. These outliers are likely measurement artifacts, as evident in the raw
data (see Sect. S1; BH3#6 and BH4#5). Only BH2 reached the bed (see Fig. 1b and Table 1), and it also generally shows lower measurement noise, likely due to better ice coupling (see Sect. S1). For both of these reasons, from now on, we only focus on the deformation rate profile recorded at BH2, which we divide into three parts.

The upper part (uppermost 120 m) of the deformation rate profile at BH2 shows a low shear deformation rate ($\approx 0.02 \text{ a}^{-1}$), except at BH2#19, which exhibits very noisy records and is thus disregarded from our analysis (see Sect. S1). The middle part, from -120 m to -219 m (BH2#6), is characterized by much higher deformation rates, which increase nonlinearly toward the bed, ranging from less than $0.02 \text{ a}^{-1}$ at -120 m to a local maximum of $0.63 \text{ a}^{-1}$ at -219 m. Below -219 m begins the section we refer to as the boundary layer, which includes a 40% decrease in deformation rates over the first 10 m below -219 m, followed by more than a 3-fold increase in du/dz from -230 to -235 m, where the deformation rate reaches a maximum of $1.5 \text{ a}^{-1}$.

## 4.2 Ice flow model

Ice flow as predicted by our model is mainly dominated by along-flow shear ($\tau_{xz}$) and lateral shear ($\tau_{xy}$) (Fig. 3). Lateral shear is significant on the sides of the glacier, where it greatly affects the effective stress ($\tau_e$) and, consequently, the deformation rate (Fig. 3). However, along the centerline near BH2, the magnitude of $\tau_{xy}$ is very small ($< 0.02 \text{ MPa}$) compared to that of along-flow shear ($\tau_{xz}$) ($> 0.08 \text{ MPa}$, where most deformation occurs), meaning our measurements are mainly influenced by $\tau_{xz}$. We find that other components of the stress tensor are small compared to $\tau_{xz}$ ($< 10 \%$), except near the surface, where they can be of a similar amplitude (see Sect. S5). The modeled flow gradient tensor components (du/dx and dw/dz) are very small (see Sect. S5), which suggests that the use of Eq. (1) to calculate du/dz is appropriate (see Sect. 3.2). We also find that the stress tensor is rather insensitive to the choice of rheological parameters and is mainly controlled by glacier geometry (Sect. S5). Comparing our results with a simplified plane-strain model, commonly referred to as the shallow-ice approximation (SIA), we show that shear stress estimated using the SIA is overestimated (see Fig. 4b), leading to an incorrect quantification of the rheological parameters. Following Nye (1965), the use of the shape factor $f = 0.646$ in the SIA formulation accounts for the effect of the parabolic valley shape (considering a half-width-to-thickness ratio of 2, a reasonable approximation of the Argentière Glacier cross section at the study site; see Sect. S6). Applying this correction improves the stress representation obtained using the SIA, but large discrepancies remain (see the dotted green line in Fig. 4b), with the corrected SIA underestimating the effective stress along most of the profile, except at the bottom, where the basal shear stress matches that of the full-Stokes model.

We find that simulations with a constant creep factor (A) and a varying Glen flow law exponent ($n = 3, 4$, or $5$) yield deformation rate profiles exhibiting much lower nonlinearity with depth than that observed (see Fig. 4a). The use of a constant creep factor thus cannot explain the observed profile, regardless of the value of $n$, which turns out to be poorly constrained in this context due to stresses not varying strongly with depth (Fig. 4b). To remain consistent with the commonly used value of $n$ and infer a relevant creep factor for ice flow modeling in general, we thus assume $n = 3$ and quantify A by following the method described in Sect. 3.2 (Fig. 4c).

Using a linear piecewise approximation of A (black line in Fig. 4c), the simulation provides a good match with observations (Fig. 4a), except for the deformation rates recovered in the boundary layer, which are poorly reproduced by the numerical model. This discrepancy is likely due to the bed topography being too low in resolution to capture stress variations near unresolved local bed irregularities, as well as to our observed du/dz values potentially being strongly biased as a result of the neglect of compressive/extensive strains in the present analysis, which are likely significant within the boundary layer (see Sect. 3.2). For this reason, we consider the values of A inferred from the boundary layer to not be relevant. In Fig. 4c, we see that the inferred creep factor, A(z), in the upper half of the glacier is compatible with the value proposed by Cuffey and Paterson (2010) for temperate ice but increases by up to a factor of $\sim 4$ from -140 m down to the top of the boundary layer (-219 m). We discuss this depth-increasing creep factor in Sect. 5.2.

## 4.3 Temporal evolution of velocities

On a seasonal timescale, the surface velocity exhibits an annual cycle, with a gradual increase between December and May, a period of stagnation until mid-September, and a decrease until December (Fig. 5). The period of stagnation coincides with the melting period, as indicated by the observed high discharges. Surprisingly, however, the deformation velocity also shows seasonal variation, with a phase and amplitude similar to those of the surface velocity. As a result, basal velocity remains roughly constant throughout the year. This is in contrast to in situ observations of basal velocity at the cavitometer, where sliding velocities show strong seasonal variability, marked by the highest sliding speeds in July and the lowest in February and March (Fig. 5b; see also Gimbert et al., 2021a). These results are further discussed in Sect. 5.3.

On shorter timescales, the relationships between surface, deformation, and sliding velocities are different from those observed at the seasonal timescale. Peaks in surface velocity that occur at the beginning of the melt season (Togaibekov et al., 2024) coincide with a decrease in deformation velocity and a strong increase in sliding velocity (highlighted as vertical dashed lines in Fig. 5). These peaks are also visible at the cavitometer, particularly during the May events. At the end of the melt season (October), peaks in surface velocities are
<figure><img src="image_4.png" /><figcaption>Figure 3. (a) Modeled shear stress ($\tau_{xz}$), (b) effective stress ($\tau_e$), (c) deformation rate (du/dz), and (d) shear stress ($\tau_{xy}$) along a transversal cross section at the BH2 location. The inclinometers from BH2 are indicated by red dots.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 4. (a) Observed and modeled deformation profiles for BH2, using a uniform creep factor with $n = 3$, $n = 4$, and $n = 5$, as well as a depth-variable creep factor ($A = A(z)$) with $n = 3$. (b) Vertical profiles of different modeled stress components (negative values correspond to compressive stress). The shallow-ice approximation solution, without (SIA) and with (fSIA) a shape factor, is provided for comparison (dotted and solid green lines, respectively). (c) An inferred creep factor calculated using the measured deformation rate and modeled effective stress at BH2, based on Eq. (6) (bottom horizontal axis). The corresponding water content, according to Duval (1977), is shown along the top horizontal axis. The dashed blue line marks the value of $A$ proposed by Cuffey and Paterson (2010) (C. & P. (2010) in the figure). Note that “const” stands for constant.</figcaption></figure>

also mainly explained by changes in sliding speed but do not necessarily coincide with a decrease in deformation velocity, as observed for the early-melt-season peaks.

To better understand the origin of temporal changes in deformation velocities, we investigate the extent to which these changes occur in proportion to changes in du/dz with depth. For this, we define $s_n$ as the standard deviation of the monthly
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Computed time series of daily subglacial discharge (a) and velocities (b) at the Argentière Glacier. Panel (b) shows the surface velocity computed with the linear model for the GNSS station ARG1 ($u_s$), the deformation velocity at BH2 ($u_d$), the inferred basal velocity at BH2 ($u_b$), and the sliding velocity at the cavitometer ($u_{cav}$). The horizontal dashed lines correspond to the average velocities over the studied period. In the shaded areas, an estimate of the sliding velocity (brown line) is shown, assuming a linear trend in the deformation velocity (blue line). The vertical dashed lines indicate speed-up events, which are referred to in the text.</figcaption></figure>

averaged values of $du/dz$, normalized by $du/dz$ values averaged between 15 February and 15 October 2020. We show that $s_n$ is consistently on the order of 20%, with no preferential localization in any specific section of the ice column (Fig. 6b). A more detailed analysis of $du/dz$ variations at each sensor shows that most of these variations exhibit a similar seasonal cycle, with maximum deformation occurring in summer (see Sect. S7). Temporal variations in $du/dz$ at any given depth thus generally occur in proportion to the long-term average value of $du/dz$. Only four sensors (BH2#6, BH2#11, BH2#12, and BH2#14) significantly depart from this general picture, either by exhibiting a different phasing (BH2#6, BH2#12, and BH2#14) or by showing a different amplitude (BH2#11) (Fig. S17), which may be due to poor coupling with the ice (the noise levels for BH2#12 and BH2#14 are particularly high).

# 5 Discussion

## 5.1 Identification and interpretation of the boundary layer

The shape of the deformation rate profile and the retrieved values of $du/dz$ close to the bed suggest intense deformation due to local stress concentrations around bedrock irregularities within the boundary layer (Maier et al., 2019; Weertman, 1957; Gudmundsson, 1997a). We qualitatively explored the viability of this explanation with a simulation of tilt evolution close to the bed, using the deformation rates provided by Gudmundsson (1997b) and the model of tilt evolution for a given velocity field provided by Gudmundsson et al. (1999), as explained in Sect. S8. In Fig. S18, we show that significant rates of compressive or extensive horizontal strain ($du/dx$) localize close to the bed and dominate the flow gradient in a boundary layer with a thickness equal to ~6 times the vertical amplitude of the bedrock bump. We find that under the hypothesis of Eq. (1), and using synthetic tilt time series produced from the estimated strain rates, the inferred apparent $du/dz$ values have a zigzag shape similar to that observed in Fig. 2. This suggests that the zigzag shape in the recon-
<figure><img src="image_7.png" /><figcaption>Figure 6. (a) Observed monthly and yearly du/dz profiles for four periods (see legend). (b) Standard deviation of the monthly averaged du/dz profile, normalized by the du/dz value averaged between 15 February and 15 October 2020. (c) Observed monthly and yearly horizontal-velocity ($u_h$) profiles for four periods (see legend).</figcaption></figure>

structed du/dz values is an artifact resulting from neglecting du/dx rather than a real variation in du/dz. A comparison between our reconstructed du/dz values from the tilt data and the inferred apparent du/dz values using the models of Gudmundsson (1997b) and Gudmundsson et al. (1999) suggests that we drilled on the summit of a bedrock bump with a vertical amplitude of ~2 m.

## 5.2 Interpretation of the depth-increasing creep factor

The Argentière Glacier is fully temperate – i.e., ice lies at the pressure melting point throughout the entire ice column, meaning that the observed increase in the creep factor with depth must result from factors other than changing temperature. We test the hypothesis that the depth increase in the creep factor is due to an associated increase in interstitial water content (W) by linking these two quantities using the formula proposed by Duval (1977). Adapting the formulation from Duval (1977) by considering $A = 50 \text{ MPa}^{-3} \text{ a}^{-1}$ as a reference value when $W = 0$ (i.e., assuming no water content in the upper half of the glacier), we obtain

$$W = \frac{1}{2.34} \left( \frac{A}{50} - 1 \right), \qquad (7)$$

where W is expressed as a percentage and A in terms of MPa$^{-3}$ a$^{-1}$. The inferred water content values are shown along the top horizontal axis of Fig. 4c. Discarding the negative values as artifacts of our chosen parametrization of the Duval (1977) model, we find that the expected water content above –219 m ranges between 0% and 1.3%, increasing toward the bed. Below this depth, the deformation profile is likely influenced by enhanced stress and underestimated velocity gradients due to local basal roughness and cannot be interpreted in terms of an enhanced creep factor (see Sect. 5.1). These values of W and this type of spatial distribution are comparable to those observed for temperate ice (see Fig. 7). The depth-increasing creep factor is thus compatible with the effect of depth-increasing water content.

Our suggestion of increased water content with depth might seem to contradict the findings of Lliboutry and Duval (1985), who reported no relationship between W and depth based on their analysis of an ice core obtained from the Argentière Glacier at a location close to our boreholes (Hantz and Lliboutry, 1983). However, the researchers who performed the measurements considered them unreliable as the water content was primarily correlated with the air temperature.
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Synthesis of water content measurements taken in temperate glaciers, with a distinction made between the top layer (close to the surface) and the bottom layer (close to the bedrock). The colors indicate the methods (ground-penetrating radar (GPR)/seismic methods, calorimeter measurements, and ion concentrations). Note that “pers. comm.” stands for personal communication.</figcaption></figure>

perature at which the measurements were taken at the glacier surface, with higher temperatures increasing ice core melt and thus artificially increasing the measured water content (Michel Vallon, personal communication). For these reasons, we do not include the results from Lliboutry and Duval (1985) in Fig. 7, although they have previously appeared in other compilations of observations of interstitial water content (e.g., Pettersson et al., 2004; Cuffey and Paterson, 2010). We also note that the absolute values of water content, as inferred using the relationship proposed by Duval (1977), may be associated with high uncertainty since this empirical relationship was derived in the laboratory by shearing temperature ice in tertiary creep with a water content of up to 0.8% and has not been validated for higher water content values. Recently, Adams et al. (2021) found in similar experiments that ice viscosity under secondary creep is not sensitive to water content when $W > 0.6\%$. However, we do not expect these results to apply to the present case since tertiary creep is expected under the significant cumulative deformations of the Argentière Glacier (Lliboutry and Duval, 1985; Budd and Jacka, 1989).

Although water content is our best identified candidate for explaining the increase in the creep factor with depth, we note that other factors might also be at play, such as depth-decreasing grain sizes due to the increased rate of mechanical work (Behn et al., 2021) or depth-increasing ice anisotropy due to increased shearing (Cuffey and Paterson, 2010; Montagnat and Duval, 2004). The role of these factors would need

to be specifically analyzed by quantifying the evolution of grain size and ice texture with depth at the measurement site in order to determine whether their respective influences on ice viscosity can be discarded.

## 5.3 Temporal changes in deformation velocity

We show that deformation velocity varies seasonally by ~30% (Fig. 5). Using Glen's flow law with $n = 3$, this implies either a stress change of ~9% or a change in the creep factor of ~30%, of which the latter would require the water content to be multiplied by 1.6 throughout the entire ice column. Although microcracking in response to changes in surface melting and subglacial water pressure can provide significant amounts of water at depth (Gajek et al., 2021), the extent to which this mechanism could alter interstitial water content at the source of the ice softening remains unknown. In addition, an increase in microcracking at depth during summer is inconsistent with the observed decrease in basal water pressure. Aside from microcracking, a 30% increase in deformation heat (≈ +3 × 10⁻⁴ W m⁻³) would produce an excess water content of only ~0.17% over 1 month (the time frame within which the deformation change occurs), which would have a much smaller effect on the creep factor than what was observed. Finally, with regard to the stress variation hypothesis, a 9% stress change cannot be explained by the seasonal evolution of ice thickness, which is at its maximum in May and decreases until October (Vincent et al., 2022), making
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Modeled changes in effective stress (a) and horizontal velocity (b) resulting from nonuniform variations in basal friction. We assume that, during the melting period, friction increases in the central part (thick blue line in panels a and b) and decreases elsewhere due to channelized efficient drainage in the central line contrasting with the inefficient cavity network elsewhere. The inclinometers for BH2 are indicated by black dots in panels (a) and (b). Panel (c) compares the observed velocity profiles at BH2 for August and October with the modeled velocity profiles with nonuniform (channel) and uniform (reference) basal friction.</figcaption></figure>

it out of phase with the observed deformation rate. Furthermore, the amplitude of the thickness change at the BH2 site is about 6 m (Vincent et al., 2022), which represents only a 2 % stress change.

The hypothesis that we find most plausible for explaining the seasonal variability in deformation rates observed at BH2 is that it arises from changes in stress distribution due to nonuniform spatial changes in basal drag during the melting season, as previously suggested in Hooke et al. (1992) and Willis et al. (2003). The loss of drag in response to melt-water input in areas surrounding the observation site may locally increase stress through stress transfer in the form of lateral shear. This may be caused by differences in subglacial hydrological conditions between the glacier's centerline and its sides. The lack of change in basal velocity during summer at BH2 (Fig. 5), despite increases in deformation velocity and stress, suggests that basal friction also increased in the central part of the glacier to accommodate more stress while remaining at a constant velocity. This may have resulted from the development of efficient drainage in the central part of the glacier, as previously identified from seismic observations (Nanni et al., 2021). This efficient drainage reduces water pressure and thus promotes high friction along the central line, while the sides remain dominated by higher-

pressure, inefficient drainage, promoting lower basal friction. To investigate this hypothesis, we perform a complementary numerical simulation in which friction is increased around a central ~50 m wide area (assumed to be affected by efficient drainage) and decreased elsewhere (see Fig. 8a and b). We find that, relative to the reference state from October 2020, a 20 % decrease in $A_s$ in the central part of the glacier, combined with an increase of 70 % in $A_s$ elsewhere, is sufficient for producing the stress change needed to explain the enhanced internal deformation observed in summer (Fig. 8c). Such a contrast in $A_s$ (i.e., a factor of 2.1) is compatible with that expected when cavitation does or does not occur (Gimbert et al., 2021a; Gilbert et al., 2022; Maier et al., 2022). This view is also consistent with the findings of Vincent et al. (2022), which show that bed separation by means of cavitation increases between January and July, with greater amplitudes along the glacier margin than along the midline. This leads to an increase in basal sliding during this period, as observed in our data (Figs. 5 and 6c). An evolving drag contrast between the center and sides of the glacier would also explain why strong peaks in surface velocity are associated with a decrease in deformation velocity during the early melt season. Early water input into the not-yet-developed efficient drainage system along the centerline of the glacier may lead
to pressurization of the central channel, reducing drag along the centerline relative to the sides of the glacier. This effect may be less pronounced in fall because the channels are still of significant size after summer, and this would explain why the decrease in deformation velocity during surface velocity peaks is observed not in fall but only in spring (Fig. 5).

## 6 Conclusions

Using borehole inclinometry, we were able to reconstruct the deformation profile along the central line of the ablation area of the Argentière Glacier and its evolution over 8 months, including the entire melt season. We quantified the Glen flow law creep factor by combining our observations with modeling of the local stress field using a three-dimensional full-Stokes ice flow model, which appears to be different from what the SIA would predict, making its use inappropriate at our study site. We show that surface values (above 100 m depth) of the creep factor are consistent with the standard value for temperate ice (Cuffey and Paterson, 2010) but increase progressively from below this depth up to $200 \text{ MPa}^{-3} \text{ a}^{-1}$. We interpret the depth-increasing creep factor as being associated with an increase in water content from 0 % to 1.3 %, which are reasonable values given previous in situ water content measurements found in the literature (Fig. 7). We also find seasonally evolving deformation rates, with higher deformation rates occurring during the melt season, meaning that the local deformation velocity explains most of the evolution of the surface velocity on a seasonal timescale. We show that this variability is due to changes in stress distribution within the ice body in response to the evolving contrast in basal drag between the centerline and the rest of the glacier during the melt season. We interpret the difference in drag at the center of the glacier as the effect of efficient drainage developing at the deepest point of the subglacial valley.

This study demonstrates the significant value of using borehole inclinometry to infer ice rheology and local changes in basal friction that cannot be detected through surface velocity observations. The data obtained provide rare insights into how subglacial hydrology, basal friction, and surface velocity are interconnected.

*Code and data availability.* All data used in this article, as well as the code used to process the data, can be accessed via the Zenodo repository available at https://doi.org/10.5281/zenodo.13961256 (Roldán-Blasco et al., 2024).

*Supplement.* The supplement related to this article is available online at: https://doi.org/10.5194/tc-19-267-2025-supplement.

*Author contributions.* JPRB processed the tilt data, and AG designed and performed the numerical simulations. JPRB, AG, and FG wrote the paper. LP designed the tiltmeters, directed the field campaign, and performed an early analysis of the data. LP, AG, FG, CV, and OG participated in the inclinometry field campaign and, together with JPRB, analyzed the data, while AT and AW provided the GNSS measurements and velocity time series. All authors contributed to finalizing the paper.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Acknowledgements.* The authors would like to thank Bruno Jourdain, Olivier Laarmann, Maël Richard, and everybody else who participated in the 2019 field campaign and assisted with the fabrication of the instruments. This work has been supported by the French ANR project SAUSSURE (grant no. ANR-18-CE01-0015-01; https://saussure.osug.fr, last access: 15 January 2025). We also thank Luc Moreau for providing the cavitometer data and Electricité d'Emosson SA for the discharge measurements. Glacier surface elevation and precipitation data were acquired within the framework of the GLACIOCLIM program (https://glacioclim.osug.fr, last access: 15 January 2025). The authors would like to thank Dominik Gräff and Manuela Köpfli for their careful reading and suggestions during the review process, which significantly improved the quality of the paper.

*Financial support.* This research has been supported by the Agence nationale de la recherche (grant no. ANR-18-CE01-0015-01). Olivier Gagliardini received support from the project FricFrac, funded by the Centre for Advanced Study (CAS) at the Norwegian Academy of Science and Letters, provided during the academic year 2023–2024.

*Review statement.* This paper was edited by Adam Booth and reviewed by Dominik Gräff and Manuela Köpfli.

## References

Adams, C. J., Iverson, N. R., Helanow, C., Zoet, L. K., and Bate, C. E.: Softening of Temperate Ice by Interstitial Water, Front. Earth Sci., 9, 1–11, https://doi.org/10.3389/feart.2021.702761, 2021.

Amundson, J. M., Truffer, M., and Lüthi, M. P.: Time-dependent basal stress conditions beneath Black Rapids
Glacier, Alaska, USA, inferred from measurements of ice deformation and surface motion, J. Glaciol., 52, 347–357, https://doi.org/10.3189/172756506781828593, 2006.

Arthern, R. J. and Gudmundsson, G. H.: Initialization of ice-sheet forecasts viewed as an inverse Robin problem, J. Glaciol., 56, 527–533, https://doi.org/10.3189/002214310792447699, 2010.

Barnes, P., Tabor, D., and Walker, J. C. F.: The friction and creep of polycrystalline ice, P. Roy. Soc. Lond. A, 324, 127–155, https://doi.org/10.1098/rspa.1971.0132, 1971.

Behn, M. D., Goldsby, D. L., and Hirth, G.: The role of grain size evolution in the rheology of ice: implications for reconciling laboratory creep data and the Glen flow law, The Cryosphere, 15, 4589–4605, https://doi.org/10.5194/tc-15-4589-2021, 2021.

Benjumea, B., Macheret, Y. Y., Navarro, F. J., and Teixidó, T.: Estimation of water content in a temperate glacier from radar and seismic sounding data, Ann. Glaciol., 37, 317–324, https://doi.org/10.3189/172756403781815924, 2003.

Beraud, L., Cusicanqui, D., Rabatel, A., Brun, F., Vincent, C., and Six, D.: Glacier-wide seasonal and annual geodetic mass balances from Pléiades stereo images: application to the Glacier d'Argentière, French Alps, J. Glaciol., 69, 525–537, https://doi.org/10.1017/jog.2022.79, 2022.

Bock, Y., Gourevitch, S. A., Counselman, III, C. C., King, R. W., and Abbot, R. I.: Interferometric analysis of GPS phase observations, Manuscripta Geodaetica, 11, 282–288, 1986.

Booth, A. D., Christoffersen, P., Schoonman, C., Clarke, A., Hubbard, B., Law, R., Doyle, S. H., Chudley, T. R., and Chalari, A.: Distributed Acoustic Sensing of Seismic Properties in a Borehole Drilled on a Fast-Flowing Greenlandic Outlet Glacier, Geophys. Res. Lett., 47, e2020GL088148, https://doi.org/10.1029/2020GL088148, 2020.

Budd, W. F. and Jacka, T. H.: A review of ice rheology for ice sheet modelling, Cold Reg. Sci. Technol., 16, 107–144, https://doi.org/10.1016/0165-232X(89)90014-1, 1989.

Chandler, D., Hubbard, B., Hubbard, A., Murray, T., and Rippin, D.: Optimising ice flow law parameters using borehole deformation measurements and numerical modelling, Geophys. Res. Lett., 35, L12502, https://doi.org/10.1029/2008GL033801, 2008.

Chauve, T., Montagnat, M., Dansereau, V., Saramito, P., Fourteau, K., and Tommasi, A.: A physically-based formulation for texture evolution during dynamic recrystallization. A case study of ice, Comptes Rendus. Mécanique, 352, 99–134, https://doi.org/10.5802/crmeca.243, 2024.

Cohen, D.: Rheology of ice at the bed of engabreen, Norway, J. Glaciol., 46, 611–621, https://doi.org/10.3189/172756500781832620, 2000.

Cuffey, K. M. and Paterson, W. S. B.: The physics of glaciers, Academic Press, Amsterdam, 4th edn., ISBN 9781493300761, 2010.

Doyle, S. H., Hubbard, B., Christoffersen, P., Young, T. J., Hofstede, C., Bougamont, M., Box, J. E., and Hubbard, A.: Physical Conditions of Fast Glacier Flow: 1. Measurements From Boreholes Drilled to the Bed of Store Glacier, West Greenland, J. Geophys. Res.-Earth Surf., 123, 324–348, https://doi.org/10.1002/2017JF004529, 2018.

Duval, P.: The role of the water content on the creep rate of polycrystalline ice, IAHS Publ, 118, 29–33, 1977.

Endres, A. L., Murray, T., Booth, A. D., and West, L. J.: A new framework for estimating englacial water content and pore geometry using combined radar and seismic wave velocities, Geophys.

Res. Lett., 36, L04501, https://doi.org/10.1029/2008GL036876, 2009.

Fichtner, A., Hofstede, C., Gebraad, L., Zunino, A., Zigone, D., and Eisen, O.: Borehole fibre-optic seismology inside the Northeast Greenland Ice Stream, Geophys. J. Int., 235, 2430–2441, https://doi.org/10.1093/gji/gga344, 2023.

Fürst, J. J., Durand, G., Gillet-Chaulet, F., Merino, N., Tavard, L., Mouginot, J., Gourmelen, N., and Gagliardini, O.: Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points, The Cryosphere, 9, 1427–1443, https://doi.org/10.5194/tc-9-1427-2015, 2015.

Gagliardini, O., Zwinger, T., Gillet-Chaulet, F., Durand, G., Favier, L., de Fleurian, B., Greve, R., Malinen, M., Martín, C., Råback, P., Ruokolainen, J., Sacchettini, M., Schäfer, M., Seddik, H., and Thies, J.: Capabilities and performance of Elmer/Ice, a new-generation ice sheet model, Geosci. Model Dev., 6, 1299–1318, https://doi.org/10.5194/gmd-6-1299-2013, 2013.

Gajek, W., Gräff, D., Hellmann, S., Rempel, A. W., and Walter, F.: Diurnal expansion and contraction of englacial fracture networks revealed by seismic shear wave splitting, Commun. Earth Environment, 2, 1–8, https://doi.org/10.1038/s43247-021-00279-4, 2021.

Gilbert, A., Gimbert, F., Thøgersen, K., Schuler, T. V., and Kääb, A.: A Consistent Framework for Coupling Basal Friction With Subglacial Hydrology on Hard-Bedded Glaciers, Geophys. Res. Lett., 49, e2021GL097507, https://doi.org/10.1029/2021GL097507, 2022.

Gilbert, A., Gimbert, F., Gagliardini, O., and Vincent, C.: Inferring the Basal Friction Law From Long Term Changes of Glacier Length, Thickness and Velocity on an Alpine Glacier, Geophys. Res. Lett., 50, e2023GL104503, https://doi.org/10.1029/2023GL104503, 2023.

Gimbert, F., Gilbert, A., Gagliardini, O., Vincent, C., and Moreau, L.: Do Existing Theories Explain Seasonal to Multi-Decadal Changes in Glacier Basal Sliding Speed?, Geophys. Res. Lett., 48, 1–10, https://doi.org/10.1029/2021GL092858, 2021a.

Gimbert, F., Nanni, U., Roux, P., Helmstetter, A., Garambois, S., Lecointre, A., Walpersdorf, A., Jourdain, B., Langlais, M., Laarman, O., Lindner, F., Sergeant, A., Vincent, C., and Walter, F.: A multi-physics experiment with a temporary dense seismic array on the argentière Glacier, French Alps: The RESOLVE project, Seismol. Res. Lett., 92, 1185–1201, https://doi.org/10.1785/0220200280, 2021b.

Glen, J. W.: The creep of polycrystalline ice, P. Roy. Soc. Lond. A, 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955.

Goldsby, D. L. and Kohlstedt, D. L.: Superplastic deformation of ice: Experimental observations, J. Geophys. Res.-Sol. Ea., 106, 11017–11030, https://doi.org/10.1029/2000JB900336, 2001.

Gudmundsson, G. H.: Basal-flow characteristics of a non-linear flow sliding frictionless over strongly undulating bedrock, J. Glaciol., 43, 80–89, https://doi.org/10.1017/s0022143000002835, 1997a.

Gudmundsson, G. H.: Basal-flow characteristics of a linear medium sliding frictionless over small bedrock undulations, J. Glaciol., 43, 71–79, https://doi.org/10.1017/s0022143000002823, 1997b.

Gudmundsson, G. H., Bauder, A., Lüthi, M., Fischer, U. H., and Funk, M.: Estimating rates of basal motion and internal ice deformation from continuous tilt measurements, Ann. Glaciol., 28, 247–252, https://doi.org/10.3189/172756499781821751, 1999.
Hantz, D. and Lliboutry, L.: Waterways, Ice Permeability at Depth, and Water Pressures at Glacier D'Argentière, French Alps, J. Glaciol., 29, 227–239, https://doi.org/10.3189/S0022143000008285, 1983.

Harper, J. T., Humphrey, N. F., Pfeffer, W. T., Huzurbazar, S. V., Bahr, D. B., and Welch, B. C.: Spatial variability in the flow of a valley glacier: Deformation of a large array of boreholes, J. Geophys. Res.-Sol. Ea., 106, 8547–8562, https://doi.org/10.1029/2000jb900440, 2001.

Herring, T. A., King, R. W., and McClusky, S. C.: GPS analysis at MIT, GAMIT reference manual, Massachusetts Institute of Technology, 2018.

Hooke, R. L.: Structure and Flow in the Margin of the Barnes Ice Cap, Baffin Island, N.W.T., Canada, J. Glaciol., 12, 423–438, https://doi.org/10.3189/s0022143000031841, 1973.

Hooke, R. L., Pohjola, V. A., Jansson, P., and Kohler, J.: Intra-seasonal changes in deformation profiles revealed by borehole studies, Storglaciaren, Sweden, J. Glaciol., 38, 348–358, https://doi.org/10.1017/S0022143000002239, 1992.

Hooke, R. L. B. and Hanson, B.: Borehole deformation experiments, Barnes Ice Cap, Canada, Cold Reg. Sci. Technol., 12, 261–276, https://doi.org/10.1016/0165-232X(86)90039-X, 1986.

Hubbard, B. P., Hubbard, A., Mader, H. M., Tison, J. L., Grust, K., and Nienow, P. W.: Spatial variability in the water content and rheology of temperate glaciers: Glacier de Tsanfleuron, Switzerland, Ann. Glaciol., 37, 1–6, https://doi.org/10.3189/172756403781815474, 2003.

Jones, S. J. and Glen, J. W.: The effect of dissolved impurities on the mechanical properties of ice crystals, Philos. Mag., 19, 13–24, https://doi.org/10.1080/14786436908217758, 1969.

Joubert, J.-L.: Stratigraphie de la glace tempérée à l'aide de la teneur en eau liquide, Comptes rendus académie des sciences, p. 3638, 1963.

Keller, A. and Blatter, H.: Measurement of strain-rate components in a glacier with embedded inclinometers, J. Glaciol., 58, 692–698, https://doi.org/10.3189/2012JoG11J234, 2012.

Law, R., Christoffersen, P., MacKie, E., Cook, S., Haseloff, M., and Gagliardini, O.: Complex motion of Greenland Ice Sheet outlet glaciers with basal temperate ice, Sci. Adv., 9, eabq5180, https://doi.org/10.1126/sciadv.abq5180, 2023.

Lee, I. R., Hawley, R. L., Bernsen, S., Campbell, S. W., Clemens-Sewall, D., Gerbi, C. C., and Hruby, K.: A novel tilt sensor for studying ice deformation: Application to streaming ice on Jarvis Glacier, Alaska, J. Glaciol., 66, 74–82, https://doi.org/10.1017/jog.2019.84, 2019.

Legchenko, A., Vincent, C., Baltassat, J. M., Girard, J. F., Thibert, E., Gagliardini, O., Desclaitres, M., Gilbert, A., Garambois, S., Chevalier, A., and Guyard, H.: Monitoring water accumulation in a glacier using magnetic resonance imaging, The Cryosphere, 8, 155–166, https://doi.org/10.5194/tc-8-155-2014, 2014.

Lliboutry, L.: Une théorie du frottement du glacier sur son lit, Annales de Geophysique, 15, 250, 1959.

Lliboutry, L.: General theory of subglacial cavitation and sliding of temperate glaciers, J. Glaciol., 7, 21–58, 1968.

Lliboutry, L.: Permeability, Brine Content and Temperature of Temperate Ice, J. Glaciol., 10, 15–29, https://doi.org/10.3189/s002214300001296x, 1971

Lliboutry, L.: Multivariate Statistical Analysis of Glacier Annual Balances, J. Glaciol., 13, 371–392, https://doi.org/10.3189/s0022143000023169, 1974.

Lliboutry, L. and Duval, P.: Various isotropic and anisotropic ices found in glaciers and polar ice caps and their corresponding rheologies, Int. J. Rock Mech. Min., 22, 198, https://doi.org/10.1016/0148-9062(85)90267-0, 1985.

Lüthi, M., Funk, M., Iken, A., Gogineni, S., and Truffer, M.: Mechanisms of fast flow in Jakobshavn Isbræ, West Greenland: Part III. Measurements of ice deformation, temperature and cross-borehole conductivity in boreholes to the bedrock, J. Glaciol., 48, 369–385, https://doi.org/10.3189/172756502781831322, 2002.

Maier, N., Humphrey, N., Harper, J., and Meierbachtol, T.: Sliding dominates slow-flowing margin regions, Greenland Ice Sheet, Sci. Adv., 5, eaaw5406, https://doi.org/10.1126/sciadv.aaw5406, 2019.

Maier, N., Humphrey, N., Meierbachtol, T., and Harper, J.: Deformation motion tracks sliding changes through summer, western Greenland, J. Glaciol., 68, 187–196, https://doi.org/10.1017/jog.2021.87, 2021.

Maier, N., Gimbert, F., and Gillet-Chaulet, F.: Threshold response to melt drives large-scale bed weakening in Greenland, Nature, 607, 714–720, https://doi.org/10.1038/s41586-022-04927-3, 2022.

Marshall, H. P., Harper, J. T., Pfeffer, W. T., and Humphrey, N. F.: Depth-varying constitutive properties observed in an isothermal glacier, Geophys. Res. Lett., 29, 61-1-61-4, https://doi.org/10.1029/2002GL015412, 2002.

Montagnat, M. and Duval, P.: The viscoplastic behaviour of ice in polar ice sheets: experimental results and modelling, Comptes Rendus Physique, 5, 699–708, https://doi.org/10.1016/j.crhy.2004.06.002, 2004.

Mosbeux, C., Gillet-Chaulet, F., and Gagliardini, O.: Comparison of adjoint and nudging methods to initialise ice sheet model basal conditions, Geosci. Model Dev., 9, 2549–2562, https://doi.org/10.5194/gmd-9-2549-2016, 2016.

Murray, T., Stuart, G. W., Fry, M., Gamble, N. H., and Crabtree, M. D.: Englacial water distribution in a temperate glacier from surface and borehole radar velocity analysis, J. Glaciol., 46, 389–398, https://doi.org/10.3189/172756500781833188, 2000

Murray, T., Booth, A., and Rippin, D. M.: Water-content of Glaciers ice: Limitations on estimates from velocity analysis of surface ground-penetrating radar surveys, J. Environ. Eng. Geoph., 12, 87–99, https://doi.org/10.2113/JEEG12.1.87, 2007.

Nanni, U., Gimbert, F., Roux, P., and Lecointre, A.: Observing the subglacial hydrology network and its dynamics with a dense seismic array, P. Natl. Acad. Sci. USA, 118, e2023757118, https://doi.org/10.1073/pnas.2023757118, 2021.

Nye, J. F.: The Flow of a Glacier in a Channel of Rectangular, Elliptic or Parabolic Cross-Section, J. Glaciol., 5, 661–690, https://doi.org/10.3189/s0022143000018670, 1965.

Ogier, C., Manen, D.-J. V., Maurer, H., Räss, L., Hertrich, M., Bauder, A., and Farinotti, D.: Ground penetrating radar in temperate ice: englacial water inclusions as limiting factor for data interpretation, J. Glaciol., 69, 1874–1885, https://doi.org/10.1017/jog.2023.68, 2023.

Perutz, M. F.: Direct Measurement of the Velocity Distribution in a Vertical Profile Through a Glacier, J. Glaciol., 1, 382–383, https://doi.org/10.3189/s0022143000012594, 1949.
Pettersson, R., Jansson, P., and Blatter, H.: Spatial variability in water content at the cold-temperate transition surface of the polythermal Storglaciären, Sweden, J. Geophys. Res.-Earth Surf., 109, F02009, https://doi.org/10.1029/2003jf000110, 2004.

Rabatel, A., Sanchez, O., Vincent, C., and Six, D.: Estimation of Glacier Thickness From Surface Mass Balance and Ice Flow Velocities: A Case Study on Argentière Glacier, France, Front. Earth Sci., 6, 112, https://doi.org/10.3389/feart.2018.00112, 2018.

Rathmann, N. M. and Lilien, D. A.: Inferred basal friction and mass flux affected by crystal-orientation fabrics, J. Glaciol., 68, 236–252, https://doi.org/10.1017/jog.2021.88, 2022.

Raymond, C.: Flow in a Transverse Section of Athabasca Glacier, Alberta, Canada, J. Glaciol., 10, 55–84, https://doi.org/10.3189/s0022143000012995, 1971.

Roldan-Blasco, J. P., Gilbert, A., Piard, L., Gimbert, F., Christian, V., Gagliardini, O., Togaibekov, A., Walpersdorf, A., and Maier, N.: Data for “Creep enhancement and sliding in a temperate, hard-bedded alpine glacier”, Zenodo [code and data set], https://doi.org/10.5281/zenodo.13961256, 2024.

Röthlisberger, H.: Water Pressure in Intra- and Subglacial Channels, J. Glaciol., 11, 177–203, https://doi.org/10.3189/S0022143000022188, 1972.

Ryser, C., Lüthi, M. P., Andrews, L. C., Hoffman, M. J., Cataisia, G. A., Hawley, R. L., Neumann, T. A., and Kristensen, S. S.: Sustained high basal motion of the Greenland ice sheet revealed by borehole deformation, J. Glaciol., 60, 647–660, https://doi.org/10.3189/2014JoG13J196, 2014.

Schoof, C.: Ice-sheet acceleration driven by melt supply variability, Nature, 468, 803–806, https://doi.org/10.1038/nature09618, 2010.

Sergeant, A., Chmiel, M., Lindner, F., Walter, F., Roux, P., Chaput, J., Gimbert, F., and Mordret, A.: On the Green's function emergence from interferometry of seismic wave fields generated in high-melt glaciers: implications for passive imaging and monitoring, The Cryosphere, 14, 1139–1171, https://doi.org/10.5194/tc-14-1139-2020, 2020.

Shreve, R. and Sharp, R.: Internal Deformation and Thermal Anomalies in Lower Blue Glacier, Mount Olympus, Washington, U.S.A., J. Glaciol., 9, 65–86, https://doi.org/10.3189/S0022143000026800, 1970.

Togaibekov, A., Gimbert, F., Gilbert, A., and Walpersdorf, A.: Observing and Modeling Short-Term Changes in Basal Friction During Rain-Induced Speed-Ups on an Alpine Glacier, Geophys. Res. Lett., 51, e2023GL107999, https://doi.org/10.1029/2023GL107999, 2024.

Vallon, M., Petit, J.-R., and Fabre, B.: Study of an Ice Core to the Bedrock in the Accumulation zone of an Alpine Glacier, J. Glaciol., 17, 13–28, https://doi.org/10.3189/S0022143000030677, 1976.

Vincent, C. and Moreau, L.: Sliding velocity fluctuations and subglacial hydrology over the last two decades on Argentière glacier, Mont Blanc area, J. Glaciol., 62, 805–815, https://doi.org/10.1017/jog.2016.35, 2016.

Vincent, C., Fischer, A., Mayer, C., Bauder, A., Galos, S. P., Funk, M., Thibert, E., Six, D., Braun, L., and Huss, M.: Common climatic signal from glaciers in the European Alps over the last 50 years, Geophys. Res. Lett., 44, 1376–1383, https://doi.org/10.1002/2016GL072094, 2017.

Vincent, C., Gilbert, A., Walpersdorf, A., Gimbert, F., Gagliardini, O., Jourdain, B., Roldan Blasco, J. P., Laarman, O., Piard, L., Six, D., Moreau, L., Cusicanqui, D., and Thibert, E.: Evidence of Seasonal Uplift in the Argentière Glacier (Mont Blanc Area, France), J. Geophys. Res.-Earth Surf., 127, e2021JF006454, https://doi.org/10.1029/2021JF006454, 2022.

Vivian, R. and Bocquet, G.: Subglacial Cavitation Phenomena Under the Glacier D'Argentière, Mont Blanc, France, J. Glaciol., 12, 439–451, https://doi.org/10.3189/S0022143000031853, 1973.

Weertman, J.: On the Sliding of Glaciers, J. Glaciol., 3, 33–38, https://doi.org/10.3189/s0022143000024709, 1957.

Weertman, J.: CREEP DEFORMATION OF ICE, Annu. Rev. Earth Planet. Sci., 11, 215–240, https://doi.org/10.1146/annurev.ea.11.050183.001243, 1983.

Willis, I., Mair, D., Hubbard, B., Nienow, P., Fischer, U. H., and Hubbard, A.: Seasonal variations in ice deformation and basal motion across the tongue of Haut Glacier d'Arolla, Switzerland, Ann. Glaciol., 36, 157–167, https://doi.org/10.3189/172756403781816455, 2003.

Young, T. J., Martín, C., Christoffersen, P., Schroeder, D. M., Tullaczyk, S. M., and Dawson, E. J.: Rapid and accurate polarimetric radar measurements of ice crystal fabric orientation at the Western Antarctic Ice Sheet (WAIS) Divide ice core site, The Cryosphere, 15, 4117–4133, https://doi.org/10.5194/tc-15-4117-2021, 2021.

Zryd, A.: Conditions dans la couche basale des glaciers tempérés: contraintes, teneur en eau et frottement intérieur, Ph.D. thesis, ETH, 1991.