Analysis Image Service generated from Extract Raster Data

This dataset is available for download from: Wetlands (File Geodatabase).

Wetlands in California are protected by several federal and state laws, regulations, and policies. This layer was extracted from the broader land cover raster from the CA Nature project which was recently enhanced to include a more comprehensive definition of wetland. This wetlands dataset is used as an exclusion as part of the biological planning priorities in the CEC 2023 Land-Use Screens.

This layer is featured in the CEC 2023 Land-Use Screens for Electric System Planning data viewer.

For more information about this layer and its use in electric system planning, please refer to the Land Use Screens Staff Report in the CEC Energy Planning Library.

Change Log

Version 1.1 (January 26, 2023)

• Full resolution of wetlands replaced a coarser resolution version that was previously shared. Also, file type changed from polygon to raster (feature service to tile layer service).

Spatial analysis and statistical summaries of the Protected Areas Database of the United States (PAD-US) provide land managers and decision makers with a general assessment of management intent for biodiversity protection, natural resource management, and recreation access across the nation. The PAD-US 3.0 Combined Fee, Designation, Easement feature class in the full geodatabase inventory (with Military Lands and Tribal Areas from the Proclamation and Other Planning Boundaries feature class) was modified to prioritize overlapping designations, avoiding massive overestimation in protected area statistics, and simplified by the following PAD-US attributes to support user needs for raster analysis data: Manager Type, Manager Name, Designation Type, GAP Status Code, Public Access, and State Name. The rasterization process (see processing steps below) prioritized overlapping designations previously identified (GAP_Prity field) in the Vector Analysis File (e.g. Wilderness within a National Forest) based upon their relative biodiversity conservation (e.g. GAP Status Code 1 over 2). The 30-meter Image (IMG) grid Raster Analysis Files area extents were defined by the Census state boundary file used to clip the Vector Analysis File, the data source for rasterization ("PADUS3_0VectorAnalysis_State_Clip_CENSUS2020" feature class from ("PADUS3_0VectorAnalysisFileOtherExtents_Clip_Census.gdb"). Alaska (AK) and Hawaii (HI) raster data are separated from the contiguous U.S. (CONUS) to facilitate analyses at manageable scales. Note, the PAD-US inventory is now considered functionally complete with the vast majority of land protection types (with a legal protection mechanism) represented in some manner, while work continues to maintain updates, improve data quality, and integrate new data as it becomes available (see inventory completeness estimates at: http://www.protectedlands.net/data-stewards/ ). In addition, protection status represents a point-in-time and changes in status between versions of PAD-US may be attributed to improving the completeness and accuracy of the spatial data more than actual management actions or new acquisitions. USGS provides no legal warranty for the use of this data. While PAD-US is the official aggregation of protected areas ( https://www.fgdc.gov/ngda-reports/NGDA_Datasets.html ), agencies are the best source of their lands data.

A clip of impervious surfaces to only include what is within the Yukon River Drainage. Original data was pulled from the impervious surface index, and was clipped to the extent of the drainage, then converted from a raster to polygon.

Notice: this is the latest Heat Island Severity image service.This layer contains the relative heat severity for every pixel for every city in the United States, including Alaska, Hawaii, and Puerto Rico. This 30-meter raster was derived from Landsat 8 imagery band 10 (ground-level thermal sensor) from the summer of 2023.To explore previous versions of the data, visit the links below:Heat Severity - USA 2022Heat Severity - USA 2021Heat Severity - USA 2020Heat Severity - USA 2019Federal statistics over a 30-year period show extreme heat is the leading cause of weather-related deaths in the United States. Extreme heat exacerbated by urban heat islands can lead to increased respiratory difficulties, heat exhaustion, and heat stroke. These heat impacts significantly affect the most vulnerable—children, the elderly, and those with preexisting conditions.The purpose of this layer is to show where certain areas of cities are hotter than the average temperature for that same city as a whole. Severity is measured on a scale of 1 to 5, with 1 being a relatively mild heat area (slightly above the mean for the city), and 5 being a severe heat area (significantly above the mean for the city). The absolute heat above mean values are classified into these 5 classes using the Jenks Natural Breaks classification method, which seeks to reduce the variance within classes and maximize the variance between classes. Knowing where areas of high heat are located can help a city government plan for mitigation strategies.This dataset represents a snapshot in time. It will be updated yearly, but is static between updates. It does not take into account changes in heat during a single day, for example, from building shadows moving. The thermal readings detected by the Landsat 8 sensor are surface-level, whether that surface is the ground or the top of a building. Although there is strong correlation between surface temperature and air temperature, they are not the same. We believe that this is useful at the national level, and for cities that don’t have the ability to conduct their own hyper local temperature survey. Where local data is available, it may be more accurate than this dataset. Dataset SummaryThis dataset was developed using proprietary Python code developed at Trust for Public Land, running on the Descartes Labs platform through the Descartes Labs API for Python. The Descartes Labs platform allows for extremely fast retrieval and processing of imagery, which makes it possible to produce heat island data for all cities in the United States in a relatively short amount of time.What can you do with this layer?This layer has query, identify, and export image services available. Since it is served as an image service, it is not necessary to download the data; the service itself is data that can be used directly in any Esri geoprocessing tool that accepts raster data as input.In order to click on the image service and see the raw pixel values in a map viewer, you must be signed in to ArcGIS Online, then Enable Pop-Ups and Configure Pop-Ups.Using the Urban Heat Island (UHI) Image ServicesThe data is made available as an image service. There is a processing template applied that supplies the yellow-to-red or blue-to-red color ramp, but once this processing template is removed (you can do this in ArcGIS Pro or ArcGIS Desktop, or in QGIS), the actual data values come through the service and can be used directly in a geoprocessing tool (for example, to extract an area of interest). Following are instructions for doing this in Pro.In ArcGIS Pro, in a Map view, in the Catalog window, click on Portal. In the Portal window, click on the far-right icon representing Living Atlas. Search on the acronyms “tpl” and “uhi”. The results returned will be the UHI image services. Right click on a result and select “Add to current map” from the context menu. When the image service is added to the map, right-click on it in the map view, and select Properties. In the Properties window, select Processing Templates. On the drop-down menu at the top of the window, the default Processing Template is either a yellow-to-red ramp or a blue-to-red ramp. Click the drop-down, and select “None”, then “OK”. Now you will have the actual pixel values displayed in the map, and available to any geoprocessing tool that takes a raster as input. Below is a screenshot of ArcGIS Pro with a UHI image service loaded, color ramp removed, and symbology changed back to a yellow-to-red ramp (a classified renderer can also be used): A typical operation at this point is to clip out your area of interest. To do this, add your polygon shapefile or feature class to the map view, and use the Clip Raster tool to export your area of interest as a geoTIFF raster (file extension ".tif"). In the environments tab for the Clip Raster tool, click the dropdown for "Extent" and select "Same as Layer:", and select the name of your polygon. If you then need to convert the output raster to a polygon shapefile or feature class, run the Raster to Polygon tool, and select "Value" as the field.Other Sources of Heat Island InformationPlease see these websites for valuable information on heat islands and to learn about exciting new heat island research being led by scientists across the country:EPA’s Heat Island Resource CenterDr. Ladd Keith, University of ArizonaDr. Ben McMahan, University of Arizona Dr. Jeremy Hoffman, Science Museum of Virginia Dr. Hunter Jones, NOAA Daphne Lundi, Senior Policy Advisor, NYC Mayor's Office of Recovery and ResiliencyDisclaimer/FeedbackWith nearly 14,000 cities represented, checking each city's heat island raster for quality assurance would be prohibitively time-consuming, so Trust for Public Land checked a statistically significant sample size for data quality. The sample passed all quality checks, with about 98.5% of the output cities error-free, but there could be instances where the user finds errors in the data. These errors will most likely take the form of a line of discontinuity where there is no city boundary; this type of error is caused by large temperature differences in two adjacent Landsat scenes, so the discontinuity occurs along scene boundaries (see figure below). Trust for Public Land would appreciate feedback on these errors so that version 2 of the national UHI dataset can be improved. Contact Dale.Watt@tpl.org with feedback.

Spatial analysis and statistical summaries of the Protected Areas Database of the United States (PAD-US) provide land managers and decision makers with a general assessment of management intent for biodiversity protection, natural resource management, and recreation access across the nation. The PAD-US 4.0 Combined Fee, Designation, Easement feature class in the full geodatabase inventory (with Military Lands and Tribal Areas from the Proclamation and Other Planning Boundaries feature class) was modified to prioritize overlapping designations, avoiding massive overestimation in protected area statistics, and simplified by the following PAD-US attributes to support user needs for raster analysis data: Manager Type, Manager Name, Designation Type, GAP Status Code, Public Access, and State Name. The rasterization process prioritized overlapping designations previously identified (GAP_Prity field) in the Vector Analysis file (e.g. Wilderness within a National Forest) based upon their relative biodiversity conservation (e.g. GAP Status Code 1 over 2).The 30-meter Image (IMG) grid Raster Analysis Files area extents were defined by the Census state boundary file used to clip the Vector Analysis File, the data source for rasterization ("PADUS4_0VectorAnalysis_State_Clip_CENSUS2022") feature class from ("PADUS4_0VectorAnalysisFile_OtherExtents_ClipCENSUS2022.gdb"). Alaska (AK) and Hawaii (HI) raster data are separated from the contiguous U.S. (CONUS) to facilitate analyses at manageable scales. Note, the PAD-US inventory is now considered functionally complete with the vast majority of land protection types (with a legal protection mechanism) represented in some manner, while work continues to maintain updates, improve data quality, and integrate new data as it becomes available (see inventory completeness estimates at: http://www.protectedlands.net/data-stewards/ ). In addition, protection status represents a point-in-time and changes in status between versions of PAD-US may be attributed to improving the completeness and accuracy of the spatial data more than actual management actions or new acquisitions. USGS provides no legal warranty for the use of this data. While PAD-US is the official aggregation of protected areas ( https://ngda-portfolio-community-geoplatform.hub.arcgis.com/pages/portfolio ), agencies are the best source of their lands data.

Wetlands in California are protected by several federal and state laws, regulations, and policies. This layer was extracted from the broader vegetation raster from the CA Nature project which was recently enhanced to include a more comprehensive definition of wetland. This wetlands dataset is used as an exclusion as part of the biological planning priorities in the CEC 2023 Land-Use Screens.

This layer is featured in the CEC 2023 Land-Use Screens for Electric System Planning data viewer.

For more information about this layer and its use in electric system planning, please refer to the Land Use Screens Staff Report in the CEC Energy Planning Library.

Notice: this is not the latest Heat Island Severity image service. For 2023 data, visit https://tpl.maps.arcgis.com/home/item.html?id=db5bdb0f0c8c4b85b8270ec67448a0b6. This layer contains the relative heat severity for every pixel for every city in the United States, including Alaska, Hawaii, and Puerto Rico. This 30-meter raster was derived from Landsat 8 imagery band 10 (ground-level thermal sensor) from the summer of 2022, patched with data from 2021 where necessary.Federal statistics over a 30-year period show extreme heat is the leading cause of weather-related deaths in the United States. Extreme heat exacerbated by urban heat islands can lead to increased respiratory difficulties, heat exhaustion, and heat stroke. These heat impacts significantly affect the most vulnerable—children, the elderly, and those with preexisting conditions.The purpose of this layer is to show where certain areas of cities are hotter than the average temperature for that same city as a whole. Severity is measured on a scale of 1 to 5, with 1 being a relatively mild heat area (slightly above the mean for the city), and 5 being a severe heat area (significantly above the mean for the city). The absolute heat above mean values are classified into these 5 classes using the Jenks Natural Breaks classification method, which seeks to reduce the variance within classes and maximize the variance between classes. Knowing where areas of high heat are located can help a city government plan for mitigation strategies.This dataset represents a snapshot in time. It will be updated yearly, but is static between updates. It does not take into account changes in heat during a single day, for example, from building shadows moving. The thermal readings detected by the Landsat 8 sensor are surface-level, whether that surface is the ground or the top of a building. Although there is strong correlation between surface temperature and air temperature, they are not the same. We believe that this is useful at the national level, and for cities that don’t have the ability to conduct their own hyper local temperature survey. Where local data is available, it may be more accurate than this dataset. Dataset SummaryThis dataset was developed using proprietary Python code developed at The Trust for Public Land, running on the Descartes Labs platform through the Descartes Labs API for Python. The Descartes Labs platform allows for extremely fast retrieval and processing of imagery, which makes it possible to produce heat island data for all cities in the United States in a relatively short amount of time.What can you do with this layer?This layer has query, identify, and export image services available. Since it is served as an image service, it is not necessary to download the data; the service itself is data that can be used directly in any Esri geoprocessing tool that accepts raster data as input.In order to click on the image service and see the raw pixel values in a map viewer, you must be signed in to ArcGIS Online, then Enable Pop-Ups and Configure Pop-Ups.Using the Urban Heat Island (UHI) Image ServicesThe data is made available as an image service. There is a processing template applied that supplies the yellow-to-red or blue-to-red color ramp, but once this processing template is removed (you can do this in ArcGIS Pro or ArcGIS Desktop, or in QGIS), the actual data values come through the service and can be used directly in a geoprocessing tool (for example, to extract an area of interest). Following are instructions for doing this in Pro.In ArcGIS Pro, in a Map view, in the Catalog window, click on Portal. In the Portal window, click on the far-right icon representing Living Atlas. Search on the acronyms “tpl” and “uhi”. The results returned will be the UHI image services. Right click on a result and select “Add to current map” from the context menu. When the image service is added to the map, right-click on it in the map view, and select Properties. In the Properties window, select Processing Templates. On the drop-down menu at the top of the window, the default Processing Template is either a yellow-to-red ramp or a blue-to-red ramp. Click the drop-down, and select “None”, then “OK”. Now you will have the actual pixel values displayed in the map, and available to any geoprocessing tool that takes a raster as input. Below is a screenshot of ArcGIS Pro with a UHI image service loaded, color ramp removed, and symbology changed back to a yellow-to-red ramp (a classified renderer can also be used): A typical operation at this point is to clip out your area of interest. To do this, add your polygon shapefile or feature class to the map view, and use the Clip Raster tool to export your area of interest as a geoTIFF raster (file extension ".tif"). In the environments tab for the Clip Raster tool, click the dropdown for "Extent" and select "Same as Layer:", and select the name of your polygon. If you then need to convert the output raster to a polygon shapefile or feature class, run the Raster to Polygon tool, and select "Value" as the field.Other Sources of Heat Island InformationPlease see these websites for valuable information on heat islands and to learn about exciting new heat island research being led by scientists across the country:EPA’s Heat Island Resource CenterDr. Ladd Keith, University of ArizonaDr. Ben McMahan, University of Arizona Dr. Jeremy Hoffman, Science Museum of Virginia Dr. Hunter Jones, NOAA Daphne Lundi, Senior Policy Advisor, NYC Mayor's Office of Recovery and ResiliencyDisclaimer/FeedbackWith nearly 14,000 cities represented, checking each city's heat island raster for quality assurance would be prohibitively time-consuming, so The Trust for Public Land checked a statistically significant sample size for data quality. The sample passed all quality checks, with about 98.5% of the output cities error-free, but there could be instances where the user finds errors in the data. These errors will most likely take the form of a line of discontinuity where there is no city boundary; this type of error is caused by large temperature differences in two adjacent Landsat scenes, so the discontinuity occurs along scene boundaries (see figure below). The Trust for Public Land would appreciate feedback on these errors so that version 2 of the national UHI dataset can be improved. Contact Dale.Watt@tpl.org with feedback.

Crop Storage Final Location: Aggregated Production (Kenya - ~1km) consists of a 0.01 decimal degree grid produced under the scope of the Covid-19 sub-Saharan African Corridor project pilot case, using a Geographical Information Systems - Multicriteria Decision Analysis (GIS-MCDA) methodology for the identification and definition of mobile storage locations (movable warehouses). The Top Score Locations were computed using the Location Score raster and adding as criteria, access to finance (distance to bank agency) and linear distance from major roads. This output step filters the location score grid (aggregated production) into a top location score raster. Applying the following functions: • Buffering: o Banks - 10km (0.09 degree) buffer radius; o Major roads - 2km (0.18 degree) buffer radius; • Intersection - extracts the overlapping portions of Banks_Buffer and Roads_Buffer. • Dissolve - Takes the intersection vector layer and combines the features into a new feature, a single polygon; • Clip Raster by Mask Layer - Extracts the raster-production aggregate location score using the polygon obtained by the dissolved intersection of the banks and roads buffers; • The raster pixel cell in the final areas are classified into the final location map.

Digital surfaces and thicknesses of selected hydrogeologic units of the Floridan aquifer system were developed to define an updated hydrogeologic framework as part of the U.S. Geological Survey Groundwater Resources Program. This feature class contains a gridded surface depicting the top of the lower Avon Park permeable zone in feet relative to NGVD29.

Crop Storage Final Location: Aggregated Production, consists of a 0.01 decimal degree grid produced under the scope of the Covid-19 sub-Saharan African Corridor project pilot case, using a Geographical Information Systems - Multicriteria Decision Analysis (GIS-MCDA) methodology for the identification and definition of mobile storage locations (movable warehouses). The Top Score Locations were computed using the Location Score raster and adding as criteria, access to finance (distance to bank agency) and linear distance from major roads. This output step filters the location score grid (aggregated production) into a top location score raster. Applying the following functions: • Buffering: o Banks - 10km (0.09 degree) buffer radius; o Major roads - 2km (0.18 degree) buffer radius; • Intersection - extracts the overlapping portions of Banks_Buffer and Roads_Buffer. • Dissolve - Takes the intersection vector layer and combines the features into a new feature, a single polygon; • Clip Raster by Mask Layer - Extracts the raster-production aggregate location score using the polygon obtained by the dissolved intersection of the banks and roads buffers; • The raster pixel cell in the final areas are classified into the final location map.

This collection of files are part of a larger dataset uploaded in support of Low Temperature Geothermal Play Fairway Analysis for the Appalachian Basin (GPFA-AB, DOE Project DE-EE0006726). Phase 1 of the GPFA-AB project identified potential Geothermal Play Fairways within the Appalachian basin of Pennsylvania, West Virginia and New York. This was accomplished through analysis of 4 key criteria: thermal quality, natural reservoir productivity, risk of seismicity, and heat utilization. Each of these analyses represent a distinct project task, with the fifth task encompassing combination of the 4 risks factors. Supporting data for all five tasks has been uploaded into the Geothermal Data Repository node of the National Geothermal Data System (NGDS).

This submission comprises the data for Thermal Quality Analysis (project task 1) and includes all of the necessary shapefiles, rasters, datasets, code, and references to code repositories that were used to create the thermal resource and risk factor maps as part of the GPFA-AB project. The identified Geothermal Play Fairways are also provided with the larger dataset. Figures (.png) are provided as examples of the shapefiles and rasters. The regional standardized 1 square km grid used in the project is also provided as points (cell centers), polygons, and as a raster. Two ArcGIS toolboxes are available: 1) RegionalGridModels.tbx for creating resource and risk factor maps on the standardized grid, and 2) ThermalRiskFactorModels.tbx for use in making the thermal resource maps and cross sections. These toolboxes contain item description documentation for each model within the toolbox, and for the toolbox itself. This submission also contains three R scripts: 1) AddNewSeisFields.R to add seismic risk data to attribute tables of seismic risk, 2) StratifiedKrigingInterpolation.R for the interpolations used in the thermal resource analysis, and 3) LeaveOneOutCrossValidation.R for the cross validations used in the thermal interpolations.

Some file descriptions make reference to various 'memos'. These are contained within the final report submitted October 16, 2015.

Each zipped file in the submission contains an 'about' document describing the full Thermal Quality Analysis content available, along with key sources, authors, citation, use guidelines, and assumptions, with the specific file(s) contained within the .zip file highlighted.

UPDATE: Newer version of the Thermal Quality Analysis has been added here: https://gdr.openei.org/submissions/879 (Also linked below) Newer version of the Combined Risk Factor Analysis has been added here: https://gdr.openei.org/submissions/880 (Also linked below) This is one of sixteen associated .zip files relating to thermal resource interpolation results within the Thermal Quality Analysis task of the Low Temperature Geothermal Play Fairway Analysis for the Appalachian Basin. This file contains 6 images (.png) including predicted and associated error for surface heat flow, depth to 80 degrees C, depth to 100 degrees C, temperature at 1.5 km, temperature at 2.5 km and temperature at 3.5 km.

The sixteen files contain the results of the thermal resource interpolation as binary grid (raster) files, images (.png) of the rasters, and toolbox of ArcGIS Models used. Note that raster files ending in “pred” are the predicted mean for that resource, and files ending in “err” are the standard error of the predicted mean for that resource. Leave one out cross validation results are provided for each thermal resource.

Several models were built in order to process the well database with outliers removed. ArcGIS toolbox ThermalRiskFactorModels contains the ArcGIS processing tools used. First, the WellClipsToWormSections model was used to clip the wells to the worm sections (interpolation regions). Then, the 1 square km gridded regions (see series of 14 Worm Based Interpolation Boundaries .zip files) along with the wells in those regions were loaded into R using the rgdal package. Then, a stratified kriging algorithm implemented in the R gstat package was used to create rasters of the predicted mean and the standard error of the predicted mean. The code used to make these rasters is called StratifiedKrigingInterpolation.R Details about the interpolation, and exploratory data analysis on the well data is provided in 9_GPFA-AB_InterpolationThermalFieldEstimation.pdf (Smith, 2015), contained within the final report.

The output rasters from R are brought into ArcGIS for further spatial processing. First, the BufferedRasterToClippedRaster tool is used to clip the interpolations back to the Worm Sections. Then, the Mosaic tool in ArcGIS is used to merge all predicted mean rasters into a single raster, and all error rasters into a single raster for each thermal resource.

A leave one out cross validation was performed on each of the thermal resources. The code used to implement the cross validation is provided in the R script LeaveOneOutCrossValidation.R. The results of the cross validation are given for each thermal resource.

Other tools provided in this toolbox are useful for creating cross sections of the thermal resource. ExtractThermalPropertiesToCrossSection model extracts the predicted mean and the standard error of predicted mean to the attribute table of a line of cross section. The AddExtraInfoToCrossSection model is then used to add any other desired information, such as state and county boundaries, to the cross section attribute table. These two functions can be combined as a single function, as provided by the CrossSectionExtraction model.

The ISC2010_FRAGMENTATION polygon features represent gap areas in vegetation cover in the 40m Riparian Zone (defined as 40m from ISC2010_STREAMBED_WIDTH) . The gap crietria is defined as any area where vegetation cover is less than 20% for at least 10m x 10m. This data set is derived from source Fractional Cover raster data and has been clipped to a 40m buffer of Stream Bed. Small polygons (<10m2) resulting from the clip process have been excluded from the Fragmentation feature class.

River condition in Victoria is assessed every 5 years using the Index of Stream Condition (ISC). The Department of Environment and Primary Industries (DEPI) developed a methodology to assess the Physical Form and Riparian Vegetation components of the ISC using remote sensing data, specifically LIDAR and aerial photography.

A State Wide mapping project was undertaken in 2010-13 to accurately map the Physical Form and Riparian Vegetation metrics of the ISC . Other ISC metrics were not assessed in the project and were derived from other sources.

The Physical Form and Riparian Vegetation Metric products are a combination of mapped Vector and Raster data as well as Tabular Summary Statistics about the mapped features. In the context of the project, the term Metrics is used to refer to both the mapped features and the summary statistics.

Remote sensing data used includes 15cm true colour and infra-red aerial photography and four return multi-pulse LiDAR data. This source data was used to derive a variety of Raster data sets including Digital Terrain Models, Slope, Vegetation Height and Vegetation Cover. The Digital Terrain and Slope rasters were used to map Physical Form metrics including Stream Bed, Top of Bank and River Centre Lines while the Vegetation Height and Cover rasters were used to map the Riparian Vegetation metrics. The Project Report "Aerial Remote Sensing for Physical Channel Form and Riparian Vegetation Mapping" describes the remote sensing and mapping approach used to create this data set.

These files include a derived 50 meter spatial resolution Compound Topographic (or Wetness) Index (CTI or TWI) and Flow Accumulation (as represented by specific catchment area, SCA) calculated from a continuous, transboundary DEM developed across the Alaska perhumid coastal temperate rainforest (AKPCTR). The extent of this dataset includes all of the Alaska and Canada watersheds that discharge into southeast Alaska coastal waters, which covers essentially the northern half of the full PCTR.The transboundary DEM used to calculate the CTI can be here: linkFlow accumulation is represented by a grid of specific catchment area (SCA), which is the contributing area per unit contour length using the multiple flow direction D-infinity approach. Unit contour length is equal to the DEM resolution of 50 meters.There are two versions of both the CTI and SCA provided:CTI_AKPCTR_NoFlowMask.tif and SCA_AKPCTR_NoFlowMask.tif This is the resulting CTI file when no masks are applied to the study area before the CTI procedure is run. After the CTI procedure is run, all glacial areas are masked out, as CTI is not meaningful over glaciers.CTI_AKPCTR_FlowMask.tif and SCA_AKPCTR_FlowMask.tif This is the resulting CTI file where prior to running the CTI procedure, we apply a mask across all active glaciers and all downslope cells receiving flow from glacially affected cells. These masked cells are excluded from the CTI procedure. This results in fewer CTI values on the landscape. This dataset is provided to identify cells in the direct downslope path of glaciers because CTI values for cells receiving upslope accumulation from glaciers may not be reliable due to uncertainties in surface water flowpaths in glaciated areas.Processing Steps: 1. Using the pitremove function from TAUDEM, filled sinks in the DEM following the method of Planchon and Darboux (2001). 2. (For *_FlowMask versions only) Mask out glacial areas and associated downslope cells. Prior to calculating flow direction and flow accumulation, change cells in the filled DEM to NoData where glaciers exist, using the Randolph Glacier Inventory version 5.0 raster dataset to identify the presence of glaciers. Cells downslope from the masked glaciers will then also be identified as NoData in subsequent processing and the final CTI raster. Compute D-infinity slopes and flow direction and flow accumulation (as specific catchment area) rasters from the filled DEM using TAUDEM D-infinity method. As part of the D-infinity routine, TauDEM uses the method of Garbrecht & Martz (1997) to resolve flats.3. The D-infinity slope raster was modified in the following way: if slope = 0, then change slope value to 0.0001; otherwise leave slope value unchanged. This was done to avoid dividing by zero when calculating CTI.4. Compute Compound Topographic Index (CTI) using these steps: Dsca = D-infinity specific catchment area raster Dslp = D-infinity slope raster Then, CTI = Ln(Dsca/Dslp)5. Clip the CTI raster to the Alaska perhumid coastal temperate rainforest (AKPCTR) watershed boundary.

Notice: this is not the latest Heat Island Anomalies image service. For 2023 data visit https://tpl.maps.arcgis.com/home/item.html?id=e89a556263e04cb9b0b4638253ca8d10.This layer contains the relative degrees Fahrenheit difference between any given pixel and the mean heat value for the city in which it is located, for every city in the contiguous United States, Alaska, Hawaii, and Puerto Rico. This 30-meter raster was derived from Landsat 8 imagery band 10 (ground-level thermal sensor) from the summer of 2022, with patching from summer of 2021 where necessary.Federal statistics over a 30-year period show extreme heat is the leading cause of weather-related deaths in the United States. Extreme heat exacerbated by urban heat islands can lead to increased respiratory difficulties, heat exhaustion, and heat stroke. These heat impacts significantly affect the most vulnerable—children, the elderly, and those with preexisting conditions.The purpose of this layer is to show where certain areas of cities are hotter or cooler than the average temperature for that same city as a whole. This dataset represents a snapshot in time. It will be updated yearly, but is static between updates. It does not take into account changes in heat during a single day, for example, from building shadows moving. The thermal readings detected by the Landsat 8 sensor are surface-level, whether that surface is the ground or the top of a building. Although there is strong correlation between surface temperature and air temperature, they are not the same. We believe that this is useful at the national level, and for cities that don’t have the ability to conduct their own hyper local temperature survey. Where local data is available, it may be more accurate than this dataset. Dataset SummaryThis dataset was developed using proprietary Python code developed at The Trust for Public Land, running on the Descartes Labs platform through the Descartes Labs API for Python. The Descartes Labs platform allows for extremely fast retrieval and processing of imagery, which makes it possible to produce heat island data for all cities in the United States in a relatively short amount of time.In order to click on the image service and see the raw pixel values in a map viewer, you must be signed in to ArcGIS Online, then Enable Pop-Ups and Configure Pop-Ups.Using the Urban Heat Island (UHI) Image ServicesThe data is made available as an image service. There is a processing template applied that supplies the yellow-to-red or blue-to-red color ramp, but once this processing template is removed (you can do this in ArcGIS Pro or ArcGIS Desktop, or in QGIS), the actual data values come through the service and can be used directly in a geoprocessing tool (for example, to extract an area of interest). Following are instructions for doing this in Pro.In ArcGIS Pro, in a Map view, in the Catalog window, click on Portal. In the Portal window, click on the far-right icon representing Living Atlas. Search on the acronyms “tpl” and “uhi”. The results returned will be the UHI image services. Right click on a result and select “Add to current map” from the context menu. When the image service is added to the map, right-click on it in the map view, and select Properties. In the Properties window, select Processing Templates. On the drop-down menu at the top of the window, the default Processing Template is either a yellow-to-red ramp or a blue-to-red ramp. Click the drop-down, and select “None”, then “OK”. Now you will have the actual pixel values displayed in the map, and available to any geoprocessing tool that takes a raster as input. Below is a screenshot of ArcGIS Pro with a UHI image service loaded, color ramp removed, and symbology changed back to a yellow-to-red ramp (a classified renderer can also be used): A typical operation at this point is to clip out your area of interest. To do this, add your polygon shapefile or feature class to the map view, and use the Clip Raster tool to export your area of interest as a geoTIFF raster (file extension ".tif"). In the environments tab for the Clip Raster tool, click the dropdown for "Extent" and select "Same as Layer:", and select the name of your polygon. If you then need to convert the output raster to a polygon shapefile or feature class, run the Raster to Polygon tool, and select "Value" as the field.Other Sources of Heat Island InformationPlease see these websites for valuable information on heat islands and to learn about exciting new heat island research being led by scientists across the country:EPA’s Heat Island Resource CenterDr. Ladd Keith, University of ArizonaDr. Ben McMahan, University of Arizona Dr. Jeremy Hoffman, Science Museum of Virginia Dr. Hunter Jones, NOAA Daphne Lundi, Senior Policy Advisor, NYC Mayor's Office of Recovery and ResiliencyDisclaimer/FeedbackWith nearly 14,000 cities represented, checking each city's heat island raster for quality assurance would be prohibitively time-consuming, so The Trust for Public Land checked a statistically significant sample size for data quality. The sample passed all quality checks, with about 98.5% of the output cities error-free, but there could be instances where the user finds errors in the data. These errors will most likely take the form of a line of discontinuity where there is no city boundary; this type of error is caused by large temperature differences in two adjacent Landsat scenes, so the discontinuity occurs along scene boundaries (see figure below). The Trust for Public Land would appreciate feedback on these errors so that version 2 of the national UHI dataset can be improved. Contact Dale.Watt@tpl.org with feedback.

Carte de l'occupation du sol en 25 classes issue de capteurs multiples, disponible sur la période 1992 à 2015 (CLIP_OMVS_ESACCI-LC-L4-LCCS-Map-300m-P1Y-2015-v.2.0.7.tif)

Ce modèle numérique d'altitude (MNT) a été créé à partir du jeu de données de terrain principal (MTD) de Hakai au moyen de l'outil « MNT to raster » dans ArcGIS for Desktop d'ESRI à l'aide d'une méthode d'échantillonnage Natural Neighbour. Le DEM a été créé en mode natif à une résolution de 3 m. Ce DEM a été fixé à une zone tampon à 10 m du rivage. Une combinaison de différentes altitudes autour de l'île a été utilisée pour créer le rivage.

Le MNT qui en résulte est un modèle d'élévation hydroaplati en terre nue et donc considéré comme « topographiquement complet ». Chaque pixel représente l'altitude en mètres au-dessus du niveau moyen de la mer de la terre nue à cet endroit. Le système de référence vertical est le « Système de référence géodésique vertical canadien 1928 » (CGVD28).

Hakai a produit des DEM à différentes résolutions de manière native directement à partir du MTD des données LiDAR. Pour vos recherches, veuillez utiliser le produit de résolution approprié parmi ceux produits par Hakai. Afin de maintenir l'homogénéité, il n'est pas recommandé de procéder à un suréchantillonnage ou à une mise à l'échelle supérieure à partir de produits de résolution supérieure car cela pourrait introduire et propager des erreurs de différentes grandeurs dans les analyses en cours ; veuillez utiliser des produits déjà disponibles, et si vous avez besoin d'une résolution non disponible, contactez data@hakai.org afin d'obtenir un DEM produit directement à partir du MTD.

Les DEM topographiquement complets suivants ont été produits en mode natif à partir du DTM par Hakai :

MNE topographiquement complète de 3 m. Ce produit a été utilisé pour produire les ensembles de données hydrologiques de Hakai (cours d'eau et bassins versants) DEM Topographiquement complet de 20 m. Compatible avec les mesures du couvert végétal de Hakai et les rasters associés. MNT topographiquement complet de 25 m. Compatible avec les produits de données TRIM BCGov. DEM Topographiquement complet de 30 m. Compatible avec les produits STRM.

Création du jeu de données de terrain principal

Nuages de points LiDAR issus de missions effectuées en 2012 et 2014 au-dessus de l'île Calvert où ils ont été chargés (XYZ uniquement) dans une classe d'entités ponctuelles d'une géodatabase ESRI.

Seul le sol (classe 2) renvoie l'endroit où il est chargé dans la géodatabase.

Le « jeu de données de MNT » ESRI a été créé dans la même géodatabase à l'aide des points LiDAR en tant que points de masse intégrés.

Les lacs et les étangs TEM Plus avec des valeurs d'altitude moyennes au-dessus des miroirs des plans d'eau ont été utilisés comme lignes de rupture de remplacement dur pour obtenir un hydroaplatissement.

La géométrie d'emprise minimale de toutes les étendues de fichiers LAS contigus a été utilisée comme masque de découpe souple lors de la création du jeu de données de MNT en tant que limite de projet.

Le système de coordonnées horizontales et le datum utilisés pour le jeu de données de MNT sont : UTM Zone 9 NAD1983 ; le système de référence vertical a été défini sur CGVD28. Les deux systèmes de référence correspondent au système de référence natif des nuages de points LiDAR.

L'espacement minimal des points défini pendant la création du jeu de données de MNT a été défini sur 1.

Cet ensemble de données fournit les limites des bassins versants dérivés du LiDAR pour toutes les îles Calvert et Hecate, en Colombie-Britannique. Les bassins versants ont été délimités à partir d'un modèle altimétrique numérique de 3 m. Pour chaque polygone de bassin versant, le jeu de données comprend un identificateur unique et des statistiques sommaires simples pour décrire la topographie et l'hydrologie. Polygones de bassin versant Cet ensemble de données a été produit à partir des résultats de la modélisation hydrologique « traditionnelle » menée à l'aide du MNT de terre nue complet topographiquement complet basé sur lidar de 2012 + 2014 avec une zone tampon de 10 m autour de la côte pour s'assurer que tous les bassins versants modélisés atteignent l'océan. Les bassins versants ont été délimités à l'aide de points d'coulée créés à l'intersection des cours d'eau modélisés et du littoral. Après la délimitation du bassin versant, ceux-ci ont été coupés sur le rivage de l'île.

Carte de l'occupation du sol en 10 classes issue de l'instrument imageur de Sentinel 2 (Clip_ESACCI-LC-L4-LC10-Map-20m-P1Y-2016-v.1.0.tif)