Radar-Detected Englacial Stratigraphy in the Pensacola Mountains, Antarctica: Implications for Recent Changes in Ice Flow and Accumulation
Total Page:16
File Type:pdf, Size:1020Kb
Annals of Glaciology 54(63) 2013 doi: 10.3189/2013AoG63A371 91 Radar-detected englacial stratigraphy in the Pensacola Mountains, Antarctica: implications for recent changes in ice flow and accumulation Seth CAMPBELL,1,2 Greg BALCO,3 Claire TODD,4 Howard CONWAY,5 Kathleen HUYBERS,5 Christopher SIMMONS,6 Michael VERMEULEN4 1University of Maine, Orono, ME, USA E-mail: [email protected] 2US Army Cold Regions Research and Engineering Laboratory (CRREL), Hanover, NH, USA 3Berkeley Geochronology Center, Berkeley, CA, USA 4Pacific Lutheran University, Tacoma, WA, USA 5University of Washington, Seattle, WA, USA 6International Federation of Mountain Guide Associations (IFMGA)/American Mountain Guides Association (AMGA) Mountain Guide, Seattle, WA, USA ABSTRACT. We used measurements of radar-detected stratigraphy, surface ice-flow velocities and accumulation rates to investigate relationships between local valley-glacier and regional ice-sheet dynamics in and around the Schmidt Hills, Pensacola Mountains, Antarctica. Ground-penetrating radar profiles were collected perpendicular to the long axis of the Schmidt Hills and the margin of Foundation Ice Stream (FIS). Within the valley confines, the glacier consists of blue ice, and profiles show internal stratigraphy dipping steeply toward the nunataks and truncated at the present-day ablation surface. Below the valley confines, the blue ice is overlain by firn. Data show that upward-progressing overlap of actively accumulating firn onto valley-glacier ice is slightly less than ice flow out of the valleys over the past 1200 years. The apparent slightly negative mass balance (–0.25 cm a–1) suggests that ice-margin elevations in the Schmidt Hills may have lowered over this time period, even without a change in the surface elevation of FIS. Results suggest that (1) mass-balance gradients between local valley glaciers and regional ice sheets should be considered when using local information to estimate regional ice surface elevation changes; and (2) interpretation of shallow ice structures imaged with radar can provide information about local ice elevation changes and stability. INTRODUCTION combined with surface ice-flow velocities and accumulation Reconstructions of past changes in the thickness and extent rates, estimate ages of buried unconformities (e.g. Campbell of the Antarctic ice sheet are important for understanding and others, 2012) and infer differences in accumulation rates past and present sea-level change, and for validating numer- (e.g. Arcone and others, 2012). ical models, necessary to make realistic predictions of In recent years, cosmogenic-nuclide exposure dating changes under future possible environmental conditions. methods (Gosse and Phillips, 2001; Dunai, 2010; Balco, Direct observations of ice-sheet history from satellite data 2011) have been used to date glacial deposits in Antarctica extend back until the 1970s; obtaining longer records (e.g. Stone and others, 2003; Todd and others, 2010). Meas- requires the use of geologic and glaciological data. For urements of trace nuclides produced by cosmic-ray bom- example, radar-detected layers and dated ice cores have bardment of rock surfaces are used to determine the exposure been used to infer changes in ice thickness through the age of the surface. Although this technique has been used to Holocene (e.g. Conway and others, 1999; Waddington and reconstruct the overall pattern of LGM-to-present ice-sheet others, 2005; Price and others, 2007). Marine geologic data thinning in Antarctica, it has several limitations: (1) measure- have been used to infer ice-margin positions and the ment precision decreases at low nuclide concentrations, so presence or absence of ice shelves during the Last Glacial exposure ages from sites that have been exposed for less than Maximum (LGM) (e.g. Domack and others, 2005; Hillen- 1000 years can be subject to large uncertainty; (2) when ice brand and others, 2012; Stolldorf and others, 2012). surface-elevation changes are not monotonic, i.e. ice-sheet Terrestrial geologic data, including drift sheets and moraines thinning is interrupted or followed by thickening, exposure- deposited by outlet glaciers, have also provided important age records are ambiguous; (3) exposure ages of samples constraints on past ice thickness. Mapping of these deposits is collected from outcrops may reflect changes in local rather especially valuable when they can be dated by radiocarbon than regional flow, because the pattern of flow around or cosmogenic-nuclide exposure methods (e.g. Bockheim outcrops is typically complex. and others, 1989; Hall, 2009; Bentley and others, 2010; Here we evaluate the hypothesis that the internal stratig- Bromley and others, 2012; Whitehouse and others, 2012). raphy of ice adjacent to nunataks contains information about Radar has also been used to locate past ice-stream margin changes in ice surface elevation and flow dynamics over the positions (Catania and others, 2006), define boundaries past hundreds to thousands of years. Subsurface and between flow elements (e.g. King, 2009) and, when exposed nunataks affect ice flow and surface mass balance 92 Campbell and others: Englacial stratigraphy in the Pensacola Mountains Fig. 1. Locator maps of the Schmidt Hills and Foundation Ice Stream showing: (a) elevations from the RADARSAT-1 Antarctic Mapping Project (RAMP) digital elevation model (DEM), the current grounding-line interpretation from the Antarctic Digital Database (ADD) and the location of the Schmidt Hills in the Pensacola Mountains relative to FIS and Ronne–Filchner Ice Shelf; (b) QuickBird 0.6 m pixel resolution satellite image of the Schmidt Hills showing location of the Childs Glacier/FIS flow relative to GPR profiles 1–4 (white solid lines, L1–L4), GPR profiles collected but not imaged (white dashed lines), surface ice-flow velocities (black arrows, underlined lettering; m a–1), the meltwater pond near line 2 (MW), blue-ice areas (BIA), accumulation rates (circles; cm a–1) and the location of the photo shown as (c) (white arrow, P); and (c) photo with approximate radar profile locations (dotted lines; L1–L4), the shear margin and FIS with dotted arrow pointing towards the Ronne–Filchner Ice Shelf. Imagery #2012 DigitalGlobe, Inc. as the ice thickness changes. For example, ice will flow FOUNDATION ICE STREAM AND THE SCHMIDT around rather than over a topographic obstacle as it emerges HILLS through the ice sheet, and protruding nunataks influence patterns of snow transport by wind, resulting in localized Foundation Ice Stream (FIS) is the southernmost ice stream areas of accumulation and ablation. Such changes would be that flows into the Ronne–Filchner Ice Shelf in the Weddell 2 evident in the radar-detected stratigraphy. Following this Sea Embayment. The catchment of FIS ( 566 000 km ) reasoning, information from patterns of shallow stratigraphy consists of ice originating from both East and West near nunataks should help interpret exposure-age records of Antarctica (Joughin and Bamber, 2005; Rignot and others, ice surface-elevation change over the past millennia. Also, 2011). Thus, thickness and velocity changes in the lower FIS interpretation of local shallow stratigraphy may provide will contain a record of changes in both regions. The insight into similarities or differences in local vs regional Schmidt Hills are a group of nunataks situated near the accumulation and flow patterns, which is useful for regional eastern edge of FIS as it enters the Filchner Trough extrapolation of exposure-age data. embayment (Fig. 1) 40 km downstream of the currently We pursue this idea by conducting an ice-penetrating estimated grounding line (SCAR, 2012). The hills are radar study of the shallow (0–400 m) firn and ice surround- approximately 27 km long and 12 km wide, with nunatak ing the Schmidt Hills, Pensacola Mountains, Antarctica. We elevations reaching >700 m above the current ice surface. also measured surface ice-flow velocity and accumulation Childs Glacier drains from the southern boundary of the rate to help interpret the radar data. We show that the Schmidt Hills and turns north to merge into the eastern edge shallow-ice stratigraphy around these nunataks: (1) is of FIS. Other small valley glaciers originating within the complex; (2) reflects the interplay of through-flowing ice Schmidt Hills flow westward into Childs Glacier and FIS. associated with the adjacent Childs Glacier/Foundation Ice Graf and others (1999) collected firn cores along a Stream and locally sourced ice derived from alpine glaciers; traverse up FIS; these cores show a regional pattern of and (3) displays features that may reflect change in the increasing accumulation rate from the FIS grounding line to surface elevation of the ice stream. the north toward the ice-shelf edge. They estimated a Campbell and others: Englacial stratigraphy in the Pensacola Mountains 93 Fig. 2. Line 1 (L1) 100 MHz GPR profile. (a) Elevation-corrected full depth profile with labels for bedrock (BR), blue-ice areas (BIA), an ablation surface (AS), surface tensional ice fractures (Fx) and the approximate division between confined (valley) and unconfined (Open) ice flow. Depth is based on a dielectric permittivity (DI) of 3.0 and two-way travel time (TWTT). (b) Zoom from (a) showing detail of dipping stratigraphy in the BIA, AS and hyperbolas caused by the Fx. (c) Zoom from (a) showing detail of the firn overlap described as wind- deposited bedding sequences (BS), buried-firn–blue-ice unconformity (UC) and the ice beneath the UC which lacks internal stratigraphy. 10.1 Æ 3.0 cm w.e. a–1 average accumulation rate between lying the scoured and exposed blue ice within the Schmidt 1947 and 1994 at their core closest to the Schmidt Hills, Hills valleys, where winds are more confined (e.g. Bintanja, 20 km from the end of our ground-penetrating radar (GPR) 1999). In contrast, accumulation zones are situated near the profile 1 (L1 in Fig. 1b). Joughin and Bamber (2005) used mouth of each valley, where winds are not as confined. regionally derived elevation and velocity data to estimate regional mass balance through a flux-gate method.