We discuss current advances in forward and inverse modeling for local helioseismology. We report theoretical uniqueness results, in particular the Novikov-Agaltsov reconstruction algorithm, which is relevant to solving the non-linear inverse problem of time-distance helioseismology (finite amplitude pertubations to the medium). Numerical experiments were conducted to determine the number of frequencies required to reconstruct density and sound speed in the solar interior.

Local helioseismology has so far relied on semi-analytical methods to compute the spatial sensitivity of wave travel times to perturbations in the solar interior. These methods are cumbersome and lack flexibility. Here we propose a convenient framework for numerically solving the forward problem of time-distance helioseismology in the frequency domain. The fundamental quantity to be computed is the cross-covariance of the seismic wavefield. We choose sources of wave excitation that enable us to relate the cross-covariance of the oscillations to the Green's function in a straightforward manner. We illustrate the method by considering the 3D acoustic wave equation in an axisymmetric reference solar model, ignoring the effects of gravity on the waves. The symmetry of the background model around the rotation axis implies that the Green's function can be written as a sum of longitudinal Fourier modes, leading to a set of independent 2D problems. We use a high-order finite-element method to solve the 2D wave equation in frequency space. The computation is `embarrassingly parallel', with each frequency and each azimuthal order solved independently on a computer cluster. We compute travel-time sensitivity kernels in spherical geometry for flows, sound speed, and density perturbations under the first Born approximation. Convergence tests show that travel times can be computed with a numerical precision better than one millisecond, as required by the most precise travel-time measurements. The method presented here is computationally efficient and will be used to interpret travel-time measurements in order to infer, e.g., the large-scale meridional flow in the solar convection zone. It allows the implementation of (full-waveform) iterative inversions, whereby the axisymmetric background model is updated at each iteration.

In time-distance helioseismology, information about the solar interior is encoded in measurements of travel times between pairs of points on the solar surface. Travel times are deduced from the cross-covariance of the random wave field. Here we consider travel times and also products of travel times as observables. They contain information about e.g. the statistical properties of convection in the Sun. The basic assumption of the model is that noise is the result of the stochastic excitation of solar waves, a random process which is stationary and Gaussian. We generalize the existing noise model (Gizon and Birch 2004) by dropping the assumption of horizontal spatial homogeneity. Using a recurrence relation, we calculate the noise covariance matrices for the moments of order 4, 6, and 8 of the observed wave field, for the moments of order 2, 3 and 4 of the cross-covariance, and for the moments of order 2, 3 and 4 of the travel times. All noise covariance matrices depend only on the expectation value of the cross-covariance of the observed wave field. For products of travel times, the noise covariance matrix consists of three terms proportional to $1/T$, $1/T^2$, and $1/T^3$, where $T$ is the duration of the observations. For typical observation times of a few hours, the term proportional to $1/T^2$ dominates and $Cov[\tau_1 \tau_2, \tau_3 \tau_4] \approx Cov[\tau_1, \tau_3] Cov[\tau_2, \tau_4] + Cov[\tau_1, \tau_4] Cov[\tau_2, \tau_3]$, where the $\tau_i$ are arbitrary travel times. This result is confirmed for $p_1$ travel times by Monte Carlo simulations and comparisons with SDO/HMI observations. General and accurate formulae have been derived to model the noise covariance matrix of helioseismic travel times and products of travel times. These results could easily be generalized to other methods of local helioseismology, such as helioseismic holography and ring diagram analysis.

Inversions for local helioseismology are an important and necessary step for obtaining three-dimensional maps of various physical quantities in the solar interior. Frequently, the full inverse problems that one would like to solve prove intractable because of computational constraints. Due to the enormous seismic data sets that already exist and those forthcoming, this is a problem that needs to be addressed. To this end, we present a very efficient linear inversion algorithm for local helioseismology. It is based on a subtractive optimally localized averaging (SOLA) scheme in the Fourier domain, utilizing the horizontal-translation invariance of the sensitivity kernels. In Fourier space the problem decouples into many small problems, one for each horizontal wave vector. This multi-channel SOLA method is demonstrated for an example problem in time-distance helioseismology that is small enough to be solved both in real and Fourier space. We find that both approaches are successful in solving the inverse problem. However, the multi-channel SOLA algorithm is much faster and can easily be parallelized.