CC BY
100
ГЕОГРАФИЯ
УДК 551.590.2
ПЕРИОДИЧЕСКИЕ ИЗМЕНЕНИЯ ПРИЗЕМНОЙ ТЕМПЕРАТУРЫ ВОЗДУХА, ВЫЗВАННЫЕ ВЛИЯНИЕМ КОСМИЧЕСКИХ ФАКТОРОВ
М.Б. Богданов, А.В. Катрущенко
Саратовский государственный университет, кафедра метеорологии и климатологии E-mail: [email protected]
Проведен анализ спектра мощности временного ряда среднесуточных значений приземной температуры воздуха в г.Саратове на интервале времени с 1945 по 1999 г. В спектре ряда обнаружены две приливные гармоники, гармоника с периодом синодического месяца, связанная с изменением инсоляции при движении системы Земля - Луна вокруг общего центра масс, а также гармоника, соответствующая среднему периоду обращения активных областей Солнца. Проявлений 11-летнего цикла солнечной активности и периодических изменений инсоляции, вызванных планетными возмущениями орбиты Земли, не обнаружено. Ключевые слова: воздух, приземная температура, г. Саратов, гармоники, космические факторы.
Periodical Changes in the Surface air Temperature, Caused by the Influence of the Space Factors
M.B. Bogdanov, A.V. Katrushchenko
We carried out the analysis of the power spectrum of the time series of the daily mean values of the surface air temperature in Saratov city (Russia) in the time interval from 1945 through 1999. In the power spectrum are discovered two tidal harmonics, the harmonic with the period of synodic month, connected with a change in the insolation during the motion of system the Earth - Moon around the common center of masses, and also harmonic, which corresponds to the average period of rotation of the active regions of the Sun. The manifestations of the 11-year solar cycle and changes in the insolation, caused by the planetary perturbations of the Earth's orbit, are not discovered. Key words: ground air temperature, Saratov, harmonic components, space factors.
Введение
В последние годы многими представителями наук о Земле исследуется проблема влияния космических факторов на процессы, протекающие в глубинах Земли, ее атмосфере и гидросфере. В число этих факторов входят разнообразные проявления солнечной активности [1], особенности движения Земли в Солнечной системе, ее вращения вокруг собственной оси и движения полюса [2, 3], а также приливные воздействия Луны и Солнца, способные вызывать изменения океанической циркуляции и влиять на характеристики регионального и глобального климата [4-6].
В ходе предыдущих исследований мы изучили короткопериоди-ческие изменения инсоляции, вызванные планетными возмущениями орбиты Земли [7], совместным действием этих возмущений и вариацией солнечной постоянной [8], а также движением системы Земля - Луна вокруг общего центра масс [9]. Было показано, что наиболее заметными короткопериодическими гармониками изменения инсоляции является гармоника L0 с периодом синодического месяца (промежуток времени между одинаковыми фазами Луны - 29d.53) и амплитудой 82.54 ± ± 0.07 мВт/м2, а также гармоника J0, соответствующая синодическому периоду (промежуток времени между одноименными конфигурациями
© М.Б. Богданов, А.В. Катрущенко, 2008
планеты) Юпитера 399d.0 с амплитудой 44.06 ± ± 3.19 мВт/м2, и первая гармоника синодического периода Венеры V1 с периодом 291d.9 и амплитудой 42.72 ± 1.12 мВт/м2. Представляет большой интерес попытка обнаружить проявление этих гармоник во временных рядах измерений метеорологических величин.
Целью настоящей работы является поиск возможных периодических составляющих, связанных с изменениями инсоляции и влиянием других космических факторов, в ряде среднесуточных значений приземной температуры воздуха в Саратове, измеренных в период с 1945 по 1999 г.
Исходные данные и их предварительная обработка
Необходимо отметить, что поставленная нами задача достаточно сложна. Это связано не только с малыми амплитудами гармоник инсоляции, но и с тем, что они действуют на фоне больших изменений температуры, вызванных сменой сезонов года. Кроме того, приземная температура воздуха определяется главным образом процессами циркуляции атмосферы и сильно зависит от состояния облачности. Однако периодический характер гармоник позволяет обнаружить их проявление при анализе достаточно продолжительных временных рядов.
Мы проанализировали ряд среднесуточных значений приземной температуры воздуха, измеренных на метеостанции 34172 Саратов (51 °34' N, 46°02' E, H = 156 м) с 1 октября 1945 г. по 31 августа 1999 г. Информация была получена из ВНИИГМИ - МЦД в рамках проекта European Climate Assessment [10] и в настоящее время доступна по сети Интернет (http:eca.knmi.nl). В нескольких случаях отсутствующие значения были рассчитаны с помощью линей-
lgS(/)
2,0
1,0
0,0
-1,0
-2,0
ной интерполяции данных. Полное число отсчетов ряда составляет 19693.
Для устранения влияния сезонных изменений температуры из ряда был вычтен вклад синусоиды с периодом, равным тропическому году (промежутку времени между двумя последовательными прохождениями Солнца точки весеннего равноденствия - 365.ё2422). Амплитуда 16.680 К и фаза 2.887 рад. этой синусоиды определялись по способу наименьших квадратов. Аналогичным образом были определены амплитуда 1.028 К и фаза 6.109 рад. первой гармоники тропического года, имеющей период 182.ё6211, которая также была вычтена из ряда среднесуточных значений температуры.
Стандартное отклонение полученного ряда характеризует изменение среднесуточных значений температуры в Саратове, не связанное с сезонными вариациями, и равно 4.93 К.
Анализ спектра мощности временного ряда
Нами использовался классический вариант цифрового спектрального анализа с расчетом спектра мощности Б/ путем фурье-преобразования автокорреляционной функции [11]. При выбранном максимальном сдвиге автокорреляции в 5000 отсчетов для использованного корреляционного окна Ханна ширина спектрального окна составляет Д/ = 2.667 х 10-4 ё-1. Число степеней свободы спектральной оценки равно 11, а длина 90% - го логарифмического доверительного интервала - 0.634. График логарифма спектра мощности в диапазоне частот / от нуля до максимальной частоты Найквиста 0.5 ё-1, построенный по 5001 расчетной точке, приведен на рис. 1. Вертикальным отрезком на этом рисунке показан логарифмический доверительный интервал.
0,0
0,1
0,2
0,3
0,4
14
л 0,5
Рис. 1. График логарифма спектра мощности временного ряда изменения приземной температуры воздуха во всем диапазоне частот
В целом характер спектра мощности соответствует чисто случайному процессу, близкому к «красному» шуму, при котором спектральная плотность возрастает с уменьшением частоты. Единственный хорошо заметный пик на частоте / = 0.00821 ё-1, соответствующей периоду 121. ё7474, представляет собой вторую гармонику тропического года, не исключавшуюся при предварительной обработке ряда.
Участок спектра мощности в диапазоне частот от нуля до / = 0.01 ё-1, содержащий гармоники с периодами Т, превышающими 100ё, приведен на рис. 2. Вертикальным отрезком на этом рисунке показан логарифмический доверительный интервал, а горизонтальным - ширина спектрального окна, характеризующая разрешающую способность по частоте. Единственный статистически значимый пик соответствует второй гармонике тропического года. Анализ этого участка спектра свидетельствует об отсутствии пиков на частотах, соответствующих 11-летнему циклу солнечной активности, синоди-
ческому периоду Юпитера /0 и первой гармонике синодического периода Венеры У1. Если для обнаружения проявления солнечного цикла длина нашего ряда может считаться недостаточной, то отсутствие гармоник планетных возмущений свидетельствует либо о том, что их влияние слишком мало, либо оно в значительной степени подавляется шумовыми флуктуациями.
На рис. 2 отсутствуют также пики, соответствующие квазидвухлетнему циклу (РВО, Т = = 2.40 года) и чандлеровскому периоду движения полюса Земли (Т = 1.19 года). Квазидвухлетний цикл был обнаружен по изменениям скорости ветра в экваториальной стратосфере и часто проявляется в различных геофизических явлениях [3]. С чандлеровским периодом связан так называемый полюсный прилив, влияющий также на изменение приземного атмосферного давления [12]. Данный результат свидетельствует о слабости влияния этих факторов на температуру нижних слоев атмосферы.
2,5
2,0
1,5
1,0
0,5
0,000
0,002
0,004
0,006
0,008
М"1
0,010
Рис. 2. График логарифма спектра мощности в диапазоне гармоник с периодами, превышающими 100ё. Единственный значимый пик соответствует второй гармонике тропического года
Область спектра в диапазоне частот f от 0.03 ё-1 до 0.04 ё-1, соответствующих периодам около одного месяца, приведена на рис. 3. Совмещая середину доверительного интервала с основанием пиков в спектре мощности можно убедиться, что ряд наблюдаемых гармоник являются статистически значимыми. Если принять, что погрешность оценки частоты гармоники т соответствующей максимуму пика, равна ширине спектрального окна А/ то погрешность определения соответствующего периода Т = И^ будет равна АТ = Т^А/ Считая, что найденные периоды гармоник должны отличаться от точных значений
не более чем на величину АТ, мы можем отождествить в данной области спектра три периодических составляющих: гармонику синодического месяца Ь0 (Т = 29ё.53), приливную гармонику Мт (Т = 27ё.55) и гармонику, связанную с солнечной активностью, Т @ (Т = 27ё.0). Обозначения этих гармоник показаны у вершин соответствующих пиков на рис. 3.
Данные об отождествленных гармониках приведены в таблице. В первом столбце дано обозначение гармоники, во втором - частота т соответствующая максимуму спектрального пика, в третьем - найденное значение периода Т
Рис. 3. График логарифма спектра мощности в диапазоне гармоник с периодами около месяца. Отмечены пики, соответствующие гармонике синодического месяца £0, приливной гармонике Мт, а также гармонике Т д, совпадающей со средним периодом обращения активных областей Солнца
с его стандартным уклонением ДТ, а в четвертом - амплитуда, измеренная в кельвинах. В качестве оценки амплитуды гармоники нами бралась величина
А(/т) = ^2 А/(5(/т ) -%), где а - среднеквадратичное отклонение временного ряда, равное 4.93 К; Д/- ширина спектрального окна; Б(/т) - значение спектральной плотности в максимуме пика; а £0 - величина спектральной
плотности у основания пика, оцениваемая как полусумма значений Б(/), измеренных слева и справа от пика. Известно, что при цифровом спектральном анализе частоты гармоник определяются гораздо точнее, чем их амплитуды. С учетом величины доверительного интервала спектральной оценки можно утверждать, что с вероятностью 90% амплитуды гармоник должны отличаться от найденных значений не более чем в два раза.
Характеристики гармоник, найденных в спектре мощности временного ряда изменения приземной температуры воздуха
Гармоника m d-1 T, d A f), К
L0 0.03364 29.73 ± 0.24 0.29
M 0.03616 27.65 ± 0.20 0.36
T ® 0.03700 27.03 ± 0.19 0.24
Mf 0.07332 13.64 ± 0.05 0.16
Область спектра в диапазоне частот / от 0.07 ё-1 до 0.08 ё-1, соответствующих периодам около полумесяца приведена на рис. 4. Данная область содержит пик с частотой 0.07332 ё-1, соответствующий приливной гармонике М (Т = = 13ё.66). Характеристики этой гармоники также приведены в таблице.
Обсуждение результатов
Несмотря на то что не все спектральные пики обнаруженных гармоник заметно выходят за границу принятого доверительного интервала, хоро-
шее совпадение найденных периодов с точными значениями служит дополнительным аргументом в пользу их достоверности.
Очевидным объяснением присутствия во флуктуациях приземной температуры воздуха гармоники £0 с периодом синодического месяца является изменение инсоляции при движении системы Земля - Луна вокруг общего центра масс с амплитудой 82.54 ± 0.07 мВт/м2 [9]. Существует еще один космический фактор, способный повлиять на температуру атмосферы, который действует синхронно с данным изменением инсоляции. Речь идет об излучении Луны, зависящем от ее фазы и
Рис. 4. График логарифма спектра мощности в диапазоне гармоник с периодами около полумесяца. Отмечен пик, соответствующий приливной гармонике Ы^
достигающего максимума в момент полнолуния. Однако влияние этого фактора достаточно мало. В видимом диапазоне спектра (фотометрическая полоса V) звездные величины полной Луны и Солнца равны соответственно VL = - 12.73 и V @ = - 26.74 [13], что дает отношение потоков излучения от этих небесных тел 0.25 х 10-5. Следует однако учесть, что альбедо Луны мало (« 0.07) и основная часть освещающего ее солнечного излучения поглощается поверхностью нашего спутника. Нагретая поверхность Луны излучает в ИК-диапазоне спектра. Поэтому отношение болометрических потоков должно быть больше. По данным радиометрических измерений [14], для среднего полнолуния оно равно 0.96 х 10-5, что с учетом величины солнечной постоянной 1366.22 Вт/м2 дает поток от Луны 13 мВт/м2. В фазе новолуния поток от Луны очень мал. Поэтому в качестве верхнего предела амплитуды изменения потока ее излучения можно принять значение 6.5 мВт/м2, что более чем на порядок меньше амплитуды изменения инсоляции при движении системы Земля - Луна вокруг общего центра масс.
Присутствие в спектре мощности нашего временного ряда гармоник лунно-солнечных приливов Ыт и Ы ^ не является неожиданным. Атмосферные приливы вызывают изменение давления, что должно сопровождаться и изменением температуры. Теория этих приливов достаточно хорошо разработана. Однако в известной монографии [15] и более поздних исследованиях [16] рассматривались только наиболее мощные гармоники гравитационных и термических приливов с суточным и меньшими периодами. Было бы ин-
тересно провести аналогичные расчеты амплитуд приземных температур воздуха, связанных с месячной и полумесячной приливными гармониками для сравнения с нашими результатами.
Вызванная солнечной активностью гармоника Т @ имеет период Т = 27а.0, совпадающий со средним периодом обращения активных областей Солнца [17]. Дать однозначное объяснение появлению этой гармоники в ряде приземных температур воздуха достаточно трудно. С одной стороны, здесь, возможно, проявляет себя блокировка излучения фотосферы, вызванная наличием групп пятен. С другой - в районе этих активных областей происходят хромосферные вспышки, сопровождающиеся резким усилением жесткого УФ-излучения, а иногда и рентгеновского. Излучение этих диапазонов спектра поглощается верхними слоями земной атмосферы, вызывая их заметный нагрев. Однако до настоящего времени не известен сколь-либо эффективный механизм передачи этой энергии к поверхности Земли. Оба указанных фактора солнечной активности могут действовать и совместно.
Заключение
Проведенный анализ спектра мощности временного ряда среднесуточных значений приземной температуры воздуха в Саратове в период с 1945 по 1999 г. позволил достаточно надежно обнаружить присутствие в изменениях температуры периодических составляющих. В спектре ряда найдены две приливные гармоники Ыт и Ы , гармоника с периодом синодического месяца П0, связанная с изменением инсоляции при движении системы
Земля - Луна вокруг общего центра масс, а также гармоника T @, соответствующая среднему периоду обращения активных областей Солнца.
Проявлений 11-летнего цикла солнечной активности и периодических изменений инсоляции, вызванных планетными возмущениями орбиты Земли, нами не обнаружено. В ходе солнечного цикла, наряду с большими вариациями потоков жесткого УФ-излучения и рентгеновского, наблюдается также изменение инсоляции с амплитудой около 1 Вт/м2. Отсутствие проявлений этого цикла можно было бы объяснить сравнительно малой длиной нашего ряда. Однако анализ более продолжительных рядов среднегодовых температур различных пунктов также не дает однозначных свидетельств в пользу связи их изменений с солнечной активностью [17].
Две главные гармоники изменения инсоляции J0 и V1, вызванные планетными возмущениями, имеют периоды около года и вполне могли бы проявить себя. Отсутствие в спектре мощности соответствующих пиков свидетельствует либо о том, что их влияние слишком мало, либо оно в значительной степени подавляется шумовыми флуктуациями. То же самое можно сказать в отношении гармоник, соответствующих квазидвухлетнему циклу и чандле-ровскому периоду движения полюса Земли. Хотя эти циклы часто проявляются в различных геофизических явлениях, отсутствие соответствующих пиков в спектре свидетельствует о слабости их влияния на температуру нижних слоев атмосферы.
Библиографический список
1. Авдюшин С.И., Данилов А.Д. Солнце, погода и климат: сегодняшний взгляд на проблему. (Обзор) // Геомагнетизм и аэрономия. 2000. Т. 40, № 5. С. 3-14.
2. Монин А.С., ШишковЮ.А. Климат как проблема физики // Успехи физ. наук. 2000. Т. 170, № 4. С. 419-445.
3. Сидоренков Н.С. Атмосферные процессы и вращение Земли. СПб.: Гидрометеоиздат, 2002. 200 с.
4. Авсюк Ю.Н. Приливные силы и природные процессы. М.: Объединенный ин-т физики Земли им. О.Ю. Шмидта, 1996. 188 с.
5. Wunsch C. Moon, tides and climate // Nature. 2000. Vol. 405, P. 743-744.
6. Treloar N.C. Luni - solar tidal influences on climate variability // Intern. J. of Climatology. 2002. Vol. 22, № 12. P. 1527-1542.
7. Богданов М.Б., Сурков А.Н. Короткопериодные изменения инсоляции, вызванные планетными возмущениями орбиты Земли // Метеорология и гидрология. 2006. № 1. С. 48-54.
8. Богданов М.Б., Катрущенко А.В., Сурков А.Н. Изменения инсоляции, вызванные планетными возмущениями орбиты Земли и вариацией солнечной постоянной // Изв. Сарат. ун-та. Сер. Науки о Земле. 2006. Т. 6, вып. 1. С. 3-9.
9. Богданов М.Б., Катрущенко А.В. Изменения инсоляции, вызванные влиянием Луны // Изв. Сарат. ун-та. Сер. Науки о Земле. 2008. Т. 8, вып. 1. С. 3-6.
10. Klein Tank A.M.G., Wijngaard J.B., Konnen G.P. et al. Daily dataset of 20th-century surface air temperature and precipitation series for the European Climate Assessment // Intern. J. of Climatology. 2002. Vol. 22. P. 1441-1453.
11. Дженкинс Г., Ваттс Д. Спектральный анализ и его приложения. М.: Мир, 1972. 603 с.
12. МаксимовИ.В. Геофизические силы и воды океана. Л.: Гидрометеоиздат, 1970. 447 с.
13. АлленК.У. Астрофизические величины. М.: Мир, 1977. 446 с.
14. Pettit E. Lunar radiation as related to phase // Astrophys. J. 1935. Vol. 81. P. 17-36.
15. Чепмен С., Линдзен Р. Атмосферные приливы: термические и гравитационные. М.: Мир, 1972. 296 с.
16. ZurekR.W. The form of Newtonian cooling in atmospheric tidal theory // Pure and Applied Geophys. 1985. Vol. 123. P. 921-929.
17. Герман Дж.Р., Голдберг Р.А. Солнце, погода и климат. Л.: Гидрометеоиздат, 1981. 319 с.
УДК 502.3
СПОСОБНОСТЬ АТМОСФЕРЫ РАЗЛИЧНЫХ РАЙОНОВ САРАТОВСКОЙ ОБЛАСТИ К САМООЧИЩЕНИЮ
С.Н. Лапина, Е.А. Полянская, Л.М. Фетисова, H.A. Фетисова
Саратовский государственный университет, кафедра метеорологии и климатологии E-mail: [email protected]
В статье приводятся результаты исследования способности атмосферы к самоочищению в Саратове и Саратовской области. Проанализирована динамика коэффициента самоочищения атмосферы в различных районах Саратовской области в теплый и холодный сезоны за многолетний период и 2000-2004 гг. Формула расчета коэффициента самоочищения атмосферы адап-
тирована к региональным климатическим условиям. Построены карты-схемы районирования территории Саратова по коэффициенту самоочищения атмосферы. Исследование метеорологического потенциала самоочищения атмосферы имеет практическое значение при определении оптимального распределения техногенных нагрузок на окружающую среду.
© С.Н. Лапина, Е.А. Полянская, Л.М. Фетисова, НА Фетисова, 2008
