Guardado en:
Detalles Bibliográficos
Autores principales: Yu, Ziheng, Al-Attar, David, Syvret, Frank, Lloyd, Andrew J.
Formato: Preprint
Publicado: 2024
Materias:
Acceso en línea:https://arxiv.org/abs/2408.13564
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866913576256536576
author Yu, Ziheng
Al-Attar, David
Syvret, Frank
Lloyd, Andrew J.
author_facet Yu, Ziheng
Al-Attar, David
Syvret, Frank
Lloyd, Andrew J.
contents This paper revisits and extends the adjoint theory for glacial isostatic adjustment (GIA) of Crawford et al. (2018). Rotational feedbacks are now incorporated, and the application of the second-order adjoint method is described for the first time. The first-order adjoint method provides an efficient means for computing sensitivity kernels for a chosen objective functional, while the second-order adjoint method provides second-derivative information in the form of Hessian kernels. These latter kernels are required by efficient Newton-type optimisation schemes and within methods for quantifying uncertainty for non-linear inverse problems. Most importantly, the entire theory has been reformulated so as to simplify its implementation by others within the GIA community. In particular, the rate-formulation for the GIA forward problem introduced by Crawford et al. (2018) has been replaced with the conventional equations for modelling GIA in laterally heterogeneous earth models. The implementation of the first- and second-order adjoint problems should be relatively easy within both existing and new GIA codes, with only the inclusions of more general force terms being required.
format Preprint
id arxiv_https___arxiv_org_abs_2408_13564
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Application of first- and second-order adjoint methods to glacial isostatic adjustment incorporating rotational feedbacks
Yu, Ziheng
Al-Attar, David
Syvret, Frank
Lloyd, Andrew J.
Geophysics
This paper revisits and extends the adjoint theory for glacial isostatic adjustment (GIA) of Crawford et al. (2018). Rotational feedbacks are now incorporated, and the application of the second-order adjoint method is described for the first time. The first-order adjoint method provides an efficient means for computing sensitivity kernels for a chosen objective functional, while the second-order adjoint method provides second-derivative information in the form of Hessian kernels. These latter kernels are required by efficient Newton-type optimisation schemes and within methods for quantifying uncertainty for non-linear inverse problems. Most importantly, the entire theory has been reformulated so as to simplify its implementation by others within the GIA community. In particular, the rate-formulation for the GIA forward problem introduced by Crawford et al. (2018) has been replaced with the conventional equations for modelling GIA in laterally heterogeneous earth models. The implementation of the first- and second-order adjoint problems should be relatively easy within both existing and new GIA codes, with only the inclusions of more general force terms being required.
title Application of first- and second-order adjoint methods to glacial isostatic adjustment incorporating rotational feedbacks
topic Geophysics
url https://arxiv.org/abs/2408.13564