Численный прогностический расчет течений в Черном море с высоким горизонтальным разрешением Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

Математическое моделирование

морских систем

УДК 551.465 С.Г. Демышев

Численный прогностический расчет течений в Черном море с высоким горизонтальным разрешением

Представлены результаты численного прогностического эксперимента с разрешением 1,64 км по горизонтальным координатам. В отличие от расчетов с более грубым разрешением получены следующие динамические особенности. Основное Черноморское течение как непрерывная струя, направленная вдоль свала глубин, прослеживается до 400 м. Вдоль восточной части Анатолийского побережья регулярно образуются мезомасштабные антициклонические вихри, которые влияют на формирование Батумского антициклона. Севастопольский, Синоп-ский, Кизилирмакский и Кавказский антициклонические вихри являются квазипериодическими. Между вновь образованным Севастопольским антициклоном и предыдущим вихрем формируется область циклонического вращения вод.

Ключевые слова: Черное море, численная модель, синоптические и мезомасштабные вихри.

Введение

К настоящему времени общепринято, что циркуляция в Черном море представляет собой обширный циклонический круговорот, по периферии которого периодически возникают и эволюционируют синоптические (с диаметром от нескольких десятков до сотен километров) и мезомасштабные (от нескольких до десятков километров) вихри (см., например, [1]). Синоптические вихри могут быть квазистационарными (Батумский антициклон) и квазипериодическими (Севастопольский, Синопский и др.). По спутниковым данным поверхностной температуры неоднократно фиксировались вихри с радиусом порядка десяти километров. Они наблюдались в области Анатолийского и Кавказского побережий. Насколько такие вихри важны в обмене энергией и перераспределении импульса, тепла и соли между движениями различных масштабов?

Для ответа на этот вопрос, помимо натурных исследований, необходимы численные эксперименты с пространственным разрешением, позволяющим аккуратно аппроксимировать мезомасштабные структуры - вихри и течения. В данной работе приводятся первые результаты прогностического расчета с разрешением 1,64 км по горизонтальным координатам. Получены особенности в поведении Основного Черноморского течения (ОЧТ), Батумского, Севастопольского антициклонов и прибрежных мезомасштабных вихрей.

© С.Г. Демышев, 2011 36

ШБЫ 0233-7584. Мор. гидрофиз. журн., 2011, № 1

Уравнения модели

Система уравнений модели в приближении Буссинеска, гидростатики и несжимаемости морской воды [2] имеет вид

ыг - (Х + /> + ^ыг =-gp0gx - —(Р1 + Е)х + (ьуыг )г-инУ4ы , (1)

А

0

V, + (X + f)u + wv, = -gPoVy - — (P1 + E)y + (u Vz) z - vH V4v , (2)

А

0

2

Р = gpZ + g I рЛт= gPoZ+Р1, (3)

0

Ых + уу + = 0, (4)

н

V + I(Ых + Уу)й2 = (Рг -Еу)/ р, (5)

0

Т+(ЫТ) х + (УТ ) у + (*Т) 2 = -к V 4Т + (кТ) 2, (6)

5,+(Ы5) х + (У^)^ + 2 =-к V4,? + (кЧ) 2, (7)

р= <Р(Т, (8)

где X =----+1, Е = р0-. Остальные обозначения общепринятые

Эх Эу 2

[2]. Заметим, что уравнение (5) получено в предположении выполнения линеаризованного кинематического условия

Рг - ЕУ

w = -v +-

Pi

Для системы уравнений (1) - (8) краевые условия следующие. На поверхности z = 0

Ev — Pr

vvuz = -tx, nvVz = -ty, kTz = QT, kvSz =-S0 + ß(Scl - S0), (9)

P

где tx, ty - касательное напряжение трения ветра; QT - поток тепла; Ev - испарение морской воды; Pr - осадки; S0 - поверхностная соленость; Scl - климатическая соленость; р1 - плотность морской воды; ß - параметр релаксации.

На дне z = H(x, y)

u = v = w = 0, Tz = 0, Sz = 0. (10)

На твердых боковых стенках:

ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 1 37

для меридиональных участков границы

и = 0, V 2и = 0, V = 0, V V = 0,

х

(11)

Тх = 0, (V Т) х = 0, ^х = 0, ^) х = 0,

для зональных участков границы

V = 0, V2v = 0, uv = 0, V X = 0,

* * (12) Ту = 0, ^Т) у = 0, Я = 0, (V 2£) у = 0.

На участках границы, где вода втекает, используются условия Дирихле: для меридиональных участков

и = ир, V2u = 0, V = 0, V2v = 0,

х

Т = Тр, £ = £р, (У2Г)х = 0, (V2S)х = 0,

(13)

для зональных участков

V = ур, V2v = 0, иу = 0, V2uv = 0,

* * (14) Т = Тр, £ = £р, ^2Т)у = 0, (V2 £ )у = 0.

Для верхнебосфорского течения и Керченского пролива, когда течение направлено из Черного моря в Азовское,

V = у', V2v = 0, иу = 0, V2uv = 0,

* * (15) Ту = 0, = 0, (V2Т)* = 0, (V2£)* = 0.

В качестве начальных при ? = ?0 потребуем выполнение следующих условий:

и = и0(х, г), V = v0(х, г), д = д°(х, *), Т = Т0(х, г), £ = £0(х, г). (16)

Система уравнений (1) - (8) с соответствующими краевыми (9) - (15) и начальными (16) условиями решалась численно.

Условия проведения эксперимента

За начальные поля принимались климатические температура, соленость, уровень и скорость течений, полученные ранее в работе [3]. Начальный момент соответствовал 1 ноября.

Для расчета на поверхности моря использовалось касательное напряжение трения ветра на каждые сутки года, полученное по распределениям приземного давления [4]. Потоки тепла на поверхности моря задавались в соответствии с работой [5], поля осадков и испарений - согласно [6]. Сток при-

38

0233-7584. Мор. гидрофиз. журн., 2011, № 1

черноморских рек (Дунай, Днестр, Днепр, Южный Буг, Риони, более мелкие кавказские реки и реки Турции - Ешиль-Ирмак, Кизил-Ирмак, Сакарья) и годовой ход температуры воды в устьях рек, кроме рек Турции, рассчитывался по [7]. В силу отсутствия данных по рекам Турции предполагалось, что температура воды в них равна температуре прибрежных вод моря.

В нижнебосфорском течении соленость принималась равной 35%о, температура 16°С, что соответствует характеристикам вод Мраморного моря.

Расчеты проводились на горизонтальной сетке 1,64 х 1,64 км, по вертикали использовалось 27 горизонтов с глубинами 2,5, 5, 7,5, 10, 12,5, 15, 17,5, 20, 22,5, 25, 30, 35, 40, 55, 75, 100, 125, 150, 200, 250, 300, 500, 700, 1000, 1400, 1800, 2100 м. Шаг по времени равнялся 1,5 мин.

Коэффициенты горизонтального турбулентного обмена импульсом и горизонтальной турбулентной диффузии тепла и соли принимались равными

пн = 1016ш4 • й-1, к = 5- 1016ш4 • й-1.

Коэффициенты турбулентности по вертикали вычислялись в соответствии с параметризацией Филандера - Пакановского [8, 9].

Интегрирование уравнений модели (1) - (8) с учетом краевых (9) - (15) и начальных (16) условий проводилось на полтора года модельного времени.

Результаты численного эксперимента

Приведенный уровень моря характеризует пространственную структуру и временную эволюцию течений в верхнем слое для средних широт, где геострофический баланс выполняется довольно точно. В расчете четко проявилось несколько особенностей, которые в экспериментах с разрешением 5 км по горизонтали либо не наблюдались, либо проявлялись в существенно сглаженном виде.

Сопоставление рис. 1 и 2 показывает, что с улучшением разрешения ОЧТ принимает вид узкого струйного течения, направленного строго вдоль свала глубин. Из анализа расчета следует, что градиент поперек струи в поле уровня увеличился примерно в три раза.

В поле уровня наблюдаются вихри диаметром в несколько десятков километров. Наиболее интенсивно такие вихри образуются вдоль Кавказского и в восточной части Анатолийского побережья (рис. 2, а, б). Они формируются квазипериодически и передвигаются вдоль свала глубин по направлению ОЧТ. Анализ их изменчивости показал, что основной причиной образования мезомасштабных вихрей между берегом и периферией ОЧТ является влияние на ОЧТ особенностей орографии берега Черного моря. Такие вихри неоднократно фиксировались на спутниковых снимках поверхностной температуры [10, 11].

0233-7584. Мор. гидрофиз. журн., 2011, № 1

39

с.ш. —I-1-1-1-1-1-1-1-1-г

46° 45° 44°

43°-

42°-

41°Ь

28° 30° 32° 34° 36° 38° 40° в.д.

Р и с. 1. Уровенная поверхность (см) для начального момента времени

44° 43° 42° 41°

29° 32° 35° 38° в.д. 29° 32° 35° 38° в г

в.д.

Р и с. 2. Рассчитанная уровенная поверхность (см) для различных моментов времени

с.ш

с.ш

о

о

о

о

40

0233-7584. Мор. гидрофиз. журн, 2011, № 1

Вихри, которые передвигаются вдоль Кавказского побережья, влияют на формирование и эволюцию Керченского антициклона. Его спонтанное образование объясняется нерегулярным поведением и меняющейся интенсивностью мезомасштабных вихрей (рис. 3, а), образующихся в различные сезоны.

Вдоль Анатолийского побережья мезомасштабные круговороты формируются регулярно. Они передвигаются на восток и не просто взаимодействуют с Батумским антициклоном, но и оказывают существенное влияние на его формирование. Ранее предполагалось, что Батумский антициклон образуется в результате действия ветра, который в Кавказском регионе имеет положительную завихренность. Анализ результатов эксперимента показал, что в отсутствие положительной завихренности в поле ветра образование Батумского круговорота в значительной мере обусловлено именно мезомасштабными антициклонами (рис. 3, б). Важная особенность, проявляющаяся на всем периоде интегрирования, - сопутствующий Батумскому антициклону мезомас-штабный циклонический вихрь, который формируется по его периферии. На рис. 2, в, г он расположен к западу от центра антициклона. В качестве примера на рис. 4 приведено поле скорости в районе Батумского антициклона. Зимой циклонический вихрь наблюдается на западной периферии антициклона (рис. 4, а), весной - на восточной (рис. 4, б). Анализ спутниковых данных показал замедленное движение вод по восточной и западной периферии Батум-ского антициклона по сравнению с его северной и южной частями, что косвенно свидетельствует о наличии циклона [11].

В районе Крымского побережья периодически формируется Севастопольский антициклон (рис. 2), который скатывается вниз по течению, и на его месте вновь образуется новый вихрь. Как показывает анализ поля уровня и течений, он также формируется за счет влияния на ОЧТ орографии берега. О том, что Севастопольский антициклон не стоит на месте, свидетельствуют данные наблюдений [12]. В проведенном эксперименте получено, что в период формирования нового антициклона между «старым» и новым вихрем образуется область циклонической завихренности. В этом состоит отличие от ранее проведенных расчетов. В поле скорости (рис. 5) указанная изменчивость течений подтверждается траекторией дрейфующих буев в районе Севастопольского антициклона: два из пяти дрифтеров описали кривую, которая соответствует циклоническому вращению вод [11, 12].

Анализ временной изменчивости показал, что Севастопольский, Синоп-ский, Кизилирмакский и Кавказский антициклонические вихри не являются квазистационарными. Они формируются и смещаются вниз по течению, образуя цепочку антициклонических вихрей между ОЧТ и шельфом Черного моря.

0233-7584. Мор. гидрофиз. журн., 2011, № 1

41

а

б

Р и с. 3. Поле уровня (см) в районе Кавказского (а) и Анатолийского (б) побережий 15 мая

42 0233-7584. Мор. гидрофиз. журн, 2011, № 1

44° с.ш.

-'-Г

40°

1-1-1-1-1-Г"

41° в.д.

с.ш.

1-1-'-11° в.д.

Р и с. 4. Поле скорости в районе Батумского антициклона на глубине 2,5 м 20 декабря (а) и 1 апреля (б)

0233-7584. Мор. гидрофиз. журн., 2011, № 1

43

О

а

О

О

б

44°

45°

46°

с.ш

29°

30°

31°

32°

33° вд34°

Р и с. 5. Поле скорости на северо-западном шельфе 20 декабря на глубине 2,5 м

В отличие от предыдущих расчетов с горизонтальным разрешением порядка 5 км в данном эксперименте ОЧТ в виде непрерывной струи, направленной вдоль свала глубин, наблюдается до 400 м. Внутри циклонического круговорота формируются вихри и струи синоптического и мезомасштабного характера. В верхнем 400-метровом слое западной части циклонического круговорота образуются в зависимости от сезона два-три синоптических вихря и струйное течение шириной ~ 30 км (рис. 6, а). В его восточной части наблюдаются несколько вихрей диаметрами 20 - 50 км (рис. 6, б) с различным знаком завихренности. Ниже 500 м поле течений имеет сложную структуру, которая характеризуется наличием вихрей и узких струйных течений.

Уровень моря может быть представлен в виде д = да + дй, где да - ади-

характеризует придонные геострофические течения. В анализируемом расчете на протяжении всего периода интегрирования придонный слой охвачен циклоническим круговоротом, периферия которого строго соответствует свалу глубин. Во внутренней области этого круговорота возникают и эволюционируют вихревые образования (рис. 7) различного знака завихренности. Их изменчивость по времени составляет порядка нескольких десятков суток. Причем на востоке бассейна процессы вихреобразования и эволюции выражены заметно ярче, чем на западе.

44 ISSN 0233-7584. Мор. гидрофиз. журн., 2011, № 1

намическая

[13]. Нетрудно убедиться, что V

а

44.0° с.ш.

43.5°-

43.0°-^

42.5°-',

42.0°

29°

30°

"1--1-1-1-'

31° 32° 33°

34° в.д. 35°

с.ш.

44.0°-;

43.5°-

43.0°-

42.5°-

42.0°-

41.5°

34°

35°

36°

37°

38°

39° в.д. 40°

б

Р и с. 6. Поле скорости на глубине 400 м в западной части циклонического круговорота 20 декабря (а) и в его восточной части 1 апреля (б)

а

0233-7584. Мор. гидрофиз. журн., 2011, № 1

45

с.ш

47

46'

45'

44'

43'

42'

41'

28°

30°

32°

34°

36°

38°

40°

в.д.

Р и с. 7. Поле адинамической поправки (см) 20 декабря

Заключение

Проведенный эксперимент показал, что улучшение горизонтального разрешения до 1-2 км при моделировании течений в Черном море может иметь принципиальное значение. Основные динамические эффекты, которые проявились в расчете, следующие.

ОЧТ как непрерывная струя, направленная вдоль свала глубин, прослеживается до 400 м. Скорость в ней достигает 90 см-с-1, что соответствует наблюдениям. Между ОЧТ и берегом образуются и эволюционируют синоптические и мезомасштабные вихри. Вдоль восточной части Анатолийского побережья в результате влияния береговой орографии непрерывно формируются мезомас-штабные антициклонические вихри. Они передвигаются вдоль берега, достигают юго-восточной части моря и образуют Батумский антициклон. На его восточной и западной периферии периодически возникает мезомасштабный циклонический вихрь, что соответствует спутниковым наблюдениям.

Следующий важный элемент циркуляции, который ранее не воспроизводился, касается динамики вод на северо-западном шельфе. В процессе квазипериодического формирования Севастопольского антициклона и смещения его вниз по течению формируется область циклонического вращения вод, что также подтверждается данными наблюдений.

Расчет продемонстрировал, что Севастопольский, Босфорский, Синоп-ский, Кизилирмакский и Кавказский антициклонические вихри являются квазипериодическими.

Необходимо отметить, что указанные особенности могут быть связаны с перестройкой циркуляции, которая обусловлена приспособлением поля плотности к рельефу дна и орографии берега. Этот процесс длительный. Например, первая бароклинная мода (в зависимости от длины волны) имеет декремент затухания до трех лет [14]. Для подтверждения полученных результатов требуются дальнейшие исследования.

СПИСОК ЛИТЕРАТУРЫ

1. Блатов А.С., Булгаков Н.П., Иванов В.А. и др. Изменчивость гидрофизических полей в Черном море. - Л.: Гидрометеоиздат, 1984. - 271 с.

2. Демышев С.Г., Коротаев Г.К. Численная энергосбалансированная модель бароклинных течений океана на сетке С // Численные модели и результаты калибровочных расчетов течений в Атлантическом океане. -М.: ИВМ РАН, 1992. - С. 163 -231.

46

0233-7584. Мор. гидрофиз. журн, 2011, № 1

3. Демышев С.Г., Иванов В.А., Маркова Н.В., Черкесов Л.В. Построение поля течений в Черном море на основе вихреразрешающей модели с ассимиляцией климатических полей температуры и солености // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофи-зика, 2007. - Вып. 15. - C. 215 - 226.

4. Дорофеев В.Л., Коротаев Г.К. Ассимиляция данных спутниковой альтиметрии в вихреразрешающей модели циркуляции Черного моря // Морской гидрофизический журнал. - 2004. - № 1. - С. 52 - 68.

5. Staneva J.V., Stanev E.V. Oceanic response to atmospheric forcing derived from different climatic data sets. Intercomparison study for the Black Sea // Oceanologia Acta. - 1998. - 21, № 3. - Р. 383 - 417.

6. Ефимов В.В., Тимофеев Н.А. Теплобалансовые исследования Черного и Азовского морей. - Обнинск: ВНИИГМИ-МЦД, 1990. - 237 с.

7. Гидрометеорология и гидрохимия морей СССР. Т. IV. Черное море. Выпуск 1. Гидрометеорологические условия. - СПб.: Гидрометеоиздат, 1991. - 428 с.

8. Pacanowski R.C., Philander S.G.H. Parameterization of vertical mixing in numerical models of tropical oceans // J. Phys. Oceanogr. - 1981. - 11, № 3. - P. 1443 - 1451.

9. Демышев С.Г., Коротаев Г.К., Кныш В.В. Моделирование сезонной изменчивости температурного режима деятельного слоя Черного моря // Изв. РАН. ФАО. - 2004. - 40, № 2. - С. 259 - 270.

10. Sur H.I., Ozsoy E., Ilyin Y.P., Unluata U. Coastal/deep ocean interactions in the Black Sea and their ecological/environmental impacts // J. Mar. Syst. - 1996. - 7. - Р. 293 - 320.

11. Zatsepin A.G., Ginzburg A.I., Kostianoy A.G. et al. Observation of Black Sea mesoscale eddies and associated horizontal mixing // J. Geophys. Res. - 2003. - 108, № C8. - Р. 1 - 27.

12. Ginzburg A.I., Kostianoy A.G., Nezlin N.P. et al. Anticyclonic eddies in the northwestern Black Sea // J. Mar. Syst. - 2002. - 32. - Р. 91 - 106.

13. Саркисян А.С., Передерей А.И. Динамический метод как первое приближение при расчете уровенной поверхности бароклинного океана // Метеорология и гидрология. -1974. - № 3. - С. 45 - 54.

14. Демышев С.Г., Коротаев Г.К., Куфтарков А.Ю. Анализ начального периода приспособления геофизических полей Черного моря к новой численной модели // Изв. РАН. ФАО. - 1996. - 32, № 5. - С. 635 - 644.

Морской гидрофизический институт НАН Украины, Материал поступил

Севастополь в редакцию 28.09.09

E-mail: sgdem@stel.sebastopol.ua После доработки 25.11.09

АНОТАЦ1Я Представлен результати чисельного прогностичного експерименту з просторо-вою роздшьною здатшстю 1,64 км по горизонтальних координатах. На вщмшу вщ розрахунюв з бшьш грубою просторовою роздшьною здатшстю отримаш наступш динамiчнi особливоси. Основна Чорноморська течiя як безперервний струмшь, направлений уздовж звалювання гли-бин, простежуеться до 400 м. Уздовж схщно! частини Анатолшського побережжя регулярно утворюються мезомасштабш антициклонш вихори, яю впливають на формування Батумського антициклону. Севастопольський, Сшопський, Юзшмарський i Кавказький антициклонш вихори е квазшерюдичними. Мiж знов утвореним Севастопольським антициклоном i попередшм вихором формуеться область циклонного обертання вод.

Ключовi слова: Чорне море, чисельна модель, синоптичш i мезомасштабш вихори.

ABSTRACT Results of a numerical prognostic experiment with 1.64 km resolution on the horizontal coordinates are represented. In contrast to the calculations with more rough resolution, the following dynamical peculiarities are obtained. The Rim Current as a continuous flow directed along the bottom slope is observed up to 400 m. Along the eastern part of the Anatolian coast meso-scale anticyclonic vortices effecting formation of the Batumi anticyclone are regularly formed. Sevastopol, Sinop, Kizilirmak and Caucasian anticyclonic vortices are quasi-periodical ones. Between the newly formed Sevastopol anticyclone and the previous vortex, the area of the cyclonic vortex rotation is formed.

Keywords: Black Sea, numerical model, synoptical and meso-scale vortices.

ISSN 0233-7584. Мор. гидрофиз. журн, 2011, № 1

47