The interpretations of relevant interfaces (i.e. the surface and bed) in radar sounding datasets over glaciers and ice sheets are primary boundary conditions in a variety of climate studies and particularly subglacial water routing models. It is therefore necessary to ensure these interpretations are consistent and not affected by cross-track clutter. For the surface interface, interferometry and a family of methods relying on digital elevation models have been used to successfully discriminate cross-track surface clutter. Here we present how interferometry can be applied to the problem of basal clutter from cross-track bed topography. Our approach is based on a comparison of the differential phases of ambiguous reflectors that may represent bed clutter and the differential phase of a reflector in an adjacent area that appears unaffected by basal clutter. The reflector yielding the smallest interferometric phase difference relative to the unambiguous bed reflector is considered to represent its consistent continuation. We successfully demonstrate our approach using 60 MHz center frequency MARFA data collected over Devon Ice Cap in the Canadian Arctic. Finally, we investigate the effects of clutter-affected and interferometry-corrected bed interpretations on ice layer thickness estimates, basal hydraulic head gradients and the potential extent of inferred subglacial water bodies.

Sea ice is composed of columnar-shaped grains. To investigate the influence of the loading direction on the uniaxial compressive strength and failure processes of sea ice, field experiments were performed with first-year level ice. Loads were applied both horizontally (parallel to the grain columns) and vertically (across the grain columns) with various nominal strain rates. Two failure modes have been observed: a ductile failure mode at low nominal strain rates, and a brittle failure mode at high nominal strain rates. However, the failure pattern of sea ice was clearly dependent on the loading direction. At low nominal strain rates (ductile failure mode), the sea-ice samples yielded due to the development of wing cracks under horizontal loading and due to splaying out at one end under vertical loading. When sea ice fails in the ductile mode, the deformation is driven by grain boundary sliding under horizontal loading and by grain decohesion and crystal deflection under vertical loading. At high nominal strain rates (brittle failure mode), the sea-ice samples failed in shear faulting under horizontal loading and in cross-column buckling under vertical loading. The nominal strain rate at the brittle–ductile transition zone is about ten times higher under vertical loading.

Basal melt of ice shelves is not only an important part of Antarctica's ice sheet mass budget, but it is also the origin of platelet ice, one of the most distinctive types of sea ice. In many coastal Antarctic regions, ice crystals form and grow in supercooled plumes of Ice Shelf Water. They usually rise towards the surface, becoming trapped under an ice shelf as marine ice or forming a semi-consolidated layer, known as the sub-ice platelet layer, below an overlying sea ice cover. In the latter, sea ice growth consolidates loose crystals to form incorporated platelet ice. These phenomena have numerous and profound impacts on the physical properties, biological processes and biogeochemical cycles associated with Antarctic fast ice: platelet ice contributes to sea ice mass balance and may indicate the extent of ice-shelf basal melting. It can also host a highly productive and uniquely adapted ecosystem. This paper clarifies the terminology and reviews platelet ice formation, observational methods as well as the geographical and seasonal occurrence of this ice type. The physical properties and ecological implications are presented in a way understandable for physicists and biologists alike, thereby providing the background for much needed interdisciplinary research on this topic.

Here we use polarimetric measurements from an Autonomous phase-sensitive Radio-Echo Sounder (ApRES) to investigate ice fabric within Whillans Ice Stream, West Antarctica. The survey traverse is bounded at one end by the suture zone with the Mercer Ice Stream and at the other end by a basal ‘sticky spot’. Our data analysis employs a phase-based polarimetric coherence method to estimate horizontal ice fabric properties: the fabric orientation and the magnitude of the horizontal fabric asymmetry. We infer an azimuthal rotation in the prevailing horizontal c-axis between the near-surface (z ≈ 10–50 m) and deeper ice (z ≈ 170–360 m), with the near-surface orientated closer to perpendicular to flow and deeper ice closer to parallel. In the near-surface, the fabric asymmetry increases toward the center of Whillans Ice Stream which is consistent with the surface compression direction. By contrast, the fabric orientation in deeper ice is not aligned with the surface compression direction but is consistent with englacial ice reacting to longitudinal compression associated with basal resistance from the nearby sticky spot.

The widespread occurrence of snow-ice formation on the pack ice plays a critical role in the mass balance of Antarctic sea ice. The stable isotope composition, ice texture and salinity of eight ice cores, obtained from the Amundsen Sea during the Oden Southern Ocean 2010/11 expedition from late December 2010 to January 2011, were investigated to illustrate the snow-ice growth process and its contribution to sea-ice development. Most previous research has utilized δ18O as an index tracer to determine the percentages of core length that contain meteoric water, i.e. snow ice. However, this standard practice of snow-ice identification might be biased due to normally low-resolution isotopic measurements and mixing/diffusion processes between the snow ice and underlying ice layers. Snow-ice contributions in these ice cores based instead on an updated isotope mixing model are also presented. Depth profiles of ice texture and salinity are described to serve as representations of the structures of these ice cores. Our isotope mixing model produced an average of 15.9% snow-ice contribution for pack ice in the Amundsen Sea, and meteoric water occupying 40% of snow-ice mass for all ice stations. These results are compared to previous investigations of snow-ice occurrence around Antarctica.

The impact of April sea-ice thickness (SIT) initialization on the predictability of September sea-ice extent (SIE) is investigated based on a series of perfect model ensemble experiments using the MIROC5.2 climate model. Ensembles with April SIT initialization can accurately predict the September SIE for greater lead times than in cases without the initialization – up to 2 years ahead. The persistence of SIT correctly initialized in April contributes to the skilful prediction of SIE in the first September. On the other hand, errors in the initialization of SIT in April cause errors in the predicted sea-ice concentration and thickness in the Pacific sector from July to September and consequently influence the predictive skill with respect to SIE in September. The present study suggests that initialization of the April SIT in the Pacific sector significantly improves the accuracy of the September SIE forecasts by decreasing the errors in sea-ice fields from July to September.

The Askaryan Radio Array (ARA) experiment at the South Pole is designed to detect high-energy neutrinos which, via in-ice interactions, produce coherent radiation at frequencies up to 1000 MHz. Characterization of ice birefringence, and its effect upon wave polarization, is proposed to enable range estimation to a neutrino interaction and hence aid in neutrino energy reconstruction. Using radio transmitter calibration sources, the ARA collaboration recently measured polarization-dependent time delay variations and reported significant time delays for trajectories perpendicular to ice flow, but not parallel. To explain these observations, and assess the capability for range estimation, we use fabric data from the SPICE ice core to model ice birefringence and construct a bounding radio propagation model that predicts polarization time delays. We compare the model with new data from December 2018 and demonstrate that the measurements are consistent with the prevailing horizontal crystallographic axis aligned near-perpendicular to ice flow. The study supports the notion that range estimation can be performed for near flow-perpendicular trajectories, although tighter constraints on fabric orientation are desirable for improving the accuracy of estimates.

The Antarctic Impulsive Transient Antenna (ANITA) balloon experiment was designed to detect radio signals initiated by high-energy neutrinos and cosmic ray (CR) air showers. These signals are typically discriminated by the polarization and phase inversions of the radio signal. The reflected signal from CRs suffer phase inversion compared to a direct ‘tau neutrino’ event. In this paper, we study subsurface reflection, which can occur without phase inversion, in the context of the two anomalous up-going events reported by ANITA. It is found that subsurface layers and firn density inversions may plausibly account for the events, while ice fabric layers and wind ablation crusts could also play a role. This hypothesis can be tested with radar surveying of the Antarctic region in the vicinity of the anomalous ANITA events. Future experiments should not use phase inversion as a sole criterion to discriminate between down-going and up-going events, unless the subsurface reflection properties are well understood.

Airborne electromagnetic induction sensors have demonstrated their extensive capacities to measure sea-ice thickness distributions. However, biases can emerge when comparing these 1-D measurements to a broader 2-D regional scale due to the spatial anisotropy inherent to sea-ice cover. Automated processing of available sea-ice maps could significantly ease the decision on how to set up an optimised flight pattern, which would result in representative ice thickness numbers for the region. In this study, first we investigate the extent to which the sea-ice anisotropy can influence the representativeness of an airborne survey compared to the regional situation. Second, we propose a method to process sea-ice maps prior to flights to help preparing the most representative flight plan possible for the local area. The method is based on automated segmentation of radar satellite images and extensive simulation of flight transects over the image. The spatial analysis of these transects enables for the identification of the most representative survey trajectories for the area. The method was applied for seven different synthetic aperture radar satellite images over Arctic sea ice north of Svalbard. The results indicate that the proposed method improved the representativeness of the airborne survey by identifying the most suitable transect over the ice pack.

Our understanding of ice algal responses to the recent changes in Arctic sea ice is impeded by limited field observations. In the present study, environmental characteristics of the landfast sea-ice zone as well as primary production and macromolecular composition of ice algae and phytoplankton were studied in the Kitikmeot Sea near Cambridge Bay in spring 2017. Averaged total chlorophyll-a (Chl-a) concentration was within the lower range reported previously for the same region, while daily carbon uptake rates of bottom-ice algae were significantly lower in this study than previously reported for the Arctic. Based on various indicators, the region's low nutrient concentrations appear to limit carbon uptake rates and associated accumulation of bottom-ice algal biomass. Furthermore, the lipids-dominant biochemical composition of bottom-ice algae suggests strong nutrient limitation relative to the distinctly different carbohydrates-dominant composition of phytoplankton. Together, the results confirm strong nitrate limitation of the local marine system.

Clouds play an important role in the Arctic surface radiative budget, impacting the seasonal evolution of Arctic sea-ice cover. We explore the large-scale impacts of springtime and early summer (March through July) cloud and radiative fluxes on sea ice by comparing these fluxes to seasonal ice volume losses over the central Arctic basin, calculated for available observational years 2004–2007 (ICESat) and 2011–2017 (CryoSat-2). We also supplement observation data with sea-ice volume computed from the Pan-Arctic Ice–Ocean Modeling and Assimilation System (PIOMAS) during summer months. We find that the volume of sea ice lost over the melt season is most closely related to observed downwelling longwave radiation in March and early summer (June and July) longwave cloud radiative forcing, which together explain a large fraction of interannual variability in seasonal sea-ice volume loss (R2 = 0.71, p = 0.007). We show that downwelling longwave fluxes likely impact the timing of melt onset near the sea-ice edge, and can limit the magnitude of ice thickening from March to April. Radiative fluxes in June and July are likely critical to seasonal volume loss because modeled data show the greatest ice volume reductions occur during these months.

Stationary ice-penetrating radar (sIPR) systems can be used to monitor temporal changes in electromagnetically sensitive properties of glaciers and ice sheets. We describe a system intended for autonomous operation in remote glacial environments, and document its performance during deployments in cold and temperate settings. The design is patterned after an existing impulse radar system, with the addition of a fibre-optic link and timing module to control transmitter pulses, a micro-UPS (uninterruptable power supply) to prevent uncontrolled system shutdown and a customized satellite telemetry scheme. Various implementations of the sIPR were deployed on the Kaskawulsh Glacier near an ice-marginal lake in Yukon, Canada, for 44–77 days in summers 2014, 2015 and 2017. Pronounced perturbations to englacial radiostratigraphy were observed commensurate with lake filling and drainage, and are interpreted as changes in englacial water storage. Another sIPR was deployed in 2015–2016 on ice island PII-A-1-f, which originated from the Petermann Glacier in northwest Greenland. This system operated autonomously for almost a year during which changes in thickness of the ice column were clearly detected.

The surface flooding of Antarctic sea ice in summer covers 50% or more of the sea-ice area in the major summer ice packs, the western Weddell and the Bellingshausen-Amundsen Seas. Two CRREL ice mass-balance buoys were deployed on the Amundsen Sea pack in late December 2010 from the icebreaker Oden, bridging the summer period (January–February 2011). Temperature records from thermistors embedded vertically in the snow and ice showed progressive increases in the depth of the flooded layer (up to 0.3–0.35 m) on the ice cover during January and February. While the snow depth was relatively unchanged from accumulation (<10 cm), ice thickness decreased by up to a meter from bottom melting during this period. Contemporaneous with the high bottom melting, under-ice water temperatures up to 1°C above the freezing point were found. The high temperature arises from solar heating of the upper mixed layer which can occur when ice concentration in the local area falls and lower albedo ocean water is exposed to radiative heating. The higher proportion of snow ice found in the Amundsen Sea pack ice therefore results from both winter snowfall and summer ice bottom melt found here that can lead to extensive surface flooding.

The Arctic sea-ice cover has undergone a significant decline in recent decades. The melt season is starting earlier, ice is thinner and seasonal ice dominates. Here we examine the effects of these changes on the solar heat input to the upper ocean in ice-covered Arctic waters from 1985 to 2014. Satellite observations of ice concentration, onset dates of melt and freeze-up and ice age, are combined with ice thicknesses from the PIOMAS model and incident solar irradiance from reanalysis products to calculate the contributions of open ocean and ice to the solar heat input in the upper ocean. Of the total, 86% of the area has positive trends for solar heat input to the ocean through leads due to decreases in ice concentration. Only 62% of the area shows positive trends of solar heat input to the ocean explicitly through the ice. Positive trends are due to thinning ice, while negative trends occur in regions where the ice-free season has lengthened. The annual total solar heat input to the ocean exhibits positive trends in 82% of the area. The spatial pattern of the cumulative annual total solar heat input is similar to the pattern of solar heat input directly to leads.

The Polar Sea Ice Topography REconstruction System, or PSITRES, is a 3D camera system designed to continuously monitor an area of ice and water adjacent to an ice-going vessel. Camera systems aboard ships in the polar regions are common; however, the application of computer vision techniques to extract high-level information from the imagery is infrequent. Many of the existing systems are built for human involvement throughout the process and lack automation necessary for round the clock use. The PSITRES was designed with computer vision in mind. It can capture images continuously for days on end with limited oversight. We have applied the system in different ice observing scenarios. The PSITRES was deployed on three research expeditions in the Arctic and Subarctic, and we present applications in measuring ice concentration, melt pond fraction and presence of algae. Systems like PSITRES and the computer vision algorithms applied represent steps toward automatically observing, evaluating and analyzing ice and the environment around ships in ice-covered waters.

Several autonomous phase-sensitive radio-echo sounders (ApRES) were deployed at Greenland glaciers to investigate ice deformation. Different attenuation settings were tested and it was observed that, in the presence of clipping of the deramped ApRES signal, each setting produced a different result. Specifically, higher levels of clipping associated with lower attenuation produced an apparent linear increase of diurnal vertical cumulative displacement with depth, and obscured the visibility of the basal reflector in the return amplitude. An example with a synthetic deramped signal confirmed that these types of artifacts result from the introduction of harmonics from square-wave-like features introduced by clipping. Apparent linear increase of vertical displacement with depth occurs when the vertical position of a near-surface internal reflector changes in time. Artifacts in the return amplitude may obscure returns from internal reflectors and the basal reflector, making it difficult to detect thickness evolution of the ice and to correctly estimate vertical velocities. Variations in surface melt during ApRES deployments can substantially modulate the received signal strength on short timescales, and we therefore recommend using higher attenuator settings for deployments in such locations.

Columnar sea ice grows with an interface of tiny parallel ice plates, the distance of which is known as plate spacing. While it has been proposed as a fundamental microstructure scale of sea ice, the physics behind its formation has not been fully understood. Here the problem is analysed on the basis of morphological stability theory to propose a model that results in a physically consistent prediction of the relationship between the plate spacing a0 and growth velocity V. The relationship may be divided into two regimes. In the diffusive regime, for V above ≈2 × 10−4 cm s−1 one finds a0 ~ V−2/3 to first order. In the convective regime, the extent of diffusive boundary layer is controlled by solutal convection near the interface, which leads to the proportionality a0 ~ V−1/3. From a comparison to observations it is evident that the plate spacing is predictable over 5 orders of magnitude in the growth velocity, covering the range from fast laboratory ice growth to slow accretion at the bottom of marine ice shelves. The predictability opens new paths towards concise modelling of marine and sea-ice microstructure and physical properties.

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.

The role of synoptic-scale cyclones in the trends and variability of Arctic sea ice conditions has remained uncertain. In recognition, we conduct a systematic investigation of how sea-ice concentration (SIC) changes with cyclone passage, including all individual storms that pass over any part of the region's ice pack. For all seasons, especially summer and autumn, we find a pattern of higher ice concentration after a region is influenced by a cyclone compared to when it is not, primarily due to thermodynamic effects. During warm months, cyclones appear to slow the general day-to-day decline in concentration; in cold months, cyclones augment the day-to-day increase. These relationships are changing over time, with cyclone-associated concentration changes becoming less distinct from overall changes. Cyclone effects on ice divergence are spatially variable; computed fields are noisy. In summer, these dynamic effects of cyclone passage generally decrease SIC, but are outweighed by the thermodynamic effects (e.g., reductions in air temperature, shortwave radiation). In autumn, cyclone-associated concentration changes are not as easily explained by observed cyclone conditions. Key questions remain regarding the extent to which our findings are influenced by artifacts of surface melt and weather effects on the passive microwave retrievals.

Radio-echo sounding (RES) can be used to understand ice-sheet processes, englacial flow structures and bed properties, making it one of the most popular tools in glaciological exploration. However, RES data are often subject to ‘strip noise’, caused by internal instrument noise and interference, and/or external environmental interference, which can hamper measurement and interpretation. For example, strip noise can result in reduced power from the bed, affecting the quality of ice thickness measurements and the characterization of subglacial conditions. Here, we present a method for removing strip noise based on combined wavelet and two-dimensional (2-D) Fourier filtering. First, we implement discrete wavelet decomposition on RES data to obtain multi-scale wavelet components. Then, 2-D discrete Fourier transform (DFT) spectral analysis is performed on components containing the noise. In the Fourier domain, the 2-D DFT spectrum of strip noise keeps its linear features and can be removed with a ‘targeted masking’ operation. Finally, inverse wavelet transforms are performed on all wavelet components, including strip-removed components, to restore the data with enhanced fidelity. Model tests and field-data processing demonstrate the method removes strip noise well and, incidentally, can remove the strong first reflector from the ice surface, thus improving the general quality of radar data.