Application of linear dynamical mode decomposition to surface air temperature in 20th century

According to recent study [1], the current state-of-the-art climate models lack the substantial part of internal multidecadal climate signal which is observed in the 20th century surface air temperature reanalysis data as a global stadium wave (GSW). In the presented work we further investigate this phenomenon using the recently developed method [2] of empirical spatio-temporal data decomposition into linear dynamical modes (LDMs). The important property of LDMs is their ability to take into account the time scales of the system evolution (they are extracted from observed dataset by the Bayesian optimization technique) better than some other linear techniques, e.g. traditional empirical orthogonal function decomposition. Like any linear decomposition, it provides the time series of principal components and corresponding spatial patterns.
We modify the initially developed LDM decomposition to make it possible to take into account a prescribed external forcing (like CO2 emissions, sun activity etc.) and then find part of variability which may be considered as an internal climate dynamics decomposed into set of modes with different time scales, and hence may be helpful in GSW interpretation. The results of applying the method to the 20th century surface air temperature with different ways of forcing inclusion will be presented and discussed.

1. Kravtsov, S., Grimm, C., & Gu, S. (2018). Global-scale multidecadal variability missing in state-of-the-art climate models. Npj Climate and Atmospheric Science, 1(1), 34. https://doi.org/10.1038/s41612-018-0044-6
2. Gavrilov, A., Seleznev, A., Mukhin, D., Loskutov, E., Feigin, A., & Kurths, J. (2018). Linear dynamical modes as new variables for data-driven ENSO forecast. Climate Dynamics. https://doi.org/10.1007/s00382-018-4255-7