We carry out model studies of turbulence quantities for flow over two-dimensional hills using anon-hydrostatic version of the Regional Atmospheric Modeling System (RAMS). We test two turbulence closure models: the first one is an explicit Algebraic Reynolds Stress Model (ARSM) and the second one is a combination of the ARSM and a transport equation for the shear stress <(uw)over bar>. Model predictions for the turbulent stresses are compared with data from a wind-tunnel experiment containing isolated two-dimensional hills of varying slope. From the comparison, it is concluded that the first model can only predict the normal stresses adequately while the second model provides satisfactory predictions for the normal stresses as well as giving an improved result for the shear stress <(uw)over bar>.