  We compute the contribution of twist-2 and twist-3 parton distribution
functions to the small-$b$ expansion for transverse momentum dependent (TMD)
distributions at all powers of $b$. The computation is done by the
twist-decomposition method based on the spinor formalism for all eight quark
TMD distributions. The newly computed terms are accompanied by the prefactor
$(M^2b^2)^n$ and represent the target-mass corrections to the resummed
cross-section. For the first time, a non-trivial expression for the
pretzelosity distribution is derived.
