We formulate and study a low-order nonlinear coupled ocean-atmosphere model with an emphasis on the impact of radiative and heat fluxes and of the frictional coupling between the two components. This model version extends a previous 24-variable version by adding a dynamical equation for the passive advection of temperature in the ocean, together with an energy balance model.The bifurcation analysis and the numerical integration of the model reveal the presence of low-frequency variability (LFV) concentrated on and near a longperiodic, attracting orbit. This orbit combines atmospheric and oceanic modes, and it arises for large values of the meridional gradient of radiative input and of frictional coupling. Chaotic behavior develops around this orbit as it loses its stability; this behavior is still dominated by the LFV on decadal and multidecadal time scales that is typical of oceanic processes. Atmospheric diagnostics also reveals the presence of predominant low-and high-pressure zones, as well as of a subtropical jet; these features recall realistic climatological properties of the oceanic atmosphere.Finally, a predictability analysis is performed. Once the decadal-scale periodic orbits develop, the coupled system's short-term instabilities -as measured by its Lyapunov exponents -are drastically reduced, indicating the ocean's stabilizing role on the atmospheric dynamics. On decadal time scales, the recurrence of the solution in a certain region of the invariant subspace associated with slow modes displays some extended predictability, as reflected by the oscillatory behavior of the error for the atmospheric variables at long lead times.

CapsuleState-of-the-Art statistical postprocessing techniques for ensemble forecasts are reviewed, together with the challenges posed by a demand for timely, high-resolution and reliable probabilistic information. Possible research avenues are also discussed.

Abstract. This paper describes a reduced-order quasi-geostrophic coupled ocean–atmosphere model that allows for an arbitrary number of atmospheric and oceanic modes to be retained in the spectral decomposition. The modularity of this new model allows one to easily modify the model physics. Using this new model, coined the "Modular Arbitrary-Order Ocean-Atmosphere Model" (MAOOAM), we analyse the dependence of the model dynamics on the truncation level of the spectral expansion, and unveil spurious behaviour that may exist at low resolution by a comparison with the higher-resolution configurations. In particular, we assess the robustness of the coupled low-frequency variability when the number of modes is increased. An "optimal" configuration is proposed for which the ocean resolution is sufficiently high, while the total number of modes is small enough to allow for a tractable and extensive analysis of the dynamics.

A stochastic subgrid-scale parameterization based on the Ruelle's response theory and proposed in Wouters and Lucarini [2012] is tested in the context of a low-order coupled ocean-atmosphere model for which a part of the atmospheric modes are considered as unresolved. A natural separation of the phase-space into an invariant set and its complement allows for an analytical derivation of the different terms involved in the parameterization, namely the average, the fluctuation and the long memory terms. In this case, the fluctuation term is an additive stochastic noise. Its application to the low-order system reveals that a considerable correction of the low-frequency variability along the invariant subset can be obtained, provided that the coupling is sufficiently weak. This new approach of scale separation opens new avenues of subgrid-scale parameterizations in multiscale systems used for climate forecasts.