Modeling helioseismic modes in the magnetic atmosphere of the Sun

The pulsation of the solar surface is caused by acoustic waves traveling in the solar interior. Thorough analyses of observational data indicate that these f and p helioseismic oscillation modes are not bounced back completely at the surface but they partially penetrate into the atmosphere. Atmospheric effects and their possible observational application are investigated in one‐dimensional magnetohydrodynamic models. It is found that f and p mode frequencies are shifted of the order of μHz due to the presence of an atmospheric magnetic field. This shift varies with the direction of the wave propagation.Resonant coupling of global helioseismic modes to local Alfvén and slow waves reduce the life time of the global modes. The resulting line width of the frequency line is of the order of nHz, and it also varies with propagation angle. These features enable us to use helioseismic observations in magnetic diagnostics of the lower atmosphere. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling p-Mode Interaction with a Spreading Sunspot Field

Sunspots absorb and scatter incident p modes. The dominant mechanism is still uncertain. One possibility, mode conversion to slow magneto-acoustic waves, has been shown to yield results in agreement with observations for the f mode only. Absorption of p modes in simple vertical magnetic field models is too weak by an order of magnitude or more. Here we report on numerical calculations of p modes encountering a simple sunspot model with field which spreads with height. It is found that p-mode absorption is greatly enhanced by field spread, to a level consistent with observations, and it appears that it occurs preferentially in the outer regions of the spot, in line with recent results from acoustic holography.     

Imaging bright-spots in the accretion flow near the black hole horizon of Sgr A*

Solar active regions are distinguished by their strong magnetic fields. Modern local helioseismology seeks to probe them by observing waves which emerge at the solar surface having passed through their interiors. We address the question of how an acoustic wave from below is partially converted to magnetic waves as it passes through a vertical magnetic field layer where the sound and Alfvén speeds coincide (the equipartition level), and find that (i) there is no associated reflection at this depth, either acoustic or magnetic, only transmission and conversion to an ongoing magnetic wave; and (ii) conversion in active regions is likely to be strong, though not total, at frequencies typically used in local helioseismology, with lower frequencies less strongly converted. A simple analytical formula is presented for the acoustic-to-magnetic conversion coefficient.     

Mode Conversion of Solar <Emphasis Type="Italic">p</Emphasis>-Modes in Non-Vertical Magnetic Fields

Sunspots absorb and scatter incident f- and p-modes. Until recently, the responsible absorption mechanism was uncertain. The most promising explanation appears to be mode conversion to slow magnetoacoustic-gravity waves, which carry energy down the magnetic field lines into the interior. In vertical magnetic field, mode conversion can adequately explain the observed f-mode absorption, but is too inefficient to account for the absorption of p-modes. In the first paper of the present series we calculated the efficiency of fast-to-slow magnetoacoustic-gravity wave conversion in uniform non-vertical magnetic fields. We assumed two-dimensional propagation, where the Alfvén waves decouple. In comparison to vertical field, it was found that mode conversion is significantly enhanced in moderately inclined fields, especially at higher frequencies. Using those results, Cally, Crouch, and Braun showed that the resultant p-mode absorption produced by simple sunspot models with non-vertical magnetic fields is ample to explain the observations. In this paper, we further examine mode conversion in non-vertical magnetic fields. In particular, we consider three-dimensional propagation, where the fast and slow magnetoacoustic-gravity waves and the Alfvén waves are coupled. Broadly speaking, the p-mode damping rates are not substantially different to the two-dimensional case. However, we do find that the Alfvén waves can remove similar quantities of energy to the slow MAG waves.     

On mode conversion and wave reflection in magnetic Ap stars

We investigate the effect of a strong large-scale magnetic field on the reflection of high-frequency acoustic modes in rapidly oscillating Ap stars. To that end, we consider a toy model composed of an isothermal atmosphere matched on to a polytropic interior and determine the numerical solution to the set of ideal magnetohydrodynamic equations in a local plane-parallel approximation with constant gravity. Using the numerical solution in combination with approximate analytical solutions that are valid in the limits where the magnetic and acoustic components are decoupled, we calculate the relative fraction of energy flux that is carried away in each oscillation cycle by running acoustic waves in the atmosphere and running magnetic waves in the interior. For oscillation frequencies above the acoustic cut-off, we show that most energy losses associated with the presence of running waves occur in regions where the magnetic field is close to vertical. Moreover, by considering the depth dependence of the energy associated with the magnetic component of the wave in the atmosphere we show that a fraction of the wave energy is kept in the oscillation every cycle. For frequencies above the acoustic cut-off frequency, such energy is concentrated in regions where the magnetic field is significantly inclined in relation to the local vertical. Even though our calculations were aimed at studying oscillations with frequencies above the acoustic cut-off frequency, based on our results we discuss what results may be expected for oscillations of lower frequency.     

Distribution of magnetic flux on the solar surface and low-degree p-modes

The frequencies of solar p-modes are known to change over the solar cycle. There is also recent evidence that the relation between frequency shift of low-degree modes and magnetic flux or other activity indicators differs between the rising and falling phases of the solar cycle, leading to a hysteresis in such diagrams. We consider the influence of the changing large-scale surface distribution of the magnetic flux on low-degree ( l ≤3) p-mode frequencies. To that end, we use time-dependent models of the magnetic flux distribution and study the ensuing frequency shifts of modes with different order and degree as a function of time. The resulting curves are periodic functions (in simple cases just sine curves) shifted in time by different amounts for the different modes. We show how this may easily lead to hysteresis cycles comparable to those observed. Our models suggest that high-latitude fields are necessary to produce a significant difference in hysteresis between odd- and even-degree modes. Only magnetic field distributions within a small parameter range are consistent with the observations by Jiménez-Reyes et al. Observations of p-mode frequency shifts are therefore capable of providing an additional diagnostic of the magnetic field near the solar poles. The magnetic distribution that is consistent with the p-mode observations also appears reasonable compared with direct measurements of the magnetic field.     

Physical Properties of Wave Motion in Inclined Magnetic Fields within Sunspot Penumbrae

At the surface of the Sun, acoustic waves appear to be affected by the presence of strong magnetic fields in active regions. We explore the possibility that the inclined magnetic field in sunspot penumbrae may convert primarily vertically-propagating acoustic waves into elliptical motion. We use helioseismic holography to measure the modulus and phase of the correlation between incoming acoustic waves and the local surface motion within two sunspots. These correlations are modeled by assuming the surface motion to be elliptical, and we explore the properties of the elliptical motion on the magnetic-field inclination. We also demonstrate that the phase shift of the outward-propagating waves is opposite to the phase shift of the inward-propagating waves in stronger, more vertical fields, but similar to the inward phase shifts in weaker, more-inclined fields.     

Three-Dimensional MHD Wave Propagation and Conversion to Alfvén Waves near the Solar Surface. I. Direct Numerical Solution

The efficacy of fast?–?slow MHD mode conversion in the surface layers of sunspots has been demonstrated over recent years using a number of modelling techniques, including ray theory, perturbation theory, differential eigensystem analysis, and direct numerical simulation. These show that significant energy may be transferred between the fast and slow modes in the neighbourhood of the equipartition layer where the Alfvén and sound speeds coincide. However, most of the models so far have been two dimensional. In three dimensions the Alfvén wave may couple to the magnetoacoustic waves with important implications for energy loss from helioseismic modes and for oscillations in the atmosphere above the spot. In this paper, we carry out a numerical “scattering experiment,” placing an acoustic driver 4 Mm below the solar surface and monitoring the acoustic and Alfvénic wave energy flux high in an isothermal atmosphere placed above it. These calculations indeed show that energy conversion to upward travelling Alfvén waves can be substantial, in many cases exceeding loss to slow (acoustic) waves. Typically, at penumbral magnetic field strengths, the strongest Alfvén fluxes are produced when the field is inclined 30°?–?40° from the vertical, with the vertical plane of wave propagation offset from the vertical plane containing field lines by some 60°?–?80°.     

The Effect of Abnormal Granulation on Acoustic Wave Travel Times and Mode Frequencies

Observations indicate that in plage areas (i.e. in active regions outside sunspots) acoustic waves travel faster than in the quiet Sun, leading to shortened travel times and higher p-mode frequencies. Coupled with the 11-year variation of solar activity, this may also explain the solar cycle variation of oscillation frequencies. While it is clear that the ultimate cause of any difference between the quiet Sun and plage is the presence of magnetic fields of order 100 G in the latter, the mechanism by which the magnetic field exerts its influence has not yet been conclusively identified. One possible such mechanism is suggested by the observation that granular motions in plage areas tend to be slightly “abnormal”, dampened compared to the quiet Sun. In this paper we consider the effect that abnormal granulation observed in active regions should have on the propagation of acoustic waves. Any such effect is found to be limited to a shallow surface layer where sound waves propagate nearly vertically. The magnetically suppressed turbulence implies higher sound speeds, leading to shorter travel times. This time shift Δ τ is independent of the travel distance, while it shows a characteristic dependence on the assumed plage field strength. As a consequence of the variation of the acoustic cutoff with height, Δ τ is expected to be significantly higher for higher frequency waves within the observed regime of 3 – 5 mHz. The lower group velocity near the upper reflection point further leads to an increased envelope time shift, as compared to the phase shift. p-mode frequencies in plage areas are increased by a corresponding amount, Δ ν/ν=ν Δ τ. These characteristics of the time and frequency shifts are in accordance with observations. The calculated overall amplitudes of the time and frequency shifts are comparable to, but still significantly less than (by a factor of 2 to 5), those suggested by measurements.     

Oscillatory Modes of a Prominence – PCTR – Corona Slab Model

Oscillations of magnetic structures in the solar corona have often been interpreted in terms of magnetohydrodynamic waves. We study the adiabatic magnetoacoustic modes of a prominence plasma slab with a uniform longitudinal magnetic field, surrounded by a prominence – corona transition region (PCTR) and a coronal medium. Considering linear small-amplitude oscillations, we deduce the dispersion relation for the magnetoacoustic slow and fast modes by assuming evanescentlike perturbations in the coronal medium. In the system without PCTR, a classification of the oscillatory modes according to the polarisation of their eigenfunctions is made to distinguish modes with fastlike or slowlike properties. Internal and external slow modes are governed by the prominence and coronal properties, respectively, and fast modes are mostly dominated by prominence conditions for the observed wavelengths. In addition, the inclusion of an isothermal PCTR does not substantially influence the mode frequencies, but new solutions (PCTR slow modes) are present.     

Diagnostics of the Internal Structure of Stars using the Differential Response Technique

We address the problem of the diagnosing the deep interior structure of stars using acoustic p-modes, and investigate the diagnostic capabilities of two complementary approaches both based on the differential response technique (Vorontsov, 1998): (a) direct calibration using a grid of evolutionary stellar models, and (b) linear and non-linear (with consecutive linearisations) inversion of low-degree frequencies. We apply this analysis to the frequencies of a model of an old 0.8M star, and to the solar frequencies obtained from BiSON measurements, using a 2-D grid of reference models of different mass and age. We explore the convergence and stability of the asteroseismic inversion, performed with the adaptive regularisation technique of Strakhov and Vorontsov (2001).     

Mode Conversion of Solar p Modes in non-Vertical Magnetic Fields – i. two-Dimensional Model

Sunspots absorb incident p modes. The responsible mechanism is uncertain. One possibility is mode conversion to slow magnetoacoustic–gravity waves. In vertical field mode conversion can adequately explain the observed f-mode absorption, but is too inefficient to explain the absorption of p modes. In this investigation we calculate the efficiency of fast-to-slow magnetoacoustic–gravity wave conversion in non-vertical field. We assume two-dimensional propagation where the Alfvén waves decouple. It is found that resultant p-mode absorption is significantly enhanced for moderate inclinations at higher frequencies, whereas for p modes at lower frequencies, and the f mode in general, there is no useful enhancement. However, the enhancement is insufficient to explain the observed p-mode absorption by sunspots. Paper II considers the efficiency of mode conversion in non-vertical field with three-dimensional propagation, where fast and slow magnetoacoustic–gravity waves and Alfvén waves are coupled.     

A New Model for Propagating Parts of EIT Waves: A Current Shell in a CME

EIT waves are observed in EUV as bright fronts. Some of these bright fronts propagate across the solar disk. EIT waves are all associated with a flare and a CME and are commonly interpreted as fast-mode magnetosonic waves. Propagating EIT waves could also be the direct signature of the gradual opening of magnetic field lines during a CME. We quantitatively addressed this alternative interpretation. Using two independent 3D MHD codes, we performed nondimensional numerical simulations of a slowly rotating magnetic bipole, which progressively result in the formation of a twisted magnetic flux tube and its fast expansion, as during a CME. We analyse the origins, the development, and the observability in EUV of the narrow electric currents sheets that appear in the simulations. Both codes give similar results, which we confront with two well-known SOHO/EIT observations of propagating EIT waves (7 April and 12 May 1997), by scaling the vertical magnetic field components of the simulated bipole to the line of sight magnetic field observed by SOHO/MDI and the sign of helicity to the orientation of the soft X-ray sigmoids observed by Yohkoh/SXT. A large-scale and narrow current shell appears around the twisted flux tube in the dynamic phase of its expansion. This current shell is formed by the return currents of the system, which separate the twisted flux tube from the surrounding fields. It intensifies as the flux tube accelerates and it is co-spatial with weak plasma compression. The current density integrated over the altitude has the shape of an ellipse, which expands and rotates when viewed from above, reproducing the generic properties of propagating EIT waves. The timing, orientation, and location of bright and faint patches observed in the two EIT waves are remarkably well reproduced. We conjecture that propagating EIT waves are the observational signature of Joule heating in electric current shells, which separate expanding flux tubes from their surrounding fields during CMEs or plasma compression inside this current shell. We also conjecture that the bright edges of halo CMEs show the plasma compression in these current shells.     

A New Look at Mode Conversion in a Stratified Isothermal Atmosphere

Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641, 2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-β plasma to a high-β plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-β to high-β environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.     

On transition region convection cells in simulations of p‐mode propagation into the solar atmosphere

This paper examines the way that transition region surface waves, generated in 2‐D numerical simulations of the nonmagnetic solar atmosphere when various synthetic photospheric drivers are applied, drive the granulation of the transition region/lower coronal region into convection cells. It is shown that these cells are generated by both synthetic point drivers and synthetic horizontally coherent p‐mode drivers. These cells cause the conversion of driven signals in vertical velocity into coronal signals predominantly in horizontal velocity, which if carried over to a case with a magnetic field included could cause mode conversion. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic field inclination and atmospheric oscillations above solar active regions

Recent observational evidence for magnetic field direction effects on helioseismic signals in sunspot penumbrae is suggestive of magnetohydrodynamic (MHD) mode conversion occurring at lower levels. This possibility is explored using wave mechanical and ray theory in a model of the Sun's surface layers permeated by uniform inclined magnetic field. It is found that fast-to-slow conversion near the equipartition depth at which the sound and Alfvén speeds coincide can indeed greatly enhance the atmospheric acoustic signal at heights observed by Solar and Heliospheric Observatory/Michelson Doppler Imager and other helioseismic instruments, but that this effect depends crucially on the wave attack angle , i.e. the angle between the wavevector and the magnetic field at the conversion/transmission depth. A major consequence of this insight is that the magnetic field acts as a filter, preferentially allowing through acoustic signal from a narrow range of incident directions. This is potentially testable by observation.     

Dissipation of Longitudinal Oscillations in Stratified Nonisothermal Hot Coronal Loops

We investigate the damping of longitudinal (i.e., slow or acoustic) waves in nonisothermal, hot (T≥ 5.0 MK), gravitationally stratified coronal loops. Motivated by SOHO/SUMER and Yohkoh/SXT observations, and by taking into account a range of dissipative mechanisms such as thermal conduction, compressive viscosity, radiative cooling, and heating, the nonlinear governing equations of one-dimensional hydrodynamics are solved numerically for standing-wave oscillations along a magnetic field line. A semicircular shape is chosen to represent the geometry of the coronal loop. It was found that the decay time of standing waves decreases with the increase of the initial temperature, and the periods of oscillations are affected by the different initial footpoint temperatures and loop lengths studied by the numerical experiments. In general, the period of oscillation of standing waves increases and the damping time decreases when the parameter that characterises the temperature at the apex of the loop increases for a fixed footpoint temperature and loop length. A relatively simple second-order scaling polynomial between the damping time and the parameter determining the apex temperature is found. This scaling relation is proposed to be tested observationally. Because of the lack of a larger, statistically relevant number of observational studies of the damping of longitudinal (slow) standing oscillations, it can only be concluded that the numerically predicted decay times are well within the range of values inferred from Doppler shifts observed by SUMER in hot coronal loops.     

The influence of turbulence on the solar p-mode frequencies I. Free-mode approximation

The influence of turbulence on the frequencies of free acoustic modes in convection zones is considered. The frequencies are modified via the speed of sound by the turbulence-induced alterations of the effective pressure: (i) by the correlated fluctuations of temperature and density and (ii) the pressure part of the Reynolds stress. The two effects shift the frequency of low l p-modes in opposite directions. In addition, the correlation of the density fluctuations with the random velocity — the eddy-mass flow — is also relevant. It is, in a steady state, balanced by a vertical mean velocity. The balance results in a rather small net effect completely disappearing for highly nonradial oscillations. Both effects of the density fluctuations produce a redshift of the low l p-mode frequencies. The Reynolds stress, however, makes a blueshift of the frequencies relative to that computed for a laminar gas. This effect dominates for subsonic turbulences. The applied second-order correlation-approximation, however, only holds for the lowest frequencies, where the KORONAS (solar minimum) data are indicating a blueshift. Of particular importance for the present concept is the expected cycle-variations of the lineshifts, i.e. the consideration of the magnetic modification of the various contributions. Observations may show whether the suggested modifications of the solar oscillation theory are correct.     

太阳5分钟振荡的激发机制

1 引言 太阳5分钟振荡是上世纪1个重要的发现[1],它使得人们可以通过观测太阳表面的振动来探测其内部的结构,日震学已取得了巨大的进展,然而我们至今仍不了解其脉动的激发机制,它依然是1个存在争议的问题.太阳位于造父变星脉动不稳定区之外,所以大多数人都相信,由于对流的阻尼,太阳是脉动稳定的,太阳和太阳型恒星的振荡都是由所谓的湍流随机激发机制所激发[2-8].     

Dissipation of Standing Slow Magnetoacoustic Waves in Hot Coronal Loops

The damping of standing slow waves in hot (T>6 MK) coronal loops of semicircular shape is revisited in both the linear and nonlinear regimes. Dissipation by thermal conduction, compressive viscosity, radiative cooling, and heating are examined for nonstratified and stratified loops. We find that for typical conditions of hot SUMER loops, thermal conduction increases the period of damped oscillations over the sound-crossing time, whereas the decay times are mostly shaped by compressive viscosity. Damping from optically thin radiation is negligible. We also find that thermal conduction alone results in slower damping of the density and velocity waves compared to the observations. Only when compressive viscosity is added do these waves damp out at the same rate as the observed rapidly decaying modes of hot SUMER loop oscillations, in contrast to most current work, which has pointed to thermal conduction as the dominant mechanism. We compare the linear predictions with numerical hydrodynamic calculations. Under the effects of gravity, nonlinear viscous dissipation leads to a reduction of the decay time compared to the homogeneous case. In contrast, the linear results predict that the damping rates are barely affected by gravity.