Modeling the (upper) solar atmosphere including the magnetic field

The atmosphere of the Sun is highly structured and dynamic in nature. From the photosphere and chromosphere into the transition region and the corona plasma-β changes from above to below one, i.e., while in the lower atmosphere the energy density of the plasma dominates, in the upper atmosphere the magnetic field plays the governing role – one might speak of a “magnetic transition”. Therefore the dynamics of the overshooting convection in the photosphere, the granulation, is shuffling the magnetic field around in the photosphere. This leads not only to a (re-)structuring of the magnetic field in the upper atmosphere, but induces also the dynamic reaction of the coronal plasma, e.g., due to reconnection events. Therefore the (complex) structure and the interaction of various magnetic patches is crucial to understand the structure, dynamics and heating of coronal plasma as well as its acceleration into the solar wind.

The present article will emphasize the need for three-dimensional modeling accounting for the complexity of the solar atmosphere to understand these processes. Some advances on 3D modeling of the upper solar atmosphere in magnetically closed as well as open regions will be presented together with diagnostic tools to compare these models to observations. This highlights the recent success of these models which in many respects closely match the observations.     

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).     

Local-area helioseismology as a diagnostic tool for solar variability

Dynamical and thermal variations of the internal structure of the Sun can affect the energy flow and result in variations in irradiance at the surface. Studying variations in the interior is crucial for understanding the mechanisms of the irradiance variations. “Global” helioseismology based on analysis of normal mode frequencies, has helped to reveal radial and latitudinal variations of the solar structure and dynamics associated with the solar cycle in the deep interior. A new technique, - “local-area” helioseismology or heliotomography, offers additional potentially important diagnostics by providing three-dimensional maps of the sound speed and flows in the upper convection zone. These diagnostics are based on inversion of travel times of acoustic waves which propagate between different points on the solar surface through the interior. The most significant variations in the thermodynamic structure found by this method are associated with sunspots and complexes of solar activity. The inversion results provide evidence for areas of higher sound speed beneath sunspot regions located at depths of 4–20 Mm, which may be due to accumulated heat or magnetic field concentrations. However, the physics of these structures is not yet understood. Heliotomography also provides information about large-scale stable longitudinal structures in the solar interior, which can be used in irradiance models. This new diagnostic tool for solar variability is currently under development. It will require both a substantial theoretical and modeling effort and high-resolution data to develop new capabilities for understanding mechanisms of solar variability.     

Acceleration and transport of energeticparticles at CME-driven shocks

Historically, solar energetic particle (SEP) events are classified in two classes as “impulsive” and “gradual”. Whether there is a clear distinction between the two classes is still a matter of debate, but it is now commonly accepted that in large “gradual” SEP events, Fermi acceleration, also known as diffusive shock acceleration, is the underlying acceleration mechanism. At shock waves driven by coronal mass ejections (CMEs), particles are accelerated diffusively at the shock and often reach > MeV energies (and perhaps up to GeV energies). As a CME-driven shock propagates, expands and weakens, the accelerated particles can escape ahead of the shock into the interplanetary medium. These escaping energized particles then propagate along the interplanetary magnetic field, experiencing only weak scattering from fluctuations in the interplanetary magnetic field (IMF). In this paper, we use a Monte-Carlo approach to study the transport of energetic particles escaping from a CME-driven shock. We present particle spectra observed at 1 AU. We also discuss the particle “crossing number” at 1AU and its implication to particle anisotropy. Based on previous models of particle acceleration at CME-driven shocks, our simulation allows us to investigate various characteristics of energetic particles arriving at various distances from the sun. This provides us an excellent basis for understanding the observations of high-energy particles made at 1 AU by ACE and WIND.     

Solar wind dependence of the Pc1 wave activity

There are a host of factors influencing the excitation of Pc1 geomagnetic pulsations, which are ULF waves in the frequency range between 0.2 and 5 Hz. We have studied carefully the dependence of the pearl-type Pc1 activity at Sodankylä, Finland (L = 5.1) on the plasma density N in front of the magnetosphere, the bulk velocity V of the solar wind, and the intensity B of the IMF. The result is as follows: high values of N and reduced values of V are favorable to appearance of Pc1, whereas the dependence of Pc1 activity on B is practically absent. We also show that the probability of Pc1 occurrence decreases with the interplanetary electric field, and increases with solar wind impact pressure and with the plasma to magnetic pressure ratio “beta”.     

Investigation of magnetospheric processes with the use of a source of strong magnetic field in the ionosphere

More than 20 years ago V.P. Shabansky suggested that the magnetic system installed aboard the satellite, could be used as a physical instrument for studying the processes which occur in the near Earth space. The corresponding space scales of an artificial “magnetosphere”—“magnisphere”—are 10 m in the experiment with relatively small magnets in the ionosphere and 100 m in the solar wind. The corresponding similarity criteria are estimated. The possible scheme of the experiment with a superconducting magnet (magnetic moment 10⁵ A · m²) installed aboard the satellite is considered. The experimental complex includes a number of systems for measuring the fluxes of charged particles in a wide energy range, DC electric and magnetic fields, the electromagnetic fields in different frequency bands (from X-rays to radio). The scientific objectives are discussed in detail.     

Total solar and spectral irradiance variations from solar cycles 21 to 23

Total solar and UV irradiances have been measured from various space platforms for more than two decades. More recently, observations of the “Variability of solar IRradiance and Gravity Oscillations” (VIRGO) experiment on SOHO provided information about spectral irradiance variations in the near-UV at 402 nm, visible at 500 nm, and near-IR at 862 nm. Analyses based on these space-borne irradiance measurements have convinced the skeptics that solar irradiance at various wavelengths and in the entire spectrum is changing with the waxing and waning solar activity. The main goal of this paper is to review the short- and long-term variations in total solar and spectral irradiances and their relation to the evolution of magnetic fields from solar cycles 21 to 23.     

Ion signatures of reconnection at the magnetopause

For a rotational MHD discontinuity the bulk flow is Alfvénic in the de Hoffmann-Teller frame. Using AMPTE/IRM data, we present ion distribution functions during three crossings of the dayside low-latitude magnetopause. For these crossings a well defined de Hoffmann-Teller frame can be found, but the field-aligned bulk flow is always slower than the Alfvén speed. Nevertheless, we find signatures in the distribution functions that provide clear evidence for magnetic reconnection: solar wind ions reflected off the magnetopause, “D-shaped” solar wind ion distributions in the boundary layer, and counterstreaming of solar wind ions and ionospheric ions in the boundary layer.     

膨胀磁场和太阳风粒子相互作用三维数值模拟

采用三维模型,使用混合网格质点法HPIC(Hybrid Particle-in-Cell)对膨胀的磁场和太阳风相互作用过程进行数值模拟.研究了线圈产生的偶极子磁场在注入等离子体后和太阳风粒子的相互作用过程,并对以不同速度入射的等离子体引起的太阳风粒子的变化和磁场变化进行了比较.研究结果表明,偶极子磁场和太阳风作用时会产生弓形激波,此时磁压等于太阳风粒子的动压,当向线圈产生的偶极子磁场中注入高能等离子体时引起磁场膨胀,膨胀的磁场将会排斥太阳风粒子向外运动,从而引起弓形激波的变化,增大与太阳风相互作用的面积,并且粒子入射速度越大,磁场膨胀越明显,与太阳风相互作用愈强.      

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

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

Synoptic magnetic field in cycle 23 and in the beginning of the cycle 24

The SOHO/MDI data provide the uniform time series of the synoptic magnetic maps which cover the period of the cycle 23 and the beginning of the cycle 24. It is very interesting period because of the long and deep solar minimum between the cycles 23 and 24. Synoptic structure of the solar magnetic field shows variability during solar cycles. It is known that the magnetic activity contributes to the solar irradiance. The axisymmetrical distribution of the magnetic flux (Fig. 3c) is closely associated with the ‘butterfly’ diagram in the EUV emission (Benevolenskaya et al., 2001). And, also, the magnetic field (B_∥) shows the non-uniform distributions of the solar activity with longitude, so-called ‘active zones’, and ‘coronal holes’ in the mid-latitude. Polar coronal holes are forming after the solar maxima and they persist during the solar minima. SOHO/EIT data in the emission of Fe XII (195 Å) could be a proxy for the coronal holes tracking. The active longitudinal zones or active longitude exist due to the reappearance of the activity and it is clearly seen in the synoptic structure of the solar cycle. On the descending branch of the solar cycle 23 active zones are less pronounced comparing with previous cycles 20, 21 and 22. Moreover, the weak polar magnetic field precedes the long and deep solar minimum. In this paper we have discussed the development of solar cycles 23 and 24 in details.     

The role of the solar magnetic field systems in modulating the solar irradiance

Variations of indices that characterize various systems of the large-scale solar magnetic field (LSSMF) - magnetic field multipoles of different order, LSSMF energy index, index of the effective solar multipole, etc.- are compared with variations of the solar irradiance in different frequency ranges during 1978–1996. The role of the local and global magnetic fields in modulating the solar irradiance is investigated in various time intervals, in particular, in different phases of the 11-year solar cycle.     

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

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

The effects of the solar magnetic polarity and the solar wind velocity on Bz-component of the interplanetary magnetic field

We have studied the effect of both solar magnetic polarity and the solar wind velocity on the Bz-component of the interplanetary magnetic field, IMFBz, for the minimum activity of the solar cycles 21, 22, 23 and 24. We made a statistical study of IMFBz in the first section which is considered as an extension of Lyatsky et al. (2003). They made a statistical study of IMFBz for two periods of minimum solar activity 22 and 23 related to 1985–1987 and 1995–1997 when the solar magnetic field had opposite polarity. Our results seem to be consistent with the results obtained by Lyatsky et al. (2003). We found that there is a dependence of IMFBz on the IMFBx and the solar magnetic polarity for the minimum periods of the selected four solar cycles. In addition, we found that there is a dependence of IMFBz on the solar wind velocity.     

Solar wind interaction with lunar crustal magnetic anomalies

Using Lunar Prospector data, we review the magnetic field and electron signatures of solar wind interaction with lunar crustal magnetic sources. Magnetic field amplifications, too large to represent direct measurements of crustal fields, appear in the solar wind over strong crustal sources, with the chance of observing these amplifications depending on upstream solar wind parameters. We often observe increases in low-energy (?100 eV) electron energy fluxes simultaneously with large magnetic field amplifications, consistent with an increase in plasma density across a shock surface. We also often observe low frequency wave activity in the magnetic field data (both broadband turbulence and monochromatic waves), often associated with electron energization, sometimes up to keV energies. Electron energization appears to be correlated more closely with wave activity than with magnetic amplifications. Detailed studies of the interaction region will be necessary in order to understand the physics of the Moon–solar wind interaction. At present, the Moon represents the only natural laboratory available to us to study solar wind interaction with small-scale crustal magnetic fields, though simulation results and theoretical work can also help us understand the physical processes at work.     

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

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

太阳大气中有限振幅磁场扰动的动力学过程

本文从完整的磁流体动力学方程组出发,通过太阳大气中磁力线管根部有限振幅磁场的扰动,研究了非线性磁场的动力学演化。假设初始磁场位形足β<<1的势场,根部磁力线管磁场扰动,驱动等离子体运动,一部分磁能转换为等离子体动能。等离子体压缩运动具有快磁声波的特性。计算结果给出非线性磁场演化的定量关系,可以解释太阳大气中日冕活动过程。也可用于模拟实验室里高β实验装置中的等离子体的持性。      

The local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

How much of the solar irradiance variations is caused by the magnetic field at the solar surface?

The contribution to total solar irradiance variations by the magnetic field at the solar surface is estimated. Detailed models of the irradiance changes on the basis of magnetograms show that magnetic features at the solar surface account for over 90% of the irradiance variations on a solar rotation time scale and at least 70% on a solar cycle time scale. If the correction to the VIRGO record proposed by Fröhlich & Finsterle (2001) is accepted, then magnetic features at the solar surface are responsible for over 90% of the solar cycle irradiance variations as well.     

月球尾迹的二维MHD模拟

采用理想的二维单流体MHD方程,对太阳风通过月球时所形成的尾迹结构进行数值模拟,得到了太阳风尾迹的粒子分布及磁场分布.模拟结果表明,在月球背阳面的本影区,太阳风粒子密度急剧下降,行星际磁场增强.当行星际磁场与太阳风流动方向平行时,尾迹被拖得很长,而磁场与太阳风流动方向垂直时,尾迹较短.