SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Comparison between numerical model ice-shelf flow simulations and synthetic aperture radar (SAR) interferograms is used to study the dynamics at the Hemmen Ice Rise (HIR) and Lassiter Coast (LC) corners of the iceberg-calving front of the Filchner-Ronne Ice Shelf (FRIS).
Three stations near the calving front of the Ross Ice Shelf, Antarctica, recorded GPS data through a full spring-neap tidal cycle in November 2005. The data revealed a diurnal horizontal motion that varied both along and transverse to the long-term average velocity direction, similar to tidal signals observed in other ice shelves and ice streams. Based on its periodicity, it was hypothesized that the signal represents a flow response of the Ross Ice Shelf to the diurnal tides of the Ross Sea. To assess the influence of the tide on the ice-shelf motion, two hypotheses were developed. The first addressed the direct response of the ice shelf to tidal forcing, such as forces due to sea-surface slopes or forces due to sub-ice-shelf currents. The second involved the indirect response of ice-shelf flow to the tidal signals observed in the ice streams that source the ice shelf. A finite-element model, based on viscous creep flow, was developed to test these hypotheses, but succeeded only in falsifying both hypotheses, i.e. showing that direct tidal effects produce too small a response, and indirect tidal effects produce a response that is not smooth in time. This nullification suggests that a combination of viscous and elastic deformation is required to explain the observations.
Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.
Fifteen synthetic-aperture radar (SAR) images of the Ronne Ice Shelf, Antarctica, obtained by the European Space Agency (ESA)'s Earth Remote Sensing satellites (ERS) 1 & 2 are used to study ice-shelf dynamics near two ends of the iceberg-calving front.
We compare Earth Remote Sensing (ERS) satellite synthetic-aperture radar (SAR) interferograms with artificial interferograms constructed using output of a finite-element ice-shelf flow model to study the dynamics of Filchner-Ronne Ice Shelf (FRIS) near Hemmen Ice Rise (HIR) where the iceberg-calving front itersects Berkener Island (BI).
Multiple radar inteferograms from the ERS satellite, may be used to eliminate the influence of the long-term velocity of an ice shelf to isolate the deformation signals associated with the cyclic, vertical motion of the ice shelf under the influence of tidal forcing from the ocean.
We use radar interferometric images collected by ERS-1 and Radarsat-1 to observe the rupture tip of rifts along Hemmen Ice Rise on the Ronne Ice Shelf, Antarctica. Interferograms generated in 1992 and 1997 allow us to observe the deformation of ice accumulated over respectively 9 and24 days. The results are an important step towards developing a better model of ice-shelf calving mechanism.
Viscoelastic models of ice-shelf flexure and ice-stream velocity perturbations are combined into a single efficient flowline model to study tidal forcing of grounded ice. The magnitude and timing of icestream response to tidally driven changes in hydrostatic pressure and/or basal drag are found to depend significantly on bed rheology, with only a perfectly plastic bed allowing instantaneous velocity response at the grounding line. The model can reasonably reproduce GPS observations near the grounding zone of Bindschadler Ice Stream (formerly Ice Stream D) on semidiurnal time scales; however, other forcings such as tidally driven ice-shelf slope transverse to the flowline and flexurally driven till deformation must also be considered if diurnal motion is to be matched
ERS images of the two ends of the ice front of the Filchner-Ronne ice shelf (Antarctica) were utilized interferometrically to study the deformation rate of the ice shelf in response to viscous creep. On the western bank of Berkner Island (BI), near Hemmen Ice Rise (HIR), a time series of ERS data acquired in February 1992, in both ascending and descending mode, and with a three-day time interval, were utilized to map the ice velocity in two dimensions. Finite-element ice-shelf flow simulations are compared with the ERS interferograms to interpret the ice motion in terms of the physical constraints on ice-shelf flow. The efforts to fit artificial interferograms generated with model velocity output suggest that the flow regime is strongly influenced by three processes. First, a void-creation process responsible for rifts at coastal margins tends to uncouple the ice shelf from the ice rise and neighboring coast of BI. Secondly, sea ice within the void space appears to act as a binding agent between discrete ice-shelf fragments, allowing rigid-body rotations of these fragments. Third, strain rates appear to be enhanced in a narrow zone adjacent to HIR, implying significant strain softening along the boundary. It is believed that synthetic aperture radar (SAR) interferogram/model intercomparison represents a powerful impetus towards the development of better, more physically realistic ice-shelf flow models.
Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf–ocean modelling, the second Ice Shelf–Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0–2 (Asay-Davis et al., 2016), which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice–ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf–ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice–ocean boundary layer and ice shelf cavity circulation.
The Antarctic Ice Sheet is surrounded by large floating ice shelves that spread under their own weight into the ocean. Ice shelf rigidity depends on ice temperature and fabrics, and is influenced by ice flow and the delicate balance between bottom and surface accumulation. Here, we use an inverse control method to infer the rigidity of the Ronne Ice Shelf that best matches observations of ice velocity from satellite radar interferometry. Ice rigidity, or flow law parameter B, is shown to vary between 300 and 900 kPa a(sup 1/3). Ice is softer along the side margins due to frictional heating, and harder along the outflow of large glaciers, which advect cold continental ice. Melting at the bottom surface of the ice shelf increases its rigidity, while freezing decreases it. Accurate numerical modelling of ice shelf flow must account for this spatial variability in mechanical characteristics.
As the global climate warms, increased surface meltwater production on ice shelves may trigger ice-shelf collapse and enhance global sea-level rise. The formation of surface rivers could help prevent ice-shelf collapse if they can efficiently evacuate meltwater. Here we present observations of the evolution of a surface river into an ice-shelf estuary atop the Petermann Ice Shelf in northwest Greenland and identify a second estuary at the nearby Ryder Ice Shelf. This surface-hydrology process can foster fracturing and enhance calving. At the Petermann estuary, sea ice was observed converging at the river mouth upstream, indicating a flow reversal. Seawater persists in the estuary after the surrounding icescape is frozen. Along the base of Petermann estuary, linear fractures were initiated at the calving front and propagated upstream along the channel. Similar fractures along estuary channels shaped past large rectilinear calving events at the Petermann and Ryder ice shelves. Increased surface melting in a warming world will enhance fluvial incision, promoting estuary development and longitudinal fracturing orthogonal to ice-shelf fronts, and increase rectilinear calving. Estuaries could develop in Antarctica within the next half-century, resulting in increased calving and accelerating both ice loss and global sea-level rise.
We develop a model of an ice shelf-ice stream system as a viscoelastic beam partially supported by an elastic foundation. When bed rock near the grounding line acts as a fulcrum, leverage from the ice shelf dropping at low tide can cause significant (approx 1 cm) uplift in the first few kilometers of grounded ice.This uplift and the corresponding depression at high tide lead to basal pressure variations of sufficient magnitude to influence subglacial hydrology.Tidal flexure may thus affect basal lubrication, sediment flow, and till strength, all of which are significant factors in ice-stream dynamics and grounding-line stability. Under certain circumstances, our results suggest the possibility of seawater being drawn into the subglacial water system. The presence of sea water beneath grounded ice would significantly change the radar reflectivity of the grounding zone and complicate the interpretation of grounded versus floating ice based on ice-penetrating radar observations.
Observations from satellite and airborne platforms are combined with model calculations to infer the nature and efficiency of basal melting of the Pine Island Glacier ice shelf, West Antarctica, by ocean waters. Satellite imagery shows surface features that suggest ice-shelf-wide changes to the ocean s influence on the ice shelf as the grounding line retreated. Longitudinal profiles of ice surface and bottom elevations are analyzed to reveal a spatially dependent pattern of basal melt with an annual melt flux of 40.5 Gt/a. One profile captures a persistent set of surface waves that correlates with quasi-annual variations of atmospheric forcing of Amundsen Sea circulation patterns, establishing a direct connection between atmospheric variability and sub-ice-shelf melting. Ice surface troughs are hydrostatically compensated by ice-bottom voids up to 150m deep. Voids form dynamically at the grounding line, triggered by enhanced melting when warmer-than-average water arrives. Subsequent enlargement of the voids is thermally inefficient (4% or less) compared with an overall melting efficiency beneath the ice shelf of 22%. Residual warm water is believed to cause three persistent polynyas at the ice-shelf front seen in Landsat imagery. Landsat thermal imagery confirms the occurrence of warm water at the same locations.
The warmest water below parts of the Ross Ice Shelf resides in the lowest portion of the water column because of its high salinity. Vertical mixing caused by tidal stirring can thus induce ablation by lifting the warm but dense water into contact with the ice shelf. A numerical tidal simulation indicates that vertically well-mixed conditions predominate in the southeastern part of the sub-ice shelf cavity, where the water column thickness is small. Basal melting in this region is expected to be between 0.05 and 0.5 m/yr and will drive a thermohaline circulation having the following characteristics: high salinity shelf water (at - 1.8 C), formed by winter sea ice production in the open Ross Sea, flows along the seabed toward the tidal mixing fronts below the ice shelf; and meltwater (at -2.2 C), produced in the well-mixed region, flows out of the sub-ice shelf cavity along the ice shelf bottom. Sensitivity of this ablation process to climatic change is expected to be small because high salinity shelf water is constrained to have the sea surface freezing temperature.
Meltwater stored in ponds and crevasses can weaken and fracture ice shelves, triggering their rapid disintegration. This ice-shelf collapse results in an increased flux of ice from adjacent glaciers and ice streams, thereby raising sea level globally. However, surface rivers forming on ice shelves could potentially export stored meltwater and prevent its destructive effects. Here we present evidence for persistent active drainage networks-interconnected streams, ponds and rivers-on the Nansen Ice Shelf in Antarctica that export a large fraction of the ice shelf's meltwater into the ocean. We find that active drainage has exported water off the ice surface through waterfalls and dolines for more than a century. The surface river terminates in a 130-metre-wide waterfall that can export the entire annual surface melt over the course of seven days. During warmer melt seasons, these drainage networks adapt to changing environmental conditions by remaining active for longer and exporting more water. Similar networks are present on the ice shelf in front of Petermann Glacier, Greenland, but other systems, such as on the Larsen C and Amery Ice Shelves, retain surface water at present. The underlying reasons for export versus retention remain unclear. Nonetheless our results suggest that, in a future warming climate, surface rivers could export melt off the large ice shelves surrounding Antarctica-contrary to present Antarctic ice-sheet models, which assume that meltwater is stored on the ice surface where it triggers ice-shelf disintegration.
Beginning in February 2003, NASA's Ice, Cloud, and land Elevation Satellite (ICESat) has determined surface elevations from approx. 86degN to 86degS latitude. To date, altimetry data have been acquired in a series of observation periods in repeated track patterns using all three Geoscience Laser Altimeter System (GLAS) lasers. This paper will focus on ice shelf elevation data that were obtained in 2003 across the Brunt Ice Shelf and the Stancomb-Wills Ice Tongue. Integrating the altimetry with the recently available MODIS Mosaic of Antarctica (MOA), quantifies the relative accuracy and precision of the resulting ice shelf elevations. Furthermore, the elevation data was processed onto an elevation grid, by regional interpolation across the area s complex glacial features only. Ice thickness estimation from the altimetry of the floating ice is discussed. ICESat operates at 40Hz and its elevation data is obtained every 172m along track. These elevations have a relative accuracy of about 14cm based on the standard deviation of low-slope crossover differences and a precision of close to 2cm for the Laser 2a, Release 21, GLA12 data used here.
We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.