Geographic Information System (GIS) analyses are an essential part of natural resource management and research. Calculating and summarizing data within intersecting GIS layers is common practice for analysts and researchers. However, the various tools and steps required to complete this process are slow and tedious, requiring many tools iterating over hundreds, or even thousands of datasets. USGS scientists will combine a series of ArcGIS geoprocessing capabilities with custom scripts to create tools that will calculate, summarize, and organize large amounts of data that can span many temporal and spatial scales with minimal user input. The tools work with polygons, lines, points, and rasters to calculate relevant summary data and combine them into a single output table that can be easily incorporated into statistical analyses. These tools are useful for anyone interested in using an automated script to quickly compile summary information within all areas of interest in a GIS dataset.

Toolbox Use
License
Creative Commons-PDDC
Recommended Citation
Welty JL, Jeffries MI, Arkle RS, Pilliod DS, Kemp SK. 2021. GIS Clipping and Summarization Toolbox: U.S. Geological Survey Software Release. https://doi.org/10.5066/P99X8558

Analysis Image Service generated from Extract Raster Data

Lesson 1. An Introduction to working with multispectral satellite data in ArcGIS Pro In which we learn: • How to unpack tar and gz files from USGS EROS • The basic map interface in ArcGIS • How to add image files • What each individual band of Landsat spectral data looks like • The difference between: o Analysis-ready data: surface reflectance and surface temperature o Landsat Collection 1 Level 3 data: burned area and dynamic surface water o Sentinel2data o ISRO AWiFS and LISS-3 data Lesson 2. Basic image preprocessing In which we learn: • How to composite using the composite band tool • How to represent composite images • All about band combinations • How to composite using raster functions • How to subset data into a rectangle • How to clip to a polygon Lesson 3. Working with mosaic datasets In which we learn: o How to prepare an empty mosaic dataset o How to add images to a mosaic dataset o How to change symbology in a mosaic dataset o How to add a time attribute o How to add a time dimension to the mosaic dataset o How to view time series data in a mosaic dataset Lesson 4. Working with and creating derived datasets In which we learn: • How to visualize Landsat ARD surface temperature • How to calculate F° from K° using ARD surface temperature • How to generate and apply .lyrx files • How to calculate an NDVI raster using ISRO LISS-3 data • How to visualize burned areas using Landsat Level 3 data • How to visualize dynamic surface water extent using Landsat Level 3 data

Analysis Image Service generated from Extract Raster Data

Notice: this is not 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. Heat Severity is a reclassified version of Heat Anomalies raster which is also published on this site. This data is generated from 30-meter 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.

Notice: this is not 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 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.

Dataset description-br /- This dataset is a recalculation of the Copernicus 2015 high resolution layer (HRL) of imperviousness density data (IMD) at different spatial/territorial scales for the case studies of Barcelona and Milan. The selected spatial/territorial scales are the following: * a) Barcelona city boundaries * b) Barcelona metropolitan area, Àrea Metropolitana de Barcelona (AMB) * c) Barcelona greater city (Urban Atlas) * d) Barcelona functional urban area (Urban Atlas) * e) Milan city boundaries * f) Milan metropolitan area, Piano Intercomunale Milanese (PIM) * g) Milan greater city (Urban Atlas) * h) Milan functional urban area (Urban Atlas)-br /- In each of the spatial/territorial scales listed above, the number of 20x20mt cells corresponding to each of the 101 values of imperviousness (0-100% soil sealing: 0% means fully non-sealed area; 100% means fully sealed area) is provided, as well as the converted measure into squared kilometres (km2). -br /- -br /- -br /- Dataset composition-br /- The dataset is provided in .csv format and is composed of: -br /- _IMD15_BCN_MI_Sources.csv_: Information on data sources -br /- _IMD15_BCN.csv_: This file refers to the 2015 high resolution layer of imperviousness density (IMD) for the selected territorial/spatial scales in Barcelona: * a) Barcelona city boundaries (label: bcn_city) * b) Barcelona metropolitan area, Àrea metropolitana de Barcelona (AMB) (label: bcn_amb) * c) Barcelona greater city (Urban Atlas) (label: bcn_grc) * d) Barcelona functional urban area (Urban Atlas) (label: bcn_fua)-br /- _IMD15_MI.csv_: This file refers to the 2015 high resolution layer of imperviousness density (IMD) for the selected territorial/spatial scales in Milan: * e) Milan city boundaries (label: mi_city) * f) Milan metropolitan area, Piano intercomunale milanese (PIM) (label: mi_pim) * g) Milan greater city (Urban Atlas) (label: mi_grc) * h) Milan functional urban area (Urban Atlas) (label: mi_fua)-br /- _IMD15_BCN_MI.mpk_: the shareable project in Esri ArcGIS format including the HRL IMD data in raster format for each of the territorial boundaries as specified in letter a)-h). -br /- Regarding the territorial scale as per letter f), the list of municipalities included in the Milan metropolitan area in 2016 was provided to me in 2016 from a person working at the PIM. -br /- In the IMD15_BCN.csv and IMD15_MI.csv, the following columns are included: * Level: the territorial level as defined above (a)-d) for Barcelona and e)-h) for Milan); * Value: the 101 values of imperviousness density expressed as a percentage of soil sealing (0-100%: 0% means fully non-sealed area; 100% means fully sealed area); * Count: the number of 20x20mt cells corresponding to a certain percentage of soil sealing or imperviousness; * Km2: the conversion of the 20x20mt cells into squared kilometres (km2) to facilitate the use of the dataset.-br /- -br /- -br /- Further information on the Dataset-br /- This dataset is the result of a combination between different databases of different types and that have been downloaded from different sources. Below, I describe the main steps in data management that resulted in the production of the dataset in an Esri ArcGIS (ArcMap, Version 10.7) project.-br /- 1. The high resolution layer (HRL) of the imperviousness density data (IMD) for 2015 has been downloaded from the official website of Copernicus. At the time of producing the dataset (April/May 2021), the 2018 version of the IMD HRL database was not yet validated, so the 2015 version was chosen instead. The type of this dataset is raster. 2. For both Barcelona and Milan, shapefiles of their administrative boundaries have been downloaded from official sources, i.e. the ISTAT (Italian National Statistical Institute) and the ICGC (Catalan Institute for Cartography and Geology). These files have been reprojected to match the IMD HRL projection, i.e. ETRS 1989 LAEA. 3. Urban Atlas (UA) boundaries for the Greater Cities (GRC) and Functional Urban Areas (FUA) of Barcelona and Milan have been checked and reconstructed in Esri ArcGIS from the administrative boundaries files by using a Eurostat correspondence table. This is because at the time of the dataset creation (April/May 2021), the 2018 Urban Atlas shapefiles for these two cities were not fully updated or validated on the Copernicus Urban Atlas website. Therefore, I had to re-create the GRC and FUA boundaries by using the Eurostat correspondence table as an alternative (but still official) data source. The use of the Eurostat correspondence table with the codes and names of municipalities was also useful to detect discrepancies, basically stemming from changes in municipality names and codes and that created inconsistent spatial features. When detected, these discrepancies have been checked with the ISTAT and ICGC offices in charge of producing Urban Atlas data before the final GRC and FUA boundaries were defined.-br /- Steps 2) and 3) were the most time consuming, because they required other tools to be used in Esri ArcGIS, like spatial joins and geoprocessing tools for shapefiles (in particular dissolve and area re-calculator in editing sessions) for each of the spatial/territorial scales as indicated in letters a)-h). -br /- Once the databases for both Barcelona and Milan as described in points 2) and 3) were ready (uploaded in Esri ArcGIS, reprojected and their correctness checked), they have been ‘crossed’ (i.e. clipped) with the IMD HRL as described in point 1) and a specific raster for each territorial level has been calculated. The procedure in Esri ArcGIS was the following: * Clipping: Arctoolbox - Data management tools - Raster - Raster Processing - Clip. The ‘input’ file is the HRL IMD raster file as described in point 1) and the ‘output’ file is each of the spatial/territorial files. The option "Use Input Features for Clipping Geometry (optional)” was selected for each of the clipping. * Delete and create raster attribute table: Once the clipping has been done, the raster has to be recalculated first through Arctoolbox - Data management tools - Raster - Raster properties - Delete Raster Attribute Table and then through Arctoolbox - Data management tools - Raster - Raster properties - Build Raster Attribute Table; the "overwrite" option has been selected. -br /- -br /- Other tools used for the raster files in Esri ArcGIS have been the spatial analyst tools (in particular, Zonal - Zonal Statistics). As an additional check, the colour scheme of each of the newly created raster for each of the spatial/territorial attributes as per letters a)-h) above has been changed to check the consistency of its overlay with the original HRL IMD file. However, a perfect match between the shapefiles as per letters a)-h) and the raster files could not be achieved since the raster files are composed of 20x20mt cells.-br /- The newly created attribute tables of each of the raster files have been exported and saved as .txt files. These .txt files have then been copied in the excel corresponding to the final published dataset.

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).

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.

The Grid Garage Toolbox is designed to help you undertake the Geographic Information System (GIS) tasks required to process GIS data (geodata) into a standard, spatially aligned format. This format is required by most, grid or raster, spatial modelling tools such as the Multi-criteria Analysis Shell for Spatial Decision Support (MCAS-S). Grid Garage contains 36 tools designed to save you time by batch processing repetitive GIS tasks as well diagnosing problems with data and capturing a record of processing step and any errors encountered.

Grid Garage provides tools that function using a list based approach to batch processing where both inputs and outputs are specified in tables to enable selective batch processing and detailed result reporting. In many cases the tools simply extend the functionality of standard ArcGIS tools, providing some or all of the inputs required by these tools via the input table to enable batch processing on a 'per item' basis. This approach differs slightly from normal batch processing in ArcGIS, instead of manually selecting single items or a folder on which to apply a tool or model you provide a table listing target datasets. In summary the Grid Garage allows you to:

• List, describe and manage very large volumes of geodata.
• Batch process repetitive GIS tasks such as managing (renaming, describing etc.) or processing (clipping, resampling, reprojecting etc.) many geodata inputs such as time-series geodata derived from satellite imagery or climate models.
• Record any errors when batch processing and diagnose errors by interrogating the input geodata that failed.
• Develop your own models in ArcGIS ModelBuilder that allow you to automate any GIS workflow utilising one or more of the Grid Garage tools that can process an unlimited number of inputs.
• Automate the process of generating MCAS-S TIP metadata files for any number of input raster datasets.

The Grid Garage is intended for use by anyone with an understanding of GIS principles and an intermediate to advanced level of GIS skills. Using the Grid Garage tools in ArcGIS ModelBuilder requires skills in the use of the ArcGIS ModelBuilder tool.

Download Instructions: Create a new folder on your computer or network and then download and unzip the zip file from the GitHub Release page for each of the following items in the 'Data and Resources' section below. There is a folder in each zip file that contains all the files. See the Grid Garage User Guide for instructions on how to install and use the Grid Garage Toolbox with the sample data provided.

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.

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.

ArcGIS Pro/QGIS to modify layers R for scripts

TIFF raster image format with TFW files. It represents the entire municipal territory divided into 29 sheets. The information content is on a scale of 1:5000. The technical paper is made up of: 1) geometric elements. 2) constituent elements of the anthropic landscape such as: buildings, technical artifacts, roads, railways, canals, trees and rows, etc; 3) constitutive elements of the natural landscape such as: hydrography, vegetation, etc.; 4) administrative limits; 5) toponymy. WMS 1.3.0 service available at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1946_TIF_CTC_GB/ImageServer/WMSServer and download service for portions of cartography at http://geodata.sitmilano .opendata.arcgis.com/ - Raster image TIFF file format with TFW file. It represents the entire territory of the Municipality of Milan divided into 29 tiles. The scale is 1:5000. The Municipal Base Map is composed by: 1) geometric elements. 2) elements of the anthropic landscape such as buildings, technical structures, roads, railroads, canals, trees and rows of plants, etc.; 3) elements of the natural landscape such as hydrography, vegetation, etc.; 4) administrative boundaries; 5) toponyms. Available WMS service 1.3.0 at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1946_TIF_CTC_GB/ImageServer/WMSServer and Clip and ship service of map areas at http://geodata.sitmilano.opendata.arcgis .com/

Terminal lakes are lakes with no hydrologic surface outflows and with losses of water occurring only through surface evaporation and groundwater discharge. We quantified the extent of the littoral zones (areas where 1% or more of surface irradiation reaches the lake bottom) and open water zones (areas where less than 1% of surface irradiation reaches the lake bottom) in 18 terminal lakes. Additionally, we quantified habitat usage and diets of the fish species inhabiting these lakes. This dataset contains includes seven lakes from North America (Atitlan, Crater, Eagle, Mann, Pyramid, Summit, Walker), one from South America (Titicaca), five from Eurasia (Caspian, Issyk-Kul, Neusiedl, Qinghai, Van), and five from Africa (Abijatta, Manyara, Nakuru, Shala, Turkana). Methods Measurements of the surface areas of the littoral and open water zones were performed using ArcGIS Pro Version 2.9. First, we generated year-specific digital elevation models (DEMs) of the lake’s bathymetry by a) using existing bathymetry raster data or b) by digitizing published depth contours of the lake’s bathymetry and interpolating a bathymetry raster using a natural neighbor interpolation. For several lakes that showed significant changes in lake level and where data regarding lake level change were available, we were able to produce a second year closer to the present by using the Raster Calculator function in ArcGIS Pro and then clipping the bathymetry raster to the lower lake level. This was possible for 5 of the 18 lakes (Mann Lake, Eagle Lake, Lake Abijatta, Walker Lake, and Lake Turkana), allowing us to map changes in the littoral zone size between the two years. For the lakes containing two years of data, we used only the most recent year in all subsequent analyses. We defined the portions of the littoral zone of the lake as the portions where the intensity of photosynthetically active radiation (PAR) reaching the lake bottom is 1% or greater relative to the intensity at the surface. For lakes where 1% PAR depth was not published, we calculated 1% PAR depth from published light profiles using the Lambert-Beer Law: 0.01 = e-u*z where µ is the light attenuation coefficient (meters-1) and z is 1% PAR depth (meters). For lakes where neither 1% PAR depth nor light profiles were published, we approximated the 1% PAR depth by multiplying the Secchi depth of the lake by a coefficient of 2.5. We sought the most recently collected Secchi depth to make these calculations. We then used the Raster Calculator function in ArcGIS PRO 2.9 to determine the portions of the lake where depth was less than or greater than the 1% PAR depth to map the open water and littoral zones, respectively. Fish species inventories and information regarding each species’ habitat and diet was compiled from 1) published peer-reviewed primary literature, 2) non-peer-reviewed literature (books, reports by government agencies or private firms), 3) online databases (i.e., FishBase (https://www.fishbase.de/home.htm), California Fish Website (www.calfish.ucdavis.edu)), and/or 4) experts studying the ecology of the species or lake ecosystem. We employed a conservative view regarding species taxonomy (i.e., ‘lumping’ rather than ‘splitting’). We classified species’ habitats with respect to three categories: 1) littoral zone (occurring in parts of the lake where 1% or more of the surface radiation reaches the lake bottom), 2) open water zone (occurring in parts of the lake where less than 1% of the surface radiation reaches the lake bottom), and 3) littoral & open water zone (occurring in both lake zones). These habitat classifications were based on adult habitat use only, and habitat use during larval and juvenile stages was not considered. We classified diets with respect to seven categories: 1) plankton only, 2) periphyton only, 3) periphyton and macroinvertebrates, 4) periphyton, macroinvertebrates, and plankton, 5) periphyton, macroinvertebrates, and fish, 6) fish OR fish and plankton, and 7) fish, plankton, periphyton, and macroinvertebrates.

JPEG raster image format with JGW files. It represents the entire municipal territory divided into 341 sheets. The information content is on a scale of 1:1,000. The Technical Map is made up of: 1) elements and entities of a geometric type: cartographic grid, quoted points. The altimetry, expressed in metres, relating to both the ground and the buildings, refers to the average sea level (Mareograph of Genoa 1942). The equidistance between the contour lines is 1 metre; for auxiliary curves (partly) it is 0.5 metres; 2) constituent elements of the anthropic landscape such as: buildings destination indication for public buildings, technical artifacts, roads, railways, canals, trees and rows, etc; 3) constitutive elements of the natural landscape such as: hydrography, vegetation, etc.; 4) administrative limits; 5) toponymy. WMS service available at http://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_2000_JPG_CTC_GB/ImageServer/WMSServer and download service for portions of cartography at http://geodata.sitmilano.opendata. arcgis.com/ - Raster image JPEG file format with JGW file. It represents the entire territory of the Municipality of Milan divided into 341 tiles. The scale is 1:1000. The Municipal Base Map is composed by: 1) geometric elements and entities: metric grid, elevation points. The elevation points, in meters, are both of ground and buildings and are referred to the average sea level (tide gauge of Genova 1942). The interval for contour lines is 1 meter and for auxiliary lines (dashed) is 0.5 meters; 2) elements of the anthropic landscape such as buildings (with building function for public buildings), technical structures, roads, railroads, canals, trees and rows of plants, etc.; 3) elements of the natural landscape such as hydrography, vegetation, etc.; 4) administrative boundaries; 5) toponyms. Available WMS service at http://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_2000_JPG_CTC_GB/ImageServer/WMSServer and Clip and ship service of map areas at http://geodata.sitmilano.opendata.arcgis.com/

JPEG raster image format with JGW files. It represents the entire municipal territory divided into 343 sheets. The information content is on a scale of 1:1,000. The Technical Map is made up of: 1) elements and entities of a geometric type: cartographic grid, quoted points. The altimetry, expressed in metres, relating to both the ground and the buildings, refers to the average sea level (Mareograph of Genoa 1942). The equidistance between the contour lines is 1 metre; for auxiliary curves (partly) it is 0.5 metres; 2) constituent elements of the anthropic landscape such as: buildings destination indication for public buildings, technical artifacts, roads, railways, canals, trees and rows, etc; 3) constitutive elements of the natural landscape such as: hydrography, vegetation, etc.; 4) administrative limits; 5) toponymy. Download of the entire cartography available through the dedicated service https://geoportale.comune.milano.it/ATOM/SIT/CTC2012/CTC2012_Service.xml (compressed file .zip - 3.19 Gb). Also available WMS service at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_2012_JPG_CTC_1000_RDN2008_TM32/ImageServer/WMSServer and download service for portions of cartography at http://geodata.sitmilano.opendata .arcgis.com/ - Raster image JPEG file format with JGW file. It represents the entire territory of the Municipality of Milan divided into 343 tiles. The scale is 1:1000. The Municipal Base Map is composed by: 1) geometric elements and entities: metric grid, elevation points. The elevation points, in meters, are both of ground and buildings and are referred to the average sea level (tide gauge of Genova 1942). The interval for contour lines is 1 meter and for auxiliary lines (dashed) is 0.5 meters; 2) elements of the anthropic landscape such as buildings (with building function for public buildings), technical structures, roads, railroads, canals, trees and rows of plants, etc.; 3) elements of the natural landscape such as hydrography, vegetation, etc.; 4) administrative boundaries; 5) toponyms. Download of the entire map available at https://geoportale.comune.milano.it/ATOM/SIT/CTC2012/CTC2012_Service.xml (.zip file - 3,19 Gb). Also available WMS service at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_2012_JPG_CTC_1000_RDN2008_TM32/ImageServer/WMSServer and Clip and ship service of map areas at http://geodata.sitmilano.opendata.arcgis.com /

TIFF raster image format with TFW files. It represents the entire municipal territory divided into 116 sheets. The information content is on a scale of 1:2000. The Technical Map is made up of: 1) elements and entities of a geometric type: cartographic grid, quoted points. 2) constituent elements of the anthropic landscape such as: buildings not covered with indication of destination for public buildings, technical artifacts, roads, railways, canals, trees and rows, etc. 3) constitutive elements of the natural landscape such as: hydrography, vegetation, etc.; 4) administrative limits; 5) toponymy. WMS service available at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1972_PNG_CTC_GB/ImageServer/WMSServer and download service for portions of cartography at http://geodata.sitmilano.opendata. arcgis.com/ - Raster image TIFF file format with TFW file. It represents the entire territory of the Municipality of Milan divided into 116 tiles. The scale is 1:2000. The Municipal Base Map is composed by: 1) geometric elements and entities: metric grid, elevation points. 2) elements of the anthropic landscape such as buildings (drawn without graphic pattern; function is specified for public buildings), technical structures, roads, railroads, canals, trees and rows of plants, etc.; 3) elements of the natural landscape such as hydrography, vegetation, etc. 4) administrative boundaries; 5) toponyms. Available WMS service at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1972_PNG_CTC_GB/ImageServer/WMSServer and Clip and ship service of map areas at http://geodata.sitmilano.opendata.arcgis.com/

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.

TIFF raster image format with TFW files. It represents the entire municipal territory divided into 57 sheets. The information content is on a scale of 1:5000. The technical paper is made up of: 1) geometric elements. 2) constituent elements of the anthropic landscape such as: buildings, technical artifacts, roads, railways, canals, trees and rows, etc; 3) constitutive elements of the natural landscape such as: hydrography, vegetation, etc.; 4) administrative limits; 5) toponymy. WMS 1.3.0 service available at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1930_TIF_CTC_GB/ImageServer/WMSServer and download service for portions of cartography at http://geodata.sitmilano .opendata.arcgis.com/ PLEASE NOTE: The map derives from the georeferencing of paper map sheets acquired via scanner. The paper support available at the time for the acquisition was found to be missing sheets 8,14,37 which are therefore not present in the digital format. - Raster image TIFF file format with TFW file. It represents the entire territory of the Municipality of Milan divided into 57 tiles. The scale is 1:5000. The Municipal Base Map is composed by: 1) geometric elements. 2) elements of the anthropic landscape such as buildings, technical structures, roads, railroads, canals, trees and rows of plants, etc.; 3) elements of the natural landscape such as hydrography, vegetation, etc.; 4) administrative boundaries; 5) toponyms. Available WMS service 1.3.0 at https://geoportale.comune.milano.it/arcgis/services/Cartografie_Raster/Milano_1930_TIF_CTC_GB/ImageServer/WMSServer and Clip and ship service of map areas at http://geodata.sitmilano.opendata.arcgis .com/ PLEASE NOTE: Digital map as result of scan process over a paper map. Original map had a lack of tiles number 8,14,37 during scan process, so that mentioned tiles are not available in digital copy.