Vol 40, No 2 (2024)

Thermohydrodynamics of the ocean and the atmosphere

Finite-Difference Approximation of the Potential Vorticity Equation for a Stratified Incompressible Fluid and an Example of its Application for Modeling the Black Sea Circulation Part I. Finite-Difference Equation of Potential Vorticity of Ideal Fluid

Demyshev S.G.

Abstract

Purpose. The study is purposed at deriving the discrete equations of absolute vorticity and potential vorticity for a three-dimensional stratified incompressible fluid as an exact consequence of the finite-difference equations of sea dynamics in the field of a potential mass force in the adiabatic approximation provided that viscosity and diffusion are absent. The purpose also consists in analyzing both the features of two-dimensional projections of the absolute vorticty equation onto the coordinate planes and the three-dimensional equation of potential vorticity.

Methods and Results. In order to determine the discrete analogues of absolute and potential vorticity, introduced is the additional grid on which the finite-difference equations for the components both of absolute vorticty and potential vorticity are written down. Two-dimensional analogues of the three-dimensional equation of absolute vorticity on the planes (x, y), (y, z) and (x, z) are obtained; they possess the feature of preserving vorticity, energy and enstrophy (square of vorticiy). A discrete equation for potential vorticity of a stratified incompressible fluid is derived from the finite-difference system of three-dimensional equations of sea dynamics in the adiabatic approximation at the absence of viscosity and diffusion.

Conclusions. In the case of a linear equation of state, the discrete equations of absolute vorticity and potential vorticity which are the exact consequence of finite-difference formulation are obtained. The equation of potential vorticity is of a divergent form, and two-dimensional analogues of the absolute vorticity equation on the planes (x, y), (y, z) and (x, z) have two quadratic invariants that provide preservation of the average wave number.

Morskoj gidrofizičeskij žurnal. 2024;40(2):165-179
pages 165-179 views

Generation of Vertical Fine Structure by Internal Waves on a Shear Flow

Slepyshev A.A., Ankudinov N.O.

Abstract

Purpose. The work is purposed at studying the wave mechanism of fine structure generation as well as at determining the vertical wave fluxes of mass.

Methods and Results. In contrast to the previously used mechanism of forming a fine structure by internal waves due to breaking, a new approach based on determining the vertical wave fluxes of mass in the field of an inertia-gravity internal wave without breaking is proposed. The inertia-gravity internal waves on the Black Sea northwestern shelf are considered on a current with a vertical velocity shift. The flow is assumed to be geostrophically balanced, with the vertical velocity shifts compensated by a horizontal density gradient. The f-plane approximation is used. Thus, the classical scheme for describing a wave field by the hydrodynamic equations is applicable with the regard for the nonlinear effects. A weakly nonlinear approach is used. In the linear approximation, the eigenfunction and the dispersion relation are found by solving numerically the boundary value problem which determines the vertical structure of a mode in the presence of an average flow. In this case, the wave frequency is of a complex character since the coefficients in the differential equation of the specified boundary value problem are complex. Depending on the wave period and the mode number, either weak attenuation of a wave or its weak amplification is possible. The eigenfunction of internal waves is also complex. Therefore, the vertical wave fluxes of mass and the vertical component of the Stokes drift velocity are nonzero and lead to the generation of vertical fine structure which is irreversible.

Conclusions. In the presence of a two-dimensional shear flow, taking into account the horizontal inhomogeneity of the average density field enhances the effect of generation of a vertical fine structure by the inertia-gravity internal waves. The vertical wave fluxes of mass also increase. The indicated fluxes and the generated fine structure for the waves of different frequencies are close, and the effect is enhanced in the presence of the waves of different frequencies.

Morskoj gidrofizičeskij žurnal. 2024;40(2):180-197
pages 180-197 views

Analysis of observations and methods of calculating Hydrophysical fields in the ocean

New Method for Determining Spectral Absorption of Light in the Sea

Lee M.E., Shybanov E.B.

Abstract

Purpose. The study is purposed at presenting and analyzing a new method for determining light absorption in the sea which for the first time made it possible to redirect almost all the scattered rays from the studied light beam to the photodetector along the path of its propagation in a weakly absorbing medium, as well as at showing that application of the new method providing such an efficient collecting of the scattered rays, permits not only to avoid significant errors from the strong influence of scattering upon the results of determining the light absorption, but also to give up the necessity in correcting the data by theoretical modeling.

Methods and Results. It is known that sea water is a weakly absorbing light-scattering medium in which the light propagation is accompanied by its attenuation that is many times stronger due to scattering than due to absorption. Therefore, the determination of light absorption by seawater at a receiving device requires collecting of not only the light that has traveled a certain distance in the absorption medium, but also all the light scattered along this path. Previously, proposed was the method for measuring the light absorption in a cylindrical mirror cuvette with a light source at the input and a collector with a photodetector at the output (reflective-tube absorption meter). Somewhat later, a similar method based on the phenomenon of total internal reflection was applied. Since these methods do not provide the sufficiently complete collection of scattered rays, the data are to be corrected by theoretical modeling. The authors propose a new method for determining the spectral absorption of light in a quartz glass conical cuvette with an external mirror cone. It is shown that the cone cuvette permits to collect most of the scattered rays in the beam passing through the water medium by means of more efficient redirection of these rays from the place of light scattering to the receiver. The rest part of the scattered rays that has left the cuvette reaches the receiver in the air space between the cuvette and the cone mirror due to multiple reflections from it. As a result, the new method makes it possible to redirect almost all of the scattered light to the receiver and thus to minimize the errors in determining the light absorption in a weakly absorbing medium. To quantify the advantages of the new method, the authors have calculated the geometric parameters of the scattered light propagation for a conical quartz cuvette in air and for the same cuvette placed inside the external cone mirror.

Conclusions. The combination of a conical quartz cuvette and an external mirror cone in the new method made it possible to collect all the rays scattered in a weakly absorbing medium in the receiver. Thus, it permitted not only to exclude their strong influence upon the determination of light absorption in the sea, but also to give up the necessity in correcting the data by theoretical modeling.

Morskoj gidrofizičeskij žurnal. 2024;40(2):198-214
pages 198-214 views

Interannual Salinity Changes in the Upper 1000-Meter Layer of Extratropical Zone in the Northwestern Pacific Ocean under Conditions of the Intensification of Global Hydrological Cycle

Rostov I.D., Dmitrieva E.V.

Abstract

Purpose. The study is purposed at determining the trends and the regional features of interannual changes in salinity and salt content in the upper 1000-m layer of extratropical zone in the northwestern Pacific Ocean, and at analyzing their possible cause-and-effect relations with large-scale and regional processes in the ocean and atmosphere over the last two decades of the current period of global warming.

Methods and Results. The NOAA climate data sets on salinity and current velocity in the nodes of the 0.3 × 1º regular grid at 31 levels for 2000–2022, derived from the oceanographic observation system GODAS (NCEP Global Ocean Data Assimilation System), as well as the series both of climate indices from the NOAA websites and routine salinity observations performed at the Rosgidromet coastal hydrometeorological stations were used in the study. Besides, average monthly ERA5 reanalysis data on precipitation (Reanalysis Data ERA5 monthly 2d Surface) and evaporation from the underlying surface (WHOI OAFlux version3 monthly evapr oaflux) obtained from the NOAA Oceanographic Data Access ERDDAP server were also applied. The process of analyzing involved the methods of cluster, correlation and regression analysis, as well as the apparatus of empirical orthogonal functions. The conducted research made it possible to identify and characterize the regional spatial and temporal features of the accelerated changes in salinity and salt content in the upper 1000-m water column of the study area under conditions of the current warming phase accompanied by the intensification of global and local hydrological cycles. The quantitative characteristics of the noted trends and their statistical significance were assessed.

Conclusions. On the whole, over the water area under study there is a tendency towards a gradual decrease in average salinity near the sea surface and water desalination in the upper 1000-m layer, which is accompanied by an increase of water heat content in this layer by 3 %.

Morskoj gidrofizičeskij žurnal. 2024;40(2):215-230
pages 215-230 views

Wave Nature and Modulation of Annual Fluctuations in the Level of the Baltic Sea

Zakharchuk E.A., Sukhachev V.N., Tikhonova N.A.

Abstract

Purpose. The study is purposed at assessing the features of spatial-temporal variability of the characteristics of annual fluctuations of the Baltic Sea level based on satellite and contact measurements, their comparing with theoretical dispersion relations of the low-frequency waves of different types, as well as at investigating possible mechanisms of the amplitude modulation of annual fluctuations of the Baltic Sea level.

Methods and Results. The hypothesis on a wave nature of annual fluctuations of the Baltic Sea level is tested and the reasons for their amplitude modulation are investigated based on the harmonic analysis of satellite altimetry data and the 132-year series of tide gauge sea level measurements in Stockholm. It is shown that the wave-like annual perturbations in the sea level field propagate from the southwest to the northeast at the velocity from 0.06 to 0.36 m/s. Having been compared, the estimated characteristics of annual waves and the theoretical dispersion ratios of the low-frequency waves of different types has shown that they are identified as the internal Kelvin waves over the most of the sea area, and only in rare cases in the southwest of the sea, their characteristics agree with the theoretical dispersion relations of the baroclinic topographic Rossby waves. The perceptible interdecadal changes in the annual wave parameters in the sea level field were noted. As compared to the 1993–2021period, the decade 1993–2002 is characterized by a decrease of the Sa harmonic amplitude by 1.5–3 times, by later onset of the maximum of the sea level annual variation (about 1 month later), and also by a noticeable slowdown of the annual wave phase velocity in the sea southwest.

Conclusions. The reasons for the amplitude modulation of annual waves in the sea level field are related to the impact of the oscillations with periods 352, 374, and 379 days, which are identified in a form of small but significant amplitude maxima in the Fourier series spectra of sea level, wind speed and atmospheric pressure. One more mechanism of the amplitude modulation of annual waves is assumed to be related to the changes in frequency of the natural baroclinic oscillations of the Baltic Sea due to the interannual variations of its stratification.

Morskoj gidrofizičeskij žurnal. 2024;40(2):231-254
pages 231-254 views

Experimental and field research

Influence of Coastal Upwelling on Chlorophyll a Distribution in the Coastal Zone of the Southeastern Baltic Sea in Summer Periods, 2000–2019

Kapustina M.V., Zimin A.V.

Abstract

Purpose. The study is purposed at obtaining the quantitative estimates of coastal upwelling influence on the distribution of chlorophyll a in the coastal zone of the southeastern Baltic Sea during the summer seasons in 2000–2019.

Methods and Results. Based on the data both on frequency and duration of upwelling events for June – August 2000–2019 and the chlorophyll a concentrations derived from multi-sensor satellite observations in the coastal zone of the southeastern Baltic Sea, the long-term and monthly average values of the studied parameters are obtained, and the influence of upwelling events on the chlorophyll a concentration in the sea surface layer is assessed. The spatial variability of chlorophyll a in the coastal areas is found to be related to the influence of upwelling events. On the average, the chlorophyll a concentration decreases by more than 1 mg/m3 after an upwelling of any duration and in all summer months. The concentration drop is most significant after the upwelling events lasting more than 6 days.

Conclusions. The reduced chlorophyll a concentrations (as compared to the pre-upwelling values) are observed in course of a week after a coastal upwelling event in the southeastern Baltic Sea.

Morskoj gidrofizičeskij žurnal. 2024;40(2):255-270
pages 255-270 views

Temporal Variability of the Beryllium-7 (7Ве) Scavenging Ratio in the Sevastopol Region

Kremenchutskii D.A.

Abstract

Purpose. The study is purposed at identifying the features of temporal variability in the relation of 7Ве scavenging ratio at the seasonal and interannual time intervals depending on the amount and frequency of precipitation.

Methods and Results. The scavenging ratio was assessed based on the field data on the 7Ве concentration values in the atmosphere and precipitation in 2012–2020. Application of ANOVA made it possible to examine the relationship between the scavenging ratio and the variations in 7Be concentrations in the atmosphere and precipitation. The relationship between the scavenging ratio and the precipitation amount and frequency was studied by the correlation method of analysis.

Conclusions. The scavenging ratio values averaged over a season and a year varied within the ranges 423–1286 and 508–919, respectively. The geometric mean value of the scavenging ratio was . At the 95% confidence level, the intrayear variability of the scavenging ratio is absent. The variability in average annual values of the scavenging ratio demonstrates a decrease in 2013 (as compared to 2012) from 664 to 508, an increase to 919 in 2016 and again a decrease to 516 in 2020. The ANOVA results indicate that variability of the scavenging ratio values averaged over a season and a year is due to the variation in 7Ве concentration in precipitation by 90 and 74%, respectively. In its turn, the long-term seasonal variability of 7Ве concentration in precipitation is conditioned by the variations in 7Ве concentration in the atmosphere (r = 0.64) and in the precipitation amount (r = –0.50). The relationship between the annual values of 7Ве concentration in precipitation, on the one hand, and its concentration in the atmosphere or the precipitation parameters, on the other, has not been revealed. The results of the correlation analysis indicate that the variability of precipitation parameters (amount and frequency) produces no a statistically significant effect (at the 95 % confidence level) upon the fluctuations of the scavenging ratio values at the seasonal and annual time intervals.

Morskoj gidrofizičeskij žurnal. 2024;40(2):271-283
pages 271-283 views

Vertical Distribution of Oxygen and Hydrogen Sulfide in the Deep Part of the Black Sea Based on the 2017–2019 Expedition Data

Kondratev S.I., Masevich A.V.

Abstract

Purpose. The purpose of the work is to analyze the features of vertical distribution of the dissolved oxygen and hydrogen sulfide in the deep part of the Black Sea in the modern period.

Methods and Results. The data obtained in 11 expeditions of Marine Hydrophysical Institute (MHI) RAS in the Black Sea within the economic zone of Russia in 2017–2019 were used. These surveys included more than 200 deep-sea stations at which by means of a cassette of 12 bathometers of the Sea-Bird 911 plus CTD Seabird-Electronics INC device the hydrochemical samples were taken at specific isopycnal surfaces; as a rule, it was σt = 16.30; 16.20; 16.10; 16.00; 15.95; 15.90; 15.80; 15.60; 15.40; 15.20; 15.00 and 14.60 kg/m3. Such a scheme made it possible to determine both the general vertical distribution of oxygen in the oxycline and the depth of hydrogen sulfide occurrence with the accuracy 0.05 kg/m3 in the conditional density scale.

Conclusions. In all the surveys, a decrease in oxygen content with depth (and occurrence of oxycline, respectively) began below the isopycnal surface σt = 14.5 kg/m3. The position of the sub-oxygen zone upper boundary defined by the isooxygen 10 µmol/l was not strictly isopycnal, but fell on the range of isopycns – σt = 15.7 – 15.85 kg/m3. However, it was not possible to identify a relation between the change in the position of the upper boundary and a certain hydrological season. For example, the deepest occurrence of the upper boundary below σt = 15.8 kg/m3 was observed both in November and December 2017, and August 2018. The lowering of isooxygen 10 µmol/l to σt = 15.9 kg/m3 in the Kerch shelf area is related to a more voluminous and colder intermediate layer over the shelf in December 2017. The position of hydrogen sulfide upper boundary determined by isosulfide 3 µmol/l was found raised almost to σt = 16.0 kg/m3 in only one of 11 surveys in April 2017. And in all other cases (including the one in August 2017, i. e. six months after it was raised) it was invariably within the range σt = 16.10–16.15 kg/m3. Over the past 25 years, the concentration of hydrogen sulfide at the depths 1750–2000 m has remained unchanged at the level 383 ± 2 µmol/l.

Morskoj gidrofizičeskij žurnal. 2024;40(2):284-297
pages 284-297 views

Underwater Ridge Impact on the Motion of Anticyclonic Eddies over a Sloping Bottom as a Result of the Topographic Beta-Effect: Laboratory Experiment

Zatsepin A.G., Elkin D.N.

Abstract

Purpose. The work is purposed at investigating the underwater ridge impact on the motion of anticyclonic eddies over a sloping bottom as a result of the topographic beta-effect in the laboratory conditions.

Methods and Results. The experiments are carried out in a cylindrical tank located on a rotating platform. A cone with its apex up is inscribed into the tank so that its base coincides with the cylinder lower base. The cone height is less than the base radius. Before the start of each experiment, the tank was filled with fresh or salt water of certain salinity. The fluid layer height exceeds that of the cone in the tank. The anticyclonic eddies were generated using a local constant source of a blue-colored fresh water flow. The source was located directly below the water layer surface at a distance equal to half the tank radius from its center. Having achieved the critical diameter, the generated eddies, due to the topographic beta-effect, drifted along the isobaths in the “western” direction (“north” is at the cone top in the tank center). The experiments were carried out over the cone with a smooth surface, and over the cone with a ridge on its side whose height was significantly smaller than that of the cone located at the path of the eddy drift. In the experimental runs with the ridge, the drift both of barotropic (fresh water in the tank) and baroclinic (salt water in the tank) eddies slowed down as compared to the eddy drift velocities in the absence of the ridge. After crossing the ridge, the orbital velocity of eddies also decreased significantly.

Conclusions. Field observations and numerical modeling of the Sevastopol anticyclonic eddy in the Black Sea moving over the continental slope along the isobaths in the southwestern direction showed that the eddy motion slowed down in the area of the underwater ridge formed by a local rise in the bottom relief between two canyons – the Danube and the Western Dnieper paleochannels. The results of the laboratory experiment have confirmed the data of field observations and numerical modeling on a slowdown of the Sevastopol eddy motion and a decrease in its orbital velocity while crossing the underwater ridge due to the topographic beta-effect.

Morskoj gidrofizičeskij žurnal. 2024;40(2):298-311
pages 298-311 views

Спутниковая гидрофизика

Evaluation of Sea Ice Drift in the Arctic Marginal Ice Zone based on Sentinel-1A/B Satellite Radar Measurements

Plotnikov E.V., Kozlov I.E., Zhuk E.V., Marchenko A.V.

Abstract

Purpose. The object of the work is to construct an automated system for calculating sea ice drift velocity fields using Sentinel-1A/B radar measurements based on the normalized maximum cross-correlation approach. The conditions and results of a numerical experiment aimed at evaluating the effectiveness of this technique for 63 pairs of radar images of the Fram Strait region for the summer-autumn periods in 2017 and 2018 are presented. Both the calculation algorithm and the qualitative and quantitative characteristics of the results are described in details. The effectiveness of the approach being applied to regular monitoring of ice drift is considered.

Methods and Results. The maximum cross-correlation (MCC) method was used for calculations. It is based on an automated finding of photographically similar fragments in the pairs of images, for which the sensing time interval is known. The Pearson correlation coefficient was applied as a proximity metric. As a result, 63 sea ice drift velocity fields were constructed for the Fram Strait region, each of which has a spatial scale of approximately several hundred thousand square kilometers. The method for filtering false correlations is proposed.

Conclusions. The approach applied in the study makes it possible to obtain automatically the sea ice drift velocity fields from the satellite data with high spatial resolution (40 m). The reconstructed velocity fields cover significant areas of the ocean surface. The method proposed for filtering false correlations permits to extract effectively the fragments with distortions resulting from the MCС algorithm limitations, from the calculation results.

Morskoj gidrofizičeskij žurnal. 2024;40(2):312-324
pages 312-324 views

Mathematical modeling of marine systems

Extreme Black Sea Storm in November, 2023

Dulov V.A., Yurovskaya M.V., Fomin V.V., Shokurov M.V., Yurovsky Y.Y., Barabanov V.S., Garmashov A.V.

Abstract

Purpose. The purpose of the study is to describe comprehensively the extreme storm in the Black Sea in November, 2023 in the terms of characteristics of the wind and wave fields, based on model calculations, satellite data and field measurements.

Methods and Results. The atmospheric fields are calculated using the WRF model, and the wave fields – by the SWAN model. The wind and wave fields, as well as their development during the storm are described in details. The phenomenon of wave shading by the Crimean Peninsula is studied. The data available for the storm period were involved in comparing the calculation results to the data of satellite altimeters, the CFOSAT SWIM wave scatterometer, and the synthetic aperture radars. The data of contact measurements carried out during the storm by the standard equipment installed at the oceanographic platform at the Black Sea hydrophysical subsatellite polygon of Marine Hydrophysical Institute of RAS (the coastal zone of the southern coast of Crimea) are presented. The wave characteristics near the oceanographic platform were calculated using the nested grid method.

Conclusions. It is found that during the storm in the Black Sea in November 2023, the maximum wave heights and the maximum wave periods exceeded 9 m and 13 s, respectively. The calculation results are confirmed by a large amount of satellite data. The results of calculating the characteristics of waves near the oceanographic platform are consistent with the in situ measurements performed from the oceanographic platform. Since the applied configuration of models and their settings permitted to retrieve the fields of wave physical characteristics with a high degree of reliability, they can be used for an authentic forecast of extreme storms in the Black Sea. The shading of waves by the Crimean Peninsula has led to the two and more times decrease in the heights of extreme waves in the coastal waters from the southern tip of the peninsula to Cape Chauda (35.8 оE).

Morskoj gidrofizičeskij žurnal. 2024;40(2):325-347
pages 325-347 views

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

8. Субъект персональных данных вправе запретить своему оборудованию прием этих данных или ограничить прием этих данных. При отказе от получения таких данных или при ограничении приема данных некоторые функции Сайта могут работать некорректно. Субъект персональных данных обязуется сам настроить свое оборудование таким способом, чтобы оно обеспечивало адекватный его желаниям режим работы и уровень защиты данных файлов «cookie», Оператор не предоставляет технологических и правовых консультаций на темы подобного характера.

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».