Full waveform inversion of solar interior flows

Research output: Contribution to journalArticle

Abstract

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 2 norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p 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).

Original languageEnglish (US)
Article number23
JournalAstrophysical Journal
Volume797
Issue number1
DOIs
StatePublished - Dec 10 2014

Fingerprint

solar interior
waveforms
inversions
seismic noise
helioseismology
stratification
inference
norms
flow velocity
travel time
travel
inversion
sun
spatial resolution
education
deviation
formulations
monitoring
profiles
prediction

Keywords

  • hydrodynamics
  • Sun: helioseismology
  • Sun: interior
  • Sun: oscillations
  • waves

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Full waveform inversion of solar interior flows. / Hanasoge, Shravan.

In: Astrophysical Journal, Vol. 797, No. 1, 23, 10.12.2014.

Research output: Contribution to journalArticle

@article{6ccadda6eec947b29071468c6da93a97,
title = "Full waveform inversion of solar interior flows",
abstract = "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 2 norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p 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).",
keywords = "hydrodynamics, Sun: helioseismology, Sun: interior, Sun: oscillations, waves",
author = "Shravan Hanasoge",
year = "2014",
month = "12",
day = "10",
doi = "10.1088/0004-637X/797/1/23",
language = "English (US)",
volume = "797",
journal = "Astrophysical Journal",
issn = "0004-637X",
publisher = "IOP Publishing Ltd.",
number = "1",

}

TY - JOUR

T1 - Full waveform inversion of solar interior flows

AU - Hanasoge, Shravan

PY - 2014/12/10

Y1 - 2014/12/10

N2 - 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 2 norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p 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).

AB - 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 2 norm deviations in travel times between prediction and observation, as measured using short-distance filtered f and p 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).

KW - hydrodynamics

KW - Sun: helioseismology

KW - Sun: interior

KW - Sun: oscillations

KW - waves

UR - http://www.scopus.com/inward/record.url?scp=84914094631&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84914094631&partnerID=8YFLogxK

U2 - 10.1088/0004-637X/797/1/23

DO - 10.1088/0004-637X/797/1/23

M3 - Article

AN - SCOPUS:84914094631

VL - 797

JO - Astrophysical Journal

JF - Astrophysical Journal

SN - 0004-637X

IS - 1

M1 - 23

ER -