Search

One of the last remaining floating tongues of the Greenland ice sheet (GrIS), the Petermann Glacier ice shelf (PGIS), is seasonally shielded from warm Atlantic water (AW) by the formation of sea ice arches in the Nares Strait. However, continued decline of the Arctic sea ice extent and thickness suggests that arch formation is likely to become anomalous, necessitating an investigation into the response of PGIS to a year-round mobile and thin sea ice cover. We use a high-resolution unstructured grid 3-D ocean–sea ice–ice shelf setup, featuring an improved sub-ice-shelf bathymetry and a realistic PGIS geometry, to investigate in unprecedented detail the implications of transitions in the Nares Strait sea ice regime, that is, from a thick and landfast sea ice regime to a mobile, and further, a thin and mobile sea ice regime, with regard to PGIS basal melt. In all three sea ice regimes, basal melt near the grounding line (GL) presents a seasonal increase during summer, driven by a higher thermal driving. The stronger melt overturning increases the friction velocity slightly downstream, where enhanced friction-driven turbulent mixing further increases the thermal driving, substantially increasing the local melt. As the sea ice cover becomes mobile and thin, wind and (additionally in winter) convectively upwelled AW from the Nares Strait enter the PGIS cavity. While its effect on basal melting is largely limited to the shallower (<200 m) drafts during winter, in summer it extends to the GL (ca. 600 m) depth. In the absence of an increase in thermal driving, increased melting under the deeper (>200 m) drafts in winter is solely driven by the increased vertical shear of a more energetic boundary layer current. A similar behaviour is noted when transitioning from a mobile to a thin mobile sea ice cover in summer, when increases in thermal driving are negligible and increases in melt are congruent with increases in friction velocity. These results suggest that the projected continuation of the warming of the Arctic Ocean until the end of the 21st century and the accompanying decline in the Arctic sea ice extent and thickness will amplify the basal melt of PGIS, impacting the long-term stability of the Petermann Glacier and its contribution to the future GrIS mass loss and sea level rise. [ABSTRACT FROM AUTHOR]

Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve due to the nonlinearity introduced by the rheology. Therefore, approximations to the FS equations are commonly used, especially when modeling a marine ice sheet (ice sheet, ice shelf, and/or ice stream) for 103 years or longer. The shallow ice approximation (SIA) and shallow shelf approximation (SSA) are commonly used but are accurate only for certain parts of an ice sheet. Here, we report a novel way of iteratively coupling FS and SSA that has been implemented in Elmer/Ice and applied to conceptual marine ice sheets. The FS–SSA coupling appears to be very accurate; the relative error in velocity compared to FS is below 0.5 % for diagnostic runs and below 5 % for prognostic runs. Results for grounding line dynamics obtained with the FS–SSA coupling are similar to those obtained from an FS model in an experiment with a periodical temperature forcing over 3000 years that induces grounding line advance and retreat. The rapid convergence of the FS–SSA coupling shows a large potential for reducing computation time, such that modeling a marine ice sheet for thousands of years should become feasible in the near future. Despite inefficient matrix assembly in the current implementation, computation time is reduced by 32 %, when the coupling is applied to a 3-D ice shelf. [ABSTRACT FROM AUTHOR]

Recent geological and geophysical data suggest that a 1 km thick ice shelf extended over the glacial Arctic Ocean during Marine Isotope Stage 6, about 140 000 years ago. Here, we theoretically analyse the development and equilibrium features of such an ice shelf, using scaling analyses and a one-dimensional ice-sheet-ice-shelf model. We find that the dynamically most consistent scenario is an ice shelf with a nearly uniform thickness that covers the entire Arctic Ocean. Further, the ice shelf has two regions with distinctly different dynamics: a vast interior region covering the central Arctic Ocean and an exit region towards the Fram Strait. In the interior region, which is effectively dammed by the Fram Strait constriction, there are strong back stresses and the mean ice-shelf thickness is controlled primarily by the horizontally integrated mass balance. A narrow transition zone is found near the continental grounding line, in which the ice-shelf thickness decreases offshore and approaches the mean basin thickness. If the surface accumulation and mass flow from the continental ice masses are sufficiently large, the ice-shelf thickness grows to the point where the ice shelf grounds on the Lomonosov Ridge. As this occurs, the back stress increases in the Amerasian Basin and the ice-shelf thickness becomes larger there than in the Eurasian Basin towards the Fram Strait. Using a one-dimensional ice-dynamic model, the stability of equilibrium ice-shelf configurations without and with grounding on the Lomonosov Ridge are examined. We find that the grounded ice-shelf configuration should be stable if the two Lomonosov Ridge grounding lines are located on the opposites sides of the ridge crest, implying that the downstream grounding line is located on a downward sloping bed. This result shares similarities with the classical result on marine ice-sheet stability of Weertman, but due to interactions between the Amerasian and Eurasian ice-shelf segments the mass flux at the downstream grounding line decreases rather than increases with ice thickness. [ABSTRACT FROM AUTHOR]

Abstract: The simulation of dynamically coupled ice sheet, ice stream, and ice shelf-systems poses a challenge to most numerical ice sheet models. Here we review present ice sheet model limitations targeting a broader audience within Earth Sciences, also those with no specific background in numerical modeling, in order to facilitate cross-disciplinary communication between especially paleoglaciologists, marine and terrestrial geologists, and numerical modelers. The ‘zero order’ (Shallow Ice Approximation, SIA)-, ‘higher order’-, and ‘full Stokes’ ice sheet models are described conceptually and complemented by an outline of their derivations. We demonstrate that higher order models are required to simulate coupled ice sheet-ice shelf and ice sheet-ice stream systems, in particular if the results are aimed to complement spatial ice flow reconstructions based on higher resolution geological and geophysical data. The zero order SIA model limitations in capturing ice stream behavior are here illustrated by conceptual simulations of a glaciation on Svalbard. The limitations are obvious from the equations comprising a zero order model. However, under certain circumstances, simulation results may falsely give the impression that ice streams indeed are simulated with a zero order SIA model. [Copyright &y& Elsevier]

Abstract: While there are numerous hypotheses concerning glacial–interglacial environmental and climatic regime shifts in the Arctic Ocean, a holistic view on the Northern Hemisphere's late Quaternary ice-sheet extent and their impact on ocean and sea-ice dynamics remains to be established. Here we aim to provide a step in this direction by presenting an overview of Arctic Ocean glacial history, based on the present state-of-the-art knowledge gained from field work and chronological studies, and with a specific focus on ice-sheet extent and environmental conditions during the Last Glacial Maximum (LGM). The maximum Quaternary extension of ice sheets is discussed and compared to LGM. We bring together recent results from the circum-Arctic continental margins and the deep central basin; extent of ice sheets and ice streams bordering the Arctic Ocean as well as evidence for ice shelves extending into the central deep basin. Discrepancies between new results and published LGM ice-sheet reconstructions in the high Arctic are highlighted and outstanding questions are identified. Finally, we address the ability to simulate the Arctic Ocean ice sheet complexes and their dynamics, including ice streams and ice shelves, using presently available ice-sheet models. Our review shows that while we are able to firmly reject some of the earlier hypotheses formulated to describe Arctic Ocean glacial conditions, we still lack information from key areas to compile the holistic Arctic Ocean glacial history. [Copyright &y& Elsevier]

The catastrophic break-ups of the floating Larsen A and B ice shelves (Antarctica) in 1995 and 2002 and associated acceleration of glaciers that flowed into these ice shelves were among the most dramatic glaciological events observed in historical time. This raises a question about the larger West Antarctic ice shelves. Do these shelves, with their much greater glacial discharge, have a history of collapse? Here we describe features from the seafloor in Pine Island Bay, West Antarctica, which we interpret as having been formed during a massive ice shelf break-up and associated grounding line retreat. This evidence exists in the form of seafloor landforms that we argue were produced daily as a consequence of tidally influenced motion of mega-icebergs maintained upright in an iceberg armada produced from the disintegrating ice shelf and retreating grounding line. The break-up occurred prior to ca. 12 ka and was likely a response to rapid sea-level rise or ocean warming at that time. [ABSTRACT FROM AUTHOR]