Understanding the controls on seasonal velocity change for tidewater glaciers may provide insight into long-term retreat and acceleration. Leveraging recent high-resolution satellite data, we examine changes in surface elevation, velocity, and terminus position for four glaciers in Central Western Greenland over 2015–2021. Our approach uses a simplified force balance focused at the terminus to model the expected response in upstream velocity caused by the observed terminus changes. We find that seasonal velocities are strongly controlled by terminus advance/retreat for two glaciers. Residuals between modeled and observed velocities reveal two distinct signals: summertime pulses coincident with peak runoff and wintertime speedup extending several kilometers inland of the terminus. We evaluate the sensitivity of terminus-driven velocity to elevation change by incorporating observed seasonally varying surface topography and applying controlled modifications to the profile, specifically uniform vertical shifts and variations in surface slope. We find surface slope changes impact velocity response to terminus changes more than spatially uniform changes in elevation. Increased surface slope amplifies velocity response to terminus changes. While simplified, our model could be applied to other glaciers to assess the importance of terminus position change as a driver of seasonal velocity.
NASA HeadquartersNNH20ZDA001N-ICESAT2Jackson School of Geosciences, University of Texas at AustinPostdoctoral FellowshipIntroduction
The Greenland Ice Sheet (GrIS) is currently the largest land ice contributor to present-day rising sea level with an acceleration in mass loss over the past few decades in part attributed to ice discharge through outlet glaciers . This acceleration underscores the importance of comprehending the intricate mechanisms that govern glacier dynamics. Despite the prevalence of glacier acceleration in Greenland, there exists notable spatio-temporal variability in glacier velocity change at a range of time scales . This is likely because glacier velocity is influenced by local conditions, like topography, and regionally by environmental factors. For example, runoff and ocean thermal forcing can influence velocity by changing basal friction , subaqueous melt rates , and through terminus fluctuations .
Outlet glaciers in Greenland present seasonal velocity changes , but it is not clear what exact processes are involved. In part, this is because many of the processes that influence outlet glacier velocity are synchronized, making it difficult to disentangle cause and effect. Some studies have suggested increased air temperature leads to summertime melting and terminus retreat, which reduces ice contact with the bed and/or fjord walls and increases net force at the calving cliff, both driving glacier acceleration . Numerical simulations also suggest a non-linear feedback between the rate of terminus change and ice discharge . Others have indicated that glacier acceleration is due to seasonal changes in basal lubrication related to changes in the subglacial hydrological system , which can cause complex responses in glacier velocities as subglacial conduits grow more efficient .
Understanding seasonal drivers of glacier change may help researchers to understand longer term glacier changes. Moreover, numerical simulations suggest that seasonal changes can induce systematic bias in mass loss estimates at the multi-decadal time scale, compared to simulations that do not include seasonal forcing, because the magnitude of mass loss is higher during retreat than mass gain for an equal amount of advance . Thus, investigating the patterns of glacier seasonality may be critical to unveiling the factors that control glacier dynamics into the future.
Previous efforts have classified seasonal velocity variations into types based on the observed timing of velocity change and their correlation to runoff and terminus change . This classification enhances our understanding of what drives these seasonal variations. Three categories have been generally described in : type-1 has positive correlation of glacier velocity to glacier terminus retreat; type-2 has positive correlation of glacier velocity to summer runoff, and; type-3 shows glacier velocity that slows in late summer and speeds up in winter. As noted in these studies, these velocity patterns connect to multiple underlying processes, and the types are used to describe patterns rather than to assign a single process. In a subsequent study, used high temporal resolution velocity data, terminus changes, and radar-backscatter intensity values to classify 45 glaciers into the same three types. extended their analysis to 221 glaciers across Greenland, focusing on type-2 and type-3 behaviors. combined velocity and runoff data using a machine-learning approach to classify glaciers throughout Greenland and discussed how subglacial drainage development, runoff volume, ocean interaction, and local conditions can produce similar seasonal patterns. They also noted that multiple processes (e.g., terminus position and runoff) jointly control seasonal flow for many glaciers. applied principal component analysis to decompose glacier seasonality in Sermilik Fjord from a statistical perspective, quantifying the prevalence of multiple seasonality types at a single glacier and emphasizing the importance of spatial patterns across the glacier.
These studies have made clear that seasonal velocity patterns are driven by multiple processes. For instance, reported that some glaciers exhibited multiple velocity types within a single year, and explicitly quantified the coexistence of multiple seasonal patterns at individual glaciers. However, these studies were primarily designed to describe dominant or statistically separable patterns. Due in part to the temporal resolution of the data (e.g., monthly or bi-monthly), it has been challenging to resolve the sequential occurrence of multiple processes within a single season. Consequently, the compounding effects of processes that act at sub-seasonal timescales have received less direct attention, especially in physics-based modeling.
Building on the foundation of previous studies, we re-examine glacier seasonality using high-frequency terminus , velocity , and surface elevation change observations for four glaciers in Greenland. We interpret these observations with an analytical model of velocity response to terminus position change . While simple, this model provides a physics-based method to help isolate the myriad processes that influence ice velocity by removing one potential driver, terminus change. By isolating the velocity component driven by terminus changes, we find that glaciers more typically exhibit multiple velocity patterns simultaneously or sequentially within a year. Using high resolution (subweekly) data in this study, we can more accurately decompose these patterns throughout a season, providing a sub-seasonal characterization that may be more appropriate for glacier velocity classification.
Study region
We investigate four glaciers in central-west Greenland (Fig. ): Sermeq Avannarleq (AVA), Sermeq Kujalleq (KUJ), Kangilernata Sermia (KAN), and Eqip Sermia (EQP) over the time period 2015–2021. This time span is specifically chosen to take advantage of the increased sample frequency available in both velocity and terminus position data due to the launch of Sentinel-1/2 in April 2014, which offers full year coverage of velocity data. We choose these 4 glaciers because we want to understand how our simple model might aid interpretation of seasonal velocity change for glaciers that have been extensively studied in the past. These glaciers are selected because they (1) are grounded; (2) exhibit regular seasonal changes in both terminus position and velocity with minimal long-term variations over our study period and; (3) exhibit a range of sub-seasonal behavior in both the terminus and velocity variability. For example, all the four glaciers advance in winter and retreat in summer and yet their seasonal velocity behavior differs over time and space. EQP and KUJ speed up during summertime terminus retreat while AVA and KAN slow down during summertime terminus retreat . In addition, AVA, KUJ, KAN, and EQP are located close to one another, suggesting that they likely experience the same regional climate forcing. All flowline samples for these glaciers are located within the ablation zone , implying that meltwater influence is expected across all sampled sites. They also exhibit distinct local effects, such as geometry that impact their dynamic behaviors. For instance, EQP's steep slope results in higher glacier velocities and more frequent calving events , whereas the shallow sills of AVA and KUJ contribute to their remarkable stability over the past century .
These glaciers have been studied extensively by previous authors. All glaciers have been relatively stable (they have not retreated more than the amplitude of seasonal fluctuations) since 2015 (Table ). All four glaciers have experienced terminus retreat over the more extensive satellite era; AVA 0.64 km, KUJ 1.80 km, KAN 2.65 km, and EQP 2.54 km from 1999 to 2016 . Studies of calving style for these glaciers indicate that they all primarily calve via serac failure due to undercutting related to terminus melt , although KUJ has been reported to have a mix of serac failure and buoyant flexure calving styles . This suggests that AVA, KAN, and EQP are grounded and that KUJ may approach floatation at times or in smaller regions of the terminus. This is consistent with an earlier study finding that these four glaciers experience terminus change that is predominantly driven by changes in subglacial runoff . Additional details on each glacier can be found in Table .
Glacier and data locations. Glaciers examined in this study include Sermeq Avannarleq (AVA), Sermeq Kujalleq (KUJ), Kangilernata Sermia (KAN), and Eqip Sermia (EQP). The terminus traces in 2018 from AutoTerm are colored by day of year. The red points show the locations where ITS_LIVE velocities were extracted for the model comparison in the results section. Blue points are the locations where we tested the sensitivity of the location of the velocity observation to the results. Blue and red points mark the locations for velocities shown in Figs. S2–S7 in the Supplement.
Summary of statistics for observations used, glacier geometry, and simulation results. Glacier width changes negligibly with terminus position, so fixed values are used. Frontal depth is averaged over the area where the terminus varies. Distance from the sampling location to the terminus is averaged across all terminus positions. Seasonal changes and model misfits are both averaged over the study period. Note that “n/a” stands for “not applicable”.
Glacier NameKUJEQPKANAVAGlacier Width(km)5346Glacier Frontal Depth(m)421134403399Average velocity over the study period(m yr−1)3699 ± 5791560 ± 1881177 ± 1701906 ± 332Terminus observations per month14111314Velocity observations per month18191919Velocity Uncertainty(m yr−1)12410910899Terminus Uncertainty(m)23161317Model Velocity Uncertainty(m yr−1)6250459Distance of the Glacier Terminus to the Sampling Location(km)1.5 ± 0.42 ± 0.23.4 ± 0.41.7 ± 0.1Glacier retreat rate over the study period(m yr−1)10241110Seasonal Terminus Change(m)400217370138Seasonal Velocity Change(m yr−1)853359404862Seasonal Simulated Velocity Change(m yr−1)80933126065Upper Limit of Stress Coupling(km)2010n/an/aAverage Misfit between Simulated and Observed Velocity5.0 %6.20 %14.2 %63.1 %Data
We use velocity time series data generated using auto-RIFT and provided by the NASA MEaSUREs ITS_LIVE project . ITS_LIVE combines velocity products derived from Landsat-8, Sentinel-1, and Sentinel-2 and produces velocity data with a spatial resolution of 240 m. The number of velocity observations per month ranges from 18 to 19 for the four glaciers (Table ). For each glacier, we use six flowlines across the glacier from to extract velocities at six points along each flowline (red and blue points in Fig. ). We then average these velocities across all flowlines at each cross section to produce mean (across-flow) velocity time series from downstream to upstream for each glacier. The chosen flowlines of AVA predominantly concentrate on the western side of the basin because AVA is formed by the confluence of two upstream tributaries, and we focus on the main tributaries with higher velocities on the western side. Blue points in Fig. are used to test the model sensitivity to distance from the terminus as described in the Supplement (Sect. S2). We applied a 40 d trailing rolling average to improve the visualization of the time series, as the original data were noisy. To compensate for the phase shift introduced by the rolling average, we shifted the time series back by 20 d. The rolling average was used only for visualization; all statistical analyses, including correlations and residuals, were conducted without it. To ensure the quality of the ITS_LIVE velocity data, we filter out measurements with dt greater than 45 d, where dt represents the interval between two scenes used to obtain glacier velocities. In addition to ITS_LIVE velocity data, we also extract along-flow velocity profiles from monthly velocity mosaics provided by the Greenland Ice Sheet Mapping Project (GrIMP) with a spatial resolution of 30 m. These along-flow profiles are exclusively used to examine the spatial pattern of velocity changes over seasons (Fig. ), while all velocity time series observations are sourced from ITS_LIVE velocity data.
Time series for EQP. (a) Terminus position change time series from AutoTerm, with light yellow shading showing uncertainty. (b) Surface elevation from ICESat-2 data fusion. The dynamic elevation change is calculated by subtracting elevation change due to surface mass balance from total elevation change. (c) Observed and simulated velocity. Light blue and pink shading show the uncertainties of observations and simulations, respectively. (d) Velocity residual between observed and simulated velocity, along with the daily runoff from GSFC-FDMv1.2.1 (values are per-unit-area runoff, in m yr−1). In all panes, gray vertical bars mark melt seasons using the 1st percentile of each year’s runoff as the threshold. The vertical black lines mark the onset of annual glacier retreat. The orange dashed lines in (c) marks the two periods in Fig. a and b, respectively.
Terminus position data come from AutoTerm , a machine learning pipeline that automatically produces terminus traces with an average sampling frequency of 10 per month since 2014. We derive a time series of terminus change by calculating the sequential area changes between termini, accumulating these over time, and then normalizing this by a static glacier width of 6 km for AVA, 5 km for KUJ, 4 km for KAN, and 3 km for EQP. The area-based metrics account for both terminus shape and fjord geometry , and unlike the curvilinear box method, they do not require a predefined box to represent the fjord geometry. We apply the same 40 d rolling average to smooth terminus data as described above. The uncertainties for ITS_LIVE, AutoTerm data, and simulated velocity from the terminus-driven model are represented as shaded areas in our time series figures (Figs. to ) and we present the average uncertainties in Table .
Time series for KUJ. The figure follows the same design as Fig. . The orange dashed lines in (c) marks the two periods in Fig. c and d, respectively.
The full record of KAN. The figure follows the same design as Fig. . The orange dashed lines in (c) marks the two periods in Fig. e and f, respectively.
Time series for AVA. The figure follows the same design as Fig. . The orange dashed lines in (c) marks the two periods in Fig. g and h, respectively.
We generate surface elevation data through a novel fusion of ICESat-2 data with DigitalGlobe high-resolution digital elevation models (DEMs), termed “DG-IS2-DEM” producing four DEMs per year from Fall 2018. The algorithms that generate the DG-IS2-DEMs are described in the Supplementary Information (Sect. S1). We also use ArcticDEM as supplementary elevation data in locations where the DG-IS2-DEMs do not extend to the most advanced terminus position found in AutoTerm. To determine ice thickness, we subtract surface elevation data from bed elevation data from BedMachineV5 , which assimilates seafloor bathymetry and ice thickness data through a mass conservation approach . BedMachine is an estimate of the bed topography, not a measurement and carries a spatial resolution of 150 m and uncertainties of ∼ 60 m for AVA, ∼ 70 m for KUJ, ∼ 20 m for KAN, and ∼ 90 m for EQP. These low uncertainties result because each of these glaciers has radar-determined bed elevation data to within ∼ 1 km of the current terminus position . For both of these data sets we extract the elevation profiles along each flowline individually.
We use the runoff component in GSFC-FDMv1.2.1 , which provide the evolution of surface mass balance and its individual components from 1 January 1980 to 30 June 2022 at 5 d temporal resolution and 12.5 km spatial resolution on North Polar Stereographic Projection. We extract 5 d runoff time series at points in Fig. . The onset and cessation of runoff served as proxies for the melt season timing. Following , we define runoff values below the 1st percentile of each year as wintertime, and values above that threshold as summertime (Figs. –).
Terminus-Driven model
Force balance methods can be used to understand the dynamic evolution of glaciers by examining the balance of stresses on them . However, the force balance method requires a second-order derivative of surface velocity data, which can result in large uncertainties when using satellite products. Therefore, we adopt a different approach that employs a modified force balance termed the “terminus-driven model” described by . This model explicitly considers the influence of dynamic changes at the glacier terminus on upstream velocity. By using the terminus-driven model and near daily velocity data, we are able to isolate the contribution of terminus variations to the observed variations in the velocity time series. The terminus-driven model is a 1-D model along the flow direction, assuming ice mélange-free conditions, a constant glacier geometry over time, and a grounded glacier. The terminus-driven model first quantify the driving stress (τd) at the points where we simulate velocity, expressed as
τd=-ρigH∂h∂x
where ρi is the density of ice (910 kg m−3), g is the gravitational acceleration (9.81 m s−2), H is ice thickness, h is ice surface elevation, x is the distance along the flowline with the direction from downstream to upstream. The driving stress is only calculated at the points where we simulate velocity. Secondly, the terminus-driven model examines an additional force at the water-terminating calving face, determined by the height above the fjord surface at the calving front and the density of seawater. The difference between these two forces at the terminus is expressed as
F=12ρigH2-12ρwg(H-h)2
where ρw is the density of seawater (1028 kg m−3). Here, we call F the “frontal force” following the naming convention found in , although it has units of N m−1. The force balance at the terminus requires the frontal force to be balanced upstream by the longitudinal stress, which redistributes much of the frontal force to the glacier upstream. We term the longitudinal stress in the upstream region that originates from the frontal force as τF(X), which pulls/opposes the glacier and enhances/reduces the original driving stress (e.g., τd+τF(X)). X is the distance between terminus and the point where we simulate velocity.
The integration of τF(X) along the flowline equals F and is assumed to linearly decrease upstream of the terminus to zero at the upper limit of stress coupling following :
τF(X)=2Fλ1-Xλ
where λ is the upper limit of stress coupling. Terminus variations cause changes in the geometry at the free calving face, consequently influencing the frontal force (F). These changes in the frontal force, subsequently, lead to modifications in the enhanced driving stress (τd+τF(X)) in the upstream region. Thus, the calculation of F and τF(X) over time requires elevation data along the profile.
In the laminar flow model or the shallow ice approximation, when basal sliding is zero, the surface velocity has a linear relationship with the cube of driving stress (n=3):
V=12AHτd3
where A is a constant from Glen's flow law and V is surface velocity. Such a linear relationship also applies when the driving stress is mainly balanced by lateral drag :
V=12AτdH3W4
where W is half of the glacier width. Based on the above two models and following the method designed by , we assume a linear relationship between velocity and the cube of the enhanced driving stress: τd+τF(X). Finally, the predicted velocity from terminus changes is thus given by:
VX,tV0=τd+τF(X)τd+τF03
where V0 is a reference velocity at the same location as velocity observations and τF0 is the τF(X) corresponding to the reference velocity. Following , for each year, we use the minimum velocity observation as the reference velocity and calculate τF0 based on the terminus position nearest the date of the reference velocity. Using Eq. (), we simulate velocity observations at both the blue and red points in Fig. using corresponding geometry profiles along each flowline. Then, we average the simulated velocity across all the flowlines at each cross-section to produce mean velocity time series in a manner consistent with the observed velocity.
We vary λ for each glacier and choose the one that produces the lowest mean difference between observations and simulated velocity (Table S1). The mean difference is determined by:
∑Nn∑6mabs(Vm,nmodel-Vm,nobs)Vm,nobs÷(N×6)×100%
where 6 is the total number of velocity points on the same cross section and N is the total number of time steps in the entire time series. We use the red points in Fig. to find an optimal upper limit of stress coupling (λ). The optimized λ values are then used to simulate velocity time series at the blue points in Fig. , which allow us to test the sensitivity of the velocity position to the results.
We also calculate the velocity residual by subtracting the simulated velocity with fixed geometry from the observed velocity. The velocity residual refers to the portion of the velocity that is not associated to terminus changes but is affected by other factors, such as runoff. We then plot residuals overlapping with runoff for ease of comparison (Figs. –). We assign the uncertainty of the terminus change as the error of X in Eq. (), and propagate this error to obtain the uncertainty in the simulated velocity.The elevation uncertainties are not directly included, as our model assumes invariant geometry (following ).
Although the terminus-driven model assumes invariant geometry, we analytically examine the impact of seasonal variations in surface elevation on velocity simulations, offering a complementary sensitivity test for elevation-related effects. Specifically, we leverage the new time-varying DG-IS2-DEM and periodically update the elevation profiles each quarter from Fall 2018, maintaining profile fixed within each quarter. We produce a simulated velocity for all glaciers with and without time-varying surface elevation to evaluate the impact of seasonally-varying surface elevation change on velocity. For the fixed geometry simulations, we choose a time step from DG-IS2-DEM with an extent that aligns best with the position of the terminus when it is most advanced. This provides the most complete elevation profile across the terminus region. For EQP, KAN, and AVA we choose the October 2019 DG-IS2-DEM and for KUJ, we use the April 2019 DG-IS2-DEM time step. To further investigate the impact of changing elevation on the simulation results, we consider only KUJ as an example and conduct two experiments using artificially modified surface elevation; (1) we shift the entire elevation profile vertically by ± 10–20 m and; (2) we alter the surface slope by ± 2 % within the 2 km-frontal region.
Results
Our comparative analysis of model outputs, runoff data, and velocity observations reveals that the seasonal velocity variations of all four glaciers exhibit three distinct patterns, which align with those identified by . We find that each glacier exhibits two of these patterns, occurring either alternately or simultaneously each year. Following the velocity-pattern classification of , we describe each glacier’s seasonal behavior as a combination of patterns: EQP and KUJ exhibit patterns consistent with type-1 plus type-2, whereas KAN and AVA exhibit patterns consistent with type-2 plus type-3.
Comparison between velocity simulation and observations
We compare the simulated velocity time series with velocities from satellite observations to determine whether seasonal velocity variations are influenced primarily by terminus change, co-influenced by subglacial hydrological system, or entirely independent of terminus change. We find that the terminus-driven model, after tuning λ, produces a good match to velocity observations for KUJ and EQP (Figs. and ) but not for KAN and AVA (Figs. and ). Note that λ only scales the amplitude of the simulated velocity without affecting its phase or temporal pattern. This indicates that both EQP and KUJ have type-1 velocity behavior. This is supported by the coincident timing of the end of terminus retreat and the peak summertime velocity (vertical black lines in Figs. and ), even in instances when retreat continues beyond the end of the melt season (Fig. ). This is also supported by the percentage difference between the simulated and observed velocities over all years, which are 5.0 % for KUJ and 6.2 % for EQP (Table S1), with correlations of 0.73 and 0.67, respectively (Fig. ) and small residuals (Figs. and ). The λ chosen based on minimizing the difference between observed and modeled velocity are EQP = 10 km, KUJ = 20 km, KAN = 25 km, and AVA = 25 km (Table S1). However, KAN and AVA’s seasonal velocities are not primarily driven by terminus variation, so we do not consider the estimation of λ valid for these glaciers.
Scatter plot showing the comparison between simulated and observed velocity for KUJ, EQP, KAN, and AVA. The 1:1 line is as shown for each.
For AVA and KAN, we find that simulated velocities differ substantially from the observed velocities suggesting that the terminus-driven model is insufficient to explain the observed velocities and that velocity change is not primarily driven by terminus fluctuations. For AVA, we find a poor correlation between observed and simulated velocity (-0.36; Fig. ) and residuals are just slightly smaller than the velocity observations themselves. For KAN, the simulated and observed velocities are out of phase but of similar magnitude, producing small residuals (Fig. ) but a poor correlation, R=0.36 (Fig. ). Although our simulation results indicate no significant correlation between seasonal-scale glacier velocity changes and terminus variations at KAN, found a large, full-thickness calving event induced speed up in the frontal region of KAN on an hourly time-scale.
Non-terminus velocity processes
Although the terminus-driven model adequately resolves the majority of seasonal velocity variability for KUJ and EQP, the velocity residuals reveal additional sub-seasonal velocity changes that are not related to terminus change. One such signal is an additional pulse in velocity (acceleration and deceleration) that occurs near the middle of the melt season for all glaciers (Figs. –), meaning that all glaciers have type-2 behavior. For KUJ, these melt-season pulses are not evident in the original velocity time series, but only become apparent in the velocity residuals, which further validates the effectiveness of our method in decomposing the coupled seasonal velocity variations of glaciers. We find that these mid-summer residual speed ups account for more than 80% of the observed summertime velocity signal across all glaciers, with KAN and AVA exhibiting a higher percentage (91 % and 86 % respectively) than EQP and KUJ (84 % and 82 % respectively). We typically observe a lag time between the peak of the velocity residual followed by the peak of summer runoff (Figs. –). We conduct cross correlation analysis between the velocity residual and runoff to assess the time lag between the two and quantify the maximum correlation coefficient, thereby investigating the extent to which velocity is influenced by runoff. These lag times differ among glaciers, measuring 24 d for KUJ, 22 d for EQP, 47 d for KAN, and 51 d for AVA. The maximum correlation coefficients are 0.49 for KUJ, 0.73 for EQP, 0.67 for KAN, and 0.67 for AVA (Fig. ).
Scatter plot showing the comparison between velocity residual and runoff for KUJ, EQP, KAN, and AVA. The phases of velocity residual have been adjusted to achieve the maximum correlation coefficient.
For KAN and AVA, we observe accelerations during winter (during terminus advance); velocities then plateau before the onset of the melt season in the following year, and early melt season accelerations with the annual maximum velocity reached in the middle of the melt season (Figs. and ). The terminus-driven model does not capture wintertime acceleration because across all glaciers the terminus is advancing in winter. For KAN, the model predicts slight deceleration in winter (Fig. ). For AVA, the simulated velocity lacks significant seasonality because (1) the seasonal terminus variations are small (Fig. ), and (2) the frontal region has a very low surface slope (Fig. S6). Together, these two reasons lead to small change in frontal force F (Eq. 2), and thus little seasonality in simulated velocity. For reference, the averaged seasonal range in terminus position is 138 m for AVA, while EQP is 217 m, KAN is 370 m, and KUJ is 400 m.
Experiments with seasonally varying surface elevation
We investigate the influence of changing surface topography by comparing the velocity simulated using a fixed geometry against the velocity simulated using a seasonally varying surface elevation from 2018–2022. We find minimal differences between these results for all glaciers (black versus red lines in Figs. –). Using KUJ as an example, the results from two analytical experiments involving artificially modified elevation profiles suggest that terminus-driven velocities are relatively insensitive to spatially uniform along-flow changes in surface elevation. In contrast, they show a high sensitivity to changes in surface slope (Fig. ): a 2 % absolute change in slope leads to a 2 % change in average velocity and a 6 % change in peak velocity. This result is important for providing context for interpreting the results of this model and our assumptions.
Experiment results for KUJ using artificially modified surface elevations. Panels (a) and (c) correspond to each other, as do (b) and (d). The blue shaded area represents observation uncertainties. Vertical black lines indicate the onset of the annual glacier retreat, and gray shaded vertical bars denote the melt season. In panel (d), the surface slope is altered by ±2 % within the 2 km frontal region.
Discussion
Using observations with high temporal resolution and a terminus-driven model to simulate velocity variations resulting solely from terminus change, we investigate the controls on sub-seasonal velocity changes for four GrIS outlet glaciers and find that in all cases, glacier velocity responds to multiple compounding processes. The seasonal velocity changes of two glaciers (KUJ and EQP) show a strong correspondence to seasonal terminus variation (type-1). All study glaciers experience sub-seasonal peaks in velocity that occur in the middle of the melt season (type-2) and, AVA and KAN exhibit wintertime speedup that occurs when their termini are advancing (type-3) (Figs. and ). Our combined pattern descriptions offer a different perspective from , which classified EQP and KAN as type-3, identified AVA as type-2 in 2018 and type-3 in 2017 and 2019, and did not classify KUJ. This difference reflects our use of high-temporal-resolution observations together with a simplified physical model, which allows us to identify the presence of multiple coexisting patterns and link them to interacting processes. We emphasize that these different classifications are not contradictory; rather, they arise from different methodological focuses (pattern-based annual classification vs. process-informed sub-seasonal analysis) and together enrich our understanding of glacier seasonal dynamics.
We hypothesize that the peaks in velocity residuals observed for all glaciers during the middle of the melt season (Figs. –) are the result of meltwater-driven acceleration and subsequent evolution of the subglacial drainage system . Early in the melt season, the subglacial drainage system is not fully channelized , thus as meltwater availability begins to increase (marked by increasing runoff in early summer), subglacial water pressure increases, enhancing basal sliding by reducing friction between the ice and the bed . As the melt season progresses, the drainage system becomes more channelized, allowing for more efficient water drainage . Consequently, the available meltwater decreases, producing a reduction in glacier speed. This change in efficiency of the subglacial drainage system may explain the lag between peak velocity residual and peak runoff as has been shown by previous authors . The more readily the subglacial drainage system channelizes, the earlier the peak velocity residual occurs relative to peak runoff. Beyond the melt season, the impact of terminus retreat on seasonal velocities can become more pronounced. For example, EQP typically has a melt season that ends in October, but the terminus continues to retreat until December/January (the vertical black lines are usually behind the end of the gray bar in the same year in Fig. ). This produces a wintertime velocity peak that is coincident with the most retreated terminus of EQP and is distinct from the melt-season peak.
To further explore the velocity increases that occur in summer, we examine the along-flow variability in velocity to determine how far upstream velocity changes occur (Fig. ). Our hypothesis is that changes in basal lubrication will be widespread over the region where water enters the bed, but downstream of the equilibrium line . To confirm that runoff drives the summertime speed up, we compare the along-flow spatial pattern of velocity change that occurs in the summer melt season of 2017 for EQP where summer velocity peaks are strongest (Apr 2017–Sep 2017; Fig. a) and the subsequent time period when runoff is low and the velocity is primarily influenced by terminus changes (Oct 2017–Mar 2018; Fig. b). The idea is that runoff-driven expansion of the subglacial drainage system would produce velocity changes that extend far inland from the terminus, while terminus-driven increases in velocity would decay rapidly from the terminus . The rapid decline in speed with distance from the terminus is expected when a glacier is terminus-driven because of the reduction in the terminus force with distance from the terminus . Conversely, elevated inland velocities are typical for melt-driven acceleration because meltwater percolates throughout a large portion of the ablation zone , which extends about 70 km inland of the terminus for EQP .
Velocity profiles over time for EQP, KUJ, KAN, and AVA. The average velocity profile has been subtracted for a better display of changes over time. The original velocity profiles are shown in Fig. S7. The shaded areas indicate regions where we obtain velocity variations in the frontal (left) and upstream (right) sections. Left column (a, c, e, g) Melt-season profiles for EQP (a), KUJ (c), KAN (e), and AVA (g). Right column (b, d, f, h) Post-melt-season profiles for EQP (b), KUJ (d), KAN (f), and AVA (h), when velocity is primarily influenced by terminus changes.
We choose a distance of 10 km upstream from the terminus to represent the inland distance and compare the range in velocity in this location to the range in velocity found at a location that is 2 km from the terminus. A larger ratio of velocity range between the upstream and the frontal region, in conjunction with the correlation between velocity residuals and runoff data, indicate that the velocity is likely driven by runoff. For EQP, we find that in the summer (Apr–Sep) the upstream velocity range is 64 % of that observed at the terminus (Fig. a) while in winter (Oct–Mar) the upstream velocity is just 12 % of what is observed at the terminus (Fig. b). In addition to EQP, we also find that in summer (Apr–Sep) the upstream velocity ranges for both KAN and AVA are both 56 % of that observed at the terminus (Fig. e and g). Thus, we conclude that summer accelerations of EQP, KAN, and AVA are strongly caused by runoff, as they are more noticeable upstream compared to the velocity changes in EQP during the winter, which are driven by terminus changes. For KUJ, the ratio of velocity range between the upstream and frontal region is 13 % in summer (Apr–Sep) (Fig. c), slightly larger than the 10 % ratio observed in winter (Oct–Mar) (Fig. d). Based on this observation, along with the correlation between KUJ's velocity residual and runoff (Fig. a), we conclude that runoff also influences the velocity changes of the KUJ glacier during the summer, although its impact is not as pronounced as that on the other three glaciers.
While AVA and KAN experience mid-summer velocity pulses similar to EQP, they do not exhibit notable terminus-driven seasonal acceleration (Figs. and ) and instead accelerate in winter. This is suggested by the poor correlation between velocity changes and terminus changes for these glaciers. Although our simulation results indicate no significant correlation between seasonal-scale glacier velocity changes and terminus variations at KAN, found a large, full-thickness calving event induced speed up in the frontal region of KAN on an hourly time-scale. To examine wintertime acceleration, we investigate the along-flow pattern of velocity change similar to above (Fig. f and h). For both glaciers we find significant inland acceleration in winter. For KAN, the range in upstream (farther from the terminus) velocity during winter is 72 % of the range in frontal velocity and for AVA, the range in upstream velocity is 36 % of the range of frontal velocity, which is larger than terminus-driven upstream velocity range that was observed for EQP (Figs. f and h versus b). We hypothesize that the elevated range of KAN and AVA's upstream velocities in winter suggests that winter acceleration is due to enhanced extensive basal slip.
For many parts of Greenland, basal water is generated year-round because of friction at the ice-bed interface and geothermal heat . During the onset of winter, refreezing of available basal water and viscous deformation over subglacial conduits can obstruct the drainage system. Consequently, water becomes trapped within an inefficient drainage network, leading to increased water pressure and winter acceleration . The winter acceleration phase ends when the melt season begins to supply significant volumes of additional water to the subglacial system, which increases basal water pressures causing summertime acceleration .
Our examination of how our model performs further upstream of the terminus (Supplementary Information) informs the transmission distances for terminus perturbations. and quantified the Péclet number (Pe) for the glaciers in our study. Most of the dynamic mass loss occurs downstream of Pe = 3 , a threshold that aligns closely with the upper limit of stress coupling we derived from our analysis. The Pe = 3 locations for these glaciers are: KAN = 32 km, AVA = 27 km, KUJ = 28 km, and EQP = 15 km . The values of λ chosen by minimizing the difference between observed and modeled velocities are 10 km for EQP and 20 km for KUJ. For KAN and AVA, seasonal velocities are not primarily driven by terminus variation, so we do not consider the estimation of λ valid for these glaciers. Notably, the λ values for EQP and KUJ are shorter than their respective Pe = 3 locations. We assume a constant upper limit of stress coupling along the profile and use a simplified one-dimensional terminus-driven model, both of which may influence our estimations.
Impact of observed seasonal elevation changes
The availability of seasonally-resolved elevation change allows us to investigate the degree to which outlet glacier terminus velocity is sensitive to changing surface elevation. We find that the seasonal elevation changes for EQP, KAN, KUJ, and AVA show limited spatial variation along the flow (Fig. ). Quantitatively, for each profile, the standard deviation along the profile (i.e., spatial variation with distance) averages only 22 % of the profile’s mean value. This low ratio indicates that the profiles are relatively uniform along flow. Thus, the seasonal elevation change is largely uniform across the DEM and these elevation changes do not significantly alter terminus-driven velocity (black lines in Figs. –). This aligns with our experimental results that suggest that vertical shifts in elevation have a limited contribution to velocity seasonality (Fig. a) but that steepening the surface elevation will cause stronger velocity responses to terminus variations (Fig. b).
Surface elevation difference profiles for AVA, KUJ, KAN, and EQP in 2019 relative to the October profile. Profile lengths are limited by the overlapping area of the DG-IS2-DEM.
Limitations and Outlook
We recognize that the choice of the terminus-driven model to explore how outlet glacier velocity might respond to terminus change is simplified over that from numerical results. There are several important limitations that merit further discussion. Seasonal elevation change is not yet available widely and may be worth including for some glaciers, particularly ones that are experiencing changes in slope near the terminus, which might influence ice velocity. While our experiments indicate insensitivity to uniform vertical elevation changes under the assumptions of this model, we cannot rule out that other physical models (e.g., including basal sliding or full stress coupling) might produce different results. Future work with higher-order models could test the robustness of this finding. The model only looks at the impact of terminus change on ice velocity and thus many important processes that control ice velocity are lumped into the examination of the velocity residuals. Here we have considered processes such as melt-driven changes in velocity and wintertime acceleration related to meltwater at the bed. Other factors such as air temperatures may play a secondary role, but terminus change and meltwater are likely the dominant drivers of seasonal and inter-annual variations.
In addition to these limitations, the model only captures longitudinal stresses, not lateral stresses. Since floating ice can experience velocity variations that are driven by nearby velocities on the grounded ice through transmission of lateral stress, the model is not appropriate for use on floating ice. The terminus-driven model importantly only applies close to the terminus since by definition, the frontal force F linearly decreases upstream of the terminus to zero at the upper limit of stress coupling.
Despite these shortcomings, the terminus-driven model provides a means to isolate the impact of seasonal terminus variations on outlet glacier velocity without numerical methods. Thus, it offers a quick assessment of how terminus variations impact glacier velocity across different glaciers, and how this impact changes over time. Future work could examine how seasonal velocity patterns vary from glacier to glacier and how those patterns change during periods of sustained year-over-year terminus retreat. Furthermore, the terminus-driven model helps illustrate that glaciers exhibit multiple, coexisting seasonal velocity patterns, suggesting that categorization of glacier velocity resulting from a single pattern may be insufficient to capture the full seasonal behavior.
Conclusions
We apply a terminus-driven model to elucidate the degree to which seasonal terminus change can explain seasonal velocity variations for four glaciers in Central West Greenland. The comparison between simulated and observed velocity suggests that all glaciers experience multiple drivers of velocity change. Some glaciers experience seasonal velocity changes that are dominated by terminus change and some glaciers do not respond to terminus change because either the terminus change is too small to be impactful and/or other processes dominate velocity, notably wintertime acceleration. Importantly, velocity residuals allow for us to see that runoff-driven acceleration occurs for all glaciers but its impact is limited to the early melt season. We present a new fusion of elevation products that resolve seasonal elevation changes for our study glaciers. When integrated into our model we do not find a significant difference between simulated velocities driven by terminus change assuming a fixed surface elevation versus those with an evolving surface elevation. We show that this is because the observed seasonal elevation change is largely the result of a vertical shift in surface elevation with no noticeable slope change, which would drive a velocity response. Our study highlights how a simple model can be used to examine the impact of terminus change on near-terminus velocity and allows a method for distinguishing the controls on seasonal velocity change of outlet glaciers. Moreover, the same framework could be applied to investigate the long-term changes in glacier dynamics with adequate historical data.
Code and data availability
Velocity time series data are available from . Velocity monthly mosaics are available from . ArcticDEMs are available from . BedMachineV5 is available from . QGreenland package is available from . The code, terminus data, flowlines, and DG-IS2-DEMs on which this article is based are available from.
The supplement related to this article is available online at https://doi.org/10.5194/tc-20-3875-2026-supplement.
Author contributions
EZ developed the code, performed the data processing, and wrote the manuscript. BS developed “DG-IS2-DEM”. GC, DF, BC, and DT advised EZ and revised the manuscript.
Competing interests
At least one of the (co-)authors is a member of the editorial board of The Cryosphere. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
We acknowledge the National Snow and Ice Data Center QGreenland package . We acknowledge DEMs provided by the Polar Geospatial Center under NSF-OPP awards 1043681, 1559691, 1542736, 1810976, and 2129685.
Financial support
This research has been supported by the NASA Headquarters (grant no. NNH20ZDA001N-ICESAT2) and the Jackson School of Geosciences, University of Texas at Austin (Postdoctoral Fellowship).
Review statement
This paper was edited by Sainan Sun and reviewed by Twila Moon and two anonymous referees.
ReferencesAn, L., Rignot, E., Mouginot, J., and Millan, R.: A Century of Stability of Avannarleq and Kujalleq Glaciers, West Greenland, Explained Using High-Resolution Airborne Gravity and Other Data, Geophys. Res. Lett., 45, 3156–3163, 10.1002/2018GL077204, 2018.Andrews, L. C., Catania, G. A., Hoffman, M. J., Gulley, J. D., Lüthi, M. P., Ryser, C., Hawley, R. L., and Neumann, T. A.: Direct observations of evolving subglacial drainage beneath the Greenland Ice Sheet, Nature, 514, 80–83, 10.1038/nature13796, 2014.Bartholomaus, T. C., Anderson, R. S., and Anderson, S. P.: Response of glacier basal motion to transient water storage, Nat. Geosci., 1, 33–37, 10.1038/ngeo.2007.52, 2008.Bartholomaus, T. C., Anderson, R. S., and Anderson, S. P.: Growth and collapse of the distributed subglacial hydrologic system of Kennicott Glacier, Alaska, USA, and its effects on basal motion, J. Glaciol., 57, 9850–1002, 10.3189/002214311798843269, 2011.Bartholomew, I., Nienow, P., Mair, D., Hubbard, A., King, M. A., and Sole, A.: Seasonal evolution of subglacial drainage and acceleration in a Greenland outlet glacier, Nat. Geosci., 3, 408–411, 10.1038/ngeo863, 2010.Bézu, C. and Bartholomaus, T. C.: Greenland Ice Sheet's Distinct Calving Styles Are Identified in Terminus Change Timeseries, Geophys. Res. Lett., 51, 10.1029/2024gl110224, 2024.Black, T. E. and Joughin, I.: Weekly to monthly terminus variability of Greenland's marine-terminating outlet glaciers, The Cryosphere, 17, 1–13, 10.5194/tc-17-1-2023, 2023.Boon, S. and Sharp, M.: The role of hydrologically-driven ice fracture in drainage system evolution on an Arctic glacier, Geophys. Res. Lett., 30, 10.1029/2003GL018034, 2003.Carnahan, E., Catania, G., and Bartholomaus, T. C.: Observed mechanism for sustained glacier retreat and acceleration in response to ocean warming around Greenland, The Cryosphere, 16, 4305–4317, 10.5194/tc-16-4305-2022, 2022.Carr, J. R., Vieli, A., and Stokes, C.: Influence of sea ice decline, atmospheric warming, and glacier width on marine-terminating outlet glacier behavior in northwest Greenland at seasonal to interannual timescales, J. Geophys. Res.: Earth Surf., 118, 1210–1226, 10.1002/jgrf.20088, 2013.Catania, G. A., Neumann, T. A., and Price, S. F.: Characterizing englacial drainage in the ablation zone of the Greenland ice sheet, J. Glaciol., 54, 567–578, 10.3189/002214308786570854, 2017.Catania, G. A., Stearns, L. A., Sutherland, D. A., Fried, M. J., Bartholomaus, T. C., Morlighem, M., Shroyer, E., and Nash, J.: Geometric controls on tidewater glacier retreat in central western Greenland, J. Geophys. Res.: Earth Surf., 123, 2024–2038, 10.1029/2017JF004499, 2018.Davison, B. J., Sole, A. J., Cowton, T. R., Lea, J. M., Slater, D. A., Fahrner, D., and Nienow, P. W.: Subglacial Drainage Evolution Modulates Seasonal Ice Flow Variability of Three Tidewater Glaciers in Southwest Greenland, J. Geophys. Res.: Earth Surf., 125, e2019JF005492, 10.1029/2019JF005492, 2020.Enderlin, E. M., Howat, I. M., Jeong, S., Noh, M.-J., van Angelen, J. H., and van den Broeke, M. R.: An improved mass budget for the Greenland ice sheet, Geophys. Res. Lett., 41, 866–872, 10.1002/2013GL059010, 2014.Felikson, D., Bartholomaus, T. C., Catania, G. A., Korsgaard, N. J., Kjær, K. H., Morlighem, M., Noël, B., Van Den Broeke, M., Stearns, L. A., Shroyer, E. L., and Sutherland, D. A.: Inland thinning on the Greenland ice sheet controlled by outlet glacier geometry, Nat. Geosci., 10, 366–369, 10.1038/ngeo2934, 2017.Felikson, D., A. Catania, G., Bartholomaus, T. C., Morlighem, M., and Noël, B. P. Y.: Steep Glacier Bed Knickpoints Mitigate Inland Thinning in Greenland, Geophys. Res. Lett., 48, e2020GL090112, 10.1029/2020GL090112, 2021.Felikson, D., Nowicki, S., Nias, I., Morlighem, M., and Seroussi, H.: Seasonal Tidewater Glacier Terminus Oscillations Bias Multi-Decadal Projections of Ice Mass Change, J. Geophys. Res.: Earth Surf., 127, e2021JF006249, 10.1029/2021JF006249, 2022.Fried, M. J., Catania, G. A., Stearns, L. A., Sutherland, D. A., Bartholomaus, T. C., Shroyer, E., and Nash, J.: Reconciling Drivers of Seasonal Terminus Advance and Retreat at 13 Central West Greenland Tidewater Glaciers, J. Geophys. Res.: Earth Surf., 123, 1590–1607, 10.1029/2018JF004628, 2018.Gardner, A. S., Moholdt, G., Scambos, T., Fahnstock, M., Ligtenberg, S., Van Den Broeke, M., and Nilsson, J.: Increased West Antarctic and unchanged East Antarctic ice discharge over the last 7 years, The Cryosphere, 12, 521–547, 10.5194/tc-12-521-2018, 2018.Gardner, A. S., Fahnestock, M. A., and Scambos, T. A.: ITS_LIVE Regional Glacier and Ice Sheet Surface Velocities (Version 1), National Snow and Ice Data Center [data set], 10.5067/6II6VW8LLWJ7, 2023.Goliber, S. A. and Catania, G. A.: Glacier Terminus Morphology Informs Calving Style, Geophys. Res. Lett., 51, 10.1029/2024gl108530, 2024.Holland, D. M., Thomas, R. H., De Young, B., Ribergaard, M. H., and Lyberth, B.: Acceleration of Jakobshavn Isbræ triggered by warm subsurface ocean waters, Nat. Geosci., 1, 659–664, 10.1038/ngeo316, 2008.Howat, I. M., Joughin, I., Tulaczyk, S., and Gogineni, S.: Rapid retreat and acceleration of Helheim Glacier, east Greenland, Geophys. Res. Lett., 32, L22502, 10.1029/2005GL024737, 2005.Howat, I. M., Joughin, I., Fahnestock, M., Smith, B. E., and Scambos, T. A.: Synchronous retreat and acceleration of southeast Greenland outlet glaciers 2000–06: ice dynamics and coupling to climate, J. Glaciol., 54, 646–660, 10.3189/002214308786570908, 2008.IPCC: Climate Change 2022: Mitigation of Climate Change, Cambridge University Press, Cambridge, UK, and New York, NY, USA, 10.1017/9781009157926, 2022.Joughin, I.: MEaSUREs Greenland Monthly Ice Sheet Velocity Mosaics from SAR and Landsat (Version 5), National Snow and Ice Data Center [data set], https://nsidc.org/data/NSIDC-0731/versions/5 (last access: 11 July 2026), 2023.Joughin, I., Das, S. B., King, M. A., Smith, B. E., Howat, I. M., and Moon, T.: Seasonal speedup along the Western flank of the Greenland Ice Sheet, Science, 320, 781–783, 10.1126/science.1153288, 2008.Joughin, I., Smith, B. E., Howat, I. M., Floricioiu, D., Alley, R. B., Truffer, M., and Fahnestock, M.: Seasonal to decadal scale variations in the surface velocity of Jakobshavn Isbrae, Greenland: Observation and model-based analysis, J. Geophys. Res.: Earth Surf., 117, F02030, 10.1029/2011JF002110, 2012.Kane, E., Rignot, E., Mouginot, J., Millan, R., Li, X., Scheuchl, B., and Fahnestock, M.: Impact of Calving Dynamics on Kangilernata Sermia, Greenland, Geophys. Res. Lett., 47, e2020GL088524, 10.1029/2020GL088524, 2020.Karlsson, N. B., Solgaard, A. M., Mankoff, K. D., Gillet-Chaulet, F., MacGregor, J. A., Box, J. E., Citterio, M., Colgan, W. T., Larsen, S. H., Kjeldsen, K. K., Korsgaard, N. J., Benn, D. I., Hewitt, I. J., and Fausto, R. S.: A first constraint on basal melt-water production of the Greenland ice sheet, Nat. Commun., 12, 3461, 10.1038/s41467-021-23739-z, 2021.King, M. D., Howat, I. M., Jeong, S., Noh, M. J., Wouters, B., Noël, B., and van den Broeke, M. R.: Seasonal to decadal variability in ice discharge from the Greenland Ice Sheet, The Cryosphere, 12, 3813–3825, 10.5194/tc-12-3813-2018, 2018.King, M. D., Howat, I. M., Candela, S. G., Noh, M. J., Jeong, S., Noël, B. P., van den Broeke, M. R., Wouters, B., and Negrete, A.: Dynamic ice loss from the Greenland Ice Sheet driven by sustained glacier retreat, Commun. Earth Environ., 1, 1, 10.1038/s43247-020-0001-2, 2020.Kneib-Walter, A., Lüthi, M. P., Moreau, L., and Vieli, A.: Drivers of Recurring Seasonal Cycle of Glacier Calving Styles and Patterns, Front. Earth Sci., 9, 10.3389/feart.2021.667717, 2021.Kneib-Walter, A., Lüthi, M. P., Funk, M., Jouvet, G., and Vieli, A.: Observational constraints on the sensitivity of two calving glaciers to external forcings, J. Glaciol., 69, 459–474, 10.1017/jog.2022.74, 2023.Lüthi, M. P., Vieli, A., Moreau, L., Joughin, I., Reisser, M., Small, D., and Stober, M.: A century of geometry and velocity evolution at Eqip Sermia, West Greenland, J. Glaciol., 62, 640–654, 10.1017/jog.2016.38, 2016.Medley, B., Neumann, T. A., Zwally, H. J., Smith, B. E., and Stevens, C. M.: Simulations of firn processes over the Greenland and Antarctic ice sheets: 1980–2021, The Cryosphere, 16, 3971–4011, 10.5194/tc-16-3971-2022, 2022.Moon, T., Joughin, I., Smith, B., and Howat, I.: 21st-Century evolution of Greenland outlet glacier velocities, Science, 336, 576–578, 10.1126/science.1219985, 2012.Moon, T., Joughin, I., Smith, B., van den Broeke, M. R., van de Berg, W. J., Noël, B., and Usher, M.: Distinct patterns of seasonal Greenland glacier velocity, Geophys. Res. Lett., 41, 7209–7216, 10.1002/2014GL061836, 2014.Moon, T., Joughin, I., and Smith, B.: Seasonal to multiyear variability of glacier surface velocity, terminus position, and sea ice/ice mélange in northwest Greenland, J. Geophys. Res.: Earth Surf., 120, 818–833, 10.1002/2015JF003494, 2015.Moon, T., Gardner, A., Csatho, B., Parmuzin, I., and Fahnestock, M.: Rapid reconfiguration of the Greenland Ice Sheet coastal margin, J. Geophys. Res.: Earth Surf., 125, 10.1029/2020jf005585, 2020.Moon, T., Fisher, M., Stafford, T., and Thurber, A.: QGreenland (v3), National Snow and Ice Data Center [data set], 10.5281/zenodo.8326507, 2023.Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: BedMachine v3: Complete bed topography and ocean bathymetry mapping of Greenland from multibeam echo dounding combined with mass conservation, Geophys. Res. Lett., 44, 11,051–11,061, 10.1002/2017GL074954, 2017.Morlighem, M., Williams, C. N., Rignot, E., An, L., Arndt, J. E., Bamber, J. L., Catania, G., Chauché, N., Dowdeswell, J. A., Dorschel, B., Fenty, I., Hogan, K., Howat, I., Hubbard, A., Jakobsson, M., Jordan, T. M., Kjeldsen, K. K., Millan, R., Mayer, L., Mouginot, J., Noël, B. P. Y., O'Cofaigh, C., Palmer, S., Rysgaard, S., Seroussi, H., Siegert, M. J., Slabon, P., Straneo, F., van den Broeke, M. R., Weinrebe, W., Wood, M., and Zinglersen, K. B.: IceBridge BedMachine Greenland, Version 5, National Snow and Ice Data Center [data set], 10.5067/GMEVBWFLWA7X, 2022.Noël, B., van de Berg, W. J., Lhermitte, S., and Van Den Broeke, M.: Rapid ablation zone expansion amplifies north Greenland mass loss, Sci. Adv., 5, eaaw0123, 10.1126/sciadv.aaw0123, 2019.Poinar, K.: Seasonal Flow Types of Glaciers in Sermilik Fjord, Greenland, Over 2016–2021, J. Geophys. Res.: Earth Surf., 128, e2022JF006901, 10.1029/2022JF006901, 2023.Porter, C., Morin, P., Howat, I. M., Noh, M.-J., Bates, B., Peterman, K., Keesey, S., Schlenk, M., Gardiner, J., Tomko, K., Willis, M. J., Kelleher, C., Cloutier, M., Husby, E., Foga, S., Nakamura, H., Platson, M., Wethington, M. J., Williamson, C. J., Bauer, G., Enos, J., Arnold, G., Kramer, W., Becker, P., Doshi, A., D'Souza, C., Cummens, P., Laurier, F., and Bojesen, M.: ArcticDEM – Strips (Version 4.1), Harvard Dataverse [data set], 10.7910/DVN/C98DVS, 2022.Robel, A. A., Roe, G. H., and Haseloff, M.: Response of marine-terminating glaciers to forcing: time scales, sensitivities, instabilities and stochastic dynamics, J. Geophys. Res.: Earth Surf., 1–62, 10.1029/2018jf004709, 2018.Schoof, C. G.: Ice-sheet acceleration driven by melt supply variability, Nature, 468, 803 806, 10.1038/nature09618, 2010.Sergienko, O. V.: Marine outlet glacier dynamics, steady states and steady-state stability, J. Glaciol., 68, 946–960, 10.1017/jog.2022.13, 2022.Shepherd, A., Ivins, E. R., A, G., Barletta, V. R., Bentley, M. J., Bettadpur, S., Briggs, K. H., Bromwich, D. H., Forsberg, R., Galin, N., Horwath, M., Jacobs, S., Joughin, I., King, M. A., Lenaerts, J. T. M., Li, J., Ligtenberg, S. R. M., Luckman, A., Luthcke, S. B., McMillan, M., Meister, R., Milne, G., Mouginot, J., Muir, A., Nicolas, J. P., Paden, J., Payne, A. J., Pritchard, H., Rignot, E., Rott, H., Sørensen, L. S., Scambos, T. A., Scheuchl, B., Schrama, E. J. O., Smith, B., Sundal, A. V., van Angelen, J. H., van de Berg, W. J., van den Broeke, M. R., Vaughan, D. G., Velicogna, I., Wahr, J., Whitehouse, P. L., Wingham, D. J., Yi, D., Young, D., and Zwally, H. J.: A reconciled estimate of ice-sheet mass balance, Science, 338, 1183–1189, 10.1126/science.1228102, 2012.Smith, L. C., Yang, K., Pitcher, L. H., Overstreet, B. T., Chu, V. W., Åsa K. Rennermalm, Ryan, J. C., Cooper, M. G., Gleason, C. J., Tedesco, M., Jeyaratnam, J., van As, D., van den Broeke, M. R., van de Berg, W. J., Noël, B., Langen, P. L., Cullather, R. I., Zhao, B., Willis, M. J., Hubbard, A., Box, J. E., Jenner, B. A., and Behar, A. E.: Direct measurements of meltwater runoff on the Greenland ice sheet surface, Proc. Natl. Aca. Sci. USA, 114, E10622–E10631, 10.1073/pnas.1707743114, 2017.Solgaard, A. M., Rapp, D., Noël, B. P. Y., and Hvidberg, C. S.: Seasonal Patterns of Greenland Ice Velocity From Sentinel-1 SAR Data Linked to Runoff, Geophys. Res. Lett., 49, e2022GL100343, 10.1029/2022GL100343, 2022.Stevens, L. A., Behn, M. D., Das, S. B., Joughin, I., Noël, B. P. Y., van den Broeke, M. R., and Herring, T.: Greenland Ice Sheet flow response to runoff variability, Geophys. Res. Lett., 43, 11295–11303, 10.1002/2016GL070414, 2016.Stevens, L. A., Nettles, M., Davis, J. L., Creyts, T. T., Kingslake, J., Hewitt, I. J., and Stubblefield, A.: Tidewater-glacier response to supraglacial lake drainage, Nat. Commun., 13, 6065, 10.1038/s41467-022-33763-2, 2022.Sundal, A. V., Shepherd, A., Nienow, P., Hanna, E., Palmer, S., and Huybrechts, P.: Melt-induced speed-up of Greenland ice sheet offset by efficient subglacial drainage, Nature, 469, 521–524, 10.1038/nature09740, 2011.The GlaMBIE Team: Community estimate of global glacier mass changes from 2000 to 2023, Nature, 639, 382–388, 10.1038/s41586-024-08545-z, 2025.Ultee, L., Felikson, D., Minchew, B., Stearns, L. A., and Riel, B.: Helheim Glacier ice velocity variability responds to runoff and terminus position change at different timescales, Nat. Commun., 13, 6022, 10.1038/s41467-022-33292-y, 2022.van den Broeke, M. R., Enderlin, E. M., Howat, I. M., Kuipers Munneke, P., Noël, B. P. Y., van de Berg, W. J., van Meijgaard, E., and Wouters, B.: On the recent contribution of the Greenland ice sheet to sea level change, The Cryosphere, 10, 1933–1946, 10.5194/tc-10-1933-2016, 2016. van der Veen, C. J.: Fundamentals of glacier dynamics, 2nd ed., CRC Press, New York, 10.1201/b14059, 2013. van der Veen, C. J., Plummer, J., and Stearns, L. A.: Controls on the recent speed-up of Jakobshavn Isbrae, West Greenland, J. Glaciol., 57, 770–782, 2011.Vieli, A., Jania, J., Blatter, H., and Funk, M.: Short-term velocity variations on Hansbreen, a tidewater glacier in Spitsbergen, J. Glaciol., 50, 389–398, 10.3189/172756504781829963, 2004.Vijay, S., Khan, S. A., Kusk, A., Solgaard, A. M., Moon, T., and Bjørk, A. A.: Resolving Seasonal Ice Velocity of 45 Greenlandic Glaciers With Very High Temporal Details, Geophys. Res. Lett., 46, 1485–1495, 10.1029/2018GL081503, 2019.Vijay, S., King, M. D., Howat, I. M., Solgaard, A. M., Khan, S. A., and Noël, B.: Greenland ice-sheet wide glacier classification based on two distinct seasonal ice velocity behaviors, J. Glaciol., 67, 1241–1248, 10.1017/jog.2021.89, 2021.Walter, J. I., Box, J. E., Tulaczyk, S. M., Brodsky, E. E., Howat, I. M., Ahn, Y., and Brown, A. A.: Oceanic mechanical forcing of a marine-terminating Greenland glacier, Ann. Glaciol., 53, 181–192, 10.3189/2012aog60a083, 2012.Werder, M. A., Hewitt, I. J., Schoof, C. G., and Flowers, G. E.: Modeling channelized and distributed subglacial drainage in two dimensions, J. Geophys. Res.: Earth Surf., 118, 2140–2158, 10.1002/jgrf.20146, 2013.Wood, M., Rignot, E., Fenty, I., An, L., Bjørk, A., Van Den Broeke, M., Cai, C., Kane, E., Menemenlis, D., Millan, R., Morlighem, M., Mouginot, J., Noël, B., Scheuchl, B., Velicogna, I., Willis, J. K., and Zhang, H.: Ocean forcing drives glacier retreat in Greenland, Sci. Adv., 7, eaba7282, 10.1126/sciadv.aba7282, 2021.Zhang, E.: Compounding seasonal variations in outlet glacier dynamics revealed by high-resolution observations, Zenodo [data set], 10.5281/zenodo.8428196, 2023.Zhang, E., Liu, L., and Huang, L.: Automatically delineating the calving front of Jakobshavn Isbræ from multitemporal TerraSAR-X images: a deep learning approach, The Cryosphere, 13, 1729–1741, 10.5194/tc-13-1729-2019, 2019.Zhang, E., Catania, G., and Trugman, D. T.: AutoTerm: an automated pipeline for glacier terminus extraction using machine learning and a “big data” repository of Greenland glacier termini, The Cryosphere, 17, 3485–3503, 10.5194/tc-17-3485-2023, 2023.