Ghost Images in Helioseismic Holography? Toy Models in a Uniform Medium

Helioseismic holography is a powerful technique used to probe the solar interior based on estimations of the 3D wavefield. The Porter–Bojarski holography, which is a well-established method used in acoustics to recover sources and scatterers in 3D, is also an estimation of the wavefield, and hence it has the potential of being applied to helioseismology. Here we present a proof-of-concept study, where we compare helioseismic holography and Porter–Bojarski holography under the assumption that the waves propagate in a homogeneous medium. We consider the problem of locating a point source of wave excitation inside a sphere. Under these assumptions, we find that the two imaging methods have the same capability of locating the source, with the exception that helioseismic holography suffers from “ghost images” (i.e. artificial peaks away from the source location). We conclude that Porter–Bojarski holography may improve the method currently used in helioseismology.


Introduction
Local helioseismology is a powerful tool used to probe the 3D interior of the Sun by exploiting the information contained within the acoustic and surface-gravity waves observed at the surface (see, e.g., Gizon and Birch, 2005;Gizon, Birch, and Spruit, 2010). Helioseismic holography is one branch of local helioseismology that aims at imaging the subsurface structure by estimating the wavefield inside the Sun Braun, 1997, 2000a). One significant achievement of helioseismic holography has been the detection of active regions on the far side of the Sun (far-side imaging: Lindsey and Braun, 2000b). The technique used in far-side imaging, known as phase-sensitive holography, has been validated with synthetic data (see, e.g., Hartlep et al., 2008;Birch et al., 2011;Braun, 2014), and it is extensively used in studying active regions in the near hemisphere (e.g. Braun and Birch, 2008;Braun, 2016). B D. Yang yangd@mps.mpg.de their respective sensitivity kernels. This allows us to refine the scope of the previous question; specifically, which method is more accurate at locating acoustic sources? Comparisons between these methods have been done in ocean acoustics, where Jackson and Dowling (1991) showed that both methods can locate the source. Hence, the authors concluded that the egression is the simplified version of the PB hologram, and chose to use the egression as their preferred method since it is easier to implement. In the case of helioseismology, however, the question of the optimal method has yet to be answered. This is the goal of this article. Additionally, preliminary work by Lindsey and Braun (2004) showed that helioseismic holography suffers from unintended mirror-like images ("ghost images") that are caused by the use of only a monopole source. An examination of ghost images in the egression and the PB hologram is also a focus of this study.
In this article, we present a proof-of-concept study, where we compare the sourcesensitivity kernels of helioseismic holography and PB holography by assuming waves are propagating in a homogeneous medium. Specifically, we examine which method is more accurate at locating acoustic sources, and hence at estimating the wavefield. Additionally, we examine the effect that the observational coverage area has on the kernels. This will provide an opportunity to improve the method currently used in helioseismology. This article is organized as follows: Section 2 states the derivation of source-sensitivity kernels of helioseismic holography (egression) and PB holography, Section 3 states the toy model used in this article, Section 4 compares the two methods with a discussion, and conclusions are provided in Section 5.

Source-Sensitivity Kernels
In this section, we present the relation of helioseismic holography and PB holography to acoustic sources. We work entirely in the temporal Fourier domain, using the convention where f (ω) is the Fourier transform of a given function F (t).
The egression [ LB ], as described by Lindsey and Braun (1997: LB), is one of the basic quantities used in helioseismic holography, where r denotes the focal point, A is the coverage of the wavefield (r , ω) at any point r on the solar surface, and G(r, r , ω) is the Green's function associated with a wave operator defined below, with the asterisk denoting the complex conjugate. For simplicity, we drop the ω within the function's arguments for the remainder of this study. Further definitions and explanations concerning the Green's function are given below. The PB hologram PB is defined by Devaney and Porter (1985), where ∂ n denotes the outward-normal derivative with respect to r , and Im G is the imaginary part of the Green's function.
In order to relate LB and PB to the acoustic sources, we first need to determine the impulse response function [G, the Green's function] of the wave equation. Here we assume that the wavefield [ ] is related to the sources through the application of a linear wave operator [L], where S(r) is the source function. For generality, we choose not to explicitly state L here. The Green's function is the impulse response of Equation 4, and it is defined as the solution to LG(r, r s ) = δ(r − r s ), where δ(r − r s ) is the Dirac delta function and r s is the location of the source. One property of the Green's function, which is crucial to deriving source-sensitivity kernels, is that it can be used to solve Equation 4 through Through expansion of in Equation 2 with the definition in Equation 6, the egression becomes LB and through a change in the order of integration, one obtains the definition for the source sensitivity kernel [K LB ] for the egression LB where We note that the 4π factor is included here such that K LB possesses a desired near-unitary amplitude for this study.
The same procedure is repeated for the derivation of the PB hologram by expanding in Equation 3 with Equation 6 and changing the order of integration: where K PB A is the source-sensitivity kernel of the PB hologram. A comparison of the egression and the PB hologram requires only the knowledge of their respective source-sensitivity kernels, whereas details of the source function are not needed. Therefore, we examine the source-sensitivity kernels of the two imaging methods in this study. We compare the two source kernels under simplifying assumptions about the medium in which the waves propagate. We note that in practice, the egression power | LB | 2 is used to estimate the location of acoustic sources, since the wavefield in the Sun is stochastically excited (see, e.g., Lindsey and Braun, 1997;Hanson, Donea, and Leka, 2015). Therefore, we also compare the squared modulus of the source kernels.

Toy Model: Waves in a Homogeneous Medium
With Equations 9 and 11 in hand, we require the computation of the Green's functions in order to determine the source-sensitivity kernels. In general, a Green's function can be obtained numerically for any given linear operator [L]. However, as stated in the introduction, we examine a homogeneous medium, and therefore the Green's function can be computed analytically.
We consider this homogeneous medium in R 3 space with a constant sound speed c = 10 5 m s −1 , and we adopt the linear wave operator, where k is the wavenumber, ω is the angular frequency and γ is the damping rate. The solution of Equation 5 with the above wave operator and free boundary condition is given by which is also known as the outgoing free-space Green's function (Born and Wolf, 1999). In this study, we set the damping rate [γ ] to be 0.1% of the angular frequency, and we use a Cartesian coordinate system (x, y, z) with its origin at the center of a sphere V with the radius R = 696 Mm. Current observational capabilities mean that we can only observe the wavefield on a fraction of the solar surface. To study the consequence of this limitation on observations, we examine both the case where the entire surface is observed and the case where only a fraction of the solar surface is observed. In these cases, we assume that the sources are located along the z-axis and that the coverage is symmetric with respect to the z-axis and is centered above the North Pole (0, 0, R ). Under these assumptions, K LB and K PB are axisymmetric about the z-axis.

Source-Sensitivity Kernels at 3 mHz
Here we examine the source-sensitivity kernels [K LB and K PB ] at a frequency of ω/2π = 3 mHz. Figure 1 shows 2D slices of both the real and imaginary parts of K LB and cK PB /ω through the z-axis, for the source locations at z s = 0.9R (panels a to h) and −0.9R (panels i to l) on the z-axis. The factor c/ω is added to K PB so that the kernel is dimensionless like the egression. The first row of panels shows the source-sensitivity kernels under the assumption that the entire surface is observed. The remaining rows show the kernels assuming a coverage of 60 degrees from the North Pole. Considering all panels, both Re[K LB ] and Re[K PB ] peak at the source, while in comparison, Im[K LB ] and Im[K PB ] are negligible if the source is located at the near side. These results demonstrate that Re[K LB ] and Re[K PB ] can both locate the source in either of these coverage geometries. In the case of sources located on the far side, all of the kernels have become less localized. While these kernels can locate the sources, we also see that the egression kernels (both the real and imaginary parts) have "ghost images" above the surface, while the PB holograms do not. These ghost images appear as peaks at points away from the source location. We note that we here also D. Yang observed ghost images below the surface in the egression, when the source is above the surface. This suggests that the egression cannot distinguish sources from below and above the surface, since one cannot differentiate sources and ghosts. The PB holograms do not suffer from this problem. Further explanations and discussions concerning the ghost images are given in Section 5.1.
For a more focused comparison of the kernels in Figure 1, Figure 2 shows 1D slices of K LB and K PB along the z-axis with the real parts shown in panels a, c, and e, and the imaginary parts in panels b, d, and f. Re[K LB ] and Re[K PB ] again have peaks at the source location. Here we see that despite the coverage geometry, Im[K LB ] and Im[K PB ] are always zero at the source location, with the peaks seen in Figure 1 surrounding the source location. This suggests that Im[K LB ] and Im[K PB ] cannot pinpoint the exact source location. In the case of sources located on the far side, K LB and K PB are both highly oscillatory and nonlocalized, and hence recovering the source location may be problematic with observations at a single frequency. Figure 3 shows 1D slices of K LB and K PB along a line that is perpendicular to the z-axis and passes through the source. Owing to axial symmetry, only half of the slice is plotted.
Here the imaginary parts of K LB and K PB are not shown since they are negligible compared to the real parts. Unlike the vertical slices, neither K LB nor K PB have ghost images, and they are less oscillatory when the source is located at the far side.

Kernels Averaged over Frequency
Specifically, the two methods are highly oscillatory and non-localized for the source located at the far side for observations made at a single frequency, and the egression suffers from ghost images for the source located at the near side. One possible solution to these problems is to average kernels over a number of frequencies, since the ghost images for the source located at the near side and the side-lobes for the source located at the far side may peak at different locations for different frequencies. Figure 4 shows 1D vertical slices of K LB and K PB averaged from 41 frequencies equally distributed from 1 to 5 mHz. A Gaussian weight function centered at 3 mHz and with a standard deviation of 1 mHz was applied for the averaging. From these results, it is clear that averaging the kernels over frequency reduces the amplitude of the ghost images for the near-side source and the side-lobes for the far-side source. The averaged kernels along the horizontal direction are similar to the kernels with a single frequency, and as such they are not shown here.
As mentioned in Section 2, the egression power, which is related to the source covariance via |K LB | 2 , has been used in observations as estimates of the acoustic sources. Therefore, the effect of averaging |K LB | 2 and |K PB | 2 over different frequencies is also of great interest. Figure 5 shows a comparison of |K LB | 2 and |K PB | 2 with or without averaging over different frequencies. Only the slices along the vertical direction are shown, as the difference between |K LB | 2 (|K PB | 2 ) from a single frequency at 3 mHz and averaged from 1 to 5 mHz along the horizontal direction is small. In the vertical direction, averaging |K LB | 2 and |K PB | 2 over different frequencies reduces the amplitude of the ghosts when the source is located at the near side, and improves the spatial resolution when the source is located at the far side.

Dependence of the Spatial Resolution on the Coverage
The results thus far have shown that both of the methods can locate the source, although the egression has the complication of ghost peaks. The question then arises of how well these methods resolve the sources with differing observational coverages. We define the spatial resolution the egression and the PB hologram as the full-width at half-maximum (FWHM) of |K LB | 2 and |K PB | 2 , respectively. Figures 2 and 3 show that both of the methods behave differently in the vertical and horizontal directions, and the FWHMs in these two directions are therefore considered separately. Additionally, to quantify the effect of averaging kernels over different frequencies, we considered the FWHM for kernels from a single frequency at 3 mHz and the value averaged over frequencies from 1 to 5 mHz. Figure 6 shows the FWHM of |K LB | 2 and |K PB | 2 as a function of the angle [θ h ], which defines the observational coverage (cap of area A). The FWHMs along the vertical (top row) and horizontal (bottom row) directions are considered in the case of kernels at 3 mHz (left column) and averaged from 1 to 5 mHz (averaged kernel, right column). At 3 mHz, the difference between the two imaging methods is small, which implies that either method has the same capability of resolving the source. Furthermore, the FWHM is close to the resolution limit despite the size of the coverage area when the source is located at the near side (z s = 0.9R ). When the source is located at the far-side (z s = −0.9R ), the resolution improves (FWHM decreases) with increasing coverage. When averaging over different frequencies, a clear improvement of the spatial resolution can be found along the vertical direction when the source is located on the far side, while the spatial resolution is almost the same as before for the other cases. Figure 4 1D slices of kernels averaged over frequencies using a Gaussian weight function centered at 3 mHz with a standard deviation of 1 mHz. Here we use a bar to denote the averaged quantity. The observational coverage A is depicted in the left panels (thick arcs). The frequency averaging reduces the amplitudes of the ghost images for the near-side source and the side-lobes for the source on the far side.

Ghost Images in the Egression
The appearance of ghost images in the egression can be understood by Huygens' principle, whereby each arbitrarily small section of the observed wavefield is regarded as a point source, and the egression as a superposition of the back-propagated (in time) waves generated from all of the point sources. Furthermore, each newly created wave is spherically symmetric with respect to its source location in a homogeneous medium, and thus it will propagate in all directions with the same behavior. This is the cause of the ghost images. A clear example of this is evident when the wavefield is recorded on a plane, where all of the newly created waves are symmetric with respect to the recording surface, and hence the egression will focus on both the source location and its counterpart on the other side of the surface (see Figure 7 for an example). When the wavefield is observed on a sphere, however, the newly created waves are no longer symmetric with respect to the surface, and the ghost images show a complicated diffraction pattern that is due to the interference among the newly created waves (Figure 1). This provides a simple explanation for the ghost images seen in the egression above the surface (see also Lindsey and Braun, 2004).
The PB hologram does not suffer from ghost images like the egression, since it includes not only a monopole source, but also a dipole source, which is not symmetric with respect D. Yang Figure 5 1D slices of |K LB | 2 and |K PB | 2 in a plane containing the z-axis. The left panels show a single frequency of 3 mHz. The right panels show averages over frequencies using a Gaussian weight function centered at 3 mHz with a standard deviation of 1 mHz. Averaging the kernels over frequencies reduces the amplitude of the ghosts when the source is located on the near side, and it improves the spatial resolution when the source is located on the far side.
to the source location. Additionally, the amplitudes of the monopole and dipole sources are chosen such that the PB hologram only focuses on the source location.
Future work should include a solar-like density stratification to confirm this simple explanation. The sharp drop in density at the solar surface leads to a reflection of the waves below 5.3 mHz, which is not captured in our toy model.
We also note that Lindsey and Braun (2005a,b) proposed that ghost images may explain the presence of phase anomalies observed around active regions in phase-sensitive holography. For further implications and discussions of the ghost images in helioseismic holography, we refer to Lindsey and Braun (2004).

Ingression and the PB Hologram
We did not consider the ingression in our analysis. The ingression is an equally important quantity used in helioseismic holography, which is an estimation of the location at which the wavefield converges by propagating the wavefield forward in time (Lindsey and Braun, 1997). So far, we have considered the wavefield as diverging away from the source. However, to compare the ingression and the PB hologram, a wavefield that converges from infinity to the source location is desired. Such a wavefield can be achieved by considering the wavefield diverging from the source as before, but with the reversed sign in time, i.e. The FWHM along two directions, horizontal and vertical, is shown in the case of a near-side source (z s = 0.9R , blue) and a far-side source (z s = −0.9R , red). The theoretical resolution limit of λ/2 is also shown with a horizontal black line. Additionally, the results for the kernels at a single frequency ω/2π = 3 mHz and the frequency-averaged kernels are shown in the left and right columns, respectively.

Figure 7
2D slice of the real (left panel) and imaginary (right panel) part of K LB through the z-axis when the wavefield is observed at the z = 0 plane. Here the coverage is a circle that is centered at the origin and has a radius R . The source is located along the z-axis at z s = −0.1R and is indicated by the focus of the crosshairs in each plot.
(r, −t). In the frequency domain, this time-reversal corresponds to taking the complex conjugate. Additionally, the wave number [k] in the wave equation is also conjugated and thus the wavefield decays when propagating toward the source location (Devaney, 2012). In this case, the ingression is and the PB hologram becomes PB A,− (r, ω) = A d 2 r * r , ω ∂ n Im G r, r , ω − Im G r, r , ω ∂ n * r , ω . (15) We can see that LB − and PB − are simply the complex conjugates of LB and PB . Since we have discussed the real and imaginary parts and the power of PB and LB separately in the results, those of LB ( PB ) will be the same as LB − ( PB − ). In particular, LB − will also have ghost images while PB − will not, and | LB − | 2 and | PB − | 2 will have the same spatial resolution when imaging acoustic sources.

Application to Stereoscopic Helioseismology
The results in Section 4.3 showed that the spatial resolution of the hologram on the far side increases as the coverage area increases. Thus the resolution has a fundamental limit when observing from a single vantage point. It has been suggested to combine observations from two or several vantage points to increase the observation coverage and therefore to improve the spatial resolution (and signal-to-noise ratio) of holography. Stereoscopic helioseismology is believed to be our best chance to probe the subsurface structure in the polar regions and the deep convection zone, which is crucial for understanding the 11-year solar cycle (see, e.g., Ruzmaikin and Lindsey, 2003). This conjecture, however, has not been studied in detail. Figure 8 shows the PB hologram on the solar surface with a Dirac delta source located 0.7 Mm below the surface at 270 • longitude along the Equator. Different coverage geometries are considered in the case of a single spacecraft (top row), two spacecraft in the Ecliptic (middle row), and two spacecraft with one in the Ecliptic and the other at 45 • inclination (bottom row). Here only the real parts of the PB hologram are shown, since the imaginary parts are negligible. The results show a clear improvement of the spatial resolution when a second spacecraft is added, whereas the FWHM along the Equator is about twice smaller than that of a single spacecraft. Additionally, the spatial resolution is increased along the great circle (black-dashed line) at the intersection of the plane that contains the center of the sphere and the two spacecraft.
Stereoscopic helioseismology might be implemented in future space missions such as Solar Orbiter and Solar Activity Far Side Investigation (see, e.g., Sekii et al., 2015 and references therein) together with observations collected from the ground (Global Oscillation Network Group) or from near-Earth orbit (Solar Dynamics Observatory). In particular, the Polarimetric and Helioseismic Imager onboard Solar Orbiter is to be launched soon and will provide high-resolution line-of-sight velocity and continuum intensity at the photosphere, which are suitable for helioseismic studies (Woch and Gizon, 2007;Müller et al., 2013;Löptien et al., 2015). The orbit of Solar Orbiter will have a period of 168 days during the nominal mission and reach a heliographic latitude of up to 25 • (35 • during an extended mission) (Müller et al., 2013). This means that Solar Orbiter will cover a wide range of spacecraft-Sun-Earth angles to test stereoscopic helioseismology.

Outlook
We found that helioseismic holography and PB holography are similar techniques, with the exception that the egression and the ingression suffer from ghost images. In principle, Here the source is located 0.7 Mm below the surface at 270 • longitude along the Equator, and the plots are shown after divided by the maximum value of Re[K PB ]. We note that Im[K PB ] is negligible, and as such is not shown here. We plot Re[K PB ] on the entire surface in the left column, where the boundary of the coverage is marked by a red curve and the point below the spacecraft by a cross. Additionally, a zoom of images around the source location is shown in the right column, where a circle centered above the source location with a diameter of the wavelength is added in each plot. The spatial resolution is clearly improved when a second spacecraft is added. Furthermore, a preferred direction that possesses higher spatial resolution is found along the great circle (black-dashed line) that intersects the points below the two spacecraft.
we could apply the PB holograms to phase-sensitive holography by replacing the egression and the ingression with the appropriate PB and PB − . Our toy model suggests that the PB holograms will improve current helioseismic holography since they do not suffer from ghost images. However, additional modeling work is needed. Future studies should consider random acoustic sources and scatterers. Furthermore, the computations must be carried out in a solar-like stratified background medium. Finally, in order to implement PB holography, a method for determining the normal derivative of the wavefield needs to be developed.

D. Yang
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.