Anisotropy of Glaciers and Ice sheets

By Tamara Gerber

Aniso… what?

Glaciers and ice sheets move by sliding over their bed and by deforming internally under their own weight. How easily the ice deforms depends not only on temperature and stress, but also on the way the ice crystals are arranged. When these crystals become aligned in a preferred orientation, the ice behaves differently depending on direction—a property known as anisotropy.

A simple example of anisotropy is a cat’s fur. Stroke it from head to tail and it feels smooth; stroke it the other way and it feels rougher (and the cat might be less happy). Ice can behave in a similar way: it deforms more easily in some directions than in others.

Because ice flows to a large degree by deforming internally, this directional behaviour can affect both the speed and pattern of glacier flow, and therefore how we understand and predict the behaviour of glaciers and ice sheets.

From crystals to flowing ice

Ice crystals are inherently anisotropic because of their molecular structure. Water molecules form hexagonal rings stacked in layers. The plane of these rings is called the basal plane, and the direction perpendicular to it the c-axis (Fig. 1). Ice deforms most easily by sliding along the basal planes, like a deck of cards. In addition to this mechanical anisotropy, ice crystals are also anisotropic in their thermal and electrical properties [1, 2, 3]. In particular the dielectric anisotropy plays a crucial role in how we can measure anisotropy in glacier ice.

Figure 1: In natural ice on Earth, water molecules are arranged in layers of hexagonal rings. Because of this crystal structure, deformation occurs a lot more easily by shearing along the so-called basal plane compared to other directions, so ice is mechanically anisotropic. Additionally the crystal arrangement also explains other forms of anisotropy, for example in the dielectric permittivity which affects the propagation of electromagnetic waves (e.g. light or radar). After Libbrecht [4] and Cuffey and Paterson [5].

When snow first falls, ice crystals are randomly oriented and their anisotropic properties average out, so the bulk ice behaves almost isotropically (equally in all directions). As snow is transformed to ice and transported by glacier flow, crystals gradually rotate and re-crystallize [6, 7]. Over time, many become aligned, producing an anisotropic crystal orientation fabric (COF), or simply fabric, much like a cat’s fur becomes directional due to the alignment of individual hairs.

The fabric describes how crystal c-axes are distributed in space and is often visualized using stereoplots (Fig. 2). Because the crystal orientation is a result of past deformation it records ice-flow history [8, 9], but also influences how ice deforms today[10]. Measuring crystal orientation in ice sheets therefore helps scientists reconstruct how glaciers and ice sheets flowed in the past and evolved to their present state. At the same time, these observations are used to improve and validate computer models that predict how ice sheets will respond to future climate change and how much they will contribute to sea-level rise.

Figure 2: Whether glacier ice is anisotropic or not depends on how the individual crystals are arranged in the bulk: for a random orientation, the anisotropic effects of individual crystals cancel each other out, while if crystals are aligned in preferred directions, the bulk ice becomes anisotropic. Stereoplots are a helpful means of visualizing the crystals orientation in space.

Measuring fabric in ice cores

Ice cores provide the most direct observations of crystal fabric. Thin sections of ice are analyzed under polarized light, where each crystal appears differently depending on its orientation. Modern fabric analyzers can automatically measure thousands of crystals, revealing how fabric changes with depth [11]. Most ice cores show a transition from randomly oriented crystals near the surface to strongly aligned fabrics at depth.

The resulting patterns depend on the local stress regime: vertical compression produces vertical alignment, while shear leads to c-axis alignment perpendicular to the shear plane, and flow acceleration produces girdle-like patterns with c-axes aligned on a plane perpendicular to ice flow [14, 15, 16, 17].

While ice cores are crucial to understand the ice fabric in detail, they can only provide point measurements and are typically drilled in slow-flowing regions, limiting their ability to capture spatial variability.

Figure 3: The crystal orientation fabric can be directly measured from thin section of ice cores and using a fabric analyzer [12] that measures the orientation of the crystals with polarimetric light. This type of analysis produces images and stereoplots like shown from the Neem ice core [13, https://doi.org/10.5194/tc-8-1129-2014].

Mapping fabric across ice sheets with radar

To extend observations beyond ice cores, scientists use a method called radio-echo sounding (or ice-penetrating radar)[19]. Radar systems send pulses of polarized radio waves into the ice and record echoes from internal layers and the bed. Mounted on aircrafts or sledges, they can survey large regions rapidly.

Because ice crystals are dielectrically anisotropic, radar waves travel at different speeds depending on their polarization (the direction in which the electromagnetic field oscillates) relative to the crystal fabric [2]. Even small differences in wave speed accumulate over the long distances through km-thick ice, causing reflections from the same interface to return with a measurable time delay for different wave polarizations (Fig. 4a,c). This time delay can be used to infer the strength and orientation of the crystal fabric when measurements are made with multiple wave polarizations.        

More advanced phase-sensitive radars (pRES) can also track subtle phase changes (i.e. small shifts in the timing of the wave’s oscillation cycle) in the returning signal, in addition to the travel-time delay, allowing fabric to be detected even in thinner ice or where it is weaker [20]. Those phase shifts sometimes also produce  birefringence or beat patterns in radargrams, where the signal oscillates due to the interference between two wave components drifting out of phase (Fig. 4b,d), causing distinct stripy patterns that can be used to estimate the fabric strength [21].

Figure 4: Radar waves polarized in different directions travel at slightly different velocities in anisotropic ice. At radar profile cross-over points (a), this causes reflections from the same internal layer to arrive at different times (c). Anisotropy also produces characteristic birefringence, or ”beat”, signatures in radargrams (b), which arise because the radar wave splits into two orthogonal components that gradually drift out of phase as they propagate through the ice. After Gerber et al. [18, https://doi.org/10.1038/s41467-023-38139-8].

From past ice deformation to future change

Radar observations have shown that the fabric in fast-flowing regions such as Thwaites Glacier [22], Whillans Ice Stream[23], and the Northeast Greenland Ice Stream[18]  are strongly horizontally anisotropic, and seems to play a critical role in the dynamics of these regions. The fabric tends to stiffen ice in flow direction while facilitating shear along margins (Fig. 5), helping sustain the fast flow that is difficult to explain without anisotropy [18].

Figure 5: In the Northeast Greenland Ice Stream (NEGIS), radar-derived (a) and modeled (d) crystal fabrics show good agreement with observations from ice cores. This convergence between methods increases confidence in our ability to infer ice fabric and assess its influence on deformation. Within the ice stream, crystal orientation fabric facilitates horizontal shearing, particularly in the shear margins (b,e), helping to sustain fast ice flow. At the same time, it makes along-flow compression and extension more difficult than in isotropic ice, with potential implications for how far inland dynamic perturbations can propagate. Adapted from Gerber et al.  [18, https://doi.org/10.1038/s41467-023-38139-8].

Ice-flow models are now increasingly able to simulate fabric evolution explicitly [24, 25]. Together with growing computational power and better observational constraints, this allows anisotropic ice behavior to be included in large-scale simulations used to study ice flow. However, uncertainties remain because radar and field data only partially constrain the full three-dimensional fabric, and converting fabric into ice stiffness still relies on simplifying assumptions [21, 26, 27, 28, 29].

Despite these limitations, the combination of ice cores, radar observations, and modeling is rapidly improving our understanding of the internal structure and dynamics of ice sheets. This is reshaping how ice flow is represented in models and will lead to better representation of fast-flowing regions such as ice streams, and ultimately to better predictions of future ice-sheet evolution and sea-level rise.

Further reading

Two EOS articles on radar polarimetry and how this method can be used to measure fabric in ice sheets:

https://eos.org/science-updates/new-directions-in-mapping-ice-sheet-fabrics-and-flow

https://eos.org/editors-vox/how-radar-reveals-the-hidden-fabric-of-ice-sheets

About the Author

Dr Tamara Gerber (she/her) is a post-doc at the Université de Lausanne, Switzerland.

I am Tamara Gerber, a geophysicist and glaciologist specialized in the use of radio-echo sounding to investigate glaciers and ice sheets. My research focuses on understanding their past, present, and future dynamics, and their role in the Earth system. I am particularly interested in addressing these questions through an interdisciplinary approach that combines field observations, remote sensing, and numerical modeling. I am currently a Postdoc at the University of Lausanne.

Beyond research, I am passionate about science communication and enjoy explaining complex scientific concepts in creative and accessible ways, making glaciology engaging and approachable for diverse audiences.

References

[1] J. Klinger and G. Rochas. “Anisotropic heat conduction of fresh hexagonal ice single crystals at low temperature”. In: Journal of Physics C: Solid State Physics 15.21 (July 1982), p. 4503. doi: 10.1088/0022-3719/15/21/014.

[2] S. Fujita and S. Mae. “Relation between ice sheet internal radio-echo reflections and ice fabric at Mizuho Station, Antarctica”. In: Annals of Glaciology 17 (1993), pp. 269–275. doi: 10.3189/S0260305500012957.

[3] T. Matsuoka et al. “Precise measurement of dielectric anisotropy in ice Ih at 39 GHz”. In: Journal of Applied Physics 81.5 (Mar. 1997), pp. 2344–2348. issn: 0021-8979. doi: 10.1063/1.364238.

[4] K. G. Libbrecht. “The physics of snow crystals”. In: Reports on Progress in Physics 68.4 (Mar. 2005), p. 855. doi: 10.1088/0034-4885/68/4/R03.

[5] K. M. Cuffey and W. S. B. Paterson. The Physics of Glaciers. 4th ed. Burlington, MA: Academic Press, 2010. isbn: 9780123694614. doi: 10.1016/C2009-0-14802-X.

[6] A. J. GOW and T. WILLIAMSON. “Rheological implications of the internal structure and crystal fabrics of the West Antarctic ice sheet as revealed by deep core drilling at Byrd Station”. In: GSA Bulletin 87.12 (Dec. 1976), pp. 1665–1677. issn: 0016-7606. doi: 10.1130/0016-7606(1976)87<1665:RIOTIS>2.0.CO;2.

[7] R. B. Alley. “Fabrics in Polar Ice Sheets: Development and Prediction”. In: Science 240.4851 (1988), pp. 493–495. doi: 10.1126/science.240.4851.493.

[8] D. A. Lilien et al. “Modeling Ice-Crystal Fabric as a Proxy for Ice-Stream Stability”. In: Journal of Geophysical Research: Earth Surface 126.9 (2021), e2021JF006306. doi: 10.1029/2021JF006306.

[9] M.-G. Llorens et al. “Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?” In: The Cryosphere 16.5 (2022), pp. 2009–2024. doi: 10.5194/tc-16-2009-2022.

[10] P. Duval, M. Ashby, and I. Anderman. “Rate-controlling processes in the creep of polycrystalline ice”. In: The Journal of Physical Chemistry 87.21 (1983), pp. 4066–4074. doi: 10.1021/j100244a014.

[11] C. J. L. Wilson, D. S. Russell-Head, and H. M. Sim. “The application of an automated fabric analyzer system to the textural evolution of folded ice layers in shear zones”. In: Annals of Glaciology 37 (2003), pp. 7–17. doi: 10.3189/172756403781815401.

[12] M. Peternell et al. “Quantification of the microstructural evolution of polycrystalline fabrics using FAME: Application to in situ deformation of ice”. In: Journal of Structural Geology 61 (2014). Microdynamics of Ice, pp. 109–122. issn: 0191-8141. doi: 10.1016/j.jsg.2013.05.005.

[13] M. Montagnat et al. “Fabric along the NEEM ice core, Greenland, and its comparison with GRIP and NGRIP ice cores”. In: The Cryosphere 8.4 (2014), pp. 1129–1138. doi: 10.5194/tc-8-1129-2014.

[14] T. Thorsteinsson, J. Kipfstuhl, and H. Miller. “Textures and fabrics in the GRIP ice core”. In: Journal of Geophysical Research: Oceans 102.C12 (1997), pp. 26583–26599. doi: 10.1029/97JC00161.

[15] Y. Wang et al. “A vertical girdle fabric in the NorthGRIP deep ice core, North Greenland”. In: Annals of Glaciology 35 (2002), pp. 515–520. doi: 10.3189/172756402781817301.

[16] R. E. Thomas et al. “Microstructure and Crystallographic Preferred Orientations of an Azimuthally Oriented Ice Core from a Lateral Shear Margin: Priestley Glacier, Antarctica”. In: Frontiers in Earth Science Volume 9 – 2021 (2021). issn: 2296-6463. doi: 10.3389/feart.2021.702213.

[17] N. Stoll et al. “Linking crystallographic orientation and ice stream dynamics: evidence from the EastGRIP ice core”. In: The Cryosphere 19.9 (2025), pp. 3805–3830. doi: 10.5194/tc-19-3805-2025.

[18] T. A. Gerber et al. “Crystal orientation fabric anisotropy causes directional hardening of the Northeast Greenland Ice Stream”. In: Nature Communications 14.1 (2023), p. 2653. doi: 10 . 1038 /s41467-023-38139-8.

[19] B. H. Hills et al. “Radar Polarimetry in Glaciology: Theory, Measurement Techniques, and Scientific Applications for Investigating the Anisotropy of Ice Masses”. In: Reviews of Geophysics 63.4 (2025), e2024RG000842. doi: 10.1029/2024RG000842.

[21] O. Zeising et al. “Improved estimation of the bulk ice crystal fabric asymmetry from polarimetric phase co-registration”. In: The Cryosphere 17.3 (2023), pp. 1097–1105. doi: 10.5194/tc-17-1097-2023.

[21] S. Fujita, H. Maeno, and K. Matsuoka. “Radio-wave depolarization and scattering within ice sheets: a matrix-based model to link radar and ice-core measurements and its application”. In: Journal of Glaciology 52.178 (2006), pp. 407–424. doi: 10.3189/172756506781828548.

[22] T. J. Young et al. “Inferring Ice Fabric From Birefringence Loss in Airborne Radargrams: Application to the Eastern Shear Margin of Thwaites Glacier, West Antarctica”. In: Journal of Geophysical Research: Earth Surface 126.5 (2021), e2020JF006023. doi: 10.1029/2020JF006023.

[23] T. M. Jordan et al. “Estimation of ice fabric within Whillans Ice Stream using polarimetric phase-sensitive radar sounding”. In: Annals of Glaciology 61.81 (2020), pp. 74–83. doi: 10.1017/aog.2020.6.

[24] D. A. Lilien et al. “Simulating higher-order fabric structure in a coupled, anisotropic ice-flow model: application to Dome C”. In: Journal of Glaciology 69.278 (2023), pp. 2007–2026. doi: 10.1017/jog.2023.78.

[25] D. H. Richards et al. “Bridging the Gap Between Experimental and Natural Fabrics: Modeling Ice Stream Fabric Evolution and its Comparison With Ice-Core Data”. In: Journal of Geophysical Research: Solid Earth 128.11 (2023), e2023JB027245. doi: 10.1029/2023JB027245.

[26] M. R. Ershadi et al. “Polarimetric radar reveals the spatial distribution of ice fabric at domes and divides in East Antarctica”. In: The Cryosphere 16.5 (2022), pp. 1719–1739. doi: 10.5194/tc-16-1719-2022.

[27] N. M. Rathmann et al. “On the Limitations of Using Polarimetric Radar Sounding to Infer the Crystal Orientation Fabric of Ice Masses”. In: Geophysical Research Letters 49.1 (2022), e2021GL096244.doi: 10.1029/2021GL096244.

[28] F. Gillet-Chaulet et al. “A user-friendly anisotropic flow law for ice-sheet modeling”. In: Journal of Glaciology 51.172 (2005), pp. 3–14. doi: 10.3189/172756505781829584.

[29] T. M. Jordan et al. “Radar Characterization of Ice Crystal Orientation Fabric and Anisotropic Viscosity Within an Antarctic Ice Stream”. In: Journal of Geophysical Research: Earth Surface 127.6 (2022), e2022JF006673. doi: 10.1029/2022JF006673.

Leave a Reply

Your email address will not be published. Required fields are marked *

This site uses Akismet to reduce spam. Learn how your comment data is processed.

This site uses cookies. Find out more about this site’s cookies.