Главная Другое
Экономика Финансы Маркетинг Астрономия География Туризм Биология История Информатика Культура Математика Физика Философия Химия Банк Право Военное дело Бухгалтерия Журналистика Спорт Психология Литература Музыка Медицина |
страница 1 На правах рукописи Шульгина Тамара Михайловна Разработка и применение комплекса программ для анализа региональных изменений климата на основе данных моделирования Специальность 05.13.18 – математическое моделирование, численные методы и комплексы программ АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Томск – 2012 Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте мониторинга климатических и экологических систем Сибирского отделения РАН (ИМКЭС СО РАН)
Защита состоится 20 марта 2012 г. в 15-00 на заседании диссертационного совета Д 003.061.02 при Учреждении Российской академии наук Институте вычислительной математики и математической геофизики Сибирского отделения РАН по адресу: 630090, г. Новосибирск, пр. академика Лаврентьева, 6, Тел. (383)330-71-59. С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Института вычислительной математики и математической геофизики Сибирского отделения РАН. Автореферат разослан «___» _____________ 20__ г. Ученый секретарь диссертационного совета Д 003.061.02, д.ф.-м.н. Сорокин С.Б. ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность работы. Происходящие в последние десятилетия изменения климата [1] влекут стремительное развитие технологий их мониторинга, моделирования и прогнозирования [2]. Накапливаемые архивы метеорологических данных и разрабатываемые методики их обработки открывают перед исследователями широкие возможности для решения климатических задач различного уровня сложности. Однако, разнородность имеющейся метеорологической информации и многообразие методов ее обработки заметно усложняют процесс анализа и сопоставление получаемых результатов. Для систематизации и стандартизации процедур исследования на сегодняшний день помимо инструментариев рабочего стола (Excel, Statistica, ArcGIS, MatLab и др.) активно используются тематические веб-порталы (KNMI Climate Explorer, ECA&D EUMETNET, Giovanni GES DISC NASA, Climvis NNDC Climate Data Explorer, NCAR Command Language, ECMWF Eraclim и др.), созданные на основе современных информационно-вычислительных технологий. Однако, большинство стандартных инструментариев требуют наличия архивов данных на рабочем месте, что создает ряд трудностей их использования ввиду больших объемов архивов метеорологической информации. Архивы метеорологических данных, предлагаемые для работы специализированными веб-порталами, решают эти вопросы, но набор типовых процедур обработки этих данных дает возможность изучить лишь основные закономерности атмосферных процессов. Кроме того, задачи описания региональных климатических процессов, где необходима детализация изменений метеорологической величины в пространстве с учетом природно-климатических особенностей региона [3], порождают ряд вопросов о репрезентативности используемых метеорологических данных и накладывают требования на используемые методики анализа. Получение достоверной информации о пространственно-временном поведении метеорологических величин на территории исследуемого региона является актуальной задачей и достигается созданием и применением комплекса программ, который позволяет проводить:
Детальное описание поведения метеорологических величин выбранного региона создает надежную основу для валидации климатических моделей и разработки новых инструментов моделирования регионального климата. Кроме того, реализация такого подхода в виде комплекса вычислительных процедур, расширяет круг его применения до возможности использования при решении разного рода междисциплинарных практических задач, вызванных к жизни наблюдаемыми и прогнозируемыми изменениями климата. Рассмотренный в работе регион Сибири, в силу своих природно-климатических особенностей [4], представляет огромный интерес для проведения такого рода исследований. Таким образом, целью диссертационной работы явилась разработка комплекса программ для исследования региональных климатических изменений на основе статистического анализа метеорологических данных. Для достижения поставленной цели необходимо было решить следующие задачи:
На защиту выносится следующее:
Научная новизна, полученных автором результатов, заключается в следующем:
Практическая значимость работы определяется возможностью решения междисциплинарных задач изучения наблюдаемых и прогнозируемых изменений регионального климата. Программный комплекс может быть использован для оценки динамики метеорологических величин для любого региона. Результаты проведенного исследования дают достоверную основу для валидации существующих и создания новых региональных климатических моделей. Работа выполнялась в рамках научных проектов, поддержанных грантами РФФИ (№10-07-00547, 11-05-01190), проектом EC FP6 Enviro-RISKS (INCO-CT-2004-013427), Программой фундаментальных исследований СО РАН (проекты 4.31.1.5, 4.31.2.7) и интеграционными проектами СО РАН 4, 8, 9, 50 и 66. Результаты диссертации использованы при выполнении НИР в проектах Федеральных целевых программ (гос. контракт Минобрнауки РФ 07.514.11.4044 и гос. контракт Минобрнауки РФ 14.740.11.0409).
Объем и структура диссертационной работы. Диссертационная работа состоит из введения, трех глав, заключения и приложений. Общий объем диссертации составил 111 страниц; работа содержит 74 рисунков и 16 таблиц; список цитируемой литературы включает 94 наименования. ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность работы, сформулированы основные цели и задачи диссертационного исследования. Подчеркнута научная новизна, раскрыты теоретическая значимость и практическая ценность полученных результатов исследования. Перечислены положения, выносимые на защиту. В первой главе диссертационной работы представлено описание архивов метеорологических данных, методики выбора данных климатического моделирования, на основе их сравнения с данными инструментальных наблюдений, а также набора климатических характеристик и апробированных статистических методов для обработки и анализа этих архивов данных. Основным источником инструментальных наблюдений послужил специализированный массив данных суточной температуры воздуха и количества осадков по 223 станциям [5], расположенным на территории 35,28° – 80,60° с.ш. и 14,30° - 190,17° в.д. Максимальный временной период покрытия наблюдениями составил интервал 1.01.1881 – 31.12.2010. Также в работе был рассмотрен архив 9092с Global Synoptic Network [6], содержащий ряды почасовых измерений метеорологических величин с 2095 наземных синоптических станций на той же территории, покрывая наблюдениями период 01.01.1881 – 01.01.2000 в зависимости от срока работы станции. Данные моделирования, использованные в работе, представлены архивами полей метеорологических величин, полученными в рамках нескольких крупнейших проектов Реанализа (данные реанализа) (Таб. 1) [7]. При формировании таких наборов данных используются глобальные климатические модели, содержащие комплексные системы усвоения многолетних архивов инструментальных наблюдений. Указанные поля метеорологических величин, как правило, заданы на равномерной горизонтальной сетке и имеют глобальное покрытие по всему земному шару. Таблица 1. Характеристики данных реанализов
Кроме того, был рассмотрен набор данных APHRODITE JMA [8], содержащий поля суточных данных количества осадков за период 01.01.1979 – 31.12.2007, полученные путем интерполяции наблюдений метеорологических величин сети наземных станций на горизонтальную равномерную сетку с разрешением 0,5°×0,5°. Поскольку рассмотренные архивы данных, представляют собой результаты компьютерных расчетов, они, безусловно, имеют отличия от натурных измерений, а также различаться между собой в виду влияния на получаемые значения метеорологических переменных численного моделирования соответствующих физических процессов. Для определения полей климатического моделирования наиболее близких к натурным измерениям метеорологических величин в работе предложена методика сравнения этих полей с рядами инструментальных наблюдений. Алгоритм, реализующий эту методику, включает подготовительный этап формирования однородных рядов наблюдений, этап восстановления данных реанализа в координатах положения станций с использованием выбранного метода интерполяции и этап проверки статистических гипотез об однородности рядов интерполированных данных реанализов и инструментальных наблюдений. Ввиду того, что инструментальные наблюдения, используемые в данной работе, имеют пропуски измерений и различия в длине рядов, автором сформулирован набор условий, выполнение которых обеспечивает формирование набора рядов наблюдений, обладающих полной информацией о поведении метеорологической величины на заданном временном промежутке исследования. Согласно поставленным условиям нижняя и верхняя временные граница рядов устанавливаются в соответствии с временным интервалом наличия данных наблюдений и реанализов. Анализ информативности рядов в пределах выбранного временного диапазона включает процедуру проверки полноты измерений метеорологической величины в сутках, месяце, годе, выборке. Наблюдение метеорологической величины в пределах суток должно содержать полный набор ее измерений в соответствии со стандартами установленной системы наблюдений. В противном случае наблюдения за эти сутки исключаются из рассмотрения. В случае, если на выбранном временном отрезке произошла смена системы наблюдений, это должно быть учтено при анализе пропусков измерений метеорологической величины. В рамках месяца наблюдения могут отсутствовать не более чем за трое изолированных друг от друга суток. В противном случае наблюдения этого месяца также исключаются из рассмотрения. Каждый год выбранного временного участка должен содержать измерения метеорологической величины за все 12 месяцев. Кроме того, полученная выборка (временной ряд наблюдений станции), элементами которой являются наблюдения метеовеличины для каждого года, должна содержать не более 20% пропусков ее элементов, а приоритет отдан выборкам, имеющим не более 5% пропусков. Полученные ряды наблюдений приводятся к единому формату представления данных: единый шаг измерений по времени (суточные значения), файлы данных в формате .dbf с названием согласно идентификационному номеру станции, установленному ВМО. Этот набор условий позволяет выбрать ряды с достаточно полным объемом метеорологической информации для достоверного воспроизведения картины поведения метеорологической величины на выбранном промежутке исследования. На следующем шаге схемы выполняется интерполяция данных реанализа в координатах местоположения станций. С этой целью был произведен выбор метода интерполяции, по результатам тестирования методов билинейной интерполяции, интерполяции полиномом третьей степени, метода взвешенных обратных расстояний, модифицированного метода Шепарда и базового геостатистического крайгинг-метода на основе данных реанализов разных версий при различных размерах сетки. Наименьшие значения среднего и среднеквадратического отклонений были получены при реализации модифицированного метода Шепарда [9], основная идея которого заключается в построении сглаженной двумерной функции где точечная функция ![]() ![]() где Выбранным методом значения данных реанализа интерполированы в координатах положения станций. Полученные ряды интерполированных данных сопоставляются с данными инструментальных наблюдений в терминах различных климатических показателей, рассчитанных на основе этих рядов. Этап сравнения рядов климатических характеристик включает расчет среднего и среднеквадратического отклонений, оценку однородности рядов по критериям Хи-квадрат и Вилкоксона, а также оценку степени корреляционной связи сравниваемых рядов [10]. Полученные в результате вычислений сравнительные оценки позволяют выбрать архив данных реанализа, поля метеорологических величин которого наиболее близки к их натурным измерениям, полученным с наземных метеостанций. Для анализа пространственно-временной динамики выбранных метеорологических данных в работе сформирован и представлен перечень климатических характеристик и апробированных статистических методов, учитывающих специфику используемого климатологического материала и природно-климатические особенности региона. В рассмотрение были включены основные статистические характеристики, показатели временной структуры метеорологических рядов и показатели экстремальных явлений. Статистические характеристики представлены средним арифметическим, дисперсией, асимметрией и эксцессом. Показатели временной структуры метеорологических рядов выражены характеристиками периодических и непериодических изменений метеорологических величин во времени. К показателям периодических изменений отнесены климатические нормы (средние значения за период 1961-1990 гг.), отклонения от этих норм, суточный и годовой ход метеорологических величин. Показатели непериодических изменений представлены характеристиками непрерывной продолжительности и повторяемости атмосферных явлений со значениями метеорологических величин выше или ниже заданных пределов для разных моментов времени: Таким образом, предложенный автором подход комплексной обработки метеорологической информации обеспечивают достоверное детальное описание поведения исследуемой метеорологической величины в заданных пространственно временных масштабах. Во второй главе диссертационной работы приведено описание созданного программного комплекса, реализующего методы выбора, обработки и анализа метеорологических данных, предложенные автором и описанные в первой главе текста диссертации. Разработанный программный комплекс представляет собой модульную структуру (Рис. 1). Методика ![]() сравнения данных наблюдений и реанализа реализована в виде ряда программных модулей для формирования однородных рядов наблюдений, интерполирования значений данных реанализа в координатах отобранных метеостанций и оценки согласованности этих значений с величинами натурных измерений. Программные модули формирования рядов наблюдений включают процедуру компоновки рядов в результате анализа пропусков измерений метеовеличин на различных временных сегментах (сутки, месяц, год, период исследования) по сформулированным условиям Рис.1. Блок-схема комплекса программ. отбора, процедуру приведения полученных рядов измерений к единому формату представления данных (структура и формат файлов данных). Программные модули разработаны с использованием высокоуровневого языка программирования Python 2.6, выбор которого обусловлен простотой написания кода и широкими возможностями стандартной библиотеки функций. Программы, реализующие процедуры выбора метода интерполяции, с последующим его применением для восстановления данных реанализа в координатах метеостанций, оценкой согласованности интерполированных и инструментально измеренных метеовеличин, разработаны с использованием языка программирования высокого уровня ITTVIS Interactive Data Language 7.0 (IDL), ввиду наличия в нем обширной библиотеки процедур для математической и графической обработки данных. Результаты сравнения данных представлены числовыми значениями отклонений. Программные модули обработки и анализа выбранных метеорологических данных обеспечивают расчет климатических характеристик и применение выбранных статистических методов, сформированных автором и представленных в первой главе текста диссертации. Визуализация результатов вычислений производится на каждом этапе последовательности математических и/или статистических процедур обработки. Для каждого программного модуля в работе приводятся блок-схемы, демонстрирующие принципы их работы, а также подробное описание входных и выходных данных. Разработанные модули являются независимыми относительно обрабатываемой метеорологической информации и пространственно-временных границ. Гибкая структура полученного комплекса программ позволяет добавлять в него новые модули. В заключительной части главы представлено описание роли разработанного комплекса программ в информационно-вычислительной веб-ГИС системе, создаваемой в ИМКЭС СО РАН для проведения научных исследований, связанных с обработкой и анализом архивов пространственно-привязанных геофизических данных, полученных в результате наблюдений и моделирования. В третьей главе диссертационной работы представлены результаты исследования динамики приземной температуры воздуха и количества осадков на территории Сибири в период с 1979 по 2007 годы, полученные с использованием созданного комплекса программ. Анализ состояния данных метеорологических величин (наблюдения и реанализы) для территории Сибири дал следующие результаты. В результате применения процедуры выбора однородных рядов наблюдений, для описания приземной температуры воздуха и количества осадков на территории Сибири в период 1961-2010 гг. только 59 временных рядов предоставляют необходимый для проведения исследований объем метеорологической информации (не более 5% пропусков измерений во временном ряде). Различие данных приземной температуры из различных версий реанализов составило в среднем 0,3 °С между данными версий ECMWF (ERA-40, ERA INTERIM) и NCEP (NCEP/NCAR, NCEP/DOE), с наименьшей величиной различий в теплый сезон года (0,06 °С), и наибольшей – в зимний сезон (1,5 – 2 °С). Анализ этих различий в широтном диапазоне показал заметное их увеличение с ростом широты в зимний период года (5,4 °С в Западной Сибири, 7,7 °С – в Восточной Сибири). В летний сезон эта разница существенно меньше. Анализ данных количества осадков из архивов ERA INTERIM и APHRODITE JMA выявил отличия данных друг от друга в среднем в два раза. Изучение пространственного распределения значений количества осадков не выявило широтной зависимости различий. В результате сравнения данных реанализов с выбранными 59 рядами инструментальных наблюдений по предложенной схеме сравнения, для описания поведения температуры воздуха на территории Сибири был выбран архив реанализа ECMWF ERA INTERIM (Рис.2а). Полученные величины различий в данных (станции с номерами 44 - 47 и 55) связаны с орографией положения станций (наибольшие различия – в горных районах Восточной Сибири, наименьшие – на Западно-Сибирской равнине). Для описания поведения количества осадков автором выбран архив данных APHRODITE JMA (Рис.2.б). Однако, результаты исследования осадков, с использованием этого набора данных следует принимать с осторожностью, поскольку при сравнении их с рядами натурных наблюдений однородность на половине станций не была достигнута.
широтах Восточной Сибири (0,6 – 1,0 °С/10 лет) и слабой тенденцией к похолоданию в ее центральной части – до -0,2 °С/10 лет. Динамика температур весеннего и осеннего сезонов, сопровождающаяся ростом температур весной до 1,2 °С/10 лет (в апреле) на востоке, до 1,1 °С/10 лет (в мае) на западе Сибири и осенью – до 1,4 °С/10 лет на востоке (в ноябре) и 1,1 °С/10 лет (в октябре) на западе Сибири, нашла отклик в колебаниях продолжительности вегетационного периода года (Рис. 3.а). На севере Сибири и к югу Западной Сибири наблюдается увеличения числа дней со среднесуточной температурой выше 5 °С до 15 дней/10 лет. Сумма эффективных температур ( Динамика количества осадков холодного (ноябрь – апрель) и теплого (май – октябрь) сезонов на территории Сибири не выражена. Количество осадков холодного периода года характеризуется их увеличением на большей части региона Сибири в среднем на 2 мм/10 лет. Отрицательные величины тренда наблюдаются лишь в верхних широтах Западный Сибири, достигая -2,5 мм/10 лет в марте и -7 мм/10 лет в апреле. Для теплого периода года отмечен рост количества осадков на территории Средней и Восточной Сибири, достигая наибольших значений, 18 мм/10 лет, в июле и августе. Однако территории средних широт Западной Сибири характеризуется уменьшением количества осадков до -6 мм/10 лет (в июне). Анализ интенсивности осадков (Рис.3.б) не выявил значимых изменений числа дней
Рис.3. Тренд продолжительности вегетационного периода года (а), тренд интенсивности количества осадков (б). с осадками выше 1 мм, за исключением северного региона, где наблюдается их рост до 5 дней/10 лет. Сопоставление наблюдающихся изменений температуры воздуха и количества осадков позволило сделать предварительные выводы о том, что на юге территории Сибири формируются благоприятные условия для увеличения продуктивности растительности. В Заключении диссертационной работы сформулированы основные результаты работы, которые заключаются в следующем:
ПУБЛИКАЦИИ АВТОРА ПО ТЕМЕ ДИССЕРТАЦИИ
1 Анализ изменений климата и их последствий / Труды ГУ «ВНИИГМИ-МЦД», под ред. Разуваева В.Н., Шерстюкова Б.Г. – Обнинск, 2010. – 278 с. 2 Дымников В.П., Лыкосов В.Н.и др. Моделирование климата и его изменений // «Современные проблемы вычислительной математики и математического моделирования». – М.: Наука, 2005. – Т. 2. – С. 36-175. 3 IPCC: 2007. Climate change. The physical science basis. Contribution of Working Group I to the Forth Assessment Report of the Intergovernmental Panel on Climate Change. 4 Кабанов М.В., Лыкосов В.Н. Мониторинг и моделирование природно-климатических изменений в Сибири // Оптика атмосферы и океана, 2006. – Т. 19, № 6. – С. 753-764. 5 Булыгина О.Н., Разуваев В.Н. Описание массива данных суточной температуры воздуха и количества осадков на 223 метеорологических станциях, http://www.meteo.ru/climate/descrip1.htm 6 Data documentation for dataset 9290c. Global Synoptic Climatology Network. The firmer USSR. National Climatic Data Center, 2005. 7 Jenne R. Reanalysis; The observations and analysis // Large scale observations: a SEARCH workshop. – Seattle, WA, 2001. – 7 p. 8 APHRODITE JMA, http://www.chikyu.ac.jp/precip/data/APHRO_V1003R1_readme.txt 9 Renka, R. J. Algorithm 790 - CSHEP2D: Cubic Shepard Method for Bivariate Interpolation of Scattered Data // ACM Trans. Math Softw. – V. 25. – P. 70-73. 10 Storch H. von, Zwiers F.W. Statistical analysis in climate research / Cambridge University Press, Cambridge, 1999. – 484 p. 11 Дюкарев Е.А. Изменчивость климата Западной Сибири во второй половине ХХ века // География и природные ресурсы, 2005. – № 4. – С. 70 - 75. Смотрите также: Разработка и применение комплекса программ для анализа региональных изменений климата на основе данных моделирования
226kb.
1 стр.
Разработка методик защиты программ от анализа и модификации на основе запутывания кода и данных
262.98kb.
1 стр.
Разработка программных средств для интерактивного анализа публикаций на основе olap-технологии
27.73kb.
1 стр.
Применение графов для анализа сложных систем на основе исходного кода программ
96.04kb.
1 стр.
Отчет о лаботарорной работе методы и средства анализа данных по теме: «Система анализа данных weka»
383.87kb.
2 стр.
Вычислительная среда r как платформа для моделирования и обучения школьников и студентов
51.54kb.
1 стр.
Природа квазирегулярности катастрофических изменений климата и их влияние на биоразнообразие земли
292.52kb.
1 стр.
Разработка и применение лекарственных средств на основе ароксиалканкарбоновой кислоты при болезнях животных, вызываемых условно патогенной микрофлорой 16. 00. 03 ветеринарная микробиология, вирусология, эпизоотология
729.2kb.
6 стр.
Лабораторные работы по дисциплине "Теория экономических информационных систем"
95.98kb.
1 стр.
Программный комплекс для математического моделирования ядерно-физических приборов и устройств
19.01kb.
1 стр.
Отчет о лаботарорной работе по дисциплине Методы и средства анализа данных по теме: «Система анализа данных weka»
229.16kb.
1 стр.
Численное моделирование и разработка комплекса программ исследования теплообмена и ламинарного течения в регулярных продольнооребренных коридорных структурах 465.83kb.
2 стр.
|