Near-Earth Interplanetary Coronal Mass Ejections During Solar Cycle 23 (1996 – 2009): Catalog and Summary of Properties

In a previous study (Cane and Richardson, J. Geophys. Res. 108(A4), SSH6-1, 2003), we investigated the occurrence of interplanetary coronal mass ejections in the near-Earth solar wind during 1996 – 2002, corresponding to the increasing and maximum phases of solar cycle 23, and provided a “comprehensive” catalog of these events. In this paper, we present a revised and updated catalog of the ≈300 near-Earth ICMEs in 1996 – 2009, encompassing the complete cycle 23, and summarize their basic properties and geomagnetic effects. In particular, solar wind composition and charge state observations are now considered when identifying the ICMEs. In general, these additional data confirm the earlier identifications based predominantly on other solar wind plasma and magnetic field parameters. However, the boundaries of ICME-like plasma based on charge state/composition data may deviate significantly from those based on conventional plasma/magnetic field parameters. Furthermore, the much studied “magnetic clouds”, with flux-rope-like magnetic field configurations, may form just a substructure of the total ICME interval.     

Modeling Magnetic Tower Jets in the Laboratory

The twisting of magnetic fields threading an accretion system can lead to the generation on axis of toroidal field loops. As the magnetic pressure increases, the toroidal field inflates, producing a flow. Collimation is due to a background corona, which radially confines this axially growing “magnetic tower”. We investigate the possibility of studying in the laboratory the dynamics, confinement and stability of magnetic tower jets. We present two-dimensional resistive magnetohydrodynamic simulations of radial arrays, which consist of two concentric electrodes connected radially by thin metallic wires. In the laboratory, a radial wire array is driven by a 1 MA current which produces a hot, low density background plasma. During the current discharge a low plasma beta (β < 1), magnetic cavity develops in the background plasma (β is the ratio of thermal to magnetic pressure). This laboratory magnetic tower is driven by the magnetic pressure of the toroidal field and it is surrounded by a shock envelope. On axis, a high density column is produced by the pinch effect. The background plasma has >rsim1, and in the radial direction the magnetic tower is confined mostly by the thermal pressure. In contrast, in the axial direction the pressure rapidly decays and an elongated, well collimated magnetic-jet develops. This is later disrupted by the development of m = 0 instabilities arising in the axial column.     

Enhancement of Solar Energetic Particles During a Shock – Magnetic Cloud Interacting Complex Structure

The behavior of solar energetic particles (SEPs) in a shock – magnetic cloud interacting complex structure observed by the Advanced Composition Explorer (ACE) spacecraft on 5 November 2001 is analyzed. A strong shock causing magnetic field strength and solar wind speed increases of about 41 nT and 300 km s⁻¹, respectively, propagated within a preceding magnetic cloud (MC). It is found that an extraordinary SEP enhancement appeared at the high-energy (≥10 MeV) proton intensities and extended over and only over the entire period of the shock – MC structure passing through the spacecraft. Such SEP behavior is much different from the usual picture that the SEPs are depressed in MCs. The comparison of this event with other top SEP events of solar cycle 23 (2000 Bastille Day and 2003 Halloween events) shows that such an enhancement resulted from the effects of the shock – MC complex structure leading to the highest ≥10 MeV proton intensity of solar cycle 23. Our analysis suggests that the relatively isolated magnetic field configuration of MCs combined with an embedded strong shock could significantly enhance the SEP intensity; SEPs are accelerated by the shock and confined into the MC. Further, we find that the SEP enhancement at lower energies happened not only within the shock – MC structure but also after it, probably owing to the presence of a following MC-like structure. This is consistent with the picture that SEP fluxes could be enhanced in the magnetic topology between two MCs, which was proposed based on numerical simulations by Kallenrode and Cliver (Proc. 27th ICRC 8, 3318, 2001b).     

Laboratory Experiment of Plasma Flow Around Magnetic Sail

To propel a spacecraft in the direction leaving the Sun, a magnetic sail (MagSail) blocks the hypersonic solar wind plasma flow by an artificial magnetic field. In order to simulate the interaction between the solar wind and the artificially deployed magnetic field produced around a magnetic sail spacecraft, a laboratory simulator was designed and constructed inside a space chamber. As a solar wind simulator, a high-power magnetoplasmadynamic arcjet is operated in a quasisteady mode of 0.8 ms duration. It can generate a simulated solar wind that is a high-speed (above 20 km/s), high-density (10¹⁸ m⁻³) hydrogen plasma plume of ∼0.7 m in diameter. A small coil (2 cm in diameter), which is to simulate a magnetic sail spacecraft and can obtain 1.9-T magnetic field strength at its center, was immersed inside the simulated solar wind. Using these devices, the formation of a magnetic cavity (∼8 cm in radius) was observed around the coil, which indicates successful simulation of the plasma flow of a MagSail in the laboratory.     

Hybrid Simulation of Collisionless Shock Formation in Support of Laboratory Experiments at Unr

The problem of producing collisionless shocks in the laboratory is of great interest for space and astrophysical plasmas. One approach is based on the idea of combining strong magnetic field (up to 100 Tesla) created during a Z-pinch discharge with a plasma flow produced in the process of the interaction of a laser pulse with a solid target. In support of laboratory experiments we present hybrid simulations of the interaction of the plasma flow with frozen in it magnetic field, with the spherical obstacle. Parameters of the flow correspond to a laser plasma ablation produced in the laboratory during irradiation of the target by a 3 J laser. Magnetic fields in the plasma flow and around the obstacle are created by the currents produced by the pulse power ZEBRA voltage generator. With the appropriate set of initial conditions imposed on the flow collisionless shocks can be created in such a system. Using independent generators for plasma flow and magnetic field allows for the exploration of a wide range of shock parameters. We present simulations of the formation of supercritical collisionless shock relevant to the experiment, performed with the 2D version of the hybrid code based on the CAM-CL algorithm [Planet. Space Sci. 51, 649, 2003].     

Polar Coronal Structures and the Global Magnetic Field Evolution Through the Cycle

We compare the shape and position of some plasma formations visible in the polar corona with the cyclic evolution of the global magnetic field. The first type of object is polar crown prominences. A two-fold decrease of the height of polar crown prominences was found during their poleward migration from the middle latitudes to the poles before a polar magnetic field reversal. The effect could be assigned to a decrease of the magnetic field scale. The second type of object is the polar plumes, ray like structures that follow magnetic field lines. Tangents to polar ray structures are usually crossed near some point, “a magnetic focus,” below the surface. The distance q between the focus and the center of the solar disk changes from the maximum value about 0.65 R _⊙ at solar minimum activity to the minimum value about 0.45 R _⊙ at solar maximum. At first glance this behaviour seems to be contrary to the dynamics of spherical harmonics of the global magnetic field throughout a cycle. We believe that the problem could be resolved if one takes into account not only scale changes in the global magnetic field but also the phase difference in the cyclic variations of large-scale and small-scale components of the global field.     

A Study of the Frequency of Occurrence of Large-Fluence Solar Proton Events and the Strength of the Interplanetary Magnetic Field

It has been shown previously that the number of very-large-fluence solar proton events inferred for the period since 1561 were more frequent at times of low solar activity (e.g., following the recovery from the Maunder minimum), than in the present epoch of high solar activity. An inverse dependence is demonstrated between the probability of observation of the very large-fluence solar proton events and the strength of the interplanetary magnetic field derived from empirical predictions. Using the observed dependence, it is predicted and demonstrated that large-fluence solar proton events have been observed at Earth more frequently near the recurrent minima of the solar activity cycle in the past than during the present epoch. We show that these results are explicable in terms of the linear dependence of the Alfvén velocity upon the strength of the interplanetary magnetic field, leading to higher shock compression ratios in the past. These results indicate that this aspect of “solar weather” will be significantly influenced by the prevailing strength of the interplanetary magnetic field, and that recurrence of solar conditions similar to those of the solar activity minimum of solar cycles 12–14 (1878.9–1913.6) would be accompanied by a factor of ∼4 increase in the occurrence of large-fluence solar proton events.     

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.     

The Bastille day Magnetic Clouds and Upstream Shocks: Near-Earth Interplanetary Observations

The energetic charged particle, interplanetary magnetic field, and plasma characteristics of the `Bastille Day' shock and ejecta/magnetic cloud events at 1 AU occurring over the days 14–16 July 2000 are described. Profiles of MeV (WIND/LEMT) energetic ions help to organize the overall sequence of events from the solar source to 1 AU. Stressed are analyses of an outstanding magnetic cloud (MC2) starting late on 15 July and its upstream shock about 4 hours earlier in WIND magnetic field and plasma data. Also analyzed is a less certain, but likely, magnetic cloud (MC1) occurring early on 15 July; this was separated from MC2 by its upstream shock and many heliospheric current sheet (HCS) crossings. Other HCS crossings occurred throughout the 3-day period. Overall this dramatic series of interplanetary events caused a large multi-phase magnetic storm with min Dst lower than −300 nT. The very fast solar wind speed (≥ 1100 km s⁻¹) in and around the front of MC2 (for near average densities) was responsible for a very high solar wind ram pressure driving in the front of the magnetosphere to geocentric distances estimated to be as low as ≈ 5 R _E, much lower than the geosynchronous orbit radius. This was consistent with magnetic field observations from two GOES satellites which indicated they were in the magnetosheath for extended times. A static force-free field model is used to fit the two magnetic cloud profiles providing estimates of the clouds' physical and geometrical properties. MC2 was much larger than MC1, but their axes were nearly antiparallel, and their magnetic fields had the same left-handed helicity. MC2's axis and its upstream shock normal were very close to being perpendicular to each other, as might be expected if the cloud were driving the shock at the time of observation. The estimated axial magnetic flux carried by MC2 was 52×10²⁰ Mx, which is about 5 times the typical magnetic flux estimated for other magnetic clouds in the WIND data over its first 4 years and is 17 times the flux of MC1. This large flux is due to both the strong axially-directed field of MC2 (46.8 nT on the axis) and the large radius (R ₀=0.189 AU) of the flux tube. MC2's average speed is consistent with the expected transit time from a halo-CME to which it is apparently related.     

Excitation of Electromagnetic Flute Modes in the Process of Interaction of Plasma Flow with Inhomogeneous Magnetic Field

Laboratory experiments on the interaction of a plasma flow, produced by laser ablation of a solid target with the inhomogeneous magnetic field from the Zebra pulsed power generator demonstrated the presence of strong wave activity in the region of the flow deceleration. The deceleration of the plasma flow can be interpreted as the appearance of a gravity-like force. The drift due to this force can lead to the excitation of flute modes. In this paper a linear dispersion equation for the excitation of electromagnetic flute-type modes with plasma and magnetic field parameters, corresponding to the ongoing experiments is examined. Results indicate that the wavelength of the excited flute modes strongly depends on the strength of the external magnetic field. For magnetic field strengths ∼0.1 MG the excited wavelengths are larger than the width of the laser ablated plasma plume and cannot be observed during the experiment. But for magnetic field strengths ∼1 MG the excited wavelengths are much smaller and can then be detected.     

Vortex structure of neutron stars with triplet neutron superfluidity

The vortex structure of the “npe” phase of neutron stars with a ³P₂ superfluid neutron condensate of Cooper pairs is discussed. It is shown that, as the star rotates, superfluid neutron vortex filaments described by a unitary ordering parameter develop in the “npe” phase. The entrainment of superconducting protons by the rotating superfluid neutrons is examined. The entrainment effect leads to the appearance of clusters of proton vortices around each neutron vortex and generates a magnetic field on the order of 10¹² G. ³P₂ neutron vortex filaments combine with quark semi-superfluid vortex filaments at the boundary of the “npe” and “CFL” phases. At the boundary of the “Aen” and “npe” phases, they combine with ¹S₀ neutron vortex filaments. In this way, a unified vortex structure is formed. The existence of this structure and its collective elastic oscillations explain the observed oscillations in the angular rotation velocity of pulsars.     

Streamer Belt and Chains as the Main Sources of Quasi-Stationary Slow Solar Wind

Synoptic maps of white-light coronal brightness from SOHO/LASCO C2 and distributions of solar wind velocity obtained from interplanetary scintillation are studied. Regions with velocity V≈300 – 450 km s⁻¹ and increased density N>10 cm⁻³, typical of the “slow” solar wind originating from the belt and chains of streamers, are shown to exist at Earth’s orbit, between the fast solar wind flows (with a maximum velocity V _max≈450 – 800 km s⁻¹). The belt and chains of streamers are the main sources of the “slow” solar wind. As the sources of “slow” solar wind, the contribution from the chains of streamers may be comparable to that from the streamer belt.     

The Sun’s Large-Scale Magnetic Field and Its Long-Term Evolution

The Sun’s large-scale external field is formed through the emergence of magnetic flux in active regions and its subsequent dispersal over the solar surface by differential rotation, supergranular convection, and meridional flow. The observed evolution of the polar fields and open flux (or interplanetary field) during recent solar cycles can be reproduced by assuming a supergranular diffusion rate of 500 – 600 km² s⁻¹ and a poleward flow speed of 10 –20 m s⁻¹. The nonaxisymmetric component of the large-scale field decays on the flow timescale of ∼1 yr and must be continually regenerated by new sunspot activity. Stochastic fluctuations in the longitudinal distribution of active regions can produce large peaks in the Sun’s equatorial dipole moment and in the interplanetary field strength during the declining phase of the cycle; by the same token, they can lead to sudden weakenings of the large-scale field near sunspot maximum (Gnevyshev gaps). Flux transport simulations over many solar cycles suggest that the meridional flow speed is correlated with cycle amplitude, with the flow being slower during less active cycles.     

From a Weak to a Strong Comet – 3d Global Hybrid Simulation Studies

Plasma structures resulting from the solar wind interaction with weakcomets are discussed. Numericalsimulations using a newly developed hybrid code are presented. The simulations are primarily applied to quantitative data for cometWirtanen, which will be the target of the Rosetta mission. It isexpected that Wirtanen is very weak during the first encounter. The main purpose is the discussion of the different features of the plasmaenvironment, such as the structured cycloidal plasma tail and non-linear Mach cones typical for weak comets and their relation tostructures like shocklets, bow shock, diamagnetic cavity and the “classical”magnetotail found at stronger comets. Furthermore, the sensitivity ofthese various features in dependence on the plasma parameters is investigated.     

The solar wind interaction with Venus through the eyes of the Pioneer Venus Orbiter

Venus has no internal magnetic dynamo and thus its ionosphere and hot oxygen exosphere dominate the interaction with the solar wind. The solar wind at 0.72 AU has a dynamic pressure that ranges from 4.5 nPa (at solar max) to 6.6 nPa (at solar min), and its flow past the planet produces a shock of typical magnetosonic Mach number 5 at the subsolar point. At solar maximum the pressure in the ionospheric plasma is sufficient to hold off the solar wind at an altitude of 400 km above the surface at the subsolar point, and 1000 km above the terminators. The deflection of the solar wind occurs through the formation of a magnetic barrier on the inner edge of the magnetosheath, or shocked solar wind. Under typical solar wind conditions the time scale for diffusion of the magnetic field into the ionosphere is so long that the ionosphere remains field free and the barrier deflects almost all the incoming solar wind. Any neutral atoms of the hot oxygen exosphere that reach the altitude of the magnetosheath are accelerated by the electric field of the flowing magnetized plasma and swept along cycloidal paths in the antisolar direction. This pickup process, while important for the loss of the Venus atmosphere, plays a minor role in the deceleration and deflection of the solar wind. Like at magnetized planets, the Venus shock and magnetosheath generate hot electrons and ions that flow back along magnetic field lines into the solar wind to form a foreshock. A magnetic tail is created by the magnetic flux that is slowed in the interaction and becomes mass-loaded with thermal ions.The structure of the ionosphere is very much dependent on solar activity and the dynamic pressure of the solar wind. At solar maximum under typical solar wind conditions, the ionosphere is unmagnetized except for the presence of thin magnetic flux ropes. The ionospheric plasma flows freely to the nightside forming a well-developed night ionosphere. When the solar wind pressure dominates over the ionospheric pressure the ionosphere becomes completely magnetized, the flow to the nightside diminishes, and the night ionosphere weakens. Even at solar maximum the night ionosphere has a very irregular density structure. The electromagnetic environment of Venus has not been well surveyed. At ELF and VLF frequencies there is noise generated in the foreshock and shock. At low altitude in the night ionosphere noise, presumably generated by lightning, can be detected. This paper reviews the plasma environment at Venus and the physics of the solar wind interaction on the threshold of a new series of Venus exploration missions.     

Interstellar matter and the location of the shock front

The initially supersonic flow of the solar wind passes through a magnetic shock front where its velocity is supposed to be reduced to subsonic values. The location of this shock front is primarily determined by the energy density of the external interstellar magnetic field and the momentum density of the solar wind plasma. Interstellar hydrogen penetrating into the heliosphere undergoes charge exchange processes with the solar wind protons and ionization processes by the solar EUV radiation. This results in an extraction of momentum from the solar wind plasma. Changes of the geometry and the location of the shock front due to this interaction are studied in detail and it is shown that the distance of the magnetic shock front from the Sun decreases from 200 to 80 AU for an increase of the interstellar hydrogen density from 0.1 to 1.0 cm⁻³. The geometry of the shock front is essentially spherical with a pronounced embayment in the direction opposite to the approach of interstellar matter which depends very much on the temperature of the interstellar gas. Due to the energy loss by the interaction with neutral matter the solar wind plasma reduces its velocity with increasing distance from the Sun. This modifies Parker's solution of a constant solar wind velocity.     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.     

Simulating Astrophysical Jets in Laboratory Experiments

Pulsed-power technology and appropriate boundary conditions have been used to create simulations of magnetically driven astrophysical jets in a laboratory experiment. The experiments are quite reproducible and involve a distinct sequence. Eight initial flux tubes, corresponding to eight gas injection locations, merge to form the jet, which lengthens, collimates, and eventually kinks. A model developed to explain the collimation process predicts that collimation is intimately related to convection and pile-up of frozen-in toroidal flux convected with the jet. The pile-up occurs when there is an axial non-uniformity in the jet velocity so that in the frame of the jet there appears to be a converging flow of plasma carrying frozen-in toroidal magnetic flux. The pile-up of convected flux at this “stagnation region” amplifies the toroidal magnetic field and increases the pinch force, thereby collimating the jet.     

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.     

Gamma-ray lines and neutrons from solar flares

An analysis, with a representative (canonical) example of solar-flare-generated equatorial disturbances, is presented for the temporal and spatial changes in the solar wind plasma and magnetic field environment between the Sun and one astronomical unit (AU). Our objective is to search for first order global consequences rather than to make a parametric study. The analysis - an extension of earlier planar studies - considers all three plasma velocity and magnetic field components (V _r, V_φ, V₀, and B _r, B₀, B_φ) in any convenient heliospheric plane of symmetry such as the ecliptic plane, the solar equatorial plane, or the heliospheric equatorial plane chosen for its ability (in a tilted coordinate system) to order northern and southern hemispheric magnetic topology and latitudinal solar wind flows. Latitudinal velocity and magnetic field gradients in and near the plane of symmetry are considered to provide higher-order corrections of a specialized nature and, accordingly, are neglected, as is dissipation, except at shock waves. The representative disturbance is examined for the canonical case in which one describes the temporal and spatial changes in a homogeneous solar wind caused by a solar-flare-generated shock wave. The ‘canonical’ solar flare is assumed to produce a shock wave that has a velocity of 1000 km s^#X2212;1 at 0.08 AU. We have examined all plasma and field parameters at three radial locations: central meridian and 33° W and 90° W of the flare's central meridian. A higher shock velocity (3000 km s^#X2212;1) was also used to demonstrate the model's ability to simulate a strongly-kinked interplanetary field. Among the global (first-order) results are the following: (i) incorporation of a small meridional magnetic field in the ambient magnetic spiral field has negligible effect on the results; (ii) the magnetic field demonstrates strong kinking within the interplanetary shocked flow, even reversed polarity that - coupled with low temperature and low density - suggests a viable explanation for observed ‘magnetic clouds’ with accompanying double-streaming of electrons at directions ～ 90° to the heliocentric radius.