scholarly journals Feedback between ice dynamics and bedrock deformation with 3D viscosity in Antarctica

Author(s):  
Wouter van der Wal ◽  
Caroline van Calcar ◽  
Bas de Boer ◽  
Bas Blank

<p>Over glacial-interglacial cycles, the evolution of an ice sheet is influenced by Glacial isostatic adjustment (GIA) via two negative feedback loops. Firstly, vertical bedrock deformation due to a changing ice load alters ice-sheet surface elevation. For example, an increasing ice load leads to a lower bedrock elevation that lowers ice-sheet surface elevation. This will increase surface melting of the ice sheet, following an increase of atmospheric temperature at lower elevations. Secondly, bedrock deformation will change the height of the grounding line of the ice sheet. For example, a lowering bedrock height following ice-sheet advance increases the melt due to ocean water that in turn leads to a retreat of the grounding line and a slow-down of ice-sheet advance.      <br>               GIA is mainly determined by the viscosity of the interior of the solid Earth which is radially and laterally varying. Underneath the Antarctic ice sheet, there are relatively low viscosities in West Antarctica and higher viscosities in East Antarctica, in turn affecting the response time of the above mentioned feedbacks. However, most ice-dynamical models do not consider the lateral variations of the viscosity in the GIA feedback loops when simulating the evolution of the Antarctic ice sheet. The method developed by Gomez et al. (2018) includes the feedback between GIA and ice-sheet evolution and alternates between simulations of the two models where each simulation covers the full time period. We presents a different method to couple ANICE, a 3-D ice-sheet model, to a 3-D GIA finite element model. In this method the model computations alternates between the ice-sheet and GIA model until convergence of the result occurs at each timestep. We simulate the evolution of the Antarctic ice sheet from 120 000 years ago to the present. The results of the coupled simulation will be discussed and compared to results of the uncoupled ice-sheet model (using an ELRA GIA model) and the method developed by Gomez et al. (2018).</p>

2020 ◽  
Author(s):  
Stephen Chuter ◽  
Jonathan Rougier ◽  
Geoffrey Dawson ◽  
Jonathan Bamber

<p>Long-term continuous monitoring of Antarctic Ice Sheet mass balance is imperative to better understand its multi-decadal response to changes in climate and ocean forcing. Additionally, more accurate knowledge of contemporaneous mass balance is key for improved parameterisations in ice sheet models. The Antarctic Peninsula has undergone rapid changes in mass balance and ice dynamics over the last two decades, with satellite observations showing the presence of grounding line retreat and increases in ice sheet velocity. This is particularly the case after the collapse of the Larsen A and B ice shelves in 1995 and 2002, and more recently the glaciers draining the southern Antarctic Peninsula. As a result, this region provides analogues for future ice sheet response to ice shelf collapse in other regions of Antarctica. </p><p>Despite the region’s importance to understanding ice sheet dynamics, it is challenging to accurately assess mass balance due its geometry and mountainous topography. Conventional pulse-limited altimetry suffers from poor coverage and data loss over steep mountainous terrain, particularly before the launch of CryoSat-2 in 2010. In the case of gravimetry, the geometry of the region means the coarse spatial resolution of the GRACE mission (~300 km) cannot resolve small spatial scale glacier changes (particularly over northern Antarctic Peninsula) and suffers from signal leakage into the ocean. For the mass budget approach, the challenge of accurately modelling surface mass balance over the region’s mountainous topography coupled with the sparsity of ice thickness observations at the grounding line for many sectors can result in large uncertainties. As a result, it can be difficult to reconcile the results from different conventional approaches in this region. </p><p>To resolve this, we have developed and optimised the BHM framework used previously over the Antarctic Ice Sheet to specifically investigate the Antarctic Peninsula. This enables each latent process driving ice sheet mass change to be resolved at a higher spatial resolution compared to previous implementations across Antarctica as a whole. The new regional solution also incorporates more recent and higher resolution observations including: CryoSat-2 swath altimetry, stereo-image DEM differencing and NASA Operation Ice Bridge laser altimetry elevation rates. This is the first time such a range of observations of varying spatio-temporal resolutions will be combined into one assessment for the region. We will present results from the regionally optimised model from 2003 until present, including basin-scale mass trends and changes in spatial latent processes at an annual resolution. Additionally, we will discuss future opportunities, such as extending the record from this approach into the next decade and further understanding of the GIA response in this region. </p>


2018 ◽  
Vol 12 (5) ◽  
pp. 1767-1778 ◽  
Author(s):  
Fifi Ibrahime Adodo ◽  
Frédérique Remy ◽  
Ghislain Picard

Abstract. Spaceborne radar altimeters are a valuable tool for observing the Antarctic Ice Sheet. The radar wave interaction with the snow provides information on both the surface and the subsurface of the snowpack due to its dependence on the snow properties. However, the penetration of the radar wave within the snowpack also induces a negative bias on the estimated surface elevation. Empirical corrections of this space- and time-varying bias are usually based on the backscattering coefficient variability. We investigate the spatial and seasonal variations of the backscattering coefficient at the S (3.2 GHz ∼ 9.4 cm), Ku (13.6 GHz ∼ 2.3 cm) and Ka (37 GHz ∼ 0.8 cm) bands. We identified that the backscattering coefficient at Ku band reaches a maximum in winter in part of the continent (Region 1) and in the summer in the remaining (Region 2), while the evolution at other frequencies is relatively uniform over the whole continent. To explain this contrasting behavior between frequencies and between regions, we studied the sensitivity of the backscattering coefficient at three frequencies to several parameters (surface snow density, snow temperature and snow grain size) using an electromagnetic model. The results show that the seasonal cycle of the backscattering coefficient at Ka frequency is dominated by the volume echo and is mainly driven by snow temperature evolution everywhere. In contrast, at S band, the cycle is dominated by the surface echo. At Ku band, the seasonal cycle is dominated by the volume echo in Region 1 and by the surface echo in Region 2. This investigation provides new information on the seasonal dynamics of the Antarctic Ice Sheet surface and provides new clues to build more accurate corrections of the radar altimeter surface elevation signal in the future.


2020 ◽  
Author(s):  
Mariel Dirscherl ◽  
Andreas Dietz ◽  
Celia Baumhoer ◽  
Christof Kneisel ◽  
Claudia Kuenzer

<p>Antarctica stores ~91 % of the global ice mass making it the biggest potential contributor to global sea-level-rise. With increased surface air temperatures during austral summer as well as in consequence of global climate change, the ice sheet is subject to surface melting resulting in the formation of supraglacial lakes in local surface depressions. Supraglacial meltwater features may impact Antarctic ice dynamics and mass balance through three main processes. First of all, it may cause enhanced ice thinning thus a potentially negative Antarctic Surface Mass Balance (SMB). Second, the temporary injection of meltwater to the glacier bed may cause transient ice speed accelerations and increased ice discharge. The last mechanism involves a process called hydrofracturing i.e. meltwater-induced ice shelf collapse caused by the downward propagation of surface meltwater into crevasses or fractures, as observed along large coastal sections of the northern Antarctic Peninsula. Despite the known impact of supraglacial meltwater features on ice dynamics and mass balance, the Antarctic surface hydrological network remains largely understudied with an automated method for supraglacial lake and stream detection still missing. Spaceborne remote sensing and data of the Sentinel missions in particular provide an excellent basis for the monitoring of the Antarctic surface hydrological network at unprecedented spatial and temporal coverage.</p><p>In this study, we employ state-of-the-art machine learning for automated supraglacial lake and stream mapping on basis of optical Sentinel-2 satellite data. With more detail, we use a total of 72 Sentinel-2 acquisitions distributed across the Antarctic Ice Sheet together with topographic information to train and test the selected machine learning algorithm. In general, our machine learning workflow is designed to discriminate between surface water, ice/snow, rock and shadow being further supported by several automated post-processing steps. In order to ensure the algorithm’s transferability in space and time, the acquisitions used for training the machine learning model are chosen to cover the full circle of the 2019 melt season and the data selected for testing the algorithm span the 2017 and 2018 melt seasons. Supraglacial lake predictions are presented for several regions of interest on the East and West Antarctic Ice Sheet as well as along the Antarctic Peninsula and are validated against randomly sampled points in the underlying Sentinel-2 RGB images. To highlight the performance of our model, we specifically focus on the example of the Amery Ice Shelf in East Antarctica, where we applied our algorithm on Sentinel-2 data in order to present the temporal evolution of maximum lake extent during three consecutive melt seasons (2017, 2018 and 2019).</p>


2018 ◽  
Vol 12 (1) ◽  
pp. 49-70 ◽  
Author(s):  
Werner M. J. Lazeroms ◽  
Adrian Jenkins ◽  
G. Hilmar Gudmundsson ◽  
Roderik S. W. van de Wal

Abstract. Basal melting below ice shelves is a major factor in mass loss from the Antarctic Ice Sheet, which can contribute significantly to possible future sea-level rise. Therefore, it is important to have an adequate description of the basal melt rates for use in ice-dynamical models. Most current ice models use rather simple parametrizations based on the local balance of heat between ice and ocean. In this work, however, we use a recently derived parametrization of the melt rates based on a buoyant meltwater plume travelling upward beneath an ice shelf. This plume parametrization combines a non-linear ocean temperature sensitivity with an inherent geometry dependence, which is mainly described by the grounding-line depth and the local slope of the ice-shelf base. For the first time, this type of parametrization is evaluated on a two-dimensional grid covering the entire Antarctic continent. In order to apply the essentially one-dimensional parametrization to realistic ice-shelf geometries, we present an algorithm that determines effective values for the grounding-line depth and basal slope in any point beneath an ice shelf. Furthermore, since detailed knowledge of temperatures and circulation patterns in the ice-shelf cavities is sparse or absent, we construct an effective ocean temperature field from observational data with the purpose of matching (area-averaged) melt rates from the model with observed present-day melt rates. Our results qualitatively replicate large-scale observed features in basal melt rates around Antarctica, not only in terms of average values, but also in terms of the spatial pattern, with high melt rates typically occurring near the grounding line. The plume parametrization and the effective temperature field presented here are therefore promising tools for future simulations of the Antarctic Ice Sheet requiring a more realistic oceanic forcing.


Author(s):  
Eric Rignot

The concept that the Antarctic ice sheet changes with eternal slowness has been challenged by recent observations from satellites. Pronounced regional warming in the Antarctic Peninsula triggered ice shelf collapse, which led to a 10-fold increase in glacier flow and rapid ice sheet retreat. This chain of events illustrated the vulnerability of ice shelves to climate warming and their buffering role on the mass balance of Antarctica. In West Antarctica, the Pine Island Bay sector is draining far more ice into the ocean than is stored upstream from snow accumulation. This sector could raise sea level by 1 m and trigger widespread retreat of ice in West Antarctica. Pine Island Glacier accelerated 38% since 1975, and most of the speed up took place over the last decade. Its neighbour Thwaites Glacier is widening up and may double its width when its weakened eastern ice shelf breaks up. Widespread acceleration in this sector may be caused by glacier ungrounding from ice shelf melting by an ocean that has recently warmed by 0.3 °C. In contrast, glaciers buffered from oceanic change by large ice shelves have only small contributions to sea level. In East Antarctica, many glaciers are close to a state of mass balance, but sectors grounded well below sea level, such as Cook Ice Shelf, Ninnis/Mertz, Frost and Totten glaciers, are thinning and losing mass. Hence, East Antarctica is not immune to changes.


Nature ◽  
2009 ◽  
Vol 460 (7256) ◽  
pp. 766-766 ◽  
Author(s):  
Eric J. Steig ◽  
David P. Schneider ◽  
Scott D. Rutherford ◽  
Michael E. Mann ◽  
Josefino C. Comiso ◽  
...  

2021 ◽  
Author(s):  
Yijing Lin ◽  
Yan Liu ◽  
Zhitong Yu ◽  
Xiao Cheng ◽  
Qiang Shen ◽  
...  

Abstract. The input-output method (IOM) is one of the most popular methods of estimating the ice sheet mass balance (MB), with a significant advantage in presenting the dynamics response of ice to climate change. Assessing the uncertainties of the MB estimation using the IOM is crucial to gaining a clear understanding of the Antarctic ice-sheet mass budget. Here, we introduce a framework for assessing the uncertainties in the MB estimation due to the methodological differences in the IOM, the impact of the parameterization and scale effect on the modeled surface mass balance (SMB, input), and the impact of the uncertainties of ice thickness, ice velocity, and grounding line data on ice discharge (D, output). For the assessment of the D’s uncertainty, we present D at a fine scale. Compared with the goal of determining the Antarctic MB within an uncertainty of 15 Gt yr−1, we found that the different strategies employed in the methods cause considerable uncertainties in the annual MB estimation. The uncertainty of the RACMO2.3 SMB caused by its parameterization can reach 20.4 Gt yr−1, while that due to the scale effect is up to 216.7 Gt yr−1. The observation precisions of the MEaSUREs InSAR-based velocity (1–17 m yr−1), the airborne radio-echo sounder thickness (±100 m), and the MEaSUREs InSAR-based grounding line (±100 m) contribute uncertainties of 17.1 Gt yr−1, 10.5 ± 2.7 Gt yr−1 and 8.0~27.8 Gt yr−1 to the D, respectively. However, the D’s uncertainty due to the remarkable ice thickness data gap, which is represented by the thickness difference between the BEDMAP2 and the BedMachine reaches 101.7 Gt yr−1, which indicates its dominant cause of the future D’s uncertainty. In addition, the interannual variability of D caused by the annual changes in the ice velocity and ice thickness are considerable compared with the target uncertainty of 15 Gt yr−1, which cannot be ignored in annual MB estimations.


2018 ◽  
Author(s):  
Qiang Shen ◽  
Hansheng Wang ◽  
C. K. Shum ◽  
Liming Jiang ◽  
Hou Tse Hsu ◽  
...  

Abstract. Ice velocity constitutes a key parameter for estimating ice-sheet discharge rates and is crucial for improving coupled models of the Antarctic ice sheet to accurately predict its future fate and contribution to sea-level change. Here, we present a new Antarctic ice velocity map at a 100-m grid spacing inferred from Landsat 8 imagery data collected from December 2013 through March 2016 and robustly processed using the feature tracking method. These maps were assembled from over 73,000 displacement vector scenes inferred from over 32,800 optical images. Our maps cover nearly all the ice shelves, landfast ice, ice streams, and most of the ice sheet. The maps have an estimated uncertainty of less than 10 m yr-1 based on robust internal and external validations. These datasets will allow for a comprehensive continent-wide investigation of ice dynamics and mass balance combined with the existing and future ice velocity measurements and provide researchers access to better information for monitoring local changes in ice glaciers. Other uses of these datasets include control and calibration of ice-sheet modelling, developments in our understanding of Antarctic ice-sheet evolution, and improvements in the fidelity of projects investigating sea-level rise (https://doi.pangaea.de/10.1594/PANGAEA.895738).


2015 ◽  
Vol 56 (70) ◽  
pp. 63-69 ◽  
Author(s):  
Denis Callens ◽  
Nicolas Thonnard ◽  
Jan T.M. Lenaerts ◽  
Jan M. Van Wessem ◽  
Willem Jan Van de Berg ◽  
...  

AbstractMass changes of polar ice sheets have an important societal impact, because they affect global sea level. Estimating the current mass budget of ice sheets is equivalent to determining the balance between surface mass gain through precipitation and outflow across the grounding line. For the Antarctic ice sheet, grounding line outflow is governed by oceanic processes and outlet glacier dynamics. In this study, we compute the mass budget of major outlet glaciers in the eastern Dronning Maud Land sector of the Antarctic ice sheet using the input/output method. Input is given by recent surface accumulation estimates (SMB) of the whole drainage basin. The outflow at the grounding line is determined from the radar data of a recent airborne survey and satellite-based velocities using a flow model of combined plug flow and simple shear. This approach is an improvement on previous studies, as the ice thickness is measured, rather than being estimated from hydrostatic equilibrium. In line with the general thickening of the ice sheet over this sector, we estimate the regional mass balance in this area at 3.15 ± 8.23 Gt a−1 according to the most recent SMB model results.


Sign in / Sign up

Export Citation Format

Share Document