The thermohaline structure of the North Atlantic waters in different phases of the Atlantic multidecadal oscillation

Full Text

1. ВВЕДЕНИЕ

В циркуляции Северной Атлантики (СА) хорошо выделяются естественные долгопериодные квазициклические колебания с типичным периодом 50–70 лет [1, 2]. Это явление получило название Атлантической мультидекадной осцилляции (АМО) [3]. Временной ход индекса АМО рассчитывается как аномалия температуры поверхности океана (ТПО), осредненной по акватории СА от экватора до 70° N [4], из которой удален линейный тренд.

В ряде работ выявлено воздействие АМО на климатические условия в Северном полушарии. Так АМО может оказывать влияние на изменения температуры воздуха, осадков и стока рек в Северной Америке, Евразии и Арктике [4–8], на активность ураганов [9] в СА, а также на атмосферные переносы тепла и влаги в Атлантико-Европейском регионе [7].

В настоящее время механизм формирования мультидекадной изменчивости в Северной Атлантике является предметом дискуссии. Это обусловлено тем, что роли атмосферы и океана в наблюдаемой изменчивости климата поняты не полностью [10]. В настоящее время предложено несколько возможных механизмов для объяснения естественных мультидекадных колебаний климата в Арктике и СA. В работе [11] предполагается, что колебания климата связаны с колебаниями интенсивности меридиональной циркуляции в СА — так называемой Атлантической меридиональной опрокидывающейся циркуляции (АМОЦ). При повышении интенсивности меридиональной циркуляции в верхних слоях происходит перенос более теплых вод в высокие широты СА. Однако эти более теплые и, следовательно, легкие воды в высоких широтах СА замедляют меридиональную циркуляцию, что, в свою очередь, приводит к охлаждению и переход в противоположную фазу колебаний.

В работе [12], авторы указывают на то, что механизм генерации мультидекадной изменчивости связан с взаимодействиями на границе океан–атмосфера и заключается в реакции атмосферы на аномалии ТПО. Одним из возможных механизмов поддержания АМО является сдвиг фаз между долгопериодными изменениями потоков тепла на границе раздела океан–атмосфера в высоких широтах СА и величинами меридиональных переносов тепла в Субтропической Атлантике, определяемый временем адвекции термохалинных аномалий из Тропической Атлантики в субарктические широты [13–18].

В работах [19, 20] на основе результатов длительных расчетов модели земной системы мультидекадные колебания связываются с колебаниями в Арктике. Альтернативная гипотеза состоит в том, что колебательный процесс поддерживается распространением волн типа Россби (солености или плотности волн Россби) вдоль береговой линии [21, 22]. Предлагаются также и другие механизмы, которые дополняют упомянутые выше.

Таким образом, АМО представляет собой ярко выраженный климатический сигнал мультидекадного и декадного масштабов, проявляющийся в ряде климатических характеристик системы океан–атмосфера. Природа собственной мультидекадной изменчивости термохалинного состояния в СА до сих пор до конца не изучена. Поэтому цель настоящей работы заключается в том, чтобы проанализировать характер связей АМО и зонально-осредненных характеристик термохалинного состояния СА по современным, длительным данным объективных океанских анализов EN4 [23, 24] и WOA2013 [25, 26].

2. ВРЕМЕННОЙ ХОД ИНДЕКСА АТЛАНТИЧЕСКОЙ МУЛЬТИДЕКАДНОЙ ОСЦИЛЛЯЦИИ

В эволюции индекса АМО (рис. 1а) с 1940 г. хорошо выделяется мультидекадное квази-шестидесятилетнее колебание [6], на фоне которого развивается декадная изменчивость, скорее всего имеющие разную физическую природу. Гипотезы о формировании мультидекадной изменчивости приведены во введении. Природа декадной изменчивости в бо́льшей степени связана с процессами взаимодействия атмосферы и океана и волнами Россби в Северной Атлантике и Северном Ледовитом океане [21, 22]. В работе [27] показывается, что роль АМОЦ в формировании декадных колебаний несущественна. В нашей работе мы будем исследовать мультидекадную изменчивость термохалинной циркуляции в СА.

 

Рис. 1. (а) ступенчатая линия — значения среднегодового индекса АМО, рассчитанные в NOAA [4] и плавная кривая — результат его низкочастотной фильтрации с отсечением периодов до 8 лет. Штриховыми линиями показаны величины ±0.1 °С, отсекающие периоды, когда индекс АМО принимает значения по модулю больше 0.1 °С. (б) средние тренды приповерхностной температуры Северной Атлантики на глубине 5 м по данным EN4 (сплошная линия) и WOA2013 (штриховая линия) (см. ниже).

 

Будем рассматривать термохалинное состояние СА в различные фазы временного хода индекса АМО только с 1951 г. Это объясняется тем, что наиболее точно трехмерное термохалинное состояние океана по данным наблюдений можно определить только с конца 1940-х гг., поскольку именно с этого времени происходит быстрый рост количества наблюдений в океане [28].

Для изучения климатической изменчивости термохалинной циркуляции СА выделим следующие периоды в соответствии с эволюцией АМО (рис. 1а), по которым будем проводить осреднение по времени для исследования соответствующих климатических состояний СА. Три теплых периода, когда наблюдаются положительные значения индекса АМО больше 0.1 °С с 1951 по 1959 гг., с 2000 по 2008 гг. и с 2009 по 2017 гг.; переходный период от теплой к холодной фазе — с 1961 по 1969 гг.; два холодных периода, когда наблюдаются отрицательные значения индекса АМО меньше –0.1 °С с 1973 по 1981 гг. и с 1982 по 1990 гг.; переходный период от холодной к теплой фазе — с 1991 по 1999 гг. Интервалы этих периодов составляют примерно 10 лет для того, чтобы при осреднении исключить влияние декадных колебаний, оставляя только фазы хорошо выраженного квази-шестидесятилетнего колебания, хорошо заметного на рис. 1а. При этом размах этих колебаний АМО и средних изменений ТПО за 70 лет (рис. 1б) сопоставимы.

3. ИСПОЛЬЗУЕМЫЕ ДАННЫЕ ПО ТЕМПЕРАТУРЕ И СОЛЕНОСТИ НА АКВАТОРИИ СЕВЕРНОЙ АТЛАНТИКИ

Для определения и сравнения термохалинного состояния Северной Атлантики нами были выбраны глобальные поля температуры и солености из двух объективных анализов данных наблюдений EN4 [23] (с поправками В.В. Гурецкого [24]), созданным в MetOffice (Великобритания) и World Ocean Atlas 2013 (WOA2013) [25, 26] из National Oceanic and Atmospheric Administration (NOAA), США. EN4 явился продуктом двух Европейских проектов — Enhanced Ocean Data Assimilation and Climate Prediction (ENACT) и ENSEMBLES [29]. Описание рассматриваемых массивов данных приведено в табл. 1. Следует отметить, что основным их источником служит World Ocean Database (WOD, который получен NOAA), в случае EN4 это WOD09 [30], а в случае WOA13–WOD13 [31].

 

Таблица 1. Описание объективных анализов EN4 и WOA2013

	EN	WOA
Версия	4.2.0	2013 v2
Разработчик	MetOffice (UK)	NODC/NOAA (USA)
Периоды	1900–2017	1955–2012
Шаг по времени	1 месяц	10 лет (для каждого месяца года)
Пространственное разрешение	1° × 1°	0.25° × 0.25°
Число уровней	42	102
Нижний уровень, м	5350	5500
Источник данных	WOD09, GTSPP, ARGO, ASBO	WOD13

 

Объективные анализы EN4 и WOA2013 рассчитаны с использованием пространственно-временной оптимальной интерполяции данных наблюдений. Указанные массивы данных также получены с использованием комбинированного подхода: учитывалось спрогнозированное состояние океана, основанное на более ранних данных объективного анализа и профили данных наблюдений за доступный временной период (см. табл. 1).

Для анализа будем использовать потенциальную температуру (θ), как более наглядно, чем температура in situ (T), отражающую тепловое состояние океана по глубине. По массивам среднемесячных данных θ и солености (S) из EN4 с 1951 по 2017 гг. были рассчитаны их среднегодовые трехмерные поля для каждого рассматриваемого периода (рис. 1а). Массив данных WOA2013 содержит среднемесячные, среднесезонные и среднегодовые трехмерные поля по температуре in situ и солености уже осредненные для следующих 6 периодов: с 1955 по 1964 гг., с 1965 по 1974 гг., с 1975 по 1984 гг., с 1985 по 1994 гг., с 1995 по 2004 гг., с 2005 по 2012 гг. Для анализа были использованы среднегодовые значения для этих периодов. По температуре in situ данных WOA2013 были рассчитаны трехмерные поля потенциальной температуры с помощью методики [32].

На рис. 1б приведены линейные тренды приповерхностной температуры в СА на глубине 5 м. Положительный тренд наблюдается как по данным EN4, так и по данным WOA2013, с коэффициентами 0.101 °C/10 лет и 0.120 °C/10 лет соответственно. Небольшая разность между ними примерно в 0.02 °C/10 лет связана, по-видимому, с разной длиной временных рядов, выбранных для анализа. Массив WOA2013 имеет более короткую продолжительность и, следовательно, вклад периода наиболее интенсивного потепления с конца 1970-х по начало 2000-х гг. в этих данных сказывается сильнее.

4. ХАРАКТЕРИСТИКИ КЛИМАТИЧЕСКИХ ТРЕНДОВ И КОМПОЗИТНЫЙ АНАЛИЗ АНОМАЛИЙ ТЕМПЕРАТУРЫ И СОЛЕНОСТИ

В данной работе применялся композитный анализ [33] аномалий среднегодовых среднезональных характеристик θ, S и потенциальной плотности, характеризующих термохалинный сигнал климатических изменений в различные фазы АМО. Композитный анализ или метод эмпирических условных статистик — это общепринятая методика при анализе изменений климата, он используется для характеристики изменчивости какого-либо параметра в противоположных фазах. С помощью этой методики можно выделить чистый сигнал, соответствующий влиянию отдельного климатического сигнала на гидрофизические поля. Композитный анализ используется во многих работах, в качестве примера можно привести работы [4, 33, 34], в которых исследуется связь АМО с изменчивостью приповерхностных характеристик атмосферы и океана.

Аномалии вычислялись путем вычитания из полученных трехмерных полей θ и S линейных трендов, определяемых за рассматриваемые периоды EN4 и WOA2013. Вычитание линейных трендов проводилось для того, чтобы исключить более долгопериодные климатические изменения климата океана, связываемые в основном с внешним, главным образом антропогенным, влиянием [35], и оставить только собственные изменения климата.

Для того, чтобы оценить сигнал климатической изменчивости в СА проводился композитный анализ картины аномалий θ и S по отношению к фазам индекса АМО. Композитный анализ аномалий какой-либо характеристики основан на расчете разности между композитом аномалий этой характеристики в положительную фазу индекса АМО и композитом аномалий в отрицательную фазу индекса АМО (так называемая разность композитов). При этом композит аномалий рассматриваемой характеристики в положительную фазу индекса АМО, определяется как усредненное за все периоды состояние этих аномалий, когда фаза индекса АМО положительна (для EN4 периоды — с 1951 по 1959 гг., с 2000 по 2008 гг., с 2009 по 2017 гг., в случае WOA2013 — с 2005 по 2012 гг.), а композит аномалий в отрицательную фазу индекса АМО, как усредненное за все периоды состояние аномалий, когда фаза индекса АМО отрицательна (с 1973 по 1981 гг., с 1982 по 1990 гг. для EN4, в случае WOA2013 — с 1975 по 1984 гг.). Разность композитов рассматриваемой характеристики рассчитывалась в каждом узле пространственной сетки и затем зонально осреднялась. Таким образом, полученная зонально-осредненная разность композитов аномалий рассматриваемой характеристики отражает среднезональный сигнал естественных мультидекадных колебаний, показывая характер изменений в теплой фазе по отношению к холодной фазе АМО.

Прежде, чем перейти, собственно, к анализу разностей композитов, рассмотрим среднезональные тренды θ и S, которые также представляют важный интерес с точки зрения климатических изменений, происходящих в СА. На рис. 2 представлены эти зонально-осредненные тренды, рассчитанные по данным EN4 и WOA2013, которые показывают практически одну и ту же пространственную структуру долгопериодных изменений термохалинного состояния СА. Это подтверждается высокими значениями пространственных коэффициентов корреляции между данными EN4 и WOA2013, равными 0.92 для трендов температуры (рис. 2а, б) и 0.77 для трендов солености (рис. 2в, г). При этом данные EN4 показывают более сглаженную структуру долгопериодных изменений, чем WOA2013 (рис. 2), что связано, скорее всего, с более высоким пространственным разрешением данных WOA2013 (см. табл. 1).

 

Рис. 2. Зонально-осредненные линейные тренды потенциальной температуры θ (а, б),  солености S (в, г),  и потенциальной плотности в величинах σ_Т (д, е) (σ_Т = ρ(θ, S)–1000),  рассчитанные по объективным анализам EN4 (слева), WOA2013 (справа). Изолиниями показаны средние этих величин за рассматриваемые периоды.

 

Согласно рис. 2 в верхнем 1-км слое, в основном, наблюдается потепление, особенно сильное около 40° N, где Гольфстрим переходит в Северо-Атлантическое течение. Это потепление в области 40° N связывается со сдвигом положения начала Северо-Атлантического течения к северу, что приводит к сильной аномалии температуры и солености в этих широтах. Ниже 1-км слоя наблюдается значительная область похолодания и распреснения, образование которой связывается с таянием льдов Гренландии, выносом распресненных вод из Северного Ледовитого океана и затягиванием этих холодных и более пресных вод в глубинные слои. С большой долей уверенности можно полагать, что основной вклад в таяние ледников Гренландии, так и в сдвиг Северо-Атлантического течения на север, делает парниковый эффект [35, 36], вызванный, главным образом, антропогенным воздействием [35].

Также были рассчитаны коэффициенты корреляции между трендами температуры и солености, которые составили 0.79 как для данных EN4, так и для WOA2013, что говорит о высокой связи между температурой и соленостью. Высокое значение этого коэффициента говорит и о том, что аномалии θ и S делают противоположный, компенсирующий друг друга вклад в формирование аномалий поля потенциальной плотности.

Рис. 2д, е показывают, что, когда тренды θ способствуют уменьшению/увеличению плотности, тогда же тренды S направлены на ее увеличение/уменьшение. При этом в верхнем 1-км слое отрицательные тренды плотности определяются в основном положительными трендами температуры, а в более глубинных слоях отрицательные тренды плотности формируются главным образом отрицательными трендами солености. Суммарный вклад трендов температуры и солености дает общий тренд к уменьшению плотности по всей глубине. Это говорит о том, что потенциальная энергия вод СА за счет изменения их термохалинного состояния увеличилась за счет расширения вод, в основном в поверхностном слое, что приводит к увеличению стерического уровня океана, что, в свою очередь, согласуется с климатическими расчетами и оценками по данным наблюдений IPCC (Intergovernmental Panel on Climate Change, или в русской транскрипции МГЭИК — Межправительственная группа экспертов по изменению климата) [35].

Распределение температурных трендов (рис. 2а, б) показывает, что в СА в верхнем одно-километровом слое за последние 70 лет накоплено колоссальное количество дополнительного тепла, оцениваемое по используемым данным EN4 в 6.89 ∙ 10⁷ ПДж, а по WOA2013 в 5.63 ∙ 10⁷ ПДж. Ясно, что этот теплозапас заключает в себе такую инерцию климатического потепления, что даже если остановить антропогенные выбросы парниковых газов потепление климата Земли не прекратится, о чем также свидетельствуют климатические оценки IPCC [35]. Так, чтобы скомпенсировать это потепление вод в верхних слоях СА, например, за 50-летний промежуток времени, необходимо охлаждать поверхность океана со средним потоком 0.76 вт/м² согласно данным EN4 и 0.65 вт/м² согласно WOA2013. Однако даже при самом умеренном сценарии IPCC будущего развития климатических изменений, радиационный поток парникового эффекта, направленный на разогрев поверхности Земли, составляет 2.5 вт/м².

Показанные на рис. 2д, е распределения трендов плотности означают, что вертикальная среднезональная стратификация в СА за 70 лет стала более устойчивой. Это должно приводить к уменьшению интенсивности атлантической меридиональной опрокидывающей циркуляции (АМОЦ) в СА, что соответствует мультимодельным расчетам по климатическими моделям, выполненными в рамках международной программы исследования климата IPCC [35, 36].

При анализе разностей композитов аномалий θ (рис. 3а, б) и S (рис. 3в, г), которые отражают сигнал естественных мультидекадных колебаний, была обнаружена очень интересная структура вертикальной стратификации. Она проявляется практически одинаково как для данных EN4, так и для WOA2013. Коэффициенты пространственной корреляции между данными EN4 и WOA2013 равны 0.73 и 0.62 для разностей композитов аномалий θ и S соответственно. Установлено, что в теплые периоды индекса АМО (рис. 1а) по сравнению с холодными наблюдается потепление и осолонение в верхнем 1-км слое Северной Атлантики (рис. 3). Ниже 1-км слоя наблюдаются значительные области похолодания и распреснения, которые более ярко выражены в данных EN4. Косвенно этот эффект был также обнаружен в работе [37] при изучении изменчивости термохалинного состояния в отдельных районах СА. По-видимому, такая структура изменчивости связана с характером изменчивости термохалинной циркуляции СА и требует своего исследования.

 

Рис. 3. Зонально-осредненные разности композитов (показаны цветом) аномалий потенциальной температуры θ (а, б), °С, солености S (в, г), ШПС и потенциальной плотности в величинах σ_Т (д, е) (σ_Т = ρ(θ, S)–1000),  рассчитанные по объективным анализам EN4 (слева), WOA2013 (справа). Изолиниями показаны средние этих величин за рассматриваемые периоды.

 

Также были рассчитаны коэффициенты корреляции между разностями композитов аномалий θ и S, которые равны 0.72 и 0.78 для данных EN4 и для WOA2013 соответственно. Высокие значения этих коэффициентов корреляции, так же как было отмечено и для трендов, означают, что аномалии θ и S вносят противоположный, компенсирующий друг друга вклад в формирование аномалий поля плотности. При этом вклад температурных изменений в изменения потенциальной плотности в верхних слоях превалирует над соленостным, а в глубинных слоях наоборот, наблюдается превалирующий вклад солености, что хорошо согласуется с работой [18].

Разности композитов аномалий потенциальной плотности демонстрируют несколько отличную структуру пространственных изменений для данных EN4 (рис. 3д) и WOA2013 (рис. 3е), особенно проявляющуюся в глубинных слоях. Согласно данным WOA2013 в глубинных слоях больше областей с положительными аномалиями плотности. Структура аномалий плотности по данным WOA2013 отличается и бо́льшей перемежаемостью.

Обнаруженные отличия в термохалинном состоянии СА в периоды положительного и отрицательного индекса АМО привели авторов статьи к гипотезе, что основную роль в этом процессе может играть Атлантическая меридиональная опрокидывающаяся циркуляция (АМОЦ). Пространственная структура АМОЦ, оцениваемая на основе данных наблюдений и воспроизводимая в различных моделях общей циркуляции океана, показана, например, в работах [38–40]. В целом АМОЦ является характеристикой термохалинной циркуляции океана в СА, когда в верхнем 1-км слое воды медленно переносятся на север, отдавая накопленное тепло в атмосферу в средних и высоких широтах и, вследствие этого, опускаются (опрокидываются), возвращаясь в глубинных слоях на юг. Именно АМОЦ и определяет меридиональный перенос тепла на север водами СА. Таким образом движение вод в АМОЦ можно косвенно связать с атлантическим сектором океанского конвейера Брокера–Лаппо [41, 42].

Для подтверждения нашей гипотезы мы будем использовать структуру АМОЦ, рассчитанную нами с помощью хорошо апробированной модели океана INMOM (Institute of Numerical Mathematics Ocean Model) [43]. Эта модель уже на протяжении более полутора десятков лет используется в качестве океанического блока совместной модели земной системы INMCM (Institute of Numerical Mathematic Climate Model) различных версий [44–49]. Именно эта совместная модель является единственным представителем от России в международном проекте сравнения климатических моделей под эгидой IPCC в международных проектах CMIP (Coupled Model Intercomparison Project, http://cmip-pcmdi.llnl.gov). Последнее такое сравнение CMIP5 (CMIP, Phase 5) проходило в 2010–2013 гг. [50]. Кроме того, глобальная версия INMOM, так же единственная от России, участвовала в международном проекте CORE (Coordinated Ocean-ice Reference Experiment) [51], вторая фаза которого была посвящена изучению изменчивости циркуляции в СА [38, 52].

Модель INMOM была настроена на акваторию Северной Атлантики от 100° W до 15° E и от экватора до 80° N с пространственным разрешением 0.5° по долготе и широте. Были подготовлены данные о батиметрии дна по массиву ETOPO5 для этой области и начальные состояния по температуре и солености из EN4 осредненные за период с 1948 по 2017 гг., причем для лета и зимы отдельно. Последнее связано с тем, что невозможно корректно восстановить среднегодовую циркуляцию океана по его среднегодовому термохалинному состоянию. Атмосферное воздействие вычислялось путем осреднения в соответствии с выделенными климатическим периодом приповерхностных атмосферных характеристик из данных CORE (Coordinated Ocean-ice Reference Experiment) [51].

Расчет циркуляции для климатического периода с 1948 по 2017 гг. проводился по методу диагноза–адаптации, предложенный академиком А.С. Саркисяном с соавторами [53]. Запуск модели осуществлялся на 75 суток, из которых первые 30 суток проводился разгон модели в режиме диагноза при предписанных полях температуры и солености. В последующие 45 суток эксперимент проводился в режиме адаптации, при котором рассчитывались итоговые поля термохалинного состояния и циркуляции вод для всех летних и зимних выделенных периодов АМО. Для анализа среднегодового состояния этих периодов результаты расчета осреднялись для зимы и лета и за последние 30 суток расчетов. При таком подходе роль модели вторична, поскольку она главным образом используется для восстановления циркуляции океана по оценке его термохалинного состояния по данным наблюдений, а также устранения на этапе адаптации заведомых ошибок в полях θ и S. Расчеты показали, что эта версия INMOM по предложенному сценарию диагноз–адаптация адекватно воспроизводит среднеклиматическую наблюдаемую структуру основных гидрофизических полей. По результатам расчетов определялась АМОЦ, по стандартной методике, описанной, например, в [43]. На рис. 4 и 5 она показана в виде изолиний. Ее характерные особенности согласуются с климатическими расчетами [38–40]. Как правило в этих расчетах она проявляет индивидуальные черты, свойственные отдельным моделям океана. В нашем же случае АМОЦ, рассчитанная по методу диагноз–адаптация, в первую очередь определяется данными наблюдений. Модель здесь играет вторичную роль [53]. Следует отметить, что наши расчеты АМОЦ очень хорошо согласованы с реанализом, выполненным по модели GFDL, показанным в работе [18].

 

Рис. 4. Зонально осредненные аномалии потенциальной температуры θ в различные периоды АМО, (°С, показана цветом) и наложенная на нее средняя за 1948–2017 гг. функция тока АМОЦ, рассчитанная по модели INMOM (Св, показана изолиниями, стрелки указывают на генерализованное направление движения вод меридиональной циркуляции).

 

Рис. 5. То же, что и на рис. 4, но для аномалий солености S (ШПС).

 

Поскольку композитный анализ выявил более отчетливую вышеописанную характерную структуру изменчивости термохалинных полей в данных EN4, то покажем связь ее формирования с АМОЦ по этим данным. На рис. 4 и 5 показан процесс перемещения среднезональных аномалий θ и S по данным объективного анализа EN4 для выделенных фаз АМО на фоне функции тока АМОЦ, средней за 1948–2017 гг. Она показана изолиниями начиная от 20° N, так как метод диагноза–адаптации работает не достаточно корректно в приэкваториальной области из-за особенностей океанической циркуляции, приводящей здесь к дополнительному шуму в динамике океана [53]. Генерализованные направления движения вод в АМОЦ показаны дополнительными стрелочками. Из анализа зонально-осредненных детрендированных аномалий температуры и солености, показанных на рис. 4 и 5 в чередующиеся последовательные периоды индекса АМО, было обнаружено, что положительные и отрицательные аномалии как температуры, так и солености циркулируют по ходу движения вод в АМОЦ вокруг ее ядра, опускаясь в глубинные слои океана примерно на 60° N и поднимаясь на поверхность на 25° N. Это вполне согласуется с теорией мультидекадных колебаний, предложенной еще в работе [11].

Из этих рисунков видно, что в теплые периоды АМО (рис. 4а, е, ж, рис. 5а, е, ж), а также при переходе от теплого к холодному периоду (рис. 4б, рис. 5б), воды с положительными аномалиями температуры и солености опускаются в глубинные слои ниже 1-км, а воды с отрицательными аномалиями поднимаются из глубинных слоев на поверхность в соответствии с перемещениями вод в АМОЦ [38–40]. В холодные же периоды индекса (рис. 4в, г, рис. 5в, г) и при переходе от холодного к теплому периоду (рис. 4д, рис. 5д), наблюдается обратная картина, когда воды уже с отрицательными аномалиями температуры и солености опускаются, а холодные воды с положительными температурными и соленостными аномалиями поднимаются. И этот процесс происходит приблизительно с периодом 60 лет, также в соответствии с оценкой времен перемещения вод в АМОЦ [38–40]. Можно предположить, что благодаря именно этому процессу наблюдаются теплые и холодные фазы индекса АМО.

При этом нельзя утверждать, что в поверхностном слое аномалии θ и S пассивно переносятся движением вод за счет АМОЦ, так как, во-первых, здесь они вступают во взаимодействие с атмосферой, а, во-вторых, в приблизительно 100-метровом слое в средних широтах 40°–60° N АМОЦ меняет свое направление на обратное.

Обнаруженный нами процесс формирования зонально-осредненных детрендированных аномалий температуры и солености в чередующиеся последовательные периоды индекса АМО косвенно подтверждает положение о том, что одним из возможных механизмов поддержания АМО является сдвиг фаз между долгопериодными изменениями потоков тепла на границе раздела океан–атмосфера в высоких широтах СА и величинами меридиональных переносов тепла в Субтропической Атлантике, определяемый временем адвекции термохалинных аномалий из Тропической Атлантики в субарктические широты [13–18].

5. ЗАКЛЮЧЕНИЕ

В настоящее время актуальны задачи диагноза и прогноза изменений климата и выявления причин, которые эти изменения обуславливают. Атлантический океан является важным звеном климатической системы Земли. Имеются основания полагать, что мультидекадные изменения климата во многом порождаются долгопериодными собственными изменениями в термохалинной циркуляции СА. Эти мультидекадные (главным образом квази-шестидесятилетние) колебания климата хорошо проявляются во временно́м ходе индекса АМО [4].

Для определения и сравнения термохалинного состояния СА в различные фазы АМО по различным массивам данных были выбраны глобальные поля температуры и солености по двум объективным анализам EN4 [23, 24] и WOA2013 [25, 26]. По данным объективного анализа EN4 были рассчитаны средние климатические состояния для трех теплых периодов АМО, для переходного периода от теплой к холодной фазе, двух холодных периодов и для переходного периода от холодной к теплой фазе (рис. 1а). Массив данных WOA2013 содержит данные по температуре и солености уже осредненные для 6 десятилетних периодов, которые также можно отнести к разным фазам АМО.

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

Из анализа зонально-осредненных трендов потенциальной температуры θ и солености S (рис. 2) установлено, что в обоих рассматриваемых массивах данных EN4 и WOA2013, в верхнем 1-км слое, в основном, наблюдается потепление и осолонение. Ниже 1-км слоя наблюдается значительная область похолодания и распреснения, что связывается с таянием льдов Гренландии, выносом распресненных вод из Северного Ледовитого океана и затягиванием этих холодных и более пресных вод в глубинные слои. Такие структуры трендов температуры и солености, наиболее вероятно, являются следствием парникового эффекта, вызванного, согласно парадигме IPCC, в основном антропогенным воздействиям [35, 36]. Показано, что θ и S вносят противоположный, компенсирующий вклад в формирование поля плотности как в долгопериодных трендах, так и мультидекадных колебаниях. Вклад температурных изменений превалирует над соленостным в верхних слоях СА, а в глубинных наоборот, что хорошо согласуется с работой [18]. При этом тренд в потенциальной плотности показывает, что за рассматриваемый промежуток времени потенциальная энергия вод СА увеличилась за счет расширения при изменении их термохалинного состояния.

Распределение температурных трендов также показывает, что в СА в верхнем одно-километровом слое за последние 70 лет накоплено значительное количество тепла и, очевидно, что этот теплозапас заключает в себе инерцию климатического потепления, так, что даже если остановить антропогенные выбросы парниковых газов потепление климата Земли не прекратится [35], о чем также свидетельствуют климатические расчеты по международной программе IPCC [35]. Распределение трендов плотности указывает на то, что вертикальная среднезональная стратификация в СА за 70 лет стала более устойчивой. Это также должно приводить к уменьшению интенсивности меридиональной опрокидывающейся циркуляции (МОЦ) в СА, что соответствует модельным расчетам [35, 36].

Анализируя разности композитов аномалий потенциальной температуры и солености (рис. 3) была обнаружена важная особенность структуры пространственной термохалинной изменчивости, которая наблюдается как для данных EN4, так и для WOA2013. Было установлено, что в теплые периоды индекса АМО (рис. 1а) по сравнению с холодными наблюдаются потепление и осолонение в верхнем 1-км слое Северной Атлантики. Ниже 1-км слоя наблюдаются значительные области похолодания и распреснения, которые более ярко выражены в данных EN4.

Было показано, что в чередующиеся последовательные периоды индекса АМО положительные и отрицательные аномалии как температуры, так и солености циркулируют по ходу движения вод в АМОЦ вокруг ее ядра, опускаясь в глубинные слои океана примерно на 60° N и поднимаясь на поверхность на 25° N. И этот процесс происходит приблизительно с периодом 60 лет. Можно предположить, что благодаря этому наблюдаются теплые и холодные фазы индекса АМО. По-видимому, такая структура изменчивости связана с характером формирования термохалинной циркуляции в СА и требует своего дальнейшего исследования.

Обнаруженный в настоящей работе процесс формирования зонально-осредненных детрендированных аномалий температуры и солености в чередующиеся последовательные периоды индекса АМО косвенно подтверждает положение о том, что одним из возможных механизмов поддержания АМО является сдвиг фаз между долгопериодными изменениями потоков тепла на границе раздела океан–атмосфера в высоких широтах СА и величинами меридиональных переносов тепла в Субтропической Атлантике, определяемый временем адвекции термохалинных аномалий из Тропической Атлантики в субарктические широты [13–18]. Следует отметить, что этот процесс до конца не изучен. Поэтому в дальнейшем для более глубокого исследования предполагается привлечь данные океанских реанализов и климатические расчеты циркуляции океана с помощью российской модели INMOM.

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

 

About the authors

N. A. Diansky

Lomonosov Moscow State University; Marchuk Institute of Numerical Mathematics of the RAS; Zubov State Oceanographic Institute

Author for correspondence.
Email: nikolay.diansky@gmail.com
Russian Federation, Leninskiye Gory, 1, 119991, Moscow; Gubkina st., 8, 119333, Moscow; Kropotkinsky lane, 6, 119034, Moscow

V. A. Bagatinsky

Lomonosov Moscow State University

Email: nikolay.diansky@gmail.com
Russian Federation, Leninskiye Gory, 1, 119991, Moscow

References

Supplementary files