Reproduction of the results of Fig. 3 from Mellon (1997). Top panel: temperature over a Martian year for the Viking 2 landing site at a depth of 0.5 m. The solid curve is a digitisation of the Mellon (1997) data with the grey shading representing a ± 5 K uncertainty. Bottom panel: induced thermal stresses from Mellon (1997; solid) and from our code (dashed), including the uncertainty from the temperature digitalisation. RMS residuals between the Mellon (1997) results and ours are 746 kPa.