Calculating Basal Thermal Zones Beneath the Antarctic Ice Sheet

A procedure is presented for using a simple flowline model to calculate the fraction of the bed that is thawed beneath present-day ice sheets, and therefore for mapping thawed, frozen, melting and freezing basal thermal zones. The procedure is based on the proposition, easily demonstrated, that variations in surface slope along ice flowlines are due primarily to variations in bed topography and ice-bed coupling, where ice-bed coupling for sheet flow is represented by the basal thawed fraction. This procedure is then applied to the central flowlines of flow bands on the Antarctic ice sheet where accumulation rates, surface elevations and bed topography are mapped with sufficient accuracy, and where sheet flow rather than stream flow prevails. In East Antarctica, the usual condition is a low thawed fraction in subglacial highlands, but a high thawed fraction in subglacial basins and where ice converges on ice streams. This is consistent with a greater depression of the basal melting temperature and a slower rate of conducting basal heat to the surface where ice is thick, and greater basal frictional heat production where ice flow is fast, as expected for steady-state flow. This correlation is reduced or even reversed where steady-state flow has been disrupted recently, notably where ice-stream surges produced the Dibble and Dalton Iceberg Tongues, both of which are now stagnating. In West Antarctica, for ice draining into the Pine Island Bay polynya of the Amundsen Sea, the basal thawed fraction is consistent with a prolonged and ongoing surge of Pine Island Glacier and with a recently initiated surge of Thwaites Glacier. For ice draining into the Ross Ice Shelf, long ice streams extend nearly to the West Antarctic ice divide. Over the rugged bed topography near the ice divide, no correlation consistent with steady-state sheet flow exists between ice thickness and the basal thawed fraction. The bed is wholly thawed beneath ice streams, even where stream flow is slow. This is consistent with ongoing gravitational collapse of ice entering the Ross Sea embayment and with unstable flow in the ice streams.

Antarctic Temperatures Over the Past Two Centuries from Ice Cores

We present a reconstruction of Antarctic mean surface temperatures over the past two centuries based on water stable isotope records from high-resolution, precisely dated ice cores. Both instrumental and reconstructed temperatures indicate large interannual to decadal scale variability, with the dominant pattern being anti-phase anomalies between the main Antarctic continent and the Antarctic Peninsula region. Comparative analysis of the instrumental Southern Hemisphere (SH) mean temperature record and the reconstruction suggests that at longer timescales, temperatures over the Antarctic continent vary in phase with the SH mean. Our reconstruction suggests that Antarctic temperatures have increased by about 0.2 degrees C since the late nineteenth century. The variability and the long-term trends are strongly modulated by the SH Annular Mode in the atmospheric circulation.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Sliding of Ice Past an Obstacle at Engabreen, Norway

At Engabreen, Norway, an instrumented panel containing a decimetric obstacle was mounted flush With the bed surface beneath 210 m of ice. Simultaneous measurements of normal and shear stresses, ice velocity and temperature were obtained as dirty basal ice flowed past the obstacle. Our measurements were broadly consistent with ice thickness, flow conditions and bedrock topography near the site of the experiment. Ice speed 0.45 m above the bed was about 130 mm d(-1), much less than the surface velocity of 800 mm d(-1) Average normal stress on the panel was 1.0-1.6 MPa, smaller than the expected ice overburden pressure. Normal stress was larger and temperature was lower on the stoss side than on the lee side, in accord with flow dynamics and equilibrium thermodynamics. Annual differences in normal stresses were correlated with changes in sliding speed and ice-flow direction. These temporal variations may have been caused by changes in ice rheology associated with changes in sediment concentration, water content or both. Temperature and normal stress were generally correlated, except when clasts presumably collided with the panel. Temperature gradients in the obstacle indicated that regelation was negligible, consistent with the obstacle size. Melt rate was about 10% of the sliding speed. Despite high sliding speed, no significant ice/bed separation was observed in the lee of the obstacle. Frictional forces between sediment particles in the ice and the panel, estimated from Hallet's (1981) model, indicated that friction accounted for about 5% of the measured bed-parallel force. This value is uncertain, as friction theories are largely untested. Some of these findings agree with sliding theories, others do not.

Ice Core Evidence for a Second Volcanic Eruption Around 1809 in the Northern Hemisphere

A volcanic signal observed in ice cores from both polar regions six years prior to Tambora is attributed to an unknown tropical eruption in 1809. Recovery of dacitic tephra from the 1809 horizon in a Yukon ice core ( Eclipse) that is chemically distinct from andesitic 1809 tephra found in Antarctic ice cores indicates a second eruption in the Northern Hemisphere at this time. Together with the similar magnitude and timing of the 1809 volcanic signal in the Arctic and Antarctic, this could suggest a large tropical eruption produced the sulfate and Antarctic tephra and a minor Northern Hemisphere eruption produced the Eclipse tephra. Nonetheless, the possibility that there were coincidental eruptions of similar magnitude in both hemispheres, rather than a single tropical eruption, should not be discounted. Correctly attributing the source of the 1809 volcanic signal has important implications for modeling the magnitude and latitudinal distribution of volcanic radiative forcing.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.

Polar MM5 Simulations of the Winter Climate of the Laurentide Ice Sheet at the LGM

Optimized regional climate simulations are conducted using the Polar MM5, a version of the fifth-generation Pennsylvania State University-NCAR Mesoscale Model (MM5), with a 60-km horizontal resolution domain over North America during the Last Glacial Maximum (LGM, 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). The objective is to describe the LGM annual cycle at high spatial resolution with an emphasis on the winter atmospheric circulation. Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. Polar MM5 produces a substantially different atmospheric response to the LGM boundary conditions than CCM3 and other recent GCM simulations. In particular, from November to April the upper-level flow is split around a blocking anticyclone over the LIS, with a northern branch over the Canadian Arctic and a southern branch impacting southern North America. The split flow pattern is most pronounced in January and transitions into a single, consolidated jet stream that migrates northward over the LIS during summer. Sensitivity experiments indicate that the winter split flow in Polar MM5 is primarily due to mechanical forcing by LIS, although model physics and resolution also contribute to the simulated flow configuration. Polar MM5 LGM results are generally consistent with proxy climate estimates in the western United States, Alaska, and the Canadian Arctic and may help resolve some long-standing discrepancies between proxy data and previous simulations of the LGM climate.

Thermal Convection and the Origin of Ice Streams

Ice streams are a fact of ice-sheet dynamics, draining up to 90% of the ice. Thermal convection in ice below the density inversion is a speculation. An attempt is made to meld the two in such a way that the speculation becomes an explanation for the fact.

What Limits Bacterial Production in the Suboxic Region of Permanently Ice-Covered Lake Bonney, Antarctica?

Bacterial production assays (thymidine incorporation rates) were used to evaluate the activity of heterotrophic bacteria at the chemocline region in both the East (ELB) and West (WLB) Lobes of permanently ice-covered Lake Bonney, in the Taylor Valley of Antarctica. The magnitude of activity varied dramatically within the depth interval of 1 to 2 m from moderate to very low levels below the chemocline, especially in the East Lobe, where chemical distributions indicate the absence of a normally functioning nitrogen cycle. Several parameters (e.g. addition of nutrients or chelators, dilution) were manipulated in incubation experiments in order to identify factors that would enhance activity in the suboxic deep waters of the East Lobe. Activity, in terms of thymidine incorporation, was consistently detected in the deep-water communities, implying that, although the water may be 'toxic', the cells remain viable. None of the treatments resulted in consistent enhancement of thymidine incorporation rates in samples from below the chemocline. Bacterial populations above the chemocline appear to be phosphorus-limited. The nature of the limitation, toxicity or inhibition that limits bacterial activity in the suboxic waters has not been identified.

Denitrification in the Hypolimnion of Permanently Ice-Covered Lake Bonney, Antarctica

The distribution of denitrification was investigated in the hypolimnion of the east and west lobes of permanently ice-covered Lake Bonney, Taylor Valley, Antarctica. Anomalously high concentrations of dissolved inorganic nitrogen (DIN; nitrate, nitrite, ammonium and nitrous oxide) in the oxygen-depleted hypolimnion of the east lobe of the Lake implied that denitrification is or was active in the west, but not in the east lobe. While previous investigations reported no detectable denitrification in the east lobe, we measured active denitrification in samples from both the east and west lobes. In the west lobe, measured denitrification rates exhibited a maximum at the depth of the chemocline and denitrification was not detectable in either the oxic surface waters or in the deep water where nitrate was absent. In the east lobe, denitrification was detected below the chemocline, at the depths where ammonium, nitrate, nitrite and nitrous oxide are all present at anomalously high levels, Trace metal availability was manipulated in incubation experiments in order to determine whether trace metal toxicity in the east lobe could explain the difference in nitrogen cycling between the 2 lobes. There were no consistent stimulatory effects of metal chelators or nutrient addition on the rate of denitrification in either lobe, so that the mechanisms underlying the unusual N cycle of the east lobe remain unknown. We conclude that all the ingredients necessary to allow denitrification to occur are present in the east lobe. However, even though denitrification could be detected under certain conditions in incubations, denitrification is inhibited under the in situ conditions of the lake.

Crevassing and Calving of Glacial Ice

Calving of ice is a relatively new area of research in the still young field of glaciology. In the short time that calving has been studied, it has been mainly treated as an afterthought, with the predominant mode of thinking being that it will happen so to concern oneself with why is not important. Many studies dealt with observations of calving front positions over time vs. ice velocity in an attempt to quantify the calving rate as the difference between the two, while others have attempted to deduce some empirical relationship between calving rate and variables such as water depth or temperature. This study instead addresses the question of why, where, and when ice will first become crevassed, which is an obviously necessary condition for a later calving event to occur. Previous work examining the causes of calving used ideas put forth from a variety of fields, including civil engineering, materials science, and results from basic physics and mechanics. These theories are re-examined here and presented as part of a larger whole. Important results from the field of fracture mechanics are utilized frequently, and these results can be used as a predictor of ice behavior and intrinsic properties of ice, as well as properties like back stresses induced by local pinning points and resistive shears along glacial ice boundaries. A theory of fracture for a material experiencing creep is also presented with applications to ice shelves and crevasse penetration. Finally, a speculative theory regarding large scale iceberg formation is presented. It is meant mainly as an impetus to further discussion on the topic, with the hope that a model relating crevasse geometries to flow parameters can result in crevasse spacings that could produce the tabular icebergs which are so newsworthy. The primary focus of this thesis is to move away from the "after the fact" studies that are so common in calving research, and instead devote energy to determining what creates the conditions that drive the calving of ice in the first place.