Two-dimensional vibrational and electronic spectroscopic observables of isotropically oriented molecular samples in solution are sensitive to laser field intensities and polarization. The third-order response function formalism predicts a signal that grows linearly with the field strength of each laser pulse, thus lacking a way of accounting for non-trivial intensity-dependent effects, such as saturation and finite bleaching. An analytical expression to describe the orientational part of the molecular response, which, in the weak-field limit, becomes equivalent to a four-point correlation function, is presented. Such an expression is evaluated for Liouville-space pathways accounting for diagonal and cross peaks for all-parallel and cross-polarized pulse sequences, in both the weak- and strong-field conditions, via truncation of a Taylor series expansion at different orders. The results obtained in the strong-field conditions suggest how a careful analysis of two-dimensional spectroscopic experimental data should include laser pulse intensity considerations when determining molecular internal coordinates.