The Velocity Field of Our Milky Way Outer Stellar Halo Based on DESI DR2

Abstract Using 64,000 halo K giants from the Dark Energy Spectroscopic Instrument second Data Release (DR2), we decompose the Milky Way (MW) stellar halo between 3 and 160 kpc into metal-rich (MR) and metal-poor (MP) components via a Gaussian mixture model. The two populations are nearly equal in number but chemically and kinematically distinct: MR stars occupy highly radial orbits with velocity anisotropy of β ≈ 0.94 and metallicity dispersion σ [Fe/H] ≈ 0.17 dex, without obvious dependence on distance, and are mainly contributed by Gaia-Sausage/Enceladus (GSE) debris. The MR component dominates the inner 30 kpc and reemerges beyond 50 kpc, implying GSE debris can extend to ∼70–80 kpc. MP stars exhibit a weaker radial bias of β ≈ 0.46, decreasing to −0.5 beyond 80 kpc, and with a larger metallicity dispersion of σ [Fe/H] ≈ 0.46 dex, showing signatures of multiple minor mergers. Both components exhibit net prograde rotation at ∼10–30 kpc with a stronger azimuthal signal in the MP population. The nonequilibrium motions of the outer halo (>50 kpc) are quantified with a dipole-plus-contraction velocity field. We find that the outer halo is simultaneously contracting ( μ compr = −19 km s −1 , distance-independent) and subject to reflex motions ( μ dipole increases from −19 to −44 km s −1 with radius), reflecting the perturbation from the Large Magellanic Cloud (LMC). We also confirm a linear dependence of mean polar velocity for the outer stellar halo on μ dipole , a direct consequence of the LMC and MW interaction. Our results provide a quantitative distance-resolved description of the MW’s last major accretion event and its ongoing response to the first infall of the LMC.