Linear Inverse Modeling for Coupled Atmosphere‐Ocean Ensemble Climate Prediction

Abstract

Abstract Paleoclimate data assimilation (PDA) experiments reconstruct climate fields by objectively blending information from climate models and proxy observations. Due to high computational cost and relatively low forecast skill, most reconstruction experiments omit the prediction step, where a forecast is made from the previously reconstructed state to the next time proxy data is available. In order to enable this critical aspect of PDA, we propose an efficient method of generating forecast ensembles of coupled climate fields using a linear inverse model (LIM). We describe the general calibration of a LIM on multiple fields using a two‐step empirical orthogonal function field compression to efficiently represent the state. We find that a LIM calibrated on global climate model (GCM) data yields skillful forecasts, including for out‐of‐sample tests on data from a different GCM. The deterministic forecast skill tests for scalar indices show that the LIM outperforms damped persistence at leads up to 3 years and has skill up to 10 years for global average sea surface temperature. Analysis of 1‐year forecasts reveals that the LIM captures dynamic climate features with local and remote predictability related to teleconnections. The forecast ensemble characteristics of the LIM, which in part determine the weighting of information for PDA experiments, show that the LIM generally produces ensemble forecast errors that are 10% to 70% larger than ensemble variance for 1‐year forecasts on data representative of the last millennium. These results show that the LIM produces ensembles with reasonable calibration but also that LIMs for PDA may require some variance tuning to work optimally for data assimilation experiments.

Keywords

• linear inverse model
• climate reanalysis
• coupled climate model
• ensemble forecasting
• empirical climate model