Abstract. In recent years a vast amount of glacier surface velocity data from satellite imagery has emerged based on correlation between repeat images. Thereby, much emphasis has been put on fast processing of large data volumes. The metadata of such measurements are often highly simplified when the measurement precision is lumped into a single number for the whole dataset, although the error budget of image matching is in reality not isotropic and constant over the whole velocity field. The spread of the correlation peak of individual image offset measurements is dependent on the image content and the non-uniform flow of the ice. Precise dispersion estimates for each individual velocity measurement can be important for inversion of, for instance, rheology, ice thickness and/or bedrock friction. Errors in the velocity data can propagate into derived results in a complex and exaggerating way, making the outcomes very sensitive to velocity noise and errors. Here, we present a computationally fast method to estimate the matching precision of individual displacement measurements from repeat imaging data, focussing on satellite data. The approach is based upon Gaussian fitting directly on the correlation peak and is formulated as a linear least squares estimation, making its implementation into current pipelines straightforward. The methodology is demonstrated for Sermeq Kujalleq, Greenland, a glacier with regions of strong shear flow and with clearly oriented crevasses, and Malaspina Glacier, Alaska. Directionality within an image seems to be dominant factor influencing the correlation dispersion. In our cases these are crevasses and moraine bands, while a relation to differential flow, such as shear, is less pronounced.
Glacier surface velocity measurements from radar interferometry and the principle of mass conservation
