scholarly journals Inferring geothermal heat flux from an ice-borehole temperature profile at Law Dome, East Antarctica

2020 ◽  
Vol 66 (257) ◽  
pp. 509-519 ◽  
Author(s):  
Laura Mony ◽  
Jason L. Roberts ◽  
Jacqueline A. Halpin

AbstractGeothermal heat flux (GHF) is an important control on the dynamics of Antarctica's ice sheet because it controls basal melt and internal deformation. However, it is hard to estimate because of a lack of in-situ measurements. Estimating GHF from ice-borehole temperature profiles is possible by combining a heat-transfer equation and the physical properties of the ice sheet in a numerical model. In this study, we truncate ice-borehole temperature profiles to determine the minimum ratio of temperature profile depth to ice-sheet thickness required to produce acceptable GHF estimations. For Law Dome, a temperature profile that is within 60% of the local ice thickness is sufficient for an estimation that is within approximately one median absolute deviation of the whole-profile GHF estimation. This result is compared with the temperature profiles at Dome Fuji and the West Antarctic Ice Sheet divide which require a temperature profile that is 80% and more than 91% of the ice thickness, respectively, for comparable accuracy. In deriving GHF median estimations from truncated temperature profiles, it is possible to discriminate between available GHF models. This is valuable for assessing and constraining future GHF models.

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.


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.


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.


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.


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>


2019 ◽  
Vol 13 (7) ◽  
pp. 2023-2041 ◽  
Author(s):  
Johannes Sutter ◽  
Hubertus Fischer ◽  
Klaus Grosfeld ◽  
Nanna B. Karlsson ◽  
Thomas Kleiner ◽  
...  

Abstract. The international endeavour to retrieve a continuous ice core, which spans the middle Pleistocene climate transition ca. 1.2–0.9 Myr ago, encompasses a multitude of field and model-based pre-site surveys. We expand on the current efforts to locate a suitable drilling site for the oldest Antarctic ice core by means of 3-D continental ice-sheet modelling. To this end, we present an ensemble of ice-sheet simulations spanning the last 2 Myr, employing transient boundary conditions derived from climate modelling and climate proxy records. We discuss the imprint of changing climate conditions, sea level and geothermal heat flux on the ice thickness, and basal conditions around previously identified sites with continuous records of old ice. Our modelling results show a range of configurational ice-sheet changes across the middle Pleistocene transition, suggesting a potential shift of the West Antarctic Ice Sheet to a marine-based configuration. Despite the middle Pleistocene climate reorganisation and associated ice-dynamic changes, we identify several regions conducive to conditions maintaining 1.5 Myr (million years) old ice, particularly around Dome Fuji, Dome C and Ridge B, which is in agreement with previous studies. This finding strengthens the notion that continuous records with such old ice do exist in previously identified regions, while we are also providing a dynamic continental ice-sheet context.


2015 ◽  
Vol 11 (9) ◽  
pp. 1165-1180 ◽  
Author(s):  
M. Willeit ◽  
A. Ganopolski

Abstract. Permafrost influences a number of processes which are relevant for local and global climate. For example, it is well known that permafrost plays an important role in global carbon and methane cycles. Less is known about the interaction between permafrost and ice sheets. In this study a permafrost module is included in the Earth system model CLIMBER-2, and the coupled Northern Hemisphere (NH) permafrost–ice-sheet evolution over the last glacial cycle is explored. The model performs generally well at reproducing present-day permafrost extent and thickness. Modeled permafrost thickness is sensitive to the values of ground porosity, thermal conductivity and geothermal heat flux. Permafrost extent at the Last Glacial Maximum (LGM) agrees well with reconstructions and previous modeling estimates. Present-day permafrost thickness is far from equilibrium over deep permafrost regions. Over central Siberia and the Arctic Archipelago permafrost is presently up to 200–500 m thicker than it would be at equilibrium. In these areas, present-day permafrost depth strongly depends on the past climate history and simulations indicate that deep permafrost has a memory of surface temperature variations going back to at least 800 ka. Over the last glacial cycle permafrost has a relatively modest impact on simulated NH ice sheet volume except at LGM, when including permafrost increases ice volume by about 15 m sea level equivalent in our model. This is explained by a delayed melting of the ice base from below by the geothermal heat flux when the ice sheet sits on a porous sediment layer and permafrost has to be melted first. Permafrost affects ice sheet dynamics only when ice extends over areas covered by thick sediments, which is the case at LGM.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document