scholarly journals Inversion of Geothermal Heat Flux under the Ice Sheet of Princess Elizabeth Land, East Antarctica

2021 ◽  
Vol 13 (14) ◽  
pp. 2760
Author(s):  
Lin Li ◽  
Xueyuan Tang ◽  
Jingxue Guo ◽  
Xiangbin Cui ◽  
Enzhao Xiao ◽  
...  

Antarctic geothermal heat flux is a basic input variable for ice sheet dynamics simulation. It greatly affects the temperature and mechanical properties at the bottom of the ice sheet, influencing sliding, melting, and internal deformation. Due to the fact that the Antarctica is covered by a thick ice sheet, direct measurements of heat flux are very limited. This study was carried out to estimate the regional heat flux in the Antarctic continent through geophysical inversion. Princess Elizabeth Land, East Antarctica is one of the areas in which we have a weak understanding of geothermal heat flux. Through the latest airborne geomagnetic data, we inverted the Curie depth, obtaining the heat flux of bedrock based on the one-dimensional steady-state heat conduction equation. The results indicated that the Curie depth of the Princess Elizabeth Land is shallower than previously estimated, and the heat flux is consequently higher. Thus, the contribution of subglacial heat flux to the melting at the bottom of the ice sheet is likely greater than previously expected in this region. It further provides research clues for the formation of the developed subglacial water system in Princess Elizabeth Land.

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

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


2020 ◽  
Author(s):  
Francesco Salvini ◽  
Paola Cianfarra ◽  
Giovanni Capponi ◽  
Laura Crispini ◽  
Laura Federico ◽  
...  

<p>Estimation of subglacial Geothermal Heat Flux (GHF) is of paramount importance to better understand the dynamics  of cryosphere and ice flow of the East Antarctica Ice Sheet (EAIS). Unfortunately, the GHF of East Antarctica is still poorly known and constrained, and direct measurements are still challenging. The EIAS is underlain by major subglacial mountain ranges and basins resulting from distinct geodynamic domains. These include Northern Victoria Land-Ross Sea, the Transantarctic Mountains, the Wilkes Subglacial Basin, the Gamburtsev Subglacial Mountains, the East Antarctic System and a major transpressional fault zone in between (e.g. Cianfarra & Maggi, 2017), which hosts clusters of subglacial lakes. The distribution of sedimentary basins and tectonic structures may affect the GHF in that it exhibits strong regional variations as testified by the presence of subglacial lakes at bedrock topographic elevation/depth with a range exceeding 1500 m, from deep subglacial basins to the flanking highlands. In the framework of the G-IDEA (Geo Ice Dynamics of East Antarctica) project, heat flow from the basement is quantified in key areas of East Antarctica between 60°E and 180°E, by an innovative application of the HCA (Hybrid Cellular Automata) method: the description of stationary conditions of the temperature field is used to replicate the observed distribution of wet vs dry ice-rock contacts in an ice-flowing environment. Evaluation of the geothermal flux is performed in key areas based on the numeric modeling of the ice-rock interaction, which can replicate the spatial distribution of wet contacts and subglacial lakes and is related to local dynamics of the ice sheet and its interaction with the atmosphere. The model takes into account the spatial distribution of the Curie temperature depth as derived from literature. The heat flux is estimated by modeling the stationary state of the ice-rock system with the HCA numerical method, and by its discretization into a large number of cells. Each cell is characterized by physical parameters such as density, enthalpy, thermal capacity and conductivity. By their interaction it is possible to compute their temperature evolution and to replicate the heat diffusion by conduction and convection (the ice movement) in the interfaces ice-rock and ice-atmosphere. The final resolution of the model is about 100 m. The presence of possible anomalous heath flow in the bedrock are identified by a stochastic approach that allow the estimation of the error in the computed heath flow values.</p>


Author(s):  
T. M Kyrke-Smith ◽  
R. F Katz ◽  
A. C Fowler

Antarctic ice streams are associated with pressurized subglacial meltwater but the role this water plays in the dynamics of the streams is not known. To address this, we present a model of subglacial water flow below ice sheets, and particularly below ice streams. The base-level flow is fed by subglacial melting and is presumed to take the form of a rough-bedded film, in which the ice is supported by larger clasts, but there is a millimetric water film which submerges the smaller particles. A model for the film is given by two coupled partial differential equations, representing mass conservation of water and ice closure. We assume that there is no sediment transport and solve for water film depth and effective pressure. This is coupled to a vertically integrated, higher order model for ice-sheet dynamics. If there is a sufficiently small amount of meltwater produced (e.g. if ice flux is low), the distributed film and ice sheet are stable, whereas for larger amounts of melt the ice–water system can become unstable, and ice streams form spontaneously as a consequence. We show that this can be explained in terms of a multi-valued sliding law, which arises from a simplified, one-dimensional analysis of the coupled model.


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>


2021 ◽  
Author(s):  
Agnes Wansing ◽  
Jörg Ebbing ◽  
Mareen Lösing ◽  
Sergei Lebedev ◽  
Nicolas Celli ◽  
...  

<p>The lithospheric structure of Greenland is still poorly known due to its thick ice sheet, the sparseness of seismological stations, and the limitation of geological outcrops near coastal areas. As only a few geothermal measurements are available for Greenland, one must rely on geophysical models. Such models of Moho and LAB depths and sub-ice geothermal heat-flow vary largely.</p><p>Our approach is to model the lithospheric architecture by geophysical-petrological modelling with LitMod3D. The model is built to reproduce gravity observations, the observed elevation with isostasy assumptions and the velocities from a tomography model. Furthermore, we adjust the thermal parameters and the temperature structure of the model to agree with different geothermal heat flow models. We use three different heat flow models, one from machine learning, one from a spectral analysis of magnetic data and another one which is compiled from a similarity study with tomography data.</p><p>For the latter, a new shear wave tomography model of Greenland is used. Vs-depth profiles from Greenland are compared with velocity profiles from the US Array, where a statistical link between Vs profiles and surface heat flow has been established. A similarity function determines the most similar areas in the U.S. and assigns the mean heat-flow from these areas to the corresponding area in Greenland.</p><p>The geothermal heat flow models will be further used to discuss the influence on ice sheet dynamics by comparison to friction heat and viscous heat dissipation from surface meltwater.</p>


2021 ◽  
Author(s):  
Haoran Kang ◽  
Liyun Zhao ◽  
Michael Wolovick ◽  
John C. Moore

Abstract. Basal thermal conditions play an important role in ice sheet dynamics, and they are sensitive to geothermal heat flux (GHF). Here we estimate the basal thermal conditions, including basal temperature, basal melt rate, and friction heat underneath the Lambert-Amery glacier system in east Antarctica, using a combination of a forward model and an inversion from a 3D ice flow model. We assess the sensitivity and uncertainty of basal thermal conditions using six different GHFs. We evaluate the modelled results using all available observed subglacial lakes. There are very large differences in modelled spatial pattern of temperate basal conditions using the different GHFs. The two most-recent GHF fields inverted from aerial geomagnetic observations have higher values of GHF in the region, produce a larger warm-based area, and match the observed subglacial lakes better than the other GHFs. The fast flowing glacier region has a lower modelled basal friction coefficient, faster basal velocity, with higher basal frictional heating in the range of 50–2000 mW m−2 than the base under slower flowing glaciated areas. The modelled basal melt rate reaches ten to hundreds of mm per year locally in Lambert, Lepekhin and Kronshtadtskiy glaciers feeding the Amery ice shelf, and ranges from 0–5 mm yr−1 on the temperate base of the vast inland region.


2007 ◽  
Vol 53 (181) ◽  
pp. 289-297 ◽  
Author(s):  
Sergey V. Popov ◽  
Valery N. Masolov

AbstractDuring the summer field seasons of 1987–91, studies of central East Antarctica by airborne radio-echo sounding commenced. This scientific work continued in the 1990s in the Vostok Subglacial Lake area and along the traverse route from Mirny, and led to the discovery of 16 new subglacial water cavities in the areas of Domes Fuji and Argus and the Prince Charles Mountains. Twenty-nine subglacial water cavities were revealed in the area near Vostok, along with a feature we believe to be a subglacial river. Two subglacial lakes were discovered along the Mirny–Vostok traverse route. These are located 50 km north of Komsomolskaya station and under Pionerskaya station. We find high geothermal heat flux in the vicinity of the largest of the subglacial lakes, and suggest this may be due to their location over deep faults where additional mantle heat is available.


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.


2000 ◽  
Vol 46 (153) ◽  
pp. 197-205 ◽  
Author(s):  
Christoph Mayer ◽  
Martin J. Siegert

AbstractA numerical model of the ice-sheet/ice-shelf transition was used to investigate ice-sheet dynamics across the large subglacial lake beneath Vostok station, central East Antarctica. European Remote-sensing Satellite (ERS-1) altimetry of the ice surface and 60 MHz radio-echo sounding (RES) of the ice-sheet base and internal ice-sheet layering were used to develop a conceptual flowline across the ice sheet, which the model used as input. The model calculates horizontal and vertical velocities and stresses, from which particle flow paths can be obtained, and the ice-sheet temperature distribution. An inverse approach to modelling was adopted, where particle flow paths were forced to match those identified from internal RES layering. Results show that ice dynamics across the inflow grounding line are similar to an ice-sheet/ice-shelf transition. Model particle flow paths match internal RES layering when ice is (a) taken away from the ice base across the first 2 km of the flowline over the lake and (b) added to the base across the remainder of the lake. We contend that the process causing this transfer of ice is likely to be melting of ice and freezing of water at the ice–water interface. Other explanations, such as enhanced rates of accumulation over the grounding line, or three-dimensional convergent/divergent flow of ice are inconsistent with available measurements. Such melting and refreezing would be responsible for circulation and mixing of at least the surface layers of the lake water. Our model suggests that several tens of metres of refrozen “basal ice” would accrete from lake water to the ice sheet before the ice regrounds.


Sign in / Sign up

Export Citation Format

Share Document