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)     

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.     

Changes in convective properties over the solar cycle: effect on p-mode damping rates

Measurements of both solar irradiance and p-mode oscillation frequencies indicate that the structure of the Sun changes with the solar cycle. Balmforth, Gough & Merryfield investigated the effect of symmetrical thermal disturbances on the solar structure and the resulting pulsation frequency changes. They concluded that thermal perturbations alone cannot account for the variations in both irradiance and p-mode frequencies, and that the presence of a magnetic field affecting acoustical propagation is the most likely explanation of the frequency change, in the manner suggested earlier by Gough & Thompson and by Goldreich et al. Numerical simulations of Boussinesq convection in a magnetic field have shown that at high Rayleigh number the magnetic field can modify the preferred horizontal length scale of the convective flow.
Here, we investigate the effect of changing the horizontal length scale of convective eddies on the linewidths of the acoustic resonant mode peaks observed in helioseismic power spectra. The turbulent fluxes in these model computations are obtained from a time-dependent, non-local generalization of the mixing-length formalism. The modelled variations are compared with p-mode linewidth changes revealed by the analysis of helioseismic data collected by the Birmingham Solar-Oscillations Network (BiSON); these low-degree (low- l ) observations cover the complete falling phase of solar activity cycle 22. The results are also discussed in the light of observations of solar-cycle variations of the horizontal size of granules and with results from 2D simulations by Steffen of convective granules.     

On helioseismic tests of basic physics

An important goal of helioseismology is to provide information about the basic physics and parameters that determine the structure of the solar interior. Here we discuss the procedures applied in such analyses, using as an example attempts to obtain significant constraints on the value of Newton's gravitational constant G from helioseismology. The analysis is based on complete direct and inverse helioseismic analysis of a set of accurate observed acoustic frequencies. We confirm, as found by previous investigations based on different approaches, that the actual level of precision of the helioseismic inferences does not allow us to constrain G with a precision better than that which can be reached with direct experimental measurements. The conclusion emphasizes the importance in helioseismic inferences of considering not only the accuracy with which solar oscillations are measured, but also the effect of uncertainties in other aspects of the model computation and helioseismic analysis.     

On the variations in the parameters of high-degree acoustic modes with solar cycle

We analyze the pattern of behavior of p-mode wave packets with solar cycle using TON one-day helioseismic data with a high spatial resolution. The time—distance method is used to perform this task. We make an attempt to determine the variations in the travel time of acoustic waves at maximum and minimum solar activity; at maximum activity, this time decreases by 2 s compared to that at minimum activity to a depth of 0.8R_⊙. In addition, the correlation amplitudes of acoustic wave packets from minimum to maximum solar activity were found to decrease by 10–20% for all angular distances.     

Turbulent spectrum of Alfvén waves excited by a kinetic instability for explaining the modulations with multi-timescales in solar flares

A low-frequency wave is treated as a local oscillation to modulate the guiding center of electrons beam, which is considered as free energy to excite Alfvén waves by a kinetic plasma instability under low-frequency approximation. The nonlinearity of the model is shown by a critical value of the amplitude of the low-frequency wave, and Alfvén waves are growing in a broad turbulent spectrum with fractional harmonics, which strongly depend on the criterion. The instability is limited in the direction nearly perpendicular to the ambient magnetic field. The growth rates are very sensitive to the beam speed that perpendicular to the magnetic field, the propagational angle, and the magnetic field strength, but not sensitive to the beam speed parallel to the magnetic field. This model is used to explain the modulations with multiple timescales in the flare light curves at radio, hard X-ray and H-alpha bands.     

Magnetic and thermal phase shifts in the local helioseismology of sunspots

Phase perturbations due to inclined surface magnetic field of active region strength are calculated numerically in quiet Sun and simple sunspot models in order to estimate and compare the direct and indirect (thermal) effects of the fields on helioseismic waves. It is found that the largest direct effects occur in highly inclined field characteristic of penumbrae, and scale roughly linearly with magnetic field strength. The combined effects of sunspot magnetic and thermal anomalies typically yield negative travel-time perturbations in penumbrae. Travel-time shifts in umbrae depend on details of how the thermal and density structure differs from the quiet Sun. The combined shifts are generally not well approximated by the sum of the thermal and magnetic effects applied separately, except at low field strengths of around 1 kG or less, or if the thermal shift is small. A useful rule-of-thumb appears to be that travel-time perturbations in umbrae are predominantly thermal, whereas in penumbrae they are mostly magnetic.     

Meridional flow profile measurements with SOHO/MDI

We present meridional flow measurements of the Sun using a novel helioseismic approach for analyzing SOHO/MDI data in order to push the current limits in radial depth. Analyzing three consecutive months of data during solar minimum, we find that the meridional flow is as expected poleward in the upper convection zone, turns equatorward at a depth of around 40 Mm (∼ 0.95 R_⊙), and possibly changes direction again in the lower convection zone. This may indicate two meridional circulation cells in each hemisphere, one beneath the other. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Probing sunspot magnetic fields with p-mode absorption and phase shift data

Long-standing observations of incoming and outgoing f- and p-modes in annuli around sunspots reveal that the spots partially absorb and substantially shift the phase of waves incident upon them. The commonly favoured absorption mechanism is partial conversion to slow magneto-acoustic waves that disappear into the solar interior channelled by the magnetic field of the sunspot. However, up until now, only f-mode absorption could be accounted for quantitatively by this means. Based on vertical magnetic field models, the absorption of p-modes was insufficient. In this paper, we use the new calculations of Crouch & Cally for inclined fields, and a simplified model of the interaction between spot interior and exterior. We find excellent agreement with phase shift data assuming field angles from the vertical in excess of 30° and Alfvén/acoustic equipartition depths of around 600–800 km. The absorption of f-modes produced by such models is considerably larger than is observed, but consistent with numerical simulations. On the other hand, p-mode absorption is generally consistent with observed values, up to some moderate frequency dependent on radial order. Thereafter, it is too large, assuming absorbing regions comparable in size to the inferred phase-shifting region. The excess absorption produced by the models is in stark contrast with previous calculations based on a vertical magnetic field, and is probably due to finite mode lifetimes and excess emission in acoustic glories. The excellent agreement of phase shift predictions with observational data allows some degree of probing of subsurface field strengths, and opens up the possibility of more accurate inversions using improved models. Most importantly, though, we have confirmed that slow mode conversion is a viable, and indeed the likely, cause of the observed absorption and phase shifts.     

Properties of filaments in solar cycle 20-23 from McIntosh Archive

A filament is a cool, dense structure suspended in the solar corona. The eruption of a filament is often associated with a coronal mass ejection(CME), which has an adverse effect on space weather. Hence,research on filaments has attracted much attention in the recent past. The tilt angle of active region(AR)magnetic bipoles is a crucial parameter in the context of the solar dynamo, which governs the conversion efficiency of the toroidal magnetic field to poloidal magnetic field. Filaments always form over polarity inversion lines(PILs), so the study of tilt angles for these filaments can provide valuable information about generation of a magnetic field in the Sun. We investigate the tilt angles of filaments and other properties using McIntosh Archive data. We fit a straight line to each filament to estimate its tilt angle. We examine the variation of mean tilt angle with time. The latitude distribution of positive tilt angle filaments and negative tilt angle filaments reveals that there is a dominance of positive tilt angle filaments in the southern hemisphere and negative tilt angle filaments dominate in the northern hemisphere. We study the variation of the mean tilt angle for low and high latitudes separately. Investigations of temporal variation with filament number indicate that total filament number and low latitude filament number vary cyclically, in phase with the solar cycle. There are fewer filaments at high latitudes and they also show a cyclic pattern in temporal variation. We also study the north-south asymmetry of filaments with different latitude criteria.     

Detection and temporal coherence of p‐modes below 1.4 mHz

Data collected recently by the helioseismic experiments aboard the SOHO spacecraft have allowed the detection of low degree p‐modes with increasingly lower order n. In particular, the GOLF experiment is currently able to unambiguously identify low degree modes with frequencies as low as 1.3 mHz. The detection of p‐modes with very low frequency (i.e., low n), is difficult due to the low signal‐to‐noise ratio in this spectral region and its contamination by solar signals that are not of acoustic origin. To address this problem without using any theoretical a priory, we propose a methodology that relies only on the inversion of observed values to define a spectral window for the expected locations of these low frequency modes. The application of this method to 2920‐day‐long GOLF observations is presented and its results discussed. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Constraints on angular momentum transport in the Sun from simulations of the tachocline

I present a new numerical tool for studying the interaction of meridional flows and magnetic fields, and study their role in establishing angular-momentum balance in the solar radiative zone. Quantitative comparisons with helioseismic observations provide stringent constraints on existing models of the dynamics of the solar interior. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A study of possible temporal and latitudinal variations in the properties of the solar tachocline

Temporal variations of the structure and the rotation rate of the solar tachocline region are studied using helioseismic data from the Global Oscillation Network Group (GONG) and the Michelson Doppler Imager (MDI) obtained during the period 1995–2000. We do not find any significant temporal variation in the depth of the convection zone, the position of the tachocline or the extent of overshoot below the convection zone. No systematic variation in any other properties of the tachocline, like width, etc., is found either. The possibility of periodic variations in these properties is also investigated. Time-averaged results show that the tachocline is prolate with a variation of about 0.02 R_⊙ in its position. Neither the depth of the convection zone nor the extent of overshoot shows any significant variation with latitude.     

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°.     

Development of structure in pores and sunspots: flows around axisymmetric magnetic flux tubes

Flux elements, pores and sunspots form a family of magnetic features observed at the solar surface. As a first step towards developing a fully non-linear model of the structure of these features and of the dynamics of their interaction with solar convection, we conduct numerical experiments on idealized axisymmetric flux tubes in a compressible convecting atmosphere in cylindrical boxes of radius up to 8 times their depth. We find that the magnetic field strength of the flux tubes is roughly independent of both distance from the centre and the total flux content of the flux tube, but that the angle of inclination from the vertical of the field at the edge of the tube increases with flux content. In all our calculations, fluid motion converges on the flux tube at the surface. The results compare favourably with observations of pores; in contrast, large sunspots lie at the centre of an out-flowing moat cell. We conjecture that there is an inflow hidden beneath the penumbrae of large spots, and that this inflow is responsible for the remarkable longevity of such features.     

关于磁声重力波传播性质的探讨

本文对充满垂直均匀磁场的等温大气内的磁声重力波做了严格的解析分析,并将其通解表述成广义超几何函数的形式。该解可用于对磁大气内振荡现象的进一步数值模拟研究。对解的分析澄清了若干磁声重力波的传播性质。     

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.     

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.     

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.