Two approaches reveal nonlinear, regime-like behavior of the atmosphere. The probability density function (PDF) of projections onto planes defined by empirical orthogonal functions (EOFs) of the GCM's 500-hPa height field indicate significant departures from Gaussianity though no more than one maximum is found. Regime-like behavior is even more evident in the mean phase space tendencies where in some planes signatures of more than one equilibrium point are distinct (see figure).
The nonlinear structure of the mean tendencies suggests a piecewise linear behavior. It is examined to what degree a local linear decomposition describes this structure and thus justifies a multinormal approximation. The PDF of the geopotential height is decomposed into a sum of component Gaussians which will be related to the piecewise linear fit. Finally, the Fokker-Planck equation is used to relate the PDFs to the mean tendencies of the GCM.
Figure: Mean tendencies with damping removed projected onto the phase space plane spanned by the 1st and 4th EOFs.