Quaternary uplift of northern England

Upland flats, attributable to erosion, have long been recognised in the landscape of the Lake District region of NW England, at altitudes of up to ~ 800 m O.D. Extrapolation using uplift rates derived from dated Pleistocene sites (karstic caves and other features) in the adjacent Pennine uplands suggests that if this succession of flats formed close to sea-level they date from the Middle Pliocene onwards, indicating a subsequent time-averaged uplift rate of almost 0.3 mm a^− 1. Numerical modelling indicates that erosion of surrounding areas at a typical rate of 0.2 mm a^− 1 since 3.1 Ma could have caused this uplift, as well as constraining the local effective viscosity of the lower crust as ~ 4 × 10¹⁸ Pa s and the typical local Moho temperature as ~ 650 °C. It is thus feasible that most of the topography of northern England has developed since the Middle Pliocene, as a consequence of coupling between erosion and the resulting induced flow in the lower continental crust. The much faster vertical crustal motions indicated in this part of northern England, compared with SE England, are thus mainly a consequence of much greater mobility of the lower crust in the north, due to its younger thermal age and the heating effect of radioactive Palaeozoic granites. Uplift of this magnitude, which has previously gone unrecognised, may have affected post-Pliocene global climate.     

Late Cenozoic uplift of western Turkey: Improved dating of the Kula Quaternary volcanic field and numerical modelling of the Gediz River terrace staircase

A set of 13 new unspiked K–Ar dates has been obtained for the Quaternary basaltic volcanism in the Kula area of western Turkey, providing improved age control for the fluvial deposits of the Gediz River that underlie these basalt flows. This dating is able, for the first time, to resolve different ages for the oldest basalts, assigned to category β2, that cap the earliest Gediz deposits recognised in this area, at altitudes of 140 to 210 m above present river level. In particular, the β2 basalt capping the Sarnıç Plateau is dated to 1215 ± 16 ka (± 2σ), suggesting that the youngest underlying fluvial deposits, 185 m above present river level, are no younger than marine oxygen isotope stage (MIS) 38. In contrast, the β2 basalt capping the adjacent Burgaz Plateau is dated to 1014 ± 23 ka, suggesting that the youngest underlying fluvial deposits, 140 m above present river level, date from MIS 28. The staircase of 11 high Gediz terraces capping the latter plateau is thus dated to MIS 48-28, assuming they represent consecutive 40 ka Milankovitch cycles, although it is possible that as many as two cycles are missing from this sequence such that the highest terrace is correspondingly older. Basalt flows assigned to the β3 category, capping Gediz terraces 35 and 25 m above the present river level, have been dated to 236 ± 6 ka and 180 ± 5 ka, indicating incision rates of 0.15 mm a^− 1, similar to the time-averaged rates since the eruptions of the β2 basalts. The youngest basalts, assigned to category β4, are Late Holocene; our K–Ar results for them range from zero age to a maximum of 7 ± 2 ka.This fluvial incision is interpreted using numerical modelling as a consequence of uplift caused by a regional-scale increase in spatial average erosion rates to 0.1 mm a^− 1, starting at 3100 ka, caused by climate deterioration, since when a total of 410 m of uplift has occurred. Parameters deduced on this basis from the observed disposition of the Early Pleistocene Gediz terraces include the local effective viscosity of the lower crust, which is 2 × 10¹⁸ Pa s, the Moho temperature of 660 °C, and the depth of the base of the brittle upper crust, which is 13 km. The thin lithosphere in this area results in high heat flow, causing this relatively shallow base of the brittle upper crust and the associated relatively thick lower-crustal layer, situated between depths of 13 and 30 km. It estimated that around 900 ka, at the start of the 100 ka Milankovitch forcing, the spatial average erosion rate increased slightly, to 0.12 mm a^− 1; the associated relatively sluggish variations in uplift rates are as expected given the relatively thick lower-crustal layer.This modelling indicates that the growth of topography since the Pliocene in this study region has not involved a steady state. The landscape was significantly perturbed by the Middle Pliocene increase in erosion rates, and has subsequently adjusted towards—but not reached—a new steady state consistent with these increased erosion rates. It would not be possible to constrain what has been occurring from the Middle to Late Pleistocene or even the Early Pleistocene uplift response alone; information regarding the starting conditions is also essential, this being available in this region from the older geological record of stacked fluvial and lacustrine deposition. This result has major implications for the rigorous modelling of uplift histories in regions of rapid erosion, where preservation of information to constrain the starting conditions is unlikely.     

Climatic implications of correlated Upper Pleistocene glacial and fluvial deposits on the Cinca and Gállego Rivers (NE Spain) based on OSL dating and soil stratigraphy

We correlate Upper Pleistocene glacial and fluvial deposits of the Cinca and Gállego River valleys (south central Pyrenees and Ebro basin, Spain) using geomorphic position, luminescence dates, and time-related trends in soil development. The ages obtained from glacial deposits indicate glacial periods at 85 ± 5 ka, 64 ± 11 ka, and 36 ± 3 ka (from glacial till) and 20 ± 3 ka (from loess). The fluvial drainage system, fed by glaciers in the headwaters, developed extensive terrace systems in the Cinca River valley at 178 ± 21 ka, 97 ± 16 ka, 61 ± 4 ka, 47 ± 4 ka, and 11 ± 1 ka, and in the Gállego River valley at 151 ± 11 ka, 68 ± 7 ka, and 45 ± 3 ka. The times of maximum geomorphic activity related to cold phases coincide with Late Pleistocene marine isotope stages and Heinrich events. The maximum extent of glaciers during the last glacial occurred at 64 ± 11 ka, and the terraces correlated with this glacial phase are the most extensive in both the Cinca (61 ± 4 ka) and Gállego (68 ± 7 ka) valleys, indicating a strong increase in fluvial discharge and availability of sediments related to the transition to deglaciation. The global Last Glacial Maximum is scarcely represented in the south central Pyrenees owing to dominantly dry conditions at that time. Precipitation must be controlled by the position of the Iberian Peninsula with respect to the North Atlantic atmospheric circulation system. The glacial systems and the associated fluvial dynamic seem sensitive to 1) global climate changes controlled by insolation, 2) North Atlantic thermohaline circulation influenced by freshwater pulses into the North Atlantic, and 3) anomalies in atmospheric circulation in the North Atlantic controlling precipitation on the Iberian Peninsula. Our scenario of glacial and fluvial evolution during the Late Pleistocene in northern Spain could be extrapolated to other glaciated mountainous areas in southern Europe.     

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.     

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.     

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.     

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.     

Reconstruction of Mediterranean sea level fields for the period 1945–2000

The distribution of sea level in the Mediterranean Sea is recovered for the period 1945–2000 by using a reduced space optimal interpolation analysis. The method involves estimating empirical orthogonal functions from satellite altimeter data spanning the period 1993–2005 that are then combined with tide gauge data to recover sea level fields over the period 1945–2000. The reconstruction technique is discussed and its robustness is checked through different tests. For the altimetric period (1993–2000) the prediction skill is quantified over the whole domain by comparing the reconstructed fields with satellite altimeter observations. For past times the skill can only be tested locally, by validating the reconstruction against independent tide gauge records. The reconstructed distribution of sea level trends for the period 1945–2000 shows a positive peak in the Ionian Sea (up to 1.5 mm yr^− 1) and a negative peak of − 0.5 mm yr^− 1 in a small area to the south-east of Crete. Positive trends are found nearly everywhere, being larger in the western Mediterranean (between 0.5 and 1 mm yr^− 1) than in the eastern Mediterranean (between 0 and 0.5 mm yr^− 1). The estimated rate of mean sea level rise for the period 1945–2000 is 0.7 ± 0.2 mm yr^− 1, i.e. about a half of the rate estimated for global mean sea level. These overall results do not appear to be very sensitive to the distribution of tide gauges. The poorest results are obtained in open-sea regions with intense mesoscale variability not correlated with any tide gauge station, such as the Algerian Basin.     

Growth patterns of the last ice age coral terraces at Huon Peninsula

At Huon Peninsula, Papua New Guinea, prolific coral growth during the last-glacial was episodic and in response to a series of sea-level rises. The resultant step-like coral terraces are currently situated from 20 m up to 140 m above sea-level due to continuous tectonic uplift of the Peninsula. The sea-level rises were in response to periodic partial disintegration of Northern Hemisphere ice sheets associated with severe climate swings and occurred within decadal timescales. The relatively rapid 15 m to 35 m rise in sea-levels exposed new head-room for corals to colonize. The resulting terrace structures contain individual corals that do not appear to have grown sequentially in time and with elevation. Additionally, following the peak, sea level fell relatively slowly over several thousand years and corals grew and filled in the flanks of the terrace such that younger corals now occupy lower elevations. We have labeled these structures “pack-up” reefs. This is in contrast to coral terraces formed during major sea-level rises from glacial to interglacial or glacial to interstadial transitions where the rate of sea level rise is commensurate with coral growth rates and corals can keep up with sea-level rise by growing on top of each other in a time orderly sequence. Deriving sea-level information from pack-up terraces is difficult and is likely to be ambiguous. The periodic fluctuations in climate were associated with atmospheric radiocarbon swings that seem to have varied smoothly with time. The same corals that show a scatter in stratigraphic temporal ordering appear regularly distributed in time and with radiocarbon content attesting to the veracity of the age measurements and at the same time confirm the disordered distribution of corals in “pack-up” type reefs.     

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.     

Preservation by basalt of a staircase of latest Pliocene terraces of the River Murat in eastern Turkey: Evidence for rapid uplift of the eastern Anatolian Plateau

Unspiked K–Ar dating makes the age of the Çakmaközü basalt in eastern Turkey 1818 ± 39 ka (± 2σ). This basalt overlies a staircase of four terraces of the River Murat, a Euphrates tributary, each separated vertically by  20 m. We deduce from the relationship with the basalt that these fluvial deposits aggraded during successive  40 ka climate cycles around the Pliocene–Pleistocene boundary (probably MIS 72-66). The incision and rock uplift at  0.5 mm a^− 1, thus indicated, are roughly consistent with the  500 m of entrenchment of this  1.8 Ma Murat palaeo-valley into a former lake basin since the Mid-Pliocene climatic optimum. We infer that the  130 m of incision in this locality since  1.8 Ma dramatically underestimates the associated rock uplift, estimated as  600 m. The  1100 m of rock uplift and  800 m of surface uplift thus estimated since the Mid-Pliocene indicate (assuming Airy isostatic equilibrium)  5 km of thickening of the continental crust, from  37 km to the present 42 km. Eastern Anatolia was thus at a much lower altitude in the Mid-Pliocene than at present, consistent with the low-relief lacustrine palaeo-environment. We infer that the subsequent development of topography and excess crustal thickness are being caused by coupling between surface processes and induced flow in the lower crust: climate change following the Mid-Pliocene climatic optimum resulted in faster erosion that has drawn mobile lower crust beneath the study region.     

Recent organic matter accumulation in relation to some climatic factors in ombrotrophic peat bogs near heavy metal emission sources in Finland

Accumulation of organic matter (OM) was studied in four ombrotrophic peat bogs in Finland: Harjavalta (vicinity of a Cu–Ni smelter), Outokumpu (near a closed Cu–Ni mine), Alkkia (Ni-treated site) and Hietajärvi (a pristine site). At each sampling site, two peat cores (15 × 15 × 100 cm) were taken. Age-dating of peat was determined using ²¹⁰Pb method (CRS model). The local annual temperature sum and precipitation for the past 125 years were modeled. The objective was to compare recent net accumulation rates of heavy metal polluted ombrotrophic peat bogs with those of a pristine bog, and to study the relationship between weather and net accumulation rates. Based on ²¹⁰Pb age-dating, the upper 16-cm peat layer at Harjavalta, 35 cm at Outokumpu and 25 cm at Hietajärvi represents 125 years of peat formation, yielding the following average peat accumulation rates: Harjavalta 1.3 mm year^− 1, Outokumpu 2.8 mm year^− 1 and Hietajärvi 2.0 mm year^− 1. At the Alkkia site, the Ni treatment in 1962 had completely stopped the peat accumulation. Net accumulation rates were related to precipitation at Outokumpu, Harjavalta and Hietajärvi sites. In addition, emissions released from the nearby located Cu–Ni smelter could have affected negatively net OM accumulation rate at Harjavalta site.     

A peat core based estimate of Late-glacial and Holocene methane emissions from northern peatlands

We estimate the intensity of Late-glacial and Holocene methane emissions from peatlands based on their paleo net primary production (PNPP). The PNPP is derived from the carbon accumulation rates of the studied bog profile (Etang de la Gruère, Switzerland), which are corrected for the degree of peat degradation. The obtained PNPP curve is taken as a proxy for methane emissions. It shows relatively high values (90 g C m^− 2 yr^− 1) early in the Bolling/Allerod and drops to low values (40 g C m^− 2 yr^− 1) during the Younger Dryas cold period. With the onset of the Holocene the PNPP increases strongly up to 150 g C m^− 2 yr^− 1 around ca. 10,000 Cal. yr bp. This is followed by a decline to minimum values (30 to 40 g C m^− 2 yr^− 1) between 6500 and 4000 Cal. yr bp. Thereafter, the PNPP starts to increase again to reach its highest value (175 g C m^− 2 yr^− 1) around 1000 Cal. yr bp.The PNPP curve correlates well with the evolution of the atmospheric methane concentrations as derived from Greenland ice-cores. For example, minima in atmospheric methane reported during the Younger Dryas and around 5200 Cal. yr bp are coinciding with the lowest values of PNPP and the negative atmospheric methane peak at 8200 Cal. yr bp corresponds to a marked decrease in PNPP.Our PNPP curve suggests that the methane emissions from northern peatlands evolved similar to those of low latitude wetlands and together they largely determined the evolution of atmospheric methane throughout the Late-glacial and the Holocene. The abruptness of the rise of atmospheric methane at the end of the Younger Dryas probably points to an additional source (e.g. marine gas hydrates), but very early in the Holocene the peatlands have likely become the dominant source of atmospheric methane.     

Source, sea level and circulation effects on the sediment flux to the deep ocean over the past 15 ka off eastern New Zealand

The last post-glacial transgression and present highstand of sea level were accompanied by a reduction in the terrigenous flux to the deep ocean bordering the active convergent margin off the eastern North Island of New Zealand. Although in accord with long-established models of highstand shelf deposition, new data from giant piston core MD97 2121 (2314 m depth) reveal that the flux also varied with terrigenous supply and palaeocirculation. Between 15 and 9.5 ka, the flux reduced from 33 to 20 g/cm²/ka as supply declined with an expanding vegetation cover, and mud depocentres became established on the continental shelf. An increase from 20 to 27 g/cm²/ka during 9.5–3.5 ka coincided with a strengthened East Cape Current which probably introduced sediment from fluvial and shelf sources in the north. The flux profile shows no immediate response to the establishment of modern sea level 7 ka. However, accumulation decreased from 3.5 to 1 ka as more sediments were retained on the shelf, possibly under wind-strengthened, along-shelf currents. Over the last 1 ka, the flux decline halted under increased terrigenous supply during anthropogenic development of the land.Despite the proximity of the North Island's Central Volcanic Region, major eruptions caused only brief increases (centuries duration) in the terrigenous flux through direct deposition of airfall and possibly fluvial redistribution of onshore volcanic deposits. Frequent earthquakes also had little short-term effect on accumulation although such events, along with volcanism, probably contribute to the long-term high flux of the region.The other measured flux component, biogenic carbonate, reached maxima of 6 g/cm²/ka between 11 and 8.5 ka when nutrient-bearing waters of the East Cape Current dominated the palaeoceanography. After these peaks, carbonate accumulation declined gradually to modern levels of 3 g/cm²/ka.     

Implications of new Quaternary uplift models for correlation between the Middle and Upper Thames terrace sequences, UK

Modelling of uplift histories in the Upper and Middle Thames valleys has revealed an important difference, in the form of additional early post-Anglian uplift in the Middle Thames, attributed to an Anglian (~ 440 ka) glacio-isostatic effect. Terraces in the Upper Thames around Oxford seem unaffected by glacio-isostasy and their heights show regional uplift of ~ 35–40 m since the Anglian. The result of the glacio-isostasy is that Anglian terraces are significantly higher above the valley floor in the Middle Thames (up to 55 m) than in the Upper Thames. Recognition of this displacement of Middle Thames terraces has solved long-standing problems of correlation between this area and the Upper Thames: the pre-Anglian (Cromerian Complex) age of the Sugworth Channel deposits, indicated by biostratigraphy, is no longer a difficulty, whereas the Hanborough Terrace is now thought to be of Anglian age, albeit incorporating pre-Anglian faunal remains and perhaps with a significant early post-Anglian component. These findings have implications for the understanding of the effects of Middle Pleistocene glacio-isostasy and of landscape evolution on the periphery of glaciated regions.     

Sea-level curves: Perplexities of an end-user in morphotectonic applications

Several sea-level curves, from the classical ones to some of the most recently published, are analysed and compared. Although this short paper is not a review of the original data and analytical methods used by the different authors, similarities and mainly differences between published sea-level curves are emphasised especially as concern the suggested height and age of the highstand sea-level peaks during the last 450 ka. The perplexities of an end-user when attempting to exploit different curves for investigating uplifted marine terraces are discussed and recommendations are suggested.     

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.     

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.     

Late Pleistocene and Holocene environments in the Nile basin

Owing to the very gently sloping nature of the flood plain in the lower White Nile valley, which is underlain by a former lake-bed, the depositional record in that area is unusually well preserved. In Egypt and along the Blue Nile phases of erosion have destroyed segments of the sedimentary record, but the White Nile sequence is a good proxy for both the main Nile and the Blue Nile. During the last 15 ka, at least, times of high flow in the Blue Nile and main Nile were synchronous with those in the White Nile.Not all the White Nile flood deposits have been preserved but calibrated radiocarbon dates obtained on fossil freshwater and amphibious Pila shells and fish bones indicate that White Nile levels were high around 14.7–13.1 ka, 9.7–9.0 ka, 7.9–7.6 ka, 6.3 ka and 3.2–2.8 ka. The Blue Nile record is more fragmentary and that of the main Nile even more so except for the Holocene Nile delta. Calibrated radiocarbon ages for high Blue Nile flows indicate very high flood levels towards 13.9–13.2 ka, 8.6 ka, 7.7 ka and 6.3 ka.Incision by the Blue Nile and main Nile has caused progressive incision in the White Nile amounting to at least 4 m since the terminal Pleistocene  15 ka ago and at least 2 m over the past 9 ka. The Blue Nile seems to have cut down at least 10 m since  15 ka and at least 4 m since 9 ka. The time-transgressive and relatively late inception of plant domestication in the Nile valley may partly reflect this history of incision. Nile incision would propagate upstream into the White Nile valley, draining previously swampy areas along the valley floor, which would then become accessible to cultivation.     

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.