  The behavior at bifurcation from global synchronization to partial
synchronization in finite networks of coupled oscillators is a complex
phenomenon, involving the intricate dynamics of one or more oscillators with
the remaining synchronized oscillators. This is not captured well by standard
macroscopic model reduction techniques which capture only the collective
behavior of synchronized oscillators in the thermodynamic limit. We introduce
two mesoscopic model reductions for finite sparse networks of coupled
oscillators to quantitatively capture the dynamics close to bifurcation from
global to partial synchronization. Our model reduction builds upon the method
of collective coordinates. We first show that standard collective coordinate
reduction has difficulties capturing this bifurcation. We identify a particular
topological structure at bifurcation consisting of a main synchronized cluster,
the oscillator that desynchronizes at bifurcation, and an intermediary node
connecting them. Utilizing this structure and ensemble averages we derive an
analytic expression for the mismatch between the true bifurcation from global
to partial synchronization and its estimate calculated via the collective
coordinate approach. This allows to calibrate the standard collective
coordinate approach without prior knowledge of which node will desynchronize.
We introduce a second mesoscopic reduction, utilizing the same particular
topological structure, which allows for a quantitative dynamical description of
the phases near bifurcation. The mesoscopic reductions significantly reduce the
computational complexity of the collective coordinate approach, reducing from
$\mathcal{O}(N^2)$ to $\mathcal{O}(1)$. We perform numerical simulations for
Erd\H{o}s-R\'enyi networks and for modified Barab\'asi-Albert networks
demonstrating excellent quantitative agreement at and close to bifurcation.
