FULL WAVEFORM INVERSION OF SOLAR INTERIOR FLOWS

• Published in: The Astrophysical journal Vol. 797; no. 1; pp. 23 - 9
• Main Author: Hanasoge, Shravan M.
• Format: Journal Article
• Language: English
• Published: United States 10.12.2014
• Subjects: ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; Categories; CONVERGENCE; DENSITY; ERRORS; FILTERS; FLOW MODELS; FORECASTING; HYDRODYNAMICS; Inversions; ITERATIVE METHODS; Mathematical models; MONTE CARLO METHOD; Noise; OSCILLATIONS; PROBABILISTIC ESTIMATION; Searching; SOLAR GRANULATION; SPATIAL RESOLUTION; STAR EVOLUTION; STRATIFICATION; SUN; Uncertainty; VELOCITY; Waveforms
• ISSN: 1538-4357; 0004-637X; 1538-4357
• DOI: 10.1088/0004-637X/797/1/23

The inference of flows of material in the interior of the Sun is a subject of major interest in helioseismology. Here, we apply techniques of full waveform inversion (FWI) to synthetic data to test flow inversions. In this idealized setup, we do not model seismic realization noise, training the focus entirely on the problem of whether a chosen supergranulation flow model can be seismically recovered. We define the misfit functional as a sum of L sub(2) norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p sub(1) and large-distance unfiltered p modes. FWI allows for the introduction of measurements of choice and iteratively improving the background model, while monitoring the evolution of the misfit in all desired categories. Although the misfit is seen to uniformly reduce in all categories, convergence to the true model is very slow, possibly because it is trapped in a local minimum. The primary source of error is inaccurate depth localization, which, due to density stratification, leads to wrong ratios of horizontal and vertical flow velocities ("cross talk"). In the present formulation, the lack of sufficient temporal frequency and spatial resolution makes it difficult to accurately localize flow profiles at depth. We therefore suggest that the most efficient way to discover the global minimum is to perform a probabilistic forward search, involving calculating the misfit associated with a broad range of models (generated, for instance, by a Monte Carlo algorithm) and locating the deepest minimum. Such techniques possess the added advantage of being able to quantify model uncertainty as well as realization noise (data uncertainty).