Judith T. Karpen

THE MECHANISMS FOR THE ONSET AND EXPLOSIVE ERUPTION OF CORONAL MASS EJECTIONS AND ERUPTIVE FLARES

We have investigated the onset and acceleration of coronal mass ejections (CMEs) and eruptive flares. To isolate the eruption physics, our study uses the breakout model, which is insensitive to the energy buildup process leading to the eruption. We performed 2.5D simulations with adaptive mesh refinement that achieved the highest overall spatial resolution to date in a CME/eruptive flare simulation. The ultra-high resolution allows us to separate clearly the timing of the various phases of the eruption. Using new computational tools, we have determined the number and evolution of all X- and O-type nulls in the system, thereby tracking both the progress and the products of reconnection throughout the computational domain. Our results show definitively that CME onset is due to the start of fast reconnection at the breakout current sheet. Once this reconnection begins, eruption is inevitable; if this is the only reconnection in the system, however, the eruption will be slow. The explosive CME acceleration is triggered by fast reconnection at the flare current sheet. Our results indicate that the explosive eruption is caused by a resistive instability, not an ideal process. Moreover, both breakout and flare reconnections begin first as a form of weak tearing characterized by slowly evolving plasmoids, but eventually transition to a fast form with well-defined Alfvenic reconnection jets and rapid flux transfer. This transition to fast reconnection is required for both CME onset and explosive acceleration. We discuss the key implications of our results for CME/flare observations and for theories of magnetic reconnection.

Formation of the slow solar wind in a coronal streamer

We have investigated a magnetohydrodynamic mechanism that accounts for several fundamental properties of the slow solar wind, in particular its variability, latitudinal extent, and bulk acceleration. In view of the well-established association between the streamer belt and the slow wind, our model begins with a simplified representation of a streamer beyond the underlying coronal helmet: a neutral sheet embedded in a plane fluid wake. This wake-neutral sheet configuration is characterized by two parameters that vary with distance from the Sun: the ratio of the cross-stream velocity scale to the neutral sheet width, and the ratio of the typical Alfven velocity to the typical flow speed far from the neutral sheet. Depending on the values of these parameters, our linear theory predicts that three kinds of instability can develop when this system is perturbed: a tearing instability and two ideal fluid instabilities with different cross-stream symmetries (varicose and sinuous). In the innermost, magnetically dominated region beyond the helmet cusp, we find that the streamer is resistively and ideally unstable, evolving from tearing-type reconnection in the linear regime to a nonlinear varicose fluid instability. Traveling magnetic islands are formed which are similar to features recently revealed by the large-angle spectroscopic coronagraph on the joint European Space Agency/NASA Solar and Heliospheric Observatory (SOHO) [Brueckner et al., 1995]. During this process, the center of the wake is accelerated and broadened slightly. Past the Alfven point, where the kinetic energy exceeds the magnetic energy, the tearing mode is suppressed, but an ideal sinuous fluid mode can develop, producing additional acceleration up to typical slow wind speeds and substantial broadening of the wake. Farther from the Sun, the system becomes highly turbulent as a result of the development of ideal secondary instabilities, thus halting the acceleration and producing strong filamentation throughout the core of the wake. We discuss the implications of this model for the origin and evolution of the slow solar wind, and compare the predicted properties with current observations from SOHO.

Coronal current-sheet formation - The effect of asymmetric and symmetric shears

A 2.5D numerical code is used to investigate the results of an asymmetric shear imposed on a potential quadrupolar magnetic field under two sets of atmospheric boundary conditions - a low-beta plasma with line tying at the base, similar to the line-tied analytic model, and a hydrostatic-equilibrium atmosphere with solar gravity, typical of the observed photosphere-chromosphere interface. The low-beta simulation confirms the crucial role of the line-tying assumption in producting current sheets. The effects of a symmetric shear on the same hydrostatic-equilibrium atmosphere is examined, using more grid points to improve the resolution of the current structures which form along the flux surfaces. It is found that true current sheets do not form in the corona when a more realistic model is considered. The amount of Ohmic dissipation in the thick currents is estimated to be two to four orders of magnitude below that required to heat the corona. It is concluded that magnetic topologies of the type examined here do not contribute significantly to coronal heating.

Large-amplitude Longitudinal Oscillations in a Solar Filament

We have developed the first self-consistent model for the observed large-amplitude oscillations along filament axes that explains the restoring force and damping mechanism. We have investigated the oscillations of multiple threads formed in long, dipped flux tubes through the thermal nonequilibrium process, and found that the oscillation properties predicted by our simulations agree with the observed behavior. We then constructed a model for the large-amplitude longitudinal oscillations that demonstrates that the restoring force is the projected gravity in the tube where the threads oscillate. Although the period is independent of the tube length and the constantly growing mass, the motions are strongly damped by the steady accretion of mass onto the threads by thermal nonequilibrium. The observations and our model suggest that a nearby impulsive event drives the existing prominence threads along their supporting tubes, away from the heating deposition site, without destroying them. The subsequent oscillations occur because the displaced threads reside in magnetic concavities with large radii of curvature. Our model yields a powerful seismological method for constraining the coronal magnetic field and radius of curvature of dips. Furthermore, these results indicate that the magnetic structure is most consistent with the sheared-arcade model for filament channels.

The Kelvin-Helmholtz instability in photospheric flows - Effects of coronal heating and structure

A series of hydrodynamic numerical simulations has been used to investigate the nonlinear evolution of driven, subsonic velocity shears under a range of typical photospheric conditions. These calculations show that typical photospheric flows are susceptible to the Kelvin-Helmholtz instability (KHI), with rapid nonlinear growth times that are approximately half of a typical granule lifetime. The KHI produces vortical structures in intergranule lanes comparable to a typical fluxule radius; this is precisely the correct scale for maximum power transfer to the corona.

Nonlinear evolution of radiation-driven thermally unstable fluids

The nonlinear evolution of a radiation-driven thermally unstable planar fluid is simulated numerically using a semiimplicit finite-difference algorithm. When the equilibrium state of the fluid is perturbed by random initial excitation of the velocity field, dense, cool, two-dimensional structures are found to form in a rarer, warmer surrounding medium. The nonlinear phase of evolution is characterized by the turbulent contraction of the condensed region, accompanied by a significant increase in the amount of energy radiated. It is found that, if the random velocity perturbation has a sufficiently large amplitude, the fluid will not form condensed structures. Finally, the relationship of these results to observations of the solar chromosphere, transition region, and corona is discussed. 30 references.

Transition to turbulence in solar surges

Transition to turbulence in magnetohydrodynamic (MHD) tearing jets has been invoked as a mechanism underlying some of the complex behavior observed in solar surges, including deceleration of the upflowing plasma and temporal correlations with types I and III radio bursts. In this paper we investigate a possible mechanism for this transition: three-dimensional secondary instabilities on two-dimensional saturated states. We find through linear analysis that these MHD configurations -- in particular, the tearing jet -- are secondarily unstable, with the dominant energy transfer from the one-dimensional field into the 3-dimensional fields. Using nonlinear simulations, we also investigate the system evolution after the secondary modes attain finite amplitude. When the tearing jet transitions to turbulence, the total kinetic energy drops rapidly corresponding to the deceleration of the jet. The electric field grows rapidly as the primary mode saturates and the three-dimensional secondary mode develops, and then decays quickly as the tearing jet becomes turbulent, providing a possible explanation for the finite duration of the associated meter-wave bursts. The electric field decays as the magnetic and velocity fields both decay. The system is dominated at late times by spanwise modes, which strongly resemble the magnetic field-aligned filamentary flows characteristic of many surges.

The effects of Kelvin-Helmholtz instability on resonance absorption layers in coronal loops

One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Nonlocal thermal transport in solar flares

A flaring solar atmosphere is modeled assuming classical thermal transport, locally limited thermal transport, and nonlocal thermal transport. The classical, local, and nonlocal expressions for the heat flux yield significantly different temperature, density, and velocity profiles throughout the rise phase of the flare. Evaporation of chromospheric material begins earlier in the nonlocal case than in the classical or local calculations, but reaches much lower upward velocities. Much higher coronal temperatures are achieved in the nonlocal calculations owing to the combined effects of delocalization and flux limiting. The peak velocity and momentum are roughly the same in all three cases. A more impulsive energy release influences the evolution of the nonlocal model more than the classical and locally limited cases. 23 references.

Dynamic spectral characteristics of thermal models for solar hard X-ray bursts

The dynamic spectral characteristics of the thermal model for solar hard X-ray bursts recently proposed by Brown et al. (1979) (BMS) are investigated. It is pointed out that this model, in which a single source is heated impulsively and cooled by anomalous conduction across an ion-acoustic turbulent thermal front, predicts that the total source emission measure should rise as the temperature falls. This prediction, which is common to all conductively cooled single-source models, is contrary to observations of many simple spike bursts. It is proposed, therefore, that the hard X-ray source may consist of a distribution of many small impulsively-heated kernels, each cooled by anomalous conduction, with lifetimes shorter than current burst data temporal resolution. In this case the dynamic spectra of bursts are governed by the dynamic evolution of the kernel production process, such as magnetic-field dissipation in the tearing mode. An integral equation is formulated, the solution of which yields information on this kernel production process, from dynamic burst spectra, for any kernel model.With a BMS-type kernel model in one-dimensional form, the derived instantaneous spectra are limited in hardness to spectral indices γ ≳ 4 for any kernel production process, due to the nature of the conductive cooling. Ion-acoustic conductive cooling in three dimensions, however, increases the limiting spectral hardness to γ ≳ 3. Other forms of anomalous conduction yield similar results but could permit bursts as hard as γ ≳ 2, consistent with the hardest observed.The contribution to the X-ray spectrum from the escaping tail of high-energy kernel electrons in the BMS model is calculated in various limits. If this tail dissipates purely collisionally, for example, its thick-target bremsstrahlung can significantly modify the kernel spectrum at the high-energy end. The energetics of this dynamic dissipation model for thermal hard X-ray bursts also are briefly discussed.