Metadata: Identification_Information: Citation: Citation_Information: Originator: Zachery W. Fehr Originator: David G. Zawada Originator: Kimberly K. Yates Publication_Date: 20210825 Title: Seafloor Elevation and Volume Change Analyses from 2016 to 2019 Along the Florida Reef Tract, USA Geospatial_Data_Presentation_Form: Tabular, vector, and raster digital data Larger_Work_Citation: Citation_Information: Originator: Zachery W. Fehr Originator: David G. Zawada Originator: Kimberly K. Yates Publication_Date: 20210825 Title: Seafloor Elevation and Volume Change Analyses from 2016 to 2019 Along the Florida Reef Tract, USA Series_Information: Series_Name: U.S. Geological Survey data release Issue_Identification: doi:10.5066/P9CHC95D Publication_Information: Publication_Place: St. Petersburg, Florida Publisher: U.S. Geological Survey - St. Petersburg Coastal and Marine Science Center Online_Linkage: https://doi.org/10.5066/P9CHC95D Description: Abstract: The U.S. Geological Survey (USGS) St. Petersburg Coastal and Marine Science Center conducted research to quantify bathymetric changes along the Florida Reef Tract (FRT) from Miami to Marquesas Keys within a 939.4 square-kilometer area between 2016 and 2019. USGS staff used light detection and ranging (lidar)-derived data acquired by the National Oceanic and Atmospheric Administration (NOAA) during two separate lidar surveys. The first is dataset is referenced as "2016 lidar" data and was collected between July 21, 2016 and February 20, 2017 (NOAA, 2017a-c and 2018) and the second is denoted as "2019 lidar" data and was collected between November 20, 2018 and March 23, 2019 (NOAA, 2020). The 2016 and 2019 NOAA lidar datasets were used to assess changes in elevation and structure that occurred during this period along the FRT. An elevation change analysis between the 2016 and 2019 lidar data was performed to quantify and map impacts to seafloor elevation and determine elevation and volume change statistics for individual habitats found within the FRT. Data were collected under Florida Keys National Marine Sanctuary permit FKNMS-2016-068. Purpose: These data were used to determine seafloor elevation and volume changes from 2016 to 2019, throughout the Florida Reef Tract. Supplemental_Information: The 2016 lidar data were collected by the NOAA National Geodetic Survey (NGS) Remote Sensing Division using a Riegl VQ-820-G System. The 2019 lidar data were collected by Quantum Spatial, Inc. (QSI) using three Riegl systems: a Riegl VQ-880-G+, a Riegl VQ-880-GII, and a Riegl VQ-880-GH. The lidar data are an ancillary product of NOAA's Coastal Mapping Program (CMP), created through a wider Integrated Ocean and Coastal Mapping initiative to increase support for multiple uses of the data. Time_Period_of_Content: Time_Period_Information: Range_of_Dates/Times: Beginning_Date: 20160721 Ending_Date: 20190323 Currentness_Reference: ground condition Status: Progress: Complete Maintenance_and_Update_Frequency: None planned Spatial_Domain: Bounding_Coordinates: West_Bounding_Coordinate: -81.958437 East_Bounding_Coordinate: -80.086515 North_Bounding_Coordinate: 25.799164 South_Bounding_Coordinate: 24.442791 Keywords: Theme: Theme_Keyword_Thesaurus: USGS Metadata Identifier Theme_Keyword: USGS:e27b8547-ca2f-4d8c-8558-d9103ab12e36 Theme: Theme_Keyword_Thesaurus: USGS Thesaurus Theme_Keyword: marine geology Theme_Keyword: sea-floor characteristics Theme_Keyword: sea-level change Theme_Keyword: reef ecosystems Theme_Keyword: coelenterates Theme_Keyword: lidar Theme_Keyword: digital elevation models Theme: Theme_Keyword_Thesaurus: ISO 19115 Topic Category Theme_Keyword: geoscientificInformation Theme_Keyword: elevation Theme_Keyword: oceans Theme: Theme_Keyword_Thesaurus: None Theme_Keyword: seafloor elevation Theme_Keyword: sea level rise Theme_Keyword: seafloor accretion Theme_Keyword: seafloor erosion Theme_Keyword: elevation change Theme_Keyword: altimetry Theme_Keyword: submerged topography Theme: Theme_Keyword_Thesaurus: Global Change Master Directory (GCMD) Science Theme_Keyword: OCEAN > BATHYMETRY/SEAFLOOR TOPOGRAPHY > WATER DEPTH Theme_Keyword: OCEAN > COASTAL PROCESSES > COASTAL ELEVATION Theme_Keyword: OCEAN > COASTAL PROCESSES > CORAL REEFS Theme_Keyword: OCEAN > COASTAL PROCESSES > EROSION Theme_Keyword: OCEAN > COASTAL PROCESSES > SEA LEVEL RISE Theme: Theme_Keyword_Thesaurus: GCMD Providers Theme_Keyword: DOI/USGS/CMG > COASTAL AND MARINE GEOLOGY, U.S. GEOLOGICAL SURVEY, U.S. DEPARTMENT OF INTERIOR Theme: Theme_Keyword_Thesaurus: GCMD Instrument Theme_Keyword: LIDAR > LIGHT DETECTION AND RANGING Place: Place_Keyword_Thesaurus: None Place_Keyword: Florida Reef Tract Place_Keyword: Florida Keys National Marine Sanctuary Place: Place_Keyword_Thesaurus: Geographic Names Information System Place_Keyword: Florida Keys Place_Keyword: Florida Stratum: Stratum_Keyword_Thesaurus: None Stratum_Keyword: submerged Stratum_Keyword: seafloor Temporal: Temporal_Keyword_Thesaurus: None Temporal_Keyword: 2016-2019 Access_Constraints: Please note the very large files and ensure you have the adequate computing compacity to access these files. See process step 9 for more information. Use_Constraints: Public domain data from the U.S. Government are freely redistributable with proper metadata and source attribution. The U.S. Geological Survey requests to be acknowledged as originator of these data in future products or derivative research. Point_of_Contact: Contact_Information: Contact_Person_Primary: Contact_Person: Kimberly K. Yates Contact_Organization: Southeast Region: ST. PETE COASTAL & MARINE SC Contact_Position: Research Oceanographer Contact_Address: Address_Type: mailing and physical Address: 600 4Th Street South City: St. Petersburg State_or_Province: FL Postal_Code: 33701 Country: US Contact_Voice_Telephone: 727-502-8059 Contact_Electronic_Mail_Address: kyates@usgs.gov Cross_Reference: Citation_Information: Originator: Yates, Kimberly K. Originator: Zawada, David G. Originator: Smiley, Nathan A. Originator: Tiling-Range, Ginger Publication_Date: 20170420 Title: Divergence of seafloor elevation and sea level rise in coral reef ecosystems Publication_Information: Publication_Place: Munich, Germany Publisher: Biogeosciences Online_Linkage: https://doi.org/10.5194/bg-14-1739-2017 Cross_Reference: Citation_Information: Originator: National Oceanic and Atmospheric Administration Office for Coastal Management Publication_Date: 20170914 Title: 2016 NOAA NGS Topobathy Lidar DEM: Florida Keys Outer Reef Block 01 Publication_Information: Publication_Place: Charleston, South Carolina Publisher: NOAA Other_Citation_Details: 2017a Online_Linkage: https://www.fisheries.noaa.gov/inport/item/48373 Cross_Reference: Citation_Information: Originator: National Oceanic and Atmospheric Administration Office for Coastal Management Publication_Date: 20170614 Title: 2016 NOAA NGS Topobathy Lidar DEM: Florida Keys Outer Reef Block 02 Publication_Information: Publication_Place: Charleston, South Carolina Publisher: NOAA Other_Citation_Details: 2017b Online_Linkage: https://www.fisheries.noaa.gov/inport/item/49423 Cross_Reference: Citation_Information: Originator: National Oceanic and Atmospheric Administration Office for Coastal Management Publication_Date: 20170712 Title: 2016 NOAA NGS Topobathy Lidar DEM: Florida Keys Outer Reef Block 03 Publication_Information: Publication_Place: Charleston, South Carolina Publisher: NOAA Other_Citation_Details: 2017c Online_Linkage: https://www.fisheries.noaa.gov/inport/item/48374 Cross_Reference: Citation_Information: Originator: National Oceanic and Atmospheric Administration Office for Coastal Management Publication_Date: 20180219 Title: 2017 NOAA NGS Topobathy Lidar DEM: Florida Keys Outer Reef Block 04 Publication_Information: Publication_Place: Charleston, South Carolina Publisher: NOAA Online_Linkage: https://www.fisheries.noaa.gov/inport/item/51529 Cross_Reference: Citation_Information: Originator: National Oceanic and Atmospheric Administration Office for Coastal Management Publication_Date: 20201002 Title: 2018-2019 NOAA NGS Topobathy Lidar DEM Hurricane Irma: Miami to Marquesas Keys, FL Publication_Information: Publication_Place: Charleston, South Carolina Publisher: NOAA Online_Linkage: https://www.fisheries.noaa.gov/inport/item/63018 Data_Quality_Information: Attribute_Accuracy: Attribute_Accuracy_Report: Datasets were visually inspected by USGS staff in Esri ArcGIS Desktop Advanced version 10.7, ArcGIS Pro version 2.4.0, and Global Mapper version 20.1 for identification of anomalous elevations or data inconsistencies. Logical_Consistency_Report: Data cover area specified for this project without known issues. Completeness_Report: Dataset is considered complete for the information presented, as described in the abstract. Users are advised to carefully read the rest of the metadata record and refer to Yates and others (2017) for additional details. Positional_Accuracy: Horizontal_Positional_Accuracy: Horizontal_Positional_Accuracy_Report: For the 2016 lidar data, the elevation positions were obtained using post-processed kinematic global positioning system (KGPS) methods. The horizontal accuracy of the 2016 lidar data is better than plus or minus 1.0 meters (m); Quantitative Value: 1.0 m. The 2019 lidar data have been reported horizontal accuracy to be 0.115 m at the 95% confidence level using the National Standards for Spatial Data Accuracy (NSSDA) reporting standard. Vertical_Positional_Accuracy: Vertical_Positional_Accuracy_Report: For the 2016 lidar, the data positions were obtained using post-processed KGPS methods. Data used to validate the lidar were collected with static GPS observational equipment and compared against the published data. For the 2019 lidar data, submerged topography accuracy tested 0.110 m vertical accuracy at the 95% confidence level against the classified points cloud using NSSDA reporting standards. Lineage: Process_Step: Process_Description: Step 1: Seafloor elevation and volume change analyses were performed using methods from Yates and others (2017). Original lidar data were downloaded from NOAA’s DigitalCoast website using the site’s 'Customized Download' capability, at https://coast.noaa.gov/dataviewer/#/lidar/search/. The full extent of the 2016 data in the form of four blocks were downloaded as well as the full extent of the 2019 data which is available as a single file. The data were downloaded using the following parameters: UTM zone: Zone 17 Range 084W-078W; Horizontal Datum: NAD83; Horizontal Units: Meters; Vertical Datum: NAVD88; Vertical Units: Meters; File Format: Tiff 32-bit Float; Bin Method: TIN; Bin Size: 1.0; Bin Units: Meters; Data Classification: Bathymetric Lidar Points; Data Returns: Any Points; Ancillary Data: No Ancillary Data; and Geoid Name: GEOID12B. Using VDatum version 3.9, a publicly available software from NOAA (https://vdatum.noaa.gov/), the lidar tagged image file format (TIFF) files for both 2016 and 2019 datasets were transformed from the North American Horizontal Datum of 1983 (NAD83) Universal Transverse Mercator (UTM) Zone 17 North to NAD83 National Spatial Reference System (NSRS2011) horizontal datum, and the vertical datum (North American Vertical Datum of 1988 [NAVD88]) and geoid model (GEOID12B) were kept the same. Process_Date: 2020 Process_Step: Process_Description: Step 2: The 2016 lidar data (blocks 01-04) were merged into a single TIFF using the “Mosaic To New Raster (Data Management)” tool in ArcMap 10.7. Blocks were merged by prioritizing the most recently collected block in areas of overlap. This was accomplished by utilizing the following parameters: blocks 01-04 were used as input rasters and were loaded in numerical order; Output Location was set to a personal workspace; Raster Dataset Name with Extension: 2016_Data_Mosaic_Last.tif; Pixel Type (optional): 32_BIT_FLOAT; Number of Bands: 1; and Mosaic Operator (optional); LAST. All other parameters were left as default. Process_Date: 2020 Process_Step: Process_Description: Step 3: Footprints of both the 2016 and 2019 lidar data were then created by opening the Digital Elevation Models (DEMs) in ArcMap. The “Reclassify (Spatial Analyst)” tool was used on both DEMs by replacing all old values with ‘1’ and keeping ‘No Data’ values as ‘No Data.’ The resulting TIFF was converted into a footprint using the “Raster to Polygon (Conversion)” tool to convert the raster files to polygon shapefiles. An intersect footprint was created between the 2016 and 2019 lidar data using the “Intersect (Analysis)” tool. The output polygon (2016_2019_lidar_footprint.shp) represents the shared coverage of 2016 and 2019 DEMs. Both DEMs were then clipped to the extent of the intersect footprint using the “Clip (Data Management) tool. Process_Date: 2020 Process_Step: Process_Description: Step 4: The clipped DEMs were then edited for removal of reduced quality data. Using Global Mapper 20.1, each DEM was visually inspected using the methods of Yates and others (2017) to identify areas characterized by coarse interpolation relative to surrounding areas and appearing as small groups of large triangles created during the development of each DEM. Polygons were manually drawn to encompass the areas for removal using the “Digitizer” tool and the ‘Create Area/Polygon Features’ function. By default, these polygons are generated in the ‘User Created Features’ layer within the Control Center. The ‘User Created Features’ layer was then exported as a polygon shapefile. The 2016 data has four polygon shapefiles (Block01_2016lidar_RemovedPolygons.shp, Block02_2016lidar_RemovedPolygons.shp, Block03_2016lidar_RemovedPolygons.shp, and Block04_2016lidar_RemovedPolygons.shp) and the 2019 data has one polygon shapefile (RemovedPolygons_2019lidar.shp). The “Erase (Analysis)” tool in ArcMap was used to remove the areas of reduced quality from the 2016_2019_Intersect_Footprint.shp shapefile creating the 2016_2019_edited_lidar_footprint_intersect.shp shapefile. Using the “Clip (Data Management)” tool, areas of reduced quality data were removed from the 2016 and 2019 lidar DEMs by clipping each DEM to the 2016_2019_edited_lidar_footprint_intersect.shp shapefile. Process_Date: 2020 Process_Step: Process_Description: Step 5: Data limitations exist in ArcMap 10.7 for elevation change analyses performed at this geographic extent (around 1 billion square meters). Approximately 250 million data points were required to sample the 2016 and 2019 DEMs using a 2-m grid spacing. However, point shapefiles were limited to approximately 80 million data points in ArcMap 10.7. DEMs were divided into five separate blocks that were used to generate five sets of 2-m grid shapefiles containing approximately 50 million points each to reduce the time required for elevation change computations. A polygon shapefile was created using the ‘Create Features’ function in ArcMap to cover roughly one fifth of the DEM beginning at the northeastern extent and was named Block01_footprint.shp. The DEMs were then clipped to the extent of the Block01_footprint.shp shapefile using the “Clip (Data Management)” tool. The symbology panel was opened for the resulting TIFF and classification statistics were inspected to determine if the resulting coverage was between 170 and 230 million cells, the Block01_footprint.shp polygon shapefile redrawn and the classification statistics were reinspected until the target range was met. A buffer was added to the block using a clipped portion of the intersect footprint shapefile that covered the extent of the block boundary with a few hundred meters of coverage on either side of the boundary. The clipped portion was added to the Block01_footprint.shp shapefile using the “Merge (Data Management)” tool. For Block 01, this buffer was added to the southern boundary, all other blocks had their buffer added to the northern boundary. The resulting polygon was clipped to the 2016_2019_edited_lidar_footprint_intersect.shp shapefile creating the Block01_TINextensionArea.shp shapefile. The buffer was required to prevent block edge artifacts in final volume change analysis steps requiring merging of elevation change TIN models for adjacent data blocks. Moving southwest along the DEM, another arbitrary polygon was drawn beginning at the boundary of the Block01_footprint.shp polygon covering the next fifth of the DEM as well as a slight overlap with the Block01_footprint.shp shapefile. The overlap extent of the Block01_footprint.shp shapefile was removed from the newly created polygon using the “Erase (Analysis)” tool in ArcMap creating the Block02.shp shapefile. The DEM was then clipped to the extent of the Block02.shp shapefile and classification statistics were inspected for the resulting TIFF to confirm the target range cell count. Following cell count target verification, the Block02.shp shapefile was clipped to the extent of the intersect footprint shapefile using the “Clip (Analysis)” tool resulting in the Block02_footprint.shp shapefile. A second version with a buffer was created using the same method used with the Block01 shapefiles, with the buffer being added to the northern boundary, creating the Block02_TINextensionArea.shp shapefile. This process was repeated until full coverage of the DEM was achieved, and 10 shapefiles were created: five polygons without a buffer, and five polygons with a buffer. This method resulted in a polygon overlap at the intersection of the polygons with buffers and full coverage with no overlap with the non-buffered polygon version. Process_Date: 2020 Process_Step: Process_Description: Step 6: 2-m grids were generated for every block’s ‘TINextensionArea.shp’ shapefile in Global Mapper Version 20.1. Each 2-m grid was generated for each block shapefile by selecting the polygon with the “Digitizer” tool, and right clicking on the polygon within the frame and selecting “GRID – Create Regular Grid of User-Specified Size/Orientation…” from the “Advanced Feature Creation Options.” In the Grid Setup window, the following parameters were used: Grid Cell Width and Grid Cell Height were set to 2; the Calculate Grid Cell Counts to Fill Rectangle Option was selected; Crop to Selected Feature(s) was selected within the Grid Bounds window; Crop Generated Grid to Selected Area Feature was enabled; Keep Area if Any Part in Crop Area was enabled; Create Points at Grid Cell Centers was enabled; and all other options were left as default. The generated 2-m grid was then manually shifted to center grid points on the DEM pixels. The 2-m grid was selected using the “Digitizer” tool in Global Mapper, and in the Move/Reshape Feature(s) window, the Shift (Offset) Selected Feature(s) tab was selected to open the Specify Offset to Apply to Point(s) menu and the following parameters were set: Units: meters; X/Longitude Offset: +0.5; and Y/Latitude Offset: +0.5. Grid points that were shifted into no-data areas were removed to avoid extraction of no-data values. The relevant ‘TINextensionArea.shp’ polygon shapefile was selected using the Digitizer tool and the 2-m grid “Options” menu was opened by right clicking the 2-m grid layer in the “Control Center,” and “Export Layers to New File(s)” was selected from the “Layer” window. Within the “Select Export Format” menu, “Data Format” was set to Shapefile. The “Export Bounds” tab was selected, and the “Crop to Selected Area Features” parameter was enabled in the “Shapefile Export Options” window to remove grid points that were in no-data areas during shapefile export. The resulting shapefile was reimported into Global Mapper for the addition of horizontal coordinates. Horizontal coordinates were added to each grid point (based on the projection of the data frame) by selecting the shapefile with the “Digitizer” tool and “Add Coordinates/Bounds Attributes to Selected Features” was selected from the “Attribute/Style Functions” window. Process_Date: 2020 Process_Step: Process_Description: Step 7: Elevation values were extracted from the 2016 and 2019 DEMs at the location of the 2-m grid points generated within each block’s ‘TINextensionArea.shp’ shapefile. The 2016 DEM was opened in Global Mapper and the resampling method was changed by opening the “Options” menu of the DEM and “Resampling” was set to No Resampling (Nearest Neighbor). The block’s 2-m grid was selected using the “Digitizer” tool. Within the “Attributes/Style Functions” window, “Apply Elevations to Selected Feature(s)” was selected, with “Terrain Layers” as the only selected parameter. An attribute named ELEVATION was created containing the elevation values extracted from the 2016 DEM. Using the “Digitizer” tool, the “Edit Selected Features” window was selected to rename the ELEVATION attribute to ELEV2016. The same steps were repeated to extract elevation values from the 2019 DEM at the location of grid points in each block’s shapefile and to rename the ELEVATION attribute containing the 2019 extracted elevation values to ELEV2019. Extracted elevation values are based on the source elevation so these data are in NAD83 (NSRS2011) UTM17N and meters. Process_Date: 2020 Process_Step: Process_Description: Step 8: The elevation-difference between 2016 and 2019 extracted elevation values in each block was calculated using the “Digitizer” tool. The corresponding 2-m grid was selected for each block and “Calculate/Copy Attributes for Feature Selection” was selected from the “Attribute/Style Functions” window. The elevation-difference was calculated using the following parameters: “Select Existing or Create New Attribute to Assign Calculated Values”: Diff_m; “Source Attribute”: ELEV2019; “Operation”: Subtract; and “Use Attribute Value”: ELEV2016. ELEV2019 represents the modern elevation values and ELEV2016 represents the historical elevation values. The final 2-m grid for each block was exported in two separate formats including one shapefile using the block’s buffered footprint as an export bound (for example, Block01_TINextensionArea.shp polygon), and a second shapefile using the block’s original footprint as an export bound (for example, Block01_footprint.shp polygon). The “Export Layers to New File(s)” function with “Data Format” set to Shapefile and “Crop to Selected Area Feature(s)” was selected in the Export Bounds tab. Each block’s shapefile exported using the buffer export bound was named with respect to the Block it was generated in, for example, Block01_ElevationChangePoints_TINExtension.shp. Each block’s shapefile exported using the original block footprint polygon was named in the same format, for example, Block01_ElevationChangePoints.shp. Steps 6 through 8 were repeated for each block, with one notable exception. When shifting the points for Blocks 02 through 05, the shift values were determined by matching any point generated within the buffer zone to its nearest corresponding point in the adjacent block’s buffer zone in the buffered version of the grid shapefiles. Point generation in these blocks were based on extent bounds and therefore did not use a simple half-meter shift. Latitude and longitude values were noted for two points. For example, one point was chosen within the buffer zone of the Block01_TINextensionArea.shp shapefile and the coordinates of that point were documented, the nearest point that results from the Step 6 process step for Block02 also had its coordinates documented. The difference in X and Y coordinates were used to determine the horizontal and vertical shift to be used for that dataset. Process_Date: 2020 Process_Step: Process_Description: Step 9: Elevation change surface models were created in ArcGIS Pro version 2.4.0 using the calculated elevation-difference (Diff_m) points from the various ‘TINExtension.shp’ point shapefiles. The five shapefiles used for elevation surface model generation are as follows: Block01_ElevationChangePoints_TINExtension.shp, Block02_ElevationChangePoints_TINExtension.shp, Block03_ElevationChangePoints_TINExtension.shp, Block04_ElevationChangePoitns_TINExtension.shp, and Block05_ElevationChangePoints_TINExtension.shp. Due to data size limitations in ArcMap, Triangulated Irregular Network (TIN) generation is only possible using ArcGIS Pro 2.1.3 or lRWE at this time. Any attempt to generate a TIN in ArcMap version 10.7 or earlier will result in a CPU memory error. Beginning with Block 01, the Block01_ElevationChangePoints.shp shapefile was loaded into ArcGIS Pro version 2.4.0. A TIN was created from the Diff_m points using the "Create TIN (3D Analyst)" tool by specifying the Block01_ElevationChangePoints.shp shapefile as the "Input Feature Class", Diff_m as the "Height Field" and Mass_Points as the "Type", creating the Block01_Extension_TIN file. The Block01_Extension_TIN was then delineated using the “Delineate TIN Data Area (3D Analyst)” tool by specifying the Block01_Extension_TIN as the “Input TIN”, a "Maximum Edge Length" of 2.828428 (hypotenuse of a triangle with 2-m legs) and the "Method" set to ALL EDGES. It was noted that the resultant delineated TIN needed to be reloaded in a new application for changes to be applied. The tool was allowed a few minutes post completion and a new ArcGIS Pro application was used to import the delineated TIN. The delineated Block01_Extension_TIN file was then clipped to the extent of the 2016 and 2019 DEMs using the “Edit TIN (3D Analyst)” tool with the following parameters: "Input TIN": Block01_Extension_TIN file; "Input Features Class": Block01_TINextensionArea.shp; "Height Field": None; "Tag Field": None; and "Type": Hard clip. As with the delineation step, the file needed to be reloaded in a new application after tool completion for changes to be applied. The Block01_Extension_TIN file was then clipped again using clipped versions of the habitat shapefile (outlined in Step 10) under the following parameters: "Input TIN": Block01_Extension_TIN file; "Input Features Class": Block01_HabMapClip.shp; "Height Field": None; "Tag Field": None; and "Type": Hard clip. An application reset was used after this tool as well. The final TIN was exported as the Block01_ProcessedTIN file. This process was repeated for each ‘TINExtension.shp’ point shapefile, creating five TINs covering the extent of the study area. Please note the very large files and ensure you have the adequate computing compacity to access these files. These TIN files were generated and can be viewed using ArcGIS Pro version 2.1.3 or later to accommodate the large file sizes. They may not be accessible in other geoprocessing software with file size limits. Smaller TIN files can be generated for viewing in other programs using the following methods: generate a polygon shapefile covering a data extent that does not exceed available computing capabilities. Generate a point shapefile that covers the extent of polygon shapefile using the “Clip (Analysis)” tool and relevant elevation-change point shapefile (ElevationChangeShapefiles.zip). Follow process steps 7 through 9 for developing a TIN at this smaller data extent Process_Date: 2020 Process_Step: Process_Description: Step 10: The original Unified Florida Reef Tract Map version 2.0 shapefile was downloaded from https://ocean.floridamarine.org/IntegratedReefMap/UnifiedReefTract.htm. Using ArcMap, the original habitat shapefile was modified using the "Clip (Analysis)" tool to clip the habitat shapefile to the extent of the 2016 and 2019 DEMs by specifying the habitat shapefile as the "Input Features" and the 2016_2019_edited_lidar_footprint_intersect.shp as the "Clip Features", creating the UR_HabMap_Clip.shp shapefile. The shapefile was then loaded into Global Mapper and exported using the 2011 realization of NAD83 as UR_Habmap_2011_GM_Export.shp. The DEMs used in this analysis cover areas that the Unified Reef map do not, so an additional polygon was integrated into the shapefile. Using the "Erase (Analysis)" tool in ArcMap, the 2016_2019_edited_lidar_footprint_intersect.shp shapefile was used as the "Input Features" and the UR_HabMap_2011_GM_Export.shp shapefile was used as the "Erase Features" creating the Unclassified_Polygons.shp shapefile. The Unclassified_Polygons.shp shapefile was edited in an edit session to add the "ClassLv2" attribute and assign its value as "Unclassified.” The edit session was closed, and the changes were saved. Using the "Merge (Data Management)" tool, the Unclassifed_Polygons.shp shapefile and the UR_HabMap_2011_GM_Export.shp shapefile were used as the "Input Datasets" creating the UR_HabMap_2011_GM_Export_UnclassifiedAdded.shp shapefile. Using the "Select by Attribute" tool, 17 individual habitat shapefiles were created from the UR_HabMap_2011_GM_Export_UnclassifiedAdded.shp shapefile by selecting one ClassLv2 habitat and exporting it as a separate shapefile. For volume analyses, this habitat file was split into five sections by using the “Clip (Analysis)” tool in ArcMap and clipping the UR_Habmap_2011_GM_Export_UnclassifiedAdded.shp file to each of the individual block footprints (Block01_footprint.shp, Block02_footprint.shp, Block03_footprint.shp, Block04_footprint.shp, and Block05_footprint.shp) creating the following files: Block01_HabMapClip.shp, Block02_HabMapClip.shp, Block03_HabMapClip.shp, Block04_HabMapClip.shp, and Block05_HabMapClip.shp. Process_Date: 2020 Process_Step: Process_Description: Step 11: Five sets of elevation change statistics were determined by habitat type for each block using the elevation-difference points from the following five shapefiles: Block01_ElevationChangePoints.shp, Block02_ElevationChangePoints.shp, Block03_ElevationChangePoints.shp, Block04_ElevationChangePoints.shp, and Block05_ChangeElevationPoints.shp. In ArcGIS Pro, the "Select Layer by Location (Data Management)" tool was used to extract points within or on the boundary of a specific habitat type by using the following parameters: "Input Feature Layer": Block01_ElevationChangePoints.shp; "Relationship": INTERSECT; "Selecting Features": Clipped Habitat shapefile; "Search Distance": left blank; and "Selection Type": NEW_SELECTION. An ArcGIS Pro model named Seafloor Elevation Change Analysis Tool (SECAT) was created to automate the process as these steps had to be repeated for 17 habitat types across each block's shapefile. Elevation change statistics were compiled by habitat type into a comma separated values (CSV) file using Microsoft Excel 2016, see Block01_Elevation.csv, Block02_Elevation.csv, Block03_Elevation.csv, Block04_Elevation.csv, Block05_Elevation, and 2016_2019_FRT_CollectiveElevationStatistics.csv included in SECAT_Analyses.zip. Excel sheet versions (XSLX) of the collective statistics are also included for convenience. Process_Date: 2020 Process_Step: Process_Description: Step 12: Using ArcGIS Pro, volume change statistics per habitat type were calculated utilizing each of the TINs generated in Step 9. Surface volume changes were calculated for six cases using the "Surface Volume (3D Analyst)" tool. To calculate the net erosion lower limit (case 1) the "Reference Plane" was set to BELOW and the "Plane Height" was set to -0.19 m. For the net erosion upper limit (case 2) the "Reference Plane" was set to BELOW and the "Plane Height" set to 0 m. For the net accretion lower limit (case 3) the "Reference Plane" was set to ABOVE and the "Plane Height" was set to 0.19 m. For the net accretion upper limit (case 4) the "Reference Plane" was set to ABOVE and the "Plane Height" was set to 0 m. A 0.19 m threshold was determined by vertical error analysis using the uncertainties reported in the metadata of the original 2016 (0.15 m) and 2019 (0.11 m) DEMs to calculate the Root Mean Square Error (RMSE) of 0.19 m (Yates and others, 2017). An additional threshold of 0.31 m (1.65 x RMSE) was used as a conservative threshold. The "Reference Plane" was set to BELOW and the "Plane Height" was set to -0.31 m for the net erosion lower limit at this threshold (case 5). The "Reference Plane" was set to ABOVE and the "Plane Height" was set to 0.31 m for the net accretion lower limit at this threshold (case 6). Minimum net volume change was calculated by summing results from cases 1 and 3 or 1 and 5 depending on lower limit threshold. Maximum net volume change was calculated by summing results from cases 2 and 4 or 2 and 6 depending on lower limit threshold. The area normalized volume change lower limit was calculated by dividing the minimum net volume change for each habitat by the habitat's total area. The area normalized volume change upper limit was calculated by dividing the maximum net volume for each habitat by the habitat's total area. An ArcGIS Pro model was created to automate the process, as these steps had to be repeated for 17 habitat types. Volume change statistics were generated for each individual block and the results were then compiled by habitat type in CSV format using Excel, see 2016_2019_FRT_CollectiveVolumeStatistics_019.csv and 2016_2019_FRT_CollectiveVolumeStatistics_031.csv as well as the individual block statistics at each threshold which are also included in SECAT_Analyses.zip. Excel sheet versions (XLSX) of the collective statistics are also included for convenience. Process_Date: 2020 Spatial_Data_Organization_Information: Direct_Spatial_Reference_Method: Point Spatial_Reference_Information: Horizontal_Coordinate_System_Definition: Planar: Grid_Coordinate_System: Grid_Coordinate_System_Name: Universal Transverse Mercator Universal_Transverse_Mercator: UTM_Zone_Number: 17 Transverse_Mercator: Scale_Factor_at_Central_Meridian: 0.9996 Longitude_of_Central_Meridian: -81.0 Latitude_of_Projection_Origin: 0.0 False_Easting: 500000.0 False_Northing: 0.0 Planar_Coordinate_Information: Planar_Coordinate_Encoding_Method: coordinate pair Coordinate_Representation: Abscissa_Resolution: 0.6096 Ordinate_Resolution: 0.6906 Planar_Distance_Units: METERS Geodetic_Model: Horizontal_Datum_Name: North American Datum of 1983 National Spatial Reference System (2011) Ellipsoid_Name: GRS_1980 Semi-major_Axis: 6378137.0 Denominator_of_Flattening_Ratio: 298.257222 Vertical_Coordinate_System_Definition: Depth_System_Definition: Depth_Datum_Name: North American Vertical Datum of 1988 (NAVD88) GEOID12B Depth_Resolution: 0.2 Depth_Distance_Units: meters Depth_Encoding_Method: Explicit depth coordinate included with horizontal coordinates Entity_and_Attribute_Information: Overview_Description: Entity_and_Attribute_Overview: The detailed attribute descriptions for the elevation and volume change statistics are provided in the included data dictionaries (DataDictionary_ElevationChange.pdf, DataDictionary_HabitatTypes.pdf, and DataDictionary_VolumeChange.pdf). These metadata are not complete without these files. Entity_and_Attribute_Detail_Citation: The entity and attribute information were generated by the individual and/or agency identified as the originator of the dataset. Please review the rest of the metadata record for additional details and information. Distribution_Information: Distributor: Contact_Information: Contact_Person_Primary: Contact_Person: Kimberly K. Yates Contact_Organization: Southeast Region: ST. PETE COASTAL & MARINE SC Contact_Position: Research Oceanographer Contact_Address: Address_Type: mailing and physical Address: 600 4Th Street South City: St. Petersburg State_or_Province: FL Postal_Code: 33701 Country: US Contact_Voice_Telephone: 727-502-8059 Contact_Electronic_Mail_Address: kyates@usgs.gov Distribution_Liability: Although these data have been processed successfully on a computer system at the U.S. Geological Survey (USGS), no warranty expressed or implied is made regarding the display or utility of the data on any other system, or for general or scientific purposes, nor shall the act of distribution constitute any such warranty. The USGS shall not be held liable for improper or incorrect use of the data described or contained herein. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Standard_Order_Process: Digital_Form: Digital_Transfer_Information: Format_Name: JPG, SHP, CSV, XLSX, TIFF Format_Version_Number: None, ArcGIS 10.7, RFC 4180, ArcGIS Pro 2.4.0, Global Mapper 20.1 Format_Specification: Joint photographic group, Esri point and polygon shapefiles, comma-separated values, Microsoft Excel spreadsheet, and tagged image file format Digital_Transfer_Option: Online_Option: Computer_Contact_Information: Network_Address: Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/ElevationChangeShapefiles.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/VolumeChangeShapefiles.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/2016_2019_LidarFiles.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/HabitatClips.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/2016_2019_ProcessedTINs.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/ElevationChangeFigures.zip Network_Resource_Name: https://coastal.er.usgs.gov/data-release/doi-P9CHC95D/data/SECAT_Analyses.zip Fees: None Technical_Prerequisites: Users are recommended to utilize ArcGIS Pro version 2.1.3 or later and Global Mapper version 20.1 or later when attempting to access geospatial files in this data release. Metadata_Reference_Information: Metadata_Date: 20210825 Metadata_Contact: Contact_Information: Contact_Person_Primary: Contact_Person: Kimberly K. Yates Contact_Organization: Southeast Region: ST. PETE COASTAL & MARINE SC Contact_Position: Research Oceanographer Contact_Address: Address_Type: mailing and physical Address: 600 4Th Street South City: St. Petersburg State_or_Province: FL Postal_Code: 33701 Country: US Contact_Voice_Telephone: 727-502-8059 Contact_Electronic_Mail_Address: kyates@usgs.gov Metadata_Standard_Name: Content Standard for Digital Geospatial Metadata Metadata_Standard_Version: FGDC-STD-001-1998