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.     

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.     

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.     

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.     

Possible physical properties of the flow in the vicinity of the heliopause

An interface between the fully ionized hydrogen plasma of the solar wind (SW) and the partially ionized hydrogen gas flow of the local interstellar medium (LISM) is formed as a region where there is a strong interaction between these two flows. The interface is bounded by the solar wind termination shock (TS) and the LISM bow shock (BS) and is separated on two regions by the heliopause (HP) separating the solar wind and charged component of the LISM (plasma component below). The BS is formed due to the deceleration of the supersonic LISM flow relative to the solar system. Regions of the interface between the TS and HP and between the HP and BS were in literature named as the inner and outer heliosheaths, respectively. An investigation of the structure and physical properties of the heliosheath is at present especially interested due to the fact that Voyager-1 and Voyager-2 have crossed the TS in December 2004 (Burlaga, L.F., Ness, N.F., Acuna, M.Y., et al. Crossing the termination shock into the the heliosheath. Magnetic fields. Science 309, 2027–2029, 2005; Fisk, L.A. Journey into the unknown beyond. Science 309, 2016–2017, 2005; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Voyager 1 in the foreshock, termination shock and heliosheath. Science 309, 2020–2024, 2005; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. Voyager 1 explores the termination shock region and the heliosheath beyond. Science 309, 2017–2020, 2005) and in September 2007 (Jokipii, J.R. A shock for Voyager 2. Nature 454, 38–39, 2008; Gurnett, D.A., Kurth, W.S. Intense plasma waves at and near the solar wind termination shock. Nature 454, 78–80, 2008. doi: 10.1038/nature07023; Wang, L., Lin, R.P., Larson, D.E., Luhmann, J.G. Domination of heliosheath pressure by shock-accelerated pickup ions from observations of neutral atoms. Nature 454, 81–83, 2008. doi: 10.1038/nature07068.14; Burlaga, L.F., Ness, N.F., Acuna, M.H., et al. Magnetic fields at the solar wind termination shock. Nature 454, 75–77, 2008. doi: 10.1038/nature07029; Richardson, J.D., Kasper, J.C., Wang, C., et al. Cool heliosheath plasma and deceleration of the upstream solar wind at the termination shock. Nature 454, 63–66, 2008. doi: 10.1038/nature07024; Stone, E.C., Cummings, A.C., McDonald, F.B., et al. An asymmetric solar wind termination shock. Nature 454, 71–74, 2008. doi: 10.1038/nature07022; Decker, R.B., Krimigis, S.M., Roelof, E.C., et al. Mediation of the solar wind termination shock by non-thermal ions. Nature 454, 67–70, 2008. doi: 10.1038/nature 07030), respectively, and entered to the inner heliosheath.     

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.     

Heliospheric asymmetries due to the action of the interstellar magnetic field

We discuss the asymmetry of the heliospheric discontinuities obtained from the analysis of 3D modeling of the solar wind (SW) interaction with local interstellar medium (LISM). The flow of charged particles is governed by the ideal MHD equations and the flow of neutral particles is described by the Boltzmann equation. The emphasis is made on the asymmetries of the termination shock (TS) and the heliopause under the combined action of the interstellar and interplanetary magnetic fields (ISMF and IMF) in the presence of neutral hydrogen atoms whose transport through the heliosphere is modeled kinetically, using a Monte Carlo approach. We show that the deflection of neutral hydrogen flow from its original direction in the unperturbed LISM is highly anisotropic and evaluate a possible angle between the hydrogen deflection plane measured in the SOHO SWAN experiment and the plane containing the ISMF and LISM velocity vectors for different ISMF strengths. It is shown that the ISMF of a strength greater than 4 μG can account for the 10 AU difference in the TS heliocentric difference observed during its crossing by the Voyager 1 and Voyager 2 spacecraft, which however results in a larger discrepancy between the calculated and observed velocity distributions. The effect of a strong ISMF on the distribution of plasma quantities in the inner heliosheath and on 2–3 kHz radio emission is discussed.     

Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

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.     

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.     

太阳风速度的源表面二维结构

以IPS观测数据为基础，在考虑了加速剖面对太阳风源位置和速度大小的影响下，探讨了1984年1745-1755Carrington周太阳风速度在源表面上的二维分布.结果显示：观测和理论结果均显示出此期间北半球太阳风速度明显高于南半球，特别是经度0一30°，300-360°N的高速流均有不同程度的从极区向南侵的现象.低速带与磁中性线有一定程度的类似即均为单峰结构.但理论和观测给出的低速带的峰值位置与磁中性线的峰值位置略有差异，而且磁中性线的报幅大于观测和理论的结果.     

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.     

Coupling of the interstellar and interplanetary magnetic fields at the heliospheric interface: The effect of neutral hydrogen atoms

We present numerical results showing the effect of neutral hydrogen atoms on the solar wind (SW) interaction with the local interstellar medium (LISM), where the interstellar magnetic field (ISMF) is coupled to the interplanetary magnetic field (IMF) at the surface of the heliopause. The IMF on the inner boundary surrounding the Sun is specified in the form of a Parker spiral and self-consistently develops in accordance with the SW motion inside the heliopause. The model of the SW–LISM interaction involves both plasma and neutral components which are treated as fluids. The configuration is chosen where the ISMF is orthogonal to the LISM velocity and tilted 60° to the ecliptic plane. This orientation of the magnetic field is a possible explanation of the 2–3 kHz emission data which is believed to originate ahead of the heliopause. It is shown that the topology of the heliospheric current sheet is substantially affected by the ISMF. The interaction pattern dependence on the neutral hydrogen density is analyzed.     

Spectral properties of backscattered solar Ly-α radiation in the heliosphere: A search for heliospheric boundary effects

We explore the imprints of the solar wind interaction with the local interstellar medium on spectral properties of backscattered solar Ly-α radiation. We employ a newly developed effective model for the interstellar H atom velocity distribution function in the heliosphere ( Katushkina and Izmodenov, 2011). The model takes into account both global effects of H atom perturbations at the heliospheric boundaries and local (i.e. within 10–20 AU from the Sun) effects of the solar ionization, charge exchange, solar gravitation and radiation pressure.     

Spatial variation of the supersonic thermal plasma flow downstream of the termination shock

In this paper we start from the most recently observed fact that the solar wind plasma after passage over the termination shock is still supersonic with a Mach number of about 2. To explain this unexpected phenomenon and to predict the evolution of properties of the downstream plasma flow we here consider a two-fluid proton plasma with pick-up protons as a separate suprathermal, second proton fluid. We then formulate a self-consistent system of hydrodynamical conservation equations coupling the two fluids by dynamical and thermodynamical coupling terms and taking into account the effects of newly incorporated protons due to charge exchange with the H-atoms in the heliosheath. This then allows us to predict that in the most probable case the solar wind protons will become subsonic over a distance of about 30 AU downstream of the shock. As we can also show, it may, however, happen that the plasma mixture later again reconverts to a supersonic signature and has to undergo a second shock before meeting the heliopause.     

Interplanetary medium – A dusty plasma

The average mass of dust per volume in space equals that of the solar wind so that the interplanetary medium should provide an obvious region to study dust plasma interactions. While dust collective behavior is typically not observed in the interplanetary medium, the dust component rather consists of isolated grains screened by and interacting with the plasma. Space measurements have revealed several phenomena possibly resulting from dust plasma interactions, but most of the dust plasma interactions are at present not quantified. Examples are the production of neutrals and pick-up ions from the dust, dust impact generated field variations at spacecraft and magnetic field variations possibly caused by solar wind interacting with dust trails. Since dust particles carry a surface charge, they are exposed to the Lorentz force in the interplanetary magnetic field and for grains of sub-micrometer sizes acceleration can be substantial.     

Correlations between neutral and ionized solar wind

We report results of a statistical study correlating ionized solar wind (ISW) fluxes observed by ACE during late 2000 and throughout 2001 with neutral solar wind (NSW) fluxes observed by IMAGE/LENA over the same period. The average correlation coefficient between the neutral and ionized solar wind is 0.66 with correlations greater than 0.80 occurring about 29% of the time. Correlations appear to be driven by high solar wind flux variability, similar to results obtained by in situ multi-spacecraft correlation studies. In this study, however, IMAGE remains inside the magnetosphere on over 95% of its orbits. As a function of day of year, or equivalently ecliptic longitude, the slope of the relationship between the neutral solar wind flux and the ionized solar wind flux shows an enhancement near the upstream direction, but the symmetry point appears shifted toward higher ecliptic longitudes than the interstellar neutral (ISN) flow direction by about 20°. The estimated peak interstellar neutral upstream density inside of 1 AU is about 7 × 10⁻³ cm⁻³.     

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.     

Temporal and spatial variation of equatorial ionization anomaly by using multistation ionosonde data for the 19th solar cycle over the Indian region

The paper deals with the study of temporal and spatial variation of equatorial ionization anomaly (EIA) phenomenon along with its dependence on solar activity and season during the 19th solar cycle by using seven Indian ionosonde stations. Present study is an attempt to carry out the comprehensive study of EIA by using the limited number of ground based instruments. This has been achieved by performing the Gaussian fitting over the latitudinal distribution of F2-region critical frequency (foF2) data. Results reveal that the phenomenon of EIA has a strong dependence on solar activity and seasons. The EIA crest exhibits the feature of latitudinal shifting and expansion with increasing solar activity. It is found out that the effect of solar cycle and seasons on EIA is local time dependent. The observations were also compared with the IRI-2001 model predictions and results reveal that the model values are in general agreement with the observed values with some discrepancies, particularly during the high solar activity period and morning sector. The results have been discussed in the light of relative contribution from transequatorial interhemispheric neutral wind and strength of equatorial fountain process during different local time, season and solar activity levels. Furthermore, an attempt is made to parameterize the location and foF2 of the EIA crest by using the regression analysis. These results can be used to predict the latitudinal position and foF2 of the EIA crest for any given 12-month running average sunspot number (R12).