In this thesis, parameterizations of non-linear interactions in quasigeostrophic (QG) flows for severely truncated models (STM) and Large Eddy Simulations (LES) are studied. Firstly, using Direct Numerical Simulations (DNS), atmospheric barotropic flows over topography are examined, and it is established that such flows exhibit multiple equilibrium states for a wide range of parameters. A STM is then constructed, consisting of the large scale zonal flow and a topographic mode. It is shown that, qualitatively, this system behaves similarly to the DNS as far as the interaction between the zonal flow and topography is concerned, and, in particular, exhibits multiple equilibrium states. By fitting the analytical form of the topographic stationary wave amplitude, obtained from the STM, to the results obtained from DNS, renormalized dissipation and rotation parameters are obtained. The usage of renormalized parameters in the STM results in better quantitative agreement with the DNS.¶ In the second type of problem, subgrid-scale parameterizations in LES are investigated with both atmospheric and oceanic parameters. ...