应用一阶与二阶 adjoint 方法至含旋转反馈的冰床等离子调整
地球物理
2024-11-13 v2
摘要
本文复现并扩展了 Crawford 等人(2018)关于冰床等离子调整(GIA)的 adjoint 理论。现在已包含旋转反馈,并首次描述二阶 adjoint 方法的应用。一阶 adjoint 方法提供了计算特定目标函数灵敏度核的有效方法,而二阶 adjoint 方法则以 Hessian 核形式提供二阶导数信息。这些后者核是高效 Newton 类型优化方案以及非线性反问题不确定性定量方法所必需的。最重要的是,整个理论已重新表述,以简化其他 GIA 社区成员的实现。特别是,Crawford 等人(2018)引入的 GIA 前向问题的速率形式公式已被取代为用于 modeling GIA 的常规方程式。一阶和二阶 adjoint 问题的实现应该在现有的 GIA 代码以及新的代码中相对容易,仅需包含更一般的力项。
引用
@article{arxiv.2408.13564,
title = {Application of first- and second-order adjoint methods to glacial isostatic adjustment incorporating rotational feedbacks},
author = {Ziheng Yu and David Al-Attar and Frank Syvret and Andrew J. Lloyd},
journal= {arXiv preprint arXiv:2408.13564},
year = {2024}
}