Glacier inventory of the Gran Campo Nevado Ice Cap in the Southern Andes and glacier changes observed during recent decades

The Gran Campo Nevado (GCN) forms an isolated ice cap on the Península Muñoz Gamero (PMG) located 200 km to the south of the Southern Patagonia Icefield (SPI). We present a glacier inventory of the GCN made up by 27 drainage basins (in total 199.5 km²) and other small cirque and valley glaciers of the southern part of PMG (in total 53 km²). The glacier inventory is based on a digital elevation model (DEM) and ortho-photos. Contour lines from maps, relief information derived from Landsat TM satellite imagery from 1986 and 2002 and stereoscopic data from aerial photos were combined in a knowledge-based scheme to obtain a DEM of the area. A digital ortho-photo map based on aerial photos from 1998 and several ortho-photos based on aerial photos from 1942 and 1984 could be produced from the initial DEM. A geographical information system (GIS) served to outline the extent of the present glaciation. All major glaciers of the GCN show a significant glacier retreat during the last 60 yr. Some of the outlet glaciers lost more than 20% of their total area during this period. Overall glacier retreat amounts to 2.8% of glacier length per decade and the glacier area loss is 2.4% per decade in the period from 1942 to 2002. We hypothesise that GCN glaciers may have reacted faster and more synchronously with the observed warming trend during recent decades when compared with the SPI.     

Sedimentary effects on the expansion of a Himalayan supraglacial lake

Supraglacial Tsho Rolpa Lake in the Nepal Himalaya has been increasing rapidly in size since the 1950s, corresponding to the mountain-glacier shrinkage after the Little Ice Age. The lake basin expansion results from the subsidence by dead-ice melt below the bottom of the lake, and the retreat of the glacier terminus. Field observations of Tsho Rolpa in 1996 revealed that the retreat of glacier terminus is connected to a wind-induced vertical circulation of surface water heated by solar radiation. In order to clarify the mechanism of the lake expansion associated with sedimentary processes, we measured bottom sedimentation rate with some sediment traps, and vertical suspended sediment concentration (SSC) and water temperature, and analyzed the grain size of suspended and trapped sediments. The sediments, mostly composed of clay-sized grains, are dominantly supplied by glacier-melt water inflow at the glacier terminus. Sedimentary processes of such fine sediment comprise: (1) suspended-sediment fallout from intrusion of horizontal currents; (2) sediment sorting by sediment-laden underflows; and (3) the debris supply from the ice collapse at the glacier terminus. The (1) and (2) processes produce the density stratification of the lake, accompanied by a pycnocline at a depth of about 27 m. The existence of the pycnocline builds up the vertical water circulation in the surface layer to enhance the glacier-melt at the terminus. With respect to the subsidence of the lake bottom, nearly molecular thermal diffusion is probably dominant near the bottom of the deepest point, which results from the kinetic-energy dissipation of sediment-laden underflows. The stable existence of the bottom turbid water throughout the year could cause continuous dead-ice melt below the lake bottom.     

Palaeoecological constraints on late Glacial and Holocene ice retreat in the Southern Andes (53°S)

Late Glacial to Holocene ice retreat was investigated along a 120 km long fjord system, reaching from Gran Campo Nevado (GCN) to Seno Skyring in the southernmost Andes (53°S). The aim was to improve the knowledge on regional and global control on glacier recession with special emphasis on latitudinal shifting of the westerlies. The timing of ice retreat was derived from peat and sediment cores, using mineralogical and chemical characteristics, and pollen as proxies. Stratigraphy was based on ¹⁴C-AMS ages and tephrochronology. The ice retreat of the Seno Skyring Glacier lobe is marked by an ice rafted debris layer which was formed around 18,300 to 17,500 cal. yr B.P. Subsequently, fast glacier retreat occurred until around 15,000 to 14,000 cal. yr B.P. during which around 84% of Skyring Glacier were lost. This fast recession was probably also triggered by an increase of the Equilibrium Line Altitude (ELA) from 200 to 300 m. Subsequently, the ice surface was lowered below the ELA in an area that previously made up more than 50% of the accumulation area. Much slower retreat and glacier fluctuations of limited extent in the fjord channel system northeast of GCN occurred between around 14,000 to 11,000 cal. yr B.P. during both the Antarctic Cold Reversal and the Younger Dryas. This slow down of retreat indicates a decline in the general warming trend and/or increased precipitation, due to a southward migration of the westerlies. After around 11,000 cal. yr B.P. pollen distribution shows evolved Magellanic Rainforest and similar climate as at present, which lasted throughout most of the Holocene. Only Late Neoglacial moraine systems were formed in the period 1220–1460 AD, and subsequently in the 1620s AD, and between 1870 and 1910 AD. The results indicate that the Gran Campo Nevado ice cap has reacted more sensitive and partly distinct to climate change, compared to the Patagonian Ice Field.     

Climate changes and recent glacier behaviour in the Chilean Lake District

Atmospheric temperatures measured at the Chilean Lake District (38°–42°S) showed contrasting trends during the second half of the 20th century. The surface cooling detected at several meteorological stations ranged from − 0.014 to − 0.021 °C a^− 1, whilst upper troposphere (850–300 gpm) records at radiosonde of Puerto Montt (41°26′S/73°07′W) revealed warming between 0.019 and 0.031 °C a^− 1. Regional rainfall data collected from 1961 to 2000 showed the overall decrease with a maximum rate of − 15 mm a^− 2 at Valdivia st. (39°38′S/73°05′W). These ongoing climatic changes, especially the precipitation reduction, seem to be related to El Niño–Southern Oscillation (ENSO) phenomena which has been more frequent after 1976. Glaciers within the Chilean Lake District have significantly retreated during recent decades, in an apparent out-of-phase response to the regional surface cooling. Moreover, very little is known about upper troposphere changes and how they can enhance the glacier responses. In order to analyse their behaviour in the context of the observed climate changes, Casa Pangue glacier (41°08′S/71°52′W) has been selected and studied by comparing Digital Elevation Models (DEMs) computed at three different dates throughout the last four decades. This approach allowed the determination of ice elevation changes between 1961 and 1998, yielding a mean thinning rate of − 2.3 ± 0.6 m a^− 1. Strikingly, when ice thinning is computed for the period between 1981 and 1998, the resulting rate is 50% higher (− 3.6 ± 0.6 m a^− 1). This enhanced trend and the related area loss and frontal retreat suggests that Casa Pangue might currently be suffering negative mass balances in response to the upper troposphere warming and decreased precipitation of the last 25–30 yr, as well as debris cover would not prevent the glacier from a fast reaction to climate forcing. Most of recent glaciological studies regarding Andean glaciers have concentrated on low altitude changes, namely frontal variations, however, in order to better understand the regional glacier changes, new data are necessary, especially from the accumulation areas.     

Recent evolution and mass balance of Cordón Martial glaciers, Cordillera Fueguina Oriental

Past and present glacier changes have been studied at Cordón Martial, Cordillera Fueguina Oriental, Tierra del Fuego, providing novel data for the Holocene deglaciation history of southern South America and extrapolating as well its future behavior based on predicted climatic changes. Regional geomorphologic and stratigraphic correlations indicate that the last glacier advance deposited the ice-proximal (“internal”) moraines of Cordón Martial, around 330 ¹⁴C yr BP, during the Late Little Ice Age (LLIA). Since then glaciers have receded slowly, until 60 years ago, when major glacier retreat started. There is a good correspondence for the past 100 years between the surface area variation of four small cirque glaciers at Cordón Martial and the annual temperature and precipitation data of Ushuaia. Between 1984 and 1998, Martial Este Glacier lost 0.64 ± 0.02 × 10⁶ m³ of ice mass (0.59 ± 0.02 × 10⁶ m³ w.e.), corresponding to an average ice thinning of 7.0 ± 0.2 m (6.4 ± 0.2 m w.e), according to repeated topographic mapping. More detailed climatic data have been obtained since 1998 at the Martial Este Glacier, including air temperature, humidity and solar radiation. These records, together with the monthly mass balance measured since March 2000, document the annual response of the Martial Este Glacier to the climate variation. Mass balances during hydrological years were positive in 2000, negative in 2001 and near equilibrium in 2002. Finally, using these data and the regional temperature trend projections, modeled for different future scenarios by the Atmosphere-Ocean Model (GISS-NASA/GSFC), potential climatic-change effects on this mountain glacier were extrapolated. The analysis shows that only the Martial Este Glacier may survive this century.     

Morphological and ice-dynamical changes on the Tasman Glacier, New Zealand, 1990–2007

This paper presents data concerning recent (1990–2007) surface morphological and ice-dynamical changes on the Tasman Glacier, New Zealand. We use remote-sensing data to derive rates of lake growth, glacier velocities and rates of glacier surface lowering. Between 1990 and 2007, the glacier terminus receded ~ 3.5 km and a large ice-contact proglacial lake developed behind the outwash head. By 2007 the lake area was ~ 6 km² and had replaced the majority of the lowermost 4 km of the glacier tongue. There is evidence that lake growth is proceeding at increasing rates — the lake area doubled between 2000 and 2007 alone. Measured horizontal glacier velocities decline from 150 m a^− 1 in the upper glacier catchment to almost zero at the glacier terminus and there is a consequent down-glacier increase in surface debris cover. Surface debris mapping shows that a large catastrophic rockfall onto the glacier surface in 1991 is still evident as a series of arcuate debris ridges below the Hochstetter icefall. Calculated glacier surface lowering is most clearly pronounced around the terminal area of the glacier tongue, with down-wasting rates of 4.2 ± 1.4 m a^− 1 in areas adjacent to the lateral moraine ridges outside of the current lake extent. Surface lowering rates of approximately 1.9 ± 1.4 m a^− 1 are common in the upper areas of the glacier. Calculations of future lake expansion are dependent on accurate bathymetric and bed topography surveys, but published data indicate that a further 8–10 km of the glacier is susceptible to calving and further lake development in the future.     

Calculation and visualisation of future glacier extent in the Swiss Alps by means of hypsographic modelling

The observed rapid glacier wastage in the European Alps during the past 20 years already has strong impacts on the natural environment (rock fall, lake formation) as well as on human activities (tourism, hydro-power production, etc.) and poses several new challenges also for glacier monitoring. With a further increase of global mean temperature in the future, it is likely that Alpine glaciers and the high-mountain environment as an entire system will further develop into a state of imbalance. Hence, the assessment of future glacier geometries is a valuable prerequisite for various impact studies. In order to calculate and visualize in a consistent manner future glacier extent for a large number of individual glaciers (> 100) according to a given climate change scenario, we have developed an automated and simple but robust approach that is based on an empirical relationship between glacier size and the steady-state accumulation area ratio (AAR₀) in the Alps. The model requires digital glacier outlines and a digital elevation model (DEM) only and calculates new glacier geometries from a given shift of the steady-state equilibrium line altitude (ELA₀) by means of hypsographic modelling. We have calculated changes in number, area and volume for 3062 individual glacier units in Switzerland and applied six step changes in ELA₀ (from + 100 to + 600 m) combined with four different values of the AAR₀ (0.5, 0.6, 0.67, 0.75). For an AAR₀ of 0.6 and an ELA₀ rise of 200 m (400 m) we calculate a total area loss of − 54% (− 80%) and a corresponding volume loss of − 50% (− 78%) compared to the 1973 glacier extent. In combination with a geocoded satellite image, the future glacier outlines are also used for automated rendering of perspective visualisations. This is a very attractive tool for communicating research results to the general public. Our study is illustrated for a test site in the Upper Engadine (Switzerland), where landscape changes above timberline play an important role for the local economy. The model is seen as a first-step approach, where several parts can be (and should be) further developed.     

Glacier change in Garibaldi Provincial Park, southern Coast Mountains, British Columbia, since the Little Ice Age

Fluctuations of glaciers during the 20th century in Garibaldi Provincial Park, in the southern Coast Mountains of British Columbia, were reconstructed from historical documents, aerial photographs, and fieldwork. Over 505 km², or 26%, of the park, was covered by glacier ice at the beginning of the 18th century. Ice cover decreased to 297 km² by 1987–1988 and to 245 km² (49% of the early 18th century value) by 2005. Glacier recession was greatest between the 1920s and 1950s, with typical frontal retreat rates of 30 m/a. Many glaciers advanced between the 1960s and 1970s, but all glaciers retreated over the last 20 years. Times of glacier recession coincide with warm and relatively dry periods, whereas advances occurred during relatively cold periods. Rapid recession between 1925 and 1946, and since 1977, coincided with the positive phase of the Pacific Decadal Oscillation (PDO), whereas glaciers advanced during its negative phase (1890–1924 and 1947–1976). The record of 20th century glacier fluctuations in Garibaldi Park is similar to that in southern Europe, South America, and New Zealand, suggesting a common, global climatic cause. We conclude that global temperature change in the 20th century explains much of the behaviour of glaciers in Garibaldi Park and elsewhere.     

20 years of mass balances on the Piloto glacier, Las Cuevas river basin, Mendoza, Argentina

Climatic changes of the 20th century have altered the water cycle in the Andean basins of central Argentina. The most visible change is seen in the mountain glaciers, with loss of part of their mass due to decreasing thickness and a substantial recession in the last 100 years. This paper briefly describes the results of glacier mass balance research since 1979 in the Piloto Glacier at the Cajón del Rubio, in the headwaters of Las Cuevas River, presenting new results for the period 1997–2003. Very large interannual variability of net annual specific balance is evident, due largely to variations in winter snow accumulation, with a maximum net annual value of + 151 cm w.e. and a minimum value of - 230 cm w.e. Wet El Niño years are normally associated with positive net annual balances, while dry La Niña years generally result in negative balances. Within the 24-year period, 67% of the years show negative net annual specific balances, with a cumulative mass balance loss of - 10.50 m water equivalent (w.e.). Except for exceptions normally related to El Niño events, a general decreasing trend of winter snow accumulation is evident in the record, particularly after 1992, which has a strong effect in the overall negative mass balance values. The glacier contribution to Las Cuevas River runoff is analysed based on the Punta de Vacas River gauge station for a hypothetical year without snow precipitation (YWSP), when the snowmelt component is zero. Extremely dry years similar to a YWSP have occurred in 1968–1969, 1969–1970 and 1996–1997. The Punta de Vacas gauge station is located 62 km downstream from Piloto Glacier, and the basin contains 3.0% of uncovered glacier ice and 3.7% of debris-covered ice. The total glacier contribution to Las Cuevas River discharge is calculated as 82 ± 8% during extremely dry years. If glacier wastage continues at the present trend as observed during the last 2 decades, it will severely affect the water resources in the arid central Andes of Argentina.     

Glacier ice mass fluctuations and fault instability in tectonically active Southern Alaska

Across the plate boundary zone in south central Alaska, tectonic strain rates are high in a region that includes large glaciers undergoing wastage (glacier retreat and thinning) and surges. For the coastal region between the Bering and Malaspina Glaciers, the average ice mass thickness changes between 1995 and 2000 range from 1 to 5 m/year. These ice changes caused solid Earth displacements in our study region with predicted values of −10 to 50 mm in the vertical and predicted horizontal displacements of 0–10 mm at variable orientations. Relative to stable North America, observed horizontal rates of tectonic deformation range from 10 to 40 mm/year to the north–northwest and the predicted tectonic uplift rates range from approximately 0 mm/year near the Gulf of Alaska coast to 12 mm/year further inland. The ice mass changes between 1995 and 2000 resulted in discernible changes in the Global Positioning System (GPS) measured station positions of one site (ISLE) located adjacent to the Bagley Ice Valley and at one site, DON, located south of the Bering Glacier terminus. In addition to modifying the surface displacements rates, we evaluated the influence ice changes during the Bering glacier surge cycle had on the background seismic rate. We found an increase in the number of earthquakes (M_L≥2.5) and seismic rate associated with ice thinning and a decrease in the number of earthquakes and seismic rate associated with ice thickening. These results support the hypothesis that ice mass changes can modulate the background seismic rate.During the last century, wastage of the coastal glaciers in the Icy Bay and Malaspina region indicates thinning of hundreds of meters and in areas of major retreat, maximum losses of ice thickness approaching 1 km. Between the 1899 Yakataga and Yakutat earthquakes (M_w=8.1, 8.1) and prior to the 1979 St. Elias earthquake (M_s=7.2), the plate interface below Icy Bay was locked and tectonic strain accumulated. We used estimated ice mass change during the 1899–1979 time period to calculate the change in the fault stability margin (FSM) prior to the 1979 St. Elias earthquake. Our results suggest that a cumulative decrease in the fault stability margin at seismogenic depths, due to ice wastage over 80 years, was large, up to 2 MPa. Ice wastage would promote thrust faulting in events such as the 1979 earthquake and subsequent aftershocks.     

Rapid tidewater glacier retreat: a comparison between Columbia Glacier, Alaska and Patagonian calving glaciers

Time-changes in the terminus positions of Patagonian grounded calving glaciers are studied. The framework for the study is a model of terminus retreat based on observations of Columbia Glacier, Alaska. The interpretation favored in this work is that rapid retreat is caused by the terminus thinning to near flotation and weakening of the ice by bottom crevasses. Both thinning and bottom crevasse formation are caused by the increase in near-terminus stretching during rapid retreat. The model is consistent with the behavior of calving Patagonian glaciers regardless of the salinity of the water body calved into. Continuity calculations on Glaciar San Rafael, Chile indicate internal voids, possibly in the form of bottom crevasses, as on Columbia Glacier.     

Tide-induced perturbations of glacier velocities

Recent observations showing substantial diurnal changes in velocities of glaciers flowing into the ocean, measured at locations far inland of glacier grounding lines, add fuel to the ongoing debate concerning the ability of glaciers to transmit longitudinal-stress perturbations over large distances. Resolution of this debate has major implications for the prediction of glacier mass balance, because it determines how rapidly a glacier can respond dynamically to changes such as weakening or removal of an ice shelf. Current IPCC assessment of sea-level rise takes little account of such changes, on the assumption that dynamic responses would be too slow to have any appreciable effect on ice discharge fluxes. However, this assumption must be questioned in view of observations showing massive increases in glacier velocities following removal of parts of the Larsen Ice Shelf, Antarctica, and of others showing diurnal velocity changes apparently linked to the tides.Here, I use a simple force-perturbation model to calculate the response of glacier strain rates to tidal rise and fall, assuming associated longitudinal-force perturbations are transmitted swiftly far inland of the glacier grounding line. Results show reasonable agreement with observations from an Alaskan glacier, where the velocity changes extended only a short distance up-glacier. However, for larger Antarctic glaciers, big velocity changes extending far upstream cannot be explained by this mechanism, unless ice-shelf “back forces” change substantially with the tides.Additional insight will require continuous measurement of velocity and strain-rate profiles along flow lines of glaciers and ice shelves. An example is suggested, involving continuous GPS measurements at a series of locations along the centre line of Glaciar San Rafael, Chile, extending from near the calving front to perhaps 20 km inland. Tidal range here is about ± 0.8 m, which should be sufficient to cause a variation in ice-front velocity of ± 2 cm h^− 1 about its average value of 75 cm h^− 1, assuming local seawater depth of 150 m and glacier thickness of 200–400 m.     

Glacier retreat and climatic variability in the eastern Terskey–Alatoo, inner Tien Shan between the middle of the 19th century and beginning of the 21st century

Changes in the extent of glaciers and rates of glacier termini retreat in the eastern Terskey–Alatoo Range, the Tien Shan Mountains, Central Asia have been evaluated using the remote sensing techniques. Changes in the extent of 335 glaciers between the end of the Little Ice Age (LIA; mid-19th century), 1990 and 2003 have been estimated through the delineation of glacier outlines and the LIA moraine positions on the Landsat TM and ASTER imagery for 1990 and 2003 respectively. By 2003, the glacier surface area had decreased by 19% of the LIA value, which constitutes a 76 km² reduction in glacier surface area. Mapping of 109 glaciers using the 1965 1:25,000 maps revealed that glacier surface area decreased by 12.6% of the 1965 value between 1965 and 2003. Detailed mapping of 10 glaciers using historical maps and aerial photographs from the 1943–1977 period, has enabled glacier extent variations over the 20th century to be identified with a higher temporal resolution. Glacial retreat was slow in the early 20th century but increased considerably between 1943 and 1956 and then again after 1977. The post-1990 period has been marked by the most rapid glacier retreat since the end of the LIA. The observed changes in the extent of glaciers are in line with the observed climatic warming. The regional weather stations have revealed a strong climatic warming during the ablation season since the 1950s at a rate of 0.02–0.03 °C a^− 1. At the higher elevations in the study area represented by the Tien Shan meteorological station, the summer warming was accompanied by negative anomalies in annual precipitation in the 1990s enhancing glacier retreat. However, trends in precipitation in the post-1997 period cannot be evaluated due to the change in observational practices at this station. Neither station in the study area exhibits significant long-term trends in precipitation.     

Cryospheric change in China

This paper provides an overview of the current status of the cryosphere in China and its changes. Up-to-date statistics of the cryosphere in China are summarized based on the latest available data. There are 46,377 glaciers in China, covering an area of 59,425 km². The glacier ice reserve is estimated to be about 5600 km³ and the annual glacier runoff is about 61.6 × 10⁹ m³. The continuous snow cover extent (> 60 days) in China is about 3.4 × 10⁶ km² and the maximum water equivalent is 95.9 × 10⁹ m³ yr^− 1. The permafrost area in China is about 1.72 × 10⁶ km². The total ground ice reserve on the Qinghai–Tibetan Plateau is estimated to be about 10,923 km³. Recent investigations indicated that glacier areas in China have shrunk about 2–10% over the past 45 yr. Total glacier area has receded by about 5.5%. Snow mass has increased slightly. Permafrost is clearly degrading, as indicated by shrinking areas of permafrost, increasing depth of the active layer, rising of lower limit of permafrost, and thinning of the seasonal frost depth. Some models predict that glacier area shrinkage could be as high as 26.7% in 2050, with glacier runoff increasing until its maximum in about 2030. Although snow mass shows an increasing trend in western China, in eastern China the trend is toward decreasing snow mass, with increasing interannual fluctuations. Permafrost degradation is likely to continue, with one-third to one-half of the permafrost on the Qinghai–Tibetan Plateau anticipated to degrade by 2100. Most of the high-temperature permafrost will disappear by then. The permafrost in northeastern China will retreat further northward.     

Twentieth century thinning of Mendenhall Glacier,Alaska, and its relationship to climate,lake calving,and glacier run-off

Mendenhall Glacier is a dynamic maritime glacier in southeast Alaska that is undergoing substantial recession and thinning. The terminus has retreated 3 km during the 20th century and the lower part of the glacier has thinned 200 m or more since 1909. Glacier-wide volume loss between 1948 and 2000 is estimated at 5.5 km³. Wastage has been the strongest in the glacier's lower reaches, but the glacier has also thinned at higher elevations. The shrinkage of Mendenhall Glacier appears to be due primarily to surface melting and secondarily to lake calving. The change in the average rate of thinning on the lower glacier, <1 m a⁻¹ between 1948 and 1982 and >2 m a⁻¹ since 1982, agrees qualitatively with observed warming trends in the region. Mean annual temperatures in Juneau decreased slightly from 1947 to 1976; they then began to increase, leading to an overall warming of ∼1.6 °C since 1943. Lake calving losses have periodically been a small but significant fraction of glacier ablation. The portion of the terminus that ends in the lake is becoming increasingly vulnerable to calving because of a deep pro-glacial lake basin. If current climatic trends persist, the glacier will continue to shrink and the terminus will recede onto land at a position about 500 m inland within one to two decades. The glacier and the meltwaters that flow from it are integral components of the Mendenhall Valley hydrologic system. Approximately 13% of the recent average annual discharge of the Mendenhall River is attributable to glacier shrinkage. Glacier melt contributes 50% of the total river discharge in summer.     

The timing and nature of recession of outlet glaciers of Hielo Patagónico Norte, Chile, from their Neoglacial IV (Little Ice Age) maximum positions

The dates of recession of eleven outlet glaciers of the Hielo Patagónico Norte (Northern Patagonian Icefield) from their recent maximum positions have been inferred from dendrochronology, lichenometry, radiocarbon dating and historical sources. We have refined the dating for part of the Little Ice Age period in this area placing a glacial advance to between AD 1650 and 1766 with the latter date favoured as conformable with historical records and an uncalibrated radiocarbon determination. Glacier recession from maximal positions began in the early 1860s–1870s. Recession was largely synchronous on the western and eastern sides of the Icefield. This synchronicity suggests that climate forcing over-rides second-order controls on glacier behaviour such as the nature of the terminal environment (e.g. calving/non-calving) or differences in glacier drainage basin area. We argue that this icefield-wide glacier recession represents a response to post-Little Ice Age warming, and provides further evidence for the global extent and near synchronous termination of the Little Ice Age.     

Airborne laser altimetry survey of Glaciar Tyndall, Patagonia

The first airborne laser altimetry measurements of a glacier in South America are presented. Data were collected in November of 2001 over Glaciar Tyndall, Torres del Paine National Park, Chilean Patagonia, onboard a Twin Otter airplane of the Chilean Air Force. A laser scanner with a rotating polygon-mirror system together with an Inertial Navigation System (INS) were fixed to the floor of the aircraft, and used in combination with two dual-frequency GPS receivers. Together, the laser–INS–GPS system had a nominal accuracy of 30 cm after data processing. On November 23rd, a total of 235 km were flown over the ablation area of Glaciar Tyndall, with 5 longitudinal tracks with a mean swath width of 300 m, which results in a point spacing of approximately 2 m both along and across track. A digital elevation model (DEM) generated using the laser altimetry data was compared with a DEM produced from a 1975 map (1:50,000 scale — Instituto Geográfico Militar (IGM), Chile). A mean thinning of − 3.1 ± 1.0 m a^− 1 was calculated for the ablation area of Glaciar Tyndall, with a maximum value of − 7.7 ± 1.0 m a^− 1 at the calving front at 50 m a.s.l. and minimum values of between − 1.0 and − 2.0 ± 1.0 m a^− 1 at altitudes close to the equilibrium line altitude (900 m a.s.l.). The thinning rates derived from the airborne survey were similar to the results obtained by means of ground survey carried out at  600 m of altitude on Glaciar Tyndall between 1975 and 2002, yielding a mean thinning of − 3.2 m a^− 1 [Raymond, C., Neumann, T.A., Rignot, E., Echelmeyer, K.A., Rivera, A., Casassa, G., 2005. Retreat of Tyndall Glacier, Patagonia, over the last half century. Journal of Glaciology 173 (51), 239–247.]. A good agreement was also found between ice elevation changes measured with laser data and previous results obtained with Shuttle Radar Topography Mission (SRTM) data. We conclude that airborne laser altimetry is an effective means for accurately detecting glacier elevation changes in Patagonia, where an ice thinning acceleration trend has been observed during recent years, presumably in response to warming and possibly also drier conditions.     

Teleconnections between Andean and New Zealand glaciers

Retreat and advance of glaciers in the Southern Alps of New Zealand have occurred over two distinct 20-yr climate periods (1954–1974) and (1974–1994). Changes in tropical and southern Andean glaciers are compared over these same periods. Behaviour of glaciers in the tropical Andes are out of phase with the Southern Alps glaciers, but some glaciers in Patagonia appear to be in phase. Southern Hemisphere atmospheric circulation using 700 hPa geopotential height anomalies and sea surface temperature patterns are examined for these periods. Glacier response on inter-decadal timescales is linked with distinctive shifts in atmospheric circulation patterns around the Southern Hemisphere. Retreat (advance) of glaciers in the Southern Alps and southern Andean glacier and advance (retreat) of glaciers in the tropical Andes are all associated with weaker (stronger) westerlies, blocking events in the South-east Pacific, negative (positive) geopotential height anomalies over Southern Africa and higher latitudes of the Southern Hemisphere. These glacier changes are also linked with the negative (positive) phase of the Inter-decadal Pacific Oscillation, a higher frequency of La Niña (El Niño) events, and warm (cool) sea surface temperatures in the New Zealand region and cool (warm) sea surface temperatures in the equatorial eastern region of the Pacific Ocean off the coast of Peru.     

Glacier mass balance of tropical Zongo glacier, Bolivia, comparing hydrological and glaciological methods

A glaciological program has been undertaken since 1991 on Zongo glacier in Bolivia (6000–4850 m asl, 2.4 km², 16°S). This program involves mass balance measurements, hydrological studies and energy balance investigations. On outer-tropical glaciers, melting and snow accumulation are both maximum in the wet season (austral summer), whereas the dry season (winter) is a period of low ablation. Errors on each term of the glaciological (stakes, snow-pits and integration method of the measurements) and hydrological (precipitation, discharge and runoff coefficient of free ice areas) methods are investigated to estimate the overall accuracy of the mass balance measurements. The hydrological budget is less than the glaciological one (mean difference: 60 cm w.e. per year), but both methods reproduce similar inter-annual variations. Errors in assessment of evaporation or water storage inside the glacier cannot explain the discrepancy. Errors using the glaciological method are large (around ± 40 cm w.e. per year), but no bias can explain the departure from the hydrological balance. Errors on discharge measurements are small and the uncertainty on the runoff coefficient has a minor effect on the mass balance. We concluded that hydrological budgets are too low due to the catch deficiency of rain gauges and absence of precipitation measurements at high altitudes, emphasizing the difficulty to assess snowfall distribution in high mountainous basins.     

Modelling observed and future runoff from a glacierized tropical catchment (Cordillera Blanca, Perú)

Monthly runoff from the 34.3% glacierized tropical catchment of Llanganuco in the tropical Cordillera Blanca, Perú, is successfully simulated and compared with a measured 44 year time series. In the investigation area, the climate is characterized by all-year round homogenous temperature conditions and a strong variability in air humidity and moisture content of the atmosphere. Thus, contrary to the mid latitudes, the seasonal variation in glacier melt strongly depends on moisture-related variables, rather than on air temperature. The here presented ITGG-2.0-R model aims for these requirements. The lack of moisture-related input data other than precipitation demands for an intermediate calibration step. Net shortwave radiation, the emissivity of the atmosphere and a sublimation/melt ratio are related to precipitation amounts. Runoff is well simulated and correlates with the measured record with r² = 0.76. Seasonally obtained r² are only slightly smaller. On a long-term, the cumulative deviation is minor, and the mean annual cycle of runoff is reproduced rather well (r² = 0.99). Based on four different IPCC climate change scenarios, future runoff is simulated. All runoff scenarios are modelled for the respective steady-state glacier extent. This leads to a reduction in the glacier size and a decreased amount of glacier melt. On the other hand, direct runoff increases due to larger glacier free areas. Consequently, mean annual runoff remains almost unchanged, but the seasonality intensifies considerably with more runoff during the wet and less runoff during the dry season.