МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МОРСКИХ СИСТЕМ
Научная статья
УДК 551.465 БОГ 10.29039/0233-7584-2023-3-384-398
ЕБ№ РОЬСКБ
Особенности изменчивости плотности доступной потенциальной энергии и работы силы плавучести в верхнем 300-метровом слое Черного моря по результатам численного моделирования
С. Г. Демышев, О. А. Дымова Н. А. Миклашевская
Морской гидрофизический институт РАН, Севастополь, Россия н [email protected]
Поступила в редакцию 09.08.2022; одобрена после рецензирования 04.10.2022; принята к публикации 09.03.2023.
Аннотация
Цель. Основная цель исследования - выявление возможных физических механизмов изменчивости плотности доступной потенциальной энергии и потоков плавучести в верхнем деятельном слое Черного моря.
Методы и результаты. Исследование пространственного распределения плотности доступной потенциальной энергии и потоков плавучести проводилось на основе анализа термохалинных характеристик циркуляции Черного моря в 2011 и 2016 гг., полученных в результате численных экспериментов. Расчеты выполнены с помощью модели динамики Черного моря, разработанной в Морском гидрофизическом институте. В модели использовались батиметрия EMODNet и данные системы SKIRON о скорости ветра, потоках тепла, осадков, испарения и температуры поверхности моря. В результате численных экспериментов получены ежесуточные поля скорости течений, температуры и солености, по которым вычислялись плотность доступной потенциальной энергии и работа силы плавучести. Показано, что пространственно-временная изменчивость плотности доступной потенциальной энергии в Черном море формируется различными механизмами для верхнего 30-метрового слоя и слоя основного галоклина. Выявлена сезонная изменчивость работы силы плавучести.
Выводы. В верхнем слое изменчивость плотности доступной потенциальной энергии связана прежде всего с распространением распресненных речных вод, тогда как в слое основного галоклина (75-150 м) структура поля определяется мезомасштабной динамикой. В первом случае повышенные значения плотности доступной потенциальной энергии в течение года наблюдаются в районе северо-западного шельфа и на периферии бассейна; в центральной части моря распределение плотности доступной потенциальной энергии определяется атмосферными условиями. В слое ниже горизонта 75 м максимальные значения плотности доступной потенциальной энергии соответствуют антициклоническим круговоротам. Вследствие интенсивного перемешивания вод верхнего деятельного слоя в холодный период года работа силы плавучести определяется вертикальной скоростью. В весенне-летний период наблюдается двухслойная структура поля, определяемая знаком аномалий плотности. Толщина верхнего слоя составляет 20-30 м и соответствует глубине сезонного термоклина. В основном галоклине наибольшие абсолютные значения работы силы плавучести наблюдаются в зонах интенсивных мезомасштабных антициклонов.
Ключевые слова: Черное море, численное моделирование, вихри, доступная потенциальная энергия, работа силы плавучести, скорости течений, плотность, сезонная изменчивость
© Демышев С. Г., Дымова О. А., Миклашевская Н. А., 2023
384 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 39 №3 2023
Благодарности: проведение численных экспериментов и анализ термохалинных характеристик выполнены О. А. Дымовой при финансовой поддержке гранта РНФ № 22-77-10056. Анализ энергетических характеристик выполнен в рамках государственного задания по теме № FNNN-2021-0004 (шифр «Океанологические процессы»).
Для цитирования: Демышев С. Г., Дымова О. А., Миклашевская Н. А. Особенности изменчивости плотности доступной потенциальной энергии и работы силы плавучести в верхнем 300-метровом слое Черного моря по результатам численного моделирования II Морской гидрофизический журнал. 2023. Т. 39, № 3. С. 384-398. EDNPGLCKB. doi:10.29039/0233-7584-2023-3-384-398
Original article
Variability of Available Potential Energy Density and Buoyancy Work in the Upper 300-m Layer of the Black Sea Based on the Simulation Results
S. G. Demyshev, O. A. Dymova N. A. Miklashevskaya
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia H [email protected]
Abstract
Purpose. The study is aimed at identifying possible physical mechanisms for the variability of available potential energy density and buoyancy fluxes in the upper active layer of the Black Sea. Methods and Results. Spatial distribution of the available potential energy density and buoyancy fluxes was studied based on analyzing the thermohaline characteristics of the Black Sea circulation in 2011 and 2016 resulted from the numerical experiments performed using the Black Sea dynamics model developed in the Marine Hydrophysical Institute. The model included the EMODNet bathymetry and the SKIRON system data on wind velocity, heat fluxes, precipitation, evaporation, and sea surface temperature. The numerical experiments provided the daily fields of current velocities, temperature and salinity based on which the density of available potential energy and the buoyancy work were calculated. It is shown that the spatial-temporal variability of the available potential energy density in the Black Sea was formed by the mechanisms different for the upper 30-m layer and for the main halocline layer. The buoyancy work was revealed to be of seasonal variability.
Conclusions. In the upper layer, the variability of available potential energy density is related primarily to the propagation of freshened river waters, whereas in the main halocline layer (75-150 m), the field structure is conditioned by mesoscale dynamics. In the first case, the increased values of the available potential energy density are observed during a year on the northwestern shelf and on the basin periphery; in the central part of the sea, the distribution of available potential energy density is determined by the atmospheric conditions. In the layer below 75 m, the maximum values of the available potential energy density correspond to the anticyclonic eddies. In consequence of the intensive water mixing in the upper active layer during the cold period of a year, the buoyancy work is conditioned by vertical velocity. In a spring-summer period, a two-layer structure of the field is observed which is governed by the sign of density anomalies. The upper layer thickness constitutes 20-30 m and corresponds to the depth of seasonal thermocline. In the main halocline, the highest absolute values of the buoyancy work are observed in the zones of intense mesoscale anticyclones.
Keywords: Black Sea, numerical modeling, eddies, available potential energy, buoyancy work, current velocities, density, seasonal variability
Acknowledgments: Numerical experiments and analysis of thermohaline characteristics were performed by O. A. Dymova with the financial support of the RSF grant No. 22-77-10056. The analysis of energy characteristics was carried out within the framework of the state task on theme No. FNNN-2021-0004 (code "Oceanological processes").
For citation: Demyshev, S.G., Dymova, O.A. and Miklashevskaya, N.A., 2023. Variability of Available Potential Energy Density and Buoyancy Work in the Upper 300-m Layer of the Black Sea Based on the Simulation Results. Physical Oceanography, 30(3), pp. 355-368. doi: 10.29039/0233-75842023-3-355-368
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 39 №3 2023 385
Введение
Доступная потенциальная энергия (ДПЭ) - одно из ключевых понятий при анализе трансформации механической энергии в атмосфере и океане. По определению ДПЭ представляет собой часть полной потенциальной энергии, которая может переходить в кинетическую энергию (КЭ). Основные физические процессы, формирующие бюджет ДПЭ, - это диссипация энергии и работа, совершаемая силой плавучести при перемещении водных масс в вертикальном направлении. Потоки плавучести на границах области (свободная поверхность, устья рек и проливы) могут быть как источниками, так и стоками ДПЭ. Согласно классическим представлениям [1-3], высвобождение ДПЭ в океане и ее превращение в КЭ связано с вихревой изменчивостью. Современные исследователи подтверждают эту концепцию на основе результатов реалистичных численных экспериментов и данных наблюдений [4-6]. Помимо оценок ДПЭ как основного источника энергии мезомасштабных вихрей в океане, диссипация ДПЭ играет важную роль в турбулентном перемешивании [7, 8]. Вклад ДПЭ океана в глобальное взаимодействие океана и атмосферы обсуждается в [9]. Из вышесказанного следует, что анализ механизмов изменчивости ДПЭ представляет актуальную океанологическую задачу.
Для Черного моря оценки величины запаса ДПЭ по данным наблюдений представлены, например, в монографии 1 и работах [10, 11]. По сведениям авторов, средняя плотность ДПЭ в верхнем деятельном слое составляет около 20-30 Дж/м3, в [10, 11] также указано, что наибольшая плотность ДПЭ наблюдается в слое постоянного галоклина (75-150 м). Работы по численному анализу бюджета ДПЭ в Черном море немногочисленны. Первые количественные оценки вклада работы силы плавучести можно найти в [12]. В [13] впервые для Черного моря предложена методика численного анализа всех компонентов бюджета энергии и представлены результаты для климатической циркуляции. Энергетические переходы между ДПЭ и КЭ в идеализированной (циркуляция возбуждается стационарным ветром) двухслойной модели исследованы в [14].
По результатам численных экспериментов с реалистичными атмосферными условиями [15-17] было показано, что среднегодовой интегральный поток энергии, формируемый работой силы плавучести, направлен из средней ДПЭ в КЭ среднего течения. Следовательно, средняя циркуляция в Черном море поддерживается как ветровой накачкой, так и бароклинной неустойчивостью среднего течения. В случае ослабленного ветрового воздействия энергетические вклады потоков плавучести могут быть соизмеримы с вкладом работы силы ветра (см. рис. 6, в [17, с. 268]).
Анализ сезонной изменчивости энергетических характеристик [18] показал, что величина плотности ДПЭ существенным образом зависит от сезона года: ее наибольшие значения наблюдаются летом и расположены в слое от поверхности до верхней границы сезонного термоклина (верхний 20-метровый слой). Вертикальное распределение работы силы плавучести носит более сложный и нерегулярный характер. В верхнем 20-метровом слое работа силы
1 Изменчивость гидрофизических полей Черного моря / А. С. Блатов [и др.]. Л. : Гидроме-
теоиздат, 1984. 240 с.
386 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 39 №3 2023
плавучести положительна, что указывает на трансформацию энергии из ДПЭ в КЭ, в слое 20-40 м - отрицательна, т. е. ДПЭ увеличивается за счет КЭ. Ниже преобладают положительные значения.
Цель настоящей работы - детально исследовать пространственное распределение плотности ДПЭ и потоков плавучести на основе анализа термохалин-ных характеристик циркуляции и выявить возможные физические механизмы изменчивости ДПЭ и потоков плавучести в верхнем деятельном слое Черного моря.
Материалы и методы
Данная работа является продолжением исследования, начатого в [17, 18], поэтому для достижения поставленной цели использованы результаты двух численных экспериментов по моделированию циркуляции Черного моря в 2011 и 2016 гг. Эти годы выбраны исходя из тех соображений, что динамика вод отражала «бассейновый» (2011) и «вихревой» (2016) режимы циркуляции [19]. В первом случае в структуре поля течений доминировало Основное Черноморское течение (ОЧТ), во втором - в центральной части преобладали мезо-масштабные вихри.
Расчеты выполнены с помощью модели Морского гидрофизического института (МГИ) [20]. Модель построена на основе полной системы уравнений термогидродинамики океана в приближении Буссинеска, гидростатики и несжимаемости морской воды. Плотность нелинейно зависит от температуры и солености. Возвышение уровня рассчитывается с использованием линеаризованного кинематического условия на свободной поверхности. Потоки тепла, осадки, испарение и касательные напряжения трения ветра задаются в качестве граничных полей на поверхности моря. На дне ставятся условия прилипания и отсутствия потоков тепла и соли, трение о дно не учитывается. На твердых боковых участках границы выполняются условия скольжения, для температуры и солености задано равенство нулю нормальных производных. На жидких участках границы учитывается сток рек и водообмен через проливы. Вертикальное турбулентное перемешивание параметризовано с помощью модели турбулентного замыкания Меллора - Ямады. Горизонтальная турбулентная вязкость и диффузия представлены бигармоническими операторами с коэффициентом 1016 см4-с-1.
Модель МГИ реализована на сетке С с равномерным шагом по горизонтальным координатам 1,6 км, по вертикали заданы 27 неравномерных z-гори-зонтов со сгущением в верхнем деятельном слое. Батиметрия бассейна построена по данным EMODnet (URL: http://portal.emodnet-bathymetry.eu) с разрешением (1/8)'. В качестве атмосферного воздействия использованы 6-часовые данные о скорости ветра на высоте 10 м, испарении, температуре поверхности моря, потоках тепла и осадков, предоставляемые системой SKIRON с пространственным разрешением 0,1° [21]. Математическая формулировка модели МГИ, параметры численных экспериментов и валидация результатов расчетов подробно представлены в [17]. В целом структура модельной циркуляции соответствует общепринятым представлениям для Черного моря [22]: в осенне-
зимний сезон наблюдается усиление скорости течений и интенсивное вертикальное перемешивание; в теплое время года - ослабление бассейновой динамики, развитие мезомасштабных вихрей различного знака, нагрев и распрес-нение верхних слоев моря, формирование холодного промежуточного слоя.
В результате численных экспериментов получены ежесуточные поля скорости течений, температуры и солености, по которым вычислялись ДПЭ и работа силы плавучести. Основные гидрофизические поля для оценки исследуемых энергетических характеристик - это плотность морской воды и вертикальная компонента скорости течения. В модели МГИ локальная плотность р рассчитывается по формуле Мамаева [23] как нелинейная зависимость от температуры и солености:
где ро - плотность пресной воды; Т - температура; - соленость; а - коэффициенты термического расширения и соленостного сжатия [23].
Вертикальная скорость вычисляется из уравнения неразрывности
где и, V - горизонтальные компоненты вектора скорости; ^ - вертикальная компонента. Модель МГИ построена в правой декартовой системе координат: ось х направлена на восток, ось у - на север, ось г - вниз от поверхности ко дну. Положительное или отрицательное значение вертикальной скорости определяет опускание или подъем вод соответственно. Точность выполнения уравнения неразрывности составляет 10~12 см/с.
Для океана не существует точной формулы для расчета ДПЭ, поэтому мы использовали формулу, предложенную в [24]:
где Рт - средняя ДПЭ; V- объем слоя; g - ускорение свободного падения; р -локальная плотность; р* - аномалия локальной плотности. Черта над символом обозначает осреднение по времени; < > - осреднение по площади слоя с учетом маски «суша - море». Параметр (р) рассчитан как осредненная по соответствующему модельному горизонту локальная плотность и является константой для каждого слоя. В формуле (1), в отличие от формулы, используемой в [24], отсутствует знак «минус» перед интегралом. Это связано с выбором направления оси г (вниз в модели МГИ), что приводит к положительным значениям вертикального градиента плотности, совпадающего с направлением ускорения силы тяжести. Работа силы плавучести ЕЖ вычисляется по формуле
р = р0(1 + а[Т + а^Б + а^Т2 + а ^ТБ),
,Т5
(2)
Формула (2) строго выводится из уравнений скорости изменения КЭ и ДПЭ. Как показано в [3, 12, 13, 16], работа силы плавучести входит в виде отдельного слагаемого в оба уравнения, отличается только ее знак в соответствующем уравнении. Следовательно, величина и знак работы силы плавучести указывают на интенсивность и направление преобразования энергии между КЭ и ДПЭ. В модели МГИ положительное значение ЕЖ соответствует преобразованию из ДПЭ в КЭ. Интервал осреднения по времени в формулах (1) и (2) выбран равным одному месяцу, что позволяет учесть сезонные вариации потоков энергии.
Результаты
Для анализа были рассчитаны среднемесячные значения плотности ДПЭ и плотности работы силы плавучести для каждого модельного горизонта и проведено сопоставление полей Т, Б, м>, р*,Рт и ЕЖ на горизонтах 5, 30, 50 и 100 м, а также вертикальных разрезов всех указанных характеристик на зональном разрезе вдоль 43° с. ш. Далее в тексте для краткости мы опускаем термин «плотность» для обозначения ДПЭ и работы силы плавучести в единице объема.
28 30 32 34 36 3 8 40 °в.д. 28 30 32 34 36 38 40 °В,Д-
Р и с. 1. Средняя за июль плотность ДПЭ (a, c) и соленость (b, d) на горизонте 5 м в 2011 (a, b) и 2016 гг. (c, d)
F i g. 1. Monthly average APE density (a, c) and salinity (b, d) at the 5 m horizon in July 2011 (a, b) and July 2016 (c, d)
Для обоих рассматриваемых периодов пространственные структуры полей ДПЭ, аномалии плотности и солености на всех горизонтах качественно сходны. Зоны максимальных значений ДПЭ совпадают с локализацией пониженных значений солености. В верхнем 20-метровом слое пространственно-временная изменчивость ДПЭ отражает влияние речного стока и атмосферных
условий. В течение всего года наибольшие значения ДПЭ наблюдаются в районе устьев рек Днепр и Дунай, в центральной части моря ДПЭ изменяется мало. На рис. 1 показаны средние за июль поля ДПЭ и солености на горизонте 5мв2011 и 2016 гг. Видно, что наибольшие значения ДПЭ в северо-западной части моря привязаны к зонам распреснения прибрежных вод (рис. 1, Ь, с1), формирующимся за счет речного стока.
С увеличением глубины ДПЭ уменьшается на 1-2 порядка. На горизонтах ниже 50 м пространственная структура ДПЭ определяется вихревой активностью. На рис. 2 представлены распределения ДПЭ и солености на горизонте 100 м в 2011 и 2016 гг. Области повышенных значений ДПЭ (рис. 2, а, с) совпадают с зонами больших по модулю аномалий плотности, которые формируются в результате изменения солености (рис. 2, Ь, при подъеме и опускании вод в циклонических и антициклонических вихрях соответственно. Отметим, что запас ДПЭ в антициклонах (области пониженных значений солености на рис. 2, Ь, с) выше, чем в циклонах, примерно в 2-3 раза.
Сезонная изменчивость атмосферных потоков мало влияет на ДПЭ в основном галоклине. Для исследуемых периодов максимум ДПЭ на горизонте 100 м выявлен в июле 2011 г. в зоне Батумского антициклона, а в 2016 г. -в марте и октябре в зоне Севастопольского антициклона.
28 30 32 54 И 3» 40 °В.Д. 2S 30 32 34 36 38 40 °В.Д.
28 30 32 34 36 3R 40 °в.д. 2а 30 32 34 36 38 40 "пд
Рис. 2. Средняя за июнь плотность ДПЭ (a, c) и соленость (Ъ, d) на горизонте 100 м в 2011 (a, Ъ) и 2016 гг. (c, d)
F i g. 2. Monthly average APE density (a, c) and salinity (b, d) at the 100 m horizon in June 2011 (a, b) and June 2016 (c, d)
По результатам анализа пространственно-временную изменчивость работы силы плавучести можно условно разделить на два типа: в верхнем 30-метровом слое и в слое основного галоклина (50-150 м согласно [22]). Рассмотрим подробнее причины такого распределения по глубине. В верхнем слое для обоих периодов в холодное время года аномалии плотности положительны,
в теплое - отрицательны. Поэтому пространственное распределение положительных и отрицательных областей величины ЕЖ определяется только сменой знака вертикальной скорости. Как видно из рис. 3, а, Ъ, осенью и зимой вертикальная скорость (а следовательно, и поле работы силы плавучести) характеризуется сильной перемежаемостью зон положительных и отрицательных значений с максимальными абсолютными величинами на периферии бассейна. В июне - сентябре (рис. 3, с, ф в восточной части моря преобладают области подъема вод ^ < 0), в западной части - опускания ^ > 0). Таким образом, в летний период в верхнем слое западной части моря энергия преобразуется из КЭ в ДПЭ (ЕЖ< 0), а в восточной части наоборот - из ДПЭ в КЭ (ЕЖ>0) независимо от режима циркуляции.
F i g. 3. Vertical velocity averaged over the upper 30 m layer in February 2011 (a), February 2016 (b), August 2011 (c), and August 2016 (d)
Ниже горизонта 30 м не наблюдается горизонтальной однородности аномалий плотности по пространству. Это связано, прежде всего, с бассейновой циклонической схемой циркуляции вод Черного моря (рис. 4, а), которая обусловливает опускание изопикн на периферии бассейна и их подъем в центральной части, формируя таким образом положительную аномалию плотности в центре (рис. 4, с). При слабой циклонической циркуляции (рис. 4, b) опускание изопикн у материкового склона выражено слабо.
Однако в обоих экспериментах развивающиеся вблизи склона интенсивные мезомасштабные антициклоны (рис. 4, a, b) способствуют дополнительному поступлению распресненных вод с вышележащих горизонтов, что четко проявляется в виде областей с минимальными аномалиями плотности (рис. 4, с, d)rn максимальными значениями вертикальной скорости (рис. 4, e, f).
Рис. 4. Средние за июнь горизонтальная (a, Ъ) и вертикальная (в, f) скорость, аномалия плотности (c, d) на горизонте 100 м в 2011 (a, c, в) и 2016 гг. (Ъ, d, f)
F i g. 4. Monthly average horizontal (a, b) and vertical (в, f) velocities, and density anomaly (c, d) at the 100 m horizon in June 2011 (a, c, в) and June 2016 (b, d,f)
При этом на периферии антициклонов происходит подъем вод, что можно продемонстрировать на примере Севастопольского антициклона. На рис. 5 представлено распределение горизонтальной (рис. 5, a) и вертикальной (рис. 5, b) скорости в Севастопольском антициклоне в середине июня 2011 г. Видно, что в ядре вихря (примерно область с координатами 44,23° с. ш., 31,5° в. д.) происходит опускание вод (w > 0) и прогиб изохалин. На периферии вихря вода поднимается (градация синего цвета на рис. 5, b), причем скорость подъема выше на правой границе, где наблюдаются повышенные значения орбитальной скорости. В циклонических вихрях ситуация обратная: в центре вода поднимается, на периферии - опускается. Как видно из рис. 5, с, ядру циклонического вихря в юго-восточной части моря в июне 2016 г. соответствует подъем изохалин (рис. 5, d).
Рис. 5. Горизонтальная скорость на горизонте 100 м (а, с), вертикальная скорость на зональных разрезах вдоль 44,23° с. ш. (Ъ) и вдоль 42° с. ш. (d) на 15.06.2011 (a, Ъ) и 18.06.2016 (с, d). Черные линии - изохалины (интервал 0,5 %о)
F i g. 5. Horizontal velocity at the 100 m horizon (a, c), vertical velocity on the zonal sections along 44.23° N (b) and 42° N (d) on June 15, 2011 (a, b) and June 18, 2016 (c, d). Black lines denote isohalines (shading is 0.5 %o)
Положительные и отрицательные аномалии плотности и чередование зон подъема - опускания вод приводят к сложной структуре поля работы силы плавучести ниже горизонта 30 м. На рис. 6 показаны средние за июнь карты BW для двух расчетов. В силу вида формулы (2) нулевые изолинии BW, р* и w совпадают. Вдоль изолинии, разделяющей положительные и отрицательные аномалии плотности (рис. 4, с, d), наблюдается смена знака работы силы плавучести (рис. 6, а, Ъ) для обоих экспериментов. Эта граница пространственно соответствует стрежню ОЧТ в 2011 г. (рис. 4, а) и областям мезомасштабных вихревых образований в 2016 г. (рис. 4, Ъ). Минимальные значения BWнаблюдаются в ядрах антициклонических вихрей с максимальными значениями орбитальной и вертикальной скорости, максимумы BW - на периферии интенсивных антициклонов. В центральной глубоководной части моря диапазон изменчивости BW в несколько раз меньше, чем в вихрях, а смена знака вертикальной скорости может быть обусловлена быстрыми динамическими и тер-мохалинными процессами, трудно идентифицируемыми на среднемесячных картах исследуемых параметров. Таким образом, наиболее энергоактивные зоны с точки зрения преобразования энергии между кинетической и доступной потенциальной - это мезомасштабные антициклоны.
Рис. 6. Средняя за нюнь плотность работы силы плавучести на горизонте 100 м в 2011 (а) и 2016 гг. (b)
F i g. 6. Monthly average buoyancy work per unit volume at the 100 m horizon in June 2011 (a) and June 2016 (b)
4 4 -i
0,2
0,16
ro 2
"Й
-0,12 4
m " О
-0,08 ¿E
0,04
0,0
""" 28 32 36 40 °в.д. JW 28 32 . 36 40 °в.д. a D
Рис. 7. Средние плотность ДПЭ (цвет) и соленость (черные линии) на зональном разрезе вдоль 43° с. ш. в апреле 2011 г. (а) и ноябре 2016 г. (b)
F i g. 7. Monthly average APE density (color) and salinity (black lines) on the zonal section along 43° N in April 2011 (a) and November 2016 (b)
Анализ вертикальных сечений энергетических и гидрофизических полей в 2011 и 2016 гг. показал, что максимальные значения ДПЭ локализованы в прибрежной северо-западной части моря, где распространяются речные воды (рис. 1, а, с). По вертикали влияние речного стока прослеживается до глубин 25-30 м. С глубиной плотность ДПЭ уменьшается, а в пространственном распределении наблюдаются области повышенных значений (по сравнению с окружающими водами), которые можно идентифицировать как проявление вихревой активности в слое галоклина. На рис. 7 показаны зональные разрезы полей ДПЭ и солености. Видно, что в слое 75-125 м в центральной части моря повышенные значения ДПЭ соответствуют зонам изменения наклонов изоха-лин. Как показано на рис. 5 прогиб изохалин соответствует опусканию вод в ядре антициклона, подъем - движению к поверхности в ядре циклона. Таким
образом, по аналогии можно заключить, что наибольшие значения ДПЭ в области опускания изохалин (указано красной стрелкой на рис. 7, Ъ) находятся в антициклоне, а наибольшие значения ДПЭ в области подъема изохалин (синие стрелки на рис. 7, а) в циклоне. Полученный результат подтверждается тем фактом, что антициклоны в Черном море более интенсивны [25]. Этот вывод согласуется с рис. 2, с, где повышенные значения ДПЭ на горизонте 100 м пространственно соответствуют зонам с максимальными и минимальными аномалиями плотности (рис. 4, в циклонах и антициклонах (рис. 4, Ъ) соответственно. Следует отметить, что временная изменчивость ДПЭ в верхнем слое прибрежной зоны определяется сезонным увеличением речного стока, тогда как в слое галоклина связи со сменой сезонов не выявлено. Максимумы ДПЭ ниже горизонта 30 м обнаружены в апреле 2011 г. (рис. 7, а) и в ноябре 2016 г. (рис. 7, Ъ).
d е /
Рис. 8. Среднемесячные работа силы плавучести в единице объема (a, d), аномалия плотности (b, e) и вертикальная скорость (с, f на зональном разрезе вдоль 43° с. ш. в марте (a - с) и августе (d -f 2011 г.
F i g. 8. Monthly average buoyancy work per unit volume (a, d), density anomaly (b, e) and vertical velocity (с,f on the zonal section along 43° N in March (a - c) and August (d -f), 2011
Согласно формуле (2), знак работы силы плавучести зависит от произведения вертикальной скорости и аномалии плотности. На рис. 8 представлены вертикальные разрезы BW, р* и w в конце зимнего и летнего гидрологических сезонов. В холодный сезон, когда аномалия плотности достаточно однородна по вертикали (рис. 8, b), направление преобразования энергии между ДПЭ и КЭ определяется зонами подъема и опускания вод. Там, где вода поднимается (w < 0, градация голубого цвета на рис. 8, с), работа силы плавучести отрицательна, следовательно, ДПЭ увеличивается за счет КЭ. В теплый сезон,
когда аномалия плотности меняет знак с глубиной (рис. 8, в), работа силы плавучести (рис. 8, й) также меняет знак в зависимости от знаков р* и w. На рис. 8, й наблюдается некая условная граница примерно на глубине 30 м, пересекая которую работа силы плавучести меняет знак. Слой 20-30 м для данных экспериментов соответствует глубине залегания сезонного термоклина (рис. 10 в работе [18, с. 15]). Таким образом в теплый период года, когда вследствие формирования сезонного термоклина меняется структура поля аномалии плотности, работа силы плавучести может менять знак с глубиной. При этом в распределении вертикальной скорости не выявлено связи с глубиной залегания сезонного термоклина или глубиной постоянного галоклина. Качественно эта ситуация повторяется и в 2016 г. с тем отличием, что наблюдается большее количество областей ЕЖ и ^ различного знака вследствие интенсивной вихревой изменчивости циркуляции в2016г.по сравнению с 2011 г.
Выводы
Согласно результатам проведенных экспериментов, пространственно-временная изменчивость плотности ДПЭ в Черном море формируется различными механизмами для верхнего 30-метрового слоя и слоя основного галоклина. В верхнем слое изменчивость ДПЭ связана прежде всего с распространением распресненных речных вод. Так, в течение года максимальные значения плотности ДПЭ наблюдаются в районе северо-западного шельфа, повышенные значения - на периферии бассейна. Верхний слой в центральной части моря подвержен слабой пространственной и временной изменчивости, и плотность ДПЭ здесь определяется атмосферными условиями.
В слое основного галоклина 75-150 м структура поля плотности ДПЭ определяется мезомасштабной динамикой. Максимальные значения плотности ДПЭ соответствуют зонам отрицательных аномалий плотности, которые формируются в антициклонических круговоротах. Анализ горизонтального пространственного распределения плотности ДПЭ показал, что наиболее энергоактивные зоны Черного моря - это области Севастопольского и Батумского антициклонов, причем ДПЭ тем выше, чем выше орбитальная скорость вихря. Отмечено, что плотность ДПЭ в циклонах в 2-3 раза меньше, чем в антициклонах. Такая структура соответствует тому факту, что антициклоны в Черном море более интенсивны и, следовательно, наклоны изопикн в них более крутые, что приводит к росту аномалий плотности.
Изменчивость работы силы плавучести, которая определяет скорость и направление трансформации энергии между ДПЭ и КЭ, проявляет сезонный характер. В холодный период года вследствие интенсивного перемешивания формируется положительная аномалия плотности во всем верхнем деятельном слое Черного моря. Поэтому работа силы плавучести в большей степени определяется знаком вертикальной скорости. В центральной части, где происходит подъем вод, поток энергии направлен из КЭ в ДПЭ, на периферии бассейна -наоборот. Весной и летом в поле работы силы плавучести наблюдается двухслойная структура, определяемая зонами положительных и отрицательных аномалий плотности. Толщина верхнего слоя составляет 20-30 м и, согласно
нашей работе, опубликованной в журнале Water в 2022 г., соответствует глубине сезонного термоклина.
Анализ структуры поля работы силы плавучести в основном галоклине показал, что наиболее интенсивное преобразование энергии между ДПЭ и КЭ происходит в зонах мезомасштабных антициклонических вихрей с максимальными орбитальными скоростями. Причем знак работы силы плавучести определяется преимущественно знаком вертикальной скорости, которая разнона-правлена в ядре и на периферии вихря.
СПИСОК ЛИТЕРАТУРЫ
1. Каменкович В.М., Кошляков М. Н.,Монин А. С. Синоптические вихри в океане. Л. : Гид-рометеоиздат, 1982. 264 с.
2. ГиллА. Динамика атмосферы и океана. В2т. Пер. с англ. М. : Мир, 1986. Т. 2. 415 с.
3. Holland W. R. Energetics of baroclinic oceans // Numerical models of ocean circulation : proceedings of a symposium held at Durham, New Hampshire, October 17-20, 1972. Washington : National Academy Press, 1975. P. 168-177.
4. The global mesoscale eddy available potential energy field in models and observations // C. Luecke [et al.] // Journal of Geophysical Research: Oceans. 2017. Vol. 122, iss. 11. P. 91269143. https://doi.org/10.1002/2017JC013136
5. Li Q., Zhou L., Xie L. Seasonal and interannual variability of EAPE inthe South China Sea derived from ECCO2 data from 1997 to 2019 // Water. 2021. Vol. 13, iss. 7. 926. https://doi.org/10.3390/w13070926
6. Травкин В. С., Белоненко Т. В. Исследование вихревой изменчивости в Лофотенской котловине на основе анализа доступной потенциальной и кинетической энергии II Морской гидрофизический журнал. 2021. Т. 37, № 3. С. 318-332. doi:10.22449/0233-7584-2021-3-318-332
7. Available potential energy and mixing in density-stratified fluids / K. B. Winters [etal.] // Journal of Fluid Mechanics. 1995. Vol. 289. P. 115-128. https: //doi.org/10.1017/S002211209500125X
8. Tailleux R. Irreversible compressible work and available potential energy dissipation in turbulent stratified fluid // Physica Scripta. 2013. Vol. 2013, no. T155. 014033. https://doi.org/10.1088/0031-8949/2013/T155/014033
9. BishopS. P., SmallR. J., BryanF. O. The global sink of available potential energy by mesoscale air-sea interaction // Journal of Advances in Modeling Earth Systems. 2020. Vol. 12, iss. 10. e2020MS002118. https://doi.org/10.1029/2020MS002118
10. Гидрометеорология и гидрохимия морей СССР. Т. IV. Черное море. Вып. 1. Гидрометеорологические условия / Иод ред. А. И. Симонова, Э. Н. Альтмана. СПб. : Гидрометеоиздат, 1991. 428 с.
11. СуворовА. М., Шокурова И. Г. Годовая и междесятилетняя изменчивость доступной потенциальной энергии в Черном море II Морской гидрофизический журнал. 2004. № 2. С. 29-41.
12. Stanev E. V. On the mechanisms of the Black Sea circulation // Earth-Science Reviews. 1990. Vol. 28, iss. 4. P. 285-319. https://doi.org/10.1016/0012-8252(90)90052-W
13. Демышев С. Г. Энергетика климатической циркуляции Черного моря. Ч. I. Дискретные уравнения скорости изменения кинетической и потенциальной энергий II Метеорология и гидрология. 2004. № 9. С. 65-80. EDN PGCNXF.
14. ПавлушинА. А., Шапиро Н. Б., Михайлова Э. Н. Энергетические переходы в двухслойной вихреразрешающей модели Черного моря II Морской гидрофизический журнал. 2019. Т. 35, № 3. С. 201-219. https://doi.org/10.22449/0233-7584-2019-3-201-219
15. Демышев С. Г., Дымова О. А. Численный анализ мезомасштабных особенностей циркуляции в прибрежной зоне Черного моря II Известия РАН. Физика атмосферы и океана. 2013. Т. 49, № 6. С. 655-663. EDNRFWNFH. https://doi.org/10.7868/S0002351513060035
16. Demyshev S. G., Dymova O. A. Numerical analysis of the Black Sea currents and mesoscale eddies in 2006 and 2011 // Ocean Dynamics. 2018. Vol. 68, iss. 10. P. 1335-1352. https: //doi.org/10.1007/s10236-018-1200-6
17. Demyshev S. G., Dymova O. A. Analysis of the annual mean energy cycle of the Black Sea circulation for the climatic, basin-scale and eddy regimes II Ocean Dynamics. 2022. Vol. 72, iss. 3-4. P. 259-278. https://doi.org/10.1007/s10236-022-01504-0
18. Demyshev S., Dymova O., Miklashevskaya N. Seasonal variability of the dynamics and energy transport in the Black Sea by simulation data // Water. 2022. Vol. 14, iss. 3. 338. https://doi.org/10.3390/w14030338
19. Stanev E. V., Staneva J. V. The impact of the baroclinic eddies and basin oscillations on the transitions between different quasi-stable states of the Black Sea circulation // Journal of Marine Systems. 2000. Vol. 24, iss. 1-2. P. 3-26. https://doi.org/10.1016/S0924-7963(99)00076-7
20. Демышев С. Г. Численная модель оперативного прогноза течений в Черном море II Известия Российской академии наук. Физика атмосферы и океана. 2012. Т. 48, № 1.С. 137149. EDN OOWHLL
21. The regional weather forecasting system SKIRON: an overview / G. Kallos [et al.] // Proceedings of the International Symposium on Regional Weather Prediction on Parallel Computer Environments (Athens, Greece, 15-17 October 1997). Athens, 1997. P. 109-122.
22. Иванов В. А., Белокопытов В. H. Океанография Черного моря. Севастополь, 2011. 212 c.
23. Мамаев О. И. T, S-анализ вод Мирового океана. Л. : Гидрометеоиздат, 1970. 364 с.
24. New estimates of the available potential energy in the world ocean / A. H. Oort [et al.] // Journal of Geophysical Research: Oceans. 1989. Vol. 94, iss. C3. P. 3187-3200. https://doi.org/10.1029/JC094iC03p03187
25. Thermohaline structure, transport and evolution of the Black Sea eddies from hydrological and satellite data / A. A. Kubryakov [et al.] // Progress in Oceanography. 2018. Vol. 167. P. 44-63. https://doi.org/10.1016/j.pocean.2018.07.007
Об авторах:
Демышев Сергей Германович, главный научный сотрудник, заведующий отделом теории волн, ФГБУН ФИЦ МТИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), доктор физико-математических наук, SPIN-код: 1848-2350, ORCID ID: 0000-0002-5405-2282, ResearcherID: C-1729-2016, Scopus Author ID: 57862712800, [email protected]
Дымова Ольга Алексеевна, старший научный сотрудник, отдел теории волн, ФГБУН ФИЦ МТИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), кандидат физико-математических наук, SPIN-код: 7565-1082, ORCID ID: 0000-0003-4036-2447, ResearcherID: P-9669-2015, [email protected]
Миклашевская Надежда Андриановна, младший научный сотрудник, отдел теории волн, ФГБУН ФИЦ МТИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), SPIN-код: 8476-2604, ORCID ID: 0000-0003-2619-343X, ResearcherID: P-2167-2017, [email protected]