No abstract

ABSTRACT. Polycrystalline ice near an ice divide typically shows a crystal fabric (crystal preferred orientation) with c axes clustered vertically. We explore the effect of this fabric on the large-scale flow pattern near an ice divide. We incorporate an analytical formulation for anisotropy into a non-linear flow law within a finite-element ice-sheet flow model. With four different depth profiles of crystal fabric, we find that the effect of fabric is significant only when a profile has a minimum cone angle of less than $ $258. For a steady-state divide, the shape and size of the isochrone arch can depend as much on the crystal fabric as it does on the non-linearity of ice flow. A vertically oriented fabric tends to increase the size of the isochrone arch, never to reduce it. Also, non-random fabric has little effect on the ice-divide-flow pattern when ice is modeled as a linear (Newtonian) fluid. Finally, when we use a crystal-fabric profile that closely approximates the measured profile for Siple Dome, West Antarctica, the model predicts concentrated bed-parallel shearing 300 m above the bed.

Geophysical data are used with ice flow models and generalized accumulation histories to estimate age and annual layer thickness versus depth for two anticipated ice core sites in West Antarctica: Siple Dome (81.65°S, 148.81°W) and an inland site on the West Antarctic Ice Sheet (WAIS). This modeling experiment predicts that 104 year‐old ice is at ∼50% depth and 105 year‐old ice is at ∼90% depth at both sites. Both of these cores could contain climate information through the last glacial cycle with annual resolution through the Holocene. The predicted similarity in resolution and record length between the two cores suggests that they could be compared in detail to obtain both spatial and temporal information about the paleoclimate and history of the West Antarctic ice sheet.

If an ice sheet is frozen to its bed, deep ice directly under a divide experiences low deviatoric stress and is relatively hard, because the rheology of polar ice is described by a power-law constitutive relation. In steady state, stratigraphic layers tend to form an arch (“Raymond bump”) in this region. However, when the basal ice can slide, the stresses are redistributed, and longitudinal extension due to sliding is associated with increased deviatoric stress in the deep ice under the divide. This increased deviatoric stress weakens the tendency to form a Raymond bump. To find a realistic spatial distribution of sliding under an ice divide, we incorporate a thin layer of viscous till in a finite-element plane-strain flow model. The resulting basal “sliding” velocity varies approximately linearly with distance from the ice divide. By varying the till viscosity, we can adjust the amount of basal motion. We find that the Raymond bump decays exponentially with the fraction of total ice flux carried by sliding: the arch is 50% smaller when the sliding flux is only 7% of the total ice flux. This implies that the possibility of a wet bed must be considered when inferring past ice-divide locations from radar internal layering.