Министерство образования и наук
и Российской Федерации
Российский Государственный Гидрометеорологический Университет
На правах рукописи
УДК 551.509.314(215 – 17)
Борисова Алла Семеновна
СТАТИСТИЧЕСКИЙ АНАЛИЗ И ПРОГНОЗ
ЕСТЕСТВЕННЫХ КОЛЕБАНИЙ ПОЛЕЙ ГЕОПОТЕНЦИАЛА
ПОВЕРХНОСТИ 500 ГПА В СЕВЕРНОМ ПОЛУШАРИИ
Специальность 25.00.30 – метеорология, климатология, агрометеорологияАвтореферат диссертации
на соискание ученой степени кандидата географических наук Санкт – Петербург 2007 2 Диссертация выполнена на кафедре метеорологических прогнозов Российского Государственного Гидрометеорологического УниверситетаНаучный руководитель: Доктор физико – математических наук, профессор Репинская Раиса Петровна
Официальные оппоненты: Доктор географических наук, профессор, Хандожко Леонид Андреевич Кандидат физ.-мат. наук, старший научный сотрудник Научного Фонда “Нансен-центр” Кузьмина Светлана Игоревна
Ведущая организация: Санкт – Петербургский Государственный Университет, факультет географии и геоэкологии
Защита состоится “25” октября 2007 г. в часов минут на заседании диссертационного совета Д 212.197.01 Российского Государственного Гидрометеорологического Университета Адрес: 195196, Санкт – Петербург, Малоохтинский пр., 98.
С диссертацией можно ознакомиться в библиотеке Российского Государственного Гидрометеорологического Университета.
Автореферат разослан “ ” сентября 2007 г.
Ученый секретарь диссертационного совета Доктор физико – математических наук, профессор А.Д. Кузнецов
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы Изменение атмосферных процессов, и, следовательно, погоды и климата, определяют динамические и физические процессы, протекающие в атмосфере, а также взаимодействия с окружающей средой: поверхности суши, океана и льда.
С практической точки зрения особый интерес представляют волны синоптического масштаба и связанные с ними синоптические вихри, играющие важную роль в формировании погодных условий на больших территориях. В системе общей циркуляции атмосферы синоптические вихри выполняют функцию основного механизма междуширотного воздухообмена.
Происходит меридиональный перенос тепла, водяного пара и импульса в планетарном масштабе, что приводит к сглаживанию температурных контрастов между высокими и низкими широтами и вызывает значительные изменения погодных условий. Анализ естественных колебаний полей геопотенциала и его прогноз играет важную роль в исследовании циркуляционных процессов в Северном полушарии. Для разработки как гидродинамических, так и физико – статистических моделей прогноза, а также для статистической параметризации физических процессов необходима детальная информация о макро – и мезоструктуре полей различных метеорологических величин. В значительной мере это “белое пятно” в исследовании циркуляционного режима атмосферы.
Цель и задачи диссертационной работы Цель диссертационного исследования: дать всестороннюю характеристику полей геопотенциала поверхности 500 гПа и разработать статистическую модель для прогноза естественных колебаний полей геопотенциала в Северном полушарии.
Для достижения этой цели необходимо было решить следующие задачи:
- Описать пространственно – временную структуру полей геопотенциала при помощи метода эмпирических ортогональных составляющих (ЭОС) и выявить климатические особенности исследуемых полей.
- Исследовать изменчивость рассматриваемого временного ряда при помощи вейвлет анализа;
- Исследовать наличие трендов в рядах полей геопотенциала поверхности гПа;
- Разработать схему прогноза временных рядов синоптически важных ЕК;
- Исследовать статистическую предсказуемость коэффициентов разложения полей АТ – 500 гПа по ЭОС, описывающих эволюцию наиболее мощных естественных колебаний (ЕК) при переходе от месяца к месяцу и от года к году.
Научная новизна работы - Использован метод ЭОС для сжатия исходной информации с целью компактного описания полей геопотенциала и фильтрации метеорологических шумов для Северного полушария;
Оценена роль процессов глобального, синоптического и мезометеорологического масштабов в формировании полей геопотенциала поверхности 500 гПа;
- Выполнен анализ статистически значимых линейных трендов в ходе суточных, месячных и годовых полей геопотенциала;
- Разработаны авторегрессионные модели различных порядков для разных широтных зон Северного полушария на основе нелинейного подхода по Боксу и Дженкинсу;
- Исследована статистическая предсказуемость рядов ЕК полей геопотенциала поверхности 500 гПа.
Практическая значимость работы Практическая значимость состоит в том, что 1) на основе работы можно создать методику и провести расчеты изменений циркуляционных режимов Северного полушария 2) разработанная модель может быть внедрена в практику оперативной службы долгосрочных прогнозов Госкомгидромета.
3) полученные результаты могут быть использованы в учебном процессе при преподавании таких дисциплин, как долгосрочные методы прогноза погоды, методы статистической обработки метеорологической информации и принятия статистических решений и климатологии.
Апробация результатов Наиболее значимые результаты докладывались на II-ой Всероссийской молодежной научной конференции “История науки и техники” (2002г.), межкафедральных научных семинарах (март 2006 г., октябрь 2006 г.).
По теме диссертации опубликованы 4 работы, 1 принята к публикации.
На защиту выносятся - Целесообразность использования недетерминированных базисных функций для аналитического представления совокупности ежедневных и среднемесячных полей геопотенциала;
- Оценка сходимости разложения полей геопотенциала при переходе от месяца к месяцу;
- Результаты разложения полей геопотенциала и анализа их сходимости;
- Доминирование крупномасштабных ЕК геопотенциала, позволяющее выделить процессы различных масштабов и оценить дисперсионную нагрузку «климатического шума»;
- Результаты разработки авторегрессионных моделей наиболее мощных ЕК геопотенциала нелинейным методом Бокса и Дженкинса;
- Предсказуемость на основе «модель плюс шум» ЕК; в средних и высоких широтах ЕК прогнозируются лучше, чем это дает климатический прогноз, в отличие от тропической зоны, где на ЕК воздействуют более мелкомасштабные процессы;
- Оценка линейного тренда в ходе наиболее крупных ЕК свидетельствует об отсутствии четкой тенденции во временной эволюции поля геопотенциала;
- Оценка прогноза по АР - моделям ЕК полей геопотенциала поверхности 500 гПа.
Структура и объем диссертации Диссертация состоит из введения, шести глав, заключения, приложений и списка использованных источников, включающего 123 наименования. Объем диссертации составляет страниц.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертации, сформулированы цель и задачи исследования, приведены основные положения и результаты, выносимые на защиту.
В первой главе рассматриваются особенности волновых движений в различных широтных зонах Северного полушария.
Во второй главе дано описание пространственно – временной структуры полей геопотенциала поверхности 500 гПа при помощи метода эмпирических ортогональных составляющих (ЭОС) и выявлены климатические особенности исследуемых полей.
Представление полей геопотенциала при помощи ЭОС было осуществлено с целью получения более сглаженных и, следовательно, более предсказуемых характеристик временных рядов. Суть метода эмпирических составляющих (функций) заключается в том, что рассчитываются временные ковариационные матрицы с различным временным сдвигом и затем находятся собственные числа и собственные векторы этих матриц. Кроме того, важно подчеркнуть, что в отличие от классических ортогональных многочленов в этом случае требование минимума ошибки аппроксимации налагается не на коэффициенты разложения, а на собственные векторы матрицы ковариации. Иначе говоря, разложение временного ряда будет производиться по собственному базису матрицы. Таким образом, базис можно рассматривать как систему детерминированных ортогональных функций, порожденных климатической системой. Основная особенность метода ЭОС заключается в экономном представлении наиболее существенной информации малым числом членов ряда, для чего вводится система ортонормированных функций, значения которых зависят от статистической структуры исследуемых полей.
Метеорологическое поле задается значениями своего элемента в ряде произвольно выбранных дискретных точках х=1, 2,…, п, т. е. совокупность метеорологических полей можно описать некоторой функцией F(t,x);
t=1,2,…,m- фиксированное для каждого индивидуального поля время, x=1, 2,…, n - параметр - координата.
Функция F(t,x) представлена в виде:
где каждая строка описывает в п заданных точках состояние исследуемого поля в некоторый момент времени.
Необходимо найти разложение любого поля из совокупности по функциям Хh(x) (h=1,2,…) с коэффициентами Т(t), изменяющимися от одного поля к другому, то есть должно иметь место разложение при условии, что сумма квадратов ошибок разложения по всем точкам данной совокупности полей достигает минимума при данном фиксированном h.
Функции Хh(x) являются собственными функциями данной совокупности полей.
Исследуемые поля геопотенциала поверхности 500 гПа были аппроксимированы однократными рядами по ЭОС и был проведен анализ скорости сходимости рядов разложений в различные месяцы года. При этом определялось число спектральных компонентов, отображающих основную часть дисперсии исходных временных последовательностей.
В процессе анализа результатов разложений рассматриваются два параметра: часть выборочной дисперсии j%, приходящуюся на каждую волновую моду (ЕК), и кумулятивная дисперсия j%.
Формулы для вычисления параметров j и j имеют вид:
где j – соответствующее собственное число корреляционной матрицы, а j = i=1,S i - след корреляционной матрицы, который при j=S интерпретируется как 100% - я выборочная дисперсия.
Детальному анализу подвергались колебания, несущие 5% выборочной дисперсии, так как именно наиболее крупные из них описывают влияние зонального и меридионального переносов, а также наиболее существенных региональных циркуляционных механизмов на различные процессы, протекающие в атмосфере и на подстилающей поверхности.
Рассматривался годовой ход ЕК, несущих более 5% выборочной дисперсии, где четко видна доминирующая роль первого ЕК, на долю которого приходится наибольшая часть дисперсии (таблица 1 – приведена одна из таблиц).
Основное влияние на формирование и эволюцию первого энергетически наиболее мощного ЕК оказывают: наличие зонального переноса; ВЗК и примыкающие к ней зоны пассатов; климатические и сезонные центры действия атмосферы; муссонная циркуляция и меридиональный перенос, обусловленный подвижными циклонами и региональными тропическими циркуляциями; распределение источников тепла над океанами и континентами; конфигурация континентов и рельеф поверхности.
Таблица 1 - Оценка сходимости рядов разложений естественных колебаний полей геопотенциала поверхности 500 гПа по ЭОС на 5 с.ш. Северного полушария 15.22442.29 42.29 15.427 42.85 42.85 12.34134.28 34.28 10.62 29.50 29. 4.702 13.06 55.35 4.544 12.62 55.48 5.579 15.50 49.78 5.103 14.18 43. 3.336 9.27 64.62 3.068 8.52 64.00 3.525 9.79 59.57 3.699 10.28 53. 10.58429.40 29.40 12.28 34.11 34.11 10.98530.51 30.51 10.713 29.76 29. 4.762 13.23 42.63 4.758 13.22 47.33 4.79 13.31 43.82 4.807 13.35 43. 3.719 10.33 52.96 3.952 10.98 58.31 4.071 11.31 55.13 4.028 11.19 54. 11.20331.12 31.12 10.886 30.24 30.24 9.789 27.19 27.19 13.723 38.12 38. 4.71 13.08 44.20 4.83 13.42 43.66 5.47 15.19 42.39 5.133 14.26 52. 3.933 10.93 55.13 3.993 11.09 54.75 3.75 10.42 52.80 3.261 9.06 61. Анализ полей геопотенциала поверхности 500 гПа, описываемых естественными колебаниями, показал, что существует четко выраженная граница между процессами различных масштабов; в средних и тропических широтах и непосредственно на экваторе указанная граница выражена более резко, чем в высоких широтах. На мезомасштабном участке спектра атмосферных волновых мод имеет место монотонная сходимость аппроксимирующих рядов, что свидетельствует о бесконечном многообразии флуктуаций мод подсиноптического и мезометеорологического масштабов.
Исключением является экваториальный пояс, где наблюдается четкая граница между процессами синоптического и мезомететорологического масштаба. Почти на всех широтах в переходные периоды от холодного полугодия к теплому (апрель – май) и от теплого к холодному (октябрь ноябрь) наблюдается наименее четкая граница между процессами различных масштабов.
Показано, что в экваториальной зоне, тропиках и полярных широтах (вплоть до =80 с.ш.) достаточно первых пяти ЕК в течение всего года, чтобы описать доминирующую часть дисперсии. В средних широтах необходимое количество ЕК увеличивается до восьми в зимнее полугодие и девяти в летнее. Указанная разница свидетельствует о значительной флуктуации ЕК разных масштабов на сезонном и годовом интервалах времени.
В таблице 2 представлены значения выборочной дисперсии первых двух ЕК в различные сезоны. Следует отметить, что в умеренных широтах в теплое полугодие доля дисперсии, приходящаяся на первое ЕК, уменьшается.
Таблица 2 - Значения выборочной дисперсии первых двух ЕК в различные сезоны При рассмотрении выборочной доли дисперсии d1, описываемой первым ЕК, обнаружено, что максимальная часть дисперсии, описываемая первым ЕК, в течение года наблюдается в июне на 20-45 с.ш., зимой максимум дисперсии d1 наблюдается на экваторе и в тропиках.
В большинстве случаев начиная с третьего – четвертого членов имеет место монотонная сходимость аппроксимирующих рядов, что свидетельствует о бесконечном многообразии флуктуаций атмосферных мод более мелкого масштаба, что затрудняет усечение спектра зональных естественных колебаний с целью фильтрации метеорологического шума.
Исключением является экваториальный пояс, где наблюдается четкая граница между процессами различных масштабов. Здесь, как известно, существенное влияние на атмосферные процессы различных масштабов оказывает положение ВЗК и ее сильная сезонная миграция.
В третьей главе рассмотрены аналогии и различия вейвлет – анализа и анализа Фурье, приведены различные типы вейвлет – преобразования (ВП). С целью выделения интервала периодов сильных колебаний во временных рядах на различных широтах был проведен вейвлет – анализ коэффициентов разложения естественных колебаний полей геопотенциала на различных широтных кругах Северного полушария при использовании базисного вейвлета Morlet.
При проведении вейвлет – анализа были рассмотрены пять первых коэффициентов разложения по ЭОС полей геопотенциала за период с 1964 по 1986 годы на следующих широтах Северного полушария: 10с.ш., 35с.ш., 55с.ш., 75с.ш.
Максимально возможный масштаб по времени вычисляется по формуле где п – объем выборки, dt - шаг по времени, s0 – минимальный масштаб по времени. Для вейвлета Morlet s0 = dt, dj – шаг по времени (обычно dj = 0.25, но для повышения точности можно выбрать dj = 0.125). В данной работе максимально возможный масштаб равен 66 месяцам. Рассматривались следующие масштабы: jtot =47; для выявления высокочастотных колебаний масштаб уменьшали до 37; для выявления низкочастотных колебаний масштаб увеличивали до 57. На рисунке 1 представлены амплитуды ВП первого коэффициента разложения (Т1) по ЭОС полей геопотенциала на 10с.ш. при использовании вейвлета Morlet.
Рисунок 1 – Амплитуды вейвлет – преобразования первого коэффициента разложения среднемесячных полей геопотенциала поверхности 500 гПа на 10 с.ш. на масштабах по времени а) jtot = 47, б) jtot = 57, в) jtot = Анализ портретов вейвлет - преобразования временных рядов коэффициентов разложения полей геопотенциала поверхности 500 гПа на различных широтах Северного полушария показал, что основные волновые моды, описывающие наибольшую дисперсию, испытывают колебания: в тропиках - с периодом от 25 до 50 месяцев (рисунок 1), обусловленные явлением Эль – Ниньо и квазидвухлетней цикличностью. На широте можно выделить колебания с периодом до 8 месяцев в течение всего интервала лет, возможно, эти колебания связаны с тропическими циркуляциями синоптического масштаба (ВЗК, муссоны). В средних и полярных широтах Северного полушария в течение всего интервала наблюдаются колебания с периодом до одного года.
Показано, что в тропиках до 1975г. существуют колебания с периодом от 25 до 130 месяцев, связанные с двухлетним и 11 – летним солнечным циклами. В высоких широтах с 1968 года ярко выражены колебания с периодом 25 месяцев, что также связано с двухлетним солнечным циклом.
В четвертой главе описано построение авторегрессионной модели временного ряда и дана оценка параметров модели. Более предпочтительным для климатических временных рядов является использование моделей авторегрессии, параметры которых оцениваются методом Берга-Левинсона, а оптимальный порядок модели выбирается по критерию Парзена. Это метод идентификации и оценивания реализован для вычисления параметров моделей коэффициентов разложения зональных полей геопотенциала поверхности 500 гПа в Северном полушарии (широтных кругов 10 0, 350, 550, 750) за период 1964-1986 гг. Фрагмент результатов идентификации представлен в таблице 2.
Таблица 2 - Порядки АР – моделей временных рядов коэффициентов разложения по ЭОС (Т1 – Т5) для различных широт Северного При анализе таблицы видно, что временные ряды коэффициентов разложения описываются моделями авторегрессии второго и более высокого порядков, авторегрессионные модели для младших компонентов разложения имеют такой же высокий порядок, как и первые компоненты, таким образом, все рассмотренные временные ряды допускают прогнозирование.
При решении задачи идентификации моделей климатических временных рядов использовался метод максимального правдоподобия, разработанный Боксом и Дженкинсом. При подгонке параметрических моделей к реальным климатическим временным рядам суммарный порядок оптимальной смешанной модели оказывается невысоким, что позволяет выбирать алгоритм Бокса-Дженкинса в качестве основного средства для оценивания параметров смешанных моделей. Однако надежность оценок снижается с ростом порядка модели. Высокий порядок построенной модели для тропических широт свидетельствует о том, что надежность оценок в этих широтах ниже, чем в средних и высоких широтах Северного полушария. Для средних и высоких широт порядок модели варьируется от 2 до 7.
В пятой главе решался вопрос о наличии трендов в рядах полей геопотенциала поверхности 500 гПа с целью выявления тенденций в ходе различных ЕК.
Анализ и статистическое описание динамики какого – либо физического процесса начинают с выявления тренда и оценки его параметров. Форма тренда отражает закономерности эволюции изучаемого процесса. При анализе тренда, в первую очередь, выявляют форму тренда, затем выбирают то уравнение, которое наилучшим образом аппроксимирует тренд.
Результаты анализа трендов циркуляционных волновых мод в различных климатических зонах Северного полушария представлены в таблице 3.
Наблюдаются различные по знаку и по величине углы наклона линейного тренда, причем наибольший по величине статистически значимый положительный тренд для К1 имеет место в экваториальных и полярных широтах (=0.0009 и 0.001 соответственно), что говорит о значительных изменениях в поле геопотенциала (с течением времени значения временного ряда имеют тенденцию возрастать). В субтропиках (35 с.ш.) тренд отрицателен для всех К, что означает понижение давления и, скорее всего, температуры в этом широтном поясе Северного полушария. Для К3 в субтропиках статистически значимого тренда не наблюдается, т. е. тренд не является выраженным, и использование роста или убывания этого тренда при прогнозировании динамки временного ряда на будущее невозможно.
Таблица 3 – Оценка тренда в ходе сигналонесущих ЕК В шестой главе дано описание методики определения статистической предсказуемости временных метеорологических последовательностей.
Представлены результаты анализа статистической предсказуемости сигналонесущих естественных колебаний полей геопотенциала поверхности 500 гПа в различных широтных зонах Северного полушария. Дана оценка прогноза естественных колебаний полей геопотенциала поверхности 500 гПа.
Статистическую предсказуемость процесса характеризуют ошибкой его оптимального линейного прогноза. Для ее количественного описания используют два критерия: относительную погрешность прогноза с одношаговым упреждением d p,q(1) и предел статистической предсказуемости p,q. Индексы p и q указывают, что эти критерии отвечают параметрическим моделям авторегрессии – скользящего среднего АРСС(p, q) порядка (p, q). В соответствии с дисперсией ошибки D ( ) = a2 j2 относительная ошибка прогноза при = и удовлетворяет неравенству: 0 dp,q(1) 1. Случаи dp,q(1) = 0 и dp,q(1) = соответствуют детерминированному процессу и белому шуму.
Относительная погрешность прогноза dp,q() и коэффициент корреляции r ( ) между истинным значением процесса и его прогнозом с упреждением связаны соотношением Это соотношение использовалось для оценки предсказуемости процесса.
Предел статистической предсказуемости p,q определяется как упреждение прогноза, при котором его средняя квадратическая ошибка меньше дисперсии прогнозируемого процесса. Значит, предел предсказуемости достигается, если относительная ошибка прогноза не меньше принятого уровня : dp,q() при p,q. Если ошибка прогноза распределена нормально, то доверительный интервал для экстраполяционной функции определяется в виде где t+;, t; – верхняя и нижняя граница доверительного интервала соответственно; µ1 – квантиль нормального распределения при уровне значимости.
Результаты оценки статистической предсказуемости первых пяти естественных колебаний приведены в таблице 4 со следующими обозначениями:
Таблица 4 – Оценка статистической предсказуемости ЕК полей геопотенциала на различных широтах Северного полушария Е2 – нормированная дисперсия ошибки прогноза, Е – нормированный стандарт ошибки прогноза, L0,95, L0,90, L0,85 – предел предсказуемости исследуемой последовательности при ошибке прогноза Е=0.95, 0.90 и 0.85.
Анализ таблицы показывает, что предсказуемость в различных широтных зонах не превышает 6 – 8 суток.
На основе модели «сигнал плюс шум» самая слабая предсказуемость сигнала (3 суток при Е=0.95, при уменьшении ошибки - 2 и 1.5 суток) наблюдается в экваториальных широтах.
В умеренных и полярных широтах предсказуемость может достигать 8 суток при Е=0.95, (6 суток при Е=0.90 и 4.5 суток при Е=0.85) и является максимальной предсказуемостью естественных колебаний из всех рассматриваемых, что лучше климатического прогноза L, который в этих широтах составляет 6 суток.
Такой же вывод можно сделать при анализе значений ошибки прогноза на один шаг вперед рядов коэффициентов по АР – модели. Лучшее качество прогнозов достигается на 55 и 75 с.ш., где = 0.49. Предел предсказуемости сигнала при = 0.85 в этих широтах оказывается максимальным и превышает 4 суток.
При анализе прогноза коэффициентов разложения полей геопотенциала поверхности 500 гПа по разработанным авторегрессионным моделям можно отметить следующее: в средних и высоких широтах наблюдается хорошее соответствие фактических и вычисленных по моделям значений коэффициентов разложений полей геопотенциала поверхности 500 гПа. В низких широтах ход вычисленных по моделям значений сильно сглажен, и, следовательно, использование такого прогноза не целесообразно.
Рисунок 2 –Вычисленные по моделям и фактические значения первого коэффициента разложения Т1 полей геопотенциала поверхности 500 гПа Январь 1986 г. ( временной ход Т1, ---- прогноз Т1) а) 10 с.ш. b) 35 с.ш. c) 55 с.ш. d) 75 с.ш.
В заключении сформулированы основные выводы диссертации.
• Для аппроксимации полей геопотенциала в различные месяцы года целесообразно применять метод ЭОС, позволяющих избавиться от шума;
• Для описания погодообразующих естественных колебаний полей геопотенциала в различных климатических зонах Северного полушария необходимо учитывать первые пять ЕК в экваториальной зоне, тропиках и полярных широтах (вплоть до =80 с.ш.) в течение всего года. В средних широтах необходимое количество ЕК увеличивается до восьми в зимнее полугодие и девяти в летнее. Указанная разница свидетельствует о значительной флуктуации ЕК разных масштабов на сезонном и годовом интервалах времени;
• Вследствие достаточно резкого скачка в значениях части суммарной дисперсии, которую несут различные ЕК во все сезоны и на всех параллелях полушария, можно достаточно уверенно выделить процессы различных масштабов. Летом (июнь) максимум дисперсии, описываемый первым ЕК, наблюдается на 20-45 с.ш. В этот сезон происходит усиление азорского и гононульского максимумов и их миграция к северу, а также их приближение к материкам.
Зимой максимум дисперсии d1 наблюдается на экваторе и в тропиках.
Граница между естественными колебаниями синоптического и мезометеорологического масштабов прослеживается, как правило, менее четко, вследствие монотонной сходимости рядов разложения во всех месяцах года и на всех широтах, что усложняет усечение спектра естественных колебаний с целью фильтрации "метеорологического шума". Исключением является экваториальный пояс;
• Анализ портретов ВП временных рядов коэффициентов разложения полей АТ – 500 гПа на различных широтах Северного полушария показал, что основные волновые моды, описывающие наибольшую дисперсию, испытывают колебания: в тропиках - с периодом от 25 до 50 месяцев, обусловленные явлением Эль – Ниньо и квазидвухлетней цикличностью. На широте 35 можно выделить колебания с периодом до 8 месяцев в течение всего интервала лет. В средних и полярных широтах Северного полушария в течение всего интервала наблюдаются колебания с периодом до одного года;
• Показано, что в тропиках до 1975г. существуют колебания с периодом от 25 до 130 месяцев, связанные с двухлетним и 11 – летним солнечным циклами. В высоких широтах с 1968 года ярко выражены колебания с периодом 25 месяцев, что также связано с двухлетним солнечным циклом;
• Для первых пяти естественных колебаний полей АТ – 500 гПа на различных широтах Северного полушария были синтезированы авторегрессионные модели различных порядков. Все исследуемые волновые моды описаны моделями второго и более высоких порядков, что свидетельствует об отсутствии белошумовых и марковских процессов первого порядка. В тропических широтах порядок построенной модели выше, чем в средних и высоких широтах, что говорит о ненадежности оценок и сильной зашумленности процессов в этом широтном поясе;
• Анализ трендов циркуляционных волновых мод в различных климатических зонах Северного полушария показал, что в субтропиках тренд отрицателен для всех рассматриваемых коэффициентов разложения, что означает понижение оценок геопотенциала в этом широтном поясе Северного полушария;
• Установлено, что в средних и высоких широтах ЕК прогнозируются лучше, чем это дает климатический прогноз, в отличие от тропической зоны, где на ЕК воздействуют более мелкомасштабные процессы;
• При оценке прогноза по разработанным авторегрессионным моделям видно, что на широтах 55 и 75 с.ш. наблюдается хорошее соответствие фактических и вычисленных по моделям значений коэффициентов разложений полей геопотенциала поверхности 500 гПа.
По теме диссертации опубликованы следующие работы:
1. Третьякова А.С., Кузнецова Е.С. Использование вейвлет – анализа для использования временных метеорологических рядов // Сборник трудов II – ой Всероссийской молодежной научной конференции “История науки и техники” 4-5 декабря 2002 г. – СПб.: изд. СПбГУ – 2003. – 225 с.
2. Репинская Р.П., Третьякова А.С. Вейвлет - анализ естественных колебаний среднемесячной температуры воздуха на различных параллелях Северного полушария // Материалы итоговой сессии Ученого совета 27 – января 2003 г. Информационные материалы. Часть I. Секции метеорологии и гидрологии. – СПб.: изд. РГГМУ, 2003. – С.58 – 59.
3. Repinskaya R.P., Tretyakova A.S. Spectrum analysis of natural fluctuations of the monthly-mean air temperature over the northern hemisphere //Mon.
Geograph. J.: Dijon.: Universit de Bourgogne (принята к опубликованию).
4. Борисова А.С. Некоторые методы построения и идентификации математических моделей временных климатических рядов // Научнотехнические ведомости СПБГПУ. - №4.- 2007. – С. 138-141.
5. Борисова А.С. Методы анализа и прогноза климатических временных рядов // Известия вузов. Северо – Кавказский регион. Естественные науки. - №6. – 2007.