The Importance of Ice Thermodynamics
A thorough understanding of ice thermodynamics is essential for accurately describing the behavior of glaciers, ice sheets, and ice shelves. Ice temperatures control both the rate of ice deformation and the occurrence of sliding when the base reaches the melting point. Precisely, ice softens by an order of magnitude as the temperature increases from -10°C to the melting point, and velocities can increase by 2-3 orders of magnitude over a temperate base that yields rapid sliding.
However, accurately estimating ice temperatures is challenging, as the heat transfer balance is the result of a complex interplay between advection, diffusion, and various heat sources. Only an accurate representation of these processes will allow for a robust assessment of ice flow, mass balance, and overall stability.
In this context, the development of analytical solutions for ice thermodynamics can provide deeper comprehension of the fundamental physics of ice. Analytical solutions are intuitively interpretable, reveal hidden symmetries, and serve as a verification tool or benchmark for numerical models.
Existing Analytical Solutions and their Limitations
The pioneering work of Robin (1955) and Lliboutry (1963) laid the groundwork for understanding ice-column thermodynamics in the presence of vertical advection and diffusion by providing analytical solutions for stationary scenarios. These seminal studies offered valuable insights into the steady-state behavior of ice columns subject to advective-diffusive processes. However, they did not consider the time-dependent evolution of ice temperatures, limiting their applicability to situations involving steady-state ice flow and fixed environmental conditions.
While subsequent studies have aimed to address the time-dependent nature of the problem, they have often relied on strong assumptions regarding the particular vertical velocity profile, such as linear (Robin, 1955) or quadratic (Raymond, 1983) profiles. This has resulted in a lack of independent analytical descriptions of the temperatures, as noted by Huybrechts and Payne (1996).
Furthermore, traditional approaches from both numerical and analytical perspectives have assumed the simplest heat-flux boundary condition at the ice surface: the imposition of the air temperature at the uppermost ice layer. However, knowing that glacial ice forms through snow densification, this imposition appears to be an oversimplification, given that thermal conductivity increases with density.
Advancing the Analytical Understanding of Ice Thermodynamics
To address these limitations, the current study presents an analytical formulation of the transient ice temperature equation that accounts for the temporal evolution of the temperature profile, rather than assuming an equilibrated state. This approach takes a step towards a more accurate representation of the ice thermal behavior.
The formulation of the problem considers a 1D ice column with diffusive heat transport, vertical advection, strain heat, and depth-integrated horizontal advection. Importantly, a refined top boundary condition is introduced, allowing for potential non-equilibrium temperature states across the ice-air interface by applying Newton’s law of cooling.
The analytical solution is expressed in terms of confluent hypergeometric functions following a separation of variables approach. Non-dimensionalization reduces the parameter space to five numbers that fully determine the shape of the solution at equilibrium: surface insulation, effective geothermal heat flow, the Péclet number, the Brinkman number, and a parameter capturing the depth-integrated horizontal advection.
The transient component of the solution is shown to exponentially converge to the stationary profile, with a decay time that solely depends on vertical advection and surface insulation. This provides valuable insights into the fundamental physics of heat transfer in ice, as well as a way of benchmarking numerical solvers for heat transfer in ice sheet models.
Benchmark Experiments for Numerical Solvers
To test the reliability of numerical solvers, the study presents a suite of benchmark experiments with gradually increasing complexity. These experiments capture the main physical processes for heat propagation, including diffusion, vertical advection, strain heating, and depth-integrated horizontal advection.
The analytical solutions developed in this work are used to assess the accuracy of numerical discretization schemes and the required spatial resolution. The results show that a symmetric scheme for the advective term and a three-point asymmetric scheme for the basal boundary condition best match the analytical solutions. Furthermore, a convergence study indicates that a minimum of 15 vertical points is sufficient to accurately reproduce the temperature profile.
These benchmark experiments provide a valuable tool for testing and validating numerical solvers, ensuring their reliability in representing the complex thermodynamics of ice sheets. The analytical solutions presented in this work can serve as a reference for future studies, contributing to the improvement of ice sheet modeling efforts.
Practical Implications and Applications
The analytical solutions developed in this study have several practical implications and applications:
-
Transient Temperature Modeling: The time-dependent description of ice temperatures can significantly reduce the computational demands in modeling exercises, as it overcomes the need for computationally expensive transient optimization approaches.
-
Ice Sheet Initialization: Poorly known parameter fields, such as the ice temperature, are often estimated by minimizing the mismatch between observations and model output variables. The analytical solutions can provide a more accurate initial condition for these initialization procedures, improving the overall reliability of ice sheet models.
-
Troubleshooting and Performance Optimization: The insights gained from the analytical solutions, such as the role of different physical processes and the dependency on key parameters, can inform troubleshooting and performance optimization efforts for air-cooled heat exchangers used in cold climate applications.
-
Benchmarking and Validation: The benchmark experiments proposed in this study can serve as a valuable tool for testing and validating numerical solvers employed in various ice-related applications, ensuring the accuracy and reliability of the results.
-
Educational and Training Resources: The analytical solutions and their interpretation can be used as educational and training resources, fostering a deeper understanding of the fundamental principles governing heat transfer in ice and their practical implications.
By providing these analytical solutions and benchmark experiments, the current study contributes to the advancement of our understanding of ice thermodynamics and the improvement of modeling capabilities, ultimately enhancing our ability to predict and manage the impacts of climate change on glaciers, ice sheets, and related systems.
Conclusion
This study presents an analytical formulation of the transient ice temperature equation, accounting for the combined effects of diffusion, vertical advection, strain heat, and depth-integrated horizontal advection. The analytical solutions are expressed in terms of confluent hypergeometric functions and are validated through a suite of benchmark experiments.
The key findings of this work include:
- The transient component of the solution exponentially converges to the stationary profile, with a decay time that solely depends on vertical advection and surface insulation.
- The equilibrium temperature profiles are mostly independent of the surface insulation parameter, but the transient regime shows a strongly distinct behavior.
- A symmetric scheme for the advective term and a three-point asymmetric scheme for the basal boundary condition best match the analytical solutions.
- A minimum of 15 vertical points is sufficient to accurately reproduce the temperature profile.
These analytical solutions and benchmark experiments provide valuable insights into the fundamental physics of heat transfer in ice, serve as a verification tool for numerical solvers, and have practical implications for ice sheet modeling, air-cooled heat exchanger applications, and educational resources.
By advancing the analytical understanding of ice thermodynamics, this study contributes to improving our ability to predict and manage the impacts of climate change on glaciers, ice sheets, and related systems.
