Land Surface Temperature - 2022

The Council initiated a Climate Vulnerability Assessment in response to the increased frequency of extreme heat and flooding events associated with climate change. The Land Surface Temperature raster dataset is intended for planning and alternatives analysis. This work supports regional climate mitigation and resiliency planning efforts. The Land Surface Temperature (LST) raster dataset is a dataset derived from a satellite image from Landsat 8 to show the land surface temperature for the seven-county metropolitan region of the Twin Cities, in degrees Fahrenheit. In essence, this map shows how hot the ground was to the touch at any given location. This satellite image was taken at 11:59 am CDT on July 22, 2016. At this time, the air temperature was 90F with a heat index of 90.3F, as taken at the Minneapolis-St. Paul International Airport (Minneapolis-St Paul International Airport, 2016). This was the third day of a regional heat wave, and while temperatures overnight had dipped down to around 74F, temperatures had climbed up to a maximum temperature of 97F by around 5-6pm, resulting in the hottest day in roughly three years (Midwestern Regional Climate Center, 2016). The original thermal image was taken at a 100 x 100 meter resolution, but was re-scaled and processed at the 30 x 30 meter scale. Landsat 8 collects image data at a spatial resolution of 30 meters (visible, NIR, SWIR); 100 meters (thermal); and 15 meters (panchromatic). This Landsat 8 raster image was first used to create a NDVI raster dataset. Using that NDVI dataset with Bands 10 and 11 of the Landsat 8 raster dataset, a Land Surface Temperature raster can be derived. That raster is further processed from celsius to fahrenheit, then clipped to the 7-county metropolitan area, and all water bodies removed from the final image.  The LST values without regional water bodies is the primary basis for the Metropolitan Council Climate Vulnerability Assessment Report. Part of this report specifically considers the urban heat island effect, or the effect of human activity and the built environment on increases in urban temperature on human life, and as water has different heat retention properties than most land surfaces, we chose to do our analysis without regional water bodies. It should be noted that the lowest original LST values were water bodies, and so by removing the water bodies from the map the minimum LST value raised 0.6F. This is not particularly significant, except that it would minimally affect which values would be highlighted in the third layer of this package.

The Harmonized Landsat and Sentinel-2 (HLS) project takes input data from the joint NASA/USGS Landsat 8 and Landsat 9 satellites and the European Space Agency's Sentinel-2A, Sentinel-2B, and Sentinel-2C satellites to generate a harmonized, analysis-ready surface reflectance data product with observations every 1.6 days.  Intended Use: The frequent revisit times of HLS allow for time series of land surface applications at the field/plot scale, which supports a wide range of applications, including agricultural health, insect infestations, and natural hazard impacts. The HLSS30 and HLSL30 products are gridded to the same resolution and tiling system and thus are “stackable” for time series analysis. Key Properties: Geographic Coverage: Global Land Temporal Coverage: March 2026 - Present Spatial Resolution: 30-meter Revisit Time: ~1 Day Update Frequency: Daily Product Level: Level-3 Source Data Coordinate System: Universal Transverse Mercator (UTM) WGS84 Service Coordinate System: Web Mercator Auxiliary Sphere (EPSG:3857) Multispectral Bands NameDescriptionUnitsData TypeFill ValueScale FactorBand 1Coastal AerosolN/Aint16-99990.0001Band 2BlueN/Aint16-99990.0001Band 3GreenN/Aint16-99990.0001Band 4RedN/Aint16-99990.0001Band 5NIRN/Aint16-99990.0001Band 6SWIR 1N/Aint16-99990.0001Band 7SWIR 2N/Aint16-99990.0001Band 9CirrusN/Aint16-99990.0001Band 10TIRS 1N/Aint16-99990.01Band 11TIRS 2N/Aint16-99990.01FmaskQuality BitsBit Fielduint8255N/ASAASun Azimuth AngleDegreeuint16400000.01SZASun Zenith AngleDegreeuint16400000.01VAAView Azimuth AngleDegreeuint16400000.01VZAView Zenith AngleDegreeuint16400000.01Processing Templates Server-side renderings available on this layer (select via the Rendering / Image Display options). The default rendering is Natural Color. Band-combination templates are for visualization; index templates (NDVI, NDWI, NBR, NDBI) return values suitable for analysis. "DRA" applies dynamic range adjustment; "Cloud Mask" variants suppress pixels flagged by the Fmask quality band. NDVI Colormap — Normalized difference vegetation index (NDVI) with colormap. Dark green represents vigorous vegetation and brown represents sparse vegetation. It is computed as (b5 - b4) / (b5 + b4) and is suitable for vegetation, land cover and plant health monitoring. (also available: Cloud Mask) NDVI Raw — Normalized difference vegetation index (NDVI) computed as (b5 - b4) / (b5 + b4) and is suitable for vegetation, land cover and plant health monitoring. (also available: Cloud Mask) NDWI Raw — The NDWI index is computed as Green(Band03)-NIR(Band05)/ Green(Band03)+NIR(Band05) and is most suitable for mapping of water bodies, droughts, boundary evaluation. Water bodies have strong absorbability and low radiation within the visible to infrared spectral ranges and this index is based on this phenomenon. (also available: Cloud Mask) Natural Color — Natural Color bands red, green, blue(4, 3, 2) displayed with fixed stretch applied. (also available: Cloud Mask, DRA, DRA and Cloud Mask) Normalized Burn Ratio — The NBR index is computed as NIR(Band05)-SWIR(Band7)/NIR(Band05)+SWIR(Band7) and appropriate for highlighting burned areas and estimating fire severity. It highlights burnt areas in fire zones greater than 500 acres. You will need to generate a differenced NBR image by subtracting a post-fire image from a pre-fire image. The resulting darker pixels represent burned areas. (also available: Cloud Mask) Normalized Difference Built-Up Index — The NDBI index is computed as SWIR(Band6)-NIR(Band05)/ SWIR(Band6)+NIR(Band05) and is suitable for watershed runoff predictions, tracking urbananisation and land-use planning. It highlights urban areas with higher reflectance in the SWIR region, compared to the NIR region. This index also contains a cloud mask. (also available: Cloud Mask) Limitations of Use: HLSL30 v2.0 data contain several known limitations related to atmospheric correction, cloud masking, geolocation, and metadata that users should carefully consider before analysis. A full description of these known issues is available here. Latency: Imagery tiles are configured to render at a 150 meter scale. Dynamic Mosaic rendering can take an extended period of time.  Resource Links: Download Data User Guide General Documentation For inquiries about this service, please contact support@earthdata.nasa.gov or post/view questions on the Earthdata Forum Citation: Neigh, C., Ju, J., Roger, J.-C., Skakun, S., Vermote, E., Claverie, M., Dungan, J., Yin, Z., Freitag, B., & Justice, C. (2021). HLS Operational Land Imager Surface Reflectance and TOA Brightness Daily Global 30m v2.0 [Dataset]. NASA Land Processes Distributed Active Archive Center. https://doi.org/10.5067/HLS/HLSL30.002

The Harmonized Landsat and Sentinel-2 (HLS) project takes input data from the joint NASA/USGS Landsat 8 and Landsat 9 satellites and the European Space Agency's Sentinel-2A, Sentinel-2B, and Sentinel-2C satellites to generate a harmonized, analysis-ready surface reflectance data product with observations every 1.6 days.  Intended Use: The frequent revisit times of HLS allow for time series of land surface applications at the field/plot scale, which supports a wide range of applications, including agricultural health, insect infestations, and natural hazard impacts. The HLSS30 and HLSL30 products are gridded to the same resolution and tiling system and thus are “stackable” for time series analysis. Key Properties: Geographic Coverage: Global Land Temporal Coverage: May 2026 - Present Spatial Resolution: 30-meter Revisit Time: ~1 Day Update Frequency: Daily Product Level: Level-3 Source Data Coordinate System: Universal Transverse Mercator (UTM) WGS84 Service Coordinate System: Web Mercator Auxiliary Sphere (EPSG:3857) Multispectral Bands Name Description Units Data Type Fill Value Scale Factor Band 1Coastal AerosolN/Aint16-99990.0001 Band 2BlueN/Aint16-99990.0001 Band 3GreenN/Aint16-99990.0001 Band 4RedN/Aint16-99990.0001 Band 5Red Edge 1N/Aint16-99990.0001 Band 6Red Edge 2N/Aint16-99990.0001 Band 7Red Edge 3N/Aint16-99990.0001 Band 8NIR BroadN/Aint16-99990.0001 Band 8ANIR NarrowN/Aint16-99990.0001 Band 9Water VaporN/Aint16-99990.0001 Band 10CirrusN/Aint16-99990.0001 Band 11SWIR 1N/Aint16-99990.0001 Band 12SWIR 2N/Aint16-99990.0001 FmaskQuality BitsBit Fielduint8255N/A SAASun Azimuth AngleDegreeuint16400000.01 Processing Templates Server-side renderings available on this layer (select via the Rendering / Image Display options). The default rendering is Natural Color. Band-combination templates are for visualization; index templates (NDVI, NDWI, NBR, NDBI) return values suitable for analysis. "DRA" applies dynamic range adjustment; "Cloud Mask" variants suppress pixels flagged by the Fmask quality band. NDVI Colormap — Normalized difference vegetation index (NDVI) with colormap. Dark green represents vigorous vegetation and brown represents sparse vegetation. It is computed as (b5 - b4) / (b5 + b4) and is suitable for vegetation, land cover and plant health monitoring. (also available: Cloud Mask) NDVI Raw — Normalized difference vegetation index (NDVI) computed as (b5 - b4) / (b5 + b4) and is suitable for vegetation, land cover and plant health monitoring. (also available: Cloud Mask) NDWI Raw — The NDWI index is computed as Green(Band03)-NIR(Band05)/ Green(Band03)+NIR(Band05) and is most suitable for mapping of water bodies, droughts, boundary evaluation. Water bodies have strong absorbability and low radiation within the visible to infrared spectral ranges and this index is based on this phenomenon. (also available: Cloud Mask) Natural Color — Natural Color bands red, green, blue(4, 3, 2) displayed with fixed stretch applied. (also available: Cloud Mask, DRA, DRA and Cloud Mask) Normalized Burn Ratio — The NBR index is computed as NIR(Band05)-SWIR(Band7)/NIR(Band05)+SWIR(Band7) and appropriate for highlighting burned areas and estimating fire severity. It highlights burnt areas in fire zones greater than 500 acres. You will need to generate a differenced NBR image by subtracting a post-fire image from a pre-fire image. The resulting darker pixels represent burned areas. (also available: Cloud Mask) Normalized Difference Built-Up Index — The NDBI index is computed as SWIR(Band6)-NIR(Band05)/ SWIR(Band6)+NIR(Band05) and is suitable for watershed runoff predictions, tracking urbananisation and land-use planning. It highlights urban areas with higher reflectance in the SWIR region, compared to the NIR region. This index also contains a cloud mask. (also available: Cloud Mask) Limitations of Use: HLSS30 v2.0 data contain several known limitations related to atmospheric correction, cloud masking, geolocation, and metadata that users should carefully consider before analysis. A full description of these known issues is available here. Latency: Imagery tiles are configured to render at a 150 meter scale. Dynamic Mosaic rendering can take an extended period of time.  Resource Links: Download Data User Guide General Documentation For inquiries about this service, please contact support@earthdata.nasa.gov or post/view questions on the Earthdata Forum Citation: Neigh, C., Ju, J., Roger, J.-C., Skakun, S., Vermote, E., Claverie, M., Dungan, J., Yin, Z., Freitag, B., & Justice, C. (2021). HLS Sentinel-2 Multi-spectral Instrument Surface Reflectance Daily Global 30m v2.0 [Dataset]. NASA Land Processes Distributed Active Archive Center. https://doi.org/10.5067/HLS/HLSS30.002.

The NASA Making Earth System Data Records for Use in Research Environments (MEaSUREs) LEOLSTCMG30 version 2 Climate Modeling Grid (CMG) product provides Land Surface Temperature (LST) derived from the Low Earth Orbit (LEO) satellite data record from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS) instruments as well as LST error estimates for both day and night. The MEaSUREs LEOLST product is generated by regridding the monthly CMG products from Aqua MODIS (MYD21C3) and VIIRS (VNP21C3 and VJ121). The product is available on 0.25, 0.5, and 1 degree optimized climate grids with well characterized per-pixel uncertainties. A low-resolution browse is also available showing LST as an RGB (red, green, blue) image in PNG format.

This Collection provides Sentinel-3 SLSTR Level-2 Land Surface Temperature products containing data on land surface temperature measurements on a 1km grid. Radiance is measured in two channels to determine the temperature of the Earth's surface skin in the instrument field of view, where the term "skin" refers to the top surface of bare soil or the effective emitting temperature of vegetation canopies as viewed from above.

The Landsat Surface Temperature (ST) product represents the temperature of the Earth's surface in Kelvin (K).