Меридиональный перенос тепла в Северной Атлантике и тенденции его изменений во второй половине XX века

УДК 465.551.7

А.Б. Полонский, С.Б. Крашенинникова

Меридиональный перенос тепла в Северной Атлантике и тенденции его изменений во второй половине XX века

Посвящается 75-летнему юбилею заслуженного деятеля науки Украины доктора географических наук НА. Тимофеева

По данным зональных разрезов WOCE ( World Océan Circulation Experiment\ выполненных в Северной Атлантике в 1992 - 1998 гг., рассчитан меридиональный перенос тепла (МПТ) в океане. Проведен обобщенный анализ оценок МПТ, полученных прямыми методами различными авторами по данным разрезов, выполненных во второй половине XX в. между 7,5 и 48° с.ш. МПТ, осредненный за весь период наблюдений, достигает максимума (1,38 ± 0,19 ПВт) в Субтропической Атлантике. МПТ подвержен достаточно интенсивной сезонной изменчивости. В конце лета в Субтропической Атлантике наблюдается максимум МПТ (около 1,9 ПВт), в конце зимы - минимум (около 0,8 ПВт). Выявлена значимая тенденция увеличения МПТ в окрестности 36° с.ш. между 1959 и 1993 гг. (от 0,75 до 1,1 ПВт), что свидетельствует об интенсификации меридиональной океанической циркуляции в Северной Атлантике.

Введение. Одним из главных механизмов, посредством которых океан влияет на формирование глобальной климатической системы, является перенос тепла океаническими течениями. Особенно важен МПТ в тропических и субтропических широтах, поскольку за счет МПТ в системе океан - атмосфера происходит уменьшение термических контрастов между высокими и низкими широтами, а на океанический МПТ приходится основная доля суммарного переноса в этой широтной полосе. Атлантический океан играет особую роль в формировании глобального МПТ. В нем существует квазистационарный перенос тепла на север [1 - 8]. Он обусловлен глобальной термохалИнной циркуляционной ячейкой, которая образуется за счет конвекции холодных глубинных вод в очагах их формирования в Северной Атлантике и компенсирующего переноса относительно теплых вод термоклина на север [9, 10]. Изменчивость интенсивности меридиональной циркуляции в Атлантическом океане в значительной степени регулирует низкочастотную изменчивость системы океан - атмосфера. Таким образом, корректное моделирование глобального климата и его изменений требует точного знания средних величин МПТ в океане и тенденций их изменений в процессе наблюдаемого глобального потепления. К сожалению, опубликованные данные об МПТ в океане достаточно противоречивы. Они основываются на двух основных методиках оценки МПТ:

- прямые океанографические измерения [И];

© А.Б. Полонский, С.Б. Крашенинникова, 2007 ISSN 0233-7584. Мор. гидрофиз. журн., 2007, № 1

- косвенные оценки с использованием уравнения бюджета тепла для океана [3, 5, 6] либо для системы океан - атмосфера [4, 7, 8].

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

По оценкам разных авторов, кроссэкваториальный перенос тепла в Северной Атлантике составляет порядка 1 ПВт и направлен на север. МПТ в океане достигает максимального значения (до -1,5 ПВт) в Субтропической Атлантике. Причем косвенные оценки, выполненные в последние годы, а также оценки, полученные инверсными методами, в среднем дают меньшие величины МПТ в Северной Атлантике по сравнению с прямыми оценками. Анализ последних с использованием наиболее точных,океанографических данных, полученных по программе 1¥ОСЕ в 90-е годы XX в., подтвердил справедливость большинства выполненных ранее прямых оценок океанического МПТ. Таким образом, между прямыми и косвенными оценками существуют систематические расхождения [17, 19]. Поэтому для анализа тенденции изменения МПТ в Северной Атлантике лучше использовать оценки, выполненные по одной методике, с тем, чтобы исключить эти расхождения.

С целью уточнения оценок МПТ и обнаружения тенденции изменения МПТ в Северной Атлантике в настоящей работе был рассчитан МПТ с использованием массива данных \VOCE и проанализированы все доступные оценки МПТ, полученные прямыми методами по данным зональных разрезов в широтной полосе между 7,5 и 48° с.ш. с 1957 по 1998 гг.

Методика расчета меридионального переноса тепла по прямым океанографическим измерениям. Анализ погрешностей. Квазистационарный среднегодовой МПТ в океане (Н) через круг широты можно рассчитать непосредственно как

Н = Ср^\\

о о

где Ни Ь- глубина и ширина океана; V - меридиональная компонента скорости; Г-температура воды; Срр& 4,18 Дж/(град-см3), Ср - удельная теплоемкость морской воды при постоянном давлении, - плотность морской воды. Черта сверху означает осреднение по времени. Ясно, что для расчета МПТ прямым методом необходимо иметь данные о меридиональной компоненте скорости и температуре на каждом круге широты.

В принципе квазистационарный среднегодовой перенос тепла может быть представлен как сумма четырех следующих МПТ: 1) обусловленный среднегодовой меридиональной циркуляцией; 2) связанный с наличием корреляций скорости и температуры на сезонном масштабе; 3) связанный с наличием корреляций скорости и температуры на синоптическом масштабе;

4) связанный с различными корреляциями между сезонными и синоптическими флюктуациями [14, 19]. В соответствии с результатами, полученными в работах [14, 19], главный вклад в квазистационарный МПТ в Северной Атлантике вносит перенос, формирующийся за счет среднегодовой меридиональной циркуляции, структура которой достаточно устойчива, несмотря на наличие интенсивной сезонной изменчивости. Это дает возможность оценить квазистационарный МПТ, даже не имея достаточного для надежного определения сезонных и синоптических вариаций Уи Г количества разрезов. На это впервые было указано в работе [11], в которой МПТ в океане вообще оценивался по данным одного разреза. Вместе с тем привлечение информации по нескольким разрезам для каждого круга широты позволяет получить оценку МПТ для разных сезонов, уменьшить погрешность вычисления квазистационарного МПТ и оценить долговременную тенденцию его изменения, что и делается в настоящей работе.

Если представить меридиональную скорость как сумму дрейфовой, ба-роклинной и баротропной скоростей, то выражение (1) запишется в виде

Я = Срр]рдрГ +F6kГ +УбтГ \xdz,

о о

Удр — существенна в верхнем экмановском слое, V в к - в пределах главного

термоклина, Vqt - в глубинных слоях океана, Т* =T-TD, TD - температура глубинных вод океана. Дрейфовые переносы оценивались с использованием экмановских соотношений. Бароклинные скорости течений рассчитывались по геострофическим соотношениям относительно нулевой поверхности 2 ООО м. Баротропный перенос вычислялся с использованием закона сохранения массы. Более подробно методика расчета МПТ изложена в работах [14, 19, 20].

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

^общ = ^инст +СГсез +

где <тобщ - суммарная погрешность оценки МПТ, атст - инструментальная погрешность, <7сез - погрешность, связанная с недостаточным разрешением сезонной изменчивости, асин - погрешность, обусловленная синоптическими флюктуациями меридиональной компоненты скорости и температуры, ^межгод ~ погрешность, вызванная межгодовой изменчивостью МПТ. Вместе

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

ISSN 0233-7584. Мор. гидрофиз. журн., 2007, № 1 47

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

Инструментальная погрешность расчета МПТ по данным зональных гидрологических разрезов, выполненных до начала реализации программы WOCE, была много меньше физических погрешностей, несмотря на относительно невысокую точность стандартных океанографических измерений [11, 17]. Это связано с интегральным характером оценки МПТ, что существенно уменьшает случайные погрешности. Разрезы WOCE выполнялись более совершенными зондирующими комплексами, точность измерений которых на порядок выше по сравнению с точностью измерений в 50 - 70-е годы. Поэтому инструментальной погрешностью расчета МПТ по этим данным можно пренебречь по сравнению с физическими погрешностями.

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

деленной на л/2 , так как

где А - амплитуда сезонного хода.

Чтобы описать синоптические возмущения и оценить возможную погрешность, связанную с их недостаточным разрешением на гидрологических разрезах, а также несинхронностью выполнения разреза, мы применяли следующую методику. Для каждой станции на разрезе рассчитываются единичные интегральные переносы массы и тепла, которые затем интерполируются в узлы регулярной сетки с небольшим шагом по долготе, величина которого много меньше типичной длины синоптического возмущения. С помощью спектрального анализа определяются длина волны и амплитуда синоптического возмущения для каждого разреза. Полученные длины волн (порядка сотен километров) характеризуют планетарные волны (волны Россби). Для того чтобы рассчитать период и фазовую скорость синоптического возмущения, отвечающего наиболее часто реализующейся в океане низшей моде, определяется частота первой бароклинной моды (п = 1) из дисперсионного соотношения для планетарных волн:

-А

г

*>,=-'

к1+-

ёК

к 1

ГДе кк = {к2 + &22)2 ~ волновое число для горизонтального разреза (к2 =0); Нп - эквивалентная глубина; / = /0 + ¡}у - параметр Кориолиса,

/0 = (<р); Р - Юсо; £2 = 7,29 • 10~5 с-1 - частота вращения Земли; (р - широта; Л « 6300 км - радиус Земли; g « 10 м/с2 - ускорение свободного падения.

Для случая постоянной стратификации (который мы будем использовать в дальнейшем для наших оценок) эквивалентные глубины определяются как

и ^

к=—гт>

К

где N2 =- — — - частота плавучести; к - типичная глубина океана (для Ро &

нашего случая А = 5 ООО м); п - номер моды.

Тогда период синоптического возмущения вычисляется по формуле

2л

Т = — .

Фазовая скорость определяется из следующего соотношения:

С„ =-

К

ёК

Поскольку рассчитанная для первой моды волн Россби фазовая скорость в рассматриваемой области Северной Атлантики (от 7,5 до 48° с.ш.) много меньше скорости выполнения разреза (см. ниже), несинхронностью наблюдений можно пренебречь. Остается погрешность, обусловленная недостаточным разрешением сетки гидрологических станций на разрезах. Если расстояние между станциями много меньше длин основных мод планетарных волн, то соответствующая погрешность невелика. Она может быть оценена исходя из следующих соображений. Зная длину разреза и основную длину волны синоптического возмущения, можно определить, сколько целых волн синоптического масштаба укладывается на разрезе. Остаток синоптического возмущения и будет вносить искомую погрешность, которую можно оценить как

Ега-а/

.

где (). - остаточные единичные переносы тепла на каждой станции после

осреднения данных на разрезе по длине волны синоптического масштаба; ЛГ-число станций.

Тогда погрешность, связанная с межгодовой изменчивостью, рассчитывается как разность между общей погрешностью и суммой инструментальной, сезонной и синоптической погрешностей:

^межгод = ^общ ~

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

Характеристика использованных данных. Для расчета МПТ в Северной Атлантике использовались данные, полученные в 'период проведения эксперимента 1¥ОСЕ за 1992 - 1998 гг. в окрестности стандартных разрезов ЛОЗ (вдоль 36° с.ш.), А05 и АЯ01А (вдоль 24,5° с.ш.), А06 (вдоль 7,5° с.ш.) со средним расстоянием между станциями порядка 40 миль. Распределение данных во времени и количество станций на каждом разрезе приведены в табл. 1. Зондирования выполнялись С7Х>-комплексами от поверхности до дна океана. Данные о температуре и солености интерполировались линейно на равномерную по вертикали сетку с шагом 10 м до глубины 2 ООО м.

Таблица!

Распределение данных во времени и количество станций на каждом _разрезе ТУОСЕ__

Разрез

с.ш.

Период выполнения

Количество станций

лоз 36 09.1993 - 10.1993 133

А05 24,5 07.1992-08.1992 • 113

ARQIA 24,5 01.1998-02.1998 131

А06 7,5 02.1993 -03.1993 83

Данные о поле ветра, полученные с помощью скаттерометров со спутников ERS-1,-2 в период выполнения гидрологических наблюдений на разрезах, использовались для оценок дрейфовых переносов тепла. В стандартном информационном продукте WOCE [21] приведены с недельным разрешением касательные напряжения, рассчитанные по этому полю ветра. По этим напряжениям оценивались полные потоки в экмановском слое для каждой станции, которые затем интегрировались вдоль разреза методом трапеции.

Для анализа тенденции изменения МПТ в Северной Атлантике использовались все доступные прямые оценки МПТ, выполненные по данным наблюдений 1957 - 1998 гг. и опубликованные в работах [11, 20, 22 - 27] .

Результаты, их анализ и обсуждение. На рис. 1 показаны квазистационарные МПТ в Северной Атлантике через 6 кругов широты (7,5; 24,5; 26,5; 32; 36; 48° с.ш.), полученные путем осреднения всех прямых оценок, приведенных в табл. 2, и их суммарные погрешности ± а^. Подчеркнем, что

большинство прямых оценок МПТ было получено в Субтропической Атлантике между 24 и 36° с.ш. (табл. 2). Из рис. 1 видно, что максимум МПТ приходится на 26,5 - 32° с.ш. Меридиональный перенос тепла на этих широтах достигает 1,38 ПВт и направлен на север. Этот результат получен в работах [16, 27] путем осреднения большого количества данных с 1960 по 1992 гг., что дало возможность надежно выделить сезонный ход МПТ на этих широтах (рис.2, б, в).

И

2.0

1.2

0.8 0.4 0.0

10°

1 - 4 /

1 Ф

- Л.

20е

30°

40°

50°

с.ш.

Рис.1. Средний МПТ в Северной Атлантике по данным прямых оценок, проанализированных в работе. МПТ на 48° с.ш. обозначен перечеркнутым черным квадратом; на 36° с.ш. - ромбами (зачерненный ромб - оценка, полученная в настоящей работе, незачерненный - по [20]); на 32° с.ш. - незачерненным квадратом (по [27]); на 26,5° с.ш. - зачерненным кружком (по [26]) и незачерненным треугольником (по [16]); на 24,5° с.ш. - зачерненным квадратом (по нашим расчетам); на 7,5° с.ш. - зачерненным треугольником (по нашим расчетам). Вертикальные тонкие линии обозначают доверительные интервалы (±

Меридиональный перенос тепла на 24,5° с.ш. и средний МПТ на 26,5° с.ш. близки между собой и к МПТ на 32° с.ш. (рис. 1). Он составляет 1,27 ± 0,18 (на 24,5° с.ш.) и 1,33 ± 0,33 ПВт (на 26,5° с.ш.; средняя между двумя, приведенными на рис. 1, оценка). Представленная оценка для 24,5° с.ш. существенно превышает оценку, рассчитанную в настоящей работе по данным повторного разреза АЯ0Ы, выполненного зимой 1998 г. (0,85 ПВт). Максимальные величины МПТ на 24,5° с.ш. (1,38; 1,48 и 1,54 ПВт) были получены в работе [24] по данным 1957, 1981 и 1992 гг. соответственно. Оценки МПТ, сделанные в настоящей работе по данным стандартного разреза Л05 (июль - август 1992 г.), а также в работах [11, 18, 22, 25], указывают на то, что МПТ подвержен интенсивной сезонной изменчивости (см. рис. 2 и обсуждение этого вопроса ниже).

На 36° с.ш. осредненная по всем данным величина МПТ равна 1,02 ± 0,33 ПВт (рис. 1). Вместе с тем по оценкам, полученным в настоящей работе с использованием только данных №ОСЕ за сентябрь - октябрь 1993 г.,

МПТ достигает здесь 1,3 ПВт. По данным Колтермана и др. [24] МПТ на 36° с.ш. изменялся в различные годы и сезоны от 0,47 до 1,29 ПВт (табл. 2). Подробный анализ МПТ на 36° с.ш. был проведен ранее в работе Сато, Росс-би [20]. Авторы проанализировали данные нескольких трансатлантических разрезов, а также около 400 разрезов, выполненных в западной части Субтропической Атлантики. В результате ими был получен средний МПТ, направленный на север и равный 1,2 ± 0,3 ПВт. Оценка Тэлли [18] близка к оценкам, сделанным Роеммичем и Вуншем [17]. По их расчетам МПТ составляет на 36° с.ш. от 0,50 до 0,86 ПВт. Заметим, что большие из этих оценок получены по данным разрезов, выполненных в 1981 г., а меньшие - по данным за 1957 г., причем без учета сезонных вариаций.

Таблица2

Меридиональный перенос тепла (Н, ПВт) по-прямым оценкам в _Северной Атлантике

с.ш.

Авторы, годы выполнения разрезов

Н, ПВт

48

48

36

36

36

36

36

32

26,5

26,5

24

24

24

24,5

24,5

24,5

7,5

Колтерман и др. [24]

Колтерман и др. [24]

Роеммич, Вунш [17]

Наши оценки

Сато, Россби [20]

Paro, Россби [27]

Молинари и др. [26]

Филленбаум и др. [16]

Тэлли [18]

Лэвин и др. [25]

Колтерман и др. [24]

Холл, Брайден [11]

Роеммич, Вунш [17]

Наши оценки

Тэлли [18]

Фридрихе, Холл [23]

Наши оценки

0,27; 0,62; 0,53 0,62 0,86 0,47; 1,29; 0,7 0,5; 0,8

1,3 1,2 ±0,3 1,38 ±0,19 1,21 ±0,34 1,44 ±0,32 1,28 1,33; 1,2 1,38; 1,48; 1,54 1,22 1,2; 1,2 1,25; 0,85 0,73 0,3 0,79

Расчет характеристик осредненного сезонного цикла МПТ в Субтропической Атлантике выполнен в нескольких работах. По данным Paro, Россби [27] и Филленбаума и др. [16], МПТ в Субтропической Атлантике достигает своего максимума в летний период, а минимума - зимой. Типичная амплитуда годового хода составляет 0,4 0,5 ПВт (рис.2, б, в). Поэтому квазистацио-

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

Чтобы оценить амплитуду сезонной изменчивости МПТ в окрестности 36° с.ш., мы осреднили все оценки МПТ, полученные с использованием только прямых методов, по сезонам. После фильтрации сезонной изменчивости мы получили средний МПТ на 36° с.ш., приблизительно равный 1,07 ± 0,27 ПВт, и годовой ход МПТ, амплитуда которого составляет около 0,3 ПВт. Максимум МПТ наблюдается здесь в конце лета, а минимум - весной. Отметим, что зимние месяцы плохо обеспечены наблюдениями, что приводит, в частности, к завышению амплитуды полугодовой гармоники при выделении среднего сезонного хода методом наименьших квадратов (рис. 2, а). Неравномерность распределения данных по сезонам приводит также к тому, что их формальное осреднение вносит дополнительную погрешность из-за плохого разрешения сезонного хода. Из полученных результатов ясно, что эта погрешность составляет около 5 % среднегодовой величины МПТ. Действительно, при среднеарифметическом осреднении всех оценок среднегодовая величина МПТ оказалась равной около 1,02 ПВт, а после фильтрации сезонного хода МПТ изменился приблизительно на 0,05 ПВт. Аналогичная (возможно, несколько большая из-за худшей обеспеченности сезонных средних и близкого по амплитуде сезонного хода МПТ на анализируемых широтах) погрешность отмечается и по данным в окрестности 24,5° с.ш. Если осреднить все оценки МПТ вначале по сезонам, а затем за год, то среднегодовой МПТ здесь составит около 1,23 ± 0,14 ПВт.

Для зональных разрезов по 7,5 и 48° с.ш. выполнено значительно меньше прямых оценок МПТ, чем для Субтропической Атлантики (табл. 2). В настоящей работе по данным стандартного разреза И^ОСЕА06 МПТ на 7,5° с.ш. оценен в 0,79 ПВт. Это соответствует оценке, опубликованной в работе [18]. Вместе с тем, по данным Фридрихса и Холл [23], МПТ на 8° с.ш. равен 0,3 ПВт. В результате средний МПТ оценивается в 0,6 ПВт, а суммарная погрешность оценки МПТ на 7,5° с.ш. достигает 0,27 ПВт .(рис. 1). На 48° с.ш. осредненный МПТ оценивается в 0,51 ПВт.

Таким образом, можно заключить, что по данным прямых расчетов МПТ достигает максимума в Субтропической Атлантике. По данным балансовых расчетов квазистационарный МПТ достигает максимума в Тропической Атлантике. Однако различие МПТ в тропиках и субтропиках по данным этих расчетов не является статистически значимым [2, 28]. На основании проведенного анализа прямых оценок можно заключить, что квазистационарный МПТ характеризуется максимумом именно в Субтропической Атлантике.

Оцененная в настоящей работе погрешность МПТ, связанная с синоптическими флюктуациями, оказалась меньше, чем проанализированная выше погрешность из-за неадекватности описания сезонного цикла. Применение описанной выше методики, позволило получить характерные длины волн синоптических возмущений на разрезах. Они вычислялись по значимым пикам в спектрах волновых чисел и оказались равными для разных разрезов от 633 до 739 км. Типичные фазовые скорости, оцененные по формуле (8) с использованием дисперсионного соотношения (5), изменяются в пределах 1,6 - 11,4 см/с. Чем ближе к экватору, тем скорость синоптического возму-

щения выше (табл. 3). Однако даже для самого южного разреза эта скорость на порядок меньше средней скорости судна при выполнении гидрологических работ, которую можно оценить по данным, приведенным в табл. 1. Поэтому несинхронностью выполнения разрезов можно пренебречь, на что указано в предыдущем разделе. Поскольку на разрезе, как правило, не укладывается целое число длин волн синоптического масштаба, то в оценку МПТ вносится некоторая дополнительная погрешность (асин)- Она оценивалась по описанной выше методике и составила 0,04; 0,03 и 0,01 ПВт на 36; 24,5 и 7,5° с.ш. соответственно. Полученные результаты указывают на то, что погрешности оценки интегральных МПТ за счет синоптических шумов малы и составляют 1 - 4%.

2.0 -т------

1.6-----а---

д

н 12____^А ¿^-^^А__|

« 0.8---—Н1--А-4*"--

0.4----Л----

0.0 -I------

а

2.0 1.6

д 1.2

И

И 0.8 0.4 0.0

Месяцы

Р и с.2. Внутригодовая изменчивость МПТ: а - на 36° с.ш. по табл.2, б - на 32° с.ш. по [27], в - на 26,5° с.ш. по [16]

А

Л £1 д Л ^- -А д

N Л А

д А £ д

Таким образом, типичная суммарная погрешность оценки среднегодового МПТ в Субтропической Атлантике по данным выполненных зональных разрезов за счет неполной фильтрации сезонной и синоптической изменчивости не превышает 10 %.

ТаблицаЗ

Характеристики синоптических возмущений на зональных разрезах_

Длина волны, км Амплитуда, ПВт Фазовая скорость, см/с

36 739 0,86 1,6

24,5 633 0,55 2,7

7,5 738 0,11 11,4

Примечание. Длина волны выделена по пикам в спектрах волновых чисел, значимых на уровне 95 %. Фазовая скорость определялась из дисперсионного соотношения (8).

Принимая во внимание величины приведенных МПТ и их суммарных погрешностей можно сделать вывод, что разница между полученными оценками стационарного МПТ через 24,5; 26,5; 32 и 36° с.ш. статистически незначима. Другими словами, при существующей наблюдательной базе и погрешностях оценок МПТ невозможно достоверно определить широту в Субтропической Атлантике, на которую приходится максимум МПТ.

Межгодовая изменчивость МПТ, оцененная для 36° с.ш. по формуле (10) с учетом полученных выше результатов, привюдит к среднеквадратичной погрешности оценки, равной ~ 0,27 ПВт. Рассмотрим долговременную тенденцию изменения МПТ на этой широте, после фильтрации которой <7„ежгод

уменьшается на 25%. На рис. 3, а можно отметить тенденцию роста МПТ в окрестности 36° с.ш. за период с конца 50-х до середины 90-х гг. Видно, что с 1959 по 1993 гг. МПТ увеличился от 0,6 до 1,1 ПВт. Полученный линейный тренд значим на 95%-ном доверительном уровне. Оценки, полученные Коль-терманом и др. на этой широте, также свидетельствуют об увеличении МПТ за период с 1959 по 1981 гг. Причем, по их данным, это увеличение составило 0,82 ПВт (от 0,47 ПВт до 1,29 ПВт [24]). Однако существенная часть этого тренда обусловлена тем обстоятельством, что в 50-е гг. наблюдения в основном выполнялись весной, а в 80 - 90-е - летом. После фильтрации сезонной изменчивости с использованием данных об осреднённом сезонном цикле МПТ (см. рис. 2, а) тренд МПТ на 36° с.ш. уменьшился. Из сравнения рис. 3, а и 3, б хорошо видно, что в действительности увеличение МПТ между 1959 и 1993 гг. составило 0,35 ПВт (от 0,75 до 1,1 ПВт). Подчеркнем вместе с тем, что тренд остался значимым на 95%-ном доверительном уровне. Таким образом, можно сделать вывод об интенсификации МПТ в Северной Атлантике во второй половине XX в. Это свидетельствует о том, что океан активно участвует в процессе глобального изменения климата, уменьшая термические контрасты между низкими и высокими широтами, поглощая часть ISSN 0233-7584. Мор. гидрофиз. журн., 2007, Nb 1 55

тепла и ослабляя, таким образом, парниковый эффект. Вполне вероятно, что обнаруженная тенденция, по крайней мере частично, представляет собой проявление низкочастотных вариаций климатической системы, обусловленных изменчивостью интенсивности термохалинной циркуляционной ячейки с типичным периодом, превышающим 50 лет. Эти вариации не могут быть надежно описаны с помощью прямых оценок МПТ из-за их низкого разрешения [1,29].

1.6

1950 1960 1970 1980 1990 2000

Годы 6

Р и с.З. Межгодовая изменчивость МПТ в Северной Атлантике на 36° с.ш. без фильтрации сезонной изменчивости

Заключение. В работе проведен обобщенный анализ оценок МПТ, полученных прямыми методами различными авторами, а также рассчитанных нами по данным стандартных разрезов WOCE в окрестности 7,5; 24,5; 26,5; 32; 36 и 48° с.ш. Анализ показал, что максимального значения суммарный МПТ достигает в Субтропической Атлантике. Его величина на 32° с.ш. оценивается в 1,38 ± 0,19 ПВт. Средний МПТ на 7,5; 24,5; 26,5; 36 и 48° с.ш. составляет соответственно 0,6 ± 0,27; 1,23 ± 0,14; 1,33 ± 0,33; 1,07 ± 0,27 и 0,51 ± 0,17 ПВт. Меридиональный перенос тепла в Субтропической Атлантике подвержен достаточно интенсивной сезонной изменчивости. В конце лета наблюдается максимум МПТ (около 1,9 ПВт), весной - минимум (около 0,8 ПВт). Анализ погрешностей разного типа позволил сделать вывод, что в настоящее время невозможно достоверно определить широту в Субтропической Атлантике, на которую приходится максимум квазистационарного 56 ISSN 0233-7584. Мор. гидрофиз. журн., 2007, № 1

МПТ. Выявлена значимая тенденция увеличения МПТ в Субтропической Атлантике во второй половине XX в. В окрестности 36° с.ш. МПТ увеличился за 1959 - 1993 гг. от 0,75 до 1,1 ПВт, что свидетельствует об интенсификации меридиональной океанической циркуляции.

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

СПИСОК ЛИТЕРАТУРЫ

1. Полонский А.Б. Роль океана в современных изменениях климата // Морской гидрофизический журнал. - 2001. - № 6. - С. 28 - 49.

2. Тимофеев Н.А., Юровский А.В. Радиационные, тепло- и водобалансовые режимы океанов. Климат и изменчивость. - Севастополь: ЭКОСИ-Гидрофизика, 2004. - 256 с.

3. Biyan К. Poleward heat transport by the oceans: observations and models // Ann. Rev. Rarth Planet. Sci. - 1982. - 10, № 1. - P. 15 - 38.

4. Hastenrath S. Heat budget of ocean // J. Phys. Oceanogr. - 1980. - № 2. - P. 159 - 170.

5. Hsuing J., Newell R.E., Houghtby T. The annual cycle of oceanic heat storage and oceanic meridional heat transport // Quart. J. Royal Met. Soc. - 1989. - P. A. - П5, № 485. - P. 1 -28.

6. Lamb P.J. Estimate of annual variation of Atlantic Ocean heat transport // Nature. - 1981. -290, № 4. - P. 766-768.

7. Oort A.H., Haar Т.Н. V. On the observed annual cycle in the ocean - atmosphere heat balance over the northern hemisphere // J. Phys. Oceanogr. - 1976. - 6, № 6. - P. 781 - 800.

8. Trenberth K.E., Caron J.M., Sepaniak D.P. The atmospheric energy budget and implications for surfaces fluxes and ocean heat transports // Clim. Dyn. - 2001. - 17, № 5. - P. 259 - 276.

9. Стоммел Г. Обзор теории морских течений // Проблемы океанической циркуляции. - М.: Мир, 1965.-С. 5-14.

10. Полонский А.Б. Циркуляция вод Тропической Атлантики и меридиональный перенос тепла // Морской гидрофизический журнал. - 1985. - № 1. - С. 58 - 62.

11. Hall М.М., Bryden H.L. Direct estimates and mechanisms of ocean heat transport // Deep-Sea Res. - 1982. - 29, № 3 A. - P. 339 - 359.

12. Roemmich D. Estimation of meridional heat flux in the North Atlantic by inverse methods // J. Phys. Oceanogr. - 1980. - JjO, № 12. - P. 1972 - 1983.

13. Wunsch C. The North Atlantic circulation west of 50° W determined by inverse methods // Rev. Geophys. Space Phys. - 1978. - 16.- P. 583 - 620.

14. Полонский А.Б. Меридиональные переносы тепла в Тропической Атлантике и механизмы их формирования / Дис. ... д-ра геогр. наук. - Севастополь: МГИ АН УССР, 1990. -С. 320.

15. Stammer D., Wunsch С., Giering R. et al. Volume, heat, and freshwater transports of the global ocean circulation in 1993 - 2000, estimated from a general circulation model constrained by World Ocean Circulation Experiment (WOCE) data // J. Geophys. Res. - 2003. - 108, № CI.-3007 doi.: 10.1029 / 2001JC001115.

16. Fillenbaum E.R, Lee T.N., Jouhns W.E. et al. Meridional heat transport variability at 26.5°N in the North Atlantic // J. Phys. Oceanogr. - 1997. - 27, № 1. - P. 153 - 174.

17. Roemmich D., Wunsch C. Two transatlantic sections: meridional circulation and heat flux in the subtropical North Atlantic Ocean // Deep-Sea Res. - 1985. - 32* № 6 A. - P. 619 - 664.

18. Talley L.D. Shallow, intermediate, and deep overturning components of the global heat budget // J. Phys. Oceanogr. - 2003.-33, № 3. - P. 530 - 560.

19. Полонский А.Б., Крашенинникова С.Б. Оценка океанического меридионального переноса тепла в Субтропической Атлантике // Морской гидрофизический журнал. - 2006. - № 1. -С. 3-15.

20. Sato О.Т., Rossby Т.Н. Seasonal and low-frequency variability of the meridional heat flux at 36° N in the North Atlantic // J. Phys. Oceanogr. - 2000. - 30, № 3. - P. 606 - 621.

21. WOCE Global Data. Version 3.0 2002. WOCE International Project Office // WOCE Report № 180 / 02 Southampton, UK.

22. Bryden H.L., Hall M.M. Heat transport by currents across 25°N latitude in the Atlantic Ocean // Science. - 1980. - 207, № 4433. - P. 884 - 886.

23. Friedrichs M.A., Hall M.M. Deep circulation in the tropical North Atlantic // J. Mar. Res. -1993. - 5JL - P. 697 - 736.

24. Koltermann K.P., Sokov A.V., Tereschenkov V.P. et al. Decadal changes in the thermohaline circulation of the North Atlantic // Deep-Sea Res. - 1999. - 46, pt 2. - P. 109 - 138.

25. Lavin A., Bryden H.L., Parrilla G. Meridional transport and heat flux variations in the subtropical North Atlantic I I Global Atmos. Ocean. Syst. - 1998. - 6, - P. 269 - 293.

26. Molinari R.L., Johns E., Eesta J.F. The annual cycle of meridional heat flux in the Atlantic Ocean at 26.5°N // J. Phys. Oceanogr. - 1990. - 20, № 3. - P. 476 - 482.

27. Rago T.A., Rossby Т.Н. Heat transport into the North Atlantic ocean north of 32°N latitude // J. Phys. Oceanogr. - 1987. - Г7, № 7. - P. 854 - 871.

28. Trenberth K.E., Caron J.M. Estimates of Meridional Atmosphere and Ocean Heat Transports I I J. Climate. - 2001. - 14, № 16. - P. 3433 - 3443.

29. Born K. Changes in baroclinicity and synoptic activity on the northern hemisphere from the 1960s to the present as seen by distinct data sets // Ann. Geophys. - 1998. - 16. № 3. -P. 318.

Морской гидрофизический институт HAH Украины. Материал поступил

Севастополь в редакцию 05.09.05

После доработки 07.10.05

ABSTRACT WOCE (World Ocean Circulation Experiment) performed in the North Atlantic in 1992 - 1998 are used to calculate the meridian heat transport (MHT) in the ocean. The generalized analysis of MHT estimates obtained by direct methods using the data of the sections carried out between 7.5° and 36°N in the second part of the 20th century is done. MHT values averaged for the whole period of observations achieve their maximum (1.38 ± 0.19 PW) in the Subtropical Atlantic. Season-to season MHT variability is rather intensive. Maximum MHT (about 1.9 PW) is observed in the Subtropical Atlantic in the end of summer, and minimum MHT (about 0.8 PW) - in the end of winter. Significant MHT increase (from 0.75 to 1.1 PW) at 36°N in 1959 - 1993 is revealed. It testifies to the intensification of the meridian oceanic circulation in the North Atlantic.