2018
Известия ТИНРО
Том 193
УДК 551.465.62(265.51)
И.Д. Ростов1, Е.В. Дмитриева1, А.А. Воронцов2*
1 Тихоокеанский океанологический институт им. В.И. Ильичева ДВО РАН, 690041, г. Владивосток, ул. Балтийская, 43;
2 Всероссийский научно-исследовательский институт гидрометеорологической информации — Мировой центр данных, 249035, г. Обнинск Калужской обл., ул. Королева, 6
ТЕНДЕНЦИИ КЛИМАТИЧЕСКИХ ИЗМЕНЕНИЙ ТЕРМИЧЕСКИХ УСЛОВИЙ В ПРИБРЕЖНЫХ АКВАТОРИЯХ ЗАПАДНОЙ ЧАСТИ БЕРИНГОВА МОРЯ И ПРИЛЕГАЮЩИХ РАЙОНАХ ЗА ПОСЛЕДНИЕ ДЕСЯТИЛЕТИЯ
По данным срочных наблюдений на сети гидрометеорологических станций Росгидромета за 1980-2016 гг. и реанализа (по данным NOAA) выполнен анализ межгодовой изменчивости поверхностной температуры воды в прибрежных акваториях и прилегающих районах западной части Берингова моря и Камчатки на современном этапе «потепления». С использованием аппарата эмпирических ортогональных функций, методов кластерного и корреляционного анализа проведен анализ пространственно-временной структуры колебаний температуры и районирование акваторий по особенностям современных климатических изменений. Исследованы возможные причинно-следственные связи этих изменений с климатическими индексами. Проведенные исследования позволили выявить, уточнить и дать количественную оценку современных тенденций и региональных особенностей межгодовой изменчивости термических условий в выделенных районах.
Ключевые слова: Берингово море, западная часть, прибрежные районы, термические условия, межгодовая изменчивость, тренды температуры, климатические индексы, корреляционные связи, разложение полей по ЭОФ, региональные особенности.
DOI: 10.26428/1606-9919-2018-193-167-182.
Rostov I.D., Dmitrieva E.V., Vorontsov A.A. Tendencies of climate changes for thermal conditions in the coastal waters of the western Bering Sea and adjacent areas in the last decades // Izv. TINRO. — 2018. — Vol. 193. — P. 167-182.
The sea surface temperature (SST) variability in recent conditions of climate warming is analyzed for the coastal waters in the western Bering Sea and adjacent areas on the data of regular observations on coastal hydrometeorological stations in 1980-2016 and the data of NOAA_OI_SST_V2 reanalysis, taking into account the climate indices variations. The data series are analyzed using the methods of Empirical Orthogonal Functions, cluster analysis, and correlation analysis. Mean climatic trends are estimated for the air temperature as +0.27 °C/decade
* Ростов Игорь Дмитриевич, кандидат географических наук, заведующий лабораторией, e-mail: [email protected]; Дмитриева Елена Витальевна, кандидат технических наук, старший научный сотрудник, e-mail: [email protected]; Воронцов Александр Анатольевич, кандидат физико-математических наук, заведующий лабораторией, e-mail: [email protected].
Rostov Igor D., Ph.D., head of laboratory, e-mail: [email protected]; Dmitrieva Elena V., Ph.D., senior researcher, e-mail: [email protected]; Vorontsov Alexander A. Ph.D., head of laboratory, e-mail: [email protected].
and for SST as +0.23 0C/decade; they explain 27 % and 35 % of the year-to-year variability, respectively. The warming patterns are not identical for all seasons and areas. The most rapid growth of the air temperature and SST is observed for warm season (on average +0.36 oC/decade and +0.49 oC/decade). Both year-to-year fluctuations of SST and its climatic trend increase from south to north, the SST trend has increased from 0.11 to 0.32 oC/decade. Statistically significant climatic trends are found also in the changes of zonal and meridional components of wind velocity and certain indices of atmosphere circulation.
Key words: Bering Sea, western Bering Sea, coastal waters, thermal conditions, year-to-year variability, SST trend, index of atmosphere circulation, correlation, EOF analysis, regional features.
Введение
Процессы «глобального потепления», широко обсуждаемого в последние годы, оказывают значительное влияние на региональный климат Дальнего Востока. Происходящие изменения климатического режима отчетливо проявляются в последние 4 десятилетия в виде устойчивых тенденций и статистически значимых трендов межгодовых колебаний температуры воды, воздуха и других гидрометеорологических характеристик (Глебова и др., 2009; Ростов и др., 2017; WMO Statement..., 2017*). Эти изменения определяются крупномасштабными планетарными процессами, приводят к перестройке атмосферной циркуляции, смене типов атмосферных процессов, положения и выраженности сезонных центров действия атмосферы (Глебова, 2003; Пономарев и др., 2007; Шатилина и др., 2016; и др.).
Цель настоящей работы — уточнение тенденций, характеристик и возможных причинно-следственных связей межгодовых изменений термических условий в при-брежно-морских районах западной части Берингова моря и прилегающих районах за последние 4 десятилетия (в период общего «потепления» климатических условий в дальневосточном бассейне).
Материалы и методы
Для анализа межгодовых изменений термических условий использовались проконтролированные данные 4 срочных измерений температуры воды и воздуха на сети из 12 ГМС Росгидромета (рис. 1, табл. 1) за период c 1980 по 2016 г. из архивов Госфонда (http://meteo.ru/). Номер станции на рис. 1 соответствует последним значениям 5-значного общероссийского кода данного пункта наблюдений. Выбранные для анализа ГМС охватывают систематическими наблюдениями большую часть прибрежной зоны западной части Берингова моря, юго-восточного побережья Камчатки и о. Беринга и позволяют получить сопоставимые оценки тенденций и особенностей межгодовой изменчивости этих климатических переменных. Отобранные ряды наблюдений на этих станциях содержали минимальное количество пропусков, которые были восстановлены с помощью статистических методов (Ростов и др., 2017).
По данным срочных наблюдений рассчитывались и корректировались величины средних месячных, сезонных и средних годовых значений, что обеспечивало получение сглаженных, относительно однородных рядов межгодовых изменений температуры воды (Tw) и приземной температуры воздуха (Т). С учетом продолжительности ледового периода и особенностей внутригодового хода** были выделены условно теплый (июнь-октябрь) и условно холодный (ноябрь-апрель) сезоны (периоды года).
Также в работе были использованы данные реанализа поверхностной температуры (Twr) на регулярной сетке, поля ветра и климатические (циркуляционные) индексы, позволяющие оценить влияние и вклад крупномасштабных процессов в межгодовую изменчивость термических характеристик исследуемого района: индекс тихоокеанского
* WMO Statement on the state of the global climate in 2016. Geneva: World Meteorological Organization, 2017. WMO-N 1189. 26 p.
** Гидрометеорология и гидрохимия морей. Т. 10: Берингово море, вып. 1: Гидрометеорологические условия. СПб.: Гидрометеоиздат, 1999. 300 с.
Рис. 1. Схема расположения станций наблюдений и выделенных районов: СЗ — северозападный; Ц — центральный; ЮЗ — юго-западный
Fig. 1. Scheme of location for the coastal hydrometeorological stations and areas in the western Bering Sea: СЗ — northern area; Ц — central area; Ю3 — southern area
Таблица 1
Характер и тенденции межгодовых изменений температуры воздуха на прибрежных ГМС
за период наблюдений 1980-2016 гг.
Table 1
Year-to-year changes and long-term trends of annual mean air temperature at coastal stations
for 1980-2016
№ ст. Название T a T м a а2 D T± a КЛТ КЛТм КЛТт КЛТх
10 Эгвекинот (1980-2003 гг.) -5,6 -4,1 -7,6 0,98 2 0,4 0,19 - 0,27 0,23
12 Анадырь (1980-2005 гг.) -7,1 -4,4 -9,3 1,23 10 1,2 0,46 - 0,72 0,54
18 Апука (Олюторка) -1,9 03 -4,1 0,91 16 1,3 0,36 0,68 7 0,54 0,17
27 Оссора (Карага) -1,6 04 -3,0 0,63 10 0,8 0,23 0,49 8 0,27 0,19
31 Мыс Озерной -0,2 13 -1,5 0,54 19 1,2 0,29 0,57 9 0,31 0,25
32 Мыс Африка 1,5 0,3 0,39 10 0,7 0,18 0,42 9 0,25 0,11
34 О. Беринга (Никольское) 2,8 41 1,0 0,32 28 1,0 0,28 0,70 8 0,44 0,13
37 Семлячики 2,5 33 1,2 0,25 28 0,9 0,24 0,55 8 0,37 0,15
41 Петропавловский маяк 1,6 25 0,3 0,30 43 1,2 0,33 0,62 8 0,42 0,28
43 Петропавловск-Камчатский 2,8 37 1,7 0,26 31 1,0 0,27 0,56 8 0,37 0,18
48 Бухта Водопадная 2,4 34 1,3 0,21 24 0,8 0,21 0,48 9 0,31 0,14
50 Мыс Лопатка - - - - 1,3 0,37 0,94 7 0,68 0,07
Все станции - - - 27 1,0 0,27 0,62 8 0,36 0,18
Примечания. Та — среднее значение температуры, "С; Хм — экстремальные величины среднегодовых изменений (над чертой — максимальное значение, под чертой — минимальное); с2 — дисперсия Т; Б — вклад тренда в суммарную дисперсию, %; Т± — тренд за 37 лет; КЛТ — коэффициент наклона линейного тренда, оС/10 лет; КЛТм — максимальная величина этого коэффициента (над чертой — значение, под чертой — месяц); КЛТт, КЛТх — коэффициент тренда соответственно для теплого и холодного сезонов. Жирным шрифтом выделены оценки значимого тренда, прочерк — отсутствие сопоставимых данных.
декадного колебания (PDO), индекс южного (SOI) и арктического колебаний (AO), региональный индекс летнего охотоморского максимума (OKH), западно-тихоокеанский (WP) и северотихоокеанский (NPI) индексы (сайты NOAA http://www.cpc.ncep.noaa. gov и https://www.esrl.noaa.gov/psd/cgi-bin/data/getpage.pl).
Тренды оценивались для временных рядов, включающих среднегодовые значения наблюдаемых параметров и их аномалий, определяемых как отклонение от средней величины за все годы 30-летнего периода 1984-2013 гг. В качестве критического значения коэффициента корреляции на 5 %-ном уровне значимости линейного тренда была принята величина 0,33 (коэффициент детерминации R2 > 0,11).
Исследование структуры межгодовых изменений температуры воды проводилось с использованием метода разложения множества рядов (полей) на главные компоненты, или эмпирические ортогональные функции (ЭОФ). Районирование и выделение акваторий, относительно обособленных по характеру климатических изменений, проводилось по трем главным компонентам ЭОФ методами кластерного анализа (Ростов и др., 2017).
Результаты и их обсуждение
Общие характеристики межгодовой изменчивости температуры воздуха и воды в прибрежной зоне западной части Берингова моря и юго-востока Камчатки
Общее представление о характере и тенденциях межгодовых изменений температуры воздуха в исследуемом районе дают количественные оценки, приведенные в табл. 1.
Среднегодовые значения температуры воздуха на станциях возрастают вдоль побережья с севера на юг от минус 7,1 оС в Анадырском заливе до плюс 2,8 оС в районе Петропавловска-Камчатского и о. Беринга, а размах и дисперсия межгодовых колебаний, напротив, увеличиваются в 2-3 раза в противоположном направлении. На всех станциях тренды межгодовых изменений среднегодовых значений температуры положительны. Величины коэффициентов наклона линейного тренда изменялись от 0,18 до 0,46 оС/10 лет, что говорит о повсеместном возрастании температуры воздуха за последние 37 лет на величину от 0,4 до 1,4 оС (в среднем на 1,0 оС). Величина вклада тренда в суммарную дисперсию Ta у южной части побережья (24-43 %) больше, чем у северной (2-19 %). Обобщенно по всем ГМС оценки тренда Ta в теплый сезон в 1,3 раза превосходят среднегодовые, а наибольшие величины коэффициентов наклона линейного тренда (до 0,9 оС/10 лет) в сопоставимых по продолжительности рядах наблюдений наблюдались в июле-августе (табл. 1). В холодный сезон значимый положительный тренд отмечался только на одной ГМС — Петропавловский маяк.
На рис. 2 приведены обобщенные кривые временного хода среднегодовых величин аномалий и накопленных аномалий температуры воздуха и воды по данным всех прибрежных ГМС. Колебания аномалий температуры отражают чередование теплых и холодных лет в рядах наблюдений с интервалом 2-7 лет, а интегральные кривые накопленных аномалий — низкочастотные составляющие и отдельные периоды в межгодовом ходе температуры.
В начальный период наблюдений преобладали отрицательные аномалии температуры воздуха. В середине периода, до 2004 г., наблюдалось резкое возрастание амплитуды колебаний Ta с абсолютным максимумом в 1996 г., обусловленным аномальным потеплением в холодный период года. В последующий период происходил устойчивый рост суммы положительных аномалий температуры воздуха (рис. 2, а, б). Обобщенные по данным ГМС экстремально высокие среднегодовые величины аномалий температуры (AT > 1,5 оС) отмечались в 1995, 1996, 2003 и 2014 гг., а наименьшие (AT < 1,0 оС) — в а1985-1988, 1993 и 1998-2000 гг. Совокупно по всем ГМС величина статистически значимого линейного тренда изменений температуры воздуха составила 0,27 оС/10 лет (в теплый период 0,36 оС/10 лет), а величина вклада тренда в суммарную дисперсию D = 27 % (в теплый период 38 %). В прибрежной зоне Охотского
в г
Рис. 2. Межгодовая изменчивость аномалий и накопленных аномалий температуры воздуха (а, б) (среднегодовые) и температуры воды (теплый период) (в, г) на ГМС. Показан диапазон колебаний и линия тренда
Fig. 2. Year-to-year variability of anomalies and accumulated anomalies at coastal stations: а, б — annual mean air temperature; в, г — SST in the warm season. Range of fluctuations and linear climatic trends are shown
и северной части Японского морей эти величины (среднегодовые) несколько выше: соответственно 0,34 оС/10 лет, D = 43 % и 0,32 оС/10, D = 34 % (Ростов и др., 2017). В теплый период максимум 1996 г. в ходе аномалий температуры воздуха выражен слабо, а в ходе аномалий температуры воды он совсем отсутствует (рис. 2, в, г). Максимальная величина отрицательных аномалий (ATw = -2,4 оС) наблюдалась в 2001 г.
Количественные оценки межгодовых изменений температуры воды в исследуемом районе в теплый период года по данным 9 ГМС представлены в табл. 2.
Как следует из данных табл. 2, на большинстве станций отмечаются положительные тренды величиной от 0,15 до 0,40 °C/10 лет. На двух станциях, одна из которых расположена в мелководной бухте Карагинского залива (Оссора), а другая — у его южной границы (мыс Африка), обнаружены статистически незначимые отрицательные тренды. Максимальные положительные тренды наблюдались в августе и сентябре, а в холодный период значимые тренды отсутствовали. В среднегодовом ходе по имеющимся сопоставимым данным отмечались положительные тренды от 0,16 °C/10 лет (Петропавловск-Камчатский) до 0,25 °C/10 лет (о. Беринга).
Межгодовые флюктуации температуры воды и воздуха на большинстве ГМС характеризуются сильными корреляционными связями (r = 0,6-0,8). На некоторых участках побережья из-за нарушения естественного процесса прогрева и охлаждения водной поверхности они слабее (табл. 2). При осреднении данных по всем ГМС оценка статистически значимого КЛТ изменений T в теплый сезон составила 0,16 оС/10 лет, а
w ' '
величина вклада тренда в суммарную дисперсию — 10 %. В прибрежной зоне Охотского моря эти величины значительно выше: 0,34 оС/10 лет, D = 50 % (Ростов и др., 2017). В отличие от температуры воздуха данные наблюдений над температурой воды на отдельных ГМС отражают в основном локальные особенности термического режима отдельных участков прибрежных акваторий, включая мелководные бухты, которые зачастую изолированы от прилегающих районов моря.
Таблица 2
Характер и тенденции межгодовых изменений температуры воды в теплый период на прибрежных ГМС за период наблюдений 1980-2016 гг.
Table 2
Year-to-year changes and long-term trends of SST in warm season at coastal stations
in 1980-2016
№ ст. Название а2 D T ± w КЛТт КЛТм г
10 Эгвекинот 0,89 6 0,7 0,20 0,34 9 0,3
12 Анадырь (1980-2000 гг.) 0,57 - 0,2 0,15 0,48 9 0,8
18 Апука (Олюторка) 0,54 7 0,7 0,18 0,31 8 0,7
27 Оссора (Карага) 0,75 0 -0,1 -0,02 0,19 9 0,5
32 Мыс Африка 0,46 -1 -0,2 -0,05 0,31 8 0,4
34 О. Беринга (Никольское) 0,71 25 1,4 0,40 0,67 8 0,8
37 Семлячики 0,80 10 1,0 0,27 0,52 8 0,6
41 Петропавловский маяк 0,55 5 0,6 0,17 0,36 8 0,4
43 Петропавловск-Камчатский - - 1,0 0,28 0,68 8 0,7
Примечания. с2 — дисперсия Тш; Б — вклад тренда в суммарную дисперсию, %; Т^ — тренд за 37 лет; КЛТт — коэффициент наклона линейного тренда в теплый сезон, оС/10 лет; КЛТм — максимальная величина этого коэффициента (над чертой — значение, под чертой — месяц); г — коэффициент корреляции температуры воды и воздуха на станции. Жирным шрифтом выделены оценки значимого тренда, прочерк — отсутствие сопоставимых данных.
Тенденции межгодовой изменчивости температуры воды в акваториях, прилегающих к прибрежной зоне
По данным реанализа полей среднегодовых аномалий температуры воды ДТ^ в западной части Берингова моря и сопредельной области северо-западной части Тихого океана (СЗТО) методами кластерного анализа выделены 7 районов, три из которых — северо-западный (СЗ), центральный (Ц) и юго-западный (ЮЗ) — расположены в акваториях, прилегающих к прибрежной зоне (см. рис. 1). В северо-западный район входят Анадырский залив, наваринский и корякский шельф, западная часть Алеутской котловины; в центральный — шельф Карагинского и Олюторского заливов и Командорская котловина, расположенные в Беринговом море; в юго-западный, расположенный в СЗТО, — Камчатская впадина и Командорские острова.
На рис. 3 сопоставлены кривые временного хода среднегодовых аномалий и накопленных аномалий температуры воды в центрах кластеров выделенных районов. В каждом районе наблюдается постепенное увеличение температуры по сравнению с нормой базового периода. Корреляционная связь межгодовых колебаний температуры в удаленных районах (СЗ и ЮЗ) статистически значима (г = 0,40), но для соседних (СЗ/Ц и ЮЗ/Ц) она значительно теснее (г = 0,80 и 0,65). Величина коэффициента корреляции между изменениями Тот в теплый и холодный периоды значима только для районов СЗ и Ц.
В табл. 3 представлены количественные оценки тенденций и амплитудных характеристик вариаций температуры воды в центрах кластеров исследуемых районов. Как следует из данных, в направлении с юга на север наблюдается возрастание дисперсии и трехкратное увеличение скорости роста среднегодовой температуры воды:
Д е
Рис. 3. Межгодовая изменчивость аномалий AT и накопленных аномалий YAT темпе-
wr wr
ратуры воды по районам СЗ (а, б), Ц (в, г) и ЮВ (д, е): 1 — среднегодовые значения; 2, 3 — значения за теплый и холодный сезоны. Показан линейный тренд
Fig. 3. Year-to-year variability of SST anomalies and accumulated anomalies, by areas: а, б — northern part of the western Bering Sea; в, г — central part of the western Bering Sea; д, е — southern part of the western Bering Sea; 1 — annual means; 2 — means for warm season; 3 — means for cold season. Linear trend is shown
Таблица 3
Характер и тенденции межгодовых изменений температуры воды в выделенных районах
за период наблюдений 1982-2016 гг.
Table 3
Year-to-year changes and long-term trends of SST in certain areas for the 1982-2016
Район а2 D T ± КЛТ КЛТм КЛТт КЛТх
СЗ 0,29 38 1,1 0,32 0,78 9 0,64 0,08
Ц 0,25 28 0,9 0,25 0,76 8 0,52 0,06
ЮЗ 0,15 8 0,4 0,11 0,56 8 0,31 0,01
Вся акватория 0,17 35 0,8 0,23 0,70 8 0,49 0,04
Примечания. с2 — дисперсия Т; D — вклад тренда в суммарную дисперсию, %; Т^ — тренд за 35 лет; КЛТ — коэффициент наклона линейного тренда, оС/10 лет; КЛТм — максимальная величина этого коэффициента (над чертой — значение, под чертой — месяц); КЛТт, КЛТх — коэффициент тренда для теплого и холодного сезонов. Жирным шрифтом выделены оценки значимого тренда.
от 0,11 оС/10 лет (ЮЗ) до 0,32 оС/10 лет (СЗ). При этом величина вклада трендовой составляющей в суммарную дисперсию возрастала в том же направлении от 8 до 38 %. За исследуемый период среднегодовая температура воды в этих районах возросла от 0,4 до 1,1 оС.
В теплый период года во всех районах коэффициенты линейного тренда в 2-3 раза больше по сравнению со среднегодовыми и достигают максимальных положительных значений (до 0,6-0,8 оС/10 лет) в августе-сентябре. В холодный сезон тренды малы и статистически незначимы. В среднем по всему району исследований температура воды увеличивалась на величину 0,23 оС/10 лет, а в теплый период — на 0,49 оС/10 лет. Можно предположить, что быстрое потепление на севере исследуемой области во многом обусловлено поступлением в северо-западную часть Берингова моря более теплых, чем в прежние годы, воздушных масс из арктических районов. Предварительный анализ данных наблюдений на ГМС в восточном секторе Арктики показал, что величина положительных трендов температуры воздуха на станциях, расположенных на побережье Восточно-Сибирского и Чукотского морей, примерно в 2-3 раза превосходит соответствующие оценки, приведенные в табл. 1. Другим фактором потепления может быть увеличение переноса более теплых тихоокеанских вод Центрально-Беринговоморским (Склоновым) течением при усилении отдельных звеньев циклонического круговорота Берингова моря (Андреев, Жабин, 2015).
Анализ интегральных кривых аномалий £AT накопленных за исследуемый период (рис. 3, б, г, е), позволяет выделить и сопоставить низкочастотные составляющие межгодового хода температуры на различных временных интервалах в выделенных районах. Интервалы спада этих кривых отражают процесс выхолаживания вод, стабилизации — варьирования около средних значений, а интервалы роста — тренд накопления тепла в результате совокупного воздействия разнообразных адвективных и климатических факторов. Хорошо заметны различия в интенсивности (размах колебаний) и характере протекания (синхронности колебаний) этого процесса в выделенных районах. Процесс накопления отрицательных аномалий завершается в 2002-2003 гг., после чего во всех районах начинается фаза потепления.
Исходя из величины соотношения среднегодовых аномалий ATwr и их стандартных отклонений (Карпова, Шатилина, 2000) были выделены годы, соответствующие 5 режимам теплового состояния вод: аномально холодные, холодные, средние (нормальные), теплые и аномально теплые. Результаты расчетов приведены в табл. 4.
Таблица 4
Годы с наблюдениями аномальных режимов теплового состояния вод в выделенных районах
Table 4
Years with abnormal regime of thermal conditions in certain areas of the western Bering Sea
Режим вод Северо-западный район Центральный район Юго-западный район
Аномально холодный - 2001 1992,2001
Холодный 1982, 1985, 1987, 1988, 1992,1998-2000, 2012 1982, 1985, 1987, 1988, 1992, 1995, 1998-2000, 2002, 2012 1982, 1985, 1987, 1988, 1999, 2002,2012
Теплый 1996, 1997, 2005-2008, 2013 1991, 1996, 1997, 2003-2006, 2010, 2011, 2015, 2016 1997, 1998, 2005, 2006, 2009, 2011-2013
Аномально теплый 2003, 2004, 2014-2016 2008, 2014 1986,1996,2014
В многолетнем ходе температуры воды в последнее десятилетие рассматриваемого периода на фоне чередования «холодных» и «теплых» лет хорошо выражена общая тенденция потепления вод.
Анализ структуры межгодовых изменений температуры воды
Отмеченные особенности пространственно-временной изменчивости поля температуры воды отражают результат влияния и взаимодействия большого числа разно-
масштабных процессов, происходящих в море и атмосфере. Для выявления наиболее важной информации о структуре этого поля по данным реанализа были вычислены ЭОФ разложения среднегодовых и сезонных колебаний аномалий ATwr в западной части Берингова моря и прилегающих районов СЗТО. Выявлено, что первые три моды ЭОФ аккумулируют 72-76 % дисперсии межгодовых колебаний температуры воды и дают наиболее полное представление об особенностях структуры межгодовой изменчивости термических полей исследуемого района (табл. 5).
Таблица 5
Вклад первых трех мод ЭОФ разложения полей аномалий температуры в дисперсию межгодовых изменений, %
Table 5
Contribution of the first three EOF modes in dispersion of SST year-to-year variability, %
ЭОФ C-1 C-2 C-3 Xе
Среднегодовые 51 19 5 75
Теплый сезон 41 28 7 76
Холодный сезон 52 12 9 73
Примечание. Жирным шрифтом выделены оценки значимого тренда.
Первая мода (С-1) описывает 51 % общей дисперсии вариаций среднегодовых значений ДТ а в теплый и холодный сезоны — 41 и 52 %. Статистически значимые положительные тренды в ходе временных рядов коэффициентов разложения выделялись только для второй моды (в теплый сезон и для среднегодовых значений ДТ^) и третьей моды (холодный сезон).
В целом для временных колебаний температуры характерна слабая пространственная связанность, более тесная в зональном направлении. Поле пространственного распределения коэффициентов первых трех мод ЭОФ представлено как положительными, так и отрицательными значениями, что свидетельствует о несинхронности колебаний температуры и различной интенсивности процессов теплоотдачи, аккумуляции и адвекции тепла в различных участках акватории моря (рис. 4). Они обусловлены различиями климатических условий, гидрологического режима выделенных районов и определяются множеством факторов, одним из которых является несинхронность и даже противофазность гидрометеорологических процессов, протекающих в восточной и западной частях Берингова моря*.
В теплый период года в распределении коэффициентов первой моды (41 % общей дисперсии) выделяются две области противофазных колебаний — на северо-востоке и юго-западе рассматриваемой акватории. Первая из них полностью включает район СЗ с входящим в него участком беринговоморского шельфа, а вторая — районы Ц и ЮЗ, расположенные в глубоководной части акватории. Граница между этими областями проходит южнее беринговоморского шельфа примерно вдоль границы максимального распространения дрейфующего льда. Можно предположить, что положение этих областей в основном связано с различием скорости процессов теплоотдачи (выхолаживания) и аккумуляции тепла в шельфовом и глубоководном участках акватории. Это обусловливает существование очагов аномалий температуры воды и воздуха различного знака (Устинова и др., 2008; Хен и др., 2008) и отражает региональные особенности термического и динамического воздействия муссонной атмосферной циркуляции и различия гидрологических условий северо-восточной и юго-западной частей моря.
В пространственном распределении коэффициентов второй моды в этот же сезон (28 % дисперсии) прослеживается синхронность вариаций температуры в пределах всего побережья и большей части прилегающих районов. За пределами этой области, в открытой части акватории моря, одновременные изменения температуры происходят в противофазе. Эта мода может быть связана с влиянием локальных синоптических процессов и формированием зон преобладающих траекторий и максимальной по-
в. д. 160° 165° 170° 175° 180° 185° 1980 1990 2000 2010 2020
Д е
Рис. 4. Пространственное распределение первой (а), второй (в) и третьей (д) мод эмпирических ортогональных функций и временной ход (б, г, е) соответствующих коэффициентов разложения по эмпирическим ортогональным функциям аномалий поля температуры воды в теплый период года
Fig. 4. Spatial distribution of the first (а), second (в) and third (д) EOF modes for SST anomalies in warm season and temporal variability of their temporal coefficients (б, г, е)
вторяемости циклонов, приуроченных к крупномасштабным атмосферным фронтам (Гидрометеорология..., 1999*; Глебова, 2003; Устинова и др., 2008). Третья мода (7 % дисперсии) характеризует синфазные колебания в поле температуры небольшой амплитуды в зональном и меридиональном направлениях по отношению к центральной части Алеутской котловины и может отражать региональные особенности динамики вод моря и перераспределения трансформирующихся тихоокеанских вод по акватории (Hughes et al., 1974).
Общие особенности структуры поля температуры, описываемые первой и второй модами в теплый период года, характерны и для холодного сезона. В это время происходит смещение к северу границы, разделяющей области положительных и отрицательных значений коэффициентов ЭОФ и уменьшение дисперсии колебаний температуры воды. Поле же третьей моды в холодный период (9 % общей дисперсии
колебаний) существенно перестраивается. В это время года область положительных значений коэффициентов ЭОФ распространяется на всю акваторию Алеутской котловины и южнее, а отрицательных, противофазных — охватывает все прибрежные районы Чукотки и восточной Камчатки. Межгодовые вариации температуры воды и воздуха в выделенных районах значимо коррелируют с колебаниями временных коэффициентов главным образом первой и второй мод ЭОФ (табл. 6).
Таблица 6
Коэффициенты корреляции между рядами межгодовых колебаний аномалий температуры воздуха АТа, воды ATwr в выделенных районах и временных коэффициентов первых мод эмпирических ортогональных функций полей аномалий ATwr в среднем за год,
в теплый и холодный сезоны
Table 6
Coefficients of correlation between year-to-year fluctuations of air temperature and SST in the northern, central, and southern parts of the western Bering Sea and temporary coefficients of the first EOF modes of SST (annual mean and average for warm and cold seasons)
Показатель Северо-западный район Цент ральный район Юго-западный район
Ср. Тепл. Хол. Ср. Тепл. Хол. Ср. Тепл. Хол.
АТ wr
С-1 0,4 0,5 0,5 -0,1 0,0 0,0 -0,5 -0,3 -0,5
С-2 0,7 0,6 0,5 0,7 0,8 0,6 0,2 0,5 0,0
С-3 0,1 0,0 0,3 0,0 -0,2 0,1 -0,1 0,0 -0,1
АТ а
С-1 0,2 0,4 0,0 0,0 0,0 -0,1 -0,2 -0,2 -0,3
С-2 0,6 0,6 0,2 0,6 0,5 0,3 0,6 0,5 0,3
С-3 0,0 -0,1 0,0 0,0 -0,2 -0,1 0,0 -0,2 0,0
Примечание. Значимые величины выделены жирным шрифтом.
В этих связях во временном ходе аномалий температуры воды и коэффициентов первой моды во все сезоны года отчетливо проявляется противофазность колебаний в районе СЗ по отношению к району ЮЗ (рис. 4, а). Корреляционные связи между изменениями аномалий АТ^, и коэффициентов второй моды характеризуются только положительным знаком. В холодный период они ослабевают, а на юге района полностью отсутствуют. Временная изменчивость коэффициентов третьей моды слабо коррелирована с колебаниями АТот только в районе СЗ в холодный период года. Как отмечалось выше, статистически значимые положительные тренды температуры воды (среднегодовые значения и в теплый сезон) присутствовали во всех выделенных районах. В вариациях коэффициентов ЭОФ такие тренды выделялись только в изменениях второй (среднегодовые и в теплый сезон) и третьей мод (холодный сезон).
В отличие от температуры воды корреляционные связи между временными изменениями коэффициентов мод С1-С3 и колебаниями температуры воздуха на прибрежных ГМС, расположенных в выделенных районах, выражены главным образом только для второй моды в теплый период года, что проявляется и в среднегодовых значениях коэффициентов корреляции (табл. 6).
Корреляционные связи между характеристиками термических условий и атмосферными процессами
Отмеченные особенности пространственно-временной изменчивости поля температуры — результат комбинации различных факторов и сложного взаимодействия процессов глобального, регионального и локального масштабов. В условиях глобального потепления причинами формирования областей положительных аномалий приземной температуры воздуха и воды в дальневосточном регионе являются изменения крупномасштабной атмосферной циркуляции, расположения и выраженности центров действия атмосферы, а также интенсивности и соотношения зональных и меридиональных
переносов в атмосфере. Эти изменения характеризуются рядом признаков и отражаются в отдельных циркуляционных показателях, в частности в колебаниях климатических индексов (Смирнов, Воробьев, 2002; Хен и др., 2008; Котляков и др., 2015; Ростов и др., 2017; и др.). Был проведен взаимный корреляционный анализ временных рядов температуры воды и воздуха, коэффициентов мод ЭОФ разложения полей аномалий ДТжг и ветра, а также индексов атмосферной циркуляции, характеризующих динамику климатической системы северо-западной части Тихого океана.
Произошедшие в результате климатического сдвига нарушения барического поля и сезонного режима региональной муссонной циркуляции выражены в изменениях ветрового режима в холодный и теплый сезоны года (рис. 5, а-г). Поле приводного ветра существенно неоднородно и характеризуется значительной завихренностью. Приведенные на рис. 5 (д, е) кривые указывают на тенденции некоторого ослабления зимних муссонных ветров северных румбов и усиления летнего муссона при статистически значимых трендах возрастания скорости ветра западных и южных румбов. Величина тренда скорости ветра западных направлений за теплый период года составила 0,6 м/с, а его вклад в суммарную дисперсию — 15 %. В ряде работ причина потепления в Северной Евразии связывается с усилением межширотного обмена при ослаблении зонального переноса в средней тропосфере (Попова, Шмакин, 2010; Попова, 2014; и др.).
Рис. 5. Поля приземного ветра в холодный (а, в) и теплый (б, г) сезоны года в начале (1983 г.) и второй половине (2011 и 2016 гг.) периода наблюдений и сглаженные 3-летним скользящим средним меридиональная V (1) и зональная U (2) компоненты скорости ветра в холодный (д) и теплый (е) сезоны
Fig. 5. Wind fields at the sea surface in cold (а, в) and warm (б, г) seasons in 1983 and 2011, 2016 and meridional (1) and zonal (2) components of wind velocity in cold (д) and warm seasons (е) smoothed by 3-years averaging
Ветровое воздействие оказывает влияние как на процессы переноса тепла в системе течений всего Субарктического круговорота, так и на характер теплообмена между морем и атмосферой, что приводит к формированию очагов аномалий температуры воды и воздуха в различных районах моря. Межгодовая изменчивость полей атмосферного давления и ветра связана с состоянием основных сезонных центров действия атмосферы, их взаимодействиями и региональными особенностями проявления, а также с характером циклонической деятельности, приуроченной к крупномасштабным атмосферным фронтам*. При этом циклические колебания скорости ветра (рис. 5, д, е) примерно соответствуют периоду межгодовых изменений повторяемости холодных и теплых синоптических типов (Глебова, 2005). Из-за существенной неоднородности структуры поля ветра коэффициенты корреляции между колебаниями температуры воздуха и воды в отдельных районах с меридиальной и зональной компонентами V и U, осредненными по всему морю, относительно невелики и значимы только для отдельных сезонов (табл. 7).
Таблица 7
Коэффициенты корреляции между рядами временных колебаний меридиональной (V) и зональной (U) компонент поля ветра и межгодовых аномалий температуры воды и воздуха в выделенных районах
Table 7
Coefficients of correlation between temporary changes of meridional (V) and zonal (U) components of wind velocity and year-to-year anomalies of SST and air temperature, by areas
Показатель Северо-западный район Цент ральный район Юго-западный район
Ср. Тепл. Хол. Ср. Тепл. Хол. Ср. Тепл. Хол.
Т а
V 0,2 0,0 0,3 0,3 0,1 0,4 0,2 0,1 0,4
U -0,2 0,1 -0,6 -0,2 0,3 -0,6 0,1 0,5 -0,3
Т w
V 0,1 0,1 0,0 0,0 -0,1 0,0 0,1 0,3 0,1
U -0,1 0,2 -0,3 -0,1 0,3 -0,5 -0,2 0,3 -0,3
Примечание. Значимые величины выделены жирным шрифтом.
В поле Та корреляционные связи проявляются во всех районах с колебаниями как меридиональной, так и зональной компонент скорости ветра. В период зимнего муссона ослабление ветров северной четверти сопровождается ростом температуры воздуха. Колебания температуры воды коррелированы с изменениями зональной компоненты скорости в двух районах — СЗ и Ц — только в холодный период. При усилении западного переноса в холодный период года, когда Аляскинское течение ослабевает и поступление более теплых океанических вод сокращается, происходит понижение температуры воды в районах СЗ и Ц Берингова моря. В теплый период года усиление ветров западной четверти (см. рис. 4, е) сопровождается ростом температуры воздуха в районах Ц и ЮЗ (табл. 7).
Для уточнения характеристик пространственно-временной структуры поля ветра были вычислены ЭОФ вариаций компонент V и и скорости. Первые три моды характеризуют пространственное положение основных энергоактивных зон в пределах акватории, где вклад изменчивости поля ветра в суммарную дисперсию достигает 71-78 %. Отмеченные выше тенденции межгодовых изменений поля ветра (рис. 5, д, е) проявляются в наличии статистически значимых положительных трендов во временном ходе коэффициентов ЭОФ разложения в теплый период года. Анализ корреляционных связей межгодовых колебаний полей ветра и температуры воды показал, что они статистически значимы только для первой (все районы), а также первой и второй (районы Ц и ЮЗ) мод ЭОФ разложения V и и компонент скорости
ветра. Интерпретация этих связей сложна и неоднозначна, она требует рассмотрения конкретных синоптических ситуаций и самостоятельного анализа.
Рассмотрим результаты корреляционного анализа связей межгодовых колебаний температуры воды и воздуха, характеристик поля ветра с вариациями региональных климатических (циркуляционных) индексов. В рассматриваемый период в ходе межгодовых изменений климатических индексов присутствовали статистически значимые положительные и отрицательные трендовые составляющие (табл. 8), которые указывают на произошедшую в конце 90-х — начале 2000-х гг. смену режима и ослабление крупномасштабной атмосферной циркуляции в северотихоокеанском регионе (Ростов и др., 2017).
Таблица 8
Коэффициенты корреляции межгодовых колебаний климатических индексов c вариациями компонент скорости ветра, аномалий температуры воздуха и воды в выделенных районах и коэффициентов временных рядов первых мод ЭОФ
поля аномалий Т
wr
Table 8
Coefficients of correlation between year-to-year changes of climatic indices, components of wind velocity, SST and air temperature anomalies and EOF temporal coefficients of SST anomalies. Significant values are shown by bold numerals
Показатель PDO NPI WP SOI AO OKH
Теплый сезон (июль)
КЛТ -0,5 - -0,3 0,2 -0,3 -0,8
V -0,2 - 0,6 0,2 0,3 -0,3
U -0,5 - 0,4 0,4 0,3 -0,5
Та (СЗ) -0,4 - -0,5 0,4 0,2 -0,3
Та (Ц) -0,3 - -0,4 0,4 -0,2 -0,4
Та (ЮЗ) -0,4 - -0,4 0,4 -0,2 -0,3
Тт (СЗ) 0,0 - -0,4 0,3 -0,3 -0,4
Т„г (Ц) -0,2 - -0,3 0,3 -0,1 -0,4
Т„г (ЮЗ) -0,2 - -0,4 0,2 -0,1 -0,3
С-1 0,4 -0,3 -0,2 -0,1 0,0
С-2 -0,4 -0,3 0,4 0,0 -0,6
С-3 -0,3 0,1 0,0 -0,2 -0,3
Холодный сезон (февраль)
КЛТ -0,3 0,8 0,5 0,3 0,1 -
V -0,4 0,7 0,3 0,4 0,4 -
U -0,3 0,6 0,6 0,1 0,3 -
Та (СЗ) 0,1 0,1 -0,5 0,1 0,0 -
Та (Ц) 0,0 0,2 -0,6 0,2 0,0 -
Та (ЮЗ) 0,0 0,3 -0,5 0,2 0,1 -
Т„г (СЗ) 0,6 -0,3 -0,2 -0,1 -0,2 -
Т„г (Ц) 0,5 -0,3 -0,2 -0,2 -0,1 -
Т„г (ЮЗ) 0,0 0,0 -0,2 0,0 -0,1 -
С-1 0,6 -0,4 0,1 -0,1 -0,1
С-2 0,0 0,0 -0,1 0,1 0,1
С-3 0,2 0,0 0,1 -0,3 -0,3
Примечания. КЛТ — коэффициент линейного тренда, отн. ед./10 лет. Прочерк — отсутствие оценок из-за сезонности центров действия атмосферы. СЗ, Ц, ЮЗ — обозначения выделенных районов. Жирным шрифтом выделены статистически значимые величины.
В теплый сезон, когда наблюдаются статистически значимые положительные тренды межгодовых вариаций температуры воды и воздуха, в ходе рассматриваемых климатических индексов произошли следующие изменения.
— Смена фазы колебаний PDO c положительной на преимущественно отрицательную, с чем обычно связывается начало процесса накопления положительных аномалий в северной и западной частях Тихого океана (Mantua et al., 1997; Данова, 2011).
—Переход значений индекса WP к отрицательной фазе, который можно связать с уменьшением градиента давления между субтропиками и субполярными широтами, блокированием западного переноса, ослаблением зональной атмосферной циркуляции в средней тропосфере на уровне АТ-500, усилением межширотного обмена и потеплением в Северо-Тихоокеанском секторе и на севере Евразии (Попова, Шмакин, 2010; Попова, 2014; Котляков и др., 2015).
—Уменьшение значений OKH, свидетельствующее об ослаблении охотского антициклона. При слабовыраженном охотском антициклоне наблюдается усиление циклонической активности над СЗТО, выход циклонов и вынос теплых воздушных масс на акваторию Охотского моря (Шатилина и др., 2016).
Эти процессы происходят синхронно с межгодовыми изменениями температуры воды и воздуха и являются следствием изменения давления в климатических центрах действия атмосферы, их смещения и тесных связей с крупномасштабными колебаниями в системе океан-атмосфера всей северной части Тихого океана (Пономарев и др., 2007; и др.). В течение холодного сезона наблюдается только значимый положительный тренд в ходе NPI, что указывает на ослабление влияния на циркуляционную деятельность алеутской депрессии. Между изменениями климатических индексов и вариациями климатических параметров, характеризующих состояние термического и барического полей, имеются статистически значимые прямые и обратные корреляционные связи (табл. 8). Как было показано выше, благодаря этим воздействиям в приземном слое и на поверхности моря формируется сложная пространственная многомодовая структура распределения амплитуд и дисперсии этих колебаний, а также статистически значимые тренды межгодовых изменений полей температуры воды, воздуха и скорости ветра. Во всех районах корреляционные связи лучше всего выражены в начале теплого сезона (июль) — периода наибольших трендов в межгодовых изменениях климатических параметров (см. табл. 1-3, 8). Так, в поле температуры воздуха в теплый сезон для большинства выделенных районов коэффициенты корреляции разного знака значимы для всех климатических индексов, кроме АО, а в холодный — только для WP. В поле температуры воды в теплый сезон для большинства районов коэффициенты корреляции разного знака значимы для WP, OKH и SOI, а в холодный — только для PDO. Эти корреляционные связи проявляются в ходе временных коэффициентов различных мод ЭОФ как в теплый, так и в холодный сезоны года.
Заключение
По данным наблюдений на ГМС и реанализа были определены общие тенденции и региональные особенности межгодовых изменений характеристик термического режима прибрежных акваторий западной части Берингова моря и прилегающих районов в 1980-2016 гг. на фоне глобального «потепления» климата.
Климатические изменения термического режима в исследуемом районе проявляются в увеличении температуры воздуха в среднем со скоростью 0,27 оС/10 лет, а воды — 0,23 оС/10 лет, что объясняет соответственно 27 и 35 % их межгодовой изменчивости. Для поля температуры воздуха эта величина больше глобальных оценок (0,17 оС/10 лет), но меньше значений, характерных для северной части Японского (0,32 оС/10 лет) и Охотского морей (0,34 оС/10 лет) за тот же период. Для поля температуры воды в районе исследований эта величина несколько выше соответствующих оценок для названных морей (0,19 и 0,17 оС/10 лет).
Характер потепления неодинаков в различные сезоны и в отдельных районах. Наиболее быстрый рост температуры воздуха и воды наблюдался в теплый сезон (в среднем 0,36 и 0,49 оС/10 лет). В направлении с юга на север в выделенных районах наблюдается возрастание дисперсии колебаний и трехкратное увеличение скорости роста среднегодовой температуры воды от 0,11 до 0,32 оС/10 лет. Величина вклада трендовой составляющей в суммарную дисперсию изменялась в том же направлении от 8 до 38 %, а температура воды за исследуемый период в этих районах возросла от 0,4 до 1,1 оС.
Исследование особенностей пространственной структуры поля межгодовых колебаний температуры воды и корреляционных связей между характеристиками термических условий и атмосферными процессами показало, что в результате атмосферных
воздействий на поверхности моря формируется сложная многомодовая структура распределения амплитуд и дисперсии колебаний температуры. Произошедшие в результате климатического сдвига нарушения барического поля и сезонного режима региональной муссонной циркуляции выражены в изменениях ветрового режима как в холодный, так и в теплый сезоны года. Во временном ходе межгодовых изменений зональной, меридиональной компонент скорости ветра и климатических (циркуляционных) индексов присутствуют статистически значимые трендовые составляющие. В поле температуры воздуха в теплый сезон для большинства выделенных районов коэффициенты корреляции разного знака значимы для всех индексов, кроме АO, а в холодный — только для WP. В поле температуры воды в теплый сезон для большинства районов коэффициенты корреляции разного знака значимы для WP, OKH и SOI, а в холодный — только для PDO.
Список литeрaтуры
Андреев А.Г., Ж;1бш i И.А. Влияние продолжения Аляскинского течения на динамику вод восточной части Oхотского моря // Вестн. flPO РАИ. — 2015. — № 2. — C. S7-92.
Глебов! С.Ю. Изменения атмосферного режима над Дальневосточным регионом в 2000-2004 гг. и предполагаемые тенденции развития на ближайшие годы // Изв. l^OTO. — 2005. — Т. i42. — C. 214-222.
ГЛeбoвa С.Ю. Типы атмосферных процессов над дальневосточными морями, межгодовая изменчивость их повторяемости и сопряженность // Изв. 'MOTO. — 2003. — Т. i34. — C. 209-257.
Глебов! С.Ю., Устинов! Е.И., Сорокин Ю.Д. Долгопериодные тенденции в ходе атмосферных процессов и термического режима дальневосточных морей за последний 30-летний период // Изв. TOOTO. — 2009. — Т. 159. — C. 2S5-29S.
Динови Т.Е. Взаимосвязь региональных климатических индексов и динамики морского льда Арктического бассейна // Региональные проблемы. — 2011. — Т. i4, № i. — C. 42-47.
Кирпови И.П., Шитмлмни Т.А. Долгопериодная изменчивость температуры воды и воздуха у юго-западного побережья Cахалина // Изв. ТИHРO. — 2000. — Т. 127. — C. 50-60.
Котляков В.М., Хромов! Т.Е., Носенко Г.А. и др. ^временные изменения ледников горных районов России : моногр. — М. : Тов-во научных изданий КМК, 2015. — 2SS с.
Пoнoмaрeв В.И., Kaплунeнкo Д.Д., Дмитров! Е.В. и др. Климатические изменения в северо-западной части Азиатско-Тихоокеанского региона // Дальневосточные моря России. Кн. 1 : Oкеанологические исследования. — М. : Hаука, 2007. — C. 17-4S.
Понови В.В. Летнее потепление на европейской территории России и экстремальная жара 2010 г. как проявление тенденций крупномасштабной атмосферной циркуляции в конце XX в. — начале XXI в. // Метеорол. и гидрол. — 2014. — № 3. — C. 37-49.
Попов! В.В., Шишкин А.Б. Региональная структура колебаний температуры приземного воздуха в северной Евразии во второй половине XX — начале XXI веков // Изв. РАЙ. Физика атмосферы и океана. — 2010. — Т. 46, № 2. — C. 161-175.
Ростов И.Д., Дмитриев! Е.В., Воронцов A.A. Тенденции климатических изменений термических условий прибрежных районов Oхотского моря за последние десятилетия // Изв. TOOTO. — 2017. — Т. 191. — C. 176-195.
Смирнов Н.П., Воробьев В.Н. Cеверо-тихоокеанское колебание и динамика климата в северной части Тихого океана : моногр. — CM. : РГГМУ, 2002. — 122 с.
Устинов! Е.И., Глебов! С.Ю., Сорокин Ю.Д. Гидрометеорологические условия дальневосточных морей и северо-западной части Тихого океана в 200S г. // Вопр. промысл. океанологии. — 200S. — Вып. 5, № 2. — C. 4S-67.
Хен Г.В., Бискж' Е.О., Сорокин Ю.Д. и др. Термические условия на поверхности Берингова и Oхотского морей в начале 21-го века на фоне полувековой изменчивости // Изв. TOOTO. — 200S. — Т. 153. — C. 254-263.
Шитилини Т.А., Ц||||111111ми1. ni Г.Ш., Paдчeнкoвa Т.В. Oценка тенденций изменчивости центров действия атмосферы над Азиатско-Тихоокеанским регионом в летние периоды 1950-1979 и 19S0-2012 гг. // Метеорол. и гидрол. — 2016. — № 1. — C. 17-2S.
Hughes F.W., Coachman L.K., Aagaard K. Circulation, transport and water exchange in the Bering Sea // Oceanography of the Bering Sea with emphasis on renewable resources. Int. Symp. for Bering Sea study. — Fairbanks, Alaska, 1974. — P. 59-9S.
Mantua N.J., Hare S.R., Zhang Y. et al. A Pacific interdecadal climate oscillation with impacts on salmon production // Bull. Am. Meteorol. Soc. — 1997. — Vol. 7S, № 6. — P. 1069-1079.
Поступила в редакцию 28.03.18 г.
Принята в печать 12.04.18 г.