(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

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 under¬standing 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 ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

Surface velocity and ice discharge of the ice cap on King George Island, Antarctica

Glaciers on King George Island, Antarctica, have shown retreat and surface lowering in recent decades, concurrent with increasing air temperatures. A large portion of the glacier perimeter is ocean-terminating, suggesting possible large mass losses due to calving and submarine melting. Here we estimate the ice discharge into the ocean for the King George Island ice cap. L-band synthetic aperture radar images covering the time-span January 2008 to January 2011 over King George Island are processed using an intensity-tracking algorithm to obtain surface velocity measurements. Pixel offsets from 40 pairs of radar images are analysed and inverted to estimate a weighted average surface velocity field. Ice thicknesses are derived from simple principles of ice flow mechanics using the computed surface velocity fields and in situ thickness data. The maximum ice surface speeds reach mayor que 225 m/yr, and the total ice discharge for the analysed flux gates of King George Island is estimated to be 0.720+/-0.428 Gt/yr, corresponding to a specific mass loss of 0.64+/-0.38 m w.e./yr over the area of the entire ice cap (1127 km2).

Sedimentological and geophysical investigation of sediment cores from the Amundsen Sea

Ice loss from the marine-based, potentially unstable West Antarctic Ice Sheet (WAIS) contributes to current sea-level rise and may raise sea level by up to 3.3 to 5 meters in the future. Over the past few decades, glaciers draining the WAIS into the Amundsen Sea Embayment (ASE) have shown accelerated ice flow, rapid thinning and grounding-line retreat. However, the long-term context of this ice-sheet retreat is poorly constrained, limiting our ability to accurately predict future WAIS behaviour. Here we present a new chronology for WAIS retreat from the inner continental shelf of the eastern ASE based on radiocarbon dates from three marine sediment cores. The ages document a retreat of the grounding line to within ~93 km of its modern position before 11.7±0.7 kyr BP (thousand years before present). This early deglaciation is consistent with ages for grounding-line retreat from the western ASE. Our new data demonstrate that, other than in the Ross Sea, WAIS retreat in the ASE has not continued progressively since the Last Glacial Maximum. Furthermore, our results suggest that the grounding-line position in the ASE was predominantly stable throughout the Holocene, and that any episodes of fast retreat similar to that observed today must have been short-lived. Alternatively, today's rapid retreat was unprecedented during the Holocene. Therefore, the current ice loss must originate in recent changes in regional climate, ocean circulation or ice-sheet dynamics. Incorporation of these results into models is essential to produce robust predictions of future ice-sheet change and its contribution to sea-level rise.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Quantifying ice-shelf basal melt and ice-stream dynamics using high-resolution DEM and GPS time series

Thesis (Ph.D.)--University of Washington, 2016-06

Methane and stable isotope record of ice core EDML from Dronning Maud Land, Antarctica

Precise knowledge of the phase relationship between climate changes in the two hemispheres is a key for understanding the Earth's climate dynamics. For the last glacial period, ice core studies have revealed strong coupling of the largest millennial-scale warm events in Antarctica with the longest Dansgaard-Oeschger events in Greenland through the Atlantic meridional overturning circulation. It has been unclear, however, whether the shorter Dansgaard-Oeschger events have counterparts in the shorter and less prominent Antarctic temperature variations, and whether these events are linked by the same mechanism. Here we present a glacial climate record derived from an ice core from Dronning Maud Land, Antarctica, which represents South Atlantic climate at a resolution comparable with the Greenland ice core records. After methane synchronization with an ice core from North Greenland, the oxygen isotope record from the Dronning Maud Land ice core shows a one-to-one coupling between all Antarctic warm events and Greenland Dansgaard-Oeschger events by the bipolar seesaw. The amplitude of the Antarctic warm events is found to be linearly dependent on the duration of the concurrent stadial in the North, suggesting that they all result from a similar reduction in the meridional overturning circulation.

Relative effect of slope and equilibrium line altitude on the retreat of Himalayan glaciers

It has been observed that a majority of glaciers in the Himalayas have been retreating. In this paper, we show that there are two major factors which control the advance/retreat of the Himalayan glaciers. They are the slope of the glacier and changes in the equilibrium line altitude. While it is well known, that these factors are important, we propose a new way of combining them and use it to predict retreat. The functional form of this model has been derived from numerical simulations using an ice-flow code. The model has been successfully applied to the movement of eight Himalayan glaciers during the past 25 years. It explains why the Gangotri glacier is retreating while Zemu of nearly the same length is stationary, even if they are subject to similar environmental changes. The model has also been applied to a larger set of glaciers in the Parbati basin, for which retreat based on satellite data is available, though over a shorter time period.

Mechanics of deformable glacier beds

My focus in this thesis is to contribute to a more thorough understanding of the mechanics of ice and deformable glacier beds. Glaciers flow under their own weight through a combination of deformation within the ice column and basal slip, which involves both sliding along and deformation within the bed. Deformable beds, which are made up of unfrozen sediment, are prevalent in nature and are often the primary contributors to ice flow wherever they are found. Their granular nature imbues them with unique mechanical properties that depend on the granular structure and hydrological properties of the bed. Despite their importance for understanding glacier flow and the response of glaciers to changing climate, the mechanics of deformable glacier beds are not well understood.

Our general approach to understanding the mechanics of bed deformation and their effect on glacier flow is to acquire synoptic observations of ice surface velocities and their changes over time and to use those observations to infer the mechanical properties of the bed. We focus on areas where changes in ice flow over time are due to known environmental forcings and where the processes of interest are largely isolated from other effects. To make this approach viable, we further develop observational methods that involve the use of mapping radar systems. Chapters 2 and 5 focus largely on the development of these methods and analysis of results from ice caps in central Iceland and an ice stream in West Antarctica. In Chapter 3, we use these observations to constrain numerical ice flow models in order to study the mechanics of the bed and the ice itself. We show that the bed in an Iceland ice cap deforms plastically and we derive an original mechanistic model of ice flow over plastically deforming beds that incorporates changes in bed strength caused by meltwater flux from the surface. Expanding on this work in Chapter 4, we develop a more detailed mechanistic model for till-covered beds that helps explain the mechanisms that cause some glaciers to surge quasi-periodically. In Antarctica, we observe and analyze the mechanisms that allow ocean tidal variations to modulate ice stream flow tens of kilometers inland. We find that the ice stream margins are significantly weakened immediately upstream of the area where ice begins to float and that this weakening likely allows changes in stress over the floating ice to propagate through the ice column.

Displacement and mixing ventilation driven by opposing wind and buoyancy

The effect of an opposing wind on the stratification and flow produced by a buoyant plume rising from a heat source on the floor of a ventilated enclosure is investigated. Ventilation openings located at high level on the windward side of the enclosure and at low level on the leeward side allow a wind-driven flow from high to low level, opposite to the buoyancy-driven flow. One of two stable steady flow regimes is established depending on a dimensionless parameter F that characterizes the relative magnitudes of the wind-driven and buoyancy-driven velocities within the enclosure, and on the time history of the flow. A third, unstable steady flow solution is identified. For small opposing winds (small F) a steady, two-layer stratification and displacement ventilation is established. Exterior fluid enters through the lower leeward openings and buoyant interior fluid leaves through the upper windward openings. As the wind speed increases, the opposing wind may cause a reversal in the flow direction. In this case, cool exterior fluid enters through the high windward openings and mixes the interior fluid, which exits through the leeward openings. There are now two possibilities. If the rate of heat input by the source exceeds the rate of heat loss through the leeward openings, the temperature of the interior increases and this flow reversal is only maintained temporarily. The buoyancy force increases with time, the flow reverts to its original direction, and steady two-layer displacement ventilation is re-established and maintained. In this regime, the increase in wind speed increases the depth and temperature of the warm upper layer, and reduces the ventilation flow rate. If, on the other hand, the heat loss exceeds the heat input, the interior cools and the buoyancy-driven flow decreases. The reversed flow is maintained, the stratification is destroyed and mixing ventilation occurs. Further increases in wind speed increase the ventilation rate and decrease the interior temperature. The transitions between the two ventilation flow patterns exhibit hysteresis. The change from displacement ventilation to mixing ventilation occurs at a higher F than the transition from mixing to displacement. Further, we find that the transition from mixing to displacement ventilation occurs at a fixed value of F, whereas the transition from displacement to mixing flow is dependent on the details of the time history of the flow and the geometry of the openings, and is not determined solely by the value of F. Theoretical models that predic t the steady stratification profiles and flow rates for the displacement and mixing ventilation, and the transitions between them, are presented and compared with measurements from laboratory experiments. The transition between these ventilation patterns completely changes the internal environment, and we discuss some of the implications for the natural ventilation of buildings. © 2004 Cambridge University Press.

东南极格罗夫山地区的陨石回收与研究

Chinese National Antarctic Research Expedition (CHTNARE) has collected 4480 meteorite specimens in the Grove Mountains, East Antarctica, from 1998 to 2003. According to the location characteristics and the diversity of the classification, the paper concludes that the Grove Mountains is another important meteorite concentration area in the Antarctica. The Concentration mechanisms at the site could be related to the last glacier activity and katabatic wind. An empirical model was proposed: 1) Probably during the Last Glacial Maximum, ice flow overrided the Gale Escarpment range in the area. Formerly concentrated meteorites were carried by the new glacier and stayed in the terminal moraine when the glacier retreated. 2) Blown by strong katabatic wind, Newly exposed meteorites on the ablation zone were scattered on the blue ice at the lee side of the Gale escarpment. Some of them would be buried when they were moved further onto the firn snow zone. Many floating meteorites stopped and mustered at the fringe of the moraine. The chemical-petrographic of 31 meteorites were assigned based on electron probe microanalyses, petrography and mineralogy, including 1 martian lherzolitic shergottite, 1 eucrite, 1 extreme fine grain octahedron iron meteorite, and 28 ordinary chondrites (the chemical groups: 7 H-group, 13 L-group, 6 LL-group, 2 L/LL group; the petrographic types: 6 unequilibrated type 3 and 22 equilibrated type 4-6). GRV99028 meteorite has the komatiite-like spinifex texture consisting of acicular olivine crystals and some hornblende-family minerals in the interstitial region. Possibly it has crystallized from a supercooled, impact-generated, ultramafic melt of the host chondrite, then experienced the retrogressive metamorphism. Four typical chondrule textures were studied: porphyritic texture, radiative texture, barred texture and glass texture. The minerals are characteristically enriched in MgO content.

In vivo relevance for platelet glycoprotein Ib alpha residue Tyr276 in thrombus formation

Background: Platelet glycoprotein (GP) Ib-IX-V supports platelet adhesion on damaged vascular walls by binding to von Willebrand factor (VWF). For several decades it has been recognized that the alpha-subunit of GP (GPIb alpha) also binds thrombin but the physiological relevance, if any, of this interaction was unknown. Previous studies have shown that a sulfated tyrosine 276 (Tyr276) is essential for thrombin binding to GPIb alpha.Objectives: This study investigated the in vivo relevance of GPIb alpha residue Tyr276 in hemostasis and thrombosis.Methods: Transgenic mouse colonies expressing the normal human GPIb alpha subunit or a mutant human GPIb alpha containing a Phe substitution for Tyr276 (hTg(Y276F)) were generated. Both colonies were bred to mice devoid of murine GPIb alpha.Results: Surface-expressed GPIb alpha levels and platelet counts were similar in both colonies. hTg(Y276F) platelets were significantly impaired in binding alpha-thrombin but displayed normal binding to type I fibrillar collagen and human VWF in the presence of ristocetin. In vivo thrombus formation as a result of chemical damage (FeCl3) demonstrated that hTg(Y276F) mice have a delayed time to occlusion followed by unstable blood flow indicative of embolization. In models of laser-induced injury, thrombi developing in hTg(Y276F) animals were also less stable.Conclusions: The results demonstrate that GPIb alpha residue Tyr276 is physiologically important, supporting stable thrombus formation in vivo.

The surficial geology, sedimentology and geochemistry of the late glacial sediments and Paleozoic bedrock in the Campbellford area, Ontario, with special reference to the Dummer Complex /

The Dummer Complex extends 180 km along the Precambrian - Paleozoic contact from Tamworth to Lake Simcoe. It is composed of coarse, angular Paleozoic clasts in discontinuous, pitted, hummocky deposits. Deposits are usually separated by bare or boulder strewn bedrock, but have been found in the southern drumlinized till sheet. Dummer Complex deposits show rough alignment with ice-flow. Eskers cross-cut many of the deposits. Dummer sediment subfacies are defined on the basis of dominant coarse grain size and lithology, which relate directly to the underlying Paleozoic formation. Three subglacial tills are identified based on the degree of comminution and distance of transport; the immature facies of the Dummer Complex; the mature facies of the drumlinized till sheet and; the submature facies which is transitional. Carbonate geochemistry was used for till-bedrock correlation in various grain sizes. Of the 3 Paleozoic formations underlying the Dummer Complex, the Gull River Fm. is geochemically distinctive from the Bobcaygeon and Verulam Formations using Ca, Mg, Sr, Cu, Mn, Fe and Na. The Bobcaygeon Fm. and Verulam Fm. can be differentiated using Ca and the Sr/Ca ratio. The immature facies from 1.0 phi and finer is dominated by the non-carbonate, long distance transported component which decreases slightly downice. The submature till facies contains more long distance material than the immature facies. Sr and Mn can be used to correlate the Gull River immature till facies to the underlying bedrock the other subfacies could not be distinguished from each other or their respective source formation. This method proved to be ineffective for sediments with greater than 35% non-carbonate component, due to leaching of elements by the dissolving acid.The Dummer Complex is produced subglacially , as the compressional ice encounters the permeable Paleozoic carbonates. The increased shear strength of the ice and pore pressures in the carbonates results in the basal ice zones becoming debris ladden. Cleaner ice overrides the basal debris . laden dead ice which then acts as the glacier bed. During retreat, the Simcoe lobe stagnates as flow is cut-off by the Algonquin Highlands.