Ice-crystal fabric can induce mechanical anisotropy that significantly affects flow, but ice-flow models generally do not include fabric development or its effect upon flow. Here, we incorporate a new spectral expansion of fabric, and more complete description of its evolution, into the ice-flow model Elmer/Ice. This approach allows us to model the effect of both lattice rotation and migration recrystallization on large-scale ice flow. The fabric evolution is coupled to flow using an unapproximated non-linear orthotropic rheology that better describes deformation when the stress and fabric states are misaligned. These improvements are most relevant for simulating dynamically interesting areas, where recrystallization can be important, tuning data are scarce and rapid flow can lead to misalignment between stress and fabric. We validate the model by comparing simulated fabric to ice-core and phase-sensitive radar measurements on a transect across Dome C, East Antarctica. With appropriately tuned rates for recrystallization, the model is able to reproduce observations of fabric. However, these tuned rates differ from those previously derived from laboratory experiments, suggesting a need to better understand how recrystallization acts differently in the laboratory compared to natural settings.

We compare different ways the bulk flow nonlinearity of glacier ice can be captured in an orthotropic rheology. Specifically, we compare the unapproximated orthotropic rheology, derived from plastic potential theory, to existing approximations that assume either the nonlinear viscosity or fluidity is identical to that of Glen's isotropic flow law. We find, overall, a reasonable agreement between the three orthotropic rheologies, and with existing Dye 3 ice-core deformation tests, although assuming Glen's viscosity provides the best approximation to the unapproximated rheology. Our results therefore suggest that previous studies based on either approximation to the orthotropic rheology are on relatively safe ground in the sense that both approximations generally agree with the unapproximated rheology and experimental data. Finally, we provide the forward and inverse analytical forms of all three rheologies for use in future numerical ice-flow modelling.

We investigate the errors caused by neglecting the crystal-orientation fabric when inferring the basal friction coefficient field, and whether such errors can be alleviated by inferring an isotropic enhancement factor field to compensate for missing fabric information. We calculate the steady states that arise from ice flowing over a sticky spot and a bedrock bump using a vertical-slab numerical ice-flow model, consisting of a Weertman sliding law and the anisotropic Johnson flow law, coupled to a spectral fabric model of lattice rotation and dynamic recrystallisation. Given the steady or transient states as input for a canonical adjoint-based inversion, we find that Glen's isotropic flow law cannot necessarily be used to infer the true basal drag or friction coefficient field, which are obscured by the orientation fabric, thus potentially affecting vertically integrated mass fluxes. By inverting for an equivalent isotropic enhancement factor, a more accurate mass flux can be recovered, suggesting that joint inversions for basal friction and the isotropic flow-rate factor may be able to compensate for mechanical anisotropies caused by the fabric. Thus, in addition to other sources of rheological uncertainty, fabric might complicate attempts to relate subglacial conditions to basal properties inferred from an inversion relying on Glen's law.

Streaming ice accounts for a major fraction of global ice flux, yet we cannot yet fully explain the dominant controls on its kinematics. In this contribution, we use an anisotropic full-Stokes thermomechanical flow solver to characterize how mechanical anisotropy and temperature distribution affect ice flux. For the ice stream and glacier geometries we explored, we found that the ice flux increases 1–3% per °C temperature increase in the margin. Glaciers and ice streams with crystallographic fabric oriented approximately normal to the shear plane increase by comparable amounts: an otherwise isotropic ice stream containing a concentrated transverse single maximum fabric in the margin flows 15% faster than the reference case. Fabric and temperature variations independently impact ice flux, with slightly nonlinear interactions. We find that realistic variations in temperature and crystallographic fabric both affect ice flux to similar degrees, with the exact effect a function of the local fabric and temperature distributions. Given this sensitivity, direct field-based measurements and models incorporating additional factors, such as water content and temporal evolution, are essential for explaining and predicting streaming ice dynamics.

We developed a tilt sensor for studying ice deformation and installed our tilt sensor systems in two boreholes drilled close to the shear margin of Jarvis Glacier, Alaska to obtain kinematic measurements of streaming ice. We used the collected tilt data to calculate borehole deformation by tracking the orientation of the sensors over time. The sensors' tilts generally trended down-glacier, with an element of cross-glacier flow in the borehole closer to the shear margin. We also evaluated our results against flow dynamic parameters derived from Glen's exponential flow law and explored the parameter space of the stress exponent n and enhancement factor E. Comparison with values from ice deformation experiments shows that the ice on Jarvis is characterized by higher n values than that is expected in regions of low stress, particularly at the shear margin (~3.4). The higher n values could be attributed to the observed high total strains coupled with potential dynamic recrystallization, causing anisotropic development and consequently sped up ice flow. Jarvis' n values place the creep regime of the ice between basal slip and dislocation creep. Tuning E towards a theoretical upper limit of 10 for anisotropic ice with single-maximum fabric reduces the n values by 0.2.

Ice rheology governs how glaciers flow and respond to environmental change. The rheology of glacier ice evolves in response to a variety of mechanisms, including damage, heating, melting and the development of crystalline fabric. The relative contributions of these rheological mechanisms are not well understood. Using remotely sensed data and physical models, we decouple the influence of each of the aforementioned mechanisms along the margins of Rutford Ice Stream, a laterally confined outlet glacier in West Antarctica. We show that fabric is an important control on ice rheology in the shear margins, with an inferred softening effect consistent with a single-maximum fabric. Fabric evolves to steady state near the onset of streaming flow, and ice progressively softens downstream almost exclusively due to shear heating. The rate of heating is sensitive to local shear strain rates, which respond to local changes in bed topography as ice is squeezed through the basal trough. The impact of shear heating on the downstream evolution of ice rheology in a laterally confined glacier suggests that the thermoviscous feedback – wherein faster ice flow leads to higher rates of shear heating, further softening the ice – is a fundamental control on glacier dynamics.

Flow at ice divides, their shape, size and internal structure depend not only on local conditions, but also on the flow regimes and past histories of the surrounding ice masses. Here we use field data from Kealey Ice Rise, Ellsworth Land, West Antarctica, in combination with flow modelling to investigate any possible signs of transients in the flow of the surrounding ice masses. Kealey Ice Rise shows linear surface features running parallel to its ridge in satellite imagery and a conspicuous layering in the ground-penetrating radar data known as double-peaked Raymond bumps. Through numerical modelling, by using an anisotropic full-Stokes thermomechanical flow solver, we analyse the evolution of Kealey Ice Rise and the timescales involved. We conclude that the features observed in the stratigraphy of Kealey Ice Rise require at least 3 ka of near-stationary flow conditions. However, we cannot exclude the possibility of a recent flow reorganization in the last century. We stress that the signs of stationary flow in radar stratigraphy observed in Kealey Ice Rise have been observed in other ice divides in the East Ellsworth Land area, suggesting stationary flow conditions over a millennial timescale in the region.

In this study we compare the anisotropic flow relations for polycrystalline ice of Azuma and Goto-Azuma (1996), Thorsteinsson (2002), Placidi and others (2010) and Budd and others (2013). Observations from the Dome Summit South (DSS) ice-coring site at Law Dome, East Antarctica, are used to model the vertical distribution of deviatoric stress components at the borehole site. The flow relations in which the anisotropic rheology is parameterized by a scalar function, so that the strain-rate and deviatoric stress tensor components are collinear, provide simple shear and vertical compression deviatoric stress profiles that are most consistent with laboratory observations of tertiary creep in combined stress configurations. Those flow relations where (1) the anisotropy is derived from the magnitude of applied stresses resolved onto the basal planes of individual grains and (2) the macroscopic deformation is obtained via homogenization of individual grain responses provide stress estimates less consistent with laboratory observations. This is most evident in combined simple shear and vertical compression flow regimes where shear is dominant. Our results highlight the difficulties associated with developing flow relations which incorporate a physically based description of microdeformation processes. In particular, this requires that all relevant microdeformation, recrystallization and recovery processes are adequately parameterized.

The orientations of individual crystals within a polycrystalline aggregate subjected to stress have a strong influence on its bulk strain rate and flow behavior. The ability to include the effect of crystal fabric and recrystallization processes in an ice flow law, especially at the bottom of glaciers and ice sheets where temperature is close to the pressure-melting point, is important because the stratigraphy of the ice body may be affected and the paleoclimate reconstruction hampered. We present herein three newly developed deformation apparatuses offering the possibility, from single experiments, of investigating different finite strain stages and their corresponding c-axis fabric and grain texture patterns in various confined, shear flow configurations (simple shear, pure shear and compression/extension bending). The technical set-ups and major advantages compared to classical methods are explained, and results from experiments are discussed in order to illustrate the functioning and purposes of the methods. In all experiments, significant variations in the microstructural development have been observed that reflect the varying orientations of the anisotropy and its relationship to the stress pattern. In monocrystalline ice-bending experiments, the pre-existing c-axis fabric is shown to have a profound influence on the response to stress and possibly to the type of slip system activated.