We present the results of a calorimetric study on iron-garnet polycrystalline samples of Y, Eu, Tm, Sm, and Yb between 0.4 and 4.2\ifmmode^\circ\else\textdegree\fi{}K in zero magnetic field. Results for another sample of YIG between 0.7 and 4.2\ifmmode^\circ\else\textdegree\fi{}K in magnetic fields up to 13 kOe are also presented. In YIG it is shown that in addition to the lattice vibration and acoustic spin-wave-mode contributions, ${C}_{L}$ and ${C}_{\mathrm{acous}}$, there is a residual specific heat $\mathcal{R}(T)$ that is linear with temperature and somewhat dependent on magnetic field. This new term $\mathcal{R}(T)$ is tentatively attributed to surface effects although its magnitude is much larger than that predicted for surface spin-wave modes in a ferromagnetic sintered sample with the same dispersion constant $D$. The data in YIG are compared with previous work, and it is shown that the analysis of the calorimetric data on this compound is more complex than was thought previously, and in particular that data on ${C}_{v}$ in zero field are not suited for obtaining the value of $D$ from the usual $\frac{{C}_{v}}{{T}^{\frac{3}{2}}}$-versus-${T}^{\frac{3}{2}}$ plot. The same comments apply to EuIG, where a residual term, linear in temperature, was also obtained. From the nuclear specific heat of this compound, we obtain an rms magnetic field at the nucleus ${H}_{N}=597\ifmmode\pm\else\textpm\fi{}5$ kOe, in excellent agreement with results from M\ossbauer experiments. In TmIG, the effect of acoustical spin waves on ${C}_{v}$ could be detected, and (as expected) there was no evidence of population of the optical modes. However, no firm conclusions could be reached about a residue, since it was difficult to estimate ${C}_{\mathrm{acous}}$ theoretically. The nuclear contribution $\frac{{C}_{N}{T}^{2}}{R}=(7.8\ifmmode\pm\else\textpm\fi{}0.4)\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4}$ was found to be in good agrement with that predicted from M\ossbauer data. In SmIG, the first excited electronic level was deduced to be $\frac{E}{{k}_{B}}=(38\ifmmode\pm\else\textpm\fi{}3)\ifmmode^\circ\else\textdegree\fi{}$K, in good agreement with predictions from White's theory. The nuclear specific heat of SmIG was analyzed as corresponding to an rms field at the nucleus ${H}_{N}=(2.2\ifmmode\pm\else\textpm\fi{}0.2)\ifmmode\times\else\texttimes\fi{}{10}^{6}$ Oe, while an analysis of M\ossbauer data had given (2.8\ifmmode\pm\else\textpm\fi{}0.15)\ifmmode\times\else\texttimes\fi{}${10}^{6}$ Oe. A broad anomaly in ${C}_{v}$ centered at 2\ifmmode^\circ\else\textdegree\fi{}K remains unexplained, but is possibly due to an unreacted samarium oxide impurity. It is responsible for the uncertainty in $E$ and ${H}_{N}$. In YbIG, a sizable discrepancy remains below 3\ifmmode^\circ\else\textdegree\fi{}K between the theory using the 12 low optical spin-wave modes and the experimental data. Apparently the dispersion curve of the exchange mode is significantly different from that calculated by Tinkham using a simple model. The nuclear specific heat in YbIG, $\frac{{C}_{N}{T}^{2}}{R}=(1.6\ifmmode\pm\else\textpm\fi{}0.2)\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}3}$, was in fair agreement with that predicted from the electron-nucleus coupling parameter and the ground-state electronic moment.