Hostname: page-component-7c8c6479df-8mjnm Total loading time: 0 Render date: 2024-03-29T07:10:40.969Z Has data issue: false hasContentIssue false

Premelting increases the rate of regelation by an order of magnitude

Published online by Cambridge University Press:  23 May 2019

ALAN W. REMPEL*
Affiliation:
Department of Earth Sciences, University of Oregon, Eugene, OR, USA
COLIN R. MEYER
Affiliation:
Department of Earth Sciences, University of Oregon, Eugene, OR, USA
*
Correspondence: Alan W. Rempel <rempel@uoregon.edu>
Rights & Permissions [Opens in a new window]

Abstract

Glacier sliding over small obstacles relies on melting on their upstream sides and refreezing downstream. Previous treatments have appealed to ‘pressure melting’ as the cause of the spatial variations in melting temperature that drive this regelation process. However, we show that typical liquid pressure variations across small obstacles are negligible and therefore variations in ice pressure closely approximate variations in effective stress. For a given change in effective stress, the equilibrium melting temperature changes by an order of magnitude more than when the pressure of ice and liquid both change by an equal amount. In consequence, the temperature gradients that drive heat flow across small obstacles are larger than previously recognized and the rate of regelation is faster. Under typical conditions, the transition wavelength at which ice deformation and regelation contribute equally is of m-scale, several times longer than previous predictions, which have been reported to underestimate field inferences.

Type
Letter
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2019

1. INTRODUCTION AND SCALING ARGUMENTS

Textbook treatments of glacier sliding based on pioneering work by Weertman (Reference Weertman1957); Nye (Reference Nye1969); Kamb (Reference Kamb1970) describe how glacier motion is facilitated through a combination of ice deformation over larger bed obstacles, and the process of melting and refreezing across smaller obstacles that is referred to as regelation (Cuffey and Paterson, Reference Cuffey and Paterson2010). In these treatments, the rate of regelation is limited by the conduction of latent heat across the obstacle in response to changes in melting temperature that accompany variations in ice pressure. Here, we follow the classical scaling arguments for the rate of regelation (e.g. Fowler, Reference Fowler2011) to demonstrate that the permeability between the ice and the obstacle, which limits the rate at which water can flow, requires that significant local differences arise between the ice and liquid pressures. When the pressure in a liquid differs from that imposed on its adjacent solid (more precisely, the first invariant of the stress tensor), the theory of premelting describes how equilibrium is achieved at a temperature that differs from the bulk melting temperature that holds when the pressures in the two phases are the same (Dash and others, Reference Dash, Rempel and Wettlaufer2006). Taking premelting into account, we modify the Kamb (Reference Kamb1970) treatment of sliding and demonstrate both that the rate of regelation is faster than previously thought, and that the effects of pressure melting and refreezing are typically of negligible importance in comparison to melting and refreezing sourced to variations in effective stress.

We begin with the scaling for conduction of the latent heat due to the refreezing of water. The rate of regelation v R is linked to the time required to transport volumetric latent heat $\rho_{i} {\cal L}$ by conduction at effective thermal conductivity K from the downstream side of small obstacles, where the melting temperature is elevated slightly, to the upstream side, where higher ice pressures cause the melting temperature to be depressed. For obstacles of wavelength λ, scaling considerations suggest that energy transport requires

(1)$$v_R\sim -{{2\,K \Delta T} \over {\rho_{i} {\cal L} \lambda}}\;,$$

where ΔT is the (negative) difference in melting temperature between the upstream and downstream sides. Conventionally, this melting temperature difference is linked to the (positive) difference in ice pressure ΔP i between the upstream and downstream sides through the Clapeyron slope, i.e. ΔT = −C 0ΔP i, with

(2)$$C_0= {{T_{\rm m}} \over {\rho_{i} {\cal L}}} {{\rho_{l} - \rho_{i}} \over {\rho_{l}}}\;,$$

where T m ≈ 273 K is a reference melting temperature, and ρ i and ρ l are the ice and liquid densities; C 0 is sensitive to the air-content dissolved in the liquid (e.g. Röthlisberger and Lang, Reference Röthlisberger, Lang, Gurnell and Clark1987; Clarke, Reference Clarke2005), with reported values ranging between 9.8 × 10−8 K Pa−1 (saturated) and 7.4 × 10−8 K Pa−1 (pure). Where this description of phase equilibrium is valid, it follows that the regelation velocity scales as

(3)$$v_R\sim {{2\,K C_0 \Delta P_i} \over {\rho_{i} {\cal L} \lambda}}\;.$$

As noted above, the equilibrium between liquid water and ice tracks along the Clapeyron slope of the bulk phase diagram only when the ice and liquid pressures are equal to each other along their common interface so that ΔP i = ΔP l. Considerations of mass conservation constrain the difference in liquid pressure between the upstream and downstream sides that is required to drive transport of liquid with viscosity μ, whereupon

(4)$$\Delta P_l \sim {{\rho_{i} \mu v_R \lambda } \over {\rho_{l}\Pi}}\;.$$

Here, the resistance to flow at low Reynolds number is characterized by effective permeability Π; for example, Darcy flow through a porous medium with the permeability Π or an equivalent Poiseuille flow through a liquid film of thickness δ, which can be described with Π = δ2/12. Substituting in for the regelation velocity scaling from the energy transport argument that led to Eqn (1) suggests that

(5)$$\Delta P_l\sim - {{ 2\mu K \Delta T} \over {\rho_{l} {\cal L} \Pi}}\;.$$

Since the melting temperature difference between the upstream and downstream sides falls along the Clapeyron slope so that ΔT = −C 0ΔP i only when ΔP l = ΔP i, the special case in which regelation is driven by pressure melting alone holds when the permeability to liquid transport is

(6)$$\Pi\sim\Pi_0= {{ 2\mu K C_0} \over {\rho_{l} {\cal L}}}\;.$$

Taking nominal values of μ = 1.8 × 10−3 Pa s, K = 2 W m−1 K−1, ${\cal L}= 3.3\times 10^5\,{\rm J}\,{\rm kg}^{-1}$, and treating the case of pure water so that ρ l = 103 kg m−3 and C 0 = 7.4 × 10−8 K Pa−1, implies that Π0 ≈ 1.6 × 10−18 m2, which in the case of thin film flow corresponds with δ0 ≈ 4.4 nm.

Permeabilities as low as Π0 are uncommon in surficial geological materials. Moreover, when the proximity between ice and mineral substrates δ is reduced below the micron scale, the intermolecular forces that cause premelting (e.g. Dash and others, Reference Dash, Fu and Wettlaufer1995, Reference Dash, Rempel and Wettlaufer2006; Wettlaufer and Worster, Reference Wettlaufer and Worster2006) enable significant load transfers between the ice and mineral surfaces so that the ice pressure exceeds the liquid pressure, with their difference identified as the effective stress N = P i − P l. Hence, for the idealized case of sliding over a completely impermeable bed with all liquid transport confined within an interfacial film, we infer that the film thickness δ0 implied by mass-balance constraints is inconsistent with the central assumption of regelation by pressure melting alone, which holds that P i = P l and N=0. Instead, the generalized Clapeyron equation (see Fig. 1) must be invoked to describe the changes in melting temperature that result from variations in both liquid pressure and effective stress as (cf. Rempel, Reference Rempel2008, Eqn (1))

(7)$$\Delta T \approx -C_0\left(\Delta P_l + {{\rho_{l}} \over {\rho_{l}-\rho_{i}}}\Delta N\right)\;,$$

where each of these quantities (i.e. ΔT, ΔP l, ΔN = ΔP i − ΔP l) represents the difference between values on the upstream and downstream sides of bed obstacles. Substituting Eqn (7) into Eqn (1) gives a regelation velocity of

(8)$$v_R\sim {{\rho_{l}} \over {\rho_{l}-\rho_{i}}} {{2\,K C_0 \Delta P_i} \over {\rho_{i} {\cal L} \lambda}} \left[1- {{\rho_{i}\Delta P_l} \over {\rho_{l}\Delta P_i}} \right]\;.$$

Fig. 1. Generalized Clapeyron conditions for regelation. (a) Schematic of regelation past a bump. (b) Gibbs free energy per molecule of ice and liquid water. The ice lines are more shallow because ice has a lower specific entropy than liquid. Separate lines for a given substance correspond to different pressures. Intersections of ice and liquid lines correspond to phase equilibria. Consider the intersection of the two green lines (e.g. at the pressure conditions on the downstream side of a bump): at temperatures below the phase equilibrium temperature, ice has the lower free energy, and above the melting temperature, liquid prevails due to its lower free energy. The red lines show the free energy in each phase when pressure is elevated as labeled. If the ice and liquid pressures both increase by the same amount, the melting temperature decreases along the Clapeyron slope with the intersection between the dashed and the dot-dashed red lines. If, however, the ice and liquid pressures are perturbed by different amounts (e.g. at the pressure conditions on the upstream side of a bump), the melting temperature decrease can be much larger, as shown by the intersection of the two dot-dashed red lines.

In the typical case where Π ≫ Π0 from Eq. (6) efficient liquid drainage is expected to diminish liquid pressure variations so that ΔP i ~ (Π/Π0P l ≫ ΔP l, and in this limit Eqn (8) implies a regelation velocity of

(9)$$v_R\sim {{\rho_{l}} \over {\rho_{l}-\rho_{i}}} {{2\,K C_0 \Delta P_i} \over {\rho_{i} {\cal L} \lambda}}\;.$$

For pure water, Eqn (9) suggests that v R is ρ l/(ρ l − ρ i) ≈ 12 times faster than predicted by conventional treatments, which rely upon temperature differences that result from pressure melting alone (cf. Eqn (3)). The enhanced regelation speed predicted by Eqn (9) is the primary finding of this work.

2. LINEARIZED ICE FLOW OVER WAVY TOPOGRAPHY

Further insight can be gained by revisiting the detailed analysis provided by Kamb (Reference Kamb1970) for regelation (at velocity v R) and ‘plastic’ deformation (at velocity v p) over a rough glacier bed, while making the minor extensions required to describe how the water pressure P l varies independently from the ice pressure P i. The details of this calculation are provided in the Supplementary Information. The interface we consider is described by

(10)$$z_0(x,y)= {{A} \over {4\pi^2}} \int\!\!\!\!\int_{-\infty}^{\infty} a(h,k){\rm e}^{i(hx+ky)}\,{\rm d}h{\rm d}k\;,$$

for area element A, Fourier spectrum a, and horizontal wavenumbers h (along flow) and k (perpendicular to flow). We account for the transport of heat and ice of viscosity η as in Kamb (Reference Kamb1970), while assuming that liquid flow follows Darcy's law and that the interface is at its melting temperature everywhere, as given by Eqn (7). This results in a modified version of Kamb's expression relating the shear stress τ and sliding velocity v = v R + v P, given by

(11)$$\tau= {{A} \over {4\pi^2}} 2\eta v \int\!\!\!\!\int_{-\infty}^{\infty} h^2 \ell a^2(h,k) {{\ell_0^2+\ell_\Pi^2} \over {\ell^2+\ell_0^2+\ell_\Pi^2}} \,{\rm d}h{\rm d}k\;.$$

where ℓ2 = h 2 + k 2, and the characteristic wavenumbers ℓ0 and $\ell _\Pi $ are defined by

(12)$$\ell_0^2 = {{\rho_{l}-\rho_{i}} \over {\rho_{l}}} {{\rho_{i} {\cal L}} \over {4\eta C_0K}} \;\;{\rm and}\;\;\ell_\Pi^2= {{\rho_{i}^{2} \mu} \over {\rho_{l}^{2} 2 \eta \Pi }}\;.$$

The length scale 2π/ℓ0 is the transition wavelength at which the conduction of latent heat required by regelation and the deformation of ice over an obstacle offer equal resistance, whereas $2\pi /\ell _\Pi $ marks distances for which the liquid transport needed to accommodate regelation imposes a significant impediment. Since reasonable permeability estimates satisfy Π ≫ ρ iΠ0/[(ρ l − ρ i)], a comparison of the limiting wave numbers under such conditions yields $\ell _\Pi ^2\ll \ell _0^2$, implying that liquid transport plays a negligible role in restricting regelation. Hence, as in conventional treatments of regelation due to pressure melting (Weertman, Reference Weertman1957; Nye, Reference Nye1969; Kamb, Reference Kamb1970; Gudmundsson, Reference Gudmundsson1997), the relative importance of plastic deformation to regelation in controlling sliding speed depends on a comparison between the topographic wavelength 2π/ℓ and the transition wavelength 2π/ℓ0. Importantly, the transition wavelength is governed by considerations of heat flow, and because liquid transport is relatively efficient, variations in ice pressure produce commensurate variations in effective stress. The resultant perturbations to the melting temperature (see Fig. 1) cause the transition wavelength to be $\sqrt {(\rho _l-\rho _i)/\rho _l}\approx 3.5$ times larger than would be expected if the phase behavior were governed by pressure melting alone. For example, adopting the nominal parameter values given above and assigning an ice viscosity of η = 10−13 Pa s leads to 2π/ℓ0 ≈ 3 m.

Of particular interest here, for interfaces with roughness only at wavelengths small enough that $\ell ^2\gg \ell _0^2 + \ell _\Pi ^2$ so that plastic deformation is negligible (i.e. v ≈ v R), Eqn (11) reduces to

(13)$$\eqalign{\tau\approx {{A} \over {4\pi^2}} {{\rho_{i} {\cal L}v} \over {2C_0K}} \left({{\rho_{l}-\rho_{i}} \over {\rho_{l}}} + {{\rho_{i}} \over {\rho_{l}}} {{\Pi_0 } \over {\Pi }}\right) \int\!\!\!\!\int_{-\infty}^{\infty} {{h^2a^2(h,k) } \over {\ell}} \,{\rm d}h{\rm d}k,}$$

where Π0 is defined as in Eqn (6). In the special case where Π = Π0 we note that the term in parentheses on the right tends to unity so that the relationship between τ and v predicted by treatments that set P i = P l and consider pressure melting alone (e.g. Weertman, Reference Weertman1957; Nye, Reference Nye1969; Kamb, Reference Kamb1970) is recovered; with Π ≫ Π0, for a given basal shear stress τ the regelation velocity predicted by Eqn (13) is a factor of ρ l/(ρ l − ρ i) ≈ 12 greater.

3. IMPLICATIONS

In attempting to validate aspects of his treatment of glacier sliding Kamb (Reference Kamb1970, see §17) compared his predicted transition wavelength to field observations and found evidence suggesting that conventional theory underestimated 2π/ℓ0. In revising expectations of the transition wavelength upwards by a factor of 3.5, the current treatment may yield an explanation for this discrepancy and is promising for efforts to map basal roughness at the resolution of the Weertman (Reference Weertman1957) ‘controlling obstacles’. Perhaps more important than such quantitative improvements is the enhanced recognition of the importance of the effective stress in basal dynamics, and particularly basal thermodynamics. Though the role of subglacial water pressure in controlling sliding behavior through its effects on subglacial mechanics is well established in a variety of contexts (Lliboutry, Reference Lliboutry1968; Iken and Bindschadler, Reference Iken and Bindschadler1986; Tulaczyk and others, Reference Tulaczyk, Kamb and Engelhardt2000; Schoof, Reference Schoof2010), its importance in affecting the ice–liquid equilibrium temperature and premelting (Dash and others, Reference Dash, Rempel and Wettlaufer2006) has arguably been under-appreciated (Alley and others, Reference Alley1997).

Here, we have demonstrated that liquid transport near the basal interface (primarily as a consequence of premelting rather than the density difference that determines the sign of the Clapeyron slope) causes an order of magnitude increase to the speed of regelation. Indeed, because effective stress variations drive melting and refreezing in the same qualitative manner in any system where a solid is wetted by its own melt (i.e. premelting takes place), this suggests a more general role for regelation as a viable mechanism for accommodating slip over small obstacles in diverse geological systems where pressure melting cannot occur because the melt is less dense than its solid. In the glacial context, we anticipate that effective-stress-driven changes in melting temperature are likely to be common causes of basal melting and refreezing, possibly even supercooling (e.g. Röthlisberger and Lang, Reference Röthlisberger, Lang, Gurnell and Clark1987; Alley and others, Reference Alley, Lawson, Evenson, Strasser and Larson1998; Lawson and others, Reference Lawson1998; Creyts and Clarke, Reference Creyts and Clarke2010). Evidence has been given that the thermo-mechanical feedbacks made possible by premelting and variations in effective stress are instrumental in limiting both the range of observed glacier driving stresses (Meyer and others, Reference Meyer, Downey and Rempel2018) and the occurrence of regions amenable to stick-slip behavior (Lipovsky and others, Reference Lipovsky2019). This adds further emphasis to the importance of basal liquid transport for understanding glacier dynamics (e.g. Clarke, Reference Clarke2005; Iverson, Reference Iverson2010).

ACKNOWLEDGEMENTS

We thank Tim Creyts, Neal Iverson, Brad Lipovsky, Elisa Mantelli, Brent Minchew, Kiya Riverman, Olga Sergienko and Christian Schoof for discussions; we appreciate the guidance of scientific editor Ian Hewitt, and the constructive comments of Martin Lüthi and an anonymous reviewer; we also acknowledge financial support from NSF-1603907.

SUPPLEMENTARY MATERIAL

The supplementary material for this article can be found at https://doi.org/10.1017/jog.2019.33

References

Alley, RB and 5 others (1997) How glaciers entrain and transport basal sediment: physical constraints. Quat. Sci. Rev., 16(9), 10171038Google Scholar
Alley, RB, Lawson, DE, Evenson, EB, Strasser, JC and Larson, GJ (1998) Glaciohydraulic supercooling: a freeze-on mechanism to create stratified, debris-rich basal ice: II. Theory. J. Glaciol., 44(148), 563569Google Scholar
Clarke, GKC (2005) Subglacial processes. Annu. Rev. Earth Planet. Sci., 33, 247276Google Scholar
Creyts, TT and Clarke, GKC (2010) Hydraulics of subglacial supercooling: theory and simulations for clear water flows. J. Geophys. Res., 115, F03021 (doi: 10.1029/2009JF001417)Google Scholar
Cuffey, KM and Paterson, WSB (2010) The physics of glaciers, 4th edn. Elsevier, San FranciscoGoogle Scholar
Dash, JG, Fu, H and Wettlaufer, JS (1995) The premelting of ice and its environmental consequences. Rep. Prog. Phys., 58(1), 115Google Scholar
Dash, JG, Rempel, AW and Wettlaufer, JS (2006) The physics of premelted ice and its geophysical consequences. Rev. Mod. Phys., 78(3), 695Google Scholar
Fowler, A (2011) Mathematical geoscience, vol. 36. Springer Science & Business Media, LondonGoogle Scholar
Gudmundsson, GH (1997) Basal-flow characteristics of a non-linear flow sliding frictionless over strongly undulating bedrock. J. Glaciol., 43(143), 8089Google Scholar
Iken, A and Bindschadler, RA (1986) Combined measurements of subglacial water pressure and surface velocity of Findelengletscher, Switzerland: conclusions about drainage system and sliding mechanism. J. Glaciol., 32(110), 101119Google Scholar
Iverson, NR (2010) Shear resistance and continuity of subglacial till: hydrology rules. J. Glaciol., 56(200), 11041114Google Scholar
Kamb, B (1970) Sliding motion of glaciers: theory and observation. Rev. Geophys., 8(4), 673728Google Scholar
Lawson, DE and 5 others (1998) Glaciohydraulic supercooling: a freeze-on mechanism to create stratified, debris-rich basal ice: I. Field evidence. J. Glaciol., 44(148), 547562Google Scholar
Lipovsky, BP and 5 others (2019) Glacier sliding, seismicity, and sediment entrainment. Ann. Glaciol. (doi: 10.1017/aog.2019.24)Google Scholar
Lliboutry, L (1968) General theory of subglacial cavitation and sliding of temperate glaciers. J. Glaciol., 7, 2158Google Scholar
Meyer, CR, Downey, AS and Rempel, AW (2018) Freeze-on limits bed strength beneath sliding glaciers. Nat. Commun., 9(1), 3242Google Scholar
Nye, JF (1969) A calculation on the sliding of ice over a wavy surface using a Newtonian viscous approximation. Proc. R. Soc. Lond. A, 311(1506), 445467Google Scholar
Rempel, AW (2008) A theory for ice-till interactions and sediment entrainment beneath glaciers. J. Geophys. Res., 113(F1), f01013Google Scholar
Röthlisberger, H and Lang, H (1987) Glacial hydrology. In Gurnell, AM and Clark, MJ, eds. Glacio-fluvial sediment transfer: an alpine perspective. John Wiley and Sons, New York, 207284Google Scholar
Schoof, C (2010) Ice sheet acceleration driven by melt supply variability. Nature, 468(7325), 803806Google Scholar
Tulaczyk, S, Kamb, WB and Engelhardt, HF (2000) Basal mechanics of Ice Stream B, West Antarctica: 1. Till mechanics. J. Geophys. Res., 105(B1), 463481Google Scholar
Weertman, J (1957) On the sliding of glaciers. J. Glaciol, 3(21), 3338Google Scholar
Wettlaufer, JS and Worster, MG (2006) Premelting dynamics. Ann. Rev. Fluid Mech., 38(1), 427452Google Scholar
Figure 0

Fig. 1. Generalized Clapeyron conditions for regelation. (a) Schematic of regelation past a bump. (b) Gibbs free energy per molecule of ice and liquid water. The ice lines are more shallow because ice has a lower specific entropy than liquid. Separate lines for a given substance correspond to different pressures. Intersections of ice and liquid lines correspond to phase equilibria. Consider the intersection of the two green lines (e.g. at the pressure conditions on the downstream side of a bump): at temperatures below the phase equilibrium temperature, ice has the lower free energy, and above the melting temperature, liquid prevails due to its lower free energy. The red lines show the free energy in each phase when pressure is elevated as labeled. If the ice and liquid pressures both increase by the same amount, the melting temperature decreases along the Clapeyron slope with the intersection between the dashed and the dot-dashed red lines. If, however, the ice and liquid pressures are perturbed by different amounts (e.g. at the pressure conditions on the upstream side of a bump), the melting temperature decrease can be much larger, as shown by the intersection of the two dot-dashed red lines.

Supplementary material: PDF

Rempel and Meyer supplementary material

Rempel and Meyer supplementary material 1

Download Rempel and Meyer supplementary material(PDF)
PDF 193 KB