Alfvén waves in the solar atmosphere

We examine the propagation of Alfvén waves in the solar atmosphere. The principal theoretical virtues of this work are: (i) The full wave equation is solved without recourse to the small-wavelength eikonal approximation (ii) The background solar atmosphere is realistic, consisting of an HSRA/VAL representation of the photosphere and chromosphere, a 200 km thick transition region, a model for the upper transition region below a coronal hole (provided by R. Munro), and the Munro-Jackson model of a polar coronal hole. The principal results are:

1. If the wave source is taken to be near the top of the convection zone, where n _H = 5.2 × 10¹⁶ cm^?3, and if B _⊙ = 10.5 G, then the wave Poynting flux exhibits a series of strong resonant peaks at periods downwards from 1.6 hr. The resonant frequencies are in the ratios of the zeroes of J ₀, but depend on B _⊙, and on the density and scale height at the wave source. The longest period peaks may be the most important, because they are nearest to the supergranular periods and to the observed periods near 1 AU, and because they are the broadest in frequency.
2. The Poynting flux in the resonant peaks can be large enough, i.e. P _⊙ ≈ 10⁴–10⁵ erg cm^?2s^?1, to strongly affect the solar wind.
3. ¦δv¦ and ¦δB¦ also display resonant peaks.
4. In the chromosphere and low corona, ¦δv ≈ 7–25 kms^?1 and ¦δB¦ ≈0.3–1.0 G if P _⊙≈10⁴-10⁵ erg cm^?2s^?1.
5. The dependences of ¦δv¦ and ¦δB¦ on height are reduced by finite wavelength effects, except near the wave source where they are enhanced.
6. Near the base, ¦δB¦ ≈ 350–1200 G if P _⊙ ～- 10⁴–10⁵. This means that nonlinear effects may be important, and that some density and vertical velocity fluctuations may be associated with the Alfvén waves.
7. Below the low corona most wave energy is kinetic, except near the base where it becomes mostly magnetic at the resonances.
8. ?₀ < δv ² > v _A or < δB ² > v _A/4π are not good estimators of the energy flux.
9. The Alfvén wave pressure tensor will be important in the transition region only if the magnetic field diverges rapidly. But the Alfvén wave pressure can be important in the coronal hole.

     

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.     

Corotating Shock Waves and the Solar-wind Source of Energetic Ion bundances: Power Laws in $$/Q$Q$

We find that element abundances in energetic ions accelerated by shock waves formed at corotating interaction regions (CIRs) mirror the abundances of the solar wind modified by a decreasing power-law dependence on the mass-to-charge ratio \(A\)/\(Q\) of the ions. This behavior is similar in character to the well-known power-law dependence on \(A\)/\(Q\) of abundances in large gradual solar energetic particles (SEP). The CIR ions reflect the pattern of \(A\)/\(Q\), with \(Q\) values of the source plasma temperature or freezing-in temperature of 1.0?–?1.2 MK typical of the fast solar wind in this case. Thus the relative ion abundances in CIRs are of the form \((A\mbox{/}Q)^{a}\), where \(a\) is nearly always negative and evidently decreases with distance from the shocks, which usually begin beyond 1 AU. For one unusual historic CIR event where \(a \approx 0\), the reverse shock wave of the CIR seems to occur at 1 AU, and these abundances of the energetic ions become a direct proxy for the abundances of the fast solar wind.     

The multiple modes of interaction of the solar wind with a comet as it approaches the Sun

A quasi-steady 1-D hydrodynamic model, with mass addition, has been used to study the various modes of interaction of the solar wind with a medium-bright, H₂O-dominated comet (such as P/Halley) approaching the Sun.At large heliocentric distances (d 5 AU) the solar wind penetrates unimpeded on to the surface. As the comet moves further in, mass loading of the solar wind by heavy ions from the fledgling cometary atmosphere causes it to slow down, thereby causing a significant enhancement of the interplanetary field. Still further in at d 3.14 AU, as the mass loading reaches a critical value, a collision-less standing shock is formed in the solar wind upstream of the nucleus. As d decreases further, the distance of this shock from the nucleus increases. The cometary atmosphere becomes dense enough to stand off the solar wind ahead of the nuclear surface and form a well defined tangential discontinuity surface (or ionopause) only when d reaches the value 2.65 AU. When d 2.65 AU an inner shock could, in principle, also form within the cometary ionosphere, although its existence would depend on the detailed thermodynamics of the cometary ionosphere. Resolution of this question is beyond the scope of the present analysis.The conclusions of the present study would be qualitatively valid for other comets having sizes, surface optical properties and chemical compositions, different from those adopted here. The helio-centric distances at which the various transitions take place from the one mode of solar wind interaction to another, would, of course, be different, with all these distances being smaller for less active comets.     

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.     

The motion of the Martian satellites

An extensive analysis of the motion of Phobos and Deimos from 1877 to 1973 has been fulfilled. The new values of the parameters of the orbital model first developed by Struve have been determined for both satellites. The new sets of the orbital parameters compete with the solutions of similar accuracy found by Wilkins and Sinclair. A secular acceleration in longitude of Phobos is found to be equal to +(0.107±0.011)×10^?7 deg day^?2. The value of the acceleration is little affected when one or another group of oppositions is omitted. The acceleration of Deimos is determined with great uncertainty: +(0.06±0.34)×10^?9 deg day^?2. Values found for the orbital parameters seem to be in good agreement since the mass, oblateness and coordinates of the pole of Mars inferred from the motion of each satellite have similar values in both cases.     

Upper limits on the total radiant energy of solar flares

We establish limits on the total radiant energy of solar flares during the period 1980 February – November, using the solar-constant monitor (ACRIM) on board the Solar Maximum Mission. Typical limits amount to 6 × 10²⁹ erg/s for a 32-second integration time, with 5σ statistical significance, for an impulsive emission; for a gradual component, about 4 × 10³² ergs total radiant energy. The limits lie about an order of magnitude higher than the total radiant energy estimated from the various known emission components, suggesting that no heretofore unknown dominant component of flare radiation exists.     

Ion acceleration and Alfvén wave generation at the Earth’s bow shock

The processes of ion acceleration and Alfvén wave generation by accelerated particles at the Earth’s bow shock are studied within a quasi-linear approach. Steady-state ion and wave spectra are shown to be established in a time of 0.3–4 h, depending on the background level of Alfvénic turbulence in the solar wind. The Alfvén waves produced by accelerated ions are confined within the frequency range 10^?2–1 Hz and their spectral peak with a wave amplitude βB ～ B comparable to the interplanetary magnetic field strength B corresponds to the frequency v = (2–3) × 10^?2 Hz. The high-frequency part of the wave spectrum (v > 0.2 Hz) undergoes damping by thermal ions. The calculated spectra of the accelerated ions and the Alfvén waves generated by them reproduce the main features observed in experiments.     

On Solar-Wind Electron Heating at Large Solar Distances

We study the temperature of electrons advected with the solar wind to large solar distances far beyond 1 AU. Almost nothing is known about the thermodynamics of these electrons from in-situ plasma observations at these distances, and usually it is tacitly assumed that electrons, due to adiabatic behaviour and vanishing heat conduction, rapidly cool off to very low temperatures at larger distances. In this article we show, however, that electrons on their way to large distances undergo non-adiabatic interactions with travelling shocks and solar-wind bulk-velocity jumps and thereby are appreciably heated. Examining this heating process on an average statistical basis, we find that solar-wind electrons first cool down to a temperature minimum, which depending on the occurrence frequency of bulk velocity jumps is located between 3 and 6 AU, but beyond this the lowest electron temperature again starts to increase with increasing solar distance, finally achieving temperatures of about 7×10⁴ K to 7×10⁵ K at the location of the termination shock. Hence these electrons are unexpectedly shown to play an important dynamical role in structuring this shock and in determining the downstream plasma properties.     

Non-thermal processes in large solar flares

We analyze particle acceleration processes in large solar flares, using observations of the August, 1972, series of large events. The energetic particle populations are estimated from the hard X-ray and γ-ray emission, and from direct interplanetary particle observations. The collisional energy losses of these particles are computed as a function of height, assuming that the particles are accelerated high in the solar atmosphere and then precipitate down into denser layers. We compare the computed energy input with the flare energy output in radiation, heating, and mass ejection, and find for large proton event flares that:

1. The ～10–10² keV electrons accelerated during the flash phase constitute the bulk of the total flare energy.
2. The flare can be divided into two regions depending on whether the electron energy input goes into radiation or explosive heating. The computed energy input to the radiative quasi-equilibrium region agrees with the observed flare energy output in optical, UV, and EUV radiation.
3. The electron energy input to the explosive heating region can produce evaporation of the upper chromosphere needed to form the soft X-ray flare plasma.
4. Very intense energetic electron fluxes can provide the energy and mass for interplanetary shock wave by heating the atmospheric gas to energies sufficient to escape the solar gravitational and magnetic fields. The threshold for shock formation appears to be ～10³¹ ergs total energy in >20 keV electrons, and all of the shock energy can be supplied by electrons if their spectrum extends down to 5–10 keV.
5. High energy protons are accelerated later than the 10–10² keV electrons and most of them escape to the interplanetary medium. The energetic protons are not a significant contributor to the energization of flare phenomena. The observations are consistent with shock-wave acceleration of the protons and other nuclei, and also of electrons to relativistic energies.
6. The flare white-light continuum emission is consistent with a model of free-bound transitions in a plasma with strong non-thermal ionization produced in the lower solar chromosphere by energetic electrons. The white-light continuum is inconsistent with models of photospheric heating by the energetic particles. A threshold energy of ～5×10³⁰ ergs in >20 keV electrons is required for detectable white-light emission.

The highly efficient electron energization required in these flares suggests that the flare mechanism consists of rapid dissipation of chromospheric and coronal field-aligned or sheet currents, due to the onset of current-driven Buneman anomalous resistivity. Large proton flares then result when the energy input from accelerated electrons is sufficient to form a shock wave.     

Revivals of a coronal arch

The giant post-flare arch of 6 November 1980 revived 11 hr and 25 hr after its formation. Both these revivals were caused by two-ribbon flares with growing systems of loops. The first two brightenings of the arch were homologous events with brightness maxima moving upwards through the corona with rather constant speed; during all three brightenings the arch showed a velocity pattern with two components: a slow one (8–12 km^?1), related to the moving maxima of brightness, and a fast one (～ 35 km s^?1), the source of which is unknown. During the first revival, at an altitude of 100000 km, temperature in the arch peaked ～ 1 hr, brightness ～ 2 hr, and emission measure ～ 3.5 hr after the onset of the brightening. Thus the arch looks like a magnified flare, with the scales both in size and time increased by an order of magnitude. At ～ 100000 km altitude the maximum temperature was ?14 × 10⁶K, max.n _e? 2.5 × 10⁹cm^?3, and max. energy density ? 11.2 erg cm^?3. The volume of the whole arch can be estimated to 1.1 × 10³⁰ cm³, total energy ?1.2 × 10³¹ erg, and total mass ?4.4 × 10¹⁵g. The density decreased with the increasing altitude and remained below 7 × 10⁹ cm^?3 anywhere in the arch. The arch cooled very slowly through radiation whereas conductive cooling was inhibited. Since its onset the revived arch was subject to energy input within the whole extent of the preexisting arch while a thermal disturbance (a new arch?) propagated slowly from below. We suggest that the first heating of the revived arch was due to reconnection of some of the distended flare loops with the magnetic field of the old preexisting arch. The formation of the ‘post’-flare loop system was delayed and started only some 30–40 min later. Since that time a new arch began to be formed above the loops and the velocities we found reflect this formation.     

A simple shock wave model of supernova remnant IC 443

This paper contains a model of supernova remnant IC 443 and the interstellar gas surrounding it. The basis of this model is the analysis of the motion of non-spherical adiabatic shock waves due to Kompaneets (1960). Observations of adjacent Hi and Hii regions have been used by several authors to determine the density of gas in the neighbourhood of IC 443. The model gives for the explosion energy and age of IC 443, 1.8×10⁵⁰ erg and 13 000 yr, respectively. The expansion velocity of IC 443 using the present model is in reasonabel agreement with observations of the remnant.     

A dynamo theory of solar flares

It is proposed that the solar flare phenomenon can be understood as a manifestation of the electrodynamic coupling process of the photosphere-chromosphere-corona system as a whole. The system is coupled by electric currents, flowing along (both upward and downward) and across the magnetic field lines, powered by the dynamo process driven by the neutral wind in the photosphere and the lower chromosphere. A self-consistent formulation of the proposed coupling system is given. It is shown in particular that the coupling system can generate and dissipate the power of 10²⁹ erg s^#X2212;1 and the total energy of 10³² erg during a typical life time (10³ s) of solar flares. The energy consumptions include Joule heat production, acceleration of current-carrying particles along field lines, magnetic energy storage and kinetic energy of plasma convection. The particle acceleration arises from the development of field-aligned potential drops of 10–150 kV due to the loss-cone constriction effect along the upward field-aligned currents, causing optical, X-ray and radio emissions. The total number of precipitating electrons during a flare is shown to be of order 10³⁷–10³⁸.     

Partial analysis of the flare-prominence of 30 April 1974

A portion of an east limb flare-prominence observed in Hα by NOAA/Boulder and NASA/ MSFC patrol facilities on 30 April 1974 is analyzed. Following a rapid (～2 min) achievement of a maximum mass ejection velocity of about 375 km s^?1, the ascending prominence reached a height of, at least, 2 × 10⁵ km. We use a one-dimensional, time-dependent hydrodynamic theory (Nakagawa et al., 1975) to compute the total mass (～2 × 10¹¹ g) and energy (～4 × 10²⁶erg) ejected during this part of this event. Theoretical aspects of the coronal response are discussed. We conclude that a moderate temperature and density pulse (factors of ten and two, respectively), for a duration of only 3 min, is sufficient for an acceptable simulation of the Hα observations and the likely coronal response to the ascending prominence and flare-related ejections. No attempt was made to simulate the additionally-important spray and surge features which probably contributed a higher level of mass and energy efflux.     

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.     

Numerical evaluation of statistical acceleration for energetic particles in the interplanetary medium at 2.5 and 5.0 AU

Numerical integration of particle trajectories is performed to evaluate the statistical acceleration coefficients D _TT for 1 to 100 MeV protons in a solar wind corotating interaction region (CIR) seen at 2.5 and 5.0 AU. Acceleration is followed in the solar wind reference frame and is due to random wave-particle interactions and to random drift motion in moderate scale field gradients. D _TT due to the first effect reaches a peak value of 4 × 10 ^–7 MeV² s^–1 post shock at 10 MeV at 2.5 AU consistent with estimates based both upon cyclotron resonance and transit time damping theory. D _TT from the second effect is less well established but is of the order of 10^–7 MeV² s^–1 at 10 MeV, 5 AU. A comparison is made between the time constant for statistical acceleration within this CIR and estimates for diffuse shock acceleration and adiabatic deceleration. All three time constants are of the same order, but deceleration is faster than shock acceleration which in turn is faster than statistical acceleration.     

On the velocity of jets powering hot spots similar to those in Cygnus A

Hot spots similar to those in the radio galaxy Cygnus A can be explained by the strong shock produced by a supersonic but classical jet \(\left( {u_{jet}< c/\sqrt 3 } \right)\) . The high integrated radio luminosity (L?2×10⁴⁴ erg s^?1) and the strength of mean magnetic field (B?2×10^?4 G) suggest the hot spots are the downstream flow of a very strong shock which generates the ultrarelativistic electrons of energy ?≥20 MeV. The fully-developed subsonic turbulence amplifies the magnetic field of the jet up to 1.6×10^?4 G by the dynamo effect. If we assume that the post-shock pressure is dominated by relativistic particles, the ratio between the magnetic energy density to the energy density in relativistic particles is found to be ?2×10^?2, showing that the generally accepted hypothesis of equipartition is not valid for hot spots. The current analysis allows the determination of physical parameters inside hot spots. It is found that:

1. The velocity of the upstream flow in the frame of reference of the shock isu ₁?0.2c. Radio observations indicate that the velocity of separation of hot spots isu _sep?0.05c, so that the velocity of the jet isu _jet=u ₁+u _sep?0.25c.
2. The density of the thermal electrons inside the hot spot isn ₂?5×10^?3 e ^? cm^?3 and the mass ejected per year to power the hot spot is ?4M ₀yr^?1.
3. The relativistic electron density is less than 20% of the thermal electron density inside the hot spot and the spectrum is a power law which continues to energies as low as 30 MeV.
4. The energy density of relativistic protons is lower than the energy density of relativistic electrons unlike the situation for cosmic rays in the Galaxy.

     

TWO-FLUID MODEL FOR HEATING OF THE SOLAR CORONA AND ACCELERATION OF THE SOLAR WIND BY HIGH-FREQUENCY ALFVÉN WAVES

A model of the solar corona and wind is developed which includes for the first time the heating and acceleration effects of high-frequency Alfvén waves in the frequency range between 1 Hz and 1 kHz. The waves are assumed to be generated by the small-scale magnetic activity in the chromospheric network. The wave dissipation near the gyro-frequency, which decreases with increasing solar distance, leads to strong coronal heating. The resulting heating function is different from other artificial heating functions used in previous model calculations. The associated thermal pressure-gradient force and wave pressure-gradient force together can accelerate the wind to high velocities, such as those observed by Helios and Ulysses. Classical Coulomb heat conduction is also considered and turns out to play a role in shaping the temperature profiles of the heated protons. The time-dependent two-fluid (electrons and protons) model equations and the time-dependent wave-spectrum equation are numerically integrated versus solar distance out to about 0.3 AU. The solutions finally converge and settle on time-stationary profiles which are discussed in detail. The model computations can be made to fit the observed density profiles of a polar coronal hole and polar plume with the sonic point occurring at 2.4 R and 3.2 R , respectively. The solar wind speeds obtained at 63 R are 740 km s^-1 and 540 km s^-1; the mass flux is 2.1 and 2.2 × 10⁸ cm^-2 s^-1 (normalized to 1 AU), respectively. The proton temperature increases from a value of 4 × 10⁵ K at the lower boundary to 2 × 10⁶ K in the corona near 2 R .     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).