Abstract A hybrid coupled ocean–atmosphere model is used to investigate low-frequency variability in the climate system. The model's atmospheric component is a Budyko-Sellers-North, two-dimensional energy-balance model; the oceanic component is a simplified general circulation model. The coupled model is confined to an idealized, rectangular North Atlantic basin. In the present model version, the ocean density depends exclusively on temperature. An interdecadal oscillation with a period of 40–50 years is found in the hybrid coupled model when model parameters are within the climatological range, even though density does not depend on salinity. This interdecadal oscillation is characterized by a pair of vortices of opposite signs, that grow and decay in quadrature with each other in the ocean's upper layer; their centers follow each other anticlockwise through the northwestern quadrant of the model domain. The interdecadal oscillation's physical mechanism resembles that of the interdecadal oscillation anal...