Calculates zonal statistics on polygons from many categorical rasters for multiple attributes

The National Insect and Disease Risk map identifies areas with risk of significant tree mortality due to insects and plant diseases. The layer identifies lands in three classes: areas with risk of tree mortality from insects and disease between 2013 and 2027, areas with lower tree mortality risk, and areas that were formerly at risk but are no longer at risk due to disturbance (human or natural) between 2012 and 2018. Areas with risk of tree mortality are defined as places where at least 25% of standing live basal area greater than one inch in diameter will die over a 15-year time frame (2013 to 2027) due to insects and diseases.The National Insect and Disease Risk map, produced by the US Forest Service FHAAST, is part of a nationwide strategic assessment of potential hazard for tree mortality due to major forest insects and diseases. Dataset Summary Phenomenon Mapped: Risk of tree mortality due to insects and diseaseUnits: MetersCell Size: 30 meters in Hawaii and 240 meters in Alaska and the Contiguous USSource Type: DiscretePixel Type: 2-bit unsigned integerData Coordinate System: NAD 1983 Albers (Contiguous US), WGS 1984 Albers (Alaska), Hawaii Albers (Hawaii)Mosaic Projection: North America Albers Equal Area ConicExtent: Alaska, Hawaii, and the Contiguous United States Source: National Insect Disease Risk MapPublication Date: 2018ArcGIS Server URL: https://landscape11.arcgis.com/arcgis/This layer was created from the 2018 version of the National Insect Disease Risk Map.What can you do with this Layer? This layer is suitable for both visualization and analysis across the ArcGIS system. This layer can be combined with your data and other layers from the ArcGIS Living Atlas of the World in ArcGIS Online and ArcGIS Pro to create powerful web maps that can be used alone or in a story map or other application.Because this layer is part of the ArcGIS Living Atlas of the World it is easy to add to your map:In ArcGIS Online you can add this layer to a map by selecting Add then Browse Living Atlas Layers. A window will open. Type "insects and disease" in the search box and browse to the layer. Select the layer then click Add to Map.In ArcGIS Pro open a map and select Add Data from the Map Tab. Select Data at the top of the drop down menu. The Add Data dialog box will open on the left side of the box, expand Portal if necessary, then select Living Atlas. Type "insects and disease" in the search box, browse to the layer then click OK.In ArcGIS Pro you can use raster functions to create your own custom extracts of the data. Imagery layers provide fast, powerful inputs to geoprocessing tools, models, or Python scripts in Pro. For example, Zonal Statistics as Table tool can be used to summarize risk of tree mortality across several watersheds, counties, or other areas that you may be interested in such as areas near homes.In ArcGIS Online you can change then layer's symbology in the image display control, set the layer's transparency, and control the visible scale range.The ArcGIS Living Atlas of the World provides an easy way to explore many other beautiful and authoritative maps on hundreds of topics like this one.

Important Note: This item is in mature support as of February 2023 and will be retired in December 2025. A new version of this item is available for your use. Esri recommends updating your maps and apps to use the new version. This layer displays change in pixels of the Sentinel-2 10m Land Use/Land Cover product developed by Esri, Impact Observatory, and Microsoft. Available years to compare with 2021 are 2018, 2019 and 2020. By default, the layer shows all comparisons together, in effect showing what changed 2018-2021. But the layer may be changed to show one of three specific pairs of years, 2018-2021, 2019-2021, or 2020-2021.Showing just one pair of years in ArcGIS Online Map ViewerTo show just one pair of years in ArcGIS Online Map viewer, create a filter. 1. Click the filter button. 2. Next, click add expression. 3. In the expression dialogue, specify a pair of years with the ProductName attribute. Use the following example in your expression dialogue to show only places that changed between 2020 and 2021:ProductNameis2020-2021By default, places that do not change appear as a transparent symbol in ArcGIS Pro. But in ArcGIS Online Map Viewer, a transparent symbol may need to be set for these places after a filter is chosen. To do this:4. Click the styles button. 5. Under unique values click style options. 6. Click the symbol next to No Change at the bottom of the legend. 7. Click the slider next to "enable fill" to turn the symbol off.Showing just one pair of years in ArcGIS ProTo show just one pair of years in ArcGIS Pro, choose one of the layer's processing templates to single out a particular pair of years. The processing template applies a definition query that works in ArcGIS Pro. 1. To choose a processing template, right click the layer in the table of contents for ArcGIS Pro and choose properties. 2. In the dialogue that comes up, choose the tab that says processing templates. 3. On the right where it says processing template, choose the pair of years you would like to display. The processing template will stay applied for any analysis you may want to perform as well.How the change layer was created, combining LULC classes from two yearsImpact Observatory, Esri, and Microsoft used artificial intelligence to classify the world in 10 Land Use/Land Cover (LULC) classes for the years 2017-2021. Mosaics serve the following sets of change rasters in a single global layer: Change between 2018 and 2021Change between 2019 and 2021Change between 2020 and 2021To make this change layer, Esri used an arithmetic operation combining the cells from a source year and 2021 to make a change index value. ((from year * 16) + to year) In the example of the change between 2020 and 2021, the from year (2020) was multiplied by 16, then added to the to year (2021). Then the combined number is served as an index in an 8 bit unsigned mosaic with an attribute table which describes what changed or did not change in that timeframe. Variable mapped: Change in land cover between 2018, 2019, or 2020 and 2021 Data Projection: Universal Transverse Mercator (UTM)Mosaic Projection: WGS84Extent: GlobalSource imagery: Sentinel-2Cell Size: 10m (0.00008983152098239751 degrees)Type: ThematicSource: Esri Inc.Publication date: January 2022What can you do with this layer?Global LULC maps provide information on conservation planning, food security, and hydrologic modeling, among other things. This dataset can be used to visualize land cover anywhere on Earth. This layer can also be used in analyses that require land cover input. For example, the Zonal Statistics tools allow a user to understand the composition of a specified area by reporting the total estimates for each of the classes. Land Cover processingThis map was produced by a deep learning model trained using over 5 billion hand-labeled Sentinel-2 pixels, sampled from over 20,000 sites distributed across all major biomes of the world. The underlying deep learning model uses 6 bands of Sentinel-2 surface reflectance data: visible blue, green, red, near infrared, and two shortwave infrared bands. To create the final map, the model is run on multiple dates of imagery throughout the year, and the outputs are composited into a final representative map. Processing platformSentinel-2 L2A/B data was accessed via Microsoft’s Planetary Computer and scaled using Microsoft Azure Batch.Class definitions1. WaterAreas where water was predominantly present throughout the year; may not cover areas with sporadic or ephemeral water; contains little to no sparse vegetation, no rock outcrop nor built up features like docks; examples: rivers, ponds, lakes, oceans, flooded salt plains.2. TreesAny significant clustering of tall (~15-m or higher) dense vegetation, typically with a closed or dense canopy; examples: wooded vegetation,
clusters of dense tall vegetation within savannas, plantations, swamp or mangroves (dense/tall vegetation with ephemeral water or canopy too thick to detect water underneath).4. Flooded vegetationAreas of any type of vegetation with obvious intermixing of water throughout a majority of the year; seasonally flooded area that is a mix of grass/shrub/trees/bare ground; examples: flooded mangroves, emergent vegetation, rice paddies and other heavily irrigated and inundated agriculture.5. CropsHuman planted/plotted cereals, grasses, and crops not at tree height; examples: corn, wheat, soy, fallow plots of structured land.7. Built AreaHuman made structures; major road and rail networks; large homogenous impervious surfaces including parking structures, office buildings and residential housing; examples: houses, dense villages / towns / cities, paved roads, asphalt.8. Bare groundAreas of rock or soil with very sparse to no vegetation for the entire year; large areas of sand and deserts with no to little vegetation; examples: exposed rock or soil, desert and sand dunes, dry salt flats/pans, dried lake beds, mines.9. Snow/IceLarge homogenous areas of permanent snow or ice, typically only in mountain areas or highest latitudes; examples: glaciers, permanent snowpack, snow fields. 10. CloudsNo land cover information due to persistent cloud cover.11. Rangeland Open areas covered in homogenous grasses with little to no taller vegetation; wild cereals and grasses with no obvious human plotting (i.e., not a plotted field); examples: natural meadows and fields with sparse to no tree cover, open savanna with few to no trees, parks/golf courses/lawns, pastures. Mix of small clusters of plants or single plants dispersed on a landscape that shows exposed soil or rock; scrub-filled clearings within dense forests that are clearly not taller than trees; examples: moderate to sparse cover of bushes, shrubs and tufts of grass, savannas with very sparse grasses, trees or other plants.CitationKarra, Kontgis, et al. “Global land use/land cover with Sentinel-2 and deep learning.” IGARSS 2021-2021 IEEE International Geoscience and Remote Sensing Symposium. IEEE, 2021.AcknowledgementsTraining data for this project makes use of the National Geographic Society Dynamic World training dataset, produced for the Dynamic World Project by National Geographic Society in partnership with Google and the World Resources Institute.For questions please email environment@esri.com

WorldClim 2.1 provides downscaled estimates of climate variables as monthly means over the period of 1970-2000 based on interpolated station measurements. Here we provide analytical image services of precipitation for each month along with an annual mean. Each time step is accessible from a processing template.Time Extent: Monthly/Annual 1970-2000Units: mm/monthCell Size: 2.5 minutes (~5 km)Source Type: StretchedPixel Type: 16 Bit IntegerData Projection: GCS WGS84Mosaic Projection: GCS WGS84Extent: GlobalSource: WorldClim v2.1Using Processing Templates to Access TimeThere are 13 processing templates applied to this service, each providing access to the 12 monthly and 1 annual mean precipitation layers. To apply these in ArcGIS Online, select the Image Display options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left-hand menu. From the Processing Template pull down menu, select the version to display.What can you do with this layer?This layer may be added to maps to visualize and quickly interrogate each pixel value. The pop-up provides a graph of the time series along with the calculated annual mean value.This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro and an area count of precipitation may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from month to month to show seasonal patterns.To calculate precipitation by land area, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Source Data: The datasets behind this layer were extracted from GeoTIF files produced by WorldClim at 2.5 minutes resolution. The mean of the 12 GeoTIFs was calculated (annual), and the 13 rasters were converted to Cloud Optimized GeoTIFF format and added to a mosaic dataset.Citation: Fick, S.E. and R.J. Hijmans, 2017. WorldClim 2: new 1km spatial resolution climate surfaces for global land areas. International Journal of Climatology 37 (12): 4302-4315.

Pockmarks are defined as depressions on the seabed and are usually formed by fluid expulsions. Recently discovered, pockmarks along the Aquitaine slope within the French EEZ, were manually mapped although two semi-automated methods were tested without convincing results. In order to potentially highlight different groups and possibly discriminate the nature of the fluids involved in their formation and evolution, a morphological study was conducted, mainly based on multibeam data and in particular bathymetry from the marine expedition GAZCOGNE1, 2013. Bathymetry and seafloor backscatter data, covering more than 3200 km², were acquired with the Kongsberg EM302 ship-borne multibeam echosounder of the R/V Le Suroît at a speed of ~8 knots, operated at a frequency of 30 kHz and calibrated with ©Sippican shots. Precision of seafloor backscatter amplitude is +/- 1 dB. Multibeam data, processed using Caraibes (©IFREMER), were gridded at 15x15 m and down to 10x10 m cells, for bathymetry and seafloor backscatter, respectively. The present table includes 11 morphological attributes extracted from a Geographical Information System project (Mercator 44°N conserved latitude in WGS84 Datum) and additional parameters related to seafloor backscatter amplitudes. Pockmark occurrence with regards to the different morphological domains is derived from a morphological analysis manually performed and based on GAZCOGNE1 and BOBGEO2 bathymetric datasets. The pockmark area and its perimeter were calculated with the “Calculate Geometry” tool of Arcmap 10.2 (©ESRI) (https://desktop.arcgis.com/en/arcmap/10.3/manage-data/tables/calculating-area-length-and-other-geometric-properties.htm). A first method to calculate pockmark internal depth developed by Gafeira et al. was tested (Gafeira J, Long D, Diaz-Doce D (2012) Semi-automated characterisation of seabed pockmarks in the central North Sea. Near Surface Geophysics 10 (4):303-315, doi:10.3997/1873-0604.2012018). This method is based on the “Fill” function from the Hydrology toolset in Spatial Analyst Toolbox Arcmap 10.2 (©ESRI), (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/fill.htm) which fills the closed depressions. The difference between filled bathymetry and initial bathymetry produces a raster grid only highlighting filled depressions. Thus, only the maximum filling values which correspond to the internal depths at the apex of the pockmark were extracted. For the second method, the internal pockmark depth was calculated with the difference between minimum and maximum bathymetry within the pockmark. Latitude and longitude of the pockmark centroid, minor and major axis lengths and major axis direction of the pockmarks were calculated inside each depression with the “Zonal Geometry as Table” tool from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI) (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/zonal-statistics.htm). Pockmark elongation was calculated as the ratio between the major and minor axis length. Cell count is the number of cells used inside each pockmark to calculate statistics (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/zonal-geometry.htm). Cell count and minimum, maximum and mean bathymetry, slope and seafloor backscatter values were calculated within each pockmark with “Zonal Statistics as Table” tool from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI). Slope was calculated from bathymetry with “Slope” function from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI) and preserves its 15 m grid size (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/slope.htm). Seafloor backscatter amplitudes (minimum, maximum and mean values) of the surrounding sediments were calculated within a 100 m buffer around the pockmark rim.

This dataset contains 50-ft contours for the Hot Springs shallowest unit of the Ouachita Mountains aquifer system potentiometric-surface map. The potentiometric-surface shows altitude at which the water level would have risen in tightly-cased wells and represents synoptic conditions during the summer of 2017. Contours were constructed from 59 water-level measurements measured in selected wells (locations in the well point dataset). Major streams and creeks were selected in the study area from the USGS National Hydrography Dataset (U.S. Geological Survey, 2017), and the spring point dataset with 18 spring altitudes calculated from 10-meter digital elevation model (DEM) data (U.S. Geological Survey, 2015; U.S. Geological Survey, 2016). After collecting, processing, and plotting the data, a potentiometric surface was generated using the interpolation method Topo to Raster in ArcMap 10.5 (Esri, 2017a). This tool is specifically designed for the creation of digital elevation models and imposes constraints that ensure a connected drainage structure and a correct representation of the surface from the provided contour data (Esri, 2017a). Once the raster surface was created, 50-ft contour interval were generated using Contour (Spatial Analyst), a spatial analyst tool (available through ArcGIS 3D Analyst toolbox) that creates a line-feature class of contours (isolines) from the raster surface (Esri, 2017b). The Topo to Raster and contouring done by ArcMap 10.5 is a rapid way to interpolate data, but computer programs do not account for hydrologic connections between groundwater and surface water. For this reason, some contours were manually adjusted based on topographical influence, a comparison with the potentiometric surface of Kresse and Hays (2009), and data-point water-level altitudes to more accurately represent the potentiometric surface. Select References: Esri, 2017a, How Topo to Raster works—Help | ArcGIS Desktop, accessed December 5, 2017, at ArcGIS Pro at http://pro.arcgis.com/en/pro-app/tool-reference/3d-analyst/how-topo-to-raster-works.htm. Esri, 2017b, Contour—Help | ArcGIS Desktop, accessed December 5, 2017, at ArcGIS Pro Raster Surface toolset at http://pro.arcgis.com/en/pro-app/tool-reference/3d-analyst/contour.htm. Kresse, T.M., and Hays, P.D., 2009, Geochemistry, Comparative Analysis, and Physical and Chemical Characteristics of the Thermal Waters East of Hot Springs National Park, Arkansas, 2006-09: U.S. Geological Survey 2009–5263, 48 p., accessed November 28, 2017, at https://pubs.usgs.gov/sir/2009/5263/. U.S. Geological Survey, 2015, USGS NED 1 arc-second n35w094 1 x 1 degree ArcGrid 2015, accessed December 5, 2017, at The National Map: Elevation at https://nationalmap.gov/elevation.html. U.S. Geological Survey, 2016, USGS NED 1 arc-second n35w093 1 x 1 degree ArcGrid 2016, accessed December 5, 2017, at The National Map: Elevation at https://nationalmap.gov/elevation.html.

This map shows the change in particulate matter 2.5 (PM 2.5) air quality data for the US between 2010 and 2016 based on NASA SEDAC gridded data. The color indicates better or worse air quality, and the size of the symbol indicates population growth.This map shows particulate matter in the air sized 2.5 micrometers of smaller (PM 2.5). The data is aggregated from NASA Socioeconomic Data and Applications Center (SEDAC) gridded data into state, county, congressional district (116th) and 50 km hex bins. The unit of measurement is micrograms per cubic meter.The data is averaged for each year and over the the 19 years to provide an overall picture of air quality in the United States, including Puerto Rico. A space time cube was performed on a multidimensional mosaic version of the data in order to derive an emerging hot spot analysis. The county and state layers provide a population-weighted PM 2.5 value to emphasize which areas have a higher human impact. Each layer has been enriched with a set of 2019 US demographic attributes (excluding Puerto Rico) apportioned to the geography in order to map patterns alongside each other. Citations:van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2018. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998-2016. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4ZK5DQS. Accessed 1 April 2020van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2016. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites. Environmental Science & Technology 50 (7): 3762-3772. https://doi.org/10.1021/acs.est.5b05833.Boundaries:50km hex bins generated using the Generate Tessellation toolStates and counties come from 2018 TIGER boundaries with coastlines clipped116th Congressional Districts come from this ArcGIS Living Atlas layerData processing notes:NASA's GeoTIFF files for 19 years (1998-2016) were first brought into ArcGIS Pro 2.5.0 and put into a multidimensional mosaic dataset.For each geography level, the following was performed: Zonal Statistics were run against the mosaic as a multidimensional layer.A Space Time Cube was created to compare the 19 years of PM 2.5 values and detect hot/cold spot patterns. To learn more about Space Time Cubes, visit this page.The Space Time Cube is processed for Emerging Hot Spots where we gain the trends and hot spot results.The Enrich tool was run to add 2019 Esri demographic and 2014-2018 ACS attributes to the geographies. Attributes such as population, poverty, minority population, and others were added to the layer.To create the population-weighted attributes on the state and county layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average PM 2.5 were multiplied.The hex bins were converted into centroids and summarized within the state and county boundaries.The summation of these values were then divided by the total population of each state/county.

This webmap is a subset of Global Landcover 1992 - 2020 Image Layer. You can access the source data from here. This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years 1992-2020. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2020Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: Annual until 2020, no updates thereafterWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

This layer is a subset of Global Landcover 1992- 2020 Layer. This layer is a time series of the annual ESA CCI (Climate Change Initiative) land cover maps of the world. ESA has produced land cover maps for the years 1992-2020. These are available at the European Space Agency Climate Change Initiative website.Time Extent: 1992-2020Cell Size: 300 meterSource Type: ThematicPixel Type: 8 Bit UnsignedData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: GlobalSource: ESA Climate Change InitiativeUpdate Cycle: Annual until 2020, no updates thereafterWhat can you do with this layer?This layer may be added to ArcGIS Online maps and applications and shown in a time series to watch a "time lapse" view of land cover change since 1992 for any part of the world. The same behavior exists when the layer is added to ArcGIS Pro.In addition to displaying all layers in a series, this layer may be queried so that only one year is displayed in a map. This layer can be used in analysis. For example, the layer may be added to ArcGIS Pro with a query set to display just one year. Then, an area count of land cover types may be produced for a feature dataset using the zonal statistics tool. Statistics may be compared with the statistics from other years to show a trend.To sum up area by land cover using this service, or any other analysis, be sure to use an equal area projection, such as Albers or Equal Earth.Different Classifications Available to MapFive processing templates are included in this layer. The processing templates may be used to display a smaller set of land cover classes.Cartographic Renderer (Default Template)Displays all ESA CCI land cover classes.*Forested lands TemplateThe forested lands template shows only forested lands (classes 50-90).Urban Lands TemplateThe urban lands template shows only urban areas (class 190).Converted Lands TemplateThe converted lands template shows only urban lands and lands converted to agriculture (classes 10-40 and 190).Simplified RendererDisplays the map in ten simple classes which match the ten simplified classes used in 2050 Land Cover projections from Clark University.Any of these variables can be displayed or analyzed by selecting their processing template. In ArcGIS Online, select the Image Display Options on the layer. Then pull down the list of variables from the Renderer options. Click Apply and Close. In ArcGIS Pro, go into the Layer Properties. Select Processing Templates from the left hand menu. From the Processing Template pull down menu, select the variable to display.Using TimeBy default, the map will display as a time series animation, one year per frame. A time slider will appear when you add this layer to your map. To see the most current data, move the time slider until you see the most current year.In addition to displaying the past quarter century of land cover maps as an animation, this time series can also display just one year of data by use of a definition query. For a step by step example using ArcGIS Pro on how to display just one year of this layer, as well as to compare one year to another, see the blog called Calculating Impervious Surface Change.Hierarchical ClassificationLand cover types are defined using the land cover classification (LCCS) developed by the United Nations, FAO. It is designed to be as compatible as possible with other products, namely GLCC2000, GlobCover 2005 and 2009.This is a heirarchical classification system. For example, class 60 means "closed to open" canopy broadleaved deciduous tree cover. But in some places a more specific type of broadleaved deciduous tree cover may be available. In that case, a more specific code 61 or 62 may be used which specifies "open" (61) or "closed" (62) cover.Land Cover ProcessingTo provide consistency over time, these maps are produced from baseline land cover maps, and are revised for changes each year depending on the best available satellite data from each period in time. These revisions were made from AVHRR 1km time series from 1992 to 1999, SPOT-VGT time series between 1999 and 2013, and PROBA-V data for years 2013, 2014 and 2015. When MERIS FR or PROBA-V time series are available, changes detected at 1 km are re-mapped at 300 m. The last step consists in back- and up-dating the 10-year baseline LC map to produce the 24 annual LC maps from 1992 to 2015.Source dataThe datasets behind this layer were extracted from NetCDF files and TIFF files produced by ESA. Years 1992-2015 were acquired from ESA CCI LC version 2.0.7 in TIFF format, and years 2016-2018 were acquired from version 2.1.1 in NetCDF format. These are downloadable from ESA with an account, after agreeing to their terms of use. https://maps.elie.ucl.ac.be/CCI/viewer/download.phpCitationESA. Land Cover CCI Product User Guide Version 2. Tech. Rep. (2017). Available at: maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdfMore technical documentation on the source datasets is available here:https://cds.climate.copernicus.eu/cdsapp#!/dataset/satellite-land-cover?tab=doc*Index of all classes in this layer:10 Cropland, rainfed11 Herbaceous cover12 Tree or shrub cover20 Cropland, irrigated or post-flooding30 Mosaic cropland (>50%) / natural vegetation (tree, shrub, herbaceous cover) (<50%)40 Mosaic natural vegetation (tree, shrub, herbaceous cover) (>50%) / cropland (<50%)50 Tree cover, broadleaved, evergreen, closed to open (>15%)60 Tree cover, broadleaved, deciduous, closed to open (>15%)61 Tree cover, broadleaved, deciduous, closed (>40%)62 Tree cover, broadleaved, deciduous, open (15-40%)70 Tree cover, needleleaved, evergreen, closed to open (>15%)71 Tree cover, needleleaved, evergreen, closed (>40%)72 Tree cover, needleleaved, evergreen, open (15-40%)80 Tree cover, needleleaved, deciduous, closed to open (>15%)81 Tree cover, needleleaved, deciduous, closed (>40%)82 Tree cover, needleleaved, deciduous, open (15-40%)90 Tree cover, mixed leaf type (broadleaved and needleleaved)100 Mosaic tree and shrub (>50%) / herbaceous cover (<50%)110 Mosaic herbaceous cover (>50%) / tree and shrub (<50%)120 Shrubland121 Shrubland evergreen122 Shrubland deciduous130 Grassland140 Lichens and mosses150 Sparse vegetation (tree, shrub, herbaceous cover) (<15%)151 Sparse tree (<15%)152 Sparse shrub (<15%)153 Sparse herbaceous cover (<15%)160 Tree cover, flooded, fresh or brakish water170 Tree cover, flooded, saline water180 Shrub or herbaceous cover, flooded, fresh/saline/brakish water190 Urban areas200 Bare areas201 Consolidated bare areas202 Unconsolidated bare areas210 Water bodies

This layer shows particulate matter in the air sized 2.5 micrometers of smaller (PM 2.5). The data is aggregated from NASA Socioeconomic Data and Applications Center (SEDAC) gridded data into country boundaries, administrative 1 boundaries, and 50 km hex bins. The unit of measurement is micrograms per cubic meter.The layer shows the annual average PM 2.5 from 1998 to 2016, highlighting if the overall mean for an area meets the World Health Organization guideline of 10 micrograms per cubic meter annually. Areas that don't meet the guideline and are above the threshold are shown in red, and areas that are lower than the guideline are in grey.The data is averaged for each year and over the the 19 years to provide an overall picture of air quality globally. Some of the things we can learn from this layer:What is the average annual PM 2.5 value over 19 years? (1998-2016)What is the annual average PM 2.5 value for each year from 1998 to 2016?What is the statistical trend for PM 2.5 over the 19 years? (downward or upward)Are there hot spots (or cold spots) of PM 2.5 over the 19 years?How many people are impacted by the air quality in an area?What is the death rate caused by the joint effects of air pollution?Choose a different attribute to symbolize in order to reveal any of the patterns above.A space time cube was performed on a multidimensional mosaic version of the data in order to derive an emerging hot spot analysis, trends, and a 19-year average. The country and administrative 1 layers provide a population-weighted PM 2.5 value to emphasize which areas have a higher human impact. Citations:van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2018. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998-2016. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4ZK5DQS. Accessed 1 April 2020van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2016. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites. Environmental Science & Technology 50 (7): 3762-3772. https://doi.org/10.1021/acs.est.5b05833.Boundaries and population figures:Antarctica is excluded from all maps because it was not included in the original NASA grids.50km hex bins generated using the Generate Tessellation tool - projected to Behrmann Equal Area projection for analysesPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Administrative boundaries from World Administrative Divisions layer from ArcGIS Living Atlas - projected to Behrmann Equal Area projection for analyses and hosted in Web MercatorSources: Garmin, CIA World FactbookPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Country boundaries from Esri 2019 10.8 Data and Maps - projected to Behrmann Equal Area projection for analyses and hosted in Web Mercator. Sources: Garmin, Factbook, CIAPopulation figures attached to the country boundaries come from the World Population Estimate 2016 Sources Living Atlas layer Data processing notes:NASA's GeoTIFF files for 19 years (1998-2016) were first brought into ArcGIS Pro 2.5.0 and put into a multidimensional mosaic dataset.For each geography level, the following was performed: Zonal Statistics were run against the mosaic as a multidimensional layer.A Space Time Cube was created to compare the 19 years of PM 2.5 values and detect hot/cold spot patterns. To learn more about Space Time Cubes, visit this page.The Space Time Cube is processed for Emerging Hot Spots where we gain the trends and hot spot results.The layers are hosted in Web Mercator Auxillary Sphere projection, but were processed using an equal area projection: Behrmann. If using this layer for analysis, it is recommended to start by projecting the data back to Behrmann.The country and administrative layer were dissolved and joined with population figures in order to visualize human impact.The dissolve tool ensures that each geographic area is only symbolized once within the map.Country boundaries were generalized post-analysis for visualization purposes. The tolerance used was 700m. If performing analysis with this layer, find detailed country boundaries in ArcGIS Living Atlas. To create the population-weighted attributes on the country and Admin 1 layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average PM 2.5 were multiplied.The hex bins were converted into centroids and the PM2.5 and population figures were summarized within the country and Admin 1 boundaries.The summation of the PM 2.5 values were then divided by the total population of each geography. This population value was determined by summarizing the population values from the hex bins within each geography.Some artifacts in the hex bin layer as a result of the input NASA rasters. Because the gridded surface is created from multiple satellites, there are strips within some areas that are a result of satellite paths. Some areas also have more of a continuous pattern between hex bins as a result of the input rasters.Within the country layer, an air pollution attributable death rate is included. 2016 figures are offered by the World Health Organization (WHO). Values are offered as a mean, upper value, lower value, and also offered as age standardized. Values are for deaths caused by all possible air pollution related diseases, for both sexes, and all age groups. For more information visit this page, and here for methodology. According to WHO, the world average was 95 deaths per 100,000 people.To learn the techniques used in this analysis, visit the Learn ArcGIS lesson Investigate Pollution Patterns with Space-Time Analysis by Esri's Kevin Bulter and Lynne Buie.

Although there are a large number of software products available for calculating landscape metrics (e.g. FRAGSTATS, landscapemetrics package in R) no tools are currently available (to my knowledge) that calculate landscape metrics directly in ArcGIS Pro. Moreover, many, if not most, landscape metrics were designed with vector data in mind, but most software calculates landscape metrics from raster data due to processing time and complexity. Scaling landscape metrics can also be tedious in some instances. This toolbox was designed to calculate attributes of patches that are easily calculated on polygons in ArcGIS (i.e. area, number of patches, Landscape Shape Index, edge density, patch size, distance to the nearest patch) and scales those calculations to coarser resolutions using Block Statistics. The tool also summarizes the relationships among metrics by using Principal Component Analysis and correlation matrices to assess relationships among variables. All variables are output to a single folder.

Pockmarks are defined as depressions on the seabed and are usually formed by fluid expulsions. Recently discovered, pockmarks along the Aquitaine slope within the French EEZ, were manually mapped although two semi-automated methods were tested without convincing results. In order to potentially highlight different groups and possibly discriminate the nature of the fluids involved in their formation and evolution, a morphological study was conducted, mainly based on multibeam data and in particular bathymetry from the marine expedition GAZCOGNE1, 2013. Bathymetry and seafloor backscatter data, covering more than 3200 km², were acquired with the Kongsberg EM302 ship-borne multibeam echosounder of the R/V Le Suroît at a speed of ~8 knots, operated at a frequency of 30 kHz and calibrated with ©Sippican shots. Precision of seafloor backscatter amplitude is +/- 1 dB. Multibeam data, processed using Caraibes (©IFREMER), were gridded at 15x15 m and down to 10x10 m cells, for bathymetry and seafloor backscatter, respectively. The present table includes 11 morphological attributes extracted from a Geographical Information System project (Mercator 44°N conserved latitude in WGS84 Datum) and additional parameters related to seafloor backscatter amplitudes. Pockmark occurrence with regards to the different morphological domains is derived from a morphological analysis manually performed and based on GAZCOGNE1 and BOBGEO2 bathymetric datasets. The pockmark area and its perimeter were calculated with the “Calculate Geometry” tool of Arcmap 10.2 (©ESRI) (https://desktop.arcgis.com/en/arcmap/10.3/manage-data/tables/calculating-area-length-and-other-geometric-properties.htm). A first method to calculate pockmark internal depth developed by Gafeira et al. was tested (Gafeira J, Long D, Diaz-Doce D (2012) Semi-automated characterisation of seabed pockmarks in the central North Sea. Near Surface Geophysics 10 (4):303-315, doi:10.3997/1873-0604.2012018). This method is based on the “Fill” function from the Hydrology toolset in Spatial Analyst Toolbox Arcmap 10.2 (©ESRI), (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/fill.htm) which fills the closed depressions. The difference between filled bathymetry and initial bathymetry produces a raster grid only highlighting filled depressions. Thus, only the maximum filling values which correspond to the internal depths at the apex of the pockmark were extracted. For the second method, the internal pockmark depth was calculated with the difference between minimum and maximum bathymetry within the pockmark. Latitude and longitude of the pockmark centroid, minor and major axis lengths and major axis direction of the pockmarks were calculated inside each depression with the “Zonal Geometry as Table” tool from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI) (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/zonal-statistics.htm). Pockmark elongation was calculated as the ratio between the major and minor axis length. Cell count is the number of cells used inside each pockmark to calculate statistics (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/zonal-geometry.htm). Cell count and minimum, maximum and mean bathymetry, slope and seafloor backscatter values were calculated within each pockmark with “Zonal Statistics as Table” tool from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI). Slope was calculated from bathymetry with “Slope” function from Spatial Analyst Toolbox in ArcGIS 10.2 (©ESRI) and preserves its 15 m grid size (https://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/slope.htm). Seafloor backscatter amplitudes (minimum, maximum and mean values) of the surrounding sediments were calculated within a 100 m buffer around the pockmark rim.

This layer shows particulate matter in the air sized 2.5 micrometers of smaller (PM 2.5). The data is aggregated from NASA Socioeconomic Data and Applications Center (SEDAC) gridded data into country boundaries, administrative 1 boundaries, and 50 km hex bins. The unit of measurement is micrograms per cubic meter.The layer shows the annual average PM 2.5 from 1998 to 2016, highlighting if the overall mean for an area meets the World Health Organization guideline of 10 micrograms per cubic meter annually. Areas that don't meet the guideline and are above the threshold are shown in red, and areas that are lower than the guideline are in grey.The data is averaged for each year and over the the 19 years to provide an overall picture of air quality globally. Some of the things we can learn from this layer:What is the average annual PM 2.5 value over 19 years? (1998-2016)What is the annual average PM 2.5 value for each year from 1998 to 2016?What is the statistical trend for PM 2.5 over the 19 years? (downward or upward)Are there hot spots (or cold spots) of PM 2.5 over the 19 years?How many people are impacted by the air quality in an area?What is the death rate caused by the joint effects of air pollution?Choose a different attribute to symbolize in order to reveal any of the patterns above.A space time cube was performed on a multidimensional mosaic version of the data in order to derive an emerging hot spot analysis, trends, and a 19-year average. The country and administrative 1 layers provide a population-weighted PM 2.5 value to emphasize which areas have a higher human impact. Citations:van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2018. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998-2016. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4ZK5DQS. Accessed 1 April 2020van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2016. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites. Environmental Science & Technology 50 (7): 3762-3772. https://doi.org/10.1021/acs.est.5b05833.Boundaries and population figures:Antarctica is excluded from all maps because it was not included in the original NASA grids.50km hex bins generated using the Generate Tessellation tool - projected to Behrmann Equal Area projection for analysesPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Administrative boundaries from World Administrative Divisions layer from ArcGIS Living Atlas - projected to Behrmann Equal Area projection for analyses and hosted in Web MercatorSources: Garmin, CIA World FactbookPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Country boundaries from Esri 2019 10.8 Data and Maps - projected to Behrmann Equal Area projection for analyses and hosted in Web Mercator. Sources: Garmin, Factbook, CIAPopulation figures attached to the country boundaries come from the World Population Estimate 2016 Sources Living Atlas layer Data processing notes:NASA's GeoTIFF files for 19 years (1998-2016) were first brought into ArcGIS Pro 2.5.0 and put into a multidimensional mosaic dataset.For each geography level, the following was performed: Zonal Statistics were run against the mosaic as a multidimensional layer.A Space Time Cube was created to compare the 19 years of PM 2.5 values and detect hot/cold spot patterns. To learn more about Space Time Cubes, visit this page.The Space Time Cube is processed for Emerging Hot Spots where we gain the trends and hot spot results.The layers are hosted in Web Mercator Auxillary Sphere projection, but were processed using an equal area projection: Behrmann. If using this layer for analysis, it is recommended to start by projecting the data back to Behrmann.The country and administrative layer were dissolved and joined with population figures in order to visualize human impact.The dissolve tool ensures that each geographic area is only symbolized once within the map.Country boundaries were generalized post-analysis for visualization purposes. The tolerance used was 700m. If performing analysis with this layer, find detailed country boundaries in ArcGIS Living Atlas. To create the population-weighted attributes on the country and Admin 1 layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average PM 2.5 were multiplied.The hex bins were converted into centroids and the PM2.5 and population figures were summarized within the country and Admin 1 boundaries.The summation of the PM 2.5 values were then divided by the total population of each geography. This population value was determined by summarizing the population values from the hex bins within each geography.Some artifacts in the hex bin layer as a result of the input NASA rasters. Because the gridded surface is created from multiple satellites, there are strips within some areas that are a result of satellite paths. Some areas also have more of a continuous pattern between hex bins as a result of the input rasters.Within the country layer, an air pollution attributable death rate is included. 2016 figures are offered by the World Health Organization (WHO). Values are offered as a mean, upper value, lower value, and also offered as age standardized. Values are for deaths caused by all possible air pollution related diseases, for both sexes, and all age groups. For more information visit this page, and here for methodology. According to WHO, the world average was 95 deaths per 100,000 people.To learn the techniques used in this analysis, visit the Learn ArcGIS lesson Investigate Pollution Patterns with Space-Time Analysis by Esri's Kevin Bulter and Lynne Buie.

This is a collection of all GPS- and computer-generated geospatial data specific to the Alpine Treeline Warming Experiment (ATWE), located on Niwot Ridge, Colorado, USA. The experiment ran between 2008 and 2016, and consisted of three sites spread across an elevation gradient. Geospatial data for all three experimental sites and cone/seed collection locations are included in this package. ––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––––– Geospatial files include cone collection, experimental site, seed trap, and other GPS location/terrain data. File types include ESRI shapefiles, ESRI grid files or Arc/Info binary grids, TIFFs (.tif), and keyhole markup language (.kml) files. Trimble-imported data include plain text files (.txt), Trimble COR (CorelDRAW) files, and Trimble SSF (Standard Storage Format) files. Microsoft Excel (.xlsx) and comma-separated values (.csv) files corresponding to the attribute tables of many files within this package are also included. A complete list of files can be found in this document in the “Data File Organization” section in the included Data User's Guide. Maps are also included in this data package for reference and use. These maps are separated into two categories, 2021 maps and legacy maps, which were made in 2010. Each 2021 map has one copy in portable network graphics (.png) format, and the other in .pdf format. All legacy maps are in .pdf format. .png image files can be opened with any compatible programs, such as Preview (Mac OS) and Photos (Windows). All GIS files were imported into geopackages (.gpkg) using QGIS, and double-checked for compatibility and data/attribute integrity using ESRI ArcGIS Pro. Note that files packaged within geopackages will open in ArcGIS Pro with “main.” preceding each file name, and an extra column named “geom” defining geometry type in the attribute table. The contents of each geospatial file remain intact, unless otherwise stated in “niwot_geospatial_data_list_07012021.pdf/.xlsx”. This list of files can be found as an .xlsx and a .pdf in this archive. As an open-source file format, files within gpkgs (TIFF, shapefiles, ESRI grid or “Arc/Info Binary”) can be read using both QGIS and ArcGIS Pro, and any other geospatial softwares. Text and .csv files can be read using TextEdit/Notepad/any simple text-editing software; .csv’s can also be opened using Microsoft Excel and R. .kml files can be opened using Google Maps or Google Earth, and Trimble files are most compatible with Trimble’s GPS Pathfinder Office software. .xlsx files can be opened using Microsoft Excel. PDFs can be opened using Adobe Acrobat Reader, and any other compatible programs. A selection of original shapefiles within this archive were generated using ArcMap with associated FGDC-standardized metadata (xml file format). We are including these original files because they contain metadata only accessible using ESRI programs at this time, and so that the relationship between shapefiles and xml files is maintained. Individual xml files can be opened (without a GIS-specific program) using TextEdit or Notepad. Since ESRI’s compatibility with FGDC metadata has changed since the generation of these files, many shapefiles will require upgrading to be compatible with ESRI’s latest versions of geospatial software. These details are also noted in the “niwot_geospatial_data_list_07012021” file.

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.

This layer shows particulate matter in the air sized 2.5 micrometers of smaller (PM 2.5). The data is aggregated from NASA Socioeconomic Data and Applications Center (SEDAC) gridded data into country boundaries, administrative 1 boundaries, and 50 km hex bins. The unit of measurement is micrograms per cubic meter.The layer shows the annual average PM 2.5 from 1998 to 2016, highlighting if the overall mean for an area meets the World Health Organization guideline of 10 micrograms per cubic meter annually. Areas that don't meet the guideline and are above the threshold are shown in red, and areas that are lower than the guideline are in grey.The data is averaged for each year and over the the 19 years to provide an overall picture of air quality globally. Some of the things we can learn from this layer:What is the average annual PM 2.5 value over 19 years? (1998-2016)What is the annual average PM 2.5 value for each year from 1998 to 2016?What is the statistical trend for PM 2.5 over the 19 years? (downward or upward)Are there hot spots (or cold spots) of PM 2.5 over the 19 years?How many people are impacted by the air quality in an area?What is the death rate caused by the joint effects of air pollution?Choose a different attribute to symbolize in order to reveal any of the patterns above.A space time cube was performed on a multidimensional mosaic version of the data in order to derive an emerging hot spot analysis, trends, and a 19-year average. The country and administrative 1 layers provide a population-weighted PM 2.5 value to emphasize which areas have a higher human impact. Citations:van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2018. Global Annual PM2.5 Grids from MODIS, MISR and SeaWiFS Aerosol Optical Depth (AOD) with GWR, 1998-2016. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC). https://doi.org/10.7927/H4ZK5DQS. Accessed 1 April 2020van Donkelaar, A., R. V. Martin, M. Brauer, N. C. Hsu, R. A. Kahn, R. C. Levy, A. Lyapustin, A. M. Sayer, and D. M. Winker. 2016. Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites. Environmental Science & Technology 50 (7): 3762-3772. https://doi.org/10.1021/acs.est.5b05833.Boundaries and population figures:Antarctica is excluded from all maps because it was not included in the original NASA grids.50km hex bins generated using the Generate Tessellation tool - projected to Behrmann Equal Area projection for analysesPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Administrative boundaries from World Administrative Divisions layer from ArcGIS Living Atlas - projected to Behrmann Equal Area projection for analyses and hosted in Web MercatorSources: Garmin, CIA World FactbookPopulation figures generated using Zonal Statistics from the World Population Estimate 2016 layer from ArcGIS Living Atlas.Country boundaries from Esri 2019 10.8 Data and Maps - projected to Behrmann Equal Area projection for analyses and hosted in Web Mercator. Sources: Garmin, Factbook, CIAPopulation figures attached to the country boundaries come from the World Population Estimate 2016 Sources Living Atlas layer Data processing notes:NASA's GeoTIFF files for 19 years (1998-2016) were first brought into ArcGIS Pro 2.5.0 and put into a multidimensional mosaic dataset.For each geography level, the following was performed: Zonal Statistics were run against the mosaic as a multidimensional layer.A Space Time Cube was created to compare the 19 years of PM 2.5 values and detect hot/cold spot patterns. To learn more about Space Time Cubes, visit this page.The Space Time Cube is processed for Emerging Hot Spots where we gain the trends and hot spot results.The layers are hosted in Web Mercator Auxillary Sphere projection, but were processed using an equal area projection: Behrmann. If using this layer for analysis, it is recommended to start by projecting the data back to Behrmann.The country and administrative layer were dissolved and joined with population figures in order to visualize human impact.The dissolve tool ensures that each geographic area is only symbolized once within the map.Country boundaries were generalized post-analysis for visualization purposes. The tolerance used was 700m. If performing analysis with this layer, find detailed country boundaries in ArcGIS Living Atlas. To create the population-weighted attributes on the country and Admin 1 layers, the hex value population values were used to create the weighting. Within each hex bin, the total population figure and average PM 2.5 were multiplied.The hex bins were converted into centroids and the PM2.5 and population figures were summarized within the country and Admin 1 boundaries.The summation of the PM 2.5 values were then divided by the total population of each geography. This population value was determined by summarizing the population values from the hex bins within each geography.Some artifacts in the hex bin layer as a result of the input NASA rasters. Because the gridded surface is created from multiple satellites, there are strips within some areas that are a result of satellite paths. Some areas also have more of a continuous pattern between hex bins as a result of the input rasters.Within the country layer, an air pollution attributable death rate is included. 2016 figures are offered by the World Health Organization (WHO). Values are offered as a mean, upper value, lower value, and also offered as age standardized. Values are for deaths caused by all possible air pollution related diseases, for both sexes, and all age groups. For more information visit this page, and here for methodology. According to WHO, the world average was 95 deaths per 100,000 people.To learn the techniques used in this analysis, visit the Learn ArcGIS lesson Investigate Pollution Patterns with Space-Time Analysis by Esri's Kevin Bulter and Lynne Buie.

The World Terrestrial Ecosystems map classifies the world into areas of similar climate, landform, and land cover, which form the basic components of any terrestrial ecosystem structure. This map is important because it uses objectively derived and globally consistent data to characterize the ecosystems at a much finer spatial resolution (250-m) than existing ecoregionalizations, and a much finer thematic resolution (431 classes) than existing global land cover products. This item was updated on Apr 14, 2023 to distinguish between Boreal and Polar climate regions in the terrestrial ecosystems. Cell Size: 250-meter Source Type: ThematicPixel Type: 16 Bit UnsignedData Projection: GCS WGS84Extent: GlobalSource: USGS, The Nature Conservancy, EsriUpdate Cycle: NoneWhat can you do with this layer?This map allows you to query the land surface pixels and returns the values of all the input parameters (landform type, landcover/vegetation type, climate region) and the name of the terrestrial ecosystem at that location.This layer can be used in analysis at global and local regions. However, for large scale spatial analysis, we have also provided an ArcGIS Pro Package that contains the original raster data with multiple table attributes. For simple mapping applications, there is also a raster tile layer. This layer can be combined with the World Protected Areas Database to assess the types of ecosystems that are protected, and progress towards meeting conservation goals. The WDPA layer updates monthly from the United Nations Environment Programme.Developing the World Terrestrial EcosystemsWorld Terrestrial Ecosystems map was produced by adopting and modifying the Intergovernmental Panel on Climate Change (IPCC) approach on the definition of Terrestrial Ecosystems and development of standardized global climate regions using the values of environmental moisture regime and temperature regime. We then combined the values of Global Climate Regions, Landforms and matrix-forming vegetation assemblage or land use, using the ArcGIS Combine tool (Spatial Analyst) to produce World Ecosystems Dataset. This combination resulted of 431 World Ecosystems classes.Each combination was assigned a color using an algorithm that blended traditional color schemes for each of the three components. Every pixel in this map is symbolized by a combination of values for each of these fields.The work from this collaboration is documented in the publication:Sayre et al. 2020. An assessment of the representation of ecosystems in global protected areas using new maps of World Climate Regions and World Ecosystems - Global Ecology and Conservation More information about World Terrestrial Ecosystems can be found in this Story Map.

The data here are large wood data collected in the Summer of 2022 in the Frying Pan River Basin in the Sawatch Mountains of Central Colorado. There are six datasets included or referenced here. The first is general geomorphic and watershed characteristics of the stream reaches surveyed. The second is data from field reports to the CAIC. The third is topographic data for the studied avalanche pathways. The fourth is summary data of the wood volumes within each surveyed reach. The fifth is the unprocessed raw data for all wood jams and individual pieces surveyed. The sixth is a table of literature-derived annual recruitment rates for mechanisms common to mountain streams. Data may also be accessed via the Dryad data repository as linked in the data accessibility statement. Raw data relate several wood jam and individual piece properties, including length, width, and depth of the former, and length and diameter of the latter. Data also indicate several other wood characteristics, such as piece orientation, stability and decay class, and the presence of a rootwad. Finally, data include information about the geomorphic impact of each surveyed piece and jam. Data were collected to examine research questions related to in-stream wood load volumes supplied by snow avalanches and the resultant geomorphic impacts. Methods Raw data were collected in the field by two trained observers. Wood loads were measured in the field using a census approach for all wood jams and individual wood pieces within the bankfull channel. Jams were identified as accumulations with three or more contiguous pieces; jam volume was quantified by measuring the length, width, and height of a rectangular prism fit to the dimensions of the jam and visually estimating porosity. Porosity was consistently estimated by two independent observers to minimize systematic bias. For individual wood pieces, diameter and length were measured and then used to compute volume via the formula of a cylinder. All measurements were made for wood at least partially contained within the bankfull channel, which was visually estimated in the field using topography (e.g., slope breaks). For wood pieces or jams that extended laterally beyond the bankfull dimensions, those portions outside were excluded from measurements. Measured wood loads were normalized by stream surface area (in ha) for comparisons between reaches with varied bankfull widths and lengths. These data have not been processed other than the above related volume calculations, which were then summarized across sites and watersheds (which have yielded the wood_volume_summary.csv data also presented here). Data regarding topographic and vegetation characteristics of the studied avalanche pathways were obtained via freely available remote datasets. These include a 1-m DEM for the study area available from USGS EarthExplorer (https://earthexplorer.usgs.gov/) and vegetation data from the LANDFIRE program (https://landfire.gov/). Canopy cover comes from the Existing Vegetation Cover raster in the LANDFIRE 2016 dataset. Planform curvature and slope angle were dervied from a 1-m DEM using the Curvature and Slope functions in ArcGIS Pro Version 2.8. Median was calculated for each pathway area using the Zonal Statistics function in ArcGIS Pro. Avalanche report data was gathered from Colorado Avalanche Information Center field reports (https://avalanche.state.co.us/observations/view-field-reports), restricting the "area" field to Sawatch and Aspen and the dates from 2019-03-01 to 2019-03-15. Literature data were gathered using a publication database/search engine (scholar.google.com) to find relevant sources via keyword sources. Data were then processed using the information in each publication to determine rates in units of m3/ha/yr. NAs in all datasets mean that the measure in question was inapplicable to the parameter under consideration. An open-source R code is provided to re-create data processing and figure creation.

Please note that this dataset is not an official City of Toronto land use dataset. It was created for personal and academic use using City of Toronto Land Use Maps (2019) found on the City of Toronto Official Plan website at https://www.toronto.ca/city-government/planning-development/official-plan-guidelines/official-plan/official-plan-maps-copy, along with the City of Toronto parcel fabric (Property Boundaries) found at https://open.toronto.ca/dataset/property-boundaries/ and Statistics Canada Census Dissemination Blocks level boundary files (2016). The property boundaries used were dated November 11, 2021. Further detail about the City of Toronto's Official Plan, consolidation of the information presented in its online form, and considerations for its interpretation can be found at https://www.toronto.ca/city-government/planning-development/official-plan-guidelines/official-plan/ Data Creation Documentation and Procedures Software Used The spatial vector data were created using ArcGIS Pro 2.9.0 in December 2021. PDF File Conversions Using Adobe Acrobat Pro DC software, the following downloaded PDF map images were converted to TIF format. 9028-cp-official-plan-Map-14_LandUse_AODA.pdf 9042-cp-official-plan-Map-22_LandUse_AODA.pdf 9070-cp-official-plan-Map-20_LandUse_AODA.pdf 908a-cp-official-plan-Map-13_LandUse_AODA.pdf 978e-cp-official-plan-Map-17_LandUse_AODA.pdf 97cc-cp-official-plan-Map-15_LandUse_AODA.pdf 97d4-cp-official-plan-Map-23_LandUse_AODA.pdf 97f2-cp-official-plan-Map-19_LandUse_AODA.pdf 97fe-cp-official-plan-Map-18_LandUse_AODA.pdf 9811-cp-official-plan-Map-16_LandUse_AODA.pdf 982d-cp-official-plan-Map-21_LandUse_AODA.pdf Georeferencing and Reprojecting Data Files The original projection of the PDF maps is unknown but were most likely published using MTM Zone 10 EPSG 2019 as per many of the City of Toronto's many datasets. They could also have possibly been published in UTM Zone 17 EPSG 26917 The TIF images were georeferenced in ArcGIS Pro using this projection with very good results. The images were matched against the City of Toronto's Centreline dataset found here The resulting TIF files and their supporting spatial files include: TOLandUseMap13.tfwx TOLandUseMap13.tif TOLandUseMap13.tif.aux.xml TOLandUseMap13.tif.ovr TOLandUseMap14.tfwx TOLandUseMap14.tif TOLandUseMap14.tif.aux.xml TOLandUseMap14.tif.ovr TOLandUseMap15.tfwx TOLandUseMap15.tif TOLandUseMap15.tif.aux.xml TOLandUseMap15.tif.ovr TOLandUseMap16.tfwx TOLandUseMap16.tif TOLandUseMap16.tif.aux.xml TOLandUseMap16.tif.ovr TOLandUseMap17.tfwx TOLandUseMap17.tif TOLandUseMap17.tif.aux.xml TOLandUseMap17.tif.ovr TOLandUseMap18.tfwx TOLandUseMap18.tif TOLandUseMap18.tif.aux.xml TOLandUseMap18.tif.ovr TOLandUseMap19.tif TOLandUseMap19.tif.aux.xml TOLandUseMap19.tif.ovr TOLandUseMap20.tfwx TOLandUseMap20.tif TOLandUseMap20.tif.aux.xml TOLandUseMap20.tif.ovr TOLandUseMap21.tfwx TOLandUseMap21.tif TOLandUseMap21.tif.aux.xml TOLandUseMap21.tif.ovr TOLandUseMap22.tfwx TOLandUseMap22.tif TOLandUseMap22.tif.aux.xml TOLandUseMap22.tif.ovr TOLandUseMap23.tfwx TOLandUseMap23.tif TOLandUseMap23.tif.aux.xml TOLandUseMap23.tif.ov Ground control points were saved for all georeferenced images. The files are the following: map13.txt map14.txt map15.txt map16.txt map17.txt map18.txt map19.txt map21.txt map22.txt map23.txt The City of Toronto's Property Boundaries shapefile, "property_bnds_gcc_wgs84.zip" were unzipped and also reprojected to EPSG 26917 (UTM Zone 17) into a new shapefile, "Property_Boundaries_UTM.shp" Mosaicing Images Once georeferenced, all images were then mosaiced into one image file, "LandUseMosaic20211220v01", within the project-generated Geodatabase, "Landuse.gdb" and exported TIF, "LandUseMosaic20211220.tif" Reclassifying Images Because the original images were of low quality and the conversion to TIF made the image colours even more inconsistent, a method was required to reclassify the images so that different land use classes could be identified. Using Deep learning Objects, the images were re-classified into useful consistent colours. Deep Learning Objects and Training The resulting mosaic was then prepared for reclassification using the Label Objects for Deep Learning tool in ArcGIS Pro. A training sample, "LandUseTrainingSamples20211220", was created in the geodatabase for all land use types as follows: Neighbourhoods Insitutional Natural Areas Core Employment Areas Mixed Use Areas Apartment Neighbourhoods Parks Roads Utility Corridors Other Open Spaces General Employment Areas Regeneration Areas Lettering (not a land use type, but an image colour (black), used to label streets). By identifying the letters, it then made the reclassification and vectorization results easier to clean up of unnecessary clutter caused by the labels of streets. Reclassification Once the training samples were created and saved, the raster was then reclassified using the Image Classification Wizard tool in ArcGIS Pro, using the Support...

Calculating the total volume of water stored in a landscape can be challenging. In addition to lakes and reservoirs, water can be stored in soil, snowpack, or even inside plants and animals, and tracking the all these different mediums is not generally possible. However, calculating the change in storage is easy - just subtract the water output from the water input. Using the GLDAS layers we can do this calculation for every month from January 2000 to the present day. The precipitation layer tells us the input to each cell and runoff plus evapotranspiration is the output. When the input is higher than the output during a given month, it means water was stored. When output is higher than input, storage is being depleted. Generally the change in storage should be close to the change in soil moisture content plus the change in snowpack, but it will not match up exactly because of the other storage mediums discussed above.Dataset SummaryThe GLDAS Change in Storage layer is a time-enabled image service that shows net monthly change in storage from 2000 to the present, measured in millimeters of water. It is calculated by NASA using the Noah land surface model, run at 0.25 degree spatial resolution using satellite and ground-based observational data from the Global Land Data Assimilation System (GLDAS-2.1). The model is run with 3-hourly time steps and aggregated into monthly averages. Review the complete list of model inputs, explore the output data (in GRIB format), and see the full Hydrology Catalog for all related data and information!Phenomenon Mapped: Change in Water StorageUnits: MillimetersTime Interval: MonthlyTime Extent: 2000/01/01 to presentCell Size: 28 kmSource Type: ScientificPixel Type: Signed IntegerData Projection: GCS WGS84Mosaic Projection: Web Mercator Auxiliary SphereExtent: Global Land SurfaceSource: NASAUpdate Cycle: SporadicWhat can you do with this layer?This layer is suitable for both visualization and analysis. It can be used in ArcGIS Online in web maps and applications and can be used in ArcGIS for Desktop. It is useful for scientific modeling, but only at global scales.In ArcGIS Pro you can use the built-in raster functions or create your own to create custom extracts of the data. Imagery layers provide fast, powerful inputs to geoprocessing tools, models, or Python scripts in Pro.Online you can filter the layer to show subsets of the data using the filter button and the layer's built-in raster functions.By applying the "Calculate Anomaly" raster function, it is possible to view these data in terms of deviation from the mean, instead of total change in storage. Mean change in storage for a given month is calculated over the entire period of record - 2000 to present.Time: This is a time-enabled layer. By default, it will show the first month from the map's time extent. Or, if time animation is disabled, a time range can be set using the layer's multidimensional settings. If you wish to calculate the average, sum, or min/max change in storage over the time extent, change the mosaic operator used to resolve overlapping pixels. In ArcGIS Online, you do this in the "Image Display Order" tab. In ArcGIS Pro, use the "Data" ribbon. In ArcMap, it is in the 'Mosaic' tab of the layer properties window. The minimum time extent is one month, and the maximum is 8 years. Important: You must switch from the cartographic renderer to the analytic renderer in the processing template tab in the layer properties window before using this layer as an input to geoprocessing tools.