Current methods for deriving thermal inertia from spacecraft observations of planetary brightness temperature generally assume that surface properties are uniform for any given observation or co-located set of observations. As a result of this assumption and the nonlinear relationship between temperature and thermal inertia, sub-pixel horizontal heterogeneity may yield different apparent thermal inertia at different times of day or seasons. We examine the effects of horizontal heterogeneity on Mars by modeling the thermal behavior of various idealized mixed surfaces containing differing proportions of either dust, sand, duricrust, and rock or slope facets at different angles and azimuths. Latitudinal effects on mixed-surface thermal behavior are also investigated. We find large (several 100 J m−2 K−1 s−1/2) diurnal and seasonal variations in apparent thermal inertia even for small (∼10%) admixtures of materials with moderately contrasting thermal properties or slope angles. Together with similar results for layered surfaces [Mellon, M.T., Putzig, N.E., 2007. Lunar Planet. Sci. XXXVIII. Abstract 2184], this work shows that the effects of heterogeneity on the thermal behavior of the martian surface are substantial and may be expected to result in large variations in apparent thermal inertia as derived from spacecraft instruments. While our results caution against the over-interpretation of thermal inertia taken from median or average maps or derived from single temperature measurements, they also suggest the possibility of using a suite of apparent thermal inertia values derived from single observations over a range of times of day and seasons to constrain the heterogeneity of the martian surface.