Abstract:Second order derivatives of model outputs with respect to input parameters are key to several applications in ice sheet modelling. For example, the ability to compute Hessian-vector products broadens the list of available optimisation methods, and facilitates certain kinds of parametric uncertainty quantification. Some modern ice sheet models are built on frameworks supporting algorithmic differentiation (AD), allowing for the computation of higher order derivatives with relative ease. However, many of our most widely-used models are not. A natural alternative might be to follow common practise in first order gradient computation and construct an approximate second-order adjoint model at the PDE level, which neglects the nonlinear dependence of ice viscosity on velocity. Here, we present such a model for the shallow-stream approximation allowing one to compute approximate second-order derivatives, and compare with full second-order derivates found using AD. We find that this produces Hessian-vector products that are superficially similar to those computed via AD. However, an analysis of the spectral decomposition of the Hessians calculated in each way reveals that the subspaces spanned by their eigenvectors diverge after the leading 4 modes, though divergence does not accelerate after this. We conclude that the utility of the approximate Hessian is case-dependent, and a full Hessian, likely computed using AD, should be used where high fidelity is required above very low rank.
From: Trystan Surawy-Stepney [view email]
[v1]
Mon, 20 Apr 2026 19:16:42 UTC (1,196 KB)
此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。