Subglacial hydrologic systems in Antarctica and Greenland play a fundamental role in ice-sheet dynamics, yet critical aspects of these systems remain poorly understood due to a lack of observations. Ground-based electromagnetic (EM) geophysical methods are established for mapping groundwater in many environments, but have never been applied to imaging lakes beneath ice sheets. Here, we study the feasibility of passive- and active-source EM imaging for quantifying the nature of subglacial water systems beneath ice streams, with an emphasis on the interfaces between ice and basal meltwater, as well as deeper groundwater in the underlying sediments. We describe a suite of model studies that exam the data sensitivity as a function of ice thickness, water conductivity and hydrologic system geometry for models representative of a subglacial lake and a grounding zone estuary. We show that EM data are directly sensitive to groundwater and can image its lateral and depth extent. By combining the conductivity obtained from EM data with ice thickness and geological structure from conventional geophysical techniques, such as ground-penetrating radar and active seismic surveying, EM data have the potential to provide new insights on the interaction between ice, rock and water at critical ice-sheet boundaries.

Grounding zones are vital to ice-sheet mass balance and its coupling to the global ocean circulation. Processes here determine the mass discharge from the grounded ice sheet, to the floating ice shelves. The response of this transition zone to tidal forcing has been described by both elastic and viscoelastic models. Here we examine the validity of these models for grounding zone flexure over tidal timescales using field data from the Southern McMurdo Ice Shelf (78° 15′S, 167° 7′E). Observations of tidal movement were carried out by simultaneous tiltmeter and GPS measurements along a profile across the grounding zone. Finite-element simulations covering a 64 d period reveal that the viscoelastic model fits best the observations using a Young's modulus of 1.6 GPa and a viscosity of 1013.7 Pa s (≈ 50.1 TPa s). We conclude that the elastic model is only well-constrained for tidal displacements >35% of the spring-tidal amplitude using a Young's modulus of 1.62 ± 0.69 GPa, but that a viscoelastic model is necessary to adequately capture tidal bending at amplitudes below this threshold. In grounding zones where bending stresses are greater than at the Southern McMurdo Ice Shelf or ice viscosity is lower, the threshold would be even higher.

Fine resolution topographic data derived from methods such as Structure from Motion (SfM) and Multi-View Stereo (MVS) have the potential to provide detailed observations of geomorphological change, but have thus far been limited by the logistical constraints of conducting repeat surveys in the field. Here, we present the results from an automated time-lapse camera array, deployed around an ice-marginal lake on the western margin of the Greenland ice sheet. Fifteen cameras acquired imagery three-times per day over a 426 day period, yielding a dataset of ~19 000 images. From these data we derived 18 point clouds of the ice-margin across a range of seasons and successfully identified calving events (ranging from 234 to 1475 m2 in area and 815–8725 m3 in volume) induced by ice cliff undercutting at the waterline and the collapse of spalling flakes. Low ambient light levels, locally reflective surfaces and the large survey range hindered analysis of smaller scale ice-margin dynamics. Nevertheless, this study demonstrates that an integrated SfM-MVS and time-lapse approach can be employed to generate long-term 3-D topographic datasets and thus quantify ice-margin dynamics at a fine spatio-temporal scale. This approach provides a template for future studies of geomorphological change.

Surface meltwater can refreeze within firn layers and crevasses to warm ice through latent-heat transfer on decadal to millennial timescales. Earlier work posited that the consequent softening of the ice might accelerate ice flow, potentially increasing ice-sheet mass loss. Here, we calculate the effect of meltwater refreezing on ice temperature and softness in the Pâkitsoq (near Swiss Camp) and Jakobshavn Isbræ regions of western Greenland using a numeric model and existing borehole measurements. We show that in the Jakobshavn catchment, meltwater percolation within the firn warms the ice at depth by 3–5°C. By contrast, meltwater refreezing in crevasses (cryo-hydrologic warming) at depths of ~300 m warms the ice in Pâkitsoq by up to 10°C, but this causes minimal increase in ice motion (<10 m a−1). Pâkitsoq is representative of western Greenland's land-terminating ice, where the slow movement of ice through a wide ablation zone provides ideal conditions for cryo-hydrologic warming to occur. We find that only ~37% of the western Greenland ice flux, however, travels through such areas. Overall, our findings suggest that cryo-hydrologic warming will likely have only a limited effect on the dynamic evolution of the Greenland ice sheet.

Taylor Glacier hosts an active englacial hydrologic system that feeds Blood Falls, a supraglacial outflow of iron-rich subglacial brine at the terminus, despite mean annual air temperatures of −17°C and limited surface melt. Taylor Glacier is an outlet glacier of the East Antarctic ice sheet that terminates in Lake Bonney, McMurdo Dry Valleys. To image and map the brine feeding Blood Falls, we used radio echo sounding to delineate a subhorizontal zone of englacial brine upstream from Blood Falls and elongated in the ice flow direction. We estimate volumetric brine content in excess of 13% within 2 m of the central axis of this zone, and likely much higher at its center. Brine content decreases, but remains detectable, up to 45 m away along some transects. Hence, we infer a network of subparallel basal crevasses allowing injection of pressurized subglacial brine into the ice. Subglacial brine is routed towards Blood Falls by hydraulic potential gradients associated with deeply incised supraglacial valleys. The brine remains liquid within the subglacial and englacial environments through latent heat of freezing coupled with elevated salt content. Our findings suggest that cold glaciers could support freshwater hydrologic systems through localized warming by latent heat alone.

Recognising the scarcity of glacier mass-balance data in the Southern Hemisphere, a mass-balance measurement programme was started at Brewster Glacier in the Southern Alps of New Zealand in 2004. Evolution of the measurement regime over the 11 years of data recorded means there are differences in the spatial density of data obtained. To ensure the temporal integrity of the dataset a new geostatistical approach is developed to calculate mass balance. Spatial co-variance between elevation and snow depth allows a digital elevation model to be used in a co-kriging approach to develop a snow depth index (SDI). By capturing the observed spatial variability in snow depth, the SDI is a more reliable predictor than elevation and is used to adjust each year of measurements consistently despite variability in sampling spatial density. The SDI also resolves the spatial structure of summer balance better than elevation. Co-kriging is used again to spatially interpolate a derived mean summer balance index using SDI as a co-variate, which yields a spatial predictor for summer balance. The average glacier-wide surface winter, summer and annual balances over the period 2005–15 are 2484, −2586 and −102 mm w.e., respectively, with changes in summer balance explaining most of the variability in annual balance.

Retreating glaciers are a consequence of a warming climate. Thus, numerous monitoring campaigns are being carried out to increase understanding of this on-going process. One phenomenon related to dynamic glacial changes is glacier-induced seismicity; however, weak seismic events are difficult to record due to the sparse seismological network in arctic areas. We have developed an automatic procedure capable of detecting glacier-induced seismic events using records from a single permanent seismological station. To distinguish between glacial and non-glacial signals, we developed a fuzzy logic algorithm based on the signal frequency and energy flow analysis. We studied the long-term changes in glacier-induced seismicity in Hornsund (southern Spitsbergen) and in Kongsfjorden (western Spitsbergen). We found that the number of detected glacial-origin events in the Hornsund dataset over the years 2013-14 has doubled. In the Kongsfjorden dataset, we observed a steady increase in the number of glacier-induced events with each year. We also observed that the seasonal event distribution correlates best with 1 month lagged temperatures, and that extreme rain events can intensify seismic emissions. Our study demonstrates the possibility of using long-term seismological observations from a single permanent station to automatically monitor the dynamic activity of nearby glaciers and retrieve its characteristic features.

We present climate data, direct surface mass balance (SMB) observations and model results for Mocho Glacier in the Chilean Lake District. Mean annual temperature on a nunatak of Mocho Glacier at an elevation of ~2000 m was +2.6°C in 2006–15 and mean annual precipitation in Puerto Fuy (13 km from the glacier, at an elevation of 600 m) was 4000 mm for the same period. High interannual variations in the SMB of Mocho Glacier were observed. A simple SMB model is able to reproduce the observed annual variations in SMB, but fails to predict the steep observed mass-balance gradient. The average of the measured annual glacier mass balances in the four hydrological years 2009/10–2012/13 was −0.90 m w.e. a−1 and the average modelled annual glacier mass balance 2006/07–2014/15 was −1.05 m w.e. a−1. The observed distributed ablation shows a clear altitudinal dependency, while accumulation is determined by patterns of snow drift as well. These patterns are only poorly represented in the model and have to be included in order to be able to reproduce a realistic SMB map of the glacier.

The extensive dataset of Variegated Glacier's 1982/83 surge and antecedent quiescent phase is fragmentary and predominantly of flow line nature. Applying the raw centre line data in a flow line model to conduct research into the mechanism behind Variegated Glacier's surge behaviour inevitably entails problems. In this study, the incomplete dataset is extrapolated into the upper and lower parts of the glacier. Furthermore, the centre line data are adapted to account for differences between width-averaged and centre line surface mass balance and ice thickness change, inflow from tributaries, and changes in surface width over time. The modifications to the dataset are backed-up by observations and clear, plausible physical explanations. Moreover, the revised dataset meets the imposed constraints regarding ice volume flux, specific mass balance and net volume change. Hence, the final dataset is considered a satisfying revision that makes the dataset more valuable for future research. Subsequent application of the revised dataset corroborates the idea that glacier evolution during quiescence is basically the growth back towards steady state after the glacier was brought out of balance by the preceding surge.

A new instrument for high-resolution optical logging has been built and tested in Antarctica. Its purpose is to obtain records of volcanic products and other scattering features, such as bubbles and impurities, preserved in polar ice sheets, and it achieves this by using long wavelength near-infrared light that is absorbed by the ice before many scattering events occur. Longer wavelengths ensure that the return signal is composed primarily of a single or few backscattering event(s) that limit its spatial spread. The compact optical logger features no components on its body that draw power, which minimizes its size and weight. A prototype of the logger was built and tested at Siple Dome A borehole, and the results were correlated with prior optical logging profiles and records of volcanic products from collected ice core samples.

Evolution of cold dry snow and firn plays important roles in glaciology; however, the physical formulation of a densification law is still an active research topic. We forced eight firn-densification models and one seasonal-snow model in six different experiments by imposing step changes in temperature and accumulation-rate boundary conditions; all of the boundary conditions were chosen to simulate firn densification in cold, dry environments. While the intended application of the participating models varies, they are describing the same physical system and should in principle yield the same solutions. The firn models all produce plausible depth-density profiles, but the model outputs in both steady state and transient modes differ for quantities that are of interest in ice core and altimetry research. These differences demonstrate that firn-densification models are incorrectly or incompletely representing physical processes. We quantitatively characterize the differences among the results from the various models. For example, we find depth-integrated porosity is unlikely to be inferred with confidence from a firn model to better than 2 m in steady state at a specific site with known accumulation rate and temperature. Firn Model Intercomparison Experiment can provide a benchmark of results for future models, provide a basis to quantify model uncertainties and guide future directions of firn-densification modeling.

Wind-packed snow in the form of slabs or crusts is an important part of alpine and polar snow covers. Yet, the formation process of such layers is poorly understood. For example, it remains unclear whether drifting snow is necessary for wind-packing or not. A better understanding of wind-packing could improve snow-cover models and avalanche danger forecasts and contribute to the assessment of mass balances in polar regions. We designed a closed-circuit, obround wind tunnel to study the process of wind crust formation. A SnowMicroPen was used to measure how the hardness of the snow evolved. The results show that no crust forms without saltation. Drifting snow is a necessary but not sufficient condition for wind-packing. The dynamics of erosion and deposition appear to be equally important.

Tidal flexure in ice shelf grounding zones has been used extensively in the past to determine grounding line position and ice properties. Although the rheology of ice is viscoelastic at tidal loading frequencies, most modelling studies have assumed some form of linear elastic beam approximation to match observed flexure profiles. Here we use density, radar and DInSAR measurements in combination with full-Stokes viscoelastic modelling to investigate a range of additional controls on the flexure of the Southern McMurdo Ice Shelf. We find that inclusion of observed basal crevasses and density dependent ice stiffness can greatly alter the flexure profile and yet fitting a simple elastic beam model to that profile will still produce an excellent fit. Estimates of the effective Young's modulus derived by fitting flexure profiles are shown to vary by over 200% depending on whether these factors are included, even when the local thickness is well constrained. Conversely, estimates of the grounding line position are relatively insensitive to these considerations for the case of a steep bed slope in our study region. By fitting tidal amplitudes only, and ignoring phase information, elastic beam theory can provide a good fit to observations in a wide variety of situations. This should, however, not be taken as an indication that the underlying rheological assumptions are correct.

Microtremor measurements and the horizontal-to-vertical spectral ratio (HVSR) technique, generally used for site effect studies as well as to determine the thickness of soft sedimentary layers, can effectively be applied to map the thickness of glaciers. In this work the radio-echo sounding, geoelectric and active seismic methods, widely employed to image the earth interior, are applied to verify the reliability of the HVSR technique in Alpine and Antarctic glacial environments. The technique has been used to analyze passive seismic data from glaciers of the Adamello and Ortles-Cevedale massifs (Italy), the Bernese Oberland Alps (Switzerland) and from the Whillans Ice Stream (West Antarctica). Comparing with the results obtained from the different geophysical imaging methods, we show that the resonance frequency in the HVSR spectra correlates well with the ice thickness at the site, in a wide range from a few tens of meters to more than 800 m. The reliability of the method mainly depends on the coupling of sensors at the glacier surface and on the basal impedance contrast. This passive seismic technique offers a logistically efficient and cost effective method to map glacier and ice-sheet thicknesses. Moreover, under certain conditions, it allows reliable estimations of the basal seismic properties.

We investigate the temporal evolution and spatial distribution of mass balance on the glacier Midtre Lovénbreen, Svalbard. Running a diagnostic high-resolution full-stress ice flow model with geometries obtained from five digital elevation models (DEMs) in the period 1962–2005, we compute velocity fields and linearly interpolated volume change of the glacier. We evaluate the kinematic free surface equation using these model outputs to solve the surface mass balance (SMB). Monitoring data on Midtre Lovénbreen allows model results to be compared with point measurements from the glacier over several decades. This method allows us to estimate the mass balance over the entire glacier surface, beyond the spatially limited field measurements, and to derive past SMB over an extended time period.

Hydraulic roughness exerts an important but poorly understood control on water pressure in subglacial conduits. Where relative roughness values are <5%, hydraulic roughness can be related to relative roughness using empirically-derived equations such as the Colebrook–White equation. General relationships between hydraulic roughness and relative roughness do not exist for relative roughness >5%. Here we report the first quantitative assessment of roughness heights and hydraulic diameters in a subglacial conduit. We measured roughness heights in a 125 m long section of a subglacial conduit using structure-from-motion to produce a digital surface model, and hand-measurements of the b-axis of rocks. We found roughness heights from 0.07 to 0.22 m and cross-sectional areas of 1–2 m2, resulting in relative roughness of 3–12% and >5% for most locations. A simple geometric model of varying conduit diameter shows that when the conduit is small relative roughness is >30% and has large variability. Our results suggest that parameterizations of conduit hydraulic roughness in subglacial hydrological models will remain challenging until hydraulic diameters exceed roughness heights by a factor of 20, or the conduit radius is >1 m for the roughness elements observed here.

The full history of ice sheet and climate interactions is recorded in the vertical profiles of geochemical tracers in polar ice sheets. Numerical simulations of these archives promise great advances both in the interpretation of these reconstructions and the validation of the models themselves. However, fundamental mathematical shortcomings of existing models subject tracers to spurious diffusion, thwarting straightforward solutions. Here, I propose a new vertical discretization for ice-sheet models that eliminates numerical diffusion entirely. Vertical motion through the model mesh is avoided by mimicking the real-world flow of ice as a thinning of underlying layers. A new layer is added to the surface at equidistant time intervals, isochronally, thus identifying each layer uniquely by its time of deposition and age. This new approach is implemented for a two-dimensional section through the summit of the Greenland ice sheet. The ability to directly compare simulations of vertical ice cores with reconstructed data is used to find optimal model parameters from a large ensemble of simulations. It is shown that because this tuning method uses information from all times included in the ice core, it constrains ice-sheet sensitivity more robustly than a realistic reproduction of the modern ice-sheet surface.

The Darwin–Hatherton Glacial system (DHGS) connects the East Antarctic Ice Sheet (EAIS) with the Ross Ice Shelf and is a key area for understanding past variations in ice thickness of surrounding ice masses. Here we present the first detailed measurements of ice thickness and grounding zone characteristics of the DHGS as well as new measurements of ice velocity. The results illustrate the changes that occur in glacier geometry and ice flux as ice flows from the polar plateau and into the Ross Ice Shelf. The ice discharge and the mean basal ice shelf melt for the first 8.5 km downstream of the grounding line amount to 0.24 ± 0.05 km3 a−1 and 0.3 ± 0.1 m a−1, respectively. As the ice begins to float, ice thickness decreases rapidly and basal terraces develop. Constructed maps of glacier geometry suggest that ice drainage from the EAIS into the Darwin Glacier occurs primarily through a deep subglacial canyon. By contrast, ice thins to <200 m at the head of the much slower flowing Hatherton Glacier. The glaciological field study establishes an improved basis for the interpretation of glacial drift sheets at the link between the EAIS and the Ross Ice Sheet.

We describe methods for measuring crystal orientation fabric with sonic waves in an ice core borehole, with special attention paid to vertical-girdle fabrics that are prevalent at the WAIS Divide. The speed of vertically propagating compressional waves in ice is influenced by vertical clustering of the ice crystal c-axes. Shear-wave speeds – particularly the speed separation between fast and slow shear polarizations – are sensitive to azimuthal anisotropy. Sonic data from the WAIS Divide complement thin-section measurements of fabric. Thin sections show a steady transition to strong girdle fabrics in the upper 2000 m of ice, followed by a transition to vertical-pole fabrics below 2500 m depth. Compressional-wave sonic data are inconclusive in the upper ice, due to noise, as well as the method's inherent insensitivity to girdle fabrics. Compared with available thin sections, sonic data provide better resolution of the transition to pole fabrics below 2500 m, notably including an abrupt increase in vertical clustering near 3000 m. Our compressional-wave measurements resolve fabric changes occurring over depth ranges of a few meters that cannot be inferred from available thin sections, but are sensitive only to zenithal anisotropy. Future logging tools should be designed to measure shear waves in addition to compressional waves, especially for logging in regions where ice flow patterns favor the development of girdle fabrics.

We describe the use of a terrestrial laser scanner (TLS) to monitor the net mass balance of Urumqi Glacier No. 1, eastern Tien Shan. We used an ultra-long-range Riegl VZ®-6000 TLS, which is specially designed for surveying snow- and ice-covered terrain, to create repeated high spatiotemporal resolution DEMs, focusing on the monthly-scale (25 April–28 May 2015) net mass balance. According to the TLS-derived DEMs, the area of Urumqi Glacier No. 1 was 1.558 km2 on 25 April 2015 and the average surface elevation change was 0.225 m. By comparing the results from the use of TLS with the conventional glaciological mass-balance method, the correlation coefficient (R2) between glaciological elevation changes of individual stakes and the TLS-derived geodetic elevation change of corresponding points was 0.85. Considering the uncertainty of both methods, this is a promising result. Using the in situ measured snow densities (snow pits) of the glacier surface, the geodetic net mass balance was 0.074 m w.e., which is slightly positive. The mean uncertainty in the TLS-derived monthly net mass balance was 0.018 m w.e., showing that the TLS surveying system presented accurate and relevant results and is therefore suitable to monitor mass-balance evolution of mountain glaciers.