CC BY
21
УДК 551.465
Роль рельефа дна и ^-эффекта в динамике Черного моря © 2017 А. А. Павлушин*, Н. Б. Шапиро, Э. Н. Михайлова
Морской гидрофизический институт РАН, Севастополь, Россия *E-mail: pavlushin@mhi-ras.ru
Поступила в редакцию 01.06.2017 г.
Обсуждаются результаты численных экспериментов, выполненных в рамках двухслойной вих-реразрешающей модели Черного моря. В модели учитываются рельеф дна, ß-эффект, придонное трение и горизонтальная турбулентная вязкость, параметризуемая бигармоническим оператором. Расчеты проведены на длительный срок до выхода решения на статистически-равновесный режим, который характеризуется наличием интенсивных течений, волн и вихрей. Показано, что в режиме статистического равновесия в море существует круговая циркуляция циклонической направленности, в верхнем слое образуется струйное меандрирующее течение - аналог ОЧТ. В нижнем слое выявлены достаточно интенсивные волны, которые накладываются на поток, распространяющийся вдоль изобат. Эти волны можно характеризовать как топографические волны Россби, захваченные материковым склоном. Предложена методика анализа таких волн. Показано, что волновые возмущения течений в нижнем слое существенно влияют на течения в верхнем слое, способствуя их неустойчивости и меандрированию.
Ключевые слова: Черное море, вихреразрешающая модель, численный эксперимент, ß-эф-фект, топографические волны Россби, рельеф дна, захваченные волны.
DOI: 10.22449/023 3-7584-2017-6-27-39
Введение. В работе в рамках двухслойной вихреразрешающей модели приводятся результаты исследования влияния рельефа дна и ß-эффекта на формирование циркуляции в Черном море при воздействии стационарного циклонического ветра с постоянной завихренностью. Она является продолжением серии численных экспериментов [1, 2], посвященных моделированию ветровой циркуляции в Черном море. В опубликованных ранее работах обсуждаются результаты численных экспериментов по исследованию роли таких факторов, как завихренность ветра, форма бассейна, орография береговой линии, интенсивность придонного трения и горизонтальной турбулентной вязкости, ß-эффект в случаях, когда рельеф дна не учитывался (дно предполагалось горизонтальным).
Модель основана на примитивных уравнениях гидродинамики океана [3]. Уравнения модели представляют собой проинтегрированные по вертикали в пределах каждого слоя уравнения движения и неразрывности
U ) + («U )х + (v,U, )y - fVi = ghCx + rx - RX + ABVfrV(A« )),
(V ) + (uV )x + (vV )y + fU 1 = gh+ ry - Ry + ABV(hV(äVj )), (U2) +(u2U2)x +(v2U2)y - V = gh2^x + gh(h)x + RX -RX + AbV(H2V(M2)) ,
(V )( +(«2V2 )x +(v2V2 )y + fU2 = glh£y + g h (h )y + Ry - RX + AbV^V(AV2 )) ,
(Л ) + (и1 + (У )у = 0, (Л) + (и2 )х +(У )у = 0,
где индексы 1 и 2 указывают номер слоя; х, у, t — дифференцирование по соответствующей переменной; ик, ук — горизонтальные компоненты скорости течений в к-м слое (к = 1, 2), направленные вдоль осей X (на восток) и У (на север) соответственно; Л1, к2 — толщины слоев; ик = икНк , ¥к = УкЪк — компоненты потоков в слоях; Щ, ЯУ — компоненты силы трения на поверхности раздела слоев; Щ, Щ — компоненты придонного трения; / = /0 + /Зу — параметр Кориолиса, /0 = 10~41/с, / = 2 -1013 1/смх; g = 980 г-см/с2 — ускорение свободного падения; g'= g(p~p)/Рт.; Т>Х, тУ — компоненты тангенциального напряжения ветра; Ав — коэффициент горизонтальной турбулентной вязкости.
Для замыкания системы уравнений используется интегральное уравнение неразрывности в приближении «твердой крышки»
их + Уу = 0,
что позволяет ввести интегральную функцию тока и = ~Уу, У = фх, где и = и + и2, У = У + У — составляющие полных потоков.
Сток рек и водообмен через проливы не учитывался, на боковых границах бассейна задавались условия прилипания и условие Ли^ = 0. В начальный момент времени вода находилась в состоянии покоя.
Для численной аппроксимации использовалась двухслойная по времени полунеявная (экономичная явная [4]) схема с неявным учетом силы Кориолиса, придонного трения и трения на поверхности раздела. Уравнения движения и неразрывности аппроксимировались на сетке В со вторым порядком точности по схеме центральных разностей. Отметим, что в предыдущих экспериментах с горизонтальным дном нелинейные члены в уравнениях неразрывности аппроксимировались схемой первого порядка точности (направленными разностями), а в уравнениях движения - схемой Лакса - Вендроффа второго порядка точности [4]. Аппроксимация центральными разностями обеспечивает соблюдение энергетического баланса в модели.
Было показано [1], что вследствие влияния нелинейных факторов даже при воздействии стационарного ветра решение задачи со временем выходит не на стационарный, а на квазипериодический статистически-равновесный режим. При этом приток энергии от ветра в море не является постоянным и регулируется изменчивостью поля поверхностных течений.
Под притоком энергии подразумевается суммарная по площади бассейна работа Жт, совершаемая тангенциальным напряжением ветра, которая равняется скалярному произведению векторов т и и. Поскольку напряжение ветра стационарно, то изменчивость Жт зависит от изменчивости и •
Другой важный момент, на котором стоит заострить внимание, - учет придонного и горизонтального трения. Горизонтальная вязкость в данной версии модели параметризуется бигармоническим оператором, в отличие от предыдущих экспериментов [1, 2], в которых использовался гармонический оператор. Бигармоническая вязкость в меньшей степени препятствует возникновению в модели синоптических и мезомасштабных вихрей [5].
Придонное трение в рассматриваемых экспериментах параметризовалось
формулой = (г + г2 |и2 |)и2 . Коэффициент г1 задавался для того, чтобы придонное трение действовало при слабых течениях. Трение на поверхности раздела слоев не учитывалось (Ка = 0) в случаях, когда Н2 > 0, чтобы в чистом виде описать формирование вихревых структур (и вообще движения) в нижнем слое. В тех случаях, когда ^2 = 0 и верхний слой контактировал
с дном, И я = 01 + Г2 |и1|)и1.
Численные эксперименты. Ниже приводятся результаты двух численных экспериментов R1 и Я2, в которых рассматривается движение в двухслойном море с реальным рельефом дна Н(х, у), возбуждаемое стационарным напряжением ветра т(х, у), имеющим постоянную циклоническую завихренность 1^(т) = 2,53-10-7 Н/м3. В эксперименте R1 учитывается ^-эффект, в Я2 в-эффект не учитывается (в = 0). Составляющие тангенциального напряжения
ветра тх,ту, как и в [1, 2], вычислялись по формулам тх = тX + у(тХ ~ТХ)/В,
ту =ту + х(ту -ту)/Ь , где тХх = туу = 0,5 см2/с2, тхм = ту = -0,5 см2/с2 - заданные значения составляющих т на северной, восточной, южной и западной границах области (0 < х < L, 0 < у < B), в которую вписано Черное море.
В начальный момент времени толщина верхнего слоя к0 равнялась 175 м, или глубине моря Н при Н < 175 м. Коэффициенты придонного трения г, Г
на основании предыдущих расчетов были выбраны равными 0,001 и 0,002 см/с соответственно, коэффициент бигармонической вязкости Ав = 41016 см4/с4 Использовалась прямоугольная сетка с шагами Ах = Ду = 3 км, шаг по времени Д( = 1,5 мин.
Продолжительность расчетов в обоих экспериментах составила 20 лет (один год состоит из 12 мес по 30 сут). За это время оба решения вышли на квазипериодический статистически-равновесный режим. При этом средние по площади характеристики, а именно доступная потенциальная энергия ВРЕ, кинетическая энергия в слоях КЕ1, КЕ2, работа сил тангенциального напряжения ветра Ж, придонного трения Жв, горизонтальной вязкости Жав
БРЕ = (Р1 g / 2) - Р %о2 / 2), КЕ = (М К2 + V2) / 2), КЕ2 =(р2И2(ы2 + v22)/2), Жт=(р1(и1тх + vlтy)),
&кв = 12 (Рк (Г + г2д/и1 + VI )(и1 + VIЛ,
и—л » '
Жав = § ( ркАш (ик V 2(Нк V \) + Ук V2 (Нк V \ )\
к=1
изменялись лишь в ограниченных пределах относительно некоторых средних значений (рис. 1). Угловые скобки означают осреднение по площади.
Рис. 1. Графики БРЕ, КЕ\, КЕг в экспериментах R1 (а) и (в); графики Wг, Шлв, Wдв, 3 в экспериментах R1 (б) и К2 (г)
Энергия и работа сил удовлетворяют уравнению энергетического баланса
ЯЕ
= Жт+ Жкв + Жав +$ ,
от
где Е = БРЕ + КЕХ + КЕ2 ; д - невязка, обусловленная точностью вычисления составляющих баланса энергии. Как видно на рис. 1, б, г, в обоих экспериментах R1 и в течение всего расчетного периода с хорошей точностью выполняется энергетический баланс - д практически равняется нулю.
Эксперимент К1, р > 0. В этом эксперименте в течение первого года ВРЕ и КЕ1 достигают значений, в пределах которых в дальнейшем и происходят их колебания (рис. 1, а). Причем амплитуда этих колебаний увеличивается после 10-го года. Энергия КЕ2 постепенно растет в течение первых 10 лет, затем, после выхода решения на статистически-равновесный режим, также как ВРЕ и КЕ1 , колеблется относительно некоторого среднего значения.
Работа Жт (рис. 1, б) изменяется подобно КЕ1 - растет в течение первого года, затем колеблется относительно среднего значения, характерного для статистически-равновесного режима. Характеристики Жав и Жив увеличиваются с ростом КЕ2 и уменьшаются при ее падении.
Обратим внимание на наличие на графиках рис. 1, а, б колебаний двух видов. Хорошо видны (особенно на графике Ж) осцилляции с периодом ~50 сут, которые связаны с процессами гидродинамической неустойчивости течений в верхнем слое. Эти более высокочастотные колебания накладываются на колебания с периодами от полугода до двух лет. Можно предположить, что такие долгопериодные колебания обусловлены изменчивостью интегральной циркуляции.
Далее перейдем к анализу пространственных полей, полученных в результате расчетов. Вначале рассмотрим осредненные за последние 15 лет поля М, и1, и2, и1, и2, щ (рис. 2).
В средних полях и1 и и (рис. 2, а, в) наблюдается струйное циклоническое течение, расположенное над материковым склоном (аналог Основного Черноморского течения - ОЧТ [6]), скорости в котором достигают 80 см/с.
Поле Ь имеет хорошее соответствие с полями течений и1 и потоков и1. В стрежне струи градиенты \ максимальны, слева от стрежня к центру бассейна \ уменьшается, правее по направлению к берегу - увеличивается. Западнее Крыма и в юго-восточной части моря находятся локальные области с большими значениями \ - Севастопольский и Батумский антициклоны [7].
Средние течения в нижнем слое также имеют циклонический характер, но, в отличие от течений в верхнем слое, они привязаны к особенностям рельефа дна и направлены преимущественно вдоль изобат Н (рис. 2, з). Наибольшие скорости течений (до 17 см/с) наблюдаются на материковом склоне в местах резкого свала глубин и вблизи мест пересечения поверхности раздела с дном, где толщина нижнего слоя мала. Несмотря на то что скорости течений в нижнем слое существенно меньше, чем в верхнем, потоки в нижнем слое, наоборот, превышают потоки в верхнем слое. Причиной этому является большая разница в толщинах слоев. Это приводит к тому, что интегральная функция тока щ (рис. 2, е) в большей мере отражает потоки в нижнем слое, особенно в глубоководной части моря, где толщина нижнего слоя велика.
О 200 400 600 800 1000Л', кмО 200 400 600 800 1000 X, км
Рис. 2. Осредненные за 15 лет поля и1 (а), и 2 (б), и1 (в), и 2 (г), ^ (д), р (е) в эксперименте R1; модуль средних потоков в верхнем слое и изолиния ^ = 200 м с контрольными точками
(ж); рельеф дна H(x, у) (з). На картах и1, и2, и1, и2 под стрелками указаны максимальные значения, цветовые шкалы соответствуют абсолютным значениям
Перейдем к рассмотрению мгновенных полей, полученных в эксперименте Я1. На рис. 3, а - е поля М, и, и2, и, и2 , V приводятся для одного момента времени. Они являются достаточно типичными для статистически-равновесного режима. В верхнем слое в полях и, и отчетливо проявляется струйное меандрирующее течение, расположенное над материковым склоном вдоль всего периметра бассейна. Из анализа результатов следует, что меандры распространяются в направлении струи в виде волн с фазовой скоростью, которая меньше скорости течения. Ширина струи составляет 30 - 50 км, ско-
рость в стрежне 80 - 120 см/с. Справа и слева от потока наблюдаются антициклонические и циклонические вихри соответственно. Их образование является следствием меандрирования основного течения, а также обтекания особенностей береговой линии. Кроме того, слева от струи периодически происходит формирование антициклонических рингов в результате отрыва крупных меандров. На рис. 3, а, в можно видеть образование такого ринга в северо-восточной части бассейна.
0 200 400 600 600 1000Л", кмО 200 400 600 800 ЮООЛ", км
Рис. 3. Мгновенные поля и (а), и2 (б), и (®), и2 (г), \ (д), щ(е) в эксперименте R1 для момента времени 15.04.0008 г.; поле ^ для момента времени 20.04.0008 г. (ж); поле ^ для момента
времени 25.04.0008 г. (з). На картах щ , и2, и1, и2 под стрелками указаны максимальные значения, цветовые шкалы соответствуют абсолютным значениям
Как уже отмечалось, с полем течений и потоков в верхнем слое тесно связано поле М. В антициклонах поверхность раздела слоев заглубляется (толщина верхнего слоя увеличивается), в циклонах поднимается (М уменьшается), градиенты М обостряются в области струйных течений. Благодаря этому все особенности циркуляции верхнего слоя (струйные течения, меандры, вихри) отчетливо видны в поле М (рис. 3, д). По последовательным распределениям М с интервалом 5 сут (рис. 3, д, ж, з) можно проследить динамику верхнего слоя, перемещение меандров вдоль основного потока и образование ринга.
В нижнем слое, при преимущественно циклонической направленности течений и потоков, картина циркуляции иная, чем в верхнем слое. В ней отсутствует круговое струйное течение. Наблюдаемые потоки имеют нерегулярный характер, представляют собой как бы отдельные струи, вытянутые вдоль изобат и к тому же перемещающиеся вдоль материкового склона. Это происходит в результате наложения интенсивных волн и/или вихрей на средний поток.
Скорости течений в нижнем слое в большей его части не превышают 5 см/с, в струях они составляют ~10 - 15 см/с, а вблизи мест пересечения поверхности раздела с дном, где толщина нижнего слоя h2 мала, они могут достигать 35 см/с. Возникновение нереально больших скоростей в тонких слоях связано, по-видимому, с не совсем корректным учетом придонного трения. Таким течениям не стоит придавать большого значения, поскольку в потоках подобные особенности в тонких слоях не проявляются (рис. 3, г).
Вследствие большой разницы в толщинах слоев потоки в нижнем слое оказываются по величине больше, чем в верхнем (рис. 3, в, г). Это приводит к тому, что интегральная функция тока щ (рис. 3, е) в большей степени отражает характер циркуляции нижнего слоя, особенно в глубоководной части моря.
Чтобы идентифицировать волновые процессы в области материкового склона, были построены временные диаграммы ) = \(?) -\ и у'(X) = у(?) -у вдоль изолинии ^ = 200 м, практически совпадающей со стрежнем осредненного
потока в верхнем слое (рис. 2, ж) и расположенной в области свала глубин примерно над изобатой 1600 м. На этой линии выбраны контрольные точки с интервалом 120 км. На временных диаграммах (рис. 4) им соответствуют вертикальные пунктирные линии, номера точек указаны на верхней границе. На обеих диаграммах хорошо видно, что вдоль материкового склона в направлении против часовой стрелки распространяются возмущения в виде волн. Скорости распространения этих волн отличаются для разных участков траектории. Так, вдоль южного побережья волны движутся быстрее, чем вдоль северного. Вероятнее всего, это связано с различной крутизной материкового склона. На участке 15 - 13 - 10 (рис. 4) кроме отмеченных волн наблюдаются более медленные возмущения, связанные с образованием в районе Анатолийского побережья антициклонических вихрей.
Существует большое внешнее сходство между двумя приведенными диаграммами, что указывает на тесную взаимосвязь волновых процессов в различных слоях. При этом диаграмма И[ (() отражает волновые процессы, происходящие в верхнем слое, а диаграмма у' ) демонстрирует динамику интегральной циркуляции, больший вклад в которую вносит нижний слой. Отме-
тим, что на диаграмме у'(?) на фоне более длинных волн видны пульсации, связанные с топографическими эффектами.
Рис. 4. Временные диаграммы к[ (?), м, (слева) и у'(?), х 106 м3/с, (справа) вдоль изолинии \ = 200 м в течение 10-го года
По наклону фазовых линий (рис. 4) можно рассчитать скорость распространения волн, которая составляет ~19 см/с на участке 1 - 15 - 10 и ~16 см/с на участке 10 - 5 - 1. Полученные величины фазовых скоростей сравнимы со средними значениями скорости течений в нижнем слое, но они гораздо меньше средних скоростей течений в верхнем слое в области ОЧТ.
Для определения частотных характеристик волн временные ряды hi(t) и y(t) в контрольных точках были исследованы методом спектрального анализа [8]. На построенных графиках спектральной плотности SA (у), Sw (у) для
всех рассмотренных точек выделяются значимые пики на частотах vi « 0,021 и V2 « 0,036 1/сут, что указывает на дискретный характер наблюдаемых волн. Эти частоты соответствуют периодам колебаний Ti « 48 и T2 « 27 сут. Полученные результаты согласуются с аналитическими оценками, приведенными в работах [9, 10].
Кроме волн, распространяющихся над материковым склоном, в центральной глубоководной части моря в поле y(t) были выделены волны Рос-сби, перемещающиеся в западном направлении. На рис. 5 приведена временная диаграмма y(t ), построенная для разреза Y = 250 км. Видно, что в течение двух лет (10 - 11-й годы) такая волна четко проявлялась два раза, в апреле -
мае 10-го года и в октябре - ноябре 11-го года. Фазовая скорость волны, рассчитанная по наклону фазовых линий, составляет ~7,5 см/с и соответствует фазовой скорости первой нормальной моды баротропной волны Россби для замкнутого бассейна с горизонтальным дном [11 - 13]. Распространение этой волны показано на последовательных (с интервалом 20 сут) распределениях интегральной функции тока 30 сентября - 30 ноября 11-го года.
I, мес
^^ТТ............^^ТТЛ .....1.....
-50 -40 -30 -20 -10 0 Х106м3/С -50 -40 -30 -20 -10 0 х106м3/с
Рис. 5. Временные диаграммы ¡^у), х 106 м3/с, для разреза У = 250 км с 10-го по 11-й год
в эксперименте R1 (слева) и последовательные с интервалом 20 сут поля щ, х 106 м3/с, (справа). Разрез показан штриховой линией
Эксперимент Я2, р = 0. При отсутствии в-эффекта процесс формирования циркуляции в течение первых 5 лет подобен процессу в эксперименте R1. В обоих слоях в области материкового склона образуются и развиваются течения, распространяющиеся по всему периметру Черного моря. На эти тече-
ния накладывались волновые возмущения, связанные с обтеканием рельефа дна (топографические волны Россби [9, 14, 15]). Графики доступной потенциальной энергии, кинетической энергии, работы сил тангенциального напряжения ветра, придонного трения и горизонтальной вязкости в эксперименте И2 (рис. 1, в, г) в первые 5 лет мало отличаются от аналогичных графиков в эксперименте Я1 (рис. 1, а, б). Также в течение этого периода пространственные распределения течений, толщины верхнего слоя и интегральной функции тока качественно и даже количественно близки к распределениям, полученным с учетом ^-эффекта.
О 200 400 600 800 ЮООЛГ, кмО 200 400 600 800 1000 X, км
Рис. 6. Мгновенные (на 30.12.0020 г.) и осредненные за 15 лет поля в эксперименте К2: и (а);
и (б); и 2 (в); и 2 (г); \ (д); ^ (е); щ(ж); у (з). На картах и , и , и 2, и 2 под стрелками указаны максимальные значения, цветовые шкалы соответствуют абсолютным значениям
Начиная с 6-го года картина циркуляции в эксперименте R2 меняется и принимает вид, показанный на рис. 6. В отличие от эксперимента R1 основной циклонический круговорот не распространяется на весь бассейн, а сжимается к его центральной части. Это видно как по мгновенным (рис. 6,
а, в, д, ж), так и по осредненным (рис. 6, б, г, е, з) полям.
Заключение. Сравнивая результаты экспериментов R1 и R2, можно сделать следующие выводы.
Длительное воздействие циклонического ветра при учете реального рельефа дна и ^-эффекта приводит к формированию в бассейне меандрирующе-го струйного циклонического течения (аналог ОЧТ) со стрежнем, расположенным в области материкового склона. На это течение накладываются длинноволновые колебания (топографические волны Россби), причем как в верхнем, так и в нижнем слое. По своим параметрам и характеру проявления эти колебания соответствуют волнам, захваченным материковым склоном и распространяющимся в циклоническом направлении.
Изменение параметра Кориолиса с широтой (S-ффект) в Черном море, кроме влияния на топографические волны Россби, проявляется также в образовании в центральной глубоководной части моря баротропных волн Россби для замкнутого бассейна, распространяющихся с востока на запад. Отсутствие этих волн приводит к тому, что круговое течение не распространяется на все море, а локализуется в его центральной части.
Работа выполнена в рамках государственного задания по теме № 08272014-0011 «Исследования закономерностей изменений состояния морской среды на основе оперативных наблюдений и данных системы диагноза, прогноза и реанализа состояния морских акваторий» (шифр «Оперативная океанография»).
СПИСОК ЛИТЕРАТУРЫ
1. Павлушин А.А., Шапиро Н.Б., Михайлова Э.Н., Коротаев Г.К. Двухслойная вихрераз-решающая модель ветровых течений в Черном море // Морской гидрофизический журнал. - 2015. - № 5. - С. 3 - 22. - doi:10.22449/0233-7584-2015-5-3-22
2. Павлушин А.А., Шапиро Н.Б., Михайлова Э.Н. Проявление ^-эффекта в двухслойной модели Черного моря // Морской гидрофизический журнал. - 2016. - № 5. - С. 3 - 24. -doi:10.22449/0233-7584-2016-5-3-24
3. 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 // J. Phys. Oceanogr. - 1975. - 5, Iss. 4. - P. 642 - 657. - doi: 10.1175/1520-0485(1975)005<0642:0TG0ME>2.0.C0;2
4. Мезингер Ф., Аракава А. Численные методы, используемые в атмосферных моделях. -Л.: Гидрометеоиздат, 1979. - 136 с.
5. Каменкович В.М., Кошляков М.Н., Монин А.С. Синоптические вихри в океане. - Л.: Гидрометеоиздат, 1987. - 512 с.
б. Блатов А.С., Булгаков Н.П., Иванов В.А. и др. Изменчивость гидрофизических полей Черного моря. - Л.: Гидрометеоиздат, 1984. - 239 с.
7. Иванов В.А., Белокопытов В.Н. Океанография Черного моря. - Севастополь: МГИ НАН Украины, 2011. - 212 c.
8. Нуссбаумер Г. Быстрое преобразование Фурье и алгоритмы вычисления сверток. - М.: Радио и связь, 1985. - 248 с.
9. Ефимов В.В., Куликов Е.А., Рабинович А.Б. и др. Волны в пограничных областях океана. - Л.: Гидрометеоиздат, 1985. - 280 c. - URL: http://elib.rshu.ru/files_books/pdf/img-503182759.pdf (дата обращения 17.05.2017).
10. Иванов В.А., Янковский А.Е. Длинноволновые движения в Черном море. - Киев: Науко-ва думка, 1992. - 110 с.
11. Rachev N.H., Stanev E. V. Eddy Processes in Semienclosed Seas: A Case Study for the Black Sea // J. Phys. Oceanogr. - 1997. - 27. - P. 1581 - 1601. - doi:10.1175/1520-0485(1997)027<1581:EPISSA>2.0.C0;2
12. Stanev E. V., Rachev N.H. Numerical study on the planetary Rossby modes in the Black Sea // J. Mar. Sys. - 1999. - 21, Iss. 1 - 4. - P. 283 - 306.
13. Longuet-Higgins M.S. Planetary Waves on a Rotating Sphere // Proc. Roy. Soc. London. Ser. A. -1964. - 279, No. 1379. - P. 446 - 473. - URL: http://www.jstor.org/stable/2414690?origin=JST0R-pdf&seq=1#page_scan_tab_contents (дата обращения 26.05.2017).
14. Физика океана. Т. 2. Гидродинамика океана / Отв. ред. В.М. Каменкович, А.С. Мо-нин. - М.: Наука, 1978. - 456 с.
15. Монин А.С., Жихарев Г.М. Океанские вихри // Успехи физических наук. - 1990. - 160, вып. 5. - P. 1 - 46.
The role of the bottom relief and the ^-effect in the Black Sea
dynamics
A. A. Pavlushin*, N. B. Shapiro, E. N. Mikhailova
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia *e-mail:pavlushin@mhi-ras.ru
The results of the numerical experiments carried out within the two-layer eddy-resolving Black Sea model are discussed. Bottom relief, ^-effect, bottom friction and horizontal turbulent viscosity parameterized by the bi-harmonic operator are taken into account in the model. The calculations cover a long period up to the moment when the solution achieves the statistically equilibrium mode which is characterized by presence of intense currents, waves and vortices. It is shown that under the statistically equilibrium mode, a cyclonic circulation is formed in the sea: in the upper layer - a meandering flow (the Rim Current analog); in the lower layer - rather intensive waves which are imposed on the flow propagating along the isobaths. These waves can be characterized as the topographic Rossby waves trapped by the continental slope. The technique for analyzing such waves is proposed. It is shown that the currents' wave disturbances in the lower layer significantly influence the flows in the upper layer contributing to their instability and meandering.
Keywords: Black Sea, eddy resolving model, numerical experiment, ^-effect, topographic Rossby waves, bottom relief, trapped waves.
