## 1. Introduction

The Antarctic Circumpolar Current (ACC) deflects steadily poleward while circumnavigating Antarctica, returning abruptly equatorward along the Argentine continental slope. This behavior is consistent with Ekman upwelling over the Southern Ocean north of Cape Horn, causing fluid columns to stretch and drift poleward, with an equatorward return flow confined to a narrow western boundary current. This is the paradigm of “Sverdrup balance” applied to an ocean gyre (Sverdrup 1947; Stommel 1957), or, in this case, to the Southern Hemisphere “supergyre” (Ridgway and Dunn 2007). Extending this paradigm to the ACC, including the latitudes of Drake Passage, is challenged by a lack of continental barriers, without which there can be no equatorward boundary current to balance the otherwise poleward flow induced by the wind stress curl at these latitudes.

One proposed extension is an “eddy Sverdrup balance,” where Sverdrup balance is augmented by the eddy-induced downwelling, added alongside the Eulerian-mean upwelling that it largely negates. In an idealized reduced-gravity channel model, Marshall et al. (2016) witnessed such a three-way balance involving the Eulerian-mean northward velocity, the Eulerian-mean vertical velocity, and the eddy-induced vertical velocity. This result may be particular to the reduced-gravity formulation, but it cannot be correct in the real ACC: the eddy-induced velocity is only relevant in the context of fluxing tracers or momentum, arising as it does from correlations between the instantaneous velocity and tracers or momentum. In Sverdrup balance, no such eddy correlations arise: the planetary vorticity is time invariant, so its time-mean advection is simply its advection by the time-mean velocity, not by any residual (time-mean plus eddy-induced) velocity.

Instead, in the Drake Passage latitudes, one must consider the more general form of Sverdrup balance, namely, the barotropic vorticity equation, wherein all torques are retained, not just the wind stress curl. Which additional torque is most important? It is well known that bottom form stress removes the vast majority of zonal momentum input by the westerly winds in the Drake Passage latitudes (Munk and Palmén 1951; Masich et al. 2015). Therefore, the analog of bottom form stress in the barotropic vorticity equation, namely, the bottom pressure torque, is of prime importance in these latitudes (Hughes and De Cuevas 2001; Lu and Stammer 2004). Indeed, the ACC extends to the seafloor and must navigate considerable topographic relief, inducing a vertical bottom velocity that may be an order of magnitude larger than the Ekman upwelling velocity, locally.

Moving from barotropic dynamics to layer-wise dynamics, it is also well known that eddy form stress (also called interfacial form stress) is largely responsible for transferring eastward momentum imparted by the surface winds downward to the seafloor, where it is removed by bottom form stress (Johnson and Bryden 1989). The same logic as above should apply: the curl of the eddy form stress divergence should be a dominant term, locally and in the net, in the layer-wise residual vorticity budget in the Drake Passage latitudes.

Now, dividing the residual vorticity budget by the mean layer thickness (the vertical separation between buoyancy surfaces) yields a residual potential vorticity (PV) equation. It stands to reason that a dominant term in the residual PV budget is the curl of the eddy form stress divergence divided by the mean layer thickness. With such a dominant forcing term present in the Drake Passage latitudes, we should not expect the mean PV to be materially conserved, following the residual velocity, in the eddying Southern Ocean. Rather, varying mean PV along streamlines of the ACC is part-and-parcel of how the ocean navigates these circumpolar latitudes.

The variation of PV along streamlines in the ACC has been studied, most notably by Gille (1997). Our results agree with Gille’s key result, that mean PV rises dramatically along streamlines transiting Drake Passage and the Argentine slope, balanced by PV decreasing along streamlines elsewhere. However, different model data and numerical methods alter the finer details and interpretation. In particular, Gille attributes PV decreasing along streamlines to large-scale patterns in Ekman upwelling. However, subsequent advances in understanding the residual overturning circulation of the Southern Ocean (Danabasoglu et al. 1994; Gent et al. 1995; Marshall 1997; Marshall and Radko 2003) indicate that eddy-induced downwelling largely compensates the Ekman upwelling of the Eulerian-mean Deacon cell. Instead, we argue that eddies are the agents of such PV forcing along deep ACC streamlines, and this eddy PV forcing is concomitant with downward eddy form stress and poleward eddy buoyancy fluxes. In summary, this paper attempts to update Gille’s (1997) important findings in light of theoretical and modeling development over the past two decades.

A word on averaging formalisms is in order, since eddy forcings are defined by the choice of averaging operator. From a theoretical perspective, the thickness-weighted average (TWA) formalism (Young 2012) has emerged as the ideal averaging technique, since it effectively captures the residual overturning of the Southern Ocean. Under TWA, averages are taken following buoyancy surfaces as they heave vertically, and are weighted by the layer thickness to account for the bolus transport. Unfortunately, we found the TWA PV budget using 5-day Eulerian-averaged outputs from a *z*-level model cannot be closed. The primary numerical issue, detailed in appendix C, is that the buoyancy and continuity equations are extremely delicate, and do not survive the transformation off their native *z*-level grids. Nonetheless, the mean PV, whether from TWA or Eulerian averages, exhibits the same qualitative structure along ACC streamlines (appendix B). Hence, we employ Eulerian averaging, which enables closure of the PV budget in a *z*-level ocean model.

Section 2 pedagogically describes why PV cannot be conserved in the eddying Southern Ocean, using large-scale approximations. These approximations are removed in section 3, which describes the Ertel PV appropriate for a numerical ocean model, with numerical methods detailed in appendix A. The structure and forcing of PV along geostrophic streamlines in the ACC are studied in sections 4 and 5, before concluding in section 6.

## 2. From vorticity to potential vorticity

Eddies play a key role in ocean dynamics through eddy form stress, eddy buoyancy fluxes, eddy-induced velocity, and PV diffusion. Since all of these are different manifestations of the same physical process, there is great scope for double-counting the effects of eddies. This confusion is evidenced by the discussion in Warren et al. (1996); see subsequent replies by Hughes (1997) and Olbers (1998). As such, the averaged equations should be derived carefully and systematically. This section provides a pedagogical derivation, deriving the linear vorticity budget, with or without some common approximations, then linking this to the large-scale PV budget.

### a. Notation and setup

For simplicity (of presentation only), we use a Cartesian coordinate system with eastward, northward, and upward coordinates (*x*, *y*, *z*) and orthonormal basis vectors *t* the time coordinate. Let an overbar denote the arithmetic time-mean—which when applied a function of Cartesian coordinates (*x*, *y*, *z*, *t*) gives the Eulerian time-mean. Let *p* be the pressure, *ρ* the in situ density, *ρ*_{0}(*z*) a reference density profile, *ρ _{c}* the Boussinesq reference density,

*g*the (constant) gravitational acceleration,

*a*is

*a*at constant height is denoted

### b. Vorticity balance

**v**= 0 by the Coriolis parameter

*f*and applying the product rule yields one form of the vorticity balance,

*β*=

*f*and

_{y}*f*= 0. Only the ageostrophic flow contributes to ∇ · (

_{z}*f*

**u**), but we temporarily leave this expressed with the full flow. No eddy terms arise when time-averaging (1), hence

**v**

^{*}yields the eddy vorticity balance,

This equation, postulated as relevant to the ACC, is merely two balances combined as one. The mean flow must satisfy its own vorticity balance (2) independently of that for the eddying flow (3), and vice versa.

Whereas both the mean and eddy-induced velocities pierce mean buoyancy surfaces, the residual velocity is closely aligned with these surfaces; this will help to frame the PV equation. But first, the above vorticity budgets are cast in a more familiar form.

### c. Approximate vorticity balance

**is the vertical stress. Neglected terms (lateral viscous stresses, tendency, and nonlinear terms) can be carried alongside**

*τ*

*τ**if desired. The geostrophic velocity*

_{z}*f*

**u**

*) = 0. The ageostrophic velocity*

_{g}**u**

_{ag}=

**u**−

**u**

*satisfies*

_{g}*b*(see section 3a) is assumed to be quasi-materially conserved,

*b*′ and

**u**′ are fluctuations from their Eulerian time-means

*ξ*′, are related to fluctuations in local buoyancy

*b*′ by

*ξ*′ = −

*b*′/

*N*

^{2}, and (ii) fluctuations of the horizontal velocity

**u**′ are geostrophic so that

^{1}Thus,

### d. Approximate residual vorticity balance

Under nonacceleration conditions, the residual overturning circulation vanishes, i.e.,

The stretching terms (*fw _{z}*) will be annihilated when forming a PV equation, but the stress divergences will remain, translating into potentially significant PV forcings. Before moving to the PV budget, we first derive and discuss Sverdrup balance.

### e. Barotropic vorticity balance and Sverdrup balance

*z*= −

*H*to the sea surface at

*z*= 0 (assuming a rigid lid for simplicity of presentation), apply the Leibniz integral rule to the integrated pressure gradient, and finally cross differentiate to get the (approximate) barotropic vorticity balance,

*p*=

_{b}*p*|

_{z}_{=−}

*is the bottom pressure, and*

_{H}The wind stress curl is generally negative in the Southern Ocean, inducing a southward Sverdrup flow. But the zonal integral of *βV* must be zero to maintain mass balance (ignoring evaporation, precipitation, and river runoff). Thus, other terms in (12) must be important—particularly the bottom pressure torque, ∇* _{z}* × (

*p*∇

_{b}*H*). In the oceanic gyres, the Sverdrup flow is also generally single-signed, so the bottom pressure torque—which is nonlocally linked to the wind stress curl (Stewart et al. 2021)—establishes a western boundary current whose meridional mass flux balances that of the interior Sverdrup flow. The Drake Passage latitudes lack boundaries to support such a boundary current, so the bottom pressure torque becomes geographically distributed; it alternates in sign and locally swamps the wind stress curl, inducing a depth-integrated meridional flow with a complicated spatial structure, though zonally integrating to zero (Hughes and De Cuevas 2001). As the bottom pressure torque is linked to the topographic structure and largely balances the residual between the wind stress curl and the net meridional flow (i.e., the failure of Sverdrup balance), its rich structure is also seen in non-eddy-resolving models (Lu and Stammer 2004).

This patchy nature of the barotropic vorticity balance is concomitant with a similarly patchy nature of the vorticity balance on buoyancy surfaces. The eddy form stress *τ*^{*} that transfers momentum down from the sea surface to the seafloor exhibits a similarly complicated spatial structure, typically enhanced downstream of major topographic obstacles (MacCready and Rhines 2001; Tansley and Marshall 2001). This eddy form stress, as we show next, appears in the PV budget, thereby inducing large and spatially complicated material changes in PV on interior buoyancy surfaces.

### f. Large-scale potential vorticity budget

The preceding residual vorticity balance equations were intrinsically the sum of mean and eddy balances, because time-invariant quantities like *f* and *β* multiplied the velocity components. However, the large-scale PV, *fb _{z}*, is time varying, so the residual equation for its material derivative cannot be similarly split.

*b*as a function only of

*ρ*and

*z*annihilates the bracketed (solenoidal) term. The eddy part is simplified similarly, but using (8) to give

The LHS is the material derivative of the large-scale, mean PV following the residual velocity. The RHS provides its material forcings, namely, the tracer forcing gradient projected onto the large-scale vorticity, and the vorticity forcing projected onto the tracer gradient.

### g. Discussion

Since the eddy form stress divergence is a large forcing of the residual momentum budget (and the eddy buoyancy flux divergence is a large forcing of the Eulerian-mean buoyancy budget), we can reasonably expect that the corresponding *τ*^{*} (or

Mean streamlines in the ACC’s northern extent, in particular, exhibit first-order variations of mean PV (Gille 1997, and our section 4). Moving downstream along these streamlines, PV rises markedly along the South American continental slope, balanced by a general decrease elsewhere. Gille (1997) suggested that wind forcing was primarily responsible for the latter effect, arguing that Ekman upwelling squeezes isopycnals toward the surface, thereby decreasing layer thickness and decreasing PV (toward larger negative values). However, the direct effect of wind on the mean PV, captured by the vertical stress divergence term in (20) or (21), is only significant in the upper few tens of meters. Beneath that, the wind stress only influences the PV budget indirectly, by impacting PV advection and eddy form stress, which are the dominant terms in the steady-state version of (20)—or for (21), by impacting PV advection and eddy buoyancy fluxes. We therefore argue that eddies, not winds, are the primary drivers of PV variations along ACC streamlines.

Beyond driving variations in general, can we also explain the sign of that eddy PV forcing, and is it consistent with the robust feature, mentioned above, that PV decreases moving downstream along northern ACC streamlines?

For simplicity, suppose the eddy form stress is zonal, *τ*^{*} > 0, transporting eastward momentum downward, whereas *τ*^{*} ≈ 0 north and south of the ACC (where momentum input by the winds can be handled by leaning against topography—a type of bottom form stress). The vertical structure of *τ*^{*} has a subsurface maximum, increasing from near zero in the Ekman layer where viscous stress handles the downward momentum transport (

A consistent conclusion is reached from the Eulerian perspective, wherein the northern flank of the ACC experiences Ekman convergence and a mean downwelling, largely compensated for by an eddy-induced upwelling. This eddy-induced upwelling advects high PV (low in magnitude) upward, increasing PV at the point in question. The reverse holds on the southern flank, with eddy-induced downwelling causing PV to decrease.^{2}

The prediction is thus that eddy PV forcing increases PV moving downstream along northern ACC streamlines, away from boundaries. Yet, the opposite is observed: PV decreases. By investigating a numerical model, we shall see that the ACC adjusts its geography to avoid this apparent contradiction.

## 3. Ertel PV in a numerical ocean model

We use the Southern Ocean State Estimate (SOSE; Mazloff et al. 2010) spanning 2005–10 to diagnose the time-mean Ertel PV and its budget along time-mean geostrophic streamlines. The necessary theory is derived here, using the same approximations as made in SOSE: the Boussinesq approximation, hydrostatic balance, the traditional approximation, and its partner the thin-film approximation.

### a. Buoyancy and its surfaces

Buoyancy is generally defined as (−*g*/*ρ _{c}*) multiplied by a density variable that would ideally (i) be quasi-materially conserved, (ii) have isosurfaces that are everywhere tangent to the neutral tangent plane (McDougall 1987), and (iii) possess an exact geostrophic streamfunction on its isosurfaces (McDougall 1989; Stanley 2019b). Unfortunately, property (ii) is prevented by nonlinearities in the equation of state (McDougall and Jackett 1988).

Many density variables exist that approximate these properties, but only specific volume anomaly (Montgomery 1937) and orthobaric density (de Szoeke et al. 2000) satisfy property (iii), which is equivalent to annihilating the solenoidal term forcing an Ertel PV defined using that density variable (de Szoeke 2000)—essential for our purposes. We use specific volume anomaly, for computational simplicity and the ready availability of its partial derivatives. Future analyses might use (modified) topobaric surfaces (Stanley 2019a,b); however, these are 2D surfaces not 3D density fields, hindering the PV budget’s calculation on the model’s native *z*-level grid.

*ρ*is determined by the equation of state

*S*is the salinity,

*θ*the potential temperature, and

*p*the pressure. Under the Boussinesq approximation, however, the equation of state takes a hydrostatic pressure determined from the height

*z*:

*S*

_{0}and

*θ*

_{0}are to-be-defined constants.

Let *ξ*(*x*, *y*, *b*, *t*) denote height as a function of buoyancy in each water column, i.e., the inverse function of *b*(*x*, *y*, *z*, *t*) for a given (*x*, *y*, *t*). Let *ξ*^{#}(*x*, *y*, *b*) denote height as a function of Eulerian time-mean buoyancy, i.e., the inverse function of *x*, *y*). The sharp notation, borrowed from the transformed residual mean (TRM; McDougall and McIntosh 2001) and TWA (Young 2012) theories, indicates a quantity constructed from mean fields, rather than a mean of some quantity. Indeed,

### b. Geostrophic streamfunction

*z*

_{1}is a constant and

*p*

_{1}(

*x*,

*y*,

*t*) =

*p*(

*x*,

*y*,

*z*

_{1},

*t*). To evaluate

*m*on the chosen buoyancy surface, simply set

*z*=

*ξ*(

*x*,

*y*,

*b*

_{0},

*t*), which also replaces

*b*(

*x*,

*y*,

*z*,

*t*) by the constant

*b*

_{0}.

### c. Instantaneous budgets

#### 1) Momentum budget

*ζ*=

*υ*−

_{x}*u*is the relative vorticity’s vertical component,

_{y}*K*= (

*u*

^{2}+

*υ*

^{2})/2 is the specific kinetic energy,

*B*=

*K*+

*ϕ*is the Bernoulli potential, and

**= 0 in the thin-film approximation (Dellar 2011).**

*ω*#### 2) Buoyancy budget

*S*,

*θ*, and

*z*(the latter trivial),

*b*,

We use ∂ notation to denote partial derivatives with respect to thermodynamic variables *S*, *θ*, and *z*; we use subscript notation to denote partial derivatives with respect to coordinate variables *x*, *y*, *z*, *t*. This helps distinguish the dual roles of *z* as a thermodynamic variable and a coordinate, which are related by

#### 3) Ertel PV budget

*Q*is derived via

**= 0 and ∇ × ∇**

*ω**b*= 0 to obtain the divergence of

The baroclinic production term *b* is a function of *ρ* and *z*.

### d. Time-mean: The Ertel–Euler PV

*z*-level model like SOSE, the natural averaging operator is the Eulerian time-mean. We define the “Ertel–Euler PV”—the Ertel PV formed from Eulerian averaged variables—as

*Q*

^{#}is that its budget automatically closes upon closing the Eulerian time-mean momentum and buoyancy budgets. These are readily found by averaging (25) to get

*Q*

^{#}:

*b*.

The material derivative in (39) uses the Eulerian-mean velocity, but PV is best conserved following the residual velocity (section 2f). Ideally, we would employ the TRM-II theory of McDougall and McIntosh (2001), which is designed to approximate buoyancy-following averages using Eulerian averages. However, TRM-II requires the variance of a quasi-conservative density variable, which is lacking from SOSE’s standard outputs—understandably so, since there is no universally best quasi-conservative density variable (section 3a). Computing this variance from SOSE’s 5-day Eulerian averages may produce similar errors in TRM-II as arise when using TWA with SOSE’s 5-day Eulerian averages (appendix C), so we have not pursued this further. Nonetheless, much can still be learned from the Ertel–Euler PV approach.

## 4. The structure of PV along streamlines

We now assess the structure of the Ertel–Euler PV along time-mean geostrophic streamlines in the ACC on the buoyancy surface defined by *x*_{0} = 133.25°E, *y*_{0} = 54.0416°S, and *z*_{0} = −700 m, which further defines *σ*_{0} = 27.46 kg m^{−3}, near the bottom of Antarctic Intermediate Water.

Figure 1 shows the geography of geostrophic streamlines on the *Q*^{#} along these streamlines (Fig. 1a) and *Q*^{#} itself (Fig. 1b). There are large *Q*^{#} variations along southerly streamlines where the surface shoals to near the base of the mixed layer, though *Q*^{#} is relatively well conserved through southern Drake Passage. In contrast, northern streamlines exhibit a dramatic rise in *Q*^{#} in the Drake Passage region, in agreement with Gille (1997), and a similar but smaller rise along Campbell Plateau; these are balanced by a steady decrease of *Q*^{#} elsewhere.

These *Q*^{#} variations are seen more clearly in Fig. 2, which shows *Q*^{#} (in black) along four streamlines: A, the northernmost circumpolar streamline; B, another northerly streamline, passing slightly further off the Argentine slope; C, the most southerly streamline that passes north of Kerguelen Plateau; and D, a southerly streamline passing between Kerguelen Plateau and Heard Island. Along streamline A, *Q*^{#} rises by ∼50% through Drake Passage (∼60°W), rises by ∼15% along Campbell Plateau (∼170°E), and drops steadily everywhere else. The steady rate of PV decrease is remarkable, given the richly heterogeneous eddy field of PV forcings. Along streamline C, *Q*^{#} decreases between Kerguelen Plateau (∼70°E) and Drake Passage, and actually increases slightly on the approach to Kerguelen to pass north of the plateau. The structure of *Q*^{#} along streamline B is a mixture of that along streamlines A and C, with *Q*^{#} roughly constant between Drake Passage and Kerguelen Plateau. The largest fluctuations of *Q*^{#} along streamline D occur upstream and in the vicinity of Kerguelen Plateau. None of these streamlines exhibit the PV structure found by Gille (1997), who shows that PV is essentially constant in the Atlantic and Indian sectors, and that all PV decreases occur between Campbell Plateau and Drake Passage (her Fig. 7a). This is likely due to differences between SOSE and the Semtner–Chervin model that Gille (1997) analyzed (appendix B).

The Ertel–Euler PV *Q*^{#} (thick black) along four streamlines A–D corresponding to the first, second, fourth, and seventh streamlines in Fig. 1 from north to south. The latitude and height at which they cross 133.25°E are indicated at the bottom left of each panel. Thin lines show the thin-film PV, *f* varies (red), stratification

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

The Ertel–Euler PV *Q*^{#} (thick black) along four streamlines A–D corresponding to the first, second, fourth, and seventh streamlines in Fig. 1 from north to south. The latitude and height at which they cross 133.25°E are indicated at the bottom left of each panel. Thin lines show the thin-film PV, *f* varies (red), stratification

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

The Ertel–Euler PV *Q*^{#} (thick black) along four streamlines A–D corresponding to the first, second, fourth, and seventh streamlines in Fig. 1 from north to south. The latitude and height at which they cross 133.25°E are indicated at the bottom left of each panel. Thin lines show the thin-film PV, *f* varies (red), stratification

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Figure 2 also shows how the thin-film PV approximation, *f*, *f* contribution is fairly stable—as expected to remain north of Kerguelen. Downstream of Kerguelen on streamline A, PV decreases are mostly achieved by decreasing *f* as streamlines slowly migrate southward to meet Drake Passage, whereas *f* variations as the streamline wends southward toward Drake Passage, while *f* variations on short scales. Streamline C shows large compensation between *f* variations, while streamline D remains consistently far south and is dominated by

Gille (1997) proposed a mechanism wherein Ekman upwelling squeezes isopycnals toward the surface, thereby increasing |*Q ^{#}*| and decreasing

*Q*. However, Figs. 1a and 2a provide evidence against this hypothesis, as the fastest drop of

^{#}*Q*on streamline A due to stratification changes occurs between ∼10° and ∼45°E, which is generally a region of Ekman downwelling, not upwelling. We shall discuss this further in the next section, after establishing a framework for the material forcing of

^{#}*Q*along streamlines.

^{#}## 5. The forcing of Ertel–Euler PV along streamlines

### a. Theory

*Q*

^{#}over a distance

*s*along such a streamline is

*is the lateral gradient operator acting within the time-mean buoyancy surface.*

_{b}*Φ = ∇*

_{b}*Φ + Φ*

_{z}*∇*

_{z}*. This yields*

_{b}z### b. Numerics

Computing the RHS of (42) and (45) requires care (the LHS was described in section 3). Regarding the integrated geostrophic advection of PV, the RHS of (42), the mean geostrophic velocity is staggered, with its zonal and meridional components living on the *υ* grid (south edge) and *u* grid (west edge), respectively. This directly enables computing *Q*^{#} on the surface. Then, the PV forcing *ζ* grid to the streamline’s contour points, then integrated as follows. The geostrophic speed *υ* grid and *u* grid, respectively. The distance, d*s*, is computed as the great circle distance between consecutive points on the streamline. Finally,

The result, however, cannot exactly match the LHS of (42), the true PV. To see why, consider a streamline passing south and east of the stencil’s central point: the north and west *Q*^{#} values never enter the interpolation for *Q*^{#}(*s*) (the LHS), but they do affect

*Q*^{#} (black) and its integrated forcing (colors) along the same four streamlines as in Fig. 2, starting from 40°W. The integrated forcings derive from ageostrophic advection (purple), buoyancy drift (dotted blue), eddy buoyancy fluxes (solid blue), buoyancy forcing (dashed blue), *w*^{#} advection (green), momentum drift (not shown, being negligible), eddy momentum fluxes (solid red), and momentum forcing (dashed red), in order from top to bottom then left to right on the RHS of (45). Also shown is the integrated geostrophic advection of *Q*^{#} [light gray, RHS of (42)], which should match *Q*^{#} (black), apart from advection errors caused by numerical discretization.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

*Q*^{#} (black) and its integrated forcing (colors) along the same four streamlines as in Fig. 2, starting from 40°W. The integrated forcings derive from ageostrophic advection (purple), buoyancy drift (dotted blue), eddy buoyancy fluxes (solid blue), buoyancy forcing (dashed blue), *w*^{#} advection (green), momentum drift (not shown, being negligible), eddy momentum fluxes (solid red), and momentum forcing (dashed red), in order from top to bottom then left to right on the RHS of (45). Also shown is the integrated geostrophic advection of *Q*^{#} [light gray, RHS of (42)], which should match *Q*^{#} (black), apart from advection errors caused by numerical discretization.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

*Q*^{#} (black) and its integrated forcing (colors) along the same four streamlines as in Fig. 2, starting from 40°W. The integrated forcings derive from ageostrophic advection (purple), buoyancy drift (dotted blue), eddy buoyancy fluxes (solid blue), buoyancy forcing (dashed blue), *w*^{#} advection (green), momentum drift (not shown, being negligible), eddy momentum fluxes (solid red), and momentum forcing (dashed red), in order from top to bottom then left to right on the RHS of (45). Also shown is the integrated geostrophic advection of *Q*^{#} [light gray, RHS of (42)], which should match *Q*^{#} (black), apart from advection errors caused by numerical discretization.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

For the RHS of (45), most PV forcing terms *G* live natively on the same grid as *Q*^{#}, so are handled exactly as above. The two exceptions are the ageostrophic and *w*^{#} advections, discussed next.

The zonal ageostrophic velocity *υ* grid, obtained by averaging *u* onto the *T* grid, then linearly interpolating vertically onto the surface, then averaging horizontally onto the *υ* grid. The meridional component is handled similarly. Then,

*w*

^{#}advection integral directly using (43) and (40b), which yields

^{3}

Almost indistinguishable results were obtained using other integration methods, e.g., assuming *G* and *G* over the time required for the geostrophic flow to travel the distance *ds*.

### c. Results

Figure 3 presents the results of the above integrations, forming the basis for the following discussion. Figures 4–6 provide supporting information, mapping *Q*^{#} forcings across the Southern Ocean, around Drake Passage, and on meridional sections, respectively.

Maps of selected PV forcing terms (s^{−4}) deriving from (a) the buoyancy tendency, (b) the *w*^{#} advection, (c) the eddy buoyancy flux, and (d) the eddy momentum flux. Their signs are as in (45), so *Q*^{#} increases following the geostrophic flow through red regions. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 2/3° zonal standard deviation and 1/3° meridional standard deviation. Light gray indicates where the mean buoyancy surface has incropped or outcropped. The four streamlines in Fig. 2 and Fig. 3 are shown in green.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Maps of selected PV forcing terms (s^{−4}) deriving from (a) the buoyancy tendency, (b) the *w*^{#} advection, (c) the eddy buoyancy flux, and (d) the eddy momentum flux. Their signs are as in (45), so *Q*^{#} increases following the geostrophic flow through red regions. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 2/3° zonal standard deviation and 1/3° meridional standard deviation. Light gray indicates where the mean buoyancy surface has incropped or outcropped. The four streamlines in Fig. 2 and Fig. 3 are shown in green.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Maps of selected PV forcing terms (s^{−4}) deriving from (a) the buoyancy tendency, (b) the *w*^{#} advection, (c) the eddy buoyancy flux, and (d) the eddy momentum flux. Their signs are as in (45), so *Q*^{#} increases following the geostrophic flow through red regions. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 2/3° zonal standard deviation and 1/3° meridional standard deviation. Light gray indicates where the mean buoyancy surface has incropped or outcropped. The four streamlines in Fig. 2 and Fig. 3 are shown in green.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. 4, but zoomed into the Drake Passage region and showing all PV forcings, deriving from (a) buoyancy tendency, (b) buoyancy forcing, (c) eddy buoyancy flux, (d) *w*^{#} advection, (e) momentum forcing, (f) eddy momentum flux, (g) ageostophic advection, and (h) geostrophic advection. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 1/3° zonal standard deviation and 1/6° meridional standard deviation.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. 4, but zoomed into the Drake Passage region and showing all PV forcings, deriving from (a) buoyancy tendency, (b) buoyancy forcing, (c) eddy buoyancy flux, (d) *w*^{#} advection, (e) momentum forcing, (f) eddy momentum flux, (g) ageostophic advection, and (h) geostrophic advection. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 1/3° zonal standard deviation and 1/6° meridional standard deviation.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. 4, but zoomed into the Drake Passage region and showing all PV forcings, deriving from (a) buoyancy tendency, (b) buoyancy forcing, (c) eddy buoyancy flux, (d) *w*^{#} advection, (e) momentum forcing, (f) eddy momentum flux, (g) ageostophic advection, and (h) geostrophic advection. For presentation, these maps are smoothed by a mean-preserving Gaussian filter with 1/3° zonal standard deviation and 1/6° meridional standard deviation.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Meridional sections of (top) the eddy buoyancy flux and (bottom) its contribution to the *Q*^{#} budget, centered at (left) 55°E and (right) 90°W and zonally averaged over 0.5°. Green dots indicate where the nine streamlines from Fig. 1 intersect the meridional section. Buoyancy surfaces are contoured every 10^{−3} m s^{−2}, up to the buoyancy value at the surface of the northernmost streamline.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Meridional sections of (top) the eddy buoyancy flux and (bottom) its contribution to the *Q*^{#} budget, centered at (left) 55°E and (right) 90°W and zonally averaged over 0.5°. Green dots indicate where the nine streamlines from Fig. 1 intersect the meridional section. Buoyancy surfaces are contoured every 10^{−3} m s^{−2}, up to the buoyancy value at the surface of the northernmost streamline.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Meridional sections of (top) the eddy buoyancy flux and (bottom) its contribution to the *Q*^{#} budget, centered at (left) 55°E and (right) 90°W and zonally averaged over 0.5°. Green dots indicate where the nine streamlines from Fig. 1 intersect the meridional section. Buoyancy surfaces are contoured every 10^{−3} m s^{−2}, up to the buoyancy value at the surface of the northernmost streamline.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Integrating the geostrophic *Q*^{#} advection along streamlines gives a reconstruction of *Q*^{#}(*s*) (gray), which closely matches the true *Q*^{#}(*s*), capturing its small-scale features and large-scale structure—except for near topography such as Kerguelen Plateau (∼70°E), Campbell Plateau (∼170°E), and Drake Passage (∼60°W), where one extreme *Q*^{#} value can spoil the reconstruction. Generally though, advection errors caused by numerical discretization are small.

On the southernmost streamline D, *Q*^{#}(*s*) is fairly well conserved relative to its integrated forcings. The dominant forcings are the *w*^{#} advection (green) that increases *Q*^{#}, and eddy buoyancy fluxes (solid blue) that decrease *Q*^{#}. This streamline is consistently in a region of Ekman upwelling (Fig. 1a), which increases *Q*^{#} by advecting deeper, higher *Q*^{#} (lower |*Q*^{#}|) water up onto the buoyancy surface. Though this heuristic argument is helpful, the relevant vertical advection is not the Ekman upwelling velocity, but rather is *w*^{#} which is heavily modified by the eddy field. Nonetheless, *w*^{#} is more positive than negative in this region (not shown), so *w*^{#} advection forces *Q*^{#} to increase along this and other southerly streamlines. This is opposite to the mechanism proposed by Gille (1997) that Ekman upwelling squeezes isopycnals toward the surface and decreases *Q*^{#} (increases |*Q*^{#}|). Gille’s (1997) proposed mechanism is rather Lagrangian in nature: *Q*^{#} supposedly decreases as isopycnals rise because it is measured rising with the isopycnals. Since our (and Gille’s) analysis framework is Eulerian, so must be our interpretation: *Q*^{#} at a fixed point will increase as upwelling causes isopycnals to rise past this point, which now resides in less stratified water originating from below. Returning to the dominant balance, the *w*^{#} upwelling is largely counterbalanced by an eddy-induced downwelling, bringing low *Q*^{#} water from above down onto the buoyancy surface, and predominantly captured by the eddy buoyancy flux divergence in our Eulerian framework.

On streamline C, *w*^{#} advection causes *Q*^{#} to decrease, at least between Drake Passage and Kerguelen. Here, there is notable Ekman downwelling even on this reasonably southerly streamline, heuristically consistent with *w*^{#} < 0 and *Q*^{#} decreasing. As above, eddies act to oppose the *w*^{#} forcing, though here eddy momentum fluxes (solid red) also contribute strongly. This may be surprising, given that *Q*^{#} is very well approximated by the thin-film PV,

Downstream of Kerguelen, streamline C is more clearly in a region of Ekman upwelling, heuristically consistent with *w*^{#} forcing *Q*^{#} to increase, although interpretation in this region is challenged by the fact that buoyancy drift (dotted blue) becomes a major driver of *Q*^{#} loss downstream of Kerguelen. Indeed, buoyancy drift does not generally dominate the *Q*^{#} budget in any single water column (except in the far south, where the buoyancy surface is near the sea surface), but it is single-signed over vast regions, so can integrate to a considerable forcing (Fig. 4a). Finally, though momentum forcing (dashed red) is generally insignificant as a driver of *Q*^{#}, close encounters with topography can lead to massive and rapid *Q*^{#} forcing by frictional drag, as at ∼70°E where streamline C passes just north of Kerguelen Plateau, and sidewall drag exerts a negative torque on the fluid, decreasing *Q*^{#}. The frictional forcing of PV is a robust feature along topography (Fig. 5e); it is not due to one malicious grid point, though it appears so from the zoomed-out view of Fig. 3c.^{4}

Along streamline B, eddy forcing—of both buoyancy and momentum—force *Q*^{#} to decrease in the Indian sector from ∼20° to ∼170°E, opposed by *w*^{#} advection. This region contains the Dynamical Subtropical Front, which (unlike the Subtropical Front) is also a deep water mass boundary (Graham and De Boer 2013). This front is seen clearly as a dipolar structure in the Indian sector of Figs. 4b and 4c, where eddy buoyancy fluxes cause *Q*^{#} to decrease on the southern side of the dipole, and to increase on the northern side—as predicted in section 2g—while *w*^{#} advection does the opposite. In fact, *Q*^{#} actually decreases northward here (Fig. 1b), suggesting that the *w*^{#} advection establishes this northward *Q*^{#} decrease, while eddy forcing erodes it, homogenizing PV as eddies that flux PV downgradient do. Indeed, this region exhibits relatively homogeneous PV (Stanley et al. 2020, their Fig. 8). Along this front, mean buoyancy surfaces rise steeply poleward (Figs. 6a,c), lifted in the south and depressed in the north by Ekman suction/pumping (Fig. 1a). Through baroclinic instability, mesoscale eddies oppose this, fluxing buoyancy poleward across the front, creating a dipolar structure to the eddy buoyancy flux divergence. Moreover, this divergence increases in magnitude moving up the water column until ∼300-m depth (Fig. 6a), broadly consistent with the vertical structure of the fastest growing baroclinic modes (Smith 2007). Thus, on the southern flank of this front at middepths, *Q*^{#} to decrease. Put another way, mesoscale eddies tend to flatten isopycnals, largely by means of an eddy-induced velocity **v**^{*} that is downward on the southern flank, advecting low *Q*^{#} water onto the buoyancy surface and materially decreasing *Q*^{#} (with the reverse occurring on the northern flank). Streamline B transits through the south flank of the Dynamical Subtropical Front (Figs. 4b,c and 6a,c), and so encounters eddy forcing that materially decreases *Q*^{#}.

This balance for streamline B reverses in the Pacific sector, where eddy forcing increases *Q*^{#} and *w*^{#} advection decreases *Q*^{#}. While the dipolar structure of the Dynamical Subtropical Front is less clear in the Pacific sector, there is evidently still a net signal of eddy-induced upwelling that increases *Q*^{#}, even in the eastern Pacific where eddy buoyancy forcing of *Q*^{#} appears more noisy and eddy-like (Fig. 4c). This is typical of the Drake Passage latitudes, with Ekman pumping counterbalanced by eddy-induced upwelling. Streamline B is northerly enough that it escapes this regime in the Indian sector, but must inhabit this regime in the Pacific sector in order to pass through Drake Passage.

Where streamline B flows northward along Campbell Plateau and the Argentine slope, eddy buoyancy and momentum forcing both cause *Q*^{#} to increase. Usually, the ACC’s eastward flow is sufficiently strong that eddies are advected eastward (Hughes 1995; Hughes et al. 1998; Klocker and Marshall 2014). With the flow deflected northward here, eddies can propagate west [see Fig. 3.1 of Zajaczkovski (2017), or Fig. 1c of Klocker and Marshall (2014)], impinging on the continental slope where they break—largely by transferring energy into high wavenumber vertical modes (Dewar and Hogg 2010)—and form an “eddy graveyard” (Zhai et al. 2010). These breaking eddies strongly force the PV, generally causing *Q*^{#} to increase and enabling the northward flow in a self-supporting cycle.

Streamline A is the northernmost circumpolar streamline, passing very close to the Argentine slope. Here, the dramatic *Q*^{#} gain is driven overwhelmingly by momentum forcing, where the northward flow along a no-slip western boundary induces a positive torque on the fluid, increasing *f* + *ζ*, and *Q*^{#}. Again, this frictional forcing is a robust signal, with large and generally positive *Q*^{#} forcings along the Argentine slope (Fig. 5e).

Away from Drake Passage and Campbell Plateau, the steady decrease of *Q*^{#} along streamline A is driven by a complicated combination of eddy buoyancy and momentum fluxes and *w*^{#} advection, not to mention buoyancy drift. These signals are subtler than for streamline B because streamline A threads between the large dipolar structure of *Q*^{#} forcings by eddy buoyancy fluxes and *w*^{#} advection in the Indian sector (Figs. 4b,c and 6c). Noting that the net *Q*^{#} eddy forcing in the Pacific sector is not as strong and coherent as in the Indian sector (Figs. 4b,c and 6c,d), we hypothesize that streamline A cannot pass any further north than it does while remaining circumpolar: if streamline A inhabited the northern flank of the Subtropical Dynamical Front, eddy buoyancy fluxes would strongly force *Q*^{#} to increase there, opposite to the *Q*^{#} decrease that is required to counterbalance the rise of *Q*^{#} that robustly occurs along the Argentine slope. The usual coarse-grained view of the Southern Ocean’s residual overturning circulation involves mean (Ekman) upwelling and eddy-induced downwelling in the southern reaches of the ACC, and the reverse in the northern reaches. However, finer-resolution inspection reveals that isopycnals are steeper in the Dynamical Subtropical Front than in the rest of the ACC further south (Fig. 6a), leading to an eddy buoyancy flux convergence and eddy-induced downwelling on the southern flank of this front, even though this lies in the northern reaches of the ACC. If we define the northern limit of the ACC to be the northernmost circumpolar (geostrophic) streamline, then the northern limit of the ACC—at least in the Indian sector—is actually the southern flank of the Dynamical Subtropical Front, where mean upwelling and eddy-induced downwelling occur.

For all streamlines, diapycnal forcing of *Q*^{#} (dashed blue) is generally small, in agreement with Gille (1997), though not insignificant for streamline D, it being relatively shallow. Our choice of in situ density anomaly to define buoyancy and our numerics ensure the solenoidal term is precisely zero, whereas Gille (1997) used potential density and found the solenoidal term a major PV forcing (her Fig. 7b). We repeated our numerical methods using potential density and found the solenoidal PV forcing becomes nonzero but remains negligible, meaning the large solenoidal forcing Gille (1997) found is due to her numerics, not inherent issues with potential density. Gille (1997) also found buoyancy and momentum forcing is a dominant PV forcing even away from topographic boundaries, whereas we find these terms are negligible in the ocean interior; this difference is likely an artifact of the lower-resolution Semtner–Chervin model analyzed by Gille (1997).

Ageostrophic advection is generally insignificant at forcing *Q*^{#}, except near topographic obstacles. The dramatic rise of *Q*^{#} along streamline A at Drake Passage, for example, is driven partly by ageostrophic advection (purple). Indeed, some geostrophic streamlines intersect South America, therefore requiring relatively strong ageostrophic flow to navigate Drake Passage. High *Q*^{#} water originating north of the ACC is ageostrophically advected onto streamline A, increasing *Q*^{#}. Real fluid parcels do not follow the mean streamlines, geostrophic or otherwise. If a fluid parcel transiting Drake Passage is better at conserving its PV than fluid following the streamlines, its *Q*^{#} will rise less at Drake Passage and consequently will migrate southward across streamlines. This is consistent with the residual meridional overturning circulation taking water, particularly in the Atlantic sector, across the ACC from north to south at these middepths, as water “spirals in” toward Antarctica (Tamsitt et al. 2017). However, this effect of ageostrophic advection forcing *Q*^{#} to increase at Drake Passage is not robust, as other streamlines exhibit the opposite tendency (Fig. 7a). Indeed, the ageostrophic forcing of *Q*^{#} in the Drake Passage region is rich and complex (Fig. 5g). In the southern ACC however, ageostrophic forcing of *Q*^{#} appears to be more placid, creating a small but steady rise of *Q*^{#} along streamline D, again consistent with poleward ageostrophic flow on this surface.

As in Fig. 3, but for streamlines on different buoyancy surfaces, all intersecting 133.25°E at −700 m.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. 3, but for streamlines on different buoyancy surfaces, all intersecting 133.25°E at −700 m.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. 3, but for streamlines on different buoyancy surfaces, all intersecting 133.25°E at −700 m.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

To examine the sensitivity of these results to the particular depth of the surface, we repeat the analysis along four new streamlines—labeled A′, B′, C′, and D′—that intersect 133.25°E at −700 m height, regardless of latitude (Fig. 7). Each such streamline exists on a different buoyancy surface, so we choose reference *S*_{0} and *θ*_{0} local to where each streamline intersects (133.25°E, −700 m), and recalculate all quantities. The dominant balance along streamline D′ is the same as that for streamline D, despite the latter being far shallower (generally at 150 to 250 m depth). Likewise, the results along streamline C′ are qualitatively unchanged, as expected since streamline C intersected (133.25°E, −54.87°S) at −577 m, not far from −700 m. Further analyses of *Q*^{#} forcing along streamlines A′ and B′ are withheld, owing to the dominant effect of buoyancy drift around 700 m depth in the ACC during the SOSE record. This is not particular to Eulerian averaging, but also arises under TWA averaging (Fig. C1). This dominant buoyancy drift is also seen in gridded Argo data (not shown, but expected since SOSE incorporates Argo).

## 6. Conclusions

The wind stress curl over the Southern Ocean lifts up deep, dense water through Ekman suction and drives poleward Sverdrup flow that, in the latitudes north of Cape Horn, is balanced by an equatorward-flowing western boundary current along the Argentine slope. However, this is not a classic subpolar gyre, since the boundary current is fed not by westward drift, but by the circumpolar current itself (Stommel 1957). In the Drake Passage latitudes where there are no continental barriers, the poleward Sverdrup flow is instead balanced by a richly heterogeneous but net equatorward flow supported by bottom pressure torque. That is, the bottom pressure torque is a dominant term in the barotropic vorticity budget, balancing wind stress curl in the net. This is linked, as Hughes and De Cuevas (2001) made clear, to bottom form drag (the curl of which is the bottom pressure torque) being the dominant sink of easterly momentum input by the wind stress in the barotropic momentum budget.

Moving from barotropic to layerwise dynamics, eddies become of prime importance. The downward transfer of mean momentum from the surface to the bottom is overwhelmingly achieved by eddy form stress on sloping isopycnals, enabled by hundreds of mesoscale eddies that populate the ACC. It is well known that this process is accompanied by a mean poleward buoyancy transport (Johnson and Bryden 1989). We have shown that these related eddy processes are also accompanied by a material forcing of mean potential vorticity, driven by the curl of the eddy form stress. This is a dominant term in the mean potential vorticity budget, just as the eddy form stress is a dominant term in the mean layerwise momentum budget—analogous to the aforementioned link established by Hughes and De Cuevas (2001).

For a jet like the ACC (or its filaments) in the open ocean, this eddy PV forcing ideally takes a dipolar structure, acting at middepths to increase PV on the northern flank, and decrease PV on the southern flank. Meanwhile, as the ACC transits along the Argentine slope, PV is materially increased (toward zero), as shown by Gille (1997) and verified here. This material PV increase is driven by drag spinning down the current or, farther offshore, breaking eddies having the same effect.

Together, these two regimes create an apparent contradiction for fluid parcels recirculating along the northern flank of the ACC, as their PV would continually increase. One might argue that fluid parcels do not recirculate in the ACC, but rather migrate across closed streamlines. Given the material PV increase, however, such a migration would be equatorward, in contrast to the poleward migration at these middepths required to balance the water mass transformation into Antarctic Bottom Water. Regardless, such migration is slight, as ageostrophic forcing of mean PV along mean geostrophic streamlines is almost negligible compared to eddy or near-boundary PV forcing. Instead, the ACC resolves this contradiction by merging with the subtropical gyres to the north, creating an eddy-rich jet for which the northernmost Antarctic Circumpolar Current’s streamlines inhabit the jet’s southern flank, where PV decreases. The PV decrease thus accumulated in the open ocean is able to balance the PV increase as circumpolar streamlines transit along the Argentine slope.

The length of the ACC’s northward excursion along the Argentine slope may therefore be dynamically controlled by this requirement for the ACC to merge with the subtropical gyres. Since the length of this northward ACC excursion (as the Malvinas Current) determines the latitude of the Brazil/Malvinas Confluence (BMC), this suggests a new mechanism by which remote factors can influence the BMC latitude: more (or less) eddy-driven PV loss along ACC streamlines in the open ocean must be balanced by longer (or shorter) northward ACC excursions along the Argentine slope. Though a complete theory for the BMC latitude is still lacking, this mechanism may play a role alongside other factors such as the relative transports of the Brazil and Malvinas Currents (Matano 1993) and the latitude of the maximum wind stress curl (Lumpkin and Garzoli 2011). Such dynamic constraints may control the mean BMC latitude as well as its fluctuations on decadal time scales—for which there appears to be a recent southward trend (Lumpkin and Garzoli 2011; Combes and Matano 2014)—while fluctuations of the BMC latitude on shorter time scales appear to be driven by Rossby waves (Bodnariuk et al. 2021). In these matters, care must be taken to distinguish between the latitude of separation from the continental slope and the more southerly latitude at which the BMC settles offshore; between these latitudes the Malvinas Current is overshooting, enabling further PV gain along the continental slope. This hypothesis, by which remote eddy PV forcing on the northern flank of the ACC influences the BMC latitude, remains speculative and requires further research. One implication of this hypothesis is that the BMC latitude may be different between coarse-resolution and eddy-resolving models. An intriguing and possibly related question regards why the eddy PV forcing appears, in SOSE, to most coherently exhibit the aforementioned dipolar structure with axis parallel to the jet in the Indian sector, and least in the Atlantic sector.

These two regimes of the mean PV structure along mean streamlines of the northern ACC—decreasing in the open ocean and increasing downstream of Drake Passage—are robust signals, regardless of whether the averaging is Eulerian or semi-Lagrangian (following buoyancy surfaces as in TWA).

In the Eulerian perspective, the open-ocean PV increase is framed as the near but not exact cancelation of the mean Ekman overturning circulation by an eddy-induced overturning cell. On the northern flank of a jet, the mean and eddy-induced overturning circulations respectively advect low (large and negative) and high (near zero) PV from above and below, acting to materially decrease or increase PV. The reverse holds on the southern flank. The eddy PV forcing slightly outcompetes the mean PV advection, producing a remarkably steady decrease of PV along mean northern circumpolar streamlines in the open ocean, particularly in the Indian sector.

Under a TWA analysis, wherein the residual circulation of the Southern Ocean is naturally captured, the two dominant but counterbalancing PV forcings of the Eulerian PV budget would, we expect, coalesce into a single dominant PV forcing due to eddy form stress, at least in the ocean interior. Close to boundaries, other TWA PV forcings such as due to friction and ageostrophic advection will be significant, as in the Eulerian PV budget.

As Gille (1997) pointed out, this nonconservation of mean PV along mean streamlines bears on many theoretical models of the ACC. For instance, the equivalent barotropic structure of Killworth (1992) relies on a functional relationship between PV and the mass-weighted streamfunction on layers, which is broken by this nonconservation of mean PV. Going one step further, Marshall (1995) assumed that PV is uniform on layers. This assumption is one (probably major) reason why his modeled ACC does not veer northward downstream of Drake Passage, as the real ACC does: with PV uniform, his flow can only move northward by shrinking the layer thickness, upon which mass conservation puts strong limits. Such models are intended only to capture the first-order dynamics of the ACC, and hence can be quite tractable and illuminating. However, given the dominant role of eddy PV forcing and that mean PV can vary by ∼50% along a mean streamline—not to mention variations across streamlines (Stanley et al. 2020)—we argue that PV-conserving flow should not be considered as the ACC’s first-order dynamics. Rather, nonconservation of PV is a first-order effect.

Several multibasin idealized models have recently coupled global oceanic processes—notably convective downwelling in the North Atlantic and diffusive upwelling predominantly in the Pacific—through the Southern Ocean with simplified ACC dynamics (Thompson et al. 2016; Jones and Cessi 2016). To geostrophically support a mass transfer of intermediate waters from the Pacific to Atlantic (and vice versa of deep waters), these models require a difference of PV between the Atlantic and Pacific sectors. As this PV difference is established by eddy and boundary PV forcings, further work with such models that brings PV forcing to the forefront of ACC dynamics may prove enlightening.

A similar result can be obtained more carefully from the TRM-II theory of McDougall and McIntosh (2001). Also note the equivalence between *E ^{u}* in (C13): simply use

*p*=

_{x}*ϕ*=

_{x}The major difference between (20) and (21) is that **v**^{†} in (20) is nearly aligned with mean buoyancy surfaces, so the PV variation along mean buoyancy surfaces is clearly driven by the eddy form stress term. In contrast, *ϖ*^{†}, the residual between the mean PV advection and the eddy buoyancy flux PV forcing is the eddy-induced PV advection.

To be clear, this does not mean that the

The discretization of appendix A was crucial to obtain accurate analysis this close to topography.

This second-order centered advection scheme for buoyancy is simpler than SOSE’s third-order upwind advection scheme for salinity and potential temperature, but is sufficient to illustrate the grid problems.

## Acknowledgments.

We thank Ed Doddridge, Ryan Holmes, Chris Hughes, Matt Mazloff, Trevor McDougall, Jan Zika, and two reviewers for helpful discussions. GJS acknowledges partial support from the Clarendon Fund and the Canadian Alumni Scholarship at Linacre College, Oxford, and from the Australian Research Council though Grant FL150100090. DPM acknowledges the partial support of the Natural Environment Research Council, NE/R000999/1.

## Data availability statement.

SOSE output is available at http://sose.ucsd.edu.

## APPENDIX A

### Numerical Methods

Careful numerical methods are required to study the structure and forcing of PV along streamlines. SOSE uses an Arakawa C-grid (Arakawa and Lamb 1977): *u*, *υ*, and *w* live on the west, south, and upper faces of the tracer cell, respectively, while vertical vorticity (*f* and *ζ*) lives on the southwest edge of the tracer cell.

#### a. Buoyancy

*ρ*from the known pressure

_{k}*p*and height

_{k}*z*at each model level

_{k}*k*, starting with

*η*is the sea surface height, and continuing for

*k*> 1 with

Next, we select the reference values *S*_{0} and *θ*_{0}, which should be chosen close to *S* and *θ* found in one’s study region to make the buoyancy surface close to neutral. We simply choose a reference location (*x*_{0}, *y*_{0}, *z*_{0}) and define *ρ* replaced by

This buoyancy variable can exhibit overturns, which we simply and effectively deal with by sorting *b* is monotonic in *z* by choosing their reference density profile as the vertical integral of a virtual compressibility, itself chosen as a function of height to be the minimum compressibility in their domain at each height. Their buoyancy variable therefore has *b _{z}* ≥

*N*

^{2}(with equality at one point for each height, where their surface is neutral) leading to surfaces that are systematically inclined away from the neutral tangent plane and toward geopotentials, and hence are less neutral than our prescription. Essentially, the contribution of vertical variations of

*S*and

*θ*to vertical variations of in situ density—which Uchida et al. (2021) missed in their Eq. (4)—means the virtual compressibility should be allowed to exceed the true compressibility in roughly half the ocean, such that

*b*≈

_{z}*N*

^{2}.

The isovalue of our buoyancy surface is *z*. Likewise, the height of our surface, *ξ*^{#}(*x*, *y*, *b*_{0}), is obtained by linear interpolation of *z* as a function of *x*, *y*) to *b*_{0}. Hence, *ξ*^{#}(*x*_{0}, *y*_{0}, *b*_{0}) = *z*_{0}. Note *b*_{0} ≈ 0 m s^{−2} but not exactly because of nonlinearities in

#### b. Geostrophic streamfunction

*ξ*

^{#}:

Numerically, we use trapezoidal integration in (41), consistent with obtaining *ξ*^{#} by linear interpolation of *b*_{0}. We use *z*_{1} = −5 m, the height of the shallowest tracer cell center.

We desire the streamline of *x*_{0}, *y*_{0}), but in general the level set of *x*_{0}, *y*_{0})—this is guaranteed to be circumpolar—then remove all filled data points from this contour. This yields a circumpolar streamline of *x*_{0}, *y*_{0}) that may consist of multiple disjoint contours.

#### c. Ertel–Euler PV budget

To close the time-mean momentum budget, we use model outputs for

The time-mean buoyancy budget cannot be determined purely from time-mean quantities, because of nonlinearity in *b _{t}*,

**v**· ∇

*b*, and

Our numerical discretization is carefully designed to satisfy the vector identities *ζ* = *υ _{x}* −

*u*, is computed using Stokes’ theorem to exchange an area integral of

_{y}*ζ*for a line integral of the circulation around the perimeter of the vertical vorticity cell centered on the southwest corner of the tracer cell (matching the model’s numerics). We extend this analogously to compute the horizontal components of the vorticity. We average

*b*onto these three vorticity grids, multiply by the vorticity components, integrate across the six faces of the vorticity cell centered on the upper southwest corner of the tracer cell, and divide by the vorticity cell’s volume to obtain

The above procedure calculates (∇ × **a**) · ∇*T* as ∇ · (*T*∇ × **a**) for any “staggered” vector field **a** (having components on the *u*, *υ*, and *w* grids) and any scalar field (tracer) *T*. Hence, the same numerics that calculate *Q*^{#} naturally calculate its forcing terms in (39). For example, **a** = **v** yields the eddy buoyancy forcing of PV, **a** is trivially augmented by *Q*^{#} will also have valid *Q*^{#} forcings. This discretization ensures the solenoidal term

We now have 3D data for *Q*^{#} and its forcings on the upper southwest corner of the tracer cells. We average *ξ*^{#}(*x*, *y*, *b*_{0}) four ways onto the southwest corner of the horizontal tracer grid, onto which we vertically interpolate *Q*^{#} and its forcings from the height of the *w* grid. This interpolation is linear, so the PV budget (39) that holds at each grid point also holds on the surface. Finally, *Q*^{#} and its forcings are horizontally interpolated to the contour points of each geostrophic streamline; these points live directly between two adjacent grid points of the geostrophic streamfunction, so this interpolation is linear (not bilinear).

## APPENDIX B

### Euler–Ertel and Thickness-Weighted Average PV

Here, we verify that the structure of the Ertel–Euler PV along mean geostrophic streamlines is a robust signature, insofar as it qualitatively matches with results using alternative definitions of PV and geostrophic streamlines.

#### a. Euler–Ertel PV

##### 1) Definition

The time-mean PV used by Gille (1997) is

*σ*

_{1}(potential density referenced to 1000 dbar) is used instead of

*b*from (22).

##### 2) Geostrophic streamfunction

No exact geostrophic streamfunction exists for potential density surfaces (McDougall 1989). Gille (1997) used the Zhang and Hogg (1992) generalization of the Montgomery (1937) potential, which references the surface’s height to its mean value, thereby reducing errors in estimating the geostrophic velocity caused by variations of *b* on the *σ*_{1} surface. Instead, we use the orthobaric Montgomery potential (Stanley 2019b), which references the surface’s height to a value that depends on the local value of *b* on the *σ*_{1} surface. This further generalizes the Zhang and Hogg (1992) streamfunction and further reduces errors.

##### 3) Numerics

In contrast to *Q*^{#}, we assemble *σ*_{1} from 5-day means of *S* and *θ*, from which we build *σ*_{1})* _{z}* by finite difference on the

*w*grid, then average back to the tracer grid. We calculate

*ζ*from 5-day means of

*u*and

*υ*as in appendix A, section c, then average

*ζ*onto the tracer grid. This gives a 5-day mean

^{−3}), which intersects the reference point (

*x*

_{0},

*y*

_{0},

*z*

_{0}) = (133.25°E, 54.0416°S, −700 m). Finally, we linearly interpolate

#### b. Thickness-weighted average PV

##### 1) Definition

*x*,

*y*,

*z*,

*t*) by

*ξ*is a dependent variable. The layer thickness is

*F*is

*F*and

*h*are functions of

Note the TWA does not commute with partial derivatives:

##### 2) Geostrophic streamfunction

*z*

_{1}—implicitly defined by

##### 3) Numerics

To calculate Π^{#}, we use piecewise cubic Hermite interpolating polynomials (PCHIPs) because they preserve monotonicity and are continuously differentiable at knots. Specifically, *z* as a function of *b* at each *u*, averaged onto the tracer grid, onto *b* = *b*_{0}, and similarly for *υ*. We repeat this for each 5-day average to assemble *m*, themselves obtained from (B5) with the integral evaluated analytically.

For reference, the time-mean height of the *b* = *b*_{0} surface at (*x*_{0}, *y*_{0}) = (133.25°E, 54.0416°S) is

#### c. Relation between three time-mean PVs

*Q*

^{#}somewhat resembles the TWA PV Π

^{#}, at least more so than does the Euler–Ertel PV. Consider the thin-film limit, which neglects horizontal contributions to the PV:

*Q*

^{#}resembles the TWA of the thin-film PV (

*f*+

*ζ*)

*b*, albeit with two differences. First, Π

_{z}^{#}involves division by

*b*, but (B6) involves multiplication by the arithmetic time-mean of

_{z}*b*=

_{z}*h*

^{−1}. Second and more importantly, for

*Q*

^{#}the time-mean is Eulerian, sampling from a constant depth, whereas for Π

^{#}the time-mean is quasi-Lagrangian, sampling from a buoyancy surface as it heaves vertically.

#### d. Results

Figure B1 verifies that *Q*^{#}, Π^{#}, and *σ*_{1} and

Comparison of different definitions for PV along geostrophic streamlines: the Ertel–Euler PV *Q*^{#} from (35) on Montgomery streamlines on ^{#} from (B4) on Montgomery streamlines on the surface *b* = *b*_{0} (dashed red); the Euler–Ertel PV *σ*_{1} surface (dotted blue). The *σ*_{1} surfaces intersect (133.25°E, 54.04°S) at −700 m, while the mean height of the *b* surface is very near −700 m here. The streamlines cross 133.25°E at the indicated latitudes, as in Fig. 2.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Comparison of different definitions for PV along geostrophic streamlines: the Ertel–Euler PV *Q*^{#} from (35) on Montgomery streamlines on ^{#} from (B4) on Montgomery streamlines on the surface *b* = *b*_{0} (dashed red); the Euler–Ertel PV *σ*_{1} surface (dotted blue). The *σ*_{1} surfaces intersect (133.25°E, 54.04°S) at −700 m, while the mean height of the *b* surface is very near −700 m here. The streamlines cross 133.25°E at the indicated latitudes, as in Fig. 2.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Comparison of different definitions for PV along geostrophic streamlines: the Ertel–Euler PV *Q*^{#} from (35) on Montgomery streamlines on ^{#} from (B4) on Montgomery streamlines on the surface *b* = *b*_{0} (dashed red); the Euler–Ertel PV *σ*_{1} surface (dotted blue). The *σ*_{1} surfaces intersect (133.25°E, 54.04°S) at −700 m, while the mean height of the *b* surface is very near −700 m here. The streamlines cross 133.25°E at the indicated latitudes, as in Fig. 2.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

## APPENDIX C

### Thickness-Weighted Average Momentum, Thickness, and PV Budgets

Here, we review the necessary TWA theory to reach the PV budget (for full details, see Young 2012) interspersed with our numerical methods. Then, we illuminate and quantify the numerical problems of using TWA with archived *z*-level model diagnostics.

#### a. Theory and methods

##### 1) Coordinate-transformations and averages

*z*coordinates and buoyancy coordinates employs the following transformations of the partial derivatives:

##### 2) Interpolation

We evaluate *ξ*, *u*, *υ*, and *m* using PCHIPs [appendix B, section b(3)], but evaluate terms in the *z*-level buoyancy and momentum equations on a buoyancy surface using linear interpolation. Linear interpolants, unlike PCHIPs, preserve equations: terms that sum to zero on a grid also sum to zero after each term is linearly interpolated between grid points. (Cubic spline interpolation also preserves equations, but can overshoot the data—particularly problematic near the pycnocline.)

##### 3) Buoyancy-coordinate thickness equation

*w*stands alone before taking the

*z*derivative. This enables continuity, (33), to annihilate the second column of (C5), while the first column sums to zero by (29) and linearity of

Numerically, we implement *W* grid, then divide by *b _{z}* calculated by first difference. Take the

*z*derivative by first difference, returning to the tracer grid. Linearly interpolate this onto the buoyancy surface. Multiply by

##### 4) TWA thickness equation

This has no eddy terms; they are confined to the momentum equation, shown next. Numerically, we form (C7) simply by averaging (C6).

##### 5) Buoyancy coordinate momentum equation

*z*-coordinate momentum Eq. (24) transforms term by term into buoyancy coordinates as

Note (C8e) and its meridional version *b* surfaces possess an exact geostrophic streamfunction, *m* (section 3a). Numerically, the first column of the RHS of (C8) is obtained by interpolation of model outputs onto *ξ* averaged onto the *u* grid, while other RHS terms are calculated directly on the buoyancy surface with averages onto the *u* grid as necessary, and *hb _{x}* and

*hb*taken as

_{y}##### 6) TWA momentum equation

*h*, time average, then divide by

*h*times (C9); the second column is

*u*times (C6). Hence, dividing (C10) by

*E*are Reynolds stresses, while the last is the eddy form stress, which derives from and dominates the RHS of

^{u}*E*transport momentum within (across) buoyancy surfaces.

^{u}Numerically, the Reynolds terms in the second column of (C10) are taken as residuals, with all other terms directly computed. The vorticity and kinetic energy terms (C10c) and (C10d) are added, so *C* is not computed and

For (C10b) we require

##### 7) TWA potential vorticity equation

^{#}following the TWA velocity. Hence, Π

^{#}is materially conserved, except for the effects of mechanical forcing and dissipation, eddies, and nonconservation of buoyancy.

#### b. Numerical errors

On the RHS, we can add and subtract the stretching term *z*-level buoyancy and continuity equations are both delicately balanced and cannot be transformed into the TWA thickness term

*b*onto the

*u*grid, and

*u*onto the

*υ*grid. Any two variables

*α*and

*β*on the same grid satisfy

^{5}

Applying *b _{z}*, which does not isolate

*w*for the continuity equation to annihilate, as needed in (C5). Thus,

*u*,

*υ*, and

*h*on the buoyancy surface. This numerical imbalance is severe, because the buoyancy and continuity equations are both extremely finely balanced.

^{#}on the tracer grid:

On the RHS, vorticity is advected by *υ* and *u* grids, correctly producing terms in the discretized version of (C17). The LHS, however, contains vorticity averaged onto the tracer grid being advected by *u* and *υ* grids. This difference is captured by *ϵ*^{(Advec)}, taken as a residual. Provided the model resolution is such that the horizontal velocity varies smoothly laterally, this advection error is small relative to other terms (Fig. C1).

Maps of thickness-weighted average potential vorticity forcing terms (s^{−4}) on the buoyancy surface *b* = *b *(133.25°E, 54.04°S, −1000 m, 18–22 Sep 2008). From left to right and then from top to bottom, panels show tendency, rotational advection, divergent advection, diapycnal advection, diapycnal shear, momentum forcing, isopycnal Reynolds stresses, diapycnal Reynolds stresses, eddy form stress, stretching error, thickness error, and advection error. All maps are smoothed with a Gaussian filter of 5 × 5 grid points, and masked out (light gray) where the surface incropped or outcropped at any time during the record.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Maps of thickness-weighted average potential vorticity forcing terms (s^{−4}) on the buoyancy surface *b* = *b *(133.25°E, 54.04°S, −1000 m, 18–22 Sep 2008). From left to right and then from top to bottom, panels show tendency, rotational advection, divergent advection, diapycnal advection, diapycnal shear, momentum forcing, isopycnal Reynolds stresses, diapycnal Reynolds stresses, eddy form stress, stretching error, thickness error, and advection error. All maps are smoothed with a Gaussian filter of 5 × 5 grid points, and masked out (light gray) where the surface incropped or outcropped at any time during the record.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Maps of thickness-weighted average potential vorticity forcing terms (s^{−4}) on the buoyancy surface *b* = *b *(133.25°E, 54.04°S, −1000 m, 18–22 Sep 2008). From left to right and then from top to bottom, panels show tendency, rotational advection, divergent advection, diapycnal advection, diapycnal shear, momentum forcing, isopycnal Reynolds stresses, diapycnal Reynolds stresses, eddy form stress, stretching error, thickness error, and advection error. All maps are smoothed with a Gaussian filter of 5 × 5 grid points, and masked out (light gray) where the surface incropped or outcropped at any time during the record.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

*ϵ*terms on the RHS of (C24) are (rescaled) terms in the TWA vorticity equation (C17) and in the version of the TWA thickness equation (C7) that can be numerically balanced, i.e.,

#### c. Results: TWA PV budget in the Southern Ocean

Figure C1 maps terms in the TWA PV budget (C18), with the in-buoyancy-surface advective terms calculated from the LHS of (C24), but further separated into rotational and divergent components found by a Helmholtz decomposition of *ϵ*^{(Thick)} is among the largest terms in the PV budget, and the stretching error *ϵ*^{(Stretch)} is also sizeable (and substantially worse without smoothing—not shown).

To study the source of the thickness error, we repeat the PV budget with an alternate prescription of the thickness advection and stretching term, **v** · ∇*b* obtained from (30b) using model diagnostics for salt and heat fluxes. Now, we continue using **v** · ∇*b* using SOSE’s third-order upwind advection scheme with **v** the 5-day Eulerian-averaged velocity field and *S* and *θ*. This means the 5-day buoyancy budget (29) will not close, hence the buoyancy forcing is now taken as the residual, *ϖ* ≡ *b _{t}* +

**v**· ∇

*b*. Having marred the fine balance between large terms in the buoyancy equation, this new

*ϖ*is over an order of magnitude larger than the true

*ϖ*, and this will appear as fictitious diapycnal forcing in the PV budget. The goal here is not to fix the PV budget, but merely to locate the source of

*ϵ*

^{(Thick)}.

Figure C2 shows the PV budget under this alternative prescription. The thickness error is considerably reduced compared to Fig. C1, suggesting our numerical treatment of the coordinate transformation works decently (though still not perfectly). However, the errors are shunted to diapycnal forcings: predominantly the diapycnal shear, with smaller adjustments to the diapycnal advection and diapycnal Reynolds stress. This implies that a proper TWA analysis of PV cannot be done with SOSE’s Eulerian 5-day averaged data.

As in Fig. C1, but having instead computed the *z*-coordinate buoyancy advection offline, using the 5-day Eulerian-mean velocity and buoyancy computed from the 5-day Eulerian mean salinity and potential temperature with a centered second-order advection scheme. Only the subset of panels that differ from those in Fig. C1 are shown, namely, the diapycnal advection, diapycnal shear, diapycnal Reynolds stresses, and thickness error.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. C1, but having instead computed the *z*-coordinate buoyancy advection offline, using the 5-day Eulerian-mean velocity and buoyancy computed from the 5-day Eulerian mean salinity and potential temperature with a centered second-order advection scheme. Only the subset of panels that differ from those in Fig. C1 are shown, namely, the diapycnal advection, diapycnal shear, diapycnal Reynolds stresses, and thickness error.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

As in Fig. C1, but having instead computed the *z*-coordinate buoyancy advection offline, using the 5-day Eulerian-mean velocity and buoyancy computed from the 5-day Eulerian mean salinity and potential temperature with a centered second-order advection scheme. Only the subset of panels that differ from those in Fig. C1 are shown, namely, the diapycnal advection, diapycnal shear, diapycnal Reynolds stresses, and thickness error.

Citation: Journal of Physical Oceanography 52, 8; 10.1175/JPO-D-21-0195.1

Higher-frequency output is required, but most likely TWA should be done online. Indeed, Uchida et al. (2021) “found the difference between the eddy–eddy correlation terms diagnosed from instantaneous and temporally averaged fields to be up to same order of the total variance of the instantaneous field.” The buoyancy surface used in Figs. C1 and C2 can change by 200 m from one 5-day window to the next, in SOSE’s ACC. Even if TWA is done online, though, the not-insignificant stretching and advective errors would remain.

#### d. Discussion

Some recent studies have used the TWA formalism to study ocean dynamics using *z*-level models (Cessi and Wolfe 2013; Aoki et al. 2016), without simultaneously and consistently balancing the thickness and momentum budget. Instead, these studies construct the 3D TWA velocity, **v**^{#}, to be divergence-free: the horizontal TWA velocity is computed, its horizontal divergence calculated, and *w*^{#} obtained by vertical integration. This does not translate into buoyancy coordinates to give a diapycnal velocity *ϖ* as small as the model’s buoyancy forcing. Cessi and Wolfe (2013) recognize a “striking cancellation” between the mean and eddy-induced diapycnal velocity in their model eastern boundary current. While this may be true, their residual diapycnal velocity may well be an order of magnitude larger than the actual diapycnal velocity *ϖ*. Nonetheless, Cessi and Wolfe (2013) develop a TWA PV budget, but only by ignoring the diapycnal velocity (among other terms), which could be far larger in their budgets than they believe, and not negligible. TWA manipulations of *z*-level model data that do not accurately balance the thickness equation should be regarded with skepticism.

## REFERENCES

Aoki, K., A. Kubokawa, R. Furue, and H. Sasaki, 2016: Influence of eddy momentum fluxes on the mean flow of the Kuroshio extension in a 1/10° ocean general circulation model.

,*J. Phys. Oceanogr.***46**, 2769–2784, https://doi.org/10.1175/JPO-D-16-0021.1.Arakawa, A., and V. R. Lamb, 1977: Computational design of the basic dynamical processes of the UCLA general circulation model.

,*Methods Comput. Phys.***17**, 173–265, https://doi.org/10.1016/B978-0-12-460817-7.50009-4.Bodnariuk, N., C. G. Simionato, M. Saraceno, M. Osman, and L. B. Diaz, 2021: Interannual variability of the latitude of separation of the Brazil current: Teleconnections and oceanic Rossby waves propagation.

,*J. Geophys. Res. Oceans***126**, e2021JC017557, https://doi.org/10.1029/2021JC017557.Cessi, P., and C. L. Wolfe, 2013: Adiabatic eastern boundary currents.

,*J. Phys. Oceanogr.***43**, 1127–1149, https://doi.org/10.1175/JPO-D-12-0211.1.Combes, V., and R. P. Matano, 2014: Trends in the Brazil/Malvinas Confluence region.

,*Geophys. Res. Lett.***41**, 8971–8977, https://doi.org/10.1002/2014GL062523.Danabasoglu, G., J. C. McWilliams, and P. R. Gent, 1994: The role of mesoscale tracer transports in the global ocean circulation.

,*Science***264**, 1123–1126, https://doi.org/10.1126/science.264.5162.1123.Dellar, P. J., 2011: Variations on a beta-plane: Derivation of non-traditional beta-plane equations from Hamilton’s principle on a sphere.

,*J. Fluid Mech.***674**, 174–195, https://doi.org/10.1017/S0022112010006464.de Szoeke, R. A., 2000: Equations of motion using thermodynamic coordinates.

,*J. Phys. Oceanogr.***30**, 2814–2829, https://doi.org/10.1175/1520-0485(2001)031<2814:>2.0.CO;2.de Szoeke, R. A., S. R. Springer, and D. M. Oxilia, 2000: Orthobaric density: A thermodynamic variable for ocean circulation studies.

,*J. Phys. Oceanogr.***30**, 2830–2852, https://doi.org/10.1175/1520-0485(2001)031<2830:>2.0.CO;2.Dewar, W. K., and A. M. Hogg, 2010: Topographic inviscid dissipation of balanced flow.

,*Ocean Modell.***32**, 1–13, https://doi.org/10.1016/j.ocemod.2009.03.007.Gent, P. R., J. Willebrand, T. J. McDougall, and J. C. McWilliams, 1995: Parameterizing eddy-induced tracer transports in ocean circulation models.

,*J. Phys. Oceanogr.***25**, 463–474, https://doi.org/10.1175/1520-0485(1995)025<0463:PEITTI>2.0.CO;2.Gille, S. T., 1997: Why potential vorticity is not conserved along mean streamlines in a numerical Southern Ocean.

,*J. Phys. Oceanogr.***27**, 1286–1299, https://doi.org/10.1175/1520-0485(1997)027<1286:WPVINC>2.0.CO;2.Graham, R. M., and A. M. De Boer, 2013: The dynamical subtropical front.

,*J. Geophys. Res. Oceans***118**, 5676–5685, https://doi.org/10.1002/jgrc.20408.Hughes, C. W., 1995: Rossby waves in the Southern Ocean: A comparison of TOPEX/POSEIDON altimetry with model predictions.

,*J. Geophys. Res.***100**, 15 933–15 950, https://doi.org/10.1029/95JC01380.Hughes, C. W., 1997: Comments on “On the obscurantist physics of ‘form drag’ in theorizing about the Circumpolar Current.”

,*J. Phys. Oceanogr.***27**, 209–210, https://doi.org/10.1175/1520-0485(1997)027<0209:COOTOP>2.0.CO;2.Hughes, C. W., and B. A. De Cuevas, 2001: Why western boundary currents in realistic oceans are inviscid: A link between form stress and bottom pressure torques.

,*J. Phys. Oceanogr.***31**, 2871–2885, https://doi.org/10.1175/1520-0485(2001)031<2871:WWBCIR>2.0.CO;2.Hughes, C. W., M. S. Jones, and S. Carnochan, 1998: Use of transient features to identify eastward currents in the Southern Ocean.

,*J. Geophys. Res.***103**, 2929–2943, https://doi.org/10.1029/97JC02442.Johnson, G. C., and H. L. Bryden, 1989: On the size of the Antarctic Circumpolar Current.

,*Deep-Sea Res.***36A**, 39–53, https://doi.org/10.1016/0198-0149(89)90017-4.Jones, C. S., and P. Cessi, 2016: Interbasin transport of the meridional overturning circulation.

,*J. Phys. Oceanogr.***46**, 1157–1169, https://doi.org/10.1175/JPO-D-15-0197.1.Killworth, P. D., 1992: An equivalent-barotropic mode in the fine resolution Antarctic model.

,*J. Phys. Oceanogr.***22**, 1379–1387, https://doi.org/10.1175/1520-0485(1992)022<1379:AEBMIT>2.0.CO;2.Klocker, A., and D. P. Marshall, 2014: Advection of baroclinic eddies by depth mean flow.

,*Geophys. Res. Lett.***41**, 3517–3521, https://doi.org/10.1002/2014GL060001.Lu, Y., and D. Stammer, 2004: Vorticity balance in coarse-resolution global ocean simulations.

,*J. Phys. Oceanogr.***34**, 605–622, https://doi.org/10.1175/2504.1.Lumpkin, R., and S. Garzoli, 2011: Interannual to decadal changes in the western South Atlantic’s surface circulation.

*J. Geophys. Res.*,**116**, C01014, https://doi.org/10.1029/2010JC006285.MacCready, P., and P. B. Rhines, 2001: Meridional transport across a zonal channel: Topographic localization.

,*J. Phys. Oceanogr.***31**, 1427–1439, https://doi.org/10.1175/1520-0485(2001)031<1427:MTAAZC>2.0.CO;2.Marshall, D. P., 1995: Topographic steering of the Antarctic circumpolar current.

,*J. Phys. Oceanogr.***25**, 1636–1650, https://doi.org/10.1175/1520-0485(1995)025<1636:TSOTAC>2.0.CO;2.Marshall, D. P., 1997: Subduction of water masses in an eddying ocean.

,*J. Mar. Res.***55**, 201–222, https://doi.org/10.1357/0022240973224373.Marshall, D. P., D. R. Munday, L. C. Allison, R. J. Hay, and H. L. Johnson, 2016: Gill’s model of the Antarctic Circumpolar Current, revisited: The role of latitudinal variations in wind stress.

,*Ocean Modell.***97**, 37–51, https://doi.org/10.1016/j.ocemod.2015.11.010.Marshall, J., and T. Radko, 2003: Residual-mean solutions for the Antarctic Circumpolar Current and its associated overturning circulation.

,*J. Phys. Oceanogr.***33**, 2341–2354, https://doi.org/10.1175/1520-0485(2003)033<2341:RSFTAC>2.0.CO;2.Masich, J., T. K. Chereskin, and M. R. Mazloff, 2015: Topographic form stress in the Southern Ocean state estimate.

,*J. Geophys. Res. Oceans***120**, 7919–7933, https://doi.org/10.1002/2015JC011143.Matano, R. P., 1993: On the separation of the Brazil current from the coast.

,*J. Phys. Oceanogr.***23**, 79–90, https://doi.org/10.1175/1520-0485(1993)023<0079:OTSOTB>2.0.CO;2.Mazloff, M. R., P. Heimbach, and C. Wunsch, 2010: An eddy-permitting Southern Ocean state estimate.

,*J. Phys. Oceanogr.***40**, 880–899, https://doi.org/10.1175/2009JPO4236.1.McDougall, T. J., 1987: Neutral surfaces.

,*J. Phys. Oceanogr.***17**, 1950–1964, https://doi.org/10.1175/1520-0485(1987)017<1950:NS>2.0.CO;2.McDougall, T. J., 1989: Streamfunctions for the lateral velocity vector in a compressible ocean.

,*J. Mar. Res.***47**, 267–284, https://doi.org/10.1357/002224089785076271.McDougall, T. J., and D. R. Jackett, 1988: On the helical nature of neutral trajectories in the ocean.

,*Prog. Oceanogr.***20**, 153–183, https://doi.org/10.1016/0079-6611(88)90001-8.McDougall, T. J., and P. C. McIntosh, 1996: The temporal-residual-mean velocity. Part I: Derivation and the scalar conservation equations.

,*J. Phys. Oceanogr.***26**, 2653–2665, https://doi.org/10.1175/1520-0485(1996)026<2653:TTRMVP>2.0.CO;2.McDougall, T. J., and P. C. McIntosh, 2001: The temporal-residual-mean velocity. Part II: Isopycnal interpretation and the tracer and momentum equations.

,*J. Phys. Oceanogr.***31**, 1222–1246, https://doi.org/10.1175/1520-0485(2001)031<1222:TTRMVP>2.0.CO;2.Montgomery, R., 1937: A suggested method for representing gradient flow in isentropic surfaces.

,*Bull. Amer. Meteor. Soc.***18**, 210–212, https://doi.org/10.1175/1520-0477-18.6-7.210.Munk, W. H., and E. Palmén, 1951: Note on the dynamics of the Antarctic Circumpolar Current.

,*Tellus***3**, 53–55, https://doi.org/10.3402/tellusa.v3i1.8609.Olbers, D., 1998: Comments on “On the obscurantist physics of ‘form drag’ in theorizing about the Circumpolar Current.”

,*J. Phys. Oceanogr.***28**, 1647–1654, https://doi.org/10.1175/1520-0485(1998)028<1647:COOTOP>2.0.CO;2.Ridgway, K. R., and J. R. Dunn, 2007: Observational evidence for a Southern Hemisphere oceanic supergyre.

,*Geophys. Res. Lett.***34**, L13612, https://doi.org/10.1029/2007GL030392.Smith, K. S., 2007: The geography of linear baroclinic instability in Earth’s oceans.

,*J. Mar. Res.***65**, 655–683, https://doi.org/10.1357/002224007783649484.Stanley, G. J., 2019a: Neutral surface topology.

,*Ocean Modell.***138**, 88–106, https://doi.org/10.1016/j.ocemod.2019.01.008.Stanley, G. J., 2019b: The exact geostrophic streamfunction for neutral surfaces.

,*Ocean Modell.***138**, 107–121, https://doi.org/10.1016/j.ocemod.2019.04.002.Stanley, G. J., T. E. Dowling, M. E. Bradley, and D. P. Marshall, 2020: Ertel potential vorticity versus Bernoulli potential on approximately neutral surfaces in the Antarctic Circumpolar Current.

,*J. Phys. Oceanogr.***50**, 2621–2648, https://doi.org/10.1175/JPO-D-19-0140.1.Stewart, A. L., J. C. McWilliams, and A. Solodoch, 2021: On the role of bottom pressure torques in wind-driven gyres.

,*J. Phys. Oceanogr.***51**, 1441–1464, https://doi.org/10.1175/JPO-D-20-0147.1.Stommel, H., 1957: A survey of ocean current theory.

,*Deep-Sea Res.***4**, 149–184, https://doi.org/10.1016/0146-6313(56)90048-X.Sverdrup, H. U., 1947: Wind-driven currents in a baroclinic ocean; with application to the equatorial currents of the eastern Pacific.

,*Proc. Natl. Acad. Sci. USA***33**, 318–326, https://doi.org/10.1073/pnas.33.11.318.Tamsitt, V., and Coauthors, 2017: Spiraling pathways of global deep waters to the surface of the Southern Ocean.

,*Nat. Commun.***8**, 172, https://doi.org/10.1038/s41467-017-00197-0.Tansley, C. E., and D. P. Marshall, 2001: On the dynamics of wind-driven circumpolar currents.

,*J. Phys. Oceanogr.***31**, 3258–3273, https://doi.org/10.1175/1520-0485(2001)031<3258:OTDOWD>2.0.CO;2.Thompson, A. F., A. L. Stewart, and T. Bischoff, 2016: A multibasin residual-mean model for the global overturning circulation.

,*J. Phys. Oceanogr.***46**, 2583–2604, https://doi.org/10.1175/JPO-D-15-0204.1.Uchida, T., Q. Jamet, W. K. Dewar, D. Balwada, J. Lesommer, and T. Penduff, 2021: Diagnosing the thickness-weighted averaged eddy-mean flow interaction in an eddying North Atlantic ensemble, version 2. Earth and Space Science Open Archive, https://doi.org/10.1002/essoar.10504082.2.

Warren, B. A., J. H. LaCasce, and P. E. Robbins, 1996: On the obscurantist physics of “form drag” in theorizing about the Circumpolar Current.

,*J. Phys. Oceanogr.***26**, 2297–2301, https://doi.org/10.1175/1520-0485(1996)026<2297:OTOPOD>2.0.CO;2.Young, W. R., 2010: Dynamic enthalpy, conservative temperature, and the seawater Boussinesq approximation.

,*J. Phys. Oceanogr.***40**, 394–400, https://doi.org/10.1175/2009JPO4294.1.Young, W. R., 2012: An exact thickness-weighted average formulation of the Boussinesq equations.

,*J. Phys. Oceanogr.***42**, 692–707, https://doi.org/10.1175/JPO-D-11-0102.1.Zajaczkovski, U., 2017: A study of the Southern Ocean: Mean state, eddy genesis & demise, and energy pathways. Ph.D. thesis, University of California, San Diego, 125 pp.

Zhai, X., H. L. Johnson, and D. P. Marshall, 2010: Significant sink of ocean-eddy energy near western boundaries.

,*Nat. Geosci.***3**, 608–612, https://doi.org/10.1038/ngeo943.Zhang, H.-M., and N. G. Hogg, 1992: Circulation and water mass balance in the Brazil Basin.

,*J. Mar. Res.***50**, 385–420, https://doi.org/10.1357/002224092784797629.