The factorial moments of any Markov branching process describe the behaviour of its probability generating function $F(t,s)$ in the neighbourhood of the point $s=1$. They are applied to solve the forward Kolmogorov equation for the critical Markov branching process with geometric reproduction of particles. The solution includes quickly convergent recurrent iterations of polynomials. The obtained results on factorial moments enable computation of statistical measures as shape and skewness. They are also applicable to the comparison between critical geometric branching and linear birth-death processes.