Hostname: page-component-78c5997874-8bhkd Total loading time: 0 Render date: 2024-11-05T04:49:13.128Z Has data issue: false hasContentIssue false
  • English
  • Français

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

Published online by Cambridge University Press:  22 April 2020

Laura Mony Open the ORCID record for Laura Mony [Opens in a new window] ,
Jason L. Roberts Open the ORCID record for Jason L. Roberts [Opens in a new window]  and
Jacqueline A. Halpin Open the ORCID record for Jacqueline A. Halpin [Opens in a new window]
Laura Mony*
Affiliation:
Institute for Marine and Antarctic Studies, University of Tasmania, Private Bag 129, Hobart, Tasmania7001, Australia
Jason L. Roberts
Affiliation:
Australian Antarctic Division, Channel Highway, Kingston, Tasmania7050, Australia Antarctic Climate & Ecosystems Cooperative Research Centre, University of Tasmania, Private Bag 80, Hobart, Tasmania7001, Australia
Jacqueline A. Halpin
Affiliation:
Institute for Marine and Antarctic Studies, University of Tasmania, Private Bag 129, Hobart, Tasmania7001, Australia
*
Author for correspondence: Laura Mony, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Geothermal 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.

Keywords

Glacier geophysicsice coreice dynamicsice temperature
Type
Papers
Information
Journal of Glaciology , Volume 66 , Issue 257 , June 2020 , pp. 509 - 519
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s), 2020

Introduction

The potential of the Antarctic Ice Sheet (AIS) to raise global sea levels by up to 58 m (Fretwell and others, Reference Fretwell2013) means that understanding the impacts of a warming climate on the southern polar region is one of the grand scientific challenges of our time (Bell, Reference Bell2008). The extent of the AIS is sensitive to change in both the atmosphere and the ocean (Whitehouse and others, Reference Whitehouse, Gomez, King and Wiens2019). The vulnerability of the AIS underscores the need for more accurate ice-sheet models, these are crucial for understanding ice dynamics and future ice flow variations. Geothermal heat flux (GHF) is an important thermal boundary condition in these models because it affects the temperature at the base of the ice sheet and is the dominant control on basal temperatures in slow-flowing areas (Larour and others, Reference Larour, Morlighem, Seroussi, Schiermeier and Rignot2012; Pittard and others, Reference Pittard, Roberts, Galton-Fenzi and Watson2016). Ice flow is sensitive to the local spatial variation of underlying GHF, especially near ice divides and at the edge of ice streams (Bell and others, Reference Bell, Studinger, Shuman, Fahnestock and Joughin2007; Pittard and others, Reference Pittard, Roberts, Galton-Fenzi and Watson2016). In turn, the heat supplied to the base of the ice sheet influences the rate of internal ice deformation (Budd and Jacka, Reference Budd and Jacka1989; Durand and others, Reference Durand2007). The GHF also modulates the production of basal meltwater that can lubricate the bedrock–ice interface and lead to faster ice velocities and increase discharge (Näslund and others, Reference Näslund, Jansson, Fastook, Johnson and Andersson2005; Pattyn, Reference Pattyn2010). Local GHF is also a first-order consideration in site selection for climate records, especially in the search for ‘old’ (million-year) ice core sites. Locations where ice is at pressure melting point can lead to melt-driven erosion of the glacial stratigraphy, reducing the temporal range of the climate record (Fischer and others, Reference Fischer2013; Liefferinge and Pattyn, Reference Liefferinge and Pattyn2013; Parrenin and others, Reference Parrenin2017; Karlsson and others, Reference Karlsson2018).

GHF is the sum of heat generated by mantle convection and conducted through the lithosphere and the heat generated by radiogenic decay of heat-producing elements within the lithosphere. These elements are mainly uranium, thorium and potassium which together produce 98% of the lithospheric heat from radiogenic decay (Pollack, Reference Pollack1982). With direct access to the bedrock, GHF (q(0)) is estimated by multiplying the vertical temperature gradient of the subsurface (∂T/∂z) and the thermal conductivity of the rock (K) (Turcotte and Schubert, Reference Turcotte and Schubert2014):

(1)$$q\lpar 0\rpar = -K {\partial T\over \partial z}.$$

In Antarctica, GHF is one of the most difficult parameters to constrain because accessing subglacial bedrock and continental shelves is logistically challenging due to kilometres-thick ice cover. Few in-situ GHF measurements have been attempted, and all of them are located under the West Antarctic Ice Sheet (WAIS) in subglacial sediments (e.g. Morin and others, Reference Morin2010; Fisher and others, Reference Fisher, Mankoff, Tulaczyk, Tyler and Foley2015; Begeman and others, Reference Begeman, Tulaczyk and Fisher2017) or oceanic sediments (e.g. Dziabek and others, Reference Dziadek, Gohl, Diehl and Kaul2017) (Fig. 1). For East Antarctica, there are GHF estimates derived from ice-borehole temperature profiles (Fig. 1), but no measurements have yet been reported from bedrock or subglacial sediments.

Fig. 1. Location of deep temperature profile boreholes and different local GHF estimations. Ice-sheet thickness data from BEPMAP 2 (Fretwell and others, Reference Fretwell2013) and Basemap from Quantarctica and the Norwegian Polar Institute 1: Gow and others (Reference Gow, Ueda and Garfield1968); 2: Risk and Hochstein (Reference Risk and Hochstein1974); 3: Foster (Reference Foster1978); 4: Decker and Bucher (Reference Decker and Bucher1982); 5: Nicholls and Paren (Reference Nicholls and Paren1993); 6: Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999); 7: Hondoh and others (Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002); 8: Price and others (Reference Price2002); 9: Engelhardt (Reference Engelhardt2004); 10: Carter and others (Reference Carter, Blankenship, Young and Holt2009); 11: Morin and others (Reference Morin2010); 12: Schröder and others (Reference Schröder, Paulsen and Wonik2011); 13: Clow and others (Reference Clow, Cuffey and Waddington2012); 14: Zagorodnov and others (Reference Zagorodnov2012); 15: Mulvaney and others (Reference Mulvaney2012); 16 and 17: Carson and others (Reference Carson, McLaren, Roberts, Boger and Blankenship2014); 18: Fisher and others (Reference Fisher, Mankoff, Tulaczyk, Tyler and Foley2015); 19: Dmitriev and others (Reference Dmitriev, Bolshunov and Podoliak2016); 20: Parrenin and others (Reference Parrenin2017); 21: Begeman and others (Reference Begeman, Tulaczyk and Fisher2017); 22: Dziabek and others (Reference Dziadek, Gohl, Diehl and Kaul2017) (26 measurements).

Several studies have attempted to map Antarctic GHF using different methods, including: (1) interpolated geophysical models using either seismic velocities (e.g. Shapiro and Ritzwoller, Reference Shapiro and Ritzwoller2004; An and others, Reference An2015) or magnetic data (e.g. Maule and others, Reference Maule, Purucker, Olsen and Mosegaard2005; Purucker, Reference Purucker2013; Martos and others, Reference Martos2017), (2) geological and structural properties of the bedrock (e.g. Carson and others, Reference Carson, McLaren, Roberts, Boger and Blankenship2014; Burton-Johnson and others, Reference Burton-Johnson, Halpin, Whittaker, Graham and Watson2017), (3) the properties within and at the base of the ice sheet (e.g. Siegert and Dowdeswell, Reference Siegert and Dowdeswell1996; Shroeder and others, Reference Schroeder, Blankenship, Young and Quartini2014) and (4) a combination of available models (e.g. Liefferinge and Pattyn, Reference Liefferinge and Pattyn2013). The maps produced by these models allow us to understand the regional variation in GHF, especially the differences between East and West Antarctica, but they are unable to predict local variations due to a limited resolution typically tens to hundreds of kilometres. Indeed, the actual spatial variation of GHF can be on a scale of few kilometres due to local geological contrasts, hydrothermal circulation and volcanic activity (Carson and others, Reference Carson, McLaren, Roberts, Boger and Blankenship2014; Begeman and others, Reference Begeman, Tulaczyk and Fisher2017). Furthermore, continent-wide geophysical models use simplified lithospheric parameters including laterally homogeneous radiogenic heat production, although this is known to vary with rock type (e.g. Hasterok and others, Reference Hasterok, Gard and Webb2017) and has been estimated to contribute up to 70% of total surface GHF on the Antarctic Peninsula (Burton-Johnson and others, Reference Burton-Johnson, Halpin, Whittaker, Graham and Watson2017). The fact that GHF is currently not well-constrained is evident in the inconsistencies between different GHF models.

GHF estimates from the ice sheet

The temperature of the ice sheet at different depths is a good indicator of past surface temperature variation because, over the years, the heat from the surface is conducted and advected into the ice sheet and this leaves an imprint on the ice temperature (e.g. Dahl-Jensen and others, Reference Dahl-Jensen1998). GHF will also affect the temperature of the ice sheet by providing a source of heat at the bottom of the ice. By measuring the temperature profile of the ice sheet and compensating for heat from the bedrock, it is possible to reconstruct a smoothed version of surface temperature history and estimate GHF. However, using ice-borehole temperature profiles to estimate GHF is challenging because the heat generated at the bottom of the ice sheet can be produced from several sources, especially when the ice sheet is moving horizontally. Ice flow can produce shear heating, heat advection and basal frictional heating which could affect the temperature of the ice sheet (Engelhardt, Reference Engelhardt2004; Pittard and others, Reference Pittard, Roberts, Galton-Fenzi and Watson2016). Bedrock topography also has an impact on GHF magnitude at a specific location by concentrating the heat in an area where there is small-scale sharp variation in the slope of the bedrock, like a narrow deep valley, which focusses and increases the basal temperature (Van der Veen and others, Reference Van der Veen, Leftwich, Von Frese, Csatho and Li2007). In an environment with basal melting, this estimation would be a minimum value because the heat used to melt the ice is lost in our model (see section Assumption). These factors should be considered when GHF is modelled from borehole temperature data.

Deep ice-borehole temperature profiles have been measured across Antarctica (Fig. 1) and GHF has been derived from some of these profiles (e.g. Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999; Hondoh and others, Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002; Engelhardt, Reference Engelhardt2004; Clow and others, Reference Clow, Cuffey and Waddington2012). The differences between these ice-borehole estimations of GHF and available geophysically-derived GHF models are significant. For example, Clow and others (Reference Clow, Cuffey and Waddington2012) calculated a value of 240 mW m−2 (Fig. 1) in the central portion of the WAIS, this is higher than the values modelled by any geophysically-derived models. Martos and others (Reference Martos2017) estimated the highest value for this area at approximately 120 mW m−2. However, Clow and others (Reference Clow, Cuffey and Waddington2012) did not account for horizontal flow effects related to bedrock topography and the GHF estimation could have been lower than the one they produced. Figure 2 shows the differences between estimated GHF from various studies at locations where deep ice-borehole temperature profiles are available (Fig. 1). GHF estimations for the East Antarctic Ice Sheet (EAIS) display a relatively narrow range, standard deviations are between 5.9 mW m−2 for South Pole/Vostok and 7.9 mW m−2 for Law Dome/Dome Fuji. The WAIS displays a larger range of GHD estimates, especially the WAIS divide and Byrd Station, where the standard deviation are 24.5 and 25.1 mW m−2, respectively.

Fig. 2. GHF values from available continent-wide GHF maps for deep ice borehole locations in Antarctica (see Fig. 1).

This study provides an understanding of the variation of the GHF estimation as a function of the depth of ice-borehole temperature profile. Here we analyse a 1177 m ice-borehole temperature profile from Law Dome and derive a revised GHF estimate using a heat diffusion-advection equation in the ice sheet integrated into a numerical model which includes a forward and an inverse model. The use of a two-step numerical model to estimate surface temperature history from an ice-borehole temperature profile has been performed previously (e.g. Cuffey and others, Reference Cuffey1995; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999; Muto and others, Reference Muto, Scambos, Steffen, Slater and Clow2011; Orsi and others, Reference Orsi, Cornuelle and Severinghaus2012; Roberts and others, Reference Roberts2013). We then systematically truncate the temperature profile from the base upwards, sequentially calculating a new GHF and its associated uncertainties, until the shallowest hundred metres. We explore how the accuracy of modelled GHF changes with depth, including the minimum depth required to discriminate among available geophysical GHF models at the Law Dome borehole location (Dome Summit South; DSS) (Fig. 2). We then explore the regional applicability of this analysis with comparisons to the Dome Fuji and WAIS divide ice core sites (WDC) (Fig. 1). Finally, the impact of the different assumptions used in the model on the final GHF estimation is assessed.

Measurements

Law Dome

To mitigate the impact of horizontal flow, the Dome Summit South (DSS) profile has been measured 4.6 km from the summit of the Dome (Fig. 1), where the horizontal velocity of the ice sheet is assumed to be minimal and the environment approximates a steady state (Morgan and others, Reference Morgan1997; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999). Furthermore, the DSS borehole has been drilled in an area of relatively level bedrock which limits the impact of topography on GHF measurement (Morgan and others, Reference Morgan1997). In 1996, when the DSS temperature profile was measured, the mean annual air temperature was -21.8oC, this precludes melting of the snow at the surface during summer. The recent average accumulation rate was 0.679 m ice a−1 (Morgan and others, Reference Morgan1997). The thickness of the ice sheet is 1220 ± 25 m (Morgan and others, Reference Morgan1997) and the temperature profile stopped around 20 m above the bedrock (Fig. 3). The ice equivalent thickness of the profile is 1176.8 m (Morgan and others, Reference Morgan1997; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999).

Fig. 3. Ice-borehole temperature profiles for Law Dome (DSS site) (Morgan and others, Reference Morgan1997), Dome Fuji (Hondoh and others, Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002) and West Antarctic Ice Sheet Divide (Clow and others, Reference Clow, Cuffey and Waddington2012).

Dome Fuji

At Dome Fuji, a temperature profile has also been measured on a topographic high and this profile reached the bedrock at an ice thickness of 3090 m. Since 2002, the average accumulation rate is 0.032 m ice a−1 (Hondoh and others, Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002).

West Antarctic Ice Sheet divide

The WDC temperature profile is located 24 km away from an ice divide where it is also considered to have minimal horizontal velocity. A horizontal velocity of 3 m a−1 has been measured by Koutnik and others (Reference Koutnik2016) and the impact of this on the GHF estimation is discussed in the Assumptions section. The deepest temperature recorded was at 3327 m and ice was recovered up to 3405 m (Buizert and others, Reference Buizert2015). The accumulation rate is 0.3 m ice a−1 (Koutnik and others, Reference Koutnik2016).

Steady-state environment

The assumption of a steady-state flow environment has been tested after the modelling at the three locations by comparing the observed surface to basal temperature difference with a steady-state theoretical thermal estimate (Cuffey and Paterson, Reference Cuffey and Paterson2010),

(2)$$\Delta T = {{\rm GHF} \times h \times \theta _B\over K}\comma\; $$

where ΔT is the temperature difference between the top and the bottom of the temperature profile, h is the vertical distance in metres and K is the thermal conductivity of the ice. θB is a dimensionless parameter for temperature which gives the temperature difference between the surface and the base of the ice sheet using the following equation:

(3)$$\theta_{B}=\left[{\pi\over 2\gamma}\right]^{{1}/{2}}{\rm erf}\left({\gamma\over 2}\right)^{{1}/{2}}\comma\; $$

where γ is the advection parameter (Cuffey and Paterson, Reference Cuffey and Paterson2010). K and θB have been calculated using the average temperature (T) measured by the different temperature profiles. The temperatures measured at the base and the surface are also derived from the same profiles. The GHF used in this equation is the one that has been estimated by our model.

At Law Dome, for a perfect steady-state environment, the temperature difference between the beginning and the end of the profile must be 10.42 °C, which is lower than the measured temperature difference of 14.3 °C. It shows that the environment is not in a perfect steady state, but the two values are similar and we use this to justify an assumption of a steady-state flow environment within our model. For Dome Fuji the calculated temperature difference is 49.16 °C which is close to the measured difference of 55.3 °C. For the WDC site, the calculated temperature difference is 32.62 °C and the measured difference is 21.29 °C, clearly the assumption of steady state is less valid at the WDC site. Even if this assumption is not entirely valid, we use it to simplify our model. We assess the impact of this assumption on the model results in the section Assumptions.

Method

Heat flow model

To estimate GHF values, a numerical model has been developed which combines: (1) a forward numerical model that reconstructs an ice temperature profile from a given surface temperature history, GHF, accumulation rate and (2) an inverse model which updates the surface temperature history and GHF to minimise the error between the measured and estimated temperature profile. To analyse how shallow a temperature profile can be while still allowing a sufficiently accurate estimation of the GHF value, we use a temperature profile systematically truncated at increasingly shallower intervals. The top 100 m of the observed temperature profile is not included because it is dominated by recent changes in climate and accounting for this introduces too many additional free parameters in the inverse model with associated significant increases in computational cost.

To solve how heat flows through and in the ice, a series of prescribed surface temperature combinations and GHF values are used as boundary conditions in a heat flow model. A one-dimensional finite difference approximation to the heat diffusion-advection equation is used to estimate the temperature profile in the ice. To minimise the difference between the simulated temperature profile and observed one, the inverse model uses particle swarm optimisation with polishing (see section Inverse model) to minimise the value of the regularised root mean squared error (E rms):

(4)$$E_{\rm rms} = \sqrt{\sum {\lpar T_{\rm m} \lpar i\rpar -T_{{\rm o}} \lpar i\rpar \rpar ^2\over n}}+\lambda \sum\left({\partial ^2 T\over \partial t ^2}\right)^2\comma\; $$

where λ is the regularisation parameter, T m is the temperature simulated by the model and T o the temperature measured in the borehole in kelvin, n is the number of data points and ∂2T/∂t 2 is the second derivative of surface temperature with respect to time. The first term is the usual RMS error and the second term is the regularisation term which enforces a degree of smoothness in the model surface temperature function by reducing the impact of overfitting values. Overfitting values could cause problems because they tend to capture local signal variations and noise in the dataset obscuring a general trend.

The change of temperature over time in the ice sheet is controlled by diffusion, advection, firn densification and a correction to the diffusive term due to variations in the thermal conductivity (Cuffey and Paterson, Reference Cuffey and Paterson2010). The heat flow equation used to simulate the variation in temperature of the ice through time,

(5)$${\partial T\over \partial t}={K\over \rho C}{\partial ^2 T\over \partial z^2}- w_{z\comma t} \left({\partial T\over \partial z}\right)+{\,f\over \rho C}+ {1\over \rho C}{\partial K \over \partial z}{\partial T\over \partial z}\comma\; $$

does not include horizontal advection because sites chosen have minimal horizontal velocity (Johnsen and others, Reference Johnsen, Dahl-Jensen, Dansgaard and Gundestrup1995; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999).

The specific heat capacity (C) of the ice in J kg−1 K−1,

(6)$$C=152.5+7.122T\comma\; $$

is dependent on the temperature (T) in kelvin. The firn densification energy term (f) (J m−3 s−1) is given by:

(7)$$f={wP\over \rho}\left({{\rm d}\rho\over {\rm d}z}\right).$$

The thermal conductivity, in W m−1 K−1, is also dependent on the temperature. It follows Eqns (8) and (9) for the firn (Cuffey and Paterson, Reference Cuffey and Paterson2010), where K i is the thermal conductivity of the pure ice and ρ is the density in kg m−3,

(8)$$K=2.1\times10^{-2} + \lpar 4.2\times 10^{-4}\rpar \rho + 2.2 \times 10^{-9}\rpar \rho ^3$$
(9)$$K = {2K_{\rm i} \rho\over \lpar 3\rho _{\rm i}-\rho\rpar }\comma\; $$

and Eqn 10 for the ice (Cuffey and Paterson, Reference Cuffey and Paterson2010):

(10)$$K=9.828\exp\lpar -5.7 \times 10 ^{-3} T\rpar.$$

The vertical flow of the ice for DSS site equation (w z,t) is from Dansgaard and Johnsen (Reference Dansgaard and Johnsen1969) and is based on a vertical strain rate which is constant in the upper part of the ice sheet, then decreases linearly towards the bedrock (Morgan and others, Reference Morgan1997; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999):

(11)$$w_{z\comma t}= \left\{\matrix{ -a\lpar t\rpar {z - h/2\over H-h/2} \hfill & z \in \lsqb h\comma\; H\rsqb \hfill \cr -a\lpar t\rpar {z^2\over 2h\lpar H-h/2\rpar } \hfill & z \in \lsqb 0\comma\; h\rsqb \comma\; \hfill }\right.$$

where H is the thickness of the ice sheet and is assumed to be constant in time. The resulting vertical velocity profile is only dependent on the accumulation rate (a(t)). The accumulation rate for the DSS site, which varies in time depending on the surface temperature (T sur(t)), is based on Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999) model:

(12)$$a\lpar t\rpar = a_0 \exp\lpar \alpha\lsqb T_{\rm sur}\lpar t\rpar - T_0\rsqb \rpar \comma\; $$

where the recent average accumulation rate (a 0) and the temperature at the surface (T 0) at the time the temperature profile was measured are 0.69 m ice a−1 and − 21.8 °C, respectively. α is the change of accumulation per °C. According to Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999) this value is around 7% for the DSS site, which falls in between the α range in general global circulation models (3–4%) and the one for polar models (10%). We used the same α for Dome Fuji and WDC, because we considered both sites to be in a similar environment.

Initial and boundary conditions

The thermal boundary conditions of the forward model are the GHF and prescribed surface temperature history for a period of 50 ky at the surface of the ice sheet. We use a piecewise linear surface temperature history function composed of five temperature values to produce a general trend of the surface temperature history and to mitigate the impact of non-unique solutions due to the presence of high-frequency variation of surface temperature through time. During a certain time period, the surface temperature diffuses in the ice sheet and affects the temperature profile. At deeper locations, the effect of surface temperatures on englacial temperatures is less intense and the initial high-frequency variations spread out and can be produced by a wide range of values (Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999; Roberts and others, Reference Roberts2013).

The bedrock is considered as a granite and its properties are listed in Table 2. The type of rock used in the model and its impact on the final estimation of GHF at Law Dome are assessed in the discussion in section Assumptions. The GHF is a boundary condition at the bottom of a layer of bedrock 1300 m thick where it is assumed to be the only source of heat and the surface temperature variation does not have an impact. To calculate the temperature profile in the 1300 m thick layer of bedrock (Table 2), Eqn (1) is used. GHF of the rock underneath the ice sheet (q(0)) is going to be updated by the inverse model to produce the best temperature profile.

Table 1. Parameters used for testing the steady-state assumption

Table 2. Input values for the bedrock

To produce a robust result without being computationally expensive, the model is run several times to optimise the time step, the grid cell size in the ice, the regularisation term (λ) and the number of repetitions during each simulation (Table 3). We picked the value that gave robust results with an acceptable Erms. For λ, the best value is obtained when there is a balance between the regularised and the standard error.

Table 3. Optimised input values in the ice

Inverse model

In the inverse model, the optimisation is divided into two main steps which are particle swarm optimisations and subsequent steepest descent polish. Particle swarm optimisation improves a set of solution according to a list of criteria by moving particles towards the best solution possible by iteration. It is used to update prescribed temperature history and the GHF value, and to minimise the E rms (Eqn (4)). Particle swarm optimisation allows a good investigation of the search space with a focus on areas that produce the best result possible. For this study, we use 8 swarm members and 250 iterations where each member is assigned a velocity vector which has a random and a momentum component. After evaluating all the criteria for the actual data set, the velocity is adjusted for the next iteration according to how far the particle is from the best solution. The temperature history and GHF that produce the lowest E rms is kept if the E rms is below the error threshold of 0.1 °C, if not the whole iteration process restarts with the same initial boundary conditions. To further minimise E rms we polish the results using direct descent optimisation from the particle swarm optimisation solution results. The steepest descent polish finds the nearest local minimum of a function by computing a new gradient at every iteration to help the model converge. This optimisation is done 100 times, or until the difference between the previous and actual E rms is sufficiently small.

Discussion

For each ice core site, we truncated ice-borehole temperature profiles at different depths. This allowed us to produce several probability distributions of GHF values using a one-dimensional finite difference approximation of the heat diffusion-advection equation. We explore the impact of the depth of a measured ice-sheet temperature profile on the median GHF estimation at specific locations.

Dome Summit South, Law Dome

At DSS, using the whole temperature profile (LD100), the median GHF estimation is 72.0 mW m−2 with a median absolute deviation (MAD) of ± 3.2 mW m−2 (Fig. 4). MAD is a robust estimation of the variability of a probability distribution because it is less impacted by outliers than the standard deviation (SD), which is 1.4826 × MAD if the population follows a normal distribution. We also calculate the percentage of difference between values by using the following equation:

(13)$$\percnt _{\rm diff}={\vert x_1-x_0\vert \over \lpar x_1+x_0\rpar /2}.$$

The LD100 value is 4.21% lower than the one estimated by Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999) of 75.1 mW m−2, but the difference between the two estimations is within ±1 MAD of LD100 GHF estimation. The differences between the two models are the optimisation technique and the surface boundary conditions. The surface temperature history used in our model follows a piecewise linear function with five points which, combined with the regularisation term added in the E rms equation, decreases the impact of sharp, but short, variation of temperature by smoothing the surface temperature history.

Fig. 4. GHF estimation for DSS site from 260 iterations of our model by using a temperature profile of 1176.8 m ice equivalent depth.

To explore the effect of temperature profile depth on acceptable GHF median values we compare the different probability distributions from truncated temperature profiles with the LD100. With a temperature profile depth as shallow as 57% of the thickness of the ice sheet (LD57), the median GHF is 75.5 mW m−2 (Fig. 5). We consider this value to be the shallowest acceptable because it is at the limit of ±1 MAD of the LD100 estimation and more than 60% of the two density distributions overlap. A temperature profile which is half the thickness of the ice sheet gives an estimation of GHF that is 8.13% of LD100 estimation and is within ±2 MAD. The difference increases with shallower temperature profiles, reaching ±9.5 MAD for a temperature profile of 20% total ice thickness (LD20). With LD20, the difference with the GHF median value using the full temperature profile (LD100) is around 35.5% and only 7% of the two probability distributions overlap, indicating that few GHF values from the LD20 models are considered acceptable (Fig. 6).

Fig. 5. (A) GHF median values with their MAD for the DSS site using temperature profile truncated at different depths. (B) GHF estimation from Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999), Martos and others (Reference Martos2017), Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013), Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004).

Fig. 6. Comparison between the probability distribution at DSS for the entire temperature profile (LD100) and 20% of the temperature profile (LD20).

Regional applicability

To analyse the regional applicability of the DSS model, we use a modified model at Dome Fuji and WDC to estimate GHF. The ice-borehole temperature profile of Dome Fuji and WDC is located at the summit of a dome and near an ice divide, respectively, where the horizontal velocity is considered to be zero (Fig. 1). The main difference between the DSS model and the one used for the two other locations is the vertical velocity equation:

(14)$$w_{z\comma t}={-a\lpar t\rpar \ast z\over H}.$$

Also, the E rms thresholds for both locations are higher than DSS. The Dome Fuji threshold is 0.35 °C because the uncertainty of the measured ice-borehole temperature data is approximately 20 times higher for Dome Fuji than for the DSS site. For WDC, it is between 0.9 and 0.2oC, depending on the length of the profile. The WDC site is far enough from the ice divide that some of the assumptions used in the model are not entirely valid, especially the assumption of a steady-state environment.

For DSS, as the temperature profile is truncated, the median GHF shifts to a higher value and the MAD increases; this makes the value less reliable (Fig. 5). At Dome Fuji, GHF estimations decrease with increasing MAD values for shallower temperature profiles (Fig. 7). At WDC, GHF median estimations shift to lower values until approximately half of the temperature profile is used, and then shift to higher GHF values with an MAD following the same behaviour (Fig. 8). The median GHF estimated by our model using the whole temperature profile at Dome Fuji (DF100) is 50.4 mW m−2 with an MAD of 0.9 mW m−2, this is ≈16% lower than that estimated by Hondoh and others (Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002) of 59 mW m−2. The difference between the two values could be explained by the difference in the method used in our model and that of Hondoh and others (Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002). In fact, Hondoh and others (Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002) used a local δ18O record to estimate a GHF value, instead of using the ice-borehole temperature profile. For Dome Fuji, the temperature profile must be at least 90% of the thickness of the ice sheet (DF90) to give an acceptable GHF median value (within ±1 MAD of the DF100 value). For WDC, the median GHF for the whole profile (WDC100) is 90.5 mW m−2 with an MAD of 1.2 mW m−2 and 91% of the whole temperature profile is required in order to be within ±1 MAD of WDC100 value. However, the MAD for DF100 and WDC100 is smaller than LD100. Using the LD100 MAD of 3.2 mW m−2, Dome Fuji requires a temperature profile of at least 80% of the total ice thickness and more than 91% for WDC.

Fig. 7. (A) GHF median values with their MAD for Dome Fuji by using temperature profile truncated at different depths. (B) GHF estimation from Martos2017: Martos and others (Reference Martos2017), Liefferinge2013: Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013), Shapiro2004: Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004), Hondoh2002: Hondoh and others (Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002).

Fig. 8. (A) GHF median values with their MAD for WAIS divide by using temperature profile truncated at different depths. (B) GHF estimation from Martos2017: Martos and others (Reference Martos2017), Liefferinge2013: Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013), Shapiro2004: Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004).

The median GHF behaviour is different at the three locations. The change in surface temperature and accumulation through time has an impact on the ice-sheet temperature profile which might not have reached equilibrium yet (Engelhardt, Reference Engelhardt2004). For example, in an environment with no horizontal velocity, a constant thickness and an accumulation rate of 0.07 m a −1, Cuffey and Paterson (Reference Cuffey and Paterson2010) show than 10 ka after an abrupt but persistent increases of the surface temperature of 15 °C, the ice sheet still contains a signature of the temperature variation at 1000 m below the surface. However, an increase of the accumulation rate to 0.3 m a −1 increases the speed of the propagation of the surface temperature in the ice sheet and the imprint is visible ≈2000 m below surface (Cuffey and Paterson, Reference Cuffey and Paterson2010). At Law Dome, the ice sheet is ≈1200 m thick and the accumulation rate is higher than 0.3 m a −1 (Morgan and others, Reference Morgan1997). The DSS results are considered to be robust because, according to Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999), there is no remnant of cold from the past glacial cycle at the bottom of the ice sheet. In several studies δ18O in the ice is used as a proxy for past surface temperature (e.g. Hondoh and others, Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002). This is achieved by measuring the isotopic ratio oxygen/hydrogen in the ice core and comparing it with a standard value (e.g. Alley and Bentley, Reference Alley and Bentley1988; Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999; Watanabe and others, Reference Watanabe2003). Using an ice-borehole temperature profile, instead of δ18O, to analyse past temperature history, limits error to the differing attenuation/suppression of high-frequency information between observations and modelling associated with the calibration of the observed data and past temperature. In contrast, the ice sheets at Dome Fuji and WDC are greater than 3000 m thick (>2.5 times thickness of Law Dome) and likely contain the imprint of multiple glacial cycles that still affects the temperature of the ice at depth. The cold temperatures of the last glacial cycle might still have an impact on the temperature profile at these two locations which would serve to lower temperature values below those assumed in a steady-state environment. In low accumulation rate environments like Dome Fuji (Hondoh and others, Reference Hondoh, Shoji, Watanabe, Salamatin and Lipenkov2002), the vertical strain rate is minimal, and this lowers the advection velocity in the ice sheet. In contrast, the WDC temperature advection velocity is faster, due to a higher accumulation rate (Koutnik and others, Reference Koutnik2016). Consequently, the impact of past glacial cycles decreases more rapidly. At Dome Fuji and WDC, the ice-sheet temperature profiles are more affected by downwelling surface perturbation than at Law Dome, which in turn affects the steady-state thermal regime of the ice sheet at these locations. In the absence of glacial cycles, the reduced advection of cold from the surface at Dome Fuji compared to Law Dome and WDC would result in the GHF dominating the vertical temperature distribution at relatively shallower depths. However, the large glacial cycle variations in surface temperature leave an imprint at Dome Fuji, complicating this interpretation. These factors account for the differences between the three locations and explain why Dome Fuji and WDC require a deeper temperature profile to produce an acceptable GHF value. The estimation of GHF from truncated temperature profiles should therefore be undertaken with caution. Each profile, depending on the maximum depth reached, does not contain the same information. Especially, the surface temperature history and the initial conditions due to the impact of past glacial cycles and intense variation of temperature that could have affected the temperature profile in the ice (Beltrami and others, Reference Beltrami, Smerdon, Matharoo and Nickerson2011).

Assumptions

Our model is a transient model with six main assumptions, these are: (1) thickness does not change over time; (2) no change of ice/snow density over time; (3) vertical velocity is dependent on the accumulation rate only; (4) there is no melt at the bottom of the ice sheet; (5) the environment is at a steady state with no horizontal velocity; and (6) the bedrock is only a granite. To assess the robustness of the results, it is important that the assumptions are approximately valid for each location. Firstly, the ice-sheet thickness might have changed over time, especially in coastal areas which are more sensitive to climate variation (Huybrechts, Reference Huybrechts1990; Levy and others, Reference Levy2016). The whole WAIS thickness is more affected by variations in the climate than the EAIS because the bedrock is mainly below sea level which makes it more vulnerable to change in temperature (Huybrechts, Reference Huybrechts1990; Ritz and others, Reference Ritz, Rommelaere and Dumas2001; Rintoul, Reference Rintoul2018). Law Dome and WDC are in areas that might have been more impacted by climate variation than Dome Fuji because of their locations. However, because the thickness variation of the ice sheet in the past, especially for Law Dome, is in the order of tens of metres, this impacted the model outputs minimally (Huybrechts, Reference Huybrechts1990; Ritz and others, Reference Ritz, Rommelaere and Dumas2001). Secondly, the impact of the ice/snow density variation is minimal in our model because, even if the firn densification is considered in our equation, the temperature profiles we use start below the firn-ice boundary. Thirdly, the rate at which the temperature will sink in the ice sheet is dependent on the accumulation rate, which in turn is affected by the surface temperature. Like the ice thickness, the impact of climate variability, which will affect the accumulation rate, is more important at the coast and for the WAIS. Fourthly, we assume that the temperature at the bottom of the ice sheet is below the pressure melting point. Where the basal ice is at or above pressure melting point, the GHF estimation is only a minimum value, since the latent heat used to melt the ice is lost in the assumption of our model. The effect of each millimetre per year of melting at the base of the ice sheet a year on the GHF estimation is 9.7 mW m−2 a−1. This value is higher than the MAD of any of our locations, which could explain the shifting in some of our values, especially at WDC. At Law Dome, the basal ice is frozen to bedrock (Dahl-Jensen and others, Reference Dahl-Jensen, Morgan and Elcheikh1999). However, this is probably not the case for Dome Fuji and WDC. Koutnik and others (Reference Koutnik2016) report that ice sheet below WDC is probably melting at a rate of few millimetres a year and according to Motoyama and others (Reference Motoyama2008), the ice at the base of Dome Fuji seems to be at pressure melting point, which is also shown in Fig. 3. Fifth, we are assuming that the ice boreholes are drilled in areas where there is minimal horizontal velocity. This is an important assumption because a high horizontal velocity produces shear heating which could affect the temperature profile in the ice, and thereby, the GHF estimation. This assumption has also been made by Dahl-Jensen and others (Reference Dahl-Jensen, Morgan and Elcheikh1999) for the DSS site, even if the surface ice velocity is 2.9 m a−1 (Morgan and others, Reference Morgan1997). It seems representative of the ice state at DSS site because of the frozen basal ice which limits the production of shear heating at the interface ice/rock. Like DSS, WDC site has a horizontal velocity of 3 m a−1, but with an ice sheet at pressure melting point. Dome F has a horizontal velocity in between 0 and 0.2 m a−1, which is considered minimal. Finally, we are assuming that the rock underneath the ice boreholes is granite. The estimated temperature profile in the bedrock, which gives the temperature at the interface ice–bedrock, is affected by the thermal conductivity. After running the DSS model with a gabbro, which has a higher heat capacity (980 J kg−1 K−1) and density (3000 kg m−3), but a lower thermal conductivity (2.57 W m−41 K−1) (Schön, Reference Schön2011) than the granite, we note that the GHF estimation increases with a difference with the granite estimations between 0.6 and 0.8 mW m−2, which is within an MAD of the GHF estimation using granite.

Comparison with geophysically-derived maps

We compare our GHF estimations with Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004), Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013) and Martos and others (Reference Martos2017) who modelled GHF and produced maps of the whole Antarctic continent. In their maps, the standard deviation of the GHF can be as high as the GHF estimation itself, which makes some of the estimations unreliable. The WAIS is a critical region with a standard deviation that could be as high as 80 mW m−2 (Shapiro and Ritzwoller, Reference Shapiro and Ritzwoller2004; Liefferinge and Pattyn, Reference Liefferinge and Pattyn2013; Martos and others, Reference Martos2017) but, with an average ice thickness of 2.5 km (Fretwell and others, Reference Fretwell2013), it is almost impossible to get in-situ GHF measurements in the bedrock, except where there are outcrops. Our model allows us to assess the borehole depth required to discriminate between the different GHF maps from available interpolated geophysical models. For Law Dome, the model of Martos and others (Reference Martos2017) yields a GHF estimation that is associated with the smallest MAD (percentage difference of 16%, Table 4). With our model, a temperature profile that is around 42% of the whole ice thickness (LD42) is enough to produce a GHF estimation that is within the same range of error (Fig. 5). However, the value from our model is more than 3 MADs away from LD100 GHF median value. Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004) and Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013) accept a higher error (Table 4) and, with our model, only ≈35% of the whole temperature profile is necessary to estimate a GHF which has the same difference (Fig. 5). For Dome Fuji, the lowest accepted error is 21% for the Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013) model, which corresponds to ≈60% of the entire temperature profile using our model (Fig. 7). Finally, for WDC, the percentage difference from the three maps are higher than 30% with a minimum value of 31% for Martos and others (Reference Martos2017) which, with our model, corresponds to a temperature profile of approximately 45% of the whole profile (Fig. 8). As explained above, the WDC temperature profile is more difficult to analyse because of the shift in the variation of GHF depending on the depth of the temperature profile. Like Pattyn (Reference Pattyn2010) and Martos and others (Reference Martos2017), it could be valuable to integrate ice-borehole temperature profiles which will produce localised GHF estimations, to optimise the regional-scale model.

Table 4. Mean GHF estimations and Median Absolute Deviation (MAD) for Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004), Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013) and Martos and others (Reference Martos2017)

aMedian Absolute Deviation.

bPercentage of difference between the mean GHF and the MAD from the map using Eqn (12).

Notes For this study, the standard deviation values from the models are converted into MAD.

Implications

We have shown that ice-borehole temperature profiles can be used to estimate local GHF values and to discriminate between available GHF maps produced by geophysical models. Our approach could also be useful in assessing the location of potential ‘old’ ice core sites by verifying the required assumption of no melting at the bottom of the ice sheet. In the aim of finding sites that contain million-year-old ice, it is important that melting has not occurred at the base of the ice sheet in the past million years (Fischer and others, Reference Fischer2013; Liefferinge and Pattyn, Reference Liefferinge and Pattyn2013; Parrenin and others, Reference Parrenin2017; Karlsson and others, Reference Karlsson2018). Even if it is possible to identify the current presence of melt water at the ice/bedrock interface by using radio echo sounding (e.g. Pattyn, Reference Pattyn2010), understanding past subglacial conditions of the AIS is possible only by modelling past variations in ice-sheet basal temperature. In regions with negligible dissipative heating, which is a common characteristic of most ice core sites, the main factors that influence the basal temperature are accumulation rate, surface temperature, thickness of the ice sheet and GHF. Even with a temperature profile that is not complete, our model gives a range of likely GHF values. With a range of acceptable GHF produced by the model, the surface temperature history, the variation of ice thickness and the variation of accumulation rate through time, it is possible to produce a range of likely basal temperatures that a specific location could have experienced in the past. Then, it would be possible to assess if the basal temperature had reached the pressure melting point in the past.

Conclusion

Our model is able to assess the variation in GHF estimations from ice-borehole temperature profiles truncated at different depths. At our focus location near the top of Law Dome (DSS), our models suggest that only 57% of the available temperature profile is required to produce an acceptable GHF value (i.e. within 1 MAD of LD100 GHF median value). However, by comparing with Dome Fuji and WDC, which require as much as 80 and 90% of the available temperature profiles, respectively, for the same level of accuracy, we notice that several external factors have an impact on our estimations. First, the six assumptions used in the model, which was designed primarily for the DSS site, could have an impact on the model results for Dome Fuji and WDC, especially the assumption of no horizontal velocity. It is important to build a suitable model for each location. Additionally, depending on the thickness of the ice sheet and the accumulation rate, the temperature of the ice sheet could be affected by past glacial cycles. The WDC ice core site, which arguably does not meet the different assumptions of our model and where the ice sheet is >3000 m thick, produces GHF estimations with a more erratic behaviour when the profile is truncated than the other two locations because WDC seems to have the most complex relationship between temperature and depth.

Depending on the required application of the GHF estimation, we show that it is possible to use relatively shallow temperature profiles to produce satisfying results. At Law Dome, for a GHF value that is less than 10% different from the median GHF using the whole 1177 m, only 50% of the temperature profile is necessary. Comparison of our modelled GHF median and MAD values with those from Shapiro and Ritzwoller (Reference Shapiro and Ritzwoller2004), Liefferinge and Pattyn (Reference Liefferinge and Pattyn2013) and Martos and others (Reference Martos2017) shows that the percentage of the individual temperature profiles (compared to total profile depth available) required to discriminate between the models at the borehole locations are 42% at Law Dome, 60% at Dome Fuji and ≈45% at the WDC.

Finally, by producing a range of likelihood GHF values, our model could be used for other purposes like confirming the assumption that there was no melting at the base of the ice sheet in the past. The values in our method open up the opportunity for fairly accessible GHF estimations across wide areas of ice cover to enable more representative estimates and assess spatial variation, and for targeted use in regions of interest (e.g. old ice sites). Future direct observations of GHF are necessary to produce maps that will be able to represent the actual spatial variability of GHF and minimise the error on the estimations. The use of shallow ice-borehole temperature profiles, and the modelling approach outlined here could aid this process, avoiding the logistical challenges of drilling deep boreholes in ice and bedrock.

Acknowledgments

We thank the members of the ice-sheet group at IMAS, especially Adam Treverrow and Felicity McCormack, for reviewing and improving an earlier version of the manuscript. This research was supported under Australian Research Council's Special Research Initiative for Antarctic Gateway Partnership (Project ID SR140300001).

References

Alley, RB and Bentley, CR (1988) Ice-core analysis on the Siple Coast of West Antarctica. Annals of Glaciology 11, 1–7.CrossRefGoogle Scholar
An, Mand 8 others (2015) Temperature, lithosphere-asthenosphere boundary, and heat flux beneath the Antarctic Plate inferred from seismic velocities. Journal of Geophysical Research: Solid Earth 120(12), 8720–8742. doi: 10.1002/2015JB011917Google Scholar
Begeman, CB, Tulaczyk, SM and Fisher, AT (2017) Spatially variable geothermal heat flux in West Antarctica: evidence and implications. Geophysical Research Letters 44, 9823–9832. doi: 10.1002/2017GL075579.CrossRefGoogle Scholar
Bell, RE (2008) The role of subglacial water in ice-sheet mass balance. Nature Geoscience 1(5), 297–304. doi: 10.1038/ngeo186CrossRefGoogle Scholar
Bell, RE, Studinger, M, Shuman, CA, Fahnestock, MA and Joughin, I (2007) Large subglacial lakes in East Antarctica at the onset of fast-flowing ice streams. Nature 445(7130), 904–907. doi: 10.1038/nature05554CrossRefGoogle ScholarPubMed
Beltrami, H, Smerdon, JE, Matharoo, GS and Nickerson, N (2011) Impact of maximum borehole depths on inverted temperature histories in borehole paleoclimatology. Climate of the Past 7(3), 745–756. doi: 10.5194/cp-7-745-2011CrossRefGoogle Scholar
Budd, W and Jacka, T (1989) A review of ice rheology for ice sheet modelling. Cold Regions Science and Technology 16(2), 107–144. doi: 10.1016/0165-232X(89)90014-1CrossRefGoogle Scholar
Buizert, Cand 9 others (2015) Precise interpolar phasing of abrupt climate change during the last ice age. Nature 520(7549), 661–665. doi: 10.1038/nature14401Google Scholar
Burton-Johnson, A, Halpin, J, Whittaker, JM, Graham, FS and Watson, SJ (2017) A new heat flux model for the antarctic peninsula incorporating spatially variable upper crustal radiogenic heat production. Geophysical Research Letters 44(11), 5436–5446. doi: 10.1002/2017GL073596CrossRefGoogle Scholar
Carson, CJ, McLaren, S, Roberts, JL, Boger, SD and Blankenship, DD (2014) Hot rocks in a cold place: high sub-glacial heat flow in East Antarctica. Journal of the Geological Society 171(1), 9–12. doi: 10.1144/jgs2013-030CrossRefGoogle Scholar
Carter, SP, Blankenship, DD, Young, DA and Holt, JW (2009) Using radar-sounding data to identify the distribution and sources of subglacial water: application to Dome C, East Antarctica. Journal of Glaciology 55(194), 1025–1040. doi: 10.3189/002214309790794931CrossRefGoogle Scholar
Clow, G, Cuffey, K and Waddington, E (2012) High Heat-Flow Beneath the Central Portion of the West Antarctic Ice Sheet. Abstract C31A-0577, Presented at the 2012 AGU Fall Meeting, San Francisco, CA.Google Scholar
Cuffey, KMand 5 others (1995) Large Arctic temperature change at the Wisconsin-Holocene glacial transition. Science (New York, N.Y.) 270(5235), 455–458. doi: 10.1126/science.270.5235.455CrossRefGoogle Scholar
Cuffey, KM and Paterson, WSB (2010) The physics of glaciers, 4th ed., Amsterdam, Academic Press, pp. 399–428.Google Scholar
Dahl-Jensen, Dand 6 others (1998) Past temperatures directly from the Greenland ice sheet. Science (New York, N.Y.) 282(5387), 268–271. doi: 10.1126/science.282.5387.268CrossRefGoogle ScholarPubMed
Dahl-Jensen, D, Morgan, VI and Elcheikh, A (1999) Monte Carlo inverse modelling of the Law Dome (Antarctica) temperature profile. Annals of Glaciology 29(1), 145–150. doi: 10.3189/172756499781821102CrossRefGoogle Scholar
Dansgaard, W and Johnsen, S (1969) A flow model and a time scale for the ice core from Camp Century, Greenland. Journal of Glaciology 8(53), 215–223.CrossRefGoogle Scholar
Decker, E and Bucher, G (1982) Geothermal studies in the Ross Island-Dry Valley region. Antarc Geoscience 4, 887–894.Google Scholar
Dmitriev, AN, Bolshunov, AV and Podoliak, AV (2016) Assessment of ice borehole temperature conditions at interface with subglacial Lake Vostok (Antarctica). International Journal of Applied Engineering Research 11(11), 7230–7233.Google Scholar
Durand, Gand 8 others (2007) Change in ice rheology during climate variations–implications for ice flow modelling and dating of the EPICA Dome C core. Climate of the Past 3(1), 155–167.CrossRefGoogle Scholar
Dziadek, R, Gohl, K, Diehl, A and Kaul, N (2017) Geothermal heat flux in the Amundsen Sea sector of West Antarctica: new insights from temperature measurements, depth to the bottom of the magnetic source estimation and thermal modeling. Geochemistry, Geophysics, Geosystems 18(7), 2657–2672. doi: 10.1002/2016GC006755CrossRefGoogle Scholar
Engelhardt, H (2004) Ice temperature and high geothermal flux at Siple Dome, West Antarctica, from borehole measurements. Journal of Glaciology 50(169), 251–256. doi: 10.3189/172756504781830105CrossRefGoogle Scholar
Fischer, Hand 9 others (2013) Where to find 1.5 million yr old ice for the IPICS oldest-ice ice core. Climate of the Past 9(6), 2489–2505. doi: 10.5194/cp-9-2489-2013CrossRefGoogle Scholar
Fisher, AT, Mankoff, KD, Tulaczyk, SM, Tyler, SW and Foley, N (2015) High geothermal heat flux measured below the West Antarctic ice sheet. Science Advances 1(6), 1–9. doi: 10.1126/sciadv.1500093.CrossRefGoogle ScholarPubMed
Foster, TD (1978) Temperature and salinity fields under the Ross Ice Shelf. Antarctic Journal of the United States 13, 81–82.Google Scholar
Fretwell, Pand 34 others (2013) Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. Cyrosphere 7, 375–393. doi: 10.5194/tc-7-375-2013CrossRefGoogle Scholar
Gow, AJ, Ueda, HT and Garfield, DE (1968) Antarctic ice sheet: preliminary results of first core hole to bedrock. Science (New York, N.Y.) 161(3845), 1011–1013.CrossRefGoogle Scholar
Hasterok, D, Gard, M and Webb, J (2017) On the radiogenic heat production of metamorphic, igneous, and sedimentary rocks. Geoscience Frontiers 9, 1777–1794. doi: 10.1016/j.gsf.2017.10.012.CrossRefGoogle Scholar
Hondoh, T, Shoji, H, Watanabe, O, Salamatin, AN and Lipenkov, VY (2002) Depth-age and temperature prediction at Dome Fuji station, East Antarctica. Annals of Glaciology 35(1), 384–390. doi: 10.3189/172756402781817013CrossRefGoogle Scholar
Huybrechts, P (1990) A 3-D model for the Antarctic ice sheet: a sensitivity study on the glacial-interglacial contrast. Climate Dynamics 5(2), 79–92.CrossRefGoogle Scholar
Johnsen, SJ, Dahl-Jensen, D, Dansgaard, W and Gundestrup, N (1995) Greenland palaeotemperatures derived from GRIP bore hole temperature and ice core isotope profiles. Tellus B: Chemical and Physical Meteorology 47(5), 624–629. doi: 10.3402/tellusb.v47i5.16077CrossRefGoogle Scholar
Karlsson, NBand 6 others (2018) Glaciological characteristics in the Dome Fuji region and new assessment for ‘oldest ice’. Cryosphere 12(7), 2413–2424. doi: 10.5194/tc-12-2413-2018CrossRefGoogle Scholar
Koutnik, MRand 7 others (2016) Holocene accumulation and ice flow near the West Antarctic Ice Sheet Divide ice core site. Journal of Geophysical Research: Earth Surface 121(5), 907–924. doi: 10.1002/2015JF003668Google Scholar
Larour, E, Morlighem, M, Seroussi, H, Schiermeier, J and Rignot, E (2012) Ice flow sensitivity to geothermal heat flux of Pine Island Glacier, Antarctica. Journal of Geophysical Research F: Earth Surface 117(4). doi: 10.1029/2012JF002371Google Scholar
Levy, Rand 10 others (2016) Antarctic ice sheet sensitivity to atmospheric CO2 variations in the early to mid-Miocene. Proceedings of the National Academy of Sciences 113(13), 3453–3458. doi: 10.1073/pnas.1516030113CrossRefGoogle ScholarPubMed
Liefferinge, BV and Pattyn, F (2013) Using ice-flow models to evaluate potential sites of million year-old ice in Antarctica. Climate of the Past 9(5), 2335–2345. doi: 10.5194/cp-9-2335-2013CrossRefGoogle Scholar
Martos, YMand 6 others (2017) Heat flux distribution of Antarctica unveiled. Geophysical Research Letters 44(22), 11417–11426. doi: 10.1002/2017GL075609.CrossRefGoogle Scholar
Maule, CF, Purucker, ME, Olsen, N and Mosegaard, K (2005) Geophysics: heat flux anomalies in Antarctica revealed by satellite magnetic data. Science (New York, N.Y.) 309(5733), 464–467. doi: 10.1126/science.1106888CrossRefGoogle Scholar
Morgan, Vand 5 others (1997) Site information and initial results from deep ice drilling on Law Dome, Antarctica. Journal of Glaciology 43(143), 3–10.CrossRefGoogle Scholar
Morin, RHand 5 others (2010) Heat flow and hydrologic characteristics at the AND-1B borehole, ANDRILL McMurdo Ice Shelf Project, Antarctica. Geosphere 6(4), 370–378. doi: 10.1130/GES00512.1CrossRefGoogle Scholar
Motoyama, Hand 6 others (2008) Characteristics of basal ice and subglacial water at Dome Fuji, Antarctica ice sheet. Presented at the 2008 AGU Fall Meeting Abstracts, San Francisco, CA.Google Scholar
Mulvaney, Rand 8 others (2012) Recent Antarctic Peninsula warming relative to Holocene climate and ice-shelf history. Nature 489(7414), 141–144. doi: 10.1038/nature11391.CrossRefGoogle ScholarPubMed
Muto, A, Scambos, TA, Steffen, K, Slater, AG and Clow, GD (2011) Recent surface temperature trends in the interior of East Antarctica from borehole firn temperature measurements and geophysical inverse methods. Geophysical Research Letters 38(15). doi: 10.1029/2011GL048086CrossRefGoogle Scholar
Näslund, JO, Jansson, P, Fastook, JL, Johnson, J and Andersson, L (2005) Detailed spatially distributed geothermal heat-flow data for modeling of basal temperatures and meltwater production beneath the Fennoscandian ice sheet. Annals of Glaciology 40, 95–101. doi: 10.3189/172756405781813582CrossRefGoogle Scholar
Nicholls, K and Paren, J (1993) Extending the Antarctic meteorological record using ice-sheet temperature profiles. Journal of Climate 6(1), 141–150. doi: 10.1175/1520-0442(1993)006%3C0141:ETAMRU%3E2.0.CO2.0.CO;2>CrossRefGoogle Scholar
Orsi, AJ, Cornuelle, BD and Severinghaus, JP (2012) Little ice age cold interval in West Antarctica: evidence from borehole temperature at the West Antarctic Ice Sheet (WAIS) divide. Geophysical Research Letters 39. doi: 10.1029/2012GL051260CrossRefGoogle Scholar
Parrenin, Fand 10 others (2017) Is there 1.5-million-year-old ice near Dome C, Antarctica?. Cryosphere 11(6), 2427–2437. doi: 10.5194/tc-11-2427-2017CrossRefGoogle Scholar
Pattyn, F (2010) Antarctic subglacial conditions inferred from a hybrid ice sheet/ice stream model. Earth and Planetary Science Letters 295(3), 451–461. doi: 10.1016/j.epsl.2010.04.025CrossRefGoogle Scholar
Pittard, ML, Roberts, JL, Galton-Fenzi, BK and Watson, CS (2016) Sensitivity of the Lambert-Amery glacial system to geothermal heat flux. Annals of Glaciology 57(73), 56–68. doi: 10.1017/aog.2016.26CrossRefGoogle Scholar
Pollack, HN (1982) The heat flow from the continents. Annual Review of Earth and Planetary Sciences 10(1), 459–481.CrossRefGoogle Scholar
Price, PBand 9 others (2002) Temperature profile for glacial ice at the South Pole: implications for life in a nearby subglacial lake. Proceedings of the National Academy of Sciences 99(12), 7844–7847. doi: 10.1073/pnas.082238999CrossRefGoogle Scholar
Purucker, M (2013) Geothermal heat flux data set based on low resolution observations collected by the CHAMP satellite between 2000 and 2010, and produced from the MF-6 model following the technique described in Fox Maule et al. (2005). Available at: http://websrv.cs.umt.edu/isis/index.php (last access: 16 July 2018).Google Scholar
Rintoul, SR (2018) The global influence of localized dynamics in the southern ocean. Nature 558(7709), 209–218. doi: 10.1038/s41586-018-0182-3CrossRefGoogle ScholarPubMed
Risk, GF and Hochstein, MP (1974) Heat flow at arrival heights, Ross Island, Antarctica. New Zealand Journal of Geology and Geophysics 17(3), 629–644. doi: 10.1080/00288306.1973.10421586CrossRefGoogle Scholar
Ritz, C, Rommelaere, V and Dumas, C (2001) Modeling the evolution of Antarctic ice sheet over the last 420,000 years: implications for altitude changes in the Vostok region. Journal of Geophysical Research: Atmospheres 106(D23), 31943–31964. doi: 10.1029/2001JD900232CrossRefGoogle Scholar
Roberts, Jand 6 others (2013) Borehole temperatures reveal a changed energy budget at Mill Island, East Antarctica, over recent decades. Cryosphere 7(1), 263–273. doi: 10.5194/tc-7-263-2013CrossRefGoogle Scholar
Schön, J (2011) Physical properties of rocks: A workbook, volume 8. Amsterdam,Elsevier.Google Scholar
Schröder, H, Paulsen, T and Wonik, T (2011) Thermal properties of the AND-2A borehole in the southern Victoria Land Basin, McMurdo sound, Antarctica. Geosphere 7(6), 1324–1330. doi: 10.1130/GES00690.1CrossRefGoogle Scholar
Schroeder, DM, Blankenship, DD, Young, DA and Quartini, E (2014) Evidence for elevated and spatially variable geothermal flux beneath the West Antarctic Ice Sheet. Proceedings of the National Academy of Sciences 111(25), 9070–9072. doi: 10.1073/pnas.1405184111CrossRefGoogle ScholarPubMed
Shapiro, NM and Ritzwoller, MH (2004) Inferring surface heat flux distributions guided by a global seismic model: particular application to Antarctica. Earth and Planetary Science Letters 223(1), 213–224. doi: 10.1016/j.epsl.2004.04.011CrossRefGoogle Scholar
Siegert, MJ and Dowdeswell, JA (1996) Spatial variations in heat at the base of the Antarctic ice sheet from analysis of the thermal regime above subglacial lakes. Journal of Glaciology 42(142), 501–509. doi: 10.1017/S0022143000003488CrossRefGoogle Scholar
Turcotte, D and Schubert, G (2014) Geodynamics. New York,Cambridge University Press.CrossRefGoogle Scholar
Van der Veen, CJ, Leftwich, T, Von Frese, R, Csatho, B and Li, J (2007) Subglacial topography and geothermal heat flux: potential interactions with drainage of the Greenland ice sheet. Geophysical Research Letters 34(12). doi: 10.1029/2007GL030046CrossRefGoogle Scholar
Watanabe, Oand 6 others (2003) Dating of the Dome Fuji, Antarctica deep ice core. Memoirs of National Institute of Polar Research, Special issue 57, 25–37.Google Scholar
Whitehouse, PL, Gomez, N, King, MA and Wiens, DA (2019) Solid earth change and the evolution of the Antarctic ice sheet. Nature Communications 10(1), 503. doi: 10.1038/s41467-018-08068-y.CrossRefGoogle ScholarPubMed
Zagorodnov, Vand 6 others (2012) Borehole temperatures reveal details of 20th century warming at Bruce Plateau, Antarctic Peninsula. Cryosphere 6(3), 675–686. doi: 10.5194/tc-6-675-2012CrossRefGoogle Scholar
Figure 0

Fig. 1. Location of deep temperature profile boreholes and different local GHF estimations. Ice-sheet thickness data from BEPMAP 2 (Fretwell and others, 2013) and Basemap from Quantarctica and the Norwegian Polar Institute 1: Gow and others (1968); 2: Risk and Hochstein (1974); 3: Foster (1978); 4: Decker and Bucher (1982); 5: Nicholls and Paren (1993); 6: Dahl-Jensen and others (1999); 7: Hondoh and others (2002); 8: Price and others (2002); 9: Engelhardt (2004); 10: Carter and others (2009); 11: Morin and others (2010); 12: Schröder and others (2011); 13: Clow and others (2012); 14: Zagorodnov and others (2012); 15: Mulvaney and others (2012); 16 and 17: Carson and others (2014); 18: Fisher and others (2015); 19: Dmitriev and others (2016); 20: Parrenin and others (2017); 21: Begeman and others (2017); 22: Dziabek and others (2017) (26 measurements).

Figure 1

Fig. 2. GHF values from available continent-wide GHF maps for deep ice borehole locations in Antarctica (see Fig. 1).

Figure 2

Fig. 3. Ice-borehole temperature profiles for Law Dome (DSS site) (Morgan and others, 1997), Dome Fuji (Hondoh and others, 2002) and West Antarctic Ice Sheet Divide (Clow and others, 2012).

Figure 3

Table 1. Parameters used for testing the steady-state assumption

Figure 4

Table 2. Input values for the bedrock

Figure 5

Table 3. Optimised input values in the ice

Figure 6

Fig. 4. GHF estimation for DSS site from 260 iterations of our model by using a temperature profile of 1176.8 m ice equivalent depth.

Figure 7

Fig. 5. (A) GHF median values with their MAD for the DSS site using temperature profile truncated at different depths. (B) GHF estimation from Dahl-Jensen and others (1999), Martos and others (2017), Liefferinge and Pattyn (2013), Shapiro and Ritzwoller (2004).

Figure 8

Fig. 6. Comparison between the probability distribution at DSS for the entire temperature profile (LD100) and 20% of the temperature profile (LD20).

Figure 9

Fig. 7. (A) GHF median values with their MAD for Dome Fuji by using temperature profile truncated at different depths. (B) GHF estimation from Martos2017: Martos and others (2017), Liefferinge2013: Liefferinge and Pattyn (2013), Shapiro2004: Shapiro and Ritzwoller (2004), Hondoh2002: Hondoh and others (2002).

Figure 10

Fig. 8. (A) GHF median values with their MAD for WAIS divide by using temperature profile truncated at different depths. (B) GHF estimation from Martos2017: Martos and others (2017), Liefferinge2013: Liefferinge and Pattyn (2013), Shapiro2004: Shapiro and Ritzwoller (2004).

Figure 11

Table 4. Mean GHF estimations and Median Absolute Deviation (MAD) for Shapiro and Ritzwoller (2004), Liefferinge and Pattyn (2013) and Martos and others (2017)