Рассматривается небесно-механическая модель колебательного процесса земного полюса на основе пространственного варианта задачи «деформируемая Земля – Луна» в поле притягивающего центра (Солнца). В модели учитываются как основные возмущения с большими амплитудами, так и более сложные мелкомасштабные свойства движения, обусловленные короткопериодическими возмущениями Луны с комбинационными частотами, подтверждаемыми наблюдениями Международной службы вращения Земли (МСВЗ). На базе астрометрических данных МСВЗ проведено численное моделирование колебаний земного полюса – выбор опорных функций и оценивание неизвестных параметров модели на различных интервалах времени.