We derive a two-dimensional (2D) extension of a recently developed formalism for slow-fast quasilinear (QL) systems subject to fast instabilities. The emergent dynamics of these systems is characterized by a slow evolution of (suitably defined) mean fields coupled to marginally stable, fast fluctuation fields. By exploiting this scale separation, an efficient hybrid fast-eigenvalue/slow-initial-value solution algorithm can be developed in which the amplitude of the fast fluctuations is slaved to the slowly evolving mean fields to ensure marginal stability-and temporal scale separation-is maintained. For 2D systems, the fluctuation eigenfunctions are labeled by their Fourier wave numbers characterizing spatial variability in that extended spatial direction, and the marginal mode(s) must coincide with the fastest-growing mode(s) over all admissible Fourier wave numbers. Here we derive an ordinary differential equation governing the slow evolution of the wave number of the fastest-growing fluctuation mode that simultaneously must be slaved to the mean dynamics to ensure the mode has zero growth rate. We illustrate the procedure in the context of a 2D model partial differential equation that shares certain attributes with the equations governing strongly stratified shear flows and other strongly constrained forms of geophysical turbulence in extreme parameter regimes. The slaved evolution follows one or more marginal stability manifolds, which constitute select state-space structures that are not invariant under the full flow dynamics yet capture quasicoherent structures in physical space in a manner analogous to invariant solutions identified in, e.g., transitionally turbulent shear flows. Accordingly, we propose that marginal stability manifolds are central organizing structures in a dynamical systems description of certain classes of multiscale flows in which scale separation justifies a QL approximation of the dynamics.