We conducted a series of experiments that revealed the formation of mm-scale penitente structures in ice illuminated by broadband light under moderate vacuum conditions between 50 and 2000 Pa. The experimental apparatus consists of a 0.3 m diameter cylindrical vacuum chamber with a cooling jacket surrounding the outer radius and bottom surface. Light shines in through an optical window at the top to illuminate most of the ice surface. We observe penitente-like structures at temperatures between −15$^\circ$C and $-2^\circ$C and pressures close to the equilibrium vapor pressure at the ice surface temperature. The formation of these structures is very sensitive to slight changes in background pressure, and the structures tend to vanish with significant deviations away from the equilibrium curve, resulting in a smooth sublimated crater formation instead of penitentes. Application of the physical model by Claudin and others (2015, doi: 10.1103/PhysRevE.92.033015) at experimental conditions generally agrees with observations for penitente spacing.

New experiments have revealed that the addition of a small amount of wood ash doubles the angle of repose of granulated ice at −10°C, but has little effect at −30°C. It is hypothesized that this behavior may be explained in terms of the freezing of water that is formed at the higher temperature through the reaction of salt within the ash.

Perennial ice can be studied for many purposes, including paleoclimate records or rheological properties. For most of those purposes, the ice microstructure must be studied, often through optical microscopy. The aim of this work is to assess the viability of immersion microscopy for the study of ice microstructures. It consists of using an oil between the objective lens and the specimen, to increase image resolution. Immersion microscopy is a technique well-developed for the investigation of diverse materials, but it has so far not been explored for ice research. Here we investigate the challenges and advantages of that technique. The main challenge is related to the selection of the immersion oil itself, which must satisfy a number of criteria, ranging from refractive index and viscosity to toxicity and reactivity. We identify pure silicone oil (dimethicone) as a simple and safe option for immersion microscopy of inner ice structures. Among its advantages, it provides higher resolution (compared to standard ‘dry’ microscopy) and it can be simultaneously used as a long-term coating to prevent undesired sublimation of the ice-sample surfaces. For the observation of surface structures, however, another type of oil with higher refractive index should be used.

To better constrain meltwater transport and ice viscosity in temperate glaciers, particularly in ice stream shear margins, we use a custom permeameter to study the untested model relationship between the permeability of temperate ice and its liquid water content. The permeability of lab-made ice of two mean grain diameters (1.8 and 4.2 mm) is measured, and water content is controlled with the ice salinity and measured calorimetrically. Fluorescein dye is added to through-flowing, chilled water to highlight flow pathways through the ice after experiments. As predicted by a simple model, permeability increases with approximately the square of the water content and by about three orders of magnitude across water contents of 0.1–4.4%. However, permeability values are less than those of the model by average factors of 2.6 and 4.1 for the finer and coarser ice, respectively. This discrepancy is likely due to tortuous, truncated or air-clogged veins. The order-of-magnitude agreement between measured and modeled values may indicate that reduced permeability from these factors is nearly compensated by preferential flow in oversized veins that are isolated or arborescent. Both kinds of preferred flow pathways are observed but the latter only in fine-grained ice at water contents > 2%.

We recently reported on the radio-frequency attenuation length of cold polar ice at Summit Station, Greenland, based on bi-static radar measurements of radio-frequency bedrock echo strengths taken during the summer of 2021. Those data also allow studies of (a) the relative contributions of coherent (such as discrete internal conducting layers with sub-centimeter transverse scale) vs incoherent (e.g. bulk volumetric) scattering, (b) the magnitude of internal layer reflection coefficients, (c) limits on signal propagation velocity asymmetries (‘birefringence’) and (d) limits on signal dispersion in-ice over a bandwidth of ~100 MHz. We find that (1) attenuation lengths approach 1 km in our band, (2) after averaging 10 000 echo triggers, reflected signals observable over the thermal floor (to depths of ~1500 m) are consistent with being entirely coherent, (3) internal layer reflectivities are ≈–60$\to$–70 dB, (4) birefringent effects for vertically propagating signals are smaller by an order of magnitude relative to South Pole and (5) within our experimental limits, glacial ice is non-dispersive over the frequency band relevant for neutrino detection experiments.

Small quantities of liquid water lining triple junctions in polycrystalline glacier ice form connected vein networks that enable material exchange with underlying basal environments. Diffuse debris concentrations commonly observed in ice marginal regions might be attributed to this mechanism. Following recent cryogenic ring-shear experiments, we observed emplacement along grain boundaries of loess particles several tens of microns in size. Here, we describe an idealized model of vein liquid flow to elucidate conditions favoring such particle transport. Gradients in liquid potential drive flow toward colder temperatures and lower solute concentrations, while deviations of the ice stress state from hydrostatic balance produce additional suction toward anomalously low ice pressures. Our model predicts particle entrainment following both modest warming along the basal interface resulting from anticipated natural changes in effective stress, and the interior relaxation of temperature and solute concentration imposed by our experimental protocols. Comparisons with experimental observations are encouraging, but suggest that liquid flow rates are somewhat higher and/or more effective at dragging larger particles than predicted by our idealized model with nominal parameter choices. Diffuse debris entrainment extending several meters above the glacier bed likely requires a more sophisticated treatment that incorporates effects of ice deformation or other processes.

Most theoretical descriptions of stresses induced by freezing are rooted in the (generalized) Clapeyron equation, which predicts the pressure that a solid can exert as it cools below its melting temperature. This equation is central for topics ranging beyond glaciology to geomorphology, civil engineering, food storage and cryopreservation. However, it has inherent limitations, requiring isotropic solid stresses and conditions near bulk equilibrium. Here, we examine when the Clapeyron equation is applicable by providing a rigorous derivation that details all assumptions. We demonstrate the natural extension for anisotropic stress states, and we show how the temperature and pressure ranges for validity depend on well-defined material properties. Finally, we demonstrate how the range of applicability of the (linear) Clapeyron equation can be extended by adding higher-order terms, yielding results that are in good agreement with experimental data for the pressure melting of ice.

Over the last 25 years, radiowave detection of neutrino-generated signals, using cold polar ice as the neutrino target, has emerged as perhaps the most promising technique for detection of extragalactic ultra-high energy neutrinos (corresponding to neutrino energies in excess of 0.01 Joules, or 1017 electron volts). During the summer of 2021 and in tandem with the initial deployment of the Radio Neutrino Observatory in Greenland (RNO-G), we conducted radioglaciological measurements at Summit Station, Greenland to refine our understanding of the ice target. We report the result of one such measurement, the radio-frequency electric field attenuation length $L_\alpha$. We find an approximately linear dependence of $L_\alpha$ on frequency with the best fit of the average field attenuation for the upper 1500 m of ice: $\langle L_\alpha \rangle = ( ( 1154 \pm 121) - ( 0.81 \pm 0.14) \, ( \nu /{\rm MHz}) ) \,{\rm m}$ for frequencies ν ∈ [145 − 350] MHz.

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.

The production of meltwater from glacier ice, which is exposed at the margins of land ice during the summer, is responsible for a large proportion of glacier mass loss. The rate of meltwater production from glacier ice is especially sensitive to its physical structure and chemical composition which combine to determine the albedo of glacier ice. However, the optical properties of near-surface glacier ice are not well known since most prior work has focused on laboratory-grown ice or deep cores. Here, we demonstrate a measurement technique based on diffuse propagation of nanosecond-duration laser pulses in near-surface glacier ice that enables the independent measurement of the scattering and absorption coefficients, allowing for a complete description of the processes governing radiative transfer. We employ a photon-counting detector to overcome the high losses associated with diffuse optics. The instrument is highly portable and rugged, making it optimally suited for deployment in remote regions. A set of measurements taken on Crook and Collier Glaciers, Oregon, serves as a demonstration of the technique. These measurements provide insight into both physical structure and composition of near-surface glacier ice and open new avenues for the analysis of light-absorbing impurities and remote sensing of the cryosphere.

Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

Sea-ice pore microstructure constrains ice transport properties, affecting fluid flow relevant to oil-in-ice transport and biogeochemical processes. Motivated by a lack of pore microstructural data, in particular for granular ice and across the seasonal cycle, throat size, tortuosity, connectivity, and other microstructural variables were derived from X-ray computed tomography for brine-filled pores in seasonal landfast ice off northern Alaska. Data were obtained for granular and columnar ice during the ice growth, transition, and melt season. While granular ice exhibits a more heterogeneous pore space than columnar ice, pore and throat size distributions are comparable. The greater tortuosity of pores in granular (1.2 < τg < 1.7) compared to columnar ice (1.0 < τc < 1.1) compounded with a less interconnected pore space translates into lower permeability for granular ice during the growth season for a given porosity. The microstructural data explain findings of granular ice hindering vertical oil-in-ice transport during ice growth and transition stage. With granular ice more frequent in the changing Arctic, data from studies such as this are needed to inform improved modeling of porosity-permeability relationships.

Results of ice-stream models that treat temperate ice deformation as a two-phase flow are sensitive to the ice permeability. We have constructed and begun using a custom, falling-head permeameter for measuring the permeability of temperate, polycrystalline ice. Chilled water is passed through an ice disk that is kept at the pressure-melting temperature while the rate of head decrease indicates the permeability. Fluorescein dye in the water allows water-vein geometry to be studied using fluorescence microscopy. Water flow over durations of seconds to hours is Darcian, and for grain diameter d increasing from 1.7 to 8.9 mm, average permeability decreases from 2 × 10−12 to 4 × 10−15 m2. In tests with dye on fine (d = 2 mm) and coarse (d = 7 mm) ice, average area-weighted vein radii are nearly equal, 41 and 34 μm, respectively. These average radii, if included in a theory slightly modified from Nye and Frank (1973), yield permeability values within a factor of 2.0 of best-fit values based on regression of the data. Permeability values depend on d−3.4, rather than d−2 as predicted by models if vein radii are considered independent of d. In future experiments, the dependence of permeability on liquid water content will be measured.

Ice flow models typically assume that ice is incompressible, a reasonable assumption because ice density changes are indeed small and have a correspondingly small effect on the overall mass balance of glaciers and ice sheets. Given the immense volume of the ice sheets, however, even relatively small changes may influence global mean sea level to a degree that severely impacts humanity (Hauer and others, 2020). Here, we quantify the role of gravitational compression and thermal contraction in estimating ice sheet mass.

Brash ice forms in harbours and ship channels from frequent ship passages and the resulting freezing–breaking cycles create a unique ice formation. The brash ice accumulation over the winter season is a result of meteorological, thermodynamical and mechanical processes. A reliable brash ice growth model is an important asset when determining navigation routes through ice conditions and when establishing port ice management solutions. This review aims to describe the brash ice development and its modelling as well as the key parameters that influence the brash ice growth and its estimation. This paper summarises the brash ice growth models and the fundamental theories of level ice growth upon which these models are based, and outlines the main knowledge gaps. The results highlight the importance of porosity and piece size distribution and their effect on the consolidation process. The inclusion of the brash ice lateral movement and the side ridge formation would improve the accuracy of forecast models. Furthermore, the findings of the study identify the effect of omitting meteorological parameters such as snow and radiation, from the brash ice growth models. Their contribution to the level ice thickness suggests a significant influence on the brash ice consolidation process.

We analysed the three-dimensional microstructure of sea ice by means of X-ray-micro computed tomography. Microscopic (brine- and air- pore sizes, numbers and connectivity) and macroscopic (salinity, density, porosity) properties of young Arctic sea ice were analysed. The analysis is based on ice cores obtained during spring 2016. Centrifuging of brine prior to CT imaging has allowed us to derive confident relationships between the open, vertically connected and total porosity of young sea ice at relatively high temperatures. We analysed the dependence of the microscopic properties on vertical position and total brine porosity. Most bulk properties (salinity, density) and pore space properties (pore sizes and their distribution) show a strong dependence on total brine porosity, but did not change significantly over the course of the field work. However, significant changes were observed for pore numbers (decreasing over time) and pore connectivity (increasing over time). CT-based salinity determinations are subject to larger than standard uncertainties (from conductivity), while the CT method yields important information about the salinity contributions from closed and open pores. We also performed a comparison of CT-based air porosity with calculations based on density from hydrostatic weighing. The consistency is encouraging and gives confidence to our CT-based results.

The importance of glacier sliding has motivated a rich literature describing the thermomechanical interactions between ice, liquid water and bed materials. Early recognition of the gradient in melting temperature across small bed obstacles led to focused studies of regelation. An appreciation for the limits on ice deformation rates downstream of larger obstacles highlighted a role for cavitation, which has subsequently gained prominence in descriptions of subglacial drainage. Here, we show that the changes in melting temperature that accompany changes in normal stress along a sliding ice interface near cavities and other macroscopic drainage elements cause appreciable supercooling and basal mass exchange. This provides the basis of a novel formation mechanism for widely observed laminated debris-rich basal ice layers.

Polycrystalline ice is known to exhibit macroscopic anisotropy in relative permittivity (ɛ) depending on the crystal orientation fabric (COF). Using a new system designed to measure the tensorial components of ɛ, we investigated the dielectric anisotropy (Δɛ) of a deep ice core sample obtained from Dome Fuji, East Antarctica. This technique permits the continuous nondestructive assessment of the COF in thick ice sections. Measurements of vertical prism sections along the core showed that the Δɛ values in the vertical direction increased with increasing depth, supporting previous findings of c-axis clustering around the vertical direction. Analyses of horizontal disk sections demonstrated that the magnitude of Δɛ in the horizontal plane was 10–15% of that in the vertical plane. In addition, the directions of the principal axes of tensorial ɛ in the horizontal plane corresponded to the long or short axis of the elliptically elongated single-pole maximum COF. The data confirmed that Δɛ in the vertical and horizontal planes adequately indicated the preferred orientations of the c-axes, and that Δɛ can be considered to represent a direct substitute for the normalized COF eigenvalues. This new method could be extremely useful as a means of investigating continuous and depth-dependent variations in COF.

In this study, we explore small-scale (~1 to 20 km) thermal-refractive effects on basal geothermal heat flux (BGHF) at subglacial boundaries resulting from lateral thermal conductivity contrasts associated with subglacial topography and geologic contacts. We construct a series of two-dimensional, conductive, steady-state models that exclude many of the complexities of ice sheets in order to demonstrate the effect of thermal refraction. We show that heat can preferentially flow into or around a subglacial valley depending on the thermal conductivity contrast with underlying bedrock, with anomalies of local BGHF at the ice–bedrock interface between 80 and 120% of regional BGHF and temperature anomalies on the order of ±15% for the typical range of bedrock conductivities. In the absence of bed topography, subglacial contacts can produce significant heat flux and temperature anomalies that are locally extensive (>10 km). Thermal refraction can result in either an increase or decrease in the likelihood of melting and ice-sheet stability depending on the conductivity contrast and bed topography. While our models exclude many of the physical complexities of ice behavior, they illustrate the need to include refractive effects created by realistic geology into future glacial models to improve the prediction of subglacial melting and ice viscosity.

Further to systematic experiments on the flexural strength of laboratory-grown, fresh water ice loaded cyclically, this paper describes results from new experiments of the same kind on lake ice harvested in Svalbard. The experiments were conducted at −12 °C, 0.1 Hz frequency and outer-fiber stress in the range from ~ 0.1 to ~ 0.7 MPa. The results suggest that the flexural strength increases linearly with stress amplitude, similar to the behavior of laboratory-grown ice.