Besides soil hydrology and snow processes, the NASA Catchment Land Surface Model (CLSM) simulates soil temperature in six layers from the surface down to 13m depth. In this study, to examine CLSM's treatment of subsurface thermodynamics, a baseline simulation produced subsurface temperatures for 1980-2014 across Alaska at 9-km resolution. The results were evaluated using in situ observations from permafrost sites across Alaska. The baseline simulation was found to capture the broad features of inter- and intra-annual variations in soil temperature. Additional model experiments revealed that: (i) the representativeness of local meteorological forcing limits the model's ability to accurately reproduce soil temperature, and (ii) vegetation heterogeneity has a profound influence on subsurface thermodynamics via impacts on the snow physics and energy exchange at surface. Specifically, the profile-average RMSE for soil temperature was reduced from 2.96°C to 2.10°C at one site and from 2.38°C to 2.25°C at another by using local forcing and land cover, respectively. Moreover, accounting for the influence of soil organic carbon on the soil thermal properties in CLSM leads to further improvements in profile-average soil temperature RMSE, with reductions of 16% to 56% across the different study sites. The mean bias of climatological ALT is reduced by 36% to 89%, and the RMSE is reduced by 11% to 47%. Finally, results reveal that at some sites it may be essential to include a purely organic soil layer to obtain, in conjunction with vegetation and snow effects, a realistic "buffer zone" between the atmospheric forcing and soil thermal processes.