CC BY
11Original Russian Text O A. S. Samodurov, A. M. Chukharev, D. A. Kazakov, 2021, published in MORSKOY GIDROFIZICHESKIY ZHURNAL, Vol. 37, Iss. 4 (2021)
Basic Regularities of Vertical Turbulent Exchange in the Mixed and Stratified Layers of the Black Sea
A. S. Samodurov A. M. Chukharev, D. A. Kazakov
Marine Hydrophysical Institute of RAS, Sevastopol, Russian Federation H [email protected]
Purpose. The purpose of the study is to assess the coefficient of vertical turbulent exchange for different layers of the Black Sea basin based on the experimental data on micro structure of the physical fields obtained for the period 2004-2019 in the Black Sea and using the semi-empirical models.
Methods and Results. For the upper mixed layer, the turbulent energy dissipation rate e and the exchange coefficient were calculated using the velocity fluctuation spectra based on the Kolmogorov hypotheses on the turbulence spectrum inertial range. In the stratified layers, the turbulence coefficient and the dissipation rate were experimentally determined both from the spectra of the velocity horizontal fluctuations' gradients and the vertical spectra of temperature fluctuations using the concept of the effective scale of turbulent patches. Depending on the features of the hydrological regime and the prevailing energy contributors to turbulence generation, five layers were identified and described (including their characteristic power dependences of the vertical turbulent diffusion coefficients K on the buoyancy frequency N) using the 1.5D-model of vertical turbulent exchange for the basin under study. For the stratified layers, the 1.5D-model results were comparatively analyzed with those of the other semi-empirical and theoretical models describing the most probable hydrophysical processes in each specific layer; the relations for the vertical turbulent exchange coefficient were obtained depending on the buoyancy frequency. Conclusions. Comparison of the experimental data collected under different hydrometeorological conditions with the simulations resulted from the known turbulence models for the sea upper layer showed that the best agreement between the simulation and measurement data was provided by a multiscale model taking into account three basic mechanisms of turbulence generation: current velocity shear, instability of wave motions, and wave breaking. The turbulent exchange coefficient dependencies on depth are conditioned by the effect of the turbulence dominant source at a given level. In the stratified layers, the exchange coefficient dependence on buoyancy frequency is determined by the hydrophysical processes in each layer; the relations obtained for individual layers indicate intensity of the contributions of vertical advection, internal wave breakings, turbulence diffusion and geothermal flux.
Keywords: Black Sea, energy dissipation, stratified layer, vertical turbulent exchange, buoyancy frequency, measuring complex, field measurements, turbulent exchange models, 1.5D-model
Acknowledgements: the investigation was carried out within the framework of state task No. 05552021-0004 "Oceanological processes".
For citation: Samodurov, A.S., Chukharev, A.M. and Kazakov, D.A., 2021. Basic Regularities of Vertical Turbulent Exchange in the Mixed and Stratified Layers of the Black Sea. Physical Oceanography, [e-journal] 28(4), pp. 376-391. doi:10.22449/1573-160X-2021-4-376-391
DOI: 10.22449/1573-160X-2021-4-376-391
O A. S. Samodurov, A. M. Chukharev, D. A. Kazakov, 2021
O Physical Oceanography, 2021
Introduction
Over the past decades, the need for predictive assessments of the functioning of the oceans and seas ecosystems and the solution of a wide range of scientific and practical problems associated with this has increased. For this purpose, it is necessary to use the data on the nature and intensity of physical exchange
processes in the sea active layer to create adequate practical calculations models of the nutrients transfer, pollutants from submarine wastewater discharges and accidental emissions, to assess the deep-sea layers ventilation activity, etc. All of these processes depend in the most essential way on turbulent exchange, which, in its turn, is determined by the hydrological regime and can vary widely in time and space. Revealing the regularities that determine the distribution of vertical turbulent exchange coefficients is a necessary condition for objective assessments of the intensity of various substances vertical flows. These regularities can be further used to solve a wide range of oceanological problems.
The purpose of this work is to generalize the long-term results of experimental and theoretical studies of turbulent regime in various layers of the Black Sea in order to give an integral picture of turbulent exchange intensity dependence on environmental conditions and physical processes affecting the vertical mixing.
Experimental data
The collection of experimental data was carried out using a specialized measuring complex "Sigma-1", developed at Marine Hydrophysical Institute and designed for studying the microstructure of hydrophysical fields [1]. The complex allows registering a wide range of physical characteristics: current velocity fluctuations, temperature and electrical conductivity (including their fluctuation values), pressure. The appearance of the complex various modifications is demonstrated in Fig. 1. The studies in the upper stratified layers were carried out during expeditionary work at R/V in 2004-2019. When performing the probing, the instrument freely moves downward at a constant velocity of approximately 0.75 m/s. The nature of movement is recorded by the built-in accelerometers system, which makes it possible to estimate the natural oscillations of the carrier. Such control is necessary to ultimately achieve "pure" data on the velocity fluctuations in the medium. In this case, the carrier's own movements influence on the measurement results is eliminated by appropriate data processing.
Fig. 1. Measuring facilities "Sigma-1": on the left - for studying the hydrophysical fields' microstructure (probing version); on the right - for studying turbulent processes in the near-surface layer (positional version)
When studying the turbulent regime in the homogeneous upper layer, the positional version of "Sigma-1" is installed on a stationary oceanographic
PHYSICAL OCEANOGRAPHY VOL. 28 ISS. 4 (2021) 377
platform using a specially designed system 1. In this case, measurements of turbulent fluctuations are supplemented with information on the mean current, surface waves and meteorological data.
As an objective assessment of the turbulence intensity, as a rule, we used the turbulent energy dissipation rate values, which were determined differently for the positional and probing measurement options. In the first case, the rate of turbulent energy s dissipation was calculated by the method proposed in [2] and then described in [3]. At the same time, signal distortions caused by suspention device swell and vibrations do not significantly affect the result. The method is based on Kolmogorov's hypothesis, according to which the velocity fluctuations spectral density E can be represented as
E (k) = s1/4v5/4F (X),
where k is a wave number; v is a kinematic viscosity; F(X) is a universal function (model spectrum X = k / (s14v-34) is a dimensionless wave number.
Energy dissipation rate was determined by combining the velocity fluctuations experimental spectrum with the model one, which was the Nasmith spectrum [4, 5]. Our data comparison, namely, the calculated values of the dissipation rate when normalizing to the energy flux from the atmosphere to the waves, with the results of in situ measurements carried out using the WAVES research programs [6] on Lake Ontario and SWADE [7] in the North Atlantic, where the characteristics of turbulence near the air - water surface were studied, showed their good agreement.
In the second case, the dissipation rate was determined using the technique described in [8] with regard to MST (Microstructure-Turbulence) Profiler instrument equipped with PNS 93 gauge for measuring the velocity shifts.
The main difference of "Sigma-1" complex from the specified instrument is that the device directly measures the velocity fluctuations, separates the horizontal component and calculates the smoothed values of the speed fluctuation gradient Au / Az in a certain layer, according to which the turbulent energy dissipation rate is determined. In this case, vertical fluctuations are not quite informative due to the high vertical velocity of the instrument; therefore, a vertical gradient of horizontal fluctuations is used. For this purpose we first remove outliers and noise emissions due to natural oscillations of the instrument, and then use a band-pass filter to restrict the measured inertia frequency fluctuations by the subband of the turbulence spectrum from the low-frequency side (kl) and the Kolmogorov wavenumber
kc = (1/2 n)(s/v3 )1/4 at the upper boundary where v is a kinematic viscosity.
1 Barabash, V.A., Samodurov, A.S., Chukharev, A.M., 2015. [Measuring System for the Study of Small-Scale Turbulence in the Near-Surface Layer of the Sea]. Patent RF, no. 2014151917/93.
The dissipation rate is calculated by the iterative method: firstly, the boundary wave numbers ki and kmax are determined, then the spectrum of du/dz value is calculated by the Welch method (P. D. Welch). In this case, the records are divided into overlapping segments, which are multiplied by the Hann time window, then we perform the Fourier transform, followed by averaging the spectral functions
over all segments. The dispersion
(du / dz)2
is determined by integrating
the spectrum values in the selected range of wavenumbers. The dissipation rate of turbulent energy is determined by the relation
15
s = — V 2
(du / dz )2
This s value is used to calculate the Kolmogorov wave number. If the stop conditions are not met, the cycle is repeated starting from the spectrum calculation. The stopping criteria are a small variation of kc value (less than a step in M spectrum) and an excess of kc over kmax. An example of the calculated dimensionless spectrum in 72-75 m depth range is represented in Fig. 2. For comparison, the same figure demonstrates the Nasmyth Universal Spectrum [4].
F i g. 2. Comparison of the Nasmyth Universal Spectrum (blue line) and the experimental spectrum of the horizontal fluctuations' gradients (green line) at the 72-75 m depth. The fluctuations were preliminary smoothed by the median filter and subjected to a band-pass filtering. The dissipation rate is 5.3-10-8 m2-s-3
In order to estimate the coefficients of vertical turbulent diffusion and turbulent energy dissipation rate, a method based on the energy analysis of the evolution of turbulent (mixed) patches was also used. The effective scale of turbulent patches was estimated from measurements of temperature fluctuations and the vertical spectra calculated from them. The method is convenient when the density gradient is mainly determined by the temperature gradient contribution.
In [9, 10], it was developed a method for calculating the dissipation rate e dependence on the buoyancy frequency N using the effective scale of turbulent patches L and the dependence of vertical turbulent diffusion coefficient K on N:
8 = 0.1L2 N3, K = —'f— (1)
1 - R N2
where n = — — is a buoyancy frequency; g is a gravitational acceleration; p is
\p dz
density; Rf is flux Richardson number (the ratio of potential energy increase rate in the system to the rate of energy input spent on mixing) in acts of a stratified flux shear instability and breaking of wave disturbances. In [11-13], using various approaches, it was determined that Rf is constant for the phenomena under consideration. For application in calculations, several approximate Rf values were proposed (1/3 in [11] and 1/4 in [12]), as well as a value equal to 0.2 to estimate the common factor of the right-hand side of the second of formulas (1) -in [14-15]. Function e is the energy dissipation rate in the act of breaking per unit area.
An approach to determining L is based on the structure analysis of the spectra
AT'
of the first temperature fluctuation differences -measured in the ocean [16].
Az
It was previously determined that the effective vertical scale of natural turbulent patches corresponds to the vertical scale of a stable minimum in the vertical spectrum small-scale region of the first differences of fluctuations [17].
Turbulent mixing in the upper mixed layer
Vertical turbulent exchange in the upper homogeneous layer of the sea is determined by a large number of natural factors. The main mechanisms for turbulence generation in this layer are instability of wave motions induced by surface waves, wave breaking, drift current velocity shear, convection, Langmuir circulation, and a number of others of lesser importance. Most models describing turbulent exchange in a layer adjacent to the atmosphere are limited to two or three of the mechanisms listed above. However, in recent years, researchers have been paying more and more attention to Langmuir circulations, suggesting their importance [18-21].
Convective mixing can be considered as a separate mechanism of vertical exchange, which, under appropriate conditions, can dominate over turbulent, but this type of mixing is not taken into account in the models considered in this paper.
Experimental data on the turbulent regime in the upper mixed layer (with neutral stratification) collected under various hydrometeorological conditions provided the verification of the most well-known models of turbulent exchange for this layer. The analyzed models included the model of the near-wall (logarithmic) layer, the possibility of which was proposed in [22], the models described in [23-25], as well as the multiscale model [26] developed at MHI Turbulence Department.
At this stage of research, a number of conclusions can be drawn from the comparative analysis of models and experiments.
1. At weak winds, the calculations based on none of the models provide satisfactory agreement with the experimental data: all the calculation results are obtained significantly lower.
2. The curve calculated according to Craig - Banner model [23], at moderate and strong winds, in many cases quite well falls on the experimental points but the roughness parameter z0 has to be changed within very large limits. In some cases, in order to achieve agreement between the results of calculations and experiment, z0 parameter must be many times higher than the wave height. It turns out a paradoxical situation when at a greater wave height the roughness parameter should be much less than in cases with small wave amplitude.
3. In addition to the current velocity shift as a turbulence generation mechanism, Benilov and Ly model [24] also takes into account the energy inflow from the surface waves. However, the proposed calculation method does not provide the desired effect - the calculated curve just slightly differs from the logarithmic curve (wall-bounded turbulence models for rough walls), and only in the uppermost part of the considered 1-2 m thick layer, below it these curves practically coincide. In most cases, these simulation results did not fit well with the measurement data 2.
4. In the model of Kudryavtsev et al. [25], in most cases, the calculation results were less than the measurement data, which could be associated with the use of only one set of coefficients in the calculations, regardless of hydrometeorological conditions.
5. The results of calculations using the multiscale Chukharev model [26], as a rule, were in agreement with the experimental data, but in a number of cases a noticeable discrepancy was observed. The introduction of directly measured values of the wave characteristics and the current velocity profile into the model instead of parameterizations (especially in the upper 3 m layer) clearly improves the agreement between the results of calculations and experiments.
The results of various models verification indicate the need to include in the turbulent exchange models for the near-surface layer at least three main mechanisms of turbulence generation: drift current velocity shear, nonlinear effects of surface waves and their breaking. It is also important to note that the parameterizations of both the current velocity profile and the surface wave spectrum require further improvement, possibly taking into account the regional characteristics of the basin.
Thus, of the considered models, Chukharev model [26] turned out to be the most preferable for assessing the intensity of turbulent exchange in the near-surface layer of the sea. However, the model dependences used in this model (taken from the literature) of wave characteristics from wind velocity and current velocity profile do not always correspond to the actual values, which is reflected in the results.
2 Chukharev, A.M., 2014. [Contribution of the Main Mechanisms of Turbulence Generation to Vertical Exchange in the Sea Stirring Layer. Dr. Phys.-Math. Diss.]. Sevastopol: MHI RAS, 274 p.
An analysis of the relative contribution of various mechanisms to the overall generation of turbulence demonstrates that at different depths and under different hydrometeorological conditions, one or another mechanism can dominate. Thus, we again come to the conclusion that it is necessary to take into account all the main mechanisms, since neglecting any of them will distort the real picture.
The upper mixed layer turbulent exchange coefficient does not have a simple dependence on any one parameter and is a complex function of wind velocity, the degree of wave development, and others. In Fig. 3 an example of calculating K using a multiscale model with the input parameters shown in the figure are demonstrated. The calculation of K in the model is carried out through the turbulence energy and dissipation rate:
E2
K = C E^ , C = 0.09. ^ e
The coefficient decrease in the uppermost layer occurs very rapidly, i.e. the direct effect of surface waves and their breaking is relatively shallow. K(z) dependence is well described here by the power function K(z)= azn, where a = 0.005-0.04; n = -0.6...-0.9 at zo < z < 2HS; Hs is the height of significant waves. The values of numerical coefficients vary depending on the hydrometeorological conditions. For the lower layer (2HS < z < 20HS), a polynomial dependence of the form K(z) = a0 + aiz + a2z2 + a3z3 suits better, which is consistent with the KPP parameterizations usually used for the upper layer in global models. The turbulence generation by the velocity shear prevails here but the turbulence diffusion from the upper layer also affects.
10
l?
20
2?
30
________
\ K(z) = 0.03Z-0-94 /
R2 = 0.98 / 271,
V10= 10 M-c-1
z0 = 0.1 m
A = 0.25 M f =0.18Tu F
*
K(z) = 2.2e- 6z3-1.4e-4z2+1.9e-3z+0.011
j if2 = 0.9
0.02
0.04
0.06 0.08 K. uV1
0.1
0.12
0.14
F i g. 3. Dependence of the turbulent viscosity coefficient on the upper homogeneous layer depth: - calculation by the Chukharev model [24]; — — - coefficient approximation by a power function for the upper layer; ■ ■ ■ ■ - polynomial approximation for the lower layer. The figure shows the input parameter values of the model: wind speed at the height 10 m (F10), sea surface roughness parameter (z0), amplitude (A) and frequency of the wave spectral peak f)
Thus, the value of turbulent viscosity coefficient near the sea surface is a complex function of depth and significantly depends on the wind velocity, wave parameters, and other physical factors. In this case, the approximation of K below the wave layer is well described by a third-degree polynomial in z.
It is considered that the turbulence generated by surface waves and breakings is "faster", but the scale of eddies is smaller than in shear turbulence, where they can reach the dimensions comparable to the mixed layer thickness.
Turbulent exchange in the Black Sea stratified layers
The main structural characteristic of layers in natural basins is well known - it is the presence of vertical density gradients that prevent the development of large-scale turbulence, but at the same time create favorable conditions for the development of internal waves (the shear instability formation dominant mechanism). In this case, the value of local Richardson number at the boundary of interlayers moving with multidirectional velocity becomes less than the critical one, which leads to the occurrence of small-scale turbulence, also due to the overturning of internal waves with the formation of turbulent patches. Inside such patches, turbulence is developed and, as a consequence, leads to vertical diffusion. The paper also discusses other mechanisms that generate turbulence in various layers of the Black Sea basin.
The purpose of this part of work is to identify the hydrophysical mechanisms of the layers formation with different density gradients observed in the Black Sea general stratified layer. As a solution to this problem, the construction of a 1.5.D-model of vertical turbulent exchange for the studied stratified sea layer is presented. It includes the role of turbulent diffusion and vertical advection in different ways for various depth ranges. Based on this model, it is possible to identify differences in the vertical structure formation mechanisms for each layer. The first such model of vertical exchange for the Black Sea, including vertical advection and vertical turbulent diffusion, was constructed in [27, 28]. Later, some changes were made to the model [27] based on the results of the monograph [29]: the depth-averaged area of the studied basin and the mean temperature and salinity profiles in the upper stratified layer were refined. In this case, the method for calculating the depth distribution of the vertical diffusion coefficient for a long period from [27, 28] was used. Note that, in contrast to the previously mentioned works, with this approach, the sought functions are already calculated within 501750 m depth range [15].
In the last decade, using the mentioned models and based on rich experimental data, we have obtained a number of important results characterizing the features of vertical turbulent diffusion in the Black Sea upper stratified layers [30-33]. Experimental data collected during the expeditions at R/V Professor Vodyanitsky from 2016 to 2018 provided, on the basis of a theoretical model [34], the development of ideas about the regularities of vertical turbulent diffusion in the main pycnocline.
F i g. 4. Dependence of the vertical turbulent diffusion coefficient K on the buoyancy frequency N in the stratified layers 1-5 at 50-1750 m depths
As indicated, the buoyancy frequency N (Vaisala - Brunt frequency) is used as the main parameter characterizing the density stability. Determination of the vertical turbulent exchange coefficient dependence K(N) is the main task of the described analysis. Calculation results comparison by semi-empirical and theoretical models revealed the fact that K(N) coefficient dependence is best approximated by a power function, the exponent of which is determined by the nature of internal waves in a particular medium [15]. A visual representation of K on N dependence variability is given in Fig. 4 and 5. The first of them shows (black line) the averaged dependence of the vertical diffusion coefficient K(N) within the framework of 1.5.-model. White lines indicate the approximating power-law dependences we have determined for various stratified layers. In the intervals between the calculated dependences K ~ Na, there are black curves connecting two different selected dependences K(N). It should be noted that such a transition gap is practically absent between the dependences for layers 1 and 2 with negative a powers. K(z) dependences in the layers 1-5 we have selected within the framework of 1.5.-model are shown in the left part of Fig. 5. On the right, the figure shows the averaged depth distribution of buoyancy frequency N. In each separate layer, the power dependence K(N) is indicated (Fig. 4) with a close integer degree a at N. The following research results are presented in the table.
O.OE+OO 5.0E-06 1.0E-05 1.5E-05 2.0E-05 O.OE+OO 5.0E-03 1.0E-02 1.5E-02 2.0E-02
F i g. 5. Model distribution of the vertical turbulent diffusion coefficient K in the Black Sea stratified layers 1-5 (on the left), and dependence of the buoyancy frequency on depth (on the right) based on the average multi-year data (o - the extrema)
Dependences K(N), m2s in different models for the stratified area
of the Black Sea
Layers (according to Fig. 5) 1.5D-model Theoretical models Semi-empirical models
1. Lower stirring layer 3T0-9N-202 2T0-8N~1-8 [33] 2T0-8N-2 [9, 10]
2. Upper layer of the main pycnocline 2T0-7N_a96 5.6T0-5N-1 [15]
3. Lower layer of the main pycnocline 5T0-4№'88 1.6T0"2N [15]
4. The layer affected by a geothermal flow from the inclined bottom 8T0-14N~22 2^10-14N-2 [35]
5. Stratified boundary of the bottom layer 4T0N23
This table compares the dependences of vertical turbulent energy coefficient K on the buoyancy frequency N in various stratified layers of the Black Sea, obtained using 1.5.-model, as well as theoretical and semi-empirical models in earlier studies. Below we discuss the analysis results of five identified stratified layers of the basin under study.
The lower stirring layer 1 (cold intermediate) is located close to the upper stirring layer, which serves as a source of various disturbances for it. With a constantly maintained stratification in the layer under consideration, sources of disturbances from above form internal waves in it, which undergo breakings and maintain turbulent exchange. In addition, quasi-inertial internal waves (quasihorizontal stratified unstable currents) also contribute to the maintenance of vertical turbulent exchange in this layer due to local wave breaking and the formation of turbulent patches [34, 36]. In our work [33], the power function K(N) in the indicated layer was calculated based on the data analysis of high-resolution measuring complex "Sigma-1" (the probing version). As a result, it was revealed that the power-law dependences of the buoyancy frequency in this layer, calculated from the measurement data (N-18) and from the modeling results using 1.5.-model (N-2-02) (see table), are close. For the rest of the studied layers, there are no measurement data due to the fact that the working immersion depth of the measuring probe complex is limited to 300 m. At the same time, there is another semi-empirical method that confirms K(N) dependence in the specified layer: e xLN, Kx e • N_2(1). Hence, using the relation [9, 10], from formulas (1) we obtain the required dependence K = 2 10-8N2 m2^s-1. This value is given in the table together with the results of calculations for 1.5. and semi-empirical models in the layer under study.
The next two stratified layers (2 and 3) constitute the main pycnocline of the studied basin. As was determined in [15, 34], vertical turbulent exchange is formed here due to shear instability of ray quasi-inertial internal waves (quasihorizontal currents) in a stratified fluid. In contrast to the lower, weakly stratified layer, in the upper, strongly stratified layer, the characteristic scales of waves transferring their energy to the turbulence (as shown) also depend on the derivative of the buoyancy frequency function N(z). As a result, the following results were obtained for the upper and lower layers:
e x N31ON/dz\, K x N\dN/dz\(layer 2); e x N3, K x N (layer 3).
The variation in depth dependence is exclusively due to the significant difference in the values of the buoyancy frequency N (vertical density gradients) in the upper and lower parts of the main pycnocline. More precisely, the pycnocline model is a structure of two layers described by two different power-law dependences K(N) and the function N(z) that unites them (see Fig. 4 and the Table). The power-law dependence N(z) in the main pycnocline (sum of layers) has the form Nx f1. In this case, the relations for the upper layer 2 are constructed using the real dependence N(z): s x N, Kx N. The obtained result demonstrated
close power-law dependences of buoyancy frequencies N with the 1.5. model for both layers (see Fig. 5, Table).
In layer 4, the dependence of vertical turbulent diffusion coefficient K on the buoyancy frequency N, calculated within the framework of 1.5.-model, has the form K <x N~22 (see Figs. 4 and 5). The probable hydrophysical mechanism of turbulent exchange in this deep-water layer is presented in [35]. In this paper, a model of vertical exchange has been developed due to geothermal heat flux from the inclined bottom, causing intrusive layering in the lower stratified layer of the Black Sea. The basin shape in the model is assumed to be conical, locally varying in depth, and the circle radius is determined in accordance with the change in the natural local area of the basin. The general expression for the vertical diffusion coefficient is
K = 8 PN-2,
where ep is an energy dissipation rate (with regard to coefficient from formula (1)), here the bottom topography, bottom geothermal flux and the efficiency of heat engine are taken into account. As a result, the model dependence of vertical turbulent exchange coefficient has the form K~ 2-10~14N"2 m2-s-1 (see the Table). It should be noted that, in contrast to the model [35], which was constructed in the depth range from about 800 m to the upper boundary of the near-bottom mixed layer, the dependence in the model used in this experimental approach (K <x N~2) actually manifested itself only in a "layer" of about 200 m thickness (see Fig. 4 and 5). This indicates that in the studied depth range we observe the layers in which a joint formation mechanism of vertical turbulent exchange, which manifests itself as a "gap" in the depth between the layers where the approximating dependence is not determined. For example, the same "gap" is also noted between layers 3 and 4. In this case, it should be noted that the mechanism under consideration manifests itself in the upper part of layer 4 more noticeably than was observed in the main pycnocline (aggregate of layers 2-3). In the upper layer adjacent to layer 3, the function K(N) begins to increase (in contrast to layer 3, in which it decreased) from a depth of approximately 1000 m (the point is marked as an extremum in Fig. 5), but with a different dependence on N.
The last of the considered stratified layers - layer 5 - is located at 1600-1750 m depths, below the stratified layer 4, and precisely adjoins the bottom homogeneous layer, not shown in the figures. The vertical structure of the stratification of two layers, 4 and 5, is shown in Fig. 5. In [36], the stationary state of the bottom layer is considered. The heat, continuously coming from the bottom, maintains the stationarity of the homogeneous layer and, together with the dissolved salt constantly getting from the flowing lower Bosporus current, penetrates upward through a thin stratified layer 5 adjacent from above. The bottom layer can be represented as a heat engine, in which the only source of energy for maintaining the vertical exchange in the system is the bottom heat flux [37]. For such a system, it is fair to use the dependences of e and K on the buoyancy frequency N in accordance with relations (1). However, since currently there are no developed exchange models, L(N) dependence in layer 5 can be found using the vertical
turbulent diffusion coefficient from 1.5.D-model (Fig. 4). Then the desired dependence of turbulent patches scale takes the form
It should be noted that the expression L(N) is applicable to this situation. Indeed, the lower boundary of layer 5 is linked to the homogeneous layer upper boundary where the gradients should vanish. Consequently, on the lower boundary N = 0, and the solution obtained here Lb = 0, Kb = 0 satisfies this requirement. In this case, the presence of "turbulent patches" in the upper layers, naturally, does not turn the exchange at the boundaries to zero:
Summarizing the foregoing, we can say that the performed experimental and theoretical studies make it possible to make practical estimates of vertical mixing intensity in different layers of the Black Sea.
For the upper mixed layer, various models of turbulent exchange were verified and the limits of their applicability were determined for different hydrometeorological conditions. The dependence of the exchange coefficient on depth is determined by the effect of dominant turbulence source at a given horizon.
For stratified sea layers, theoretical and semi-empirical relations for the coefficient of vertical turbulent exchange depending on the buoyancy frequency are obtained. It should be noted that the results comparison of calculating
the dependence K(N), which has the form K = ANa, according to 1.5.-model and hydrophysical models, shows a significant difference in the factor A values. Differences in the values of the coefficients (factors) may be due to the peculiarities of the Black Sea bottom topography and unevenness in the exchange processes distribution of various origins. So, for layer 1, the factor differs by an order of magnitude, which is most likely due to the close layer location to the upper 50-m mixed layer, while for layers 2 and 3 the discrepancy in the coefficients increases by two orders of magnitude due to the limited region of quasi-horizontal internal waves' generation.
Based on the analysis of 1.5.-model, the dependence of vertical turbulent diffusion coefficient K on the buoyancy frequency N in five identified stratified layers of the Black Sea basin at depths from 50 to 1750 m was revealed. Based on the analysis of a number of previously constructed semi-empirical vertical exchange models, various physical mechanisms for maintaining vertical turbulent exchange in different layers of the basin have been identified. A comparison of the calculation results using 1.5.-model and semi-empirical models demonstrates that the powers of a at N in the dependence K ~ Na have similar values. When analyzing the results, the layers thickness was limited by the condition of power constancy in the layer, i.e. most layers did not touch each other. As it was found out, this is explained by the fact that different physical mechanisms of vertical exchange often act simultaneously in one stratified layer.
(I)L x N_1; (II)L x N~2/3; (III)L x №; (IV)L x N1.
Conclusions
In order to clarify the causes for the difference between the model and
experimental coefficients in K = AN dependences, it is necessary to increase the number of measurements with the determination of local and seasonal
differences in A coefficient and a power.
REFERENCES
1. Samodurov, A.S., Dykman, V.Z., Barabash, V.A., Efremov, O.I., Zubov, A.G., Pavlenko, O.I. and Chukharev, A.M., 2005. "Sigma-1" Measuring Complex for the Investigation of Small-Scale Characteristics of Hydrophysical Fields in the Upper Layer of the Sea. Physical Oceanography, 15(5), pp. 311-322. https://doi.org/10.1007/s11110-006-0005-1
2. Stewart, R.W. and Grant, H.L., 1962. Determination of the Rate of Dissipation of Turbulent Energy near the Sea Surface in the Presence of Waves. Journal of Geophysical Research, 67(8), pp. 3177-3180. https://doi.org/10.1029/JZ067i008p03177
3. Monin, A.S. and Ozmidov, R.V., 1985. Turbulence in the Ocean. D. Reidel Publishing Company, 248 p. doi:10.1007/978-94-009-5217-1
4. Oakey, N.S., 1982. Determination of the Rate of Dissipation of Turbulent Energy from Simultaneous Temperature and Velocity Shear Microstructure Measurements. Journal of Physical Oceanography, 12(3), pp. 256-271. https://doi.org/10.1175/1520-0485(1982)012<0256: DOTROD>2.0.CO;2
5. Chukharev, A.M. and Repina, I.A., 2012. Interaction of Atmosphere and Sea Boundary Layers on Small and Meso-Scales in a Coastal Zone. Morskoy Gidrofizicheskiy Zhurnal, (2), pp. 60-78 (in Russian).
6. Terray, E.A., Donelan, M.A., Agrawal, Y.C., Drennan, W.M., Kahma, K.K., Williams, A.J., Hwang, P.A. and Kitaigorodskii, S.A., 1996. Estimates of Kinetic Energy Dissipation under Breaking Waves. Journal of Physical Oceanography, 26(5), pp. 792-807. https://doi.org/10.1175/1520-0485(1996)026<0792:EOKEDU>2.0.CO;2
7. Drennan, W.M., Donelan, M.A., Terray, E.A. and Katsaros, K.B., 1996. Oceanic Turbulence Dissipation Measurements in SWADE. Journal of Physical Oceanography, 26(5), pp. 808815. https://doi.org/10.1175/1520-0485(1996)026<0808:OTDMIS>2.0.CO;2
8. Prandke, H. and Stips, A., 1998. Test Measurements with an Operational Micro structure-Turbulence Profiler: Detection Limit of Dissipation Rates. Aquatic Sciences, 60(3), pp. 191209. http://dx.doi.org/10.1007/s000270050036
9. Samodurov, A.S., Lubitsky, A.A. and Panteleev, N.A., 1995. Contribution of Breaking Internal Waves to Structure Formation, Energy Dissipation, and Vertical Diffusion in the Ocean. Physical Oceanography, 6(3), pp. 177-190. https://doi.org/10.1007/BF02197516
10. Samodurov, A.S. and Ivanov, L.I., 2003. Mixing and Energy Dissipation Rate in Mediterranean Seas: an Intercomparison of Existing Models. In: A. Yilmaz, ed., 2003. Oceanography of the Eastern Mediterranean and Black Sea: Similarities and Differences of Two Interconnected Basins. Ankara: Tubitak Publishers, pp. 369-375.
11. Thorpe, S.A., 1973. Experiments of Instability and Turbulence in a Stratified Shear Flow. Journal of Fluid Mechanics, 61(4), pp. 731-751. doi:10.1017/S0022112073000911
12. Mcewan, A.D., 1983. The Kinematics of Stratified Mixing through Internal Wavebreaking. Journal of Fluid Mechanics, 128, pp. 47-57. doi:10.1017/S0022112083000373
13. Mcewan, A.D., 1983. Internal Mixing in Stratified Fluids. Journal of Fluid Mechanics, 128, pp. 59-80. doi:10.1017/S0022112083000385
14. Osborn, T.R., 1980. Estimates of the Local Rate of Vertical Diffusion from Dissipation Measurements. Journal of Physical Oceanography, 10(1), pp. 83-89. doi: 10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2
15. Samodurov, A.S., 2016. Complementarity of Different Approaches for Assessing Vertical Turbulent Exchange Intensity in Natural Stratified Basins. Physical Oceanography, (6), pp. 32-42. doi: 10.22449/1573-160X-2016-6-32-42
16.
17.
18.
19.
20
21
22.
23
24
25.
26
27
28
29.
30
31
32
33
39
Samodurov, A.S., Ivanov, V.A., Belokopytov, V.N. and Kulsha, O.E., 2016. Modeling Average Annual Vertical Exchange in the Black Sea Cold Intermediate Layer. Proceses in Geomedia, (2), pp. 141-147 (in Russian).
Gregg, M.C., 1989. Scaling Turbulent Dissipation in the Thermocline. Journal of Geophysical Research: Oceans, 94(C7), pp. 9686-9698. doi:10.1029/JC094iC07p09686
McWilliams, J.C., Huckle, E., Liang, J.-H. and Sullivan, P.P., 2012. The Wavy Ekman Layer: Langmuir Circulations, Breaking Waves, and Reynolds Stress. Journal of Physical Oceanography, 42(11), pp. 1793-1816. doi:10.1175/JPO-D-12-07.1
Sullivan, P.P., McWilliams, J.C. and Patton, E.G., 2014. Large-Eddy Simulation of Marine Atmospheric Boundary Layers above a Spectrum of Moving Waves. Journal of the Atmospheric Sciences, 71(11), pp. 4001-4027. https://doi.org/10.1175/JAS-D-14-0095.1
Sullivan, P.P., Romero, L., McWilliams, J.C. and Melville, W.K., 2012. Transient Evolution of Langmuir Turbulence in Ocean Boundary Layers Driven by Hurricane Winds and Waves. Journal of Physical Oceanography, 42(11), pp. 1959-1980. https://doi.org/10.1175/JPO-D-12-025.1
Li, S., Li, M., Gerbi, G.P. and Song, J.-B., 2013. Roles of Breaking Waves and Langmuir Circulation in the Surface Boundary Layer of a Coastal Ocean. Journal of Geophysical Research: Oceans, 118(10), pp. 5173-5187. https://doi.org/10.1002/jgrc.20387 Csanady, G.T., 1984. The Free Surface Turbulent Shear Layer. Journal of Physical Oceanography, 14(2), pp. 402-411. https://doi.org/10.1175/1520-0485(1984)014<0402:TFSTSL>2.0.CO;2
Craig, P.D. and Banner, M.L., 1994. Modeling Wave-Enhanced Turbulence in the Ocean Surface Layer. Journal of Physical Oceanography, 24(12), pp. 2546-2559. https://doi.org/10.1175/1520-0485(1994)024<2546:MWETIT>2.0.CO;2
Benilov, A.Yu. and Ly, L.N., 2002. Modelling of Surface Waves Breaking Effects in the Ocean Upper Layer. Mathematical and Computer Modelling, 35(1-2), pp. 191-213. https://doi.org/10.1016/S0895-7177(01)00159-5
Kudryavtsev, V., Shrira, V., Dulov, V. and Malinovsky, V., 2008. On the Vertical Structure of Wind-Driven Sea Currents. Journal of Physical Oceanography, 38(10), pp. 2121-2144. https://doi.org/10.1175/2008JPO3883.1
Chukharev, A.M., 2013. Multitime Scale Model of Turbulence in the Sea Surface Layer. Izvestiya, Atmospheric and Oceanic Physics, 49(4), pp. 439-449. https://doi.org/10.1134/S0001433813040026
Ivanov, L.I. and Samodurov, A.S., 2001. The Role of Lateral Fluxes in Ventilation of the Black Sea. Journal of Marine Systems, 31(1-3), pp. 159-174. https://doi.org/10.1016/S0924-7963(01)00051-3
Samodurov, A.S. and Ivanov, L.I., 2002. Balance Model for Calculating Average Vertical Fluxes of Fluid, Heat, Salt and Dissolved Chemical Substances in the Black Sea Thermohaline. Morskoy Gidrofizicheskiy Zhurnal, (1), pp. 7-24 (in Russian). Ivanov, V.A. and Belokopytov, V.N., 2013. Oceanography ofthe Black Sea. Sevastopol, 210 p.
Samodurov, A.S., Chukharev, A.M., Nosova, A.V. and Globina, L.V., 2013. Internal Waves Intensification in the Shelf Break Region as a Factor of the Vertical Exchange Intensification. Fundamentalnaya i Prikladnaya Gidrofizika, 6(2), pp. 12-24 (in Russian).
Samodurov, A.S., Chukharev, A.M., Zubov, A.G. and Pavlenko, O.I., 2015. Structure-Formation and Vertical Turbulent Exchange in the Coastal Area of the Sevastopol Region. Physical Oceanography, (6), pp. 3-14. doi:10.22449/1573-160X-2015-6-3-14
Samodurov, A.S. and Chukharev, A.M., 2018. Vertical Turbulent Exchange Features in the Layer of Seasonal Pycnocline in the Northwestern Part of the Black Sea. Journal of Physics: Conference Series, 1128, 012148. doi:10.1088/1742-6596/1128/1/012148 Samodurov, A.S. and Chukharev, A.M., 2018. The Vertical Turbulent Exchange Features in the Black Sea Active Layer. In: V. Karev, D. Klimov, K. Pokazeev, eds., 2018. Physical and Mathematical Modeling of Earth and Environment Processes. Cham: Springer, pp. 148-156. doi:10.1007/978-3-319-77788-7_16
34. Samodurov, A.S. and Globina, L.V., 2012. Dissipation of Energy and Vertical Exchange in Stratified Basins Caused by the Shear Instability in the Field of Quasiinertial Internal Waves. Physical Oceanography, 21(6), pp. 383-393. https://doi.org/10.1007/s11110-012-9130-1
35. Samodurov, A.S., 2003. Intrusive Lamination and Vertical Exchange in the Black Sea because of Geothermal Heat Flow on a Sloping Bottom. In: MHI, 2003. Ekologicheskaya Bezopasnost' Pribrezhnykh i Shelfovykh Zon i Kompleksnoe Ispol'zovanie Resursov Shel'fa [Ecological Safety of Coastal and Shelf Zones and Comprehensive Use of Shelf Resources]. Sevastopol: ECOSI-Gidrofizika. Iss. 8, pp. 152-156 (in Russian).
36. Wunsch, C. and Ferrari, R., 2004. Vertical Mixing, Energy, and the General Circulation of the Oceans. Annual Review of Fluid Mechanics, 36, pp. 281-314. https://doi.org/10.1146/annurev.fluid.36.050802.122121
37. Samodurov, A.S., 2009. Bottom Boundary Layer in the Black Sea: Formation of the Stationary State. Physical Oceanography, 19(1), pp. 13-21. https://doi.org/10.1007/s11110-009-9035-9
About the authors:
Anatoliy S. Samodurov, Chief Research Associate, Head of Turbulence Department, Marine Hydrophysical Institute of RAS (2 Kapitanskaya St., Sevastopol, 299011, Russian Federation), Dr. Sci. (Phys.-Math.), ORCID ID: 0000-0002-9910-5325, [email protected]
Aleksandr M. Chukharev, Leading Research Associate, Marine Hydrophysical Institute of RAS (2 Kapitanskaya St., Sevastopol, 299011, Russian Federation), Dr. Sci. (Phys.-Math.), ORCID ID: 0000-0003-1078-6425, [email protected]
Dmitry A. Kazakov, Junior Research Associate, Marine Hydrophysical Institute of RAS (2 Kapitanskaya St., Sevastopol, 299011, Russian Federation), Postgraduate, ORCID ID: 0000-00015083-4968, [email protected]
Contribution of the authors:
Anatoly S. Samodurov - research task setting, mathematical model development, review of literature on the vertical turbulent exchange models for stratified layers of the Ocean, choice of data analysis methods, formulation and analysis of the results obtained in the study, their interpretation, discussion of the work results, formulation of conclusions, article writing
Alexander M. Chucharev - research task setting, mathematical model development, review of literature on the vertical turbulent exchange models for near-surface quasi-homogeneous layer of the Ocean, software development and calculations, plotting, qualitative analysis of the research results and their description, formulation of conclusions, article writing
Dmitrii A. Kazakov - data processing and analysis, verification of vertical turbulent exchange models in stratified ocean layers, calculation results analysis and their interpretation, plotting, discussion of the work results, formulation of conclusions, article writing and editing
All the authors have read and approved the final manuscript.
The authors declare that they have no conflict of interest.
