Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet
The majority of research on efficient and scalable algorithms in computational science and engineering has focused on the forward problem: given parameter inputs, solve the governing equations to determine output quantities of interest. In contrast, here we consider the broader question: given a (la...
Published in: | Journal of Computational Physics |
---|---|
Main Authors: | , , , |
Format: | Text |
Language: | unknown |
Published: |
2014
|
Subjects: | |
Online Access: | http://arxiv.org/abs/1410.1221 https://doi.org/10.1016/j.jcp.2015.04.047 |
id |
ftarxivpreprints:oai:arXiv.org:1410.1221 |
---|---|
record_format |
openpolar |
spelling |
ftarxivpreprints:oai:arXiv.org:1410.1221 2023-09-05T13:12:52+02:00 Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet Isaac, Tobin Petra, Noemi Stadler, Georg Ghattas, Omar 2014-10-05 http://arxiv.org/abs/1410.1221 https://doi.org/10.1016/j.jcp.2015.04.047 unknown http://arxiv.org/abs/1410.1221 doi:10.1016/j.jcp.2015.04.047 Mathematics - Optimization and Control Mathematics - Numerical Analysis Statistics - Computation Statistics - Methodology 35Q62 62F15 35R30 35Q93 65C60 49M15 86A40 text 2014 ftarxivpreprints https://doi.org/10.1016/j.jcp.2015.04.047 2023-08-16T13:27:21Z The majority of research on efficient and scalable algorithms in computational science and engineering has focused on the forward problem: given parameter inputs, solve the governing equations to determine output quantities of interest. In contrast, here we consider the broader question: given a (large-scale) model containing uncertain parameters, (possibly) noisy observational data, and a prediction quantity of interest, how do we construct efficient and scalable algorithms to (1) infer the model parameters from the data (the deterministic inverse problem), (2) quantify the uncertainty in the inferred parameters (the Bayesian inference problem), and (3) propagate the resulting uncertain parameters through the model to issue predictions with quantified uncertainties (the forward uncertainty propagation problem)? We present efficient and scalable algorithms for this end-to-end, data-to-prediction process under the Gaussian approximation and in the context of modeling the flow of the Antarctic ice sheet and its effect on sea level. The ice is modeled as a viscous, incompressible, creeping, shear-thinning fluid. The observational data come from InSAR satellite measurements of surface ice flow velocity, and the uncertain parameter field to be inferred is the basal sliding parameter. The prediction quantity of interest is the present-day ice mass flux from the Antarctic continent to the ocean. We show that the work required for executing this data-to-prediction process is independent of the state dimension, parameter dimension, data dimension, and number of processor cores. The key to achieving this dimension independence is to exploit the fact that the observational data typically provide only sparse information on model parameters. This property can be exploited to construct a low rank approximation of the linearized parameter-to-observable map. Text Antarc* Antarctic Ice Sheet ArXiv.org (Cornell University Library) Antarctic The Antarctic Journal of Computational Physics 296 348 368 |
institution |
Open Polar |
collection |
ArXiv.org (Cornell University Library) |
op_collection_id |
ftarxivpreprints |
language |
unknown |
topic |
Mathematics - Optimization and Control Mathematics - Numerical Analysis Statistics - Computation Statistics - Methodology 35Q62 62F15 35R30 35Q93 65C60 49M15 86A40 |
spellingShingle |
Mathematics - Optimization and Control Mathematics - Numerical Analysis Statistics - Computation Statistics - Methodology 35Q62 62F15 35R30 35Q93 65C60 49M15 86A40 Isaac, Tobin Petra, Noemi Stadler, Georg Ghattas, Omar Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
topic_facet |
Mathematics - Optimization and Control Mathematics - Numerical Analysis Statistics - Computation Statistics - Methodology 35Q62 62F15 35R30 35Q93 65C60 49M15 86A40 |
description |
The majority of research on efficient and scalable algorithms in computational science and engineering has focused on the forward problem: given parameter inputs, solve the governing equations to determine output quantities of interest. In contrast, here we consider the broader question: given a (large-scale) model containing uncertain parameters, (possibly) noisy observational data, and a prediction quantity of interest, how do we construct efficient and scalable algorithms to (1) infer the model parameters from the data (the deterministic inverse problem), (2) quantify the uncertainty in the inferred parameters (the Bayesian inference problem), and (3) propagate the resulting uncertain parameters through the model to issue predictions with quantified uncertainties (the forward uncertainty propagation problem)? We present efficient and scalable algorithms for this end-to-end, data-to-prediction process under the Gaussian approximation and in the context of modeling the flow of the Antarctic ice sheet and its effect on sea level. The ice is modeled as a viscous, incompressible, creeping, shear-thinning fluid. The observational data come from InSAR satellite measurements of surface ice flow velocity, and the uncertain parameter field to be inferred is the basal sliding parameter. The prediction quantity of interest is the present-day ice mass flux from the Antarctic continent to the ocean. We show that the work required for executing this data-to-prediction process is independent of the state dimension, parameter dimension, data dimension, and number of processor cores. The key to achieving this dimension independence is to exploit the fact that the observational data typically provide only sparse information on model parameters. This property can be exploited to construct a low rank approximation of the linearized parameter-to-observable map. |
format |
Text |
author |
Isaac, Tobin Petra, Noemi Stadler, Georg Ghattas, Omar |
author_facet |
Isaac, Tobin Petra, Noemi Stadler, Georg Ghattas, Omar |
author_sort |
Isaac, Tobin |
title |
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
title_short |
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
title_full |
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
title_fullStr |
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
title_full_unstemmed |
Scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the Antarctic ice sheet |
title_sort |
scalable and efficient algorithms for the propagation of uncertainty from data through inference to prediction for large-scale problems, with application to flow of the antarctic ice sheet |
publishDate |
2014 |
url |
http://arxiv.org/abs/1410.1221 https://doi.org/10.1016/j.jcp.2015.04.047 |
geographic |
Antarctic The Antarctic |
geographic_facet |
Antarctic The Antarctic |
genre |
Antarc* Antarctic Ice Sheet |
genre_facet |
Antarc* Antarctic Ice Sheet |
op_relation |
http://arxiv.org/abs/1410.1221 doi:10.1016/j.jcp.2015.04.047 |
op_doi |
https://doi.org/10.1016/j.jcp.2015.04.047 |
container_title |
Journal of Computational Physics |
container_volume |
296 |
container_start_page |
348 |
op_container_end_page |
368 |
_version_ |
1776202108969680896 |