CC BY
22
УДК 551.465
П.Н. Лишаев, Г.К. Коротаев, В.В. Кныш, А.И. Мизюк, О.А. Дымова
Восстановление синоптической изменчивости гидрофизических полей Черного моря на основе реанализа за 1980 - 1993 годы
Синоптическая изменчивость гидрофизических полей Черного моря воспроизведена по результатам реанализа, выполненного за период 1980 - 1993 гг. посредством вихреразрешающей модели МГИ НАН Украины с использованием атмосферных полей из реанализа ERA-Interim и с ассимиляцией в модели специально подготовленных среднегодовых профилей невозмущенных температуры и солености. Диаграммы межгодовой и сезонной изменчивости средней по горизонтам температуры подтверждают зависимость характеристик холодного промежуточного слоя моря от зимних термических условий. Получено, что тенденции изменчивости средних по слоям температуры (0 - 300 м) и солености (100 - 250 м) из реанализа с вихреразрешением и реанализа крупномасштабных гидрофизических полей за период 1971 - 1993 гг., наиболее обеспеченный данными измерений, качественно совпадают. Выявленная при реанализе синоптическая изменчивость в 1980 - 1993 гг. характеризуется наличием известных антициклонических вихрей справа от Основного Черноморского течения и интенсивных циклонических круговоротов. Синоптические вихревые структуры присутствуют в полях и рассчитанного и аль-тиметрического уровней, но различаются по расположению.
Ключевые слова: реанализ ERA-Interim, профили невозмущенных температуры и солености, прогностический расчет, ассимиляция профилей, синоптические вихри, сравнение с данными альтиметрии.
Введение
Реанализ крупномасштабных гидрофизических полей Черного моря, выполненный по о- и z-координатным моделям за период 1971 - 1993 гг. с ассимиляцией среднемесячных данных измерений температуры и солености и с использованием атмосферных полей из глобального реанализа ERA-40 [1], показал, что тенденции межгодовой изменчивости температуры, солености и кинетической энергии течений получены вполне реалистичными [2 - 6]. Однако восстановленные гидрофизические поля моря характеризуются слабой синоптической изменчивостью. Причиной этого являются использовавшиеся при ассимиляции поля данных измерений температуры и солености, сглаженные вследствие применения специальной методики их получения [7]. Кроме того, пространственное разрешение атмосферных полей, задававшихся в качестве граничных условий, было достаточно грубым (2,5 х 2,5°). Поэтому задача проведения реанализа за последние 33 года (1980 - 2012 гг.) с учетом модельных расчетов и наблюдающихся в структурах гидрофизических полей моря синоптических вихрей [8 - 14] становится весьма актуальной. Адекватно восстановленные гидрофизические поля, их сезонная и межгодовая изменчивость необходимы для изучения климата моря и создания трехмерных экологических моделей.
Для решения задачи реанализа с ассимиляцией в численной модели данных измерений температуры и солености наиболее сложными являются вре-
© П.Н. Лишаев, Г.К. Коротаев, В.В. Кныш, А.И. Мизюк, О.А. Дымова, 2014
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
49
менной интервал 1993 - 2001 гг. и современный период 2002 - 2012 гг. Они характеризуются редкими нерегулярными по пространству и времени измерениями температуры и солености на гидрологических станциях и буях-профилемерах Argo. Отметим, что с октября 1992 г. и по настоящее время имеются данные об аномалиях альтиметрического уровня Черного моря, которые позволили разработать методику формирования полей псевдонаблюдений температуры и солености и выполнить реанализ гидрофизических полей с использованием измерений буев Argo в 2007 - 2008 гг. [15]. Однако построенные поля псевдонаблюдений температуры и солености оказались сглаженными, поскольку при оптимальной интерполяции данных измерений в узлы сеточной области численной модели использовались автокорреляционные функции, оцененные по климатическим полям. В синоптической изменчивости этих полей присутствовали лишь некоторые антициклонические вихри справа от Основного Черноморского течения (ОЧТ). Воспроизведение синоптической изменчивости гидрофизических полей Черного моря следует получить посредством более высокого пространственного разрешения атмосферных параметров и численной модели с одновременной коррекцией стратификации вод моря данными измерений для обеспечения реалистичной глубины главного пикноклина.
В данной работе синоптическая изменчивость гидрофизических полей Черного моря воспроизводится в реанализе за период 1980 - 1993 гг. с использованием вихреразрешающей нелинейной модели МГИ НАН Украины. Атмосферные поля брались из реанализа ЕВА-1Швпт с пространственным разрешением 0,75 х 0,75°. Модельные значения температуры и солености корректировались среднегодовыми профилями этих параметров, полученными по данным реанализа за 1971 - 1993 гг. [16], что позволило воспроизвести межгодовую изменчивость полей моря и реалистичное распределение термо-халинных параметров в постоянном галоклине и главном пикноклине. Среднегодовые профили температуры и солености получены на основе модифицированной однопараметрической связи полей плотности (солености, температуры) и полей динамического уровня. Выполнен анализ сезонной и межгодовой изменчивости средних по слоям температуры, солености, кинетической энергии, а также средней по горизонтам температуры. Проведено сопоставление полей рассчитанного и альтиметрического уровней.
Используемая модель динамики морских течений и сведения о параметрах атмосферного воздействия
Для проведения реанализа гидрофизических полей за 1980 - 1993 гг. использовалась вихреразрешающая трехмерная нелинейная модель с учетом воздействия атмосферных полей из реанализа ЕЕА-1п(епт и с ассимиляцией в модели специально подготовленных среднегодовых профилей температуры и солености.
Система уравнений модели в приближении Буссинеска, гидростатики и несжимаемости морской воды записана в декартовой системе координат. Ось х направлена на восток, ось у - на север, ось г - вертикально вниз. Система уравнений имеет вид:
50
ISSN 0233-7584. Мор. гидрофиз. журн, 2014, № 5
Щ - (% + /> + wuz = — §Ро£х —— (Р' + Е)х + (нуыг)г - няV4и , (1)
Ро
+ (С + [Щ + WVz = — §Ро^у — — (Р' + Е)у + (н^)z — ня V4v , (2)
Ро
Щх + Уу + Wz = 0, (3)
н
£ + | (Щх + Уу ^ = (Рг— БУ)/Р1, (4)
Р = ЯР0С + Я = яр0С + Р', (5)
о
Т+(иТ) х + (УТ) у + (ъТ) z = —кн V 4Т + {ктТ2) z — ^, (6)
дz
+(и5) х + (у5 ) у + ) z = —кн V4 5 + (кSSz) z, (7)
Р = РТ, 5), (8)
где и, V, ъ - компоненты вектора скорости, С, - возвышение свободной поверхности, [ - параметр Кориолиса, я - ускорение свободного падения, , -время, Н - глубина моря, Р - давление, Рг - осадки, Еу - испарение, Т - температура, 5 - соленость, р - плотность морской воды, р 1 - средняя плотность в поверхностном слое, р0 = 1 г/см3, н^, нн - коэффициенты вертикальной и горизонтальной турбулентной вязкости соответственно, кН - коэффициент горизонтальной турбулентной диффузии, к - коэффициент вертикальной турбулентной диффузии тепла, к5 - коэффициент вертикальной турбулентной
диффузии соли, % = ду — дщ, е = р0 -—. Уравнение (4) получено в пред-
дх ду 2
положении выполнения линеаризованного кинематического условия в виде
ъ = + —Г—— . Последнее слагаемое в правой части (6) описывает изме-Р1
нение температуры с учетом влияния коротковолновой радиации [17]. Здесь I(z) = ехр(—z /Ь1) + (1 — а) ехр(—z /Ь2)], - коротковолновая радиация на поверхности моря, а, Ь ь Ь 2 - эмпирические константы, которые подбираются на основе численных расчетов.
Для описания верхнего перемешанного слоя используется приближение Пакановски - Филандера [18]. Коэффициенты турбулентного обмена импульсом и диффузии тепла и соли по вертикали вычисляются следующим образом:
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
51
ну = Но(1 + И!)-2 + нГ,
К = [н,(1 + И1)-2 + н1]/(1 + И1) + К , (9)
кт = [Но(1 + Ш)-2 + н1]/(1 + Ш) + К ,
£ др 1 V Т 5 где И =------ - число Ричардсона; н„, н, , н, , к , к
р0 дz (ди / дг)2 + (^ / дг)2 -01111
- заданные константы.
Для системы уравнений (1) - (8) краевые условия на поверхности (г = 0) имеют вид:
X V
нуиг =-Т , =-Т ,
т т 1 с ЕУ - Рг 1 (10)
КТТг = ОТ + в(Т ° - То), =-5о + в(5С - 5>) .
р1
Здесь использованы следующие обозначения: тх, ту - компоненты вектора касательного напряжения трения ветра, О - поток тепла за вычетом коротковолновой радиации, Т0, 50 - поверхностные температура и соленость, Тс1, 5с1 - климатические температура и соленость, в - параметр релаксации. На дне (г = Н(х, у)):
и = 0, V = 0, М> = 0, Тг = 0, 5 = 0. (1 1)
На твердых боковых стенках для меридиональных участков границы:
и = 0, V 2и = 0, vx = 0, = 0,
Тх = 0, ^2Т) х = 0, Бх = 0, (V2?) х = 0,
(12)
для зональных участков границы:
V = 0, V2v = 0, иу = 0, V 2иу = 0, Ту = 0, ^2Т) у = 0, 5у = 0, (V2 5) у = 0.
(13)
На участках границы, где вода втекает, используются условия Дирихле для меридиональных участков:
и = ир, V2u = 0, vx = 0, V2vx = 0,
Т = Тр, 5 = 5р, (V2Т)х = 0, (V25)х = 0, для зональных участков:
V = V13, V^ = 0, иу = 0, V2иу = 0, Т = Тр, 5 = 5р, (V2Т)у = 0, (V25)у = 0.
(14)
(15)
52 0233-7584. Мор. гидрофиз. журн, 2014, № 5
Для верхнебосфорского течения и Керченского пролива, когда течение направлено из Черного моря в Азовское:
V = vn, V2v = 0, uv = 0, V2uv = 0,
' 2 ' 2 (16)
Ty = 0, Sy = 0, (V T)y = 0, (V2S)y = 0.
В соотношениях (14) - (16) введены следующие обозначения: up, vp, vn -
скорость в устьях рек и проливах соответственно; Tp, Sp - температура и соленость в реках.
В качестве начальных при t = t0 потребуем выполнения таких условий:
u = u0 ( x, y, z), v = v0(x, y, z), Z = Z0(x, y), (17)
T = T 0( x, y, z), S = S 0( x, y, z).
Таким образом, решается система уравнений (1) - (8) с соответствующими краевыми (10) - (16) и начальными (17) условиями. Конечно-разностная дискретизация уравнений модели, начальных и граничных условий проведена на сетке С [19]. Дискретная модель при равномерном шаге обладает вторым порядком аппроксимации по пространству и первым (в результате периодического применения схемы Мацуно [20]) - по времени. Используемые разностные операторы и особенности аппроксимации уравнений модели подробно описаны в работе [21].
Расчеты проводились на горизонтальной сетке 5 х 5 км (238 х 132 точки). По вертикали использовались 38 горизонтов: 2,5, 5, 10, ..., 30, 40, 50, 63, 75, 88, 100, 113, 125, 150, ..., 200, 250, 300, 400, ..., 2100 м. Шаг по времени равнялся 5 мин.
Значения эмпирических констант для учета коротковолновой радиации следующие: а = 0,85; b 1 = 15,5; b2 = 2. В соотношениях (9) задано, что
2 2 V 2 T 2 S 2
н0 = 30 см /с; н1 = 1 см /с; н1 = 5 см /с; к1 = 5 см /с; к1 = 1 см /с. В формуле (10) параметр релаксации в = 1,15 х 10-3см/с. Для задания температуры воды и скорости течения в устьях рек и проливах использовались литературные данные [22]. В зоне смешения речных и морских вод соленость считалась равной 7%о. В расчете полагалось, что в верхнебосфорском течении температура и соленость те же, что и в море. В нижнебосфорском потоке соленость принималась равной 22%о, а температура 16°С, что соответствует характеристикам вод Мраморного моря. Принималось также, что верхнее течение в Босфоре распространяется до глубины 27,5 м, нижнее - до 68,75 м. Коэффициенты горизонтального турбулентного обмена импульсом и горизонтальной турбулентной диффузии тепла и соли брались следующими: нн = 5х1017 см4/с, кH = 1016 см4/с.
На поверхности моря задавались потоки тепла, осадков и испарения, полученные по данным атмосферного реанализа ERA-Interim с разрешением по пространству 0,75° через каждые 12 ч. Поля тангенциального напряжения трения ветра имелись на каждые 6 ч. Граничные поля линейно интерполировались на каждый шаг по времени.
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5 53
Анализ результатов реанализа гидрофизических полей за период 1980 - 1993 гг. на основе вихреразрешающей модели
В качестве начальных полей в вихреразрешающей модели принимались климатические поля температуры, солености, уровня и скоростей течений, полученные в статье [23] и интерполированные на сеточную область модели. Начальный момент времени соответствовал 1 января 1980 г. С использованием атмосферных полей из реанализа ЕЕЛ-1п(впт было проведено три численных эксперимента. В первом из них выполнялся прогностический расчет с учетом реального атмосферного воздействия. Второй эксперимент (далее реанализ за период 1980 - 1993 гг.) отличался от первого тем, что на каждом шаге по времени для области моря с глубинами, равными и большими 500 м, производилась коррекция модельных значений, например, температуры по формуле
Т + (х, у, ^ 0 = Т~ (х, у, ^ 0 + ср(^ ,г) — Т (г, 0], (18)
320 п(z)
где Т (z,,) - среднее по горизонту значение модельной температуры;
Тср(z, ,г) - среднегодовое значение невозмущенной температуры; п(z) - количество точек на горизонте. Верхние индексы «-», «+» соответствуют значениям температуры до и после коррекции. Заметим, что аналогичная формула использовалась для коррекции солености.
Весовой коэффициент в последнем слагаемом правой части (18) подобран в результате численных экспериментов. Среднегодовые профили невозмущенных температуры и солености для каждого года периода 1980 -1993 гг. рассчитывались по данным реанализа за 1971 - 1993 гг. [16]. С этой целью были построены диаграммы эмпирических зависимостей «среднегодовая соленость (температура) - среднегодовой уровень» в градациях уровня с дискретностью 1 см. В пределах каждой градации вычислялись средние значения температуры и солености за год. Расчеты проводились для области моря, ограниченной изобатой 500 м. Среднегодовые профили невозмущенных температуры и солености (плотности) соответствовали «нулевой» градации среднегодового уровня моря (-1, 0, +1), т. е. адиабатическому процессу в море.
В третьем численном эксперименте восстановление гидрофизических полей моря выполнялось только для 1993 г. таким же способом, как во втором эксперименте, но с увеличением значений компонентов вектора тангенциального напряжения трения ветра в 1,5 раза.
Данные реанализов по вихреразрешающей модели получены с дискретностью по времени 1 сут и являются исходными для исследования сезонной и межгодовой изменчивости среднесуточных температуры, солености, скорости течений и кинетической энергии (КЭ), а также синоптических вихрей.
Обратимся к рассмотрению некоторых интегральных характеристик. Профили солености, полученные посредством осреднения полей по горизонтам для зимы и лета, приведены на рис. 1. Зимой хорошо выражен верхний квазиоднородный слой (ВКС), особенно в прогностическом расчете (рис. 1, а). Положение по глубине постоянного галоклина несколько
54
ISSN 0233-7584. Мор. гидрофиз. журн, 2014, № 5
выше по данным реанализа за 1971 - 1993 гг. В прогностическом расчете зимой и летом слой постоянного галоклина расположен на большей глубине. Галоклин менее обострен по сравнению с данными второго эксперимента и реанализа за 1971 - 1993 гг. Рис. 1, б свидетельствует о том, что летом глубина ВКС по данным численных экспериментов равна ~10 м. Обострение постоянного галоклина максимально по результатам реанализа за 1971 - 1993 гг. Отметим, что профили солености зимой и летом по данным второго численного эксперимента и реанализа за 1971 - 1993 гг. качественно близки между собой.
□
1984 1985 Ассимиляция среднегодовых профилей Т, 5 (1980- 1993)
Р и с. 1. Профили солености 15 февраля (а) и 15 августа (б) в 1984 и 1985 гг.
Показательными являются восстановленные для 1984 и 1985 гг. профили температуры зимой (рис. 2, а) и летом (рис. 2, б, в). Зимой в слое 0 -50 м температура воды по данным реанализа за период 1980 - 1993 гг. существенно выше, чем в реанализе за 1971 - 1993 гг. (рис. 2, а). Таким образом, ассимиляция в модели среднегодовых профилей невозмущенной температуры завышает температуру верхнего слоя моря. Профили температуры, полученные по данным прогностического расчета, расположены ближе к профилям по данным реанализа за 1971 - 1993 гг. В холодную зиму 1985 г., в отличие от теплой зимы 1984 г., температура в слое 0 - 50 м меньше 8°С во всех трех рассматриваемых случаях.
В августе в те же годы сезонный термоклин наиболее ярко выражен по данным реанализа за 1971 - 1993 гг. (рис. 2, б). Ассимиляция в модели среднегодовых профилей невозмущенной температуры приводит к заниже-ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5 55
нию глубины ВКС летом. На рис. 2, в видно, что значения температуры воды в холодном промежуточном слое (ХПС), восстановленные во втором эксперименте, ближе к таковым по данным реанализа за 1971 - 1993 гг. Распределение профилей температуры по годам (рис. 2, в) подтверждает зависимость характеристик ХПС от зимних термических условий в соответствии с работами [24 - 26].
6 6.5 7 7.5 8 8.5 9 5 10 15 20 25 6.5 7 7.5 8 8.5 9
Температура, "С Температура, "С Температура, °С
С1 б в
О • □ ■ ДА
1984 1985 1984 1985 1984 1985
I Грогностический Ассимиляция Реанализ
расчет среднегодовых (1971 - 1993) (1980 1993) профилей Т, Б (1980 1993)
Р и с. 2. Профили температуры 15 февраля (а) и 15 августа (б, в) в 1984 и 1985 гг.
Поведение со временем средней в слое 0 - 300 м температуры подтверждает, что она выше по данным реанализа с коррекцией модельных значений температуры (рис. 3). Важно отметить, что тенденции изменчивости средней в слое температуры в обоих реанализах качественно совпадают. Линейные тренды являются отрицательными (рис. 3). Характеристики трендов следующие: -5,51Е-6 °С/сут в реанализе за 1980 - 1993 гг. и -1,47Е-5 °С/сут в реанализе крупномасштабных полей моря за 1971 - 1993 гг.
Диаграмма межгодовой и сезонной изменчивости средней по горизонтам температуры в слое 0 - 300 м (рис. 4) подтверждает завышение температуры в слое 0 - 40 м в осенне-зимний период (см. рис. 7 из статьи [16]). На ней также четко видна межгодовая изменчивость. Выделяются аномально теплая (1981 г.), теплая (1984 г.), нормальные (1987, 1989, 1991 гг.) и холодные (1985, 1992, 1993 гг.) зимы [24, 26]. На основе полученных результатов можно сделать вывод о том, что при выполнении реанализа по предложенной методике необходим учет сезонной изменчивости ассимилируемых в модели профилей температуры.
56
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
9.5
1980 1982 1984 1986 1988 1990 1992 1994
Р и с. 3. Сезонная и межгодовая изменчивость средней в слое 0 - 300 м температуры: кривая 1 - по данным реанализа за 1980 - 1993 гг. с ассимиляцией среднегодовых профилей температуры и солености; кривая 2 - по данным реанализа крупномасштабных гидрофизических полей за 1971 - 1993 гг. [1] (прямые линии здесь и далее - линейные тренды)
Р и с. 4. Диаграмма межгодовой и сезонной изменчивости средней по горизонтам температуры в слое 0 - 300 м по данным реанализа за 1980 - 1993 гг.
Анализ тенденций сезонной и межгодовой изменчивости средней в слое 100 - 250 м солености выявил следующее. Линейные тренды по данным обоих реанализов качественно одинаковы (положительные) (рис. 5, а). По реана-лизу за 1980 - 1993 гг. увеличение средней солености равно 7,25Е-7 %о/сут (кривая 1, рис. 5, а), по результатам реанализа за 1971 - 1993 гг. характеристика тренда равна 4,27Е-6 %о/сут (кривая 2, рис. 5, а). В верхнем слое 0 -40 м отрицательные линейные тренды по данным обоих реанализов качественно соответствуют друг другу (рис. 5, б). В расчете с коррекцией солености характеристика тренда равна -3,19Е-6 %о/сут (кривая 1, рис. 5, б), в реанализе крупномасштабных полей уменьшение солености равно -1,73Е-6 %о/сут (кривая 2, рис. 5, б). Тем самым по данным реанализа за 1980 - 1993 гг. подтверждается распреснение вод верхнего слоя моря.
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5 57
Р и с. 5. Сезонная и межгодовая изменчивость средней в слоях 100 - 250 м (а) и 0 - 40 м (б) солености: кривая 1 - по данным реанализа за 1980 - 1993 гг. с ассимиляцией среднегодовых профилей температуры и солености; кривая 2 - по данным реанализа крупномасштабных гидрофизических полей за 1971 - 1993 гг.
Пониженные значения средней по слоям солености в реанализе за 1980 -1993 гг. объясняются ассимилируемыми в модели профилями среднегодовых значений солености «нулевой» градации, заниженных по сравнению со значениями солености отрицательных градаций среднегодового уровня моря. Необходимы дальнейшие исследования, направленные на выбор оптимальных профилей невозмущенной солености.
На приведенном на рис. 6, б графике (кривая 1) четко прослеживается сезонная изменчивость КЭ в течение года: летом и осенью она меньше, чем зимой. Наибольшая интенсивность циркуляции вод в слое 0 - 300 м наблюдается в феврале - марте, наименьшая — в сентябре - октябре. Особенности сезонной и межгодовой изменчивости среднего по поверхности моря модуля вектора касательного напряжения трения ветра (кривая 1, рис. 6, а) обуславливают сезонную и межгодовую изменчивость величин средней по горизонтам КЭ. На рис.6, б видно, что повышенные значения средней в слое 0 - 300 м КЭ наблюдаются в 1988, 1992 и 1993 гг. Подчеркнем, что существенный ее рост виден в 1992 и 1993 гг. Межгодовая изменчивость средней в слое 0 - 300 м КЭ характеризуется хорошо выраженным положительным линейным трендом.
Сезонная и межгодовая изменчивость средней в слое 0 - 300 м КЭ (кривая 2, рис. 6, б) заметно отличается от полученной нами. Она характеризуется, во-первых, весьма слабым положительным линейным трендом и, во-вторых, меньшей интенсивностью циркуляции вод. Причина этого состоит в том, что при выполнении реанализа крупномасштабных гидрофизических полей [16] использовались поля из реанализа ERA-40, в котором компоненты вектора касательного напряжения трения ветра оказались заниженными (кривая 2, рис. 6, а). Данные по атмосферным параметрам реанализа ERA-Interim в этом отношении являются более предпочтительными.
58 ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
Р и с. 6. Сезонная и межгодовая изменчивость: а - модуля вектора касательного напряжения трения ветра (кривая 1 - по данным ЕКЛ-1п1епт, кривая 2 - по данным ЕЯЛ-40); б - средней в слое 0 - 300 м кинетической энергии (кривая 1 - по данным реанализа за 1980 - 1993 гг. с ассимиляцией среднегодовых профилей температуры и солености, кривая 2 - по данным реанализа крупномасштабных гидрофизических полей за 1971 - 1993 гг.)
Проведем анализ воспроизводимости вихреразрешающей моделью наблюдаемых синоптических вихреобразований. Основное внимание обратим на вихри, формирующиеся как справа от ОЧТ, так и в восточном и западном циклонических круговоротах. Для анализа используем данные восстановления гидрофизических полей в отмеченных выше трех численных экспериментах в нормальном (1987) и холодном (1993) годах по шкале зимних термических условий [26].
Топография уровня моря 25 апреля 1987 г. в прогностическом расчете отличается от полученной в численном эксперименте с коррекцией модельных профилей соответствующими среднегодовыми профилями невозмущенных температуры и солености (рис. 7, а, в). В прогностическом расчете (рис. 7, а) вследствие взаимодействия двух антициклонических вихрей - Си-нопского и Кавказского - возникают два циклонических круговорота. Центрально-западный циклонический круговорот включает в себя два циклони-
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
59
ческих вихря. Присутствуют также Батумский и Севастопольский антициклоны. Перепад уровня моря в прогностическом расчете составляет 35 см.
Р и с. 7. Топография уровня моря (см) и течения (см/с) на глубине 50 м, восстановленные на 25.04.1987 г. в прогностическом расчете (а, б) и с ассимиляцией в модели среднегодовых профилей температуры и солености (в, г)
60
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
На рис. 7, в, построенном по результатам реанализа, видно, что на указанную дату охватывающий все море циклонический круговорот состоит из интенсивных западного и восточного циклонических круговоротов и юго-восточного циклонического вихря меньшей интенсивности. Справа от ОЧТ наблюдается большее количество синоптических антициклонических вихрей: Босфорский, Сакарьякский, Синопский, Кизилирмакский, Батумский, два Кавказских и Севастопольский. Перепад уровня в этом расчете составляет 31 см, однако синоптическая изменчивость гидрофизических полей Черного моря справа от ОЧТ более выражена в сравнении с прогностическим расчетом.
Вихреобразования, проявляющиеся в поле уровня, четко прослеживаются на картах течений (рис. 7, б, г). Течения на глубине 50 м в центрально-западном циклоническом круговороте прогностического расчета более интенсивны, чем в ОЧТ по данным реанализа за 1980 - 1993 гг. Максимальные значения модуля вектора скорости течений равны 53,9 и 48,5 см/с соответственно. В то же время интенсивность течений в большинстве антициклонов справа от ОЧТ на рис. 7, г выше, чем в антициклонах на рис. 7, б.
Топография уровня моря, восстановленная на 05.05.1993 г. (рис. 8), показывает, что ОЧТ, восточный и западный циклонические круговороты и синоптическая изменчивость справа от ОЧТ ярко выражены как по данным первого, так и второго численных экспериментов. Однако количество вихрей больше и интенсивность антициклонической циркуляции выше в реанализе с коррекцией модельных значений температуры и солености их среднегодовыми профилями. Перепад уровня составляет 46 см в прогностическом расчете и 53 см - в реанализе. Обращает на себя внимание полученный в этом реана-лизе циклон, находящийся на периферии Батумского антициклона в юго-восточном районе Черного моря. Отметим, что Синопский, Батумский, Кавказский и Севастопольский антициклоны, воспроизводимые в обоих экспериментах, образуются в различные моменты времени.
Интенсивность ОЧТ, восстановленная в реанализе, выше, чем в прогностическом расчете (рис. 8, б, г). Максимальные значения модуля вектора скорости течений, полученные в этих численных экспериментах, равны 51,2 и 45,1 см/с соответственно. Интенсивность циркуляции в антициклонических структурах справа от ОЧТ выше на карте рис. 8, г, чем на карте рис. 8, б.
Ценной информацией для качественной и количественной оценки восстановленных в численных экспериментах полей уровня моря являются карты динамического уровня, построенные по данным альтиметрии. Отметим, что поля альтиметрического динамического уровня Черного моря впервые были получены в работе [27]. В данной работе для построения таких полей использовался региональный массив аномалий уровня Черного моря AVISO (Archivage Validation Interprétation des données des Satellites Océanographiques) для периода с октября 1992 г. по 2011 г. (http://www.aviso.oceanobs.com/). Эти данные были получены по комбинированным измерениям вдольтрековых аномалий уровня с нескольких спутников с помощью методов, предложенных в статьях [28, 29]. Пространственное разрешение карт аномалий уровня составляет 1/8° (~12,5 км), что в два раза выше, чем для глобального массива данных (1/4°). Временное разрешение составляет 7 дней (delayed-time product). Поля альтиметрического динамическо-
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5 61
го уровня моря получены посредством суммирования значений аномалий уровня со значениями среднегодового климатического уровня, восстановленного по модели с ассимиляцией климатических значений температуры и солености [30].
■20 -10 0 10 20 см
■20 -10 0 10 20 см
Р и с. 8. Топография уровня моря (см) и течения (см/с) на глубине 50 м, восстановленные на 05.05.1993 г. в прогностическом расчете (а, б) и с ассимиляцией в модели среднегодовых профилей температуры и солености (в, г)
62 ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
20 -10 0 10 20 см
20 -10 0 10 20 см
Р и с. 9. Топография уровня моря (см) (а) и течения (см/с) на глубине 50 м (б), восстановленные на 05.05.1993 г. в реанализе за 1980 - 1993 гг. с увеличением компонентов вектора касательного напряжения трения ветра в 1,5 раза в 1993 г., а также динамический уровень моря (см) на 05.05.1993 г. по данным альтиметрических измерений (в)
Сопоставление полученной в реанализе карты рельефа свободной поверхности на 05.05.1993 г. (рис. 8, в) с картой альтиметрического динамического уровня (рис. 9, в) свидетельствует о невысокой корреляции между ними. Анализ полей уровня и течений, восстановленных в расчете с увеличенными в 1,5 раза компонентами вектора тангенциального напряжения трения ветра (рис. 9, а, б), показал следующее. Интенсивность ОЧТ, а также циркуляции в антициклонах справа от него оказались заметно выше, чем без увеличения компонентов вектора (рис. 8, в, г). Перепад уровня равен 55 см. Количество антициклонов увеличилось. Структуры полей альтиметрического уровня и уровня, воспроизведенного в реанализе, близки между собой (рис. 9, а, в). Между антициклоническими вихрями наблюдается расхождение
lSSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
63
по их расположению. Интенсивность поверхностных геострофических течений градиентной природы в некоторых антициклонах альтиметрического уровня выше, чем интенсивность аналогичной циркуляции в антициклонах уровня, восстановленного в модели (рис. 9, а, в).
Р и с. 10. Топография уровня моря (см), восстановленного в реанализе за 1980 - 1993 гг. с увеличением компонентов вектора касательного напряжения трения ветра в 1,5 раза в 1993 г. (а, в), и альтиметрический динамический уровень моря (см) (б, г)
64
ISSN 0233-7584. Мор. гидрофиз. журн, 2014, № 5
Как показывает анализ, синоптическая изменчивость гидрофизических полей ярко проявляется тогда, когда интенсифицируется ОЧТ. К лету ОЧТ имеет тенденцию к разрушению (рис. 10, а, б). Видно, что интенсивные синоптические вихри справа от ОЧТ, восстановленные в модели, присутствуют также в поле альтиметрического динамического уровня (рис. 10, а, б). Понижение уровня в северо-западной области моря коррелирует с таковым на карте альтиметрического уровня. Различие состоит в том, что ОЧТ, восстановленное в модели, не разрушилось. Интенсивность течений в Батумском антициклоне существенно выше по альтиметрическим данным. По этим данным наблюдается много мезомасштабных вихревых структур в западном и восточном циклонических круговоротах (рис. 10, б).
Общие структуры в полях восстановленного в модели и альтиметриче-ского уровней моря наблюдаются в зимний период (рис. 10, в, г). Весьма близки структуры полей в юго-восточной области. ОЧТ ярко выражено в обоих случаях. Севастопольский антициклон на карте альтиметрического уровня выражен слабо.
Выводы
Основные результаты численных экспериментов состоят в следующем. Синоптическая изменчивость гидрофизических полей Черного моря восстановлена в реанализе за 1980 - 1993 гг., который выполнен с использованием вихреразрешающей модели МГИ НАН Украины, атмосферных полей из реа-нализа ERЛ-Interim с пространственным разрешением 0,75° и с ассимиляцией среднегодовых профилей невозмущенных температуры и солености.
В слое 0 - 300 м четко проявляется межгодовая изменчивость средней по горизонтам температуры. Тенденции изменчивости средней в слое 0 - 300 м температуры качественно совпадают с таковыми по результатам реанализа за 1971 - 1993 гг., наиболее обеспеченные данными наблюдений.
Линейные тренды средней солености в слое 100 - 250 м по данным реа-нализа с ассимиляцией среднегодовых профилей и реанализа за 1971 -1993 гг. положительны и качественно одинаковы. Данные реанализа за 1980 -1993 гг. подтверждают распреснение вод верхнего слоя моря.
Синоптическая изменчивость гидрофизических полей ярко выражена тогда, когда интенсифицируется ОЧТ. Летом и осенью ОЧТ имеет тенденцию к разрушению.
Структуры полей альтиметрического уровня и уровня моря, воспроизведенного в реанализе за 1980 - 1993 гг., близки между собой, однако положения синоптических вихрей справа от ОЧТ не совпадают. В период разрушения ОЧТ поверхностная геострофическая циркуляция вод моря, полученная
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
65
по альтиметрическим данным, и циркуляция, восстановленная в реанализе, различаются между собой в наибольшей степени.
При воспроизведении синоптической изменчивости гидрофизических полей по предложенной методике необходимо учитывать сезонную изменчивость ассимилируемых в модели профилей невозмущенных температуры и солености.
Авторы выражают благодарность С.Г. Демышеву за консультации по вопросам модификации численной модели, А. А. Кубрякову, любезно предоставившему картированные аномалии уровня Черного моря, а также В.Л. Дорофееву, предложившему соотношение для коррекции модельных температуры и солености.
СПИСОК ЛИТЕРАТУРЫ
1. Кныш В.В., Коротаев Г.К., Моисеенко В.А. и др. Сезонная и межгодовая изменчивость гидрофизических полей Черного моря, восстановленных на основе реанализа за период 1971 - 1993 гг. // Изв. РАН. Физика атмосферы и океана. - 2011. - 47, № 3. - С. 433 -446.
2. Полонский А.Б., Ловенкова Е.А. Тренд температуры и солености деятельного слоя в Черном море во второй половине XX века и его возможные причины // Там же. - 2004. -40, № 6. - С. 832 - 841.
3. Белокопытов В.Н., Шокурова И.Г. Оценки междесятилетней изменчивости температуры и солености в Черном море в период 1951 - 1995 гг. // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: НПЦ «ЭКОСИ-Гидрофизика». - 2005. - Вып. 12. - С. 12 - 21.
4. Oguz T., Dippner J.W., Kaymaz Z. Climatic regulation of the Black Sea hydrometeorological and properties at interannual-to decadal time scales // J. Mar. Syst. - 2006. - 60, № 3 - 4. -P. 235 - 254.
5. Полонский А.Б., Шокурова И.Г, Белокопытов В.Н. Десятилетняя изменчивость температуры и солености в Черном море // Морской гидрофизический журнал. - 2013.-№ 6. - С. 27 - 40.
6. Полонский А.Б., Шокурова И.Г. Декадная изменчивость характеристик пикноклина и геострофической циркуляции вод Черного моря в зимний период // Метеорология и гидрология. - 2009. - № 4. - С. 75 - 92.
7. Моисеенко В.А., Белокопытов В.Н. Оценка качества массива данных гидрологических измерений, подготовленного для решения задачи реанализа состояния Черного моря за период 1985 - 1994 гг. // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: НПЦ «ЭКОСИ-Гидро-физика». - 2008. - Вып. 16. - С. 184 - 189.
8. Oguz T., Besiktepe S. Observations on the Rim Current structure, CIW formation and transport in the western Black Sea // Deep-Sea Res. Part I. - 1999. - 46, № 10. - P. 1733 -1753.
9. Sokolova E., Stanev E. V., Yakubenko V. et al. Synoptic variability in the Black Sea. Analysis of hydrographic survey and altimeter data //J. Mar. Syst. - 2001. - 31, № 1 - 3. - P. 45 - 63.
10. Staneva J.V., Dietrich D.E., Stanev E.V. et al. Rim current and coastal eddy mechanisms in an eddy-resolving Black Sea general circulation model // Ibid. - 2001. - 31, № 1 - 3. - P. 137 -157.
66 ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
11. Zatsepin A.G., Ginzburg A.I., Kostianoy A.G. et al. Observations of Black Sea mesoscale eddies and associated horizontal mixing // J. Geophys. Res. - 2003. - 108, № C8. - P. 3246 -3272.
12. Журбас В.М., Зацепин А.Г., Григорьева Ю.В. и др. Циркуляция вод и характеристики разномасштабных течений в верхнем слое Черного моря по дрифтерным данным // Океанология. - 2004. - 44, № 1. - С. 34 - 48.
13. Enriquez C.E., Shapiro G.I., Souza A.J. et al. Hydrodynamic modeling of mesoscale eddies in the Black Sea // Ocean Dyn. - 2005. - 55. - P. 476 - 489.
14. Джиганшин Г.Ф., Полонский А.Б., Толстошеее А.П. Исследование синоптических и инерционных движений в поле Основного Черноморского течения по данным дриф-терных измерений // Морской гидрофизический журнал. - 2013.- № 2. - С. 3 - 25.
15. Мизюк А.И., Кныш В.В., Коротаее Г.К. Методика восстановления гидрофизических полей Черного моря в период 2007 - 2008 гг. на основе ассимиляции в модели данных измерений буев-профилемеров // Экологическая безопасность прибрежной и шельфо-вой зон и комплексное использование ресурсов шельфа. - Севастополь: НПЦ «ЭКОСИ-Гидрофизика». - 2011. - Вып. 24. - С. 38 - 55.
16. Кныш В.В., Коротаее Г.К., Мизюк А.И. и др. Усвоение гидрологических наблюдений для расчета течений в морях и океанах // Изв. РАН. Физика атмосферы и океана. -2012. - 48, № 1. - С. 67 - 85.
17. Paulson C.A., Simpson J.J. Irradiance measurements in the upper ocean // J. Phys. Oceanogr. - 1977. - 7, № 6. - P. 952 - 956.
18. Pacanowsci R.C., Philander S.G.H. Parameterization of vertical mixing in numerical models of tropical oceans // Ibid. - 1981. - 11, № 11. - P. 1443 - 1451.
19. Arakawa A. Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow // J. Comput. Phys. - 1966. - № 1. -P. 119 - 143.
20. Matsuno T. Numerical integrations of the primitive equations by a simulated backward difference method // J. Meteor. Soc. Japan. Ser. 2. - 1966. - 44, № 1. - P. 76 - 84.
21. Демышее С.Г. Численная модель оперативного прогноза течений в Черном море // Изв. РАН. Физика атмосферы и океана. - 2012. - 48, № 1. - С. 137 - 149.
22. Гидрометеорология и гидрохимия морей СССР. Т. IV. Черное море. Вып. 1. Гидрометеорологические условия. - СПб.: Гидрометеоиздат, 1991. - 429 с.
23. Кныш В.В., Коротаее Г.К., Демышее С.Г. и др. Долговременные изменения термоха-линных и динамических характеристик Черного моря по климатическим данным температуры и солености и их ассимиляции в модели // Морской гидрофизический журнал. - 2005. - № 3. - С. 11 - 30.
24. Титое В.Б. Формирование зимней гидрологической структуры Черного моря в зависимости от суровости зим // Океанология. - 2000. - 40, № 6. - С. 826 - 832.
25. Титое В.Б. Об оценке температурного режима атмосферы, формирующего гидрологическую структуру Черного моря // Метеорология и гидрология. - 2000. - № 10. - С. 78 -84.
26. Титое В.Б. Влияние многолетней изменчивости климатических условий на гидрологическую структуру и межгодовое обновление холодного промежуточного слоя в Черном море // Океанология. - 2003. - 43, № 2. - С. 176 - 184.
27. Коротаее Г.К., Саенко О.А., Коблински Ч. Дж. и др. Оценка точности, методика и некоторые результаты усвоения альтиметрических данных TOPEX/POSEIDON в модели общей циркуляции Черного моря // Исследование Земли из космоса. - 1998. - № 3. -С. 3 - 17.
ISSN 0233-7584. Мор. гидрофиз. журн., 2014, № 5
67
28. Le Traon P.-Y., Dibarboure G., Ducet N. Use of a high-resolution model to analyze the mapping capabilities of multiple-altimeter missions // J. Atmos. Ocean. Technol. - 2001. -№ 18. - P. 1277 - 1288.
29. Pascual A., Faugere Y., Larnicol G. et al. Improved description of the ocean mesoscale variability by combining four satellite altimeters // Geophys. Res. Lett. - 2006. - № 33. - L02611. - P. 1 - 16.
30. Кныш В.В., Кубряков А.И., Инюшина Н.В. и др. Восстановление климатической сезонной циркуляции Черного моря на основе модели в с-координатах с использованием ассимиляции данных о температуре и солености // Экологическая безопасность прибрежной и шельфовой зон и комплексное использование ресурсов шельфа. - Севастополь: НПЦ «ЭКОСИ-Гидрофизика». - 2008. - Вып. 16. - С. 243 - 265.
Морской гидрофизический институт НАН Украины, Материал поступил
Севастополь в редакцию 29.01.14
E-mail: [email protected] После доработки 03.03.14
АНОТАЦ1Я Синоптична мшлишсть гiдрофiзичних полш Чорного моря вщгворена за результатами реаналiзу, виконаного за перюд 1980 - 1993 рр. за допомогою вихоророздшювально! мо-делi МП НАН Украши з використанням атмосферних полiв з реаналiзу ЕКЛ-1п1ег1т i з асимь лящею в моделi спещально тдготовлених середньорiчних профг™ незбурених температури та солоностг Дiаграми мiжрiчноi i сезонно! мшливосп середньо! за горизонтами температури тдтверджують залежнiсть характеристик холодного промiжного шару моря вiд зимових тер-мiчних умов. Отримано, що тенденци мiнливостi середнiх за шарами температури (0 - 300 м) та солоност (100 - 250 м) з реаналiзу з вихоророздшьшстю та реаналiзу великомасштабних гiдрофiзичних полiв за перюд 1971 - 1993 рр., найбшьш забезпечений даними вимiрювань, яюсно спiвпадають. Виявлена при реаналiзi синоптична мшливють в 1980 - 1993 рр. характеризуешься наявнiстю ввдомих антициклонiчних вихорiв праворуч вiд Основно! Чорноморсько! течii та iнтенсивних циклошчних кругообiгiв. Синоптичнi вихровi структури присутт в полях розрахованого та альтиметричного рiвнiв моря, але вiдрiзняються за розташуванням.
Ключовi слова: реаналiз ЕКЛ-1п1ег1т, профш незбурених температури та солоносл, про-гностичний розрахунок, асимiляцiя профЫв, синоптичнi вихори, порiвняння з даними альти-метрii.
ABSTRACT Mesoscale variability of the Black Sea hydrophysical fields is reproduced based on the results of reanalysis for 1980 - 1993 carried out by the eddy-resolving model of MHI NAS of Ukraine. The model is run using the atmospheric fields from ERA-Interim reanalysis and assimilation of specially prepared annual-average profiles of undisturbed temperature and salinity. Diagrams of inter-annual and seasonal variability of average over horizons temperature confirm dependence of cold intermediate layer characteristics on winter thermal conditions. It is found that variability trends of average over layers temperature (0 - 300 m) and salinity (100 - 250 m) derived from the eddy-resolving reanalysis and those obtained from reanalysis of large-scale hydrophysical fields for 1971 -1993 (mostly provided with measurements) are in qualitative agreement. Mesoscale variability revealed in reanalysis for 1980 - 1993 is characterized by presence of known anticyclone eddies to the right of the Rim current and intense cyclonic gyres. Mesoscale vortex structures are present both in the simulated sea levels and the altimetry ones, but they differ in location.
Keywords: ERA-Interim reanalysis, profiles of undisturbed temperature and salinity, prognostic calculation, assimilation of profiles, mesoscale eddies, comparison with altimetry data.
68
ISSN 0233-7584. Мор. гидрофиз. журн, 2014, № 5
