Off-shell, transverse-momentum dependent splitting functions can be defined from the high-energy limit of partonic decay amplitudes. Based on these splitting functions, we construct Sudakov form factors and formulate a new parton branching algorithm. We present a first Monte Carlo implementation of the algorithm. We use the numerical results to verify explicitly momentum sum rules for TMD parton distributions.