Ice Rise Research Articles

Subglacial controls on the flow of Institute Ice Stream, West Antarctica

ABSTRACTThe Institute Ice Stream (IIS) rests on a reverse-sloping bed, extending >150 km upstream into the ~1.8 km deep Robin Subglacial Basin, placing it at the threshold of marine ice-sheet instability. Understanding IIS vulnerability has focused on the effect of grounding-line melting, which is forecast to increase significantly this century. Changes to ice-flow dynamics are also important to IIS stability, yet little is known about them. Here we reveal that the trunk of the IIS occurs downstream of the intersection of three discrete subglacial features; a large ‘active’ subglacial lake, a newly-discovered sharp transition to a zone of weak basal sediments and a major tectonic rift. The border of IIS trunk flow is confined by the sediment on one side, and by a transition between basal melting and freezing at the border with the Bungenstock Ice Rise. By showing how basal sediment and water dictate present-day flow of IIS, we reveal that ice-sheet stability here is dependent on this unusual arrangement.

Temporally stable surface mass balance asymmetry across an ice rise derived from radar internal reflection horizons through inverse modeling

ABSTRACTIce rises are locally grounded parts of Antarctic ice shelves that play an important role in regulating ice flow from the continent towards the ocean. Because they protrude out of the otherwise horizontal ice shelves, ice rises induce an orographic uplift of the atmospheric flow, resulting in an asymmetric distribution of the surface mass balance (SMB). Here, we combine younger and older internal reflection horizons (IRHs) from radar to quantify this distribution in time and space across Derwael Ice Rise (DIR), Dronning Maud Land, Antarctica. We employ two methods depending on the age of the IRHs, i.e. the shallow layer approximation for the younger IRHs near the surface and an optimization technique based on an ice flow model for the older IRHs. We identify an SMB ratio of 2.5 between the flanks and the ice divide with the SMB ranging between 300 and 750 kg m−2 a−1. The SMB maximum is located on the upwind side, ~4 km offset to today's topographic divide. The large-scale asymmetry is consistently observed in time until 1966. The SMB from older IRHs is less-well constrained, but the asymmetry has likely persisted for >ka, indicating that DIR has been a stable features over long time spans.

Airborne radar evidence for tributary flow switching in Institute Ice Stream, West Antarctica: Implications for ice sheet configuration and dynamics

AbstractDespite the importance of ice streaming to the evaluation of West Antarctic Ice Sheet (WAIS) stability we know little about mid‐ to long‐term dynamic changes within the Institute Ice Stream (IIS) catchment. Here we use airborne radio echo sounding to investigate the subglacial topography, internal stratigraphy, and Holocene flow regime of the upper IIS catchment near the Ellsworth Mountains. Internal layer buckling within three discrete, topographically confined tributaries, through Ellsworth, Independence, and Horseshoe Valley Troughs, provides evidence for former enhanced ice sheet flow. We suggest that enhanced ice flow through Independence and Ellsworth Troughs, during the mid‐Holocene to late Holocene, was the source of ice streaming over the region now occupied by the slow‐flowing Bungenstock Ice Rise. Although buckled layers also exist within the slow‐flowing ice of Horseshoe Valley Trough, a thicker sequence of surface‐conformable layers in the upper ice column suggests slowdown more than ~4000 years ago, so we do not attribute enhanced flow switch‐off here, to the late Holocene ice‐flow reorganization. Intensely buckled englacial layers within Horseshoe Valley and Independence Troughs cannot be accounted for under present‐day flow speeds. The dynamic nature of ice flow in IIS and its tributaries suggests that recent ice stream switching and mass changes in the Siple Coast and Amundsen Sea sectors are not unique to these sectors, that they may have been regular during the Holocene and may characterize the decline of the WAIS.

Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points

Abstract. In ice flow modelling, the use of control methods to assimilate the dynamic and geometric state of an ice body has become common practice. These methods have primarily focussed on inverting for one of the two least known properties in glaciology, namely the basal friction coefficient or the ice viscosity parameter. Here, we present an approach to infer both properties simultaneously for the whole of the Antarctic ice sheet. After the assimilation, the root-mean-square deviation between modelled and observed surface velocities attains 8.7 m a−1 for the entire domain, with a slightly higher value of 14.0 m a−1 for the ice shelves. An exception in terms of the velocity mismatch is the Thwaites Glacier Ice Shelf, where the RMS value is almost 70 m a−1. The reason is that the underlying Bedmap2 geometry ignores the presence of an ice rise, which exerts major control on the dynamics of the eastern part of the ice shelf. On these grounds, we suggest an approach to account for pinning points not included in Bedmap2 by locally allowing an optimisation of basal friction during the inversion. In this way, the velocity mismatch on the ice shelf of Thwaites Glacier is more than halved. A characteristic velocity mismatch pattern emerges for unaccounted pinning points close to the marine shelf front. This pattern is exploited to manually identify seven uncharted features around Antarctica that exert significant resistance to the shelf flow. Potential pinning points are detected on Fimbul, West, Shackleton, Nickerson and Venable ice shelves. As pinning points can provide substantial resistance to shelf flow, with considerable consequences if they became ungrounded in the future, the model community is in need of detailed bathymetry there. Our data assimilation points to some of these dynamically important features not present in Bedmap2 and implicitly quantifies their relevance.

Antarctic ice rise formation, evolution, and stability

AbstractAntarctic ice rises originate from the contact between ice shelves and one of the numerous topographic highs emerging from the edge of the continental shelf. While investigations of the Raymond effect indicate their millennial‐scale stability, little is known about their formation and their role in ice shelf stability. Here we present for the first time the simulation of an ice rise using the BISICLES model. The numerical results successfully reproduce several field‐observable features, such as the substantial thinning downstream of the ice rise and the successive formation of a promontory and ice rise with stable radial ice flow center, showing that ice rises are formed during the ice sheet deglaciation. We quantify the ice rise buttressing effect, found to be mostly transient, delaying grounding line retreat significantly but resulting in comparable steady state positions. We demonstrate that ice rises are key in controlling simulations of Antarctic deglaciation.

Oceanic and atmospheric forcing of Larsen C Ice-Shelf thinning

Abstract. The catastrophic collapses of Larsen A and B ice shelves on the eastern Antarctic Peninsula have caused their tributary glaciers to accelerate, contributing to sea-level rise and freshening the Antarctic Bottom Water formed nearby. The surface of Larsen C Ice Shelf (LCIS), the largest ice shelf on the peninsula, is lowering. This could be caused by unbalanced ocean melting (ice loss) or enhanced firn melting and compaction (englacial air loss). Using a novel method to analyse eight radar surveys, this study derives separate estimates of ice and air thickness changes during a 15-year period. The uncertainties are considerable, but the primary estimate is that the surveyed lowering (0.066 ± 0.017 m yr−1) is caused by both ice loss (0.28 ± 0.18 m yr−1) and firn-air loss (0.037 ± 0.026 m yr−1). The ice loss is much larger than the air loss, but both contribute approximately equally to the lowering because the ice is floating. The ice loss could be explained by high basal melting and/or ice divergence, and the air loss by low surface accumulation or high surface melting and/or compaction. The primary estimate therefore requires that at least two forcings caused the surveyed lowering. Mechanisms are discussed by which LCIS stability could be compromised in the future. The most rapid pathways to collapse are offered by the ungrounding of LCIS from Bawden Ice Rise or ice-front retreat past a "compressive arch" in strain rates. Recent evidence suggests that either mechanism could pose an imminent risk.

Ice‐flow structure and ice dynamic changes in the Weddell Sea sector of West Antarctica from radar‐imaged internal layering

Recent studies have aroused concerns over the potential for ice draining the Weddell Sea sector of West Antarctica to figure more prominently in sea-level contributions should buttressing from the Filchner-Ronne Ice Shelf diminish. To improve understanding of how ice-stream dynamics there evolved through the Holocene, we interrogate Radio-Echo Sounding (RES) data from across the catchments of Institute and Moller Ice Streams (IIS and MIS), focusing especially on the use of internal layering to investigate ice-flow change. As an important component of this work, we investigate the influence that the orientation of the RES acquisition-track with respect to ice flow exerts on internal layering, and find that this influence is minimal unless a RES flight track parallels ice flow. We also investigate potential changes to internal layering characteristics with depth to search for important temporal transitions in ice-flow regime. Our findings suggest that ice in northern IIS, draining the Ellsworth Subglacial Highlands, has retained its present ice-flow configuration throughout the Holocene. This contrasts with less topographically-constrained ice in southern IIS and much of MIS, whose internal layering evinces spatial changes to the configuration of ice flow over the past ~10,000 years. Our findings confirm Siegert et al.’s (2013) inference that fast flow was diverted from Bungenstock Ice Rise during the Late Holocene, and suggest that this may have represented just one component of wider regional changes to ice flow occurring across the IIS and MIS catchments as the West Antarctic Ice Sheet has thinned since the Last Glacial Maximum.

Evolution of Derwael Ice Rise in Dronning Maud Land, Antarctica, over the last millennia

AbstractIce rises situated in the ice‐shelf belt around Antarctica have a spatially confined flow regime with local ice divides. Beneath the divides, ice stratigraphy often develops arches with amplitudes that record the divide's horizontal residence time and surface elevation changes. To investigate the evolution of Derwael Ice Rise, Dronning Maud Land, Antarctica, we combine radar and GPS data from three consecutive surveys, with a two‐dimensional, full Stokes, thermomechanically coupled, transient ice‐flow model. We find that the surface mass balance (SMB) is higher on the upwind and lower on the downwind slopes. Near the crest, the SMB is anomalously low and causes arches to form in the shallow stratigraphy, observable by radar. In deeper ice, arches are consequently imprinted by both SMB and ice rheology (Raymond effect). The data show how arch amplitudes decrease as along‐ridge slope increases, emphasizing that the lateral positioning of radar cross sections is important for the arch interpretation. Using the model with three rheologies (isotropic with n=3,4.5 and anisotropic with n=3), we show that Derwael Ice Rise is close to steady state but is best explained using ice anisotropy and moderate thinning. Our preferred, albeit not unique, scenario suggests that the ice divide has existed for at least 5000 years and lowered at approximately 0.03 m a−1 over the last 3400 years. Independent of the specific thinning scenario, our modeling suggests that Derwael Ice Rise has exhibited a local flow regime at least since the Mid‐Holocene.

Grounding-line dynamics and margin lakes

AbstractAt both corners of the now stagnant Kamb Ice Stream (KIS, West Antarctica) outlet, shear margins of different ages confine wedge-shaped areas with relatively flat, smooth surfaces that stagnated before the main trunk of the ice stream. We identify these features as lakes or past lakes, and consider scenarios for their development in a regional history of ongoing adjustment to grounding-line position. We focus here on the centuries leading up to the recent stagnation of KIS, a time when its grounding line appears to have advanced >100km from an earlier upstream location. Starting from stagnation of Crary Ice Rise and changes in the grounding zone of Whillans Ice Stream, we trace feedbacks associated with local thickening, ice grounding and thickness transients that both advance the grounding line and leave remnant lakes in their wake. These lakes in turn promote the development of secondary margins that may appear as ‘margin jumps’ in the ice record.

Modelling of Kealey Ice Rise, Antarctica, reveals stable ice-flow conditions in East Ellsworth Land over millennia

AbstractFlow at ice divides, their shape, size and internal structure depend not only on local conditions, but also on the flow regimes and past histories of the surrounding ice masses. Here we use field data from Kealey Ice Rise, Ellsworth Land, West Antarctica, in combination with flow modelling to investigate any possible signs of transients in the flow of the surrounding ice masses. Kealey Ice Rise shows linear surface features running parallel to its ridge in satellite imagery and a conspicuous layering in the ground-penetrating radar data known as double-peaked Raymond bumps. Through numerical modelling, by using an anisotropic full-Stokes thermomechanical flow solver, we analyse the evolution of Kealey Ice Rise and the timescales involved. We conclude that the features observed in the stratigraphy of Kealey Ice Rise require at least 3 ka of near-stationary flow conditions. However, we cannot exclude the possibility of a recent flow reorganization in the last century. We stress that the signs of stationary flow in radar stratigraphy observed in Kealey Ice Rise have been observed in other ice divides in the East Ellsworth Land area, suggesting stationary flow conditions over a millennial timescale in the region.

Creep deformation and buttressing capacity of damaged ice shelves: theory and application to Larsen C ice shelf

Abstract. Around the perimeter of Antarctica, much of the ice sheet discharges to the ocean through floating ice shelves. The buttressing provided by ice shelves is critical for modulating the flux of ice into the ocean, and the presently observed thinning of ice shelves is believed to be reducing their buttressing capacity and contributing to the acceleration and thinning of the grounded ice sheet. However, relatively little attention has been paid to the role that fractures play in the ability of ice shelves to sustain and transmit buttressing stresses. Here, we present a new framework for quantifying the role that fractures play in the creep deformation and buttressing capacity of ice shelves. We apply principles of continuum damage mechanics to derive a new analytical relation for the creep of an ice shelf that accounts for the softening influence of fractures on longitudinal deformation using a state damage variable. We use this new analytical relation, combined with a temperature calculation for the ice, to partition an inverse method solution for ice shelf rigidity into independent solutions for softening damage and stabilizing backstress. Using this new approach, field and remote sensing data can be utilized to monitor the structural integrity of ice shelves, their ability to buttress the flow of ice at the grounding line, and thus their indirect contribution to ice sheet mass balance and global sea level. We apply this technique to the Larsen C ice shelf using remote sensing and Operation IceBridge data, finding damage in areas with known crevasses and rifts. Backstress is highest near the grounding line and upstream of ice rises, in agreement with patterns observed on other ice shelves. The ice in contact with the Bawden ice rise is weakened by fractures, and additional damage or thinning in this area could diminish the backstress transmitted upstream. We model the consequences for the ice shelf if it loses contact with this small ice rise, finding that flow speeds would increase by 25% or more over an area the size of the former Larsen B ice shelf. Such a perturbation could potentially destabilize the northern part of Larsen C along pre-existing lines of weakness, highlighting the importance of the feedback between buttressing and fracturing in an ice shelf.

Real‐Time Changes in Arctic Ice Presented in Online Portal

The contemporary Arctic cryosphere is in a highly transient state, responding to a persistent and accelerating climate forcing. Over the past decade, the frequency with which annual ice losses in Greenland reached record highs and annual sea ice area in the Arctic reached record lows has been exceptional: Each year, the Greenland ice sheet has lost roughly as much ice, if melted, to fill most of Lake Erie, and since 2007, all sea ice minima have been smaller than those seen between 1979 and 2006. The anticipated consequences of cryospheric change, including more solar radiation absorbed by newly exposed patches of seawater where there once was ice and substantial sea level rise thanks to melting glaciers, provide a strong motivation to understand the state of Arctic ice.

Late Holocene ice-flow reconfiguration in the Weddell Sea sector of West Antarctica

Here we report Late Holocene ice sheet and grounding-line changes to the Weddell Sea sector of West Antarctica. Internal radio-echo layering within the Bungenstock Ice Rise, which comprises very slow-flowing ice separating the fast-flowing Institute and Möller ice streams, reveals ice deformed by former enhanced flow, overlain by un-deformed ice. The ice-rise surface is traversed by surface lineations explicable as diffuse ice-flow generated stripes, which thus capture the direction of flow immediately prior to the creation of the ice rise. The arrangement of internal layers can be explained by adjustment to the flow path of the Institute Ice Stream, during either a phase of ice sheet retreat not longer than ∼4000 years ago or by wholesale expansion of the grounding-line from an already retreated situation not sooner than ∼400 years ago. Some combination of these events, involving uplift of the ice rise bed during ice stream retreat and reorganisation, is also possible. Whichever the case, the implication is that the ice sheet upstream of the Bungenstock Ice Rise, which currently grounds over a >1.5 km deep basin has been, and therefore may be, susceptible to significant change.

Recent changes in the flow of the Ross Ice Shelf, West Antarctica

Comparison of surface velocities measured during the Ross Ice Shelf Geophysical and Glaciological Survey (RIGGS, 1973 to 1978) and velocities measured via feature tracking between two Moderate-resolution Imaging Spectroradiometer (MODIS) mosaics (compiled from 2003/4 and 2008/9 images) reveals widespread slowing and minor areas of acceleration in the Ross Ice Shelf (RIS) over the approximately 30 year interval. The largest changes (−13 ma−2) occur near the Whillans and Mercer Ice Streams grounding line in the southernmost part of the ice shelf. Speed has increased over the interval (up to 5 ma−2) between the MacAyeal Ice Stream grounding line and the shelf front, and along the eastern shelf front. Changes in ice thickness computed using ICESat laser altimetry are used together with a well-tested model of the ice shelf to investigate underlying causes of change in the flow of the ice shelf over time. The observed transients represent a combination of recent forcings and ongoing response to ice stream discharge variations over the past millennium. While evidence of older events may be present, the modern signal is dominated by shorter time scale events, including the stagnation of Kamb Ice Stream about 160 years ago, recent changes in basal drag on the Whillans Ice Stream ice plain and, perhaps, iceberg calving. Details in embayment geometry, for example the shallow sea floor below Crary Ice Rise, modulate the spatial pattern of ice shelf response to boundary condition perturbations.

Glacial geological implications of overconsolidated sediments on the Lomonosov Ridge and Yermak Plateau

With the coupled use of multibeam swath bathymetry, high-resolution subbottom profiling and sediment coring from icebreakers in the Arctic Ocean, there is a growing awareness of the prevalence of Quaternary ice-grounding events on many of the topographic highs found in present water depths of <1000 m. In some regions, such as the Lomonosov Ridge and Yermak Plateau, overconsolidated sediments sampled through either drilling or coring are found beneath seismically imaged unconformities of glacigenic origin. However, there exists no comprehensive analysis of the geotechnical properties of these sediments, or how their inferred stress state may be related to different glacigenic processes or types of ice-loading. Here we combine geophysical, stratigraphic and geotechnical measurements from the Lomonosov Ridge and Yermak Plateau and discuss the glacial geological implications of overconsolidated sediments. The degree of overconsolidation, determined from measurements of porosity and shear strength, is shown to result from consolidation and/or deformation below grounded ice and, with the exception of a single region on the Lomonosov Ridge, cannot be explained by erosion of overlying sediments. We demonstrate that the amount and depth of porosity loss associated with a middle Quaternary (∼790–950 thousand years ago – ka) grounding on the Yermak Plateau is compatible with sediment consolidation under an ice sheet or ice rise. Conversely, geotechnical properties of sediments from beneath late Quaternary ice-groundings in both regions, independently dated to Marine Isotope Stage (MIS) 6, indicate a more transient event commensurate with a passing tabular iceberg calved from an ice shelf.

Grounding line movement and ice shelf buttressing in marine ice sheets

Understanding the dynamics of marine ice sheets is integral to studying the evolution of the Antarctic ice sheet in both the short and long terms. An important component of the dynamics, grounding line migration, has proved difficult to represent in numerical models, and most successful attempts have made use of techniques that are only readily applicable to flow line models. However, to capture the stress transmission involved in another important component, the buttressing of a marine ice sheet by its ice shelf, the transverse direction must also be resolved. We introduce a model that solves the time‐dependent shelfy stream equations and makes use of mesh adaption techniques to overcome the difficulties typically associated with the numerics of grounding line migration. We compare the model output with a recent benchmark for flow line models and show that our model yields an accurate solution while using far less resources than would be required without mesh adaption. We also show that the mesh adapting techniques extend to two horizontal dimensions. Experiments are carried out to determine how both ice shelf buttressing and ice rises affect the marine instability predicted for an ice sheet on a foredeepened bed. We find that buttressing is not always sufficient to stabilize such a sheet but collapse of the grounded portion is still greatly delayed. We also find that the effect of an ice rise is similar to that of narrowing the ice shelf.

On the effects of anisotropic rheology on ice flow, internal structure, and the age‐depth relationship at ice divides

We use numerical modeling with a full‐system Stokes solver to elucidate the effects of nonlinear rheology and strain‐induced anisotropy on ice flow at ice divides. We find that anisotropic rheology profoundly affects the shape of both isochrone layering and surface topography. Anisotropic effects cause the formation of a downward curving fold, i.e., a syncline, in isochrones in the lower central area beneath the ice divide. When the resulting syncline is superimposed on the well‐known Raymond anticline, a double‐peaked Raymond bump is formed. Furthermore, to each side of the Raymond bump, flanking synclines are formed. In addition, anisotropic effects are found to give rise to a subtle concavity in the surface profile to both sides of the summit. The lower center syncline, the flanking synclines, and the near‐summit surface concavity have all previously been observed in nature, but hitherto no explanation for the genesis of these features has been given. We compare modeling results with radiograms collected from Fuchs Ice Piedmont and Kealey Ice Rise, Antarctica. Good overall agreement is found. In particular, we are able to reproduce all observed qualitative features of surface geometry and internal layering by including, and only by including, the effects of induced nonlinear rheological anisotropy on flow. Rheological anisotropy has the potential to profoundly affect the age distribution with depth, and caution must be exercised when estimating age of ice from ice cores with an isotropic model. The occurrence of linear features parallel to the ridge of ice divides, often seen in satellite imagery, is indicative of long‐term stability rather than signs of ongoing ice divide migration as previously suggested. Such ice divides are ideal locations for extracting ice cores.

Parameter sensitivity studies for the ice flow of the Ross Ice Shelf, Antarctica

The diagnostic, dynamic/thermodynamic ice shelf model Finite Element Shallow Shelf Approximation Code (FESSACODE) is applied to the Ross Ice Shelf. We simulate the present ice flow which results from the ice thickness distribution, the inflow at the grounding line, and the surface and bottom temperatures and compare results with measured flow velocities. Our reference simulation reproduces the general flow pattern and the magnitudes of the flow velocities reasonably well. The ice flow is found to be very sensitive to the flow enhancement factor, the ice thickness, and the ice temperature but robust against inflow velocities from ice streams, glaciers, and ice rises. The ice rises (Roosevelt Island and Crary Ice Rise) stabilize the ice shelf by significantly decreasing the flow velocities for the entire ice shelf area. The ice shelf is susceptible to global warming in that a 2°C surface warming entails an increase of the flow velocities by a factor 1.25, whereas a 10°C warming leads to an increase by a factor 3.1.

Changes in the ice plain of Whillans Ice Stream, West Antarctica

AbstractData from the mouth of the decelerating Whillans Ice Stream (WIS), West Antarctica, spanning 42 years are reviewed. Deceleration has continued, with local areas of both thinning and thickening occurring. The mean thinning rate is 0.48 ± 0.77 ma–1. No consistent overall pattern is observed. Ice thickens immediately upstream of Crary Ice Rise where deceleration and divergence are strongest, suggesting expanded upstream influence of the ice rise. Thinning is prevalent on the Ross Ice Shelf. Grounding-line advance at a rate of 0.3 km a–1 is detected in a few locations. Basal stresses vary across an ice-stream transect with a zone of enhanced flow at the margin. Marginal shear is felt at the ice-stream center. Mass-balance values are less negative, but larger errors of earlier measurements mask any possible temporal pattern. Comparisons of the recent flow field with flow stripes suggest WIS contributes less ice to the deep subglacial channel carved by Mercer Ice Stream and now flows straighter. The general lack of geometric changes suggests that the regional velocity decrease is due to changing basal conditions.

Surface mass balance of the Ward Hunt Ice Rise and Ward Hunt Ice Shelf, Ellesmere Island, Nunavut, Canada

The Ward Hunt Ice Rise and Ward Hunt Ice Shelf, located on Ellesmere Island, Canada, are two of the northernmost land ice masses on the North American continent. Surface mass balance measurements (excluding calving and subice processes) began in 1959 on the ice rise and in 1966 on the ice shelf but were frequently interrupted, most recently between 1986 and 2002. The surface balance of the ice rise and ice shelf follows the temporal pattern seen on other measured High Arctic glaciers. The overall surface mass losses over the last 45 years have been comparatively low (1.68 m water equivalent (w eq) for the ice rise and 3.1 m w eq for the ice shelf), which reflects their proximity to the Arctic Ocean. Nevertheless, the ice shelf appears to have weakened sufficiently in recent years to raise concerns about its possible disintegration in the near future. The 2002/2003 balance year was the most negative year on record (−0.33 m w eq for the ice rise and −0.54 m w eq for the ice shelf). Dynamical stresses related to wind, wave, and tidal action may further accelerate this process, as open water conditions on the Arctic Ocean become more prevalent. The Ward Hunt Ice Rise has so far remained in a reasonably healthy state in terms of its overall surface mass balance, although its long‐term survival is also threatened by current and predicted future climatic conditions.