This dynamic World Elevation Terrain layer returns float values representing ground heights in meters and compiles multi-resolution data from many authoritative data providers from across the globe. Heights are orthometric (sea level = 0), and water bodies that are above sea level have approximated nominal water heights.Height units: MetersUpdate Frequency: QuarterlyCoverage: World/GlobalData Sources: This layer is compiled from a variety of best available sources from several data providers. To see the coverage and extents of various datasets comprising this service in an interactive map, see World Elevation Coverage Map.What can you do with this layer?Use for Visualization: This layer is generally not optimal for direct visualization. By default, 32 bit floating point values are returned, resulting in higher bandwidth requirements. Therefore, usage should be limited to applications requiring elevation data values. Alternatively, client applications can select from numerous additional functions, applied on the server, that return rendered data. For visualizations such as multi-directional hillshade, hillshade, elevation tinted hillshade, and slope, consider using the appropriate server-side function defined on this service.Use for Analysis: Yes. This layer provides data as floating point elevation values suitable for use in analysis. There is a limit of 5000 rows x 5000 columns.Note: This layer combine data from different sources and resamples the data dynamically to the requested projection, extent and pixel size. For analyses using ArcGIS Desktop, it is recommended to filter a dataset, specify the projection, extent and cell size using the Make Image Server Layer geoprocessing tool. The extent is factor of cell size and rows/columns limit. e.g. if cell size is 10 m, the extent for analysis would be less than 50,000 m x 50,000 m.Server Functions: This layer has server functions defined for the following elevation derivatives. In ArcGIS Pro, server function can be invoked from Layer Properties - Processing Templates.

Slope Degrees Slope Percent Aspect Ellipsoidal height Hillshade Multi-Directional Hillshade Dark Multi-Directional Hillshade Elevation Tinted Hillshade Slope Map Aspect Map Mosaic Method: This image service uses a default mosaic method of "By Attribute”, using Field 'Best' and target of 0. Each of the rasters has been attributed with ‘Best’ field value that is generally a function of the pixel size such that higher resolution datasets are displayed at higher priority. Other mosaic methods can be set, but care should be taken as the order of the rasters may change. Where required, queries can also be set to display only specific datasets such as only NED or the lock raster mosaic rule used to lock to a specific dataset.Accuracy: Accuracy will vary as a function of location and data source. Please refer to the metadata available in the layer, and follow the links to the original sources for further details. An estimate of CE90 and LE90 are included as attributes, where available.This layer allows query, identify, and export image requests. The layer is restricted to a 5,000 x 5,000 pixel limit in a single request.This layer is part of a larger collection of elevation layers that you can use to perform a variety of mapping analysis tasks.

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

Maps of California's Wildland Urban Interface (WUI) generated using the Time Step Moving Window (TSMW) method outlined in the paper "Remapping California's Wildland Urban Interface: A Property-Level Time-Space Framework, 2000-2020".

Please cite the original paper:

Berg, Aleksander K, Dylan S. Connor, Peter Kedron, and Amy E. Frazier. 2024. “Remapping California’s Wildland Urban Interface: A Property-Level Time-Space Framework, 2000–2020.” Applied Geography 167 (June): 103271. https://doi.org/10.1016/j.apgeog.2024.103271.

WUI maps were generated using Zillow ZTRAX parcel level attributes joined with FEMA USA Structures building footprints and the National Land Cover Database (NLCD).

All files are geotiff rasters with WUI areas mapped at a ~30m resolution. A raster value of null indicates not WUI, raster value of 1 indicates intermix WUI, and a raster value of 2 indicates interface WUI.

Three WUI maps were generated using structures built on of before the years indicated below:

2000 - "CA_WUI_2000.tif"

2010 - "CA_WUI_2010.tif"

2020 - "CA_WUI_2020.tif"

Acknowledgments -

We thank our reviewers and editors for helping us to improve the manuscript. We gratefully acknowledge access to the Zillow Transaction and Assessment Dataset (ZTRAX) through a data use agreement between the University of Colorado Boulder, Arizona State University, and Zillow Group, Inc. More information on accessing the data can be found at http://www.zillow.com/ztrax. The results and opinions are those of the author(s) and do not reflect the position of Zillow Group. Support by Zillow Group Inc. is acknowledged. We thank Johannes Uhl and Stefan Leyk for their great work in preparing the original dataset. For feedback and comments, we also thank Billie Lee Turner II, Sharmistha Bagchi-Sen, and participants at the 2022 Global Conference on Economic Geography, the 2022 Young Economic Geographers Network meeting, and the 2023 annual meeting of the American Association of Geographers. Funding for our work has been provided by Arizona State University's Institute of Social Science Research (ISSR) Seed Grant Initiative. Additional funding was provided through the Humans, Disasters, and the Built Environment program of the National Science Foundation, Award Number 1924670 to the University of Colorado Boulder, the Institute of Behavioral Science, Earth Lab, the Cooperative Institute for Research in Environmental Sciences, the Grand Challenge Initiative and the Innovative Seed Grant program at the University of Colorado Boulder as well as the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health under Award Numbers R21 HD098717 01A1 and P2CHD066613.

The Raster Based GIS Coverage of Mexican Population is a gridded coverage (1 x 1 km) of Mexican population. The data were converted from vector into raster. The population figures were derived based on available point data (the population of known localities - 30,000 in all). Cell values were derived using a weighted moving average function (Burrough, 1986), and then calculated based on known population by state. The result from this conversion is a coverage whose population data is based on square grid cells rather than a series of vectors. This data set is produced by the Columbia University Center for International Earth Science Information Network (CIESIN) in collaboration with the Instituto Nacional de Estadistica Geografia e Informatica (INEGI).

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

The Terrain Ruggedness Index (TRI) is used to express the amount of elevation difference between adjacent cells of a DEM. This raster function template is used to generate a visual representation of the TRI with your elevation data. The results are interpreted as follows:0-80m is considered to represent a level terrain surface81-116m represents a nearly level surface117-161m represents a slightly rugged surface162-239m represents an intermediately rugged surface240-497m represents a moderately rugged surface498-958m represents a highly rugged surface959-4367m represents an extremely rugged surfaceWhen to use this raster function templateThe main value of this measurement is that it gives a relatively accurate view of the vertical change taking place in the terrain model from cell to cell. The TRI provides data on the relative change in height of the hillslope (rise), such as the side of a canyon.How to use this raster function templateIn ArcGIS Pro, search ArcGIS Living Atlas for raster function templates to apply them to your imagery layer. You can also download the raster function template, attach it to a mosaic dataset, and publish it as an image service. The output is a visual TRI representation of your imagery. This index supports elevation data.References:Raster functionsApplicable geographiesThe index is a standard index which is designed to work globally.

This map shows the global vulnerability with respect to natural hazards. To derive the vulnerability using the above described resilience the population density is used as the main impact factor for the vulnerability. So the more people live at a certain place, the higher the vulnerability.To calculate the vulnerability the value gained from the population reclassification is divided by the value from the global resilience map.The population density is evaluated by the following reclass table:• < 0.1 person/km²: not being used• 0.1 - 10 person/km²: 1• 10 - 100 person/km²: 2• 100 - 1 000 person/km²: 3• 1 000 - 10 000 person/km²: 4• > 10 000 person/km²: 5The resulting raster values have a range between 0 (low vulnerability) and 4 (high vulnerability).DataPopulation density © Socioeconomic Data and Applications Center (SEDAC)Airports © ourairports.comWorld Port Index © msi.nga.milTerrain Data (250 grid resolution) © SRTMGlobal Needs Assessment - Vulnerability Index and Crisis Index 2013 © European Commission ECHOMedical Care © www.laenderdaten.deIncome © World BankHuman rights and state stability © The Global Slavery Index 2013Education © World BankBackgroundBasemap © Esri, DeLorme, GEBCO, NOAA NGDC, and other contributors

From the site: "This raster dataset has been created using the original data "Pennsylvania conservation gap fish habitat models" as originated by the Environmental Resources Research Institute of The Pennsylvania State University. Conservation Values were then assigned to species as determined by SmartConservation® methodology and combined to create an overall conservation value raster for fish. The resulting raster was then reclassified into 10 quantiles as follows: Old Value New Value 0 0 1-59 1 60-67 2 68-81 3 82-92 4 93 5 94-122 6 123-125 7 126 8 127-177 9 178-202 10 Conservation values were determined by experts gathered by Natural Lands Trust through SmartConservation®. This data set is one of several that have been combined to create an overall aquatic resources conservation value raster for the Central Appalachian Forest Ecoregion. Therefore the values were determined as a relative rank, comparable in value only to the other input aquatic resources data. Conservation value ranges from 1 - 10 with 10 being the highest value."

This dataset represents a unique compiled environmental data set for the circumpolar Arctic ocean region 45N to 90N region. It consists of 170 layers (mostly marine, some terrestrial) in ArcGIS 10 format to be used with a Geographic Information System (GIS) and which are listed below in detail. Most layers are long-term average raster GRIDs for the summer season, often by ocean depth, and represent value-added products easy to use. The sources of the data are manifold such as the World Ocean Atlas 2009 (WOA09), International Bathimetric Chart of the Arctic Ocean (IBCAO), Canadian Earth System Model 2 (CanESM2) data (the newest generation of models available) and data sources such as plankton databases and OBIS. Ocean layers were modeled and predicted into the future and zooplankton species were modeled based on future data: Calanus hyperboreus (AphiaID104467), Metridia longa (AphiaID 104632), M. pacifica (AphiaID 196784) and Thysanoessa raschii (AphiaID 110711). Some layers are derived within ArcGIS. Layers have pixel sizes between 1215.819573 meters and 25257.72929 meters for the best pooled model, and between 224881.2644 and 672240.4095 meters for future climate data. Data was then reprojected into North Pole Stereographic projection in meters (WGS84 as the geographic datum). Also, future layers are included as a selected subset of proposed future climate layers from the Canadian CanESM2 for the next 100 years (scenario runs rcp26 and rcp85). The following layer groups are available: bathymetry (depth, derived slope and aspect); proximity layers (to,glaciers,sea ice, protected areas, wetlands, shelf edge); dissolved oxygen, apparent oxygen, percent oxygen, nitrogen, phosphate, salinity, silicate (all for August and for 9 depth classes); runoff (proximity, annual and August); sea surface temperature; waterbody temperature (12 depth classes); modeled ocean boundary layers (H1, H2, H3 and Wx).This dataset is used for a M.Sc. thesis by the author, and freely available upon request. For questions and details we suggest contacting the authors. Process_Description: Please contact Moritz Schmid for the thesis and detailed explanations. Short version: We model predicted here for the first time ocean layers in the Arctic Ocean based on a unique dataset of physical oceanography. Moreover, we developed presence/random absence models that indicate where the studied zooplankton species are most likely to be present in the Arctic Ocean. Apart from that, we develop the first spatially explicit models known to science that describe the depth in which the studied zooplankton species are most likely to be at, as well as their distribution of life stages. We do not only do this for one present day scenario. We modeled five different scenarios and for future climate data. First, we model predicted ocean layers using the most up to date data from various open access sources, referred here as best-pooled model data. We decided to model this set of stratification layers after discussions and input of expert knowledge by Professor Igor Polyakov from the International Arctic Research Center at the University of Alaska Fairbanks. We predicted those stratification layers because those are the boundaries and layers that the plankton has to cross for diel vertical migration and a change in those would most likely affect the migration. I assigned 4 variables to the stratification layers. H1, H2, H3 and Wx. H1 is the lower boundary of the mixed layer depth. Above this layer a lot of atmospheric disturbance is causing mixing of the water, giving the mixed layer its name. H2, the middle of the halocline is important because in this part of the ocean a strong gradient in salinity and temperature separates water layers. H3, the isotherm is important, because beneath it flows denser and colder Atlantic water. Wx summarizes the overall width of the described water column. Ocean layers were predicted using machine learning algorithms (TreeNet, Salford Systems). Second, ocean layers were included as predictors and used to predict the presence/random absence, most likely depth and life stage layers for the zooplankton species: Calanus hyperboreus, Metridia longa, Metridia pacifica and Thysanoessa raschii, This process was repeated for future predictions based on the CanESM2 data (see in the data section). For zooplankton species the following layers were developed and for the future. C. hyperboreus: Best-pooled model as well as future predictions (rcp26 including ocean layer(also excluding), rcp85 including oocean layers (also excluding) for 2010 and 2100.For parameters: Presence/random absence, most likely depth and life stage layers M. longa: Best-pooled model as well as future predictions (rcp26 including ocean layer(also excluding), rcp85 including oocean layers (also excluding) for 2010 and 2100. For parameters: Presence/rand... Visit https://dataone.org/datasets/f63d0f6c-7d53-46ce-b755-42a368007601 for complete metadata about this dataset.

This raster GIS dataset contains 5-meter-resolution cells depicting the areas of total marsh gain (value=1), lost (value=-1) and remaining (no change; value=0). Total marsh (TM) was defined as the sum of low marsh and high marsh [SLAMM category 8 + SLAMM category 7 + SLAMM category 20]. Based on SLAMM simulation outputs, we generated the gain and loss map by using the “Raster Calculator” tool under “Spatial Analyst Tools” in ArcGIS software. The methodology consists of the three steps listed below (where we use low marsh [LM] as an example). The same process can be applied to other SLAMM land cover categories.

1) Open ArcMap, add SLAMM simulation raster outputs (all SLAMM categories) for baseline year and future years.

2) In Raster Calculator, set the SLAMM codeequal to8 (low marsh = SLAMM category 8) to generate a new raster. Each individual cell in the new raster is assigned a value of “0” or “1”. “1” is low marsh and “0” is any other SLAMM land cover category. Perform this step for both the baseline year and future year.

3) In Raster Calculator, subtract the new raster for the baseline year from the new raster for the future year (formula = new future year raster - new baseline year raster). The calculation generates a new raster, in which each individual cell is assigned a value of “-1”, “0”, or “1”. Based on the calculation, “-1” means low marsh loss in the future (the cell has converted from low marsh to a different SLAMM category), “0” means low marsh is remaining (the cell stays the same), and “1” means low marsh gain in the future (the cell has converted from a different SLAMM category to low marsh).

Prior SLAMM work has been performed in the Delaware Bay, but our methods differ in that we derive results for specific marsh areas and utilize more recent, higher resolution elevation data (2015 USGS CoNED Topobathy Model: New Jersey and Delaware), the most recent SLR projections, and site-specific accretion data (through 2016). These SLAMM simulations were performed as part of a larger project by the USEPA on frameworks and methods for characterizing relative wetland vulnerabilities.

Note: additional raster files from this project are available upon request. These include files from low and high SLR scenarios and different model protection scenarios. For more information, contact Jordan West (West.Jordan@epa.gov).

Impervious surfaces are surfaces that do not allow water to pass through. Examples of these surfaces include highways, parking lots, rooftops, and airport runways. Instead of allowing rain to pass into the soil, impervious surfaces cause water to collect at the surface, then run off. An increase in impervious surface area causes an increase of water volume which needs to be managed by stormwater systems. With the flow come pollutants, which collect on impervious surfaces then discharge with the runoff into streams and the ocean. Runoff water does not enter the water table, and that can cause other management issues, such as interruptions in baseline stream flow.The NLCD imperviousness layer represents urban impervious surfaces as a percentage of developed surface over every 30-meter pixel in the United States. Phenomenon Mapped: The proportion of the landscape that is impervious to water.Time Extent: 2001, 2004, 2006, 2008, 2011, 2013, 2016, 2019, and 2021 for the lower 48 conterminous US states. A small portion of Alaska around Anchorage displays a time series of 2001, 2011, and 2016. Hawaii, Puerto Rico, and the US Virgin Islands unfortunately only have data for 2001 so there is only one image in the series. This information may be used in conjunction with the USA NLCD Land Cover layer.Units: PercentCell Size: 30 metersSource Type: DiscretePixel Type: Unsigned integerData Coordinate System: North America Albers Equal Area Conic (102008)Mosaic Projection: North America Albers Equal Area Conic (102008)Extent: CONUS, Hawaii, A portion of Alaska around Anchorage, District of Columbia, Puerto RicoNoData Value: 127Source: Multi-Resolution Land Characteristics ConsortiumPublication Date: June 30, 2023ArcGIS Server URL: https://landscape10.arcgis.com/arcgis/Time SeriesBy default, this layer will appear in your client with a time slider which allows you to play the series as an animation. The animation will advance year by year, but the layer only changes appearance every few years in the lower 48 states, in 2001, 2004, 2006, 2008, 2011, 2013, 2016, 2019, and 2021. To select just one year in the series, first turn the time series off on the time slider, then create a definition query on the layer which selects only the desired year.Time Series DescriptorMRLC issued a set of companion rasters with this impervious surface layer showing the reason why each pixel is impervious. This companion layer, called the Developed Imperviousness Descriptor, is not currently available in this map service. The descriptor layer identifies types of roads, core urban areas, and energy production sites for each impervious pixel to allow deeper analysis of developed features. The descriptor layer may be downloaded directly from MRLC and added to ArcGIS Pro.Alaska, Hawaii, and Puerto RicoAt this time Alaska, Hawaii, and Puerto Rico are produced with a different methodology, and are not set up to be directly compared the way the CONUS time series is. To analyze change between the latest two data years for this portion of Alaska, be sure to use the NLCD 2011 to 2016 Developed Impervious Change raster. For Hawaii and Puerto Rico, only the year 2001 is available for download at the MRLC.North America Albers ProjectionAll NLCD layers in the Living Atlas are projected into the North America Albers Projection before serving in the Living Atlas. This allows the coterminous USA, Puerto Rico, Hawaii, and Alaska to be served from a common projection and analyzed together. In tests performed by esri, the NLCD land cover classes after projection to North America Albers had the exact same number of pixels in input as output, but pixels had been slightly rearranged after projection. Processing TemplatesThis layer comes with two color schemes, cool and warm. The default is a cool gray color scheme, designed to look good on light and dark gray web maps. To choose a warm color scheme which was the default until 2021, change the processing template to the Impervious Surface Warm Renderer in your map client.Dataset SummaryThe National Land Cover Database products are created through a cooperative project conducted by the Multi-Resolution Land Characteristics Consortium (MRLC). The MRLC Consortium is a partnership of federal agencies, consisting of the U.S. Geological Survey, the National Oceanic and Atmospheric Administration, the U.S. Environmental Protection Agency, the U.S. Department of Agriculture, the U.S. Forest Service, the National Park Service, the U.S. Fish and Wildlife Service, the Bureau of Land Management and the USDA Natural Resources Conservation Service.What can you do with this layer?This layer can be used to create maps and to visualize the underlying data. This layer can be used as an analytic input in ArcGIS Desktop.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.

From the site: "This raster dataset has been created using the "Floodplains from the PA Explorer CD-ROM edition" for each county in the study area as originated by the Office of Remote Sensing for Earth Resources, Penn State University (see metadata entitled "Bucks FP.doc"). All areas designated in the shapefile were assigned a conservation value of 5. Conservation values were determined by experts gathered by Natural Lands Trust through SmartConservation®. This data set is one of several that have been combined to create an overall aquatic resources conservation value raster for the expanded piedmont ecoregion. Therefore the values were determined as a relative rank, comparable in value only to the other input aquatic resources data. Conservation value ranges from 1 - 10 with 10 being the highest value."

This raster represents a continuous surface of sage-grouse habitat suitability index (HSI, created using ArcGIS 10.2.2) values for Nevada during spring, which is a surrogate for habitat conditions during the sage-grouse breeding and nesting period. Summary of steps to create Habitat Categories: HABITAT SUITABILITY INDEX: The HSI was derived from a generalized linear mixed model (specified by binomial distribution) that contrasted data from multiple environmental factors at used sites (telemetry locations) and available sites (random locations). Predictor variables for the model represented vegetation communities at multiple spatial scales, water resources, habitat configuration, urbanization, roads, elevation, ruggedness, and slope. Vegetation data was derived from various mapping products, which included NV SynthMap (Petersen 2008, SageStitch (Comer et al. 2002, LANDFIRE (Landfire 2010), and the CA Fire and Resource Assessment Program (CFRAP 2006). The analysis was updated to include high resolution percent cover within 30 x 30 m pixels for Sagebrush, non-sagebrush, herbaceous vegetation, and bare ground (C. Homer, unpublished; based on the methods of Homer et al. 2014, Xian et al. 2015 ) and conifer (primarily pinyon-juniper, P. Coates, unpublished). The pool of telemetry data included the same data from 1998 - 2013 used by Coates et al. (2014); additional telemetry location data from field sites in 2014 were added to the dataset. The dataset was then split according calendar date into three seasons (spring, summer, winter). Summer included telemetry locations (n = 14,058) from mid-March to June. All age and sex classes of marked grouse were used in the analysis. Sufficient data (i.e., a minimum of 100 locations from at least 20 marked Sage-grouse) for modeling existed in 10 subregions for spring and summer, and seven subregions in winter, using all age and sex classes of marked grouse. It is important to note that although this map is composed of HSI values derived from the seasonal data, it does not explicitly represent habitat suitability for reproductive females (i.e., nesting and with broods). Insufficient data were available to allow for estimation of this habitat type for all seasons throughout the study area extent. A Resource Selection Function (RSF) was calculated using R Software (v 3.13) for each subregion and using generalized linear models to derive model-averaged parameter estimates for each covariate across a set of additive models. Subregional RSFs were transformed into Habitat Suitability Indices, and averaged together to produce an overall statewide HSI whereby a relative probability of occurrence was calculated for each raster cell during the spring. In order to account for discrepancies in HSI values caused by varying ecoregions within Nevada, the HSI was divided into north and south extents using a slightly modified flood region boundary (Mason 1999) that was designed to represent respective mesic and xeric regions of the state. North and south HSI rasters were each relativized according to their maximum value to rescale between zero and one, then mosaicked once more into a state-wide extent. REFERENCES: California Forest and Resource Assessment Program (CFRAP). 2006. Statewide Land Use / Land Cover Mosaic. [Geospatial data.] California Department of Forestry and Fire Protection, http://frap.cdf.ca.gov/data/frapgisdata-sw-rangeland-assessment_data.php Census 2010. TIGER/Line Shapefiles. Urban Areas [Geospatial data.] U.S. Census Bureau, Washington D.C., https://www.census.gov/geo/maps-data/data/tiger-line.html Census 2014. TIGER/Line Shapefiles. Roads [Geospatial data.] U.S. Census Bureau, Washington D.C., https://www.census.gov/geo/maps-data/data/tiger-line.html Census 2015. TIGER/Line Shapefiles. Blocks [Geospatial data.] U.S. Census Bureau, Washington D.C., https://www.census.gov/geo/maps-data/data/tiger-line.html Coates, P.S., Casazza, M.L., Brussee, B.E., Ricca, M.A., Gustafson, K.B., Overton, C.T., Sanchez-Chopitea, E., Kroger, T., Mauch, K., Niell, L., Howe, K., Gardner, S., Espinosa, S., and Delehanty, D.J. 2014, Spatially explicit modeling of greater sage-grouse (Centrocercus urophasianus) habitat in Nevada and northeastern California—A decision-support tool for management: U.S. Geological Survey Open-File Report 2014-1163, 83 p., http://dx.doi.org/10.3133/ofr20141163. ISSN 2331-1258 (online) Comer, P., Kagen, J., Heiner, M., and Tobalske, C. 2002. Current distribution of sagebrush and associated vegetation in the western United States (excluding NM). [Geospatial data.] Interagency Sagebrush Working Group, http://sagemap.wr.usgs.gov Homer, C.G., Aldridge, C.L., Meyer, D.K., and Schell, S.J. 2014. Multi-Scale Remote Sensing Sagebrush Characterization with Regression Trees over Wyoming, USA; Laying a Foundation for Monitoring. International Journal of Applied Earth Observation and Geoinformation 14, Elsevier, US. LANDFIRE. 2010. 1.2.0 Existing Vegetation Type Layer. [Geospatial data.] U.S. Department of the Interior, Geological Survey, http://landfire.cr.usgs.gov/viewer/ Mason, R.R. 1999. The National Flood-Frequency Program—Methods For Estimating Flood Magnitude And Frequency In Rural Areas In Nevada U.S. Geological Survey Fact Sheet 123-98 September, 1999, Prepared by Robert R. Mason, Jr. and Kernell G. Ries III, of the U.S. Geological Survey; and Jeffrey N. King and Wilbert O. Thomas, Jr., of Michael Baker, Jr., Inc. http://pubs.usgs.gov/fs/fs-123-98/ Peterson, E. B. 2008. A Synthesis of Vegetation Maps for Nevada (Initiating a 'Living' Vegetation Map). Documentation and geospatial data, Nevada Natural Heritage Program, Carson City, Nevada, http://www.heritage.nv.gov/gis Xian, G., Homer, C., Rigge, M., Shi, H., and Meyer, D. 2015. Characterization of shrubland ecosystem components as continuous fields in the northwest United States. Remote Sensing of Environment 168:286-300. NOTE: This file does not include habitat areas for the Bi-State management area and the spatial extent is modified in comparison to Coates et al. 2014

This data set represents a 2-meter resolution LiDAR first return surface or Digital Surface Model (DSM) for New Hampshire. It was derived from a statewide LAS Dataset which comprised 8 separate LiDAR collections that covered the state as of January, 2020. The LAS Dataset was used as input to the ArcGIS "LAS Dataset to Raster" geoprocessing tool which converted the LAS first return values to raster values in the output data set. In some areas, users may notice unusual linear edges which appear unlikely or anomalous. The LiDAR vendor explained that these anomalies may be the result of changes in the degrees of tree canopy closure that occurred between the times adjacent flight lines were completed. Although leaf-off conditions were specified for data collection, strict adherence to the project specifications was not possible in all locations and exceptions occurred in order to complete data acquisition in a timely manner. As a result, abrupt discontinuities may be noticeable where data were collected on different dates. Eamples of these anomalies can be found in the areas of Cave Mountain in Bartlett and to the west of Woodstock.

The Federal Emergency Management Agency (FEMA) produces Flood Insurance Rate maps and identifies Special Flood Hazard Areas as part of the National Flood Insurance Program's floodplain management. Special Flood Hazard Areas have regulations that include the mandatory purchase of flood insurance for holders of federally regulated mortgages. In addition, this layer can help planners and firms avoid areas of flood risk and also avoid additional cost to carry insurance for certain planned activities. Dataset SummaryPhenomenon Mapped: Flood Hazard AreasGeographic Extent: Contiguous United States, Alaska, Hawaii, Puerto Rico, Guam, US Virgin Islands, Northern Mariana Islands and American Samoa.Projection: Web Mercator Auxiliary SphereData Coordinate System: USA Contiguous Albers Equal Area Conic USGS version (contiguous US, Puerto Rico, US Virgin Islands), WGS 1984 Albers (Alaska), Hawaii Albers Equal Area Conic (Hawaii), Western Pacific Albers Equal Area Conic (Guam, Northern Mariana Islands, and American Samoa)Cell Sizes: 10 meters (default), 30 meters, and 90 metersUnits: NoneSource Type: ThematicPixel Type: Unsigned integerSource: Federal Emergency Management Agency (FEMA)Update Frequency: AnnualPublication Date: May 7, 2025 This layer is derived from the May 7, 2025 version Flood Insurance Rate Map feature class S_FLD_HAZ_AR. The vector data were then flagged with an index of 94 classes, representing a unique combination of values displayed by three renderers. (In three resolutions the three renderers make nine processing templates.) Repair Geometry was run on the set of features, then the features were rasterized using the 94 class index at a resolutions of 10, 30, and 90 meters, using the Polygon to Raster tool and the "MAXIMUM_COMBINED_AREA" option. Not every part of the United States is covered by flood rate maps. This layer compiles all the flood insurance maps available at the time of publication. To make analysis easier, areas that were NOT mapped by FEMA for flood insurance rates no longer are served as NODATA but are filled in with a value of 250, representing any unmapped areas which appear in the US Census boundary of the USA states and territories. The attribute table corresponding to value 250 will indicate that the area was not mapped.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 "flood hazard areas" 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 "flood hazard areas" in the search box, browse to the layer then click OK. In ArcGIS Pro you can use the built-in raster functions to create custom extracts of the data. Imagery layers provide fast, powerful inputs to geoprocessing tools, models, or Python scripts in Pro. 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. Processing TemplatesCartographic Renderer - The default. These are meaningful classes grouped by FEMA which group its own Flood Zone Type and Subtype fields. This renderer uses FEMA's own cartographic interpretations of its flood zone and zone subtype fields to help you identify and assess risk. Flood Zone Type Renderer - Specifically renders FEMA FLD_ZONE (flood zone) attribute, which distinguishes the original, broadest categories of flood zones. This renderer displays high level categories of flood zones, and is less nuanced than the Cartographic Renderer. For example, a fld_zone value of X can either have moderate or low risk depending on location. This renderer will simply render fld_zone X as its own color without identifying "500 year" flood zones within that category.Flood Insurance Requirement Renderer - Shows Special Flood Hazard Area (SFHA) true-false status. This may be helpful if you want to show just the places where flood insurance is required. A value of True means flood insurance is mandatory in a majority of the area covered by each 10m pixel. Each of these three renderers have templates at three different raster resolutions depending on your analysis needs. To include the layer in web maps to serve maps and queries, the 10 meter renderers are the preferred option. These are served with overviews and render at all resolutions. However, when doing analysis of larger areas, we now offer two coarser resolutions of 30 and 90 meters in processing templates for added convenience and time savings.Data DictionaryMaking a copy of your area of interest using copyraster in arcgis pro will copy the layer's attribute table to your network alongside the local output raster. The raster attribute table in the copied raster will contain the flood zone, zone subtype, and special flood hazard area true/false flag which corresponds to each value in the layer for your area of interest. For your convienence, we also included a table in CSV format in the box below as a data dictionary you can use as an index to every value in the layer. Value,FLD_ZONE,ZONE_SUBTY,SFHA_TF 2,A,, 3,A,,F 4,A,,T 5,A,,T 6,A,,T 7,A,1 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN CHANNEL,T 8,A,1 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN STRUCTURE,T 9,A,ADMINISTRATIVE FLOODWAY,T 10,A,COASTAL FLOODPLAIN,T 11,A,FLOWAGE EASEMENT AREA,T 12,A99,,T 13,A99,AREA WITH REDUCED FLOOD RISK DUE TO LEVEE,T 14,AE,,F 15,AE,,T 16,AE,,T 17,AE,,T 18,AE,1 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN CHANNEL,T 19,AE,1 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN STRUCTURE,T 20,AE,"1 PCT CONTAINED IN STRUCTURE, COMMUNITY ENCROACHMENT",T 21,AE,"1 PCT CONTAINED IN STRUCTURE, FLOODWAY",T 22,AE,ADMINISTRATIVE FLOODWAY,T 23,AE,AREA OF SPECIAL CONSIDERATION,T 24,AE,COASTAL FLOODPLAIN,T 25,AE,COLORADO RIVER FLOODWAY,T 26,AE,COMBINED RIVERINE AND COASTAL FLOODPLAIN,T 27,AE,COMMUNITY ENCROACHMENT,T 28,AE,COMMUNITY ENCROACHMENT AREA,T 29,AE,DENSITY FRINGE AREA,T 30,AE,FLOODWAY,T 31,AE,FLOODWAY CONTAINED IN CHANNEL,T 32,AE,FLOODWAY CONTAINED IN STRUCTURE,T 33,AE,FLOWAGE EASEMENT AREA,T 34,AE,RIVERINE FLOODWAY IN COMBINED RIVERINE AND COASTAL ZONE,T 35,AE,RIVERINE FLOODWAY SHOWN IN COASTAL ZONE,T 36,AE,STATE ENCROACHMENT AREA,T 37,AH,,T 38,AH,,T 39,AH,FLOODWAY,T 40,AO,,T 41,AO,COASTAL FLOODPLAIN,T 42,AO,FLOODWAY,T 43,AREA NOT INCLUDED,,F 44,AREA NOT INCLUDED,,T 45,AREA NOT INCLUDED,,U 46,D,,F 47,D,,T 48,D,AREA WITH FLOOD RISK DUE TO LEVEE,F 49,OPEN WATER,,F 50,OPEN WATER,,T 51,OPEN WATER,,U 52,V,,T 53,V,COASTAL FLOODPLAIN,T 54,VE,,T 55,VE,,T 56,VE,COASTAL FLOODPLAIN,T 57,VE,RIVERINE FLOODWAY SHOWN IN COASTAL ZONE,T 58,X,,F 59,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD,F 60,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD,T 61,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD,U 62,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN CHANNEL,F 63,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD CONTAINED IN STRUCTURE,F 64,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD IN COASTAL ZONE,F 65,X,0.2 PCT ANNUAL CHANCE FLOOD HAZARD IN COMBINED RIVERINE AND COASTAL ZONE,F 66,X,"1 PCT CONTAINED IN STRUCTURE, COMMUNITY ENCROACHMENT",F 67,X,"1 PCT CONTAINED IN STRUCTURE, FLOODWAY",F 68,X,1 PCT DEPTH LESS THAN 1 FOOT,F 69,X,1 PCT DRAINAGE AREA LESS THAN 1 SQUARE MILE,F 70,X,1 PCT FUTURE CONDITIONS,F 71,X,1 PCT FUTURE CONDITIONS CONTAINED IN STRUCTURE,F 72,X,"1 PCT FUTURE CONDITIONS, COMMUNITY ENCROACHMENT",F 73,X,"1 PCT FUTURE CONDITIONS, FLOODWAY",F 74,X,"1 PCT FUTURE IN STRUCTURE, COMMUNITY ENCROACHMENT",F 75,X,"1 PCT FUTURE IN STRUCTURE, FLOODWAY",F 76,X,AREA OF MINIMAL FLOOD HAZARD, 77,X,AREA OF MINIMAL FLOOD HAZARD,F 78,X,AREA OF MINIMAL FLOOD HAZARD,T 79,X,AREA OF MINIMAL FLOOD HAZARD,U 80,X,AREA OF SPECIAL CONSIDERATION,F 81,X,AREA WITH REDUCED FLOOD RISK DUE TO LEVEE,F 82,X,AREA WITH REDUCED FLOOD RISK DUE TO LEVEE,T 83,X,FLOWAGE EASEMENT AREA,F 84,X,1 PCT FUTURE CONDITIONS,T 85,AH,COASTAL FLOODPLAIN,T 86,AE,,U 87,AE,FLOODWAY,F 88,X,AREA WITH REDUCED FLOOD HAZARD DUE TO ACCREDITED LEVEE SYSTEM,F 89,X,530,F 90,VE,100,T 91,AE,100,T 92,A99,AREA WITH REDUCED FLOOD HAZARD DUE TO LEVEE SYSTEM,T 93,A99,AREA WITH REDUCED FLOOD HAZARD DUE TO NON-ACCREDITED LEVEE SYSTEM,T 94,A,COMBINED RIVERINE AND COASTAL FLOODPLAIN,T 250,AREA NOT INCLUDED,Not Mapped by FEMA, Questions?Please leave a comment below if you have a question about this layer, and we will get back to you as soon as possible.

Supplementary materials used in the following studies:Kwong, I. H. Y., Lai, D. Y. F., Wong, F. K. K., & Fung, T. (Manuscript submitted for publication). Integrating five decades of Landsat imagery for territory-wide habitat mapping and change detection in a subtropical metropolitan city.Kwong, I. H. Y. (2025). Spatio-Temporal Changes in Habitat Type and Quality in Hong Kong Using a 50-Year Archive of Remote Sensing Imagery [Doctoral thesis, Department of Geography and Resource Management, The Chinese University of Hong Kong].Kwong, I. H. Y., Lai, D. Y. F., Wong, F. K. K., & Fung, T. (2025). Spatial variations in forest succession rates revealed from multi-temporal habitat maps using Landsat imagery in subtropical Hong Kong. European Geosciences Union (EGU) General Assembly 2025, Vienna, Austria, 27 Apr–2 May 2025. https://doi.org/10.5194/egusphere-egu25-2667. [Poster Presentation: https://presentations.copernicus.org/EGU25/EGU25-2667_presentation-h291057.pdf]Disclaimer: All datasets described here are for reference only. No express or implied warranty or representation is given to the accuracy or completeness of the data or its appropriateness for use in any particular circumstances.GIS mapping results:All raster layers (GeoTiff format) have a pixel size of 30 m covering the 1117-km2 terrestrial area in Hong Kong in this study (Hong Kong 1980 Grid coordinate system). The time period of 1973–2022 was divided into 10 five-year periods in the mapping process.HabitatMapHK_6class_yyyy-yyyy.tif: Raster data showing the 6 habitat classes mapped in this study. Pixel values range from 1 to 6 representing woodland, shrubland, grassland, barren land, built-up area, and water respectively.HabitatMapHK_EstimatedArea.csv: Area coverage (km2) of different habitat classes, as well as their confidence intervals, as mapped in this study.HabitatMapHK_6class_ArcGISsymbology.lyrx: Used to apply the suggested symbology in ArcGIS Pro.ClassificationProbability_yyyy-yyyy.tif: The probability values belonging to each class for every pixel. They were the intermediate products generated from the classification workflow and used to determine the final class with the highest probability and compute the forest index in this study. The sum of probabilities for all six classes is equal to 1. A scale factor of 10000 was applied to the GeoTiff files for storage convenience.HabitatMapHK_8class_yyyy-yyyy.tif: Based on the 6-class outputs, two more classes are added in this product, including wetland (pixel value 7) and plantation (pixel value 8), to serve as inputs for the habitat quality model.HabitatQualityHK_yyyy-yyyy.tif: Habitat quality maps produced in this study. The pixel value is a continuous variable ranging from 0 to 1, with 1 meaning the highest habitat quality.GIS supplementary data:All datasets were collected and compiled from January to June 2024 and represent the conditions at that time.Environmental Raster:DistanceFromCoast.tif: Geometric distance (m) from the coastline.Elevation.tif: Terrain height (m) from a LiDAR-based digital terrain model.Hillfire_10periods.tif: Hill fires occurred in each five-year period, based on burn-area products by Chan et al. (2023) and manual digitisation for early years.Insolation.tif: Annual amount of incoming solar radiation (kWh/m2) computed using SAGA GIS.Landslide_10periods.tif: Landslides occurred in each five-year period, based on the Enhanced Natural Terrain Landslide Inventory (Dias et al., 2009).Northness.tif: Terrain aspect from 1 (due north) to -1 (due south) computed from the DTM.Precipitation.tif: Annual precipitation (mm) (average between 1991-2020) from Hong Kong Observatory.Slope.tif: Steepness (°) of the ground surface computed from the DTM.SoilCEC.tif: Cation exchange capacity (CEC) (mmol/kg) of topsoil from Luo et al. (2007).SoilOrganicMatter.tif: Organic matter content (%) of topsoil from Luo et al. (2007).Temperature.tif: Annual mean temperature (°C) from Morgan and Guénard (2019).TopographicWetnessIndex.tif: Amount of water accumulation due to topographic effects computed using SAGA GIS.Typhoon_10periods.tif: Wind speed (km/h) estimated from WindNinja based on maximum hourly mean wind records associated with typhoon events in each five-year period.WindSpeed.tif: Mean wind speed (km/h) estimated from WindNinja based on monthly prevailing wind records.Human Activities:BuiltupAreas_10periods_shp.zip: Shapefile (polygons) of built-up areas, with attributes on the years of construction (estimated from topographic maps) and density (high and low). It was used as a threat factor in habitat quality mapping and variables in habitat changes.CountryParksProtectedAreas_shp.zip: Shapefile (polygons) of protected areas (Country Parks, Special Areas, etc.), with attributes on the years of designation and revision. It was used as a protection factor in habitat quality mapping and variables in habitat changes.PollutionSource_shp.zip: Shapefile (polygons) of pollution sources (landfills, power stations, and incineration plants), with attributes on the years of construction and closure. It was used as a threat factor in habitat quality mapping.Roads_10periods_shp.zip: Shapefile (polylines) of roads, with attributes on the years of construction (estimated from topographic maps) and type (main and secondary). It was used as a threat factor in habitat quality mapping.Mapping Reference:ForestIndex_FieldCollectedReferenceData.csv: Field survey records of habitat types which were used to evaluate the forest index variable in this study.HabitatMapHK_FieldCollectedReferenceData.csv: Field survey records of habitat types which were used to assess the habitat mapping results in this study.HabitatMapHK_OfficeInterpretedReferenceData.csv: Reference points where the habitat class in each period was determined through visual interpretation of the aerial photographs and other historical records. The points were used for both training and validation of the habitat maps in this study.HabitatQualityHK_FieldSurveyedEcologicalValue2008.csv: Field survey records of ecological values in 2008 which were used to evaluate the habitat quality maps in this study.LandsatHK_CrossSensorCalibrationPoints.csv: Selected points that were assumed to remain unchanged over time and used to cross-calibrate different Landsat sensors in this study.LandsatHK_ImageMetadata.csv: Metadata of the Landsat imagery (1,100 downloaded scenes and 607 valid scenes after pre-processing) acquired and processed in this study.Plantation_1975_1990_2008_2019.tif: Pixels that were identified as plantations on four existing maps in different years (1975, 1990, 2008, 2019), as represented by the four layers contained in this raster file respectively. These pixels were used to help extract plantation class on the habitat map (when producing habitat quality) and denote areas with plantation activities (when modelling habitat changes) in this study.SpeciesObsHK_SpeciesChecklist.csv: A species checklist of 7 taxa in Hong Kong (Plants, Butterflies, Birds, Reptiles, Dragonflies, Amphibians, Mammals) compiled from AFCD, Hong Kong Biodiversity Information Hub, and other secondary sources. Species of conservation concern are identified based on local assessments (Corlett et al., 2000; Fellowes et al., 2002), environmental protection laws, and national and global assessments. The checklist was used to match with the iNaturalist observation data to compute biodiversity metrics at grid levels and evaluate habitat quality maps in this study.SpeciesObsHK_SynonymList.csv: A list of species name synonyms for matching names used in iNaturalist and other secondary sources with the species checklist. It was used to pre-process the iNaturalist observation data and unify the species names from different records in this study.Analysis scripts:Part 1: Mapping Vegetation Habitats from a Satellite Image Time-SeriesP1_01_SearchAndDownloadFromGEE.ipynb: Query and download all available Landsat 1-9 imagery covering the study area using Google Earth Engine. Atmospheric correction is performed if necessary.P1_02_Preprocess_part1.py: Some basic pre-processing steps after downloading the images from cloud platform to local computer, such as mosaicking adjacent scenes and reprojecting to local coordinate system.P1_03_TopographicCorrection.R: SCS+C topographic correction based on terrain slope, aspect, sun azimuth and sun elevation angles.P1_04_CrossSensorCal.R: Cross-calibration of different Landsat sensors based on pseudo-invariant features, followed by computing variables for image classification.P1_05_ImageComposite.R: Create image composites (median and standard deviation statistics) by combining all imagery acquired in the same period.P1_06_ExtractPixelValue.R: Extract pixel values at the locations of reference points.P1_07_TrainingDataStat.R: Summarise the characteristics of pixel values (e.g., spectral reflectance) of each habitat class and Landsat sensor.P1_08_TrainRFModel.R: Train the Random Forest model, fuse probability outputs from each image, evaluate the model accuracies with cross-validation, and create the final model for classifying the entire dataset.P1_09_TestProcedures.R: Modify the classification procedures and re-run the Random Forest models to evaluate their impacts on the classification accuracies.P1_10_ApplyModel.R: Apply the Random Forest model and fusion steps to all images to create the habitat map for each period.P1_11_AreaCoverage.R: Obtain the area coverage of each class on the habitat map as well as the confidence interval of the area estimates.P1_12_CompareFieldData.R: Assess the accuracies of the habitat maps by overlaying with field-collected points and LiDAR height information at different times.P1_13_SurvivalAnalysis.R: Analyse the number of years required for transitioning between vegetation classes as well as the correlations between transition times and environmental variables.Part 2: Computing Habitat Quality Maps with Reference to

Wadi Hasa Sample Dataset — GRASS GIS Location
Version 1.0 (2025-09-19)

Overview
--------
This archive contains a complete GRASS GIS *Location* for the Wadi Hasa region (Jordan), including base data and exemplar analyses used in the Geomorphometry chapter. It is intended for teaching and reproducible research in archaeological GIS.

How to use
----------
1) Unzip the archive into your GRASSDATA directory (or a working folder) and add the Location to your GRASS session.
2) Start GRASS and open the included workspace (Workspace.gxw) or choose a Mapset to work in.
3) Set the computational region to the default extent/resolution for reproducibility:
g.region n=3444220 s=3405490 e=796210 w=733450 nsres=30 ewres=30 -p
4) Inspect layers as needed:
g.list type=rast,vector
r.info

Citation & License
------------------
Please cite this dataset as:

Isaac I. Ullah. 2025. *Wadi Hasa Sample Dataset (GRASS GIS Location)*. Zenodo. https://doi.org/10.5281/zenodo.17162040

All contents are released under the Creative Commons Attribution 4.0 International (CC BY 4.0) license. The original Wadi Hasa survey dataset is available at: https://figshare.com/articles/dataset/Wadi_Hasa_Ancient_Pastoralism_Project/1404216 The original Wadi Hasa survey dataset is available at: https://figshare.com/articles/dataset/Wadi_Hasa_Ancient_Pastoralism_Project/1404216

Coordinate Reference System
---------------------------
- Projection: UTM, Zone 36N
- Datum/Ellipsoid: WGS84
- Units: meter
- Coordinate system and units are defined in the GRASS Location (PROJ_INFO/UNITS).

Default Region (computational extent & resolution)
--------------------------------------------------
- North: 3444220
- South: 3405490
- East: 796210
- West: 733450
- Resolution: 30 (NS), 30 (EW)
- Rows x Cols: 1291 x 2092 (cells: 2700772)

Directory / Mapset Structure
----------------------------
This Location contains the following Mapsets (data subprojects), each with its own raster/vector layers and attribute tables (SQLite):
- Boolean_Predictive_Modeling: 8 raster(s), 4 vector(s)
- ISRIC_soilgrid: 31 raster(s), 0 vector(s)
- Landsat_Imagery: 3 raster(s), 0 vector(s)
- Landscape_Evolution_Modeling: 41 raster(s), 0 vector(s)
- Least_Cost_Analysis: 13 raster(s), 4 vector(s)
- Machine_Learning_Predictive_Modeling: 70 raster(s), 11 vector(s)
- PERMANENT: 4 raster(s), 2 vector(s)
- Sentinel2_Imagery: 4 raster(s), 0 vector(s)
- Site_Buffer_Analysis: 0 raster(s), 2 vector(s)
- Terrain_Analysis: 27 raster(s), 2 vector(s)
- Territory_Modeling: 14 raster(s), 2 vector(s)
- Trace21k_Paleoclimate_Downscale_Example: 4 raster(s), 2 vector(s)
- Visibility_Analysis: 11 raster(s), 5 vector(s)

Data Content (summary)
----------------------
- Total raster maps: 230
- Total vector maps: 34

Raster resolutions present:
- 10 m: 13 raster(s)
- 30 m: 183 raster(s)
- 208.01 m: 2 raster(s)
- 232.42 m: 30 raster(s)
- 1000 m: 2 raster(s)

Major content themes include:
- Base elevation surfaces and terrain derivatives (e.g., DEMs, slope, aspect, curvature, flow accumulation, prominence).
- Hydrology, watershed, and stream-related layers.
- Visibility analyses (viewsheds; cumulative viewshed analyses for Nabataean and Roman towers).
- Movement and cost-surface analyses (isotropic/anisotropic costs, least-cost paths, time-to-travel surfaces).
- Predictive modeling outputs (boolean/inductive/deductive; regression/classification surfaces; training/test rasters).
- Satellite imagery products (Landsat NIR/RED/NDVI; Sentinel‑2 bands and RGB composite).
- Soil and surficial properties (ISRIC SoilGrids 250 m products).
- Paleoclimate downscaling examples (CHELSA TraCE21k MAT/AP).

Vectors include:
- Archaeological point datasets (e.g., WHS_sites, WHNBS_sites, Nabatean_Towers, Roman_Towers).
- Derived training/testing samples and buffer polygons for modeling.
- Stream network and paths from least-cost analyses.

Important notes & caveats
-------------------------
- Mixed resolutions: Analyses span 10 m (e.g., Sentinel‑2 composites, some derived surfaces), 30 m (majority of terrain and modeling rasters), ~232 m (SoilGrids products), and 1 km (CHELSA paleoclimate). Set the computational region appropriately (g.region) before processing or visualization.
- NoData handling: The raw SRTM import (Hasa_30m_SRTM) reports extreme min/max values caused by nodata placeholders. Use the clipped/processed DEMs (e.g., Hasa_30m_clipped_wshed*) and/or set nodata with r.null as needed.
- Masks: MASK rasters are provided for analysis subdomains where relevant.
- Attribute tables: Vector attribute data are stored in per‑Mapset SQLite databases (sqlite/sqlite.db) and connected via layer=1.

Provenance (brief)
------------------
- Primary survey points and site datasets derive from the Wadi Hasa projects (see Figshare record above).
- Base elevation and terrain derivatives are built from SRTM and subsequently processed/clipped for the watershed.
- Soil variables originate from ISRIC SoilGrids (~250 m).
- Paleoclimate examples use CHELSA TraCE21k surfaces (1 km) that are interpolated to higher resolutions for demonstration.
- Satellite imagery layers are derived from Landsat and Sentinel‑2 scenes.

Reproducibility & quick commands
--------------------------------
- Restore default region: g.region n=3444220 s=3405490 e=796210 w=733450 nsres=30 ewres=30 -p
- Set region to a raster: g.region raster=

Change log
----------
- v1.0: Initial public release of the teaching Location on Zenodo (CC BY 4.0).

Contact
-------
For questions, corrections, or suggestions, please contact Isaac I. Ullah

This data set contains reduced-resolution QuickBird imagery and geospatial data for the entire Barrow QuickBird image area 156.15° W - 157.07° W, 71.15° N - 71.41° N) and the Barrow B4 Quadrangle (156.29° W - 156.89° W, 71.25° N - 71.40° N), for use in Geographic Information Systems (GIS) and remote sensing software. The original QuickBird data sets were acquired by DigitialGlobe from 1 to 2 August 2002, and consist of orthorectified satellite imagery. Federal Geographic Data Committee (FGDC)-compliant metadata for all value-added data sets are provided in text, HTML, and XML formats.

Accessory layers include: 1:250,000- and 1:63,360-scale USGS Digital Raster Graphic (DRG) mosaic images (GeoTIFF format); 1:250,000- and 1:63,360-scale USGS quadrangle index maps (ESRI Shapefile format); an index map for the 62 QuickBird tiles (ESRI Shapefile format); and a simple polygon layer of the extent of the Barrow QuickBird image area and the Barrow B4 quadrangle area (ESRI Shapefile format).

The baseline geospatial data support education, outreach, and multi-disciplinary research of environmental change in Barrow, which is an area of focused scientific interest.

Data are available either via FTP or on CD-ROM.

For a series of studies on the ecosystem service values of chaparral in Southern California, we developed a raster data layer providing an ecological unit classification of the Southern California landscape. This raster dataset is at a 30 meter pixel resolution and partitions the landscape into 37 different ecological unit types. This dataset was derived through a GIS-based cluster analysis of 10 different physiographic variables, namely soil suborder type, terrain geomorphon type, flow accumulation, slope, solar irradiation, annual precipitation, annual minimum temperature, actual evapotranspiration, and climatic water deficit. This partitioning was based on physiographic variables rather than vegetation types because of the wish to have the ecological units reflect biophysical characteristics rather than the historical land use patterns that may influence vegetation. The cluster analysis was performed across a set of 10,000 points randomly placed on a GIS layer stack for the 10 variables. These random points were grouped into 37 discrete clusters using an algorithm called partitioning around medoids. This assignment of points to clusters was then used to train a random forest classifier, which in turn was run across the GIS stack to produce the output raster layer.

This dataset is described in the following book chapter publication:

Underwood, Emma C., Allan D. Hollander, Patrick R. Huber, and Charlie Schrader-Patton. 2018. “Mapping the Value of National Forest Landscapes for Ecosystem Service Provision.” In Valuing Chaparral, 245–70. Springer Series on Environmental Management. Springer, Cham. https://doi.org/10.1007/978-3-319-68303-4_9.

Methods Summary of Methods for Developing Ecological Units in Southern California

Allan Hollander and Emma Underwood, University of California Davis.

1) Compiling GIS layers. These data were compiled from a variety of sources and resolutions (Table 1) for the southern California study area (see Methods_figure_1.png for the study area). The original resolution of these raster layers ran from 10 meters to 270 meters, and resampling was conducted so all analyses were performed at a 30 meter raster resolution. We decided not to include vegetation in the data stack as the aim was to capture biophysical characteristics and vegetation will reflect current landscape history and land use patterns (e.g. fire history, type conversion from shrubland, or agricultural use). Lakes and reservoirs were omitted from the subsequent analysis. Data compiled:

a) Soil suborders. This was a discretely-classified raster layer with 22 soil suborder classes included in the southern California region. This was derived from the gridded Soil Survey Geographic Database (gSSURGO, available at http://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/geo/?cid=nrcs142p2_053628). This product is a rasterization at a 10-meter resolution of the county-scale SSURGO data published by the USDA Natural Resources Conservation Service.

b) Terrain geomorphons. This raster layer derives from a DEM surface and classifies the landscape into 10 discrete landform types, examples being ridges, slopes, hollows, and valleys. The algorithm for geomorphon classification uses a pattern recognition approach based on line of sight analysis (Jasiewisc and Stepinski 2013). This layer was created from a 30 meter DEM in GRASS 7.0.0, using the extension r.geomorphon (https://grass.osgeo.org/grass70/manuals/addons/r.geomorphon.html).

c) Annualized solar irradiation. This layer uses the r.sun model available in GRASS 7.0.0 (https://grass.osgeo.org/grass70/manuals/r.sun.html) which calculates direct, diffuse, and reflected solar irradiation for a given day, location, topography, and atmospheric conditions. This layer was created from a 30 meter DEM and assumes clear-sky conditions. To estimate the total annual irradiation, the model was run for every 15th day and these values were integrated over the year.

d) Flow accumulation. This layer is another product of 30 meter DEM data and measures the upslope area in pixel count that conceivably drains into a given pixel. This was calculated using the accumulation option in the GRASS 7.0.0 command r.watershed (https://grass.osgeo.org/grass70/manuals/r.watershed.html)

e) Slope. This was derived from 30 meter DEM data using the GRASS 7.0.0 command r.slope.aspect, and is measured in degrees.

f) Annual precipitation. This layer came from the 2014 Basin Characterization Model (BCM) for California (Flint et al. 2013) and gives the average annual precipitation between 1981 and 2010 at a 270-meter resolution.

g) Annual minimum temperature. This layer also came from BCM (Flint et al. 2013) and gives the average annual minimum temperature between 1981 and 2010 at a 270-meter resolution. Minimum temperature was included in the set of climate variables to represent montane winter conditions.

h) Climatic water deficit. This layer also came from the BCM (Flint et al. 2013) and gives the average climatic water deficit between 1981 and 2010 at a 270-meter resolution. The two evapotranspiration variables (climatic water deficit and actual evapotranspiration) are included in this set because they are strong drivers of vegetation distribution (Stephenson 1998).

i) Actual evapotranspiration. This layer also came from the BCM (Flint et al. 2013) and gives the average actual evapotranspiration between 1981 and 2010 at a 270-meter resolution.

Table 1. Summary of GIS data stack

    LAYER


    ORIGINAL SOURCE


    ORIGINAL RESOLUTION


    THEME








    Soil suborders


    gSSURGO


    10 meters


    Soil type




    Terrain geomorphons


    Digital elevation model


    30 meters


    Geomorphometry




    Solar irradiation


    Digital elevation model


    30 meters


    Energy balance




    Flow accumulation


    Digital elevation model


    30 meters


    Geomorphometry




    Slope


    Digital elevation model


    30 meters


    Geomorphometry




    Annual precipitation


    Basin Characterization Model


    270 meters


    Climate




    Annual min temperature


    Basin Characterization Model


    270 meters


    Climate




    Climatic water deficit


    Basin Characterization Model


    270 meters


    Climate




    Actual evapotranspiration


    Basin Characterization Model


    270 meters


    Climate

2) Generating 10,000 random points. A mask was imposed to limit analyses to the 35,158 square study area and 10,000 random points were generated to create a data table of the values of each GIS layer at each of the random points. This data table was the basis for sorting the random points into a limited number of clustered types. The first step in doing this is calculating in multivariate space the distance with respect to these environmental variables each random point is from every other point, in other words creating a dissimilarity matrix.

3) Assigning weights to variables. Because the 9 environmental variables use completely different metrics and are a combination of numerical and categorical types, calculating an environmental distance between any two of these random points requires some weighting to be assigned to each of the environmental variables to sum up their relative distances. A subanalysis to determine these weightings used a subset of the study area, the Santa Clara River watershed. Since these ecological units are intended to summarize a diverse set of ecological services, we chose three different proxy variables from the GIS data available for this area to represent biomass, hydrological response, and biodiversity. These proxies included mean annual MODIS Enhanced Vegetation Index (EVI) value for biomass, recharge for hydrological response, and habitat type in the California Wildlife Habitat Relations (CWHR) classification for biodiversity.

The MODIS EVI data was derived by averaging over the 2000-2014 period the maximum EVI value in a single year. The MODIS index used was MOD13Q1 (https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod13q1) at a 250 meter resolution, available at 16-day intervals.


The hydrological recharge data were extracted from the 2014 Basin Characterization Model (Flint et al. 2013) at 270 meter resolution.


The CWHR habitat type came from the 2015 FRAP vegetation layer (FVEG15_1, from http://frap.fire.ca.gov/data/frapgisdata-sw-fveg_download), available at a 30 meter resolution.

a) We used random forest regression and classification (Hastie et al. 2009) to determine a ranking of importance values of these predictor variables using random forest regression for EVI and recharge and random forest classification for the habitat type. These were calculated using the randomForest package in R (Liaw and Wiener 2002).

b) We then averaged these three sets of importance values to create an overall set of weightings to enter into the dissimilarity matrix (Table 2).

Table 2. Weightings for each variable to reflect their relative importance to the ecological units

    VARIABLE NAME


    WEIGHT




    Precipitation


    1.00




    Annual minimum temperature


    0.600




    Slope


    0.507




    Climatic water deficit


    0.413




    Annualized solar radiation


    0.404




    Soil suborder