the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Hydraulic suppression of basal glacier melt in sill fjords
Eef van Dongen
Martin Jakobsson
Matt O'Regan
Christian Stranne
Using a conceptual model, we examine how hydraulically controlled exchange flows in silled fjords affect the relationship between the basal glacier melt and the features of warm intermediate Atlantic Water (AW) outside the fjords. We show that an exchange flow can be forced to transit into the hydraulic regime if the AW interface height decreases, the AW temperature increases, or the production of glacially modified water is boosted by subglacial discharge. In the hydraulic regime, the heat transport across the sill becomes a ratelimiting factor for the basal melt, which is suppressed. An interplay between processes near the ice–ocean boundary and the hydraulically controlled exchange flow determines the melt dynamics, and the sensitivity of the basal melt to changes in the AW temperature is reduced. The model results are discussed in relation to observations from the Petermann, Ryder, and 79^{∘} N glaciers in northern Greenland.
In response to global warming, the Greenland Ice Sheet (GIS) has lost mass over the past decades and its marine outlet glaciers are retreating (Mouginot et al., 2019; Straneo and Heimbach, 2013; Wood et al., 2021). The GIS holds an ice volume equivalent to 7.4 m of sea level (Morlighem et al., 2017), and it may contribute up to 0.3 m global mean sealevel rise by the end of this century (Aschwanden et al., 2019). Over half of the recent mass loss from the GIS is from increased discharge of ice into the ocean from marine outlet glaciers (Mouginot et al., 2019), where calving and oceanic melt of marine ice have increased. Mass loss from marineterminating glaciers can cause a positive feedback: resistive stresses in grounded or floating marine glaciers buttress ice inland, and icestream flow speed and ice export across the grounding line can increase when marine glaciers retreat (Schoof, 2007; Gudmundsson, 2013; Nick et al., 2013; Schoof et al., 2017). This accelerated ice loss contributes directly to sealevel increase.
The subsurface melt on marine glaciers is primarily controlled by the excess ocean temperature over the (pressuredependent) freezing temperature at the grounding line (Holland and Jenkins, 1999), the point where the ice begins to float (or for tidewater glaciers, the water depth at their essentially vertical fronts). The melt depends also on factors such as basal slope, subglacial discharge, tidal currents, and water column stratification (Jenkins, 2011; Truffer and Motyka, 2016; De Andrés et al., 2020). The ice melt mixes with ocean water, which creates a buoyant meltwater plume that rises along the base of the ice tongue (Lewis and Perkin, 1986). Turbulence in the meltwater plume transports heat to the ice–ocean boundary and sustains melt in the rising plume. Marine glaciers in Greenland terminate in fjords, and basal melt is chiefly driven by heat supplied in subsurface Atlantic Water (AW) that enters the fjords (Straneo et al., 2012). In Greenland, basal melt is sensitive to the AW temperature (Straneo and Heimbach, 2013), and increases in AW temperature and subglacial discharge have been the major drivers of the retreat of outlet glaciers in deep Greenlandic fjords since the mid 1990s (Wood et al., 2021; Slater and Straneo, 2022). However, local features, such as fjord geometry and wind conditions, affect the sensitivity of the basal melt to changes in the AW temperature in the open ocean (Straneo and Cenedese, 2015; Khazendar et al., 2019; Wood et al., 2021).
The present study is motivated by recent observations of hydraulically controlled exchange flows at sills in the Greenlandic fjords that host the ice tongues of the Ryder and 79^{∘} N glaciers (Jakobsson et al., 2020; Schaffer et al., 2020). The hydraulic control sets an upper limit on the exchange flow that depends on sill geometry and upstream stratification (Pratt and Whitehead, 2007). Accordingly, hydraulic control limits the heat transport that sustains the basal melt and has the potential to stabilize marine glaciers. Numerous observations of sill flows demonstrate that the vertical mixing increases strongly when the flow becomes hydraulically controlled (Pratt and Whitehead, 2007), and Jakobsson et al. (2020) and Schaffer et al. (2020) show that as inflowing AW passes over the sills and descends on the landward slopes, it mixes with overlying cold, glacially modified water. As a result, the waters reaching these glaciers' grounding lines are colder than the AW outside the fjords. This reduces the basal melt compared to the case when unmodified AW reaches the glacier. Thus, hydraulic control can reduce basal melt by limiting the exchange flow as well as by decreasing the water temperature at the grounding line of the glacier.
The feature that hydraulic control (and/or fjord geometry) can limit the heat transport suggests that there are two different regimes of icetongue basal melt in fjords. First, one where the basal melt is controlled locally by turbulent processes near the ice–ocean boundary, which determines the heat flux from the fjord water to the ice (Fig. 1a). In this scenario, the fjordscale circulation adjusts to deliver the heat required for the basal melt, and AW reaches the grounding line of the ice tongue. Second, in sill fjords hydraulic control may be established, which constrains the exchange circulation in the fjord and its associated heat transport (Fig. 1b). In this case, the basal melt is not solely controlled by local processes near the ice–ocean boundary, as the fjordscale heat transport towards the glacier enters as a ratelimiting factor. In the relatively narrow Greenlandic fjords, sill geometry is likely to be a major factor constraining ocean heat transport towards marine glaciers (Zhao et al., 2021; Bao and Moffat, 2023). For large Antarctic ice shelves, effects due to Earth's rotation are important, and the oceanic heat flux available for basal melt may be controlled by mesoscale ocean eddies or largescale flows constrained by conservation of potential vorticity (Little et al., 2008; Hattermann et al., 2014; Zhao et al., 2019).
The observations from the Ryder and 79^{∘} N glaciers (Jakobsson et al., 2020; Schaffer et al., 2020) raise the question of how strongly hydraulic control limits basal melt and how it affects the dependence of basal melt to the temperature and height of the AW layer outside the fjords. Here, we examine this question using a conceptual twolayer fjord model that includes ocean–glacier interactions. The model results are discussed in relation to observations from the Greenlandic ice tongues of the Petermann, Ryder, and 79^{∘} N glaciers. However, with some modifications the model can be applied also to fjords with tidewater glaciers. Before the model is presented, we give a brief overview of the oceanographic conditions at the Petermann and Ryder glaciers.
The model result will be discussed in relation to the ice tongues of the 79^{∘} N, Petermann, and Ryder glaciers. These glaciers have the three largest ice tongues in Greenland (Wilson et al., 2017; Hill et al., 2018) and are located in the northern part of the island (Fig. 2). The geometries of the fjords in which these glaciers drain have some general features in common, including relatively large sill depths: about 500 m for 79^{∘} N and about 400 m for the Petermann and Ryder glaciers. Here, we will describe fjord geometries and oceanographic conditions for Petermann and Ryder, which are located relatively close to each other (∼ 200 km apart) and drain into fjords that terminate in the Lincoln Sea. The oceanographic conditions in the fjord of the 79^{∘} N ice tongue, which is Greenland's largest and about 80 km long, are described by, for example, Lindeman et al. (2020) and Schaffer et al. (2020).
Figures 2 and 3 show bathymetric and temperature conditions in the Sherard Osborn and Petermann fjords, where the Ryder and Petermann glaciers drain. In Petermann Fjord, which has a ∼ 400 m deep and ∼ 12 km wide sill, AW with similar features is encountered inside as well as outside the fjord, and there are no indications of hydraulic control at the sill (Johnson et al., 2011; Jakobsson et al., 2020).
Sherard Osborn Fjord has a more constrictive fjord topography, with an outer and an inner sill. The temperature in the AW depth range decreases across the sills, with the coldest temperature in the fjord basin landward of the inner sill that is largely capped by the ice tongue. The largest temperature drop occurs over the inner sill. Here, a strong nearbottom inflow was observed, occurring in a ∼ 400 m deep and ∼ 1 km wide channel on the eastern sill, demonstrating that the inflow is hydraulically controlled (see Fig. 4 in Jakobsson et al., 2020). Accordingly, the inner sill provides the main geometrical constraint on the exchange flow and heat transport to the glacier.
During the last 6 decades the ice tongues of Petermann and Ryder have evolved differently: the former retreating significantly (∼ 300 m yr^{−1}) and the latter advancing modestly (∼ 40 m yr^{−1}) (Hill et al., 2018), and in 2010 and 2012 Petermann lost ∼ 35 km of its ice tongue in major calving events (Johannessen et al., 2013). The two glaciers are relatively closely located and are expected to experience similar atmospheric conditions and to have AW of similar temperatures outside the fjords. Jakobsson et al. (2020) proposed that the differences in fjord bathymetry are the reason for different behavior of the two glaciers: Ryder Glacier has a more restrictive sill geometry, partly protecting the ice tongue from the inflow of warmer subsurface AW (Fig. 2).
In summary, these results suggest tentatively that the basal melt on Petermann is chiefly rate limited by processes near the ice–ocean boundary, whereas the basal melt on Ryder is partly rate limited by largescale heat transport towards the glacier. In addition to differences in sill geometry, differences in summer seaice conditions influence the efficiency in transporting AW into the two fjords. Perennial landfast sea ice outside Sherard Osborne Fjord curtails winddriven water exchange between the fjord and the open ocean, whereas icefree conditions in and around Petermann Fjord allow for a more vigorous winddriven water exchange during summer (Shroyer et al., 2017; Jackson et al., 2018; Stranne et al., 2021). We will now go on to describe a twolayer model that will be used to examine the interplay between basal melt dynamics and hydraulic control.
We consider a twolayer model of glacier–ocean interaction in a fjord, with AW (T_{A}, S_{A}) and glacially modified water (T, S) (see Straneo and Cenedese, 2015; Jackson and Straneo, 2016, for background). Figure 1 shows the model geometry for two different circulation regimes that will be examined. The model represents near steadystate conditions, and we assume that the timemean exchange flow in the fjord is primarily driven by basal melting of the ice tongue, which creates a buoyant meltwater plume rising along its base. Higher up, the plume becomes neutrally buoyant and feeds the outflow of glacially modified water. A fresh, lowdensity layer of Polar Surface Water caps the two water masses represented in the model, insulating them from surface runoff and contact with sea ice and the atmosphere. The Polar Surface Water is not explicitly represented in the model.
Although subglacial discharge can have a strong impact on subsurface melt rates, we will for simplicity neglect subglacial discharge in the model's freshwater budget. The reason is twofold. First, the resulting model becomes simpler and more tractable analytically. In Appendix A, we describe a more complex model version that includes subglacial discharge in the conservation relations: this shows that the results remain qualitatively similar even when the subglacial discharge is significantly greater than the subsurface melt. Second, observations indicate that freshwater input due to basal melt exceeds subglacial discharge for large ice tongues such as 79^{∘} N, Ryder, and Petermann: Schaffer et al. (2020) estimated that in the annual mean the subglacial discharge constitutes only about 10 % of the freshwater exported from 79^{∘} N Glacier, and the summer measurements from Petermann of Washam et al. (2019) indicate that the freshwater fraction due to subglacial discharge in the glacially modified water column below the ice tongue is less than 30 % (see their Fig. 5). This suggests that for large ice tongues, subglacial discharge may, as a leading order approximation, be neglected in the model's freshwater budget, but subglacial discharge will be allowed to affect the model's melt rates.
3.1 Conservation relations
In the twolayer model of the fjord, the melting of the ice tongue is the only local source/sink of freshwater/heat. This can be used to formulate conservation relations for volume, salt, and heat (Jackson and Straneo, 2016; Truffer and Motyka, 2016). At the sill, conservation of volume is given by
where Q and Q_{A} are the volume outflow and inflow of glacially modified and Atlantic waters, respectively, and M the freshwater input due to melting. The ice consists of pure freshwater, implying that the meltwater input M does not affect the salinity budget. Hence conservation of salt is given as
Combining Eqs. (1) and (2) yields Knudsen's relation for salt conservation:
The heat budget involves a balance between advective heat transport towards the glacier and basal melt. The advective heat flux is
where c is the heat capacity (per unit volume) of seawater and the melt freshwater input M is assumed to have the salinitydependent freezing temperature T_{f}. Note that T_{f} will be taken as a constant set by the groundingline pressure and S_{A}. Using Eqs. (1) and (4), we obtain
The heat flux is related to the ice melt (M) as
where L is the latent heat of freezing, c_{i} the heat capacity of ice, and T_{i} the ice temperature.
By combining Eqs. (5) and (6), we obtain
Here, $L/c\approx \mathrm{75}$ ^{∘}C and in northern Greenlandic fjords T_{A}−T_{f} is typically 3 ^{∘}C. Hence, the terms involving T_{A}−T_{f} can to a good approximation be neglected in Eqs. (5) and (7). It is convenient to define a generalized Gade temperature (Gade, 1979):
which would be the decrease in temperature of a unit volume of water from which sensible heat is extracted to melt ice corresponding to a unit volume of liquid water.^{1} By using the definition of T_{G}, Eq. (7) can be written as
which gives the relation between the advective heat flux and the melt. Note that unless the ice temperature is very cold, ${T}_{\mathrm{G}}\approx L/c$.
The conservation relations of the model can be summarized as follows.

Volume. The meltwater input M is small, implying that Q≈Q_{A}. Thus, the inflow of AW approximately equals the volume outflow of glacially modified water. In what follows, we will denote the exchange flow simply by Q.

Salt. The salt balance is given by Eq. (3). Here, M cannot be neglected since it is multiplied by S_{A}, which is larger than ΔS: Eq. (3) states that $M/Q=\mathrm{\Delta}S/{S}_{\mathrm{A}}$.

Heat. The heat budget is specified by Eq. (9), which in combination with Eq. (3) yields
$$\begin{array}{}\text{(10)}& {\displaystyle \frac{\mathrm{\Delta}S}{{S}_{\mathrm{A}}}}={\displaystyle \frac{\mathrm{\Delta}T}{{T}_{\mathrm{G}}}}.\end{array}$$For the melting of ice in seawater, heat and salt conservation yields a linear relationship between the salinity and temperature differences (Gade, 1979). This allows us to either use ΔT or ΔS in our analyses; we will use ΔT.
The layer density difference is calculated using a linear equation of state
where ρ_{0} is a constant seawater reference density and where $\mathit{\alpha}=\mathrm{4}\times {\mathrm{10}}^{\mathrm{5}}$ K^{−1} and $\mathit{\beta}=\mathrm{8}\times {\mathrm{10}}^{\mathrm{4}}$ are the thermal and haline expansion coefficients, respectively. Equation (10) allows the density difference to be written as
Here, $\left(\mathit{\alpha}{T}_{\mathrm{G}}\right)/\left(\mathit{\beta}{S}_{\mathrm{A}}\right)\approx \mathrm{0.1}$ showing that the salinity dominates the density difference.
3.2 Basal melt parameterization
We will use a parametrization of the basal melt (M), which depends on the difference between the ocean water temperature (T_{C}) and the freezing point temperature of seawater (T_{f}) at the grounding line (Holland et al., 2008; Jenkins, 2011; Xu et al., 2013; Favier et al., 2019). We denote the thermal forcing as
If AW reaches the grounding line, then T_{C}=T_{A}, and the thermal forcing is denoted as
However, mixing between in and outflowing waters over a sill can lower T_{C} relative to T_{A}, which implies that 𝒯 can be lower than 𝒯_{A}.
We model the areaintegrated basal melt as (Xu et al., 2013)
where γ_{1} is a coefficient that depends on features such as the icetongue geometry and subglacial discharge and n_{1}>0 is an exponent. We assume that the plume volume flow Q_{P} is also a function of the thermal forcing 𝒯 and given by
where γ_{2} is a constant and n_{2}>0 an exponent.
Several studies have used models of varying complexity to examine the relationship between thermal forcing and the areaintegrated melt on Greenlandic ice tongues and Antarctic ice shelves (e.g., Jenkins, 1991, 2011; Little et al., 2009; Lazeroms et al., 2018, 2019; Holland et al., 2008; Cai et al., 2017; Favier et al., 2019). These investigations report values of n_{1} in the range from 1 to 2, with a preference for n_{1} values of around 1.5 (Xu et al., 2013; Cai et al., 2017) to 2.0 (Holland et al., 2008; Little et al., 2009). The reported range of n_{1} is likely to reflect both different ice–ocean interaction regimes and model assumptions on the boundary conditions at the ice–ocean boundary. There are fewer studies that specifically comment on the relationship between the volume transport in the plume and the thermal forcing, but Holland et al. (2008) reported a linear dependence of Q_{P} on the thermal forcing; i.e., n_{2}≈1.
Primarily, the heat flux to the ice and the associated melting depend on the product of the thermal forcing and the plume velocity (say u), which in turn is related to the plume buoyancy (Holland and Jenkins, 1999; Favier et al., 2019). If the plume buoyancy is proportional to 𝒯 and the buoyancy force is balanced by a linear basal friction, then u∝𝒯. This gives n_{1}=2, corresponding to a quadratic relation between melt and thermal forcing (Holland et al., 2008; Little et al., 2009). If the basal friction is quadratic, i.e., proportional to u^{2}, on the other hand, scaling analyses suggest that $u\propto {\mathcal{T}}^{\mathrm{1}/\mathrm{2}}$ (Lazeroms et al., 2018), which gives n_{1}=1.5. Jakobsson et al. (2020) applied the plume model of Jenkins (1991) to the Ryder ice tongue, and their results suggest that n_{1}≈1.7 and n_{2}≈0.7. Notably, if M∝u𝒯 and Q_{P}∝u, then $M/{Q}_{\mathrm{P}}\propto \mathcal{T}$. In view of Eqs. (15) and (16) this implies that
This constraint on the exponents leads to some attractive simplifications of the dynamics, which will be used in the analyses.
It is worth noting that theoretical considerations (e.g., Straneo and Cenedese, 2015; Jenkins, 2011, and references therein) indicate that on marine glaciers where the buoyancy source is dominated by the subglacial discharge near the grounding line rather than by the distributed melt along the submerged glacier, the melt is approximately proportional to the thermal forcing and the plume volume transport is essentially independent of the thermal forcing. This limit of high subglacial discharge, characteristic of summer conditions at Greenlandic tidewater glaciers (e.g., Straneo and Cenedese, 2015), is described by the case of n_{1}=1 and n_{2}=0, which satisfies the constraint of Eq. (17). Some aspects of the case with high subglacial discharge are discussed in Appendix A.
In summary, the literature reports a range of values for n_{1} and n_{2}. However, n_{1}=2 and n_{2}=1 appear as one reasonable choice for the exponents for qualitatively examining the dynamics of large Greenlandic ice tongues such as the 79^{∘} N, Petermann, and Ryder glaciers. This will be the baseline case when we examine the interplay between melt dynamics and hydraulic control in Sect. 3. In Sect. 4.2.5, we will consider how variations in the values of n_{1} and n_{2} affect the results, including the case of n_{1}=1 and n_{2}=0. When we derive general results below, however, we will allow n_{1} and n_{2} to be arbitrary positive numbers but subject to the constraint n_{1}>n_{2}.
3.3 The meltcontrolled exchange flow regime
Consider a situation in which the fjord geometry, via frictional resistance or hydraulic control, does not limit the exchange flow and its associated heat transport towards the ice tongue. We assume that unmodified AW reaches the glacier (T_{C}=T_{A}), and the melt processes create a plume volume flow Q_{P} that sets the exchange flow: if T_{A} increases, the exchange flow also increases at the rate given by Eq. (16). This regime, in which the strength of the exchange flow is controlled locally by the glacier basal melt, will be referred to as the meltcontrolled regime, and a contrasting hydraulically controlled exchange flow regime will be presented in Sect. 3.4.
In the meltcontrolled regime, where Q=Q_{P}, we can use the heat conservation relation (Eq. 9) together with Eqs. (15) and (16) to obtain
This relationship, which is equal to $M/Q$, shows that ΔT as well as the ratio $M/Q$ increase with 𝒯; i.e., the meltwater fraction in the plume increases with thermal forcing. Note that the condition ${n}_{\mathrm{1}}{n}_{\mathrm{2}}=\mathrm{1}$ yields a linear relation between ΔT as well as $M/Q$ and the thermal forcing.
By dividing Eq. (18) with $\mathcal{T}/{T}_{\mathrm{G}}$, we obtain
Since 𝒯=𝒯_{A} in the meltcontrolled regime, the lefthand side in this expression equals $({T}_{\mathrm{A}}T)/({T}_{\mathrm{A}}{T}_{\mathrm{f}})$, which is less or equal to 1 since T≥T_{f}. When ${n}_{\mathrm{1}}{n}_{\mathrm{2}}=\mathrm{1}$, the righthand side becomes independent of 𝒯 and equals
This parameter is a nondimensional measure of the temperature of the outflowing glacially modified water (T): when σ=1, T=T_{f}, and when σ=0, T=T_{A}. (This interpretation of the Eq. (19) applies also when ${n}_{\mathrm{1}}{n}_{\mathrm{2}}\ne \mathrm{1}$, but then $\mathrm{\Delta}T/\mathcal{T}$ is no longer constant in the meltcontrolled regime.) As will be shown below, σ influences aspects of the hydraulically controlled regime.
To summarize, the flow in the meltcontrolled regime is specified by a knowledge of the AW properties T_{A} and S_{A}, which determine the thermal forcing 𝒯=𝒯_{A}. In turn, this yields M, Q, and ΔT (Eqs. 15, 16, 18), and ΔS is obtained from Eq. (10).
We will now go on to examine how a hydraulically controlled exchange flow affects the melt dynamics. For this purpose, it is useful to write Q_{P} as a function of the temperature difference. By using Eq. (18), we obtain
Since n_{1}>0 and n_{2}>0, this shows that the meltcontrolled exchange flow increases with ΔT when ${n}_{\mathrm{1}}{n}_{\mathrm{2}}>\mathrm{1}$.
3.4 Hydraulic control
Fjord and sill geometries may impose limits on the exchange flow, which in turn can potentially alter the basal melt dynamics. In particular, hydraulic control of a twolayer exchange flow over a sill sets an upper bound for the exchange flow (say Q_{H}), which is determined by the upstream height of the AW layer above the sill (h) and the layer density difference (Pratt and Whitehead, 2007; Zhao et al., 2021). Exchange flow strengths below the critical value Q_{H} are unconstrained by the geometry and are referred to as subcritical flows. Thus, it is conceivable that a sufficiently strong meltdriven exchange flow or a high sill can cause a transition from a subcritical flow to a critical, hydraulically controlled flow (Pratt and Whitehead, 2007).
How the flow evolves as the meltdriven exchange flow (or the sill height) is gradually increased and approaches the hydraulically controlled limit is complex and depends on fjord and sill geometry (Armi, 1986; Pratt and Whitehead, 2007; Nycander et al., 2008). Observations from the Ryder and 79^{∘} N glaciers show that the inflow at the sills in front of the ice tongues are hydraulically controlled and that the thickness of the inflow layer is thin compared to the upper outflowing layer (Jakobsson et al., 2020; Schaffer et al., 2020). This implies that the flow can be approximated by a onelayer hydraulic model representing the inflowing AW layer. Two additional features allow for simplifications of the hydraulic model. First, the depth of the AW layer on the seaward side (upstream) of the sill is much larger than at the sill, which implies that the upstream AW inflow velocity is negligible. Second, the inflow over the sill is confined in a channel that is small compared to the internal Rossby radius, which allows the Earth's rotation to be neglected. In this situation, the maximum hydraulically controlled volume flow is given by (Pratt and Whitehead, 2007)
where W is the crosssectional width of the lower layer on the sill,^{2} h the height of the lower layer above the sill upstream (Fig. 1), and g the gravitational acceleration. By using Eq. (12), Q_{H} can be written as
where we have introduced
In the hydraulically controlled regime, the exchange flow is given by Eq. (23); i.e., Q=Q_{H}. By using this in the heat conservation relation Eq. (9), we obtain
3.5 Steadystate regimes: meltcontrolled and hydraulically controlled exchange flows
The results presented above suggest that there can exist two different flow regimes in a fjord with basal icetongue melting: one where the exchange flow Q is determined locally by the basal melt processes on the glacier and one where it is hydraulically controlled. The two regimes have the following characteristics.

In the meltcontrolled regime, the volume flow of the meltwater plume is smaller than the upper limit set by hydraulic control; i.e., Q_{P}<Q_{H}. Accordingly, the sill does not constrict the exchange flow, which is specified by Eq. (16) (or Eq. 21). The flow at the sill is subcritical (Pratt and Whitehead, 2007), and as result there is limited mixing between inflowing and outflowing waters. Essentially unmodified AW reaches the grounding line of the ice tongue (Fig. 1a), which implies that the thermal forcing is given by 𝒯=𝒯_{A}. Since the freezing temperature (T_{f}) is set by the groundingline depth (and S_{A}, which to a good approximation can be taken as constant here), the thermal forcing at a specific glacier is externally determined by the AW temperature (T_{A}).

In the hydraulically controlled regime, the plume volume flow exceeds the hydraulic limit; i.e., Q_{P}>Q_{H}. The exchange flow is now determined by Eq. (23), and the flow at the sill crest is critical, and it accelerates down the landward slope of the sill (Pratt and Whitehead, 2007). Here the flow becomes supercritical, and inflowing AW mixes with colder outflowing water (Price and O'Neil Baringer, 1994; Pratt and Whitehead, 2007; Jakobsson et al., 2020; Schaffer et al., 2020). This lowers the temperature of the water reaching the grounding line (Fig. 1b); i.e., 𝒯<𝒯_{A}. Importantly, this implies that the thermal forcing is no longer directly set by 𝒯_{A}: 𝒯_{A} is the external forcing, but the local thermal forcing 𝒯 is determined by dynamics in the fjord. The relationship between 𝒯 and 𝒯_{A} can be expressed as
$$\begin{array}{}\text{(26)}& \mathcal{T}=R\cdot {\mathcal{T}}_{\mathrm{A}},\end{array}$$where $R=R({\mathcal{T}}_{\mathrm{A}},h)$ is a reduction factor which arises when the exchange flow is hydraulically controlled. Note that R<1 in the hydraulic regime; in the meltcontrolled regime R=1.
At the transition between the two flow regimes 𝒯=𝒯_{A} and Q_{P}=Q_{H}, which implies that the transports given by Eqs. (21) and (23) are equal. By using this and Eq. (18), which applies at the transition, after some rearrangements we obtain
This represents the condition ${Q}_{\mathrm{P}}/{Q}_{\mathrm{H}}=\mathrm{1}$ and gives the relationship between 𝒯 and h at the regime transition. The height of the AW above the sill at the transition (say h_{L}) as a function of 𝒯_{A} is given by
When h>h_{L}, the flow is in the meltcontrolled regime, and when h<h_{L}, the exchange flow becomes hydraulically controlled. Alternatively, Eq. (27) gives the thermal forcing at the regime transition as a function of h
Using this and Eq. (18), the temperature difference at the regime transition can be written as
If ${n}_{\mathrm{1}}{n}_{\mathrm{2}}=\mathrm{1}$, then ΔT_{L}=σ𝒯_{L}; see Eq. (20).
Lowering the AW height always brings the flow towards the hydraulic regime. If h is fixed, then 𝒯_{L} is also fixed, which implies that changes in the AW temperature can cause a transition between the meltcontrolled and the hydraulically controlled regimes. As illustrated in Fig. 4, the nature of the transition depends on the value of the exponent 3n_{2}−n_{1} in Eq. (27). If $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}>\mathrm{0}$, increasing 𝒯_{A} will bring the flow towards the hydraulically controlled regime. On the other hand if $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}<\mathrm{0}$, decreasing 𝒯_{A} will brings the flow towards the hydraulically controlled regime. This behavior follows from the fact that ${Q}_{\mathrm{P}}\propto \mathrm{\Delta}{T}^{\frac{{n}_{\mathrm{2}}}{{n}_{\mathrm{1}}{n}_{\mathrm{2}}}}$ and ${Q}_{\mathrm{H}}\propto \mathrm{\Delta}{T}^{\mathrm{1}/\mathrm{2}}$. In the case of $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}=\mathrm{0}$, the regime transition height h_{L} is independent of 𝒯_{A}. This is because Q_{P} and Q_{H} then have the same dependence on ΔT.
Figure 5 illustrates the relation between exchange flow and temperature difference (Eqs. 21 and 23) for the case of n_{1}=2 and n_{2}=1, which implies that Q_{P}∝ΔT. Here, the flow is in the meltcontrolled regime if ΔT<ΔT_{L} or from Eq. (18) equivalently if 𝒯_{A}<𝒯_{L}. By increasing ΔT, the flow increases and is shifted towards and into hydraulic control.
In the case where $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}<\mathrm{0}$, hydraulic control – for a fixed h – occurs for weak thermal forcing and exchange flow. In a real fjord, additional exchange flows driven by winds and tides may be larger than a modelpredicted weak hydraulic flow (Jackson and Straneo, 2016). This can prevent the establishment of hydraulic control and in effect yield an exchange flow in the meltcontrolled regime. We will briefly discuss this case in Sect. 4.2.5 .
It is possible that transitions between meltcontrolled and hydraulically controlled regimes can also be caused by seasonal variations in subglacial discharge, even if the AW features remain unchanged. The reason is that the basal melt M increases with the subglacial discharge (Jenkins, 2011; Xu et al., 2013), which has a pronounced seasonal cycle that tracks the surface melt on the glaciers (e.g., Truffer and Motyka, 2016; Cai et al., 2017; Slater and Straneo, 2022). Thus, it is possible that in some fjords the exchange circulation can be stronger and hydraulically controlled in summer when the subglacial discharge peaks and weaker and meltcontrolled in winter. The details of such seasonal regime transitions will depend, among other things, on how the ratio between exchange flow and melt ($M/Q$) depends on the subglacial discharge. However, we will not pursue this topic further.
In the hydraulic regime, the volume flow of the melt plume is predicted to exceed the exchange flow at the sill (Fig. 5). This would cause an imbalance in the production and export of glacially modified waters in the fjord, preventing a steady state to be established. Therefore, changes in the flow are expected when hydraulic control is established.
To examine some of the new features emerging when the flow becomes hydraulically controlled, it is instructive to consider a thought experiment in which the sill height is suddenly increased and hydraulic control is established. Initially, the production of glacially modified water will be larger than the exchange flow across the sill. As a result, the layer of glacial water inside the sill will thicken and possibly extend below the sill crest (Fig. 1). This has two important consequences for the basal melt. First, the inflowing AW will entrain glacial water, causing the temperature of water reaching the grounding line to decrease. Second, the meltwater plume will rise partly through ambient waters that are colder and lighter than the displaced AW. This reduces the buoyancy and speed of the plume, which will now also entrain colder water. Theses changes in the stratification in the ice cavity act to reduce the basal melt. Thus, we expect that the temperature and salinity distributions inside the sill evolve such that a new steady state, compatible with the hydraulically constrained exchange flow, is established.
The reasoning above suggests that, in the hydraulic regime, the interface height of (pure or modified) AW is no longer the same on the seaward and landward side of the sill. Thus additional variables, such as a fjord interface height, may be needed to model flow features and melt features in the hydraulic regime. However, we will not introduce additional model variables. Instead, we consider two idealized scenarios for how the interplay between melt dynamics and hydraulic control can determine the steadystate flow. In these scenarios, the features of the fjord stratification (i.e., interface height and layer difference in temperature and salinity) can be viewed as hidden model variables that influence the flow. In scenario 1, we implicitly assume that no glacially modified water is entrained into the inflowing AW near the sill. This is an extreme and less likely scenario as observations and modeling show that entrainment generally occurs (Schaffer et al., 2020; Jakobsson et al., 2020; Hager et al., 2022; Bao and Moffat, 2023). In scenario 2, on the other hand, entrainment plays a key role for closing the volume budget in the ice cavity.
4.1 Scenario 1: hydraulically constrained plume volume flow
4.1.1 Physical assumptions
Here, we assume that the steadystate flow and stratification inside the sill adjust such that

the plume volume flow (Eq. 16) and the exchange flow (Eq. 23) are equal, which implies that ${Q}_{\mathrm{P}}\left(\mathcal{T}\right)={Q}_{\mathrm{H}}(h,\mathrm{\Delta}T)$;

the relationship between basal melt M and thermal forcing (Eq. 15) still applies and equals the formula for M in the hydraulically controlled regime (Eq. 25); this gives a relationship of the form $\mathcal{T}=\mathcal{T}(h,\mathrm{\Delta}T)$.
From these assumptions, we obtain the following expressions for the thermal forcing and temperature difference:
Note that 𝒯 and ΔT both depend on the features of the melt representation and the hydraulic constant k_{H}. In the reference case (n_{1}=2 and n_{2}=1) the exponents in the expressions above simplify, and the hydraulic exchange flow (Eq. 23) becomes
Notably, the flow is independent of the AW temperature T_{A}: the strength of the hydraulically controlled flow is determined by h and the parameters γ_{1}, γ_{2}, and k_{H}, which control ΔT that is proportional to layer density difference.
4.1.2 Dynamical features
Figure 6 shows the dependence of the thermal forcing and the reduction factor on the AW forcing (𝒯_{A} and h) in the case of n_{1} and n_{2}. The flow is in the hydraulic regime when 𝒯_{A}>𝒯_{L} and h<h_{L}, and the opposite applies in the meltcontrolled regime. In the meltcontrolled regime, 𝒯 is equal to 𝒯_{A} and is independent of h. An increase in 𝒯_{A} or a decrease in h brings the flow towards the hydraulically controlled regime. In this regime, the flow features become (in this scenario) independent of the AW temperature. As a result, the R factor decreases with increasing 𝒯_{A} at a fixed h. This provides negative feedback on the basal melt. The flow features are sensitive to changes in h: when n_{1}=2 and n_{2}=1, one finds that 𝒯∝h^{3} and M∝h^{6}. Accordingly, the basal melt drops sharply with decreasing AW height.
In this hypothetical scenario, one should view 𝒯 as an effective thermal forcing, which can have an implicit dependence on features such as the fjord stratification rather than the actual thermal forcing near the grounding line (i.e., T_{C}−T_{f}). We expect that a sudden increase in the inflowing AW temperature at the sill would initially cause warmer water to reach the grounding line and increase Q_{P} and M. However, adjustments of the temperature and stratification in the ice cavity are assumed to reestablish a state with the T_{A}independent melt rate determined by Eqs. (31) and (15). Since $\mathrm{\Delta}T={T}_{\mathrm{A}}T$ is constant, the outflow temperature T mirrors T_{A}.
4.2 Scenario 2: unconstrained plume volume flow
4.2.1 Physical assumptions
We will now consider another hypothetical scenario in which the plume volume transport is determined by the thermal forcing via Eq. (16) and is not set directly by the hydraulic constraints. This implies that Q_{P}>Q_{H}, and additional physical processes need to be invoked to balance the production and export over the sill of glacially modified water. Entrainment of some glacially modified water into the inflowing AW will help to achieve this balance: the entrained glacially modified water is, in effect, recirculating in the basin inside the sill. This can allow a steady state to develop even if Q_{P} exceeds Q_{H}. Specifically, in this scenario, we assume the following.

The excess volume flow in the plume is supplied by entrainment of glacially modified water (Q_{E}) into the inflow of AW on the landward side of the sill:
$$\begin{array}{}\text{(34)}& {Q}_{\mathrm{P}}={Q}_{\mathrm{E}}+{Q}_{\mathrm{H}}.\end{array}$$By introducing the entrainment fraction (Price and O'Neil Baringer, 1994)
$$\begin{array}{}\text{(35)}& \mathrm{\Phi}\stackrel{\mathrm{def}}{=}{\displaystyle \frac{{Q}_{\mathrm{E}}}{{Q}_{\mathrm{H}}+{Q}_{\mathrm{E}}}}=\mathrm{1}{\displaystyle \frac{{Q}_{\mathrm{H}}}{{Q}_{\mathrm{P}}}},\end{array}$$we can relate the plume volume flow and the exchange flow as
$$\begin{array}{}\text{(36)}& {Q}_{\mathrm{P}}={\displaystyle \frac{{Q}_{\mathrm{H}}}{\mathrm{1}\mathrm{\Phi}}}.\end{array}$$The parameter Φ ranges from 0 (no entrainment) to 1 (in the limit of strong entrainment). We note that Φ is essentially equal to the reflux factor used by Hager et al. (2022).

The relationship between basal melt M and thermal forcing (Eq. 15) still applies and equals the formula for M in the hydraulically controlled regime (Eq. 25).
The second assumption yields the following relationship:
By using this result in Eq. (23), we obtain
These formulas, which depend on features of the basal melt as well as h and k_{H}, specify the flow dependence on the thermal forcing in this hydraulicregime scenario. However, 𝒯 is a function of the AW forcing 𝒯_{A} and h that remains to be determined. This is obtained by considering how the groundingline temperature is affected by the entrainment of glacially modified waters as outlined below.
The entrainment is controlled by local conditions on the landward side of the sill, where the denser inflowing AW accelerates down the sill slope (Price and O'Neil Baringer, 1994; Pratt and Whitehead, 2007). However, we assume for simplicity that Q_{E} adjusts to satisfy Eq. (34). By using Eqs. (16) and (38), after some manipulation we can express the entrainment rate Φ as
where
Here Φ is given by Eq. (39) when Z≤1, and when Z>1 Φ=0. Note that Z=1 yields the condition (Eq. 27) that defines the flowregime transition.
Next, we consider the relationship between the AW (S_{A}, T_{A}) and the water properties in the ice cavity (S_{C}, T_{C}), which are affected by the entrainment: conservation of heat in the cavity (see Fig. 1) yields
By using the entrainment parameter Φ (Eq. 35), this yields the following expression for the temperature in the ice cavity:
Conservation of salt yields an analogous expression for S_{C}. Finally by using Eqs. (13) and (42), we obtain
This relation and Eqs. (37) and (39) determine the functional relationship $\mathcal{T}=\mathcal{T}({\mathcal{T}}_{\mathrm{A}},h)$ in this scenario. Note that since ΔT and Φ are specified as functions of 𝒯 and h, Eq. (43) also yields the function ${\mathcal{T}}_{\mathrm{A}}(\mathcal{T},h)=\mathcal{T}+\mathrm{\Delta}T\mathrm{\Phi}$, which is algebraically easier to use when constructing graphical solutions. This is because the function 𝒯(𝒯_{A},h) cannot, generally, be obtained on a closed analytical form; see Appendix B.
4.2.2 Dynamical features: general aspects
Figure 7 shows, for scenario 2, the dependence of the thermal forcing and the reduction factor on the AW forcing (𝒯_{A} and h). Again, the case of n_{1}=2 and n_{2}=1 is illustrated, which is qualitatively representative of melt representations satisfying 3n_{2}>n_{1}. Qualitatively, the behavior is similar to that of scenario 1 (Fig. 6). A difference is that the thermal forcing now increases with 𝒯_{A} in the hydraulic regime. However, the rate of increase is weaker than linear ($\frac{\mathrm{d}\mathcal{T}}{\mathrm{d}{\mathcal{T}}_{\mathrm{A}}}<\mathrm{1}$).
For the hydraulically controlled flow, entrainment of colder glacially modified waters into the inflowing water lowers the icecavity temperature relative to T_{A}. This decreases the thermal forcing and thereby the melt rate: 𝒯 decreases with decreasing h and increases more slowly with 𝒯_{A} than in the meltcontrolled regime. Figure 7b shows the reduction factor R (Eq. 26). By definition R=1 in the meltcontrolled regime. In the hydraulic regime, R<1, and the sensitivity of the thermal forcing to changes in 𝒯_{A} is reduced: the isolines of constant R become shallower with increasing 𝒯_{A}. In the present scenario 2, the flow response also depends on the parameter σ, a nondimensional measure of the outflow temperature T at the regime transition, which will be discussed below.
4.2.3 Dynamical features: dependence on AW height
Here, we examine the flow and melt response to changes in the AW height h for a fixed 𝒯_{A}, i.e., moving vertically in Fig. 7. Specifically, we consider how the parameter σ affects the response. Recall that σ is a nondimensional measure of the outflow temperature T in meltcontrolled regime where $\mathrm{\Delta}T/{\mathcal{T}}_{\mathrm{A}}$= σ; see Eq. (20). We will consider the whole range of possible σ values ($\mathrm{0}\le \mathit{\sigma}\le \mathrm{1}$), but our observationally based estimates indicate that σ is about 0.1 (Table 2).
Figure 8 illustrates how Q, entrainment fraction Φ, ΔT, 𝒯, and M vary with the AW height h. (In the figures, we have normalized Q, 𝒯, and M to be unity at the regime transition; but ΔT is normalized equal to σ.) If h is decreased, either by an increase in the sill height or by lowering the upper boundary of the AW, the flow is unchanged until h=h_{L}, the point at which the flow becomes hydraulically controlled. By further reducing h, Q, 𝒯, and M decline but ΔT grows (implying decreasing T). In the limiting case of σ=1, where ΔT is constant, Q is proportional to ${h}^{\mathrm{3}/\mathrm{2}}$; see Eq. (23). When σ<1, Q falls less steeply with h because the layer density difference (proportional to ΔT) increases with decreasing h.
In the hydraulic regime, outflowing glacially modified water is entrained into inflowing AW, thereby reducing 𝒯 (Fig. 8c). This effect is most pronounced for larger values of σ, which correspond to colder outflow temperatures of the glacially modified water. In the limiting case of σ=1, $M\propto {h}^{\frac{\mathrm{3}}{\mathrm{2}}}$ because ΔT=𝒯_{A} is constant; see Eq. (25). The dependence of M on h in this limiting case describes the response for all values of σ when $h/{h}_{L}$ becomes small. Here, $\mathrm{\Delta}T/{\mathcal{T}}_{\mathrm{A}}$ approaches its maximum value of 1. By using that ΔT≈𝒯_{A} in Eq. (25), the melt rate becomes
Notably, the melt rate is independent of the features of the melt representation in this limit. From the melt rate formula (Eq. 15), it follows that the thermal forcing is approximately given by
Thus, the thermal forcing still depends on the melt parameterization, but in such a way that the melt itself only depends on k_{H}, h, and 𝒯_{A}.
To summarize, in scenario 2 hydraulic control constrains the exchange flow and induces entrainment, which acts to lower the water temperature at the grounding line relative to T_{A}. Notably, when $h/{h}_{L}$ becomes sufficiently small, the flow enters a regime where the exchange flow and the basal melt become independent of the physical processes near the ice–ocean interface that govern the local melt rates. A partly analogous situation is an overmixed estuary, where hydraulic control at a sill or a fjord mouth sets the exchange flow rate and vertical salinity difference independently of the nature of the mixing processes in the estuary (Stommel and Farmer, 1953; Timmermans, 1998).
4.2.4 Dynamical features: dependence on AW temperature
Next, we consider how the melt and flow features depend on the AW temperature when h is fixed. Figure 9a and b show how the exchange flow Q (normalized to be unity at the regime transition), temperature difference ΔT, and the entrainment fraction Φ vary with 𝒯_{A}. In the meltcontrolled regime, Q and ΔT depend linearly on 𝒯_{A} (when n_{1}=2 and n_{2}=1), and a larger σ is associated with a larger ΔT. If 𝒯_{A} is increased, the flow enters the hydraulically controlled regime, where Q is proportional to $\mathrm{\Delta}{T}^{\mathrm{1}/\mathrm{2}}$. This constrains the exchange flow, and increasing entrainment lowers the temperature at the grounding line. In response, the outflow temperature decreases, which causes ΔT to increase with 𝒯_{A} at a rate that is slightly higher than linear. In the limiting case of σ=1, ΔT=𝒯_{A} in both regimes, and hence Q is proportional to ${\mathcal{T}}_{\mathrm{A}}^{\mathrm{1}/\mathrm{2}}$ in the hydraulic regime. The entrainment rate depends only weakly on σ and increases relatively slowly with 𝒯_{A}.
Figure 9c and d show the dependence of the thermal forcing and basal melt on 𝒯_{A}, which are both normalized to be unity at the regime transition. In the hydraulic regime, 𝒯 and M are lower for a given 𝒯_{A} than they would have been in a meltcontrolled regime. Since M∝ΔTQ, the hydraulicregime melt rate is proportional to ${\mathcal{T}}_{\mathrm{A}}^{\mathrm{3}/\mathrm{2}}$ when σ=1. Note that in Fig. 9, the melt rates have been normalized to be unity at the regime transition for visual clarity. The actual melt rates are proportional to ${\mathit{\sigma}}^{\mathrm{3}/\mathrm{2}}$, implying that σ=1 corresponds to the highest melt rate of a hydraulically controlled exchange flow for a given 𝒯_{A}.
4.2.5 Dynamical features: dependence on the melt parameterization exponents n_{1} and n_{2}
So far we have considered the case of $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}>\mathrm{0}$, where increasing 𝒯_{A} brings the flow towards the hydraulic regime (Fig. 4). The large 79^{∘}, Petermann, and Ryder ice tongues should be described by this case. However, some qualitatively different flow features emerge if $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}<\mathrm{0}$, and this case may be relevant for tidewater glaciers with high glacial discharge: theoretical considerations suggest that the exponents n_{1}=1 and n_{2}=0 describe the melt processes in this limit (Jenkins, 2011; Straneo and Cenedese, 2015). Therefore, we consider briefly two cases for which $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}\le \mathrm{0}$ in the context of scenario 2.
Figure 10 shows flow features for the cases $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}=\mathrm{0}$ and $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}=\mathrm{1}$. In the former case, the flow has the same dependence on 𝒯 and ΔT in both regimes; see Eqs. (16) and (38). As a result, the boundary between the flow regimes depends only on h. Further the R factor, the suppression of the thermal forcing due to hydraulic control, is independent of 𝒯_{A}.
Also in the case where $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}<\mathrm{0}$, a hydraulically controlled exchange flow suppresses the thermal forcing and basal melt. However, here the R factor – for a fixed h – increases with increasing 𝒯_{A}. Thus, as the AW temperature increases, the hydraulic suppression of the melt decreases. The reason is that the hydraulically determined upper bound on the exchange flow Q_{H} (Eq. 38) now increases faster with the thermal forcing than the plume volume transport Q_{P} (Eq. 16). Figure 10b and d illustrate the case of n_{1}=1 and n_{2}=0, where Q_{P} is independent of the thermal forcing but captures the qualitative features for the general case of $\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}}<\mathrm{0}$.
4.3 Applications to the Ryder, 79^{∘} N, and Petermann glaciers
To examine some concrete aspects of the model, we will now apply it in a qualitative way to the ice tongues of the Ryder, 79^{∘} N, and Petermann glaciers. We consider scenario 2, assuming n_{1}=2 and n_{2}=1 and try to crudely estimate model parameters characterizing the melt–flow dynamics. We use observations of flow rates and melt rates (Q and M) and hydrography from the three glaciers reported in the literature (Johnson et al., 2011; Wilson et al., 2017; Jakobsson et al., 2020; Schaffer et al., 2020). We recall that observations show that the sill exchange flows of the Ryder and 79^{∘} N glaciers are hydraulically controlled but that the exchange flow over the relatively deep and wide sill in Petermann Fjord is not. Note that there are uncertainties in the observations and in model assumptions, and the present exercise is primarily an illustration of how the model can be applied.
Guided by the model physics, we estimate model variables and parameters as follows: the hydrographic observations give T_{C} as the nearbottom temperatures inside the sill and T_{A} as the silldepth temperature outside the sill. Note that outside the sills, the temperatures are nearly constant below the sill depths; see Fig. 3. The outflow model temperature T, which represents a flowweighted mean over of an outflow distributed vertically over a range of temperatures, is less straightforward to determine from hydrography. Here, we use the relation $M/Q=\mathrm{\Delta}T/{T}_{\mathrm{G}}=\mathrm{\Delta}S/{S}_{\mathrm{A}}$ (which follows from Eqs. (3) and (9)) to find values of T and S that roughly satisfy these conditions and at the same time characterize outflowing water. This allows Φ to be determined from Eq. (42).
The model parameters γ_{1}, γ_{2}, and k_{H} are estimated as follows: from Eq. (15), we obtain ${\mathit{\gamma}}_{\mathrm{1}}\approx M/{\mathcal{T}}^{\mathrm{2}}$. By using Eqs. (15), (16), and (36), we obtain ${\mathit{\gamma}}_{\mathrm{2}}\approx Q/\left[\mathcal{T}\right(\mathrm{1}\mathrm{\Phi}\left)\right]$. This provides the estimate
By using Eqs. (23) and (25), we obtain ${k}_{\mathrm{H}}{h}^{\mathrm{3}/\mathrm{2}}\approx {Q}^{\mathrm{3}/\mathrm{2}}/{M}^{\mathrm{1}/\mathrm{2}}$, and an estimate of h then gives k_{H}. Note that k_{H} can also be determined from a knowledge of the crosssectional sill width W.
The estimates of Q from Ryder and Petermann are more uncertain than the ones from 79^{∘} N reported by Schaffer et al. (2020), which are based on a 1year moored time series of velocity. The Ryder estimate of Q is based on a single instantaneous current measurement on the inner sill in Sherard Osborn Fjord (Jakobsson et al., 2020). Our Petermann estimate of Q is based on hydrography and geostrophic velocities presented by Johnson et al. (2011): using their Fig. 7, we estimate the outflow of glacially modified water (in the depth range of 150 to 250 m) to be on the order of 50×10^{3} m^{3} s^{−1}.
Table 2 summarizes observational features and estimated model parameters. Notably, the estimates of σ, around 0.1–0.2, are similar for the three glaciers. This reflects that the cooling due to ice melt of the meltwater plume is small compared to the upper limit (σ=1), in which the outflow temperature (T) approaches the freezing point. For the Ryder and 79^{∘} N glaciers, which have hydraulically controlled exchange flows, the estimated $h/{h}_{L}$ is 0.7 and 0.4, respectively. The lower value of $h/{h}_{L}$ at 79^{∘} N Glacier implies a higher sensitivity of the melt to changes in $h/{h}_{L}$; see Fig. 8. Further, the reduction factors ($R=\mathcal{T}/{\mathcal{T}}_{\mathrm{A}}$), which are directly inferred from the hydrography, are only slightly below unity: R is 0.96 and 0.86 at Ryder and 79^{∘} N, respectively. If as assumed here M∝𝒯^{2}, this implies that, relative to the situation where unmodified AW reaches the grounding line, the melt rates are reduced by about 10 % and 30 % at Ryder and 79^{∘} N, respectively. Taken together, this suggests that hydraulic control and associated entrainment reduce the basal melt on both ice tongues but that currently this effect is more pronounced at the 79^{∘} N Glacier.
Figure 11 shows modelpredicted melt rates M and R factors as a function of 𝒯_{A} for the Ryder and 79^{∘} N glaciers when h is kept at observed values. The figure also shows the melt rates that would occur if unmodified AW reached the grounding lines: the melt rates in the hydraulically controlled regime are lower and their dependence on 𝒯_{A} is weaker. The model results show that a 1 ^{∘}C increase in 𝒯_{A} from present values increases M by about 75 % and 50 % at Ryder and 79^{∘} N, respectively. The basal melt is more sensitive to the same change in T_{A} at Ryder than at 79^{∘} N simply because 𝒯_{A} is larger at 79^{∘} N. For fractional changes in 𝒯_{A}, the response in basal melt is more similar: a 10 % increase in 𝒯_{A} yields an increase in M of about 15 % at both glaciers. (Note that a 10 % increase in 𝒯_{A} corresponds to an increase in ∼0.3 and ∼0.4 ^{∘}C of T_{A} at Ryder and 79^{∘} N, respectively). In the absence of hydraulic control, where $M\propto {\mathcal{T}}_{\mathrm{A}}^{\mathrm{2}}$, the corresponding increase in M would be 20 %. Figure 12 shows the modelpredicted dependence of Ryder basal melt on T_{A} and h. Currently, the AW interface is about 30 m below the transition depth h_{L} at which hydraulic control ceases; i.e., a lifting of the AW interface by about 30 m would bring the flow into the meltcontrolled regime.
At Petermann Glacier, the exchange flow is not hydraulically controlled because the height of the AW above the sill is larger than the transition height h_{L} (Eq. 28). Note that h_{L} decreases with the crosssectional width of the sill (${h}_{L}\propto {W}^{\mathrm{2}/\mathrm{3}}$). Ryder and Petermann have similar values of 𝒯_{A}, and it is primarily the wide sill in Petermann Fjord that yields a small value of h_{L} (see Fig. 5 in Jakobsson et al., 2020).
To analyze the impact of hydraulic control on basal glacier melt, we have developed a twolayer fjord model that includes simple representations of melt and exchange flow dynamics. Despite model idealizations, we believe that Figs. 6 and 7 qualitatively illustrate how the interplay between nearice melt processes and hydraulic control affects the relationship between basal icetongue melt and AW features.^{3} Our results suggest that there are two flow regimes, with different relationships between basal melt and AW features. To begin with, there is a meltcontrolled flow regime, in which the fjord geometry and sills do not restrict the exchange flow, and unmodified AW reaches the grounding line (Fig. 1). In this regime, the basal melt is rate limited by nearice processes rather than by the heat flux carried by the horizontal exchange flow in the fjord. Accordingly, the thermal forcing is set by the AW temperature and the basal melt processes determine the strength of the exchange flow. In this regime, the dependence of the basal melt on the AW height is weak and neglected here.
In a sill fjord, hydraulic control can set an upper bound on the exchange flow, which depends on the height of the AW above the sill crest and the density difference between the in and outflowing waters. If hydraulic control is established, the outflow of glacially modified water created by the basal melt must be compatible with the hydraulically determined exchange flow at the sill to ensure volume conservation. In a hydraulically controlled flow regime, accordingly, the heat transport supplied by the fjord circulation enters as a ratelimiting factor for the basal melt (Fig. 1). The flow can transit from a meltcontrolled to a hydraulically controlled regime if the AW height is decreased or the AW temperature is increased (Figs. 6, 7). Such transitions can also occur when increasing subglacial discharge enhances the basal melt and the production of glacially modified water.
In the hydraulic regime, decreasing AW height causes the exchange flow and its associated heat flux to decrease. As a result, the basal melt decreases. The hydraulic constraint also reduces the sensitivity of the basal melt to changes in the AW temperature. We have examined this effect using a simple representation of the melt processes (Eqs. 15 and 16), and considered two idealized scenarios for how the flow adjusts to satisfy the hydraulic constraint; see Sect. 4.1 and 4.2. In these scenarios changes in the fjord stratification or entrainment of glacially modified water into the inflowing AW are assumed to regulate the thermal forcing such that the hydraulic constraint is satisfied. In scenario 1 (Fig. 6), the thermal forcing is completely blind to the AW temperature but sensitive to the AW height. Scenario 2 (Fig. 7) is less extreme, and here the thermal forcing has a muted response to changes in the AW temperature; i.e., $\frac{\mathrm{d}\mathcal{T}}{\mathrm{d}{\mathcal{T}}_{\mathrm{A}}}<\mathrm{1}$. These scenarios involve some fairly ad hoc assumptions, and further studies are needed to more accurately quantify the suppression of the thermal forcing in hydraulic flow regimes. Nevertheless, the qualitative features shown in Figs. 6 and 7 are expected to be robust. We note that scenario 2, which assumes entrainment and recirculation in the ice cavity, is more consistent with observations and modeling (Schaffer et al., 2020; Jakobsson et al., 2020; Bao and Moffat, 2023) than scenario 1.
The suppression of basal melt due to hydraulic control can be quantified by the reduction factor R (Eq. 26): the melt relative to the case when AW reaches the grounding line is proportional to ${R}^{{n}_{\mathrm{1}}}$. At an ice tongue or tidewater glacier in a specific fjord, R is a function of the AW features; i.e., $R=R({\mathcal{T}}_{\mathrm{A}},h)$. This feature could be used to parameterize effects of hydraulic control in simulations of marineglacier response to changes in AW forcing, which is crucial for the evolution of Greenlandic marine glaciers on decadal and centennial timescales (Straneo and Heimbach, 2013; Aschwanden et al., 2019; Wood et al., 2021).
We have considered a situation where the sill is seaward of the icetongue front, which is presently the case for the Ryder and 79^{∘} N glaciers. However, if an ice tongue extends above the sill, the ice draft will contribute to the geometrical constraints that determine the hydraulic exchange flow: the ice reduces the water column depth over the sill. We will not explore this problem here. However, we note that in the early 1900s the Ryder ice tongue was some 40 km longer than today, and covered the inner sill; its front reached roughly the 82^{∘}12^{′} N mark in Fig. 2c; see Jakobsson et al. (2020) and O'Regan et al. (2021) for additional information. This should have strongly restricted the water exchange over the inner sill, resulting in very low basal melt on the inner part of the ice tongue.
Ryder Glacier has been relatively stable in recent decades (Hill et al., 2018). In contrast, Petermann Glacier, located ∼ 200 km southwest of Ryder, has been retreating and lost 35 km of its ice tongue in 2010 and 2012 (Johannessen et al., 2013; Hill et al., 2018). Jakobsson et al. (2020) proposed that Ryder Glacier has been stable because of its more restrictive sill geometry, which partly protects the ice tongue from the inflow of warmer subsurface AW (Fig. 2). The present study suggests that Ryder has a relatively high R value (R≈0.9), implying that despite the double sill geometry in Sherard Osborne Fjord, the modified AW reaching the grounding line is currently weakly cooled as its flows through the fjord. However, the sensitivity of basal melt to thermal forcing depends on local conditions such as icetongue geometry, subglacial discharge, and tidal currents. Notably, our simple model fit (Table 2) suggests that the thermal sensitivity of basal melt per unit area (${\mathit{\gamma}}_{\mathrm{1}}/A$) is 40 % higher for Petermann than for the Ryder and 79^{∘} N glaciers, which have comparable sensitivities. Further, remotesensing analyses show that the basal melt per unit area ($M/A$) is about 50 % higher on Petermann than on Ryder (Wilson et al., 2017). Even if our model fit is quite uncertain, this indicates that Ryder and 79^{∘} N glaciers, which are shielded by hydraulically controlled sill flows, also have basal melt processes characterized by lower thermal sensitivity coefficients (${\mathit{\gamma}}_{\mathrm{1}}/A$) than Petermann.
We emphasize that the oceanic conditions at Ryder have only been observed a single time in the summer of 2019 and may not give a representative view of the melt–flow dynamic. As documented in the observations of Schaffer et al. (2020) at 79^{∘} N, variations in the AW height on monthly to annual timescales cause significant variations in exchange flow and basal melt. Thus, in hydraulically controlled fjords, longterm melt variations may be strongly controlled by the evolution of the AW height, a quantity that has received less attention than the AW temperature for the evolution of marine glaciers in Greenland (Straneo and Heimbach, 2013; Wood et al., 2021). Central Arctic Ocean observations document changes, on decadal timescales, of the AW height that are up to 100 m (Polyakov et al., 2004). If similar height changes occurred along the Arctic coast of Greenland, significant changes in basal glacial melt would result: the present model (scenario 2) suggests that a lowering of the AW height by ∼ 40 m would halve the basal melt on the Ryder ice tongue.
In the conservation relations of Sect. 3.1, subglacial discharge (say D) is neglected, whereas it is allowed to affect the melt rate; see Eq. (15). Neglecting D in the conservation relations is generally not a valid approximation for tidewater glaciers, and we show here how to generalize the results to cases where the subglacial discharge is not small compared to the freshwater input due to subsurface ice melt (M). Essentially, this is accomplished by replacing M by M+D in the derivations presented in Sect. 3.1, and this is outlined below.
When including D, conservation of volume is given by
but salt conservation is still given by Eq. (2). This yields the modified Knudsen's relation:
The advective heat flux (Eq. 5) becomes
The advective heat flux determines the melt rate (Eq. 6), which in combination with Eq. (A3) yields
where the Gade temperature T_{G} is defined in Eq. (8). This equation is the modified form of Eq. (7). For conditions in northern Greenland $({T}_{\mathrm{A}}{T}_{\mathrm{f}})/{T}_{\mathrm{G}}\approx \mathrm{0.05}$, implying that the last term in Eq. (A4) can be neglected unless D≫M. This approximation is made in Sect. 2.1, where D is taken to be 0.
By combining Eqs. (A2) and (A4) and eliminating Q, we obtain the counterpart of Eq. (10):
where we have introduced
Since $({T}_{\mathrm{A}}{T}_{\mathrm{f}})/{T}_{\mathrm{G}}<\mathrm{1}$, it follows that Γ≥1. The density difference (Eq. 12) as function of ΔT becomes
Thus for a given ΔT, the primary effect of subglacial discharge is to increase the associated ΔS and Δρ.
Figure A1 shows that depending on the value of $D/M$, there are two limiting regimes.

When $D/M\ll \mathrm{1}$, Γ≈1. Formally, this is the limit considered in section 2.1, where $D/M$ is taken to be 0. However, Fig. A1 indicates that this limit may also serve as a leading order approximation when $D/M\approx \mathrm{1}$.

When $D/M\gg \mathrm{1}$, $\mathrm{\Gamma}\approx {T}_{\mathrm{G}}/({T}_{\mathrm{A}}{T}_{\mathrm{f}})$ (≈20 for conditions in northern Greenland). This implies that $\mathrm{\Delta}S/{S}_{\mathrm{A}}\approx \mathrm{\Delta}T/({T}_{\mathrm{A}}{T}_{\mathrm{f}})$, which is the relationship between salinity and temperature changes when freshwater at the freezing temperature is mixed with AW. Here, βS_{A}Γ≫αT_{G}, and from Eqs. (A5) and (A7) it follows that the density difference becomes approximately controlled by the salinity difference alone: $\mathrm{\Delta}\mathit{\rho}/{\mathit{\rho}}_{\mathrm{0}}\approx \mathit{\beta}\mathrm{\Delta}S$. This limit is approached when $D/M$ is large compared to ${T}_{\mathrm{G}}/({T}_{\mathrm{A}}{T}_{\mathrm{f}})$ and can be appropriate for tidewater glaciers with high subglacial discharge.
As long as the volume flux of subglacial discharge and basal melt is small compared to the fjord exchange flow (the generally valid case for Greenlandic fjords where $Q\gg M+D$), the inclusion of the factor Γ in the generalized Eqs. (A5) and (A7) is the only model modification needed for treating cases where the subglacial discharge is not small compared to the melt. The flow in the meltcontrolled regime (Sect. 3.3) does not depend on ΔS and Δρ and is therefore not dependent on the value of Γ. To describe the hydraulically controlled flow regime (Sect. 3.4), it is convenient to define a modified hydraulic coefficient:
where k_{H} is defined in Eq. (24), and we note that ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}\ge {k}_{\mathrm{H}}$. Replacing k_{H} with ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}$ in Sect. 3.4 allows us to describe cases with high subglacial discharge.
For given AW features associated with a specific ΔT, the primary effect of subglacial discharge (besides increasing the subsurface melt) is to enhance the layer density difference. Essentially, this causes the transition into the hydraulically controlled regime to occur for somewhat greater sills heights (or a lower height of the AW layer) than when $D/M$ is taken to be 0. The regime transition is still defined from the condition that Q_{P}=Q_{H} (see Eq. (27)), which yields the equivalent of Eq. (28):
where k_{H} is replaced ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}$. This shows that h_{L} (the height of the AW layer above the sill for which the flow becomes hydraulically controlled) decreases when Γ increases ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}$. In the limits when $D/M\ll \mathrm{1}$ or $D/M\gg \mathrm{1}$, Γ and therefore ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}$ are constants independent of M and D. Accordingly, the results for the limit of small subglacial discharge ($D/M\ll \mathrm{1}$) in the present paper are qualitatively similar to those for the limit of high subglacial discharge ($D/M\gg \mathrm{1}$).
The situation is slightly more complicated between these two limits, where D∼M. This is because Γ and ${\stackrel{\mathrm{\u0303}}{k}}_{\mathrm{H}}$ in this regime both depend on M and D, which yields a model that is more complex algebraically. However, the main effect of finite values of $D/M$ is still to decrease the transition height h_{L}. Thus, the model results should qualitatively also describe cases where D∼M.
Here, we derive a few mathematical relationships that can be used to construct graphs for scenario 2. We assume that ${n}_{\mathrm{1}}{n}_{\mathrm{2}}=\mathrm{1}$, which simplifies the algebra but is not strictly necessary.
When 𝒯_{A} is fixed, it is convenient to put Eq. (43) in nondimensional form using the variables
where h_{L} is defined in Eq. (28). Note that $\stackrel{\mathrm{\u0303}}{\mathcal{T}}=R$; see Eq. (26). By using these nondimensional variables and Eqs. (37) and (35), we obtain
Here, the term in the square brackets on the righthand side is Φ, and σ is defined by Eq. (20). From Eq. (B2), we can obtain $\stackrel{\mathrm{\u0303}}{h}$ as a function of $\stackrel{\mathrm{\u0303}}{\mathcal{T}}$:
It is generally not possible to find the inverse function $\stackrel{\mathrm{\u0303}}{\mathcal{T}}=\stackrel{\mathrm{\u0303}}{\mathcal{T}}\left(\stackrel{\mathrm{\u0303}}{h}\right)$ in a closed analytical form, but Eq. (B3) allows us to examine it graphically by plotting 𝒯 versus h. By using Eqs. (37) and (B3), we can after some manipulations find the dependence of ΔT and Φ on $\stackrel{\mathrm{\u0303}}{\mathcal{T}}$:
To examine the melt dynamics when h is fixed, it is useful to rewrite Eq. (43) as
where the terms on the righthand side are now known functions of 𝒯 and h specified by Eqs. (37) and (35). We put Eq. (B6) in nondimensional form using the definitions of 𝒯_{L}, ΔT_{L}, and h_{L}:
By using Eq. (B7), after some straightforward calculations we obtain
which applies at the regime transition in the hydraulically controlled regime. In the meltcontrolled regime, $\frac{\mathrm{d}\mathcal{T}}{\mathrm{d}{T}_{\mathrm{A}}}=\mathrm{1}$. Thus, the suppression of 𝒯 relative to the AW thermal forcing (𝒯_{A}) is governed by the exponents n_{1} and n_{2} and σ. An inspection of Eq. (35) shows that $(\mathrm{3}{n}_{\mathrm{2}}{n}_{\mathrm{1}})/\mathrm{3}$ determines how fast the entrainment increases with 𝒯.
Data presented in the paper (multibeam bathymetry and oceanographic stations) are available in the Bolin Centre for Climate Research database. Multibeam bathymetry: https://doi.org/10.17043/odenryder2019bathymetry1 (Calder et al., 2020). LADCP (current measurements): https://doi.org/10.17043/odenryder2019ladcp1 (Stranne et al., 2020a). CTD stations: https://doi.org/10.17043/odenryder2019ctd1 (Stranne et al., 2020b).
JN lead the work on the paper. All authors contributed to the writing.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We thank Jonas Nycander, Jonathan Wiskandt, and Inga Koszalka for valuable comments on our work. We also thank two anonymous reviewers, who offered many insightful suggestions for revising the paper.
This research has been supported by the Vetenskapsrådet (grant nos. 201804350, 202005076, 202104512, and 202204018) and the Swedish National Space Agency (grant no. 202000171).
The article processing charges for this openaccess publication were covered by Stockholm University.
This paper was edited by Jan De Rydt and reviewed by two anonymous referees.
Armi, L.: The hydraulics of two flowing layers with different densities, J. Fluid Mech., 163, 27–58, https://doi.org/10.1017/S0022112086002197, 1986. a
Aschwanden, A., Fahnestock, M. A., Truffer, M., Brinkerhoff, D. J., Hock, R., Khroulev, C., Mottram, R., and Khan, S. A.: Contribution of the Greenland Ice Sheet to sea level over the next millennium, Science Advances, 5, eaav9396, https://doi.org/10.1126/sciadv.aav9396, 2019. a, b
Bao, W. and Moffat, C.: Impact of shallow sills on heat transport and stratification regimes in proglacial fjords, The Cryosphere Discuss. [preprint], https://doi.org/10.5194/tc202332, in review, 2023. a, b, c
Cai, C., Rignot, E., Menemenlis, D., and Nakayama, Y.: Observations and modeling of oceaninduced melt beneath Petermann Glacier Ice Shelf in northwestern Greenland, Geophys. Res. Lett., 44, 8396–8403, https://doi.org/10.1002/2017GL073711, 2017. a, b, c
Calder, B., Eriksson, B., Jerram, K., Weidner, E., Holmes, F., Muchowski, J., Prakash, A., Handl, T., Ståhl, E., Mayer, L., and Jakobsson, M.: Highresolution bathymetry from the Ryder 2019 expedition to Northwest Greenland. Dataset version 1, Bolin Centre Database [data set], https://doi.org/10.17043/odenryder2019bathymetry1, 2020. a
De Andrés, E., Slater, D. A., Straneo, F., Otero, J., Das, S., and Navarro, F.: Surface emergence of glacial plumes determined by fjord stratification, The Cryosphere, 14, 1951–1969, https://doi.org/10.5194/tc1419512020, 2020. a
Favier, L., Jourdain, N. C., Jenkins, A., Merino, N., Durand, G., Gagliardini, O., GilletChaulet, F., and Mathiot, P.: Assessment of subshelf melting parameterisations using the ocean–icesheet coupled model NEMO(v3.6)–Elmer/Ice(v8.3) , Geosci. Model Dev., 12, 2255–2283, https://doi.org/10.5194/gmd1222552019, 2019. a, b, c
Gade, H.: Melting of ice in sea water: a primitive model with application to the Antarctic Shelf Ice and icebergs, J. Phys. Oceanogr., 9, 189–198, https://doi.org/10.1175/15200485(1979)009<0189:MOIISW>2.0.CO;2, 1979. a, b
Gudmundsson, G. H.: Iceshelf buttressing and the stability of marine ice sheets, The Cryosphere, 7, 647–655, https://doi.org/10.5194/tc76472013, 2013. a
Hager, A. O., Sutherland, D. A., Amundson, J. M., Jackson, R. H., Kienholz, C., Motyka, R. J., and Nash, J. D.: Subglacial Discharge Reflux and Buoyancy Forcing Drive Seasonality in a Silled Glacial Fjord, J. Geophys. Res.Oceans, 127, e2021JC018355, https://doi.org/10.1029/2021JC018355, 2022. a, b
Hattermann, T., Smedsrud, L. H., Nøst, O. A., Lilly, J. M., and GaltonFenzi, B.: Eddyresolving simulations of the Fimbul Ice Shelf cavity circulation: Basal melting and exchange with open ocean, Ocean Model., 82, 28–44, 2014. a
Hill, E. A., Carr, J. R., Stokes, C. R., and Gudmundsson, G. H.: Dynamic changes in outlet glaciers in northern Greenland from 1948 to 2015, The Cryosphere, 12, 3243–3263, https://doi.org/10.5194/tc1232432018, 2018. a, b, c, d
Holland, D. M. and Jenkins, A.: Modeling Thermodynamic Ice–Ocean Interactions at the Base of an Ice Shelf, J. Phys. Oceanogr., 29, 1787–1800, https://doi.org/10.1175/15200485(1999)029<1787:MTIOIA>2.0.CO;2, 1999. a, b
Holland, P. R., Jenkins, A., and Holland, D. M.: The Response of Ice Shelf Basal Melting to Variations in Ocean Temperature, J. Climate, 21, 2558–2572, https://doi.org/10.1175/2007JCLI1909.1, 2008. a, b, c, d, e
Jackson, R. H. and Straneo, F.: Heat, Salt, and Freshwater Budgets for a Glacial Fjord in Greenland, J. Phys. Oceanogr., 46, 2735–2768, https://doi.org/10.1175/JPOD150134.1, 2016. a, b, c
Jackson, R. H., Lentz, S. J., and Straneo, F.: The Dynamics of Shelf Forcing in Greenlandic Fjords, J. Phys. Oceanogr., 48, 2799–2827, https://doi.org/10.1175/JPOD180057.1, 2018. a
Jakobsson, M., Mayer, L. A., Nilsson, J., Stranne, C., Calder, B., O'Regan, M., Farrell, J. W., Cronin, T. M., Brüchert, V., Chawarski, J., Eriksson, B., Fredriksson, J., Gemery, L., Glueder, A., Holmes, F. A., Jerram, K., Kirchner, N., Mix, A., Muchowski, J., Prakash, A., Reilly, B., Thornton, B., Ulfsbo, A., Weidner, E., Åkesson, H., Handl, T., Ståhl, E., Boze, L.G., Reed, S., West, G., and Padman, J.: Ryder Glacier in northwest Greenland is shielded from warm Atlantic water by a bathymetric sill, Communications Earth & Environment, 1, 45, https://doi.org/10.1038/s43247020000430, 2020. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r
Jenkins, A.: A onedimensional model of ice shelfocean interaction, J. Geophys. Res.Oceans, 96, 20671–20677, https://doi.org/10.1029/91JC01842, 1991. a, b
Jenkins, A.: The Impact of Melting Ice on Ocean Waters, J. Phys. Oceanogr., 29, 2370–2381, https://doi.org/10.1175/15200485(1999)029<2370:TIOMIO>2.0.CO;2, 1999. a
Jenkins, A.: ConvectionDriven Melting near the Grounding Lines of Ice Shelves and Tidewater Glaciers, J. Phys. Oceanogr., 41, 2279–2294, https://doi.org/10.1175/JPOD1103.1, 2011. a, b, c, d, e, f
Johannessen, O. M., Babiker, M., and Miles, M. W.: Unprecedented retreat in a 50year observational record for Petermann Glacier, North Greenland, Atmospheric and Oceanic Science Letters, 6, 259–265, https://doi.org/10.3878/j.issn.16742834.13.0021, 2013. a, b
Johnson, H. L., Münchow, A., Falkner, K. K., and Melling, H.: Ocean circulation and properties in Petermann Fjord, Greenland, J. Geophys. Res.Oceans, 116, C01003, https://doi.org/10.1029/2010JC006519, 2011. a, b, c
Khazendar, A., Fenty, I. G., Carroll, D., Gardner, A., Lee, C. M., Fukumori, I., Wang, I., Zhang, H., Seroussi, H., Moller, D., Noël, B. P. Y., van den Broeke, M. R.,Dinardo, S., and Willis, J.: Interruption of two decades of Jakobshavn Isbrae acceleration and thinning as regional ocean cools, Nat. Geosci., 12, 277–283, 2019. a
Lazeroms, W. M. J., Jenkins, A., Gudmundsson, G. H., and van de Wal, R. S. W.: Modelling presentday basal melt rates for Antarctic ice shelves using a parametrization of buoyant meltwater plumes, The Cryosphere, 12, 49–70, https://doi.org/10.5194/tc12492018, 2018. a, b
Lazeroms, W. M. J., Jenkins, A., Rienstra, S. W., and van de Wal, R. S. W.: An Analytical Derivation of IceShelf Basal Melt Based on the Dynamics of Meltwater Plumes, J. Phys. Oceanogr., 49, 917–939, https://doi.org/10.1175/JPOD180131.1, 2019. a
Lewis, E. L. and Perkin, R. G.: Ice pumps and their rates, J. Geophys. Res., 91, 11756–11762, https://doi.org/10.1029/JC091iC10p11756, 1986. a
Lindeman, M. R., Straneo, F., Wilson, N. J., Toole, J. M., Krishfield, R. A., Beaird, N. L., Kanzow, T., and Schaffer, J.: Ocean Circulation and Variability Beneath Nioghalvfjerdsbræ (79 North Glacier) Ice Tongue, J. Geophys. Res.Oceans, 125, e2020JC016091, https://doi.org/10.1029/2020JC016091, 2020. a
Little, C. M., Gnanadesikan, A., and Hallberg, R.: LargeScale Oceanographic Constraints on the Distribution of Melting and Freezing under Ice Shelves, J. Phys. Oceanogr., 38, 2242–2255, https://doi.org/10.1175/2008JPO3928.1, 2008. a
Little, C. M., Gnanadesikan, A., and Oppenheimer, M.: How ice shelf morphology controls basal melting, J. Geophys. Res., 114, C12007, https://doi.org/10.1029/2008JC005197, 2009. a, b, c
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, A., 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 Sounding Combined With Mass Conservation, Geophys. Res. Lett., 44, 11051–11061, https://doi.org/10.1002/2017GL074954, 2017. a
Mouginot, J., Rignot, E., Bjørk, A. A., van den Broeke, M., Millan, R., Morlighem, M., Noël, B., Scheuchl, B., and Wood, M.: Fortysix years of Greenland Ice Sheet mass balance from 1972 to 2018, P. Natl. Acad. Sci. USA, 116, 9239–9244, https://doi.org/10.1073/pnas.1904242116, 2019. a, b
Nick, F. M., Vieli, A., Andersen, M. L., Joughin, I., Payne, A., Edwards, T. L., Pattyn, F., and van de Wal, R. S. W.: Future sealevel rise from Greenland's main outlet glaciers in a warming climate, Nature, 497, 235–238, https://doi.org/10.1038/nature12068, 2013. a
Nycander, J., Hogg, A. M., and Frankcombe, L. M.: Open boundary conditions for nonlinear channel flow, Ocean Model., 24, 108–121, https://doi.org/10.1016/j.ocemod.2008.06.003, 2008. a
O'Regan, M., Cronin, T. M., Reilly, B., Alstrup, A. K. O., Gemery, L., Golub, A., Mayer, L. A., Morlighem, M., Moros, M., Munk, O. L., Nilsson, J., Pearce, C., Detlef, H., Stranne, C., Vermassen, F., West, G., and Jakobsson, M.: The Holocene dynamics of Ryder Glacier and ice tongue in north Greenland, The Cryosphere, 15, 4073–4097, https://doi.org/10.5194/tc1540732021, 2021. a
Polyakov, I. V., Alekseev, G. V., Timokhov, L. A., Bhatt, U. S., Colony, R. L., Simmons, H. L., Walsh, D., Walsh, J. E., and Zakharov, V. F.: Variability of the Intermediate Atlantic Water of the Arctic Ocean over the Last 100 Years, J. Climate, 17, 4485–4497, https://doi.org/10.1175/JCLI3224.1, 2004. a
Pratt, L. J. and Whitehead, J. A.: Rotating Hydraulics: Nonlinear Topographic Effects in the Ocean and Atmosphere, first edn., edited by: Mysak, L. A., Springer Verlag, ISBN 9780387495729, https://doi.org/10.1007/9780387495729, 2007. a, b, c, d, e, f, g, h, i, j
Price, J. F. and O'Neil Baringer, M.: Outflows and deep water production by marginal seas, Prog. Oceanogr., 33, 161–200, https://doi.org/10.1016/00796611(94)900272, 1994. a, b, c
Schaffer, J., T. Kanzow, W. J. v., L. von Albedyll, J. E. A., and Roberts, D. H.: Bathymetry constrains ocean heat supply to Greenland's largest glacier tongue, Nat. Geosci., 13, 227–231, https://doi.org/10.1038/s415610190529x, 2020. a, b, c, d, e, f, g, h, i, j, k, l
Schoof, C.: Ice sheet grounding line dynamics: Steady states, stability, and hysteresis, J. Geophys. Res., 112, F03S28, https://doi.org/10.1029/2006JF000664, f03S28, 2007. a
Schoof, C., Davis, A. D., and Popa, T. V.: Boundary layer models for calving marine outlet glaciers, The Cryosphere, 11, 2283–2303, https://doi.org/10.5194/tc1122832017, 2017. a
Shroyer, E. L., Padman, L., Samelson, R. M., Münchow, A., and Stearns, L. A.: Seasonal control of Petermann Gletscher iceshelf melt by the ocean's response to seaice cover in Nares Strait, J. Glaciol., 63, 324–330, https://doi.org/10.1017/jog.2016.140, 2017. a
Slater, D. A. and Straneo, F.: Submarine melting of glaciers in Greenland amplified by atmospheric warming, Nat. Geosci., 15, 794–799, https://doi.org/10.1038/s41561022010359, 2022. a, b
Stommel, H. and Farmer, H. G.: Control of salinity in an estuary by a transition, J. Mar. Res, 12, 13–20, 1953. a
Straneo, F. and Cenedese, C.: The Dynamics of Greenland's Glacial Fjords and Their Role in Climate, Annu. Rev. Mar. Sci., 7, 89–112, https://doi.org/10.1146/annurevmarine010213135133, 2015. a, b, c, d, e
Straneo, F. and Heimbach, P.: North Atlantic warming and the retreat of Greenland's outlet glaciers, Nature, 504, 36–43, https://doi.org/10.1038/nature12854, 2013. a, b, c, d
Straneo, F., Sutherland, D., Holland, D., Gladish, C., Hamilton, G., Johnson, H., and Koppes, M.: Characteristics of ocean waters reaching Greenland's glaciers, Ann. Glaciol., 53, 202–210, https://doi.org/10.3189/2012AoG60A059, 2012. a
Stranne, C., Nilsson, J., and Muchowski, J.: Oceanographic LADCP data from the Ryder 2019 expedition, Bolin Centre Database [data set], https://doi.org/10.17043/odenryder2019ladcp1, 2020a. a
Stranne, C., Nilsson, J., Muchowski, J., and Chawarski, J.: Oceanographic CTD data from the Ryder 2019 expedition, Dataset version 1, Bolin Centre Database [data set], https://doi.org/10.17043/odenryder2019ctd1, 2020b. a
Stranne, C., Nilsson, J., Ulfsbo, A., O'Regan, M., Coxall, H. K., Meire, L., Muchowski, J., Mayer, L. A., Brüchert, V., Fredriksson, J., Thornton, B., Chawarski, J., West, G., Weidner, E., and Jakobsson, M.: The climate sensitivity of northern Greenland fjords is amplified through seaice damming, Communications Earth & Environment, 2, 70, https://doi.org/10.1038/s43247021001408, 2021. a, b
Timmermans, M.L.: Hydraulic Control and Mixing in a SemiEnclosed Reservoir, Tech. rep., Geophys. Fluid Dyn. Summer Study Program Tech. Rep. WHOI9809, 1998. a
Truffer, M. and Motyka, R. J.: Where glaciers meet water: Subaqueous melt and its relevance to glaciers in various settings, Rev. Geophys., 54, 220–239, https://doi.org/10.1002/2015RG000494, 2016. a, b, c
Washam, P., Nicholls, K. W., Münchow, A., and Padman, L.: Summer surface melt thins Petermann Gletscher Ice Shelf by enhancing channelized basal melt, J. Glaciol., 65, 662–674, https://doi.org/10.1017/jog.2019.43, 2019. a
Wilson, N., Straneo, F., and Heimbach, P.: Satellitederived submarine melt rates and mass balance (2011–2015) for Greenland's largest remaining ice tongues, The Cryosphere, 11, 2773–2782, https://doi.org/10.5194/tc1127732017, 2017. a, b, c
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, Science Advances, 7, eaba7282, https://doi.org/10.1126/sciadv.aba7282, 2021. a, b, c, d, e
Xu, Y., Rignot, E., Fenty, I., Menemenlis, D., and Flexas, M. M.: Subaqueous melting of Store Glacier, west Greenland from threedimensional, highresolution numerical modeling and ocean observations, Geophys. Res. Lett., 40, 4648–4653, https://doi.org/10.1002/grl.50825, 2013. a, b, c, d
Zhao, K. X., Stewart, A. L., and McWilliams, J. C.: SillInfluenced Exchange Flows in Ice Shelf Cavities, J. Phys. Oceanogr., 49, 163–191, https://doi.org/10.1175/JPOD180076.1, 2019. a
Zhao, K. X., Stewart, A. L., and McWilliams, J. C.: Geometric Constraints on Glacial Fjord–Shelf Exchange, J. Phys. Oceanogr., 51, 1223–1246, https://doi.org/10.1175/JPOD200091.1, 2021. a, b
Note that the equivalent ice temperature defined by Jenkins (1999) is approximately ${T}_{\mathrm{G}}\cdot \frac{{\mathit{\rho}}_{w}}{{\mathit{\rho}}_{i}}$, where $\frac{{\mathit{\rho}}_{w}}{{\mathit{\rho}}_{i}}\approx \mathrm{1.1}$ is the density ratio between water and ice.
For simplicity, we assume a rectangular cross section, which implies that W does not depend on h. Note that the lowerlayer width W may be smaller than the fjord width if the inflow is confined in a deeper channel crossing the sill, which is the case for Ryder; see Fig. 1 in Jakobsson et al. (2020).
 Abstract
 Introduction
 Ryder and Petermann glaciers
 A twolayer model
 The dynamics in the hydraulic regime
 Conclusions
 Appendix A: Subglacial discharge and conservations relations
 Appendix B: Mathematical relationships for scenario 2
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Ryder and Petermann glaciers
 A twolayer model
 The dynamics in the hydraulic regime
 Conclusions
 Appendix A: Subglacial discharge and conservations relations
 Appendix B: Mathematical relationships for scenario 2
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References