Significant groundwater contribution to Antarctic ice streams hydrologic budget

Satellite observations have revealed active hydrologic systems beneath Antarctic ice streams, but sources and sinks of water within these systems are uncertain. Here we use numerical simulations of ice streams to estimate the generation, flux, and budget of water beneath five ice streams on the Siple Coast. We estimate that 47% of the total hydrologic input (0.98 km3 yr−1) to Whillans (WIS), Mercer (MIS), and Kamb (KIS) ice streams comes from the ice sheet interior and that only 8% forms by local basal melting. The remaining 45% comes from a groundwater reservoir, an overlooked source in which depletion significantly exceeds recharge. Of the total input to Bindschadler (BIS) and MacAyeal (MacIS) ice streams (0.56 km3 yr−1), 72% comes from the interior, 19% from groundwater, and 9% from local melting. This contrasting hydrologic setting modulates the ice streams flow and has important implications for the search for life in subglacial lakes.This work was carried out with support from the Isaac Newton Trust, Cecil H. and Ida M. Green Foundation and Natural Environment Research Council (grant NE/E005950/1 and NE/J005800/1).This is the final version of the article. It was originally published in Geophysical Research Letters and is also available from the Wiley website at http://onlinelibrary.wiley.com/doi/10.1002/2014GL059250/abstract. © American Geophysical Union 201

Reactivation of Kamb Ice Stream tributaries triggers century-scale reorganization of Siple Coast ice flow in West Antarctica

Ongoing, centennial-scale flow variability within the Ross ice streams of West Antarctica suggests that the present-day positive mass balance in this region may reverse in the future. Here, we use a three-dimensional ice-sheet model to simulate ice flow in this region over 250 years. The flow responds to changing basal properties, as a subglacial till layer interacts with water transported in an active subglacial hydrological system. We show that a persistent weak bed beneath the tributaries of the dormant Kamb Ice Stream is a source of internal ice-flow instability, which reorganizes all ice streams in this region, leading to a reduced (positive) mass balance within decades and a net loss of ice within two centuries. This hitherto unaccounted for flow variability could raise sea-level by 5mm this century. Better constraints on future sea-level change from this region will require improved estimates of geothermal heat flux and subglacial water transport.This work was carried out with support from the Isaac Newton trust, Cecil H. and Ida M. Green Foundation and Natural Environment Research Council (grants NE/E005950/1 and NE/J005800/1). SFP was supported by the U.S. Department of Energy Office of Science, Biological and Environmental Research program. ST acknowledges support from National Science Foundation (grant #0338295). SPC was supported by funding from the Cryospheric Sciences program of NASA and HAF was supported by funding from NSF (grant ANT-0838885 (Fricker)). The source code for the results presented can be obtained by contacting the corresponding author directlyThis is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1002/2015GL06578

Crater Lakes on Mars: Development of Quantitative Thermal and Geomorphic Models

Impact craters on Mars have served as catchments for channel-eroding surface fluids, and hundreds of examples of candidate paleolakes are documented [1,2] (see Figure 1). Because these features show similarity to terrestrial shorelines, wave action has been hypothesized as the geomorphic agent responsible for the generation of these features [3]. Recent efforts have examined the potential for shoreline formation by wind-driven waves, in order to turn an important but controversial idea into a quantitative, falsifiable hypothesis. These studies have concluded that significant wave-action shorelines are unlikely to have formed commonly within craters on Mars, barring Earth-like weather for approx.1000 years [4,5,6]

A new methodology to simulate subglacial deformation of water-saturated granular material

The dynamics of glaciers are to a large degree governed by processes operating at the ice–bed interface, and one of the primary mechanisms of glacier flow over soft unconsolidated sediments is subglacial deformation. However, it has proven difficult to constrain the mechanical response of subglacial sediment to the shear stress of an overriding glacier. In this study, we present a new methodology designed to simulate subglacial deformation using a coupled numerical model for computational experiments on grain-fluid mixtures. The granular phase is simulated on a per-grain basis by the discrete element method. The pore water is modeled as a compressible Newtonian fluid without inertia. The numerical approach allows close monitoring of the internal behavior under a range of conditions. <br><br> Our computational experiments support the findings of previous studies where the rheology of a slowly deforming water-saturated granular bed in the steady state generally conforms to the rate-independent plastic rheology. Before this so-called critical state, deformation is in many cases accompanied by volumetric changes as grain rearrangement in active shear zones changes the local porosity. For previously consolidated beds porosity increases can cause local pore-pressure decline, dependent on till permeability and shear rate. We observe that the pore-water pressure reduction strengthens inter-granular contacts, which results in increased shear strength of the granular material. In contrast, weakening takes place when shear deformation causes consolidation of dilated sediments or during rapid fabric development. Both processes of strengthening and weakening depend inversely on the sediment permeability and are transient phenomena tied to the porosity changes during the early stages of shear. <br><br> We find that the transient strengthening and weakening in turn influences the distribution of shear strain in the granular bed. Dilatant strengthening has the ability to distribute strain during early deformation to large depths, if sediment dilatancy causes the water pressure at the ice–bed interface to decline. Oppositely, if the ice–bed interface is hydrologically stable the strengthening process is minimal and instead causes shallow deformation. The depth of deformation in subglacial beds thus seems to be governed by not only local grain and pore-water feedbacks but also larger-scale hydrological properties at the ice base

Inferring ice fabric from birefringence loss in airborne radargrams : application to the eastern shear margin of Thwaites Glacier, West Antarctica

This work is ITGC Contribution No. ITGC-036 and is an output from the Thwaites Interdisciplinary Margin Evolution (TIME) project as part of the International Thwaites Glacier Collaboration (ITGC), supported by Natural Environment Research Council (NERC) research grant #NE/S006788/1 supporting T. J. Young and P. Christoffersen, and National Science Foundation (NSF) research grant #1739027 supporting S. M. Tulaczyk and D. M. Schroeder. Logistics for this project were provided by the NSF-U.S. Antarctic Program and NERC-British Antarctic Survey. R. Culberg is supported by a USA Department of Defense NDSEG Fellowship and N. L. Bienert is supported by an NSF Graduate Research Fellowship.In airborne radargrams, undulating periodic patterns in amplitude that overprint traditional radiostratigraphic layering are occasionally observed, however, they have yet to be analyzed from a geophysical or glaciological perspective. We present evidence supported by theory that these depth-periodic patterns are consistent with a modulation of the received radar power due to the birefringence of polar ice, and therefore indicate the presence of bulk fabric anisotropy. Here, we investigate the periodic component of birefringence-induced radar power recorded in airborne radar data at the eastern shear margin of Thwaites Glacier and quantify the lateral variation in azimuthal fabric strength across this margin. We find the depth variability of birefringence periodicity crossing the shear margin to be a visual expression of its shear state and its development, which appears consistent with present-day ice deformation. The morphology of the birefringent patterns is centered at the location of maximum shear and observed in all cross-margin profiles, consistent with predictions of ice fabric when subjected to simple shear. The englacial fabric appears stronger inside the ice stream than outward of the shear margin. The detection of birefringent periodicity from non-polarimetric radargrams presents a novel use of subsurface radar to constrain lateral variations in fabric strength, locate present and past shear margins, and characterize the deformation history of polar ice sheets.Publisher PDFPeer reviewe

Rapid and accurate polarimetric radar measurements of ice crystal fabric orientation at the Western Antarctic Ice Sheet (WAIS) Divide ice core site

The crystal orientation fabric (COF) of ice sheets records the past history of ice sheet deformation and influences present-day ice flow dynamics. Though not widely implemented, coherent ice-penetrating radar is able to detect bulk anisotropic fabric patterns by exploiting the birefringence of ice crystals at radar frequencies, with the assumption that one of the crystallographic axes is aligned in the vertical direction. In this study, we conduct a suite of quad-polarimetric measurements consisting of four orthogonal antenna orientation combinations near the Western Antarctic Ice Sheet (WAIS) Divide Ice Core site. From these measurements, we are able to quantify the azimuthal fabric asymmetry at this site to a depth of 1400 m at a bulk-averaged resolution of up to 15 m. Our estimates of fabric asymmetry closely match corresponding fabric estimates directly measured from the WAIS Divide Ice Core. While ice core studies are often unable to determine the absolute fabric orientation due to core rotation during extraction, we are able to identify and conclude that the fabric orientation is depth-invariant to at least 1400 m, equivalent to 6700 years BP (years before 1950), and aligns closely with the modern surface strain direction at WAIS Divide. Our results support the claim that the deformation regime at WAIS Divide has not changed substantially through the majority of the Holocene. Rapid polarimetric determination of bulk fabric asymmetry and orientation compares well with much more laborious sample-based COF measurements from thin ice sections. Because it is the bulk-averaged fabric that ultimately influences ice flow, polarimetric radar methods provide an opportunity for its accurate and widespread mapping and its incorporation into ice flow models.Natural Environment Research Council (NERC) research grant NE/S006788/1 National Science Foundation (NSF) research grant #173902

Petrology and Geochemistry of the AND-1B Core, ANDRILL McMurdo Ice Shelf Project, Antarctica

This section reports preliminary data and results on petrology and geochemistry of AND-1B core

Roughness of a subglacial conduit under Hansbreen, Svalbard

K.M., J.G., X.L. and Y.C. were supported by the National Science Foundation (NSF) under Grant No. #1503928. Thefieldwork team (K.M., J.G., M.C.) were supported by the Norwegian Arctic Research Council and Svalbard Science Forum, RiS #6106. K.M. was also supported by the National Aeronautics and Space Administration (NASA)Headquarters under the NASA Earth and Space Science Fellowship Program – Grant NNX10AN83H, the University of California, Santa Cruz, and the Woods Hole Oceanographic Institution Ocean and Climate Change Institute post-graduate fellowship. Portions of this work were conducted while J.G. was supported by the NSF EAR Postdoctoral Fellowship (#0946767). S.T. was funded by NASA grant NNX11AH61G.Hydraulic roughness exerts an important but poorly understood control on water pressure in subglacial conduits. Where relative roughness values are 5%. Here we report the first quantitative assessment of roughness heights and hydraulic diameters in a subglacial conduit. We measured roughness heights in a 125 m long section of a subglacial conduit using structure-from-motion to produce a digital surface model, and hand-measurements of the b-axis of rocks. We found roughness heights from 0.07 to 0.22 m and cross-sectional areas of 1-2 m2, resulting in relative roughness of 3-12% and >5% for most locations. A simple geometric model of varying conduit diameter shows that when the conduit is small relative roughness is >30% and has large variability. Our results suggest that parameterizations of conduit hydraulic roughness in subglacial hydrological models will remain challenging until hydraulic diameters exceed roughness heights by a factor of 20, or the conduit radius is >1 m for the roughness elements observed here.Publisher PDFPeer reviewe