Abstract

We investigate the Oneshot Optimization strategy introduced by Hamdi and Griewank for the applicability and efficiency to identify parameters in models of the earth's climate system. Parameters of a box model of the North Atlantic Thermohaline Circulation are optimized with respect to the fit of model output to data given by another model of intermediate complexity. Since the model is run into a steady state by a pseudo time-stepping, efficient techniques are necessary to avoid extensive recomputations or storing when using gradient-based local optimization algorithms. The Oneshot approach simultaneously updates state, adjoint and parameter values. For the required partial derivatives, the algorithmic/automatic differentiation tool TAF was used. Numerical results are compared to results obtained by the BFGS-quasi-Newton method.

Highlights

  • Parameter optimization is an important task in all kinds of climate models or models that simulate parts of the climate system, as for example ocean or atmospheric models

  • As mentioned in the introduction, in climate modeling an optimization is at first performed for steady states, which means in this example for temperatures and salinities which do not change in time anymore

  • We compare results to values obtained by the Limited-memory BFGS (L-BFGS) algorithm implemented by Zhu, Byrd, Nocedal and Morales, see [18], version 3.0 from 2011, without and with box constraints on the control parameters (L-BFGS-B) because we find that computed optimal parameter values of the BFGS and limited-memory BFGS (L-BFGS) method are far away from actual real world values

Read more

Summary

Introduction

Parameter optimization is an important task in all kinds of climate models or models that simulate parts of the climate system, as for example ocean or atmospheric models. (and this is the second point where engineering and geophysical flow problems are similar), the computation of steady states is often performed by running a transient model into the steady state This strategy is called pseudo time-stepping, since the time variable may be regarded as a kind of iteration counter. In the One-shot approach used here, the idea is that for fixed parameters there is a given (not necessarily (pseudo-) time-stepping) strategy to solve the state equations. We apply the One-shot approach to a box model of the North Atlantic This problem is different from the application in [6] in that the parameters enter in a nonlinear fashion resulting in so-called non-separable adjoints where the adjoint is no longer only the sum of a term on the state and a term on design.

One-shot Optimization Strategy
Problem formulation
One-shot iteration and its properties
Preconditioner B and the doubly augmented Lagrangian
BGFS update to avoid computation of full Jacobians and 2nd order derivatives
Application in Earth System Modeling
The optimization problem
Numerical results and discussion
Effect of the weighting factor α on the numerical results
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call