Search

Arctic lakes are exposed to warming during increasingly longer ice-free periods and, if located in glaciated areas, to increased inflow of meltwater and sediments. However, direct monitoring of how such lakes respond to changing environmental conditions is challenging not only because of their remoteness but also because of the scarcity of present and previously observed lake states. At the glacier-proximal Lake Tarfala in the Kebnekaise Mountains, northern Sweden, temperatures throughout the water column at its deepest part (50 m) were acquired between 2016 and 2019. This three-year record shows that Lake Tarfala is dimictic and is overturning during spring and fall, respectively. Timing, duration, and intensity of mixing processes, as well as of summer and winter stratification, vary between years. Glacial meltwater may play an important role regarding not only mixing processes but also cooling of the lake. Attribution of external environmental factors to (changes in) lake mixing processes and thermal states remains challenging owing to for example, timing of ice-on and ice-off but also reflection and absorption of light, both known to play a decisive role for lake mixing processes, are not (yet) monitored in situ at Lake Tarfala. publishedVersion

Arctic lakes are exposed to warming during increasingly longer ice-free periods and, if located in glaciated areas, to increased inflow of meltwater and sediments. However, direct monitoring of how such lakes respond to changing environmental conditions is challenging not only because of their remoteness but also because of the scarcity of present and previously observed lake states. At the glacier-proximal Lake Tarfala in the Kebnekaise Mountains, northern Sweden, temperatures throughout the water column at its deepest part (50 m) were acquired between 2016 and 2019. This three-year record shows that Lake Tarfala is dimictic and is overturning during spring and fall, respectively. Timing, duration, and intensity of mixing processes, as well as of summer and winter stratification, vary between years. Glacial meltwater may play an important role regarding not only mixing processes but also cooling of the lake. Attribution of external environmental factors to (changes in) lake mixing processes and thermal states remains challenging owing to for example, timing of ice-on and ice-off but also reflection and absorption of light, both known to play a decisive role for lake mixing processes, are not (yet) monitored in situ at Lake Tarfala.

Predictions of future mass loss from ice sheets are afflicted with uncertainty, caused, among others, by insufficient understanding of spatiotemporally variable processes at the inaccessible base of ice sheets for which few direct observations exist and of which basal friction is a prime example. Here, we present a general numerical framework for studying the relationship between bed and surface properties of ice sheets and glaciers. Specifically, we use an inverse modeling approach and the associated time-dependent adjoint equations, derived in the framework of a full Stokes model and a shallow-shelf/shelfy-stream approximation model, respectively, to determine the sensitivity of grounded ice sheet surface velocities and elevation to time-dependent perturbations in basal friction and basal topography. Analytical and numerical examples are presented showing the importance of including the time-dependent kinematic free surface equation for the elevation and its adjoint, in particular for observations of the elevation. A closed form of the analytical solutions to the adjoint equations is given for a two-dimensional vertical ice in steady state under the shallow-shelf approximation. There is a delay in time between a seasonal perturbation at the ice base and the observation of the change in elevation. A perturbation at the base in the topography has a direct effect in space at the surface above the perturbation, and a perturbation in the friction is propagated directly to the surface in time.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve due to the nonlinearity introduced by the rheology. Therefore, approximations to the FS equations are commonly used, especially when modeling a marine ice sheet (ice sheet, ice shelf, and/or ice stream) for 103 years or longer. The shallow ice approximation (SIA) and shallow shelf approximation (SSA) are commonly used but are accurate only for certain parts of an ice sheet. Here, we report a novel way of iteratively coupling FS and SSA that has been implemented in Elmer/Ice and applied to conceptual marine ice sheets. The FS-SSA coupling appears to be very accurate; the relative error in velocity compared to FS is below 0.5 % for diagnostic runs and below 5 % for prognostic runs. Results for grounding line dynamics obtained with the FS-SSA coupling are similar to those obtained from an FS model in an experiment with a periodical temperature forcing over 3000 years that induces grounding line advance and retreat. The rapid convergence of the FS-SSA coupling shows a large potential for reducing computation time, such that modeling a marine ice sheet for thousands of years should become feasible in the near future. Despite inefficient matrix assembly in the current implementation, computation time is reduced by 32 %, when the coupling is applied to a 3-D ice shelf..

GRANTISM (GReenland and ANTarctic Ice Sheet Model) is an educational Excel™ model introduced by Pattyn (2006). Here, GRANTISM is amended to simulate the Svalbard-Barents-Sea Ice Sheet during the Last Glacial Maximum, an analogue for the contemporary West Antarctic Ice Sheet. A new name, “GRANTSISM,” is suggested; the added S represents Svalbard. GRANTSISM introduces students of bachelor's or master's programs in Earth sciences (first or second cycle program in the Bologna system for higher education), but with little or no background in numerical modeling, to basic ice sheet modeling. GRANTSISM provides hands-on learning experiences related to ice sheet dynamics in response to climate forcing, and fosters understanding of processes and feedbacks. GRANTSISM was successfully used in noncompulsory courses in which students have been able to reproduce paleo-ice sheet evolution scenarios discussed here as examples. Students progressed further by designing, developing, and analyzing their own modeling scenarios. Here, we describe GRANTSISM and report on how learning activities with GRANTSISM were assessed by students who had no prior experience in ice sheet modeling. The response rate for a noncompulsory survey of the learning activity was less than 40%. A subsequent control experiment with a compulsory survey, however, showed the same patterns of answers, so the student response is considered representative. First, GRANTSISM is concluded to be a highly attractive tool to introduce learners with an interest in ice sheet behavior to ice sheet modeling. Second, it triggers an interest for more in-depth learning experiences related to numerical ice sheet modeling.

Geophysical mapping and coring of the central Arctic Ocean seafloor provide evidence for repeated occurrences of ice sheet/ice shelf complexes during previous glacial periods. Several ridges and bathymetric highs shallower than present water depths of ~1000?m show signs of erosion from deep drafting (armadas of) icebergs which originated from thick outlet glaciers and ice shelves. Mapped glacigenic landforms and dates of cored sediments suggest that the largest ice shelf complex was confined to the Amerasian sector of the Arctic Ocean during Marine Isotope Stage (MIS) 6. However the spatial extent of ice shelves can not be well reconstructed from occasional groundings on bathymetric highs. Therefore we apply a statistical approach to provide independent support for an extensive MIS 6 ice shelf complex which previously was inferred only from interpretation of geophysical and geological data. Specifically we assess whether this ice shelf complex comprises a likely source of the deep draft icebergs responsible for the mapped scour marks. The statistical modeling is based on exploiting relations between contemporary Antarctic ice shelves and their local physical environments and the assumption that Arctic Ocean MIS?6 ice shelves scale similarly. Analyzing ice thickness data along the calving front of contemporary ice shelves a peak over threshold method is applied to determine sources of deep drafting icebergs in the Arctic Ocean MIS?6 ice shelf complex. This approach is novel to modeling Arctic paleoglacial configurations. Predicted extreme calving front drafts match observed deep draft iceberg scours if the ice shelf complex is sufficiently large.

Ice sheets extending over parts of the East Siberian continental shelf have been proposed for the last glacial period and during the larger Pleistocene glaciations. The sparse data available over this sector of the Arctic Ocean have left the timing, extent and even existence of these ice sheets largely unresolved. Here we present new geophysical mapping and sediment coring data from the East Siberian shelf and slope collected during the 2014 SWERUS-C3 expedition (SWERUS-C3: Swedish-Russian-US Arctic Ocean Investigation of Climate-Cryosphere-Carbon Interactions). The multibeam bathymetry and chirp sub-bottom profiles reveal a set of glacial landforms that include grounding zone formations along the outer continental shelf, seaward of which lies a > 65-m thick sequence of glacio-genic debris flows. The glacial landforms are interpreted to lie at the seaward end of a glacial trough-the first to be reported on the East Siberian margin, here referred to as the De Long Trough because of its location due north of the De Long Islands. Stratigraphy and dating of sediment cores show that a drape of acoustically laminated sediments covering the glacial deposits is older than 50 kyr BP. This provides direct evidence for extensive glacial activity on the Siberian shelf that predates the Last Glacial Maximum and most likely occurred during the Saalian (Marine Isotope Stage (MIS) 6).

I. Microphysical properties, deformation, texture and grain growth, Polar ice masses, including grounded (inland) ice, ice shelves, icebergs and sea ice, constitute an essential part of the global climate system. They interact with the hydrosphere and the atmosphere, and, if persisting through millennia, serve as archives of the past climate of Earth. Currently. significant efforts in climate reconstructions are being undertaken, serving as a basis for climate predictions into the future. It has been argued in Part I of this work that more attention should be given to interactions between structures on distinct scales (multiscale structural interactions), which seem to play a decisive role in the structure-form-environment interplay (SFEI) in grounded polar ice. As ice shelves play an equally important role in the global climate, the present article is devoted to the description of the multiscale structure of ice shelves.

The Tibetan Plateau is a topographic feature of extraordinary dimension and has an important impact on regional and global climate. Yet, the glacial history of the Tibetan Plateau is less constrained than the history of some other formerly glaciated regions, especially in the Northern Hemisphere (e.g. Laurentide Ice Sheet, Fennoscandian Ice Sheet). Nevertheless, field evidence for extensive valley glaciation indicates that ice sheet glaciation on the Tibetan Plateau did not evolve during the Last Glacial maximum (LGM). This is an important and robust result that has not been widely investigated using numerical ice sheet models, despite potentially important climate ramifications. Perhaps this is because reconstructions of the LGM glacial configurations of the Tibetan Plateau in the framework of numerical simulations covering an entire glacial cycle exhibit a pronounced variability then entire range of which is not supported by field evidence. Using the 3d thermomechanical ice sheet model SICOPOLIS, we simulated the evolution of Tibetan Plateau ice configurations during the last 125.000 years. Temperature and precipitation data driving the simulations have been applied in the form of a large ensemble of glacial/interglacial climate scenarios. It is observed that variations in ice sheet configuration resulting from the prescription of different present-day precipitation- and temperature data sets, on the one hand, and different paleoclimates as obtained from reconstructions based on different GCM-model outputs, on the other hand, include as extreme end members an entirely ice free Tibetan Plateau during the last glacial cycle as well as a plateau-scale Tibetan Ice sheet during the LGM. Comparison of such numerical results with available field data indicates that further refinements in the numerical simulations are required, and that these must include atmosphere-ice sheet feedback mechanisms. However, because mapped and simulated glacial extents are represented at different spatial scales, this task is not straightforward.

Arctic lakes are exposed to warming during increasingly longer ice-free periods and, if located in glaciated areas, to increased inflow of meltwater and sediments. However, direct monitoring of how such lakes respond to changing environmental conditions is challenging not only because of their remoteness but also because of the scarcity of present and previously observed lake states. At the glacier-proximal Lake Tarfala in the Kebnekaise Mountains, northern Sweden, temperatures throughout the water column at its deepest part (50 m) were acquired between 2016 and 2019. This three-year record shows that Lake Tarfala is dimictic and is overturning during spring and fall, respectively. Timing, duration, and intensity of mixing processes, as well as of summer and winter stratification, vary between years. Glacial meltwater may play an important role regarding not only mixing processes but also cooling of the lake. Attribution of external environmental factors to (changes in) lake mixing processes and thermal states remains challenging owing to for example, timing of ice-on and ice-off but also reflection and absorption of light, both known to play a decisive role for lake mixing processes, are not (yet) monitored in situ at Lake Tarfala.

Arctic lakes are exposed to warming during increasingly longer ice-free periods and, if located in glaciated areas, to increased inflow of meltwater and sediments. However, direct monitoring of how such lakes respond to changing environmental conditions is challenging not only because of their remoteness but also because of the scarcity of present and previously observed lake states. At the glacier-proximal Lake Tarfala in the Kebnekaise Mountains, northern Sweden, temperatures throughout the water column at its deepest part (50 m) were acquired between 2016 and 2019. This three-year record shows that Lake Tarfala is dimictic and is overturning during spring and fall, respectively. Timing, duration, and intensity of mixing processes, as well as of summer and winter stratification, vary between years. Glacial meltwater may play an important role regarding not only mixing processes but also cooling of the lake. Attribution of external environmental factors to (changes in) lake mixing processes and thermal states remains challenging owing to for example, timing of ice-on and ice-off but also reflection and absorption of light, both known to play a decisive role for lake mixing processes, are not (yet) monitored in situ at Lake Tarfala.

This work is affiliated with the CRIOS project (Calving Rates and Impact On Sea Level), which was supported by the Conoco Phillips-Lundin Northern Area Program. PH is funded by a NERC PhD studentship (reference number 1396698). We present a highly detailed study of calving dynamics at Tunabreen, a tidewater glacier in Svalbard. A time-lapse camera was trained on the terminus and programmed to capture images every 3 seconds over a 28-hour period in August 2015, producing a highly detailed record of 34 117 images from which 358 individual calving events were distinguished. Calving activity is characterised by frequent events (12.8 events h-1) that are small relative to the spectrum of calving events observed, demonstrating the prevalence of small-scale calving mechanisms. Five calving styles were observed, with a high proportion of calving events (82%) originating at, or above, the waterline. The tidal cycle plays a key role in the timing of calving events, with 68% occurring on the falling limb of the tide. Calving activity is concentrated where meltwater plumes surface at the glacier front, and a ∼ 5 m undercut at the base of the glacier suggests that meltwater plumes encourage melt-under-cutting. We conclude that frontal ablation at Tunabreen may be paced by submarine melt rates, as suggested from similar observations at glaciers in Svalbard and Alaska. Using submarine melt rate to calculate frontal ablation would greatly simplify estimations of tidewater glacier losses in prognostic models. Publisher PDF

Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve due to the nonlinearity introduced by the rheology. Therefore, approximations to the FS equations are commonly used, especially when modeling a marine ice sheet (ice sheet, ice shelf, and/or ice stream) for 103 years or longer. The shallow ice approximation (SIA) and shallow shelf approximation (SSA) are commonly used but are accurate only for certain parts of an ice sheet. Here, we report a novel way of iteratively coupling FS and SSA that has been implemented in Elmer/Ice and applied to conceptual marine ice sheets. The FS–SSA coupling appears to be very accurate; the relative error in velocity compared to FS is below 0.5% for diagnostic runs and below 5% for prognostic runs. Results for grounding line dynamics obtained with the FS–SSA coupling are similar to those obtained from an FS model in an experiment with a periodical temperature forcing over 3000 years that induces grounding line advance and retreat. The rapid convergence of the FS–SSA coupling shows a large potential for reducing computation time, such that modeling a marine ice sheet for thousands of years should become feasible in the near future. Despite inefficient matrix assembly in the current implementation, computation time is reduced by 32%, when the coupling is applied to a 3-D ice shelf., Geoscientific Model Development, 11 (11), ISSN:1991-9603, ISSN:1991-959X

The Barents Sea Ice Sheet was a marine-based ice sheet, i.e., it rested on the Barents Sea floor during the Last Glacial Maximum (21 ky BP). The Bjørnøyrenna Ice Stream was the largest ice stream draining the Barents Sea Ice Sheet and is regarded as an analogue for contemporary ice streams in West Antarctica. Here, the retreat of the Bjørnøyrenna Ice Stream is simulated by means of two numerical ice sheet models and results assessed against geological data. We investigate the sensitivity of the ice stream to changes in ocean temperature and the impact of grounding-line physics on ice stream retreat. Our results suggest that the role played by sub-shelf melting depends on how the grounding-line physics is represented in the models. When an analytic constraint on the ice flux across the grounding line is applied, the retreat of Bjørnøyrenna Ice Stream is primarily driven by internal ice dynamics rather than by oceanic forcing. This suggests that implementations of grounding-line physics need to be carefully assessed when evaluating and predicting the response of contemporary marine-based ice sheets and individual ice streams to ongoing and future ocean warming.