CC BY
33
УДК 551.465
Расчет и анализ энергетики циркуляции вод в прибрежных районах Черного моря
© 2017 С.Г. Демышев, О.А. Дымова*
Морской гидрофизический институт РАН, Севастополь, Россия * E-mail: olgadym@yahoo.com
Поступила в редакцию 25.07.2016 г. После доработки 13.10.2016 г.
Представлены результаты анализа компонентов бюджета кинетической и потенциальной энергии Черного моря, рассчитанных с помощью вихреразрешающей модели Морского гидрофизического института. Численные эксперименты выполнены с пространственным разрешением 1,6 км при учете реального атмосферного воздействия за 2006 г. Рассмотрена пространственно-временная изменчивость слагаемых в уравнениях бюджета кинетической и потенциальной энергии на среднегодовом и сезонном масштабах для отдельных районов моря.
Показано, что энергетически наиболее значимыми компонентами бюджета были работа силы плавучести в северной и северо-восточной частях моря, работа силы ветра в южной его части, работа силы давления и сил трения над свалом глубин в зоне Основного Черноморского течения. Формирование мезомасштабных вихрей у Черноморских побережий Крыма и Северного Кавказа происходило преимущественно в весенне-летний сезон в результате бароклинной неустойчивости. Основной фактор, определяющий мезомасштабную динамику у южного и юго-восточного побережий во все сезоны при слабом ветре, - влияние особенностей береговой черты на периферию Основного Черноморского течения.
Ключевые слова: Черное море, численное моделирование, кинетическая энергия, потенциальная энергия, мезомасштабные вихри.
DOI: 10.22449/0233-7584-2017-3-49-62
Введение. С развитием современных методов контактного и дистанционного зондирования вод Мирового океана стали доступны данные о вихревых структурах с пространственными размерами порядка нескольких километров и временной изменчивостью с масштабом от нескольких суток до нескольких недель. Относительно локального бароклинного радиуса деформации Россби Rd эти структуры разделяют на мезомасштабные квазигеострофи-ческие вихри с радиусом больше Rd, для которых число Россби Ro гораздо меньше единицы, и субмезомасштабные агеострофические вихри с радиусом меньше Rd и значением числа Ro ~ 1 [1]. В прибрежной зоне Черного моря, где Rd в среднем составляет 7,5 км, наблюдаются оба типа вихрей [2, 3]. Необходимость исследования циркуляции на масштабе менее 10 км диктуется несколькими причинами. Во-первых, вихревые структуры мезо- и субме-зомасштабного характера могут осуществлять обмен массой, теплом и солью между прибрежной и глубоководной частями и тем самым играть существенную роль в гидродинамике Черного моря. Во-вторых, итоги работ по проекту MyOcean, в котором проведено сравнение результатов численных расчетов циркуляции по трем моделям с разрешением от 5 до 8 км по горизонтали [4], продемонстрировали, что горизонтальное разрешение может существенно влиять на точность прогноза гидрофизических характеристик. Кроме того,
для решения задач мониторинга прибрежных экосистем востребована подробная информация о состоянии морской среды. Изучению и предсказанию последствий природных катастрофических явлений Азово-Черноморского региона на основе численных моделей было посвящено множество работ С.Ф. Доценко [5].
Существует несколько вихреразрешающих моделей динамики Черного моря. В Институте вычислительной математики РАН и Институте океанологии им. П.П. Ширшова РАН разработаны модели с высоким разрешением [6, 7], с помощью которых исследуется изменчивость в Черном и Азовском морях. Ряд диагностических и прогностических экспериментов по расчету гидрофизических полей Черного моря выполнен с использованием моделей РОМ [8], NEMO [9] и МГИ [10] с учетом реального воздействия атмосферы и спутниковой информации.
Эффективным инструментом для оценки роли различных физических факторов в динамике течений является анализ энергетического баланса в моделируемой системе. В статье [11] численные оценки баланса кинетической энергии используются при выборе критериев окончания адаптационных расчетов гидродинамических характеристик. В работе [12] выполнен численный анализ энергетических переходов при изучении механизмов генерации аномалий гидрофизических параметров в Северном Ледовитом океане. В книге [13] рассчитан энергетический цикл Мирового океана с акцентом на вклад мезомасштабных вихрей в энергетику глобальной циркуляции. Энергетические балансы полузамкнутых морей (в том числе и Черного моря) приведены в работе [14], в которой показано, что работа силы плавучести является одним из основных факторов формирования мезомасштабных особенностей динамики в полузамкнутых морях.
Энергетика климатической циркуляции Черного моря, рассчитанной на сетке 5 х 5 км, представлена в статье [15], в которой определены пять энергоактивных зон, расположенных вдоль свала глубин в районе северо-западного шельфа (СЗШ), у западной части Анатолийского побережья, в юго-западном углу моря, в центральной части Кавказского побережья и около восточного Крыма. В работе [16] на основе конечно-разностных уравнений изменения кинетической (КЭ) и потенциальной энергии (ПЭ) рассчитаны составляющие бюджета энергии в 2006 г. с горизонтальным разрешением 1,6 км при учете реального атмосферного воздействия. Проведен анализ осредненных за год и за каждый сезон интегральных компонентов бюджета КЭ и ПЭ. Показано, что в течение 2006 г. изменение энергии Черного моря определялось преимущественно балансом между вкладами работ силы плавучести, сил трения и силы ветра.
Статья посвящена исследованию энергетики Черного моря и является продолжением анализа результатов численного эксперимента [16]. Рассмотрена пространственно-временная изменчивость проинтегрированных по вертикали энергетических характеристик. Целью данной работы было проанализировать наиболее энергетически значимые компоненты бюджета кинетической и потенциальной энергии и определить возможные механизмы формирования и эволюции мезомасштабных вихрей.
н
Численный эксперимент. Исходя из предположения, что <ф>н = ,
о
энергетические уравнения в символьном виде можно записать следующим образом [15]:
<Е{ >н = СА "^(С2/2) - Advhor(Р)-Advhw(E) + <П^Е >н +
+ т ^ Е - От-щ. (Е) - < Ыля^ (Е) >н-< алл^г(Е) > , (1)
где Et - скорость изменения КЭ; £ - возвышение свободной поверхности; £0 = Е(х, у, 0, 0; Advhor (Р) =< (иР)^ + (уР)^ >н - работа силы давления;
Advhor (Е) =< (иЕ)х + (уЕ)у >н - адвекция КЭ; < П ^Е >н =< ^р >н - работа силы плавучести; т ^ Е = щтх + у,ту - работа касательных напряжений
трения ветра; О/ллм (Е) = уг
2 | 2=н ( X, у) 2 | 2=н ( X, у )
- изменение КЭ
V /
за счет трения о дно; < Огля^ (Е) >н(ииг + уу - изменение КЭ за счет
тт О О
вертикального внутреннего трения; < Ом^,,,. (Е) > =уя [(V и) + (V у) ] -изменение КЭ за счет горизонтального трения, куда входит внутреннее горизонтальное трение и трение о боковые стенки.
<П >н = -^уга(П-<П^ Е >н +< 0#1юг(П >н --<(П)>н +<а//^иг Ъо1(и)>н -<т^П>н, (2)
где П - скорость изменения ПЭ; Advhor(П) =< (иП)х + (уП)у >н - адвекция
ПЭ; < О1//Ъог: (П) >н =< -к н (V4П - gzQH) >н - изменение ПЭ за счет горизон-
>Ъ)Г ^ Ч--- Vй ^ Чх^Уу1-Чу
/^Г(п) > —н (» П-,
тальной турбулентной диффузии; < (П) >н =< §р(кг)2 >н - изменение
ПЭ за счет неоднородности коэффициента вертикальной турбулентной диффузии по глубине; <О/й ™_Ь(* (П >н = §(кVр)г=о -§(кVр)г=н(х,у) - изменение ПЭ за счет разницы между придонной и поверхностной плотностями; < ЩГ£П) >н = § < 2кrQr) >н - слагаемое, возникающее при интегрировании уравнения состояния вследствие его нелинейного характера. Остальные обозначения общепринятые.
В отличие от работы [15], в настоящем исследовании расчеты выполнены на сетке 1,6 км, вертикальные турбулентные процессы параметризованы на основе теории Меллора - Ямады уровня замыкания 2.5 [17], учтены ветер, потоки тепла, осадков и испарения, полученные по данным атмосферной модели ЛЬАБШ [18]. Для задания гидродинамических характеристик в устьях рек и проливах использованы климатические данные [19]. Для удобства восприятия знак интегрирования по вертикали < > н далее в тексте опущен.
Анализ результатов. На рис. 1 представлено пространственное распределение некоторых осредненных за год слагаемых уравнений (1), (2). Адвекция КЭ на среднегодовом масштабе характеризовалась близкими к нулю значениями, малыми по сравнению со значениями других слагаемых. Из рис. 1, а видно, что максимальные и минимальные значения работы силы давления наблюдаются между изобатами 200 - 900 м в западной части моря и у берегов Турции. В этих же районах над свалом глубин выявлены наибольшие среднегодовые значения скорости в стрежне Основного Черноморского течения (ОЧТ), которые приводят к росту величины адвективной производной и, соответственно, интенсификации работы силы давления. В зоне Батумского антициклона расположены взаимно компенсационные области величины Advhor поэтому суммарный вклад работы силы давления в изменение КЭ в этом районе в среднем за год незначителен. Среднегодовое распределение работы силы плавучести представлено на рис. 1, б. Работа силы плавучести зависит от вертикальной скорости и плотности: П Е = gwp. Сопоставление с картами вертикальной скорости показало, что в глубоководной части (за исключением юго-восточного угла моря) области положительных и отрицательных значений величины П Е соответствуют зонам опускания и подъема уровня воды, которые определяются знаком вертикальной скорости. На СЗШ преобладают положительные значения П Е, при этом корреляция с пространственным распределением вертикальной скорости не выявлена. Следовательно, работа силы плавучести на СЗШ в большей степени была обусловлена полем плотности, которое определялось потоками тепла, осадков и испарения на поверхности моря и втоками пресных речных вод. В среднем за год в этом районе наблюдался переход доступной потенциальной энергии (ДПЭ) в КЭ и тем самым происходило выравнивание изопикниче-ских поверхностей. Таким образом, за исследуемый период времени в районе СЗШ процессы охлаждения и зимней конвекции были энергетически более значимы по сравнению с весенне-летним прогревом вод и формирование ме-зомасштабных вихрей, по мнению авторов, происходило на фоне баротропи-зации течений. Однако данный вывод является лишь предположением, так как оценить баротропную и бароклинную составляющие кинетической энергии (как, например, в работе [20]) на основе используемого нами подхода затруднительно.
В среднем за год КЭ возрастала за счет притока от ветра: наибольший вклад отмечен в юго-западном углу моря и у восточной части Анатолийского побережья (рис. 1, в). Такое распределение обусловлено полем напряжений трения ветра. Анализ поля ветра показал, что в 2006 г. наиболее сильные и часто повторяющиеся ветра имели северо-восточное направление с максимальными значениями скорости у западного побережья. Отрицательные значения работы силы ветра около западной части Анатолийского побережья (рис. 1, в) означают, что ветер и течения на поверхности направлены в противоположные стороны: ОЧТ направлено на северо-восток, а вектор напряжения ветра - преимущественно на юго-запад. Вследствие этого развиваются процессы сдвиговой неустойчивости, которые являются одним из механизмов образования мезомасштабных вихрей [21]. Анализ карт мгновенных течений и поля ветра показал, что в указанном районе при слабых ветрах в результате набегания течения на береговые выступы вихревые структуры фор-52 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2017
мировались между ОЧТ и берегом. При сильных ветрах вихри не наблюдались. Таким образом, около западной части Анатолийского побережья в среднем за год влияние ветра на структуру циркуляции оказалось энергетически более слабым по сравнению с влиянием орографических особенностей береговой черты.
Рис. 1. Осредненные за год слагаемые (эрг/см2) в уравнениях (1), (2): а - АЛьг (Р); б - П ^ Е;
в - т ^ Е; г - Dissver (Е); д - АЛы (П); е -
Значительные потери КЭ вследствие вертикального трения выявлены в районе юго-восточного побережья Крыма и западной части Черного моря с максимальными значениями у берегов Румынии и Болгарии (рис. 1, г), что объясняется интенсификацией процессов перемешивания и вертикального обмена в зонах, где были отмечены наиболее сильные и часто повторяющиеся ветры. При расчете Dissver (Е) касательные напряжения трения ветра непосредственно учитываются в качестве краевых условий не только в уравнениях движения, но и при вычислении коэффициентов вертикального обмена с использованием параметризации Меллора - Ямады [17]. Поэтому усиление ветра приводит к увеличению вертикального градиента скорости, коэффициентов вертикального турбулентного обмена и, как следствие, к большим по модулю значениям Dissver (Е). Сток КЭ за счет горизонтальной диссипации наиболее интенсивен в струе ОЧТ. Так как коэффициенты горизонтальной турбулентной вязкости постоянны, то потеря КЭ в результате трения зависит исключительно от величины горизонтальных градиентов скорости, которые
в ОЧТ имеют большие значения, чем в остальных областях моря. При этом основной вклад в величину Disshor (Е) вносили верхние модельные горизонты в зимний период, когда ОЧТ наиболее интенсивно. Летом вследствие ослабления ОЧТ и ветров развиваются многочисленные вихри различных масштабов и уменьшаются потери энергии вследствие горизонтального трения по сравнению с холодным сезоном.
В среднем за 2006 г. самые значимые вклады в изменение ПЭ вносили потоки тепла и соли на поверхности моря и адвекция ПЭ. Наибольшие по модулю значения адвективного слагаемого в уравнении (2) локализованы над свалом глубин (рис. 1, д), где увеличение скорости в ядре ОЧТ приводит к интенсификации адвективного переноса. В глубоководной части моря к западу от 35° в. д. наблюдалось чередование зон положительных и отрицательных значений величины Аа^ьог (П). Аналогичная структура с обратным знаком прослеживалась в распределении слагаемого ЛёУь_ог (Р) на рис. 1, а. Эти компоненты баланса энергии являются взаимно-компенсационными, их поведение определяется вертикальной скоростью. В работе [22] указывается на наличие в Черном море стоячих волн Россби, структура которых по пространственно-временным масштабам соответствует распределению характеристик на рис. 1, а и 1, д.
В пространственном распределении осредненного за год члена О/ЙС(П) (рис. 1, е) выявлено соответствие среднегодовому пространственному распределению полного потока тепла из атмосферы: там, где поток тепла направлен из атмосферы в море (вода нагревается), слагаемое Д/Й™ (П) имеет положительное значение. Таким образом, на среднегодовом масштабе поток тепла наиболее значительно влияет на изменение ПЭ в верхнем слое. Максимальное значения Щй*™ (П) в районах устьев рек Дунай и Риони (рис. 1, е) обусловлены весенней интенсификацией речного стока. При увеличении разности между поверхностной и придонной плотностями (то есть при увеличении значения вертикального градиента плотности) растет отклонение изопикнических поверхностей от среднего «невозмущенного» уровня и, как следствие, повышается запас ДПЭ. Из литературы известно, что ДПЭ является источником энергии мезомасштабных вихрей [23]. Слагаемое
Й(П) мало по сравнению с другими компонентами баланса ПЭ, поэтому в дальнейшем не рассматривается.
Таким образом, из рис. 1 видно, что наиболее интенсивный энергообмен в 2006 г. происходил в районах ОЧТ, у западного побережья, в юго-восточной части моря, около Крымского и Анатолийского побережий. Ниже представлены результаты анализа сезонной изменчивости компонентов бюджета КЭ и ПЭ. Значения проинтегрированных по вертикали энергетических слагаемых в уравнениях (1) и (2) были осреднены посезонно. Для каждого сезона оценены наибольшие вклады в баланс энергии. Мы рассматривали гидрологические сезоны, то есть трехмесячные периоды, начинающиеся 1 января, 1 апреля, 1 июля и 1 октября.
12 12.5 13 13.5 14 14.5 15 15.5 16 16.5 17 17.5 18
Рис. 2. Поле солености (%о) на разрезе 45,2° с. ш.: а - 22.05.06 г., б - 23.12.06 г.
Для СЗШ и западного побережья в осенне-зимний сезон наиболее значимые компоненты баланса КЭ - вклад от ветра и потери вследствие вертикального трения. Эти величины характеризуют перестройку поля течений, вызванную штормовыми ветрами (скорости ветра порядка 10 м/с и более) над западной частью моря в холодный период. В это время в поле скорости
наблюдались струйные течения, направленные преимущественно на юго-запад. В западной части моря вклад величины Аа^ьог (Р) в КЭ имел наибольшие значения над свалом глубин во все сезоны года. Изменчивость ПЭ в районе СЗШ наиболее интенсивна весной и летом в результате усиления речного стока. На рис. 2 представлены карты поля солености на разрезе, соответствующем широте устья р. Дунай для холодного и теплого сезонов. Видно, что весной вток речных вод приводит к изгибу изохалин, в то время как зимой они практически вертикальны. Таким образом, ПЭ растет вследствие увеличения разницы между поверхностной и придонной плотностями. При этом возрастает отклонение изопикнических поверхностей от невозмущенного уровня, что приводит к увеличению запаса ДПЭ.
Сезонная изменчивость слагаемых т ^ Е и Dissver (Е) у крымского побережья определяется тем же фактором, что и на западном побережье, а именно: заданными полями напряжения трения ветра. Это значит, что в осенне-зимний сезон вклад от ветра и потеря вследствие трения - наиболее значимые компоненты баланса энергии, а течения носят преимущественно струйный характер. Во все сезоны года у юго-восточного побережья Крыма наблюдаются взаимно компенсационные зоны Л^ъ.ог (Е) и ЛЛьог (Р), в итоге доминирующей величиной в балансе КЭ становится работа силы плавучести (рис. 3), значение которой определяется знаком вертикальной скорости. Это значение большую часть времени положительно, то есть потенциальная энергия переходит в кинетическую. При отсутствии сильных ветров этот процесс завершается формированием короткоживущих мезомасштабных вихрей в поле скорости.
Рис. 3. Осредненная за сезон работа силы плавучести (эрг/см2): а - зимой, б - летом
В балансе ПЭ в холодный период года преобладает слагаемое (П),
величина которого максимальна в Каркинитском заливе. Механизм изменения
Щ(П) следующий. С увеличением глубины значение коэффициента вертикальной турбулентной диффузии к к уменьшается. А так как оно рассчитывается с использованием параметризации Меллора - Ямады [17], в которой ветер учитывается непосредственно в качестве граничного условия на поверхности, то при сильных и длительных ветрах в холодные сезоны года вертикальный 56 МОРСКОЙ ГИДРОФИЗИЧЕСКИЙ ЖУРНАЛ № 3 2017
диффузионный поток возрастает за счет увеличения вертикального градиента (кг )2. Весной прирост ПЭ определяется большими положительными значениями полного вертикального диффузионного потока около побережий Крыма и Кавказа, усиление которого соответствует заданному потоку тепла на поверхности моря.
Анатолийское побережье в течение года характеризуется пространственно-временной неоднородностью компонентов бюджета КЭ и относительно однородным распределением составляющих бюджета ПЭ (за исключением ЛОуьог (П)). Для удобства изложения побережье Турции к западу от 34° в. д. будем называть западной частью, к востоку от 34° в. д. - восточной. Западная часть отличается крутым континентальным склоном, восточная -относительно широким свалом глубин и изрезанностью береговой черты. Зимой и осенью работа силы давления уменьшает КЭ в восточной части и увеличивает в западной. В холодный период года, когда вода охлаждается, возрастает плотность воды и, соответственно, давление. Увеличение давления и интенсификация течений над узким континентальным склоном в западной части приводят к значительному росту величины Лакъог (Р). Во все сезоны у берегов Турции фиксируются области наиболее интенсивной диссипации КЭ. К усилению диссипативных процессов приводит интенсификация ОЧТ над узким континентальным склоном. При анализе распределения слагаемого т ^ Е, осредненного посезонно, выявлена следующая особенность: в западной части Анатолийского побережья работа силы ветра отрицательна, в восточной - положительна. Сопоставление карт течений на верхнем расчетном горизонте с пространственным распределением слагаемого, характеризующего мгновенную работу силы ветра, показало, что при слабых ветрах областям отрицательных значений т ^ Е соответствуют зоны локализации мезомас-штабных вихрей. На рис. 4 в качестве примера представлены поля течений и т ^ Е на 19.03.06 г. Отметим, что значение модуля вектора касательных напряжений трения ветра над изучаемым районом на указанную дату изменяется в пределах 0,2 - 0,5 дин/см2, что соответствует скоростям приводного ветра в диапазоне 1,5 - 2 м/с. Скорость ветра высчитывалась по формуле
т = рС^Г^ [24]. Из рис. 4 видно, что слагаемое т ^ Е увеличивает кинетическую энергию ОЧТ (красные зоны), а около берега формируются области стока КЭ (синие зоны). Расположение областей стока КЭ (рис. 4, б) вдоль берега четко совпадает с локализацией мезомасштабных вихрей в поле скорости (рис. 4, а). Слагаемое т ^ Е стоит в правой части выражения (1) и определяется произведением компонент векторов скорости течений и касательного напряжения трения ветра (щтх + у,ту ) . Если исследуемая величина отрицательна, то происходит торможение потока и при набегании на неровности береговой линии в шельфовых зонах формируются запрепятственные антициклонические вихри [25]. Глубина шельфа в восточной части Анатолийского побережья в среднем составляет 100 м. Следовательно, при слабом ветре орография берега оказывает существенное влияние на течения и за мысами в относительно мелководных районах возникают мезомасштабные вихри. При усилении ветра ОЧТ прижимается к берегу, площадь зон стока КЭ
уменьшается и мезомасштабные вихри в прибрежной зоне не развиваются. Таким образом, ослабление ветрового воздействия в сочетании с влиянием орографии береговой черты могут являться основными предпосылками для генерации мезомасштабных особенностей циркуляции в южной части Черного моря. Исходя из вышесказанного, мы предполагаем, что в областях, где слагаемое т ^ Е отрицательно, КЭ течений трансформируется в вихревую энергию. Однако данная гипотеза требует дальнейшего исследования.
Рис. 4. Поле скорости (см/с) (а) и поле т ^ Е (эрг/с2) (б) на 19.03.06
Несмотря на отсутствие значительных среднегодовых вкладов в изменение энергии в районах Батумского и Севастопольского антициклонов,
адвективный перенос существен в этих областях на сезонном масштабе. Так, с весны формируются взаимно компенсационные области Лёу^ог (Е) иЛёуьог(Р), соответствующие расположению антициклонов (рис. 5, а, б). В зонах соприкосновения антициклонов с основным циклоническим круговоротом наблюдаются наибольшие скорости течений, соответственно, возрастают и горизонтальные градиенты скорости. В результате этого вклад адвективных слагаемых увеличивается. В районе Батумского антициклона ЛЛьог (П) вносит главный вклад в изменение ПЭ (рис. 5, в), что обусловлено большими градиентами скорости и большой величиной ПЭ в глубоководной части моря. Летом ОЧТ ослабевает, а Севастопольский и Батумский антициклоны
интенсифицируются, что в энергетическом аспекте приводит к увеличению площади областей, в которых адвективный перенос играет главную роль в энергообмене. В теплый период года по периферии Батумского антициклона наблюдается также рост величины П Е (рис. 5, г). Многочисленные мелкие особенности изолиний в восточной части моря проявляются в результате осреднения по времени. Сопоставление мгновенных полей скорости и работы силы плавучести показало соответствие зон интенсивной изменчивости П Е, обусловленной изменением плотности, зонам локализации мезомасштабных вихрей, что свидетельствует об усилении процессов бароклинной неустойчивости. На юго-востоке в весенне-летний период за счет поступления речных вод возрастает вклад слагаемого ШГ(П)
в изменение ПЭ, при этом растет запас ДПЭ вследствие увеличения отклонения изопикнических поверхностей от «невозмущенного» уровня. Для осеннего периода характерна интенсификация ОЧТ (увеличение скорости и уменьшение ширины потока), Севастопольский и Батумский антициклоны ослабевают, поэтому площадь областей, соответствующих большим значениям адвективных слагаемых в уравнении (1), уменьшается. Охлаждение вод и увеличение толщины верхнего квазиоднородного слоя приводят к росту ПЭ, однако запас ДПЭ уменьшается.
28° 30° 32° 34° 36° 38° 40° 38. 39= 40° 41° в.Д.
Рис. 5. Осредненные за весну адвективные потоки (эрг/см2): а - АЛиог (Е), б - АЛъог (Р), в -АЛьог (П); г - П ^ Е
Зоной, где происходит интенсивная трансформация энергии, является ОЧТ, стрежень которого локализован над континентальным склоном. Прежде всего, ОЧТ является основной областью стока энергии, обусловленного горизонтальными диссипативными и диффузионными процессами. В осенне-зимний период, когда ОЧТ наиболее интенсивно, потери КЭ вследствие горизонтального трения соизмеримы с величиной адвективных вкладов. Максимальный сток КЭ наблюдался западнее оконечности Крымского п-ова и около западной части Анатолийского побережья. Анализ изменения значения
Disshoi (E) с глубиной показал, что основная часть энергии терялась зимой в верхнем 75-метровом слое. Весной и летом работа силы трения уменьшилась на порядок. Коэффициент горизонтальной турбулентной вязкости постоянный, поэтому величина потерь КЭ в результате трения определяется исключительно горизонтальным градиентом скорости течений. Уменьшение ПЭ, обусловленное горизонтальной диффузией, также наиболее велико в зоне ОЧТ, однако величина Diffhoi(П) меняется слабо на межсезонном масштабе.
Выводы. В результате численного эксперимента рассчитаны компоненты бюджета кинетической и потенциальной энергии Черного моря в 2006 г. Рассмотрена пространственно-временная изменчивость основных сил, определяющих особенности циркуляции на среднегодовом и сезонном масштабах для различных районов моря. Анализ карт осредненных по вертикали компонентов бюджета энергии показал, что в среднем за 2006 г. наиболее интенсивный энергообмен происходил у западного побережья и в юго-восточной части Черного моря, около Крымского и Анатолийского побережий, в районе ОЧТ.
В холодный период у западного побережья и на СЗШ штормовые северовосточные ветры обусловливали преимущественно струйный режим течений, а изменение КЭ года определялось балансом между вкладом работы силы ветра и потерей вследствие вертикального трения. Весной и летом в результате усиления речного стока на СЗШ увеличивались отклонения изопикниче-ских поверхностей от «невозмущенного» уровня, что привело к увеличению запаса ДПЭ. Анализ работы силы плавучести в районе СЗШ показал, что при отсутствии сильных ветров ДПЭ переходила из потенциальной в кинетическую энергию, и в результате формировались мезомасштабные вихри.
Около Крымского побережья в течение года доминирующим компонентом баланса энергии была работа силы плавучести, при этом ее величина определялась знаком вертикальной скорости. Показано, что работа силы плавучести у берегов Крыма была преимущественно положительна, то есть потенциальная энергия переходила в кинетическую.
Основной вклад в изменение КЭ осенью и зимой возле Анатолийского побережья вносила работа силы давления, большая величина которой обусловлена интенсификацией адвективного переноса вследствие увеличения скоростей течений в стрежне ОЧТ. Анализ пространственного распределения работы силы ветра показал, что во все сезоны при ослаблении ветрового воздействия под влиянием орографии береговой черты возникали зоны стока КЭ, которая, по-видимому, трансформировалась в энергию мезомасштабных вихрей. При сопоставлении карт работы сила ветра с модельным полем скорости в южной части Черного моря в 2006 г. выявлено пространственное совпадение областей локализации мезомасштабных вихрей с зонами стока КЭ.
Наибольшие вклады в баланс энергии в районах Севастопольского и Батумского антициклонов вносили адвекции КЭ и ПЭ. Пространственно -временная изменчивость этих энергетических потоков обусловлена величиной горизонтальных градиентов скорости, которые возрастали при увеличении скорости течений в зонах соприкосновения антициклонов с ОЧТ. Область ОЧТ являлась зоной стока КЭ и ПЭ вследствие интенсивных горизонтальных диссипативных и диффузионных процессов, при этом зимой в верхнем 75-метровом слое терялась основная часть КЭ.
Работа выполнена по теме госзадания 0827-2014-0010 «Комплексные междисциплинарные исследования океанологических процессов, определяющих функционирование и эволюцию экосистем Черного и Азовского морей, на основе современных методов контроля состояния морской среды и гридтехнологий» (шифр «Фундаментальная океанология»).
СПИСОК ЛИТЕРАТУРЫ
1. Thomas L.N., Tandon A., Mahadevan A. Submesoscale processes and dynamics // Ocean Modeling in an Eddying Regime. Geophysical monograph series. Oceanography-Mathematical models. - Washington, D.C.: American Geophysical Union. - 2008. - 177. -Р. 17 - 38. - doi:10.1029/177GM04
2. Karimova S. Spiral eddies in the Baltic, Black and Caspian seas as seen by satellite radar data // Advances in Space Research. - 2012. - 50, No. 8. - P. 1107 - 1124. -doi:10.1016/j.asr.2011.10.027
3. Zatsepin A.G., Baranov V.I., Kondrashov A.A. et al. Submesoscale eddies at the Cauasus Black Sea shelf and the mechanisms of their generation // Oceanology. - 2011. - 51, No. 4. -P. 554 - 567. - doi:10.1134/S0001437011040205
4. Demyshev S., Knysh V., Korotaev G. et al. The MyOcean Black Sea from a scientific point of view // Mercator Ocean Quarterly Newsletter. - 2010. - No. 39. - Р. 16 - 24.
5. Доценко С.Ф., Иванов В.А. Природные катастрофы Азово-Черноморского региона. - Севастополь, 2010. - 174 c. - URL: http://meteo.geofaq.ru/books/612759_D9-E3D_docenko_s_f_ivanov_v_a_prirodnye_katastrofy_azovo_chernomors.pdf (дата обращения: 20.08.2016).
6. Zalesny V.B., Gusev A.V., Moshonkin S.N. Numerical model of the hydrodynamics of the Black Sea and the Sea of Azov with variational initialization of temperature and salinity // Izv. Atmos. Ocean. Phys. - 2013. - 49, No. 6. - P. 642 - 658. - doi:10.1134/S0001433813060133
7. Ibraev R.A. A study of the sensitivity of the model of the Black Sea current dynamics to the surface boundary conditions // Oceanology. - 2001. - 41, No. 5. - P. 615 - 621.
8. Кубряков А.И. Применение технологии вложенных сеток при создании системы мониторинга гидрофизических полей в прибрежных районах Черного моря // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: ЭКОСИ-Гидрофизика, 2004. - Вып. 11. - С. 31 - 50.
9. Мизюк А.И., Коротаев Г.К., Ратнер Ю.Б. Адаптация модели NEMO для выполнения прогнозов в Черном море // Современное состояние и перспективы наращивания морского ресурсного потенциала юга России: межд. науч. конф., 15 - 18 сентября 2014 г.: тез. докл. - Севастополь: ЭКОСИ-Гидрофизика, 2014. - С. 181 - 184.
10. Demyshev S.G. A numerical model of online forecasting Black Sea currents // Izv. Atmos. Ocean. Phys. - 2012. - 48, No. 1. - P. 120 - 132. - doi:10.1134/S0001433812010021
11. Sarkisyan A.S. On some problems and results of ocean modeling // Oceanology. - 1996. - 36, No. 5. - P. 607 - 617.
12. Zalesny V.B., Ivchenko V.O. Modeling the global circulation response and the regional response of the Arctic Ocean to the external forcing anomalies // Oceanology. - 2010. - 50, No. 6. - P. 829 - 840. - doi:10.1134/S0001437010060020
13. Aiki H., Zhai X., Greatbatch R.J. Energetics of the global ocean: the role of mesoscale eddies // The Indo-Pacific Climate Variability and Predictability. Edition: World Scientific Series on Asia-Pacific Weather and Climate. - 2016. - 7, Ch. 4. - P. 109 - 134.
14. Cessi P., Pinardi N., Lyubartsev V. Energetics of semienclosed basins with two-layer flows at the strait // J. Phys. Oceanogr. - 2014. - 44, No. 3. - P. 967 - 979. - doi:10.1175/JPO-D-13-0129.1
15. Demyshev S. G. Energy of the Black Sea climatic circulation. 1. Discrete equations of the time rate of change of kinetic and potential energy // Meteorologiya i Gidrologiya. - 2004. - 9. -P. 65 - 80.
16. Demyshev S.G., Dymova O.A. Analyzing intraannual variations in the energy characteristics of circulation in the Black Sea // Izv. Atmos. Ocean. Phys. - 2016. - 52, No. 4. - P. 386 -393. - doi:10.1134/S0001433816040046
17. Mellor G.L., Yamada T. Development of a turbulence closure model for geophysical fluid problems // Rev. Geophys. Space Phys. - 1982. - 20, Issue 4. - Р. 851 - 875. -doi: 10.1029/RG020i004p00851
18. Farda A., Deue M., Somot S. et al. Model ALADIN as regional climate model for Central and Eastern Europe // Studia Geophysica et Geodaetica. - 2010. - 54, Issue 2. - P. 313 - 332. -doi:10.1007/s11200-010-0017-7
19. Гидрометеорология и гидрохимия морей СССР. Т. IV. Черное море. Вып. 1. Гидрометеорологические условия. - СПб.: Гидрометеоиздат, 1991. - 428 с.
20. Robinson A.R., Harrison D.E., Mintz Y. et al. Eddies and the general circulation of an idealized oceanic gyre: a wind and thermally driven primitive equation numerical experiment // J. Phys. Oceanogr. - 1977. - 7, Issue 2. - P. 182 - 207. - doi:http://dx.doi.org/10.1175/1520-0485(1977)007<0182:EATGC0>2.0.C0;2
21. Зацепин А.Г., Кременецкий В.В., Станичный С.В. и др. Бассейновая циркуляция и мезо-масштабная динамика Черного моря под ветровым воздействием // Современные проблемы динамики океана и атмосферы: сб. ст. - М.: ТРИАДА ЛТД, 2010. - С. 347 - 368.
22. Stanev E., Rachev N. Numerical study on the planetary Rossby modes in the Black Sea // J. Marine Syst. - 1999. - 21, Issues 1 - 4. - P. 283 - 306. - doi:http://dx.doi.org/10.1016/S0924-7963(99)00019-6
23. Коротаев Г.К., Коснырев В.К., Куфтарков Ю.М. О доступной потенциальной энергии и генерации мезомасштабных вихрей в океане // Морские гидрофизические исследования. - 1977. - 79, № 4. - С. 200 - 212.
24. Вольцингер Н.Е., Пясковский Р.В. Основные океанологические задачи теории мелкой воды. - Л.: Гидрометеоиздат, 1968. - 300 с.
25. Elkin D.N., Zatsepin A.G. Laboratory investigation of the mechanism of the periodic eddy formation behind capes in a coastal sea // Oceanology. - 2013. - 53, No. 1. - P. 24 - 35. -doi: 10.1134/S0001437012050062
Calculation and analysis of water circulation energetics in the Black Sea coastal regions
S.G. Demyshev, O.A. Dymova*
Marine Hydrophysical Institute, Russian Academy of Sciences, Sevastopol, Russia *e-mail: olgadym@yahoo.com
The results of analysis of the kinetic and potential energy budget components calculated using the MHI eddy-resolving model are represented. The numerical experiments are carried out with the spatial resolution 1.6 km under allowing for real atmospheric forcing in 2006. Temporal-spatial variability of the terms in the equations of kinetic and potential energy budget is considered on the average annual and seasonal scales for certain sea regions. It is shown that the most energy-significant budget components are the buoyancy force in the northern and northeastern parts of the Black Sea, the wind force in its southern part, and the pressure and friction forces above the continental slope in the RIM current zone. The mesoscale eddies nearby the Crimean and the North Caucasus coasts are formed mainly in a spring-summer season resulting from baroclinic instability. The basic factor conditioning mesoscale dynamics along the southern and the southeastern Black Sea coast in all the seasons and at a weak wind is the influence of the coastal features upon the RIM current periphery.
Keywords: Black Sea, numerical modeling, kinetic energy, potential energy, mesoscale eddies.
