1. Introduction

Conceptual models have clarified scientific understanding of the ocean’s deep stratification and meridional overturning circulation (MOC), while making predictions that may be tested via observations and numerical modeling. Our conceptual understanding has evolved, with two end points being the uniform vertical advective–diffusive balance assumed by Munk (1966) and the theory posited by Nikurashin and Vallis (2011, hereafter NV11), sketched in Fig. 1. This more recent view recognizes that the diapycnal diffusivity is an order of magnitude smaller than the estimate of Munk (1966) (e.g., Ledwell et al. 1993), and so transport is typically directed along isopycnals in the ocean interior (Lumpkin and Speer 2007). The purpose of this note is to highlight that, while the NV11 conceptual model illustrates the importance of the Southern Ocean in the setting, the global stratification, and MOC, it may change dramatically with the with improved representation of the dynamics in this region, especially at the Antarctic margins. As an example, modification of the surface boundary condition over the Southern Ocean dramatically alters the sensitivities of the deep stratification and MOC to diapycnal mixing.

A schematic of the residual-mean model of the deep overturning circulation. In the Southern Ocean channel (y < 0), the deep circulation is directed along isopycnals, and is constructed as a residual of a wind-driven mean circulation that tends to steepen isopycnals and the opposing “eddy” circulation that relaxes isopycnals by releasing potential energy through baroclinic instability (Karsten and Marshall 2002). The circulation is closed by buoyancy loss to the atmosphere, which supports a southward diabatic transport in the surface mixed layer, and by downward diapycnal diffusion of buoyancy in the northern basin (y > 0), which supports upwelling across isopycnal surfaces.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

View Full Size

A schematic of the residual-mean model of the deep overturning circulation. In the Southern Ocean channel (y < 0), the deep circulation is directed along isopycnals, and is constructed as a residual of a wind-driven mean circulation that tends to steepen isopycnals and the opposing “eddy” circulation that relaxes isopycnals by releasing potential energy through baroclinic instability (Karsten and Marshall 2002). The circulation is closed by buoyancy loss to the atmosphere, which supports a southward diabatic transport in the surface mixed layer, and by downward diapycnal diffusion of buoyancy in the northern basin (y > 0), which supports upwelling across isopycnal surfaces.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

A schematic of the residual-mean model of the deep overturning circulation. In the Southern Ocean channel (y < 0), the deep circulation is directed along isopycnals, and is constructed as a residual of a wind-driven mean circulation that tends to steepen isopycnals and the opposing “eddy” circulation that relaxes isopycnals by releasing potential energy through baroclinic instability (Karsten and Marshall 2002). The circulation is closed by buoyancy loss to the atmosphere, which supports a southward diabatic transport in the surface mixed layer, and by downward diapycnal diffusion of buoyancy in the northern basin (y > 0), which supports upwelling across isopycnal surfaces.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

The deep cell of the MOC is supplied by the outflow of Antarctic Bottom Water (AABW; Gordon 2009), upwells via the action of diapycnal mixing, and returns to the surface via outcropping isopycnal surfaces in the Southern Ocean (Lumpkin and Speer 2007). This circulation is of interest because it may control the ocean–atmosphere CO₂ partitioning over millennial time scales (Skinner et al. 2010). This has motivated a series of investigations of the sensitivity of the deep MOC, for example, to changes in surface forcing (Stewart and Thompson 2012; Meredith et al. 2012).

The role of conceptual models is to summarize in a clear and accessible way the mechanisms that control the deep stratification and MOC. The landmark study of Munk (1966) assumed the deep circulation to consist of uniform upwelling supported by downward diffusive fluxes, and established that a diapycnal diffusivity of κ ≈ 1 × 10⁻⁴ m² s⁻¹ was required to balance the export of AABW. This model dominated thinking for decades until direct measurements of mixing rates in the upper ocean established that the diffusivity is closer to κ ≈ 1 × 10⁻⁵ m² s⁻¹ (Ledwell et al. 1993). The implication that transport in the ocean interior is largely along isopycnals emphasizes the role of the Southern Ocean in closing water mass pathways, and is reflected in models like that of Gnanadesikan (1999), later extended to include the deep cell by Shakespeare and Hogg (2012). Ito and Marshall (2008) proposed a residual-mean model in which lateral mixing by mesoscale eddies and enhanced diapycnal mixing at depth constrain water mass transformation in the deep ocean, and thus the strength of the deep MOC.

The most recent dynamically consistent model of the deep cell is that of NV11, which agrees qualitatively with observations of the deep MOC (Lumpkin and Speer 2007). Their scalings for the sensitivity of the deep stratification and MOC to diapycnal mixing are corroborated by coarse-resolution simulations with a Gent and McWilliams (1990) parameterization, and agree qualitatively with the eddy-resolving simulations conducted by Munday et al. (2013). Here we identify the Southern Ocean surface buoyancy condition as an unexplored, but critical, component of this model that has received little attention and is poorly understood.

We demonstrate that a straightforward but plausible change in the surface boundary condition in a NV11-like configuration qualitatively changes the properties of the model. NV11 apply a fixed-buoyancy profile at the ocean surface, which is a suitable approximation for the strong restoring to atmospheric temperature acting along the Antarctic Circumpolar Current latitudes (Haney 1971). However, this boundary condition is not appropriate for the coastal regions of Antarctica in which AABW is formed. For example, brine rejection in nearshore polynyas may be more closely approximated as a fixed surface buoyancy flux (Chapman 1999; Wilchinsky and Feltham 2008). Here, we contrast the sensitivity of the model’s deep stratification and MOC to the diapycnal diffusivity κ in the limiting cases of a prescribed surface buoyancy versus a prescribed surface flux.

2. Scalings for deep stratification and overturning

a. Residual-mean model

We employ the conceptual model of NV11, illustrated in Fig. 1. For simplicity we restrict our attention to the deepest cell of the MOC, which is permitted to occupy the entire water column. We could restrict the deep cell to below 2–3-km depth, as in the real ocean, but our results would not change qualitatively.

Following NV11, the stratification and overturning in the Southern Ocean “channel” portion −l < y < 0 of our model domain are governed by

View Expanded

View Expanded

View Expanded

Here y and z are Cartesian coordinates corresponding to latitude and depth, and J(ψ, b) = ψ_yb_z − ψ_zb_y denotes the Jacobian operator. Equation (1a) balances advective transport of the time- and zonal-mean buoyancy b by a streamfunction ψ with diapycnal diffusion κ across buoyancy surfaces. The advecting streamfunction is a small residual in Eq. (1b) of a wind-driven mean overturning, proportional to the surface wind stress τ, and an opposing eddy overturning proportional to the isopycnal slope s = −b_y/b_z (Marshall and Radko 2003; Plumb and Ferrari 2005). We denote the reference density as ρ₀, the Coriolis parameter as f, and the isopycnal and diapycnal diffusivities as K and κ, respectively. We obtain the boundary condition in Eq. (1c) at y = 0 by assuming zero isopycnal slope in the northern basin (b_y ≡ 0), and integrating Eq. (1a) from y = 0 to the end of the basin y = L. To simplify the presentation we have taken τ, K, κ, and f to be constants; in the real ocean these quantities vary by orders of magnitude.

In the limit of a small residual overturning streamfunction, the two terms on the right-hand side of Eq. (1b) balance at leading order and prescribe a uniform isopycnal slope in the channel:

View Expanded

This approximation holds under the assumption that the diffusively driven residual overturning is much weaker than the wind-driven mean overturning:

View Expanded

Here we scaled ψ using the right-hand side of Eq. (1c), using h ~ l · s as a vertical scale for the stratification.

We also assume that the northern basin is much wider than the channel, as is the case in the real ocean,

View Expanded

Under this assumption we may neglect the right-hand side of Eq. (1a),

View Expanded

which states that the residual streamfunction is constant along isopycnals in the channel. It follows from Eqs. (5) and (2) that the buoyancy and streamfunction can be mapped between z = 0 and y = 0 via

View Expanded

b. Fixed surface buoyancy

NV11 studied the case of a prescribed surface buoyancy profile:

View Expanded

They argue that this is a suitable approximation in regions where temperature dominates the buoyancy variations and restoring to the atmosphere is fast.

From Eq. (6) it follows that the buoyancy in the northern basin is prescribed by the surface buoyancy profile:

View Expanded

The stratification in the northern basin is therefore prescribed by the surface buoyancy profile:

View Expanded

so it is independent of the diapycnal diffusivity κ.

We may similarly combine Eq. (7) with the northern boundary condition in Eq. (1c) to obtain an expression for the surface residual streamfunction:

View Expanded

At any fixed latitude y, the overturning ψ scales linearly with κ. The latitude of the surface streamfunction maximum ψ_max is determined by ∂_yψ(y_max, 0) = 0, so y_max is independent of κ. Therefore, the deep MOC strength should scale linearly with κ, in agreement with the findings of NV11. Intuitively this is because the surface buoyancy profile sets the deep stratification, which in turn sets the structure of the overturning streamfunction by Eq. (1c).

c. Fixed surface flux

We now consider the case of a fixed surface buoyancy flux at the ocean surface. This is arguably a more appropriate boundary condition where salinity plays an important role in setting buoyancy variations, like under ice.

Following Marshall and Radko (2003), we assume zero stratification in the mixed layer (b_z ≡ 0) and integrate Eq. (1a) from the base of the mixed layer (z = 0) to the ocean surface:

View Expanded

Here B is the downward buoyancy flux into the ocean surface. We will restrict our attention to the case of surface buoyancy loss, B ≤ 0. Using Eq. (6) along with Eq. (1c) we obtain the following:

View Expanded

Integrating once with respect to y yields the surface lateral stratification:

View Expanded

where M² = b_y|_y=−l,z=0 is the minimum lateral buoyancy gradient. Substituting Eq. (13) into Eq. (11) yields the surface streamfunction:

View Expanded

In Fig. 2, we contrast the solutions produced from a fixed surface buoyancy profile:

View Expanded

versus a fixed surface buoyancy flux:

View Expanded

Here, γ is the maximum of the second derivative of b_s(y) with respect to y. We could modify b_s(y) and B(y) such that the fixed-buoyancy and fixed-flux cases in Figs. 2a and 2b agree even more closely. We have chosen Eqs. (15) and (16) because they are simple, produce overturning streamfunctions of similar strength that vanish at y = −l and y = 0, and yield robust numerical results in section 3. Our parameter choices are listed in Table 1.

Plots of surface properties in the Southern Ocean channel: (a) the buoyancy profile b_s, (b) the downward buoyancy flux converted to an equivalent energy flux Q = ρ₀C_pB/αg, and (c) the streamfunction ψ as derived in section 2. In each case the diapycnal diffusivity is κ = 5 × 10⁻⁵ m² s⁻¹.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

View Full Size

Plots of surface properties in the Southern Ocean channel: (a) the buoyancy profile b_s, (b) the downward buoyancy flux converted to an equivalent energy flux Q = ρ₀C_pB/αg, and (c) the streamfunction ψ as derived in section 2. In each case the diapycnal diffusivity is κ = 5 × 10⁻⁵ m² s⁻¹.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

Plots of surface properties in the Southern Ocean channel: (a) the buoyancy profile b_s, (b) the downward buoyancy flux converted to an equivalent energy flux Q = ρ₀C_pB/αg, and (c) the streamfunction ψ as derived in section 2. In each case the diapycnal diffusivity is κ = 5 × 10⁻⁵ m² s⁻¹.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

Table 1.

List of parameters used in our analytical scaling and numerical solutions.

View Table

As the buoyancy flux is always directed out of the ocean (B ≤ 0), it follows from Eq. (12) that the stratification maximum lies at y = z = 0. Substituting y = 0 into Eq. (13) and dividing by −s we obtain the maximum stratification:

View Expanded

Intuitively, the net diffusive buoyancy flux convergence in the northern basin must balance the net buoyancy loss at the ocean surface.

As in the prescribed buoyancy case, at any fixed latitude y the streamfunction in Eq. (14) scales linearly with diffusivity κ when κ is sufficiently small. However, the form of Eq. (14) implies that the latitude at which the streamfunction maximum lies itself depends on the diffusivity, y_max = y_max(κ) (here the “maximum” of ψ refers to its most negative value). For example, a maximum of ψ always exists close to y = −l: for a surface buoyancy flux that satisfies B(y) → C(y + l)^p as y → −l for some constant C, the streamfunction in Eq. (14) is maximized at

View Expanded

so there exists an maximum of ψ that approaches the southern boundary y = −l as κ → 0. Substituting Eq. (18) into Eq. (14), we find that the overturning scales as ψ_max ~ κ^p/(p+1) for small κ. In general ψ may have local maxima at other latitudes, but it follows from Eq. (14) that any other maximum satisfies ψ_max ~ κ ≪ κ^p/(p+1) as κ → 0, so for sufficiently small κ the maximum close to y = −l is the global maximum. Thus, the dependence of the deep MOC strength on diapycnal mixing depends not only upon whether the surface buoyancy fluxes are fixed, but also upon the latitudinal structure of the fluxes. We illustrate this point in Fig. 3, which shows the sensitivity of the MOC to diapycnal mixing for different buoyancy flux profiles.

Dependence of the max overturning ψ_max = max_y,z|ψ| on the diapycnal diffusivity κ for different surface buoyancy flux profiles B(y). The total buoyancy loss at the ocean surface is the same for each profile. The parameters otherwise match those listed in Table 1.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

View Full Size

Dependence of the max overturning ψ_max = max_y,z|ψ| on the diapycnal diffusivity κ for different surface buoyancy flux profiles B(y). The total buoyancy loss at the ocean surface is the same for each profile. The parameters otherwise match those listed in Table 1.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

Dependence of the max overturning ψ_max = max_y,z|ψ| on the diapycnal diffusivity κ for different surface buoyancy flux profiles B(y). The total buoyancy loss at the ocean surface is the same for each profile. The parameters otherwise match those listed in Table 1.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

3. Numerical sensitivity to diapycnal mixing

a. Numerical configuration

The scalings derived in section 2 are only valid for small ε ≪ 1, or equivalently for small κ. To obtain a more general picture of the dependence of the deep MOC on diapycnal mixing, we now solve the residual-mean equations in Eqs. (1a)–(1c) and (7), (11) numerically. We obtain a steady solution of Eq. (1a) by integrating its time-dependent equivalent to steady state:

View Expanded

Our numerical approach follows that of Stewart and Thompson (2013).

We apply boundary conditions of no-normal-flow (ψ = 0) at y = −l and z = −H. An intuitive additional boundary condition at the ocean bed is that there should be no-normal buoyancy flux (i.e., b_z = 0 at z = −H). However, this boundary condition is ill posed when combined with the assumption of flat isopycnals in the northern basin. Integrating Eq. (1c) vertically from z = −H, it may be shown that if the stratification vanishes at the bottom boundary, b_z|_z=−H = 0, then it must also vanish throughout the water column, b_z ≡ 0. Numerical simulations employing this boundary condition develop a spuriously large streamfunction adjacent to the ocean bed, and do not converge under refinement of the numerical grid. We circumvent this problem by instead prescribing the stratification at the ocean bed,

View Expanded

In reality isopycnals should outcrop from the ocean bed throughout the northern basin, and dense water should cross these isopycnals laterally in a bottom boundary layer whose thickness scales with the diapycnal diffusivity κ. Our results may therefore be sensitive to the choice of bottom boundary condition, particularly when κ is large and the bottom boundary layer is thick. This sensitivity needs to be quantified using a three-dimensional general circulation model, but such a study is beyond the scope of this paper. Here we focus on the sensitivity of the deep circulation to the surface boundary condition, which remains regardless of the treatment of the bottom boundary.

For the case of fixed surface fluxes in Eq. (11), we may simply choose , which corresponds to M² = 5 × 10⁻⁹ s⁻² using the parameters listed in Table 1. For the case of fixed surface buoyancy, we prescribe to balance the steady-state buoyancy budget:

View Expanded

where square brackets denote that the expressions within should be evaluated between limits, indicated by the subscript and superscript following the brackets. We approximate the surface terms in Eq. (21) using a uniform slope in Eq. (2) and the analytical streamfunction in Eq. (10), using continuity of b_z at y = 0 to determine the surface stratification in the northern basin 0 < y < L:

View Expanded

Note that the buoyancy profile in Eq. (15) ensures that the streamfunction ψ vanishes at y = 0, so there is no advective buoyancy exchange with the mixed layer in the northern basin.

Rather than impose Eqs. (7) or (11) directly on z = 0, we have found that the stability and accuracy of the solution are improved by enforcing these conditions over a layer of finite depth H_s = 100 m close to the surface. For the fixed-flux case in Eq. (11) we impose a downward buoyancy flux B that decreases linearly in magnitude from z = 0 to z = −H_s. For the fixed-buoyancy case in Eq. (7) we replace the fixed-buoyancy flux B with a relaxative flux W_s(b − b_s). The surface piston velocity W_s = 5 × 10⁻⁵ m s⁻¹ restores the surface buoyancy with a time scale of around 3 weeks.

b. Sensitivity to diapycnal mixing

In Fig. 4 we plot the numerically computed stratification and overturning in the fixed-buoyancy and fixed-flux cases, for κ ranging over two orders of magnitude. Over this range the overturning at y = 0 varies by 0.68 m² s⁻¹ in the fixed-buoyancy case, and only by 0.26 m² s⁻¹ in the fixed-flux case. This is consistent with our scalings in section 2, which suggest that the overturning should scale linearly with κ in the fixed-buoyancy case, and as κ^1/2 in the fixed-flux case (see Fig. 3). In contrast the stratification exhibits almost no change in the fixed-buoyancy case, but varies by an order of magnitude in the fixed-flux case.

Plots of the numerically computed buoyancy b (10⁻² m s⁻²; solid contours) and overturning streamfunction ψ (m² s⁻¹; dashed contours) for diapycnal diffusivities of (a),(d) κ = 5 × 10⁻⁶ m² s⁻¹; (b),(e) κ = 5 × 10⁻⁵ m² s⁻¹; and (c),(f) κ = 5 × 10⁻⁴ m² s⁻¹. In (a)–(c), the fixed surface buoyancy profile in Eq. (15) is employed, while in (d)–(f) the fixed surface flux condition in Eq. (16) is employed. All other parameters are listed in Table 1. Note the larger buoyancy contour intervals in (d). The dashed contour interval in each panel is one-tenth of the max overturning strength.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

View Full Size

Plots of the numerically computed buoyancy b (10⁻² m s⁻²; solid contours) and overturning streamfunction ψ (m² s⁻¹; dashed contours) for diapycnal diffusivities of (a),(d) κ = 5 × 10⁻⁶ m² s⁻¹; (b),(e) κ = 5 × 10⁻⁵ m² s⁻¹; and (c),(f) κ = 5 × 10⁻⁴ m² s⁻¹. In (a)–(c), the fixed surface buoyancy profile in Eq. (15) is employed, while in (d)–(f) the fixed surface flux condition in Eq. (16) is employed. All other parameters are listed in Table 1. Note the larger buoyancy contour intervals in (d). The dashed contour interval in each panel is one-tenth of the max overturning strength.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

Plots of the numerically computed buoyancy b (10⁻² m s⁻²; solid contours) and overturning streamfunction ψ (m² s⁻¹; dashed contours) for diapycnal diffusivities of (a),(d) κ = 5 × 10⁻⁶ m² s⁻¹; (b),(e) κ = 5 × 10⁻⁵ m² s⁻¹; and (c),(f) κ = 5 × 10⁻⁴ m² s⁻¹. In (a)–(c), the fixed surface buoyancy profile in Eq. (15) is employed, while in (d)–(f) the fixed surface flux condition in Eq. (16) is employed. All other parameters are listed in Table 1. Note the larger buoyancy contour intervals in (d). The dashed contour interval in each panel is one-tenth of the max overturning strength.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

In Fig. 5 we compare the numerically calculated maximum overturning and stratification with the scalings derived in section 2, for a range of κ between 10⁻⁶ and 10⁻³ m² s⁻¹. For κ larger than this, the overturning cell begins to interact with the bottom boundary. The scalings are quantitatively accurate as long as κ is sufficiently small, as one would expect from the small-ε approximation in Eq. (3). The discrepancy between our scaling and the numerical results in Fig. 5b is due to diffusive modification of the overturning streamfunction ψ along isopycnals in the channel. Some discrepancy is present for all values of κ, no matter how small, because the maximum stratification scales as , and so the diffusive term on the right-hand side of Eq. (1a) is nonzero as κ → 0.

Sensitivity of (a),(b) the max overturning strength ψ_max and (c),(d) the max stratification in the northern basin to the rate of diapycnal mixing. In (a),(c) the fixed surface buoyancy profile in Eq. (15) is employed, while in (b),(d) the fixed surface flux condition in Eq. (16) is employed. For our numerical calculations we use the max overturning at the boundary between the channel and the northern basin, ψ_max = max_y=0|ψ|.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

View Full Size

Sensitivity of (a),(b) the max overturning strength ψ_max and (c),(d) the max stratification in the northern basin to the rate of diapycnal mixing. In (a),(c) the fixed surface buoyancy profile in Eq. (15) is employed, while in (b),(d) the fixed surface flux condition in Eq. (16) is employed. For our numerical calculations we use the max overturning at the boundary between the channel and the northern basin, ψ_max = max_y=0|ψ|.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

Sensitivity of (a),(b) the max overturning strength ψ_max and (c),(d) the max stratification in the northern basin to the rate of diapycnal mixing. In (a),(c) the fixed surface buoyancy profile in Eq. (15) is employed, while in (b),(d) the fixed surface flux condition in Eq. (16) is employed. For our numerical calculations we use the max overturning at the boundary between the channel and the northern basin, ψ_max = max_y=0|ψ|.

Citation: Journal of Physical Oceanography 44, 3; 10.1175/JPO-D-13-0206.1

• Download Figure
• Download figure as PowerPoint slide

In the fixed-buoyancy case the overturning appears to undergo transition to a large-κ regime in which ψ_max ~ κ^1/2, as in NV11. The stratification also undergoes a transition, scaling approximately as for large κ. This transition is due to the fact that lateral transports in the channel are supported by eddy thickness fluxes (i.e., that the vertical isopycnal spacing must narrow in the direction of the residual transport). As κ and ψ increase, this effect begins to modify the stratification in the northern basin, which in turn impacts the overturning circulation. In the fixed-flux case, both ψ_max and agree qualitatively with our scalings from section 2. For very strong diapycnal mixing κ the stratification becomes almost uniform and equal to . This fixes the surface buoyancy gradient, which in turn fixes the overturning circulation via Eq. (11).

4. Discussion

We have extended the conceptual model of NV11 to demonstrate that the sensitivities of the deep stratification and MOC ψ_max to diapycnal mixing are themselves sensitive to the surface boundary condition. If the surface buoyancy is prescribed then we recover the results of NV11, in which is effectively prescribed and ψ_max scales as κ¹ for small κ and κ^1/2 for large κ. In contrast, if the surface fluxes are prescribed via Eq. (16) then scales as κ⁻¹ for small κ, so varying κ over the range observed in the ocean changes the stratification by orders of magnitude. The scaling for the overturning also changes, with ψ_max scaling as κ^1/2 for small κ and as κ⁰ for large κ. However, in the fixed-flux case the scaling of ψ_max with κ also depends on the latitudinal structure of the flux, as shown in Fig. 3. This complicates prediction of the deep MOC sensitivity to diapycnal mixing, as surface buoyancy fluxes are poorly constrained (Cerovecki et al. 2011).

For the purpose of illustration we have focused on one particular property of the NV11 model: its sensitivity to diapycnal mixing. A complete study would covary many of the model parameters (e.g., the wind stress τ, eddy diffusivity K, and diapycnal diffusivity κ). NV11 argued that their model suggested that the deep ocean overturning is sensitive to diapycnal mixing rates, while the stratification is set by the surface boundary conditions in the Southern Ocean independently of diapycnal mixing. This result overturned the traditional view pioneered by Munk (1966), who argued that the ocean stratification is set by diapycnal mixing alone. We have shown that the result of NV11 depends on the choice of surface boundary condition. NV11 used a fixed-buoyancy boundary condition, appropriate for waters that outcrop equatorward of the sea ice line and are strongly restored to atmospheric temperature. Switching to fixed-flux boundary conditions more appropriate to describe the waters that outcrop around the Antarctic continent, we find that both the stratification and overturning become sensitive to diapycnal mixing, but at rates depending on the specific latitudinal distribution of the fluxes. In reality the ocean surface is likely subject to a combination of these two boundary conditions. M. Nikurashin (2013, personal communication) plans to extend our work in an upcoming paper and quantify the relative importance of these two forcings in the present ocean.