Бесплатный автореферат и диссертация по наукам о земле на тему
Межгодовая изменчивость уровня мирового океана в современных климатических условиях
ВАК РФ 25.00.28, Океанология
Автореферат диссертации по теме "Межгодовая изменчивость уровня мирового океана в современных климатических условиях"
□□3432572
РОССИЙСКИЙ ГОСУДАРСТВЕННЫЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ
На правах рукописи
УДК 551.461.22: 551.583
ШЕВЧУК ОЛЕГ ИГОРЕВИЧ
МЕЖГОДОВАЯ ИЗМЕНЧИВОСТЬ УРОВНЯ МИРОВОГО ОКЕАНА В СОВРЕМЕННЫХ КЛИМАТИЧЕСКИХ УСЛОВИЯХ
Специальность 25.00.28- «Океанология»
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата географических наук
Санкт-Петербург 2009
003492572
Работа выполнена в Российском государственном гидрометеорологическом университете (РГГМУ)
Научный руководитель:
доктор географических наук, Малинин Валерий Николаевич профессор
Официальные оппоненты:
доктор географических наук, Смирнов Николай Павлович профессор
доктор географических наук Добровольский Сергей Гавриилович
Ведущая организация: Арктический и антарктический научно-исследовательский институт (ААНИИ)
Защита состоится " 24 " декабря 2009 г. в 15 час. 30 мин. на заседании диссертационного совета Д 212.197.02 в Российском государственном гидрометеорологическом университете по адресу: 195196, Санкт-Петербург, Малоохтинский проспект, 98, актовый зал.
С диссертацией можно ознакомиться в библиотеке Российского государственного гидрометеорологического университета.
Автореферат разослан "£Ъ " ноября 2009 г.
Ученый секретарь диссертационного совета, кандидат географических наук
В.Н. Воробьев
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Уровень Мирового океана является важнейшей характеристикой, определяющей как гидродинамические процессы в самом океане, так и характер взаимодействия океана и атмосферы. Изменения уровня Мирового океана тесно связаны с современными колебаниями климата, которые в свою очередь определяют испарение, осадки, стерические колебания. Эвстати-ческие колебания можно рассматривать как индикатор глобального водообмена, перераспределяющего воды гидросферы между отдельными оболочками Земли и, прежде всего, между водами литосферы и криосферы.
Как следует из докладов Межправительственной группы экспертов по изменениям климата (МГЭИК) (1РСС, 2001, 2007), в соответствии с некоторыми сценариями изменений климата, основанными на результатах его численного моделирования, уровень Мирового океана может повыситься на 10—30 см к 2030 г. и в пределах 30-100 см (наиболее вероятно 65 см) к концу столетия. Естественно, если такое развитие изменений климата станет реальностью, то это грозит катастрофическим ущербом для инфраструктуры прибрежных территорий, где проживает около миллиарда жителей Земли. Возможно подтопление крупнейших городов мира: Лондона, Нью-Йорка, Токио, С-Петербурга и др. В зависимости от разных причин темпы и масштабы подъема уровня моря могут значительно превзойти способность антропогенных и естественных систем к адаптации.
Поэтому задача изучения долговременных колебаний уровня Мирового океана и особенно построения систем прогноза уровня Мирового океана в системе океан-атмосфера представляется чрезвычайно актуальной. Эксперты МГЭИК отнесли данный вопрос к числу важнейших проблем современной гидрометеорологии.
Цель работы: изучение межгодовых колебаний уровня Мирового океана и эвстатических факторов его формирования.
Для достижения указанной цели решались следующие задачи:
1. Анализ межгодовых колебаний уровня океана за период инструментальных наблюдений.
2. Выявление станций с многолетними наблюдениями за морским уровнем, обеспечивающих минимальную ошибку восстановления уровня Мирового океана (УМО) с достаточной для практических целей точностью и построение статистической модели оценки межгодовых колебаний УМО.
3. Анализ межгодовой изменчивости эвстатических факторов (испарения и осадков) формирования уровня Мирового океана. Построение диагностической статистической модели оценки межгодовых колебаний УМО на основе эвстатических факторов.
4. Сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XX столетия.
5. Оценка связи колебаний уровня Мирового океана с приповерхностной температурой воздуха над океаном.
6. Построение физико-статистической модели прогноза межгодовых колебаний уровня Мирового океана на ближайшие десятилетия.
7. Построение адаптивной модели прогноза сезонной изменчивости уровня Мирового океана по альтиметрическим данным.
Материалы и методы. Основой работы являлись несколько глобальных архивов гидрометеорологических данных находящихся в свободном доступе в Интернет:
- футшточные данные об уровне Мирового океана полученные из архива PSMSL (Permanent Service for Mean Sea Level);
- альтиметрические данные центра по хранению, верификации и интерпретации спутниковой океанографической информации (AVISO) и Колорадского центра астродинамических исследований (CCAR) при Колорадском университете со спутников Topex/Poseidon и Jason-1, с пространственным разрешением 1/3°х 1/3°;
- среднемесячные значения осадков и испарения в узлах первичной широтно-долготной сетки 1.875°х1.875° из глобального гидрометеорологического архива CDAS (Climate Data Assimilation System);
- аномалии поверхностной температуры воздуха из архива HadCrut3 с пространственным разрешением 5°х5°;
Для выполнения поставленных задач использовался комплекс методов многомерного статистического анализа (MCA), а также теория проверки статистических гипотез и анализ временных рядов.
Научная новизна:
- Впервые построена статистическая модель, обеспечивающая минимальную ошибку восстановления межгодовых колебаний уровня Мирового океана по данным 5 береговых станций с достаточной для практических целей точностью.
- Установлены закономерности межгодовых колебаний УМО за период инструментальных наблюдений.
- Впервые выявлены особенности межгодовой изменчивости испарения, осадков и эффективного испарения над отдельными океанами и Мировым океаном в целом в период интенсивного глобального потепления климата.
- Установлена степень взаимосвязи колебаний УМО с температурой воздуха над его акваторией.
- Впервые построена физико-статистическая модель прогноза УМО на ближайшие десятилетия на основе данных о температуре воздуха над его акваторией.
- Впервые построена адаптивная модель прогноза сезонной изменчивости УМО по альтиметрическим данным.
На защиту выносятся:
1. Закономерности межгодовых колебаний уровня Мирового океана за период инструментальных наблюдений.
2. Особенности межгодовой изменчивости эвстатических факторов (осадки, испарение, эффективное испарение) над океанами.
3. Результаты прогнозирования уровня Мирового океана по данным о температуре воздуха над его акваторией на ближайшие десятилетия.
4. Адаптивная модель прогноза сезонной изменчивости уровня Мирового океана по альтиметрическим данным.
Практическая значимость работы. Полученные в работе результаты могут быть использованы при решении задач мониторинга УМО и моделирования климата, а также при разработке долгосрочных прогнозов межгодовых колебаний и сезонного хода УМО.
Личный вклад автора. Подготовка архивов информации, проведение статистических расчетов, выполнение необходимых графических построений, обобщение и анализ результатов.
Апробация работы и публикации. Основные положения и отдельные результаты обсуждались на итоговых сессиях Ученого Совета РГГМУ, 2005 и 2006; на международной школе-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда», РГТМУ, 2005; научном семинаре «Проблемы изменения климата в полярных широтах» в рамках проекта BRIDGE «Keeping it Cool», Университет г. Шеффилда, Великобритания, 2007; международной научно-практической конференции «Географическое образование и наука в России: история и современное состояние», СПБГУ, 2008; всероссийской научной конференции «Селиверстовские чтения», 2009.
По теме диссертации опубликовано 9 печатных работ, из них 2 статьи в рецензируемых журналах (по перечню ВАК).
Структура и объем работы. Диссертация состоит из введения, 3 глав, заключения и библиографического списка использованной литературы. Объем диссертации составляет 113 страниц, включая 21 таблицу и 34 рисунка. Список литературы составляет 84 наименования, в том числе 43 на иностранных языках.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность темы, формулируются цель и задачи исследования.
Первая глава «Колебания уровня Мирового океана в современных климатических условиях по данным наблюдений» посвящена реконструкции и исследованию межгодовой изменчивости уровня Мирового океана.
В разделе 1.1 «Краткий обзор исследований колебаний уровня океана на основе футшточных наблюдений» рассматриваются исследования колебаний уровня океана на основе футшточных наблюдений по литературным данным
(Gutenberg, 1941; Fairbridge, Krebs, 1962; Калинин, Книге, 1972; Дуванин и др., 1957; Клиге 1985; Church et al., 2006; Jevrejeva et. al., 2006 и др.).Показано, что в современную эпоху наблюдается повсеместное повышение УМО.
В общем, в XX в. рост уровня согласно обобщенным оценкам многих исследователей, приведенных в докладе Межправительственной группой экспертов по изменениям климата в 2001 г. составил 1,5-2 мм/год. Темпы роста в последние десятилетия заметно увеличились. В частности за время интенсивного потепления климата с 1985 по 2001 гг. по данным Church et al. (1985 -2001 гг.) - 3,27 мм/год, Jevrejeva et. al. - 3,68 мм/год.
Раздел 1.2 «Межгодовая изменчивость уровня Мирового океана по данным футшточных наблюдений» посвящен реконструкции и анализу УМО с 1861 - 2005 гг. Разными авторами (Клиге, 1985; Barnett, 1984; Douglas, 1997; Fairbridge et al., 1962; Gornits et al., 1982) расчет УМО осуществлялся за различные, но не одинаковые промежутки времени. Поэтому было выполнено арифметическое осреднение указанных рядов за общий период с 1901 по 1958 г. Полученный таким образом ряд УМО имеющий довольно короткую длину (рис. 1, линия 1) считался базовым. После этого решалась задача его удлинения: вперед с 1958 г. и назад до 1901 г. С этой целью использовался глобальный архив футшточных (береговых) наблюдений уровня моря - Permanent Service for Mean Sea Level (PSMSL), находящийся в свободном доступе в сети Интернет. Реконструкция УМО в обе стороны осуществлялась методом пошаговой множественной регрессии, который позволяет последовательно отбирать наиболее значимые факторы и подробно исследовать их статистическую значимость. При этом линейная модель для расчета УМО (Ьумо) имеет вид:
где А/ - уровень моря на j станции; т — число станций в модели.
В состав предикторов была включена 21 станция, длина наблюдений на которых за уровнем превышает 90 лет, причем в качестве зависимой выборки использовался период с 1909 по 1958 гг., т.е. 50 лет. Коэффициенты модели множественной линейной регрессии находились с помощью метода включения переменных. Порядок модели (ее оптимальность) оценивался по минимуму среднеквадратической ошибки при условии значимости всех коэффициентов регрессии по критерию Стьюдента и значимости модели по критерию Фишера (при уровне значимости 5%). Было получено, что оптимальной является модель, для которой ш=5. В порядке значимости, то есть по вкладу в описание изменчивости исходного ряда hyMo, в данную модель входят: New York, Seattle, Galveston, Swinoujscie и Ratan. При этом уже первый коэффициент корреляции УМО со ст. New York достигает 0,91. Коэффициент детерминации данной мо-
О)
дели равен Л2 = 0,94, то есть модель описывает 94 % дисперсии функции отклика, а стандартная ошибка модели равна ау(х)= 6,6 мм.
Рисунок 1 - Межгодовая изменчивость уровня Мирового океана (1 - базовый уровень; 2 - реконструкция уровня вперед; 3 - реконструкция уровня назад). В качестве доверительного интервала принято значение ±а
На зависимых данных модель (I) с достаточно высокой точностью описывает глобальную изменчивость исходного ряда УМО. Принимая во внимание данное обстоятельство, нами был произведен расчет УМО за независимый промежуток времени с 1959 по 2005 гг. В результате этого общая продолжительность временного ряда УМО возросла до 105 лет (с 1901 по 2005 гг.). Его межгодовой ход также приведен на рисунке 1 (линия 2).
Далее решалась задача реконструкции УМО назад до 1901 г. (рис. 1, линия 3). С этой целью было отобрано 8 станций, для которых общее начало наблюдений за уровнем начинается с 1861 г., а общий конец наблюдений датируется 1934 г. Таким образом, зависимая выборка имеет длину 34 года (с 1901 по 1934 гг.). Использование аппарата пошаговой регрессии позволило установить, что оптимальной также следует считать модель с пятью предикторами. В порядке значимости в нее входят следующие станции: Brest, Lyokki, Swinoujscie, San Francisco, Lypyrtti. Первый коэффициент корреляции УМО со ст. Brest равен г = 0,63, коэффициент детерминации данной модели равен R = 0,75, стандартная ошибка модели cryW=6,5 мм. При этом коэффициенты регрессии значимы по критерию Стьюдента, а сама модель значима по критерию Фишера (при уровне значимости 5%). Заметим, что хотя доля описания дисперсии исходного УМО по данной модели существенно меньше по сравне-
нию с моделью реконструкции УМО вперед, однако оценки случайной ошибки стУ(Х) практически одинаковы.
В таблице 1 приводятся более подробные сведения об оценках линейного тренда и статистических характеристиках рассматриваемых рядов УМО в XX в. по независимым данным разных авторов. Отметим, что в работах Church et. al., (2006) и Jevrejeva et. al., (2006) УМО определялся осреднением большого числа береговых станций из архива PSMSL методами пространственной интерполяции (например, в работе Jevrejeva et. al., использовалось 1023 станции). Установлено, что расхождения между средними значениями и дисперсиями по критериям Стьюдента и Фишера являются незначимыми. Оценки трендов тоже очень близки друг к другу.
Таблица 1 - Оценки линейного тренда и статистических характеристик рядов
УМО за 1901- 1999 гг.
Авторы R2 тренд Среднее значение o£/Xcp
Малинин, Гордеева, Шевчук, 2007 0,96 1,79 10,40 170,17 0,06
Church et. al., 2006 0,97 1,65 9,09 166,93 0,05
Jevrejeva et. al., 2006 0,94 1,89 13,77 167,73 0,08
Таким образом, выявлены станции с многолетними наблюдениями за морским уровнем, обеспечивающие минимальную ошибку восстановления УМО с достаточной для практических целей точностью. А результаты сравнения сделанных выше расчетов с расчетами других авторов свидетельствуют, что выполненная автором реконструкция УМО является весьма точной, отражает общие закономерности его межгодовых колебаний и совместима с аналогичными данными других авторов. Однако принципиальное ее отличие заключается в том, что оценки УМО могут быть получены по регрессионной модели с малым числом станций.
В разделе 1.3 «Изменчивость уровня Мирового океана по алътиметри-ческим данным» рассмотрен межгодовой ход УМО за период 1993-2008 гг., построенный по среднемесячным альтиметрическим данным с изомаршрутных спутников Торех/Розе1с1оп и 1азоп-1 с учетом поправки обратного барометра в университете Колорадо. Главной его особенностью является наличие ярко выраженного тренда, равного 7У=3,2+0,4 мм/год, который после исключения поправки на гляциоизостатический эффект поверхности океанов, становится равным 7>=2,9±0,4 мм/год.
Выполнено совмещение УМО по альтиметрическим и береговым наблюдениям, которое возможно двумя простыми способами. Первый - по совмещению средних значений двух рядов. Второй способ заключается в построении модели регрессии между альтиметрическим и футшточным рядами УМО.
В итоге показано, что второй способ по существу представляет собой тренд, свойственный фухшточному ряду УМО за период 1993-2005 гг., в то время как первый способ и соотношение Х0^ = 258 + Х111Т, связывающее альтимет-рические и футшточные данные, можно рассматривать в виде основы для построения мониторинга УМО по спутниковой информации.
Глава 2 «Факторы формирования уровня Мирового океана и их анализ».
В разделе 2.1 «Эвстатические изменения уровня океана как интегральный индикатор глобального водообмена» показано, что УМО отражает процессы перераспределения природных вод в гидросфере, представляющей сплошную оболочку Земли, содержащую воду во всех агрегатных состояниях в пределах Мирового океана, криосферы, литосферы и атмосферы и принимающую непосредственное участие в планетарном круговороте влаги. Поэтому эвстатические колебания УМО можно рассматривать как интегральный индикатор глобального водообмена.
Раздел 2.2 «Эвстатические колебания уровня Мирового океана» посвящен исследованию эвстатических колебаний УМО. В частности в подразделе 2.2.1 «Вертикальный влагообмен в системе океан-атмосфера» рассматривается глобальный гидрометеорологический архив СБАБ, являющегося частью системы «ретроспективного» анализа» (реанализа) >ЮАА ИСЕР/ЫСАК Яеапа1ук15 (Ка1пау ег а!., 1996) и находящий в свободном доступе в сети Интернет, из которого заимствованы поля среднемесячных значений испарения, осадков и их разности ЕЭ^=Е-Р (эффективное испарение) над акваторией Мирового океана (МО), исключая Северный Ледовитый океан. Затем было осуществлено осреднение исходных данных по пятиградусным трапециям акватории МО, что позволило сократить число узлов с 15965 до 1305.
В качестве основного периода исследования выбран промежуток с 1980 по 2005 гг., характеризующийся особенно резким потеплением глобального климата. Для каждого из трех океанов рассчитаны зонально-осредненные значения испарения, осадков и эффективного испарения, которые затем дополнительно усреднены для отдельных океанов и МО в целом. Установлено, что в межгодовом ходе испарения и осадков присутствует сильный положительный тренд, который над МО в осадках составляет 4,2 мм/год, а для испарения - 3,6 мм/год. Главной причиной их формирования является интенсивное потепление климата в рассматриваемый период времени, которое характеризуется трендом в глобальной температуре воздуха равным Тг = 0.18 °С АТ/10лет. Показано, что корреляция глобальной температуры воздуха с рядами испарения и осадков для МО, составляет соответственно Гв,т~ 0,82 и гр г= 0,80.
Поскольку количество выпавших осадков над МО растет быстрее испарения, то тренд в эффективном испарении оказывается отрицательным (рис. 2). Наибольшая величина его отмечается для Атлантического и Индийского океанов. Для Тихого океана тренд незначим по критерию Стьюдента. Величина тренда в межгодовом ходе эффективного испарения для МО в целом составляет Тг =-0,6 мм/год. Итак, уменьшению эффективного испарения на 0,6 мм/год со-
ответствует увеличение УМО на 0,6 мм/год или на 217 км3/год. Поскольку общий тренд УМО за рассматриваемый промежуток времени составляет около 2,4 мм/год, то вклад эффективного испарения в изменения УМО достигает 25 %.
см
----
• Л Л
* \ . -<
Л / •
ЕЭ(Ь =-0,061 + 142 Я2 = 0,15
годы
1980
1985
1990
1995 2000
2005
Рисунок 2 - Межгодовой ход эффективного испарения над отдельными океанами за период 1980-2006 гг. в см/год. (1 - Атлантический, 2 — Индийский, 3 -Тихий, 4 - Мировой океан)
Межгодовая изменчивость испарения и осадков над отдельными океанами и МО в целом примерно одинакова и невелика, коэффициент вариации составляет С=0,2-0,4. Существенно более высокой оказывается изменчивость Еэф, так как его средние значения существенно меньше средних значений Е и Р. Самый большой коэффициент вариации отмечается для Тихого океана (С=0,42). Для МО в целом С=ОДО. Заметим, что он значительно больше коэффициента вариации притока пресных вод к МО (С=0,02).
Далее построена диагностическая модель оценки внутригодовых изменений уровня АЛм, определяемых по центральным разностям, в зависимости от Еъф в узлах сеточной области Мирового океана
= а0 + а, Езф[ + а2Езф2 +... + а10Еэфю (2)
С помощью данной модели легко вычислить оценки А/гм, а затем значения среднего годового УМО как /¡м(Н1)= Лм(И) + 2Л/гм„ межгодовой ход которого приведен на рисунке 3.
330 г мм
310
290
270
250
230
210
190
' годы
1980 1985 1990 1995 2000 2005 2010
Рисунок 3 - Изменения УМО за период 1981-2006 гг., рассчитанные по различным вариантам (5, 8 и 10 предикторов) модели (2) по данным об эффективном испарении с поверхности океанов
В подразделе 2.2.2 «Приток пресных вод к океану. Айсберговый сток» по литературным данным рассмотрен приток пресных вод к океану и айсберговый сток.
В разделе 2.3 «Стерические колебания уровня Мирового океана» показано, что за период с 1955 по 2003 год, подъем уровня из-за термического расширения в слое 0-700 м Мирового Океана составляет приблизительно 0.33 мм/год. Причем наибольшее влияние на термическое расширение оказывает нагрев Атлантического океана, который описывает примерно половину данной величины.
В разделе 2.4 «Сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XXстолетия» представлены последние оценки вклада различных факторов в колебания УМО по данным четвертого доклада МГЭИК за 1961-2003 и 1993-2003 гг. Показано, что точность многих оценок является весьма низкой. В последнее десятилетие отмечается резкое увеличение вклада всех факторов. В частности, в 4 раза возросли вклады стерической компоненты и гренландского ледяного щита, причем по абсолютной величине преобладает вклад термического расширения вод океана, который составляет половину роста УМО по данным альтиметрических наблюдений. Изменения массы льда Антарктиды за период 1961-2003 гг. определены с ошибкой в 3 раза превышающей саму оценку. Поэтому величина невязки (дисбаланса), равная 1,4 мм/год, превышает вклад любого из факторов. Еще больше дисбаланс за период 1993-2003 гг. (2 мм/год). Рост УМО эксперты объясняют плотностными изменениями толщи вод океана, колебаниями таяния как горных ледников, так и ледниковых щитов Антарктиды и Гренландии. Однако, таяние горных ледников, которое действительно весьма существенно, может давать непосредствен-
ный вклад в изменения УМО только с ледников, расположенных на островах в Северном Ледовитом океане и территории Аляски. Горные ледники, находящиеся в Европе, Азии, Африке и Ю. Америке, могут влиять на УМО, очевидно, только через приток речных вод к океану. Поэтому вряд ли имеет смысл непосредственно учитывать их вклад в тренд УМО (Малинин, 2009). Как известно, тренд в притоке материковых вод к МО отсутствует. Отсюда следует, что оценки вклада горных ледников в колебания УМО являются многократно преувеличенными. Кроме того, экспертами МГЭИК не учитывается вклад в УМО изменения вод литосферы.
Поэтому был выполнен сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XX в. (1980-1999 и 1990-2000 гг.) на основе уравнения пресноводного баланса океана (табл. 2). Установлено, что полученные невязки существенно меньше аналогичных ошибок, полученных экспертами МГЭИК на основе изменений содержания воды в криосфере и литосфере.
Таблица 2 - Оценка вклада разных факторов в изменения УМО за различные
промежутки времени XX столетия
Период времени 1980-1999 гг. 1990-2000 гг.
§ й По данным наблюдений 1,9 2,3
ш к о й Вертикальный влагообмен (Р-Е) 0,5 0,2
^О § Стерические колебания 0,2 0,6
я Р-Ег £ о 3 8 § я Твердый сток Антарктиды 0,6 1,1
Донное таяние шельфовых ледников 0,3 0,4
I & Суммарный сток Гренландии ОД 0,2
¡3 3 Невязка (дисбаланс) 0,2 -0,3
Глава 3 «Возможные изменения уровня Мирового океана в ближайшие десятилетия». В разделе 3.1 «Сценарии возможных изменений климата в XXI в. и их влияние на уровень Мирового океана» показано, что УМО можно рассматривать в качестве индикатора изменений климата.
В разделе 3.2 «Влияние межгодовых колебаний температуры воздуха на уровень океана» по литературным данным показаны связи между значениями УМО и аномалиями температуры воздуха, сделанные разными авторами в разное время. Выявленные связи подтверждаются высокими коэффициентами корреляции 0,80-0,95 при запаздывании температуры воздуха относительно УМО на 18, 19, 20, 30 лет. Очевидно, в этих случаях рост УМО может быть связан с «медленными» изменениями эвстатической и стерической компонент, т.е. с изменениями массы ледниковых щитов Антарктиды и Гренландии, а также повышением температуры воды в толще вод океана. В разделе 3.3 «Физико-статистическая модель прогноза межгодовых колебаний уровня Мирового океана» на основе аномалий приповерхностной температуры воздуха (ПТВ) из архива НасЮШТЗ за период с 1850 года и обновляемый каждый месяц по на-
стоящее время построена физико-статистическая модель прогноза межгодовых колебаний УМО.
С этой целью вначале был сформирован новый архив, который содержит только узлы, находящиеся непосредственно над акваторией океана и имеющие минимальное число пропусков. В результате из 2592 точек в новый архив вошло 196. Длина исходной матрицы значений температуры воздуха составила 107 (1901-2007 гг.) лет, а длина зависимой выборки п=99 (1901-1999 гг.). Большая их часть охватывает прибрежные районы вдоль всех материков земного шара и северную часть Атлантического океана, причем в южном полушарии их менее 20 %. В качестве функции отклика использовался рассчитанный УМО на основе футшточных (береговых) наблюдений за 1901-1999 гг..
Прогностическая регрессионная модель средних годовых значений УМО в зависимости от аномалий среднегодовой температуры воздуха над океаном имеет вид
йумо! ЛГго-т),--, ЛУш^Х (3)
где х — сдвиг в годах, ш - число исходных точек (ш=196).
Непосредственное применение к исходной матрице ПТВ метода включения переменных позволило осуществить расчет полного комплекса пошаговых моделей для сдвигов от т=0 до т=33 лет до последнего шага, на котором все включаемые в модель предикторы еще являются значимыми. Таким образом, общее число моделей составило более 500. Наиболее сложным моментом сравнения моделей является оценка точности их по стандартной ошибке с"у(х). Это связано с явной нестационарностью ряда УМО. Действительно, к последнему сдвигу длина ряда уменьшается от п=99 до «=66 значений. Из-за мощного тренда выборочные оценки средней и стандартного отклонения сильно изменяются. Так, для исходного ряда стандартное отклонение составляет Оумо = 53,05 мм/год, а при сдвиге т=33 лет стум0 = 38,29 мм/год. Для устранения влияния нестационарности функции отклика введен коэффициент эффективности модели у=сгУ(Х/оумо. Самый высокий коэффициент эффективности отмечается на 15 шаге для сдвига т=21 г., на котором он равен у=0.25. Поэтому примем его за основу при анализе прогностических оценок УМО.
Из сравнения моделей для различных сдвигов следует, что уже на первом шаге включаемая переменная описывает 52-70% дисперсии функции отклика. При этом почти на всех сдвигах (в 30 случаях из 34) первой переменной является точка с координатами 32,5° с.ш. и 137,5° в.д., расположенная вблизи побережья Японии.
В разделе 3.4 «Прогностические оценки уровня Мирового океана на ближайшие десятилетия» представлены результатам построения прогностических моделей для основного сдвига т=21 год. Приводится распределение статистических параметров до 19 шага включительно. Исходные данные ПТВ позволяют рассчитать прогностические значения УМО на 29 лет вперед, т.е. на 2000-2028
гг. На рисунке 4 представлены его значения для всех 19 шагов за независимый период 2000-2028 гг. Весьма важно, что разброс прогностических оценок уровня на каждый год относительно невелик и составляет в среднем около 4,5 см. При этом максимальный разброс оценок отмечается в 2003 г. (разброс 7 см), а минимальный - в 2011 г. (разброс 2 см). Принимая во внимание данное обстоятельство, очевидно, вряд ли имеет смысл придавать значение конкретным оценкам прогностических значений УМО и вполне достаточно ограничиться лишь величиной тренда (7>).
Рисунок 4 - Прогностические значения УМО за период 2000-2028 гг., рассчитанные по модели (4) при запаздывании уровня на 21 год для различных шагов от т=1 до «1=19. Жирная линия - осредненная кривая по 19 моделям УМО
На первом шаге Тг=2,6 мм/год, затем вплоть до 13 шага его величина Тг>3,0 мм/год с максимумом на 3 шаге (Тг=3,8 мм/год). После этого до 19 шага тренд медленно убывает до Тг-2,8 мм/год. Величина тренда по осредненной кривой для 19 моделей равна Тг=3,1 мм/год. Таким образом, результаты расчетов свидетельствуют, что в ближайшие годы следует ожидать сравнительно быстрого роста УМО.
В разделе 3.5 «Адаптивная модель прогноза сезонной изменчивости УМО по алътиметрическим данными обсуждаются их закономерности на основе данных университета Колорадо (США) за период 1993-2008 гг. При этом рассмотрены две прогностические модели сезонных колебаний УМО: детерминированная модель в виде суммы тренда и годовой гармоники УМО и адаптивная модель линейного роста с аддитивной сезонностью.
После приведения ряда УМО к стационарному виду путем исключения трендовой компоненты, был выполнен гармонический анализ, который позволил выделить значимую годовую гармонику. Остальные гармоники по критерию Стьюдента оказались незначимыми. Вклад годовой гармоники в дисперсию ряда УМО составляет 3,6 %, а ее амплитуда — 4,13 мм. Таким образом, на вклад остаточной компоненты в колебания УМО, характеризующей случайные колебания, приходится 4,2 %, а стандартное отклонение ее равно 3,1 мм. Отсюда следует, что вклад в сезонно-межгодовые колебания УМО космогеофизиче-ских и геолого-геодинамических факторов в современных климатических условиях ничтожно мал, они практически полностью обусловлены соответствующими изменениями климата.
Для адаптивного прогноза использовалась модель Тейла-Вейджа, описывающая сезонные явления с линейным ростом. Отличительной чертой таких моделей служит возможность «подстраиваться» под изменения временного ряда, придавая больший вес (более высокую информационную ценность) тем наблюдениям, которые ближе к текущему моменту прогнозирования. В результате возможно построение самокорректирующихся моделей, способных учитывать результат прогноза, сделанного на предыдущем шаге. Скорость реакции модели на изменения в динамике процесса характеризует так называемый параметр адаптации. Этот параметр адаптации должен быть выбран таким образом, чтобы обеспечивалось адекватное отображение тенденции при одновременной фильтрации случайных отклонений. Подбор оптимального параметра адаптации осуществлялся по минимизации среднеквадратической ошибки модели на ретроспективном статистическом материале.
С целью удобства использования исходных альтиметрических данных в практических расчетах предварительно было осуществлено их среднемесячное осреднение. Вся выборка за 1993-2008 гг. (N=168) была разделена на зависимую (1993-2004 гг.) и независимую (2005-2008 гг.) части. По зависимой части осуществлялось построение самой модели и ее верификация, а независимые данные использовались для проверки прогнозов УМО. Результаты тестирования и опытных прогнозов среднемесячных значений УМО на независимых данных показали возможность прогноза сезонного хода УМО на 4 года, т.е. на 48 месяцев. В результате прогноза по адаптивной модели коэффициент детерминации 0,75 и ошибка 2,62 мм, а для детерминированного прогноза Я2=0,73 и Су(Х)=2,86 мм.
Рассмотрены возможные изменения УМО в течение 2009-2012 гг. (рис. 5), рассчитанные по обеим моделям. В этом случае верификация моделей осуществлялась для выборки 1993-2008 гг. Нетрудно видеть, что обе модели дают согласованные оценки возможных изменений УМО. При этом даже скорость роста уровня примерно одинакова. Для детерминированной модели он соответствует зависимой выборке и составляет 3,26 мм/год, а для адаптивной модели -3,54 мм/год. Итак, к 2013 г. УМО может повыситься на 13-14 мм по сравнению с 2009 г.
47
42
37
32
27
; годы
2009
2010
2011
2012
2013
Рисунок 5 — Независимый прогноз УМО по адаптивной и детерминированной моделям и их линейные тренды
В заключении сформулированы основные результаты и выводы диссертационной работы:
1. Впервые построена статистическая модель УМО, обеспечивающая минимальную ошибку его восстановления по данным 5 береговых станций с достаточной для практических целей точностью и рассчитан межгодовой ход УМО с 1861 по 2005 г., т.е. за период инструментальных наблюдений. Показано, что полученные оценки УМО не уступают в точности результатам расчетов Church et al. и Jevreeva et al., в которых УМО вычислялся пространственным осреднением по данным более 1000 береговых станций. В межгодовом ходе УМО довольно отчетливо выделяется несколько периодов с различным характером изменения УМО, для которых локальные тренды существенно различаются. Это довольно быстрый рост УМО в период с 1861 по 1878 г. (Тг = 2,0 мм/год), затем в течение периода с 1879 по 1923 г. уровень практически не менялся, т.е. наблюдалась фаза относительного стояния уровня (Тг = 0,4 мм/год), затем вплоть до 2005 г. он опять начал быстро расти (Тг = 2,0 мм/год). В период интенсивного потепления климата (1976-2005 гг.) тренд УМО не изменялся (2,0 мм/год). В течение XX в. величина тренда составляла 1,8 мм/год.
2. Впервые выполнена оценка межгодовой изменчивости эвстатических факторов формирования УМО (испарение, осадки, эффективное испарение) на основе данных архива CDAS (Climate Data Assimilation System) за период с 1980 по 2006 г. Для временных рядов осадков и испарения характерен ярко выраженный положительный тренд, описывающий более 50 % дисперсии. Самый высокий тренд, как в осадках, так и в испарении отмечается над Индийским океаном. Величина тренда для осадков над Мировым океаном составляет Тг =
4.2 мм/год, а для испарения Тг = 3,6 мм/год, соответственно тренд в эффективном испарении оказывается отрицательным. Вследствие этого в течение 19802005 гг. происходил рост УМО со скоростью 0,6 мм/год, что соответствует 217 км3/год. Корреляция рядов испарения и осадков с глобальной температурой воздуха составляет соответственно г=0,82 и г=0,80.
3. Впервые построена диагностическая модель средних годовых значений УМО по данным о компонентах вертикального влагообмена над океанами. Оптимальная модель включает 8 точек по эффективному испарению Еэф, относительно равномерно расположенных на акватории Мирового океана, причем в зоне АЦТ находится сразу три точки. На основе независимых данных Еэф за 1998-2005 гг. рассчитаны значения УМО до 2006 г. Сопоставление трендов в УМО, рассчитанного по модели и по альтиметрическим данным за период 1993-2006 гг. показало их хорошее соответствие (2,7 и 2,9 мм/год). Кроме того, вычисленные по модели значения УМО в период 1981-2000 гг. согласуются с оценками УМО Church et al. и Jevreeva et al., коэффициенты корреляции между ними соответственно равны 0,85 и 0,91.
4. Выполнен сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XX в. (1980-1999 и 1990-2000 гг.) на основе уравнения пресноводного баланса океана. С учетом оценок трендов в эффективном испарении, материковом притоке вод, твердом и жидком стоке с Антарктиды и Гренландии, в таянии шельфовых ледников в Антарктиде и стерических колебаний уровня океана получены суммарные оценки тренда в УМО. Сравнение с оценками тренда УМО по альтиметрическим данным (1,9 мм/год за 1980-1999 и
2.3 мм/год за 1990-2000 гг.) показало, что невязка в первом случае составляет 0,2 мм/год, а во втором -0,3 мм/год, т.е. существенно меньше, чем аналогичные ошибки, полученные экспертами МГЭИК на основе изменений содержания воды в криофере и литосфере.
5. Установлено наличие высокой корреляции приповерхностной температуры воздуха (ПТВ) над океаном с УМО при его запаздывании относительно ПТВ на основе архива HadCrut3. В частности, в точке, расположенной вблизи побережья Японии (ф = 32,5°с.ш., X = 137,5° в.д.), корреляция температуры воздуха с УМО при его запаздывании относительно ПТВ не опускается ниже г=0,70 на всех сдвигах от 0 до 35 лет.
6. Впервые построена физико-статистическая модель прогноза УМО при запаздывании последнего относительно ПТВ. В результате расчета комплекса моделей для сдвигов от т=0 до т=33 лет установлено, что оптимальной в статистическом смысле является модель при запаздывании УМО на 21 год. Оценка прогностического тренда УМО до 2029 г. показала его небольшое уменьшение (Тг = 2,8 мм/год) по сравнению с трендом в начале XXI в. (Тг = 3,1 мм/год). Можно полагать, что в ближайшие годы повышение УМО будет продолжаться с некоторым замедлением темпов роста по сравнению с началом XXI в.
7. Впервые построена адаптивная модель прогноза сезонной изменчивости УМО на 2009-2012 гг. на основе альтиметрических данных. С этой целью ис-
пользована модель Тейла-Вейджа, описывающая сезонные явления с линейным ростом. Одновременно для сравнения рассматривалась «детерминированная» модель, состоящая из линейного тренда и годовой гармоники, описывающая 94 % дисперсии исходного ряда УМО. Показано, что обе модели дают согласованные оценки возможных изменений УМО. Для детерминированной модели скорость роста УМО соответствует зависимой выборке и составляет 3,23 мм/год, а для адаптивной модели - 3,54 мм/год. Установлено, что адаптивная модель обладает более высокой точностью, так как дополнительно учитывает определенную часть дисперсии, сосредоточенной в случайной компоненте разложения. В результате расчетов выявлено, что по сравнению с 2009 г. УМО к 2013 г. может повыситься на 13-14 мм.
Список публикаций по теме диссертации
Статьи в изданиях, входящих в Перечень ВАК
1. В.Н. Малинин, О.И. Шевчук. Эвстатические колебания уровня Мирового океана в современных климатических условиях И Изв. Рус. геогр. общества. 2008. Т.140. Вып.4. С. 20-30.
2. Клиге Р.К., Малинин В.Н., Шевчук О.И. Колебания уровня Мирового океана в XX столетии // Вести. Моск. ун-та. Сер. 5. География. 2009. №1. С. 7-14.
Материалы конференций
3. О.И. Шевчук. О проблеме мониторинга уровня Мирового океана // Сборник трудов международной школы-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда». - СПб.: Гранд,
2005. - С. 235-239.
4. О.И. Шевчук. К анализу трендов уровня Мирового океана // Материалы итоговой сессии Ученого совета 25-26 января 2006 года. - СПб.: изд-во РГТМУ,
2006. - С. 63-64.
5. В.Н. Малинин, О.И. Шевчук. О прогнозе колебаний уровня Мирового океана на основе физико-статистической модели // Географическое образование и наука в России: история и современное состояние. Материалы международной научно-практической конференции, Санкт-Петербург, 3-4 декабря 2008 г. -СПб.: изд-во СПбГУ, 2009. - (в печати).
6. Малинин В.Н., Гордеева С.М., Шевчук О.И. Статистический анализ изменчивости уровня Мирового океана по альтиметическим данным // География и геоэкология на современном этапе взаимодействия природы и общества: мат-лы Всерос. научн. конф. «СЕЛИВЕРСТОВСКИЕ ЧТЕНИЯ», Санкт-Петербург, 19-20 ноября 2009 г. - СПб.: СПбГУ, ВВМ, 2009. - С.56-61.
Статьи в периодических изданиях
7. Малинин В.Н., Гордеева С.М., Шевчук О.И. Изменчивость уровня Мирового океана за последние 140 лет // Учен, записки РГГМУ. 2007. Вып. 4. С. 125-131.
8. Малинин В.Н., Шевчук О.И. О возможных изменениях уровня мирового океана в ближайшие десятилетия // Общество. Среда. Развитие. 2009. №2. С. 172-180.
9. В.Н. Малинин, О.И. Шевчук. О прогнозе сезонного хода уровня Мирового океана // Фундаментальная и прикладная гидрофизика. 2009. № 3 (5). С. 25-32.
Отпечатано с готового оригинал-макета в ЦНИТ «АСТЕРИОН» Заказ № 322. Подписано в печать 20.11.2009 г. Бумага офсетная. Формат 60х84'/1б.0бъем 1,25 п. л. Тираж 100 экз. Санкт-Петербург, 191015, а/я 83, (812) 275-73-00,970-35-70 asterion@asterion.ru
Содержание диссертации, кандидата географических наук, Шевчук, Олег Игоревич
Введение
Глава 1 Колебания уровня Мирового океана в современных климата- 11 ческих условиях по данным наблюдений
1.1 Краткий обзор исследований колебаний уровня океана на ос- 11 нове футшточных наблюдений
1.2 Межгодовая изменчивость уровня Мирового океана по дан- 17 ным футшточных наблюдений
1.3 Изменчивость уровня Мирового океана по альтиметрическим 34 данным
Глава 2 Факторы формирования уровня Мирового океана и их ана
2.1 Эвстатические изиенения уровня океана как интегральный 40 индикатор глобального водообмена
2.2 Эвстатические колебания уровня Мирового океана.
2.2.1 Вертикальный влагообмен в системе океан-атмосфера.
2.2.2 Приток пресных вод к океану. Айсберговый сток.
2.3 Стерические колебания уровня Мирового океана
2.4. Сравнительный анализ вклада разных факторов в рост уров- 67 ня Мирового океана в XX столетии
Глава 3 Возможные изменения уровня Мирового океана на ближай- 70 шую перспективу
3.1 Сценарии возможных изменений климата в XXI веке и их 70 влияние на уровень Мирового океана.
3.2. Влияние межгодовых колебаний температуры воздуха на 73 уровень океана
3.3. Физико-статистическая модель прогноза межгодовых коле- 77 баний уровня Мирового океана
3.4. Прогностические оценки уровня Мирового океана в бли- 86 жайшие десятилетия
3.5. Адаптивная модель прогноза сезонной изменчивости УМО 89 по альтиметрическим данным
Введение Диссертация по наукам о земле, на тему "Межгодовая изменчивость уровня мирового океана в современных климатических условиях"
Актуальность темы. Уровень Мирового океана является важнейшей характеристикой, определяющей как гидродинамические процессы в самом океане, так и характер взаимодействия океана и атмосферы.
Изменения уровня Мирового океана тесно связаны с современными колебаниями климата, которые в свою очередь определяют испарение, осадки, стерические колебания. Эвстатические колебания можно рассматривать как индикатор глобального водообмена, перераспределяющего воды гидросферы между отдельными оболочками Земли и, прежде всего, между водами литосферы и криосферы.
Как следует из докладов Межправительственной группы экспертов по изменениям климата (МГЭИК) [61, 63], в соответствии с некоторыми сценариями изменений климата, основанными на результатах его численного моделирования, уровень Мирового океана может повыситься на 10-3 0'см к 2030 г. и в пределах 30-100 см (наиболее вероятно 65 см) к концу столетия. Естественно, если такое развитие изменений климата станет реальностью, то это грозит катастрофическим ущербом для инфраструктуры прибрежных территорий, где проживает около миллиарда жителей Земли. Возможно подс* топление крупнейших городов мира: Лондона, Нью-Йорка, Токио, С-Петербурга и др. Повышение уровня океана опасно для человечества не только сокращением доступной для жизни территории суши с затоплением множества заселенных территорий. Прибрежные пресноводные водоносные слои могут быть засолены при повышении уровня соленых грунтовых вод. Перемещение фронта соленых вод вверх по устьям рек негативно скажется на работе сооружений по забору пресной воды, расположенных вверх по течению, на рыбный промысел в солоноватых водах и сельское хозяйство. В зависимости от разных причин темпы и масштабы подъема уровня моря могут значительно превзойти способность антропогенных и естественных систем к адаптации без существенных последствий.
Поэтому задача изучения долговременных колебаний уровня Мирового океана и особенно построения систем прогноза уровня Мирового океана в системе океан атмосфера представляется чрезвычайно актуальной. Эксперты МГЭИК отнесли данный вопрос к числу важнейших проблем современной гидрометеорологии.
Цель работы: изучение межгодовых колебаний уровня Мирового океана и эвстатических факторов его формирования.
Для достижения указанной цели решались следующие задачи:
1. Анализ межгодовых колебаний уровня океана за период инструментальных наблюдений.
2. Выявление станций с многолетними наблюдениями за морским уровнем, обеспечивающих минимальную ошибку восстановления уровня Мирового океана (УМО) с достаточной для практических целей точностью и построение статистической модели оценки межгодовых колебаний УМО.
3. Анализ межгодовой изменчивости эвстатических факторов (испарения и осадков) формирования уровня Мирового океана. Построение диагностической статистической модели оценки межгодовых колебаний УМО на основе эвстатических факторов.
4. Сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XX столетия.
5. Оценка связи колебаний уровня Мирового океана с приповерхностной температурой воздуха над океаном.
6. Построение физико-статистической модели прогноза межгодовых колебаний уровня Мирового океана на ближайшие десятилетия.
7. Построение адаптивной модели прогноза сезонной изменчивости уровня Мирового океана по альтиметрическим данным.
Материалы и методы. Основой работы являлись несколько глобальных архивов гидрометеорологических данных находящихся в свободном доступе в Интернет:
- футшточные данные об уровне Мирового океана полученные из архива PSMSL (Permanent Service for Mean Sea Level);
- альтиметрические данные центра по хранению, верификации и интерпретации спутниковой океанографической информации (AVISO) и Колорадского центра астродинамических исследований (CCAR) при Колорадском университете со спутников Topex/Poseidon и Jason-1, с пространственным разрешением 1/3°х1/3°;
- среднемесячные значения осадков и испарения в узлах первичной широтно-долготной сетки 1.875°х 1.875° из глобального гидрометеорологического архива CD AS (Climate Data Assimilation System);
- аномалии поверхностной температуры воздуха из архива HadCrut3 г с пространственным разрешением 5°х5°;
Для выполнения поставленных задач использовался комплекс методов многомерного статистического анализа (MCA), а также теория проверки статистических гипотез и анализ временных рядов.
Научная новизна:
- Впервые построена статистическая модель, обеспечивающая минимальную ошибку восстановления межгодовых колебаний уровня Мирового' океана по данным 5' береговых станций с достаточной для практических целей точностью.
- Установлены закономерности межгодовых колебаний УМО за период инструментальных наблюдений.
- Впервые выявлены особенности межгодовой изменчивости испарения, осадков и эффективного испарения над отдельными океанами и Мировым океаном в целом в период интенсивного глобального потепления климата,
- Установлена степень взаимосвязи колебаний УМО с температурой воздуха над его акваторией.
- Впервые построена физико-статистическая модель прогноза УМО на ближайшие десятилетия на основе данных о температуре воздуха над его акваторией.
- Впервые построена адаптивная модель прогноза сезонной изменчивости УМО по альтиметрическим данным.
На защиту выносятся:
1. Закономерности межгодовых колебаний уровня Мирового океана за период инструментальных наблюдений.
2. Особенности межгодовой изменчивости эвстатических факторов (осадки, испарение, эффективное испарение) над океанами.
3. Результаты прогнозирования уровня Мирового океана по данным о температуре воздуха над его акваторией на ближайшие десятилетия.
4. Адаптивная модель прогноза сезонной изменчивости уровня Мирового океана по альтиметрическим данным.
Практическая значимость работы. Полученные в работе результаты могут быть использованы при решении задач мониторинга УМО и моделирования климата, а также при разработке долгосрочных прогнозов межгодовых колебаний и сезонного хода УМО.
Личный вклад автора. Подготовка архивов информации, проведение статистических расчетов, выполнение необходимых графических построений, обобщение и анализ результатов.
Апробация работы и публикации. Основные положения и отдельные результаты обсуждались на итоговых сессиях Ученого Совета РГГМУ, 2005 и 2006; на международной школе-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда», РГГМУ, 2005; научном семинаре «Проблемы изменения климата в полярных широтах» в рамках проекта BRIDGE «Keeping it Cool», Университет г. Шеффилда, Великобритания, 2007; международной научно-практической конференции «Географическое образование и наука в России: история и современное состояние», СПБГУ, 2008; всероссийской научной конференции «Селиверстовские чтения», 2009.
По теме диссертации опубликовано 9 печатных работ, из них 2 статьи в рецензируемых журналах (по перечню ВАК):
Статьи в изданиях, входящих в Перечень ВАК
1. В.Н. Малинин, О.И. Шевчук. Эвстатические колебания уровня Мирового океана в современных климатических условиях // Изв. Рус. геогр. общества. 2008. Т. 140. Вып.4. С. 20-30.
2. Клиге Р.К., Малинин В.Н., Шевчук О.И. Колебания уровня Мирово- * го океана в XX столетии // Вестн. Моск. ун-та. Сер. 5. География. 2009. №1. С. 7-14.
Материалы конференций
3. О.И. Шевчук. О проблеме мониторинга уровня Мирового океана // Сборник трудов международной школы-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда». - СПб.: Гранд, 2005. - С. 235-239.
4. О.И. Шевчук. К анализу трендов уровня Мирового океана // Материалы итоговой сессии Ученого совета 25-26 января 2006 года. - СПб.: изд-во РГГМУ, 2006. - С. 63-64.
5. В.Н. Малинин, О.И. Шевчук. О прогнозе колебаний уровня Мирового океана на основе физико-статистической модели // Географическое образование и наука в России: история и современное состояние. Материалы международной научно-практической конференции, Санкт-Петербург, 3-4 декабря 2008 г. - СПб.: изд-во СПбГУ, 2009. - (в печати).
6. Малинин В.Н., Гордеева С.М., Шевчук О.И. Статистический анализ изменчивости уровня Мирового океана по альтиметическим данным // География и геоэкология на современном этапе взаимодействия природы и общества: мат-лы Всерос. научн. конф. «СЕЛИВЕРСТОВСКИЕ ЧТЕНИЯ», Санкт-Петербург, 19-20 ноября 2009 г. - СПб.: СПбГУ, ВВМ, 2009. -С.56-61.
Статьи в периодических изданиях
7. Малинин В.Н., Гордеева С.М., Шевчук О.И. Изменчивость уровня Мирового океана за последние 140 лет // Учен, записки РГГМУ. 2007. Вып. 4. С. 125-131. ^
8. Малинин В.Н., Шевчук О.И. О возможных изменениях уровня мирового океана в ближайшие десятилетия // Общество. Среда. Развитие. 2009. №2. С. 172-180.
9. В.Н. Малинин, О.И. Шевчук. О прогнозе сезонного хода уровня Мирового океана // Фундаментальная и прикладная гидрофизика. 2009. № 3 (5). С. 25-32.
Структура и объем работы. Диссертация состоит из введения, 3 глав, заключения и библиографического списка использованной литературы. Объем диссертации составляет 113 страниц, включая 23 таблицы и 34 рисунка. Список литературы составляет 84 наименования, в том числе 42 на иностранных языках.
Заключение Диссертация по теме "Океанология", Шевчук, Олег Игоревич
Заключение
В результате выполненной работы могут быть сделаны следующие основные выводы:
1. Впервые построена статистическая модель УМО, обеспечивающая минимальную ошибку его восстановления по данным 5 береговых станций с достаточной для практических целей точностью и рассчитан межгодовой ход УМО с 1861 по 2005 г., т.е. за период инструментальных наблюдений. Показано, что полученные оценки УМО не уступают в точности результатам расчетов Church et al. и Jevreeva et al., в которых УМО вычислялся пространственным осреднением по данным более 1000 береговых станций. В межгодовом ходе УМО довольно отчетливо выделяется несколько периодов с различным характером изменения УМО, для которых локальные тренды существенно различаются. Это довольно быстрый рост УМО в период с 1861 по 1878 г. (Тг = 2,0 мм/год), затем в течение периода с 1879 по 1923 г. уровень практически не менялся, т.е. наблюдалась фаза относительного стояния уровня (Тг = 0,4 мм/год), затем вплоть до 2005 г. он опять начал быстро расти (Тг = 2,0 мм/год). В период интенсивного потепления климата (1976-2005 гг.) тренд УМО не изменялся (2,0 мм/год). В течение XX в. величина тренда составляла 1,8 мм/год.
2. Впервые выполнена оценка межгодовой изменчивости эвстатических факторов формирования УМО (испарение, осадки, эффективное испарение) на основе данных архива CDAS (Climate Data Assimilation System) за период с 1980 по 2006 г. Для временных рядов осадков и испарения характерен ярко выраженный положительный тренд, описывающий более 50 % дисперсии. Самый высокий тренд, как в осадках, так и в испарении отмечается над Индийским океаном. Величина тренда для осадков над Мировым океаном составляет Тг = 4,2 мм/год, а для испарения Тг = 3,6 мм/год, соответственно тренд в эффективном испарении оказывается отрицательным. Вследствие этого в течение 1980-2005 гг. происходил рост УМО со скоростью 0,6 мм/год, что соответствует 217 км /год. Корреляция рядов испарения и осадков с глобальной температурой воздуха составляет соответственно г=0,82 и г=0,80.
3. Впервые построена диагностическая модель средних годовых значений УМО по данным о компонентах вертикального влагообмена над океанами. Оптимальная модель включает 8 точек по эффективному испарению ЕЭф, относительно равномерно расположенных на акватории Мирового океана, причем в зоне АЦТ находится сразу три точки. На основе независимых данных ЕЭф за 1998-2005 гг. рассчитаны значения УМО до 2006 г. Сопоставление трендов в УМО, рассчитанного по модели и по альтиметрическим данным за период 1993-2006 гг. показало их хорошее соответствие (2,7 и 2,9 мм/год). Кроме "того, вычисленные по модели значения УМО в период 1981-2000 гг. согласуются с оценками УМО Church et al. и Jevreeva et al., коэффициенты корреляции между ними соответственно равны 0,85 и 0,91.
4. Выполнен сравнительный анализ вклада разных факторов в рост уровня Мирового океана в конце XX в. (1980-1999 и 1990-2000 гг.) на основе уравнения пресноводного баланса океана. С учетом оценок трендов в эффективном испарении, материковом притоке вод, твердом и жидком стоке с Антарктиды и Гренландии, в таянии шельфовых ледников в Антарктиде и сте-рических колебаний уровня океана получены суммарные оценки тренда в УМО. Сравнение с оценками тренда УМО по альтиметрическим данным (1,9 мм/год за 1980-1999 и 2,3 мм/год за 1990-2000 гг.) показало, что невязка в первом случае составляет 0,2 мм/год, а во втором -0,3 мм/год, т.е. существенно меньше, чем аналогичные ошибки, полученные экспертами МГЭИК на основе изменений содержания воды в криофере и литосфере.
5. Установлено наличие высокой корреляции приповерхностной температуры воздуха (ПТВ) над океаном с УМО при его запаздывании относительно
ПТВ на основе архива НасЮпйЗ. В частности, в точке, расположенной вблизи побережья Японии (ф = 32,5°с.ш., Х= 137,5° в.д.), корреляция температуры воздуха с УМО при его запаздывании относительно ПТВ не опускается ниже 1-0,70 на всех сдвигах от 0 до 35 лет.
6. Впервые построена физико-статистическая модель прогноза УМО при запаздывании последнего относительно ПТВ. В результате расчета комплекса моделей для сдвигов от т=0 до т=33 лет установлено, что оптимальной в статистическом смысле является модель при запаздывании УМО на 21 год. Оценка прогностического тренда УМО до 2028 г. показала его небольшое уменьшение (Тг = 2,6-2,8 мм/год) по сравнению с трендом в начале XXI в. (Тг = 3,1 мм/год). Можно полагать, что в ближайшие годы повышение УМО будет продолжаться с некоторым замедлением темпов роста по сравнению с началом XXI в.
7. Впервые построена адаптивная модель прогноза сезонной изменчивости УМО на 2009-2012 гг. на основе альтиметрических данных. С этой целью использована модель Тейла-Вейджа, описывающая сезонные явления с линейным ростом. Одновременно для сравнения рассматривалась «детерминированная» модель, состоящая из линейного тренда и годовой гармоники, описывающая 94 % дисперсии исходного ряда УМО. Показано, что обе модели дают согласованные оценки возможных изменений УМО. Для детерминированной модели скорость роста УМО соответствует зависимой выборке и составляет 3,23 мм/год, а для адаптивной модели - 3,54 мм/год. Установлено, что адаптивная модель обладает более высокой точностью, так как дополнительно учитывает определенную часть дисперсии, сосредоточенной в случайной компоненте разложения. В результате расчетов выявлено, что по сравнению с 2009 г. УМО к 2013 г. может повыситься на 13-14 мм.
Библиография Диссертация по наукам о земле, кандидата географических наук, Шевчук, Олег Игоревич, Санкт-Петербург
1. Айвазян С.А., Мхитарян B.C. Прикладная статистика и основыэконометрики. -М.: ЮНИТИ, 1998. 1022 с.
2. Афифи А., Эйзен С. Статистический анализ: подход с использованием
3. ЭВМ. М.: Мир, 1982. - 488 с.
4. Буторин Н. В. Вековые изменения среднего уровня Атлантического океанаи их связь с циркуляцией атмосферы. М.- JL: Изд-во АН СССР, 1960. -334 с.
5. Вайновский П. А. Методы обработки и анализа океанологическойинформации. 4.2. Многомерный анализ. СПб.: Изд. РГМИ, 1992. - 96 с.
6. Гордеева С.М., Малинин В.Н. О межгодовой изменчивости компонентовводного баланса Балтийского моря // Современные проблемы гидрометеорологии. Сборник научных трудов РГГМУ. Вып. 122. СПб: Изд. РГГМУ, 1999. - С. 64-80.
7. Дженнрич Р.И. Пошаговая регрессия //Статистические методы для ЭВМ /
8. Под ред. К. Энслейна и др. М.: Наука, 1986. - С. 77-93.
9. Добровольский С.Г. Климатические изменения в системе «гидросфераатмосфера». М.: Геос, 2002. - 232 с.
10. Добровольский С.Г. Изменения климата и составляющих водного баланса
11. Раткович Д.Я. Актуальные проблемы водообеспечения. М.: Наука, 2003.-С. 119-130.
12. Догановский A.M., В.Н. Малинин. Гидросфера Земли. СПб.:
13. Гидрометеоиздат, 2004. 630 с.
14. Дрейпер Н., Смит Г. Прикладной регрессионный анализ: В 2 кн. М.:
15. Финансы и статистика, 1986. — 366 е.; 1987. — 351 с.
16. Дуброва Т.А. Прогнозирование социально-экономических процессов. —
17. М.: Маркет ДС, 2007. 192 с.
18. Дуванин А. И., Калинин Г. П., Книге Р. К. О многолетних колебанияхуровня океанов, некоторых морей и озер // Вестн. Моск. ун-та. Сер. 5. География. 1975. № 6. С. 3-15.
19. Зайцева И. С. Антропогенные воздействия на водные ресурсыконтинентов. М.: Научный Мир, 2000. - С. 183-194
20. Калинин Г. П., Клиге Р. К. К вопросу о вековых колебаниях уровня
21. Мирового океана. / В сб.: Формирование ресурсов вод суши / Под ред. Г.П. Калинина. -М.: «Наука», 1972. 132 с.
22. Клиге Р. К. Изменения глобального водообмена М.: Наука. 1985 - 247 с.
23. Клиге Р. К., Захаров В. Г. Изменения снежно-ледового режима
24. Антарктиды // В кн.: Современные глобальные изменения природной среды: в 2 т. /Отв. ред.: Чл.-кор. РАН Н.С. Касимов, проф. Р.К. Клиге. -М.: Научный Мир, 2006. Т. 2. С. 577-606.
25. Клиге Р.К., Малинин В.Н., Шевчук О.И. Колебания уровня Мировогоокеана в XX столетии // Вестн. Моск. ун-та. Сер. 5. География. 2009. №1. С.7-14.
26. Котляков В. М. Гляциология Антарктиды. / Избранные сочинения: В 6 кн.-М.: Наука, 2000. Кн. 1.-430 с.
27. Лебедев С.А., Костяной А.Г. Спутниковая альтиметрия Каспийскогоморя. М.: Издательский центр «МОРЕ» Международного института океана, 2005. - 366 с.
28. Лукашин Ю.П. Адаптивные методы краткосрочного прогнозированиявременных рядов. М.: Финансы и статистика. 2003. - 415 с.
29. Максимов И.В., Саруханян Э.И., Смирнов Н.П. Океан и космос. Л.:
30. Гидрометеоиздат, 1970.-215 с.
31. Малинин, В.Н. Влагообмен в системе океан-атмосфера. СПб,
32. Гидрометеоиздат, 1994. 198 с.
33. Малинин В.Н. О современном состоянии проблемы изменчивости водгидросферы. Учен. Записки РГГМУ. Вып.1. - СПб.: Изд. РГГМУ, 2005. — С.54-75.
34. Малинин В.Н. Статистические методы анализа гидрометеорологическойинформации: Учебник. СПб., Изд. РГГМУ, 2008. -408с.
35. Малинин В.Н. Изменчивость глобального водообмена в условияхменяющегося климата // Водные ресурсы. 2009. Т. 36. №1. С. 1-14.
36. Малинин В.Н., Гордеева С.М. Физико-статистический метод прогнозаокеанологических характеристик. Мурманск: Изд. ПИНРО, 2003 .- 164 с.
37. Малинин В.Н., Гордеева С.М., Шевчук О.И. Изменчивость уровня
38. Мирового океана за последние 140 лет // Учен, записки РГГМУ. 2007. Вып. 4. С. 125-131.
39. Малинин В.Н., Радикевич В.М. Гордеева С.М., Куликова Л.А.
40. Изменчивость вихревой активности атмосферы над Северной Атлантикой. СПб, Гидрометеоиздат, 2003. - 171 с.
41. Малинин В.Н., Шевчук О.И. Эвстатические колебания уровня Мировогоокеана в современных климатических условиях. // Изв. Рус. геогр. общества. 2008. Т. 140. Вып.4. С. 20-30.
42. Малинин В.Н., Шевчук О.И. О прогнозе сезонного хода уровня Мирового океана // Фундаментальная и прикладная гидрофизика. 2009. № 3 (5). С. 25-32.
43. Мелешко В.П., Катцов В.М., Мирвис В.М., Говоркова В.А., Павлова Т.В.
44. Климат России в XXI веке. Часть 1. Новые свидетельства антропогенного изменения климата и современные возможности его расчета // Метеорология и гидрология. 2008. № 6. С.5-19.
45. Мировой водный баланс и водные ресурсы Земли. JL: Гидрометеоиздат,1974.-638 с.
46. Оценочный доклад об изменениях климата и их последствиях натерритории Российской Федерации. Техническое резюме. М.: Изд. Росгидромет, 2008. - 89 с. - Режим доступа: http://climate2008.igce.ru/ v2008/pdf7resumeteh.pdf)
47. Провоторов П. П. Стерические колебания уровня моря. // Колебанияуровня в морях. Сборник научных трудов РГГМУ. СПб.: Изд. РГГМУ, 2003.-С. 129-138.
48. Романов A.A. Ледовые условия плавания в Южном океане. СПб.:
49. Гидрометеоиздат, 1996. 119 с. - URL: http://www.aari.nw.ru/gdsidb/ pub/WMO783rus.pdf
50. Степанов В.Н. Океаносфера. М.: Мысль, 1983. - 270 с.
51. Фукс В.Р. Уровень Мирового океана как индикатор глобальногопотепления. II География и современность. 2005. Вып. 10. С. 73-93.
52. Шевчук О.И. О проблеме мониторинга уровня Мирового океана //
53. Сборник трудов международной школы-конференции молодых ученых, аспирантов и студентов «Изменение климата и окружающая среда». -СПб.: Гранд, 2005. С. 235-239.
54. Шнитников А. В. Внутривековая изменчивость компонентов общейувлажненности. JL: Наука, 1969. - 246 с.
55. Alley R.B., Clark P.U., Huybrechts P., Joughin I. Ice-Sheet and Sea-Level
56. Change // Science. 2005. 21 October. Vol. 310. P. 456-460.
57. Antonov J. I., Levitus S., Boyer T.B. Thermosteric Sea Level Rise, 1955-2003
58. Geophysical Research Letters. 2005. Vol. 32. L12602, doi: 10.1029/2005GL023112.
59. Ballantyne J.A. Multidecadal Study of the Number of Atlantic Icebergs Using
60. Scatterometr Data. Brigham: Brigham Young Universit. 2002. 232 p.
61. Barnett T.P. The Estimation of Global Sea Level Change: A Problem of
62. Uniqueness // Journal of Geophysical Research. 1984. Vol. 89. C.5, P. 79807988.
63. Brohan, P., Kennedy J.J., Harris I., Tett S.F.B., Jones P.D. Uncertainty
64. Estimates in Regional and Global Observed Temperature Changes: A New Dataset from 1850 // Journal of Geophysical Research. 2006. Vol. Ill, D12106, doi:10.1029/2005JD006548.
65. Brown R. G. Smoothing, forecasting and prediction of discrete time series,
66. New Jersey, Prentice-Hall. 1963. 125 p
67. Cazenave A., Nerem R.S. Present-day Sea Level Change: Observations and
68. Causes // Reviews of Geophysics. 2004. Vol. 42. № 3. P. 1-20.
69. Church J.A., White N.J. A 20th Century Acceleration in Global Sea-level Rise
70. Geophysical Research Letters. 2006. Vol. 33. № 1. L01602. doi: 10.1029/2005GL024826. URL: http://www.pol.ac.uk/psmsl/author archive/churchwhite/GRLChurchWhite2006024826.pdf
71. Climate Change 2007: The Physical Science Basis. Intergovernmental Panel on
72. Climate Change WGI Fourth Assessment Report. Geneva, 2007. 18 p.
73. Dietrich C. Ozeanographisch-meteorologische Einflüsse and Wasserstandsanderungen des Meersam Beispiel der Regelbeo-Bachtungen von Esbjerg. Die Küste, Bd2. H. 2, 1954.
74. Douglas B.C. Global Sea Rise: A Redetermination // Survey in Geophysics. 1997. Vol. 18. P. 279-292.
75. Fairbridge R.W., Krebs O.A. Sea Level and the Southern Oscillation // Journalof Geophysical Research. 1962. V. 6. P. 32-545
76. Fairbridge, R.W. Mean Sea Level Related to Solar Radiation During the Last20,000 Years. // Changes of Climate. 1963. № 20. P. 229-241
77. Gornits V., Lebedeff L., Hansen J. Global Sea Level Trend in the Past Century
78. Science. 1982. Vol. 215. P. 1611-1614.
79. Gutenberg B. Changes in Sea Level Postglacial Uplift and Mobility of the
80. Earth Interior// Bull. Geol. Soc. Amer. 1941. Vol. 52. № 5. P. 721- 772.
81. Hanna H., Huybrechts P., Janssens I., Cappelen L., Steffen K., Stephens A.
82. Runoff and Mass Balance of the Greenland Ice Sheet: 1958-2003 // Journal Geophysical Research. 2005. V. 110. D13108. doi: 10.1029/2004JD005641
83. Hansen J. E., R. Ruedy, Mki. Sato, M. Imhoff, W. Lawrence, D. Easterling, T.
84. Peterson, and T. Karl. A Closer Look at United States and Global Surface Temperature Change // Journal of Geophysical Research. 2001. Vol. 106. P. 23947-23963. doi: 10.1029/2001JD000354
85. Intergovernmental Panel on Climate Change (IPCC) Working Group I Third
86. Assessment Report. Chapter 11. Changes in Sea Level / Ed. by Douglas B.C., Ramirez A. 2000. 76 p.
87. IPCC. Climate Change 2001: The Scientific Basis. Contribution of Working
88. Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change // Ed. by Houghton J.T. Cambridge: Cambridge University Press, 2001. 881 p.
89. IPCC. Climate Change 2007: Synthesis Report. Contribution of Working
90. Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change // Ed. by Pachauri R.K, Reisinger A. Geneva, 2007. 104 p.
91. IPCC. Climate Change 2007: The Physical Science Basis. Intergovernmental
92. Panel on Climate Change Fourth Assessment Report Climate Change 2007 // Ed. by Bernstein L. Cambridge: Cambridge University Press, 2007. 940 p.
93. IRI/LDEO Climate Data Library / Lamont-Doherty Earth Observatory. New York: Columbia University Press. URL:http://iridl.ldeo.columbia.edu
94. Ishii M., Kimoto M., Sakamoto K., Iwasaki S.I. Steric Sea Level Changes
95. Estimated from Historical Ocean Subsurface Temperature and Salinity Analyses//Journal of Oceanography, 2006. Vol. 62 (No. 2). P.155-170.
96. Jevrejeva S., Grinsted A., Moore J.C., Holgate S. Nonlinear Trends and
97. Multiyear Cycles in Sea Level Records // Journal of Geophysical Research. 2006. Vol. 111. C09012. doi: 10.1029/2005JC003229
98. Kalnay E., Kanamitsu M., Kistler R., Collins W., Deaven D., Gandin L., Iredell
99. Kuenen P. Marine Geology. New York: John Wiley and Sons, 1950. p. 568.
100. Leuliette E. W., Nerem R.S., Mitchum G.T. Calibration of TOPEX/Poseidonand Jason Altimeter Data to Construct a Continuous Record of Mean Sea Level Change // Marine Geodesy. 2004. Vol. 27(1-2). P. 79-94.
101. Lisitzen E. Bull Inf. C.O.E.C, 1958. Vol. 10. № 5.
102. Meade R. H., Emery K. O. Sea Level as Affected by River Runoff // eastern
103. United States: Science. 1971. Vol. 173. № 3995. P. 425-427. 73 Nerem R.S., Chambers D.P., Leuliette E.W., Mitchum G.T. Satellite
104. Measurements of Sea Level Change: Where have been and here are we going // ESA Portal. Earthnet online. URL: http://earth.esa.int/cgi-bin/confaltl5y.pl? abstract=1092J
105. WOCE). URL: http://woce.nodc.noaa.gov/wocev3/wocedata2/sleveldm/ psmsl /psmsl.htm
106. Sea Level Change // University of Colorado at Boulder. URL:http://sealevel.colorado.edu
107. Smith T.M., Reynolds R.W. A Global Merged Land and Sea Surface
108. Temperature Reconstruction Based on Historical Observations (1880-1997) // Journal of Climate. 2005. Vol. 18, P. 2021-2036.
109. Theil H., Wage S. Some Observations on Adaptive Forecasting // Management
110. Science. 1964. Vol. 10. P. 198—206
111. Valentin H. Die Kusten der Erde. Petermanns Geogr. Mitt. Erganzungsheff.
112. Berlin: Erganzungsheft, Justus Perthes, Gotha. 246 p.
113. Trupin A., Wahr J. Spectroscopic Analysis of Global Tide Gauge Sea Level
114. Data // Geophysical Journal International. 1996. Vol. 100. P. 441-453.
115. Tushingham A.M., Peltier W.R. Ice-3G: A New Global Model of Late
116. Pleistocene Deglaciation Based upon Geophysical Predictions of Post-glacial Relative Sea Level Change // Journal of Geophysical Research. 1991. Vol. 96. P. 4497-4523.
117. Unal Y. S., Ghil M. Interannual an Interdecadal Oscillation Patterns in Sea Level // Climate Dynamics Vol. 11. P. 255- 278.
118. Walsh J.E., Portis D.H. Variations of Precipitations and Evaporation over the
119. North Atlantic Ocean, 1958-1997 //-Journal of Geophysical Research. 1999. Vol. 104. P.16613-16631.
120. Willis J. K., Roemmich D., Cornuelle B. Interannual Variability in Upperocean Heat Content, Temperature and Thermosteric Expansion on Global Scales // Journal of Geophysical Research. 2004. Vol. 109. C12036, doi: 10.1029/2003JC002260.
121. World Water Resources of the Earth on the Border XXI Century. London:
122. Cambridge Press, 2002. 465 p.
- Шевчук, Олег Игоревич
- кандидата географических наук
- Санкт-Петербург, 2009
- ВАК 25.00.28
- Статистический анализ межгодовой изменчивости теплосодержания и стерических колебаний уровня Мирового океана
- Структура и изменчивость океанических фронтальных зон: анализ глобальной спутниковой информации
- Структура и динамика промежуточных водных масс Антарктического происхождения в южной части Тихого океана
- Изменчивость гидрометеорологических полей Атлантико-Евразийского региона под влиянием Североатлантического и Южного колебаний
- Крупномасштабный термохалинный режим вод и формирование зон биологической продуктивности Атлантического океана