The purpose of this paper is to investigate the close relation between Okounkov bodies and Zariski decompositions of pseudo effective divisors on smooth projective surfaces. Firstly, we completely determine the limiting Okounkov bodies on such surfaces, and give applications to Nakayama constants and Seshadri constants. Secondly, we study how the shapes of Okounkov bodies change as we vary the divisors in the big cone.