Решается задача повышения точности моделирования движения земного полюса. Известно, что наблюдаемые вариации параметров основных колебательных составляющих (чандлеровской и годичной) содержат синфазную с прецессией лунной орбиты (∼18.61 лет) компоненту, не объясняемую в рамках стандартных моделей с геофизическими возмущениями. Для учета этого эффекта предложена уточненная динамическая модель, уравнения которой представляют собой систему дифференциальных уравнений с периодическими коэффициентами, зависящими от долготы восходящего узла орбиты Луны. С помощью численных расчетов на основе данных Международной службы вращения Земли (МСВЗ, IERS) за период 1976–2025 гг. определены оптимальные параметры модели: коэффициент влияния лунного узла χ=0.07 и добротность Q=63. Показано, что учет долгопериодического лунного возмущения позволяет снизить среднеквадратичное отклонение модели от наблюдаемых данных. На тестовых расчетах точность определения положения полюса повышается на величину, соответствующую 3.7 см на поверхности Земли, а максимальный эффект может достигать 5 см. Полученные результаты подтверждают необходимость явного учета долгопериодических вариаций, связанных с движением Луны, в высокоточных моделях движения полюса.