scholarly journals Polythermal modelling of steady states of the Antarctic ice sheet in comparison with the real world

1996 ◽  
Vol 23 ◽  
pp. 382-387 ◽  
Author(s):  
I. Hansen ◽  
R. Greve

An approach to simulate the present Antarctic ice sheet with reaped to its thermomechanical behaviour and the resulting features is made with the three-dimensional polythermal ice-sheet model designed by Greve and Hutter. It treats zones of cold and temperate ice as different materials with their own properties and dynamics. This is important because an underlying layer of temperate ice can influence the ice sheet as a whole, e.g. the cold ice may slide upon the less viscous binary ice water mixture.Measurements indicate that the geothermal heat flux below the Antarctic ice sheet appears to be remarkably higher than the standard value of 42 m W m−2 that is usually applied for Precambrian shields in ice-sheet modelling. Since the extent of temperate ice at the base is highly dependent on this heat input from the lithosphere, an adequate choice is crucial for realistic simulations. We shall present a series of steady-state results with varied geothermal heat flux and demonstrate that the real ice-sheet topography can be reproduced fairly well with a value in the range 50–60 m W m−2. Thus, the physical parameters of ice (especially the enhancement factor in Glen’s flow law) as used by Greve (1995) for polythermal Greenland ice-sheet simulations can be adopted without any change. The remaining disagreements may he explained by the neglected influence of the ice shelves, the rather coarse horizontal resolution (100 km), the steady-state assumption and possible shortcomings in the parameterization of the surface mass balance.

1996 ◽  
Vol 23 ◽  
pp. 382-387 ◽  
Author(s):  
I. Hansen ◽  
R. Greve

An approach to simulate the present Antarctic ice sheet with reaped to its thermomechanical behaviour and the resulting features is made with the three-dimensional polythermal ice-sheet model designed by Greve and Hutter. It treats zones of cold and temperate ice as different materials with their own properties and dynamics. This is important because an underlying layer of temperate ice can influence the ice sheet as a whole, e.g. the cold ice may slide upon the less viscous binary ice water mixture. Measurements indicate that the geothermal heat flux below the Antarctic ice sheet appears to be remarkably higher than the standard value of 42 m W m−2 that is usually applied for Precambrian shields in ice-sheet modelling. Since the extent of temperate ice at the base is highly dependent on this heat input from the lithosphere, an adequate choice is crucial for realistic simulations. We shall present a series of steady-state results with varied geothermal heat flux and demonstrate that the real ice-sheet topography can be reproduced fairly well with a value in the range 50–60 m W m−2. Thus, the physical parameters of ice (especially the enhancement factor in Glen’s flow law) as used by Greve (1995) for polythermal Greenland ice-sheet simulations can be adopted without any change. The remaining disagreements may he explained by the neglected influence of the ice shelves, the rather coarse horizontal resolution (100 km), the steady-state assumption and possible shortcomings in the parameterization of the surface mass balance.


1995 ◽  
Vol 21 ◽  
pp. 144-148
Author(s):  
Garth W. Paltridge ◽  
Christopher M. Zweck

A simple steady-state energy and mass-balance model of the Antarctic ice sheet is developed. Basically it is a set of two equations with two unknowns of steady-state height h and potential basal temperature Tb. Tb determines whether, and to what extent, there is liquid water at the base of the ice which in turn affects the values of h and Tb. Simultaneous changes of sea-level temperature and precipitation (changes related to each other as might be expected from global climate models) indicate a maximum in the field of possible steady-state ice volumes which may not be far from the presently observed conditions. The possibility of cyclical variation in ground heat flux associated with convection of water and heat in the continental crust is discussed. The mechanism might be capable of generating cycles of ice-sheet volume with relatively short periods similar to those of Milankovitch forcing.


1989 ◽  
Vol 12 ◽  
pp. 16-22 ◽  
Author(s):  
W.F. Budd ◽  
D. Jenssen

A three-dimensional dynamic, thermodynamic ice-sheet model has been developed to simulate the past, present, and future behaviour of the Antarctic ice sheet. The present ice velocities depend on the deep ice temperatures which in turn depend on the past changes of the ice sheet, including surface temperature, accumulation rate, and ice thickness. The basal temperatures are also strongly dependent on the geothermal heat flux. The model has therefore been used to study the effect on the basal temperatures, of changes to the geothermal heat flux, as well as the past changes of surface temperature and accumulation rate based on results obtained from the Vostok deep ice core. The model is also used to compute the distribution of surface velocity required to balance the present accumulation rate and the dynamics velocity based on the stress, temperature, and flow properties of ice, for the internal deformation, plus a component due to ice sliding. These velocities are compared to observed surface velocities in East Antarctica to assess the state of balance and the performance of the dynamics formulation.


1995 ◽  
Vol 21 ◽  
pp. 144-148 ◽  
Author(s):  
Garth W. Paltridge ◽  
Christopher M. Zweck

A simple steady-state energy and mass-balance model of the Antarctic ice sheet is developed. Basically it is a set of two equations with two unknowns of steady-state heighthand potential basal temperatureTb.Tbdetermines whether, and to what extent, there is liquid water at the base of the ice which in turn affects the values ofhandTb. Simultaneous changes of sea-level temperature and precipitation (changes related to each other as might be expected from global climate models) indicate a maximum in the field of possible steady-state ice volumes which may not be far from the presently observed conditions. The possibility of cyclical variation in ground heat flux associated with convection of water and heat in the continental crust is discussed. The mechanism might be capable of generating cycles of ice-sheet volume with relatively short periods similar to those of Milankovitch forcing.


1989 ◽  
Vol 12 ◽  
pp. 16-22 ◽  
Author(s):  
W.F. Budd ◽  
D. Jenssen

A three-dimensional dynamic, thermodynamic ice-sheet model has been developed to simulate the past, present, and future behaviour of the Antarctic ice sheet. The present ice velocities depend on the deep ice temperatures which in turn depend on the past changes of the ice sheet, including surface temperature, accumulation rate, and ice thickness. The basal temperatures are also strongly dependent on the geothermal heat flux. The model has therefore been used to study the effect on the basal temperatures, of changes to the geothermal heat flux, as well as the past changes of surface temperature and accumulation rate based on results obtained from the Vostok deep ice core. The model is also used to compute the distribution of surface velocity required to balance the present accumulation rate and the dynamics velocity based on the stress, temperature, and flow properties of ice, for the internal deformation, plus a component due to ice sliding. These velocities are compared to observed surface velocities in East Antarctica to assess the state of balance and the performance of the dynamics formulation.


2020 ◽  
Author(s):  
Pavel Talalay ◽  
Yazhou Li ◽  
Laurent Augustin ◽  
Gary Clow ◽  
Jialin Hong ◽  
...  

Abstract. The temperature at the Antarctic ice sheet bed and the temperature gradient in subglacial rocks have been directly measured only a few times, although extensive thermodynamic modelling has been used to estimate geothermal heat flux under ice sheet. During the last five decades, deep ice-core drilling projects at six sites – Byrd, WAIS Divide, Dome C, Kohnen, Dome F, and Vostok – have succeeded in reaching to, or nearly to, the bed in inland locations in Antarctica. When temperature profiles in these boreholes and heat flow model are combined with estimations of vertical velocity, the heat flow at ice sheet base is translated to a geothermal heat flux of 117.8 ± 3.3 mW m−2 at Byrd, 67.3 ± 8.6 mW m−2 at Dome C, 79.0 ± 5.0 mW m−2 at Dome F, and −3.3 ± 5.6 mW m−2 at Vostok, close to predicted values. However, estimations at Kohnen and WAIS Divide gave flux of 161.5 ± 10.2 mW m−2 and 251.3 ± 24.1 mW m−2, respectively, far higher than that predicted by existing heat flow models. The question arises as to whether this high heat flow represents regional values, or if the Kohnen and WAIS Divide boreholes were drilled over local hot spots.


2020 ◽  
Vol 14 (11) ◽  
pp. 4021-4037
Author(s):  
Pavel Talalay ◽  
Yazhou Li ◽  
Laurent Augustin ◽  
Gary D. Clow ◽  
Jialin Hong ◽  
...  

Abstract. The temperature at the Antarctic Ice Sheet bed and the temperature gradient in subglacial rocks have been directly measured only a few times, although extensive thermodynamic modeling has been used to estimate the geothermal heat flux (GHF) under the ice sheet. During the last 5 decades, deep ice-core drilling projects at six sites – Byrd, WAIS Divide, Dome C, Kohnen, Dome F, and Vostok – have succeeded in reaching or nearly reaching the bed at inland locations in Antarctica. When temperature profiles in these boreholes and steady-state heat flow modeling are combined with estimates of vertical velocity, the heat flow at the ice-sheet base is translated to a geothermal heat flux of 57.9 ± 6.4 mW m−2 at Dome C, 78.9 ± 5.0 mW m−2 at Dome F, and 86.9 ± 16.6 mW m−2 at Kohnen, all higher than the predicted values at these sites. This warm base under the East Antarctic Ice Sheet (EAIS) could be caused by radiogenic heat effects or hydrothermal circulation not accounted for by the models. The GHF at the base of the ice sheet at Vostok has a negative value of −3.6 ± 5.3 mW m−2, indicating that water from Lake Vostok is freezing onto the ice-sheet base. Correlation analyses between modeled and measured depth–age scales at the EAIS sites indicate that all of them can be adequately approximated by a steady-state model. Horizontal velocities and their variation over ice-age cycles are much greater for the West Antarctic Ice Sheet than for the interior EAIS sites; a steady-state model cannot precisely describe the temperature distribution here. Even if the correlation factors for the best fitting age–depth curve are only moderate for the West Antarctic sites, the GHF values estimated here of 88.4 ± 7.6 mW m−2 at Byrd and 113.3 ± 16.9 mW m−2 at WAIS Divide can be used as references before more precise estimates are made on the subject.


2021 ◽  
Author(s):  
Maximilian Lowe ◽  
Fausto Ferraccioli ◽  
Duncan Young ◽  
Donald Blankenship ◽  
Egidio Armadillo ◽  
...  

<p>The Wilkes Subglacial Basin in East Antarctica hosts one of the largest marine-based and hence potentially more unstable sectors of the East Antarctic Ice Sheet (EAIS). Predicting the past, present and future behaviour of this key sector of the EAIS requires that we also improve our understanding of the lithospheric cradle on which it flows. This is particularly important in order to quantify geothermal heat flux heterogeneity in the region.   </p><p>The WSB stretches for almost 1600 km from the Southern Ocean towards South Pole. Like many intracratonic basins, it is a long-lived geological feature, which originated and evolved in different tectonic settings. A wide basin formed in the WSB in a distal back arc basin setting, likely in response to a retreating West Antarctic Paleo-Pacific active margin from Permo-Triassic times. Jurassic extension then led to the emplacement of part of a huge flood basalt province that extends from South Africa to Australia.  The region was then affected by relatively minor upper crustal Mesozoic to Cenozoic(?) extension and transtension, producing narrow graben-like features that were glacially overdeepened, and presently steer enhanced glacial flow of the Matusevich, Cook and Ninnis glaciers.</p><p>Here we present the results of our enhanced geophysical imaging and modelling in the WSB region performed within the 4D Antarctica project of ESA, which aims to help quantify the spatial variability in subglacial Antarctic geothermal heat flux (GHF), one of the least well constrained parameters of the entire continent.</p><p>We exploit a combination of airborne radar and aeromagnetic data compilations and crustal and lithosphere thickness estimates from both satellite and airborne gravity and independent passive seismic constraints to develop new geophysical models for the region. To help constrain the starting models, including depth to basement beneath the Permian to Jurassic cover rocks, we applied a variety of depth to magnetic and gravity source estimation approaches from both line and gridded datasets. Given the huge differences between recent satellite gravity estimates of crustal thickness (Pappa et al., 2019, JGR) and sparse seismological constraints, we examine different scenarios for isostatic compensation of Rock Equivalent Topography and intracrustal loads, as a function of variable effective elastic thickness (Te) across the WSB and its flanks.  </p><p>Our models reveal a major lithospheric-scale boundary along the northeastern margin of the WSB, separating the Ross Orogen from a cryptic and composite Precambrian Wilkes Terrane. At the onset of enhanced flow for the central Cook ice stream, we image a Precambrian basement high with a felsic bulk composition. We suggest based on the similarity in potential field signatures that it represents late Paleoproterozoic to Mesoproterozoic igneous basement as exposed in South Australia, where it also associated with high GHF (80-120 mW/m<sup>2</sup>), primarily caused by anomalously radiogenic granitoids.</p><p>We hypothesise that the differences in basement depth and metasediment/sediment thickness, coupled with differences in intracrustal heat production give rise to significantly greater heterogeneity in GHF beneath different sectors of the WSB than previously recognised. To help quantify such heterogeneity we develop a suite of new probabilistic thermal models for the study region.</p>


2013 ◽  
Vol 9 (6) ◽  
pp. 2489-2505 ◽  
Author(s):  
H. Fischer ◽  
J. Severinghaus ◽  
E. Brook ◽  
E. Wolff ◽  
M. Albert ◽  
...  

Abstract. The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.


2010 ◽  
Vol 3 (1) ◽  
pp. 195-230 ◽  
Author(s):  
A. M. Le Brocq ◽  
A. J. Payne ◽  
A. Vieli

Abstract. The dataset described in this paper (ALBMAP) has been created for the purposes of high-resolution numerical ice sheet modelling of the Antarctic Ice Sheet. It brings together data on the ice sheet configuration (e.g. ice surface and ice thickness) and boundary conditions, such as the surface air temperature, accumulation and geothermal heat flux. The ice thickness and basal topography is based on the BEDMAP dataset (Lythe et al., 2001), however, there are a number of inconsistencies within BEDMAP and, since its release, more data has become available. The dataset described here addresses these inconsistencies, including some novel interpolation schemes for sub ice-shelf cavities, and incorporates some major new datasets. The inclusion of new datasets is not exhaustive, this considerable task is left for the next release of BEDMAP, however, the data and procedure documented here provides another step forward and demonstrates the issues that need addressing in a continental scale dataset useful for high resolution ice sheet modelling. The dataset provides an initial condition that is as close as possible to present-day ice sheet configuration, aiding modelling of the response of the Antarctic Ice Sheet to various forcings, which are, at present, not fully understood.


Sign in / Sign up

Export Citation Format

Share Document