We identify three major areas of ignorance which limit predictability in current ocean GCMs. One is the very crude representation of subgrid-scale mixing processes. These processes are parameterized with coefficients whose values and variations in space and time are poorly known. A second problem derives from the fact that ocean models generally contain multiple equilibria and bifurcations, but there is no agreement as to where the current ocean sits with respect to the bifurcations. A third problem arises from the fact that ocean circulations are highly nonlinear, but only weakly dissipative, and therefore are potentially chaotic. The few studies that have looked at this kind of behavior have not answered fundamental questions, such as what are the major sources of error growth in model projections, and how large is the chaotic behavior relative to realistic changes in climate forcings. Advances in computers will help alleviate some of these problems, for example by making it more practical to explore to what extent the evolution of the oceans is chaotic. However models will have to rely on parameterizations of key small-scale processes such as diapycnal mixing for a long time. To make more immediate progress here requires the development of physically based prognostic parameterizations and coupling the mixing to its energy sources. Another possibly fruitful area of investigation is the use of paleoclimate data on changes in the ocean circulation to constrain more tightly the stability characteristics of the ocean circulation. (© 2004 International Union of Geodesy and Geophysics and the American Geophysical Union)