We first focus on the most prominent feature at a beating energy of 4.40 meV, which corresponds to the mode M2. Qualitatively, it can be noted that as the temperature increases the mode amplitude decays at a faster rate. This can be qualitatively explained by the increased probability of elastic many-body interactions with a thermally generated incoherent phonon population. We also observe that the relative amplitude of M2 decreases at higher temperatures in comparison to the other modes, which can also be observed in the temperature-dependent resonant Raman spectra (Fig. ). On the other hand, mode M1 (2.61 meV) presents a very different scenario in which its dynamics are not sensitive to the temperature, at least within the temperature range of the measurement. This implies a significantly different anharmonic interaction potential for M1 compared to M2. To further highlight the differences between M1 and M2, kinetic traces at their energies, at selected temperatures, are shown at Figs. ) and (b). The amplitude of M2 decays exponentially with clearly temperature-dependent dynamics. The effective lifetime of this mode can be estimated via mono-exponential fits; the inverse of the decay time constant is plotted in Fig. ) for all the temperatures. The amplitude of M1, on the contrary, exhibits a non-exponential, temperature-independent decay. An effective lifetime can be estimated, for the sake of comparison, as the time taken for the signal to reach 1/e of its maximum value and the corresponding values of dephasing rate is shown in Fig. ). The difference in the temperature dependence of M1 and M2 is evident here. Additionally, in Fig. , we show the time-integrated spectrum between 5-10 and 10-15 ps. From this data it is clear that M1 dominates at longer times and at higher temperatures.