Joël Mayor, Marta Spinelli, Gabriella De Lucia, Robert Yates, Alexandre Refregier, Fabio Fontanot, Lizhi Xie, Michaela Hirschmann

Upcoming galaxy surveys with the SKA Observatory will detect neutral hydrogen (HI) across unprecedented volumes, and their scientific return will crucially depend on predictive models for HI observables. In this work, we present a framework to simulate the neutral hydrogen 21cm emission line in such large-scale HI galaxy surveys. This framework is developed as a modular layer that builds on semi-analytical models. In particular we use as bases the Galaxy Evolution and Assembly (GAEA) and L-Galaxies semi-analytical models, coupled to merger trees from the Millennium Simulation suite. We validate our framework against local Universe observations, demonstrating consistency with velocity functions, and generalised Tully-Fisher relations. Predictions based on GAEA and L-Galaxies exhibit mutual consistency despite the distinct underlying physical prescriptions. We construct mock galaxy catalogues that incorporate forward-modelled selection functions, inclination effects, and redshift broadening, reproducing the statistical distributions of HI-selected galaxies in the ALFALFA survey. Finally, we present redshift distribution forecasts for future SKA Observatory HI galaxy surveys. This framework offers a flexible tool for interpreting upcoming large-scale radio surveys and studying HI line observables as cosmological probes.