Thermohaline staircase

Last updated

Thermohaline staircases are patterns that form in oceans and other bodies of salt water, characterised by step-like structures observed in vertical temperature and salinity profiles; the patterns are formed and maintained by double diffusion of heat and salt. The ocean phenomenon consists of well-mixed layers of ocean water stacked on top of each other. The well-mixed layers are separated by high-gradient interfaces, which can be several meters thick. The total thickness of staircases ranges typically from tens to hundreds of meters. [1]

Contents

Two types of staircases are distinguished. Salt-fingering staircases can be found at locations where relatively warm, salty water overlies relatively colder, fresher water. Here, large-scale temperature and salinity both increase upward, making the mixing process of salt fingering possible. [2] Locations where you can find these type of staircases are for example beneath the Mediterranean outflow, in the Tyrrhenian Sea, and northeast Caribbean. [1] Diffusive staircases can be found at locations where both temperature and salinity increase downward, for example in the Arctic Ocean and in the Weddell Sea. [2] An important feature of thermohaline staircases is their extreme stability in space and time. They can persist several years or more and can extend for hundreds of kilometers. [1] The interest in thermohaline staircases is partly due to the fact that the staircases represent mixing hot spots in the main thermocline. [2]

Extensive definition and detection

Schematic of a typical vertical temperature profile with staircases, indicating the definitions of the quantities used to detect the thermohaline staircases (green: mixed layer; orange: interface, the part of the water column between two mixed layers). The black line is used to verify the result using the Turner angle
(
T
u
)
{\displaystyle (Tu)}
and density ratio
(
R
r
)
{\displaystyle (R_{\rho })}
. Here,
D
T
M
L
,
D
S
M
L
{\displaystyle \Delta T_{ML},\Delta S_{ML}}
, and
D
s
M
L
{\displaystyle \Delta \sigma _{ML}}
represent the conservative temperature, salinity, and density variations in the mixed layer.
D
T
I
F
,
D
S
I
F
{\displaystyle \Delta T_{IF},\Delta S_{IF}}
and
D
s
I
F
{\displaystyle \Delta \sigma _{IF}}
represent the conservative temperature, salinity, and density variations in the interface.
h
M
L
{\displaystyle h_{ML}}
, and
h
I
F
{\displaystyle h_{IF}}
denote the height of the mixed layer and of the interface respectively.
T
-
M
L
,
S
-
M
L
,
p
-
M
L
{\displaystyle {\overline {T}}_{ML},{\overline {S}}_{ML},{\overline {p}}_{ML}}
represent the averages of conservative temperature, salinity and pressure in the mixed layer. Figure 1 Typical Profile Staircase.png
Schematic of a typical vertical temperature profile with staircases, indicating the definitions of the quantities used to detect the thermohaline staircases (green: mixed layer; orange: interface, the part of the water column between two mixed layers). The black line is used to verify the result using the Turner angle and density ratio . Here, , and represent the conservative temperature, salinity, and density variations in the mixed layer. and represent the conservative temperature, salinity, and density variations in the interface. , and denote the height of the mixed layer and of the interface respectively. represent the averages of conservative temperature, salinity and pressure in the mixed layer.

To determine the presence of thermohaline staircases, the following steps can be taken according to the algorithm designed by Van der Boog. [3]

The first step of the algorithm is to identify the mixed layers by locating weak vertical density gradients in conservative temperature and absolute salinity. To do so, the threshold gradient method is used with a threshold of , with the pressure and the reference pressure. The vertical conservative temperature, absolute salinity, and potential density gradients are all below the threshold value by meeting these three conditions:

with the thermal expansion coefficient, the haline contraction coefficient, the reference density, the conservative temperature, and the salinity.

The second step is to define the interface, which is the part of the water column in the middle of two mixed layers. It is required that the conservative temperature, absolute salinity, and potential density variations in the interface should be larger than the variations within each mixed layer to ensure a stepped structure. Therefore the following conditions should be met:

where subscript 1 corresponds to the mixed layer above the interface and subscript 2 corresponds to the mixed layer below the interface.

The third step is to limit the interface height . The interface height should be smaller than the height of the mixed layers directly above and below the interface . This condition has to be met in order to ensure that the interface is relatively thin compared to the mixed layers surrounding it. Furthermore, the algorithm removes all interfaces with conservative temperature or absolute salinity inversions to make sure that it only detects step-like structures that are associated with the presence of thermohaline staircases.

The fourth step is to determine the double-diffusive regime (salt-fingering or diffusive) of each interface. When both conservative temperature and absolute salinity of the mixed layers above and below the interface increase downward, the interface belongs to the diffusive regime. When both conservative temperature and absolute salinity of the mixed layers above and below the interface both increase upward, the interface is classified as the salt-fingering regime.

Finally, only vertical sequences of at least two interfaces in the same double-diffusive regime are selected, where the interfaces should be separated from each other by only one mixed layer. This way, most thermohaline intrusions are removed, as these are characterised by alternating mixed layers in the diffusive and salt-finger regimes. Furthermore, the algorithm removes salt-fingering interfaces and diffusive-convective interfaces outside their favourable Turner angle , a parameter used to describe the local stability of an inviscid water column. Interfaces with salt-fingering characteristics should correspond to Turner angles of and interfaces with diffusive-convective characteristics should correspond to Turner angles of .

Staircase origin

The origin of thermohaline staircases relies on double diffusive convection, and specifically on the fact that heated water diffuses more readily than salty water. However, there is still much debate on which specific mechanism of layering plays a role. [2] Six possible mechanisms are described below. [4]

Collective instability mechanism

This mechanism, involving collective instability, relies on the idea that after a period of active internal wave motion, layers appear. This hypothesis was motivated by laboratory experiments in which staircases formed from the initially uniform temperature and salinity gradients. [5] Growing waves might overturn and generate the stepped structure of thermohaline staircases. [6]

Thermohaline intrusion mechanism

This hypothesis states that staircases represent the final stage in the evolution of thermohaline intrusions. [1] Intrusions can evolve either to a state consisting of alternating salt-finger and diffusive interfaces separated by convecting layers, which is common at high density ratio , or to a series of salt-finger interfaces when the density ratio is low . This proposition relies on the presence of lateral property gradients to drive interleaving. This mechanism, where thermohaline intrusions are transformed into staircases, are likely to exist in strong temperature-salinity fronts. [4]

Metastable equilibria mechanism

A different theory states that staircases represent distinct metastable equilibria. It is suggested that finite amplitude perturbations to the gradient state force the system into a layered regime where it can remain for long periods of time. [5] Large initial perturbations to the gradient state make the transition to the staircase more likely and accelerate the process. Once the staircase is created, the system becomes resilient to further structural changes. [4]

Applied flux mechanism

The applied flux mechanism was mainly tested in laboratory experiments, and is most likely at work in cases when layering is caused by geothermal heating. [7] [8] When a stable salinity gradient is heated from below, top-heavy convection will take place in the lower part of the water column. The well-mixed convecting layer is bounded from above by a thin high-gradient interface. By a combination of molecular diffusion and entrainment across the interface, heat is transferred upward from the convecting layer. The molecular transfer of heat exceeds that of salt, resulting in a supply of buoyancy to the region immediately above the interface. This leads to the formation of a second convecting layer. The process can repeat itself over and over, which results in a sequence of mixed layers separated by sharp interfaces, a thermohaline staircase. [4]

Negative density diffusion

In salt-fingering staircases, vertical temperature and salinity fluxes are downgradient, while the vertical density flux is upgradient. This is explained by the fact that the potential energy released in transporting salt downward must exceed that expended in transporting heat upward, resulting in a net downward transport of mass. This negative diffusion sharpens the fluctuations and therefore suggests a means for generating and maintaining staircases. [1]

Instability of flux-gradient laws

This mechanism is based on negative density diffusion as well. However, instead of combining temperature and salinity into a single density term, it treats both density components individually. In a publication by Radko, [9] it is shown that formation of steps in numerical models is caused by the parametric variation of the flux ratio as a function of the density ratio , leading to an instability of equilibrium with uniform stratification. These unstable perturbations continuously grow in time until well-defined layers are formed. [4]

Observations

Salt-fingering staircases in the Tyrrhenian Sea: (a) The Tyrrhenian Sea; the white dot indicates the station where data have been collected. (b, c) Thermohaline staircases envelope from 2003 to 2016, both for conservative temperature (degC) and absolute salinity (g/kg), in the central Tyrrhenian station. Figure 2, Tyrrhenian Sea staircases.png
Salt-fingering staircases in the Tyrrhenian Sea: (a) The Tyrrhenian Sea; the white dot indicates the station where data have been collected. (b, c) Thermohaline staircases envelope from 2003 to 2016, both for conservative temperature (°C) and absolute salinity (g/kg), in the central Tyrrhenian station.

Two types of staircases exist: salt-fingering staircases, where both temperature and salinity of the mixed layers decrease with pressure (and therefore with depth); and diffusive staircases, where both temperature and salinity of the mixed layers increase with pressure (so with depth). [3]

Salt-fingering staircases

Most observations of salt-fingering staircases have come from three locations: the western Tropical Atlantic, the Tyrrhenian Sea, and the Mediterranean outflow. In these regions the density ratio has a very low value, which appears to be a condition for sufficient staircase formation. No staircases have been reported for values below 2. For values below 1.7, the step-like structures in vertical temperature and salinity profiles becomes apparent. [4] Moreover, the spatial pattern of staircases is very sensitive to . With decreasing , the height of steps sharply increases and the staircases become more pronounced. The importance of the density ratio for the formation is a sign that staircases are a product of double diffusive convection. [4]

Acoustic observations of a diffusive staircase in the Arctic Ocean. (a) Processed echogram with layer depths derived from the echogram scatter strength (white circles). (b) Vertical potential temperature with reference at the surface
(
Th
)
{\displaystyle (\Theta )}
and salinity profiles with black horizontal lines indicating the depth of the individual layers identified in the echogram. FIGWIKI2.png
Acoustic observations of a diffusive staircase in the Arctic Ocean. (a) Processed echogram with layer depths derived from the echogram scatter strength (white circles). (b) Vertical potential temperature with reference at the surface and salinity profiles with black horizontal lines indicating the depth of the individual layers identified in the echogram.

In the Tyrrhenian Sea, thermohaline staircases due to salt fingers are observed. [10] The step-like shape is visible in the vertical temperature and salinity profiles. Staircases in the Tyrrhenian Sea show a very high stability in space and time. The weak deep circulation in this area might be an explanation for this stability. [12]

Diffusive staircases

Diffusive staircases are found at higher latitudes. In the Arctic Ocean, warm and salty water from the Atlantic enters the Arctic basin and subducts beneath the colder and fresher waters of the upper Arctic. In some regions, also Pacific waters sit below the mixed layer and above the Atlantic layer. A thermocline is found at the top of the Atlantic Water layer. In that region, temperature and salinity increases with depth and step-like patterns are observed in vertical temperature and salinity profiles. [13] [11] These staircases mediate the heat transport from the warm water of Atlantic origin to the Arctic halocline and therefore serve as an important process in determining the heat flux from the Atlantic Water upward to the sea ice. Staircases in the Arctic are characterised by much smaller steps than in salt-fingering staircases. [4]

On a much smaller scale, diffusive staircases have also been observed in low- and mid-latitudes. For example, Lake Kivu and Lake Nyos show characteristic staircase patterns. [14] [15] In these salt-water lakes, geothermal springs supply heat at the bottom resulting in the diffusive background stratification. [4]

See also

Related Research Articles

<span class="mw-page-title-main">Fick's laws of diffusion</span> Mathematical descriptions of molecular diffusion

Fick's laws of diffusion describe diffusion and were first posited by Adolf Fick in 1855 on the basis of largely experimental results. They can be used to solve for the diffusion coefficient, D. Fick's first law can be used to derive his second law which in turn is identical to the diffusion equation.

<span class="mw-page-title-main">Convection</span> Fluid flow that occurs due to heterogeneous fluid properties and body forces.

Convection is single or multiphase fluid flow that occurs spontaneously due to the combined effects of material property heterogeneity and body forces on a fluid, most commonly density and gravity. When the cause of the convection is unspecified, convection due to the effects of thermal expansion and buoyancy can be assumed. Convection may also take place in soft solids or mixtures where particles can flow.

The Prandtl number (Pr) or Prandtl group is a dimensionless number, named after the German physicist Ludwig Prandtl, defined as the ratio of momentum diffusivity to thermal diffusivity. The Prandtl number is given as:

In fluid mechanics, the Grashof number is a dimensionless number which approximates the ratio of the buoyancy to viscous forces acting on a fluid. It frequently arises in the study of situations involving natural convection and is analogous to the Reynolds number.

In fluid mechanics, the Rayleigh number (Ra, after Lord Rayleigh) for a fluid is a dimensionless number associated with buoyancy-driven flow, also known as free (or natural) convection. It characterises the fluid's flow regime: a value in a certain lower range denotes laminar flow; a value in a higher range, turbulent flow. Below a certain critical value, there is no fluid motion and heat transfer is by conduction rather than convection. For most engineering purposes, the Rayleigh number is large, somewhere around 106 to 108.

<span class="mw-page-title-main">Boundary layer</span> Layer of fluid in the immediate vicinity of a bounding surface

In physics and fluid mechanics, a boundary layer is the thin layer of fluid in the immediate vicinity of a bounding surface formed by the fluid flowing along the surface. The fluid's interaction with the wall induces a no-slip boundary condition. The flow velocity then monotonically increases above the surface until it returns to the bulk flow velocity. The thin layer consisting of fluid whose velocity has not yet returned to the bulk flow velocity is called the velocity boundary layer.

The Richardson number (Ri) is named after Lewis Fry Richardson (1881–1953). It is the dimensionless number that expresses the ratio of the buoyancy term to the flow shear term:

The Hagen number (Hg) is a dimensionless number used in forced flow calculations. It is the forced flow equivalent of the Grashof number and was named after the German hydraulic engineer G. H. L. Hagen.

Spiciness (τ) is a term in oceanography that defines the salinity and potential temperature variation, often at constant density. Here, a temperature change offsets a salinity change; an increase in temperature decreases density whereas an increase in salinity increases density. Warmer and more saline water is spicier whereas cooler and less saline water is mintier.

There are two different Bejan numbers (Be) used in the scientific domains of thermodynamics and fluid mechanics. Bejan numbers are named after Adrian Bejan.

Ocean stratification is the natural separation of an ocean's water into horizontal layers by density, which is generally stable because warm water floats on top of cold water, and heating is mostly from the sun, which reinforces that arrangement. Stratification is reduced by wind-forced mechanical mixing, but reinforced by convection. Stratification occurs in all ocean basins and also in other water bodies. Stratified layers are a barrier to the mixing of water, which impacts the exchange of heat, carbon, oxygen and other nutrients. The surface mixed layer is the uppermost layer in the ocean and is well mixed by mechanical (wind) and thermal (convection) effects. Climate change is causing the upper ocean stratification to increase.

<span class="mw-page-title-main">Thermal effusivity</span> Ability of a material to exchange thermal energy with surroundings

In thermodynamics, a material's thermal effusivity, also known as thermal responsivity, is a measure of its ability to exchange thermal energy with its surroundings. It is defined as the square root of the product of the material's thermal conductivity and its volumetric heat capacity or as the ratio of thermal conductivity to the square root of thermal diffusivity.

Salt fingering is a mixing process, example of double diffusive instability, that occurs when relatively warm, salty water overlies relatively colder, fresher water. It is driven by the fact that heated water diffuses more readily than salty water. A small parcel of warm, salty water sinking downwards into a colder, fresher region will lose its heat before losing its salt, making the parcel of water increasingly denser than the water around it and sinking further. Likewise, a small parcel of colder, fresher water will be displaced upwards and gain heat by diffusion from surrounding water, which will then make it lighter than the surrounding waters, and cause it to rise further. Paradoxically, the fact that salinity diffuses less readily than temperature means that salinity mixes more efficiently than temperature due to the turbulence caused by salt fingers.

<span class="mw-page-title-main">Double diffusive convection</span> Convection with two density gradients

Double diffusive convection is a fluid dynamics phenomenon that describes a form of convection driven by two different density gradients, which have different rates of diffusion.

Paleosalinity is the salinity of the global ocean or of an ocean basin at a point in geological history.

The density ratio of a column of seawater is a measure of the relative contributions of temperature and salinity in determining the density gradient. At a density ratio of 1, temperature and salinity are said to be compensated: their density signatures cancel, leaving a density gradient of zero. The formula for the density ratio, R, is:

<span class="mw-page-title-main">Turner angle</span>

The Turner angleTu, introduced by Ruddick(1983) and named after J. Stewart Turner, is a parameter used to describe the local stability of an inviscid water column as it undergoes double-diffusive convection. The temperature and salinity attributes, which generally determine the water density, both respond to the water vertical structure. By putting these two variables in orthogonal coordinates, the angle with the axis can indicate the importance of the two in stability. Turner angle is defined as:

The Atlantic meridional overturning circulation (AMOC) is a large system of ocean currents, like a conveyor belt. It is driven by differences in temperature and salt content and it is an important component of the climate system. However, the AMOC is not a static feature of global circulation. It is sensitive to changes in temperature, salinity and atmospheric forcings. Climate reconstructions from δ18O proxies from Greenland reveal an abrupt transition in global temperature about every 1470 years. These changes may be due to changes in ocean circulation, which suggests that there are two equilibria possible in the AMOC. Stommel made a two-box model in 1961 which showed two different states of the AMOC are possible on a single hemisphere. Stommel’s result with an ocean box model has initiated studies using three dimensional ocean circulation models, confirming the existence of multiple equilibria in the AMOC.

A baroclinic instability is a fluid dynamical instability of fundamental importance in the atmosphere and ocean. It can lead to the formation of transient mesoscale eddies, with a horizontal scale of 10-100 km. In contrast, flows on the largest scale in the ocean are described as ocean currents, the largest scale eddies are mostly created by shearing of two ocean currents and static mesoscale eddies are formed by the flow around an obstacle (as seen in the animation on eddy. Mesoscale eddies are circular currents with swirling motion and account for approximately 90% of the ocean's total kinetic energy. Therefore, they are key in mixing and transport of for example heat, salt and nutrients.

<span class="mw-page-title-main">Faroe-Bank Channel overflow</span> Overflow current from Nordic Seas towards North Atlantic Ocean

Cold and dense water from the Nordic Seas is transported southwards as Faroe-Bank Channel overflow. This water flows from the Arctic Ocean into the North Atlantic through the Faroe-Bank Channel between the Faroe Islands and Scotland. The overflow transport is estimated to contribute to one-third of the total overflow over the Greenland-Scotland Ridge. The remaining two-third of overflow water passes through Denmark Strait, the Wyville Thomson Ridge (0.3 Sv), and the Iceland-Faroe Ridge (1.1 Sv).

References

  1. 1 2 3 4 5 Merryfield, William J. (2000-05-01). "Origin of Thermohaline Staircases". Journal of Physical Oceanography. 30 (5): 1046–1068. Bibcode:2000JPO....30.1046M. doi: 10.1175/1520-0485(2000)030<1046:OOTS>2.0.CO;2 . ISSN   0022-3670.
  2. 1 2 3 4 Radko, T.; Bulters, A.; Flanagan, J. D.; Campin, J.-M. (2014-05-01). "Double-Diffusive Recipes. Part I: Large-Scale Dynamics of Thermohaline Staircases". Journal of Physical Oceanography. 44 (5): 1269–1284. Bibcode:2014JPO....44.1269R. doi:10.1175/JPO-D-13-0155.1. hdl: 1721.1/91958 . ISSN   0022-3670. S2CID   66801727.
  3. 1 2 3 van der Boog, Carine G.; Koetsier, J. Otto; Dijkstra, Henk A.; Pietrzak, Julie D.; Katsman, Caroline A. (2021-01-13). "Global dataset of thermohaline staircases obtained from Argo floats and Ice-Tethered Profilers". Earth System Science Data. 13 (1): 43–61. Bibcode:2021ESSD...13...43V. doi: 10.5194/essd-13-43-2021 . ISSN   1866-3508. S2CID   234356971. CC-BY icon.svg Text was copied from this source, which is available under a Creative Commons Attribution 4.0 International License.
  4. 1 2 3 4 5 6 7 8 9 Radko, Timour (2013). Double-Diffusive Convection. Cambridge: Cambridge University Press. doi:10.1017/cbo9781139034173. ISBN   978-0-521-88074-9. S2CID   67802029.
  5. 1 2 Stern, Melvin E.; Stewart Turner, J. (1969-11-01). "Salt fingers and convecting layers". Deep Sea Research and Oceanographic Abstracts. 16 (5): 497–511. Bibcode:1969DSRA...16..497S. doi:10.1016/0011-7471(69)90038-2. ISSN   0011-7471.
  6. Stern, Melvin E. (January 1969). "Collective instability of salt fingers". Journal of Fluid Mechanics. 35 (2): 209–218. Bibcode:1969JFM....35..209S. doi:10.1017/S0022112069001066. ISSN   1469-7645. S2CID   121945515.
  7. Turner, J. S.; Stommel, Henry (1964-07-01). "A New Case of Convection in the Presence of Combined Vertical Salinity and Temperature Gradients". Proceedings of the National Academy of Sciences. 52 (1): 49–53. Bibcode:1964PNAS...52...49T. doi: 10.1073/pnas.52.1.49 . ISSN   0027-8424. PMC   300570 . PMID   16591193.
  8. Huppert, Herbert E.; Linden, P. F. (December 1979). "On heating a stable salinity gradient from below". Journal of Fluid Mechanics. 95 (3): 431–464. Bibcode:1979JFM....95..431H. doi:10.1017/S0022112079001543. ISSN   1469-7645. S2CID   123813620.
  9. Radko, Timour (2003-12-25). "A mechanism for layer formation in a double-diffusive fluid". Journal of Fluid Mechanics. 497: 365–380. Bibcode:2003JFM...497..365R. doi:10.1017/S0022112003006785. S2CID   119737327.
  10. 1 2 Durante, S.; Schroeder, K.; Mazzei, L.; Pierini, S.; Borghini, M.; Sparnocchia, S. (2019). "Permanent Thermohaline Staircases in the Tyrrhenian Sea". Geophysical Research Letters. 46 (3): 1562–1570. Bibcode:2019GeoRL..46.1562D. doi: 10.1029/2018GL081747 . ISSN   1944-8007. S2CID   134607170.
  11. 1 2 Stranne, Christian; Mayer, Larry; Weber, Thomas C.; Ruddick, Barry R.; Jakobsson, Martin; Jerram, Kevin; Weidner, Elizabeth; Nilsson, Johan; Gårdfeldt, Katarina (2017-11-09). "Acoustic Mapping of Thermohaline Staircases in the Arctic Ocean". Scientific Reports. 7 (1): 15192. Bibcode:2017NatSR...715192S. doi:10.1038/s41598-017-15486-3. ISSN   2045-2322. PMC   5680214 . PMID   29123176.
  12. Tait, R.I.; Howe, M.R. (1968-06-01). "Some observations of thermo-haline stratification in the deep ocean". Deep Sea Research and Oceanographic Abstracts. 15 (3): 275–280. Bibcode:1968DSRA...15..275T. doi:10.1016/0011-7471(68)90005-3. ISSN   0011-7471.
  13. Brown, Justin M.; Radko, Timour (2021-06-01). "Diffusive Staircases in Shear: Dynamics and Heat Transport". Journal of Physical Oceanography. 51 (6): 1915–1928. Bibcode:2021JPO....51.1915B. doi: 10.1175/JPO-D-20-0193.1 . ISSN   0022-3670. S2CID   233632601.
  14. Newman, Fred C. (1976-03-01). "Temperature Steps in Lake Kivu: A Bottom Heated Saline Lake". Journal of Physical Oceanography. 6 (2): 157–163. Bibcode:1976JPO.....6..157N. doi: 10.1175/1520-0485(1976)006<0157:TSILKA>2.0.CO;2 . ISSN   0022-3670.
  15. Schmid, Martin; Lorke, Andreas; Dinkel, Christian; Tanyileke, Gregory; Wüest, Alfred (August 2004). "Double-diffusive convection in Lake Nyos, Cameroon". Deep Sea Research Part I: Oceanographic Research Papers. 51 (8): 1097–1111. Bibcode:2004DSRI...51.1097S. doi:10.1016/j.dsr.2004.02.010.