Flow laws for ice predict rates of deformation (strain) and are fundamental to modelling glacier and ice-sheet dynamics. Here we apply Bayesian inference to laboratory measurements accumulated over 70 years to constrain flow laws for ice-sheet modelling. At low strains, commonly used flow laws-derived from individual experimental datasets with narrow stress, temperature and grain-size ranges-fail to capture the full complexity of ice behaviour. We show that a multicomponent flow law that sums strain rates from different deformation mechanisms is needed to capture grain-size and temperature sensitivities observed in the full set of experiments. This multicomponent flow law is applicable to natural scenarios where the anisotropy of ice is weak or where the deformation kinematics differ from those that formed the crystallographic preferred orientation, making the ice more viscous. Low-strain flow laws, including this multicomponent flow law, have limited validity at high strain, where viscosity evolves and anisotropy develops, making ice less viscous. A one-component, grain-size insensitive flow law gives a reasonable fit to high-strain experimental data and is better suited to modelling the large-scale flow behaviour of ice sheets.

The anelastic response of planetary materials to dynamic forcing and the associated dissipation of elastic strain energy influence a wide range of geophysical processes across seismic to tidal to isostatic rebound timescales, such as the attenuation of seismic waves in Earth’s mantle and tidal heating of icy satellites orbiting large host planets. Most attenuation experiments on ice and other planetary materials are conducted at relatively low stresses, to mimic seismic attenuation and to prevent fracturing of experimental samples in the absence of confining pressure. Typically, these experiments explore a linear, amplitude-independent regime, primarily governed by mechanisms associated with grain boundaries. At sufficiently high stresses, nonlinear attenuation in planetary materials is predicted to occur above some critical strain amplitude from the motion of dislocations and becomes important, for example, in exoplanetary systems wherein highly eccentric orbits of planets with close proximity to host stars can yield enormous tidal stresses. To study the onset of nonlinear dissipation behavior of ice and the associated amplitude dependence, we conducted cyclic deformation experiments in a high-pressure gas apparatus at elevated confining pressures. Cylindrical ice samples with a uniform initial grain size of ~0.3 mm were sealed in indium jackets, pressurized to 30 MPa and loaded to a median differential stress in the range 4 to 8 MPa. The samples were then subjected to sinusoidal stress oscillations from 0.5 to 3 MPa about the median stress at frequencies of 10-3 to 5x10-3 Hz. The resulting sample strain was measured with a strain gauge sandwiched between the outer surface of the ice cylinder and the 0.4 mm-thick indium jacket. The phase lag between the applied stress and the resulting strain response was used to calculate the inverse quality factor, Q-1. For experiments conducted at 260 K, we observe a marked dependence of Q-1 on stress amplitude at three different frequencies. We also observed that Q-1 decreases, and Young’s modulus increases, with increasing forcing frequency. Our experimental results will be evaluated in the context of existing theoretical models of nonlinear attenuation.

The conformer distribution of normal-propyl cyanide is investigated using broadband chirped pulse rotational spectroscopy in the millimeter-wave regime coupled with buffer gas cooling. Here we explore the relative abundances of the anti and gauche conformers following room-temperature gas-phase injection into a 25 K buffer gas cell and compare to that which is observed following temperature programmed desorption from an ice surface, similar to the slow warm-up experienced by ice grains as they approach warmer regions within the interstellar medium. The conformer distributions observed in the gas phase from room temperature injection are then used to determine their relative energies, an important parameter needed to interpret the isomer and conformer abundances derived from astronomical observations. We find the gauche conformer to be the most stable species by similar to 97 +/- 21 cm(-1). We further examine the relative conformer abundances following ice desorption, which are distinct from those following the gas-phase introduction. The ratios measured off the ice correspond to a conformer temperature of similar to 56 K, which is much lower than their sublimation temperature of 170 K.

We present a new technique for the detection of molecules desorbed from an ice surface using broad-band millimeter-wave rotational spectroscopy. The approach permits interrogation of molecules that have undergone the slow warmup process of temperature-programmed desorption (TPD), analogous to the warmup phase of icy grains in the interstellar medium as they approach the central protostar. The detection is conformer- and isomer-specific and quantitative, as afforded by chirped-pulse rotational spectroscopy. To achieve this, we combine ice TPD with buffer gas cooling, followed by detection in the millimeter-wave regime. In this report we examine the TPD profiles of n- and i-propanol, the former of which may be in five different conformational isomeric forms, and which display distinct desorption profiles. The limited conformational isomerization and temperature-dependent relative yields of n-propanol conformers observed show that the desorption is highly conformer-specific.

Grain growth can modify the microstructure of natural ice, including the grain size and crystallographic preferred orientation (CPO). To better understand grain-growth processes and kinetics, we compared microstructural data from synthetic and natural ice samples of similar starting grain sizes that were annealed at the solidus temperature (0 ∘C) for durations of a few hours to 33 d. The synthetic ice has a homogeneous initial microstructure characterized by polygonal grains, little intragranular distortion, few bubbles, and a near-random CPO. The natural ice samples were subsampled from ice cores acquired from the Priestley Glacier, Antarctica. This natural ice has a heterogeneous microstructure characterized by a considerable number of air bubbles, widespread intragranular distortion, and a CPO. During annealing, the average grain size of the natural ice barely changes, whereas the average grain size of the synthetic ice gradually increases. These observations demonstrate that grain growth in natural ice can be much slower than in synthetic ice and therefore that the grain-growth law derived from synthetic ice cannot be directly applied to estimate the grain-size evolution in natural ice with a different microstructure. The microstructure of natural ice is characterized by many bubbles that pin grain boundaries. Previous studies suggest that bubble pinning provides a resisting force that reduces the effective driving force of grain-boundary migration and is therefore linked to the inhibition of grain growth observed in natural ice. As annealing progresses, the number density (number per unit area) of bubbles on grain boundaries in the natural ice decreases, whilst the number density of bubbles in the grain interiors increases. This observation indicates that some grain boundaries sweep through bubbles, which should weaken the pinning effect and thus reduce the resisting force for grain-boundary migration. Some of the Priestley ice grains become abnormally large during annealing. We speculate that the contrast of dislocation density amongst neighbouring grains, which favours the selected growth of grains with low dislocation densities, and bubble pinning, which inhibits grain growth, are tightly associated with abnormal grain growth. The upper 10 m of the Priestley ice core has a weaker CPO and better-developed second maximum than deeper samples. The similarity of this difference to the changes observed in annealing experiments suggests that abnormal grain growth may have occurred in the upper 10 m of the Priestley Glacier during summer warming.

During plastic deformation, strain weakening can be achieved, in part, via strain energy reduction associated with intragranular boundary development and grain boundary formation. Grain boundaries (in 2D) are segments between triple junctions, that connect to encircle grains; every boundary segment in the encircling loop has a high (>10 degrees) misorientation angle. Intragranular boundaries terminate within grains or dissect grains, usually containing boundary segments with a low (<10 degrees) misorientation angle. We analyze electron backscatter diffraction (EBSD) data from ice deformed at -30 degrees C (Th approximate to ${T}_{h}\approx $ 0.9). Misorientation and weighted Burgers vector (WBV) statistics are calculated along planar intragranular boundaries. Misorientation angles change markedly along each intragranular boundary, linking low- (<10 degrees) and high-angle (10-38 degrees) segments that exhibit distinct misorientation axes and WBV directions. We suggest that these boundaries might be produced by the growth and intersection of individual intragranular boundary segments comprising dislocations with distinct slip systems. There is a fundamental difference between misorientation axis distributions of intragranular boundaries (misorientation axes mostly confined to ice basal plane) and grain boundaries (no preferred misorientation axis). These observations suggest during progressive subgrain rotation, intragranular boundaries remain crystallographically controlled up to large misorientation angles (>>10 degrees). In contrast, the apparent lack of crystallographic control for grain boundaries suggests misorientation axes become randomized, likely due to the activation of additional mechanisms (such as grain boundary sliding) after grain boundary formation, linking boundary segments to encircle a grain. Our findings on ice intragranular boundary development and grain boundary formation may apply more broadly to other rock-forming minerals (e.g., olivine, quartz).

Kinking can accommodate significant amounts of strain during crystal plastic deformation under relatively large stresses and may influence the mechanical properties of cold planetary cryosphere. To better understand the origins, mechanisms, and microstructural effects of kinking, we present detailed microstructural analyses of coarse-grained ice (~1300 µm) deformed under uniaxial compression at -30°C. Microstructural data are generated using cryogenic electron backscattered diffraction (cryo-EBSD). Deformed samples have bimodal grain size distributions, with thin and elongated (aspect ratio ≥ 4) kink domains that develop within, or at the tips of, remnant original grains (≥ 300 µm, aspect ratio < 4). Small, equiaxed subgrains also develop along margins of remnant grains. Moreover, many remnant grains are surrounded by fine-grained mantles of small, recrystallized grains (< 300 µm, aspect ratio < 4). Together, these observations indicate that grain nucleation is facilitated by both kinking and dynamic recrystallization (via subgrain rotation). Low- (< 10°) and high-angle (mostly > 10°, many > 20°) kink bands within remnant grains have misorientation axes that lie predominantly within the basal plane. Moreover, previous studies suggest the kinematics of kinking and subgrain rotation should be fundamentally the same. Therefore, progressive kinking and subgrain rotation should be crystallographically controlled, with rotation around fixed misorientation axes. Furthermore, the c-axes of most kink domains are oriented sub-perpendicular to the sample compression axis, indicating a tight correlation between kinking and the development of crystallographic preferred orientation. Kink band densities are the highest within remnant grains that have basal planes sub-parallel to the compression axis (i.e., c-axes perpendicular to the compression axis)—these data are inconsistent with models suggesting that, if kinking is the only strain-accommodating process, there should be higher kink band densities within grains that have basal planes oblique to the compression axis (for low kink-host misorientation angles, e.g., ≤ 20°, as in this study). One way to rationalize this inconsistency between kink models and experimental observations is that kinking and dynamic recrystallization are both active during deformation, but their relative activities depend on the crystallographic orientations of grains. For grains with basal planes sub-parallel to the compression axis, strain-induced GBM is inhibited, and large intragranular strain incompatibilities can be relaxed via kinking, when other processes such as subgrain rotation recrystallization are insufficient. For grains with basal planes oblique to the compression axis, strain-induced grain boundary migration (GBM) might be efficient enough to relax the strain incompatibility via selective growth of these grains, and kinking is therefore less important. For grains with basal planes sub-perpendicular to the compression axis, kink bands are seldom observed—for these grains, the minimum shear stress required for kinking exceeds the applied compressive stress, such that kinks cannot nucleate.

Dynamic recrystallization is an important mechanical weakening mechanism during the deformation of ice, yet we currently lack robust quantitative tools for identifying recrystallized grains in the “migration” recrystallization regime that dominates ice deformation at temperatures close to the ice melting point. Here, we propose grain boundary irregularity as a quantitative means for discriminating between recrystallized (high sphericity, low irregularity) and remnant (low sphericity, high irregularity) grains. To this end, we analysed cryogenic electron backscatter diffraction (cryo-EBSD) data of deformed polycrystalline ice, to quantify dynamic recrystallization using grain boundary irregularity statistics. Grain boundary irregularity has an inverse relationship with a sphericity parameter, Ψ, defined as the ratio of grain area and grain perimeter, divided by grain radius in 2-D so that the measurement is grain size independent. Sphericity (Ψ) typically decreases with increasing grain size, up to a threshold grain size, above which Ψ either plateaus (at temperature, T< -10°C) or increases much more gradually (at T≥ -10°C). There is no apparent relationship between grain boundary sphericity and grain c-axis orientation even at very high temperatures (-4 and -5°C), where GBM dominants, suggesting little crystallographic control on the activity of grain boundary migration (GBM). Decreasing sphericity up to the threshold grain size can be explained by newly-formed, small, spherical recrystallized grains growing via strain-induced GBM and thereby developing increasingly irregular grain boundaries. We suggest that the plateau (or gradual decrease) in sphericity at larger sizes represents a population of original grains (i.e., remnant grains) that becomes increasingly irregular (at similar rates) due to the balance between GBM and nucleation. In this interpretation, the threshold grain size represents the largest grain size reached by a growing recrystallized grain by the end of an experiment. Thus, grain boundary irregularity provides a means for discriminating between recrystallized and remnant grains—a capability that is potentially useful for evaluating dynamic recrystallization processes in ice deformed at temperatures close to the melting point. The threshold grain size and experiment duration can be used to calculate the rates of recrystallization and grain size evolution associated with GBM. Grain size evolution rates are similar at high and low temperatures, suggesting similar GBM rates. Previous studies show that grain boundary mobility decreases with decreasing temperature. The driving force of GBM, on the other hand, has a positive correlation with stress, which increases with a decreasing temperature if strain rate remains unchanged. The balance between boundary mobility and driving force is likely the cause of similar GBM rates between high and low temperatures.

Strain weakening leads to the formation of high‐strain shear zones and strongly influences terrestrial ice discharge. In glacial flow models, strain weakening is assumed to arise from the alignment of weak basal planes—the development of a crystallographic preferred orientation, CPO—during flow. However, in experiments, ice strain weakening also coincides with grain size reduction, which has been invoked as a weakening mechanism in other minerals. To interrogate the relative contributions of CPO development and grain size reduction toward ice strain weakening, we deformed initially isotropic polycrystalline ice samples to progressively higher strains between −4 and −30°C. Microstructural measurements were subsequently combined with flow laws to separately model the mechanical response expected to arise from CPO development and grain size reduction. Magnitudes of strain weakening predicted by the constitutive flow laws were then compared with the experimental measurements. Flow laws that only consider grain size do not predict weakening with strain despite grain size reduction. In contrast, flow laws solely considering CPO effects can reproduce the measured strain weakening. Thus, it is reasonable to assume that strain weakening in ice is dominated by CPO development, at least under high temperature ( Th ${T}_{h}$ ≥ 0.9) and high stress (>1 MPa), like those in our experiments. We speculate that at high homologous temperatures ( Th ${T}_{h}$ ≥ 0.9), CPO development will also govern the strain weakening behavior of other viscously anisotropic minerals, like olivine and quartz. Overall, we emphasize that geodynamic and glaciological models should incorporate CPOs to account for strain weakening, especially at high homologous temperatures.

In order to better understand ice deformation mechanisms, we document the microstructural evolution of ice with increasing strain. We include data from experiments at relatively low temperatures (−20 and −30 ∘C), where the microstructural evolution with axial strain has never before been documented. Polycrystalline pure water ice was deformed under a constant displacement rate (strain rate ∼1.0×10-5 s−1) to progressively higher strains (∼ 3 %, 5 %, 8 %, 12 % and 20 %) at temperatures of −10, −20 and −30 ∘C. Microstructural data were generated from cryogenic electron backscattered diffraction (cryo-EBSD) analyses. All deformed samples contain subgrain (low-angle misorientations) structures with misorientation axes that lie dominantly in the basal plane, suggesting the activity of dislocation creep (glide primarily on the basal plane), recovery and subgrain rotation. Grain boundaries are lobate in all experiments, suggesting the operation of strain-induced grain boundary migration (GBM). Deformed ice samples are characterized by interlocking big and small grains and are, on average, finer grained than undeformed samples. Misorientation analyses between nearby grains in 2-D EBSD maps are consistent with some 2-D grains being different limbs of the same irregular grain in the 3-D volume. The proportion of repeated (i.e. interconnected) grains is greater in the higher-temperature experiments suggesting that grains have more irregular shapes, probably because GBM is more widespread at higher temperatures. The number of grains per unit area (accounting for multiple occurrences of the same 3-D grain) is higher in deformed samples than undeformed samples, and it increases with strain, suggesting that nucleation is involved in recrystallization. “Core-and-mantle” structures (rings of small grains surrounding big grains) occur in −20 and −30 ∘C experiments, suggesting that subgrain rotation recrystallization is active. At temperatures warmer than −20 ∘C, c axes develop a crystallographic preferred orientation (CPO) characterized by a cone (i.e. small circle) around the compression axis. We suggest the c-axis cone forms via the selective growth of grains in easy slip orientations (i.e. ∼ 45∘ to shortening direction) by GBM. The opening angle of the c-axis cone decreases with strain, suggesting strain-induced GBM is balanced by grain rotation. Furthermore, the opening angle of the c-axis cone decreases with temperature. At −30 ∘C, the c-axis CPO changes from a narrow cone to a cluster, parallel to compression, with increasing strain. This closure of the c-axis cone is interpreted as the result of a more active grain rotation together with a less effective GBM. We suggest that lattice rotation, facilitated by intracrystalline dislocation glide on the basal plane, is the dominant mechanism controlling grain rotation. Low-angle neighbour-pair misorientations, relating to subgrain boundaries, are more extensive and extend to higher misorientation angles at lower temperatures and higher strains supporting a relative increase in the importance of dislocation activity. As the temperature decreases, the overall CPO intensity decreases, primarily because the CPO of small grains is weaker. High-angle grain boundaries between small grains have misorientation axes that have distributed crystallographic orientations. This implies that, in contrast to subgrain boundaries, grain boundary misorientation is not controlled by crystallography. Nucleation during recrystallization cannot be explained by subgrain rotation recrystallization alone. Grain boundary sliding of finer grains or a different nucleation mechanism that generates grains with random orientations could explain the weaker CPO of the fine-grained fraction and the lack of crystallographic control on high-angle grain boundaries.

Vast quantities of solid CO2 reside in topographic basins of the south polar layered deposits (SPLD) on Mars and exhibit morphological features indicative of glacial flow. Previous experimental studies showed that CO2 ice is 1–2 orders of magnitude weaker than water ice under Martian polar conditions. Here we present data from deformation experiments on pure, fine‐grained CO2 ice, over a broader range of temperatures than previously explored (158–213 K). The experiments confirm previous observations of highly nonlinear power law creep at larger stresses, but also show a transition to a previously unseen linear‐viscous creep regime at lower stresses. We examine the viscosity of CO2 within the SPLD and predict that the CO2‐rich deposits are modestly stronger than previously thought. Nevertheless, CO2 ice flows much more readily than H2O ice, particularly on the steep flanks of SPLD topographic basins, allowing the CO2 to pond as observed.