Recently significant attention has been paid to magnetic field dependent photoluminescence (PL) features of the negatively charged nitrogen-vacancy (NV) centers in diamond. These features are used for microwavefree sensing and are indicative of the spin-bath properties in the diamond sample. Examining the temperature dependence of the PL features allows identifying both temperature-dependent and -independent features, of potential utility in diamond-based quantum sensing and dynamic nuclear polarization applications. Here we study the thermal variability of many different features visible in a wide range of magnetic fields. To this end, we first discuss the origin of the features and tentatively assign the previously unidentified features to cross relaxation of NV center containing multispin systems. The experimental results are compared with theoretically predicted temperature shifts deduced from a combination of thermal expansion and electron-phonon interactions. A deeper insight into the thermal behavior of a wide array of the features may come with important consequences for various applications in high-precision NV thermometry, gyroscopes, solid-state clocks, and biomagnetic measurements.