Introduction
Although glaciers have long been recognized as significant agents of erosion and deposition (e.g. Reference PlayfairPlayfair, 1802; Reference AgassizAgassiz, 1840; Reference RamsayRamsay, 1862; Reference GilbertGilbert, 1905), little is known about the rates at which debris is transported and released at their termini. Modern and ancient moraine comlexes, widespread till sheets, glacial outwash plains and massive offshore deposits attest to the effective mechanisms of glacial sediment transport. Modern process studies demonstrate that the accumulation of ice-proximal landforms varies according to the balance between the glacier sediment yield and the processes that control sediment deposition and redistribution near the terminus (e.g. Reference LawsonLawson, 1979; Reference Powellpowell, 1981; Reference HunterHunter 1994). It is now clear that this sediment balance influences ice-cliff stability in the marine environment (Reference MayoMayo, 1988; Reference AlleyAlley, 1991b; Reference PowellPowell, 1991; Reference Hunter, Powell and EngstromHunter and Powell, in press b) by moderating grounding-line water depth (Reference Brown, Meter and PostBrown and others, 1982; Reference Pelto and Warrenpelto and Warren, 1991). Therefore, an understanding of sediment transport and depositional processes not only helps us understand the formation of sedimentary facies, but should also be considered when studying the dynamics of tidewater glaciers.
Variables affecting glacial debris flux include source-rock type, relief, weathering state and tectonic setting. Climate influences precipitation and temperature, which in turn influences the ice flux, thermal conditions and debris entrainment processes. Our study was located in the tectonically active region of Glaciers Bay, Alaska, U.S.A. (Fig. 1), where friable bedrock and unstable valley walls are the result of regional metamorphism, local metasomatism, deforamtion and fracturing (Reference PowellPowell, 1984b). Mountains in Glacier Bay reach elevations of over 4660 m a.s.l., where orographic precipitation and winds develop thick, often drifted, snow covers that nourish numerous glaciers (Reference HunterHunter, 1994). Local bedrock conditions combine with rapidly flowing glaciers to produce high sediment yields. Reference Hallet, Hunter and BogenHallet and others (in press) demonstrated recently that sediment yields iin Glacier Bay correspond to effective denudation rates ranging from about 30 to 60 mm a −1. Such trates are the highest yet reported in the literautre for glacierized and non-glacierized basins (Reference Milliman and SyvitskiMilliman and Syvitski, 1992; Reference Hallet, Hunter and BogenHallet and other, in press). Glacier Bay erosion rates are about an order of magnitude greater than those reported from southern Iceland (3.8 mm a1: Reference Lawler, Dolan, Tómasson and ZóphóníassonLawler and others, 1992), Punjab Himalaya (1.4–2.1 mm a−1: Gardner and Jones, 1984), Central Asia (3–5mm a−1: Reference ChernovaChernova, 1981) and the Ivory Basin, New Zealand (5.6 mm a−1: Reference Hicks, McSaveney and ChinnHicks and others, 1990), and two orders of magnitude grater than Scandinavian glaciers (∼0.1−0.5mm a−1: Reference BogenBogen, 1989).
Map of Glacier Bay National Park and Preserve, showing the location of (1) Muir Glacier and (2) Grand Pacific and Margerie Glaciers.
Sediment released at the terminus is transported in the glacier by turbulent englacial and subglacial streams (e.g. Reference Gustavson and BoothroydGustavson and Boothroyd, 1987; Reference LawsonLawson, 1993) and as debris in the ice. Debris transported in the basal and englacial zones is released through ablation at the terminus or grounding line, while supraglacial debris is mainly dumped during calving events (Reference Powellpowell, 1981). Some basal and englacial debris may be deposited away from the grounding line following transport by icebergs downfjord (e.g. Reference OvenshineOvenshine, 1970; Reference Dowdeswell and DowdeswellDowdeswell and Dowdeswell, 1989; Reference Dowdeswell, Murray, Dowdeswell and ScourceDowdeswell and Murray, 1990; Reference Gottler, Powell, Milner and WoodGottler and Powell, 1990). Additional sediment transport may occur subglacially in a deforming layer (e.g. Reference Boulton, Menzies and RoseBoulton, 1987; Reference Boulton and HindmarshBoulton and Hindmarsh, 1987; Reference AlleyAlley, 1991a; Reference Humphrey, Kamb, Fahnestock and EngelhardtHumphery and others, 1993); but monitoring of this layer and measurement of its flux are difficult due to its inaccessibilty.
A tidewater terminus is ideal for debris-flux analysis since an ice cliff represents a vertical cross-section through the glacier. Icebergs calving from above and below sea level (Fig. 2) eventually drift away from the ice cliff where they can be sampled. Ice from each stratigraphic horizon is collected to develop a cmposite stratigraphic section through the glacier. Selective sampling allows us to concentrate on defining the spatial distribution of debris so that the debris flux can be calculated. Our analysis differs from previous debris studies that have focused on demonstrating processes of debris entrainment and ice-facies formation (e.g. Reference LawsonLawson, 1979; Reference KnightKnight, 1989; Hubbard, 1991; Reference GottlerGottler, 1992; Reference Ronnert and MickelsonRonnert and Mickelson, 1992; Reference GoodwinGoodwin, 1993; Reference Sharp, Jouzel, Hubbard and LawsonSharp and others, 1994; Reference SouchezSouchez and others, 1994).
Glacier calving process in Glacier Bay: (a) subaerial calving at Johns Hopkins Glacier, in which ice shatters upon impact with the water surface, and (b) submarine calving at McBride Glacier. Submarine icebergs lend to remain more-or-less intact on reaching the water surface (photo from McBride Glacier by K. Seramur).
In this paper we describe the debris distribution in three tidewater glaciers (Grand Pacific, Margerie and Muir) in Glacier Bay, Alaska (Fig. 3), and evaluate their sediment yields, emphasizing the distribution of debris in and on the glaciers. Debris-flux estimates are useful for interpreting glacimarine records of the Gulf of Alaska (e.g. Reference Plafker, Addicott and MillerPlafker and Addicott, 1976; Reference Eyles, Eyles and LagoeEyles and others, 1991) and perhaps for evaluating the minimum volume of ice required to distribute the large quantities of iceberg-rafted debris associated with Heinrich events in the North Atlantic (e.g. Reference Dowdeswell, Maslin, Andrews and McCaveDowdeswell and others, 1995).
UIL high-altitude aerial photographs of (a) Muir Glacier and (b) Grand Pacific, Ferris, and Margerie Glaciers, taken in Angust l979. The number of moraines reflects tributary confluence where debris reaches englacial-transport pathways in medial moraines that project into the glacier (courtesy of the U.S. Geological Survey ICe and Climate Project, Tacoma, Washington).
Terminology
Terms used in this paper have been chosen for the purpose of modeling sedimentary processes along the marine grounding line of a tidewater terminus. A tidewater terminus is a grounded vertical ice cliff of a glacier exposed to tidal activity (Reference PowellPowell, 1988). We use grounding line to describe the boundary at the base of the ice cliff where it is in contact with the sea floor. Since sedimentary processes at the grounding line form a complex depositional system that extends away from the grounding line (cf. Reference Powellpowell, 1981, Reference Powell1984b), it is appropriate to describe the depositional environment as the grounding-line system (Reference PowellPowell, 1988). Sediment piles that accumulate in a bank geometry as part of grounding-line systems are called morainal banks (Reference Powellpowell, 1981, Reference Powell1984a).
The glacifluvial and glacial components of the sediment budget are differentiated by reserving the term sediment for lithic material transported by fluvial or marine processes, the latter including debris once it has been released by melt-out or calving. Lithic material in transport in and on top of a glacier are referred to as debris (Reference LawsonLawson, 1979). Therefore we describe glacial sediment flues as glacifluvial sediment flux (sediment flux) and glacial debris flux debris flux). Combined, the two fluxes represent the total glacier sediment yield of the basin, which we term total glacial sediment flux.
Sampling Procedures
Debris concentrations were calculated for 282 iceberg samples and 139 glacier-ice samples (Fig. 4). The percentage of each facies exposed at each terminus was estimated in the field using a comparison chart for percentage composition (AGI, 1982) and subsequently by point-count analyses of photographs, The thickness of basal stratigraphic sequences was determined from large icebergs that originate at the grounding line and from measurements at terrestrial exposures along these margins. Calving events originating below the water-line (Fig. 2b) produce large icebergs that do not fragment like those that topple or free-fall from subaerial parts of the ice cliff (Fig. 2a). These large icebergs often contain thick basal ice sequences that can be measured and sampled.
Sites of iceberg and glacier-ice samples from (a) Grand Pacific and Margerie Glaciers, and (b) Muir Glacier. B, iceberg sample; S, glacier ice sample at terminus; T, glacier-ice sample collected along transects up-glacier from terminus.
Englacial and basal ice samples were collected at exposures along the terminus and supraglacially along transects near the terminus (Fig. 4). Weathered surface ice was removed using an ice axe to expose fresh ice prior to the collection of each sample. All iceberg and glacier-ice samples were then described, crushed and melted before their volumes were measured (for details, see Reference HunterHunter, 1994). Samples were then allowed to settle before the meltwater was decanted and the sediments shipped to the laboratory for drying and weighing.
Debris concentration (C; in volume fraction) of each ice sample was calculated by:
where the ice volume (V i) is
using 1.089 to adjust for the reduction in volume when ice melts, V t is the total sample volume, and V d is the debris volume defined by
where m is the mass and ρd is the density of the glacial debris as measured in the laboratory (2560 ± 109 kg m−3 at Grand Pacific Glacier. 2870 ± 511 kg m −3 at Margerie Glacier and 2800 ± 264 kg m −3 at Muir Glacier).
Ice-Facies Analyses
Debris-flux calculations depend on the identification of the characteristic ice facies that Compose the glacier, measurements of debris concentrations, and appraisal of the percentage of each ice Facies in the glacier. We use an ice-facies classification modified from Reference LawsonLawson (1979) that includes additional ice facies exposed in tidewater ice cliffs associated with medial moraines (Fig. 5). Lawson’s original classification distinguishes ice facies by styles of debris suspension, concentration and zones of transport (englacial or basal). His classification consists of two englacial facies, the diffuse and banded, and two basal facies, the dispersed and Stratified. The diffuse ice facies forms the greatest volume of the glacier and consists of white to light-blue, bubble-rich ice with an extremely low debris content (<1%). Debris is suspended in the ice as isolated mud aggregates or silt and sand particles. Thin horizons with greater debris concentrations are called the banded ice facies and form a minor constituent of the glacier. We add a diffuse-stratified facies to describe pervasively foliated ice observed near the lateral margins of tidewater termini or in association with narrow bands of englacial ice confined by medial moraines.
General ice-facies distribution at a tidewater terminus (modified after Reference Gottler, Powell, Milner and WoodGotter and Powell (1990)).
For the debris-flux analyses, we assume that ice facies in the basal and morainal zones share similar processes of formation in the subglacial environment. Therefore, our distinction is based on glacier stratigraphy, where ice in the morainal zone reaches an englacial-transporl path-way by tributary confluence and the formation of medial moraines. The dispersed facies is generally transitional between englacial diffuse ice and the debris-rich horizons of the basal and morainal facies, and consists of blue ice that is generally bubble-free to bubble-poor, with debris concentrations slightly higher (mean about 4–5%) than englacial ice (Reference Lawson and KullaLawson and Kulla, 1978; Reference LawsonLawson. 1979). The stratified ice facies is debris-rich and forms the lowermost horizons in the basal zone and the core of moraines. This facies is divided into three subfacies based on the distribution of sediment in the ice. The suspended subfacies consists of relatively clear, bubble-free ice within which are suspended sand, silt, pebble-size particles and aggregates. The discontinuous subfacies consists of thin, discontinuous debris-rich lenses and laminae separated by thin clear-ice laminae. The solid subfacies consists of frozen debris held together by interstitial ice, and may occasionally exhibit primary sedimentary structures.
Englacial Ice Facies
The diffuse facies is dominated by white, bubbly ice that comprises 65–83% of the ice in the termini of Grand Pacific, Margerie and Muir Glaciers (Fig. 6a; Table 1). Debris occurs as isolated sand and silt particles, silty mud-filled bubbles and rare cobbles and boulders with concentrations in the diffuse facies ranging from trace values to 0.9% (Reference HunterHunter, 1994). The highest average debris concentration of diffuse facies occurs at Muir Glacier (0.6% by volume), whereas average concentrations at Mangerie and Grand Pacific Glaciers are about 0.2% and 0.1% by volume, respectively (Table 1).
Ice facies observed at tidewater termini in Glacier Bay: (a) bubbly englacial diffuse, (b) englacial diffuse-stratified. (c) basal dispersed, (d) basal discontinuous subfacies, (e) basal solid subfacies. and (f) alternating layers of debris-rich discontinuous subfacies and debris-poor suspended subfacies.
Summary of average debris calculations. x , average debris concentration with 1 σ in kg m 3; %, representative percentage of each ice facies in ice cliff; Vol, volume fraction of debris in concentration
Diffuse-stratified facies occurs in zones of pervasively foliated englacial ice (Fig. 6b). Bulk sampling yielded debris concentrations of 0.4–14.4% by volume. The diffuse-stratified facies is most abundant at Grand Pacific Glacier where it comprises up to 7.1% of the terminus and has an average debris concentration of 2.7% by volume (Table 1). At Margerie and Muir Glaciers, this facies is a minor component (<2%) of the ice cliff and contains up to 1.4% debris by volume.
Banded facies constitutes a small percentage (<%; Table 1) of ice cliffs and is characterized by thin debris-rich horizons. This facies is poorly developed at Margerie and Muir Glaciers where thin (up to a few millimeters) horizons have debris concentrations of 0.2–0.4% by volume. At Grand Pacific Glacier, similar horizons are up to 3cm thick with concentrations ranging from 1.7% to 27.1 % by volume. Debris horizons are commonly stratified, consisting of debris-rich bands separated by clear-ice layers. Debris in these hands is in the form of suspended sand-particle and mud clots (Reference HeinyHeiny, 1983). At higher concentrations, the debris occurs in lenses and more continuous strata with only interstitial ice.
Basal-Zone and Morainal Ice Facies
Morainal and basal-ice stratigraphic sequences are composed of dispersed and stratified facies. Morainal sequences contain ice facies similar to those in the basal zone, but the facies are often repeated on both sides of a central debris-rich horizon, and the dispersed facies may be absent in morainal sequences. Morainal sequences appear to represent condensed basal sequences (Reference LawsonLawson, 1979) that range from a few meters to much thicker sequences, as at Grand Pacific Glacier where the sequence marking the boundary of the Ferris tributary is 30 m wide and includes a 1 m thick debris-rich horizon of Reference LawsonLawson’s (1979) solid subfacies.
The dispersed facies is characterized by low average debris concentrations ranging from about 0.1% to 1.4% by volume (Table 1). A sharp contact between the englacial diffuse and basal dispersed ice (Fig. 6c) is generally discordant with glacial sedimentary layering and ice-banding. Debris is dispersed in the form of isolated particles and aggregates that occur predominantly at grain boundaries.
Basal and morainal stratified facies constitute 4.3 10.4% of the ice observed (Table 1). The discontinuous subfacies is characterized by alternating thinly laminated debris-rich and clear ice (Fig. 6d). Debris horizons range from a few sand grains to several centimeters thick, and form lenses of various lengths in each debris-rich horizon. Clear-ice laminae are bubble-free, with minor amounts of suspended silt and sand particles.
Solid subfacies are generally devoid of clear-ice laminae and are composed of debris horizons greater than 1 cm thick. Interstitial ice-crystal growth is confined in the fine pore structure of silty muds that exhibit a frosty (ice-cement) luster on freshly broken surfaces. Reference BoultonBoulton (1970) and Reference LawsonLawson (1979) have described sedimentary structures in solid layers that indicate possible incorporation at the glacier sole by freeze-on and vertical accretion (Reference Lawson and KullaLawson and Kulla, 1978; Reference LawsonLawson, 1979) or perhaps bv fluvial sediments deposited in basal-zone conduits.
Thick horizons of solid subfacies have been observed at the stratigraphic bases of icebergs that originated below sea level in contact with the sea floor (Fig. 6e).
These horizons are composed of diamictons similar to those described by Reference Powell and MolniaPowell and Molnia (1989) and Reference Gottler, Powell, Milner and WoodGottler and Powell (1990), but they melt rapidly or break off icebergs and sink. Where solid layers have been observed in an iceberg, they form anisotropies along which icebergs tend to break. At Margerie Glacier, subglacial debris fills flutes (up to a few tens of centimeters deep) on the surface of icebergs that had previously been in contact with the sea floor (Fig. 7). Solid layers at Muir Glacier are composed of well-sorted and well-rounded sand and gravel apparently derived from fluvial sediments in the morainal bank (Reference HunterHunter, 1994). Solid subfacies resulting from seasonal freeze-on of saturated sediments have been observed in basal ice sequences ol Gulkana, Matanuska and Kennicott Glaciers in southern Alaska, and Findeln and Z’Mutt Gletchers in the Swiss Alps (Lawson, unpublished data). Solid subfacies contain up to 70% debris by volume (Table 1).
The suspended subfacies (Fig. 6f) consists of ice layers about 1cm to over 3 m thick. This subfacies was a minor constituent (2.6–5.7%: Table 1) of the three tidewater termini, with debris concentrations less than 0.1–12.6% by volume (Reference HunterHunter, 1994). Debris is suspended as isolated silt, sand and gravel particles, and commonly forms sandy-mud aggregates.
Glacier-Sole morphology preserved on a submarine iceberg that calved from Margerie Glacier. (a) Note linear fluting across face of berg, (b) Close-up of right edge of berg shown in (a), highlighting surface texture and row of cobbles that are aligned clast to clast.
Debris Flux of Grand Pacific, Margerie and Muir Glaciers
Compilation of the ice-facies data (Table 1) gives us the opportunity to calculate the total volume of debris transported to the terminus. The debris concentrations of each facies and their percentage in an ice cliff (Tables 1 and 2) are used to calculate the debris flux (D t) for the different glaciers In the simple relationship:
where D s is the supraglacial debris flux (i.e. surface transport) and D g is the debris flux within basal and englacial ice facies calculated using:
so that D g is the product of the ice flux (Q i) and the sum of the debris concentrations (C) of each ice facies weighted by their relative volume (Vj in the ice cliff. The ice flux is:
where ν t is the averaged surface velocity near the terminus as estimated from aerial photographs, ω is the glacier width, and h t is the height of the ice cliff (Table 3). This calculation assumes that ice flow near the glacier terminus occurs as a plug with little internal deformation. We assume that the observed surface velocity is due mostly to rapid basal sliding, as observed at another calving glacier in southern Alaska, Columbia Glacier (Reference Humphrey, Kamb, Fahnestock and EngelhardtHumphrey and others, 1993).
Calculated debris flux values *. TR, trace
Glacier parameters during study (1988–91)
Calculations of Qi for Margerie and Muir Glaciers yielded ice fluxes of 1.38 × 108 and 1.35 × 108 m 3a−1, respectively. The ice flux at Grand Pacific Glacier (0.62 × 108 m3 a−1) is half as large because of the glacier’s slower flow (Table 3). By solving Equation (5) for each glacier we estimate debris fluxes of 1.3 × 106 m3a−1 (Muir Glacier), 9.7 × 105m3a−1 (Margerie Glacier) and 1.1 × 106m3a−1 (Grand Pacific Glacier; Table 2).
Comparison of these estimates demonstrates the influence of debris distribution and ice-flow rates on the debris flux. Muir and Margerie Glaciers transport less debris per unit volume than Grand Pacific Glacier (Table 1);1 however, Muir Glacier has the highest debris flux, due to its large ice flux (Table 2). Alternatively, despite the low ice flux of Grand Pacific Glacier, its debris flux is about 30% greater than that of Margerie Glacier because of its higher debris concentrations. Margerie Glacier’s debris concentrations are about 60% per unit volume of Muir Glacier’s. The concentration and distribution of debris alone are not enough to produce a high debris flux, but when combined with rapid glacier flow, exceptionally high rates of debris transfer to the grounding line can occur.
Supraglacial debris flux was calculated as the product of glacier surface velocity, moraine widths measured photogrammetrically and surficial debris thickness measured in the field. Moraine thicknesses observed during transects of Margerie and Grand Pacific Glaciers (Fig. 4) ranged from less than 1 mm to 1.5 m (Reference HunterHunter, 1994), and were generally less than the 0.08 m estimate of Reference GottlerGottler (1992) at McBride Glacier. A debris cover of about 1 mm is sufficient to discolor the surface while a thickness of 1–2 cm appears to form a complete cover on aerial photographs (e.g. Fig. 3). Thicker moraine cover (0.5– 1.0 m) fills surface crevasses and forms debris ridges and rounded gravel surfaces. Despite the conspicuous appearance of supraglacial moraines (Fig. 8), their fluxes are relatively low: 3.4 × 104m3a−1 at Grand Pacific Glacier, 4.1 × 104 m3 a−1 at Margerie Glacier and 1.4 × 104m3a ×1 at Muir Glacier (Table 2). Supraglacial debris is assumed to be released entirely by various gravitational processes at tidewater ice cliffs; however, some is also released on the heads of fan deltas along terrestrial parts of the margins (e.g. the northern edge of Margerie Glacier and the eastern margin of Grand Pacific Glacier).
Oblique aerial photos of (a) Muir Glacier.
(b) Ferris tributary to Grand Pacific Glacier.
(c) Margerie Glacier. The number and complexity of tributaries, local relief and valley width and area-elevation distribution relative to zones of maximum snowfall influence debris distribution (see text for description: photographs by A. Post, U.S. Geological Survey).
Discussion
Debris-flux estimates provide some indication of the role various factors play in determining the transport of debris by tidewater glaciers. The differences in debris distributions and fluxes between Muir, Grand Pacific and Margerie Glaciers indicate that valley geometry, especially width, glacier flow speed and the number of tributary glaciers flowing into the main trunk glacier, strongly influence debris flux. Grand Pacific and Margerie Glaciers are larger and wider than Muir Glacier (Table 3; Figs 3 and 8); the apparent result is that debris is more easily distributed across Muir Glacier by slope mass movement and eolian processes. While 52% of Muir Glacier’s debris is transported englacially. only 18% and 19% of the debris is transported by the equivalent ice facies at Grand Pacific and Margerie Glaciers, respectively. Debris concentration in the diffuse facies of Muir Glacier is more than four times those of Grand Pacific and Margerie Glaciers (Table 1). We believe that the debris distribution at Muir Glacier reflects the relative importance of lateral debris sources (e.g. rockfalls, slides and avalanches) in the narrow valleys through which Muir Glacier flows (Figs 3a and 8a). These soures introduce supraglacial debris to the surface of the glacier where it becomes entrained by snow accumulation and ice metamorphism.
The relative abundance of debris in basal and morainal facies at Margerie and Grand Pacific Glaciers documents the influence of tributary confluences. Dominant zones of debris transport in Grand Pacific are about equally subdivided between morainal (39%) and basal ice facies (43%). The relatively high contribution of morainal facies must reflect the complex interactions of its multiple tributaries, especially Ferris Glacier (Figs 3b and 8b). In sharp contrast, Margerie Glacier (Figs 3a and 8c) transports nearly 74% of its debris in the basal ice. Its low englacial-debris flux (Table 2) reflects the influence of fewer tributaries and the wide valley of its upper accumulation area (Reference Hunter, Powell and EngstromHunter and Powell, in press a).
Sediment Introduction into the Fjord Environment
Debris transported to the terminus of tidewater glaciers is released by a variety of processes (Reference PowellPowell, 1980, Reference Powell1988; Reference HunterHunter, 1994) including grain-by-grain melt-out at the ice cliff (ice-cliff melt-out, D m), dumping of supraglacial debris during calving (calve dumping, D d) and melt-out from icebergs (iceberg-rafting. D i), or ablation and melt-out at terrestrial margins (D ∆).
The debris flux can be partitioned as follows:
where D d and D m add debris to the morainal bank along the grounding line. D ∆ is the debris dumped on to terrestrial parts of the morainal bank (e.g. deltas) and valley sides, and D i represents debris transport away from the terminus.
From our observations at both terrestrial and tide-Water termini, we feel that the lowermost solid facies represents localized seasonal freezing of subglacial sediment, or rapid freezing due to thermodynamic disequilibrium near calving margins. Such layers would experience only local transport from the subglacial slope of a morainal bank to the grounding line. In this regard, we do not treat these layers as adding sediment to the morainal bank, but simply as a form of recycling (freeze-recycling: R f). The flux of sediment recycled by this mechanism is estimated by sampling the lowermost solid layer along terrestrial sections of termini and similar layers carried to the fjord-water surface b) icebergs (Fig. 3b) that were in contact with the fjord floor prior to calving (Fig. 7: Reference HunterHunter, 1994). Solid layers stratigraphically higher in the basal ice are, however, added to the calculation of the debris flux (Equation (7)).
Using each glacier’s attributes (Table 3), we can calculate the volume of debris released into the fjord environment by these processes. The calving speed (νc) is calculated using a modified continuity equation (Reference Brown, Meter and PostBrown and others, 1982):
where X is the rate of change of glacier length (positive for advance; Reference MeierMeier and others, 1980) and R is the melting rate (see Table 3) calculated using the Reference Weeks and CampbellWeeks and Campbell (1973) equation:
where ν is the boundary-layer water velocity. ΔT is the temperature difference between water and ice, and l is a linear measure of the ic e cliff in contact with water, either its height in contact with the sea (i.e. water depth, h w; Reference Powell and MolniaPowell and Molnia, 1989: Reference SyvitskiSyvitski, 1989) or width (W c; Reference Powell and MolniaPowell, 1983, Reference Powell1988). An average temperature difference (ΔT) of 2.95°C was measured by a remotely operated submersible that collected data in front of tidewater termini in Glacier Bay (personal communication from S. Bograd, 1994).
The choice of characteristic length (l in Equation (9)) depends upon whether the dominant current at the ice cliff is from buoyant upwelling of meltwater or horizontal fjord currents. Melting rates were calculated by substituting h w for l for Grand Pacific and Muir Glaciers where upwelling along the ice cliffs is important and wc for Margerie Glacier where tidal currents across the face prevail. Open water (2–5 m wide) is common along the ice cliffs of Grand Pacific and Muir Glaciers that indicate the spreading of upwelling fresh water that displaces longitudinal currents. In addition, fjord circulation at the terminus of Muir Glacier was restricted by a shallow morainal-bank sill left behind by grounding-line retreat (Reference HunterHunter, 1994). In contrast, Margerie Glacier terminates in a hanging valley along the western margin of Tarr Inlet (Fig. 1) where its ice cliff is exposed to tidal currents moving up- and downfjord (e.g. Reference Syvitski, Burrell and skeiSyvitski and others, 1987). Turbid surface water flows rapidly across the face of Margerie Glacier and forms strong near-surface (<10 m depth) currents while tidal exchange produces currents at depth.
On this basis, melting rates of 21 and 20 m a−1 have been calculated for Grand Pacific and Muir Glaciers, respectively, using a boundary-layer velocity of 0.03 ms−1 to represent buoyant upwelling following the arguments of Reference Cowan and PowellCowan and Powell (1986) and Reference Powell and MolniaPowell and Molnia (1989), and oceanographic measurements by Reference Mathews and QuinlanMathews and Quinlan (1975). Using terminus width exposed to tidewater (ωc) and a velocity of 0.25 ms−1 (estimated From iceberg drifting), a melting rate (R m) of 31 m a−1 has been calculated for Margerie Glacier.
We can calculate the melting flux (Q m) and debris volume released due in melting Dm by substituting melting rate R for velocity νt in Equation (6), and Qm for Qi in Equation (5) (Table 2). Dm ranges from 3.7×104m3a1 at Margerie Glacier to 5.0 × 104m3a 1 at Grand Pacific Glacier where most of each terminus is exposed to tidewater. At Muir Glacier, the width of the ice cliff exposed to tidewater diminished from 900 m in 1978 (Reference Brown, Meter and PostBrown and others, 1982) to a perimeter of 475 m around a calving embayment in 1991. The debris released by ice-cliff melting is estimated at 5.8 × 104m3a 1 in the early 1990s.
The morainal banks of Grand Pacific, Margerie and Muir Glaciers receive a net addition of 8.4 × 104, 7.8 × 104 and 7.2 × 104m3a−1, respectively, through calve dumping and ice-cliff melt-out. Local recycling of morainal-bank sediment near the grounding line delivers up to 6.2 × 105, 1.0 × 104 and 7.2 × 105m3a−1 to the grounding lines of these glaciers (Table 2). Iceberg-rafting (D i) accounts for the removal of 1.0 × 106m3a −1 (Grand Pacific), 9.3 × 105 m3 a −1 (Margerie) and 3.2 × 105 m3a −1 (Muir), or about 25–96% of the debris fluxes.
The flux of sediment transported in a deformable bed has been roughly estimated assuming a 60 cm thick deforming layer (e.g. Reference Humphrey, Kamb, Fahnestock and EngelhardtHumphrey and others, 1993) and a linear velocity profile (Reference AlleyAlley, 1991a). Assuming plug-How conditions, the average velocity of the deforming layer is about half of the surface velocity, or about 262, 405 and 850 ma−1 for the Ferris tributary of Grand Pacific Glacier, Margerie and Muir Glaciers, respectively. We use the Ferris tributary here because it is the source of nearly 100% of the ice discharged into sea water by Grand Pacific Glacier (Reference HunterHunter, 1994). Based on these assumptions, we estimate that 1.3 × 105, 2.3 × 105 and 2.2 × 105 m3 a1could be transported, respectively, to the grounding lines of Grand Pacific, Margerie and Muir Glaciers by deforming layers. However, if soft-bed deformation is more localized, the subglacial sediment flux could be considerably less.
In comparison, the glacifluvial system supplies 7.7 × 107, 2.6 × 107, and 6.6 × 106m3 a−1 of sediment to the grounding line at Grand Pacific, Margerie and Muir Glaciers, respectively (Reference HunterHunter, 1994). The glacifluvial flux in Glacier Bay is, therefore, one to two orders of magnitude greater than the total glacial debris flux. The significant contribution of fluvial sediment to the marine environment in Glacier Bay reflects the region’s maritime climate and heavy precipitation, snow and ice melt and temperate ice conditions, each of which produce prodigious amounts of water for the glacifluvial system. These factors, along with friable bedrock, produce conditions favoring high sediment yields.
In other regions of the world. Reference Andrews, Milliman, Jennings, Rynes and DwyerAndrews and others (1994) and Reference ElverhøiElverhøi and others (1995) suggest that in Greenland and Svalbard the glacifluvial input is secondary, as the relative contribution of the debris flux into the total sediment budget is greater at latitudes higher than Glacier Bay (>59°N). Clearly, the role of glacifluvial transport varies as the rate of ablation increases (Reference Hallet, Hunter and BogenHallet and others, in press). Thus, both the glacifluvial sediment flux and the debris flux must be considered in regional comparisons.
Conclusions
In this paper, we have described the flux of debris transported in and on the surfaces of the Grand Pacific, Margerie and Muir Glaciers. We have isolated the debris flux from the flux of sediment discharged by outwash streams into the fjord environment to provide constraints on debris transport. The debris flux is governed by several attributes including glacier size, number of tributaries, glacier flow speed, the bedrock geology, substrate and character, thermal regime, local relief and accumulation patterns. Muir Glacier’s diffuse-facies debris concentration ( x =0.6% by volume) is more than four times greater than those of Grand Pacific and Margerie Glaciers (∼ 0.1% by volume) and reflects the importance of lateral inputs (e.g. rockfalls, slides, avalanches, eolian transport) in narrow valleys. The proportion of Muir Glacier’s debris transported in englacial ice is 52%, with the remaining 48% being transported in basal and morainal ice facies. In contrast. 18% of Grand Pacific Glacier’s and 19% of Margerie Glacier’s debris flux is transported by englacial ice. Debris transport in Grand Pacific Glacier is about equally divided between morainal (39%) and basal ice facies (43%), indicating the importance of medial moraines in debris transport for glaciers with multiple tributaries. Basal ice facies transport nearly 74% of the total debris llux at Margerie Glacier where englacial ice is relatively clean and there are few medial moraines (Figs 2b and 8c).
Data presented in this paper represent a first attempt to quantify the volume of debris transported by ice to tidewater termini as part of the total sediment flux of glacial systems. In Glacier Bay, glaciers with tidewater termini lose as much as 96% of their debris to iceberg-rafting. At Grand Pacific and Margerie Glaciers, less than 10% of debris is released at the terminus and deposited near the grounding line. This volume is relatively insignificant when compared to the volume of glacifluvial sediment discharged from outwash conduits. Glacifluvial sediment production is as much as two orders of magnitude greater than the debris flux and constitutes 83.6–98.6% of total sediment yields. This relationship, however, may vary in cooler, drier subpolar and polar climates. In addition, our estimates of deforming-bed transport indicate that it may deliver 2–8 times as much sediment to the grounding line as is released by melting of debris-rich ice and dumped during calving events. By developing a better understanding of debris flux and the mechanisms that control grounding-line deposition, our data should be useful in constraining models of glacier sensitivity to sediment dynamics (Reference AlleyAlley. 1991b) and evaluating process variations under different climatic regimes.
Acknowledgements
The authors thank B. Hallet, J. Strasser, J. Hardy and two anonymous reviewers for critically reviewing previous versions of the manuscript. Their comments were insightful and informative. Thanks are also forwarded to the National Park Service at Glacier Bay National Park and Preserve for logistical support, and especially to J. Luthy, captain of the M/V Nunatak. Special thanks are given to L. Hunter for her assistance in the field, drafting and editing the numerous versions of the manuscript. Field assistance was also provided by J. Cai. The lead author thanks D. Hirst and E. Josherger for assistance and access to the aerial photograph archives at the U.S. Geological Survey Ice and Climate Project, Tacoma, Washington. Funding was provided by the U.S. National Science Foundation under grant DPP 88-22098 to R.D.P.; the Geological Society of America, Sigma Xi and the Department of Geology at Northern Illinois University to L.E.H.; and U.S. Army CRREL for D.E.L. under the Civil Works Cold Regions Water Resources Program, “Predicting Runoff and Sediment Yield from Partly Glacierized Basins”.
