Missing higher order uncertainties (MHOU) in perturbative computations are usually estimated by varying the unphysical scales present in the process. However, it is known that scale variation prescriptions often underestimate the actual uncertainty. In this talk, I will present an approach to approximate the unknown next-to-next-to-leading order (NNLO) transverse momentum distribution of colourless final states. The approximation relies on the combination of various resummation formalisms, namely threshold, small-pt and high energy resummations, by exploiting the singularity structure of the large logarithms in Mellin space. I will show that for the case of the Higgs boson production, for instance, the approximate NNLO transverse momentum distribution amounts to a correction of a few percent with respect to the NLO result with a reduction in the scale dependence.