This paper shows that the most predictable components of internal variability in coupled atmosphere-ocean models are remarkably similar to the most predictable components of climate models without interactive ocean dynamics (i.e., models whose ocean is represented by a 50-m-deep slab ocean mixed layer with no interactive currents). Furthermore, a linear regression model derived solely from dynamical model output can skillfully predict observed anomalies in these components at least a year or two in advance, indicating that these model-derived components and associated linear dynamics are realistic. These results suggest that interactive ocean circulation is not essential for the existence of multiyear predictability previously identified in coupled models and observations. decadal prediction | decadal predictability | average predictability time | CMIP I n recent years, climate prediction on decadal time scales has gained increased attention because of its growing scientific, geopolitical, and societal importance (1-5). Sources of decadal predictability often are divided into two kinds: predictability caused by external forcing, such as changes in solar insolation, volcanic aerosols, and anthropogenic greenhouse gases, and predictability due to internal variability arising naturally from the coupled atmosphere-ocean-land-ice climate system (6). Climate models suggest that certain structures of internal variability in sea surface temperatures (SST), such as the Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO), are predictable on decadal/multidecadal time scales (5, 7-9). The precise mechanisms for this decadal predictability are not clear, although the dynamics of ocean circulation are widely believed to play a major role (10-17).Recently, a few studies have challenged the notion that interactive ocean dynamics play a dominant role in decadal predictability (18,19). These studies are based on integrations of atmospheric global circulation models coupled to a slab ocean mixed layer, in which ocean circulation is a prescribed function of time and the ocean interacts with the atmosphere only thermodynamically, through radiative, sensible, and latent heat fluxes. Despite the absence of interactive ocean dynamics, these models can produce realistic variability that is predictable on interannual and longer time scales. For example, El Niño-Southern Oscillation (ENSO)-like variability can arise on interannual and decadal time scales from such models (20)(21)(22). Also, the main features of the observed AMO (e.g., spatial pattern, power spectra, and associated atmospheric circulation) have been reproduced in models without interactive ocean dynamics (23). However, not all details of these simulations are perfect: In some locations, the lag correlation between AMO and surface heat flux has the opposite sign relative to coupled models and observations (14,15). On the other hand, these heat fluxes tend to be nearly canceled by the ocean heat transport convergence on long time scales (24), raisin...