CC BY
28
Том 149, кн. 4
Естественные науки
2007
УДК 551.501
ОПЫТ АНАЛИЗА МНОГОЛЕТНЕЙ ДИНАМИКИ АНОМАЛИЙ СРЕДНИХ ГОДОВЫХ ТЕМПЕРАТУР ВОЗДУХА В КАЗАНИ С ИСПОЛЬЗОВАНИЕМ НАЧАЛЬНЫХ ПОЛОЖЕНИЙ ТЕОРИИ СЛУЧАЙНЫХ ФУНКЦИЙ
М.А. Верещагин, Ю.П. Переведенцев, Э.П. Наумов, К.М. Шанталинский, В.Д. Тудрий, Ф.В. Гоголь
Аннотация
Многолетняя (1900-2004 гг.) динамика аномалий средних годовых температур воздуха, регистрируемых в метеорологической обсерватории Казанского государственного университета, рассматривается как поведение некоторой системы, проходящей через определенные ее состояния. С использованием начальных положений теории случайных функций исследуются некоторые особенности функционирования этой системы, имеющие научное и прикладное значение.
Введение
Город Казань относится к числу немногих пунктов регулярных метеорологических наблюдений на территории России, где их непрерывная продолжительность (и в их числе измерения температуры воздуха) насчитывает уже почти два столетия. История становления этих наблюдений неразрывно связана с деятельностью метеорологической обсерватории Казанского университета, основанной в 1812 г., вскоре после его открытия (1804 г.).
Попытки наладить систематические измерения температуры воздуха, атмосферного давления, характеристик ветра в Казани предпринимались и много раньше (в 1733, 1805 гг.) [1]. К сожалению, большая часть архива этих наблюдений оказалась утраченной, да и сами наблюдения в те годы не были регулярными.
До наших дней в наиболее полном и качественном виде дошли лишь непрерывные ряды (с 1828 г.) средних месячных и годовых температур воздуха [2], явившиеся в последующем уникальной базой данных для изучения исторических особенностей колебаний термического режима города [3-5 и др.].
Наблюдательная площадка метеорологической обсерватории с момента начала ее функционирования неизменно находилась во дворе главного корпуса университета (центральная часть города), будучи окруженной со всех сторон плотной городской застройкой. Таким образом, по условиям местонахождения метеорологическую обсерваторию следует относить к категории типично городских пунктов наблюдений.
Начиная со середины XIX столетия и позднее в ближайших окрестностях Казани для производства метеорологических наблюдений было открыто еще несколько станций. Тем самым возникли благоприятные предпосылки для получения объективных оценок воздействия городской среды на многолетнюю динамику термического режима города [3, 6, 7].
Наиболее важной особенностью указанной динамики является тенденция длительного потепления климата города [3, 5, 6, 8], которую в ходе последующего изложения нам придется неоднократно упоминать. В связи с этим представляется целесообразным привести небольшую сводку данных, характеризующих процесс потепления климата города. Ограничимся при этом показателями, отражающими лишь динамику средних годовых температур воздуха (СГТВ).
На протяжении последних 176 лет (1828-2004 гг.) в межгодовых изменениях температуры воздуха присутствовала линейная составляющая тренда потепления. Средняя за этот период времени скорость потепления составила 1.39°С/100 лет, что существенно больше соответствующих показателей потепления (за период 1856-1997 гг.) на Северном полушарии (0.38°С/100 лет) и Земле в целом (0.42°С/100 лет) [8].
В итоге за упомянутые 176 лет СГТВ внутри города увеличилась на 2.4°С [3]. Результаты работы [7] показывают, что в указанном повышении СГТВ 58.3% (~1.4°С) объясняется влиянием «городской составляющей» потепления и лишь 41.7% (~1°С) - действием факторов глобального масштаба.
Из результатов работ [3, 8] следует также, что процесс потепления климата города развивался нелинейно (с ускорением). Наиболее высокими темпами потепление шло в последние десятилетия (1941-2004 гг.). В 1995 г. был достигнут абсолютный максимум СГТВ (6.8°С) за всю историю измерений температуры воздуха.
Объект изучения, методика и результаты исследования
Настоящую работу следует рассматривать как дальнейшее развитие традиционного направления научных исследований кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета, связанного с изучением динамики климата города.
Однако, в отличие от предшествующих работ, основное внимание будет направлено на изложение результатов анализа марковских свойств многолетней (1900-2004 гг.) хронологической последовательности аномалий средних готовых температур воздуха (АСГТВ).
Ежегодные значения АСГТВ (At) вычислялись как разности между текущими величинами СГТВ, регистрируемыми в метеорологической обсерватории, и их средней величиной (t = 4.4°С), полученной из осреднения за период 1961-1990 гг. (рис. 1).
Согласно принятого в практике метеорологического информационного обеспечения стандарта [9], АСГТВ при ее величине At < -1°С квалифицировалась как температура воздуха «ниже нормы» (НН), а при |At| < 1 С и At > 1 С -«в норме» (Н) и «выше нормы» (ВН) соответственно. В итоге многолетняя по-
следовательность АСГТВ рассматривалась как поведение некоторой системы, проходящей указанные три состояния.
При самом беглом рассмотрении рис. 1 может сложиться впечатление о хаотичности функционирования рассматриваемой системы. При более внимательном изучении того же рисунка в сильно «зашумленном» поведении системы все же удается обнаружить некоторые устойчивые закономерности. Так, например, видно, что в начале исследуемого периода (1900-1925 гг.) система в состоянии Н наблюдалась наиболее часто, да и сама интенсивность отрицательных АСГТВ была заметно большей, нежели в конце исследуемого периода (1975-2004 гг.). В конце исследуемого периода чаще, чем в его начале, система находилась в состоянии ВН. В этих особенностях функционирования уже усматриваются первые признаки нестационарного характера ее поведения, обусловленные тенденцией длительного потепления климата города.
Существенно важными для последующего анализа являются оценки безусловных вероятностей обнаружения р7 аномалий А( в каждой из трех градаций (НН (7 = 1), Н (7 = 2), ВН (7 = 3)):
Р, = П / N. (1)
Здесь N = 105 - объем изучаемой выборки, п7 - число обнаружений А7 в 1-й градации. Полученные таким образом оценки р7 (рНН = 0.2115, рН = 0.6630, рВн = 0.1250) далее определяются как компоненты базисного вектора вероятностей состояний системы Р(70):
Р('0) =[ Р1С0), Р2С0), РзМ, (2)
где 70 - символ начального момента времени.
Из условий рН > рНН, рН > рВН следует, что рассматриваемая система при всей сложности ее поведения и наличии в ней систематической составляющей (тенденции потепления) обладает еще значительным потенциалом саморегулирования: после выхода термического режима из «нормального» (Н) в его экстремальные (НН или ВН) состояния она чаще всего опять входит в «норму», чтобы затем из этого состояния снова перейти в состояние НН или ВН.
Изучаемую хронологическую последовательность АСГТВ (или поведение выше указанной системы, рис. 1, б) будем одновременно трактовать как реализацию случайной функции Х(7) с дискретным временем (с интервалом дискретизации в один год) и дискретным множеством (НН, Н, ВН) ее значений. Начальное (0 состояние системы задается базисным вектором вероятностей Р(?о).
Существенный интерес может представлять при этом выяснение вопроса о том, в какой градации (из 3-х возможных) следует ожидать чаще (реже) всего значения АСГТВ при числе переходов системы в п лет (п > 1) вне зависимости от того, в какой градации она находилась в начальный момент времени. Очевидно, что ответ на этот вопрос предполагает оценки величин безусловных вероятностей обнаружения АСГТВ в каждой из трех градаций:
р, (п) ^^ (7 = 1,2,3). (3)
N - п к
Здесь т7(п) - число попаданий А( в 7-ю градацию, к - объем выборки, учитываемый при оценке р7(п).
3
2
1
0
о
о
«-1
-2
-3
-4
0 0 0 0 0 0 0 0 0 0 о
0 1 2 3 4 ю со СО О) С 3
с» О) О) О) О) О) О) О) О) О) о
1 1 1 1 1 1 1 1 1 1 см
Рис. 1. Многолетняя динамика АСГТВ (А/, °С) на ст. Казань, университет в ее натуральном (а) и градуированном (б) представлениях и результаты ее сглаживания (скользящее осреднение при т = 31 год)
При фиксированном n
I Pi (n) = 1,0,
(4)
а вектор-строку
P( n) =[ Pi(n), p2( n), Рз(п)]
(5)
следует рассматривать как вектор вероятностей состояний системы при числе ее переходов п.
Оценки частот р(п) выполнялись в интервале 1 < п < 90 (рис. 2) с использованием «метода скольжения». При заданном значении п в качестве начальных ординат функции Х(0 рассматривались 1, 2, ..., (Ы- п - 1), (Ы- п)-й члены многолетнего ряда градуированных величин АСГТВ, а в качестве последующих ординат соответственно (1 + п), (1 + п + 1), ..., (Ы - п), Ы-й члены того же ряда.
С увеличением числа переходов п наблюдается уменьшение объема выборки к, используемого при оценке р(п) (3) и некоторое увеличение средней квад-ратической погрешности
определения р(п). Одновременно увеличивается и ширина зоны неопределенности поведения частот р(п)
и ухудшаются условия их распознавания.
При n = const признаком относительно надежного (при уровне доверия
> 68%) распознавания, например, частот pH(n) и Рнн(п) будут являться неперекрывающиеся зоны sH(n) и sBH(n).
Из поведения эмпирических кривых pi(n) (рис. 2) следует, что на всем интервале переходов системы 1 < n < 90 ординаты pH(n) по отношению к ординатам кривых pBH(n), pHH(n) занимают самое высокое положение, а ординаты pHH(n) (при n > 25) - наиболее низкое положение. Можно также легко убедиться в том, что при всех n > 1 нижняя граница зоны sH(n) нигде не перекрывается верхними границами зон sBH(n), sHH(n). Это позволяет утверждать о том, что, указывая в прогнозах (с их заблаговременностью n > 1 и вне зависимости от начального положения АСГТВ) будущее значение At в градации «норма», в 6075% случаев мы можем рассчитывать на «успех». Но если в тех же прогнозах будущее положение At указывать в классе «ниже нормы» мы рискуем ошибиться не менее, чем в 90% случаев.
При n > 60 возрастает надежность прогнозов, если будущее положение At указывать в классе «выше нормы» (рис. 2). Однако, и в этом случае pH(n) >
> pBH(n). Обратим также внимание на то, что с ростом числа переходов n в поведении ординат pH(n) и pHH(n) наблюдаются устойчивые тенденции к их уменьшению, а в поведении ординат pBH(n) - тенденция к их быстрому увеличению.
(6)
st (n) = Рi (n) ±аг (n)
(7)
Рис. 2. Безусловные вероятности попадания АСГТВ в градации «НН» (а), «Н» (б), «ВН.» (в): 1 - фактические, 2 - виртуальные величины, следующие из воспроизводства многолетней динамики АСГТВ простой однородной цепью Маркова
Это свидетельствует о том, что с ростом п наблюдается систематическое изменение закона распределения АСГТВ, что является верным признаком нестационарности многолетней динамики термического режима города [10-11].
Очевидно, что главным ее источником является уже упоминавшаяся длительная тенденция потепления климата города.
Указанная нестационарность проявляется в изменчивости векторов Р(и), в возникновении «ползучей неоднородности» в исходном ряду А( и, как следствие, в становлении проблемы определения оптимальной длины периода осреднения [12-14] для расчета средних многолетних показателей термического режима.
Прямым следствием систематических изменений безусловных вероятностей р(п) является все возрастающая надежность
/(и) = Рн (и) + Рвн (и) (8)
попаданий А( в укрупненную градацию «норма - выше нормы», особенно заметная при и > 50 (табл. 1).
Табл. 1
Вероятности /(и), %
и, годы 50 60 70 80 90
/(и) 90.7 90.9 94.0 95.8 100
Можно при этом полагать, что анализ функций р^и) может быть полезным для выявления скрытых неоднородностей и в других, помимо рядов АСГТВ, временных рядах.
Опыт показывает [12-14], что временные метеорологические ряды в большинстве заключают в себе в той или иной мере систематическую составляющую (тренд) и, строго говоря, нестационарны. Однако в ряде случаев на определенных интервалах времени они могут рассматриваться как почти стационарные. В связи с этим возникает необходимость сравнения динамики реально наблюдаемого процесса Х(() с динамикой строго стационарного процесса 7(0, имитируемого посредством какого-либо алгоритма.
В нашем случае в качестве эталона стационарного процесса рассматривалась простая однородная цепь Маркова. Условия ее функционирования предполагают такое поведение системы, когда условная вероятность перехода из начального (/) состояния в последующее (/) зависит (наряду с учетом /) лишь от длины интервала времени между указанными состояниями (в нашем случае -от числа переходов и) и не зависит от начала его отсчета [10, 11]. Вектор вероятностей состояний такой системы задается равенством
Е(и) = Р(0-Г и, (9)
где Г = ||ргу|| - матрица перехода системы.
Простая однородная модель Маркова стационарна и эргодична, и для нее можно определить вид предельного вектора вероятностей
Е = [Енн,^,^ ]. (10)
При числе переходов и ^ да компоненты вектора Е определяются при этом из решения системы линейных уравнений
[ ^НН , ^Н , ^ВН ] = [ ^НН , ^Н , ^ВН ]-Г. (11)
Применительно к изучаемому ряду вид матрицы Г определился следующим образом
Г р„ Р12 Р13 ' Г 0.32 0.64 0.04 >
Г= Р21 Р22 Р23 = = 0.17 0.71 0.12 . (12)
V Р31 Р32 Р33 у ч 0.17 0.50 0.33 у
Графический вид функций Fi(n) для строго стационарного процесса представлен на рис. 2. Там же показаны и значения компонент предельного вектора FНН, FН, FВН. Разности ординат функций рг(п) и Fi(n) позволяют составить качественное представление о степени согласия (или отличия) наблюдаемой динамики АСГТВ [Р(п)] с их возможной динамикой |Т(п)], предполагающей отсутствие в ряду А( систематической составляющей.
В качестве объективной меры сходства между двумя видами динамики теплового режима может служить комплексный показатель сходства векторов
С(п) = --Р(П)ТГ(П)-"1/2 + (1 + [[) - Г(п)]Т [Р(п) - Щп)]] , (13)
(Р(п)т Р(п) |>(п)т Е(п)])
где Р(п)Т - транспонированный вектор Р(п), п > 1.
При полной аналогии векторов Р(п) и Е(п) - С(п) = 2.0, при полной контраналогии С(п) = -1.0. Вид функции С(п) воспроизводится на рис. 3.
Из совместного рассмотрения рис. 2, 3 следует, что на интервале времени 1 < п < 15(20) ординаты р;(п) и F1(n) почти совпадают, а функция С(п) близка к своему верхнему пределу. Таким образом, при числе переходов п < 15 20 динамика градуированных АСГТВ почти стационарна и достаточно полно воспроизводится простой однородной цепью Маркова. Указанное обстоятельство, очевидно, объясняется тем, что на коротких интервалах времени влияние тренда потепления на колебания А( еще почти не проявляется. Отсюда следует, что на фиксированных отрезках времени продолжительностью до 15 20 лет средние показатели термического режима обладают значительной устойчивостью. Указанный период времени в связи с этим следует рассматривать как наиболее оптимальный (для условий Казани) для получения упомянутых средних величин.
При числе переходов п > 15 20 согласие между наблюдаемой [Р(п)] и теоретически возможной [Р(п)] динамикой А( утрачивается: при 15 < п < 40 -сравнительно медленно, при 40 < п < 65 - заметно быстрее, а при п > 60 - наиболее быстро (рис. 3).
Таким образом, на больших интервалах времени динамика АСГТВ, рассматриваемая условно как реализация случайной функции Х(0, много сложнее динамики рассматриваемой цепи Маркова. Последняя предполагает [10, 11], что закон распределения ординаты Х(^), вполне определяется значением ординаты Х(^-_1), зависит от длины промежутка времени ^ - ^ (числа переходов п) и не зависит от значений предшествующих ординат Х^^), Х(^г-з) и др., т. е. от более ранней предыстории процесса. Последнее противоречит ранее установленному присутствию в Х(0 систематической составляющей.
п, годы
Рис. 3. Динамика комплексного показателя сходства векторов С(п)
Поэтому для описания поведения АСГТВ на длительных интервалах времени можно рассчитывать на наиболее полный успех, привлекая для этого более сложные модели (аппарат марковских процессов высоких порядков). Так, например, в работе [15] для экстраполяции сглаженной динамики СГТВ проверялась возможность использования для этой цели стационарных марковских моделей 10 15-го порядков. При этом были получены положительные результаты.
Из всего сказанного следует также, что при числе переходов п > 15 20 используемая последовательность АСГТВ неэргодична, что приводит, в частности, к неустойчивости средних показателей термического режима и бессмысленности попыток их уточнения путем удлинения имеющихся рядов.
Присутствие в многолетней динамике А( систематической составляющей влечет за собою некоторую предсказуемость их поведения [15]. Об этом свидетельствует, в частности, данные о динамике условных вероятностей р^(п) переходов At из начальной () в последующую (/) градацию (табл. 2). Здесь же наряду с эмпирическими оценками р^(п) показаны также (в круглых скобках) вероятности Fy{n), следующие из результатов имитации динамики А( простой однородной цепью Маркова. Вероятности F¡¡(n) вычислялись с использованием матрицы переходов Г:
(п)|| =ГП. (14)
При п ^ да матрица Гп стремится к предельной матрице с одинаковыми строками [10, 11]. Это свойство начинает проявляться при п > 5 и уже при п = 10 мы имеем (с небольшими округлениями)
(0.200 0.668 0.132^
^ (п) =
0.200 0.668 0.132 0.200 0.668 0.132
(15)
Табл. 2
Условные вероятности р,/п) и Еу(п) (в круглых скобках), %
Градации АСГТВ Ри(п) и Fil(n) при различных п (годы)
начальная, / последующая, ] 1 5 10 20 30
НН НН 31.8 (31.0) 13.6 (20.0) 18.2 (20.0) 14.3 (20.0) 10.0 (20.0)
Н 63.6 (64.0) 77.2 (66.9) 81.8 (66.8) 81.0 (66.8) 80.0 (66.8)
ВН 4.6 (4.0) 9.2 (13.1) 0 (13.2) 4.7 (13.2) 10.0 (13.2)
Н НН 17.4 (17.0) 22.4 (19.9) 20.6 (20.0) 15.8 (20.0) 11.6 (20.0)
Н 71.0 (71.0) 64.2 (66.8) 63.5 (66.8) 68.4 (66.8) 69.2 (66.8)
ВН 11.6 (12.0) 13.4 (13.1) 15.9 (13.2) 15.8 (13.2) 19.2 (13.2)
ВН НН 16.7 (17.0) 0 (19.9) 0 (20.0) 16.7 (20.0) 50.0 (20.0)
Н 50.0 (50.0) 80.0 (66.8) 66.7 (66.8) 33.3 (66.8) 50.0 (66.8)
ВН 33.3 (33.0) 20.0 (13.2) 33.3 (13.2) 50.0 (13.2) 0.0 (13.2)
Из сравнения частот р^(п) и Fi](n) (табл. 2) следует, что увеличение предсказуемости поведения А/ за счет присутствия в многолетнем ряду АСГТВ тенденции потепления наблюдается при этом не во всех возможных случаях, а лишь в определенных.
Наибольшее увеличение предсказуемости (на 10 20%) достигается, в первую очередь, тогда, когда начальное значение А/ относится к градации НН, а в прогнозах при их заблаговременности п > 5 будущее положение А/ указывается в градации Н.
Надежность таких прогнозов значительна: 77.2 <р^(п) < 81.8%. Заметное увеличение предсказуемости наблюдается также в прогнозах, когда будущее положение предиктанта (А/) указывается в градации ВН (до 7 27%) и Н (до 13.2% при п = 5) при начальном положении А/ в градации ВН и в некоторых других случаях.
Выводы
1. Многолетняя динамика АСГТВ в Казани, условно рассматриваемая как поведение некоторой системы, на больших интервалах времени (более 15^20 лет) нестационарна и существенно отличается от поведения простой однородной цепи Маркова. Основным источником указанной нестационарности является присутствие в многолетнем ряду А/ длительной тенденции потепления. Отмеченные обстоятельства влекут за собою неустойчивость средних показателей термического режима, бессмысленность попыток их уточнения за счет удлинения соответствующих рядов. А в отдельных случаях это приводит к заметному возрастанию предсказуемости его поведения.
2. На интервалах времени, не больших 15 20 лет, влияние тренда потепления практически не проявляется, и многолетняя динамика АСГТВ достаточно полно воспроизводится поведением простой однородной цепи Маркова. В рамках указанного периода времени средние показатели термического режима обнаруживают значительную устойчивость, а сам промежуток времени 15 20 лет для условий Казани можно рассматривать как наиболее оптимальный для их определения.
Summary
M.A. Vereschagin, Y.P. Perevedentsev, E.P. Naumov, K.M. Shantalinsky, V.D. Tudry, F.V. Gogol. Experience of the analysis of long-term dynamics averages annual air temperatures anomalies in Kazan using initial positions of the theory of stochastic functions.
Long-term (1900-2004) dynamics of averages annual anomalies of air temperatures, registered in a meteorological observatory of the Kazan State University, is considered as behaviour of some system, which are passing through its certain conditions. Using initial positions of the theory of stochastic functions some features of this system functioning having scientific and applied importance are investigated.
Литература
1. Переведенцев Ю.П. Метеорология в Казанском университете: становление, развитие, основные достижения. - Казань: Изд-во Казан. ун-та, 2001. - 128 с.
2. Климатологический справочник СССР. Вып. 12. Ч. 1. Метеорологические данные за отдельные годы. Температура воздуха. - Л.: Гидрометеоиздат, 1989. - 431 с.
3. Верещагин М.А. Переведенцев Ю.П., Шанталинский К.М. О вероятных изменениях средней годовой температуры воздуха на востоке Русской равнины в конце текущего и в начале XXI столетия // Метеорология и гидрология. - 1996. - № 8. -С. 38-44.
4. Кушников П.Д. Об изменении температуры воздуха в Казани за 97 лет (18281925 гг.) // Учен. зап. Казан. ун-та. - 1925. - Т. 85, Кн. 1. - С. 68-69.
5. СмоляковП.Т. Температурный режим г. Казани // Учен. зап. Казан. ун-та. - 1956. -Т. 116, Кн. 5. - С. 232-235.
6. Колобов Н.В., Верещагин М.А., Переведенцев Ю.П., Шанталинский К.М. Оценка влияния роста Казани на изменения термического режима внутри города // Тр. Зап. СибНИИ. - 1983. - Вып. 57. - С. 37-41.
7. Переведенцев Ю.П., Верещагин М.А., Шанталинский К.М. О многолетних колебаниях температуры воздуха по данным метеорологической обсерватории Казанского университета // Метеорология и гидрология. - 1994. - № 7. - С. 59-67.
8. Переведенцев Ю.П., Верещагин М.А., Шанталинский К.М., Тудрий В.Д. Современное глобальные и региональные изменения окружающей среды и климата. - Казань: УНИПРЕСС, 1999. - 96 с.
9. Багров Н.А., Кондратович К.В., Педь Д.А., Угрюмов А.И. Долгосрочные метеорологические прогнозы. - Л.: Гидрометеоиздат, 1985. - 248 с.
10. Вентцель Е.С. Теория вероятностей. - М.: Физматгиз, 1982. - 564 с.
11. Иванова В.М., Калинина В.Н., Нешумова Л.А., Решетникова И.О. Математическая статистика. - М.: Высш. шк., 1975. - 398 с.
12. Кобышева Н.В., Наровлянский Г.Я. Климатологическая обработка метеорологической информации. - Л.: Гидрометеоиздат, 1978. - 295 с.
13. Дроздов О.А., Васильев В.А., Кобышева Н.В., Раевский А.Н., Смекалова Л.К. Климатология. - Л.: Гидрометеоиздат, 1989. - 568 с.
14. Рубинштейн Е.С. К вопросу о периоде осреднения в климатологии // Тр. ГГО. -1965. - Вып. 181. - С. 46-51.
15. Верещагин М.А., Переведенцев Ю.П., Наумов Э.П., Шанталинский К.М., Гоголь Ф.В. Многолетние изменения температуры воздуха и атмосферных осадков в Казани и их некоторые проявления // Учен. зап. Казан. ун-та. Сер. Естеств. науки. -2005. - Т. 147, Кн. 3. - С. 151-166.
Поступила в редакцию 12.03.07
Верещагин Михаил Алексеевич - кандидат географических наук, доцент кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета.
Переведенцев Юрий Петрович - доктор географических наук, профессор, заведующий кафедрой метеорологии, климатологии и экологии Казанского государственного университета.
E-mail: [email protected]
Наумов Эдуард Петрович - кандидат географических наук, доцент кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета.
Шанталинский Константин Михайлович - кандидат географических наук, доцент кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета.
E-mail: [email protected]
Тудрий Вадим Дмитриевич - кандидат географических наук, доцент кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета.
E-mail: [email protected]
Гоголь Феликс Витальевич - ассистент кафедры метеорологии, климатологии и экологии атмосферы Казанского государственного университета.
E-mail: [email protected]
