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.     

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.     

Large subglacial lakes in East Antarctica at the onset of fast-flowing ice streams

Water plays a crucial role in ice-sheet stability and the onset of ice streams. Subglacial lake water moves between lakes and rapidly drains, causing catastrophic floods. The exact mechanisms by which subglacial lakes influence ice-sheet dynamics are unknown, however, and large subglacial lakes have not been closely associated with rapidly flowing ice streams. Here we use satellite imagery and ice-surface elevations to identify a region of subglacial lakes, similar in total area to Lake Vostok, at the onset region of the Recovery Glacier ice stream in East Antarctica and predicted by ice-sheet models. We define four lakes through extensive, flat, featureless regions of ice surface bounded by upstream troughs and downstream ridges. Using ice velocities determined using interferometric synthetic aperture radar (InSAR), we find the onset of rapid flow (moving at 20 to 30 m yr(-1)) of the tributaries to the Recovery Glacier ice stream in a 280-km-wide segment at the downslope margins of these four subglacial lakes. We conclude that the subglacial lakes initiate and maintain rapid ice flow through either active modification of the basal thermal regime of the ice sheet by lake accretion or through scouring bedrock channels in periodic drainage events. We suggest that the role of subglacial lakes needs to be considered in ice-sheet mass balance assessments.     

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-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.     

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.     

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.     

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.     

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.     

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.     

Similar meltwater contributions to glacial sea level changes from Antarctic and northern ice sheets

The period between 75,000 and 20,000 years ago was characterized by high variability in climate and sea level. Southern Ocean records of ice-rafted debris suggest a significant contribution to the sea level changes from melt water of Antarctic origin, in addition to likely contributions from northern ice sheets, but the relative volumes of melt water from northern and southern sources have yet to be established. Here we simulate the first-order impact of a range of relative meltwater releases from the two polar regions on the distribution of marine oxygen isotopes, using an intermediate complexity model. By comparing our simulations with oxygen isotope data from sediment cores, we infer that the contributions from Antarctica and the northern ice sheets to the documented sea level rises between 65,000 and 35,000 years ago were approximately equal, each accounting for a rise of about 15 m. The reductions in Antarctic ice volume implied by our analysis are comparable to that inferred previously for the Antarctic contribution to meltwater pulse 1A (refs 16, 17), which occurred about 14,200 years ago, during the last deglaciation.     

Fractures as the main pathways of water flow in temperate glaciers

Understanding the flow of water through the body of a glacier is important, because the spatial distribution of water and the rate of infiltration to the glacier bottom is one control on water storage and pressure, glacier sliding and surging, and the release of glacial outburst floods. According to the prevailing hypothesis, this water flow takes place in a network of tubular conduits. Here we analyse video images from 48 boreholes drilled into the small Swedish glacier Storglaci?ren, showing that the glacier's hydrological system is instead dominated by fractures that convey water at slow speeds. We detected hydraulically connected fractures at all depths, including near the glacier bottom. Our observations indicate that fractures provide the main pathways for surface water to reach deep within the glacier, whereas tubular conduits probably form only in special circumstances. A network of hydraulically linked fractures offers a simple explanation for the origin and evolution of the englacial water flow system and its seasonal regeneration. Such a fracture network also explains radar observations that reveal a complex pattern of echoes rather than a system of conduits. Our findings may be important in understanding the catastrophic collapse of ice shelves and rapid hydraulic connection between the surface and bed of an ice sheet.     

Active formation of 'chaos terrain' over shallow subsurface water on Europa

Europa, the innermost icy satellite of Jupiter, has a tortured young surface and sustains a liquid water ocean below an ice shell of highly debated thickness. Quasi-circular areas of ice disruption called chaos terrains are unique to Europa, and both their formation and the ice-shell thickness depend on Europa's thermal state. No model so far has been able to explain why features such as Conamara Chaos stand above surrounding terrain and contain matrix domes. Melt-through of a thin (few-kilometre) shell is thermodynamically improbable and cannot raise the ice. The buoyancy of material rising as either plumes of warm, pure ice called diapirs or convective cells in a thick (>10 kilometres) shell is insufficient to produce the observed chaos heights, and no single plume can create matrix domes. Here we report an analysis of archival data from Europa, guided by processes observed within Earth's subglacial volcanoes and ice shelves. The data suggest that chaos terrains form above liquid water lenses perched within the ice shell as shallow as 3?kilometres. Our results suggest that ice-water interactions and freeze-out give rise to the diverse morphologies and topography of chaos terrains. The sunken topography of Thera Macula indicates that Europa is actively resurfacing over a lens comparable in volume to the Great Lakes in North America.     

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.     

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.     

Timing of the Last Glacial Maximum from observed sea-level minima

During the Last Glacial Maximum, ice sheets covered large areas in northern latitudes and global temperatures were significantly lower than today. But few direct estimates exist of the volume of the ice sheets, or the timing and rates of change during their advance and retreat. Here we analyse four distinct sediment facies in the shallow, tectonically stable Bonaparte Gulf, Australia--each of which is characteristic of a distinct range in sea level--to estimate the maximum volume of land-based ice during the last glaciation and the timing of the initial melting phase. We use faunal assemblages and preservation status of the sediments to distinguish open marine, shallow marine, marginal marine and brackish conditions, and estimate the timing and the mass of the ice sheets using radiocarbon dating and glacio-hydro-isostatic modelling. Our results indicate that from at least 22,000 to 19,000 (calendar) years before present, land-based ice volume was at its maximum, exceeding today's grounded ice sheets by 52.5 x 10(6) km. A rapid decrease in ice volume by about 10% within a few hundred years terminated the Last Glacial Maximum at 19,000 +/- 250 years.     

黄河流域水量调度的自校正控制模型

为克服单纯人工调水的局限,满足现代黄河水调工作的需要,建立了黄河流域水量调度的自校正控制模型。以防凌、防断流以及"公平分水"为模型控制目标,以黄河干流骨干水库为调节器,及时跟踪、掌握水调系统控制流量、引水流量、水库泄量、水位等状态参数,调节水库泄流、调整引水计划。运用该模型编制了黄河流域2001至2002水文年的水调方案并与当年实际情况对比。结果显示与实际水调情况下的控制断面流量平均偏差1.05%,小浪底水库下泄流量平均偏差2.77%,各省区耗水量平均偏差4.40%,省区耗水过程偏差14.37%。该模型目前已经在黄河水量调度中心投入运行。     

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.     

Recent mass balance of polar ice sheets inferred from patterns of global sea-level change

Global sea level is an indicator of climate change, as it is sensitive to both thermal expansion of the oceans and a reduction of land-based glaciers. Global sea-level rise has been estimated by correcting observations from tide gauges for glacial isostatic adjustment--the continuing sea-level response due to melting of Late Pleistocene ice--and by computing the global mean of these residual trends. In such analyses, spatial patterns of sea-level rise are assumed to be signals that will average out over geographically distributed tide-gauge data. But a long history of modelling studies has demonstrated that non-uniform--that is, non-eustatic--sea-level redistributions can be produced by variations in the volume of the polar ice sheets. Here we present numerical predictions of gravitationally consistent patterns of sea-level change following variations in either the Antarctic or Greenland ice sheets or the melting of a suite of small mountain glaciers. These predictions are characterized by geometrically distinct patterns that reconcile spatial variations in previously published sea-level records. Under the--albeit coarse--assumption of a globally uniform thermal expansion of the oceans, our approach suggests melting of the Greenland ice complex over the last century equivalent to -0.6 mm yr(-1) of sea-level rise.     

DUT-1非冻结合成模型冰物模技术及应用

基于渤海海水的物理和力学特性、海冰物理模型试验的相似理论以及国际上已有模型冰的经验,较成功地研制出一种用于模拟我国渤海海冰的非冻结可破碎功能材料－－DUT－1模型冰。该模型冰采用聚丙烯粉、水泥、塑料微珠等材料,加入一定量的水、由特制的搅拌机搅拌均匀倒入模具后成型制成。利用DUT－1模型冰材料,在实验室常温波流水槽中进行了冰排在波浪作用下的断裂、冰排对直桩及锥体结构物结构作用的物理模型试验。通过对试验现象的观察,对实验结果的综合分析、对相关理论及经验公式计算结果的比较,证明DUT－1非冻结合成模型冰材料物模实验技术有较好可靠性。