There are a variety of reasons why in geodesy an improved formula for nutation is needed; related topics of interest are the determination of time-dependences in low degree zonals of the earth gravity field, ocean tide modeling, determination of odd harmonics of gravity field etc. in satellite geodesy. A combined model of deterministic and stochastic components is used in order to evaluate two new nutation series where, in an adjustment, mainly VLBI data (IRIS, GSFC, IERS) have been applied. Contrary to earlier revisions of the present nutation formula, not only the five significantly affected waves (annual, semi-annual, FCN etc.) are corrected but rather all constituents are revised in such a way that white noise residuals result from the adjusted observations, based on the new formula. Still remaining problems (such as the separation of long-period terms from precession etc.) are outlined.