We develop a cosmological parameter estimation code for (tomographic) angular power spectra analyses of galaxy number counts, for which we include, for the first time, redshift-space distortions (RSD) in the Limber approximation. This allows for a speed-up in computation time, and we emphasise that only angular scales where the Limber approximation is valid are included in our analysis. Our main result shows that a correct modelling of RSD is crucial not to bias cosmological parameter estimation. This happens not only for spectroscopy-detected galaxies, but even in the case of galaxy surveys with photometric redshift estimates. Moreover, a correct implementation of RSD is especially valuable in alleviating the degeneracy between the amplitude of the underlying matter power spectrum and the galaxy bias. We argue that our findings are particularly relevant for present and planned observational campaigns, such as the Euclid satellite or the Square Kilometre Array, which aim at studying the cosmic large-scale structure and trace its growth over a wide range of redshifts and scales.

Following on our purpose of developing a unified pipeline for large-scale structure data analysis with angular (i.e. harmonic-space) power spectra, we now include the weak lensing effect of magnification bias on galaxy clustering in a publicly available, modular parameter estimation code. We thus forecast constraints on the parameters of the concordance cosmological model, dark energy, and modified gravity theories from galaxy clustering tomographic angular power spectra. We find that a correct modelling of magnification is crucial in order not to bias the estimation of cosmological parameters, especially in the case of deep galaxy surveys. Our case study adopts specifications of the Evolutionary Map of the Universe (EMU), which is a full-sky, deep radio-continuum survey, and is expected to probe the Universe up to redshift $z\sim6$. We assume the Limber approximation, and include magnification bias on top of density fluctuations and redshift-space distortions. By restricting our analysis to the regime where the Limber approximation holds true, we significantly minimise the computational time needed, compared to that of the exact calculation. We also show that there is a trend for more biased parameter estimates from neglecting magnification when the redshift bins are very wide. We conclude that this result implies a strong dependence on the lensing contribution, which is an integrated effect and becomes dominant when wide redshift bins are considered. Finally, we note that instead of being considered a contaminant, magnification bias encodes important cosmological information, and its inclusion leads leading to an alleviation of its degeneracy between the galaxy bias and the amplitude normalisation of the matter fluctuations.