Computational helioseismology in the frequency domain: acoustic waves in axisymmetric solar models with flows

Abstract

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.

Accepted on 3 Nov 2016 for publication in Astronomy and Astrophysics

References in corpus (6)

• ADIPLS -- the Aarhus adiabatic oscillation package
• Helioseismology of Sunspots: Confronting Observations with Three-Dimensional MHD Simulations of Wave Propagation
• Linear Sensitivity of Helioseismic Travel Times to Local Flows
• High-resolution mapping of flows in the solar interior: Fully consistent OLA inversion of helioseismic travel times
• Time-distance helioseismology: Sensitivity of f-mode travel times to flows
• Generalization of the noise model for time-distance helioseismology

Cited by in corpus (17)

• Waves in the lower solar atmosphere: the dawn of next-generation solar telescopes
• A Comprehensive Method to Measure Solar Meridional Circulation and Center-to-Limb Effect Using Time-Distance Helioseismology
• Inversions for Deep Solar Meridional Flow Using Spherical Born Kernels
• Imaging individual active regions on the Sun's far side with improved helioseismic holography
• Finite-frequency sensitivity kernels in spherical geometry for time-distance helioseismology
• Atmospheric radiation boundary conditions for high frequency waves in time-distance helioseismology
• Direct assessment of SDO/HMI helioseismology of active regions on the Sun's far side using SO/PHI magnetograms
• Predicting frequency changes of global-scale solar Rossby modes due to solar cycle changes in internal rotation
• Helioseismic Modeling of Background Flows
• Contribution of flows around active regions to the north-south helioseismic travel-time measurements
• An Optimized Ly$\alpha$ Forest Inversion Tool Based on a Quantitative Comparison of Existing Reconstruction Methods
• Quantitative passive imaging by iterative holography: The example of helioseismic holography
• A general formulation for computing spherical helioseismic sensitivity kernels while incorporating systematical effects
• Computing Helioseismic Sensitivity Kernels for the Sun's Large-Scale Internal Flows using Global-Scale Wave-Propagation Simulations
• Helioseismic finite-frequency sensitivity kernels for flows in spherical geometry including systematic effects
• C2 representations of the solar background coefficients for the model S-AtmoI
• Numerical evaluation of time-distance helioseismic sensitivity kernels in spherical geometry