Abstract. An abstract framework for the theory of statistical solutions is developed for general evolution equations, extending the theory initially developed for the three-dimensional incompressible Navier-Stokes equations. The motivation for this concept is to model the evolution of uncertainties on the initial conditions for systems which have global solutions that are not known to be unique. Both concepts of statistical solution in trajectory space and in phase space are given, and the corresponding results of existence of statistical solution for the associated initial value problems are proved. The wide applicability of the theory is illustrated with the very incompressible Navier-Stokes equations, a reaction-diffusion equation, and a nonlinear wave equation, all displaying the property of global existence of weak solutions without a known result of global uniqueness.