In a thermal gradient, surface nanostructures have been experimentally observed to move due to thermomigration. However, analytical models that describe the thermomigration force acting on surfaces are still controversial. In this work, we start from a thermodynamic approach based on the Massieu function, which is used to describe thermomigration of single adatoms, to develop an expression for the velocity of thermomigrating 2D holes. The model can be simplified in two limiting cases: (i) When the hole motion is limited by adspecies diffusion, the velocity is independent from the hole size (as in our experiments). (ii) If the hole motion is limited by the attachment or detachment of species to or from steps, then the velocity is proportional to the hole width. We have studied the thermomigration of 2D monatomic deep holes on Si(100) using low energy electron microscopy. From the velocity measurements taken at different temperatures, we find, using our model, that the sum of the migration energy and the adatom creation energy is 1.95 ± 0.16 eV. This value is consistent with those found by other authors, reinforcing the validity of our thermomigration model.