The parabolic wave equation in local helioseismology

In recent years methods of time-distance helioseismology have been used to produce maps of local flows in the surface layers of the Sun. Usually, these studies rely on ray theory to describe the propagation of sound waves. Ray theory, however, is a poor approximation of the acoustic wavefield near the surface of the Sun. In particular, it is inappropriate for the study of scattering and diffraction by inhomogeneities. But an exact solution of the acoustic wave equation in the Sun is not trivial. In this paper I present an approximation to the full wave equation, which transforms it into a parabolic equation. The parabolic equation is commonly used in ocean acoustics and geoseismology because it is much simpler to solve numerically. Here I discuss the parabolic approximation, its limitations and potential applications in helioseismology. Finally, I present some numerical results to demonstrate the capabilities of this method.