24 Feb 2023
 | 24 Feb 2023

A microstructure-based parameterization of the effective, anisotropic elasticity tensor of snow, firn, and bubbly ice

Kavitha Sundu, Johannes Freitag, Kévin Fourteau, and Henning Löwe

Abstract. Quantifying the link between microstructure and effective elastic properties of snow, firn, and bubbly ice is essential for many applications in cryospheric sciences. The microstructure of snow and ice can be characterized by different types of fabrics (crystallographic, geometrical) that gives rise to macroscopically anisotropic elastic behavior. While the impact of the crystallographic fabric has been extensively studied in deep firn, the present work investigates the influence of the geometrical fabric over the entire range of possible volume fractions. To this end we have computed the effective elasticity tensor of snow, firn, and ice by finite element simulations based on 395 X-ray tomography images comprising samples from the laboratory, Alps, Greenland, and Antarctica. We employed a variant of the Eshelby tensor that has been previously utilized for the parametrization of thermal and dielectric properties of snow and utilized Hashin-Shtrikman bounds to capture the nonlinear interplay between density and geometrical anisotropy. From that we derive a closed-form parametrization for all components of the (transverse isotropic) elastic tensor for all volume fractions using 2 fit parameters per tensor component. Finally we used the Thomsen parameter to compare the geometrical anisotropy to the crystallographic anisotropy in bubbly ice. While the geometrical anisotropy is clearly dominating up to ice volume fractions of Ø ≈ 0.7, a thorough understanding of elasticity in bubbly ice may require a coupled elastic theory that includes geometrical and crystallographic anisotropy.

Kavitha Sundu et al.

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on egusphere-2023-220', Pascal Hagenmuller, 04 Apr 2023
  • RC2: 'Comment on egusphere-2023-220', Antoine Wautier, 05 Jun 2023
  • RC3: 'Comment on egusphere-2023-220', Kris Houdyshell, 26 Jul 2023
  • RC4: 'Comment on egusphere-2023-220', Anonymous Referee #4, 26 Jul 2023

Kavitha Sundu et al.

Kavitha Sundu et al.


Total article views: 517 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
350 141 26 517 9 9
  • HTML: 350
  • PDF: 141
  • XML: 26
  • Total: 517
  • BibTeX: 9
  • EndNote: 9
Views and downloads (calculated since 24 Feb 2023)
Cumulative views and downloads (calculated since 24 Feb 2023)

Viewed (geographical distribution)

Total article views: 501 (including HTML, PDF, and XML) Thereof 501 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 01 Oct 2023
Short summary
Ice crystals often show a rod-like, vertical orientation in snow and firn, they are said to be anisotropic. The stiffness in vertical direction therefore differs from the horizontal which e.g. impacts the propagation of seismic waves. To quantify this anisotropy we conducted finite element simulations of 395 snow, firn, and ice core microstructures obtained from X-ray tomography. From the results we derived a parametrization that may be employed for advanced seismic studies in polar regions.