Changes in the firn structure of the Greenland Ice Sheet caused by recent warming

Introduction Conclusions References Tables Figures


Conclusions References
Tables Figures

Introduction
Investigations in the percolation zone of the Greenland Ice Sheet (GrIS) have revealed a highly variable snowpack structure characterized by the presence of ice lenses, pipes, and layers (Benson, 1962;Scott et al., 2006a;Parry et al., 2007;Harper et al., 2012).
The heterogeneous snowpack characteristic of this region results from periods of rel-Introduction

Conclusions References
Tables Figures

Back Close
Full atively high snow accumulation followed by short melt events during summer.During melt episodes, surface meltwater percolates through the snowpack and may refreeze at depth.The atmosphere has warmed considerably in the last decade over the GrIS (van den Broeke et al., 2009;Box et al., 2012;Bennartz et al., 2013), with 2010 and 2012 being the warmest years in western Greenland since records began (Tedesco et al., 2011(Tedesco et al., , 2013;;Bennartz et al., 2013;Tingley and Huybers, 2013).As a consequence, the area of the ice sheet covered by percolation facies has grown, and the amount of surface melt and subsequently refrozen meltwater retained in the firn has increased (Tedesco et al., 2008;Fettweis et al., 2011Fettweis et al., , 2013;;Harper et al., 2012;van Angelen et al., 2014).Extreme warming events such as in July 2012 (Nghiem et al., 2012;Bennartz et al., 2013) have further intensified melt, but it is unknown whether increasing meltwater production and subsequent percolation and refreezing at high elevations has shifted the equilibrium line higher and affected the buffering and transport of meltwater over the interior of the ice sheet.Furthermore, these processes and their high spatial and temporal variability have implications for altimetry-derived mass balance estimates.
Thus, continuous monitoring of densification, refreezing, and accumulation across the percolation zone (i.e.areas affected by significant melt and refreezing, but with little or no runoff so that all meltwater is refrozen and retained in the snowpack) is important for improving altimetry-derived mass balance estimates and for better understanding of the effects that increased seasonal melt across the ice sheet interior is having on the supra-glacial hydrology of the GrIS.
A model-based study in the early nineties (Pfeffer et al., 1991) showed how predictions of runoff-induced sea level rise from Greenland that did not consider meltwater refreezing within the firn could be overestimating sea level rise by as much as 5 cm over the next 150 years.The importance of meltwater retention was further highlighted by a study (Harper et al., 2012)  as a buffer between surface melt and runoff.Firn ice content was greater than 50 % by volume below 1600 m elevation, decreasing steadily with height.At the time, the presence of ice at an elevation of 2000 m was found to be relatively uncommon.That work added to a series of studies conducted between 2003 and 2006 as part of the first CryoSat Validation Experiment that helped assess near-surface snowpack and firn conditions at the higher end of the percolation zone (1950-2350 m).The CryoSat validation work focused on the region's spatially variable stratigraphy, as characterized by the presence of thin ice layers that were the main source of backscatter of Ku-band altimeter signals (Scott et al., 2006a;Parry et al., 2007;Helm et al., 2007).Together, these studies can provide a decadal record of ice content and can be linked with earlier data from NASA's Program for Arctic Regional Climate Assessment (PARCA, Abdalati et al., 1998).In the late 1990s, PARCA collected scores of shallow firn and ice cores to quantify spatial and temporal variability of annual accumulation rates over the GrIS (e.g.McConnell et al., 2000;Bales et al., 2001;Mosley-Thompson et al., 2001).
Here, we use field and remote sensing observations from the percolation zone of western Greenland in conjunction with output from regional climate model to (a) quantify changes in percolation conditions after the unusually warm years of 2010 and 2012, (b) identify areas where widespread percolation layers are found, and (c) assess the state and extent of the percolation zone of the GrIS given current melt trends over the ice sheet interior.For this, we characterized near-surface snowpack structure across a wide elevation range on the western slope of the GrIS, and estimated the total ice content and area covered by percolation facies resulting from melting and refreezing patterns.Our field sites are all located within one degree of latitude (69.2-71.1 • N) on the western slope of the ice sheet, spanning an elevation range from 1900 to 2500 m (Fig. 1).At each site, snowpits were excavated in April of 2011April of , 2012April of , 2013April of , and 2014 in order to characterize regional near-surface snowpack conditions and percolation facies in the region following the melt season from the previous year.Extensive melt layers are identified and traced using airborne radar to identify areas covered by percolation facies, and airborne laser altimetry data are used to estimate annual eleva-Introduction

Conclusions References
Tables Figures

Back Close
Full tion changes in the region.Additionally, output from the Regional Atmospheric Climate Model (RACMO2.3/GR,van Meijgaard et al., 2008;Ettema et al., 2009;van den Broeke, 2009) are used to assess, at a regional scale, the intensity and extent of the abnormally strong 2010 and 2012 melt seasons in Greenland's accumulation zone.

Observations
We measured firn ice content in snow pits during April (pre-summer melt) field campaigns between 2011 and 2014.The 2011 campaign was part of CryoSat-2 validation activities continuing efforts initiated in 2004 and 2006 (Parry et al., 2007;Scott et al., 2006a, b;Helm et al., 2007).Near-surface density and stratigraphy measurements were obtained from snow pits at 4 sites in the upper end of the percolation zone, located between 2350 and 2490 m along the EGIG line (named T12-T15 hereafter for historical reasons).The April 2012, 2013 and 2014 campaigns were part of a NASAfunded study of outlet glacier dynamics and also included snow pit surveys of winter accumulation and ice layer thickness at four sites in the catchment of Jakobshavn Glacier (hereafter referred as J1-J4 sites) between elevations of 1935 and 2350 m (Fig. 1).
Additionally, we extended the firn-ice content record in this area back in time to 1977 using 3 shallow firn cores drilled in 1998 at 2000, 2500 m, and 2795 m elevation (69 We used remote sensing observations acquired by NASA's Operation IceBridge Snow Radar and the Airborne Topographic Mapper (ATM).The Snow Radar is a wideband radar that operates in the S-C bands and can map polar firn internal layers with high spatial resolution (Leuschen et al., 2014), and is used to trace the ice layers beneath the winter accumulation in 2011 by identifying internal radar reflections caused by melt layers.The ATM is a scanning laser altimeter used in this study to estimate annual elevation change between 2010 and 2013 along a frequently repeated transect in western Greenland (Fig. 1).Introduction

Conclusions References
Tables Figures

Back Close
Full Output data from the RACMO2.3/GRmodel with a horizontal resolution of ∼ 11 km (Ettema et al., 2009;van den Broeke, 2009) are used in combination with the remote sensing observations to map the percolation zone extent.RACMO2.3/GR is coupled to a physical snow model that treats surface albedo as a function of melt as well as percolation and refreezing, making the model results more realistic (Bougamont et al., 2005;Van Angelen et al., 2014).A recent study comparing RACMO2.3/GRresults with QuickScat melt rates shows that discrepancies between the model and satellite data occur mainly in areas of high meltwater runoff (Fettweis et al., 2011), due in part to the lack of remote sensing capability of assessing melt and runoff rates.There is good agreement at higher elevations, which suggests that assumptions about snowpack heat transfer and energy balance in the model are well parameterized for the percolation zone of the GrIS.Thus, the model results used here are limited to regions of no meltwater runoff.

Snowpack structure
Total annual ice content found beneath the wintertime snow accumulation during each of the campaigns described in the previous section is summarized in Fig. 2, revealing an increase in the later years, where the annual ice layers were an order of magnitude thicker in total than the 1977-1997 average from PARCA core 6945.Figure 3a et al., 2001).Although the thickness of each annual layer (after being converted to water equivalent using density) can be ascertained using any one of these parameters, for this study the layer thicknesses (net accumulation) were calculated using the winter minima in dust concentration.The thickness of the different annual layers is affected Introduction

Conclusions References
Tables Figures

Back Close
Full by annual snowfall, deflation and re-deposition by wind, compaction, and melting, and varies on the order of tens of cm from year to year in western Greenland (McConnell et al., 2000;de la Peña et al., 2010;Burgess et al., 2010).Small amounts of ice are commonly found in the core, but are not present on an annual basis.While no ice was observed in core 6945 over the periods of 1977-1979 and 1996-1997, between 1987 and 1991 several thin ice layers were found separated by a few cm.Most layers were found to be 3-4 cm thick, and the thickest (1989) was only 9 cm thick.Cores 6943 and 6941 show similar patterns, with occasional ice layers 1-2 cm thick.
Current conditions are represented in the schematic (Fig. 3b), made from measurements of winter snow depth and ice layer thickness for each site visited during spring in 2013 and 2014 and illustrating the near-surface stratigraphy as found in spring 2014.Ice content measured in 2011 at T12, a site located at the same altitude as J4 but roughly 10 km north is included as well.The stratigraphy of the shallow firn cores shows that at least between 1977 and 1997, melt events at higher elevations were more rare and much less intensive relative to those observed since.The extreme melt events of 2012 created conditions that facilitated the formation of impermeable ice layers several times thicker than previously observed across this elevation in the percolation zone of the GrIS.These layers were found at all sites visited and appear to be continuous throughout the area surveyed.In 2004, total ice thickness at an elevation of 1950 m averaged 10 cm (Parry et al., 2007), significantly less than the ice content in 2013 at J4, located 400 m higher on the ice sheet.The total mass of the ice layers found at the J sites averaged 441±12 kg m −2 and 306±118 kg m the T12-T15 sites was 1.10 ± 0.12 m.Winter snow density was measured in 2011 and exhibited little variability (247 ± 8 kg m −3 ), consistent with previous observations (Parry et al., 2007).Measurements in April 2011 by the NASA Snow Radar show the extent of the type of percolation features as described above that formed during the 2010 melt season (Fig. 4).The ∼ 420 km-long transect extends from the ablation zone to the ice divide in the dry snow zone.Radar signals are partially backscattered from within a stratified snowpack by abrupt changes encountered in snow density and/or ice structure, such as ice layers in the percolation zone, or "autumn hoar" in the dry snow zone.Near surface layering is observed in the 1600-2200 m elevation range and winter accumulation over the previous summer melt layer is clearly resolved while deeper layers are obscured by infiltration ice that limits radar signal penetration.The topmost reflection under the observed winter surface is continuous over the percolation zone, confirming that ice layers observed at the J field sites are widespread over an elevation range of 1600 to 2200 m, about 220 km inland from the ice margin.The ice layer was retraced with a custom-made threshold algorithm configured to identify strong reflections underneath the surface (shown with a black line in Fig. 4).The algorithm tracks continuity between horizontally adjacent pixels, so that the retraced layer in one individual radar acquisition is not separated vertically by more than ∼ 20 cm from the next measurement.In some sections the buried signal could not be differentiated from the surface, and at some points there is discontinuity in the signal.However, for most of the percolation zone the reflection appears continuous, tracked at a depth of 0.81 ± 0.29 m, slightly lower than the field measurements at T sites.The underestimation is likely the result of tracking on the leading edge of the signal and does not affect our analysis.
At higher elevations, annual accumulation layers are clearly seen to depths of at least 15 m, consistent with previous observations from the dry snow zone (Hawley et al., 2005;de la Peña et al., 2010;Simonsen et al., 2013).While melt has intensified over the last decade, the 2012 melt episodes had notably bigger impacts on the firn structure, with significant melting and infiltration extending to the ice divide (Tedesco Introduction

Conclusions References
Tables Figures

Back Close
Full ), but as discussed in the previous section, total ice content increased more than 5 times over the same period.Ice layer thickness does not increase as a simple function of surface melt.Infiltration of meltwater into firn is complicated by factors including the initial thermal state and structure of the firn and the timing and duration of melt.For example, more frequent extreme melt events, such as occurred in 2012, may cause rapid snow saturation and refreezing near the surface, resulting in the formation of thicker and more extensive ice layers.The formation of ice reservoirs of this magnitude is thus not only dependent on melt increase, but also on how melt is distributed during summer.Uncertainties in melt and refreezing rates in RACMO2.3/GRadd to this complexity, and as a result, it is difficult to confidently predict firn changes using only the model results.For instance, in 2003, the ice layers formed during the summer at 2000 m elevation accounted for 29 % of total modeled melt, while Introduction

Conclusions References
Tables Figures

Back Close
Full in 2012 discernible ice content was 79 % of the modeled melt (52 % at the highest site at 2350 m).The 2012 melt season anomaly was driven partly by extreme, short-lived melt episodes (i.e.Nghiem et al., 2012;Bennartz et al., 2013;Tedesco et al., 2013), which as discussed may result in more ice content.Increases in modeled melt rates at higher sites are even more pronounced than at J1, with rates at J4 increasing from less than 100 kg m −2 over the period from 1958 to 1999, to a rate of 305 and 580 kg m −2 in 2010 and 2012 respectively (Fig. 6).Prior to 2012, the largest modeled melt rates occurred in 2010 at all sites except at the highest regions (above 2200 m), which registered the largest melt rate in 2007 (second highest for the lower sites).Large modeled melt rates at all sites occurred during the warm years of 1995, 1998, 2007, and 2010.
Melt at J1 increased from less than 300 kg m −2 between 1958 and 1999, to 340 kg m −2 in 2003 and to 755 kg m −2 in 2010.At J4 there was only an increase of 30 % in 2010 with respect to 2003 (from 235 to 305 kg m −2 ), but this increased to 580 kg m −2 in 2012.
Before 2012, the modeled melt rate at this site was highest in 2007, in contrast to the 2010 peak melt rates estimated for most of the western Greenland margin (Tedesco et al., 2011).The mass of ice found in 2013 at sites J1 and J2 exceeds the mass of the total annual accumulation for 2012, a trend that according to RACMO2.3/GRsuggests has been occurring at J1 since 2007.Modeled melt rates in 2013 are significantly lower than in 2012, but still larger than the 1958-1999 average.We found no evidence of meltwater pathways or pond formation at our field sites.If, however, melt continues to exceed total accumulation, it is likely that lateral supraglacial transport and eventual runoff of meltwater will reach to 2000 m elevation.in 2010 the percolation zone covered a total area of 486 000 km 2 , only a 5 % increase with respect to 2003, the intensity of melt increased drastically after 1999, especially in west and southwest Greenland.To quantify the magnitude of refrozen meltwater in 2012, we estimate total annual ice content in western Greenland by identifying areas where snowfall and melt rates are sufficiently high to form the stratified snowpack structure described above.We estimate the total mass refrozen in the snowpack in western Greenland by relating melt rates from RACMO2/GR to the measured ice content found during 2011 and 2013 (a percentage of the total melt), and subtract the mass of the snow that occupied the volume where meltwater refroze (total ice content cannot exceed the total winter accumulation).Based on the field measurements, we assume an average winter snow density of 250 kg m −3 .Our estimates indicate that the total 2012 infiltration ice content above the equilibrium line and west of the ice divide of the GrIS was 74 ± 25 Gt, compared to an average of only 13 ± 2 Gt between 1958 and 1999.
Most of this ice is distributed near the surface in the form of thick layers that appear to be widespread.

Discussion and conclusions
Our field observations reveal substantial, recent changes in the structure of highelevation percolation facies that differ from those recorded in previous decades.Observations in 2004 and 2006 showed an increase from the 1977-1997 average ice content, but the lack of significant ice above 2000 m was evidence that the transition between the percolation and dry snow zones was near this elevation (Scott et al., 2006b).In contrast, the Snow Radar data obtained in spring 2011 reveals that the percolation zone has risen with the transition in the previous few years occurring at elevations between 2200 and 2400 m.More importantly, the concentration of ice at higher elevations has increased dramatically, and the formation of thick, widespread continuous ice layers is increasingly common over extensive areas of the percolation zone causing large changes in the air content of the firn.Furthermore, the increase in

Conclusions References
Tables Figures

Back Close
Full firn ice content raises questions regarding the permeability of the firn column across the percolation zone, and the effects it could have on meltwater retention and transport and on firn compaction rates.The formation of percolation facies in the last few decades of the 20th Century above 2000 m was limited, and during some years melt was not even strong enough to form identifiable ice layers.In this context, the percolation features described here have no recent precedent.Considering, as measured in the field, an average surface snow density before melt of 250 kg m −3 , and that meltwater refreezes in a volume occupied by snow of this density, each meter of winter snow would equate to layers of ice totaling ∼ 0.35 m thick after refreezing if little or no air remains trapped within.This process will create the observed seasonal change in the volume of the firn governed by accumulation, compaction, and melting and refreezing.Even without thermal snow densification, commonly used steady-state snow density and firn compaction assumptions used for altimetry-derived mass balance estimates would need to be reconsidered for regions where percolation features become more impermeable.These ice layers not only significantly limit the total meltwater buffering capacity of the percolation zone estimated by Harper et al. (2012), but also introduce variability in compaction rates if air is trapped within ice layers.Moreover, and as observed at sites J1 and J2, the winter snow in these regions of the accumulation zone of the ice sheet completely melted and refroze in 2012, meaning mass gain in these areas was in the form of superimposed ice.
Melt affecting the GrIS has increased over recent decades, with pronounced departures from the 1958-1999 mean melt rate during the last few years in regions well within the accumulation zone of the GrIS.High (and highly variable) accumulation rates play an important role in the surface mass balance of the ice sheet at higher elevations, but with melt rates increasing faster than snow accumulation in Greenland (Fettweis, 2007;Van Angelen et al., 2014), meltwater excess during warm years may saturate the already limited firn meltwater storage capacity over large areas (Van Angelen et al., 2013).Although no evidence of extensive lateral hydrological pathways was found at any of our field sites, it is likely that if melt continues to exceed total accumulation, Introduction

Conclusions References
Tables Figures

Back Close
Full water will soon be transported to lower elevations by supraglacial flow, especially during years experiencing extreme melt.Recent studies show that lakes are increasingly found at higher elevations (Banwell et al., 2012;Liang et al., 2012;Howat et al., 2013), and it is unclear how the percolation zone will evolve if more meltwater ponds are being formed over firn.Furthermore, if current warming trends continue, percolation facies in the future will cover an even greater extent potentially extending all areas above the equilibrium line.Regardless of the role they will play in the future, these huge ice reservoirs are the result of an intense melting process of the same order of magnitude as the total mass imbalance of the GrIS, and the consequences of their formation described here underlines the importance of monitoring the evolution of Greenland's firn layer in the coming years.Introduction

Conclusions References
Tables Figures

Back Close
Full  Tedesco, M., Fettweis, X., Mote, T., Wahr, J., Alexander, P., Box, J. E., and Wouters, B.: Ev-idence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data, The Cryosphere, 7, 615-630, doi:10.5194/tc-7-615-2013Cryosphere, 7, 615-630, doi:10.5194/tc-7-615- , 2013.Tingley, M. P. and Huybers, P.: Recent temperature extremes at high northern latitudes unprecedented in the past 600 years, Nature, 496, 201-205, doi:10.1038Nature, 496, 201-205, doi:10. /nature11969, 2013. van . van  Full  Full Discussion Paper | Discussion Paper | Discussion Paper | based on field measurements obtained in 2007, 2008 and 2009 along the Expédition Glaciologique Internationale au Groenland (EGIG) line in western Greenland.Harper et al. (2012) estimated that Greenland's firn has the potential to store between 322 and 1289 Gt of meltwater, confirming its importantce Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | illustrates the ice layers and lenses present in each core's visible stratigraphy.Cores 6945 and 6943 (named for the coordinates where they were acquired) are ∼ 18 m long and extend back to 1976 while core 6941 is 11.7 m long and extends to a depth equating to 1985's accumulation.The cores were dated using a combination of the winter minima in the seasonal variations of dust concentration, δ 18 O, and H 2 O 2 (Mosley-Thompson Discussion Paper | Discussion Paper | Discussion Paper | −2 in 2013 (related to 2012 melt season) and 2014 (related to 2013 melt season) respectively.In April 2011 at the higher elevation sites, the total ice layer mass was 285, 255 and 160 kg m −2 for T12, T13 and T14 respectively.Ice content decreased at T15 (h = 2490 m), where an ice layer just 4 cm thick (36 kg m −2 ) was found, suggesting this site was close to the boundary of the dry snow zone that year.At the J-sites, winter accumulation was measured at 1.11 ± 0.17 m, 0.885 ± 0.08 m, and 1.35 ± 0.07 m in the winters of 2011-2012, 2012-2013 and 2013-2014 respectively.The 2010-2011 winter accumulation measured at Discussion Paper | Discussion Paper | Discussion Paper | Figure6presents results from RACMO2.3/GR between 1958 and 2013 for each of the sites visited between 2011 and 2014, and Table1shows at each site the modeled melt and accumulation rates at each site for the years 2003, 2010, and 2012.Melt at J1 increased from 340 kg m −2 in 2003 to 810 kg m −2 in 2012 (3 times the 1958-1999 average of 285 kg m −2 ), but as discussed in the previous section, total ice content increased more than 5 times over the same period.Ice layer thickness does not increase as a simple function of surface melt.Infiltration of meltwater into firn is complicated by factors including the initial thermal state and structure of the firn and the timing and duration of melt.For example, more frequent extreme melt events, such as Discussion Paper | Discussion Paper | Discussion Paper |

Figure 7
presents annual melt predicted by RACMO2.3/GR in areas of the GrIS with predicted positive surface mass balance (i.e. the accumulation zone).The left figure shows the 1958-1999 average, while the right figure is the annual melt rates in 2012, illustrating the changes in melt intensity and the size of the areas affected by melt over Greenland.The predicted average area of the accumulation zone that experienced melt between 1958 and 1999 was about 405 000 km 2 .The extent of this area then increased from 465 000 km 2 in 2003 to a maximum of 1 155 000 km 2 in 2012.Although Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 1 .Figure 2 .Figure 4 .Figure 5 .Figure 6 .Figure 7 .
Figure 1.Map of Greenland showing the research area and research sites (all T-sites are along the EGIG line).Black line is the flightpath of the NASA Operation IceBridge data shown in Figs. 4 and 5. Snowpit sites are shown with an X, while core sites are shown with a blue circle.

Table 1 .
Melt and accumulation at each field site as simulated by RACMO2.3/GR.