scholarly journals Exceptionally High Geothermal Heat Flux Needed to Sustain the Northeast Greenland Ice Stream

2019 ◽  
Author(s):  
Silje Smith-Johnsen ◽  
Basile de Fleurian ◽  
Nicole Schlegel ◽  
Helene Seroussi ◽  
Kerim Nisanciolgu

Abstract. The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area, and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet, but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting hotspots with various locations and intensities. We find that a minimum geothermal heat flux value of 970 mW/m2 located close to EastGRIP is required locally to reproduce the observed NEGIS velocities, consistent with previous estimates. By including high geothermal heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.

2020 ◽  
Vol 14 (3) ◽  
pp. 841-854 ◽  
Author(s):  
Silje Smith-Johnsen ◽  
Basile de Fleurian ◽  
Nicole Schlegel ◽  
Helene Seroussi ◽  
Kerim Nisancioglu

Abstract. The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting heat flux anomalies with various locations and intensities. In our model experiment, a minimum heat flux value of 970 mW m−2 located close to the East Greenland Ice-core Project (EGRIP) is required locally to reproduce the observed NEGIS velocities, giving basal melt rates consistent with previous estimates. The value cannot be attributed to geothermal heat flux alone and we suggest hydrothermal circulation as a potential explanation for the high local heat flux. By including high heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.


2019 ◽  
Vol 65 (254) ◽  
pp. 1023-1034 ◽  
Author(s):  
Soroush Rezvanbehbahani ◽  
Leigh A. Stearns ◽  
C. J. van der Veen ◽  
Gordon K. A. Oswald ◽  
Ralf Greve

AbstractThe spatial distribution of basal water critically impacts the evolution of ice sheets. Current estimates of basal water distribution beneath the Greenland Ice Sheet (GrIS) contain large uncertainties due to poorly constrained boundary conditions, primarily from geothermal heat flux (GHF). The existing GHF models often contradict each other and implementing them in numerical ice-sheet models cannot reproduce the measured temperatures at ice core locations. Here we utilize two datasets of radar-detected basal water in Greenland to constrain the GHF at regions with a thawed bed. Using the three-dimensional ice-sheet model SICOPOLIS, we iteratively adjust the GHF to find the minimum GHF required to reach the bed to the pressure melting point, GHFpmp, at locations of radar-detected basal water. We identify parts of the central-east, south and northwest Greenland with significantly high GHFpmp. Conversely, we find that the majority of low-elevation regions of west Greenland and parts of northeast have very low GHFpmp. We compare the estimated constraints with the available GHF models for Greenland and show that GHF models often do not honor the estimated constraints. Our results highlight the need for community effort to reconcile the discrepancies between radar data, GHF models, and ice core information.


2021 ◽  
Author(s):  
Paul D. Bons ◽  
Tamara de Riese ◽  
Steven Franke ◽  
Maria-Gema Llorens ◽  
Till Sachau ◽  
...  

<p>The prominent North East Greenland Ice Stream (NEGIS) is an exceptionally large ice stream in the Greenland Ice sheet. It is over 500 km long, originates almost at the central ice divide, and contributes significantly to overall ice drainage from the Greenland Ice sheet. Surface velocities in the inland part of the ice stream are several times higher inside NEGIS than in the adjacent ice sheet. Modelling NEGIS is still a challenge as it remains unclear what actually causes and controls the ice stream.</p><p>An elevated geothermal heat flux is one of the factors that are being considered to trigger or drive the fast flow inside NEGIS. Unfortunately, the geothermal heat flux below NEGIS and its upstream area is poorly constrained and estimates vary from close to the global average for continental crust (ca. 60 mW/m<sup>2</sup>) to values as high as almost 1000 mW/m<sup>2</sup>. The latter would cause about 10 cm/yr of melting at the base of the ice sheet.</p><p>We present a brief survey of global geothermal heat flux data, especially from known hotspots, such as Iceland and Yellowstone. Heat fluxes in these areas that are known to be among the hottest on Earth rarely, if ever, exceed 300 mW/m<sup>2</sup>. A plume hotspot or its trail can therefore not cause heat fluxes at the high end of the suggested range. Other potential factors, such as hydrothermal fluid flow and radiogenic heat, also cannot raise the heat flux significantly. We conclude that the heat flux at NEGIS is very unlikely to exceed 100-150 mW/m<sup>2</sup>, and future modelling studies on NEGIS should thus be mindful of implementing realistic geothermal heat flux values. If NEGIS is not the result of an exceptionally high heat flux, we are left with the exciting challenge to find the true trigger of this fascinating structure.</p>


2005 ◽  
Vol 42 ◽  
pp. 424-432 ◽  
Author(s):  
Ralf Greve

AbstractThe thermomechanical, three-dimensional ice-sheet model SICOPOLIS is applied to the Greenland ice sheet. Simulations over two glacial–interglacial cycles are carried out, driven by a climatic forcing interpolated between present conditions and Last Glacial Maximum anomalies. Based on the global heat-flow representation by Pollack and others (1993), we attempt to constrain the spatial pattern of the geothermal heat flux by comparing simulation results to direct measurements of basal temperatures at the GRIP, NorthGRIP, Camp Century and Dye 3 ice-core locations. The obtained heat-flux map shows an increasing trend from west to east, a high-heat-flux anomaly around NorthGRIP with values up to 135 mWm–2 and a low-heat-flux anomaly around Dye 3 with values down to 20 mW m–2. Validation is provided by the generally good fit between observed and measured ice thicknesses. Residual discrepancies are most likely due to deficiencies of the input precipitation rate and further variability of the geothermal heat flux not captured here.


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

<p>The origin of Greenland’s largest ice stream – the Northeast Greenland Ice Stream (NEGIS) – is located far inland of the Greenland Ice Sheet. High surface flow velocities in the center of NEGIS are attributed to the lubrication of the ice sheet base facilitated by basal melt water. In order to derive basal melt rates at the EastGRIP drill site (~2668 m thick ice), we performed in-situ measurements with an autonomous phase-sensitive radar (ApRES; Brennan et al., 2014; Nicholls et al., 2015) in two consecutive years. The precise processing method (Stewart et al., 2019 and Vankova et al., 2020) detects englacial and basal vertical displacements, but it is limited due to noisy data in the lower half of the ice column. Thus, we made assumptions for the vertical strain in the lower half and adapted simulation results (Rückamp et al., 2020). We found melt rates ranging from 0.16 to 0.22 m/a, which is extremely large for inland ice. However, our results are only slightly above melt rates from previous studies (Fahnestock et al., 2001 and MacGregor et al., 2016) which found melt rates of 0.10 m/a and more through airborne radar measurements (evaluated using radiostratigraphy methods) in the vicinity of EastGRIP. Melt rates of >0.16 m/a require a heat flux into the ice of >1.55 W/m<sup>2</sup> which is not exclusively the geothermal heat flux, as also the subglacial hydrological system may supply a significant heat flux into the ice.</p>


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>


2011 ◽  
Vol 52 (59) ◽  
pp. 43-50 ◽  
Author(s):  
Douglas J. Brinkerhoff ◽  
Toby W. Meierbachtol ◽  
Jesse V. Johnson ◽  
Joel T. Harper

AbstractA full-stress, thermomechanically coupled, numerical model is used to explore the interaction between basal thermal conditions and motion of a terrestrially terminating section of the west Greenland ice sheet. The model domain is a two-dimensional flowline profile extending from the ice divide to the margin. We use data-assimilation techniques based on the adjoint model in order to optimize the basal traction field, minimizing the difference between modeled and observed surface velocities. We monitor the sensitivity of the frozen/melted boundary (FMB) to changes in prescribed geothermal heat flux and sliding speed by applying perturbations to each of these parameters. The FMB shows sensitivity to the prescribed geothermal heat flux below an upper threshold where a maximum portion of the bed is already melted. The position of the FMB is insensitive to perturbations applied to the basal traction field. This insensitivity is due to the short distances over which longitudinal stresses act in an ice sheet.


2020 ◽  
Author(s):  
Ilka Weikusat ◽  
David Wallis ◽  
Steven Franke ◽  
Nicolas Stoll ◽  
Julien Westhoff ◽  
...  

<p>Drilling an ice core through an ice sheet (typically 2000 to 3000 m thick) is a technical challenge that nonetheless generates valuable and unique information on palaeo-climate and ice dynamics. As technically the drilling cannot be done in one run, the core has to be fractured approximately every 3 m to retrieve core sections from the bore hole. This fracture process is initiated by breaking the core with core-catchers which also clamp the engaged core in the drill head while the whole drill is then pulled up with the winch motor.</p><p> </p><p>This standard procedure is known to become difficult and requires extremely high pulling forces (Wilhelms et al. 2007), in the very deep part of the drill procedure, close to the bedrock of the ice sheet, especially when the ice material becomes warm (approximately -2°C) due to the geothermal heat released from the bedrock. Recently, during the EastGRIP (East Greenland Ice coring Project) drilling we observed a similar issue with breaking off cored sections only with extremely high pulling forces, but started from approximately 1800 m of depth, where the temperature is still very cold (approximately -20°C). This has not been observed at other ice drilling sites. As dependencies of fracture behaviour on crystal orientation and grain size are known (Schulson & Duval 2009) for ice, we thus examined the microstructure in the ice samples close to and at the core breaks.</p><p> </p><p>First preliminary results suggest that these so far unexperienced difficulties are due to the profoundly different c-axes orientation distribution (CPO) in the EastGRIP ice core. In contrast to other deep ice cores which have been drilled on ice domes or ice divides, EastGRIP is located in an ice stream. This location means that the deformation geometry (kinematics) is completely different, resulting in a different CPO (girdle pattern instead of single maximum pattern). Evidence regarding additional grain-size dependence will hopefully help to refine the fracturing procedure, which is possible due to a rather strong grain size layering observed in natural ice formed by snow precipitation.</p><p> </p><p>---------------------</p><p>Wilhelms, F.; Sheldon, S. G.; Hamann, I. & Kipfstuhl, S. Implications for and findings from deep ice core drillings - An example: The ultimate tensile strength of ice at high strain rates. Physics and Chemistry of Ice (The proceedings of the International Conference on the Physics and Chemistry of Ice held at Bremerhaven, Germany on 23-28 July 2006), <strong>2007</strong>, 635-639</p><p>Schulson, E. M. & Duval, P. Creep and Fracture of Ice. Cambridge University Press, <strong>2009</strong>, 401</p>


1995 ◽  
Vol 21 ◽  
pp. 8-12 ◽  
Author(s):  
R. Greve ◽  
K. Hutter

Computations over 50 000 years into steady state with Greve’s polythermal ice-sheet model and its numerical code are performed for the Greenland ice sheet with today’s climatological input (surface temperature and accumulation function) and three values of the geothermal heat flux: (42, 54.6, 29.4) mW m−2. It is shown that through the thermomechanical coupling the geometry as well as the thermal regime, in particular that close to the bed, respond surprisingly strongly to the basal thermal heat input. The most sensitive variable is the basal temperature field, but the maximum height of the summit also varies by more than ±100m. Furthermore, some intercomparison of the model outputs with the real ice sheet is carried out, showing that the model provides reasonable results for the ice-sheet geometry as well as for the englacial temperatures.


Sign in / Sign up

Export Citation Format

Share Document