We study the impact of the QCD DGLAP evolution on the geometric scaling of the gluon distributions which is expected to hold at small x within the saturation models. To this aim we solve the DGLAP evolution equations with the initial conditions provided along the critical line Q²=Q_s²(x) with Q_s²(x) \sim x^-\lambda and satisfying geometric scaling. Both fixed and running coupling cases are studied. We show that in the fixed coupling case the geometric scaling at low x is stable against the DGLAP evolution for sufficiently large values of the parameter \lambda and in the double logarithmic approximation of the DGLAP evolution this happens for \lambda \ge 4N_c\alpha _s/\pi . In the running coupling case geometric scaling is found to be approximately preserved at very small x. The residual geometric scaling violation in this case can be approximately factored out and the corresponding form-factor controlling this violation is found.

PACS numbers: 12.38.Cy