  The differential cross section of coherent photo-production of heavy
quarkonia on nuclear targets is calculated within the QCD color dipole
formalism. The higher-twist nuclear shadowing corresponding to the $\bar QQ$
Fock component of the photon, is calculated including the correlation between
dipole orientation $\vec r$ and impact parameter of a collision $\vec b$, which
is related to the transverse momentum transfer via Fourier transform. We also
included the leading twist gluon shadowing corresponding to higher Fock
components of the photon containing gluons, which have specifically short
coherence time, especially for multi-gluon components, even at very high
energies. The contribution of such fluctuating gluonic dipole is calculated
employing the path-integral technique. Our results are in good agreement with
recent ALICE data on charmonium production in ultra-peripheral nuclear
collisions.
