This map shows areas of anomalous surface temperature in northern Saguache Counties identified from ASTER and LANDSAT thermal data and spatial based insolation model. The temperature for the ASTER data was calculated using the Emissivity Normalization Algorithm that separate temperature from emissivity. The incoming solar radiation was calculated using spatial based insolation model developed by Fu and Rich (1999). Then the temperature due to solar radiation was calculated using emissivity derived from ASTER data. The residual temperature, i.e. temperature due to solar radiation subtracted from ASTER temperature was used to identify thermally anomalous areas. Areas having anomalous temperature in the ASTER data are shown in blue diagonal hatch, while areas having anomalous temperature in the LANDSAT data are shown in magenta on the map. Thermal springs and areas with favorable geochemistry are also shown. Springs or wells having non-favorable geochemistry are shown as blue dots.