CC BY
67
УДК 533.933; 533.932
В. М. Журавлев, С. В. Летуновский
АНАЛИЗ ДОЛГОВРЕМЕННОЙ ЭВОЛЮЦИИ АКТИВНОСТИ СОЛНЦА НА ОСНОВЕ РЯДА ЧИСЕЛ ВОЛЬФА (I. Методика)1
Аннотация. Представлены результаты исследования эволюции статистических характеристик ряда чисел Вольфа и ряда групп пятен на масштабах их изменчивости порядка 100 лет. Строится полуэмпирическая модель вероятностного распределения чисел Вольфа. Излагается метод моментов в применении к задаче вычисления эволюции параметров распределения чисел Вольфа.
Ключевые слова: солнечная активность, числа Вольфа, статистика, эволюция распределения, метод моментов.
Abstract. Statistical characteristics of Volf number and number of groups sequence distributions and their evolution in time on scales about 100 years are investigated.
The statistical semi-empirical model this distributions are represented. The statistical moment’s method in application to a problem of evaluation Volf numbers distribution parameters is discussed.
Keywords: solar activity, Volf numbers, statistics, distribution evolution, moments method.
Введение
Задача исследования солнечной активности, кроме важности в выяснении физических механизмов, управляющих динамикой процессов, происходящих на Солнце, имеет большое значение для выяснения и прогноза влияния этой изменчивости на изменения климата на Земле. Такие исследования особенно важны, например, в связи с наблюдаемым сдвигом средней температуры на Земле в сторону ее увеличения, что не находит пока однозначного объяснения. Решение общей задачи предсказания солнечной активности осуществляется в настоящее время множеством различных способов, ориентированных на выделение и анализ различных периодов ее изменчивости. Поскольку одной из наиболее ярко выраженных квазипериодических составляющих изменчивости солнечной активности является 11-летний цикл, то имеется множество работ, связанных с решением задачи предсказания очередного солнечного квази-одиннадцатилетнего цикла на основе информации
о предыдущих. Но в изменчивости характеристик самого 11-летнего цикла обнаруживаются как более короткие периодичности (например, 2-летняя составляющая), так и длинно-периодические изменения, которым посвящено также множество работ [1-22] (см. также библиографию там). Цель большинства работ состоит в отыскании прогностических параметров, исследуя эволюцию которых можно с той или иной степенью надежности предсказать изменение солнечной активности в очередном цикле. Большинство проведенных исследований опирается на данные в виде чисел Вольфа и использует различные типы параметров, которые могут быть оценены из самого ряда.
1 Работа выполнена при частичной поддержке РФФИ, проект 08-01-97013-р_поволжье_а.
Это и авторегрессионные модели [11], различные типы индексов, например: индекс, связанный с числом групп пятен [13]; параметр экспоненты Харста [16] и т.д. В значительно меньшей степени в современных работах уделяется внимание исследованию изменчивости непосредственно параметров статистического распределения чисел Вольфа. В работе [23] приведены гистограммы распределения ряда чисел Вольфа, восстановленного по палеоданным для больших периодов времени. Однако на таких больших масштабах времени, как тысячи лет, значения индексов солнечной активности являются малонадежными и подвержены значительному статистическому разбросу. Поэтому такую работу имело бы смысл проводить для данных, которые установлены с достаточной надежностью. В известной монографии [10] исследован целый набор статистических свойств различных характеристик параметров солнечных пятен на Солнце, однако вероятностному распределению чисел Вольфа практически внимания не уделялось. Вместе с тем можно ожидать, что исследование самого вероятностного распределения ряда чисел Вольфа может дать полезную информацию об эволюции солнечной активности. Эти надежды основываются на том, что вероятностные распределения содержат в себе всю детерминированную информацию о случайном процессе.
Вместе с тем при исследовании вероятностного распределения такого случайного процесса, как ряд чисел Вольфа, необходимо отдавать себе отчет, что сам этот процесс не стационарен и параметры распределения могут зависеть от времени. Именно зависимость от времени параметров распределения может дать полезную информацию об изменениях в состоянии термодинамической машины Солнца. Однако для этого необходимо иметь возможность связывать параметры распределения, например, его моменты, с какими-либо параметрами физической модели, описывающей такую изменчивость.
В настоящей работе проводится анализ ряда чисел Вольфа с помощью относительно простой полуэмпирической модели, основанной на гипотезе о существовании двух несовместных механизмов возникновения солнечных пятен. Один из них соответствует некоторому равновесному процессу, а второй интерпретируется как «взрывной». Обоснование такого представления приводится в данной статье. Исследования предпринимаются с целью показать возможность извлекать из такой модели существенную информацию о характеристиках солнечной термодинамической машины как целого. Это позволяет также выявить ряд особенностей в эволюции солнечной активности, которые связаны с глобальными изменениями в функционировании этой машины на больших масштабах времени, значительно превышающей характерный временной масштаб солнечного цикла. В отличие от плавных квази-периодических изменений самих чисел Вольфа, обнаруживаемые изменения статистических характеристик вероятностного распределения имеют отрезки времени, когда параметры меняются очень быстро - практически скачкообразно. Быстро изменяющиеся процессы на Солнце обычно связывают со вспышками, характерный временной масштаб которых оценивается, как правило, минутами или десятками минут. Изменения в статистических характеристиках среднемесячных чисел Вольфа, как кажется, не могут быть связаны со вспышечной активностью. Однако можно предположить, что вспышки типа вспышки 1859 г. могут быть связаны с существенными изменениями термодинамической системы Солнца, поскольку энергия, выбрасываемая в таких
вспышках, уже может быть сравнимой с энергией глобальных процессов на Солнце. В настоящей работе проводится исследование параметров предложенной модели в течение последних 250 лет. Основным набором данных, который используется для анализа изменчивости активности Солнца, выбран ряд ежемесячных чисел Вольфа с 1749 по 2009 г., имеющийся на Интернет-сайте [1]. Вспомогательным набором данных, который позволяет увеличить интервал, на котором можно получить оценки исследуемых параметров, является ряд ежемесячных чисел групп пятен. Этот ряд имеется на том же Интернет-сайте, но, по всей видимости, менее надежен. Однако в силу того, что в работе анализируются не отдельные значения рядов, а их усредненные за 100 лет статистические характеристики, флуктуации оценок оказываются не очень существенными. Задачей этого исследования является выявление скрытых механизмов, управляющих динамикой солнечной активности на больших масштабах времени, и отыскание новых прогностических параметров, позволивших бы сделать предсказания о характере возможных изменений в функционировании солнечной термодинамической системы.
1. Скользящие ряды чисел Вольфа
Для выявления изменчивости статистических характеристик ряда чисел Вольфа необходимо иметь множество рядов, сдвинутых во времени, для каждого из которых строится гистограмма распределения. Для этого в данной работе используется метод скользящих рядов. Метод скользящих рядов сводится к сравнительному исследованию свойств рядов чисел Вольфа
Wa ={wa г- }/=!, полученных из полного ряда W = {^ }2=Г749 1 общей длиной
N =3115 с помощью процедуры выделения части ряда заданной длины N, начиная с некоторого элемента ряда с номером a :
wa,i = Was+(г -1), г' = 1,2,---, N.
В этом соотношении целое число a определяет номер каждого ряда в создаваемом многомерном наборе данных: a = 0,1,2,...,M . Фактически номер a связан со временем ta = aДt между начальными значениями рядов, которое можно рассматривать в качестве реального параметра, позволяющего упорядочить во времени изменения средних характеристик отдельных рядов Wa . Целое число 5 определяет величину временного сдвига между началами отдельных рядов. Числа М, 5 и N связаны между собой и общей длиной ряда N0 соотношением 5 = [(N0 - Щ/М], [к] - целая часть числа к . Обозначим через W N (а) и а N (а) среднее значение ряда с номером сдвига а длиной N.
2. Распределение вероятностей среднемесячных чисел Вольфа
Гистограмма распределения вероятности появления определенного числа пятен (рассчитанных с помощью формулы Вольфа) в течение месяца, вычисленная по всему ряду (рис. 1), достаточно хорошо описывается показательным распределением следующего вида:
Р0(и) = -ке-Хп . (1)
1750 1800 1850 1900 1950 2000
І
а)
0 50 100 150 200 250
W
б)
Рис. 1. Числа Вольфа: а - ежемесячный ряд; б - гистограмма
Это устанавливается с помощью критерия у2 и, собственно, наглядно подтверждается близостью среднего значения полного набора ежемесячных чисел Вольфа W к стандартному отклонению этого же полного набора а: W = 51,9, а = 44,3. Для теоретического распределения (1) эти величины
должны совпадать точно: wo(t) = а0 = ^-1. Некоторое различие в величинах W, а указывает на существование дополнительной составляющей в истинном распределении, которое несколько искажает распределение (1) . Длина иссле-
дуемого ряда в целом и отдельных скользящих рядов вполне достаточна для того, чтобы выявить искажения в структуре распределения со временем и установить форму дополнительной составляющей с требуемой для прогноза степенью надежности.
Для учета этого искажения в данной работе предлагается модель в форме смеси двух распределений, одним из которых является (1), а вторым -некоторое распределение, которое можно выбрать, изучая результаты моделирования.
3. Особенности эволюции моментов распределения чисел Вольфа
Расчет параметров распределения будет проводиться в данной работе на основе метода моментов. Прежде чем переходить к анализу эволюции распределений чисел Вольфа, рассмотрим особенности в долговременной эволюции самих моментов распределений, на основе которых и строятся оценки параметров распределения р(п,Г). На рис. 2 и 3 представлены графики изменения оценок моментов распределения р(п, t) скользящих рядов. Наиболее заметными изменениями со временем оценок математического ожидания распределения р(п, t) (рис. 2) является монотонный рост этой величины. Если определять момент времени, к которому относится оценка среднего по началу скользящего ряда, то рост среднего значения начинается с середины XIX в. и остается таковым и для последнего скользящего ряда с началом в 1909 г.
г
Рис. 2. Эволюция среднего значения скользящих рядов чисел Вольфа.
По оси абсцисс обозначен номер месяца, начиная с января 1749 г.
Аналогичный, но с некоторого момента более крутой рост наблюдается и во втором моменте (рис. 3,а). Основной вывод, который можно сделать, исходя из наблюдаемого одновременного роста среднего значения чисел Вольфа и второго момента, состоит в том, что со временем возрастает не только вероятность наблюдения все большего числа пятен, но и вероятность флуктуаций этой величины от среднего.
а)
б)
ю
в) г)
Рис. 3. Эволюция центральных моментов скользящих рядов чисел Вольфа: а - т = 2; б - т = 3; в - т = 4; г - т = 5
№ 4 (16), 2010 Физико-математические науки. Физика
Поскольку основной составляющей в распределении р(и, г) является показательное распределение, то это означает, что величина стандартного отклонения, характеризующая среднестатистическое отклонение от среднего, почти равна самому среднему. В свою очередь это означает, что с достаточно большой вероятностью в ежемесячном исчислении могут наблюдаться не только периоды большого числа пятен, но и периоды с числом пятен, практически равным нулю. Это формально объясняет со статистической точки зрения наблюдаемое в начале текущего цикла солнечной активности аномально низкое число пятен, которое в феврале 2010 г. сменилось резким увеличением числа пятен и вспышек. Такое поведение характерно для процесса с показательным распределением.
Эволюция старших моментов представлена на рис. 3. Наиболее существенным элементом в эволюции старших моментов является наличие скачка в XIX в. Этот скачек прослеживается на всех моментах порядка к > 1.
Поскольку моменты оценивались на интервале времени 100 лет, то абсолютное положение скачка на графиках моментов (рис. 3) установить можно лишь, связав какой-либо известный факт из истории наблюдения солнечной активности с характеристиками эволюции исследуемых параметров. В качестве такого временного маркера можно избрать известную мощную вспышку 1859 г. Эта вспышка хорошо прослеживается по данным химических маркеров (нитратный след [24, 25]) во льду Антарктики и Гренландии. Согласно графикам на рис. 3 наиболее значительный скачок в моментах распределения чисел Вольфа наблюдается спустя почти сто лет от начала ряда, датируемого 1749 г. Скачок в моментах распределений как раз приходится на время вблизи вспышки 1859 г., если в качестве момента, к которому привязываются значения моментов, полученные на отрезках длиной 100 лет, выбрать начало этого отрезка. Сам факт того, что в параметрах распределения изменения происходят, начиная с момента вспышки, указывает на то, что само распределение формируется в результате происходящих изменений на Солнце и является их следствием и индикатором. В результате мы получаем в руки новый тип индексов, характеризующих скачкообразные изменения в состоянии Солнца. Однако сами моменты распределения чисел Вольфа с физической точки зрения мало информативны. Для получения более информативных параметров далее предлагается модель вероятностного распределения чисел Вольфа, на основании которой можно уже строить гипотезы именно о физических механизмах, управляющих вспышками.
4. Модель распределения чисел Вольфа
Модель распределения, используемая далее, имеет следующий вид:
р„ (и) = (1 - р(Г)) ^0 (г)е(г)п + р(г)р1 (и, Г). (2)
Это распределение можно интерпретировать как полную вероятность появления определенного числа пятен под действием двух несовместных механизмов. Вероятность срабатывания одного из этих механизмов (основного) равна 1 - р(Г), а второго р(Г). Плотность вероятности появления ^ чисел Вольфа при срабатывании первого механизма описывается показательным
^ (Г) п
распределением р0 (м>, Г) = А0 (Г)е 0 , а второго - р1 (и, Г), функциональ-
ный вид которого еще следует выбрать на основе анализа эмпирического распределения ряда чисел Вольфа.
Показательное распределение рд(м>, Г) можно интерпретировать как распределение по энергиям, если при действии этого механизма каждое появившееся пятно имеет некоторое слабо меняющееся среднее значение энергии Е. В этом случае показательное распределение есть распределение Больцмана по энергиям. Формально мы можем при такой интерпретации считать, что ^0 (Г) = Е0 (Г)/кТ* (Г), где Т* (Г) - некоторая условная температура системы. Эта интерпретация не является достаточно обоснованной и приводится здесь лишь качестве иллюстрации возможных типов объяснения этой части распределения.
Как видно из рис. 1,б, наиболее существенные отклонения от показательного распределения наблюдаются в области чисел Вольфа порядка 50-80 и практически исчезают при больших и малых (вблизи нуля) значениях ^ . Это означает, что второй механизм срабатывает реже, но порождает сразу множество пятен, число которых находится вблизи значения 50-80, которое должно соответствовать модели распределения р1 (и, Г). В силу этого мы можем предполагать, что такой механизм носит «взрывной» характер, а само распределение р1(и, Г) можно аппроксимировать либо распределением типа максвелловского:
с параметрами д, ^ (Г) и а(Г), которые должны вычисляться из экспериментальных данных. К этим параметрам добавляется еще и параметр р(Г), который наиболее важен с физической точки зрения, поскольку определяет вероятность срабатывания взрывного механизма в солнечной активности. Для оценивания этих параметров воспользуемся методом моментов [26, 27].
5. Метод моментов
Метод моментов базируется на вычислении параметров теоретических распределений на основе оценок моментов случайной величины по эмпирическим данным. Распределения р0(и, Г) и р(1,11)(и, Г; д) содержат по одному параметру (д считается заданным). Еще один параметр р(Г) принадлежит смеси распределений р(и, Г). Таким образом, для построения распределения р(и, Г) требуется вычислить три параметра. Для этого согласно методу моментов достаточно вычислять три момента эмпирического распределения.
Предполагая, что р1(и,Г) нормировано, получаем следующие соотношения для параметров полного распределения:
«2 (і) = °2(() = (1 - Р(і Жа0 (і) + (м(і) - мо(і ))2 ] + р(і )[СТ! (і) + (м(і) - м\(і ))2 ]; (6)
(3)
либо распределением следующего вида:
р(П)( п, і) = 1 )п,
(4)
т (і) = м(і) = (1 - р(і))мо + р(і )м>\;
(5)
m3(t) = (1 -p(t))[2^03(t) + (w(t)-wo(t))3] + p(t)[A(t) + (w(t)-wi(t))3], (7) где
wi(t) = J«Pi(«, t)dn, a2(t) =
o
= J(n - wi (t))2 P1 (n, t)dn, A1 (t) = J(n - wi (t))3 P1 (n, t)dn .
0 0
В этих уравнениях в силу того, что распределение Po(n, t) имеет вид (1) , wo (t) и ao (t) выражаются через одну функцию ^ (t):
wo(t ) = ao (t ) = V(t).
Эти уравнения содержат пять неизвестных функций: p(t), ^o(t),
wi(t), ai(t), Ai(t). Для случая выбора распределения Pi(n, t) в форме (3) и (4), параметры wi (t), a (t), Ai (t) связаны между собой так, что только один из них оказывается произвольным. Для случая в качестве свободного параметра удобно выбирать ^ (t) = n-1, а в случае (4) - параметр no (t).
Поскольку в этом случае число параметров оказывается равным 3, то система (5)-(7) из трех уравнений позволяет рассчитать все параметры двух распределений и вероятность взрыва p(t). Вычисления проводились для каждого ряда из набора скользящих рядов с шагом в один месяц. При этом численно решалась система из трех алгебраических уравнений суммарного порядка, равного 7. Из корней этого уравнения отбирались корни, удовлетворяющие определенным требованиям.
Заключение
Разработан метод анализа эволюции параметров вероятностнорго распределения чисел Вольфа с целю выяснения долговременной эволюции солнечной активности. Предлагаемый подход строится на методе моментов [26, 27] в сочетании со специальным выбором смеси распределений, которые наилучшим образом характеризуют форму распределения чисел Вольфа. Важным элементом такого подхода является интерпретация смеси распределений как полной вероятности наблюдения пятен. Это позволяет характеризовать отдельные элементы распределения как вероятности несовместных механизмов образования пятен. Одно из распределений, как вероятность появления пятен в «равновесии», а второе, как неравновесную, или взрывную, составляющую. Метод предназначен для проведения исследований ряда чисел Вольфа, результаты которых изложены в следующей работе.
Список литературы
1. [Электронный ресурс] FTP сервер NASA. - URL: ftp://itp.ngdc.noaa.gov/STP/ SOLAR_DATA/SUNSPOT_ NUMBERS/
2. [Электронный ресурс] Basu S., Antia H. M. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] oooi294v1
3. [Электронный ресурс] Petrovay K. - URL: http://arxiv.org;articl-id: //arXiv:[astro-ph] 0010096v2
4. [Электронный ресурс] Stefano Sello. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0010106v1
5. Ziеeba, S. Cover illustration: First Doppler image of a solar-type G dwarf / S. Ziеeba, J. Maslowski, A. Michalec, A. Kulak // Astronomy & Astrophysics. -2001. - V. 377. - № 1. - Р. 297-311.
6. [Электронный ресурс] Dean-Yi Chou, Alexander Serebryanskiy. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0405175v1
7. Atac, T. Flare Index During the Rising Phase of Solar Cycle 23 / T. Atac, A. Ozguc // Solar Physics. - 2004. - V. 198. - № 2. - Р. 399-407.
8. [Электронный ресурс] R. Cameron and M. SchEussler. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0612693v1
9. [Электронный ресурс] Arnab Rai Choudhuri, Piyali Chatterjee, Jie Jiang. - URL: http://arxiv.org;articl-id: arXiv:astro-ph/0701527v1
10. Витинский, Ю. И. Статистика пятнообразовательной деятельности Солнца / Ю. И. Витинский, М. Копецкий, Г. В. Куклин. - М. : Наука, 1986. - 294 с.
11. [Электронный ресурс] K. M. Hiremath. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph]0704.1346v1
12. [Электронный ресурс] J. Bushby, Steven M. Tobias. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0704.2345v1
13. [Электронный ресурс] A. G. Tlatov. - URL: http://arxiv.org;articl-id: //arXiv:[astro-ph]0706.1624;arXiv:[astro-ph] 0703681
14. [Электронный ресурс] M. Scffiussler. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0712.1917v1
15. [Электронный ресурс] R. Cameron, M. Scffiussler. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0806.2833v1
16. [Электронный ресурс] A. Kilcik1, C.N.K. Anderson, J.P. Rozelot, A. Ozguc. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0811.1708v5
17. [Электронный ресурс] P.A.Semi. - URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0903.5009
18. [Электронный ресурс] D. Salabert, R. A. Garcm, P. L. Pallre1,S. J. Jimrenez-Reyes. -URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0907.3888v1
19. [Электронный ресурс] Leonid V. Didkovsky, Darrell L. Judge, Seth R. Wieman. -URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0911.0870v1
20. Joshi, B. Periodicities in sunspot activity during solar cycle 23 / Bhuwan Joshi, P. Pant1, P. K. Manoharan // Astronomy & Astrophysics. - 2006. - V. 452. - № 2. -P. 647-650.
21. [Электронный ресурс] L. Rogers, Mercedes T. Richards, Donald St. P. Richards. -URL: http://arxiv.org;articl-id: arXiv:[astro-ph] 0606426v3
22. Shea, M. A. History of Energetic Solar Protons for the Past Three Solar Cycles Including Cycle 22 Update / M. A. Shea, D. F. Smart // Solar Physics. - 1990. - June. -V. 127. - P. 297-320.
23. [Электронный ресурс] I. G. Usoskin. A History of Solar Activity over Millennia. -URL: http://www.livingreviews.org/lrsp-2008-3
24. Кочаров, Г. Е. Естественные архивы солнечной активности и термоядерной истории Солнца за последние миллионы лет / Г. Е. Кочаров // Соросовский образовательный журнал. Физика. - 2000. - № 1. - С. 91-95.
25. Gisela, A. M. Dreschhoff, Edward J. Zeller Evidence Of Individual Solar Proton Events in Antarctic Snow / A. M. Gisela // Solar Physics. - 1990. - V. 127. -P.333-346.
26. Митропольский, А. К. Техника статистических вычислений / А. К. Митропольский. - М. : Наука, 1971. - 576 с.
27. Крамер, Г. Математические методы статистики / Г. Крамер. - М. : Мир, 1975. -625 с.
Журавлев Виктор Михайлович
доктор физико-математических наук, профессор, кафедра теоретической физики, Мордовский государственный педагогический институт им. М. Е. Евсевьева
E-mail: zhvictorm@mail.ru
Летуновский Сергей Владимирович старший преподаватель, кафедра физики, Ульяновский государственный университет (филиал в г. Димитровграде)
E-mail: grayser@bk.ru
Zhuravlev Viktor Mikhaylovich Doctor of physical and mathematical sciences, professor, sub-department of theoretical physics, Mordovia State Pedagogical University named after M. E. Evsevyev
Letunovsky Sergey Vladimirovich Senior lecturer, sub-department of physics, Ulyanovsk State University (affiliated branch in Dimitrovgrad)
УДК 533.933; 533.932 Журавлев, В. М.
Анализ долговременной эволюции активности солнца на основе ряда чисел Вольфа (I. Методика) / В. М. Журавлев, С. В. Летуновский // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2010. - № 4 (16). - С. 120-130.
