利用EMD方法提取太阳活动周期成分

EMD(经验模态分解)方法在处理非线性及非平稳时间序列时表现出了很大的优势和应用潜力.利用EMD方法研究太阳活动周期,对110年(1894-2003)和55年(1949-2003)的太阳黑子数月均值进行分解,分别得到一系列模式和一个趋势项,其中都可能包含有1.3至1.4年周期分量,25至30个月QBO(准双年振荡)分量,11年太阳周分量和22年Hale周分量.其中11年周期分量幅度最大,变化特征与太阳黑子数原始数据具有很高的相似性.不同于传统方法,EMD方法给出了太阳活动在不同时间尺度上各自分离的变化特征.      

太阳黑子数及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小波变换图中体现出来。     

武汉站第23周电离层TEC与太阳及地磁活动的相关性分析

利用武汉电离层观测站1997-2007年电离层TEC资料, 采用连续小波变换和交叉小波以及小波相干方法, 分析了该站电离层TEC的周期变化特征以及与太阳和地磁活动之间的关系. 分析结果表明, 武汉站TEC变化的长期趋势主要随太阳活动的强弱而变化; 在局部时域上分别存在128～256d, 256～512d和512d～1024d的周期尺度, 且与同时期的太阳黑子数和地磁Dst指数的周期特征存在很好的对应关系; 太阳黑子数在512～1024d周期尺度上超前TEC变化约1/6个周期; 在准半年的周期尺度上武汉站TEC与地磁Dst指数几乎呈反相位变化, 但TEC对$Dst$指数的这种响应仅在太阳活动高年存在, 具体机理尚需进一步分析研究.      

地磁Ap指数滞后太阳周循环分析

把1932-2006年地磁Ap指数12个月流动均值分解成为(Ap)R和(Ap)I.其中(Ap)R为太阳黑子数R的线性函数,与太阳黑子数R相位相同,可能对应于日冕物质抛射(CME)等地磁控制因素. (Ap)I分量与太阳黑子数R相位相差约180°,该分量可能对应于极冕洞变化(从太阳峰年开始,由日面极区逐渐向赤道延伸).以地磁Ap指数与太阳黑子数R滞后非常严重的第20太阳周为例,证实了(Ap)I分量与极冕洞向赤道延伸循环变化相对应.因此极冕洞循环变化可能是导致地磁扰动指数与太阳周循环相位不一致,出现滞后现象的一个十分重要原因.      

基于小波与交叉小波分析的太阳黑子与宇宙线相关性研究

利用小波分析和交叉小波分析方法, 根据太阳黑子数以及Huancayo和Climax两个测站的月均宇宙线数据, 分析了两个测站的月均宇宙线周期变化, 同时利用太阳黑子数R₁₂对Climax站宇宙线流量进行预测研究. 小波分析结果表明, 太阳黑子与宇宙线除存在显著的11年周期外, 太阳活动高年期间还存在1～6个月尺度的周期特性, 在第22太阳周活动高年时还出现了6～8和1～22个月的变化周期; 交叉小波分析结果表明, 在130个月左右的周期上宇宙线与太阳黑子具有显著的负相关性, 并且宇宙线的变化滞后太阳黑子约8个月; 分别采用预测时刻和8个月前的太阳黑子数, 预测相对误差为3.8912%和3.2386%. 本文方法同样适用于估算其他空间天气参量之间的周期和相关性, 提高各种空间天气参量的预测或预报精度.      

第23太阳活动周武汉站电离层TEC特征分析

利用武汉站(30.5°N, 114.4°E)1997年1月1日至2007年12月31日电离层TEC、太阳黑子数及地磁指数等资料, 分析了第23周武汉站TEC的周日变化、季节变化、半年变化以及与太阳活动的相关性等特征; 以2006年4月13-17日发生的磁暴为例, 讨论了武汉站TEC对磁暴的响应以及可能的机理. 结果表明,武汉站电离层TEC在太阳活动高、低年均呈典型的周日变化特征; 冬季异常和半年异常特征明显, 且受太阳活动强弱影响; TEC和太阳黑子数年均值相关系数为0.9611; TEC对磁暴的响应可能是由磁层穿透电场和中性风共同作用导致的, 具体影响机制有待深入研究.      

基于经验模态分解的地球同步轨道高能电子通量预报

在磁暴恢复相期间,大量相对论(高能)电子从磁层的外辐射带渗透到地球同步轨道区.其中> 2 MeV的高能电子能够穿透卫星表面并聚积在材料内部,导致卫星无法正常运行或完全损坏.磁暴期间的高能电子通量变化的非平稳与非线性特征十分明显.通过实验发现,经验模态分解法能够极大地降低高能电子通量非平稳性问题造成的预报影响.以2008-2009年的数据作为训练集,2010-2013年数据作为测试集.结果表明:2010-2013年的预报率约为0.84;在太阳活动较为复杂的2013年,预报率达到0.81.引入经验模态分解后预报效率得到显著提高.     

北京和东北地区雷暴与太阳活动的关系

太阳活动对于大气电性能可能会产生影响,而大气电性能在雷暴形成过程中可能会起着重要作用.因此,太阳活动对于大气电性能的影响,最终将有可能影响到雷暴的形成.本文用1957年到1978年北京地区和东北地区的10个气象站的13000多个雷暴资料,用时序迭加法分析了雷暴数与太阳磁扇形界面通过、太阳耀斑及太阳黑子11年周期变化等的关系.结果表明,有些年分的雷暴数当行星际磁场方向由指向太阳变为离开太阳的磁扇形界面扫过地球时有相关影响,在界面通过前后3天内或7天后雷暴数明显增大.从季节来看,上半年的界面通过比下半年的界面通过对雷暴的影响大.在各种太阳耀斑分类统计分析中,看到出现在日面上西区特别是西一区(0°-30°)的耀斑对雷暴的影响显著.雷暴数与太阳黑子11年周期变化则没有很显著的相关关系.      

电离层N_mF_2的准27天变化

利用1969-1980年期间东亚和澳大利亚扇区不同地磁纬度11个台站的电离层垂直探测数据,采用带通滤波的方法,分析了白天NmF2的准27天变化特征及其与太阳活动指数F10.7的相关关系.结果表明,在大多数年里,太阳活动指数F10.7和NmF2的短周期(2～70天)频谱中,27天周期附近都出现明显的极大值;在27天波段,F10.7和NmF2的标准偏差有逐年变化特征,F10.7的平均标准偏差为10.9%,NmF2 的标准偏差随地磁纬度变化,赤道地区最小,纬度越高标准偏差越大,11个站的平均标准偏差为8.2%;在27天波段,NmF2与F10.7存在显著的相关,在0.05的显著水平下,显著相关的概率在90%以上;NmF2相对于F10.7的准27天变化平均滞后2天左右;从总体上看,太阳EUV辐射的准27天变化是造成NmF2准27天变化的主要原因.     

太阳活动对电离层TEC变化影响分析

为研究太阳活动对电离层TEC变化的影响,从整体到局部分析了2000-2016年的太阳黑子数、太阳射电流量F_(10.7)指数日均值与电离层TEC的关系,并重点分析了2017年9月6日太阳爆发X9.3级特大耀斑前后15天太阳活动与电离层TEC变化的相关性.结果表明:由2000-2016年的数据整体看来,太阳黑子数、太阳F_(10.7)指数、TEC两两之间具有很强的整体相关性,但局部相关性强弱不均;此次耀斑爆发前后太阳黑子数、太阳F_(10.7)指数和TEC具有很强的正相关特性,太阳活动对TEC的影响时延约为2天;太阳活动对全球电离层TEC的影响不同步,从高纬至低纬约有1天的延迟,且对低纬度的影响远大于中高纬度.太阳活动是影响电离层TEC变化的主要原因,但局部也可能存在其他重要影响因素.     

Phase asynchrony between coronal index and sunspot numbers

In this paper, the phase asynchrony between coronal index and sunspot numbers is investigated. It is found that, (1) the sunspot numbers begin one month earlier than coronal index, which should mathematically lead to phase asynchrony between them but with a slight effect; (2) the 11-year Schwabe cycle is the only one period with statistical significance for coronal index and sunspot numbers, and the difference between the length of the Schwabe cycle of them should also lead to phase asynchrony between them; (3) although coronal index and sunspot numbers are coherent in low-frequency components corresponding to the 11-year Schwabe cycle, they are asynchronous in phase in high-frequency components; (4) their different definitions and physical meanings may be a major reason why there is a phase asynchrony between them.     

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.     

A 22-year cycle in sunspot activity

We study the recently presented group sunspot number series and show that a persistent 22-year periodicity exists in sunspot activity throughout the entire period of about 400 years of direct sunspot observations. The amplitude of this periodicity in total cycle intensity is about 20% of the present intensity level. A 22-year periodicity in sunspot activity is naturally produced by the 22-year magnetic dynamo cycle in the presence of a relic magnetic field. Accordingly, a persistent 22-year periodicity in sunspot activity gives strong evidence for the existence of such a relic magnetic field in the Sun. The stable phase and the roughly constant amplitude of this periodicity during times of very different sunspot activity level strongly support this interpretation.     

Forecasting the 23 and 24 solar cycles on the basis of MGM spectrum

A new method of nonlinear spectral analysis (called the method of global minimum: MGM), based on the best presentation (in sense of minimal squares) of a given time data set as a sum of sinusoids whose frequencies, amplitudes and phases are to be determined, has been used to find periodicities in annual Wolf sunspot numbers (W) during the period 1700–1995. The possible future behaviour of the 11-year solar cycle (based on an extrapolation of the calculated model) is also presented. The main characteristics of the 23^rd solar cycle are as follows: the W maximum occurs about 2004, with a peak of nearly 220. An unusually large value of W will occur during the 23^rd cycle, which should be characterised by the longest maximum, specifically, W will be greater than 100 during the 11-year period from 1997 to 2007. The first sharp rise will occur during the period 1996–1998, the second sharp during 2002–2004. The main features of the 24 year cycle are as follows: the next minimum in W, associated with the 24^th solar cycle, should occur in the year 2008 and the maximum in 2014. W is expected to peak at about 180. The minimum value for the 25^th year cycle is expected to occur in the year 2019. It is shown that the accuracy of these predictions depends, first of all, on the extrapolation of the hyperlong harmonic of the calculated polyharmonic model fit of observed annual sunspot numbers during the period 1700–1995. The error bars in the definition of the maximum and minimum epochs can be as large as two years.     

Evolution of the solar corona from solar maximum to minimum

A study on solar coronal activities related to the 11-year activity cycle is presented from the Yohkoh soft X-ray observations. Yohkoh was launched in August 1991, just after the solar maximum of the cycle 22 and continues to observe the Sun in the declining phase of the magnetic activity cycle toward the solar minimum. The soft X-ray flux from the whole Sun in the declining phase essentially decreases with the size of active regions. The X-ray intensity in quiet regions in the declining phase decreases with the magnetic flux observed at the photosphere. The whole-Sun soft X-ray flux does not monotonically decrease, but there are periodic enhancements of the flux with about a one-year interval. The activity appears as bright clusters in the butterfly diagram of the soft X-ray intensity and corresponds to the emergence of complexes of activity in the sunspot zones. The high-latitude activity is also studied, and we find that the X-ray intensity of high-latitude regions fluctuates with time scale of about one year.     

高层大气原子氧的半年周变

地球高层大气成分的长期变化受太阳黑子周、太阳活动程度和地磁活动程度等诸多因子所控制．本文利用国外有关高层大气成分的数据,分别讨论受上述控制因子影响的高层大气成分长期变化,讨论范围仅限原子氧半年周期变化．选用LDEF在轨飞行器1984年4月—1990年1月高度470km附近的长期资料进行统计分析,结果表明,高度470km附近原子氧在年平均太阳黑子数＜20、太阳活动程度相对低而平稳期间,半年周期的变化尤为明显,相对变幅约为40％—60％、井随平均太阳黑子数增加而增大．而年平均太阳黑子数峰值的1989年期间（＞120）,半年变化的相对变幅可达87％左右．     

Trends in foF2 and the 24th solar activity cycle

The trends in foF2 are analyzed based on the data of Juliusruh and Boulder ionospheric stations. It is shown that using the traditional solar activity index F10.7 leads to an impossible trend in foF2 when the data for the 24th solar activity cycle are included into the analysis. It is assumed that the F10.7 index does not describe correctly the solar ultraviolet radiation variations in that cycle. A correction of this index using the Rz (sunspot number) and Ly (intensity of the Lyman-α line in the solar spectrum) is performed, and it is shown that in that case reasonable values of the foF2 trends are obtained.