In the presence of background rotation, conventional two-dimensional models of geostrophic flow in a rotating system usually include Ekman friction – associated with the no-slip condition at the bottom – by adding a linear term in the vorticity evolution equation. This term is proportional to E1/2 (where E is the Ekman number), and arises from the linear Ekman theory, which yields an expression for the vertical velocity produced by the thin Ekman layer at the flat bottom. In this paper, a two- dimensional model with Ekman damping is proposed using again the linear Ekman theory, but now including nonlinear Ekman terms in the vorticity equation. These terms represent nonlinear advection of relative vorticity as well as stretching effects. It ...