Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model

Hongyu Zhu, Noemi Petra, Georg Stadler, Tobin Isaac, Thomas J R Hughes, Omar Ghattas

Research output: Contribution to journalArticle

Abstract

We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two-and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. We show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

Original languageEnglish (US)
Pages (from-to)1477-1494
Number of pages18
JournalCryosphere
Volume10
Issue number4
DOIs
StatePublished - Jul 13 2016

Fingerprint

heat flux
ice sheet
inverse problem
ice flow
cost
advection-diffusion equation
flow velocity
energy
boundary condition
inversion
wavelength
prediction

ASJC Scopus subject areas

  • Water Science and Technology
  • Earth-Surface Processes

Cite this

Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model. / Zhu, Hongyu; Petra, Noemi; Stadler, Georg; Isaac, Tobin; Hughes, Thomas J R; Ghattas, Omar.

In: Cryosphere, Vol. 10, No. 4, 13.07.2016, p. 1477-1494.

Research output: Contribution to journalArticle

Zhu, Hongyu ; Petra, Noemi ; Stadler, Georg ; Isaac, Tobin ; Hughes, Thomas J R ; Ghattas, Omar. / Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model. In: Cryosphere. 2016 ; Vol. 10, No. 4. pp. 1477-1494.
@article{afe99004a4a44086b18b217c32deddcc,
title = "Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model",
abstract = "We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two-and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. We show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.",
author = "Hongyu Zhu and Noemi Petra and Georg Stadler and Tobin Isaac and Hughes, {Thomas J R} and Omar Ghattas",
year = "2016",
month = "7",
day = "13",
doi = "10.5194/tc-10-1477-2016",
language = "English (US)",
volume = "10",
pages = "1477--1494",
journal = "Cryosphere",
issn = "1994-0416",
publisher = "Copernicus Group",
number = "4",

}

TY - JOUR

T1 - Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model

AU - Zhu, Hongyu

AU - Petra, Noemi

AU - Stadler, Georg

AU - Isaac, Tobin

AU - Hughes, Thomas J R

AU - Ghattas, Omar

PY - 2016/7/13

Y1 - 2016/7/13

N2 - We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two-and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. We show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

AB - We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two-and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. We show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

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

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

U2 - 10.5194/tc-10-1477-2016

DO - 10.5194/tc-10-1477-2016

M3 - Article

VL - 10

SP - 1477

EP - 1494

JO - Cryosphere

JF - Cryosphere

SN - 1994-0416

IS - 4

ER -