INTRODUCTION
Ocean tides cause variations in observed glacier flow in Antarctica, Greenland and other glaciated regions (e.g. Walters, Reference Walters1989; Harrison and others, Reference Harrison, Echelmeyer and Engelhardt1993; Reeh and others, Reference Reeh, Mayer, Olesen, Christensen and Thomsen2000; Podrasky and others, Reference Podrasky, Truffer, Lüthi and Fahnestock2014). These flow variations have different amplitudes, periods and phases at different glaciers, indicating that there are multiple physical mechanisms through which tides influence ice-sheet flow. At Whillans Ice Stream, almost all flow occurs through tidally-modulated stick-slip events over basal ‘sticky spots’ (Bindschadler and others, Reference Bindschadler, King, Alley, Anandakrishnan and Padman2003; Wiens and others, Reference Wiens, Anandakrishnan, Winberry and King2008; Lipovsky and Dunham, Reference Lipovsky and Dunham2017). At Bindschadler Ice Stream (hereafter BIS), ice flow primarily varies at the same period as the dominant local ocean tides and may be detected over 100 km upstream of the grounding line (Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997; Anandakrishnan and others, Reference Anandakrishnan, Voigt, Alley and King2003). At Rutford Ice Stream (hereafter RIS), fortnightly ocean tides (known as M f and M sf tides) are weak compared with the dominant high-frequency diurnal and semidiurnal tides. Ice flow velocities in RIS, however, vary almost entirely at the M sf period (14.77 days), which corresponds to the beat frequency of the two primary semidiurnal tidal constituents. These M sf velocity variations in RIS are nearly 25% of the time-averaged flow rate within the ice shelf, and can be detected more than 100 km upstream of the grounding line (Gudmundsson, Reference Gudmundsson2006; Murray and others, Reference Murray, Smith, King and Weedon2007; Brunt and others, Reference Brunt, King, Fricker and MacAyeal2010; King and others, Reference King, Makinson and Gudmundsson2011; Rosier and others, Reference Rosier, Gudmundsson and Green2015; Minchew and others, Reference Minchew, Simons, Riel and Milillo2017).
Many previous studies have suggested that tides can cause flow variations in an ice shelf through changes in hydrostatic and flexural stresses (e.g. Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997; Gudmundsson, Reference Gudmundsson2007; Goldberg and others, Reference Goldberg, Schoof and Sergienko2014; Thompson and others, Reference Thompson, Simons and Tsai2014; Rosier and Gudmundsson, Reference Rosier and Gudmundsson2016). Variations in basal stress from tidal currents and variations in driving stress due to changing ice shelf surface slope (Makinson and others, Reference Makinson, King, Nicholls and Hilmar Gudmundsson2012) have also been put forward as potential tidal mechanisms for producing ice shelf flow variability, but have had limited success in reproducing observations in model investigations (Brunt and MacAyeal, Reference Brunt and MacAyeal2014). Once the tidal signal enters the grounded ice stream, models suggest that nonlinearity in the basal sliding law is capable of producing flow variations at periods longer than the tidal forcing in the grounded ice stream, which then propagate up and downstream (Gudmundsson, Reference Gudmundsson2007, Reference Gudmundsson2011; Williams and others, Reference Williams, Hindmarsh and Arthern2012; Rosier and others, Reference Rosier, Gudmundsson and Green2014; Walker and others, Reference Walker, Parizek, Alley, Brunt and Anandakrishnan2014). Yet further studies have suggested that ocean tides lead to variations in ice shelf flow through a decrease in contact with sub-shelf pinning points and sidewalls (Schmeltz and others, Reference Schmeltz, Rignot and MacAyeal2001; Doake and others, Reference Doake2002; Heinert and Riedel, Reference Heinert and Riedel2007; Thomas, Reference Thomas2007; Minchew and others, Reference Minchew, Simons, Riel and Milillo2017). However, the extent to which such tidally-modulated changes in ice shelf buttressing combine with other potential tidal forcing mechanisms to produce the period, amplitude and phase of observed changes in ice shelf flow has not been tested in a systematic fashion.
Notwithstanding such modeling advances, the full range of mechanisms that cause ice flow to vary with tides, and their interaction with each other, is still unknown. As GPS stations are deployed for longer durations, more glaciers and ice shelves are being found to exhibit low-frequency variations in ice flow (Rosier and others, Reference Rosier, Gudmundsson and Green2015). Existing models have proposed explanations for some characteristics of tidal ice flow variations on some glaciers, but not others. Spatially-resolved observations are also calling into question existing theories of low-frequency ice flow variations. For example, new observations from interferometric synthetic aperture radar (InSAR) (Minchew and others, Reference Minchew, Simons, Riel and Milillo2017) show that the low-frequency tidal signal in ice flow at RIS originates in the ice shelf, rather than in the grounded ice stream. These observations also indicate that tidal variations in ice shelf strain rate are not uniform, but exhibit significant small-scale variability, suggesting that ice flow variability originates with processes within the ice shelf, rather than being passively advected downstream from the grounded ice stream. Here, we propose explanations for the period, amplitude and phase of tidal signals in ice shelf and ice stream flow. We start by using observations to argue that tides modulate the contact stresses near the grounding line and between the ice shelf and sub-shelf bathymetric pinning points. Using a lumped Maxwell model, we show that nonlinearities in ice shelf buttressing with respect to tidal height are the most plausible explanations of low-frequency variations in ice shelf flow. We also find that the phase of high-frequency variations in displacement relative to tidal height indicates whether variations in buttressing or hydrostatic stresses are more important for the modulation of ice flow.
OBSERVATIONAL CONSTRAINTS AND MODEL MOTIVATION
Observations indicate that a rising ocean tide causes grounding line retreat and lifts the ice shelf away from sub-shelf bathymetric pinning points, decreasing resistance to the flow of the ice shelf (Schmeltz and others, Reference Schmeltz, Rignot and MacAyeal2001; Doake and others, Reference Doake2002; Minchew and others, Reference Minchew, Simons, Riel and Milillo2017). Conversely, a rising tide also increases hydrostatic stresses in the ice shelf, which increases the resistance to ice flow (e.g. Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997; Gudmundsson, Reference Gudmundsson2007; Rosier and Gudmundsson, Reference Rosier and Gudmundsson2016). In this section, we analyze observations of tidally-modulated ice shelf flow variations at two West Antarctic ice streams and discuss which physical processes might play a role in causing these variations.
At RIS, horizontal ice shelf displacement (blue line in Fig. 1a) varies primarily at fortnightly periods, with smaller semidiurnal variations in displacement that are nearly in phase with tidal height (calculated to be 15° phase lag, though not necessarily discernible in Fig. 1a). At BIS (Fig. 1b), horizontal ice displacement varies at fortnightly and diurnal periods with nearly equal amplitude, with diurnal variations completely out of phase with tidal height (210° phase lag). At both ice streams, there is a significant low-frequency signal in horizontal displacement downstream of the grounding line. However, the variations in displacement occurring at high tidal frequencies have different phase with respect to tidal height. We aim to explain the causes of these similarities and differences between the tidally-induced ice shelf flow variations at these two ice streams.
Data from Gudmundsson (Reference Gudmundsson2006) and personal communication from R. Bindschadler and K. Brunt.
We first discuss different ways in which tides cause variations in the stress applied on an ice shelf. To do so, we analyze observations of the relationship between horizontal strain and tidal height at RIS and BIS (Figs. 1c and d, respectively). Strain is measured over a baseline distance between stations near and downstream of the grounding lines of RIS and BIS (with the understanding that ice stream grounding line locations vary with tidal height and measurements of grounding line position are subject to significant errors, see Rignot and others, Reference Rignot, Mouginot and Scheuchl2011). Observations are indicated by grey circles, and the least squares fit of the relationship between tidal height and ice shelf strain is indicated by the black solid line (black dashed lines indicate the ±2σ least squares fits), with the functional form written explicitly in the inset boxes. For elastic deformations, strain is proportional to the imposed stress. If viscous deformation is small compared with elastic deformation, then the phase of high-frequency variations in observed strain with respect to tidal height are expected to be close to either 0° or 180° (see the section ‘Parameter exploration of tidal influence on ice shelf flow’, below). In such a scenario, measured variations in strain can be used as indicators of variations in stress. In the floating portion of RIS (Fig. 1c), strain is positively correlated (close to 0° in phase) and nonlinear with respect to tidal height. This nonlinearity indicates that strain is more sensitive to changing tides at high tide than at low tide. In contrast, strain in the floating portion of BIS (Fig. 1d) is negatively correlated (close to 180° out of phase) and nonlinear with respect to tidal height, with higher sensitivity to changing tides at low tide.
Variations in hydrostatic (and flexural) stresses within an ice shelf are linearly proportional to local tidal height (Holdsworth, Reference Holdsworth1969; Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997; Gudmundsson, Reference Gudmundsson2007; Rosier and Gudmundsson, Reference Rosier and Gudmundsson2016). Thus, our primary goal in this study is to determine mechanisms which can plausibly produce the observed nonlinear relationships between ice shelf strain and tidal height in the ice shelves at RIS and BIS. Spatially-resolved InSAR observations indicate that these tidal strain variations are localized in nature (Minchew and others, Reference Minchew, Simons, Riel and Milillo2017), suggesting that the mechanism causing ice shelf flow variations is unrelated to large-scale processes including either changes in ice shelf slopes and driving stresses (Makinson and others, Reference Makinson, King, Nicholls and Hilmar Gudmundsson2012) or changes in the velocity of the grounded ice stream, which advect downstream to the ice shelf (Gudmundsson, Reference Gudmundsson2011). We argue that beyond hydrostatic, flexural and driving stresses, the only way for tides to cause localized variations in ice shelf strain is through variations in resistance from localized contact at the lower and lateral ice shelf surfaces. The total of the stresses exerted back on the grounded ice stream by all different contact points on these pinned ice shelf surfaces is typically referred to as the ‘buttressing stress’ (Goldberg and others, Reference Goldberg, Holland and Schoof2009).
We consider the locations of ‘final contact points’ where ice goes afloat and loses contact with the bed because the weight of ice is approximately equal (or exactly equal in the hydrostatic limit) to the hydrostatic pressure of ocean water. These final contact points occur at (what is conventionally called) the grounding line, and also at pinning points within the ice shelf where local maxima in bathymetry come into contact with the lower ice shelf surface (see schematic example in Fig. 2). Tidal variations in hydrostatic pressure lead to migration of these final contact points and change the total area where ice-bed contact occurs (as previously argued by Anandakrishnan and Alley, Reference Anandakrishnan and Alley1997). By simply considering where hydrostatic equilibrium occurs near the grounding line and in ice shelves, it has been shown (Tsai and Gudmundsson, Reference Tsai and Gudmundsson2015) that final contact points (either a grounding line or pinning point) migrate ~10 times further upstream from mean to high tide than from low to mean tide. This asymmetry arises because spatial variations in ice thickness are set by slow viscous deformation of the ice sheet. However, on the typical timescale of tidal hydrostatic pressure variations (12–24 h), viscous deformation is too slow to adjust ice thickness in the vicinity of the final contact points in response to the accumulation of hydrostatic stresses.
An explanatory schematic (Fig. 2) shows the asymmetric migration of the grounding line and a pinning point from low to mean to high tide. The precise functional dependence of the total contact area on tidal height will depend on the particulars of local bathymetry. However, if there is non-zero basal shear stress that is approximately spatially uniform over all locations of ice contact, the decrease in the spatially-integrated basal shear stress associated with grounding line retreat will be greater in going from mean to high tide than from low to mean tide. The consequence of this asymmetry is a nonlinear relationship between tidal height and the total buttressing stress within the ice shelf. Indeed, at at both RIS and BIS, ice shelf strain is a nonlinear function of tidal height (Figs. 1c and d). To capture this asymmetry, we formulate (in the next section) a very general empirical parameterization for total ice shelf buttressing stress as an increasing function of tidal height with arbitrary shape. We then use observations of the relationship between tidal height and ice shelf strain (e.g. Figs. 1c and d) to constrain the shape of this function, and explore the consequences of this nonlinearity on the observed ice shelf strain using a lumped viscoelastic ice shelf model.
VISCOELASTIC ICE SHELF MODEL
To determine the relative importance of variations in hydrostatic and buttressing stresses in generating the period, phase and amplitude of observed ice shelf flow variability, we use a lumped viscoelastic model. For simplicity, the normalized tidal height is assumed to be composed of two high-frequency components of equal amplitude
By only forcing the model at high frequencies, we can isolate the effect of nonlinearities in the ice shelf dynamics in producing a low-frequency ice shelf response. Also, in reality, high-frequency tidal constituents may have different amplitudes. Regardless, simulations forced directly by tidal heights observed in ice shelves of interest (e.g. red lines in Figs.1a and 1b) do not produce a qualitatively different ice flow response than simulations forced by the simpler tidal height in Eqn (1).
Considering the ice shelf as a simple elastic beam, tidal variations in water level will cause stress variations within the ice shelf through hydrostatic pressure and ice flexure. Changes in hydrostatic stresses within the ice shelf are proportional to tidal height
where Δh is the tidal range, ρ w is the density of ocean water, g is the acceleration due to gravity and negative values denote horizontal compression of the ice shelf (Rosier and others, Reference Rosier, Gudmundsson and Green2014). Changes in longitudinal stress associated with ice shelf flexure at the lower ice surface around a non-migrating grounding line are proportional to tidal height
where $L_x^{ - 1} $ is a flexure length scale dependent on ice shelf rheology (Holdsworth, Reference Holdsworth1969). However, flexural stresses quickly drop off within the ice shelf (Rosier and Gudmundsson, Reference Rosier and Gudmundsson2016), and average to zero when integrated over the ice shelf thickness. Since we consider such a depth-integrated, lumped representation of an ice shelf in this study, flexural stresses are negligible and thus we only consider hydrostatic stresses,
where σ h0 is the magnitude of tidal variations in hydrostatic stress within the ice shelf, which cause compression of the ice shelf at high tide and extension of the ice shelf at low tide.
Downstream of the grounding lines of RIS and BIS, measured ice shelf strain varies nonlinearly with tidal height (Fig. 1), likely as a consequence of tidal asymmetries in grounding line migration and the contact between the ice shelf and sub-shelf bathymetry (as discussed in the section ‘Observational constraints and model motivation’). Thus, we can use the shape of the measured relationship between ice shelf strain and tidal height, captured by the (black line) nonlinear least squares estimates in Figs. 1c and d, as an empirical measure of the shape of this asymmetry. The equivalent buttressing stress (which as we argue above, is proportional to strain for elastic deformations) is parameterized as an increasing function with arbitrary shape
where α is a shape parameter indicating the strength of the asymmetry in buttressing stress with respect to tidal height (the same shape parameter that is fit in Figs. 1c and d), and
is a non-dimensional parameter measuring the relative importance of the magnitude of changes in buttressing stress (σ b0) to the magnitude of changes in hydrostatic stress (σ h0). The value of β cannot be determined strictly from the strain/height relationship in Figs. 1c and d, since measured strains include the effect of hydrostatic stress variations. The bracketed term in Eqn (5) varies between −1 to 1 regardless of the value chosen for α. It should be noted that this parameterized buttressing equation could just as easily have taken other forms, such as a polynomial. Here, we choose this particular functional form of Eqn (5) in order to easily parameterize the magnitude (β) and shape (α) of buttressing stress variations with a minimum of parameters and without introducing other potential sources of nonlinearity or asymmetry with respect to tidal height. It is not necessarily the case that the buttressing relationship will be nonlinear (as this form permits the possibility that α = 1), but rather, we use the observations of nonlinearities in ice shelf strain as a guide for setting the value of α. In the section ‘Parameter exploration of tidal influence on ice shelf flow’, we further explore the sensitivity of the modeled ice flow response to the values of α and β.
We simulate the tidal variations in local ice shelf deformation using a lumped Maxwell model (Tsai and others, Reference Tsai, Rice and Fahnestock2008) for a viscoelastic beam unconfined in the y and z direction and driven by tidal variations in stress in the x direction
where ε xx is the longitudinal strain, E is the Young's modulus of ice and η e is the effective viscosity within the ice shelf (MacAyeal and others, Reference MacAyeal, Sergienko and Banwell2015). The model is ‘lumped’ in that it is forced by the average stresses in a section of the ice stream that includes the grounding line and some part of the ice shelf. The effective viscosity in this section of the ice stream is a function of the time-averaged (indicated by overbar) second invariant of the deviatoric stress tensor
where we expect that the longitudinal and lateral shear stresses will dominate for ice shelves (as in Goldberg and others, Reference Goldberg, Schoof and Sergienko2014). We separate tidally-driven variations in hydrostatic (σ h) and buttressing stresses (σ b) from the sum of the squares of the time-averaged deviatoric stresses ( $\bar I_2 $ ,) giving an effective viscosity
where A and n are the rate factor and exponent in the Nye–Glen constitutive law (Cuffey and Paterson, Reference Cuffey and Paterson2010), respectively, and
is a non-dimensional parameter measuring the relative importance of the magnitude of total time-averaged stresses in the ice shelf (such as driving and lateral shear stresses, which determine the time-averaged ice flow) to the magnitude of tidal variations in hydrostatic stress. The driving and other background stresses within the ice shelf are only included in the effective viscosity and nowhere else in the Maxwell model, as we intend to simulate the tidal anomalies in ice shelf deformation, but not the time-averaged ice flow. We choose to separate tidally-varying from time-averaged terms in the effective viscosity in order to simplify the effective viscosity terms, but this choice has no significant quantitative effect on the results discussed hereafter. The parameterizations and the lumped nature of this viscoelastic model enable straightforward parameter exploration, but would be made unnecessary by a fully 3-D viscoelastic ice-stream model with realistic sub-shelf bathymetry. However, the bathymetry and many other parameters that would be required for such a model are generally poorly constrained. We therefore leave such a modeling approach to future studies (Table 2).
MODELED TIDAL ICE FLOW VARIATIONS AT RUTFORD AND BINDSCHADLER ICE STREAMS
The lumped Maxwell model we have described can explain why glaciers respond differently to tides. In this section, we focus on the two cases of RIS and BIS where the phase, amplitude and period of tidally-driven ice flow variations are considerably different. We use the lumped Maxwell model to calculate the detrended displacement calculated for parameter regimes similar to these two cases. We have plotted the results in Figs. 1e and f for easy comparison with observations (Figs. 1a and b).
In a Rutford-like case (Fig. 1e, with α = 1.54, β = 1.4 and γ = 2.9), the ice shelf is forced by tidal variation at M 2 and S 2 periods, with hydrostatic stress variations of 55 kPa amplitude (consistent with combined M 2 and S 2 tidal height variations of ~6 m). In this case, β > 1 and γ is large, corresponding to a situation where changes in hydrostatic stresses are smaller than changes in both buttressing stresses and background stresses in the ice shelf. Variations in simulated displacement are nearly in phase with semidiurnal variations in tidal height. These semidiurnal variations in displacement are smaller in amplitude than the dominant low-frequency signal at the M sf tidal period, which arises from the buttressing nonlinearity (as α > 1 in this case). We can understand this low-frequency ice-stream response by noting that a nonlinear power of the tidal height as it appears in Eqn (5) can be expanded as
which can be rewritten using the double angle and product to sum of sinusoidal functions
which includes power at a low beat frequency ω L = ω h1 − ω h2. The ‘frequency mixing’ of Eqn (13) is a mathematical property of any system with nonlinearities, which is forced at more than one frequency. In this case, the nonlinearity arises through processes occurring in the ice shelf. By providing a nonlinear mechanism that produces a low-frequency ice flow response in the ice shelf, we can explain InSAR observations that show that the low-frequency signal at RIS originates in the ice shelf (Minchew and others, Reference Minchew, Simons, Riel and Milillo2017). More broadly, with these parameters, the model quantitatively reproduces the characteristic phase, amplitude and frequency of observed flow variations in the ice shelf downstream of RIS (Fig. 1a).
In a Bindschadler-like case (Fig. 1f, with α = 2.53, β = 0.3 and γ = 4), the ice shelf is forced by tidal variation at O 1 (25.92 h) and K 1 (24.00 h) periods, with hydrostatic stress variations of 30 kPa amplitude (consistent with combined O 1 and K 1 tidal height variations of ~3.2 meters). In this case, β < 1 and γ is large, indicating a scenario where changes in hydrostatic stresses are large compared with changes in buttressing stresses and small compared with background stresses in the ice shelf. Consequently, diurnal variations in horizontal displacement are close to 180° out of phase with the prescribed diurnal variations in tidal height. Although the dominant variations in hydrostatic stress are linear with tidal height, nonlinear changes in buttressing stresses (α > 1 also in this case) and effective viscosity are sufficiently important to introduce a significant low-frequency component into the tidal response that is comparable with the amplitude of the diurnal response. This case quantitatively reproduces the phase, amplitude and frequency of ice flow variations at BIS, which are close to 180° out of phase with tidal height and comparable in amplitude with the low-frequency response (Fig. 1b).
Comparison of the two cases of RIS and BIS shows that even though the dominant tidal periods in the Weddell and Ross Seas are quite different, the existence of an asymmetry in buttressing stresses with respect to tidal height produces an ice stream response that includes variation on long spring-neap timescales for both ice streams. However, there are also important differences in the flow response of BIS and RIS. In particular, there are significant differences in the phasing of the high-frequency ice shelf flow response relative to the dominant tidal forcing and the amplitude of the high-frequency ice shelf flow response relative to the low-frequency flow response.
PARAMETER EXPLORATION OF TIDAL INFLUENCE ON ICE SHELF FLOW
In this section, we systematically explore the parameter space of the lumped Maxwell model. We explore how the high-frequency phase and dominant period of the modeled variations in ice shelf displacement are functions of parameters α, β and γ. We show that it is the relative importance of tidal variations in hydrostatic and buttressing stresses (β in Eqn (5)), and to a lesser extent the strength of the viscous ice shelf response to tides (γ in Eqn (9)), that cause differences in the tidally-driven ice flow variation at different ice streams. We start by analyzing the phase of high-frequency variations in horizontal ice displacement with respect to tidal height, as a function of β and γ (Figs. 3a and b) for two values of α that are similar to those for RIS and BIS (see Figs. 1c and d). When β < 1, changes in hydrostatic stresses dominate changes in buttressing stresses and displacement is close to 180° out of phase with tidal height. Conversely, when β > 1, changes in buttressing stresses dominate and the displacement is nearly in phase with tidal height. There is a sharp transition between these two regimes at β = 1, suggesting that the observed phase lag between high-frequency variations in ice displacement and tidal height (and whether it is near 0° or 180°) is a strong indicator of whether hydrostatic stress variations or buttressing stress variations dominate the tidal modulation of ice shelf flow. Furthermore, increasing γ decreases the effective viscosity, therefore increasing the amplitude of the viscous deformation response to forcing, leading to an increase in the phase lag of horizontal displacement with respect to tidal forcing.
To determine the dominant timescale of the ice shelf response to tides, we examine (Figs. 3c and d) the difference between the logarithm of power spectral density (S) for detrended horizontal ice displacement at low frequencies (M sf) and the high frequencies of tidal forcing (M 2 and S 2 in Fig. 3), as a function of β and γ for the same two values of α as in Figs. 3a and b. Where values are large, the low-frequency response of ice flow to tidal forcing dominates over the high-frequency ice flow response (at the same frequency as the forcing). Changes in γ cause two opposite effects in the ice flow response to tidal forcing. When γ increases, the time-averaged deviatoric stresses in the nonlinear effective viscosity ( $\bar I_2 $ in Eqn (9)) become larger relative to the tidally-varying terms (σ h and σ b). The result is that tidal forcing causes proportionately less change in effective viscosity and the nonlinearity of effective viscosity in the Nye-Glen flow law for ice becomes less important. In addition, increasing γ decreases the overall magnitude of effective viscosity (through the dependence of viscosity on deviatoric stresses in ice, see Eqn (9)). When the ice shelf strain varies on periods longer than the Maxwell time (the characteristic relaxation timescale for a Maxwell material: τ m = η e/E), the corresponding deformation will primarily be viscous. At lower viscosity, the amplitude of the resulting variations in viscous deformation will be larger, while the high-frequency elastic deformations are largely unchanged (Lakes, Reference Lakes1998). Consequently, this second effect leads to an increase in the amplitude of low-frequency variations in ice flow relative to high-frequency variations with increasing γ. Figs. 3c and d show that in a range of parameters relevant to ice shelves, increasing γ generally leads to an increase in the amplitude of low-frequency components relative to high-frequency components. Thus, we conclude that the second effect is generally more important than the first, except when α and γ are low and β is very high (like in the bottom right corner of Fig. 3c). We also find that, even when the magnitude of linear changes in hydrostatic stresses dominate nonlinear changes in buttressing stresses (β < 1), low-frequency variations in ice flow can be significant (as long as γ is sufficiently large). Even in the Bindschadler-like case shown in Fig. 1f, there is a significant component of the variation in ice shelf flow occurring at M sf periods.
As argued in the section ‘Modeled tidal ice flow variations at Rutford and Bindschadler ice streams’, when α ≠ 1, the nonlinearity in buttressing stress with respect to tidal height causes mixing of high-frequency forcings to produce a low-frequency response in ice flow. Thus, increasing α leads to an expansion in the range of parameter space over which the low-frequency ice-shelf response dominates (in Fig. 3, white regions demarcate where this ratio is near one). In the marginal case where variations in buttressing stresses are exactly linear with respect to tidal height (α = 1; not shown), we find that no amount of viscous nonlinearity is able to produce power at low frequencies. This lack of low-frequency response is similar to previous studies (Gudmundsson, Reference Gudmundsson2007, Reference Gudmundsson2011), which found that the nonlinearity of the Nye-Glen flow law is insufficient on its own to reproduce the low-frequency response of RIS to high-frequency tidal forcing. Other nonlinearities, which come from asymmetries in buttressing in our model (and sliding relations for grounded ice in other studies) are necessary, though not strictly sufficient, to produce a dominantly low-frequency response in horizontal ice displacement. Given that γ is large in most ice streams, we expect there generally to be some detectable ice flow variation at the tidal beat frequency if tides cause variations in ice shelf buttressing stress with magnitude greater than ~25% of the magnitude of hydrostatic stress variations (with the specific number depending on precise values of γ and α).
DISCUSSION
A combination of tidal variations in hydrostatic and buttressing stresses can simultaneously explain the phase of high-frequency variations in ice-shelf flow, the origin of low-frequency variations in the ice shelf, and the likely cause of small-scale spatial variability in strain rate in the ice shelf. Other studies have proposed alternative mechanisms through which tides influence ice flow, including variations in ice shelf driving stress, basal friction from sub-shelf tidal currents, subglacial water pressure and nonlinear sliding laws. However, each of these mechanisms has difficulty in explaining particular aspects of existing observations. In one modeling study of the Ross Ice Shelf, Brunt and MacAyeal (Reference Brunt and MacAyeal2014) were unable to reproduce observed changes in flow through tidally-driven changes in ice shelf driving stress and sub-shelf ocean currents. In other modeling studies of RIS (Gudmundsson, Reference Gudmundsson2007, Reference Gudmundsson2011; Rosier and others, Reference Rosier, Gudmundsson and Green2014; Thompson and others, Reference Thompson, Simons and Tsai2014; Rosier and Gudmundsson, Reference Rosier and Gudmundsson2016), tidal modulation of ice flow through variations in subglacial water pressure or nonlinearity in the basal sliding law create a low-frequency signal first in the grounded ice stream, which is then passively propagated uniformly through the ice shelf. Such an M sf signal in the grounded ice stream would precede (in time) the M sf signal in the floating ice shelf and also cause uniform changes in ice shelf strain rates. However, Minchew and others (Reference Minchew, Simons, Riel and Milillo2017) have shown, using spatially-resolved InSAR observations from RIS, that the largest amplitude and earliest appearance of the low-frequency tidal signal in horizontal flow occurs in the ice shelf, not in the grounded ice stream. That same study showed that, in the ice shelf, strain rates do not vary uniformly in space over the M sf tidal cycle, but rather exhibit significant spatiotemporal variability, presumably due to spatial variations in contact with sub-shelf pinning points. Both InSAR and GPS observations, together with the viscoelastic model of this study, show that processes within the ice shelf are capable of producing variations in ice flow at the spring-neap beat frequency.
The tidal height variations in the sub-shelf ocean cavities of the Ross and Weddell Seas are not spatially uniform (Robertson and others, Reference Robertson, Padman and Egbert1998; Padman and others, Reference Padman, Erofeeva and Joughin2003). Additionally, other tidal constituents, such as the M f tide (13.66 days), also cause vertical deflection and horizontal flow variations in the Ross and Filchner-Ronne ice shelves. Further yet, there may also be far-field tidal variations at other ice streams, which cause interactions between the tidal responses of ice streams within an ice shelf. However, these other tidal constituents and interactions seem to be present at much lower amplitudes than the primary semidiurnal and diurnal tides (Padman and others, Reference Padman, Erofeeva and Joughin2003; Murray and others, Reference Murray, Smith, King and Weedon2007) and do not play a significant role in producing the largest amplitude tidal ice flow variations. To account for many of these detailed (but secondary) aspects of observations, future studies may include the processes highlighted in this study in conjunction with fully 3-D viscoelastic ice stream/shelf models with realistic tidal and bathymetric boundary conditions. However, such an approach is beyond the scope of this study.
CONCLUSIONS
We have shown that asymmetries in buttressing stress with respect to tidal height can produce low-frequency ice flow variations and, in combination with tidal variations in hydrostatic stresses, can also explain the phase of high-frequency ice flow variations. While previous studies have suggested a potential role for tidal variations in buttressing stress (Schmeltz and others, Reference Schmeltz, Rignot and MacAyeal2001; Doake and others, Reference Doake2002; Heinert and Riedel, Reference Heinert and Riedel2007; Thomas, Reference Thomas2007; Minchew and others, Reference Minchew, Simons, Riel and Milillo2017), we describe a simple viscoelastic model that can quantitatively reproduce the observed temporal variations in ice shelf horizontal displacement at different ice streams. We show that nonlinearities in ice flow dynamics with respect to tidal height cause frequency mixing and ice flow variations at low frequencies corresponding to the beat of the primary tidal frequencies. Viscoelasticity then acts to increase the amplitude of low-frequency ice flow variations relative to the amplitude of high-frequency ice flow variations. We show that the M 2 and S 2 tides produce an ice flow response primarily at their beat frequency (the M sf frequency).
A new generation of observations shows that the ice sheet response to ocean tides is more complex than previously thought. We have proposed new mechanisms of tidal influence on ice flow and a simple model to explain the wide range of tidally-modulated ice flow variability exhibited by glaciers, including the characteristics of tidally-driven ice flow variability at Rutford and Bindschadler Ice Streams. Pinning points beneath ice shelves are difficult to observe, especially when they only make ephemeral contact with the ice shelf. However, as the rate of basal ice shelf melting in West Antarctica increases due to warming of the sub-shelf ocean cavity (Pritchard and others, Reference Pritchard2012), these marginal pinning points are likely to be the most immediate locations of decreasing ice shelf buttressing stress. Thus, observations of tidally-modulated ice shelf flow can potentially reveal the previously hidden influence of pinning points on ice stream flow under changes to the lower ice shelf surface. In this fashion, this study lays the theoretical and conceptual groundwork for pushing further in our understanding of the ice sheet response to tides and future climate-driven changes to ice shelves.
ACKNOWLEDGMENTS
We thank Doug MacAyeal, Kelly Brunt, Hilmar Gudmundsson and Olga Sergienko for helpful conversations during the completion of this work. We also thank Bob Bindschadler and Hilmar Gudmundsson for providing GPS data. Source code and documentation of the models used in this study are available as public repositories on GitHub: https://github.com/aarobel/. AR has been supported by the NOAA Climate & Global Change and Caltech Stanback postdoctoral fellowships during the completion of this work. BM was supported by a NSF Earth Sciences Postdoctoral Fellowship award EAR-1452587. This work was also supported by NSF grant EAR-1453263.