We present a novel method of improving signal-to-noise ratio in radargrams. The method uses singular spectrum analysis (SSA) to separate each individual radar trace into orthogonal components. The components that explain most of the original trace variance contain mainly physically meaningful signal, while the components with little variance tend to be noise. Adding the largest-magnitude components together until the sum of components accounts for the variance above the noise level (typically 60–80%) of the original trace variance results in a much cleaner radargram with more easily seen internal features than in traditionally filtered data. The radargrams can be further enhanced by envelope-detecting the SSA-filtered data, as this measure of instantaneous energy minimizes the deleterious effects of innumerable phase changes at dielectric boundaries. Subsequent incoherent stacking results in far more structured radargrams than are achieved with traditionally processed radar data and amplitude stacking.

A minimal model of a tidewater glacier based solely on mass conservation is compared with two one-dimensional numerical flowline models, one with the calving rate proportional to water depth, and the other with the flotation criterion as a boundary condition at the glacier terminus. The models were run with two simplified bed geometries and two mass-balance formulations. The models simulate the full cycle of length variations and the equilibrium states for a tidewater glacier. This study shows that the branching of the equilibrium states depends significantly on the bed geometry. The similarity between the results of the three models indicates that if there is a submarine undulation at the terminus of a tidewater glacier, any model in which the frontal ice loss is related to the water depth yields qualitatively the same non-linear behaviour. For large glaciers extending into deep water, the flotation model causes unrealistic behaviour.

Precise measurements of basal melting have been made at a series of 14 sites lying within a few kilometres of the grounding line of the Ronne Ice Shelf, Antarctica, where the ice thickness ranges from 1570 to 1940 m. The study was conducted over the course of 1 year and included a detailed survey of the horizontal deformation, as well as phase-sensitive radar measurements of the vertical displacement of both internal reflecting horizons and the ice-shelf base. Results from the surface survey show that the long-term viscous strain rate is modulated at tidal frequencies by (probably) elastic strains of order 10−5 per metre of tidal elevation. The radar measurements show a similar modulation of the long-term thinning/thickening of the ice shelf, with thickness oscillations up to a few centimetres in range. The long-term trends in ice thickness determined at points moving with the ice-shelf flow are consistent with a steady-state thickness profile. Vertical strain rates within the ice shelf were determined from the relative motion of internal reflectors. At two sites the observations were sufficient to discern the effect of tidal bending about a neutral surface 60% of the way down the ice column. Coincident measurements of horizontal and vertical strain imply a Poisson’s ratio of 0.5, and this combined with the asymmetric bending gives rise to the observed oscillations in thickness. At a number of sites the longterm viscous strain rates were found to be a linear function of depth. For an ice shelf this is an unexpected result. It can be attributed to the presence of significant vertical shear stresses set up close to the grounding line where the ice is still adjusting to flotation. Additional vertical motion arising from firn compaction was observed within the upper layers of the ice shelf. The additional motion was consistent with the assumption that firn density is a function only of the time since burial by steady surface accumulation. With both spatial and temporal fluctuations in the vertical strain rate accurately quantified it was possible to estimate the vertical motion of the ice-shelf base in response. Differences between the calculated and observed motion of the basal reflector arise because of basal melting. Derived melt rates at the 14 sites ranged from −0.11 ±0.31 to 2.51 ±0.10 m a−1, with a mean of 0. 85 m a−1 and a standard deviation of 0.69 m a−1, and showed no signs of significant sub-annual temporal variability. There was no obvious global correlation with either ice thickness or distance from the grounding line, although melt rates tended to decrease downstream along each of the flowlines studied. Previous estimates of basal melting in this region have been obtained indirectly from an assumption that the ice shelf is locally in equilibrium and have included a broad range of values. Only those at the lower end of the published range are consistent with the directly measured melt rates reported here.

We present an overview and the first systematic results on the backscatter properties of snow samples with different grain properties. During a 2 year study we investigated the effect of the apparent characteristics of the samples on the intensity enhancement in the direction of exact backscatter (retroreflection), using a specific instrument designed for backscatter measurements. We observed that a sharp peak in intensity (the hot spot) occurs for most types of snow, which is not observable with traditional goniometers. The key factors in the peak properties are the temperature (related to the changes in grain structure) and grain shape and size. These results form the basis of a larger backscatter data archive which is being applied in the systematic study and exploitation of ‘hot spots’ in remote sensing of natural targets, as well as in the development of airborne laser intensity measurement.

Slab avalanche release requires fracture initiation and propagation in a weak snowpack layer. While field tests of weak-layer strength are useful for fracture initiation, the challenge remains to find a verified field test for fracture propagation. We introduce the two current versions of a field test for fracture propagation propensity, and report results of testing conducted in the Columbia Mountains of British Columbia, Canada, during the winter of 2005. By extending the column of a stability test approximately 3 m in the downslope direction, the test method allows for the development of a flexural wave in the slab, and thereby maintains the contribution of this wave and the associated weak-layer collapse to the fracture process. Fracture lengths collected on a day and location where the propagation propensity of the snowpack was locally high show a bimodal distribution, with approximately 50% of observed fractures similar to those collected in stable snowpacks, and approximately 50% with much longer fracture lengths.

Generalized numerical models of sub-debris ice ablation are preferable to empirical approaches for predicting runoff and glacier response to climate change, as empirical methods are site-specific and strongly dependent upon the conditions prevailing during the measurement period. We present a modified surface energy-balance model to calculate melt beneath a surface debris layer from daily mean meteorological variables. Despite numerous simplifications, the model performs well and modelled melt rates give a good match to observed melt rates, suggesting that this model can produce reliable estimates of ablation rate beneath debris layers several decimetres thick. This is a useful improvement on previous models which are inappropriate for thick debris cover.

Controls on glacier calving rates are receiving increased scientific interest. At fresh-watercalving glaciers, limnological factors might be more important than glaciological ones. Measurements of thermo-erosional notch development at the calving ice cliff of Tasman Glacier, New Zealand, suggest that the calving rates at this glacier are directly controlled by the rate of thermal undercutting. Notch formation rates typically vary between 10 and 30 cm d–1 (maximum rate 65cmd–1) in summer, corresponding to an average calving rate of 34 m a–1. Notch formation is slower than waterline melt and is controlled by water temperatures and circulation, cliff geometry, debris supply and water-level fluctuations. The latter shift the position of undercutting, resetting the level of the notch formation process and thereby slowing it. The geometry of the notch and the debris supply determine the extent of influence of the lake on notch water temperatures and circulation. Hence, water temperatures in the lake are not necessarily indicative of the rate of notch formation. The prediction of rate of notch formation from far-field variables is hampered by the complex interaction of the influencing factors. The significance of thermal undercutting as a calving rate-controlling process decreases with increasing ice velocities, calving rates and surface gradients.