Christopher Nuth

Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas

Glaciers are among the best indicators of terrestrial climate variability, contribute importantly to water resources in many mountainous regions and are a major contributor to global sea level rise. In the Hindu Kush–Karakoram–Himalaya region (HKKH), a paucity of appropriate glacier data has prevented a comprehensive assessment of current regional mass balance. There is, however, indirect evidence of a complex pattern of glacial responses in reaction to heterogeneous climate change signals. Here we use satellite laser altimetry and a global elevation model to show widespread glacier wastage in the eastern, central and south-western parts of the HKKH during 2003–08. Maximal regional thinning rates were 0.66 ± 0.09 metres per year in the Jammu–Kashmir region. Conversely, in the Karakoram, glaciers thinned only slightly by a few centimetres per year. Contrary to expectations, regionally averaged thinning rates under debris-mantled ice were similar to those of clean ice despite insulation by debris covers. The 2003–08 specific mass balance for our entire HKKH study region was −0.21 ± 0.05 m yr−1 water equivalent, significantly less negative than the estimated global average for glaciers and ice caps. This difference is mainly an effect of the balanced glacier mass budget in the Karakoram. The HKKH sea level contribution amounts to one per cent of the present-day sea level rise. Our 2003–08 mass budget of −12.8 ± 3.5 gigatonnes (Gt) per year is more negative than recent satellite-gravimetry-based estimates of −5 ± 3 Gt yr−1 over 2003–10 (ref. 12). For the mountain catchments of the Indus and Ganges basins, the glacier imbalance contributed about 3.5% and about 2.0%, respectively, to the annual average river discharge, and up to 10% for the Upper Indus basin.

On the accuracy of glacier outlines derived from remote-sensing data

Abstract Deriving glacier outlines from satellite data has become increasingly popular in the past decade. In particular when glacier outlines are used as a base for change assessment, it is important to know how accurate they are. Calculating the accuracy correctly is challenging, as appropriate reference data (e.g. from higher-resolution sensors) are seldom available. Moreover, after the required manual correction of the raw outlines (e.g. for debris cover), such a comparison would only reveal the accuracy of the analyst rather than of the algorithm applied. Here we compare outlines for clean and debris-covered glaciers, as derived from single and multiple digitizing by different or the same analysts on very high- (1 m) and medium-resolution (30 m) remote-sensing data, against each other and to glacier outlines derived from automated classification of Landsat Thematic Mapper data. Results show a high variability in the interpretation of debris-covered glacier parts, largely independent of the spatial resolution (area differences were up to 30%), and an overall good agreement for clean ice with sufficient contrast to the surrounding terrain (differences ∼5%). The differences of the automatically derived outlines from a reference value are as small as the standard deviation of the manual digitizations from several analysts. Based on these results, we conclude that automated mapping of clean ice is preferable to manual digitization and recommend using the latter method only for required corrections of incorrectly mapped glacier parts (e.g. debris cover, shadow).

Glacier geometry and elevation changes on Svalbard (1936–90): a baseline dataset

Abstract This study uses older topographic maps made from high-oblique aerial photographs for glacier elevation change studies. We compare the 1936/38 topographic map series of Svalbard (Norwegian Polar Institute) to a modern digital elevation model from 1990. Both systematic and random components of elevation error are examined by analyzing non-glacier elevation difference points. The 1936/38 photographic aerial survey is examined to identify areas with poor data coverage over glaciers. Elevation changes are analyzed for seven regions in Svalbard (~5000 km2), where significant thinning was found at glacier fronts, and elevation increases in the upper parts of the accumulation areas. All regions experience volume losses and negative geodetic balances, although regional variability exists relating to both climate and topography. Many surges are apparent within the elevation change maps. Estimated volume change for the regions is –1.59±0.07km3 a–1 (ice equivalent) for a geodetic annual balance of –0.30ma–1w.e., and the glaciated area has decreased by 16% in the 54 year time interval. The 1936–90 data are compared to modern elevation change estimates in the southern regions, to show that the rate of thinning has increased dramatically since 1990.

Spatial and temporal distribution of mass loss from the Greenland Ice Sheet since AD 1900

The response of the Greenland Ice Sheet (GIS) to changes in temperature during the twentieth century remains contentious, largely owing to difficulties in estimating the spatial and temporal distribution of ice mass changes before 1992, when Greenland-wide observations first became available. The only previous estimates of change during the twentieth century are based on empirical modelling and energy balance modelling. Consequently, no observation-based estimates of the contribution from the GIS to the global-mean sea level budget before 1990 are included in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Here we calculate spatial ice mass loss around the entire GIS from 1900 to the present using aerial imagery from the 1980s. This allows accurate high-resolution mapping of geomorphic features related to the maximum extent of the GIS during the Little Ice Age at the end of the nineteenth century. We estimate the total ice mass loss and its spatial distribution for three periods: 1900–1983 (75.1 ± 29.4 gigatonnes per year), 1983–2003 (73.8 ± 40.5 gigatonnes per year), and 2003–2010 (186.4 ± 18.9 gigatonnes per year). Furthermore, using two surface mass balance models we partition the mass balance into a term for surface mass balance (that is, total precipitation minus total sublimation minus runoff) and a dynamic term. We find that many areas currently undergoing change are identical to those that experienced considerable thinning throughout the twentieth century. We also reveal that the surface mass balance term shows a considerable decrease since 2003, whereas the dynamic term is constant over the past 110 years. Overall, our observation-based findings show that during the twentieth century the GIS contributed at least 25.0 ± 9.4 millimetres of global-mean sea level rise. Our result will help to close the twentieth-century sea level budget, which remains crucial for evaluating the reliability of models used to predict global sea level rise.

Glacier Remote Sensing Using Sentinel-2. Part I: Radiometric and Geometric Performance, and Application to Ice Velocity

With its temporal resolution of 10 days (five days with two satellites, and significantly more at high latitudes), its swath width of 290 km, and its 10 m and 20 m spatial resolution bands from the visible to the shortwave infrared, the European Sentinel-2 satellites have significant potential for glacier remote sensing, in particular mapping of glacier outlines and facies, and velocity measurements. Testing Level 1C commissioning and ramp-up phase data for initial sensor quality experiences, we find a high radiometric performance, but with slight striping effects under certain conditions. Through co-registration of repeat Sentinal-2 data we also find lateral offset patterns and noise on the order of a few metres. Neither of these issues will complicate most typical glaciological applications. Absolute geo-location of the data investigated was on the order of one pixel at the time of writing. The most severe geometric problem stems from vertical errors of the DEM used for ortho-rectifying Sentinel-2 data. These errors propagate into locally varying lateral offsets in the images, up to several pixels with respect to other georeferenced data, or between Sentinel-2 data from different orbits. Finally, we characterize the potential and limitations of tracking glacier flow from repeat Sentinel-2 data using a set of typical glaciers in different environments: Aletsch Glacier, Swiss Alps; Fox Glacier, New Zealand; Jakobshavn Isbree, Greenland; Antarctic Peninsula at the Larsen C ice shelf.

Digital elevation model and orthophotographs of Greenland based on aerial photographs from 1978-1987.

Digital Elevation Models (DEMs) play a prominent role in glaciological studies for the mass balance of glaciers and ice sheets. By providing a time snapshot of glacier geometry, DEMs are crucial for most glacier evolution modelling studies, but are also important for cryospheric modelling in general. We present a historical medium-resolution DEM and orthophotographs that consistently cover the entire surroundings and margins of the Greenland Ice Sheet 1978–1987. About 3,500 aerial photographs of Greenland are combined with field surveyed geodetic ground control to produce a 25 m gridded DEM and a 2 m black-and-white digital orthophotograph. Supporting data consist of a reliability mask and a photo footprint coverage with recording dates. Through one internal and two external validation tests, this DEM shows an accuracy better than 10 m horizontally and 6 m vertically while the precision is better than 4 m. This dataset proved successful for topographical mapping and geodetic mass balance. Other uses include control and calibration of remotely sensed data such as imagery or InSAR velocity maps.

Dynamic perennial firn aquifer on an Arctic glacier

Ice-penetrating radar and GPS observations reveal a perennial firn aquifer (PFA) on a Svalbard ice field, similar to those recently discovered in southeastern Greenland. A bright, widespread radar reflector separates relatively dry and water-saturated firn. This surface, the phreatic firn water table, is deeper beneath local surface elevation maxima, shallower in surface lows, and steeper where the surface is steep. The reflector crosscuts snow stratigraphy; we use the apparent deflection of accumulation layers due to the higher dielectric permittivity below the water table to infer that the firn pore space becomes progressively more saturated as depth increases. Our observations indicate that PFAs respond rapidly (subannually) to surface forcing, and are capable of providing significant input to the englacial hydrology system.

A digital glacier database for svalbard

The archipelago of Svalbard presently contains approximately 33,200 km2 of glaciers, with a large number of small valley glaciers as well as large areas of contiguous icefields and ice caps. While the first glacier inventory was compiled in 1993, there has not been a readily available digital version. Here we present a new digital glacier database, which is being made available through the GLIMS project (Global Land Ice Measurements from Space). Glacier outlines have been created for 1936, 1966–1971, 1990, and 2001–2010. For most glaciers, outlines are available from more than one of these years. Complete coverage of Svalbard is available for the 2001–2010 dataset. Glacier outlines were created using cartographic data from the original Norwegian Polar Institute topographic map series of Svalbard as the basis by delineating individual glaciers and ice streams, assigning unique identification codes relating to hydrological watersheds, digitizing centerlines, and providing a number of attributes for each glacier mask. The 2001–2010 glacier outlines are derived from orthorectified satellite images acquired from SPOT-5 and ASTER satellite sensors. In areas where coverage for all time periods is available, the overwhelming majority of glaciers are observed to be in sustained retreat over the period from 1936 to 2010.

Regional Glacier Mapping Using Optical Satellite Data Time Series

The first of two Sentinel-2 satellites, launched mid2015, has similar characteristics as the Landsat TM/ETM + /OLI satellites. Together, these satellites will produce a tremendous quantity of optical images worldwide for glacier mapping, with increasing temporal coverage toward the more glacierized higher latitudes due to convergence of near-polar orbits. To exploit the potential of such near-future dense time series, methods for mapping glaciers from space should be revisited. Currently, snow and ice are typically classified from an optical satellite image using a multispectral band ratio. For each scene, mapping conditions will vary (e.g., snow, ice, and clouds) and not be equally optimal over the entire scene. The increasing amount of images makes it difficult to manually select the best glacier mapping scene as is the current practice. This work is based on the above robust image ratio method for exploiting the dense temporal image coverage. Four application scenarios using time series of Landsat type data for glacier mapping are presented. First, we synthesize an optimal band ratio image from a stack of images within one season to compensate for regional differences. The second application scenario introduces robust methods to improve automatic glacier mapping by exploiting the seasonal variation in spectral properties of snow. Third, we explore the spatio-temporal variation of glacier surface types. Finally, we show how the synthesized band ratio images from the first application scenario can be used for automatic glacier change detection. In summary, we explore automatic algorithms for glacier mapping applications that exploit the temporal signatures in the satellite data time series.

Decadal Scale Changes in Glacier Area in the Hohe Tauern National Park (Austria) Determined by Object-Based Image Analysis

In this paper, we semi-automatically classify clean and debris-covered ice for 145 glaciers within Hohe Tauern National Park in the Austrian Alps for the years 1985, 2003, and 2013. We also map the end-summer transient snowline (TSL), which approximates the annual Equilibrium Line Altitude (ELA). By comparing our results with the Austrian Glacier Inventories from 1969 and 1998, we calculate a mean reduction in glacier area of 33% between 1969 and 2013. The total ice area reduced at a mean rate of 1.4 km2 per year. This TSL rose by 92 m between 1985 and 2013 to an altitude of 3005 m. Despite some limitations, such as some seasonal snow being present at higher elevations, as well as uncertainties related to the range of years that the LiDAR DEM was collected, our results show that the glaciers within Hohe Tauern National Park conform to the heavy shrinkage experienced in other areas of the European Alps. Moreover, we believe that Object-Based Image Analysis (OBIA) is a promising methodology for future glacier mapping.