  We develop a full microscopic theory for a nonlinear phase space filling
(NPSF) in strongly coupled two-dimensional polaritonic lattices. Ubiquitous in
polaritonic experiments, the theoretical description of NPSF, remains limited
to perturbative treatment and homogeneous samples. In this study, we go beyond
the existing theoretical description and discover the broad scope of regimes
where NPSF crucially modifies the optical response. Studying the quantum
effects of non-bosonicity, cooperative light-matter coupling, and Coulomb
blockade, we reveal several regimes for observing the nonlinear Rabi splitting
quench due to the phase space filling. Unlike prior studies, we derive
nonlinear Rabi frequency scaling all the way to the saturation limit and show
that the presence of a lattice potential leads to qualitatively distinct
nonlinearity. We concentrate on three regimes of NPSF: 1) planar; 2) fractured;
and 3) ultralocalized. In planar saturation, the Rabi frequency decreases
exponentially as a function of exciton density. For the fractured case, where
excitons form a lattice with sites exceeding the exciton size, we discover fast
NPSF at low occupation in the lattice. This is followed by slower NPSF as the
medium becomes fully saturated. This behavior is particularly pronounced in the
presence of Coulomb (or Rydberg) blockade, where regions of fast and slow NPSF
depend on the strength of repulsion. For the ultralocalized NPSF, we observe
the square-root saturation typical to the collection of two-level systems. Our
findings can help describing recent observations of strong nonlinearity in
heterobilayers of transition metal dichalcogenides where Moir{\'e} lattices
emerge naturally [Nature \textbf{591}, 61 (2021)]. The theory also opens the
prospects for engineering strongly nonlinear responses of polaritonic lattices
with patterned samples, driving polaritonics into the quantum regime.
