Ice streams are fast flowing outlet glaciers through which over 90% of the ice stored within the Antarctic Ice Sheet drains. The dynamic behaviour of ice streams is therefore crucial in controlling the mass balance of the ice sheet. Over the past few decades, Antarctica has been losing mass. Much of this mass loss has been focussed around coastal regions of the Antarctic Ice Sheet. Some of the most dramatic changes such as grounding-line retreat, acceleration and surface elevation change have been observed in Pine Island Glacier (PIG) and its neighbouring ice streams. This is of particular concern because these ice streams account for 10% of the discharge from the west Antarctic Ice Sheet and therefore have the potential to contribute significantly to global sea-level rise. One of the key challenges in accurately forecasting this future sea-level rise is improving understanding of processes occurring at the beds of ice streams. This requires detailed knowledge of the properties and dynamics of the bed. This thesis aims to address this knowledge gap by investigating the spatial and temporal characteristics of the bed of PIG using high-resolution geophysical data acquired in its trunk and tributaries and beneath the ice shelf. The thesis begins by analysing radar-derived high-resolution maps of subglacial topography. These data show a contrasting topography across the ice-bed interface. These diverse subglacial landscapes have an impact on ice flow through form drag, controlled by the size and orientation of bedrock undulations and protuberances. The next chapter provides a quantitative analysis of these landscapes using Fast Fourier analysis of subglacial roughness. This analysis investigates the roughness signature of subglacial bedforms and the how the orientation and wavelength of roughness elements determine their correlation with ice dynamic parameters. The slow-flowing inter-tributary site is found to have a distinct signature comparable to 'ribbed' patterns of modelled basal shear stress and transverse 'mega rib' bedforms. Roughness oriented parallel to ice flow with wavelengths approaching mean ice thickness are found to have the highest correlation with ice dynamic parameters. The temporal stability of PIG is analysed using repeat radar measurements. No significant change is observed over a period of 3-6 years with no evidence of rapid erosion or the evolution of subglacial bedforms as observed in previous repeat measurements of ice-stream beds elsewhere in Antarctica. This suggests that the widespread deforming till layer detected in extensive seismic reflection surveys is in steady state. Lastly, the thesis explores geomorphological evidence of twentieth-century grounding-line retreat beneath PIG Ice Shelf using high-resolution geophysical data acquired from autonomous underwater vehicle surveys. Evidence of erosion, deposition, meltwater flow and post-glacial modification is observed in fine detail. The observed distribution of sediment supported previous surveys indicating a geological transition coinciding with the ridge that acted as a former stable grounding-line location. Metre-scale resolution images of recently deglaciated ice stream beds were found to reveal bedforms that are not detectable with traditional offshore bathymetric surveys. Together these findings reveal the role of short wavelength topography as both an influence on, and product of fast ice stream flow. It also highlights the spatial diversity of subglacial environments and the need to focus future research on tying detailed observations of ice-stream beds with knowledge of basal properties over time.

Mountain centered glaciers have played a major role throughout the last three million years in the Scandinavian mountains. The climatic extremes, like the present warm interglacial or cold glacial maxima, are very short-lived compared to the periods of intermediate climate conditions, characterized by the persistence of mountain based glaciers and ice fields of regional size. These have persisted in the Scandinavian mountains for about 65% of the Quaternary. Mountain based glaciers thus had a profound impact on large-scale geomorphology, which is manifested in large-scale glacial landforms such as fjords, glacial lakes and U-shaped valleys in and close to the mountain range. Through a mapping of glacial landforms in the northern Scandinavian mountain range, in particular a striking set of lateral moraines, this thesis offers new insights into Weichselian stages predating the last glacial maximum. The aerial photograph mapping and field evidence yield evidence that these lateral moraines were overridden by glacier ice subsequent to their formation. The lateral moraines were dated using terrestrial cosmogenic nuclide techniques. Although the terrestrial cosmogenic nuclide signature of the moraines is inconclusive, an early Weichselian age is tentatively suggested through correlations with other landforms and stratigraphical archives in the region. The abundance and coherent spatial pattern of the lateral moraines also allow a spatial reconstruction of this ice field. The ice field was controlled by topography and had nunataks protruding also where it was thickest close to the elevation axis of the Scandinavian mountain range. Outlet glaciers discharged into the Norwegian fjords and major valleys in Sweden. The process by which mountain based glaciers grow into an ice sheet is a matter of debate. In this thesis, a feedback mechanism between debris on the ice surface and ice sheet growth is presented. In essence, the growth of glaciers and ice sheets may be accelerated by an abundance of debris in their ablation areas. This may occur when the debris cover on the glacier surface inhibits ablation, effectively increasing the glaciers mass balance. It is thus possible that a dirty ablation area may cause the glacier to advance further than a clean glacier under similar conditions. An ice free period of significant length allows soil production through weathering, frost shattering, and slope processes. As glaciers advance through this assemblage of sediments, significant amounts of debris end up on the surface due to both mass wastage and subglacial entrainment. Evidence that this chain of events may occur, is given by large expanses of hummocky moraine (local name Veiki moraine) in the northern Swedish lowlands. Because the Veiki moraine has been correlated with the first Weichselian advance following the Eemian, it implies a heavily debris charged ice sheet emanating from the mountain range and terminating in a stagnant fashion in the lowlands.

Ice sheet modeling is an important tool both for reconstructing past ice sheets and predicting their future evolution, but is complex and computationally costly. It involves modeling a system including the ice sheet, ice shelves and ice streams, which all have different dynamical behavior. The governing equations are non-linear, and to capture a full glacial cycle more than 100,000 years need to be simulated. To reduce the problem size, approximations of the equations are introduced. The most common approximation, the Shallow Ice Approximation (SIA), works well in the ice bulk but fails in e.g. the modeling of ice streams and the ice sheet/ice shelf coupling. In recent years more accurate models, so-called higher order models, have been constructed to address these problems. However, these models are generally constructed in an ad hoc fashion, lacking rigor. In this thesis, so-called Second Order Shallow Ice Approximation (SOSIA) equations for pressure, vertical shear stress and velocity are implemented into the ice sheet model SICOPOLIS. The SOSIA is a rigorous model derived by Baral in 1999 [3]. The numerical solution for a simple model problem is compared to an analytical solution, and benchmark experiments, comparing the model to other higher order models, are carried out. The numerical and analytical solution agree well, but the results regarding vertical shear stress and velocity differ from other models. It is concluded that there are problems with the model implemented, most likely in the treatment of the relation between stress and strain rate.

The Antarctic Ice Sheets responded significantly to climatic conditions during the Last Glacial Maximum (LGM) and the subsequent warming that followed. Therefore, an understanding of how Antarctica reacted to past climates is necessary to predict the response of its ice sheets to current and future climate change. This thesis presents new evidence about the timing and magnitude of East and West Antarctic ice sheet (EAIS & WAIS) changes during the Quaternary Period, from the Darwin Hatherton glacial system (DHGS, 79.5S, 158E). The DHGS drains the EAIS through the Transantarctic Mountains into the Ross Ice Shelf and glacial deposits have been used to constrain ice sheet thicknesses in this sector of the Ross Sea Embayment. At four sites along the length of the system, glacial deposits were mapped and 73 erratic and bedrock samples collected for ¹⁰Be and ²⁶Al surface exposure dating (SED). The exposure ages range from 0.01 to 2.2 Ma and generally show a trend of oldest ages at the highest elevations, thus suggesting an overall decrease of ice volume within the DHGS over the Quaternary. The older ages suggest that during the Plio-Pleistocene, DHGS ice was at least 800-1000 metres thicker than present, while in the mid to late-Holocene thickening was less than 50-80 metres. Four glacial advance and retreat events were described and mapped previously from the DHGS by Bockheim et al (1989). The Isca and Danum drifts, are ~1-2 and 0.6 Ma respectively. The Britannia-II Drift, previously assumed to mark the maximum extent of the Last Glacial Maximum advance is more complex, with clusters of ages at ~6.5, ~36 and ~125 ka. The youngest drift, the Hatherton is mid to late-Holocene (<4.5 ka) and suggests that the DHGS has been near its equilibrium position during this period. Throughout the DHGS no unequivocal evidence of the LGM was observed and therefore poses questions about the past thickness of the Antarctic ice sheets during the LGM. Exposure ages from sites near the head of the Hatherton Glacier (Dubris Valley & Lake Wellman) suggest that at the LGM, the East Antarctic Ice Sheet may have been of similar size, or slightly smaller, than present. In stark contrast, at the confluence of the Darwin Glacier and the Ross Ice Shelf, a WAIS ~400-900 metres above the modern ice surface is tentatively suggested; A value in agreement with that proposed by modern glaciological models. Additionally, while the results from the Dubris and Bibra valleys show that the EAIS thins during glacial climates (i.e. the LGM), it also suggests thickening during interglacials. The Britannia-I and II drifts representing retreats at~6.5 and ~125 ka from glacial highstands. A number of key findings related to the application of SED in Antarctic settings are also presented. The use of dual-nuclides (¹⁰Be & ²⁶Al) show that within the DHGS, the proportion of samples displaying a prior burial history increases with distance from the catchment. The spread of exposure ages observed in the dataset also show the complexity of the depositional processes occurring at cold-based glacial margins and therefore judicious sample selection is required to obtain exposure ages that are representative of the true deglaciation age.

The continental margin in the area west of Shetland was subjected to repeated and extensive ice sheet advances during the Late Pleistocene. Seabed imagery, seismic survey and borehole core data show the Late Devensian ice sheets expanded across the continental shelf three times, two of these advances reaching the shelf edge. On the inner shelf, where present-day water depths are generally less than 100m, only thin sediments from the last retreat phase and exposed rock surfaces remain, all other deposits from earlier phases having been removed by the last advance. On the mid to outer shelf elements of all three phases are preserved, including lodgement and deformation tills, melt-out and water-lain till sheets, in-filled hollows left by stagnant ice decaying in situ and a series of large recessional and terminal moraines. In addition, there is evidence of shallow troughs and overdeepend basins which indicate preferential ice-drainage pathways across the shelf which were formerly occupied by ice streams. At the shelf edge, a thick wedge of glacigenic sediment forms a transition from the till sheets and moraines of the shelf to debris flows composed of glacigenic sediments on the upper slope. Shelf-edge moraines show an architecture indicating floating ice in modern water depths over approximately 180m, suggesting the West Shetland ice sheet was no more than about 250m thick. The upper and middle slope is dominated by glacigenic debris flows which are focused in the slope areas below the proposed ice stream discharges at the shelf edge. The mid-to-lower slope has been subjected to contour current activity which has re-worked much of the glacigenic sediment in this position. The lower slope and floor of the Faroe-Shetland Channel are marked by either large debris flow lobes of glacigenic sediment or thin glacimarine muds deposited from suspension. A conceptual model of the glacigenic development of a passive continental margin based upon the West Shetland example shows the deposited sequence for both advance and retreat phases of a glacial cycle, and the actual preserved sequence which might be expected in the rock record. The model also shows that ice sheet buoyancy, thickness, and to a lesser extent, basin subsidence, are the most important factors in the deposition and preservation of a glacially-influenced marine sequence.

Glaciers and ice sheets covered extensive areas in the Northern Hemisphere during the last glacial period. Subsequently to the Last Glacial Maximum (LGM), they retreated rapidly and, except for Greenland and some other ice caps and glaciers, they vanished after the last glacial termination. This thesis examines the dynamics of deglacial environments by analysing the glacial geomorphological record with focus on the landforms created by glacial meltwater. The aims are (i) to evaluate the data available for mapping glacial meltwater features at the regional scale, and (ii) to demonstrate the potential of such features for regional ice retreat reconstructions in high-relief landscapes. Meltwater landforms such as ice-marginal meltwater channels, eskers, deltas and fossil glacial lake shorelines are used to infer former ice surface slope directions and successive positions of retreating ice margins. Evaluated high-resolution satellite imagery and digital elevation models reveal their potential to replace aerial photographs as the primary data for mapping glacial meltwater landforms. Following a methods study, reconstructions of the deglacial dynamics are carried out for central Transbaikalia, Siberia, Russia, and for the Cordilleran Ice Sheet (CIS) in central British Columbia, Canada, using regional geomorphological mapping surveys. Mapped glacial landforms in central Transbaikalia show evidence of a significant glaciation that possibly extended beyond the high mountain areas. Large glacial lakes were formed as advancing glaciers blocked rivers, and of these, Glacial Lake Vitim was the most prominent. Deglacial dynamics of the CIS reveals that the ice divide shifted to the Coast Mountains in north-central British Columbia and the eastern ice margin retreated towards the ice divide in late glacial time. This thesis demonstrates the potential to reconstruct ice retreat patterns and deglacial dynamics at regional scales by interpretation of the meltwater landform record. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 5: Submitted. Paper 6: Manuscript.

In this thesis, I explore the flow dynamics associated with ice shelves confined within channels and the buttressing they provide to grounded ice. Ice shelves are the floating extensions of ice sheets and act as the interface between the ice sheet and the ocean. They form when ice flows out from the interior of the ice sheet towards the coast and begins to float as the ice thins. Ice shelves are often found within a channel or pinned in place by stationary bedrock outcrops. The interest in their dynamics is motivated by the buttressing effect they provide to the grounded ice, which strongly controls the rate of ice discharge and thereby the contribution to sea-level rise. I use a combination of mathematical modeling, fluid-mechanical laboratory experiments and geophysical data analysis to develop an improved understanding of ice-shelf flow dynamics. Initially, geophysical data in the form of Antarctic ice-surface velocity data is analysed, producing maps of strain rate, shear rate and strain orientation for Antarctic ice shelves. This allows the geophysical setting and flow processes to be explored, particularly by identifying areas where resistance to ice flow is generated and regions of the shelf that make no contribution to buttressing. Using the geophysical data, I find good agreement between a theoretical scaling relationship for ice flow at the ice-shelf calving front and data from Antarctic ice shelves. I proceed to develop an idealized mathematical model of an ice shelf confined to flow in a channel. By assuming shear-dominated dynamics within the shelf, analytical solutions are obtained for steady-state ice-shelf thickness profiles in parallel and diverging channels. This model is developed further to include both shear and extensional stresses, from which numerical solutions for steady-state shelves are calculated. The results from these two models are then compared. It is found that shear stresses dominate the dynamics throughout the majority of the shelf, with adjustment regions at the upstream and downstream boundaries where extensional dynamics become important. Output from these models is also compared with geophysical data and it is observed that there is good agreement between several features of the thickness profiles and velocity fields. In addition to the geophysical data, comparisons are made with fluid-mechanical laboratory experiments designed to simulate the flow of an ice shelf in a channel. The advantage of performing experiments of this kind is that parameters such as the fluid rheology can be varied, allowing for direct comparison with a range of parameters in the mathematical models. From these experiments, surface velocity fields and thickness profiles are collected, which are used to make comparisons with the models. Clear differences are observed in the velocity and strain-rate fields produced using fluids with different rheologies, for which there is qualitative agreement with the output from the mathematical models.

Despite their central importance to the Earth system, nowhere within the literatures of Polar Studies or the Humanities does there exist a comprehensive cultural and scientific history of ice sheets that takes into consideration both Greenland and Antarctica, or that is not constrained to a particular exploratory, technological, or geopolitical period or framing. My thesis addresses this lacunae by contributing a bi-polar, empirical history and analysis of the different scientific and cultural processes, transformations, and discontinuities through which ice sheets have been transformed from unexplored realms of the imagination, into tangible, material objects of investigation and meaning. Its scope extends from early Greek mapping to 1970. Within this timeframe, it identifies three broad phases in the perception of ice sheets. The first, preceding their earliest physical exploration, corresponds to the perception of ice sheets as one-dimensional realms defined and bounded by the human imagination. The second phase, associated with their early surface exploration between 1870 and 1930, corresponds to the perception of ice sheets as undifferentiated, two-dimensional 'topographies of absence', best characterized by their horizontal desolation. The third phase, triggered by the deployment of new technologies of sub-surface investigation such as seismic sounding, radio echo sounding (RES), and the practice of ice coring, corresponds to the perception of ice sheets as three-dimensional, super-massive, and interdependent objects of internal and material complexity. Although primarily rooted in archival research and the study of first hand textual and visual materials, my arguments and observations also draw on secondary literatures from the history of science and technology, geopolitics, visual culture, and the geography of space and place. These literatures allow me to contextualize and substantiate my analysis of historical processes within broader perspectives, notably Humboldtian science, Romanticism, visual abstraction, scientific imagery, and the Cold War.

There is increasing recognition that the hydrology of the Greenland Ice Sheet plays an important role in the dynamics and therefore mass balance of the ice sheet. Understanding the hydrology of the ice sheet and being able to predict its future behaviour is therefore a key aspect of glaciological research. To date, the ice sheet’s hydrology has tended to be inferred from the analysis of surface velocity measurements, or modelled in a theoretical, idealised way. This study focuses on the development of a high spatial (100 m) and temporal (1 hour) resolution, physically based, time-dependent hydrological model which is applied to the ~2,300 km2 Paakitsoq region, West Greenland, and is driven, calibrated, and evaluated using measured data. The model consists of three components. First, net runoff is calculated across the ice sheet from a distributed, surface energy- balance melt model coupled to a subsurface model, which calculates changes in temperature, density and water content in the snow, firn and upper-ice layers, and hence refreezing. The model is calibrated by adjusting key parameter values to minimize the error between modelled output and surface height and albedo measurements from the three Greenland Climate Network (GC-Net) stations, JAR 1, JAR 2 and Swiss Camp. Model performance is evaluated in two ways by comparing: i) modelled snow and ice distribution with that derived from Landsat-7 ETM+ satellite imagery using Normalised Difference Snow Index (NDSI) classification and supervised image thresholding; and ii) modelled albedo with that retrieved from the Moderate- resolution Imaging Spectroradiometer (MODIS) sensor MOD10A1 product. Second, a surface routing / lake filling model takes the time-series of calculated net runoff over the ice sheet and calculates flow paths and water velocities over the snow / ice covered surface, routing the water into ‘open’ moulins or into topographic depressions which can fill to form supraglacial lakes. This model component is calibrated against field measurements of a filling lake in the study area made during June 2011. Supraglacial lakes are able to drain by a simulated hydrofracture mechanism if they reach a critical volume. Once water is at the ice / bed interface, discharge and hydraulic head within subglacial drainage pathways are modelled using the third model component. This consists of an adaptation of a component (EXTRAN) of the U.S. Environmental Protection Agency Storm Water Management Model (SWMM), modified to allow for enlargement and closure of ice-walled conduits. The model is used to identify how the subglacial hydrological system evolves in space and time in response to varying surface water inputs due to melt and lake drainage events, driven ultimately by climate data. A key output from the model is the spatially and temporally varying water pressures which are of interest in helping to explain patterns of surface velocity and uplift found by others, and will ultimately be of interest for driving ice dynamics models.

Indiana University-Purdue University Indianapolis (IUPUI) / A broad range of samples were collected from the Huron-Erie Lobe, Lake Michigan Lobe, Saginaw Lobe, and Tipton Till Plain of northern Indiana to determine the provenance of Laurentide Ice Sheet till in the Midwest U.S. during the Illinoian and Wisconsinan glaciations. U-Pb age distributions from approximately 300 detrital zircons (DZ) were used as provenance indicators for each till sample. Till from the Lake Michigan Lobe and was found to be largely homogenized. The distinct lobe DZ age distributions are the Lake Michigan Lobe till with a dominant ~1465 Ma peak, the northern Huron-Erie Lobe till with a dominant ~1060 Ma and a secondary peak at ~1450 Ma, the southern Huron-Erie Lobe till with nearly equal peaks at ~1435 Ma, ~1175 Ma, and ~1065 Ma, and the southern Saginaw Lobe till with a dominant peak at ~1095 Ma. Those four DZ age distributions were treated as endmembers in a nonlinear least-squares mixing model to calculate the contribution of each lobe to till in the Tipton Till Plain. Huron-Erie and Saginaw lobe tills were found to be the primary components of the Tipton Till Plain, and Lake Michigan Lobe till was only found in the western Tipton Till Plain. Zircons from the Saginaw Lobe till increased 39 % in the eastern Tipton Till Plain between the Illinoisan and Wisconsinan glaciations. The mixing model was also applied to relate the DZ age distributions of the lobes to bedrock within and near their flow paths. When comparing nearby bedrock to each lobe’s till, mixing model results, yield an approximate maximum transport distance between 500 and 630 kilometers for the matrix vii fraction of till in the Lake Michigan, Huron-Erie, and Saginaw lobes. Samples for the southern Huron-Erie Lobe indicate that the most of the zircon ages within the southern Huron-Erie Lobe till in Indiana were specifically entrained between Niagara County, New York and east-central Indiana. Within the model’s error, 93 – 100 % of the detrital zircons in each of the three lobes are relatable to nearby Paleozoic and Precambrian sedimentary and metamorphic bedrock formations.