Ice-stream initiation, duration and thinning on James Ross Island, northern Antarctic Peninsula

Predicting the future response of the Antarctic Ice Sheet to climate change requires an understanding of the ice streams that dominate its dynamics. Here we use cosmogenic isotope exposure-age dating (26Al, 10Be and 36Cl) of erratic boulders on ice-free land on James Ross Island, north-eastern Antarctic Peninsula, to define the evolution of Last Glacial Maximum (LGM) ice in the adjacent Prince Gustav Channel. These data include ice-sheet extent, thickness and dynamical behaviour. Prior to ∼18 ka, the LGM Antarctic Peninsula Ice Sheet extended to the continental shelf-edge and transported erratic boulders onto high-elevation mesas on James Ross Island. After ∼18 ka there was a period of rapid ice-sheet surface-lowering, coincident with the initiation of the Prince Gustav Ice Stream. This timing coincided with rapid increases in atmospheric temperature and eustatic sea-level rise around the Antarctic Peninsula. Collectively, these data provide evidence for a transition from a thick, cold-based LGM Antarctic Peninsula Ice Sheet to a thinner, partially warm-based ice sheet during deglaciation.

Glaciodynamics of the central sector of the last British–Irish ice sheet in Northern England

The central sector of the last British–Irish Ice Sheet (BIIS) was characterised by considerable complexity, both in terms of its glacial stratigraphy and geomorphological signature. This complexity is reflected by the large number and long history of papers that have attempted to decipher the glaciodynamic history of the region. Despite significant advances in our understanding, reconstructions remain hotly debated and relatively local, thereby hindering attempts to piece together BIIS dynamics. This paper seeks to address these issues by reviewing geomorphological mapping evidence of palimpsest flow signatures and providing an up-to-date stratigraphy of the region. Reconciling geomorphological and sedimentological evidence with relative and absolute dating constraints has allowed us to develop a new six-stage glacial model of ice-flow history and behaviour in the central sector of the last BIIS, with three major phases of glacial advance. This includes: I. Eastwards ice flow through prominent topographic corridors of the north Pennines; II. Cessation of the Stainmore ice flow pathway and northwards migration of the North Irish Sea Basin ice divide; III. Stagnation and retreat of the Tyne Gap Ice Stream; IV. Blackhall Wood–Gosforth Oscillation; V. Deglaciation of the Solway Lowlands; and VI. Scottish Re-advance and subsequent final retreat of ice out of the central sector of the last BIIS. The ice sheet was characterised by considerable dynamism, with flow switches, initiation (and termination) of ice streams, draw-down of ice into marine ice streams, repeated ice-marginal fluctuations and the production of large volumes of meltwater, locally impounded to form ice-dammed glacial lakes. Significantly, we tie this reconstruction to work carried out and models developed for the entire ice sheet. This therefore situates research in the central sector within contemporary understanding of how the last BIIS evolved over time.

Landscape evolution and ice-sheet behaviour in a semi-arid polar environment: James Ross Island, NE Antarctic Peninsula

This study of landscape evolution presents both new modern and palaeo process-landform data, and analyses the behaviour of the Antarctic Peninsula Ice Sheet through the Last Glacial Maximum (LGM), the Holocene and to the present day. Six sediment-landform assemblages are described and interpreted for Ulu Peninsula, James Ross Island, NE Antarctic Peninsula: (1) the Glacier Ice and Snow Assemblage; (2) the Glacigenic Assemblage, which relates to LGM sediments and comprises both erratic-poor and erratic-rich drift, deposited by cold-based and wet-based ice and ice streams respectively; (3) the Boulder Train Assemblage, deposited during a Mid-Holocene glacier readvance; (4) the Ice-cored Moraine Assemblage, found in front of small cirque glaciers; (5) the Paraglacial Assemblage including scree, pebble-boulder lags, and littoral and fluvial processes; and (6) the Periglacial Assemblage including rock glaciers, protalus ramparts, blockfields, solifluction lobes and extensive patterned ground. The interplay between glacial, paraglacial and periglacial processes in this semi-arid polar environment is important in understanding polygenetic landforms. Crucially, cold-based ice was capable of sediment and landform genesis and modification. This landsystem model can aid the interpretation of past environments, but also provides new data to aid the reconstruction of the last ice sheet to overrun James Ross Island.

Timing and depositional environments of a Middle Pleistocene glaciation of northeast England: New evidence from Warren House Gill, County Durham

At various times during the Quaternary, north-eastern England was a zone of confluence between dynamic ice lobes sourced from the Pennines, northern Scotland, the Cheviots, and Scandinavia. The region thus has some of the most complex exposures of Middle to Late Pleistocene sediments in Britain, with both interglacial and glacial sediments deposited in terrestrial and marine settings. We investigated sedimentary sequences exposed on the coastline of County Durham at Warren House Gill, and present a new model of British and Fennoscandian Ice Sheet interaction in the North Sea Basin during the Middle Pleistocene. The stratigraphy at Warren House Gill consists of a lower diamicton and upper estuarine sediments, both part of the Warren House Formation. They are separated from the overlying Weichselian Blackhall and Horden tills by a substantial unconformity. The lower diamicton of the Warren House Formation is re-interpreted here as an MIS 8 to 12 glaciomarine deposit containing ice-rafted lithics from north-eastern Scotland and the northeast North Sea, and is renamed the ‘Ash Gill Member’. It is dated by lithological comparison to the Easington Raised Beach, Middle Pleistocene Amino Acid Racemisation values, and indirectly by optically stimulated luminescence. The overlying shallow subaqueous sediments were deposited in an estuarine environment by suspension settling and bottom current activity. They are named the ‘Whitesides Member’, and form the uppermost member of the Warren House Formation. During glaciation, ice-rafted material was deposited in a marine embayment. There is no evidence of a grounded, onshore Scandinavian ice sheet in County Durham during MIS 6, which has long been held as the accepted stratigraphy. This has major implications for the currently accepted British Quaternary Stratigraphy. Combined with recent work on the Middle Pleistocene North Sea Drift from Norfolk, which is now suggested to have been deposited by a Scottish ice sheet, the presence of a Scandinavian ice sheet in eastern England at any time during the Quaternary is becoming increasingly doubtful.

Dynamic Devensian ice flow in NE England: a sedimentological reconstruction

This study reconstructs the depositional environments that accompanied both ice advance and ice retreat of the last British–Irish Ice Sheet in NE England during the Last Glacial Maximum, and proposes three regional ice-flow phases. The Late Devensian (29–22 cal. ka BP) Tyne Gap Ice Stream initially deposited the Blackhall Till Formation during shelf-edge glaciation (Phase I). This subglacial traction till comprises several related facies, including stratified and laminated diamictons, tectonites, and sand and gravel beds deposited both in subglacial canals and in proglacial streams. Eventually, stagnation of the Tyne Gap Ice Stream led to ice-marginal sedimentation in County Durham (Phase II). During the Dimlington Stadial (21 cal. ka BP), the North Sea Lobe advanced towards the coastline of N Norfolk. This resulted initially in sandur deposition (widespread, tabular sand and gravel; the Peterlee Sand and Gravel Formation; Phase II) and ultimately in deposition of the Horden Till Formation (Phase III), a massive subglacial till. As the North Sea Lobe overrode previous formations, it thrusted and stacked sediments in County Durham, and dammed proglacial lakes between the east-coast ice, the Pennine uplands and the remaining Pennine ice. The North Sea Lobe retreated after Heinrich Event 1 (16 ka). This study highlights the complexity of ice flow during the Late Devensian glaciation of NE England, with changing environmental and oceanic conditions forcing a mobile and sensitive ice sheet.

Antarctic Peninsula ice sheet evolution during the Cenozoic Era

The Antarctic Peninsula region is currently undergoing rapid environmental change, resulting in the thinning, acceleration and recession of glaciers and the sequential collapse of ice shelves. It is important to view these changes in the context of long-term palaeoenvironmental complexity and to understand the key processes controlling ice sheet growth and recession. In addition, numerical ice sheet models require detailed geological data for tuning and testing. Therefore, this paper systematically and holistically reviews published geological evidence for Antarctic Peninsula Ice Sheet variability for each key locality throughout the Cenozoic, and brings together the prevailing consensus of the extent, character and behaviour of the glaciations of the Antarctic Peninsula region. Major contributions include a downloadable database of 186 terrestrial and marine calibrated dates; an original reconstruction of the LGM ice sheet; and a new series of isochrones detailing ice sheet retreat following the LGM. Glaciation of Antarctica was initiated around the Eocene/Oligocene transition in East Antarctica. Palaeogene records of Antarctic Peninsula glaciation are primarily restricted to King George Island, where glacigenic sediments provide a record of early East Antarctic glaciations, but with modification of far-travelled erratics by local South Shetland Island ice caps. Evidence for Neogene glaciation is derived primarily from King George Island and James Ross Island, where glaciovolcanic strata indicate that ice thicknesses reached 500–850 m during glacials. This suggests that the Antarctic Peninsula Ice Sheet draped, rather than drowned, the topography. Marine geophysical investigations indicate multiple ice sheet advances during this time. Seismic profiling of continental shelf-slope deposits indicates up to ten large advances of the Antarctic Peninsula Ice Sheet during the Early Pleistocene, when the ice sheet was dominated by 40 kyr cycles. Glacials became more pronounced, reaching the continental shelf edge, and of longer duration during the Middle Pleistocene. During the Late Pleistocene, repeated glacials reached the shelf edge, but ice shelves inhibited iceberg rafting. The Last Glacial Maximum (LGM) occurred at 18 ka BP, after which transitional glaciomarine sediments on the continental shelf indicate ice-sheet retreat. The continental shelf contains large bathymetric troughs, which were repeatedly occupied by large ice streams during Pleistocene glaciations. Retreat after the LGM was episodic in the Weddell Sea, with multiple readvances and changes in ice-flow direction, but rapid in the Bellingshausen Sea. The late Holocene Epoch was characterised by repeated fluctuations in palaeoenvironmental conditions, with associated glacial readvances. However, this has been subsumed by rapid warming and ice-shelf collapse during the twentieth century.

Late-Holocene changes in character and behaviour of land-terminating glaciers on James Ross Island, Antarctica

Virtually no information is available on the response of land-terminating Antarctic Peninsula glaciers to climate change on a centennial timescale. This paper analyses the topography, geomorphology and sedimentology of prominent moraines on James Ross Island, Antarctica, to determine geometric changes and to interpret glacier behaviour. The moraines are very likely due to a late-Holocene phase of advance and featured (1) shearing and thrusting within the snout, (2) shearing and deformation of basal sediment, (3) more supraglacial debris than at present and (4) short distances of sediment transport. Retreat of ∼100 m and thinning of 15–20 m has produced a loss of 0.1 km3 of ice. The pattern of surface lowering is asymmetric. These geometrical changes are suggested most simply to be due to a net negative mass balance caused by a drier climate. Comparisons of the moraines with the current glaciological surface structure of the glaciers permits speculation of a transition from a polythermal to a cold-based thermal regime. Small land-terminating glaciers in the northern Antarctic Peninsula region could be cooling despite a warming climate.

Palaeovegetation (Communications arising): diversity of temperate plants in east Asia

The exceptionally broad species diversity of vascular plant genera in east Asian temperate forests, compared with their sister taxa in North America, has been attributed to the greater climatic diversity of east Asia, combined with opportunities for allopatric speciation afforded by repeated fragmentation and coalescence of populations through Late Cenozoic ice-age cycles1. According to Qian and Ricklefs1, these opportunities occurred in east Asia because temperate forests extended across the continental shelf to link populations in China, Korea and Japan during glacial periods, whereas higher sea levels during interglacial periods isolated these regions and warmer temperatures restricted temperate taxa to disjunct refuges. However, palaeovegetation data from east Asia2, 3, 4, 5, 6 show that temperate forests were considerably less extensive than today during the Last Glacial Maximum, calling into question the coalescence of tree populations required by the hypothesis of Qian and Ricklefs1.

How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

Simulated climate and biomes of Africa during the Late Quaternary: comparison with pollen and lake status data

New compilations of African pollen and lake data are compared with climate (CCM1, NCAR, Boulder) and vegetation (BIOME 1.2, GSG, Lund) simulations for the last glacial maximum (LGM) and early to mid-Holocene (EMH). The simulated LGM climate was ca 4°C colder and drier than present, with maximum reduction in precipitation in semi-arid regions. Biome simulations show lowering of montane vegetation belts and expansion of southern xerophytic associations, but no change in the distribution of deserts and tropical rain forests. The lakes show LGM conditions similar or drier than present throughout northern and tropical Africa. Pollen data indicate lowering of montane vegetation belts, the stability of the Sahara, and a reduction of rain forest. The paleoenvironmental data are consistent with the simulated changes in temperature and moisture budgets, although they suggest the climate model underestimates equatorial aridity. EMH simulations show temperatures slightly less than present and increased monsoonal precipitation in the eastern Sahara and East Africa. Biome simulations show an upward shift of montane vegetation belts, fragmentation of xerophytic vegetation in southern Africa, and a major northward shift of the southern margin of the eastern Sahara. The lakes indicate conditions wetter than present across northern Africa. Pollen data show an upward shift of the montane forests, the northward shift of the southern margin of the Sahara, and a major extension of tropical rain forest. The lake and pollen data confirm monsoon expansion in eastern Africa, but the climate model fails to simulate the wet conditions in western Africa.

Climate and biome simulations for the past 21,000 years

This paper reports on a set of paleoclimate simulations for 21, 16, 14, 11 and 6 ka (thousands of years ago) carried out with the Community Climate Model, Version 1 (CCM1) of the National Center for Atmospheric Research (NCAR). This climate model uses four interactive components that were not available in our previous simulations with the NCAR CCM0 (COHMAP, 1988Science, 241, 1043–1052; Wright et al., 1993Global Climate Since the Last Glocial Maximum, University of Minnesota Press, MN): soil moisture, snow hydrology, sea-ice, and mixed-layer ocean temperature. The new simulations also use new estimates of ice sheet height and size from ( Peltier 1994, Science, 265, 195–201), and synchronize the astronomically dated orbital forcing with the ice sheet and atmospheric CO2 levels corrected from radiocarbon years to calendar years. The CCM1 simulations agree with the previous simulations in their most general characteristics. The 21 ka climate is cold and dry, in response to the presence of the ice sheets and lowered CO2 levels. The period 14–6 ka has strengthened northern summer monsoons and warm mid-latitude continental interiors in response to orbital changes. Regional differences between the CCM1 and CCM0 simulations can be traced to the effects of either the new interactive model components or the new boundary conditions. CCM1 simulates climate processes more realistically, but has additional degrees of freedom that can allow the model to ‘drift’ toward less realistic solutions in some instances. The CCM1 simulations are expressed in terms of equilibrium vegetation using BIOME 1, and indicate large shifts in biomes. Northern tundra and forest biomes are displaced southward at glacial maximum and subtropical deserts contract in the mid-Holocene when monsoons strengthen. These vegetation changes could, if simulated interactively, introduce additional climate feedbacks. The total area of vegetated land remains nearly constant through time because the exposure of continental shelves with lowered sea level largely compensates for the land covered by the expanded ice sheets.

Deglacial changes of the southern margin of the southern westerly winds revealed by terrestrial records from SW Patagonia (52 degrees S)

Much of the ongoing discussion regarding synchrony or bipolar asynchrony of paleoclimate events has centered on the timing and structure of the last glacial termination in the southern mid- latitudes, in particular the southwestern Patagonian region (50�e55�S). Its location adjacent to the Drake Passage andnear the southern margin of the southern westerly winds (SWW) allows examining the postulated links between the Southern Oceane SWW coupled system and tmospheric CO2 variations through the last glacial termination. Results from two sites located in the Última Esperanza area (52�S) allow us to infer SWW-driven changes in hydrologic balance during this critical time interval. These findings indicate peatland development under temperate/wet conditions between 14,600 and 14,900 cal yr BP, followed by cooling and a lake transgressive phase that led to a shallow lake during the early part of the Antarctic Cold Reversal (ACR, 13,600-14,600 cal yr BP), followed in turn by a deeper lake and modest warming during Younger Dryas time (YD, 11,800-13,000 cal yr BP), superseded by terrestrialization and forest expansion at the beginning of the Holocene. We propose that the SWW (i) strengthened and shifted northward during ACR time causing a precipitation rise in northwestern and southwestern Patagonia coeval with mid- and high-latitude cooling and a halt in the deglacial atmospheric CO2 rise; (ii) shifted southward during YD time causing a precipitation decline/increase in NW/SW Patagonia, respectively, high-latitude warming, and invigorated CO2 release from the Southern Ocean; (iii) became weaker between 10,000 and 11,500 cal yr BP causing a precipitation decline throughout Patagonia, concurrent with peak mid- and high-latitude temperatures and atmospheric CO2 concentrations.

The response of vegetation on the Andean flank in western Amazonia to Pleistocene climate change

A reconstruction of past environmental change from Ecuador reveals the response of lower montane forest on the Andean flank in western Amazonia to glacial-interglacial global climate change. Radiometric dating of volcanic ash indicates that deposition occurred ~324,000 to 193,000 years ago during parts of Marine Isotope Stages 9, 7, and 6. Fossil pollen and wood preserved within organic sediments suggest that the composition of the forest altered radically in response to glacial-interglacial climate change. The presence of Podocarpus macrofossils ~1000 meters below the lower limit of their modern distribution indicates a relative cooling of at least 5°C during glacials and persistence of wet conditions. Interglacial deposits contain thermophilic palms suggesting warm and wet climates. Hence, global temperature change can radically alter vegetation communities and biodiversity in this region.

Hydrology and climatology at Laguna La Gaiba, lowland Bolivia: complex responses to climatic forcings over the last 25,000 years

Diatom, geochemical and isotopic data provide a record of environmental change in Laguna La Gaiba, lowland Bolivia (17°450S, 57°350W), over the last ca. 25 000 years. High-resolution diatom analysis around the Last Glacial–Interglacial Transition provides new insights into this period of change. The full and late glacial lake was generally quite shallow, but with evidence of periodic flooding. At about 13 100 cal a BP, just before the start of the Younger Dryas chronozone, the diatoms indicate shallower water conditions, but there is a marked change at about 12 200 cal a BP indicating the onset of a period of high variability, with rising water levels punctuated by periodic drying. From ca. 11 800 to 10 000 cal a BP, stable, deeper water conditions persisted. There is evidence for drying in the early to middle Holocene, but not as pronounced as that reported from elsewhere in the southern hemisphere tropics of South America. This was followed by the onset of wetter conditions in the late Holocene consistent with insolation forcing. Conditions very similar to present were established about 2100 cal a BP. A complex response to both insolation forcing and millennial-scale events originating in the North Atlantic is noted.

Terrestrial biosphere changes over the last 120 kyr and their impact on ocean δ 13C

A new global synthesis and biomization of long (>40 kyr) pollen-data records is presented, and used with simulations from the HadCM3 and FAMOUS climate models to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Global modelled (BIOME4) biome distributions over time generally agree well with those inferred from pollen data. The two climate models show good agreement in global net primary productivity (NPP). NPP is strongly influenced by atmospheric carbon dioxide (CO2) concentrations through CO2 fertilization. The combined effects of modelled changes in vegetation and (via a simple model) soil carbon result in a global terrestrial carbon storage at the Last Glacial Maximum that is 210–470 Pg C less than in pre-industrial time. Without the contribution from exposed glacial continental shelves the reduction would be larger, 330–960 Pg C. Other intervals of low terrestrial carbon storage include stadial intervals at 108 and 85 ka BP, and between 60 and 65 ka BP during Marine Isotope Stage 4. Terrestrial carbon storage, determined by the balance of global NPP and decomposition, influences the stable carbon isotope composition (δ13C) of seawater because terrestrial organic carbon is depleted in 13C. Using a simple carbon-isotope mass balance equation we find agreement in trends between modelled ocean δ13C based on modelled land carbon storage, and palaeo-archives of ocean δ13C, confirming that terrestrial carbon storage variations may be important drivers of ocean δ13C changes.