scholarly journals Sensitivity of the Lambert-Amery glacial system to geothermal heat flux

2016 ◽  
Vol 57 (73) ◽  
pp. 56-68 ◽  
Author(s):  
M. L. Pittard ◽  
J. L. Roberts ◽  
B. K. Galton-Fenzi ◽  
C. S. Watson

ABSTRACTGeothermal heat flux (GHF) is one of the key thermal boundary conditions for ice-sheet models. We assess the sensitivity of the Lambert-Amery glacial system in East Antarctica to four different GHF datasets using a regional ice-sheet model. A control solution of the regional model is initialised by minimising the misfit to observations through an optimisation process. The Lambert-Amery glacial system simulation contains temperate ice up to 150 m thick and has an average basal melt of 1.3 mm a−1, with maximum basal melting of 504 mm a−1. The simulations which use a relatively high GHF compared to the control solution increase the volume and area of temperate ice, which causes higher surface velocities at higher elevations, which leads to the advance of the grounding line. The grounding line advance leads to changes in the local flow configuration, which dominates the changes within the glacial system. To investigate the difference in spatial patterns within the geothermal datasets, they were scaled to have the same median value. These scaled GHF simulations showed that the ice flow was most sensitive to the spatial variation in the underlying GHF near the ice divides and on the edges of the ice streams.

2021 ◽  
Author(s):  
Ole Zeising ◽  
Angelika Humbert

Abstract. The accelerated ice flow of ice streams that reach far into the interior of the ice sheet, is associate with lubrication of the ice sheet base by basal melt water. However, the amount of basal melting under the large ice streams – such as the Northeast Greenland Ice Stream (NEGIS) – are largely unknown. In-situ measurements of basal melt rates are important from various perspectives as they indicate the heat budget, the hydrological regime and the role of sliding in glacier motion. The few previous estimates of basal melt rates in the NEGIS region were 0.1 m a−1 and more, based on radiostratigraphy methods. These finding raised the question of the heat source, since even an increased geothermal heat flux could not deliver the necessary amount of heat. Here, we present basal melt rates at the recent deep drill site EastGRIP, located in the center of NEGIS. Within two subsequent years, we found basal melt rates of (0.16–0.22) ± 0.01 m a−1, that are based on analysis of repeated phase-sensitive radar measurements. In order to quantify the contribution of processes that cause a heat flux into the ice, we carried out an assessment of the energy sources and found the subglacial water system to play a key role in facilitating such high melt rates.


2011 ◽  
Vol 52 (59) ◽  
pp. 43-50 ◽  
Author(s):  
Douglas J. Brinkerhoff ◽  
Toby W. Meierbachtol ◽  
Jesse V. Johnson ◽  
Joel T. Harper

AbstractA full-stress, thermomechanically coupled, numerical model is used to explore the interaction between basal thermal conditions and motion of a terrestrially terminating section of the west Greenland ice sheet. The model domain is a two-dimensional flowline profile extending from the ice divide to the margin. We use data-assimilation techniques based on the adjoint model in order to optimize the basal traction field, minimizing the difference between modeled and observed surface velocities. We monitor the sensitivity of the frozen/melted boundary (FMB) to changes in prescribed geothermal heat flux and sliding speed by applying perturbations to each of these parameters. The FMB shows sensitivity to the prescribed geothermal heat flux below an upper threshold where a maximum portion of the bed is already melted. The position of the FMB is insensitive to perturbations applied to the basal traction field. This insensitivity is due to the short distances over which longitudinal stresses act in an ice sheet.


1996 ◽  
Vol 42 (140) ◽  
pp. 37-42 ◽  
Author(s):  
Jan T. Heine ◽  
David F. Mctigue

AbstractA finite-difference numerical model is used to simulate the temperature profile at the center of an ice sheet throughout the course of a glaciation. The ice sheet is gradually built to a thickness of 3000 m over about 10 000 years, starting on permafrost. A geothermal heat flux is applied at large depth. For an initial surface temperature of –12.5°C, our model shows that basal melting occurs 72000 years after the onset of the glaciation. The important parameters determining the basal temperatures are the initial temperature of the ice and substrate, the rate of downward advection of cold ice and, to a lesser extent, the thickness of the ice sheet. The growth history of the ice sheet does not significantly influence the time at which basal melting occurs. Our results show the possibility that the central parts of the continental ice sheets were cold-based for a significant part of their existence. Heating due to the geothermal heat flux cannot account for basal melting during most or all of a glacial cycle. These results may help to explain the existence of preserved land forms under the ice sheets.


2015 ◽  
Vol 1 (6) ◽  
pp. e1500093 ◽  
Author(s):  
Andrew T. Fisher ◽  
Kenneth D. Mankoff ◽  
Slawek M. Tulaczyk ◽  
Scott W. Tyler ◽  
Neil Foley ◽  
...  

The geothermal heat flux is a critical thermal boundary condition that influences the melting, flow, and mass balance of ice sheets, but measurements of this parameter are difficult to make in ice-covered regions. We report the first direct measurement of geothermal heat flux into the base of the West Antarctic Ice Sheet (WAIS), below Subglacial Lake Whillans, determined from the thermal gradient and the thermal conductivity of sediment under the lake. The heat flux at this site is 285 ± 80 mW/m2, significantly higher than the continental and regional averages estimated for this site using regional geophysical and glaciological models. Independent temperature measurements in the ice indicate an upward heat flux through the WAIS of 105 ± 13 mW/m2. The difference between these heat flux values could contribute to basal melting and/or be advected from Subglacial Lake Whillans by flowing water. The high geothermal heat flux may help to explain why ice streams and subglacial lakes are so abundant and dynamic in this region.


1996 ◽  
Vol 42 (140) ◽  
pp. 37-42 ◽  
Author(s):  
Jan T. Heine ◽  
David F. Mctigue

AbstractA finite-difference numerical model is used to simulate the temperature profile at the center of an ice sheet throughout the course of a glaciation. The ice sheet is gradually built to a thickness of 3000 m over about 10 000 years, starting on permafrost. A geothermal heat flux is applied at large depth. For an initial surface temperature of –12.5°C, our model shows that basal melting occurs 72000 years after the onset of the glaciation. The important parameters determining the basal temperatures are the initial temperature of the ice and substrate, the rate of downward advection of cold ice and, to a lesser extent, the thickness of the ice sheet. The growth history of the ice sheet does not significantly influence the time at which basal melting occurs. Our results show the possibility that the central parts of the continental ice sheets were cold-based for a significant part of their existence. Heating due to the geothermal heat flux cannot account for basal melting during most or all of a glacial cycle. These results may help to explain the existence of preserved land forms under the ice sheets.


2019 ◽  
Vol 13 (3) ◽  
pp. 1043-1049 ◽  
Author(s):  
Tyler Pelle ◽  
Mathieu Morlighem ◽  
Johannes H. Bondzio

Abstract. Basal melting at the bottom of Antarctic ice shelves is a major control on glacier dynamics, as it modulates the amount of buttressing that floating ice shelves exert onto the ice streams feeding them. Three-dimensional ocean circulation numerical models provide reliable estimates of basal melt rates but remain too computationally expensive for century-scale projections. Ice sheet modelers therefore routinely rely on simplified parameterizations based on either ice shelf depth or more sophisticated box models. However, existing parameterizations do not accurately resolve the complex spatial patterns of sub-shelf melt rates that have been observed over Antarctica's ice shelves, especially in the vicinity of the grounding line, where basal melting is one of the primary drivers of grounding line migration. In this study, we couple the Potsdam Ice-shelf Cavity mOdel (PICO, Reese et al., 2018) to a buoyant plume melt rate parameterization (Lazeroms et al., 2018) to create PICOP, a novel basal melt rate parameterization that is easy to implement in transient ice sheet numerical models and produces a melt rate field that is in excellent agreement with the spatial distribution and magnitude of observations for several ocean basins. We test PICOP on the Amundsen Sea sector of West Antarctica, Totten, and Moscow University ice shelves in East Antarctica and the Filchner-Ronne Ice Shelf and compare the results to PICO. We find that PICOP is able to reproduce inferred high melt rates beneath Pine Island, Thwaites, and Totten glaciers (on the order of 100 m yr−1) and removes the “banding” pattern observed in melt rates produced by PICO over the Filchner-Ronne Ice Shelf. PICOP resolves many of the issues contemporary basal melt rate parameterizations face and is therefore a valuable tool for those looking to make future projections of Antarctic glaciers.


2020 ◽  
Vol 14 (7) ◽  
pp. 2283-2301 ◽  
Author(s):  
Stephen L. Cornford ◽  
Helene Seroussi ◽  
Xylar S. Asay-Davis ◽  
G. Hilmar Gudmundsson ◽  
Rob Arthern ◽  
...  

Abstract. We present the result of the third Marine Ice Sheet Model Intercomparison Project, MISMIP+. MISMIP+ is intended to be a benchmark for ice-flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient to model buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments first tests that models are able to maintain a steady state with the grounding line located on a retrograde slope due to buttressing and then explore scenarios where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. The majority of participating models passed the first test and then produced similar responses to the loss of buttressing. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress but also the differences between numerical methods – taking a secondary role.


2020 ◽  
Vol 66 (256) ◽  
pp. 303-312
Author(s):  
Silje Smith-Johnsen ◽  
Basile de Fleurian ◽  
Kerim H. Nisancioglu

AbstractThe spatial distribution of geothermal heat flux (GHF) under ice sheets is largely unknown. Nonetheless, it is an important boundary condition in ice-sheet models, and suggested to control part of the complex surface velocity patterns observed in some regions. Here we investigate the effect of including subglacial hydrology when modelling ice streams with elevated GHF. We use an idealised ice stream geometry and a thermomechanical ice flow model coupled to subglacial hydrology in the Ice Sheet System Model (ISSM). Our results show that the dynamic response of the ice stream to elevated GHF is greatly enhanced when including the interactive subglacial hydrology. On the other hand, the impact of GHF on ice temperature is reduced when subglacial hydrology is included. In conclusion, the sensitivity of ice stream dynamics to GHF is likely to be underestimated in studies neglecting subglacial hydrology.


2001 ◽  
Vol 47 (157) ◽  
pp. 271-282 ◽  
Author(s):  
Richard C.A. Hindmarsh ◽  
E. Le Meur

AbstractMarine ice sheets with mechanics described by the shallow-ice approximation by definition do not couple mechanically with the shelf. Such ice sheets are known to have neutral equilibria. We consider the implications of this for their dynamics and in particular for mechanisms which promote marine ice-sheet retreat. The removal of ice-shelf buttressing leading to enhanced flow in grounded ice is discounted as a significant influence on mechanical grounds. Sea-level rise leading to reduced effective pressures under ice streams is shown to be a feasible mechanism for producing postglacial West Antarctic ice-sheet retreat but is inconsistent with borehole evidence. Warming thins the ice sheet by reducing the average viscosity but does not lead to grounding-line retreat. Internal oscillations either specified or generated via a MacAyeal–Payne thermal mechanism promote migration. This is a noise-induced drift phenomenon stemming from the neutral equilibrium property of marine ice sheets. This migration occurs at quite slow rates, but these are sufficiently large to have possibly played a role in the dynamics of the West Antarctic ice sheet after the glacial maximum. Numerical experiments suggest that it is generally true that while significant changes in thickness can be caused by spatially uniform changes, spatial variability coupled with dynamical variability is needed to cause margin movement.


2021 ◽  
Author(s):  
Parviz Ajourlou ◽  
François PH Lapointe ◽  
Glenn A Milne ◽  
Yasmina Martos

<p>Geothermal heat flux (GHF) is known to be an important control on the basal thermal state of an ice sheet which, in turn, is a key factor in governing how the ice sheet will evolve in response to a given climate forcing. In recent years, several studies have estimated GHF beneath the Greenland ice sheet using different approaches (e.g. Rezvanbehbahani et al., Geophysical Research Letters, 2017; Martos et al., Geophysical Research Letters, 2018; Greve, Polar Data Journal, 2019). Comparing these different estimates indicates poor agreement and thus large uncertainty in our knowledge of this important boundary condition for modelling the ice sheet. The primary aim of this study is to quantify the influence of this uncertainty on modelling the past evolution of the ice sheet with a focus on the most recent deglaciation. We build on past work that considered three GHF models (Rogozhina et al., 2011) by considering over 100 different realizations of this input field. We use the uncertainty estimates from Martos et al. (Geophysical Research Letters, 2018) to generate GHF realisations via a statistical sampling procedure. A sensitivity analysis using these realisations and the Parallel Ice Sheet Model (PISM, Bueler and Brown, Journal of Geophysical Research, 2009) indicates that uncertainty in GHF has a dramatic impact on both the volume and spatial distribution of ice since the last glacial maximum, indicating that more precise constraints on this boundary condition are required to improve our understanding of past ice sheet evolution and, consequently, reduce uncertainty in future projections.</p>


Sign in / Sign up

Export Citation Format

Share Document