scholarly journals Surface velocity and ice thickness of the Müller ice cap, Axel Heiberg Island

2021 ◽  
Author(s):  
Ann-Sofie Priergaard Zinck

Muller ice cap is situated on Axel Heiberg Island in Arctic Canada. It is characterised by a mountanious region separating the ice cap in the east from the outlet glaciers in the west. Research has taken place on the outlet glaciers of the ice cap since 1959, but only limited research has been conducted on the main ice cap, and no full depth ice cores have ever been drilled. The interesting location of the ice cap facing the Arctic Ocean, the chance of finding ice dating back to the Innuitian Ice Sheet, and the fact that no full depth ice cores have been drilled, makes it an obvious place to do so. In order to achieve a long and undisturbed chronology of the ice core, one needs to find a location for the drilling site, where there is a great ice thickness, low surface velocity and little melt. In this project the aim is to make surface velocity maps of the ice cap and estimate the ice thickness to be able to come up with suggestions of possible drill site areas.Surface velocities are calculated using feature tracking of optical satellite images from the Landsat satellites in the period of 2004-2019. A median velocity map of all Landsat 8 velocity maps is used as validation in modelling the ice thickness and in the investigation of possible drill site areas.To estimate the ice thickness various methods are used and are being compared to the ice thickness measured by Operation IceBridge. The first method is an iterative inverse method where the ice sheet model PISM works as a forward model. The model is found to work rather well on the ice cap, with a root mean squared error (RMS) of 138.9 m, but overestimates the ice thickness on the outlet glaciers. The second model uses a simple inversion of the shallow ice approximation. It overestimates the ice thickness in areas with low surface slope, but has a RMS of 131.4 m on the ice cap. The third and and fourth models uses Monte Carlo sampling methods of the shallow ice approximation without and with sliding, respectively. The latter uses an initial ice thickness guess, and the modelled ice thickness was proofed not to differ from that initial guess at all. The RMSs on the icecap of the two models were found to be 132.1 m and 129.9 m, respectively. Finally, the fifth model uses the PISM setup but with an initial geometry defined by the SIA inversion. The RMS on the ice cap is found to be 135.4 m.Based on the median Landsat 8 surface velocity map, the modelled ice thicknesses and the surface elevation from the Arctic Digital Elevation Model, a map of suggested drill site areas is made. The site which fulfilled the criteria the most is located at 526629 m easting and 8866463 m northing in UTM zone 15N. In this site the surface velocity is 1.2 m yr−1, the surface elevation is 1804 m and the modelled ice thicknesses varies from 535-579 m. The melt in this area is estimated to be less than 20 melt days per year based on the backscatter from Sentinel-1.

2013 ◽  
Vol 7 (2) ◽  
pp. 499-510 ◽  
Author(s):  
J. L. Bamber ◽  
J. A. Griggs ◽  
R. T. W. L. Hurkmans ◽  
J. A. Dowdeswell ◽  
S. P. Gogineni ◽  
...  

Abstract. We present a new bed elevation dataset for Greenland derived from a combination of multiple airborne ice thickness surveys undertaken between the 1970s and 2012. Around 420 000 line kilometres of airborne data were used, with roughly 70% of this having been collected since the year 2000, when the last comprehensive compilation was undertaken. The airborne data were combined with satellite-derived elevations for non-glaciated terrain to produce a consistent bed digital elevation model (DEM) over the entire island including across the glaciated–ice free boundary. The DEM was extended to the continental margin with the aid of bathymetric data, primarily from a compilation for the Arctic. Ice thickness was determined where an ice shelf exists from a combination of surface elevation and radar soundings. The across-track spacing between flight lines warranted interpolation at 1 km postings for significant sectors of the ice sheet. Grids of ice surface elevation, error estimates for the DEM, ice thickness and data sampling density were also produced alongside a mask of land/ocean/grounded ice/floating ice. Errors in bed elevation range from a minimum of ±10 m to about ±300 m, as a function of distance from an observation and local topographic variability. A comparison with the compilation published in 2001 highlights the improvement in resolution afforded by the new datasets, particularly along the ice sheet margin, where ice velocity is highest and changes in ice dynamics most marked. We estimate that the volume of ice included in our land-ice mask would raise mean sea level by 7.36 m, excluding any solid earth effects that would take place during ice sheet decay.


1996 ◽  
Vol 23 ◽  
pp. 364-373 ◽  
Author(s):  
Jonathan L. Bamber ◽  
Philippe Huybrechts

This paper presents improved geometric boundary conditions (surface elevation and ice thickness) required as inputs to calculations of the surface-velocity field for the Antarctic ice sheet. A comparison of the two-dimensional horizontal velocity field obtained on the basis of conservation of mass (balance velocity) with the diagnostic velocity field calculated with an ice-sheet model (dynamic velocity) may yield information on shortcomings in the way the ice-sheet model describes the ice flow. Here, the surface-elevation grid is described in detail, as it has been generated specifically for such a study and represents a new standard in accuracy and resolution for calculating surface slopes. The digital-elevation model was generated on a 10 km grid size from over 20 000 000 height estimates obtained from eight 35 d repeat cycles of ERS-1 radar-altimeter data. For surface slopes less than 0.4°, the accuracy is better than 1.5 m. In areas of high surface slope (coastal and mountainous regions), the altimeter measurements have been supplemented with data taken from the Antarctic Digital Database. South of 81.5°, data from the SPRI folio map have been used. The ice-thickness grid was produced from a combination of a redigitization of the SPRI folio and the original radio-echo-sounding flight lines. For areas of grounded ice, the elevation of the bed was estimated from surface elevation and ice thickness. Significant differences (in excess of 25% of ice thickness) were obtained between an earlier digitization of the folio bed-elevation map and the data set derived here. Furthermore, a new value of 25.6 × 106 km3 was obtained for the total volume of the ice sheet and ice shelves, which is a reduction of 12% compared with the original estimate derived during the compilation of the SPRI folio. These differences will have an important influence on the results obtained by numerical ice-sheet models.


1996 ◽  
Vol 23 ◽  
pp. 364-373 ◽  
Author(s):  
Jonathan L. Bamber ◽  
Philippe Huybrechts

This paper presents improved geometric boundary conditions (surface elevation and ice thickness) required as inputs to calculations of the surface-velocity field for the Antarctic ice sheet. A comparison of the two-dimensional horizontal velocity field obtained on the basis of conservation of mass (balance velocity) with the diagnostic velocity field calculated with an ice-sheet model (dynamic velocity) may yield information on shortcomings in the way the ice-sheet model describes the ice flow. Here, the surface-elevation grid is described in detail, as it has been generated specifically for such a study and represents a new standard in accuracy and resolution for calculating surface slopes. The digital-elevation model was generated on a 10 km grid size from over 20 000 000 height estimates obtained from eight 35 d repeat cycles of ERS-1 radar-altimeter data. For surface slopes less than 0.4°, the accuracy is better than 1.5 m. In areas of high surface slope (coastal and mountainous regions), the altimeter measurements have been supplemented with data taken from the Antarctic Digital Database. South of 81.5°, data from the SPRI folio map have been used. The ice-thickness grid was produced from a combination of a redigitization of the SPRI folio and the original radio-echo-sounding flight lines. For areas of grounded ice, the elevation of the bed was estimated from surface elevation and ice thickness. Significant differences (in excess of 25% of ice thickness) were obtained between an earlier digitization of the folio bed-elevation map and the data set derived here. Furthermore, a new value of 25.6 × 106 km3 was obtained for the total volume of the ice sheet and ice shelves, which is a reduction of 12% compared with the original estimate derived during the compilation of the SPRI folio. These differences will have an important influence on the results obtained by numerical ice-sheet models.


2021 ◽  
pp. 1-13
Author(s):  
Iben Koldtoft ◽  
Aslak Grinsted ◽  
Bo M. Vinther ◽  
Christine S. Hvidberg

Abstract To assess the amount of ice volume stored in glaciers or ice caps, a method to estimate ice thickness distribution is required for glaciers where no direct observations are available. In this study, we use an existing inverse method to estimate the bedrock topography and ice thickness of the Renland Ice Cap, East Greenland, using satellite-based observations of the surface topography. The inverse approach involves a procedure in which an ice dynamical model is used to build-up an ice cap in steady state with climate forcing from a regional climate model, and the bedrock is iteratively adjusted until the modelled and observed surface topography match. We validate our model results against information from airborne radar data and satellite observed surface velocity, and we find that the inferred ice thickness and thereby the stored total volume of the ice cap is sensitive to the assumed ice softness and basal slipperiness. The best basal model parameters for the Renland Ice Cap are determined and the best estimated total ice volume of 384 km3 is found. The Renland Ice Cap is particularly interesting because of its location at a high elevation plateau and hence assumed low sensitivity to climate change.


2021 ◽  
Author(s):  
Maria Hoerhold ◽  
Thomas Münch ◽  
Stefanie Weißbach ◽  
Sepp Kipfstuhl ◽  
Bo Vinther ◽  
...  

<p>Climate variability of the Arctic region has been investigated by means of temperature reconstructions based on proxies from various climate archives around the Arctic, compiled over the last 2000a in the so called Arctic2k record. However, the representativeness of the Arctic2k reconstruction for central Greenland remains unclear, since only a few ice cores have been included in the reconstruction, and observations from the Greenland Ice Sheet (GIC) report ambiguous warming trends for the end of the 20th and the beginning of the 21st century which are not displayed by Arctic2k. Today, the GIC experiences periods with temperatures close to or above the freezing point at high elevations, area-wide melting and mass loss. In order to assess the recent warming as signature of global climate change, records of past climate changes with appropriate temporal and spatial coverage can serve as a benchmark for naturally driven climate variability. Instrumental records for Greenland are short and geographically sparse, and existing temperature reconstructions from single ice cores are noisy, leading to an inconclusive assessment of the recent warming for Greenland.</p><p>Here, we provide a Greenland firn-core stack covering the time span of the last millennium until the first decade of the 21<sup>st </sup>century in unprecedented quality by re-drilling as well as analyzing 16 existing firn core sites. We find a strong decadal to bi-decadal natural variability in the record, and, while the record exhibits several warming events with trends that show a similar amplitude as the recent one, we find that the recent absolute values of stable oxygen isotope composition are unprecedented for the last 1000 years.</p><p> </p><p>Comparing our Greenland record with the Arctic 2k temperature reconstruction shows that the correlation between the two records changes throughout the last millennium. While in the periods of 1200-1300 and 1400-1650 CE the records correlate positively, between 1300 and 1400 and 1650-1700 CE shorter periods with negative correlation are found. Since then the correlation is characterized by alternation between positive and zero correlation, with a drop towards negative values at the end of the 20<sup>th</sup> century. Including re-analysis data, we hypothesize that the climate on top of the GIC was decoupled from the surrounding Arctic for the last decades, leading to the observed mismatch in observations of warming trends.</p><p>We suggest that the recently observed Greenland temperatures are a superposition of a strong natural variability with an anthropogenic long-term trend. Our findings illustrate that global warming has reached the interior of the Greenland ice sheet, which will have implications for its surface mass balance and Greenland’s future contribution to sea level rise.</p><p>Our record complements the Arctic 2k record to a profound view on the Arctic climate variability, where regional compilations may not be representative for specific areas.</p>


2021 ◽  
Author(s):  
Ann-Sofie Priergaard Zinck ◽  
Aslak Grinsted

<p><span>The ice thickness of the Müller Ice Cap, Arctic Canada, is estimated using regression parameters obtained from an inversion of the shallow ice approximation by the use of a single Operation IceBridge flight line in combination with the glacier outline, surface slope, and elevation. The model is compared with an iterative inverse method of estimating the bedrock topography using PISM as a forward model. In both models the surface elevation is given by the Arctic Digital Elevation Model. The root mean squared errors of the ice thickness on the ice cap is 131 m and 139 m for the shallow ice inversion and the PISM model, respectively. Including the outlet glaciers increases the root mean squared errors to 136 m and 396 m, respectively. </span></p><p><span>The simplicity of the shallow ice inversion model, combined with the good results and the fact that only remote sensing data is needed, means that there is a possibility of applying this model in a global glacier thickness estimate by using the Randolph Glacier Inventory. Most global glacier estimates only provide the volume and not the ice thickness of the glaciers. Hence, global ice thickness models is of great importance in quantifying the potential contribution of sea level rise from the glaciers and ice caps around the globe. </span></p>


1996 ◽  
Vol 23 ◽  
pp. 226-236 ◽  
Author(s):  
Philippe Huybrechts

A high-resolution, three-dimensional thermomechanical ice-sheet model, which includes isostasy, the possibility of ice-sheet expansion on the continental shelf and refined climatic parameterizations, was used to investigate the basal thermal regime of the Greenland ice sheet. The thermodynamic calculations take into account the usual terms of heat flow within the ice, a thermally active bedrock layer and all of the effects associated with changes in ice thickness and flow pattern. Basal temperature conditions are documented with respect to glacial–interracial shifts in climatic boundary conditions, both in steady state as during simulations over the last two glacial cycles using the GRIP δ180 record. It is found that the basal temperature field shows a large sensitivity in steady-state experiments but that, during a glacial cycle, basal temperature variations are strongly damped, in particular in central areas. A comparison has been made with measured data from deep ice cores and the implications are discussed.


2013 ◽  
Vol 59 (213) ◽  
pp. 9-20 ◽  
Author(s):  
Reinhard Drews ◽  
Carlos Martín ◽  
Daniel Steinhage ◽  
Olaf Eisen

AbstractWe present a comprehensive approach (including field data, remote sensing and an anisotropic ice-flow model) to characterize Halvfarryggen ice dome in coastal Dronning Maud Land, Antarctica. This is a potential drill site for the International Partnerships in Ice Core Sciences, which has identified the need for ice cores covering atmospheric conditions during the last few millennia. We derive the surface topography, the ice stratigraphy from radar data, and accumulation rates which vary from 400 to 1670 kg m−2 a−1 due to preferred wind directions and changing surface slope. The stratigraphy shows anticlines and synclines beneath the divides. We transfer Dansgaard–Johnsen age–depth scales from the flanks along isochrones to the divide in the upper 20–50% of the ice thickness and show that they compare well with the results of a full-Stokes, anisotropic ice-flow model which predicts (1) 11 ka BP ice at 90% of the ice thickness, (2) a temporally stable divide for at least 2700–4500 years, (3) basal temperatures below the melting point (−12°C to −5°C) and (4) a highly developed crystal orientation fabric (COF). We suggest drilling into the apices of the deep anticlines, providing a good compromise between record length and temporal resolution and also facilitating studies of the interplay of anisotropic COF and ice flow.


1993 ◽  
Vol 5 (1) ◽  
pp. 63-75 ◽  
Author(s):  
M. O. Jeffries ◽  
W. F. Weeks

The internal structure of ice cores from western Ross Sea pack ice floes showed considerable diversity. Snow-ice formation made a small, but significant contribution to ice growth. Frazil ice was common and its growth clearly occurred during both the pancake cycle and deformation events. Congelation ice was also common, in both its crystallographically aligned and non-aligned varieties. Platelet ice was found in only one core next to the Drygalski Ice Tongue, an observation adding to the increasing evidence that this unusual ice type occurs primarily in coastal pack ice near ice tongues and ice shelves. The diverse internal structure of the floes indicates that sea ice development in the Ross Sea is as complex as that in the Weddell Sea and more complex than in the Arctic. The mean ice thickness at the ice core sites varied between 0.71 m and 1.52 m. The thinnest ice generally occurred in the outer pack ice zone. Regardless of latitude, the ice thickness data are further evidence that Antarctic sea ice is thinner than Arctic sea ice.


1989 ◽  
Vol 12 ◽  
pp. 46-50 ◽  
Author(s):  
D.M. Etheridge

The internal dynamics of the Law Dome ice cap have been investigated by measuring the deformation of three bore holes located on an approximate flow line. Bore holes BHC1 (300 m deep) and BHC2 (344 m) were drilled in the coastal area to within several metres of bedrock and BHQ (418 m) was drilled about half-way towards the dome centre to about 50% of the ice thickness. Detailed measurements of orientation (inclination and azimuth), diameter, and temperature were taken through each bore hole over a 1 year span for BHC1 and BHC2 and a 10 year span for BHQ. The orientation data were reduced to obtain ∂u/∂z, a measure of the shear strain-rate. Changes in the depth of features located by bore-hole diameter measurements were used to obtain vertical velocity profiles. Other measurements discussed are temperatures, oxygen isotopes, crystal structure, surface velocities, and surface and bedrock topography.At the coastal sites, the ∂u/∂z profiles show two maxima in the lower third of the ice sheet. Flow due to the measured deformation accounts for about 55% of the surface velocity, the remainder being due to deformation and sliding in the basal zone. The vertical velocity profiles show mostly firn compression. The deeper ∂u/∂z maximum occurs in ice from the Wisconsin period which appears to deform more rapidly than the Holocene ice immediately above. The upper ∂u/∂z maximum may be related to the stress history of the ice, which can also explain the presence of significant shear strain and crystal-fabric development at only half the ice thickness at the BHQ site.


Sign in / Sign up

Export Citation Format

Share Document