To summarise, we have shown that a Wigner sampling approach representing the ensemble at different temperatures is indeed able to reproduce the vibrational temperature dependence of the quadrupole splitting. We identified as the root cause of the vibrational component of its temperature dependence that distortions lead to a swap in sign. A quantitative comparison with experimental data is likely hampered by correctly representing different electronic spin ground states upon distortion with DFT. Using a linear combination of normal modes, we have shown which collective nuclear movements are responsible for the observed changes. This provides a deeper understanding of the quadrupole splitting and its temperature dependence in general: the calculated movement explains how the charge distribution is converted from an oblate to a prolate shape and hence rationalises the sign swap for the quadrupole splitting.