Nonlinear Numerical Simulations of Magneto-Acoustic Wave Propagation in Small-Scale Flux Tubes

We present results of nonlinear, two-dimensional, numerical simulations of magneto-acoustic wave propagation in the photosphere and chromosphere of small-scale flux tubes with internal structure. Waves with realistic periods of three to five minutes are studied, after horizontal and vertical oscillatory perturbations are applied to the equilibrium model. Spurious reflections of shock waves from the upper boundary are minimized by a special boundary condition. This has allowed us to increase the duration of the simulations and to make it long enough to perform a statistical analysis of oscillations. The simulations show that deep horizontal motions of the flux tube generate a slow (magnetic) mode and a surface mode. These modes are efficiently transformed into a slow (acoustic) mode in the v _A<c _S atmosphere. The slow (acoustic) mode propagates vertically along the field lines, forms shocks, and remains always within the flux tube. It might effectively deposit the energy of the driver into the chromosphere. When the driver oscillates with a high frequency, above the cutoff, nonlinear wave propagation occurs with the same dominant driver period at all heights. At low frequencies, below the cutoff, the dominant period of oscillations changes with height from that of the driver in the photosphere to its first harmonic (half period) in the chromosphere. Depending on the period and on the type of the driver, different shock patterns are observed.     

Physical properties of solar chromospheric plages

We compute a new grid of plage models to determine the difference in temperature versus mass column density structure T(m) between plage regions and the quiet solar chromosphere, and to test whether the solar chromosphere is geometrically thinner in plages. We compare partial redistribution calculations of Mg ii h and k and Ca ii K to NRL Skylab observations of Mg ii h and k in six active regions and Ca ii K intensities obtained from spectroheliograms taken at approximately the same time as the Mg ii observations. We find that the plage observations are better matched by models with linear (in log m) temperature distributions and larger values of m ₀ (the mass column density at the 8000 K layer in the chromosphere), than by models with larger low chromosphere temperature gradients but values of m ₀ similar to the quiet Sun. Our derived temperature structures are in agreement with the grid originally proposed by Shine and Linsky, but our analysis is in contrast to the study by Kelch which implies that stellar chromospheric geometrical thickness is not affected by chromospheric activity. We conclude that either the stellar Mg ii observations upon which the Kelch study was based are of poorer quality than had been assumed, or that the spatial averaging of inhomogeneous structures, which is inherent in the stellar data, does not lead to a best fit one-component model similar in detail to that of a stellar or a solar plage.Visiting Astronomer at Kitt Peak National Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.Staff member, Quantum Physics Division, National Bureau of Standards.     

Interpretation of XUV spectroheliograms

Some parameters of chromospheric structure are drawn from recently published XUV spectroheliograms. The HeII emission above the limb arises from the small amount of He⁺ still existing at 10⁶°. The larger amounts of He⁺ in the cooler corona at the poles explain the polar cap absorption in 304. The flat distribution of emission in Oiv and Ov, with a sharp spike at the limb, is caused by the rough structure of the chromosphere and the variable excitation in the emitting spicules. The intensity of the Nevii lines shows that the transition zone between chromosphere and corona is very sharp.This research was supported by the National Aeronautics and Space Administration under Grant NASA NGR 05 002 034.     

Mechanics of Seismic Emission from Solar Flares

Instances of seismic transients emitted into the solar interior in the impulsive phases of some solar flares offer a promising diagnostic tool, both for understanding the physics of solar flares and for the general development of local helioseismology. Among the prospective contributors to flare acoustic emission that have been considered are: i) chromospheric shocks propelled by pressure transients caused by impulsive thick-target heating of the upper and middle chromosphere by high-energy particles, ii) heating of the photosphere by continuum radiation from the chromosphere or possibly by high-energy protons, and iii) magnetic-force transients caused by magnetic reconnection. Hydrodynamic modeling of chromospheric shocks suggests that radiative losses deplete all but a small fraction of the energy initially deposited into them before they penetrate the photosphere. Comparisons between the spatial distribution of acoustic sources, derived from seismic holography of the surface signatures of flare acoustic emission, and the spatial distributions of sudden changes both in visible-light emission and in magnetic signatures offer a possible means of discriminating between contributions to flare acoustic emission from photospheric heating and magnetic-force transients. In this study we develop and test a means for estimating the seismic intensity and spatial distribution of flare acoustic emission from photospheric heating associated with visible-light emission and compare this with the helioseismic signatures of seismic emission. Similar techniques are applicable to transient magnetic signatures.     

On the convective instability of finitely conducting viscous model chromosphere

The convective stability of a simple model chromosphere, consisting of protons, electrons and hydrogen atoms in the ground state, has been studied in the presence of a vertically upward uniform magnetic field to include the effects of FLR, Hall-currents, finite conductivity and ionization. The ionization in the chromosphere is collisional and the recombination is radiative. It is found that the Schwarzchild criterion should necessarily be satisfied for the stability together with the condition thatv > 2v ₀, where is kinematic viscosity andv ₀ is gyroviscosity. Some special cases have also been investigated.     

A magnetodynamic mechanism for the heating of emerging magnetic flux tubes and loop flares

A new magnetodynamic model for loop flares is proposed to explain the following observational facts obtained from space during the last solar activity maximum: (i) Blueshifted lines of Ca xix and Fe xxv appear in some cases a minute or so before the initiation of impulsive bursts and relax into the unshifted lines with large width by the time of the onset of impulsive bursts, (ii) the hot source is formed by that time at the top of a loop-like structure, and confined there for a considerable time, and (iii) -ray line enhancement occurs at about the same time as hard X-ray spikes.In our model, the supply of energy to the loop top comes from below the chromosphere immediately before the flare (30 s-1 min before the hard X-ray impulsive bursts) in the form of the relaxing fronts of magnetic twist of opposite sign. These packets are thought to be built up in the process of loop emergence, stored at the footpoints of the loop below the photosphere, and released when the part of the feet floats up further. These released packets of magnetic twist drive the mass in the high chromosphere and transition zone into helical flows with pinch heating, and when these collide at the top of the loop, a very hot region appears there with a violent unwinding of the twists, resulting in the rapid dynamical annihilation of the magnetic energy, . Electrons and ions, raised to medium energies in the pinch at the incidence of the packets to the loop, are accelerated further by the Fermi-I mechanism between the approaching fronts of magnetic twist, and when B is weakened by unwinding they are released towards the chromosphere, and cause simultaneous -ray and hard X-ray bursts.     

Solar oscillations in strong and weak fraunhofer lines over a quiet region

We have analysed a 35-min-long time sequence of spectra in the Caii H line, Nai D1 and D2 lines, and in a large number of strong and weak Fei lines taken over a quiet region at the center of the solar disk. The time series of these spectra have been observed simultaneously in these lines under high spatial, spectral, and temporal resolution at the Vacuum Tower Telescope (VTT) of the Sacramento Peak Observatory. We have derived the line profiles and their central intensity values at the sites of the chromospheric bright points, which are visible in the H line for easy identification. We have done a power spectrum analysis for all the lines, using their central intensity values to determine the period of oscillations. It is shown that the 3 Fei lines, present 23 Å away from the core of the H line representing the pure photospheric lines, Nai D1 and D2 lines, 6 Fei lines at the wings of H line, and Can H line exhibit 5-min, 4.05-min, 3.96-min, and 3.2-min periodicity in their intensity oscillations, respectively. Since all these lines form at different heights in the solar atmosphere from low photosphere to middle chromosphere and show different periodicities in their intensity oscillations, these studies may give an idea about the spatial and temporal relation between the photospheric and chromospheric intensities. Therefore these studies will help to better understand the physical mechanisms of solar oscillations. It is clearly seen that the period of intensity oscillations decreases outward from the low photosphere to the middle chromosphere. Since we have studied a single feature at a time on the Sun (i.e., bright points seen in the H line) in all these spectral lines simultaneously, this may explain about the footpoints of the bright points, the origin of 3-min oscillations, and the relation to other oscillations pertaining to these locations on the Sun. We have concluded that 80% of the bright points are associated with dark elements in the true continuum, and they may seem to have a relationship with the dark intergranular lanes of the photosphere, after carefully examining the brightness (bright threads) extending from the core to the far wings of the H line at the locations of a large number of bright points, using their time sequence of spectra.NRC Resident Research Associate, on leave from Indian Institute of Astrophysics, Bangalore 560034, India.     

A simple approximation forJ v (σ v ,κ v ,τ v,B v ), and its application for temperature-correction procedures

For the case of isotropic coherent scattering plus absorption a simple expression is given (Equation 12) to compute the mean intensity of the radiationJ _v (as a function of optical depth _v ) if the scattering coefficient _v , the absorption coefficientK _v and the Planck functionB _v are given as a function of depth. In general the accuracy of this approximation is of the order of a few percent.A fairly simple temperature-correction procedure for the case when scattering is important is described.     

Force-free magnetic fields in the solar atmosphere

Reliable measurements of the solar magnetic field are restricted to the level of the photosphere. For about half a century attempts have been made to calculate the field in the layers above the photosphere, i.e. in the chromosphere and in the corona, from the measured photospheric field. The procedure is known as magnetic field extrapolation. In the superphotospheric parts of active regions the magnetic field is approximately force-free, i.e. electric currents are aligned with the magnetic field. The practical application to solar active regions has been largely confined to constant-α or linear force-free fields, with a spatially constant ratio, α, between the electric current and the magnetic field. We review results obtained from extrapolations with constant-α force-free fields, in particular on magnetic topologies favourable for flares and on magnetic and current helicities. Presently, different methods are being developed to calculate non-constant-α or nonlinear force-free fields from photospheric vector magnetograms. We also briefly discuss these methods and present a comparison of a linear and a nonlinear force-free magnetic field extrapolation applied to the same photospheric boundary data. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the occurrence of convective motions in the upper photosphere

We have examined whether the motion field in the photosphere in the range of optical depths 0.25< ₀< 0.6 is dominated by thermal convection or by vibrations. The observed asymmetries of infrared Fraunhofer lines indicate the presence of motions, and the fact that the asymmetry is zero for lines of low excitation and increases with the excitation potential shows that these motions are chiefly convective in this part of the photosphere: upward moving elements appear to be hotter than downward moving ones.Assuming furthermore that the photosphere can be described by a three-column model, with temperature differences as given by Edmonds (1967), we find that in the range of optical depths given above, where T seems to vary between 80 and 160 °K, average convective velocities of 2.3 to 3.2 km/sec should occur. This result is in numerical agreement with (a) a previous one by the present authors (1967) derived from the variation of line asymmetry with depth in lines of one multiplet, (b) a finding by Lambert and Mallia (1968) deduced from absolute wavelength measurements of Fraunhofer lines, and (c) a recent result of Beckers (1968) found from a comparison of two granulation pictures obtained simultaneously with a narrow-band filter centred on the two wings of a faint line.     

The sensitivity of the Ca II lines to the possible stratification of calcium abundance

Chemical composition anomalies of the corona and solar wind (FIP effect) have their origins in the chromosphere. We study the influence of the possible stratification of abundances in the chromosphere on the shape of corresponding strong Fraunhofer line profiles. A rough estimate is made of this influence for some variants of the variation in elemental abundance a _el(m) (m being the column mass). Calculations are made for the K and 8542Å Caii lines in two semi-empirical models of the chromosphere above sunspot umbra. In practice it is difficult to realize a direct search of the FIP effect in the chromosphere for lack of self-consistent theoretical models of the chromosphere required to match observations to them. An attempt is therefore proposed to find indirect signatures of the effect using a statistically significant sample of spectra for different objects, say, sunspots with a different magnetic structure and a different age.     

Granular-Scale Elementary Flux Emergence Episodes in a Solar Active Region

We analyse data from Hinode spacecraft taken over two 54-minute periods during the emergence of AR 11024. We focus on small-scale portions within the observed solar active region and discover the appearance of very distinctive small-scale and short-lived dark features in Ca ii H chromospheric filtergrams and Stokes I images. The features appear in regions with close-to-zero longitudinal magnetic field, and are observed to increase in length before they eventually disappear. Energy release in the low chromospheric line is detected while the dark features are fading. Three complete series of these events are detected with remarkably similar properties, i.e. lifetime of ≈ 12 min, maximum length and area of 2 – 4 Mm and 1.6 – 4 Mm², respectively, and all with associated brightenings. In time series of magnetograms a diverging bipolar configuration is observed accompanying the appearance of the dark features and the brightenings. The observed phenomena are explained as evidencing elementary flux emergence in the solar atmosphere, i.e. small-scale arch filament systems rising up from the photosphere to the lower chromosphere with a length scale of a few solar granules. Brightenings are explained as being the signatures of chromospheric heating triggered by reconnection of the rising loops (once they have reached chromospheric heights) with pre-existing magnetic fields, as well as being due to reconnection/cancellation events in U-loop segments of emerging serpentine fields. The characteristic length scale, area and lifetime of these elementary flux emergence events agree well with those of the serpentine field observed in emerging active regions. We study the temporal evolution and dynamics of the events and compare them with the emergence of magnetic loops detected in quiet Sun regions and serpentine flux emergence signatures in active regions. The physical processes of the emergence of granular-scale magnetic loops seem to be the same in the quiet Sun and active regions. The difference is the reduced chromospheric emission in the quiet Sun attributed to the fact that loops are emerging in a region of lower ambient magnetic field density, making interactions and reconnection less likely to occur. Incorporating the novel features of granular-scale flux emergence presented in this study, we advance the scenario for serpentine flux emergence.     

The production of lithium in the solar chromosphere and photosphere during white light flares

Recently new values of the lithium formation rate in low energy flares have been reported in the literature. These values are applied to the white light flare phenomenon on the Sun. It is found that the formation rate in the chromosphere is much larger than in the upper photosphere and that the ratio between the time integrated flare created abundance and the initial photospheric abundance is modest in the chromosphere and small in the upper photosphere. The yield of Li⁶ in the upper photosphere is, however, comparable to the upper limit of Li⁶ there.     

Characteristics of the Ca ii K-line profiles in the quiet Sun

Feature-to-feature identification is made on simultaneous Ca ii K-line spectrograms (SG) and K_2v spectroheliograms (SHG). The line profiles in plages and in the network boundary nearly always have double-peaked reversal in the core, while those inside the cells present all possibilities: double-peaked, single-peaked on violet side, single-peaked on red side, and unreversed absorption. Statistics of the profiles in the quiet chromosphere show that 50% are K₂₃₂ double-peaked, 20% are K_2v single-peaked, 10% are K_2r single-peaked, and 20% show only incipient reversal or even totally lack any reversal. We call attention to the nontrivial contribution of these absorption profiles which are formed in dark regions shown on SHG's.The physical conditions inferred from different kinds of profiles are briefly discussed.In part to be included in a dissertation to be submitted to the Graduate School, University of Maryland, by S. Y. Liu in partial fulfillment of the requirements for the Ph. D. degree in Astronomy.Visiting Graduate Student, Solar Division, Kitt Peak National Observatory.Visiting Astronomer, Solar Division, Kitt Peak National Observatory.     

The impulsive and gradual phases of a solar limb flare as observed from the solar maximum mission satellite

Simultaneous observations of a solar limb flare in the X-ray and ultraviolet regions of the spectrum are presented. Temporal and spectral X-ray observations were obtained for the 25–300 keV range while temporal, spectral, and spatial X-ray observations were obtained for the 30–0.3 keV range. The ultraviolet observations were images with a 10 spatial resolution in the lines of O v (T _e 2.5 × 10⁵ K) and Fe xxi (T _e 1.1 × 10⁷ K). The hard X-ray and O v data indicate that the impulsive phase began in the photosphere or chromosphere and continued for several minutes as material was ejected into the corona. Impulsive excitation was observed up to 30 000 km above the solar surface at specific points in the flare loop. The Fe xxi observations indicate a preheating before the impulsive phase and showed the formation of hot post-flare loops. This later formation was confirmed by soft X-ray observations. These observations provide limitations for current flare models and will provide the data needed for initial conditions in modeling the concurrent coronal transient.     

Photospheric and chromospheric activity on EY Dra

Magnetic activity in the photosphere and chromosphere of the M dwarf EY Dra is studied and possible correlations between the two are investigated using photometric observations in the V and R bands and optical and near infrared spectroscopy. The longitudinal spot configuration in the photosphere is obtained from the V band photometry, and the chromospheric structures are investigated using variations in the Hα line profile and observations of the Paschen β line. The shape of the V band light‐curve indicates two active regions on the stellar surface, about 0.4 in phase apart. The spectroscopic observations show enhanced Hα emission observed close to the phases of the photometrically detected starspots. This could indicate chromospheric plages associated with the photospheric starspots. Some indications of prominence structures are also seen. The chromospheric pressure is limited to log m_TR < –4 based on the non‐detection of emission in the Paschen β wavelength region. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of Wave Motions in the Active Region of the Solar Surface on Convection

The results of the observations of the active region (facula) near the center of the solar disk obtained with the German Vacuum Tower Telescope (VTT; Tenerife, Spain), are discussed. We have determined that the decrease in the contrast (brightness) of the facula with the magnetic field increasing from 130 to 160 mT is due to the fact that the V_V phase shift of waves in this range of magnetic field densities is close to zero (Φ _VV ≈ 0), i.e., the wave becomes stationary and does not transfer energy from the photosphere to the chromosphere. The sound waves that propagate from the chromosphere towards the photosphere significantly affect the temperature characteristics of turbulent vortices at the level of formation of the continuous spectrum. In particular, the contrast of granules under the influence of these waves can increase by 25%.     

Direct Propagation of Photospheric Acoustic <Emphasis Type="Italic">p</Emphasis> Modes into Nonmagnetic Solar Atmosphere

Solar p modes are one of the dominant types of coherent signals in Doppler velocity in the solar photosphere, with periods showing a power peak at five minutes. The propagation (or leakage) of these p-mode signals into the higher solar atmosphere is one of the key drivers of oscillatory motions in the higher solar chromosphere and corona. This paper examines numerically the direct propagation of acoustic waves driven harmonically at the photosphere, into the nonmagnetic solar atmosphere. Erdélyi et al. (Astron. Astrophys. 467, 1299, 2007) investigated the acoustic response to a single point-source driver. In the follow-up work here we generalise this previous study to more structured, coherent, photospheric drivers mimicking solar global oscillations. When our atmosphere is driven with a pair of point drivers separated in space, reflection at the transition region causes cavity oscillations in the lower chromosphere, and amplification and cavity resonance of waves at the transition region generate strong surface oscillations. When driven with a widely horizontally coherent velocity signal, cavity modes are caused in the chromosphere, surface waves occur at the transition region, and fine structures are generated extending from a dynamic transition region into the lower corona, even in the absence of a magnetic field.     

Dynamical fine structure of a quiescent filament

A quiescent filament was observed near the center of the disk (N5, W5) with the MSDP spectrograph of the 50 cm refractor of the Pic-du-Midi Observatory on June 17, 1986. We focus our study on the statistical moments of the Dopplershift,V ₁, and the intensity,I ₁, at the center of a chord of the Hα profile (±0.256 Å), versus the minimum intensityI ₀. We use a statistical model simulating a numbern _max of threads (of optical thicknessτ ₀ and source functionS ₀), seen over the chromosphere. The threadsj along the same line-of-sighti are identical except for the velocityv _j (gaussian distributionv ₀,σ _v). We search for the best fit between the observed and simulated quantities:V ₁,σ (V ₁),I ₁,σ(I ₁), and the histogram of theI ₀ values over the field of view. A good fit is obtained with: (a) threads characterized byτ ₀ = 0.2,S ₀ = 0.06 (unit of the continuum at disk center), mean upward velocityv ₀ = 1.7 km s⁻¹ and gaussian-type velocity distributionσ _v = 3.5 km s⁻¹. Other possible values ofτ ₀ andσ _v are discussed; (b) underlying chromosphere deduced from observed quiet Sun (outside the filament) by modifying the chromospheric velocities: additional mean upward velocity 0.7 km s⁻¹, standard deviation reduced by a factorF _c ∼ 0.7. The results are discussed in connection with the values deduced from prominence observations.     

Response of an optically thin,isothermal atmosphere to a convective overshoot

Radiation is believed to be hostile to the generation of gravity waves by granulation at the base of photosphere where the radiation is effective. A convective overshoot from subphotosphere seems able to penetrate to a height where the solar temperature is minimum and to excite the gravity waves in a stable region there.The response of the solar atmosphere to a Gaussian disturbance characterizing such a convective overshoot is studied in an unbounded isothermal atmosphere. Radiative effects are included, but only in regions which are optically thin. The response is measured in terms of mean vertical kinetic energy density (E _z) and mean vertical external energy flux (Q _z). E _z and Q _z were calculated for a wide range of frequencies centered at the observed 5-min velocity oscillation period. The computed sharp and broad power spectra at the lower chromosphere and the upper photosphere, respectively, are attributed to the combined effects of space damping and source function. Low-frequency waves (2000 s or longer) are found to be not responsible for depositing energy in the upper solar atmosphere.The National Center for Atmospheric Research is sponsored by the National Science Foundation.