Off-resonance rotating frame technique offers a novel tool to explore the dynamics of paramagnetic agents at high magnetic fields (B0 > 3T). Based on the effect of paramagnetic relaxation enhancement in the off-resonance rotating frame, a new method is described here for determining the dynamics of paramagnetic ion chelates from the residual z-magnetizations of water protons. In this method, the dynamics of the chelates are identified by the difference magnetization profiles, which are the subtraction of the residual z-magnetization as a function of frequency offset obtained at two sets of RF amplitude omega(1) and pulse duration tau. The choices of omega(1) and tau are guided by a 2-D magnetization map that is created numerically by plotting the residual z-magnetization as a function of effective field angle theta and off-resonance pulse duration tau. From the region of magnetization map that is the most sensitive to the alteration of the paramagnetic relaxation enhancement efficiency R(1rho)/R1, the ratio of the off-resonance rotating frame relaxation rate constant R(1rho) verse the laboratory frame relaxation rate constant R(1), three types of difference magnetization profiles can be generated. The magnetization map and the difference magnetization profiles are correlated with the rotational correlation time tauR of Gd-DTPA through numerical simulations, and further validated by the experimental data for a series of macromolecule conjugated Gd-DTPA in aqueous solutions. Effects of hydration water number q, diffusion coefficient D, magnetic field strength B0 and multiple rotational correlation times are explored with the simulations of the magnetization map. This method not only provides a simple and reliable approach to determine the dynamics of paramagnetic labeling of molecular/cellular events at high magnetic fields, but also a new strategy for spectral editing in NMR/MRI based on the dynamics of paramagnetic labeling in vivo.