Some of the radio sources used by VLBI to materialize the celestial reference frame are known to have apparent motions at the sub - milliarcsecond level , particularly those observed in the early years . On the other hand , state of the art precession - nutation models match the observations at this same level . We investigate to what extent the source instability may contaminate the VLBI determination of precession and nutation corrections in several frequency domains .