1998年5月2日日冕亮度观测图的数值研究

采用数值手段模拟了1998年5月2日日冕亮度观测图，计算模式改进为球坐标下系下特殊的二维理想磁流体（MHD）模型，即把（r,φ）坐标建立在SOHO观测日冕亮度的子午面上，消除了子午面极区的几何奇异性，根据SOHD日冕观测布置磁极子得到初始磁场位形，内边界条件采用自治的投影特征边界条件，计算迭代出稳态的多 磁场结构，得得到了与观测基本一致的亮度图，计算结果表明在太阳表面附近磁场位形对太阳凤等离子参数分布起控制作用。     

WSA太阳风经验模型及其应用

Wang-Sheely-Arge （WSA）模型是对准稳态太阳风的经验和物理相结合的描述,其利用观测的日面磁图作为输入,可以提前3到4天预测L1点处的太阳风速度和行星际磁场极性.WSA模型是在WS模型的基础上经过若干改进形成的实时预报模式,之后又借鉴Distance from the Coronal Hole Boundary （DCHB）模型的参数,进一步改进了太阳风速度关系式,形成了目前常见的形式.WSA经验模型由日冕磁场模型、太阳风速度关系式和一维运动学模型三部分组成.在实际应用过程中,基本步骤包括观测磁图预处理、日冕三维磁场反演、计算日冕磁场参数、计算太阳风的速度分布和将太阳风映射到1AU等环节.在发展过程中,WSA模型经历了一些细节上的调整变化,例如观测磁图数据的来源、日冕磁场模型的类型、经验速度关系中自由系数的取值等.许多研究对如何改善模型的预报效果进行了探索.      

多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

基于CPU/GPU的日冕偏振亮度并行计算模型

三维磁流体力学(MHD)数值模拟是用来研究日冕和太阳风最常用的方法之一, 其中将计算得到的日冕电子数密度转化为日冕偏振亮度(Polarization Brightness, PB)是与观测对比的重要方法. 由于待转换电子数据网格密度、PB数据网格密度和计算模型的复杂度, 使得日冕偏振亮度的计算比较耗时, 利用单CPU计算无法达到近实时转换日冕偏振亮度的要求, 从而影响了数值模拟的验证效率. 本文在CPU/GPU环境下, 利用CUDA编程技术, 提出了一个日冕偏振亮度并行计算模型. 实验结果表明, 该模型比CPU上的串行模型计算速度提高了31.86倍, 达到了近实时模拟与观测数据比对的计算要求.      

1998年5月日冕三维结构的数值模拟

采用三维理想磙流体力学（MHD）模式，内边界条件把二维投影特征线边界方法推广应用到三维计算，有效地稳定了数值计算并保证稳态解的自洽性；初始猜解磁场由1935卡林顿周光球磁场观测数据得到，这样计算得到的1998年5月份期间日冕三维结构比较符合实际，计算结果表明：（1）计算得到的源表面非径向磁场量值在磁中性线附近不超过2μT，表明源表面磁场基本径向。（2）模拟得到的源表面径向磁场量值除了在磁中性线附近的区域外变化不大，这和观测一致。（3）由源表面磁场按平方反比的规律计算出1AU处磁场量值更接近观测值。（4）计算得到的日冕结构和观测定性一致，三维数值模拟结果表明，日冕的三维大尺度背景结构主要是由磁场决定的，在闭磁场处或者电流片附近，太阳风的密度高，速度低；在开场区，太阳风的密度低，速度高。     

二维太阳风速度结构与日球电流片

本文分析了第1733—1742Carrington周的太阳耀斑、行星际闪烁(IPS)、太阳风及K-日冕资料,对二维太阳风速度结构进行了综合研究。结果表明:速度的大尺度分布、高速区、速度梯度、经度不均匀性和边界形状等方面相对于太阳赤道面都是不对称的;通过对理想的和实际的电流片以及飞船测量的磁场和速度比较,可以看到电流片附近大体上为低速区,但不存在简单的对应关系;扰动太阳风对背景太阳风的二维速度结构有显著的影响;理想电流片与实际电流片之间也存在较大的差距。      

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

利用多卫星观测CIR事件演变的统计分析

根据2007-2009年STEREO-BEHIND (STB)和ACE卫星的行星际磁场和太阳风数据, 基于冕洞高速流从太阳向外匀速径向传输假设, 讨论了随着STB和ACE卫星与太阳之间的夹角从0°增大至70°时, 冕洞发出的高速太阳风形成的相互作用区(CIR)依次扫过STB和ACE卫星的时间差特性, 并统计分析了两颗卫星观测到的CIR参数的变化特征. 结果显示, 可以利用STB对CIR事件的观测来预测这个CIR事件到达ACE的理论时间, 时间误差均值和最大值分别为0.217d和0.952d, 时间误差的产生与STB和ACE卫星观测到的CIR速度大小的不同有关, 用速度差异矫正后, 时间误差的平均值和最大值可分别减小为0.194d 和0.489d; STB和ACE卫星观测的CIR事件太阳风速度最大值的线性相关系数达到了0.84, STB和ACE卫星观测到的CIR事件对特征物理量中速度、质子温度的变化最小, 而质子密度及总压力的变化最大. 分析结果表明, STB和ACE卫星观测到的CIR事件有很强的相似性, STB卫星的CIR观测可以作为ACE卫星观测CIR事件特征的参考, 从而为地球空间环境扰动预报提供依据.      

太阳风速度及日球电流片对CME渡越时间的影响

基于1996-2005年88个引起重大地磁暴的CME(日冕物质抛射)事件、1996-2000年的47个CME事件以及1997-2002年的29个全晕状CME事件,结合ACE卫星在1AU处的太阳风和行星际磁场观测资料以及Wilcox Solar Observatory(WSO)天文台的太阳光球层磁图,分析了背景太阳风速度和日球电流片对CME到达1AU处渡越时间预报误差的影响.结果表明,背景太阳风速度与CME渡越时间误差并没有明显的相关性,在考虑了磁云通量管轴相对黄道面夹角的影响后相关性依然不明显.然而日球电流片对CME渡越时间却有明显的影响,对于初速度较小的异侧CME事件,其渡越时间大于同侧事件;而对于具有较大初速度的CME事件,异侧事件的渡越时间明显小于同侧事件.研究结果表明,CME与太阳风以及日球电流片的相互作用并不是简单的对流相互作用,造成高速CME异侧事件快于同侧事件到达地球的因素非常复杂,有待深入研究.      

ICME期间磁暴的太阳风参数和极光沉降能量

基于1995-2004年ICME驱动的强烈磁暴(SA型)、强磁暴(SB型)和延迟型主相暴(SC型)三种磁暴类型,对1AU处太阳风动压、太阳风速度、行星际磁场、EK-L电场以及极光沉降能量进行时序叠加分析,并分别与-vBz耦合函数和Newell耦合函数进行对比.结果表明,三种磁暴在ICME到达前期的太阳风动压较稳定,背景太阳风、极光沉降能量、行星际磁场和磁层存在相对平静期. ICME到达前期SA型磁暴的背景太阳风速度、行星际磁场南向分量以及极光沉降能量的均值高于另外两种磁暴类型,这说明大型日冕物质抛射在ICME到达前就对行星际磁场、背景太阳风和HP产生了影响.磁暴急始后,SC型磁暴的EK-L电场斜率小,峰值延后且行星际磁场北向分量增强,这些都是磁暴主相延迟的表现,极光沉降能量随着行星际磁场转为南向而增加.     

Dynamic processes in the outer heliosphere: voyager observations and models

We review recent Voyager 2 observations in the vicinity of 70 AU. The character of the solar wind plasma data between 2002 and 2003 changed to a regime in which the speed, density and magnetic field magnitude are positively correlated. The average speed of the solar wind at Voyager 2 increased between early 2003 and mid-2004, which we attribute to a return of fast coronal hole flow. We use solar wind data at Earth as input to numerical models which include the effect of pickup ions to model the radial evolution of the solar wind. The model reproduces the basic features of the observations. As a specific example, we investigate the propagation of the Halloween (Oct.–Nov.), 2003 storms in the outer heliosphere. The model predictions are in reasonable agreement with Voyager 2 observations.     

Relation between solar wind velocity and properties of its source region

Different kinds of coronal holes are sources of different kind of solar winds. A successful solar wind acceleration model should be able to explain all those solar winds. For the modeling it is important to find a universal relation between the solar wind physical parameters, such as velocity, and coronal physical parameters such as magnetic field energy. To clarify the physical parameters which control the solar wind velocity, we have studied the relation between solar wind velocity and properties of its source region such as photospheric/coronal magnetic field and the size of each coronal hole during the solar minimum. The solar wind velocity structures were derived by using interplanetary scintillation tomography obtained at Solar-Terrestrial Environment Laboratory, Japan. Potential magnetic fields were calculated to identify the source region of the solar wind. HeI 1083 nm absorption line maps obtained at Kitt Peak National Solar Observatory were used to identify coronal holes. As a result, we found a relation during solar minimum between the solar wind velocity and the coronal magnetic condition which is applicable to different kind of solar winds from different kind of coronal holes.     

Two-velocity structure observed in the inner solar wind

Measurements of the motion of plasma density inhomogeneities in the inner solar wind are presented. The speeds were estimated using a cross-correlation analysis of radio frequency fluctuations of the Galileo spacecraft measured simultaneously at widely spaced ground stations. The radial projections of the correlation baselines on the pattern plane were of the order of several thousand kilometers. For cross-correlation functions calculated with comparatively short averaging times, we find that a pronounced two-velocity configuration is occasionally observed over the range of heliocentric distances 20 R < R < 40 R. The typical mean speed for such observations is about 300–400 km/s and the difference between the two predominant speeds is about 150–200 km/s. These results may indicate that the density fluctuations are associated with slow magnetosonic waves propagating in opposite directions at the local speed of sound in the reference frame moving with the mean solar wind speed. Quite reasonable estimates of the solar wind speed and speed of sound are obtained from this model. Another possible explanation of the two-velocity structures is that two independent solar wind streams are present simultaneously along different segments of the radio ray path.     

The characteristics of sharp (small-scale) boundaries of solar wind plasma and magnetic field structures

We investigate properties of large (>20%) and sharp (<10 min) solar wind ion flux changes using INTERBALL-1 and WIND plasma and magnetic field measurements from 1996 to 1999. These ion flux changes are the boundaries of small-scale and middle-scale solar wind structures. We describe the behavior of the solar wind velocity, temperature and interplanetary magnetic field (IMF) during these sudden flux changes. Many of the largest ion flux changes occur during periods when the solar wind velocity is nearly constant, so these are mainly plasma density changes. The IMF magnitude and direction changes at these events can be either large or small. For about 55% of the ion flux changes, the sum of the thermal and magnetic pressure are in balance across the boundary. In many of the other cases, the thermal pressure change is significantly more than the magnetic pressure change. We also attempted to classify the types of discontinuities observed.     

Solar wind structure in the outer heliosphere

A solar wind parcel evolves as it moves outward, interacting with the solar wind plasma ahead of and behind it and with the interstellar neutrals. This structure varies over a solar cycle as the latitudinal speed profile and current sheet tilt change. We model the evolution of the solar wind with distance, using inner heliosphere data to predict plasma parameters at Voyager. The shocks which pass Voyager 2 often have different structure than expected; changes in the plasma and/or magnetic field do not always occur simultaneously. We use the recent latitudinal alignment of Ulysses and Voyager 2 to determine the solar wind slowdown due to interstellar neutrals at 80 AU and estimate the interstellar neutral density. We use Voyager data to predict the termination shock motion and location as a function of time.     

Prompt effects of solar wind variations on the inner magnetosphere and midlatitude ionosphere

It is well known that the solar wind can significantly affect high-latitude ionospheric dynamics. However, the effects of the solar wind on the middle- and low-latitude ionosphere are much less studied. In this paper, we report observations that large perturbations in the middle- and low-latitude ionosphere are well correlated with solar wind variations. In one event, a significant (20–30%) decrease of the midlatitude ionospheric electron density over a large latitudinal range was related to a sudden drop in the solar wind pressure and a northward turning of the interplanetary magnetic field, and the density decrease became larger at lower latitudes. In another event, periodic perturbations in the dayside equatorial ionospheric E × B drift and electrojet were closely associated with variations in the interplanetary electric field. Since the solar wind is always changing with time, it can be a very important and common source of ionospheric perturbations at middle- and low-latitudes. The relationship between solar wind variations and significant ionospheric perturbations has important applications in space weather.     

Comparison of heliospheric conditions near the earth during four recent solar maxima

Statistical properties of the daily averaged values of the solar activity (sunspot numbers, total solar irradiance and 10.7 cm radio emission indices), the solar wind plasma and the interplanetary magnetic field parameters near the Earth’s orbit are investigated for a period from 1964 to 2002 covering the maxima of four solar cycles from 20th to 23rd. Running half-year averages show significant solar cycle variations in the solar activity indices but only marginal and insignificant changes in comparison with background fluctuations for heliospheric bulk plasma and magnetic field parameters. The current 23rd cycle maximum is weaker than 21st and 22nd maxima, but slightly stronger than 20th cycle in most of solar and heliospheric manifestations.     

嵌套闭磁场结构内CME产生和传播的数值模拟

给出了特殊类型的日冕物质抛射(CME)数值模拟定性结果,这种CME核心闭磁场结构前半部分磁力线的方向与太阳整体偶极场磁力线的方向相反．计算结果表明,这种CME核心闭磁场结构磁力线与太阳整体偶极场反向磁力线之间存在过渡磁场结构,在向外传播时过渡磁场结构所占的面积逐渐增大．这一结果可以用来解释飞船为什么能够观测到一类双极磁云,这类磁云前半部分磁场方向与太阳整体偶极场方向相反．为了模拟这一数值结果,强调需要采用包含嵌套闭磁场的冕流背景结构,并在合适的位置触发CME．