Interplanetary magnetic field structure and solar wind parameters as inferred from solar magnetic field observations and by using a numerical 2-D MHD model

An attempt is made to infer parameters of the solar corona and the solar wind by means of a numerical, self-consistent MHD simulation. Boundary conditions for the magnetic field are given from the observations of the large-scale magnetic field at the Sun. A two-region, planar (the ecliptic plane is assumed) model for the solar wind flow is considered. Region I of transonic flow is assumed to cover the distances from the solar surface up to 10R _S (R _S is the radius of the Sun). Region II of supersonic, super-Alfvénic flow extends between 10R _S and the Earth's orbit. Treatment for region I is that for a mixed initial-boundary value problem. The solution procedure is similar to that discussed by Endler (1971) and Steinolfson, Suess, and Wu (1982): a steady-state solution is sought as a relaxation to the dynamic equilibrium of an initial state. To obtain a solution to the initial value problem in region II with the initial distribution of dependent variables at 10R _S (deduced from the solution for region I), a numerical scheme similar to that used by Pizzo (1978, 1982) is applied. Solar rotation is taken into account for region II; hence, the interaction between fast and slow solar wind streams is self-consistently treated. As a test example for the proposed formulation and numerical technique, a solution for the problem similar to that discussed by Steinolfson, Suess, and Wu (1982) is obtained. To demonstrate the applicability of our scheme to experimental data, solar magnetic field observations at Stanford University for Carrington rotation 1682 are used to prescribe boundary conditions for the magnetic field at the solar surface. The steady-state solution appropriate for the given boundary conditions was obtained for region I and then traced to the Earth's orbit through region II. We compare the calculated and spacecraft-observed solar wind velocity, radial magnetic field, and number density and find that general trends during the solar rotation are reproduced fairly well although the magnitudes of the density in comparison are vastly different.     

Understanding the Behavior of the Heliospheric Magnetic Field and the Solar Wind During the Unusual Solar Minimum Between Cycles 23 and 24

The properties of the heliospheric magnetic field and the solar wind were substantially different in the unusual solar minimum between Cycles 23 and 24: the magnetic-field strength was substantially reduced, as were the flow properties of the solar wind, such as the mass flux. Explanations for these changes are offered that do not require any substantial reconsideration of the general understandings of the behavior of the heliospheric magnetic field and the solar wind that were developed in the minimum of Cycle 22?–?23. Solar-wind composition data are used to demonstrate that there are two distinct regions of solar wind: solar wind likely to originate from the stalk of the streamer belt (the highly elongated loops that underlie the heliospheric current sheet), and solar wind from outside this region. The region outside the streamer-stalk region is noticeably larger in the minimum of Cycle 23?–?24; however, the increased area can account for the reduction in the heliospheric magnetic-field strength in this minimum. Thus, the total magnetic flux contained in this region is the same in the two minima. Various correlations among the solar-wind mass flux and coronal electron temperature inferred from solar-wind charge states were developed for the Cycle 22?–?23 solar minimum. The data for the minimum of Cycle 23?–?24 suggest that the correlations still hold, and thus the basic acceleration mechanism is unchanged in this minimum.     

RADIO MAPS OF THE SOLAR WIND TRANSITION REGION

Occultation studies of near-Sun plasmas using several natural sources simultaneously result in large-scale patterns, radio maps of the solar wind flow. Large radio telescopes of the P.N. Lebedev Physical Institute, Pushino, were used. Previously it had been shown that the plasma acceleration and traverse of the sound barrier proceed in an extended region, the transition region of the solar wind, located at radial distances of about 10–40 solar radii from the Sun. The 1989–1994 experiments showed that the evolution of the transition region geometry is very close to that of the optical corona. On the other hand, the plasma flow structures characteristic of the transition region persist in the course of the 11-year cycle, which demonstrates the existence of some specific mechanism of the solar wind acceleration, independent of wide variations of the general solar activity state. These experimental facts are discussed in connection with the existing theoretical approaches.     

Diagnostics of solar wind flows

The formation of the solar wind, the plasma flows from the Sun, is studied by new methods that have been developed in recent years. Experiments on circumsolar plasma sounding at radial heliocentric distances of ～2.5–60R _⊙ form their basis. Experimental data are used to construct the correlation diagrams-the location of the boundary of the transonic solar wind transition region versus the magnetic field strength in the region of the flow sources. The 2000–2004 correlation diagrams reveal flows of six types that differ by the magnetic field structure in their sources. During the decline of solar activity in 2003–2004, the evolution of the slow solar wind flows has been found to be determined not by the Wolf numbers, but by the total strength of the global magnetic field in the solar corona.     

North-south asymmetry in solar, interplanetary, and geomagnetic indices

Data of hourly interplanetary plasma (field magnitude, solar wind speed, and ion density), solar (sunspot number, solar radio flux), and geomagnetic indices (Kp, Ap) over the period 1970-2010, have been used to examine the asymmetry between the solar field north and south of the heliospheric current sheet (HCS). A persistent yearly north-south asymmetry of the field magnitude is clear over the considered period, and there is no magnetic solar cycle dependence. There is a weak N-S asymmetry in the averaged solar wind speed, exhibited well at times of maximum solar activities. The solar plasma is more dense north of the current sheet than south of it during the second negative solar polarity epoch (qA < 0). Moreover, the N - S asymmetry in solar activity (Rz) can be statistically highly significant. The sign of the average N - S asymmetry depends upon the solar magnetic polarity. The annual magnitudes of N - S asymmetry depend positively on the solar magnetic cycle. Most of the solar radio flux asymmetries occurred during the period of positive IMF polarity.     

Ion escape from Mars as a function of solar wind conditions: A statistical study

The influence of solar EUV and solar wind conditions on ion escape at Mars is investigated using ion data from the Aspera-3 instrument on Mars Express, combined with solar wind proxy data obtained from the Mars Global Surveyor (MGS) spacecraft. A solar EUV flux proxy based on data from the Earth position, scaled and shifted in time for Mars, is used to study relatively long time scale changes related to solar EUV variability. Data from May 2004 until November 2005 has been used. A clear dependence on the strength of the subsolar magnetic field as inferred from MGS measurements is seen in the ion data. The region of significant heavy ion flows is compressed and the heavy ion flux density is higher for high subsolar magnetic field strength. Because of the difference in outflow area, the difference in estimated total outflow is somewhat less than the difference in average flux density. We confirm previous findings that escaping planetary ions are mainly seen in the hemisphere into which the solar wind electric field is pointed. The effect is more pronounced for the high subsolar magnetic field case.The average ion motion has a consistent bias towards the direction of the solar wind electric field, but the main motion is in the antisunward direction. The antisunward flow velocity increases with tailward distance, reaching above at 2 to 3 martian radii downtail from Mars for O⁺ ions. Different ion species reach approximately the same bulk flow energy. We did not find any clear correlation between the solar EUV flux and the ion escape distribution or rate, probably because the variation of the solar EUV flux over our study interval was too small. The results indicate that the solar wind and its magnetic field directly interacts with the ionosphere of Mars, removing more ions for high subsolar magnetic field strength. The interaction region and the tail heavy ion flow region are not perfectly shielded from the solar wind electric field, which accelerates particles over relatively large tail distances.     

Comparison of single-site interplanetary scintillation solar wind speed structure with coronal features

Interplanetary Scintillation (IPS) Observations were made during the period 1984–1990 using a single radio telescope at 103 MHz situated at Thaltej (Ahmedabad), India. Solar wind speeds were estimated using a recently developed method based on matching the observed IPS spectra with model solar wind spectra for Kolmogorov turbulence. The best-fit speeds derived are traced back to a source surface, and average velocity maps are made for each year, averaging over a number of Carrington rotations. It is found that the resulting single-site, large-scale IPS speed structure agrees well with that derived from 3-site observations from earlier workers. The IPS speed structure during this period was compared with other coronal features. Nearly 85% of the observed high-speed regions were associated with coronal holes. At solar minimum, in 1986, a quasi-sinusoidal, narrow belt of slow solar wind was observed which matched well with the neutral line structure of the solar magnetic field and the belt of active centers. Near solar maximum, in 1990, the speed structure was chaotic, similar to that of the neutral line, with low speed regions appearing all over the source surface.     

The influence of momentum addition on the acceleration of the solar and stellar winds

The influence of the momentum addition, which may be associated with the average or fluctuation transverse component of the magnetic field or others, on the acceleration the solar wind or stellar wind is studied in a local streamtube. The results show that the larger the momentum addition the stronger the acceleration of the wind. For example, if the typical transverse magnetic field is about 0.1 of the longitudinal field, the velocity of the solar wind at 1 AU may be increased by 40%. The coronal hole may be considered as a streamtube, the presence of a high stream from the coronal hole may be explained by the existence of an average or fluctuation transverse magnetic field in the streamtube. A similar conclusion may be applied to the polar region, where the velocity of the solar wind will be larger than elsewhere as if there is a transverse component of magnetic field, as well as to the stellar wind. The influence of other parameters on the acceleration of the solar wind is also discussed. From the viewpoint of the solar wind mechanism, the present paper shows that the momentum addition in the subsonic flow region can increase the velocity of the solar wind at 1 AU.     

Substorm recovery phase

The recovery phase of the magnetospheric substorm is studied numerically by means of a two-dimensional time-dependent nonlinear resistive MHD code. The initial configuration was chosen from the earlier numerical model in which the magnetospheric substorm was driven by the solar wind plasmas. In order to study the recovery phase, the entering solar wind energy flux was reduced when the magnetospheric substorm was in its expansive phase. The system was found to respond instantly to this change and the result showed many characteristic features related to the recovery phase including the tailward motion of thex-point of the reconnected magnetic field lines and the restoration of a tail-like configuration of the magnetic field. Thex-point moved at almost the same speed of the plasma flow in the upstream region, which was considerably smaller than the speed of the plasma jetting or the speed of the plasmoid. As the recovery phase progressed, the plasma jetting across thex-point was reduced very much in the Earthside region. Although the plasma flow was generally in the Earthward direction in the Earthside region of thex-point, the tailward flow was also found near thex-point. The current density was reduced near thex-point and the neutral sheet was broadened in the recovered region. The plasma sheet also became thick in this region. During the recovery of the substorm, the energy conversion rate, both in the form of plasma acceleration and the Joule heating, was reduced. These results on the recovery phase together with the earlier simulation result on the expansive phase indicate that driven reconnection can be a viable mechanism for the magnetospheric substorm including the recovery phase.     

The space environment of Mercury at the times of the second and third MESSENGER flybys

The second and third flybys of Mercury by the MESSENGER spacecraft occurred, respectively, on 6 October 2008 and on 29 September 2009. In order to provide contextual information about the solar wind properties and the interplanetary magnetic field (IMF) near the planet at those times, we have used an empirical modeling technique combined with a numerical physics-based solar wind model. The Wang–Sheeley–Arge (WSA) method uses solar photospheric magnetic field observations (from Earth-based instruments) in order to estimate the inner heliospheric radial flow speed and radial magnetic field out to 21.5 solar radii from the Sun. This information is then used as input to the global numerical magnetohydrodynamic model, ENLIL, which calculates solar wind velocity, density, temperature, and magnetic field strength and polarity throughout the inner heliosphere. WSA-ENLIL calculations are presented for the several-week period encompassing the second and third flybys. This information, in conjunction with available MESSENGER data, aid in understanding the Mercury flyby observations and provide a basis for global magnetospheric modeling. We find that during both flybys, the solar wind conditions were very quiescent and would have provided only modest dynamic driving forces for Mercury's magnetospheric system.     

Sharp Changes of Solar Wind Ion Flux and Density Within and Outside Current Sheets

Analysis of the Interball-1 spacecraft data (1995 – 2000) has shown that the solar wind ion flux sometimes increases or decreases abruptly by more than 20% over a time period of several seconds or minutes. Typically, the amplitude of such sharp changes in the solar wind ion flux (SCIFs) is larger than 0.5×10⁸ cm⁻² s⁻¹. These sudden changes of the ion flux were also observed by the Solar Wind Experiment (SWE), on board the Wind spacecraft, as the solar wind density increases and decreases with negligible changes in the solar wind velocity. SCIFs occur irregularly at 1 AU, when plasma flows with specific properties come to the Earth’s orbit. SCIFs are usually observed in slow, turbulent solar wind with increased density and interplanetary magnetic field strength. The number of times SCIFs occur during a day is simulated using the solar wind density, magnetic field, and their standard deviations as input parameters for a period of five years. A correlation coefficient of ∼0.7 is obtained between the modelled and the experimental data. It is found that SCIFs are not associated with coronal mass ejections (CMEs), corotating interaction regions (CIRs), or interplanetary shocks; however, 85% of the sector boundaries are surrounded by SCIFs. The properties of the solar wind plasma for days with five or more SCIF observations are the same as those of the solar wind plasma at the sector boundaries. One possible explanation for the occurrence of SCIFs (near sector boundaries) is magnetic reconnection at the heliospheric current sheet or local current sheets. Other probable causes of SCIFs (inside sectors) are turbulent processes in the slow solar wind and at the crossings of flux tubes.     

Magnetic Field Strengths and Structures from Radio Observations of Solar Active Regions

Radio observations of some active regions (ARs) obtained with the Nobeyama radioheliograph at λ=1.76cm are used for estimating the magnetic field strength in the upper chromosphere, based on thermal bremsstrahlung. The results are compared with the magnetic field strength in the photosphere from observations with the Solar Magnetic Field Telescope (SMFT) at Huairou Solar Observing Station of Beijing Astronomical Observatory. The difference in the magnetic field strength between the two layers seems reasonable. The solar radio maps of active regions obtained with the Nobeyama radioheliograph, both in total intensity (I-map) and in circular polarizations (V-map), are compared with the optical magnetograms obtained with the SMFT. The comparison between the radio map in circular polarization and the longitudinal photospheric magnetogram of a plage region suggest that the radio map in circular polarization is a kind of magnetogram of the upper chromosphere. The comparison of the radio map in total intensity with the photospheric vector magnetogram of an AR shows that the radio map in total intensity gives indications of magnetic loops in the corona, thus we have a method of defining the coronal magnetic structure from the radio I-maps at λ=1.76 cm. Analysing the I-maps, we identified three components: (a) a compact bright source; (b) a narrow elongated structure connecting two main magnetic islands of opposite polarities (observed in both the optical and radio magnetograms); (c) a wide, diffuse, weak component that corresponds to a wide structure in the solar active region which shows in most cases an S or a reversed S contour, which is probably due to the differential rotation of the Sun. The last two components suggest coronal loops on different spatial scales above the neutral line of the longitudinal photospheric magnetic field.     

The acceleration of minor ion species in the solar wind

This paper provides a comprehensive analysis of the dynamics of the flow of minor ion species in the solar wind under the combined influences of gravity, Coulomb friction (with protons), rotational forces (arising from the Sun's rotation and the interplanetary spiral magnetic field) and wave forces (induced in the minor ion flow by Alfvén waves propagating in the solar wind). It is assumed that the solar wind can be considered as a proton-electron plasma which is, to a first approximation, unaffected by the presence of minor ions. In the dense hot region near the Sun Coulomb friction accelerates minor ions outwards against the gravitational force, part of which is cancelled by the charge-separation electric field. Once the initial acceleration has been achieved, wave and rotational forces assist Coulomb friction in further increasing the minor ion speed so that it becomes comparable with, or perhaps even exceeds, the solar wind speed. A characteristic feature of the non-resonant wave force is that it tends to bring the minor ion flow into an equilibrium where the radial speed matches the Alfvén speed relative to the solar wind speed, whereas Coulomb friction and rotational forces tend to bring the flow into an equilibrium where the radial speed of the minor ions equals the solar wind speed. Therefore, provided that there is sufficient wave energy and Coulomb friction is weak, the minor ion speed can be trapped between these two speeds. This inteststing result is in qualitative agreement with observational findings to the effect that the differential flow speed between helium ions and protons is controlled by the ratio of the solar wind expansion time to the ion-proton collision time. If the thermal speeds of the protons and minor ions are small compared to the Alfvén speed, two stable equilibrium speeds can exist because the rapid decrease in the Coulomb cross-section with increasing differential flow speed allows the non-resonant wave force to balance Coulomb friction at more than one ion speed. However, it must be emphasized that resonant wave acceleration and/or strong ion partial pressure gradients are required to achieve radial speeds of minor ions in excess of the proton speed, since, as is shown in Section 4, the non-resonant wave acceleration on protons and minor ions are identical when their radial speeds are the same, with the result that, in the solar wind, non-resonant wave acceleration tends (asymptotically) to equalize minor ion and proton speeds.     

Solar wind speed and its acceleration inferred using the interplanetary scintillation method in carrington rotation 1753

Solar wind speeds (SWSs) estimated by interplanetary scintillation (IPS) observations during Carrington rotation 1753 are projected onto the so-called source-surface of 2.5 solar radii along magnetic field lines in interplanetary space. The following two working hypotheses are examined from different points of view: (1) The SWS is a weighted mean along the line of sight to a radio source; the weight for the SWS depends on the distance from theP-point, the closest approach to the Sun on the line of sight. (2) The weighting function has a very sharp peak at theP-point, so that the SWS shows a real solar wind speed at theP-point. In both the two cases, the SWSs projected onto the source surface are further projected onto the photosphere along magnetic field lines in the corona. Footpoints of these field lines are inferred as photospheric source regions of the solar wind. The intensity of the Hei (1083 nm) absorption line (HEI) in the chromosphere corresponding to these photospheric sources is interpolated from observational data. The weighted mean of the HEI is calculated in case 1. The HEI corresponding to theP-point is used in case 2. The SWS is compared with the HEI in both the two cases. It is found that the correlation between the SWS and the HEI is better in case 2 than in case 1. It is further inferred by correlation analysis between the SWS and the HEI that the solar wind is accelerated within 27 solar radii on average. Although the data examined in this paper were limited to just one solar rotation, these results suggest that the SWS estimated by the IPS technique corresponds to the solar wind speed near theP-point and the weighting function along the line of sight may have a very sharp peak near theP-point.     

A global 3-D simulation of interplanetary dynamics in June 1991

The global dynamics of the solar wind and interplanetary magnetic field in June 1991 is simulated based on a fully three-dimensional, time-dependent numerical MHD model. The numerical simulation includes eight transient disturbances associated with the major solar flares of June 1991. The unique features of the present simulation are: (i) the disturbances are originated at the coronal base (1R _s) and their propagation through inhomogeneous ambient solar wind is simulated out to 1.5 AU; (ii) as a background for the transients, the global steady-state solar wind structure inferred from the 3-D steady-state model (Usmanov, 1993c) is used. The parameters of the initial pulses are prescribed in terms of the near-Sun shock velocities (as inferred from the metric Type II radio burst observations) relative to the preshock steady-state flow parameters at the flare sites. The computed parameters at the Earth's location for the period 1–18 June, 1991 are compared with the available observations of the interplanetary magnetic field, solar wind velocity, density, and with variation of the geomagnetic activityK _pindex.     

The Bastille day Shocks and Merged Interaction Region

Near 1 AU the solar wind structure associated with the solar flare of 14 July 2000 (Bastille Day) consisted of a large high-speed stream of 15 July and five nearby small streams during a 10-day period. At the leading edge of the large high-speed stream, in less than 6 hours, the flow speed increased from 600 km s⁻¹ to 1100 km s⁻¹, the magnetic field intensity increased from 10 nT to 60 nT, and an interaction region was identified. The interaction region was bounded between the pair of a forward shock F and a reverse shock R. Additional forward shocks were also identified at the leading edge of each of the five smaller streams. This paper presents a magnetohydrodynamics (MHD) simulation using ACE plasma and magnetic field data near 1 AU as input to study the radial evolution of the Bastille Day solar wind event. The two shocks, F and R, propagated in opposite directions away from each other in the solar wind frame and interacted with neighboring shocks and streams; the spatial and temporal extent of the interaction region continued to increase with the heliocentric distance. The solar wind was restructured from a series of streams at 1 AU to a huge merged interaction region (MIR) extending over a period of 12 days at 5.5 AU. Throughout the interior of the MIR bounded by the shock pair F and R the magnetic field intensity was a few times stronger than that outside the MIR. The simulation shows how merging of shocks, collision of shocks, and formation of new shocks contributed to the evolution process.     

Preferred bartels days of high-speed plasma streams in the solar wind

An analysis of 346 high-speed solar wind streams observed at 1 AU during 1964–75 is presented. The analysis shows that a two-sector structure was the dominant feature of the interplanetary magnetic field associated with the high-speed solar wind plasma. The high-speed streams occurred at preferred Bartels days: Positive polarity streams were most frequent near Bartels day 4, negative polarity streams were most frequent near Bartels day 17. Since the solar wind carries with it the photospheric magnetic polarity of the solar source region, the observed distribution of Bartels days must indicate a fundamental property of the distribution of the solar sources of high-speed plasma streams. The observations are explained in terms of a tilted dipole model of the solar-interplanetary field.Proceedings of the 14th ESLAB Symposium on Physics of Solar Variations, 16–19 September 1980, Scheveningen, The Netherlands.     

Phase mixing of standing Alfvén waves with shear flows in solar spicules

Alfvénic waves are thought to play an important role in coronal heating and solar wind acceleration. Here we investigate the dissipation of such waves due to phase mixing at the presence of shear flow and field in the stratified atmosphere of solar spicules. The initial flow is assumed to be directed along spicule axis and to vary linearly in the x direction and the equilibrium magnetic field is taken 2-dimensional and divergence-free. It is determined that the shear flow and field can fasten the damping of standing Alfvén waves. In spite of propagating Alfvén waves, standing Alfvén waves in Solar spicules dissipate in a few periods. As height increases, the perturbed velocity amplitude does increase in contrast to the behavior of perturbed magnetic field. Moreover, it should be emphasized that the stratification due to gravity, shear flow and field are the facts that should be considered in MHD models in spicules.     

Solar Wind Forecasting with Coronal Holes

An empirical model for forecasting solar wind speed related geomagnetic events is presented here. The model is based on the estimated location and size of solar coronal holes. This method differs from models that are based on photospheric magnetograms (e.g., Wang–Sheeley model) to estimate the open field line configuration. Rather than requiring the use of a full magnetic synoptic map, the method presented here can be used to forecast solar wind velocities and magnetic polarity from a single coronal hole image, along with a single magnetic full-disk image. The coronal hole parameters used in this study are estimated with Kitt Peak Vacuum Telescope He I 1083 nm spectrograms and photospheric magnetograms. Solar wind and coronal hole data for the period between May 1992 and September 2003 are investigated. The new model is found to be accurate to within 10% of observed solar wind measurements for its best 1-month period, and it has a linear correlation coefficient of ∼0.38 for the full 11 years studied. Using a single estimated coronal hole map, the model can forecast the Earth directed solar wind velocity up to 8.5 days in advance. In addition, this method can be used with any source of coronal hole area and location data.     

Energization of positive ions in the cometary foreshock region

The energization of positive ions in front of a cometary bow shock is investigated. Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce, among other waves, large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading) ; hence, they can energize the suprathermal ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting ion energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind, such that the cometary bow shock of Halley-type comet may be regarded as a “cosmic ray shock”.