BLOG

We have a new open-access study about geothermal heat flow beneath the Greenland and Antarctic ice sheets in the Journal of Geophysical Research: Earth Surface. 

Presently, there’s a lot of uncertainty about the magnitude and pattern of geothermal heat flow beneath both ice sheets. That’s because it has only been sampled at a handful of widely spaced deep ice cores (Figure 1). While the average value of geothermal heat flow is relatively small, getting it right is essential for ice-flow models. If you run a computer simulation of an ice sheet with a severe over- or under-estimation of the geothermal heat flow, you can easily end up generating an ice sheet that is either too warm or too cold. Ice flow is very sensitive to temperature – particularly near the bed – so geothermal heat flow is a critical variable for simulating the form and flow of Earth’s ice sheets.

Our study took a fresh look at changes in geothermal heat flow across space, but not those due to the subtle variations in Earth’s crust and mantle properties over tens of km. Instead, we examined the effect of the ice-sheet bed’s topographic relief on geothermal heat flow to generate the first comprehensive snapshot of changes in geothermal heat flow at scales of hundreds of meters due to that relief. It’s been known for over a century that geothermal heat flow is greater in valleys and smaller on ridges. Basically, if the heat escaping Earth’s interior is looking for the quickest way to radiate into the atmosphere, a deeply incised valley provides the fastest exit. This effect is readily observable from the fact that geotherms – surfaces of constant temperature – are packed more closely together beneath valleys, indicating a stronger temperature gradient there (and hence heat flow) in comparison to ridges.

We created a simple statistical model to estimate this topographic influence on geothermal heat flow. This model essentially uses a digital elevation model of the bedrock topography to assess local topographic relief and then converts this local relief into a fractional correction for geothermal heat flow. It produces a positive correction – an increase – for valleys, and vice versa for ridges. Our approach is admittedly simple and empirical – a literal “first-order” approximation – but it seems to reliably reproduce the topographic variability in geothermal heat flow in all the settings for which we could find previous studies. So, we applied this statistical model to digital elevation models for Greenland and Antarctica. This revealed much more detail in a geothermal heat flow map than we are used to seeing.

Across both Greenland and Antarctica, we see patterns of increased geothermal heat flow within deeply incised glacier valleys and decreased geothermal heat flow along ridges and mountains. In many regions, most notably the Antarctic Peninsula (Figure 2) and Central East Greenland (Figure 3), we find that local topography routinely modifies regional geothermal heat flow by more than ~50%.

In Greenland, we estimate that there are ~100 outlet glaciers that are both sufficiently narrow and deeply incised to more than double local geothermal heat flow relative to that of the regional average value. The model also suggests that – especially deep within the interior of the Greenland Ice Sheet – local geothermal heat flow may be sufficiently suppressed along prominent subglacial ridges to cause subglacial water to refreeze. (At least, in ice-sheet areas where the ice-bed interface is near the freezing point.) 

The topographic correction for geothermal heat flow that we model is only as a good as the topographic relief that we derive from subglacial digital elevation models. Generally, in areas where subglacial topography is best resolved, the effect of the topography on the local geothermal heat flow is greater than uncertainty in the underlying regional geothermal heat flow model (Figure 4). There are still large swaths of the ice sheets where subglacial topography remains poorly resolved. In these areas, our model is not tremendously useful; the existing uncertainties between regional geothermal heat flow models are still larger than any local topographic correction that we can estimate.

The topographic corrections for geothermal heat flow in Greenland and Antarctica that we have calculated are now available as dimensionless fields in NetCDF format, with grids that are the same resolution as BedMachine, for each region via the PROMICE data portal (www.promice.dk). This means that they can be anonymously downloaded and applied to any regional geothermal heat flow model of the user’s choice. We hope that these topographic corrections for geothermal heat flow will be adopted into ice-flow models to improve both present-day ice-sheet simulations, as well as our understanding of the role of geothermal heat flow in the feedback between ice flow and topography on geologic timescales.

It has certainly been a long and winding road to this publication – the original draft of this article was first submitted in June 2019 – and we are grateful for Noah Finnegan (University of California Santa Cruz) and Olga Sergienko (Princeton University) for serving as editors to four very helpful peer-reviewers. Interdisciplinary projects can clearly provide a bumpier ride than staying in your own lane, but – in this instance – the journey seems to have taken us to a very different view of geothermal heat flow in Greenland and Antarctica.

Development of this data product was funded by the award “HOTROD: Prototype for Rapid Sampling of Ice-Sheet Basal Temperatures” provided by the Experiment Programme of the Villum Foundation. Improved understanding of the spatial variability in subglacial geothermal heat flow helps us optimize drill site selection and analyze ice temperature measurements. This data product was also supported by the Danish Ministry for Climate, Energy and Supplies through the Programme for Monitoring of the Greenland Ice Sheet (PROMICE).

Colgan, W., J. MacGregor, K. Mankoff, R. Haagenson, H. Rajaram, Y. Martos, M. Morlighem, M. Fahnestock and K. Kjeldsen. 2021. Topographic Correction of Geothermal Heat Flux in Greenland and Antarctica. Journal of Geophysical Research. 125: e2020JF005598. doi:10.1029/2020JF005598.

We have a new open-access study about Greenland Ice Sheet mass balance – or health – in the current issue of Geophysical Research Letters. In this study, we present a new 28-year record of ice-sheet mass balance. This record is relatively unique for two reasons.

Firstly, because of its length. The most recent ice-sheet mass balance inter-comparison exercise (IMBIE2) clearly highlighted how the availability of ice-sheet mass balance estimates has changed through time. During the GRACE satellite gravimetry era (2003-2017), there are usually more than twenty independent estimates of annual Greenland Ice Sheet mass balance. Prior to 2003, however, there are just two independent estimates. Our new 1992-2020 mass balance record will therefore provide especially welcomed additional insight on ice-sheet mass balance during the 1990s.

Secondly, because of its consistency. This new mass balance record has been constructed by merging radar altimetry measurements from four ESA satellites (ERS-1/2, ENVISAT, CryoSat-2 and Sentinel-3A/B) over nearly three decades into one consistent framework. While all four of these satellites use the same type of Ku-band radar altimeter, to date, their measurements have usually been analyzed independently of each other. This time, however, we use machine learning to merge the elevation changes measured by these similar-but-different satellites into a common mass balance signal through space and time. This makes our new record the only satellite altimetry record that spans the entire IMBIE period.

When we compare our new radar altimetry record of mass balance to two existing input-output records of mass balance, we find good agreement in the capture of Greenland’s high and low mass balance years. These other two multi-decade records are derived from the input-output method, in which estimated iceberg calving into the oceans is differenced from estimated surface mass balance (or net snow accumulation) over the ice sheet. While the input-output method often has limited spatial (and temporal) resolution, our radar altimetry derived record can resolve spatial variability in mass balance across the ice sheet every month since 1992.

While our new long-term record provides a new overview of the health of the Greenland ice sheet, it can also be helpful to understand the processes that influence ice-sheet health. For example, we see a sharp increase in mass balance between 2016 and 2017. When we look at this event in detail, we can attribute it to unusually high snowfall in fall 2016, especially in East Greenland, and unusually little surface melting in summer 2017, throughout the ice-sheet ablation area. We estimate that the 2017 hydrological year was likely the first year during the 21^st Century during which the ice sheet was actually in a state of true “mass balance” – or equilibrium – as opposed to mass loss.

The development of this new dataset was primarily funded by the European Space Agency (ESA), with a little help from the Programme for Monitoring of the Greenland Ice Sheet (www.promice.dk). Our multi-satellite Ku-band altimetry mass balance record is now available as tabulated data – both for the ice sheet, as well as the eight major ice-sheet drainage sectors – at https://doi.org/10.11583/DTU.13353062. Within the next two years, the ongoing Sentinel-3A/B satellite missions are clearly poised to extend Greenland’s radar altimetry record to three decades. This will allow us to start assessing ice-sheet health using the statistics of a 30-year climatology record. This keeps us excited at the prospect of updating this record in the near future. Stay tuned!

Simonsen, S., V. Barletta, W. Colgan and L. Sørensen. 2021. Greenland Ice Sheet mass balance (1992-2020) from calibrated radar altimetry. Geophysical Research Letters. L61865. doi:10.1029/2020GL091216.