New data on the variation of the thermal expansion of ice Ih with temperature at ambient pressure together with new evaluations of the bulk modulus and earlier data for the heat capacity provide the basis for a coherent thermodynamic modeling of the main thermophysical properties of ice Ih over its whole range of stability. The quasi-harmonic approximation with one Debye term and seven Einstein terms, together with explicit anharmonicity, represents the dominant contribution next to minor “anomalies” from hydrogen ordering and lattice defects. The model accurately fits the main features of all experimental data and provides a basis for the comparison with earlier determinations of the phonon density of states and the Grüneisen parameters.