By carrying out random walk simulations we systematically study the effects of various types of complex molecular dynamics on (2)H NMR experiments in solids. More precisely, we calculate one-dimensional (1D) (2)H NMR spectra and the results of two dimensional (2D) (2)H NMR experiments in time domain, taking into account isotropic as well as highly restricted motions which involve rotational jumps about different finite angles. Although the dynamical models are chosen to mimic the primary and secondary relaxation in supercooled liquids and glasses, we do not intend to describe experimental results quantitatively but rather to show general effects appearing for complex reorientations. We carefully investigate whether 2D (2)H NMR in time domain, which was originally designed to measure correlation times of ultraslow motions (tau >/= 1 ms), can be used to obtain shorter tau, too. It is demonstrated that an extension of the time window to tau >/= 10 &mgrs is possible when dealing with exponential relaxation, but that it will fail if there is a distribution of correlation times G(lgtau). Vice versa, we show that 1D (2)H NMR spectra, usually recorded to look at dynamics with tau in the microsecond regime, are also applicable for studying ultraslow motions provided that the loss of correlation is achieved step by step. Therefore, it is useful to carry out 1D and 2D NMR experiments simultaneously in order to reveal the mechanism of complex molecular motions. In addition, we demonstrate that highly restricted dynamics can be clearly observed in 1D spectra and in 2D NMR in time domain if long solid-echo delays and large evolution times are applied, respectively. Finally, unexpected observations are described which appear in the latter experiment when considering very broad distributions G(lgtau). Because of these effects, time scale and geometry of a considered motion cannot be extracted from a straightforward analysis of experimental results. Copyright 2000 Academic Press.