ALBERT

Description: Tidally induced velocity variations of the Beardmore Glacier, Antarctica, and their representation in satellite measurements of ice velocity The Cryosphere, 7, 1375-1384, 2013 Author(s): O. J. Marsh, W. Rack, D. Floricioiu, N. R. Golledge, and W. Lawson Ocean tides close to the grounding line of outlet glaciers around Antarctica have been shown to directly influence ice velocity, both linearly and non-linearly. These fluctuations can be significant and have the potential to affect satellite measurements of ice discharge, which assume displacement between satellite passes to be consistent and representative of annual means. Satellite observations of horizontal velocity variation in the grounding zone are also contaminated by vertical tidal effects, the importance of which is highlighted here in speckle tracking measurements. Eight TerraSAR-X scenes from the grounding zone of the Beardmore Glacier are analysed in conjunction with GPS measurements to determine short-term and decadal trends in ice velocity. Diurnal tides produce horizontal velocity fluctuations of 〉50% on the ice shelf, recorded in the GPS data 4 km downstream of the grounding line. This variability decreases rapidly to 〈5% only 15 km upstream of the grounding line. Daily fluctuations are smoothed to 〈1% in the 11-day repeat pass TerraSAR-X imagery, but fortnightly variations over this period are still visible and show that satellite-velocity measurements can be affected by tides over longer periods. The measured tidal displacement observed in radar look direction over floating ice also allows the grounding line to be identified, using differential speckle tracking where phase information cannot be easily unwrapped.

Description: Snow on the Ross Ice Shelf: comparison of reanalyses and observations from automatic weather stations The Cryosphere, 7, 1399-1410, 2013 Author(s): L. Cohen and S. Dean Snow accumulation measurements from automatic weather stations (AWS) around the Ross Ice Shelf (RIS), Antarctica, are used to provide a new set of ground-based observations which are compared to precipitation from the ECMWF ERA-Interim and NCEP/NCAR Reanalysis-2 datasets. The high temporal resolution of the AWS snow accumulation measurements allow for an event-based comparison of reanalyses precipitation to the in situ observations. Snow accumulation records from nine AWS provide multiple years of accumulation data between 2008 and 2012 over a relatively large, homogeneous region of Antarctica, and also provide the basis for a statistical evaluation of accumulation and precipitation events. The complex effects of wind on snow accumulation (which can both limit and enhance accumulation) complicate the use of the accumulation measurements, but this analysis shows that they can provide a valuable source of ground-based observations for comparisons to modelled precipitation on synoptic timescales. The analysis shows that ERA-Interim reproduces more precipitation events than NCEP-2, and these events correspond to an average 8.2% more precipitation. Significant correlations between reanalyses and AWS event sizes are seen at several stations and show that ERA-Interim consistently produces larger precipitation events than NCEP-2.

Description: A general treatment of snow microstructure exemplified by an improved relation for thermal conductivity The Cryosphere, 7, 1473-1480, 2013 Author(s): H. Löwe, F. Riche, and M. Schneebeli Finding relevant microstructural parameters beyond density is a longstanding problem which hinders the formulation of accurate parameterizations of physical properties of snow. Towards a remedy, we address the effective thermal conductivity tensor of snow via anisotropic, second-order bounds. The bound provides an explicit expression for the thermal conductivity and predicts the relevance of a microstructural anisotropy parameter Q , which is given by an integral over the two-point correlation function and unambiguously defined for arbitrary snow structures. For validation we compiled a comprehensive data set of 167 snow samples. The set comprises individual samples of various snow types and entire time series of metamorphism experiments under isothermal and temperature gradient conditions. All samples were digitally reconstructed by micro-computed tomography to perform microstructure-based simulations of heat transport. The incorporation of anisotropy via Q considerably reduces the root mean square error over the usual density-based parameterization. The systematic quantification of anisotropy via the two-point correlation function suggests a generalizable route to incorporate microstructure into snowpack models. We indicate the inter-relation of the conductivity to other properties and outline a potential impact of Q on dielectric constant, permeability and adsorption rate of diffusing species in the pore space.

Description: Influence of high-order mechanics on simulation of glacier response to climate change: insights from Haig Glacier, Canadian Rocky Mountains The Cryosphere, 7, 1527-1541, 2013 Author(s): S. Adhikari and S. J. Marshall Evolution of glaciers in response to climate change has mostly been simulated using simplified dynamical models. Because these models do not account for the influence of high-order physics, corresponding results may exhibit some biases. For Haig Glacier in the Canadian Rocky Mountains, we test this hypothesis by comparing simulation results obtained from 3-D numerical models that deal with different assumptions concerning physics, ranging from simple shear deformation to comprehensive Stokes flow. In glacier retreat scenarios, we find a minimal role of high-order mechanics in glacier evolution, as geometric effects at our site (the presence of an overdeepened bed) result in limited horizontal movement of ice (flow speed on the order of a few meters per year). Consequently, high-order and reduced models all predict that Haig Glacier ceases to exist by ca. 2080 under ongoing climate warming. The influence of high-order mechanics is evident, however, in glacier advance scenarios, where ice speeds are greater and ice dynamical effects become more important. Although similar studies on other glaciers are essential to generalize such findings, we advise that high-order mechanics are important and therefore should be considered while modeling the evolution of active glaciers. Reduced model predictions may be adequate for other glaciologic and topographic settings, particularly where flow speeds are low and where mass balance changes dominate over ice dynamics in determining glacier geometry.

Description: Meteorological drivers of ablation processes on a cold glacier in the semi-arid Andes of Chile The Cryosphere, 7, 1513-1526, 2013 Author(s): S. MacDonell, C. Kinnard, T. Mölg, L. Nicholson, and J. Abermann Meteorological and surface change measurements collected during a 2.5 yr period are used to calculate surface mass and energy balances at 5324 m a.s.l. on Guanaco Glacier, a cold-based glacier in the semi-arid Andes of Chile. Meteorological conditions are marked by extremely low vapour pressures (annual mean of 1.1 hPa), strong winds (annual mean of 10 m s −1 ), shortwave radiation receipt persistently close to the theoretical site maximum during cloud-free days (mean annual 295 W m −2 ; summer hourly maximum 1354 W m −2 ) and low precipitation rates (mean annual 45 mm w.e.). Snowfall occurs sporadically throughout the year and is related to frontal events in the winter and convective storms during the summer months. Net shortwave radiation provides the greatest source of energy to the glacier surface, and net longwave radiation dominates energy losses. The turbulent latent heat flux is always negative, which means that the surface is always losing mass via sublimation, which is the main form of ablation at the site. Sublimation rates are most strongly correlated with net shortwave radiation, incoming shortwave radiation, albedo and vapour pressure. Low glacier surface temperatures restrict melting for much of the period, however episodic melting occurs during the austral summer, when warm, humid, calm and high pressure conditions restrict sublimation and make more energy available for melting. Low accumulation (131 mm w.e. over the period) and relatively high ablation (1435 mm w.e.) means that mass change over the period was negative (−1304 mm w.e.), which continued the negative trend recorded in the region over the last few decades.

Description: Antarctic ice-mass balance 2003 to 2012: regional reanalysis of GRACE satellite gravimetry measurements with improved estimate of glacial-isostatic adjustment based on GPS uplift rates The Cryosphere, 7, 1499-1512, 2013 Author(s): I. Sasgen, H. Konrad, E. R. Ivins, M. R. Van den Broeke, J. L. Bamber, Z. Martinec, and V. Klemann We present regional-scale mass balances for 25 drainage basins of the Antarctic Ice Sheet (AIS) from satellite observations of the Gravity and Climate Experiment (GRACE) for time period January 2003 to September 2012. Satellite gravimetry estimates of the AIS mass balance are strongly influenced by mass movement in the Earth interior caused by ice advance and retreat during the last glacial cycle. Here, we develop an improved glacial-isostatic adjustment (GIA) estimate for Antarctica using newly available GPS uplift rates, allowing us to more accurately separate GIA-induced trends in the GRACE gravity fields from those caused by current imbalances of the AIS. Our revised GIA estimate is considerably lower than previous predictions, yielding an estimate of apparent mass change of 53 ± 18 Gt yr −1 . Therefore, our AIS mass balance of −114 ± 23 Gt yr −1 is less negative than previous GRACE estimates. The northern Antarctic Peninsula and the Amundsen Sea sector exhibit the largest mass loss (−26 ± 3 Gt yr −1 and −127 ± 7 Gt yr −1 , respectively). In contrast, East Antarctica exhibits a slightly positive mass balance (26 ± 13 Gt yr −1 ), which is, however, mostly the consequence of compensating mass anomalies in Dronning Maud and Enderby Land (positive) and Wilkes and George V Land (negative) due to interannual accumulation variations. In total, 6% of the area constitutes about half the AIS imbalance, contributing 151 ± 7 Gt yr −1 (ca. 0.4 mm yr −1 ) to global mean sea-level change. Most of this imbalance is caused by ice-dynamic speed-up expected to prevail in the near future.

Description: Evidence of meltwater retention within the Greenland ice sheet The Cryosphere, 7, 1433-1445, 2013 Author(s): A. K. Rennermalm, L. C. Smith, V. W. Chu, J. E. Box, R. R. Forster, M. R. Van den Broeke, D. Van As, and S. E. Moustafa Greenland ice sheet mass losses have increased in recent decades with more than half of these attributed to surface meltwater runoff. However, the magnitudes of englacial storage, firn retention, internal refreezing and other hydrologic processes that delay or reduce true water export to the global ocean remain less understood, partly due to a scarcity of in situ measurements. Here, ice sheet surface meltwater runoff and proglacial river discharge between 2008 and 2010 near Kangerlussuaq, southwestern Greenland were used to establish sub- and englacial meltwater storage for a small ice sheet watershed (36–64 km 2 ). This watershed lacks significant potential meltwater storage in firn, surface lakes on the ice sheet and in the proglacial area, and receives limited proglacial precipitation. Thus, ice sheet surface runoff not accounted for by river discharge can reasonably be attributed to retention in sub- and englacial storage. Evidence for meltwater storage within the ice sheet includes (1) characteristic dampened daily river discharge amplitudes relative to ice sheet runoff; (2) three cold-season river discharge anomalies at times with limited ice sheet surface melt, demonstrating that meltwater may be retained up to 1–6 months; (3) annual ice sheet watershed runoff is not balanced by river discharge, and while near water budget closure is possible as much as 54% of melting season ice sheet runoff may not escape to downstream rivers; (4) even the large meltwater retention estimate (54%) is equivalent to less than 1% of the ice sheet volume, which suggests that storage in en- and subglacial cavities and till is plausible. While this study is the first to provide evidence for meltwater retention and delayed release within the Greenland ice sheet, more information is needed to establish how widespread this is along the Greenland ice sheet perimeter.

Description: High-resolution provenance of desert dust deposited on Mt. Elbrus, Caucasus in 2009–2012 using snow pit and firn core records The Cryosphere, 7, 1481-1498, 2013 Author(s): S. Kutuzov, M. Shahgedanova, V. Mikhalenko, P. Ginot, I. Lavrentiev, and S. Kemp The first record of dust deposition events on Mt. Elbrus, Caucasus Mountains derived from a snow pit and a shallow firn core is presented for the 2009–2012 period. A combination of isotopic analysis, SEVIRI red-green-blue composite imagery, MODIS atmospheric optical depth fields derived using the Deep Blue algorithm, air mass trajectories derived using the HYSPLIT model and analyses of meteorological data enabled identification of dust source regions with high temporal (hours) and spatial (ca. 20–100 km) resolution. Seventeen dust deposition events were detected; fourteen occurred in March–June, one in February and two in October. Four events originated in the Sahara, predominantly in northeastern Libya and eastern Algeria. Thirteen events originated in the Middle East, in the Syrian Desert and northern Mesopotamia, from a mixture of natural and anthropogenic sources. Dust transportation from Sahara was associated with vigorous Saharan depressions, strong surface winds in the source region and mid-tropospheric southwesterly flow with daily winds speeds of 20–30 m s −1 at 700 hPa level. Although these events were less frequent than those originating in the Middle East, they resulted in higher dust concentrations in snow. Dust transportation from the Middle East was associated with weaker depressions forming over the source region, high pressure centred over or extending towards the Caspian Sea and a weaker southerly or southeasterly flow towards the Caucasus Mountains with daily wind speeds of 12–18 m s −1 at 700 hPa level. Higher concentrations of nitrates and ammonium characterised dust from the Middle East deposited on Mt. Elbrus in 2009 indicating contribution of anthropogenic sources. The modal values of particle size distributions ranged between 1.98 μm and 4.16 μm. Most samples were characterised by modal values of 2.0–2.8 μm with an average of 2.6 μm and there was no significant difference between dust from the Sahara and the Middle East.

Description: Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: evaluation by remote sensing and modeling of the morpho-topographic and climate controls The Cryosphere, 7, 1455-1471, 2013 Author(s): A. Rabatel, A. Letréguilly, J.-P. Dedieu, and N. Eckert We present time series of equilibrium-line altitude (ELA) measured from the end-of-summer snow line altitude computed using satellite images, for 43 glaciers in the western Alps over the 1984–2010 period. More than 120 satellite images acquired by Landsat, SPOT and ASTER were used. In parallel, changes in climate variables, summer cumulative positive degree days (CPDD) and winter precipitation, were analyzed over the same time period using 22 weather stations located inside and around the study area. Assuming a continuous linear trend over the study period: (1) the average ELA of the 43 glaciers increased by about 170 m; (2) summer CPDD increased by about 150 PDD at 3000 m a.s.l.; and (3) winter precipitation remained rather stationary. Summer CPDD showed homogeneous spatial and temporal variability; winter precipitation showed homogeneous temporal variability, but some stations showed a slightly different spatial pattern. Regarding ELAs, temporal variability between the 43 glaciers was also homogeneous, but spatially, glaciers in the southern part of the study area differed from glaciers in the northern part, mainly due to a different precipitation pattern. A sensitivity analysis of the ELAs to climate and morpho-topographic variables (elevation, aspect, latitude) highlighted the following: (1) the average ELA over the study period of each glacier is strongly controlled by morpho-topographic variables; and (2) the interannual variability of the ELA is strongly controlled by climate variables, with the observed increasing trend mainly driven by increasing temperatures, even if significant nonlinear, low-frequency fluctuations appear to be driven by winter precipitation anomalies. Finally, we used an expansion of Lliboutry's approach to reconstruct fluctuations in the ELA of any glacier of the study area with respect to morpho-topographic and climate variables, by quantifying their respective weight and the related uncertainties in a consistent manner within a hierarchical Bayesian framework. This method was tested and validated using the ELA measured on the satellite images.

Description: A note on the water budget of temperate glaciers The Cryosphere, 7, 1557-1564, 2013 Author(s): J. Oerlemans In this note, the total dissipative melting in temperate glaciers is studied. The analysis is based on the notion that the dissipation is determined by the loss of potential energy due to the downward motion of mass (ice, snow, meltwater and rain). A mathematical formulation of the dissipation is developed and applied to a simple glacier geometry. In the next step, meltwater production resulting from enhanced ice motion during a glacier surge is calculated. The amount of melt energy available follows directly from the lowering of the centre of gravity of the glacier. To illustrate the concept, schematic calculations are presented for a number of glaciers with different geometric characteristics. Typical dissipative melt rates, expressed as water-layer depth averaged over the glacier, range from a few centimetres per year for smaller glaciers to half a metre per year for Franz Josef Glacier, one of the most active glaciers in the world (in terms of mass turnover). The total generation of meltwater during a surge is typically half a metre. For Variegated Glacier a value of 70 cm is found, for Kongsvegen 20 cm. These values refer to water layer depth averaged over the entire glacier. The melt \textit{rate} depends on the duration of the surge. It is generally an order of magnitude greater than water production by `normal' dissipation. On the other hand, the additional basal melt rate during a surge is comparable in magnitude with the water input from meltwater and precipitation. This suggests that enhanced melting during a surge does not grossly change the total water budget of a glacier. Basal water generated by enhanced sliding is an important ingredient in many theories of glacier surges. It provides a positive feedback mechanism that actually makes the surge happen. The results found here suggest that this can only work if water generated by enhanced sliding accumulates in a part of the glacier base where surface meltwater and rain have no or very limited access. This finding seems compatible with the fact that, on many glaciers, surges are initiated in the lower accumulation zone.

Description: Global glacier changes: a revised assessment of committed mass losses and sampling uncertainties The Cryosphere, 7, 1565-1577, 2013 Author(s): S. H. Mernild, W. H. Lipscomb, D. B. Bahr, V. Radić, and M. Zemp Most glaciers and ice caps (GIC) are out of balance with the current climate. To return to equilibrium, GIC must thin and retreat, losing additional mass and raising sea level. Because glacier observations are sparse and geographically biased, there is an undersampling problem common to all global assessments. Here, we further develop an assessment approach based on accumulation-area ratios (AAR) to estimate committed mass losses and analyze the undersampling problem. We compiled all available AAR observations for 144 GIC from 1971 to 2010, and found that most glaciers and ice caps are farther from balance than previously believed. Accounting for regional and global undersampling errors, our model suggests that GIC are committed to additional losses of 32 ± 12% of their area and 38 ± 16% of their volume if the future climate resembles the climate of the past decade. These losses imply global mean sea-level rise of 163 ± 69 mm, assuming total glacier volume of 430 mm sea-level equivalent. To reduce the large uncertainties in these projections, more long-term glacier measurements are needed in poorly sampled regions.

Description: Recent changes in spring snowmelt timing in the Yukon River basin detected by passive microwave satellite data The Cryosphere, 7, 905-916, 2013 Author(s): K. A. Semmens and J. M. Ramage Spring melt is a significant feature of high latitude snowmelt dominated drainage basins influencing hydrological and ecological processes such as snowmelt runoff and green-up. Melt duration, defined as the transition period from snowmelt onset until the end of the melt refreeze, is characterized by high diurnal amplitude variations (DAV) where the snowpack is melting during the day and refreezing at night, after which the snowpack melts constantly until depletion. Determining trends for this critical period is necessary for understanding how the Arctic is changing with rising temperatures and provides a baseline from which to assess future change. To study this dynamic period, brightness temperature ( T b ) data from the Special Sensor Microwave Imager (SSM/I) 37 V-GHz frequency from 1988 to 2010 were used to assess snowmelt timing trends for the Yukon River basin, Alaska/Canada. Annual T b and DAV for 1434 Equal-Area Scalable Earth (EASE)-Grid pixels (25 km resolution) were processed to determine melt onset and melt refreeze dates from T b and DAV thresholds previously established in the region. Temporal and spatial trends in the timing of melt onset and melt refreeze, and the duration of melt were analyzed for the 13 sub-basins of the Yukon River basin with three different time interval approaches. Results show a lengthening of the melt period for the majority of the sub-basins with a significant trend toward later end of melt refreeze after which the snowpack melts day and night leading to snow clearance, peak discharge, and green-up. Earlier melt onset trends were also found in the higher elevations and northernmost sub-basins (Porcupine, Chandalar, and Koyukuk rivers). Latitude and elevation displayed the dominant controls on melt timing variability and spring solar flux was highly correlated with melt timing in middle (∼600–1600 m) elevations.

Description: Seasonal controls on snow distribution and aerial ablation at the snow-patch and landscape scales, McMurdo Dry Valleys, Antarctica The Cryosphere, 7, 917-931, 2013 Author(s): J. W. Eveland, M. N. Gooseff, D. J. Lampkin, J. E. Barrett, and C. D. Takacs-Vesbach Accumulated snow in the McMurdo Dry Valleys, while limited, has great ecological significance to subnivian soil environments. Though sublimation dominates the ablation process in this region, measurable increases in soil moisture and insulation from temperature extremes provide more favorable conditions with respect to subnivian soil communities. While precipitation is not substantial, significant amounts of snow can accumulate, via wind transport, in topographic lees along the valley bottoms, forming thousands of discontinuous snow patches. These patches have the potential to act as significant sources of local meltwater, controlling biogeochemical cycling and the landscape distribution of microbial communities. Therefore, determining the spatial and temporal dynamics of snow at multiple scales is imperative to understanding the broader ecological role of snow in this region. High-resolution satellite imagery acquired during the 2009–2010 and 2010–2011 austral summers was used to quantify the distribution of snow across Taylor and Wright valleys. Extracted snow-covered area from the imagery was used as the basis for assessing inter-annual variability and seasonal controls on accumulation and ablation of snow at multiple scales. In addition to landscape analyses, fifteen 1 km 2 plots (3 in each of 5 study regions) were selected to assess the prevalence of snow cover at finer spatial scales, referred to herein as the snow-patch scale. Results confirm that snow patches tend to form in the same locations each year with some minor deviations observed. At the snow-patch scale, neighboring patches often exhibit considerable differences in aerial ablation rates, and particular snow patches do not reflect trends for snow-covered area observed at the landscape scale. These differences are presumably related to microtopographic influences acting on individual snow patches, such as wind sheltering and differences in snow depth due to the underlying topography. This highlights the importance of both the landscape and snow-patch scales in assessing the effects of snow cover on biogeochemical cycling and microbial communities.

Description: Ice-shelf buttressing and the stability of marine ice sheets The Cryosphere, 7, 647-655, 2013 Author(s): G. H. Gudmundsson Ice-shelf buttressing and the stability of marine-type ice sheets are investigated numerically. Buttressing effects are analysed for a situation where a stable grounding line is located on a bed sloping upwards in the direction of flow. Such grounding-line positions are known to be unconditionally unstable in the absence of transverse flow variations. It is shown that ice-shelf buttressing can restore stability under these conditions. Ice flux at the grounding line is, in general, not a monotonically increasing function of ice thickness. This, possibly at first somewhat counterintuitive result, is found to be fully consistent with recent theoretical work. Grounding lines on retrograde slopes are conditionally stable, and the stability regime is a non-trivial function of bed and ice-shelf geometry. The stability of grounding lines cannot be assessed from considerations of local bed slope only.

Description: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data The Cryosphere, 7, 615-630, 2013 Author(s): M. Tedesco, X. Fettweis, T. Mote, J. Wahr, P. Alexander, J. E. Box, and B. Wouters A combined analysis of remote sensing observations, regional climate model (RCM) outputs and reanalysis data over the Greenland ice sheet provides evidence that multiple records were set during summer 2012. Melt extent was the largest in the satellite era (extending up to ∼97% of the ice sheet) and melting lasted up to ∼2 months longer than the 1979–2011 mean. Model results indicate that near surface temperature was ∼3 standard deviations (σ) above the 1958–2011 mean, while surface mass balance (SMB) was ∼3σ below the mean and runoff was 3.9σ above the mean over the same period. Albedo, exposure of bare ice and surface mass balance also set new records, as did the total mass balance with summer and annual mass changes of, respectively, −627 Gt and −574 Gt, 2σ below the 2003–2012 mean. We identify persistent anticyclonic conditions over Greenland associated with anomalies in the North Atlantic Oscillation (NAO), changes in surface conditions (e.g., albedo, surface temperature) and preconditioning of surface properties from recent extreme melting as major driving mechanisms for the 2012 records. Less positive if not increasingly negative SMB will likely occur should these characteristics persist.

Description: Balanced conditions or slight mass gain of glaciers in the Lahaul and Spiti region (northern India, Himalaya) during the nineties preceded recent mass loss The Cryosphere, 7, 569-582, 2013 Author(s): C. Vincent, Al. Ramanathan, P. Wagnon, D. P. Dobhal, A. Linda, E. Berthier, P. Sharma, Y. Arnaud, M. F. Azam, P. G. Jose, and J. Gardelle The volume change of the Chhota Shigri Glacier (India, 32° 20 N, 77° 30' E) between 1988 and 2010 has been determined using in situ geodetic measurements. This glacier has experienced only a slight mass loss between 1988 and 2010 (−3.8 ± 2.0 m w.e. (water equivalent) corresponding to −0.17 ± 0.09 m w.e. yr −1 ). Using satellite digital elevation models (DEM) differencing and field measurements, we measure a negative mass balance (MB) between 1999 and 2010 (−4.8 ± 1.8 m w.e. corresponding to −0.44 ± 0.16 m w.e. yr −1 ). Thus, we deduce a slightly positive or near-zero MB between 1988 and 1999 (+1.0 ± 2.7 m w.e. corresponding to +0.09 ± 0.24 m w.e. yr −1 ). Furthermore, satellite DEM differencing reveals that the MB of the Chhota Shigri Glacier (−0.39 ± 0.15 m w.e. yr −1 ) has been only slightly less negative than the MB of a 2110 km 2 glaciarized area in the Lahaul and Spiti region (−0.44 ± 0.09 m w.e. yr −1 ) during 1999−2011. Hence, we conclude that the ice wastage is probably moderate in this region over the last 22 yr, with near equilibrium conditions during the nineties, and an ice mass loss after. The turning point from balanced to negative mass budget is not known but lies probably in the late nineties and at the latest in 1999. This positive or near-zero MB for Chhota Shigri Glacier (and probably for the surrounding glaciers of the Lahaul and Spiti region) during at least part of the 1990s contrasts with a recent compilation of MB data in the Himalayan range that indicated ice wastage since 1975. However, in agreement with this compilation, we confirm more negative balances since the beginning of the 21st century.

Description: The geomorphological effect of cornice fall avalanches in the Longyeardalen valley, Svalbard The Cryosphere, 7, 1361-1374, 2013 Author(s): M. Eckerstorfer, H. H. Christiansen, L. Rubensdotter, and S. Vogel The study of snow avalanches and their geomorphological effect in the periglacial parts of the cryosphere is important for enhanced geomorphological process understanding as well as hazard-related studies. Only a few field studies, and particularly few in the High Arctic, have quantified avalanche sedimentation. Snow avalanches are traditionally ranked behind rockfall in terms of their significance for mass-wasting processes of rockslopes. Cornice fall avalanches are at present the most dominant snow avalanche type at two slope systems, called Nybyen and Larsbreen, in the valley Longyeardalen in central Svalbard. Both slope systems are on northwest-facing lee slopes underneath a large summit plateau, with annual cornices forming on the top. High-frequency and magnitude cornice fall avalanching is observed by daily automatic time-lapse photography. In addition, rock debris sedimentation by cornice fall avalanches was measured directly in permanent sediment traps or by snow inventories. The results from a maximum of seven years of measurements in a total of 13 catchments show maximum mean rock debris sedimentation rates ranging from 8.2 to 38.7 kg m −2 at Nybyen, and from 0.8 to 55.4 kg m −2 at Larsbreen. Correspondingly, avalanche fan surfaces accreted from 2.6 to 8.8 mm yr −1 at Nybyen, and from 0.2 to 13.9 mm yr −1 at Larsbreen. This comparably efficient rockslope mass wasting is due to collapsing cornices producing cornice fall avalanches containing large amounts of rock debris throughout the entire winter. The rock debris of different origin stems from the plateau crests, the adjacent free rock face and the transport pathway, accumulating distinct avalanche fans at both slope systems. Cornice fall avalanche sedimentation also contributed to the development of a rock glacier at the Larsbreen site during the Holocene. We have recorded present maximum rockwall retreat rates of 0.9 mm yr −1 at Nybyen, but as much as 6.7 mm yr −1 at Larsbreen, while average Holocene rockwall retreat rates of 1.1 mm yr −1 at Nybyen have been determined earlier. As cornice fall avalanches are the dominant type of avalanche in central Svalbard, the related geomorphological effect is assumed to be of significance at periglacial landscape scale. A climate-induced shift in prevailing winter wind direction could change the rockslope sedimentation effectively by changing the snow avalanche activity.

Description: Heterogeneity in glacier response in the upper Shyok valley, northeast Karakoram The Cryosphere, 7, 1385-1398, 2013 Author(s): R. Bhambri, T. Bolch, P. Kawishwar, D. P. Dobhal, D. Srivastava, and B. Pratap Glaciers in the Karakoram show long-term irregular behaviour with comparatively frequent and sudden advances. A glacier inventory of the upper Shyok valley situated in northeast Karakoram has been generated for the year 2002 using Landsat ETM+ and SRTM3 DEM as baseline data for the investigations and subsequent change analysis. The upper Shyok valley contained 2123 glaciers (larger than 0.02 km 2 in size) with an area of 2977.9 ± 95.3 km 2 in 2002. Out of these, 18 glaciers with an area of 1004.1 ± 32.1 km 2 showed surge-type behaviour. Change analysis based on Hexagon KH-9 (years 1973 and 1974) and Landsat TM/ETM+ (years 1989, 2002 and 2011) images had to be restricted to a subset of 136 glaciers (covering an area of 1609.7 ± 51.5 km 2 in 2002) due to adverse snow conditions. The area of the investigated glaciers, including the 18 surge-type glaciers identified, showed no significant changes during all studied periods. However, the analysis provides a hint that the overall glacier area slightly decreased until about 1989 (area 1973: 1613.6 ± 43.6 km 2 ; area 1989: 1602.0 ± 33.6 km 2 ) followed by an increase (area 2002: 1609.7 ± 51.5; area 2011: 1615.8 ± 35.5 km 2 ). Although the overall change in area is insignificant, advances in glacier tongues since the end of the 1980s are clearly visible. Detailed estimations of length changes for individual glaciers since the 1970s and for Central Rimo Glacier since the 1930s confirm the irregular retreat and advance.

Description: Weak layer fracture: facets and depth hoar The Cryosphere, 7, 1447-1453, 2013 Author(s): I. Reiweger and J. Schweizer Understanding failure initiation within weak snow layers is essential for modeling and predicting dry-snow slab avalanches. We therefore performed laboratory experiments with snow samples containing a weak layer consisting of either faceted crystals or depth hoar. During these experiments the samples were loaded with different loading rates and at various tilt angles until fracture. The strength of the samples decreased with increasing loading rate and increasing tilt angle. Additionally, we took pictures of the side of four samples with a high-speed video camera and calculated the displacement using a particle image velocimetry (PIV) algorithm. The fracture process within the weak layer could thus be observed in detail. Catastrophic failure started due to a shear fracture just above the interface between the depth hoar layer and the underlying crust.

Description: Scatter of mass changes estimates at basin scale for Greenland and Antarctica The Cryosphere, 7, 1411-1432, 2013 Author(s): V. R. Barletta, L. S. Sørensen, and R. Forsberg During the last decade, the GRACE mission has provided valuable data for determining the mass changes of the Greenland and Antarctic ice sheets. Yet, discrepancies still exist in the published mass balance results, and comprehensive analyses on the sources of errors and discrepancies are lacking. Here, we present monthly mass changes together with trends derived from GRACE data at basin scale for both the Greenland and Antarctic ice sheets, and we assess the variability and errors for each of the possible sources of discrepancies, and we do this in an unprecedented systematic way, taking into account mass inference methods, data sets and background models. We find a very good agreement between the monthly mass change results derived from two independent methods, which represents a cross validation. For the monthly solutions, we find that most of the scatter is caused by the use of the two different data sets rather than the two different methods applied. Besides the well-known GIA trend uncertainty, we find that the geocenter motion and the recent de-aliasing corrections significantly impact the trends, with contributions of +13.2 Gt yr −1 and −20 Gt yr −1 , respectively, for Antarctica, which is more affected by these than Greenland. We show differences between the use of release RL04 and the new RL05 and confirm a lower noise content in the new release. The overall scatter of the solutions well exceeds the uncertainties propagated from the data errors and the leakage (as done in the past); hence we calculate new sound total errors for the monthly solutions and the trends. We find that the scatter in the monthly solutions caused by applying different estimates of geocenter motion time series (degree-1 corrections) is significant – contributing with up to 40% of the total error. For the whole GRACE period (2003–2011) our trend estimate for Greenland is −234 ± 20 Gt yr −1 and −83 ± 36 Gt yr −1 for Antarctica (−111 ± 15 Gt yr −1 in the western part). We also find a clear (with respect to our errors) increase of mass loss in the last four years.

Description: Are seasonal calving dynamics forced by buttressing from ice mélange or undercutting by melting? Outcomes from full-Stokes simulations of Store Glacier, West Greenland The Cryosphere, 8, 2353-2365, 2014 Author(s): J. Todd and P. Christoffersen We use a full-Stokes 2-D model (Elmer/Ice) to investigate the flow and calving dynamics of Store Glacier, a fast-flowing outlet glacier in West Greenland. Based on a new, subgrid-scale implementation of the crevasse depth calving criterion, we perform two sets of simulations: one to identify the primary forcing mechanisms and another to constrain future stability. We find that the mixture of icebergs and sea ice, known as ice mélange or sikussak, is principally responsible for the observed seasonal advance of the ice front. On the other hand, the effect of submarine melting on the calving rate of Store Glacier appears to be limited. Sensitivity analysis demonstrates that the glacier's calving dynamics are sensitive to seasonal perturbation, but are stable on interannual timescales due to the strong topographic control on the flow regime. Our results shed light on the dynamics of calving glaciers and may help explain why neighbouring glaciers do not necessarily respond synchronously to changes in atmospheric and oceanic forcing.

Description: Deglaciation of the Caucasus Mountains, Russia/Georgia, in the 21st century observed with ASTER satellite imagery and aerial photography The Cryosphere, 8, 2367-2379, 2014 Author(s): M. Shahgedanova, G. Nosenko, S. Kutuzov, O. Rototaeva, and T. Khromova Changes in the map area of 498 glaciers located on the Main Caucasus ridge (MCR) and on Mt. Elbrus in the Greater Caucasus Mountains (Russia and Georgia) were assessed using multispectral ASTER and panchromatic Landsat imagery with 15 m spatial resolution in 1999/2001 and 2010/2012. Changes in recession rates of glacier snouts between 1987–2001 and 2001–2010 were investigated using aerial photography and ASTER imagery for a sub-sample of 44 glaciers. In total, glacier area decreased by 4.7 ± 2.1% or 19.2 ± 8.7 km 2 from 407.3 ± 5.4 km 2 to 388.1 ± 5.2 km 2 . Glaciers located in the central and western MCR lost 13.4 ± 7.3 km 2 (4.7 ± 2.5%) in total or 8.5 km 2 (5.0 ± 2.4%) and 4.9 km 2 (4.1 ± 2.7%) respectively. Glaciers on Mt. Elbrus, although located at higher elevations, lost 5.8 ± 1.4 km 2 (4.9 ± 1.2%) of their total area. The recession rates of valley glacier termini increased between 1987–2000/01 and 2000/01–2010 (2000 for the western MCR and 2001 for the central MCR and Mt.~Elbrus) from 3.8 ± 0.8, 3.2 ± 0.9 and 8.3 ± 0.8 m yr −1 to 11.9 ± 1.1, 8.7 ± 1.1 and 14.1 ± 1.1 m yr −1 in the central and western MCR and on Mt. Elbrus respectively. The highest rate of increase in glacier termini retreat was registered on the southern slope of the central MCR where it has tripled. A positive trend in summer temperatures forced glacier recession, and strong positive temperature anomalies in 1998, 2006, and 2010 contributed to the enhanced loss of ice. An increase in accumulation season precipitation observed in the northern MCR since the mid-1980s has not compensated for the effects of summer warming while the negative precipitation anomalies, observed on the southern slope of the central MCR in the 1990s, resulted in stronger glacier wastage.

Description: Glacier-like forms on Mars The Cryosphere, 8, 2047-2061, 2014 Author(s): B. Hubbard, C. Souness, and S. Brough More than 1300 glacier-like forms (GLFs) are located in Mars' mid-latitudes. These GLFs are predominantly composed of ice–dust mixtures and are visually similar to terrestrial valley glaciers, showing signs of downhill viscous deformation and an expanded former extent. However, several fundamental aspects of their behavior are virtually unknown, including temporal and spatial variations in mass balance, ice motion, landscape erosion and deposition, and hydrology. Here, we investigate the physical glaciology of martian GLFs. We use satellite images of specific examples and case studies to build on existing knowledge relating to (i) GLF current and former extent, exemplified via a GLF located in Phlegra Montes; (ii) indicators of GLF motion, focusing on the presence of surface crevasses on several GLFs; (iii) processes of GLF debris transfer, focusing on mapping and interpreting boulder trains on one GLF located in Protonilus Mensae, the analysis of which suggests a best-estimate mean GLF flow speed of 7.5 mm a −1 ; and (iv) GLF hydrology, focusing on supra-GLF gulley networks. On the basis of this information, we summarize the current state of knowledge of the glaciology of martian GLFs and identify future research avenues.

Description: Modeling the elastic transmission of tidal stresses to great distances inland in channelized ice streams The Cryosphere, 8, 2007-2029, 2014 Author(s): J. Thompson, M. Simons, and V. C. Tsai Geodetic surveys suggest that ocean tides can modulate the motion of Antarctic ice streams, even at stations many tens of kilometers inland from the grounding line. These surveys suggest that ocean tidal stresses can perturb ice stream motion at distances about an order of magnitude farther inland than tidal flexure of the ice stream alone. Recent models exploring the role of tidal perturbations in basal shear stress are primarily one- or two-dimensional, with the impact of the ice stream margins either ignored or parameterized. Here, we use two- and three-dimensional finite-element modeling to investigate transmission of tidal stresses in ice streams and the impact of considering more realistic, three-dimensional ice stream geometries. Using Rutford Ice Stream as a real-world comparison, we demonstrate that the assumption that elastic tidal stresses in ice streams propagate large distances inland fails for channelized glaciers due to an intrinsic, exponential decay in the stress caused by resistance at the ice stream margins. This behavior is independent of basal conditions beneath the ice stream and cannot be fit to observations using either elastic or nonlinear viscoelastic rheologies without nearly complete decoupling of the ice stream from its lateral margins. Our results suggest that a mechanism external to the ice stream is necessary to explain the tidal modulation of stresses far upstream of the grounding line for narrow ice streams. We propose a hydrologic model based on time-dependent variability in till strength to explain transmission of tidal stresses inland of the grounding line. This conceptual model can reproduce observations from Rutford Ice Stream.

Description: Fluctuations of a Greenlandic tidewater glacier driven by changes in atmospheric forcing: observations and modelling of Kangiata Nunaata Sermia, 1859–present The Cryosphere, 8, 2031-2045, 2014 Author(s): J. M. Lea, D. W. F. Mair, F. M. Nick, B. R. Rea, D. van As, M. Morlighem, P. W. Nienow, and A. Weidick Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Description: The surface energy balance of a polygonal tundra site in northern Siberia – Part 2: Winter The Cryosphere, 5, 509-524, 2011 Author(s): M. Langer, S. Westermann, S. Muster, K. Piel, and J. Boike In this study, we present the winter time surface energy balance at a polygonal tundra site in northern Siberia based on independent measurements of the net radiation, the sensible heat flux and the ground heat flux from two winter seasons. The latent heat flux is inferred from measurements of the atmospheric turbulence characteristics and a model approach. The long-wave radiation is found to be the dominant factor in the surface energy balance. The radiative losses are balanced to about 60 % by the ground heat flux and almost 40 % by the sensible heat fluxes, whereas the contribution of the latent heat flux is small. The main controlling factors of the surface energy budget are the snow cover, the cloudiness and the soil temperature gradient. Large spatial differences in the surface energy balance are observed between tundra soils and a small pond. The ground heat flux released at a freezing pond is by a factor of two higher compared to the freezing soil, whereas large differences in net radiation between the pond and soil are only observed at the end of the winter period. Differences in the surface energy balance between the two winter seasons are found to be related to differences in snow depth and cloud cover which strongly affect the temperature evolution and the freeze-up at the investigated pond.

Description: A full Stokes ice flow model for the vicinity of Dome Fuji, Antarctica, with induced anisotropy and fabric evolution The Cryosphere, 5, 495-508, 2011 Author(s): H. Seddik, R. Greve, T. Zwinger, and L. Placidi A three-dimensional, thermo-mechanically coupled ice flow model with induced anisotropy has been applied to a ~200 × 200 km domain around the Dome Fuji drill site, Antarctica. The model ("Elmer/Ice") is based on the open-source multi-physics package Elmer ( http://www.csc.fi/elmer/ ) and solves the full Stokes equations. Flow-induced anisotropy in ice is accounted for by an implementation of the C ontinuum-mechanical, A nisotropic F low model, based on an anisotropic F low E nhancement factor ("CAFFE model"). Steady-state simulations for present-day climate conditions are conducted. The main findings are: (i) the flow regime at Dome Fuji is a complex superposition of vertical compression, horizontal extension and bed-parallel shear; (ii) for an assumed geothermal heat flux of 60 mW m −2 the basal temperature at Dome Fuji reaches the pressure melting point and the basal melting rate is ~0.35 mm a −1 ; (iii) in agreement with observational data, the fabric shows a strong single maximum at Dome Fuji, and the age of the ice is decreased compared to an isotropic scenario; (iv) as a consequence of spatially variable basal melting conditions, the basal age tends to be smaller where the ice is thicker and larger where the ice is thinner. The latter result is of great relevance for the consideration of a future drill site in the area.

Description: Scale-dependent measurement and analysis of ground surface temperature variability in alpine terrain The Cryosphere, 5, 431-443, 2011 Author(s): S. Gubler, J. Fiddes, M. Keller, and S. Gruber Measurements of environmental variables are often used to validate and calibrate physically-based models. Depending on their application, the models are used at different scales, ranging from few meters to tens of kilometers. Environmental variables can vary strongly within the grid cells of these models. Validating a model with a single measurement is therefore delicate and susceptible to induce bias in further model applications. To address the question of uncertainty associated with scale in permafrost models, we present data of 390 spatially-distributed ground surface temperature measurements recorded in terrain of high topographic variability in the Swiss Alps. We illustrate a way to program, deploy and refind a large number of measurement devices efficiently, and present a strategy to reduce data loss reported in earlier studies. Data after the first year of deployment is presented. The measurements represent the variability of ground surface temperatures at two different scales ranging from few meters to some kilometers. On the coarser scale, the dependence of mean annual ground surface temperature on elevation, slope, aspect and ground cover type is modelled with a multiple linear regression model. Sampled mean annual ground surface temperatures vary from −4 °C to 5 °C within an area of approximately 16 km 2 subject to elevational differences of approximately 1000 m. The measurements also indicate that mean annual ground surface temperatures vary up to 6 °C (i.e., from −2 °C to 4 °C) even within an elevational band of 300 m. Furthermore, fine-scale variations can be high (up to 2.5 °C) at distances of less than 14 m in homogeneous terrain. The effect of this high variability of an environmental variable on model validation and applications in alpine regions is discussed.

Description: A new approach to mapping permafrost and change incorporating uncertainties in ground conditions and climate projections The Cryosphere, 8, 2177-2194, 2014 Author(s): Y. Zhang, I. Olthof, R. Fraser, and S. A. Wolfe Spatially detailed information on permafrost distribution and change with climate is important for land use planning, infrastructure development, and environmental assessments. However, the required soil and surficial geology maps in the North are coarse, and projected climate scenarios vary widely. Considering these uncertainties, we propose a new approach to mapping permafrost distribution and change by integrating remote sensing data, field measurements, and a process-based model. Land cover types from satellite imagery are used to capture the general land conditions and to improve the resolution of existing permafrost maps. For each land cover type, field observations are used to estimate the probabilities of different ground conditions. A process-based model is used to quantify the evolution of permafrost for each ground condition under three representative climate scenarios (low, medium, and high warming). From the model results, the probability of permafrost occurrence and the most likely permafrost conditions are determined. We apply this approach at 20 m resolution to a large area in Northwest Territories, Canada. Mapped permafrost conditions are in agreement with field observations and other studies. The data requirements, model robustness, and computation time are reasonable, and this approach may serve as a practical means to mapping permafrost and changes at high resolution in other regions.

Description: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice The Cryosphere, 8, 2219-2233, 2014 Author(s): S. Arndt and M. Nicolaus Arctic sea ice has not only decreased in volume during the last decades, but has also changed in its physical properties towards a thinner and more seasonal ice cover. These changes strongly impact the energy budget, and might affect the ice-associated ecosystems. In this study, we quantify solar shortwave fluxes through sea ice for the entire Arctic during all seasons. To focus on sea-ice-related processes, we exclude fluxes through open water, scaling linearly with sea ice concentration. We present a new parameterization of light transmittance through sea ice for all seasons as a function of variable sea ice properties. The maximum monthly mean solar heat flux under the ice of 30 × 10 5 Jm −2 occurs in June, enough heat to melt 0.3 m of sea ice. Furthermore, our results suggest that 96% of the annual solar heat input through sea ice occurs during only a 4-month period from May to August. Applying the new parameterization to remote sensing and reanalysis data from 1979 to 2011, we find an increase in transmitted light of 1.5% yr −1 for all regions. This corresponds to an increase in potential sea ice bottom melt of 63% over the 33-year study period. Sensitivity studies reveal that the results depend strongly on the timing of melt onset and the correct classification of ice types. Assuming 2 weeks earlier melt onset, the annual transmitted solar radiation to the upper ocean increases by 20%. Continuing the observed transition from a mixed multi-year/first-year sea ice cover to a seasonal ice cover results in an increase in light transmittance by an additional 18%.

Description: Post-LIA glacier changes along a latitudinal transect in the Central Italian Alps The Cryosphere, 8, 2235-2252, 2014 Author(s): R. Scotti, F. Brardinoni, and G. B. Crosta The variability of glacier response to atmospheric temperature rise in different topo-climatic settings is still a matter of debate. To address this question in the Central Italian Alps, we compile a post-LIA (Little Ice Age) multitemporal glacier inventory (1860–1954–1990–2003–2007) along a latitudinal transect that originates north of the continental divide in the Livigno Mountains and extends south through the Disgrazia and Orobie ranges, encompassing continental-to-maritime climatic settings. In these sub-regions, we examine the area change of 111 glaciers. Overall, the total glacierized area has declined from 34.1 to 10.1 km 2 , with a substantial increase in the number of small glaciers due to fragmentation. The average annual decrease (AAD) in glacier area has risen by about 1 order of magnitude from 1860–1990 (Livigno: 0.45; Orobie: 0.42; and Disgrazia: 0.39 % a −1 ) to 1990–2007 (Livigno: 3.08; Orobie: 2.44; and Disgrazia: 2.27 % a −1 ). This ranking changes when considering glaciers smaller than 0.5 km 2 only (i.e., we remove the confounding caused by large glaciers in Disgrazia), so that post-1990 AAD follows the latitudinal gradient and Orobie glaciers stand out (Livigno: 4.07; Disgrazia: 3.57; and Orobie: 2.47 % a −1 ). More recent (2007–2013) field-based mass balances in three selected small glaciers confirm post-1990 trends showing the consistently highest retreat in continental Livigno and minimal area loss in maritime Orobie, with Disgrazia displaying transitional behavior. We argue that the recent resilience of glaciers in Orobie is a consequence of their decoupling from synoptic atmospheric temperature trends, a decoupling that arises from the combination of local topographic configuration (i.e., deep, north-facing cirques) and high winter precipitation, which ensures high snow-avalanche supply, as well as high summer shading and sheltering. Our hypothesis is further supported by the lack of correlations between glacier change and glacier attributes in Orobie, as well as by the higher variability in ELA,sub〉0 positioning, post-LIA glacier change, and interannual mass balances, as we move southward along the transect.

Description: Melting trends over the Greenland ice sheet (1958–2009) from spaceborne microwave data and regional climate models The Cryosphere, 5, 359-375, 2011 Author(s): X. Fettweis, M. Tedesco, M. van den Broeke, and J. Ettema To study near-surface melt changes over the Greenland ice sheet (GrIS) since 1979, melt extent estimates from two regional climate models were compared with those obtained from spaceborne microwave brightness temperatures using two different remote sensing algorithms. The results from the two models were consistent with those obtained with the remote sensing algorithms at both daily and yearly time scales, encouraging the use of the models for analyzing melting trends before the satellite era (1958–1979), when forcing data is available. Differences between satellite-derived and model-simulated results still occur and are used here to identify (i) biases in the snow models (notably in the albedo parametrization, in the thickness of a snow layer, in the maximum liquid water content within the snowpack and in the snowfall impacting the bare ice appearance in summer) and (ii) limitations in the use of passive microwave data for snowmelt detection at the edge of the ice sheet due to mixed pixel effect (e.g., tundra or rock nearby the ice sheet). The results from models and spaceborne microwave sensors confirm a significant (p-value = 0.01) increase in GrIS surface melting since 1979. The melt extent recorded over the last years (1998, 2003, 2005 and 2007) is unprecedented in the last 50 yr with the cumulated melt area in the 2000's being, on the average, twice that of the 1980's.

Description: Isotope hydrological studies of the perennial ice deposit of Saarhalle, Mammuthöhle, Dachstein Mts, Austria The Cryosphere, 5, 291-298, 2011 Author(s): Z. Kern, I. Fórizs, R. Pavuza, M. Molnár, and B. Nagy A 5.28 m-long ice core was extracted from a major cave ice body in the Mammuthöhle cave system. The upper ~1.2 m of ice most likely originate from precipitation fallen before the 1960s (based on

Description: Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change The Cryosphere, 5, 271-290, 2011 Author(s): C. Nuth and A. Kääb There are an increasing number of digital elevation models (DEMs) available worldwide for deriving elevation differences over time, including vertical changes on glaciers. Most of these DEMs are heavily post-processed or merged, so that physical error modelling becomes difficult and statistical error modelling is required instead. We propose a three-step methodological framework for assessing and correcting DEMs to quantify glacier elevation changes: (i) remove DEM shifts, (ii) check for elevation-dependent biases, and (iii) check for higher-order, sensor-specific biases. A simple, analytic and robust method to co-register elevation data is presented in regions where stable terrain is either plentiful (case study New Zealand) or limited (case study Svalbard). The method is demonstrated using the three global elevation data sets available to date, SRTM, ICESat and the ASTER GDEM, and with automatically generated DEMs from satellite stereo instruments of ASTER and SPOT5-HRS. After 3-D co-registration, significant biases related to elevation were found in some of the stereoscopic DEMs. Biases related to the satellite acquisition geometry (along/cross track) were detected at two frequencies in the automatically generated ASTER DEMs. The higher frequency bias seems to be related to satellite jitter , most apparent in the back-looking pass of the satellite. The origins of the more significant lower frequency bias is uncertain. ICESat-derived elevations are found to be the most consistent globally available elevation data set available so far. Before performing regional-scale glacier elevation change studies or mosaicking DEMs from multiple individual tiles (e.g. ASTER GDEM), we recommend to co-register all elevation data to ICESat as a global vertical reference system.

Description: Present dynamics and future prognosis of a slowly surging glacier The Cryosphere, 5, 299-313, 2011 Author(s): G. E. Flowers, N. Roux, S. Pimentel, and C. G. Schoof Glacier surges are a well-known example of an internal dynamic oscillation whose occurrence is not a direct response to the external climate forcing, but whose character (i.e. period, amplitude, mechanism) may depend on the glacier's environmental or climate setting. We examine the dynamics of a small (∼5 km 2 ) valley glacier in Yukon, Canada, where two previous surges have been photographically documented and an unusually slow surge is currently underway. To characterize the dynamics of the present surge, and to speculate on the future of this glacier, we employ a higher-order flowband model of ice dynamics with a regularized Coulomb-friction sliding law in both diagnostic and prognostic simulations. Diagnostic (force balance) calculations capture the measured ice-surface velocity profile only when non-zero basal water pressures are prescribed over the central region of the glacier, coincident with where evidence of the surge has been identified. This leads to sliding accounting for 50–100% of the total surface motion in this region. Prognostic simulations, where the glacier geometry evolves in response to a prescribed surface mass balance, reveal a significant role played by a bedrock ridge beneath the current equilibrium line of the glacier. Ice thickening occurs above the ridge in our simulations, until the net mass balance reaches sufficiently negative values. We suggest that the bedrock ridge may contribute to the propensity for surges in this glacier by promoting the development of the reservoir area during quiescence, and may permit surges to occur under more negative balance conditions than would otherwise be possible. Collectively, these results corroborate our interpretation of the current glacier flow regime as indicative of a slow surge that has been ongoing for some time, and support a relationship between surge incidence or character and the net mass balance. Our results also highlight the importance of glacier bed topography in controlling ice dynamics, as observed in many other glacier systems.

Description: In-situ multispectral and bathymetric measurements over a supraglacial lake in western Greenland using a remotely controlled watercraft The Cryosphere, 5, 445-452, 2011 Author(s): M. Tedesco and N. Steiner Supraglacial lakes form from meltwater on the Greenland ice sheet in topographic depressions on the surface, affecting both surface and sub-glacial processes. As the reflectance in the visible and near-infrared regions of a column of water is modulated by its height, retrieval techniques using spaceborne remote sensing data (e.g. Landsat, MODIS) have been proposed in the literature for the detection of lakes and estimation of their volume. These techniques require basic assumptions on the spectral properties of the water as well as the bottom of the lake, among other things. In this study, we report results obtained from the analysis of concurrent in-situ multi-spectral and depth measurements collected over a supraglacial lake during early July 2010 in West Greenland (Lake Olivia, 69°36'35" N, 49°29'40" W) and aim to assess some of the underlying hypotheses in remote sensing based bathymetric approaches. In particular, we focus our attention on the analysis of the lake bottom albedo and of the water attenuation coefficient. The analysis of in-situ data (collected by means of a remotely controlled boat equipped with a GPS, a sonar and a spectrometer) highlights the exponential trend of the water-leaving reflectance with lake depth. The values of the attenuation factor obtained from in-situ data are compared with those computed using approaches proposed in the literature. Also, the values of the lake bottom albedo from in-situ measurements are compared with those obtained from the analysis of reflectance of shallow waters. Finally, we quantify the error between in-situ measured and satellite-estimated lake depth values for the lake under study.

Description: Temperature variability and offset in steep alpine rock and ice faces The Cryosphere, 5, 977-988, 2011 Author(s): A. Hasler, S. Gruber, and W. Haeberli The thermal condition of high-alpine mountain flanks can be an important determinant of climate change impact on slope stability and correspondingly down-slope hazard regimes. In this study we analyze time-series from 17 shallow temperature-depth profiles at two field sites in steep bedrock and ice. Extending earlier studies that revealed the topographic variations in temperatures, we demonstrate considerable differences of annual mean temperatures for variable surface characteristics and depths within the measured profiles. This implies that measurements and model related to compact and near-vertical bedrock temperatures may deviate considerably from conditions in the majority of bedrock slopes in mountain ranges that are usually non-vertical and fractured. For radiation-exposed faces mean annual temperatures at depth are up to 3 °C lower and permafrost is likely to exist at lower elevations than reflected by estimates based on near-vertical homogeneous cases. Retention of a thin snow cover and ventilation effects in open clefts are most likely responsible for this cooling. The measurements presented or similar data could be used in the future to support the development and testing of models related to the thermal effect of snow-cover and fractures in steep bedrock.

Description: On the influence of model physics on simulations of Arctic and Antarctic sea ice The Cryosphere, 5, 687-699, 2011 Author(s): F. Massonnet, T. Fichefet, H. Goosse, M. Vancoppenolle, P. Mathiot, and C. König Beatty Two hindcast (1983–2007) simulations are performed with the global, ocean-sea ice models NEMO-LIM2 and NEMO-LIM3 driven by atmospheric reanalyses and climatologies. The two simulations differ only in their sea ice component, while all other elements of experimental design (resolution, initial conditions, atmospheric forcing) are kept identical. The main differences in the sea ice models lie in the formulation of the subgrid-scale ice thickness distribution, of the thermodynamic processes, of the sea ice salinity and of the sea ice rheology. To assess the differences in model skill over the period of investigation, we develop a set of metrics for both hemispheres, comparing the main sea ice variables (concentration, thickness and drift) to available observations and focusing on both mean state and seasonal to interannual variability. Based upon these metrics, we discuss the physical processes potentially responsible for the differences in model skill. In particular, we suggest that (i) a detailed representation of the ice thickness distribution increases the seasonal to interannual variability of ice extent, with spectacular improvement for the simulation of the recent observed summer Arctic sea ice retreats, (ii) the elastic-viscous-plastic rheology enhances the response of ice to wind stress, compared to the classical viscous-plastic approach, (iii) the grid formulation and the air-sea ice drag coefficient affect the simulated ice export through Fram Strait and the ice accumulation along the Canadian Archipelago, and (iv) both models show less skill in the Southern Ocean, probably due to the low quality of the reanalyses in this region and to the absence of important small-scale oceanic processes at the models' resolution (~1°).

Description: Glacier contribution to streamflow in two headwaters of the Huasco River, Dry Andes of Chile The Cryosphere, 5, 1099-1113, 2011 Author(s): S. Gascoin, C. Kinnard, R. Ponce, S. Lhermitte, S. MacDonell, and A. Rabatel Quantitative assessment of glacier contribution to present-day streamflow is a prerequisite to the anticipation of climate change impact on water resources in the Dry Andes. In this paper we focus on two glaciated headwater catchments of the Huasco Basin (Chile, 29° S). The combination of glacier monitoring data for five glaciers (Toro 1, Toro 2, Esperanza, Guanaco, Estrecho and Ortigas) with five automatic streamflow records at sites with glacier coverage of 0.4 to 11 % allows the estimation of the mean annual glacier contribution to discharge between 2003/2004 and 2007/2008 hydrological years. In addition, direct manual measurements of glacier runoff were conducted in summer at the snouts of four glaciers, which provide the instantaneous contribution of glacier meltwater to stream runoff during summer. The results show that the mean annual glacier contribution to streamflow ranges between 3.3 and 23 %, which is greater than the glaciated fraction of the catchments. We argue that glacier contribution is partly enhanced by the effect of snowdrift from the non-glacier area to the glacier surface. Glacier mass loss is evident over the study period, with a mean of −0.84 m w.e. yr −1 for the period 2003/2004–2007/2008, and also contributes to increase glacier runoff. An El Niño episode in 2002 resulted in high snow accumulation, modifying the hydrological regime and probably reducing the glacier contribution in favor of seasonal snowmelt during the subsequent 2002/2003 hydrological year. At the hourly timescale, summertime glacier contributions are highly variable in space and time, revealing large differences in effective melting rates between glaciers and glacierets (from 1 mm w.e. h −1 to 6 mm w.e. h −1 ).

Description: Coupling of climate models and ice sheet models by surface mass balance gradients: application to the Greenland Ice Sheet The Cryosphere, 6, 255-272, 2012 Author(s): M. M. Helsen, R. S. W. van de Wal, M. R. van den Broeke, W. J. van de Berg, and J. Oerlemans It is notoriously difficult to couple surface mass balance (SMB) results from climate models to the changing geometry of an ice sheet model. This problem is traditionally avoided by using only accumulation from a climate model, and parameterizing the meltwater run-off as a function of temperature, which is often related to surface elevation ( H s ). In this study, we propose a new strategy to calculate SMB, to allow a direct adjustment of SMB to a change in ice sheet topography and/or a change in climate forcing. This method is based on elevational gradients in the SMB field as computed by a regional climate model. Separate linear relations are derived for ablation and accumulation, using pairs of H s and SMB within a minimum search radius. The continuously adjusting SMB forcing is consistent with climate model forcing fields, also for initially non-glaciated areas in the peripheral areas of an ice sheet. When applied to an asynchronous coupled ice sheet – climate model setup, this method circumvents traditional temperature lapse rate assumptions. Here we apply it to the Greenland Ice Sheet (GrIS). Experiments using both steady-state forcing and glacial-interglacial forcing result in realistic ice sheet reconstructions.

Description: Estimating ice phenology on large northern lakes from AMSR-E: algorithm development and application to Great Bear Lake and Great Slave Lake, Canada The Cryosphere, 6, 235-254, 2012 Author(s): K.-K. Kang, C. R. Duguay, and S. E. L. Howell Time series of brightness temperatures ( T B ) from the Advanced Microwave Scanning Radiometer–Earth Observing System (AMSR-E) are examined to determine ice phenology variables on the two largest lakes of northern Canada: Great Bear Lake (GBL) and Great Slave Lake (GSL). T B measurements from the 18.7, 23.8, 36.5, and 89.0 GHz channels (H- and V- polarization) are compared to assess their potential for detecting freeze-onset/melt-onset and ice-on/ice-off dates on both lakes. The 18.7 GHz (H-pol) channel is found to be the most suitable for estimating these ice dates as well as the duration of the ice cover and ice-free seasons. A new algorithm is proposed using this channel and applied to map all ice phenology variables on GBL and GSL over seven ice seasons (2002–2009). Analysis of the spatio-temporal patterns of each variable at the pixel level reveals that: (1) both freeze-onset and ice-on dates occur on average about one week earlier on GBL than on GSL (Day of Year (DY) 318 and 333 for GBL; DY 328 and 343 for GSL); (2) the freeze-up process or freeze duration (freeze-onset to ice-on) takes a slightly longer amount of time on GBL than on GSL (about 1 week on average); (3) melt-onset and ice-off dates occur on average one week and approximately four weeks later, respectively, on GBL (DY 143 and 183 for GBL; DY 135 and 157 for GSL); (4) the break-up process or melt duration (melt-onset to ice-off) lasts on average about three weeks longer on GBL; and (5) ice cover duration estimated from each individual pixel is on average about three weeks longer on GBL compared to its more southern counterpart, GSL. A comparison of dates for several ice phenology variables derived from other satellite remote sensing products (e.g. NOAA Interactive Multisensor Snow and Ice Mapping System (IMS), QuikSCAT, and Canadian Ice Service Database) show that, despite its relatively coarse spatial resolution, AMSR-E 18.7 GHz provides a viable means for monitoring of ice phenology on large northern lakes.

Description: Large surface meltwater discharge from the Kangerlussuaq sector of the Greenland ice sheet during the record-warm year 2010 explained by detailed energy balance observations The Cryosphere, 6, 199-209, 2012 Author(s): D. van As, A. L. Hubbard, B. Hasholt, A. B. Mikkelsen, M. R. van den Broeke, and R. S. Fausto This study uses data from six on-ice weather stations, calibrated MODIS-derived albedo and proglacial river gauging measurements to drive and validate an energy balance model. We aim to quantify the record-setting positive temperature anomaly in 2010 and its effect on mass balance and runoff from the Kangerlussuaq sector of the Greenland ice sheet. In 2010, the average temperature was 4.9 °C (2.7 standard deviations) above the 1974–2010 average in Kangerlussuaq. High temperatures were also observed over the ice sheet, with the magnitude of the positive anomaly increasing with altitude, particularly in August. Simultaneously, surface albedo was anomalously low in 2010, predominantly in the upper ablation zone. The low albedo was caused by high ablation, which in turn profited from high temperatures and low winter snowfall. Surface energy balance calculations show that the largest melt excess (∼170%) occurred in the upper ablation zone (above 1000 m), where higher temperatures and lower albedo contributed equally to the melt anomaly. At lower elevations the melt excess can be attributed to high atmospheric temperatures alone. In total, we calculate that 6.6 ± 1.0 km 3 of surface meltwater ran off the ice sheet in the Kangerlussuaq catchment in 2010, exceeding the reference year 2009 (based on atmospheric temperature measurements) by ∼150%. During future warm episodes we can expect a melt response of at least the same magnitude, unless a larger wintertime snow accumulation delays and moderates the melt-albedo feedback. Due to the hypsometry of the ice sheet, yielding an increasing surface area with elevation, meltwater runoff will be further amplified by increases in melt forcings such as atmospheric heat.

Description: Climatic drivers of seasonal glacier mass balances: an analysis of 6 decades at Glacier de Sarennes (French Alps) The Cryosphere, 7, 47-66, 2013 Author(s): E. Thibert, N. Eckert, and C. Vincent Refined temporal signals extracted from a winter and summer mass balance series recorded at Glacier de Sarennes (French Alps) using variance decomposition are related to local meteorological data and large-scale North Atlantic Oscillation (NAO) anomalies in terms of interannual variability, trends of the low-frequency signals, and breaks in the time series. The winter balance has increased by +23% since 1976 due to more precipitation in early and late winter. The summer balance has decreased since 1982 due to a 43% increase in snow and ice melt. A 24-day lengthening of the ablation period – mainly due to longer ice ablation – is the main component in the overall increase in ablation. In addition, the last 25 yr have seen increases in ablation rates of 14 and 10% for snow and ice, respectively. A simple degree-day analysis can account for both the snow/ice melt rate rise and the lengthening of the ablation period as a function of higher air temperatures. From the same analysis, the equilibrium-line altitude of this 45° N latitude south-facing glacier has a sensitivity to temperature of +93 m °C −1 around its mean elevation of 3100 m a.s.l. over 6 decades. The sensitivity of summer balance to temperature is −0.62 m w.e. yr −1 °C −1 for a typical 125-day-long ablation season. Finally, the correlation of winter and summer mass balance terms with NAO anomalies is investigated. Singularly, highest values are obtained between winter NAO anomalies and summer balance. Winter NAO anomalies and winter balance and precipitation are almost disconnected. However, these results strongly depend on how the NAO signal is smoothed, so that the link between Sarennes mass balance seasonal terms and NAO signal remains tenuous and hard to interpret.

Description: A synthesis of the Antarctic surface mass balance during the last 800 yr The Cryosphere, 7, 303-319, 2013 Author(s): M. Frezzotti, C. Scarchilli, S. Becagli, M. Proposito, and S. Urbini Global climate models suggest that Antarctic snowfall should increase in a warming climate and mitigate rises in the sea level. Several processes affect surface mass balance (SMB), introducing large uncertainties in past, present and future ice sheet mass balance. To provide an extended perspective on the past SMB of Antarctica, we used 67 firn/ice core records to reconstruct the temporal variability in the SMB over the past 800 yr and, in greater detail, over the last 200 yr. Our SMB reconstructions indicate that the SMB changes over most of Antarctica are statistically negligible and that the current SMB is not exceptionally high compared to the last 800 yr. High-accumulation periods have occurred in the past, specifically during the 1370s and 1610s. However, a clear increase in accumulation of more than 10% has occurred in high SMB coastal regions and over the highest part of the East Antarctic ice divide since the 1960s. To explain the differences in behaviour between the coastal/ice divide sites and the rest of Antarctica, we suggest that a higher frequency of blocking anticyclones increases the precipitation at coastal sites, leading to the advection of moist air in the highest areas, whereas blowing snow and/or erosion have significant negative impacts on the SMB at windy sites. Eight hundred years of stacked records of the SMB mimic the total solar irradiance during the 13th and 18th centuries. The link between those two variables is probably indirect and linked to a teleconnection in atmospheric circulation that forces complex feedback between the tropical Pacific and Antarctica via the generation and propagation of a large-scale atmospheric wave train.

Description: Future Arctic marine access: analysis and evaluation of observations, models, and projections of sea ice The Cryosphere, 7, 321-332, 2013 Author(s): T. S. Rogers, J. E. Walsh, T. S. Rupp, L. W. Brigham, and M. Sfraga There is an emerging need for regional applications of sea ice projections to provide more accuracy and greater detail to scientists, national, state and local planners, and other stakeholders. The present study offers a prototype for a comprehensive, interdisciplinary study to bridge observational data, climate model simulations, and user needs. The study's first component is an observationally based evaluation of Arctic sea ice trends during 1980–2008, with an emphasis on seasonal and regional differences relative to the overall pan-Arctic trend. Regional sea ice loss has varied, with a significantly larger decline of winter maximum (January–March) extent in the Atlantic region than in other sectors. A lead–lag regression analysis of Atlantic sea ice extent and ocean temperatures indicates that reduced sea ice extent is associated with increased Atlantic Ocean temperatures. Correlations between the two variables are greater when ocean temperatures lag rather than lead sea ice. The performance of 13 global climate models is evaluated using three metrics to compare sea ice simulations with the observed record. We rank models over the pan-Arctic domain and regional quadrants and synthesize model performance across several different studies. The best performing models project reduced ice cover across key access routes in the Arctic through 2100, with a lengthening of seasons for marine operations by 1–3 months. This assessment suggests that the Northwest and Northeast Passages hold potential for enhanced marine access to the Arctic in the future, including shipping and resource development opportunities.

Description: Greenland ice sheet surface mass balance: evaluating simulations and making projections with regional climate models The Cryosphere, 6, 1275-1294, 2012 Author(s): J. G. L. Rae, G. Aðalgeirsdóttir, T. L. Edwards, X. Fettweis, J. M. Gregory, H. T. Hewitt, J. A. Lowe, P. Lucas-Picher, R. H. Mottram, A. J. Payne, J. K. Ridley, S. R. Shannon, W. J. van de Berg, R. S. W. van de Wal, and M. R. van den Broeke Four high-resolution regional climate models (RCMs) have been set up for the area of Greenland, with the aim of providing future projections of Greenland ice sheet surface mass balance (SMB), and its contribution to sea level rise, with greater accuracy than is possible from coarser-resolution general circulation models (GCMs). This is the first time an intercomparison has been carried out of RCM results for Greenland climate and SMB. Output from RCM simulations for the recent past with the four RCMs is evaluated against available observations. The evaluation highlights the importance of using a detailed snow physics scheme, especially regarding the representations of albedo and meltwater refreezing. Simulations with three of the RCMs for the 21st century using SRES scenario A1B from two GCMs produce trends of between −5.5 and −1.1 Gt yr −2 in SMB (equivalent to +0.015 and +0.003 mm sea level equivalent yr −2 ), with trends of smaller magnitude for scenario E1, in which emissions are mitigated. Results from one of the RCMs whose present-day simulation is most realistic indicate that an annual mean near-surface air temperature increase over Greenland of ~ 2°C would be required for the mass loss to increase such that it exceeds accumulation, thereby causing the SMB to become negative, which has been suggested as a threshold beyond which the ice sheet would eventually be eliminated.

Description: Simulating snow maps for Norway: description and statistical evaluation of the seNorge snow model The Cryosphere, 6, 1323-1337, 2012 Author(s): T. M. Saloranta Daily maps of snow conditions have been produced in Norway with the seNorge snow model since 2004. The seNorge snow model operates with 1 × 1 km resolution, uses gridded observations of daily temperature and precipitation as its input forcing, and simulates, among others, snow water equivalent (SWE), snow depth (SD), and the snow bulk density (ρ). In this paper the set of equations contained in the seNorge model code is described and a thorough spatiotemporal statistical evaluation of the model performance from 1957–2011 is made using the two major sets of extensive in situ snow measurements that exist for Norway. The evaluation results show that the seNorge model generally overestimates both SWE and ρ, and that the overestimation of SWE increases with elevation throughout the snow season. However, the R 2 -values for model fit are 0.60 for (log-transformed) SWE and 0.45 for ρ, indicating that after removal of the detected systematic model biases (e.g. by recalibrating the model or expressing snow conditions in relative units) the model performs rather well. The seNorge model provides a relatively simple, not very data-demanding, yet nonetheless process-based method to construct snow maps of high spatiotemporal resolution. It is an especially well suited alternative for operational snow mapping in regions with rugged topography and large spatiotemporal variability in snow conditions, as is the case in the mountainous Norway.

Description: Ground penetrating radar detection of subsnow slush on ice-covered lakes in interior Alaska The Cryosphere, 6, 1435-1443, 2012 Author(s): A. Gusmeroli and G. Grosse Lakes are abundant throughout the pan-Arctic region. For many of these lakes ice cover lasts for up to two thirds of the year. The frozen cover allows human access to these lakes, which are therefore used for many subsistence and recreational activities, including water harvesting, fishing, and skiing. Safe traveling condition onto lakes may be compromised, however, when, after significant snowfall, the weight of the snow acts on the ice and causes liquid water to spill through weak spots and overflow at the snow-ice interface. Since visual detection of subsnow slush is almost impossible our understanding on overflow processes is still very limited and geophysical methods that allow water and slush detection are desirable. In this study we demonstrate that a commercially available, lightweight 1 GHz, ground penetrating radar system can detect and map extent and intensity of overflow. The strength of radar reflections from wet snow-ice interfaces are at least twice as much in strength than returns from dry snow-ice interface. The presence of overflow also affects the quality of radar returns from the base of the lake ice. During dry conditions we were able to profile ice thickness of up to 1 m, conversely, we did not retrieve any ice-water returns in areas affected by overflow.

Description: The footprint of Asian monsoon dynamics in the mass and energy balance of a Tibetan glacier The Cryosphere, 6, 1445-1461, 2012 Author(s): T. Mölg, F. Maussion, W. Yang, and D. Scherer Determinations of glacier-wide mass and energy balance are still scarce for the remote mountains of the Tibetan Plateau, where field measurements are challenging. Here we run and evaluate a physical, distributed mass balance model for Zhadang Glacier (central Tibet, 30° N) based on in-situ measurements over 2009–2011 and an uncertainty estimate by Monte Carlo and ensemble strategies. The model application aims to provide the first quantification of how the Indian Summer Monsoon (ISM) impacts an entire glacier over the various stages of the monsoon's annual cycle. We find a strong and systematic ISM footprint on the interannual scale. Early (late) monsoon onset causes higher (lower) accumulation, and reduces (increases) the available energy for ablation primarily through changes in absorbed shortwave radiation. By contrast, only a weak footprint exists in the ISM cessation phase. Most striking though is the core monsoon season: local mass and energy balance variability is fully decoupled from the active/break cycle that defines large-scale atmospheric variability during the ISM. Our results demonstrate quantitatively that monsoon onset strongly affects the ablation season of glaciers in Tibet. However, we find no direct ISM impact on the glacier in the main monsoon season, which has not been acknowledged so far. This result also adds cryospheric evidence that, once the monsoon is in full swing, regional atmospheric variability prevails on the Tibetan Plateau in summer.

Description: The first complete inventory of the local \newline glaciers and ice caps on Greenland The Cryosphere, 6, 1483-1495, 2012 Author(s): P. Rastner, T. Bolch, N. Mölg, H. Machguth, R. Le Bris, and F. Paul Glacier inventories provide essential baseline information for the determination of water resources, glacier-specific changes in area and volume, climate change impacts as well as past, potential and future contribution of glaciers to sea-level rise. Although Greenland is heavily glacierised and thus highly relevant for all of the above points, a complete inventory of its glaciers was not available so far. Here we present the results and details of a new and complete inventory that has been compiled from more than 70 Landsat scenes (mostly acquired between 1999 and 2002) using semi-automated glacier mapping techniques. A digital elevation model (DEM) was used to derive drainage divides from watershed analysis and topographic attributes for each glacier entity. To serve the needs of different user communities, we assigned to each glacier one of three connectivity levels with the ice sheet (CL0, CL1, CL2; i.e. no, weak, and strong connection) to clearly, but still flexibly, distinguish the local glaciers and ice caps (GIC) from the ice sheet and its outlet glaciers. In total, we mapped ~ 20 300 glaciers larger than 0.05 km 2 (of which ~ 900 are marine terminating), covering an area of 130 076 ± 4032 km 2 , or 89 720 ± 2781 km 2 without the CL2 GIC. The latter value is about 50% higher than the mean value of more recent previous estimates. Glaciers smaller than 0.5 km 2 contribute only 1.5% to the total area but more than 50% (11 000) to the total number. In contrast, the 25 largest GIC (〉 500 km 2 ) contribute 28% to the total area, but only 0.1% to the total number. The mean elevation of the GIC is 1700 m in the eastern sector and around 1000 m otherwise. The median elevation increases with distance from the coast, but has only a weak dependence on mean glacier aspect.

Description: Bedmap2: improved ice bed, surface and thickness datasets for Antarctica The Cryosphere, 7, 375-393, 2013 Author(s): P. Fretwell, H. D. Pritchard, D. G. Vaughan, J. L. Bamber, N. E. Barrand, R. Bell, C. Bianchi, R. G. Bingham, D. D. Blankenship, G. Casassa, G. Catania, D. Callens, H. Conway, A. J. Cook, H. F. J. Corr, D. Damaske, V. Damm, F. Ferraccioli, R. Forsberg, S. Fujita, Y. Gim, P. Gogineni, J. A. Griggs, R. C. A. Hindmarsh, P. Holmlund, J. W. Holt, R. W. Jacobel, A. Jenkins, W. Jokat, T. Jordan, E. C. King, J. Kohler, W. Krabill, M. Riger-Kusk, K. A. Langley, G. Leitchenkov, C. Leuschen, B. P. Luyendyk, K. Matsuoka, J. Mouginot, F. O. Nitsche, Y. Nogi, O. A. Nost, S. V. Popov, E. Rignot, D. M. Rippin, A. Rivera, J. Roberts, N. Ross, M. J. Siegert, A. M. Smith, D. Steinhage, M. Studinger, B. Sun, B. K. Tinto, B. C. Welch, D. Wilson, D. A. Young, C. Xiangbin, and A. Zirizzotti We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km 3 ) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.

Description: Variability and trends in Laptev Sea ice outflow between 1992–2011 The Cryosphere, 7, 349-363, 2013 Author(s): T. Krumpen, M. Janout, K. I. Hodges, R. Gerdes, F. Girard-Ardhuin, J. A. Hölemann, and S. Willmes Variability and trends in seasonal and interannual ice area export out of the Laptev Sea between 1992 and 2011 are investigated using satellite-based sea ice drift and concentration data. We found an average total winter (October to May) ice area transport across the northern and eastern Laptev Sea boundaries (NB and EB) of 3.48 × 10 5 km 2 . The average transport across the NB (2.87 × 10 5 km 2 ) is thereby higher than across the EB (0.61 × 10 5 km 2 ), with a less pronounced seasonal cycle. The total Laptev Sea ice area flux significantly increased over the last decades (0.85 × 10 5 km 2 decade −1 , p 〉 0.95), dominated by increasing export through the EB (0.55 × 10 5 km 2 decade −1 , p 〉 0.90), while the increase in export across the NB is smaller (0.3 × 10 5 km 2 decade −1 ) and statistically not significant. The strong coupling between across-boundary SLP gradient and ice drift velocity indicates that monthly variations in ice area flux are primarily controlled by changes in geostrophic wind velocities, although the Laptev Sea ice circulation shows no clear relationship with large-scale atmospheric indices. Also there is no evidence of increasing wind velocities that could explain the overall positive trends in ice export. The increased transport rates are rather the consequence of a changing ice cover such as thinning and/or a decrease in concentration. The use of a back-propagation method revealed that most of the ice that is incorporated into the Transpolar Drift is formed during freeze-up and originates from the central and western part of the Laptev Sea, while the exchange with the East Siberian Sea is dominated by ice coming from the central and southeastern Laptev Sea. Furthermore, our results imply that years of high ice export in late winter (February to May) have a thinning effect on the ice cover, which in turn preconditions the occurence of negative sea ice extent anomalies in summer.

Description: Retention and radiative forcing of black carbon in eastern Sierra Nevada snow The Cryosphere, 7, 365-374, 2013 Author(s): K. M. Sterle, J. R. McConnell, J. Dozier, R. Edwards, and M. G. Flanner When contaminated by absorbing particles, such as refractory black carbon (rBC) and continental dust, snow's albedo decreases and thus its absorption of solar radiation increases, thereby hastening snowmelt. For this reason, an understanding of rBC's affect on snow albedo, melt processes, and radiation balance is critical for water management, especially in a changing climate. Measurements of rBC in a sequence of snow pits and surface snow samples in the eastern Sierra Nevada of California during the snow accumulation and ablation seasons of 2009 show that concentrations of rBC were enhanced sevenfold in surface snow (~25 ng g –1 ) compared to bulk values in the snowpack (~3 ng g –1 ). Unlike major ions, which were preferentially released during the initial melt, rBC and continental dust were retained in the snow, enhancing concentrations well into late spring, until a final flush occurred during the ablation period. We estimate a combined rBC and continental dust surface radiative forcing of 20 to 40 W m −2 during April and May, with dust likely contributing a greater share of the forcing.

Description: Surface undulations of Antarctic ice streams tightly controlled by bedrock topography The Cryosphere, 7, 407-417, 2013 Author(s): J. De Rydt, G. H. Gudmundsson, H. F. J. Corr, and P. Christoffersen Full Stokes flow-line models predict that fast-flowing ice streams transmit information about their bedrock topography most efficiently to the surface for basal undulations with length scales between 1 and 20 times the mean ice thickness. This typical behaviour is independent of the precise values of the flow law and sliding law exponents, and should be universally observable. However, no experimental evidence for this important theoretical prediction has been obtained so far, hence ignoring an important test for the physical validity of current-day ice flow models. In our work we use recently acquired airborne radar data for the Rutford Ice Stream and Evans Ice Stream, and we show that the surface response of fast-flowing ice is highly sensitive to bedrock irregularities with wavelengths of several ice thicknesses. The sensitivity depends on the slip ratio, i.e. the ratio between mean basal sliding velocity and mean deformational velocity. We find that higher values of the slip ratio generally lead to a more efficient transfer, whereas the transfer is significantly dampened for ice that attains most of its surface velocity by creep. Our findings underline the importance of bedrock topography for ice stream dynamics on spatial scales up to 20 times the mean ice thickness. Our results also suggest that local variations in the flow regime and surface topography at this spatial scale cannot be explained by variations in basal slipperiness.

Description: Grounding line transient response in marine ice sheet models The Cryosphere, 7, 395-406, 2013 Author(s): A. S. Drouet, D. Docquier, G. Durand, R. Hindmarsh, F. Pattyn, O. Gagliardini, and T. Zwinger Marine ice-sheet stability is mostly controlled by the dynamics of the grounding line, i.e. the junction between the grounded ice sheet and the floating ice shelf. Grounding line migration has been investigated within the framework of MISMIP (Marine Ice Sheet Model Intercomparison Project), which mainly aimed at investigating steady state solutions. Here we focus on transient behaviour, executing short-term simulations (200 yr) of a steady ice sheet perturbed by the release of the buttressing restraint exerted by the ice shelf on the grounded ice upstream. The transient grounding line behaviour of four different flowline ice-sheet models has been compared. The models differ in the physics implemented (full Stokes and shallow shelf approximation), the numerical approach, as well as the grounding line treatment. Their overall response to the loss of buttressing is found to be broadly consistent in terms of grounding line position, rate of surface elevation change and surface velocity. However, still small differences appear for these latter variables, and they can lead to large discrepancies (〉 100%) observed in terms of ice sheet contribution to sea level when cumulated over time. Despite the recent important improvements of marine ice-sheet models in their ability to compute steady state configurations, our results question the capacity of these models to compute short-term reliable sea-level rise projections.

Description: Glacier changes from 1966–2009 in the Gongga Mountains, on the south-eastern margin of the Qinghai-Tibetan Plateau and their climatic forcing The Cryosphere, 6, 1087-1101, 2012 Author(s): B. T. Pan, G. L. Zhang, J. Wang, B. Cao, H. P. Geng, J. Wang, C. Zhang, and Y. P. Ji In order to monitor the changes of the glaciers in the Gongga Mountain region on the south-eastern margin of the Qinghai-Tibetan Plateau, 74 monsoonal temperate glaciers were investigated by comparing the Chinese Glacier Inventory (CGI), recorded in the 1960s, with Landsat MSS in 1974, Landsat TM in 1989, 1994, 2005, and ASTER data in 2009. The remote sensing data have been applied to map the glacier outline by threshold ratio images (TM4/TM5). Moreover, the glacier outlines were verified by GPS survey on four large glaciers (Hailuogou (HLG), Mozigou (MZG), Yanzigou (YZG), and Dagongba (DGB)) in 2009. The results show that the area dominated by the 74 glaciers has shrunk by 11.3% (29.2 km 2 ) from 1966 to 2009. Glacier area on the eastern and western slopes of the Gongga Mountains decreased by 9.8% and 14.6% since 1966, respectively. The loss in glacier area and length is, respectively, 0.8 km 2 and 1146.4 m for the HLG Glacier, 2.1 km 2 and 501.8 m for the MZG Glacier, 0.8 km 2 and 724.8 m for the YZG Glacier, and 2.4 km 2 and 1002.3 m for the DGB Glacier. Decades of climate records obtained from three meteorological stations in the Gongga Mountains were analyzed to evaluate the impact of the temperature and precipitation on glacier retreat. The mean annual temperatures over the eastern and western slopes of the Gongga Mountains have been increasing by 0.34 K decade −1 and 0.24 K decade −1 (1988–2009), respectively. Moreover, mean annual precipitation has only increased by 1% in the past 50 yr. The increasing amount of precipitation could not compensate for the glacier mass loss due to the temperature increase in the Gongga Mountains. This suggests that the warming of the climate is probably also responsible for the glacier retreat in the study region. At the region scale, glacier changes were also controlled by local topographical factors.

Description: Accelerated contributions of Canada's Baffin and Bylot Island glaciers to sea level rise over the past half century The Cryosphere, 6, 1103-1125, 2012 Author(s): A. Gardner, G. Moholdt, A. Arendt, and B. Wouters Canadian Arctic glaciers have recently contributed large volumes of meltwater to the world's oceans. To place recently observed glacier wastage into a historical perspective and to determine the region's longer-term (~50 years) contribution to sea level, we estimate mass and volume changes for the glaciers of Baffin and Bylot Islands using digital elevation models generated from airborne and satellite stereoscopic imagery and elevation postings from repeat airborne and satellite laser altimetry. In addition, we update existing glacier mass change records from GRACE satellite gravimetry to cover the period from 2003 to 2011. Using this integrated approach, we find that the rate of mass loss from the region's glaciers increased from 11.1 ± 3.4 Gt a −1 (271 ± 84 kg m −2 a −1 ) for the period 1963–2006 to 23.8 ± 6.1 Gt a −1 (581 ± 149 kg m −2 a −1 ) for the period 2003–2011. The doubling of the rate of mass loss is attributed to higher temperatures in summer with little change in annual precipitation. Through both direct and indirect effects, changes in summer temperatures accounted for 70–98% of the variance in the rate of mass loss, to which the Barnes Ice Cap was found to be 1.7 times more sensitive than either the Penny Ice Cap or the region's glaciers as a whole. This heightened sensitivity is the result of a glacier hypsometry that is skewed to lower elevations, which are shown to have a higher mass change sensitive to temperature compared to glacier surfaces at higher elevations. Between 2003 and 2011 the glaciers of Baffin and Bylot Islands contributed 0.07 ± 0.02 mm a −1 to sea level rise accounting for 16% of the total contribution from glaciers outside of Greenland and Antarctica, a rate much higher than the longer-term average of 0.03 ± 0.01 mm a −1 (1963 to 2006).

Description: Ice velocity changes in the Ross and Ronne sectors observed using satellite radar data from 1997 and 2009 The Cryosphere, 6, 1019-1030, 2012 Author(s): B. Scheuchl, J. Mouginot, and E. Rignot We report changes in ice velocity of a 6.5 million km 2 region around South Pole encompassing the Filchner-Ronne and Ross Ice Shelves and a significant portion of the ice streams and glaciers that constitute their catchment areas. Using the first full interferometric synthetic aperture radar (InSAR) coverage of the region completed in 2009 and partial coverage acquired in 1997, we processed the data to assemble a comprehensive map of ice speed changes between those two years. On the Ross Ice Shelf, our results confirm a continued deceleration of Mercer and Whillans Ice Streams with a 12-yr velocity difference of −50 m yr −1 (−16.7%) and −100 m yr −1 (−25.3%) at their grounding lines. The deceleration spreads 450 km upstream of the grounding line and more than 500 km onto the shelf, beyond what was previously known. Ross and Filchner Ice Shelves exhibit signs of pre-calving events, representing the largest observed changes, with an increase in speed in excess of +100 m yr −1 in 12 yr. Other changes in the Ross Ice Shelf region are less significant. The observed changes in glacier speed extend on the Ross Ice Shelf along the ice streams' flow lines. Most tributaries of the Filchner-Ronne Ice Shelf show a modest deceleration or no change between 1997 and 2009. Slessor Glacier shows a small deceleration over a large sector. No change is detected on the Bailey, Rutford, and Institute Ice Streams. On the Filchner Ice Shelf itself, ice decelerated rather uniformly with a 12-yr difference in speed of −50 m yr −1 , or −5% of its ice front speed, which we attribute to a 12 km advance in its ice front position. Our results show that dynamic changes are present in the region. They highlight the need for continued observation of the area with a primary focus on the Siple Coast. The dynamic changes in Central Antarctica between 1997 and 2009 are generally second-order effects in comparison to losses on glaciers in the Bellingshausen and Amundsen Seas region and on the Antarctic Peninsula. We therefore conclude that the dynamic changes shown here do not have a strong impact on the mass budget of the Antarctic continent.

Description: Uncertainties in the global temperature change caused by carbon release from permafrost thawing The Cryosphere, 6, 1063-1076, 2012 Author(s): E. J. Burke, I. P. Hartley, and C. D. Jones Under climate change thawing permafrost will cause old carbon which is currently frozen and inert to become vulnerable to decomposition and release into the climate system. This paper develops a simple framework for estimating the impact of this permafrost carbon release on the global mean temperature (P-GMT). The analysis is based on simulations made with the Hadley Centre climate model (HadGEM2-ES) for a range of representative CO 2 concentration pathways. Results using the high concentration pathway (RCP 8.5) suggest that by 2100 the annual methane (CH 4 ) emission rate is 2–59 Tg CH 4 yr −1 and 50–270 Pg C has been released as CO 2 with an associated P-GMT of 0.08–0.36 °C (all 5th–95th percentile ranges). P-GMT is considerably lower – between 0.02 and 0.11 °C – for the low concentration pathway (RCP2.6). The uncertainty in climate model scenario causes about 50% of the spread in P-GMT by the end of the 21st century. The distribution of soil carbon, in particular how it varies with depth, contributes to about half of the remaining spread, with quality of soil carbon and decomposition processes contributing a further quarter each. These latter uncertainties could be reduced through additional observations. Over the next 20–30 yr, whilst scenario uncertainty is small, improving our knowledge of the quality of soil carbon will contribute significantly to reducing the spread in the, albeit relatively small, P-GMT.

Description: Refreezing on the Greenland ice sheet: a comparison of parameterizations The Cryosphere, 6, 743-762, 2012 Author(s): C. H. Reijmer, M. R. van den Broeke, X. Fettweis, J. Ettema, and L. B. Stap Retention and refreezing of meltwater are acknowledged to be important processes for the mass budget of polar glaciers and ice sheets. Several parameterizations of these processes exist for use in energy and mass balance models. Due to a lack of direct observations, validation of these parameterizations is difficult. In this study we compare a set of 6 refreezing parameterizations against output of two Regional Climate Models (RCMs) coupled to an energy balance snow model, the Regional Atmospheric Climate Model (RACMO2) and the Modèle Atmosphérique Régional (MAR), applied to the Greenland ice sheet. In both RCMs, refreezing is explicitly calculated in a snow model that calculates vertical profiles of temperature, density and liquid water content. Between RACMO2 and MAR, the ice sheet-integrated amount of refreezing differs by only 4.9 mm w.e yr −1 (4.5 %), and the temporal and spatial variability are very similar. For consistency, the parameterizations are forced with output (surface temperature, precipitation and melt) of the RCMs. For the ice sheet-integrated amount of refreezing and its inter-annual variations, all parameterizations give similar results, especially after some tuning. However, the spatial distributions differ significantly and the spatial correspondence between the RCMs is better than with any of the parameterizations. Results are especially sensitive to the choice of the depth of the thermally active layer, which determines the cold content of the snow in most parameterizations. These results are independent of which RCM is used to force the parameterizations.

Description: Vapor flux and recrystallization during dry snow metamorphism under a steady temperature gradient as observed by time-lapse micro-tomography The Cryosphere, 6, 1141-1155, 2012 Author(s): B. R. Pinzer, M. Schneebeli, and T. U. Kaempfer Dry snow metamorphism under an external temperature gradient is the most common type of recrystallization of snow on the ground. The changes in snow microstructure modify the physical properties of snow, and therefore an understanding of this process is essential for many disciplines, from modeling the effects of snow on climate to assessing avalanche risk. We directly imaged the microstructural changes in snow during temperature gradient metamorphism (TGM) under a constant gradient of 50 K m −1 , using in situ time-lapse X-ray micro-tomography. This novel and non-destructive technique directly reveals the amount of ice that sublimates and is deposited during metamorphism, in addition to the exact locations of these phase changes. We calculated the average time that an ice volume stayed in place before it sublimated and found a characteristic residence time of 2–3 days. This means that most of the ice changes its phase from solid to vapor and back many times in a seasonal snowpack where similar temperature conditions can be found. Consistent with such a short timescale, we observed a mass turnover of up to 60% of the total ice mass per day. The concept of hand-to-hand transport for the water vapor flux describes the observed changes very well. However, we did not find evidence for a macroscopic vapor diffusion enhancement. The picture of {temperature gradient metamorphism} produced by directly observing the changing microstructure sheds light on the micro-physical processes and could help to improve models that predict the physical properties of snow.

Description: Transition in the fractal geometry of Arctic melt ponds The Cryosphere, 6, 1157-1162, 2012 Author(s): C. Hohenegger, B. Alali, K. R. Steffen, D. K. Perovich, and K. M. Golden During the Arctic melt season, the sea ice surface undergoes a remarkable transformation from vast expanses of snow covered ice to complex mosaics of ice and melt ponds. Sea ice albedo, a key parameter in climate modeling, is determined by the complex evolution of melt pond configurations. In fact, ice–albedo feedback has played a major role in the recent declines of the summer Arctic sea ice pack. However, understanding melt pond evolution remains a significant challenge to improving climate projections. By analyzing area–perimeter data from hundreds of thousands of melt ponds, we find here an unexpected separation of scales, where pond fractal dimension D transitions from 1 to 2 around a critical length scale of 100 m 2 in area. Pond complexity increases rapidly through the transition as smaller ponds coalesce to form large connected regions, and reaches a maximum for ponds larger than 1000 m 2 , whose boundaries resemble space-filling curves, with D ≈ 2. These universal features of Arctic melt pond evolution are similar to phase transitions in statistical physics. The results impact sea ice albedo, the transmitted radiation fields under melting sea ice, the heat balance of sea ice and the upper ocean, and biological productivity such as under ice phytoplankton blooms.

Description: Albedo of the ice covered Weddell and Bellingshausen Seas The Cryosphere, 6, 479-491, 2012 Author(s): A. I. Weiss, J. C. King, T. A. Lachlan-Cope, and R. S. Ladkin This study investigates the surface albedo of the sea ice areas adjacent to the Antarctic Peninsula during the austral summer. Aircraft measurements of the surface albedo, which were conducted in the sea ice areas of the Weddell and Bellingshausen Seas show significant differences between these two regions. The averaged surface albedo varied between 0.13 and 0.81. The ice cover of the Bellingshausen Sea consisted mainly of first year ice and the sea surface showed an averaged sea ice albedo of α i = 0.64 ± 0.2 (± standard deviation). The mean sea ice albedo of the pack ice area in the western Weddell Sea was α i = 0.75 ± 0.05. In the southern Weddell Sea, where new, young sea ice prevailed, a mean albedo value of α i = 0.38 ± 0.08 was observed. Relatively warm open water and thin, newly formed ice had the lowest albedo values, whereas relatively cold and snow covered pack ice had the highest albedo values. All sea ice areas consisted of a mixture of a large range of different sea ice types. An investigation of commonly used parameterizations of albedo as a function of surface temperature in the Weddell and Bellingshausen Sea ice areas showed that the albedo parameterizations do not work well for areas with new, young ice.

Description: Seasonal speed-up of two outlet glaciers of Austfonna, Svalbard, inferred from continuous GPS measurements The Cryosphere, 6, 453-466, 2012 Author(s): T. Dunse, T. V. Schuler, J. O. Hagen, and C. H. Reijmer A large part of the ice discharge from ice caps and ice sheets occurs through spatially limited flow units that may operate in a mode of steady flow or cyclic surge behaviour. Changes in the dynamics of distinct flow units play a key role in the mass balance of Austfonna, the largest ice cap on Svalbard. The recent net mass loss of Austfonna was dominated by calving from marine terminating outlet glaciers. Previous ice-surface velocity maps of the ice cap were derived by satellite radar interferometry (InSAR) and rely on data acquired in the mid-1990s with limited information concerning the temporal variability. Here, we present continuous Global Positioning System (GPS) observations along the central flowlines of two fast flowing outlet glaciers over 2008–2010. The data show prominent summer speed-ups with ice-surface velocities as high as 240% of the pre-summer mean. Acceleration follows the onset of the summer melt period, indicating enhanced basal motion due to input of surface meltwater into the subglacial drainage system. In 2008, multiple velocity peaks coincide with successive melt periods. In 2009, the major melt was of higher amplitude than in 2008. Flow velocities appear unaffected by subsequent melt periods, suggesting a transition towards a hydraulically more efficient drainage system. The observed annual mean velocities of Duvebreen and Basin-3 exceed those from the mid-1990s by factors two and four, respectively, implying increased ice discharge at the calving front. Measured summer velocities up to 2 m d −1 for Basin-3 are close to those of Kronebreen, often referred to as the fastest glacier on Svalbard.

Description: Repeat optical satellite images reveal widespread and long term decrease in land-terminating glacier speeds The Cryosphere, 6, 467-478, 2012 Author(s): T. Heid and A. Kääb By matching of repeat optical satellite images it is now possible to investigate glacier dynamics within large regions of the world and also between regions to improve knowledge about glacier dynamics in space and time. In this study we investigate whether the negative glacier mass balance seen over large parts of the world has caused the glaciers to change their speeds. The studied regions are Pamir, Caucasus, Penny Ice Cap, Alaska Range and Patagonia. In addition we derive speed changes for Karakoram, a region assumed to have positive mass balance and that contains many surge-type glaciers. We find that the mapped glaciers in the five regions with negative mass balance have over the last decades decreased their velocity at an average rate per decade of: 43 % in the Pamir, 8 % in the Caucasus, 25 % on Penny Ice Cap, 11 % in the Alaska Range and 20 % in Patagonia. Glaciers in Karakoram have generally increased their speeds, but surging glaciers and glaciers with flow instabilities are most prominent in this area. Therefore the calculated average speed change is not representative for this area.

Description: Thermal conductivity of snow measured by three independent methods and anisotropy considerations The Cryosphere, 7, 217-227, 2013 Author(s): F. Riche and M. Schneebeli The thermal conductivity of snow determines the temperature gradient, and by this, it has a direct effect on the rate of snow metamorphism. It is therefore a key property of snow. However, thermal conductivities measured with the transient needle probe and the steady-state, heat flux plate differ. In addition, the anisotropy of thermal conductivity plays an important role in the accuracy of thermal conductivity measurements. In this study, we investigated three independent methods to measure snow thermal conductivity and its anisotropy: a needle probe with a long heating time, a guarded heat flux plate, and direct numerical simulation at the microstructural level of the pore and ice structure. The three methods were applied to identical snow samples. We analyzed the consistency and the difference between these methods. As already shown in former studies, we observed a distinct difference between the anisotropy of thermal conductivity in small rounded grains and in depth hoar. Indeed, the anisotropy between vertical and horizontal thermal conductivity components ranges between 0.5–2. This can cause a difference in thermal conductivity measurements carried out with needle probes of up to –25 % to +25 % if the thermal conductivity is calculated only from a horizontally inserted needle probe. Based on our measurements and the comparison of the three methods studied here, the direct numerical simulation is the most reliable method, as the tensorial components of the thermal conductivity can be calculated and the corresponding microstructure is precisely known.

Description: Brief communication "Important role of the mid-tropospheric atmospheric circulation in the recent surface melt increase over the Greenland ice sheet" The Cryosphere, 7, 241-248, 2013 Author(s): X. Fettweis, E. Hanna, C. Lang, A. Belleflamme, M. Erpicum, and H. Gallée Since 2007, there has been a series of surface melt records over the Greenland ice sheet (GrIS), continuing the trend towards increased melt observed since the end of the 1990's. The last two decades are characterized by an increase of negative phases of the North Atlantic Oscillation (NAO) favouring warmer and drier summers than normal over GrIS. In this context, we use a circulation type classification based on daily 500 hPa geopotential height to evaluate the role of atmospheric dynamics in this surface melt acceleration for the last two decades. Due to the lack of direct observations, the interannual melt variability is gauged here by the summer (June–July–August) mean temperature from reanalyses at 700 hPa over Greenland; analogous atmospheric circulations in the past show that ~70% of the 1993–2012 warming at 700 hPa over Greenland has been driven by changes in the atmospheric flow frequencies. Indeed, the occurrence of anticyclones centred over the GrIS at the surface and at 500 hPa has doubled since the end of 1990's, which induces more frequent southerly warm air advection along the western Greenland coast and over the neighbouring Canadian Arctic Archipelago (CAA). These changes in the NAO modes explain also why no significant warming has been observed these last summers over Svalbard, where northerly atmospheric flows are twice as frequent as before. Therefore, the recent warmer summers over GrIS and CAA cannot be considered as a long-term climate warming but are more a consequence of NAO variability affecting atmospheric heat transport. Although no global model from the CMIP5 database projects subsequent significant changes in NAO through this century, we cannot exclude the possibility that the observed NAO changes are due to global warming.

Description: Restoring mass conservation to shallow ice flow models over complex terrain The Cryosphere, 7, 229-240, 2013 Author(s): A. H. Jarosch, C. G. Schoof, and F. S. Anslow Numerical simulation of glacier dynamics in mountainous regions using zero-order, shallow ice models is desirable for computational efficiency so as to allow broad coverage. However, these models present several difficulties when applied to complex terrain. One such problem arises where steep terrain can spuriously lead to large ice fluxes that remove more mass from a grid cell than it originally contains, leading to mass conservation being violated. This paper describes a vertically integrated, shallow ice model using a second-order flux-limiting spatial discretization scheme that enforces mass conservation. An exact solution to ice flow over a bedrock step is derived for a given mass balance forcing as a benchmark to evaluate the model performance in such a difficult setting. This benchmark should serve as a useful test for modellers interested in simulating glaciers over complex terrain.

Description: Borehole temperatures reveal a changed energy budget at Mill Island, East Antarctica, over recent decades The Cryosphere, 7, 263-273, 2013 Author(s): J. L. Roberts, A. D. Moy, T. D. van Ommen, M. A. J. Curran, A. P. Worby, I. D. Goodwin, and M. Inoue A borehole temperature record from the Mill Island (East Antarctica) icecap reveals a large surface warming signal manifested as a 0.75 K temperature difference over the approximate 100 m depth in the zone of zero annual amplitude below the seasonally varying zone. The temperature profile shows a break in gradient around 49 m depth, which we model with inverse numerical simulations, indicating that surface warming started around the austral summer of 1980/81 AD ±5 yr. This warming of approximately 0.37 K per decade is consistent with trends seen in both instrumental and other reconstructions for Antarctica and, therefore, suggests that regional- rather than local-scale processes are largely responsible. Alteration of the surface energy budget arising from changes in radiation balances due to local cloud, the amount of liquid deposition and local air temperatures associated with altered air/sea exchanges also potentially plays a role at this location due to the proximity of the Shackleton Ice Shelf and sea-ice zone.

Description: Analysis of ice phenology of lakes on the Tibetan Plateau from MODIS data The Cryosphere, 7, 287-301, 2013 Author(s): J. Kropáček, F. Maussion, F. Chen, S. Hoerz, and V. Hochschild The Tibetan Plateau includes a large system of endorheic (closed basin) lakes. Lake ice phenology, i.e. the timing of freeze-up and break-up and the duration of the ice cover may provide valuable information about climate variations in this region. The ice phenology of 59 large lakes on the Tibetan Plateau was derived from Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day composite data for the period from 2001 to 2010. Ice cover duration appears to have a high variability in the studied region due to both climatic and local factors. Mean values for the duration of ice cover were calculated for three groups of lakes defined by clustering, resulting in relatively compact geographic regions. In each group several lakes showed anomalies in ice cover duration in the studied period. Possible reasons for such anomalous behaviour are discussed. Furthermore, many lakes do not freeze up completely during some seasons. This was confirmed by inspection of high resolution optical data. Mild winter seasons, large water volume and/or high salinity are the most likely explanations. Trends in the ice cover duration derived by linear regression for all the studied lakes show a high variation in space. A correlation of ice phenology variables with parameters describing climatic and local conditions showed a high thermal dependency of the ice regime. It appears that the freeze-up tends to be more thermally determined than break-up for the studied lakes.

Description: Paleo ice flow and subglacial meltwater dynamics in Pine Island Bay, West Antarctica The Cryosphere, 7, 249-262, 2013 Author(s): F. O. Nitsche, K. Gohl, R. D. Larter, C.-D. Hillenbrand, G. Kuhn, J. A. Smith, S. Jacobs, J. B. Anderson, and M. Jakobsson Increasing evidence for an elaborate subglacial drainage network underneath modern Antarctic ice sheets suggests that basal meltwater has an important influence on ice stream flow. Swath bathymetry surveys from previously glaciated continental margins display morphological features indicative of subglacial meltwater flow in inner shelf areas of some paleo ice stream troughs. Over the last few years several expeditions to the eastern Amundsen Sea embayment (West Antarctica) have investigated the paleo ice streams that extended from the Pine Island and Thwaites glaciers. A compilation of high-resolution swath bathymetry data from inner Pine Island Bay reveals details of a rough seabed topography including several deep channels that connect a series of basins. This complex basin and channel network is indicative of meltwater flow beneath the paleo-Pine Island and Thwaites ice streams, along with substantial subglacial water inflow from the east. This meltwater could have enhanced ice flow over the rough bedrock topography. Meltwater features diminish with the onset of linear features north of the basins. Similar features have previously been observed in several other areas, including the Dotson-Getz Trough (western Amundsen Sea embayment) and Marguerite Bay (SW Antarctic Peninsula), suggesting that these features may be widespread around the Antarctic margin and that subglacial meltwater drainage played a major role in past ice-sheet dynamics.

Description: A recent tipping point in the Arctic sea-ice cover: abrupt and persistent increase in the seasonal cycle since 2007 The Cryosphere, 7, 275-286, 2013 Author(s): V. N. Livina and T. M. Lenton There is ongoing debate over whether Arctic sea ice has already passed a "tipping point", or whether it will do so in the future. Several recent studies argue that the loss of summer sea ice does not involve an irreversible bifurcation, because it is highly reversible in models. However, a broader definition of a "tipping point" also includes other abrupt, non-linear changes that are neither bifurcations nor necessarily irreversible. Examination of satellite data for Arctic sea-ice area reveals an abrupt increase in the amplitude of seasonal variability in 2007 that has persisted since then. We identified this abrupt transition using recently developed methods that can detect multi-modality in time-series data and sometimes forewarn of bifurcations. When removing the mean seasonal cycle (up to 2008) from the satellite data, the residual sea-ice fluctuations switch from uni-modal to multi-modal behaviour around 2007. We originally interpreted this as a bifurcation in which a new lower ice cover attractor appears in deseasonalised fluctuations and is sampled in every summer–autumn from 2007 onwards. However, this interpretation is clearly sensitive to how the seasonal cycle is removed from the raw data, and to the presence of continental land masses restricting winter–spring ice fluctuations. Furthermore, there was no robust early warning signal of critical slowing down prior to the hypothesized bifurcation. Early warning indicators do however show destabilization of the summer–autumn sea-ice cover since 2007. Thus, the bifurcation hypothesis lacks consistent support, but there was an abrupt and persistent increase in the amplitude of the seasonal cycle of Arctic sea-ice cover in 2007, which we describe as a (non-bifurcation) "tipping point". Our statistical methods detect this "tipping point" and its time of onset. We discuss potential geophysical mechanisms behind it, which should be the subject of further work with process-based models.

Description: Analysis of the snow-atmosphere energy balance during wet-snow instabilities and implications for avalanche prediction The Cryosphere, 7, 205-216, 2013 Author(s): C. Mitterer and J. Schweizer Wet-snow avalanches are notoriously difficult to predict; their formation mechanism is poorly understood since in situ measurements representing the thermal and mechanical evolution are difficult to perform. Instead, air temperature is commonly used as a predictor variable for days with high wet-snow avalanche danger – often with limited success. As melt water is a major driver of wet-snow instability and snow melt depends on the energy input into the snow cover, we computed the energy balance for predicting periods with high wet-snow avalanche activity. The energy balance was partly measured and partly modelled for virtual slopes at different elevations for the aspects south and north using the 1-D snow cover model SNOWPACK. We used measured meteorological variables and computed energy balance and its components to compare wet-snow avalanche days to non-avalanche days for four consecutive winter seasons in the surroundings of Davos, Switzerland. Air temperature, the net shortwave radiation and the energy input integrated over 3 or 5 days showed best results in discriminating event from non-event days. Multivariate statistics, however, revealed that for better predicting avalanche days, information on the cold content of the snowpack is necessary. Wet-snow avalanche activity was closely related to periods when large parts of the snowpack reached an isothermal state (0 °C) and energy input exceeded a maximum value of 200 kJ m −2 in one day, or the 3-day sum of positive energy input was larger than 1.2 MJ m −2 . Prediction accuracy with measured meteorological variables was as good as with computed energy balance parameters, but simulated energy balance variables accounted better for different aspects, slopes and elevations than meteorological data.

Description: Effect of higher-order stress gradients on the centennial mass evolution of the Greenland ice sheet The Cryosphere, 7, 183-199, 2013 Author(s): J. J. Fürst, H. Goelzer, and P. Huybrechts We use a three-dimensional thermo-mechanically coupled model of the Greenland ice sheet to assess the effects of marginal perturbations on volume changes on centennial timescales. The model is designed to allow for five ice dynamic formulations using different approximations to the force balance. The standard model is based on the shallow ice approximation for both ice deformation and basal sliding. A second model version relies on a higher-order Blatter/Pattyn type of core that resolves effects from gradients in longitudinal stresses and transverse horizontal shearing, i.e. membrane-like stresses. Together with three intermediate model versions, these five versions allow for gradually more dynamic feedbacks from membrane stresses. Idealised experiments are conducted on various resolutions to compare the time-dependent response to imposed accelerations at the marine ice front. If such marginal accelerations are to have an appreciable effect on total mass loss on a century timescale, a fast mechanism to transmit such perturbations inland is required. While the forcing is independent of the model version, inclusion of direct horizontal coupling allows the initial speed-up to reach several tens of kilometres inland. Within one century, effects from gradients in membrane stress alter the inland signal propagation and transmit additional dynamic thinning to the ice sheet interior. But the centennial overall volume loss differs only by some percents from the standard model, as the dominant response is a diffusive inland propagation of geometric changes. For the experiments considered, this volume response is even attenuated by direct horizontal coupling. The reason is a faster adjustment of the sliding regime by instant stress transmission in models that account for the effect of membrane stresses. Ultimately, horizontal coupling decreases the reaction time to perturbations at the ice sheet margin. These findings suggest that for modelling the mass evolution of a large-scale ice sheet, effects from diffusive geometric adjustments dominate effects from successively more complete dynamic approaches.

Description: Remote sensing of sea ice: advances during the DAMOCLES project The Cryosphere, 6, 1411-1434, 2012 Author(s): G. Heygster, V. Alexandrov, G. Dybkjær, W. von Hoyningen-Huene, F. Girard-Ardhuin, I. L. Katsev, A. Kokhanovsky, T. Lavergne, A. V. Malinka, C. Melsheimer, L. Toudal Pedersen, A. S. Prikhach, R. Saldo, R. Tonboe, H. Wiebe, and E. P. Zege In the Arctic, global warming is particularly pronounced so that we need to monitor its development continuously. On the other hand, the vast and hostile conditions make in situ observation difficult, so that available satellite observations should be exploited in the best possible way to extract geophysical information. Here, we give a résumé of the sea ice remote sensing efforts of the European Union's (EU) project DAMOCLES (Developing Arctic Modeling and Observing Capabilities for Long-term Environmental Studies). In order to better understand the seasonal variation of the microwave emission of sea ice observed from space, the monthly variations of the microwave emissivity of first-year and multi-year sea ice have been derived for the frequencies of the microwave imagers like AMSR-E (Advanced Microwave Scanning Radiometer on EOS) and sounding frequencies of AMSU (Advanced Microwave Sounding Unit), and have been used to develop an optimal estimation method to retrieve sea ice and atmospheric parameters simultaneously. In addition, a sea ice microwave emissivity model has been used together with a thermodynamic model to establish relations between the emissivities from 6 GHz to 50 GHz. At the latter frequency, the emissivity is needed for assimilation into atmospheric circulation models, but is more difficult to observe directly. The size of the snow grains on top of the sea ice influences both its albedo and the microwave emission. A method to determine the effective size of the snow grains from observations in the visible range (MODIS) is developed and demonstrated in an application on the Ross ice shelf. The bidirectional reflectivity distribution function (BRDF) of snow, which is an essential input parameter to the retrieval, has been measured in situ on Svalbard during the DAMOCLES campaign, and a BRDF model assuming aspherical particles is developed. Sea ice drift and deformation is derived from satellite observations with the scatterometer ASCAT (62.5 km grid spacing), with visible AVHRR observations (20 km), with the synthetic aperture radar sensor ASAR (10 km), and a multi-sensor product (62.5 km) with improved angular resolution (Continuous Maximum Cross Correlation, CMCC method) is presented. CMCC is also used to derive the sea ice deformation, important for formation of sea ice leads (diverging deformation) and pressure ridges (converging). The indirect determination of sea ice thickness from altimeter freeboard data requires knowledge of the ice density and snow load on sea ice. The relation between freeboard and ice thickness is investigated based on the airborne Sever expeditions conducted between 1928 and 1993.

Description: Environmental controls on the thermal structure of alpine glaciers The Cryosphere, 7, 167-182, 2013 Author(s): N. J. Wilson and G. E. Flowers Water entrapped in glacier accumulation zones represents a significant latent heat contribution to the development of thermal structure. It also provides a direct link between glacier environments and thermal regimes. We apply a two-dimensional mechanically-coupled model of heat flow to synthetic glacier geometries in order to explore the environmental controls on flowband thermal structure. We use this model to test the sensitivity of thermal structure to physical and environmental variables and to explore glacier thermal response to environmental changes. In different conditions consistent with a warming climate, mean glacier temperature and the volume of temperate ice may either increase or decrease, depending on the competing effects of elevated meltwater production, reduced accumulation zone extent and thinning firn. For two model reference states that exhibit commonly-observed thermal structures, the fraction of temperate ice is shown to decline with warming air temperatures. Mass balance and aquifer sensitivities play an important role in determining how the englacial thermal regimes of alpine glaciers will adjust in the future.

Description: Brief Communication "Expansion of meltwater lakes on the Greenland Ice Sheet" The Cryosphere, 7, 201-204, 2013 Author(s): I. M. Howat, S. de la Peña, J. H. van Angelen, J. T. M. Lenaerts, and M. R. van den Broeke Forty years of satellite imagery reveal that meltwater lakes on the margin of the Greenland Ice Sheet have expanded substantially inland to higher elevations with warming. These lakes are important because they provide a mechanism for bringing water to the ice bed, causing sliding. Inland expansion of lakes could accelerate ice flow by bringing water to previously frozen bed, potentially increasing future rates of mass loss. Increasing lake elevations closely follow the rise of the mass balance equilibrium line over much of the ice sheet, suggesting no physical limit on lake expansion. Data are not yet available to detect a corresponding change in ice flow, and the potential effects of lake expansion on ice sheet dynamics are not included in ice sheet models.

Description: Ikaite crystals in melting sea ice – implications for p CO 2 and pH levels in Arctic surface waters The Cryosphere, 6, 901-908, 2012 Author(s): S. Rysgaard, R. N. Glud, K. Lennert, M. Cooper, N. Halden, R. J. G. Leakey, F. C. Hawthorne, and D. Barber A major issue of Arctic marine science is to understand whether the Arctic Ocean is, or will be, a source or sink for air–sea CO 2 exchange. This has been complicated by the recent discoveries of ikaite (a polymorph of CaCO 3 ·6H 2 O) in Arctic and Antarctic sea ice, which indicate that multiple chemical transformations occur in sea ice with a possible effect on CO 2 and pH conditions in surface waters. Here, we report on biogeochemical conditions, microscopic examinations and x-ray diffraction analysis of single crystals from a melting 1.7 km 2 (0.5–1 m thick) drifting ice floe in the Fram Strait during summer. Our findings show that ikaite crystals are present throughout the sea ice but with larger crystals appearing in the upper ice layers. Ikaite crystals placed at elevated temperatures disintegrated into smaller crystallites and dissolved. During our field campaign in late June, melt reduced the ice floe thickness by 0.2 m per week and resulted in an estimated 3.8 ppm decrease of p CO 2 in the ocean surface mixed layer. This corresponds to an air–sea CO 2 uptake of 10.6 mmol m −2 sea ice d −1 or to 3.3 ton km −2 ice floe week −1 . This is markedly higher than the estimated primary production within the ice floe of 0.3–1.3 mmol m −2 sea ice d −1 . Finally, the presence of ikaite in sea ice and the dissolution of the mineral during melting of the sea ice and mixing of the melt water into the surface oceanic mixed layer accounted for half of the estimated p CO 2 uptake.

Description: Drifting snow climate of the Greenland ice sheet: a study with a regional climate model The Cryosphere, 6, 891-899, 2012 Author(s): J. T. M. Lenaerts, M. R. van den Broeke, J. H. van Angelen, E. van Meijgaard, and S. J. Déry This paper presents the drifting snow climate of the Greenland ice sheet, using output from a high-resolution (∼11 km) regional climate model. Because reliable direct observations of drifting snow do not exist, we evaluate the modeled near-surface climate instead, using automatic weather station (AWS) observations from the K-transect and find that RACMO2 realistically simulates near-surface wind speed and relative humidity, two variables that are important for drifting snow. Integrated over the ice sheet, drifting snow sublimation (SU ds ) equals 24 ± 3 Gt yr −1 , and is significantly larger than surface sublimation (SU s , 16 ± 2 Gt yr −1 ). SU ds strongly varies between seasons, and is only important in winter, when surface sublimation and runoff are small. A rapid transition exists between the winter season, when snowfall and SU ds are important, and the summer season, when snowmelt is significant, which increases surface snow density and thereby limits drifting snow processes. Drifting snow erosion (ER ds ) is only important on a regional scale. In recent decades, following decreasing wind speed and rising near-surface temperatures, SU ds exhibits a negative trend (0.1 ± 0.1 Gt yr −1 ), which is compensated by an increase in SU s of similar magnitude.

Description: Antarctic sea ice variability and trends, 1979–2010 The Cryosphere, 6, 871-880, 2012 Author(s): C. L. Parkinson and D. J. Cavalieri In sharp contrast to the decreasing sea ice coverage of the Arctic, in the Antarctic the sea ice cover has, on average, expanded since the late 1970s. More specifically, satellite passive-microwave data for the period November 1978–December 2010 reveal an overall positive trend in ice extents of 17 100 ± 2300 km 2 yr −1 . Much of the increase, at 13 700 ± 1500 km 2 yr −1 , has occurred in the region of the Ross Sea, with lesser contributions from the Weddell Sea and Indian Ocean. One region, that of the Bellingshausen/Amundsen Seas, has (like the Arctic) instead experienced significant sea ice decreases, with an overall ice extent trend of −8200 ± 1200 km 2 yr −1 . When examined through the annual cycle over the 32-yr period 1979–2010, the Southern Hemisphere sea ice cover as a whole experienced positive ice extent trends in every month, ranging in magnitude from a low of 9100 ± 6300 km 2 yr −1 in February to a high of 24 700 ± 10 000 km 2 yr −1 in May. The Ross Sea and Indian Ocean also had positive trends in each month, while the Bellingshausen/Amundsen Seas had negative trends in each month, and the Weddell Sea and western Pacific Ocean had a mixture of positive and negative trends. Comparing ice-area results to ice-extent results, in each case the ice-area trend has the same sign as the ice-extent trend, but the magnitudes of the two trends differ, and in some cases these differences allow inferences about the corresponding changes in sea ice concentrations. The strong pattern of decreasing ice coverage in the Bellingshausen/Amundsen Seas region and increasing ice coverage in the Ross Sea region is suggestive of changes in atmospheric circulation. This is a key topic for future research.

Description: Improved modelling of Siberian river flow through the use of an alternative frozen soil hydrology scheme in a land surface model The Cryosphere, 6, 859-870, 2012 Author(s): D. L. Finney, E. Blyth, and R. Ellis A parameterisation to incorporate the effects of frozen soil on modelled hydrology is described and implemented within a land surface model, the Joint UK Land Surface Environment Simulator. It is shown to generally improve the modelled flow of Siberian rivers compared to observations, specifically in seasons of freezing or thawing soil. Most noticeably, the revised model increases the snowmelt flow peak by 26–100% compared to the control model, thereby better matching observed flows. The model physics resulting in the changes to river flow are discussed and attention is given to the effect of inaccuracies in snowfall driving data which can hinder the comparison of new model processes.

Description: Arctic sea ice variability and trends, 1979–2010 The Cryosphere, 6, 881-889, 2012 Author(s): D. J. Cavalieri and C. L. Parkinson Analyses of 32 yr (1979–2010) of Arctic sea ice extents and areas derived from satellite passive microwave radiometers are presented for the Northern Hemisphere as a whole and for nine Arctic regions. There is an overall negative yearly trend of −51.5 ± 4.1 × 10 3 km 2 yr −1 (−4.1 ± 0.3% decade −1 ) in sea ice extent for the hemisphere. The yearly sea ice extent trends for the individual Arctic regions are all negative except for the Bering Sea: −3.9 ± 1.1 × 10 3 km 2 yr −1 (−8.7 ± 2.5% decade −1 ) for the Seas of Okhotsk and Japan, +0.3 ± 0.8 × 10 3 km 2 yr −1 (+1.2 ± 2.7% decade −1 ) for the Bering Sea, −4.4 ± 0.7 × 10 3 km 2 yr −1 (−5.1 ± 0.9% decade −1 ) for Hudson Bay, −7.6 ± 1.6 × 10 3 km 2 yr −1 (−8.5 ± 1.8% decade −1 ) for Baffin Bay/Labrador Sea, −0.5 ± 0.3 × 10 3 km 2 yr −1 (−5.9 ± 3.5% decade −1 ) for the Gulf of St. Lawrence, −6.5 ± 1.1 × 10 3 km 2 yr −1 (−8.6 ± 1.5% decade −1 ) for the Greenland Sea, −13.5 ± 2.3 × 10 3 km 2 yr −1 (−9.2 ± 1.6% decade −1 ) for the Kara and Barents Seas, −14.6 ± 2.3 × 10 3 km 2 yr −1 (−2.1 ± 0.3% decade −1 ) for the Arctic Ocean, and −0.9 ± 0.4 × 10 3 km 2 yr −1 (−1.3 ± 0.5% decade −1 ) for the Canadian Archipelago. Similarly, the yearly trends for sea ice areas are all negative except for the Bering Sea. On a seasonal basis for both sea ice extents and areas, the largest negative trend is observed for summer with the next largest negative trend being for autumn. Both the sea ice extent and area trends vary widely by month depending on region and season. For the Northern Hemisphere as a whole, all 12 months show negative sea ice extent trends with a minimum magnitude in May and a maximum magnitude in September, whereas the corresponding sea ice area trends are smaller in magnitude and reach minimum and maximum values in March and September.

Description: Borehole temperatures reveal details of 20th century warming at Bruce Plateau, Antarctic Peninsula The Cryosphere, 6, 675-686, 2012 Author(s): V. Zagorodnov, O. Nagornov, T. A. Scambos, A. Muto, E. Mosley-Thompson, E. C. Pettit, and S. Tyuflin Two ice core boreholes of 143.18 m and 447.73 m (bedrock) were drilled during the 2009–2010 austral summer on the Bruce Plateau at a location named LARISSA Site Beta (66°02' S, 64°04' W, 1975.5 m a.s.l.). Both boreholes were logged with thermistors shortly after drilling. The shallow borehole was instrumented for 4 months with a series of resistance thermometers with satellite uplink. Surface temperature proxy data derived from an inversion of the borehole temperature profiles are compared to available multi-decadal records from weather stations and ice cores located along a latitudinal transect of the Antarctic Peninsula to West Antarctica. The LARISSA Site Beta profiles show temperatures decreasing from the surface downward through the upper third of the ice, and warming thereafter to the bed. The average temperature for the most recent year is −14.78°C (measured at 15 m depth, abbreviated T 15 ). A minimum temperature of −15.8°C is measured at 173 m depth, and basal temperature is estimated to be −10.2°C. Current mean annual temperature and the gradient in the lower part of the measured temperature profile have a best fit with an accumulation rate of 1.9×10 3 kg m −2 a −1 and basal heat flux ( q ) of 88 mW m −2 , if steady-state conditions are assumed. However, the mid-level temperature variations show that recent temperature has varied significantly. Reconstructed surface temperatures ( T s = T 15 ) over the last 200 yr are derived by an inversion technique (Tikhonov and Samarskii, 1990). From this, we find that cold temperatures (minimum T s =−16.2°C) prevailed from ~1920 to ~1940, followed by a gradual rise of temperature to −14.2°C around 1995, then cooling over the following decade and warming in the last few years. The coldest period was preceded by a relatively warm 19th century at T 15 ≥−15°C. To facilitate regional comparisons of the surface temperature history, we use our T 15 data and nearby weather station records to refine estimates of lapse rates (altitudinal, adjusted for latitude: Γ a(l) ). Good temporal and spatial consistency of Γ a(l) over the last 35 yr are observed, implying that the climate trends observed here are regional and consistent over a broad altitude range.

Description: Statistical adaptation of ALADIN RCM outputs over the French Alps – application to future climate and snow cover The Cryosphere, 6, 785-805, 2012 Author(s): M. Rousselot, Y. Durand, G. Giraud, L. Mérindol, I. Dombrowski-Etchevers, M. Déqué, and H. Castebrunet In this study, snowpack scenarios are modelled across the French Alps using dynamically downscaled variables from the ALADIN Regional Climate Model (RCM) for the control period (1961–1990) and three emission scenarios (SRES B1, A1B and A2) for the mid- and late 21st century (2021–2050 and 2071–2100). These variables are statistically adapted to the different elevations, aspects and slopes of the Alpine massifs. For this purpose, we use a simple analogue criterion with ERA40 series as well as an existing detailed climatology of the French Alps (Durand et al., 2009a) that provides complete meteorological fields from the SAFRAN analysis model. The resulting scenarios of precipitation, temperature, wind, cloudiness, longwave and shortwave radiation, and humidity are used to run the physical snow model CROCUS and simulate snowpack evolution over the massifs studied. The seasonal and regional characteristics of the simulated climate and snow cover changes are explored, as is the influence of the scenarios on these changes. Preliminary results suggest that the snow water equivalent (SWE) of the snowpack will decrease dramatically in the next century, especially in the Southern and Extreme Southern parts of the Alps. This decrease seems to result primarily from a general warming throughout the year, and possibly a deficit of precipitation in the autumn. The magnitude of the snow cover decline follows a marked altitudinal gradient, with the highest altitudes being less exposed to climate change. Scenario A2, with its high concentrations of greenhouse gases, results in a SWE reduction roughly twice as large as in the low-emission scenario B1 by the end of the century. This study needs to be completed using simulations from other RCMs, since a multi-model approach is essential for uncertainty analysis.

Description: Sensitivity of basal conditions in an inverse model: Vestfonna ice cap, Nordaustlandet/Svalbard The Cryosphere, 6, 771-783, 2012 Author(s): M. Schäfer, T. Zwinger, P. Christoffersen, F. Gillet-Chaulet, K. Laakso, R. Pettersson, V. A. Pohjola, T. Strozzi, and J. C. Moore The dynamics of Vestfonna ice cap (Svalbard) are dominated by fast-flowing outlet glaciers. Its mass balance is poorly known and affected dynamically by these fast-flowing outlet glaciers. Hence, it is a challenging target for ice flow modeling. Precise knowledge of the basal conditions and implementation of a good sliding law are crucial for the modeling of this ice cap. Here we use the full-Stokes finite element code Elmer/Ice to model the 3-D flow over the whole ice cap. We use a Robin inverse method to infer the basal friction from the surface velocities observed in 1995. Our results illustrate the importance of the basal friction parameter in reproducing observed velocity fields. We also show the importance of having variable basal friction as given by the inverse method to reproduce the velocity fields of each outlet glacier – a simple parametrization of basal friction cannot give realistic velocities in a forward model. We study the robustness and sensitivity of this method with respect to different parameters (mesh characteristics, ice temperature, errors in topographic and velocity data). The uncertainty in the observational parameters and input data proved to be sufficiently small as not to adversely affect the fidelity of the model.

Description: Glacier dynamics over the last quarter of a century at Helheim, Kangerdlugssuaq and 14 other major Greenland outlet glaciers The Cryosphere, 6, 923-937, 2012 Author(s): S. L. Bevan, A. J. Luckman, and T. Murray The Greenland ice sheet is experiencing increasing rates of mass loss, the majority of which results from changes in discharge from tidewater glaciers. Both atmospheric and ocean drivers have been implicated in these dynamic changes, but understanding the nature of the response has been hampered by the lack of measurements of glacier flow rates predating the recent period of warming. Here, using Landsat-5 data from 1985 onwards, we extend back in time the record of surface velocities and ice-front position for 16 of Greenland's fastest-flowing tidewater glaciers, and compare these to more recent data from Landsat-7 and satellite-borne synthetic-aperture radar. Climate re-analysis data and sea surface temperatures from 1982 show that since 1995 most of Greenland and its surrounding oceans have experienced significant overall warming, and a switch to a warming trend. During the period from 1985 to 1995 when Greenland and the surrounding oceans were not warming, major tidewater outlet glaciers around Greenland, including Kangerdlugssuaq and Helheim, were dynamically stable. Since the mid-1990s, glacier discharge has consistently been both greater and more variable. Together, these observations support the hypothesis that recent dynamic change is a rapid response to climate forcing. Both air and ocean temperatures in this region are predicted to continue to warm, and will therefore likely drive further change in outlet glacier discharge.

Description: A simple inverse method for the distribution of basal sliding coefficients under ice sheets, applied to Antarctica The Cryosphere, 6, 953-971, 2012 Author(s): D. Pollard and R. M. DeConto Variations in intrinsic bed conditions that affect basal sliding, such as the distribution of deformable sediment versus hard bedrock, are important boundary conditions for large-scale ice-sheet models, but are hard to observe and remain largely uncertain below the modern Greenland and Antarctic ice sheets. Here a very simple model-based method is described for deducing the modern spatial distribution of basal sliding coefficients. The model is run forward in time, and the basal sliding coefficient at each grid point is periodically increased or decreased depending on whether the local ice surface elevation is too high or too low compared to observed in areas of unfrozen bed. The method considerably reduces large-scale errors in Antarctic ice elevation, from several 100s to several 10s of meters in most regions. Remaining ice elevation errors over mountain ranges such as the Transantarctics are further improved by parameterizing the possible effect of sub-grid topography in the basal sliding law, representing sliding in deep valleys. Results are compared with modern velocity data, and various sensitivity tests are described in Appendices.

Description: 3-D image-based numerical computations of snow permeability: links to specific surface area, density, and microstructural anisotropy The Cryosphere, 6, 939-951, 2012 Author(s): N. Calonne, C. Geindreau, F. Flin, S. Morin, B. Lesaffre, S. Rolland du Roscoat, and P. Charrier We used three-dimensional (3-D) images of snow microstructure to carry out numerical estimations of the full tensor of the intrinsic permeability of snow ( K ). This study was performed on 35 snow samples, spanning a wide range of seasonal snow types. For several snow samples, a significant anisotropy of permeability was detected and is consistent with that observed for the effective thermal conductivity obtained from the same samples. The anisotropy coefficient, defined as the ratio of the vertical over the horizontal components of K , ranges from 0.74 for a sample of decomposing precipitation particles collected in the field to 1.66 for a depth hoar specimen. Because the permeability is related to a characteristic length, we introduced a dimensionless tensor K *= K / r es 2 , where the equivalent sphere radius of ice grains ( r es ) is computed from the specific surface area of snow (SSA) and the ice density (ρ i ) as follows: r es =3/(SSA×ρ i . We define K and K * as the average of the diagonal components of K and K *, respectively. The 35 values of K * were fitted to snow density (ρ s ) and provide the following regression: K = (3.0 ± 0.3) r es 2 exp((−0.0130 ± 0.0003)ρ s ). We noted that the anisotropy of permeability does not affect significantly the proposed equation. This regression curve was applied to several independent datasets from the literature and compared to other existing regression curves or analytical models. The results show that it is probably the best currently available simple relationship linking the average value of permeability, K , to snow density and specific surface area.

Description: Evaluation of the criticality of cracks in ice shelves using finite element simulations The Cryosphere, 6, 973-984, 2012 Author(s): C. Plate, R. Müller, A. Humbert, and D. Gross The ongoing disintegration of large ice shelf parts in Antarctica raise the need for a better understanding of the physical processes that trigger critical crack growth in ice shelves. Finite elements in combination with configurational forces facilitate the analysis of single surface fractures in ice under various boundary conditions and material parameters. The principles of linear elastic fracture mechanics are applied to show the strong influence of different depth dependent functions for the density and the Young's modulus on the stress intensity factor K I at the crack tip. Ice, for this purpose, is treated as an elastically compressible solid and the consequences of this choice in comparison to the predominant incompressible approaches are discussed. The computed stress intensity factors K I for dry and water filled cracks are compared to critical values K Ic from measurements that can be found in literature.

Description: Limitations of a coupled regional climate model in the reproduction of the observed Arctic sea-ice retreat The Cryosphere, 6, 985-998, 2012 Author(s): W. Dorn, K. Dethloff, and A. Rinke The effects of internal model variability on the simulation of Arctic sea-ice extent and volume have been examined with the aid of a seven-member ensemble with a coupled regional climate model for the period 1948–2008. Beyond general weaknesses related to insufficient representation of feedback processes, it is found that the model's ability to reproduce observed summer sea-ice retreat depends mainly on two factors: the correct simulation of the atmospheric circulation during the summer months and the sea-ice volume at the beginning of the melting period. Since internal model variability shows its maximum during the summer months, the ability to reproduce the observed atmospheric summer circulation is limited. In addition, the atmospheric circulation during summer also significantly affects the sea-ice volume over the years leading to a limited ability to start with reasonable sea-ice volume into the melting period. Furthermore, the sea-ice volume pathway shows notable decadal variability that varies in amplitude among the ensemble members. The scatter is particularly large in periods when the ice volume increases, indicating limited skill in reproducing high-ice years.

Description: Thin-layer effects in glaciological seismic amplitude-versus-angle (AVA) analysis: implications for characterising a subglacial till unit, Russell Glacier, West Greenland The Cryosphere, 6, 909-922, 2012 Author(s): A. D. Booth, R. A. Clark, B. Kulessa, T. Murray, J. Carter, S. Doyle, and A. Hubbard Seismic amplitude-versus-angle (AVA) methods are a powerful means of quantifying the physical properties of subglacial material, but serious interpretative errors can arise when AVA is measured over a thinly-layered substrate. A substrate layer with a thickness less than 1/4 of the seismic wavelength, λ, is considered "thin", and reflections from its bounding interfaces superpose and appear in seismic data as a single reflection event. AVA interpretation of subglacial till can be vulnerable to such thin-layer effects, since a lodged (non-deforming) till can be overlain by a thin (metre-scale) cap of dilatant (deforming) till. We assess the potential for misinterpretation by simulating seismic data for a stratified subglacial till unit, with an upper dilatant layer between 0.1–5.0 m thick (λ / 120 to 〉 λ / 4, with λ = 12 m). For dilatant layers less than λ / 6 thick, conventional AVA analysis yields acoustic impedance and Poisson's ratio that indicate contradictory water saturation. A thin-layer interpretation strategy is proposed, that accurately characterises the model properties of the till unit. The method is applied to example seismic AVA data from Russell Glacier, West Greenland, in which characteristics of thin-layer responses are evident. A subglacial till deposit is interpreted, having lodged till (acoustic impedance = 4.26±0.59 × 10 6 kg m −2 s −1 ) underlying a water-saturated dilatant till layer (thickness 〈 2 m, Poisson's ratio ~ 0.5). Since thin-layer considerations offer a greater degree of complexity in an AVA interpretation, and potentially avoid misinterpretations, they are a valuable aspect of quantitative seismic analysis, particularly for characterising till units.

Description: Laboratory study of initial sea-ice growth: properties of grease ice and nilas The Cryosphere, 6, 729-741, 2012 Author(s): A. K. Naumann, D. Notz, L. Håvik, and A. Sirevaag We investigate initial sea-ice growth in an ice-tank study by freezing an NaCl solution of about 29 g kg −1 in three different setups: grease ice grew in experiments with waves and in experiments with a current and wind, while nilas formed in a quiescent experimental setup. In this paper we focus on the differences in bulk salinity, solid fraction and thickness between these two ice types. The bulk salinity of the grease-ice layer in our experiments remained almost constant until the ice began to consolidate. In contrast, the initial bulk-salinity evolution of the nilas is well described by a linear decrease of about 2.1 g kg −1 h −1 independent of air temperature. This rapid decrease can be qualitatively understood by considering a Rayleigh number that became maximum while the nilas was still less than 1 cm thick. Comparing three different methods to measure solid fraction in grease ice based on (a) salt conservation, (b) mass conservation and (c) energy conservation, we find that the method based on salt conservation does not give reliable results if the salinity of the interstitial water is approximated as being equal to the salinity of the underlying water. Instead the increase in salinity of the interstitial water during grease-ice formation must be taken into account. In our experiments, the solid fraction of grease ice was relatively constant with values of 0.25, whereas it increased to values as high as 0.50 as soon as the grease ice consolidated at its surface. In contrast, the solid fraction of the nilas increased continuously in the first hours of ice formation and reached an average value of 0.55 after 4.5 h. The spatially averaged ice thickness was twice as large in the first 24 h of ice formation in the setup with a current and wind compared to the other two setups, since the wind kept parts of the water surface ice free and therefore allowed for a higher heat loss from the water. The development of the ice thickness can be reproduced well with simple, one dimensional models that only require air temperature or ice surface temperature as input.

Description: Extrapolating glacier mass balance to the mountain-range scale: the European Alps 1900–2100 The Cryosphere, 6, 713-727, 2012 Author(s): M. Huss This study addresses the extrapolation of in-situ glacier mass balance measurements to the mountain-range scale and aims at deriving time series of area-averaged mass balance and ice volume change for all glaciers in the European Alps for the period 1900–2100. Long-term mass balance series for 50 Swiss glaciers based on a combination of field data and modelling, and WGMS data for glaciers in Austria, France and Italy are used. A complete glacier inventory is available for the year 2003. Mass balance extrapolation is performed based on (1) arithmetic averaging, (2) glacier hypsometry, and (3) multiple regression. Given a sufficient number of data series, multiple regression with variables describing glacier geometry performs best in reproducing observed spatial mass balance variability. Future mass changes are calculated by driving a combined model for mass balance and glacier geometry with GCM ensembles based on four emission scenarios. Mean glacier mass balance in the European Alps is −0.31 ± 0.04 m w.e. a −1 in 1900–2011, and −1 m w.e. a −1 over the last decade. Total ice volume change since 1900 is −96 ± 13 km 3 ; annual values vary between −5.9 km 3 (1947) and +3.9 km 3 (1977). Mean mass balances are expected to be around −1.3 m w.e. a −1 by 2050. Model results indicate a glacier area reduction of 4–18% relative to 2003 for the end of the 21st century.

Description: Spatial patterns of North Atlantic Oscillation influence on mass balance variability of European glaciers The Cryosphere, 6, 661-673, 2012 Author(s): B. Marzeion and A. Nesje We present and validate a set of minimal models of glacier mass balance variability. The most skillful model is then applied to reconstruct 7735 individual time series of mass balance variability for all glaciers in the European Alps and Scandinavia. Subsequently, we investigate the influence of atmospheric variability associated with the North Atlantic Oscillation (NAO) on the glaciers' mass balances. We find a spatial coherence in the glaciers' sensitivity to NAO forcing which is caused by regionally similar mechanisms relating the NAO forcing to the mass balance: in southwestern Scandinavia, winter precipitation causes a correlation of mass balances with the NAO. In northern Scandinavia, temperature anomalies outside the core winter season cause an anti-correlation between NAO and mass balances. In the western Alps, both temperature and winter precipitation anomalies lead to a weak anti-correlation of mass balances with the NAO, while in the eastern Alps, the influences of winter precipitation and temperature anomalies tend to cancel each other, and only on the southern side a slight anti-correlation of mass balances with the NAO prevails.

Description: Brief communication "Can recent ice discharges following the Larsen-B ice-shelf collapse be used to infer the driving mechanisms of millennial-scale variations of the Laurentide ice sheet?" The Cryosphere, 6, 687-693, 2012 Author(s): J. Alvarez-Solas, A. Robinson, and C. Ritz The effects of an ice-shelf collapse on inland glacier dynamics have recently been widely studied, especially since the breakup of the Antarctic Peninsula's Larsen-B ice shelf in 2002. Several studies have documented acceleration of the ice streams that were flowing into the former ice shelf. The mechanism responsible for such a speed-up lies with the removal of the ice-shelf backforce. Independently, it is also well documented that during the last glacial period, the Northern Hemisphere ice sheets experienced large discharges into the ocean, likely reflecting ice flow acceleration episodes on the millennial time scale. The classic interpretation of the latter is based on the existence of an internal thermo-mechanical feedback with the potential to generate oscillatory behavior in the ice sheets. Here we would like to widen the debate by considering that Larsen-B-like glacial analog episodes could have contributed significantly to the registered millennial-scale variablity.

Description: Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR The Cryosphere, 6, 695-711, 2012 Author(s): B. Franco, X. Fettweis, C. Lang, and M. Erpicum With the aim to force an ice dynamical model, the Greenland ice sheet (GrIS) surface mass balance (SMB) was modelled at different spatial resolutions (15–50 km) for the period 1990–2010, using the regional climate model MAR (Modèle Atmosphérique Régional) forced by the ERA-INTERIM reanalysis. This comparison revealed that (i) the inter-annual variability of the SMB components is consistent within the different spatial resolutions investigated, (ii) the MAR model simulates heavier precipitation on average over the GrIS with decreasing spatial resolution, and (iii) the SMB components (except precipitation) can be derived from a simulation at lower resolution with an "intelligent" interpolation. This interpolation can also be used to approximate the SMB components over another topography/ice sheet mask of the GrIS. These results are important for the forcing of an ice dynamical model needed to enable future projections of the GrIS contribution to sea level rise over the coming centuries.

Description: Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP The Cryosphere, 6, 573-588, 2012 Author(s): F. Pattyn, C. Schoof, L. Perichon, R. C. A. Hindmarsh, E. Bueler, B. de Fleurian, G. Durand, O. Gagliardini, R. Gladstone, D. Goldberg, G. H. Gudmundsson, P. Huybrechts, V. Lee, F. M. Nick, A. J. Payne, D. Pollard, O. Rybak, F. Saito, and A. Vieli Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.

Description: An assessment of key model parametric uncertainties in projections of Greenland Ice Sheet behavior The Cryosphere, 6, 589-606, 2012 Author(s): P. J. Applegate, N. Kirchner, E. J. Stone, K. Keller, and R. Greve Lack of knowledge about the values of ice sheet model input parameters introduces substantial uncertainty into projections of Greenland Ice Sheet contributions to future sea level rise. Computer models of ice sheet behavior provide one of several means of estimating future sea level rise due to mass loss from ice sheets. Such models have many input parameters whose values are not well known. Recent studies have investigated the effects of these parameters on model output, but the range of potential future sea level increases due to model parametric uncertainty has not been characterized. Here, we demonstrate that this range is large, using a 100-member perturbed-physics ensemble with the SICOPOLIS ice sheet model. Each model run is spun up over 125 000 yr using geological forcings and subsequently driven into the future using an asymptotically increasing air temperature anomaly curve. All modeled ice sheets lose mass after 2005 AD. Parameters controlling surface melt dominate the model response to temperature change. After culling the ensemble to include only members that give reasonable ice volumes in 2005 AD, the range of projected sea level rise values in 2100 AD is ~40 % or more of the median. Data on past ice sheet behavior can help reduce this uncertainty, but none of our ensemble members produces a reasonable ice volume change during the mid-Holocene, relative to the present. This problem suggests that the model's exponential relation between temperature and precipitation does not hold during the Holocene, or that the central-Greenland temperature forcing curve used to drive the model is not representative of conditions around the ice margin at this time (among other possibilities). Our simulations also lack certain observed physical processes that may tend to enhance the real ice sheet's response. Regardless, this work has implications for other studies that use ice sheet models to project or hindcast the behavior of the Greenland Ice Sheet.

Description: Multi-decadal marine- and land-terminating glacier recession in the Ammassalik region, southeast Greenland The Cryosphere, 6, 625-639, 2012 Author(s): S. H. Mernild, J. K. Malmros, J. C. Yde, and N. T. Knudsen Landsat imagery was applied to elucidate glacier fluctuations of land- and marine-terminating outlet glaciers from the Greenland Ice Sheet (GrIS) and local land-terminating glaciers and ice caps (GIC) peripheral to the GrIS in the Ammassalik region, Southeast Greenland, during the period 1972–2011. Data from 21 marine-terminating glaciers (including the glaciers Helheim, Midgaard, and Fenris), the GrIS land-terminating margin, and 35 GIC were examined and compared to observed atmospheric air temperatures, precipitation, and reconstructed ocean water temperatures (at 400 m depth in the Irminger Sea). Here, we document that net glacier recession has occurred since 1972 in the Ammassalik region for all glacier types and sizes, except for three GIC. The land-terminating GrIS and GIC reflect lower marginal and areal changes than the marine-terminating outlet glaciers. The mean annual land-terminating GrIS and GIC margin recessions were about three to five times lower than the GrIS marine-terminating recession. The marine-terminating outlet glaciers had an average net frontal retreat for 1999–2011 of 0.098 km yr −1 , which was significantly higher than in previous sub-periods 1972–1986 and 1986–1999. For the marine-terminating GrIS, the annual areal recession rate has been decreasing since 1972, while increasing for the land-terminating GrIS since 1986. On average for all the observed GIC, a mean net frontal retreat for 1986–2011 of 0.010 ± 0.006 km yr −1 and a mean areal recession of around 1% per year occurred; overall for all observed GIC, a mean recession rate of 27 ± 24% occurred based on the 1986 GIC area. Since 1986, five GIC melted away in the Ammassalik area.

Description: Simulating melt, runoff and refreezing on Nordenskiöldbreen, Svalbard, using a coupled snow and energy balance model The Cryosphere, 6, 641-659, 2012 Author(s): W. J. J. van Pelt, J. Oerlemans, C. H. Reijmer, V. A. Pohjola, R. Pettersson, and J. H. van Angelen A distributed energy balance model is coupled to a multi-layer snow model in order to study the mass balance evolution and the impact of refreezing on the mass budget of Nordenskiöldbreen, Svalbard. The model is forced with output from the regional climate model RACMO and meteorological data from Svalbard Airport. Extensive calibration and initialisation are performed to increase the model accuracy. For the period 1989–2010, we find a mean net mass balance of −0.39 m w.e. a −1 . Refreezing contributes on average 0.27 m w.e. a −1 to the mass budget and is most pronounced in the accumulation zone. The simulated mass balance, radiative fluxes and subsurface profiles are validated against observations and are generally in good agreement. Climate sensitivity experiments reveal a non-linear, seasonally dependent response of the mass balance, refreezing and runoff to changes in temperature and precipitation. It is shown that including seasonality in climate change, with less pronounced summer warming, reduces the sensitivity of the mass balance and equilibrium line altitude (ELA) estimates in a future climate. The amount of refreezing is shown to be rather insensitive to changes in climate.