Североатлантическое колебание и увлажнение европейской территории России

ПРИРОДНАЯ СРЕДА

УДК 551.468.6 ББК 26

В.Н. Малинин, С.М. Гордеева

СЕВЕРОАТЛАНТИЧЕСКОЕ КОЛЕБАНИЕ И УВЛАЖНЕНИЕ ЕВРОПЕЙСКОЙ ТЕРРИТОРИИ РОССИИ*

Рассматриваются межгодовые закономерности североатлантического колебания -важнейшего индикатора взаимодействия океана и атмосферы в Северной Атлантике. Выделены 4 фазы развития САК с 1875 года: по две усиления и ослабления. Показано, что при положительной фазе САК происходит значительное усиление зональной циркуляции на уровне 500 гПа в зоне 50-70° с.ш., что приводит к увеличению переноса водяного пара в системе средней циркуляции и синоптических вихрей на ЕТР. Установлено наличие высокой положительной корреляции зимних осадков с зональным ветром в северной и центральной части ЕТР.

Ключевые слова:

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

Североатлантическое колебание (САК) является одной из важнейших характеристик крупномасштабной циркуляции атмосферы в северном полушарии и представляет собой разность давления между Азорским максимумом и Исландским минимумом. Сущность САК заключается в перераспределении атмосферных масс между Арктикой и субтропической Атлантикой, при этом переход из одной фазы САК в другую вызывает большие изменения в поле ветра, переносах тепла и влаги, в интенсивности, количестве и траекториях штормов и т.д. При этом, чем выше индекс САК, тем более интенсивным является геострофический зональный перенос воздушных масс с акватории Северной Атлантики на Европейский континент и тем более ослабленной меридиональная циркуляция.

Довольно подробная сводка современного состояния изученности САК и его взаимосвязей с процессами в океане и атмосфере дается в монографии Е.С. Нестерова [9]. При этом закономерности в развитии САК и его связей с гидрометеорологическими процессами в Северной Атлантике установлены преимущественно на качественном уровне и нуждаются в дополнительной проверке и уточнении. Данная работа посвящена выявлению частотной структуры САК, уточнению некоторых

* Работа выполнена при финансовой поддержке Р<

связей САК с характеристиками климата Северной Атлантики и влиянию САК на межгодовую изменчивость увлажнения Европейской территории России.

Существует несколько подходов к расчету индекса САК. Очевидно, наиболее точным является предложенное в работе [11] двухпараметрическое описание САК через модуль градиента давления в Азорском и Исландском центрах действия атмосферы и угол отклонения его от меридиана, которое позволяет более точно определить не только интенсивность, но и направление циркуляции (геострофического ветра). К сожалению, данный индекс рассчитан только до 1984 г. и в настоящее время не поддерживается. Остальные индексы САК являются однопараметрическими. Это разность давления между Азорскими о-вами (Понта-Де-льгада) и Исландией (Акюйрейри или Ре-кьявик), между Лиссабоном (Португалия) и Стиккисхоульмюр (Исландия) [9; 17; 18], а также первые главные компоненты поля атмосферного давления над Северной Атлантикой [12] и высоты геопотенциальной поверхности 500 гПа [15]. В связи с этим рассмотрим соответствие между различными вариантами задания индекса САК. Были взяты 3 индекса: первая главная компонен- ™ та разложения поля высоты геопотенци- ю альной поверхности 500 гПа над Северной го

о

И (грант 14-05-00837). ¿^

Атлантикой (САКГК) [15], разность давления между Понта-Дельгада и Рекьявиком (САКПР) [18] и Лиссабоном и Стик-кисхоульмюром (САКЛС) [17]. При этом рассматривались как средние годовые, так и зимние (декабрь-февраль) значения индекса САК, ибо, по мнению, многих исследователей именно его зимние значения являются наиболее репрезентативными.

В табл. 1 приводятся выборочные оценки коэффициентов корреляции для различных вариантов индекса САК. Прежде всего, отметим, что все коэффициенты являются значимыми на уровне 1%. Для годовых значений связь, близкая к линейной, отмечается между САКПРгод и САКЛСгод. Существенно слабее выражена корреляЦия между САКГКгод и САКЛСгод. (Г = 0,77). Для зимних значений высокая связь отме-

Таблица 1

Распределение выборочных коэффициентов корреляции между различными вариантами индекса САК для средних годовых и зимних значений за период 1950-2012 гг.

САКГК год САКПР год САКЛС год САКГКзим САКПРзим

САКГК год 1,00

САКПР год 0,80 1,00

САКЛС год 0,77 0,95 1,00

САКГК 0,41 0,54 0,39 1,00

САКПРзим 0,38 0,43 0,41 0,39 1,00

САКЛСзим 0,38 0,56 0,41 0,92 0,33

чается только между Между зимними и годовыми значениями САК корреляция колеблется в пределах г = 0,38-0,56.

Так как априори ни один вариант САК предпочтения не заслуживает, то исследователь может манипулировать ими по своему усмотрению и тем самым добиваться нужного результата. Довольно часто в практических расчетах зимние значения САК сравниваются с годовыми значениями гидрометеорологических характеристик (например, с температурой воздуха или осадками). Однако даже получение в этом случае высоких корреляций еще не означает существование зависимости между переменными, ибо коэффициент корреляции -

САК

ПРгод

САК

ЛСгод*

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

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

Рис. 1. Особенности развития гидрометеорологических процессов над Северной Атлантикой для положительной (верхняя панель) и отрицательной (нижняя панель) фаз североатлантического колебания [16].

САИдс5-лвгос (1880 - 2012) = -0.0004 I + 0.0688 = 0.0045

- • - • - • - основной тренд,

----локальные тренды для промежутков времени, указанных в таблице 3.

Рис. 2. Среднегодовые значения (1) североатлантического колебания (САКЛС) и его 5-летнее скользящее среднее (2) за период 1880-2012 гг.

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

На рис. 1 приводится схема развития гидрометеорологических процессов над Северной Атлантикой для положительной и отрицательной фаз САК согласно работе [16]. При положительной фазе САК мы видим, что сравнительно теплый воздух выносится с американского континента на еще более теплую поверхность океана. Далее он в системе мощного зонального переноса достигает европейского континента, вследствие чего там формируются положительные аномалии температуры воздуха. Одновременно с этим над восточной Канадой и Лабрадором отмечаются значительные отрицательные аномалии температуры воды и воздуха и создаются тяжелые ледовые условия в море Лабрадор.

В работах [8; 9 и др.] представлено более подробное обобщение гидрометеорологических процессов в атмосфере и океане, характерных для разных фаз САК. Однако в большинстве случаев связи между фазами САК и фазами развития гидрометеорологических процессов носят качественный характер и не всегда со статистической

точки зрения достоверны. Кроме того, обобщение этих связей выполняется для всех случаев САК > 0 и САК < 0.

На наш взгляд, следует исключать нормальное состояние САК, когда индекс принимает значения, близкие к нулю. В результате возникает задача разбиения временного ряда на однородные градации (норма, выше и ниже нормы). В качестве нормы в данной работе были приняты значения САК в диапазоне -0,5а < САК < 0,5а, где а - стандартное отклонение исходного временного ряда. В табл. 2 приводятся выборочные оценки а для различных вариантов задания САК для годовых и зимних условий. Итак, если значения САК окажутся в пределах ±0,5а, они должны быть исключены при оценке фаз САК.

Таблица 2 Выборочные оценки стандартного отклонения для различных вариантов задания САК для годовых и зимних условий за период 1950-2012 гг.

Индекс САК Годовые значения САК Зимние (декабрь-март) значения САК

САК(ГК) 0,36 0,67

САК(ПР), гПа 0,54 1,17

САК(ЛС), гПа 0,55 2,18

о О

194 Таблица 3

Статистические характеристики фаз усиления (ослабления) индекса САК, сглаженного 5-летним скользящим средним. Данные приведены к середине окна

Характеристика фаза усиления САК фаза ослабления САК

1882-1924 гг. 1969-1992 гг. 1925-1968 гг. 1993-2010 гг.

Величина тренда 0,012 0,022 -0,010 -0,006

Коэффициент детерминации 0,30 0,39 0,38 0,06

На

рис. 2 приводится межгодовой ход

индекса САКЛС с 1880 по 2012 гг. и его 5-летнее скользящее среднее. Нетрудно видеть, что «вековой» тренд практически отсутствует. Это связано с характерными особенностями изменчивости поля атмосферного давления, которая быстро уменьшается с увеличением масштаба временного осреднения. Она весьма значительна для синоптических процессов, уменьшается для сезонных колебаний и становится малой для климатических колебаний. Именно

[1; 7], который обладает целым рядом преимуществ перед традиционными спектральным и гармоническим анализом. На рис. 3 представлено вейвлет-разложение среднемесячных значений индекса САКЛС за период 1875-2010 гг., из которого следу-

ет, что частотная структура

САК

на мас-

штабах менее 10 лет в основном состоит из случайных колебаний.

Из цикличностей наиболее сильным является квазидвадцатилетнее колебание,

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

За рассматриваемый период можно выделить 4 фазы: две - усиление САК и две - его ослабление (табл. 3). Для первых трех из них отмечаются значимые тренды с высоким коэффициентом детерминации. Усиление САК происходило в 1880-1924 и 1969-1992 гг., а ослабление - в 19251968 и 1993-2010 гг. Возможно, начавшаяся в 1993 г. фаза ослабления САК еще не закончилась.

Для выделения и оценивания циклических колебаний в гидрометеорологических рядах целесообразно использовать вейвлет-анализ

120->

110

100-^^В

90 Г

80

70-

60

50-

40-

30-

I Ш

га)«

тИ11ШШ}\

ш.......

ж

........,,Л

_I Щ [т?м11( //А

(Г////

Рис. 3. Распределение коэффициентов вейвлет-разложения среднемесячных значений САКЛС за период 1875-2010 гг. на больших

масштабах.

амплитуда и период которого заметно изменяются во времени. В начале прошлого столетия его амплитуда была максимальной, а период уменьшался до 14-15 лет. Можно выделить также слабое квазисорокалетнее колебание с примерно одинаковой амплитудой и уменьшением периода к середине XX в. Наконец, явственно прослеживается 1 цикл 80-летнего колебания, которое на рис. 2 отождествляется с фазами усиления и ослабления САК. Однако длина временного ряда САК не позволяет проверить длину следующего долгопериодного колебания.

Принято считать, что при положительной фазе САК происходит усиление ветров в тропосфере в зоне 50-60° с.ш., а при отрицательной - наоборот, ослабление ветров в тропосфере в зоне 50-60° с.ш. и усиление меридионального типа циркуляции [9]. Очевидно, что увеличение индекса САК должно сопровождаться углублением Исландского минимума, усилением процессов циклогенеза, особенно в Норвежском море, вследствие чего более мощные циклоны в системе зонального переноса будут приходить на ЕТР и оставлять здесь свою влагу.

В работе [5] было показано, что осадки на ЕТР имеют высокую положительную корреляцию с зональным вертикально-интегрированным _горизонтальным потоком водяного пара _РХ. Одновременно с этим было установлено [4], что над акваторией океана величина вертикально-интегрированного горизонтального потока водяного пара может быть представлена в виде простой формулы:

^ = иЭф \qclp = иэф¥

где и - горизонтальный вектор эффективной скорости переноса водяного пара, W - интегральное по высоте влагосодер-жание атмосферы, q - удельная влажность воздуха. Эффективная скорость ветра с достаточно высокой точностью аппроксимируется скоростью ветра на уровне 500 гПа: иА = (20.536и_п„ + 145.06)0-5 - 10.37,

эф у 500 ' 7

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

В связи с этим для характеристики зонального переноса мы выбрали зональный ветер на уровне 500 гПа на меридиональном разрезе 5-6о в.д. от 40о до 75о с.ш. Исходными данными послужили результаты реа-нализа NCEP/DOE Reanalysis (Reanalysis-2) (http://www.cpc.ncep.noaa.gov/products/wesley/

reanalysis2/index.html). Среднемесячные значения зональной скорости ветра усреднялись для прямоугольников 1о по долготе и 5о по широте. В результате были получены временные ряды зональной составляющей скорости ветра с 1979 по 2013 гг. для 7 центральных точек (42.5о, 47.5о с.ш. и т.д.). Отметим, что среднемесячные значения скорости ветра характеризуют лишь адвективный поток, т.е. за счет средней циркуляции атмосферы и не учитывают его макротурбулентную составляющую, обусловленную перемещением крупномасштабных синоптических вихрей (циклонов и антициклонов). Однако в рассматриваемом регионе движение циклонов направлено преимущественно с запада на восток, причем с усилением зональной циркуляции интенсивность циклонической активности возрастает. Действительно, в соответствии с правилом ведущего потока центры синоптических вихрей перемещаются в направлении устойчивого воздушного течения на уровне 4-6 км над ними [2]. Следовательно, адвективный перенос воздушных масс косвенно учитывает и циклоническую активность.

В табл. 4 представлены коэффициенты корреляции годовых и зимних (декабрь-март) значений индекса САК с соответствующими временными рядами зональной скорости ветра на этом меридиональном разрезе, а также даются оценки скорости ветра для положительной и отрицательной фазы состояния САК. Как видно из табл. 4 индекс САК имеет высокую корреляцию как с зимними, так и с годовыми значениями зональной скорости ветра на уровне 500 гПа в зоне 50-70° с.ш., причем зимой в зоне 50-65о с.ш. она существенно выше. Усреднение скоростей ветра для положительной и отрицательной фаз САК показало, что в первом случае они значительно выше, особенно в зимний период (например, в зоне 55-60о с.ш. более чем в 2 раза). При годовом усреднении различия между скоростями ветра для рассматриваемых фаз САК существенно уменьшаются, но, тем не менее, остаются значимыми по критерию Стьюдента (при а = 0,05) в пределах широтной зоны 55-70о с.ш.

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

о

в

196 Таблица 4

Распределение коэффициентов корреляции индекса САКгк с зональной скоростью ветра на меридиональном разрезе 40-75° с.ш. и оценки усредненной скорости ветра для положительной и отрицательной фазы состояния САК за период 1979-2012 гг. для годовых и зимних (декабрь-март) условий*

Широтная зона, с.ш. Зима (декабрь-март) Год

R + фаза САК (18 случаев) Средняя скорость ветра, м/с - фаза САК (5 случаев) Средняя скорость ветра, м/с R + фаза САК (13 случаев) Средняя скорость ветра, м/с - фаза САК (10 случаев) Средняя скорость ветра, м/с

40-45 -0,51 6,67 9,36 -0,41 8,16 8,88

45-50 0,18 8,59 6,63 -0,16 7,93 8,16

50-55 0,69 12,15 5,98 0,31 9,66 8,81

55-60 0,76 12,79 5,53 0,55 9,76 8,28

60-65 0,67 12,13 6,73 0,61 9,39 7,63

65-70 0,43 10,24 7,92 0,51 8,30 7,09

70-75 -0,03 7,82 8,21 0,31 6,99 6,44

* Значимые (при а = 0,05) различия в скорости ветра между положительной и отрицательной фазами САК выделены полужирным шрифтом.

вает лимитирующее воздействие на процесс испарения. В практических расчетах очень часто используется различные безразмерные индексы, например, коэффициент увлажнения Р/Е , обратная ему величина - коэффициент испаряемости Е0/Р, коэффициент испарения Е/Р, а также относительное испарение Е0-Е и др. Наиболее универсальным параметром является индекс Е0/Р, который посредством различных полуэмпирических уравнений связан практически со всеми остальными характеристиками увлажнения [3; 13]. Анализ межгодовой изменчивости коэффициента испаряемости для ЕТР показал, что в первом приближении можно ограничиться учетом лишь величины осадков, изменчивость которых существенно превосходит изменчивость испаряемости [5].

Источником среднемесячных данных по осадкам послужил метеорологический архив ВНИИГМИ-МЦД [14], содержащий информацию о срочных инструментальных наблюдениях почти на 500 метеорологических станциях России в основном с 1950 года, когда сеть станций резко выросла. Временные ряды по осадкам в этом архиве приводятся с 1966 года, когда на метеорологических станциях был закончен переход к осадкомеру Третьякова и с января 1966 г. начала вводиться поправка на смачивание. После 1966 года никаких изменений в методиках измерений и обработки не происходило, поэтому ряды сумм осадков можно считать однородными [10]. Всего в работе использовались данные 112 станций, расположенные на ЕТР севернее 48о с.ш. Наиболее низкая плотность сети отмечается для северо-восточной части ЕТР.

На рис. 4 представлено распределение изокоррелят между зональной скоростью ветра на уровне 500 гПа в широтной зоне 55-60° с.ш. по меридиану 5о в.д. и значениями осадков для зимних (декабрь-март) условий за период 1979-2012 гг. Значимые коэффициенты корреляции соответствуют условию |г| > 0,30. Нетрудно видеть, что вся северная и центральная часть ЕТР охвачена высокой положительной корреляцией, т.е. при усилении САК количество выпавших здесь осадков увеличивается. Максимальный очаг корреляции (г > 0,70) простирается вдоль меридиана 40о в.д. на территории Архангельской и Вологодской областей, а также отмечается в восточной части Ленинградской области. Как и следовало ожидать, южные районы ЕТР малочувствительны к изменениям САК, т.к. формирование осадков здесь в значительной степени определяется средиземноморскими циклонами.

Что касается корреляции зимних осадков непосредственно с индексом САК, то она ожидаемо оказалось несколько хуже. Максимальные значения корреляции не превысили r = 0,55, причем они сосредоточены в зоне 60-62о с.ш. в центральной части ЕТР.

Кроме того, было выполнено разбиение зимних осадков по фазам САК. При этом число случаев осадков, соответствующих положительной фазе САК, оказалось равным 23, а отрицательной фазе - 11. На рис. 5 дается пространственное распределение критерия Стьюдента (t), показывающего степень расхождения между средними значениями зимних осадков двух выборок [6]. Если t > t (при уровне значимости а = 0,05 t = 2,(5*4), то расхождения между средними значениями

Рис. 4. Распределение изокоррелят между зональной скоростью ветра на уровне 500 гПа в широтной зоне 55-60° с.ш. по меридиану 5о в.д. и значениями осадков для зимних условий (декабрь-март) за период 1979-2012 гг. Значимые коэффициенты корреляции при \г\ > 0.30 (жирная линия).

Рис. 5. Пространственное распределение критерия Стьюдента (), показывающего степень расхождения между средними значениями зимних (декабрь-март) осадков двух выборок, соответствующих положительной и отрицательной фазе развития САК. Критическое значение статистики Стьюдента при уровне значимости а = 0,05 ^ = 2,0.

о

считаются значимыми. Примем для и квазисорокалетнее колебание. Однако простоты t = 2,0. Как видно из рис. 5, вклад их в дисперсию исходного ряда неизолиния t = 2,0, исключая северо- велик, т.е. в значительной степени колеба-западный район ЕТР, имеет в основном ния САК носят случайный характер. зональную протяженность и проходит Установлено, что индекс САК имеет высо-вблизи 59о с.ш. Севернее ее расхождения кую корреляцию, как с зимними, так и с годо-между средними значениями зимних выми значениями зональной скорости ветра осадков значимы, т.е. влияние на них САК на уровне 500 гПа на меридиане 5° в.д. в зоне можно считать установленным. При t = 2,0 50-70° с.ш., причем зимой в зоне 50-65° с.ш. разность осадков между положительной она существенно выше. При положительной и отрицательной фазой САК составляет фазе САК происходит значительное усиле-20-22 мм или 20 % от средней величины ние зональной циркуляции на уровне 500 зимних осадков, а для максимального зна- гПа в зоне 50-70° с.ш., что приводит к уве-чения t = 4.03 (Няндома) - 69 мм или 36%. личению переноса водяного пара в системе Выполненные расчеты позволили вы- средней циркуляции и синоптических вих-явить закономерности межгодовой измен- рей на ЕТР. Вследствие этого вся северная чивости САК. Показано, что при отсутс- и центральная часть ЕТР охвачена высокой твии основного тренда можно выделить положительной корреляцией зимних осад-4 фазы: две - усиление САК и две - его ков с зональным ветром, т.е. при усилении ослабление, для которых отмечаются зна- САК количество выпавших здесь осадков чимые тренды с высоким коэффициентом увеличивается. Максимальный очаг корре-детерминации. При этом начавшаяся в ляции (г > 0.70) простирается вдоль мери-1993 г. фаза ослабления САК еще не за- диана 40° в.д. на территории Архангельской кончилась. В частотной структуре САК ос- и Вологодской областей, а также отмечается новными являются квазидвадцатилетнее в восточной части Ленинградской области.

Список литературы:

[1] Белоненко Т.В. Северо-западная часть Тихого океана. Исследование изменчивости уровня океана в системе вод Куросио-Ойясио на основе спутниковой альтиметрической информации.- Lambert Academic Publishing, 2010. - 218 p.

[2] Зверев А.С. Синоптическая метеорология.- Л.: Гидрометеоиздат, 1977. - 711 с.

[3] Зубенок Л.И. Испарение на континентах.- Л.: Гидрометеоиздат, 1976. - 264 с.

[4] Малинин В.Н. Влагообмен в системе океан-атмосфера.- Л.: Гидрометеоиздат, 1994. - 197 с.

[5] Малинин В.Н. Проблема прогноза уровня Каспийского моря.- СПб.: Изд. РГГМИ, 1994. - 160 с.

[6] Малинин В.Н. Статистические методы анализа гидрометеорологической информации.- СПб, Изд. РГГМУ, 2008. - 407 с.

[7] Малинин В.Н., Гордеева С.М. О современных изменениях глобальной температуры воздуха // Общество. Среда. Развитие. - 2011, № 2. - С. 215-221.

[8] Нестеров Е.С. О влиянии североатлантического колебания на температуру поверхности океана // Метеорология и гидрология. - 1992, № 5. - С. 62-68.

[9] Нестеров Е.С. Североатлантическое колебание: атмосфера и океан.- М.: Триада ЛТД, 2013. - 144 с.

[10] Описание массива данных месячных сумм осадков на станциях России // Российский гидрометеорологический портал: ВНИИГМИ-МЦД. - Интернет-ресурс. Режим доступа: http://meteo.ru/data/158-total-pгecipitation#описание-массива-данных

[11] Радикевич В. М., Энок Ийамуремье. Типизация барического поля для Северной Атлантики и описание Северо-Атлантического колебания (САК) // Современные проблемы гидрометеорологии: сборник научных трудов. Вып. 123. - СПб.: Изд-во РГГМУ, 1999. - С. 43-56.

[12] Смирнов Н.П., Воробьев В.Н., Качанов С.Ю. Североатлантическое колебание и климат.- СПб.: Изд-во РГГМУ, 1998. - 121 с.

[13] Смирнов Н.П., Малинин В.Н. Водный баланс атмосферы как гидрологическая задача.- Л.: Изд. ЛГУ, 1988. - 200 с.

[14] Специализированные массивы для климатических исследований // Российский гидрометеорологический портал: ВНИИГМИ-МЦД. - Интернет-ресурс. Режим доступа: http://aisori.meteo.ru/ClimateR

[15] Barnston A. G., Livezey R. E. Classification, seasonality and persistence of low-frequency atmospheric circulation patterns // Mon. Wea. Rev. - 1987, № 115. - P. 1083-1126.

[16] Bradbury J. A. Cameron P. W. El Nko, the North Atlantic Oscillation and New England Climate. Winter Season Teleconnections and Climate Prediction // Background Information on New England's Climate and Air Quality. AIRMAP Project. Electronic text. - Интернет-ресурс. Режим доступа: http://airmap.unh. edu/background/nao.html, free.

[17] Hurrell J.W., Kushnir Y., Ottersen G., Visbeck M. An Overview of the North Atlantic Oscillation // The North Atlantic Oscillation: Climatic Significance and Environmental Impact. Geophysical Monograph, 134. - AGU, 2003. 10.1029/13 4GM01

[18] Jones P. D., Jonsson T., Wheeler D. Extension to the North Atlantic Oscillation using early instrumental pressure observations from Gibraltar and South-West Iceland // Int.J.Climatol. - 1997, № 17. - P. 1433-1450.