Alun Hubbard

Deglacial history of the West Antarctic Ice Sheet in the Weddell Sea embayment: Constraints on past ice volume change

The retreat history of the West Antarctic Ice Sheet (WAIS) since the Last Glacial Maximum is important for understanding the process of rapid deglaciation, constraining models that seek to predict the future trajectory of the ice sheet, and for estimating rates of sea-level change. Here we report new glacial geologic data from the southwestern Weddell Sea embayment that demonstrate that this part of the WAIS was thinner than previously suggested, and that there was progressive thinning of the ice sheet by 230–480 m since ca. 15 ka. We use geomorphological data and a numerical ice sheet model to reconstruct the ice sheet in the Weddell Sea at the Last Glacial Maximum. The volume of this ice would have added between 1.4 and 2.0 m to postglacial sea-level rise and would not have been sufficient to contribute significantly to meltwater pulse 1A, a rapid rise in sea level ∼14,200 yr ago.

Hydrological controls on patterns of surface, internal and basal motion during three "spring events" : Haut Glacier d'Arolla, Switzerland

Three early-melt-season high-velocity events (or “spring events”) occurred on Haut Glacier d’Arolla, Switzerland, during the melt seasons of 1998 and 1999. The events involve enhanced glacier velocity during periods of rapidly increasing bulk discharge in the proglacial stream and high subglacial water pressures. However, differences in spatial patterns of surface velocity, internal ice deformation rates, the spatial extent of high subglacial water pressures and in rates of subglacial sediment deformation suggest different hydrological and mechanical controls. The data from two of the events suggest widespread ice–bed decoupling, particularly along a subglacial drainage axis creating the highest rates of basal motion and “plug flow” in the overlying ice. The other event showed evidence of less extensive ice–bed decoupling and sliding along the drainage axis with more mechanical support for ice overburden transferred to areas adjacent to decoupled areas. We suggest that: (1) plug flow may be a common feature on glaciers experiencing locally induced reductions in basal drag; (2) under certain circumstances, enhanced surface motion may be due in part to non-locally forced enhanced bed deformation; and (3) subglacial sediment deformation is confined to a depth of the order of centimetres to decimetres.

Persistent flow acceleration within the interior of the Greenland ice sheet

We present surface velocity measurements from a high-elevation site located 140 km from the western margin of the Greenland ice sheet, and ~ 50 km into its accumulation area. Annual velocity increased each year from 51.78 ± 0.01 m yr−1 in 2009 to 52.92 ± 0.01 m yr−1 in 2012—a net increase of 2.2%. These data also reveal a strong seasonal velocity cycle of up to 8.1% above the winter mean, driven by seasonal melt and supraglacial lake drainage. Sole et al. (2013) recently argued that ice motion in the ablation area is mediated by reduced winter flow following the development of efficient subglacial drainage during warmer, faster, summers. Our data extend this analysis and reveal a year-on-year increase in annual velocity above the equilibrium line altitude, where despite surface melt increasing, it is still sufficiently low to hinder the development of efficient drainage under thick ice.

A modelling reconstruction of the last glacial maximum ice sheet and its deglaciation in the vicinity of the Northern Patagonian Icefield, South America

ABSTRACT. A time‐dependent model is used to investigate the interaction between climate, extent and fluctuations of Patagonian ice sheet between 45° and 48°S during the last glacial maximum (LGM) and its subsequent deglaciation. The model is applied at 2 km resolution and enables ice thickness, lithospheric response and ice deformation and sliding to interact freely and is perturbed from present day by relative changes in sea level and equilibrium line altitude (ELA). Experiments implemented to identify an LGM configuration compatible with the available empirical record, indicate that a stepped ELA lowering of 750 to 950 m is required over 15000 years to bracket the Fenix I‐V suite of moraines at Lago Buenos Aires. However, 900 m of ELA lowering yields an ice sheet which best matches the Fenix V moraine (c. 23000 a BP) and Caldenius’ reconstructed LGM limit for the entire modelled area. This optimum LGM experiment yields a highly dynamic, low aspect ice sheet, with a mean ice thickness of c. 1130 m drained by numerous large ice streams to the western, seaward margin and two large, fast‐flowing outlet lobes to the east. Forcing this scenario into deglaciation using a re‐scaled Vostok ice core record results in an ice sheet that slowly shrinks by 25% to c. 14500 a bp, after which it experiences a rapid collapse, loosing some 85% of its volume in c. 800 years. Its margins stabilize during the Antarctic Cold Reversal after which it shrinks to near present‐day limits by 11 000 a bp.

Sensitive response of the Greenland Ice Sheet to surface melt drainage over a soft bed

The dynamic response of the Greenland Ice Sheet (GrIS) depends on feedbacks between surface meltwater delivery to the subglacial environment and ice flow. Recent work has highlighted an important role of hydrological processes in regulating the ice flow, but models have so far overlooked the mechanical effect of soft basal sediment. Here we use a three-dimensional model to investigate hydrological controls on a GrIS soft-bedded region. Our results demonstrate that weakening and strengthening of subglacial sediment, associated with the seasonal delivery of surface meltwater to the bed, modulates ice flow consistent with observations. We propose that sedimentary control on ice flow is a viable alternative to existing models of evolving hydrological systems, and find a strong link between the annual flow stability, and the frequency of high meltwater discharge events. Consequently, the observed GrIS resilience to enhanced melt could be compromised if runoff variability increases further with future climate warming.

Glacier mass-balance determination by remote sensing and high-resolution modelling

An indirect methodology for determining the distribution of mass balance at high spatial resolution using remote sensing and ice-flow modelling is presented. The method, based on the mass-continuity equation, requires two datasets collected over the desired monitoring interval: (i) the spatial pattern of glacier surface-elevation change, and (ii) the mass-flux divergence field. At Haut Glacier d’Arolla, Valais, Switzerland, the mass-balance distribution between September 1992 and September 1993 is calculated at 20 m resolution from the difference between the pattern of surface-elevation change derived from analytical photogrammetry and the mass-flux divergence field determined from three-dimensional, numerical flow modelling constrained by surface-velocity measurements. The resultant pattern of mass balance is almost totally negative, showing a strong dependence on elevation, but with large localized departures. The computed distribution of mass balance compares well ( R 2 = 0.91) with mass-balance measurements made at stakes installed along the glacier centre line over the same period. Despite the highly optimized nature of the flow-modelling effort employed in this study, the good agreement indicates the potential this method has as a strategy for deriving high spatial and temporal-resolution estimates of mass balance.

BedMachine v3: Complete Bed Topography and Ocean Bathymetry Mapping of Greenland From Multibeam Echo Sounding Combined With Mass Conservation

Abstract Greenlands bed topography is a primary control on ice flow, grounding line migration, calving dynamics, and subglacial drainage. Moreover, fjord bathymetry regulates the penetration of warm Atlantic water (AW) that rapidly melts and undercuts Greenlands marine‐terminating glaciers. Here we present a new compilation of Greenland bed topography that assimilates seafloor bathymetry and ice thickness data through a mass conservation approach. A new 150 m horizontal resolution bed topography/bathymetric map of Greenland is constructed with seamless transitions at the ice/ocean interface, yielding major improvements over previous data sets, particularly in the marine‐terminating sectors of northwest and southeast Greenland. Our map reveals that the total sea level potential of the Greenland ice sheet is 7.42 ± 0.05 m, which is 7 cm greater than previous estimates. Furthermore, it explains recent calving front response of numerous outlet glaciers and reveals new pathways by which AW can access glaciers with marine‐based basins, thereby highlighting sectors of Greenland that are most vulnerable to future oceanic forcing.

Subglacial water drainage, storage, and piracy beneath the Greenland ice sheet

Meltwater drainage across the surface of the Greenland Ice Sheet (GrIS) is well constrained by measurements and modeling, yet despite its critical role, knowledge of its transit through the subglac ...

Evidence for subglacial ponding across Taylor Glacier, Dry Valleys, Antarctica

Abstract Ice-penetrating radar and modelling data are presented suggesting the presence of a zone of temperate ice, water ponding or saturated sediment beneath the tongue of Taylor Glacier, Dry Valleys, Antarctica. The proposed subglacial zone lies 3–6 km up-glacier of the terminus and is 400– 1000m across. The zone coincides with an extensive topographic overdeepening to 80m below sea level. High values of residual bed reflective power across this zone compared to other regions and the margins of the glacier require a high dielectric contrast between the ice and the bed and are strongly indicative of the presence of basal water or saturated sediment. Analysis of the hydraulic equipotential surface also indicates strong convergence into this zone of subglacial water flow paths. However, thermodynamic modelling reveals that basal temperatures in this region could not exceed –7˚C relative to the pressure-melting point. Such a result is at odds with the radar observations unless the subglacial water is a hypersaline brine.

Seismic evidence of mechanically weak sediments underlying Russell Glacier, West Greenland

Abstract Amplitude-versus-angle (AVA) analysis of a seismic reflection line, imaged 13 km from Russell Glacier terminus, near the western margin of the Greenland ice sheet (GrIS), suggests the presence of sediment at the bed. The analysis was complicated by the lack of identifiable multiples in the data due to a highly irregular and crevassed ice surface, rendering deeper seismic returns noisy. A modified technique for AVA processing of glacial seismic data using forward modelling with primary reflection amplitudes and simulated multiple amplitudes is presented here. Our analysis demonstrates that AVA analysis can be applied to areas with noisy seismic returns and indicates that sediment underlies the seismic study site. Our data are inconsistent with the common assumption that the GrIS is underlain only by hard bedrock, but consistent with the presence of subglacial sediment with porosity between 30% and 40%. As analysis and modelling of ice-sheet dynamics requires a sound knowledge of the underlying basal materials, subglacial sediment should be taken into account when considering ice dynamics in this region of the GrIS.