Time-resolved two-dimensional (2D) infrared spectra of the asymmetric stretch mode of solvated CO2 show distinct features corresponding to ground and excited state thermal populations of the bend modes. The time-dependence of these peaks arises in part from solvent-driven thermal fluctuations in populations of the lower frequency bend modes through their coupling to the higher frequency asymmetric stretch. This observation illustrates the capacity of multidimensional vibrational spectroscopy to reveal details of the interactions among vibrational modes in condensed phases. The optimized mean-trajectory (OMT) method is a trajectory-based semiclassical approach to computing the vibrational response functions of multidimensional spectroscopy from a classical Hamiltonian. We perform an OMT calculation of the 2D vibrational spectrum for two coupled anharmonic modes, with the lower frequency mode undergoing stochastic transitions in energy to mimic solvent-induced fluctuations in quantum populations. The semiclassical calculation reproduces the influence of thermal fluctuations in the low-frequency mode on the 2D spectrum of the high-frequency mode, as in measured spectra of solvated CO2.