Glen's Flow Law

How can a seemingly solid, brittle substance like ice flow like a river, carving continents and shaping planetary landscapes? This apparent contradiction lies at the heart of glaciology and is explained by a fundamental principle known as Glen's Flow Law. This law moves beyond our everyday intuition, revealing that under the immense pressure of its own weight, ice deforms and flows in a process called creep. This article addresses the physics behind this large-scale motion, bridging the gap between a small ice crystal and a continent-sized ice sheet.

This exploration is divided into two parts. In the first chapter, ​​Principles and Mechanisms​​, we will dissect Glen's Flow Law itself. You will learn why ice is a "non-Newtonian" fluid, how its flow is profoundly sensitive to both stress and temperature, and how these factors create powerful feedback loops that can cause flow to accelerate. Following this, the chapter on ​​Applications and Interdisciplinary Connections​​ will demonstrate the law's remarkable power. We will journey from the complex computer models used to predict sea-level rise to the practical engineering challenges of building on ice, and even to the frozen plains of distant planets, showing how this single rule unifies a vast range of scientific fields.

To understand how a colossal river of ice can carve its way through solid rock, we must first abandon our everyday intuition about solids. We think of a block of ice as rigid, brittle, and unyielding. And over the course of a minute, or an hour, it is. But over years, decades, and centuries, under the immense pressure of its own weight, ice behaves in a most peculiar way: it flows. This slow, steady deformation is a process called ​​creep​​, and it transforms a seemingly static ice sheet into one of the most powerful and dynamic forces on the planet.

But ice does not flow like a simple liquid such as water or honey. For these familiar fluids, the relationship between the force applied (the ​​stress​​) and the rate of flow (the ​​strain rate​​) is straightforwardly linear. Double the push, and you double the speed. Physicists call these ​​Newtonian fluids​​. Ice, however, plays by a different set of rules—a rulebook discovered through painstaking experiments by the glaciologist John W. Glen in the 1950s. This rulebook, now known as ​​Glen's Flow Law​​, is the key to unlocking the secret life of glaciers.

At its heart, Glen's Flow Law is a deceptively simple-looking power law. In its scalar form, it states:

Here, $\dot{\epsilon}_{e}$ represents the ​​effective strain rate​​, a measure of how quickly the ice is deforming. $\tau_{e}$ is the ​​effective stress​​, a measure of the internal forces compelling it to move. And holding the secret to the strange behavior of ice are the two parameters, $n$ and $A$.

The Exponent $n$: The Non-Newtonian Heart of the Matter

The exponent $n$ governs the relationship between stress and strain rate. If ice were a simple Newtonian fluid, $n$ would be exactly 1. The flow rate would be directly proportional to the stress, and the material would have a constant, well-defined viscosity. In a simplified glacier, this would produce a smooth, almost parabolic velocity profile, with the fastest flow at the surface and gradually decreasing to zero at the bed.

But Glen's experiments revealed that for polycrystalline ice, $n$ is not 1; it's approximately 3. This seemingly small change has profound consequences. With $n=3$, the law becomes $\dot{\epsilon}_{e} \propto \tau_{e}^3$. This means that doubling the stress on the ice does not merely double its rate of flow—it increases it by a factor of $2^3$, or eight times! This extreme sensitivity to stress is the signature of a ​​non-Newtonian​​ material.

This non-linearity forces us to rethink the very idea of viscosity. For a Newtonian fluid, viscosity is a fixed property. For ice, it's a dynamic variable. We can define an ​​effective viscosity​​, $\eta_{\text{eff}}$, that relates stress and strain rate through $\tau_{ij} = 2\eta_{\text{eff}}\dot{\epsilon}_{ij}$. By combining this with Glen's Law, we find that the effective viscosity of ice is not constant but depends on the stress itself:

Since $n \approx 3$, this means $\eta_{\text{eff}} \propto \tau_e^{-2}$. This phenomenon is known as ​​shear-thinning​​. Where the stress is high—for example, deep within a glacier or near its rocky bed—the effective viscosity plummets, and the ice becomes "softer" and more fluid. Conversely, where the stress is low, such as near the surface, the effective viscosity is enormous, and the ice behaves almost like a rigid solid.

This single fact explains one of the most striking features of glacier motion: ​​plug flow​​. Because the ice is so much softer where stresses are highest (at the base), most of the deformation is concentrated in a relatively thin layer of rapidly shearing ice at the bottom. The vast majority of the glacier's thickness above this layer, being under lower internal shear stress, remains incredibly stiff and is simply carried along for the ride, moving almost as a single, rigid "plug".

The Rate Factor $A$: The "Softness" of Ice

If $n$ describes how ice responds to stress, the rate factor $A$ describes its intrinsic willingness to flow. It's a measure of the ice's "softness" or "fluidity." A dimensional analysis reveals that the units of $A$ actually depend on the value of $n$, confirming that it isn't a simple viscosity. So what controls this softness?

The overwhelming factor is ​​temperature​​. Creep is a thermally activated process, where molecules in the ice crystal lattice jostle and rearrange themselves. The warmer the ice, the more energetic this jostling becomes, and the more easily the crystals can deform. This relationship is not linear; it is exponential, following an ​​Arrhenius-type law​​ similar to that governing chemical reaction rates.

Here, $Q$ is an activation energy, $R$ is the gas constant, and $T$ is the absolute temperature. The crucial part is the exponential dependence. A small change in temperature can cause a huge change in the softness $A$, and thus in the flow rate.

The numbers are staggering. A seemingly minor warming of just $5\,^{\circ}\text{C}$ (from $-10\,^{\circ}\text{C}$ to $-5\,^{\circ}\text{C}$) can reduce the effective viscosity of ice by about 40%, making it flow significantly faster under the same load. A temperature increase of $20\,^{\circ}\text{C}$ (from a very cold $240\,\text{K}$ to a near-melting $260\,\text{K}$) can increase the strain rate by a factor of ten. This extreme temperature sensitivity is the fundamental physical mechanism linking climate change directly to the accelerated flow of glaciers and ice sheets.

We have seen that the flow of ice is exquisitely sensitive to its temperature. But the connection is a two-way street. The very act of deformation generates heat. As layers of ice shear past one another, the work done by the internal stresses is converted into thermal energy, a process called ​​viscous dissipation​​ or ​​strain heating​​. The rate of heat generated per unit volume, $\dot{q}$, is the product of stress and strain rate: $\dot{q} = \tau_e \dot{\epsilon}_e$.

By substituting Glen's Law into this expression, we uncover another power law:

With $n \approx 3$, the heat generation scales with stress to the fourth power. This means that regions of high stress not only deform much faster, but they also generate heat at an even more prodigious rate.

This creates a powerful positive feedback loop, a ​​thermoviscous feedback​​:

1. High stress causes rapid deformation.
2. Rapid deformation generates significant heat.
3. The generated heat warms the ice.
4. Warmer ice is softer (its rate factor $A$ increases).
5. Softer ice deforms even faster under the same stress, which generates even more heat.

This loop can, under certain conditions, spiral out of control in a phenomenon known as ​​thermal runaway​​. Imagine a localized zone of high shear within a glacier. It's a constant battle: the feedback loop generates heat, while thermal conduction tries to diffuse that heat away into the colder surrounding ice. If the shear zone is wide enough, heat cannot escape as fast as it is being produced. The temperature in the zone will rise, making the ice softer and softer, localizing the deformation further and concentrating the heating. A stability analysis shows that instability occurs when the temperature sensitivity of heating becomes greater than the rate at which heat can be diffused away, a threshold that depends critically on the width of the shear zone. This self-reinforcing process is thought to be responsible for creating the intensely localized, fast-flowing, and relatively warm shear margins that flank major ice streams.

What happens when the ice temperature reaches the pressure-melting point? At this stage, the ice is called ​​temperate ice​​, and the physics of its flow changes once again. The temperature is now buffered—any extra heat goes into melting a small fraction of the ice rather than raising the temperature further. This shuts down the strong Arrhenius temperature dependence of the softness parameter $A$.

A new factor now takes center stage: the presence of liquid water at the boundaries between ice crystals. This water acts as a lubricant, allowing grains to slide past each other far more easily, dramatically increasing the ice's softness. The rate factor $A$ is no longer primarily a function of temperature, but of water content and pressure. The higher the water pressure within the ice, the more it pushes back against the immense weight of the overlying ice, reducing the friction between grains. This is captured by the concept of ​​effective pressure​​—the difference between the ice overburden pressure and the internal water pressure. In temperate ice, the softness $A$ increases as the effective pressure decreases. This provides a direct link between the glacier's internal plumbing and its large-scale dynamics, and it is fundamental to explaining the behavior of fast-flowing, warm-based glaciers all over the world.

Now that we have acquainted ourselves with the principles behind Glen's flow law, we are ready for the real fun. Like any profound physical law, its true beauty is revealed not in isolation, but in its power to connect disparate ideas and to solve problems in the real world. We are about to embark on a journey that will take us from the heart of a glacier to the cold, distant plains of Pluto, all with this simple power law as our guide.

How does one even begin to model a continent-spanning ice sheet? The most complete, "gold-standard" approach is to treat the ice as a very thick, non-Newtonian fluid and write down the full equations of motion—a system known in fluid dynamics as the ​​Full Stokes equations​​. These equations represent a perfect balance of forces at every single point within the ice: the crushing force of gravity is balanced by internal pressure gradients and the resistance from viscous stresses. And what is the crucial ingredient that tells the equations how ice behaves under stress? It is, of course, Glen's flow law, which provides the ​​constitutive relation​​ linking the strain rate to the stress, defining the ice's unique, non-linear viscosity.

Solving the full-Stokes system is a formidable task, a bit like trying to track the motion of every molecule in a flowing river. For many purposes, especially when looking at the vast, slow-moving interiors of Greenland and Antarctica, physicists and glaciologists have developed brilliant simplifications. The most famous of these is the ​​Shallow Ice Approximation (SIA)​​. The SIA recognizes that for a very wide and relatively thin (in a geometric sense) ice sheet, the dominant action is vertical shear—layers of ice sliding over one another. It wisely neglects other, smaller stresses. By doing so, the complex Stokes equations collapse into a much more manageable form. With the SIA, we can do amazing things, like derive an explicit formula for the total volume of ice flowing past a certain line per unit time, known as the ice flux. This flux is the lifeblood of an ice sheet, and its calculation is a cornerstone of the Earth System Models used to predict our planet's climate future.

These models also benefit from another powerful tool in the physicist's arsenal: ​​scaling analysis​​. Before diving into a brute-force calculation, we can ask: what are the essential ingredients controlling the flow? By non-dimensionalizing the equations, we can group various physical parameters—like density $\rho$, gravity $g$, ice thickness $H$, and slope $\alpha$—into a single characteristic velocity scale. This process reveals the deep mathematical structure of the problem, showing, for instance, how the surface velocity scales with the $(n+1)$-th power of the ice thickness. This gives us profound physical intuition without solving the equations for every possible glacier on Earth.

Writing down beautiful equations is one thing; solving them for a real-world glacier with a bumpy bed and complex boundaries is another entirely. This is where the power of computational science comes in. The modern approach to solving the full-Stokes equations is a numerical technique called the ​​Finite Element Method (FEM)​​.

To make the problem digestible for a computer, mathematicians and engineers reformulate the equations from their "strong" form (which must hold at every point) into a "weak" form. This new form expresses the force balance in an averaged sense over small volumes, or "elements". This clever trick not only makes the problem solvable but also elegantly incorporates the physical conditions at the glacier's boundaries. Deriving this weak form is a crucial step that directly translates the physics of Glen's flow law into a language that computational algorithms can understand, forming the very foundation of modern ice sheet simulation software.

The influence of Glen's law is not confined to the grand scales of planetary ice sheets. It is just as critical for human-scale endeavors. Imagine the challenge of building a scientific research station in Antarctica, where the entire foundation is a sheet of moving, deforming ice. Or consider excavating a tunnel to access a subglacial lake. What happens to the tunnel after you've dug it?

The immense pressure from the overlying ice will cause it to slowly close. The rate of this closure is dictated precisely by Glen's flow law. Engineers must calculate this creep to ensure the tunnel remains open and safe for its intended lifespan. This application connects glaciology to geotechnical engineering, using concepts like the ​​Ground Response Curve (GRC)​​ to relate the tunnel's closure rate to the internal support pressure provided. Furthermore, a fascinating feedback loop emerges: the viscous deformation itself generates heat, warming the ice wall. This warming softens the ice (the parameter $A$ in the law is temperature-dependent), causing it to deform even faster, which generates more heat! If this process runs away, it can even cause melting at the tunnel wall, creating a fully coupled thermo-viscoplastic problem of remarkable complexity and practical importance.

A glacier is not a static block of ice. At its base, there is often a complex network of conduits and cavities, a subglacial plumbing system that transports meltwater. The behavior of this drainage system is a key controller of how fast a glacier moves, as high water pressure can lubricate the bed and accelerate sliding.

Here, Glen's law plays a starring, and perhaps surprising, role. It governs the battle between the opening and closing of these water channels. On one side, the friction of the flowing water generates heat, melting the channel walls and enlarging the passage. On the other, the crushing weight of the glacier causes the ice to creep inwards, squeezing the channel shut. The rate of this closure is, once again, determined by Glen's flow law. The equilibrium size of a subglacial channel is a delicate balance between these two competing processes. The transition from an inefficient, distributed system of cavities to an efficient, channelized one can dramatically lower water pressure, increase the effective pressure on the bed, and ironically, slow the glacier down by increasing friction. This beautiful interplay connects continuum mechanics, hydrology, and thermodynamics in the dark, hidden world beneath the ice.

All these models are wonderful, but how do we know they are right? And how can we use them to predict the future, particularly in the context of climate change and sea-level rise?

The answer is data. Satellites now continuously monitor the Earth's ice sheets, providing us with unprecedented maps of surface velocity. These observations are a reality check for our models. Scientists use sophisticated statistical methods, such as ​​Bayesian inference​​, to find the optimal values for the parameters in Glen's flow law—most notably the softness factor $A$—that make the model's output best match the satellite data. This process, known as data assimilation, closes the loop between theory and observation, allowing us to calibrate our models and make them truly predictive tools.

This predictive power is nowhere more critical than at the "grounding line," the zone where a glacier lifts off the bedrock and begins to float on the ocean. The stability of this zone is the Achilles' heel of marine-terminating ice sheets and a primary source of uncertainty in future sea-level projections. Simplified models, built on the foundations of Glen's law and principles of hydrostatic flotation, allow us to simulate how the grounding line advances and retreats in response to changes in climate, such as warming ocean waters or changing snowfall patterns. These models demonstrate the profound sensitivity of ice discharge to temperature, a direct consequence of the Arrhenius relation embedded within the flow law's softness parameter $A(T)$.

Finally, let us lift our gaze from our own world. The laws of physics are universal. While on Earth we are concerned with water ice, other worlds in our solar system have glaciers made of entirely different materials. On Pluto, vast mountains of solid nitrogen flow to fill basins. On Saturn's moon Titan, river valleys may have been carved by glaciers of solid methane. Planetary scientists apply the very same principles, using Glen's flow law with parameters appropriate for these "exotic" ices, to understand the landscapes and geological processes on these distant worlds. They even model complex feedbacks, where latent heat absorbed or released during phase transitions can dramatically alter the ice's softness and its erosive power. What began as a law to describe terrestrial glaciers has become a universal tool for exploring the geodynamics of the cosmos.

From a simple empirical rule, we have built a framework that spans engineering, climate science, and planetary geology. We have seen how it lives inside the most fundamental equations of motion and how it breathes life into the most advanced computer simulations. This is the hallmark of a truly great physical law: its ability to unify, to explain, and to empower our quest for understanding.