Aim

To understand the representativeness and accuracy of expert range maps, and explore alternate methods for accurately mapping species distributions.

Location

Global

Time period

Contemporary

Major taxa studied

Terrestrial vertebrates, and Odonata

Methods

We analyzed the biases in 50,768 animal IUCN, GARD and BirdLife species maps, assessed the links between these maps and existing political and various non-ecological boundaries to assess their accuracy for certain types of analysis. We cross-referenced each species map with data from GBIF to assess if maps captured the whole range of a species, and what percentage of occurrence points fall within the species’ assessed ranges. In addition, we use a number of alternate methods to map diversity patterns and compare these to high resolution models of distribution patterns.

Results

On average 20-30% of species’ non-coastal range boundaries overlapped with administrative national boundaries. In total, 60% of areas with the highest spatial turnover in species (high densities of species range boundaries marking high levels of shift in the community of species present) occurred at political boundaries, especially commonly in Southeast Asia. Different biases existed for different taxa, with gridded analysis in reptiles, river-basins in Odonata (except the Americas) and county-boundaries for Amphibians in the US. On average, up to half (25-46%) species recorded range points fall outside their mapped distributions. Filtered Minimum-convex polygons performed better than expert range maps in reproducing modeled diversity patterns.

Main conclusions

Expert range maps showed high bias at administrative borders in all taxa, but this was highest at the transition from tropical to subtropical regions. Methods used were inconsistent across space, time and taxa, and ranges mapped did not match species distribution data. Alternate approaches can better reconstruct patterns of distribution than expert maps, and data driven approaches are needed to provide reliable alternatives to better understand species distributions.

Methods Materials and methods

We use a combination of approaches to explore the relationship between species range maps and geopolitical boundaries and a subset of geographic features. In some cases we used the density of species range boundaries to explore the relationship between these and various features (i.e. administrative boundaries, river basin boundaries etc.). Additionally, species richness and spatial turnover are used to explore changes in richness over short geographic distances. Analyses were conducted in R statistical software unless noted otherwise. All code scripts are available at https://github.com/qiaohj/iucn_fix. Workflows are shown in Figure S1a-c with associated scripts listed.

Species ranges and boundary density maps

ERMs (Expert range maps) were downloaded from the IUCN RedList website for mammals (5,709 species), odonates (2,239 species) and amphibians (6,684 species; https://www.iucnredlist.org/resources/grid/spatial-data). Shapefile maps for birds were downloaded from BirdLife (10,423 species, http://datazone.birdlife.org/species/requestdis), and for reptiles from the Global Assessment of Reptile Distributions (GARD) (10,064 species; Roll et al., 2017). Each species’ polygon boundaries were converted to a polylines to show the boundary of each species range (Figure S1a-II; codes are lines 7 – 18 in line2raster_xxxx.r ; xxxx varies based on the taxa). The associated shapefile was then split to produce independent polyline files for each species within each taxon (see Figure S1a-I, codes are lines 29 to 83 in the same file above.).

To generate species boundary density maps, species range boundaries were rasterized at 1km spatial resolution with an equal area projection (Eckert-IV), and stacked to form a single raster for each taxon (at the level of amphibians, odonates, etc.). This represented the number of species in each group and their overlapping range boundaries (Figure S1b-II, codes are in line2raster_all.r). Each cell value indicated the number of species whose distribution boundaries overlapped with each cell, enabling us to overlay this rasterized information with other features (i.e. administrative boundaries) so that the overlaps between them can be calculated in R. These species boundary density maps underlie most subsequent analyses. R code and caveats are given in the supplements, links are provided in text and Figure S1.

Geographic boundaries

Spatial exploration of species range boundaries in ArcGIS suggested that numerous geographic datasets (i.e. political and in few cases geographic features such as river basins) were used to delineate the species ranges for different regions and taxa (this is sometimes part of the methodology in developing ERMs as detailed by Ficetola et al., 2014). Thus in addition to analyzing the administrative bias and the percentage of occurrence records within each species’ ERM for all taxa, additional analyses were conducted when other biases were evident in any given taxa or region (detailed later in methods on a case-by-case basis).

For all taxa, we assessed the percentage of overlap between species range boundaries and national and provincial boundaries by digitizing each to 1km (equivalent to buffering thie polyline by 500m), both with and without coastal boundaries. An international map was used because international (Western) assessors use them, and does not necessarily denote agreed country boundaries (https://gadm.org/). The different buffers (500m, 1000m, 2500m, 5000m) were added to these administrative boundaries in ArcMap to account for potential, insignificant deviations from political boundaries (Figure S1b). An R script for the same function is provided in “country_line_buffer.r”.

To establish where multiple species shared range boundaries we reclassified the species range boundary density rasters for each taxa into richness classes using the ArcMap quartile function (Figure S1). From these ten classes the percentage of the top-two, and top-three quartiles of range densities within different buffers (500m, 1000m, 2500m, 5000m) was calculated per country to determine what percentage of highest range boundary density approximately followed administrative borders. This was done because people drawing ERMs may use detailed administrative maps or generalize near political borders, or may use political shapefiles that deviate slightly. It is consequently useful to include varying distances from administrative features to assess how range boundary densities vary in relation to administrative boundaries. Analyses of relationships between individual species range boundaries and administrative boundaries (coastal, non-coastal) were made in R and scripts provided (quantile_country_buffer_overlap.r).

Spatial turnover and administrative boundaries

Heatmaps of species richness were generated by summing entire sets of compiled species ranges for each taxon in polygonal form (Figure 1; Figure S1b-I). To assess abrupt diversity changes, standard deviations for 10km blocks were calculated using the block statistics function in ArcMap. Abrupt changes in diversity were signified by high standard deviations based on the cell statistics function in ArcGIS, which represented rapid changes in the number of species present. Maps were then classified into ten categories using the quartile function. Given the high variation in maximum diversity and taxonomic representation, only the top two –three richness categories were retained per taxon. This was then extracted using 1km buffers of national administrative boundaries to assess percentages of administrative boundaries overlapping turnover hotspots by assessing what proportion of political boundaries were covered by these turnover hotspots.

Taxon-specific analyses

Data exploration and mapping exposed taxon and regional-specific biases requiring additional analysis. Where other biases and irregularities were clear from visual inspection of the range boundary density maps for each taxa, the possible causes of biases were assessed by comparing range boundary density maps to high-resolution imagery and administrative maps via the ArcGIS server (AGOL). Standardized overlay of the taxon boundary sets with administrative or geophysical features from the image-server revealed three types of bias which were either spatially or taxonomically limited between: 1) amphibians with county borders in the United States, 2) dragonflies and river basins globally and 3) gridding of distributions of reptiles. In these cases, species boundary density maps were used as a basis to identify potential biases which were then explored empirically using appropriate methods.

For amphibians, counties in the United States (US) were digitized using a county map from the US (https://gadm.org/), then buffered by with 2.5km either side. Amphibian species range boundary density maps were reclassified showing where species range boundaries existed (with other non-range boundary areas reclassified as “no data,”) and all species boundaries numerically indicated (i.e. values of 1 indicates one species range boundary, values of 10 indicates ten species range boundaries). Percentages of species boundary areas falling on county and in the buffers, in addition to species range boundaries which did not overlap with county boundaries were calculated to give measures of what percentage of the species boundaries fell within 2.5km of county boundaries.

For Odonata, many species were mapped to river basin borders. We used river basins of levels 6-8 (sub-basin to basin) in the river hierarchy (https://hydrosheds.org) to assess the relationship between Odonata boundaries and river boundaries. Two IUCN datasets exist for Odonata; the IUCN Odonata specialist group spatial dataset

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map California’s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of “lands” from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of California’s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bay’s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Service–Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Offshore of Coal Oil Point map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these “ground-truth” surveying data are available from the CSMP Video and Photograph Portal at https://doi.org/10.5066/F7J1015K. The “seafloor character” data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The “potential habitats” polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Offshore Coal Oil Point map area data layers. Data layers are symbolized as shown on the associated map sheets.

The LANDFIRE vegetation layers describe the following elements of existing and potential vegetation for each LANDFIRE mapping zone: environmental site potentials, biophysical settings, existing vegetation types, canopy cover, and vegetation height. Vegetation is mapped using predictive landscape models based on extensive field reference data, satellite imagery, biophysical gradient layers, and classification and regression trees. The environmental site potential (ESP) data layer represents the vegetation that could be supported at a given site based on the biophysical environment. Map units are named according to NatureServe's Ecological Systems classification, which is a nationally consistent set of mid-scale ecological units (Comer and others 2003). Usage of these classification units to describe environmental site potential, however, differs from the original intent of Ecological Systems as units of existing vegetation. As used in LANDFIRE, map unit names represent the natural plant communities that would become established at late or climax stages of successional development in the absence of disturbance. They reflect the current climate and physical environment, as well as the competitive potential of native plant species. The ESP layer is similar in concept to other approaches to classifying potential vegetation in the western United States, including habitat types (for example, Daubenmire 1968 and Pfister and others 1977) and plant associations (for example, Henderson and others 1989). It is important to note that ESP is an abstract concept and represents neither current nor historical vegetation. To create the ESP data layer, we first assign field plots to one of the ESP map unit classes. Go to http://www.landfire.gov/participate_acknowledgements.php for more information regarding contributors of field plot data. Assignments are based on presence and abundance of indicator plant species recorded on the plots and on the ecological amplitude and competitive potential of these species. We then intersect plot locations with a series of 30-meter spatially explicit gradient layers. Most of the gradient layers used in the predictive modeling of ESP are derived using the WX-BGC simulation model (Keane and Holsinger, in preparation; Keane and others 2002). WX-BGC simulations are based largely on spatially extrapolated weather data from DAYMET (Thornton and others 1997; Thornton and Running 1999; http://www.daymet.org/) and on soils data in STATSGO (NRCS 1994). Additional indirect gradient layers, such as elevation, slope, and indices of topographic position, are also used. We use data from plot locations to develop predictive classification tree models, using See5 data mining software (Quinlan 1993; Rulequest Research 1997), for each LANDFIRE map zone. These decision trees are applied spatially to predict the ESP for every pixel across the landscape. Finally, ESP pixel values are, in some cases, modified based on a comparison with the LANDFIRE existing vegetation type (EVT) layer created with the use of 30-meter Landsat ETM satellite imagery. We make such modifications only in non-vegetated areas (such as water, rock, snow, or ice) and where information in the EVT layer clearly enables a better depiction of the environmental site potential concept. Although the ESP data layer is intended to represent current site potential, the actual time period for this data set is variable. The weather data used in DAYMET were compiled from 1980 to 1997. Refer to spatial metadata for date ranges of field plot data and satellite imagery for each LANDFIRE map zone. A number of changes were implemented for the LF2010 ESP product that worked with this original data. LF2010 updates to mapping EVT map units for Barren, Snow-Ice, and Water were translated to the LF2010 ESP product so those map units will coincide with the EVT. Subsequent to that, each ESP map unit was stratified spatially two different ways. First, each ESP map unit was stratified by LANDFIRE map zone. Second, each ESP map unit was stratified by an ESP life form classification layer that incorporated NLCD 2001 data, LF2001 EVC data, a Vegetation Change Tracker (VCT) dataset (Huang, 2010), and the National Wetlands Inventory (NWI) data. Each layer was leveraged against each other to determine areas of stable Sparse, Upland Herb, Upland Shrub, Upland Woodland, Upland Forest, Wetland Shrub-herb, Wetland Forest, Wetland Shrub, and Wetland Herb. Areas mapped as agriculture, urban, barren, snow-ice, and water were described as “Undetermined”.

Geospatial Data Pack for Visualization 🗺️

Learn Geographic Mapping with Altair, Vega-Lite and Vega using Curated Datasets

Complete geographic and geophysical data collection for mapping and visualization. This consolidation includes 18 complementary datasets used by 31+ Vega, Vega-Lite, and Altair examples 📊. Perfect for learning geographic visualization techniques including projections, choropleths, point maps, vector fields, and interactive displays.

Source data lives on GitHub and can also be accessed via CDN. The vega-datasets project serves as a common repository for example datasets used across these visualization libraries and related projects.

Why Use This Dataset? 🤔

• Comprehensive Geospatial Types: Explore a variety of core geospatial data models:
  • Vector Data: Includes points (like airports.csv), lines (like londonTubeLines.json), and polygons (like us-10m.json).
  • Raster-like Data: Work with gridded datasets (like windvectors.csv, annual-precip.json).
• Diverse Formats: Gain experience with standard and efficient geospatial formats like GeoJSON (see Table 1, 2, 4), compressed TopoJSON (see Table 1), and plain CSV/TSV (see Table 2, 3, 4) for point data and attribute tables ready for joining.
• Multi-Scale Coverage: Practice visualization across different geographic scales, from global and national (Table 1, 4) down to the city level (Table 1).
• Rich Thematic Mapping: Includes multiple datasets (Table 3) specifically designed for joining attributes to geographic boundaries (like states or counties from Table 1) to create insightful choropleth maps.
• Ready-to-Use & Example-Driven: Cleaned datasets tightly integrated with 31+ official examples (see Appendix) from Altair, Vega-Lite, and Vega, allowing you to immediately practice techniques like projections, point maps, network maps, and interactive displays.
• Python Friendly: Works seamlessly with essential Python libraries like Altair (which can directly read TopoJSON/GeoJSON), Pandas, and GeoPandas, fitting perfectly into the Kaggle notebook environment.

Table of Contents

• Why Use This Dataset?
• Dataset Inventory
• Usage Guide
• Appendix: Geospatial Examples by Technique
• Licenses
• Sources
• Learning Resources

Dataset Inventory 🗂️

This pack includes 18 datasets covering base maps, reference points, statistical data for choropleths, and geophysical data.

1. BASE MAP BOUNDARIES (Topological Data)

2. GEOGRAPHIC REFERENCE POINTS (Point Data) 📍

A global data set of soil types is available at 0.5-degree latitude by 0.5-degree longitude resolution. There are 106 soil units, based on Zobler?s (1986) assessment of the FAO/UNESCO Soil Map of the World. This data set is a conversion of the Zobler 1-degree resolution version to a 0.5-degree resolution. The resolution of the data set was not actually increased. Rather, the 1-degree squares were divided into four 0.5-degree squares with the necessary adjustment of continental boundaries and islands. The computer code used to convert the original 1-degree data to 0.5-degree is provided as a companion file. A JPG image of the data is provided in this document. The Zobler data (1-degree resolution) as distributed by Webb et al. (1993) [http://www.ngdc.noaa.gov/seg/eco/cdroms/gedii_a/datasets/a12/wr.htm#top] contains two columns, one column for continent and one column for soil type. The Soil Map of the World consists of 9 maps that represent parts of the world. The texture data that Webb et al.(1993) provided allowed for the fact that a soil type in one part of the world may have different properties than the same soil in a different part of the world. This continent-specific information is retained in this 0.5-degree resolution data set, as well as the soil type information which is the second column. A code was written (one2half.c) to take the file CONTIZOB.LER distributed by Webb et al. (1993) [http://www.ngdc.noaa.gov/seg/eco/cdroms/gedii_a/datasets/a12/wr.htm#top] and simply divide the 1-degree cells into quarters. This code also reads in a land/water file (land.wave) that specifies the cells that are land at 0.5 degrees. The code checks for consistency between the newly quartered map and the land/water map to which the quartered map is to be registered. If there is a discrepancy between the two, an attempt was made to make the two consistent using the following logic. If the cell is supposed to be water, it is forced to be water. If it is supposed to be land but was resolved to water at 1 degree, the code looks at the surrounding 8 cells and picks the most frequent soil type and assigns it to the cell. If there are no surrounding land cells then it is kept as water in the hopes that on the next pass one or more of the surrounding cells might be converted from water to a soil type. The whole map is iterated 5 times. The remaining cells that should be land but couldn't be determined from surrounding cells (mostly islands that are resolved at 0.5 degree but not at 1 degree) are printed out with coordinate information. A temporary map is output with -9 indicating where data is required. This is repeated for the continent code in CONTIZOB.LER as well. A separate map of the temporary continent codes is produced with -9 indicating required data. A nearly identical code (one2half.c) does the same for the continent codes. The printout allows one to consult the printed versions of the soil map and look up the soil type with the largest coverage in the 0.5-degree cell. The program manfix.c then will go through the temporary map and prompt for input to correct both the soil codes and the continent codes for the map. This can be done manually or by preparing a file of changes (new_fix.dat) and redirecting stdin. A new complete version of the map is outputted. This is in the form of the original CONTIZOB.LER file (contizob.half) but four times larger. Original documentation and computer codes prepared by Post et al. (1996) are provided as companion files with this data set. Image of 106 global soil types available at 0.5-degree by 0.5-degree resolution. Additional documentation from Zobler?s assessment of FAO soil units is available from the NASA Center for Scientific Information.

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map California’s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of “lands” from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of California’s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bay’s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Service–Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Offshore of Tomales Point map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these “ground-truth” surveying data are available from the CSMP Video and Photograph Portal at https://doi.org/10.5066/F7J1015K. The “seafloor character” data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The “potential habitats” polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Offshore of Tomales Point map area data layers. Data layers are symbolized as shown on the associated map sheets.

Connecticut Quaternary Geology Geologic Basin Divides includes the drainage basins divides appearing on Sheet 1 of the The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005) The Connecticut Quaternary Geology digital spatial data combines the information portrayed on the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others 2005) with the information portrayed on its sister map, the Surficial Materials Map of Connecticut (Stone and others, 1992). When used together, these maps provide a three dimensional context for understanding and predicting the internal composition, resource potential and hydrologic character of Connecticut's glacial and postglacial deposits. Both were compiled at 1:24,000 scale, and published at 1:125,000 scale. The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005) portrays the glacial and postglacial deposits of Connecticut (including Long Island Sound) with an emphasis on where and how they were emplaced. Glacial Ice-Laid Deposits (thin till, thick till, and deposits of individual end moraines), Early Postglacial Deposits (Late Wisconsinan to Early Holocene stream terrace and inland dune deposits) and Holocene Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. This mapping is based on the concept of systematic northward retreat of the Late Wisconsinan glacier. Meltwater deposits are divided into six depositional system categories (Deposits of Major Ice-Dammed Lakes, Deposits of Major Sediment-Dammed Lakes, Deposits of Related Series of Ice-Dammed Ponds, Deposits of Related Series of Sediment-Dammed Ponds, Deposits of Proximal Meltwater Streams, and Deposits of Distal Meltwater Streams) based on the arrangement and character of the groupings of sedimentary facies (morphosequences). The Surficial Materials Map of Connecticut (Stone and others, 1992) portrays the glacial and postglacial deposits of Connecticut in terms of their aerial extent and subsurface textural relationships. Glacial Ice-Laid Deposits (thin till, thick till, end moraine deposits) and Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. The meltwater deposits are further characterized using four texturally-based map units (g = gravel, sg = sand and gravel, s = sand, and f = fines). In many places a single map unit (e.g. sand) is sufficient to describe the entire meltwater section. Where more complex stratigraphic relationships exist, "stacked" map units are used to characterize the subsurface (e.g. sg/s/f - sand and gravel overlying sand overlying fines). Where postglacial deposits overlie meltwater deposits, this relationship is also described (e.g. alluvium overlying sand). Map unit definitions (Surficial Materials Polygon Code definitions, found in the metadata) provide a short description of the inferred depositional environment for each of the glacial meltwater map units. The geologic contacts between till and meltwater deposits coincide on both the Quaternary and Surficial Materials maps, as do the boundaries of polygons that define areas of thick till, alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill. Within the meltwater deposits, a Quaternary map unit (deposit) may contain several Surficial Materials textural units (akin to facies within a delta, for example). Combining the textural and vertical stacking information from the Surficial Materials map with the orderly portrayal of morphosequence relationships, up and down valley, that can be gleaned from the Quaternary map provides a three dimensional predictive context for relating the geologic setting of Connecticut's glacial meltwater deposits to their behavior as aquifers and/or transmitters of contaminants. Since this data layer is a polygon and line feature representation of the two maps combined, each map unit's depiction and description could provide information as to its aerial extent, subsurface textural characteristics, depositional and paleogeographic settings, and facies composition in a morphosequence context. Therefore, a typical meltwater polygon would have a combination of Quaternary (e.g. Deposit of Major Sediment-Dammed Lake; Glacial Lake Middletown Cromwell Deltaic Deposit) and Surficial Materials (e.g. sand and gravel overlying sand overlying fine) map attributes. Additional polygon features are incorporated to define surface water areas for streams, lakes, ponds, bays, and estuaries greater than 5 acres in size. Line features describe the type of boundary between individual geologic or textural units such as a geologic contact line between two different geologic units or a linear shoreline feature between a textural unit and an adjacent waterbody. The data have been updated to reflect minor changes in map unit name (QUPOLY_COD) for consistency with the 2005 publication of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin. Previously distributed versions of CTQSGEOM were consistent with the 1998 Open-file Report for the same map. It is important to note that this data layer represents only the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005). The off-shore geologic units are organized in separate data layers (LISQMOR, LISQFAN, LISQLAKE, LISQCHAN, LISQMARD) which can be used in conjunction with this data layer. These Long Island Sound layers have been mapped at 1:80,000 scale using seismic reflection data. The CTQSGEOM data layer should be used as the geologic base for Connecticut Quaternary Geology / Surficial Materials Features (CTQSFEAT) data layer which represents features such as eskers, meltwater channels, spillways, and locations of radio-carbon dated samples.

Connecticut Quaternary Geology Geologic Cross Section includes the cross sections appearing on Sheets 1 and 2 of the The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005). The Connecticut Quaternary Geology digital spatial data combines the information portrayed on the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others 2005) with the information portrayed on its sister map, the Surficial Materials Map of Connecticut (Stone and others, 1992). When used together, these maps provide a three dimensional context for understanding and predicting the internal composition, resource potential and hydrologic character of Connecticut's glacial and postglacial deposits. Both were compiled at 1:24,000 scale, and published at 1:125,000 scale. The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005) portrays the glacial and postglacial deposits of Connecticut (including Long Island Sound) with an emphasis on where and how they were emplaced. Glacial Ice-Laid Deposits (thin till, thick till, and deposits of individual end moraines), Early Postglacial Deposits (Late Wisconsinan to Early Holocene stream terrace and inland dune deposits) and Holocene Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. This mapping is based on the concept of systematic northward retreat of the Late Wisconsinan glacier. Meltwater deposits are divided into six depositional system categories (Deposits of Major Ice-Dammed Lakes, Deposits of Major Sediment-Dammed Lakes, Deposits of Related Series of Ice-Dammed Ponds, Deposits of Related Series of Sediment-Dammed Ponds, Deposits of Proximal Meltwater Streams, and Deposits of Distal Meltwater Streams) based on the arrangement and character of the groupings of sedimentary facies (morphosequences). The Surficial Materials Map of Connecticut (Stone and others, 1992) portrays the glacial and postglacial deposits of Connecticut in terms of their aerial extent and subsurface textural relationships. Glacial Ice-Laid Deposits (thin till, thick till, end moraine deposits) and Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. The meltwater deposits are further characterized using four texturally-based map units (g = gravel, sg = sand and gravel, s = sand, and f = fines). In many places a single map unit (e.g. sand) is sufficient to describe the entire meltwater section. Where more complex stratigraphic relationships exist, "stacked" map units are used to characterize the subsurface (e.g. sg/s/f - sand and gravel overlying sand overlying fines). Where postglacial deposits overlie meltwater deposits, this relationship is also described (e.g. alluvium overlying sand). Map unit definitions (Surficial Materials Polygon Code definitions, found in the metadata) provide a short description of the inferred depositional environment for each of the glacial meltwater map units. The geologic contacts between till and meltwater deposits coincide on both the Quaternary and Surficial Materials maps, as do the boundaries of polygons that define areas of thick till, alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill. Within the meltwater deposits, a Quaternary map unit (deposit) may contain several Surficial Materials textural units (akin to facies within a delta, for example). Combining the textural and vertical stacking information from the Surficial Materials map with the orderly portrayal of morphosequence relationships, up and down valley, that can be gleaned from the Quaternary map provides a three dimensional predictive context for relating the geologic setting of Connecticut's glacial meltwater deposits to

Connecticut Quaternary Geology Long Island Submerged Marine Fluvial-Estuarine, Channel-Fill Deposits identifies early postglacial, channel-fill deposits submerged in Long Island Sound and Fishers Island Sound. This information appears on Sheet 1 of the The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005). The Connecticut Quaternary Geology digital spatial data combines the information portrayed on the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others 2005) with the information portrayed on its sister map, the Surficial Materials Map of Connecticut (Stone and others, 1992). When used together, these maps provide a three dimensional context for understanding and predicting the internal composition, resource potential and hydrologic character of Connecticut's glacial and postglacial deposits. Both were compiled at 1:24,000 scale, and published at 1:125,000 scale. The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005) portrays the glacial and postglacial deposits of Connecticut (including Long Island Sound) with an emphasis on where and how they were emplaced. Glacial Ice-Laid Deposits (thin till, thick till, and deposits of individual end moraines), Early Postglacial Deposits (Late Wisconsinan to Early Holocene stream terrace and inland dune deposits) and Holocene Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. This mapping is based on the concept of systematic northward retreat of the Late Wisconsinan glacier. Meltwater deposits are divided into six depositional system categories (Deposits of Major Ice-Dammed Lakes, Deposits of Major Sediment-Dammed Lakes, Deposits of Related Series of Ice-Dammed Ponds, Deposits of Related Series of Sediment-Dammed Ponds, Deposits of Proximal Meltwater Streams, and Deposits of Distal Meltwater Streams) based on the arrangement and character of the groupings of sedimentary facies (morphosequences). The Surficial Materials Map of Connecticut (Stone and others, 1992) portrays the glacial and postglacial deposits of Connecticut in terms of their aerial extent and subsurface textural relationships. Glacial Ice-Laid Deposits (thin till, thick till, end moraine deposits) and Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. The meltwater deposits are further characterized using four texturally-based map units (g = gravel, sg = sand and gravel, s = sand, and f = fines). In many places a single map unit (e.g. sand) is sufficient to describe the entire meltwater section. Where more complex stratigraphic relationships exist, "stacked" map units are used to characterize the subsurface (e.g. sg/s/f - sand and gravel overlying sand overlying fines). Where postglacial deposits overlie meltwater deposits, this relationship is also described (e.g. alluvium overlying sand). Map unit definitions (Surficial Materials Polygon Code definitions, found in the metadata) provide a short description of the inferred depositional environment for each of the glacial meltwater map units. The geologic contacts between till and meltwater deposits coincide on both the Quaternary and Surficial Materials maps, as do the boundaries of polygons that define areas of thick till, alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill. Within the meltwater deposits, a Quaternary map unit (deposit) may contain several Surficial Materials textural units (akin to facies within a delta, for example). Combining the textural and vertical stacking information from the Surficial Materials map with the orderly portrayal of morphosequence relationships, up and down valley, that can be gleaned from the Quaternary map provides a three dimensional predictive context for relating the geologic setting of Connecticut's glacial meltwater deposits to their behavior as aquifers and/or transmitters of contaminants. Since this data layer is a polygon and line feature representation of the two maps combined, each map unit's depiction and description could provide information as to its aerial extent, subsurface textural characteristics, depositional and paleogeographic settings, and facies composition in a morphosequence context. Therefore, a typical meltwater polygon would have a combination of Quaternary (e.g. Deposit of Major Sediment-Dammed Lake; Glacial Lake Middletown Cromwell Deltaic Deposit) and Surficial Materials (e.g. sand and gravel overlying sand overlying fine) map attributes. Additional polygon features are incorporated to define surface water areas for streams, lakes, ponds, bays, and estuaries greater than 5 acres in size. Line features describe the type of boundary between individual geologic or textural units such as a geologic contact line between two different geologic units or a linear shoreline feature between a textural unit and an adjacent waterbody. The data have been updated to reflect minor changes in map unit name (QUPOLY_COD) for consistency with the 2005 publication of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin. Previously distributed versions of CTQSGEOM were consistent with the 1998 Open-file Report for the same map. It is important to note that this data layer represents only the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005). The off-shore geologic units are organized in separate data layers (LISQMOR, LISQFAN, LISQLAKE, LISQCHAN, LISQMARD) which can be used in conjunction with this data layer. These Long Island Sound layers have been mapped at 1:80,000 scale using seismic reflection data. The CTQSGEOM data layer should be used as the geologic base for Connecticut Quaternary Geology / Surficial Materials Features (CTQSFEAT) data layer which represents features such as eskers, meltwater channels, spillways, and locations of radio-carbon dated samples.

Connecticut Quaternary Geology Long Island Submerged Marine Deltaic Deposits identifies early postglacial, marine deltaic deposits submerged in Long Island Sound. This information appears on Sheet 1 of the The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005). The Connecticut Quaternary Geology digital spatial data combines the information portrayed on the on-land portion of the Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others 2005) with the information portrayed on its sister map, the Surficial Materials Map of Connecticut (Stone and others, 1992). When used together, these maps provide a three dimensional context for understanding and predicting the internal composition, resource potential and hydrologic character of Connecticut's glacial and postglacial deposits. Both were compiled at 1:24,000 scale, and published at 1:125,000 scale. The Quaternary Geologic Map of Connecticut and Long Island Sound Basin (Stone and others, 2005) portrays the glacial and postglacial deposits of Connecticut (including Long Island Sound) with an emphasis on where and how they were emplaced. Glacial Ice-Laid Deposits (thin till, thick till, and deposits of individual end moraines), Early Postglacial Deposits (Late Wisconsinan to Early Holocene stream terrace and inland dune deposits) and Holocene Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. This mapping is based on the concept of systematic northward retreat of the Late Wisconsinan glacier. Meltwater deposits are divided into six depositional system categories (Deposits of Major Ice-Dammed Lakes, Deposits of Major Sediment-Dammed Lakes, Deposits of Related Series of Ice-Dammed Ponds, Deposits of Related Series of Sediment-Dammed Ponds, Deposits of Proximal Meltwater Streams, and Deposits of Distal Meltwater Streams) based on the arrangement and character of the groupings of sedimentary facies (morphosequences). The Surficial Materials Map of Connecticut (Stone and others, 1992) portrays the glacial and postglacial deposits of Connecticut in terms of their aerial extent and subsurface textural relationships. Glacial Ice-Laid Deposits (thin till, thick till, end moraine deposits) and Postglacial Deposits (alluvium, swamp deposits, marsh deposits, beach deposits, talus, and artificial fill) are differentiated from Glacial Meltwater Deposits. The meltwater deposits are further characterized using four texturally-based map units (g = gravel, sg = sand and gravel, s = sand, and f = fines). In many places a single map unit (e.g. sand) is sufficient to describe the entire meltwater section. Where more complex stratigraphic relationships exist, "stacked" map units are used to characterize the subsurface (e.g. sg/s/f - sand and gravel overlying sand overlying fines). Where postglacial deposits overlie meltwater deposits, this relationship is also described (e.g. alluvium overlying sand). Map unit definitions (Surficial Materials Polygon Code definitions, found in the metadata) provide a short description of the inferred depositional environment for each of the glacial meltwater map units. The geologic contacts between till and meltwater deposits coincide on both the Quaternary and Surficial Materials maps, as do the boundaries of polygons that define areas of thick till, alluvium, swamp deposits, marsh deposits, beach and dune deposits, talus, and artificial fill. Within the meltwater deposits, a Quaternary map unit (deposit) may contain several Surficial Materials textural units (akin to facies within a delta, for example). Combining the textural and vertical stacking information from the Surficial Materials map with the orderly portrayal of

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map California’s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of “lands” from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of California’s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bay’s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Service–Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Pigeon Point to Monterey map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these “ground-truth” surveying data are available from the CSMP Video and Photograph Portal at https://doi.org/10.5066/F7J1015K. The “seafloor character” data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The “potential habitats” polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Pigeon Point to Monterey map area data layers. Data layers are symbolized as shown on the associated map sheets.

2D Maps from CAMELS IllustrisTNG Simulations (Subset)

Author’s Note

This is my first dataset publication on Kaggle, and I’m very excited to share a small, manageable subset of the CAMELS Multifield (CMD) data to help the ML community practice image classification and astrophysical data analysis! Make sure to upvote, comment, and share if you enjoy or have suggestions for me. This subset is distributed under the MIT License with permission from the original author, Francisco Villaescusa-Navarro, and the CAMELS collaboration.

1. Overview

The CAMELS Multifield Dataset (CMD) is a massive collection of 2D maps and 3D grids derived from cosmological simulations that track the evolution of gas, dark matter, stars, black holes, and (in some suites) magnetic fields. These simulations vary important cosmological and astrophysical parameters, allowing researchers to explore and train machine learning models that could help us understand the universe’s fundamental properties.

Because the original CMD is extremely large, I’m providing a small subset from the IllustrisTNG suite, specifically from the LH (Latin Hypercube) set. This subset focuses on 2D maps at redshift z = 0.00, chosen randomly for demonstration and educational purposes.

2. About CAMELS & CMD

2.1 CAMELS Project

• CAMELS stands for Cosmology and Astrophysics with MachinE Learning Simulations.
• It aims to produce large datasets of cosmological simulations where diverse parameters (e.g., Ω_m, σ_8, supernova feedback, black-hole feedback) are systematically varied.

2.2 CMD (CAMELS Multifield Dataset)

• CMD contains hundreds of thousands of 2D maps and 3D grids from thousands of state-of-the-art hydrodynamic and N-body simulations.
• Data is organized by suite (IllustrisTNG, SIMBA, Astrid, or Nbody) and set (LH, CV, 1P, EX, BE, SB, etc.).
• Each 2D map or 3D grid is associated with a specific set of simulation parameters (the “labels”).

2.3 Suites

1. IllustrisTNG (MHD simulations): Gas, dark matter, stars, black holes, and magnetic fields.
2. SIMBA (Hydrodynamic): Similar to IllustrisTNG, but uses the SIMBA code.
3. Astrid (Hydrodynamic): Another code with its own feedback physics.
4. Nbody (Gravity-only): Follows only dark matter without astrophysical processes.

2.4 Sets

• LH (Latin Hypercube): Each simulation has unique values for all parameters, covering a broad range in parameter space.
• CV, 1P, EX, BE, SB: Other sets that vary parameters differently or keep some fixed.

3. This Subset

1. Suite & Set: IllustrisTNG – LH (Latin Hypercube sampling).
2. Field Example: Dark matter density (Mcdm), though you could encounter other fields if you download more from the official CMD resource.
3. Format: .npy files, each containing multiple 2D slices (maps).
4. Size: A small fraction of the full dataset (~5–10 maps or however many you decide to provide).
5. Coordinates: Each 2D map is a “slice” of the cosmological volume at redshift z = 0.00.
6. Parameter Labels:
   • Ω_m (matter density fraction),
   • σ_8 (root-mean-square amplitude of matter fluctuations),
   • A_SN1, A_SN2 (supernova feedback parameters),
   • A_AGN1, A_AGN2 (black-hole/AGN feedback parameters).

Because these maps are from the IllustrisTNG suite, they have non-zero values for all six parameters. If any user references the corresponding Nbody subset, note that only Ω_m and σ_8 apply there.

4. Why This Matters

The CAMELS data—especially 2D projections—are excellent for: - Machine Learning & Computer Vision: Classification, segmentation, or anomaly detection tasks. - Cosmology Research: Investigating how changes in Ω_m, σ_8, or feedback physics affect large-scale structure formation. - Educational Purposes: Students and newcomers can learn how real cosmological simulation data is structured and experiment with analysis or ML pipelines.

5. License & Attribution

This subset is shared under the MIT License granted by the original author, Francisco Villaescusa-Navarro, and the broader CAMELS collaboration. Please see the “License” section for the full text.

1. I am not the original author; I only provide a subset for educational and ML practice purposes.
2. All credits go to the CAMELS collaboration and the authors of the CMD.
3. If you publish results using this data, please cite the original sources appropriately and mention the CAMELS project.

6. References & Further Reading

-**CAMELS Project Overview**
-**IllustrisTNG Official Site**
-**[CMD Official Documentation](https://camels-multifield-dataset.read...

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map California’s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of “lands” from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of California’s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bay’s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Service–Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Offshore of Santa Cruz map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these “ground-truth” surveying data are available from the CSMP Video and Photograph Portal at https://doi.org/10.5066/F7J1015K. The “seafloor character” data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The “potential habitats” polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Offshore of Santa Cruz map area data layers. Data layers are symbolized as shown on the associated map sheets.

In this upload we share processed crop type datasets from both France and Kenya. These datasets can be helpful for testing and comparing various domain adaptation methods. The datasets are processed, used, and described in this paper: https://doi.org/10.1016/j.rse.2021.112488 (arXiv version: https://arxiv.org/pdf/2109.01246.pdf).

In summary, each point in the uploaded datasets corresponds to a particular location. The label is the crop type grown at that location in 2017. The 70 processed features are based on Sentinel-2 satellite measurements at that location in 2017. The points in the France dataset come from 11 different departments (regions) in Occitanie, France, and the points in the Kenya dataset come from 3 different regions in Western Province, Kenya. Within each dataset there are notable shifts in the distribution of the labels and in the distribution of the features between regions. Therefore, these datasets can be helpful for testing for testing and comparing methods that are designed to address such distributional shifts.

More details on the dataset and processing steps can be found in Kluger et. al. (2021). Much of the processing steps were taken to deal with Sentinel-2 measurements that were corrupted by cloud cover. For users interested in the raw multi-spectral time series data and dealing with cloud cover issues on their own (rather than using the 70 processed features provided here), the raw dataset from Kenya can be found in Yeh et. al. (2021), and the raw dataset from France can be made available upon request from the authors of this Zenodo upload.

All of the data uploaded here can be found in "CropTypeDatasetProcessed.RData". We also post the dataframes and tables within that .RData file as separate .csv files for users who do not have R. The contents of each R object (or .csv file) is described in the file "Metadata.rtf".

Preferred Citation:

-Kluger, D.M., Wang, S., Lobell, D.B., 2021. Two shifts for crop mapping: Leveraging aggregate crop statistics to improve satellite-based maps in new regions. Remote Sens. Environ. 262, 112488. https://doi.org/10.1016/j.rse.2021.112488.

-URL to this Zenodo post https://zenodo.org/record/6376160

The USGS National Hydrography Dataset (NHD) service from The National Map (TNM) is a comprehensive set of digital spatial data that encodes information about naturally occurring and constructed bodies of surface water (lakes, ponds, and reservoirs), paths through which water flows (canals, ditches, streams, and rivers), and related entities such as point features (springs, wells, stream gages, and dams). The information encoded about these features includes classification and other characteristics, delineation, geographic name, position and related measures, a "reach code" through which other information can be related to the NHD, and the direction of water flow. The network of reach codes delineating water and transported material flow allows users to trace movement in upstream and downstream directions. In addition to this geographic information, the dataset contains metadata that supports the exchange of future updates and improvements to the data. The NHD is available nationwide in two seamless datasets, one based on 1:24,000-scale maps and referred to as high resolution NHD, and the other based on 1:100,000-scale maps and referred to as medium resolution NHD. Additional selected areas in the United States are available based on larger scales, such as 1:5,000-scale or greater, and referred to as local resolution NHD. The NHD from The National Map supports many applications, such as making maps, geocoding observations, flow modeling, data maintenance, and stewardship. For additional information on the NHD, go to http://nhd.usgs.gov/index.html. The Watershed Boundary Dataset (WBD) is a companion dataset to the NHD. It defines the perimeter of drainage areas formed by the terrain and other landscape characteristics. The drainage areas are nested within each other so that a large drainage area, such as the Upper Mississippi River, will be composed of multiple smaller drainage areas, such as the Wisconsin River. Each of these smaller areas can further be subdivided into smaller and smaller drainage areas. The WBD uses six different levels in this hierarchy, with the smallest averaging about 30,000 acres. The WBD is made up of polygons nested into six levels of data respectively defined by Regions, Subregions, Basins, Subbasins, Watersheds, and Subwatersheds. For additional information on the WBD, go to http://nhd.usgs.gov/wbd.html. The National Map hydrography data is commonly combined with other data themes, such as boundaries, elevation, structures, and transportation, to produce general reference base maps. The National Map viewer allows free downloads of public domain NHD and WBD data in either Esri File or Personal Geodatabase, or Shapefile formats.

Digital data from VG2015-3 Springston, G, Thomas, E, and Kim, J, 2015,�Surficial Geologic Map of the Southern Two-Thirds of the Woodbury Quadrangle, Vermont, Washington County, Vermont: Vermont Geological Survey Open File Report VG2015-3. Data may include surficial geologic contacts, isopach contours lines, bedrock outcrop polygons, bedrock geologic contacts, hydrogeologic units and more. The surficial geologic materials data at a scale of 1:24,000 depict types of unconsolidated surficial and glacial materials overlying bedrock in Vermont. Data is created by mapping on the ground using standard geologic pace and compass techniques and/or GPS on a USGS 1:24000 topographic base map. The materials data is selected from the Vermont Geological Survey Open File Report (OFR) publication (https://dec.vermont.gov/geological-survey/publication-gis/ofr). The OFR contains more complete descriptions of map units, cross-sections, isopach maps and other information that may not be included in this digital data set.

The Watershed Boundary Dataset (WBD) from The National Map (TNM) defines the perimeter of drainage areas formed by the terrain and other landscape characteristics. The drainage areas are nested within each other so that a large drainage area, such as the Upper Mississippi River, will be composed of multiple smaller drainage areas, such as the Wisconsin River. Each of these smaller areas can further be subdivided into smaller and smaller drainage areas. The WBD uses six different levels in this hierarchy, with the smallest averaging about 30,000 acres. The WBD is made up of polygons nested into six levels of data respectively defined by Regions, Subregions, Basins, Subbasins, Watersheds, and Subwatersheds. For additional information on the WBD, go to http://nhd.usgs.gov/wbd.html. The USGS National Hydrography Dataset (NHD) service is a companion dataset to the WBD. The NHD is a comprehensive set of digital spatial data that encodes information about naturally occurring and constructed bodies of surface water (lakes, ponds, and reservoirs), paths through which water flows (canals, ditches, streams, and rivers), and related entities such as point features (springs, wells, stream gages, and dams). The information encoded about these features includes classification and other characteristics, delineation, geographic name, position and related measures, a "reach code" through which other information can be related to the NHD, and the direction of water flow. The network of reach codes delineating water and transported material flow allows users to trace movement in upstream and downstream directions. In addition to this geographic information, the dataset contains metadata that supports the exchange of future updates and improvements to the data. The NHD is available nationwide in two seamless datasets, one based on 1:24,000-scale maps and referred to as high resolution NHD, and the other based on 1:100,000-scale maps and referred to as medium resolution NHD. Additional selected areas in the United States are available based on larger scales, such as 1:5,000-scale or greater, and referred to as local resolution NHD. For more information on the NHD, go to http://nhd.usgs.gov/index.html. Hydrography data from The National Map supports many applications, such as making maps, geocoding observations, flow modeling, data maintenance, and stewardship. Hydrography data is commonly combined with other data themes, such as boundaries, elevation, structures, and transportation, to produce general reference base maps. The National Map viewer allows free downloads of public domain WBD and NHD data in either Esri File or Personal Geodatabase, or Shapefile formats.

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within Californiaâ s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map Californiaâ s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of â œlandsâ from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of Californiaâ s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bayâ s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Serviceâ Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Salt Point to Drakes Bay map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these â œground-truthâ surveying data are available from the CSMP Video and Photograph Portal at http://dx.doi.org/10.5066/F7J1015K. The â œseafloor characterâ data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The â œpotential habitatsâ polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Salt Point to Drakes Bay map area data layers. Data layers are symbolized as shown on the associated map sheets.

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The program supports a large number of coastal-zone- and ocean-management issues, including the California Marine Life Protection Act (MLPA) (California Department of Fish and Wildlife, 2008), which requires information about the distribution of ecosystems as part of the design and proposal process for the establishment of Marine Protected Areas. A focus of CSMP is to map California’s State Waters with consistent methods at a consistent scale. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. It is emphasized that the more interpretive habitat and geology data rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. This approach and CSMP planning is based in part on recommendations of the Marine Mapping Planning Workshop (Kvitek and others, 2006), attended by coastal and marine managers and scientists from around the state. That workshop established geographic priorities for a coastal mapping project and identified the need for coverage of “lands” from the shore strand line (defined as Mean Higher High Water; MHHW) out to the 3-nautical-mile (5.6-km) limit of California’s State Waters. Unfortunately, surveying the zone from MHHW out to 10-m water depth is not consistently possible using ship-based surveying methods, owing to sea state (for example, waves, wind, or currents), kelp coverage, and shallow rock outcrops. Accordingly, some of the data presented in this series commonly do not cover the zone from the shore out to 10-m depth. This data is part of a series of online U.S. Geological Survey (USGS) publications, each of which includes several map sheets, some explanatory text, and a descriptive pamphlet. Each map sheet is published as a PDF file. Geographic information system (GIS) files that contain both ESRI ArcGIS raster grids (for example, bathymetry, seafloor character) and geotiffs (for example, shaded relief) are also included for each publication. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html (last accessed September 20, 2013). The California Seafloor Mapping Program is a collaborative venture between numerous different federal and state agencies, academia, and the private sector. CSMP partners include the California Coastal Conservancy, the California Ocean Protection Council, the California Department of Fish and Wildlife, the California Geological Survey, California State University at Monterey Bay’s Seafloor Mapping Lab, Moss Landing Marine Laboratories Center for Habitat Studies, Fugro Pelagos, Pacific Gas and Electric Company, National Oceanic and Atmospheric Administration (NOAA, including National Ocean Service–Office of Coast Surveys, National Marine Sanctuaries, and National Marine Fisheries Service), U.S. Army Corps of Engineers, the Bureau of Ocean Energy Management, the National Park Service, and the U.S. Geological Survey. These web services for the Hueneme Canyon map area includes data layers that are associated to GIS and map sheets available from the USGS CSMP web page at https://walrus.wr.usgs.gov/mapping/csmp/index.html. Each published CSMP map area includes a data catalog of geographic information system (GIS) files; map sheets that contain explanatory text; and an associated descriptive pamphlet. This web service represents the available data layers for this map area. Data was combined from different sonar surveys to generate a comprehensive high-resolution bathymetry and acoustic-backscatter coverage of the map area. These data reveal a range of physiographic including exposed bedrock outcrops, large fields of sand waves, as well as many human impacts on the seafloor. To validate geological and biological interpretations of the sonar data, the U.S. Geological Survey towed a camera sled over specific offshore locations, collecting both video and photographic imagery; these “ground-truth” surveying data are available from the CSMP Video and Photograph Portal at https://doi.org/10.5066/F7J1015K. The “seafloor character” data layer shows classifications of the seafloor on the basis of depth, slope, rugosity (ruggedness), and backscatter intensity and which is further informed by the ground-truth-survey imagery. The “potential habitats” polygons are delineated on the basis of substrate type, geomorphology, seafloor process, or other attributes that may provide a habitat for a specific species or assemblage of organisms. Representative seismic-reflection profile data from the map area is also include and provides information on the subsurface stratigraphy and structure of the map area. The distribution and thickness of young sediment (deposited over the past about 21,000 years, during the most recent sea-level rise) is interpreted on the basis of the seismic-reflection data. The geologic polygons merge onshore geologic mapping (compiled from existing maps by the California Geological Survey) and new offshore geologic mapping that is based on integration of high-resolution bathymetry and backscatter imagery seafloor-sediment and rock samplesdigital camera and video imagery, and high-resolution seismic-reflection profiles. The information provided by the map sheets, pamphlet, and data catalog has a broad range of applications. High-resolution bathymetry, acoustic backscatter, ground-truth-surveying imagery, and habitat mapping all contribute to habitat characterization and ecosystem-based management by providing essential data for delineation of marine protected areas and ecosystem restoration. Many of the maps provide high-resolution baselines that will be critical for monitoring environmental change associated with climate change, coastal development, or other forcings. High-resolution bathymetry is a critical component for modeling coastal flooding caused by storms and tsunamis, as well as inundation associated with longer term sea-level rise. Seismic-reflection and bathymetric data help characterize earthquake and tsunami sources, critical for natural-hazard assessments of coastal zones. Information on sediment distribution and thickness is essential to the understanding of local and regional sediment transport, as well as the development of regional sediment-management plans. In addition, siting of any new offshore infrastructure (for example, pipelines, cables, or renewable-energy facilities) will depend on high-resolution mapping. Finally, this mapping will both stimulate and enable new scientific research and also raise public awareness of, and education about, coastal environments and issues. Web services were created using an ArcGIS service definition file. The ArcGIS REST service and OGC WMS service include all Hueneme Canyon map area data layers. Data layers are symbolized as shown on the associated map sheets.

This data package includes two related data files that can be used as input for habitat network analyses on amphibians using a specific habitat network analysis tool (HNAT; v0.1.2-alpha):

1. AmphibianHabitatNetwork_Parameters.xlsx
2. BiotopeMap_GothenburgRegion_withPondsRoadsAndBuildings.tif

HNAT is a plugin for the open-source Geographic Information System QGIS (https://qgis.org/en/site/). HNAT can be downloaded at https://github.com/SMoG-Chalmers/hnat/releases/tag/v0.1.2-alpha. To run the habitat network analyses based on the input data provided in this package one must install the plugin HNAT into QGIS. This software has been created by Chalmers within a research project financed by the Swedish government research council for sustainable development, Formas (FR -2021/0004), within the framework of the national research program "From research to implementation for a sustainable society 2021". The Excel-file contains the parameters for amphibians and the GeoTiff-file is representing a biotope raster map covering the Gothenburg region in western Sweden. SRID=3006 (Sweref99 TM). Pixel size =10x10 metres. The pixel values of the biotope map correspond to the biotope codes listed in the in the parameter file (see column “BiotopeCode”). For each biotope the parameter file holds biotope specific parameter values for two alternative amphibian models denoted “Amphibians_NMDWater_ponds” and Amphibians_NMDWater_ponds_NoFriction”. The two alternative parameter settings can be used to demonstrate the difference in model prediction with or without the assumption that amphibian movements are affected by barrier effects caused by roads, buildings and certain biotopes biotope types. The “NoFriction” version assumes that amphibian dispersal probability declines exponentially with increasing Euclidian distance whereas the other set assumes dispersal to be affected by barriers. Read the readme file for details on each parameter provided in the parameter file.

The GeoTiff-file is a biotope mape which has been created by combining a couple of publicly available geodata sets. As a base for the biotope map the Swedish land cover map NMD was used (https://geodata.naturvardsverket.se/nedladdning/marktacke/NMD2018/NMD2018_basskikt_ogeneraliserad_Sverige_v1_1.zip). To achieve a greater cartographic representation of small ponds, streams, buildings and transport infrastructure relevant for amphibian dispersal, reproduction and foraging, NMD was complemented by information from a number of vector layers. In total, 20 new biotope classes representing buildings of different height ranging from less than 5 m up to 100 m, were added to the basic land cover map. The heights were obtained by analyzing the LiDAR data provided by Swedish Land Survey (for details see Berghauser Pont et al., 2019). The data was rasterized and added on top of existing pixels representing buildings in the Swedish land cover map. The roads were separated into 101 new biotope classes with different expected number of vehicles per day. Instead of using statistics from the Swedish Transport Administration on observed number of vehicles per day relative traffic volumes were predicted based on angular betweenness centrality values calculated from the road network using PST (Place Syntax Tool, Stavroulaki et al. 2023). PST is an open-source plugin for QGIS (https://www.smog.chalmers.se/pst). Traffic volumes are expected to be correlated to the centrality values (Serra and Hillier, 2019). The vector layer with the centrality values was buffered by 15 m prior to rasterization. After that the new pixel values were added to the basic Land cover raster in sequence following the order of centrality values. Information on small streams with a maximum width of 6 m was added from a vector layer of Swedish streams (https://www.lantmateriet.se/en/geodata/geodata-products/product-list/topography-50-download-vector/). These lines where rasterized and added to the land cover raster by replacing the underlaying pixel values with new class specific pixel values. Small pondlike waterbodies was identified from the NMD data selecting contiguous fragments of the original NMD biotope class 61 with a smaller area than 1 hectare. Pixels representing the smaller water bodies was then changed to 201.

References Berghauser Pont M, Stavroulaki G, Bobkova E, et al. (2019). The spatial distribution and frequency of street, plot and building types across five European cities. Environment and Planning B: Urban analytics and city science 46(7): 1226-1242. Serra M and Hillier B (2019) Angular and Metric Distance in Road Network Analysis: A nationwide correlation study. Computers, Environment and Urban Systems 74: 194-207. Stavroulaki I, Berghauser Pont M, Fitger M, et al. (2023) PST Documentation_v.3.2.5_20231128, DOI:10.13140/RG.2.2.32984.67845.

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.

Suppose that you are provided with two items, a mapping file describing your study site of interest and a zipped folder of files that you will need to use. The mapping file describes the gridded cell centroids for each 1/16th longitude-latitude gridded cell within the study site, and the files within the folder follow the same gridded cell centroid schema in their naming conventions. The folder may contain more files than is necessary, so the mapping file will need to be used to isolate the files of interest after they have been unzipped.

These two files are intended as a tutorial materials. The Observatory for Gridded Hydrometeorology (OGH) python library contains a function called 'remapCatalog', which will be used to construct a new column in the mapping file. The column will be filled with the file path for the files that correspond to each gridded cell centroid. Files that are missing will be filled with NA.