ОЦЕНКА ДОЛИ И ГОДОВОЙ АМПЛИТУДЫ АДДИТИВНЫХ СЕЗОННЫХ КОМПОНЕНТ ВРЕМЕННЫХ РЯДОВ ПРИЗЕМНОЙ ТЕМПЕРАТУРЫ ВОЗДУХА 30-ТИ ПУНКТОВ НАБЛЮДЕНИЙ РФ ЗА 120-ЛЕТНИЙ ПЕРИОД

Научная статья
Выпуск: № 7 (38), 2015
Опубликована:
2015/08/15
PDF

Седов В.Е.

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

ОЦЕНКА ДОЛИ И ГОДОВОЙ АМПЛИТУДЫ АДДИТИВНЫХ СЕЗОННЫХ КОМПОНЕНТ ВРЕМЕННЫХ РЯДОВ ПРИЗЕМНОЙ ТЕМПЕРАТУРЫ ВОЗДУХА 30-ТИ ПУНКТОВ НАБЛЮДЕНИЙ РФ ЗА 120-ЛЕТНИЙ ПЕРИОД

Аннотация

Анализируется годовой цикл приземной температуры воздуха на основе аддитивных  сезонных компонент, восстановленных из  временных рядов 30-ти пунктов наблюдений РФ  за  120 летний период. Рассмотрена доля компонент,  годовая амплитуда температуры воздуха  и временной ход температуры  сезонных  компонент в среднем по ЕТР, РФ, Сибири и Дальнего Востока с учетом прогноза.

Ключевые слова: Амплитуда, доля, климат, компонента, температура

Sedov V.E.

PhD in Geography, Senior  researcher.

SHARE AND ANNUAL AMPLITUDE OF ADDITIVE SEASONAL COMPONENT SURFACE AIR TEMPERATURE TIME SERIES  30  OBSERVATION POIONTS OF RUSSIA 120-YEAR PERIOD

Abstract

Analyses the annual cycle in surface air temperature based on additive seasonal component, recovered from a 30-minute time series observation of Russia for 120 years. Considered the share component, annual amplitude of temperature and temporal temperature average component of the ETR, Russia, Siberia and the far East, taking into account the forecast.

Keywords: Amplitude, share, climate, component, temperature.

В Российской Федерации (РФ), как и везде на Земном шаре севернее или южнее тропического пояса,  основным процессом климатической системы (КС)  является, как известно,  12-ти месячный сезонный цикл, обусловленный годовым ходом поступающей инсоляции, зависящий от широты местности.  Разные гидрометеорологические элементы, в том числе атмосферная циркуляция, циклоническая деятельность, тепловой и водный обмен с Мировым океаном и др. в годовом цикле имеют свои особенности.  А такой элемент, как температура воздуха, считается ключевым для характеристики климата. На сезонный цикл накладываются слабо предсказуемые квазициклические процессы, обусловленные различными естественными и антропогенными факторами, детально рассмотренными во многих работах, например [1,2,3,4, 6,12],  гидрометеорологические параметры разных сезонов не такие, как если бы они зависели только от солнечной радиации.

В данной работе годовая сезонная периодичность была восстановлена из исходных временных рядов (ВР) 30-ти станций РФ за 1893-2012 годы, а затем исследовалась  на основе двух взаимозависимых переменных: доли (%) и годовой амплитуды температуры воздуха (0С).

Вычисление переменных  выполнялось  современным методом  сингулярного спектрального анализа (Singular Spectrum Analysis Catepillar SSA, в России  называется “Гусеница”) [10,11], по алгоритму которого разработана довольно эффективная  программа[13], которой автор пользовался ранее при исследовании климата Москвы.[8].

Можно кратко напомнить, что  базовsq алгоритм состоит из двух этапов:

1) разложения и 2) восстановления, включающих, в свою очередь, несколько   шагов.     На первом шаге этапа разложения, называемым «вложение», ВР из   N значений раскладывается в последовательность скользящих векторов (временных рядов) с длиной окна, равной L, в результате чего образуется

K = N – L + 1. векторов вложения, которые становятся столбцами траекторной матрицы Х размером L х K.    На втором шаге – сингулярного разложения траекторной матрицы происходит транспонирование матрицы Х, а затем перемножение их, которое образует новую матрицу.

S = X XT               (1)

В матрице S, с помощью специального алгоритма, находятся собственные числа λ, взятые в неубывающем  порядке,

≥ λ≥  ………λL  ≥  0),               (2)

ортонормированная система собственных векторов  U1,…..,UL   матрицы  S и факторные векторы  V1, ….., Vd,

где              17-08-2015 11-17-57               (3)

i = 1, 2,……,d,  образуя, собственные тройки,  каждая, из которых включает собственный  вектор U, собственное число λ и факторный вектор V. Количество троек равно К. Тогда сингулярное разложение  матрицы Х  может быть записано  как

17-08-2015 11-18-13               (4)

Каждая их матриц Хi    имеет ранг 1, т.е. один собственный вектор, одно собственное число и один факторный вектор (эти три переменные матриц называются условно “тройка”). На этапе разложения производится также вычисление периодограммы    в каждом   скользящем векторе и последующее их усреднение по “К” векторам.

Этап восстановления  состоит из двух шагов: шаг 3 – группировка или отбор троек с целью восстановления заданных аддитивных компонент (составляющих) для их самостоятельного исследования и, если нужно, продолжения (прогноза); – 4-й шаг – это диагональное усреднение значений сингулярной матрицы, заканчивающийся восстановлением ряда на основе отобранных  троек  и вычислением аппроксимирующего  ряда.

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

Сезонная составляющая, как 12-ти месячная периодичность вычисляется в виде суммы константы и шести гармоник из выражения 5

17-08-2015 11-18-53               (5)

k=1

где T(n) значение температуры СK в точке ряда n,

π –число пи, Ck и фк  - произвольные константы. При условии Ck  ≠ 0,  для к = 0….5  ряд   имеет размерность 12.

В качестве исходных данных были использованы ряды средних месячных значений приземной температуры воздуха (в ºС) инструментальных наблюдений 30-ти репрезентативных метеорологических станций РФ за 120 ле  Название станций с их географическими координатами приводится в табл.1. При восстановлении сезонных компонент двумерные матрицы данных автором статьи были транспонированы в одномерные последовательные временные ряды Т(Y), где Y = N x 12.

Для общей характеристики сезонных компонент, отметим следующее. 

В среднем по Земному Шару доля и амплитуда СК определяются широтным кругом с максимумом на полюсе и минимумом на экваторе [5]. Важными факторами изменчивости доли и амплитуды являются удаленность

места наблюдений  от океана (моря), центра действия атмосферы. Влияет  на изменчивость инерционность  КС[10].

Восстановление СК, если она зафиксирована в файле ВР, выполняется для любого значения окна при L ≥ 12. Но для более точного восстановления СК значение L должно быть кратно 12-ти.  Вид СК зависит от величины “L”: чем больше значение L, тем большую сглаженность имеет ряд СК.

В данной работе анализируются средние параметры упомянутых переменных за рассматриваемый период наблюдений, рассчитанные  с окнами   L=N/2.

Значение доли СК (далее просто доля) определялось на шаге группировки до восстановления компоненты. Как показали расчеты,  в зависимости от “L”  доля СК изменяется незначительно. Например,  по станции Москва  для

L =  60, 120. 360, 720  доля, соответственно, равна  77.21,  77.30,  77.62  и 77.86%. т.е. вычисления статистически надежны.

Полученные, с окном L=N/2,   средние  доли СК  по 30 пунктам наблюдений РФ, в порядке возрастания их географической долготы, приводятся в табл.1 для трех временных интервалов: -за весь 120-летний период наблюдений (колонка 3.1), -за первое 30-летие ряда  1893-1922 (к. 4.1) и  -30-летие  за 1983-2012 (к. 5.1). Данные табл.1 показывают, что доли в пунктах наблюдений существенно различается в довольно широком диапазоне. Можно отметить, что  расчетные  данные неплохо согласуются с общими закономерностями, отмеченными выше, т.е. наблюдается уменьшение доли от высоких широт к низким широтам и увеличение с запада на восток. Например, на станциях Онега и Сочи,  расположенных примерно на одинаковой географической долготе,  но  на 20.1о  с разницей в широте, значение доли на станции Онега на  80% больше, чем в Сочи. Правда, в Сочи более “мористый”  климат, за счет чего доля также уменьшилась. Но все равно, это подтверждает правильность расчетов доли с помощью метода “Гусеница”.

Средняя доля  по 30 станциям РФ составляет 77, 7%. В 21 пунктах (18 из них ЕТР)  доля уменьшается, а в 9 пунктах  наоборот, возрастает. В широтном поясе от 50 до 60О с. ш. в среднем градиент доли  равен 5.8%/10Oд. Южнее 50О с.ш. доля СК уменьшается до 15.4% (Сочи). В Западной Сибири доля СК увеличивается до 92,7% (Екатеринбург), достигая максимума на востоке Сибири (в Иркутске.) до 96.8%. Далее на восток  доля уменьшается до 82.8% (Владивосток).  В субарктических широтах Сибири доля  СК также высокая.                        

Рассчитанные доли в первом и последнем 30-летии показывают, что во временном ходе доля СК под влиянием разных факторов изменяется.

За 120-летний период доля убывает на двадцати одной станции, в среднем, на 4%, в том по ЕТР – 4.7%,  Сибири и Дальнего Востока – 2.9%. Но на  семи станциях Сибири и Дальнего Востока доля наоборот возрастает на 2.1%. Различие в географической широте в пределах РФ не всегда является определяющим.  Например, ст. Махачкала  расположена южнее станции Сочи на 0.48 ш., но доля СК в Сочи на 19.5% меньше, чем в Махачкале.

Годовая амплитуда приземной  температуры воздуха является важной характеристикой климата. По годовой амплитуде температуры с помощью  разных  алгоритмов [7] могут вычисляться  индексы “континентальности” климата в пункте наблюдения. Это обстоятельство и учитывалось при выборе переменной для анализа.

СК вычислена,   как разность средних месячных значений температуры июля и января. Рассчитанные значения   амплитуды (далее амплитуды) приведены в табл.1. Данные показывают,  что амплитуда температуры СК, по 30-ти пунктам наблюдений, вносит основной вклад в общий годовой ход параметра. Также как и  доля, амплитуда СК вычислялась за весь период наблюдений, и для определения тенденции – за первое и последнее 30-летие.

Сопоставление годовой амплитуды СК с амплитудой фактических данных (колонке 3.3), показывает, что амплитуда СК является определяющей характеристикой. В среднем отношение амплитуды СК к общей амплитуде  составляет  99.17%.  Но в 12 случаях из 30,  амплитуда СК превышает общую амплитуду. Это означает, что амплитуда СК  не всегда обусловлена значением доли.

Таблица 1 - Значения доли (%) и амплитуды (О С) сезонных компонент по 30 станциям РФ

1-Название метеорологической станции, 2.1- географическая широта и  2.2- долгота в градусах;  3.1-доля СК и 3.2- амплитуда всего ряда  (1893-2012гг, 3.3- отношение амплитуды  СК к амплитуде фактической;  4.1-доля и амплитуда  первого 30-летия (1893- 1922 гг);  5.1- доля и 5.2- амплитуда последнего  30-летия (1983-2012 гг; 6.1- разность значений колонок 4.1 и 5.1;  6.2-разность значений колонок 4.2 и 5.2

17-08-2015 11-21-17

Например, на станциях Калининград и Махачкала величина доли небольшая, а амплитуда СК превышает амплитуду фактических данных. Вероятно,  это объясняется наличием во  временных рядах больших положительных и отрицательных аномалий.

На ЕТР  амплитуды изменяются   от 16.7 О С (Сочи) до 33.2 ОС (Волгоград). На станциях Сибири и Дальнего Востока, за исключением Корсаково,  амплитуды  превышают  30 ОС с максимумом 64.7 ОС в Верхоянске.

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

30-летия к последнему 30-летию от -0.2 до -3.7 ОС,  Это означает, в этих регионах степень “континентальности” за рассматриваемый период в среднем уменьшилась.  Исключение составляют лишь три пункта наблюдения (пункты 6,7,19), где годовая амплитуда температуры СК увеличивается и, следовательно, в этих регионах уровень климатической континентальности возрос.

Для общей характеристики временного хода температуры СК на рис.1 приводятся временные ряды за июль и январь в среднем по ЕТР, РФ, Сибири и Дальнего Востока с учетом прогноза на 120 лет. Чтобы исключить высокочастотные климатические квазициклы, ряды рассчитаны для L= 120.

17-08-2015 11-23-07

Рис. 1 – Временные ряды СК за июль и январь

По ЕТР (1 и 4), РФ (2 и 5), Сибири и Дальнего востока (3 и 6)

В среднем  за период 1940-1992гг.,как видно на рис.1, температура сезонных компонент в июле падает, в январе растет т.е. годовая амплитуда CК уменьшается и, следовательно, степень “континентальности”, этом случае, убывает. В среднем  за период 1940-1992гг.,как видно на рис.1, температура сезонных компонент в июле падает, в январе растет т.е. годовая амплитуда CК уменьшается и, следовательно, степень “континентальности”, этом случае, убывает.

Далее, в отличие от среднего значения, годовая амплитуда температуры сезонной компонентs увеличивается, т.е. степень “континентальности” стала возрастать.  Этот процесс согласно прогнозу будет продолжаться примерно до 2032 года.

В заключение отметим, восстановление сезонной компоненты в гидрометеорологических временных рядах, в особенности приземной температуры воздуха, как ключевого элемента, весьма важно для дополнительной характеристики климата, в том числе для определения тенденции степени “континентальности”   климата региона или даже “геоклимата”.

Сезонная компонента  годового климатического цикла под влиянием естественных и антропогенных факторов, так же, как и климатический процесс в целом, подвержена различной степени изменчивости, а наибольшую долю в эту изменчивость вносят естественные (внутренние и внешние) факторы. Вклад радиационных фактов в сезонном цикле, в зависимости от региона существенно различается. Наибольшая доля    инсоляции в сезонном цикле наблюдается в центральной части Сибирского региона. Судя по прогнозу, настоящий период в РФ  увеличивается годовая амплитуда температуры СК. А это свидетельствует о возрастании степени “континентальности” и этот процесс может продолжаться примерно до 2030 года.

Литература

  1. Будыко М.И. Причины естественной изменчивости климата. –В книге Достижения1 в области гидрометеорологии и контроля природной среды. Л. Гидрометеоиздат, 1987. С 201-229.
  2. Груза Г.В., Ранькова Э.Я. Оценка предстоящих изменений климата на Территории Российской  Федерации. Метеорология и гидрология, 2009, № 11 с.15-29.
  3. Гусакова М.А, Л.Н. Карлин. Оценка вклада парниковых газов, водяного пара и облачности в изменении глобальной приповерхностной температуры. Метеорология и Гидрология, -2013, № 3
  4. Зверев И.И., С.К.Гулёв. Сезонность и нестационарность изменчивости Европейского климата в двадцатом веке. Доклады Академии Наук, 2007, том 416, № 5, с.1-4
  5. Кислов А.В. Климат в прошлом, настоящем и будущем. -М.: Наука/Интерпериодика. 2010. -352 с.
  6. Монин А.С. Д.М. Сонечкин. Колебания климата. –М. Наука. 2005.-192с
  7. Российский Гидрометеорологический Энциклопедический словарь/Под ред. А.И. Бедрицкого. -СПб., Москва, Летний сад, 2008.-Т1 А-И.-336с.
  8. Седов В.Е. О климатических колебаниях и тенденции климата современной Москвы. Метеорология и Гидрология. -2012, № 8.
  9. Шерстюков Б.Г. Региональные и зональные закономерности изменений современного климата. –Обнинск: ГУ ВНИИГМИ-МЦД. 208. -246 с.
  10. Elsner J., Tsonis A. Singular Spectrum Analysis. A New Tool in Time Series Analysis. New York: Plenum Press, 1996. 163 p
  11. Golyandina N., Nekrutkin V., Zhigljavsky A. Analysis of Time Series Structure: SSA and Related Techniques. Boca Raton: Chapman & Hall/CRC. 2001. 305 p.
  12. Oliver J.E.,Hidore J.J. Climatology. An atmomospheric science. Prentice Hall. New Jersey.2002.-410 p.
  13. gistatgoup.com/gus/

References

  1. Budyko M.I. Prichiny estestvennoj izmenchivosti klimata. –V knige Dostizhenija1 v oblasti gidrometeorologii i kontrolja prirodnoj sredy. L. Gidrometeoizdat, 1987. S 201-229.
  2. Gruza G.V., Ran'kova Je.Ja. Ocenka predstojashhih izmenenij klimata na Territorii Rossijskoj  Federacii. Meteorologija i gidrologija, 2009, № 11 s.15-29.
  3. Gusakova M.A, L.N. Karlin. Ocenka vklada parnikovyh gazov, vodjanogo para i oblachnosti v izmenenii global'noj pripoverhnostnoj temperatury. Meteorologija i Gidrologija, -2013, № 3
  4. Zverev I.I., S.K.Guljov. Sezonnost' i nestacionarnost' izmenchivosti Evropejskogo klimata v dvadcatom veke. Doklady Akademii Nauk, 2007, tom 416, № 5, s.1-4
  5. Kislov A.V. Klimat v proshlom, nastojashhem i budushhem. -M.: Nauka/Interperiodika. 2010. -352 s.
  6. Monin A.S. D.M. Sonechkin. Kolebanija klimata. –M. Nauka. 2005.-192s
  7. Rossijskij Gidrometeorologicheskij Jenciklopedicheskij slovar'/Pod red. A.I. Bedrickogo. -SPb., Moskva, Letnij sad, 2008.-T1 A-I.-336s.
  8. Sedov V.E. O klimaticheskih kolebanijah i tendencii klimata sovremennoj Moskvy. Meteorologija i Gidrologija. -2012, № 8.
  9. Sherstjukov B.G. Regional'nye i zonal'nye zakonomernosti izmenenij sovremennogo klimata. –Obninsk: GU VNIIGMI-MCD. 208. -246 s.
  10. Elsner J., Tsonis A. Singular Spectrum Analysis. A New Tool in Time Series Analysis. New York: Plenum Press, 1996. 163 p
  11. Golyandina N., Nekrutkin V., Zhigljavsky A. Analysis of Time Series Structure: SSA and Related Techniques. Boca Raton: Chapman & Hall/CRC. 2001. 305 p.
  12. Oliver J.E.,Hidore J.J. Climatology. An atmomospheric science. Prentice Hall. New Jersey.2002.-410 p.
  13. www.gistatgoup.com/gus/