Ice Dynamics Research Articles

A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling

Abstract. The ability of numerical sea ice models to reproduce localized deformation features associated with fracture processes is key for an accurate representation of the ice dynamics and of dynamically coupled physical processes in the Arctic and Antarctic. Equally key is the capacity of these models to minimize the numerical diffusion stemming from the advection of these features to ensure that the associated strong gradients persist in time, without the need to unphysically re-inject energy for re-localization. To control diffusion and improve the approximation quality, we present a new numerical core for the dynamics of sea ice that is based on higher-order finite-element discretizations for the momentum equation and higher-order discontinuous Galerkin methods for the advection. The mathematical properties of this core are discussed, and a detailed description of an efficient shared-memory parallel implementation is given. In addition, we present different numerical tests and apply the new framework to a benchmark problem to quantify the advantages of the higher-order discretization. These tests are based on Hibler's viscous–plastic sea ice model, but the implementation of the developed framework in the context of other physical models reproducing a strong localization of the deformation is possible.

Prediction of nonlinear interface dynamics in the unidirectional freezing of particle suspensions with a rigid compacted layer

Ice growth in aqueous suspensions widely exists in natural and industrial settings where freezing of water occurs with a porous medium and results in the spatiotemporal evolution of a solid/liquid interface. Physical models have been proposed in previous efforts to describe the dynamic behaviours of interfaces in the unidirectional freezing of particle suspensions. Most previous models dealt with this process at the single-particle level, and the growth of a particle compacted layer coupled with ice growth was hardly addressed. Here, based on the fundamental momentum theorem, we propose a consistent theoretical framework to address the unidirectional freezing process in particle suspensions free of solute, which focuses exclusively on the effect of water permeation. We first propose a constant compacted layer model, in which an interface undercooling-dependent pushing force exerted on the compacted layer is derived based on surface tension. Subsequently, a dynamically growing compacted layer is established and analysed. Numerical solutions of the nonlinear models reveal the dependence of the system dynamics on some typical physical parameters, such as particle radius, initial particle concentration in the suspensions, and freezing velocity. The system dynamics are characterized by interface velocity, interface undercooling, and interface recoil as functions of time. The models allow us to reconsider the dynamic ice growth during the freezing of particle suspensions in a simple but novel way, with potential implications for both understanding and controlling not only ice formation in porous media but also crystallization processes in other complex systems.

A High Arctic inner shelf–fjord system from the Last Glacial Maximum to the present: Bessel Fjord and southwest Dove Bugt, northeastern Greenland

Abstract. The Greenland Ice Sheet (GrIS) responds rapidly to the present climate; therefore, its response to the predicted future warming is of concern. To learn more about the impact of future climatic warming on the ice sheet, decoding its behavior during past periods of warmer than present climate is important. However, due to the scarcity of marine studies reconstructing ice sheet conditions on the Northeast Greenland shelf and adjacent fjords, the timing of the deglaciation over marine regions and its connection to forcing factors remain poorly constrained. This includes data collected in fjords that encompass the Holocene thermal maximum (HTM), a period in which the climate was warmer than it is at present. This paper aims to use new bathymetric data and the analysis of sediment gravity cores to enhance our understanding of ice dynamics of the GrIS in a fjord and inner shelf environment as well as give insight into the timing of deglaciation and provide a palaeoenvironmental reconstruction of southwestern Dove Bugt and Bessel Fjord since the Last Glacial Maximum (LGM). North–south-oriented glacial lineations and the absence of pronounced moraines in southwest Dove Bugt, an inner continental shelf embayment (trough), suggest the southwards and offshore flow of Storstrømmen, the southern branch of the Northeast Greenland Ice Stream (NEGIS). Sedimentological data suggest that an ice body, theorized to be the NEGIS, may have retreated from the region slightly before ∼ 11.4 cal ka BP. The seabed morphology of Bessel Fjord, a fjord terminating in southern Dove Bugt, includes numerous basins separated by thresholds. The position of basin thresholds, which include some recessional moraines, suggest that the GrIS had undergone multiple halts or readvances during deglaciation, likely during one of the cold events identified in the Greenland Summit temperature records. A minimum age of 7.1 cal ka BP is proposed for the retreat of ice through the fjord to or west of its present-day position in the Bessel Fjord catchment area. This suggests that the GrIS retreated from the marine realm in Early Holocene, around the onset of the HTM in this region, a period when the mean July temperature was at least 2–3 ∘C higher than at present and remained at or west of this onshore position for the remainder of the Holocene. The transition from predominantly mud to muddy sand layers in a mid-fjord core at ∼ 4 cal ka BP may be the result of increased sediment input from nearby and growing ice caps. This shift may suggest that in the Late Holocene (Meghalayan), a period characterized by a temperature drop to modern values, ice caps in Bessel Fjord probably fluctuated with greater sensitivity to climatic conditions than the northeastern sector of the GrIS.

Satellite record reveals 1960s acceleration of Totten Ice Shelf in East Antarctica

Wilkes Land and Totten Glacier (TG) in East Antarctica (EA) have been losing ice mass significantly since 1989. There is a lack of knowledge of long-term mass balance in the region which hinders the estimation of its contribution to global sea level rise. Here we show that this acceleration trend in TG has occurred since the 1960s. We reconstruct ice flow velocity fields of 1963–1989 in TG from the first-generation satellite images of ARGON and Landsat-1&4, and build a five decade-long record of ice dynamics. We find a persistent long-term ice discharge rate of 68 ± 1 Gt/y and an acceleration of 0.17 ± 0.02 Gt/y2 from 1963 to 2018, making TG the greatest contributor to global sea level rise in EA. We attribute the long-term acceleration near grounding line from 1963 to 2018 to basal melting likely induced by warm modified Circumpolar Deep Water. The speed up in shelf front during 1973–1989 was caused by a large calving front retreat. As the current trend continues, intensified monitoring in the TG region is recommended in the next decades.

DynIceData: a gridded ice–water classification dataset at short-time intervals based on observations from multiple satellites over the marginal ice zone

ABSTRACT High-resolution observations of short-term changes in sea ice are critical to understanding ice dynamics and also provide important information used in advice to shipping, especially in the Arctic. Although individual satellite sensors provide periodic sea ice observations with spatial resolutions of tens of meters, information regarding changes that occur over short time intervals of minutes or hours is limited. In this study, a gridded ice–water classification dataset with a high temporal resolution was developed based on observations acquired by multiple satellite sensors in the Marginal Ice Zone (MIZ). This dataset – DynIceData – which combines Sentinel-1 Synthetic Aperture Radar (SAR) data with Gaofen-3 (GF-3) SAR and SDGSAT-1 thermal infrared imagery was used to obtain observations of the MIZ with a range of temporal resolutions ranging from minutes to tens of hours. The areas of the Arctic covered include the Kara Sea, Beaufort Sea, and Greenland Sea during the period from August 2021 to August 2022. Object-oriented segmentation and thresholding were used to obtain the ice–water classification map from Sentinel-1 and GF-3 SAR image pairs and Sentinel-1 SAR and SDGSAT-1 thermal image pairs. The time interval between the images in each pair ranged from 1 minute to 68 hours. Ten-kilometer grid sample granules with a spatial resolution of 25 m for the GF-3 SAR data and 30 m for the SDGSAT-1 thermal data were used. The classification was verified as having an overall accuracy of at least 95.58%. The DynIceData dataset consists of 7338 samples, which could be used as reference data for further research on rapid changes in sea ice patterns at different short time scales and provide support for research on thermodynamic and dynamic models of sea ice in combination with other environmental data, thus potentially improving the accuracy of sea ice forecasting using Artificial Intelligence. The dataset can be accessed at https://doi.org/10.57760/sciencedb.j00001.00784.

Dynamics of water temperature and interfacial heat flux in lake ice and under-ice layers

Abstract The purpose of this study is to understand the temperature dynamics of ice and lake water under ice during ice growth. A vertical temperature monitoring system for ice and water under ice was employed for the first time in Wuliangsuhai Lake from 20 December 2016 to 15 March 2017 to obtain continuous monitoring temperature data necessary for calculating the heat flux density at the ice–water interface. During the ice growth period, the ice body temperature and the water body temperature under the ice show an obvious rising trend from the upper layer to the lower layer. As the ice body enters the melting period, the ice temperature and the water temperature under the ice begin to change rapidly. In the presence of snow over the ice, the fluctuations of ice temperature and water temperature under the ice remain stable. According to the experimental data, the analytical equation derived from the Stefan condition is used to calculate the heat flux at the intersection boundary of ice body and water body under ice, which varies from 1.73 to 19.6 W/m2. The dynamic change of ice thickness is significantly influenced by the temperature of ice and water and the heat flux density at the ice–water interface.

Polarimetric airborne scientific instrument, mark 2, an ice‐sounding airborne synthetic aperture radar for subglacial 3D imagery

AbstractPolarimetric Airborne Scientific INstrument, mark 2 (PASIN2) is a 150 MHz coherent pulsed radar with the purpose of deep ice sounding for bedrock, subglacial channels and ice‐water interface detection in Antarctica. It is designed and operated by the British Antarctic Survey from 2014. With multiple antennas, oriented along and across‐track, for transmission and reception, it enables polarimetric 3D estimation of the ice base with a single pass, reducing the gridding density of the survey paths. The off‐line data processing stream consists of channel calibration; 2D synthetic aperture radar (SAR) imaging based on back‐projection, for along‐track and range dimensions; and finally, a direction of arrival estimation (DoA) of the remaining across‐track angle, by modifying the non‐linear MUSIC algorithm. Calibration flights, during the Antarctic Summer campaigns in 16/17 and 19/20 seasons, assessed and validated the instrument and processing performances. Imaging flights over ice streams and ice shelves close to grounding lines demonstrate the 3D sensing capabilities. By resolving directional ambiguities and accounting for reflector across‐track location, the true ice thickness and bed elevation are obtained, thereby removing the error of the usual assumption of vertical DoA, that greatly influence the output of flow models of ice dynamics.

Sea ice drift tracks from autonomous buoys in the MOSAiC Distributed Network

A network of autonomous, ice-tethered buoys was deployed around the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) experiment in late September 2019 for a year-long drift in the Arctic Transpolar Drift Stream. The buoys were deployed as part of the MOSAiC distributed network (DN) which included 12 multi-instrumented ice stations and an additional 116 GPS buoys distributed primarily within a 40 km radius of the MOSAiC Central Observatory. Buoy coverage within the DN was maintained with additional deployments throughout the year-long drift allowing for collection of data over a full sea ice growth and melt cycle. All GPS position data from buoys deployed within the DN have been assembled and processed into the collection of 216 quality-controlled buoy drift tracks presented in this dataset covering the period 26 September 2019 – 23 May 2021. The drift tracks in this collection are ideal for studies of dynamic sea ice motion around the MOSAiC experiment at cascading spatial scales ranging from 100s of meters to 100s of km.

Development of a dynamic ice penetration model for runway de-icing products through performance experimental characterization

During winter, the clear runway concept is crucial to airports for maintaining passenger and employee safety. Research can optimize de-icing product application rates, thereby reducing costs and environmental impact. Currently, three aerospace standards exist for evaluating the performance of runway de-icing products. These comparative tests—ice melting, ice undercutting, and ice penetration tests—are carried out jointly in a systematic way and indicate whether a runway de-icing product is effective in actively reducing the presence of ice. Here we present a numerical model for the ice penetration method, for fluid runway de-icing product, with the goal of deepening the knowledge in this field. We apply a series of laboratory-based tests to measure the performance of runway de-icing products, and develop a new model to analyze the behavior of each product. The principal observation is that the product is more effective at 1 min and effectiveness decreases following an inverse logarithmic curve trend to 30 min. Through modeling, we observe that the dominant mode of heat transfer is natural laminar convection, and the model predicts a uniform melt length for the entire tube section. Thus, simplified ice penetration modeling using AS6211 can acceptably predict the rate of ice penetration as a function of time (mm/min), with the modeling confirmed by our experimental results. The future application of our simplified model to other existing de-icing fluid characterization test methods could identify the effects of diluted or undiluted solid runway de-icing products on ice penetration rates and ultimately on spreading rates.

SIMUG – finite element model of sea ice dynamics on triangular grid in local Cartesian basis

Abstract The paper presents the dynamical core of the new sea ice model SIMUG (Sea Ice Model on Unstructured Grid) on the A- and CD-types of unstructured triangular grids in the local-element basis on sphere. Three standardized box tests to reproduce the Linear Kinematic Features (LKFs), and the short-term forecast in the real Arctic Ocean geometry with the realistic atmosphere and ocean forcing demonstrate the model quality compared to other sea ice models like CICE, FESOM, MITgcm, and ICON-O. The distinctive features of the model presented are a wide choice of transport schemes, and the new numerical implementation with the serial and parallel C++ coding and INMOST, Ani2D, and Ani3D packages to deal with unstructured grids. Code profiling and scalability assessment are carried out. In general, the A-version of the ice drift model works faster, but has fewer degrees of freedom on the same grid. Due to the increase in the degrees of freedom, the model on the CD grid gives ultra-resolution of LKFs, but requires more strict conditions for stability.

Melting domain size and recrystallization dynamics of ice revealed by time-resolved x-ray scattering

The phase transition between water and ice is ubiquitous and one of the most important phenomena in nature. Here, we performed time-resolved x-ray scattering experiments capturing the melting and recrystallization dynamics of ice. The ultrafast heating of ice I is induced by an IR laser pulse and probed with an intense x-ray pulse which provided us with direct structural information on different length scales. From the wide-angle x-ray scattering (WAXS) patterns, the molten fraction, as well as the corresponding temperature at each delay, were determined. The small-angle x-ray scattering (SAXS) patterns, together with the information extracted from the WAXS analysis, provided the time-dependent change of the size and the number of liquid domains. The results show partial melting (~13%) and superheating of ice occurring at around 20 ns. After 100 ns, the average size of the liquid domains grows from about 2.5 nm to 4.5 nm by the coalescence of approximately six adjacent domains. Subsequently, we capture the recrystallization of the liquid domains, which occurs on microsecond timescales due to the cooling by heat dissipation and results to a decrease of the average liquid domain size.

Stability and energy consumption analysis of arctic fleet: modeling and simulation based on future motion of multi-ship.

Ensuring the safety of Arctic shipping and preserving the Arctic ecological environment are emerging as key challenges in the shipping sector. Ship collisions and getting trapped in ice are frequently occurring under dynamic ice conditions due to the Arctic environment, making research on ship navigation in Arctic routes significant. Leveraging ship networking technology, we developed an intelligent microscopic model which considered factors such as the future motion trends of multi-ships in front and the influence of pack ice, and carried out a stability analysis of the model utilizing linear and nonlinear methods. Additionally, the accuracy of the theoretical results was further validated through simulation experiments with diverse scenarios. The conclusions manifest that the model can magnify the anti-disturbance ability of traffic flow. Additionally, the problem of energy consumption due to ship speed is explored, and it is determined that the model has a positive intention in reducing speed fluctuations and ship energy consumption. This paper highlights the potential of intelligent microscopic models in studying the safety and sustainability of Arctic shipping routes, providing targeted initiatives to improve safety, efficiency, and sustainability in Arctic shipping.

Ice Velocity Variations of the Cook Ice Shelf, East Antarctica, from 2017 to 2022 from Sentinel-1 SAR Time-Series Offset Tracking

The Cook Glacier drains a significant portion of the Wilkes Subglacial Basin, the largest subglacial basin in East Antarctica—which feeds the Cook Ice Shelf. The ice velocity of the Cook Ice Shelf needs to be monitored precisely and accurately, as it plays a critical role in determining the ice discharge from the Wilkes Subglacial Basin. In this study, we measured the annual ice velocities of the Cook Ice Shelf using the offset tracking technique on Sentinel-1 synthetic aperture radar images obtained from 2017 to 2022. Time-series offsets in the range and azimuth directions were determined from the offset tracking pairs with a temporal baseline of 36 days obtained from January to December of each year. Statistical evaluations of the spatiotemporal variations of the time-series offsets effectively eliminated the erroneous offsets in the original offset fields; the remaining offsets were then used to produce two-dimensional annual ice velocities. The direction of the ice flow of the Cook Ice Shelf was almost constant during the period 2017–2022, and the variations in the magnitude of annual ice velocities were investigated. The annual ice velocities of the Cook East Ice Shelf (CEIS) stayed constant and showed a gradual increase from the grounding line to the ice front, except in the western part. Ice velocities of the western part of the CEIS have not changed much at the grounding line during the 6-year period, while in the dynamic shelf ice zone, ice velocities accelerated by up to 22% because of the development of numerous crevasses and fractures. The ice velocities of the Cook West Ice Shelf (CWIS) were about two times higher than those of the CEIS and tended to increase rapidly from the grounding line to the ice front. The annual ice velocities at the grounding line of CWIS increased rapidly from 1330 to 1450 m/a over 6 years, with 70% of this acceleration observed after 2021. This was attributed to a reduction in the ice shelf volume because of the evolution of surface crevasses and rifts, leading to a decrease in the ice shelf’s buttressing potential. In particular, the loss of a portion of the dynamic shelf ice zone due to a series of ice front collapses in February 2022 likely caused the rapid speed-up of the ice shelf. The results of this study indicate that the buttressing potential of the CWIS and the western part of the CEIS has been significantly reduced, which could mean serious instability of the marine ice sheet in this region.

Frazil ice events: Assessing what to expect in the future

Abstract This article addresses the question: What is expected from frazil ice activity in rivers, taking into account the changing climate? It begins with an overview of what frazil ice is and what is required for the occurrence of frazil ice events, namely a supercooled water column. Methodologies to anticipate frazil ice events in the short term are based on air temperature and water discharge, underlining the significance of these two parameters for any predictive methods. Longer-term approaches, calibrated against past events (hindcasting), are used to anticipate frazil ice activity into the future, with indicators such as frazil ice risk, water temperature and frazil volume. Any of these approaches could conceivably be applied to frazil-prone river stretches. To assess climate impact, each location should be treated separately. River ice dynamics can lead to the formation of a hanging dam, a frequent outcome of frazil ice generation in the early winter, causing flow restriction. Flood modeling and forecasting capabilities have been developed and implemented for operational use. More frequent mid-winter breakups are expected to extend the occurrence of frazil ice events into the winter months – the prediction of these will require climate model output to adequately capture month-to-month variability.

Evaluation of Heat Pump Operation in a Single-Family House

Abstract Since decarbonization of the energy sector is among the top priorities in the EU, with ambitious targets of carbon neutrality until 2050, the energy efficiency of the building stock and the use of renewables are those directions, which can bring the most considerable input towards the achievement of these goals. However, it is not always obvious how to efficiently combine different aspects of low energy demand, availability of renewables, part load operation conditions etc. The results of dynamic IDA ICE simulations highlight that the introduction of renewable low-carbon energy sources should be thoroughly coupled with building energy systems and only their full compatibility can give a high efficiency of the entire energy supply system of the building. Analysing simulation data, it was concluded that for low-energy buildings, heat pumps would not always show higher COP values, compared to buildings with higher energy demand. Ground source heat pump (GSHP) will not always be more efficient than air source heat pump (ASHP). If the ground/groundwater temperature is lower than 10 °C, for GSHP it will be almost impossible to over-compete ASHP within the same system. While COP with radiators and underfloor heating differs only by 10 %, both for ASHP and GSHP.

Consistency in the fingerprints of projected sea level change 2015–2100CE

SUMMARY A large ensemble of ice sheet projections to the end of the 21st century have been compiled within community-based initiatives. These ensembles allow for assessment of uncertainties in projections associated with climate forcing and a wide range of parameters governing ice sheet and shelf dynamics, including ice-ocean interactions. Herein, we compute geographically variable sea level ‘fingerprints’ associated with ∼320 simulations of polar ice sheet projections included in the Ice Sheet Model Intercomparison Project for CMIP6 and ∼180 projections of glacier mass changes from the Glacier Model Intercomparison Project. We find a strong correlation (coefficient > 0.97) between all fingerprints of Greenland Ice Sheet projections when considering a global region outside the near field of the ice sheet. Consistency in the fingerprints for the Antarctic Ice Sheet (AIS) projections is much weaker, though correlation coefficients > 0.80 were found for all projections with global mean sea level (GMSL) greater than 10 cm. The far-field variability in the fingerprints associated with the AIS is due in large part to the sea level change driven by Earth rotation changes. The size and position of the AIS on the south pole makes the rotational component of the sea level fingerprint highly sensitive to the geometry of the ice mass flux, a geometry that becomes more consistent as the GMSL associated with the ice sheet projection increases. Finally, the fingerprints of glacier mass flux show an intermediate level of consistency, with contributions from Antarctic glaciers being the primary driver of decorrelation.

Simulating the Laurentide Ice Sheet of the Last Glacial Maximum

Abstract. In the last decades, great effort has been made to reconstruct the Laurentide Ice Sheet (LIS) during the Last Glacial Maximum (LGM; ca. 21 000 years before present, 21 kyr ago). Uncertainties underlying its modelling have led to notable differences in fundamental features such as its maximum elevation, extent and total volume. As a result, the uncertainty in ice dynamics and thus in ice extent, volume and ice stream stability remains large. We herein use a higher-order three-dimensional ice sheet model to simulate the LIS under LGM boundary conditions for a number of basal friction formulations of varying complexity. Their consequences for the Laurentide ice streams, configuration, extent and volume are explicitly quantified. Total volume and ice extent generally reach a constant equilibrium value that falls close to prior LIS reconstructions. Simulations exhibit high sensitivity to the dependency of the basal shear stress on the sliding velocity. In particular, a regularised Coulomb friction formulation appears to be the best choice in terms of ice volume and ice stream realism. Pronounced differences are found when the basal friction stress is thermomechanically coupled: the base remains colder, and the LIS volume is lower than in the purely mechanical friction scenario counterpart. Thermomechanical coupling is fundamental for producing rapid ice streaming, yet it leads to a similar ice distribution overall.

An experimental study of rime ice accretion on bundled conductors

The dynamic ice accretion process on bundled electrical power conductors and icing-induced effects on the aerodynamic forces acting on the conductors were examined experimentally. An Icing Research Tunnel was used to generate a typical atmospheric rime icing condition experienced by power transmission cables. Two pieces of Aluminum-Conductor-Steel-Reinforced (ACSR) cables were mounted in the icing wind tunnel with different spacing and angular displacements with respect to the incoming airflow. While ice accretion was found to take place primarily on the frontal surfaces of the bundled conductors, the accreted ice layers resembled well the twisted outer strands of the ACSR conductors. Averaged outer profiles of the iced conductors were elliptical with the semi-minor axis being the radius of the conductors and the semi-major axis passing through the foremost frontal points of the accreted ice layers. The icing process on the windward conductor was almost unaffected by the existence of the leeward conductor. The ice accretion and resultant aerodynamic force acting on the leeward conductor were found to vary significantly, depending on the spacing and angular displacements between the bundled conductors. The windward conductor was found to induce a “shadowing effect” by intercepting airborne supercooled water droplets to prevent them from impinging onto the leeward conductor, resulting in much less ice accretion and smaller aerodynamic drag acting on the leeward conductor. The aerodynamic force measurements were correlated with the acquired ice accretion images and wake flow field around the bundled conductors to elucidate the underlying physics.

Ensemble Modeling of Ice Sheet Dynamics in the Last Glacial Cycle

Ensemble simulations (taking into account to the uncertainty of paleoclimate reconstructions) with a model of ice sheet dynamics for the last glacial cycle (128 kyr) are carried out. In general, the model realistically reproduces the spatial structure of the ice sheets and the heights of their domes in the Northern Hemisphere, as well as the associated changes in the ocean level. Perturbations with a sufficiently large amplitude of paleoclimate data in the model show significant differences in the results of modeling the ice sheets of the Northern Hemisphere from the data obtained for the initial paleoreconstruction, including for the Last Glacial Maximum and the time interval of 58‒51 ka (the initial part of MIS3). According to the simulation results, the uncertainty of global reconstructions for the Last Glacial Maximum is 2°С, which is consistent with the existing estimates.

Reservoir Ice Conditions from Multi-Sensor Remote Sensing and ERA5-Land: The Manicouagan Hydroelectric Reservoir Case Study

Reservoir ice can have an important impact on the watershed scale and influence hydraulic operations. On the other hand, hydropower generation can also impact the ice regime. In this study, multi-source satellite and ERA5-land data are used to evaluate ice conditions. Specifically, ice-controlling variables (temperature, water levels), ice regime (cover/deformation, thickness) and their interrelations are assessed for a 5-year period from 2017 to 2021. The methodology is applied to the Manicouagan reservoir, one of the largest hydropower reservoirs in Quebec, Canada. The satellite-based land surface temperatures (LSTs) suggest that winter 2021 was the hottest one. Overall, MODIS and Landsat LSTs agree with the ERA5-land temperatures. Ice backscatter from Sentinel-1 indicates that, in general, the reservoir is completely covered by ice from January to March. A correlation of 0.6 and 0.8 is observed between C- and Ku-band Synthetic Aperture Radar (SAR) signal and ice thickness, respectively. Important ice changes inferred from Differential Interferometric SAR (D-InSAR) occur approximately at the position where the largest ERA5-land ice thickness differences are observed. Winter water levels are also evaluated using satellite altimetric data to verify their influence on the ice dynamics. They show a decreasing tendency as the winter advances.