Neutrino masses and mixing: Entering the era of subpercent precision

We perform an updated global analysis of the known and unknown parameters of the standard three-neutrino (3ν) framework, using data available at the beginning of 2025. The known oscillation parameters include three mixing angles (θ12,θ23,θ13) and two squared mass gaps, chosen as δm2=m22-m12>0 and Δm2=m32-12(m12+m22), where the discrete parameter α=sign(Δm2) distinguishes normal ordering (NO, α=+1) from inverted ordering (IO, α=-1). With respect to our previous 2021 update, the combination of accelerator, reactor, and atmospheric neutrino data leads to appreciably reduced uncertainties for θ23, θ13, and |Δm2|. In particular, |Δm2| is the first 3ν parameter to enter the domain of subpercent precision (0.8% at 1σ). We underline some issues about common systematics in combined fits that might affect (and possibly weaken) this error estimate. Concerning oscillation unknowns, we find a relatively weak preference for NO versus IO (at 2.2σ), for CP violation versus conservation in NO (1.3σ), and for the first θ23 octant versus the second in NO (1.1σ). We discuss the current status and qualitative prospects of the mass ordering hint in the plane charted by the mass parameters (δm2,Δmee2), where Δmee2=|Δm2|+12α(cos2θ12-sin2θ12)δm2, to be jointly measured by the JUNO experiment with subpercent precision. We also discuss upper bounds on nonoscillation observables, including the effective νe mass mβ in β decay, the effective Majorana mass mββ in 0νββ decay, and the sum of neutrino masses ς in cosmology. We adopt mβ<0.50 eV (2σ) from current H3 data and report mββ<0.086 eV (2σ) from a combined Ge76, Te130, and Xe136 data analysis, accounting for parametrized nuclear matrix element covariances. Concerning ς, current results show tensions within the standard Λ cold dark matter (ΛCDM) cosmological model, pulling ς toward unphysical values and suggesting possible model extensions. We discuss representative combinations of data, with or without augmenting the ΛCDM model with extra parameters accounting for possible systematics (lensing anomaly) or new physics (dynamical dark energy). The resulting 2σ upper limits are roughly spread around the bound ς<0.2 eV within a factor of 3 (both up-and downward), with different implications for NO and IO scenarios. Bounds from oscillation and nonoscillation data are also discussed in the planes charted by pairs of (mβ,mββ,ς) parameters.