The high-temperature thermodynamic properties of forsterite were reviewed in the light of a new determination of the isobaric heat capacity (C(p)), up to 1850 K, and Raman spectroscopic measurements, up to 1150 K and 10 GPa. The C(p) measurements and available data on thermal expansion (alpha) and bulk modulus (K) show that the isochoric specific heat (C-nu) exceeds the harmonic limit of Dulong and Petit above 1300 K. This intrinsic anharmonic behavior of C-nu) can be modeled by introducing anharmonic parameters a(i) = (partial-lnv(i)/partial-T)V which are calculated from the measured pressure and temperature shifts of the vibrational frequencies. These parameters are all negative, with absolute values lower for the stretching modes of the SiO4 tetrahedra (a(i) almost-equal-to - 1 x 10(-5) K-1) than for the lattice modes (a(i) almost-equal-to 2 x 10(-5) K-1). Through the relation C(p) = C-nu) + alpha-2K(T)VT, the calculated anharmonic C-nu) and the measured C(p) are then used to determine the temperature dependences of the thermal expansion and bulk modulus of forsterite, up to 2000 K, in agreement with recent experimental results. Finally, all these data point to an inconsistency for the Gruneisen parameter of forsterite, whereby the macroscopic parameter gamma = alpha-VK(T)/C-nu) cannot be evaluated simply at high temperature by summation of the individual isothermal mode Gruneisen parameters gamma-iT = K(T) (partial-lnv(i)/partial-P).