Search

A realistic initialization of ice flow models is critical for predicting future changes in ice sheet mass balance and their associated contribution to sea level rise. The initial thermal state of an ice sheet is particularly important, as it controls ice viscosity and basal conditions, thereby influencing the overall ice velocity. Englacial and subglacial conditions, however, remain poorly understood due to insufficient direct measurements, which complicate the initialization and validation of thermal models. Here, we investigate the impact of using different geothermal heat flux (GHF) datasets and vertical velocity profiles on the thermal state of the Antarctic ice sheet and compare our modeled temperatures to in situ measurements from 15 boreholes. We find that the temperature profile is more sensitive to vertical velocity than to GHF. The basal temperature of grounded ice and the amount of basal melting are influenced by both selection of GHF and vertical velocity. More importantly, we find that the standard approach, which consists of combining basal sliding speed and incompressibility to derive vertical velocities, provides reasonably good results in fast-flow regions (ice velocity >50 m yr−1) but performs poorly in slow-flow regions (ice velocity m yr−1). Furthermore, the modeled temperature profiles in ice streams, where bed geometry is generated using a mass conservation approach, show better agreement with observed borehole temperatures compared to kriging-based bed geometry.

Ice melting beneath Antarctic ice shelves is caused by heat transfer through the ice-shelf–ocean boundary layer (IOBL). However, our understanding of the fluid dynamics and thermohaline physics of the IOBL flow is poor. In this study, we utilize a large-eddy simulation (LES) model to investigate ocean dynamics and the role of turbulence within the IOBL flow near the ice front. To simulate the varying turbulence intensities, we imposed different theoretical profiles of the velocity. Far-field ocean conditions for the melting at the ice-shelf base and freezing at the sea surface were derived based on in situ observations of temperature and salinity near the ice front of the Nansen Ice Shelf. In terms of overturning features near the ice front, we validated the LES simulation results by comparing them with the in situ observational data. In the comparison of the velocity profiles to shipborne lowered acoustic Doppler current profiler (LADCP) data, the LES-derived strength of the overturning cells is similar to that obtained from the observational data. Moreover, the vertical distribution of the simulated temperature and salinity, which were mainly determined by the positively buoyant meltwater and sea-ice formation, was also comparable to that of the observations. We conclude that the IOBL flow near the ice front and its contribution to the ocean dynamics can be realistically resolved using our proposed method. Based on validated 3D-LES results, we revealed that the main forces of ocean dynamics near the ice front are driven by positively buoyant meltwater, concentrated salinity at the sea surface, and outflowing momentum of the sub-ice-shelf plume. Moreover, in the strong-turbulence case, distinct features such as a higher basal melt rate (0.153 m yr−1), weak upwelling of the positively buoyant ice-shelf water, and a higher sea-ice formation were observed, suggesting a relatively high speed current within the IOBL because of highly turbulent mixing. The findings of this study will contribute toward a deeper understanding of the complex IOBL-flow physics and its impact on the ocean dynamics near the ice front.