Addendum to Global constraints on absolute neutrino masses and their ordering

We revisit our previous work [Capozzi et al., 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) similar or equal to 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) similar or equal to Sigma/3 similar or equal to 0.23 eV, at the sensitivity frontier of the KATRIN experiment.