Addendum to: Global constraints on absolute neutrino masses and their ordering

We revisit our previous work [Phys. Rev. D 95, 096014 (2017)] where neutrino oscillation and nonoscillation data were analyzed in the standard framework with three neutrino families, in order to constrain their absolute masses and to probe their ordering (either normal, NO, or inverted, IO). We include updated oscillation results to discuss best fits and allowed ranges for the two squared mass differences $\delta m^2$ and $\Delta m^2$, the three mixing angles $\theta_{12}$, $\theta_{23}$ and $\theta_{13}$, as well as constraints on the CP-violating phase $\delta$, plus significant indications in favor of NO vs IO at the level of $\Delta\chi^2=10.0$. We then consider nonoscillation data from beta decay, from neutrinoless double beta decay (if neutrinos are Majorana), and from various cosmological input variants (in the data or the model) leading to results dubbed as default, aggressive, and conservative. In the default option, we obtain from nonoscillation data an extra contribution $\Delta\chi^2 = 2.2$ in favor of NO, and an upper bound on the sum of neutrino masses $\Sigma < 0.15$ eV at $2\sigma$; both results - dominated by cosmology - can be strengthened or weakened by using more aggressive or conservative options, respectively. Taking into account such variations, we find that the combination of all (oscillation and nonoscillation) neutrino data favors NO at the level of $3.2-3.7\sigma$, and that $\Sigma$ is constrained at the $2\sigma$ level within $\Sigma < 0.12-0.69$ eV. The upper edge of this allowed range corresponds to an effective $\beta$-decay neutrino mass $m_\beta = \Sigma/3 = 0.23$ eV, at the sensitivity frontier of the KATRIN experiment.