🔢 mathematics

Diffuse Adiabatic Flows in Thermally Coupled Grounded Shallow Ice Sheets: Modelling and Analysis

This paper proposes a novel thermodynamical model coupling the evolution of grounded shallow ice sheet thickness and internal temperature under physical constraints, and establishes the existence of solutions for this formal model using the penalty method, despite the challenges posed by degenerate equations and low solution regularity.

Original authors: Paolo Piersanti

Published 2026-06-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Paolo Piersanti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a massive, slow-moving river of ice, like the Greenland Ice Sheet, resting on a flat, frozen floor. This paper is a mathematical recipe for predicting two things simultaneously: how thick the ice gets (or shrinks) over time, and how hot or cold the ice is inside.

The author, Paolo Piersanti, has built a new mathematical model to describe this dance between the ice's shape and its temperature. Here is the breakdown of what he did, using simple analogies.

1. The Two Main Characters: The Ice Sheet and Its Temperature

Think of the ice sheet as a giant, flexible blanket lying on the floor.

The Shape (Thickness): The blanket can pile up high in some spots (accumulation) or melt away in others (ablation). However, it has a hard rule: it can never go below the floor. If the math tries to make the ice thickness negative, the model stops it at zero. This is like a "floor constraint."
The Temperature: Inside the blanket, heat moves around. It comes from the ground below, the air above, and the friction of the ice sliding. The ice has a "melting point" limit. The model assumes the ice stays frozen, so the temperature never gets hot enough to turn it into water. This is another "ceiling constraint."

2. The Big Problem: A Moving Room

Usually, when mathematicians study heat moving through a material, they assume the material is in a fixed box. But here, the "box" is the ice itself.

As the ice grows or shrinks, the room where the temperature lives changes size and shape every second.
It's like trying to measure the temperature inside a balloon while someone is constantly blowing air into it or letting air out, changing the balloon's size in real-time. This makes the math incredibly difficult because the boundaries are moving.

3. The Author's Solution: The "Diffuse" Trick

To solve this moving-room problem, the author uses a clever trick inspired by "diffuse-interface" models (think of it like blurring the sharp edge of a photo).

The Cut-off: Instead of trying to perfectly track the exact, jagged edge of the ice, the model introduces a "corrector" or a "buffer zone" near the surface.
The Analogy: Imagine the ice sheet is covered in a thin layer of dust or debris. The model assumes the temperature inside this layer behaves in a specific, simplified way. This allows the mathematician to pretend the "room" (the domain) is a fixed, solid cylinder, even though the ice inside it is changing shape.
The Result: This turns a nightmare of a moving boundary problem into a more manageable problem on a fixed stage, at the cost of slightly approximating the very top layer of the ice.

4. The "Obstacle" Course

The paper treats the ice thickness and temperature as an "obstacle problem."

The Ice Thickness: Imagine a ball rolling on a surface. It wants to go down, but the floor (the bedrock) stops it. The math has to figure out where the ball touches the floor and where it floats above it.
The Temperature: Similarly, the temperature wants to rise, but the "melting point" acts as a ceiling it cannot break through.
The model couples these two: the shape of the floor affects how the temperature moves, and the temperature affects how "soft" or "stiff" the ice is, which changes how the shape evolves.

5. The Mathematical Magic: Penalty and Limits

To prove that a solution actually exists (i.e., that the math doesn't break), the author uses a method called the Penalty Method.

The Analogy: Imagine you are trying to keep a car on a road. Instead of building a wall, you put a very strong spring (a penalty) on the side of the road. If the car tries to drive off, the spring pulls it back hard.
In the math, the author adds these "springs" to force the ice thickness to stay above zero and the temperature to stay below the melting point.
The Limit: He first solves the problem with the springs. Then, he makes the springs infinitely stiff (the "limit" process). He proves that even with infinitely stiff springs, a stable solution still exists.

6. The Main Achievement

The paper doesn't just propose a model; it rigorously proves that solutions exist.

Because the equations are so complex and "degenerate" (meaning they get messy near the edges where the ice disappears), the author cannot find a "perfect" smooth solution.
Instead, he proves the existence of "very weak solutions."
The Metaphor: Think of a "smooth solution" as a perfectly polished marble statue. A "very weak solution" is like a sculpture made of clay that holds its shape but might have rough edges. It's not perfect, but it is a valid, stable object that obeys the laws of physics described in the model.

Summary

Paolo Piersanti has created a new mathematical framework that links the changing shape of a grounded ice sheet with its internal temperature. By using a "diffuse" trick to fix the moving boundaries and a "penalty" method to handle the physical constraints, he proved that this complex, coupled system has a valid mathematical solution, even if that solution isn't perfectly smooth. This provides a solid theoretical foundation for understanding how ice sheets evolve when they are frozen to the ground.

Technical Summary: Diffuse Adiabatic Flows in Thermally Coupled Grounded Shallow Ice Sheets

Problem Statement
This paper addresses the mathematical modeling of the simultaneous evolution of the thickness and internal temperature of a grounded shallow ice sheet. The core challenge lies in the coupling of two distinct physical phenomena: the mechanical evolution of the ice surface (governed by mass balance and flow dynamics) and the thermodynamic evolution of the internal temperature (governed by advection-diffusion).

The problem is formulated under specific physical constraints:

Geometric Constraint: The ice thickness must remain non-negative, resting on a flat, frozen bedrock. This transforms the thickness evolution into a time-dependent obstacle problem.
Thermal Constraint: The internal temperature is constrained to remain below the pressure melting point (frozen regime), preventing latent heat generation and phase change. This introduces a second obstacle problem for the temperature field.
Moving Domain: The temperature equation is posed on a domain whose vertical extent is determined by the unknown ice thickness, creating a moving boundary problem.

The author notes that previous rigorous studies (e.g., [63]) treated the thickness evolution separately or assumed flat bedrock without thermal coupling. This work extends the analysis to the fully coupled system where the ice softness (viscosity) depends on the internal temperature, which in turn depends on the ice thickness.

Methodology
The paper employs a rigorous analytical framework combining formal derivation, variational inequalities, and compactness arguments.

Formal Derivation and Reduction:
- The author derives the governing equations for the ice surface elevation and internal temperature based on the Shallow Ice Approximation and Glen's flow law.
- To handle the moving domain and the degeneracy of the $p$ -Laplacian near the ice margin (where thickness vanishes), a change of variables is introduced ( $h = u^{(p-1)/2p}$ ), following strategies by Raviart.
- A "diffuse-interface" strategy is adopted to fix the computational domain. An auxiliary temperature field is constructed to extend the solution from the ice region to a fixed cylinder $\Omega \times (0, L)$ . This involves a cut-off function near the ice surface, effectively modeling the system as if the ice were covered by a debris layer where temperature vanishes near the surface. This allows the problem to be posed on a fixed domain, albeit with an approximation near the boundary.
Penalization and Discretization:
- The coupled obstacle problems are reformulated as a "penalized" system ( $P_\ell$ ) using penalty terms to enforce the constraints ( $u \ge 0$ and $0 \le \tilde{\Theta} \le \Theta_m$ ).
- Existence of solutions for the penalized problem is established via a semi-discrete time-stepping scheme (finite differences). The proof relies on the Minty-Browder theorem, utilizing the monotonicity and coercivity of the associated operators.
A Priori Estimates and Compactness:
- A comprehensive set of a priori estimates is derived for the discrete solutions, independent of the time step and the penalty parameter. These estimates bound the solutions in appropriate Sobolev and Bochner spaces.
- Compactness is achieved using the Aubin-Lions-Simon theorem and a discrete version of Dubinskii's lemma. Special care is taken to handle the non-linearities arising from the $p$ -Laplacian and the exponential dependence of ice softness on temperature.
Limit Passage:
- The existence of solutions for the penalized problem is used to pass to the limit as the time step tends to zero, recovering a solution for the penalized continuous problem.
- Finally, the penalty parameter $\ell$ is sent to zero. The author establishes the existence of "very weak" solutions for the original coupled system. The proof avoids extra regularity assumptions on the weak limit that were present in prior works, aligning the solution concept with the framework of evolutionary variational inequalities governed by pseudo-monotone operators.

Key Contributions and Novelties
The paper identifies three primary novelties:

Reduction to a Fixed Domain: The introduction of reduction techniques (Section 2) transforms the original model, which couples thickness and temperature on a moving domain, into an equivalent model on a fixed domain. This involves an approximation where the temperature evolution vanishes near the ice surface, characterizing ice sheets with a debris cover (adiabatic flows).
Existence of Very Weak Solutions: The paper establishes the existence of solutions for the initial boundary-value problem where both the ice thickness and internal temperature are subject to physical constraints (obstacle problems). The solution concept is "very weak," necessitated by the low regularity of the problem which prevents standard solution definitions.
Analytical Framework: The work provides a new analytical framework for a class of coupled obstacle problems posed on moving domains, extending beyond specific glaciological applications.

Main Results
The central result is Theorem 8.1, which proves the existence of very weak solutions for the initial boundary-value problem coupling the evolution of the shallow ice sheet thickness with its internal temperature, assuming a flat and frozen bedrock.

Key aspects of the result include:

The removal of the "extra decrease assumption" (4.6) required in previous analyses (e.g., [63]).
The derivation of an abstract limit model that does not rely on extra regularity assumptions for the weak limit.
The solution concept aligns with non-standard frameworks for evolutionary variational inequalities involving pseudo-monotone operators.

Significance and Claims
The paper claims to extend the rigorous mathematical study of ice sheet models initiated in [63] to the thermally coupled case. The author emphasizes that the coupling introduces significant challenges due to the moving domain and the non-linearities associated with the free boundary.

The significance of the work is framed as follows:

Modeling: It offers a novel thermodynamical model that rigorously couples surface elevation and internal temperature under physical constraints.
Analytical: It overcomes the difficulty of defining solutions on moving domains by employing a diffuse-interface approach and a change of variables to handle degeneracy.
Robustness: By dropping extra regularity assumptions, the result provides a more robust existence theory consistent with the low regularity inherent in shallow ice sheet models.

The paper concludes that the proposed model and the existence proof represent a step forward in the mathematical analysis of grounded shallow ice sheets, particularly in contexts where thermal coupling and physical constraints (frozen bedrock) are critical. The author explicitly notes that the "adiabatic" nature of the title refers to the specific thermal model choice justified by the presence of debris layers, rather than a general claim about all ice sheet dynamics.