This paper presents the changes in the thermal structure of the polythermal glacier Storglaciären, northern Sweden, over the 20 year period 1989-2009 derived by comparing maps of the depth of the englacial transition between cold ice (permanently frozen) and temperate ice (which contains water inclusions). The maps are based on interpreted ice-penetrating radar surveys from 1989, 2001 and 2009.

Complex thinning of the cold layer, first identified between 1989 and 2001, is still ongoing. A volume calculation shows that Storglaciären has lost one-third of its cold surface layer volume in 20 years, with a mean thinning rate of 0.80 ± 0.24 m a-1. We suggest that the thinning of the cold layer at Storglaciären is connected to the climatic warming experienced by sub-Arctic Scandinavia since the 1980s and we argue that repeated ice-penetrating radar surveys over the ablation area of polythermal glaciers offer a useful proxy for evaluating glacier responses to changes in climate.

Layered and polycrystalline ice was experimentally deformed in general shear involving axial compression (strain magnitude 0.5-17%) and simple shear (strain magnitude γ = 0.1-1.4). As the temperature is increased from -20°C to -2°C, there is at least a twofold enhancement in octahedral shear strain rate, which coincides with the onset of extensive dynamic recrystallization and a change in grain-size distribution at -15°C. Between -150C and -10°C the c-axis preferred orientation rapidly evolves with the initiation of two-maxima fabrics in shear zones. From -10°C to -2°C there is progressive evolution of a final c-axis pattern that is asymmetric with respect to the direction of shortening, with a strong maximum at ~5° to the pole of the shear zone, a sense of asymmetry in the direction of the shear, and a secondary maximum inclined at ~45° to the plane of shearing. An initial c-axis preferred orientation plays a critical role in the initial mechanical evolution. In contrast to established ideas, a strong alignment of basal planes parallel to the plane of easy glide inhibited deformation and there was an increased component of strain hardening until recrystallization processes become dominant.

During the early melt season, when sea-ice extent begins to retreat due to increasing heat absorption from the atmosphere and ocean, the sea ice becomes porous and weakens. This induces the break-up of ice floes and accelerates the melting process substantially. Determining how sea-ice structure evolves during this season is important for understanding the melting process on a global scale. We investigated this issue using field observations, laboratory experiments and numerical modeling, focusing on the effect of the C-shaped temperature profile on the internal structure of sea ice. Field observations were conducted to examine this effect on ~0.3 m thick ice on Lake Saroma, located at the coast of Hokkaido, Japan. To test the generalization from these observational case studies, laboratory experiments were conducted using a tank with 0.12 m thick ice. Additionally, evolution of ice structure from winter to early spring was numerically analyzed using meteorological data. From these studies, we find that the heat convergence, caused by the C-shape temperature profile, is essential to the internal melting particularly in the upper layer. This increases the porosity of the sea ice, leading to a rapid decrease in its flexural strength during the season.

Near ice shelves around Antarctica the ocean becomes supercooled and has been observed to carry small suspended ice crystals. Our measurements demonstrate that these small crystals are persistently present in the water column beneath the winter fast ice, and when incorporated in sea ice they reduce the mean grain size of the sea-ice cover. By midwinter, larger ice crystals below the ice/water interface are observed to form a porous sub-ice platelet layer with an ice volume fraction of 0.25 ± 0.06. The magnitude and direction of the oceanic heat flux varied between (5 ± 6) Wm-2 (upwards) and (-15 ± 10) Wm-2 (downwards) in May, but by September it settled between (-6 ± 2) and (-11 ± 2) W m-2. The negative values imply that the ocean acts as a heat sink which is responsible for the growth of 12% of the ice thickness between June and September. This oceanic contribution should not be ignored in models of Antarctic sea-ice thickness close to an ice shelf.

Ice velocities observed in 2005/06 at three GPS stations along the Sermeq Avannarleq flowline, West Greenland, are used to characterize an observed annual velocity cycle. We attempt to reproduce this annual ice velocity cycle using a 1-D ice-flow model with longitudinal stresses coupled to a 1-D hydrology model that governs an empirical basal sliding rule. Seasonal basal sliding velocity is parameterized as a perturbation of prescribed winter sliding velocity that is proportional to the rate of change of glacier water storage. The coupled model reproduces the broad features of the annual basal sliding cycle observed along this flowline, namely a summer speed-up event followed by a fall slowdown event. We also evaluate the hypothesis that the observed annual velocity cycle is due to the annual calving cycle at the terminus. We demonstrate that the ice acceleration due to a catastrophic calving event takes an order of magnitude longer to reach CU/ETH (‘Swiss’) Camp (46 km upstream of the terminus) than is observed. The seasonal acceleration observed at Swiss Camp is therefore unlikely to be the result of velocity perturbations propagated upstream via longitudinal coupling. Instead we interpret this velocity cycle to reflect the local history of glacier water balance.

A detailed history of volcanism covering the last 2840 years is reconstructed from the top 100.42 m of a 109.91 m ice core from Dome A (DA2005 ice core), East Antarctica. Using two known volcanic stratigraphic markers, the mean accumulation rate during the period AD 1260-1964 is found to be 23.2 mmw.e. a-1, consistent with the previously reported accumulation rate at Dome A. This mean accumulation rate is used to date the entire core. Volcanic eruptions in the period 840 BC-AD1998 are detected as outstanding sulphate events. Seventy-eight eruptions are identified, with a mean of 2.7 eruptions per century. Comparisons with previous Antarctic ice-core volcanic records are made to assess the quality of this new DA2005 record. In terms of dates for volcanic events, the DA2005 record is in good agreement with previous records in the second millennium ad (ad 1000-1998). A series of volcanic signatures found in both the DA2005 record and several other Antarctic ice-core records in the first millennium ad (ad 1-1000) appear to validate the DA2005 record during this time period. For the older periods, direct comparisons are difficult between the DA2005 record and other Antarctic ice-core records due to the lack of well-dated stratigraphic horizons.

Higher-order models represent a computationally less expensive alternative to the Stokes model for ice-sheet modeling. In this work, we develop linear and quadratic finite-element methods, implemented on parallel architectures, for the three-dimensional first-order model of Dukowicz and others (2010) that is based on the Blatter-Pattyn model, and for the depth-integrated model of Schoof and Hindmarsh (2010). We then apply our computational models to three of the ISMIP-HOM benchmark test cases (Pattyn and others, 2008). We compare results obtained from our models with those obtained using a reliable Stokes computational model, showing that our first-order model implementation produces reliable and accurate solutions for almost all characteristic length scales of the test geometries considered. Good agreement with the reference Stokes solution is also obtained by our depth-integrated model implementation in fast-sliding regimes and for medium to large length scales. We also provide a comprehensive comparison between results obtained from our first-order model implementation and implementations developed by ISMIP-HOM participants; this study shows that our implementation is at least as good as the previous ones. Finally, a comparison between linear and quadratic finite- element approximations is carried out, showing, as expected, the better accuracy of the quadratic finite-element method.

A laboratory device was built to measure the forces that ice exerts on a 0.05 m diameter rigid plastic sphere in two different configurations: in contact with a flat bed or isolated from the bed. Measurements indicated that bed-normal contact forces were 1.8 times larger than drag forces due to creeping flow past a slippery sphere isolated from the bed. Measurements of forces as a function of the bed-normal ice velocity, estimations of the ice viscosity parameter and observations of markers in the ice indicate ice is Newtonian with a viscosity of ~1.3x 1011 Pas. Newtonian behavior is expected due to small and transient stresses. A model of regelation indicates that it had a negligible (<5%) influence on forces. Water pressure in the cavity beneath the sphere in contact with the bed had a likewise negligible influence on contact forces. When no cavity is present, drag forces can be correctly estimated using Stokes’s law (Newtonian viscosity) for a slippery sphere. The same law with a bed-enhancement factor of 1.8 is appropriate for estimating bed-normal contact forces. These results reinforce previous laboratory measurements and theories but provide no support for explanations of high debris/bed friction or rates of abrasion that depend on high contact forces.

We used ground-penetrating radar (GPR), GPS and glaciochemistry to evaluate melt regimes and ice depths, important variables for mass-balance and ice-volume studies, of Upper Yentna Glacier, Upper Kahiltna Glacier and the Mount Hunter ice divide, Alaska. We show the wet, percolation and dry snow zones located below ~2700ma.s.l., at ~2700 to 3900ma.s.l. and above 3900ma.s.l., respectively. We successfully imaged glacier ice depths upwards of 480 m using 40-100 MHz GPR frequencies. This depth is nearly double previous depth measurements reached using mid-frequency GPR systems on temperate glaciers. Few Holocene-length climate records are available in Alaska, hence we also assess stratigraphy and flow dynamics at each study site as a potential ice-core location. Ice layers in shallow firn cores and attenuated glaciochemical signals or lacking strata in GPR profiles collected on Upper Yentna Glacier suggest that regions below 2800ma.s.l. are inappropriate for paleoclimate studies because of chemical diffusion, through melt. Flow complexities on Kahiltna Glacier preclude ice-core climate studies. Minimal signs of melt or deformation, and depth-age model estimates suggesting ~4815 years of ice on the Mount Hunter ice divide (3912ma.s.l.) make it a suitable Holocene-age ice-core location.

The evolution of ice dynamics through surges of four tidewater-terminating glaciers in northwest Svalbard is investigated by remote sensing. A 20 year time series of glacier surface flow speeds and frontal positions is presented covering the recent surges of Monacobreen, Comfortlessbreen, Blomstrandbreen and Fjortende Julibreen. Surface flow speeds were derived using feature tracking between pairs of ERS SAR and ALOS PALSAR images, while frontal positions were taken from the same imagery, as well as more frequent but lower-spatial-resolution Envisat Wide Swath Mode images. During all four surges, increased ice flow caused the tidewater margin to advance while the calving flux was initially reduced to near zero due to compressive stresses limiting crevasse propagation. As ice speed decreased, the terminus continued to advance, until the glacier’s speed had returned to its pre-surge flow rate. Only at this time did the terminus start to retreat and peak iceberg calving flux was established. We conclude that terminus advance closely tracks glacier speed-up, that there is little mass loss through calving during the most active phase of the surge, and that seasonal cycles of terminus positions diminish during the active surge phase.

This study independently quantifies geodetic elevation change and models surface mass balance to solve the continuity equation. The approach is tested on two dynamically different glaciers, Kongsvegen and Kronebreen in northwest Svalbard, through two time epochs, 1966-1990/95 (I) and 1990/95-2007 (II). On Kongsvegen, a dynamically inactive glacier, the residual term represents an error associated with determining elevation changes and surface mass balance. It is apparent that centerline mass-balance estimates are not representative of the entire glacier, which we relate to center-line accumulation being larger than the elevation bin average. On Kronebreen, a fast-flowing and actively calving glacier, a significant part of the residual is identified with the long-term calving flux. For both glaciers, the cumulative surface mass balance remained close to zero during the first epoch but became increasingly negative in the second epoch. The long-term calving flux of Kronebreen is estimated to be -0.14 ± 0.03 km3 w.e.a-1 during epoch I and-0.20 ± 0.05 km3 w.e.a-1 in epoch II.

The tropical glaciers of the Cordillera Blanca, Peru, are rapidly retreating, resulting in complex impacts on the hydrology of the upper Río Santa watershed. The effect of this retreat on water resources is evaluated by analyzing historical and recent time series of daily discharge at nine measurement points. Using the Mann-Kendall nonparametric statistical test, the significance of trends in three hydrograph parameters was studied. Results are interpreted using synthetic time series generated from a hydrologic model that calculates hydrographs based on glacier retreat sequences. The results suggest that seven of the nine study watersheds have probably crossed a critical transition point, and now exhibit decreasing dry-season discharge. Our results suggest also that once the glaciers completely melt, annual discharge will be lower than present by 2-30% depending on the watershed. The retreat influence on discharge will be more pronounced during the dry season than at other periods of the year. At La Balsa, which measures discharge from the upper Río Santa, the glacier retreat could lead to a decrease in dry-season average discharge of 30%.

Accurate knowledge of the spatial distribution of the mass balance of temperate glaciers is essential for a better understanding of the physical processes controlling the mass balance and for the monitoring of water resources. In relation to albedo variations, the shortwave radiation budget is a controlling variable of the surface energy balance of glaciers. Remotely sensed albedo observations are here assimilated in a snowpack model to improve the modeling of the spatial distribution of the glacier mass balance. The albedo observations are integrated in the snowpack simulation using a variational data assimilation scheme that modifies the surface grain conditions. The study shows that mesoscale meteorological variables and MODIS-derived albedo maps can be used to obtain a good reconstruction of the annual mass balance on Glacier de Saint-Sorlin, French Alps, on a 100 m × 100m grid. Five hydrological years within the 2000-10 decade are tested. The accuracy of the method is estimated from comparison with field measurements. Sensitivity to roughness lengths and winter precipitation fields is investigated. Results demonstrate the potential contribution of remote-sensing data and variational data assimilation to further improve the understanding and monitoring of the mass balance of snowpacks and temperate glaciers.

Recent observations and modeling studies emphasize the crucial role of fracture mechanics for the stability of ice shelves and thereby the evolution of ice sheets. Here we introduce a macroscopic fracture-density field into a prognostic continuum ice-flow model and compute its evolution incorporating the initiation and growth of fractures as well as their advection with two dimensional ice flow. To a first approximation, fracture growth is assumed to depend on the spreading rate only, while fracture initiation is defined in terms of principal stresses. The inferred fracture-density fields compare well with observed elongate surface structures. Since crevasses and other deep-reaching fracture structures have been shown to influence the overall ice-shelf dynamics, we propose the fracture- density field introduced here be used as a measure for ice softening and decoupling of the ice flow in fracture-weakened zones. This may yield more accurate and realistic velocity patterns in prognostic simulations. Additionally, the memory of past fracture events links the calving front to the upstream dynamics. Thus the fracture-density field proposed here may be employed in fracture-based calving parameterizations. The aim of this study is to introduce the field and investigate which of the observed surface structures can be reproduced by the simplest physically motivated fracture source terms.

Remote-sensing and GIS techniques in conjunction with field investigations show how glacier mass loss has led to the rapid growth of Linggo Co, a glacier-fed lake on the central Tibetan Plateau, which has expanded by 21.3% in area between 1974 and 2010, with a lake-level rise of ˜11.2m. The lake volume of Linggo Co increased at a rate of 0.02 × 106, 42.67 × 106 and 65.8 × I06m3a-1 during the periods 1974-92, 1992-99 and 1999-2010, respectively. Other nonglacier-fed lakes in the vicinity (i.e. Longwei Co, Amur Co and Darngo Co Ngion) shrank considerably from the early 1970s to 1992 and then expanded from 1992 to 2010. Despite being in the same climate region, Linggo Co and the non-glacier-fed lakes have differed in response to climate change. The glaciers in the catchment of Linggo Co retreated by 2.4% in area between 1974 and 2007, and their mean thickness decreased by 6.19 ± 1.91 m between 1974 and 2000, with an associated glacier meltwater runoff of (7.52 ± 2.32) × 108 m3. The results indicate that glacier mass loss had a significant impact on the growth of Linggo Co over the past 40 years.

We interpreted flow dynamics of the Kahiltna Pass Basin accumulation zone on Mount McKinley, Alaska, USA, using 40,100 and 900 MHz ground-penetrating radar profiles and GPS surface velocity measurements. We found dipping, englacial surface-conformable strata that experienced vertical thickening as the glacier flowed westward from a steep, higher-velocity (60 m a-1) region into flat terrain associated with a 90° bend in the glacier and lower velocities (15 m a-1) to the south. Stratigraphy near the western side of the basin was surface-conformable to ˜170m depth and thinned as flow diverged southward, down-glacier. We found complex strata beneath the conformable stratigraphy and interpret these features as buried crevasses, avalanche debris and deformed ice caused by up-glacier events. We also suggest that basin dimensions, bed topography and the sharp bend each cause flow extension and compression, significantly contributing to conformable and complex strata thickness variations. Our findings show that surface-conformable stratigraphy continuous with depth and consistent strata thicknesses cannot be assumed in accumulation basins, because local and up- glacier terrain and flow dynamics can cause structural complexities to occur under and within surface- conformable layers.