### Abstract

A major goal of helioseismology is the three-dimensional reconstruction of the three velocity components of convective flows in the solar interior from sets of wave travel-time measurements. For small amplitude flows, the forward problem is described in good approximation by a large system of convolution equations. The input observations are highly noisy random vectors with a known dense covariance matrix. This leads to a large statistical linear inverse problem. Whereas for deterministic linear inverse problems several computationally efficient minimax optimal regularization methods exist, only one minimax-optimal linear estimator exists for statistical linear inverse problems: the Pinsker estimator. However, it is often computationally inefficient because it requires a singular value decomposition of the forward operator or it is not applicable because of an unknown noise covariance matrix, so it is rarely used for real-world problems. These limitations do not apply in helioseismology. We present a simplified proof of the optimality properties of the Pinsker estimator and show that it yields significantly better reconstructions than traditional inversion methods used in helioseismology, i.e. regularized least squares (Tikhonov regularization) and SOLA (approximate inverse) methods. Moreover, we discuss the incorporation of the mass conservation constraint in the Pinsker scheme using staggered grids. With this improvement we can reconstruct not only horizontal, but also vertical velocity components that are much smaller in amplitude.

Original language | English (US) |
---|---|

Article number | 105002 |

Journal | Inverse Problems |

Volume | 32 |

Issue number | 10 |

DOIs | |

State | Published - Aug 5 2016 |

### Fingerprint

### Keywords

- helioseismology
- pinsker estimator
- statistical inverse problem

### ASJC Scopus subject areas

- Theoretical Computer Science
- Signal Processing
- Mathematical Physics
- Computer Science Applications
- Applied Mathematics

### Cite this

*Inverse Problems*,

*32*(10), [105002]. https://doi.org/10.1088/0266-5611/32/10/105002

**Pinsker estimators for local helioseismology : Inversion of travel times for mass-conserving flows.** / Fournier, Damien; Gizon, Laurent; Holzke, Martin; Hohage, Thorsten.

Research output: Contribution to journal › Article

*Inverse Problems*, vol. 32, no. 10, 105002. https://doi.org/10.1088/0266-5611/32/10/105002

}

TY - JOUR

T1 - Pinsker estimators for local helioseismology

T2 - Inversion of travel times for mass-conserving flows

AU - Fournier, Damien

AU - Gizon, Laurent

AU - Holzke, Martin

AU - Hohage, Thorsten

PY - 2016/8/5

Y1 - 2016/8/5

N2 - A major goal of helioseismology is the three-dimensional reconstruction of the three velocity components of convective flows in the solar interior from sets of wave travel-time measurements. For small amplitude flows, the forward problem is described in good approximation by a large system of convolution equations. The input observations are highly noisy random vectors with a known dense covariance matrix. This leads to a large statistical linear inverse problem. Whereas for deterministic linear inverse problems several computationally efficient minimax optimal regularization methods exist, only one minimax-optimal linear estimator exists for statistical linear inverse problems: the Pinsker estimator. However, it is often computationally inefficient because it requires a singular value decomposition of the forward operator or it is not applicable because of an unknown noise covariance matrix, so it is rarely used for real-world problems. These limitations do not apply in helioseismology. We present a simplified proof of the optimality properties of the Pinsker estimator and show that it yields significantly better reconstructions than traditional inversion methods used in helioseismology, i.e. regularized least squares (Tikhonov regularization) and SOLA (approximate inverse) methods. Moreover, we discuss the incorporation of the mass conservation constraint in the Pinsker scheme using staggered grids. With this improvement we can reconstruct not only horizontal, but also vertical velocity components that are much smaller in amplitude.

AB - A major goal of helioseismology is the three-dimensional reconstruction of the three velocity components of convective flows in the solar interior from sets of wave travel-time measurements. For small amplitude flows, the forward problem is described in good approximation by a large system of convolution equations. The input observations are highly noisy random vectors with a known dense covariance matrix. This leads to a large statistical linear inverse problem. Whereas for deterministic linear inverse problems several computationally efficient minimax optimal regularization methods exist, only one minimax-optimal linear estimator exists for statistical linear inverse problems: the Pinsker estimator. However, it is often computationally inefficient because it requires a singular value decomposition of the forward operator or it is not applicable because of an unknown noise covariance matrix, so it is rarely used for real-world problems. These limitations do not apply in helioseismology. We present a simplified proof of the optimality properties of the Pinsker estimator and show that it yields significantly better reconstructions than traditional inversion methods used in helioseismology, i.e. regularized least squares (Tikhonov regularization) and SOLA (approximate inverse) methods. Moreover, we discuss the incorporation of the mass conservation constraint in the Pinsker scheme using staggered grids. With this improvement we can reconstruct not only horizontal, but also vertical velocity components that are much smaller in amplitude.

KW - helioseismology

KW - pinsker estimator

KW - statistical inverse problem

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

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

U2 - 10.1088/0266-5611/32/10/105002

DO - 10.1088/0266-5611/32/10/105002

M3 - Article

AN - SCOPUS:85009228768

VL - 32

JO - Inverse Problems

JF - Inverse Problems

SN - 0266-5611

IS - 10

M1 - 105002

ER -