行星际磁场的螺旋扇形过渡区

从行星际磁场的大尺度螺旋形构型和扇形边界附近太阳风流动与冕旒的可能相关,本文得到一个推论,即在行星际空间可能存在一种较厚的螺旋扇形过渡区。行星际磁场和太阳风的实地观测证实了这种较厚的螺旋扇形过渡区的存在。在所分析的45个螺旋扇形过渡区中,磁场强度都不为零;大部分大于或小于周围平均场强。本文进一步分析了磁增大和磁减小两类过渡区中的物理性质和可能成因。      

二维磁流体动力学平衡日冕(Ⅱ)

基于解析和数值相结合的方法,进一步讨论了非均匀引力场中日冕的二维磁流体动力学平衡。对临界点进行了比较仔细的处理。得到了包含闭场区、中性片和开场区的大尺度日冕磁场位形,闭场区和中性片构成冕流结构。在高纬和低纬地区几个太阳半径之外,等离子体径向流动速度超过了局地声速和局地Alfvén速度。在1AU处,太阳风速度可达到400kms_-1以上      

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.     

日冕三维太阳风模拟

COIN-TVD MHD模型是近年发展起来的能有效实现日冕–行星际三维太阳风模拟的模型.本文利用此模型针对日冕区三维太阳风进行研究,为了模拟日冕太阳风的加热加速,对模型中的体积加热项做了调整.在磁流体模拟中,减小磁场散度的误差是关键问题之一,在调整体积加热项后应用扩散法、八波法、扩散八波法,对2199卡林顿周的背景太阳...     

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

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

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

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

部分开放多极磁场中的日冕磁绳灾变

采用球坐标下二维三分量理想MHD模型，研究部分开放多极背景磁场中日冕磁绳的灾变现象．背景磁场由含3个闭合双极场的冕流和带赤道中性电流片的开放场构成，磁绳位于中心双极场的下方，其特性由环向磁通和轴向磁通表征．对给定的环向磁通，存在轴向磁通的一个临界值；对给定的轴向磁通，也存在环向磁通的一个临界值．在该临界值以下，磁绳附着于太阳表面，系统处于平衡状态；该临界值一旦被超越，磁绳将脱离太阳表面向上喷发，说明部分开放多极背景磁场中的日冕磁绳系统存在灾变现象．本文算例表明，灾变点对应的磁能阈值超过对应部分开放场(中心双极场开放，两侧的双极场仍维持闭合)能量约15％，其超过部分可为日冕物质抛射一类太阳爆发提供能源．     

嵌套闭磁场日冕背景结构对触发CME特征的影响

在子午面内,研究具有嵌套闭磁场结构冕流背景对触发日冕物质抛射(CME)特征的影响.在冕流背景磁场结构内有三个小尺度的闭磁场结构,其中中间小尺度闭磁场结构的磁场方向和冕流整体偶极磁场方向相反.CME触发模型在这个小尺度闭磁场结构下方浮出,它具有同心圆形磁场结构,半径为a=0.1Rs(Rs为太阳半径);CME触发模型前半部分磁场方向和触发位置处小尺度闭磁场结构的磁场方向相反,与冕流整体偶极磁场方向相同.数值模拟结果表明,当CME触发模型中心等离子压力与边界压力之比m≥2时,上浮CME触发模型可以触发CME;当m＜2时,上浮CME触发模型不能触发CME,计算结论的误差小于1%.      

一个基于观测的太阳风背景模型

建立由太阳光球磁场和日冕偏振亮度等观测约束的单流体太阳风模型,包括日冕和太阳风的等离子体密度、速度和磁场,温度还有待于以后处理.这里采用高山观测台(HAO)MKⅢ的日冕偏振亮度(pB)在1.36Rs上的观测概图,根据Guhathakurta在1996年发展的日冕电子密度反演模型确定日冕的电子密度分布.同时采用Wilcox太阳观测台(WSO)的光球磁场视向分量的观测概图作为底部边界,根据Zhao等在1994年发展的水平电流-电流片(HCCS)模型得到全球磁场.Phillips在1995年及McComas在2003年分别用Ulysses第一次和第二次跨极飞行的观测发现,归一化到1 AU的太阳风动量流密度除了在10°～30°的纬度范围内略低以外几乎不变.根据这一结论,结合已经得到的密度数据,就可以得到日冕和太阳风的速度.将上面的模型应用于1918卡林顿自转周稳态太阳风的研究,结果与太阳活动极小期的观测基本相符,但是与观测相比较低速高密度区偏大,因此密度模型还有待改进.      

太阳宁静区磁漏斗结构的扩展特征

通过对观测的光球纵向磁场进行势场外推, 重构了一个太阳宁静区光球以上的磁场结构. 结果显示, 在20 Mm 以下, 开放磁力线呈现多个明显的小磁漏斗结构, 这些小磁漏斗结构随高度扩展, 并且在20 Mm 左右融合成大的漏斗结构. 通过系统地研究磁漏斗结构横截面积随高度的变化趋势, 发现太阳宁静区磁漏斗结构的截面积随高度近似线性扩展, 磁漏斗结构在较低高度上(<20Mm) 扩展的速度比在较高高度上 (>20Mm) 扩展的速度要快. 这一结果对太阳风起源和磁环中物质流动的二维数值模拟具有重要的意义. 同时还发现, 闭合磁力线的数目随高度以指数函数的形式减少.      

Ray Structure of the Coronal Streamer Belt and Its Manifestation as Sharp Large Peaks of Solar Wind Plasma Density at the Earth''''s Orbit

The white-light corona calibrated data with processing level L1 from the LASCO-C2/SOHO instrument, and data from the Wind spacecraft with one-hour and one-minute time resolution on quasi-stationary slow (v between 300-450 km/s at the Earth's orbit) the Solar Wind (SW) parameters in the absence of sporadic SW streams are examined. Within distances from the Sun's center less than R in the range of 20-30 Rs, (Rs, the solar radius), slow wind is known as the streamer belt, and at larger distances it is called the He-liospheric Plasma Sheet (HPS). It is shown that the streamer belt comprises a sequence of pairs of rays. In general, ray brightnesses in each pair can differ, and the magnetic field is oppositely directed in them. The neutral line of the radial magnetic field of the Sun runs along the belt between the rays of each of the pairs. The area in which the streamer belt intersects the ecliptic plane and which lies at the central meridian, will be recorded at the earth's orbit with a time delay of 5-6 days, in the form of one or several peaks with Nmax> 10cm-3. Furthermore, the simplest density profile of the portion of the HCS has the form of two peaks of a different or identical amplitude . The such a profile is observed in cases where the angle of intersection of the streamer belt with the ecliptic plane near the Sun is sufficiently large, i.e. close to 90°. The two-ray structure of the cross-section of the streamer-belt moves from the Sun to the Earth, it retains not only the angular size of the peaks but also the relative density variations, and the position of the neutral line (sector boundary) in between. At the Earth's orbit the ray structure of the streamer belt provides the source for sharp (i.e. with steep fronts of a duration of a few minutes or shorter) solar wind plasma density peaks (of a duration of several hours) with maximum values Nmax> 10cm-3.     

Numerical MHD modeling of the solar driver for the destabilization of a coronal helmet streamer

We investigate the forms of the solar driver which cause the destabilization of helmet streamers. Two forms of solar drivers are considered; (i) emergence of a flux-rope from sub-photospheric levels and (ii) application of a photospheric shear motion to a streamer-flux rope system. Numerical results showed that both cases exhibit the characteristics of commonly observed coronal mass ejections (CMEs), but the propagation speed of the CME is higher than the background solar wind speed when the solar driver is the emerging magnetic flux and is the same as the solar wind speed when the photospheric shear is used as the solar driver. The energy constraint allowing the magnetic field transition from a closed to an open configuration is also addressed.     

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

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.     

冕洞特征参数与重现型地磁暴关系的统计研究

在提取冕洞特征参数的基础上,利用1996年到2005年8月近十年来对地磁扰动有影响的356个冕洞事例,定量分析了冕洞特征参数(包括冕洞的面积比、经纬度跨度等)与冕洞高速流特征、重现型地磁扰动特征(包括扰动大小和持续时间等)之间的相关性,研究发现,从引起地磁扰动的冕洞在整个太阳活动周的分布来看,在地磁扰动峰年中冕洞影响同样具有重要的贡献;冕洞高速流太阳风速度与地磁扰动强度之间存在较强的相关性,而高速流中太阳风速度与冕洞面积比关系不大,与冕洞亮度存在一定相关性;冕洞的经度跨度与地磁扰动持续时间存在很强的正相关性.      

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.     

太阳活动下降期的日冕大尺度磁场和冕洞分布

本文应用十八、十九周下降期的地磁观测资料导出了对应时期太阳风速度的卡林顿经度变化,并结合由极盖区地磁观测导出的行星际磁扇形结构确定十八、十九周下降期太阳南北磁极极性、类偶极点的太阳余纬和卡林顿经度,从而确定了该时期的冕洞分布,并与二十周冕洞分布进行比较.      

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

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

磁场强度对日冕定态结构的影响

以二维ＭＨＤ模型及时变方法为基础,内外边界完整的设影特征线边界条件,考察了太阳日冕大气的定态结构随偶极场强度的变化情况。模拟结果表明：随着偶极场强度的增加,磁场对太阳风的约束增强,低纬闭磁场打开程度减少,高纬与低纬区速度差增加,并且在阿尔文马赫数为１的点附近达到最大,速度过度区变陡;随着日心距离增加,低纬区宽度减小,速度过渡区变陡,可定性解释Ｕｌｙｓｓｅｓ飞船的新观测事实。