Search

The Ross Ice Shelf (RIS) is currently stable but recent observations have indicated that basal melt rates beneath the ice shelf are expected to increase. It is important to know which areas of the RIS are more sensitive to enhanced basal melting as well as other external forcings or internal material properties of the ice to understand how climate change will influence RIS mass balance. In this paper, we use automatic differentiation and the Ice Sheet and Sea-level System Model to quantify the sensitivity of the RIS to changes in basal friction, ice rigidity, surface mass balance, and basal melting. Using volume above flotation (VAF) as our quantity of interest, we find that the RIS is most sensitive to changes in basal friction and ice rigidity close to grounding lines and along shear margins of the Siple Coast Ice Streams and Transantarctic Mountains Outlet Glaciers. The RIS sensitivity to surface mass balance is uniform over grounded ice, while the sensitivity to basal melting is more spatially variable. Changes in basal melting close to the grounding lines of the Siple Coast Ice Streams and Transantarctic Mountains outlet glaciers have a larger impact on the final VAF compared to elsewhere. Additionally, the pinning points and ice shelf shear margins are highly sensitive to changes in basal melt. Our sensitivity maps allow areas of greatest future vulnerability to be identified. [ABSTRACT FROM AUTHOR]

Energy exchanges between the atmosphere and the glacier surface control the net energy available for snow and ice melt. This paper explores the response of a midlatitude glacier in the Canadian Rocky Mountains to daily and interannual variations in the meteorological parameters that govern the surface energy balance. We use an energy balance model to run sensitivity tests to perturbations in temperature, specific humidity, wind speed, incoming shortwave radiation, glacier surface albedo, and winter snowpack depth. Variables are perturbed (i) in isolation, (ii) including internal feedbacks, and (iii) with co-evolution of meteorological perturbations, derived from the North American regional climate reanalysis (NARR) over the period 1979-2014. Summer melt at this site has the strongest sensitivity to interannual variations in temperature, albedo, and specific humidity, while fluctuations in cloud cover, wind speed, and winter snowpack depth have less influence. Feedbacks to temperature forcing, in particular summer albedo evolution, double the melt sensitivity to a temperature change. When meteorological perturbations covary through the NARR forcing, summer temperature anomalies remain important in driving interannual summer energy balance and melt variability, but they are reduced in importance relative to an isolated temperature forcing. Covariation of other variables (e.g., clear skies, giving reduced incoming longwave radiation) may be partially compensating for the increase in temperature. The methods introduced in this paper provide a framework that can be extended to compare the sensitivity of glaciers in different climate regimes, e.g., polar, maritime, or tropical environments, and to assess the importance of different meteorological parameters in different regions. [ABSTRACT FROM AUTHOR]

Due to recent atmospheric and oceanic warming, the Antarctic Peninsula is one of the most challenging regions of Antarctica to understand in terms of both local- and regional-scale climate signals. Steep topography and a lack of long-term and in situ meteorological observations complicate the extrapolation of existing climate models to the sub-regional scale. Therefore, new techniques must be developed to better understand processes operating in the region. Isotope signals are traditionally related mainly to atmospheric conditions, but a detailed analysis of individual components can give new insight into oceanic and atmospheric processes. This paper aims to use new isotopic records collected from snow and firn cores in conjunction with existing meteorological and oceanic datasets to determine changes at the climatic scale in the northern extent of the Antarctic Peninsula. In particular, a discernible effect of sea ice cover on local temperatures and the expression of climatic modes, especially the Southern Annular Mode (SAM), is demonstrated. In years with a large sea ice extension in winter (negative SAM anomaly), an inversion layer in the lower troposphere develops at the coastal zone. Therefore, an isotope-temperature relationship (δ-T) valid for all periods cannot be obtained, and instead the δ-T depends on the seasonal variability of oceanic conditions. Comparatively, transitional seasons (autumn and spring) have a consistent isotope-temperature gradient of +0.69 ‴ ϐC-1. As shown by firn core analysis, the near-surface temperature in the northern-most portion of the Antarctic Peninsula shows a decreasing trend (-0.33 ϐC year-1) between 2008 and 2014. In addition, the deuterium excess (dexcess) is demonstrated to be a reliable indicator of seasonal oceanic conditions, and therefore suitable to improve a firn age model based on seasonal dexcess variability. The annual accumulation rate in this region is highly variable, ranging between 1060 and 2470 kg m2 year-1 from 2008 to 2014. The combination of isotopic and meteorological data in areas where data exist is key to reconstruct climatic conditions with a high temporal resolution in polar regions where no direct observations exist. [ABSTRACT FROM AUTHOR]

The mountain cryosphere of mainland Europe is recognized to have important impacts on a range of environmental processes. In this paper, we provide an overview on the current knowledge on snow, glacier, and permafrost processes, as well as their past, current, and future evolution. We additionally provide an assessment of current cryosphere research in Europe and point to the different domains requiring further research. Emphasis is given to our understanding of climate-cryosphere interactions, cryosphere controls on physical and biological mountain systems, and related impacts. By the end of the century, Europe's mountain cryosphere will have changed to an extent that will impact the landscape, the hydrological regimes, the water resources, and the infrastructure. The impacts will not remain confined to the mountain area but also affect the downstream lowlands, entailing a wide range of socioeconomical consequences. European mountains will have a completely different visual appearance, in which low- and mid-range-altitude glaciers will have disappeared and even large valley glaciers will have experienced significant retreat and mass loss. Due to increased air temperatures and related shifts from solid to liquid precipitation, seasonal snow lines will be found at much higher altitudes, and the snow season will be much shorter than today. These changes in snow and ice melt will cause a shift in the timing of discharge maxima, as well as a transition of runoff regimes from glacial to nival and from nival to pluvial. This will entail significant impacts on the seasonality of high-altitude water availability, with consequences for water storage and management in reservoirs for drinking water, irrigation, and hydropower production. Whereas an upward shift of the tree line and expansion of vegetation can be expected into current periglacial areas, the disappearance of permafrost at lower altitudes and its warming at higher elevations will likely result in mass movements and process chains beyond historical experience. Future cryospheric research has the responsibility not only to foster awareness of these expected changes and to develop targeted strategies to precisely quantify their magnitude and rate of occurrence but also to help in the development of approaches to adapt to these changes and to mitigate their consequences. Major joint efforts are required in the domain of cryospheric monitoring, which will require coordination in terms of data availability and quality. In particular, we recognize the quantification of high-altitude precipitation as a key source of uncertainty in projections of future changes. Improvements in numerical modeling and a better understanding of process chains affecting high-altitude mass movements are the two further fields that - in our view - future cryospheric research should focus on. [ABSTRACT FROM AUTHOR]

Characterization of permafrost, particularly warm and near-surface permafrost which can contain significant liquid water, is critical to understanding complex interrelationships with climate change, ecosystems, and disturbances such as wildfires. Understanding the vulnerability and resilience of permafrost requires an interdisciplinary approach, relying on (for example) geophysical investigations, ecological characterization, direct observations, remote sensing, and more. As part of a multiyear investigation into the impacts of wildfires on permafrost, we have collected in situ measurements of the nuclear magnetic resonance (NMR) response of the active layer and permafrost in a variety of soil conditions, types, and saturations. In this paper, we summarize the NMR data and present quantitative relationships between active layer and permafrost liquid water content and pore sizes and show the efficacy of borehole NMR (bNMR) to permafrost studies. Through statistical analyses and synthetic freezing simulations, we also demonstrate that borehole NMR is sensitive to the nucleation of ice within soil pore spaces. [ABSTRACT FROM AUTHOR]

The impacts of the accelerated glacier retreat in recent decades on glacier runoff changes are still unknown in most Andean catchments, increasing uncertainties in estimating water availability. This particularly affects the outer tropics and Dry Andes, heavily impacted by prolonged droughts. Current global estimates overlook climatic and morphometric disparities, which significantly influence model parameters, among Andean glaciers. Meanwhile, local studies have used different approaches to estimate glacier runoff in a few catchments. Improving 21st-century glacier runoff projections relies on calibrating and validating models using corrected historical climate inputs and calibrated parameters across diverse glaciological zones. Here, we simulate glacier evolution and related runoff changes between the periods 2000–2009 and 2010–2019 across 786 Andean catchments (11 282 km2 of glacierized area, 11° N to 55° S) using the Open Global Glacier Model (OGGM). TerraClimate atmospheric variables were corrected using in situ data, getting a mean temperature bias by up to 2.1 °C and enhanced monthly precipitation. Glacier mass balance and volume were calibrated, where melt factor and the Glen A parameter exhibited significant alignment with varying environmental conditions. Simulation outcomes were validated against in situ data in three documented catchments (with a glacierized area > 8 %) and monitored glaciers. Our results at the Andes scale reveal an average reduction of 8.3 % in glacier volume and a decrease of 2.2 % in surface area between the periods 2000–2009 and 2010–2019. Comparing these two periods, glacier and climate variations have led to a 12 % increase in mean annual glacier melt (86.5 m3s-1) and a decrease in rainfall on glaciers of - 2 % (- 7.6 m3s-1) across the Andes, with both variables comprising the glacier runoff. We confirmed the utility of our corrected regional simulations of glacier runoff contribution at the catchment scale, where our estimations align with previous studies (e.g., Maipo 34° S, Chile) as well as provide new insights on the seasonal glaciers' largest contribution (e.g., La Paz 16° S, Bolivia) and new estimates of glacier runoff contribution (e.g., Baker 47° S, Chile). [ABSTRACT FROM AUTHOR]

A detailed understanding of snow cover and its possible feedback on climate change on the Tibetan Plateau (TP) is of great importance. However, spatiotemporal variability in snow phenology (SP) and its influencing factors on the TP remain unclear. Based on the daily gap-free snow cover product (HMRFS-TP) with 500 m resolution, this study investigated the spatiotemporal variability in snow cover days (SCDs), snow onset date (SOD), and snow end date (SED) on the TP from 2002 to 2022. A structural equation model was used to quantify the direct and indirect effects of meteorological factors, geographical location, topography, and vegetation greenness on SP. The results indicate that the spatial distribution of SP on the TP was extremely uneven and exhibited temporal heterogeneity. SP showed vertical zonality influenced by elevation (longer SCD, earlier SOD, and later SED at higher elevations). A total of 4.62 % of the TP area had a significant decrease in SCDs, at a rate of -1.74 d yr -1. The SOD of 2.34 % of the TP area showed a significant delayed trend, at a rate of 2.90 d yr -1 , while the SED of 1.52 % of the TP area had a significant advanced trend, at a rate of at -2.49 d yr -1. We also found a strong elevation dependence for the trend in SCDs (R=-0.73). Air temperature, precipitation, wind speed, and shortwave radiation can directly affect SP as well as indirectly affect it by influencing the growth of vegetation, whereas the direct effect was much greater than the indirect effect. Geographical location (latitude and longitude) and topographic conditions (elevation and slope) indirectly affected SP by modulating meteorological conditions and the growth of vegetation. Vegetation primarily influences SP by intercepting the snow and regulating the balance of the solar radiation budget. Regarding the total effect, air temperature was found to be the dominant factor. This study contributes to the understanding of snow variation in response to global warming over the past 2 decades by providing a basis for predicting future environmental and climate changes and their impacts on the TP. [ABSTRACT FROM AUTHOR]

The Torngat Mountains National Park, northern Labrador, Canada, contains more than 120 small glaciers: the only remaining glaciers in continental northeast North America. These small cirque glaciers exist in a unique topoclimatic setting, experiencing temperate maritime summer conditions yet very cold and dry winters, and may provide insights into the deglaciation dynamics of similar small glaciers in temperate mountain settings. Due to their size and remote location, very little information exists regarding the health of these glaciers. Just a single study has been published on the contemporary glaciology of the Torngat Mountains, focusing on net mass balances from 1981 to 1984. This paper addresses the extent to which glaciologically relevant climate variables have changed in northern Labrador in concert with 20th-century Arctic warming, and how these changes have affected Torngat Mountain glaciers. Field surveys and remote-sensing analyses were used to measure regional glacier area loss of 27% from 1950 to 2005, substantial rates of ice surface thinning (up to 6 m yr-1) and volume losses at Abraham, Hidden, and Minaret glaciers, between 2005 and 2011. Glacier mass balances appear to be controlled by variations in winter precipitation and, increasingly, by strong summer and autumn atmospheric warming since the early 1990s, though further observations are required to fully understand mass balance sensitivities. This study provides the first comprehensive contemporary assessment of Labrador glaciers and will inform both regional impact assessments and syntheses of global glacier mass balance. [ABSTRACT FROM AUTHOR]

Sea ice is crucial in regulating the heat balance between the ocean and atmosphere and quintessential for supporting the prevailing Arctic food web. Due to limited and often local data availability back in time, the sensitivity of sea-ice proxies to long-term climate changes is not well constrained, which renders any comparison with palaeoclimate model simulations difficult. Here we compiled a set of marine sea-ice proxy records with a relatively high temporal resolution of at least 100 years, covering the Common Era (past 2k years) in the Greenland–North Atlantic sector of the Arctic to explore the presence of coherent long-term trends and common low-frequency variability, and we compared those data with transient climate model simulations. We used cluster analysis and empirical orthogonal functions to extract leading modes of sea-ice variability, which efficiently filtered out local variations and improved comparison between proxy records and model simulations. We find that a compilation of multiple proxy-based sea-ice reconstructions accurately reflects general long-term changes in sea-ice history, consistent with simulations from two transient climate models. Although sea-ice proxies have varying mechanistic relationships to sea-ice cover, typically differing in habitat or seasonal representation, the long-term trend recorded by proxy-based reconstructions showed a good agreement with summer minimum sea-ice area from the model simulations. The short-term variability was not as coherent between proxy-based reconstructions and model simulations. The leading mode of simulated sea ice associated with the multidecadal to centennial timescale presented a relatively low explained variance and might be explained by changes in solar radiation and/or inflow of warm Atlantic waters to the Arctic Ocean. Short variations in proxy-based reconstructions, however, are mainly associated with local factors and the ecological nature of the proxies. Therefore, a regional or large-scale view of sea-ice trends necessitates multiple spatially spread sea-ice proxy-based reconstructions, avoiding confusion between long-term regional trends and short-term local variability. Local-scale sea-ice studies, in turn, benefit from reconstructions from well-understood individual research sites. [ABSTRACT FROM AUTHOR]

The Svalbard archipelago is particularly sensitive to climate change due to the relatively low altitude of its main ice fields and its geographical location in the higher North Atlantic, where the effect of Arctic amplification is more significant. The largest temperature increases have been observed during winter, but increasing summer temperatures, above the melting point, have led to increased glacier melt. Here, we evaluate the impact of this increased melt on the preservation of the oxygen isotope (δ18 O) signal in firn records. δ18 O is commonly used as a proxy for past atmospheric temperature reconstructions, and, when preserved, it is a crucial parameter to date and align ice cores. By comparing four different firn cores collected in 2012, 2015, 2017 and 2019 at the top of the Holtedahlfonna ice field (1100 m a.s.l.), we show a progressive deterioration of the isotope signal, and we link its degradation to the increased occurrence and intensity of melt events. Our findings indicate that, starting from 2015, there has been an escalation in melting and percolation resulting from changes in the overall atmospheric conditions. This has led to the deterioration of the climate signal preserved within the firn or ice. Our observations correspond with the model's calculations, demonstrating an increase in water percolation since 2014, potentially reaching deeper layers of the firn. Although the δ18 O signal still reflects the interannual temperature trend, more frequent melting events may in the future affect the interpretation of the isotopic signal, compromising the use of Svalbard ice cores. Our findings highlight the impact and the speed at which Arctic amplification is affecting Svalbard's cryosphere. [ABSTRACT FROM AUTHOR]

Glaciological and geodetic methods provide independent observations of glacier mass balance. The glaciological method measures the surface mass balance, on a seasonal or annual basis, whereas the geodetic method measures surface, internal, and basal mass balances, over a period of years or decades. In this paper, we reanalyse the 10 glaciers with long-term mass-balance series in Norway. The reanalysis includes (i) homogenisation of both glaciological and geodetic observation series, (ii) uncertainty assessment, (iii) estimates of generic differences including estimates of internal and basal melt, (iv) validation, and, if needed, (v) calibration of mass-balance series. This study comprises an extensive set of data (484 mass-balance years, 34 geodetic surveys, and large volumes of supporting data, such as metadata and field notes). In total, 21 periods of data were compared and the results show discrepancies between the glaciological and geodetic methods for some glaciers, which are attributed in part to internal and basal ablation and in part to inhomogeneity in the data processing. Deviations were smaller than 0.2 m w.e.a-1 for 12 out of 21 periods. Calibration was applied to 7 out of 21 periods, as the deviations were larger than the uncertainty. The reanalysed glaciological series shows a more consistent signal of glacier change over the period of observations than previously reported: six glaciers had a significant mass loss (14-22 m w.e.) and four glaciers were nearly in balance. All glaciers have lost mass after the year 2000. More research is needed on the sources of uncertainty to reduce uncertainties and adjust the observation programmes accordingly. The study confirms the value of carrying out independent high-quality geodetic surveys to check and correct field observations. [ABSTRACT FROM AUTHOR]

Sea ice concentration has been retrieved in polar regions with satellite microwave radiometers for over 30 years. However, the question remains as to what is an optimal sea ice concentration retrieval method for climate monitoring. This paper presents some of the key results of an extensive algorithm inter-comparison and evaluation experiment. The skills of 30 sea ice algorithms were evaluated sysematically over low and high sea ice concentrations. Evaluation criteria included standard deviation relative to independent validation data, performance in the presence of thin ice and melt ponds, and sensitivity to error sources with seasonal to inter-annual variations and potential climatic trends, such as atmospheric water vapour and water-surface roughening by wind. A selection of 13 algorithms is shown in the article to demonstrate the results. Based on the findings, a hybrid approach is suggested to retrieve sea ice concentration globally for climate monitoring purposes. This approach consists of a combination of two algorithms plus dynamic tie points implementation and atmospheric correction of input brightness temperatures. The method minimizes inter-sensor calibration discrepancies and sensitivity to the mentioned error sources. [ABSTRACT FROM AUTHOR]

With the EU-funded PACE (Permafrost and Climate in Europe) project at the turn of this century, several deep boreholes (100 m +) were drilled in European mountain sites, including in mainland Norway, Svalbard and Sweden. During other projects from 2004 and the International Polar Year (IPY) period in 2006–2007, several additional boreholes were drilled in different sites in both Norway and Iceland, measuring temperatures along both altitudinal and latitudinal gradients. At most sites, multi-temporal geophysical soundings are available using electrical resistivity tomography (ERT). Here, we study the development of permafrost and ground temperatures in mainland Norway and Iceland based on these data sets. We document that permafrost in Norway and Iceland is warming at a high rate, including the development of taliks in both Norway and Iceland in response to global climate change during the last 20 years. At most sites, ground surface temperature (GST) is apparently increasing more strongly than surface air temperature (SAT). Changing snow conditions appear to be the most important factor for the higher GST rates. Modelling exercises also indicate that the talik development can be explained by both higher air temperatures and increasing snow depth. [ABSTRACT FROM AUTHOR]

Stable water isotopologues of snow, firn and ice cores provide valuable information on past climate variations. Yet single profiles are generally not suitable for robust climate reconstructions. Stratigraphic noise, introduced by the irregular deposition, wind-driven erosion and redistribution of snow, impacts the utility of high-resolution isotope records, especially in low-accumulation areas. However, it is currently unknown how stratigraphic noise differs across the East Antarctic Plateau and how it is affected by local environmental conditions. Here, we assess the amount and structure of stratigraphic noise at seven sites along a 120 km transect on the plateau of Dronning Maud Land, East Antarctica. Replicated oxygen isotope records of 1 m length were used to estimate signal-to-noise ratios as a measure of stratigraphic noise at sites characterised by different accumulation rates (43–64 mm w.e. a -1), snow surface roughnesses and slope inclinations. While we found a high level of stratigraphic noise at all sites, there was also considerable variation between sites. At sastrugi-dominated sites, greater stratigraphic noise coincided with stronger surface roughnesses, steeper slopes and lower accumulation rates, probably related to increased wind speeds. These results provide a first step to modelling stratigraphic noise and might guide site selection and sampling strategies for future expeditions to improve high-resolution climate reconstructions from low-accumulation regions. [ABSTRACT FROM AUTHOR]

The ground ice content in cold environments influences the permafrost thermal regime and the thaw trajectories in a warming climate, especially for soils containing excess ice. Despite their importance, the amount and distribution of ground ice are often unknown due to lacking field observations. Hence, modeling the thawing of ice-rich permafrost soils and associated thermokarst is challenging as ground ice content has to be prescribed in the model setup. In this study, we present a model scheme, capable of simulating segregated ice formation during a model spinup together with associated ground heave. It provides the option to add a constant sedimentation rate throughout the simulation. Besides ice segregation, it can represent thaw consolidation processes and ground subsidence under a warming climate. The computation is based on soil mechanical processes, soil hydrology by the Richards equation and soil freezing characteristics. The code is implemented in the CryoGrid community model (version 1.0), a modular land surface model for simulations of the ground thermal regime. The simulation of ice segregation and thaw consolidation with the new model scheme allows us to analyze the evolution of ground ice content in both space and time. To do so, we use climate data from two contrasting permafrost sites to run the simulations. Several influencing factors are identified, which control the formation and thaw of segregated ice. (i) Model results show that high temperature gradients in the soil as well as moist conditions support the formation of segregated ice. (ii) We find that ice segregation increases in fine-grained soils and that especially organic-rich sediments enhance the process. (iii) Applying external loads suppresses ice segregation and speeds up thaw consolidation. (iv) Sedimentation leads to a rise of the ground surface and the formation of an ice-enriched layer whose thickness increases with sedimentation time. We conclude that the new model scheme is a step forward to improve the description of ground ice distributions in permafrost models and can contribute towards the understanding of ice segregation and thaw consolidation in permafrost environments under changing climatic conditions. [ABSTRACT FROM AUTHOR]

Volume-area scaling is the most popular method for estimating the ice volume of large glacier samples. Here, a series of resampling experiments based on different sets of synthetic data is presented in order to derive an upperbound estimate (i.e. a level achieved only within ideal conditions) for its accuracy. For real-world applications, a lower accuracy has to be expected. We also quantify the maximum accuracy expected when scaling is used for determining the glacier volume change, and area change of a given glacier population. A comprehensive set of measured glacier areas, volumes, area and volume changes is evaluated to investigate the impact of real-world data quality on the so-assessed accuracies. For populations larger than a few thousand glaciers, the total ice volume can be recovered within 30% if all data currently available worldwide are used for estimating the scaling parameters. Assuming no systematic bias in ice volume measurements, their uncertainty is of secondary importance. Knowing the individual areas of a glacier sample for two points in time allows recovering the corresponding ice volume change within 40% for populations larger than a few hundred glaciers, both for steady-state and transient geometries. If ice volume changes can be estimated without bias, glacier area changes derived from volume-area scaling show similar uncertainties to those of the volume changes. This paper does not aim at making a final judgement on the suitability of volume-area scaling as such, but provides the means for assessing the accuracy expected from its application. [ABSTRACT FROM AUTHOR]

Earth system models (ESMs) can help to improve the understanding of climate-induced cryospheric–hydrological impacts in complex mountain regions, such as High Mountain Asia (HMA). Coarse ESM grids, however, have difficulties in representing cryospheric–hydrological processes that vary over short distances in complex mountainous environments. Variable-resolution (VR) ESMs can help to overcome these limitations through targeted grid refinement. This study investigates the ability of the VR Community Earth System Model (VR-CESM) to simulate cryospheric–hydrological variables such as the glacier surface mass balance (SMB) over HMA. To this end, a new VR grid is generated, with a regional grid refinement up to 7 km over HMA. Two coupled atmosphere–land simulations are run for the period 1979–1998. The second simulation is performed with an updated glacier cover dataset and includes snow and glacier model modifications. Comparisons are made to gridded outputs derived from a globally uniform 1 ∘ CESM grid, observation-, reanalysis-, and satellite-based datasets, and a glacier model forced by a regional climate model (RCM). Climatological biases are generally reduced compared to the coarse-resolution CESM grid, but the glacier SMB is too negative relative to observation-based glaciological and geodetic mass balances, as well as the RCM-forced glacier model output. In the second simulation, the SMB is improved but is still underestimated due to cloud cover and temperature biases, missing model physics, and incomplete land–atmosphere coupling. The outcomes suggest that VR-CESM could be a useful tool to simulate cryospheric–hydrological variables and to study climate change in mountainous environments, but further developments are needed to better simulate the SMB of mountain glaciers. [ABSTRACT FROM AUTHOR]

Winter warm air intrusions entering the Arctic region can strongly modify the microwave emission of the snow-covered sea ice system due to temperature-induced snow metamorphism and ice crust formations, e.g., after melt–refreeze events. Common microwave radiometer satellite sea ice concentration retrievals are based on empirical models using the snow-covered sea ice emissivity and thus can be influenced by strong warm air intrusions. Here, we carry out a long-term study analyzing 41 years of winter sea ice concentration observations from different algorithms to investigate the impact of warm air intrusions on the retrieved ice concentration. Our results show that three out of four algorithms underestimate the sea ice concentration during (and up to 10 d after) warm air intrusions which increase the 2 m air temperature (daily maximum) above - 5 ∘ C. This can lead to sea ice area underestimations in the order of 10 4 to 10 5 km 2. If the 2 m temperature during the warm air intrusions crosses - 2 ∘ C, all algorithms are impacted. Our analysis shows that the strength of these strong warm air intrusions increased in recent years, especially in April. With a further climate change, such warm air intrusions are expected to occur more frequently and earlier in the season, and their influence on sea ice climate data records will become more important. [ABSTRACT FROM AUTHOR]

Anthropogenic climate change threatens northern permafrost environments. This compromises the existence of permafrost landforms, such as palsas and peat plateaus, which have been assessed to be critically endangered habitats. In this study, we integrated geospatial datasets and statistical methods to model the suitable environments for palsas and peat plateaus across the Northern Hemisphere permafrost region. The models were calibrated using data from years 1950–2000. The effects of climate change on the suitable environments for the landforms were assessed by using low-, moderate-, and high-emissions scenarios (Representative Concentration Pathway climate scenarios: RCP2.6, RCP4.5, and RCP8.5, respectively) for two periods (2041–2060 and 2061–2080). Hotspots for palsa and peat plateau environments occurred in northern Europe, western Siberia, and subarctic Canada. Climate change was predicted to cause an almost complete loss (decrease of 98.2 %) of suitable environmental spaces under the high-emissions scenario by 2061–2080, while under low- and moderate-emissions scenarios the predicted loss was 76.3 % and 89.3 % respectively. Our modeling results are in line with previously published thermokarst data pointing out areas of recent degradation of palsa and peat plateau environments. Our results provide new insights into the distribution of the permafrost landforms in less studied areas such as central and eastern Siberia. In addition, the predictions provide new understanding of the changing geoecological conditions of the circumpolar region with important implications for greenhouse gas emissions. [ABSTRACT FROM AUTHOR]

We have investigated the drivers behind the observed spatiotemporal mass balance heterogeneity in Tien Shan and Pamir, in High Mountain Asia. To study the consistency of the different interpretations derived from the available meteorological reanalysis and remote sensing products, we used correlation analyses between climatic and static drivers with novel estimates of region-wide annual glacier mass balance time series. These analyses were performed both spatially using different spatial classifications of glaciers and temporally for each individual glacier. Our results show that the importance of the variables studied depends strongly on the dataset used and which spatial classification of glaciers is chosen. This extends to opposing results using the different products. Even supposedly similar datasets lead to different and partly contradicting assumptions on dominant drivers of mass balance variability. The apparent but false consistencies across studies using a single dataset are related, according to our results, to the chosen dataset or spatial classification rather than to the processes or involved environmental variables. Without a glaciological, meteorological, and hydrological in situ observation network providing data that allow for the direct calibration and validation of extensive datasets, our understanding of neither the changing cryosphere at the regional scale for Tien Shan and Pamir nor glacier response to climate change or the assessment of water availability for the region's growing population can improve. [ABSTRACT FROM AUTHOR]

Areas of thin sea ice in the polar regions not only are experiencing the highest rate of sea-ice production but also are, therefore, important hot spots for ocean ventilation as well as heat and moisture exchange between the ocean and the atmosphere. Through co-location of (1) an unsupervised waveform classification (UWC) approach applied to CryoSat-2 radar waveforms with (2) Moderate Resolution Imaging Spectroradiometer-derived (MODIS) thin-ice-thickness estimates and (3) Sentinel-1A/B synthetic-aperture radar (SAR) reference data, thin-ice-based waveform shapes are identified, referenced, and discussed with regard to a manifold of waveform shape parameters. Here, strong linear dependencies are found between binned thin-ice thickness up to 25 cm from MODIS and the CryoSat-2 waveform shape parameters that show the possibility of either developing simple correction terms for altimeter ranges over thin ice or directing adjustments to current retracker algorithms specifically for very thin sea ice. This highlights the potential of CryoSat-2-based SAR altimetry to reliably discriminate between occurrences of thick sea ice, open-water leads, and thin ice within recently refrozen leads or areas of thin sea ice. Furthermore, a comparison to the ESA Climate Change Initiative's (CCI) CryoSat-2 surface type classification with classes sea ice, lead, and unknown reveals that the newly found thin-ice-related waveforms are divided up almost equally between unknown (46.3 %) and lead type (53.4 %) classifications. Overall, the UWC results in far fewer unknown classifications (1.4 % to 38.7 %). Thus, UWC provides more usable information for sea-ice freeboard and thickness retrieval and at the same time reduces range biases from thin-ice waveforms processed as regular sea ice in the CCI classification. [ABSTRACT FROM AUTHOR]

Arctic sea ice type (SITY) variation is a sensitive indicator of climate change. However, systematic inter-comparison and analysis for SITY products are lacking. This study analysed eight daily SITY products from five retrieval approaches covering the winters of 1999–2019, including purely radiometer-based (C3S-SITY), scatterometer-based (KNMI-SITY and IFREMER-SITY) and combined ones (OSISAF-SITY and Zhang-SITY). These SITY products were inter-compared against a weekly sea ice age product (i.e. NSIDC-SIA – National Snow and Ice Data Center sea ice age) and evaluated with five synthetic aperture radar (SAR) images. The average Arctic multiyear ice (MYI) extent difference between the SITY products and NSIDC-SIA varies from -1.32×106 to 0.49×106 km 2. Among them, KNMI-SITY and Zhang-SITY in the QuikSCAT (QSCAT) period (2002–2009) agree best with NSIDC-SIA and perform the best, with the smallest bias of -0.001×106 km 2 in first-year ice (FYI) extent and -0.02×106 km 2 in MYI extent. In the Advanced Scatterometer (ASCAT) period (2007–2019), KNMI-SITY tends to overestimate MYI (especially in early winter), whereas Zhang-SITY and IFREMER-SITY tend to underestimate MYI. C3S-SITY performs well in some early winter cases but exhibits large temporal variabilities like OSISAF-SITY. Factors that could impact performances of the SITY products are analysed and summarized. (1) The Ku-band scatterometer generally performs better than the C-band scatterometer for SITY discrimination, while the latter sometimes identifies FYI more accurately, especially when surface scattering dominates the backscatter signature. (2) A simple combination of scatterometer and radiometer data is not always beneficial without further rules of priority. (3) The representativeness of training data and efficiency of classification are crucial for SITY classification. Spatial and temporal variation in characteristic training datasets should be well accounted for in the SITY method. (4) Post-processing corrections play important roles and should be considered with caution. [ABSTRACT FROM AUTHOR]

The Totten Glacier in East Antarctica is of major climatic interest because of the large fluctuations in its grounding line and potential vulnerability to climate change. Here, we use a series of high-resolution, regional NEMO-LIM-based (Nucleus for European Modelling of the Ocean coupled with the Louvain-la-Neuve sea ice model) experiments, which include an explicit treatment of ocean–ice shelf interactions, as well as a representation of grounded icebergs and fast ice, to investigate the changes in ocean–ice interactions in the Totten Glacier area between the recent past (1995–2014) and the end of the 21st century (2081–2100) under SSP4–4.5 climate change conditions. By the end of the 21st century, the wide areas of multiyear fast ice simulated in the recent past are replaced by small patches of first year fast ice along the coast, which decreases the total summer sea ice extent. The Antarctic Slope Current is accelerated by about 116 %, which decreases the heat exchange across the shelf and tends to reduce the ice shelf basal melt rate, but this effect is counterbalanced by the effect of the oceanic warming. As a consequence, despite the accelerated Antarctic Slope Current, the Totten ice shelf melt rate is increased by 91 % due to the intrusion of warmer water into its cavity. The representation of fast ice dampens the ice shelf melt rate increase throughout the 21st century, as the Totten ice shelf melt rate increase reaches 136 % when fast ice is not taken into account. The Moscow University ice shelf melt rate increase is even more impacted by the representation of fast ice, with a 36 % melt rate increase with fast ice, compared to a 75 % increase without a fast ice representation. This influence of the representation of fast ice in our simulations on the basal melting rate trend over the 21st century is explained by the large impact of the fast ice for present-day conditions (∼25 % difference in m yr -1), while the impact decreases significantly at the end of the 21st century (∼4 % difference in m yr -1). As a consequence, the reduction in the fast ice extent in the future induces a decrease in the fast ice effect on the ice shelf melt rate that partly compensates for the increase due to warming of the ocean. This highlights the importance of including a representation of fast ice to simulate realistic ice shelf melt rate increase in East Antarctica under warming conditions. [ABSTRACT FROM AUTHOR]

As changes to Earth's polar climate accelerate, the need for robust, long–term sea ice thickness observation datasets for monitoring those changes and for verification of global climate models is clear. By linking an algorithm for retrieving snow–ice interface temperature from passive microwave satellite data to a thermodynamic sea ice energy balance relation known as Stefan's law, we have developed a retrieval method for estimating thermodynamic sea ice thickness growth from space: Stefan's Law Integrated Conducted Energy (SLICE). With an initial condition at the beginning of the sea ice growth season, the method can model basin-wide absolute sea ice thickness by combining the one-dimensional SLICE retrieval with an ice motion dataset. The advantages of the SLICE retrieval method include daily basin-wide coverage, lack of atmospheric reanalysis product input requirement, and a potential for use beginning in 1987. Validation of the retrieval against measurements from 10 ice mass balance buoys shows a mean correlation of 0.89 and a mean bias of 0.06 m over the course of an entire sea ice growth season. Despite its simplifications and assumptions relative to models like the Pan-Arctic Ice–Ocean Modeling and Assimilation System (PIOMAS), basin-wide SLICE performs nearly as well as PIOMAS when compared against CryoSat-2 and Operation IceBridge using a linear correlation between collocated points. [ABSTRACT FROM AUTHOR]

Ice aprons (IAs) are part of the critical components of the Alpine cryosphere. As a result of the changing climate over the past few decades, deglaciation has resulted in a surface decrease of IAs, which has not yet been documented, except for a few specific examples. In this study, we quantify the effects of climate change on IAs since the mid-20th century in the Mont Blanc massif (western European Alps). We then evaluate the role of meteorological parameters and the local topography in the behaviour of IAs. We precisely mapped the surface areas of 200 IAs using high-resolution aerial and satellite photographs from 1952, 2001, 2012 and 2019. From the latter inventory, the surface area of the present individual IAs ranges from 0.001 to 0.04 km 2. IAs have lost their surface area over the past 70 years, with an alarming increase since the early 2000s. The total area, from 7.93 km 2 in 1952, was reduced to 5.91 km 2 in 2001 (-25.5 %) before collapsing to 4.21 km 2 in 2019 (-47 % since 1952). We performed a regression analysis using temperature and precipitation proxies to better understand the effects of meteorological parameters on IA surface area variations. We found a strong correlation between both proxies and the relative area loss of IAs, indicating the significant influence of the changing climate on the evolution of IAs. We also evaluated the role of the local topographic factors in the IA area loss. At a regional scale, factors like direct solar radiation and elevation influence the behaviour of IAs, while others like curvature, slope and size of the IAs seem to be rather important on a local scale. [ABSTRACT FROM AUTHOR]

Mountain glaciers provide us a window into past climate changes and landscape evolution, but the pattern of glacier evolution at centennial or suborbital timescale remains elusive, especially in monsoonal Himalayas. We simulated the glacier evolution in Bhutanese Himalaya (BH), a typical monsoon-influenced region, during the Little Ice Age (LIA) using the Open Global Glacier Model driven by six paleoclimate datasets and their average. Compared with geomorphologically mapped glacial landforms, the model can well capture the patterns of glacier length change. Simulation results revealed four glacial substages (the 1270s, 1470s, 1710s, and 1850s) during LIA in the study area. Statistically, a positive correlation between the number of glacial substages and glacier slope was found, indicating that the occurrence of glacial substages might be a result from heterogeneous responses of glaciers to climate change. Monthly climate change analysis and sensitivity experiments indicated that the summer temperature largely dominates the regional glacier evolution during the LIA in BH. [ABSTRACT FROM AUTHOR]

Climate projection studies of future changes in snow conditions and resulting rain-on-snow (ROS) flood events are subject to large uncertainties. Typically, emission scenario uncertainties and climate model uncertainties are included. This is the first study on this topic to also include quantification of natural climate variability, which is the dominant uncertainty for precipitation at local scales with large implications for runoff projections, for example. To quantify natural climate variability, a weather generator was applied to simulate inherently consistent climate variables for multiple realizations of current and future climates at 100 m spatial and hourly temporal resolution over a 12×12 km high-altitude study area in the Swiss Alps. The output of the weather generator was used as input for subsequent simulations with an energy balance snow model. The climate change signal for snow water resources stands out as early as mid-century from the noise originating from the three sources of uncertainty investigated, namely uncertainty in emission scenarios, uncertainty in climate models, and natural climate variability. For ROS events, a climate change signal toward more frequent and intense events was found for an RCP 8.5 scenario at high elevations at the end of the century, consistently with other studies. However, for ROS events with a substantial contribution of snowmelt to runoff (> 20 %), the climate change signal was largely masked by sources of uncertainty. Only those ROS events where snowmelt does not play an important role during the event will occur considerably more frequently in the future, while ROS events with substantial snowmelt contribution will mainly occur earlier in the year but not more frequently. There are two reasons for this: first, although it will rain more frequently in midwinter, the snowpack will typically still be too cold and dry and thus cannot contribute significantly to runoff; second, the very rapid decline in snowpack toward early summer, when conditions typically prevail for substantial contributions from snowmelt, will result in a large decrease in ROS events at that time of the year. Finally, natural climate variability is the primary source of uncertainty in projections of ROS metrics until the end of the century, contributing more than 70 % of the total uncertainty. These results imply that both the inclusion of natural climate variability and the use of a snow model, which includes a physically based process representation of water retention, are important for ROS projections at the local scale. [ABSTRACT FROM AUTHOR]

At high latitudes, long-term changes in riverine ice break-ups are exemplary measures of climatic change and variation. This study compares cryophenological trends, patterns and changes for the rivers Aura (1749–2020), Torne (1693–2020) and Kokemäki (1793–2020); all sites are located in Finland. The Kokemäki River series is a new series from the city of Pori. The findings show statistically significant cross-correlations between the Aura and Kokemäki rivers but weaker cross-correlations with the Torne River. We attribute the latter to climatic differences caused by the higher latitude of the Torne River. Taken together, the many results of this study suggest that the spring climate in the south has changed more rapidly and become less predictable than in the north. Climatic extremes – warmer and wetter winters – in the 2000s resulted in the first recorded no-freeze events on the Aura and Kokemäki rivers. This was the culmination of a rapid increase in early ice break-up events and interannual variability over the last 30 years. The number of early events has increased in all three rivers since the early or mid-1900s, but the earliest recorded break-up day on the Torne River has changed only marginally in the last 100 years. Our dynamic temperature analysis shows that the ice break-up on the Torne River requires higher temperatures than in the south, and future changes in the timing of the break-up depend on April temperatures. In the south, on the other hand, future changes concerns the return period of no-freeze events, which depend on temperature and precipitation during winter. [ABSTRACT FROM AUTHOR]

Ice shelves play a key role in the dynamics of marine ice sheets by buttressing grounded ice and limiting rates of ice flux to the oceans. In response to recent climatic and oceanic change, ice shelves fringing the West Antarctic Ice Sheet (WAIS) have begun to fragment and retreat, with major implications for ice-sheet stability. Here, we focus on the Thwaites Eastern Ice Shelf (TEIS), the remaining pinned floating extension of Thwaites Glacier. We show that TEIS has undergone a process of fragmentation in the last 5 years, including brittle failure along a major shear zone, formation of tensile cracks on the main body of the shelf, and a release of tabular bergs on both the eastern and western flanks. Simulations with the Helsinki Discrete Element Model (HiDEM) show that this pattern of failure is associated with high backstress from a submarine pinning point at the distal edge of the shelf. We show that a significant zone of shear, upstream of the main pinning point, developed in response to the rapid acceleration of the shelf between 2002 and 2006, seeding damage on the shelf. Subsequently, basal melting and positive feedback between damage and strain rates weakened TEIS, allowing damage to accumulate. Thus, although backstress on TEIS has likely diminished over time as the pinning point shrunk, accumulation of damage has ensured that the ice in the shear zone remained the weakest link in the system. Experiments with the BISICLES ice-sheet model indicate that additional damage to or unpinning of TEIS is unlikely to trigger significantly increased ice loss from WAIS, but the calving response to the loss of TEIS remains highly uncertain. It is widely recognised that ice-shelf fragmentation and collapse can be triggered by hydrofracturing and/or unpinning from ice-shelf margins or grounding points. Our results indicate a third mechanism, backstress triggered failure, that can occur if and when an ice shelf is no longer able to withstand stress imposed by pinning points. In most circumstances, pinning points are essential for ice-shelf stability, but as ice shelves thin and weaken, the concentration of backstress in damaged ice upstream of a pinning point may provide the seeds of their demise. [ABSTRACT FROM AUTHOR]

As Arctic sea ice extent continues to decline, remote sensing observations are becoming even more vital for the monitoring and understanding of sea ice. Recently, the sea ice community has entered a new era of synthetic aperture radar (SAR) satellites operating at C-band with the launch of Sentinel-1A in 2014 and Sentinel-1B (S1) in 2016 and the RADARSAT Constellation Mission (RCM) in 2019. These missions represent five spaceborne SAR sensors that together routinely cover the pan-Arctic sea ice domain. Here, we describe, apply, and validate the Environment and Climate Change Canada automated sea ice tracking system (ECCC-ASITS) that routinely generates large-scale sea ice motion (SIM) over the pan-Arctic domain using SAR images from S1 and RCM. We applied the ECCC-ASITS to the incoming image streams of S1 and RCM from March 2020 to October 2021 using a total of 135 471 SAR images and generated new SIM datasets (7 d 25 km and 3 d 6.25 km) by combining the image stream outputs of S1 and RCM (S1 + RCM). Results indicate that S1 + RCM SIM provides more coverage in the Hudson Bay, Davis Strait, Beaufort Sea, Bering Sea, and directly over the North Pole compared to SIM from S1 alone. Based on the resolvable S1 + RCM SIM grid cells, the 7 d 25 km spatiotemporal scale is able to provide the most complete picture of SIM across the pan-Arctic from SAR imagery alone, but considerable spatiotemporal coverage is also available from 3 d 6.25 SIM products. S1 + RCM SIM is resolved within the narrow channels and inlets of the Canadian Arctic Archipelago, filling a major gap from coarser-resolution sensors. Validating the ECCC-ASITS using S1 and RCM imagery against buoys indicates a root-mean-square error (RMSE) of 2.78 km for dry ice conditions and 3.43 km for melt season conditions. Higher speeds are more apparent with S1 + RCM SIM as comparison with the National Snow and Ice Data Center (NSIDC) SIM product and the Ocean and Sea Ice Satellite Application Facility (OSI SAF) SIM product indicated an RMSE of u=4.6 km d -1 and v=4.7 km d -1 for the NSIDC and u=3.9 km d -1 and v=3.9 km d -1 for OSI SAF. Overall, our results demonstrate the robustness of the ECCC-ASITS for routinely generating large-scale SIM entirely from SAR imagery across the pan-Arctic domain. [ABSTRACT FROM AUTHOR]

Most earth system models (ESMs) neglect climate feedbacks arising from carbon release from thawing permafrost, especially from thawing of subsea permafrost (SSPF). To assess the fate of SSPF in the next 1000 years, we implemented SSPF into JSBACH, the land component of the Max Planck Institute Earth System Model (MPI-ESM). This is the first implementation of SSPF processes in an ESM component. We investigate three extended scenarios from the 6th phase of the Coupled Model Intercomparison Project (CMIP6). In the 21st century only small differences are found among the scenarios, but in the upper-end emission scenario SSP5–8.5 (shared socio-economic pathway), especially in the 22nd century, SSPF ice melting is more than 15 times faster than in the pre-industrial period. In this scenario about 35 % of total SSPF volume and 34 % of SSPF area are lost by the year 3000 due to climatic changes. In the more moderate scenarios, the melting rate maximally exceeds that of pre-industrial times by a factor of 4, and the climate change induced SSPF loss (volume and area) by the year 3000 does not exceed 14 %. Our results suggest that the rate of melting of SSPF ice is related to the length of the local open-water season and thus that the easily observable sea ice concentration may be used as a proxy for the change in SSPF. [ABSTRACT FROM AUTHOR]

Glaciers in the Russian High Arctic have been subject to extensive atmospheric warming due to global climate change, yet their contribution to sea level rise has been relatively small over the past decades. Here we show surface elevation change measurements and geodetic mass balances of 93 % of all glacierized areas of Novaya Zemlya, Severnaya Zemlya, and Franz Josef Land using interferometric synthetic aperture radar measurements taken between 2010 and 2017. We calculate an overall mass loss rate of -22±6 Gt a -1 , corresponding to a sea level rise contribution of 0.06±0.02 mm a -1. Compared to measurements prior to 2010, mass loss of glaciers on the Russian archipelagos has doubled in recent years. [ABSTRACT FROM AUTHOR]

Although snow depth on sea ice is a key parameter for sea ice thickness (SIT) retrieval, there currently does not exist reliable estimations. In the Arctic, nearly all SIT products use a snow depth climatology (the modified Warren-99 climatology, W99m) constructed from in situ data obtained prior to the first significant impacts of climate change. In the Antarctic, the lack of information on snow depth remains a major obstacle in the development of reliable SIT products. In this study, we present the latest version of the altimetric snow depth (ASD) product computed over both hemispheres from the difference of the radar penetration into the snow pack between the Ka-band frequency SARAL/Altika and the Ku-band frequency CryoSat-2. The ASD solution is compared against a wide range of snow depth products including model data (Pan-Arctic Ice-Ocean Modelling and Assimilation System (PIOMAS) or its equivalent in the Antarctic the Global Ice-Ocean Modeling and Assimilation System (GIOMAS), the MERCATOR model, and NASA's Eulerian Snow On Sea Ice Model (NESOSIM, only in the Arctic)), the Advanced Microwave Scanning Radiometer-2 (AMSR2) passive radiometer data, and the Dual-altimeter Snow Thickness (DuST) Ka–Ku product (only in the Arctic). The ASD product is further validated in the Arctic against the ice mass balance (IMB) buoys, the CryoSat Validation Experiment (CryoVEx) and Operation Ice Bridge's (OIB) airborne measurements. These comparisons demonstrate that ASD is a relevant snow depth solution, with spatiotemporal patterns consistent with those of the alternative Ka–Ku DuST product but with a mean bias of about 6.5 cm. We also demonstrate that ASD is consistent with the validation data: comparisons with OIB's airborne snow radar in the Arctic during the period of 2014–2018 show a correlation of 0.66 and a RMSE of about 6 cm. Furthermore, a first-guess monthly climatology has been constructed in the Arctic from the ASD product, which shows a good agreement with OIB during 2009–2012. This climatology is shown to provide a better solution than the W99m climatology when compared with validation data. Finally, we have characterised the SIT uncertainty due to the snow depth from an ensemble of SIT solutions computed for the Arctic by using the different snow depth products previously used in the comparison with the ASD product. During the period of 2013–2019, we found a spatially averaged SIT mean standard deviation of 20 cm. Deviations between SIT estimations due to snow depths can reach up to 77 cm. Using the ASD data instead of W99m to estimate SIT over this time period leads to a reduction in the average SIT of about 30 cm. [ABSTRACT FROM AUTHOR]

Permafrost is being degraded worldwide due to the change in external forcing caused by climate change. This has also been shown to affect the morphodynamics of active rock glaciers. We studied these changes, depending on the analysis, on nine or eight active rock glaciers, respectively, with different characteristics in multiple epochs between 1953 and 2017 in Kaunertal, Austria. A combination of historical aerial photographs and airborne laser scanning data and their derivatives were used to analyse surface movement and surface elevation change. In general, the studied landforms showed a significant acceleration of varying magnitude in the epoch 1997–2006 and a volume loss to variable degrees throughout the investigation period. Rock glaciers related to glacier forefields showed significantly higher rates of subsidence than talus-connected ones. Besides, we detected two rock glaciers with deviating behaviour and one that showed an inactivation of its terminal part. By analysing meteorological data (temperature, precipitation and snow cover onset and duration), we were able to identify possible links to these external forcing parameters. The catchment-wide survey further revealed that, despite the general trend, timing, magnitude and temporal peaks of morphodynamic changes indicate a slightly different sensitivity, response or response time of individual rock glaciers to fluctuations and changes in external forcing parameters. [ABSTRACT FROM AUTHOR]

Landfast sea ice (fast ice) is an important though poorly understood component of the cryosphere on the Antarctic continental shelf, where it plays a key role in atmosphere–ocean–ice-sheet interaction and coupled ecological and biogeochemical processes. Here, we present a first in-depth baseline analysis of variability and change in circum-Antarctic fast-ice distribution (including its relationship to bathymetry), based on a new high-resolution satellite-derived time series for the period 2000 to 2018. This reveals (a) an overall trend of -882±824 km 2 yr -1 (-0.19±0.18 % yr -1) and (b) eight distinct regions in terms of fast-ice coverage and modes of formation. Of these, four exhibit positive trends over the 18-year period and four negative. Positive trends are seen in East Antarctica and in the Bellingshausen Sea, with this region claiming the largest positive trend of +1198±359 km 2 yr -1 (+1.10±0.35 % yr -1). The four negative trends predominantly occur in West Antarctica, with the largest negative trend of -1206±277 km 2 yr -1 (-1.78±0.41 % yr -1) occurring in the Victoria and Oates Land region in the western Ross Sea. All trends are significant. This new baseline analysis represents a significant advance in our knowledge of the current state of both the global cryosphere and the complex Antarctic coastal system, which are vulnerable to climate variability and change. It will also inform a wide range of other studies. [ABSTRACT FROM AUTHOR]

Short-term glacier variations can be important for water supplies or hydropower production, and glaciers are important indicators of climate change. This is why the interest in near-real-time mass balance nowcasting is considerable. Here, we address this interest and provide an evaluation of continuous observations of point mass balance based on online cameras transmitting images every 20 min. The cameras were installed on three Swiss glaciers during summer 2019, provided 352 near-real-time point mass balances in total, and revealed melt rates of up to 0.12 m water equivalent per day (mw.e.d-1) and of more than 5 mw.e. in 81 d. By means of a particle filter, these observations are assimilated into an ensemble of three TI (temperature index) and one simplified energy-balance mass balance models. State augmentation with model parameters is used to assign temporally varying weights to individual models. We analyze model performance over the observation period and find that the probability for a given model to be preferred by our procedure is 39 % for an enhanced TI model, 24 % for a simple TI model, 23 %, for a simplified energy balance model, and 14 % for a model employing both air temperature and potential solar irradiation. When compared to reference forecasts produced with both mean model parameters and parameters tuned on single mass balance observations, the particle filter performs about equally well on the daily scale but outperforms predictions of cumulative mass balance by 95 %–96 %. A leave-one-out cross-validation on the individual glaciers shows that the particle filter is also able to reproduce point observations at locations not used for model calibration. Indeed, the predicted mass balances is always within 9 % of the observations. A comparison with glacier-wide annual mass balances involving additional measurements distributed over the entire glacier mostly shows very good agreement, with deviations of 0.02, 0.07, and 0.24 mw.e. [ABSTRACT FROM AUTHOR]

Icequakes are the result of processes occurring within the ice mass or between the ice and its environment. Studying icequakes provides a unique view on ice dynamics, specifically on the basal conditions. Changes in conditions due to environmental or climate changes are reflected in icequakes. Counting and characterizing icequakes is thus essential to monitor them. Most of the icequakes recorded by the seismic station at the Belgian Princess Elisabeth Antarctica Station (PE) have small amplitudes corresponding to maximal displacements of a few nanometres. Their detection threshold is highly variable because of the rapid and strong changes in the local seismic noise level. Therefore, we evaluated the influence of katabatic winds on the noise measured by the well-protected PE surface seismometer. Our purpose is to identify whether the lack of icequake detection during some periods could be associated with variations in the processes generating them or simply with a stronger seismic noise linked to stronger wind conditions. We observed that the wind mainly influences seismic noise at frequencies greater than 1 Hz. The seismic noise power exhibits a bilinear correlation with the wind velocity, with two different slopes at a wind velocity lower and greater than 6 m s -1 and with, for example at a period of 0.26 s, a respective variation of 0.4 dB (m -1 s) and 1.4 dB (m -1 s). These results allowed a synthetic frequency and wind-speed-dependent noise model to be presented that explains the behaviour of the wind-induced seismic noise at PE, which shows that seismic noise amplitude increases exponentially with increasing wind speed. This model enables us to study the influence of the wind on the original seismic dataset, which improves the observation of cryoseismic activity near the PE station. [ABSTRACT FROM AUTHOR]

Arctic snow and ice cover are vital indicators of climate variability and change, yet while the Arctic shows overall warming and dramatic changes in snow and ice cover, the response of these high-latitude regions to recent climatic change varies regionally. Although previous studies have examined changing snow and ice separately, examining phenology changes across multiple components of the cryosphere together is important for understanding how these components and their response to climate forcing are interconnected. In this work, we examine recent changes in sea ice, lake ice, and snow together at the pan-Arctic scale using the Interactive Multisensor Snow and Ice Mapping System 24 km product from 1997–2019, with a more detailed regional examination from 2004–2019 using the 4 km product. We show overall that for sea ice, trends toward earlier open water (-7.7 d per decade, p<0.05) and later final freeze (10.6 d per decade, p<0.05) are evident. Trends toward earlier first snow-off (-4.9 d per decade, p<0.05), combined with trends toward earlier first snow-on (-2.8 d per decade, p<0.05), lead to almost no change in the length of the snow-free season, despite shifting earlier in the year. Sea ice-off, lake ice-off, and snow-off parameters were significantly correlated, with stronger correlations during the snow-off and ice-off season compared to the snow-on and ice-on season. Regionally, the Bering and Chukchi seas show the most pronounced response to warming, with the strongest trends identified toward earlier ice-off and later ice-on. This is consistent with earlier snow-off and lake ice-off and later snow-on and lake ice-on in west and southwest Alaska. In contrast to this, significant clustering between sea ice, lake ice, and snow-on trends in the eastern portion of the North American Arctic shows an earlier return of snow and ice. The marked regional variability in snow and ice phenology across the pan-Arctic highlights the complex relationships between snow and ice, as well as their response to climatic change, and warrants detailed monitoring to understand how different regions of the Arctic are responding to ongoing changes. [ABSTRACT FROM AUTHOR]

Knowledge of past ice sheet configurations is useful for informing projections of future ice sheet dynamics and for calibrating ice sheet models. The topology of grounding line retreat in the Ross Sea sector of Antarctica has been much debated, but it has generally been assumed that the modern ice sheet is as small as it has been for more than 100 000 years (Conway et al., 1999; Lee et al., 2017; Lowry et al., 2019; McKay et al., 2016; Scherer et al., 1998). Recent findings suggest that the West Antarctic Ice Sheet (WAIS) grounding line retreated beyond its current location earlier in the Holocene and subsequently readvanced to reach its modern position (Bradley et al., 2015; Kingslake et al., 2018). Here, we further constrain the post-LGM (Last Glacial Maximum) grounding line retreat and readvance in the Ross Sea sector using a two-phase model of radiocarbon input and decay in subglacial sediments from six sub-ice sampling locations. In addition, we reinterpret high basal temperature gradients, measured previously at three sites in this region (Engelhardt, 2004), which we explain as resulting from recent ice shelf re-grounding accompanying grounding line readvance. At one location – Whillans Subglacial Lake (SLW) – for which a sediment porewater chemistry profile is known, we estimate the grounding line readvance by simulating ionic diffusion. Collectively, our analyses indicate that the grounding line retreated over SLW 4300-2500+1500 years ago, and over sites on Whillans Ice Stream (WIS), Kamb Ice Stream (KIS), and Bindschadler Ice Stream (BIS) 4700-2300+1500 , 1800-700+2700 , and 1700-600+2800 years ago, respectively. The grounding line only recently readvanced back over those sites 1100-100+200 , 1500-200+500 , 1000-300+200 , and 800±100 years ago for SLW, WIS, KIS, and BIS, respectively. The timing of grounding line retreat coincided with a warm period in the mid-Holocene to late Holocene. Conversely, grounding line readvance is coincident with cooling climate in the last 1000–2000 years. Our estimates for the timing of grounding line retreat and readvance are also consistent with relatively low carbon-to-nitrogen ratios measured in our subglacial sediment samples (suggesting a marine source of organic matter) and with the lack of grounding zone wedges in front of modern grounding lines. Based on these results, we propose that the Siple Coast grounding line motions in the mid-Holocene to late Holocene were primarily driven by relatively modest changes in regional climate, rather than by ice sheet dynamics and glacioisostatic rebound, as hypothesized previously (Kingslake et al., 2018). [ABSTRACT FROM AUTHOR]

Atmospheric circulation exerts an important control on a region's snow avalanche activity by broadly determining the mountain weather patterns that influence snowpack development and avalanche release. In central Spitsbergen, the largest island in the High Arctic Svalbard archipelago, avalanches are a common natural hazard throughout the winter months. Previous work has identified a unique snow climate reflecting the region's climatically dynamic environmental setting but has not specifically addressed the synoptic-scale control of atmospheric circulation on avalanche activity here. In this work, we investigate atmospheric circulation's control on snow avalanching in the Nordenskiöld Land region of central Spitsbergen by first constructing a four-season (2016/2017–2019/2020) regional avalanche activity record using observations available on a database used by the Norwegian Water Resources and Energy Directorate (NVE). We then analyze the synoptic atmospheric conditions on days with differing avalanche activity situations. Our results show atmospheric circulation conducive to elevated precipitation, wind speeds, and air temperatures near Svalbard are associated with increased avalanche activity in Nordenskiöld Land, but different synoptic signals exist for days characterized by dry, mixed, and wet avalanche activity. Differing upwind conditions help further explain differences in the frequency and nature of avalanche activity resulting from these various atmospheric circulation patterns. We further employ a daily atmospheric circulation calendar to help contextualize our results in the growing body of literature related to climate change in this location. This work helps expand our understanding of snow avalanches in Svalbard to a broader spatial scale and provides a basis for future work investigating the impacts of climate change on avalanche activity in Svalbard and other locations where avalanche regimes are impacted by changing climatic and synoptic conditions. [ABSTRACT FROM AUTHOR]

Subarctic peatlands underlain by permafrost contain significant amounts of organic carbon. Our ability to quantify the evolution of such permafrost landscapes in numerical models is critical for providing robust predictions of the environmental and climatic changes to come. Yet, the accuracy of large-scale predictions has so far been hampered by small-scale physical processes that create a high spatial variability of thermal surface conditions, affecting the ground thermal regime and thus permafrost degradation patterns. In this regard, a better understanding of the small-scale interplay between microtopography and lateral fluxes of heat, water and snow can be achieved by field monitoring and process-based numerical modeling. Here, we quantify the topographic changes of the Šuoššjávri peat plateau (northern Norway) over a three-year period using drone-based repeat high-resolution photogrammetry. Our results show thermokarst degradation is concentrated on the edges of the plateau, representing 77 % of observed subsidence, while most of the inner plateau surface exhibits no detectable subsidence. Based on detailed investigation of eight zones of the plateau edge, we show that this edge degradation corresponds to an annual volume change of 0.13±0.07 m 3 yr -1 per meter of retreating edge (orthogonal to the retreat direction). Using the CryoGrid3 land surface model, we show that these degradation patterns can be reproduced in a modeling framework that implements lateral redistribution of snow, subsurface water and heat, as well as ground subsidence due to melting of excess ice. By performing a sensitivity test for snow depths on the plateau under steady-state climate forcing, we obtain a threshold behavior for the start of edge degradation. Small snow depth variations (from 0 to 30 cm) result in highly different degradation behavior, from stability to fast degradation. For plateau snow depths in the range of field measurements, the simulated annual volume changes are broadly in agreement with the results of the drone survey. As snow depths are clearly correlated with ground surface temperatures, our results indicate that the approach can potentially be used to simulate climate-driven dynamics of edge degradation observed at our study site and other peat plateaus worldwide. Thus, the model approach represents a first step towards simulating climate-driven landscape development through thermokarst in permafrost peatlands. [ABSTRACT FROM AUTHOR]

Debris-covered glaciers are commonly found in alpine landscapes of high relief and play an increasingly important role in a warming climate. As a result of the insulating effect of supraglacial debris, their response to changes in climate is less direct and their dynamic behaviour more complex than for debris-free glaciers. Due to a lack of observations, here we use numerical modelling to explore the dynamic interactions between debris cover and geometry evolution for an idealized glacier over centennial timescales. The main goal of this study is to understand the effects of debris cover on the glacier's transient response. To do so, we use a numerical model that couples ice flow, debris transport, and its insulating effect on surface mass balance and thereby captures dynamic feedbacks that affect the volume and length evolution. In a second step we incorporate the effects of cryokarst features such as ice cliffs and supraglacial ponds on the dynamical behaviour. Our modelling indicates that thick debris cover delays both the volume response and especially the length response to a warming climate signal. Including debris dynamics therefore results in glaciers with extended debris-covered tongues and that tend to advance or stagnate in length in response to a fluctuating climate at century timescales and hence remember the cold periods more than the warm. However, when including even a relatively small amount of melt enhancing cryokarst features in the model, the length is more responsive to periods of warming and results in substantial mass loss and thinning on debris-covered tongues, as is also observed. [ABSTRACT FROM AUTHOR]

Mean sea ice thickness is a sensitive indicator of Arctic climate change and is in long-term decline despite significant interannual variability. Current thickness estimations from satellite radar altimeters employ a snow climatology for converting range measurements to sea ice thickness, but this introduces unrealistically low interannual variability and trends. When the sea ice thickness in the period 2002–2018 is calculated using new snow data with more realistic variability and trends, we find mean sea ice thickness in four of the seven marginal seas to be declining between 60 %–100 % faster than when calculated with the conventional climatology. When analysed as an aggregate area, the mean sea ice thickness in the marginal seas is in statistically significant decline for 6 of 7 winter months. This is observed despite a 76 % increase in interannual variability between the methods in the same time period. On a seasonal timescale we find that snow data exert an increasingly strong control on thickness variability over the growth season, contributing 46 % in October but 70 % by April. Higher variability and faster decline in the sea ice thickness of the marginal seas has wide implications for our understanding of the polar climate system and our predictions for its change. [ABSTRACT FROM AUTHOR]

Long-term observations of ice phenology in lakes are ideal for studying climatic variation in time and space. We used a large set of observations from 1890 to 2020 of the timing of freeze-up and break-up, and the length of ice-free season, for 101 Norwegian lakes to elucidate variation in ice phenology across time and space. The dataset of Norwegian lakes is unusual, covering considerable variation in elevation (4–1401 m a.s.l.) and climate (from oceanic to continental) within a substantial latitudinal and longitudinal gradient (58.2–69.9 ∘ N, 4.9–30.2 ∘ E). The average date of ice break-up occurred later in spring with increasing elevation, latitude and longitude. The average date of freeze-up and the length of the ice-free period decreased significantly with elevation and longitude. No correlation with distance from the ocean was detected, although the geographical gradients were related to regional climate due to adiabatic processes (elevation), radiation (latitude) and the degree of continentality (longitude). There was a significant lake surface area effect as small lakes froze up earlier due to less volume. There was also a significant trend that lakes were completely frozen over later in the autumn in recent years. After accounting for the effect of long-term trends in the large-scale North Atlantic Oscillation (NAO) index, a significant but weak trend over time for earlier ice break-up was detected. An analysis of different time periods revealed significant and accelerating trends for earlier break-up, later freeze-up and completely frozen lakes after 1991. Moreover, the trend for a longer ice-free period also accelerated during this period, although not significantly. An understanding of the relationship between ice phenology and geographical parameters is a prerequisite for predicting the potential future consequences of climate change on ice phenology. Changes in ice phenology will have consequences for the behaviour and life cycle dynamics of the aquatic biota. [ABSTRACT FROM AUTHOR]

At high latitudes and altitudes one of the main controls on hydrological and biogeochemical processes is the breakup and freeze-up of lake and river ice. This study uses 3510 time series from across 678 Northern Hemisphere lakes and rivers to explore historical patterns in lake and river ice phenology across five overlapping time periods (1931–1960, 1946–1975, 1961–1990, 1976–2005, and 1931–2005). These time series show that the number of annual open-water days increased by 0.63 d per decade from 1931–2005 across the Northern Hemisphere, with trends for breakup and, to a lesser extent, freeze-up closely correlating with regionally averaged temperature. Breakup and freeze-up trends display a spatiotemporally complex evolution and reveal considerable caveats with interpreting the implications of ice phenology changes at lake and river sites that may only have breakup or freeze-up data, rather than both. These results provide an important contribution by showing regional variation in ice phenology trends through time that can be hidden by longer-term trends. The overlapping 30-year time periods also show evidence for an acceleration in warming trends through time. Understanding the changes on both long- and short-term timescales will be important for determining the causes of this change, the underlying biogeochemical processes associated with it, and the wider climatological significance as global temperatures rise. [ABSTRACT FROM AUTHOR]

Owing to differing and complex snow geophysical properties, radar waves of different wavelengths undergo variable penetration through snow-covered sea ice. However, the mechanisms influencing radar altimeter backscatter from snow-covered sea ice, especially at Ka- and Ku-band frequencies, and the impact on the Ka- and Ku-band radar scattering horizon or the "track point" (i.e. the scattering layer depth detected by the radar re-tracker) are not well understood. In this study, we evaluate the Ka- and Ku-band radar scattering horizon with respect to radar penetration and ice floe buoyancy using a first-order scattering model and the Archimedes principle. The scattering model is forced with snow depth data from the European Space Agency (ESA) climate change initiative (CCI) round-robin data package, in which NASA's Operation IceBridge (OIB) data and climatology are included, and detailed snow geophysical property profiles from the Canadian Arctic. Our simulations demonstrate that the Ka- and Ku-band track point difference is a function of snow depth; however, the simulated track point difference is much smaller than what is reported in the literature from the Ku-band CryoSat-2 and Ka-band SARAL/AltiKa satellite radar altimeter observations. We argue that this discrepancy in the Ka- and Ku-band track point differences is sensitive to ice type and snow depth and its associated geophysical properties. Snow salinity is first increasing the Ka- and Ku-band track point difference when the snow is thin and then decreasing the difference when the snow is thick (>0.1 m). A relationship between the Ku-band radar scattering horizon and snow depth is found. This relationship has implications for (1) the use of snow climatology in the conversion of radar freeboard into sea ice thickness and (2) the impact of variability in measured snow depth on the derived ice thickness. For both (1) and (2), the impact of using a snow climatology versus the actual snow depth is relatively small on the radar freeboard, only raising the radar freeboard by 0.03 times the climatological snow depth plus 0.03 times the real snow depth. The radar freeboard is a function of both radar scattering and floe buoyancy. This study serves to enhance our understanding of microwave interactions towards improved accuracy of snow depth and sea ice thickness retrievals via the combination of the currently operational and ESA's forthcoming Ka- and Ku-band dual-frequency CRISTAL radar altimeter missions. [ABSTRACT FROM AUTHOR]

Uncertain or inaccurate parameters in sea ice models influence seasonal predictions and climate change projections in terms of both mean and trend. We explore the feasibility and benefits of applying an ensemble Kalman filter (EnKF) to estimate parameters in the Los Alamos sea ice model (CICE). Parameter estimation (PE) is applied to the highly influential dry snow grain radius and combined with state estimation in a series of perfect model observing system simulation experiments (OSSEs). Allowing the parameter to vary in space improves performance along the sea ice edge but degrades in the central Arctic compared to requiring the parameter to be uniform everywhere, suggesting that spatially varying parameters will likely improve PE performance at local scales and should be considered with caution. We compare experiments with both PE and state estimation to experiments with only the latter and have found that the benefits of PE mostly occur after the data assimilation period, when no observations are available to assimilate (i.e., the forecast period), which suggests PE's relevance for improving seasonal predictions of Arctic sea ice. [ABSTRACT FROM AUTHOR]

The present study investigates the response of the high-latitude carbon cycle to changes in atmospheric greenhouse gas (GHG) concentrations in idealized climate change scenarios. To this end we use an adapted version of JSBACH – the land surface component of the Max Planck Institute for Meteorology Earth System Model (MPI-ESM) – that accounts for the organic matter stored in the permafrost-affected soils of the high northern latitudes. The model is run under different climate scenarios that assume an increase in GHG concentrations, based on the Shared Socioeconomic Pathway 5 and the Representative Concentration Pathway 8.5, which peaks in the years 2025, 2050, 2075 or 2100, respectively. The peaks are followed by a decrease in atmospheric GHGs that returns the concentrations to the levels at the beginning of the 21st century, reversing the imposed climate change. We show that the soil CO2 emissions exhibit an almost linear dependence on the global mean surface temperatures that are simulated for the different climate scenarios. Here, each degree of warming increases the fluxes by, very roughly, 50 % of their initial value, while each degree of cooling decreases them correspondingly. However, the linear dependence does not mean that the processes governing the soil CO2 emissions are fully reversible on short timescales but rather that two strongly hysteretic factors offset each other – namely the net primary productivity and the availability of formerly frozen soil organic matter. In contrast, the soil methane emissions show a less pronounced increase with rising temperatures, and they are consistently lower after the peak in the GHG concentrations than prior to it. Here, the net fluxes could even become negative, and we find that methane emissions will play only a minor role in the northern high-latitude contribution to global warming, even when considering the high global warming potential of the gas. Finally, we find that at a global mean temperature of roughly 1.75 K (±0.5 K) above pre-industrial levels the high-latitude ecosystem turns from a CO2 sink into a source of atmospheric carbon, with the net fluxes into the atmosphere increasing substantially with rising atmospheric GHG concentrations. This is very different from scenario simulations with the standard version of the MPI-ESM, in which the region continues to take up atmospheric CO2 throughout the entire 21st century, confirming that the omission of permafrost-related processes and the organic matter stored in the frozen soils leads to a fundamental misrepresentation of the carbon dynamics in the Arctic. [ABSTRACT FROM AUTHOR]

The retreat of ice shelves and glaciers over the last century provides unequivocal evidence of recent global warming. Glacierets (miniature glaciers) and ice patches are important components of the cryosphere that highlight the global retreat of glaciers, but knowledge of their behaviour prior to the Little Ice Age is lacking. Here, we report the uranium–thorium age of subglacial carbonate deposits from a recently exposed surface previously occupied by the disappearing Triglav Glacier (southeastern European Alps) that may elucidate the glacier's presence throughout the entire Holocene. The ages suggest the deposits' possible preservation since the Last Glacial Maximum and Younger Dryas. These thin deposits, formed by regelation, are easily eroded if exposed during previous Holocene climatic optima. The age data indicate the glacier's present unprecedented level of retreat since the Last Glacial Maximum and the potential of subglacial carbonates as additional proxies to highlight the extraordinary nature of the current global climatic changes. [ABSTRACT FROM AUTHOR]

We have studied the evolution of the Greenland ice sheet under a range of constant climates typical of those projected for the end of the present century using a dynamical ice sheet model (Glimmer) coupled to an atmosphere general circulation model (FAMOUS–ice AGCM). The ice sheet surface mass balance (SMB) is simulated within the AGCM by a multilayer snow scheme from snowfall and surface energy fluxes, including refreezing and dependence on altitude within AGCM grid boxes. Over millennia under any warmer climate, the ice sheet reaches a new steady state, whose mass is correlated with the magnitude of global climate change imposed. If a climate that gives the recently observed SMB were maintained, global-mean sea level rise (GMSLR) would reach 0.5–2.5 m. For any global warming exceeding 3 K, the contribution to GMSLR exceeds 5 m. For the largest global warming considered (about +5 K), the rate of GMSLR is initially 2.7 mm yr -1 , and eventually only a small ice cap endures, resulting in over 7 m of GMSLR. Our analysis gives a qualitatively different impression from previous work in that we do not find a sharp threshold warming that divides scenarios in which the ice sheet suffers little reduction from those in which it is mostly lost. The final steady state is achieved by withdrawal from the coast in some places and a tendency for increasing SMB due to enhancement of cloudiness and snowfall over the remaining ice sheet by the effects of topographic change on atmospheric circulation, outweighing the tendency for decreasing SMB from the reduction in surface altitude. If late 20th-century climate is restored after the ice sheet mass has fallen below a threshold of about 4 m of sea level equivalent, it will not regrow to its present extent because the snowfall in the northern part of the island is reduced once the ice sheet retreats from there. In that case, about 2 m of GMSLR would become irreversible. In order to avoid this outcome, anthropogenic climate change must be reversed before the ice sheet has declined to the threshold mass, which would be reached in about 600 years at the highest rate of mass loss within the likely range of the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. [ABSTRACT FROM AUTHOR]