Broadband transient absorption and two-dimensional electronic spectroscopy (2DES) studies of methylene blue in aqueous solution are reported. By isolating the coherent oscillations of the nonlinear signal amplitude and Fourier transforming with respect to the population time, we analyzed a significant number of coherences in the frequency domain and compared them with predictions of the vibronic spectrum from density function theory (DFT) calculations. We show here that such a comparison enables reliable assignments of vibrational coherences to particular vibrational modes, with their constituent combination bands and overtones also being identified via Franck-Condon analysis aided by DFT. Evaluation of the Fourier transform (FT) spectrum of transient absorption recorded to picosecond population times, in coincidence with 2D oscillation maps that disperse the FT spectrum into the additional excitation axis, is shown to be a complementary approach toward detailed coherence determination. Using the Franck-Condon overlap integrals determined from DFT calculations, we modeled 2D oscillation maps up to two vibrational quanta in the ground and excited state (six-level model), showing agreement with experiment. This semiquantitative analysis is used to interpret the geometry change upon photoexcitation as an expansion of the central sulfur/nitrogen containing ring due to the increased antibonding character in the excited state.