Modeling the (upper) solar atmosphere including the magnetic field

The atmosphere of the Sun is highly structured and dynamic in nature. From the photosphere and chromosphere into the transition region and the corona plasma-β changes from above to below one, i.e., while in the lower atmosphere the energy density of the plasma dominates, in the upper atmosphere the magnetic field plays the governing role – one might speak of a “magnetic transition”. Therefore the dynamics of the overshooting convection in the photosphere, the granulation, is shuffling the magnetic field around in the photosphere. This leads not only to a (re-)structuring of the magnetic field in the upper atmosphere, but induces also the dynamic reaction of the coronal plasma, e.g., due to reconnection events. Therefore the (complex) structure and the interaction of various magnetic patches is crucial to understand the structure, dynamics and heating of coronal plasma as well as its acceleration into the solar wind.

The present article will emphasize the need for three-dimensional modeling accounting for the complexity of the solar atmosphere to understand these processes. Some advances on 3D modeling of the upper solar atmosphere in magnetically closed as well as open regions will be presented together with diagnostic tools to compare these models to observations. This highlights the recent success of these models which in many respects closely match the observations.     

Solar cosmic rays in the near Earth space and the atmosphere

Solar energetic particles (SEPs) constitute a distinct population of energetic charged particles, which can be often observed in the near Earth space. SEP penetration into the Earth’s magnetosphere is a complicated process depending on particle magnetic rigidity and geomagnetic field structure. Particles in the several MeV energy range can only access to periphery of the magnetosphere and the polar cap regions, while the GeV particles can arrive at equatorial latitudes. Solar protons with energies higher than 100 MeV may be observed in the atmosphere above ∼30 km, and those with energies more than 1 GeV may be recorded even at the sea level. There are some observational evidences of SEP influence on atmospheric processes. Intruding into the atmosphere, SEPs affect middle atmosphere odd-nitrogen and ozone chemistry. Since spatial and temporal variations of SEP fluxes in the near Earth space are controlled by solar activity, SEPs may present an important link between solar activity and climate. The paper outlines dynamics of SEP fluxes in the near Earth space during the last decades. This can be useful for tracing relationship between SEPs and atmospheric processes.     

太阳大气中磁场的扭转储能

本文从完整的磁流体动力学方程组出发, 研究太阳大气中磁力线管根部的扭转储存能量。通过数值方法, 研究了包含太阳表面过渡区在内的非线性无力场的动力学演化。假设初始磁场位形为势场, 根部等离子体运动使活动区磁场扭转。磁能蓬新分布, 在局部区域中储存大量的磁能。计算结果给出非线性无力场的定量关系, 可以解释太阳耀斑的储能过程。      

Dependence of sporadic-E layer and lower thermosphere dynamics on solar activity

The influence of solar activity level on frequency and stochastic sporadic-E layer parameters is considered for station Moscow (56 N, 37 E) during three solar cycles (1958–1988). Es layers with various intensity levels and times of appearance were considered separately. Some regularities of sporadic-E layer reaction to variations of solar activity are revealed. An explanation of the revealed regularities is offered in terms of dynamic processes of the neutral atmosphere.     

The thermosphere as a sink of magnetospheric energy: A review of recent observations of dynamics

The study of the dynamics and thermodynamics of the earth's upper atmosphere has made significant progress over the past few years owing to the availability of new global-scale data sets from the Dynamics Explorer satellites. The thermospheric wind and temperature fields at high latitude have been observed to depend strongly on forcing processes of magnetospheric origin. A key momentum source is due to the drag effect of ions convecting in response to electric fields mapped down on the ionosphere from magnetospheric boundary regions. Likewise, an important heat source derives from Joule or frictional dissipation due to ion/neutral difference velocities governed, in turn, by magnetospheric forcing. In this paper we discuss the progress made over the last 2–3 years initiated by the new satellite measurements and we review published data on ion and neutral motions in the context of the energy and momentum coupling between the magnetosphere and the ionosphere/neutral upper atmosphere. The observations indicate the existence of a “flywheel effect” which implies direct feedback from the neutral thermosphere to the magnetosphere via the release of energy and momentum previously “stored” in the neutral thermosphere.     

Relationship of solar flare accelerated particles to solar energetic particles (SEPs) observed in the interplanetary medium

Observations of hard X-ray (HXR)/γ-ray continuum and γ-ray lines produced by energetic electrons and ions, respectively, colliding with the solar atmosphere, have shown that large solar flares can accelerate ions up to many GeV and electrons up to hundreds of MeV. Solar energetic particles (SEPs) are observed by spacecraft near 1 AU and by ground-based instrumentation to extend up to similar energies as in large SEP events, but it appears that a different acceleration process, one associated with fast coronal mass ejections is responsible. Much weaker SEP events are observed that are generally rich in electrons, ³He, and heavy elements. The energetic particles in these events appear to be similar to those accelerated in flares. The Ramaty high energy solar spectroscopic imager (RHESSI) mission provides high-resolution spectroscopy and imaging of flare HXRs and γ-rays. Such observations can provide information on the location, energy spectra, and composition of the flare accelerated energetic particles at the Sun. Here, preliminary comparisons of the RHESSI observations with observations of both energetic electron and ion near 1 AU are reviewed, and the implications for the particle acceleration and escape processes are discussed.     

The dynamics of the solar radiative zone

The picture of the solar radiative zone is evolving quickly. This review is separated in two parts. We first recall how the two powerful probes of the solar interior, namely the neutrinos and helioseismology have scrutinized the microscopic properties of the solar radiative plasma. Recent observations stimulate today complementary activities beyond the standard stellar model through theoretical modeling of angular momentum transport by rotation, internal waves or (and) by magnetic fields to get access to the dynamical motions of this important region of the Sun. So in the second part, we summarize the first impact of such processes on the radiative zone.     

Probing heliospheric diffusion coefficients with solar energetic particles

The dynamics of solar particle events provide a direct link to the understanding of properties of wave–particle interactions, and to the nature of the solar wind fluctuations. Depending on their energy, the often simultaneously observed electrons, protons and ions interact with different wavenumber ranges of the fluctuations, and are sensitive to various aspects of the dynamical nature of the solar wind turbulence. In general, the evolution of particle events is also sensitive to the spatial variation of the transport parameters between the Sun and a few AU. Together with in situ plasma and magnetic field observations this information can be used to extrapolate the properties of transport parameters into the more distant Heliosphere. Recent developments in the theory of parallel transport of energetic particles, and examples for the modelling of solar particle events and the derivation of transport parameters are considered. A dynamical quasi-linear theory is presented which gives special emphasis to the geometry and dynamic nature of the fluctuations, and which is able to provide particle mean free paths solely from observed plasma parameters, in good agreement with those derived by the modelling. Possibilities to apply the above results to the study of other energetic particle processes in the Heliosphere are discussed.     

On the evolution of dust in the solar vicinity.

The analysis of interplanetary dust shows that the majority of particles in out-of-ecliptic regions comes from comets and also that near solar dust, in the ecliptic regions, results most probably largely from comets. The intense radiation flux in the solar vicinity is expected to cause strong modifications in the material composition and surface structure of interplanetary dust particles and hence the analysis of near solar dust provides interesting insights into the evolution of meteoritic, especially cometary materials. Because of the lack of in-situ measurements our present knowledge concerning these processes derives from remote sensing, i.e. observations of the solar F-corona. In particular these are observations of albedo, polarization and colour temperature given in terms of average particle properties. For example the analysis of near infra-red F-corona data points to the existence of a strong component of irregularly structured silicate particles, most probably of cometary origin. The data may indicate a subsequent sublimation of different particles or different constituents of the particles. Here we compare particle properties derived from F-corona observations with model calculations of single particle properties and discuss perspectives of future analysis of cometary dust in the interplanetary cloud.     

Inner radiation belt source of helium and heavy hydrogen isotopes

Nuclear interactions between inner zone protons and atoms in the upper atmosphere provide the main source of energetic H and He isotopes nuclei in the radiation belt. This paper reports on the specified calculations of these isotope intensities using various inner zone proton intensity models (AP-8 and SAMPEX/PET PSB97), the atmosphere drift-averaged composition and density model MSIS-90, and cross-sections of the interaction processes from the GNASH nuclear model code. To calculate drift-averaged densities and energy losses of secondaries, the particles were tracked in the geomagnetic field (modelled through IGRF-95) by integrating numerically the equation of the motion. The calculations take into account the kinematics of nuclear interactions along the whole trajectory of trapped proton. The comparison with new data obtained from the experiments on board RESURS-04 and MITA satellites and with data from SAMPEX and CRRES satellites taken during different phases of solar activity shows that the upper atmosphere is a sufficient source for inner zone helium and heavy hydrogen isotopes. The calculation results are energy spectra and angular distributions of light nuclear isotopes in the inner radiation belt that may be used to develop helium inner radiation belt model and to evaluate their contribution to SEU (single event upset) rates.     

Middle atmosphere electrical structure,dynamics and coupling

The ram current to ion traps and the insensitivity of ion conductivity to compressibility provide the basis of robust techniques for middle atmosphere measurements. Gerdien condensers are more difficult to implement but provide more information. Mesospheric electrical conductivity shows many orders of magnitude variability, with depressions below gas phase model values indicating dominance by aerosol particles. The mobility of these ions has been directly measured and indicates particles of thousands of AMU. Large mesospheric fields have come into question, and diagnostic measurements show that many such measurements may be artifacts. However, some measurements of V/m fields with symmetrical and redundant sensors appear to be real. These fields complicate the “mapping” picture of electrical coupling and may also modulate the transport of aerosol particles. They are probably related to neutral atmosphere dynamics and/or the aerosol particles. Lightning couples much more energy to the middle atmosphere and above than previously suspected, primarily in the ELF-ULF range. There are many important unanswered questions in this relatively unexplored frontier area which may be answered with low cost balloon and sounding rocket experiments.     

Understanding solar flares from opticalobservations: How do particle beams affect the lower atmosphere?

During the impulsive phase of solar flares, both hard X-ray (HXR) and optical emissions exhibit fast temporal fluctuations detectable down to sub-second scales. This is usually ascribed to the propagation of beams of accelerated particles and to the dissipation of their energy in lower layers of the solar atmosphere. Although it is rather difficult to prove a temporal correlation between HXR and optical intensity variations, we discuss here some previous results and recent attempts. Namely in coordination with RHESSI observations, several ground-based observatories started to detect fast optical variations in the H line. In addition to this, we also mention a possibility of using some other diagnostically important lines. The proper interpretation of coordinated HXR and optical observations further requires robust tools for radiation-hydrodynamical (RHD) forward modeling. We briefly describe a new ‘hybrid’ code which consists of RHD part and particle-simulation part. Short-duration heating due to beam pulses is modeled which allows us to predict temporal fluctuations of HXR and selected optical and UV lines formed in chromospheric layers and in the transition region. Particularly the line asymmetries originating in a highly dynamical lower atmosphere of the flare can be used to diagnose the response of these layers to particle beams.     

Viewing radiation signatures of solar energetic particles in interplanetary space

A current serious limitation on the studies of solar energetic particle (SEP) events is that their properties in the inner heliosphere are studied only through in situ spacecraft observations. Our understanding of spatial distributions and temporal variations of SEP events has come through statistical studies of many such events over several solar cycles. In contrast, flare SEPs in the solar corona can be imaged through their radiative and collisional interactions with solar fields and particles. We suggest that the heliospheric SEPs may also interact with heliospheric particles and fields to produce signatures which can be remotely observed and imaged. A challenge with any such candidate signature is to separate it from that of flare SEPs. The optimum case for imaging high-energy (E > 100 MeV) heliospheric protons may be the emission of π⁰-decay γ-rays following proton collisions with solar wind (SW) ions. In the case of E > 1 MeV electrons, gyrosynchrotron radio emission may be the most readily detectible remote signal. In both cases we may already have observed one or two such events. Another radiative signature from nonthermal particles may be resonant transition radiation, which has likely already been observed from solar flare electrons. We discuss energetic neutrons as another possible remote signature, but we rule out γ-ray line and 0.511 MeV positron annihilation emission as observable signatures of heliospheric energetic ions. We are already acquiring global signatures of large inner-heliospheric SW density features and of heliosheath interactions between the SW and interstellar neutral ions. By finding an appropriate observable signature of remote heliospheric SEPs, we could supplement the in situ observations with global maps of energetic SEP events to provide a comprehensive view of SEP events.     

Combining MHD and kinetic modelling of solar flares

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport in converging magnetic fields.     

Calculation of X-ray polarization during flares

The polarization of X-radiation emitted by electrons which are accelerated during flares is investigated in a non-thermal model. Monte Carlo calculations combining analytically treated multiple-scattering and random large-angle scattering processes yield the energy and angular distributions of flare electrons penetrating the upper solar atmosphere. The X-ray polarization derived from these distributions is computed as a function of energy, observation angle and column density traversed by the electrons. Initially monoenergetic electrons as well as electrons with a power-law spectrum are considered. The degree of polarization of the total X-radiation obtained by summing over all layers of the atmosphere is compared with observations.     

Towards multi-dimensional space weather modeling for energetic oxygen ions in the earth''s inner magnetosphere: Equilibrium configuration

We report the first 3+1 dimensional model development for energetic atomic oxygen ions in the Earth's radiation belts. Energetic Oxygen ions cans be supplied to the Earth's Inner magnetosphere from the sun (as a component of solar wind and solar energetic particles), from anomalous cosmic rays, and from acceleration processes acting on ionospheric atomic oxygen ions. We have built a multi-dimensional oxygen ion model in the following free parameters: geomagnetic L-shell, the magnetic moment, the second adiabatic invariant, and the discrete charge state number. Quiet time, steady state oxygen ion distributions have been obtained numerically from an assumed outer radiation zone boundary condition at L=7, average values of the radial diffusion coefficients, and standard values for the exospheric neutral densities due to the MSIS-86 upper atmosphere and exosphere neutral thermal particle density model. Average distributions of free electrons in the plasmasphere were also assumed with a mean plasmapause location just beyond L=4. We included the six lowest ionic charge states of atomic oxygen (¹⁶O) based on an existing charge exchange cross section compilation by Spjeldvik and Fritz (1978). Computed oxygen ion distributions include the resulting equilibrium structure of energy oxygen ions between 10 KeV and 100 MeV.     

Photochemical response of neutral and ionized middle atmosphere composition to the strong solar proton event of October 1989

One of the strongest solar proton events (SPE) occurred in October 1989. Its forcing of the middle atmosphere chemistry including ionized components in the D-region is examined. The ionization rate, and ozone, NO and OH density temporal and spatial (vertical) deviations induced by the SPE, calculated by a 1-D time-dependent photochemical model separately for daytime and nighttime (not shown here), are used in a 1-D model of the lower ionosphere to calculate the response of ionized components to combined forcing by ionization rate and neutral chemical composition disturbances. The radio wave absorption caused by electron density disturbances after the SPE is calculated and compared with observations. The computed ozone values are compared with observations, as well.     

Modeling ground and space based cosmic ray observations

After entering our local astrosphere (called the heliosphere), galactic cosmic rays, as charged particles, are affected by the Sun’s turbulent magnetic field. This causes their intensities to decrease towards the inner heliosphere, a process referred to as modulation. Over the years, cosmic ray modulation has been studied extensively at Earth, utilizing both ground and space based observations. Moreover, modelling cosmic ray modulation and comparing results with observations, insight can be gained into the transport of these particles, as well as offering explanations for observed features. We review some of the most prominent cosmic ray observations made near Earth, how these observations can be modelled and what main insights are gained from this modelling approach. Furthermore, a discussion on drifts, as one of the main modulation processes, are given as well as how drift effects manifest in near Earth observations. We conclude by discussing the contemporary challenges, fuelled by observations, which are presently being investigated. A main challenge is explaining observations made during the past unusual solar minimum.     

Atmospheric dynamics on Venus and Mars

New results from Pioneer Orbiter observations indicate a continued vortex organization of the cloud level atmosphere in either hemisphere, centered over respective poles. Significant changes in the magnitude of the cloud level zonal circulation over a period of several years have been detected. A strong signature of the solar tidal circulation has been detected in the atmospheric circulation with the lowest speeds occurring in equatorial latitudes about 20° upstream of the sub-solar point. Finally, a solar-locked persistent spatial structure has been discovered in the variance of the ultraviolet brightness measured from brightness normalized images of Venus. Vega balloons (drifting at about 53 km altitude near 7°N and 7°S latitudes) have also provided some unique observations of atmospheric circulation, significant among them being the strong vertical motions, the zonality of their drift speeds as well as a significant temperature difference between the two balloons. The temperature difference which amounts to 6.5°K on average is currently being interpreted as a temperature variation with longitude or time.

Diagnostic modelling efforts towards simulating the atmospheric circulation on Venus are continuing and have provided some clues about the processes that maintain them but have not yet been successful in explaining the superrotation of the atmosphere.

Knowledge of the Martian atmospheric dynamics on the other hand is still limited by lack of adequate observations. Numerical modelling of the Martian atmosphere continues to provide most of the information about the atmospheric circulation. The situation regarding the paucity of observations should improve with the completion of the proposed Mars Observer mission. The low circular polar orbit planned provides an excellent opportunity to study the Martian atmosphere.