skip to main content


Title: Hydrologic Properties of a Highly Permeable Firn Aquifer in the Wilkins Ice Shelf, Antarctica
Abstract

We present measurements of the density, hydraulic conductivity, and specific discharge of a widespread firn aquifer in Antarctica, within the Wilkins Ice Shelf. At the field site, the aquifer is 16.2 m thick, starting at 13.4 m from the snow surface and transitioning from water‐saturated firn to ice at 29.6 m. Hydraulic conductivity derived from slug tests show a geometric mean value of 1.4 ± 1.2 × 10−4 m s−1, equivalent to permeability of 2.6 ± 2.2 × 10−11 m2. A borehole dilution test indicates an average specific discharge value of 1.9 ± 2.8 × 10−6 m s−1. Ground‐penetrating radar profiles and a groundwater flow model show the aquifer is draining laterally into a large nearby rift. Our findings indicate that the firn aquifer in the vicinity of the field site is likely not in a steady state and its presence likely contributed to past ice shelf instability.

 
more » « less
NSF-PAR ID:
10362181
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
47
Issue:
22
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The hydrology of alpine and subalpine areas in the Colorado Front Range (USA) is evolving, driven by warming and by the alteration of precipitation patterns, the timing of snowmelt, and other components of the hydrologic budget. Field measurements of soil hydraulic conductivity and moisture along 30-m transects (n = 13) of representative soils developed in surficial deposits and falling head slug tests of shallow groundwater in till demonstrate that hydraulic conductivity in the soil is comparable to hydraulic conductivity values in the shallow aquifer. Soil hydraulic conductivity values were variable (medians ranged from 5.6 × 10−7 to 4.96 × 10−5 m s−1) and increased in alpine areas underlain by periglacial deposits. Hydraulic conductivities measured by a modified Hvorslev technique in test wells ranged from 4.86 × 10−7 to 1.77 × 10−4 m s−1 in subalpine till. The results suggest a gradient from higher hydraulic conductivity in alpine zones, where short travel paths through periglacial deposits support ephemeral streams and wetlands, to lower hydraulic conductivity in the till-mantled subalpine zone. In drier downstream areas, streambed infiltration contributes substantially to near-channel groundwater. As summer temperatures and evapotranspiration (ET) increase and snowmelt occur earlier, alpine soils are likely to become more vulnerable to drought, and groundwater levels in the critical zone may lower, affecting the connectivity between late-melting snow, meltwater streams, and the areas they affect downstream.

     
    more » « less
  2. Abstract

    A set of collocated, in situ oceanographic and glaciological measurements from Petermann Gletscher Ice Shelf, Greenland, provides insights into the dynamics of under‐ice flow driving basal melting. At a site 16 km seaward of the grounding line within a longitudinal basal channel, two conductivity‐temperature (CT) sensors beneath the ice base and a phase‐sensitive radar on the ice surface were used to monitor the coupled ice shelf‐ocean system. A 6 month time series spanning 23 August 2015 to 12 February 2016 exhibited two distinct periods of ice‐ocean interactions. Between August and December, radar‐derived basal melt rates featured fortnightly peaks of15 m yr−1which preceded the arrival of cold and fresh pulses in the ocean that had high concentrations of subglacial runoff and glacial meltwater. Estimated current speeds reached 0.20 – 0.40 m s−1during these pulses, consistent with a strengthened meltwater plume from freshwater enrichment. Such signals did not occur between December and February, when ice‐ocean interactions instead varied at principal diurnal and semidiurnal tidal frequencies, and lower melt rates and current speeds prevailed. A combination of estimated current speeds and meltwater concentrations from the two CT sensors yields estimates of subglacial runoff and glacial meltwater volume fluxes that vary between 10 and 80 m3 s−1during the ocean pulses. Area‐average upstream ice shelf melt rates from these fluxes are up to 170 m yr−1, revealing that these strengthened plumes had already driven their most intense melting before arriving at the study site.

     
    more » « less
  3. Abstract

    Seasonal formation of Dense Shelf Water (DSW) in the Ross Sea is a direct precursor to Antarctic Bottom Water, which fills the deep ocean with atmospheric gases in what composes the southern limb of the solubility pump. Measurements of seawater noble gas concentrations during katabatic wind events in two Ross Sea polynyas reveal the physical processes that determine the boundary value properties for DSW. This decomposition reveals 5–6 g kg−1of glacial meltwater in DSW and sea‐ice production rates of up to 14 m yr−1within the Terra Nova Bay polynya. Despite winds upwards of 35 m s−1during the observations, air bubble injection had a minimal contribution to gas exchange, accounting for less than 0.01 μmols kg−1of argon in seawater. This suggests the slurry of frazil ice and seawater at the polynya surface inhibits air‐sea exchange. Most noteworthy is the revelation that sea‐ice formation and glacial melt contribute significantly to the ventilation of DSW, restoring 10% of the gas deficit for krypton, 24% for argon, and 131% for neon, while diffusive gas exchange contributes the remainder. These measurements reveal a cryogenic component to the solubility pump and demonstrate that while sea ice blocks air‐sea exchange, sea ice formation and glacial melt partially offset this effect via addition of gases. While polynyas are a small surface area, they represent an important ventilation site within the southern‐overturning cell, suggesting that ice processes both enhance and hinder the solubility pump.

     
    more » « less
  4. Abstract

    The present work is the first to undertake systematicin situobservations of the microstructural changes on samples taken at ∼10‐m intervals along the length of a 80‐m firn core, extracted at Summit, Greenland (72°35’ N, 38°25’ W) in June, 2017, under interrupted load and at a strain rate of ∼8 × 10−5s−1at −10°C, using a X‐ray micro‐computed tomograph. Several noteworthy features of the densification were found: the ice particle size increases, while the specific surface area, the total porosity, the pore size, and the structure model index (a measure of convexity/concavity of ice surface) decreases. The results were used to formulate semi‐empirical models (valid in the stress range of ∼0.05–2.15 MPa) that can be used to quantitatively assess the relative contributions of lattice diffusion (LD) and grain boundary diffusion (GBD) to the densification of polar firn. We found that 0.28 and 2.15 MPa are two critical stresses, which represent the start and end of LD as the dominant deformation mechanism to the densification of polar firn under the interrupted increasing loads. This bimodality when LD dominates implies that stress is not the only factor governing the densification of polar firn. On the other hand, GBD dominates the densification of polar firn both for stresses lower than 0.28 MPa and greater than 2.15 MPa. At stresses greater than 2.41 MPa, the firn specimens either fractured or other deformation mechanisms dominated, e.g., grain boundary sliding or power‐law dislocation creep.

     
    more » « less
  5. Abstract

    Basal melting of ice shelves is a major source of mass loss from the Antarctic Ice Sheet. In situ measurements of ice shelf basal melt rates are sparse, while the more extensive estimates from satellite altimetry require precise information about firn density and characteristics of near‐surface layers. We describe a novel method for estimating multidecadal basal melt rates using airborne ice penetrating radar data acquired during a 3‐year survey of the Ross Ice Shelf. These data revealed an ice column with distinct upper and lower units whose thicknesses change as ice flows from the grounding line toward the ice front. We interpret the lower unit as continental meteoric ice that has flowed across the grounding line and the upper unit as ice formed from snowfall onto the relatively flat ice shelf. We used the ice thickness difference and strain‐induced thickness change of the lower unit between the survey lines, combined with ice velocities, to derive basal melt rates averaged over one to six decades. Our results are similar to satellite laser altimetry estimates for the period 2003–2009, suggesting that the Ross Ice Shelf melt rates have been fairly stable for several decades. We identify five sites of elevated basal melt rates, in the range 0.5–2 m a−1, near the ice shelf front. These hot spots indicate pathways into the sub‐ice‐shelf ocean cavity for warm seawater, likely a combination of summer‐warmed Antarctic Surface Water and modified Circumpolar Deep Water, and are potential areas of ice shelf weakening if the ocean warms.

     
    more » « less