Müller Ice Cap sits on Umingmat Nunaat (Axel Heiberg Island), Nunavut, Canada, ~ 80°N. Its high latitude and elevation suggest it experiences relatively little melt and preserves an undisturbed paleoclimate record. Here, we present a suite of field measurements, complemented by remote sensing, that constrain the ice thickness, accumulation rate, temperature, ice-flow velocity, and surface-elevation change of Müller Ice Cap. These measurements show that some areas near the top of the ice cap are more than 600 m thick, have nearly stable surface elevation, and flow slowly, making them good candidates for an ice core. The current mean annual surface temperature is −19.6 °C, which combined with modeling of the temperature profile indicates that the ice is frozen to the bed. Modeling of the depth-age scale indicates that Pleistocene ice is likely to exist with measurable resolution (300–1000 yr m−1) 20–90 m from the bed, assuming that Müller Ice Cap survived the Holocene Climatic Optimum with substantial ice thickness (~400 m or more). These conditions suggest that an undisturbed Holocene climate record could likely be recovered from Müller Ice Cap. We suggest 91.795°W, 79.874°N as the most promising drill site.

Upernavik Isstrøm, the largest contributor to sea-level rise in northwest Greenland, has experienced complex and contrasting ice-flow-speed changes across its five outlets over the last two decades. In this study, we present a detailed remote-sensing analysis of the ice dynamics at Upernavik's outlets from 2000 to 2021 to evaluate the details of these changes. Previous research suggested that the presence or absence of floating ice tongues strongly influences Upernavik's ice dynamics. We use several lines of evidence to document the presence of floating ice tongues, and find that, while several outlets experienced ice-tongue formation and/or loss during the study period, these changes do not explain observed fluctuations in ice-flow velocity. Further exploration of ice-dynamic forcings using a flowline model suggests that changes in basal slipperiness near the terminus have a strong impact on upstream ice dynamics and can explain the velocity variations. Our results suggest that speed fluctuations at Upernavik's outlets may be seasonally and interannually controlled by bed conditions near the terminus, and highlight the need for further research on the influence of basal conditions on complex tidewater glacier dynamics.

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.

Converting measurements of ice-sheet surface elevation change to mass change requires measurements of accumulation and knowledge of the evolution of the density profile in the firn. Most firn-densification models are tuned using measured depth–density profiles, a method which is based on an assumption that the density profile in the firn is invariant through time. Here we present continuous measurements of firn-compaction rates in 12 boreholes near the South Pole over a 2 year period. To our knowledge, these are the first continuous measurements of firn compaction on the Antarctic plateau. We use the data to derive a new firn-densification algorithm framed as a constitutive relationship. We also compare our measurements to compaction rates predicted by several existing firn-densification models. Results indicate that an activation energy of 60 kJ mol−1, a value within the range used by current models, best predicts the seasonal cycle in compaction rates on the Antarctic plateau. Our results suggest models can predict firn-compaction rates with at best 7% uncertainty and cumulative firn compaction on a 2 year timescale with at best 8% uncertainty.

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.

Bulk directional enhancement factors are determined for axisymmetric (girdle and single-maximum) orientation fabrics using a transversely isotropic grain rheology with an orientation-dependent non-linear grain fluidity. Compared to grain fluidities that are simplified as orientation independent, we find that bulk strain-rate enhancements for intermediate-to-strong axisymmetric fabrics can be up to a factor of ten larger, assuming stress homogenization over the polycrystal scale. Our work thus extends previous results based on simple basal slip (Schmid) grain rheologies to the transversely isotropic rheology, which has implications for large-scale anisotropic ice-flow modelling that relies on a transversely isotropic grain rheology. In order to derive bulk enhancement factors for arbitrary evolving fabrics, we expand the c-axis distribution in terms of a spherical harmonic series, which allows the rheology-required structure tensors through order eight to easily be calculated and provides an alternative to current structure-tensor-based modelling.

Despite widespread use of radio-echo sounding (RES) in glaciology and broad distribution of processed radar products, the glaciological community has no standard software for processing impulse RES data. Dependable, fast and collection-system/platform-independent processing flows could facilitate comparison between datasets and allow full utilization of large impulse RES data archives and new data. Here, we present ImpDAR, an open-source, cross-platform, impulse radar processor and interpreter, written primarily in Python. The utility of this software lies in its collection of established tools into a single, open-source framework. ImpDAR aims to provide a versatile standard that is accessible to radar-processing novices and useful to specialists. It can read data from common commercial ground-penetrating radars (GPRs) and some custom-built RES systems. It performs all the standard processing steps, including bandpass and horizontal filtering, time correction for antenna spacing, geolocation and migration. After processing data, ImpDAR's interpreter includes several plotting functions, digitization of reflecting horizons, calculation of reflector strength and export of interpreted layers. We demonstrate these capabilities on two datasets: deep (~3000 m depth) data collected with a custom (3 MHz) system in northeast Greenland and shallow (<100 m depth, 500 MHz) data collected with a commercial GPR on South Cascade Glacier in Washington.