Типологизация объектов циркумполярного вихря Северного полушария

Тип работы:
Реферат
Предмет:
Общие и комплексные проблемы естественных и точных наук


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 551. 513
ТИПОЛОГИЗАЦИЯ ОБЪЕКТОВ ЦИРКУМПОЛЯРНОГО ВИХРЯ СЕВЕРНОГО ПОЛУШАРИЯ
Шипко Ю. В., Шувакин Е. В.
ФГВКОУ ВПО Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н. Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, (394 064, Воронеж, ул. Старых Большевиков, 54А), e-mail: argentina_85@mail. ru_
На базе моделирования циркумполярных образований Северного полушария на изобарическом уровне 500 гПа для января проведен многомерный статистический анализ, позволивший выявить три кластера объектов. Учитывались следующие характеристики базовой модели циркумполярного вихря: положение центра условных масс- площадь, масса и плотность вихря- ориентация оси эллипса рассеяния масс, главные средние квадратические отклонения и коэффициент сжатия эллипса. Использовались срочные данные реанализа параметров атмосферы NCEP/DOEAMIP-II за период 1979—2013 гг. Построена факторная модель циркумполярного вихря, дана содержательная интерпретация факторов. Объекты выделенных кластеров геометрически проявляются на факторных плоскостях. Дана оценка статистического различия распределений факторов в кластерах. В результате проведенного дискриминант-ного анализа построены функции, дискриминирующие объекты выделенных кластеров, определены факторы, оказывающие наибольшее влияние на дискриминацию между группами.
Ключевые слова: циркумполярный вихрь Северного полушария, кластеризация, факторная модель, дискрими-нантный анализ.
TYPOLOGIZATION OF CIRCUMPOLAR VORTEX OBJECTS OF THE NORTHERN HEMISPHERE
Shipko Y.V., Shuvakin E.V.
Federal State Official Military Educational Institution of Higher Professional Education Military Educational-Research Centre of Air Force «Air Force Academy named after professor N.E. Zhukovsky and Y.A. Gagarin» (Voronezh) of the Ministry of Defense of the Russian Federation, (396 064, Voronezh, street Starich Bolshevikov, 54A), email: argentina_85@mail. ru_
On the basis of modeling of circumpolar formations of the Northern hemisphere at the isobaric level 500 hPa for January the multidimensional statistical analysis which allowed to reveal three clusters of objects is carried out. The following characteristics of basic model of circumpolar vortexwere considered: position of the center of masses- area, weight and density of vortex- orientation of dispersion ellipse axis, main standard deviations and coefficient of ellipsecompression. Data of NCEP/DOE AMIP-II reanalysis of atmosphericparameters of 19 792 013 years period were used. The factorial model of circumpolar vortex is constructed, interpretation of factors is given. Objects of the allocated clusters are shown geometrically on the factorial planes. The estimation of statistical difference of factorsdistributions in clusters is given. As a result of the carried-out discriminant analysis the functions discriminating objects of the allocated clusters are constructed, the factors having the greatest impact
on discrimination between groups are defined. _
Keywords: circumpolar vortex of the Northern hemisphere, clusterization, factorial model, discriminant analysis.
В целях обеспечения безопасного и устойчивого развития страны в условиях изменяющегося климата принята Климатическая доктрина Российской Федерации, что указывает на актуальность развития научной основы политики в области климата.
В формировании климатической изменчивости важная роль принадлежит изменчивости циркуляции атмосферы. Глобальная атмосферная циркуляция является одним из механизмов, определяющих все многообразие региональных особенностей климата. Основным звеном атмосферной циркуляции Северного полушария является циркумполярный вихрь (ЦПВ) с центром в районе полюса и западно-восточным движением воздуха. Его характеристики
существенно влияют на формирование и распределение аномальных погодных условий. Поэтому надежное метеорологическое обеспечение оперативной и плановой деятельности потребителей метеорологической информации невозможно без всестороннего изучения внутренней структуры этой компоненты общей циркуляции атмосферы.
Цель работы — совершенствование моделей климатической системы и прогностических схем за счет выявления и учета структурных особенностей ЦПВ Северного полушария.
Для достижения поставленной цели решается задача так называемого разведочного статистического анализа [1], позволяющего выявить вероятностную и геометрическую природу типового метеорологического материала — многомерных данных барической топографии. При этом используются методы факторного, кластерного, дискриминантного анализа.
Исходная информация. В качестве исходного материала использовались данные реана-лиза параметров атмосферы NCEP/DOEAMIP-II [6] по Северному полушарию. Рассматривались значения геопотенциальной высоты Н500(на изобарическом уровне 500 гПа) в узлах регулярной широтно-долготной сетки с шагом 2,5о по срокам 00, 06, 12, 18 час. (Всемирного скоординированного времени) за период 1979—2013 гг. По срочным данным геопотенциальной высоты проведен расчет среднесуточных значений. Задача решалась для центрального зимнего месяца (января), когда циклонический характер ЦПВ наиболее выражен.
Используемая модель ЦПВ. При построении математической модели ЦПВ выбран подход [2−4], определяющий аналогию между проекцией поверхности Н500 на плоскость XOY и двумерным нормальным распределением на этой плоскости случайных величин X, Yс весами pij для точек (xi, y). При этом веса проекций элементов объема «чаши», образованной топографией изобарической поверхностии горизонтальной плоскостью, проходящей через краевую изогипсу Но, определяются по формуле [4]:
p = (Н -Hj COS j
j ЕЕ (Ho -Hj cos j
i j
где Hij — значение высоты изобарической поверхности в узле (x--, y) — j — широта точки (x--, y) — учитываются только точки (узлы), для которых выполняется условие Hij& lt-H0.
В декартовой системе координат XOY центр совмещен с географическим полюсом, ось X направлена на юг вдоль гринвичского меридиана, ось Y — на восток (90о западной — 90о восточной долготы). В этой системе центр циркуляции (xc, yc) определяется как центр масс системы материальных точек.
Характеристиками модели циркумполярного вихря являются: площадь S, ограниченная с юга изогипсой H- условная «масса» М рассматриваемого объема- условная плотность ЦПВг = М/S- геометрические параметры концентрации двумерного распределения около
центра масс: угол 0 между главной осью эллипса рассеяния и осью X- «?, «V — полуоси или главные средние квадратические отклонения эквивалентного единичного эллипса рассеяния- коэффициент сжатия эллипса к = & amp- / «?[2−4].
Используемая модель ЦПВ [2−4] программно реализована для поверхности Н500. В качестве краевой изогипсы, Но выбрана характерная для января изогипса, совпадающая с осью планетарной высотной фронтальной зоны, Но= 540 дам [2].
Кластерный анализ. С целью определения однородных по структуре групп множества циркумполярных образований {Ог}(/ =1, 2, …, п- где п — объем выборки суточных значений характеристик ЦПВ, п=1085) с описательными признаками хц (] = 1, 2, …, р- к расчету приняты параметры модели вихря: Хс, уД М, 2, 0, «?, «, к), проведена процедура иерархической кластеризации объектов многомерного пространства. При этом не требуется априорной информации о распределении генеральной совокупности [1]. Мерой близости, «похожести» различных объектов О/, От выбрано хеммингово расстояние («городских кварталов») [1]:
ён (О/, От) = Ъхи — Хт. (2)
] =1
Предварительно проведена стандартизация параметров по формуле х* = (х-]- - Х])/"Х]-,
чтобы компоненты вектора признаков были однородны по своему физическому смыслу. Для меры (2) влияние отдельных больших разностей (выбросов) уменьшается, так как они не возводятся в квадрат как в евклидовой метрике.
Дерево иерархической кластеризации имеет три отчетливо выделяющиеся ветви. Таким образом, при рассмотрении среднесуточных характеристик циклонических образований можно говорить, в первом приближении, о трех группах (типах) ЦПВ. При этом объекты циркуляции 1-го типа составляют 30,1% случаев, 2-го — 29,6%, 3-го — 40,3%.
В таблице 1 сведены статистические оценки основных характеристик ЦПВ отмеченных кластеров. Из данных таблицы 1 видно, что объекты циркуляции 3-го кластера, в сравнении с объектами кластеров 1, 2, имеют в среднем большую «массу» и «плотность», более сжатый эллипс рассеяния масс. Выявлена особенность распределений угла ориентации эллипса рассеяния по кластерам: для кластера 1 (в основном) 0& gt- 0- для кластера 2 — 0& lt- 0 (т. е. объекты циркуляции кластера 1 имеют направление главной оси эллипса рассеяния восточной ориентации, а кластера 2 — западной) — для кластера 3 угол 0 имеет бимодальное распределение.
Таблица 1
Статистические оценки параметров ЦПВ по кластерам (500 гПа, январь)
Параметр ЦПВ Кластер 1 Кластер 2 Кластер 3
средняя ср. кв. откл. средняя ср. кв. откл. средняя ср. кв. откл.
Хс, км -420,02 394,43 -552,63 330,64 -503,71 284,68
Je, км -91,15 328,29 -59,09 326,82 27,47 264,89
ЛЮ& quot-6, км2 63,39 3,83 65,34 3,32 65,57 2,92
M* 10−9, км2дам 1,37 0,16 1,38 0,17 1,56 0,15
z, дам 21,72 2,54 21,19 2,71 23,86 2,42
0, град. 32,34 14,34 -31,98 17,75 2,80 31,80
Cu, км 2453,85 147,77 2411,18 127,15 2217,08 119,02
Су, км 1617,53 138,59 1688,51 143,82 1788,39 129,93
k 0,66 0,07 0,70 0,07 0,81 0,07
Факторная модель ЦПВ. С целью сжатия информации и определения структуры данных проведен факторный анализ, позволяющий обнаружить скрытые факторы, объясняющие связи между наблюдаемыми признаками объекта. В факторном анализе применяется модель описания переменных Х](] = 1, 2, …, p) линейной комбинацией r ненаблюдаемых факторов fi, fi,… fr, при этом r & lt-p:
Х'-]= ajifi + 0,2/2 + … +ajr- fr + ej, (3)
где Х*] - центрированные и нормализованные исходные данные- fi, i= 1, 2, …, r, называются общими факторами, так как они являются общими для представления всех признаков- ej -остаток или специфичный фактор для переменной Х- aji — неизвестные коэффициенты, имеют специальное название — факторные нагрузки (задача факторного анализа состоит в том, чтобы оценить их некоторым оптимальным образом).
В результате факторного анализа определены общие факторы модели ЦПВ на уровне 500 гПа, объясняющие 89,5% общей дисперсии признаков. Содержательная интерпретация факторов дана в соответствии с факторными нагрузками (после «варимакс-вращения») на параметры ЦПВ, превышающими по модулю 0,7 [1].
Первый фактор наиболее тесно связан с предикторами с, Су, k, что дает основание интерпретировать fi как фактор формы рассеяния условных масс ЦПВ (объясняет 28,9% общей дисперсии признаков).
Фактор f2, судя по нагрузкам на параметры M, z, можно назвать фактором массы ЦПВ (объясняет 24,0% общей дисперсии).
Соответственно другим факторам даны следующие интерпретации: -f — фактор положения центра вихря (15,7% общей дисперсии) — - f — фактор ориентации эллипса рассеяния (11,3% общей дисперсии) — -f — фактор площади вихря (9,6% общей дисперсии).
На рисунках 1−3 объекты выделенных кластеров представлены в разных проекциях: на плоскости факторовf, f (рисунки 1, 2) иf, f (рисунок 3).
3
Фактор 2 2
-3
-3−2-10 1 2 Фактор 1
Рис. 1. Распределение объектов ЦПВ на плоскости факторов /1, / с маркировкой групп, полученных кластеризацией (без объектов кластера 2)
3
Фактор 2 2
-1
-2
-3
-4
2 Фактор31
Рис. 2. Распределение объектов ЦПВ на плоскостифакторов /1, / с маркировкой групп, полученных кластеризацией (без объектов кластера 1) Как показано на рисунках, «сгустки» объектов ЦПВ полученных кластеров выделяются визуально (чтобы не затенять картину отличия скоплений объектов определенных групп, на каждом рисунке не показаны объекты одного из кластеров).
1
0
1
2
1
0
3
2
1
0
1
2
Фактор 4 1. 5
1
0.5 0
и» -ф---. ¦
кластер 1 кластер 2
• V"'-.
¦ •
-0.5 -1 -1.5 -2
-2
2 Фактор 1
Рис. 3. Распределение объектов ЦПВ на плоскостифакторов /1, / с маркировкой групп, полученных кластеризацией (без объектов кластера 3) Оценка однородности кластеров. В целях определения однородности групп, проведена проверка гипотезы о принадлежности факторов рассматриваемых кластеров одним и тем же распределениям. Использовались непараметрические методы (или свободные от распределения): на основе критериев Колмогорова — Смирнова и рангового Манна — Уитни[6], — не требующие предпосылки существования закона распределения.
В таблице 2 представлены результаты оценки по совокупности двух указанных критериев статистического различия распределений соответствующих факторов в рассматриваемых группах на уровне значимости, а = 0,05 (знаком «+» отмечены случаи принятия статистической гипотезы, что две совокупности имеют одно и то же распределение- знаком «-» -принятия альтернативной гипотезы).
Таблица 2
Оценки статистического различия распределений факторов в кластерах (а= 0,05)
3
1
0
1
Кластеры /1 /2 /3 /4 /5
1, 2 — + + - -
1, 3 — - - - -
2, 3 — - + - +
Как следует из данных таблицы 2, существенно различаются распределения всех факторов в кластерах 1, 3- в кластерах 1, 2 различаются только распределения факторов /1, /4, /5- в кластерах 2 и 3 — /1, /2, /4.
На рисунке 4 даны в сравнении статистические характеристики факторов по группам (мода, 50-ти процентный интервал и размах варьирования), графическое представление которых визуально подтверждает данное утверждение.
Рис. 4. Сравнение для выделенных кластеров статистических характеристик факторов: а) /- б) /2- в) /3- г) /4- д) /5 Дискриминантный анализ. Для выяснения, какие факторы вносят свой вклад в дискриминацию между полученными совокупностями объектов ЦПВ, в задаче использовался дискриминантный анализ. В результате дискриминантного анализа, во-первых, построены классификационные функции, предназначенные для определения того, к какой группе наиболее вероятно может быть отнесен отдельный объект. Во-вторых, определены дискриминантные функции, разделяющие возникающие совокупности.
Каждой группе (совокупности объектов) соответствует своя классификационная функция йс-(1 = 1, 2, 3), определяемая вектором коэффициентов и константой:
ёс1 =-2,734- 1,663/1- 0,808/2 + 0,396/3+1,779/4- 0,692/5, (4)
ёс2 = -2,726- 0,999/1- 0,986/2+ 0,056/3- 2,212/4+ 0,178 /5, (5)
ёс3 = -2,079 + 1,979/1+1,329/2- 0,337/3 + 0,294/4+ 0,387/5, (6)
Данные функции классификации (14)-(16) можно использовать для прямого вычисления показателя классификации для некоторых новых значений. Наблюдение (объект циркуляции) считается принадлежащим той совокупности, для которой получен наивысший показатель классификации.
В результате проведенного анализа определены дискриминантные функции ё1 и Функция ё1 в основном дискриминирует между группой 3 и объединением групп 1 и 2. Функция й?2 дискриминирует между группой 1 и 2. По коэффициентам полученных дискри-минантных функций, представленных в стандартизованном виде:
ё1 = - 1,004/1- 0,781/2+ 0,226/э- 0,070/4- 0,267/5, ё2 = - 0,095/1+ 0,125 /2+ 0,101/э + 1,001/4- 0,281/5,
(8)
можно судить о степени влияния каждого предиктора на дискриминацию групп. Так, согласно (7), (8), наибольшее влияние на дискриминацию между группой 3 и совокупностью групп 1, 2 оказывают факторы /1 и /2 (формы и массы ЦПВ) — на дискриминацию между группами 1 и 2 наибольшее влияние оказывает фактор /4 (ориентации эллипса рассеяния условных масс ЦПВ), что подтверждает и распределение объектов на рисунках 1−3.
Другой целью применения дискриминантного анализа является проведение классификации. Если модель установлена и получены дискриминирующие функции, возникает вопрос о том, как хорошо они могут предсказывать, к какой совокупности принадлежит конкретный объект циркуляции. Такой анализ не проводился, поскольку это цель следующего этапа -модельного дискриминантного анализа. Цель разведочного анализа — предварительная оценка внутренней структуры исходной совокупности.
Достоверность полученных результатов подтверждается большим объемом данных, привлеченных для исследования, и применением классических методов статистической обработки информации.
Полученные результаты позволяют говорить о дифференцированном использовании данных о состоянии ЦПВ в разрабатываемых погодно-климатических моделях, тем самым повысить эффективность метеорологического обеспечения.
1. Айвазян С. А., Бухштабер В. М., Енюков И. С., Мешалкин Л. Д. Прикладная статистика: Классификация и снижение размерности. Справочное издание, [под ред. С. А. Айвазяна]. -М.: Финансы и статистика, 1989. — 607 с.
2. Багров Н. А., Орлова И. И. К вопросу определения центра циркуляции атмосферы // Труды Гидрометцентра СССР. — Вып. 211. — Л.: Гидрометеоиздат, 1978. — С. 3−14.
3. Багров Н. А. О центре циркуляции // Метеорология и гидрология. — 1975. — № 2. — С. 3−11.
4. Глызь Г. А. О некоторых характеристиках циркумполярного вихря // Труды ВНИИГМИ-МЦД. — Вып. 58. — М.: Гидрометеоиздат, 1979. — С. 98−104.
5. Справочник по прикладной статистике. В 2 т. Т. 2: [под ред. Э. Ллойда, У. Ледермана- пер. с англ. под ред. С. А. Айвазяна, Ю.Н. Тюрина]. — М.: Финансы и статистика, 1990. — 526 с.
6. Капаш^и М. еЫ. КСЕР/БОЕАМХР-ПКеапа^в // Ви11. Ашепсап Мйеог. Бос., 2002. — Уо1. 83. — Р. 1631−1643.
Список литературы
Рецензенты:
Куролап С. А., д.г.н., профессор, заведующий кафедрой геоэкологии и мониторинга окружающей среды, ФГБОУ ВПО «Воронежский государственный университет», г. Воронеж- Кузнецов И. Е., д.т.н., начальник13 кафедры гидрометеорологического обеспечения, ФГБОУ ВПО «Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н. Е. Жуковского и Ю.А. Гагарина» Министерства обороны Российской Федерации, г. Воронеж.

ПоказатьСвернуть
Заполнить форму текущей работой