Validated helioseismic inversions for 3-D vector flows

According to time-distance helioseismology, information about internal fluid motions is encoded in the travel times of solar waves. The inverse problem consists of inferring 3-D vector flows from a set of travel-time measurements. Here we investigate the potential of time-distance helioseismology to infer 3-D convective velocities in the near-surface layers of the Sun. We developed a new Subtractive Optimally Localised Averaging (SOLA) code suitable for pipeline pseudo-automatic processing. Compared to its predecessor, the code was improved by accounting for additional constraints in order to get the right answer within a given noise level. The main aim of this study is to validate results obtained by our inversion code. We simulate travel-time maps using a snapshot from a numerical simulation of solar convective flows, realistic Born travel-time sensitivity kernels, and a realistic model of travel-time noise. These synthetic travel times are inverted for flows and the results compared with the known input flow field. Additional constraints are implemented in the inversion: cross-talk minimization between flow components and spatial localization of inversion coefficients. Using modes f, p1 through p4, we show that horizontal convective flow velocities can be inferred without bias, at a signal-to-noise ratio greater than one in the top 3.5 Mm, provided that observations span at least four days. The vertical component of velocity (v_z), if it were to be weak, is more difficult to infer and is seriously affected by cross-talk from horizontal velocity components. We emphasise that this cross-talk must be explicitly minimised in order to retrieve v_z in the top 1 Mm. We also show that statistical averaging over many different areas of the Sun allows for reliably measuring of average properties of all three flow components in the top 5.5 Mm of the convection zone.Comment: 14 pages main paper, 9 pages electronic supplement, 28 figures. Accepted for publication in Astronomy & Astrophysic

Time-distance helioseismology: Sensitivity of f-mode travel times to flows

Time-distance helioseismology has shown that f-mode travel times contain information about horizontal flows in the Sun. The purpose of this study is to provide a simple interpretation of these travel times. We study the interaction of surface-gravity waves with horizontal flows in an incompressible, plane-parallel solar atmosphere. We show that for uniform flows less than roughly 250 m s$^{-1}$, the travel-time shifts are linear in the flow amplitude. For stronger flows, perturbation theory up to third order is needed to model waveforms. The case of small-amplitude spatially-varying flows is treated using the first-order Born approximation. We derive two-dimensional Fr\'{e}chet kernels that give the sensitivity of travel-time shifts to local flows. We show that the effect of flows on travel times depends on wave damping and on the direction from which the observations are made. The main physical effect is the advection of the waves by the flow rather than the advection of wave sources or the effect of flows on wave damping. We compare the two-dimensional sensitivity kernels with simplified three-dimensional kernels that only account for wave advection and assume a vertical line of sight. We find that the three-dimensional f-mode kernels approximately separate in the horizontal and vertical coordinates, with the horizontal variations given by the simplified two-dimensional kernels. This consistency between quite different models gives us confidence in the usefulness of these kernels for interpreting quiet-Sun observations.Comment: 34 pages, accepted to Astrophysical Journa

A procedure for the inversion of f-mode travel times for solar flows

We perform a two-dimensional inversion of f-mode travel times to determine near-surface solar flows. The inversion is based on optimally localized averaging of travel times. We use finite-wavelength travel-time sensitivity functions and a realistic model of the data errors. We find that it is possible to obtain a spatial resolution of 2 Mm. The error in the resulting flow estimate ultimately depends on the observation time and the number of travel distances used in the inversion.Comment: 8 pages, 9 figure

Local helioseismology and the active Sun

The goal of local helioseismology is to elicit three-dimensional information about the sub-surface (or far-side) structure and dynamics of the Sun from observations of the helioseismic wave field at the surface. The physical quantities of interest include flows, sound-speed deviations and magnetic fields. However, strong surface magnetic fields induce large perturbations to the waves making inversions difficult to interpret. The purpose of this paper is to outline the methods of analysis used in local helioseismology, review discoveries associated with the magnetic Sun made using local helioseismology from the past three years, and highlight the efforts towards imaging the interior in the presence of strong magnetic fields.Comment: 6 pages, 4th HELAS International Conference, Lanzarote, Spain, 1-5 February 201

Interpretation of Helioseismic Travel Times - Sensitivity to Sound Speed, Pressure, Density, and Flows

Time-distance helioseismology uses cross-covariances of wave motions on the solar surface to determine the travel times of wave packets moving from one surface location to another. We review the methodology to interpret travel-time measurements in terms of small, localized perturbations to a horizontally homogeneous reference solar model. Using the first Born approximation, we derive and compute 3D travel-time sensitivity (Fr\'echet) kernels for perturbations in sound-speed, density, pressure, and vector flows. While kernels for sound speed and flows had been computed previously, here we extend the calculation to kernels for density and pressure, hence providing a complete description of the effects of solar dynamics and structure on travel times. We treat three thermodynamic quantities as independent and do not assume hydrostatic equilibrium. We present a convenient approach to computing damped Green's functions using a normal-mode summation. The Green's function must be computed on a wavenumber grid that has sufficient resolution to resolve the longest lived modes. The typical kernel calculations used in this paper are computer intensive and require on the order of 600 CPU hours per kernel. Kernels are validated by computing the travel-time perturbation that results from horizontally-invariant perturbations using two independent approaches. At fixed sound-speed, the density and pressure kernels are approximately related through a negative multiplicative factor, therefore implying that perturbations in density and pressure are difficult to disentangle. Mean travel-times are not only sensitive to sound-speed, density and pressure perturbations, but also to flows, especially vertical flows. Accurate sensitivity kernels are needed to interpret complex flow patterns such as convection

F-mode sensitivity kernels for flows

We compute f-mode sensitivity kernels for flows. Using a two-dimensional model, the scattered wavefield is calculated in the first Born approximation. We test the correctness of the kernels by comparing an exact solution (constant flow), a solution linearized in the flow, and the total integral of the kernel. In practice, the linear approximation is acceptable for flows as large as about 400 m/s.Comment: 4 pages, 3 figures. Proceedings of SOHO18/GONG 2006/HELAS I. Beyond the Spherical Sun: A new era of helio- and asteroseismology. Sheffield, England. August, 200

Propagating Linear Waves in Convectively Unstable Stellar Models: a Perturbative Approach

Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization, however, affects the result. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field, while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within $1~\mu$Hz for low-degree modes with frequencies near $3~\mu$Hz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection.Comment: 14 pages. Published online in Solar Physics, available at http://link.springer.com/article/10.1007/s11207-013-0457-

Seismic Sounding of Convection in the Sun

Our Sun, primarily composed of ionized hydrogen and helium, has a surface temperature of 5777~K and a radius $R_\odot \approx 696,000$ km. In the outer $R_\odot/3$, energy transport is accomplished primarily by convection. Using typical convective velocities $u\sim100\,\rm{m\,s^{-1}}$ and kinematic viscosities of order $10^{-4}$ m$^{2}$s$^{-1}$, we obtain a Reynolds number $Re \sim 10^{14}$. Convection is thus turbulent, causing a vast range of scales to be excited. The Prandtl number, $Pr$, of the convecting fluid is very low, of order $10^{-7}$\,--\,$10^{-4}$, so that the Rayleigh number ($\sim Re^2 Pr$) is on the order of $10^{21}\,-\,10^{24}$. Solar convection thus lies in extraordinary regime of dynamical parameters, highly untypical of fluid flows on Earth. Convective processes in the Sun drive global fluid circulations and magnetic fields, which in turn affect its visible outer layers ("solar activity") and, more broadly, the heliosphere ("space weather"). The precise determination of the depth of solar convection zone, departures from adiabaticity of the temperature gradient, and the internal rotation rate as a function of latitude and depth are among the seminal contributions of helioseismology towards understanding convection in the Sun. Contemporary helioseismology, which is focused on inferring the properties of three-dimensional convective features, suggests that transport velocities are substantially smaller than theoretical predictions. Furthermore, helioseismology provides important constraints on the anisotropic Reynolds stresses that control the global dynamics of the solar convection zone. This review discusses the state of our understanding of convection in the Sun, with a focus on helioseismic diagnostics. We present our considerations with the interests of fluid dynamicists in mind.Comment: 29 pages, 12 figures, in review, Annual Reviews of Fluid Mechanic

German Science Center for the Solar Dynamics Observatory

A data and computation center for helioseismology has been set up at the Max Planck Institute for Solar System Research in Germany to prepare for the SDO mission. Here we present the system infrastructure and the scientific aims of this project, which is funded through grants from the German Aerospace Center and the European Research Council

The art of fitting p-mode spectra: Part I. Maximum Likelihood Estimation

In this article we present our state of the art of fitting helioseismic p-mode spectra. We give a step by step recipe for fitting the spectra: statistics of the spectra both for spatially unresolved and resolved data, the use of Maximum Likelihood estimates, the statistics of the p-mode parameters, the use of Monte-Carlo simulation and the significance of fitted parameters. The recipe is applied to synthetic low-resolution data, similar to those of the LOI, using Monte-Carlo simulations. For such spatially resolved data, the statistics of the Fourier spectrum is assumed to be a multi-normal distribution; the statistics of the power spectrum is \emph{not} a $\chi^{2}$ with 2 degrees of freedom. Results for $l=1$ shows that all parameters describing the p modes can be obtained without bias and with minimum variance provided that the leakage matrix is known. Systematic errors due to an imperfect knowledge of the leakage matrix are derived for all the p-mode parameters.Comment: 13 pages, ps file gzipped. Submitted to A&