The MGM spectral analysis of Wolf sunspot numbers: Preliminary results

A method of nonlinear spectral analysis, named by the method of the global minimum (MGM) has been used to find periodicities in the annual Wolf sunspot numbers W, for the 295 year time period from 1700 to 1995. The normalised mean square error of the model fit to the data set . The main spectral peak is an unresolved (due to confluence of error bars) triplet at mean period 252 yr. We hypothesise that one of the triplet peaks is connected with the rotation period of a massive 10^th planet, which has to be located at one of Pluto's Lagrange points. The solar magnetic cycle (22 yr) could be the 11^th harmonic of this triplet (22.3×11=245 yr). The most powerful peak in vicinity of 11 years is an unresolved doublet with a mean period of T=10.8 yr. We have demonstrated that the magnetic cycle is non-stationary during the period 1700–1995 and described by doublet with a mean period of T=22.3 yr. Other nonstationary harmonics are discussed, also. The 11.87- yr spectral peak could be associated with the effect of Jupiter. Error bars and the statistical significances of the spectral peaks are presented. The study is ongoing, so that our results will be rendered more precise in the next future.     

太阳黑子数及Ap指数周期变化特征的小波分析

应用Morlet小波变换方法从多个变化尺度上对1932—2000年的太阳黑子数及Ap指数的变化特征进行分析．(1)太阳黑子数存在准11年、准32年的周期变化特征及Ap指数存在准32年、准11年、准6个月、准27天和准13．9天的周期变化特征；(2)太阳黑子数及Ap指数有着相似的准11年周期变化，但Ap指数极值的出现要比太阳黑子数极值出现滞后1—2年；(3)Ap指数准27天的周期变化在太阳黑子活动高、低年不同，在太阳活动低年，Ap指数有着较稳定的准27天周期变化，但在太阳活动高年，27天的周期变化几乎消失，这种周期变化的消失和出现时间可在Morlet小波变换图中体现出来。     

Solar terrestrial effects of two distinct types

The occurrence frequencies or fluxes of most of the solar phenomena show a 11-year cycle like that of sunspots. However, the average characteristics of these phenomena may not show a 11-year cycle. Among the terrestrial parameters, some related directly to the occurrence frequencies of solar phenomena (for example, ionospheric number densities related to solar EUV fluxes which show 11-year cycle like sunspots) show 11-year cycles, including the double-peak structures near sunspot maxima. Other terrestrial parameters related to average characteristics may not show 11-year sunspot cycles. For example, long-term geomagnetic activity (Ap or Dst indices) is related to the average interplanetary solar wind speed V and the total magnetic field B. The average values of V depend not on the occurrence frequency of ICMEs and/or CIRs as such, but on the relative proportion of slow and high-speed events in them. Hence, V values (and Ap values) in any year could be low, normal or high irrespective of the phase of the 11-year cycle, except that during sunspot minimum, V (and Ap) values are also low. However, 2–3 years after the solar minimum (well before sunspot maximum), V values increase, oscillate near a high level for several years, and may even increase further during the declining phase of sunspot activity, due to increased influence of high-speed CIRs (corotating interplanetary regions). Thus, Ap would have no fixed relationship with sunspot activity. If some terrestrial parameter shows a 11-year cycle, chances are that the solar connection is through the occurrence frequencies (and not average characteristics) of some solar parameter.     

太阳风速度的周期结构分析

本文应用功率谱估计中的加窗法和最大熵法对第20太阳活动周太阳风速度的周期结构进行了频谱分析;为了克服太阳活动高年瞬变扰动的干扰,本文还用互相关滤波方法与谱分析进行了比较。得到如下结果:除了个别年代外,在太阳活动周的所有相位里,太阳风速度几乎都存在13天和27天两种主要周期变化成份。在一些年里,还有9天的周期成份存在。      

Interplanetary shocks and sudden impulses during solar maximum (2000) and solar minimum (1995–1996)

In this work a study is performed on the correlation between fast forward interplanetary shock parameters at 1 Astronomical Unit and sudden impulse (SI) amplitudes in the H-component of the geomagnetic field, for periods of solar activity maximum (year 2000) and minimum (year 1995–1996). Solar wind temperature, density and speed, and total magnetic field, were taken to calculate the static pressures (thermal and magnetic) both in the upstream and downstream sides of the shocks. The variations of the solar wind parameters and pressures were then correlated with SI amplitudes. The solar wind speed variations presented good correlations with sudden impulses, with correlation coefficients larger than 0.70 both in solar maximum and solar minimum, whereas the solar wind density presented very low correlation. The parameter better correlated with SI was the square root dynamic pressure variation, showing a larger correlation during solar maximum (r = 0.82) than during solar minimum (r = 0.77). The correlations of SI with square root thermal and magnetic pressure were smaller than with the dynamic pressure, but they also present a good correlation, with r > 0.70 during both solar maximum and minimum. Multiple linear correlation analysis of SI in terms of the three pressure terms have shown that 78% and 85% of the variance in SI during solar maximum and minimum, respectively, are explained by the three pressure variations. Average sudden impulse amplitude was 25 nT during solar maximum and 21 nT during solar minimum, while average square root dynamic pressure variation is 1.20 and 0.86 nPa^1/2 during solar maximum and minimum, respectively. Thus on average, fast forward interplanetary shocks are 33% stronger during solar maximum than during solar minimum, and the magnetospheric SI response has amplitude 20% higher during solar maximum than during solar minimum. A comparison with theoretical predictions (Tsyganenko’s model corrected by Earth’s induced currents) of the coefficient of sudden impulse change with solar wind dynamic pressure variation showed excellent agreement, with values around 17 nT/nPa^1/2.     

Specific features of the high-speed plasma stream cycles

The high-speed plasma streams in the solar wind are investigated during the solar cycles nos. 20–22 (1964–1996), separately on the two types of streams according to their solar origin: the HSPS produced by coronal holes (co-rotating) and the flare-generated, in keeping with the classification made in different catalogues. The analysis is performed taking into account the following high-speed stream parameters: the durations (in days), the maximum velocities, the velocity gradients and, the importance of the streams. The time variation of these parameters and the high-speed plasma streams occurrence rate show an 11-year periodicity with some differences between the solar cycles considered. A detailed analysis of the high-speed stream 11-year cycles is made by comparison with the “standard” cycles of the sunspot relative number (Wolf number). The different behaviour of the high-speed stream parameters between even and odd solar cycles could be due to the 22-year solar magnetic cycle. The increased activity of the high-speed plasma streams on the descendant phases of the cycles, regardless of their solar sources, proves the existence of some special local conditions of the solar plasma and the magnetic field on a large scale that allow the ejection of the high energy plasma streams. This fact has led us to the analysis the stream parameters during the different phases of the solar cycles (minimum, ascendant, maximum and, descendant) as well as during the polar magnetic field reversal intervals. The differences between the phases considered are pointed out. The solar cycles 20 and 22 reveal very similar dynamics of the flare-generated and also co-rotating stream parameters during the maximum, descendant and reversal intervals. This fact could be due to their position in a Hale Cycle (the first component of the 22-year solar magnetic cycle). The 21st solar cycle dominance of all co-rotating stream parameters against the 20th and 22nd solar cycle ones, during almost all phases, could be due to the same structure of a Hale Cycle – solar cycle 21 is the second component in a 22-year SC. During the reversal intervals, all high-speed stream parameters have comparable values with the ones of the maximum phases of the cycles even if this interval contains a small part of the descendant branch (solar cycles 20 and 22).     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

利用相似周方法预测第24活动周开始时间

利用相似周方法对第24活动周的开始时间与第23活动周下降相后期的太阳黑子数进行了预报.根据第23周已经出现的特征参量和下降相的形态特征,选取9,10,11,15,17和20等六个太阳活动周作为第23周下降相的相似周,对第24周开始时间进行预报.预报结果显示,第24活动周的开始时间为2007年5±1月,黑子数平滑月均极小值为7.1±2.6,第23太阳活动周长度为11.1年.与其他研究者的预报结果相比较,本文给出的结果与文献[11]和[12]及MSFC的结果比较一致.通过对相似周方法在下降相预测太阳活动周结束时间的研究讨论,及对第23周上升阶段的太阳黑子数和F10.7平滑月均值预报结果的评估,可以看出,相似周预报方法在太阳活动周长期预报中是很有应用价值的.      

A study of the solar cycle signal on total ozone over low-latitude Brazilian observation stations

Observations of total ozone at low latitudes in Brazil have been made using Dobson spectrophotometers since 1974 for Cachoeira Paulista (23.1° S, 45° W) and since 1978 for Natal (5.8° S, 35.2° W). Annual averages, 12 months and 36 months running averages have been analyzed. Spectral analyses of the data revealed that the most important periods found (confidence level> 90%) were: for Natal, 2.5 years (93.1%, quasi-biennial oscillation-QBO) and 10 years (98,2%, possibly the solar cycle signal); for Cachoeira Paulista, 2.4 years (96.8%, QBO) and 8 years (99.6%). The difference in total ozone between maximum and minimum solar cycles were estimated, using yearly averages of total ozone. For solar cycle 21, 1.16% and 1.26% for Natal and Cachoeira Paulista were found; for solar cycle 22, a larger difference of 3.8% for Natal and 4.1% for Cachoeira Paulista were found. The corresponding variation in UV-B at 300 nm, using Beer's law, is 8–10% for C. Paulista and 4–5% for Natal, with maxima occurring during the minimum of the solar cycle.     

海口站与Huancayo站扩展F的差异

本文统计分析了我国海口站太阳黑子高年和低年的扩展F的资料,得出了该站扩展F的出现率随地方时、季节和太阳活动周期的某些变化规律并与Huancayo站扩展F的资料进行了比较,指出了两站扩展F的差异。      

Variations of the radiation dose onboard Mir station

Dose variations, associated with the 11-year solar activity cycle, seasonal variations of particle fluxes in the Earth's radiation belts at the station orbit, and solar proton events are studied, using prolonged measurements of radiation doses inside orbital station Mir. Daily averages of radiation doses during the declining phase of the 22^nd solar cycle and during transition to the 23^rd solar activity cycle reached very large values for astronauts and significantly exceed the values calculated according to existing models.     

Noise in wireless systems from solar radio bursts

Solar radio bursts were first discovered as result of their interference in early defensive radar systems during the Second World War (1942). Such bursts can still affect radar systems, as well as new wireless technologies. We have investigated a forty-year record of solar radio burst data (1960–1999) as well as several individual radio events in the 23rd solar cycle. This paper reviews the results of a portion of this research. Statistically, for frequencies f  1 GHz (near current wireless bands), there can be a burst with amplitudes >10³ solar flux units (SFU; 1 SFU = 10⁻²² W/m²) every few days during solar maximum conditions, and such burst levels can produce problems in contemporary wireless systems.     

Relationship between daily variation of cosmogenic nuclide Be-7 concentration in atmosphere and solar activities

Daily Be-7 concentrations in air at the height of 15 m are continuously observed at 38°15.2′N, 140°20.9′E, between 2000 and 2001. The average concentration and the relative standard deviation were 4.0 mBq/m³ and 50% in 2000–2001, respectively. The Be-7 concentrations increased 2.5% with the decrease in the sunspot numbers by 6.7% for the term of two years. From the power spectral analysis, the periodicity of 26 days is shown for the daily Be-7 concentrations. The folding analysis indicates that the time variation of the Be-7 concentration is similar to that of the ground-based neutron counting rate, and the phase delay for the minimum portion of Be-7 concentration was roughly 8 days to the maximum sunspot number. These results indicate that the Be-7 concentrations in the air at ground level have 26 day periodicity as a component of time variations and the time variation is caused by the solar modulation of galactic cosmic rays, which corresponds to the variation of the sunspot number due to the rotation of the sun.     

Sunspot numbers during 1736–1739 revisited

Some periods in the sunspot number reconstruction composed by Hoyt and Schatten [Hoyt, D.V., Schatten, K.H. Group Sunspot Numbers: a new solar activity reconstruction. Sol. Phys. 179, 189–219, 1998. Reprinted with figures in Sol. Phys. 181, 491–512, 1998], are based on very few records. For example, there are only a few solar observations during the years 1736–1739. In this paper we intend to improve the reliability of the sunspot numbers reconstruction developed by Hoyt and Schatten for this 4-years period based on information about solar activity published in three journals of that epoch: “Philosophical Transactions”, “Histoire de l’Académie Royale des Sciences”, and “Nova Acta Eruditorum”. We were able to identify 42 papers with solar observations, including 30 with relevant information on sunspots. Based upon this new outlook, a reconstruction of the monthly solar activity for these years is proposed.     

第21～24太阳周4类空间天气事件爆发特征统计分析

对第21～24太阳周不同等级的太阳X射线耀斑事件、太阳质子事件、地磁暴事件及高能电子增强事件的爆发频次特征进行统计,结果表明:太阳周耀斑爆发的总数量与该太阳周的黑子数峰值呈正比,耀斑总数、X级耀斑事件数与峰值的相关系数分别为0.974,0.997;太阳质子事件主要发生在峰年前后1～2年,约占总发生次数的80%,峰值通量大于10pfu （1 pfu=1 cm-2·sr-1·s-1）的质子事件中,84%伴有耀斑爆发,并且主要伴随M或X级耀斑,少量伴随C级耀斑,峰值通量大于1000pfu的质子事件中,98%伴随M或X级耀斑,并且以X级耀斑为主;第21,22,23和24太阳周发生地磁暴最频繁的时间分别在1982,1991,2003年和2015年,分别滞后黑子数峰值时间3年、2年、2年和1年;72%的高能电子增强事件发生在太阳周下降期,24%的高能电子增强事件发生在太阳周上升期.      

Temporal evolution of statistical features of the sunspot cycles

Analysis of the general statistical features of the sunspot cycles in the period 1700–1996 AD, including the Gnevyshev–Ohl rule, Waldmeier rule and an amplitude–period effect, was performed for both Wolf numbers and group sunspot numbers. It was shown that for both solar indices all the statistical effects are weaker over the time interval 1700–1855 AD than over the time interval 1856–1996 AD. Possible causes of this difference are discussed.     

The structure of the solar cycle maximum phase in the galactic cosmic ray 11-year variation

Two phenomena connected with the maximum phase of the 11-year solar cycle in the galactic cosmic ray intensity – the change in the energy dependence of the intensity variations and the double-peak structure in the intensity modulation time profile – are considered for the last five solar cycles (Nos. 19–23). The distinct 22-year cycle in the magnitude of the so called energy hysteresis is observed.The periods of the solar cycle maximum phase in the galactic cosmic ray intensity, characterized by the specific energy dependence of the intensity, are estimated. It is found that the double-peak structures belonging to the solar cycle maximum phase and those around it are very similar both in the amplitude and in its energy dependence.     

Even–odd cycle parity in complexity of hemispheric sunspot activity

We have used the Lempel–Ziv measure to describe the complexity in sunspot activity during the solar cycles 18–23. In particular, we examined the time series of daily sunspot numbers in the northern and southern hemispheres in each of the six cycles and calculated the Lempel–Ziv complexity (LZC) value for each time series. Our results indicate that in the even cycles, the LZC values of the sunspot numbers in the two hemispheres are very close to each other, whereas in the odd cycles they differ significantly between the two hemispheres. We also find that within each hemisphere the LZC varies from cycle to cycle. This even–odd cycle parity reflects the variations in inter-hemispheric strengths of the solar magnetic field leading to different temporal distributions of sunspots in the two hemispheres. The degree of complexity may influence the predictability of sunspot activity in the two hemispheres during the various cycles. Although the physical implication of the results is not clear, these results may stimulate new ideas into modeling the complex dynamics of the solar dynamo.     

Sunspot cycle 23 descent to an unusual minimum and forecasts for cycle 24 activity

The decay phase of the sunspot cycle 23 exhibited two unusual features. First, it lasted too long. Second, the interplanetary magnetic field intensity at earth orbit reached the lowest value since in situ measurements in space began in October 1963. These physical anomalies significantly altered the early forecasts for the sunspot activity parameters for cycle 24, made by several colleagues. We note that there was a significant change in the solar behavior during cycle 22. We discuss the observed trends and their effect on our empirical solar activity forecast technique, leading to our prediction for cycle 24 parameters; cycle 24 will be only half as active as cycle 23, reaching its peak in May 2013. We speculate on the possible implications of this outcome on future earth climate change and the ensuing socio-economic consequences.     

中等磁暴太阳周分布的统计研究

依据实际观测的中等磁暴数据,统计分析了中等磁暴的太阳周分布.分析结果表明,在一个太阳活动周内,每年中等磁暴随时间的变化出现多个峰值,其中,最大峰值均出现在太阳活动周的下降段,即中等磁暴的峰值比太阳黑子数平滑年均值的峰值要滞后,滞后的时间为2～3年.超过70%的中等磁暴出现在太阳活动周的下降段,这表明绝大多数中等磁暴出现在太阳活动周的下降段.通过对中等磁暴平滑月均值与太阳黑子数平滑月均值相位差的计算分析发现,中等磁暴峰值出现的时间比太阳黑子数峰值出现的时间要滞后,不同太阳活动周中等磁暴峰值出现的时间与太阳黑子数峰值时间滞后的程度不同.