TY - JOUR
T1 - Drag parameter estimation using gradients and hessian from a polynomial chaos model surrogate
AU - Sraj, Ihab
AU - Iskandarani, Mohamed
AU - Carlisle Thacker, W.
AU - Srinivasan, Ashwanth
AU - Knio, Omar M.
PY - 2014/2
Y1 - 2014/2
N2 - A variational inverse problem is solved using polynomial chaos expansions to infer several critical variables in the Hybrid Coordinate Ocean Model's (HYCOM's) wind drag parameterization. This alternative to the Bayesian inference approach in Sraj et al. avoids the complications of constructing the full posterior with Markov chain Monte Carlo sampling. It focuses instead on identifying the center and spread of the posterior distribution. The present approach leverages the polynomial chaos series to estimate, at very little extra cost, the gradients and Hessian of the cost function during minimization. The Hessian's inverse yields an estimate of the uncertainty in the solution when the latter's probability density is approximately Gaussian. The main computational burden is an ensemble of realizations to build the polynomial chaos expansion; no adjoint code or additional forward model runs are needed once the series is available. The ensuing optimal parameters are compared to those obtained in Sraj et al. where the full posterior distribution was constructed. The similarities and differences between the new methodology and a traditional adjoint-based calculation are discussed
AB - A variational inverse problem is solved using polynomial chaos expansions to infer several critical variables in the Hybrid Coordinate Ocean Model's (HYCOM's) wind drag parameterization. This alternative to the Bayesian inference approach in Sraj et al. avoids the complications of constructing the full posterior with Markov chain Monte Carlo sampling. It focuses instead on identifying the center and spread of the posterior distribution. The present approach leverages the polynomial chaos series to estimate, at very little extra cost, the gradients and Hessian of the cost function during minimization. The Hessian's inverse yields an estimate of the uncertainty in the solution when the latter's probability density is approximately Gaussian. The main computational burden is an ensemble of realizations to build the polynomial chaos expansion; no adjoint code or additional forward model runs are needed once the series is available. The ensuing optimal parameters are compared to those obtained in Sraj et al. where the full posterior distribution was constructed. The similarities and differences between the new methodology and a traditional adjoint-based calculation are discussed
UR - http://www.scopus.com/inward/record.url?scp=84893854293&partnerID=8YFLogxK
U2 - 10.1175/MWR-D-13-00087.1
DO - 10.1175/MWR-D-13-00087.1
M3 - Article
AN - SCOPUS:84893854293
SN - 0027-0644
VL - 142
SP - 933
EP - 941
JO - Monthly Weather Review
JF - Monthly Weather Review
IS - 2
ER -