The role of two skeletal (C=C and C=O stretch) and O-H vibrational motions in the internal conversion dynamics associated with the coupled S-1(pi pi*, A') -S-1(n pi*, A '') potential energy surfaces of alizarin are investigated theoretically. Quantum wavepacket dynamics simulations reveal a nonadiabatic population transfer from the "bright" S-1(pi pi*, A') to "dark" S-2(n pi*, A '') state on a time scale of 10 fs. A detailed analysis of computed structural parameters, energetics, and time-dependent observables suggest that these vibrations promote the nonadiabatic dynamics before initiating the proton transfer process. We also discuss how the simultaneous evolution of multidimensional dynamics involving several vibrational degrees of freedom would increase the complexity, while analyzing the spectral and kinetic data of time-resolved spectroscopy measurements.