Abstract. The recent influx of remote sensing data provides new opportunities for quantifying spatiotemporal variations in glacier surface velocity and elevation fields. Here, we introduce a flexible time series reconstruction and decomposition technique for forming continuous, time-dependent surface velocity and elevation fields from discontinuous data and partitioning these time series into short- and long-term variations. The time series reconstruction consists of a sparsity-regularized

more »

... squares regression for modeling time series as a linear combination of generic basis functions of multiple temporal scales, allowing us to capture complex variations in the data using simple functions. We apply this method to the multitemporal evolution of Sermeq Kujalleq (Jakobshavn Isbrae), Greenland. Using 555 ice velocity maps generated by the Greenland Ice Mapping Project and covering the period 2009–2019, we show that the amplification in seasonal velocity variations in 2012–2016 was coincident with a longer-term speedup initiating in 2012. Similarly, the reduction in post-2017 seasonal velocity variations was coincident with a longer-term slowdown initiating around 2017. To understand how these perturbations propagate through the glacier, we introduce an approach for quantifying the spatially varying and frequency-dependent phase velocities and attenuation length scales of the resulting traveling waves. We hypothesize that these traveling waves are predominantly kinematic waves based on their long periods, coincident changes in surface velocity and elevation, and connection with variations in the terminus position. This ability to quantify wave propagation enables an entirely new framework for studying glacier dynamics using remote sensing data.