Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “firn”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗
NASA NTRS2022

Mapping Firn Saturation over Greenland using NASA’s Soil Moisture Active Passive Satellite

Mapping the spatial extent of recently identified englacial hydrological features (i.e., ice slabs and perennial firn aquifers) formed by meters-thick water-saturated firn layers over the percolation facies of the Greenland Ice Sheet using L-band microwave radiometry has recently been demonstrated. However, these initial maps are binary, and do not provide a parameter to estimate the spatial variability in the thickness and volumetric fraction of meltwater stored within the firn pore space. Here, we exploit enhanced-resolution vertical-polarization L-band brightness temperature (T^(B)_(V)) imagery (2015–2019) generated using observations collected over Greenland by NASA’s Soil Moisture Active Passive (SMAP) satellite and a simple two-layer L-band brightness temperature model. We map water-saturated firn layers via a “firn saturation” parameter, and interpret our results together with ice slab and perennial firn aquifer spatial extents, estimates of snow accumulation simulated via the Regional Atmospheric Climate Model (RACMOp2.3), and airborne radar surveys collected via NASA’s Operation IceBridge (OIB) campaigns. We find that variable firn saturation parameter values are mapped in lower snow accumulation ice slab areas in western, northern, and northeastern Greenland, where firn is colder and water-saturated firn layers seasonally refreeze as solid-ice. Higher firn saturation parameter values are mapped in higher snow accumulation perennial firn aquifer areas in southeastern, southern, and northwestern Greenland, where firn is near the melting point, and meters-thick water-saturated firn layers exist. Our results have implications for identifying expansive englacial reservoirs that store significant volumes of meltwater in locations that are vulnerable to meltwater-induced hydrofracturing and accelerated outlet glacier flow year-round.

Greenland ice sheet↗
NASA NTRS2022

An Empirical Algorithm to Map Perennial Firn Aquifers and Ice Slabs Within the Greenland Ice Sheet Using Satellite L-band Microwave Radiometry

Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature () imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

Julie Z Miller↗
NASA NTRS2022

Simulations of Firn Processes Over the Greenland and Antarctic Ice Sheets: 1980-2021

Conversion of altimetry-derived ice-sheet volume change to mass requires an understanding of the evolution of the combined ice and air content within the firn column. In the absence of suitable techniques to observe the changes to the firn column across the entirety of an ice sheet, the firn column processes are typically modeled. Here, we present new simulations of firn processes over the Greenland and Antarctic ice sheets (GrIS and AIS) using the Community Firn Model and atmospheric reanalysis variables for more than four decades. A data set of more than 250 measured depth–density profiles from both ice sheets provides the basis of the calibration of the dry-snow densification scheme. The resulting scheme results in a reduction in the rate of densification, relative to a commonly used semi-empirical model, through a decreased dependence on the accumulation rate, a proxy for overburden stress. The 1980–2020 modeled firn column runoff, when combined with atmospheric variables from MERRA-2, generates realistic mean integrated surface mass balance values for the Greenland (+390 Gt yr(exp −1)) and Antarctic (+2612 Gt yr(exp −1)) ice sheets when compared to published model-ensemble means. We find that seasonal volume changes associated with firn air content are on average approximately 2.5 times larger than those associated with mass fluxes from surface processes for the AIS and 1.5 times larger for the GrIS; however, when averaged over multiple years, ice and air-volume fluctuations within the firn column are of comparable magnitudes. Between 1996 and 2019, the Greenland Ice Sheet lost nearly 5% of its firn air content, indicating a reduction in the total meltwater retention capability. Nearly all (94 %) of the meltwater produced over the Antarctic Ice Sheet is retained within the firn column through infiltration and refreezing.

Greenland↗
NASA NTRS2018

Enhanced Firn Densification in High‐Accumulation Shear Margins of the NE Greenland Ice Stream

Firn thickness across the NE Greenland Ice Stream is a function of accumulated strain, with thinner firn in the high‐strain margins of the ice stream. We present a novel technique for extracting firn density from previously collected seismic reflection profiles and apply this technique across both shear margins of NE Greenland Ice Stream. Firn is up to 30 m thinner in the vicinity of the ice stream shear margins. Snow accumulation rates across the ice stream were calculated from airborne ice‐penetrating radar data, calibrated with ground‐based firn density measurements from a shallow core. We find that accumulation is ~20% higher in the shear margins compared to the surroundings. The higher density firn adjacent to shear margins is due to high along‐flow stresses that accelerate firn densification and develops despite the higher accumulation rate favoring lower density. These firn density variations influence subglacial hydropotential by changing the ice surface slope and overburden pressure and may influence subglacial water flow. These results demonstrate the importance of high‐resolution firn surveys in studies of shear‐margin dynamics.

ice stream↗
NASA NTRS2023

A New Model of Dry Firn Densification Constrained By Continuous Strain Measurements Near South Pole

Converting measurements of ice-sheet surface elevation change to mass change requires measurements of accumulation and knowledge of the evolution of the density profile in the firn. Most firn-densification models are tuned using measured depth-density profiles, a method which is based on an assumption that the density profile in the firn is invariant through time. Here we present continuous measurements of firn-compaction rates in 12 boreholes near the South Pole over a two-year period. To our knowledge, these are the first continuous measurements of firn compaction on the Antarctic Plateau. We use the data to derive a new firn-densification algorithm framed as a constitutive relationship. We also compare our measurements to compaction rates predicted by several existing firn densification models. Results indicate that an activation energy of 60 kJ mol (−1), a value within the range used by current models, best predicts the seasonal cycle in compaction rates on the Antarctic Plateau. Our results suggest models can predict firn-compaction rates with at best 7% uncertainty and cumulative firn compaction on a two year timescale with at best 8% uncertainty.

C. Max Stevens↗
NASA NTRS2021

Physics-Based SNOWPACK Model Improves Representation of Near-Surface Antarctic Snow and Firn Density

Estimates of snow and firn density are required for satellite-altimetry-based retrievals of ice sheet mass balance that rely on volume-to-mass conversions. Therefore, biases and errors in presently used density models confound assessments of ice sheet mass balance and by extension ice sheet contribution to sea level rise. Despite this importance, most contemporary firn densification models rely on simplified semi-empirical methods, which are partially reflected by significant modeled density errors when compared to observations. In this study, we present a new drifting-snow compaction scheme that we have implemented into SNOWPACK, a physics-based land surface snow model. We show that our new scheme improves existing versions of SNOWPACK by increasing simulated near-surface (defined as the top 10 m) density to be more in line with observations (near-surface bias reduction from −44.9 to −5.4 kg m−3). Furthermore, we demonstrate high-quality simulation of near-surface Antarctic snow and firn density at 122 observed density profiles across the Antarctic ice sheet, as indicated by reduced model biases throughout most of the near-surface firn column when compared to two semi-empirical firn densification models (SNOWPACK mean bias=−9.7 kg m−3, IMAU-FDM mean bias=−32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the near surface where firn density is most variable due to accumulation and compaction variability driven by synoptic weather and seasonal climate variability. Additionally, the GSFC-FDM exhibits lower mean density bias from 7–10 m (SNOWPACK bias=−22.5 kg m−3, GSFC-FDM bias=10.6 kg m−3) and throughout the entire near surface at high-accumulation sites (SNOWPACK bias=−31.4 kg m−3, GSFC-FDM bias=−4.7 kg m−3). However, we found that the performance of SNOWPACK did not degrade when applied to sites that were not included in the calibration of semi-empirical models. This suggests that SNOWPACK may possibly better represent firn properties in locations without extensive observations and under future climate scenarios, when firn properties are expected to diverge from their present state.

Eric Keenan↗
NASA NTRS2023

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗
NASA NTRS2001

Modeled Seasonal Variations of Firn Density Induced by Steady State Surface Air Temperature Cycle

Seasonal variations of firn density in ice-sheet firn layers have been attributed to variations in deposition processes or other processes within the upper firn. A recent high-resolution (mm scale) density profile, measured along a 181 m core from Antarctica, showed small-scale density variations with a clear seasonal cycle that apparently was not-related to seasonal variations in deposition or known near-surface processes (Gerland and others 1999). A recent model of surface elevation changes (Zwally and Li, submitted) produced a seasonal variation in firn densification, and explained the seasonal surface elevation changes observed by satellite radar altimeters. In this study, we apply our 1-D time-dependent numerical model of firn densification that includes a temperature-dependent formulation of firn densification based on laboratory measurements of grain growth. The model is driven by a steady-state seasonal surface temperature and a constant accumulation rate appropriate for the measured Antarctic ice core. The modeled seasonal variations in firn density show that the layers of snow deposited during spring to mid-summer with the highest temperature history compress to the highest density, and the layers deposited during later summer to autumn with the lowest temperature history compress to the lowest density. The initial amplitude of the seasonal difference of about 0.13 reduces to about 0.09 in five years and asymptotically to 0.92 at depth, which is consistent with the core measurements.

Jun, Li↗
NASA NTRS2014

Observations of Pronounced Greenland Ice Sheet Firn Warming and Implications for Runoff Production

Field measurements of shallow borehole temperatures in firn across the northern Greenland ice sheet are collected during May 2013. Sites first measured in 19521955 are revisited, showing long-term trends in firn temperature. Results indicate a pattern of substantial firn warming (up to +5.7C) at midlevel elevations (1400-2500 m) and little temperature change at high elevations (2500 m). We find that latent heat transport into the firn due to meltwater percolation drives the observed warming. Modeling shows that heat is stored at depth for several years, and energy delivered from consecutive melt events accumulates in the firn. The observed warming is likely not yet in equilibrium with recent melt production rates but captures the progression of sites in the percolation facies toward net runoff production.

ice sheet↗
NASA NTRS2011

Modeling of Firn Compaction for Estimating Ice-Sheet Mass Change from Observed Ice-Sheet Elevation Change

Changes in ice-sheet surface elevation are caused by a combination of ice-dynamic imbalance, ablation, temporal variations in accumulation rate, firn compaction and underlying bedrock motion. Thus, deriving the rate of ice-sheet mass change from measured surface elevation change requires information on the rate of firn compaction and bedrock motion, which do not involve changes in mass, and requires an appropriate firn density to associate with elevation changes induced by recent accumulation rate variability. We use a 25 year record of surface temperature and a parameterization for accumulation change as a function of temperature to drive a firn compaction model. We apply this formulation to ICESat measurements of surface elevation change at three locations on the Greenland ice sheet in order to separate the accumulation-driven changes from the ice-dynamic/ablation-driven changes, and thus to derive the corresponding mass change. Our calculated densities for the accumulation-driven changes range from 410 to 610 kg/cu m, which along with 900 kg/cu m for the dynamic/ablation-driven changes gives average densities ranging from 680 to 790 kg/cu m. We show that using an average (or "effective") density to convert elevation change to mass change is not valid where the accumulation and the dynamic elevation changes are of opposite sign.

Li, Jun↗
NASA NTRS2023

Evaluating Greenland Surface-Mass-Balance and Firn-Densification Data Using ICESat-2 Altimetry

Surface-mass-balance (SMB) and firn-densification (FD) models are widely used in altimetry studies as a tool to separate atmospheric-driven from ice-dynamics-driven ice-sheet mass changes and to partition observed volume changes into ice-mass changes and firn-air-content changes. Until now, SMB models have been principally validated based on comparison with ice core and weather station data or comparison with widely separated flight radar-survey flight lines. Firn-densification models have been primarily validated based on their ability to match net densification over decades, as recorded in firn cores, and the short-term time-dependent component of densification has rarely been evaluated at all. The advent of systematic ice-sheet-wide repeated ice-surface-height measurements from ICESat-2 (the Ice Cloud, and land Elevation Satellite, 2) allows us to measure the net surface-height change of the Greenland ice sheet at quarterly resolution and compare the measured surface-height differences directly with those predicted by three FD–SMB models: MARv3.5.11 (Modèle Atmosphérique Régional version 3.5.11) and GSFCv1.1 and GSFCv1.2 (the Goddard Space Flight Center FD–SMB models version 1.1 and 1.2). By segregating the data by season and elevation, and based on the timing and magnitude of modelled processes in areas where we expect minimal ice-dynamics-driven height changes, we investigate the models' accuracy in predicting atmospherically driven height changes. We find that while all three models do well in predicting the large seasonal changes in the low-elevation parts of the ice sheet where melt rates are highest, two of the models (MARv3.5.11 and GSFCv1.1) systematically overpredict, by around a factor of 2, the magnitude of height changes in the high-elevation parts of the ice sheet, particularly those associated with melt events. This overprediction seems to be associated with the melt sensitivity of the models in the high-elevation part of the ice sheet. The third model, GSFCv1.2, which has an updated high-elevation melt parameterization, avoids this overprediction.

Benjamin E. Smith↗
NASA NTRS2015

Antarctic Firn Compaction Rates from Repeat-Track Airborne Radar Data: I. Methods

While measurements of ice-sheet surface elevation change are increasingly used to assess mass change, the processes that control the elevation fluctuations not related to ice-flow dynamics (e.g. firn compaction and accumulation) remain difficult to measure. Here we use radar data from the Thwaites Glacier (West Antarctica) catchment to measure the rate of thickness change between horizons of constant age over different time intervals: 2009-10, 2010-11 and 2009-11. The average compaction rate to approximately 25m depth is 0.33ma(exp -1), with largest compaction rates near the surface. Our measurements indicate that the accumulation rate controls much of the spatio-temporal variations in the compaction rate while the role of temperature is unclear due to a lack of measurements. Based on a semi-empirical, steady-state densification model, we find that surveying older firn horizons minimizes the potential bias resulting from the variable depth of the constant age horizon. Our results suggest that the spatiotemporal variations in the firn compaction rate are an important consideration when converting surface elevation change to ice mass change. Compaction rates varied by up to 0.12ma(exp -1) over distances less than 6km and were on average greater than 20% larger during the 2010-11 interval than during 2009-10.

polar firn↗
NASA NTRS2001

Interannual Variations of Shallow Firn Temperature at Greenland Summit

Firn-temperature profiles are calculated in a thermal model using continuous surface temperatures derived from Automatic Weather Station (AWS) data and passive microwave data in the Greenland Summit region during the period 1987-1999. The results show that significant interannual variations of mean summer (June to August) and annual temperatures occur in the top 15 m, in addition to the normal seasonal cycle of firn temperature. At 5 m depth, the seasonal cycle is damped to 13% of the surface seasonal amplitude, but even at 15 m about 1% or 0.6 C of the seasonal cycle persists. Both summer and mean annual temperatures decrease from 1987 to 1992, followed by a general increasing trend. Interannual variability is 5 C at the surface, but only is only dampened to 3.2 C at 10 m depth and 0.7 C at 15 m depth. Dampening of the interannual variability with depth is slower than dampening of the seasonal cycle, because of the longer time constant of the interannual variation. The warmer spring and summer temperatures experienced in the top 5 m, due to both the seasonal cycle and interannual variations, affect the rate of firn densification, which is non-linearly dependent on temperature. During the 12 year period 1987-1999, the mean annual surface temperature is -29.2 C, and the mean annual 15 m temperature is -30. 1 C, which is more than 1 C warmer than a 15-m borehole temperature representing the period of about 1959 and warmer than the best-fit temperature history by Alley and Koci back to 1500 A.D..

Jun, Li↗
NASA NTRS2003

Determination of Mercury Content in a Shallow Firn Core from Summit, Greenland by Isotope Dilution Inductively Coupled Plasma Mass Spectrometry

The total mercury Hg content was determined in 6 cm sections of a near-surface 7 m firn core and in surrounding surface snow from Summit, Greenland (elevation: 3238 m, 72.58 N, 38.53 W) in May 2001 by isotope dilution cold-vapor inductively coupled plasma mass spectrometry (ID-CV-ICP-MS). The focus of this research was to evaluate the capability of the ID-CV-ICPMS technique for measuring trace levels of Hg typical of polar snow and firn. Highly enriched Hg-201 isotopic spike is added to approximately 10 ml melted core and thoroughly mixed. The Hg(+2) in the sample is reduced on line with tin (II) chloride (SnCl2) and the elemental Hg (Hg(0)) vapor pre-concentrated on to gold gauze using a commercial amalgam system. The Hg is then thermally desorbed and introduced into a quadrupole ICP-MS. The blank corrected Hg concentrations determined for all samples ranged from 0.25 ng/L to 1.74 ng/L (ppt) (average 0.59 ng/L plus or minus 0.28 ng/L) and fall within the range of those previously determined by Boutron et al., 1998 (less than or equal to 0.05 ng/L to 2.0 ng/L) for the Summit site. The average blank value was 0.19 ng/L plus or minus 0.045 ng/L (n=6). The Hg values specifically for the firn core range from 0.25 ng/L to 0.87 ng/L (average 0.51 ng/L plus or minus 0.13 ng/L) and show both values declining with time and larger variability in concentration in the top 1.8 m.

Mann, Jacqueline L.↗
NASA NTRS2014

Initial in Situ Measurements of Perennial Meltwater Storage in the Greenland Firn Aquifer

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.

measurements↗
NASA NTRS2013

Extensive Liquid Meltwater Storage in Firn Within the Greenland Ice Sheet

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.

meltwater↗