Vertical Stability

Storyboard

The stability of the marine water column depends on both temperature and salinity.

If the temperature increases, the water expands, creating a zone of lower density, causing the volume to tend to float.

Conversely, if salinity increases, density increases, causing the volume to sink.

In this sense, there is a competition between the effects of temperature and salinity, where the volume can try to rise or sink. The latter case is key for the generation of deep currents.

>Model

ID:(1524, 0)



Mechanisms

Iframe

>Top



Code
Concept
Diffusion concept
Instability in water in case of temperature difference
Lambda factor
Lewis number
Rayleigh number for temperature and stability
Stability condition
Stability of the seawater column
Stability of the water column
Temperature and salinity variation

Mechanisms

ID:(15508, 0)



Stability of the water column

Description

>Top


Typically, the density of seawater increases with depth.



This means that the layers closer to the surface are lighter than the deeper layers. This ensures that these layers float above the deeper layers and do not tend to displace them.

However, fluctuations in temperature and salinity can cause the deeper layers to become less dense than the upper layers. This creates an unstable situation, as these layers tend to float and emerge above the upper layers.

It is only in situations where the density is constant or increases with depth that the system is stable.

On the other hand, when a system becomes unstable, it means that in the presence of a disturbance, it may collapse, but if left undisturbed, it can maintain its current state.

ID:(12045, 0)



Temperature and salinity variation

Top

>Top


The increase in the temperature variation ($\Delta T$) leads to thermal expansion, causing the volume variation by temperature ($\Delta V_T$) to increase relative to the volumen ($V$) with the thermic dilatation coefficient ($k_T$), as shown in:

$ k_T \equiv \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_T }{ \Delta T }$



Similarly, the increase in the variación de la salinidad ($\Delta s$) due to mass causes the density variation ($\Delta\rho$) to increase relative to the sea water density ($\rho$) with the salinity coefficient ($k_s$), as shown in:

$ k_s =\displaystyle\frac{ 1 }{ \rho }\displaystyle\frac{ \Delta\rho }{ \Delta s }$



This expression is equivalent to the expression in which the volume variation due to salinity ($\Delta V_s$) decreases (negative value), as shown in:

$ k_s \equiv - \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_s }{ \Delta s }$

Therefore, the role of temperature and salinity is crucial, as they can cause the oceanic water column to become unstable, leading a volume element to start floating or sinking, thereby reversing the column.

ID:(15514, 0)



Instability in water in case of temperature difference

Description

>Top


When heating water in a pot, a region of lower density forms at the bottom near the heat source. This region begins to rise, attempting to displace the cooler, denser layer above it, which in turn tends to sink.

Once the temperature difference between the surface and the bottom exceeds a critical value, true jets of hotter water start to emerge, reaching the surface and creating space for the cooler surface water to descend towards the bottom:

ID:(12046, 0)



Stability of the seawater column

Image

>Top


In the case of seawater, not only can there be variations in temperature, but also in salinity. Salinity generally increases density, so processes that reduce salinity in the depth can lead to instabilities.

In this case, there are regions where water with higher salinity sinks while water with lower concentration rises. These sinking salt zones are called salt fingers and can be seen in the following graph generated through simulation:

ID:(12051, 0)



Diffusion concept

Description

>Top


Diffusion corresponds to the random motion of molecules gradually spreading out in space. The multiple collisions frequently cause them to change their travel direction, leading to a slow expansion. To describe this movement, statistical concepts are employed, such as describing the region where the majority of particles are located using the root-mean-square deviation. In fact, this root-mean-square deviation increases linearly over time:

The proportionality constant is known as the diffusion coefficient.

The concept is also used to describe how properties of particles, such as momentum and energy, propagate within a system. In this case, the spatial distribution of particles remains unchanged, but the spatial distribution of the considered parameter is affected.

ID:(13405, 0)



Rayleigh number for temperature and stability

Top

>Top


When water is heated in a pot, the water near the bottom begins to warm up, causing it to expand by ($$) according to the thermal expansion relationship, which complies with the thermic dilatation coefficient ($k_T$), the volumen ($V$), and the temperature variation ($\Delta T$) through:

$ k_T \equiv \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_T }{ \Delta T }$



the buoyancy force ($F_b$) is proportional to the displaced volume and can be approximately expressed as:

$F_b \sim g \Delta V \sim k_T V \Delta T$



By analyzing the units, we can observe that the factor

$\Delta V g \rightarrow \displaystyle\frac{m^4}{s^2}$



is the square of a diffusion constant. Therefore, instability can be understood as the dominance of the diffusion constant of the moment ($D_p$) of convection compared to the thermal diffusion constant ($D_T$) necessary to increase temperature and the loss of momentum due to viscosity.

Therefore, if the following proportion:

$\displaystyle\frac{g \Delta V}{D_T D_p} = \displaystyle\frac{g k_T V}{D_p D_T} \Delta T$



is much greater than unity, convection will dominate. In this sense, it makes sense to define a characteristic dimensionless number known as the rayleigh number for temperature ($Ra_T$):

$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$

In the case of a system without boundaries, it has been shown that the critical limit for instability occurs when the Rayleigh number exceeds $Ra_L=657.51$. However, this limit depends on the geometry of the system, and in the case of a cylinder (like an open pot), it has been demonstrated to be unstable when $Ra_L=1,100.65$.

ID:(15510, 0)



Lambda factor

Top

>Top


The tendency for an element of oceanic water to float due to increased temperature or to sink due to increased salinity is represented in the following diagram:



To study the situation, we introduce the lambda factor ($\Lambda$) as the proportion of the rayleigh number for temperature ($Ra_T$) and the numero de Rayleigh para la salinidad ($Ra_s$):

$\Lambda = \displaystyle\frac{Ra_T}{Ra_s} = \displaystyle\frac{k_T \Delta T}{k_s \Delta s}$



Since the rayleigh number for temperature ($Ra_T$) depends on the gravitational Acceleration ($g$), the thermic dilatation coefficient ($k_T$), the temperature variation ($\Delta T$), the diffusion constant of the moment ($D_p$), and the thermal diffusion constant ($D_T$) through the equation:

$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$



and the numero de Rayleigh para la salinidad ($Ra_s$) depends on the salinity coefficient ($k_s$), the variación de la salinidad ($\Delta s$), and the particle diffusion constant ($D_N$) through the equation:

$ Ra_s \equiv\displaystyle\frac{ g k_s h ^3 }{ D_p D_N } \Delta s $



we obtain the relationship for the lambda factor ($\Lambda$) through:

$ \Lambda \equiv \displaystyle\frac{ k_T \Delta T }{ k_s \Delta s }$

ID:(15511, 0)



Lewis number

Top

>Top


The numero de Lewis ($Le$) compares the thermal diffusion constant ($D_T$), which depends on the ocean thermal conduction ($\lambda_T$), the specific heat ($c$), and the sea water density ($\rho$), as per:

$ D_T \equiv \displaystyle\frac{ \lambda_T }{ \rho c }$



with the particle diffusion constant ($D_N$), which depends on the particle mobility ($\mu$), the boltzmann constant ($k_B$), and the absolute temperature ($T$), as per:

$ D_N \equiv \mu k_B T $



Therefore, it is defined as:

$ Le \equiv \displaystyle\frac{ D_T }{ D_N }$

ID:(15512, 0)



Stability condition

Top

>Top


To keep the system stable, it's necessary that the diffusion of energy (temperature) and salinity don't generate a the buoyancy force ($F_b$) large enough to invert the column. This is achieved when the lambda factor ($\Lambda$) is greater than the numero de Lewis ($Le$).

Therefore, the system is stable if the following condition is met:

$ Le < \Lambda $

It's important to note that the number factor depends on temperature and salinity, so if these variables vary, the system may reach a point of instability.

ID:(15515, 0)



Model

Top

>Top



Parameters

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$k_B$
k_B
Boltzmann constant
J/K
$k_s$
k_s
Contracción por salinidad
=
$\rho$
rho
Density
kg/m^3
$\Delta\rho$
Drho
Density variation
kg/m^3
$\Delta T$
DT
Diferencia de temperatura
K
$D_p$
D_p
Diffusion constant of the moment
m^2/s
$g$
g
Gravitational Acceleration
m/s^2
$Le$
Le
Numero de Lewis
-
$Ra_s$
Ra_s
Numero de Rayleigh para la salinidad
=
$\lambda_T$
lambda_T
Ocean thermal conduction
W/m K
$D_N$
D_N
Particle diffusion constant
m^2/s
$\mu$
mu
Particle mobility
$Ra_T$
Ra_T
Rayleigh number for temperature
=
$\rho$
rho
Sea water density
kg/m^3
$c$
c
Specific heat
J/kgK
$D_T$
D_T
Thermal diffusion constant
m^2/s
$k_T$
k_T
Thermic dilatation coefficient
$\Delta s$
Ds
Variación de la salinidad
=
$\eta$
eta
Viscosity of ocean water
Pa s

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$T$
T
Absolute temperature
K
$h$
h
Depth
m
$\Lambda$
Lambda
Lambda factor
=
$\Delta T$
DT
Temperature variation
K
$\Delta V_T$
DV_T
Volume variation by temperature
m^3
$\Delta V_s$
DV_s
Volume variation due to salinity
m^3
$V$
V
Volumen
m^3

Calculations


First, select the equation: to , then, select the variable: to

Calculations

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

Variable Given Calculate Target : Equation To be used




Equations

#
Equation

$ D_N \equiv \mu k_B T $

D_N = mu * k_B * T


$ D_p \equiv \displaystyle\frac{ \eta }{ \rho }$

D_p = eta / rho


$ D_T \equiv \displaystyle\frac{ \lambda_T }{ \rho c }$

D_T = lambda_T /( rho * c )


$ k_s \equiv - \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_s }{ \Delta s }$

k_s = - DV_s /( Ds * V )


$ k_s =\displaystyle\frac{ 1 }{ \rho }\displaystyle\frac{ \Delta\rho }{ \Delta s }$

k_s = Drho /( Ds * rho )


$ k_T \equiv \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_T }{ \Delta T }$

k_T = DV_T /( DT * V )


$ \Lambda \equiv \displaystyle\frac{ k_T \Delta T }{ k_s \Delta s }$

Lambda = k_T * DT /( k_s * Ds )


$ Le < \Lambda $

Le < Lambda


$ Le \equiv \displaystyle\frac{ D_T }{ D_N }$

Le = D_T / D_N


$ Ra_s \equiv\displaystyle\frac{ g k_s h ^3 }{ D_p D_N } \Delta s $

Ra_s = g * k_s * Ds * h ^3/( D_p * D_N )


$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$

Ra_T = g * k_T * DT * h ^3/( D_p * D_T )

ID:(15509, 0)



Thermal expansion

Equation

>Top, >Model


To model convection, we must consider that the water near the base of the system heats up, and as a result, expands. This expansion is what ultimately leads to a decrease in density and, therefore, the tendency to float. To describe this, the thermic dilatation coefficient ($k_T$) is introduced, indicating the proportion by which the volume variation by temperature ($\Delta V_T$) expands relative to the volumen ($V$) due to the increase in the temperature variation ($\Delta T$).

Therefore, we have:

$ k_T \equiv \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_T }{ \Delta T }$

$\Delta T$
Temperature variation
$K$
8983
$k_T$
Thermic dilatation coefficient
$1/K$
9361
$\Delta V$
Volume variation by temperature
$m^3$
8982
$V$
Volumen
$m^3$
8984

ID:(12050, 0)



Density variation due to the effect of salinity

Equation

>Top, >Model


The increase in the temperature variation ($\Delta T$) results in thermal expansion, leading to an increase in the volume variation by temperature ($\Delta V_T$) with respect to the volumen ($V$) as a function of the thermic dilatation coefficient ($k_T$), as shown in:

$ k_T \equiv \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_T }{ \Delta T }$



Similarly, adding salt to water leads to an increase in the variación de la salinidad ($\Delta s$) relative to the sea water density ($\rho$) due to the increase in the variación de la salinidad ($\Delta s$) as a function of the salinity coefficient ($k_s$), as shown in:

$ k_s =\displaystyle\frac{ 1 }{ \rho }\displaystyle\frac{ \Delta\rho }{ \Delta s }$

$k_s$
Contracción por salinidad
$-$
8996
$\rho$
Density
$kg/m^3$
5342
$\Delta\rho$
Density variation
$kg/m^3$
8624
$\Delta s$
Variación de la salinidad
$-$
8993

ID:(12053, 0)



Volume variation due to salinity effect

Equation

>Top, >Model


The increase in the variación de la salinidad ($\Delta s$) leads to changes in the density variation ($\Delta\rho$) relative to the sea water density ($\rho$) with the salinity coefficient ($k_s$), as shown in:

$ k_s =\displaystyle\frac{ 1 }{ \rho }\displaystyle\frac{ \Delta\rho }{ \Delta s }$



This can be formulated in terms of the equivalent the volume variation due to salinity ($\Delta V_s$) with respect to the volumen ($V$), resulting in:

$ k_s \equiv - \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_s }{ \Delta s }$

$k_s$
Contracción por salinidad
$-$
8996
$\Delta s$
Variación de la salinidad
$-$
8993
$\Delta V_s$
Volume variation due to salinity
$m^3$
10299
$V$
Volumen
$m^3$
8984

As the sea water density ($\rho$) is equal to a mass $m$ divided by the volumen ($V$), expressed as:

$\rho =\displaystyle\frac{m}{V}$



If we differentiate this expression for a constant mass $m$, it results in a the density variation ($\Delta\rho$) as:

$\Delta\rho =-\displaystyle\frac{m}{V^2}\Delta V=-\displaystyle\frac{\rho}{V}\Delta V$



Therefore, the expression in the salinity coefficient ($k_s$) with the variación de la salinidad ($\Delta s$):

$ k_s =\displaystyle\frac{ 1 }{ \rho }\displaystyle\frac{ \Delta\rho }{ \Delta s }$



implies:

$ k_s \equiv - \displaystyle\frac{ 1 }{ V }\displaystyle\frac{ \Delta V_s }{ \Delta s }$

It's important to note that the sign is negative, meaning that the increase in salinity leads to what is effectively a reduction in volume, causing the volume to tend to sink.

ID:(15513, 0)



Diffusion constant of the moment

Equation

>Top, >Model


The movement of a system like water tends to dissipate until the system reaches rest relative to its surroundings. This phenomenon is known as viscosity and competes with the inertia of bodies to maintain movement.

The first term is associated with the viscosity of ocean water ($\eta$), while the second is related to mass, or in the case of a liquid, with the sea water density ($\rho$).

Therefore, we introduce the diffusion constant of the moment ($D_p$) with:

$ D_p \equiv \displaystyle\frac{ \eta }{ \rho }$

$D_p$
Diffusion constant of the moment
$m^2/s$
8985
$\rho$
Sea water density
$kg/m^3$
8605
$\eta$
Viscosity of ocean water
$Pa s$
8612



The units are:

$\displaystyle\frac{\eta}{\rho} \rightarrow \displaystyle\frac{Pa,s}{kg/m^3} = \displaystyle\frac{m^3 kg,m,s}{s^2m^2kg} = \displaystyle\frac{m^2}{s}$

which corresponds to a diffusion constant. The value for water is on the order of $10^{-6} , m^2/s$.

ID:(12049, 0)



Temperature diffusion constant

Equation

>Top, >Model


The temperature in a system like water tends to diffuse until it is uniform throughout the volume. This diffusion is proportional to the ocean thermal conduction ($\lambda_T$) and inversely proportional to the sea water density ($\rho$) and the specific heat ($c$), which are necessary to increase the temperature.

Therefore, we introduce the thermal diffusion constant ($D_T$) as:

$ D_T \equiv \displaystyle\frac{ \lambda_T }{ \rho c }$

$\lambda_T$
Ocean thermal conduction
$J/m s K$
8987
$\rho$
Sea water density
$kg/m^3$
8605
$c$
Specific heat
$J/kg K$
8988
$D_T$
Thermal diffusion constant
$m^2/s$
8989



The units are:

$\displaystyle\frac{\lambda_T}{\rho,c} \rightarrow \displaystyle\frac{J/m,s,K}{kg/m^3,J/kg K} = \displaystyle\frac{m^2}{s}$

which corresponds to a diffusion constant. The value for water is on the order of $10^{-6} , m^2/s$.

ID:(12048, 0)



Particle diffusion constant

Equation

>Top, >Model


The diffusion of particles, such as salt, occurs slowly due to the interaction of the particles with the medium. This process depends, on one hand, on the particle mobility ($\mu$), expressed in $(m/s)/N=kg/s$, which corresponds to the velocity a particle achieves when a force is applied. On the other hand, it depends on the absolute temperature ($T$), associated with the velocity the particle can reach.

Therefore, the particle diffusion constant ($D_N$) for the movement of molecules is:

$ D_N \equiv \mu k_B T $

$T$
Absolute temperature
$K$
8979
$k_B$
Boltzmann constant
$J/K$
8978
$D_N$
Particle diffusion constant
$m^2/s$
8977
$\mu$
Particle mobility
$s/kg$
8980



where $k_B=1.34\times 10^{-23} J/K$ is the boltzmann constant ($k_B$).

ID:(12054, 0)



Rayleigh number for temperature and stability

Equation

>Top, >Model


Stability depends on the buoyancy force ($F_b$), which is proportional to the volume variation by temperature ($\Delta V_T$), which along with the gravitational Acceleration ($g$), needs to be compared with the diffusion constant of the moment ($D_p$) and the thermal diffusion constant ($D_T$). If we rewrite the volume variation by temperature ($\Delta V_T$) in terms of thermal expansion with the thermic dilatation coefficient ($k_T$), where the volumen ($V$) is expressed as the cube of the depth ($h$), we get:

$\displaystyle\frac{g \Delta V}{D_T D_p} = \displaystyle\frac{g k_T V}{D_p D_T} \Delta T$



This way, we can define the rayleigh number for temperature ($Ra_T$) with respect to temperature:

$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$

$h$
Depth
$m$
10064
$\Delta T$
Diferencia de temperatura
$K$
8992
$D_p$
Diffusion constant of the moment
$m^2/s$
8985
$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$Ra_T$
Rayleigh number for temperature
$-$
8990
$D_T$
Thermal diffusion constant
$m^2/s$
8989
$k_T$
Thermic dilatation coefficient
$1/K$
9361

ID:(12047, 0)



Rayleigh number for salinity

Equation

>Top, >Model


The rayleigh number for temperature ($Ra_T$) represents the comparison of the volume variation by temperature ($\Delta V_T$) in terms of the temperature variation ($\Delta T$) and the thermic dilatation coefficient ($k_T$) with the thermal diffusion constant ($D_T$) and the diffusion constant of the moment ($D_p$):

$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$



with the gravitational Acceleration ($g$). Similarly, a relationship for salinity can be established by replacing the thermic dilatation coefficient ($k_T$) with the salinity coefficient ($k_s$) and the thermal diffusion constant ($D_T$) with the particle diffusion constant ($D_N$), resulting in the numero de Rayleigh para la salinidad ($Ra_s$):

$ Ra_s \equiv\displaystyle\frac{ g k_s h ^3 }{ D_p D_N } \Delta s $

$k_s$
Contracción por salinidad
$-$
8996
$h$
Depth
$m$
10064
$D_p$
Diffusion constant of the moment
$m^2/s$
8985
$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$Ra_s$
Numero de Rayleigh para la salinidad
$-$
8991
$D_N$
Particle diffusion constant
$m^2/s$
8977
$\Delta s$
Variación de la salinidad
$-$
8993

ID:(12055, 0)



Lambda factor

Equation

>Top, >Model


The key to determining whether the volume of water will tend to float or sink can be studied by comparing the relationship between the rayleigh number for temperature ($Ra_T$) and the numero de Rayleigh para la salinidad ($Ra_s$), allowing us to define a characteristic number called the lambda factor ($\Lambda$).

$\Lambda = \displaystyle\frac{Ra_T}{Ra_s}$



By using the relationships defining the Rayleigh numbers, it can be shown that the lambda factor ($\Lambda$) is a function of the thermic dilatation coefficient ($k_T$), the salinity coefficient ($k_s$), and the temperature variation ($\Delta T$) with the variación de la salinidad ($\Delta s$):

$ \Lambda \equiv \displaystyle\frac{ k_T \Delta T }{ k_s \Delta s }$

$k_s$
Contracción por salinidad
$-$
8996
$\Lambda$
Lambda factor
$-$
8976
$\Delta T$
Temperature variation
$K$
8983
$k_T$
Thermic dilatation coefficient
$1/K$
9361
$\Delta s$
Variación de la salinidad
$-$
8993

Since the rayleigh number for temperature ($Ra_T$) depends on the gravitational Acceleration ($g$), the depth ($h$), the temperature variation ($\Delta T$), the diffusion constant of the moment ($D_p$), and the thermal diffusion constant ($D_T$), as defined by:

$ Ra_T \equiv\displaystyle\frac{ g k_T h ^3 }{ D_p D_T } \Delta T$



and the numero de Rayleigh para la salinidad ($Ra_s$) depends on the salinity coefficient ($k_s$) and the variación de la salinidad ($\Delta s$), as defined by:

$ Ra_s \equiv\displaystyle\frac{ g k_s h ^3 }{ D_p D_N } \Delta s $



then, we can assert that

$\Lambda = \displaystyle\frac{Ra_T}{Ra_s}$



reduces to:

$ \Lambda \equiv \displaystyle\frac{ k_T \Delta T }{ k_s \Delta s }$

ID:(12056, 0)



Lewis number

Equation

>Top, >Model


The numero de Lewis ($Le$) compares the thermal diffusion constant ($D_T$) with the particle diffusion constant ($D_N$) through:

$ Le \equiv \displaystyle\frac{ D_T }{ D_N }$

$Le$
Numero de Lewis
$-$
8994
$D_N$
Particle diffusion constant
$m^2/s$
8977
$D_T$
Thermal diffusion constant
$m^2/s$
8989

ID:(12058, 0)



Stability condition

Equation

>Top, >Model


The system is stable as long as the lambda factor ($\Lambda$) is greater than the numero de Lewis ($Le$), as in this case, the diffusion of energy (temperature) and salinity fails to destabilize the column:

$ Le < \Lambda $

$\Lambda$
Lambda factor
$-$
8976
$Le$
Numero de Lewis
$-$
8994

ID:(12057, 0)