1. Introduction
Overturning cells of various shapes and sizes are a pervasive feature of the global ocean, ranging from the small eddy-induced circulations to the global-scale meridional overturning circulation (MOC). Overturning circulations are often analyzed in simplified frameworks like the Stommel box model, which consists of two connected, well-mixed reservoirs of different temperature and salinity and features different circulation modes depending on the relative strength of the thermal and haline forcings (Stommel 1961).
The thermohaline loop is a variation of the Stommel box model in which properties can vary continuously along a one-dimensional loop, allowing for a more realistic representation of the advective-diffusive balance (e.g., Welander 1967; Huang 1999; Wunsch 2005). It does not only find application in physical oceanography but in the form of the so-called thermosyphons also in nuclear and solar energy engineering (e.g., Zvirin 1982; Miljkovic and Wang 2011). Despite their simple setup, thermohaline loops exhibit a variety of dynamical behaviors including instability and chaos (Malkus 1972; Welander 1986; Yuan and Wunsch 2005).
Although it is unclear to what extent analogies between these conceptual representations and the real ocean can be drawn, it has been argued that by capturing the basic physics these idealized models can indeed provide an intuitive understanding of more complex, higher-dimensional processes (Bryan 1986; Ruddick and Zhang 1996; Manabe and Stouffer 1999; Held 2005).
Wunsch (2005) devised a circular loop model with point sources and sinks of salinity and temperature applied on either side of the loop, typically at an angle of 45° from the loop’s top. The case of applying sources and sinks at the same height is of particular interest for ocean studies, as it represents the case of horizontal convection best (Hughes and Griffiths 2008). This means that in the Wunsch model part of the circulation occurs at a level above the position of sources and sinks. In the real ocean, however, the thermohaline forcing (leaving aside geothermal heating) is applied directly at the ocean surface, introducing an asymmetry between an essentially horizontal surface flow and a deeper, both horizontal and vertical, one in the interior. We therefore propose a simple modification of the Wunsch model where the upper portion of the loop above the level of source and sink is folded down until it becomes purely horizontal (cf. Fig. 1).1

(a) Schematic description of the fluid loop as proposed by Wunsch (2005) and (b) the modified folded version. The loop coordinate is represented as l = aΦ, where the angle Φ is measured clockwise from the top and a denotes the radius. Sources and sinks of temperature and salinity are applied at l+ and l−, respectively, which are located at the same height Zf. The azimuthal velocity w is constant around the loop for reasons of continuity. In the modified version, the loop is folded at the level of source and sink until the part above is purely horizontal. The steady-state temperature distribution is also shown; blue and red colors denote cold and warm fluid, respectively, with dashed lines representing the strong temperature gradients within the diffusive boundary layers.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
In this study, we investigate the dynamical modifications introduced when folding the loop, considering a simplified nonlinear EOS similar to the one in Vallis (2006, p. 34). In most idealized studies, including the one by Wunsch (2005), a linear EOS is assumed. Yet, in the real ocean, the EOS is nonlinear mainly through a strong dependence of the seawater thermal expansion on temperature and pressure, giving rise to cabbeling and thermobaricity, respectively. Both effects induce vertical advection through neutral surfaces (albeit only downward in the case of cabbeling) and have been shown to considerably affect the circulation and water mass properties (e.g., McDougall 1987; Klocker and McDougall 2010; Schanze and Schmitt 2013; Hieronymus and Nycander 2013). One key aspect of this study is to explore how the impact of these two nonlinear effects changes when the loop is folded.
This paper is organized as follows: In section 2, we describe the loop model, discuss it analytically, and present the numerical implementation. We then describe the dynamics of the linear EOS scenario (section 3) before addressing the nonlinear case in section 4. An extensive discussion of the results and their applicability is given in section 5, and a brief summary is given in section 6.
2. The loop model
a. Governing equations
The model consists of a closed fluid loop with an infinitesimal section area, so that radial velocities can be neglected. The loop natural coordinate l is defined as the distance along the circle from the top of the loop, increasing clockwise. The total length of the loop is given by L = 2πa, where a denotes the radius. The focus of this study is on heat-induced horizontal convection, where heating and cooling occur at the same height and no salinity forcing is applied. Two geometries of the loop will be considered (cf. Fig. 1):
- The circular loop: Essentially the same model as in Wunsch (2005). The loop coordinate is l = aΦ, where Φ is simply the angle from the top, which is defined positive for clockwise turning. Height is thus given as z(l) = a cos(l/a). Sinks and sources of temperature and salinity are applied at positions l− = aΦ− ≤ L/2 and l+ = aΦ+ = L − l−, respectively, that is, at the same height Zf.
- The folded loop: The part above source and sink (z ≥ Zf in the circular loop) is now folded down horizontally. The loop coordinate is kept unchanged, only the height function is modified, with z(l) = min[Zf, a cos(l/a)]. By analogy, an equivalent angle Φ = l/a can still be defined, although it must be kept in mind that it is now associated with two different circular paths, one vertical for the interior flow and one horizontal for the surface flow.


- Π(l) = sin(l/a) for the circular loop, and
- Π(l) = m(l) sin(l/a) for the folded loop,












b. Nondimensionalization








Apart from geometrical parameters (i.e., circular vs folded cases and the vertical position of forcings Zf), several nondimensional parameters are defined in the loop model:
- the inverse Rayleigh number
, which corresponds to the ratio of diffusive to advective time scales and controls the size of the diffusive boundary layers near the source and sink (as will be shown later); - the momentum stress forcing τ′ as a representation of frictional forcing;
- the Grashof number
, which represents the ratio of frictional to advective time scales. Assuming inertialess flow, that is, F → 0, is equivalent to neglecting frictional oscillations. Equation (15) then simplifies to - the salt forcing strength parameter
, which scales the magnitude of the salinity forcing relative to that of the temperature forcing; when it is set to one, the temperature and salinity forcings cancel exactly when a linear EOS is assumed. Note that a circulation can be generated even if η = 1 if the EOS is nonlinear. No salt forcing will be considered in this paper, that is, η will be kept to zero; and - finally, the cabbeling and thermobaricity parameters λ′ and μ′.
The only formal difference between the circular and folded loop model concerns the buoyancy term, which is masked in the folded loop case. Another, somewhat minor difference is that the definition of height is changed with consequences for the thermobaric term in the EOS.
In the rest of the paper, we will refer to nondimensionalized quantities only, and the primes will be dropped for clarity.
c. The numerical model
Following the approach by Wunsch [2005, cf. his Eqs. (11)–(14)], the system of equations—Eqs. (16)–(19), using the inertialess assumption—can be solved analytically for the equilibrium velocity w if a linear EOS is assumed. This procedure can be extended to nonlinear scenarios including only thermobaricity; when adding the cabbeling term or folding the loop, however, the ability to compute a simple analytical solution is lost. Therefore, a numerical investigation of loop dynamics is required.
The loop model is implemented in Fortran90, using a forward differencing scheme for the diffusive part and a leapfrog scheme for the advection and forcing terms. With a time increment of Δt = 10−4 and 360 grid points around the loop (i.e., ΔΦ = 1°), the diffusion scheme is stable for R < 0.4. Since the velocity overshoot is typically below w = 2, the Courant–Friedrich–Levy criterion w ≤ ΔΦ/Δt ≈ 175 is easily fulfilled.
Our standard parameter settings are R = 0.1 (weak diffusion), τ = 0 (no wind stress forcing), and η = 0 (no salinity forcing). The simulations are started at rest (w = 0) with temperatures set to zero everywhere. The thermal forcing, applied at the same height (i.e., Φ+ + Φ− = 2π), is switched on at t = 0. Our main focus is on the scenario Zf = 0.5 (i.e., Φ− = π/3; cf. Fig. 1). With these settings, equilibrium is typically reached at t ≈ 15 in the circular loop model and at t ≈ 20 in the folded loop model. Unless otherwise noted, the simulations involve 106 iterations (i.e., tend = 100), which take less than a minute to perform on a standard computer.
The numerical model is validated based on scenarios that allow for the computation of an analytical solution for w. We find that for our standard parameter settings, the difference between numerical and analytical solutions is several orders of magnitude smaller than the respective analytical velocity. The remaining inaccuracies are therefore considered too small to influence the conclusions drawn here.
3. Loop dynamics with a linear EOS





a. The transitory regime
The simulations are started from a state of rest with temperatures set to zero everywhere. When the thermal forcing is turned on at t = 0, temperature anomalies are generated at the source and sink (cf. Fig. 2) that induce horizontal density anomalies and thus a nonzero velocity. This will in turn advect the temperature anomalies away from the forcing region and further modify velocity. Because of the influence of the curvature term, the highest velocities can be obtained when source and sink are applied at Zf = 0, where sin(Φ) = 1. The circulation reaches its maximum speed when the horizontal density anomaly, integrated around the loop, is largest in accordance with Eq. (19). For example, when source and sink are applied at Zf = 0 in the circular loop model, this is given when the warm front reaches the top and the cold front reaches the bottom of the loop. As velocity decreases, the fluid is advected more slowly past the source and sink, which increases their efficiency and leads to higher temperature amplitudes (and vice versa). There are further oscillations in temperature (cf. Figs. 2 and 3) that are associated with the arrival of warm (cold) temperature anomalies at the sink (source) and whose amplitudes are steadily decreased by diffusion. Yet, those temperature oscillations do not yield velocity oscillations because they tend to compensate almost exactly upon integration around the loop.

Initial evolution of temperature as a function of time and position for the circular and the folded loop model. (a),(c) When source and sink are applied at Zf = 0, the temperature amplitudes are the same at source and sink. (b),(d) Shifting the source and sink up or down, for example, with Zf = 0.5, leads to an asymmetric temperature distribution. (e),(f) When the loop is folded, the temperature distribution is similar but exhibits higher amplitudes because of the significant velocity decrease, which is underlined by the temperature difference ΔΘ = Θ(full) − Θ(folded). Contour lines are placed at intervals of 1 in (a)–(d) or 0.5 in (e) and (f).
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1

Variation of velocity and temperature at the source and sink with time in the linear EOS scenario for the (a) circular loop and the (b) folded loop model. The highest equilibrium velocities are obtained when source and sink are applied at Zf = 0, where the curvature term Π is maximum. The corresponding temperature variation demonstrates that faster advection results in lower temperature amplitudes. Significant changes in the variability of the transitory regime are observed when the loop is folded down.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
When the loop is folded, significant differences in magnitude are observed. The velocity is generally reduced as expected, and it takes longer for the system to reach a steady state (cf. Fig. 3). Also, the temporal variations of velocity and temperature are modified substantially. For example, in the circular loop the maximum velocity overshoot is obtained for the case Zf = 0 (Fig. 3a, red curve), while in the folded loop the velocity maximum is highest for the setting Zf = 0.5 (Fig. 3b, green curve). Moreover, in the latter case oscillations in velocity are observed that do not exist in the circular loop case. This is because once the loop is folded, buoyancy anomalies produced in the upper part of the loop can no longer compensate anomalies in its lower part (see Fig. 3).
b. Steady-state properties
The steady state is reached when a diffusive-advective temperature balance is achieved everywhere around the loop. Then, temperatures are negative in most of the lower loop between the sink and source and positive in most of the upper loop between source and sink (cf. Figs. 1 and 2). Because of diffusion, however, the temperature changes sign slightly above the sink and slightly below the source, which gives rise to a small region around Zf where Θ and thus σ change sign horizontally. This region is pivotal for maintaining the circulation in equilibrium; if the temperature and hence the density distribution were exactly symmetrical around the vertical axis, the horizontal density anomaly would be zero everywhere and there could be no circulation.


Equation (24) indicates that the relative length of the warm path determines its temperature amplitude, that is, the temperature amplitude at the source, relative to that at the sink. In the symmetric case with Zf = 0, the warm path is as long as the cold path (L+ = L−), and consequently, ΔT+ = −ΔT−. Positioning source and sink at Zf = 0.5, on the other hand, means that the warm branch is half as long as the cold branch, which requires ΔT+ = −2ΔT− in equilibrium.


Steady-state velocity and temperature values for different positions of source and sink in the circular and the folded loop model.

The diffusive boundary layer of the loop’s right branch is situated above the sink, which in the folded loop model lies within the folded part that does not contribute to the buoyancy torque. In the dynamically relevant part of the folded loop, the section where horizontal temperature anomalies are high is then approximately half as large as in the circular model (cf. Fig. 1). This partly explains the observed reduction in velocity compared to the circular model (Δw = −18.3% for Zf = 0.5 and Δw = −22.4% for Zf = 0). However, the velocity decrease is somewhat smaller than expected considering the length of loop that is folded (30% and 50%, respectively). This is because the lower velocity induces, in turn, a distinct rise in temperature amplitudes as predicted by Eq. (26), leading to a compensation through slightly increased horizontal density gradients.
c. Sensitivity to the forcing position Zf
In the circular loop, the effect of vertically shifting the level of source and sink is symmetric about z = 0, that is, at any time, the same loop velocity is obtained for Zf and −Zf (cf. Table 1). Choosing Zf = −0.5 instead of Zf = 0.5, for example, means that the cold path is half instead of twice as long as the warm path. The relative length of upper and lower paths, however, is not affected; in consequence, the temperature distribution is the same but of opposite sign, and the velocity remains unaffected. This symmetry is broken in the folded loop model; since the loop is always folded at the level of source and sink, the size of the dynamically important part of the loop is varied when source and sink are shifted vertically.
An asymmetric positioning of source and sink (i.e., Zf ≠ 0) has important consequences for dynamics: First, the timing of the various dynamical events becomes asynchronous, as the warm front reaches the top or the sink earlier/later than the cold front reaches the bottom or the source (cf. Fig. 2b). Second, because of their different amplitudes, the warm or cold temperature anomalies no longer compensate (the asymmetric temperature distribution begins to be established at t ≈ 0.6 for Zf = 0.5). This causes oscillations and even sign changes in the horizontal density anomaly (cf. Fig. 4), which approaches its equilibrium distribution smoothly in both model setups when Zf = 0 (not shown).

Time evolution of the horizontal density anomaly σa as a function of height in the (a) full circular loop and the (b) folded loop model for Zf = 0.5 and Φ ∈ (0, π). Contour lines are placed at intervals of 1. Note that the vertical integral of σa is equal to half the loop velocity.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
4. Cabbeling and thermobaricity
Introducing cabbeling only, the EOS assumes the form σ = σlin + σcabb. Unlike the linear term σlin = −Θ, the nonlinear contribution σcabb = −0.5λΘ2 has the same sign everywhere, leading to a buoyancy gain throughout the dynamically relevant part of the loop. Including only thermobaricity in the EOS leads to the form σ = σlin + σtherm, where the nonlinear contribution σtherm = μzΘ changes sign both with height and temperature. Consequently, cabbeling and thermobaricity are expected to induce different dynamical modifications compared to the linear scenario.
a. Effect on the steady state
As illustrated in Fig. 4, the region around source and sink is pivotal for loop dynamics; the impact of the nonlinear effects at that level will therefore determine how they affect the steady-state velocity of the circular loop. We will first discuss the case of introducing cabbeling only.
At the source, density is negative and reduced even further due to cabbeling (assuming λ > 0). This leads to an increased tendency to rise and hence a faster circulation (represented by the longer red arrow in the schematic of the circular loop model in Fig. 5a). At the sink, on the other hand, density is positive, so that the nonlinear buoyancy gain due to cabbeling diminishes the fluid’s tendency to sink (represented by the shorter blue arrow in Fig. 5a). This causes a deceleration of the flow. With a completely symmetric geometry, where source and sink are applied at z = 0, the buoyancy increase by cabbeling is equally strong in the warm and the cold branch. In consequence, the effects at source and sink compensate exactly, and cabbeling does not influence the steady-state velocity in the circular loop model, which is indeed observed and can even be shown analytically.

Schematic description of the effect of (a) cabbeling and (b) thermobaricity in both model versions with Zf = 0.5. Full arrows represent the linear scenario, while dashed arrows denote the nonlinear contribution. Blue and red colors describe cold and warm fluid, respectively. The effect on steady-state velocity is illustrated in black.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
For Zf ≠ 0, the temperature distribution determines which of the opposing effects of cabbeling at source and sink ultimately dominates; with an asymmetric position of source and sink, temperature amplitudes in the warm branch differ from those in the cold branch. Because cabbeling is a quadratic function of temperature, it has the highest effect at the source and within the warm branch for the standard setting of Zf = 0.5, for which it induces a velocity increase (cf. the schematic illustration in Fig. 5a). This is underlined by the variation of steady-state velocity with λ (cf. Fig. 6a): a positive cabbeling parameter leads to an acceleration of the flow, while λ < 0 decelerates the circulation. In the former case, cabbeling induces a decrease in mass

Velocity and mass as a function of μ and λ. The linear reference value is shown in red. When the amplitude of the cabbeling parameter λ is too large, an oscillatory instability is induced, which is why the full range between −1 ≤ λ ≤ 1 is not depicted for the circular loop model. In the folded loop model, higher values of λ lead to ever lower velocities; the different behavior for λ ≥ 0.3 is related to the shift from an advection- to a diffusion-dominated regime. In that regime, the circulation is very slow, which results in very high, negative values for mass.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
This compensation between these opposite effects is also illustrated in the difference in the horizontal density anomaly with respect to the linear scenario [cf. Fig. 7a for the circular loop model with λ = 0.1 and Φ ∈ (0, π)]: in the upper loop, where temperatures are positive, cabbeling acts to increase the horizontal density anomaly of the right branch, while in the lower right branch, Δσa assumes negative values, indicating a reduced density anomaly due to cabbeling. The higher amplitudes in the warm section then lead to a stronger effect of cabbeling in the upper than in the lower loop, resulting in the observed velocity increase.

Difference in the horizontal density anomaly σa between the case with λ = 0.1 and the one based on a linear EOS (λ = 0, depicted in Fig. 4) for the (a) circular loop and the (b) folded loop with Zf = 0.5. In the circular loop, the largest differences are observed in the upper part, where temperature amplitudes are highest, especially at the level Zf = 0.5. When the loop is folded, the largest differences are found at and slightly below Zf = 0.5, that is, where the diffusive boundary layer of the left branch is located. Contour lines are placed at intervals of 0.4, and again only Φ ∈ (0, π) is considered in the computation of σa.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
When the loop is folded, cabbeling (and thermobaricity) can influence the steady-state velocity even when Zf = 0, since a perfect cancellation between warm and cold loop parts is no longer possible. For better comparison, however, source and sink are applied at z = 0.5 like in the circular loop. Moreover, we again focus on the scenario λ = 0.1, which allows for linear effects to dominate loop dynamics, while the influence of cabbeling is expected to be sufficiently strong to be discernible. In that case, cabbeling induces a 5% deceleration of the circulation in contrast to the 4.6% increase in velocity in the circular loop model (cf. Table 2 and Fig. 5a for a schematic description of how cabbeling affects dynamics in the folded loop model). The variation of velocity with λ (cf. Fig. 6b) underlines that the influence of cabbeling is reversed compared to the circular loop model. This is because the upper, warm branch is no longer relevant for velocity. The horizontal density anomaly difference to the linear EOS case (cf. Fig. 7b) illustrates that in most of the dynamically important part of the loop a decrease in σa is induced when cabbeling is included in the EOS. There, temperatures are mostly negative, and the nonlinear buoyancy gain acts to retard the flow, which causes the observed velocity decrease. However, the downward diffusion of heat from the source leads to a small section where temperatures are positive and cabbeling acts to accelerate the flow (Δσa > 0); hence even in the folded loop a compensation between the retarding and the accelerating effect of cabbeling is observed.
Final results for velocity, mass, temperature, and density at the source in the folded and the circular loop model for various forms of the EOS.

The velocity decrease due to cabbeling results in higher temperature amplitudes (cf. Table 2), which is also why mass is affected more strongly than in the circular loop (cf. Fig. 6b). In consequence, the amplitude of the horizontal density anomaly difference between the nonlinear and the linear scenario is also higher than in the circular loop (cf. Fig. 7).
In both model setups, an oscillatory instability can be provoked when the cabbeling parameter is sufficiently large or in the circular loop for any choice of λ when the linear density forcings by salinity and temperature cancel (η = 1). Critical values are λ < −0.385 and λ > 0.951 for η = 0 in the circular loop model, which might vary slightly when the simulations are run arbitrarily long (here, tend = 500). Nevertheless, they suggest that stability is affected differently when the nonlinear and linear contributions to density do not act in the same way (λ < 0) and when they do (λ > 0). In the folded loop, different parameter settings induce this instability. A detailed analysis of this phenomenon will be deferred to future work.
Considering only thermobaricity in the EOS, the nonlinear contribution σtherm = μzΘ varies both with temperature and height. Throughout most of the loop, Θ and z have the same sign, so that the thermobaric term is positive for μ < 0. In that case, thermobaricity induces an increase in mass, which is illustrated in Fig. 6a. At the source, where Θ > 0, the thermobaric term is positive if μ > 0 in the standard scenario with Zf = 0.5, which leads to a denser fluid and thus a reduced tendency to rise. At the sink, where Θ < 0, the thermobaric term is negative and acts to increase buoyancy. Both processes lead to a velocity decrease, which is shown schematically in Fig. 5b. When μ < 0, the opposite effect is obtained, as the variations of mass and velocity with μ (cf. Fig. 6a) show.
When source and sink are not applied at the loop’s equator, the thermobaric term changes sign twice in the circular loop model. With μ > 0 and Zf = 0.5, it causes a deceleration in the upper part of the loop (z > 0) and an acceleration for z < 0. The overall impact of thermobaricity on steady-state velocity is governed by its effect at source and sink and in that part of the loop with the highest temperature amplitudes; for μ = 0.1 and the standard parameter settings, velocity is decreased by 1.6% in the circular loop model. When the loop is folded, there still is a compensation between the counteracting effects of thermobaricity in the upper and lower part of the loop, but the total impact on steady-state velocity is reduced: the deceleration observed in the circular loop for z > Zf = 0.5 is irrelevant for dynamics in the folded loop model. In consequence, the accelerating effect of thermobaricity below z = 0 becomes more important, leading to a velocity change of Δw = −0.9% compared to the linear case (illustrated in Fig. 5b). Contrary to the cabbeling scenario, the region around Zf still determines the effect of thermobaricity on equilibrium velocity in the folded loop—the section where Δσa < 0 extends from z = 0.5 to z = 0.2 (cf. Fig. 8b), while for cabbeling Δσa changes sign already at z = 0.4 (cf. Fig. 7b).

As in Fig. 7, but varying the thermobaricity parameter μ instead of the cabbeling parameter λ. Differences are shown between the case where μ = 0.1 and the linear case μ = 0 (depicted in Fig. 4) and contour lines are placed at intervals of 0.05.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
When source and sink are applied below z = 0, the effect of thermobaricity on velocity is reversed in both model versions (cf. the positive horizontal density difference below z = 0 shown in Fig. 8): the dynamically important region, where temperatures change sign horizontally (i.e., within the boundary layers around the level of source and sink), as well as that part of the loop where temperature amplitudes are highest are then associated with a negative value of z.
The final results for velocity, temperature, and mass enlisted in Table 2 show that in both model geometries the effects of cabbeling and thermobaricity are superimposed almost linearly in a combined scenario with λ = μ = 0.1. For example, velocity is increased by 4.6% due to cabbeling in the circular loop, decreased by 1.6% due to thermobaricity, and increased by 3.2% in a combined scenario. In the folded loop, on the other hand, both effects add up with respect to their influence in velocity, leading to a deceleration by 6.0% for λ = μ = 0.1. Mass on the other hand is still altered oppositely by the two nonlinear processes, which have a stronger effect than in the circular loop model because of the higher temperature amplitudes caused by the folding and the resultant velocity decrease.
In the parameter ranges discussed here, cabbeling has a somewhat stronger effect on dynamics and steady-state properties in both model versions. This is because the thermobaric term is only a linear, not a quadratic, function of temperature and because z is maximum where the curvature term Π is minimum and vice versa. Yet, for other parameter choices the influence of thermobaricity and cabbeling could be comparable, or their relative importance might even be reversed. Note, too, that with respect to mass both nonlinear effects have a significant influence—when a linear EOS is used, mass is zero in equilibrium for reasons of heat conservation.
b. Effects on the mass distribution





Variation of velocity, temperature, and mass with time in the (a) circular loop and the (b) folded loop model for different versions of the EOS with Zf = 0.5. For clarity, we show the scenario μ = 0.5; the behavior is qualitatively the same as in the case with μ = 0.1, which can hardly be distinguished from the scenario with a linear EOS. The amplitudes of the steady-state velocity (temperature) are decreased (increased) by less than 10%, and the final mass is increased by a factor of 5 when considering μ = 0.5 instead of μ = 0.1.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1



c. The transitory regime
To assess the influence of the nonlinear effects during the transient regime, the difference in the horizontal density anomaly with respect to the linear reference case is analyzed. As in the previous subsection, we first analyze the scenario with λ = 0.1 (cf. Fig. 7 for the difference in σa using Zf = 0.5). In both model versions, oscillations between positive and negative Δσa can be observed for t < 10, indicating that cabbeling intensifies or weakens certain dynamical events. For example, when the cold front propagates down the right branch right after the onset of the forcing, the buoyancy gain due to cabbeling leads to a lighter right side and hence a reduced horizontal density difference compared to the linear case. When the warm front arrives at the sink and travels down the right branch for t > 1.1, the right side becomes lighter; this is intensified by cabbeling so that the horizontal density anomaly is higher than when a linear EOS is considered. These modifications of the horizontal density anomaly are the main reason for the different variation of velocity with time shown in Fig. 9. The magnitude of the changes brought about by cabbeling differs in the full and folded loops because the latter exhibits higher temperature amplitudes as a consequence of the folding.
This effect of cabbeling on buoyancy has interesting consequences in the folded loop model, where velocity oscillations can be observed in the linear case when an asymmetric positioning of source and sink is chosen (cf. Fig. 9 for the scenario with Zf = 0.5). These oscillations, however, are dampened when cabbeling is accounted for in the EOS: the velocity increase for t ≥ 3.8 is brought about by the arrival of the warm minimum at the sink, causing anomalously cold and thus dense fluid, which increases the horizontal density anomaly and hence accelerates the flow. In the presence of cabbeling, this densification is counteracted by the nonlinear buoyancy gain. In consequence, σa is not increased as much as when a linear EOS is considered, and the acceleration is weaker. This is confirmed by the negative density anomaly difference observed in the lower loop for t ≥ 3.8 (cf. Fig. 7).
Considering now the variation of the horizontal density anomaly difference between the nonlinear scenario with μ = 0.1 and the linear one (cf. Fig. 8), we see that, like cabbeling, the thermobaric effect intensifies and weakens the various dynamical events. Especially during the first time steps, the additional dependence on height is imminent; once the cold front propagating down the right branch for t > 0 passes below z = 0, the thermobaric term becomes positive, resulting in a densification of the fluid and therefore an increase in Δσa and velocity compared to the linear case.
d. Effect of salinity and wind stress forcing
We now briefly discuss the influence of introducing salinity or wind stress forcing. Figure 10 shows the variation of velocity, temperature, and mass for η = 0.5, that is, with a linear temperature forcing twice as strong as the linear salinity forcing, which effectively halves the net linear density forcing compared to the standard case (η = 0). Significant changes in the transient behavior can be observed compared to the purely thermally driven circulation (cf. Fig. 9). There also is a noticeable effect on the equilibrium properties; for the settings λ = μ = 0.1, the final velocity is increased by 7.6% relative to the linear scenario in the full loop (compared to 3.2% for η = 0) and decreased by 16.8% in the folded loop (compared to 6.0% for η = 0). This stronger influence of cabbeling and thermobaricity is expected because when the linear contributions of salinity and temperature to the buoyancy torque are opposed, the relative importance of nonlinear effects is naturally increased.

As in Fig. 9, but for η = 0.5; that is, the linear forcing by salinity is half as strong as that by temperature. Since temperature and salinity influence density in a counteracting way, the relative importance of nonlinear effects is enhanced, which is clearly visible in the larger differences between linear and nonlinear scenarios compared to the purely thermally driven case shown in Fig. 9, especially during the transient response.
Citation: Journal of Physical Oceanography 45, 10; 10.1175/JPO-D-15-0022.1
The influence of the wind stress forcing τ on the system’s behavior and stability properties is discussed in detail by Wunsch (2005) and Yuan and Wunsch (2005) for a linear equation of state. Taking λ = μ = 0.1 and setting for example τ = η = 0.5 induces a circulation that is 60.8% faster than in the linear scenario without salinity or wind stress forcing in the full loop and 65.7% faster in the folded loop. A detailed investigation of the impact of wind and salinity forcing or their interaction, however, is beyond the scope of the present paper and will be deferred to future work.
5. Discussion
Either folded or not, the loop model is without doubt a simplified description of large-scale oceanic overturning cells such as the three-dimensional, temporally varying meridional overturning circulation. It features the simplest dynamics one could think of, consisting in a single time-dependent scalar for velocity and consequently no net momentum advection. Yet, it can potentially embed the full richness and complexity of ocean thermodynamics, through the ability to use any equation of state, and also allows the study of more exotic fluid thermodynamics. Apart from being interesting in its own right, the loop model thus has great educational value when it comes to the analysis of overturning flow behavior or the demonstration of the Sandström theorem (Wunsch 2005).
A natural classification of overturning models can be obtained based on the number of dimensions used to describe the tracer fields. The Stommel box model, involving only a scalar to represent a tracer field in each box, is essentially a zero-dimensional model, while the thermohaline loop is the one-dimensional equivalent. Two-dimensional overturning models are often used in the study of horizontal convection, either in numerical simulations or laboratory experiments (Hughes and Griffiths 2008), and, finally, three-dimensional models simply correspond to the real ocean simulated with general circulation models. In this perspective, it comes as evidence that the thermohaline loop features a wider (and somewhat more realistic) range of dynamical behaviors than the Stommel box model. In particular, it allows for a more transparent investigation of loop dynamics insofar as the loop does not feature discontinuous, infinitely large, or zero Rayleigh numbers (Dewar and Huang 1996; Wunsch 2005). Importantly, this also implies that the Stommel box model is a limit case of the thermohaline loop (in practice, when the inverse Rayleigh number tends toward zero), indicating that steady states observed in a box model can be reproduced with a thermohaline loop but not necessarily vice versa.
The possibility of three-dimensional analogs of the thermohaline loop is discussed by Wunsch (2005, cf. his section 7 and the references given therein), who underlines that the great parameter sensitivity and the significantly fewer degrees of freedom of the one-dimensional system might make its stability characteristics, time histories, and responses to external disturbances very different from higher-dimensional models. Without a careful investigation of the connection between three-dimensional circulations and the flow in one-dimensional loops (which is beyond the scope of the present paper), quantitative conclusions about the real ocean cannot be drawn based on the findings of this study. However, interpreting the thermohaline loop as a “metaphor for the circulation instead of the circulation itself” (Wunsch 2005, p. 97), some basic understanding about the circulation’s fundamental properties can indeed be gained. In the present case, the asymmetry between the interior and the surface flow was shown to considerably affect loop dynamics and properties, especially with respect to the influence of the nonlinear effects of cabbeling and thermobaricity.
The relevance of this result can be tested by comparing our standard parameter settings to observations. Because of the highly idealized nature of the loop model, it could symbolize various kinds of overturning cells. We will here focus on the Atlantic meridional overturning circulation but propose similar scalings for Labrador Sea convective plumes and the wind-driven shallow overturning circulation of the subtropical Atlantic in Table 3.
Suitable parameter settings for application of the thermohaline loop to Labrador Sea convective plumes (Phenomenon 1) and the shallow overturning circulation of the subtropical Atlantic (Phenomenon 2). Refer, for example, to Klinger et al. (1996) and Marshall and Schott (1999) for characteristics of convective plumes and to Schott et al. (2004) and Boccaletti (2005) for those of the North Atlantic subtropical shallow overturning cells. The notation is as in the main text; we furthermore introduce τo for the overturning time scale, Θ+ for the temperature estimate of the upper, and Θ− for that of the lower branch of the cell.

In the North Atlantic, the mean thermocline temperature is typically about Θo = 10°C, and we estimate 20°C for the temperature in the upper thermocline (i.e., Θ = Θo + 10°C) and 5°C below the thermocline (i.e., Θ = Θo − 5°C). This roughly corresponds to the temperature distribution in the loop model when the warm section is half as long as the cold section, a case achieved when Zf = 0.5, using the nondimensionalization Θ = Θo + Θ′ × 1°C. With a circulation time scale of 600 yr (Toggweiler and Key 2003) and 100–200 days as commonly accepted values for the frictional decay time scale of the deep ocean, the Grashof number is on the order of F ~ 10−3. The inertialess assumption thus appears reasonable.
Using a circulation time scale of 600 yr and a vertical length scale of aυ ~ 103 m, we obtain a vertical diffusivity of κυ ~ 10−5 m2 s−1 for R = 0.1, a number in general agreement with observations (Toole et al. 1994). Choosing instead a length scale that corresponds to the global overturning’s horizontal length scale ah ~ 107 m produces an equivalent horizontal diffusivity of κh ~ 103 m2 s−1 for R = 0.1. Hence, it is arguably a reasonable zero-order approximation to use a single nondimensional number for representing horizontal and vertical diffusivities. We hypothesize that it may be possible to slightly modify the loop model using an elliptic rather than a circular shape to represent the large aspect ratio characterizing large-scale ocean overturning cells; however, we do not attempt to develop such a model extension here. Our hypothesis is comforted by the fact that we find similar inverse Rayleigh numbers for the horizontal and vertical components of the overturning circulations under consideration (cf. Table 3). In a generalized loop model, it would be useful to be able to apply different Rayleigh numbers for horizontal and vertical circulations, similar to the approach by Hazewinkel et al. (2012) in the context of stressed horizontal convection. On the other hand, even the discrimination between (fixed) horizontal and vertical diffusivities would only yield a poor description of oceanic turbulence, which is inhomogeneous in both space and time (Polzin et al. 1997; Whalen et al. 2012). By mixing density in the vertical, this turbulence is an important contributor to driving the global overturning circulation (e.g., Munk and Wunsch 1998; Visbeck 2007; Talley 2013), and in the past years, efforts have been made to develop more elaborate parameterizations that, for example, involve the sources of this mixing [e.g., breaking internal gravity waves, as in Müller and Natarov (2003) and Olbers and Eden (2013)].
The cabbeling parameter
6. Summary and outlook
The asymmetry between the essentially horizontal surface flow and the interior circulation, which spans both horizontal and vertical scales, is a fundamental characteristic of many oceanic overturning cells. In this study, we have investigated the effect of this asymmetry on dynamics and fluid properties in the context of a thermohaline loop. To that end, we have introduced a simple modification of the loop model devised by Wunsch (2005): by folding the loop at the level where the point sources and sinks of temperature and salinity are applied, until the upper part is purely horizontal, this asymmetry can be accounted for in a simple way. In such a setting, the forcing is concentrated at the uppermost level of the loop, which is in good agreement with the observation that in reality the buoyancy fluxes mainly occur at the ocean’s surface.
In our analysis, we focused on weak diffusion (small values of the inverse Rayleigh number R) and a weakly nonlinear equation of state (relatively small nonlinear parameters). The circulation was purely thermally driven (no salinity or wind stress forcing applied), with the source and sink of heat located at the same level to best represent the case of horizontal convection. Folding the loop led to significant changes in the transitory behavior and steady-state properties, including velocity oscillations and considerably higher temperature amplitudes associated with a slower circulation.
The influence of the two nonlinear EOS effects of cabbeling and thermobaricity on loop dynamics and steady-state properties differed significantly in the two model versions; for example, cabbeling was found to dampen the velocity oscillations observed in the folded loop, and its effect on the equilibrium velocity was reversed compared to the circular loop. This is because in the folded loop there is no compensation between density torques arising above and below the level of the heat forcing—an unrealistic property of the symmetric, circular loop model since the surface flow is basically horizontal in the real ocean. Folding the loop also induced a modification of the system’s stability properties, with dramatically different dynamic responses when large cabbeling parameters were used (cf. Fig. 6b). A detailed investigation of these differences will be the subject of future work, as the focus here was on weakly nonlinear forms of the EOS.
From the study of the one-dimensional fluid loop, only a basic, qualitative understanding of real ocean dynamics can be gained. Because of the wide range of parameter choices, different interpretations are possible; a useful analogy with the MOC, however, is encouraged by the observation of the following features in the folded loop model:
- a thermocline just below the warm source, with a typical thickness of a few hundred of meters, below which cold waters are sitting;
- a deep-reaching, cold vertical column, homogeneous almost to the bottom, below the cold source, reminiscent of the cold weakly stratified polar regions; and
- a surface region where the temperature is gradually decreasing poleward in a transition zone between a large body of warm surface waters and the region of deep convection.
We thank Jonas Nycander, Carl Wunsch, and two anonymous reviewers for their helpful comments. This work has been funded by the Bolin Centre for Climate Research (RA1).
REFERENCES
Arzel, O., , T. Huck, , and A. Colin de Verdière, 2006: The different nature of the interdecadal variability of the thermohaline circulation under mixed and flux boundary conditions. J. Phys. Oceanogr., 36, 1703–1718, doi:10.1175/JPO2938.1.
Boccaletti, G., 2005: Timescales and dynamics of the formation of a thermocline. Dyn. Atmos. Oceans, 39, 21–40, doi:10.1016/j.dynatmoce.2004.10.010.
Bryan, F., 1986: High-latitude salinity effects and interhemispheric thermohaline circulations. Nature, 323, 301–304, doi:10.1038/323301a0.
Dewar, W., , and R. Huang, 1996: On the forced flow of salty water in a loop. Phys. Fluids, 8, 954–970, doi:10.1063/1.868874.
Hazewinkel, J., , F. Paparella, , and W. Young, 2012: Stressed horizontal convection. J. Fluid Mech., 692, 317–331, doi:10.1017/jfm.2011.514.
Held, I. M., 2005: The gap between simulation and understanding in climate modeling. Bull. Amer. Meteor. Soc., 86, 1609–1614, doi:10.1175/BAMS-86-11-1609.
Hieronymus, M., , and J. Nycander, 2013: The buoyancy budget with a nonlinear equation of state. J. Phys. Oceanogr., 43, 176–186, doi:10.1175/JPO-D-12-063.1.
Huang, R. X., 1999: Mixing and energetics of the oceanic thermohaline circulation. J. Phys. Oceanogr., 29, 727–746, doi:10.1175/1520-0485(1999)029<0727:MAEOTO>2.0.CO;2.
Hughes, G. O., , and R. W. Griffiths, 2008: Horizontal convection. Annu. Rev. Fluid Mech., 40, 185–208, doi:10.1146/annurev.fluid.40.111406.102148.
IOC, SCOR, and IAPSO, 2010: The International Thermodynamic Equation of Seawater—2010: Calculation and use of thermodynamic properties. Intergovernmental Oceanographic Commission, Manuals and Guides 56, 220 pp. [Available online at http://www.teos-10.org/pubs/TEOS-10_Manual.pdf.]
Klinger, B. A., , J. Marshall, , and U. Send, 1996: Representation of convective plumes by vertical adjustment. J. Geophys. Res., 101, 18 175–18 182, doi:10.1029/96JC00861.
Klocker, A., , and T. J. McDougall, 2010: Influence of the nonlinear equation of state on global estimates of dianeutral advection and diffusion. J. Phys. Oceanogr., 40, 1690–1709, doi:10.1175/2010JPO4303.1.
Malkus, W. V., 1972: Non-periodic convection at high and low Prandtl number. Mem. Soc. Roy. Sci. Liege, 4, 125–128.
Manabe, B. S., , and R. J. Stouffer, 1999: Are two modes of thermohaline circulation stable? Tellus, 51A, 400–411, doi:10.1034/j.1600-0870.1999.t01-3-00005.x.
Marotzke, J., 2000: Abrupt climate change and thermohaline circulation: Mechanisms and predictability. Proc. Natl. Acad. Sci. USA, 97, 1347–1350, doi:10.1073/pnas.97.4.1347.
Marshall, J., , and F. Schott, 1999: Open-ocean convection: Observations, theory, and models. Rev. Geophys., 37, 1–64, doi:10.1029/98RG02739.
McDougall, T. J., 1987: Thermobaricity, cabbeling, and water-mass conversion. J. Geophys. Res., 92, 5448–5464, doi:10.1029/JC092iC05p05448.
McDougall, T. J., 2003: Potential enthalpy: A conservative oceanic variable for evaluating heat content and heat fluxes. J. Phys. Oceanogr., 33, 945–963, doi:10.1175/1520-0485(2003)033<0945:PEACOV>2.0.CO;2.
McDougall, T. J., , and C. J. Garrett, 1992: Scalar conservation equations in a turbulent ocean. Deep-Sea Res., 39A, 1953–1966, doi:10.1016/0198-0149(92)90007-G.
McDougall, T. J., , D. R. Jackett, , and F. J. Millero, 2009: An algorithm for estimating Absolute Salinity in the global ocean. Ocean Sci. Discuss., 6, 215–242, doi:10.5194/osd-6-215-2009.
Miljkovic, N., , and E. N. Wang, 2011: Modeling and optimization of hybrid solar thermoelectric systems with thermosyphons. Sol. Energy, 85, 2843–2855, doi:10.1016/j.solener.2011.08.021.
Müller, P., , and A. Natarov, 2003: The internal wave action model (IWAM). Near-Boundary Processes and Their Parameterization: Proc. 13th ‘Aha Huliko‘a Winter Workshop, Honolulu, HI, University of Hawai’i at Mānoa, 95–105.
Munk, W., , and C. Wunsch, 1998: Abyssal recipes II: Energetics of tidal and wind mixing. Deep-Sea Res. I, 45, 1977–2010, doi:10.1016/S0967-0637(98)00070-3.
Olbers, D., , and C. Eden, 2013: A global model for the diapycnal diffusivity induced by internal gravity waves. J. Phys. Oceanogr., 43, 1759–1779, doi:10.1175/JPO-D-12-0207.1.
Polzin, K., , J. Toole, , J. Ledwell, , and R. Schmitt, 1997: Spatial variability of turbulent mixing in the abyssal ocean. Science, 276, 93–96, doi:10.1126/science.276.5309.93.
Roquet, F., 2013: Dynamical potential energy: A new approach to ocean energetics. J. Phys. Oceanogr., 43, 457–476, doi:10.1175/JPO-D-12-098.1.
Ruddick, B., , and L. Zhang, 1996: Qualitative behavior and nonoscillation of Stommel’s thermohaline box model. J. Climate, 9, 2768–2777, doi:10.1175/1520-0442(1996)009<2768:QBANOS>2.0.CO;2.
Schanze, J. J., , and R. W. Schmitt, 2013: Estimates of cabbeling in the global ocean. J. Phys. Oceanogr., 43, 698–705, doi:10.1175/JPO-D-12-0119.1.
Schott, F. A., , J. P. Mccreary, , and G. C. Johnson, 2004: Shallow overturning circulations of the tropical-subtropical oceans. Earth’s Climate, Geophys. Monogr., Vol. 147, Amer. Geophys. Union, 261–304.
Stommel, H., 1961: Thermohaline convection with two stable regimes of flow. Tellus, 13A, 224–230, doi:10.1111/j.2153-3490.1961.tb00079.x.
Talley, L. D., 2003: Shallow, intermediate, and deep overturning components of the global heat budget. J. Phys. Oceanogr., 33, 530–560, doi:10.1175/1520-0485(2003)033<0530:SIADOC>2.0.CO;2.
Talley, L. D., 2013: Closure of the global overturning circulation through the Indian, Pacific, and Southern Oceans: Schematics and transports. Oceanography, 26, 80–97, doi:10.5670/oceanog.2013.07.
Toggweiler, J., , and R. Key, 2003: Thermohaline circulation. Encyclopedia of Atmospheric Sciences, A. C. Judith and J. A. Pyle, Eds., Elsevier, 1549–1555.
Toole, J. M., , R. W. Schmitt, , and K. L. Polzin, 1994: Estimates of diapycnal mixing in the abyssal ocean. Science, 264, 1120–1123, doi:10.1126/science.264.5162.1120.
Vallis, G. K., 2006: Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation. Cambridge University Press, 745 pp.
Visbeck, M., 2007: Oceanography: Power of pull. Nature, 447, 383–383, doi:10.1038/447383a.
Welander, P., 1967: On the oscillatory instability of a differentially heated fluid loop. J. Fluid Mech., 29, 17–30, doi:10.1017/S0022112067000606.
Welander, P., 1986: Thermohaline effects in the ocean circulation and related simple models. Large-Scale Transport Processes in Oceans and Atmosphere, J. Willebrand and D. L. T. Anderson, Eds., Springer, 163–200.
Whalen, C., , L. Talley, , and J. MacKinnon, 2012: Spatial and temporal variability of global ocean mixing inferred from Argo profiles. Geophys. Res. Lett., 39, L18612, doi:10.1029/2012GL053196.
Wunsch, C., 2005: Thermohaline loops, Stommel box models, and the Sandström theorem. Tellus, 57A, 84–99, doi:10.1111/j.1600-0870.2005.00093.x.
Young, W. R., 2010: Dynamic enthalpy, Conservative Temperature, and the seawater Boussinesq approximation. J. Phys. Oceanogr., 40, 394–400, doi:10.1175/2009JPO4294.1.
Yuan, S., , and C. Wunsch, 2005: Stress-driven thermohaline loops. Phys. Fluids, 17, 066601, doi:10.1063/1.1927887.
Zvirin, Y., 1982: A review of natural circulation loops in pressurized water reactors and other systems. Nucl. Eng. Des., 67, 203–225, doi:10.1016/0029-5493(82)90142-X.
Note that also the lowest portion of the loop could be folded into the horizontal plane since many overturning cells feature both a horizontal surface and a horizontal bottom branch. However, as the bulk of the buoyancy forcing is confined to the ocean’s surface, little dynamical effect is expected from a folded lower branch, which is why we restrict ourselves to folding the loop’s upper part only.
Note that from here on we refer to the nondimensional density and mass anomalies simply as density and mass for the sake of brevity. This is the reason why “negative densities” are sometimes obtained.