Glaciology of the British Isles Ice Sheet during the last glacial cycle: form, flow, streams and lobes

The last ice sheet over the British Isles, together with other mid-latitude Pleistocene ice sheets, and in contrast to the modern ice sheets of Greenland and Antarctica, had a relatively low profile, low summit elevation and extensive, elongated lobes at its margin. A thermo-mechanically coupled numerical ice sheet model, driven by a proxy climate, has been used to explore the properties that would have permitted these characteristics to develop. The approach, the key to quantitative palaeoglaciology, is to determine the boundary conditions that permit the simulated ice sheet to mimic the evolution of the real ice sheet through the last glacial cycle. Simulations show how a British ice sheet may have been confluent with a Scandinavian ice sheet during some parts of its history and how unforced periodic and asynchronous oscillations could occur in different parts of its margins. Marginal lobes are a reflection of streaming within the ice sheet. Such streams can be ephemeral, dynamic streams located because of ice sheet properties, or fixed streams whose location is determined by bed properties. The simulations that best satisfy constraints of extent, elevation and relative sea levels are those with major fixed streams that strongly draw down the ice sheet surface. In these, the core upland areas of the ice sheet were cold based at the Last Glacial Maximum, basal streaming velocities were between 500 and 1000 ma⁻¹ compared with surface velocities of 10–50 ma⁻¹ in inter-stream zones, shear stresses were as low as 15–25 kPa in streams compared with 70–110 kPa in upland areas and 60–84% of the ice flux was delivered to the margin via streams.     

Paleoglaciological reconstructions for the Tibetan Plateau during the last glacial cycle: evaluating numerical ice sheet simulations driven by GCM-ensembles

The Tibetan Plateau is a topographic feature of extraordinary dimension and has an important impact on regional and global climate. However, the glacial history of the Tibetan Plateau is more poorly constrained than that of most other formerly glaciated regions such as in North America and Eurasia. On the basis of some field evidence it has been hypothesized that the Tibetan Plateau was covered by an ice sheet during the Last Glacial Maximum (LGM). Abundant field- and chronological evidence for a predominance of local valley glaciation during the past 300,000 calendar years (that is, 300 ka), coupled to an absence of glacial landforms and sediments in extensive areas of the plateau, now refute this concept. This, furthermore, calls into question previous ice sheet modeling attempts which generally arrive at ice volumes considerably larger than allowed for by field evidence. Surprisingly, the robustness of such numerical ice sheet model results has not been widely queried, despite potentially important climate ramifications. We simulated the growth and decay of ice on the Tibetan Plateau during the last 125 ka in response to a large ensemble of climate forcings (90 members) derived from Global Circulation Models (GCMs), using a similar 3D thermomechanical ice sheet model as employed in previous studies. The numerical results include as extreme end members as an ice-free Tibetan Plateau and a plateau-scale ice sheet comparable, in volume, to the contemporary Greenland ice sheet. We further demonstrate that numerical simulations that acceptably conform to published reconstructions of Quaternary ice extent on the Tibetan Plateau cannot be achieved with the employed stand-alone ice sheet model when merely forced by paleoclimates derived from currently available GCMs. Progress is, however, expected if future investigations employ ice sheet models with higher resolution, bidirectional ice sheet-atmosphere feedbacks, improved treatment of the surface mass balance, and regional climate data and climate reconstructions.     

南极冰盖研究最新进展

南极冰盖是地球系统的重要组成部分,在全球气候系统中扮演着重要角色.通过对南极冰盖的研究将有助于了解其在全球气候系统中的作用,并为探讨全球气候过去、现在以及未来的演化提供支撑.总结分析了近年来南极冰盖研究的一些重要进展,并在此基础上对南极冰盖研究领域的一些主要结果、观测事实以及未来变化展开讨论,重点介绍南极物质平衡、冰芯研究、冰下水系统、冰盖数值模拟方面最近的进展,评述未来可能的研究方向和应该关注的问题.     

A simple holistic hypothesis for the self-destruction of ice sheets

Ice sheets are the only components of Earth’s climate system that can self-destruct. This paper presents the quantitative force balance for bottom-up modeling of ice sheets, as first presented qualitatively in this journal as a way to quantify ice-bed uncoupling leading to self-destruction of ice sheets (Hughes, 2009a). Rapid changes in sea level and climate can result if a large ice-sheet self-destructs quickly, as did the former Laurentide Ice Sheet of North America between 8100 and 7900 BP, thereby terminating the last cycle of Quaternary glaciation. Ice streams discharge up to 90 percent of ice from past and present ice sheets. A hypothesis is presented in which self-destruction of an ice sheet begins when ubiquitous ice-bed decoupling, quantified as a floating fraction of ice, proceeds along ice streams. This causes ice streams to surge and reduce thickness by some 90 percent, and height above sea level by up to 99 percent for floating ice, so the ice sheet undergoes gravitational collapse. Ice collapsing over marine embayments becomes floating ice shelves that may then disintegrate rapidly. This floods the world ocean with icebergs that reduce the ocean-to-atmosphere heat exchange, thereby triggering climate change. Calving bays migrate up low stagnating ice streams and carve out the accumulation zone of the collapsed ice sheet, which prevents its recovery, decreases Earth’s albedo, and terminates the glaciation cycle. This sequence of events may coincide with a proposed life cycle of ice streams that drain the ice sheet. A first-order treatment of these life cycles is presented that depends on the longitudinal force balance along the flowbands of ice streams and gives a first approximation to ice-bed uncoupling at snapshots during gravitational collapse into ice shelves that disintegrate, thereby removing the ice sheet. The stability of the Antarctic Ice Sheet is assessed using this bottom-up approach.     

Reconstructing the last Irish Ice Sheet 2: a geomorphologically-driven model of ice sheet growth,retreat and dynamics

The ice sheet that once covered Ireland has a long history of investigation. Much prior work focussed on localised evidence-based reconstructions and ice-marginal dynamics and chronologies, with less attention paid to an ice sheet wide view of the first order properties of the ice sheet: centres of mass, ice divide structure, ice flow geometry and behaviour and changes thereof. In this paper we focus on the latter aspect and use our new, countrywide glacial geomorphological mapping of the Irish landscape (>39 000 landforms), and our analysis of the palaeo-glaciological significance of observed landform assemblages (article Part 1), to build an ice sheet reconstruction yielding these fundamental ice sheet properties. We present a seven stage model of ice sheet evolution, from initiation to demise, in the form of palaeo-geographic maps. An early incursion of ice from Scotland likely coalesced with local ice caps and spread in a south-westerly direction 200 km across Ireland. A semi-independent Irish Ice Sheet was then established during ice sheet growth, with a branching ice divide structure whose main axis migrated up to 140 km from the west coast towards the east. Ice stream systems converging on Donegal Bay in the west and funnelling through the North Channel and Irish Sea Basin in the east emerge as major flow components of the maximum stages of glaciation. Ice cover is reconstructed as extending to the continental shelf break. The Irish Ice Sheet became autonomous (i.e. separate from the British Ice Sheet) during deglaciation and fragmented into multiple ice masses, each decaying towards the west. Final sites of demise were likely over the mountains of Donegal, Leitrim and Connemara. Patterns of growth and decay of the ice sheet are shown to be radically different: asynchronous and asymmetric in both spatial and temporal domains. We implicate collapse of the ice stream system in the North Channel – Irish Sea Basin in driving such asymmetry, since rapid collapse would sever the ties between the British and Irish Ice Sheets and drive flow configuration changes in response. Enhanced calving and flow acceleration in response to rising relative sea level is speculated to have undermined the integrity of the ice stream system, precipitating its collapse and driving the reconstructed pattern of ice sheet evolution.     

Influence of glacier hydrology on the dynamics of a large Quaternary ice sheet

The influence of glacier hydrology on the time-dependent morphology and flow behaviour of the late Weichselian Scandinavian ice sheet is explored using a simple one-dimensional ice sheet model. The model is driven by orbitally induced radiation variations, ice-albedo feedback and eustatic sea-level change. The influence of hydrology is most marked during deglaciation and on the southern side of the ice sheet, where a marginal zone of rapid sliding, thin ice and low surface slopes develops. Such a zone is absent when hydrology is omitted from the model, and its formation results in earlier and more rapid deglaciation than occurs in the no-hydrology model. The final advance to the glacial maximum position results from an increase in the rate of basal sliding as climate warms after 23000 yr BP. Channelised subglacial drainage develops only episodically, and is associated with relatively low meltwater discharges and high hydraulic gradients. The predominance of iceberg calving as an ablation mechanism on the northern side of the ice sheet restricts the occurrence of surface melting. Lack of meltwater penetration to the glacier bed in this area means that ice flow is predominantly by internal deformation and the ice sheet adopts a classical parabolic surface profile.     

中国极地冰盖数值模拟研究进展北大核心CSCD

冰盖数值模拟是一种基于多源观测数据,通过构建并求解冰流动力学方程组,理解冰流运动物理机制以及诊断和预估其演化过程的方法,目前已被广泛应用于冰盖变化研究。本文简要介绍了极地冰盖数值模拟方法,归纳综述了近十余年我国学者在极地冰盖数值模拟方面的研究进展,厘清我国在冰盖数值模拟领域遇到的瓶颈和关键问题。阐述了如何与我国的极地冰盖科考优势区域深度结合,协同多源强化观测和数值模拟,研发和改进冰盖模式,提高冰盖模拟能力,对定量估算极地冰盖的物质平衡及其对未来海平面上升的影响做出实质贡献。通过逐步发展冰盖模式的研究能力,有望将来在冰盖关键动力过程和机制的科学认识上有所突破。     

天山托木尔峰科其喀尔巴西冰川表面运动速度特征分析

天山托木尔峰科其喀尔巴西冰川是典型的树枝状山谷冰川,利用3组(6期)ASTER遥感影像通过COSI-corr软件反演了该冰川表面运动速度.与花杆测量数据进行对比,反演冰川表面运动速度平均绝对误差为3.1 m·a-1,相对误差为11.9%,二者在空间上的分布基本一致,表明其反演精度符合要求.在此基础上,分析冰川表面运动速...     

Analytical estimate of the critical global-warming level for the Antarctic ice sheet mass gain-to-loss transition

We suggested a relatively simple model describing changes in the total balance of the ice sheet mass due to global climate change. Taking into account the basic mechanisms behind change in the ice sheet mass and their relations with temperature, we obtained a nonlinear analytical dependence of the ice sheet thickness on the global near-surface temperature. The behavior of the ice sheet can be split into six regimes. Implementation of some regime or another depends essentially on the initial (present-day) value of the mass balance and the sensitivity parameters of precipitation and ice melt to the temperature. Based on this model, we obtained an analytical estimate of the critical level in global warming, in excess of which the regime of the Antarctic ice sheet gain due to snow accumulation changes to sheet degradation due to more intense growth in ice melting.     

Modelling ice-sheet sensitivity to late weichselian environments in the svalbard-barents sea region

Ice-proximal sedimentological features from the northwestern Barents Sea suggest that this region was covered by a grounded ice sheet during the Late Weichselian. However, there is debate as to whether these sediments were deposited by the ice sheet at its maximum or a retreating ice sheet that had covered the whole Barents Sea. To examine the likelihood of total glaciation of the Late Weichselian Barents Sea, a numerical ice-sheet model was run using a range of environmental conditions. Total glaciation of the Barents Sea, originating solely from Svalbard and the northwestern Barents Sea, was not predicted even under extreme environmental conditions. Therefore, if the Barents Sea was completely covered by a grounded Late Weichselian ice sheet, then a mechanism (not accounted for within the glaciological model) by which grounded ice could have formed rapidly within the central Barents Sea, may have been active during the last glaciation. Such mechanisms include (i) grounded ice migration from nearby ice sheets in Scandinavia and the central Barents Sea, (ii) the processes of sea-ice-induced ice-shelf thickening and (iii) isostatic uplift of the central Barents Sea floor.     

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.     

Constraints on the Late Weichselian ice sheet over the Barents Sea from observations of raised shorelines

Direct evidence for Late Weichselian grounded glacier ice over extensive areas of the Barents Sea is based largely on indirect observations, including elevations of old shorelines on Svalbard and arguments of isostatic rebound. Such isostatic models are discussed here for two cases representing maximum and minimum ice-sheet reconstructions. In the former model the ice extends over the Kara Sea, whereas in the latter the ice is limited to the Barents Sea and island archipelagos. Comparisons of predictions with observations from a number of areas, including Spitsbergen, Nordaustlandet, Edgeøya, Kong Karls Land, Franz Josef Land, Novaya Zemlya and Finnmark, support arguments for the existence of a large ice sheet over the region at the time of the last glacial maximum. This ice sheet is likely to have had the following characteristics, conclusions that are independent of assumptions made about the Earth's rheological parameters. (i) The maximum thickness of this ice was about 1500–2000 m with the centre of the load occurring to the south and east of Kong Karls Land. (ii) The ice sheet extended out to the western edge of the continental shelf and its maximum thickness over western Spitsbergen was about 800 m. (iii) To the north of Svalberg and Frans Josef Land the ice sheet extended out to the northern shelf edge. (iv) Retreat of the grounded ice across the southern Barents Sea occurred relatively early such that this region was largely ice free by about 15,000 BP. (v) By 12,000 BP the grounded ice had retreated to the northern archipelagos and was largely gone by 10,000 BP. (vi) The ice sheet may have extended to the Kara Sea but ice thicknesses were only a fraction of those proposed in those reconstructions where the maximum ice thickness is centered on Novaya Zemlya. Models for the palaeobathymetry for the Barents Sea at the time of the last glacial maximum indicate that large parts of the Barents Sea were either very shallow or above sea level, providing the opportunity for ice growth on the emerged plateaux, as well as on the islands, but only towards the end of the period of Fennoscandian ice sheet build-up.     

基于海冰物理和力学参数的渤海冰工程分区

根据渤海海冰物理和力学参数分布建立渤海海冰工程分区.其特点是将作用到相同直径或宽度的直立结构物上冰的作用力范围相同的区域划分为同一分区内,因此不同重现期设计冰厚度和峰值压缩强度之积作为区划指标.为了获得不同重现期渤海海冰物理和力学性质参数,利用渤海历史资料建立用气温、水温、冰期和冰厚表达非变形平整冰层和冰样温度、盐度、密度的关系,进而获得冰的孔隙率;利用渤海海冰单轴压缩强度同孔隙率的试验关系计算相应峰值压缩强度;然后发展渤海不同重现期平整冰层物理和力学参数的平面分布.在此基础上以新区划指标做不同重现期的海冰工程分区.作为新区划的实例,给出25、50、100年重现期的渤海海冰工程分区.     

Late devensian ice sheet in the western Grampians,Scotland

A radial pattern of ice flow of the last ice sheet in the largest source area of ice (c. 5000 km²) in the British Isles in western Scotland is demonstrated by the dispersal of indicator erratics and by patterns of striae, friction cracks and ice-moulded landforms. Three major ice domes and the principal ice divide are identified in the western Grampians. The ice domes coincided with the highest mountain blocks while it is inferred that the alignment of many of the pre-existing valleys controlled much of the outflow of ice, forming ice streams within the ice sheet. The importance of the Rannoch Moor basin as a radial provider of ice to surrounding areas was apparently less significant than has hitherto been considered. Ice flowed into and across the southern part of the basin from the principal ice divide located to the west and south. No evidence has yet been found that would support a model of an eastward-migrating ice divide either during the build-up or during the deglaciation of the Late Devensian ice sheet.     

Dynamics of the Late Weichselian ice sheet on Svalbard inferred from high-resolution sea-floor morphology

High-resolution bathymetric mapping of the fjords and continental shelf around the Svalbard archipelago shows an extensive pattern of large- and medium-scale submarine landforms formed by differences in ice-flow regimes. Mega-scale glacial lineations, lateral moraines, transverse ridges and glaciotectonic features are superimposed on the large-scale fjord, shelf and cross-shelf trough morphology of the margin. From these landforms we have inferred the flow and dynamics of the last ice sheet on Svalbard. Major fjords and their adjacent cross-shelf troughs have been identified as the main routes for ice streams draining the ice sheet. On the west coast of Svalbard major pathways existed along Bellsund, Isfjorden and Kongsfjorden. Along the northern Svalbard margin most of the ice drained through the Woodfjorden cross-shelf trough and Wijdefjorden-Hinlopen strait. Extensive areas with trough-parallel glacial lineations in the cross-shelf troughs suggest fast ice flow by palaeo-ice streams. Lateral ice-stream moraines, several tens of kilometres in length, have been mapped along the margins of some of the cross-shelf troughs, identifying the border zone between fast ice flow and stagnant or slow-flowing ice on intervening banks. Several general implications can be drawn from the interpretation of the glacier-derived submarine landforms around Svalbard. Firstly, the Late Weichselian ice sheet was partitioned into fast-flowing ice streams separated by slower moving ice. Secondly, our submarine morphological evidence supports earlier sedimentological, stratigraphical and chronological studies in implying that a large ice sheet reached the shelf edge around almost all of western and northern Svalbard in the Late Weichselian. The idea of a relatively restricted ice sheet over Svalbard, with ice-free conditions in some areas of the west coast at the Last Glacial Maximum, is therefore unlikely to be correct. Thirdly, the ice sheet appears to have retreated more rapidly from the cross-shelf troughs and outer fjords, although sometimes this occurred in a punctuated pattern indicated by grounding-zone wedges, and more slowly from the intervening shallower banks. In addition, a grounding zone for the ice sheet has been mapped at the shelf edge 10-20 km off the northwest coast of Svalbard, suggesting that ice did not reach the adjacent Yermak Plateau during the Late Weichselian.     

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.     

东南极冰盖/冰架变化监测与预测技术研究

南北极是研究全球变化的关键区域。＂十一五＂期间我国在南极地区系统开展了东南极冰盖/冰架变化监测与预测技术研究,对认识全球气候变化具有重要作用。通过项目实施,建立和发展了一批现场观测体系,发展了冰盖观测新技术并集成应用于东南极冰盖的综合观测,拓展了对冰穹A地区的新认识和新发现;在冰穹A边缘地区钻取的一支浅冰芯恢复了过去约780年的气候记录,记录了东南极地区存在小冰期的明显证据;发展了东南极冰盖积累和等时年层流动模型,研究在冰盖浅层、中层和深部的变化特征,反演了冰穹A地区的古积累率分布。本文概要介绍该项目基本情况。     

The last British Ice Sheet: A review of the evidence utilised in the compilation of the Glacial Map of Britain

This paper reviews the evidence presently available (as at December 2003) for the compilation of the Glacial Map of Britain (see [Clark C.D., Evans D.J.A., Khatwa A., Bradwell T., Jordan C.J., Marsh S.H., Mitchell W.A., Bateman, M.D. , 2004. Map and GIS database of glacial landforms and features related to the last British Ice Sheet. Boreas 33, 359–375] and http://www.shef.ac.uk/geography/staff/clark_chris/britice.html) in an effort to stimulate further research on the last British Ice Sheet and promote a reconstruction of ice sheet behaviour based on glacial geology and geomorphology. The wide range of evidence that has been scrutinized for inclusion on the glacial map is assessed with respect to the variability of its quality and quantity and the existing controversies in ice sheet reconstructions. Landforms interpreted as being of unequivocal ice-marginal origin (moraines, ice-contact glacifluvial landforms and lateral meltwater channels) and till sheet margins are used in conjunction with available chronological control to locate former glacier and ice-sheet margins throughout the last glacial cycle. Subglacial landforms (drumlins, flutings and eskers) have been used to demarcate former flow patterns within the ice sheet. The compilation of evidence in a regional map is crucial to any future reconstructions of palaeo-ice sheet dynamics and will provide a clearer understanding of ice sheet configuration, ice divide migration and ice thickness and coverage for the British Ice Sheet as it evolved through the last glacial cycle.     

Capabilities and limitations of numerical ice sheet models: a discussion for Earth-scientists and modelers

The simulation of dynamically coupled ice sheet, ice stream, and ice shelf-systems poses a challenge to most numerical ice sheet models. Here we review present ice sheet model limitations targeting a broader audience within Earth Sciences, also those with no specific background in numerical modeling, in order to facilitate cross-disciplinary communication between especially paleoglaciologists, marine and terrestrial geologists, and numerical modelers. The ‘zero order’ (Shallow Ice Approximation, SIA)-, ‘higher order’-, and ‘full Stokes’ ice sheet models are described conceptually and complemented by an outline of their derivations. We demonstrate that higher order models are required to simulate coupled ice sheet-ice shelf and ice sheet-ice stream systems, in particular if the results are aimed to complement spatial ice flow reconstructions based on higher resolution geological and geophysical data. The zero order SIA model limitations in capturing ice stream behavior are here illustrated by conceptual simulations of a glaciation on Svalbard. The limitations are obvious from the equations comprising a zero order model. However, under certain circumstances, simulation results may falsely give the impression that ice streams indeed are simulated with a zero order SIA model.