Accurately predicting how the cosmic abundance of neutral hydrogen evolves with redshift is a challenging problem facing modellers of galaxy formation . We investigate the predictions of four currently favoured semi-analytical galaxy formation models applied to the Millennium simulation for the mass function of cold neutral gas ( atomic and molecular ) in galaxies as a function of redshift , and we use these predictions to construct number counts for the next generation of all-sky neutral atomic hydrogen ( HI ) surveys . Despite the different implementations of the physical ingredients of galaxy formation , we find that the model predictions are broadly consistent with one another ; the key differences reflect how the models treat AGN feedback and how the timescale for star formation evolves with redshift . The models produce mass functions of cold gas in galaxies that are generally in good agreement with HI surveys at z =0 . Interestingly we find that these mass functions do not evolve significantly with redshift . Adopting a simple conversion factor for cold gas mass to HI mass that we apply to all galaxies at all redshifts , we derive mass functions of HI in galaxies from the predicted mass functions of cold gas , which we use to predict the number counts of sources likely to be detected by HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders . We find the number counts peak at \sim 4 \times 10 ^ { 3 } / 4 \times 10 ^ { 4 } / 3 \times 10 ^ { 5 } galaxies per square degree at z \sim 0.1/0.2/0.5 for a year long HI hemispheric survey on a 1 % /10 % /100 % SKA with a 30 square degree field of view , corresponding to an integration time of 12 hours . On a full SKA with a 200 square degree field of view ( equivalent to an integration time of 80 hours ) the number counts peak at 5 \times 10 ^ { 5 } galaxies per square degree at z \sim 0.6 . We show also how adopting a conversion factor for cold gas mass to HI mass that varies from galaxy to galaxy impacts on number counts . In addition , we examine how the typical angular sizes of galaxies vary with redshift . These decline strongly with increasing redshift at z \mathrel { \hbox { \hbox to 0.0 pt { \hbox { \lower 4.0 pt \hbox { $ \sim$ } } } \hbox { $ < $ } } } % 0.5 and more gently at z \mathrel { \hbox { \hbox to 0.0 pt { \hbox { \lower 4.0 pt \hbox { $ \sim$ } } } \hbox { $ > $ } } } % 0.5 ; the median angular size varies between 5 ^ { \prime \prime } and 10 ^ { \prime \prime } at z =0.1 , 0.5 ^ { \prime \prime } and 3 ^ { \prime \prime } at z =1 and 0.2 ^ { \prime \prime } and 1 ^ { \prime \prime } at z =3 for galaxies with HI masses in excess of 10 ^ { 9 } h ^ { -1 } M _ { \odot } , depending on the precise model . Taken together , these results make clear that forthcoming HI surveys will provide important and powerful tests of theoretical galaxy formation models .