Subglacial water plays an important role in the regulation of an ice sheet's mass balance. It may be the dominant control on the velocities of ice streams and outlet glaciers on scales of months to millennia. Recent satellite observations of ice surface elevation changes have given researchers new insights into how subglacial water is stored and transported. Localized uplift and settling of the ice surface implies that lakes exist beneath the ice sheet that are being filled and drained on relatively short time scales. %At the base of an ice sheet water can be transported through a variety of drainage networks or stored in subglacial lakes. Here, a numerical investigation of the mechanisms of transport and storage of subglacial water and the associated time scales is presented. Experiments are carried out using a finite element model of coupled ice and water flow. The first experiment seeks to understand the relationship between the depth of a basal depression and the area over which the feature affects basal water flow. It is found that as the perturbation to a topographic depression's depth is increased, water is rerouted in response to the perturbation. Additionally it is found that the relationship between perturbation depth and the extent upstream to which its effects reach is nonlinear. The second experiment examines how the aspect ratio of bed features (prolate, oblate, or equidimensional) influences basal water flow. It is found that the systems that develop and their interactions are mediated by both the topography and the feedbacks taken into account by the coupling of the systems in the model. Features oriented parallel to ice and water flow are associated with distributed fan systems that develop branches which migrate laterally across the domain and interact with one another on monthly and yearly timescales. Laterally oriented features develop laterally extensive ponds. As the ratio of longitudinal to lateral dimension of the topography is increased, a combination of these two water distributions is seen.

The lateral "corners" where Kamb and Whillans Ice Streams (KIS and WIS) discharge into the Ross Ice Shelf share common geometries and ice mechanical settings. At both corners of the now-stagnant KIS outlet, shear margins of apparently different ages confine regions with a relatively flat, smooth surface expression. These features are called the "Duckfoot" on the northern, right-lateral side and the "Goosefoot" on the other. It has been suggested, on evidence found in ice internal layers, that the flat ice terrains on KIS were afloat in the recent past, at a time when the ice stream grounding line was upstream of its present location. The overdeepening in the bed just upstream of the KIS grounding line supports this view of the past geometry. The right-lateral margin at the outlet of the currently active WIS, the location of Subglacial Lake Englehardt (SLE), appears to have many similarities with the right lateral margin of KIS, though with a less developed looking inboard margin. This paper presents a mechanical analysis using surface and bed topography and velocity datasets comparing the Duckfoot flat ice terrain with the terrain around Subglacial Lake Englehardt. At both locations mechanical thinning along shear margins and lows in the bed topography redirects basal water routing towards the features. Here, I consider the history of these features and their role in ice stream variability by comparison of the relict and modern features and via numerical modeling of ice shelf grounding and ungrounding in response to variations in ice flow. We propose two scenarios for the development of flat ice terrains/subglacial lakes at the outlets of ice streams. In the first, development of a lake in the hydraulic potential low along a shear margin forces a margin jump as shearing develops along the inboard shore of the margin lake. This thesis presents evidence for an inboard (relative to the main outboard shear margin) zone of shear along the inboard shoreline of SLE, suggesting that subglacial lakes along shear margins are capable of facilitating shear margin jumps. In the second, grounding line advance around a relative low in the bed, creating adjacent margins along the lakeshores, forms a remnant lake. Discerning which of these scenarios is appropriate at the KIS outlet has implications for understanding the history of the ice stream grounding line. An ice flow model is used to place these local conditions in a regional context by studying the effect of internal perturbations, such as ice rise stagnation or inward margin jumps, on grounding line position. Bathymetry is important in determining ice stream flow in the ways that might not be otherwise realized in 1-D flow model studies. In the numerical modeling experiments, grounding line advance across the KIS outlet is mediated by the overdeepening in the bed and proceeds not in the direction of ice flow but transverse to flow. This finding adds complexity to both a flowline view of grounding line migration and the theory that grounding lines are unstable in the presence of inward sloping bed topography.

The cold, lightless, and high pressure aquatic environment at the base of the East Antarctic Ice Sheet is of interest to a wide range of disciplines. Stable subglacial lakes and their connecting channels remain perennially liquid three kilometers below some of the coldest places on Earth. The presence of subglacial water impacts flow of the overlying ice and provides clues to the geologic properties of the bedrock below, and may harbor unique life forms which have evolved out of contact with the atmosphere for millions of years. Periodic release of water from this system may impact ocean circulation at the margins of the ice sheet. This research uses airborne radar sounding, with its unique ability to infer properties within and at the base of the ice sheet over large spatial scales, to locate and characterize this unique environment. Subglacial lakes, the primary storage mechanism for subglacial water, have been located and classified into four categories on the basis of the radar reflection properties from the sub-ice interface: Definite lakes are brighter than their surroundings by at least two decibels (relatively bright), and are both consistently reflective (specular) and have a reflection coefficient greater than --10 decibels (absolutely bright). Dim lakes are relatively bright and specular but not absolutely bright, possibly indicating non-steady dynamics in the overlying ice. Fuzzy lakes are both relatively and absolutely bright, but not specular, and may indicate saturated sediments or high frequency spatially heterogeneous distributions of sediment and liquid water (i.e. a braided steam). Indistinct lakes are absolutely bright and specular but no brighter than their surroundings. Lakes themselves and the different classes of lakes are not arranged randomly throughout Antarctica but are clustered around ice divides, ice stream onsets and prominent bedrock troughs, with each cluster demonstrating a different characteristic lake classification distribution. In the bedrock trough of Adventure Subglacial Trench, analysis of satellite altimetry is combined with radar sounding data to calculate a mass budget and infer a flow mechanism for a two cubic kilometer discharge reported to have traveled between two lakes in the region from 1996 -1998. The volume released from the source lake exceeded the volume received by the destination lakes by one and a tenth cubic kilometers, indicating that some water must have escaped downstream from the lowest destination lake over the course of the event. Release of water from the source lake preceded arrival of the water at the destination lakes, 260 kilometers away, by about three months. Water continued draining from the destination lakes for several years after surface subsidence at the source lake had ceased. By 2003, a total of one and a half cubic km or nearly 75% of the water released by the source lake had traveled downstream from the destination lakes. Hydraulic modeling work indicates that the initial release of water from the source lake could have been accommodated by a self-enlarging semicircular channel. Subsequent evolution of the discharge and the three-month delay between release of water from the source lake and arrival of that water at the destination lakes indicates that a shallower and broader distributed water system is responsible for the transport of subglacial water in this region. Such a system would be more stable for the given icebedrock geometry and may explain the observations of intermittent flat bright bedrock reflections in radar data acquired upstream from the destination lake in 2000. For the purpose of better understanding the long-term water budget of the Dome C region, an area upstream of Adventure Trench, eleven dated isochronal internal layers within the ice penetrating radar data were tracked. An age-depth relationship, derived from the European ice core through Dome C is used to calculate strain, estimate melt, model ice temperature, and determine absolute basal reflectivity for the entire region which covers over 28,000 square kilometers. The two largest subglacial lakes within the survey, Concordia and Vincennes, are both associated with enhanced basal melting on their upstream shores at rates locally greater than two millimeters per year. Widely distributed melt rates in the major topographic valleys upstream of these lakes are generally less than one millimeter per year throughout the region with slightly higher melts in the basin draining into Vincennes Subglacial Lake. Although published estimates for geothermal flux are capable of explaining the behavior of ice and water in most of the area, an additional source of basal heat is required to explain melt anomalies and subglacial lakes along the Concordia Ridge. Lake Concordia is expected to discharge water on a similar scale and duration as that observed in Adventure Trench, with a repeat cycle of a few hundred years. / text