Volume of Antarctic Ice at the Last Glacial Maximum, and its impact on global sea level change

The volume of Antarctic ice at the Last Glacial Maximum is a key factor for calculating the past contribution of melting ice sheets to Late Pleistocene global sea level change. At present, there are large uncertainties in our knowledge of the extent and thickness of the formerly expanded Antarctic ice sheets, and in the timing of their release as meltwater into the world’s oceans. This paper reviews the four main approaches to determining former Antarctic ice volume, namely glacial geology, glacio-isostatic studies, glaciological modelling, and ice core analysis and attempts to reconcile these to give a ‘best estimate’ for ice volume. In the Ross Sea there was a major expansion of grounded ice at the Last Glacial Maximum, accounting for 2.3–3.2 m of global sea level. At some time in the Weddell Sea a large grounded ice sheet corresponding to c. 2.7 m of global sea level extended to the shelf break. However, this ice expansion has not yet been confidently dated and may not relate to the Last Glacial Maximum. Around East Antarctica there was thickening and advance offshore of ice in coastal regions. Ice core evidence suggests that the interior of East Antarctica was either close to its present elevation or thinner during the last glacial so the effect of East Antarctica on sea level depends on the net balance between marginal thickening and interior thinning. Suggested East Antarctic contributions vary from a 3–5.5 m lowering to a 0.64 m rise in global sea level. The Antarctic Peninsula ice sheet thickened and extended offshore at the Last Glacial Maximum, with a sea level equivalent contribution of c. 1.7 m. Thus, the Antarctic ice sheets accounted for between 6.1 and 13.1 m of global sea level fall at the Last Glacial Maximum. This is substantially less than has been suggested by most previous studies but the maximum figure matches well with one modelling estimate. The timing of Antarctic deglaciation is not well known. In the Ross Sea, terrestrial evidence suggests deglaciation may have begun at c. 13,000 yr BP¹ but that grounded ice persisted until c. 6,500 yr BP. Marine evidence suggests the western Ross Sea was deglaciated by c. 11,500 yr BP. Deglaciation of the Weddell Sea is poorly constrained. Grounded ice in the northern Antarctic Peninsula had retreated by c. 13,000 yr BP, and further south deglaciation occurred sometime prior to c. 6,000 yr BP. Many parts of coastal East Antarctica apparently escaped glaciation at the LGM, but in those areas that were ice-covered deglaciation was underway by 10,000 yr BP. With existing data, the timing of deglaciation shows no firm relation to northern hemisphere-driven sea level rise. This is probably due partly to lack of Antarctic dating evidence but also to the combined influence of several forcing mechanisms acting during deglaciation.     

Antarctic Glacial Retreats and the Records of Sea Level Change in Deposits at the North Bank of the Shenzhen Bay,South China

In view of the huge ice cover of 24.5 X 10⁶ km³ in Antarctica, which accounts for over 90% by volume of the ice body on earth, the movement of Antarctic glaciers is a major control on global sea-level change and climatic fluctuation. As recorded in the Quaternary deposits in King George Island, West Antarctica, three rapid ablations can be recognized at 11000, 9000 and 6100 years ago and the global climate within the past 6000 years is characterized by small-amplitude warm-cold fluctuation. Intertidal deposits at the north bank of the Shenzhen Bay suggest a periodic variation in sea level in about every 670 years over the last 6000 years with low sea levels recorded in the periods of 5500-4900, 3900-3600, 2400-2200 and 1300-1200 years ago. Between these periods the sea level rised for about 80 cm on average. The modern warming climate in the last century corresponds with a rise of sea level at the rate of 2~3 mm/a. A state key project under the Antarctic Science & Technology Program (No. 905-02-04-03) and is partially financially supported by the Natural Science Foundation of Guangdong Province.     

Glacial rebound and sea-level change in the British Isles

Observations of sea levels around the coastline of the British Isles for the past 10,000–15,000 years exhibit a major regional variation and provide an important data base for testing models of glacial rebound as well as models of the Late Devensian ice sheet. A high-resolution rebound model has been developed which is consistent with both the spatial and temporal patterns of sea-level change and which demonstrates that the observations are the result of (i) the glacio-isostatic crustal rebound in response to the unloading of the ice sheet over Britain and, to a lesser degree, of the ice sheet over Fennoscandia, and (ii) the rise in sea-level from the melting Late Pleistocene ice sheets, including the response of the crust to the water loading (the hydro-isostatic effect). The agreement between model and observations is such that there is no need to invoke vertical crustal movements for Great Britain and Ireland of other than glacio-hydro-isostatic origin. The rebound contributions are important throughout the region and nowhere is it sufficiently small for the sea-level change to approximate the eustatic sea-level rise. The observational data distribution around the periphery as well as from sites near the centre of the former ice sheet is sufficient to permit constraints to be established on both earth model parameters specifying the mantle viscosity and lithospheric thickness and the extent and volume of the ice sheet at the time of the last glaciation. Preliminary solutions are presented which indicate an upper mantle viscosity of (3–5)10²⁰ Pas, a lithospheric thickness of about 100 km or less, and an ice model that was not confluent with the Scandinavian ice sheet during the last glaciation and whose maximum thickness over Scotland is unlikely to have exceeded about 1500 m.     

Holocene relative sea‐level changes in Harris,Outer Hebrides,Scotland, UK

Evidence for relative sea‐level changes during the middle and late Holocene is examined from two locations on the Atlantic coast of Harris, Outer Hebrides, Scotland, using morphological mapping and survey, stratigraphical, grain size and diatom analysis, and radiocarbon dating. The earliest event identified is a marine flood, which occurred after 7982–8348 cal. a (7370 ± 80 ¹⁴C a) BP, when the sea crossed a threshold lying at ?0.08 m Ordnance Datum Newlyn (OD) (?2.17 m mean high water springs (MHWS)) before withdrawing. This could have been due to a storm or to the Holocene Storegga Slide tsunami. By 6407–6122 cal. a (5500 ± 60 ¹⁴C a) BP, relative sea levels had begun to fall from a sandflat surface with an indicated MHWS level of between 0.08 and ?1.96 m (?2.01 to ?4.05 m). This fall reached between ?0.30 and ?2.35 m (?2.39 to ?4.44 m) after 5841–5050 cal. a (4760 ± 130 ¹⁴C a) BP, but was succeeded by a relative sea‐level rise which reached between 0.54 and ?1.57 m (?1.55 to ?3.66 m) by 5450–4861 cal. a (4500 ± 100 ¹⁴C a) BP. This rise continued, possibly with an interruption, until a second sandflat surface was reached between 2.34 and ?0.26 m (0.25 to ?2.35 m) between 2952–3375 cal. a (3000 ± 80 ¹⁴C a) and 1948–2325 cal. a (2130 ± 70 ¹⁴C a) BP, before present levels were reached. The regressive episode from the earliest sandflat is correlated with the abandonment of the Main Postglacial Shoreline. It is maintained that the fluctuations in relative sea level recorded can be correlated with similar events elsewhere on the periphery of the glacio‐isostatic centre and may therefore reflect secular changes in nearshore sea surface levels. Despite published evidence from trim lines of differential ice sheet loading across the area, no evidence of variations in uplift between the locations concerned could be found. Copyright © 2009 John Wiley & Sons, Ltd.     

Early Telychian changes in graptoloid diversity and sea level

It is suggested that a decrease in graptoloid diversity in the utilis Subzone (early Telychian) resulted from a eustatic fall in sea level induced by the spread of ice sheets from the highlands of central Gondwanaland. Major extinctions occurred within the retiolites and the genera Glyptograptus, Petalolithus, Stimulograptus and Rastrites. The selectivity of the extinction may reflect the degree of environmental tolerance of the various graptoloid taxa, with those exhibiting water mass specificity most seriously affected. Lithological changes, here correlated with the glacial interval, are recognized from the Yangtze Platform of China and from Europe.     

Numerical simulation of the paleohydrology of glacial Lake Oshkosh, eastern Wisconsin, USA

Proglacial lakes, formed during retreat of the Laurentide ice sheet, evolved quickly as outlets became ice-free and the earth deformed through glacial isostatic adjustment. With high-resolution digital elevation models (DEMs) and GIS methods, it is possible to reconstruct the evolution of surface hydrology. When a DEM deforms through time as predicted by our model of viscoelastic earth relaxation, the entire surface hydrologic system with its lakes, outlets, shorelines and rivers also evolves without requiring assumptions of outlet position. The method is applied to proglacial Lake Oshkosh in Wisconsin (13,600 to 12,900 cal yr BP). Comparison of predicted to observed shoreline tilt indicates the ice sheet was about 400 m thick over the Great Lakes region. During ice sheet recession, each of the five outlets are predicted to uplift more than 100 m and then subside approximately 30 m. At its maximum extent, Lake Oshkosh covered 6600 km² with a volume of 111 km³. Using the Hydrologic Engineering Center-River Analysis System model, flow velocities during glacial outburst floods up to 9 m/s and peak discharge of 140,000 m³/s are predicted, which could drain 33.5 km³ of lake water in 10 days and transport boulders up to 3 m in diameter.     

Cryosphere as a temporal sink and source of microplastics in the Arctic region

Microplastics (MPs) pollution has become a serious environmental issue of growing global concern due to the increasing plastic production and usage. Under climate warming, the cryosphere, defined as the part of Earth’s layer characterized by the low temperatures and the presence of frozen water, has been experiencing significant changes. The Arctic cryosphere (e.g., sea ice, snow cover, Greenland ice sheet, permafrost) can store and release pollutants into environments, making Arctic an important temporal sink and source of MPs. Here, we summarized the distributions of MPs in Arctic snow, sea ice, seawater, rivers, and sediments, to illustrate their potential sources, transport pathways, storage and release, and possible effects in this sentinel region. Items concentrations of MPs in snow and ice varied about 1–6 orders of magnitude in different regions, which were mostly attributed to the different sampling and measurement methods, and potential sources of MPs. MPs concentrations from Arctic seawater, river/lake water, and sediments also fluctuated largely, ranging from several items of per unit to >40,000 items m^?3, 100 items m^?3, and 10,000 items kg^?1 dw, respectively. Arctic land snow cover can be a temporal storage of MPs, with MPs deposition flux of about (4.9–14.26) × 10⁸ items km^?2 yr^?1. MPs transported by rivers to Arctic ocean was estimated to be approximately 8–48 ton/yr, with discharge flux of MPs at about (1.65–9.35) × 10⁸ items/s. Average storage of MPs in sea ice was estimated to be about 6.1×10¹⁸ items, with annual release of about 5.1×10¹⁸ items. Atmospheric transport of MPs from long-distance terrestrial sources contributed significantly to MPs deposition in Arctic land snow cover, sea ice and oceanic surface waters. Arctic Great Rivers can flow MPs into the Arctic Ocean. Sea ice can temporally store, transport and then release MPs in the surrounded environment. Ocean currents from the Atlantic brought high concentrations of MPs into the Arctic. However, there existed large uncertainties of estimation on the storage and release of MPs in Arctic cryosphere owing to the hypothesis of average MPs concentrations. Meanwhile, representatives of MPs data across the large Arctic region should be mutually verified with in situ observations and modeling. Therefore, we suggested that systematic monitoring MPs in the Arctic cryosphere, potential threats on Arctic ecosystems, and the carbon cycle under increasing Arctic warming, are urgently needed to be studied in future.     

Should Quaternary sea-level changes be used to correct glacier ELAs, vegetation belt altitudes and sea level temperatures for inferring climate changes?

Changes in the altitudes of glacier snowlines (ELAs) and the altitudes of montane vegetation belts (VBAs) measure Quaternary climatic change. An accepted ‘correction’ to such changes by deducting the amount of contemporary sea level fall is wrong, since the air displaced by the ice sheets approximately fills the space left by the falling sea level and so there is no overall downward movement of the troposphere. This also causes a reduced cooling at the lowered sea level relative to that at the former inter-glacial sea level, about 1°C at the Las Glacial Maximum, which reduces the discrepancies previously noted by others between terrestrial and marine estimates of sea-level cooling. The change in temperature is indicated by the product of the ELA or VBA lowering and the environmental lapse rate (ELR). Prior estimates of ΔELA (−900 ± 135 m) and ELR (−6° ± 0.1°C km⁻¹) would indicate a cooling of −5.4°C at interglacial sea level and −4.4°C at glacial sea level, although glacial-period ELRs are not known reliably. Established ELA corrections for local epeirogenic uplift or subsidence are appropriate.     

Mid‐late Holocene sea‐level variability in eastern Australia

A re‐analysis of sea‐level data from eastern Australia based on 115 calibrated C‐14 ages is used to constrain the origin, timing and magnitude of sea‐level change over the last 7000 years. We demonstrate that the Holocene sea‐level highstand of +1.0–1.5 m was reached ～7000 cal yr bp and fell to its present position after 2000 yr bp . These findings are in contrast to most previous studies that relied on smaller datasets and did not include the now common conversion of conventional C‐14 ages to calendar years. During this ～5000 year period of high sea level, growth hiatuses in oyster beds and tubeworms and lower elevations of coral microatolls are interpreted to represent short‐lived oscillations in sea‐level of up to 1 m during two intervals, beginning c. 4800 and 3000 cal yr bp . The rates of sea‐level rise and fall (1–2 mm yr^?1) during these centennial‐scale oscillations are comparable with current rates of sea‐level rise. The origin of the oscillations is enigmatic but most likely the result of oceanographic and climatic changes, including wind strengths, ice ablation, and melt‐water contributions of both Greenland and Antarctic ice sheets.     

The origin and flux of icebergs released into the Last Glacial Maximum Northern Hemisphere oceans: the impact of ice‐sheet topography

The precipitation fields of a palaeoatmospheric general circulation model are used to derive estimates of the geographical distribution, and flux, of icebergs from the Laurentide, Fennoscandinavian and eastern Siberian ice‐sheets at the Last Glacial Maximum (LGM). The atmospheric model fields from LGM simulations using CLIMAP or Peltier (ICE‐4G) ice orography were studied, to test the sensitivity of the predicted flux. The estimated Northern Hemispheric LGM iceberg flux is 3500–4000 km³ yr^?1, of which about 60% issued directly into the North Atlantic. The iceberg flux from the St Lawrence area is of similar significance to that issuing from Hudson Strait in all estimates. Both the North Pacific and the Arctic received substantial iceberg fluxes (ca. 700 km³ yr^?1), with relatively minor differences occurring between the two ice‐sheet reconstructions. Apparent discrepancies between Arctic deep‐sea core samples of ice‐rafted debris and our estimates of mean glacial iceberg flux may be ascribed to coastal trapping of bergs, the existence of floating ice tongues or a rapid exit of icebergs from the Arctic basin into the Greenland Sea through the Fram Strait. Copyright © 2001 John Wiley & Sons, Ltd.     

On the reconstruction of pleistocene ice sheets: A review

Pleistocene ice sheets can be reconstructed through three separate approaches: (1) Evidence based on glacial geological studies, such as erratic trains, till composition, crossing striations and exposures of multiple tills/nonglacial sediments. (2) Reconstructions based on glaciological theory and observations. These can be either two- or three-dimensional models; they can be constrained by ‘known’ ice margins at specific times; or they can be ‘open-ended’ with the history of growth and retreat controlled by parameters resting entirely within the model. (3) Glacial isostatic rebound after deglaciation provides a measure of the distribution of mass (ice) across a region. A ‘best fit’ ice sheet model can be developed that closely approximates a series of relative sea level curves within an area of a former ice sheet; in addition, the model should also provide a reasonable sea level fit to relative sea level curves at sites well removed from glaciation.This paper reviews some of the results of a variety of ice sheet reconstructions and concentrates on the various attempts to reconstruct the ice sheets of the last (Wisconsin, Weischelian, Würm, Devensian) glaciation. Evidence from glacial geology suggests flow patterns at variance with simple, single-domed ice sheets over North America and Europe. In addition, reconstruction of ice sheets from glacial isostatic sea level data suggests that the ice sheets were significantly thinner than estimates based on 18 ka equilibrium ice sheets (cf. Denton and Hughes, 1981). The review indicates it is important to differentiate between ice divides, which control the directions of glacial flow, and areas of maximum ice thickness, which control the glacial isostatic rebound of the crust upon deglaciation. Recent studies from the Laurentide Ice Sheet region indicate that the center of mass was not over Hudson Bay; that a major ice divide lay east of Hudson Bay so that flow across the Hudson Bay and James Bay lowlands was from the northeast; that Hudson Bay was probably open to marine invasions two or three times during the Wisconsin Glaciation; and that the Laurentide Ice Sheet was thinner than an equilibrium reconstruction would suggest.     

Preliminary 10Be chronology for the last deglaciation of the western margin of the Greenland Ice Sheet

The now acknowledged thinning of the Greenland Ice Sheet raises concerns about its potential contribution to future sea level rise. In order to appreciate the full extent of its contribution to sea level rise, reconstruction of the ice sheet's most recent last deglaciation could provide key information on the timing and the height of the ice sheet at a time of rapid climate readjustment. We measured ¹⁰Be concentrations in 12 samples collected along longitudinal and altitudinal transects from Sisimiut to within 10 km of the Isunguata Sermia Glacier ice margin on the western coast of Greenland. Along the longitudinal transect, we collected three perched boulders and two bedrocks. In addition, we sampled seven perched boulders along a vertical transect in a valley within 10 km of the Isunguata Sermia Glacier ice margin. Our pilot dataset constrains the height of the ice sheet during the Last Glacial Maximum (LGM) between 500 m and 840 m (including the 120 m relative sea level depression at the time of the LGM, 21 ka BP). From the transect we estimate the thinning of the ice sheet at the end of the deglaciation between 12.3 ± 1.5 ¹⁰Be ka (n = 2) and 8.3 ± 1.2 ¹⁰Be ka (n = 3) to be ～6 cm a^?1 over this time period. Direct dating of the retreat of the western margin of the Greenland Ice Sheet has the potential to better constrain the retreat rate of the ice margin, the thickness of the former ice sheet as well as its response to climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

Coastal Wetlands of China: Changes from the 1970s to 2007 Based on a New Wetland Classification System

A new classification of coastal wetlands along the coast of China has been generated that is compatible with the Ramsar Convention of 1971. The coastal wetlands have been divided into two broad categories with overall nine subcategories. On this basis, a series of coastal wetland maps, together covering the coast of mainland China, have been produced based on topographic maps acquired in the 1970s and satellite images acquired in 2007. These document substantial wetland losses over this period. In the 1970s, the total coastal wetland area in China was 5.76?×?10⁴?km², whereas in 2007, it was 5.36?×?10⁴?km², indicating a loss of 7 %. Over this approximately 40-year period, the area of natural coastal wetlands decreased from 5.74?×?10⁴ to 5.09?×?10⁴?km², while that of artificial coastal wetlands increased from 240 to 2,740 km². Due to shoreline and sea-level changes, newly formed coastal wetlands amounted to 2,460 km², while coastal wetland loss amounted to 6,310 km² in the period from the 1970s to 2007. When excluding shallow coastal waters (depths between 0 and ?5 m), nearly 16 % of Chinese coastal wetlands have been lost between the 1970s and 2007.     

Late Holocene relative sea level rise and the Neoglacial history of the Greenland ice sheet

In West Greenland, early and mid Holocene relative sea level (RSL) fall was replaced by late Holocene RSL rise during the Neoglacial, after 4–3 cal. ka BP (thousand calibrated years before present). Here we present the results of an isolation basin RSL study completed near to the coastal town of Sisimiut, in central West Greenland. RSL fell from 14 m above sea level at 5.7 cal. ka BP to reach a lowstand of ?4.0 m at 2.3–1.2 cal. ka BP, before rising by an equivalent amount to present. Differences in the timing and magnitude of the RSL lowstand between this and other sites in West and South Greenland record the varied interplay of local and non‐Greenland RSL processes, notably the reloading of the Earth's crust caused by a Neoglacial expansion of the Greenland Ice Sheet (GIS) and the subsidence associated with the collapse of the Laurentide Ice Sheet forebulge. This means that the timing of the sea level lowstand cannot be used to infer directly when the GIS advanced during the Neoglacial. The rise in Late Holocene RSL is contrary to recently reported bedrock uplift in the Sisimiut area, based on repeat GPS surveys. This indicates that a belt of peripheral subsidence around the current ice sheet margin was more extensive in the late Holocene, and that there has been a switch from subsidence to uplift at some point in the last thousand years or so. Copyright © 2008 John Wiley & Sons, Ltd.     

青藏高原东部第四纪冰川问题

本文讨论青藏高原东部地区第四纪古冰川遗迹,根据野外实地调查结合卫星影象和航空照片资料,说明末次冰期及倒数第二次冰期中古冰川从未联成M.Kuhle 等所称的大冰盖。各种证据说明,倒数第二次冰期降温最大,冰川、冰缘及河流加积作用明显,相当于深海氧同位素阶段6。在此之前有过漫长的大间冰期,致使倒数第三次冰期冰碛被侵蚀凌夷,表面发育红色风化壳。     

Ice‐distal landscape and sediment signatures evidencing damming and drainage of large pro‐glacial lakes,northwest Russia

Lyså, A., Jensen, M. A., Larsen, E., Fredin, O. & Demidov, I. N.^* 2010: Ice‐distal landscape and sediment signatures evidencing damming and drainage of large pro‐glacial lakes, northwest Russia. Boreas, Vol. 40, pp. 481–497. 10.1111/j.1502‐3885.2010.00197.x. ISSN 0300‐9483. Sediments from river sections and the morphology of the upper reaches of Severnaya Dvina and Vychegda in northwest Russia show evidence of the existence of large ice‐dammed lakes in the area twice during the Weichselian. During the Late Weichselian, three separate ice‐dammed lakes (LGM lake(s)) existed, the largest one at about 135 m a.s.l. having a volume of about 1510 km³. Stepwise and rapid lake drainage is suggested to have taken place within less than 1000 years. The locations of various passpoints controlled the drainage, and when the lake was at its maximum level water spilled southeastwards into the Volga basin. Later, but before the lake water finally drained into the White Sea, water was routed northeastwards into the southeastern part of the Barents Sea. The oldest lake, the White Sea lake, existed around 67–57 ka ago, slightly in conflict with earlier palaeogeographic reconstructions regarding the chronology. The extent of the lake was constrained by, in addition to the Barents Sea ice‐sheet margin in the north, thresholds in the drainage basin. Later, one threshold was eroded and lowered during the LGM lake drainage. Given a lake level of about 115 m a.s.l., a lake area of about 2.5 × 10⁴ km³ and a water volume of about 4800 km³, the lake drainage northwards and into the ocean probably impacted the ocean circulation.     

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.     

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.     

Shingled Quaternary debris flow lenses on the north-east Newfoundland Slope

Debris flow deposits are the principal component of Quaternary continental slope sediments between the north-east Newfoundland Shelf and central Orphan Basin. In seismic profiles, these deposits occur as shingled, elongate, acoustically transparent lenses with their long axes orientated downslope. Deposits of individual flows form positive mounds on the sea floor; subsequent flows were diverted by the pre-existing topography into bathymetric lows between older debris flow deposits. These deposits show a large variation in the area of sea floor covered by individual flows (about 60–1000 km²), average thickness of deposits (9–37 m) and volume of sediment displaced (1–27 km³). The ratio of average thickness to a measure of deposit diameter, termed the aspect ratio, has a threefold variation from 0·0006 to 0·0021. Very low depositional slopes and low aspect ratios suggest relatively low viscosities, probably due to inmixing of water during downslope transport. Stratified sediments form three distinct horizons and are locally interbedded with the debris flow deposits. These are mainly hemipelagic deposits. The slope and rise to the west of the Orphan Basin are constructional in character. The apparent absence of upper slope erosional features and the abundance of debris flow deposits on the slope suggest that the supply of sediment to the continental slope occurred predominantly during times of maximum extent of Quaternary glacial ice. The ice sheet grounding line during several glacial maxima must have been situated at or near the present shelf break, supplying vast amounts of sediment directly to the upper slope. Oversteepening and subsequent slope failures fed material into deeper water.     

Modeling of pleistocene European ice sheets: Some experiments with simple mass-balance parameterizations

A vertically integrated ice-flow model suitable for use in climate studies is formulated. Large continental ice sheets may be characterized by two fundamental quantities: the height-to-width ratio, and the steepness of the edge. So it is natural to develop a model containing two parameters that can be chosen to give the right values of those characteristic quantities. The result is a model that is close to M. A. W. Mahaffy's (Journal of Geophysical Research, 81, 1059–1066 (1976)). The model is used to study glaciation in Europe. Dropping the level of zero mass balance creates small stable ice caps in the Alps and the Scandinavian mountains. If the drop exceeds 600 m (with respect to present-day conditions), the feedback between ice-sheet height and mass balance becomes dominating and the Fennoscandian Ice Sheet keeps growing. It does not reach an equilibrium state within 60,000 yr. An experiment simulating rapid onset of a glacial cycle shows that the growth of ice volume in Europe is smaller than that in northern America (J. T. Andrews and M. A. W. Mahaffy, Quaternary Research, 6, 167–183 (1976)). After 10,000 yr, the volume of the Fennoscandian Ice Sheet (2 × 10¹⁵ m³) is about half the volume of the Laurentide Ice Sheet. This leaves the “observed” sea-level lowering in the period 125,000–115,000 yr B.P. (estimates center around 50 m) unexplained.