Whether or not anomalies in the thermal conductivity from insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity κ and specific heat cv of the half-filled 2D single-band Hubbard model using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. Both cv and κ possess two peaks as a function of temperature, with scales related to the Hubbard interaction energy U and spin superexchange energy J, respectively. At low temperatures where the charge degrees of freedom are gapped-out, our results for the contribution to both cv and the Drude weight associated with κ from the kinetic energy agree well with spin-wave theory for the spin-1 2 antiferromagnetic Heisenberg model.