scholarly journals Using ice-flow models to evaluate potential sites of million year-old ice in Antarctica

2013 ◽  
Vol 9 (3) ◽  
pp. 2859-2887 ◽  
Author(s):  
B. Van Liefferinge ◽  
F. Pattyn
Keyword(s):  
Ice Sheet ◽  
Geothermal Heat ◽  
Ice Flow ◽  
Flow Models ◽  
Drill Site ◽  
Slow Moving ◽  
Basal Melting ◽  
The Antarctic ◽  
Suitable Potential

Abstract. Finding suitable potential sites for an undisturbed record of million-year old ice in Antarctica requires a slow-moving ice sheet (preferably an ice divide) and basal conditions that are not disturbed by large topographic variations. Furthermore, ice should be thick and cold basal conditions should prevail, since basal melting would destroy the bottom layers. However, thick ice (needed to resolve the signal at sufficient high resolution) increases basal temperatures, which is a conflicting condition in view of finding a suitable drill site. In addition, slow moving areas in the center of ice sheets are also low-accumulation areas, and low accumulation reduces potential cooling of the ice through vertical advection. While boundary conditions such as ice thickness and accumulation rates are relatively well constraint, the major uncertainty in determining basal conditions resides in the geothermal heat flow (GHF) underneath the ice sheet. We explore uncertainties in existing GHF datasets and their effect on basal temperatures of the Antarctic ice sheet and propose an updated method based on Pattyn (2010) to improve existing GHF datasets in agreement with known basal temperatures and their gradients to reduce this uncertainty. Both complementary methods lead to a better comprehension of basal temperature sensitivity and a characterization of potential ice coring sites within these uncertainties.

scholarly journals Using ice-flow models to evaluate potential sites of million year-old ice in Antarctica

2013 ◽  
Vol 9 (5) ◽  
pp. 2335-2345 ◽  
Author(s):  
B. Van Liefferinge ◽  
F. Pattyn
Keyword(s):  
Ice Sheet ◽  
Thermal Conditions ◽  
Ice Thickness ◽  
Data Sets ◽  
Geothermal Heat ◽  
Flow Models ◽  
Slow Moving ◽  
Basal Melting ◽  
The Antarctic

Abstract. Finding suitable potential sites for an undisturbed record of million-year old ice in Antarctica requires slow-moving ice (preferably an ice divide) and basal conditions that are not disturbed by large topographic variations. Furthermore, ice should be thick and cold basal conditions should prevail, since basal melting would destroy the bottom layers. However, thick ice (needed to resolve the signal at sufficient high resolution) increases basal temperatures, which is a conflicting condition for finding a suitable drill site. In addition, slow moving areas in the center of ice sheets are also low-accumulation areas, and low accumulation reduces potential cooling of the ice through vertical advection. While boundary conditions such as ice thickness and accumulation rates are relatively well constrained, the major uncertainty in determining basal thermal conditions resides in the geothermal heat flow (GHF) underneath the ice sheet. We explore uncertainties in existing GHF data sets and their effect on basal temperatures of the Antarctic Ice Sheet, and propose an updated method based on Pattyn (2010) to improve existing GHF data sets in agreement with known basal temperatures and their gradients to reduce this uncertainty. Both complementary methods lead to a better comprehension of basal temperature sensitivity and a characterization of potential ice coring sites within these uncertainties. The combination of both modeling approaches show that the most likely oldest ice sites are situated near the divide areas (close to existing deep drilling sites, but in areas of smaller ice thickness) and across the Gamburtsev Subglacial Mountains.

scholarly journals Where to find 1.5 million yr old ice for the IPICS "Oldest-Ice" ice core

2013 ◽  
Vol 9 (6) ◽  
pp. 2489-2505 ◽  
Author(s):  
H. Fischer ◽  
J. Severinghaus ◽  
E. Brook ◽  
E. Wolff ◽  
M. Albert ◽  
...  
Keyword(s):  
Heat Flux ◽  
Ice Core ◽  
Ice Sheet ◽  
Ice Thickness ◽  
Geothermal Heat ◽  
Ice Flow ◽  
Antarctic Ice Sheet ◽  
Quaternary Climate ◽  
Drill Site ◽  
Heat Flow Model

Abstract. The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful pre-site selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

scholarly journals Indication of high basal melting at EastGRIP drill site on the Northeast Greenland Ice Stream

2021 ◽  
Author(s):  
Ole Zeising ◽  
Angelika Humbert
Keyword(s):  
Heat Flux ◽  
Heat Budget ◽  
Ice Sheet ◽  
Water System ◽  
Hydrological Regime ◽  
Ice Streams ◽  
Geothermal Heat ◽  
Ice Stream ◽  
Drill Site ◽  
Basal Melting

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.

scholarly journals A full-Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution

2009 ◽  
Vol 3 (1) ◽  
pp. 1-31 ◽  
Author(s):  
H. Seddik ◽  
R. Greve ◽  
T. Zwinger ◽  
L. Placidi
Keyword(s):  
Flow Model ◽  
Stokes Equations ◽  
Melting Rate ◽  
Induced Anisotropy ◽  
Geothermal Heat ◽  
Ice Flow ◽  
Climate Conditions ◽  
Anisotropic Flow ◽  
Drill Site ◽  
Basal Melting

Abstract. A three-dimensional, thermo-mechanically coupled ice flow model with induced aniso-tropy has been applied to a ~200×200 km domain around the Dome Fuji drill site, Antarctica. The model ("Elmer/Ice") is based on the open-source multi-physics package Elmer (http://www.csc.fi/elmer/) and solves the full-Stokes equations. Flow-induced anisotropy in ice is accounted for by an implementation of the Continuum-mechanical, Anisotropic Flow model, based on an anisotropic Flow Enhancement factor ("CAFFE model"). Steady-state simulations for present-day climate conditions are conducted. The main findings are: (i) the flow regime at Dome Fuji is a complex superposition of vertical compression, horizontal extension and bed-parallel shear; (ii) for a geothermal heat flux of 60 mW m−2 the basal temperature at Dome Fuji reaches the pressure melting point and the basal melting rate is ~1 mm a−1; (iii) the fabric shows a weak single maximum at Dome Fuji, which increases the age of the ice compared to an isotropic scenario; (iv) as a consequence of spatially variable basal melting conditions, and contrary to intuition, the basal age is smaller where the ice is thicker and larger where the ice is thinner. The latter result is of great relevance for the consideration of a future drill site in the area.

scholarly journals A full Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution

2011 ◽  
Vol 5 (2) ◽  
pp. 495-508 ◽  
Author(s):  
H. Seddik ◽  
R. Greve ◽  
T. Zwinger ◽  
L. Placidi
Keyword(s):  
Flow Model ◽  
Stokes Equations ◽  
Melting Rate ◽  
Induced Anisotropy ◽  
Geothermal Heat ◽  
Ice Flow ◽  
Climate Conditions ◽  
Anisotropic Flow ◽  
Drill Site ◽  
Basal Melting

Abstract. A three-dimensional, thermo-mechanically coupled ice flow model with induced anisotropy has been applied to a ~200 × 200 km domain around the Dome Fuji drill site, Antarctica. The model ("Elmer/Ice") is based on the open-source multi-physics package Elmer (http://www.csc.fi/elmer/) and solves the full Stokes equations. Flow-induced anisotropy in ice is accounted for by an implementation of the Continuum-mechanical, Anisotropic Flow model, based on an anisotropic Flow Enhancement factor ("CAFFE model"). Steady-state simulations for present-day climate conditions are conducted. The main findings are: (i) the flow regime at Dome Fuji is a complex superposition of vertical compression, horizontal extension and bed-parallel shear; (ii) for an assumed geothermal heat flux of 60 mW m−2 the basal temperature at Dome Fuji reaches the pressure melting point and the basal melting rate is ~0.35 mm a−1; (iii) in agreement with observational data, the fabric shows a strong single maximum at Dome Fuji, and the age of the ice is decreased compared to an isotropic scenario; (iv) as a consequence of spatially variable basal melting conditions, the basal age tends to be smaller where the ice is thicker and larger where the ice is thinner. The latter result is of great relevance for the consideration of a future drill site in the area.

Implications for the lithospheric structure of Greenland by applying different heat flow models

2021 ◽  
Author(s):  
Agnes Wansing ◽  
Jörg Ebbing ◽  
Mareen Lösing ◽  
Sergei Lebedev ◽  
Nicolas Celli ◽  
...  
Keyword(s):  
Heat Flow ◽  
Heat Dissipation ◽  
Ice Sheet ◽  
Magnetic Data ◽  
Similarity Function ◽  
Lithospheric Structure ◽  
Temperature Structure ◽  
Geothermal Heat ◽  
Flow Models ◽  
Ice Sheet Dynamics

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

scholarly journals Melting and freezing under Antarctic ice shelves from a combination of ice-sheet modelling and observations

2017 ◽  
Vol 63 (240) ◽  
pp. 731-744 ◽  
Author(s):  
JORGE BERNALES ◽  
IRINA ROGOZHINA ◽  
MAIK THOMAS
Keyword(s):  
Mass Balance ◽  
Ice Sheet ◽  
Model Parameters ◽  
Ice Shelf ◽  
Ice Shelves ◽  
Continental Scale ◽  
Model Set ◽  
Set Up ◽  
Basal Melting ◽  
The Antarctic

ABSTRACTIce-shelf basal melting is the largest contributor to the negative mass balance of the Antarctic ice sheet. However, current implementations of ice/ocean interactions in ice-sheet models disagree with the distribution of sub-shelf melt and freezing rates revealed by recent observational studies. Here we present a novel combination of a continental-scale ice flow model and a calibration technique to derive the spatial distribution of basal melting and freezing rates for the whole Antarctic ice-shelf system. The modelled ice-sheet equilibrium state is evaluated against topographic and velocity observations. Our high-resolution (10-km spacing) simulation predicts an equilibrium ice-shelf basal mass balance of −1648.7 Gt a−1 that increases to −1917.0 Gt a−1 when the observed ice-shelf thinning rates are taken into account. Our estimates reproduce the complexity of the basal mass balance of Antarctic ice shelves, providing a reference for parameterisations of sub-shelf ocean/ice interactions in continental ice-sheet models. We perform a sensitivity analysis to assess the effects of variations in the model set-up, showing that the retrieved estimates of basal melting and freezing rates are largely insensitive to changes in the internal model parameters, but respond strongly to a reduction of model resolution and the uncertainty in the input datasets.

scholarly journals Are longitudinal ice-surface structures on the Antarctic Ice Sheet indicators of long-term ice-flow configuration?

2014 ◽  
Vol 2 (2) ◽  
pp. 911-933 ◽  
Author(s):  
N. F. Glasser ◽  
S. J. A. Jennings ◽  
M. J. Hambrey ◽  
B. Hubbard
Keyword(s):  
Flow Line ◽  
Ice Sheet ◽  
Surface Structures ◽  
Surface Velocity ◽  
Ice Flow ◽  
Antarctic Ice Sheet ◽  
Ice Surface ◽  
Ice Stream ◽  
The Antarctic

Abstract. Continent-wide mapping of longitudinal ice-surface structures on the Antarctic Ice Sheet reveals that they originate in the interior of the ice sheet and are arranged in arborescent networks fed by multiple tributaries. Longitudinal ice-surface structures can be traced continuously down-ice for distances of up to 1200 km. They are co-located with fast-flowing glaciers and ice streams that are dominated by basal sliding rates above tens of m yr-1 and are strongly guided by subglacial topography. Longitudinal ice-surface structures dominate regions of converging flow, where ice flow is subject to non-coaxial strain and simple shear. Associating these structures with the AIS' surface velocity field reveals (i) ice residence times of ~ 2500 to 18 500 years, and (ii) undeformed flow-line sets for all major flow units analysed except the Kamb Ice Stream and the Institute and Möller Ice Stream areas. Although it is unclear how long it takes for these features to form and decay, we infer that the major ice-flow and ice-velocity configuration of the ice sheet may have remained largely unchanged for several thousand years, and possibly even since the end of the last glacial cycle. This conclusion has implications for our understanding of the long-term landscape evolution of Antarctica, including large-scale patterns of glacial erosion and deposition.

scholarly journals Enthalpy benchmark experiments for numerical ice sheet models

2015 ◽  
Vol 9 (1) ◽  
pp. 217-228 ◽  
Author(s):  
T. Kleiner ◽  
M. Rückamp ◽  
J. H. Bondzio ◽  
A. Humbert
Keyword(s):  
Analytical Solutions ◽  
Ice Sheet ◽  
Jump Conditions ◽  
Cold Side ◽  
Ice Flow ◽  
Flow Models ◽  
Transient Simulations ◽  
Rate Calculation ◽  
Transition Surface ◽  
Simulation Results

Abstract. We present benchmark experiments to test the implementation of enthalpy and the corresponding boundary conditions in numerical ice sheet models. Since we impose several assumptions on the experiment design, analytical solutions can be formulated for the proposed numerical experiments. The first experiment tests the functionality of the boundary condition scheme and the basal melt rate calculation during transient simulations. The second experiment addresses the steady-state enthalpy profile and the resulting position of the cold–temperate transition surface (CTS). For both experiments we assume ice flow in a parallel-sided slab decoupled from the thermal regime. We compare simulation results achieved by three different ice flow-models with these analytical solutions. The models agree well to the analytical solutions, if the change in conductivity between cold and temperate ice is properly considered in the model. In particular, the enthalpy gradient on the cold side of the CTS goes to zero in the limit of vanishing temperate-ice conductivity, as required from the physical jump conditions at the CTS.

Sign in / Sign up

Export Citation Format

Share Document