Non-linear damping of visco-resistive Alfvén waves in solar spicules

Interaction of Alfvén waves with plasma inhomogeneities generates phase mixing which can lead to dissipate Alfvén waves and to heat the solar plasma. Here we study the dissipation of Alfvén waves by phase mixing due to viscosity and resistivity variations with height. We also consider nonlinear magnetohydrodynamic (MHD) equations in our theoretical model. Non-linear terms of MHD equations include perturbed velocity, magnetic field, and density. To investigate the damping of Alfvén waves in a stratified atmosphere of solar spicules, we solve the non-linear MHD equations in the x–z plane. Our simulations show that the damping is enhanced due to viscosity and resistivity gradients. Moreover, energy variations is influenced due to nonlinear terms in MHD equations.     

Reflection Properties of Gravito-MHD Waves in an Inhomogeneous Horizontal Magnetic Field

We derive the dispersion equation for gravito-magnetohydrodynamical (MHD) waves in an isothermal, gravitationally stratified plasma with a horizontal inhomogeneous magnetic field. Sound and Alfvén speeds are constant. Under these conditions, it is possible to derive analytically the equations for gravito-MHD waves. The high values of the viscous and magnetic Reynolds numbers in the solar atmosphere imply that the dissipative terms in the MHD equations are negligible, except in layers around the positions where the frequency of the MHD wave equals the local Alfvén or slow wave frequency. Outside these layers the MHD waves are accurately described by the equations of ideal MHD. We consider waves that propagate energy upward in the atmosphere. For the plane boundary, z=0, between two isothermal plasma regions with horizontal but different magnetic fields, we discuss the boundary conditions and derive the equations for the reflection and transmission coefficients. In the simpler case of a gravitationally stratified plasma without magnetic field, these coefficients describe the reflection and transmission properties of gravito-acoustic waves.     

The solar corona as a minimum energy system?

This paper is an exploration of the possibility that the large-scale equilibrium of plasma and magnetic fields in the solar corona is a minimum energy state. Support for this conjecture is sought by considering the simplest form of that equilibrium in a dipole solar field, as suggested by the observed structure of the corona at times of minimum solar activity. Approximate, axisymmetric solutions to the MHD equations are constructed to include both a magnetically closed, hydrostatic region and a magnetically open region where plasma flows along field lines in the form of a transonic, thermally-driven wind. Sequences of such solutions are obtained for various degrees of magnetic field opening, and the total energy of each solution is computed, including contributions from both the plasma and magnetic field. It is shown that along a sequence of increasingly closed coronal magnetic field, the total energy curve is a non-monotonic function of the parameter measuring the degree of magnetic field opening, with a minimum occurring at moderate field opening.For reasonable choices of model parameters (coronal temperature, base density, base magnetic field strength, etc.), the morphology of the minimum energy solution resembles the observed quiet, solar minimum corona. The exact location energy minimum along a given sequence depends rather sensitively on some of the adopted parameter values. It is nevertheless argued that the existence of an energy minimum along the sequences of solutions should remain a robust property of more realistic coronal wind models that incorporate the basic characteristics of the equilibrium corona- the presence of both open and closed magnetic regions.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Solar coronal heating by magnetosonic waves

Solar coronal heating by magnetohydrodynamic (MHD) waves is investigated. ultraviolet (UV) and X-ray emission lines of the corona show non-thermal broadenings. The wave rms velocities inferred from these observations are of the order of 25–60 km s⁻¹ . Assuming that these values are not negligible, we solved MHD equations in a quasi-linear approximation, by retaining the lowest order non-linear term in rms velocity. Plasma density distribution in the solar corona is assumed to be inhomogeneous. This plasma is also assumed to be permeated by dipole-like magnetic loops. Wave propagation is considered along the magnetic field lines. As dissipative processes, only the viscosity and parallel (to the local magnetic field lines) heat conduction are assumed to be important. Two wave modes emerged from the solution of the dispersion relation. The fast mode magneto-acoustic wave, if originated from the coronal base can propagate upwards into the corona and dissipate its mechanical energy as heat. The damping length-scale of the fast mode is of the order of 500 km. The wave energy flux associated with these waves turned out to be of the order of 2.5×10⁵ ergs cm⁻² s⁻¹ which is high enough to replace the energy lost by thermal conduction to the transition region and by optically thin coronal emission. The fast magneto-acoustic waves prove to be a likely candidate to heat the solar corona. The slow mode is absent, in other words cannot propagate in the solar corona.     

灾变理论在太阳物理和天体物理其他研究领域中的应用

叙述和介绍了太阳爆发的磁通量绳灾变理论和模型的发展过程,强调了建立这样的模型所需要的观测基础。讨论了由模型所预言的爆发磁结构的几个重要特征以及观测结果对这种预言的证实。在此模型的基础上,讨论了一个典型的爆发过程中所出现的不同现象及它们之间的相互关系。最后,介绍了作者的一项最新尝试:将太阳爆发的灾变理论和模型应用到对黑洞吸积盘间歇性喷流的理论研究当中,以及研究所取得的初步结果。     

太阳异极性磁区的相互入侵

我们在文[1]的启发下,计算了磁中性线附近异极性磁区相互入侵(或挤压)引起的等离子体动力学问题。气体初态取用流行的宁静太阳光球色球大气模型,即非等温的密度指数变化的重力分层大气。采用Lagrangian格式数值求解自洽的MHD方程,这可使入侵力学变得直观明显——磁场随流体而运动。我们的新结果是入侵流动在光球低层产生出强的水平磁场(即强的横向场),但光球高层和色球低层的磁结构却变化不大,有力地支持了文[13]提出的光球色球里可能出现磁流体力学间断面的概念。入侵确实在磁中性线附近建立了电流片,但这电流片主要在光球低层,其量级和观测一致。另外还显示垂直下降运动也可能导致异极磁区的入侵。尽管在MHD~1方程里包含了电阻耗散和热传导流,但计算证明它们对入侵力学影响不大,热传导的作用只是使气体温度分布逐渐趋于宁静太阳分布(尽管高度变了)。     

Linear and non-linear MHD wave propagation in steady-state magnetic cylinders

The propagation of linear and non-linear magnetohydrodynamic (MHD) waves in a straight homogeneous cylindrical magnetic flux tube embedded in a homogeneous magnetic environment is investigated. Both the tube and its environment are in steady state. Steady flows break the symmetry of forward (field-aligned) and backward (anti-parallel to magnetic field) propagating MHD wave modes because of the induced Doppler shifts. It is shown that strong enough flows change the sense of propagation of MHD waves. The flow also induces shifts in cut-off values and phase-speeds of the waves. Under photospheric conditions, if the flow is strong enough, the slow surface modes may disappear and the fast body modes may become present. The crossing of modes is also observed due to the presence of flows. The effect of steady-state background has to be considered particularly carefully when evaluating observation signatures of MHD waves for diagnostics in the solar atmosphere.     

The heating of the solar corona

The heating of the solar corona has been a fundamental astrophysical issue for over sixty years. Over the last decade in particular, space-based solar observatories (Yohkoh, SOHO and TRACE) have revealed the complex and often subtle magnetic-field and plasma interactions throughout the solar atmosphere in unprecedented detail. It is now established that any energy release mechanism is magnetic in origin - the challenge posed is to determine what specific heat input is dominating in a given coronal feature throughout the solar cycle. This review outlines a range of possible magnetohydrodynamic (MHD) coronal heating theories, including MHD wave dissipation and MHD reconnection as well as the accumulating observational evidence for quasi-periodic oscillations and small-scale energy bursts occurring in the corona. Also, we describe current attempts to interpret plasma temperature, density and velocity diagnostics in the light of specific localised energy release. The progress in these investigations expected from future solar missions (Solar-B, STEREO, SDO and Solar Orbiter) is also assessed.Received: 6 February 2003, Published online: 14 November 2003 Correspondence to: R. W. Walsh     

Radial variation of differential rotation in the solar electron corona

Long-lived brightness structures in the solar electron corona persist over many solar rotation periods and permit an observational determination of coronal magnetic tracer rotation as a function of latitude and height in the solar atmosphere. For observations over 1964–1976 spanning solar cycle 20, we compare the latitude dependence of rotation at two heights in the corona. Comparison of rotation rates from East and West limbs and from independent computational procedures is used to estimate uncertainty. Time-averaged rotation rates based on three methods of analysis demonstrate that, on average, coronal differential rotation decreases with height from 1.125 to 1.5 R _S. The observed radial variation of differential rotation implies a scale height of approximately 0.7 R _S for coronal differential rotation.Model calculations for a simple MHD loop show that magnetic connections between high and low latitudes may produce the observed radial variations of magnetic tracer rotation. If the observed tracer rotation represents the rotation of open magnetic field lines as well as that of closed loops, the small scale height for differential rotation suggests that the rotation of solar magnetic fields at the base of the solar wind may be only weakly latitude dependent. If, instead, closed loops account completely for the radial gradients of rotation, outward extrapolation of electron coronal rotation may not describe magnetic field rotation at the solar wind source. Inward extrapolations of observed rotation rates suggest that magnetic field and plasma are coupled a few hundredths of a solar radius beneath the photosphere.     

Kolmogorov dissipation scales in weakly ionized plasmas

In a weakly ionized plasma, the evolution of the magnetic field is described by a 'generalized Ohm's law' that includes the Hall effect and the ambipolar diffusion terms. These terms introduce additional spatial and time-scales which play a decisive role in the cascading and the dissipation mechanisms in magnetohydrodynamic turbulence. We determine the Kolmogorov dissipation scales for the viscous, the resistive and the ambipolar dissipation mechanisms. The plasma, depending on its properties and the energy injection rate, may preferentially select one of these dissipation scales, thus determining the shortest spatial scale of the supposedly self-similar spectral distribution of the magnetic field. The results are illustrated taking the partially ionized part of the solar atmosphere as an example. Thus, the shortest spatial scale of the supposedly self-similar spectral distribution of the solar magnetic field is determined by any of the four dissipation scales given by the viscosity, the Spitzer resistivity (electron–ion collisions), the resistivity due to electron–neutral collisions and the ambipolar diffusivity. It is found that the ambipolar diffusion dominates for reasonably large energy injection rate. The robustness of the magnetic helicity in the partially ionized solar atmosphere would facilitate the formation of self-organized vortical structures.     

Current sheet creation by a super-Alfvénic jet in a bipolar field

Energy accumulation in a current sheet (CS) can occur during the injection of a fast plasma jet in a perpendicular magnetic field. A similar situation can occur in the solar corona when a flux of plasma appears under a magnetic arch. The flare can be produced at the CS disruption. The CS creation during plasma jet interaction with the magnetic field is demonstrated by numerical MHD simulation. The choice of dimensionless parameters Re, Re_m,, II, which are suitable for simulation of coronal phenomena, is discussed. When jet injection ceases, the CS evolution produces an unstable state and fast magnetic energy dissipation is observed.     

Magnetohydrodynamic Waves in an Asymmetric Magnetic Slab

Analytical models of solar atmospheric magnetic structures have been crucial for our understanding of magnetohydrodynamic (MHD) wave behaviour and in the development of the field of solar magneto-seismology. Here, an analytical approach is used to derive the dispersion relation for MHD waves in a magnetic slab of homogeneous plasma enclosed on its two sides by non-magnetic, semi-infinite plasma with different densities and temperatures. This generalises the classic magnetic slab model, which is symmetric about the slab. The dispersion relation, unlike that governing a symmetric slab, cannot be decoupled into the well-known sausage and kink modes, i.e. the modes have mixed properties. The eigenmodes of an asymmetric magnetic slab are better labelled as quasi-sausage and quasi-kink modes. Given that the solar atmosphere is highly inhomogeneous, this has implications for MHD mode identification in a range of solar structures. A parametric analysis of how the mode properties (in particular the phase speed, eigenfrequencies, and amplitudes) vary in terms of the introduced asymmetry is conducted. In particular, avoided crossings occur between quasi-sausage and quasi-kink surface modes, allowing modes to adopt different properties for different parameters in the external region.     

Linear MHD Wave Propagation in Time-Dependent Flux Tube

The propagation of magnetohydrodynamic (MHD) waves is an area that has been thoroughly studied for idealised static and steady state magnetised plasma systems applied to numerous solar structures. By applying the generalisation of a temporally varying background density to an open magnetic flux tube, mimicking the observed slow evolution of such waveguides in the solar atmosphere, further investigations into the propagation of both fast and slow MHD waves can take place. The assumption of a zero-beta plasma (no gas pressure) was applied in Williamson and Erdélyi (Solar Phys. 2013, doi: 10.1007/s11207-013-0366-9 , Paper I) is now relaxed for further analysis here. Firstly, the introduction of a finite thermal pressure to the magnetic flux tube equilibrium modifies the existence of fast MHD waves which are directly comparable to their counterparts found in Paper I. Further, as a direct consequence of the non-zero kinetic plasma pressure, a slow MHD wave now exists, and is investigated. Analysis of the slow wave shows that, similar to the fast MHD wave, wave amplitude amplification takes place in time and height. The evolution of the wave amplitude is determined here analytically. We conclude that for a temporally slowly decreasing background density both propagating magnetosonic wave modes are amplified for over-dense magnetic flux tubes. This information can be very practical and useful for future solar magneto-seismology applications in the study of the amplitude and frequency properties of MHD waveguides, e.g. for diagnostic purposes, present in the solar atmosphere.     

太阳大气中的光球亮点

首先是对太阳光球亮点近年来研究工作的总结。光球亮点是一种发生在太阳光球上宁静区域的的小尺度和短寿命增亮现象,平均直径在100～300knm之间,平均寿命约为几分钟。光球亮点的研究对于光球辐射和磁场性质的认识具有重要意义。过去的观测显示,绝大多数光球亮点的产生和演化与磁场,特别是光球上的小尺度磁场的演化密切相关,比如,两个同极性磁场的合并,或者反极性磁场的对消,或者一个同极性磁场的分裂,均可以促使光球亮点产生或消失。基于这样的观测结果,统计研究了2722个光球亮点(1600A)与光球上偶极磁元的关系,发现大约有1/3的光球亮点出现在偶极磁元中心附近。     

Generation of rotational discontinuities by magnetic reconnection associated with microflares

Magnetic reconnection can take place between two plasma regions with antiparallel magnetic field components. In a time-dependent reconnection event, the plasma outflow region consists of a leading bulge region and a trailing reconnection layer. Magnetohydrodynamic (MHD) discontinuities, including rotational discontinuities, can be formed in both the bulge region and the trailing layer. In this paper, we suggest that the rotational discontinuities observed in the solar wind may be generated by magnetic reconnection associated with microflares in coronal holes. The structure of the reconnection layer is studied by solving the one-dimensional Riemann problem for the evolution of an initial current sheet after the onset of magnetic reconnection as well as carrying out two-dimensional MHD simulations. As the emerging magnetic flux reconnects with ambient open magnetic fields in the coronal hole, rotational discontinuities are generated in the region with open field lines. It is also found that in the solar corona with a low plasma beta ( 0.01), the magnetic energy is converted through magnetic reconnection mostly into the plasma bulk-flow energy. Since more microflares will generate more rotational discontinuities and also supply more energy to the solar wind, it is expected that the number of rotational discontinuities observed in the solar wind would be an increasing function of solar wind speed. The observation rate of rotational discontinuities generated by microflares is estimated to be dN _RD/dt - f/63 000 s (f > 1) at 1 AU. The present mechanism favors the generation of rotational discontinuities with a large shock normal angle.     

Excess heating of corona and chromosphere above magnetic regions by non-linear Alfvén waves

Excess heating of the active region solar atmosphere is interpreted by the decay of MHD slow-mode waves produced in the corona through the non-linear coupling of Alfvén waves supplied from subphotospheric layers. It is stressed that the Alfvén-mode waves may be very efficiently generated directly in the convection layer under the photosphere in magnetic regions, and that such magnetic regions, at the same time, provide the ‘transparent windows’ for Alfvén waves in regard to the Joule and frictional dissipations in the photospheric and subphotospheric layers. Though the Alfvén waves suffer considerable reflection in the chromosphere and in the transition layer, a certain fraction of this large flux is propagated out to the corona, and a large velocity amplitude exceeding the local Alfvén velocity is attained during the propagation along the magnetic tubes of force into a region of lower density and weaker magnetic field. The otherwise divergence-free velocity field in Alfvén waves gets involved in such a case with a compressional component (slow-mode waves) which again is of considerable velocity amplitude relative to the local acoustic velocity when estimated by using the formulation for non-linear coupling between MHD wave modes derived by Kaburaki and Uchida (1971). Therefore, the compressional waves thus produced through the non-linear coupling of Alvén waves will eventually be thermalized to provide a heat source. The introduction of this non-linear coupling process and the subsequent thermalization of thus produced slow-mode waves may provide means of converting the otherwise dissipation-free Alfvén mode energy into heat in the corona. The liberated heat will readily be redistributed by conduction along the magnetic lines of force, with higher density as a consequence of increased scale height, and thus the loop-like structure of the coronal condensations (or probably also the thread-like feature of the general corona) may be explained in a natural fashion.     

Interfacing MHD Single Fluid and Kinetic Exospheric Solar Wind Models and Comparing Their Energetics

An exospheric kinetic solar wind model is interfaced with an observation-driven single-fluid magnetohydrodynamic (MHD) model. Initially, a photospheric magnetogram serves as observational input in the fluid approach to extrapolate the heliospheric magnetic field. Then semi-empirical coronal models are used for estimating the plasma characteristics up to a heliocentric distance of 0.1 AU. From there on, a full MHD model that computes the three-dimensional time-dependent evolution of the solar wind macroscopic variables up to the orbit of Earth is used. After interfacing the density and velocity at the inner MHD boundary, we compare our results with those of a kinetic exospheric solar wind model based on the assumption of Maxwell and Kappa velocity distribution functions for protons and electrons, respectively, as well as with in situ observations at 1 AU. This provides insight into more physically detailed processes, such as coronal heating and solar wind acceleration, which naturally arise from including suprathermal electrons in the model. We are interested in the profile of the solar wind speed and density at 1 AU, in characterizing the slow and fast source regions of the wind, and in comparing MHD with exospheric models in similar conditions. We calculate the energetics of both models from low to high heliocentric distances.     

A Simulation Study of the Magnetic Reconnections Between Closed Loops and Open Magnetic Fields Driven by Horizontal Flows on Solar Surface

Numerous observational events in the solar atmosphere (e.g., solar ?ares and jets) are attributed to the energy conversion due to magnetic reconnec- tions. Magnetic reconnections are also involved in a new scenario of solar wind origin to play a crucial role in opening the closed magnetic loop and releasing its mass into the open magnetic funnel. In this scenario, the closed magnetic loop moves towards the supergranular boundary by the supergranular convection, and collides with the open magnetic funnel there to trigger the magnetic reconnec- tion between each other. This work aims at studying the occurrence and effect of magnetic reconnection in this scenario in detail. The magnetohydrodynamic (MHD) numerical simulation is an important approach to investigate the mag- netic reconnection process in the solar atmosphere. A two-dimensional MHD numerical model has been established, and in combination with the strati?ed temperature and density distributions in the solar atmosphere, the numerical simulation on the process of magnetic reconnection of the closed magnetic loops driven by the horizontal ?ows with the open magnetic ?elds has been performed on the scale of supergranulation. Based on a quantitative analysis of the simula- tion result, it is suggested that the process of magnetic reconnection can really realize the mass release of closed magnetic loops, and further supply to the new open magnetic structures to produce upward mass ?ows. Our results provide a basis for the further modeling of solar wind origin.     

A comparison of global models for the solar wind interaction with Mars

We present initial results from the first community-wide effort to compare global plasma interaction model results for Mars. Seven modeling groups participated in this activity, using MHD, multi-fluid, and hybrid assumptions in their simulations. Moderate solar wind and solar EUV conditions were chosen, and the conditions were implemented in the models and run to steady state. Model output was compared in three ways to determine how pressure was partitioned and conserved in each model, the location and asymmetry of plasma boundaries and pathways for planetary ion escape, and the total escape flux of planetary oxygen ions. The two participating MHD models provided similar results, while the five sets of multi-fluid and hybrid results were different in many ways. All hybrid results, however, showed two main channels for oxygen ion escape (a pickup ion ‘plume’ in the hemisphere toward which the solar wind convection electric field is directed, and a channel in the opposite hemisphere of the central magnetotail), while the MHD models showed one (a roughly symmetric channel in the central magnetotail). Most models showed a transition from an upstream region dominated by plasma dynamic pressure to a magnetosheath region dominated by thermal pressure to a low altitude region dominated by magnetic pressure. However, calculated escape rates for a single ion species varied by roughly an order of magnitude for similar input conditions, suggesting that the uncertainties in both the current and integrated escape over martian history as determined by models are large. These uncertainties are in addition to those associated with the evolution of the Sun, the martian dynamo, and the early atmosphere, highlighting the challenges we face in constructing Mars’ past using models.