Showing papers in "Journal of Glaciology in 2012"

Glacier recession and water resources in Peru's Cordillera Blanca

Abstract: The tropical glaciers of the Cordillera Blanca, Peru, are rapidly retreating, resulting in complex impacts on the hydrology of the upper Rio 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 Rio Santa, the glacier retreat could lead to a decrease in dry-season average discharge of 30%.

An enthalpy formulation for glaciers and ice sheets

Abstract: Polythermal conditions are ubiquitous among glaciers, from small valley glaciers to ice sheets. Conventional temperature-based 'cold-ice' models of such ice masses cannot account for that portion of the internal energy which is latent heat of liquid water within temperate ice, so such schemes are not energy-conserving when temperate ice is present. Temperature and liquid water fraction are, however, functions of a single enthalpy variable: a small enthalpy change in cold ice is a change in temperature, while a small enthalpy change in temperate ice is a change in liquid water fraction. The unified enthalpy formulation described here models the mass and energy balance for the three- dimensional ice fluid, for the surface runoff layer and for the subglacial hydrology layer, together in a single energy-conserving theoretical framework. It is implemented in the Parallel Ice Sheet Model. Results for the Greenland ice sheet are compared with those from a cold-ice scheme. This paper is intended to be an accessible foundation for enthalpy formulations in glaciology.

Impact of resolution and radar penetration on glacier elevation changes computed from DEM differencing

Abstract: Regional glacier mass balances can be measured by subtracting multi-temporal digital elevation models (DEMs). However, DEMs are often biased with altitude and it remains unclear whether the elevation differences observed on the ice-free terrain can be used to correct biases on ice-covered areas. IWe investigate such altitude-related biases using DEMs from three different sensors: SPOT5, SRTM C-band and SRTM X-band. The bias due to different original resolutions can be corrected using a relationship between curvature and elevation difference, calculated on ice-free terrain. The impact of C-band radar penetration into snow and ice can be evaluated for a specific region by comparing SRTM C-band and SRTM X-band DEMs. In our test area (Karakoram), the resolution-related bias has a minor influence on region-averaged elevation change. Conversely, not accounting for C-Band penetration would significantly biased the mean elevation change.

Dynamic thinning of Antarctic glaciers from along-track repeat radar altimetry

Abstract: Since 2002, the Envisat radar altimeter has measured the elevation of the Antarctic ice sheet with a repeat cycle of 35 days. This long and regular time series is processed using an along-track algorithm to depict in detail the spatial and temporal pattern of elevation change for the whole ice sheet. We use this dataset to examine the spatial and temporal pattern of Pine Island Glacier (PIG) thinning and compare it to the neighbouring glaciers. We also examine additional areas, especially in East Antarctica whose mass balance is poorly known. One advantage of the finer along-track spacing of measurements is that it reveals places of dynamic thinning in regions of rapid ice flow. We observe the acceleration of thinning on PIG. Over the entire basin, the volume loss increased from 7 km3 a-1 during 2002-06 to ∼48 km3 a-1 during 2006-10. We also observe accelerated thinning on the lower tens of kilometres of Thwaites Glacier, with a mean thinning of 0.18 m a-1 over its entire basin during our observation period. We confirm the dynamic thinning of Totten Glacier but we do not detect significantly accelerated thinning on any glacier elsewhere than on the coast of the Amundsen Sea.

From balance to imbalance: a shift in the dynamic behaviour of Chhota Shigri glacier, western Himalaya, India

Abstract: Mass-balance and dynamic behaviour of Chhota Shigri Glacier have been investigated between 2002 and 2010 and compared to data collected in 1987/1989. During the period 2002/2010, the glacier experienced a negative glacier-wide mass balance of -0.67 ± 0.40 m/a w.e. Between 2003 and 2010, elevation and ice flow velocities are slowly decreasing in the ablation area leading to a 24 to 37% reduction in ice fluxes, an expected response of the glacier dynamics to its recent negative mass balances. The reduced ice fluxes still remain far larger than the balance fluxes calculated from the year 2002 to 2010 average surface mass balances. Therefore, further slow down, thinning and terminus retreat of Chhota Shigri Glacier are expected over the next years. Conversely, the 2003/2004 ice fluxes are in good agreement with ice fluxes calculated assuming that the glacier-wide mass balance is zero. Given the limited velocity change between 1987/1989 and 2003/2004 and the small terminus change between 1988 and 2010, we suggest that the glacier has experienced a period of near zero or slightly positive mass balance in the 1990s, before shifting to a strong imbalance in the 21st century. This result challenges the generally accepted idea that glaciers of Western Himalaya have been shrinking rapidly for the last decades.

Elevation changes of glaciers revealed by multitemporal digital elevation models calibrated by GPS survey in the Khumbu region, Nepal Himalaya, 1992-2008

Abstract: Due to remoteness and high altitude, only a few ground-based glacier change studies are available in high-mountain areas in the Himalaya. However, digital elevation models based on remotely sensed data (RS-DEMs) provide feasible opportunities to evaluate how fast Himalayan glaciers are changing. Here we compute elevation changes in glacier surface (total area 183.3 km2) in the Khumbu region, Nepal Himalaya, for the period 1992-2008 using multitemporal RS-DEMs and a map-derived DEM calibrated with differential GPS survey data in 2007. Elevation change is calculated by generating a weighted least-squares linear regression model. Our method enables us to provide the distribution of uncertainty of the elevation change. Debris-covered areas show large lowering rates. The spatial distribution of elevation change shows that the different wastage features of the debris-covered glaciers depend on their scale, slope and the existence of glacial lakes. The elevation changes of glaciers in the eastern Khumbu region are in line with previous studies. The regional average mass balance of -0.40 ± 0.25 m w.e.a-1 for the period 1992-2008 is consistent with a global value of about -0.55 m w.e. a-1 for the period 1996-2005.

Accelerating shrinkage of Patagonian glaciers from the Little Ice Age (~AD 1870) to 2011

Abstract: We used Little Ice Age (LIA) trimlines and moraines to assess changes in South American glaciers over the last ∼140 years. We determined the extent and length of 640 glaciers during the LIA (∼ AD 1870) and 626 glaciers (the remainder having entirely disappeared) in 1986, 2001 and 2011. The calculated reduction in glacierized area between the LIA and 2011 is 4131 km2 (15.4%), with 660 km2 (14.2%) being lost from the Northern Patagonia Icefield (NPI), 1643 km2 (11.4%) from the Southern Patagonia Icefield (SPI) and 306 km2 (14.4%) from Cordillera Darwin. Latitude, size and terminal environment (calving or land-terminating) exert the greatest control on rates of shrinkage. Small, northerly, land-terminating glaciers shrank fastest. Annual rates of area loss increased dramatically after 2001 for mountain glaciers north of 52° S and the large icefields, with the NPI and SPI now shrinking at 9.4 km2 a–1 (0.23% a–1) and 20.5 km2 a–1 (0.15% a–1) respectively. The shrinkage of glaciers between 52° S and 54° S accelerated after 1986, and rates of shrinkage from 1986 to 2011 remained steady. Icefield outlet glaciers, isolated glaciers and ice caps south of 54° S shrank faster from 1986 to 2001 than they did from 2001 to 2011.

Using surface velocities to calculate ice thickness and bed topography: a case study at Columbia Glacier, Alaska, USA

Abstract: Information about glacier volume and ice thickness distribution is essential for many glaciological applications, but direct measurements of ice thickness can be difficult and costly. We present a new method that calculates ice thickness via an estimate of ice flux. We solve the familiar continuity equation between adjacent flowlines, which decreases the computational time required compared to a solution on the whole grid. We test the method on Columbia Glacier, a large tidewater glacier in Alaska, USA, and compare calculated and measured ice thicknesses, with favorable results. This shows the potential of this method for estimating ice thickness distribution of glaciers for which only surface data are available. We find that both the mean thickness and volume of Columbia Glacier were approximately halved over the period 1957-2007, from 281m to 143m, and from 294km 3 to 134km 3 , respectively. Using bedrock slope and considering how waves of thickness change propagate through the glacier, we conduct a brief analysis of the instability of Columbia Glacier, which leads us to conclude that the rapid portion of the retreat may be nearing an end.

Widespread rifting and retreat of ice-shelf margins in the eastern Amundsen Sea Embayment between 1972 and 2011

Abstract: The major outlet glaciers that drain the eastern sector of the Amundsen Sea Embayment (Smith, Haynes, Thwaites and Pine Island) are among the largest, fastest-flowing and fastest-thinning glaciers in West Antarctica. Their recent ice-flow acceleration is linked to ocean-induced ice-shelf thinning, but may also arise from additional losses of ice-shelf buttressing that are not well understood. Here we present a comprehensive history of coastal change in the eastern Amundsen Sea Embayment between 1972 and 2011 derived mostly from Landsat imagery. The termini of all four major outlet glaciers have retreated, but retreat is most rapid along the ice-shelf margins, where progressive rifting has occurred. This pattern of retreat coincides with the recent acceleration of grounded ice and contributed to loss of ice-shelf buttressing. The observed pattern of margin-led gradual ice-shelf disintegration appears to be common in accelerating ocean-terminating outlet glaciers. We hypothesize that this pattern is part of a positive feedback between glacier acceleration and rift growth that could drive further buttressing loss in the eastern Amundsen Sea Embayment.

Simulations of the Greenland ice sheet 100 years into the future with the full Stokes model Elmer/Ice

Abstract: It is likely that climate change will have a significant impact on the mass balance of the Greenland ice sheet, contributing to future sea-level rise. Here we present the implementation of the full Stokes model Elmer/Ice for the Greenland ice sheet, which includes a mesh refinement technique in order to resolve fast-flowing ice streams and outlet glaciers. We discuss simulations 100 years into the future, forced by scenarios defined by the SeaRISE (Sea-level Response to Ice Sheet Evolution) community effort. For comparison, the same experiments are also run with the shallow-ice model SICOPOLIS (SImulation COde for POLythermal Ice Sheets). We find that Elmer/Ice is ~43% more sensitive (exhibits a larger loss of ice-sheet volume relative to the control run) than SICOPOLIS for the ice-dynamic scenario (doubled basal sliding), but ~61 % less sensitive for the direct global warming scenario (based on the A1 B moderate-emission scenario for greenhouse gases). The scenario with combined A1B global warming and doubled basal sliding forcing produces a Greenland contribution to sea-level rise of ~15cm for Elmer/Ice and ~12cm for SICOPOLIS over the next 100 years.

The response of Petermann Glacier, Greenland, to large calving events, and its future stability in the context of atmospheric and oceanic warming

Abstract: This study assesses the impact of a large 2010 calving event on the current and future stability of Petermann Glacier, Greenland, and ascertains the glacier's interaction with different components of the climate and ocean system. We use a numerical ice-flow model that captures the major aspects of the glacier's mass budget, the resistive forces controlling glacier flow, and includes dynamic calving. Satellite observations and model results show that the recent break-off of 25% of the floating tongue did not result in a significant glacier speed-up due to the low lateral resistance of this relatively wide and thin ice tongue. We demonstrate that seasonal speed-up at Petermann Glacier is mainly driven by meltwater lubrication rather than freeze-up conditions in the fjord. Results also show that sub-shelf ocean melt may have a profound effect on the future stability of Petermann Glacier, emphasizing the urgent need for more observations, and a better understanding of fjord temperature variability and circulation.

An inexact Gauss-Newton method for inversion of basal sliding and rheology parameters in a nonlinear Stokes ice sheet model

Abstract: We propose an infinite-dimensional adjoint-based inexact Gauss-Newton method for the solution of inverse problems governed by Stokes models of ice sheet flow with nonlinear rheology and sliding law. The method is applied to infer the basal sliding coefficient and the rheological exponent parameter fields from surface velocities. The inverse problem is formulated as a nonlinear least-squares optimization problem whose cost functional is the misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to the cost functional to render the problem well-posed and account for observational error. Our findings show that the inexact Newton method is significantly more efficient than the nonlinear conjugate gradient method and that the number of Stokes solutions required to solve the inverse problem is insensitive to the number of inversion parameters. The results also show that the reconstructions of the basal sliding coefficient converge to the exact sliding coefficient as the observation error (here, the noise added to synthetic observations) decreases, and that a nonlinear rheology makes the reconstruction of the basal sliding coefficient more difficult. For the inversion of the rheology exponent field, we find that horizontally constant or smoothly varying parameter fields can be reconstructed satisfactorily from noisy observations.

Can the snowline be used as an indicator of the equilibrium line and mass balance for glaciers in the outer tropics

Abstract: Original abstract: Because the glacier snowline is easy to identify on optical satellite images and because in certain conditions it can be used as an indicator of the equilibrium line, it may be a relevant parameter for the study of the relationships between climate and glaciers. Although several studies have shown that the snowline altitude (SLA) at the end of the hydrological year is a good indicator of the equilibrium-line altitude (ELA) for mid-latitude glaciers, such a relationship remains conjectural for tropical glaciers. Indeed, unlike in mid-latitudes, tropical climate conditions result in a distinct seasonality of accumulation/ablation processes. We examine this relationship using direct field ELA and mass-balance measurements made on Glaciar Zongo, Bolivia (~16° S), vand Glaciar Artesonraju, Peru (~9° S), and the SLA retrieved from satellite images acquired in the past two decades. We show that on glaciers in the outer tropics: (1) ablation is reduced during the dry season in austral winter (May-August), the SLA does not change much, and satellite images acquired between May and August could be used to compute the SLA; and (2) the highest SLA detected on a number of satellite images acquired during the dry season provides a good estimate of the annual ELA. However, as snowfall events can occur during the dry season, the SLA detected on satellite images tends to underestimate the ELA. Thus, we recommend validating the SLA computed from satellite images with field data collected on a benchmark glacier before measuring the SLA on other glaciers in the same mountain range for which no field data are available. This study is a major step towards extending the measurement of glacier parameters (ELA and mass balance) at the scale of a whole mountain range in the outer tropics to better document the relationships between climate and glaciers.

Ice-shelf basal channels in a coupled ice/ocean model

Abstract: A numerical model for an interacting ice shelf and ocean is presented in which the iceshelf base exhibits a channelized morphology similar to that observed beneath Petermann Gletscher’s (Greenland) floating ice shelf. Channels are initiated by irregularities in the ice along the grounding line and then enlarged by ocean melting. To a first approximation, spatially variable basal melting seaward of the grounding line acts as a steel-rule die or a stencil, imparting a channelized form to the ice base as it passes by. Ocean circulation in the region of high melt is inertial in the along-channel direction and geostrophically balanced in the transverse direction. Melt rates depend on the wavelength of imposed variations in ice thickness where it enters the shelf, with shorter wavelengths reducing overall melting. Petermann Gletscher’s narrow basal channels may therefore act to preserve the ice shelf against excessive melting. Overall melting in the model increases for a warming of the subsurface water. The same sensitivity holds for very slight cooling, but for cooling of a few tenths of a degree a reorganization of the spatial pattern of melting leads, surprisingly, to catastrophic thinning of the ice shelf 12 km from the grounding line. Subglacial discharge of fresh water along the grounding line increases overall melting. The eventual steady state depends on when discharge is initiated in the transient history of the ice, showing that multiple steady states of the coupled system exist in general.

A physically based method for estimating supraglacial debris thickness from thermal band remote-sensing data

Abstract: In order to account for the effects of debris cover in model scenarios of the response of glaciers to climate change and water resource planning, it is important to know the distribution and thickness of supraglacial debris and to monitor its change over time. Previous attempts to map surface debris thickness using thermal band remote sensing have relied upon time-specific empirical relationships between surface temperature and thickness, limiting their general applicability. In this paper, we develop a physically based model that utilizes Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) thermal band remotely sensed imagery and is based on a solution of the energy balance at the debris surface. The model is used to estimate debris thickness on Miage glacier, Italy, and is validated using field debris-thickness measurements and a previously published debris-thickness map. The temporal transferability of the model is demonstrated through successful application to a separate ASTER image from a different year using reanalysis meteorological input data. This model has the potential to be used for regional-scale supraglacial debris-thickness mapping and monitoring for debris up to at least 0.50 m thickness, but improved understanding of the spatial patterns of air temperature, aerodynamic roughness length and thermal properties across debriscovered glaciers is needed.

Variability in the mass flux of the Ross ice streams, West Antarctica, over the last millennium

Abstract: We synthesize previously published remote-sensing observations, radar data and model output to obtain a 1000 year ice flow history for the Siple Coast ice-stream system in West Antarctica to investigate the timing and magnitude of changes in mass flux. The synthesized history shows significant short-term variability in ice-stream shear margin and grounding line position due to internal variability of the coupled system. The chronology highlights the interplay between adjacent ice streams, which implies that the behavior of any individual ice stream should not be examined in isolation. Furthermore, individual events cannot be fully interpreted without an understanding of the broad-scale, long-term variability in the ice sheet. In the context of this millennium-scale history, we interpret the relatively recent stagnation of Kamb Ice Stream (KIS) as just one stage in the thermodynamic cycle of an ice stream in this region. The changes in mass balance that result from the KIS stagnation may thus be viewed as century-scale 'noise' relative to the longer-term trend. Understanding and characterizing this noise is a necessary step before accurate model-based predictions of ice-sheet mass balance for the next century can be made.

Extent of Low-accumulation 'Wind Glaze' Areas on the East Antarctic Plateau: Implications for Continental Ice Mass Balance

Abstract: Persistent katabatic winds form widely distributed localized areas of near-zero net surface accumulation on the East Antarctic ice sheet (EAIS) plateau. These areas have been called 'glaze' surfaces due to their polished appearance. They are typically 2-200 square kilometers in area and are found on leeward slopes of ice-sheet undulations and megadunes. Adjacent, leeward high-accumulation regions (isolated dunes) are generally smaller and do not compensate for the local low in surface mass balance (SMB). We use a combination of satellite remote sensing and field-gathered datasets to map the extent of wind glaze in the EAIS above 1500m elevation. Mapping criteria are derived from distinctive surface and subsurface characteristics of glaze areas resulting from many years of intense annual temperature cycling without significant burial. Our results show that 11.2 plus or minus 1.7%, or 950 plus or minus 143 x 10(exp 3) square kilometers, of the EAIS above 1500m is wind glaze. Studies of SMB interpolate values across glaze regions, leading to overestimates of net mass input. Using our derived wind-glaze extent, we estimate this excess in three recent models of Antarctic SMB at 46-82 Gt. The lowest-input model appears to best match the mean in regions of extensive wind glaze.

Conventional versus reference-surface mass balance

Abstract: Glacier surface mass balance evaluated over the actual glacier geometry depends not only on climatic variations, but also on the dynamic adjustment of glacier geometry. Therefore, it has been proposed that reference-surface balances calculated over a constant glacier hypsometry are better suited for climatic interpretation. Here we present a comparison of 82 year modelled time series (1926-2008) of conventional and reference-surface balance for 36 Swiss glaciers. Over this time period the investigated glaciers have lost 22% of their area, and ice surface elevation close to the current glacier terminus has decreased by 78 m on average. Conventional balance in the last decade, at -0.91 m w.e.a(-1), is 0.14 m w.e.a(-1) less negative than the reference-surface balance. About half of the negative (stabilizing) feedback on mass balance due to glacier terminus retreat is compensated by more negative mass balances due to surface lowering. Short-term climatic variability is clearly reflected in the conventional mass-balance series; however, the magnitude of the long-term negative trend is underestimated compared to that found in the reference-surface balance series. Both conventional and reference-surface specific balances show large spatial variability among the 36 glaciers.

Relative contribution of surface mass balance and ice flux changes to the accelerated thinning of the Mer de Glace (Alps) over 1979-2008

Abstract: By subtracting surface topographies from 1979, 1994, 2000 and 2008, we measured ice thinning rates increasing from 1 m a-1 (1979-1994) to more than 4 m a-1 (2000-2008) on the tongue of the Mer de Glace The relative contributions of changes in surface mass balance and ice fluxes to this acceleration in the thinning are estimated using field and remote sensing measurements Between 1979-1994 and 2000-2008, surface mass balance diminished by 12 m we a-1 mainly because of atmospheric warming Mass balance changes induced by the growing debris-covered area and the evolving glacier hypsometry compensated each other Meanwhile, the Mer de Glace slowed down and the ice fluxes through two cross-sections at 2200 m asl and 2050 m asl decreased by 60% Between 1979-1994 and 2000-2008, two thirds of the increase in the thinning rates was caused by reduced ice fluxes and one third by rising surface ablation However, those numbers need to be interpreted cautiously given our inability to respect mass conservation below our upper cross section An important implication is that large errors would occur if one term of the continuity equation (eg, surface mass balance) was deduced from the two others (eg, elevation and ice flux changes)

Estimating the long-term calving flux of Kronebreen, Svalbard, from geodetic elevation changes and mass-balance modelling

Abstract: 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 tertiary creep of polycrystalline ice: experimental evidence for stress-dependent levels of strain-rate enhancement

Abstract: Laboratory creep deformation experiments have been conducted on initially isotropic laboratory-made samples of polycrystalline ice. Steady-state tertiary creep rates, ∈ · ter , were determined at strains exceeding 10% in either uniaxial-compression or simple-shear experiments. Isotropic minimum strain rates, ∈ · min , determined at ∼1% strain, provide a reference for comparing the relative magnitude of tertiary creep rates in shear and compression through the use of strain-rate enhancement factors, E , defined as the ratio of corresponding tertiary and isotropic minimum creep rates, i.e. E = ∈ · ter / ∈ · min . The magnitude of strain-rate enhancement in simple shear was found to exceed that in uniaxial compression by a constant factor of 2.3. Results of experiments conducted at octahedral shear stresses of τ o = 0.04–0.80 MPa indicate a creep power-law stress exponent of n = 3 for isotropic minimum creep rates and n = 3.5 for tertiary creep rates. The difference in stress exponents for minimum and tertiary creep regimes can be interpreted as a τ o stress-dependent level of strain-rate enhancement, i.e. E α τ 1/2 o . The implications of these results for deformation in complex multicomponent stress configurations and at stresses below those used in the current experiments are discussed.

Signatures of supercooling: McMurdo Sound platelet ice

Abstract: 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.

Seasonal variations of 17 O-excess and d-excess in snow precipitation at Vostok station, East Antarctica

Abstract: The use of water isotopes in polar regions is essential for reconstructing past climate over glacial-interglacial cycles. In addition to d Do rd 18 O, linearly related to condensation temperature, the second-order parameters, d-excess and 17 O-excess, provide important information on the climatic conditions of the source of precipitations. In order to best interpret the glacial-interglacial records of d-excess and 17 O-excess in polar ice cores, it is important to document their present variability, especially in remote and cold regions of East Antarctica. Indeed, the current climatic conditions encountered in these regions provide a good analogy with glacial climatic conditions in a large part of Antarctica. Here we present the first seasonal variations of 17 O-excess and d-excess at Vostok station on an event basis (i.e. samples were collected immediately after each precipitation event) over 1 year. These records show strong correlation between 17 O-excess and d 18 O over the course of the year, with an amplitude 40 per meg (10 -3 %) in the 17 O-excess seasonal cycle, and strong anticorrelation between d-excess and d 18 O, with d-excess variations up to 20%. The d-excess and 17 O-excess variations can be explained by the influence of kinetic fractionation at very low temperatures. The comparison with simple isotopic models confirms this explanation, but cannot explain the link between 17 O-excess, d-excess and temperature without (1) a particular relationship between condensation and surface temperature and/or (2) seasonal changes in the climatic conditions of the source regions.

New shortwave infrared albedo measurements for snow specific surface area retrieval

Abstract: Snow grain-size characterization, its vertical and temporal evolution is a key parameter for the improvement and validation of snow and radiative transfer models (optical and microwave) as well as for remote-sensing retrieval methods. We describe two optical methods, one active and one passive shortwave infrared, for field determination of the specific surface area (SSA) of snow grains. We present a new shortwave infrared (SWIR) camera approach. This new method is compared with a SWIR laser- based system measuring snow albedo with an integrating sphere (InfraRed Integrating Sphere (IRIS)). Good accuracy (10%) and reproducibility in SSA measurements are obtained using the IRIS system on snow samples having densities greater than 200 kg m-3, validated against X-ray microtomography measurements. The SWIRcam approach shows improved sensitivity to snow SSA when compared to a near-infrared camera, giving a better contrast of the snow stratigraphy in a snow pit.

Iceberg signatures and detection in SAR images in two test regions of the Weddell Sea, Antarctica

Abstract: A pixel-based methodology has been established for automatic identification of icebergs in satellite synthetic aperture radar (SAR) images acquired during different seasons and for different seaice conditions. This includes, in particular, smaller icebergs (longitudinal axis 100m to 18.5 km). Investigations were carried out for two test regions located in theWeddell Sea, Antarctica, using images of the Envisat Advanced SAR (ASAR) at HH polarization and of the European Remote-sensing Satellite-2 (ERS-2) SAR (VV-polarized). From the former, a sequence of Image Mode and Wide Swath Mode data are available for the whole of 2006. The ERS data were acquired around the tip of the Antarctic Peninsula in spring and summer months of the years 2000–03. The minimum size of icebergs that could be identified in the IM images was <0.02km2. Radar backscattering coefficients of icebergs, sea ice and open water were determined separately.We demonstrate that the error in separating icebergs from their surroundings (sea ice or open water) depends on meteorological, oceanographic and sea-ice conditions. Also the pre-processing of the SAR images (e.g. speckle reduction) influences iceberg recognition. Differences in detection accuracy as a function of season could not be substantiated for our test sites, but have in general to be taken into account, as results of other investigations indicate.

Volume calculation and analysis of the changes in moraine-dammed lakes in the north Himalaya: a case study of Longbasaba lake

Abstract: Glacial lake outburst flood hazards in the Himalayan region have received considerable attention in recent years. Accurate volume estimation for glacial lakes is important for calculating outburst flood peak discharge and simulating flood evolution. Longbasaba lake, a potentially dangerous moraine-dammed lake, is located on the north side of the Himalaya. Its depth was surveyed using the SyQwest Hydrobox™ high-resolution echo sounder, and 6916 measurements were collected in September 2009. The maximum and average depths of the lake were 102 ± 2 and 48 ± 2 m, respectively. The morphology of the lake basin was modeled by constructing a triangulated irregular network, and the lake was found to have a storage capacity of 0.064 ± 0.002 km 3 . Multi-source remote-sensing images from Landsat MSS, Landsat TM/ETM+ and Terra ASTER and a topographic map were digitized to delineate the outlines of the lake between 1977 and 2009. The results indicate that the length and area of the lake have increased during the past 32 years, with a drastic expansion occurring since 2000. Based on volume and area data of Longbasaba lake in different periods, we deduced an empirical equation of the lake volume―area relationship that can be used to calculate the storage capacity of similar moraine-dammed lakes in the Himalayan region.

Variational assimilation of albedo in a snowpack model and reconstruction of the spatial mass-balance distribution of an alpine glacier

Abstract: 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.

The effect of discrete recharge by moulins and heterogeneity in flow-path efficiency at glacier beds on subglacial hydrology

Abstract: Subglacial conduit systems are thought to consist of dendritic networks that exist at lower pressure than distributed systems and have locations that are determined by theoretical hydraulic potential. On glaciers with moulins, however, meltwater is delivered to glacier beds at discrete points, violating assumptions of uniform recharge needed to calculate potential. To understand how moulins affect subglacial conduit hydrology, we used speleological techniques to map 0.4 km of subglacial conduit at the base of a moulin in Hansbreen, Svalbard, and compared our observations with theoretical predictions. The conduit began in an area predicted to lack drainage, crossed equipotential contours at oblique rather than right angles and was locally anastomotic rather than dendritic. We propose moulin locations, which are determined by the locations of supraglacial streams and crevasses, control locations of subglacial recharge. Because conduits have no direct causal relationship with gradients in effective pressure, this recharge can form conduits in areas of glacier beds that may not be predicted by hydraulic potential theory to have conduits. Recharge by moulins allows hydraulic head to increase in conduits faster and to higher values than in adjacent distributed systems, resulting in an increase rather than a decrease in glacier sliding speeds above subglacial conduits.

Glacier mass loss induced the rapid growth of Linggo Co on the central Tibetan Plateau

Abstract: 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.2 m. The lake volume of Linggo Co increased at a rate of 0.02 � 10 6 , 42.67 � 10 6 and 65.8 � 10 6 m 3 a -1 during the periods 1974-92, 1992-99 and 1999-2010, respectively. Other non- glacier-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) � 10 8 m 3 . The results indicate that glacier mass loss had a significant impact on the growth of Linggo Co over the past 40 years.

Dynamics of tidewater surge-type glaciers in northwest Svalbard

Abstract: 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.