Late Glacial ice advances in southeast Tibet

The sensitivity of Tibetan glacial systems to North Atlantic climate forcing is a major issue in palaeoclimatology. In this study, we present surface exposure ages of erratic boulders from a valley system in the Hengduan Mountains, southeastern Tibet, showing evidence of an ice advance during Heinrich event 1. Cosmogenic nuclide analyses (¹⁰Be and ²¹Ne) revealed consistent exposure ages, indicating no major periods of burial or pre-exposure. Erosion-corrected (3 mm/ka) ¹⁰Be exposure ages range from 13.4 to 16.3 ka. This is in agreement with recalculated exposure ages from the same valley system by [Tschudi, S., Schäfer, J.M., Zhao, Z., Wu, X., Ivy-Ochs, S., Kubik, P.W., Schlüchter, C., 2003. Glacial advances in Tibet during the Younger Dryas? Evidence from cosmogenic ¹⁰Be, ²⁶Al, and ²¹Ne. Journal of Asian Earth Sciences 22, 301–306.]. Thus this indicates that local glaciers advanced in the investigated area as a response to Heinrich event 1 cooling and that periglacial surface adjustments during the Younger Dryas overprinted the glacial morphology, leading to deceptively young exposure ages of certain erratic boulders.     

Glacial history of the Spitsbergen archipelago and the problem of a Barents Shelf ice sheet

The popular concept of a Late Weichselian ice sheet covering the Barents Shelf and confluent with the Scandinavian and Russian ice sheets is based primarily on the 6500 B.P. isobase which rises to the east over Spitsbergen, and to the west over Franz Joseph Land. Analysis of uplift curves from the Spitsbergen archipelago shows, however, that the strongest early Holocene uplift occurs over northeastern Spitsbergen and eastern Nordaustlandet, falling both to east and west, and that the centre of uplift migrates to the southeast during the Holocene. Direct evidence of glacier fluctuation indicates an important Billefjorden Stage of glaciation at about 11,000 to 10,000 B.P., part of whose extent can be defined by moraines and by abrupt changes in the marine limit. The dominant ice masses of the Billefjorden Stage seem to have formed over eastern Spitsbergen, Edgeøya, Barentsøya and southern Hinlopenstretet, and it is the decay of this ice mass which is primarily responsible for the pattern of early Holocene uplift. Stratigraphic evidence suggests the absence of an important glacial event at 18,000–20,000 B.P., but an important phase of Spitsbergen-centred glaciation at about 40,000 B.P., and a glacial phase at 80,000–120,000 B.P. It is suggested that many raised beach sequences outside the Billefjorden readvance show an upper sequence related to deglaciation at about 40,000 B.P., and a lower, Holocene sequence related to decay of the Billefjorden ice. The anomalous pattern of late Holocene uplift may be related to restrained rebound produced by regeneration of ice on the main islands of the archipelago and unrestrained rebound on Hopen and Kong Karls Land, which were incapable of sustaining large ice masses of their own. A pattern of LateGlacial climatic circulation which may have produced ice masses on the east coast of Spitsbergen, west coast of Novaya Zemlya and north coast of Russia is suggested. It is also suggested that this pattern of glaciation produced features which have been wrongly interpreted as evidence of a Barents ice sheet.     

Gibbsitic soils on former nunataks: implications for ice sheet reconstruction

X-ray diffraction analyses of soils above and below a periglacial trimline developed across the basalts of the Trotternish Escarpment (Isle of Skye, Scotland) demonstrate that gibbsite is restricted to soils above the trimline. This suggests that the gibbsite is a relict of pre-Late Devensian weathering, and that the trimline did not develop after the last ice sheet achieved its maximum thickness. The sharpness of the boundary between frost-weathered regolith and gibbsitic soils upslope and ice-scoured bedrock associated with gibbsite-free soils downslope suggests that the trimline represents the altitude of the last ice sheet at its maximum thickness rather than a former boundary between passive cold-based ice and erosive warm-based ice. These findings illustrate how identification of high-level periglacial trimlines and associated contrasts in clay mineralogy provide a means for constraining reconstructions of the form of the last ice sheets.     

The Barents ice sheet as a relay regulator of glacial-interglacial alternation

A slight cooling can induce the formation of ice sheets in the Scandinavian mountains and in the American Arctic. The increasing albedo and the appearance of cold air masses above the glaciers cause glaciation to spread over a vast area. As a result, the sea level lowers and a large part of the Barents and Kara seabeds dries up. Ice sheets are formed there, which spread over the northeastern part of the Kola Peninsula, the Pechora River basin, and over northwestern Siberia. The glacier barrier extending nearly from the North Pole to central Europe hinders latitudinal atmospheric circulation. Precipitation decreases sharply in the areas east and southeast of the glaciers. As a consequence, glaciers in the mid-latitudes retreat and sea level rises. Increased iceberg formation is induced in the periphery of the Barents Ice Sheet, causing it to disappear. An interglacial sets in.     

Calibrating a glaciological model of the Greenland ice sheet from the Last Glacial Maximum to present-day using field observations of relative sea level and ice extent

We constrain a three-dimensional thermomechanical model of Greenland ice sheet (GrIS) evolution from the Last Glacial Maximum (LGM, 21 ka BP) to the present-day using, primarily, observations of relative sea level (RSL) as well as field data on past ice extent. Our new model (Huy2) fits a majority of the observations and is characterised by a number of key features: (i) the ice sheet had an excess volume (relative to present) of 4.1 m ice-equivalent sea level at the LGM, which increased to reach a maximum value of 4.6 m at 16.5 ka BP; (ii) retreat from the continental shelf was not continuous around the entire margin, as there was a Younger Dryas readvance in some areas. The final episode of marine retreat was rapid and relatively late (c. 12 ka BP), leaving the ice sheet land based by 10 ka BP; (iii) in response to the Holocene Thermal Maximum (HTM) the ice margin retreated behind its present-day position by up to 80 km in the southwest, 20 km in the south and 80 km in a small area of the northeast. As a result of this retreat the modelled ice sheet reaches a minimum extent between 5 and 4 ka BP, which corresponds to a deficit volume (relative to present) of 0.17 m ice-equivalent sea level. Our results suggest that remaining discrepancies between the model and the observations are likely associated with non-Greenland ice load, differences between modelled and observed present-day ice elevation around the margin, lateral variations in Earth structure and/or the pattern of ice margin retreat.     

Last Glacial Maximum ice sheet dynamics in Arctic Canada inferred from young erratics perched on ancient tors

A long-standing debate regarding the reconstruction of former ice sheets revolves around the use of relative weathering of landscapes, i.e., the assumption that highly weathered landscapes have not been recently glaciated. New cosmogenic isotope measurements from upland bedrock surfaces and erratics along the northeastern margin of the Laurentide Ice Sheet (LIS) shed light on this debate. ¹⁰Be and ²⁶Al concentrations from three perched erratics, yielding cosmogenic exposure ages of 17–11 ka, are much lower than those measured in two unmodified, highly weathered tors upon which they lie, which yield cosmogenic exposure ages of >60 ka. These findings suggest that non-erosive ice covered weathered upland surfaces along the northeastern margin of the LIS during the last glacial maximum. These data challenge the use of relative weathering to define the margins of Pleistocene ice sheets. The juxtaposition of non-erosive ice over upland plateaus and erosive ice in adjacent fiords requires strong gradients in basal thermal regimes, suggestive of an ice-stream mode of glaciation.     

A new ice sheet model validated by remote sensing of the Greenland ice sheet

Accurate prediction of future sea level rise requires models that accurately reproduce and explain the recent observed dramatic ice sheet behaviours. This study presents a new multi-phase, multiple-rheology, scalable and extensible geofluid model of the Greenland ice sheet that shows the credential of successfully reproducing the mass loss rate derived from the Gravity Recovery and Climate Experiment (GRACE), and the microwave remote sensed surface melt area over the past decade. Model simulated early 21st century surface ice flow compares satisfactorily with InSAR measurements. Accurate simulation of the three metrics simultaneously cannot be explained by fortunate model tuning and give us confidence in using this modelling system for projection of the future fate of Greenland Ice Sheet (GrIS). Based on this fully adaptable three dimensional, thermo-mechanically coupled prognostic ice model, we examined the flow sensitivity to granular basal sliding, and further identified that this leads to a positive feedback contributing to enhanced mass loss in a future warming climate. The rheological properties of ice depend sensitively on its temperature, thus we further verified modelâ?s temperature solver against in situ observations. Driven by the NCEP/NCAR reanalysis atmospheric parameters, the ice model simulated GrIS mass loss rate compares favourably with that derived from the GRACE measurements, or about ?147 km³/yr over the 2002–2008 period. Increase of the summer maximum melt area extent (SME) is indicative of expansion of the ablation zone. The modeled SME from year 1979 to 2006 compares well with the cross-polarized gradient ratio method (XPGR) observed melt area in terms of annual variabilities. A high correlation of 0.88 is found between the two time series. In the 30-year model simulation series, the surface melt exhibited large inter-annual and decadal variability, years 2002, 2005 and 2007 being three significant recent melt episodes.     

Periglacial trimlines,former nunataks and the altitude of the last ice sheet in Wester Ross,northwest Scotland

High-level weathering limits separating ice-scoured topography from frost-weathered detritus were identified on 28 mountains in Wester Ross at altitudes of 700–960 m, and a further 22 peaks support evidence of ice scouring to summit level. Weathering limits are defined most clearly on sandstone and gneiss, which have resisted frost shattering during the Late Devensian Lateglacial, but can also be distinguished on schists and quartzite. Schmidt hammer measurements and analyses of clay mineral assemblages indicate significantly more advanced rock and soil weathering above the weathering limits. The persistence of gibbsite above weathering limits indicates that they represent the upper limit of Late Devensian glacial erosion. The regular decline of weathering-limit altitudes along former flowlines eliminates the possibility that the weathering limits represent former thermal boundaries between protective cold-based and erosive warm-based ice. The weathering limits are therefore interpreted as periglacial trimlines that define the maximum surface altitude of the last ice sheet. Calculated basal shear stresses of 50–95 kPa are consistent with this interpretation. Reconstruction of ice-sheet configuration indicates that the former ice-shed lay above 900 m along the present watershed, and that the ice surface descended northwestwards, with broad depressions along major troughs and localised domes around independent centres of ice dispersal. Extrapolation of the ice surface gradient and altitude suggests that the ice sheet did not overrun the Outer Hebrides, but was confluent with the independent Outer Hebrides ice-cap in the North Minch basin. Erratics located up to 140 m above the reconstructed ice surface are inferred to have been emplaced by a pre-Late Devensian ice sheet (or ice sheets) of unknown age. © 1997 John Wiley & Sons, Ltd.     

10Be ages from central east Greenland constrain the extent of the Greenland ice sheet during the Last Glacial Maximum

Traditional ice sheet reconstructions have suggested two distinctly different ice sheet regimes along the East Greenland continental margin during the Last Glacial Maximum (LGM): ice to the shelf break south of Scoresby Sund and ice extending no further than to the inner shelf at and north of Scoresby Sund. We report new ¹⁰Be ages from erratic boulders perched at 250 m a.s.l. on the Kap Brewster peninsula at the mouth of Scoresby Sund. The average ¹⁰Be ages, calculated with an assumed maximum erosion rate of 1 cm/ka and no erosion (respectively, 17.3±2.3 ka and 15.1±1.7 ka) overlap with a period of increased sediment input to the Scoresby Sund fan (19–15 ka). The results presented here suggest that ice reached at least 250 m a.s.l. at the mouth of Scoresby Sund during the LGM and add to a growing body of evidence indicating that LGM ice extended onto the outer shelf in northeast Greenland.     

Late Devensian ice sheet characteristics: a palaeohydraulic approach

Glacial meltwater channels are incised into bedrock and diamicton along much of the length of the Mid-Cheshire Ridge. Detailed mapping of one such system near the town of Helsby reveals a dendritic channel network developed in the opposite direction to the regional ice flow during the last (Late Devensian) glaciation. The channels formed subglacially, under atmospheric and not hydrostatic pressure, presumably as the ice sheet downwasted during deglaciation. Morphological and palaeohydraulic evidence suggests that not all of the network was necessarily active contemporaneously. Former water levels in the channels can be estimated due to the presence of bar surfaces, giving a calculated palaeodischarge of at least 111 m³ s⁻¹. The ablation rates required to account for this large discharge are an order of magnitude greater than those obtained from theoretical calculations and those observed in modern glacial environments. This implies that some form of high-magnitude discharge, such as a seasonal flood event, must have taken place in this area during deglaciation. This picture of the Late Devensian ice sheet suggests that during recession the ice sheet was static, crevassed and relatively thin (<50 m). This study also shows that there is no simple relationship between meltwater channel direction and ice dynamics, and that care is required when using the former to make inferences about the latter. Copyright © 1998 John Wiley & Sons, Ltd.     

Relative sea‐level changes,glacial isostatic modelling and ice‐sheet reconstructions from the British Isles since the Last Glacial Maximum

While contributing <1 m equivalent eustatic sea‐level rise the British Isles ice sheet produced glacio‐isostatic rebound in northern Britain of similar magnitude to eustatic sea‐level change, or global meltwater influx, over the last 18 000 years. The resulting spatially variable relative sea‐level changes combine with observations from far‐field locations to produce a rigorous test for quantitative models of glacial isostatic adjustment, local ice‐sheet history and global meltwater influx. After a review of the attributes of relative sea‐level observations significant for constraining large‐scale models of the isostatic adjustment process we summarise long records of relative sea‐level change from the British Isles and far‐field locations. We give an overview of different global theoretical models of the isostatic adjustment process before presenting intercomparisons of observed and predicted relative sea levels at sites in the British Isles and far‐field for a range of Earth and ice model parameters in order to demonstrate model sensitivity and the resolving power available from using evidence from the British Isles. For the first time we show a good degree of fit between relative sea‐level observations and predictions that are based upon global Earth and ice model parameters, independently derived from analysis of far‐field data, with a terrain‐corrected model of the British Isles ice sheet that includes extensive glaciation of the North Sea and western continental shelf, that does not assume isostatic equilibrium at the Last Glacial Maximum and keeps to trimline constraints of ice surface elevation. We do not attempt to identify a unique solution for the model lithosphere thickness parameter or the local‐scale detail of the ice model in order to provide a fit for all sites, but argue that the next stage should be to incorporate an ice‐sheet model that is based on quantitative, glaciological model simulations. We hope that this paper will stimulate this debate and help to integrate research in glacial geomorphology, glaciology, sea‐level change, Earth rheology and quantitative modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Cirque glaciers as morphological evidence for a thin Younger Dryas ice sheet in east-central southern Norway

Three localities with marginal moraines deposited by former cirque glaciers are investigated in east-central southern Norway. The wet-based (erosive) cirque glaciers with aspects towards S-SW and N-NE are mapped at altitudes above 1100 m, and have a mean equilibrium-line altitude of 1275 m. With a suggested mean annual winter precipitation close to the average for the modern accumulation season (1 October-30 April) when the cirque glaciers existed, the mean air-temperature depression during the ablation season (1 May-30 September) is calculated to be 6–7°C lower than at present. The high-altitude cirques of central Rondane were still covered by ice when the low-altitude cirque glaciers developed in distal position for this massif in eastern Rondane and on isolated mountains. Hence, the cirque glaciers are suggested to have existed during the deglaciation after the Late Weichselian maximum, and most likely during the Younger Dryas (11000–10000 BP). The cirque glaciers indicate a downwasting ice-sheet surface well below an altitude of 1100 m prior to the Younger Dryas, and this supports a limited (small) vertical extent for the Late Weichselian ice sheet in this region. With the contemporaneous level for instantaneous glacierization (glaciation threshold) just below the highest elevated peaks in east-central southern Norway, this fits with the idea of a continuous downwasting of the Late Weichselian ice sheet since the 'first' nunataks appeared. The occurrence of the cirque glaciers indicates a multidomed Scandinavian ice-sheet geometry during the Late Weichselian.     

Reconstructing the Last Glacial Maximum ice sheet in the Weddell Sea embayment,Antarctica, using numerical modelling constrained by field evidence

The Weddell Sea Embayment (WSE) sector of the Antarctic ice sheet has been suggested as a potential source for a period of rapid sea-level rise – Meltwater Pulse 1a, a 20 m rise in ～500 years. Previous modelling attempts have predicted an extensive grounding line advance in the WSE, to the continental shelf break, leading to a large equivalent sea-level contribution for the sector. A range of recent field evidence suggests that the ice sheet elevation change in the WSE at the Last Glacial Maximum (LGM) is less than previously thought. This paper describes and discusses an ice flow modelling derived reconstruction of the LGM ice sheet in the WSE, constrained by the recent field evidence. The ice flow model reconstructions suggest that an ice sheet consistent with the field evidence does not support grounding line advance to the continental shelf break. A range of modelled ice sheet surfaces are instead produced, with different grounding line locations derived from a novel grounding line advance scheme. The ice sheet reconstructions which best fit the field constraints lead to a range of equivalent eustatic sea-level estimates between approximately 1.4 and 3 m for this sector. This paper describes the modelling procedure in detail, considers the assumptions and limitations associated with the modelling approach, and how the uncertainty may impact on the eustatic sea-level equivalent results for the WSE.     

Autochthonous block fields in southern Norway: Implications for the geometry,thickness, and isostatic loading of the Late Weichselian Scandinavian ice sheet

The consistent geographical and altitudinal distribution of autochthonous block fields (mantle of bedrock weathered in situ) and trimlines in southern Norway suggests a multi-domed and asymmetric Late Weichselian ice sheet. Low-gradient ice-sheet profiles in the southern Baltic region, in the North Sea, and along the outer fjord areas of southern Norway, are best explained by movement of ice on a bed of deforming sediment, although water lubricated sliding or a combination of the two, may not be excluded. The ice-thickness distribution of the Late Weichselian Scandinavian ice sheet is not in correspondence with the modern uplift pattern of Fennoscandia. Early Holocene crustal rebound was apparently determined by an exponential, glacio-isostatic rise. Later, however, crustal movements appear to have been dominated by large-scale tectonic uplift of the Fennoscandian Shield, centred on the Gulf of Bothnia, the region of maximum lithosphere thickness.     

Relative sea‐level observations in western Scotland since the Last Glacial Maximum for testing models of glacial isostatic land movements and ice‐sheet reconstructions

Observations of relative sea‐level change and local deglaciation in western Scotland provide critical constraints for modelling glacio‐isostatic rebound in northern Britain over the last 18 000 years. The longest records come from Skye, Arisaig and Knapdale with a shorter, Holocene, record from Kintail. Biostratigraphic (diatom, pollen, dinoflagellate, foraminifera and thecamoebian), lithological and radiocarbon analyses provide age and elevation parameters for each sea‐level index point. All four sites reveal relative sea‐level change that is highly non‐monotonic in time as the local vertical component of glacio‐isostatic rebound and eustasy (or global meltwater influx) dominate at different periods. Copyright © 2006 John Wiley & Sons, Ltd.     

Variations in the pattern and rate of isostatic uplift indicated by a comparison of Holocene sea-level curves from Scotland

Differences in the timing of maxima and minima and shape between sea-level curves from the western Forth valley, lower Strathearn, inner Moray Firth and eastern Solway Firth areas of Scotland are interpreted as being due to differences in isostatic uplift and different methodologies used to interpret the basic data. Possible changes in the relative rates of uplift between areas suggest that the practice of applying the western Forth ‘model’ of relative sea-level change to other parts of Scotland must be questioned.     

Response of a marine ice sheet to changes at the grounding line

A numerical model was designed to study the stability of a marine ice sheet, and used to do some basic experiments. The ice-shelf/ice-sheet interaction enters through the flow law in which the longitudinal stress is also taken into account. Instead of applying the model to some (measured) profile and showing that this is unstable (as is common practice in other studies), an attempt is made to simulate a whole cycle of growth and retreat of a marine ice sheet, although none of the model sheets is particularly sensitive to changes in environmental conditions. The question as to what might happen to the West Antarctic Ice Sheet in the near future when a climatic warming can be expecied as a result of the CO₂ effect, seems to be open for discussion again. From the results presented in this paper one can infer that a collapse, caused by increased melting on the ice shelves, is not very likely.     

Andean tectonics as a consequence of sea-floor spreading

The evolution of sea-floor spreading (directions and rates) in the Southeast Pacific, shows a close relationship with the Andean structural directions and the plutonic episodicity, as well as the succession of short compressive orogenic phases alternating with relatively longer periods of extension.The orogenic periods coincided with the beginning of periods of relatively higher rates of plate convergence, immediately after periods of relative “quiescence”, with the following periodicity: 110-85, 76-70, 63-60, 49-45, 35-33, 16-13, 7, m.y.The presence of oceanic structures as aseismic ridges or fracture zones produces “weak” Benioff zones (shallow and poorly defined), without a corresponding active volcanism. On the other hand, ancient continental structures such as the Precambrian-Paleozoic orogenic belts have active volcanism and neat, deep, and well developed Benioff zones.     

The ice lobes of the Scandinavian ice sheet during the deglaciation in Finland.

On the basis of glacial landforms interpreted by means of Landsat satellite imagery and ice-flow data obtained by other methods, the Scandinavian ice sheet has been observed to have divided at the deglaciation stage into several ice lobes. The ice lobes were more active parts of the uniform ice sheet. They represent parts that had bordered on each other in different directions or on more passive portions of the ice. The reasons for the appearance of separate ice lobes were evidently the Fennoscandian topography, the location of accumulation areas, and regional differences in the amounts of ice generated. In the boundary zones of the different ice lobes, there occur exceptionally large glaciofluvial forms and moraines (interlobate complexes). An area of passive ice was often between ice lobes, and in such areas there occur no noteworthy eskers, marginal formations or streamlined forms. In the part of Finland located on the southern side of the Arctic Circle, six different ice lobes and four major areas of passive ice are interpreted to have existed.