Климатическая изменчивость циркуляции вод Японского моря во второй половине XX века

Вестник ДВО РАН. 2015. № 6

УДК 551.465.43(265.53) Д.В. СТЕПАНОВ

Климатическая изменчивость циркуляции вод Японского моря во второй половине XX века

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

Интенсификация и ослабление циклонического круговорота имеют периодичность 3, 4 и 5 лет. На декадных масштабах выявлена значительная интенсификация циклонического круговорота в интервале с 1948 по 1965 г., после которого вплоть до начала 2000 г. наблюдается период его значительного ослабления. На межгодовых масштабах основной вклад в изменчивость вихревой составляющей циклонического круговорота вносит положительный вихрь напряжения ветра (ВНВ), а на декадных — положительный ВНВ и меридиональный градиент ВНВ над Японским морем. Ослабление плотностной стратификации над Центральной котловиной с 1958 по 1983 г. обусловлено увеличением потока явного тепла в атмосферу в районе вдольбереговой ветви Цусимского течения. Ослабление меридионального градиента потока явного тепла над морем в масштабе 12—15 лет приводит к ослаблению стратификации в северной и центральной частях моря.

Ключевые слова: Японское море, численное моделирование, циклонический круговорот, межгодовая и декадная изменчивость циркуляции, ЭОФ анализ, SVD анализ.

Climate variability of the circulation in the Sea of Japan in the second part of the XX century. D.V. STEPANOV (V.I. Il'ichev Pacific Oceanological Institute, FEB RAS, Vladivostok).

Based on the numerical simulations, the climate variability of the Sea of Japan circulation from 1948 to 2009 was investigated. We confirm that one of the specific features of the circulation is the cyclonic gyre in the northern part of the Sea of Japan. The cyclonic gyre is intensified in spring and to be decayed in autumn. Analysis of the vortex component shows strong interannual and decadal variability in the intermediate and abyssal layers of the Sea ofJapan.

Intensification and decay of the cyclonic gyre have intervals of 3, 4 and 5 years. We found that on the decade scales from 1948 to 1965 the cyclonic gyre had been significantly intensified after which to the beginning of2000 this gear had been significantly decayed. We found that on the interannual scales the cyclonic gyre vortex variability is induced by the positive wind stress curl and on the decade scales it is induced by the positive wind stress curl and its meridional gradient over the Sea of Japan. The density stratification decay above the Central basin is induced by the sensible heat flux over the nearshore branch of the Tsushima current from 1958 to 1983. The decay of the sensible heat flux meridional gradient within 12—15 years results to the density stratification decay in the northern and central parts of the Sea of Japan.

Key words: the Sea of Japan, numerical simulations, cyclonic gyre, interannual and decadal variability, empirical orthogonal function (EOF) analysis, singular value decomposition (SVD) analysis.

В последние десятилетия наблюдаются значительные изменения в климатической системе Земли, характеризующиеся квазилинейным ростом интегрального (от

СТЕПАНОВ Дмитрий Вадимович - кандидат физико-математических наук, старший научный сотрудник (Тихоокеанский океанологический институт им. В.И. Ильичёва ДВО РАН, Владивосток). E-mail: step-nov@poi.dvo.ru

Работа осуществлена при финансовой поддержке РФФИ (проект 14-05-00255). Часть работы выполнена по госбюджетной тематике ТОИ ДВО РАН «Нелинейные динамические процессы в океане и атмосфере» (№ 01201363045).

поверхности до 2000 м) запаса тепла океана [12, 13], а также его глобального уровня [19]. Эти изменения неоднородны как в пространственном, так и во временном отношении и, следовательно, по-разному могут влиять на климат прибрежных территорий.

Японское море представляет собой квазизамкнутый бассейн, связанный с Тихим океаном и окружающими морями через мелкие проливы. Его термодинамический режим оказывает значительное влияние на климатические условия прибрежных зон граничащих с этим морем государств (России, Японии, Северной и Южной Кореи).

В результате исследования полей температуры и запаса тепла в приповерхностном слое Японского моря выявлена их значительная изменчивость на межгодовых и декадных временных масштабах. По данным натурных наблюдений установлено, что наиболее сильная изменчивость температуры промежуточных вод отмечается в северной части моря и имеет периодичность 5, 8 и 13 лет [3], а ее декадная изменчивость обусловлена положением Сибирского антициклона [14]. С другой стороны, авторы работы [16] показали, что декадная изменчивость запаса тепла приповерхностного слоя Японского моря связана с изменчивостью водообмена через Корейский пролив. Водообмен через проливы и атмосферные воздействия служат причинами изменчивости термического режима вод Японского моря. В то же время процессы адвекции и диффузии ответственны за перераспределение тепла по всему бассейну.

Согласно работе [6], крупномасштабная структура приповерхностной циркуляции вод Японского моря характеризуется тремя ветвями Цусимского течения, а также циклоническим круговоротом, охватывающим северную часть моря и состоящим из Приморского и Северо-Корейского течений, а также зоны субполярного фронта. На основе данных натурных наблюдений [8] установлено, что с глубиной в поле скорости преобладает крупномасштабный циклонический круговорот с характерными скоростями течений на его периферии >10 см/с и интенсивной мезомасштабной динамикой [18, 20]. Столь интенсивная и неоднородная по пространству циркуляция с неизбежностью будет влиять на формирование как поля температуры в приповерхностном слое, так и термического режима вод всего морского бассейна.

Несмотря на то что крупномасштабная структура динамики вод Японского моря в основном изучена, особенности ее изменчивости на климатических масштабах (>30 лет) до конца не выяснены. Отсутствие продолжительных натурных наблюдений, покрывающих весь морской бассейн от поверхности до дна, выводит на передний план численное моделирование, которое позволит раскрыть особенности изменчивости циркуляции вод, вызванные внешним воздействием (атмосферный форсинг и водообмен через проливы) на климатических масштабах.

В настоящей работе проводится анализ циркуляции вод и ее изменчивости на промежуточных горизонтах и в глубоководных слоях Японского моря, обусловленных климатическими изменениями атмосферного воздействия в период с 1948 по 2009 г. Данные по полям скоростей течений в бассейне получены с помощью численного моделирования. Расчеты выполнены с помощью численной модели циркуляции океана INMOM [1], а данные по атмосферному воздействию задавались в соответствии с экспериментом CORE, фаза II [12].

Описание модельной конфигурации

Для расчета циркуляции вод Японского моря под действием изменяющихся атмосферных условий в период с 1948 по 2009 г. использовалась численная модель циркуляции вод океана INMOM, основные характеристики которой, а также возможности ее приложений представлены в монографии [1]. Главная ее особенность заключается в применении процедуры задания граничных условий, что позволяет исключить из рассмотрения влияние изменчивости Тихого океана и окружающих морей на циркуляцию вод Японского моря. Подробно эта процедура описана в работах [2, 5].

На рис. 1 представлена модельная топография дна, полученная по данным ETOPO2 [7] и отражающая основные топографические особенности бассейна Японского моря. Для учета влияния мезомасштабных вихрей на крупномасштабную структуру циркуляции вод использовалось высокое горизонтальное (1/12° по широте и долготе) разрешение. На границе воздух-вода задавались напряжение ветра, поток тепла, учитывающий вклады коротко-и длинноволновой радиации, потоки явного и скрытого тепла, поток соли с учетом вклада от испарения и осадков, а также речного стока. Температура и влажность воздуха, скорость ветра на высоте 10 м, давление на уровне моря и радиация (длинно- и коротковолновая) принимались в соответствии с данными эксперимента CORE, фаза II за период с 1948 по 2009 г. [12]. Описание параметризаций процессов подсеточ-ных масштабов подробно представлено в работах [2, 5]. Расчеты проводились с начальных полей температуры и солености для января, полученных по данным работы [10], и нулевых полей скорости.

Рис. 1. Топография дна Японского моря по данным ЕТОР02. ЦК - Центральная котловина, КЯ - котловина Ямато, ПЯ -поднятие Ямато, КУ - котловина Уллындо

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

Рассмотрим результаты численного моделирования циркуляции вод Японского моря на интервале с 1948 по 2009 г. Для начала оценим, насколько модельное поле скорости отражает особенности динамики вод, выявленные по данным натурных наблюдений.

Согласно данным модельного эксперимента (рис. 2), для приповерхностного слоя вод характерна скорость течений до 29 см/с. Зимой почти весь бассейн Японского моря охвачен циклоническим круговоротом с интенсивными течениями у западного берега о-ва Хонсю (береговая ветвь Цусимского течения), а также у берегов Приморья (Приморское течение). Весной крупномасштабный циклонический круговорот наблюдается лишь в северной части моря. Наряду с этим отчетливо проявляется вдольбереговая ветвь Цусимского течения, антициклонический вихрь над Цусимской котловиной, а напротив Сангарского пролива - интенсивный циклонический вихрь. Летом, как и весной, циклонический круговорот формируется лишь в северной части Японского моря, а осенью он интенсифицируется и расширяется, охватывая часть котловины Ямато. Кроме этого прослеживается значительная интенсификация антициклонического вихря у берегов Северной Кореи.

На промежуточных горизонтах и в глубоководных слоях Японского моря в течение всего климатического года доминирует циклонический круговорот (рис. 3). Его граница отчетливо проявляется во все сезоны и характеризуется интенсивными течениями со скоростями от 12 до 15 см/с. Интенсивные течения наблюдаются в южной части моря: зимой в районе котловины Ямато, осенью в районе Цусимской котловины. Внутри циклонического круговорота формируются как циклонические, так и антициклонические

Рис. 2. Поле скорости в приповерхностном слое Японского моря по данным модельного эксперимента: а - зима, б - весна, в - лето, г - осень

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

Сравнение результатов численного моделирования с данными натурных наблюдений [8] показало, что в северной части Японского моря доминирует циклонический круговорот, внутри которого в северо-западной части Центральной котловины наблюдается циклонический вихрь синоптического масштаба (200-250 км). Согласно оценкам, характерные масштабы скоростей на периферии циклонического круговорота по натурным данным и данным модельного эксперимента составляют 10 и 10-12 см/с соответственно.

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

Рис. 3. Поле скорости на промежуточных горизонтах Японского моря по данным модельного эксперимента: а - зима, б - весна, в - лето, г - осень

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

Одной из характерных особенностей циркуляции вод Японского моря является ее вихревой характер [18, 20], поэтому исследование изменчивости циклонического круговорота проводилось на основе анализа поля относительной завихренности, характеризующей вихревую составляющую циркуляции вод моря [17] в период с 1948 по 2009 г.

Анализ осуществлялся на основе разложений полей по эмпирическим ортогональным функциям (ЭОФ) [4]. Предварительно из среднемесячных аномалий относительной завихренности на промежуточных горизонтах и в глубоководных слоях Японского моря был удален линейный тренд, а затем с помощью процедуры скользящего среднего

с окном, равным частоте отсечки 73 мес., они были разделены на «межгодовые» и «декадные» составляющие. С целью исключения из рассмотрения межсезонной и вну-тригодовой изменчивости для межгодовых аномалий относительной завихренности применялась процедура скользящего среднего с окном 13 мес. В качестве пространственных структур анализировались гомогенные корреляционные карты (ГКК), которые представляют собой пространственные распределения коэффициента корреляции между среднемесячными аномалиями относительной завихренности и главной компонентой ведущей ЭОФ моды.

Рассмотрим пространственно-временную изменчивость межгодовых аномалий относительной завихренности. Согласно разложению по ЭОФ две первые межгодовые ЭОФ моды охватывают более 70 % всей дисперсии, причем на первую приходится более 60 %, а на вторую - около 15 % общей дисперсии. Из рис. 4а и б видно, что максимальные значения коэффициента корреляции наблюдаются над Центральной котловиной Японского моря. Положительные значения коэффициента корреляции (черный цвет) указывают на преобладание циклонической завихренности внутри циклонического круговорота,

Рис. 4. Пространственные распределения первой (а) и второй (б) мод ЭОФ-разложения межгодовых аномалий относительной завихренности. Главные компоненты: первой (в) и второй (г) ЭОФ мод

Рис. 5. Пространственные распределения первой (а) и второй (б) мод ЭОФ-разложения декадных аномалий относительной завихренности. Главные компоненты: первой (в) и второй (г) ЭОФ мод

обусловленной неустойчивостью течений вдоль континентального склона, а также циклоническим вихрем в северной части Центральной котловины. Временная изменчивость завихренности на межгодовых масштабах характеризуется периодичностями 3, 4 и 5 лет с преобладанием 4-летнего периода (рис. 4в, г). Синфазность колебаний и незначительное уменьшение их амплитуды с глубиной сигнализируют о слабом затухании межгодовой изменчивости относительной завихренности по глубине.

На декадных масштабах вклады двух первых ЭОФ мод в их общую дисперсию составляют 50 и 23 % соответственно. Наиболее интенсивная декадная изменчивость относительной завихренности наблюдается над Центральной котловиной, охваченной циклоническим круговоротом (рис. 5а, б). На интервалах с 1950 по 1968 и с 2000 по 2007 г. преобладает положительная завихренность, сопровождающаяся интенсификацией циклонического круговорота (рис. 5в). Период с 1968 по 2000 г. характеризуется преобладанием антициклонической завихренности и ослаблением циклонического круговорота (рис. 5в, г). Кроме того, периоды ослабления циклонического круговорота наблюдаются с 1955 по 1987 г., а с 1987 по 2007 г. - их частичная интенсификация.

Анализ связей между вихрем напряжения ветра и циклоническим круговоротом вод Японского моря

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

По результатам сингулярного разложения (ЗУБ-разложения) [1] кросскорреляционной матрицы аномалий относительной завихренности и вихря напряжения ветра проведен анализ связей между ними на межгодовых и декадных временных масштабах.

На межгодовых масштабах суммарный вклад двух первых межгодовых SVD мод в совместную дисперсию составил 74 %, а коэффициент корреляции между вариациями коэффициента разложения первой и второй SVD мод составил 0,86 и 0,8 соответственно. Сильные связи наблюдаются между положительными аномалиями завихренности над Центральной котловиной Японского моря, охваченной циклоническим круговоротом, и положительными аномалиями вихря напряжения ветра над всем морем. Кроме того, значительный вклад в изменчивость завихренности вносят отрицательные аномалии

Рис. 6. Временные реализации коэффициентов разложения первой (а), второй (б) межгодовых ЗУБ мод, а также первой (в), второй (г) декадных ЗУБ мод, нормированные на максимальное значение. Сплошная линия - относительная завихренность, штриховая - вихрь напряжения ветра

зонального градиента вихря напряжения ветра. На временных масштабах сильная связь между относительной завихренностью и вихрем напряжения ветра наблюдается на периодах 3, 4 и 5 лет, которая не ослабевает с глубиной (рис. 6а, б).

На декадных масштабах на первые две SVD моды приходится до 84 % совместной декадной дисперсии. С глубиной вклад первой декадной SVD моды в совместную дисперсию становится преобладающим. Тем не менее коэффициенты корреляции между коэффициентами разложения для первой и второй декадных SVD мод почти не изменяются с глубиной и составляют 0,9 и 0,8 соответственно. Наиболее сильные связи прослеживаются между положительными аномалиями относительной завихренности на западной периферии Центральной котловины Японского моря, в ее северной части и положительными аномалиями вихря напряжения ветра над центральной частью, а также положительными аномалиями зонального градиента вихря напряжения ветра. В вариациях коэффициентов разложения, соответствующих первой и второй декадным SVD модам, выявлен лаг в регистрации отдельных сигналов. Например, сигнал, связанный с минимальной вариацией коэффициента разложения, зарегистрированный в 1980 г. для относительной завихренности, для вихря напряжения ветра наблюдается на 1-2 года раньше (рис. 6в, г).

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

Анализ связей между потоком явного тепла и стратификацией вод в зимний период года

По данным натурных наблюдений, в бассейне Японского моря фиксируются заметные различия в фоновых условиях плотностной стратификации между летним и зимним сезонами [10]. Значительное понижение температуры воздуха в зимний сезон приводит к увеличению потока тепла в атмосферу, сопровождающемуся интенсивным выхолаживанием верхнего приповерхностного слоя. Это влечет за собой перестройку плотностной стратификации вод, сформировавшейся в осенний период. В свою очередь, изменения плотностной стратификации должны привести к изменениям вихревой составляющей циркуляции вод Японского моря [17].

Согласно некоторым исследованиям [15], основной вклад в суммарный поток тепла над Японским морем в зимний период вносят потоки явного и скрытого тепла. Предполагается, что изменение потока явного тепла над морем в зимний сезон является одной из причин изменения стратификации плотности на промежуточных горизонтах моря.

Анализ межгодовых аномалий потока явного тепла и частоты плавучести не выявил значимых связей между ними. Согласно SVD-разложению декадных аномалий потока явного тепла и частоты плавучести, суммарный вклад первых трех SVD мод составил более 80 % общей декадной дисперсии, с преобладающим вкладом первой SVD моды (>50 %). Получены высокие значения коэффициентов корреляции между коэффициентом разложения, соответствующие всем трем модам SVD-разложения (от 0,7 до 0,9) и сигнализирующие о сильной связи между аномалиями потока явного тепла и частоты плавучести на декадных временных масштабах.

Анализ пространственных структур ГКК показал, что наблюдается сильная связь между отрицательными аномалиями частоты плавучести в северной части Японского моря и положительными аномалиями потока явного тепла в районе западной границы о-ва Хонсю в зоне прибрежной ветви Цусимского течения. На временных масштабах эта связь характеризуется отрицательными вариациями с 1948 по 1983 г., а после 1983 г. происходит смена знака вариаций для обоих коэффициентов разложения (рис. 7а). Значимая связь прослеживается также между отрицательными аномалиями частоты плавучести в

I -'-1-1-'-'-1-1-1-1-

1948 1953 1958 1963 1968 1973 1978 1983 1988 1993 1998 2003

Годы

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

центральной части моря и аномалиями меридионального градиента потока явного тепла, которая характеризуется квазипериодическими изменениями на масштабах 12-15 лет (рис. 7б). Кроме этого значимый вклад в изменчивость отрицательных аномалий частоты плавучести вносят положительные аномалии зонального градиента потока явного тепла на масштабах около 30 лет. Пик в вариациях коэффициента разложения наблюдается в 1959 и 1993 гг. (рис. 7в).

Заключение

На основе данных численного моделирования раскрыты особенности пространственно-временной изменчивости вихревой составляющей циклонического круговорота вод Японского моря, вызванной атмосферным воздействием в период с 1948 по 2009 г.

Результаты численного моделирования показали, что климатическое поле скоростей течений отражает основные элементы традиционной схемы течений. Одним из них является циклонический круговорот вод, расположенный в северной части Японского моря. На промежуточных горизонтах и в глубоководных слоях этот круговорот доминирует в структуре течений. Наблюдается его интенсификация в осенне-зимний и ослабление в летний сезоны. Характерные масштабы скорости течений на промежуточных горизонтах достигают 10-12 см/с, что согласуется с данными натурных наблюдений.

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

котловиной, т.е. преимущественно внутри циклонического круговорота вод. Характерные временные масштабы интенсификации и ослабления циклонического круговорота связаны с периодами 3, 4 и 5 лет. Выявлена значительная интенсификация круговорота в интервале с 1948 по 1965 г., после которого вплоть до начала 2000 г. наблюдался период его значительного ослабления. После этого циклонический круговорот вновь интенсифицировался вплоть до 2009 г.

На межгодовых масштабах изменчивость вихревой составляющей циклонического круговорота обусловлена изменчивостью положительного вихря напряжения ветра над Японским морем. Интенсификация циклонического круговорота в 1965 г. связана с усилением положительного вихря напряжения ветра над центральной частью моря, а события, связанные с интенсификацией круговорота на масштабах 20-25 лет и значительным его ослаблением в 1980 г., обусловлены изменчивостью меридионального градиента вихря напряжения ветра с лагом около 2 лет.

Анализ данных численного эксперимента показал, что с 1948 по 2009 г. в бассейне Японского моря наблюдаются значительные вариации в стратификации плотности в зимний сезон на декадных масштабах. Ослабление стратификации в юго-западной части Центральной котловины с 1958 по 1983 г. обусловлено увеличением потока явного тепла в атмосферу в районе вдольбереговой ветви Цусимского течения. С другой стороны, снижение меридионального градиента потока явного тепла в северной части моря на масштабах 12-15 лет сопровождается ослаблением стратификации в северной и центральной частях моря. Рост зонального градиента потока явного тепла в масштабе около 30 лет приводит к ослаблению стратификации вод над северо-западной частью Центральной котловины. Выявленный временной масштаб совместной изменчивости потока явного тепла и частоты плавучести согласуется с декадным масштабом, связанным с интенсификацией циклонического круговорота над Центральной котловиной Японского моря в середине 60-х и 90-х годах прошлого столетия.

ЛИТЕРАТУРА

1. Дианский Н.А. Моделирование циркуляции океана и исследование его реакции на короткопериодные и долгопериодные атмосферные воздействия. М.: Физматлит, 2013. 272 с.

2. Дианский Н.А., Степанов Д.В., Гусев А.В., Новотрясов В.В. Роль ветрового и термического воздействий в формировании изменчивости циркуляции вод в Центральной котловине Японского моря с 1958 по 2006 г. // Изв. РАН. Физика атмосферы и океана. 2016. Т. 52, № 1 (в печати).

3. Лучин В.А., Новотрясов В.В., Степанов Д.В. Межгодовая и декадная изменчивость температуры промежуточных вод Японского моря во второй половине XX века // Вестн. ДВО РАН. 2010. № 6. С. 30-35.

4. Мещерская А.В., Руховец Л.В., Юдин М.И., Яковлева Н.И. Естественные составляющие метеорологических полей. Л.: Гидрометеоиздат, 1970. 199 с.

5. Степанов Д.В., Дианский Н.А., Новотрясов В.В. Численное моделирование циркуляции вод центральной части Японского моря и исследование ее долгопериодной изменчивости в период 1958-2006 гг. // Изв. РАН. Физика атмосферы и океана. 2014. Т. 50, № 1. С. 84-96.

6. Юрасов Г.И., Яричин В.Г. Течения Японского моря. Владивосток: ДВО РАН, 1991. 176 с.

7. 2-minute Gridded Global Relief Data (ETOPO2v2). U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Geophysical Data Center. 2006.

8. Choi Y.J., Yoon J.H. Structure and Seasonal Variability of the Deep Mean Circulation of the East Sea (Sea of Japan) // J. Oceanogr. 2010. Vol. 66. P. 349-361.

9. Chu P.C., Chen Y., Lu S. Temporal and spatial variabilities of Japan Sea surface temperature and atmospheric forcing // J. Oceanogr. 1998. Vol. 54. P. 273-284.

10. Conkright M.E. et al. World Ocean Atlas 2001, Objective Analysis, Data Statistics, and Figures (CD-ROM Documentation, National Oceanographic Data Center, Silver Spring. Maryland, 2002).

11. Hogan P.J., Inagawa M. Sensitivity of simulated circulation dynamics to the choice of surface wind forcing in the Japan/East Sea // Deep-Sea Res. 2005. Vol. 52. P. 1464-1489.

12. Large W.G., Yeager S. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies // NCAR Tech. Note NCAR/TN-460+STR. 2004.

13. Levitus S., Antonov J., Boyer T., Stephens C. Warming of the World Ocean // Science. 2000. Vol. 287. P. 2225-2229.

14. Minobe S., Sako A., Nakamura M. Interannual to interdecadal variability in the Japan Sea based on a new grid-ded upper water temperature dataset // J. Phys. Oceanogr. 2004. Vol. 34. P. 2382-2397.

15. Na J., Seo J., Lie H.-J. Annual and Seasonal Variations of the Sea Surface Heat Fluxes in the East Asian Marginal Seas // J. Oceanogr. 1999. Vol. 55. P. 257-270.

16. Na H., Kim K.-Y., Chang K.-I., Park J.J., Kim K., Minobe S. Decadal variability of the upper ocean heat content in the East/Japan Sea and its possible relationship to northwestern Pacific variability // J. Geophys. Res. 2012. Vol. 117. P. C02017.

17. Pedlosky J. Geophysical Fluid Dynamics. Second edition. N.Y.: Springer-Verlag, 1987. 710 p.

18. Prants S.V. et al. Lagrangian analysis of the vertical structure of eddies simulated in the Japan Basin of the Japan/East Sea // Ocean Modelling. 2015. Vol. 86. P. 128-140.

19. Russell G.L., Gornitz V., Miller J.R. Regional sea level changes projected by the NASA/GISS atmosphere-ocean model // Climate Dynamics. 2000. Vol. 16. P. 789-797.

20. Takematsu M., Ostrovskii A.G., Nagano Z. Observations of Eddies in the Japan Basin Interior // J. Oceanogr. 1999. Vol. 55. P. 237-246.