Бесплатный автореферат и диссертация по наукам о земле на тему
Статистический анализ межгодовой изменчивости теплосодержания и стерических колебаний уровня Мирового океана
ВАК РФ 25.00.28, Океанология

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

Министерство образования и науки Российской Федерации ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ

УНИВЕРСИТЕТ"

СТАТИСТИЧЕСКИЙ АНАЛИЗ МЕЖГОДОВОЙ ИЗМЕНЧИВОСТИ ТЕПЛОСОДЕРЖАНИЯ И СТЕРИЧЕСКИХ КОЛЕБАНИЙ УРОВНЯ МИРОВОГО ОКЕАНА

На правах рукописи

Глок Наталья Ивановна

Специальность 25.00.28 - океанология

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата географических наук

2 8 НОЯ 2013

РГГМЫ

Санкт-Петербург

2013

005540500

005540500

Министерство образования и науки Российской Федерации ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ

УНИВЕРСИТЕТ"

СТАТИСТИЧЕСКИЙ АНАЛИЗ МЕЖГОДОВОЙ ИЗМЕНЧИВОСТИ ТЕПЛОСОДЕРЖАНИЯ И СТЕРИЧЕСКИХ КОЛЕБАНИЙ УРОВНЯ МИРОВОГО ОКЕАНА

На правах рукописи УДК 551.465

Глок Наталья Ивановна

Специальность 25.00.28 - океанология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата географических наук

рггма

Санкт-Петербург 2013

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Российский государственный гидрометеорологический университет (РГГМУ)

Научный руководитель:

доктор географических наук, профессор

Официальные оппоненты:

доктор географических наук, профессор

кандидат географических наук, старший научный сотрудник

Малинин Валерий Николаевич

Алексеев Генрих Васильевич Белоненко Татьяна Васильевна

Ведущая организация: Санкт-Петербургское отделение

федерального государственного бюджетного учреждения «Государственный океанографический институт имени Н.Н.Зубова» (СПО ФГБУ «ГОИН»)

Защита состоится "12" декабря 2013 г. в 15 час. 30 мин. на заседании диссертационного совета Д 212.197.02 в Российском государственном гидрометеорологическом университете по адресу: 195196, Санкт-Петербург, Малоохтинский проспект, 98, ауд. 308.

С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета по адресу: Санкт-Петербург, Малоохтинский проспект, 98

Автореферат разослан " 12 " ноября 2013г.

Ученый секретарь диссертационного совета, кандидат географических наук

Воробьев Владимир Николаевич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Климат на нашей планете постоянно меняется, и, в первую очередь, это вызвано резким изменением глобальной температуры воздуха, которая в течение XX века увеличилась на 0,6-0,7 °С. Учитывая тесную взаимосвязь атмосферы и океана, процессы, происходящие в ней, не могут не сказываться на водной оболочке нашей планеты. Одновременно процессы в океане оказывают стимулирующее влияние на тепловой режим и общую циркуляцию атмосферы. В этом плане особую роль играет теплосодержание Мирового океана (ТСМО), которое является важнейшим фактором изменений глобального климата, так как он имеет высокую теплоемкость и инерционность. Кроме того, ТСМО практически полностью определяет стерические колебания уровня Мирового океана (СКУМО), вклад которых в формирование тренда глобального уровня может составлять от 15 до 50 % в зависимости от интервала временного осреднения.

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

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

Исходя из поставленной цели, были сформулированы следующие задачи:

- выявление закономерностей межгодовой изменчивости теплосодержания Мирового океана;

- построение статистических моделей межгодовой изменчивости теплосодержания Мирового океана по данным спутников и реанализа о температуре поверхностного слоя воды;

- выявление закономерностей термостерических колебаний уровня Мирового океана (УМО);

-физическое обоснование и построение физико-статистической модели термостерических колебаний уровня Мирового океана (ТСКУМО);

- построение статистической модели термостерических колебаний уровня Мирового океана по реанализу и спутниковым данным температуры его поверхности.

Материалы и методы. Основой работы послужили следующие архивы гидрометеорологических данных, находящихся в свободном доступе в сети Интернет:

-среднемесячные данные о поверхностной температуре океана в узлах широтно-долготной сетки 2 х 2° из глобального архива NOAA NCDC ERSST (National Oceanic and Atmospheric Administration National Climatic Data Center Extended Reconstructed Sea Surface Temperature);

- суточные данные о поверхностной температуре океана в узлах широтно-долготной сетки 0,25 * 0,25° из глобального архива NOAA NCDC OISST (National Oceanic and Atmospheric Administration National Climatic Data Center Optimum Interpolation Sea Surface Temperature);

- среднемесячные данные о глубоководной температуре и вертикальных скоростях воды в узлах широтно-долготной сетки 0,5 х 0,5° из глобального архива SODA (Simple Ocean Data Assimilation);

- тренды поверхностной температуры океана из архива HadlSSTl (Hadley Central Sea Ice and Sea Surface Temperature Dataset) в узлах географической ши-ротно-долготной сетки 2 х 2°.

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

Научная новизна:

- выявлены закономерности межгодовой изменчивости теплосодержания Мирового океана. Показано, что за период с 1955 по 2010 гг. трендовая составляющая является значимой и носит положительный характер, однако, начиная с 2003 г. происходит значительное замедление роста теплосодержания. В данных ТСМО обнаружены два устойчивых цикла с периодами в 23 и 15 лет;

- установлены особенности распределения трендов временного хода температуры воды в отдельных океанах в слое 0-3000 м. Для всех океанов тренды являются значимыми до глубины порядка 1500 метров. При этом на разных горизонтах отмечаются как положительные, так и отрицательные тренды. Однако, начиная с 500 м, проявляется тенденция к уменьшению величины положительного тренда с глубиной;

-разработаны статистические модели межгодовых изменений теплосодержания Мирового океана по данным спутников и реанализа о температуре поверхностного слоя воды;

- выявлены закономерности термостерических колебаний уровня Мирового океана. Показано, что в термостерических колебаниях уровня Мирового океана с 1880 по2005 гг. присутствует значимое колебание с периодом в 57 лет, которое отчетливо проявляется и для температуры поверхности океана. Обнаружены также колебания с периодами 22 - 25 лет, 15 и 7 - 8 лет. Установлено, что трендовая компонента обеспечивает максимальный вклад в дисперсию тер-мостерического ряда, который составляет не менее 77 %, причем на долю значимых циклов приходится в среднем около 10 %;

- разработана физико-статистическая модель термостерических колебаний УМО, представляющая собой сумму двух компонент: «быстрой», отвечающей

за внешний тепловой баланс океана и «медленной», определяющей перераспределение тепла в толще воды до нижней глубины главного термоклина;

-разработана статистическая модель термостерических колебаний УМО на основе спутниковых данных о температуре поверхности океана (ТПО);

На защиту выносятся:

- закономерности изменений теплосодержания и термостерических колебаний уровня Мирового океана в современных климатических условиях;

- статистические модели межгодовых изменений теплосодержания океана по данным спутников и реанализа температуры поверхностного слоя воды;

- физико-статистическая модель термостерических колебаний уровня Мирового океана;

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

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

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

Апробация работы и публикации. Данная работа была выполнена в рамках Аналитической ведомственной целевой программы «Оценить последствия повышения уровня океана в условиях современных изменений климата» №2.1.1/6690 и Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы (государственные контракты № П726 и № 02.740.11.0849), а также в рамках гранта Правительства РФ (Договор № 11.034.31.0078) для поддержки исследований под руководством ведущих ученых. Основные результаты диссертации были представлены в рамках шестнадцатой Санкт-Петербургской ассамблеи молодых ученых и специалистов и отмечены грантом губернатора Санкт-Петербурга.

По теме диссертации опубликовано 8 работ, в том числе 6 из них в рецензируемых журналах по перечню ВАК.

Структура и объем работы. Диссертация состоит из введения, 4 глав, заключения и библиографического списка использованной литературы. Объем диссертации составляет 125 страницы, включая 13 таблиц и 33 рисунка. Список литературы составляет 100 наименований.

Благодарности. Автор искренне благодарит научного руководителя профессора В.Н. Малинина за конструктивное руководство и всестороннюю поддержку в работе, а так же к.г.н. С.М. Гордееву за консультации и помощь.

СОДЕРЖАНИЕ РАБОТЫ

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

Первая глава «Современное состояние изученности теплосодержания и стерических колебаний уровня Мирового океана» посвящена исследованию и методам расчета теплосодержания и стерических колебаний уровня Мирового океана.

В разделе 1.1 «Основные факторы формирования колебаний уровня Мирового океана» рассмотрены три основных группы факторов, влияющих на колебания УМО: космогеофизические силы, геолого-геодинамические и гидрометеорологические процессы. Оценены периоды и масштабы их действия.

На основе отчетов Межправительственной группы экспертов по изменению климата проанализированы количественные оценки вкладов отдельных факторов в формирование тренда УМО. Показано, что стерические колебания в отдельные периоды времени (1993 - 2003 гг.) могут определять до 50 % изменений уровня Мирового океана.

Раздел 1.2 «Методы получения исходных данных, используемых для оценки теплосодержания и стерических колебаний морского уровня» посвящен исследованию двух подходов определения теплосодержания и СКУМО - использованию глубоководных данных на основе контактных методов измерений с последующим их пересчетом и применению бесконтактных спутниковых измерений, как разности между альтиметрическими и GRACE (Gravity Recovery And Climate Experiment) данными. Рассмотрены всемирные базы данных World Ocean Database (WOD) и атласы World Ocean Atlas (WOA), на основе которых проводятся глобальные расчеты ТСМО и СКУМО, проанализированы технические устройства, применяемые для формирования баз WOD, их основные характеристики, принципы работы, способы измерений и погрешности, а так же районы их применения. Подробно рассмотрены проблемы, возникающие при использовании глубоководных данных, связанные с пространственным покрытием Мирового океана, а так же с точностью измерений. Для бесконтактного метода приведены принципы работы и точность альтиметрических и GRACE измерений.

Для временного интервала с 2004 по 2008 гг. проведен сравнительный анализ СКУМО по данным, полученным обоими способами. Следует отметить, что ряды сильно отличаются, оценки СКУМО по данным проекта Арго имеют заниженные величины относительно бесконтактных методов и несколько отличный характер изменения. Причина такого расхождения еще до конца не определена, но может быть связана с недостаточной выборкой in situ измерений.

В разделе 1.3 «Краткий обзор исследований теплосодержания Мирового океана» проведен анализ исследования ТСМО по литературным данным.

Существенная часть тренда ТСМО определяется слоем от поверхности до 700 м, который описывает 69 % от изменения теплосодержания в слое от 0 до 3000 м.

В работе представлены методы устранения систематических ошибок архивных данных для расчета достоверных оценок ТСМО. Подробно рассмотрен ряд версий температурного анализа [Ыш е1. а1., 2006], полученных по различным базам данных с применением коррекции ошибок, как для отдельных видов измерительных приборов, так и для их совокупности, расхождения между ними, и влияние их на оценки ТСМО.

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

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

В разделе 1.4 «Краткий обзор исследований стерических колебаний уровня Мирового океана» проведен анализ СКУМО по литературным данным. Рассмотрены различные подходы к оценке ТСКУМО.

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

Показано, что большая часть трендовой составляющей ТСКУМО формируется в слое от 0 до 700 м и описывает 75% тренда временного хода термосте-рического уровня для слоя от 0 до 3000 м, а коэффициент корреляции между рядами ТСКУМО данных слоев составляет г = 0,99. Таким образом, хотя тер-мостерические колебания в слое от 0 до 700 м заметно занижают величину тренда, они вполне адекватно описывают межгодовые ТСКУМО.

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

Вторая глава «.Статистические методы анализа данных, используемые в работе» посвящена анализу статистических методов обработки данных, их характеристикам, достоинствам и недостаткам.

В разделе 2.1 «Анализ трендовой компоненты временного ряда» рассмотрены особенности трендовой компоненты и ее изменчивость относительно выбранного временного интервала. Определены особенности применения параметрических и непараметрических характеристик тренда. Описаны основные параметрические характеристики: коэффициент детерминации, величина и значимость тренда.

Раздел 2.2 «Анализ циклических колебаний временного ряда» посвящен выявлению и оценке циклических (гармонических) колебаний. Рассмотрены ограничения применения гармонического и спектрального анализа при исследовании гидрометеорологических процессов. Обосновано применение вейвлет-анализа и описаны его основные принципы, позволяющие представлять исследуемый одномерный сигнал в двухмерной плоскости частота — положение. Рассмотрено прямое и обратное интегральное вейвлет-преобразование, дискретное вейвлет-преобразование и проведено их сравнение относительно решаемой задачи. Введены условия, выполнение которых необходимо для конструирования базисной анализирующей функции. Рассмотрены разновидности вейвлетобра-зующих функций, обоснован выбор вейвлета Морле и приведено его описание, получен вейвлет-спектр, полная энергия сигнала и глобальный спектр энергии.

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

В третьей главе «Теплосодержание Мирового океана и его межгодовая изменчивость» обсуждаются закономерности колебаний теплосодержания Мирового океана и результаты построения статистических моделей ТСМО.

В разделе 3.1 «Закономерности колебаний теплосодержания Мирового океана» рассмотрены оценки межгодовой изменчивости теплосодержания, рассчитанные различными группами ученых с применением разных методик обработки исходной информации. Проведена оценка расхождений трендов, величины которых за период с 1955 по 2009 гг. изменяются от 0,16 х 1022 до 0,26 х 1022Дж/год. Корреляция между временными рядами аномалий ТСМО достаточно высока и изменяется в пределах от 0,70 до 0,94, а после исключения трендов - в пределах 0,54-0,79. Для нивелирования случайных ошибок имеющихся рядов путем осреднения отдельных членов каждого из них в фиксиро-

ванные моменты времени сформирован осредненный ряд аномалий ТСМО. Ос-редненному ряду аномалий ТСМО свойствен отчетливый линейный тренд (7> = 0,22 х 1023 Дж/год), описывающий 75 % дисперсии ряда.

Для выявления внутренней структуры межгодовой изменчивости осред-ненного ряда аномалий ТСМО применен вейвлет-анализ, благодаря которому выявлены два значимых цикла с периодами 23 и 15 лет, первый из которых определен на всем временном интервале, а второй лишь с 1975 г. по 2005 г. (рисунок 1)

Годы

Рисунок 1 - Коэффициенты вейвлет-разложения осредненного по данным разных авторов временного ряда аномалий теплосодержания Мирового океана за период с 1955 по 2009 гг. на основе вейвлета Морле

В разделе 3.2 «Статистическая модель оценки среднегодовых значений теплосодержания Мирового океана на основе реанализа температуры его поверхности» продемонстрирован согласованный характер межгодовых колебаний аномалий ТСМО с глобальными значениями аномалий ТПО, полученной из архива Нас^ЭТ, и приповерхностной температуры воздуха архива С1ШТЕМ4. Доказана принципиальная возможность построения статистической модели ТСМО по архивным данным о температуре поверхности океана.

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

По данным реанализа аномалий ТПО архива ЕЯБЗТ методом множественной линейной регрессии построены статистические модели аномалий ТСМО при различном числе включаемых переменных, и проведен комплексный анализ их основных параметров. Статистическая модель для оценки аномалий ТСМО имеет вид:

AH=f(ATuAT2,...,ATm),

(1)

где ДТ,- - аномалии ТПО в узлах выбранной сетки, гп - число узлов (8200 значений).

В качестве исходной выборки ТСМО был выбран ряд, полученный [Ьеукш е1 а1., 2009], который рассчитан по максимальному объему глубоководной информации и имеет постоянную поддержку. Указанный ряд был разделен на 2 выборки: зависимую (1968-2004 гг.), по которой определялись параметры модели, и независимую (2005-2009 гг.), используемую для проверки ее качества.

При увеличении количества предикторов в модели коэффициент детерминации плавно возрастает, а стандартная ошибка также плавно убывает. Оценки р-1еуе1 коэффициентов регрессии показали, что в моделях с 8 и 10 предикторами несколько коэффициентов регрессии являются незначимыми, поэтому в качестве оптимальной была принята модель с 7 предикторами, для которой Я2 = 0,88, а ау = 0,6 х 1022 Дж.

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

12 п

Годы

Рисунок 2 - Временной ход аномалий теплосодержания Мирового океана по данным [Levitus et al., 2009]. Толстая кривая - исходные значения аномалий ТСМО, тонкая кривая - вычисленные значения аномалий ТСМО по модели (1). Квадрат - начало расчета аномалий ТСМО по независимым данным о ТПО.

В разделе 3.3 «Статистическая модель оценки теплосодержания Мирового океана на основе спутниковых данных температуры поверхностного слоя воды» рассмотрен архив спутниковых данных ТПО Optimum Interpolation Sea Surface Temperature (OISST), его достоинства и недостатки. Произведено сравнение глобальных значений аномалий ТПО архивов OISST и ERSST, показано

10

наличие высокой корреляции (Я = 0,94) между рядами и небольшое систематическое занижение ТПО по данным архива 0188Т.

Оценка связи аномалий ТСМО с глобальными рядами ТПО по данным архива 0188Т за период с 1982 по 2009 гг. позволила определить, что максимальная корреляция отмечается при запаздывании теплосодержания на 1 год (Я = 0,90). Это позволяет построить простую линейную прогностическую модель теплосодержания Мирового океана с заблаговременностью в 1 год:

Н,= Ъ0 + Ь,ТиЛ), (2)

Среднеквадратическая ошибка данной модели составляет 0У= 1,62 х 1022 Дж, что меньше стандартного отклонения исходной выборки.

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

Оптимальная модель теплосодержания содержит 7 предикторов, расположенных в разных частях океана, имеет коэффициент детерминации равный Я2 = 0,94 и стандартную ошибку оу = 0,9 х 1022Дж, что значительно меньше стандартного отклонения исходного временного ряда аномалий ТСМО (рисунок 3).

о

и н

Годы

Рисунок 3 - Сопоставление фактических значений аномалий теплосодержания Мирового океана (линия 1) и вычисленных по моделям (2) и (1) (линии 2 и 4, соответственно). Пунктиром показаны независимые расчеты по моделям (2) и (1) (линии 3 и 5, соответственно); маркером (6) отмечено фактическое независимое значение теплосодержания.

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

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

Четвертая глава посвящена разработке «физико-статистической модели термостерического уровня Мирового океана».

В разделе 4.1 «Закономерности термостерических колебаний уровня Мирового океана» произведено сравнение четырех рядов данных временной изменчивости ТСКУМО, полученных различными учеными. Для периода с 1960 по 2003 гг. коэффициенты корреляции между исходными рядами всех авторов являются значимыми по критерию Стьюдента и имеют корреляцию в пределах 0,71 - 0,90. Для стационарных рядов (после исключения тренда), она уменьшилась до 0,46 - 0,64.

Однако, несмотря на высокую корреляцию между рядами ТСКУМО, оценки линейных трендов существенно различаются. Максимальный тренд ТСКУМО (7> = 0,45 мм/год) отмечается в оценках [Маг'сеЦа, 2010], а минимальный проявляется в данных [Ьеукиэ е! а1., 2012] и составляет Тг = 0,32 мм/год.

Проведено вейвлет-преобразование рядов ТСКУМО, полученных различными авторами. Для ряда [Маг'сеЦа, 2010] за период с 1880 по 2005 гг., отчетливо выделяется 57-летний цикл, определенный на всем временном интервале. Цикл с периодом 19-20 лет выражен значительно слабее и присутствует в данных до середины 80-х годов. Кроме того можно отметить слабые 6-8-летние колебания, наблюдаемые с середины прошлого столетия.

На более коротких рядах ТСКУМО, рассчитанных [Ьеу^Шв е! а]., 2012] и [Вогтг^иев е1 а1., 2008], выявлены постоянные циклы с периодом 23 - 25 лет, а так же два локальных цикла: первый, с периодом 15 лет, определенный с 1970 по 1975 гг. до начала 21 века, и второй с периодом 6-8 лет, определенный на всем временном интервале в данных [Ботш^ез е! а]., 2008], а в данных [Ьеу^Б е1 а1„ 2012] с 1992 по 2010 гг.

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

В разделе 4.2 дается «физическое обоснование модели». Обсуждается гипотеза, в соответствии с которой формирование термостерических колебаний уровня в конкретный момент времени представляет сумму двух компонент: «быстрой» и «медленной». «Быстрая» компонента обусловлена внешним теп-

ловым балансом океана и определяет колебания теплосодержания деятельного слоя (с периодом тепловой адаптации системы океан-атмосфера в несколько недель), а «медленная» компонента - нагревание (охлаждение) океана сверху и перераспределение тепла по всей его толще до нижней границы главного термоклина, что в основном определяется трендом.

Для широтных зон 15,5 - 20° с.ш. и 32,5 - 37° с.ш., где наблюдаются пики максимума и минимума изменения термостерического уровня для отдельных океанов и Мирового океана в целом, дана оценка вертикального распределения линейных трендов временных рядов температуры воды. Распределение трендов по глубине в разных океанах существенно отличается, что свидетельствует о неодинаковых причинах их формирования, при этом значимые тренды температуры наблюдаются до глубин порядка 1500 м. Решающее влияние на формирование трендов оказывает система океанской циркуляции, которая «размазывает» поступающий сверху «климатический сигнал» практически по всей акватории океана. Оценки вертикальных скоростей на основе глобального архива SODA показали, что климатический сигнал от поверхности до глубин 1500 - 1800 м может дойти примерно за 20 лет.

Раздел 4.3 «Модель «быстрой» компоненты термостерических колебаний уровня Мирового океана» посвящен рассмотрению случайных колебаний ТСКУМО.

Статистическую модель «быстрой» компоненты ТСКУМО можно представить в виде:

/гвыси=ДАГь ДГ2,..., ДГ„), (3)

где/- некая функция от аномалий годовых значений ТПО, т - число исходных точек с данными ТПО, взятых из глобального архива NOAA NCDC ERSST v3 в узлах широтно-долготной сетки 2 х 2°.

Для расчета «быстрой» компоненты ряда ТСКУМО все данные поверхностной температуры были приведены к стационарному виду путем удаления трендовой составляющей с целью предотвратить ее возможное влияние на формирование тренда в «быстрой» компоненте. Из 8280 значений трендов поля ТПО 6500 (78 %) оказались значимыми по критерию Стьюдента. Для уменьшения количества используемых при расчете архивных данных были выбраны значимые коэффициенты корреляции «быстрой» компоненты термостерического уровня океана с ТПО, где в качестве «быстрой» компоненты ТСКУМО выступала разность исходных значений ряда ТСКУМО и его трендовой составляющей. В результате получено 2590 значимых коэффициентов корреляции.

Для построения модели быстрой компоненты ТСКУМО применялся метод множественной линейной регрессии с алгоритмом пошагового включения переменных. Модель рассчитывалась для зависимой выборки за период с 1955 по 2007 гг. Результаты расчетов показали, что до 9 шага все включаемые в модель переменные являются значимыми при уровне значимости а = 0,05. К сожалению, выделить оптимальный шаг модели оказалось чрезвычайно сложно, по-

этому исходя из распределения статистических параметров, рассматривались модели в которые вошли от 1 до 7 предикторов. Исключив первые 2 шага которые дают максимальную ошибку, остальные 5 были усреднены, и таким образом была получена средняя модель по 5 моделям. На зависимом периоде модель описывает 83 % дисперсии исходного ряда быстрой компоненты ТСКУМО, при этом стандартная ошибка составляет 1,28 мм. Стандартная ошибка модели для независимых данных равна 2,86 мм, что меньше стандартного отклонения исходного временного ряда ТСКУМО (ст = 3,11 мм).

В разделе 4.4 «Модель «медленной» компоненты термостерических колебаний уровня Мирового океана» отдельно рассмотрена «медленная» компонента ТСКУМО, которая была представлена как сумма тренда и колебаний уровня не описанных моделью (3).

Для построения модели «медленной» составляющей ТСКУМО применялся метод сдвиговой пошаговой множественной линейной регрессии, где в качестве предикторов выступали среднегодовые значения ТПО из архива ЕИ^Т. Использование значений поверхностной температуры за предшествующие периоды позволяет получить прогноз на несколько лет вперед.

Прогностическая модель компоненты Имед1 представлена в виде:

ЬМед1,1= Тзд-ф---, Т„,(1-т)), г=1 ,/и, (4)

где т - сдвиг в годах, т — число исходных точек.

Взаимнокорреляционная функция ТСКУМО [Ьеукиз е1 а1., 2012] от ТПО выявила, что коэффициенты корреляции остаются максимальными (0,86 - 0,90) до 8 сдвига, что> позволяет построить 8 моделей с заблаговременностью от 1 года до 8 лет. При включении в модели уже первых трех предикторов точность описания /¡„еа., достигает 88 - 92 % дисперсии, при этом стандартные ошибки моделей ст><х) изменяются от 1,8 мм/год (сдвиг 2 года) до 1,5 мм/год (сдвиг 6 лет).

Модель, полученная после осреднения оценок Имед, на каждом сдвиге т показала хорошую аппроксимацию на зависимом периоде с 1955 по 2007 гг., однако проверка на независимых данных выявила занижение ТСКУМО.

Тренд за период с 2008 по 2014 гг. (Тг = -0,18 мм/год) является незначимым по критерию Стьюдента. Так как «медленная» составляющая ТСКУМО в большей степени отвечает за повышение уровня, то полученные результаты означают значительное замедление темпа роста ТСКУМО

В разделе 4.5 «Прогностическая модель термостерического уровня Мирового океана» рассмотрена полная модель ТСКУМО, представляющая собой сумму «медленной» и «быстрой» компоненты.

При наложении двух компонент наблюдается небольшое снижение уровня, что определено значительной ролью «медленной» компоненты (рисунок 4), в то время как в оценках [Ье\киз й а1., 2012], наоборот, наблюдается слабый рост ТСКУМО. Оценки термостерических колебаний, полученные на основе альти-

метрических и GRACE данных, так же несколько расходятся с оценками [Levitus et al., 2012] и с 2005 г. показывают некоторое снижение термостериче-ского уровня.

Рисунок 4 - Временной ход термостерических колебаний уровня Мирового океана по данным [Levitus et al., 2012] (черная линия), суммарная модель (серая линия), оценки [Cazenave et al., 2009] по альтиметрическим и GRACE данным

(кружки).

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

Оптимальная прогностическая модель остаточного ряда с заблаговременностью в 1 год включает в себя 7 предикторов, при этом коэффициент детерминации составляет 0,84, а стандартная ошибка модели 1,35 мм, что меньше стандартного отклонения (3,22 мм). При использовании такой модели после добавления трендовой составляющей (рисунок 5), прогнозируемый на 2013 г. уровень термостерической компоненты составил 17,72 мм, при значении ТСКУМО на 2012 г. равном 16,05 мм. Таким образом, по полученным данным можно сделать вывод о том, что в 2013 г. ТСКУМО должно несколько возрасти.

Рисунок 5 - Временной ход термостерического уровня Мирового океана по данным [Ьеукиз е1 а1., 2012] (серая линяя) и прогностическая модель на 1 год (черная линяя). Пунктиром показаны независимые расчеты

Раздел 4.6 «Оценка термостерических колебаний уровня Мирового океана на основе спутниковых данных температуры поверхности океана». Оценка ТСКУМО проводилась на основе поверхностной температуры МО по спутниковым данным из архива 0188Т. При этом временной ряд ТСКУМО был представлен суммой двух компонент - тренда и отклонений от тренда:

ААстер = Тестер + ШТь АТ2,..., АТЯ), (5)

где - некая функция, подлежащая определению, Л Г, - значения аномалий ТПО в узлах географической сетки 0,25 х 0,25° архива ЫОАА ЫСОС ОИБТу.г, т - количество предикторов.

Исходный временной ряд ТСКУМО предварительно был разделен на 2 выборки: зависимую (1982 - 2007 гг.), по которой определялись параметры модели, и независимую (2008 - 2011 гг.), используемую для проверки ее качества.

С увеличением числа предикторов возрастает коэффициент детерминации и при этом моделью описываются не только закономерности, присущие случайной компоненте ТСКУМО, но и ее ошибки. Учитывая подобный характер распределения указанных параметров, было решено отобрать 5 моделей: с 3-го шага по 7-й. В этом случае Я2 меняется от 0,82 до 0,96, а сту(х) от 0,96 до 0,50 мм/год. Так как стандартное отклонение случайной компоненты ТСКУМО составляет 2,14 мм/год, то все модели имеют высокую точность, и, в качестве оптимальной, была принята средняя из пяти моделей.

В течение всех пяти лет (с 2007 по 2011 гг.) независимой выборки расхождения между фактическими и вычисленными по модели значениями ТСКУМО (рисунок 6) малы и носят случайный характер.

Рисунок 6 - Временной ход термостерического уровня Мирового океана по данным [ЬеуНиз е1 а1., 2012] - черная кривая и модель (5) - серая кривая. Темный квадрат - начало расчета аномалий ТСКУМО по независимым данным о

ТПО

В заключении сформулированы основные результаты и выводы диссертационной работы:

1. Рост теплосодержания Мирового океана за период с 1955 по 2012 гг. не является монотонным и на фоне основного тренда можно выделить промежутки времени, где изменения аномалий ТСМО имеют иной характер. Так, в течение 1955 - 1975 гг. тренд вообще отсутствовал, а интенсивный рост ТСМО в за 1993 - 2003 гг. сменился в 2004 - 2005 гг. его резким уменьшением, вследствие чего было высказано предположение о стабилизации теплосодержания Мирового океана. На основе вейвлет-анализа показано, что наиболее мощным является 23-летнее колебание ТСМО, период которого остается неизменным на протяжении всего рассматриваемого промежутка времени, 15-летнее колебание отмечается лишь в течение 1975 - 2005 гг. При этом оба цикла превышают 5% уровень значимости статистической модели «красный шум».

2. На основе данных реанализа (архив ЕГ^БЗТ) и спутниковых измерений (архив 0155Т) о температуре поверхности океана построены статистические модели с довольно высокой точностью описывающие изменения ТСМО. Во всех случаях стандартная ошибка моделей значительно меньше стандартного отклонения исходного временного ряда аномалий ТСМО. Их важным достоинством является простота, доступность и отсутствие необходимости в использовании глубоководных данных о температуре воды. Полученные по моделям результаты подтверждают вывод о начавшемся в последние годы некотором охлаждении толщи вод Мирового океана.

3. В ТСКУМО выделено три цикла, наиболее мощный из которых (57-летний) не меняется на протяжении всего интервала с 1880 по 2005 гг. Значительно слабее выражен 19 - 20-летний цикл, который присутствует в данных до середины 1980-х годов. Кроме того, можно отметить слабые 6 - 8-летние колебания, наблюдаемые с середины прошлого столетия. Вейвлет-анализ временного ряда ТПО выявил циклы с близкими периодами и небольшим опережением их максимумов относительно ТСКУМО. В частности, для основного 57-летнего цикла уровень отстает от ТПО на 3 — 4 года. Корреляция между вейвлет-коэффициентами на этих сдвигах составляет 0,99.

4. В общем случае формирование термостерической компоненты в конкретный момент времени может быть выражено суммой двух составляющих: «быстрой» и «медленной». «Быстрая» составляющая обусловлена внешним тепловым балансом океана, представляющим собой результирующий вертикальный поток тепла, который определяет колебания теплосодержания деятельного слоя океана. «Медленная» составляющая - нагревание (охлаждение) океана сверху и перераспределение тепла по всей его толще до нижней границы главного термоклина и, в основном, определяется трендовой компонентой.

5. Модель «быстрой» составляющей ТСКУМО была получена на основе 7 предикторов ТПО. Модель на зависимом периоде описывает 83 % дисперсии исходного ряда ТСКУМО, при этом стандартная ошибка составляет 1,28 мм. Стандартная ошибка модели для независимых данных меньше стандартного отклонения исходного временного ряда ТСКУМО. Таким образом, модель, полученная с использованием только 7 предикторов, хорошо описывает «быструю» компоненту ТСКУМО.

6. Построены прогностические модели «медленной» компоненты ТСКУМО с заблаговременностью от 1 года до 8 лет, в каждую из которых вошло по три предиктора. На независимом периоде (2008 - 2014 гг.) модель ос-редненных на каждом сдвиге т оценок Ижд., показала наличие незначимого тренда (-0,18 мм/год), что свидетельствует о замедлении роста ТСКУМО.

7. Предложена прогностическая модель ТСКУМО с заблаговременностью в 1 год, представляющая собой сумму трендовой компоненты и отклонений от нее. Оптимальная прогностическая модель для отклонений от тренда включает в себя 7 предикторов, коэффициент детерминации составляет 0,84, а стандартная ошибка модели 1,35 мм, что меньше стандартного отклонения остаточного ряда (3,22 мм). После наложения трендовой составляющей на 2013 г. прогнозируемый уровень термостерической компоненты составил 17,72 мм, при значении ТСКУМО на 2012 г. равном 16,05 мм.

8. По спутниковым данным о ТПО (архив ОГБЗТ) получена статистическая модель оценки ТСКУМО, представляющая собой осредненный вариант моделей с количеством предикторов от 3 до 7. В течение всех пяти лет независимой выборки (с 2007 по 2011 гг.) расхождения между фактическими и вычисленными по модели значениями ТСКУМО малы и носят случайный характер, при

этом стандартная ошибка «спрогнозированных» значений меньше стандартной ошибки исходной выборки.

Список публикаций по теме диссертации

Работы в журналах, рекомендованных ВАК РФ для публикации основных

результатов диссертаций:

1. Малинин В.Н., Глок Н.И. К оценке стерических колебаний уровня Мирового океана // Ученые записки РГГМУ. 2009. №10. С. 53-62.

2. Малинин В.Н., Гордеева С. М., Глок Н.И. О возможных изменениях стерической компоненты уровня Мирового океана // Изв. РГО. 2010. Т. 142. Вып. 4. С. 23-32.

3. Глок Н.И., Малинин В.Н. Статистический анализ стерических колебаний уровня Мирового океана// Ученые записки РГГМУ. 2011. № 21. С. 26-36.

4. Глок Н.И., Малинин В.Н. Изменения теплосодержания Мирового океана в период 1955-2009.//Фундаментальная и прикладная гидрофизика. 2013. Том 6, №1 С.32-39.

5. Малинин В.Н., С.М.Гордеева, Н.И. Глок. К оценке теплосодержания Мирового океана по спутниковым данным о температуре поверхностного слоя воды // Современные проблемы дистанционного зондирования Земли из космоса. 2013. Т. 10. № 3. С. 201-207.

6. Малинин В. Н., Глок Н.И. Использование спутниковых данных о температуре поверхностного слоя воды для оценки стерических колебаний уровня Мирового океана // Исследование Земли из космоса. 2014. № 2 (в печати).

Публикации в прочих научных изданиях:

7. Малинин В.Н., Глок Н.И. Стерический эффект в колебаниях уровня Мирового океана // Географ, образование и наука в России: история и современное состояние. СПб, 2009, Изд. ВВМ, с. 827-836.

8. Клиге Р.К., Малинин В.Н., Гордеева С.М, Глок Н.И. Факторы изменения

уровня океана // Современные глобальные изменения природной среды: М. Научный мир, 2012, С. 302-318.

Глок Наталья Ивановна

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата географических наук

СТАТИСТИЧЕСКИЙ АНАЛИЗ МЕЖГОДОВОЙ ИЗМЕНЧИВОСТИ ТЕПЛОСОДЕРЖАНИЯ И СТЕРИЧЕСКИХ КОЛЕБАНИЙ УРОВНЯ МИРОВОГО ОКЕАНА

_ЛР№ 020309 от 30.12.96_

Подписано в печать 07.11.13. Формат 60x90 1/16. Гарнитура Times New Roman. Печать цифровая. Усл.печ.л. 1,37. Тираж 100 экз. Заказ № 227. РГГМУ, 195196, Санкт-Петербург, Малоохтинский пр., 98. _Отпечатано в ЦОП РГГМУ_

Текст научной работыДиссертация по наукам о земле, кандидата географических наук, Глок, Наталья Ивановна, Санкт-Петербург

ГОУ ВПО РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ (РГГМУ)

На правах рукописи

04201 4 53 5 <со ГЛОК НАТАЛЬЯ ИВАНОВНА

СТАТИСТИЧЕСКИЙ АНАЛИЗ МЕЖГОДОВОЙ ИЗМЕНЧИВОСТИ ТЕПЛОСОДЕРЖАНИЯ И СТЕРИЧЕСКИХ КОЛЕБАНИЙ УРОВНЯ

МИРОВОГО ОКЕАНА

Специальность 25.00.28 - океанология

Диссертация на соискание ученой степени кандидата географических наук

Научный руководитель:

доктор географических наук, профессор

Малинин Валерий Николаевич

Санкт-Петербург 2013

СОДЕРЖАНИЕ

10

10

14

35

Стр.

Сокращения 4

Введение 5

1. Современное состояние изученности теплосодержания и стерических колебаний уровня Мирового океана

1.1 Основные факторы формирования колебаний уровня Мирового океана

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

1.3 Краткий обзор исследований теплосодержания Мирового океана 27

1.4 Краткий обзор исследований стерических колебаний уровня Мирового океана

2 Статистические методы анализа данных, используемые в работе 44

2.1 Анализ трендовой компоненты временного ряда 44

2.2 Анализ циклических колебаний временного ряда 47

2.3 Множественный регрессионный анализ 57

3 Теплосодержание Мирового океана и его межгодовая

68

изменчивость

3.1 Закономерности колебаний теплосодержания Мирового океана 68

3.2 Статистическая модель оценки среднегодовых значений теплосодержания Мирового океана на основе реанализа 71 температуры его поверхности

3.3 Статистическая модель оценки теплосодержания Мирового

океана на основе спутниковых данных температуры 78

поверхностного слоя воды

4 Физико-статистическая модель термостерического уровня Мирового океана

4.1 Закономерности термостерических колебаний уровня Мирового океана

4.2 Физическое обоснование модели 89

83 83

4.3 Модель «быстрой» компоненты термостерических колебаний уровня Мирового океана

4.4 Модель «медленной» компоненты термостерических колебаний уровня Мирового океана

4.5 Прогностическая модель термостерического уровня Мирового океана

4.6 Оценка термостерических колебаний уровня Мирового океана на основе спутниковых данных температуры поверхности океана Заключение 113 Список использованной литературы 116

96 101 105 109

СОКРАЩЕНИЯ

УМО - Уровень Мирового океана;

ВМО - Всемирная Метеорологическая Организация;

ЮНЕП - Программа ООН по окружающей среде

МГЭИК - Межправительственная группа экспертов по изменению

(IPCC) климата;

МО - Мировой океан;

СКУМО - Стерические колебания уровня Мирового океана;

ТСМО - Теплосодержание Мирового океана;

ТПО - Температура поверхности океана;

ТСКУМО - Термостерические колебания уровня Мирового океана;

MJIP - Множественная линейная регрессия;

ПТВ - Приповерхностная температура воздуха;

АТПО - Аномалии температуры поверхности океана;

GRACE - Gravity Recovery And Climate Experiment;

WOD - World Ocean Database;

WOA - World Ocean Atlas.

ВВЕДЕНИЕ

Актуальность темы. Климат на нашей планете постоянно меняется, и, в первую очередь, это вызвано резким изменением глобальной температуры воздуха, которая в течение XX века увеличилась на 0,6 - 0,7 °С. Учитывая тесную взаимосвязь атмосферы и океана, процессы, происходящие в ней, не могут не сказываться на водной оболочке нашей планеты. Одновременно процессы в океане оказывают стимулирующее влияние на тепловой режим и общую циркуляцию атмосферы. В этом плане особую роль играет теплосодержание Мирового океана (ТСМО), которое является важнейшим фактором изменений глобального климата, так как он имеет высокую теплоемкость и инерционность. Кроме того, ТСМО практически полностью определяет стерические колебания уровня Мирового океана (СКУМО), вклад которых в формирование тренда глобального уровня может составлять от 15 до 50 % в зависимости от интервала временного осреднения.

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

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

Исходя из поставленной цели, были сформулированы следующие задачи:

- выявление закономерностей межгодовой изменчивости

теплосодержания Мирового океана;

- построение статистических моделей межгодовой изменчивости теплосодержания Мирового океана по данным спутников и реанализа о температуре поверхностного слоя воды;

- выявление закономерностей термостерических колебаний уровня Мирового океана (УМО);

- физическое обоснование и построение физико-статистической модели термостерических колебаний уровня Мирового океана (ТСКУМО);

- построение статистической модели термостерических колебаний уровня Мирового океана по реанализу и спутниковым данным температуры его поверхности.

Материалы и методы. Основой работы послужили следующие архивы гидрометеорологических данных, находящихся в свободном доступе в сети Интернет:

- среднемесячные данные о поверхностной температуре океана в узлах широтно-долготной сетки 2 х 2° из глобального архива NOAA NCDC ERSST (National Oceanic and Atmospheric Administration National Climatic Data Center Extended Reconstructed Sea Surface Temperature);

- суточные данные о поверхностной температуре океана в узлах широтно-долготной сетки 0,25 х 0,25° из глобального архива NOAA NCDC OISST (National Oceanic and Atmospheric Administration National Climatic Data Center Optimum Interpolation Sea Surface Temperature);

- среднемесячные данные о глубоководной температуре и вертикальных скоростях воды в узлах широтно-долготной сетки 0,5 х 0,5° из глобального архива SODA (Simple Ocean Data Assimilation);

-тренды поверхностной температуры океана из архива HadlSSTl (Hadley Central Sea Ice and Sea Surface Temperature Dataset) в узлах географической широтно-долготной сетки 2 х 2°.

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

6

анализ.

Научная новизна:

- выявлены закономерности межгодовой изменчивости теплосодержания Мирового океана. Показано, что за период с 1955 по 2010 гг. трендовая составляющая является значимой и носит положительный характер, однако, начиная с 2003 г. происходит значительное замедление роста теплосодержания. В данных ТСМО обнаружены два устойчивых цикла с периодами в 23 и 15 лет;

-установлены особенности распределения трендов временного хода температуры воды в отдельных океанах в слое 0-3000 м. Для всех океанов тренды являются значимыми до глубины порядка 1500 метров. При этом на разных горизонтах отмечаются как положительные, так и отрицательные тренды. Однако, начиная с 500 м, проявляется тенденция к уменьшению величины положительного тренда с глубиной;

- разработаны статистические модели межгодовых изменений теплосодержания Мирового океана по данным спутников и реанализа о температуре поверхностного слоя воды;

- выявлены закономерности термостерических колебаний уровня Мирового океана. Показано, что в термостерических колебаниях уровня Мирового океана с 1880 по 2005 гг. присутствует значимое колебание с периодом в 57 лет, которое отчетливо проявляется и для температуры поверхности океана. Обнаружены также колебания с периодами 22 - 25 лет, 15 и 7 - 8 лет. Установлено, что трендовая компонента обеспечивает максимальный вклад в дисперсию термостерического ряда, который составляет не менее 77 %, причем на долю значимых циклов приходится в среднем около 10 %;

- разработана физико-статистическая модель термостерических колебаний УМО, представляющая собой сумму двух компонент: «быстрой», отвечающей за внешний тепловой баланс океана и «медленной», определяющей перераспределение тепла в толще воды до нижней глубины

7

главного термоклина;

- разработана статистическая модель термостерических колебаний УМО на основе спутниковых данных о температуре поверхности океана (ТПО);

На защиту выносятся:

- закономерности изменений теплосодержания и термостерических колебаний уровня Мирового океана в современных климатических условиях;

- статистические модели межгодовых изменений теплосодержания океана по данным спутников и реанализа температуры поверхностного слоя воды;

-физико-статистическая модель термостерических колебаний уровня Мирового океана;

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

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

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

Апробация работы и публикации. Данная работа была выполнена в рамках Аналитической ведомственной целевой программы «Оценить последствия повышения уровня океана в условиях современных изменений климата» №2.1.1/6690 и Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (государственные контракты №П726 и № 02.740.11.0849), а также в рамках гранта Правительства РФ (Договор № 11.G34.31.0078) для поддержки исследований под руководством ведущих ученых. Основные результаты

8

диссертации были представлены в рамках шестнадцатой Санкт-Петербургской ассамблеи молодых ученых и специалистов и отмечены грантом губернатора Санкт-Петербурга.

По теме диссертации опубликовано 8 работ, в том числе 6 из них в рецензируемых журналах по перечню ВАК.

Структура и объем работы. Диссертация состоит из введения, 4 глав, заключения и библиографического списка использованной литературы. Объем диссертации составляет 125 страницы, включая 13 таблиц и 33 рисунка. Список литературы составляет 100 наименований.

Благодарности. Автор искренне благодарит научного руководителя профессора В.Н. Малинина за конструктивное руководство и всестороннюю поддержку в работе, а так же к.г.н. С.М. Гордееву за консультации и помощь.

1 Современное состояние изученности теплосодержания и стерических колебаний уровня Мирового океана

1.1 Основные факторы формирования колебаний уровня Мирового океана

На формирование уровня Мирового океана (УМО) воздействуют три группы различных факторов: космогеофизические, геолого-геодинамические и гидрометеорологические [27, 13].

К космогеофизическим относятся приливообразующие силы Луны и Солнца, свободные и вынужденные колебания полюсов Земли, неравномерные изменения скорости вращения Земли и астрономические факторы, связанные с изменением орбитальных параметров Земли, её положением в Солнечной системе и т.п. Приливообразующие силы формируют, прежде всего, полусуточные и суточные колебания уровня, которые, вследствие своей периодичности, легко фильтруются при осреднении уровенных данных уже за месяц. Наложение свободных (14 мес) и вынужденных (12 мес) колебаний полюсов Земли вызывает явление резонанса, в результате чего возникает так называемый «полюсный» прилив с периодом 6-7 лет, амплитуда которого может достигать нескольких сантиметров. Неравномерные изменения скорости вращения Земли обуславливают «ротационный» прилив, амплитуда которого не превышает нескольких миллиметров. Что касается астрономических факторов, то они вызывают колебания уровня с периодами от тысяч до миллионов лет.

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

10

Эти процессы приводят либо к эпизодическим импульсивным изменениям уровня в виде волн цунами и крупных сейш, либо к медленным межвековым изменениям с масштабами, измеряемыми тысячелетиями и миллионами лет [34].

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

Применительно к уровню Мирового океана космогеофизическими факторами, как правило, пренебрегают, так как из-за периодичности большинства процессов их суммарное воздействие незначительно, а временные масштабы могут значительно превышать период инструментальных наблюдений. Деформационными колебаниями так же можно пренебречь вследствие того, что тектонические движения земной коры имеют разнонаправленный характер в различных регионах Земли и при глобальном осреднении суммарный эффект их воздействия близок к нулю. Что касается донного осадконакопления, то его величина для масштабов времени не превышающих нескольких десятков лет, является пренебрежимо малой. Вклад анемоборических колебаний в общий уровень при глобальном осреднении по МО так же не является значительным. Таким образом, межгодовые колебания УМО в современных климатических условиях практически полностью обусловлены изменениями составляющих глобального водного баланса и плотности морской воды. [17, 18, 28].

11

Для оценки состояния климата в 1988 г. Всемирной Метеорологической Организацией (ВМО) совместно с Программой Организации Объединенных Наций по окружающей среде (ЮНЕП) была учреждена Межправительственная группа экспертов по изменению климата (МГЭИК), основной задачей которой является обобщение и анализ имеющейся научной и социально-экономической информации об изменениях климата. С момента создания МГЭИК было подготовлено пять отчетов, в которых специальные разделы были посвящены изменениям УМО, так как уровень Мирового океана является неотъемлемой частью глобальной климатической системы и главным индикатором глобального изменения климата [41].

В таблице 1.1 представлены оценки трендов различных факторов, влияющих на формирование тренда УМО, по данным Третьего отчета МГЭИК [65] за период с 1910 по 1990 гг., Четвертого отчета МГЭИК [92] за периоды с 1961 по 2003 гг. и с 1993 по 2003 гг., Пятого отчета МГЭИК [94] с 1993 по 2010 гг. и с 2005 по 2010 гг. и по данным [50] за период с 2003 по 2008 гг.

Суммарная оценка трендовых составляющих представленных процессов за период времени с 1910 по 1990 гг. имеет большую погрешность. Это связано с тем, что адекватность информации о Мировом океане в начале прошлого столетия была не высока. Но даже при таких грубых расчетах тренд временной изменчивости стерических колебаний среди других факторов является наибольшим. Максимальное влияние на УМО стерические колебания оказывали в период с 1993 по 2003 гг., при этом они описывали более 50 % тренда наблюденного уровня. Также немаловажную роль на формирование УМО, по данным МГЭИК, оказывает таяние ледников и ледниковых шапок, которое вместе со стерическими колебаниями для первых трех периодов (см. таблицу 1.1) описывают более 80% от суммарного тренда.

Таблица 1.1 - Тренды различных компонент уровня Мирового океана, мм/год

Факторы, влияющие на УМО Оценки трендов за различные периоды

19101990 гг. [651 19612003 гг. [92] 19932003 гг. [92] 19932010 гг. [94] 20032008 гг. [50] 20052010 гг. [94]

Стерические колебания 0,51 ± 0,20 0,42 ± 0,12 1,6 ± 0,5 1,0 ± 0,3 0,31 ± 0,15 0,9 ± 0,3

Ледники и ледниковые шапки 0,30 ± 0,09 0,5 ± 0,18 0,77 ± 0,22 1,04 ± 0,19 1Д ± 0,24 0,30 ± 0,03

Гренландский ледниковый щит 0,06 ± 0,05 0,05 ± 0,12 0,21 ± 0,07 0,34 ± 0,14 1,00 ± 0,15 0,63 ± 0,15

Антарктический ледниковый щит -0,09 ± 0,10 0,14 ± 0,41 0,21 ± 0,35 0,18 ± 0,09 0,31 ± 0,16

Другие факторы* -0,07 ± 0,78 — — 0,38 ± 0,11 0,17 ±0,10 -0,33 ± 0,22

Сумма 0,70 ± 0,82 1Д ± 0,5 2,8 ± 0,7 2,9 ± 0,4 2,58 ± 0,34 1,8 ± 0,4

Наблюденный уровень МО 1,5 ± 0,50 1,8 ± 0,5 ЗД ± 0,7 3,2 ± 0,4 2,5 ± 0,4 2,3 ± 0,6

* Примечание — в определение другие факторы: за период с 1910 по 1990 гг. вошли: вечная мерзлота, донное осадконакопление и материковый сток, для остальных периодов -только материковый сток.

В XXI веке отмечается замедление роста стерического уровня, и при этом резко возрастает роль