Features at the Sun’s surface and atmosphere are constantly changing due to its magnetic field. The McIntosh Archive provides a long-term (45 yr) record of these features, digitized from hand-drawn synoptic maps by Patrick McIntosh. Utilizing this data, we create stack plots for coronal holes, i.e., Hovmöller-type plots of latitude bands, for all longitudes, stacked in time, allowing tracking of coronal hole movement. Using a newly developed two-step method of centroid calculation, which includes a Fourier descriptor to represent a coronal hole’s boundary and calculate the centroid by the use of Green’s theorem, we calculate the centroids of 31 unique, long-lived equatorial coronal holes for successive Carrington rotations during the entire solar cycle 23, and estimate their slopes (time versus longitude) as the coronal holes evolve. We compute coronal hole centroid drift speeds from these slopes, and find an eastward (prograde) pattern that is actually retrograde with respect to the local differential rotation. By discussing the plausible physical mechanisms which could cause these long-lived equatorial coronal holes to drift retrograde, we identify either classical or magnetically modified westward-propagating solar Rossby waves, with a speed of a few tens to a few hundreds of meters per second, to be the best candidate for governing the drift of deep-rooted, long-lived equatorial coronal holes. To explore plausible physics of why long-lived equatorial coronal holes appear few in number during solar minimum/early rising phase more statistics are required, which will be studied in future.