Victoria Lee

Incorporation of particulates into accreted ice above subglacial Vostok lake, Antarctica

Abstract The nature of microscopic particulates in meteoric and accreted ice from the Vostok (Antarctica) ice core is assessed in conjunction with existing ice-core data to investigate the mechanism by which particulates are incorporated into refrozen lake water. Melted ice samples from a range of ice-core depths were filtered through 0.2 μm polycarbonate membranes, and secondary electron images were collected at ×500 magnification using a scanning electron microscope. Image analysis software was used to characterize the size and shape of particulates. Similar distributions of major-axis lengths, surface areas and shape factors (aspect ratio and compactness) for particulates in all accreted ice samples suggest that a single process may be responsible for incorporating the vast majority of particulates for all depths. Calculation of Stokes settling velocities for particulates of various sizes implies that 98% of particulates observed could ‘float’ to the ice–water interface with upward water velocities of 0.0003 ms–1 where they could be incorporated by growing ice crystals, or by rising frazil ice crystals. The presence of particulates that are expected to sink in the water column (2%) and the uneven distribution of particulates in the ice core further implies that periodic perturbations to the lake’s circulation, involving increased velocities, may have occurred in the past.

Adaptive mesh refinement versus subgrid friction interpolation in simulations of Antarctic ice dynamics

ABSTRACT At least in conventional hydrostatic ice-sheet models, the numerical error associated with grounding line dynamics can be reduced by modifications to the discretization scheme. These involve altering the integration formulae for the basal traction and/or driving stress close to the grounding line and exhibit lower – if still first-order – error in the MISMIP3d experiments. MISMIP3d may not represent the variety of real ice streams, in that it lacks strong lateral stresses, and imposes a large basal traction at the grounding line. We study resolution sensitivity in the context of extreme forcing simulations of the entire Antarctic ice sheet, using the BISICLES adaptive mesh ice-sheet model with two schemes: the original treatment, and a scheme, which modifies the discretization of the basal traction. The second scheme does indeed improve accuracy – by around a factor of two – for a given mesh spacing, but

Initialization of an ice-sheet model for present-day Greenland

\lesssim 1

Is Vostok lake in steady state

km resolution is still necessary. For example, in coarser resolution simulations Thwaites Glacier retreats so slowly that other ice streams divert its trunk. In contrast, with

Numerical simulation of three-dimensional velocity fields in pressurized and non-pressurized Nye channels

\lesssim 1

Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP

km meshes, the same glacier retreats far more quickly and triggers the final phase of West Antarctic collapse a century before any such diversion can take place.

Calibrated prediction of Pine Island Glacier retreat during the 21st and 22nd centuries with a coupled flowline model

Abstract We construct initial conditions for an ice flow model of the Greenland ice sheet (GrIS). GrIS has been losing mass at an increasing rate over the past two decades, and a significant proportion of this loss is due to dynamic thinning of narrow outlet glaciers. We solve an inverse problem to estimate poorly known basal and englacial parameters given observed geometry and surface velocities. A weighted cost function, resolved to 4 km in the interior of the ice sheet and 1 km in regions of fast-flowing ice at the margin, is minimized to find two-dimensional fields for a stiffness factor, which is a coefficient of the effective viscosity, and basal traction coefficient. Using these fields, we run the model under present-day climate to damp large-amplitude, short-wavelength fluctuations in the flux divergence. The time-dependent model uses an adaptive mesh with resolution ranging from 8 km of the base grid to 500 m in areas of fast-flowing ice to capture the behaviour of the main outlet glaciers. The ice discharge calculated from the initial conditions for GrIS and individual glaciers compares well with values calculated from observations.

Grounding line migration in an adaptive mesh ice sheet model

Abstract Stable-isotope (δD and δ18O) data from the Vostok (East Antarctica) ice core are used to explore whether or not subglacial Vostok lake is in isotopic steady state. A simple box model shows that the lake is likely to be in steady state on time-scales of the order of 104–105 years (three to four residence times of the water in the lake), given our current knowledge of north–south and east–west gradients in the stable-isotopic composition of precipitation in the vicinity of Vostok station and Ridge B. However, the lake may not be in perfect steady state depending on the precise location of the melting area, which determines the source region of inflowing ice, and on the magnitude of the east–west gradient in isotopic compositions in the vicinity of Vostok station and Ridge B.

Design and results of the ice sheet model initialisation initMIP-Greenland : An ISMIP6 intercomparison

Abstract Channels incised into bedrock, or Nye channels, often form an important component of subglacial drainage at temperate glaciers, and their structure exerts control over patterns and rates of (a) channel erosion, (b) water flow-velocity and (c) water pressure. The latter, in turn, exerts a strong control over basal traction and, thus, ice dynamics. In order to investigate these controls, it is necessary to quantify detailed flow processes in subglacial Nye channels. However, it is effectively impossible to acquire such measurements from fully pressurized, subglacial channels. To solve this problem, we here apply a three-dimensional, finite-volume solution of the Reynolds averaged Navier– Stokes (RANS) equations with a one-equation mixing-length turbulence closure to simulate flow in a 3 m long section of an active Nye channel located in the immediate foreground of Glacier de Tsanfleuron, Switzerland. Numerical model output permits high-resolution visualization of water flow through the channel reach, and enables evaluation of the experimental manipulation of the pressure field adopted across the overlying ice lid. This yields an increased theoretical understanding of the hydraulic behaviour of Nye channels, and, in the future, of their effect on glacier drainage, geomorphology and ice dynamics.