Abstract: We present robust, model-marginalized limits on both the total neutrino mass (E m1,) and abundances (Neff) to minimize the role of parametrizations, priors and models when extracting neutrino properties from cosmology. The cosmological observations we consider are cosmic microwave background temperature fluctuation and polarization measurements, supernovae Ia luminosity distances, baryon acoustic oscillation observations and determinations of the growth rate parameter from the Data Release 16 of the Sloan Digital Sky Survey IV. The degenerate neutrino mass spectrum (which implies the prior sigma m(1), > 0) is weakly or moderately preferred over the normal and inverted hierarchy possibilities, which imply the priors sigma m(1), > 0.06 and sigma m(1), > 0.1 eV respectively. Concerning the underlying cosmological model, the ACDM minimal scenario is almost always strongly preferred over the possible extensions explored here. The most constraining 95% CL bound on the total neutrino mass in the ACDM + sigma m(1), picture is sigma m(1), < 0.087 eV. The parameter N-eff is restricted to 3.08 +/- 0.17 (68% CL) in the ACDM + Neff model. These limits barely change when considering the ACDM + sigma m(1), + Neff scenario. Given the robustness and the strong constraining power of the cosmological measurements employed here, the model -marginalized posteriors obtained considering a large spectra of nonminimal cosmologies are very close to the previous bounds, obtained within the ACDM framework in the degenerate neutrino mass spectrum. Future cosmological measurements may improve the current Bayesian evidence favoring the degenerate neutrino mass spectra, challenging therefore the consistency between cosmological neutrino mass bounds and oscillation neutrino measurements, and potentially suggesting a more complicated cosmological model and/or neutrino sector.