Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology

J. Jackiewicz, A. C. Birch, Laurent Gizon, Shravan Hanasoge, T. Hohage, J. B. Ruffio, M. Švanda

Research output: Contribution to journalArticle

Abstract

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 multichannel 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 multichannel SOLA algorithm is much faster and can easily be parallelized.

Original languageEnglish (US)
Pages (from-to)19-33
Number of pages15
JournalSolar Physics
Volume276
Issue number1-2
DOIs
StatePublished - Feb 1 2012

Fingerprint

helioseismology
inverse problem
inversions
solar interior
seismic data
invariance
sensitivity
inversion
method

Keywords

  • Helioseismology
  • Inverse modeling

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

Jackiewicz, J., Birch, A. C., Gizon, L., Hanasoge, S., Hohage, T., Ruffio, J. B., & Švanda, M. (2012). Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology. Solar Physics, 276(1-2), 19-33. https://doi.org/10.1007/s11207-011-9873-8

Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology. / Jackiewicz, J.; Birch, A. C.; Gizon, Laurent; Hanasoge, Shravan; Hohage, T.; Ruffio, J. B.; Švanda, M.

In: Solar Physics, Vol. 276, No. 1-2, 01.02.2012, p. 19-33.

Research output: Contribution to journalArticle

Jackiewicz, J, Birch, AC, Gizon, L, Hanasoge, S, Hohage, T, Ruffio, JB & Švanda, M 2012, 'Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology', Solar Physics, vol. 276, no. 1-2, pp. 19-33. https://doi.org/10.1007/s11207-011-9873-8
Jackiewicz, J. ; Birch, A. C. ; Gizon, Laurent ; Hanasoge, Shravan ; Hohage, T. ; Ruffio, J. B. ; Švanda, M. / Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology. In: Solar Physics. 2012 ; Vol. 276, No. 1-2. pp. 19-33.
@article{439b3d0c88b44e78921221a584f885c5,
title = "Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology",
abstract = "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 multichannel 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 multichannel SOLA algorithm is much faster and can easily be parallelized.",
keywords = "Helioseismology, Inverse modeling",
author = "J. Jackiewicz and Birch, {A. C.} and Laurent Gizon and Shravan Hanasoge and T. Hohage and Ruffio, {J. B.} and M. Švanda",
year = "2012",
month = "2",
day = "1",
doi = "10.1007/s11207-011-9873-8",
language = "English (US)",
volume = "276",
pages = "19--33",
journal = "Solar Physics",
issn = "0038-0938",
publisher = "Springer Netherlands",
number = "1-2",

}

TY - JOUR

T1 - Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology

AU - Jackiewicz, J.

AU - Birch, A. C.

AU - Gizon, Laurent

AU - Hanasoge, Shravan

AU - Hohage, T.

AU - Ruffio, J. B.

AU - Švanda, M.

PY - 2012/2/1

Y1 - 2012/2/1

N2 - 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 multichannel 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 multichannel SOLA algorithm is much faster and can easily be parallelized.

AB - 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 multichannel 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 multichannel SOLA algorithm is much faster and can easily be parallelized.

KW - Helioseismology

KW - Inverse modeling

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

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

U2 - 10.1007/s11207-011-9873-8

DO - 10.1007/s11207-011-9873-8

M3 - Article

VL - 276

SP - 19

EP - 33

JO - Solar Physics

JF - Solar Physics

SN - 0038-0938

IS - 1-2

ER -