Experiments have been carried at magnetic-field strengths of 9.4, 14.1, and 17.6 T to explore the evolution of intermolecular multiple quantum coherences in the nonlinear regime where the system evolves for times that are much greater than the characteristic time of action of the long-range dipolar field, tau(d). The results show the expected Bessel function form of the recorded signal as a function of time of evolution, with evident zeros and sign changes. As expected, the rate of signal evolution increases at higher-field strengths as a result of the increased equilibrium magnetization. A numerical method for calculating the evolution of magnetization under the action of the distant dipolar field, relaxation, and diffusion that is based on Fourier analysis of the magnetization distribution has been applied to the correlated two-dimensional spectroscopy revamped by asymmetric z-gradient echo detection sequence in the nonlinear regime and shown to produce results that are in good agreement with experimental data acquired at different magnetic fields and rates of spatial modulation. Experiments and simulations have also been used to explore the evolution of magnetization in a mixture of two interacting spin species in the nonlinear regime.