No abstract

[1] Subglacial Lake Ellsworth has been proposed as a candidate for direct measurement and sampling, to identify microbial life and extract sedimentary climate records. We present a detailed characterization of the physiography of this subglacial lake from geophysical surveys, allowing bathymetry and geomorphic setting to be established. Lake Ellsworth is 14.7 km × 3.1 km with an area of 28.9 km 2 . Lake depth increases downlake from 52 m to 156 m, with a water body volume of 1.37 km 3 . The ice thickness suggests an unusual thermodynamic characteristic, with the critical pressure boundary intersecting the lake. Numerical modeling of water circulation has allowed accretion of basal ice to be estimated. We collate this physiographic and modeling information to confirm that Lake Ellsworth is ideal for direct access and propose an optimal drill site. The likelihood of dissolved gas exchange between the lake and the borehole is also assessed. Citation:

The West Antarctic Ice Sheet (WAIS) is considered the major contributor to global sea level rise in the Last Interglacial (LIG) and potentially in the future. Exposed fossil reef terraces suggest sea levels in excess of 7 m in the last warm era, of which probably not much more than 2 m are considered to originate from melting of the Greenland Ice Sheet. We simulate the evolution of the Antarctic Ice Sheet during the LIG with a 3‐D thermomechanical ice sheet model forced by an atmosphere‐ocean general circulation model (AOGCM). Our results show that high LIG sea levels cannot be reproduced with the atmosphere‐ocean forcing delivered by current AOGCMs. However, when taking reconstructed Southern Ocean temperature anomalies of several degrees, sensitivity studies indicate a Southern Ocean temperature anomaly threshold for total WAIS collapse of 2–3°C, accounting for a sea level rise of 3–4 m during the LIG. Potential future Antarctic Ice Sheet dynamics range from a moderate retreat to a complete collapse, depending on rate and amplitude of warming.

Abstract. Simulations of the ocean dynamics in the cavity under the Amery Ice Shelf, Antarctica, were carried out using a three-dimensional numerical ocean model. Two different boundary conditions were used to describe the open ocean barotropic exchange at the ice front. The simulations show that the circulation in the ocean cavity is predominantly barotropic and is generally steered by the cavity topography. The circulation is driven by the density gradient in the cavity, which is strongly influenced by the heat and salt fluxes from melting and freezing processes at the ice-ocean interface, and by the horizontal exchange of heat and salt across the open ocean boundary at the ice front. The interaction at the ice-ocean interface allows the basal component of the mass loss of the Amery Ice Shelf to be estimated.In the two simulations the computed losses were 5.8 Gt yr -1 and 18.0 Gt yr -1, values consistent with observations. The bulk of the melting occurred near the southern grounding line of the ice shelf, although substantial melting also occurred in areas where heat transport by horizontal circulation was large. Accretion was restricted to areas where water, from upstream melting, became supercooled as it ascended the ice shelf base.

No abstract

tion, while results in terms of transport patterns, melting rates at the ice shelf base, and hydrography are presented in the third section. We discuss these and other sensitivity experiments in light of available direct and indirect observations from the Filchner-Ronne Ice Shelf in section 4, which is followed by a summary and our conclusions in the last section. Model Description and Experimental Design Surface Boundary ConditionsOcean and ice shelf interact through heat and freshwater fluxes at the ice shelf base as a result of melting and freezing. At the boundary to the open ocean we apply conditions as described by Stevens [1991]. The internal velocity is computed from linearized momentum equations where unknown outside values enter only through the friction term. These outside values are obtained by extrapolation from the interior. The external mode must be specified in a suitable manner (see below).At outflow points the equations for temperature and salinity are simplified such that advection must be perpendicular to the boundary. The advection velocity v is corrected to allow the passage of waves through the boundary [Orlanski, 1976].

Lake Vostok, isolated from direct exchange with the atmosphere by about 4 km of ice for millions of years, provides a unique environment. This inaccessibility raises the importance of numerical models to investigate the physical conditions within the lake. Using a three-dimensional numerical model and the best available geometry, we test different parameter settings to define a standard model configuration suitable for studying flow in this subglacial lake. From our model runs we find a baroclinic circulation within the lake that splits into three different parts: Along a topographic ridge in the northern part of Lake Vostok, bottom water masses are transported eastward, diverging away from the ridge. In the lake's surface layer, the flow in these two vertical overturning cells has opposite directions. In the southern part of the lake, where freezing occurs across about 3,500 km 2 , two opposing gyres split the water column vertically. The general flow is stronger in the southern basin with horizontal velocities in the order of 1 mm/s. The strongest upwelling, found in the eastern part of this basin, is about 25 µm/s. We estimate the lower limit of the overturning timescale to be about 2.5 years vertically and 8.6 years horizontally. The basal mass loss of ice from the ice sheet floating on the lake is 5.6 mm/year (equivalent to a fresh water flux of 2.78 m 3 /s, or a basal ice loss of 0.09 km 3 /year). This imbalance indicates either a constant growth of the lake or its continuous (or periodical) discharge into a subglacial drainage system.

Abstract. Glaciers and ice caps exhibit currently the largest cryospheric contributions to sea level rise. Modelling the dynamics and mass balance of the major ice sheets is therefore an important issue to investigate the current state and the future response of the cryosphere in response to changing environmental conditions, namely global warming. This requires a powerful, easy-to-use, versatile multi-approximation ice dynamics model. Based on the well-known and established ice sheet model of Pattyn (2003) we develop the modular multi-approximation thermomechanic ice model RIMBAY, in which we improve the original version in several aspects like a shallow ice-shallow shelf coupler and a full 3D-groundingline migration scheme based on Schoof's (2007) heuristic analytical approach. We summarise the full Stokes equations and several approximations implemented within this model and we describe the different numerical discretisations. The results are cross-validated against previous publications dealing with ice modelling, and some additional artificial set-ups demonstrate the robustness of the different solvers and their internal coupling. RIMBAY is designed for an easy adaption to new scientific issues. Hence, we demonstrate in very different set-ups the applicability and functionality of RIMBAY in Earth system science in general and ice modelling in particular.