An official website of the United States government
Abstract We develop a two-dimensional, plan-view formulation of ice-shelf flow and viscoelastic ice-shelf flexure. This formulation combines, for the first time, the shallow-shelf approximation for horizontal ice-shelf flow (and shallow-stream approximation for flow on lubricated beds such as where ice rises and rumples form), with the treatment of a thin-plate flexure. We demonstrate the treatment by performing two finite-element simulations: one of the relict pedestalled lake features that exist on some debris-covered ice shelves due to strong heterogeneity in surface ablation, and the other of ice rumpling in the grounding zone of an ice rise. The proposed treatment opens new venues to investigate physical processes that require coupling between the longitudinal deformation and vertical flexure, for instance, the effects of surface melting and supraglacial lakes on ice shelves, interactions with the sea swell, and many others.
more »
« less
A new vertically integrated MOno-Layer Higher-Order (MOLHO) ice flow model
Abstract. Numerical simulations of ice sheets rely on the momentum balance to determine how ice velocities change as the geometry of the system evolves. Ice is generally assumed to follow a Stokes flow with a nonlinear viscosity. Several approximations have been proposed in order to lower the computational cost of a full-Stokes stress balance. A popular option is the Blatter–Pattyn or higher-order model (HO), which consists of a three-dimensional set of equations that solves the horizontal velocities only. However, it still remains computationally expensive for long transient simulations. Here we present a depth-integrated formulation of the HO model, which can be solved on a two-dimensional mesh in the horizontal plane. We employ a specific polynomial function to describe the vertical variation in the velocity, which allows us to integrate the vertical dimension using a semi-analytic integration. We assess the performance of this MOno-Layer Higher-Order (MOLHO) model to compute ice velocities and simulate grounding line dynamics on standard benchmarks (ISMIP-HOM and MISMIP3D). We compare MOLHO results to the ones obtained with the original three-dimensional HO model. We also compare the time performance of both models in time-dependent runs. Our results show that the ice velocities and grounding line positions obtained with MOLHO are in very good agreement with the ones from HO. In terms of computing time, MOLHO requires less than 10 % of the computational time of a typical HO model, for the same simulations. These results suggest that the MOno-Layer Higher-Order formulation provides improved computational time performance and a comparable accuracy compared to the HO formulation, which opens the door to higher-order paleo simulations.
more »
« less
- Award ID(s):
- 2152622
- PAR ID:
- 10319979
- Date Published:
- Journal Name:
- The Cryosphere
- Volume:
- 16
- Issue:
- 1
- ISSN:
- 1994-0424
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)The use of a boundary-layer parameterization of buttressing and ice flux across grounding lines in a two dimensional ice-sheet model is improved by allowing general orientations of the grounding line. This and another modification to the model’s grounding-line parameterization are assessed in three settings: rectangular fjord-like domains – the third Marine Ice Sheet Model Intercomparison Project (MISMIPC) and Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d) – and future simulations of West Antarctic ice retreat under Representative Concentration Pathway (RCP) 8.5-based climates. The new modifications are found to have significant effects on the fjord-like results, which are now within the envelopes of other models in the MISMIP+ and MISMIP3d intercomparisons. In contrast, the modifications have little effect on West Antarctic retreat, presumably because dynamics in the wider major Antarctic basins are adequately represented by the model’s previous simpler one-dimensional formulation. As future grounding lines retreat across very deep bedrock topography in the West Antarctic simulations, buttressing is weak and deviatoric stress measures exceed the ice yield stress, implying that structural failure at these grounding lines would occur. We suggest that these grounding-line quantities should be examined in similar projections by other ice models to better assess the potential for future structural failure.more » « less
-
Abstract Delineating the grounding line of marine-terminating glaciers—where ice starts to become afloat in ocean waters—is crucial for measuring and understanding ice sheet mass balance, glacier dynamics, and their contributions to sea level rise. This task has been previously done using time-consuming, mostly-manual digitizations of differential interferometric synthetic-aperture radar interferograms by human experts. This approach is no longer viable with a fast-growing set of satellite observations and the need to establish time series over entire continents with quantified uncertainties. We present a fully-convolutional neural network with parallel atrous convolutional layers and asymmetric encoder/decoder components that automatically delineates grounding lines at a large scale, efficiently, and accompanied by uncertainty estimates. Our procedure detects grounding lines within 232 m in 100-m posting interferograms, which is comparable to the performance achieved by human experts. We also find value in the machine learning approach in situations that even challenge human experts. We use this approach to map the tidal-induced variability in grounding line position around Antarctica in 22,935 interferograms from year 2018. Along the Getz Ice Shelf, in West Antarctica, we demonstrate that grounding zones are one order magnitude (13.3 ± 3.9) wider than expected from hydrostatic equilibrium, which justifies the need to map grounding lines repeatedly and comprehensively to inform numerical models.more » « less
-
Abstract Basal crevasses threaten the stability of ice shelves through the potential to form rifts and calve icebergs. Furthermore, it is important to determine the dependence of crevasse stability on temperature due to large vertical temperature variations on ice shelves. In this work, considering the vertical temperature profile through ice viscosity, we compare (1) the theoretical crack depths and (2) the threshold stress causing the transition from basal crevasses to full thickness fractures in several fracture theories. In the Zero Stress approximation, the depth-integrated force at the crevassed and non-crevassed location are unbalanced, violating the volume-integrated Stokes equation. By incorporating a Horizontal Force Balance (HFB) argument, recent work showed analytically that the threshold stress for rift initiation is only half of that predicted by the Zero Stress approximation. We generalize the HFB theory to show that while the temperature profile influences crack depths, the threshold rifting stress is insensitive to temperature. We compare with observations and find that HFB best matches observed rifts. Using HFB instead of Zero Stress for cracks in an ice-sheet model would substantially enlarge the predicted fracture depth, reduce the threshold rifting stress and potentially increase the projected rate of ice shelf mass loss.more » « less
-
Abstract. Time-dependent simulations of ice sheets require two equations to be solved:the mass transport equation, derived from the conservation of mass, and thestress balance equation, derived from the conservation of momentum. The masstransport equation controls the advection of ice from the interior of the icesheet towards its periphery, thereby changing its geometry. Because it isbased on an advection equation, a stabilization scheme needs to beemployed when solved using the finite-element method. Several stabilizationschemes exist in the finite-element method framework, but their respectiveaccuracy and robustness have not yet been systematically assessed forglaciological applications. Here, we compare classical schemes used in thecontext of the finite-element method: (i) artificial diffusion, (ii)streamline upwinding, (iii) streamline upwind Petrov–Galerkin, (iv)discontinuous Galerkin, and (v) flux-corrected transport. We also look at thestress balance equation, which is responsible for computing the ice velocitythat “advects” the ice downstream. To improve the velocity computationaccuracy, the ice-sheet modeling community employs several sub-elementparameterizations of physical processes at the grounding line, the point wherethe grounded ice starts to float onto the ocean. Here, we introduce a newsub-element parameterization for the driving stress, the force that drives theice-sheet flow. We analyze the response of each stabilization scheme byrunning transient simulations forced by ice-shelf basal melt. The simulationsare based on an idealized ice-sheet geometry for which there is no influenceof bedrock topography. We also perform transient simulations of the AmundsenSea Embayment, West Antarctica, where real bedrock and surface elevations areemployed. In both idealized and real ice-sheet experiments, stabilizationschemes based on artificial diffusion lead systematically to a bias towardsmore mass loss in comparison to the other schemes and therefore should beavoided or employed with a sufficiently high mesh resolution in the vicinityof the grounding line. We also run diagnostic simulations to assess theaccuracy of the driving stress parameterization, which, in combination with anadequate parameterization for basal stress, provides improved numericalconvergence in ice speed computations and more accurate results.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/tc-16-179-2022
