Временные изменения глобальных характеристик вращения Солнца
З використанням методів спектрального аналізу отримано глобальні характеристики обертання спостережуваних у лінії НеI 1083 λ нм сонячних структур для трьох видів аналізу: на всьому часовому інтервалі спостережень тривалістю близько 26 років; на п’ятьох часових інтервалах, що відповідають епохам ма...
Saved in:
| Published in: | Проблемы управления и информатики |
|---|---|
| Date: | 2007 |
| Main Authors: | , , |
| Format: | Article |
| Language: | Russian |
| Published: |
Інститут кібернетики ім. В.М. Глушкова НАН України
2007
|
| Subjects: | |
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/206965 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Временные изменения глобальных характеристик вращения Солнца / Я.И. Зелык, Н.Н. Степанян, О.А. Андреева // Проблемы управления и информатики. — 2007. — № 1. — С. 127-135. — Бібліогр.: 2 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859966436286398464 |
|---|---|
| author | Зелык, Я.И. Степанян, Н.Н. Андреева, О.А. |
| author_facet | Зелык, Я.И. Степанян, Н.Н. Андреева, О.А. |
| citation_txt | Временные изменения глобальных характеристик вращения Солнца / Я.И. Зелык, Н.Н. Степанян, О.А. Андреева // Проблемы управления и информатики. — 2007. — № 1. — С. 127-135. — Бібліогр.: 2 назв. — рос. |
| collection | DSpace DC |
| container_title | Проблемы управления и информатики |
| description | З використанням методів спектрального аналізу отримано глобальні характеристики обертання спостережуваних у лінії НеI 1083 λ нм сонячних структур для трьох видів аналізу: на всьому часовому інтервалі спостережень тривалістю близько 26 років; на п’ятьох часових інтервалах, що відповідають епохам максимумів і мінімумів сонячної активності; на ковзному річному часовому інтервалі зі зсувом на півроку.
By using the spectral analysis methods the global characteristics of the solar structures rotation, observed in line НеI 1083 λ nm are obtained for three types of analysis: on the whole observation time interval of duration about 26 years; on five time intervals corresponding to maxima and minima of solar activity; epoch on the sliding one year time interval with halfyear shift .
|
| first_indexed | 2025-12-07T16:21:23Z |
| format | Article |
| fulltext |
© Я.И. ЗЕЛЫК, Н.Н. СТЕПАНЯН, О.А. АНДРЕЕВА, 2007
Проблемы управления и информатики, 2007, № 1 127
УДК 523.9-327:519.246.8:517.443:519.71
Я.И. Зелык, Н.Н. Степанян, О.А. Андреева
ВРЕМЕННЫЕ ИЗМЕНЕНИЯ ГЛОБАЛЬНЫХ
ХАРАКТЕРИСТИК ВРАЩЕНИЯ СОЛНЦА
В настоящей работе рассматривается вращение Солнца и его характеристики
в верхней хромосфере. Наблюдательный материал в виде синоптических карт
в линии HeI λ 1083 нм за три последних солнечных цикла (с начала 1977 г. до
конца сентября 2003 г.) был любезно предоставлен нам доктором Харви из Китт
Пик обсерватории США. Обработка наблюдений и методы анализа описаны
в [1, 2]. Синоптические карты преобразованы во временные ряды для нескольких
типов солнечных образований в соответствии с характерными для них интервала-
ми эквивалентных ширин W. На рис. 1 приведены изменения со временем эквива-
лентных ширин N внутри пятиградусных широтных зон для активных областей
AR (рис. 1, а) и для корональных дыр CH (рис. 1, б). На рис. 1, а видны изменения
широтного распределения W в течение трех 11-летних солнечных циклов (21-, 22-
и 23-го): появление активности в начале цикла на высоких широтах и постепенное
понижение широт возникновения активных областей к концу цикла.
1650 1700 1750 1800 1850 1900 1950 2000
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
1977 1981 1984 1988 1992 1996 1999 2003
N 75
N 65
N 55
N 45
N 35
N 25
N 15
N 05
S 10
S 20
S 30
S 40
S 50
S 60
S 70
S 80
W
AR
Время, каррингтоновские обороты
1650 1700 1750 1800 1850 1900 1950 2000
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
0
80
1977 1981 1984 1988 1992 1996 1999 2003
N 75
N 65
N 55
N 45
N 35
N 25
N 15
N 05
S 10
S 20
S 30
S 40
S 50
S 60
S 70
S 80
W
CH
Время, каррингтоновские обороты
а б
Рис. 1
С применением разработанных авторами алгоритмов цифровой обработки
наблюдений и спектрального анализа (Periodogram, Welch, Multitaper) получены
спектральные плотности, значения значимых пиков в них и соответствующих им
периодов в интервале периодов 20÷35 сут для всех типов солнечных образований.
Применительно к рядам, полученным для структур, наблюдаемых в линии HeI λ
1083 нм, проведено три типа спектрального анализа: 1) по всему временнóму интер-
валу продолжительностью около 26 лет; 2) по пяти временным интервалам, со-
128 ISSN 0572-2691
ответствующим эпохам максимумов и минимумов солнечной активности; 3) сколь-
зящий спектральный анализ во временнóм окне длиной 1 год со сдвигом на полгода.
Среднее вращение солнечных структур на всем интервале наблюдений.
Остановимся на результатах, полученных в первом случае при построении спек-
тров мощности по всем 26-летним временным рядам. На рис. 2 приведены спек-
тры мощности в интервале периодов 20÷35 сут для активных областей во всех
широтных зонах N-полусферы.
0
200
400
0
50
100
0
20
40
0
20
40
0
50
100
0
10
20
0
20
40
20 22 24 26 28 30 32 34 36
0
20
40
N80
N75
N70
N65
N60
N55
N50
N45
P
S
D
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
0
1000
2000
0
2000
4000
0
2000
4000
0
1000
2000
0
1000
2000
0
200
400
20 22 24 26 28 30 32 34 36
N40
N35
N30
N25
N20
N15
N10
N5
Период (дни)
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
0
2000
4000
20 22 24 26 28 30 32 34 36
0
2000
4000
20 22 24 26 28 30 32 34 36
20 22 24 26 28 30 32 34 36
Рис. 2
Проблемы управления и информатики, 2007, № 1 129
Изменение основных периодов с широтой, встречающиеся в спектрах мощ-
ности удобно рассмотреть на трехмерных графиках или в изолиниях.
На рис. 3 показана зависимость спектральной плотности мощности PSD(φ, Р)
от широты φ и периода Р, построенная для активных областей. На рис. 4 приведе-
на зависимость PSD(φ, Р) в виде изолиний для активных областей и корональных
дыр. Значение широты и величина периода мест с наибóльшим значением PSD на
рис. 4 даны в табл. 1.
20
25
AR
35
0
– 80 – 60
– 40
– 20 0 20 40
60 80
500
1000
1500
2000
2500
3000
3500
4000
P
S
D
Период (дни)
30
Широта (в градусах)
Рис. 3
0 – 80 – 60 – 40 – 20 20 40 60 80
22
24
26
28
30
32
AR
П
ер
и
о
д
(
д
н
и
)
0 – 80 – 60 – 40 – 20 20 40 60 80
22
24
26
28
30
32
CH
П
ер
и
о
д
(
д
н
и
)
а Широта (в градусах) б
Рис. 4
Таблица 1
Активные области Корональные дыры
Наиболее мощные
периоды (сут)
Широта (градусы)
Наиболее мощные
периоды (сут)
Широта (градусы)
25,68 – 25 30,32 – 24
28,18 – 17,5 29,25 30
27,43 – 16,2 28,30 64,8
25,4 – 15,6
26,81 20
28,92 20
Самый мощный период вращения во всех полученных спектрах мощности —
26,81 сут на широте +20°. Сравнимый с ним по мощности период 28,18 сут на
широте –17,5°. Оба они относятся к активным областям.
Вращение в N- и S- полусферах существенно различается не только по вели-
чине максимального периода и его широте. Различаются и средние спектры мощ-
ности для двух полусфер. Это можно увидеть на рис. 5 для активных областей (а)
и корональных дыр (б).
130 ISSN 0572-2691
28 20 22 24 26 30 32 34 36
0
200
400
600
800
С
р
ед
н
ее
з
н
ач
ен
и
е
P
S
D
Период (дни)
18
1
2
1 — AR N-полусфера
2 — AR S-полусфера
а
28 20 22 24 26 30 32 34 36
0,0
0,5
1,0
С
р
ед
н
ее
з
н
ач
ен
и
е
P
S
D
Период (дни)
18
1
2
1 — CH N-полусфера
2 — CH S-полусфера
1,5
2,0
2,5
б
Рис. 5
Вращение в эпохи максимумов и минимумов солнечной активности. Для
эпох трех максимумов и двух минимумов солнечной активности были построены
спектры мощности для активных областей и корональных дыр. Из всего материа-
ла выделялись массивы данных, относящиеся к пяти временным интервалам,
представленным в табл. 2.
Таблица 2
№ п /п
Интервал времени, годы
(дни от начала массива данных)
Характеристика активности
1 1978,5÷1981 (600÷1700) Максимум 21-го цикла
2 1984,5÷1986,5 (2600÷3800) Минимум 21-го цикла
3 1988,5÷1991,5 (4300÷5400) Максимум 22-го цикла
4 1994÷1997 (6300÷7500) Минимум 22-го цикла
5 1998÷2001 (8200÷9300) Максимум 23-го цикла
На рис. 6 представлена зависимость PSD(φ, Р) в виде изолиний для актив-
ных областей для эпох максимумов активности. Изолинии проведены через
375 относительных единиц. Аналогичные графики для корональных дыр пред-
ставлены на рис. 7. Интервал между изолиниями в этом случае 0,7 относи-
тельных единиц.
Проблемы управления и информатики, 2007, № 1 131
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
Активные области
П
ер
и
о
д
(
д
н
и
)
Максимум 21-го цикла
34
а
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
П
ер
и
о
д
(
д
н
и
)
Максимум 22-го цикла
34
б
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
Максимум 23-го цикла
П
ер
и
о
д
(
д
н
и
)
Широта (в градусах)
34
в
Рис. 6
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
П
ер
и
о
д
(
д
н
и
)
Максимум 21-го цикла
34
– 80 80
Корональные дыры
а
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
П
ер
и
о
д
(
д
н
и
)
Максимум 22-го цикла
34
– 80 80
б
0 – 60 – 40 – 20 20 40 60
22
24
26
28
30
32
П
ер
и
о
д
(
д
н
и
)
Максимум 23-го цикла
34
– 80 80
Широта (в градусах)
в
Рис. 7
132 ISSN 0572-2691
Анализ полученных распределений для активных областей показал, что в
максимумах и минимумах трех рассмотренных циклов вращение различается по
интервалам встречающихся периодов вращения, значению периодов с макси-
мальной мощностью в каждой полусфере и величине широтной зоны, в которой
наблюдается вращение с наиболее мощным периодом. Эти характеристики при-
ведены в табл. 3.
Таблица 3
Временнóй
интервал
Интервалы периодов
с PSDe(φ, Р) >750
для максимумов
и PSDe(φ, Р) >75
для минимумов, сут
Период вращения
с максимальным
PSDe(φ, Р), сут
Широтная зона вращения
с периодом с максималь-
ным PSDe(φ, Р), градусы
S-полу-
сфера
N-полу-
сфера
S-полу-
сфера
N-полу-
сфера
S-полу-
сфера
N-полу-
сфера
Максимум
21-го цикла
20÷35 24÷33 28,2 27 – 30÷0 –10÷35
Максимум
22-го цикла
22÷34 24÷30 26 29 – 40÷– 5 10÷45
Максимум
23-го цикла
22÷30 21÷34 28,5 26,8 – 40÷– 15 10÷35
Минимум
21-го цикла
22÷34 22÷34 27,3 26,5 – 35÷0 0÷35
Минимум
22-го цикла
22÷29 22÷32 28 28 – 20÷5 10÷20
В табл. 4, аналогичной табл. 3, приведены данные о вращении корональных дыр.
Таблица 4
Временнóй
интервал
Интервалы периодов
с PSDe(φ, Р) >0,7
для максимумов
и PSDe(φ, Р) >0,5
для минимумов, сут
Период вращения
с максимальным
PSDe(φ, Р), сут
Широтная зона вращения
с периодом с максималь-
ным PSDe(φ, Р), градусы
S-полу-
сфера
N-полу-
сфера
S-полу-
сфера
N-полу-
сфера
S-полу-
сфера
N-полу-
сфера
Максимум
21-го цикла
24÷35 22÷35 28,5 27 – 80÷80 – 10÷80
Максимум
22-го цикла
29÷32 24÷35 31 29 – 75÷–30 20÷40
Максимум
23-го цикла
20÷35 22÷35 27 30 – 20÷15 50÷80
Минимум
21-го цикла
25÷34 26÷30 27 26,3 –75÷– 10 35÷75
Минимум
22-го цикла
26÷31 26÷32 27,8 28,1 – 75÷– 35 50÷70
Обширные области, в которых наблюдается вращение с наиболее мощным
периодом (см. два последние столбца табл. 3 и табл. 4), нельзя считать зоной
твердотельного вращения, так как на тех же широтах существуют и другие перио-
ды, сравнимые по мощности с наиболее мощным периодом.
Обращает на себя внимание малое количество корональных дыр в 22-м цикле.
Исследование характеристик вращения с разрешением по времени в один
год по результатам скользящего спектрального анализа. По исходным времен-
ным рядам строились спектры мощности для текущего годичного интервала спек-
трального анализа с последующим сдвигом на полгода для каждой пятиградусной
широтной зоны. Таким образом, из одного временнóго ряда было получено
52 спектра мощности, относящиеся к последовательным годичным интервалам,
каждый из которых был смещен на полгода относительно предыдущего интерва-
ла. Таким образом, была построена зависимость PSD(φ, P, t).
Проблемы управления и информатики, 2007, № 1 133
На рис. 8 дан пример спектров мощности, полученных для широты N 20°
для шести последовательных годичных интервалов из участка 912÷1825 сут 26-лет-
него временнóго ряда для активных областей. Близкие периоды в соседних вре-
менных интервалах соединены пунктирной линией.
0
3
6
18
20 22 24 26 28 30 32 34 36
Период (дни)
912 d
1095 d
1277 d
1460 d
1642 d
1825 d
Широта 20°
PSD×10
–3
0
3
6
0
3
6
0
3
6
0
3
6
0
3
6
Рис. 8
Из рис. 8 видно, что на одной широте находится несколько больших пиков
спектральной плотности. Значения периодов, соответствующие этим пикам, ме-
няются со временем.
Анализируя спектры мощности, построенные по годичным временным рядам
для отдельных пятиградусных широтных зон, было установлено, что в них можно
выделить четыре последовательных интервала периодов длиной трое суток, в ко-
торые попадает большинство пиков, представленных в табл. 5. В этих интервалах
периодов определялась так называемая суммарная мощность процесса в виде за-
висимости ,),,(PSD),(1
2
1
=
P
P
tPtS которая лишь отсутствующим множителем
f (шаг дискретизации по частоте) отличается от средней мощности процесса
в интервале периодов ].;[ 21 PP
Таблица 5
№ п /п Интервал периодов Обозначение
1 2421 P 21÷24
2 2724 P 24÷27
3 3027 P 27÷30
4 3330 P 30÷33
Для активных областей и корональных дыр проведено суммирование зависи-
мостей ),(1 tS по всем дискретным значениям широты φ. Полученные зависимо-
сти обозначены как ).(2 tS На рис. 9 приведены зависимости )(2 tS для активных
областей и корональных дыр. Просуммированные зависимости )(2 tS по всем че-
тырем приведенным в табл. 5 интервалам периодов, представленные как ),(3 tS
для активных областей и корональных дыр отображены на рис. 10.
134 ISSN 0572-2691
0
0
20000
2000 4000 6000 8000 10000
Время (дни)
S
2
(t
)
40000
Активные области
1 — P 21–24
2 — P 30–33
1
2
50000
30000
10000
0
0
20000
2000 4000 6000 8000 10000
Время (дни)
S
2
(t
)
120000
Активные области
1 — P 24–27
2 — P 27–30
1
2
40000
60000
80000
100000
0
0
200
2000 4000 6000 8000 10000
Время (дни)
S
2
(t
)
400
Корональные дыры 1 — P 21–24
2 — P 30–33
1
2
0
0
100
2000 4000 6000 8000 10000
Время (дни)
S
2
(t
)
Корональные дыры
1 — P 24–27
2 — P 27–30
1
2
200
300
400
500
600
700
800
Рис. 9
0
0
10000
20000
30000
40000
2000 4000 6000 8000 10000
Время (дни)
S
3
(t
)
50000
P 21–33
0
0
100
200
300
400
2000 4000 6000 8000 10000
Время (дни)
500
P 21–33
600
700
S
3
(t
)
Рис. 10
Обсуждение результатов. Анализ спектров мощности, полученных по вре-
менным рядам за весь 26-летний интервал наблюдений, показал, что вращение ак-
тивных областей и корональных дыр в среднем за три 11-летних цикла суще-
ственно различается как набором значимых периодов вращения, так и их широт-
ным распределением.
Переход к более мелким временным масштабам — эпохам максимумов и ми-
нимумов цикла и годичным интервалам, позволил выявить достаточно быстрые
изменения в глобальных характеристиках вращения рассматриваемых солнечных
образований.
Спектры мощности для последовательных циклов солнечной активности суще-
ственно различаются как для эпох высокой активности, так и для ее минимумов.
Проблемы управления и информатики, 2007, № 1 135
Анализируя изменение со временем и широтой суммарной мощности спек-
тров внутри выделенных четырех периодов вращения, приходим к следующим
выводам.
● Суммарная мощность во всех четырех интервалах периодов для активных
областей в рассматриваемых циклах активности достигает примерно одинакового
максимума, но характер изменения ее со временем существенно меняется от цикла
к циклу. Так в 21- и 23-м циклах наблюдаются несколько максимумов, а в 22-м кри-
вая изменения суммарной мощности имеет один максимум.
● Для корональных дыр циклы существенно различаются максимальными
величинами и видом кривых изменения суммарной мощности со временем.
● Наибольшая суммарная мощность приходится на интервалы периодов
вращения 24÷27 и 27÷30 сут. Это справедливо и для активных областей, и для ко-
рональных дыр.
● Соотношение между суммарными мощностями в названных выше интер-
валах периодов заметно меняется со временем.
● Дифференциальность вращения слабо проявляется в широтном распреде-
лении периодов вращения. Одни и те же периоды встречаются в больших широт-
ных интервалах.
● В одном и том же широтном интервале одновременно присутствует вра-
щение с разными периодами.
Таким образом, вращение солнечных образований определяется, в основном,
подфотосферными источниками конкретных образований, а не вращением основ-
ных объемов солнечной плазмы, определяющих невозмущенное Солнце.
В работе использовались данные NSO/Kitt Peak, полученные как совместный
продукт NSF/NOAA, NASA/GSFC и NOAA/SEL.
Я.І. Зєлик, Н.М. Степанян, О.А. Андрєєва
ЧАСОВІ ЗМІНИ ГЛОБАЛЬНИХ
ХАРАКТЕРИСТИК ОБЕРТАННЯ СОНЦЯ
З використанням методів спектрального аналізу отримано глобальні характери-
стики обертання спостережуваних у лінії НеI 1083 λ нм сонячних структур для
трьох видів аналізу: на всьому часовому інтервалі спостережень тривалістю
близько 26 років; на п’ятьох часових інтервалах, що відповідають епохам мак-
симумів і мінімумів сонячної активності; на ковзному річному часовому інтер-
валі зі зсувом на півроку.
Yu.I. Zyelyk, N.N. Stepanian, O.A. Andryeyeva
TIME VARIATIONS OF GLOBAL
CHARACTERISTICS OF SUN ROTATION
By using the spectral analysis methods the global characteristics of the solar struc-
tures rotation, observed in line НеI 1083 λ nm are obtained for three types of analy-
sis: on the whole observation time interval of duration about 26 years; on five time
intervals corresponding to maxima and minima of solar activity; epoch on the sliding
one year time interval with half-year shift.
1. Зелык Я.И., Степанян Н.Н. Андреева О.А. О методах оценивания вращения солнечных
структур по наблюдениям Солнца в линии HeI λ 1083 нм // Космічна наука і технологія. —
2006. — 12, № 1. — С. 85–89.
2. Зелык Я.И., Степанян Н.Н., Андреева О.А. О применении аппарата анализа временных ря-
дов для оценивания вращения солнечных структур в верхней хромосфере // Проблемы
управления и информатики. — 2006. — № 3. — С. 102–115.
Получено 19.10.2006
|
| id | nasplib_isofts_kiev_ua-123456789-206965 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 0572-2691 |
| language | Russian |
| last_indexed | 2025-12-07T16:21:23Z |
| publishDate | 2007 |
| publisher | Інститут кібернетики ім. В.М. Глушкова НАН України |
| record_format | dspace |
| spelling | Зелык, Я.И. Степанян, Н.Н. Андреева, О.А. 2025-09-27T07:54:06Z 2007 Временные изменения глобальных характеристик вращения Солнца / Я.И. Зелык, Н.Н. Степанян, О.А. Андреева // Проблемы управления и информатики. — 2007. — № 1. — С. 127-135. — Бібліогр.: 2 назв. — рос. 0572-2691 https://nasplib.isofts.kiev.ua/handle/123456789/206965 523.9-327:519.246.8:517.443:519.71 З використанням методів спектрального аналізу отримано глобальні характеристики обертання спостережуваних у лінії НеI 1083 λ нм сонячних структур для трьох видів аналізу: на всьому часовому інтервалі спостережень тривалістю близько 26 років; на п’ятьох часових інтервалах, що відповідають епохам максимумів і мінімумів сонячної активності; на ковзному річному часовому інтервалі зі зсувом на півроку. By using the spectral analysis methods the global characteristics of the solar structures rotation, observed in line НеI 1083 λ nm are obtained for three types of analysis: on the whole observation time interval of duration about 26 years; on five time intervals corresponding to maxima and minima of solar activity; epoch on the sliding one year time interval with halfyear shift . ru Інститут кібернетики ім. В.М. Глушкова НАН України Проблемы управления и информатики Космический мониторинг Временные изменения глобальных характеристик вращения Солнца Часові зміни глобальних характеристик обертання Сонця Time variations of global characteristics of Sun rotation Article published earlier |
| spellingShingle | Временные изменения глобальных характеристик вращения Солнца Зелык, Я.И. Степанян, Н.Н. Андреева, О.А. Космический мониторинг |
| title | Временные изменения глобальных характеристик вращения Солнца |
| title_alt | Часові зміни глобальних характеристик обертання Сонця Time variations of global characteristics of Sun rotation |
| title_full | Временные изменения глобальных характеристик вращения Солнца |
| title_fullStr | Временные изменения глобальных характеристик вращения Солнца |
| title_full_unstemmed | Временные изменения глобальных характеристик вращения Солнца |
| title_short | Временные изменения глобальных характеристик вращения Солнца |
| title_sort | временные изменения глобальных характеристик вращения солнца |
| topic | Космический мониторинг |
| topic_facet | Космический мониторинг |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/206965 |
| work_keys_str_mv | AT zelykâi vremennyeizmeneniâglobalʹnyhharakteristikvraŝeniâsolnca AT stepanânnn vremennyeizmeneniâglobalʹnyhharakteristikvraŝeniâsolnca AT andreevaoa vremennyeizmeneniâglobalʹnyhharakteristikvraŝeniâsolnca AT zelykâi časovízmíniglobalʹnihharakteristikobertannâsoncâ AT stepanânnn časovízmíniglobalʹnihharakteristikobertannâsoncâ AT andreevaoa časovízmíniglobalʹnihharakteristikobertannâsoncâ AT zelykâi timevariationsofglobalcharacteristicsofsunrotation AT stepanânnn timevariationsofglobalcharacteristicsofsunrotation AT andreevaoa timevariationsofglobalcharacteristicsofsunrotation |