scholarly journals Ice-shelf basal channels in a coupled ice/ocean model

2012 ◽  
Vol 58 (212) ◽  
pp. 1227-1244 ◽  
Author(s):  
Carl V. Gladish ◽  
David M. Holland ◽  
Paul R. Holland ◽  
Stephen F. Price

AbstractA numerical model for an interacting ice shelf and ocean is presented in which the ice- shelf base exhibits a channelized morphology similar to that observed beneath Petermann Gletscher’s (Greenland) floating ice shelf. Channels are initiated by irregularities in the ice along the grounding line and then enlarged by ocean melting. To a first approximation, spatially variable basal melting seaward of the grounding line acts as a steel-rule die or a stencil, imparting a channelized form to the ice base as it passes by. Ocean circulation in the region of high melt is inertial in the along-channel direction and geostrophically balanced in the transverse direction. Melt rates depend on the wavelength of imposed variations in ice thickness where it enters the shelf, with shorter wavelengths reducing overall melting. Petermann Gletscher’s narrow basal channels may therefore act to preserve the ice shelf against excessive melting. Overall melting in the model increases for a warming of the subsurface water. The same sensitivity holds for very slight cooling, but for cooling of a few tenths of a degree a reorganization of the spatial pattern of melting leads, surprisingly, to catastrophic thinning of the ice shelf 12 km from the grounding line. Subglacial discharge of fresh water along the grounding line increases overall melting. The eventual steady state depends on when discharge is initiated in the transient history of the ice, showing that multiple steady states of the coupled system exist in general.

2021 ◽  
Author(s):  
Chen Zhao ◽  
Rupert Gladstone ◽  
Ben Galton-Fenzi ◽  
David Gwyther

<p>The ocean-driven basal melting has important implications for the stability of ice shelves in Antarctic, which largely affects the ice sheet mass balance, ocean circulation, and subsequently global sea level rise. Due to the limited observations in the ice shelf cavities, the couple ice sheet ocean models have been playing a critical role in examining the processes governing basal melting. In this study we use the Framework for Ice Sheet-Ocean Coupling (FISOC) to couple the Elmer/Ice full-stokes ice sheet model and the Regional Ocean Modeling System (ROMS) ocean model to model ice shelf/ocean interactions for an idealised three-dimensional domain. Experiments followed the coupled ice sheet–ocean experiments under the first phase of the Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1). A periodic pattern in the simulated mean basal melting rates is found to be highly consistent with the maximum barotropic stream function and also the grounding line retreat row by row,  which is likely to be related with the gyre break down near the grounding line caused by some non-physical instability events from the ocean bottom. Sensitivity tests are carried out, showing that this periodic pattern is not sensitive to the choice of couple time intervals and horizontal eddy viscosities but sensitive to vertical resolution in the ocean model, the chosen critical water column thickness in the wet-dry scheme, and the tracer properties for the nudging dry cells at the ice-ocean interface boundary. Further simulations are necessary to better explain the mechanism involved in the couple ice-ocean system, which is very significant for its application on the realistic ice-ocean systems in polar regions.</p>


2018 ◽  
Vol 12 (2) ◽  
pp. 453-476 ◽  
Author(s):  
Rachael D. Mueller ◽  
Tore Hattermann ◽  
Susan L. Howard ◽  
Laurie Padman

Abstract. Recent modeling studies of ocean circulation in the southern Weddell Sea, Antarctica, project an increase over this century of ocean heat into the cavity beneath Filchner–Ronne Ice Shelf (FRIS). This increase in ocean heat would lead to more basal melting and a modification of the FRIS ice draft. The corresponding change in cavity shape will affect advective pathways and the spatial distribution of tidal currents, which play important roles in basal melting under FRIS. These feedbacks between heat flux, basal melting, and tides will affect the evolution of FRIS under the influence of a changing climate. We explore these feedbacks with a three-dimensional ocean model of the southern Weddell Sea that is forced by thermodynamic exchange beneath the ice shelf and tides along the open boundaries. Our results show regionally dependent feedbacks that, in some areas, substantially modify the melt rates near the grounding lines of buttressed ice streams that flow into FRIS. These feedbacks are introduced by variations in meltwater production as well as the circulation of this meltwater within the FRIS cavity; they are influenced locally by sensitivity of tidal currents to water column thickness (wct) and non-locally by changes in circulation pathways that transport an integrated history of mixing and meltwater entrainment along flow paths. Our results highlight the importance of including explicit tidal forcing in models of future mass loss from FRIS and from the adjacent grounded ice sheet as individual ice-stream grounding zones experience different responses to warming of the ocean inflow.


2017 ◽  
Author(s):  
Rachael D. Mueller ◽  
Tore Hattermann ◽  
Susan L. Howard ◽  
Laurence Padman

Abstract. Recent modeling studies of ocean circulation in the southern Weddell Sea, Antarctica, project an increase over this century of ocean heat into the cavity beneath Filchner-Ronne Ice Shelf (FRIS). This increase in ocean heat would lead to more basal melting and a modification of the FRIS ice draft. The corresponding change in cavity shape will affect advective pathways and the spatial distribution of tidal currents, which play important roles in basal melting under FRIS. These feedbacks between heat flux, basal melting, and tides will affect the evolution of FRIS under the influence of a changing climate. We explore these feedbacks with a three-dimensional ocean model of the southern Weddell Sea that is forced by thermodynamic exchange beneath the ice shelf and tides along the open boundaries. Our results show regionally-dependent feedbacks that, in some areas, substantially modify the melt rates near the grounding lines of buttressed ice streams that flow into FRIS. These feedbacks are introduced by variations in meltwater production as well as the circulation of this meltwater within the FRIS cavity; they are influenced locally by sensitivity of tidal currents to water column thickness and non-locally by changes in circulation pathways that transport an integrated history of mixing and meltwater entrainment along flow paths. Our results highlight the importance of including explicit tidal forcing in models of future mass loss from FRIS and from the adjacent grounded ice sheet as individual ice stream grounding zones experience different responses to warming of the ocean inflow.


2021 ◽  
Author(s):  
Vjeran Visnjevic ◽  
Reinhard Drews ◽  
Clemens Schannwell ◽  
Inka Koch

<p>Ice shelves buttress ice flow from the continent towards the ocean, and their disintegration results in increased ice discharge.  Ice-shelf evolution and integrity is influenced by surface accumulation, basal melting, and ice dynamics. We find signals of all of these processes imprinted in the ice-shelf stratigraphy that can be mapped using isochrones imaged with radar.</p><p>Our aim is to develop an inverse approach to infer ice shelf basal melt rates using radar isochrones as observational constraints. Here, we investigate the influence of basalt melt rates on the shape of isochrones using combined insights from both forward and inverse modeling. We use the 3D full Stokes model Elmer/Ice in our forward simulations, aiming to reproduce isochrone patterns observed in our data. Moreover we develop an inverse approach based on the shallow shelf approximating, aiming to constrain basal melt rates using isochronal radar data and surface velocities. Insights obtained from our simulations can also guide the collection of new radar data (e.g., profile lines along vs. across-flow) in a way that ambiguities in interpreting the ice-shelf stratigraphy can be minimized. Eventually, combining these approaches will enable us to better constrain the magnitude and history of basal melting, which will give valuable input for ocean circulation and sea level rise projections.</p>


1988 ◽  
Vol 11 ◽  
pp. 165-172 ◽  
Author(s):  
R. H. Thomas ◽  
S. N. Stephenson ◽  
R. A. Bindschadler ◽  
S. Shabtaie ◽  
C. R. Bentley

Detailed measurements of surface topography, ice motion, snow accumulation, and ice thickness were made in January 1974 and again in December 1984, along an 8 km stake network extending from the ice sheet, across the grounding line, and on to floating ice shelf in the mouth of slow-moving Ice Stream C, which flows into the eastern side of Ross Ice Shelf, Antarctica. During the 11 years between surveys, the grounding line retreated by approximately 300 m. This was caused by net thinning of the ice shelf, which we believe to be a response to the comparatively recent, major decrease in ice discharge from Ice Stream C. Farther inland, snow accumulation is not balanced by ice discharge, and the ice stream is growing progressively thicker.There is evidence that the adjacent Ice Stream B has slowed significantly over the last decade, and this may be an early indication that this fast-moving ice stream is about to enter a period of stagnation similar to that of Ice Stream C. Indeed, these large ice streams flowing from West Antarctica into Ross Ice Shelf may oscillate between periods of relative stagnation and major activity. During active periods, large areas of ice shelf thicken and run aground on seabed to form extensive “ice plains” in the mouth of the ice stream. Ultimately, these become too large to be pushed seaward by the ice stream, which then slows down and enters a period of stagnation. During this period, the grounding line of the ice plain retreats, as we observe today in the mouth of Ice Stream C, because nearby ice shelf, no longer compressed by ice-stream motion, progressively thins. At the same time, water within the deformable till beneath the ice starts to freeze on to the base of the ice stream, and snow accumulation progressively increases the ice thickness. A new phase of activity would be initiated when the increasing gravity potential of the ice stream exceeds the total resistance of the shrinking ice plain and the thinning layer of deformable till at the bed. This could occur rapidly if the effects of the shrinking ice plain outweigh those of the thinning (and therefore stiffening) till. Otherwise, the till layer would finally become completely frozen, and the ice stream would have to thicken sufficiently to initiate significant heating by internal deformation, followed by basal melting and finally saturation of an adequate thickness of till; this could take some thousands of years.


2019 ◽  
Vol 13 (10) ◽  
pp. 2633-2656 ◽  
Author(s):  
David E. Shean ◽  
Ian R. Joughin ◽  
Pierre Dutrieux ◽  
Benjamin E. Smith ◽  
Etienne Berthier

Abstract. Ocean-induced basal melting is responsible for much of the Amundsen Sea Embayment ice loss in recent decades, but the total magnitude and spatiotemporal evolution of this melt is poorly constrained. To address this problem, we generated a record of high-resolution digital elevation models (DEMs) for Pine Island Glacier (PIG) using commercial sub-meter satellite stereo imagery and integrated additional 2002–2015 DEM and altimetry data. We implemented a Lagrangian elevation change (Dh∕Dt) framework to estimate ice shelf basal melt rates at 32–256 m resolution. We describe this methodology and consider basal melt rates and elevation change over the PIG ice shelf and lower catchment from 2008 to 2015. We document the evolution of Eulerian elevation change (dh∕dt) and upstream propagation of thinning signals following the end of rapid grounding line retreat around 2010. Mean full-shelf basal melt rates for the 2008–2015 period were ∼82–93 Gt yr−1, with ∼200–250 m yr−1 basal melt rates within large channels near the grounding line, ∼10–30 m yr−1 over the main shelf, and ∼0–10 m yr−1 over the North shelf and South shelf, with the notable exception of a small area with rates of ∼50–100 m yr−1 near the grounding line of a fast-flowing tributary on the South shelf. The observed basal melt rates show excellent agreement with, and provide context for, in situ basal melt-rate observations. We also document the relative melt rates for kilometer-scale basal channels and keels at different locations on the ice shelf and consider implications for ocean circulation and heat content. These methods and results offer new indirect observations of ice–ocean interaction and constraints on the processes driving sub-shelf melting beneath vulnerable ice shelves in West Antarctica.


2002 ◽  
Vol 48 (160) ◽  
pp. 9-19 ◽  
Author(s):  
David M. Holland

AbstractThe freezing of sea water to the base of an ice shelf can give rise to large patches of accumulated ice, a phenomenon known as marine ice. In this study a numerical method is presented for calculating the thickness of the marine-ice layer using an ice- shelf-ocean model. The present-day modeling paradigm of ice-shelf–ocean interaction usually involves the fixed specification of the ice-shelf geometry while the ocean circulation in the cavity beneath the ice shelf evolves freely. This approach relies on several assumptions, such as steady-state ice-shelf thickness and ice-shelf flow fields, in order to make reasonable quantitative estimates of the thermodynamic exchange processes occurring at the ice-shelf base. This paper discusses the impact of these and other assumptions on the estimation of the thickness of the marine-ice layer. Model simulation results are presented for an idealized ice-shelf–ocean configuration as a demonstration of the feasibility of the numerical method. A sensitivity analysis is given so as to quantify the relative uncertainty in the marine-ice thickness that arises from uncertainties in the model input parameters, these being principally the ice-shelf flow field, the basal accumulation rate and the ice-shelf thickness field.


2006 ◽  
Vol 52 (178) ◽  
pp. 325-346 ◽  
Author(s):  
Adrian Jenkins ◽  
Hugh F.J. Corr ◽  
Keith W. Nicholls ◽  
Craig L. Stewart ◽  
Christopher S.M. Doake

AbstractPrecise measurements of basal melting have been made at a series of 14 sites lying within a few kilometres of the grounding line of the Ronne Ice Shelf, Antarctica, where the ice thickness ranges from 1570 to 1940 m. The study was conducted over the course of 1 year and included a detailed survey of the horizontal deformation, as well as phase-sensitive radar measurements of the vertical displacement of both internal reflecting horizons and the ice-shelf base. Results from the surface survey show that the long-term viscous strain rate is modulated at tidal frequencies by (probably) elastic strains of order 10−5 per metre of tidal elevation. The radar measurements show a similar modulation of the long-term thinning/thickening of the ice shelf, with thickness oscillations up to a few centimetres in range. The long-term trends in ice thickness determined at points moving with the ice-shelf flow are consistent with a steady-state thickness profile. Vertical strain rates within the ice shelf were determined from the relative motion of internal reflectors. At two sites the observations were sufficient to discern the effect of tidal bending about a neutral surface 60% of the way down the ice column. Coincident measurements of horizontal and vertical strain imply a Poisson’s ratio of 0.5, and this combined with the asymmetric bending gives rise to the observed oscillations in thickness. At a number of sites the longterm viscous strain rates were found to be a linear function of depth. For an ice shelf this is an unexpected result. It can be attributed to the presence of significant vertical shear stresses set up close to the grounding line where the ice is still adjusting to flotation. Additional vertical motion arising from firn compaction was observed within the upper layers of the ice shelf. The additional motion was consistent with the assumption that firn density is a function only of the time since burial by steady surface accumulation. With both spatial and temporal fluctuations in the vertical strain rate accurately quantified it was possible to estimate the vertical motion of the ice-shelf base in response. Differences between the calculated and observed motion of the basal reflector arise because of basal melting. Derived melt rates at the 14 sites ranged from −0.11 ±0.31 to 2.51 ±0.10 m a−1, with a mean of 0. 85 m a−1 and a standard deviation of 0.69 m a−1, and showed no signs of significant sub-annual temporal variability. There was no obvious global correlation with either ice thickness or distance from the grounding line, although melt rates tended to decrease downstream along each of the flowlines studied. Previous estimates of basal melting in this region have been obtained indirectly from an assumption that the ice shelf is locally in equilibrium and have included a broad range of values. Only those at the lower end of the published range are consistent with the directly measured melt rates reported here.


1998 ◽  
Vol 27 ◽  
pp. 75-80 ◽  
Author(s):  
M.J.M. Williams ◽  
R. C. Warner ◽  
W. F. Budd

Using a three-dimensional ocean model specially adapted to the ocean cavity under the Amery Ice Shelf, we investigated the present ocean circulation and pattern of ice-shelf basal melting and freezing, the differences which would result from temperature changes in the seas adjacent to the Amery Ice Shelf, and the ramifications of these changes for the mass balance of the ice shelf. Under present conditions we estimate the net loss from the Amery Ice Shelf from excess basal melting over freezing at approximately 7.8 Gt a−1. This comprises a gross loss of 11.4 Gt a−1 at a mean rate of 0.42 m a−1, which is partially offset by freezing-on of 3.6 Gt a−1, at a mean rate of 0.19 m a−1. When the adjacent seas were assumed to warm by 1°C, we found the net melt increased to 31.6 Gt a−1, comprising 34.6 Gt a−1 of gross melt and 3.0 Gt a−1 of freezing.


2019 ◽  
Vol 13 (3) ◽  
pp. 1043-1049 ◽  
Author(s):  
Tyler Pelle ◽  
Mathieu Morlighem ◽  
Johannes H. Bondzio

Abstract. Basal melting at the bottom of Antarctic ice shelves is a major control on glacier dynamics, as it modulates the amount of buttressing that floating ice shelves exert onto the ice streams feeding them. Three-dimensional ocean circulation numerical models provide reliable estimates of basal melt rates but remain too computationally expensive for century-scale projections. Ice sheet modelers therefore routinely rely on simplified parameterizations based on either ice shelf depth or more sophisticated box models. However, existing parameterizations do not accurately resolve the complex spatial patterns of sub-shelf melt rates that have been observed over Antarctica's ice shelves, especially in the vicinity of the grounding line, where basal melting is one of the primary drivers of grounding line migration. In this study, we couple the Potsdam Ice-shelf Cavity mOdel (PICO, Reese et al., 2018) to a buoyant plume melt rate parameterization (Lazeroms et al., 2018) to create PICOP, a novel basal melt rate parameterization that is easy to implement in transient ice sheet numerical models and produces a melt rate field that is in excellent agreement with the spatial distribution and magnitude of observations for several ocean basins. We test PICOP on the Amundsen Sea sector of West Antarctica, Totten, and Moscow University ice shelves in East Antarctica and the Filchner-Ronne Ice Shelf and compare the results to PICO. We find that PICOP is able to reproduce inferred high melt rates beneath Pine Island, Thwaites, and Totten glaciers (on the order of 100 m yr−1) and removes the “banding” pattern observed in melt rates produced by PICO over the Filchner-Ronne Ice Shelf. PICOP resolves many of the issues contemporary basal melt rate parameterizations face and is therefore a valuable tool for those looking to make future projections of Antarctic glaciers.


Sign in / Sign up

Export Citation Format

Share Document