The spectrum of water is ubiquitous and therefore the demand for a reliable theoretical model for it remains pressing. The treatment of nuclear motion using methods based on the variational principle has led to a major advance in the development of the model but places great emphasis on the accurate treatment of the underlying potential energy surface. Progress in constructing high-accuracy potentials both ab initio and by fitting to spectroscopic data is discussed.