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.     

The Compton–Getting effect of energetic particles with an anisotropic pitch-angle distribution

This paper presents a simulation of anisotropy measurements by the low-energy charged particle (LECP) experiment on Voyager 1 for cases when the particle pitch-angle distribution function in the solar wind plasma reference frame is not isotropic. The model includes both the Compton–Getting anisotropy and perpendicular diffusion anisotropy that possibly exists in the upstream region of the termination shock. The results show that the Voyager 1 data cannot rule out either the model with zero solar wind speed or the one with a finite speed on qualitative basis. The determination of solar wind speed using the Compton–Getting effect is affected by the assumption of the magnetic field direction and perpendicular diffusion anisotropy. Because the pitch-angle distribution anisotropy is so large, a small uncertainty in the magnetic field direction can produce very different solar wind speeds ranging from ∼0 to >400 km/s. In fact, if the magnetic field is chosen to be in the Parker spiral direction, which is consistent with the magnetometer measurement on Voyager 1, the derived solar wind speed is still close to the supersonic value. Only the two lowest-energy channels of the LECP instrument may give a definitive answer to the solar wind speed. However, because these channels contain a very high level of cosmic ray background, an uncertainty of just a few percent in the background can entirely hamper the estimate of solar wind speed.     

Cosmic ray modulation in the outer heliosphere: Predictions for cosmic ray intensities up to the heliopause along Voyager 1 and 2 trajectories

Time dependent cosmic ray modulation in the outer heliosphere is calculated and results are compared to Voyager 1 and 2 observations using a two-dimensional time-dependent cosmic ray transport model. We predict possible future 133–242 MeV proton observations along the Voyager 1 and 2 spacecraft trajectories. Recent theoretical advances in cosmic ray transport parameters are introduced in order to provide a time-dependence for the assumed transport parameters used in the model. This leads to results that are in general compatible with the spacecraft observations in the inner and outer heliosphere over multiple solar cycles. However, for the outer heliosphere, we find that the Voyager 1 and 2 spacecraft observations cannot be fitted with an identical set of parameters along both trajectories. This indicates a possible asymmetric heliosphere or a symmetric heliosphere but with different diffusion parameters in the northern and southern hemispheres, respectively. Furthermore, results indicate that Voyager 2 observations are still under the influence of solar cycle related changes because of the large modulation volume between the heliopause and spacecraft location in contrast to Voyager 1 which shows a steady increase in cosmic ray intensities.     

Modeling of the outer heliosphere with the realistic solar cycle

Time-dependent kinetic-continuum model of the solar wind interaction with the two-component local interstellar cloud (LIC) has been developed recently [Izmodenov, V., Malama, Y.G., Ruderman, M.S. Solar cycle influence on the interaction of the solar wind with local interstellar cloud. Astron. Astrophys. 429, 1069–1080, 2005a.]. Here, we adopted this model to the realistic solar cycle, when the solar wind parameters at the Earth’s orbit are taken from space data. This paper focuses on the results related to the termination shock (TS) excursion with the solar cycle that may help to understand Voyager 1 data obtained at and after the crossing of the termination shock and to predict the time of the TS crossing by Voyager 2.     

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

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

Saturn radio emission and the solar wind: Voyager-2 studies

Voyager 2 data from the Plasma Science experiment, the Magnetometer experiment and the Planetary Radio Astronomy experiment were used to analyze the relationship between parameters of the solar wind/interplanetary medium and the nonthermal Saturn radiation. Solar wind and interplanetary magnetic field properties were combined to form quantities known to be important in controlling terrestrial magnetospheric processes.The Voyager 2 data set used in this investigation consists of 237 days of Saturn preencounter measurements. However, due to the immersion of Saturn and the Voyager 2 spacecraft into the extended Jupiter magnetic tail, substantial periods of the time series were lacking solar wind data. To cope with this problem a superposed epoch method (CHREE analysis) was used. The results indicate the superiority of the quantities containing the solar wind density in stimulating the radio emission of Saturn — a result found earlier using Voyager 1 data — and the minor importance of quantities incorporating the interplanetary magnetic field.     

Saturn kilometric radiation as a monitor for the solar wind?

Since the Voyager mission it is known that Saturn Kilometric Radiation (SKR) is strongly influenced by external forces, i.e., the solar wind and in particular the solar wind ram pressure. Recent studies using Cassini data essentially confirmed these findings for particular periods during the first Cassini orbit of Saturn. The data coverage of SKR by the Cassini/RPWS experiment for the period of six months prior to Saturn Orbit Insertion (July 1, 2004) is rather continuous, whereas there are gaps in the solar wind plasma data. The strong correlation of SKR with the solar wind may provide an indication on the variations of the solar wind plasma, specifically during the gap periods. These periods lacking solar wind data are substituted by Ulysses solar wind data which have been propagated over ∼4 AU, applying magnetohydrodynamic propagation models. Cross correlation studies showed that Ulysses solar wind data can be taken as a substitute for missing Cassini data. The use of SKR as monitor for solar wind variations is discussed. With the present set of observations the SKR proxy lacks significant reliability.     

The heliospheric modulation of 3–10 MeV electrons: Modeling of changes in the solar wind speed in relation to perpendicular polar diffusion

The discrepancy between cosmic ray model predictions representing solar minimum conditions in the heliosphere and the 3–10 MeV post-1998 electrons observations by the Kiel Electron Telescope (KET) onboard Ulysses suggests the need for consistent changes in model parameters with increasing solar activity. In order to reduce this discrepancy, an effort is made to model the KET observations realistically during periods of increased solar activity by applying an advanced three-dimensional, steady-state electron modulation model based on Parker’s transport equation including the Jovian electron source. Some elements of the diffusion tensor which were not previously emphasized are revisited. A new relation is also established between the latitudinal dependence of the solar wind speed and the perpendicular polar diffusion. Based on this relation, a transition of an average solar wind speed from solar minimum to solar maximum conditions, as observed on board the Ulysses spacecraft, is modeled on the concept of the time-evolution of large polar coronal holes. These changes are correlated to different scenarios of the enhancement of perpendicular polar diffusion. Effects of these scenarios are illustrated, as a series of steady-state solutions, on the computed 7 MeV Jovian and galactic electrons in comparison with 3–10 MeV electrons observed from the period 1998 to the end of 2003. It is shown that this approach improves compatibility with the KET observations but it also points to the need for a time-dependent electron modulation model to fully describe modulation during moderate to extreme solar maximum conditions.     

Dynamics of solar cosmic ray events: Processes at large heliocentric distances (⪢1 AU)

Observations of solar cosmic ray events far from the sun (?1 AU) became possible after the launch of Pioneer 10 in 1972. Four spacecraft have now travelled beyond the orbit of Jupiter - Pioneer 10/11 and Voyager 1/2 — and are producing a growing body of distant observations of solar cosmic ray events. Initial studies using Pioneer 10/11 data out to ～6 AU interpreted flare particle observations in terms of a diffusion model, including the effects of convection and adiabatic energy loss. This model enjoyed general success in explaining the time-intensity profiles in cases where the spacecraft connection longitude at the sun did not change significantly with time. The results implied that the radial diffusion coefficient (K_r) increased slowly with distance over that radial range. More recent results at larger distances imply that K_r may begin to decrease beyond ～5 AU. It is not yet clear whether the standard diffusion model will be adequate to explain solar events well beyond 5 AU. The fact that large events at very large distances can last up to two solar rotations implies that solar wind stream structure will also play a role in the event dynamics. In general, however, observations at large distances offer perhaps the best hope of separating interplanetary propagation effects from coronal storage and propagation effects which frequently dominate observed event profiles at 1 AU.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

GEO轨道相对论电子日积分通量预报统计建模

基于辐射带相对论电子哨声波局地加速理论,将地磁AE指数作为源电子通量和通量各向异性的指标,将地磁Dst指数作为损失机制的指标,利用滑动窗口线性滤波器方法,建立了一个地球静止轨道大于2MeV相对论电子预报模型.利用该模型开展了2000-2009年地球静止轨道相对论电子通量预报试验.研究发现,这10年总预报效率为0.818,2003年的预报效率（0.633）最低,2009年的预报效率（0.856）最高.模型预报效果与持续模型相比有很大提高,略低于利用太阳风参数作为输入的同类预报模型的预报效果.这说明即使在缺少太阳风参数的情况下,该模型利用地磁扰动参数也能取得较好的预报效果.当模型输入参数增加了太阳风速度时,即综合考虑了行星际扰动和磁层扰动对辐射带粒子加速过程的影响,模型逐年的预报效率进一步提升.其中,2005年的预报效率提升了9.5%,这10年的总预报效率增加到0.848,预报值与实测值之间的线性相关系数为0.918,均方根误差为0.422.      

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.     

Modelling of galactic Carbon in an asymmetrical heliosphere: Effects of asymmetrical modulation conditions

Observations of galactic cosmic rays (GCRs) from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that in addition to a possible global asymmetry in the north–south dimensions (meridional plane) of the heliosphere, it is also possible that different modulation (turbulence) conditions could exist between the two hemispheres of the heliosphere. We focus on illustrating the effects on GCR Carbon of asymmetrical modulation conditions combined with a heliosheath thickness that has a significant dependence on heliolatitude. To reflect different modulation conditions between the two heliospheric hemispheres in our numerical model, the enhancement of both polar and radial perpendicular diffusion off the ecliptic plane is assumed to differ from heliographic pole to pole. The computed radial GCR intensities at polar angles of 55° (approximating the Voyager 1 direction) and 125° (approximating the Voyager 2 direction) are compared at different energies and for both particle drift cycles. This is done in the context of illustrating how different values of the enhancement of both polar and radial perpendicular diffusion between the two hemispheres contribute to causing differences in radial intensities during solar minimum and moderate maximum conditions. We find that in the A > 0 cycle these differences between 55° and 125° change both quantitatively and qualitatively for the assumed asymmetrical modulation condition as reflected by polar diffusion, while in the A < 0 cycle, minute quantitative differences are obtained. However, when both polar and radial perpendicular diffusion have significant latitude dependences, major differences in radial intensities between the two polar angles are obtained in both polarity cycles. Furthermore, significant differences in radial intensity gradients obtained in the heliosheath at lower energies may suggest that the solar wind turbulence at and beyond the solar wind termination shock must have a larger latitudinal dependence.     

Rapid and Asymmetric Response of the Earth's Bow Shock:Multipoint Observations

We analyze observations of three bow shock crossings which occurred during 2007, using upstream data from STEREO A/B, ACE and WIND, combined with multi-point THEMIS and Cluster data, and TC-1 data located near noon. During the crossing of 7 May 2007, we find that following a rapid reduction in solar wind ram pressure and subsequent pressure pulse seen by ACE and WIND upstream, the bow shock responds asymmetrically from dawn to dusk. Cluster data on the dawn-side suggest the bow shock is significantly flared and responds rapidly to the pulse arrival, while TC-1 at noon, and THEMIS on the dusk-side, are well matched to the model bow shock, but show a delayed response. The crossings observed on 21 May and 2 June show contrasting response matching the model boundary for northward Interplanetary Magnetic Field (IMF). The IMF and solar wind plasma data suggest that the bow shock crossing at dawn-dusk side and subsolar point were mainly caused by large and smaller scale features of the solar wind ram pressure rise rather than the influence of IMF.      

In situ observations of magnetic field fluctuations

It is well known that the irregularities of the magnetic field are intimately related to the motion of charged particles. Although transport theories need the spatial and time variations of the magnetic field as input, in situ observations are very limited. Ulysses observations have provided a major step forward by entering the unexplored high latitude regions of the heliosphere, the knowledge of which is vital to interpret particle flux measurements, even at the ecliptic. We analyze the magnetic field data of Ulysses during the mission to study the waves and discontinuities in the heliosphere at different locations, covering a total sunspot cycle. Various tools are employed, including power spectral and structure function analysis. A remarkable difference was found between the fluctuations in the fast and slow solar wind. We argue that the latitudinal extent of the high speed solar wind contributes significantly to the latitudinal variation of the transport parameters, which should also affect the 11 (and 22) year modulation cycle.     

Modeling and experimental study of the 27-day variation of galactic cosmic-ray intensity for a solar wind velocity depending on heliolongitude

We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic-ray (GCR) intensity with a spatial variation of the solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving the corresponding Maxwell equations with a variable solar wind speed, which reproduces in situ observed experimental data for the time interval to be analyzed (24 August 2007–28 February 2008). We perform model calculations for the GCR intensity using the variable solar wind and the corresponding magnetic field. Results are compatible with experimental data; the correlation coefficient between our model predictions and observed 27-day GCR variation is 0.80 ± 0.05.     

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

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

异常SSC事件的行星际原因分析

磁暴急始(SSC)是强烈太阳风动压或行星际激波与磁层相互作用的结果.通常SSC事件的上升时间在4～10 min,我们把上升时间超过15 min的SSC事件称为异常SSC事件.本文利用地磁SYM-H指数鉴别出了5个有地磁观测历史以来发生的上升时间大于15 min的异常SSC事件,并利用Wind,ACE,IMP 8,Goes,Geotail多点卫星太阳风观测数据和地磁观测数据,分析了异常SSC事件的行星际原因.结果表明,异常SSC事件通常都是强烈行星际扰动引起的,5个异常SSC事件有4个对应于行星际激波,有3个对应于多步太阳风动压跃变,有1个对应于行星际电场大幅度变化;由行星际激波产生的异常SSC事件,其上升时间依赖于行星际激波的方向,方向相对于日地连线越偏,上升时间越长;异常SSC事件上升时间与行星际磁场方向关系不明显.      

Solar wind-driven electron radiation belt response functions at 100-min time scales

We present a simple yet numerically robust technique, using autoregressive linear filters, to remove unwanted “colored noise” from solar wind and radiation belt electron data at sub-daily resolution. The remaining signal is then studied using finite impulse response linear prediction filters to represent the driven portion of the linear dynamics that describe the coupling between solar wind speed and electron flux. Sub-daily resolution response profiles covering magnetic L-shells between 1.1 and 8.0 R_E are presented which are consistent with daily resolution response functions. Namely, while there is strong global coherence governing electron flux dynamics, there are at least two distinct responses. The first response is an immediate dropout of electrons between L = 4 and L = 7 that is at least a partly adiabatic effect associated with enhancements in the ring current. This is followed by a 1–2 day delayed enhancement across the same L-shells that is likely a result of increased radial diffusion. The second response is an immediate enhancement seen between L = 3 and L = 4 with a typical duration of less than one day. Plausible explanations for this second response are briefly discussed, but neither empirical nor theoretical evidence can establish conclusively a definite physical cause. Finally, the response profiles show significant solar cycle and seasonal dependencies, indicating that better model output might be achieved with: (1) additional simultaneous solar wind inputs; (2) more sophisticated dynamical model structures capable of incorporating non-linear feedback; and/or (3) time-adaptive linear filters that can track non-stationary dynamics in time.