Vertical stratification or vertical stability of the ocean was originally quantified by the gradient of density. However, previous researchers have noted that it is not suitable to define the vertical stability of the water column by using only the vertical gradient of potential density (e.g.,
McDougall, 1987). Currently, the squared Brunt-Väisälä frequency
$ {N}^{2} $ is primarily used to describe ocean stratification, and it has been widely applied in many studies.
Helber et al. (2012) estimated the transition layer thickness as the distance between the mixed layer depth and the deepest
$ {N}^{2} $ value, which is equal to half of the maximum
$ {N}^{2} $ that occurs in the profile. In
Maes (2008), a method was proposed that makes use of the dependency of the squared Brunt-Väisälä frequency
$ {N}^{2}(T,S) $. By separating the impact of temperature from
$ {N}^{2} $, the layer where salinity stratification plays a key role in
$ {N}^{2} $ can be isolated. This method has been applied in other studies by dividing
$ {N}^{2} $ into
$ {N}_{S}^{2} $ and
$ {N}_{T}^{2} $.
Zheng and Zhang (2012) indicated that during the 2007/2008 La Niña event, the positive salinity anomaly in the western-central Pacific played a more significant role than negative temperature anomaly in destabilizing the upper layers.
Maes and O’Kane (2014) identified global seasonal variability in ocean salinity stratification (OSS as defined by the upper-300 m averaged positive
$ {N}_{S}^{2} $) and illustrated that regions without permanent barrier layer may also existing strong salinity stratification (large OSS).
Guimbard et al. (2017) explained that
$ {N}_{S}^{2} $ can better explain the shoal of pycnocline than
$ {N}_{T}^{2} $ in the 2014/2015 rainy seasons over the central equatorial Pacific. Moreover, the role of
$ {N}^{2} $ variability on the baroclinic Rossby radius (
$ \lambda $) was studied by
Chelton et al. (1998), and they found that changes of
$ {N}^{2} $ had a limited contribution toward
$ \lambda $(1%–2%) at seasonal time scale.