Go to previous page Go up Go to next page

3.5 Effects of small steady perturbations

In order to do linear inversions of helioseismic data, it is necessary to first solve the linear forward problem. The linear forward problem is to compute the first-order effect of small perturbations to the solar model. By first-order we mean first order in the strength of the perturbations. Essentially all inversions that have been done have assumed that the perturbations to the solar model are time-independent over the time duration during which the observations are made.

In general the linear forward problem can be written as

∫ ∑ 3 i δdi = ⊙d r K α(r) δqα(r). (28 ) α
The δdi are the first-order changes in the observed helioseismic parameters, for example changes in travel times, ring-fit parameters, or changes in ingression-egression correlations (which will be described in Section 4). The kernel functions Ki (r) α depend on the three-dimensional position vector r. The first-order changes to the solar model are given by the functions δqα (r ), where the index α refers to the type of the perturbation, for example we can look at the effect of flows, changes in sound speed, and changes in wave excitation and damping. For each type of perturbation α, corresponds a particular kernel function K α. Notice that we have assumed the functions qα to be constant in time. The goal of the linear forward problem is to compute the functions Ki α for any particular type of measurement di.

There have been numerous efforts to approximate the kernel functions K for time-distance travel times. These efforts will be described in Section 4.3.5. The kernel functions for ring diagrams have typically been approximated as constant within the area the ring is measured over, with a depth dependence derived from normal mode theory (see Section 4.2.3).

The linear forward problem for normal-mode frequencies is quite well known (see the upcoming Living Reviews paper on global helioseismology, Thompson (2005Jump To The Next Citation Point)). For example, Gough and Thompson (1990), Goldreich et al. (1991Jump To The Next Citation Point), and Dziembowski and Goode (2004Jump To The Next Citation Point) applied first-order perturbation theory to estimate the effect of magnetic fields on normal mode frequencies. This work may be helpful for computing the effects of magnetic fields on local helioseismic measurements.

In local helioseismology, the sensitivity functions K α may be computed using the Born approximation. The first Born approximation gives the lowest order approximation for estimating scattering amplitudes (Born, 1926); this method appears to have been first used by Strutt (1881). As explained below, the Born approximation is essentially an equivalent-source description of wave interaction. For small perturbations to a solar model, it is used to compute the first-order change in the wavefield. The first-order change in any particular helioseismic parameter can then be obtained from the first-order, as well as the zero-order, wavefields. To compute the first-order change in the wavefield we begin from Equation (7View Equation) and Equations (19View Equation, 20View Equation, 21View Equation). To first order we have

0 0 ℒ δξ = − δℒ ξ + δS. (29 )
This is the equation for the first Born approximation. The first-order correction to the observed wavefield is therefore given by
∫ 3 ′ ′ 3 i ′ ′{ j ′ ′ } 0 δ Φ(x, t) = − d r dt d rsdts 𝒢 (x,t;r ,t) δℒG (r ,t;rs,ts) iSj (rs,ts) ∫ 3 i + d rsdts 𝒢 (x, t;rs,ts)δSi(rs,ts), (30 )
where j {δℒG }i is the i− th component of the vector j δℒG and summation on repeated indices is implied. From Equation (30View Equation) we can predict the first-order change in any quantity in local helioseismology that will result from a small change in the model. We can look at, for example, first-order changes in the cross-covariance, travel times, power spectrum, or holographic measurements.

An alternative to the Born approximation is the Rytov approximation. In the Rytov approximation one computes the first order correction to the phase of the wavefield rather than the correction to the wavefield itself. For applications in the context of helioseismology, see Brüggen (2000) and Jensen and Pijpers (2003Jump To The Next Citation Point). The Born and Rytov approximations have been compared by Keller (1969).

  Go to previous page Go up Go to next page