We present a new local, analytic scheme for the subtraction of infrared singularities at next-to-next-to-leading order (NNLO) in QCD, which aims at reducing the complexity of the problem and features remarkable aspects. It works for any infrared-safe observable, it benefits from the partition of the radiative phase-space into sectors, the subtraction counterterms are local and can be analytically integrated. All these properties enable an efficient numerical implementation. Our scheme is currently designed for massless final state radiation only, but its extension to initial state radiation will allow for predictions at hadron-hadron and lepton-nucleon colliders.