ALBERT

Description: Differential interferometric synthetic aperture radar (DInSAR) is an essential tool for detecting ice-sheet motion near Antarctica's oceanic margin. These space-borne measurements have been used extensively in the past to map the location and retreat of ice-shelf grounding lines as an indicator for the onset of marine ice-sheet instability and to calculate the mass balance of ice sheets and individual catchments. The main difficulty in interpreting DInSAR is that images originate from a combination of several SAR images and do not indicate instantaneous ice deflection at the times of satellite data acquisitions. Here, we combine the sub-centimetre accuracy and spatial benefits of DInSAR with the temporal benefits of tide models to infer the spatio-temporal dynamics of ice–ocean interaction during the times of satellite overpasses. We demonstrate the potential of this synergy with TerraSAR-X data from the almost-stagnant southern McMurdo Ice Shelf (SMIS). We then validate our algorithm with GPS data from the fast-flowing Darwin Glacier, draining the Antarctic Plateau through the Transantarctic Mountains into the Ross Sea. We are able to reconstruct DInSAR-derived vertical displacements to 7 mm mean absolute residual error and generally improve traditional tide-model output by up to 39 % from 10.8 to 6.7 cm RMSE against GPS data from areas where ice is in local hydrostatic equilibrium with the ocean and by up to 74 % from 21.4 to 5.6 cm RMSE against GPS data in feature-rich coastal areas where tide models have not been applicable before. Numerical modelling then reveals Young's modulus of E=1.0±0.56 GPa and an ice viscosity of ν=10±3.65 TPa s when finite-element simulations of tidal flexure are matched to 16 d of tiltmeter data, supporting the hypothesis that strain-dependent anisotropy may significantly decrease effective viscosity compared to isotropic polycrystalline ice on large spatial scales. Applications of our method include the following: refining coarsely gridded tide models to resolve small-scale features at the spatial resolution and vertical accuracy of SAR imagery, separating elastic and viscoelastic contributions in the satellite-derived flexure measurement, and gaining information about large-scale ice heterogeneity in Antarctic ice-shelf grounding zones, the missing key to improving current ice-sheet flow models. The reconstruction of the individual components forming DInSAR images has the potential to become a standard remote-sensing method in polar tide modelling. Unlocking the algorithm's full potential to answer multi-disciplinary research questions is desired and demands collaboration within the scientific community.

Description: We examine tidal flexure in the grounding zone of the McMurdo Ice Shelf, Antarctica, using a combination of TerraSAR-X repeat-pass radar interferometry, a precise digital elevation model, and GPS ground validation data. Satellite and field data were acquired in tandem between October and December 2014. Our GPS data show a horizontal modulation of up to 60 % of the vertical displacement amplitude at tidal periods within a few kilometres of the grounding line. We ascribe the observed oscillatory horizontal motion to varying bending stresses and account for it using a simple elastic beam model. The horizontal surface strain is removed from nine differential interferograms to obtain precise bending curves. They reveal a fixed (as opposed to tidally migrating) grounding-line position and eliminate the possibility of significant upstream bending at this location. The consequence of apparent vertical motion due to uncorrected horizontal strain in interferometric data is a systematic mislocation of the interferometric grounding line by up to the order of one ice thickness, or several hundred metres. While our field site was selected due to its simple boundary conditions and low background velocity, our findings are relevant to other grounding zones studied by satellite interferometry, particularly studies looking at tidally induced velocity changes or interpreting satellite-based flexure profiles.

Description: We examine tidal flexure in the grounding zone of the McMurdo Ice Shelf, Antarctica, using a combination of a TerraSAR-X repeat-pass radar interferometry, a precise digital elevation model, and GPS ground validation data. Satellite and field data were acquired in tandem between October and December 2014. Our GPS data show a horizontal modulation of up to 60 % of the vertical amplitude at tidal periods within a few km of the grounding line. We ascribe this to bending stresses and account for it using a simple elastic beam model. The horizontal surface strain is removed from nine differential interferograms to obtain precise vertical bending curves. This processing step allows us to identify a fixed (as opposed to tidally migrating) grounding line position and eliminates the possibility of significant upstream bending at this location. The change in apparent vertical motion due to horizontal strain can lead to a systematic mis-location of the interferometric grounding line by the order of up to one ice thickness, or several hundred metres. While our field site was selected in consideration of the simple boundary conditions and low background velocity the findings are relevant to other satellite-based grounding zone studies, particularly those looking at tidally-induced velocity changes or interpreting satellite-based flexure profiles.

Description: Differential interferometric synthetic aperture radar (DInSAR) is an essential tool for detecting ice-sheet motion near Antarctica's oceanic margin. These space-borne measurements have been used extensively in the past to map the location and retreat of ice-shelf grounding lines as an indicator for the onset of marine ice-sheet instability and to calculate the mass balance of ice-sheets and individual catchments. The main difficulty in interpreting DInSAR is that images originate from a combination of several SAR images and do not indicate instantaneous ice deflection at the time of satellite data acquisition. Here, we combine the sub-centimetre accuracy and spatial benefits of DInSAR with the temporal benefits of tide models to infer the spatiotemporal dynamics of ice-ocean interaction during the times of satellite overpasses. We demonstrate the potential of this synergy with TerraSAR-X data from the almost stagnant Southern McMurdo Ice Shelf. We then validate our algorithm with GPS data from the fast-flowing Darwin Glacier, draining the Antarctic Plateau through the Transantarctic Mountains into the Ross Sea. We are able to match DInSAR to 0.84 mm; generally improve traditional tide models by up to −39 % from 10.8 cm to 6.7 cm RMSE against GPS data from areas where ice is in local hydrostatic equilibrium with the ocean; and up to −74 % from 21.4 cm to 5.6 cm RMSE against GPS data in feature-rich coastal areas where contemporary tide-models are most inaccurate. Numerical modelling then reveals a Young’s modulus of E = 1.0 GPa and an ice viscosity of 10 TPa s when finite-element simulations of tidal flexure are matched to 16 days of tiltmeter data; supporting the theory that strain dependent anisotropy may significantly decrease effective viscosity compared to isotropic polycrystalline ice on large spatial scales. Applications of our method range from (i) refining coarsly-gridded tide models to resolve small-scale features at the spatial resolution and vertical accuracy of SAR imagery, to (ii) separating elastic and viscoelastic contributions in the satellite derived flexure measurement and (iii) gaining information about large-scale ice heterogenity in Antarctic ice-shelf grounding zones, the missing key to improve current ice-sheet flow models. The reconstruction of the individual components forming DInSAR images has the potential to become a standard remote-sensing method in polar tide modelling. Unlocking the algorithm's full potential to answer multi-disciplinary research questions is desired and demands collaboration within the scientific community.

Description: Despite of their importance for sea-level rise, seasonal water availability, and as source of geohazards, mountain glaciers are one of the few remaining sub-systems of the global climate system for which no globally applicable, open source, community-driven model exists. Here we present the Open Global Glacier Model (OGGM, http://www.oggm.org), developed to provide a modular and open source numerical model framework for simulating past and future change of any glacier in the world. The modelling chain comprises data downloading tools (glacier outlines, topography, climate, validation data), a preprocessing module, a mass-balance model, a distributed ice thickness estimation model, and an ice flow model. The monthly mass-balance is obtained from gridded climate data and a temperature index melt model. To our knowledge, OGGM is the first global model explicitly simulating glacier dynamics: the model relies on the shallow ice approximation to compute the depth-integrated flux of ice along multiple connected flowlines. In this paper, we describe and illustrate each processing step by applying the model to a selection of glaciers before running global simulations under idealized climate forcings. Even without an in-depth calibration, the model shows a very realistic behaviour. We are able to reproduce earlier estimates of global glacier volume by varying the ice dynamical parameters within a range of plausible values. At the same time, the increased complexity of OGGM compared to other prevalent global glacier models comes at a reasonable computational cost: several dozens of glaciers can be simulated on a personal computer, while global simulations realized in a supercomputing environment take up to a few hours per century. Thanks to the modular framework, modules of various complexity can be added to the codebase, allowing to run new kinds of model intercomparisons in a controlled environment. Future developments will add new physical processes to the model as well as tools to calibrate the model in a more comprehensive way. OGGM spans a wide range of applications, from ice-climate interaction studies at millenial time scales to estimates of the contribution of glaciers to past and future sea-level change. It has the potential to become a self-sustained, community driven model for global and regional glacier evolution.

Description: Despite their importance for sea-level rise, seasonal water availability, and as a source of geohazards, mountain glaciers are one of the few remaining subsystems of the global climate system for which no globally applicable, open source, community-driven model exists. Here we present the Open Global Glacier Model (OGGM), developed to provide a modular and open-source numerical model framework for simulating past and future change of any glacier in the world. The modeling chain comprises data downloading tools (glacier outlines, topography, climate, validation data), a preprocessing module, a mass-balance model, a distributed ice thickness estimation model, and an ice-flow model. The monthly mass balance is obtained from gridded climate data and a temperature index melt model. To our knowledge, OGGM is the first global model to explicitly simulate glacier dynamics: the model relies on the shallow-ice approximation to compute the depth-integrated flux of ice along multiple connected flow lines. In this paper, we describe and illustrate each processing step by applying the model to a selection of glaciers before running global simulations under idealized climate forcings. Even without an in-depth calibration, the model shows very realistic behavior. We are able to reproduce earlier estimates of global glacier volume by varying the ice dynamical parameters within a range of plausible values. At the same time, the increased complexity of OGGM compared to other prevalent global glacier models comes at a reasonable computational cost: several dozen glaciers can be simulated on a personal computer, whereas global simulations realized in a supercomputing environment take up to a few hours per century. Thanks to the modular framework, modules of various complexity can be added to the code base, which allows for new kinds of model intercomparison studies in a controlled environment. Future developments will add new physical processes to the model as well as automated calibration tools. Extensions or alternative parameterizations can be easily added by the community thanks to comprehensive documentation. OGGM spans a wide range of applications, from ice–climate interaction studies at millennial timescales to estimates of the contribution of glaciers to past and future sea-level change. It has the potential to become a self-sustained community-driven model for global and regional glacier evolution.