ABSTRACT. The generalized (Glen) flow relation for ice, involving the second invariants of the stress deviator and strain-rate tensors, is only expected to hold for isotropic polycrystalline ice. Previous single-stress experiments have shown that for the steady-state flow, which develops at large strains, the tertiary strain rate is greater than the minimum (secondary creep) value by an enhancement factor which is larger for shear than compression. Previous experiments combining shear with compression normal to the shear plane have shown that enhancement of the tertiary octahedral strain rate increases monotonically from compression alone to shear alone. Additional experiments and analyses presented here were conducted to further investigate how the separate tertiary shear and compression strain-rate components are related in combined stress situations. It is found that tertiary compression rates are more strongly influenced by the addition of shear than is given by a Glen-type flow relation, whereas shear is less influenced by additional compression. A scalar function formulation of the flow relation is proposed, which fits the tertiary creep data well and is readily adapted to a generalized form that can be extended to other stress configurations and applied in ice mass modelling. BACKGROUNDIn natural ice masses the most important and common state of deformation is arguably a combination of approximately bed-parallel shear and vertical compression. For deformational flow with a stationary boundary, a region of simple shear is associated in an essential way with bulk transport of ice in glaciers, ice sheets and ice shelves, and this is generally accompanied by normal deformations associated with increasing velocities along the flow and divergence or convergence transverse to the flow.For a coordinate system with x and y horizontal and z vertical, and corresponding component velocities (u, v, w), simple shear deformation in the x direction can be characterized by du/dz = c where we note that the horizontal planes on which the forces generating shear deformation act do not rotate, while compression normal to these planes is described by dw/dz = k, where c/2 and k are the respective shear and vertical compressive strain rates. The compressive flow may be confined or unconfined, and quite generally the accompanying horizontal normal strain rates are du dx ¼ ð À 1Þk and dv dy ¼ Àk where the factors involving indicate the proportions of the deformations in the horizontal directions, relative to the rate of vertical compression. Note that = 1/2 corresponds to uniaxial compression in the z direction, while = 1 corresponds to longitudinally confined compression in the experiments reported here (Fig. 1).The generalized flow relation for ice involving the second invariants of the stress deviator and strain-rate tensors (Nye, 1953;Glen, 1958) provides a useful formulation for the interactions between the individual stress and strain-rate components for isotropic ice. This relation is not expected to apply for anisotropic i...

Laboratory creep deformation experiments have been conducted on initially isotropic laboratory-made samples of polycrystalline ice. Steady-state tertiary creep rates, , were determined at strains exceeding 10% in either uniaxial-compression or simple-shear experiments. Isotropic minimum strain rates, , determined at ˜1 % strain, provide a reference for comparing the relative magnitude of tertiary creep rates in shear and compression through the use of strain-rate enhancement factors, E, defined as the ratio of corresponding tertiary and isotropic minimum creep rates, i.e. . The magnitude of strain-rate enhancement in simple shear was found to exceed that in uniaxial compression by a constant factor of 2.3. Results of experiments conducted at octahedral shear stresses of to = 0.040.80 MPa indicate a creep power-law stress exponent of n = 3 for isotropic minimum creep rates and n = 3.5 for tertiary creep rates. The difference in stress exponents for minimum and tertiary creep regimes can be interpreted as a t0 stress-dependent level of strain-rate enhancement, i.e. .The implications of these results for deformation in complex multicomponent stress configurations and at stresses below those used in the current experiments are discussed.

A real-time video camera probe was deployed in a hot-water drilled borehole through the Amery Ice Shelf, East Antarctica, where a total ice thickness of 480 m included at least 200 m of basal marine ice. Down-looking and side-looking digital video footage showed a striking transition from white bubbly meteoric ice above to dark marine ice below, but the transition was neither microscopically sharp nor flat, indicating the uneven nature (at centimetre scale) of the ice-shelf base upstream where the marine ice first started to accrete. Marine ice features were imaged including platelet structures, cell inclusions, entrained particles, and the interface with sea water at the base. The cells are assumed to be entrained sea water, and were present throughout the lower 100-150 m of the marine ice column, becoming larger and more prevalent as the lower surface was approached until, near the base, they became channels large enough that the camera field of view could not contain them. Platelets in the marine ice at depth appeared to be as large as 1-2 cm in diameter. Particles were visible in the borehole meltwater; probably marine and mineral particles liberated by the drill, but their distribution varied with depth.

Abstract. Accurate high-resolution records of snow accumulation rates in Antarctica are crucial for estimating ice sheet mass balance and subsequent sea level change. Snowfall rates at Law Dome, East Antarctica, have been linked with regional atmospheric circulation to the mid-latitudes as well as regional Antarctic snowfall. Here, we extend the length of the Law Dome accumulation record from 750 years to 2035 years, using recent annual layer dating that extends to 22 BCE. Accumulation rates were calculated as the ratio of measured to modelled layer thicknesses, multiplied by the long-term mean accumulation rate. The modelled layer thicknesses were based on a power-law vertical strain rate profile fitted to observed annual layer thickness. The periods 380–442, 727–783 and 1970–2009 CE have above-average snow accumulation rates, while 663–704, 933–975 and 1429–1468 CE were below average, and decadal-scale snow accumulation anomalies were found to be relatively common (74 events in the 2035-year record). The calculated snow accumulation rates show good correlation with atmospheric reanalysis estimates, and significant spatial correlation over a wide expanse of East Antarctica, demonstrating that the Law Dome record captures larger-scale variability across a large region of East Antarctica well beyond the immediate vicinity of the Law Dome summit. Spectral analysis reveals periodicities in the snow accumulation record which may be related to El Niño–Southern Oscillation (ENSO) and Interdecadal Pacific Oscillation (IPO) frequencies.

The Antarctic continental shelf supports a high level of marine primary productivity and is a globally important carbon dioxide (CO2) sink through the photosynthetic fixation of CO2 via the biological pump. Sustaining such high productivity requires a large supply of the essential micronutrient iron (Fe); however, the pathways for Fe delivery to these zones vary spatially and temporally. Our study is the first to report a previously unquantified source of concentrated bioavailable Fe to Antarctic surface waters. We hypothesize that Fe derived from subglacial processes is delivered to euphotic waters through the accretion (Fe storage) and subsequent melting (Fe release) of a marine‐accreted layer of ice at the base of the Amery Ice Shelf (AIS). Using satellite‐derived Chlorophyll‐a data, we show that the soluble Fe supplied by the melting of the marine ice layer is an order of magnitude larger than the required Fe necessary to sustain the large annual phytoplankton bloom in Prydz Bay. Our finding of high concentrations of Fe in AIS marine ice and recent data on increasing rates of ice shelf basal melt in many of Antarctica's ice shelves should encourage further research into glacial and marine sediment transport beneath ice shelves and their sensitivity to current changes in basal melt. Currently, the distribution, volume, and Fe concentration of Antarctic marine ice is poorly constrained. This uncertainty, combined with variable forecasts of increased rates of ice shelf basal melt, limits our ability to predict future Fe supply to Antarctic coastal waters.

ABSTRACT. The northwestern sector of the Amery Ice Shelf, East Antarctica, has a layered structure, due to the presence of both meteoric ice and a marine ice layer resulting from sub-shelf freezing processes. Crystal orientation fabric and grain-size data are presented for ice cores obtained from two boreholes $70 km apart on approximately the same flowline. Multiple-maxima crystal orientation fabrics and large mean grain sizes in the meteoric ice are indicative of stress relaxation and subsequent grain growth in ice that has flowed into the Amery Ice Shelf. Strongly anisotropic single-maximum crystal orientation fabrics and rectangular textures near the base of the $200 m thick marine ice layer suggest accretion occurs by the accumulation of frazil ice platelets. Crystal orientation fabrics in older marine ice exhibit vertical large circle girdle patterns, influenced by the complex stress configurations that exist towards the margins of the ice shelf. Post-accumulation grain growth and fabric development in the marine ice layer are restricted by a high concentration of brine and insoluble particulate inclusions. Differences in the meteoric and marine ice crystallography are indicative of the contrasting rheological properties of these layers, which must be considered in relation to large-scale ice-shelf dynamics.

Abstract. The microstructure of polycrystalline ice evolves under prolonged deformation, leading to anisotropic patterns of crystal orientations. The response of this material to applied stresses is not adequately described by the ice flow relation most commonly used in large-scale ice sheet models – the Glen flow relation. We present a preliminary assessment of the implementation in the Ice Sheet System Model (ISSM) of a computationally efficient, empirical, scalar, constitutive relation which addresses the influence of the dynamically steady-state flow-compatible induced anisotropic crystal orientation patterns that develop when ice is subjected to the same stress regime for a prolonged period – sometimes termed tertiary flow. We call this the ESTAR flow relation. The effect on ice flow dynamics is investigated by comparing idealised simulations using ESTAR and Glen flow relations, where we include in the latter an overall flow enhancement factor. For an idealised embayed ice shelf, the Glen flow relation overestimates velocities by up to 17 % when using an enhancement factor equivalent to the maximum value prescribed in the ESTAR relation. Importantly, no single Glen enhancement factor can accurately capture the spatial variations in flow across the ice shelf generated by the ESTAR flow relation. For flow line studies of idealised grounded flow over varying topography or variable basal friction – both scenarios dominated at depth by bed-parallel shear – the differences between simulated velocities using ESTAR and Glen flow relations depend on the value of the enhancement factor used to calibrate the Glen flow relation. These results demonstrate the importance of describing the deformation of anisotropic ice in a physically realistic manner, and have implications for simulations of ice sheet evolution used to reconstruct paleo-ice sheet extent and predict future ice sheet contributions to sea level.

ABSTRACT. We present a new method for extracting the direction of surface flow for ice sheets, based on the detection of flow-induced features that are visible in satellite imagery. The orientation of linear features is determined using a Radon transform and only requires a single image. The technique is demonstrated by applying it to the RADARSAT mosaic of Antarctica, over the Lambert Glacier-Amery Ice Shelf region of East Antarctica. Comparisons with both existing flow-direction fields and traced streamlines over the same area provide an evaluation of the method. We also illustrate its application to Landsat 7 imagery.