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

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МОРСКИХ СИСТЕМ УДК 561.465

Численное моделирование водообмена через Керченский пролив для различных типов атмосферных воздействий

© 2017 В. В. Фомин1*, Д. И. Лазоренко1, И. Н. Фомина2

'Морской гидрофизический институт РАН, Севастополь, Россия 2Севастопольское отделение Государственного океанографического института им. Н.Н. Зубова, Севастополь, Россия *E-mail: fomin. [email protected]

Поступила в редакцию 02.11.2016 г. После доработки 26.12.2016 г.

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

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

DOI: 10.22449/0233-7584-2017-4-82-93

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

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

Моделированию динамики вод Керченского пролива посвящен ряд публикаций [4 - 9], в которых используется расчетная область с двумя жидкими

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

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

Численная модель. Для расчетов полей течений и уровня моря в Азово-Черноморском бассейне применяется численная гидродинамическая модель Advanced Circulation Model for Shelves Coasts and Estuaries (ADCIRC). Подробное описание модели и результаты ее использования приведены в работах [10 - 13]. Валидация ADCIRC по данным наблюдений за уровнем Азовского моря проводилась в [14, 15]. Исходные уравнения модели имеют вид

ôu ttôu ôu d

-+ U-+ V--fV = —g—

dt dx dy dx

Л + -

gPo

T —Tl

_!_ sx_bx_ _j_

Po H

Ml

H

(1)

dV dV dV d

— + U — + V — + fU = —g — dt dx dy dy

gPo

T —Tu

+ sy by +

Po H

Mv

—, (2) H

Ôl+Ôq±+д^ = 0

dt dx ôy

(3)

здесь х, у, t - горизонтальные координаты и время; U, V - средние по глубине компоненты вектора скорости течений по осям х и у соответственно; 7] - отклонение свободной поверхности моря от невозмущенного состояния; f - параметр Кориолиса; g - ускорение силы тяжести; Pа - атмосферное давление; ро - средняя плотность морской воды; H = h + 7] - динамическая глубина, h -глубина моря; Мх = Ль Д^, Му = Ль Д^ - компоненты горизонтальной турбулентной вязкости, Ah - коэффициент горизонтальной турбулентной вязкости, Д - оператор Лапласа по переменным х, у; qx = иН, qy = УН - компоненты вектора полных потоков.

Компоненты касательных напряжений в (1), (2) определяются следующими выражениями:

' ' (4)

T sx =PfWAWX- + W , Tsy =P<CWAWl + Wy

Tbx = pCJ^U2 + V2 , Tby = PoCdVVu2 + V2 ,

(5)

где ра - плотность воздуха; Wх, Wy - компоненты вектора скорости приводного ветра; Са, Cd - коэффициенты поверхностного и донного трения.

Численный алгоритм ADCIRC основан на методе конечных элементов, использующем треугольные элементы и линейные базисные функции. Для

предотвращения образования вычислительного шума при численном интегрировании системы (1) - (3) последнее уравнение записывается в виде так называемого уравнения Generalized Wave Continuity Equation (GWCE) [11]

dG , n n

--+ т G = 0,

dt 0

8ц 8qx 5qy „ „

где G =--1--— h--— ; т> - произвольный положительный параметр. После

8t 8x 8y

некоторых дифференциальных преобразований уравнение GWCE приводится к следующему эквивалентному виду:

^ + + + ПУ = 0, (6)

8 t 8t 8x ду

где

8U 8U g 8ц2 H 8Pa т -тъ 8ц 8ц

8x 8y 2 8x p0 8— p0 8t 8—

J = a 8V q 8V fa g Ц H 8Pa , Tsy-Тъу + M +т q + V 8Ц gH 8Ц

Jy = -qx^--qy^--fqx-—^----+-+ My + т 0qy + V -gH^~ ■

y 8x y 8y 2 8y po 8y p0 y y 8t 8y

На практике оптимальное значение параметра т выбирается на основе численных экспериментов. Ориентиром для его выбора является условие

1 <то/ттах ^10, в котором ттах = max CdVU2 + V 2/H .

Результаты численных экспериментов. Численное моделирование выполнялось для всего Азово-Черноморского бассейна на неструктурированной расчетной сетке из ~158 тыс. узлов. Размер сторон конечных элементов изменялся в пределах от 60 м в Керченском проливе до 4900 м в глубоководной части Черного моря. Коэффициенты трения в выражениях (4), (5) задавались

формулами С = °,°01(0,15 + 0,067^Ж? + Ж2 ^, С = ё"21 Нуъ, где п = 0,025 -

параметр шероховатости Маннинга. Коэффициент горизонтальной турбулентной вязкости Ah = 10 м2/с, шаг интегрирования по времени составлял 3 с. С учетом вышеизложенных ограничений весовой множитель ть в уравнении (6) принимался равным 0,005. Для упрощения задачи вклад речного стока при моделировании гидродинамических полей не учитывался.

Положение разреза порт Крым - порт Кавказ показано на рис. 1. Длина разреза Ь приблизительно равна 4 км. Расчеты водообмена проводились по вычисленным гидродинамическим полям следующим образом. Вначале с помощью линейной интерполяции определялись компоненты вектора полных потоков и уровень моря на разрезе. Затем интегрированием вдоль прямой Ь вычислялись мгновенный расход воды (0, средняя скорость течений и средний уровень моря (5Ь):

й($)

| Н (I, t )М

, 8Щ) = I ^(1, t)Л,

где дХ + д2у ~Ф) - проекция вектора полных потоков на нормаль

к прямой Ь, = 52° - направление положительной нормали к разрезу, ф -направление вектора полных потоков. Здесь и далее все углы отсчитываются относительно оси х, направленной на восток против часовой стрелки. Таким образом, при Q > 0 суммарный поток воды направлен в Азовское море (прямой поток), при Q < 0 - в Черное море (обратный поток). Для прямого и обратного потоков используются обозначения QA и Qв соответственно.

Рис. 1. Схема северной узости Керченского пролива. Пунктирная линия - разрез порт Крым -порт Кавказ

В первом цикле численных экспериментов анализировалось влияние направления (в ) однородного по пространству ветра на водообмен через Керченский пролив. Задача решалась с нулевыми начальными условиями. Предполагалось, что при ^ > 0 на поверхность бассейна начинает действовать ветер, скорость которого за 12 ч линейно возрастает от 0 до Жо и далее остается неизменной. Для каждого фиксированного в общее время интегрирования составляло 48 ч, что соответствует характерному времени действия устойчивых ветров над Азово-Черноморским бассейном.

На рис. 2 для различных скоростей ветра Жо показаны зависимости расхода воды от направления ветра при ^ = 48 ч. Как видно, представленные кривые имеют хорошо выраженный синусоидальный характер. Максимальным прямым потокам соответствует диапазон углов 60° < в< 105°, максимальным обратным потокам - диапазон углов 240° < в < 285°. Минимальным значени-

ь

ь

ям потоков, а также смене их направления соответствуют диапазоны углов 165° < в < 175° и 345° < в < 355°. Увеличение ЖО вызывает интенсификацию водообмена через пролив. Так, при изменении ЖО в диапазоне 5 - 20 м/с максимальные значения QA и Qв увеличиваются в 5,7 и 6,6 раза соответственно. Связь максимальных значений прямого и обратного потоков со скоростью ветра ЖО в указанном диапазоне можно аппроксимировать следующими линейными уравнениями:

тах QA = 1013 Жо - 2008, тах Qв = -1207Жо + 3324.

20000 15000 10000 5000

С]

S

Oj -5000 -10000 -15000 -20000 -25000

0 45 90 135 180 225 270 315 360 0,°

Рис. 2. Зависимости Q от впри t = 48 ч и различных значениях Wo. Черная кривая соответствует Wo = 5 м/с; синяя - Wo = 10 м/с; красная - Wo = 15 м/с; зеленая - Wo = 20 м/с

На рис. 3 при Wo =10 м/с представлены зависимости Q от времени для характерных значений в. Кривые В1и B2 соответствуют максимальным значениям Q (направление ветра близко к меридиональному). Видно, что в интервале времени 0 - 12 ч происходит монотонное увеличение Q до ±7500 м3/с с последующей стабилизацией. На этих кривых видны незначительные колебания с инерционным периодом, обусловленные вращением Земли. Кривые C1 - C4 соответствуют минимальным значениям Q, когда направление ветра близко к зональному. На этих кривых также прослеживаются инерционные колебания, но их размах существенно больше по сравнению с кривыми B1 и B2.

Таким образом, направление ветра существенно влияет на водообмен в Керченском проливе. Выделяются диапазоны направлений ветра с максимальными и минимальными значениями Qa и Qb . При направлении ветра, близком к зональному, водообмен минимален. В этом случае возрастает роль инерционных колебаний. При всех рассмотренных значениях Wo качественный вид зависимости Q от времени не меняется.

10000

7500

2500

о

2 0

См

-2500

-7500

-10000

В1

С1

¡//02\ ----

\\сз/

С4

В2

1 | | |

12

24

Л Ч

36

48

Рис. 3. Зависимости Q от времени при Жо =10 м/с. Кривая В1 соответствует в = 75°; кривая В2 - в= 240°; кривая С1 - в= 165°; кривая С2 - в= 175°; кривая С3 - в= 355°; кривая С4 - в= = 345°

Анализировались также особенности водообмена через Керченский пролив, возникающие после резкого прекращения действия ветра. Как и в предыдущих расчетах, скорость ветра за первые 12 ч линейно увеличивалась от 0 до Ж0 и далее оставалась постоянной в течение 36 ч. После этого она принималась равной нулю и интегрирование проводилось еще на 72 ч модельного времени.

На рис. 4 показаны зависимости Q от времени после прекращения действия ветра (48 ч < ( < 120 ч). Анализ кривых показывает, что в течение первых суток происходят резкие изменения Q, обусловленные внезапным прекращением действия ветра и сопровождаемые сменой знака Q. Далее возникает режим свободных колебаний с затуханием.

При меридиональном направлении ветра (кривые С1 и С2) период этих колебаний т равен 10 - 12 ч. При зональных ветрах (кривые С3 и С4) период колебаний близок к инерционному. В обоих случаях размах колебаний Q относительно невелик. В средних течениях эти колебания проявляются в виде флуктуаций Ук, равных ~0,02 - 0,03 м/с. Можно предположить, что период т индуцирован продольными сейшевыми колебаниями Азовского моря.

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

Ч211/2

Т -—,

Ра & V 2 ) 2 ,

) =

¿Р + ( /г

где г = ^г'2 + гУ - расстояние от центра циклона (хс, ус) до точки (х, у), в которой определяется Жё; Гх = х - Хс; Гу = у - ус. Атмосферное давление Ра задавалось выражением из работы [4]

Pair) =

P) -SPaCOsZ J при r < R, P0 при r > R,

где Po - фоновое значение атмосферного давления; SPa - перепад давления между центром и периферией циклона; R - радиус циклона. Предполагалось, что направление ветра в приповерхностном слое атмосферы за счет трения отклоняется от касательных к изобарам на угол у = 20° против часовой стрелки, а его скорость меньше скорости градиентного ветра. С учетом этих допущений компоненты вектора скорости ветра в приповерхностном слое атмосферы определялись при a = 90° + у, ¡ = 0,7 следующим образом:

¡Wg ¡W

Wx =--- (rx sin a + ry cos«), Wy =-- (rx cosa- ry sin a) .

r r

Рис. 4. Зависимости Q от времени при Wo =10 м/с после прекращения действия ветра. Кривая C1 соответствует в= 90°; кривая C2 - в = 270°; кривая C3 - в= 0°; кривая C4 - в= 180°

Расчеты проводились при следующих параметрах циклона: R = 300 км, SPa = 10 гПа, x0 = 22° в. д., yo = 46,75° с. ш. Циклон двигался в восточном направлении по зональной траектории (xc = xo + c t, yc = yo) со скоростью c. Общее время интегрирования составляло 10 сут.

Зависимости Q от времени при различных скоростях перемещения циклона показаны на рис. 5, а - в. На рис. 5, г приведены зависимости от времени суммарного расхода воды между бассейнами:

I

ас)=|

0

Рис. 5. Зависимости Q и Qs от времени при различных скоростях движения циклона: а - Q при с = 20 м/с; б - Q при с = 30 м/с; в - Q при с = 40 м/с; г - Qs при с = 20 м/с (черная кривая), с = = 30 м/с (красная кривая), с = 40 м/с (синяя кривая)

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

Заключительный цикл расчетов посвящен оценкам водообмена через Керченский пролив для реальных синоптических ситуаций, которые имели место в Азово-Черноморском бассейне. В качестве атмосферного форсинга использовались поля скорости ветра и атмосферного давления из мезомасштабной метеорологической модели ШЯГ (ecobase.org.ua) с пространственным разрешением 7 км и дискретностью по времени 3 ч. Расчеты проводились для 5 месячных интервалов времени с разными синоптическими условиями.

Рис. 6. Зависимости Жи, Q, Уи, SL от времени для ноября 2007 г.

В качестве примера на рис. 6 показаны результаты расчетов для ноября 2007 г. Как известно, 11 ноября 2007 г. в Азово-Черноморском бассейне произошел мощный шторм, который привел к катастрофическим последствиям в Керченском проливе. На рис. 6, а показана проекция скорости ветра на

нормаль к разрезу WN = ^+ W^ cosO^ -в) . Для вычисления WN использовались средняя скорость и среднее направление ветра на разрезе. Сопоставление кривых показывает, что изменения расхода воды (рис. 6, б) и средней скорости течений (рис. 6, в) происходят синфазно с изменениями проекции скорости ветра на нормаль к разрезу (рис. 6, а). При этом средний уровень моря изменяется в противофазе с указанными выше характеристиками. Такой ход зависимостей является общим для всех 5 интервалов времени.

В таблице приведены экстремальные значения величин Q, Vn и SL в северной узости Керченского пролива для рассматриваемых месячных интервалов времени. Видно, что наибольшие положительные значения Q (черноморский поток) достигают 12000 - 16000 м3/е, наибольшие отрицательные значения Q (азовский поток) - 11000 м3/е. Средние скорости течений из Черного моря в Азовское составляют 0,9 м/c. При этом течения из Азовского моря в Черное имеют несколько меньшие скорости (~0,6 м/c). Наибольшие спады и подъемы уровня моря находятся в пределах -0,4 ... 0,4 м. Приведенные в таблице данные согласуются с результатами расчетов водообмена с использованием натурных измерений из работы [3], где показано, что максимальные расходы воды в период штормов составляют 11000 - 15000 м3^ при средней скорости течений ~45 - 47 см/с.

Экстремумы расхода воды Q, нормальной составляющей средней скорости течений Ук и среднего уровня моря БЬ в северной узости Керченского пролива для месячных интервалов времени

Расчетные месяцы maxQ, minQ, maxVN, м/c minVN, м/c max SL, м min SL м

Ноябрь 2007 г. 16220 -10909 0,93 -0,59 0,17 -0,34

Март 2013 г. 12880 -9700 0,79 -0,52 0,23 -0,38

Июль 2013 г. 6604 -8043 0,38 -0,43 0,14 -0,13

Сентябрь 2013 г. 9543 -7586 0,56 -0,41 0,17 -0,17

Сентябрь 2014 г. 14990 -1392 0,92 -0,61 0,41 -0,39

Заключение. Выполнены расчеты водообмена через северную узость Керченского пролива для различных типов атмосферных воздействий. В качестве входных данных использовались результаты численного моделирования течений и уровня моря для всего Азово-Черноморского бассейна на неструктурированной расчетной сетке со сгущением в Керченском проливе. Анализ результатов расчетов позволил сделать следующие выводы.

Направление ветра в существенно влияет на величину расходов воды Q. Диапазоны 60° < в < 105° и 240° < в < 285° соответствуют максимальным

значениям |Q|, диапазоны 165° < 9 < 173° и 345° < 9 < 353° - минимальным значениям |Q|. После прекращения действия ветра в Керченском проливе возникают сейшеобразные осцилляции Q с периодами 10 - 12 ч.

При распространении над Керченским проливом атмосферного циклона интенсивность водообмена определяется скоростью его перемещения c. При уменьшении c увеличивается время воздействия циклона на акваторию Керченского пролива, что приводит к усилению суммарного водообмена в проливе.

Величины |Q| для реальных синоптических ситуаций могут достигать 11000 - 16000 м3/с, что согласуется с результатами расчетов Q, полученными по натурным данным.

Расчеты проводились на вычислительном кластере Морского гидрофизического института.

Работа выполнена в рамках государственного задания по теме № 08272014-0010 «Комплексные междисциплинарные исследования океанологических процессов, определяющих функционирование и эволюцию экосистем Черного и Азовского морей, на основе современных методов контроля состояния морской среды и грид-технологий».

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

1. Альтман Э.Н. К вопросу об изменчивости расходов воды в Керченском проливе (по натурным наблюдениям) // Тр. ГОИН. - 1976. - Вып. 132. - С. 17 - 28.

2. Альтман Э.Н. Динамика вод Керченского пролива // Гидрометеорология и гидрохимия морей СССР. Том 4. Черное море. Вып. 1. Гидрометеорологические условия. - Л.: Гид-рометеоиздат, 1991. - С. 291 - 328. - URL: http://www.geokniga.org/bookfiles/geokniga-t4-chernoe-more-vyp-1 -gidrometeorologicheskie-usloviya-1991.pdf (дата обращения: 10.08.2016).

3. Дьяков Н.Н., Фомина И.Н., Тимошенко Т.Ю., Полозок А.А. Особенности водообмена через Керченский пролив по данным натурных наблюдений // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. -Севастополь: МГИ, 2016. - Вып. 1. - С. 62 - 67.

4. Гидрометеорологические условия морей Украины. Том 1. Азовское море. - Севастополь: ЭКОСИ-Гидрофизика, 2009. - 402 с.

5. Иванов В.А., Шапиро Н.Б. Моделирование течений в Керченском проливе // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофизика, 2004. - Вып. 10. - С. 207 - 232.

6. Ivanov V.A., Cherkesov L.V., Shul'ga T.Y. Extreme deviations of the sea level and the velocities of currents induced by constant winds in the Azov Sea // Physical Oceanography. - 2011. - 21, Iss. 2. - P. 98 - 105. - doi:10.1007/s11110-011-9107-5

7. Матишов Г.Г., Чикин А.Л. Исследование ветровых течений в Керченском проливе с помощью математического моделирования // Вестник Южного научного центра РАН. - 2012. - 8, № 2. -С. 27 - 32. - URL: http://ssc-ras.ru/files/files/27-32_Matishov1_.pdf (дата обращения 10.08.2016).

8. Тучковенко Ю.С. Численная математическая модель циркуляции вод в Керченском проливе // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофизика, 2002. - Вып. 6. - С. 223 - 232.

9. Фомин В.В., ПолозокА.А., Фомина И.Н. Моделирование циркуляции вод Азовского моря с учетом речного стока // Морской гидрофизический журнал. - 2015. - № 1. -С. 16 - 28. - doi:10.22449/0233-7584-2015-1-16-28

10. Luettich R.A., WesterinkJ.J. Formulation and Numerical Implementation of the 2D/3D ADCIRC Finite Element Model Version 44.XX. - 2004. - URL: http://www.unc.edu/ims/adcirc/publications/2004/2004_Luettich.pdf (дата обращения: 26.06.2017).

11. Luettich Jr. R.A., Westerink J.J., Scheffner N.W. ADCIRC: an advanced three-dimensional circulation model for shelves coasts, and estuaries, report 1: theory and methodology of ADCIRC-2DDI and ADCIRC-3DL. Dredging research program technical report DRP-92-6. -Vicksburg: U.S. Army Engineers Waterways Experiment Station, 1992. - 137 p. - URL: http://www.dtic.mil/get-tr-doc/pdf?AD=ADA261608 (дата обращения: 10.08.2016).

12. Dietrich J.C., Zijlema M., Westerink J.J. et al. Modeling hurricane waves and storm surge using integrally-coupled, scalable computations // Coast. Engineer. - 2011. - 58, Iss. 1. -P. 45 - 65. - doi:10.1016/j.coastaleng.2010.08.001

13. Sebastian A.G., Proft J.M., Dietrich J.C. et al. Characterizing hurricane storm surge behavior in Galveston Bay using the SWAN + ADCIRC model // Ibid. - 2014. - № 88. - P. 171 -181. - doi:10.1016/j.coastaleng.2014.03.002

14. Фомин В.В., ПолозокА.А. Технология моделирования штормовых нагонов и ветрового волнения в Азовском море на неструктурированных сетках // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. -Севастополь: ЭКОСИ-Гидрофизика, 2013. - Вып. 27. - С. 139 - 145.

15. Fomin V.V., PolozokA.A., KamyshnikovR.V. Wave and storm surge modelling for Sea of Azov with use of SWAN +ADCIRC // Geoinformation Sciences and Environmental Development: New Approaches, Methods, Technologies. Collection of articles. - Rostov-on-Don: Publishing house SSC RAS, 2014. - P. 111 - 116.

Numerical modeling of water exchange through the Kerch Strait for various types of the atmospheric impact

V. V. Fomin1*, D. I. Lazorenko1, I. N. Fomina2

Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia 2Sevastopol Branch of the N.N. Zubov State Oceanographic Institute, Sevastopol, Russia *e-mail: fomin. dntmm@gmail. com

Features of water exchange through the Kerch Strait are examined for various types of the atmospheric impact. The results of numerical modeling of currents and sea level in the Azov - Black Sea basin on the non-uniform computational grid concentrated in the strait are used as the input data. The ranges of wind directions at which water exchange through the strait is the most intense are identified. It is revealed that the wind field generated by an atmospheric cyclone intensifies total water exchange through the strait with decrease of the cyclone movement speed. The results of calculating water flow in the strait for real synoptic situations are consistent with the estimates obtained by the other authors from the field data.

Keywords: Sea of Azov, Black Sea, Kerch Strait, wind currents, sea level, water exchange, numerical simulation.