scholarly journals Implementing an empirical scalar constitutive relation for ice with flow-induced polycrystalline anisotropy in large-scale ice sheet models

2018 ◽  
Vol 12 (3) ◽  
pp. 1047-1067 ◽  
Author(s):  
Felicity S. Graham ◽  
Mathieu Morlighem ◽  
Roland C. Warner ◽  
Adam Treverrow

Abstract. The microstructure of polycrystalline ice evolves under prolonged deformation, leading to anisotropic patterns of crystal orientations. The response of this material to applied stresses is not adequately described by the ice flow relation most commonly used in large-scale ice sheet models – the Glen flow relation. We present a preliminary assessment of the implementation in the Ice Sheet System Model (ISSM) of a computationally efficient, empirical, scalar, constitutive relation which addresses the influence of the dynamically steady-state flow-compatible induced anisotropic crystal orientation patterns that develop when ice is subjected to the same stress regime for a prolonged period – sometimes termed tertiary flow. We call this the ESTAR flow relation. The effect on ice flow dynamics is investigated by comparing idealised simulations using ESTAR and Glen flow relations, where we include in the latter an overall flow enhancement factor. For an idealised embayed ice shelf, the Glen flow relation overestimates velocities by up to 17 % when using an enhancement factor equivalent to the maximum value prescribed in the ESTAR relation. Importantly, no single Glen enhancement factor can accurately capture the spatial variations in flow across the ice shelf generated by the ESTAR flow relation. For flow line studies of idealised grounded flow over varying topography or variable basal friction – both scenarios dominated at depth by bed-parallel shear – the differences between simulated velocities using ESTAR and Glen flow relations depend on the value of the enhancement factor used to calibrate the Glen flow relation. These results demonstrate the importance of describing the deformation of anisotropic ice in a physically realistic manner, and have implications for simulations of ice sheet evolution used to reconstruct paleo-ice sheet extent and predict future ice sheet contributions to sea level.

2017 ◽  
Author(s):  
Felicity S. Graham ◽  
Mathieu Morlighem ◽  
Roland C. Warner ◽  
Adam Treverrow

Abstract. The microstructural evolution that occurs in polycrystalline ice during deformation leads to the development of anisotropic rheological properties that are not adequately described by the most common, isotropic, ice flow relation used in large-scale ice sheet models – the Glen flow relation. We present a preliminary assessment of the implementation in the Ice Sheet System Model (ISSM) of a computationally-efficient, empirical, scalar, tertiary, anisotropic rheology (ESTAR). The effect of this anisotropic rheology on ice flow dynamics is investigated by comparing idealised simulations using ESTAR with those using the isotropic Glen flow relation, where the latter includes a flow enhancement factor. For an idealised embayed ice shelf, the Glen flow relation overestimates velocities by up to 17 % when using an enhancement factor equivalent to the maximum value prescribed by ESTAR. Importantly, no single Glen enhancement factor can accurately capture the spatial variations in flow over the ice shelf. For flow-line studies of idealised grounded flow over a bumpy topography or a sticky base – both scenarios dominated at depth by bed-parallel shear – the differences between simulated velocities using ESTAR and the Glen flow relation vary according to the value of the enhancement factor used to calibrate the Glen flow relation. These results demonstrate the importance of describing the anisotropic rheology of ice in a physically realistic manner, and have implications for simulations of ice sheet evolution used to reconstruct paleo-ice sheet extent and predict future ice sheet contributions to sea level.


2017 ◽  
Vol 11 (6) ◽  
pp. 2543-2554 ◽  
Author(s):  
Sainan Sun ◽  
Stephen L. Cornford ◽  
John C. Moore ◽  
Rupert Gladstone ◽  
Liyun Zhao

Abstract. Floating ice shelves exert a stabilizing force onto the inland ice sheet. However, this buttressing effect is diminished by the fracture process, which on large scales effectively softens the ice, accelerating its flow, increasing calving, and potentially leading to ice shelf breakup. We add a continuum damage model (CDM) to the BISICLES ice sheet model, which is intended to model the localized opening of crevasses under stress, the transport of those crevasses through the ice sheet, and the coupling between crevasse depth and the ice flow field and to carry out idealized numerical experiments examining the broad impact on large-scale ice sheet and shelf dynamics. In each case we see a complex pattern of damage evolve over time, with an eventual loss of buttressing approximately equivalent to halving the thickness of the ice shelf. We find that it is possible to achieve a similar ice flow pattern using a simple rule of thumb: introducing an enhancement factor ∼ 10 everywhere in the model domain. However, spatially varying damage (or equivalently, enhancement factor) fields set at the start of prognostic calculations to match velocity observations, as is widely done in ice sheet simulations, ought to evolve in time, or grounding line retreat can be slowed by an order of magnitude.


1979 ◽  
Vol 24 (90) ◽  
pp. 63-75 ◽  
Author(s):  
K. E. Rose

AbstractExtensive radio echo-sounding has mapped the part of West Antarctica between Byrd Station, the Whitmore Mountains, the Transantarctic Mountains, and the Ross Ice Shelf. The ice sheet in this area is dominated by five major sub-parallel ice streams (A–E), which are up to 100 km wide and extend inland from the grounding line of the Ross Ice Shelf for about 400 km. Their positions have been determined by crevassing seen on radio echo-sounding records, trimetrogon photographs, and Landsat imagery. The ice streams are characterized by their flat transverse cross-sections, while the intervening ice sheet exhibits domes and ridges. Ice flow lines are defined from the ice-surface contour pattern and the trend of the ice streams. It is apparent from this work that the flow line passing through Byrd Station joins ice stream D.The bedrock of the area is relatively smooth near the Ross Ice Shelf, becoming rougher near Byrd Station and especially so near the Whitmore Mountains. Bedrock troughs, which control the positions of the ice streams, are believed to have a tectonic origin.In this paper the role of the ice streams in the glaciological regime of West Antarctica is investigated from radio-echo data and estimates of balance velocity, basal shear stress, and basal temperatures.


1979 ◽  
Vol 24 (90) ◽  
pp. 63-75 ◽  
Author(s):  
K. E. Rose

AbstractExtensive radio echo-sounding has mapped the part of West Antarctica between Byrd Station, the Whitmore Mountains, the Transantarctic Mountains, and the Ross Ice Shelf. The ice sheet in this area is dominated by five major sub-parallel ice streams (A–E), which are up to 100 km wide and extend inland from the grounding line of the Ross Ice Shelf for about 400 km. Their positions have been determined by crevassing seen on radio echo-sounding records, trimetrogon photographs, and Landsat imagery. The ice streams are characterized by their flat transverse cross-sections, while the intervening ice sheet exhibits domes and ridges. Ice flow lines are defined from the ice-surface contour pattern and the trend of the ice streams. It is apparent from this work that the flow line passing through Byrd Station joins ice stream D.The bedrock of the area is relatively smooth near the Ross Ice Shelf, becoming rougher near Byrd Station and especially so near the Whitmore Mountains. Bedrock troughs, which control the positions of the ice streams, are believed to have a tectonic origin.In this paper the role of the ice streams in the glaciological regime of West Antarctica is investigated from radio-echo data and estimates of balance velocity, basal shear stress, and basal temperatures.


2014 ◽  
Vol 2 (2) ◽  
pp. 911-933 ◽  
Author(s):  
N. F. Glasser ◽  
S. J. A. Jennings ◽  
M. J. Hambrey ◽  
B. Hubbard

Abstract. Continent-wide mapping of longitudinal ice-surface structures on the Antarctic Ice Sheet reveals that they originate in the interior of the ice sheet and are arranged in arborescent networks fed by multiple tributaries. Longitudinal ice-surface structures can be traced continuously down-ice for distances of up to 1200 km. They are co-located with fast-flowing glaciers and ice streams that are dominated by basal sliding rates above tens of m yr-1 and are strongly guided by subglacial topography. Longitudinal ice-surface structures dominate regions of converging flow, where ice flow is subject to non-coaxial strain and simple shear. Associating these structures with the AIS' surface velocity field reveals (i) ice residence times of ~ 2500 to 18 500 years, and (ii) undeformed flow-line sets for all major flow units analysed except the Kamb Ice Stream and the Institute and Möller Ice Stream areas. Although it is unclear how long it takes for these features to form and decay, we infer that the major ice-flow and ice-velocity configuration of the ice sheet may have remained largely unchanged for several thousand years, and possibly even since the end of the last glacial cycle. This conclusion has implications for our understanding of the long-term landscape evolution of Antarctica, including large-scale patterns of glacial erosion and deposition.


2020 ◽  
Vol 14 (1) ◽  
pp. 17-37 ◽  
Author(s):  
Sebastian H. R. Rosier ◽  
G. Hilmar Gudmundsson

Abstract. An extensive network of GPS sites on the Filchner–Ronne Ice Shelf and adjoining ice streams shows strong tidal modulation of horizontal ice flow at a range of frequencies. A particularly strong (horizontal) response is found at the fortnightly (Msf) frequency. Since this tidal constituent is absent in the (vertical) tidal forcing, this observation implies the action of some non-linear mechanism. Another striking aspect is the strong amplitude of the flow perturbation, causing a periodic reversal in the direction of ice shelf flow in some areas and a 10 %–20 % change in speed at grounding lines. No model has yet been able to reproduce the quantitative aspects of the observed tidal modulation across the entire Filchner–Ronne Ice Shelf. The cause of the tidal ice flow response has, therefore, remained an enigma, indicating a serious limitation in our current understanding of the mechanics of large-scale ice flow. A further limitation of previous studies is that they have all focused on isolated regions and interactions between different areas have, therefore, not been fully accounted for. Here, we conduct the first large-scale ice flow modelling study to explore these processes using a viscoelastic rheology and realistic geometry of the entire Filchner–Ronne Ice Shelf, where the best observations of tidal response are available. We evaluate all relevant mechanisms that have hitherto been put forward to explain how tides might affect ice shelf flow and compare our results with observational data. We conclude that, while some are able to generate the correct general qualitative aspects of the tidally induced perturbations in ice flow, most of these mechanisms must be ruled out as being the primary cause of the observed long-period response. We find that only tidally induced lateral migration of grounding lines can generate a sufficiently strong long-period Msf response on the ice shelf to match observations. Furthermore, we show that the observed horizontal short-period semidiurnal tidal motion, causing twice-daily flow reversals at the ice front, can be generated through a purely elastic response to basin-wide tidal perturbations in the ice shelf slope. This model also allows us to quantify the effect of tides on mean ice flow and we find that the Filchner–Ronne Ice Shelf flows, on average, ∼ 21 % faster than it would in the absence of large ocean tides.


1979 ◽  
Vol 24 (90) ◽  
pp. 89-101 ◽  
Author(s):  
W. F. Budd ◽  
N. W. Young

AbstractIn order to determine accurate velocities of the ice sheet in the interior of Antarctica, approximately along a flow line, a detailed trilateration net was established in 1973 from the summit of Law Dome (100 km inland) to about 250 km south near the 2000 m contour. The net consisted of a double line of markers approximately 10 km apart with all sides and diagonals of the quadrilaterals measured with telluro-meters. In addition, satellite doppler survey positions and astronomical azimuths were determined at about 50 km intervals to control the net on the large scale. Other measurements carried out en route included: continuous barometric levelling, radio echo-sounding, gravimetry, accumulation, and surface sampling. The route was close to an earlier traverse route which reached Vostok in 1962 and along which other data, including snow-surface temperatures and temperature–depth gradients, were determined.The trilateration net was re-surveyed in 1975 allowing velocities and strain-rates to be determined. The results indicate that the ice sheet is close to balance in this region.Therefore, the measured velocities were used together with “balance velocities”, further inland, to carry out a modelling study of a flow line, to derive particle trajectories, ages, temperature profiles, and “dynamics velocities”, from a flow law. The results provide further insight into the dynamics and flow properties of the ice sheet.


2020 ◽  
Vol 12 (8) ◽  
pp. 1327 ◽  
Author(s):  
Anna Ruth W. Halberstadt ◽  
Colin J. Gleason ◽  
Mahsa S. Moussavi ◽  
Allen Pope ◽  
Luke D. Trusel ◽  
...  

Surface meltwater generated on ice shelves fringing the Antarctic Ice Sheet can drive ice-shelf collapse, leading to ice sheet mass loss and contributing to global sea level rise. A quantitative assessment of supraglacial lake evolution is required to understand the influence of Antarctic surface meltwater on ice-sheet and ice-shelf stability. Cloud computing platforms have made the required remote sensing analysis computationally trivial, yet a careful evaluation of image processing techniques for pan-Antarctic lake mapping has yet to be performed. This work paves the way for automating lake identification at a continental scale throughout the satellite observational record via a thorough methodological analysis. We deploy a suite of different trained supervised classifiers to map and quantify supraglacial lake areas from multispectral Landsat-8 scenes, using training data generated via manual interpretation of the results from k-means clustering. Best results are obtained using training datasets that comprise spectrally diverse unsupervised clusters from multiple regions and that include rock and cloud shadow classes. We successfully apply our trained supervised classifiers across two ice shelves with different supraglacial lake characteristics above a threshold sun elevation of 20°, achieving classification accuracies of over 90% when compared to manually generated validation datasets. The application of our trained classifiers produces a seasonal pattern of lake evolution. Cloud shadowed areas hinder large-scale application of our classifiers, as in previous work. Our results show that caution is required before deploying ‘off the shelf’ algorithms for lake mapping in Antarctica, and suggest that careful scrutiny of training data and desired output classes is essential for accurate results. Our supervised classification technique provides an alternative and independent method of lake identification to inform the development of a continent-wide supraglacial lake mapping product.


2019 ◽  
Vol 13 (10) ◽  
pp. 2771-2787 ◽  
Author(s):  
Jan De Rydt ◽  
Gudmundur Hilmar Gudmundsson ◽  
Thomas Nagler ◽  
Jan Wuite

Abstract. Despite the potentially detrimental impact of large-scale calving events on the geometry and ice flow of the Antarctic Ice Sheet, little is known about the processes that drive rift formation prior to calving, or what controls the timing of these events. The Brunt Ice Shelf in East Antarctica presents a rare natural laboratory to study these processes, following the recent formation of two rifts, each now exceeding 50 km in length. Here we use 2 decades of in situ and remote sensing observations, together with numerical modelling, to reveal how slow changes in ice shelf geometry over time caused build-up of mechanical tension far upstream of the ice front, and culminated in rift formation and a significant speed-up of the ice shelf. These internal feedbacks, whereby ice shelves generate the very conditions that lead to their own (partial) disintegration, are currently missing from ice flow models, which severely limits their ability to accurately predict future sea level rise.


1982 ◽  
Vol 3 ◽  
pp. 346 ◽  
Author(s):  
N.W. Young ◽  
D. SheehY ◽  
T. Hamley

Trilateration and single line surveys have been made to about 900 km inland of Casey, Wilkes Land, to measure surface elevation, ice thickness, horizontal velocity, and other parameters. On the large scale the velocity U increases smoothly from 8 m a−1, 800 km inland, to 280 m a−1 inland of the fast outlet streams. This increase in velocity is associated with a corresponding increase in the large-scale smoothed (over about 30 ice thicknesses) basal shear stress τb from 0.4 to 1.5 bar. The mean shear strain-rate through the ice sheet U/Z = kτb4 , where Z is the ice thickness (range 4 500 to 1 700 m). At scales of one to several ice thicknesses large variations occur in surface slope and ice thickness without proportionally large velocity variations, because of the effect of the longitudinal stress. Detailed measurements made over a 30 km section indicated that the surface longitudinal strain-rate gradient varied from -1.7 to +1.3×l0−6 a−1 m−1 along with variations in surface slope of from -3.5 to +1.5%. A multilayer model, based on the solution of the biharmonic equation for the stream function, was used in a study of the ice flow associated with these surface undulations. Given the bedrock topography and large-scale flow parameters, the model closely predicted the measured surface profile when the variation of the surface accumulation rate over an undulation was also considered.


Sign in / Sign up

Export Citation Format

Share Document