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.     

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.     

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.     

Glacier mass balance variability in the Cordillera Blanca, Peru and its relationship with climate and the large-scale circulation

A 41-year-long reconstructed annual mean glacier mass balance record from the Cordillera Blanca, Peru, was investigated for its climate sensitivity toward temperature, humidity and precipitation, and its links with the large-scale atmospheric circulation. On interannual timescales precipitation variability appears to be the main driver for glacier mass balance fluctuations in the Cordillera Blanca. This is corroborated by an analysis of the relationship between mass balance variations and local- to regional-scale precipitation variability. Temperature tends to enhance precipitation in driving the mass balance signal, as dry years are often characterized by warm conditions, while wet years usually coincide with cold anomalies. In some years, however, warm and wet or cold and dry conditions coincide, under which circumstances temperature minimizes or even neutralizes the effects of precipitation. Surface energy balance studies have shown that changes in atmospheric humidity significantly affect the melt rates of tropical glaciers, but the lack of long and high-quality in-situ measurements precludes a detailed quantitative assessment of its role on interannual timescales in the Cordillera Blanca. Sea surface temperature anomalies (SSTA) in the tropical Pacific exert the dominant large-scale forcing on interannual time scales, leading to negative mass balance anomalies during El Niño and above average mass balance during La Niña episodes. In general the teleconnection mechanism linking ENSO with glacier mass balance is similar to what has previously been described for the Bolivian Altiplano region. Changes in the upper-tropospheric zonal flow aloft associated with ENSO conditions determine the amount of snowfall during the wet season and thereby significantly affect the glacier mass balance. Because this teleconnection mechanism is spatially unstable and oscillates latitudinally along the subtropical Andes, it affects the Cordillera Blanca in most, but not all years. The relationship between ENSO and glacier mass balance is therefore characterized by occasional ‘break downs’, more common since the mid-1970's, when El Niño years with above average mass balance and La Niña events with negative mass balance have been observed.     

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.     

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.     

Snow zonation on Hielo Patagónico Sur, Southern Patagonia, derived from Landsat 5 TM data

Hielo Patagónico Sur (HPS), an icefield in Southern Patagonia, is the largest temperate ice mass in the southern hemisphere. Despite continued research efforts during the last decade many glaciological variables, especially mass balance, are still poorly known. This is partly because access to the icefield is difficult due to remoteness and persistent harsh weather conditions. Therefore, remote sensing appears to be a more suitable tool for the acquisition of data. In this work we present a remote sensing study of snow zonation on HPS using Landsat 5 Thematic Mapper data acquired on 12 March 2001. By using image processing and classification techniques, proved to be useful in other glaciated regions, we map for the first time the extent and occurrence of major snow zones on the whole HPS. We separate between two classes of ice and three classes of snow. Ice facies are classified as bare or debris (i.e. dirt) covered ice, covering 2454 km² and 777 km² respectively, or 18.4% and 5.8% of the icefield on the day of image acquisition. Snow types are classified according to spectral differences in the images, following the glacier facies concept. Two of the three snow cover types are interpreted to represent differences in snow grain size within a fairly homogeneous snow pack whereas the third one is interpreted to represent the slush zone. A first order altitudinal control on the distribution of these snow facies is evident. In addition, our results show that snow accumulation on HPS is markedly controlled by the interaction of strong west–northwest snow-bearing winds and the rough mountainous terrain. In order of decreasing altitude we find that the two snow facies and the slush facies occupy 3819 km², 3292 km² and 2295 km² respectively, or 28.6%, 24.6% and 17.2% of the icefield, on the day of image acquisition. Estimates of equilibrium line altitude using our results yield values of 800–900 m above sea level for the western side and 1500–1600 m above sea level for the eastern side, with an accumulation area ratio of 0.74.     

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.     

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.     

Interannual sea level variability in the South China Sea and its response to ENSO

Sea level observed by altimeter during the 1993–2004 period, thermosteric sea level from 1945 through 2004, and tide gauge records are analyzed to investigate the interannual variability of sea level in the South China Sea (SCS) and its relationship with ENSO (El Niño and Southern Oscillation). Both the interannual variations of the observed sea level and the thermosteric sea level are closely related to ENSO. An ‘enigma’ that the SST and sea level in the SCS have inverse response to ENSO is revealed. It is found that the thermosteric sea level has an excellent correspondence to seawater temperature at 100 m depth, and their variations are unsynchronized to SST. Detailed analysis denotes that the warming of seawater occurs only in the upper 75 m during and after the mature phase of El Niño, while the cooling appears in the layers deeper than 75 m during El Niño years. The volume transports between the SCS and the adjacent oceans and the anomalous Ekman pumping contribute a lot for the sea level fall in the developing stage of El Niño, while the mass exchange, which is dominated by precipitation, plays a more significant role in the following continuous negative sea level anomalies.     

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.     

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.     

Coral Sr/Ca and Mg/Ca records in Con Dao Island off the Mekong Delta: Assessment of their potential for monitoring ENSO and East Asian monsoon

The climate of the South China Sea (SCS) is dominated by the East Asian monsoon (EAM) and can be related to the El Niño-Southern Oscillation (ENSO) owing to the interaction between ENSO and the EAM. An annually-banded coral (Porites sp.) collected from Con Dao Island in the southern SCS was measured for Sr/Ca and Mg/Ca ratios at near-monthly resolution through the annual bands of AD 1948–1999. This island is only ~ 90 km from the Mekong Delta coast and thus significantly influenced by riverine discharge, suggesting relatively severe environmental stress on corals. The Sr/Ca time series shows a clear annual cyclicity chiefly modulated by sea-surface temperature (SST), whereas the Mg/Ca time series exhibits an indistinct annual cyclicity, indicating that the previously-proposed coral Mg/Ca thermometry is greatly disturbed. An instrumental SST record in Con Dao Island (since 1980) has been compared with the Sr/Ca time series to calibrate a Sr/Ca thermometer. The Sr/Ca vs. SST comparison shows that the Sr/Ca thermometer is sometimes disturbed by some factor and that almost all of the disturbances occur around the annual-maximum SST in the warm/wet season. The Sr/Ca data around the annual-minimum SST in the cool/dry season is almost free from the disturbance and thus useful as a SST proxy. The disturbances of the Sr/Ca and Mg/Ca thermometers may be ascribed to the Mekong River discharge and its accompanying phenomena (i.e., large freshwater input, suspended-sediment loads, intense phytoplankton blooms, etc.), which are likely to disturb coral physiological processes. Applying the Sr/Ca thermometer to the whole Sr/Ca time series provides a SST reconstruction from 1948 through 1999. Reconstructed annual-minimum SSTs show a clear quasi-biennial oscillation significantly correlated with ENSO, indicating that the annual-minimum SST in the southern SCS tends to be higher (lower) in El Niño (La Niña) phases. This is compatible with previous observations that the East Asian winter monsoon is weakened (strengthened) in El Niño (La Niña) phases. The reconstructed SST record suggests a warming of 1.0 °C for the latter half of the 20th century. The Sr/Ca and Mg/Ca time series exhibit similar decadal-to-bidecadal variations, which do not seem to be primarily due to SST variability but rather due to some other factor possibly related to disturbance or fluctuation of coral physiological processes. Although both of our Sr/Ca and Mg/Ca data are affected, to a greater or lesser extent, by some non-temperature factor, a part of the Sr/Ca data provides a useful SST proxy and suggests that coral-based SST reconstruction in the southern SCS may be an effective means for monitoring the EAM and ENSO.     

Modeled impacts of changes in tundra snow thickness on ground thermal regime and heat flow to the atmosphere in Northernmost Alaska

Seasonal snow covers the tundra surface for up to nine months of each year on the Alaskan North Slope. Variations in the snow thickness could strongly influence the thermal regime of the underlying soil and permafrost, and the surface energy balance. The impacts of increases and decreases in the tundra snow thickness on the thermal regime of snow surface, active layer, and permafrost, and on the conductive heat flow to the atmosphere were investigated numerically, by using an improved surface energy balance approach based one-dimensional heat transfer model. The baseline inputs for the numerical model are mean daily meteorological data and surface albedos collected at Barrow, Alaska from 1995 through 1999. Based on a study for the long-term mean daily maximum and minimum snow thickness distributions at Barrow in the snow season of 1948 through 1997, a snow thickness factor was defined and five simulation cases were run for the snow season of 1997–1998 by changing the snow thickness factor. The modeled results indicate that changes in snow thickness have significant impacts on ground thermal regimes and conductive heat flow to the atmosphere. Decreasing the snow thickness by 50% led to the maximum ground temperature decrease of 1.48 °C at 0.29 m depth, and 0.72 °C at 3.0 m depth; the magnitude of the mean conductive heat flow to the atmosphere for December increase of 4.3 Wm^− 2. Increasing the snow thickness by 50% resulted in the maximum ground temperature increase of 1.44 °C at 0.29 m depth, and 0.66 °C at 3.0 m depth; the magnitude of the mean conductive heat flow to the atmosphere for December decrease of 1.57 W m^− 2. On an annual basis, variation in the snow thickness by 50%, the ground temperature variations of more than 0.25 °C occurred as deep as 8.0 m below the ground surface. The modeled results also show that changes in snow thickness have a relatively small influence on the snow surface temperature.     

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.     

Estimation of net accumulation rate at a Patagonian glacier by ice core analyses using snow algae

Snow algae in a 45.97-m-long ice core from the Tyndall Glacier (50°59′05″S, 73°31′12″W, 1756 m a.s.l.) in the Southern Patagonian Icefield were examined for potential use in ice core dating and estimation of the net accumulation rate. The core was subjected to visual stratigraphic observation and bulk density measurements in the field, and later to analyses of snow algal biomass, water isotopes (¹⁸O, D), and major dissolved ions. The ice core contained many algal cells that belonged to two species of snow algae growing in the snow near the surface: Chloromonas sp. and an unknown green algal species. Algal biomass and major dissolved ions (Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻, SO₄²⁻) exhibited rapid decreases in the upper 3 m, probably owing to melt water elution and/or decomposition of algal cells. However, seasonal cycles were still found for the snow algal biomass, ¹⁸O, D-excess, and major ions, although the amplitudes of the cycles decreased with depth. Supposing that the layers with almost no snow algae were the winter layers without the melt water essential to algal growth, we estimated that the net accumulation rate at this location was 12.9 m a^− 1 from winter 1998 to winter 1999, and 5.1 m from the beginning of winter to December 1999. These estimates are similar to the values estimated from the peaks of ¹⁸O (17.8 m a^− 1 from summer 1998 to summer 1999 and 11.0 m from summer to December 1999) and those of D-excess (14.7 m a^− 1 from fall 1998 to fall 1999 and 8.6 m a^− 1 from fall to December 1999). These values are much higher than those obtained by past ice core studies in Patagonia, but are of the same order of magnitude as those predicted from various observations at ablation areas of Patagonian glaciers.     

Ice elevation and areal changes of glaciers from the Northern Patagonia Icefield, Chile

High thinning rates (up to − 4.0 ± 0.97 m a^− 1) have been measured at Campo de Hielo Patagónico Norte (CHN) or Northern Patagonia Icefield, Chile between 1975 and 2001. Results have been obtained by comparing a Digital Elevation Model (DEM) derived from regular cartography compiled by Instituto Geográfico Militar of Chile (IGM) based upon 1974/1975 aerial photographs and a DEM generated from Advanced Space-borne Thermal Emission and Reflection Radiometer (ASTER) satellite images acquired in September 2001. A complete cloud-free Landsat ETM+ satellite image mosaic acquired in March 2001 was used to update the available glacier inventory of the CHN, including all glaciers larger than 0.5 km² (48 new glaciers). A new delineation of ice divides was also performed over the accumulation areas of glaciers sharing the high plateau where the existing regular cartography exhibits poor coverage of topographic information. This updated glacier inventory produced a total ice area for 2001 of 3953 km², which represents a decrease of 3.4 ± 1.5% (140 ± 61 km² of ice) with respect to the total ice area of the CHN in 1979 calculated from a Landsat MSS satellite image. Almost 62% of the total area change between 1979 and 2001 took place in glaciers located at the western margin of the CHN, where the maximum area loss was experienced by Glaciar San Quintín with 33 km². At the southern margin, Glaciar Steffen underwent the largest ice-area loss (12 km² or 2.6% of the 1979 area), whilst at the eastern margin the greatest area loss took place in Glaciares Nef (7.9 km², 5.7% of the 1979 area) and Colonia (9.1 km², 2.7% of the 1979 area). At the northern margin of the CHN the lower debris-covered ablation area of Glaciar Grosse collapsed into a new freshwater lake formed during the late 1990s. The areal changes measured at the CHN are much larger than previously estimated due to the inclusion of changes experienced in the accumulation areas. The CHN as a whole is contributing melt water to global sea level rise at rates  25% higher than previous estimates.     

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.     

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.     

Climate change and recent water level variability in Patagonian proglacial lakes, Argentina

A long series of lakes (~ 150) borders the Patagonian Andes (south of ~ 38°S), most of which are a geomorphologic relict of Pleistocene glaciations. Employing instrumental records, we inspected lake water level departures from seasonal variations in seven proglacial lakes: Lacar, Mascardi, Steffen, Escondido, Puelo, Vinter, and Argentino. Lakes north of ~ 42°S show maximum gage (water) level during austral winter months; lakes between ~ 42° and ~ 45°S appear transitional; the one lake south of ~ 50°S (Argentino) shows maximum water level in early autumn. Most lakes show moderate level fluctuation throughout yearly records and, in general, show heteroscedacity. Furthermore, Hurst exponents reveal persistent behavior (i.e., long-term memory effect) in all water level series. In most lakes there are no trends in deseasonalized mean and maximum water levels (Seasonal Kendall test). Lake Mascardi–Manso River system (mostly fed by melt water from the retreating Manso Glacier) is a contrasting example that shows a decreasing trend during summer months that we ascribe to the also declining ice volume. Harmonic analysis (Fourier and wavelet transform) of deseasonalized mean and maximum water level time series shows interannual and decadal periodicities that we link to the occurrence of El Niño and/or the Antarctic Oscillation. The associated phase spectrum indicates that there is a ~ 13-month lag between ENSO occurrences and its effect on anomalous lake water levels. Increased snow accumulation during austral winters usually follows summertime El Niño events, which normally result in increased melt water volume that occurs with about one-year delay during the following (austral) spring/summer.