Search

During the last three years we have continued the development of extensive computer programs for constructing realistic models of the solar atmosphere and for calculating detailed spectra to use in the interpretation of solar observations. This research involves two major interrelated efforts: work by Avrett and Loeser on the Pandora computer program for optically thick non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed high-resolution synthesis of the solar spectrum using data for over 58 million atomic and molecular lines. Our objective is to construct atmospheric models from which the calculated spectra agree as well as possible with high-and low-resolution observations over a wide wavelength range. Such modeling leads to an improved understanding of the physical processes responsible for the structure and behavior of the atmosphere.

This grant funded an observational and theoretical program to study the structure and dynamics of the solar photosphere and low chromosphere, and the spectral signatures that result. The overall goal is to learn about mechanisms that cause heating of the overlying atmosphere, and produce variability of solar emission in spectral regions important for astrophysics and space physics. The program exploited two new ground-based observational capabilities: one using the Swedish Solar Telescope on La Palma for very high angular resolution observations of the photospheric intensity field (granulation) and proxies of the magnetic field (G-band images); and the other using the Near Infrared Magnetograph at the McMath-Pierce Solar Facility to map the spatial variation and dynamic behavior of the solar temperature minimum region using infrared CO lines. We have interpreted these data using a variety of theoretical and modelling approaches, some developed especially for this project. Previous annual reports cover the work done up to 31 May 1997. This final report summarizes our work for the entire period, including the period of no-cost extension from 1 June 1997 through September 30 1997. In Section 2 we discuss observations and modelling of the photospheric flowfields and their consequences for heating of the overlying atmosphere, and in Section 3 we discuss imaging spectroscopy of the CO lines at 4.67 mu.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Rurucz on the detailed synthesis of the solar spectrum based on opacity data or over 58 million atomic and molecular lines. our goals are: to determine models of the various features observed on the Sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for non-radiative heating, and for solar activity in general.

One-dimensional hydrostatic models of quiet and active solar regions can be constructed that generally account for the observed intensities of lines and continue throughout the spectrum, except for the infrared CO lines. There is an apparent conflict between: (1) observations of the strongest infrared CO lines formed in LTE at low-chromospheric heights but at temperatures much cooler than the average chromospheric values; and (2) observations of Ca II, UV (ultraviolet), and microwave intensities that originate from the same chromospheric heights but at the much higher temperatures characteristic of the average chromosphere. A model M(sub CO) has been constructed which gives a good fit to the full range of mean CO line profiles (averaged over the central area of the solar disk and over time) but this model conflicts with other observations of average quiet regions. A model L(sub CO) which is approximately 100 K cooler than M(sub CO) combined with a very bright network model F in the proportions 0.6 L(sub CO) + 0.4 F is found to be generally consistent with the CO, Ca II, UV, and microwave observations. Ayres, Testerman, and Brault found that models COOLC and FLUXT in the proportions 0.925 and 0.075 account for the CO and Ca II lines, but these combined models give an average UV intensity at 140 nm about 20 times larger than observed. The 0.6 L(sub CO) + 0.4 F result may give a better description of the cool and hot components that produce the space- and time-averaged spectra. Recent observations carried out by Uitenbroek, Noyes, and Rabine with high spatial and temporal resolution indicate that the faintest intensities in the strong CO lines measured at given locations usually become much brighter within 1 to 3 minutes. The cool regions thus seem to be mostly the low-temperature portions of oscillatory waves rather than cool structures that are stationary.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed synthesis of the solar spectrum based on opacity data for over 58 million atomic and molecular lines. Our goals are to determine models of the various features observed on the sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for nonradiative heating, and for solar activity in general.

The resonance lines of singly ionized magnesium, the MgII h&k lines at about 280 nm, are two of the small number of lines in the solar spectrum that are optically thick in the chromospheric part of the solar atmosphere. Potentially these lines contain information on the initial temperature rise that occurs at the top of the photosphere. Unfortunately, few good observations of the lines exist due to their wavelength near 280 nm the ultraviolet. However, a fair number of observations (on the order of 200) are available from the data base of the UltraViolet Polarimeter and Spectrometer (UVSP) instrument that flew on board of NASA's Solar Maximum Mission (SMM) satellite. In addition, this data base contains a number of spectra that include the Mg I resonance line at (lambda)285.2nm, just longward of the h&k lines. The neutral magnesium line is not as strong as its ionic counterparts and samples slightly lower parts of the atmosphere. Its width is a sensitive diagnostic of the ionization balance between neutral and singly ionized magnesium, which determines the opacity scale (and formation height) of other diagnostically important MgI lines like the 457.1 nm intercombination line, the magnesium b lines and the infrared MgI emission lines near 12 microns. Analysis of the observed line profiles shows that it is necessary to include the effects of partial frequency redistribution (PRD) in the formation of the line as in the case of the h&k lines. This implies that the core of the line is very sensitive to the way scattering is treated in the modeling of the line, and in turn this allows us to separate the uncertain effects in the atomic data (viz. the Van der Waals broadening) from the uncertainties in the underlying atmospheric model. The main objective of this research was to compare observed spectra of the magnesium resonance lines against theoretical line profiles calculated from recent models of the solar atmosphere by Fontenla et al., hereafter called FAL. These models extend earlier solar models by Avrett et al. and provide yet a better match between observed and calculated spectra, in particular of the hydrogen Lyman line, by including the effects of ambipolar diffusion of hydrogen and helium in the mass and energy transport equations for the chromosphere-corona transition region. Although the transition region models can be built self-consistently (without additional observational input) once the proper boundary conditions are given at the top of the chromosphere, there are uncertainties in the models precisely because these boundary conditions are not well known. With the analysis of UVSP magnesium observations we strive to constrain these uncertainties.

Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Kurucz on the detailed synthesis of the solar spectrum based on opacity data for over 58 million atomic and molecular lines. Our goals are to determine models of the various features observed on the Sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for non-radiative heating, and for solar activity in general.

Our energy-balance models of the lower transit ion region were presented in a previous paper (New Models of the Chromosphere and Transition Region). Here we show the influence of particle diffusion on the calculated hydrogen and helium number densities for a given temperature-density model (model C in the preceding paper). We have also solved the statistical equilibrium and radiative transfer equations for a 13-level He I atom (22 radiative transitions) and a 6-level He II ion (15 radiative transitions) together with He III. The resulting He I and He II level-1 number densities and He III density are shown as functions of temperature. Diffusion substantially increases eta(sub HeI) for T greater than 35,000 K and decreases eta(sub HeI) (while increasing eta(sub(HeII)) between 9,000 and 25,000 K. Including the effects of diffusion also increases eta(sub HeIII) for T less than 60,000 K.

The present characterization of the solar photosphere as a radiative boundary for the solar interior gives attention to both semiempirical and a priori models for the structure of emergent radiation, and emphasizes that any discussion of photospheric structures smaller than 300 km is easily compromised by the radiant integration over such scale lengths. Plane-parallel modes have indicated that Boltzmann-Saha statistics are largely valid over the photosphere save within 100 km of the temperature minimum. An evaluation is made of observational evidence for horizontal inhomogeneities and such characteristics as granular and temperature minimum region fluctuations, as well as the recent demonstration of relatively cold CO-absorbing regions.

Models are presented for the two continuum-emitting kernels observed in the white-light flare (WLF) of June 15, 1982. They are the first semiempirical models of a WLF which are consistent not only with observations of the continuum emission level but also with a set of spectral lines having heights of formation which span the chromosphere and upper photosphere. It is shown that the models are not compatible with the hypothesis that the continuum emission is caused by enhanced Balmer and Paschen hydrogen continua, and they present strong evidence instead that the emission is of photospheric origin and that its source is due to H(-). The observation and the models derived from them show that white-light emission can occur in areas of the active region where there is no chromospheric emission and in particular no H-alpha emission. This fact seems to rule out the viability of downward transport mechanisms as the source of the energy required for the WLF.

The CO cooling rate for models of the solar atmosphere using the detailed line-by-line CO opacity in the fundamental band, and carrying out a full radiative transfer calculation for each line is computed. The importance of the different assumptions that have been made to obtain the CO cooling rate and find that when detailed optical depth effects are taken into account, the calculated CO cooling rate at line optical depths near unity can be smaller than optically thin estimates by more than an order of magnitude is studied. It is found that CO cooling does not account for the missing source of radiative cooling in the temperature minimum region of the quiet sun.

Various methods of inverting solar Mg I 4571 and 5173 spectral line observations are examined to find the best method of using these lines to calculate the vertical temperature and electron density structure around the temperature minimum region. Following a perturbation analysis by Mein (1971), a Fredholm integral equation of the first kind is obtained which can be inverted to yield these temperature and density structures as a function of time. Several inversion methods are tested and compared. The methods are used to test data as well as to a subset of observations of these absorption lines taken on February 3, 1986 before and during a solar flare. A small but significant increase is found in the temperature and a relatively large increase in the electron density during this flare. The observations are inconsistent with heating and ionization by an intense beam of electrons and with ionization by UV photoionization of Si I.

A nine-level atomic model is presented for C I lines synthesis in the sun. All the atomic parameters are based on recent experimental and theoretical data. Profiles of the multiplets at 1561 and 1657 A are computed from a model of the solar chromosphere. The sensitivity of the results to changes in the atomic parameters is examined. Partial frequency redistribution must be included in these lines to obtain agreement between calculated and observed profiles.

Consideration is given to empirical atmospheric structures that are consistent with enhanced white-light continuum emission in solar flares. Results are presented from calculations of radiative transfer in lines and continua in empirical white-light flare model atmospheres, showing that flares with strong emission in the Balmer lines and continuum must show increases at longer wavelengths due to H(-) emission from overheated photospheric levels, which the Paschen continuum contribution in the same wavelength range is neglible. Also, plausible heating mechanisms that can lead to white-light flare emission are examined.

This paper presents a complete atomic model for Mg I line synthesis, where all the atomic parameters are based on recent experimental and theoretical data. It is shown how the computed profiles at 4571 A and 5173 A are influenced by the choice of these parameters and the number of levels included in the model atom. In addition, observed profiles of the 5173 A b2 line and theoretical profiles for comparison (based on a recent atmospheric model for the average quiet sun) are presented.