Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
https://doi.org/10.5194/gmd-2018-48
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.
Development and technical paper
06 Mar 2018
Review status
This discussion paper is a preprint. It is a manuscript under review for the journal Geoscientific Model Development (GMD).
Parameter space Kalman smoothers for multi-decadal climate analysis in high resolution coupled Global Circulation Models
Javier García-Pintado and André Paul MARUM – Center for Marine environmental Sciences and Department of Geosciences, University of Bremen, Bremen, Germany
Abstract. In climate reanalyses for multi-decadal or longer scales with coupled atmosphere-ocean General Circulation models (CGCMs) it can be assumed that the growth of prediction errors arises chiefly from imprecisely known model parameters, which have a nonlinear relationship with the climate observations (paleoclimate proxies). Also, high-resolution CGCMs for climate analysis are extremely expensive to run, which constrains the applicability of assimilation schemes. In a model framework where we assume that model dynamic parameters account for (nearly) all forecast errors at observation times, we compare two computationally efficient iterative schemes for approximate nonlinear model parameter estimation and joint flux estimation (taking the specific shape of freshwater from melting in the Greenland ice sheet), and its physically consistent state. First, a trivial adaptation of the strong constraint incremental 4D-Var formulation leads to what we refer to as the parameter space iterative extended Kalman smoother (pIKS); a Gauss-Newton scheme. Second, a so-called parameter space fractional Kalman smoother (pFKS) is an alternative controlled-step line search, which can potentially be a more stable approach. While these iterative schemes have been used in data assimilation, we revisit them together within the context of parameter estimation in climate reanalysis, as compared to the more general 4D-Var formulation. Then, the two schemes are evaluated in numerical experiments with a simple 1D energy balance model (Ebm1D) and with a fully-coupled Community Earth System Model (CESM v1.2). Firstly, with Ebm1D the pFKS obtains a cost function similar to the adjoint method with highly reduced computational cost, while an ensemble transform Kalman filter with an m = 60 ensemble size (ETKF60) behaves slightly worse. The pIKS behaves worse than the ETKF60, but an ETKF10 (m = 10) is even worst. Accordingly, with CESM we evaluate the pKFS and the ETKF60 along with an ETKF with Gaussian Anamorphosis (ETKF-GA60). From all the options, the pFKS has the lowest cost function and seems the favored overall option under heavy computational restrictions, but the ETKF obtains better estimates of the flux term.
Citation: García-Pintado, J. and Paul, A.: Parameter space Kalman smoothers for multi-decadal climate analysis in high resolution coupled Global Circulation Models, Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2018-48, in review, 2018.
Javier García-Pintado and André Paul
Javier García-Pintado and André Paul
Javier García-Pintado and André Paul

Viewed

Total article views: 263 (including HTML, PDF, and XML)

HTML PDF XML Total BibTeX EndNote
209 48 6 263 10 8

Views and downloads (calculated since 06 Mar 2018)

Cumulative views and downloads (calculated since 06 Mar 2018)

Viewed (geographical distribution)

Total article views: 263 (including HTML, PDF, and XML)

Thereof 262 with geography defined and 1 with unknown origin.

Country # Views %
  • 1

Saved

Discussed

Latest update: 19 Jun 2018
Publications Copernicus
Download
Short summary
Earth system models integrate the interactions of atmosphere, ocean, land, ice, and biosphere to estimate the state of regional and global climate under a wide variety of conditions. Here, we focus on methods for online estimation of uncertain model parameters and corresponding climate reconstruction based on available paleo-climate proxies. A supercomputing experiment for the Last Glacial Maximum, shows positive results and identifies areas in the global ocean in need of further exploration.
Earth system models integrate the interactions of atmosphere, ocean, land, ice, and biosphere to...
Share