### 4.4 Helioseismic holography

Helioseismic holography was introduced in detail by Lindsey and Braun (1990), although the basic concept was first suggested by Roddier (1975). For a recent theoretical introduction to helioseismic holography see Lindsey and Braun (2000a). The central idea in helioseismic holography is that the wavefield, e.g., the line-of-sight Doppler velocity observed at the solar surface, can be used to make an estimate of the wavefield at any location in the solar interior at any instant in time. In this sense, holography is much like seismic migration, a technique in geophysics that has been in use since the 1940s (see, e.g., Hagedoorn, 1954, and references therein). In migration, the wave equation is used to determine the wavefield in the interior of the earth at past times, given the wavefield observed at the surface (see, e.g., Claerbout, 1985).

In Section 4.4.1 we introduce the basic computations of helioseismic holography, the ingression and egression. In Section 4.4.2 we review the various calculations that have been used to obtain the Green’s functions used in holography. In Section 4.4.3 we introduce the notion of the local control correlation. Acoustic power measurements are described in Section 4.4.4. Phase-sensitive holography is introduced in Section 4.4.5. Section 4.4.6 describes the technique of far-side imaging. In Section 4.4.7 we describe acoustic imaging (Chang et al., 1997), a technique closely related to holography.

#### 4.4.1 Ingression and egression

The ingression and egression are attempts to the answer the following question: Given the wavefield observed at the solar surface, what is the best estimate of the wavefield at some point in the solar interior assuming that the observed wavefield resulted entirely from waves diverging from that point (for the egression) or waves converging towards that point (for the ingression)?

The development of holography has, historically, been motived by analogy with optics (e.g., Lindsey and Braun, 2000a) and as a result the holography literature often employs optical terminology. The target point in the solar interior at which we attempt to estimate the wavefield is termed the “focus point”. In practice, only data in a restricted area on the surface above the focus point is used to compute the egression and ingression. This region is called the “pupil”. The choice of pupil geometry depends on the particular application and will be described when we discuss particular applications of holography.

As described by Lindsey and Braun (2000a), a mathematical motivation for the ingression and egression is the Kirchoff integral solution to the wave equation (e.g., Jackson, 1975). Suppose that a scalar field solves the equation

The Green’s function, , for this problem is defined as the solution to
where is the 3D Dirac delta function. The Kirchoff integral equation for is
The integration variable ranges over any surface enclosing the locations . For a derivation see for example the introduction in Jackson (1975), noting that the Kirchoff integral equation is a special case of Green’s second identity. The normal derivative, outwards, is denoted by . Equation (69) is interesting as it says that given the value of the field on the surface of a region we can determine the value of the field anywhere inside this region. Holography is an attempt to do just that, determine the wavefield in the solar interior given the wavefield on the solar surface. As discussed in detail by Lindsey and Braun (2004), helioseismic holography is much more complicated than solving the simple scalar wave Equation (67) addressed by the Kirchoff integral.

Motivated by the Kirchoff identity (Equation 69), the egression and ingression are defined as

The focus point is located at position in the solar interior. The integration variable denotes horizontal position on the solar surface and ranges over the surface region denoted by the pupil , typically an annular region, or a sector of an annular region, centered on the horizontal position of the focus point, . For examples of pupils see Lindsey and Braun (2000a). The wavefield on the surface is and is temporal angular frequency. Recall that the egression, , is an attempt to estimate the wavefield at the location and frequency based on the waves that are seen, in the pupil at the surface, diverging away . Likewise, the ingression, is an attempt to estimate the wavefield at the location and frequency based on the waves seen, within the pupil, converging towards the focus point. The holography Green’s functions, which depend on a horizontal displacement , a focus depth , and a temporal frequency are denoted by . In the time domain, the causal Green’s function , obtained by temporal inverse Fourier transform from , is zero for . In the time domain, the anti-causal Green’s function is zero for . The time-domain causal Green’s function and anti-causal Green’s function are related by . The holography Green’s function are not defined as the solution to any particular set of equations, but rather we use the symbol to denote whatever function is used in practice in the computation of the ingression and egression. There are a number of choices that have been used for these Green’s functions (see Section 4.4.2 for details).

An intuitive way to interpret Equation (70) is to view the Green’s functions as propagators. In this way, the egression can be seen as the result of using the anti-causal Green’s function to propagate the surface wavefield backwards in time. Likewise, the ingression can be seen as the result of using the causal Green’s function to propagate the observed surface wavefield forwards in time.

The egression and ingression are the essential quantities in helioseismic holography. There are many particular ways in which these quantities can be combined to learn about the solar interior. The main techniques are control-correlations (Section 4.4.3), acoustic power holography (Section 4.4.4), phase-sensitive holography (Section 4.4.5), and far-side imaging (Section 4.4.6) which is a special case of phase-sensitive holography. Of central importance to all of these methods is the choice of the holography Green’s function , which we now describe.

#### 4.4.2 Holography Green’s functions

The two approaches that have been used to construct the holography Green’s functions, , involve ray theory and wave theory.

##### 4.4.2.1 Ray theory.
Lindsey and Braun (2000a), motived by ray theory, prescribed the holography Green’s function as

where and is the horizontal distance. The role of the function () is to propagate the observed surface wavefield forwards (backwards) in time and into the interior of the solar model. Theses functions can also be seen as the “wavefield” response at a horizontal distance , depth , and time to a surface source located at the origin at time . The observed surface wavefield can be thought of as a source, in the spirit of Huygen’s principle. The index refers to the number of skips the ray path has taken off the solar surface. Each function is the amplitude of a single skip component of the Green’s function and can be estimated from simple ray theory. In particular, in ray theory the amplitude functions are determined by the conservation of energy and the increase, with distance along the ray path, of the cross-sectional area of a ray bundle emerging from the source location. In Equation (71), the functions are the ray-theoretical travel times along the -skip ray paths connecting a point on the surface with a point at depth and a horizontal distance away from the first point. The sum over skips in Equation (71) is convenient as it can easily be truncated to estimate, for example, just the “one-skip” Green’s function, i.e., the Green’s function that contains only contributions from waves that have traveled once into the solar interior. The holography Green’s functions in the time domain (Equation 71) must then be Fourier transformed to obtain referred to in Equation (70). Notice that there are a great many simplifications involved in writing the Green’s function in Equation (71). The substantial frequency dependence of the travel time (e.g., Jefferies et al., 1994) has been ignored. There is no account taken of damping, which is quite significant for high-degree modes (e.g., Duvall Jr et al., 1998). Furthermore, as the ray approximation requires the wavelength to be much smaller than the length scale of the variation of the background medium, the ray approximation is not expected to be valid near the solar surface. Further issues include the neglect of the buoyancy and cut-off frequencies in the calculation of , and also in the computation of the ray paths themselves (Lindsey and Braun 2000a). Recent work by Lindsey and Braun (2004) suggests that after empirical dispersion corrections have been applied (see Section 4.4.3), the results of ray theory calculations are quite similar to the results of wave theory calculations (see Figure 25).

##### 4.4.2.2 Wave theory.
An alternative to the ray calculation is to choose, somewhat arbitrarily, the holography Green’s functions to solve a wave equation for the vertical displacement of a fluid element. This was carried out by Lindsey and Braun (2004). The computation proceeds most simply in the 3D Fourier domain (- domain). For each horizontal wavenumber and frequency component the Green’s function can be found by solving a one-dimensional boundary value problem. This can be seen by considering the equations of motion, written in the Fourier domain (Lindsey and Braun, 2004),

Here and are the pressure perturbation and vertical displacement associated with the Green’s function. The acceleration due to gravity is and the background density and sound speed are and , respectively. Equations (72, 73) can be used to find the pressure and vertical displacement given two boundary conditions. The data, however, provide us with only one upper boundary condition. It is for this reason that the assumptions that all of the observed waves are up-going waves at the surface, when computing the egression, and down-going waves at the surface, when computing the ingression, are made (Lindsey and Braun, 2004). In this way the second boundary condition can be generated from the first. Notice that Equations (72, 73) do not include the effect of damping.

#### 4.4.3 Local control correlations

The local control correlation is the correlation between the observed signal at a particular point on the solar surface and the holographic reconstruction of the signal at that same point computed from the data in a surrounding region. In particular the local control correlations are defined as (see, e.g., Lindsey and Braun, 2004)

Here the depth denotes the coordinate of the solar surface. The angle brackets denote the average over a frequency band of width centered at frequency . Notice that we have not specified the pupil which goes into the ingression/egression calculations. The pupil that is used depends on the particular situation. Figure 26 shows quiet Sun local control correlations, measured using both ray theory and wave theory Green’s functions.

If the egression/ingressions were perfect reconstructions of the observed signal, then the phases of the would both, on average, be zero in the quiet Sun. In practice, this is not the case (Lindsey and Braun, 2004). The phases of the holography Green’s functions are typically corrected so that the phase of the average quiet Sun local control correlation is zero (e.g., Lindsey and Braun, 2000a). We refer to this correction as the empirical dispersion correction.

#### 4.4.4 Acoustic power holography

The goal of acoustic power holography is to estimate the amount of wave power emitted from a particular region, either at particular time or a particular frequency. The estimate of the power emitted at a particular frequency, , at horizontal location and depth is

and the estimate of the power emitted at a particular time is
Here is the temporal Fourier transform of . Notice that these estimates depend on the estimate of the Green’s function and the pupil that are used to compute . For further discussions of acoustic power holography see, for example, Lindsey and Braun (1997) and Donea et al. (1999). Acoustic power holography has been used in studies of wave excitation by flares (e.g., Donea et al., 1999) and around active regions (e.g., Donea et al., 2000).

#### 4.4.5 Phase-sensitivity holography

The basic computation in phase-sensitive holography is the ingression-egression correlation (see, e.g., Lindsey and Braun, 2000a2004),

where and are two pupils and the angle brackets denote averaging over a frequency range of width centered on frequency . Phase-sensitive holography has been used to look for sound-speed perturbations and mass flows. In both cases the geometry is the quadrant geometry that is also used in time-distance measurements. In particular the pupil is chosen to be a quarter of an annulus, let us call it for “left”, and is chosen to be the quarter annulus located 180° away, call it for “right”. The symmetric phase is defined as
Here the operator returns the phase. The phase is mostly sensitive to sound-speed; this can be seen from the symmetry of the definition. In particular, does not change when the pupils and are interchanged. Notice that for a horizontally uniform sound-speed perturbation . On the other hand, a uniform horizontal flow gives . Thus is zero for a horizontally uniform flow and non-zero for a horizontally uniform sound-speed perturbation. The anti-symmetric phase is defined as
The symmetry of is such that changes sign under interchange of the pupils and and, thus, the anti-symmetric phase is sensitive to horizontal flows. The phases and are analogous to the mean and difference travel times commonly employed in time-distance helioseismology.

#### 4.4.6 Far-side imaging

Far-side imaging is a special case of phase-sensitive holography (e.g., Lindsey and Braun, 2000bBraun and Lindsey, 2001). The idea is to use the wavefield on the visible disk to learn about active regions on the far-side of the Sun. Figure 27 shows the geometry that was used by Braun and Lindsey (2001). In order to obtain full coverage of the far-side of the Sun, two different geometries were employed. For regions near the antipode of the center of the visible disk a two-skip geometry and two-skip Green’s functions were employed (panel (a) of Figure 27). For focus positions near the limb, the ingression/egression were computed using a single-skip pupil and corresponding Green’s functions, and then correlated with the egression/ingression computed using a three-skip pupil and Green’s functions (the geometry is shown in panel (b) of Figure 27).

#### 4.4.7 Acoustic imaging

Acoustic imaging was first introduced by Chang et al. (1997); for a recent review see Chou et al. (2003). We have included acoustic imaging in the section on helioseismic holography as the definitions and philosophical motivation for the two techniques are quite similar.

The central computations in acoustic imaging (see, e.g., Chou et al., 19992003) are wavefield reconstructions in the solar interior, and , defined by

The focus position is . The function denotes the azimuthal average, around , of the surface wavefield measured a horizontal distance away from . The quantity represents the ray-theoretical travel times from the subsurface focus point at to surface points a horizontal distance away from the focus point. For a given , the distance satisfies the time-distance relation established between the focus point and a surface point. This relation, , is computed from a standard solar model, using the ray approximation (see Figure 28). The function is a smooth weight function of (or ), explicitly given by Chou et al. (1999). The sum in Equation (81) involves the observed wavefield inside an annulus with inner and outer radii specified by the range .

In some sense, acoustic imaging is a special case of holography. Here, the pupil is an annulus with inner radius and outer radius given by the time-distance relation . The equivalent holography Green’s function is given by . The signal corresponds to the egression, i.e., is the signal reconstructed from the observations of waves diverging from the focal point, while corresponds to the ingression, the signal estimated from the waves seen converging towards the focus position.

As with holography, both the amplitudes and phases of are used to learn about the solar interior. The square of the amplitude of is an estimate of the power contained in the waves seen diverging from the focus position (Chang et al., 1997). In the terminology of holography, the squared of the amplitude of is called the “egression power.” Chen et al. (1998) introduced the correlation

and then measured phase and group times by fitting a Gaussian wavelet to at fixed target point . Changes in the phase between and result in changes in and thus shift the travel times. A phase shift between the two reconstructions, and , is evidence for local changes in sound speed (Chen et al., 1998). Chou and Duvall Jr (2000) discuss the relationship between time-distance travel times and the travel times measured by acoustic imaging.

The forward problem that has received the most attention in acoustic imaging is the dependence of travel times on changes in the sound speed. Chou and Sun (2001) used the ray approximation to estimate the sensitivity of acoustic imaging phase travel times to changes in sound speed. The results showed that in general the horizontal resolution is greater than the vertical resolution and that the resolution decreases with increasing focus depth.

The inverse problem of determining the sound speed from a given set of travel times has been studied in the ray approximation as well. Sun and Chou (2002) tested ray-theory RLS inversions of phase-time measurements on artificial data. These tests showed that the RLS inversions were, for good choices of the regularization parameter, capable of reconstructing the model used to generate the artificial data. Inversions of travel times for sound speed were also discussed in detail by Chou and Sun (2001).