• European Marine Science
• 2019-2023close
• Other research products close
• European Commission close
• EC|FP7 close
• English close
• The Cryosphere (TC) close

The role of icebergs in narrow fjords hosting marine terminating glaciers in Greenland is poorly understood, even though icebergs provide a substantial freshwater flux that can exceed the subglacial discharge. Iceberg melt is distributed at depth, contributing to fjord stratification, thus impacting melt and dynamics of the glacier front. We model the high-silled Ilulissat Icefjord in Western Greenland with the MITgcm ocean model, using the IceBerg package to study the effect of icebergs on fjord properties, and compare our results with available observations from 2014. We find the subglacial discharge plume to be the primary driver of the seasonality of circulation, glacier melt and iceberg melt. Icebergs are necessary to include to correctly understand the properties of Ilulissat Icefjord, since they modify the fjord in three main ways: First, icebergs cool and freshen the water column within their vertical extent; Second, icebergs depress the neutral buoyancy depth of the plume and the outflow route of glacially modified water; Third, icebergs modify the deep basin, below their vertical extent, due to both increased entrainment of glacially modified water into the fjord, and iceberg modification of the incoming ambient water. Furthermore, the depressed neutral buoyancy depth of the plume limits melt to the deep section of the front of Sermeq Kujalleq (Jakobshavn Isbræ) even during peak summer, and thus promotes undercutting. We postulate that the impact of icebergs on the neutral buoyancy depth of the plume is a key mechanism connecting iceberg melange and glacier calving, independent of mechanical support.

The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting heat flux anomalies with various locations and intensities. In our model experiment, a minimum heat flux value of 970 mW m−2 located close to the East Greenland Ice-core Project (EGRIP) is required locally to reproduce the observed NEGIS velocities, giving basal melt rates consistent with previous estimates. The value cannot be attributed to geothermal heat flux alone and we suggest hydrothermal circulation as a potential explanation for the high local heat flux. By including high heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.

The thermokarst lakes of permafrost regions play a major role in the global carbon cycle. These lakes are sources of methane to the atmosphere but the methane flux is restricted by an ice cover for most of the year. We provide insights into the methane pathways in the winter ice cover on three different water bodies in a continuous permafrost region in Siberia. The first is a bay underlain by submarine permafrost (Tiksi Bay, TB), the second a shallow thermokarst lagoon (Polar Fox, PF) and the third a land-locked, freshwater thermokarst lake (Goltsovoye Lake, GL). In total, 11 ice cores were analyzed as records of the freezing process and methane pathways during the winter season. In TB, the hydrochemical parameters indicate an open system freezing. In contrast, PF was classified as a semi-closed system, where ice growth eventually cuts off exchange between the lagoon and the ocean. The GL is a closed system without connections to other water bodies. Ice on all water bodies was mostly methane-supersaturated with respect to the atmospheric equilibrium concentration, except of three cores from the lake. Generally, the TB ice had low methane concentrations (3.48–8.44 nM) compared to maximum concentrations of the PF ice (2.59–539 nM) and widely varying concentrations in the GL ice (0.02–14817 nM). Stable delta13CCH4 isotope signatures indicate that methane above the ice-water interface was oxidized to concentrations close to or below the calculated atmospheric equilibrium concentration in the ice of PF. We conclude that methane oxidation in ice may decrease methane concentrations during winter. Therefore, understanding seasonal effects to methane pathways in Arctic saline influenced or freshwater systems is critical to anticipate permafrost carbon feedbacks in course of global warming.

The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area, and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet, but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting hotspots with various locations and intensities. We find that a minimum geothermal heat flux value of 970 mW/m2 located close to EastGRIP is required locally to reproduce the observed NEGIS velocities, consistent with previous estimates. By including high geothermal heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.

The growth of frazil or granular ice is an important mode of ice formation in the cryosphere. Recent advances have improved our understanding of the microphysical processes that control the rate of ice-crystal growth when water is cooled beneath its freezing temperature. These advances suggest that crystals grow much faster than previously thought. In this paper, we consider models of a population of ice crystals with different sizes to provide insight into the treatment of frazil ice in large-scale models. We consider the role of crystal growth alongside the other physical processes that determine the dynamics of frazil ice. We apply our model to a simple mixed layer (such as at the surface of the ocean) and to a buoyant plume under a floating ice shelf. We provide numerical calculations and scaling arguments to predict the occurrence of frazil-ice explosions, which we show are controlled by crystal growth, nucleation, and gravitational removal. Faster crystal growth, higher secondary nucleation, and slower gravitational removal make frazil-ice explosions more likely. We identify steady-state crystal size distributions, which are largely insensitive to crystal growth rate but are affected by the relative importance of secondary nucleation to gravitational removal. Finally, we show that the fate of plumes underneath ice shelves is dramatically affected by frazil-ice dynamics. Differences in the parameterization of crystal growth and nucleation give rise to radically different predictions of basal accretion and plume dynamics, and can even impact whether a plume reaches the end of the ice shelf or intrudes at depth.

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.