Search

The surface mass balance (SMB) of West Antarctica is an important glaciological input to understanding polar climate and sea-level rise but with historically poor in situ data coverage. Previous studies demonstrate the utility of frequency-modulated continuous-wave radar to image subsurface layering in ice sheets, providing an additional source of data with which to estimate SMB. Traditional methods, however, require time-intensive manual oversight. Here, we present a probabilistic, fully automated approach to estimate annual SMB and uncertainties from radar echograms using successive peak-finding and weighted neighborhood search algorithms with the Monte Carlo simulations based on annual-layer likelihood scores. We apply this method to ground-based and airborne radar in a 175-km transect of the West Antarctic interior and compare the results to traditional manual methods and independent estimates from firn cores. The method demonstrates an automated estimation of SMB across a range of accumulation rates (100–450-mm water equivalent per year) and layer gradients up to 2 m/km. Based on likelihood-weighted F-scores, automated layer picks have a success rate between 64% and 84.6% compared with manually picked layers for three validation sites dispersed across the region. Comparisons between the automated SMB estimates and independent firn cores show a bias of 24 ± 70-mm water equivalent per year (12% ± 35% water equivalent of the in situ mean accumulation rate) although individual core site biases differ. This new approach permits the fully automated extraction of annual SMB rates and should be broadly and readily applicable to previously collected and ongoing radar data sets across polar regions.

As part of the NASA SnowEx 2020 campaign, we performed a time series experiment with NASA's UAVSAR, an airborne L-band InSAR, over 13 sites across 5 states. Six flights were performed (December-March), capturing a wide range of snow conditions. On Grand Mesa, Colorado, two of the InSAR overflights were well aligned with two airborne lidar surveys. We show that in a 4 km2 wind exposed area on the west end of Grand Mesa, the measured change in phase can be used to estimate change in snow depth. In this region we observed both scouring and drifting, with a dynamic range of ± 20 cm. We use the measured surface density, change in phase, and local incidence angle to estimate change in snow depth over approximately a 2 week period, with a RMSD of 4.7cm depth, and 0.9cm SWE. This technique shows promise in open regions in dry snow conditions, and may be a useful component of a global snow satellite concept.

Synthetic Aperture Radar (SAR) is the ideal instrument for observing polar regions, with its ability to see through clouds and operate during the long polar nights. Since the launch of ERS-1 in 1991, the international constellation of SARs has enabled great progress in our understanding of sea ice, ice sheets, and glaciers. With its launch in early 2023, NISAR will contribute to the fourth decade of extensive SAR observations of the cryosphere. In contrast to all other missions, it will operate in a permanently south-looking mode, providing complete routine coverage of Antarctica. Here we summarize the NISAR requirements and how the mission will greatly advance our observational capability for the Earth's remote ice-covered regions.

Perennial water storage in firn aquifers has been observed within the lower percolation zone of the southeast Greenland ice sheet. Spatially distributed seismic and radar observations, made ~50 km upstream of the Helheim Glacier terminus, reveal spatial variations of seismic velocity within a firn aquifer. From 1.65 to 1.8 km elevation, shear‐wave velocity (Vs) is 1,290 ± 180 m/s in the unsaturated firn, decreasing below the water table (~15 m depth) to 1,130 ± 250 m/s. Below 1.65 km elevation, Vs in the saturated firn is 1,270 ± 220 m/s. The compressional‐to‐shear velocity ratio decreases in the downstream saturated zone, from 2.30 ± 0.54 to 2.01 ± 0.46, closer to its value for pure ice (2.00). Consistent with colocated firn cores, these results imply an increasing concentration of ice in the downstream sites, reducing the porosity and storage potential of the firn likely caused by episodic melt and freeze during the evolution of the aquifer. Plain Language Summary An integrated geophysical analysis of seismic, radar, and borehole measurements has been completed over a firn aquifer in southeast Greenland. We show the stiffness of the aquifer increases at lower elevations, closer to sea level, which leads to a decrease in pore space for the meltwater to be stored. This corresponds to an increase in ice content within the firn at lower elevations, as observed in borehole measurements, and likely caused by the meltwater refreezing within and below the aquifer.

Warming of the polar ice sheets causes changes in the hydrological regime of surface layers of firn and ice. Surface meltwater may undergo perennial storage of liquid water above the firn-ice transition, which could slow sea level rise or cause sudden release events, when storage capacity is reached. Firn aquifers have been commonly observed within the lower percolation zone of the southeastern Greenland ice sheet during the past decade, and more recently, across some Antarctic ice shelves. Knowledge of the geographic extent and fractional liquid water content (and storage) of such aquifers will enable a better understanding of their effects on the sub- and en-glacial hydrologic system and is crucial for accurate predictions of the contribution of meltwater discharge to global sea level rise.Quantitative geophysical analysis from surface observations can be used to infer hydrological properties of the firn and ice without time intensive direct drilling, providing an efficient spatial distribution of properties along with an estimate of their uncertainty. Furthermore, by combining multiple types of geophysical observations, joint inversions allow ambiguities of one methodology to be mitigated by resolution in the other.Here, we demonstrate that this joint approach is a powerful complement to the conventional geophysical analysis of firn aquifers, by combining seismic, ground penetrating radar and borehole data to characterise aquifer properties, using the ‘MuLTI’ algorithm. In particular, we incorporate seismic shear wave velocities (Vs), derived from surface (Rayleigh) waves offering a promising means of distinguishing zones containing liquid water, into independent compressional wave velocity, density, and radar soundings of the water table. We find Vs decreases from 1600 m/s in the unsaturated firn above the water table at around 15 m depth, to 800 m/s through saturated ‘clean’ firn aquifer at around 25 m depth. However, at lower elevations, Vs increases to 1250 m/s through thicker, older firn aquifer where there are many ice lenses, which are interpreted to correspond with episodes of refreezing aquifer water as the system has evolved through time. With access to multiple seismic wave velocities (compressional and shear) through the aquifer, a more accurate estimate of liquid water content can be derived. Thus, the application of the MuLTI algorithm to this pressing new problem can deliver an accurate assessment of firn aquifer properties, and provide clear uncertainty limits which will be valuable for predictive modelling.

The Indus River is a vital resource for food security, ecosystem services, hydropower, and economy for millions of people living in Pakistan, India, China, and Afghanistan. Glacier and snowmelt from the high altitude Himalaya, Karakoram, and Hindu Kush mountain ranges are the largest drivers of discharge in the upper Indus Basin (UIB), and contribute significantly to Indus flows. Complex climatology and topography, coupled with the challenges of field study and meteorological measurement in these rugged ranges, elicit notable uncertainties in predicting seasonal runoff as well as cryospheric response to changes in climate. Here we utilize Sentinel-1 synthetic aperture radar (SAR) imagery to track ablation season development of wet snow in the Shigar Watershed of the Karakoram Mountains in Pakistan from 2015 to 2018. We exploit opportune local image acquisition times to highlight diurnal differences in radar indications of wet snow, and examine the spatial and temporal contexts of radar diurnal differences for 2015, 2017, and 2018 ablation seasons. Radar classifications for each ablation season show spatial and temporal patterns that indicate a dry winter snowpack undergoing diurnal surface melt-refreeze cycles, transitioning to surface snow that remains wet both day and night, and finally snow free conditions following melt out. Diurnally differing SAR signals may offer insights into important snowpack energy balance processes that precede melt out, which could provide useful constraints for both glacier mass balance modeling and runoff forecasting in remote alpine watersheds.

We apply the magnetic resonance sounding (MRS) method to investigate a firn aquifer in the south-east region of the Greenland ice sheet. Our study aims to delineate and estimate the volume of the recently discovered water stored within the firn (compacted snow) that remains liquid throughout the year. We develop and test successfully a methodology for joint use of MRS and ground-penetrating radar (GPR). This noninvasive geophysical approach is particularly well-adapted to glacier conditions and has a promising future for in situ investigation of water distribution in glaciers. At our field site, MRS showed an aquifer located at variable depths between 20 and 30 m beneath the ice-sheet surface. At the monitoring site, both MRS and GPR show an increase in the water volume stored between April 2015 and July 2016. MRS estimates suggest that the volume increased by approximately 28%.

Within the lower percolation zone of the southeastern Greenland ice sheet, meltwater has accumulated within the firn pore space, forming extensive firn aquifers. Previously, it was unclear if these aquifers stored or facilitated meltwater runoff. Following mixing of a saline solution into boreholes within the aquifer, we observe that specific conductance measurements decreased over time as flowing freshwater diluted the saline mixture in the borehole. These tests indicate that water flows through the aquifer with an average specific discharge of 4.3 × 10−6 m/s (σ = 2.5 × 10−6 m/s). The specific discharge decreases dramatically to 0 m/s, defining the bottom of the aquifer between 30 to 50 m depth. The observed flow indicates that the firn pore space is a short-term (

We document the existence of widespread firn aquifers in an elevation range of ~1200–2000 m, in the high snow-accumulation regions of the Greenland ice sheet. We use NASA Operation IceBridge accumulation radar data from five campaigns (2010–2014) to estimate a firn-aquifer total extent of 21,900 km2. We investigate two locations in Southeast Greenland, where repeated radar profiles allow mapping of aquifer-extent and water table variations. In the upper part of Helheim Glacier the water table rises in spring following above-average summer melt, showing the direct firn-aquifer response to surface meltwater production changes. After spring 2012, a drainage of the firn-aquifer lower margin (5 km) is inferred from both 750 MHz accumulation radar and 195 MHz multicoherent radar depth sounder data. For 2011–2014, we use a ground-penetrating radar profile located at our Ridgeline field site and find a spatially stable aquifer with a water table fluctuating less than 2.5 m vertically. When combining radar data with surface topography, we find that the upper elevation edge of firn aquifers is located directly downstream of locally high surface slopes. Using a steady state 2-D groundwater flow model, water is simulated to flow laterally in an unconfined aquifer, topographically driven by ice sheet surface undulations until the water encounters crevasses. Simulations suggest that local flow cells form within the Helheim aquifer, allowing water to discharge in the firn at the steep-to-flat transitions of surface topography. Supported by visible imagery, we infer that water drains into crevasses, but its volume and rate remain unconstrained.

We implemented an automated approach for mapping persistent ice and snow cover (PISC) across the conterminous western U.S. using all available Landsat TM and ETM+ scenes acquired during the late summer/early fall period between 2010 and 2014. Two separate validation approaches indicate this dataset provides a more accurate representation of glacial ice and perennial snow cover for the region than either the U.S. glacier database derived from US Geological Survey (USGS) Digital Raster Graphics (DRG) maps (based on aerial photography primarily from the 1960s–1980s) or the National Land Cover Database 2011 perennial ice and snow cover class. Our 2010–2014 Landsat-derived dataset indicates 28% less glacier and perennial snow cover than the USGS DRG dataset. There are larger differences between the datasets in some regions, such as the Rocky Mountains of Northwest Wyoming and Southwest Montana, where the Landsat dataset indicates 54% less PISC area. Analysis of Landsat scenes from 1987–1988 and 2008–2010 for three regions using a more conventional, semi-automated approach indicates substantial decreases in glaciers and perennial snow cover that correlate with differences between PISC mapped by the USGS DRG dataset and the automated Landsat-derived dataset. This suggests that most of the differences in PISC between the USGS DRG and the Landsat-derived dataset can be attributed to decreases in PISC, as opposed to differences between mapping techniques. While the dataset produced by the automated Landsat mapping approach is not designed to serve as a conventional glacier inventory that provides glacier outlines and attribute information, it allows for an updated estimate of PISC for the conterminous U.S. as well as for smaller regions. Additionally, the new dataset highlights areas where decreases in PISC have been most significant over the past 25–50 years.

Recent observations of the Greenland ice sheet show an increase of the area affected by progressive melt of snow and ice, thus resulting in production of the additional meltwater . In 2011, an important storage of meltwater in the firn has been observed in the S-E Greenland. This water does not freeze during the wintertime and forms a perennial firn aquifer. The aquifer spatial extent has been initially monitored with combined ground and airborne radar observations, but these geophysical techniques are not able to inform us on the amount of meltwater stored at depth. In this study, we use the magnetic resonance soundings (MRS) method for estimating the volume of water stored in the Greenland ice sheet firn and mapping its spatial variability. Our study area covers a firn aquifer along a 16-km E-W transect, ranging between elevations of 1520 and 1760 m. In July 2015 and July 2016, we performed MRS measurements that allow estimating the water volume in the studied area as well as the one-year water volume evolution. Water storage is not homogeneous, fluctuating between 0.2 and 2 m 3/m2 , and contains discontinuities in the hydrodynamic properties. We estimate an average volume of water stored in the firn in 2016 to be 0.76 m 3/m2, which corresponds to a 0.76-m-thick layer of bulk water. MRS monitoring reveals that from April 2015 to July 2016 the volume of water stored at the location of our transect increases by about 36%. We found MRS-estimated depth to water in a good agreement with that obtained with the ground penetrating radar (GPR).

The mean, trend and variability of net snow accumulation in firn cores are often used to validate model output, develop remote-sensing algorithms and quantify ice-sheet surface mass balance. Thus, accurately defining uncertainties associated with these in situ measurements is critical. In this study, we apply statistical simulation methods to quantify the uncertainty in firn-core accumulation data due to the uncertainty in depth–age scales. The methods are applied to a suite of firn cores from central West Antarctica. The results show that uncertainty in depth–age scales can give rise to spurious trends in accumulation that are the same order of magnitude as accumulation trends reported in West Antarctica. The depth–age scale uncertainties also significantly increase the apparent interannual accumulation variability, so these uncertainties must first be accounted for before using firn-core data to assess such processes as small-spatial-scale variability. Better quantification of error in accumulation will improve our ability to meaningfully compare firn-core data across different regions of the ice sheet, and provide appropriate targets for calibration and/or validation of model output and remote-sensing data.

The objectives of this research are to (1) address remote sensing strategies and requirements for estimating snow depth and snow water equivalent (SWE) using existing L-Band interferometric data sets in coordination with field-based observations and modeling frameworks and, with this information, (2) inform the Next Generation Cold Land Processes Experiment (SnowEx) toward articulating the appropriate science and research questions for a single motivating science plan. As proposed, SnowEx is a multi-year airborne snow campaign with a primary goal of exploring multimodal sensor observations in coordination with field campaigns to inform the next generation snow remote sensing satellite platform. Based on limitations of satellite-based optical and LiDAR instruments operating in regions of the globe with consistent cloud-cover, the fact that many snow-dominated regions are at more northerly latitudes (limited solar illumination in the middle of winter), and these snow-dominated regions often experience periods of prolonged cloud cover (due to synoptic precipitation events), a microwave remote sensing platform may be the most viable path to space for a dedicated snow remote sensing mission. Specifically, L-Band radar interferometry has shown some unique promise with an archive of historical and contemporary satellite collections from JAXA's PALSAR-1 and PALSAR-2 instruments, respectively. Moreover, with the expected NISAR (NASA-ISRO Synthetic Aperture Radar) mission launch in 2020 and the unprecedented availability of dedicated global interferometric L-Band products every 12-days, as well as what is in essence a NISAR airborne simulator in JPL's UAVSAR platform, the L-Band interferometric approach to estimating snow depth and snow water equivalent (SWE) requires further investigation within the context of in-situ observations and modeling frameworks.

Some regions of the Greenland ice sheet, where snow accumulation and melt rates are high, currently retain substantial volumes of liquid water within the firn pore space throughout the year. These firn aquifers, found between ~10 and 30 m below the snow surface, may significantly affect sea level rise by storing or draining surface meltwater. The hydraulic gradient and the hydraulic conductivity control flow of meltwater through the firn. Here we describe the hydraulic conductivity of the firn aquifer estimated from slug tests and aquifer tests at six sites located upstream of Helheim Glacier in southeastern Greenland. We conducted slug tests using a novel instrument, a piezometer with a heated tip that melts itself into the ice sheet. Hydraulic conductivity ranges between 2.5 × 10−5 and 1.1 × 10−3 m/s. The geometric mean of hydraulic conductivity of the aquifer is 2.7 × 10−4 m/s with a geometric standard deviation of 1.4 from both depth specific slug tests (analyzed using the Hvorslev method) and aquifer tests during the recovery period. Hydraulic conductivity is relatively consistent between boreholes and only decreases slightly with depth. The hydraulic conductivity of the firn aquifer is crucial for determining flow rates and patterns within the aquifer, which inform hydrologic models of the aquifer, its relation to the broader glacial hydrologic system, and its effect on sea level rise.

In spring of 2011, a perennial storage of water was observed in the firn of the southeastern Greenland ice sheet, a region of both high snow accumulation and high melt. This aquifer is created through percolation of surface meltwater downward through the firn, saturating the pore space above the ice-firn transition. The aquifer may play a significant role in sea level rise through storage or draining freshwater into the ocean. We carried out a series of active source seismic experiments using continuously refracted P-waves and inverted the first P-arrivals using a transdimensional Bayesian approach where the depth, velocity, and number of layers are allowed to vary to identify the seismic velocities associated with the base of the aquifer. When our seismic approach is combined with a radar sounding of the water table situated at the top of the firn aquifer, we are able to quantify the volume of water present. In our study region, the base of the aquifer lies on average 27.7±2.9 m beneath the surface, with an average thickness of 11.5±5.5 m. Using a Wyllie average for porosity, we found the aquifer has an average water content of 16±8%, with considerable variation in water storage capacity along the studied east-west flow line, 40 km upstream of the Helheim glacier terminus. Between 2015 and 2016, we observed a 1-2 km uphill expansion of the aquifer system, with a site dry in summer 2015 exhibiting a water content of 530 kg m-2 in summer 2016. We estimate the volume of water stored in the aquifer across the entire region upstream of Helheim glacier to be 4.7±3.1 Gt, approximately 3% of the total water stored in firn aquifers across the Greenland ice sheet. Elucidating the volume of water stored within these recently discovered aquifers is vital for determining the hydrological structure and stability of the southeastern Greenland ice sheet.

The Saint Elias Mountains in southern Alaska are located at a structural syntaxis where the coastal thrust and fold belt of the Fairweather plate boundary intersects thrust faults and folds generated by collision of the Yakutat Terrane. The axial trace of this syntaxis extends southeastward out of the Saint Elias Mountains and beneath Malaspina Glacier where it is hidden from view and cannot be mapped using conventional methods. Here we examine the surface morphology and flow patterns of Malaspina Glacier to infer characteristics of the bedrock topography and organization of the syntaxis. Faults and folds beneath the eastern part of the glacier trend northwest and reflect dextral transpression near the terminus of the Fairweather fault system. Those beneath the western part of the glacier trend northeast and accommodate folding and thrust faulting during collision and accretion of the Yakutat Terrane. Mapping the location and geometry of the structural syntaxis provides important constraints on spatial variations in seismicity, fault kinematics, and crustal shortening beneath Malaspina Glacier, as well as the position of the collisional deformation front within the Yakutat Terrane. We also speculate that the geometrical complexity of intersecting faults within the syntaxis formed a barrier to rupture propagation during two regional Mw 8.1 earthquakes in September 1899.

Recent observations have shown that the firn layer on the Greenland Ice Sheet features subsurface bodies of liquid water at the end of the winter season. Using a model with basic firn hydrology, thermodynamics, and compaction in one dimension, we find that a combination of moderate to strong surface melt and a high annual accumulation rate is required to form such a perennial firn aquifer. The high accumulation rate ensures that there is pore space available to store water at a depth where it is protected from the winter cold. Low-accumulation sites cannot provide sufficiently deep pore space to store liquid water. However, for even higher accumulation rates, the total cold content of the winter accumulation becomes sufficient to refreeze the total amount of liquid water. As a consequence, wintertime or springtime observations of subsurface liquid water in these specific accumulation conditions cannot distinguish between a truly perennial firn aquifer and water layers that will ultimately refreeze completely.

A perennial storage of water in a firn aquifer was discovered in southeast Greenland in 2011. We present the first in situ measurements of the aquifer, including densities and temperatures. Water was present at depths between approx. 12 and 37m and amounted to 18.7 +/- 0.9 kg in the extracted core. The water filled the firn to capacity at approx. 35m. Measurements show the aquifer temperature remained at the melting point, representing a large heat reservoir within the firn. Using model results of liquid water extent and aquifer surface depth from radar measurements, we extend our in situ measurements to the Greenland ice sheet. The estimated water volume is 140 +/- 20 Gt, representing approx. 0.4mm of sea level rise (SLR). It is unknown if the aquifer temporary buffers SLR or contributes to SLR through drainage and/or ice dynamics.

The Wind River Range in Wyoming, USA, contains the largest concentration of glacial mass in the Rocky Mountains of the contiguous USA. Despite this distinction, only a few field or remotely sensed studies providing glacier volume changes have been published. The current study focuses on Continental Glacier located on the northern end of the range and uses two field datasets (high-accuracy GPS surface elevation points and ice-penetrating radar transects of the glacier bed) to create a three-dimensional model of glacier volume. Current surface elevations are compared with historical elevation data to calculate surface elevation change over time. An average thinning rate of 13.8 ± 7.8 m (0.30 ± 0.17 m a–1) between 1966 and 2012 was found. Surface elevation change rates varied across the glacier, ranging from +0.30 to –0.98 m a–1. Taking into account variable melt rates across the glacier, along with a glacial volume of 72.1 × 106 ± 10.8 × 106 m3, we estimate that Continental Glacier will be reduced in volume by 43% over the next 100 years and will disappear completely over the next 300–400 years, if current climatic conditions persist.

A pulse is a type of unstable glacier flow intermediate between normal flow and surging. Using Landsat MSS, TM and ETM+ imagery and feature-tracking software, a time series of mostly annual velocity maps from 1973 to 2012 was produced that reveals five pulses of Ruth Glacier, Alaska. Peaks in ice velocity were found in 1981, 1989, 1997, 2003 and 2010, approximately every 7 years. During these peak years the ice velocity increased 300%, from approximately 40 m a–1 to 160 m a–1. Based on the spatio-temporal behavior of Ruth Glacier during the pulse cycles, we suggest the pulses are due to enhanced basal motion via deformation of a subglacial till. The cyclical nature of the pulses is interpreted to be due to a thin till, with low permeability, that causes incomplete drainage of the till between the pulses, followed by eventual recharge and dilation of the till. These findings suggest care is needed when attempting to correlate changes in regional climate with decadal-scale changes in velocity, because in some instances basal conditions may have a greater influence on ice dynamics than climate.

The accelerating loss of mass from the Greenland ice sheet is a major contribution to current sea level rise. Increased melt water runoff is responsible for half of Greenlands mass loss increase. Surface melt has been increasing in extent and intensity, setting a record for surface area melt and runoff in 2012. The mechanisms and timescales involved in allowing surface melt water to reach the ocean where it can contribute to sea level rise are poorly understood. The potential capacity to store this water in liquid or frozen form in the firn (multi-year snow layer) is significant, and could delay its sea-level contribution. Here we describe direct observation of water within a perennial firn aquifer persisting throughout the winter in the southern ice sheet,where snow accumulation and melt rates are high. This represents a previously unknown storagemode for water within the ice sheet. Ice cores, groundairborne radar and a regional climatemodel are used to estimate aquifer area (70 plue or minus 10 x 10(exp 3) square kilometers ) and water table depth (5-50 m). The perennial firn aquifer represents a new glacier facies to be considered 29 in future ice sheet mass 30 and energy budget calculations.

[1] Predicting how climate change will affect glacier and ice sheet flow speeds remains a large hurdle toward accurate sea level rise forecasting. Increases in surface melt rates are known to accelerate glacier flow in summer, whereas in winter, flow speeds are believed to be relatively invariant. Here we show that wintertime flow speeds on nearly all major glaciers throughout Alaska are not only variable but are inversely related to melt from preceding summers. For each additional meter of summertime melt, we observe an 11% decrease in wintertime velocity on glaciers of all sizes, geometries, climates, and bed types. This dynamic occurs because interannual differences in summertime melt affect how much water is retained in the subglacial system during winter. The ubiquity of the dynamic indicates it occurs globally on glaciers and ice sheets not frozen to their beds and thus constitutes a new mechanism affecting sea level rise projections.

[1] Observations of snow accumulation rates from five new firn cores show a negative trend that is statistically significant over the past several decades across the central West Antarctic ice sheet (WAIS). A negative temporal trend in accumulation rates is unexpected in light of rising surface temperatures as well as model simulations predicting higher accumulation rates for the region. Both the magnitude of the mean accumulation rates and the range of interannual variability observed in the new records compare favorably to older records collected from a broad area of the WAIS, suggesting that the new data may serve as a regional proxy for recent temporal trends in West Antarctic accumulation rates. The observed negative trend in accumulation is likely the result of a shift in low-pressure systems over the Amundsen Sea region, dominated by changes in the austral fall season. Regional-scale climate models and reanalysis data do not capture the negative trend in accumulation rate observed in these firn cores. Nevertheless the models and reanalyses agree well in both accumulation-rate means and interannual variability, with no single model or dataset standing out as significantly more skilled at capturing the observed magnitude of and trend in accumulation rates in this region of the WAIS.

Seasonal snow melt in the Wind River Range, Wyoming, has been ending earlier over the last several decades leaving the region to rely more on supplemental melt water from mountain glaciers. This leads to the necessity of understanding recent glacial changes. This study uses elevation data from 1966, 2006 and 2011 to calculate surface elevation and volume changes that have occurred on Continental Glacier. Results indicate a mean volume change of −0.034 ± 0.02 km3 and surface elevation change of −0.36 ± 0.19 m y−1 between 1966 and 2006. Detailed spatial analysis shows that the glacier is divided into two sections which are thinning at different rates (lower section: −0.06±0.19 m y−1; upper section: −0.51 ± 0.19 m y−1). The upper section has experienced 97% of the thinning (or 742.5 × 103 m3 of melt water equivalent per year) and increased its rate since 2006 by 27.5%.

Bering Glacier, Alaska, USA, has a ∼20 year surge cycle, with its most recent surge reaching the terminus in 2011. To study this most recent activity a time series of ice velocity maps was produced by applying optical feature-tracking methods to Landsat-7 ETM+ imagery spanning 2001-11. The velocity maps show a yearly increase in ice surface velocity associated with the down-glacier movement of a surge front. In 2008/09 the maximum ice surface velocity was 1.5 ±0.017 km a-1 in the mid-ablation zone, which decreased to 1.2 ±0.015 km a-1 in 2009/10 in the lower ablation zone, and then increased to nearly 4.4 ± 0.03 km a-1 in summer 2011 when the surge front reached the glacier terminus. The surge front propagated down-glacier as a kinematic wave at an average rate of 4.4 ±2.0 km a-1 between September 2002 and April 2009, then accelerated to 13.9 ± 2.0 km a-1 as it entered the piedmont lobe between April 2009 and September 2010. The wave seems to have initiated near the confluence of Bering Glacier and Bagley Ice Valley as early as 2001, and the surge was triggered in 2008 further down-glacier in the mid-ablation zone after the wave passed an ice reservoir area.

Despite containing only 14% of the Greenland ice sheet by area, the southeastern sector has the highest accumulation rates, and hence receives ∼30% of the total snow accumulation. We present accumulation rates obtained during our 2010 Arctic Circle Traverse derived from three 50 m firn cores dated using geochemical analysis. We tracked continuous internal reflection horizons between the firn cores using a 400 MHz ground-penetrating radar (GPR). GPR data combined with depth-age scales from the firn cores provide accumulation rates along a 70 km transect. We followed an elevation gradient from ∼2350 to ∼1830m to understand how progressive surface melt may affect the ability to chemically date the firn cores and trace the internal layers with GPR. From the firn cores, we find a 52% (∼0.43 m w.e. a-1) increase in average snow accumulation and greater interannual variability at the lower site than the upper site. The GPR profiling reveals that accumulation rates are influenced by topographic undulations on the surface, with up to 23% variability over 7 km. These measurements confirm the presence of high accumulation rates in the southeast as predicted by the calibrated regional climate model Polar MM5.

Contemporary climate warming over the Arctic is accelerating mass loss from the Greenland Ice Sheet through increasing surface melt, emphasizing the need to closely monitor its surface mass balance in order to improve sea-level rise predictions. Snow accumulation is the largest component of the ice sheet's surface mass balance, but in situ observations thereof are inherently sparse and models are difficult to evaluate at large scales. Here, we quantify recent Greenland accumulation rates using ultra-wideband (2–6.5 GHz) airborne snow radar data collected as part of NASA's Operation IceBridge between 2009 and 2012. We use a semiautomated method to trace the observed radiostratigraphy and then derive annual net accumulation rates for 2009–2012. The uncertainty in these radar-derived accumulation rates is on average 14 %. A comparison of the radar-derived accumulation rates and contemporaneous ice cores shows that snow radar captures both the annual and long-term mean accumulation rate accurately. A comparison with outputs from a regional climate model (MAR) shows that this model matches radar-derived accumulation rates in the ice sheet interior but produces higher values over southeastern Greenland. Our results demonstrate that snow radar can efficiently and accurately map patterns of snow accumulation across an ice sheet and that it is valuable for evaluating the accuracy of surface mass balance models.

A surge cycle of the Bering Glacier system, Alaska, is examined using observations of surface velocity obtained using synthetic aperture radar (SAR) offset tracking, and elevation data obtained from the University of Alaska Fairbanks LiDAR altimetry program. After 13 yr of quiescence, the Bering Glacier system began to surge in May 2008 and had two stages of accelerated flow. During the first stage, flow accelerated progressively for at least 10 months and reached peak observed velocities of ~ 7 m d−1. The second stage likely began in 2010. By 2011 velocities exceeded 9 m d−1 or ~ 18 times quiescent velocities. Fast flow continued into July 2011. Surface morphology indicated slowing by fall 2011; however, it is not entirely clear if the surge is yet over. The quiescent phase was characterized by small-scale acceleration events that increased driving stresses up to 70%. When the surge initiated, synchronous acceleration occurred throughout much of the glacier length. Results suggest that downstream propagation of the surge is closely linked to the evolution of the driving stress during the surge, because driving stress appears to be tied to the amount of resistive stress provided by the bed. In contrast, upstream acceleration and upstream surge propagation is not dependent on driving stress evolution.

Pressing scientific questions concerning the Greenland ice sheet's climatic sensitivity, hydrology, and contributions to current and future sea level rise require hydrological datasets to resolve. While direct observations of ice sheet meltwater losses can be obtained in terrestrial rivers draining the ice sheet and from lake levels, few such datasets exist. We present a new hydrologic dataset from previously unmonitored sites in the vicinity of Kangerlussuaq, Southwest Greenland. This dataset contains measurements of river stage and discharge for three sites along the Akuliarusiarsuup Kuua (Watson) River's northern tributary, with 30 min temporal resolution between June 2008 and July 2011. Additional data of water temperature, air pressure, and lake stage are also provided. Flow velocity and depth measurements were collected at sites with incised bedrock or structurally reinforced channels to maximize data quality. However, like most proglacial rivers, high turbulence and bedload transport introduce considerable uncertainty to the derived discharge estimates. Eleven propagating error sources were quantified, and reveal that largest uncertainties are associated with flow depth observations. Mean discharge uncertainties (approximately the 68% confidence interval) are two to four times larger (±19% to ±43%) than previously published estimates for Greenland rivers. Despite these uncertainties, this dataset offers a rare collection of direct measurements of ice sheet runoff to the global ocean and is freely available for scientific use at http://dx.doi.org/10.1594/PANGAEA.762818.

Rising sea levels and increased surface melting of the Greenland ice sheet have heightened the need for direct observations of meltwater release from the ice edge to ocean. Buoyant sediment plumes that develop in fjords downstream of outlet glaciers are controlled by numerous factors, including meltwater runoff. Here, Moderate Resolution Imaging Spectroradiometer (MODIS) satellite imagery is used to average surface suspended sediment concentration (SSC) in fjords around ~80% of Greenland from 2000–2009. Spatial and temporal patterns in SSC are compared with positive-degree-days (PDD), a proxy for surface melting, from the Polar MM5 regional climate model. Over this decade significant geographic covariance occurred between ice sheet PDD and fjord SSC, with outlet type (land- vs. marine-terminating glaciers) also important. In general, high SSC is associated with high PDD and/or a high proportion of land-terminating glaciers. Unlike previous site-specific studies of the Watson River plume at Kangerlussuaq, temporal covariance is low, suggesting that plume dimensions best capture interannual runoff dynamics whereas SSC allows assessment of long-term conditions across much broader fjord environments around the ice sheet. Remote sensing of both plume charactersitics thus offers a viable approach for observing spatial and temporal patterns of meltwater release exiting the Greenland ice sheet to the global ocean.