Small mountain glaciers are an important part of the cryosphere and tend to respond rapidly to climate warming. Historically, mapping very small glaciers (generally considered to be <0.5 km2) using satellite imagery has often been subjective due to the difficulty in differentiating them from perennial snowpatches. For this reason, most scientists implement minimum size-thresholds (typically 0.01–0.05 km2). Here, we compare the ability of different remote-sensing approaches to identify and map very small glaciers on imagery of varying spatial resolutions (30–0.25 m) and investigate how operator subjectivity influences the results. Based on this analysis, we support the use of a minimum size-threshold of 0.01 km2 for imagery with coarse to medium spatial resolution (30–10 m). However, when mapping on high-resolution imagery (<1 m) with minimal seasonal snow cover, glaciers <0.05 km2 and even <0.01 km2 are readily identifiable and using a minimum threshold may be inappropriate. For these cases, we develop a set of criteria to enable the identification of very small glaciers and classify them as certain, probable or possible. This should facilitate a more consistent approach to identifying and mapping very small glaciers on high-resolution imagery, helping to produce more comprehensive and accurate glacier inventories.

The Martian mid-latitudes contain numerous small water-ice deposits, collectively termed viscous flow features (VFFs). The shape and topography of the deposits contain information on their past flow history and formation process. In order to access this information, it is imperative to get information on their deformational properties. Here we use a high-resolution digital topography map and ice-penetrating radar data in combination with an inverse method to constrain the deformational properties of a lobate debris apron, a class of VFF, in the southern hemisphere of Mars. We find that while the creep parameter and accumulation rates are not well constrained in absolute values, their ratio is robust. We also find that the creep exponent is most likely n ≤ 3.

Supraglacial streams are a significant part of the glacial hydrological system and important for understanding glacial hydrology and dynamics. Here we infer factors that influence the long-term development of perennial supraglacial streams, particularly in reference to canyon, incised and surface stream formation. Orthophotos and digital elevation models generated from high-resolution aerial imagery taken with unmanned aerial vehicles or piloted helicopters between 2010 and 2017 were used to compare seven streams on Fountain Glacier, Bylot Island, Canada over time. Results show canyon formation occurs from a combination of fluvial processes and the impact of solar radiation. The greater the discharge or slope, the faster the incision and higher the sinuosity. With greater sinuosity, the distance between the top of the valley banks increases, and cutoffs cause trapezoidal canyon-like valleys to form. Solar radiation causes the backward migration of the valley walls, further enhancing canyon area. Canyons are less likely to occur in areas of low discharge and slope. Less incised channels are also more likely to have water flow jumping the channel banks, changing the channel path. The presence of medial moraines and crevasses also increases rerouting of small streams. Lastly, windblown created snow-plugs may lead to stream diversion.

The Ross Ice Shelf (RIS) is host to a broadband, multimode seismic wavefield that is excited in response to atmospheric, oceanic and solid Earth source processes. A 34-station broadband seismographic network installed on the RIS from late 2014 through early 2017 produced continuous vibrational observations of Earth's largest ice shelf at both floating and grounded locations. We characterize temporal and spatial variations in broadband ambient wavefield power, with a focus on period bands associated with primary (10–20 s) and secondary (5–10 s) microseism signals, and an oceanic source process near the ice front (0.4–4.0 s). Horizontal component signals on floating stations overwhelmingly reflect oceanic excitations year-round due to near-complete isolation from solid Earth shear waves. The spectrum at all periods is shown to be strongly modulated by the concentration of sea ice near the ice shelf front. Contiguous and extensive sea ice damps ocean wave coupling sufficiently so that wintertime background levels can approach or surpass those of land-sited stations in Antarctica.

Spectral reflectance of natural snow samples representing various stratigraphies was investigated in a controlled dark laboratory environment. Mean and Std dev. of band specific reflectance values were determined for several satellite sensor bands utilized in remote sensing of snow. The reflectance values for dry, moist, wet and wet and littered snow for different instruments varied between 0.63–0.97 in the visible and near-infrared bands at an incoming light zenith angle of θ = 55°. The results indicate that in MODIS band 4 (545–565 nm), essential to snow mapping, the reflectance of snow drops by 9% when dry snow changes to wet snow and by a further 10% when typical forest litter inclusions exist on the wet snow surface. A separate investigation of individual snow types revealed that they can be grouped either as dry or wet snow based on their spectral behavior. However, some snow types were located between these two distinct groups, such as snow with near-surface melt-freeze crusts, and could not be clearly distinguished. The reflectance statistics collected and analyzed here can be directly used to refine accuracy characterization and parametrization of snow mapping algorithms, such as the SCAmod method, used for the mapping of snow cover area.

Light-absorbing impurities (LAIs, e.g. black carbon (BC), organic carbon (OC), mineral dust (MD)) deposited on snow cover reduce albedo and accelerate its melting. Northern Xinjiang (NX) is an arid and semi-arid inland region, where snowmelt leads to frequent floods that have been a serious threat to local ecological security. There is still a lack of quantitative assessments of the effects of LAIs on snowmelt in the region. This study investigates spatial variations of LAIs in snow and its effect on snow albedo, radiative forcing (RF) and snowmelt across NX. Results showed that concentrations of BC, OC (only water-insoluble OC), MD ranged from 32 to 8841 ng g−1, 77 to 8568 ng g−1 and 0.46 to 236 µg g−1, respectively. Weather Research and Forecasting Chemistry model suggested that residential emission was the largest source of BC. Snow, Ice, and Aerosol Radiative modelling showed that the average contribution of BC and MD to snow albedo reduction was 17 and 3%, respectively. RF caused by BC significantly exceeded RF caused by MD. In different scenarios, changes in snow cover duration (SCD) caused by BC and MD decreased by 1.36 ± 0.61 to 6.12 ± 3.38 d. Compared with MD, BC was the main dominant factor in reducing snow albedo and SCD across NX.

Wind-driven snow redistribution can increase the spatial heterogeneity of snow accumulation on ice caps and ice sheets, and may prove crucial for the initiation and survival of glaciers in areas of marginal glaciation. We present a snowdrift model (Snow_Blow), which extends and improves the model of Purves, Mackaness and Sugden (1999, Journal of Quaternary Science 14, 313–321). The model calculates spatial variations in relative snow accumulation that result from variations in topography, using a digital elevation model (DEM) and wind direction as inputs. Improvements include snow redistribution using a flux routing algorithm, DEM resolution independence and the addition of a slope curvature component. This paper tests Snow_Blow in Antarctica (a modern environment) and reveals its potential for application in palaeoenvironmental settings, where input meteorological data are unavailable and difficult to estimate. Specifically, Snow_Blow is applied to the Ellsworth Mountains in West Antarctica where ablation is considered to be predominantly related to wind erosion processes. We find that Snow_Blow is able to replicate well the existing distribution of accumulating snow and snow erosion as recorded in and around Blue Ice Areas. Lastly, a variety of model parameters are tested, including depositional distance and erosion vs wind speed, to provide the most likely input parameters for palaeoenvironmental reconstructions.

Finely resolved geodetic data provide an opportunity to assess the extent and morphology of crevasses and their change over time. Crevasses have the potential to bias geodetic measurements of elevation and mass change unless they are properly accounted for. We developed a framework that automatically maps and extracts crevasse geometry and masks them where they interfere with surface mass-balance assessment. Our study examines airborne light detection and ranging digital elevation models (LiDAR DEMs) from Haig Glacier, which is experiencing a transient response in its crevassed upper regions as the glacier thins, using a self-organizing map algorithm. This method successfully extracts and characterizes ~1000 crevasses, with an overall accuracy of 94%. The resulting map provides insight into stress and flow conditions. The crevasse mask also enables refined geodetic estimates of summer mass balance. From differencing of September and April LiDAR DEMs, the raw LiDAR DEM gives a 9% overestimate in the magnitude of glacier thinning over the summer: −5.48 m compared with a mean elevation change of −5.02 m when crevasses are masked out. Without identification and removal of crevasses, the LiDAR-derived summer mass balance therefore has a negative bias relative to the glaciological surface mass balance.

Englacial layers in Antarctica and Greenland are indicators of the dynamic, rheological and subglacial configuration of the ice sheets. Airborne radar sounder data is the primary remote sensing solution for directly observing englacial layers and structures at the glacier-catchment to ice-sheet scale. However, when traditional along-track synthetic aperture radar (SAR) processing is applied, steep layers can disappear, limiting the detectability and interpretability of englacial layer geometry. This study provides a reconstruction algorithm to address the problem of destructive phase interference during the radargram formation. We develop and apply a novel SAR processor optimized for layer detection that enhances the Signal-to-Noise ratio (SNR) of specular reflectors. The algorithm also enables the automatic estimation of layer slope. We demonstrate the algorithm using data acquired at the Institute Ice Stream, West Antarctica.

While large-scale observations of intensified fracture and rifting can be observed through remote-sensing observations, understanding crevasse initiation may best be achieved with small-scale observations in which crevasses can be directly observed. Here we investigate the kinematic drivers of crevasse initiation in the McMurdo Shear Zone (MSZ), Antarctica. We delineated 420 crevasses from ~95 km of 400 MHz frequency ground-penetrating radar data and compared these data with kinematic outputs derived from remotely-sensed ice surface velocities to develop a statistical method to estimate crevasse initiation threshold strain rate values. We found the MSZ to be dominated by simple shear and that surface shear strain rates proved best for predicting crevasse features, with regions of higher shear strain rate more likely to have a greater number of crevasses. In the surveyed portion of our study region, values of shear strain rate and vorticity rate derived from the MEaSUREs2 velocity dataset range between 0.005–0.020 and 0.006–0.022 a−1, respectively, with crevasses located at ≥0.011 and ≥0.013 a−1. While threshold values from this study cannot be directly applied to other glacial environments, the method described here should allow for the study of shear margin evolution and assessment of localized damage and weakening processes in other locations where in situ data are available.

Modelling paleo-glacier networks in mountain ranges on the millennial timescales requires ice flow approximations. Hybrid models calculating ice flow by combining vertical shearing (shallow ice approximation) and longitudinal stretching (shallow shelf approximation) have been applied to model paleo-glacier networks on steep terrain, yet their validity has not yet been assessed quantitatively. Moreover, hybrid models consistently yield higher ice thicknesses than Last Glacial Maximum geomorphological reconstructions in the European Alps. Here, we compare results based on the hybrid Parallel Ice Sheet Model (PISM) and the Stokes model Elmer/Ice on the Rhine Glacier, a catchment of the former European Alpine Icefield. For PISM, we also test two magnitudes of flux limitation in a scheme that reduces shearing velocities. We find that the flux limitation typically used in PISM yields significantly reduced shearing speeds and increases ice thicknesses by up to 500 m, partly explaining previous overestimations. However, reducing the ice flux limitation allows the hybrid model to minimize this mismatch and captures sliding speeds, ice thicknesses, ice extent and basal temperatures in close agreement with those obtained with the Stokes model.

A new type of ice core drill bit, designed with a vane swirler, was developed for ice core drilling with air reverse circulation. An orthogonal experimental design method was employed to investigate the effects of the swirler structure parameters on the reverse circulation performance of the drill bit including helical angle, number of blades, blade length and blade central angle, etc. The entrainment ratio was used to evaluate the reverse circulation effectiveness of the drill bit. The results show that the helical angle is the dominant factor regardless of whether or not the flushing nozzles are part of the design of the drill bit. The number of blades is the least important factor for the drill bit designed with the flushing nozzles (referred to as drill bit I), while the outlet area of the swirling slot is the least influential factor for the drill bit without flushing nozzles (referred to as drill bit П). In addition, the appearance of the ice core has a certain effect on the air reverse circulation for both drill bits. Within the ranges of this study, the optimal structure of the drill bit was determined based on the range analysis of the orthogonal design.

The spatial distribution of basal water critically impacts the evolution of ice sheets. Current estimates of basal water distribution beneath the Greenland Ice Sheet (GrIS) contain large uncertainties due to poorly constrained boundary conditions, primarily from geothermal heat flux (GHF). The existing GHF models often contradict each other and implementing them in numerical ice-sheet models cannot reproduce the measured temperatures at ice core locations. Here we utilize two datasets of radar-detected basal water in Greenland to constrain the GHF at regions with a thawed bed. Using the three-dimensional ice-sheet model SICOPOLIS, we iteratively adjust the GHF to find the minimum GHF required to reach the bed to the pressure melting point, GHFpmp, at locations of radar-detected basal water. We identify parts of the central-east, south and northwest Greenland with significantly high GHFpmp. Conversely, we find that the majority of low-elevation regions of west Greenland and parts of northeast have very low GHFpmp. We compare the estimated constraints with the available GHF models for Greenland and show that GHF models often do not honor the estimated constraints. Our results highlight the need for community effort to reconcile the discrepancies between radar data, GHF models, and ice core information.

Predicting the isotopic modification of ice by melting processes is important for improving the accuracy in paleoclimate reconstruction. To this end, we present results from cold room laboratory observations of changes in the isotopic ratio (D/H and 18O/16O) of ice cubes by isotopic exchange between liquid water and ice in nearly isothermal conditions. A 1-D model was fit to the isotopic results by adjusting the values of two parameters, the isotopic exchange rate constant (kr) and the fraction of ice participating in the exchange (f). We found that the rate constant for hydrogen isotopic exchange between liquid water and ice may be greater (up to 40%) than that for the oxygen isotopic exchange. The range of the rate constant obtained from four melt experiments is from 0.21 to 0.82 h–1. The model results also suggest that f decreases with the increasing wetness of the ice. This is because with increasing water saturation in ice, water may be present only in the small pores or some of the water that was exchanged with ice may be bypassed, decreasing the effective surface area over which the isotopic exchange can occur. The relationship between the two water isotopes (δ18O vs δD) was observed and modeled and the slope was <8, which is significantly different from the slope of the meteoric waterline. We note that these slopes were obtained without considering the sublimation process.