The development of zonal flows on a midlatitude β-plane subject to a time-varying topographic forcing is investigated in a series of numerical integrations in which the forcing is concentrated at large scales, and in which the usual two-dimensional inverse energy cascade is absent. In contrast to the case of small-scale forcing, where mixing of potential vorticity occurs largely through the action of small-scale eddies, mixing of potential vorticity in this case occurs predominantly in latitudinally localized Rossby wave critical layer regions, whose width grows continuously in time due to the entrainment of background fluid. The potential vorticity is found to organize into a piecewise constant staircase-like profile, monotonic in latitude, provided the ratio L Rh /L f 1, where L Rh is the usual Rhines scale and L f is the scale of the forcing; this may be regarded as supplemental to the condition L Rh /L ε 6, where L ε = (ε/β 3 ) 1/5 and ε is the rate of energy input, obtained recently [R. K. Scott and D. G. Dritschel, "The structure of zonal jets in geostrophic turbulence," J. Fluid Mech. 711, 576-598 (2012)] for the case of small-scale forcing. The numerical results further suggest that the nature of the potential vorticity mixing is controlled by the ratio L ε /L f , and occurs predominantly in critical layers when L ε /L f 1/6. A combined condition for staircase formation may therefore be expressed asFinally, in a separate set of experiments it is shown that when forcing is represented by an additive source term in the evolution equation, as is common practice in numerical investigations of two-dimensional turbulence, the effect of non-conservation of potential vorticity may obscure the development of the staircase profile in the critical layer mixing dominated regime. C 2012 American Institute of Physics.[http://dx