Two-dimensional (2D) line shapes derived from four types of coherent third-order nonlinear spectroscopies as a function of two independent time or frequency variables are compared. The signals scattered into two wave-vector-matching directions in the limits of high-time or high-frequency resolution are calculated for an inhomogeneously broadened ensemble of two-level systems with phenomenological dephasing and population relaxation times. While one-dimensional (1D) line shape analysis for this model is ambiguous, the contours of these 2D line shapes are shown to be characteristic of the relative time scales of the three relevant line-broadening parameters. These broadening mechanisms are also apparent by comparing 2D spectral profiles along the diagonal and antidiagonal frequency axes. For time domain experiments, in analogy with NMR correlation spectroscopy, the radiated signal is observed as a function of the two coherence periods and is 2D Fourier transformed to obtain a line shape. For experiments based on the photon echo, the ellipticity of the absolute value 2D line shape is related to the ratio of the inhomogeneous width to the dephasing rate, whereas experiments based on the transient grating method are shown to be functionally 1D. For frequency domain experiments, the contours of the 2D line shape can yield information on static and both dynamic broadening mechanisms. While the addition of spectral diffusion to such calculations will modify these results for frequency domain experiments, characterizing these dynamics in the time domain naturally leads to a 3D experiment analogous to the spin diffusion experiment in NMR.