Search

In some areas of Antarctica, blue-colored bare ice is exposed at the surface. These blue ice areas (BIAs) can trap meteorites or old ice and are vital for understanding the climatic history. By combining multi-sensor remote sensing data (MODIS, RADARSAT-2, and TanDEM-X) in a deep learning framework, we map blue ice across the continent at 200-m resolution. We use a novel methodology for image segmentation with “noisy” labels to learn an underlying “clean” pattern with a neural network. In total, BIAs cover ca. 140,000 km2 (∼1%) of Antarctica, of which nearly 50% located within 20 km of the grounding line. There, the low albedo of blue ice enhances melt-water production and its mapping is crucial for mass balance studies that determine the stability of the ice sheet. Moreover, the map provides input for fieldwork missions and can act as constraint for other geophysical mapping efforts., Mathematical Geodesy and Positioning, Civil Engineering & Geosciences

Correction to: Nature Reviews Earth & Environment https://doi.org/10.1038/s43017-023-00507-9, published online 23 January 2024. In the version of the article initially published, in Fig. 5, under “Radar altimeter”, “O(16–160 m)” previously read “O(16–160 km)”. This has now been corrected in the HTML and PDF versions of the article., Mathematical Geodesy and Positioning, Physical and Space Geodesy

Imagery acquired by the Moderate-resolution Imaging Spectroradiometer (MODIS) provides a global archive of dailyNormalized Difference Snow Index (NDSI) at 500 m nominal resolution since the year 2000. While Sentinel-2 (S2) NDSI provides an increased spatial resolution of 20 m since the year 2015, the temporal resolution amounts to only 5 days and thus lacks the high temporal resolution of MODIS. Efforts to combine NDSI datasets for an increased temporal and spatial resolution have so far focused on the deriving binary snow cover maps or combining data from other sensors. In contrast, we produce fine scale (20 m) fractional snow cover (FSC) by downscaling MODIS NDSI to S2 resolution. Random forest regression predicts S2 NDSI based on dynamic features (MODIS NDSI, day-of-year) and static, topographic features for an alpine study site. Subsequently, FSC is derived from S2 NDSI. Cross-validation results in R2 of 0.795 and RMSE of 0.155 for FSC and outperforms common resampling methods. Multi-annual S2 NDSI metrics are able to slightly improve model accuracy. Our results suggest that combining topographical data and low-resolution NDSI allows to produce daily, high-resolution S2 NDSI and FSC and improve fine scale characterization of snow cover dynamics in mountain landscapes., Mathematical Geodesy and Positioning

More than 60% of meteorite finds on Earth originate from Antarctica. Using a data-driven analysis that identifies meteorite-rich sites in Antarctica, we show climate warming causes many extraterrestrial rocks to be lost from the surface by melting into the ice sheet. At present, approximately 5,000 meteorites become inaccessible per year (versus ~1,000 finds per year) and, independent of the emissions scenario, ~24% will be lost by 2050, potentially rising to ∼76% by 2100 under a high-emissions scenario., Mathematical Geodesy and Positioning

While the influence of surface melt on Antarctic ice shelf stability can be large, the duration and affected area of melt events are often small. Therefore, melt events are difficult to capture with remote sensing, as satellite sensors always face the trade-off between spatial and temporal resolution. To overcome this limitation, we developed UMelt: a surface melt record for all Antarctic ice shelves with a high spatial (500 m) and high temporal (12 h) resolution for the period 2016–2021. Our approach is based on a deep learning model, specifically a U-Net, which was developed in Google Earth Engine. The U-Net combines microwave remote sensing observations from three sources: Sentinel-1, Special Sensor Microwave Imager/Sounder (SSMIS), and Advanced Scatterometer (ASCAT). The U-Net was trained on the Shackleton Ice Shelf for melt seasons 2017–2021, using the fine-scale melt patterns of Sentinel-1 as reference data and SSMIS, ASCAT, a digital elevation model, and multi-year Sentinel-1 melt fraction as predictors. The trained U-Net performed well on the Shackelton Ice Shelf for test melt season 2016–2017 (accuracy: 91.3%; F1-score: 86.9%), and the Larsen C Ice Shelf, which was not considered during training (accuracy: 91.0%; F1-score: 89.3%). Using the trained U-Net model, we have successfully developed the UMelt record. UMelt allows Antarctic-wide surface melt to be detected at a small scale while preserving a high temporal resolution, which could lead to new insights into the response of ice shelves to a changing atmospheric forcing., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Muddy coasts provide ecological habitats, supply food and form a natural coastal defence. Relative sea level rise, changing wave energy and human interventions will increase the pressure on muddy coastal zones. For sustainable coastal management it is key to obtain information on the geomorphology of and historical changes along muddy areas. So far, little is known about the distribution and behaviour of muddy coasts at a global scale. In this study we present a global scale assessment of the occurrence of muddy coasts and rates of coastline change therein. We combine publicly available satellite imagery and coastal geospatial datasets, to train an automated classification method to identify muddy coasts. We find that 14% of the world’s ice-free coastline is muddy, of which 60% is located in the tropics. Furthermore, the majority of the world’s muddy coasts are eroding at rates exceeding 1 m/yr over the last three decades., Coastal Engineering, Environmental Fluid Mechanics, Mathematical Geodesy and Positioning, Civil Engineering & Geosciences

Antarctic ice sheet (AIS) mass loss is predominantly driven by increased solid ice discharge, but its variability is governed by surface processes. Snowfall fluctuations control the surface mass balance (SMB) of the grounded AIS, while meltwater ponding can trigger ice shelf collapse potentially accelerating discharge. Surface processes are essential to quantify AIS mass change, but remain poorly represented in climate models typically running at 25-100 km resolution. Here we present SMB and surface melt products statistically downscaled to 2 km resolution for the contemporary climate (1979-2021) and low, moderate and high-end warming scenarios until 2100. We show that statistical downscaling modestly enhances contemporary SMB (3%), which is sufficient to reconcile modelled and satellite mass change. Furthermore, melt strongly increases (46%), notably near the grounding line, in better agreement with in-situ and satellite records. The melt increase persists by 2100 in all warming scenarios, revealing higher surface melt rates than previously estimated., Physical and Space Geodesy, Mathematical Geodesy and Positioning

The intrusion of Circumpolar Deep Water in the Amundsen and Bellingshausen Sea embayments of Antarctica causes ice shelves in the region to melt from below, potentially putting their stability at risk. Earlier studies have shown how digital elevation models can be used to obtain ice shelf basal melt rates at a high spatial resolution. However, there has been limited availability of high-resolution elevation data, a gap the Reference Elevation Model of Antarctica (REMA) has filled. In this study we use a novel combination of REMA and CryoSat-2 elevation data to obtain high-resolution basal melt rates of the Dotson Ice Shelf in a Lagrangian framework, at a 50 m spatial posting on a 3-yearly temporal resolution. We present a novel method: Basal melt rates Using REMA and Google Earth Engine (BURGEE). The high resolution of BURGEE is supported through a sensitivity study of the Lagrangian displacement. The high-resolution basal melt rates show a good agreement with an earlier basal melt product based on CryoSat-2. Both products show a wide melt channel extending from the grounding line to the ice front, but our high-resolution product indicates that the pathway and spatial variability of this channel is influenced by a pinning point on the ice shelf. This result emphasizes the importance of high-resolution basal melt rates to expand our understanding of channel formation and melt patterns. BURGEE can be expanded to a pan-Antarctic study of high-resolution basal melt rates. This will provide a better picture of the (in)stability of Antarctic ice shelves., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Pine Island Glacier (PIG) has recently experienced increased ice loss that has mostly been attributed to basal melt and ocean ice dynamics. However, atmospheric forcing also plays a role in the ice mass budget, as besides lower-latitude warm air intrusions, the steeply sloping terrain that surrounds the glacier promotes frequent Foehn winds. An investigation of 41 years of reanalysis data reveals that Foehn occurs more frequently from June to October, with Foehn episodes typically lasting about 5 to 9 h. An analysis of the surface mass balance indicated that their largest impact is on the surface sublimation, which is increased by about 1.43 mm water equivalent (w.e.) per day with respect to no-Foehn events. Blowing snow makes roughly the same contribution as snowfall, around 0.34–0.36 mm w.e. d−1, but with the opposite sign. The melting rate is 3 orders of magnitude smaller than the surface sublimation rate. The negative phase of the Antarctic oscillation and the positive phase of the Southern Annular Mode promote the occurrence of Foehn at PIG. A particularly strong event took place on 9–11 November 2011, when 10 m winds speeds in excess of 20 m s−1 led to downward sensible heat fluxes higher than 75 W m−2 as they descended the mountainous terrain. Surface sublimation and blowing-snow sublimation dominated the surface mass balance, with magnitudes of up to 0.13 mm w.e. h−1. Satellite data indicated an hourly surface melting area exceeding 100 km2. Our results stress the importance of the atmospheric forcing on the ice mass balance at PIG., Mathematical Geodesy and Positioning

It has been argued that the −5 °C annual mean 2 m air temperature isotherm defines a limit of ice shelf viability on the Antarctic Peninsula as melt ponding increases at higher temperatures. It is, however, presently unknown whether this threshold can also be applied to other Antarctic ice shelves. Here we use two present-day and three future high-resolution Antarctic climate simulations to predict warming thresholds for Antarctic ice shelf melt pond formation on the basis of the melt-over-accumulation ratio. The associated warming thresholds match well with observed melt pond volumes and are found to be spatially highly variable and controlled by snow accumulation. For relatively wet ice shelves, the −5 °C temperature threshold was confirmed; but cold and dry ice shelves such as Amery, Ross and Filchner-Ronne are more vulnerable than previously thought, with threshold temperatures well below −15 °C. Coupled Model Intercomparison Project Phase 6 models predict that towards the end of this century these thresholds can be reached on many ice shelves, even on cold ice shelves and under moderate warming scenarios., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Meteorites provide a unique view into the origin and evolution of the Solar System. Antarctica is the most productive region for recovering meteorites, where these extraterrestrial rocks concentrate at meteorite stranding zones. To date, meteorite-bearing blue ice areas are mostly identified by serendipity and through costly reconnaissance missions. Here, we identify meteorite-rich areas by combining state-of-the-art datasets in a machine learning algorithm and provide continent-wide estimates of the probability to find meteorites at any given location. The resulting set of ca. 600 meteorite stranding zones, with an estimated accuracy of over 80%, reveals the existence of unexplored zones, some of which are located close to research stations. Our analyses suggest that less than 15% of all meteorites at the surface of the Antarctic ice sheet have been recovered to date. The data-driven approach will greatly facilitate the quest to collect the remaining meteorites in a coordinated and cost-effective manner., Mathematical Geodesy and Positioning, Pattern Recognition and Bioinformatics

Icelandic glaciers have been losing mass since the Little Ice Age in the mid-to-late 1800s, with higher mass loss rates in the early 21st century, followed by a slowdown since 2011. As of yet, it remains unclear whether this mass loss slowdown will persist in the future. By reconstructing the contemporary (1958–2019) surface mass balance of Icelandic glaciers, we show that the post-2011 mass loss slowdown coincides with the development of the Blue Blob, an area of regional cooling in the North Atlantic Ocean to the south of Greenland. This regional cooling signal mitigates atmospheric warming in Iceland since 2011, in turn decreasing glacier mass loss through reduced meltwater runoff. In a future high-end warming scenario, North Atlantic cooling is projected to mitigate mass loss of Icelandic glaciers until the mid-2050s. High mass loss rates resume thereafter as the regional cooling signal weakens., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Glaciers play a crucial role in the Earth System: they are important water suppliers to lower-lying areas during hot and dry periods, and they are major contributors to the observed present-day sea-level rise. Glaciers can also act as a source of natural hazards and have a major touristic value. Given their societal importance, there is large scientific interest in better understanding and accurately simulating the temporal evolution of glaciers, both in the past and in the future. Here, we give an overview of the state of the art of simulating the evolution of individual glaciers over decadal to centennial time scales with ice-dynamical models. We hereby highlight recent advances in the field and emphasize how these go hand-in-hand with an increasing availability of on-site and remotely sensed observations. We also focus on the gap between simplified studies that use parameterizations, typically used for regional and global projections, and detailed assessments for individual glaciers, and explain how recent advances now allow including ice dynamics when modeling glaciers at larger spatial scales. Finally, we provide concrete recommendations concerning the steps and factors to be considered when modeling the evolution of glaciers. We suggest paying particular attention to the model initialization, analyzing how related uncertainties in model input influence the modeled glacier evolution and strongly recommend evaluating the simulated glacier evolution against independent data., Mathematical Geodesy and Positioning

Satellite radar altimetry has been an important tool for cryospheric applications such as measuring ice-sheet height or assessing anomalies in snow and ice properties (e.g. the extensive melt in Greenland in 2012). Although accurate height measurements are key for such applications, slope-induced errors due to undulating topography within the kilometre-wide beam-limited footprint can cause multi-metre errors. Two main correction methods that have been developed (referred to as the slope- and point-based methods) neglect either the actual topography or the actual footprint that can be estimated by a combination of the leading edge and topography. Therefore, a leading edge point-based (LEPTA) method is presented that corrects for the slope-induced error by including the leading edge information of the radar waveform to determine the impact point. The principle of the method is that only the points on the ground that are within the range determined by the beginning and end of the leading edge are used to determine the impact point. Benchmarking of the LEPTA method against the slope- and point-based methods based on CryoSat-2 Low Resolution Mode (LRM) acquisitions over Greenland in 2019 shows that, when compared to ICESat-2 observations, the LEPTA method has a stable performance both in the flat, interior regions of Greenland and in regions with more complex topography. The median difference between the slope-corrected CryoSat-2 heights using LEPTA and the ICESat-2 heights is at the millimetre level, whereas the slope and point-based methods can have a 0.21 and 0.48 m difference, respectively, and the Level-2I (L2I) data provided by ESA have a 0.01 m difference. The median absolute deviation of height differences between CryoSat-2 and ICESat-2, which we use as an indicator of the variation in errors, is also the lowest for LEPTA (0.09 m) in comparison to the aforementioned methods (0.19 m for slope method and 0.10 m for point-based method) and ESA Level-2 data (0.14 m)., Aerospace Engineering, Mathematical Geodesy and Positioning

Satellite radar altimetry has been an important tool for cryospheric applications such as measuring ice-sheet height or assessing anomalies in snow and ice properties (e.g. the extensive melt in Greenland in 2012). Although accurate height measurements are key for such applications, slope-induced errors due to undulating topography within the kilometre-wide beam-limited footprint can cause multi-metre errors. Two main correction methods that have been developed (referred to as the slope- and point-based methods) neglect either the actual topography or the actual footprint that can be estimated by a combination of the leading edge and topography. Therefore, a leading edge point-based (LEPTA) method is presented that corrects for the slope-induced error by including the leading edge information of the radar waveform to determine the impact point. The principle of the method is that only the points on the ground that are within the range determined by the beginning and end of the leading edge are used to determine the impact point. Benchmarking of the LEPTA method against the slope- and point-based methods based on CryoSat-2 Low Resolution Mode (LRM) acquisitions over Greenland in 2019 shows that, when compared to ICESat-2 observations, the LEPTA method has a stable performance both in the flat, interior regions of Greenland and in regions with more complex topography. The median difference between the slope-corrected CryoSat-2 heights using LEPTA and the ICESat-2 heights is at the millimetre level, whereas the slope and point-based methods can have a 0.21 and 0.48 m difference, respectively, and the Level-2I (L2I) data provided by ESA have a 0.01 m difference. The median absolute deviation of height differences between CryoSat-2 and ICESat-2, which we use as an indicator of the variation in errors, is also the lowest for LEPTA (0.09 m) in comparison to the aforementioned methods (0.19 m for slope method and 0.10 m for point-based method) and ESA Level-2 data (0.14 m)., Aerospace Engineering, Mathematical Geodesy and Positioning

Satellite radar altimetry has been an important tool for cryospheric applications such as measuring ice-sheet height or assessing anomalies in snow and ice properties (e.g. the extensive melt in Greenland in 2012). Although accurate height measurements are key for such applications, slope-induced errors due to undulating topography within the kilometre-wide beam-limited footprint can cause multi-metre errors. Two main correction methods that have been developed (referred to as the slope- and point-based methods) neglect either the actual topography or the actual footprint that can be estimated by a combination of the leading edge and topography. Therefore, a leading edge point-based (LEPTA) method is presented that corrects for the slope-induced error by including the leading edge information of the radar waveform to determine the impact point. The principle of the method is that only the points on the ground that are within the range determined by the beginning and end of the leading edge are used to determine the impact point. Benchmarking of the LEPTA method against the slope- and point-based methods based on CryoSat-2 Low Resolution Mode (LRM) acquisitions over Greenland in 2019 shows that, when compared to ICESat-2 observations, the LEPTA method has a stable performance both in the flat, interior regions of Greenland and in regions with more complex topography. The median difference between the slope-corrected CryoSat-2 heights using LEPTA and the ICESat-2 heights is at the millimetre level, whereas the slope and point-based methods can have a 0.21 and 0.48 m difference, respectively, and the Level-2I (L2I) data provided by ESA have a 0.01 m difference. The median absolute deviation of height differences between CryoSat-2 and ICESat-2, which we use as an indicator of the variation in errors, is also the lowest for LEPTA (0.09 m) in comparison to the aforementioned methods (0.19 m for slope method and 0.10 m for point-based method) and ESA Level-2 data (0.14 m)., Mathematical Geodesy and Positioning, Physical and Space Geodesy

Surface melt is an important driver of ice shelf disintegration and its consequent mass loss over the Antarctic Ice Sheet. Monitoring surface melt using satellite remote sensing can enhance our understanding of ice shelf stability. However, the sensors do not measure the actual physical process of surface melt, but rather observe the presence of liquid water. Moreover, the sensor observations are influenced by the sensor characteristics and surface properties. Therefore, large inconsistencies can exist in the derived melt estimates from different sensors. In this study, we apply state-of-the-art melt detection algorithms to four frequently used remote sensing sensors, i.e., two active microwave sensors, which are Advanced Scatterometer (ASCAT) and Sentinel-1, a passive microwave sensor, i.e., Special Sensor Microwave Imager/Sounder (SSMIS), and an optical sensor, i.e., Moderate Resolution Imaging Spectroradiometer (MODIS). We intercompare the melt detection results over the entire Antarctic Ice Sheet and four selected study regions for the melt seasons 2015-2020. Our results show large spatiotemporal differences in detected melt between the sensors, with particular disagreement in blue ice areas, in aquifer regions, and during wintertime surface melt. We discuss that discrepancies between sensors are mainly due to cloud obstruction and polar darkness, frequency-dependent penetration of satellite signals, temporal resolution, and spatial resolution, as well as the applied melt detection methods. Nevertheless, we argue that different sensors can complement each other, enabling improved detection of surface melt over the Antarctic Ice Sheet., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Areas of structural damage mechanically weaken Antarctic ice shelves. This potentially preconditions ice shelves for disintegration and enhanced grounding line retreat. The development of damage and its feedback on marine ice sheet dynamics has been identified as key to future ice shelf stability and sea level contributions from Antarctica. However, it is one of the least understood processes that impact ice shelf instability since quantifying damage efficiently and accurately is a challenging task. Challenges relate to the complex surface of Antarctica, variations in viewing-illumination geometry, snow or cloud cover and variable signal-to-noise levels in satellite imagery. Therefore, automated damage assessment approaches require careful pre- and post-processing, lacking the option to be applied to wider spatiotemporal domains. Simultaneously, studies that use manual mapping are usually limited due to the effort required for extensive mapping, which either results in a limited spatial domain or the use of low resolution data. This study proposes the NormalisEd Radon transform Damage detection (NeRD) method to detect damage features and their orientations from multi-source satellite imagery. NeRD performs robust, high resolution, large-scale damage assessments. NeRD is applied to the ice shelves in the Amundsen Sea Embayment (ASE) and validated with both manually labelled and existing fracture maps. Validation shows that NeRD detects damage with high recall and provides an accurate physical representation of multi-scale damage features and their orientation. Sensitivity analyses indicate NeRD is robust to different resolution parameter settings. NeRD consistently detects damage for different data sources ranging from optical Landsat 7/8 and Sentinel-2 optical to Synthetic Aperture Radar Sentinel-1 data. Therefore, NeRD paves the way for synergistic multi-source damage detection that overcomes remaining limitations from individual sources. Results show vary, Mathematical Geodesy and Positioning

A devastating tsunami struck Palu Bay in the wake of the 28 September 2018 Mw = 7.5 Palu earthquake (Sulawesi, Indonesia). With a predominantly strike-slip mechanism, the question remains whether this unexpected tsunami was generated by the earthquake itself, or rather by earthquake-induced landslides. In this study we examine the tsunami potential of the co-seismic deformation. To this end, we present a novel geodetic data set of Global Positioning System and multiple Synthetic Aperture Radar-derived displacement fields to estimate a 3D co-seismic surface deformation field. The data reveal a number of fault bends, conforming to our interpretation of the tectonic setting as a transtensional basin. Using a Bayesian framework, we provide robust finite fault solutions of the co-seismic slip distribution, incorporating several scenarios of tectonically feasible fault orientations below the bay. These finite fault scenarios involve large co-seismic uplift (>2 m) below the bay due to thrusting on a restraining fault bend that connects the offshore continuation of two parallel onshore fault segments. With the co-seismic displacement estimates as input we simulate a number of tsunami cases. For most locations for which video-derived tsunami waveforms are available our models provide a qualitative fit to leading wave arrival times and polarity. The modeled tsunamis explain most of the observed runup. We conclude that co-seismic deformation was the main driver behind the tsunami that followed the Palu earthquake. Our unique geodetic data set constrains vertical motions of the sea floor, and sheds new light on the tsunamigenesis of strike-slip faults in transtensional basins., Astrodynamics & Space Missions, Environmental Fluid Mechanics, Mathematical Geodesy and Positioning, Space Engineering, Physical and Space Geodesy

The calving of Antarctic ice shelves remains unpredictable to date due to a lack of understanding of the role of the different climatic components in such events. In this study, the role of atmospheric forcing in the calving of the Brunt Ice Shelf (BIS) in February 2021 is investigated using a combination of observational and reanalysis data. The occurrence of a series of extreme cyclones around the time of the calving induced an oceanward sea-surface slope of >0.08° leading to the calving along a pre-existing rift. The severe storms were sustained by the development of a pressure dipole on both sides of the BIS associated with a La Niña event and the positive phase of the Southern Annular Mode. Poleward advection of warm and moist low-latitude air over the BIS area just before the calving was also observed in association with atmospheric rivers accompanying the cyclones. Immediately after the calving, strong offshore winds continued and promoted the drift of the iceberg A-74 in the Weddell Sea at a speed up to 700 m day−1. This study highlights the contribution of local atmospheric conditions to ice-shelf dynamics. The link to the larger scale circulation patterns indicates that both need to be accounted for in the projections of Antarctic ice shelf evolution., Mathematical Geodesy and Positioning

Antarctic blue ice areas are exposed due to erosion and sublimation of snow. At the same time, surface melt can form surface types that are spectrally similar to blue ice, especially at low elevations. These are termed melt-induced blue ice areas. Both types of blue ice are sensitive indicators of climate change. Satellite remote sensing is a powerful technique to retrieve the spatial extent of blue ice areas and their variation in time. Yet, existing satellite-derived blue ice area products are either mono-temporal for the entire Antarctic ice sheet, or multi-temporal for a limited area. Here, we present FABIAN, a product of blue ice fraction over Antarctica, derived from the moderate resolution imaging spectroradiometer (MODIS) archive covering the period 2000–2021. A spectral mixture analysis (SMA) in Google Earth Engine, based on a careful selection of endmember spectra, accurately reconstructs the reflectance observed by MODIS in blue ice areas. Based on a validation with contemporaneous Sentinel-2 images, FABIAN has a root mean square error in blue ice fraction of approximately 10% ∼ 20% in wind-induced blue ice areas, and 20% ∼ 30% in melt-induced blue ice areas across six selected test sites in the coastal East Antarctic ice sheet. FABIAN is challenged in regions with shallow melt streams and lakes, since their spectral profiles are similar to those from blue ice areas in MODIS bands. For further analyses and applications, FABIAN holds the potential for (1) deriving annual blue ice area maps, (2) distinguishing between wind-and melt-induced blue ice types, (3) evaluating and correcting (regional) climate models, and (4) analyzing temporal variations in blue ice abundance and exposure., Mathematical Geodesy and Positioning

Glaciers play a crucial role in the Earth System: they are important water suppliers to lower-lying areas during hot and dry periods, and they are major contributors to the observed present-day sea-level rise. Glaciers can also act as a source of natural hazards and have a major touristic value. Given their societal importance, there is large scientific interest in better understanding and accurately simulating the temporal evolution of glaciers, both in the past and in the future. Here, we give an overview of the state of the art of simulating the evolution of individual glaciers over decadal to centennial time scales with ice-dynamical models. We hereby highlight recent advances in the field and emphasize how these go hand-in-hand with an increasing availability of on-site and remotely sensed observations. We also focus on the gap between simplified studies that use parameterizations, typically used for regional and global projections, and detailed assessments for individual glaciers, and explain how recent advances now allow including ice dynamics when modeling glaciers at larger spatial scales. Finally, we provide concrete recommendations concerning the steps and factors to be considered when modeling the evolution of glaciers. We suggest paying particular attention to the model initialization, analyzing how related uncertainties in model input influence the modeled glacier evolution and strongly recommend evaluating the simulated glacier evolution against independent data., Mathematical Geodesy and Positioning

Icelandic glaciers have been losing mass since the Little Ice Age in the mid-to-late 1800s, with higher mass loss rates in the early 21st century, followed by a slowdown since 2011. As of yet, it remains unclear whether this mass loss slowdown will persist in the future. By reconstructing the contemporary (1958–2019) surface mass balance of Icelandic glaciers, we show that the post-2011 mass loss slowdown coincides with the development of the Blue Blob, an area of regional cooling in the North Atlantic Ocean to the south of Greenland. This regional cooling signal mitigates atmospheric warming in Iceland since 2011, in turn decreasing glacier mass loss through reduced meltwater runoff. In a future high-end warming scenario, North Atlantic cooling is projected to mitigate mass loss of Icelandic glaciers until the mid-2050s. High mass loss rates resume thereafter as the regional cooling signal weakens., Physical and Space Geodesy, Mathematical Geodesy and Positioning

Meteorites provide a unique view into the origin and evolution of the Solar System. Antarctica is the most productive region for recovering meteorites, where these extraterrestrial rocks concentrate at meteorite stranding zones. To date, meteorite-bearing blue ice areas are mostly identified by serendipity and through costly reconnaissance missions. Here, we identify meteorite-rich areas by combining state-of-the-art datasets in a machine learning algorithm and provide continent-wide estimates of the probability to find meteorites at any given location. The resulting set of ca. 600 meteorite stranding zones, with an estimated accuracy of over 80%, reveals the existence of unexplored zones, some of which are located close to research stations. Our analyses suggest that less than 15% of all meteorites at the surface of the Antarctic ice sheet have been recovered to date. The data-driven approach will greatly facilitate the quest to collect the remaining meteorites in a coordinated and cost-effective manner., Mathematical Geodesy and Positioning, Pattern Recognition and Bioinformatics

Accurately estimating the surface melt volume of the Antarctic Ice Sheet is challenging and has hitherto relied on climate modeling or observations from satellite remote sensing. Each of these methods has its limitations, especially in regions with high surface melt. This study aims to demonstrate the potential of improving surface melt simulations with a regional climate model by deploying a deep learning model. A deep-learning-based framework has been developed to correct surface melt from the regional atmospheric climate model version 2.3p2 (RACMO2), using meteorological observations from automatic weather stations (AWSs) and surface albedo from satellite imagery. The framework includes three steps: (1) training a deep multilayer perceptron (MLP) model using AWS observations, (2) correcting Moderate Resolution Imaging Spectroradiometer (MODIS) albedo observations, and (3) using these two to correct the RACMO2 surface melt simulations. Using observations from three AWSs at the Larsen B and C ice shelves, Antarctica, cross-validation shows a high accuracy (root-mean-square error of 0.95ĝ€¯mmĝ€¯w.e.ĝ€¯d-1, mean absolute error of 0.42ĝ€¯mmĝ€¯w.e.ĝ€¯d-1, and a coefficient of determination of 0.95). Moreover, the deep MLP model outperforms conventional machine learning models and a shallow MLP model. When applying the trained deep MLP model over the entire Larsen Ice Shelf, the resulting corrected RACMO2 surface melt shows a better correlation with the AWS observations for two out of three AWSs. However, for one location (AWS 18), the deep MLP model does not show improved agreement with AWS observations; this is likely because surface melt is largely driven by factors (e.g., air temperature, topography, katabatic wind) other than albedo within the corresponding coarse-resolution model pixels. Our study demonstrates the opportunity to improve surface melt simulations using deep learning combined with satellite albedo observations. However, more work is required to re, Mathematical Geodesy and Positioning

Surface meltwater drains on several Antarctic ice shelves, resulting in surface and sub-surface lakes that are potentially critical for the ice shelf collapse. Despite these phenomena, our understanding and assessment of the drainage and refreezing of these lakes is limited, mainly due to lack of field observations and to the limitations of optical satellite imagery during polar night and in cloudy conditions. This paper explores the potential of backscatter intensity and of interferometric coherence and phase from synthetic aperture radar (SAR) imagery as an alternative to assess the dynamics of meltwater lakes. In four case study regions over Amery and Roi Baudouin ice shelves, East Antarctica, we examine spatial and temporal variations in SAR backscatter intensity and interferometric (InSAR) coherence and phase over several lakes derived from Sentinel-1A/B C-band SAR imagery. Throughout the year, the lakes are observed in a completely frozen state, in a partially frozen state with a floating ice lid and as open-water lakes. Our analysis reveals that the meltwater lake delineation is challenging during the melting period when the contrast between melting snow and lakes is indistinguishable. Despite this finding, we show using a combination of backscatter and InSAR observations that lake dynamics can be effectively captured during other non-summertime months. Moreover, our findings highlight the utility of InSAR-based observations for discriminating between refrozen ice and sub-surface meltwater and indicate the potential for phase-based detection and monitoring of rapid meltwater drainage events. The potential of this technique to monitor these meltwater change events is, however, strongly determined by the satellite revisit interval and potential changes in scattering properties due to snowfall or melt events., Mathematical Geodesy and Positioning

Physically based snow models provide valuable information on snow cover evolution and are therefore key to provide water availability projections. Yet, uncertainties related to snow modelling remain large as a result of differences in the representation of snow physics and meteorological forcing. While many studies focus on evaluating these uncertainties, no snow model comparison has been done in environments where sublimation is the main ablation process. This study evaluates a case study in the semi-arid Andes of Chile and aims to compare two snow models with different complexities, SNOWPACK and SnowModel, at a local point over one snow season and to evaluate their sensitivity relative to parameterisation and forcing. For that purpose, the two models are forced with (i) the most ideal set of input parameters, (ii) an ensemble of different physical parameterisations, and (iii) an ensemble of biased forcing. Results indicate large uncertainties depending on forcing, the snow roughness length z0, albedo parameterisation, and fresh snow density parameterisation. The uncertainty caused by the forcing is directly related to the bias chosen. Even though the models show significant differences in their physical complexity, the snow model choice is of least importance, as the sensitivity of both models to the forcing data was on the same order of magnitude and highly influenced by the precipitation uncertainties. The sublimation ratio ranges are in agreement for the two models: 36.4 % to 80.7 % for SnowModel and 36.3 % to 86.0 % for SNOWPACK, and are related to the albedo parameterisation and snow roughness length choice for the two models., Mathematical Geodesy and Positioning

Surface albedo is an important forcing parameter that drives the radiative energy budget as it determines the fraction of the downwelling solar irradiance that the surface reflects. Here we report on ground-based measurements of the spectral albedo (350–2200 nm) carried out at 20 sites across a North–South transect of approximately 1300 km in the Atacama Desert, from latitude 18° S to latitude 30° S. These spectral measurements were used to evaluate remote sensing estimates of the albedo derived from the Moderate Resolution Imaging Spectroradiometer (MODIS). We found that the relative mean bias error (RMBE) of MODIS-derived estimates was within ± 5% of ground-based measurements in most of the Atacama Desert (18–27° S). Although the correlation between MODIS-derived estimates and ground-based measurements remained relatively high (R= 0.94), RMBE values were slightly larger in the southernmost part of the desert (27–30° S). Both MODIS-derived data and ground-based measurements show that the albedo at some bright spots in the Atacama Desert may be high enough (up to 0.25 in visible range) for considerably boosting the performance of bifacial photovoltaic technologies (6–12%)., Mathematical Geodesy and Positioning

River ice is a major contributor to flood risk in cold regions due to the physical impediment of flow caused by ice jamming. Although a variety of classifiers have been developed to distinguish ice types using HH or VV intensity of SAR data, mostly based on data from RADARSAT-1 and -2, these classifiers still experience problems with breakup classification, because meltwater development causes overlap in co-polarization backscatter intensities of open water and sheet ice pixels. In this study, we develop a Random Forest classifier based on multiple features of Sentinel-1 data for three main classes generally present during breakup: rubble ice, sheet ice and open water, in a case study over the Athabasca River in Canada. For each ice stage, intensity of the VV and VH backscatter, pseudo-polarimetric decomposition parameters and Grey Level Co-occurrence Matrix texture features were computed for 70 verified sample areas. Several classifiers were developed, based on i) solely intensity features or on ii) a combination of intensity, pseudo-polarimetric and texture features and each classifier was evaluated based on Recursive Feature Elimination with Cross-Validation and pair-wise correlation of the studied features. Results show improved classifier performance when including GLCM mean of VV intensity, and VH intensity features instead of the conventional classifier based solely on intensity. This highlights the complementary nature of texture and intensity for the classification of breaking river ice. GLCM mean incorporates spatial patterns of the co-polarized intensity and sensitivity to context, while VH intensity introduces cross-polarized surface and volume scattering signals and is less sensitive to wind than the commonly used co-polarized intensity. We conclude that the proposed method based on the combination of texture and intensity features is suitable for and performs well in physically complex situations such as breakup, which are hard to classify otherwise. T, Physical and Space Geodesy, Mathematical Geodesy and Positioning

Ice shelf instability is one of the main sources of uncertainty in Antarctica's contribution to future sea level rise. Calving events play a crucial role in ice shelf weakening but remain unpredictable, and their governing processes are still poorly understood. In this study, we analyze the unexpected September 2019 calving event from the Amery Ice Shelf, the largest since 1963 and which occurred almost a decade earlier than expected, to better understand the role of the atmosphere in calving.We find that atmospheric extremes provided a deterministic role in this event. A series of anomalously deep and stationary explosive twin polar cyclones over the Cooperation and Davis seas generated tides and winddriven ocean slope, leading to fracture amplification along the pre-existing rift and ultimately calving of the massive iceberg. The calving was triggered by high oceanward sea surface slopes produced by the storms. The observed recordanomalous atmospheric conditions were promoted by blocking ridges and Antarctic-wide anomalous poleward transport of heat and moisture. Blocking highs helped in (i) directing moist and warm air masses towards the ice shelf and (ii) maintaining the observed extreme cyclones stationary at the front of the ice shelf for several days. Accumulation of cold air over the ice sheet, due to the blocking highs, led to the formation of an intense cold high pressure over the ice sheet, which helped fuel sustained anomalously deep cyclones via increased baroclinicity. Our results stress the importance of atmospheric extremes in ice shelf dynamics via tides and sea surface slope and its need to be accounted for when considering Antarctic ice shelf variability and contribution to sea level, especially given that more of these extremes are predicted under a warmer climate., Mathematical Geodesy and Positioning

Surface meltwater drains on several Antarctic ice shelves, resulting in surface and sub-surface lakes that are potentially critical for the ice shelf collapse. Despite these phenomena, our understanding and assessment of the drainage and refreezing of these lakes is limited, mainly due to lack of field observations and to the limitations of optical satellite imagery during polar night and in cloudy conditions. This paper explores the potential of backscatter intensity and of interferometric coherence and phase from synthetic aperture radar (SAR) imagery as an alternative to assess the dynamics of meltwater lakes. In four case study regions over Amery and Roi Baudouin ice shelves, East Antarctica, we examine spatial and temporal variations in SAR backscatter intensity and interferometric (InSAR) coherence and phase over several lakes derived from Sentinel-1A/B C-band SAR imagery. Throughout the year, the lakes are observed in a completely frozen state, in a partially frozen state with a floating ice lid and as open-water lakes. Our analysis reveals that the meltwater lake delineation is challenging during the melting period when the contrast between melting snow and lakes is indistinguishable. Despite this finding, we show using a combination of backscatter and InSAR observations that lake dynamics can be effectively captured during other non-summertime months. Moreover, our findings highlight the utility of InSAR-based observations for discriminating between refrozen ice and sub-surface meltwater and indicate the potential for phase-based detection and monitoring of rapid meltwater drainage events. The potential of this technique to monitor these meltwater change events is, however, strongly determined by the satellite revisit interval and potential changes in scattering properties due to snowfall or melt events., Mathematical Geodesy and Positioning

Accurately estimating the surface melt volume of the Antarctic Ice Sheet is challenging and has hitherto relied on climate modeling or observations from satellite remote sensing. Each of these methods has its limitations, especially in regions with high surface melt. This study aims to demonstrate the potential of improving surface melt simulations with a regional climate model by deploying a deep learning model. A deep-learning-based framework has been developed to correct surface melt from the regional atmospheric climate model version 2.3p2 (RACMO2), using meteorological observations from automatic weather stations (AWSs) and surface albedo from satellite imagery. The framework includes three steps: (1) training a deep multilayer perceptron (MLP) model using AWS observations, (2) correcting Moderate Resolution Imaging Spectroradiometer (MODIS) albedo observations, and (3) using these two to correct the RACMO2 surface melt simulations. Using observations from three AWSs at the Larsen B and C ice shelves, Antarctica, cross-validation shows a high accuracy (root-mean-square error of 0.95ĝ€¯mmĝ€¯w.e.ĝ€¯d-1, mean absolute error of 0.42ĝ€¯mmĝ€¯w.e.ĝ€¯d-1, and a coefficient of determination of 0.95). Moreover, the deep MLP model outperforms conventional machine learning models and a shallow MLP model. When applying the trained deep MLP model over the entire Larsen Ice Shelf, the resulting corrected RACMO2 surface melt shows a better correlation with the AWS observations for two out of three AWSs. However, for one location (AWS 18), the deep MLP model does not show improved agreement with AWS observations; this is likely because surface melt is largely driven by factors (e.g., air temperature, topography, katabatic wind) other than albedo within the corresponding coarse-resolution model pixels. Our study demonstrates the opportunity to improve surface melt simulations using deep learning combined with satellite albedo observations. However, more work is required to re, Mathematical Geodesy and Positioning

Physically based snow models provide valuable information on snow cover evolution and are therefore key to provide water availability projections. Yet, uncertainties related to snow modelling remain large as a result of differences in the representation of snow physics and meteorological forcing. While many studies focus on evaluating these uncertainties, no snow model comparison has been done in environments where sublimation is the main ablation process. This study evaluates a case study in the semi-arid Andes of Chile and aims to compare two snow models with different complexities, SNOWPACK and SnowModel, at a local point over one snow season and to evaluate their sensitivity relative to parameterisation and forcing. For that purpose, the two models are forced with (i) the most ideal set of input parameters, (ii) an ensemble of different physical parameterisations, and (iii) an ensemble of biased forcing. Results indicate large uncertainties depending on forcing, the snow roughness length z0, albedo parameterisation, and fresh snow density parameterisation. The uncertainty caused by the forcing is directly related to the bias chosen. Even though the models show significant differences in their physical complexity, the snow model choice is of least importance, as the sensitivity of both models to the forcing data was on the same order of magnitude and highly influenced by the precipitation uncertainties. The sublimation ratio ranges are in agreement for the two models: 36.4 % to 80.7 % for SnowModel and 36.3 % to 86.0 % for SNOWPACK, and are related to the albedo parameterisation and snow roughness length choice for the two models., Mathematical Geodesy and Positioning

Surface albedo is an important forcing parameter that drives the radiative energy budget as it determines the fraction of the downwelling solar irradiance that the surface reflects. Here we report on ground-based measurements of the spectral albedo (350–2200 nm) carried out at 20 sites across a North–South transect of approximately 1300 km in the Atacama Desert, from latitude 18° S to latitude 30° S. These spectral measurements were used to evaluate remote sensing estimates of the albedo derived from the Moderate Resolution Imaging Spectroradiometer (MODIS). We found that the relative mean bias error (RMBE) of MODIS-derived estimates was within ± 5% of ground-based measurements in most of the Atacama Desert (18–27° S). Although the correlation between MODIS-derived estimates and ground-based measurements remained relatively high (R= 0.94), RMBE values were slightly larger in the southernmost part of the desert (27–30° S). Both MODIS-derived data and ground-based measurements show that the albedo at some bright spots in the Atacama Desert may be high enough (up to 0.25 in visible range) for considerably boosting the performance of bifacial photovoltaic technologies (6–12%)., Mathematical Geodesy and Positioning

Ice shelf instability is one of the main sources of uncertainty in Antarctica's contribution to future sea level rise. Calving events play a crucial role in ice shelf weakening but remain unpredictable, and their governing processes are still poorly understood. In this study, we analyze the unexpected September 2019 calving event from the Amery Ice Shelf, the largest since 1963 and which occurred almost a decade earlier than expected, to better understand the role of the atmosphere in calving.We find that atmospheric extremes provided a deterministic role in this event. A series of anomalously deep and stationary explosive twin polar cyclones over the Cooperation and Davis seas generated tides and winddriven ocean slope, leading to fracture amplification along the pre-existing rift and ultimately calving of the massive iceberg. The calving was triggered by high oceanward sea surface slopes produced by the storms. The observed recordanomalous atmospheric conditions were promoted by blocking ridges and Antarctic-wide anomalous poleward transport of heat and moisture. Blocking highs helped in (i) directing moist and warm air masses towards the ice shelf and (ii) maintaining the observed extreme cyclones stationary at the front of the ice shelf for several days. Accumulation of cold air over the ice sheet, due to the blocking highs, led to the formation of an intense cold high pressure over the ice sheet, which helped fuel sustained anomalously deep cyclones via increased baroclinicity. Our results stress the importance of atmospheric extremes in ice shelf dynamics via tides and sea surface slope and its need to be accounted for when considering Antarctic ice shelf variability and contribution to sea level, especially given that more of these extremes are predicted under a warmer climate., Mathematical Geodesy and Positioning

River ice is a major contributor to flood risk in cold regions due to the physical impediment of flow caused by ice jamming. Although a variety of classifiers have been developed to distinguish ice types using HH or VV intensity of SAR data, mostly based on data from RADARSAT-1 and -2, these classifiers still experience problems with breakup classification, because meltwater development causes overlap in co-polarization backscatter intensities of open water and sheet ice pixels. In this study, we develop a Random Forest classifier based on multiple features of Sentinel-1 data for three main classes generally present during breakup: rubble ice, sheet ice and open water, in a case study over the Athabasca River in Canada. For each ice stage, intensity of the VV and VH backscatter, pseudo-polarimetric decomposition parameters and Grey Level Co-occurrence Matrix texture features were computed for 70 verified sample areas. Several classifiers were developed, based on i) solely intensity features or on ii) a combination of intensity, pseudo-polarimetric and texture features and each classifier was evaluated based on Recursive Feature Elimination with Cross-Validation and pair-wise correlation of the studied features. Results show improved classifier performance when including GLCM mean of VV intensity, and VH intensity features instead of the conventional classifier based solely on intensity. This highlights the complementary nature of texture and intensity for the classification of breaking river ice. GLCM mean incorporates spatial patterns of the co-polarized intensity and sensitivity to context, while VH intensity introduces cross-polarized surface and volume scattering signals and is less sensitive to wind than the commonly used co-polarized intensity. We conclude that the proposed method based on the combination of texture and intensity features is suitable for and performs well in physically complex situations such as breakup, which are hard to classify otherwise. T, Physical and Space Geodesy, Mathematical Geodesy and Positioning

Mass loss from the Greenland Ice Sheet (GrIS) has accelerated over the past two decades, coincident with rapid Arctic warming and increasing moisture transport over Greenland by atmospheric rivers (ARs). Summer ARs affecting western Greenland trigger GrIS melt events, but the physical mechanisms through which ARs induce melt are not well understood. This study elucidates the coupled surface-atmosphere processes by which ARs force GrIS melt through analysis of the surface energy balance (SEB), cloud properties, and local- to synoptic-scale atmospheric conditions during strong summer AR events affecting western Greenland. ARs are identified in MERRA-2 reanalysis (1980-2017) and classified by integrated water vapor transport (IVT) intensity. SEB, cloud, and atmospheric data from regional climate model, observational, reanalysis, and satellite-based datasets are used to analyze melt-inducing physical processes during strong,.90th percentile ''AR901'' events. Near AR ''landfall,'' AR901 days feature increased cloud cover that reduces net shortwave radiation and increases net longwave radiation. As these oppositely signed radiative anomalies partly cancel during AR901 events, increased melt energy in the ablation zone is primarily provided by turbulent heat fluxes, particularly sensible heat flux. These turbulent heat fluxes are driven by enhanced barrier winds generated by a stronger synoptic pressure gradient combined with an enhanced local temperature contrast between cool over-ice air and the anomalously warm surrounding atmosphere. During AR901 events in northwest Greenland, anomalous melt is forced remotely through a clear-sky foehn regime produced by downslope flow in eastern Greenland., Mathematical Geodesy and Positioning

Snow and ice albedo schemes in present-day climate models often lack a sophisticated radiation penetration scheme and do not explicitly include spectral albedo variations. In this study, we evaluate a new snow albedo scheme in the Regional Atmospheric Climate Model (RACMO2) for the Greenland ice sheet, version 2.3p3, that includes these processes. The new albedo scheme uses the TwostreAm Radiative TransfEr in Snow (TARTES) model and the Spectral-To-NarrOWBand ALbedo (SNOWBAL) module, version 1.2. Additionally, the bare-ice albedo parameterization has been updated. The snow and ice broadband and narrowband albedo output of the updated version of RACMO2 is evaluated using the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) and Kangerlussuaq transect (K-Transect) in situ data and Moderate Resolution Imaging Spectroradiometer (MODIS) remote-sensing observations. Generally, the modeled narrowband and broadband albedo is in very good agreement with satellite observations, leading to a negligible domain-Averaged broadband albedo bias for the interior. Some discrepancies are, however, observed close to the ice margin. Compared to the previous model version, RACMO2.3p2, the broadband albedo is considerably higher in the bare-ice zone during the ablation season, as atmospheric conditions now alter the bare-ice broadband albedo. For most other regions, however, the updated broadband albedo is lower due to spectral effects, radiation penetration or enhanced snow metamorphism., Mathematical Geodesy and Positioning

Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are among the fastest changing outlet glaciers in West Antarctica with large consequences for global sea level. Yet, assessing how much and how fast both glaciers will weaken if these changes continue remains a major uncertainty as many of the processes that control their ice shelf weakening and grounding line retreat are not well understood. Here, we combine multisource satellite imagery with modeling to uncover the rapid development of damage areas in the shear zones of Pine Island and Thwaites ice shelves. These damage areas consist of highly crevassed areas and open fractures and are first signs that the shear zones of both ice shelves have structurally weakened over the past decade. Idealized model results reveal moreover that the damage initiates a feedback process where initial ice shelf weakening triggers the development of damage in their shear zones, which results in further speedup, shearing, and weakening, hence promoting additional damage development. This damage feedback potentially preconditions these ice shelves for disintegration and enhances grounding line retreat. The results of this study suggest that damage feedback processes are key to future ice shelf stability, grounding line retreat, and sea level contributions from Antarctica. Moreover, they underline the need for incorporating these feedback processes, which are currently not accounted for in most ice sheet models, to improve sea level rise projections., Mathematical Geodesy and Positioning, Physical and Space Geodesy

Compared to other Arctic ice masses, Svalbard glaciers are low-elevated with flat interior accumulation areas, resulting in a marked peak in their current hypsometry (area-elevation distribution) at ~450 m above sea level. Since summer melt consistently exceeds winter snowfall, these low-lying glaciers can only survive by refreezing a considerable fraction of surface melt and rain in the porous firn layer covering their accumulation zones. We use a high-resolution climate model to show that modest atmospheric warming in the mid-1980s forced the firn zone to retreat upward by ~100 m to coincide with the hypsometry peak. This led to a rapid areal reduction of firn cover available for refreezing, and strongly increased runoff from dark, bare ice areas, amplifying mass loss from all elevations. As the firn line fluctuates around the hypsometry peak in the current climate, Svalbard glaciers will continue to lose mass and show high sensitivity to temperature perturbations., Mathematical Geodesy and Positioning, Physical and Space Geodesy

Between 1992 and 2017, the Antarctic Ice Sheet (AIS) lost ice equivalent to 7.6 ± 3.9 mm of sea level rise. AIS mass loss is mitigated by ice shelves that provide a buttress by regulating ice flow from tributary glaciers. However, ice-shelf stability is threatened by meltwater ponding, which may initiate, or reactivate preexisting, fractures, currently poorly understood processes. Here, through ground penetrating radar (GPR) analysis over a buried lake in the grounding zone of an East Antarctic ice shelf, we present the first field observations of a lake drainage event in Antarctica via vertical fractures. Concurrent with the lake drainage event, we observe a decrease in surface elevation and an increase in Sentinel-1 backscatter. Finally, we suggest that fractures that are initiated or reactivated by lake drainage events in a grounding zone will propagate with ice flow onto the ice shelf itself, where they may have implications for its stability., Mathematical Geodesy and Positioning

About 20 % of all snow accumulation in Antarctica occurs on the ice shelves. There, ice rises control the spatial surface mass balance (SMB) distribution by inducing snowfall variability and wind erosion due to their topography. Moreover these ice rises buttress the ice flow and represent ideal drilling locations for ice cores. In this study we assess the connection between snowfall variability and wind erosion to provide a better understanding of how ice rises impact SMB variability, how well this is captured in the regional atmospheric climate model RACMO2 and the implications of this SMB variability for ice rises as an ice core drilling site. By combining ground-penetrating radar (GPR) profiles from two ice rises in Dronning Maud Land with ice core dating, we reconstruct spatial and temporal SMB variations from 1983 to 2018 and compare the observed SMB with output from RACMO2 and SnowModel. Our results show snowfall-driven differences of up to 1.5 times higher SMB on the windward side of both ice rises than on the leeward side as well as a local erosion-driven minimum at the ice divide of the ice rises. RACMO2 captures the snowfall-driven differences but overestimates their magnitude, whereas the erosion on the peak can be reproduced by SnowModel with RACMO2 forcing. Observed temporal variability of the average SMBs, retrieved from the GPR data for four time intervals in the 1983-2018 range, are low at the peak of the easternmost ice rise (∼ 0.06 m w.e. yr−1), while they are higher (∼ 0.09 m w.e. yr−1) on the windward side of the ice rise. This implies that at the peak of the ice rise, higher snowfall, driven by orographic uplift, is balanced out by local erosion. As a consequence of this, the SMB recovered from the ice core matches the SMB from the GPR at the peak of the ice rise but not at the windward side of the ice rise, suggesting that the SMB signal is damped in the ice core., Mathematical Geodesy and Positioning

Continuous measurements of blowing snow are scarce, both in time and space. Satellites now provide the opportunity to derive blowing snow occurrences, transport and sublimation rates over Antarctica. These products are extremely valuable and offer a continental-wide assessment of blowing snow, which is an important but unknown component of the surface mass balance of the Antarctic ice sheet. However, little ground truth is available to validate these retrievals. The recent application of ceilometers for detection of blowing snow frequencies provides an opportunity to validate the satellite retrievals of blowing snow. A routine to detect blowing snow occurrence from ground-based remote sensing ceilometers has been developed at two coastal locations in East Antarctica for the 2011–2016 time period. Thanks to their ground-based location, ceilometers are able to detect blowing snow events in the presence of clouds and precipitation, which can be missed by the satellite, since optically thick clouds impede the penetration of the signal. In coastal areas, more than 90% of blowing snow occurs under cloudy conditions and represent 30 to 56% of all cloudy conditions at Princess Elisabeth and Neumayer III (Neumayer hereafter) stations, respectively. For cloud-free conditions, 8% of the measurements at Princess Elisabeth (and none at Neumayer) are identified as blowing snow by the satellite but not by the ceilometer, likely due to differences in sensors, limitation of the surface identification by the satellite, or the spatial inhomogeneity of the blowing snow event. While the satellite blowing snow retrieval is a useful product, further investigation is needed to reduce the uncertainties on blowing snow frequencies associated with clouds., Mathematical Geodesy and Positioning

To better understand and quantify the impact of clouds on the Greenland Ice Sheet surface mass balance (SMB), we study the spatiotemporal variability of the cloud radiative effect (CRE). The total CRE is separated in short-term and long-term impacts by performing multiple simulations with the SNOWPACK model for 2001-+2010. The annual total CRE is 16.8 ± 4.5 W m−2, reducing the SMB with −157 ± 3.8 Gt yr−1. Summer cloud radiative cooling is −6.4 ± 5.7 W m−2 in the ablation area, increasing the SMB with 121 ± 2.2 Gt yr−1. The annual integrated impact is cloud-reduced SMB of −36 Gt yr−1. The short-term effect dominates the opposing long-term effects through the albedo-melt feedback. A long-term warming effect decreases the albedo and so preconditions the surface for enhanced (summer) melt. The impact of the CRE, determined by spatial, temporal and initial conditions, explains existing conflicted views on the role of cloud radiation and emphasizes the need for accurate cloud and albedo representations in future studies., Mathematical Geodesy and Positioning

The response of the Greenland Ice Sheet (GrIS) to a warmer climate is uncertain on long time scales. Climate models, such as those participating in the Coupled Model Intercomparison Project phase 6 (CMIP6), are used to assess this uncertainty. The Community Earth System Model version 2.1 (CESM2) is a CMIP6 model capable of running climate simulations with either one-way coupling (fixed ice sheet geometry) or two-way coupling (dynamic geometry) to the GrIS. The model features prognostic snow albedo, online downscaling using elevation classes, and a firn pack to refreeze percolating melt water. Here we evaluate the representation of the GrIS surface energy balance and surface mass balance in CESM2 at 1° resolution with fixed GrIS geometry. CESM2 agrees closely with ERA-Interim reanalysis data for key controls on GrIS SMB: surface pressure, sea ice extent, 500 hPa geopotential height, wind speed, and 700 hPa air temperature. Cloudsat-CALIPSO data show that supercooled liquid-containing clouds are adequately represented, whereas comparisons to Moderate Resolution Imaging Spectroradiometer and CM SAF Cloud, Albedo, and Surface Radiation data set from Advanced Very High Resolution Radiometer data second edition data suggest that CESM2 underestimates surface albedo. The seasonal cycle and spatial patterns of surface energy balance and surface mass balance components in CESM2 agree well with regional climate model RACMO2.3p2, with GrIS-integrated melt, refreezing, and runoff bracketed by RACMO2 counterparts at 11 and 1 km. Time series of melt, runoff, and SMB show a break point around 1990, similar to RACMO2. These results suggest that GrIS SMB is realistic in CESM2, which adds confidence to coupled ice sheet-climate experiments that aim to assess the GrIS contribution to future sea level rise., Mathematical Geodesy and Positioning, Physical and Space Geodesy

In the semiarid Andes of Chile, farmers and industry in the cordillera lowlands depend on water from snowmelt, as annual rainfall is insufficient to meet their needs. Despite the importance of snow cover for water resources in this region, understanding of snow depth distribution and snow mass balance is limited. Whilst the effect of wind on snow cover pattern distribution has been assessed, the relative importance of melt versus sublimation has only been studied at the point scale over one catchment. Analyzing relative ablation rates and evaluating uncertainties are critical for understanding snow depth sensitivity to variations in climate and simulating the evolution of the snowpack over a larger area and over time. Using a distributed snowpack model (SnowModel), this study aims to simulate melt and sublimation rates over the instrumented watershed of La Laguna (513 km2, 3150-5630 a.s.l., 30° S, 70° W), during two hydrologically contrasting years (i.e., dry vs. wet). The model is calibrated and forced with meteorological data from nine Automatic Weather Stations (AWSs) located in the watershed and atmospheric simulation outputs from the Weather Research and Forecasting (WRF) model. Results of simulations indicate first a large uncertainty in sublimation-to-melt ratios depending on the forcing as the WRF data have a cold bias and overestimate precipitation in this region. These input differences cause a doubling of the sublimation-to-melt ratio using WRF forcing inputs compared to AWS. Therefore, the use of WRF model output in such environments must be carefully adjusted so as to reduce errors caused by inherent bias in the model data. For both input datasets, the simulations indicate a similar sublimation fraction for both study years, but ratios of sublimation to melt vary with elevation as melt rates decrease with elevation due to decreasing temperatures. Finally results indicate that snow persistence during the spring period decreases the ratio of, Mathematical Geodesy and Positioning

Snow and ice albedo schemes in present-day climate models often lack a sophisticated radiation penetration scheme and do not explicitly include spectral albedo variations. In this study, we evaluate a new snow albedo scheme in the Regional Atmospheric Climate Model (RACMO2) for the Greenland ice sheet, version 2.3p3, that includes these processes. The new albedo scheme uses the TwostreAm Radiative TransfEr in Snow (TARTES) model and the Spectral-To-NarrOWBand ALbedo (SNOWBAL) module, version 1.2. Additionally, the bare-ice albedo parameterization has been updated. The snow and ice broadband and narrowband albedo output of the updated version of RACMO2 is evaluated using the Programme for Monitoring of the Greenland Ice Sheet (PROMICE) and Kangerlussuaq transect (K-Transect) in situ data and Moderate Resolution Imaging Spectroradiometer (MODIS) remote-sensing observations. Generally, the modeled narrowband and broadband albedo is in very good agreement with satellite observations, leading to a negligible domain-Averaged broadband albedo bias for the interior. Some discrepancies are, however, observed close to the ice margin. Compared to the previous model version, RACMO2.3p2, the broadband albedo is considerably higher in the bare-ice zone during the ablation season, as atmospheric conditions now alter the bare-ice broadband albedo. For most other regions, however, the updated broadband albedo is lower due to spectral effects, radiation penetration or enhanced snow metamorphism., Mathematical Geodesy and Positioning

Mass loss from the Greenland Ice Sheet (GrIS) has accelerated over the past two decades, coincident with rapid Arctic warming and increasing moisture transport over Greenland by atmospheric rivers (ARs). Summer ARs affecting western Greenland trigger GrIS melt events, but the physical mechanisms through which ARs induce melt are not well understood. This study elucidates the coupled surface-atmosphere processes by which ARs force GrIS melt through analysis of the surface energy balance (SEB), cloud properties, and local- to synoptic-scale atmospheric conditions during strong summer AR events affecting western Greenland. ARs are identified in MERRA-2 reanalysis (1980-2017) and classified by integrated water vapor transport (IVT) intensity. SEB, cloud, and atmospheric data from regional climate model, observational, reanalysis, and satellite-based datasets are used to analyze melt-inducing physical processes during strong,.90th percentile ''AR901'' events. Near AR ''landfall,'' AR901 days feature increased cloud cover that reduces net shortwave radiation and increases net longwave radiation. As these oppositely signed radiative anomalies partly cancel during AR901 events, increased melt energy in the ablation zone is primarily provided by turbulent heat fluxes, particularly sensible heat flux. These turbulent heat fluxes are driven by enhanced barrier winds generated by a stronger synoptic pressure gradient combined with an enhanced local temperature contrast between cool over-ice air and the anomalously warm surrounding atmosphere. During AR901 events in northwest Greenland, anomalous melt is forced remotely through a clear-sky foehn regime produced by downslope flow in eastern Greenland., Mathematical Geodesy and Positioning

Compared to other Arctic ice masses, Svalbard glaciers are low-elevated with flat interior accumulation areas, resulting in a marked peak in their current hypsometry (area-elevation distribution) at ~450 m above sea level. Since summer melt consistently exceeds winter snowfall, these low-lying glaciers can only survive by refreezing a considerable fraction of surface melt and rain in the porous firn layer covering their accumulation zones. We use a high-resolution climate model to show that modest atmospheric warming in the mid-1980s forced the firn zone to retreat upward by ~100 m to coincide with the hypsometry peak. This led to a rapid areal reduction of firn cover available for refreezing, and strongly increased runoff from dark, bare ice areas, amplifying mass loss from all elevations. As the firn line fluctuates around the hypsometry peak in the current climate, Svalbard glaciers will continue to lose mass and show high sensitivity to temperature perturbations., Mathematical Geodesy and Positioning, Physical and Space Geodesy

Pine Island Glacier and Thwaites Glacier in the Amundsen Sea Embayment are among the fastest changing outlet glaciers in West Antarctica with large consequences for global sea level. Yet, assessing how much and how fast both glaciers will weaken if these changes continue remains a major uncertainty as many of the processes that control their ice shelf weakening and grounding line retreat are not well understood. Here, we combine multisource satellite imagery with modeling to uncover the rapid development of damage areas in the shear zones of Pine Island and Thwaites ice shelves. These damage areas consist of highly crevassed areas and open fractures and are first signs that the shear zones of both ice shelves have structurally weakened over the past decade. Idealized model results reveal moreover that the damage initiates a feedback process where initial ice shelf weakening triggers the development of damage in their shear zones, which results in further speedup, shearing, and weakening, hence promoting additional damage development. This damage feedback potentially preconditions these ice shelves for disintegration and enhances grounding line retreat. The results of this study suggest that damage feedback processes are key to future ice shelf stability, grounding line retreat, and sea level contributions from Antarctica. Moreover, they underline the need for incorporating these feedback processes, which are currently not accounted for in most ice sheet models, to improve sea level rise projections., Mathematical Geodesy and Positioning, Physical and Space Geodesy