Large-eddy simulations of the ice-shelf-ocean boundary layer near the ice front of Nansen Ice Shelf, Antarctica

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.