Search

The Thwaites Eastern Ice Shelf buttresses a significant portion of Thwaites Glacier through contact with a pinning point 40 km offshore of the present grounding line. Predicting future rates of Thwaites Glacier's contribution to sea-level rise depends on the evolution of this pinning point and the resultant change in the ice-shelf stress field since the breakup of the Thwaites Western Glacier Tongue in 2009. Here we use Landsat-8 feature tracking of ice velocity in combination with ice-sheet model perturbation experiments to show how past changes in flow velocity have been governed in large part by changes in lateral shear and pinning point interactions with the Thwaites Western Glacier Tongue. We then use recent satellite altimetry data from ICESat-2 to show that Thwaites Glacier's grounding line has continued to retreat rapidly; in particular, the grounded area of the pinning point is greatly reduced from earlier mappings in 2014, and grounded ice elevations are continuing to decrease. This loss has created two pinned areas with ice flow now funneled between them. If current rates of surface lowering persist, the Thwaites Eastern Ice Shelf will unpin from the seafloor in less than a decade, despite our finding from airborne radar data that the seafloor underneath the pinning point is about 200 m shallower than previously reported. Advection of relatively thin and mechanically damaged ice onto the remaining portions of the pinning point and feedback mechanisms involving basal melting may further accelerate the unpinning. As a result, ice discharge will likely increase up to 10 % along a 45 km stretch of the grounding line that is currently buttressed by the Thwaites Eastern Ice Shelf. [ABSTRACT FROM AUTHOR]

The Thwaites Eastern Ice Shelf (TEIS) buttresses the eastern grounded portion of Thwaites Glacier through contact with a pinning point at its seaward limit. Loss of this ice shelf will promote further acceleration of Thwaites Glacier. Understanding the dynamic controls and structural integrity of the TEIS is therefore important to estimating Thwaites' future sea-level contribution. We present a ∼ 20-year record of change on the TEIS that reveals the dynamic controls governing the ice shelf's past behaviour and ongoing evolution. We derived ice velocities from MODIS and Sentinel-1 image data using feature tracking and speckle tracking, respectively, and we combined these records with ITS_LIVE and GOLIVE velocity products from Landsat-7 and Landsat-8. In addition, we estimated surface lowering and basal melt rates using the Reference Elevation Model of Antarctica (REMA) DEM in comparison to ICESat and ICESat-2 altimetry. Early in the record, TEIS flow dynamics were strongly controlled by the neighbouring Thwaites Western Ice Tongue (TWIT). Flow patterns on the TEIS changed following the disintegration of the TWIT around 2008, with a new divergence in ice flow developing around the pinning point at its seaward limit. Simultaneously, the TEIS developed new rifting that extends from the shear zone upstream of the ice rise and increased strain concentration within this shear zone. As these horizontal changes occurred, sustained thinning driven by basal melt reduced ice thickness, particularly near the grounding line and in the shear zone area upstream of the pinning point. This evidence of weakening at a rapid pace suggests that the TEIS is likely to fully destabilize in the next few decades, leading to further acceleration of Thwaites Glacier. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Grounding zones are vital to ice-sheet mass balance and its coupling to the global ocean circulation. Processes here determine the mass discharge from the grounded ice sheet, to the floating ice shelves. The response of this transition zone to tidal forcing has been described by both elastic and viscoelastic models. Here we examine the validity of these models for grounding zone flexure over tidal timescales using field data from the Southern McMurdo Ice Shelf (78° 15′S, 167° 7′E). Observations of tidal movement were carried out by simultaneous tiltmeter and GPS measurements along a profile across the grounding zone. Finite-element simulations covering a 64 d period reveal that the viscoelastic model fits best the observations using a Young's modulus of 1.6 GPa and a viscosity of 1013.7 Pa s (≈ 50.1 TPa s). We conclude that the elastic model is only well-constrained for tidal displacements >35% of the spring-tidal amplitude using a Young's modulus of 1.62 ± 0.69 GPa, but that a viscoelastic model is necessary to adequately capture tidal bending at amplitudes below this threshold. In grounding zones where bending stresses are greater than at the Southern McMurdo Ice Shelf or ice viscosity is lower, the threshold would be even higher. [ABSTRACT FROM PUBLISHER]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Tidal flexure in ice shelf grounding zones has been used extensively in the past to determine grounding line position and ice properties. Although the rheology of ice is viscoelastic at tidal loading frequencies, most modelling studies have assumed some form of linear elastic beam approximation to match observed flexure profiles. Here we use density, radar and DInSAR measurements in combination with full-Stokes viscoelastic modelling to investigate a range of additional controls on the flexure of the Southern McMurdo Ice Shelf. We find that inclusion of observed basal crevasses and density dependent ice stiffness can greatly alter the flexure profile and yet fitting a simple elastic beam model to that profile will still produce an excellent fit. Estimates of the effective Young's modulus derived by fitting flexure profiles are shown to vary by over 200% depending on whether these factors are included, even when the local thickness is well constrained. Conversely, estimates of the grounding line position are relatively insensitive to these considerations for the case of a steep bed slope in our study region. By fitting tidal amplitudes only, and ignoring phase information, elastic beam theory can provide a good fit to observations in a wide variety of situations. This should, however, not be taken as an indication that the underlying rheological assumptions are correct. [ABSTRACT FROM PUBLISHER]