Evolution of fast and slow shock interactions in the inner heliosphere

We use a one-dimensional, time-dependent adaptive grid MHD code to study the interaction between fast and slow shocks in the solar wind. Our results show that: (1) a forward slow shock (FSS) can be destroyed by a forward fast shock (FFS) that overtakes it from behind; (2) two propagating FSSs can merge into a stronger FSS; (3) a strong FSS can survive by following a strong forward fast shock; and (4) the strength of a FSS is decreased by following an FFS. These simulation results reproduce an important feature of the Helios observations (Richter, 1987) where transient fast shocks were more frequently followed within a few hours by slow shock ype discontinuities rather than by fast reverse shocks.     

Numerical simulation of slow shocks in the solar wind

The evolutionary state of slow forward shock waves is examined with the use of two MHD numerical codes. Our study is intended to be exploratory rather than a detailed parametric one. The first code is one-dimensional (with three components of velocity and magnetic field) which is used to follow a slow shock that propagates into a positive gradient of density versus distance. It is found that the slow shock evolves into an extraneous (intermediate) shock wave. The second code has a spherical, one-dimensional, planar geometry (with two velocity and magnetic field components) which is used to follow a spiral interplanetary magnetic field. It is found that a slow shock type perturbation can generate a forward slow shock; a fast forward shock is generated in the front of the slow shock; a contact discontinuity is formed behind the slow shock, and a compound nonlinear MHD wave is formed behind the contact discontinuity with a fast reverse shock formed further behind. Thus, we demonstrate that the evolution of a slow shock into (solely) a fast shock, as suggested by Whang (1987), is much more complicated.     

Low density drivers of strong interplanetary shocks

The theory that most, if not all, interplanetary shocks are caused by coronal mass ejections (CMEs) faces serious problems in accounting for the strongest shocks. The difficulties include (i) a remarkable absence of very strong shocks during solar maximum 1980 when CMEs were prolific, (ii) unrealistic initial speeds near the Sun for impulsive models, (iii) the absence of rarefaction zones behind the shocks and (iv) sustained high speed flows following shocks which are not easily explained as consequences of CME eruptions. Observations of the proton temperature near 1 AU indicate that strong shock drivers have properties similar to high speed streams emitted by coronal holes. Eruptions of fast solar wind from coronal holes influenced by solar activity can explain the occurrence of the strongest interplanetary shocks.     

Interaction of interplanetary shocks with the bow shock

Fast forward interplanetary (IP) shocks have been identified as a source of large geomagnetic disturbances. However, the shocks can evolve in the solar wind, they are modified by interaction with the bow shock and during their propagation through the magnetosheath. A few previous papers refer the inclination and deceleration of the IP shock front in this region. Our contribution continues this effort and presents the study of an IP shock interaction with the bow shock. Since the bow shock is a reversed fast shock, the interaction of the IP shock and bow shock is a problem of interaction of two fast MHD shocks.

We compare profiles of magnetic field and plasma parameters observed by several spacecraft in the solar wind and magnetosheath with the profiles of the same parameters resulting from the MHD numerical model. The MHD model suggests that the interaction of an IP shock with the bow shock results in an inward bow shock displacement that is followed by its outward motion. Such motion will result in an indentation propagating along the bow shock surface. This scenario is confirmed by multipoint observations. Moreover, the model confirms also previous suggestions on the IP shock deceleration in the magnetosheath.     

On the overtaking interaction of typical shock waves in the solar wind flow

The interaction of traveling fast solar shock waves with other fast shock waves generated previously is considered in terms of magnetohydrodynamics for various solar wind parameters. The shocks are not piston ones and move freely in the flow. The magnetic structure in the interplanetary magnetic field emerging after the shock interaction is shown to correspond to the well-known magnetic configuration commonly observed on spacecraft or the classical Hundhausen R model. A head-on collision of solar shock waves with the boundary of a magnetic cloud is considered. It is pointed out that a slow shockwave refracted into the magnetic cloud can appear at an oblique collision of the shock with the cloud boundary. The results clarify our understanding of the available spacecraft data.     

Structure and stationarity of quasi-perpendicular shocks: Numerical simulations

This paper presents an overview of numerical simulation studies of fast collisionless shocks and compares these simulation results with observations of the Earth's bow shock and theoretical works. Especially, we review the structure and stationarity of the supercritical quasi-perpendicular shocks. In situ observations indicate that these shocks are generally quasi-stationary whereas full particle simulations as well as hybrid simulations often present a strong nonstationary behavior, a shock self-reformation. The simulation results, along with theoretical and observational works, suggest that the classical models of the quasi-stationary structure generated by reflected protons or by dispersive whistlers are not generally applicable for the supercritical quasi-perpendicular shocks and other phenomena are to be included into the model to ensure the observed quasi-stationarity: The role of a small scale turbulence and shock ripples is investigated. The downstream turbulence and the electron dynamics in the quasi-perpendicular shocks are also discussed.     

Alfvén-magnetosonic waves interaction in the solar corona

The nonlinear propagation of the Alfvén and magnetosonic waves in the solar corona is investigated in terms of model equations. Due to viscous effects taken into account the propagation of the fast wave itself is governed by Burgers type equations possessing both expansion and compression shock solutions. Numerical simulations show that both parallely and perpendicularly propagating fast waves can steepen into shocks if their amplitudes are in excess of some sizeable fraction of the Alfvén velocity. However, if the magnetic field changes linearly in the perpendicular direction, then formation of perpendicular shocks can be hindered. The Alfvén waves exhibit a tendency to drive both the slow and fast magnetosonic waves whose propagation is described by linearized Boussinesq type equations with ponderomotive terms due to the Alfvén wave. The limits of the slow and fast waves are investigated.     

Energetic protons associated with a forward-reverse interplanetary shock pair at 1 A.U.

A forward-reverse interplanetary shock was observed on 25 March 1969 by the magnetometer and plasma detector on the HEOS-1 satellite. This relatively rare event was described by Chao et al (1972) who concluded that the shock pair was formed at a distance 0.10–0.13 A.U. upstream of the Earth as a result of the interaction between a fast and a slow solar wind streams. Simultaneous observations of 1 MeV solar proton fluxes were also performed on HEOS-1. A characteristic intensity peak was observed as the forward shock passed by the spacecraft. The evolution of the proton intensity, together with a detailed analysis of anisotropies and pitch angle distributions show a complex dynamic picture of the effect of the forward shock on the ambient proton population. Significant changes in particle fluxes are seen to be correlated with fluctuations in the magnetic field. It is suggested that simple geometrical models of shock-associated acceleration should be expanded to include the effect of magnetic fluctuations on particle fluxes. The interaction region limited by the forward and reverse shocks contained a large variety of magnetic fluctuations. Following the tangential discontinuity separating the fast solar wind stream from the preceding slow stream, a sunward flow was observed in the proton data, followed by a small but significant drop in intensity prior to the reverse shock.     

MHD study of temporal and spatial evolution of simulated interplanetary shocks in the ecliptic plane within 1 AU

We utilize a 21/2-D MHD time-dependent model to perform a parametric study of interplanetary shock propagation to 1 AU. The input conditions are represented by the following variables:(1) initial shock velocity, (2) duration of the driving pulse, and (3) width of the pulse at the near-Sun position (18 solar radii). The total net energy added to the solar wind was calculated for each pulse. The forward shock's travel time to, and the peak dynamic pressure at, 1 AU as a function of location along the shock front have been studied over a range of total input pulse energies from 10²⁹ to 10³² ergs. For input pulses with modest angular width and temporal duration, we find that the propagation of the resulting interplanetary fast forward shock waves depends primarily upon the net input energy. The dependence of the transit time upon energy is a power law with a -1/3 index which corresponds to the classical, piston driven case. Reverse shocks are also formed behind all but the lowest energy shocks. Their properties, although also a function of input energy, depend upon the specific values of the input pulse shock velocity, width and duration. We also briefly discuss the propagation of the shocks out to 1 AU, and the conditions for which the interplanetary shocks depart from being symmetric about the input pulse central meridian due to magnetic and dynamic effects.     

Shock surfing acceleration

Analytical and numerical analysis identify shock surfing acceleration as an ideal pre-energization mechanism for the slow pick-up ions at quasiperpendicular shocks. After gaining sufficient energy by shock surfing, pick-up ions undergo diffusive acceleration to reach their observed energies. Energetic ions upstream of the cometary bow shock, acceleration of solar energetic particles by magnetosonic waves in corona, ion enhancement in interplanetary shocks, generation of anomalous cosmic rays from interstellar pick-up ions at the termination shock are some of the cases where shock surfing acceleration apply. Inclusion of the lower-hybrid wave turbulence into the laminar model of shock surfing can explain the preferential acceleration of heavier particles as observed by Voyager at the termination shock. At relativistic energies, unlimited acceleration of ions is theoretically possible; because for sufficiently strong shocks main limitation of the mechanism, caused by the escape of accelerated particles downstream of the shock during acceleration no longer exists.     

Limit shocks of relativistic magnetohydrodynamics

In this paper we show that switch-on and switch-off shocks are allowed by the shock equations of relativistic MHD and have similar properties to their Newtonian counterparts. Just like in Newtonian MHD they are limits of fast and slow shock solutions and as such they may be classified as weakly evolutionary shocks.     

Geoeffectiveness of interplanetary shocks during solar minimum (1995–1996) and solar maximum (2000)

Plasma and magnetic field parameter variations across fast forward interplanetary shocks are analyzed during the last solar cycle minimum (1995–1996, 15 shocks), and maximum year 2000 (50 shocks). It was observed that the solar wind velocity and magnetic field strength variation across the shocks were the parameters better correlated with Dst. Superposed epoch analysis centered on the shock showed that, during solar minimum, B _z profiles had a southward, long-duration variation superposed with fluctuations, whereas in solar maximum the B _z profile presented 2 peaks. The first peak occurred 4 hr after the shock, and seems to be associated with the magnetic field disturbed by the shock in the sheath region. The second peak occurred 19 hr after the shock, and seems to be associated with the ejecta fields. The difference in shape and peak in solar maximum (Dst peak =−50 nT, moderate activity) and minimum (Dst peak =−30 nT, weak activity) in average Dst profiles after shocks are, probably, a consequence of the energy injection in the magnetosphere being driven by different interplanetary southward magnetic structures. A statistical distribution of geomagnetic activity levels following interplanetary shocks was also obtained. It was observed that during solar maximum, 36% of interplanetary shocks were followed by intense (Dst≤−100 nT) and 28% by moderate (−50≤Dst <−100 nT) geomagnetic activity. During solar minimum, 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity, respectively. Thus, during solar maximum a higher relative number of interplanetary shocks might be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50–60% to be followed by intense/moderate geomagnetic activity.     

The structure of radiative slow-mode shocks

We investigate the structure of slow-mode MHD shocks in a plasma where both radiation and thermal conduction are important. In such a plasma a slow shock dissociates into an extended foreshock, an isothermal subshock, and a downstream radiative cooling region. Our analysis, which is both numerical and analytical, focuses on the nearly switch-off shocks which are generated by magnetic reconnection in a strong magnetic field. These shocks convert magnetic energy into kinetic energy and heat, and we find that for typical flare conditions about f of the conversion occurs in the subshock while the remaining 1/3 occurs in the foreshock. We also find that no stable, steady-state solutions exist for radiative slow shocks unless the temperature in the radiative region downstream of the subshock falls below 10⁵ K. These results suggest that about 2/3 of the magnetic energy released in flare loops is released at the top of the loop, while the remaining 1/3 is released in the legs of the loop.     

Properties of Stream Interactions at One AU During 1995 – 2004

A stream interaction region (SIR) forms when a fast solar stream overtakes a slow stream, leading to structure that evolves as an SIR moves away from the Sun. Based on Wind (1995 – 2004) and ACE (1998 – 2004) in situ observations, we have conducted a comprehensive survey of SIRs at one AU, including a separate assessment of the longer-lasting corotating interaction regions (CIRs) that recur on more than one solar rotation. In all there are 196 CIRs, accounting for about 54% of the 365 SIRs. The largest proportion of CIRs to SIRs (64%) appears in 1999, and the smallest proportion (49%) is in 2002. Over the ten years, the annual number of SIR events varies little, from 32 up to 45. On average, the occurrence rate of shocks at SIRs at one AU is about 24%. Seventy percent of the SIRs with shocks have only forward shocks, more than twice the percentage of SIRs with only reverse shocks. This preponderance of forward shocks is consistent with the deflections of forward and reverse shocks relative to the ecliptic plane. In order to help address the effect of SIRs and CIRs on geomagnetic activity, we determine the solar-cycle variation of the event duration, scale size, the change in velocity from slow stream to fast stream, and the solar-cycle variation of the maximum magnetic field, peak total perpendicular pressure, and other properties. These statistics also provide a baseline for future studies at other heliocentric distances and for validating heliospheric models. Electronic Supplementary Material Supplementary material is available for this article at     

Synchrotron emission in the fast cooling regime: which spectra can be explained?

We consider the synchrotron emission from relativistic shocks assuming that the radiating electrons cool rapidly (either through synchrotron or any other radiation mechanism). It is shown that the theory of synchrotron emission in the fast cooling regime can account for a wide range of spectral shapes. In particular, the magnetic field, which decays behind the shock front, brings enough flexibility to the theory to explain the majority of gamma-ray burst spectra even in the parameter-free fast cooling regime. Also, we discuss whether location of the peak in observed spectral energy distributions of gamma-ray bursts and active galactic nuclei can be made consistent with predictions of diffusive shock acceleration theory, and find that the answer is negative. This result is a strong indication that a particle injection mechanism, other than the standard shock acceleration, works in relativistic shocks.     

Propagation of Fast Coronal Mass Ejections and Shock Waves Associated with Type II Radio-Burst Emission: An Analytic Study

Coronal mass ejections (CMEs) are large-scale eruptive events in the solar corona. Once they are expelled into the interplanetary (IP) medium, they propagate outwards and “evolve” interacting with the solar wind. Fast CMEs associated with IP shocks are a critical subject for space weather investigations. We present an analytic model to study the heliocentric evolution of fast CME/shock events and their association with type II radio-burst emissions. The propagation model assumes an early stage where the CME acts as a piston driving a shock wave; beyond this point the CME decelerates, tending to match the ambient solar wind speed and its shock decays. We use the shock speed evolution to reproduce type II radio-burst emissions. We analyse four fast CME halo events that were associated with kilometric type II radio bursts, and in-situ measurements of IP shock and CME signatures. The results show good agreement with the dynamic spectra of the type II frequency drifts and the in-situ measurements. This suggests that, in general, IP shocks associated with fast CMEs evolve as blast waves approaching 1 AU, implying that the CMEs do not drive their shocks any further at this heliocentric range.     

Solar Cycle Variation of the Interplanetary Forward Shock Drivers Observed at 1 AU

Forecasting space weather more accurately from solar observations requires an understanding of the variations in physical properties of interplanetary (IP) shocks as solar activity changes. We examined the characteristics (occurrence rate, physical parameters, and types of shock driver) of IP shocks. During the period of 1995 – 2001, a total of 249 forward IP shocks were observed. In calculating the shock parameters, we used the solar wind data from Wind at the solar minimum period (1995 – 1997) and from ACE since 1998 including the solar maximum period (1999 – 2001). Most of IP shocks (68%) are concentrated in the solar maximum period. The values of physical quantities of IP shocks, such as the shock speed, the sonic Mach number, and the ratio of plasma density compression, are larger at solar maximum than at solar minimum. However, the ratio of IMF compression is larger at solar minimum. The IP shock drivers are classified into four groups: magnetic clouds (MCs), ejecta, high speed streams (HSSs), and unidentified drivers. The MC is the most dominant and strong shock driver and 150 out of total 249 IP shocks are driven by MCs. The MC is a principal and very effective shock driver not only at solar maximum but also at solar minimum, in contrast to results from previous studies, where the HSS is considered as the dominant IP shock driver.     

磁通量喷发的大气响应——Ⅲ.太阳附近激波的演化

作为本系列的最后一篇,本文研究日冕大气中的激波演化.结果表明,日冕背景大气和磁场的非均匀性对激波结构和演化起着决定作用。太阳附近形成的慢激波-快磁声波系统,将在向外传播的过程中演变为以中间激波作为必要组成部分的混合激波。该混合激波在沿电流片传播时能继续维持,而在沿单极开放磁场传播时则进一步演变为纯快激波。     

The influence of the Mach number on the stability of radiative shocks

We study the stability properties of hydrodynamic shocks with finite Mach numbers. The linear analysis supplements previous analyses which took the strong shock limit. We derive the linearized equations for a general specific heat ratio as well as temperature and density power-law cooling functions, corresponding to a range of conditions relevant to interstellar atomic and molecular cooling processes. Boundary conditions corresponding to a return to the upstream temperature  ( R = 1)  and to a cold wall  ( R = 0)  are investigated. We find that for Mach number   M > 5  , the strong shock overstability limits are not significantly modified. For   M < 3  , however, shocks are considerably more stable for most cases. In general, as the shock weakens, the critical values of the temperature power-law index (below which shocks are overstable) are reduced for the overtones more than for the fundamental, which signifies a change in basic behaviour. In the   R = 0  scenario, however, we find that the overstability regime and growth rate of the fundamental mode are increased when cooling is under local thermodynamic equilibrium. We provide a possible explanation for the results in terms of a stabilizing influence provided downstream but a destabilizing effect associated with the shock front. We conclude that the regime of overstability for interstellar atomic shocks is well represented by the strong shock limit unless the upstream gas is hot. Although molecular shocks can be overstable to overtones, the magnetic field provides a significant stabilizing influence.     

Instability of isothermal stellar wind bowshocks

We present hydrodynamical simulations illustrating the instability of stellar wind bowshocks in the limit of an isothermal equation of state. In this limit, the bowshock is characterized by a thin dense shell bounded on both sides by shocks. In a time-averaged sense the shape of this bowshock shell roughly matches the steady state solution of Wilkin (1996)[ApJ, 459, L31], although the apex of the bowshock can deviate in or out by a factor of two or more. The shape of the bowshock is distorted by large amplitude kinks with a characteristic wavelength of order the standoff distance from the star. The instability is driven by a strong shear flow within the shock-bounded shell, suggesting an origin related to the nonlinear thin-shell instability. This instability occurs when both the forward bowshock and the reverse wind shock are effectively isothermal and the star is moving through the interstellar medium with a Mach number greater than a few. This work therefore suggests that ragged, clumpy bowshocks should be expected to surround stars with a slow, dense wind (which leads to rapid cooling behind the reverse wind shock), whose velocity with respect to the surrounding interstellar medium is of order 60 km s⁻¹ (leading both to rapid cooling behind the forward bowshock and sufficiently high Mach numbers to drive the instability).