Showing papers in "The Cryosphere in 2010"

Overview of areal changes of the ice shelves on the Antarctic Peninsula over the past 50 years

Abstract: . In recent decades, seven out of twelve ice shelves around the Antarctic Peninsula (AP) have either retreated significantly or have been almost entirely lost. At least some of these retreats have been shown to be unusual within the context of the Holocene and have been widely attributed to recent atmospheric and oceanic changes. To date, measurements of the area of ice shelves on the AP have either been approximated, or calculated for individual shelves over dissimilar time intervals. Here we present a new dataset containing up-to-date and consistent area calculations for each of the twelve ice shelves on the AP over the past five decades. The results reveal an overall reduction in total ice-shelf area by over 28 000 km2 since the beginning of the period. Individual ice shelves show different rates of retreat, ranging from slow but progressive retreat to abrupt collapse. We discuss the pertinent features of each ice shelf and also broad spatial and temporal patterns in the timing and rate of retreat. We believe that an understanding of this diversity and what it implies about the underlying dynamics and control will provide the best foundation for developing a reliable predictive skill for ice-shelf change.

A glacier inventory for the western Nyainqentanglha Range and the Nam Co Basin, Tibet, and glacier changes 1976-2009

Abstract: The western Nyainqentanglha Range is located in the south-eastern centre of the Tibetan Plateau. Its north- western slopes drain into Lake Nam Co. The region is of spe- cial interest for glacio-climatological research as it is influ- enced by both the continental climate of Central Asia and the Indian Monsoon system, and situated at the transition zone between temperate and subcontinental glaciers. A glacier in- ventory for the whole mountain range was generated for the year around 2001 using automated remote sensing and GIS techniques based on Landsat ETM+ and SRTM3 DEM data. Glacier change analysis was based on data from Hexagon KH-9 and Landsat MSS (both 1976), Metric Camera (1984), and Landsat TM/ETM+ (1991, 2001, 2005, 2009). Man- ual adjustment was especially necessary for delineating the debris-covered glaciers and the glaciers on the panchromatic Hexagon data. In the years around 2001 the whole moun- tain range contained about 960 glaciers covering an area of 795.6± 22.3 km 2 while the ice in the drainage basin of Nam Co covered 198.1± 5.6 km 2 . The median elevation of the glaciers was about 5800 m with the majority termi- nating around 5600 m. Five glaciers with debris-covered tongues terminated lower than 5200 m. The glacier area de- creased by 6.1± 3% between 1976 and 2001. This is less than reported in previous studies based on the 1970s topo- graphic maps and Landsat data from 2000. Glaciers con- tinued to shrink during the period 2001-2009. No advanc- ing glaciers were detected. Detailed length measurements for five glaciers indicated a retreat of around 10 m per year

Spatial and temporal variability of snow depth and ablation rates in a small mountain catchment

Abstract: . The spatio-temporal variability of the mountain snow cover determines the avalanche danger, snow water storage, permafrost distribution and the local distribution of fauna and flora. Using a new type of terrestrial laser scanner, which is particularly suited for measurements of snow covered surfaces, snow depth was monitored in a high alpine catchment during an ablation period. From these measurements snow water equivalents and ablation rates were calculated. This allowed us for the first time to obtain a high resolution (2.5 m cell size) picture of spatial variability of the snow cover and its temporal development. A very high variability of the snow cover with snow depths between 0–9 m at the end of the accumulation season was observed. This variability decreased during the ablation phase, while the dominant snow deposition features remained intact. The average daily ablation rate was between 15 mm/d snow water equivalent at the beginning of the ablation period and 30 mm/d at the end. The spatial variation of ablation rates increased during the ablation season and could not be explained in a simple manner by geographical or meteorological parameters, which suggests significant lateral energy fluxes contributing to observed melt. It is qualitatively shown that the effect of the lateral energy transport must increase as the fraction of snow free surfaces increases during the ablation period.

Role of glaciers in watershed hydrology: a preliminary study of a "Himalayan catchment"

Abstract: . A large number of Himalayan glacier catchments are under the influence of humid climate with snowfall in winter (November–April) and south-west monsoon in summer (June–September) dominating the regional hydrology. Such catchments are defined as "Himalayan catchment", where the glacier meltwater contributes to the river flow during the period of annual high flows produced by the monsoon. The winter snow dominated Alpine catchments of the Kashmir and Karakoram region and cold-arid regions of the Ladakh mountain range are the other major glacio-hydrological regimes identified in the region. Factors influencing the river flow variations in a "Himalayan catchment" were studied in a micro-scale glacier catchment in the Garhwal Himalaya, covering an area of 77.8 km2. Three hydrometric stations were established at different altitudes along the Din Gad stream and discharge was monitored during the summer ablation period from 1998 to 2004, with an exception in 2002. These data have been analysed along with winter/summer precipitation, temperature and mass balance data of the Dokriani glacier to study the role of glacier and precipitation in determining runoff variations along the stream continuum from the glacier snout to 2360 m a.s.l. The study shows that the inter-annual runoff variation in a "Himalayan catchment" is linked with precipitation rather than mass balance changes of the glacier. This study also indicates that the warming induced an initial increase of glacier runoff and subsequent decline as suggested by the IPCC (2007) is restricted to the glacier degradation-derived component in a precipitation dominant Himalayan catchment and cannot be translated as river flow response. The preliminary assessment suggests that the "Himalayan catchment" could experience higher river flows and positive glacier mass balance regime together in association with strong monsoon. The important role of glaciers in this precipitation dominant system is to augment stream runoff during the years of low summer discharge. This paper intends to highlight the importance of creating credible knowledge on the Himalayan cryospheric processes to develop a more representative global view on river flow response to cryospheric changes and locally sustainable water resources management strategies.

Climate of the Greenland ice sheet using a high-resolution climate model – Part 1: Evaluation

Abstract: . A simulation of 51 years (1957–2008) has been performed over Greenland using the regional atmospheric climate model (RACMO2/GR) at a horizontal grid spacing of 11 km and forced by ECMWF re-analysis products. To better represent processes affecting ice sheet surface mass balance, such as meltwater refreezing and penetration, an additional snow/ice surface module has been developed and implemented into the surface part of the climate model. The temporal evolution and climatology of the model is evaluated with in situ coastal and ice sheet atmospheric measurements of near-surface variables and surface energy balance components. The bias for the near-surface air temperature (−0.8 °C), specific humidity (0.1 g kg−1), wind speed (0.3 m s−1) as well as for radiative (2.5 W m−2 for net radiation) and turbulent heat fluxes shows that the model is in good accordance with available observations on and around the ice sheet. The modelled surface energy budget underestimates the downward longwave radiation and overestimates the sensible heat flux. Due to their compensating effect, the averaged 2 m temperature bias is small and the katabatic wind circulation well captured by the model.

Understanding snow-transport processes shaping the mountain snow-cover

Abstract: . Mountain snow-cover is normally heterogeneously distributed due to wind and precipitation interacting with the snow cover on various scales. The aim of this study was to investigate snow deposition and wind-induced snow-transport processes on different scales and to analyze some major drift events caused by north-west storms during two consecutive accumulation periods. In particular, we distinguish between the individual processes that cause specific drifts using a physically based model approach. Very high resolution wind fields (5 m) were computed with the atmospheric model Advanced Regional Prediction System (ARPS) and used as input for a model of snow-surface processes (Alpine3D) to calculate saltation, suspension and preferential deposition of precipitation. Several flow features during north-west storms were identified with input from a high-density network of permanent and mobile weather stations and indirect estimations of wind directions from snow-surface structures, such as snow dunes and sastrugis. We also used Terrestrial and Airborne Laser Scanning measurements to investigate snow-deposition patterns and to validate the model. The model results suggest that the in-slope deposition patterns, particularly two huge cross-slope cornice-like drifts, developed only when the prevailing wind direction was northwesterly and were formed mainly due to snow redistribution processes (saltation-driven). In contrast, more homogeneous deposition patterns on a ridge scale were formed during the same periods mainly due to preferential deposition of precipitation. The numerical analysis showed that snow-transport processes were sensitive to the changing topography due to the smoothing effect of the snow cover.

The relation between sea ice thickness and freeboard in the Arctic

Abstract: Retrieval of Arctic sea ice thickness from CryoSat-2 radar altimeter freeboard data requires observa- tional data to verify the relation between these two variables. In this study in-situ ice and snow data from 689 observation sites, obtained during the Sever expeditions in the 1980s, have been used to establish an empirical relation between thickness and freeboard of FY ice in late winter. Estimates of mean and variability of snow depth, snow density and ice density were produced on the basis of many field observa- tions. These estimates have been used in the hydrostatic equi- librium equation to retrieve ice thickness as a function of ice freeboard, snow depth and snow/ice density. The accuracy of the ice thickness retrieval has been calculated from the es- timated variability in ice and snow parameters and error of ice freeboard measurements. It is found that uncertainties of ice density and freeboard are the major sources of error in ice thickness calculation. For FY ice, retrieval of 1.0 m (2.0 m) thickness has an uncertainty of 46% (37%), and for MY ice, retrieval of 2.4 m (3.0 m) thickness has an uncertainty of 20% (18%), assuming that the freeboard error is ± 0.03 m for both ice types. For MY ice the main uncertainty is ice density er- ror, since the freeboard error is relatively smaller than that for FY ice. If the freeboard error can be reduced to 0.01 m by av- eraging measurements from CryoSat-2, the error in thickness retrieval is reduced to about 32% for a 1.0 m thick FY floe and to about 18% for a 2.4 m thick MY floe. The remaining error is dominated by uncertainty in ice density. Provision of improved ice density data is therefore important for accurate retrieval of ice thickness from CryoSat-2 data.

Reanalysis of multi-temporal aerial images of Storglaciären, Sweden (1959–99) – Part 1: Determination of length, area, and volume changes

Abstract: Storglaci¨ aren, located in the Kebnekaise massif in northern Sweden, has a long history of glaciological re- search. Early photo documentations date back to the late 19th century. Measurements of front position variations and distributed mass balance have been carried out since 1910 and 1945/46, respectively. In addition to these in-situ mea- surements, aerial photographs have been taken at decadal intervals since the beginning of the mass balance monitor- ing program and were used to produce topographic glacier maps. Inaccuracies in the maps were a challenge to early attempts to derive glacier volume changes and resulted in major differences when compared to the direct glaciological mass balances. In this study, we reanalyzed dia-positives of the original aerial photographs of 1959, -69, -80, -90 and -99 based on consistent photogrammetric processing. From the resulting digital elevation models and orthophotos, changes in length, area, and volume of Storglaci ¨ aren were computed

Geometric changes and mass balance of the Austfonna ice cap, Svalbard

Abstract: . The dynamics and mass balance regime of the Austfonna ice cap, the largest glacier on Svalbard, deviates significantly from most other glaciers in the region and is not fully understood. We have compared ICESat laser altimetry, airborne laser altimetry, GNSS surface profiles and radio echo-sounding data to estimate elevation change rates for the periods 1983–2007 and 2002–2008. The data sets indicate a pronounced interior thickening of up to 0.5 m y−1, at the same time as the margins are thinning at a rate of 1–3 m y−1. The southern basins are thickening at a higher rate than the northern basins due to a higher accumulation rate. The overall volume change in the 2002–2008 period is estimated to be −1.3±0.5 km3 w.e. y−1 (or −0.16±0.06 m w.e. y−1) where the entire net loss is due to a rapid retreat of the calving fronts. Since most of the marine ice loss occurs below sea level, Austfonna's current contribution to sea level change is close to zero. The geodetic results are compared to in-situ mass balance measurements which indicate that the 2004–2008 surface net mass balance has been slightly positive (0.05 m w.e. y−1) though with large annual variations. Similarities between local net mass balances and local elevation changes indicate that most of the ice cap is slow-moving and not in dynamic equilibrium with the current climate. More knowledge is needed about century-scale dynamic processes in order to predict the future evolution of Austfonna based on climate scenarios.

On the potential of very high-resolution repeat DEMs in glacial and periglacial environments

Abstract: . The potential of high-resolution repeat DEMs was investigated for glaciological applications including periglacial features (e.g. rock glaciers). It was shown that glacier boundaries can be delineated using airborne LIDAR-DEMs as a primary data source and that information on debris cover extent could be extracted using multi-temporal DEMs. Problems and limitations are discussed, and accuracies quantified. Absolute deviations of airborne laser scanning (ALS) derived glacier boundaries from ground-truthed ones were below 4 m for 80% of the ground-truthed values. Overall, we estimated an accuracy of +/−1.5% of the glacier area for glaciers larger than 1 km2. The errors in the case of smaller glaciers did not exceed +/−5% of the glacier area. The use of repeat DEMs in order to obtain information on the extent, characteristics and activity of rock glaciers was investigated and discussed based on examples.

Brief Communication: Ikaite (CaCO 3 ·6H 2 O) discovered in Arctic sea ice

Abstract: We report for the first time on the discovery of calcium carbonate crystals as ikaite (CaCO3·6H2O) in sea ice from the Arctic (Kongsfjorden, Svalbard) as confirmed by morphology and indirectly by X-ray diffraction as well as XANES spectroscopy of its amorophous decomposition product. This finding demonstrates that the precipitation of calcium carbonate during the freezing of sea ice is not re- stricted to the Antarctic, where it was observed for the first time in 2008. This observation is an important step in the quest to quantify its impact on the sea ice driven carbon cy- cle.

Climate of the Greenland ice sheet using a high-resolution climate model – Part 2: Near-surface climate and energy balance

Abstract: . The spatial variability of near-surface variables and surface energy balance components over the Greenland ice sheet are presented, using the output of a regional atmospheric climate model for the period 1958–2008. The model was evaluated in Part 1 of this paper. The near-surface temperature over the ice sheet is affected by surface elevation, latitude, longitude, large-scale and small-scale advection, occurrence of summer melt and mesoscale topographical features. The atmospheric boundary layer is characterised by a strong temperature inversion, due to continuous longwave cooling of the surface. In combination with a gently sloping surface the radiative loss maintains a persistent katabatic wind. This radiative heat loss is mainly balanced by turbulent sensible heat transport towards the surface. In summer, the surface is near radiative balance, resulting in lower wind speeds. Absorption of shortwave radiation and a positive subsurface heat flux due to refreezing melt water are heat sources for surface sublimation and melt. The strongest temperature deficits (>13 °C) are found on the northeastern slopes, where the strongest katabatic winds (>9 m s−1) and lowest relative humidity (

A sea-ice thickness retrieval model for 1.4 GHz radiometry and application to airborne measurements over low salinity sea-ice

Abstract: . In preparation for the European Space Agency's (ESA) Soil Moisture and Ocean Salinity (SMOS) mission, we investigated the potential of L-band (1.4 GHz) radiometry to measure sea-ice thickness. Sea-ice brightness temperature was measured at 1.4 GHz and ice thickness was measured along nearly coincident flight tracks during the SMOS Sea-Ice campaign in the Bay of Bothnia in March 2007. A research aircraft was equipped with the L-band Radiometer EMIRAD and coordinated with helicopter based electromagnetic induction (EM) ice thickness measurements. We developed a three layer (ocean-ice-atmosphere) dielectric slab model for the calculation of ice thickness from brightness temperature. The dielectric properties depend on the relative brine volume which is a function of the bulk ice salinity and temperature. The model calculations suggest a thickness sensitivity of up to 1.5 m for low-salinity (multi-year or brackish) sea-ice. For Arctic first year ice the modelled thickness sensitivity is less than half a meter. It reduces to a few centimeters for temperatures approaching the melting point. The campaign was conducted under unfavorable melting conditions and the spatial overlap between the L-band and EM-measurements was relatively small. Despite these disadvantageous conditions we demonstrate the possibility to measure the sea-ice thickness with the certain limitation up to 1.5 m. The ice thickness derived from SMOS measurements would be complementary to ESA's CryoSat-2 mission in terms of the error characteristics and the spatiotemporal coverage. The relative error for the SMOS ice thickness retrieval is expected to be not less than about 20%.

The sea level fingerprint of recent ice mass fluxes

Abstract: The sea level contribution from glacial sources has been accelerating during the first decade of the 21st Century (Meier et al., 2007; Velicogna, 2009). This contribution is not distributed uniformly across the world's oceans due to both oceanographic and gravitational effects. We compute the sea level signature for ice mass fluxes due to changes in the gravity field, Earth's rotation and related effects for the nine year period 2000–2008. Mass loss from Greenland results in a relative sea level (RSL) reduction for much of North Western Europe and Eastern Canada. RSL rise from this source is concentrated around South America. Losses in West Antarctica marginally compensate for this and produce maxima along the coastlines of North America, Australia and Oceania. The combined far-field pattern of wastage from all ice melt sources, is dominated by losses from the ice sheets and results in maxima at latitudes between 20° N and 40° S across the Pacific and Indian Oceans, affecting particularly vulnerable land masses in Oceania. The spatial pattern of RSL variations from ice mass losses used in this study is time-invariant and cumulative. Thus, sea level rise, based on the gravitational effects from the ice losses considered here, will be amplified for this sensitive region.

Investigating the sensitivity of numerical model simulations of the modern state of the Greenland ice-sheet and its future response to climate change

Abstract: Ice thickness and bedrock topography are essential boundary conditions for numerical modelling of the evolution of the Greenland ice-sheet (GrIS). The datasets currently in use by the majority of GrIS modelling studies are over two decades old and based on data collected from the 1970s and 80s. We use a newer, high-resolution Digital Elevation Model of the GrIS and new temperature and precipitation forcings to drive the Glimmer ice-sheet model offline under steady state, present day climatic conditions. Comparisons are made of ice-sheet geometry between these new datasets and older ones used in the EISMINT-3 exercise. We find that changing to the newer bedrock and ice thickness makes the greatest difference to Greenland ice volume and ice surface extent. When all boundary conditions and forcings are simultaneously changed to the newer datasets the ice-sheet is 33 larger in volume compared with observation and 17 larger than that modelled by EISMINT-3. We performed a tuning exercise to improve the modelled present day ice-sheet. Several solutions were chosen in order to represent improvement in different aspects of the GrIS geometry: ice thickness, ice volume and ice surface extent. We applied these new parameter sets for Glimmer to several future climate scenarios where atmospheric CO2 concentration was elevated to 400, 560 and 1120 ppmv (compared with 280 ppmv in the control) using a fully coupled General Circulation Model. Collapse of the ice-sheet was found to occur between 400 and 560 ppmv, a threshold substantially lower than previously modelled using the standard EISMINT-3 setup. This work highlights the need to assess carefully boundary conditions and forcings required by ice-sheet models, particularly in terms of the abstractions required for large-scale ice-sheet models, and the implications that these can have on predictions of ice-sheet geometry under past and future climate scenarios. © 2010 Author(s).

An explanation for the dark region in the western melt zone of the Greenland ice sheet

Abstract: . The western part of the Greenland ice sheet contains a region that is darker than the surrounding ice. This feature has been analysed with the help of MODIS images. The dark region appears every year during the summer season and can always be found at the same location, which makes meltwater unlikely as the only source for the low albedos. Spectral information indicates that the ice in this region contains more debris than the ice closer to the margin. ASTER images reveal a wavy pattern in the darker ice. Based on these findings we conclude that ice, containing dust from older periods, is presently outcropping near the margin, leading to albedos lower than observed for the remaining ablation area. Therefore it can be concluded that the accumulation of meltwater is a result rather than a cause of the darkening.

Reanalysis of multi-temporal aerial images of Storglaciären, Sweden (1959–99) – Part 2: Comparison of glaciological and volumetric mass balances

Abstract: Seasonal glaciological mass balances have been measured on Storglaciwithout interruption since 1945/46. In addition, aerial surveys have been carried out on a decadal basis since the beginning of the observation pro- gram. Early studies had used the resulting aerial photographs to produce topographic glacier maps with which the in-situ observations could be verified. However, these maps as well as the derived volume changes are subject to errors which resulted in major differences between the derived volumet- ric and the glaciological mass balance. As a consequence, the original photographs were re-processed using uniform photogrammetric methods, which resulted in new volumetric mass balances for 1959-69, 1969-80, 1980-90, and 1990- 99. We compared these new volumetric mass balances with mass balances obtained by standard glaciological methods including an uncertainty assessment considering all related previous studies. The absolute differences between volumet- ric and the glaciological mass balances are 0.8 m w.e. for the period of 1959-69 and 0.3 m w.e. or less for the other survey periods. These deviations are slightly reduced when consid- ering corrections for systematic uncertainties due to differ- ences in survey dates, reference areas, and internal ablation, whereas internal accumulation systematically increases the mismatch. However, the mean annual differences between glaciological and volumetric mass balance are less than the uncertainty of the in-situ stake reading and stochastic error bars of both data series overlap. Hence, no adjustment of the glaciological data series to the volumetric one is required.

An efficient regional energy-moisture balance model for simulation of the Greenland Ice Sheet response to climate change

Abstract: . In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS.

Monitoring ice shelf velocities from repeat MODIS and Landsat data - a method study on the Larsen C ice shelf, Antarctic Peninsula, and 10 other ice shelves around Antarctica

Abstract: . We investigate the velocity field of the Larsen C ice shelf, Antarctic Peninsula, over the periods 2002–2006 and 2006–2009 based on repeat optical satellite data. The velocity field of the entire ice shelf is measured using repeat low resolution MODIS data (250 m spatial resolution). The measurements are validated for two ice shelf sections against repeat medium resolution Landsat 7 ETM + pan data (15 m spatial resolution). Horizontal surface velocities are obtained through image matching using both orientation correlation operated in the frequency domain and normalized crosscorrelation operated in the spatial domain, and the two methods compared. The uncertainty in the displacement measurements turns out to be about one fourth of the pixel size for the MODIS derived data, and about one pixel for the Landsat derived data. The difference between MODIS and Landsat based speeds is −15.4 m a−1 and 13.0 m a−1, respectively, for the first period for the two different validation sections on the ice shelf, and −26.7 m a−1 and 27.9 m a−1 for the second period for the same sections. This leads us to conclude that repeat MODIS images are well suited to measure ice shelf velocity fields and monitor their changes over time. Orientation correlation seems better suited for this purpose because it produces fewer mismatches, is able to match images with regular noise and data voids, and is faster. Since it can match images with regular data voids it is possible to match Landsat 7 ETM+ images even after the 2003 failure of the Scan Line Corrector (SLC off) that leaves significant image stripes with no data. Image matching based on the original 12-bit radiometric resolution MODIS data produced slightly better results than using the 8-bit version of the same images. Streamline interpolation from the obtained surface velocity field on Larsen C indicates ice travel times of up to 450 to 550 years between the inland boundary and the ice shelf edge. In a second step of the study we test our method successfully on 10 other ice shelves around Antarctica demonstrating that the approach presented could in fact be used for large scale monitoring of ice shelf dynamics.

Forecasting temperate alpine glacier survival from accumulation zone observations

Abstract: . Temperate alpine glacier survival is dependent on the consistent presence of an accumulation zone. Frequent low accumulation area ratio values, below 30%, indicate the lack of a consistent accumulation zone, which leads to substantial thinning of the glacier in the accumulation zone. This thinning is often evident from substantial marginal recession, emergence of new rock outcrops and surface elevation decline in the accumulation zone. In the North Cascades 9 of the 12 examined glaciers exhibit characteristics of substantial accumulation zone thinning; marginal recession or emergent bedrock areas in the accumulation zone. The longitudinal profile thinning factor, f, which is a measure of the ratio of thinning in the accumulation zone to that at the terminus, is above 0.6 for all glaciers exhibiting accumulation zone thinning characteristics. The ratio of accumulation zone thinning to cumulative mass balance is above 0.5 for glacier experiencing substantial accumulation zone thinning. Without a consistent accumulation zone these glaciers are forecast not to survive the current climate or future additional warming. The results vary considerably with adjacent glaciers having a different survival forecast. This emphasizes the danger of extrapolating survival from one glacier to the next.

Modelling snowdrift sublimation on an Antarctic ice shelf

Abstract: In this paper, we estimate the contribution of snowdrift sublimation (SUds) to the surface mass balance at Neumayer, located on the Ekstr¨ om ice shelf in Eastern Antarctica. A single column version of the RACMO2-ANT model is used as a physical interpolation tool of high-quality radiosonde and surface measurements for a 15-yr period (1993-2007), and combined with a routine to calculate snow- drift sublimation and horizontal snow transport. The site is characterised by a relatively mild, wet and windy climate, so snowdrift is a common phenomenon. The modelled tim- ing and frequency of snowdrift events compares well with observations. This is further illustrated by an additional sim- ulation for Kohnen base, where the timing of snowdrift is re- alistic, although the modelled horizontal transport is overes- timated. Snowdrift sublimation is mainly dependent on wind speed, but also on relative humidity and temperature. Dur- ing high wind speeds, SUds saturates and cools the air, limit- ing its own strength. We estimate that SUds removes around 16%±8% of the accumulated snow from the surface. The total sublimation more than triples when snowdrift is con- sidered, although snowdrift sublimation limits sublimation at the surface. SUds shows a strong seasonal cycle, as well as large inter-annual variability. This variability can be related to the variability of the atmospheric conditions in the surface layer.

Monitoring of active layer dynamics at a permafrost site on Svalbard using multi-channel ground-penetrating radar

Abstract: . Multi-channel ground-penetrating radar is used to investigate the late-summer evolution of the thaw depth and the average soil water content of the thawed active layer at a high-arctic continuous permafrost site on Svalbard, Norway. Between mid of August and mid of September 2008, five surveys have been conducted in gravelly soil over transect lengths of 130 and 175 m each. The maximum thaw depths range from 1.6 m to 2.0 m, so that they are among the deepest thaw depths recorded in sediments on Svalbard so far. The thaw depths increase by approximately 0.2 m between mid of August and beginning of September and subsequently remain constant until mid of September. The thaw rates are approximately constant over the entire length of the transects within the measurement accuracy of about 5 to 10 cm. The average volumetric soil water content of the thawed soil varies between 0.18 and 0.27 along the investigated transects. While the measurements do not show significant changes in soil water content over the first four weeks of the study, strong precipitation causes an increase in average soil water content of up to 0.04 during the last week. These values are in good agreement with evapotranspiration and precipitation rates measured in the vicinity of the the study site. While we cannot provide conclusive reasons for the detected spatial variability of the thaw depth at the study site, our measurements show that thaw depth and average soil water content are not directly correlated. The study demonstrates the potential of multi-channel ground-penetrating radar for mapping thaw depth in permafrost areas. The novel non-invasive technique is particularly useful when the thaw depth exceeds 1.5 m, so that it is hardly accessible by manual probing. In addition, multi-channel ground-penetrating radar holds potential for mapping the latent heat content of the active layer and for estimating weekly to monthly averages of the ground heat flux during the thaw period.

Parameterising the grounding line in flow-line ice sheet models

Abstract: . Realistic predictions of the behaviour of marine ice sheets require that models are able to robustly simulate grounding line migration. Fixed-grid ice sheet models have been shown to exhibit inconsistent and hence unreliable grounding line migration, except at very high resolution not yet achievable in whole ice sheet simulations. In this study we present several different approaches to parameterising the grounding line. These are distinguished by choices regarding the ice thickness profile from the last grounded to the first floating grid point, and how this profile impacts the gravitational driving stress and basal drag. We demonstrate that the most obvious choice of thickness parameterisation, linear interpolation from the last grounded to the first floating grid point, is not the most effective. We show that use of a grounding line parameterisation greatly improves performance, and that choice of a better grounding line parameterisation over a simpler one can bring further improvements, in terms of both accuracy and self consistent behaviour, comparable to increasing the grid resolution by factor two (i.e. doubling the number of grid points). The approach presented here to parameterising the grounding line does not completely solve the grounding line problem, however it reduces the resolution required. The parameterisations are presented in the context of a one dimensional "shelfy-stream" flow-line model, but could be extended to cope with more than one dimension and other model formulations.

Time-lapse refraction seismic tomography for the detection of ground ice degradation

Abstract: . The ice content of the subsurface is a major factor controlling the natural hazard potential of permafrost degradation in alpine terrain. Monitoring of changes in ice content is therefore similarly important as temperature monitoring in mountain permafrost. Although electrical resistivity tomography monitoring (ERTM) proved to be a valuable tool for the observation of ice degradation, results are often ambiguous or contaminated by inversion artefacts. In theory, the sensitivity of P-wave velocity of seismic waves to phase changes between unfrozen water and ice is similar to the sensitivity of electric resistivity. Provided that the general conditions (lithology, stratigraphy, state of weathering, pore space) remain unchanged over the observation period, temporal changes in the observed travel times of repeated seismic measurements should indicate changes in the ice and water content within the pores and fractures of the subsurface material. In this paper, a time-lapse refraction seismic tomography (TLST) approach is applied as an independent method to ERTM at two test sites in the Swiss Alps. The approach was tested and validated based on a) the comparison of time-lapse seismograms and analysis of reproducibility of the seismic signal, b) the analysis of time-lapse travel time curves with respect to shifts in travel times and changes in P-wave velocities, and c) the comparison of inverted tomograms including the quantification of velocity changes. Results show a high potential of the TLST approach concerning the detection of altered subsurface conditions caused by freezing and thawing processes. For velocity changes on the order of 3000 m/s even an unambiguous identification of significant ice loss is possible.

Multi-channel ground-penetrating radar to explore spatial variations in thaw depth and moisture content in the active layer of a permafrost site

Abstract: . Multi-channel ground-penetrating radar (GPR) was applied at a permafrost site on the Tibetan Plateau to investigate the influence of surface properties and soil texture on the late-summer thaw depth and average soil moisture content of the active layer. Measurements were conducted on an approximately 85 × 60 m2 sized area with surface and soil textural properties that ranged from medium to coarse textured bare soil to finer textured, sparsely vegetated areas covered with fine, wind blown sand, and it included the bed of a gravel road. The survey allowed a clear differentiation of the various units. It showed (i) a shallow thaw depth and low average soil moisture content below the sand-covered, vegetated area, (ii) an intermediate thaw depth and high average soil moisture content along the gravel road, and (iii) an intermediate to deep thaw depth and low to intermediate average soil moisture content in the bare soil terrain. From our measurements, we found hypotheses for the permafrost processes at this site leading to the observed late-summer thaw depth and soil moisture conditions. The study clearly indicates the complicated interactions between surface and subsurface state variables and processes in this environment. Multi-channel GPR is an operational technology to efficiently study such a system at scales varying from a few meters to a few kilometers.

Response of the ice cap Hardangerjøkulen in southern Norway to the 20th and 21st century climates

Abstract: Glaciers respond to mass balance changes by ad- justing their surface elevation and area. These properties in their turn affect the local and area-averaged mass bal- ance. To incorporate this interdependence in the response of glaciers to climate change, models should include an in- teractive scheme coupling mass balance and ice dynamics. In this study, a spatially distributed mass balance model, com- prising surface energy balance calculations, was coupled to a vertically integrated ice-flow model based on the shallow ice approximation. The coupled model was applied to the ice cap Hardangerjokulen in southern Norway. The avail- able glacio-meteorological records, mass balance and glacier length change measurements were utilized for model calibra- tion and validation. Forced with meteorological data from nearby synoptic weather stations, the coupled model real- istically simulated the observed mass balance and glacier length changes during the 20th century. The mean climate for the period 1961-1990, computed from local meteorolog- ical data, was used as a basis to prescribe climate projec- tions for the 21st century at Hardangerjokulen. For a linear temperature increase of 3 C from 1961-1990 to 2071-2100, the modelled net mass balance soon becomes negative at all altitudes and Hardangerjokulen disappears around the year 2100. The projected changes in the other meteorological variables could at most partly compensate for the effect of the projected warming.

Assessing high altitude glacier thickness, volume and area changes using field, GIS and remote sensing techniques: the case of Nevado Coropuna (Peru)

Abstract: . Higher temperatures and changes in precipitation patterns have induced an acute decrease in Andean glaciers, thus leading to additional stress on water supply. To adapt to climate changes, local governments need information on the rate of glacier area and volume losses and on current ice thickness. Remote sensing analyses of Coropuna glacier (Peru) delineate an acute glaciated area decline between 1955 and 2008. We tested how volume changes can be estimated with remote sensing and GIS techniques using digital elevation models derived from both topographic maps and satellite images. Ice thickness was measured in 2004 using a Ground Penetrating Radar coupled with a Ground Positioning System during a field expedition. It provided profiles of ice thickness on different slopes, orientations and altitudes. These were used to model the current glacier volume using Geographical Information System and statistical multiple regression techniques. The results revealed a significant glacier volume loss; however the uncertainty is higher than the measured volume loss. We also provided an estimate of the remaining volume. The field study provided the scientific evidence needed by COPASA, a local Peruvian NGO, and GTZ, the German international cooperation agency, in order to alert local governments and communities and guide them in adopting new climate change adaptation policies.

The influence of changes in glacier extent and surface elevation on modeled mass balance

Abstract: . Glaciers are widely recognized as unique demonstration objects for climate change impacts, mostly due to the strong change of glacier length in response to small climatic changes. However, glacier mass balance as the direct response to the annual atmospheric conditions can be better interpreted in meteorological terms. When the climatic signal is deduced from long-term mass balance data, changes in glacier geometry (i.e. surface extent and elevation) must be considered as such adjustments form an essential part of the glacier reaction to new climatic conditions. In this study, a set of modelling experiments is performed to assess the influence of changes in glacier geometry on mass balance for constant climatic conditions. The calculations are based on a simplified distributed energy/mass balance model in combination with information on glacier extent and surface elevation for the years 1850 and 1973/1985 for about 60 glaciers in the Swiss Alps. The results reveal that over this period about 50–70% of the glacier reaction to climate change (here a one degree increase in temperature) is "hidden" in the geometric adjustment, while only 30–50% can be measured as the long-term mean mass balance. For larger glaciers, the effect of the areal change is partly reduced by a lowered surface elevation, which results in a slightly more negative balance despite a potential increase of topographic shading. In view of several additional reinforcement feedbacks that are observed in periods of strong glacier decline, it seems that the climatic interpretation of long-term mass balance data is rather complex.

Freshwater flux to Sermilik Fjord, SE Greenland

Abstract: . Terrestrial inputs of freshwater flux to Sermilik Fjord, SE Greenland, were estimated, indicating ice discharge to be the dominant source of freshwater. A freshwater flux of 40.4 ± 4.9×109 m3 y−1 was found (1999–2008), with an 85% contribution originated from ice discharge (65% alone from Helheim Glacier), 11% from terrestrial surface runoff (from melt water and rain), 3% from precipitation at the fjord surface area, and 1% from subglacial geothermal and frictional melting due to basal ice motion. The results demonstrate the dominance of ice discharge as a primary mechanism for delivering freshwater to Sermilik Fjord. Time series of ice discharge for Helheim Glacier, Midgard Glacier, and Fenris Glacier were calculated from satellite-derived average surface velocity, glacier width, and estimated ice thickness, and fluctuations in terrestrial surface freshwater runoff were simulated based on observed meteorological data. These simulations were compared and bias corrected against independent glacier catchment runoff observations. Modeled runoff to Sermilik Fjord was variable, ranging from 2.9 ± 0.4×109 m3 y−1 in 1999 to 5.9 ± 0.9×109 m3 y−1 in 2005. The sub-catchment runoff of the Helheim Glacier region accounted for 25% of the total runoff to Sermilik Fjord. The runoff distribution from the different sub-catchments suggested a strong influence from the spatial variation in glacier coverage, indicating high runoff volumes, where glacier cover was present at low elevations.

High resolution modelling of snow transport in complex terrain using downscaled MM5 wind fields

Abstract: Snow transport is one of the most dominant pro- cesses influencing the snow cover accumulation and ablation in high mountain environments. Hence, the spatial and tem- poral variability of the snow cover is significantly modified with respective consequences on the total amount of water in the snow pack, on the temporal dynamics of the runoff and on the energy balance of the surface. For the present study we used the snow transport model SnowModel in combina- tion with MM5 (Penn State University - National Center for Atmospheric Research MM5 model) generated wind fields. In a first step the MM5 wind fields were downscaled by us- ing a semi-empirical approach which accounts for the eleva- tion difference of model and real topography, and vegetation. The target resolution of 30 m corresponds to the resolution of the best available DEM and land cover map of the test site Berchtesgaden National Park. For the numerical mod- elling, data of six automatic meteorological stations were used, comprising the winter season (September-August) of 2003/04 and 2004/05. In addition we had automatic snow depth measurements and periodic manual measurements of snow courses available for the validation of the results. It could be shown that the model performance of SnowModel could be improved by using downscaled MM5 wind fields for the test site. Furthermore, it was shown that an estimation of snow transport from surrounding areas to glaciers becomes possible by using downscaled MM5 wind fields.