We present analytic formulas that describe the fully-differential production of color-singlet final states in qq̄ and gg annihilation, including all the relevant partonic channels, through NNLO QCD. We work within the nested soft-collinear scheme, which allows the fully local subtraction of infrared divergences.We demonstrate analytic cancellation of soft and collinear poles and present formulas for the finite parts of all integrated subtraction terms. These results provide an important building block for calculating NNLO QCD corrections to arbitrary processes at hadron colliders within the nested soft-collinear subtraction scheme
