scholarly journals The impact of parametric uncertainty and topographic error in ice-sheet modelling

2008 ◽  
Vol 54 (188) ◽  
pp. 899-919 ◽  
Author(s):  
Felix Hebeler ◽  
Ross S. Purves ◽  
Stewart S.R. Jamieson

AbstractIce-sheet models (ISMs) developed to simulate the behaviour of continental-scale ice sheets under past, present or future climate scenarios are subject to a number of uncertainties from various sources. These sources include the conceptualization of the ISM and the degree of abstraction and parameterizations of processes such as ice dynamics and mass balance. The assumption of spatially or temporally constant parameters (such as degree-day factor, atmospheric lapse rate or geothermal heat flux) is one example. Additionally, uncertainties in ISM input data such as topography or precipitation propagate to the model results. In order to assess and compare the impact of uncertainties from model parameters and climate on the GLIMMER ice-sheet model, a parametric uncertainty analysis (PUA) was conducted. Parameter variation was deduced from a suite of sensitivity tests, and accuracy information was deduced from input data and the literature. Recorded variation of modelled ice extent across the PUA runs was 65% for equilibrium ice sheets. Additionally, the susceptibility of ISM results to modelled uncertainty in input topography was assessed. Resulting variations in modelled ice extent in the range of 0.8–6.6% are comparable to that of ISM parameters such as flow enhancement, basal traction and geothermal heat flux.

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>


1996 ◽  
Vol 42 (140) ◽  
pp. 37-42 ◽  
Author(s):  
Jan T. Heine ◽  
David F. Mctigue

AbstractA finite-difference numerical model is used to simulate the temperature profile at the center of an ice sheet throughout the course of a glaciation. The ice sheet is gradually built to a thickness of 3000 m over about 10 000 years, starting on permafrost. A geothermal heat flux is applied at large depth. For an initial surface temperature of –12.5°C, our model shows that basal melting occurs 72000 years after the onset of the glaciation. The important parameters determining the basal temperatures are the initial temperature of the ice and substrate, the rate of downward advection of cold ice and, to a lesser extent, the thickness of the ice sheet. The growth history of the ice sheet does not significantly influence the time at which basal melting occurs. Our results show the possibility that the central parts of the continental ice sheets were cold-based for a significant part of their existence. Heating due to the geothermal heat flux cannot account for basal melting during most or all of a glacial cycle. These results may help to explain the existence of preserved land forms under the ice sheets.


1996 ◽  
Vol 42 (140) ◽  
pp. 37-42 ◽  
Author(s):  
Jan T. Heine ◽  
David F. Mctigue

AbstractA finite-difference numerical model is used to simulate the temperature profile at the center of an ice sheet throughout the course of a glaciation. The ice sheet is gradually built to a thickness of 3000 m over about 10 000 years, starting on permafrost. A geothermal heat flux is applied at large depth. For an initial surface temperature of –12.5°C, our model shows that basal melting occurs 72000 years after the onset of the glaciation. The important parameters determining the basal temperatures are the initial temperature of the ice and substrate, the rate of downward advection of cold ice and, to a lesser extent, the thickness of the ice sheet. The growth history of the ice sheet does not significantly influence the time at which basal melting occurs. Our results show the possibility that the central parts of the continental ice sheets were cold-based for a significant part of their existence. Heating due to the geothermal heat flux cannot account for basal melting during most or all of a glacial cycle. These results may help to explain the existence of preserved land forms under the ice sheets.


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.


2021 ◽  
Vol 14 (2) ◽  
pp. 889-905
Author(s):  
Rupert Gladstone ◽  
Benjamin Galton-Fenzi ◽  
David Gwyther ◽  
Qin Zhou ◽  
Tore Hattermann ◽  
...  

Abstract. A number of important questions concern processes at the margins of ice sheets where multiple components of the Earth system, most crucially ice sheets and oceans, interact. Such processes include thermodynamic interaction at the ice–ocean interface, the impact of meltwater on ice shelf cavity circulation, the impact of basal melting of ice shelves on grounded ice dynamics and ocean controls on iceberg calving. These include fundamentally coupled processes in which feedback mechanisms between ice and ocean play an important role. Some of these mechanisms have major implications for humanity, most notably the impact of retreating marine ice sheets on the global sea level. In order to better quantify these mechanisms using computer models, feedbacks need to be incorporated into the modelling system. To achieve this, ocean and ice dynamic models must be coupled, allowing runtime information sharing between components. We have developed a flexible coupling framework based on existing Earth system coupling technologies. The open-source Framework for Ice Sheet–Ocean Coupling (FISOC) provides a modular approach to coupling, facilitating switching between different ice dynamic and ocean components. FISOC allows fully synchronous coupling, in which both ice and ocean run on the same time step, or semi-synchronous coupling in which the ice dynamic model uses a longer time step. Multiple regridding options are available, and there are multiple methods for coupling the sub-ice-shelf cavity geometry. Thermodynamic coupling may also be activated. We present idealized simulations using FISOC with a Stokes flow ice dynamic model coupled to a regional ocean model. We demonstrate the modularity of FISOC by switching between two different regional ocean models and presenting outputs for both. We demonstrate conservation of mass and other verification steps during evolution of an idealized coupled ice–ocean system, both with and without grounding line movement.


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 66 (256) ◽  
pp. 303-312
Author(s):  
Silje Smith-Johnsen ◽  
Basile de Fleurian ◽  
Kerim H. Nisancioglu

AbstractThe spatial distribution of geothermal heat flux (GHF) under ice sheets is largely unknown. Nonetheless, it is an important boundary condition in ice-sheet models, and suggested to control part of the complex surface velocity patterns observed in some regions. Here we investigate the effect of including subglacial hydrology when modelling ice streams with elevated GHF. We use an idealised ice stream geometry and a thermomechanical ice flow model coupled to subglacial hydrology in the Ice Sheet System Model (ISSM). Our results show that the dynamic response of the ice stream to elevated GHF is greatly enhanced when including the interactive subglacial hydrology. On the other hand, the impact of GHF on ice temperature is reduced when subglacial hydrology is included. In conclusion, the sensitivity of ice stream dynamics to GHF is likely to be underestimated in studies neglecting subglacial hydrology.


2020 ◽  
Author(s):  
Rupert Gladstone ◽  
Benjamin Galton-Fenzi ◽  
David Gwyther ◽  
Qin Zhou ◽  
Tore Hattermann ◽  
...  

Abstract. A number of important questions concern processes at the margins of ice sheets where multiple components of the Earth System, most crucially ice sheets and oceans, interact. Such processes include thermodynamic interaction at the ice-ocean interface, the impact of melt water on ice shelf cavity circulation, the impact of basal melting of ice shelves on grounded ice dynamics, and ocean controls on iceberg calving. These include fundamentally coupled processes in which feedback mechanisms between ice and ocean play an important role. Some of these mechanisms have major implications for humanity, most notably the impact of retreating marine ice sheets on global sea level. In order to better quantify these mechanisms using computer models, feedbacks need to be incorporated into the modelling system. To achieve this ocean and ice dynamic models must be coupled, allowing run time information sharing between components. We have developed a flexible coupling framework based on existing Earth System coupling technologies. The open-source Framework for Ice Sheet – Ocean Coupling (FISOC) provides a modular approach to online coupling, facilitating switching between different ice dynamic and ocean components. FISOC allows fully synchronous coupling, in which both ice and ocean run on the same time-step, or semi-synchronous coupling in which the ice dynamic model uses a longer time step. Multiple regridding options are available, and multiple methods for coupling the sub ice shelf cavity geometry. Thermodynamic coupling may also be activated. We present idealised simulations using FISOC with a Stokes flow ice dynamic model coupled to a regional ocean model. We demonstrate the modularity of FISOC by switching between two different regional ocean models and presenting outputs for both. We demonstrate conservation of mass and other verification steps during evolution of an idealised coupled ice – ocean system, both with and without grounding line movement.


2021 ◽  
Author(s):  
Paul Halas ◽  
Jeremie Mouginot ◽  
Basile de Fleurian ◽  
Petra Langebroek

<div> <p>Ice losses from the Greenland Ice Sheet have been increasing in the last two decades, leading to a larger contribution to the global sea level rise. Roughly 40% of the contribution comes from ice-sheet dynamics, mainly regulated by basal sliding. The sliding component of glaciers has been observed to be strongly related to surface melting, as water can eventually reach the bed and impact the subglacial water pressure, affecting the basal sliding.  </p> </div><div> <p>The link between ice velocities and surface melt on multi-annual time scale is still not totally understood even though it is of major importance with expected increasing surface melting. Several studies showed some correlation between an increase in surface melt and a slowdown in velocities, but there is no consensus on those trends. Moreover those investigations only presented results in a limited area over Southwest Greenland.  </p> </div><div> <p>Here we present the ice motion over many land-terminating glaciers on the Greenland Ice Sheet for the period 2000 - 2020. This type of glacier is ideal for studying processes at the interface between the bed and the ice since they are exempted from interactions with the sea while still being relevant for all glaciers since they share the same basal friction laws. The velocity data was obtained using optical Landsat 7 & 8 imagery and feature-tracking algorithm. We attached importance keeping the starting date of our image pairs similar, and avoided stacking pairs starting before and after melt seasons, resulting in multiple velocity products for each year.  </p> </div><div> <p>Our results show similar velocity trends for previously studied areas with a slowdown until 2012 followed by an acceleration. This trend however does not seem to be observed on the whole ice sheet and is probably specific to this region’s climate forcing. </p> </div><div> <p>Moreover comparison between ice velocities from different parts of Greenland allows us to observe the impact of different climatic trends on ice dynamics.</p> </div>


Sign in / Sign up

Export Citation Format

Share Document