Depth–density profile

Density of the core was calculated using the length, diameter and weight of each piece of core measured in the field prior to packing for storage and shipment. Figure 2 shows density as a function of depth; minor data gaps at shallow depth are related to core quality. The approximate pore close-off density of 830 kg m–₃ is reached at 42 m depth. This density is reached at 33 m in the UpB core from Whillans Ice Stream, and 52 m in the BC core (labeled UpBC in Fig. 1), sites with temperature and accumulation rate similar to UpC (Reference Alley and BentleyAlley and Bentley 1988).

Fig. 1. Location sketch map of the Siple Coast of West Antarcti showing general area of ice streams and location of cores taken from UpC, UpB, BC and Siple Dome.

Fig. 2. Density–depth data from UpC core. Data points at surface and 1m depth were calculated by averaging the measured ca densities of snow samples taken from a pit. Other points represent the measured density of the core determined in the field. Depth is taken as the top of a sample interval.

The shallowest data points in Figure 2 summarize the snow-pit densities plotted in Figure 3. The density at each point was determined from the weight of a measured volume of firn sampled from the snow-pit wall with snow-density boxes. Two boxes were used, 1 and 3 cm thick, with volumes of 33 and 99 cm₃, respectively, and firn samples were removed so as to encompass visibly homogeneous layers. The calculated densities of the samples were then averaged over the intervals 0–1 and 1–1.7 m for the top two data points in Figure 2. Also shown in Figure 3 are the summer annual layers picked by a combination of observation of the snow wall in the pit and measured firn densities (see discussion below). Note that while the summers of 1994 and 1992 do not show low-density spikes, observations in the pit show these to be summer layers.

Fig. 3. Density–depth relationships in 1.7 m snow pit. Density is calculated from weight of firn sample taken with snow-density box (see text). Years indicate summer layers based visual observation in the snow pit and measured density.surface represents the summer of 1996.

Annual layers

Annual layers in the core were counted during visual examination as the core was processed at NICL. Identification of annual layers was made on a freshly cut and planed surface of the core, bisected lengthwise. Examination of the core on a light table showed alternating layers of coarse-grained, low-density firn or ice that represent summer snow, and fine-grained, more homogeneous firn or ice that represents densely packed snow laid down in winter (Reference AlleyAlley, 1988). Figure 4 shows the depth–age relationships to 29.7 m depth, corresponding to an age of 168 years (the year 1828). Estimated accuracy is ±5%, based on ongoing studies at Siple Dome, which has temperature, accumulation rate and visible strata similar to those of Ice Stream C (personal communication from K. Taylor, 2001). Beginning at 29.7 m depth, an approximately 3 m thick layer was encountered that appears to be a break in the annuallayers. This layer is rather homogeneous, contains abundant small bubbles and, as described below, maybe a crevasse fill. Below this feature, sections of ice again display countable annual layers, but the count continuity is lost, and layers become less distinct with increasing depth.

Fig. 4. Age-depth relationship ofUpC core obtained by summing annual layers as determined by visual stratigraphy to a depth where annual layers were no longer discernible. Solid lines indicate melt layers observed in the core during visual logging.Deeper melt layers (to 50 m) in ice of uncertain age are not shown.

Electrical conductivity measurement and dielectric property measurement

During processing of the core at NICL, electrical conductivity measurement (ECM) (Reference Spencer, Alley and CreytsTaylor and others, 1992) and dielectric property measurement (DEP) (Reference Moore and ParenMoore and Paren, 1987) logs were taken. The ECM record shows one large, -100 mm thick spike in conductivity, indicating high acidity at 81 m depth. The DEP, which measures total ionic concentration in the core, also shows this spike at 81 m, and another at 63 m that has no equivalent in the ECM record.

Core characterization

Layers above the probable crevasse fill (29.7–33 m) are nearly horizontal, whereas deeper layers are inclined, as shown in Figure 5a. Dips increase to a maximum of about 45° at 40 m depth, and then generally decrease toward the bottom of the core.

Fig. 5. (a) Angle of bed dip (degrees from horizontal) vs depth; (b) number of fractures per meter of core vs depth.

Features that we interpret as the traces of healed fractures occur deeper than about 50 m (Fig. 6). A fracture trace is expressed as a plane of tiny bubbles bounded by a nearly bubble-free zone typically a few millimeters across. Occasionally the central plane of the original fracture is marked by larger disk-shaped bubbles aligned parallel to the fracture. These fractures typically are nearly vertical, and as shown in Figure 6, there is a prominent and pervasive bubble elongation normal to the fracture traces and probably linked to along-flow stretching (Reference Alley and FitzpatrickAlley and Fitzpatrick, 1999). Some fractures appear to have split or joined during propagation, and lozenge-shaped pieces of ice can be isolated between fractures. These pieces are roughly 5–50mm wide and 50–500mm long. In some cases, later fractures are observed to cross-cut earlier ones, or are inferred to have exploited pre-existing weakness by reopening former fractures to produce zones with more than one “central”plane of bubbles.

Fig. 6. Photograph of a section from the UpC core at 104.32 m depth looking down the axis of the core (sample for thin section was removed after core was cut in half for archiving). Note orientation of the elongated bubbles and fracture traces. Because of the thickness of the section and the slight inclination of the feature, some subsurface bubbles appear in the bubble-free zone around the fracture trace described in the text. The photo was taken by J. J. Fitzpatrick.

The full three-dimensional nature of these fractures cannot be characterized in a single, 100 mm diameter core, but no strong depth trends were observed in dips or waviness of the fracture traces. Thus, simply counting the number of separate fracture traces in a core segment provides a useful estimate of fracture density. As shown in Figure 5b, these fracture features were not observed above 50m or below 120 m, and exhibit abroad peak in occurrence near 90 m depth.

Visible strata

The layering is apparently annual and horizontal to 29.7 m depth, approximately the year 1828 (±8 years estimated error). This corresponds well with the timing of the shutdown of Ice Stream C based on previous studies (Reference Retzlaff and BentleyRetzlaff and Bentley, 1993). The ECM and DEP peaks at 81 m probably represent volcanic fallout; however, no ash was evident in visible inspection. Based on prominence and approximate age, this 81 m peak is possibly the fallout from the eruption of El Chichon, Mexico, in 1259 (Reference Palais, Germani and ZielinskiPalais and others, 1992). The origin of the DEP peak at 63 m is unknown and may simply represent the transportation inland of sea salts by a strong coastal storm.

The 3 m zone of relatively homogeneous-appearing ice just below 29.7 m depth that shows a lack of annual layers and numerous small bubbles is consistent with rapid infilling of a crevasse by fine-grained drift snow. Finer grains, hence smaller bubbles (Reference AlleyGow, 1968), likely reflect both the shorter time that passed for burial to given depth for crevasse fill compared to normal firn, andthe shorter time that the crevasse fill spent in the very shallow zone of strong temperature gradients where grain growth is greatly accelerated (Reference GowGow, 1975). We cannot tell whether this fine-bubbled, homogeneous zone represents complete infilling of the crevasse (which then would have been approximately 6–7m deep, allowing for twofold compaction during firn densification), or whether the initial crevasse extended deeper and the drill has exited through the side of the crevasse. Non-zero dips are observed in beds at the bottom of the crevasse fill, but are not especially steep (Fig. 5a).

Density profile

The depth–density profiles of the UpB core on active Whillans Ice Stream and of the adjacent, slow-moving ridge BC (Fig.7) were indistinguishable from each other for densities less than ∼550kg m^–3, but densification proceeded faster on the ice stream at greater depths (Reference Alley and BentleyAlley and Bentley, 1988). One interpretation of this is that ice-flow stresses in the ice stream accelerate densification only at densities higher than ∼550kgm^–3. Previous work (Reference AlleyAlley, 1987) has suggested that densification at low densities is dominated by grain-boundary sliding, which would be a linear-viscous process unaffected by ice-flow stresses. At higher densities, power-law-creep densification is believed to be dominant and would be affected by ice-flow stresses (Reference DuvalDuval, 1985; Reference Ebinuma and MaenoEbinuma and Maeno, 1985; Reference AlleyAlley, 1987; Reference Alley and BentleyAlley and Bentley, 1988). Our data are fully consistent with this view, presuming that the existence of the few meters of probable crevasse fill has not greatly affected the densification profiles. The shallow part of the UpC profile is similar to that from ridge BC, except for the shallowest part of the core where there were sampling problems. Starting at a point about 160 years before coring, at a core depth of 28.6 m, the profile at UpC exhibits a density increase with depth that becomes more similar to ice-stream values such as those seen from UpB. This change appears to occur over a period of 7–14 years, as represented by 1–2 m of core. Inspection shows (and has been validated using a firn-densification model of Reference E.M., Iliescu and RenshawSpencer and others (2001)) that one can produce a profile similar to that observed at UpC by starting with the steady, active-ice-stream UpB profile, and then burying it and densifying it along a“normal”ridge-type BC profile, beginning at the time of Ice Stream C shut-down. The non-steady curve in such an experiment would overlie the UpB and BC profiles at low density (5550kgm^–3) at all times, and would evolve from the UpB to the BC profile over time at densities between 550 kgm^–3 and that of ice, giving an intermediate profile today as observed. We thus suggest that UpC had a depth–density profile similar to that of UpB prior to shut-down, and is now developing a profile similar to that of BC. Furthermore, the extra stresses on ice streams compared to ridges appear to promote rapid densification at densities above about 550 kgm^–3 but not below, suggesting a switch from linear-viscous to higher-power creep at about that density.

Fig. 7. Depth–density relationships of UpC core (points) overlain by curves from Reference Alley and BentleyAlley and Bentley (1988) for BC and UpB.

Crevassing

We interpret the fracture traces observed below 50 m as incipient crevasses that did not fill with snow before they closed and healed. The numerous fractures observed, exhibiting some cross-cutting relations and reactivation of previously formed fractures, show that opening and closing of fractures were repeated frequently. Alternate stretching and compression during “accordion-like” flow over sticky spots may have been involved (Reference Anandakrishnan, Alley, Jacobel, Conway, Alley and BindschadlerAnandakrishnan and others, 2001). If the crevasses opened to significant width at the surface (more than a few millimeters?) and to significant depth (more than perhaps 1m below the surface?), this would have enabled introduction of finer-grained drift snow to regions of coarser-grained firn or ice, and this fill now would be evidenced as a zone as wide as (or wider than, allowing for along-flow stretching) the original crevasse with anomalously numerous and small bubbles. No such zones are recognized among these fracture traces. The small bubbles along the central planes of the healed fractures are in a plane only a single bubble wide and are more consistent with mismatch of opposing sides during closure. It thus appears that the fractures either did not reach the surface, or did not open sufficiently to capture drift snow. Our data are consistent with the suggestion by Reference Nath and VaughanNath and Vaughan (2000) that crevassing initiates at depth (approximately 10 m?) and propagates upwards as well as downwards. We also support the suggestion by Reference Vaughan, Corr, Doake and WaddingtonWhillans and others (2001) that fractures in the ice stream often heal followed by re-fracturing.

The lack of healed fracture traces in the deepest part of the core (Fig.5b) may indicate that this ice had been buried below the depth of crevassing before entering the ice stream. Assuming that the fracturing occurred especially in the upper 10–20m, has been buried by the accumulation since the slow-down of the ice stream, and was transported by active flow at a few hundred meters per year before slowdown with accompanying longitudinal strain, the most abundant fractures would have formed a few hundred km up-glacier, and may correspond to intense strain at the onset of streaming flow.

Apparent splitting and rejoining of some fractures isolating lozenge-shaped regions of ice suggests at least a small component of shear along the fractures (cf. Reference Engelder, Fischer and GrossEngelder and others, 1993, section 8; Reference Retzlaff, Lord and BentleySchulson and others, 1999). The bubbles along the central planes of fractures indicate that a mismatch developed between the sides before healing, perhaps resulting from at least slight shear motion, or alternatively, vapor transported water from the exposed surfaces of the fracture, changing the shapes of the opposing faces before the fracture closed. The bubble-free, few-millimeter-thick zones flanking these bubbly planes may represent vapor-diffusion closure of bubbles, or else some sort of process zone of enhanced strain in high-stress regions around bubbles during fracture-tip propagation (e.g. Reference Engelder, Fischer and GrossEngelder and others, 1993, section 7).

Melt layers extend back through the Holocene in the Siple Dome core (personal communication from S.B.Das and R.B. Reference Alley, Bindschadler, Alley and BindschadlerAlley, 2000, http://igloo.gsfc.nasa.gov/wais/abstracts00/Das.htm). The absence of melt layers from the deeper and older parts of the UpC core is consistent with this ice having been transported from higher, colder regions farther inland.