Retirement Notice: 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 Viewer To 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-2021 By 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 Pro To 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 2022 What 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

Aquifers are underground layers of gravel, sand, or permeable rock that contain ground water. This water can be extracted using a well and provides an important source of water in many regions of the world.This layer, produced as part of the Ground Water Atlas of the United States, provides access to the areal extent of the principal aquifers of the United States. In areas where multiple aquifers exist at different depths below the surface only the shallowest aquifer is included.This layer does not display all areas where ground water exists. The U.S. Geologic Survey (USGS) mapped these aquifers by interpreting surface features and aquifers may extend beyond these features. Ground water areas along watercourses and ground water in unconsolidated glacial sand and gravel deposits are not included in this layer. Data on these areas are provided in the layer Aquifers of Alluvial and Glacial Origin from the USGS.Dataset SummaryPhenomenon Mapped: Aquifers of the United StatesCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United States, Hawaii, Puerto Rico, and the U.S. Virgin IslandsVisible Scale: All ScalesSource: Groundwater Atlas of the United StatesPublication Date: October 1, 2003Please note: "This dataset, published in 2003, contains the shallowest principal aquifers of the conterminous United States, Hawaii, Puerto Rico, and the U.S. Virgin Islands, portrayed as polygons. The map layer was developed as part of the effort to produce the maps published at 1:2,500,000 in the printed series Ground Water Atlas of the United States. The published maps contain base and cultural features not included in these data. Please note that the maps do not show the entire extent of an aquifer, only its subcrop or outcrop area. Refer to the metadata for a complete description of the files and how they were generated." (Source USGS)What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

Terms of Use:

Data Limitations Disclaimer

The MassDEP Estimated Sewer System Service Area Boundaries datalayer may not be complete, may contain errors, omissions, and other inaccuracies, and the data are subject to change. The user’s use of and/or reliance on the information contained in the Document (e.g. data) shall be at the user’s own risk and expense. MassDEP disclaims any responsibility for any loss or harm that may result to the user of this data or to any other person due to the user’s use of the Document.

All sewer service area delineations are estimates for broad planning purposes and should only be used as a guide. The data is not appropriate for site-specific or parcel-specific analysis. Not all properties within a sewer service area are necessarily served by the system, and some properties outside the mapped service areas could be served by the wastewater utility – please contact the relevant wastewater system. Not all service areas have been confirmed by the sewer system authorities.

This is an ongoing data development project. Attempts have been made to contact all sewer/wastewater systems, but not all have responded with information on their service area. MassDEP will continue to collect and verify this information. Some sewer service areas included in this datalayer have not been verified by the POTWs, privately-owned treatment works, GWDPs, or the municipality involved, but since many of those areas are based on information published online by the municipality, the utility, or in a publicly available report, they are included in the estimated sewer service area datalayer.

Please use the following citation to reference these data

MassDEP. Water Utility Resilience Program. 2025. Publicly-Owned Treatment Work and Non-Publicly-Owned Sewer Service Areas (PubV2024_12).

We want to learn about the data uses. If you use this dataset, please notify staff in the Water Resilience program (WURP@mass.gov).

Layers and Tables:

The MassDEP Estimated Sewer System Service Area data layer comprises two feature classes and a supporting table:

Publicly-Owned Treatment Works (POTW) Sewer Service Areas feature class SEWER_SERVICE_AREA_POTW_POLY includes polygon features for sewer service areas systems operated by publicly owned treatment works (POTWs)Non-Publicly Owned Treatment Works (NON-POTW) Sewer Service Areas feature class SEWER_SERVICE_AREA_NONPOTW_POLY includes polygon features for sewer service areas for operated by NON publicly owned treatment works (NON-POTWs)The Sewer Service Areas Unlocated Sites table SEWER_SERVICE_AREA_USL contains a list of known, unmapped active POTW and NON-POTW services areas at the time of publication.

ProductionData Universe

Effluent wastewater treatment plants in Massachusetts are permitted either through the Environmental Protection Agency’s (EPA) National Pollutant Discharge Elimination System (NPDES) surface water discharge permit program or the MassDEP Groundwater Discharge Permit Program. The WURP has delineated active service areas served by publicly and privately-owned effluent treatment works with a NPDES permit or a groundwater discharge permit.

National Pollutant Discharge Elimination System (NPDES) Permits

In the Commonwealth of Massachusetts, the EPA is the permitting authority for regulating point sources that discharge pollutants to surface waters. NPDES permits regulate wastewater discharge by limiting the quantities of pollutants to be discharged and imposing monitoring requirements and other conditions. NPDES permits are typically co-issued by EPA and the MassDEP. The limits and/or requirements in the permit ensure compliance with the Massachusetts Surface Water Quality Standards and Federal Regulations to protect public health and the aquatic environment. Areas served by effluent treatment plants with an active NPDES permit are included in this datalayer based on a master list developed by MassDEP using information sourced from the EPA’s Integrated Compliance Information System (ICIS).

Groundwater Discharge (GWD) Permits

In addition to surface water permittees, the WURP has delineated all active systems served by publicly and privately owned effluent treatment works with groundwater discharge (GWD) permits, and some inactive service areas. Groundwater discharge permits are required for systems discharging over 10,000 GPD sanitary wastewater – these include effluent treatment systems for public, district, or privately owned effluent treatment systems. Areas served by an effluent treatment plant with an active GWD permit are included in this datalayer based on lists received from MassDEP Wastewater staff.

Creation of Unique IDs for Each Service Area

The Sewer Service Area datalayer contains polygons that represent the service area of a particular wastewater system within a particular municipality. Every discharge permittee is assigned a unique NPDES permit number by EPA or a unique GWD permit identifier by MassDEP. MassDEP WURP creates a unique Sewer_ID for each service area by combining the municipal name of the municipality served with the permit number (NPDES or GWD) ascribed to the sewer that is serving that area. Some municipalities contain more than one sewer system, but each sewer system has a unique Sewer_ID. Occasionally the area served by a sewer system will overlap another town by a small amount – these small areas are generally not given a unique ID. The Estimated sewer Service Area datalayer, therefore, contains polygons with a unique Sewer_ID for each sewer service area. In addition, some municipalities will have multiple service areas being served by the same treatment plant – the Sewer_ID for these will contain additional identification, such as the name of the system, to uniquely identify each system.

Classifying System Service Areas

WURP staff reviewed the service areas for each system and, based on OWNER_TYPE, classified as either a publicly-owned treatment work (POTW) or a NON-POTW (see FAC_TYPE field). Each service area is further classified based on the population type served (see SECTOR field).

Methodologies and Data Sources

Several methodologies were used to create service area boundaries using various sources, including data received from the sewer system in response to requests for information from the MassDEP WURP project, information on file at MassDEP, and service area maps found online at municipal and wastewater system websites. When MassDEP received sewer line data rather than generalized areas, 300-foot buffers were created around the sewer lines to denote service areas and then edited to incorporate generalizations. Some municipalities submitted parcel data or address information to be used in delineating service areas. Many of the smaller GWD permitted sewer service areas were delineated using parcel boundaries related to the address on file.

Verification Process

Small-scale pdf file maps with roads and other infrastructure were sent to systems for corrections or verifications. If the system were small, such as a condominium complex or residential school, the relevant parcels were often used as the basis for the delineated service area. In towns where 97% or more of their population is served by the wastewater system and no other service area delineation was available, the town boundary was used as the service area boundary. Some towns responded to the request for information or verification of service areas by stating that the town boundary should be used since all, or nearly all, of the municipality is served by one wastewater system.

To ensure active systems are mapped, WURP staff developed two work flows. For NPDES-permitted systems, WURP staff reviewed available information on EPA’s ICIS database and created a master list of these systems. Staff will work to routinely update this master list by reviewing the ICIS database for new NPDES permits. The master list will serve as a method for identifying active systems, inactive systems, and unmapped systems. For GWD permittees, GIS staff established a direct linkage to the groundwater database, which allows for populating information into data fields and identifying active systems, inactive systems, and unmapped systems.

All unmapped systems are added to the Sewer Service Area Unlocated List (SEWER_SERVICE_AREAS_USL) for future mapping. Some service areas have not been mapped but their general location is represented by a small circle which serves as a placeholder - the location of these circles are estimated based on the general location of the treatment plant or the general estimated location of the service area - these do not represent the actual service area.

Percent Served Statistics The attribute table for the POTW sewer service areas (SEWER_SERVICE_AREA_POTW_POLY) has several fields relating to the percent of the town served by the particular system and one field describing the percent of town served by all systems in the town. The field ‘Percent AREA Served by System’ is strictly a calculation done dividing the area of the system by the total area of the town and multiplying by 100. In contrast, the field ‘Percent Served by System’, is not based on a particular calculation or source – it is an estimate based on various sources – these estimates are for planning purposes only. Data includes information from municipal websites and associated plans, the 1990 Municipal Priority list from CMR 310 14.17, the 2004 Pioneer Institute for Public Policy Research “percent on sewer” document, information contained on NPDES Permits and MassDEP Wastewater program staff input. Not all POTW systems have percent served statistics. Percentage may reflect the percentage of parcels served, the percent of area within a community served or the population served and should not be used for legal boundary definition or regulatory interpretation.

Sources of information for estimated wastewater service areas:

EEOA Water Assets

This is a collection of Digital Surface Models and Highest Hit rasters covering selected U.S. Forest Service and adjoining lands in the Southwest Region, encompassing Arizona and New Mexico. The data are presented in a time-enabled format, allowing the end-user to view available data year-by-year, or all available years at once, within a GIS system. The data encompass varying years, varying resolutions, and varying geographic extents, dependent upon available data as provided by the region. DSM and Highest Hit rasters represent elevation of Earth's surface, including its natural and human-made features, such as vegetation and buildings.The data contains an attribute table. Notable attributes that may be of interest to an end-user are:lowps: the pixel size of the source raster, given in meters.highps: the pixel size of the top-most pyramid for the raster, given in meters.beginyear: the first year of data acquisition for an individual dataset.endyear: the final year of data acquisition for an individual dataset.dataset_name: the name of the individual dataset within the collection.metadata: A URL link to a file on IIPP's Portal containing metadata pertaining to an individual dataset within the image service.resolution: The pixel size of the source raster, given in meters.Terrain-related imagery are primarily derived from Lidar, stereoscopic aerial imagery, or Interferometric Synthetic Aperture Radar datasets. Consequently, these derivatives inherit the limitations and uncertainties of the parent sensor and platform and the processing techniques used to produce the imagery. The terrain images are orthographic; they have been georeferenced and displacement due to sensor orientation and topography have been removed, producing data that combines the characteristics of an image with the geometric qualities of a map. The orthographic images show ground features in their proper positions, without the distortion characteristic of unrectified aerial or satellite imagery. Digital orthoimages produced and used within the Forest Service are developed from imagery acquired through various national and regional image acquisition programs. The resulting orthoimages can be directly applied in remote sensing, GIS and mapping applications. They serve a variety of purposes, from interim maps to references for Earth science investigations and analysis. Because of the orthographic property, an orthoimage can be used like a map for measurement of distances, angles, and areas with scale being constant everywhere. Also, they can be used as map layers in GIS or other computer-based manipulation, overlaying, and analysis. An orthoimage differs from a map in a manner of depiction of detail; on a map only selected detail is shown by conventional symbols whereas on an orthoimage all details appear just as in original aerial or satellite imagery.Tribal lands have been masked from this public service in accordance with Tribal agreements.

Overview

Empower your location data visualizations with our edge-matched polygons, even in difficult geographies.

Our self-hosted GIS data cover administrative and postal divisions with up to 6 precision levels: a zip code layer and up to 5 administrative levels. All levels follow a seamless hierarchical structure with no gaps or overlaps.

The geospatial data shapes are offered in high-precision and visualization resolution and are easily customized on-premise.

Use cases for the Global Boundaries Database (GIS data, Geospatial data)

• In-depth spatial analysis

• Clustering

• Geofencing

• Reverse Geocoding

• Reporting and Business Intelligence (BI)

Product Features

• Coherence and precision at every level

• Edge-matched polygons

• High-precision shapes for spatial analysis

• Fast-loading polygons for reporting and BI

• Multi-language support

For additional insights, you can combine the GIS data with:

• Population data: Historical and future trends

• UNLOCODE and IATA codes

• Time zones and Daylight Saving Time (DST)

Data export methodology

Our geospatial data packages are offered in variable formats, including - .shp - .gpkg - .kml - .shp - .gpkg - .kml - .geojson

All GIS data are optimized for seamless integration with popular systems like Esri ArcGIS, Snowflake, QGIS, and more.

Why companies choose our map data

• Precision at every level

• Coverage of difficult geographies

• No gaps, nor overlaps

Note: Custom geospatial data packages are available. Please submit a request via the above contact button for more details.

This layer displays a global map of land use/land cover (LULC) derived from ESA Sentinel-2 imagery at 10m resolution. Each year is generated with Impact Observatory’s deep learning AI land classification model, trained using billions of human-labeled image pixels from the National Geographic Society. The global maps are produced by applying this model to the Sentinel-2 Level-2A image collection on Microsoft’s Planetary Computer, processing over 400,000 Earth observations per year.The algorithm generates LULC predictions for nine classes, described in detail below. The year 2017 has a land cover class assigned for every pixel, but its class is based upon fewer images than the other years. The years 2018-2024 are based upon a more complete set of imagery. For this reason, the year 2017 may have less accurate land cover class assignments than the years 2018-2024. Key Properties Variable mapped: Land use/land cover in 2017, 2018, 2019, 2020, 2021, 2022, 2023, 2024Source Data Coordinate System: Universal Transverse Mercator (UTM) WGS84Service Coordinate System: Web Mercator Auxiliary Sphere WGS84 (EPSG:3857)Extent: GlobalSource imagery: Sentinel-2 L2ACell Size: 10-metersType: ThematicAttribution: Esri, Impact ObservatoryAnalysis: Optimized for analysisClass Definitions: ValueNameDescription1WaterAreas 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.2TreesAny significant clustering of tall (~15 feet 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).4Flooded 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.5CropsHuman planted/plotted cereals, grasses, and crops not at tree height; examples: corn, wheat, soy, fallow plots of structured land.7Built 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.8Bare 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.9Snow/IceLarge homogenous areas of permanent snow or ice, typically only in mountain areas or highest latitudes; examples: glaciers, permanent snowpack, snow fields.10CloudsNo land cover information due to persistent cloud cover.11RangelandOpen 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.NOTE: Land use focus does not provide the spatial detail of a land cover map. As such, for the built area classification, yards, parks, and groves will appear as built area rather than trees or rangeland classes.Usage Information and Best PracticesProcessing TemplatesThis layer includes a number of preconfigured processing templates (raster function templates) to provide on-the-fly data rendering and class isolation for visualization and analysis. Each processing template includes labels and descriptions to characterize the intended usage. This may include for visualization, for analysis, or for both visualization and analysis. VisualizationThe default rendering on this layer displays all classes.There are a number of on-the-fly renderings/processing templates designed specifically for data visualization.By default, the most recent year is displayed. To discover and isolate specific years for visualization in Map Viewer, try using the Image Collection Explorer. AnalysisIn order to leverage the optimization for analysis, the capability must be enabled by your ArcGIS organization administrator. More information on enabling this feature can be found in the ‘Regional data hosting’ section of this help doc.Optimized for analysis means this layer does not have size constraints for analysis and it is recommended for multisource analysis with other layers optimized for analysis. See this group for a complete list of imagery layers optimized for analysis.Prior to running analysis, users should always provide some form of data selection with either a layer filter (e.g. for a specific date range, cloud cover percent, mission, etc.) or by selecting specific images. To discover and isolate specific images for analysis in Map Viewer, try using the Image Collection Explorer.Zonal Statistics is a common tool used for understanding the composition of a specified area by reporting the total estimates for each of the classes. GeneralIf you are new to Sentinel-2 LULC, the Sentinel-2 Land Cover Explorer provides a good introductory user experience for working with this imagery layer. For more information, see this Quick Start Guide.Global land use/land cover maps provide information on conservation planning, food security, and hydrologic modeling, among other things. This dataset can be used to visualize land use/land cover anywhere on Earth. Classification ProcessThese maps include Version 003 of the global Sentinel-2 land use/land cover data product. It is produced by a deep learning model trained using over five 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 L2A 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 for each year.The input Sentinel-2 L2A data was accessed via Microsoft’s Planetary Computer and scaled using Microsoft Azure Batch. 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.

Access National Hydrography ProductsThe National Hydrography Dataset (NHD) is a feature-based database that interconnects and uniquely identifies the stream segments or reaches that make up the nation's surface water drainage system. NHD data was originally developed at 1:100,000-scale and exists at that scale for the whole country. This high-resolution NHD, generally developed at 1:24,000/1:12,000 scale, adds detail to the original 1:100,000-scale NHD. (Data for Alaska, Puerto Rico and the Virgin Islands was developed at high-resolution, not 1:100,000 scale.) Local resolution NHD is being developed where partners and data exist. The NHD contains reach codes for networked features, flow direction, names, and centerline representations for areal water bodies. Reaches are also defined on waterbodies and the approximate shorelines of the Great Lakes, the Atlantic and Pacific Oceans and the Gulf of Mexico. The NHD also incorporates the National Spatial Data Infrastructure framework criteria established by the Federal Geographic Data Committee.The NHD is a national framework for assigning reach addresses to water-related entities, such as industrial discharges, drinking water supplies, fish habitat areas, wild and scenic rivers. Reach addresses establish the locations of these entities relative to one another within the NHD surface water drainage network, much like addresses on streets. Once linked to the NHD by their reach addresses, the upstream/downstream relationships of these water-related entities--and any associated information about them--can be analyzed using software tools ranging from spreadsheets to geographic information systems (GIS). GIS can also be used to combine NHD-based network analysis with other data layers, such as soils, land use and population, to help understand and display their respective effects upon one another. Furthermore, because the NHD provides a nationally consistent framework for addressing and analysis, water-related information linked to reach addresses by one organization (national, state, local) can be shared with other organizations and easily integrated into many different types of applications to the benefit of all.Statements of attribute accuracy are based on accuracy statements made for U.S. Geological Survey Digital Line Graph (DLG) data, which is estimated to be 98.5 percent. One or more of the following methods were used to test attribute accuracy: manual comparison of the source with hardcopy plots; symbolized display of the DLG on an interactive computer graphic system; selected attributes that could not be visually verified on plots or on screen were interactively queried and verified on screen. In addition, software validated feature types and characteristics against a master set of types and characteristics, checked that combinations of types and characteristics were valid, and that types and characteristics were valid for the delineation of the feature. Feature types, characteristics, and other attributes conform to the Standards for National Hydrography Dataset (USGS, 1999) as of the date they were loaded into the database. All names were validated against a current extract from the Geographic Names Information System (GNIS). The entry and identifier for the names match those in the GNIS. The association of each name to reaches has been interactively checked, however, operator error could in some cases apply a name to a wrong reach.Points, nodes, lines, and areas conform to topological rules. Lines intersect only at nodes, and all nodes anchor the ends of lines. Lines do not overshoot or undershoot other lines where they are supposed to meet. There are no duplicate lines. Lines bound areas and lines identify the areas to the left and right of the lines. Gaps and overlaps among areas do not exist. All areas close.The completeness of the data reflects the content of the sources, which most often are the published USGS topographic quadrangle and/or the USDA Forest Service Primary Base Series (PBS) map. The USGS topographic quadrangle is usually supplemented by Digital Orthophoto Quadrangles (DOQs). Features found on the ground may have been eliminated or generalized on the source map because of scale and legibility constraints. In general, streams longer than one mile (approximately 1.6 kilometers) were collected. Most streams that flow from a lake were collected regardless of their length. Only definite channels were collected so not all swamp/marsh features have stream/rivers delineated through them. Lake/ponds having an area greater than 6 acres were collected. Note, however, that these general rules were applied unevenly among maps during compilation. Reach codes are defined on all features of type stream/river, canal/ditch, artificial path, coastline, and connector. Waterbody reach codes are defined on all lake/pond and most reservoir features. Names were applied from the GNIS database. Detailed capture conditions are provided for every feature type in the Standards for National Hydrography Dataset available online through https://prd-wret.s3-us-west-2.amazonaws.com/assets/palladium/production/atoms/files/NHD%201999%20Draft%20Standards%20-%20Capture%20conditions.PDF.Statements of horizontal positional accuracy are based on accuracy statements made for U.S. Geological Survey topographic quadrangle maps. These maps were compiled to meet National Map Accuracy Standards. For horizontal accuracy, this standard is met if at least 90 percent of points tested are within 0.02 inch (at map scale) of the true position. Additional offsets to positions may have been introduced where feature density is high to improve the legibility of map symbols. In addition, the digitizing of maps is estimated to contain a horizontal positional error of less than or equal to 0.003 inch standard error (at map scale) in the two component directions relative to the source maps. Visual comparison between the map graphic (including digital scans of the graphic) and plots or digital displays of points, lines, and areas, is used as control to assess the positional accuracy of digital data. Digital map elements along the adjoining edges of data sets are aligned if they are within a 0.02 inch tolerance (at map scale). Features with like dimensionality (for example, features that all are delineated with lines), with or without like characteristics, that are within the tolerance are aligned by moving the features equally to a common point. Features outside the tolerance are not moved; instead, a feature of type connector is added to join the features.Statements of vertical positional accuracy for elevation of water surfaces are based on accuracy statements made for U.S. Geological Survey topographic quadrangle maps. These maps were compiled to meet National Map Accuracy Standards. For vertical accuracy, this standard is met if at least 90 percent of well-defined points tested are within one-half contour interval of the correct value. Elevations of water surface printed on the published map meet this standard; the contour intervals of the maps vary. These elevations were transcribed into the digital data; the accuracy of this transcription was checked by visual comparison between the data and the map.

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes cotton production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Cotton ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United StatesVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Area Harvested in AcresOperations with Area HarvestedOperations with SalesProduction in BalesSales in US DollarsIrrigated Area Harvested in AcresOperations with Irrigated Area HarvestedAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

[Metadata] Flood Hazard Areas Line features for the State of Hawaii as of May, 2021. The Statewide GIS Program created the statewide layer by merging all county layers (downloaded on May 1, 2021). The National Flood Hazard Layer (NFHL) data incorporates all Flood Insurance Rate Map (FIRM) databases published by the Federal Emergency Management Agency (FEMA), and any Letters of Map Revision (LOMRs) that have been issued against those databases since their publication date. It is updated on a monthly basis. The FIRM Database is the digital, geospatial version of the flood hazard information shown on the published paper FIRMs. The FIRM Database depicts flood risk information and supporting data used to develop the risk data. The primary risk classifications used are the 1-percent-annual-chance flood event, the 0.2-percent-annual-chance flood event, and areas of minimal flood risk. The FIRM Database is derived from Flood Insurance Studies (FISs), previously published FIRMs, flood hazard analyses performed in support of the FISs and FIRMs, and new mapping data, where available. The FISs and FIRMs are published by FEMA. The NFHL is available as State or US Territory data sets. Each State or Territory data set consists of all FIRM Databases and corresponding LOMRs available on the publication date of the data set. The specification for the horizontal control of FIRM Databases is consistent with those required for mapping at a scale of 1:12,000. This file is georeferenced to the Earth's surface using the Geographic Coordinate System (GCS) and North American Datum of 1983. For more information, please refer to summary metadata: https://files.hawaii.gov/dbedt/op/gis/data/s_fld_haz_ar_line_state.pdf or contact Hawaii Statewide GIS Program, Office of Planning and Sustainable Development, State of Hawaii; PO Box 2359, Honolulu, HI 96804; (808) 587-2846; email: gis@hawaii.gov; Website: https://planning.hawaii.gov/gis.

This map features the World Ecophysiographic Land Units 2015 layer, focused on the Caribbean. Ecological Land Units are areas of distinct bioclimate, landform, lithology, and land cover that form the basic components of terrestrial ecosystem structure. Click on the map to learn more about these components for a given location.Ecological Land Units (ELUs) are areas of distinct bioclimate, landform, lithology, and land cover that form the basic components of terrestrial ecosystem structure. The ELU layer was produced by combining the values in four 250-m cell-sized rasters using the ArcGIS Combine tool (Spatial Analyst). In 2015 these four components resulted in 3,639 different combinations or ELUs, which is 284 fewer than 2014 which used older land cover and a different landform methodology.Note: This layer is designed for use as a geoprocessing input layer and to support pop-ups in ArcGIS Online. Because of the large number of unique values in the image service, the legend cannot be used in a meaningful way. Use the World Ecological Land Units Map 2015 tiled map layer for mapping and visualization. These four component datasets represent the most accurate, current, globally comprehensive, and finest spatial and thematic resolution data available for each of the four inputs. Values for each of the four input layers are listed in the table below. BioclimateLandformsLithologyLand CoverArcticPlainsUndefinedBare AreaCold DryHillsUnconsolidated SedimentSparse VegetationCold Semi-DryMountainsCarbonate Sedimentary RockGrassland, Shrub, or ScrubCold Moist Mixed Sedimentary RockMostly CroplandCold Wet Non-Carbonate Sedimentary RockMostly Needleleaf/Evergreen ForestCool Dry EvaporiteMostly Deciduous ForestCool Semi-Dry PyroclasticsSwampy or Often FloodedCool Moist Metamorphic RockArtificial or Urban AreaCool Wet Acidic VolcanicsSurface WaterHot Dry Acidic PlutonicsUndefinedHot Semi-Dry Non-Acidic Volcanics Hot Moist Non-Acidic Plutonics Hot Wet Warm Dry Warm Semi-Dry Warm Moist Warm Wet Dataset SummaryThis layer is suitable for analysis and can be used in ArcGIS Online to support pop-ups. It can be used in ArcGIS Desktop. Because of the large number of unique values in the image service it cannot be symbolized and displays as an all black layer. To use in pop-ups set the transparency to 100% and configure the pop-up. The pop-up from this layer can be combined with the World Ecological Land Units Map.Layers providing access to the four input layers used to create this map see the following links:World BioclimatesWorld Landforms Improved Hammond MethodWorld LithologyWorld Land Cover ESA 2010The ecophysiographic facets layer is available here and a layer summarizing the local diversity of the ecophysiographic facets is available here. A service is available to the data tables associated with this and other global layers. These data table services can be used by developers to create custom applications. For more information see the World Ecophysiographic Tables.The layer was created by the USGS and Esri in 2015.For more information see the publication:Sayre and others. 2014. A New Map of Global Ecological Land Units — An Ecophysiographic Stratification Approach. Washington, DC: Association of American Geographers. 46 pages. Available onlineWhat can you do with this layer?This layer is suitable for analysis and can be used in ArcGIS Online to support pop-ups. It can be used in ArcGIS Desktop. Because of the large number of unique values in the image service it can not be symbolized and displays as an all white layer. To use in pop-ups set the transparency to 100% and configure the pop-up.This layer has query, identify, and export image services available. This layer is restricted to a maximum area of 16,000 x 16,000 pixels - an area 4,000 kilometers on a side or an area approximately the size of Europe. This layer is part of a larger collection of landscape layers that you can use to perform a wide variety of mapping and analysis tasks.The Living Atlas of the World provides an easy way to explore the landscape layers and many other beautiful and authoritative maps on hundreds of topics.Geonet is a good resource for learning more about landscape layers and the Living Atlas of the World. To get started see the Living Atlas Discussion Group.The Esri Insider Blog provides an introduction to the Ecophysiographic Mapping project.

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes hay production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Hay ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United States, Alaska, and HawaiiVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Area Harvested in AcresOperations with Area HarvestedProduction in TonsAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.Many other ready-to-use layers derived from the Census of Agriculture can be found in the Living Atlas Agriculture of the USA group.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

A digital orthophoto is a georeferenced image prepared from aerial imagery, or other remotely-sensed data in which the displacement within the image due to sensor orientation and terrain relief has been removed. Orthophotos combine the characteristics of an image with the geometric qualities of a map. Orthoimages show ground features such as roads, buildings, and streams in their proper positions, without the distortion characteristic of unrectified aerial imagery. Digital orthoimages produced and used within the Forest Service are developed from imagery acquired through various national and regional image acquisition programs. The resulting orthoimages, also known as orthomaps, can be directly applied in remote sensing, GIS and mapping applications. They serve a variety of purposes, from interim maps to references for earth science investigations and analysis. Because of the orthographic property, an orthoimage can be used like a map for measurement of distances, angles, and areas with scale being constant everywhere. Also, they can be used as map layers in GIS or other computer-based manipulation, overlaying, and analysis. An orthoimage differs from a map in a manner of depiction of detail; on a map only selected detail is shown by conventional symbols, whereas on an orthoimage all details appear just as in original aerial or satellite imagery.This record was taken from the USDA Enterprise Data Inventory that feeds into the https://data.gov catalog. Data for this record includes the following resources: ISO-19139 metadata ArcGIS Hub Dataset ArcGIS GeoServiceFor complete information, please visit https://data.gov.

This layer displays coal fields of the Conterminous United States from the U.S. Coal Resources and Reserves Assessment Project by the U.S. Geological Survey (USGS) Energy Resources Program.Dataset SummaryPhenomenon Mapped: Coal Fields of the United StatesCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Conterminous United StatesVisible Scale: All ScalesSource: National Coal Resource AssessmentPublication Date: 2013East, J.A., 2013, Coal fields of the conterminous United States—National Coal Resource Assessment updated version: U.S. Geological Survey Open-File Report 2012–1205, one sheet, scale 1:5,000,000, available at https://pubs.usgs.gov/of/2012/1205/. What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

The nine-banded Armadillo (Dasypus novemcinctus) is the only species of Armadillo in the United States and alters ecosystems by excavating extensive burrows used by many other wildlife species. Relatively little is known about its habitat use or population densities, particularly in developed areas, which may be key to facilitating its range expansion. We evaluated Armadillo occupancy and density in relation to anthropogenic and landcover variables in the Ozark Mountains of Arkansas along an urban to rural gradient. Armadillo detection probability was best predicted by temperature (positively) and precipitation (negatively). Contrary to expectations, occupancy probability of Armadillos was best predicted by slope (negatively) and elevation (positively) rather than any landcover or anthropogenic variables. Armadillo density varied considerably between sites (ranging from a mean of 4.88 – 46.20 Armadillos per km2) but was not associated with any environmental or anthropogenic variables. Methods Site Selection Our study took place in Northwest Arkansas, USA, in the greater Fayetteville metropolitan area. We deployed trail cameras (Spypoint Force Dark (Spypoint Inc, Victoriaville, Quebec, Canada) and Browning Strikeforce XD cameras (Browning, Morgan, Utah, USA) over the course of two winter seasons, December 2020-March 2021, and November 2021-March 2022. We sampled 10 study sites in year one, and 12 study sites in year two. All study sites were located in the Ozark Mountains ecoregion in Northwest Arkansas. Sites were all Oak Hickory dominated hardwood forests at similar elevation (213.6 – 541 m). Devils Eyebrow and ONSC are public natural areas managed by the Arkansas Natural heritage Commission (ANHC). Devil’s Den and Hobbs are managed by the Arkansas state park system. Markham Woods (Markham), Ninestone Land Trust (Ninestone) and Forbes, are all privately owned, though Markham has a publicly accessible trail system throughout the property. Lake Sequoyah, Mt. Sequoyah Woods, Kessler Mountain, Lake Fayetteville, and Millsaps Mountain are all city parks and managed by the city of Fayetteville. Lastly, both Weddington and White Rock are natural areas within Ozark National Forest and managed by the U.S. Forest Service. We sampled 5 sites in both years of the study including Devils Eyebrow, Markham Hill, Sequoyah Woods, Ozark Natural Science Center (ONSC), and Kessler Mountain. We chose our study sites to represent a gradient of human development, based primarily on Anthropogenic noise values (Buxton et al. 2017, Mennitt and Fristrup 2016). We chose open spaces that were large enough to accommodate camera trap research, as well as representing an array of anthropogenic noise values. Since anthropogenic noise is able to permeate into natural areas within the urban interface, introducing human disturbance that may not be detected by other layers such as impervious surface and housing unit density (Buxton et al. 2017), we used dB values for each site as an indicator of the level of urbanization. Camera Placement We sampled ten study sites in the first winter of the study. At each of the 10 study sites, we deployed anywhere between 5 and 15 cameras. Larger study areas received more cameras than smaller sites because all cameras were deployed a minimum of 150m between one another. We avoided placing cameras on roads, trails, and water sources to artificially bias wildlife detections. We also avoided placing cameras within 15m of trails to avoid detecting humans. At each of the 12 study areas we surveyed in the second winter season, we deployed 12 to 30 cameras. At each study site, we used ArcGIS Pro (Esri Inc, Redlands, CA) to delineate the trail systems and then created a 150m buffer on each side of the trail. We then created random points within these buffered areas to decide where to deploy cameras. Each random point had to occur within the buffered areas and be a minimum of 150m from the next nearest camera point, thus the number of cameras at each site varied based upon site size. We placed all cameras within 50m of the random points to ensure that cameras were deployed on safe topography and with a clear field of view, though cameras were not set in locations that would have increased animal detections (game trails, water sources, burrows etc.). Cameras were rotated between sites after 5 or 10 week intervals to allow us to maximize camera locations with a limited number of trail cameras available to us. Sites with more than 25 cameras were active for 5 consecutive weeks while sites with fewer than 25 cameras were active for 10 consecutive weeks. We placed all cameras on trees or tripods 50cm above ground and at least 15m from trails and roads. We set cameras to take a burst of three photos when triggered. We used Timelapse 2.0 software (Greenberg et al. 2019) to extract metadata (date and time) associated with all animal detections. We manually identified all species occurring in photographs and counted the number of individuals present. Because density estimation requires the calculation of detection rates (number of Armadillo detections divided by the total sampling period), we wanted to reduce double counting individuals. Therefore, we grouped photographs of Armadillos into “episodes” of 5 minutes in length to reduce double counting individuals that repeatedly triggered cameras (DeGregorio et al. 2021, Meek et al. 2014). A 5 min threshold is relatively conservative with evidence that even 1-minute episodes adequately reduces double counting (Meek et al. 2014). Landcover Covariates To evaluate occupancy and density of Armadillos based on environmental and anthropogenic variables, we used ArcGIS Pro to extract variables from 500m buffers placed around each camera (Table 2). This spatial scale has been shown to hold biological meaning for Armadillos and similarly sized species (DeGregorio et al. 2021, Fidino et al. 2016, Gallo et al. 2017, Magle et al. 2016). At each camera, we extracted elevation, slope, and aspect from the base ArcGIS Pro map. We extracted maximum housing unit density (HUD) using the SILVIS housing layer (Radeloff et al. 2018, Table 2). We extracted anthropogenic noise from the layer created by Mennitt and Fristrup (2016, Buxton et al. 2017, Table 2) and used the “L50” anthropogenic sound level estimate, which was calculated by taking the difference between predicted environmental noise and the calculated noise level. Therefore, we assume that higher levels of L50 sound corresponded to higher human presence and activity (i.e. voices, vehicles, and other sources of anthropogenic noise; Mennitt and Fristrup 2016). We derived the area of developed open landcover, forest area, and distance to forest edge from the 2019 National Land Cover Database (NLDC, Dewitz 2021, Table 2). Developed open landcover refers to open spaces with less than 20% impervious surface such as residential lawns, cemeteries, golf courses, and parks and has been shown to be important for medium-sized mammals (Gallo et al. 2017, Poessel et al. 2012). Forest area was calculated by combing all forest types within the NLCD layer (deciduous forest, mixed forest, coniferous forest), and summarizing the total area (km2) within the 500m buffer. Distance to forest edge was derived by creating a 30m buffer on each side of all forest boundaries and calculating the distance from each camera to the nearest forest edge. We calculated distance to water by combining the waterbody and flowline features in the National Hydrogeography Dataset (U.S. Geological Survey) for the state of Arkansas to capture both permanent and ephemeral water sources that may be important to wildlife. We measured the distance to water and distance to forest edge using the geoprocessing tool “near” in ArcGIS Pro which calculates the Euclidean distance between a point and the nearest feature. We extracted Average Daily Traffic (ADT) from the Arkansas Department of Transportation database (Arkansas GIS Office). The maximum value for ADT was calculated using the Summarize Within tool in ArcGIS Pro. We tested for correlation between all covariates using a Spearman correlation matrix and removed any variable with correlation greater than 0.6. Pairwise comparisons between distance to roads and HUD and between distance to forest edge and forest area were both correlated above 0.6; therefore, we dropped distance to roads and distance to forest edge from analyses as we predicted that HUD and forest area would have larger biological impacts on our focal species (Kretser et al. 2008). Occupancy Analysis In order to better understand habitat associations while accounting for imperfect detection of Armadillos, we used occupancy modeling (Mackenzie et al. 2002). We used a single-species, single-season occupancy model (Mackenzie et al. 2002) even though we had two years of survey data at 5 of the study sites. We chose to do this rather than using a multi-season dynamic occupancy model because most sites were not sampled during both years of the study. Even for sites that were sampled in both years, cameras were not placed in the same locations each year. We therefore combined all sampling into one single-season model and created unique site by year combinations as our sampling locations and we used year as a covariate for analysis to explore changes in occupancy associated with the year of study. For each sampling location, we created a detection history with 7 day sampling periods, allowing presence/absence data to be recorded at each site for each week of the study. This allowed for 16 survey periods between 01 December 2020, and 11 March 2021 and 22 survey periods between 01 November 2021 and 24 March 2022. We treated each camera as a unique survey site, resulting in a total of 352 sites. Because not all cameras were deployed at the same time and for the same length of time, we used a staggered entry approach. We used a multi-stage fitting approach in which we

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes cattle production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Cattle ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United States, Alaska, and HawaiiVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Cattle - Operations with SalesCattle - Sales in US DollarsCattle - Sales in HeadDairy - Operations with SalesDairy - Sales in US DollarsAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes winter wheat production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Winter Wheat ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United StatesVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Area Harvested in AcresOperations with Area HarvestedProduction in BushelsIrrigated Area Harvested in AcresAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.Additional information on wheat from the Census of Agriculture is available in the USDA Census of Agriculture 2017 - Wheat Production layer.Many other ready-to-use layers derived from the Census of Agriculture can be found in the Living Atlas Agriculture of the USA group.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes wheat production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Wheat ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United StatesVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Area Harvested in AcresOperations with Area HarvestedOperations with SalesProduction in BushelsSales in US DollarsIrrigated Area Harvested in AcresOperations with Irrigated Area HarvestedAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.Many other ready-to-use layers derived from the Census of Agriculture can be found in the Living Atlas Agriculture of the USA group.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.

The Census of Agriculture, produced by the USDA National Agricultural Statistics Service (USDA), provides a complete count of America's farms, ranches and the people who grow our food. The census is conducted every five years, most recently in 2017, and provides an in-depth look at the agricultural industry.This layer summarizes dairy production from the 2017 Census of Agriculture at the county level.This layer was produced from data downloaded using the USDA's QuickStats Application. The data was transformed using the Pivot Table tool in ArcGIS Pro and joined to the county boundary file provided by the USDA. The layer was published as feature layer in ArcGIS Online. Dataset SummaryPhenomenon Mapped: 2017 Dairy ProductionCoordinate System: Web Mercator Auxiliary SphereExtent: 48 Contiguous United States, Alaska, and HawaiiVisible Scale: All ScalesSource: USDA National Agricultural Statistics Service QuickStats ApplicationPublication Date: 2017AttributesThis layer provides values for the following attributes. Note that some values are not disclosed (coded as -1 in the layer) to protect the privacy of producers in areas with limited production.Cattle - Operations with SalesCattle - Sales in US DollarsCattle - Sales in HeadDairy - Operations with SalesDairy - Sales in US DollarsAdditionally attributes of State Name, State Code, County Name and County Code are included to facilitate cartography and use with other layers.What can you do with this layer?This layer can be used throughout the ArcGIS system. Feature layers can be used just like any other vector layer. You can use feature layers as an input to geoprocessing tools in ArcGIS Pro or in Analysis in ArcGIS Online. Combine the layer with others in a map and set custom symbology or create a pop-up tailored for your users. For the details of working with feature layers the help documentation for ArcGIS Pro or the help documentation for ArcGIS Online are great places to start. The ArcGIS Blog is a great source of ideas for things you can do with feature layers. This layer is part of ArcGIS Living Atlas of the World that provides an easy way to find and explore many other beautiful and authoritative layers, maps, and applications on hundreds of topics.