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: 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: 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: 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: 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: Interferometric swath processing of Cryosat data for glacial ice topography The Cryosphere, 7, 1857-1867, 2013 Author(s): L. Gray, D. Burgess, L. Copland, R. Cullen, N. Galin, R. Hawley, and V. Helm We have derived digital elevation models (DEMs) over the western part of the Devon Ice Cap in Nunavut, Canada, using "swath processing" of interferometric data collected by Cryosat between February 2011 and January 2012. With the standard ESA (European Space Agency) SARIn (synthetic aperture radar interferometry) level 2 (L2) data product, the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time-delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~ 0.5 to ~ 2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna so that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped, allowing the creation of DEMs from a sequence of L1b SARIn Cryosat data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order 1 m with a standard deviation of about 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.

Description: Supercooled interfacial water in fine-grained soils probed by dielectric spectroscopy The Cryosphere, 7, 1839-1855, 2013 Author(s): A. Lorek and N. Wagner Water substantially affects nearly all physical, chemical and biological processes on the Earth. Recent Mars observations as well as laboratory investigations suggest that water is a key factor of current physical and chemical processes on the Martian surface, e.g. rheological phenomena. Therefore it is of particular interest to get information about the liquid-like state of water on Martian analogue soils for temperatures below 0 °C. To this end, a parallel plate capacitor has been developed to obtain isothermal dielectric spectra of fine-grained soils in the frequency range from 10 Hz to 1.1 MHz at Martian-like temperatures down to −70 °C. Two Martian analogue soils have been investigated: a Ca-bentonite (specific surface of 237 m 2 g −1 , up to 9.4% w / w gravimetric water content) and JSC Mars 1, a volcanic ash (specific surface of 146 m 2 g −1 , up to 7.4% w / w ). Three soil-specific relaxation processes are observed in the investigated frequency–temperature range: two weak high-frequency processes (bound or hydrated water as well as ice) and a strong low-frequency process due to counter-ion relaxation and the Maxwell–Wagner effect. To characterize the dielectric relaxation behaviour, a generalized fractional dielectric relaxation model was applied assuming three active relaxation processes with relaxation time of the i th process modelled with an Eyring equation. The real part of effective complex soil permittivity at 350 kHz was used to determine ice and liquid-like water content by means of the Birchak or CRIM equation. There are evidence that bentonite down to −70 °C has a liquid-like water content of 1.17 monolayers and JSC Mars 1 a liquid-like water content of 1.96 monolayers.

Description: The sensitivity of flowline models of tidewater glaciers to parameter uncertainty The Cryosphere, 7, 1579-1590, 2013 Author(s): E. M. Enderlin, I. M. Howat, and A. Vieli Depth-integrated (1-D) flowline models have been widely used to simulate fast-flowing tidewater glaciers and predict change because the continuous grounding line tracking, high horizontal resolution, and physically based calving criterion that are essential to realistic modeling of tidewater glaciers can easily be incorporated into the models while maintaining high computational efficiency. As with all models, the values for parameters describing ice rheology and basal friction must be assumed and/or tuned based on observations. For prognostic studies, these parameters are typically tuned so that the glacier matches observed thickness and speeds at an initial state, to which a perturbation is applied. While it is well know that ice flow models are sensitive to these parameters, the sensitivity of tidewater glacier models has not been systematically investigated. Here we investigate the sensitivity of such flowline models of outlet glacier dynamics to uncertainty in three key parameters that influence a glacier's resistive stress components. We find that, within typical observational uncertainty, similar initial (i.e., steady-state) glacier configurations can be produced with substantially different combinations of parameter values, leading to differing transient responses after a perturbation is applied. In cases where the glacier is initially grounded near flotation across a basal over-deepening, as typically observed for rapidly changing glaciers, these differences can be dramatic owing to the threshold of stability imposed by the flotation criterion. The simulated transient response is particularly sensitive to the parameterization of ice rheology: differences in ice temperature of ~ 2 °C can determine whether the glaciers thin to flotation and retreat unstably or remain grounded on a marine shoal. Due to the highly non-linear dependence of tidewater glaciers on model parameters, we recommend that their predictions are accompanied by sensitivity tests that take parameter uncertainty into account.

Description: A particle based simulation model for glacier dynamics The Cryosphere, 7, 1591-1602, 2013 Author(s): J. A. Åström, T. I. Riikilä, T. Tallinen, T. Zwinger, D. Benn, J. C. Moore, and J. Timonen A particle-based computer simulation model was developed for investigating the dynamics of glaciers. In the model, large ice bodies are made of discrete elastic particles which are bound together by massless elastic beams. These beams can break, which induces brittle behaviour. At loads below fracture, beams may also break and reform with small probabilities to incorporate slowly deforming viscous behaviour in the model. This model has the advantage that it can simulate important physical processes such as ice calving and fracturing in a more realistic way than traditional continuum models. For benchmarking purposes the deformation of an ice block on a slip-free surface was compared to that of a similar block simulated with a Finite Element full-Stokes continuum model. Two simulations were performed: (1) calving of an ice block partially supported in water, similar to a grounded marine glacier terminus, and (2) fracturing of an ice block on an inclined plane of varying basal friction, which could represent transition to fast flow or surging. Despite several approximations, including restriction to two-dimensions and simplified water-ice interaction, the model was able to reproduce the size distributions of the debris observed in calving, which may be approximated by universal scaling laws. On a moderate slope, a large ice block was stable and quiescent as long as there was enough of friction against the substrate. For a critical length of frictional contact, global sliding began, and the model block disintegrated in a manner suggestive of a surging glacier. In this case the fragment size distribution produced was typical of a grinding process.

Description: Pine Island glacier ice shelf melt distributed at kilometre scales The Cryosphere, 7, 1543-1555, 2013 Author(s): P. Dutrieux, D. G. Vaughan, H. F. J. Corr, A. Jenkins, P. R. Holland, I. Joughin, and A. H. Fleming By thinning and accelerating, West Antarctic ice streams are contributing about 10% of the observed global sea level rise. Much of this ice loss is from Pine Island Glacier, which has thinned since at least 1992, driven by changes in ocean heat transport beneath its ice shelf and retreat of the grounding line. Details of the processes driving this change, however, remain largely elusive, hampering our ability to predict the future behaviour of this and similar systems. Here, a Lagrangian methodology is developed to measure oceanic melting of such rapidly advecting ice. High-resolution satellite and airborne observations of ice surface velocity and elevation are used to quantify patterns of basal melt under the Pine Island Glacier ice shelf and the associated adjustments to ice flow. At the broad scale, melt rates of up to 100 m yr −1 occur near the grounding line, reducing to 30 m yr −1 just 20 km downstream. Between 2008 and 2011, basal melting was largely compensated by ice advection, allowing us to estimate an average loss of ice to the ocean of 87 km 3 yr −1 , in close agreement with 2009 oceanographically constrained estimates. At smaller scales, a network of basal channels typically 500 m to 3 km wide is sculpted by concentrated melt, with kilometre-scale anomalies reaching 50% of the broad-scale basal melt. Basal melting enlarges the channels close to the grounding line, but farther downstream melting tends to diminish them. Kilometre-scale variations in melt are a key component of the complex ice–ocean interaction beneath the ice shelf, implying that greater understanding of their effect, or very high resolution models, are required to predict the sea-level contribution of the region.

Description: A ten-year record of supraglacial lake evolution and rapid drainage in West Greenland using an automated processing algorithm for multispectral imagery The Cryosphere, 7, 1869-1877, 2013 Author(s): B. F. Morriss, R. L. Hawley, J. W. Chipman, L. C. Andrews, G. A. Catania, M. J. Hoffman, M. P. Lüthi, and T. A. Neumann The rapid drainage of supraglacial lakes introduces large pulses of meltwater to the subglacial environment and creates moulins, surface-to-bed conduits for future melt. Introduction of water to the subglacial system has been shown to affect ice flow, and modeling suggests that variability in water supply and delivery to the subsurface play an important role in the development of the subglacial hydrologic system and its ability to enhance or mitigate ice flow. We developed a fully automated method for tracking meltwater and rapid drainages in large (〉 0.125 km 2 ) perennial lakes and applied it to a 10 yr time series of ETM+ and MODIS imagery of an outlet glacier flow band in West Greenland. Results indicate interannual variability in maximum coverage and spatial evolution of total lake area. We identify 238 rapid drainage events, occurring most often at low ( 〈 900 m) and middle (900–1200 m) elevations during periods of net filling or peak lake coverage. We observe a general progression of both lake filling and draining from lower to higher elevations but note that the timing of filling onset, peak coverage, and dissipation are also variable. Lake coverage is sensitive to air temperature, and warm years exhibit greater variability in both coverage evolution and rapid drainage. Mid-elevation drainages in 2011 coincide with large surface velocity increases at nearby GPS sites, though the relationships between ice-shed-scale dynamics and meltwater input are still unclear.

Description: Influence of grain shape on light penetration in snow The Cryosphere, 7, 1803-1818, 2013 Author(s): Q. Libois, G. Picard, J. L. France, L. Arnaud, M. Dumont, C. M. Carmagnola, and M. D. King The energy budget and the photochemistry of a snowpack depend greatly on the penetration of solar radiation in snow. Below the snow surface, spectral irradiance decreases exponentially with depth with a decay constant called the asymptotic flux extinction coefficient. As with the albedo of the snowpack, the asymptotic flux extinction coefficient depends on snow grain shape. While representing snow by a collection of spherical particles has been successful in the numerical computation of albedo, such a description poorly explains the decrease of irradiance in snow with depth. Here we explore the limits of the spherical representation. Under the assumption of geometric optics and weak absorption by snow, the grain shape can be simply described by two parameters: the absorption enhancement parameter B and the geometric asymmetry factor g G . Theoretical calculations show that the albedo depends on the ratio B /(1- g G ) and the asymptotic flux extinction coefficient depends on the product B (1- g G ). To understand the influence of grain shape, the values of B and g G are calculated for a variety of simple geometric shapes using ray tracing simulations. The results show that B and (1- g G ) generally covary so that the asymptotic flux extinction coefficient exhibits larger sensitivity to the grain shape than albedo. In particular it is found that spherical grains propagate light deeper than any other investigated shape. In a second step, we developed a method to estimate B from optical measurements in snow. A multi-layer, two-stream, radiative transfer model, with explicit grain shape dependence, is used to retrieve values of the B parameter of snow by comparing the model to joint measurements of reflectance and irradiance profiles. Such measurements were performed in Antarctica and in the Alps yielding estimates of B between 0.8 and 2.0. In addition, values of B were estimated from various measurements found in the literature, leading to a wider range of values (1.0–9.9) which may be partially explained by the limited accuracy of the data. This work highlights the large variety of snow microstructure and experimentally demonstrates that spherical grains, with B = 1.25, are inappropriate to model irradiance profiles in snow, an important result that should be considered in further studies dedicated to subsurface absorption of short-wave radiation and snow photochemistry.

Description: Influence of ice-sheet geometry and supraglacial lakes on seasonal ice-flow variability The Cryosphere, 7, 1185-1192, 2013 Author(s): I. Joughin, S. B. Das, G. E. Flowers, M. D. Behn, R. B. Alley, M. A. King, B. E. Smith, J. L. Bamber, M. R. van den Broeke, and J. H. van Angelen Supraglacial lakes play an important role in establishing hydrological connections that allow lubricating seasonal meltwater to reach the base of the Greenland Ice Sheet. Here we use new surface velocity observations to examine the influence of supraglacial lake drainages and surface melt rate on ice flow. We find large, spatially extensive speedups concurrent with times of lake drainage, showing that lakes play a key role in modulating regional ice flow. While surface meltwater is supplied to the bed via a geographically sparse network of moulins, the observed ice-flow enhancement suggests that this meltwater spreads widely over the ice-sheet bed. We also find that the complex spatial pattern of speedup is strongly determined by the combined influence of bed and surface topography on subglacial water flow. Thus, modeling of ice-sheet basal hydrology likely will require knowledge of bed topography resolved at scales (sub-kilometer) far finer than existing data (several km).

Description: The effects of additional black carbon on the albedo of Arctic sea ice: variation with sea ice type and snow cover The Cryosphere, 7, 1193-1204, 2013 Author(s): A. A. Marks and M. D. King The response of the albedo of bare sea ice and snow-covered sea ice to the addition of black carbon is calculated. Visible light absorption and light-scattering cross-sections are derived for a typical first-year and multi-year sea ice with both "dry" and "wet" snow types. The cross-sections are derived using data from a 1970s field study that recorded both reflectivity and light penetration in Arctic sea ice and snow overlying sea ice. The variation of absorption cross-section over the visible wavelengths suggests black carbon is the dominating light-absorbing impurity. The response of first-year and multi-year sea ice albedo to increasing black carbon, from 1 to 1024 ng g −1 , in a top 5 cm layer of a 155 cm-thick sea ice was calculated using a radiative-transfer model. The albedo of the first-year sea ice is more sensitive to additional loadings of black carbon than the multi-year sea ice. An addition of 8 ng g −1 of black carbon causes a decrease to 98.7% of the original albedo for first-year sea ice compared to a decrease to 99.7% for the albedo of multi-year sea ice, at a wavelength of 500 nm. The albedo of sea ice is surprisingly unresponsive to additional black carbon up to 100 ng g −1 . Snow layers on sea ice may mitigate the effects of black carbon in sea ice. Wet and dry snow layers of 0.5, 1, 2, 5 and 10 cm depth were added onto the sea ice surface. The albedo of the snow surface was calculated whilst the black carbon in the underlying sea ice was increased. A layer of snow 0.5 cm thick greatly diminishes the effect of black carbon in sea ice on the surface albedo. The albedo of a 2–5 cm snow layer (less than the e -folding depth of snow) is still influenced by the underlying sea ice, but the effect of additional black carbon in the sea ice is masked.

Description: Snow spectral albedo at Summit, Greenland: measurements and numerical simulations based on physical and chemical properties of the snowpack The Cryosphere, 7, 1139-1160, 2013 Author(s): C. M. Carmagnola, F. Domine, M. Dumont, P. Wright, B. Strellis, M. Bergin, J. Dibb, G. Picard, Q. Libois, L. Arnaud, and S. Morin The broadband albedo of surface snow is determined both by the near-surface profile of the physical and chemical properties of the snowpack and by the spectral and angular characteristics of the incident solar radiation. Simultaneous measurements of the physical and chemical properties of snow were carried out at Summit Camp, Greenland (72°36´ N, 38°25´ W, 3210 m a.s.l.) in May and June 2011, along with spectral albedo measurements. One of the main objectives of the field campaign was to test our ability to predict snow spectral albedo by comparing the measured albedo to the albedo calculated with a radiative transfer model, using measured snow physical and chemical properties. To achieve this goal, we made daily measurements of the snow spectral albedo in the range 350–2200 nm and recorded snow stratigraphic information down to roughly 80 cm. The snow specific surface area (SSA) was measured using the DUFISSS instrument (DUal Frequency Integrating Sphere for Snow SSA measurement, Gallet et al., 2009). Samples were also collected for chemical analyses including black carbon (BC) and dust, to evaluate the impact of light absorbing particulate matter in snow. This is one of the most comprehensive albedo-related data sets combining chemical analysis, snow physical properties and spectral albedo measurements obtained in a polar environment. The surface albedo was calculated from density, SSA, BC and dust profiles using the DISORT model (DIScrete Ordinate Radiative Transfer, Stamnes et al., 1988) and compared to the measured values. Results indicate that the energy absorbed by the snowpack through the whole spectrum considered can be inferred within 1.10%. This accuracy is only slightly better than that which can be obtained considering pure snow, meaning that the impact of impurities on the snow albedo is small at Summit. In the near infrared, minor deviations in albedo up to 0.014 can be due to the accuracy of radiation and SSA measurements and to the surface roughness, whereas deviations up to 0.05 can be explained by the spatial heterogeneity of the snowpack at small scales, the assumption of spherical snow grains made for DISORT simulations and the vertical resolution of measurements of surface layer physical properties. At 1430 and around 1800 nm the discrepancies are larger and independent of the snow properties; we propose that they are due to errors in the ice refractive index at these wavelengths. This work contributes to the development of physically based albedo schemes in detailed snowpack models, and to the improvement of retrieval algorithms for estimating snow properties from remote sensing data.

Description: A regional climate model hindcast for Siberia: analysis of snow water equivalent The Cryosphere, 7, 1017-1034, 2013 Author(s): K. Klehmet, B. Geyer, and B. Rockel This study analyzes the added value of a regional climate model hindcast with respect to snow water equivalent (SWE) for Siberia when compared to SWE estimates from forcing NCEP-R. In addition, we examine the discrepancies of simulated SWE to several recent reanalysis products (NCEP-R2 NCEP-CFSR, ERA-Interim). We apply the regional climate model COSMO-CLM (CCLM) to a 50 km grid spacing using NCEP-R1 as driving force to obtain a 63 yr (1948 to 2010) gridded dataset of historical SWE. Simulated regional climate data is necessary because of the absence of station data in that region. To perform large-scale assessments we use the satellite-derived daily SWE product of ESA DUE GlobSnow from 1987 to 2010. Russian station SWE data is used for cross-checking the findings. In January (mid-winter), the SWE hindcast is in good agreement with GlobSnow, whereas it overestimates SWE during the melting season. CCLM shows a clear added value in providing realistic SWE information compared to the driving reanalysis. The temporal consistency of CCLM is higher than that presented by ERA-Interim and NCEP-R2.

Description: Sea ice thickness, freeboard, and snow depth products from Operation IceBridge airborne data The Cryosphere, 7, 1035-1056, 2013 Author(s): N. T. Kurtz, S. L. Farrell, M. Studinger, N. Galin, J. P. Harbeck, R. Lindsay, V. D. Onana, B. Panzer, and J. G. Sonntag The study of sea ice using airborne remote sensing platforms provides unique capabilities to measure a wide variety of sea ice properties. These measurements are useful for a variety of topics including model evaluation and improvement, assessment of satellite retrievals, and incorporation into climate data records for analysis of interannual variability and long-term trends in sea ice properties. In this paper we describe methods for the retrieval of sea ice thickness, freeboard, and snow depth using data from a multi-sensor suite of instruments on NASA's Operation IceBridge airborne campaign. We assess the consistency of the results through comparison with independent data sets that demonstrate that the IceBridge products are capable of providing a reliable record of snow depth and sea ice thickness. We explore the impact of inter-campaign instrument changes and associated algorithm adaptations as well as the applicability of the adapted algorithms to the ongoing IceBridge mission. The uncertainties associated with the retrieval methods are determined and placed in the context of their impact on the retrieved sea ice thickness. Lastly, we present results for the 2009 and 2010 IceBridge campaigns, which are currently available in product form via the National Snow and Ice Data Center.

Description: Summer sea ice characteristics and morphology in the Pacific Arctic sector as observed during the CHINARE 2010 cruise The Cryosphere, 7, 1057-1072, 2013 Author(s): H. Xie, R. Lei, C. Ke, H. Wang, Z. Li, J. Zhao, and S. F. Ackley In the summer of 2010, atmosphere–ice–ocean interaction was studied aboard the icebreaker R/V Xuelong during the Chinese National Arctic Research Expedition (CHINARE), in the sea ice zone of the Pacific Arctic sector between 150° W and 180° W up to 88.5° N. The expedition lasted from 21 July to 28 August and comprised of ice observations and measurements along the cruise track, 8 short-term stations and one 12-day drift station. Ship-based observations of ice thickness and concentration are compared with ice thickness measured by an electromagnetic induction device (EM31) mounted off the ship's side and ice concentrations obtained from AMSR-E. It is found that the modal thickness from ship-based visual observations matches well with the modal thickness from the mounted EM31. A grid of 8 profiles of ice thickness measurements (four repeats) was conducted at the 12-day drift station in the central Arctic (~ 86°50´ N–87°20´ N) and an average melt rate of 2 cm day −1 , primarily bottom melt, was found. As compared with the 2005 data from the Healy/Oden Trans-Arctic Expedition (HOTRAX) for the same sector but ~ 20 days later (9 August to 10 September), the summer 2010 was first-year ice dominant (vs. the multi-year ice dominant in 2005), 70% or less in mean ice concentration (vs. 90% in 2005), and 94–114 cm in mean ice thickness (vs. 150 cm in 2005). Those changes suggest the continuation of ice thinning, less concentration, and younger ice for the summer sea ice in the sector since 2007 when a record minimum sea ice extent was observed. Overall, the measurements provide a valuable dataset of sea ice morphological properties over the Arctic Pacific Sector in summer 2010 and can be used as a benchmark for measurements of future changes.

Description: Hindcasting to measure ice sheet model sensitivity to initial states The Cryosphere, 7, 1083-1093, 2013 Author(s): A. Aschwanden, G. Adalgeirsdóttir, and C. Khroulev Validation is a critical component of model development, yet notoriously challenging in ice sheet modeling. Here we evaluate how an ice sheet system model responds to a given forcing. We show that hindcasting, i.e. forcing a model with known or closely estimated inputs for past events to see how well the output matches observations, is a viable method of assessing model performance. By simulating the recent past of Greenland, and comparing to observations of ice thickness, ice discharge, surface speeds, mass loss and surface elevation changes for validation, we find that the short term model response is strongly influenced by the initial state. We show that the thermal and dynamical states (i.e. the distribution of internal energy and momentum) can be misrepresented despite a good agreement with some observations, stressing the importance of using multiple observations. In particular we identify rates of change of spatially dense observations as preferred validation metrics. Hindcasting enables a qualitative assessment of model performance relative to observed rates of change. It thereby reduces the number of admissible initial states more rigorously than validation efforts that do not take advantage of observed rates of change.

Description: Corrigendum to "A combined approach of remote sensing and airborne electromagnetics to determine the volume of polynya sea ice in the Laptev Sea" published in The Cryosphere, 7, 947−959, 2013 The Cryosphere, 7, 1107-1108, 2013 Author(s): L. Rabenstein, T. Krumpen, S. Hendricks, C. Koeberle, C. Haas, and J. A. Hoelemann No abstract available.

Description: Theoretical study of solar light reflectance from vertical snow surfaces The Cryosphere, 7, 657-666, 2013 Author(s): O. V. Nikolaeva and A. A. Kokhanovsky The influence of horizontal and vertical inhomogeneity of snow surfaces on solar light reflectance is studied using the radiative transfer theory (RTT). We compared 1-D RTT and 2-D RTT and found that large errors are produced if the 1-D RTT is used for the calculation of the snow reflection function (and, therefore, also in the retrievals of the snow grain radii) in 2-D measurement geometries. Such 2-D geometries are common in the procedures for the determination of the effective snow grain radii using near-infrared photography and spectroscopy of vertical snow walls. In particular, we have considered three cases for the numerical calculations: (1) the case with no black film; (2) the case with a black film at the pit's bottom; (3) the case with a black film at the pit's bottom and also at one of the vertical snow walls.

Description: The influence of climate and hydrological variables on opposite anomaly in active-layer thickness between Eurasian and North American watersheds The Cryosphere, 7, 631-645, 2013 Author(s): H. Park, J. Walsh, A. N. Fedorov, A. B. Sherstiukov, Y. Iijima, and T. Ohata This study not only examined the spatiotemporal variations of active-layer thickness (ALT) in permafrost regions during 1948–2006 over the terrestrial Arctic regions experiencing climate changes, but also identified the associated drivers based on observational data and a simulation conducted by a land surface model (CHANGE). The focus on the ALT extends previous studies that have emphasized ground temperatures in permafrost regions. The Ob, Yenisey, Lena, Yukon, and Mackenzie watersheds are foci of the study. Time series of ALT in Eurasian watersheds showed generally increasing trends, while the increase in ALT in North American watersheds was not significant. However, ALT in the North American watersheds has been negatively anomalous since 1990 when the Arctic air temperature entered into a warming phase. The warming temperatures were not simply expressed to increases in ALT. Since 1990 when the warming increased, the forcing of the ALT by the higher annual thawing index (ATI) in the Mackenzie and Yukon basins has been offset by the combined effects of less insulation caused by thinner snow depth and drier soil during summer. In contrast, the increasing ATI together with thicker snow depth and higher summer soil moisture in the Lena contributed to the increase in ALT. The results imply that the soil thermal and moisture regimes formed in the pre-thaw season(s) provide memory that manifests itself during the summer. The different ALT anomalies between Eurasian and North American watersheds highlight increased importance of the variability of hydrological variables.

Description: Estimating the Greenland ice sheet surface mass balance contribution to future sea level rise using the regional atmospheric climate model MAR The Cryosphere, 7, 469-489, 2013 Author(s): X. Fettweis, B. Franco, M. Tedesco, J. H. van Angelen, J. T. M. Lenaerts, M. R. van den Broeke, and H. Gallée To estimate the sea level rise (SLR) originating from changes in surface mass balance (SMB) of the Greenland ice sheet (GrIS), we present 21st century climate projections obtained with the regional climate model MAR (Modèle Atmosphérique Régional), forced by output of three CMIP5 (Coupled Model Intercomparison Project Phase 5) general circulation models (GCMs). Our results indicate that in a warmer climate, mass gain from increased winter snowfall over the GrIS does not compensate mass loss through increased meltwater run-off in summer. Despite the large spread in the projected near-surface warming, all the MAR projections show similar non-linear increase of GrIS surface melt volume because no change is projected in the general atmospheric circulation over Greenland. By coarsely estimating the GrIS SMB changes from GCM output, we show that the uncertainty from the GCM-based forcing represents about half of the projected SMB changes. In 2100, the CMIP5 ensemble mean projects a GrIS SMB decrease equivalent to a mean SLR of +4 ± 2 cm and +9 ± 4 cm for the RCP (Representative Concentration Pathways) 4.5 and RCP 8.5 scenarios respectively. These estimates do not consider the positive melt–elevation feedback, although sensitivity experiments using perturbed ice sheet topographies consistent with the projected SMB changes demonstrate that this is a significant feedback, and highlight the importance of coupling regional climate models to an ice sheet model. Such a coupling will allow the assessment of future response of both surface processes and ice-dynamic changes to rising temperatures, as well as their mutual feedbacks.

Description: New estimates of Arctic and Antarctic sea ice extent during September 1964 from recovered Nimbus I satellite imagery The Cryosphere, 7, 699-705, 2013 Author(s): W. N. Meier, D. Gallaher, and G. G. Campbell Visible satellite imagery from the 1964 Nimbus I satellite has been recovered, digitized, and processed to estimate Arctic and Antarctic sea ice extent for September 1964. September is the month when the Arctic sea ice reaches its minimum annual extent and the Antarctic sea ice reaches its maximum. Images from a three-week period were manually analyzed to estimate the location of the ice edge and then composited to obtain a hemispheric estimate. Uncertainties were based on limitations in the image analysis and the variation of the ice cover over the three-week period. The 1964 Antarctic extent is higher than estimates from the 1979–present passive microwave record, but is in accord with previous indications of higher extents during the 1960s. The Arctic 1964 extent is near the 1979–2000 average from the passive microwave record, suggesting relatively stable summer extents during the 1960s and 1970s preceding the downward trend since 1979 and particularly the large decrease in the last decade. These early satellite data put the recently observed record into a longer-term context.

Description: Ikaite crystal distribution in winter sea ice and implications for CO 2 system dynamics The Cryosphere, 7, 707-718, 2013 Author(s): S. Rysgaard, D. H. Søgaard, M. Cooper, M. Pućko, K. Lennert, T. N. Papakyriakou, F. Wang, N. X. Geilfus, R. N. Glud, J. Ehn, D. F. McGinnis, K. Attard, J. Sievers, J. W. Deming, and D. Barber The precipitation of ikaite (CaCO 3 ⋅ 6H 2 O) in polar sea ice is critical to the efficiency of the sea ice-driven carbon pump and potentially important to the global carbon cycle, yet the spatial and temporal occurrence of ikaite within the ice is poorly known. We report unique observations of ikaite in unmelted ice and vertical profiles of ikaite abundance and concentration in sea ice for the crucial season of winter. Ice was examined from two locations: a 1 m thick land-fast ice site and a 0.3 m thick polynya site, both in the Young Sound area (74° N, 20° W) of NE Greenland. Ikaite crystals, ranging in size from a few μm to 700 μm, were observed to concentrate in the interstices between the ice platelets in both granular and columnar sea ice. In vertical sea ice profiles from both locations, ikaite concentration determined from image analysis, decreased with depth from surface-ice values of 700–900 μmol kg −1 ice (~25 × 10 6 crystals kg −1 ) to values of 100–200 μmol kg −1 ice (1–7 × 10 6 crystals kg −1 ) near the sea ice–water interface, all of which are much higher (4–10 times) than those reported in the few previous studies. Direct measurements of total alkalinity (TA) in surface layers fell within the same range as ikaite concentration, whereas TA concentrations in the lower half of the sea ice were twice as high. This depth-related discrepancy suggests interior ice processes where ikaite crystals form in surface sea ice layers and partly dissolve in layers below. Melting of sea ice and dissolution of observed concentrations of ikaite would result in meltwater with a p CO 2 of

Description: Transient thermal modeling of permafrost conditions in Southern Norway The Cryosphere, 7, 719-739, 2013 Author(s): S. Westermann, T. V. Schuler, K. Gisnås, and B. Etzelmüller Thermal modeling is a powerful tool to infer the temperature regime of the ground in permafrost areas. We present a transient permafrost model, CryoGrid 2, that calculates ground temperatures according to conductive heat transfer in the soil and in the snowpack. CryoGrid 2 is forced by operational air temperature and snow-depth products for potential permafrost areas in Southern Norway for the period 1958 to 2009 at 1 km 2 spatial resolution. In total, an area of about 80 000 km 2 is covered. The model results are validated against borehole temperatures, permafrost probability maps from "bottom temperature of snow" measurements and inventories of landforms indicative of permafrost occurrence. The validation demonstrates that CryoGrid 2 can reproduce the observed lower permafrost limit to within 100 m at all validation sites, while the agreement between simulated and measured borehole temperatures is within 1 K for most sites. The number of grid cells with simulated permafrost does not change significantly between the 1960s and 1990s. In the 2000s, a significant reduction of about 40% of the area with average 2 m ground temperatures below 0 °C is found, which mostly corresponds to degrading permafrost with still negative temperatures in deeper ground layers. The thermal conductivity of the snow is the largest source of uncertainty in CryoGrid 2, strongly affecting the simulated permafrost area. Finally, the prospects of employing CryoGrid 2 as an operational soil-temperature product for Norway are discussed.

Description: An updated and quality controlled surface mass balance dataset for Antarctica The Cryosphere, 7, 583-597, 2013 Author(s): V. Favier, C. Agosta, S. Parouty, G. Durand, G. Delaygue, H. Gallée, A.-S. Drouet, A. Trouvilliez, and G. Krinner We present an updated and quality controlled surface mass balance (SMB) database for the Antarctic ice sheet. Importantly, the database includes formatted metadata, such as measurement technique, elevation, time covered, etc, which allows any user to filter out the data. Here, we discard data with limited spatial and temporal representativeness, too small measurement accuracy, or lack of quality control. Applied to the database, this filtering process gives four times more reliable data than when applied to previously available databases. New data with high spatial resolution are now available over long traverses, and at low elevation in some areas. However, the quality control led to a considerable reduction in the spatial density of data in several regions, particularly over West Antarctica. Over interior plateaus, where the SMB is low, the spatial density of measurements remains high. This quality controlled dataset was compared to results from ERA-Interim reanalysis to assess whether field data allow us to reconstruct an accurate description of the main SMB distribution features in Antarctica. We identified large areas where data gaps impede model validation: except for very few areas (e.g., Adelie Land), measurements in the elevation range between 200 m and 1000 m above sea level are not regularly distributed and do not allow a thorough validation of models in such regions with complex topography, where the highest scattering of SMB values is reported. Clearly, increasing the spatial density of field measurements at low elevations, in the Antarctic Peninsula and in West Antarctica is a scientific priority.

Description: Surface mass balance model intercomparison for the Greenland ice sheet The Cryosphere, 7, 599-614, 2013 Author(s): C. L. Vernon, J. L. Bamber, J. E. Box, M. R. van den Broeke, X. Fettweis, E. Hanna, and P. Huybrechts A number of high resolution reconstructions of the surface mass balance (SMB) of the Greenland ice sheet (GrIS) have been produced using global re-analyses data extending back to 1958. These reconstructions have been used in a variety of applications but little is known about their consistency with each other and the impact of the downscaling method on the result. Here, we compare four reconstructions for the period 1960–2008 to assess the consistency in regional, seasonal and integrated SMB components. Total SMB estimates for the GrIS are in agreement within 34% of the four model average when a common ice sheet mask is used. When models' native land/ice/sea masks are used this spread increases to 57%. Variation in the spread of components of SMB from their mean: runoff 42% (29% native masks), precipitation 20% (24% native masks), melt 38% (74% native masks), refreeze 83% (142% native masks) show, with the exception of refreeze, a similar level of agreement once a common mask is used. Previously noted differences in the models' estimates are partially explained by ice sheet mask differences. Regionally there is less agreement, suggesting spatially compensating errors improve the integrated estimates. Modelled SMB estimates are compared with in situ observations from the accumulation and ablation areas. Agreement is higher in the accumulation area than the ablation area suggesting relatively high uncertainty in the estimation of ablation processes. Since the mid-1990s each model estimates a decreasing annual SMB. A similar period of decreasing SMB is also estimated for the period 1960–1972. The earlier decrease is due to reduced precipitation with runoff remaining unchanged, however, the recent decrease is associated with increased precipitation, now more than compensated for by increased melt driven runoff. Additionally, in three of the four models the equilibrium line altitude has risen since the mid-1990s, reducing the accumulation area at a rate of approximately 60 000 km 2 per decade due to increased melting. Improving process representation requires further study but the use of a single accurate ice sheet mask is a logical way to reduce uncertainty among models.

Description: How does internal variability influence the ability of CMIP5 models to reproduce the recent trend in Southern Ocean sea ice extent? The Cryosphere, 7, 451-468, 2013 Author(s): V. Zunz, H. Goosse, and F. Massonnet Observations over the last 30 yr have shown that the sea ice extent in the Southern Ocean has slightly increased since 1979. Mechanisms responsible for this positive trend have not been well established yet. In this study we tackle two related issues: is the observed positive trend compatible with the internal variability of the system, and do the models agree with what we know about the observed internal variability? For that purpose, we analyse the evolution of sea ice around the Antarctic simulated by 24 different general circulation models involved in the 5th Coupled Model Intercomparison Project (CMIP5), using both historical and hindcast experiments. Our analyses show that CMIP5 models respond to the forcing, including the one induced by stratospheric ozone depletion, by reducing the sea ice cover in the Southern Ocean. Some simulations display an increase in sea ice extent similar to the observed one. According to models, the observed positive trend is compatible with internal variability. However, models strongly overestimate the variance of sea ice extent and the initialization methods currently used in models do not improve systematically the simulated trends in sea ice extent. On the basis of those results, a critical role of the internal variability in the observed increase of sea ice extent in the Southern Ocean could not be ruled out, but current models results appear inadequate to test more precisely this hypothesis.

Description: Intercomparison of retrieval algorithms for the specific surface area of snow from near-infrared satellite data in mountainous terrain, and comparison with the output of a semi-distributed snowpack model The Cryosphere, 7, 741-761, 2013 Author(s): A. Mary, M. Dumont, J.-P. Dedieu, Y. Durand, P. Sirguey, H. Milhem, O. Mestre, H. S. Negi, A. A. Kokhanovsky, M. Lafaysse, and S. Morin This study compares different methods to retrieve the specific surface area (SSA) of snow from satellite radiance measurements in mountainous terrain. It aims at addressing the effect on the retrieval of topographic corrections of reflectance, namely slope and aspect of terrain, multiple reflections on neighbouring slopes and accounting (or not) for the anisotropy of snow reflectance. Using MODerate resolution Imaging Spectrometer (MODIS) data for six different clear sky scenes spanning a wide range of snow conditions during the winter season 2008–2009 over a domain of 46 × 50 km in the French Alps, we compared SSA retrievals with and without topographic correction, with a spherical or non-spherical snow reflectance model and, in spherical case, with or without anisotropy corrections. The retrieved SSA values were compared to field measurements and to the results of the detailed snowpack model Crocus, fed by driving data from the SAFRAN meteorological analysis. It was found that the difference in terms of surface SSA between retrieved values and SAFRAN-Crocus output was minimal when the topographic correction was taken into account, when using a retrieval method assuming disconnected spherical snow grains. In this case, the root mean square deviation was 9.4 m 2 kg −1 and the mean difference was 0.1 m 2 kg −1 , based on 3170 pairs of observation and simulated values. The added-value of the anisotropy correction was not significant in our case, which may be explained by the presence of mixed pixels and surface roughness. MODIS retrieved data show SSA variations with elevation and aspect which are physically consistent and in good agreement with SAFRAN-Crocus outputs. The variability of the MODIS retrieved SSA within the topographic classes of the model was found to be relatively small (3.9 m 2 kg −1 ). This indicates that semi-distributed snowpack simulations in mountainous terrain with a sufficiently large number of classes provides a representation of the snowpack variability consistent with the scale of MODIS 500 m pixels.

Description: Mapping radiation transfer through sea ice using a remotely operated vehicle (ROV) The Cryosphere, 7, 763-777, 2013 Author(s): M. Nicolaus and C. Katlein Transmission of sunlight into and through sea ice is of critical importance for sea-ice associated organisms and photosynthesis because light is their primary energy source. The amount of visible light transferred through sea ice contributes to the energy budget of the sea ice and the uppermost ocean. However, our current knowledge on the amount and distribution of light under sea ice is still restricted to a few local observations, and our understanding of light-driven processes and interdisciplinary interactions is still sparse. The main reasons are that the under-ice environment is difficult to access and that measurements require large logistical and instrumental efforts. Hence, it has not been possible to map light conditions under sea ice over larger areas and to quantify spatial variability on different scales. Here we present a detailed methodological description for operating spectral radiometers on a remotely operated vehicle (ROV) under sea ice. Recent advances in ROV and radiation-sensor technology have allowed us to map under-ice spectral radiance and irradiance on floe scales within a few hours of station time. The ROV was operated directly from the sea ice, allowing for direct relations of optical properties to other sea-ice and surface features. The ROV was flown close to the sea ice in order to capture small-scale variability. Results from the presented data set and similar future studies will allow for better quantification of light conditions under sea ice. The presented experiences will support further developments in order to gather large data sets of under-ice radiation for different ice conditions and during different seasons.

Description: Satellite-derived volume loss rates and glacier speeds for the Cordillera Darwin Icefield, Chile The Cryosphere, 7, 823-839, 2013 Author(s): A. K. Melkonian, M. J. Willis, M. E. Pritchard, A. Rivera, F. Bown, and S. A. Bernstein We produce the first icefield-wide volume change rate and glacier velocity estimates for the Cordillera Darwin Icefield (CDI), a 2605 km 2 temperate icefield in southern Chile (69.6° W, 54.6° S). Velocities are measured from optical and radar imagery between 2001–2011. Thirty-six digital elevation models (DEMs) from ASTER and the SRTM DEM are stacked and a weighted linear regression is applied to elevations on a pixel-by-pixel basis to estimate volume change rates. The CDI lost mass at an average rate of −3.9 ± 1.5 Gt yr −1 between 2000 and 2011, equivalent to a sea level rise (SLR) of 0.01 ± 0.004 mm yr −1 and an area-averaged thinning rate of −1.5 ± 0.6 m w.e.(water equivalent) yr −1 . Thinning is widespread, with concentrations near the front of two northern glaciers (Marinelli, Darwin) and one western (CDI-08) glacier. Thickening is apparent in the south, most notably over the advancing Garibaldi Glacier. The northeastern part of the CDI has an average thinning rate of −1.9 ± 0.7 m w.e. yr −1 , while the southwestern part has an average thinning rate of −1.0 ± 0.4 m w.e. yr −1 . Velocities are obtained over many of the CDI glaciers for the first time. We provide a repeat speed time series at the Marinelli Glacier. There we measure maximum front speeds of 7.5 ± 0.2 m day −1 in 2001, 9.5 ± 0.6 m day −1 in 2003 and 10 ± 0.3 m day −1 in 2011. The maintenance of high front speeds from 2001 to 2011 supports the hypothesis that Marinelli is in the retreat phase of the tidewater cycle, with dynamic thinning governed by the fjord bathymetry.

Description: Snow cover thickness estimation using radial basis function networks The Cryosphere, 7, 841-854, 2013 Author(s): E. Binaghi, V. Pedoia, A. Guidali, and M. Guglielmin This paper reports an experimental study designed for the in-depth investigation of how the radial basis function network (RBFN) estimates snow cover thickness as a function of climate and topographic parameters. The estimation problem is modeled in terms of both function regression and classification, obtaining continuous and discrete thickness values, respectively. The model is based on a minimal set of climatic and topographic data collected from a limited number of stations located in the Italian Central Alps. Several experiments have been conceived and conducted adopting different evaluation indexes. A comparison analysis was also developed for a quantitative evaluation of the advantages of the RBFN method over to conventional widely used spatial interpolation techniques when dealing with critical situations originated by lack of data and limited n -homogeneously distributed instrumented sites. The RBFN model proved competitive behavior and a valuable tool in critical situations in which conventional techniques suffer from a lack of representative data.

Description: Radio-frequency probes of Antarctic ice at South Pole The Cryosphere, 7, 855-866, 2013 Author(s): D. Besson and I. Kravchenko Using hardware developed for the ARA (Askaryan Radio Array) particle astrophysics experiment, we herein report on the amplitude and temporal characteristics of polarized surface radar echo data collected in South Polar ice using radio sounding equipment with 0.5-ns echo-time sampling. We observe strong echoes at 6, 9.6, 13.9, 17, and 19 μs following vertical pulse emission from the surface, corresponding to reflectors in the upper half of the ice sheet. The synchronicity of those echoes for all broadcast azimuthal polarizations affirms the lack of observable birefringence over the upper half of the ice sheet. Of the five strongest echoes, three exhibit an evident amplitude correlation with the local surface ice flow direction, qualitatively consistent with measurements in East Antarctica. Combined with other radio echo sounding data, we conclude that observed birefringent asymmetries at South Pole are generated entirely in the lower half of the ice sheet. By contrast, birefringent asymmetries are observed at shallow depths in East Antarctica.

Description: Density assumptions for converting geodetic glacier volume change to mass change The Cryosphere, 7, 877-887, 2013 Author(s): M. Huss The geodetic method is widely used for assessing changes in the mass balance of mountain glaciers. However, comparison of repeated digital elevation models only provides a glacier volume change that must be converted to a change in mass using a density assumption or model. This study investigates the use of a constant factor for the volume-to-mass conversion based on a firn compaction model applied to simplified glacier geometries with idealized climate forcing, and two glaciers with long-term mass balance series. It is shown that the "density" of geodetic volume change is not a constant factor and is systematically smaller than ice density in most cases. This is explained by the accretion/removal of low-density firn layers, and changes in the firn density profile with positive/negative mass balance. Assuming a value of 850 ± 60 kg m −3 to convert volume change to mass change is appropriate for a wide range of conditions. For short time intervals (≤3 yr), periods with limited volume change, and/or changing mass balance gradients, the conversion factor can however vary from 0–2000 kg m −3 and beyond, which requires caution when interpreting glacier mass changes based on geodetic surveys.

Description: Quantifying present and future glacier melt-water contribution to runoff in a central Himalayan river basin The Cryosphere, 7, 889-904, 2013 Author(s): M. Prasch, W. Mauser, and M. Weber Water supply of most lowland cultures heavily depends on rain and melt water from the upstream mountains. Especially melt-water release of alpine mountain ranges is usually attributed a pivotal role for the water supply of large downstream regions. Water scarcity is assumed as consequence of glacier shrinkage and possible disappearance due to global climate change (GCC), in particular for large parts of Central and Southeast Asia. In this paper, the application and validation of a coupled modeling approach with regional climate model (RCM) outputs and a process-oriented glacier and hydrological model is presented for the central Himalayan Lhasa River basin despite scarce data availability. Current and possible future contributions of ice melt to runoff along the river network are spatially explicitly shown. Its role among the other water balance components is presented. Although glaciers have retreated and will continue to retreat according to the chosen climate scenarios, water availability is and will be primarily determined by monsoon precipitation and snowmelt. Ice melt from glaciers is and will be a minor runoff component in summer monsoon-dominated Himalayan river basins.

Description: Mass balance, runoff and surges of Bering Glacier, Alaska The Cryosphere, 7, 867-875, 2013 Author(s): W. Tangborn The historical net, ablation and accumulation daily balances, as well as runoff of Bering Glacier, Alaska are determined for the 1951–2011 period with the PTAA (precipitation-temperature-area-altitude) model, using daily precipitation and temperature observations collected at the Cordova and Yakutat weather stations, together with the area-altitude distribution of the glacier. The model mean annual balance for this 61 yr period is −0.6 m w.e., the accumulation balance is +1.4 and the ablation balance is −2.0 m w.e. Average annual runoff is 2.5 m w.e. Periodic surges of this glacier transport large volumes of ice to lower elevations where the ablation rate is higher, producing more negative balances and increasing runoff. Runoff from Bering Glacier (derived from simulated ablation and precipitation as rain) is highly correlated with four of the glacier surges that have occurred since 1951. Ice volume loss for the 1972–2003 period measured with the PTAA model is 2.7 km 3 w.e. a −1 and closely agrees with losses for the same period measured with the geodetic method. It is proposed that the timing and magnitude of daily snow accumulation and runoff, both of which are controlled by the glacier's area-altitude distribution and are calculated with the PTAA model, can be used to determine the probability that a glacier will surge.

Description: Corrigendum to "Influence of ablation-related processes in the build-up of simulated Northern Hemisphere ice sheets during the last glacial cycle" published in The Cryosphere, 7, 681–698, 2013 The Cryosphere, 7, 933-934, 2013 Author(s): S. Charbit, C. Dumas, M. Kageyama, D. M. Roche, and C. Ritz No abstract available.

Description: Impacts of mean annual air temperature change on a regional permafrost probability model for the southern Yukon and northern British Columbia, Canada The Cryosphere, 7, 935-946, 2013 Author(s): P. P. Bonnaventure and A. G. Lewkowicz Air temperature changes were applied to a regional model of permafrost probability under equilibrium conditions for an area of nearly 0.5 × 10 6 km 2 in the southern Yukon and northwestern British Columbia, Canada. Associated environmental changes, including snow cover and vegetation, were not considered in the modelling. Permafrost extent increases from 58% of the area (present day: 1971–2000) to 76% under a −1 K cooling scenario, whereas warming scenarios decrease the percentage of permafrost area exponentially to 38% (+ 1 K), 24% (+ 2 K), 17% (+ 3 K), 12% (+ 4 K) and 9% (+ 5 K) of the area. The morphology of permafrost gain/loss under these scenarios is controlled by the surface lapse rate (SLR, i.e. air temperature elevation gradient), which varies across the region below treeline. Areas that are maritime exhibit SLRs characteristically similar above and below treeline resulting in low probabilities of permafrost in valley bottoms. When warming scenarios are applied, a loss front moves to upper elevations (simple unidirectional spatial loss). Areas where SLRs are gently negative below treeline and normal above treeline exhibit a loss front moving up-mountain at different rates according to two separate SLRs (complex unidirectional spatial loss). Areas that display high continentally exhibit bidirectional spatial loss in which the loss front moves up-mountain above treeline and down-mountain below treeline. The parts of the region most affected by changes in MAAT (mean annual air temperature) have SLRs close to 0 K km −1 and extensive discontinuous permafrost, whereas the least sensitive in terms of areal loss are sites above the treeline where permafrost presence is strongly elevation dependent.

Description: Variability of light transmission through Arctic land-fast sea ice during spring The Cryosphere, 7, 977-986, 2013 Author(s): M. Nicolaus, C. Petrich, S. R. Hudson, and M. A. Granskog The amount of solar radiation transmitted through Arctic sea ice is determined by the thickness and physical properties of snow and sea ice. Light transmittance is highly variable in space and time since thickness and physical properties of snow and sea ice are highly heterogeneous on variable time and length scales. We present field measurements of under-ice irradiance along transects under undeformed land-fast sea ice at Barrow, Alaska (March, May, and June 2010). The measurements were performed with a spectral radiometer mounted on a floating under-ice sled. The objective was to quantify the spatial variability of light transmittance through snow and sea ice, and to compare this variability along its seasonal evolution. Along with optical measurements, snow depth, sea ice thickness, and freeboard were recorded, and ice cores were analyzed for chlorophyll a and particulate matter. Our results show that snow cover variability prior to onset of snow melt causes as much relative spatial variability of light transmittance as the contrast of ponded and white ice during summer. Both before and after melt onset, measured transmittances fell in a range from one third to three times the mean value. In addition, we found a twentyfold increase of light transmittance as a result of partial snowmelt, showing the seasonal evolution of transmittance through sea ice far exceeds the spatial variability. However, prior melt onset, light transmittance was time invariant and differences in under-ice irradiance were directly related to the spatial variability of the snow cover.

Description: Snow specific surface area simulation using the one-layer snow model in the Canadian LAnd Surface Scheme (CLASS) The Cryosphere, 7, 961-975, 2013 Author(s): A. Roy, A. Royer, B. Montpetit, P. A. Bartlett, and A. Langlois Snow grain size is a key parameter for modeling microwave snow emission properties and the surface energy balance because of its influence on the snow albedo, thermal conductivity and diffusivity. A model of the specific surface area (SSA) of snow was implemented in the one-layer snow model in the Canadian LAnd Surface Scheme (CLASS) version 3.4. This offline multilayer model (CLASS-SSA) simulates the decrease of SSA based on snow age, snow temperature and the temperature gradient under dry snow conditions, while it considers the liquid water content of the snowpack for wet snow metamorphism. We compare the model with ground-based measurements from several sites (alpine, arctic and subarctic) with different types of snow. The model provides simulated SSA in good agreement with measurements with an overall point-to-point comparison RMSE of 8.0 m 2 kg –1 , and a root mean square error (RMSE) of 5.1 m 2 kg –1 for the snowpack average SSA. The model, however, is limited under wet conditions due to the single-layer nature of the CLASS model, leading to a single liquid water content value for the whole snowpack. The SSA simulations are of great interest for satellite passive microwave brightness temperature assimilations, snow mass balance retrievals and surface energy balance calculations with associated climate feedbacks.

Description: Alpine permafrost thawing during the Medieval Warm Period identified from cryogenic cave carbonates The Cryosphere, 7, 1073-1081, 2013 Author(s): M. Luetscher, M. Borreguero, G. E. Moseley, C. Spötl, and R. L. Edwards Coarse crystalline cryogenic cave carbonates (CCC coarse ) dated to the last glacial period are common in central European caves and provide convincing evidence of palaeo-permafrost during this time. Little is known, however, about the exact nature of the environment in which CCC coarse formed as no modern analogue setting is known. Here, we report the first findings of sub-recent, albeit inactive, CCC coarse from a cave of the Western Alps which is located in the present-day permafrost zone. The globular shape and the presence of ubiquitous euhedral crystal terminations are comparable to previously reported aggregates from the last glacial period and strongly suggest that these aggregates formed subaqueously in pools lacking agitation. Furthermore, stable isotope values of mm-sized spheroids point to calcite precipitation in a closed system with respect to CO 2 , strongly supporting the hypothesis of a cryogenic origin associated with the freezing of water ponds. U-series analyses revealed three clusters of late Holocene calcite precipitation intervals between 2129 and 751 a b2k. These ages correlate with known periods of elevated summer temperatures, suggesting that warming and thawing of the frozen catchment above the cave allowed water infiltration into the karst system. The growth of CCC coarse resulted from the re-freezing of this water in the still cold karst cavities.

Description: A balanced water layer concept for subglacial hydrology in large-scale ice sheet models The Cryosphere, 7, 1095-1106, 2013 Author(s): S. Goeller, M. Thoma, K. Grosfeld, and H. Miller There is currently no doubt about the existence of a widespread hydrological network under the Antarctic Ice Sheet, which lubricates the ice base and thus leads to increased ice velocities. Consequently, ice models should incorporate basal hydrology to obtain meaningful results for future ice dynamics and their contribution to global sea level rise. Here, we introduce the balanced water layer concept, covering two prominent subglacial hydrological features for ice sheet modeling on a continental scale: the evolution of subglacial lakes and balance water fluxes. We couple it to the thermomechanical ice-flow model RIMBAY and apply it to a synthetic model domain. In our experiments we demonstrate the dynamic generation of subglacial lakes and their impact on the velocity field of the overlaying ice sheet, resulting in a negative ice mass balance. Furthermore, we introduce an elementary parametrization of the water flux–basal sliding coupling and reveal the predominance of the ice loss through the resulting ice streams against the stabilizing influence of less hydrologically active areas. We point out that established balance flux schemes quantify these effects only partially as their ability to store subglacial water is lacking.

Description: Modeling energy and mass balance of Shallap Glacier, Peru The Cryosphere, 7, 1787-1802, 2013 Author(s): W. Gurgiser, B. Marzeion, L. Nicholson, M. Ortner, and G. Kaser We calculated the distributed surface mass and energy balance of Shallap Glacier, Cordillera Blanca, Peru (9° S, 77° W, 4700–5700 m a.s.l., ~ 7 km 2 ), on hourly time steps for two years (September 2006–August 2008) using a process-based model and meteorological measurements as input. Model parameter combinations were optimized against 21 temporal readings of 20 stakes in the ablation zone of the glacier. Uncertainty caused by model input parameters and parameterization schemes was estimated using a leave-one out cross-validation scheme, which yields values of root mean square deviation (RMSD) of surface height change 〈 1 m ( 〈 10% of the measured amplitude) for all stakes. With the best parameter combination (smallest RMSD) applied, the modeled annual surface mass balance of the glacier was −0.32 ± 0.4 m w.e. (water equivalent) for September 2006–August 2007 and 0.51 ± 0.56 m w.e. for September 2007–August 2008. While the mass balance above 5000 m was similar in both years (Δ 0.33 ± 0.68 m w.e.) due to similar annual sums of solid precipitation, a difference of 1.97 ± 0.68 m w.e. was calculated for the lower parts of the glacier. This difference is associated with more frequent occurrence of higher snow line altitudes during the first year, which was mainly caused by a higher fraction of liquid precipitation due to higher mean air temperatures. As the net shortwave budget was found to be the main source for ablation throughout the year at Shallap Glacier, lower surface albedo especially caused by lower solid precipitation amounts explains most of the difference in modeled ablation and mass balance between the two years.

Description: Recent extreme light sea ice years in the Canadian Arctic Archipelago: 2011 and 2012 eclipse 1998 and 2007 The Cryosphere, 7, 1753-1768, 2013 Author(s): S. E. L. Howell, T. Wohlleben, A. Komarov, L. Pizzolato, and C. Derksen Remarkably low mean September sea ice area in the Canadian Arctic Archipelago (CAA) was observed in 2011 (146 × 10 3 km 2 ), a record-breaking level that was nearly exceeded in 2012 (150 × 10 3 km 2 ). These values were lower than previous September records set in 1998 (200 × 10 3 km 2 ) and 2007 (220 × 10 3 km 2 ), and are ∼60% lower than the 1981–2010 mean September climatology. In this study, the processes contributing to the extreme light years of 2011 and 2012 were investigated, compared to previous extreme minima of 1998 and 2007, and contrasted against historic summer seasons with above average September ice area. The 2011 minimum was associated with positive June through September (JJAS) surface air temperature (SAT) and net solar radiation ( K * ) anomalies that facilitated rapid melt, coupled with atmospheric circulation that restricted multi-year ice (MYI) inflow from the Arctic Ocean into the CAA. The 2012 minimum was also associated with positive JJAS SAT and K * anomalies with coincident rapid melt, but further ice decline was temporarily mitigated by atmospheric circulation which drove Arctic Ocean MYI inflow into the CAA. Atmospheric circulation was comparable between 2011 and 1998 (impeding Arctic Ocean MYI inflow) and 2012 and 2007 (inducing Arctic Ocean MYI inflow). However, preconditioning was more apparent leading up to 2011 and 2012 than 1998 and 2007. The rapid melt process in 2011 and 2012 was more intense than observed in 1998 and 2007 because of the thinner ice cover being more susceptible to anomalous thermodynamic forcing. The thinner sea ice cover within the CAA in recent years has also helped counteract the processes that facilitate extreme heavy ice years. The recent extreme light years within the CAA are associated with a longer navigation season within the Northwest Passage.

Description: Net accumulation rates derived from ice core stable isotope records of Pío XI glacier, Southern Patagonia Icefield The Cryosphere, 7, 1635-1644, 2013 Author(s): M. Schwikowski, M. Schläppi, P. Santibañez, A. Rivera, and G. Casassa Pío XI, the largest glacier of the Southern Patagonia Icefield, reached its neoglacial maximum extent in 1994 and is one of the few glaciers in that area which is not retreating. In view of the recent warming it is important to understand glacier responses to climate changes. Due to its remoteness and the harsh conditions in Patagonia, no systematic mass balance studies have been performed. In this study we derived net accumulation rates for the period 2000–2006 from a 50 m (33.2 4 m weq) ice core collected in the accumulation area of Pío XI (2600 m a.s.l., 49°16'40"S, 73°21'14"W). Borehole temperatures indicate near temperate ice, but the average melt percent is only 16 ± 14%. Records of stable isotopes are well preserved and were used for identification of annual layers. Net accumulation rates range from 3.4–7.1 water equivalent (m weq) with an average of 5.8 m weq, comparable to precipitation amounts at the Chilean coast, but not as high as expected for the Icefield. Ice core stable isotope data correlate well with upper air temperatures and may be used as temperature proxy.

Description: Seasonal and annual mass balances of Mera and Pokalde glaciers (Nepal Himalaya) since 2007 The Cryosphere, 7, 1769-1786, 2013 Author(s): P. Wagnon, C. Vincent, Y. Arnaud, E. Berthier, E. Vuillermoz, S. Gruber, M. Ménégoz, A. Gilbert, M. Dumont, J. M. Shea, D. Stumm, and B. K. Pokhrel In the Everest region, Nepal, ground-based monitoring programmes were started on the debris-free Mera Glacier (27.7° N, 86.9° E; 5.1 km 2 , 6420 to 4940 m a.s.l.) in 2007 and on the small Pokalde Glacier (27.9° N, 86.8° E; 0.1 km 2 , 5690 to 5430 m a.s.l., ~ 25 km north of Mera Glacier) in 2009. These glaciers lie on the southern flank of the central Himalaya under the direct influence of the Indian monsoon and receive more than 80% of their annual precipitation in summer (June to September). Despite a large inter-annual variability with glacier-wide mass balances ranging from −0.67 ± 0.28 m w.e. in 2011–2012 (Equilibrium-line altitude (ELA) at ~ 5800 m a.s.l.) to +0.46 ± 0.28 m w.e. in 2010–2011 (ELA at ~ 5340 m a.s.l.), Mera Glacier has been shrinking at a moderate mass balance rate of −0.08 ± 0.28 m w.e. yr −1 since 2007. Ice fluxes measured at two distinct transverse cross sections at ~ 5350 m a.s.l. and ~ 5520 m a.s.l. confirm that the mean state of this glacier over the last one or two decades corresponds to a limited mass loss, in agreement with remotely-sensed region-wide mass balances of the Everest area. Seasonal mass balance measurements show that ablation and accumulation are concomitant in summer which in turn is the key season controlling the annual glacier-wide mass balance. Unexpectedly, ablation occurs at all elevations in winter due to wind erosion and sublimation, with remobilised snow potentially being sublimated in the atmosphere. Between 2009 and 2012, the small Pokalde Glacier lost mass more rapidly than Mera Glacier with respective mean glacier-wide mass balances of −0.72 and −0.23 ± 0.28 m w.e. yr −1 . Low-elevation glaciers, such as Pokalde Glacier, have been usually preferred for in-situ observations in Nepal and more generally in the Himalayas, which may explain why compilations of ground-based mass balances are biased toward negative values compared with the regional mean under the present-day climate.

Description: Decadal changes from a multi-temporal glacier inventory of Svalbard The Cryosphere, 7, 1603-1621, 2013 Author(s): C. Nuth, J. Kohler, M. König, A. von Deschwanden, J. O. Hagen, A. Kääb, G. Moholdt, and R. Pettersson We present a multi-temporal digital inventory of Svalbard glaciers with the most recent from the late 2000s containing 33 775 km 2 of glaciers covering 57% of the total land area of the archipelago. At present, 68% of the glacierized area of Svalbard drains through tidewater glaciers that have a total terminus width of ~ 740 km. The glacierized area over the entire archipelago has decreased by an average of 80 km 2 a −1 over the past ~ 30 yr, representing a reduction of 7%. For a sample of ~ 400 glaciers (10 000 km 2 ) in the south and west of Spitsbergen, three digital inventories are available from the 1930/60s, 1990 and 2007 from which we calculate average changes during 2 epochs. In the more recent epoch, the terminus retreat was larger than in the earlier epoch, while area shrinkage was smaller. The contrasting pattern may be explained by the decreased lateral wastage of the glacier tongues. Retreat rates for individual glaciers show a mix of accelerating and decelerating trends, reflecting the large spatial variability of glacier types and climatic/dynamic response times in Svalbard. Lastly, retreat rates estimated by dividing glacier area changes by the tongue width are larger than centerline retreat due to a more encompassing frontal change estimate with inclusion of lateral area loss.

Description: Recent mass balance of the Purogangri Ice Cap, central Tibetan Plateau, by means of differential X-band SAR interferometry The Cryosphere, 7, 1623-1633, 2013 Author(s): N. Neckel, A. Braun, J. Kropáček, and V. Hochschild Due to their remoteness, altitude and harsh climatic conditions, little is known about the glaciological parameters of ice caps on the Tibetan Plateau. This study presents a geodetic mass balance estimate of the Purogangri Ice Cap, Tibet's largest ice field between 2000 and 2012. We utilized data from the actual TerraSAR-X mission and its add-on for digital elevation measurements and compared it with elevation data from the Shuttle Radar Topography Mission. The employed data sets are ideal for this approach as both data sets were acquired at X-band at nearly the same time of the year and are available at a fine grid spacing. In order to derive surface elevation changes we employed two different methods. The first method is based on differential synthetic radar interferometry while the second method uses common DEM differencing. Both approaches revealed a slightly negative mass budget of −44 ± 15 and −38 ± 23 mm w.eq. a −1 (millimeter water equivalent) respectively. A slightly negative trend of −0.15 ± 0.01 km 2 a −1 in glacier extent was found for the same time period employing a time series of Landsat data. Overall, our results show an almost balanced mass budget for the studied time period. Additionally, we detected one continuously advancing glacier tongue in the eastern part of the ice cap.

Description: Current state of glaciers in the tropical Andes: a multi-century perspective on glacier evolution and climate change The Cryosphere, 7, 81-102, 2013 Author(s): A. Rabatel, B. Francou, A. Soruco, J. Gomez, B. Cáceres, J. L. Ceballos, R. Basantes, M. Vuille, J.-E. Sicart, C. Huggel, M. Scheel, Y. Lejeune, Y. Arnaud, M. Collet, T. Condom, G. Consoli, V. Favier, V. Jomelli, R. Galarraga, P. Ginot, L. Maisincho, J. Mendoza, M. Ménégoz, E. Ramirez, P. Ribstein, W. Suarez, M. Villacis, and P. Wagnon The aim of this paper is to provide the community with a comprehensive overview of the studies of glaciers in the tropical Andes conducted in recent decades leading to the current status of the glaciers in the context of climate change. In terms of changes in surface area and length, we show that the glacier retreat in the tropical Andes over the last three decades is unprecedented since the maximum extension of the Little Ice Age (LIA, mid-17th–early 18th century). In terms of changes in mass balance, although there have been some sporadic gains on several glaciers, we show that the trend has been quite negative over the past 50 yr, with a mean mass balance deficit for glaciers in the tropical Andes that is slightly more negative than the one computed on a global scale. A break point in the trend appeared in the late 1970s with mean annual mass balance per year decreasing from −0.2 m w.e. in the period 1964–1975 to −0.76 m w.e. in the period 1976–2010. In addition, even if glaciers are currently retreating everywhere in the tropical Andes, it should be noted that this is much more pronounced on small glaciers at low altitudes that do not have a permanent accumulation zone, and which could disappear in the coming years/decades. Monthly mass balance measurements performed in Bolivia, Ecuador and Colombia show that variability of the surface temperature of the Pacific Ocean is the main factor governing variability of the mass balance at the decadal timescale. Precipitation did not display a significant trend in the tropical Andes in the 20th century, and consequently cannot explain the glacier recession. On the other hand, temperature increased at a significant rate of 0.10 °C decade −1 in the last 70 yr. The higher frequency of El Niño events and changes in its spatial and temporal occurrence since the late 1970s together with a warming troposphere over the tropical Andes may thus explain much of the recent dramatic shrinkage of glaciers in this part of the world.

Description: Calving on tidewater glaciers amplified by submarine frontal melting The Cryosphere, 7, 119-128, 2013 Author(s): M. O'Leary and P. Christoffersen While it has been shown repeatedly that ocean conditions exhibit an important control on the behaviour of grounded tidewater glaciers, modelling studies have focused largely on the effects of basal and surface melting. Here, a finite-element model of stresses near the front of a tidewater glacier is used to investigate the effects of frontal melting on calving, independently of the calving criterion used. Applications of the stress model to idealized scenarios reveal that undercutting of the ice front due to frontal melting can drive calving at up to ten times the mean melt rate. Factors which cause increased frontal melt-driven calving include a strong thermal gradient in the ice, and a concentration of frontal melt at the base of the glacier. These properties are typical of both Arctic and Antarctic tidewater glaciers. The finding that frontal melt near the base is a strong driver of calving leads to the conclusion that water temperatures near the bed of the glacier are critically important to the glacier front, and thus the flow of the glacier. These conclusions are robust against changes in the basal boundary condition and the choice of calving criterion, as well as variations in the glacier size or level of crevassing.

Description: An analysis of present and future seasonal Northern Hemisphere land snow cover simulated by CMIP5 coupled climate models The Cryosphere, 7, 67-80, 2013 Author(s): C. Brutel-Vuilmet, M. Ménégoz, and G. Krinner The 20th century seasonal Northern Hemisphere (NH) land snow cover as simulated by available CMIP5 model output is compared to observations. On average, the models reproduce the observed snow cover extent very well, but the significant trend towards a reduced spring snow cover extent over the 1979–2005 period is underestimated (observed: (−3.4 ± 1.1)% per decade; simulated: (−1.0 ± 0.3)% per decade). We show that this is linked to the simulated Northern Hemisphere extratropical spring land warming trend over the same period, which is also underestimated, although the models, on average, correctly capture the observed global warming trend. There is a good linear correlation between the extent of hemispheric seasonal spring snow cover and boreal large-scale spring surface air temperature in the models, supported by available observations. This relationship also persists in the future and is independent of the particular anthropogenic climate forcing scenario. Similarly, the simulated linear relationship between the hemispheric seasonal spring snow cover extent and global mean annual mean surface air temperature is stable in time. However, the slope of this relationship is underestimated at present (observed: (−11.8 ± 2.7)% °C −1 ; simulated: (−5.1 ± 3.0)% °C −1 ) because the trend towards lower snow cover extent is underestimated, while the recent global warming trend is correctly represented.

Description: An estimate of global glacier volume The Cryosphere, 7, 141-151, 2013 Author(s): A. Grinsted I assess the feasibility of using multivariate scaling relationships to estimate glacier volume from glacier inventory data. Scaling laws are calibrated against volume observations optimized for the specific purpose of estimating total global glacier ice volume. I find that adjustments for continentality and elevation range improve skill of area–volume scaling. These scaling relationships are applied to each record in the Randolph Glacier Inventory, which is the first globally complete inventory of glaciers and ice caps. I estimate that the total volume of all glaciers in the world is 0.35 ± 0.07 m sea level equivalent, including ice sheet peripheral glaciers. This is substantially less than a recent state-of-the-art estimate. Area–volume scaling bias issues for large ice masses, and incomplete inventory data are offered as explanations for the difference.

Description: Ice tectonic deformation during the rapid in situ drainage of a supraglacial lake on the Greenland Ice Sheet The Cryosphere, 7, 129-140, 2013 Author(s): S. H. Doyle, A. L. Hubbard, C. F. Dow, G. A. Jones, A. Fitzpatrick, A. Gusmeroli, B. Kulessa, K. Lindback, R. Pettersson, and J. E. Box We present detailed records of lake discharge, ice motion and passive seismicity capturing the behaviour and processes preceding, during and following the rapid drainage of a 4 km 2 supraglacial lake through 1.1-km-thick ice on the western margin of the Greenland Ice Sheet. Peak discharge of 3300 m 3 s −1 coincident with maximal rates of vertical uplift indicates that surface water accessed the ice–bed interface causing widespread hydraulic separation and enhanced basal motion. The differential motion of four global positioning system (GPS) receivers located around the lake record the opening and closure of the fractures through which the lake drained. We hypothesise that the majority of discharge occurred through a 3-km-long fracture with a peak width averaged across its wetted length of 0.4 m. We argue that the fracture's kilometre-scale length allowed rapid discharge to be achieved by combining reasonable water velocities with sub-metre fracture widths. These observations add to the currently limited knowledge of in situ supraglacial lake drainage events, which rapidly deliver large volumes of water to the ice–bed interface.

Description: Glacier changes and climate trends derived from multiple sources in the data scarce Cordillera Vilcanota region, southern Peruvian Andes The Cryosphere, 7, 103-118, 2013 Author(s): N. Salzmann, C. Huggel, M. Rohrer, W. Silverio, B. G. Mark, P. Burns, and C. Portocarrero The role of glaciers as temporal water reservoirs is particularly pronounced in the (outer) tropics because of the very distinct wet/dry seasons. Rapid glacier retreat caused by climatic changes is thus a major concern, and decision makers demand urgently for regional/local glacier evolution trends, ice mass estimates and runoff assessments. However, in remote mountain areas, spatial and temporal data coverage is typically very scarce and this is further complicated by a high spatial and temporal variability in regions with complex topography. Here, we present an approach on how to deal with these constraints. For the Cordillera Vilcanota (southern Peruvian Andes), which is the second largest glacierized cordillera in Peru (after the Cordillera Blanca) and also comprises the Quelccaya Ice Cap, we assimilate a comprehensive multi-decadal collection of available glacier and climate data from multiple sources (satellite images, meteorological station data and climate reanalysis), and analyze them for respective changes in glacier area and volume and related trends in air temperature, precipitation and in a more general manner for specific humidity. While we found only marginal glacier changes between 1962 and 1985, there has been a massive ice loss since 1985 (about 30% of area and about 45% of volume). These high numbers corroborate studies from other glacierized cordilleras in Peru. The climate data show overall a moderate increase in air temperature, mostly weak and not significant trends for precipitation sums and probably cannot in full explain the observed substantial ice loss. Therefore, the likely increase of specific humidity in the upper troposphere, where the glaciers are located, is further discussed and we conclude that it played a major role in the observed massive ice loss of the Cordillera Vilcanota over the past decades.

Description: The impact of heterogeneous surface temperatures on the 2-m air temperature over the Arctic Ocean under clear skies in spring The Cryosphere, 7, 153-166, 2013 Author(s): A. Tetzlaff, L. Kaleschke, C. Lüpkes, F. Ament, and T. Vihma The influence of spatial surface temperature changes over the Arctic Ocean on the 2-m air temperature variability is estimated using backward trajectories based on ERA-Interim and JRA25 wind fields. They are initiated at Alert, Barrow and at the Tara drifting station. Three different methods are used. The first one compares mean ice surface temperatures along the trajectories to the observed 2-m air temperatures at the stations. The second one correlates the observed temperatures to air temperatures obtained using a simple Lagrangian box model that only includes the effect of sensible heat fluxes. For the third method, mean sensible heat fluxes from the model are correlated with the difference of the air temperatures at the model starting point and the observed temperatures at the stations. The calculations are based on MODIS ice surface temperatures and four different sets of ice concentration derived from SSM/I (Special Sensor Microwave Imager) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS) data. Under nearly cloud-free conditions, up to 90% of the 2-m air temperature variance can be explained for Alert, and 70% for Barrow, using these methods. The differences are attributed to the different ice conditions, which are characterized by high ice concentration around Alert and lower ice concentration near Barrow. These results are robust for the different sets of reanalyses and ice concentration data. Trajectories based on 10-m wind fields from both reanalyses show large spatial differences in the Central Arctic, leading to differences in the correlations between modeled and observed 2-m air temperatures. They are most pronounced at Tara, where explained variances amount to 70% using JRA and 80% using ERA. The results also suggest that near-surface temperatures at a given site are influenced by the variability of surface temperatures in a domain of about 200 km radius around the site.

Description: Future projections of the Greenland ice sheet energy balance driving the surface melt The Cryosphere, 7, 1-18, 2013 Author(s): B. Franco, X. Fettweis, and M. Erpicum In this study, simulations at 25 km resolution are performed over the Greenland ice sheet (GrIS) throughout the 20th and 21st centuries, using the regional climate model MAR forced by four RCP scenarios from three CMIP5 global circulation models (GCMs), in order to investigate the projected changes of the surface energy balance (SEB) components driving the surface melt. Analysis of 2000–2100 melt anomalies compared to melt results over 1980–1999 reveals an exponential relationship of the GrIS surface melt rate simulated by MAR to the near-surface air temperature (TAS) anomalies, mainly due to the surface albedo positive feedback associated with the extension of bare ice areas in summer. On the GrIS margins, the future melt anomalies are preferentially driven by stronger sensible heat fluxes, induced by enhanced warm air advection over the ice sheet. Over the central dry snow zone, the surface albedo positive feedback induced by the increase in summer melt exceeds the negative feedback of heavier snowfall for TAS anomalies higher than 4 °C. In addition to the incoming longwave flux increase associated with the atmosphere warming, GCM-forced MAR simulations project an increase of the cloud cover decreasing the ratio of the incoming shortwave versus longwave radiation and dampening the albedo feedback. However, it should be noted that this trend in the cloud cover is contrary to that simulated by ERA-Interim–forced MAR for recent climate conditions, where the observed melt increase since the 1990s seems mainly to be a consequence of more anticyclonic atmospheric conditions. Finally, no significant change is projected in the length of the melt season, which highlights the importance of solar radiation absorbed by the ice sheet surface in the melt SEB.

Description: Manufactured solutions and the verification of three-dimensional Stokes ice-sheet models The Cryosphere, 7, 19-29, 2013 Author(s): W. Leng, L. Ju, M. Gunzburger, and S. Price The manufactured solution technique is used for the verification of computational models in many fields. In this paper, we construct manufactured solutions for the three-dimensional, isothermal, nonlinear Stokes model for flows in glaciers and ice sheets. The solution construction procedure starts with kinematic boundary conditions and is mainly based on the solution of a first-order partial differential equation for the ice velocity that satisfies the incompressibility condition. The manufactured solutions depend on the geometry of the ice sheet, basal sliding parameters, and ice softness. Initial conditions are taken from the periodic geometry of a standard problem of the ISMIP-HOM benchmark tests. The upper surface is altered through the manufactured solution procedure to generate an analytic solution for the time-dependent flow problem. We then use this manufactured solution to verify a parallel, high-order accurate, finite element Stokes ice-sheet model. Simulation results from the computational model show good convergence to the manufactured analytic solution.

Description: Stable isotope and gas properties of two climatically contrasting (Pleistocene and Holocene) ice wedges from Cape Mamontov Klyk, Laptev Sea, northern Siberia The Cryosphere, 7, 31-46, 2013 Author(s): T. Boereboom, D. Samyn, H. Meyer, and J.-L. Tison This paper presents and discusses the texture, fabric, water stable isotopes (δ 18 O, δD) and gas properties (total gas content, O 2 , N 2 , Ar, CO 2 , and CH 4 mixing ratios) of two climatically contrasted (Holocene vs. Pleistocene) ice wedges (IW-26 and IW-28) from Cape Mamontov Klyk, Laptev Sea, in northern Siberia. The two ice wedges display contrasting structures: one being of relatively "clean" ice and the other showing clean ice at its centre as well as debris-rich ice on both sides (referred to as "ice-sand wedge"). Our multiparametric approach allows discrimination between three different ice facies with specific signatures, suggesting different climatic and environmental conditions of formation and various intensities and nature of biological activity. More specifically, crystallography, total gas content and gas composition reveal variable levels of meltwater infiltration and contrasting contributions from anaerobic and aerobic conditions to the biological signatures. Stable isotope data are drawn on to discuss changes in paleoenvironmental conditions and in the temporal variation of the different moisture sources for the snow feeding into the ice wedges infillings. Our data set also supports the previous assumption that the ice wedge IW-28 was formed in Pleistocene and the ice wedge IW-26 in Holocene times. This study sheds more light on the conditions of ice wedge growth under changing environmental conditions.

Description: Area and volume loss of the glaciers in the Ortles-Cevedale group (Eastern Italian Alps): controls and imbalance of the remaining glaciers The Cryosphere, 7, 1339-1359, 2013 Author(s): L. Carturan, R. Filippi, R. Seppi, P. Gabrielli, C. Notarnicola, L. Bertoldi, F. Paul, P. Rastner, F. Cazorzi, R. Dinale, and G. Dalla Fontana A widespread loss of glacier area and volume has been observed in the European Alps since the 1980s. In addition to differences among various regions of the Alps, different responses to climate change characterize neighboring glaciers within the same region. In this study we describe the glacier changes in the Ortles-Cevedale group, the largest glacierized area in the Italian Alps. We analyze the spatial variability, the drivers, and the main factors controlling the current loss of ice in this region, by comparing mean elevation changes derived from two digital terrain models (DTMs), along with glacier extents and snow-covered areas derived from Landsat images acquired in 1987 and 2009, to various topographic factors. Glacier outlines were obtained using the band ratio method with manual corrections. Snow was classified from a near-infrared image after topographic correction. The total glacierized area shrank by 23.4 ± 3% in this period, with no significant changes in the mean altitude of the glaciers. In 2009 the snowline was 240 m higher than in the 1960s and 1970s. From the snow-covered area at the end of summer 2009, which fairly represents the extent and local variability of the accumulation areas in the 2000s, we estimate that approximately 50% of the remaining glacier surfaces have to melt away to re-establish balanced mass budgets with present climatic conditions. The average geodetic mass budget rate, calculated for 112 ice bodies by differencing two DTMs, ranged from −0.18 ± 0.04 to −1.43 ± 0.09 m w.e. a −1 , averaging −0.69 ± 0.12 m w.e. a −1 . The correlation analysis of mass budgets vs. topographic variables emphasized the important role of hypsometry in controlling the area and volume loss of larger glaciers, whereas a higher variability characterizes smaller glaciers, which is likely due to the higher importance of local topo-climatic conditions.

Description: Brief Communication "The 2013 Erebus Glacier Tongue calving event" The Cryosphere, 7, 1333-1337, 2013 Author(s): C. L. Stevens, P. Sirguey, G. H. Leonard, and T. G. Haskell The Erebus Glacier Tongue, a small floating glacier in southern McMurdo Sound, is one of the best-studied ice tongues in Antarctica. Despite this, its calving on the 27 February 2013 (UTC) was around 10 yr earlier than previously predicted. The calving was likely a result of ocean currents and the absence of fast ice. The subsequent trajectory of the newly created iceberg supports previous descriptions of the surface ocean circulation in southern McMurdo Sound.

Description: Spectral reflectance of solar light from dirty snow: a simple theoretical model and its validation The Cryosphere, 7, 1325-1331, 2013 Author(s): A. Kokhanovsky A simple analytical equation for the snow albedo as the function of snow grain size, soot concentration, and soot mass absorption coefficient is presented. This simple equation can be used in climate models to assess the influence of snow pollution on snow albedo. It is shown that the squared logarithm of the albedo (in the visible) is directly proportional to the soot concentration. A new method of the determination of the soot mass absorption coefficient in snow is proposed. The equations derived are applied to a dusty snow layer as well.

Description: Brief communication "The aerophotogrammetric map of Greenland ice masses" The Cryosphere, 7, 445-449, 2013 Author(s): M. Citterio and A. P. Ahlstrøm The PROMICE (Programme for Monitoring of the Greenland Ice Sheet) aerophotogrammetric map of Greenland ice masses is the first high resolution dataset documenting the mid-1980s areal extent of the Greenland Ice Sheet and all the local glaciers and ice caps. The total glacierized area excluding nunataks was 1 804 638 km 2 ± 2178 km 2 , of which 88 083 ± 1240 km 2 belonged to local glaciers and ice caps (GIC) substantially independent from the Greenland Ice Sheet. This new result of GIC glacierized area is higher than most previous estimates, 81% greater than Weng's (1995) measurements, but is in line with contemporary findings based on independent data sources. A comparison between our map and the recently released Rastner et al. (2012) inventory and GIMP (Greenland Ice Mapping Project) Ice-Cover Mask (Howat and Negrete, 2013) shows potential for change-assessment studies.

Description: Speedup and fracturing of George VI Ice Shelf, Antarctic Peninsula The Cryosphere, 7, 797-816, 2013 Author(s): T. O. Holt, N. F. Glasser, D. J. Quincey, and M. R. Siegfried George VI Ice Shelf (GVIIS) is located on the Antarctic Peninsula, a region where several ice shelves have undergone rapid breakup in response to atmospheric and oceanic warming. We use a combination of optical (Landsat), radar (ERS 1/2 SAR) and laser altimetry (GLAS) datasets to examine the response of GVIIS to environmental change and to offer an assessment on its future stability. The spatial and structural changes of GVIIS (ca. 1973 to ca. 2010) are mapped and surface velocities are calculated at different time periods (InSAR and optical feature tracking from 1989 to 2009) to document changes in the ice shelf's flow regime. Surface elevation changes are recorded between 2003 and 2008 using repeat track ICESat acquisitions. We note an increase in fracture extent and distribution at the south ice front, ice-shelf acceleration towards both the north and south ice fronts and spatially varied negative surface elevation change throughout, with greater variations observed towards the central and southern regions of the ice shelf. We propose that whilst GVIIS is in no imminent danger of collapse, it is vulnerable to ongoing atmospheric and oceanic warming and is more susceptible to breakup along its southern margin in ice preconditioned for further retreat.

Description: High-resolution interactive modelling of the mountain glacier–atmosphere interface: an application over the Karakoram The Cryosphere, 7, 779-795, 2013 Author(s): E. Collier, T. Mölg, F. Maussion, D. Scherer, C. Mayer, and A. B. G. Bush The traditional approach to simulations of alpine glacier mass balance (MB) has been one-way, or offline, thus precluding feedbacks from changing glacier surface conditions on the atmospheric forcing. In addition, alpine glaciers have been only simply, if at all, represented in atmospheric models to date. Here, we extend a recently presented, novel technique for simulating glacier–atmosphere interactions without the need for statistical downscaling, through the use of a coupled high-resolution mesoscale atmospheric and physically-based climatic mass balance (CMB) modelling system that includes glacier CMB feedbacks to the atmosphere. We compare the model results over the Karakoram region of the northwestern Himalaya with remote sensing data for the ablation season of 2004 as well as with in situ glaciological and meteorological measurements from the Baltoro glacier. We find that interactive coupling has a localized but appreciable impact on the near-surface meteorological forcing data and that incorporation of CMB processes improves the simulation of variables such as land surface temperature and snow albedo. Furthermore, including feedbacks from the glacier model has a non-negligible effect on simulated CMB, reducing modelled ablation, on average, by 0.1 m w.e. (−6.0%) to a total of −1.5 m w.e. between 25 June–31 August 2004. The interactively coupled model shows promise as a new, multi-scale tool for explicitly resolving atmospheric-CMB processes of mountain glaciers at the basin scale.

Description: Brief communication "Global glacier volumes and sea level – small but systematic effects of ice below the surface of the ocean and of new local lakes on land" The Cryosphere, 7, 817-821, 2013 Author(s): W. Haeberli and A. Linsbauer The potential contribution of glaciers and ice caps to sea level rise is usually calculated by comparing the estimated total ice volume with the surface area of the ocean. Part of this total ice volume, however, does not contribute to sea level rise because it is below the surface of the ocean or below the levels of future lakes on land. The present communication points to this so far overlooked phenomenon and provides a first order-of-magnitude estimate. It is shown that the effect is small (most likely about 1 to 6 cm sea level equivalent) but systematic, could primarily affect earlier stages of global glacier vanishing, and should therefore be adequately considered. Now-available techniques of slope-related high-resolution glacier bed modelling have the potential to provide more detailed assessments in the future.