The EPR spectra of nitroxide spin labels have been simulated as a function of microwave field, H1, taking into account both magnetic field modulation and molecular rotation. It is found that the saturation of the second integral, S, of the first harmonic in-phase absorption spectrum is approximated by that predicted for slow-passage conditions, that is, S approximately H1/1 + PH21, in all cases. This result is independent of the degree of inhomogeneous broadening. In general, the fitting parameter, P, depends not only on the T1 and T2 relaxation times, but also on the rate of molecular reorientation and on the modulation frequency. Calibrations for determining the relaxation times are established from the simulations. For a given modulation frequency and molecular reorientation rate, the parameter obtained by fitting the saturation curves is given by 1/P = a + 1/gamma2eT1 . Teff2, where Teff2 is the effective T2. For molecular reorientation frequencies in the range 2 x 10(7)-2 x 10(8) s-1, Teff2 is dominated by the molecular dynamics and is only weakly dependent on the intrinsic T02, allowing a direct estimation of T1. For reorientation frequencies outside this range, the (T1T2) product may be determined from the calibrations. The method is applied to determining relaxation times for spin labels undergoing different rates of rotational reorientation in a variety of environments, including those of biological relevance, and is verified experimentally by the relaxation rate enhancements induced by paramagnetic ions.