We present the first combined oscillation analysis of multiple atmospheric neutrino datasets, featuring data from Super-Kamiokande, IceCube-DeepCore, and KM3NeT/ORCA together with reactor data from Daya Bay.
Such combinations have long been considered infeasible outside experimental collaborations; we demonstrate that a unified physics model can simultaneously describe all datasets with no significant parameter tensions.
Fitting 839,048 events across 1536 bins with 91 parameters, our combined analysis yields competitive measurements of the neutrino mixing parameters, disfavors CP conservation, and prefers the Normal over the Inverted Mass Ordering.
We disfavor the absence of CP violation at ∆χ2 = 8.06 and the Inverted Ordering at ∆χ2 = 9.11.
These preferences are statistically significant at a more than 95% confidence level.
The preference for normal ordering is about three sigma (roughly a 99% confidence level).
This preference is also corroborated by an independent statistically significant preference for a normal ordering from cosmology data that strengthens that preference when cosmology data is combined with terrestrial experimental data. But quantifying the cosmology data preference is challenging because it is cosmology model dependent (see also here).
Cosmology data does, however, consistently favor a lightest neutrino mass eigenstate far less massive than the Katrin direct neutrino mass measurement experiments (by two or three orders of magnitude).
Neutrinoless beta decay experiments (which imply Majorana neutrino mass limits) aren't yet powerful enough to make meaningful statements about neutrino masses relative to other data sources for neutrino masses.

