The motion of the Earth’s geophysical pole is modeled by numerical integration of the Liouville equations. The dependence of variable parameters of these equations (the period of free nutation T and the quality factor of the shell Q) is examined as a function of time and initial data sets used for constructing geophysical perturbing functions. In particular, we used oceanic perturbing functions constructed from TOPEX/POSEIDON altimetry data. The annual and Chandler components of the pole motion were also modeled. Apart from the subtraction of quadratic trends, initial data were not filtered or averaged. The performed analysis provided model values of T = 425–440 days and Q = 20–60, and Q was found to be somewhat unstable with time.