It has long been known that the vertical motions of Galactic disk stars increase with stellar age, commonly interpreted as vertical heating through orbit scattering. Here we map the vertical actions of disk stars as a function of age ({tau}<=8Gyr) and across a large range of Galactocentric radii, R_GC_, drawing on APOGEE and Gaia data. We fit J_z_(R_GC,{tau}) as a combination of the vertical action at birth, J_z,0, and the subsequent heating {Delta}J_z,1Gyr_(R_GC_), which scales as {tau}^{gamma}(R_GC_)^. The inferred birth temperature, J_z,0_(R_GC_) is 1kpc.km/s for 3kpc<R_GC_<10kpc, consistent with the ISM velocity dispersion, but it rapidly rises outward, to 8kpc.km/s for R_GC_=14kpc, likely reflecting the stars' birth in a warped or flared gas disk. We find the heating rate {Delta}J_z,1Gyr_ to be modest and nearly constant across all radii, 1.6kpc.km/s/Gyr. The stellar age dependence {gamma} gently grows with Galactocentric radius, from {gamma}~1 for R_GC_~1) than the time dependence theoretically expected from orbit scattering, J_z_{propto}t^0.5^. We illustrate how this conundrum can be resolved if we also account for the fact that at earlier epochs, the scatterers were more common, and the restoring force from the stellar disk surface mass density was low. Our analysis may reinstate gradual orbital scattering as a plausible and viable mechanism to explain the age-dependent vertical motions of disk stars.

Cone search capability for table J/ApJ/878/21/table1 (APOGEE red clump value-added catalog)