The evolution of melt ponds on Arctic sea ice in summer is one of the main factors that affect sea-ice albedo and hence the polar climate system. Due to the different spectral properties of open water, melt pond and sea ice, the melt pond fraction (MPF) can be retrieved using a fully constrained least-squares algorithm, which shows a high accuracy with root mean square error ~0.06 based on the validation experiment using WorldView-2 image. In this study, the evolution of ponds on first-year and multiyear ice in the Canadian Arctic Archipelago was compared based on Sentinel-2 and Landsat 8 images. The relationships of pond coverage with air temperature and albedo were analysed. The results show that the pond coverage on first-year ice changed dramatically with seasonal maximum of 54%, whereas that on multiyear ice changed relatively flat with only 30% during the entire melting period. During the stage of pond formation, the ponds expanded rapidly when the temperature increased to over 0°C for three consecutive days. Sea-ice albedo shows a significantly negative correlation (R = −1) with the MPF in melt season and increases gradually with the refreezing of ponds and sea ice.

Spatial and temporal trends of remotely sensed sea-ice cover, sea surface temperatures, chlorophyll-a concentration and primary production in the Baffin Bay, Davis Strait and Labrador Sea were analyzed for the 1998–2017 period. We found spatial variability in the trends of these cryospheric, biologic and oceanographic phenomena. For example, in the northern Baffin Bay, we observed decreases in annual sea-ice persistence, yet increases along the Labrador Sea-ice edge during winter, with the latter having significant correlations with broader atmospheric patterns. In general, we observed increases in summer sea surface temperatures across the study region, except a small area of cooling along the southern Greenlandic coast. We also found significant negative trends in April chlorophyll-a and primary production followed by significant positive trends for both biological phenomena in May, owing to anomalously high values in 2014 and 2015. Notably, we found a significant positive correlation between days of monthly sea ice presence in April with May primary production quantities. Finally, we found a significant positive trend in total annual primary production over the study period. This novel finding suggests an important relationship between the timing of breakup along the sea-ice edge and peaks in biological production.

This work is a collection of radar equations for low-frequency radar sounding and radar in general that emphasize the form of the radar equation for different target and source geometries. This is meant as a handbook for scientists and engineers that work with or analyze radar sounder systems and interpret radar sounding echoes. Lookup tables summarize the results and derivations are provided for each equation.

The Arctic winter seasonal sea ice (WSSI) concentration from 1979 to 2019 is derived from passive microwave data. Based on Empirical Orthogonal Function (EOF) analysis, the WSSI time series includes regionally different trends, abrupt shifts and interannual variations. The time series of the first EOF mode (PC1) mainly represents the WSSI trend, which is characterized by an increase, particularly in the Pacific sector. PC1 confirms two abrupt shifts in WSSI in 1989 and 2007, with a variance of 31%. After 2007, the large-scale atmospheric circulation anomaly shows a strengthened wavenumber 3 structure at high latitudes associated with a mid-tropospheric low-pressure anomaly in central and western Siberia and a high-pressure anomaly in eastern Siberia in summer and autumn. These patterns have promoted the increased transport of moist static energy to the central Arctic and contributed to increased near-surface air temperatures that may enhance ice melting in summer and reduce ice growth in autumn and winter. The changes in ice melt and growth have had opposite effects in the Pacific and Atlantic sectors: WSSI has increased in the Pacific sector due to the replacement of multi-year ice by WSSI, and decreased in the Atlantic sector due to the replacement of WSSI by open water.

We examine the relative effect of warming events (storms) and snow cover on thermodynamic growth of Arctic sea ice in winter. We use a 1-D snow and ice thermodynamic model to perform sensitivity experiments. Observations from the winter period of the Norwegian young sea ICE (N-ICE2015) campaign north of Svalbard are used to initiate and force the model. The N-ICE2015 winter was characterized by frequent storm events that brought pulses of heat and moisture, and a thick snow cover atop the sea ice (0.3–0.5 m). By the end of the winter, sea-ice bottom growth was negligible. We show that the thermodynamic effect of storms to the winter sea-ice growth is controlled by the amount of snow on sea ice. For 1.3 m initial ice thickness, the decrease in ice growth caused by the warming events ranged from −1.4% (for 0.5 m of snow) to −7.5% (for snow-free conditions). The decrease in sea-ice growth caused by the thick snow (0.5 m) was more important, ranging from −17% (with storms) to −23% (without storms). The results showcase the critical role of snow on winter Arctic sea-ice growth.

The yearly paired process of slow growth and rapid melt of some 15 million square kilometers of Antarctic sea ice takes place with a regular asymmetry; the process has been linked to the relationship of the position of the ice edge with the band of low pressure that circles the continent between 60° and 70°S. In autumn, winds to the north of the low-pressure band slow the advancing ice edge. In summer, Ekman divergence created by opposing winds on either side of the low-pressure band opens up warm water regions that rapidly melt sea ice. We use the 40 ensemble members from the CESM-LENS historical run (1920–2005) to examine the relationship between the asymmetry in the annual cycle and the position and intensity of the low-pressure band. CESM-LENS reproduces the magnitude of the annual cycle of Antarctic sea ice extent with a short lag (2 weeks). Melt rate is the characteristic of the annual cycle that varies the most. Our results provide evidence that lower pressure leads to increased melt rates, which supports the importance of the role of divergence in increasing the melt rate of Antarctic sea ice. The role of winds during the growing season remains unquantified.

Sea ice, which forms in polar and nonpolar areas, transmits light to ice-associated (sympagic) algal communities. To noninvasively study the distribution of sea-ice algae, empirical relations to estimate its biomass from under-ice hyperspectral irradiance have been developed in the Arctic and Antarctica but lack for nonpolar regions. This study examines relationships between normalised difference indices (NDI) calculated from hyperspectral transmittance and sympagic algal biomass in the nonpolar Saroma-ko Lagoon. We analysed physico-biogeochemical properties of snow and land-fast sea ice supporting 27 paired bio-optical measurements along three transects covering an area of over 250 m × 250 m in February 2019. Snow depth (0.08 ± 0.01 m) and ice-bottom brine volume fraction (0.21 ± 0.02) showed low (0.06) and high (0.58) correlations with sea-ice core bottom section chlorophyll a (Chl. a), respectively. Spatial analyses unveiled the patch size of sea-ice Chl. a to be ~65 m, which is in the same range reported from previous studies. A selected NDI (669, 596 nm) explained 63% of algal biomass variability. This reflects the bio-optical properties and environmental conditions of the lagoon that favour the wavelength pair in the orange/red part of the spectrum and suggests the necessity of a specific bio-optical relationship for Saroma-ko Lagoon.

The Northeast Greenland Ice Stream (NEGIS) is an important dynamic component for the total mass balance of the Greenland ice sheet, as it reaches up to the central divide and drains 12% of the ice sheet. The geometric boundary conditions and in particular the nature of the subglacial bed of the NEGIS are essential to understand its ice flow dynamics. We present a record of more than 8000 km of radar survey lines of multi-channel, ultra-wideband radio echo sounding data covering an area of 24 000 km2, centered on the drill site for the East Greenland Ice-core Project (EGRIP), in the upper part of the NEGIS catchment. Our data yield a new detailed model of ice-thickness distribution and basal topography in the region. The enhanced resolution of our bed topography model shows features which we interpret to be caused by erosional activity, potentially over several glacial–interglacial cycles. Off-nadir reflections from the ice–bed interface in the center of the ice stream indicate a streamlined bed with elongated subglacial landforms. Our new bed topography model will help to improve the basal boundary conditions of NEGIS prescribed for ice flow models and thus foster an improved understanding of the ice-dynamic setting.

The Ross Sea is known for showing the greatest sea-ice increase, as observed globally, particularly from 1979 to 2015. However, corresponding changes in sea-ice thickness and production in the Ross Sea are not known, nor how these changes have impacted water masses, carbon fluxes, biogeochemical processes and availability of micronutrients. The PIPERS project sought to address these questions during an autumn ship campaign in 2017 and two spring airborne campaigns in 2016 and 2017. PIPERS used a multidisciplinary approach of manned and autonomous platforms to study the coupled air/ice/ocean/biogeochemical interactions during autumn and related those to spring conditions. Unexpectedly, the Ross Sea experienced record low sea ice in spring 2016 and autumn 2017. The delayed ice advance in 2017 contributed to (1) increased ice production and export in coastal polynyas, (2) thinner snow and ice cover in the central pack, (3) lower sea-ice Chl-a burdens and differences in sympagic communities, (4) sustained ocean heat flux delaying ice thickening and (5) a melting, anomalously southward ice edge persisting into winter. Despite these impacts, airborne observations in spring 2017 suggest that winter ice production over the continental shelf was likely not anomalous.

Geophysical systems are often assumed to follow Gaussian probability density functions; however, deviations from Gaussian behaviour can shed light on the underlying dynamics. For the large-scale motion of the Arctic sea ice, such deviations have been interpreted as signatures of structure in dynamic flow fields. In this study, we use higher-order moments (skewness and kurtosis) to identify spatiotemporal changes in the Beaufort Gyre (BG) and the Transpolar Drift (TPD) sea-ice drift patterns. Higher-order moments of satellite-derived ice drift speeds are examined over the winter period of 2006–2017 to describe the persistence of features like the BG and TPD, and their variation over time. Index patterns indicate that the periphery of the BG can be identified by a combination of high positive skewness and high kurtosis in the ice drift time series on an annual basis.

Recent advances in the measurement of bedforms beneath active ice streams have been made using ground-based grid profiling using impulse radar systems operating with centre frequencies in the 3–5 MHz range. Surveys of Rutford Ice Stream and Pine Island Glacier have shown that features such as mega-scale glacial lineations with topographic relief of as little as 3 m can be traced for many kilometres downstream under more than 2 km of fast-moving ice. In the discussion of these data, it is often asked ‘How is it possible to map such fine-scale topography with such a low-frequency radar’. In answering that question, the key point is the distinction between the precision of a radar range measurement to a single, isolated reflective interface and the ability to resolve the presence of two closely-spaced interfaces of similar reflectivity (commonly referred to as the vertical resolution). This paper will discuss and illustrate this distinction and use the case study of data acquired over Pine Island Glacier to examine the limits of precision of the radar range measurement.

This paper provides the first evidence for sea-ice formation in the Cordillera Darwin (CD) fjords in southern Chile, which is farther north than sea ice has previously been reported for the Southern Hemisphere. Initially observed from a passenger plane in September 2015, the presence of sea ice was then confirmed by aerial reconnaissance and subsequently identified in satellite imagery. A time series of Sentinel-1 and Landsat-8 images during austral winter 2015 was used to examine the chronology of sea-ice formation in the Cuevas fjord. A longer time series of imagery across the CD was analyzed from 2000 to 2017 and revealed that sea ice had formed in each of the 13 fjords during at least one winter and was present in some fjords during a majority of the years. Sea ice is more common in the northern end of the CD, compared to the south where sea ice is not typically present. Is suggested that surface freshening from melting glaciers and high precipitation reduces surface salinity and promotes sea-ice formation within the semi-enclosed fjord system during prolonged periods of cold air temperatures. This is a unique set of initial observations that identify questions for future research in this remote area.

Quantifying the rate of wave attenuation in sea ice is key to understanding trends in the Antarctic marginal ice zone extent. However, a paucity of observations of waves in sea ice limits progress on this front. We deployed 14 waves-in-ice observation systems (WIIOS) on Antarctic sea ice during the Polynyas, Ice Production, and seasonal Evolution in the Ross Sea expedition (PIPERS) in 2017. The WIIOS provide in situ measurement of surface wave characteristics. Two experiments were conducted, one while the ship was inbound and one outbound. The sea ice throughout the experiments generally consisted of pancake and young ice <0.5 m thick. The WIIOS survived a minimum of 4 d and a maximum of 6 weeks. Several large-wave events were captured, with the largest recorded significant wave height over 9 m. We find that the total wave energy measured by the WIIOS generally decays exponentially in the ice and the rate of decay depends on ice concentration.

Basal motion of ice sheets depends in part on the roughness and material properties of the subglacial bed and the occurrence of water. To date, basal motion represents one of the largest uncertainties in ice-flow models. It is that component of the total flow velocity that can change most rapidly and can, therefore, facilitate rapid variations in dynamic behaviour. In this study, we investigate the subglacial properties of the East Antarctic Ice Sheet by statistically analysing the roughness of the bed topography, inferred from radio-echo sounding measurements. We analyse two sets of roughness parameters, one derived in the spatial and the other in the spectral domain, with two roughness parameters each. This enables us to compare the suitability of the four roughness parameters to classify the subglacial landscapes below the ice sheet. We further investigate the relationship of the roughness parameters with observed surface flow velocity and modelled basal temperatures of the ice sheet. We find that one of the roughness parameters, the Hurst exponent derived in the spatial domain, coincides with the thermal condition at the base of the ice sheet for slow flow velocities and varies with flow velocity.

All radar power interpretations require a correction for attenuative losses. Moreover, radar attenuation is a proxy for ice-column properties, such as temperature and chemistry. Prior studies use either paired thermodynamic and conductivity models or the radar data themselves to calculate attenuation, but there is no standard method to do so; and, before now, there has been no robust methodological comparison. Here, we develop a framework meant to guide the implementation of empirical attenuation methods based on survey design and regional glaciological conditions. We divide the methods into the three main groups: (1) those that infer attenuation from a single reflector across many traces; (2) those that infer attenuation from multiple reflectors within one trace; and (3) those that infer attenuation by contrasting the measured power from primary and secondary reflections. To assess our framework, we introduce a new ground-based radar survey from South Pole Lake, comparing selected empirical methods to the expected attenuation from a temperature- and chemistry-dependent Arrhenius model. Based on the small surveyed area, lack of a sufficient calibration surface and low reflector relief, the attenuation methods that use multiple reflectors are most suitable at South Pole Lake.

Results examining variations in the ice extent along the Norwegian coastline based on the analysis of Moderate Resolution Imaging Spectroradiometer (MODIS) images from 2001 to 2019, February through May, are presented. A total of 386 fjords and coastal areas were outlined and grouped into ten regions to assess seasonal and long-term trends in ice extent. In addition, three fjords were examined to investigate how ice extent may vary over short distances (<100 km). Of the 386 outlined, 47 fjords/coastal areas held >5 km2 of ice at least once between 2001 and 2019. Over this span of time, no statistically significant trend in ice extent is found for all ten regions; however, variations between regions and years are evident. Ice extent is assessed through comparison to three weather variables – freezing degree days (FDD), daily new snowfall and daily freshwater supply from rainfall plus snowmelt. Six out of ten regions are significantly positively correlated (p < 0.05) to FDD. In addition, ice in two regions is significantly positively correlated to daily new snowfall, and in one region negatively correlated to rainfall plus snowmelt. The importance of fjord geometry and bathymetry as well as other weather variables including wind is discussed.

Although altimeters have been widely used to monitor the spatiotemporal variation of sea-ice thickness, they are unable to separate sea-ice freeboard from snow depth. We use a floating GPS deployed on sea ice to derive the freeboard and snow depth near China's Zhongshan Station. Our results show that the standalone floating GPS can monitor freeboard with a precision of 4.2 cm. If time-varying dynamic ocean topography provided by, for example, a bottom pressure gauge is available, then the precision of GPS-derived freeboard can improve to 1.3 cm. The daily snow depth inverted by GPS interferometric reflectometry captures three precipitation events during our experiment, showing that the floating GPS can monitor the variation in snow depth and observe the freeboard variation at the same time. By studying the relationship between freeboard, snow depth and sea-ice thickness, we find that sea-ice thickness will be greatly underestimated by the negative single-point freeboard under the assumption of hydrostatic equilibrium. As a supplement to existing technologies, the GPS-derived freeboard and snow depth can be used both to evaluate the altimeter observations directly and to improve our understanding of the real-time variation of freeboard and snow depth in the experimental area.

Numerical ice-sheet models are commonly matched to surface ice velocities from InSAR measurements by modifying basal drag, allowing the flow and form of the ice sheet to be simulated. Geophysical measurements of the bed are rarely used to examine if this modification is realistic, however. Here, we examine radio-echo sounding (RES) data from the Weddell Sea sector of West Antarctica to investigate how the output from a well-established ice-sheet model compares with measurements of the basal environment. We know the Weddell Sea sector contains the Institute, Möller and Foundation ice streams, each with distinct basal characteristics: Institute Ice Stream lies partly over wet unconsolidated sediments, where basal drag is very low; Möller Ice Stream lies on relatively rough bed, where basal drag is likely larger; and Foundation Ice Stream is controlled by a deep subglacial trough with flow-aligned bedrock landforms and smooth unconsolidated sediments. In general, the ice-sheet model represents each ice-stream system well. We also find that ice velocities do not match perfectly in some locations, and that adjustment of the boundaries of low basal drag, to reflect RES evidence, should improve model performance. Our work showcases the usefulness of RES in calibrating ice-sheet model results with observations of the bed.