Distinguishing between tills from Valdaian ice sheets in the Arkhangelsk region, Northwest Russia

Recent studies in the Arkhangelsk region, northwest Russia, have identified at least three consecutive tills all associated with the last Valdaian (Weichselian) glaciation. The Scandinavian ice sheet deposited a Late Valdaian till (ca. 17 ka BP), whereas two tills were deposited in the Early–Middle Valdaian by the Barents/Kara Sea ice sheet (ca. 45–60 ka BP) and an older ice sheet with an eastern centre (ca. 74 ka BP). This article expands on previous stratigraphical work on the discrimination of regional till units by a combination of compositional characteristics and directional properties. Tills associated with the Scandinavian ice sheet were deposited by a glacier advancing from west or northwest, transporting predominantly material from the Fennoscandian shield and the White Sea area. The Barents/Kara Sea ice sheet moved from the north and northeast, whereas the oldest ice advance came from the east–southeast. Although, the two oldest tills both contain material with an eastern provenance, the Viryuga Till is dominated more by local carbonate-rich material. This study demonstrates that detailed investigation of till units facilitate the distinction of glacial events imperative for the reconstructing of the last glaciation in northern Russia.     

A mass balance model for the Eurasian Ice Sheet for the last 120,000 years

We present a mass balance model for Eurasia which is based on the calculation of accumulation from a moisture balance concept. The model is forced with 500 hPa temperatures from GCM time slices at LGM and present day. The model simulates key characteristics, such as control on the size of ice sheets through the advection of moisture, asymmetric ice sheets due to advection of moisture and orography, and the drying of ice sheets when they grow. A simulation of the Eurasian Ice Sheet through a full glacial cycle shows that the model reproduces realistic ice sheets that compare well with geomorphological data. During the Middle Weichselian and the Late Weichselian, the model picks up the trend that the Scandinavian part of the ice grows towards the south and east whilst the ice sheet covering the Barents and Kara Seas remains relatively stable. However, the model seriously underestimates the observed ice extent in the Baltic area. Uncertainties in the temperature and the wind field limit the reliability of regional modelling results.     

Late Weichselian fluvial evolution on the southern Kara Sea Shelf, North Siberia

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1–12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network.These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel–levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land.The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses.Asymmetrical submarine channel–levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel–levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.     

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.     

Evaluating digital elevation models for glaciologic applications: An example from Nevado Coropuna, Peruvian Andes

This paper evaluates the suitability of readily available elevation data derived from recent sensors – the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and the Shuttle Radar Topography Mission (SRTM) – for glaciological applications. The study area is Nevado Coropuna (6426 m), situated in Cordillera Ampato of Southern Peru. The glaciated area was 82.6 km² in 1962, based on aerial photography. We estimate the glacier area to be ca. 60.8 km² in 2000, based on analysis of the ASTER L1B scene.We used two 1:50,000 topographic maps constructed from 1955 aerial photography to create a digital elevation model with 30 m resolution, which we used as a reference dataset. Of the various interpolation techniques examined, the TOPOGRID algorithm was found to be superior to other techniques, and yielded a DEM with a vertical accuracy of ± 14.7 m. The 1955 DEM was compared to the SRTM DEM (2000) and ASTER DEM (2001) on a cell-by-cell basis. Steps included: validating the DEM's against field GPS survey points on rock areas; visualization techniques such as shaded relief and contour maps; quantifying errors (bias) in each DEM; correlating vertical differences between various DEM's with topographic characteristics (elevation, slope and aspect) and subtracting DEM elevations on a cell-by-cell basis.The RMS error of the SRTM DEM with respect to GPS points on non-glaciated areas was 23 m. The ASTER DEM had a RMS error of 61 m with respect to GPS points and displayed 200–300 m horizontal offsets and elevation ‘spikes’ on the glaciated area when compared to the DEM from topographic data.Cell-by-cell comparison of SRTM and ASTER-derived elevations with topographic data showed ablation at the toes of the glaciers (− 25 m to − 75 m surface lowering) and an apparent thickening at the summits. The mean altitude difference on glaciated area (SRTM minus topographic DEM) was − 5 m, pointing towards a lowering of the glacier surface during the period 1955–2000. Spurious values on the glacier surface in the ASTER DEM affected the analysis and thus prevented us from quantifying the glacier changes based on the ASTER data.     

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.     

The chronology of a large ice-dammed lake and the Barents–Kara Ice Sheet advances, Northern Russia

Beach and shoreface sediments deposited in the more than 800-km long ice-dammed Lake Komi in northern European Russia have been investigated and dated. The lake flooded the lowland areas between the Barents–Kara Ice Sheet in the north and the continental drainage divide in the south. Shoreline facies have been dated by 18 optical stimulated luminescence (OSL) dates, most of which are closely grouped in the range 80–100 ka, with a mean of 88±3 ka. This implies that that the Barents–Kara Ice Sheet had its Late Pleistocene maximum extension during the Early Weichselian, probably in the cold interval (Rederstall) between the Brørup and Odderade interstadials of western Europe, correlated with marine isotope stage 5b. This is in strong contrast to the Scandinavian and North American ice sheets, which had their maxima in isotope stage 2, about 20 ka. Field and air photo interpretations suggest that Lake Komi was dammed by the ice advance, which formed the Harbei–Harmon–Sopkay Moraines. These has earlier been correlated with the Markhida moraine across the Pechora River Valley and its western extension. However, OSL dates on fluvial sediments below the Markhida moraine have yielded ages as young as 60 ka. This suggests that the Russian mainland was inundated by two major ice sheet advances from the Barents–Kara seas after the last interglacial: one during the Early Weichselian (about 90 ka) that dammed Lake Komi and one during the Middle Weichselian (about 60 ka). Normal fluvial drainage prevailed during the Late Weichselian, when the ice front was located offshore.     

Cenozoic basin evolution beneath the southern McMurdo Ice Shelf, Antarctica

Fifty-two kilometres of multi-channel seismic reflection data were acquired from the southern McMurdo Ice Shelf (SMIS) during potential drill site investigations for the Antarctic Drilling (ANDRILL) program. The survey was acquired atop 110 to 220 m of floating ice and extended across ablation and accumulation zones of the ice shelf. Seismic processing was tailored to the ice shelf environment, including: datum static corrections to account for changes in the thickness and average velocity of the near-surface firn layer, and changes in the surface elevation across the survey area; residual static corrections to account for near-surface ice shelf irregularities; and two-step predictive deconvolution to suppress ice and firn layer multiples. A model for the ice shelf thickness was also incorporated in the interval velocity model during depth conversion to ensure that the ice shelf structure did not impose non-static shifts on the seismic section.The depth converted CMP stacked sections reveal several N to NE trending normal faults, that offset reflective horizons by up to 150 m within the lower part of the section and form a broad east-dipping, half-graben structure. The seafloor possesses trough and arch morphology in parallel with the half-graben structure. These features are interpreted as the southern extension of the Terror Rift. The rift succession comprises a dislocated (?)early-Miocene synrift package and a relatively undeformed (?)late-Miocene post-rift package separated by an erosional unconformity. The post-rift package infills and onlaps the rift topography, and drapes over the graben system, reaching a maximum thickness of 400 m. Throughout the post-rift phase, the basin was also influenced by Neogene volcanism, evidenced by three small volcanic features within the seismic profiles, and associated successions of inferred volcanic material. An angular unconformity within the post-rift succession is interpreted as a flexural horizon related to the load of Mount Discovery and/or Mount Morning. Up to 150 m of flexural moat fill occurs above this surface at ~ 20 km from the load centres. The post-rift succession also includes several glacio-geomorphic features, the orientation and morphology of which indicate an approximately SW to NE ice flow direction during a mid-Miocene grounding event and a SE to NW ice flow direction during Quaternary ice sheet grounding events.The thickness and lower extent of the rift succession was not able to be determined because signal-to-noise ratio and vertical resolution were low at these depths. Strata from an earlier, Paleogene, rift episode may underlie the Terror Rift succession, or it may be directly underlain by acoustic basement. A Paleogene rift episode has previously been proposed based on the occurrence of Eocene fossiliferous erratics around the margin of the SMIS and the structural setting revealed by the SMIS seismic reflection profiles is consistent with this hypothesis.     

Climate forced atmospheric CO2 variability in the early Holocene: A stomatal frequency reconstruction

The dynamic climate in the Northern Hemisphere during the early Holocene could be expected to have impacted on the global carbon cycle. Ice core studies however, show little variability in atmospheric CO₂. Resolving any possible centennial to decadal CO₂ changes is limited by gas diffusion through the firn layer during bubble enclosure. Here we apply the inverse relationship between stomatal index (measured on sub-fossil leaves) and atmospheric CO₂ to complement ice core records between 11,230 and 10,330 cal. yr BP. High-resolution sampling and radiocarbon dating of lake sediments from the Faroe Islands reconstruct a distinct CO₂ decrease centred on ca. 11,050 cal. yr BP, a consistent and steady decline between ca. 10,900 and 10,600 cal. yr BP and an increased instability after ca. 10,550 cal. yr BP. The earliest decline lasting ca. 150 yr is probably associated with the Preboreal Oscillation, an abrupt climatic cooling affecting much of the Northern Hemisphere a few hundred years after the end of the Younger Dryas. In the absence of known global climatic instability, the decline to ca. 10,600 cal. yr BP is possibly due to expanding vegetation in the Northern Hemisphere. The increasing instability in CO₂ after 10,600 cal. yr BP occurs during a period of increasing cooling of surface waters in the North Atlantic and some increased variability in proxy climate indicators in the region.The reconstructed CO₂ changes also show a distinct similarity to indicators of changing solar activity. This may suggest that at least the Northern Hemisphere was particularly sensitive to changes in solar activity during this time and that atmospheric CO₂ concentrations fluctuated via rapid responses in climate.     

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.     

A lake in Uzboi Vallis and implications for Late Noachian-Early Hesperian climate on Mars

Uzboi Vallis (centered at ∼28°S, 323°E) is ∼400 km long and comprises the southernmost segment of the northward-draining Uzboi-Ladon-Morava (ULM) meso-scale outflow system that emerges from Argyre basin. Bond and Holden craters blocked the valley to the south and north, respectively, forming a Late Noachian-to-Hesperian paleolake basin that exceeded 4000 km³. Limited CRISM data suggest lake deposits in Uzboi and underlying basin floor incorporate relatively more Mg-clays and more Fe-clays, respectively. The short-lived lake overflowed and breached Holden crater’s rim at an elevation of −350 m and rapidly drained into the crater. Fan deltas in Holden extend 25 km from the breach and incorporate meter-sized blocks, and longitudinal grooves along the Uzboi basin floor are hundreds of meters long and average 60 m wide, suggesting high-discharge drainage of the lake. Precipitation-derived runoff rather than regional groundwater or overflow from Argyre dominated contributions to the Uzboi lake, although the failure of most tributaries to respond to a lowering of base level indicates their incision largely ended when the lake drained. The Uzboi lake may have coincided with alluvial and/or lacustrine activity in Holden, Eberswalde, and other craters in southern Margaritifer Terra, where fluvial/lacustrine activity may have required widespread, synoptic precipitation (rain or snow), perhaps associated with an ephemeral, global hydrologic system during the Late Noachian into the Hesperian on Mars.     

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.     

Aram Chaos and its constraints on the surface heat flux of Mars

The surface heat flux of a planet is an important parameter to characterize its internal activity and to determine its thermal evolution. Here we report on a new method to constrain the surface heat flux of Mars during the Hesperian. For this, we explore the consequences for the martian surface heat flux from a recently presented new hypothesis for the formation of Aram Chaos (Zegers, T.E., Oosthoek, J.H.P., Rossi, A.P., Blom, J.K., Schumacher, S. [2010]. Earth Planet. Sci. Lett. 297, 496-504. doi:10.1016/j.epsl.2010.06.049.). In this hypothesis the chaotic terrain is thought to have formed by melting of a buried ice sheet. The slow sedimentation and burial of the ice sheet led to an increased thermal insulation of the ice and subsequently to a temperature increase high enough to trigger melting and the formation of the subsurface lake. As these processes highly depend on the thermal properties of the subsurface and especially on the surface heat flux, it is possible to constrain the latter by using numerical simulations. Based on the hypothesis for the formation of Aram Chaos, we conducted an extensive parameter study to determine the parameter settings leading to sufficient melting of the buried ice sheet. We find that the surface heat flux in the Aram Chaos region during the Hesperian was most likely between 20 and 45 mW m⁻² with a possible maximum value of up to 60 mW m⁻².     

Dynamic topography of the East European craton: Shedding light upon lithospheric structure, composition and mantle dynamics

Most of the East European Craton lacks surface relief; however, the amplitude of topography at the top of the basement exceeds 20 km, the amplitude of topography undulations at the crustal base reaches almost 30 km with an amazing amplitude of ca. 50 km in variation in the thickness of the crystalline crust, and the amplitude of topography variations at the lithosphere–asthenosphere boundary exceeds 200 km. This paper examines the relative contributions of the crust, the subcrustal lithosphere, and the dynamic support of the sublithospheric mantle to maintain surface topography, using regional seismic data on the structure of the crystalline crust and the sedimentary cover, and thermal and large-scale P- and S-wave seismic tomography data on the structure of the lithospheric mantle. For the Precambrian lithosphere, an analysis of Vp/Vs ratio at 100, 150, 200, and 250 km depths does not show any age-dependence, suggesting that while Vp/Vs ratio can be effectively used to outline the cratonic margins, it is not sensitive to compositional variations within the cratonic lithosphere.Statistical analysis of age-dependence of velocity, density, and thermal structure of the continental crust and subcrustal lithosphere in the study area (0–62E, 45–72N) allows to link lithospheric structure with the tectonic evolution of the region since the Archean. Crustal thickness decreases systematically with age from 42–44 km in regions older than 1.6 Ga to 37–40 km in the Paleozoic–Mesoproterozoic structures, and to ca. 31 km in the Meso-Cenozoic regions. However, the isostatic contribution of the crust to the surface topography of the East European Craton is almost independent of age (ca. 4.5 km) due to an interplay of age-dependent crustal and sedimentary thicknesses and lithospheric temperatures.On the contrary, the contribution of the subcrustal lithosphere to the surface topography strongly depends on the age, being slightly positive (+ 0.3 + 0.7 km) for the regions older than 1.6 Ga and negative (− 0.5–1 km) for younger structures. This leads to age-dependent variations in the residual topography, i.e. the topography which cannot be explained by the assumed thermal and density structure of the lithosphere, and which can (at least partly) originate from the dynamic component caused by the mantle flow. Positive dynamic topography at the cratonic margins, which exceeds 2 km in the Norwegian Caledonides and in the Urals, clearly links their on-going uplift with deep mantle processes. Negative residual topography beneath the Archean-Paleoproterozoic cratons (− 1–2 km) indicates either a smaller density deficit (ca. 0.9%) in their subcrustal lithosphere than predicted by global petrologic data on mantle-derived xenoliths or the presence of a strong convective downwelling in the mantle. Such mantle downflows can effectively divert heat from the lithospheric base, leading to a long-term survival of the Archean-Paleoproterozoic lithosphere.     

Topographic control on the dynamics of the Svalbard-Barents Sea ice sheet

The land surface of what is now the Barents Sea region may have been eroded to a sub-aerial platform prior to the Quaternary, due to both tectonic uplift-induced and sea-level lowering-induced erosion processes. The Barents Sea was then further eroded into its present form by the subsequent action of ice sheets. Two bedrock configurations, representing the pre-Quaternary sub-aerial Barents Shelf topography and the largely submarine morphology of the present day, were used as input to a glaciological ice sheet model so that the dynamic evolution of the maximum-sized ice sheets, caused solely by a change in bedrock elevation, could be identified. The ice-sheet model was run under constant glacial environmental conditions, until mass balance stability was reached, over both bedrock configurations. The simple parabolic ice sheet surface, which formed on a flat sub-aerial bedrock platform, was found to be significantly different in dynamic character compared with an ice sheet developed on the present submarine bedrock topography. In this latter situation, the central ice dome is drained by ice streams in Bjørnøyrenna, Storfjordrenna and smaller outlet glaciers in the north of the ice sheet.     

Comparison of Mediterranean sea level fields for the period 1961–2000 as given by a data reconstruction and a 3D model

Two Mediterranean sea level distributions spanning the last decades are examined. The first one is a reconstruction of sea level obtained by a reduced-space optimal interpolation applied to tide gauge and altimetry data. The second distribution is obtained from a 3D (baroclinic) regional circulation model. None of the two representations includes the mechanical atmospheric forcing. Results are presented for two different periods: 1993–2000 (for which altimetry data are available) and 1961–2000 (the longest period common to both distributions).The first period is examined as a test period for the model, since the reconstruction is very similar to altimetry observations. The modelled sea level is in fair agreement with the reconstruction in the Western Mediterranean and in the Aegean Sea (except in the early nineties), but in the Ionian Sea the model departs from observations. For the whole period 1961–2000 the main feature is a marked positive trend in the Ionian Sea (up to 1.8 mm yr^− 1), observed both in the reconstruction and in the model. Also the distribution of positive trends in the Western Mediterranean (mean value of 1.1 mm yr^− 1) and the smaller trends in the Aegean Sea (0.5 mm yr^− 1) are similar in the reconstruction and in the model, despite the first implicitly accounts for sea level variations due to remote sources such as ice melting and the second does not. The interannual sea level variability associated with key regional events such as the Eastern Mediterranean Transient is apparently captured by the reconstruction but not by the model (at least in its present configuration). Hence, the reconstruction can be envisaged as a useful tool to validate further long-term numerical simulations in the region.     

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.     

Permafrost warming in the Tien Shan Mountains, Central Asia

The general features of alpine permafrost such as spatial distribution, temperatures, ice content, permafrost and active-layer thickness within the Tien Shan Mountains, Central Asia are described. The modern thermal state of permafrost reflects climatic processes during the twentieth century when the average rise in mean annual air temperature was 0.006–0.032 °C/yr for the different parts of the Tien Shan. Geothermal observations during the last 30 yr indicate an increase in permafrost temperatures from 0.3 °C up to 0.6 °C. At the same time, the average active-layer thickness increased by 23% in comparison to the early 1970s. The long-term records of air temperature and snow cover from the Tien Shan's high-mountain weather stations allow reconstruction of the thermal state of permafrost dynamics during the last century. The modeling estimation shows that the altitudinal lower boundary of permafrost distribution has shifted by about 150–200 m upward during the twentieth century. During the same period, the area of permafrost distribution within two river basins in the Northern Tien Shan decreased approximately by 18%. Both geothermal observations and modeling indicate more favorable conditions for permafrost occurrences and preservation in the coarse blocky material, where the ice-rich permafrost could still be stable even when the mean annual air temperatures exceeds 0 °C.     

Are sea-level-rise trends along the coasts of the north Indian Ocean consistent with global estimates?

Mean-sea-level data from coastal tide gauges in the north Indian Ocean were used to show that low-frequency variability is consistent among the stations in the basin. Statistically significant trends obtained from records longer than 40 years yielded sea-level-rise estimates between 1.06–1.75 mm yr^− 1, with a regional average of 1.29 mm yr^− 1, when corrected for global isostatic adjustment (GIA) using model data. These estimates are consistent with the 1–2 mm yr^− 1 global sea-level-rise estimates reported by the Intergovernmental Panel on Climate Change.     

The North Taymyr ice-marginal zone, Arctic Siberia—a preliminary overview and dating

The North Taymyr ice-marginal zone (NTZ) is a complex of glacial, glaciofluvial and glaciolacustrine deposits, laid down on the northwestern Taymyr Peninsula in northernmost Siberia, along the front of ice sheets primarily originating on the Kara Sea shelf. It was originally recognised from satellite radar images by Russian scientists; however, before the present study, it had not been investigated in any detail. The ice sheets have mainly inundated Taymyr from the northwest, and the NTZ can be followed for 700–750 km between 75°N and 77°N, mostly 80–100 km inland from the present Kara Sea coast.The ice-marginal zone is best developed in its central parts, ca. 100 km on each side of the Lower Taymyr River, and has there been studied by us in four areas. In two of these, the ice sheet ended on land, whereas in the two others, it mainly terminated into ice-dammed lakes. The base of the NTZ is a series of up to 100-m-high and 2-km-wide ridges, usually consisting of redeposited marine silts. These ridges are still to a large extent ice-cored; however, the present active layer rarely penetrates to the ice surface. Upon these main ridges, smaller ridges of till and glaciofluvial material are superimposed. Related to these are deltas corresponding to two generations of ice-dammed lakes, with shore levels at 120–140 m and ca. 80 m a.s.l. These glacial lakes drained southwards, opposite to the present-day pattern, via the Taymyr River valley into the Taymyr Lake basin and, from there, most probably westwards to the southern Kara Sea shelf.The basal parts of the NTZ have not been dated; however, OSL dates of glaciolacustrine deltas indicate an Early–Middle Weichselian age for at least the superimposed ridges. The youngest parts of the NTZ are derived from a thin ice sheet (less than 300 m thick near the present coast) inundating the lowlands adjacent to the lower reaches of the Taymyr River. The glacial ice from this youngest advance is buried under only ca. 0.5 m of melt-out till and is exposed by hundreds of shallow slides. This final glaciation is predated by glacially redeposited marine shells aged ca. 20,000 BP (¹⁴C) and postdated by terrestrial plant material from ca. 11,775 and 9500 BP (¹⁴C)–giving it a last global glacial maximum (LGM; Late Weichselian) age.