Rapid discharge connects Antarctic subglacial lakes

The existence of many subglacial lakes provides clear evidence for the widespread presence of water beneath the East Antarctic ice sheet, but the hydrology beneath this ice mass is poorly understood. Such knowledge is critical to understanding ice flow, basal water transfer to the ice margin, glacial landform development and subglacial lake habitats. Here we present ice-sheet surface elevation changes in central East Antarctica that we interpret to represent rapid discharge from a subglacial lake. Our observations indicate that during a period of 16 months, 1.8 km3 of water was transferred over 290 km to at least two other subglacial lakes. While viscous deformation of the ice roof above may moderate discharge, the intrinsic instability of such a system suggests that discharge events are a common mode of basal drainage. If large lakes, such as Lake Vostok or Lake Concordia, are pressurizing, it is possible that substantial discharges could reach the coast. Our observations conflict with expectations that subglacial lakes have long residence times and slow circulations, and we suggest that entire subglacial drainage basins may be flushed periodically. The rapid transfer of water between lakes would result in large-scale solute and microbe relocation, and drainage system contamination from in situ exploration is, therefore, a distinct risk.     

Inland thinning of West Antarctic Ice Sheet steered along subglacial rifts

Current ice loss from the West Antarctic Ice Sheet (WAIS) accounts for about ten per cent of observed global sea-level rise. Losses are dominated by dynamic thinning, in which forcings by oceanic or atmospheric perturbations to the ice margin lead to an accelerated thinning of ice along the coastline. Although central to improving projections of future ice-sheet contributions to global sea-level rise, the incorporation of dynamic thinning into models has been restricted by lack of knowledge of basal topography and subglacial geology so that the rate and ultimate extent of potential WAIS retreat remains difficult to quantify. Here we report the discovery of a subglacial basin under Ferrigno Ice Stream up to 1.5?kilometres deep that connects the ice-sheet interior to the Bellingshausen Sea margin, and whose existence profoundly affects ice loss. We use a suite of ice-penetrating radar, magnetic and gravity measurements to propose a rift origin for the basin in association with the wider development of the West Antarctic rift system. The Ferrigno rift, overdeepened by glacial erosion, is a conduit which fed a major palaeo-ice stream on the adjacent continental shelf during glacial maxima. The palaeo-ice stream, in turn, eroded the 'Belgica' trough, which today routes warm open-ocean water back to the ice front to reinforce dynamic thinning. We show that dynamic thinning from both the Bellingshausen and Amundsen Sea region is being steered back to the ice-sheet interior along rift basins. We conclude that rift basins that cut across the WAIS margin can rapidly transmit coastally perturbed change inland, thereby promoting ice-sheet instability.     

Origin and fate of Lake Vostok water frozen to the base of the East Antarctic ice sheet

The subglacial Lake Vostok may be a unique reservoir of genetic material and it may contain organisms with distinct adaptations, but it has yet to be explored directly. The lake and the overlying ice sheet are closely linked, as the ice-sheet thickness drives the lake circulation, while melting and freezing at the ice-sheet base will control the flux of water, biota and sediment through the lake. Here we present a reconstruction of the ice flow trajectories for the Vostok core site, using ice-penetrating radar data and Global Positioning System (GPS) measurements of surface ice velocity. We find that the ice sheet has a significant along-lake flow component, persistent since the Last Glacial Maximum. The rates at which ice is frozen (accreted) to the base of the ice sheet are greatest at the shorelines, and the accreted ice layer is subsequently transported out of the lake. Using these new flow field and velocity measurements, we estimate the time for ice to traverse Lake Vostok to be 16,000-20,000 years. We infer that most Vostok ice analysed to date was accreted to the ice sheet close to the western shoreline, and is therefore not representative of open lake conditions. From the amount of accreted lake water we estimate to be exported along the southern shoreline, the lake water residence time is about 13,300 years.     

Water exchange between the subglacial Lake Vostok and the overlying ice sheet

It has now been known for several years that a 200-km-long lake, called Lake Vostok, lies beneath the ice sheet on which sits Vostok Station in Antarctica. The conditions at the base of the ice sheet above this subglacial lake can provide information about the environment within the lake, including the likelihood that it supports life. Here we present an analysis of the ice-sheet structure from airborne 60-MHz radar studies, which indicates that distinct zones of basal ice loss and accretion occur at the ice-water interface. Subglacial melting and net ice loss occur in the north of the lake and across its 200-km-long western margin, whereas about 150 m of ice is gained by subglacial freezing in the south. This indicates that significant quantities of water are exchanged between the base of the ice sheet and the lake waters, which will enrich the lake with gas hydrates, cause sediment deposition and encourage circulation of the lake water.     

DINSAR measurement of glacier motion in Antarctic Grove Mountain

Grove Mountain is an important nunatak region on East Antarctic Glacier that blocks the ice flow to- ward Lambert Glacier. The existence of nunataks and subglacial mountains leads to complex ice flow patterns, which are difficult to be measured by conventional ground-based methods. In this study, several JERS-1 and ERS-1/2 SAR images covering this area are used for 3-pass and 4-pass differential interferometric processing. The ice flow field of Grove Mountain and the eastern zone are derived and validated with related knowledge. The research shows that DINSAR is an effective method for meas- uring complex ice flow in Antarctic inland glacier. L-band DINSAR is more suitable for mid or fast ice flow than C-band over this region.     

Physical, chemical and biological processes in Lake Vostok and other Antarctic subglacial lakes.

Over 70 lakes have now been identified beneath the Antarctic ice sheet. Although water from none of the lakes has been sampled directly, analysis of lake ice frozen (accreted) to the underside of the ice sheet above Lake Vostok, the largest of these lakes, has allowed inferences to be made on lake water chemistry and has revealed small quantities of microbes. These findings suggest that Lake Vostok is an extreme, yet viable, environment for life. All subglacial lakes are subject to high pressure (approximately 350 atmospheres), low temperatures (about -3 degrees C) and permanent darkness. Any microbes present must therefore use chemical sources to power biological processes. Importantly, dissolved oxygen is available at least at the lake surface, from equilibration with air hydrates released from melting basal glacier ice. Microbes found in Lake Vostok's accreted ice are relatively modern, but the probability of ancient lake-floor sediments leads to a possibility of a very old biota at the base of subglacial lakes.     

Melt-induced speed-up of Greenland ice sheet offset by efficient subglacial drainage

Fluctuations in surface melting are known to affect the speed of glaciers and ice sheets, but their impact on the Greenland ice sheet in a warming climate remains uncertain. Although some studies suggest that greater melting produces greater ice-sheet acceleration, others have identified a long-term decrease in Greenland's flow despite increased melting. Here we use satellite observations of ice motion recorded in a land-terminating sector of southwest Greenland to investigate the manner in which ice flow develops during years of markedly different melting. Although peak rates of ice speed-up are positively correlated with the degree of melting, mean summer flow rates are not, because glacier slowdown occurs, on average, when a critical run-off threshold of about 1.4?centimetres a day is exceeded. In contrast to the first half of summer, when flow is similar in all years, speed-up during the latter half is 62?±?16 per cent less in warmer years. Consequently, in warmer years, the period of fast ice flow is three times shorter and, overall, summer ice flow is slower. This behaviour is at odds with that expected from basal lubrication alone. Instead, it mirrors that of mountain glaciers, where melt-induced acceleration of flow ceases during years of high melting once subglacial drainage becomes efficient. A model of ice-sheet flow that captures switching between cavity and channel drainage modes is consistent with the run-off threshold, fast-flow periods, and later-summer speeds we have observed. Simulations of the Greenland ice-sheet flow under climate warming scenarios should account for the dynamic evolution of subglacial drainage; a simple model of basal lubrication alone misses key aspects of the ice sheet's response to climate warming.     

Ice radar investigation at Dome A, East Antarctica: Ice thickness and subglacial topography

Dome A, located in the central East Antarctic ice sheet (EAIS), is the highest summit of the Antarctic ice sheet. From ice-sheet evolution modeling results, Dome A is likely to preserve over one million years of the Earth’s paleo-climatic and -environmental records, and considered an ideal deep ice core drilling site. Ice thickness and subglacial topography are critical factors for ice-sheet models to determine the timescale and location of a deep ice core. During the 21st and 24th Chinese National Antarctic Research Expedition (CHINARE 21, 2004/05; CHINARE 24, 2007/08), ground-based ice radar systems were used to a three-dimensional investigation in the central 30 km×30 km region at Dome A. The successfully obtained high resolution and accuracy data of ice thickness and subglacial topography were then interpolated into the ice thickness distribution and subglacial topography digital elevation model (DEM) with a regular grid resolution of 140.5 m×140.5 m. The results of the ice radar investigation indicate that the average ice thickness in the Dome A central 30 km×30 km region is 2233 m, with a minimal ice thickness of 1618 m and a maximal ice thickness of 3139 m at Kunlun Station. The subglacial topography is relatively sharp, with an elevation range of 949–2445 m. The typical, clear mountain glaciation morphology is likely to reflect the early evolution of the Antarctic ice sheet. Based on the ice thickness distribution and subglacial topography characteristics, the location of Kunlun Station was suggested to carry out the first high-resolution, long time-scale deep ice core drilling. However, the internal structure and basal environments at Kunlun Station still need further research to determine.     

Ice-sheet acceleration driven by melt supply variability

Increased ice velocities in Greenland are contributing significantly to eustatic sea level rise. Faster ice flow has been associated with ice-ocean interactions in water-terminating outlet glaciers and with increased surface meltwater supply to the ice-sheet bed inland. Observed correlations between surface melt and ice acceleration have raised the possibility of a positive feedback in which surface melting and accelerated dynamic thinning reinforce one another, suggesting that overall warming could lead to accelerated mass loss. Here I show that it is not simply mean surface melt but an increase in water input variability that drives faster ice flow. Glacier sliding responds to melt indirectly through changes in basal water pressure, with observations showing that water under glaciers drains through channels at low pressure or through interconnected cavities at high pressure. Using a model that captures the dynamic switching between channel and cavity drainage modes, I show that channelization and glacier deceleration rather than acceleration occur above a critical rate of water flow. Higher rates of steady water supply can therefore suppress rather than enhance dynamic thinning, indicating that the melt/dynamic thinning feedback is not universally operational. Short-term increases in water input are, however, accommodated by the drainage system through temporary spikes in water pressure. It is these spikes that lead to ice acceleration, which is therefore driven by strong diurnal melt cycles and an increase in rain and surface lake drainage events rather than an increase in mean melt supply.     

Switch of flow direction in an Antarctic ice stream

Fast-flowing ice streams transport ice from the interior of West Antarctica to the ocean, and fluctuations in their activity control the mass balance of the ice sheet. The mass balance of the Ross Sea sector of the West Antarctic ice sheet is now positive--that is, it is growing--mainly because one of the ice streams (ice stream C) slowed down about 150 years ago. Here we present evidence from both surface measurements and remote sensing that demonstrates the highly dynamic nature of the Ross drainage system. We show that the flow in an area that once discharged into ice stream C has changed direction, now draining into the Whillans ice stream (formerly ice stream B). This switch in flow direction is a result of continuing thinning of the Whillans ice stream and recent thickening of ice stream C. Further abrupt reorganization of the activity and configuration of the ice streams over short timescales is to be expected in the future as the surface topography of the ice sheet responds to the combined effects of internal dynamics and long-term climate change. We suggest that caution is needed when using observations of short-term mass changes to draw conclusions about the large-scale mass balance of the ice sheet.     

North American ice-sheet dynamics and the onset of 100,000-year glacial cycles

The onset of major glaciations in the Northern Hemisphere about 2.7 million years ago was most probably induced by climate cooling during the late Pliocene epoch. These glaciations, during which the Northern Hemisphere ice sheets successively expanded and retreated, are superimposed on this long-term climate trend, and have been linked to variations in the Earth's orbital parameters. One intriguing problem associated with orbitally driven glacial cycles is the transition from 41,000-year to 100,000-year climatic cycles that occurred without an apparent change in insolation forcing. Several hypotheses have been proposed to explain the transition, both including and excluding ice-sheet dynamics. Difficulties in finding a conclusive answer to this palaeoclimatic problem are related to the lack of sufficiently long records of ice-sheet volume or sea level. Here we use a comprehensive ice-sheet model and a simple ocean-temperature model to extract three-million-year mutually consistent records of surface air temperature, ice volume and sea level from marine benthic oxygen isotopes. Although these records and their relative phasings are subject to considerable uncertainty owing to limited availability of palaeoclimate constraints, the results suggest that the gradual emergence of the 100,000-year cycles can be attributed to the increased ability of the merged North American ice sheets to survive insolation maxima and reach continental-scale size. The oversized, wet-based ice sheet probably responded to the subsequent insolation maximum by rapid thinning through increased basal-sliding, thereby initiating a glacial termination. Based on our assessment of the temporal changes in air temperature and ice volume during individual glacials, we demonstrate the importance of ice dynamics and ice-climate interactions in establishing the 100,000-year glacial cycles, with enhanced North American ice-sheet growth and the subsequent merging of the ice sheets being key elements.     

新疆冰湖溃决灾害风险评价与预警研究

受全球气候变暖及人类活动的影响,冰湖溃决灾害已成为当前新疆地区经济可持续发展中的突出问题之一。利用1977(MSS)～2015年Landsat影像、2015年的Google Earth高精度影像和DEM等数据解译获取新疆冰湖的时空分布信息,并选取导致冰湖溃决的主要影响因子,运用MATLAB和SPSS等软件,利用西北地区已有冰湖数据作为训练样本拟合评价方程,建立了新疆冰湖溃决风险评价模型和基于BP神经网络算法的预警模型。结果表明:(1)1977年以来新疆三大山系中冰湖的数量和面积不断增加,冰湖数量上增加了62.0%,冰湖面积增长率为0.98 km~2·a~(-1),但增长率呈现逐年减小的趋势;(2)对分布在新疆范围内的17个典型冰湖溃决风险评价的结果显示:有5个冰湖为高危风险、5个冰湖为中等风险、7个冰湖为低等风险,新疆大约29%的冰湖需要进行实时监测。     

Antarctic ice-sheet loss driven by basal melting of ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Atmospheric and oceanic forcing have the potential to reduce the thickness and extent of floating ice shelves, potentially limiting their ability to buttress the flow of grounded tributary glaciers. Indeed, recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. But the extent and magnitude of ice-shelf thickness change, the underlying causes of such change, and its link to glacier flow rate are so poorly understood that its future impact on the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary control of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet leading to accelerated glacier flow. The highest thinning rates occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen seas, and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic ice-sheet mass balance, and hence global sea level, on annual to decadal timescales.     

Preliminary results of ice radar investigation along the traverse between Zhongshan and Dome A in East Antarctic ice sheet: Ice thickness and subglacial topography

During the 24th Chinese National Antarctic Research Expedition (CHINARE 24, 2007/08), a ground-based ice radar was used to survey ice thickness and subglacial topography along the 1170 km traverse between Zhongshan and Dome A in East Antarctic ice sheet (EAIS). Ice-bedrock interface was detected along 82% of the traverse and data was collected at a horizontal resolution of <5.6 m. The data was processed to produce curves of ice thickness distribution and subglacial topography along the traverse. The results indicate that, along the traverse, the average ice thickness is 2037 m, smaller than the average ice thickness in EAIS; the thickest ice is at 730 km, and the thinnest ice (891 m) is at the edge of the ice sheet, but the slightly larger ice thickness (1078 m) in inland appears at 1020 km; the average subglacial topography elevation is 728 m, greatly larger than the average value in EAIS, and the largest elevation reaches up to 2650 m at 1034 km. The lowest terrain is located at 765 km. In further inland of 900–1170 km, the subglacial topography is relatively high due to the existence of the Gamburtsev Subglacial Mountains in the region. Generally, the influence of subglacial topography on ice surface is not significant, except at 900 km where great rise of subglacial topography causes evident uplift of ice surface. Where ice-bedrock interface was detected, the frequent and strong change of ice thickness and subglacial topography in small-scale means large bedrock roughness along the traverse, and is considered as the result of the integrated influence of ice flow, basal environments and geology. The segment where bedrock was not detected has very large ice thickness. The strong ice flow there also makes internal structure more complicated and induces serious attenuation of radar signals.     

冰碛湖库容数学模型的计算方法

全球气候的变暖致使冰湖溃决灾害频繁发生,严重威胁高原地区人民生命财产安全,随之对于冰湖溃决的研究也愈发深入,冰湖库容量作为冰湖研究基础参数,对冰湖危险性评价、冰湖溃决洪峰流量以及模拟洪水演进具有不可代替的作用。计算冰湖库容的方法主要有实地测量和冰碛湖库容计算公式,由于冰碛湖多位于高山极高山地区,对每个冰湖进行实地测量难以实现,传统冰湖库容的计算方法是基于大量冰湖参数与冰湖库容的统计关系,建立冰碛湖库容经验公式,虽然其拟合优度相对较高,但其相关性只能衡量冰碛湖面积和冰碛湖库容在统计规律上的一致性,对于精确计算冰碛湖库容仍有难以完成。本文通过文献资料整理、结合实地考察对湖盆形态模拟分析,将典型冰碛湖湖盆分为三段,通过模拟建立数学模型,以冰湖面积为约束条件,抽象成数学中的曲面立体模型,分别计算其形态体积,由此推算出计算冰湖体积的表达式V=0.0717w²l (w-冰碛湖湖宽;l-冰碛湖整体长度),并利用该模型验证了公式的有效精确性,为我国喜马拉雅山地区冰碛湖库容计算提供理论参考。     

Constraints on hydrothermal processes and water exchange in Lake Vostok from helium isotopes

Lake Vostok, the largest subglacial lake in Antarctica, is covered by the East Antarctic ice sheet, which varies in thickness between 3,750 and 4,100 m (ref. 1). At a depth of 3,539 m in the drill hole at Vostok station, sharp changes in stable isotopes and the gas content of the ice delineate the boundary between glacier ice and ice accreted through re-freezing of lake water. Unlike most gases, helium can be incorporated into the crystal structure of ice during freezing, making helium isotopes in the accreted ice a valuable source of information on lake environment. Here we present helium isotope measurements from the deep section of the Vostok ice core that encompasses the boundary between the glacier ice and accreted ice, showing that the accreted ice is enriched by a helium source with a radiogenic isotope signature typical of an old continental province. This result rules out any significant hydrothermal energy input into the lake from high-enthalpy mantle processes, which would be expected to produce a much higher 3He/4He ratio. Based on the average helium flux for continental areas, the helium budget of the lake leads to a renewal time of the lake of the order of 5,000 years.     

抚仙湖集水域径流过程的观测与模拟

抚仙湖集水域的径流过程包括河道径流、坡面散流和地下径流,各径流间发生水量交换,且同时向湖泊输送污染物。针对这一特点,提出以整体水文过程思路建立径流及物质输移模型,以完整统一的模型模拟整个集水域的径流过程及污染物入湖路径。模拟方法考虑各径流之间的水量交换,能定量区分各径流量及各径流分别向湖泊输送的物质数量,可望深层次地揭示湖泊集水域的径流系统及物质输移路径。设计了相应的野外观测方案,对各径流及降雨量作同步连续观测。观测数据显示,该集水域降雨强度分布不均匀;河道径流量与降雨相关性强,表现出暴涨暴落的特点,反映了山区河道径流特征;北岸农田区地下水埋深较浅,表层土体松散,接受良好的降雨补给。这些成果为观测方案的进一步细化和模型建设提供了可靠的依据。     

1990—2020年那曲地区冰湖变化研究

基于Landsat遥感影像,采用目视解译的方法提取了青藏高原内部那曲地区冰湖和冰川1990、2000、2010及2020年4期边界数据,并分析近30年来该地区冰湖变化的特征与原因．结果表明:那曲地区现有冰湖255个,总面积（27.829±4.62） km2,冰湖主要集中在东南部,其次是西南部;1990—2020年,研究区新增冰湖72个,面积增长6.14 km2;冰湖扩张趋势明显,具体表现为低海拔地区（≤4 700 m）原有冰湖快速扩张,较高海拔地区（>5 100~5 700 m）新增冰湖集中出现;气温与降水是冰湖变化的关键因素,由于降水与气温分布及变化存在时空差异,冰湖变化呈现“南快北慢,逐期加快”的特征;冰湖与冰川的位置关系也会影响冰湖变化,离冰川位置越近,对冰湖发育越有利,同时与冰湖接触的冰川退缩速率相较于其余冰川有更大的退缩速率,但冰川与冰湖面积变化速率并无明显相关性．      

A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ～34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5?km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved understanding of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1?km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two iceflow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.     

Vertical extension of the subglacial drainage system into basal crevasses

Water plays a first-order role in basal sliding of glaciers and ice sheets and is often a key constituent of accelerated glacier motion. Subglacial water is known to occupy systems of cavities and conduits at the interface between ice and the underlying bed surface, depending upon the history of water input and the characteristics of the substrate. Full understanding of the extent and configuration of basal water is lacking, however, because direct observation is difficult. This limits our ability to simulate ice dynamics and the subsequent impacts on sea-level rise realistically. Here we show that the subglacial hydrological system can have a large volume of water occupying basal crevasses that extend upward from the bed into the overlying ice. Radar and seismic imaging combined with in situ borehole measurements collected on Bench Glacier, Alaska, reveal numerous water-filled basal crevasses with highly transmissive connections to the bed. Some crevasses extend many tens of metres above the bed and together they hold a volume of water equivalent to at least a decimetre layer covering the bed. Our results demonstrate that the basal hydrologic system can extend high into the overlying ice mass, where basal crevasses increase water-storage capacity and could potentially modulate basal water pressure. Because basal crevasses can form under commonly observed glaciological conditions, our findings have implications for interpreting and modelling subglacial hydrologic processes and related sliding accelerations of glaciers and ice sheets.