CC BY
1
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МОРСКИХ СИСТЕМ
Научная статья
УДК 551.465 ЕБ№ НТШОА
Потоки энергии между средними течениями и мезомасштабными вихрями в восточной и западной частях Черного моря
А. А. Павлушин
Морской гидрофизический институт РАН, Севастополь, Россия н рау1изЫп@тЫ-газ. ги
Поступила в редакцию 31.12.2022; одобрена после рецензирования 22.06.2023; принята к публикации 15.11.2023.
Аннотация
Цель. Целью данной работы является анализ потоков энергии между течениями разных масштабов в восточной и западной частях Черного моря.
Методы и результаты. Для определения составляющих энергетического баланса использовались результаты численного расчета полей течений в Черном море, выполненного с помощью двухслойной вихреразрешающей модели, возбуждаемой ветром с циклонической завихренностью. Полная нестационарная крупномасштабная циркуляция моря условно представлялась в виде суммы средних течений и мезомасштабных вихрей. Были рассчитаны переходы между кинетической и доступной потенциальной энергией, а также потоки энергии между крупномасштабными течениями и мезомасштабными вихрями отдельно для восточной и западной части Черного моря. Кроме того, были определены адвективные потоки энергии через меридиональную границу между двумя частями моря.
Выводы. Основной вклад в адвективный перенос потенциальной энергии из восточной части Черного моря в западную вносят мезомасштабные вихри. Этот адвективный перенос энергии приводит к возникновению различий в направлении потоков энергии между течениями разных масштабов в восточной и западной частях моря. В восточной части энергия передается от средних течений к мезомасштабным вихрям, а в западной наблюдается обратный поток энергии. Показано, что обмен кинетической энергией между различными масштабами обеспечивается работой силы Кориолиса, которая для исходных течений равна нулю, но оказывается не равной нулю отдельно для средних течений и мезомасштабных вихрей.
Ключевые слова: Черное море, крупномасштабная циркуляция, средние течения, мезомасштабные вихри, энергетический баланс, поток энергии, Р-эффект, волны Россби
Благодарности: работа выполнена в рамках темы государственного задания ^NNN-2022-0003 «Развитие методов оперативной океанологии на основе междисциплинарных исследований процессов формирования и эволюции морской среды и математического моделирования с привлечением данных дистанционных и контактных измерений».
Для цитирования: Павлушин А. А. Потоки энергии между средними течениями и мезомасштабными вихрями в восточной и западной частях Черного моря // Морской гидрофизический журнал. 2024. Т. 40, № 1. С. 108-129. ЕБЫНТИБОА.
© Павлушин А. А., 2024
Original article
Energy Flows between the Mean Currents and the Mesoscale Eddies in the Eastern and Western Parts of the Black Sea
А. A. Pavlushin
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia H pavlushin@mhi-ras. ru
Abstract
Purpose. The study consists in analyzing the energy flows between the currents of different scales in the eastern and western parts of the Black Sea.
Methods and Results. The energy balance components were determined based on the results of numerical calculation of current fields in the Black Sea performed using a two-layer eddy-resolving model subjected to a forcing of a wind with cyclonic vorticity. A complete non-stationary large-scale sea circulation was represented conventionally as a sum of mean currents and mesoscale eddies. Conversions between the kinetic and available potential energy, as well as the energy flows between the mean currents and the mesoscale eddies were calculated separately for the eastern and western parts of the Black Sea. Besides, the advective energy flows across the meridional boundary between two parts of the sea were also calculated.
Conclusions. The main contribution to the potential energy advective transfer from the eastern part of the Black Sea to its western part is made by the mesoscale eddies. This advective energy transfer results in arising of differences in the energy flow directions in the currents of different scales in the eastern and western parts of the sea. In the Black Sea eastern part, energy is transferred from the mean currents to the mesoscale eddies, whereas in its western part, a reverse energy flow is observed. It is shown that the kinetic energy exchange between the movements of different scales is provided by the Coriolis force work, which for the initial currents is equal to zero, but turns out to be non-zero separately for the mean flows and the mesoscale eddies.
Keywords: Black Sea, large-scale circulation, mean currents, mesoscale eddies, energy balance, energy flow, P-effect, Rossby waves
Acknowledgements: The study was carried out within the framework of state assignment FNNN-2022-0003 "Development of operational oceanology methods based on interdisciplinary research of the processes of marine environment formation and evolution, and mathematical modeling using the data of remote and contact measurements".
For citation: Pavlushin, A.A., 2024. Energy Flows between the Mean Currents and the Mesoscale Eddies in the Eastern and Western Parts of the Black Sea. Physical Oceanography, 31(1), pp. 99-119.
Введение
Большую роль в формировании и изменчивости горизонтальной циркуляции вод в океанах и морях играет взаимодействие между динамическими процессами различных масштабов. Согласно классификации, предложенной А. С. Мониным \ для циркуляции Черного моря актуальной является изменчивость межгодового, сезонного и синоптического масштабов. Процессы межгодового и сезонного масштабов проявляются в формировании течений, ограниченных размерами бассейна. В Черном море это Основное Черноморское течение (ОЧТ), которое в некоторые годы распространяется по всему периметру глубоководной части моря, а в другие годы распадается на два или три
1 МонинА. С., Каменкович В. М, Корт В. Г. Изменчивость Мирового океана. Ленинград : Гидрометеоиздат, 1974. 262 с.
МОРЖОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 40 № 1 2024 109
суббассейновых циклонических круговорота - «очки Книповича» [1]. К процессам синоптического масштаба относятся меандры ОЧТ, антициклонические и циклонические вихри, образующиеся в ложбинах меандров, ринги -оторвавшиеся меандры, квазистационарные антициклоны (Севастопольский, Батумский), волны Россби с масштабами порядка бароклинного радиуса деформации [2, 3].
В научных статьях зарубежных авторов [4-6] при описании особенностей полей течений чаще применяется другая классификация, согласно которой вихревые образования, по определению А. С. Монина являющиеся синоптическими [7], отнесены к мезомасштабным объектам. Такие объекты называются мезомасштабными вихрями (mesoscale eddies), имеют размеры от десяти до нескольких сотен километров, время их жизни составляет от недели до нескольких месяцев [8]. Соответственно, явления большего масштаба, чем мезомас-штабные вихри, относятся к крупномасштабным, меньшего - к субмезомас-штабным. В данной работе субмезомасштабные динамические объекты не рассматриваются, поскольку пространственное разрешение двухслойной вихре-разрешающей модели [9], с помощью которой проводится расчет полей гидрофизических характеристик, не позволяет этого сделать. В дальнейшем планируется провести некоторую модернизацию модели и увеличить ее пространственное разрешение, чтобы исследовать влияние субмезомасштабных процессов на более крупномасштабные, так как этот вопрос представляет научный интерес [10-12].
В настоящей работе для удобства изложения используется разделение процессов на крупномасштабные и мезомасштабные, как принято в зарубежной литературе. В качестве крупномасштабных рассматриваются средние течения (СТ), полученные путем осреднения по времени мгновенных полей течений, а отклонения мгновенных течений от этого среднего состояния являются вихревыми образованиями, чаще называемыми мезомасштабными вихрями (МВ). В Черном море МВ образуются в основном вследствие гидродинамической неустойчивости крупномасштабных течений [13]. Они участвуют в передаче энергии в глубинные слои моря, влияют на приток ветровой энергии и на горизонтальное перемешивание вод [14].
Для исследования динамики Черного моря в работе применяется метод энергетического анализа, при котором вычисляются и оцениваются энергетические параметры течений: кинетическая энергия и потенциальная энергия, энергетические потоки (переходы) между различными видами энергии и масштабами движения, источники и стоки энергии. Знание энергетики исследуемых процессов необходимо для лучшего и правильного понимания их физики [15, 16].
Постановка задачи, описание модели
Предметом рассмотрения в настоящей работе являются взаимные переходы кинетической и потенциальной энергии, вызванные крупномасштабными и мезомасштабными динамическими процессами в Черном море. На данном этапе исследований ставилась цель - рассчитать и проанализировать средние за большой период времени переходы механической энергии между ее раз-
личными видами и между движениями различного масштаба отдельно для восточной и западной части Черного моря. Разделение моря на две половины осуществлялось по условной границе Б, которая совпадает с меридианом, проходящим через м. Сарыч (рис. 1, а).
Р и с. 1. Поля касательного напряжения ветра (a) и завихренности т (b) F i g. 1. Fields of tangential wind stress (a) and vorticity т (b)
В предыдущей работе [17], также посвященной анализу энергетики Черного моря, было установлено, что под влиянием ß-эффекта формируется средний адвективный перенос потенциальной энергии из восточной половины моря в западную, который приводит к пространственной неравномерности распределения энергетических характеристик по акватории моря. В работе анализировались осредненные по времени интегральные потоки и переходы энергии в западной и восточной частях моря, но их разделение по масштабам движения не проводилось.
В качестве исходных данных для расчета энергетических характеристик как в предыдущей, так и в настоящей работе использовались результаты численного моделирования черноморской циркуляции с помощью двухслойной вихреразрешающей модели [9], в которой учитываются реальный рельеф дна, трение на поверхности раздела слоев, придонное трение и горизонтальная турбулентная вязкость.
Несмотря на кажущуюся простоту, в используемой модели учитываются все основные факторы, определяющие формирование и изменчивость поля течений. При правильном выборе внешних и внутренних параметров модель вполне адекватно воспроизводит известные особенности крупномасштабной черноморской циркуляции - Основное Черноморское течение с меандрами и тесно связанными с ними мезомасштабными вихревыми образованиями, Ба-тумский и Севастопольский квазистационарные антициклоны, «очки Книпо-вича» и др. Течения в нижнем слое, полученные с помощью двухслойной модели [18], соответствуют данным наблюдений глубинных течений в Черном море [19, 20].
Если для возбуждения движения в модели использовать стационарный или периодически меняющийся ветер, то удается получить выход модельного решения на статистически равновесный режим, при котором средние значения всех рассчитываемых характеристик и параметров модели остаются постоянными.
Модель основана на системе уравнений мелкой воды для двухслойной жидкости и включает в себя уравнения движения и неразрывности для каждого слоя:
Ц)( + ЦUl )х + (УЛ )у - V = g\Сх + тх - + Aв/ (AAUl), V), + ЦV), + (VVl )у + /и, = ^И,^ + ту -Я* + Ав/ (ДД^ ),
(И ), + Ц )Я + (у )У= 0 (Ц2 )( + Ци2 )х + Уи2 )у - ^ = gИ2Сх + gИ (/¡1 )х + ДЬХ2 - Я£ + АвИ (ААи2 ), V ), ^ ^2 ), ^ ^2 X + Ц = gh2^ + gИ (¡1 )у + ^2 - + АвИ2 (ААУ2 ),
(к2 \+{и21 + (у2 Х= 0
где нижние индексы 1, 2 указывают номер слоя; нижние индексы I, х, у означают дифференцирование по соответствующей переменной; (и1, у1), (и2, у2) -горизонтальные компоненты скорости течений в верхнем и нижнем слое, направленные вдоль осейX (на восток) и У (на север); И1, И2 — толщины слоев; £ - уровень моря; Xх, ту — компоненты касательного напряжения ветра, действующего на морскую поверхность; (Ц = и1И1, У = v1h1), (Ц = и2И2, У2 = v2h2) —
компоненты потоков в слоях; (Я^, RfУl), (Я£2, Я[2) — компоненты силы трения на нижней границе верхнего слоя и верхней границе нижнего слоя соответственно (на жидкой границе между слоями Я^ =— 2, Я[1 =— Я[2); (Я^, Я^ ) — компоненты придонного трения в нижнем слое;/=/ + ву — параметр Кориолиса, / = 10-4 1/с, в = 2-10-8 1/(с^м); g = 9,81м/с2 — ускорение свободного падения; g = g(p2 - р1)/р2 = 0,032 м/с2 — редуцированное ускорение свободного падения, р1, р2 - плотность воды в верхнем и нижнем слое; Ав - коэффициент бигармонической горизонтальной турбулентной вязкости;
л 52 52 й _
Д = —- л---г - двумерный оператор Лапласа.
дх су
Для замыкания системы уравнений в модели используется уравнение неразрывности для полных потоков в приближении «твердой крышки» Ц + Уу = = 0, где и = Ц + и2, V = У1 + У2 — составляющие полных потоков. Это позволяет ввести интегральную функцию тока у, для которой и = -уу, V = ух.
Речной сток и водообмен через проливы в модели не учитываются, на боковых границах бассейна задаются условия прилипания и1 = 0, и2 = 0.
Поскольку для параметризации горизонтальной турбулентной вязкости в модели используется бигармонический оператор, возникает необходимость определения граничных условий для лапласиана скорости течений. По аналогии с условиями для и1, и2 задаем Ди1 = 0, Ди2 = 0, хотя физического смысла в этом нет.
Модель вследствие слабой диссипативности наряду с крупномасштабными течениями хорошо воспроизводит мезомасштабные вихри, что дает возможность, используя результаты проведенных экспериментов, выделить доли энергии и работы сил, связанные с СТ и МВ. Методика подобного анализа описана в [21] и позволяет рассчитать энергетические переходы при взаимодействии мезомасштабных вихрей и крупномасштабных течений.
Исходные данные для расчета энергетических характеристик представляют собой поля толщины верхнего слоя, уровня моря и скоростей течений в верхнем и нижнем слоях за 30 лет с дискретностью 1 сут. Данные были получены в численном эксперименте при следующих параметрах модели: размер горизонтальной ячейки Ах = Ау = 3000 м, шаг интегрирования по времени А^ = = 120 с, толщина верхнего слоя в состоянии покоя Но = 100 м, коэффициент горизонтальной турбулентной вязкости Ав = 2,0-108 м4/с, коэффициент придонного трения г0 = 0,002, коэффициент трения между слоями гь = 2,0-Ю-6 м/с.
Ветровое воздействие задавалось стационарным полем касательного напряжения ветра т (х, у) с переменной по площади завихренностью (рис. 1, а, Ь), при этом средняя по акватории моря завихренность была циклонической и составляла 0,3•Ю-7 Н/м3. Выбор стационарного поля т объясняется тем, что в этом случае численная модель быстрее выходит на статистически равновесный режим.
На рис. 2 приводятся полученные в результате расчетов характерные для статистически равновесного режима поля уровня моря ^ (х, у) и потенциала Монтгомери М (х, у) - функции, которая характеризует давление на поверхности раздела слоев 2. Вследствие квазигеострофичности движений рассматриваемого масштаба пространственное распределение £ и М дает хорошее представление о крупномасштабных и мезомасштабных течениях. В верхнем слое моря вода движется вдоль изолиний в нижнем слое - вдоль изолиний М. Движение воды происходит таким образом, что большее давление находится справа от направления течения. Чем гуще расположены изолинии, тем большую скорость имеет поток воды.
В мгновенном поле £ (рис. 2, а) хорошо проявляются известные особенности черноморской циркуляции. В верхнем слое моря наблюдаются меандриру-ющее ОЧТ и мезомасштабные антициклонические вихри справа от него. В нижнем слое (рис. 2, Ь) движение воды происходит в циклоническом направлении так же, как в верхнем слое, но с меньшей скоростью (это видно по размаху шкалы М). Круговое течение в нижнем слое сильно меандрирует, причем меандры вытягиваются вдоль изобат и движутся в направлении против часовой стрелки как волны, захваченные материковым склоном [9].
Следующие два поля (рис. 2, с, получены осреднением за большой период времени мгновенных полей М. Они отражают крупномасштабную циркуляцию в бассейне. Согласно рисунку, в верхнем слое СТ представляет собой круговое течение, расположенное над материковым склоном (среднее ОЧТ), а в нижнем слое - течение, у которого линии тока близки к изобатам.
2 Шакина Н. П. Лекции по динамической метеорологии. М. : ТРИАДА ЛТД, 2013. 160 с. МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 40 № 1 2024 113
Р и с. 2. Мгновенные поля Z (a) и M (b), средние поля Z (с) и M (d), отклонения мгновенных значений от средних Z - Z (e) и M - M f)
F i g. 2. Instantaneous fields Z (a) and M(b), average fields Z (c) and M (d), deviations of the instantaneous values from the average ones Z - Z (e) and M - M (f)
На рис. 2, e, f приведены отклонения мгновенных полей Z, M от средних. На этих рисунках хорошо видны вихревые структуры (МВ) в нижнем и верхнем слоях моря. Эти вихревые образования связаны между собой и перемещаются в циклоническом направлении с одинаковой фазовой скоростью [20].
Баланс механической энергии в двухслойной модели
Используемая в работе двухслойная вихреразрешающая модель является энергетически сбалансированной и позволяет в достаточно простой форме описать процессы обмена механической энергией между верхним и нижним слоями моря. Механическая энергия в двухслойной модели состоит из кинетической энергии течений в отдельных слоях и потенциальной энергии.
Кинетическая энергия столба воды единичного сечения в верхнем и нижнем слое K\, K2 определяется как
2 2 2 2 K _ h K _ h _ U + V1 _ u2 + V2
K1 _ h1eK1, K2 _ h2eK2, eK1 _ P1 2 ' eK2 _ P2 2 '
где еК1, еК2 - кинетическая энергия единицы объема воды верхнего и нижнего слоя; р1, р2 - плотность воды в слоях.
Потенциальная энергия положения столба воды единичного сечения от поверхности моря до дна равна
¡1+<;
(
РЕ = -
¡р(г^г = - | р^гЖ - | p2gzdz = р2
¡1+<;
2 2 2 Х2
где Н (х, у) - глубина моря.
Вертикальная ось координат Z направлена вниз, за начало отсчета (г = 0) принимается невозмущенная поверхность моря. В этом случае потенциальная энергия отрицательна и равна работе, которую нужно совершить, чтобы поднять всю воду на нулевой уровень.
Иногда вместо потенциальной энергии удобнее пользоваться доступной потенциальной энергией (Р), которая представляет собой разность между текущей потенциальной энергией (РЕ) и энергией наиболее устойчивого состояния жидкости (Р0). Для двухслойной жидкости наиболее устойчивым является состояние покоя, при котором поверхность моря и граница раздела слоев горизонтальны, поэтому
Р = РЕ - Р0 = р2
( С2 И2 И2 ^
7+g' ¡2-+g«-
где И0 - толщина верхнего слоя в состоянии покоя.
, (и2-Ир)
В приближении «твердой крышки» Р = Р2 g 2 .
Энергетика двухслойной модели описывается системой уравнений энергетического баланса, состоящей из двух уравнений для кинетической энергии (верхнего и нижнего слоя) и уравнения для потенциальной энергии. Уравнения энергетического баланса К1 и К2 можно получить, если для каждого слоя сложить уравнения движения, умноженные на соответствующие составляющие скорости течений. Уравнение баланса потенциальной энергии получается дифференцированием по времени уравнения для Р с последующим использованием приближений Буссинеска и «твердой крышки»:
К
1 | 8Ц1ек1
сх
8Уе
+ —= Ж + Ж + ж + ж + ж
5у
гК гШ е ЛУ е
иг^ + ии = ж + ж + Ж + Ж + Ж
С2 ^ "02 ^ "КЬ2 ^ "ко ^ "ЛВ2'
йх ду
(1)
8Р_
т
(
-рg
ди 2 И, дУ2 /
йх су
Л
г
рg
Зх ду
Л
= _ ж - Ж ''01 02'
где Жа = Ц • /V! -у • /и, = 0; ЖС2 = Ц2 ■ /V2 -У2 • /и2 = 0; Жг = р{и,тх + Vlту); жт = Р (и^х +VlЯLУl); Ж^ = р^ + V2ЯLУ2); жю = р(и2Яу +v2ЯDy);
Ж = Р^ЦСх + Щу); Ж =РЕ(и2^х +Г2^)+РЕ'{и2Нх +У2Ку);
Ж* = рАв (Ц (ААи1)+V; (ААУ1)); = рАв (и2 (Ми,) + V ( Д^ ));
р - средняя плотность воды в Черном море.
Левая часть уравнений (1) включает в себя локальные производные по времени и дивергенции адвективных потоков соответствующих энергий.
В правой части уравнений находятся члены, описывающие работы в единицу времени (мощность) сил, входящих в уравнения движения, которые определяют соответствующие потоки энергии: Жс1, ЖС - работа силы Кориолиса в верхнем и нижнем слое; Ж - работа касательного напряжения ветра, равная притоку кинетической энергии в море (ветровая накачка); Жи, Жб - работа сил трения (диссипация энергии вследствие трения на нижней границе верхнего и нижнего слоев); Жь2 - работа сил трения (поступление энергии в нижний слой из верхнего за счет трения между слоями); Ж1, Ж2 - работа градиента давления (переход между кинетической и потенциальной энергией в верхнем и нижнем слое); Жв1, Жв2 - работа сил горизонтальной турбулентной вязкости в слоях (диссипация энергии).
Потоки Ж1, Ж2 с разными знаками входят одновременно в уравнения баланса кинетической и потенциальной энергии, тем самым обеспечивая энергетическую связь между верхним и нижним слоями.
Работа силы Кориолиса Ж1 и Ж2 в уравнениях (1) равна нулю, но, как будет показано ниже, она может быть ненулевой отдельно для СТ и МВ. Поэтому эти члены в уравнениях сохранены. В этом случае работа силы Кориолиса определяет переходы кинетической энергии между СТ и МВ.
Энергетические характеристики СТ и МВ
Итак, согласно поставленной задаче, требуется провести разделение энергетических характеристик течений, полученных с помощью численной модели, по масштабам движения, а именно: рассчитать средние по времени значения энергии и ее переходов отдельно для средних течений и мезомасштаб-ных вихревых образований. Также необходимо определить, каким образом происходит обмен энергией между СТ и МВ. Расчеты проведем отдельно для восточной и западной части моря.
Рассмотрение средних характеристик дает возможность исключить в уравнениях энергетического баланса производные энергии по времени и тем самым уменьшить количество факторов, влияющих на изменчивость энергетических характеристик. Это несколько упрощает анализ.
Перепишем осредненные за большой период времени уравнения (1) в терминах переходов энергии, используя для обозначения среднего по времени перехода одного вида энергии в другой конструкцию {Е1, £2} [21]. Если {Е1, Ег} > 0, это означает, что Ег переходит в Е\, и наоборот: если {£1, Е2} < 0, то Е1 переходит в Е2:
ди,ек1 : ду1ек1
дх ду
ди2еК2 . дУ2еК2
дх ду
- РВ' ди2\ | д дх V
= жс2 + {к2, р} + {к2, к1} + {к1, лн} + {к15 (2)
ду
РВ
ГдЩ | дУ^ ^ дх ду
= ~ {к1,р}"{к2,Р},
где Жа = 0; = 0; {К1,т} = Ж ; {К1,Р} = Жв, ; {К2,Р} = Жв2;
{К1,^ь} = Ж^; {К2,К1} = {К1,ЛН} = ЖБ; {К^Н ЖДВ1;
|К1, £>Т2} = ЖАВ2; В - диссипация вследствие трения на нижней границе верхнего слоя; Пи - диссипация вследствие придонного трения; Вл, Вт2 - диссипация вследствие горизонтальной турбулентной вязкости в верхнем и нижнем слоях; черта сверху означает осреднение по времени.
Проинтегрируем уравнения (2) по пространству в пределах западной и восточной областей, после чего получим две системы уравнений, в которых для записи суммарных по площади энергетических характеристик применим квадратные скобки с верхним индексом W или Е, означающим область интегрирования:
[...Г = {{(...) <^хёу, [,..]Е ={{(...) ёхёу,
ди,е„, дУ,вг
И1 дх
ду
Я1
W
Я
дипе„п дУе
дх
су
( Гди2И1 +дУ2И1Л
Р8
V V
Я|
Е
III
Е
я
дх ду
Л
^у = [Жс1 ]" +[К ,Р]" +[К ,г]" +[К А 1 Г,
\су = [ЖС2]" +[К2,Р]" + [К2,К1]Ш +[К1,Ви]" +[К1,ДГ2]", (3) сксу =- [К^Р]" ~[К2,Р]™,
-рв
/у
аи!еК1+ау;еК1
сх
ду
дх ду
dxdy = [Жа ]Е +[К„ Р]Е +[К„т]Е +[К„ ]Е +[К„ £>„ ]Е
дОпе„п дУе
Л
дх
( (
рв
V V
Эу
СХСУ = [ЖС2 ]Е +[К2,Р]Е +[К2,К1 ]Е +[К,,Ви ]Е +[К,,ВТ2 ]Е, (4)
ди2Н1 дУ2И1 дх ду
^ ( рв
дЩ | У дх ду
\Л
скс!у = - [К1,Р]Е -[К2, Р]Е.
Используя формулу Остроградского - Гаусса и условия прилипания на твердой границе, заменим стоящие в левой части уравнений двойные интегралы по площади от дивергенции адвективных потоков энергии линейными интегралами от х-составляющей этих потоков вдоль границы Б:
f FkL( y)dy = ffi ^ + dIeL 1 dxdy = -ffi ^KL + W^L 1 dxdy, f F2( y)dy = ffi ^^ + 1 dxdy = -ffi + 1 dxdy,
D K JW l &_ay ^ J л ay J J ()
J FP( y)dy = JJÎ^,g + M1-pg te11 dxdy =
D ^ ^ ^ dy J { âx dy JJ
ccf ,(dûh dVh ) (дЩ ÔVZ)). ,
= -III-pg +—2r -pg +dxdy,
W ^ dy ) y dx dy JJ
где Fkl = UjeK1, FK2 = иею, Fp = -p(g'U2h + g(UL + U2 )Ç) - средние адвективные потоки кинетической и потенциальной энергии между западной и восточной частями моря.
Чтобы разделить энергетические характеристики по масштабам движения, воспользуемся данным выше определением для СТ и МВ. Как уже говорилось, СТ получается осреднением по времени мгновенных полей течений, а МВ -это отклонения мгновенных течений от СТ. Для идентификации СТ и его энергетических характеристик используем верхний индекс m, для МВ - индекс p.
Математически разделение циркуляции по масштабам движения можно выразить через течения (потоки воды) в слоях следующим образом:
Ui = um + Ul = Uh +(uh )', v = Vm + Vp = Vh+ {vh )',
_ _ (6)
U2 = U2m + U2p = u2h + (u2h2 )', V2 = V2m + V2P = v2h2 + (v2h2 )',
где штрих означает отклонение (пульсацию) от среднего значения; (UL,V), (U2,V2) -составляющие течений (потоков) в слоях; (ULm,Vm), (U22,V2m) - составляющие среднего течения; (UP,VP), {Up,V2p) - составляющие пульсационных течений, или МВ.
В данной работе период осреднения течений был выбран равным периоду осреднения уравнений (2).
Если подставить выражения (6) в систему уравнений (1), то после раскрытия скобок и разделения членов уравнений по масштабам движения можно получить отдельные системы уравнений энергетического баланса для СТ и МВ.
Покажем, как работает разделение по масштабам движения на примере члена, описывающего работу силы Кориолиса в первом уравнении системы (1):
wCi = u h • А - vi h • fui =(um + Uip )• & -(V m + vp )• fui , wCi =(Um ■ V - Vim • fu ) + (up • fvi - Vip • u ) = + W,
где WCi = Um • fvi - Vm • fui и WCP = Up • fvi - Vp • fui - работа силы Кориолиса, совершаемая СТ и МВ соответственно.
На рис. 3 приведен вектор мгновенного потока, разложенный на средний и пульсационный потоки, и вектор силы Кориолиса, действующей на мгновенный поток.
Р и с. 3. Возникновение ненулевых работ силы Кориолиса ( W™ = -WC ф 0 ) при разложении течений на сумму средних ulhl и пульсационных (иД)'
F i g. 3. Arising of non-zero values of the Coriolis force work (W™ = -W<P1 ф 0 ) at decomposing the flows into a sum of the average uh and pulsation (иД)' values
Как видно из рисунка, разделение течений (потоков) на средние (СТ) и пульсационные (МВ) приводит к появлению ненулевых работ силы Кориолиса, с ними связанных. В сумме эти работы равны нулю, т. е. взаимно уничтожаются. Но при разделении уравнений энергетического баланса реальных потоков на уравнения для СТ и МВ соответствующие работы силы Кориолиса должны быть учтены для согласования систем уравнений. Правильно считать эти работы мнимыми (фиктивными), поскольку разделение течений на средние и пульсационные также является мнимым, существующим только в рамках используемого метода представления.
В уравнениях баланса энергии СТ и МВ работы силы Кориолиса определяют переходы кинетической энергии между средними и пульсационными течениями. Можно сказать, что не существующие на самом деле средние и пуль-сационные течения обмениваются между собой кинетической энергией за счет мнимой работы силы Кориолиса, которая сама является мнимой (фиктивной) силой.
Учитывая, что уравнения (3) - (5) являются следствием уравнений (1), получим для западной и восточной частей Черного моря две системы уравнений, в которых члены, относящиеся к СТ и МВ, разнесены в разные уравнения:
фт = Г пЖ Г пЖ г пЖ г пЖ г пЖ [кг,крр] +[кг,рт] +[кг,г] +[кг,ЮГ] +[КГ,БГ] ,
фр = г пЖ Г г пЖ г г = [крр,кг] +[кр,рр] +[крр,г] +[крр,юг] +[кррЮтр1 ] ,
фт = К2 г пЖ г -]Ж г пЖ г -|Ж г пЖ -- [кг,кр] +[К2Г,РГ] +[кг,кг] +[кг,юг] +[кг,юг2]
Фр = К2 = [ кр, кг ]Ж + [кр, рр ]Ж + [кр, крр ]Ж + [кр, юн ]Ж +[ кр, юр2 ]Ж,
Ф г = Г пЖ г !Ж - [кг,рг] -[кг,рг] ,
Фрр = - [кр,рр]Ж -[кр,рр]Ж,
"ФГг = [кг, кр ]Е +[ кг, р г ]Е + [кг,г]Е + [ к;, ю; ]Е +[ кг,БГ; ]Е,
-ФРК1 = [ кр, к; ]Е + [кр, рр ]Е + [кр,г]Е +[ кр, юр ]Е +[ кр, ют; ]E,
-фт К2 = [ к;, кр ]Е + [ к;, р г ]Е + [ к;, к; ]Е + [ к2Г, ю; ]Е +[ к;, ю; ]Е
— фр К2 = [ кр, кг ]Е + [кр, рр ]Е + [кр, кр ]Е +[ кр, юн ]Е +[ кр, юр2 ]Е,
-Ф; = - [кг,рг]Е -[к;,Рг]Е,
"ФРр : = - [ кр, р р ]Е -[ кр, р р ]Е,
где ФК! = |; ФК2 = |Е™ёу; ФГ = {ЕрГёу; ФК = {ЕКрёу; Ф|2 = {Е£ёу;
Б Б Б Б ю
Фр - Г ЕрёУ • Ет + Ер - • Ет + Ер - • Ет + Ер - Е • Ет - Г/т • •
Б
ер - ^ - Ет • Ет - Г/т • • Ер - - Ет •
К1 К1 К1 ' 1 К2 2 СК2 ' 1 К2 1 К2 1 К2 з
Ерт = и; К + g(стт + иг)С); Ерр = ЕР -ЕРт ;
Г пШ Г пШ Г пШ
[к;кр] =[жст] = [цг• /У1 -^¡т• /й;] ;
г пШ г пШ Г пШ г г — — пШ
[кр\к;] =[жр1 ] = -[>;,кр] ; [кг,рг] = р^[ит • с, + *Рт• су] ; [кр ,,р 7 = Р£ \ЩХ + ^ 7 - [кг,,г 7;
Г пШ Г -Т -ПШ Г Г- -г пи
[КГ,т] = р|_ЦГ • т^1 + ^Рт • А ] ; [кр,т] = р[мРт* + уу ] "[^1 ;
[ кт, юр ]Ж =Р . ^+Грт. «ПК7 Г;
Г пШ Г- -Г пШ
[кр,] =4+ ^] -[к*,] ; [кг, Ютт17 =РЛВ [ ит (ААЙр) + ГГ (АДУ; )]Ж; [кр, ]Ж =рАв [и, {Мм 1)+^ (АДу1 )]Ж -Цкг,юотр ]Ш; [К2т,кр]Ж =[жт2]Ш = [цт • /У2 -Г2т • /й2]Ж ;
Г пШ г пШ г п«1
[ к2р, к 2т ] =[Жр2 ] ="[ к2т, к2р ] ;
[к2т,Рт]Ж = р[я(щ • СТ + г2т • £)+ я'(^2т • К + • )]";
[ Ц,РР ]№ =р[8 (и2 с х + у2 с,) + Я '(и2 к + г2 И у - [ кт,Рт ] ^;
[кг, кт ]" = р . ^+V? • Г;
Г Г- -ПШ Г П'М
[к>р] =Р[«х2+] -[кг,] ; [ к2т, р" ]ш = р |>2т • кх + • Г;
г Г- -Пш г ^w
[кт,] + ] -[кт,Рн] ;
[К2т, Р^ ]W = рлв и (АМ2) + V" (АА?2 )]W; [ к?,, ] W= Р4 |>2 (ААм 2) + Г2 (МУ2 [ кт,Ртт2 ] W•
Формулы для вычисления переходов энергии в восточной половине моря получаются заменой W на Е в вышеприведенных выражениях.
Выражения [к™,кр^ , [кР,к"]W , [к",кр]W , [к2Р,к"]W , [к^,кр]Е,
[кР,кр"] , [к2т ,к2 ] , [к2, кт ] не равны нулю и определяют переходы кинетической энергии между СТ и МВ в соответствующей части моря.
Результаты расчета энергетических характеристик
Прежде всего были определены средние по времени значения энергий крупномасштабных течений и мезомасштабных вихревых образований в западной и восточной частях моря, для чего использовались формулы для кр, кг, Р, выражения (6) и представление И = ИТ + ИР (где ИТ = И1):
к ии+ИУ {и?щ + утУр )+(и?Ы1 + у\)
к1 - Иек1 - р—2— ~ р-2-'
к ии + И У2 (и2ти2 + V" У2 ) + {и2РМ2 + У2рУ2 )
к2 = И2еК2 = р-2- = р-2-'
Р = р§, /¿^ = к Хит + И )2 - И2 = ^, И )2+2ИтИР + (ИР )2 - И2
2 2 2 Осреднение по времени и разделение по масштабам движения дает
к (иртщ + утц ) + (цри1 + ) (цч + урт ур) (ич + уч) к1 =Р-2-= р-2- -2-'
кр=к1т + к1Р, к1т = рКА^Л., кр =р^ •1 = к1 -к1т,
V 2 ' ' 2 ' 2 / V 2 "2 ' ' 2 ' 2 / '
к2 =Р-^- = Р"
(ЦР и + VlРVl)
2
й2 + V/ л >,р<
2 1 Р
(С/2Р и2 + ^2Р ^2 )
к 2 = к т + к 2Р, к2т = рити2+ , к 2Р =р ^ = к 2 - кт,
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 42 № 1 2224 121
- Дг )2 - Ко2 +(Ир )2 ,(ИГ )2 - Ко2 )2 р = ря—20 = рял0 + РЯ,
- К )2 _ К2 (КЛ^ _
р = рг + рр, рг = ря—-, рр = ря= р -рг.
После интегрирования по площади для западной и восточной частей моря осредненных по времени энергий СТ и МВ получим искомые значения (таблица).
Осредненные по времени кинетическая и доступная потенциальная энергии (ТДж) Time-averaged kinetic and available potential energy (TJ)
! К \ ' К 2m " [ Kp ] ! К2 ! " К 2m " [ К ] P [ P m " [ p p ]
Западная часть / Western part
313,7 238,3 75,4 48,4 1,9 46,4 881,1 756,0 125,0
Восточная часть / Eastern part
370,7 286,6 84,1 61,7 2,2 59,4 2915,0 2773,2 141,8
Средние по времени энергии стационарны и не влияют на переходы энергии, но их знание полезно для общего понимания энергетики Черного моря и сравнения результатов моделирования с данными наблюдений.
Согласно полученным данным, кинетическая энергия СТ и МВ распределена в Черном море достаточно равномерно между его восточной и западной частями. Немного более высокие значения кинетической энергии в восточной половине моря можно объяснить большей протяженностью крупномасштабных течений в этой части бассейна.
В верхнем слое моря кинетическая энергия СТ больше кинетической энергии МВ примерно в три раза, а в нижнем слое - наоборот: в МВ содержится в 25 раз больше кинетической энергии, чем в СТ.
Больше всего средней механической энергии в море (~ 60%) приходится на доступную потенциальную энергию СТ, сосредоточенную в его восточной половине. Она почти в 10 раз превосходит кинетическую энергию СТ в этой части бассейна. В западной половине моря общая доступная потенциальная энергия СТ в три раза превышает кинетическую энергию СТ. При этом она в три раза меньше, чем доступная потенциальная энергия СТ в восточной половине.
Для МВ разница в распределении Р между западной и восточной частями пропорциональна длине среднего ОЧТ в соответствующей области.
Следует обратить внимание на выявленную закономерность: в сумме кинетическая энергия МВ верхнего и нижнего слоев получилась примерно равной доступной потенциальной энергии МВ:
К Г +[ кр ]W ]W, [кр ]E +[ к p ]E ]E.
МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 40 № 1 2024
На следующем шаге вычислялись средние адвективные переносы энергии между западной и восточной частями моря, осуществляемые СТ и МВ. Именно эти потоки, возникающие вследствие Р-эффекта, как было показано в работе [17], являются причиной образования пространственной неоднородности полей энергетических характеристик в Черном море.
На рис. 4 показано распределение адвективных потоков энергии по оси У. В верхнем слое моря максимальный средний перенос кинетической энергии отмечается в стрежне среднего ОЧТ (рис. 4, а). Направление переноса энергии совпадает с направлением течения. Основной вклад в адвекцию кинетической энергии вносит СТ (рис. 4, Ь), максимум совпадает со стрежнем среднего потока.
Р и с. 4. Распределение адвективных потоков энергии вдоль границы D: средние потоки -сплошная линия, доля СТ - штриховая, доля МВ - пунктир
F i g. 4. Distribution of the energy advective flows along boundary D: mean flows are shown by a solid line, share of mean currents - by a dashed line, share of mesoscale eddies - by a dotted one
Южная ветвь СТ переносит больше энергии, чем северная, в результате суммарный поток кинетической энергии Фт1 составляет 33 МДж/с и имеет направление из западной половины моря в восточную (рис. 5, а).
Р и с. 5. Суммарные адвективные потоки энергии (МДж/с) между восточной и западной частями моря: Ф™, Ф^ (а); Ф» , Ф(b); ф£, Ф£ (с)
F i g. 5. Total advective energy flows (MJ/s) between the eastern and western parts of the sea: Ф^ ,
Ф^(а); Фm2, Ф^2 Ш ФP (с)
Мезомасштабные вихревые образования в верхнем слое моря осуществляют перенос кинетической энергии так же, как СТ, в направлении среднего ОЧТ, но максимумы потоков энергии МВ не совпадают со стрежнем среднего потока, а расположены рядом с ним в зонах наибольшего градиента скорости течений (рис. 4, с). Причем справа от стрежня СТ перенос энергии F¡pi больше, чем слева. Вероятнее всего, это происходит потому, что мезомасштабные антициклонические вихри, которые образуются между ОЧТ и берегом, содержат в себе больше энергии, чем циклонические МВ, образующиеся в основном слева от ОЧТ. Максимальный адвективный поток кинетической энергии МВ наблюдается правее северной ветви СТ, он вносит основной вклад в суммарный адвективный перенос Ф^, который составляет 14 МДж/с в западном направлении (рис. 5, а). Это меньше, чем Ф^, поэтому общий адвективный перенос кинетической энергии течениями в верхнем слое направлен из западной половины моря в восточную.
В нижнем слое моря южная ветвь СТ переносит в восточном направлении больше кинетической энергии, чем северная ветвь в западном (рис. 4, e). Но общий перенос кинетической энергии (рис. 4, d) оказывается все же направленным на запад за счет адвективного потока МВ, который максимален в северной части разреза D (рис. 4, f). Суммарный адвективный поток кинетической энергии СТ Фm2, направленный на восток, равен 0,4 МДж/с, а суммарный поток кинетической энергии МВ Ф|2 направлен на запад и равен 0,5 МДж/с (рис. 5, b).
Средний адвективный перенос потенциальной энергии СТ происходит в направлении, противоположном потоку, в отличие от переноса кинетической
энергии (рис. 4, g, h). Максимальные значения FP и FPm отмечаются в стрежне СТ нижнего слоя. Они намного превышают максимумы потоков кинетической
энергии СТ в верхнем и нижнем слоях. Но суммарный поток Ф* при этом оказывается одного порядка с ФГ и составляет 25 МДж/с в западном направлении за счет южной ветви СТ (рис. 5, с).
Наиболее интересным представляется распределение вдоль разреза Б адвективного потока ЕРр (рис. 4, /). Он имеет максимумы в области стрежня СТ в нижнем слое моря и направлен по потоку. Правее стрежня СТ поток ЕРр меняет направление на противоположное, достигает максимума и обращается в ноль на берегу. Можно предположить, что такое распределение ерр связано с образованием и движением справа от СТ мезомасштабных антициклонических вихрей (рис. 1, /). В [9] показано, что такие вихри являются элементами градиентно-вихревых волн (захваченных материковым склоном), которые относятся к классу топографических волн Россби. В целом потоки ЕРр между стрежнем СТ и берегом вследствие разнонаправленности компенсируют друг друга.
В центральной части разреза Б (рис. 4, /, между 200 и 300 км на вертикальной шкале) существует еще один максимум адвективного потока ЕРр , который и определяет основной вклад МВ в общий перенос потенциальной энергии, в результате чего суммарный поток Фр = 40 МДж/с направлен на запад (рис. 5, с).
В целом можно сказать, что МВ переносят из восточной половины Черного моря в западную и кинетическую, и потенциальную энергию, тогда как СТ переносят на запад только потенциальную энергию. Наибольший вклад в западный перенос энергии вносит адвективный поток Фр (рис. 5, с), который формируется вследствие переноса потенциальной энергии мезомасштабными вихревыми образованиями через центральную часть разреза Б. Такими МВ, согласно работам [9, 22], являются баротропные волны Россби замкнутого бассейна, образующиеся в глубоководной части моря вследствие Р-эффекта. Как показано в 3, в бассейне с двухслойной жидкостью при стационарном ветровом воздействии происходит генерация именно баротропных волн Россби. Если для возбуждения движения в модели использовать сезонно меняющийся ветер, то в бассейне кроме баротропных волн будут формироваться и бароклинные волны Россби 3 [23].
Рассмотрим средние по времени переходы механической энергии, являющиеся следствием динамических процессов в море, с учетом разделения течений на СТ и МВ. Эти переходы, рассчитанные для западной и восточной частей Черного моря, показаны в виде схемы на рис. 6. Период осреднения по времени составил 20 модельных лет, этого оказалось достаточно, чтобы временная изменчивость средних энергетических характеристик оказалась близкой к нулю.
3 Сафронов Г. Ф. Возбуждение длинных волн в океане крупномасштабными изменениями в поле касательного напряжения ветра. М. : Гидрометеоиздат,1985. 108 с. МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ том 40 № 1 2024 125
Р и с. 6. Средние переходы энергии между СТ и МВ с учетом деления моря на две части, западную и восточную
F i g. 6. Mean energy flows between the mean currents and the mesoscale eddies with the regard for the sea division into two parts - the western and eastern ones
Отметим ряд особенностей в переходах энергии.
1. Ветровая накачка энергией происходит в основном в восточной части моря (94 МДж/с), что объясняется особенностями используемого в модели поля ветра и соответствует наблюдениям [14]. В западной части приток ветровой энергии составляет 11 МДж/с, что в 9 раз меньше, чем в восточной.
2. Вся приходящая от ветра энергия идет на пополнение кинетической энергии CX в верхнем слое. На МВ ветер в среднем оказывает тормозящее воздействие, что приводит к потере кинетической энергии МВ как в восточной (6 МДж/с), так и в западной части моря (4 МДж/с).
3. В верхнем слое восточной части моря МВ получают кинетическую энергию от СТ, а в западной части наблюдается обратный переход кинетической энергии от МВ к СТ.
4. В верхнем слое восточной части моря кинетическая энергия СТ и МВ переходит в потенциальную энергию СТ и МВ, в западной части - противоположное направление переходов энергии от потенциальной к кинетической.
5. В нижний слой моря энергия поступает за счет перехода доступной потенциальной энергии МВ в кинетическую энергию МВ как в восточной, так и в западной части моря.
6. Кинетическая энергия СТ нижнего слоя в обеих частях моря пополняется за счет кинетической энергии МВ, т. е. наблюдается переход энергии от движений мелкого масштаба к более крупным. В теории турбулентности такой переход энергии называется эффектом отрицательной вязкости, но в данном случае это результат осреднения движущихся над материковым склоном захваченных волн, которые вследствие нелинейности переносят некоторую массу воды, т. е. имеют свойства вихрей.
7. В восточной половине моря часть кинетической энергии СТ нижнего слоя переходит в потенциальную энергию СТ, а в западной половине - наоборот: потенциальная энергия СТ частично переходит в кинетическую энергию СТ.
8. Несмотря на то что приток энергии в море находится преимущественно в его восточной половине, диссипация энергии на западе больше, чем на востоке (53 и 42 МДж/с соответственно). Причем больше всего энергии, 67 МДж/с, теряется за счет диссипации в верхнем слое, при этом 43 МДж/с приходится на западную часть моря.
На рис. 6 не обозначены переходы кинетической энергии из верхнего слоя в нижний вследствие их малости. Соответствующие им значения учтены в диссипации энергии.
Важной особенностью представленной схемы, связывающей воедино все ее элементы, является наличие рассмотренных выше адвективных потоков энергии, направленных из одной половины моря в другую (рис. 5). С одной стороны, можно сказать, что данные адвективные потоки компенсируют различия в размере и направлении энергетических переходов в западной и восточной частях моря. Но, с другой стороны, эти адвективные потоки энергии, являющиеся следствием Р-эффекта, и есть главная причина неравномерного распределения по акватории моря переходов энергии.
Заключение
Рассмотрение энергетических характеристик отдельно для восточной и западной части Черного моря позволило выявить ряд закономерностей в энергетике двухслойной модели, получить новую информацию о протекании динамических процессов в море. В частности, оказалось, что направление и размер осредненных за большой период времени переходов механической энергии между ее видами в различных частях моря значительно отличаются друг от друга.
Были установлены различия в направлении и размере потоков энергии в восточной и западной частях моря. Их причиной оказались адвективные потоки энергии, которые в сумме направлены из восточной половины моря в западную. В результате вычисления вклада средних течений и мезомасштабных вихревых образований в средние по времени адвективные потоки энергии получено, что основной вклад в западный перенос энергии вносят МВ, локализованные в центральной части моря вне зоны распространения крупномасштабных течений. На этом основании сделан вывод, что данные МВ являются волнами Россби, которые образуются вследствие Р-эффекта в замкнутом (ограниченном) бассейне. Этот вывод еще раз подтверждает важность Р-эффекта для формирования черноморской циркуляции.
Показана роль работы силы Кориолиса в переходах кинетической энергии между СТ и МВ. Возникает резонный вопрос: как такое возможно? Известно, что сила Кориолиса является мнимой силой и не производит работы. Чтобы разрешить данное противоречие, напомним, что разделение течений на СТ и МВ искусственное, т. е. мнимое, поэтому и переходы энергии между ними мнимые. Эти мнимые переходы определяются работой мнимой силы Корио-лиса, которая для суммарных течений равна нулю, но получилась не равной нулю отдельно для СТ и МВ.
На данном этапе исследований анализировались энергетические характеристики, осредненные за большой период времени в статистически равновесном режиме, в котором все средние характеристики и параметры модели остаются постоянными. Такой подход значительно упрощает энергетический анализ, так как позволяет исключить производные по времени в уравнениях энергетического баланса. Но при этом из рассмотрения исключаются некоторые важные процессы, которые имеют стохастический или периодический характер и происходят в различных частях моря (напр., бароклинная неустойчивость). Для исследования подобных процессов при определении средней циркуляции необходимо выбирать соответствующий период осреднения и учитывать пространственно-временную локализацию самого явления. В этом направлении и планируется проведение дальнейших исследований.
СПИСОК ЛИТЕРАТУРЫ
1. Иванов В. А., Белокопытов В. Н. Океанография Черного моря. Севастополь, 2011. 212 с.
2. Seasonal, interannual, and mesoscale variability of the Black Sea upper layer circulation derived from altimeter data / G. Korotaev [et al.] // Journal of Geophysical Research: Oceans. 2003. Vol. 108, iss. C4. 3122. doi:10.1029/2002JC001508
3. Коротенко К. А. Моделирование мезомасштабной циркуляции Черного моря // Океанология. 2015. Т. 55, № 6. С. 909-915. EDNUVEPDV. doi:10.7868/S0030157415060076
4. Universal structure of mesoscale eddies in the ocean / Z. Zhang [et al.] // Geophysical Research Letters. 2013. Vol. 40, iss. 14. P. 3677-3681. doi:10.1002/grl.50736
5. Chelton D. B., SchlaxM. G., Samelson R. M. Global observations of nonlinear mesoscale eddies// Progress in Oceanography. 2011. Vol. 91, iss. 2. P. 167-216. doi:10.1016/j.pocean.2011.01.002
6. Chen G., Hou Y., Chu X. Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure // Journal of Geophysical Research: Oceans. 2011. Vol. 116, iss. C6. C06018. https://doi.org/10.1029/2010JC006716
7. Кубряков А. А., Станичный С. В. Синоптические вихри в Черном море по данным спутниковой альтиметрии // Океанология. 2015. Т. 55, № 1. С. 65-77. EDN TGWBIX. doi:10.7868/S0030157415010104
8. Numerical Models of Oceans and Oceanic Processes / Ed. L. H. Kantha, C. A. Clayson. San-Diego : Academic Press, 2000. 940 p. (International Geophysics Series ; vol. 66). doi:10.1016/s0074-6142(00)x8001-1
9. Павлушин А. А., ШапироН. Б.,Михайлова Э. Н. Роль рельефа дна и Р-эффекта в динамике Черного моря // Морской гидрофизический журнал. 2017. № 6. С. 27-39. EDN QHQKPF. doi:10.22449/0233-75 84-2017-6-27-39
10. Субмезомасштабные вихри на кавказском шельфе Черного моря и порождающие их механизмы / А. Г. Зацепин [и др.] // Океанология. 2011. Т. 51, № 4. С. 592-605. EDN NXXEQF.
11. Mesoscale to Submesoscale Transition in the California Current System. Part II: Frontal Processes / X. Capet [et al.] // Journal of Physical Oceanography. 2008. Vol. 38, iss. 1. P. 4464. doi:10.1175/2007JPO3672.1
12. Влияние субмезомасштабных вихрей на перенос взвешенного вещества в прибрежной зоне Крыма по данным БПЛА, спутниковых и контактных измерений / А. А. Кубряков [и др.] // Океанология. 2021. Т. 61, № 2. С. 182-197. EDN VZKVBX. doi:10.31857/ S0030157421020106
13. Калашник М. В., Курганский М. В., Чхетиани О. Г. Бароклинная неустойчивость в геофизической гидродинамике // Успехи физических наук. 2022. Т. 192, № 10. С. 11101144. EDNPLRMZL. doi:10.3367/UFNr.2021.08.039046
14. Бассейновая циркуляция и мезомасштабная динамика Черного моря под ветровым воздействием / А. Г. Зацепин [и др.] // Современные проблемы динамики океана и атмосферы : сборник статей, посвященный 100-летию со дня рождения П. С. Линейкина. Москва : ТРИАДА ЛТД, 2010. С. 347-368. EDN TVXEMZ.
15. Kang D., Curchitser E. N. Energetics of Eddy-Mean Flow Interactions in the Gulf Stream Region // Journal of Physical Oceanography. 2015. Vol. 45, iss. 4. P. 1103-1120. doi:10.1175/JPO-D-14-0200.1
16. Демышев С. Г., Дымова О. А. Анализ энергетического цикла Лоренца для различных режимов циркуляции Черного моря // Труды Карельского научного центра РАН. 2022. № 6. С. 26-40. EDN QAWITM. doi:10.17076/lim1621
17. Павлушин А. А. Особенности и причины пространственной неоднородности потоков механической энергии в Черном море // Морской гидрофизический журнал. 2023. Т. 39, № 3. С. 328-341. EDN OSWDCN.
18. Павлушин А. А., Шапиро Н. Б., Михайлова Э. Н. Захваченные волны и меандрирование Основного Черноморского течения // Экологическая безопасность прибрежной и шель-фовой зон моря. 2019. № 4. С. 14-21. EDN SYNPMJ. doi:10.22449/2413-5577-2019-4-14-21
19. МарковаН. В., Багаев А. В. Оценка скоростей глубоководных течений в Черном море по данным дрейфующих буев-профилемеров Argo // Морской гидрофизический журнал. 2016. № 3. С. 26-39. EDN WNAFSV.
20. Энергетический спектр скорости течения в глубокой части Черного моря / А. А. Клю-виткин [и др.] // Доклады Академии наук. 2019. Т. 488, № 5. C. 550-554. EDN RUPYJT. doi:10.31857/S0869-56524885550-554
21. Holland W. R., Lin L. B. On the Generation of Mesoscale Eddies and their Contribution to the Oceanic General Circulation. I. A Preliminary Numerical Experiment // Journal of Physical Oceanography. 1975. Vol. 5, iss. 4. P. 642-657. doi:10.1175/1520-0485(1975)005<0642:OTGOME>2.0.CO;2
22. Stanev E. V., Rachev N. H. Numerical study on the planetary Rossby modes in the Black Sea // Journal of Marine Systems. 1999. Vol. 21, iss. 1-4. P. 283-306. https://doi.org/10.1016/S0924-7963(99)00019-6
23. Павлушин А. А. Автоколебания интенсивности крупномасштабной циркуляции в Черном море // Морской гидрофизический журнал. 2022. Т. 38, № 6. С. 605-619. EDN UFLDYX. doi:10.22449/0233-7584-2022-6-605-619
Об авторе:
Павлушин Андрей Александрович, младший научный сотрудник, отдел динамики океанических процессов, ФГБУН ФИЦ МГИ (299011, Россия, г. Севастополь, ул. Капитанская, д. 2), ORCШ ГО: 0000-0002-2098-5068, ResearcherГО: R-4908-2018, 8РШ-код: 1439-2290, pavlushin@mhi-ras.ru
