Predicting the terminus position of Greenland’s marine-terminating glaciers remains a major challenge for ice-sheet modeling. Improving model parameterizations requires understanding the dominant drivers of terminus change. Here, we employ the Liang–Kleeman information flow method to quantify the contributions of bed slope, sea-ice concentration, runoff, air temperature and ocean thermal forcing to changes in glacier terminus position and frontal ablation. We focus on three major marine-terminating glaciers, Jakobshavn Isbræ (JI), Kangerlussuaq Glacier (KG) and Helheim Glacier (HG), between 2000 and 2021. Our analysis reveals that atmospheric and oceanic forcings exert comparable influence on glacier retreat, with their causal effects demonstrating pronounced seasonal variability. These drivers dominate during the summer months while exerting minimal impact in winter. Analysis of ocean thermal forcing highlights the 0–100 m surface layer as the primary driver of terminus variability for JI and KG, underscoring the critical role of surface-ocean processes. Notably, despite theoretical relevance, bed slope at kilometer scales shows no significant causal influence, emphasizing the nonlinearity and complexity of topographic controls. Finally, we develop a parameterized statistical model that explains 53–66% of the variance in frontal ablation for the highly responsive JI and KG, but only 38% for the less sensitive HG.

The Arctic is a hotspot for climate warming, making it crucial to quantify the sea level rise contribution from its ice masses. Novaya Zemlya's ice caps are the largest glacier complex in Europe and are a major contributor to contemporary sea level rise. Here we show that Novaya Zemlya outlet glaciers on the Barents Sea coast respond rapidly and consistently to oceanic forcing at annual timescales, likely due to their exposure to Atlantic Water variability. Glaciers on the Kara Sea show more variable response, likely reflecting their reduced exposure to Atlantic Water. Data demonstrate that the pause in glacier retreat previously observed on Novaya Zemlya between 2013 and 2015 has not persisted and that these changes correspond to ocean temperature variability on the Barents Sea coast. We document a marked shift to warmer air and ocean temperatures, and reduced sea ice concentrations from 2005 onwards. Although we identify ocean warming as the primary trigger for glacier retreat, we suggest that multi-year thinning, driven by the shift towards warmer air temperatures since 2005, pre-conditioned Novaya Zemlya's glaciers to retreat. Despite commonality in the timing of outlet glacier retreat, the magnitude is highly variable during rapid retreat phases, which we attribute to glacier-specific factors.

We describe the annual pattern of frontal ablation driven by submarine melting mechanisms at the Hansbreen terminus: these are reflected in the intensity and spatial distribution of calving events. Analysis of time-lapse images of the Hansbreen front in conjunction with oceanographic and meteorological data shows that calving intensity is driven primarily by seawater temperature. Regression analysis also highlights the importance of air temperature, which we take to be a proxy for surface ablation and subglacial discharge. This, combined with seasonal changes in ice cliff tortuosity and the increasing significance of wave motion outside the ablation season, enabled us to determine seasonal changes in the mechanisms of ice cliff undercutting by submarine melting. While submarine melting controlled by estuarine circulation primarily drives frontal ablation in summer, wave-driven melting at the waterline is more important outside the ablation season. During winter, ice cliff undercutting by melting is suspended by low seawater temperature, negligible subglacial water discharge and sea-ice cover. The most intense frontal ablation, recorded in summer, was related to higher sea temperature and vigorous estuarine circulation.

To increase inclusivity, diversity, equity and accessibility in Antarctic science, we must build more positive and inclusive Antarctic field work environments. The International Thwaites Glacier Collaboration (ITGC) has engaged in efforts to contribute to that goal through a variety of activities since 2018, including creating an open-access ‘Field and Ship Best Practices’ guide, engaging in pre-field season team dynamics meetings, and surveying post-field season reflections and experiences. We report specific actions taken by ITGC and their outcomes. We found that strong and supported early career researchers brought new and important perspectives regarding strategies for transforming culture. We discovered that engaged and involved senior leadership was also critical for expanding participation and securing funding to support efforts. Pre-field discussions involving all field team members were particularly helpful for setting expectations, improving sense of belonging, describing field work best practices, and co-creating a positive work culture.

Processes driving changes in sea-ice seasonality and sea-ice thickness were explored for a ‘warm-shelf’ region along the West Antarctic Peninsula using vertically coupled sea-ice-ocean thermodynamic simulations, with and without assimilated satellite sea-ice observations and moored ocean temperature observations. Simulations with assimilated sea-ice observations permitted investigation of surface [thermodynamic and dynamic (e.g., wind-driven)] processes affecting sea-ice thickness and seasonality. Assimilation of quasi-weekly variability in the depth and temperature of the deep warm pycnocline permitted examination of subsurface processes affecting sea-ice. Simulations using assimilated sea-ice observations (and implied motion) always produced greater surface heat fluxes and overall thinner sea ice. Assimilating seasonal and quasi-weekly variability in the depth and temperature of the pycnocline modified the start of the sea-ice season by −23 to +1 d, and also modified the sea ice thickness/seasonality to be thinner/shorter or thicker/longer at sub-seasonal and seasonal timescales, highlighting a mechanism where a shoaling pycnocline enhanced upward deep-water heat fluxes as transient surface-induced turbulence had a greater effect on a reduced mixed layer volume. The observed interplay of surface, subsurface, and sea-ice modulation of ocean-atmosphere heat transfer underscores the importance of representing the interaction between sea-ice concentration and upper ocean variability in climate projections.

Arctic rain-on-snow (ROS) events can have significant impacts on Arctic wildlife and socio-economic systems. This study addresses the meteorology of two different Arctic ROS events. The first, occurring near Nuuk, Greenland, generated significant impacts, including slush avalanches. The second, less severe, event occurred within the community of Iqaluit, Nunavut, Canada. This research utilizes atmospheric reanalysis, automated surface observation station data and atmospheric soundings to determine the meteorological conditions driving these events and the differences between each case. In both cases, atmospheric blocking played a leading role in ROS initiation, with atmospheric rivers – narrow bands of high water vapor transport, typically originating from the tropics and subtropics – having both direct and indirect effects. Cyclone-induced low-level jets and resultant ‘warm noses’ of higher air temperatures and moisture transport were other key features in ROS generation. To our knowledge, our study is the first to visualize how the varying strength and manifestation of these coupled features contribute to differences in the severity of Arctic ROS events. The meteorological drivers identified here find support from other studies on Arctic ROS events and are similar to weather features associated with Arctic precipitation events of extreme magnitude.

Using data from three automatic weather stations (LGB69, Eagle and Dome A) from distinctly different climatological zones along the CHINARE (Chinese National Antarctic Research Expedition) traverse route from Zhongshan Station to Dome A, we investigated the characteristics of meteorological conditions and subsurface heat conduction. Spatial analysis indicated decreasing trends in air temperature, relative humidity and wind speed from the coastal katabatic wind zone to the inland plateau region, and air temperatures clearly showed a strong daily variability in winter, suggesting the effect from the fluctuation in the Antarctic atmospheric system. We also analyzed the optimal response time of the 1 and 3 m depth snow temperatures to the 0.1 m depth snow temperature for each site under clear/overcast and day/night situations. This showed an important enhancement to the heat transfer from shortwave radiation penetration. Using an iterative optimization method, we estimated the subsurface heat conduction variations along the transect. This was ~3–5 W m–2. Multiple maxima in daily mean subsurface fluxes were found in winter, with a typical value above 2 W m–2, while a single minimum value under –2 W m–2 was found in summer. On an annual scale, a larger mean loss of subsurface heat conduction was observed in the inland plateau compared to in the coastal katabatic area. Finally, we discussed the possible influences of turbulent and radiant transport on the vertical heat response and confirmed the wind enhancement on the growth of thermal conductivity. This preliminary study provides a brief perspective and an important reference for studying subsurface heat conduction in inland areas of Antarctica.

In situ observations of wave attenuation by sea ice are required to develop and validate wave–ice interaction parameterizations in coupled wave models. To estimate ice-induced wave attenuation in the field, the wave field is typically assumed to be stationary. In this study we investigate the validity of this assumption by creating a synthetic wave field in sea ice for different attenuation rates. We observe that errors in estimates of the wave attenuation rates are largest when attenuation rates are small or temporal averaging periods are short. Moreover, the adoption of the wave stationarity assumption can lead to negative estimates of the instantaneous wave attenuation rate. These apparent negative values should therefore not be attributed to wave growth or erroneous measurements a priori. Surprisingly, we observe that the validity of the wave stationarity assumption is irrelevant to the accuracy of estimates of wave attenuation rates as long as the temporal averaging period is taken sufficiently long. This may provide opportunities in using satellite-derived products to estimate wave attenuation rates in sea ice at global scales.

Sea-ice drift in the Antarctic marginal ice zone (MIZ) is discussed using data from a 4-month-long drift of a buoy deployed on a pancake ice floe during the winter sea-ice expansion. We demonstrate increased meandering and drift speeds, and changes in the dynamical regimes of the absolute dispersion during cyclone activity, together with high correlations between drift velocities and wind from atmospheric reanalyses. This indicates a dominant physical control of wind forcing on ice drift and the persistence of free-drift conditions. These conditions occurred despite the buoy remaining largely in >80% ice concentrations and at distances >200 km from the estimated ice edge. The drift is additionally characterised by a strong inertial signature at 13.47 h, which appears initiated by passing cyclones. A wavelet analysis of the buoy's velocity confirms that the momentum transfer from winds at the multi-day frequencies is due to atmospheric forcing, while the initiation of inertial oscillations of sea ice has been identified as the secondary effect. Propagating storm-generated waves may initiate inertial oscillations by increasing the mobility of floes and enhance the drag of the inertial current. This analysis indicates that the Antarctic MIZ in the Indian Ocean sector remains much wider and mobile, during austral winter-to-spring, than defined by sea-ice concentration.

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.

The influence of surface melt on the flow of Greenland's largest outlet glaciers remains poorly known and in situ observations are few. We use field observations to link surface meltwater forcing to glacier-wide diurnal velocity variations on East Greenland's Helheim Glacier over two summer melt seasons. We observe diurnal variations in glacier speed that peak ~6.5 h after daily maximum insolation and extend from the terminus region to the equilibrium line. Both the amplitude of the diurnal speed variation and its sensitivity to daily melt are largest at the glacier terminus and decrease up-glacier, suggesting that the magnitude of the response is controlled not only by melt input volume and temporal variability, but also by background effective pressure, which approaches zero at the terminus. Our results provide evidence that basal lubrication by meltwater drives diurnal velocity variations at Greenland's marine-terminating glaciers in a similar manner to alpine glaciers and Greenland's land-terminating outlet glaciers.

The recently-advancing Taku Glacier is excavating subglacial sediments at high rates over multi-decadal timescales. However, sediment redistribution over shorter timescales remains unquantified. We use a variety of methods to study subglacial and proglacial sediment redistribution on decadal, seasonal, and daily timescales to gain insight into sub- and proglacial landscape formation. Both excavation and deposition were observed from 2003 to 2015 (2.8 ± 0.9 m a−1 to +2.9 ± 0.9 m a−1). The observed patterns imply that a subglacial conduit has occupied the same site over the past decade. Outwash fans on the subaerial end moraine experience fluvial sediment reworking almost year-round, with net sediment gain in winter and net sediment loss in summer, and an overall mass gain between 2005 and 2015. We estimate that tens of meters of sediment still underlie the glacier terminus, sediments which can be remobilized during future activity. However, imminent retreat from the proglacial moraine will limit its sediment supply, leaving the moraine vulnerable to erosion by bordering rivers. Retreat into an over-deepened basin will leave the glacier vulnerable to increased frontal ablation and accelerating retreat.

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.

In the marginal ice zone, surface waves drive motion of sea ice floes. The motion of floes relative to each other can cause periodic collisions, and drives the formation of pancake sea ice. Additionally, the motion of floes relative to the water results in turbulence generation at the interface between the ice and ocean below. These are important processes for the formation and growth of pancakes, and likely contribute to wave energy loss. Models and laboratory studies have been used to describe these motions, but there have been no in situ observations of relative ice velocities in a natural wave field. Here, we use shipboard stereo video to measure wave motion and relative motion of pancake floes simultaneously. The relative velocities of pancake floes are typically small compared to wave orbital motion (i.e. floes mostly follow the wave orbits). We find that relative velocities are well-captured by existing phase-resolved models, and are only somewhat over-estimated by using bulk wave parameters. Under the conditions observed, estimates of wave energy loss from ice–ocean turbulence are much larger than from pancake collisions. Increased relative pancake floe velocities in steeper wave fields may then result in more wave attenuation by increasing ice–ocean shear.

Climate variability can complicate efforts to interpret any long-term glacier mass-balance trends due to anthropogenic warming. Here we examine the impact of climate variability on the seasonal mass-balance records of 14 glaciers throughout Norway, Sweden and Svalbard using dynamical adjustment, a statistical method that removes orthogonal patterns of variability shared between each mass-balance record and sea-level pressure or sea-surface temperature predictor fields. For each glacier, the two leading predictor patterns explain 27–81% of the winter mass-balance variability and 24–69% of the summer mass-balance variability. The spatial and temporal structure of these patterns indicates that accumulation variability for all of the glaciers is strongly related to the North Atlantic Oscillation (NAO), with the Atlantic Multidecadal Oscillation (AMO) also modulating accumulation variability for the northernmost glaciers. Given this result, predicting glacier change in the region may depend on NAO and AMO predictability. In the raw mass-balance records, the glaciers throughout southern Norway have significantly negative summer trends, whereas the glaciers located closer to the Arctic have negative winter trends. Removing the effects of climate variability suggests it can bias trends in mass-balance records that span a few decades, but its effects on most of the longer-term mass-balance trends are minimal.

The response of the atmospheric boundary layer (ABL) to subgrid-scale variations of sea ice properties and fracturing is poorly understood and not taken into account in mesoscale Numerical Weather Prediction (NWP) model parametrizations. In this paper we analyze three-dimensional air circulation within the ABL over fragmented sea ice. A series of idealized high-resolution simulations with the Weather Research and Forecasting (WRF) model is performed for several spatial distributions of ice floes and leads for two values of sea ice concentration (0.5 and 0.9) and several ambient wind speed profiles. The results show that the convective circulation within the ABL is sensitive to the subgrid-scale spatial distribution of sea ice. Considerable variability of several domain-averaged quantities – cloud liquid water content, surface turbulent heat flux (THF) – is found for different arrangements of floes. Moreover, the organized structure of air circulation leads to spatial covariance of variables characterizing the ABL. Based on the example of THF, it is demonstrated that this covariance may lead to substantial errors when THF values are estimated from area-averaged quantities, as it is done in mesoscale NWP models. This suggests the need for developing suitable parametrizations of ABL effects related to subgrid-scale sea ice features for these models.

Significant salinity anomalies have been observed in the Arctic Ocean surface layer during the last decade. Our study is based on an extensive gridded dataset of winter salinity in the upper 50 m layer of the Arctic Ocean for the periods 1950–1993 and 2007–2012, obtained from ~20 000 profiles. We investigate the interannual variability of the salinity fields, identify predominant patterns of anomalous behavior and leading modes of variability, and develop a statistical model for the prediction of surface-layer salinity. The statistical model is based on linear regression equations linking the principal components of surface-layer salinity obtained through empirical orthogonal function decomposition with environmental factors, such as atmospheric circulation, river runoff, ice processes and water exchange with neighboring oceans. Using this model, we obtain prognostic fields of the surface-layer salinity for the winter period 2013–2014. The prognostic fields generated by the model show tendencies of surface-layer salinification, which were also observed in previous years. Although the used data are proprietary and have gaps, they provide the most spatiotemporally detailed observational resource for studying multidecadal variations in basin-wide Arctic salinity. Thus, there is community value in the identification, dissemination and modeling of the principal modes of variability in this salinity record.

Large discrepancies have been observed between satellite-derived sea-ice concentrations(IC) from passive microwave remote sensing and those derived from optical images at several locations in the East Antarctic, between February and April 2014. These artefacts, that resemble polynyas in the IC maps, appear in areas where optical satellite data show that there is landfast sea ice. The IC datasets and the corresponding retrieval algorithms are investigated together with microwave brightness temperature, air temperature, snowfall and bathymetry to understand the failure of the IC retrieval. The artefacts are the result of the application of weather filters in retrieval algorithms. These filters use the 37 and 19 GHz channels to correct for atmospheric effects on the retrieval. These channels show significant departures from typical ranges when the artefacts occur. A melt–refreeze cycle with associated snow metamorphism is proposed as the most likely cause. Together, the areas of the artefacts account for up to 0.5% of the Antarctic sea-ice area and thus cause a bias in sea-IC time series. In addition, erroneous sea ICs can adversely affect shipping operations.

Accelerated discharge through marine-terminating outlet glaciers has been a key component of the rapid mass loss from Arctic glaciers since the 1990s. However, glacier retreat and its climatic controls have not been assessed at the pan-Arctic scale. Consequently, the spatial and temporal variability in the magnitude of retreat, and the possible drivers are uncertain. Here we use remotely sensed data acquired over 273 outlet glaciers, located across the entire Atlantic Arctic (i.e. areas potentially influenced by North Atlantic climate and/or ocean conditions, specifically: Greenland, Novaya Zemlya, Franz Josef Land and Svalbard), to demonstrate high-magnitude, accelerating and near-ubiquitous retreat between 1992 and 2010. Overall, mean retreat rates increased by a factor of 3.5 between 1992 and 2000 (−30.5 m a−1) and 2000–10 (−105.8 m a−1), with 97% of the study glaciers retreating during the latter period. The Retreat was greatest in northern, western and south-eastern Greenland and also increased substantially on the Barents Sea coast of Novaya Zemlya. Glacier retreat showed no significant or consistent relationship with summer air temperatures at decadal timescales. The rate of frontal position change showed a significant, but weak, correlation with changes in sea-ice concentrations. We highlight large variations in retreat rates within regions and suggest that fjord topography plays an important role. We conclude that marine-terminating Arctic outlet glaciers show a common response of rapid and accelerating retreat at decadal timescales.

The Arctic marine environment is undergoing a transition from thick multi-year to first-year sea-ice cover with coincident lengthening of the melt season. Such changes are evident in the Baffin Bay-Davis Strait-Labrador Sea (BDL) region where melt onset has occurred ~8 days decade−1 earlier from 1979 to 2015. A series of anomalously early events has occurred since the mid-1990s, overlapping a period of increased upper-air ridging across Greenland and the northwestern North Atlantic. We investigate an extreme early melt event observed in spring 2013. (~6σ below the 1981–2010 melt climatology), with respect to preceding sub-seasonal mid-tropospheric circulation conditions as described by a daily Greenland Blocking Index (GBI). The 40-days prior to the 2013 BDL melt onset are characterized by a persistent, strong 500 hPa anticyclone over the region (GBI >+1 on >75% of days). This circulation pattern advected warm air from northeastern Canada and the northwestern Atlantic poleward onto the thin, first-year sea ice and caused melt ~50 days earlier than normal. The episodic increase in the ridging atmospheric pattern near western Greenland as in 2013, exemplified by large positive GBI values, is an important recent process impacting the atmospheric circulation over a North Atlantic cryosphere undergoing accelerated regional climate change.