The detailed interpretation of time-resolved spectroscopic signals in terms of the molecular rearrangement during a photoreaction or a photophysical event is one of the most important challenges of both experimental and theoretical chemistry. Here we simulate a time-resolved fluorescence spectrum of a dye in aqueous solution, the N-methyl-6-oxyquinolinium betaine, and analyze it in terms of far IR and THz frequency contributions, providing a direct connection to specific molecular motions. To obtain this result, we build up an innovative and general approach based on excited state ab-initio molecular dynamics and a wavelet-based time-dependent frequency analysis of nonstationary signals. We obtain a nice agreement with key parameters of the solvent dynamics, such as the total Stokes shift and the Stokes shift relaxation times. As an important finding, we observe a strong change of specific solute-solvent interactions upon the electronic excitation, with the migration of about 1.5 water molecules from the first solvation shell toward the bulk. In spite of this event, the Stokes shift dynamics is ruled by collective solvent motions in the THz and far IR, which guide and modulate the strong rearrangement of the dye microsolvation. By the relaxation of THz and IR contributions to the emission signal, we can follow and understand in detail the molecularity of the process. The protocol presented here is, in principle, transferable to other time-resolved spectroscopic techniques.