Projected Seafloor Elevation Along the Florida Reef Tract From Port St. Lucie to Marquesas Key, Florida-50 Years From 2001 Based on Historical Rates of Mean Erosion

Step 1: The original U.S. Coastal Relief Model Vol.3 - Florida and East Gulf of Mexico DEM was downloaded from https://data.nodc.noaa.gov/cgi-bin/iso?id=gov.noaa.ngdc.mgg.dem:307. The network common data format (netCDF) DEM file was converted to a tagged image file format (TIFF) using Blue Marble Global Mapper version 19.1.0. All remaining steps were completed with Esri ArcGIS Desktop Advanced version 10.6. Using the original TIFF, a footprint of the DEM was created with the "Reclassify (Spatial Analyst)" tool in ArcToolbox by replacing all old data values with 1 and the "No Data" value with 0 to create a raster file. Then, the "Raster to Polygon (Conversion)" tool was used to create a footprint of the original Florida Reef Tract DEM by converting the raster file to a polygon shapefile (SHP).

Step 2: The original Unified Florida Reef Tract Map version 2.0 polygon SHP file was downloaded from http://ocean.floridamarine.org/IntegratedReefMap/UnifiedReefTract.htm. Using Esri ArcGIS, the Artificial, Land, and Mangrove areas (ClassLv2) were removed from the habitat SHP file using the "Select by Attribute" tool to select the three habitat types, and then the Editor Toolbox to delete them.

Step 3: The original Florida Reef Tract DEM TIFF (from Step 1), was clipped to the extent of the modified habitat SHP file (from Step 2) using the "Clip (Data Management)" tool in Esri ArcGIS by specifying the Florida Reef Tract DEM TIFF as the 'Input Raster' and the modified habitat SHP file as the 'Output Extent'. The clipped TIFF was used to extract two shoreline contours with the "Contour List (Spatial Analyst)" tool by adding 'Contour values' for 0-m, 0.05 m and -0.4 m contours. The shoreline contour SHP file was manually connected across inlets and gaps using the "Straight Segment" tool in the Editor Toolbox to draw straight lines between the seaward most points of channels or other breaks along the coastline. The backshore was removed from the contour SHP file by using the "Buffer (Analysis)" tool with a 'Distance Linear unit' of 1000 m. Then, the shoreline contour was smoothed with the "Smooth Line (Cartography)" tool using the Polynomial Approximation with Exponential Kernel (PAEK) smoothing algorithm and a 50-m smoothing tolerance. Using the extended contour SHP file, a polygon SHP file was created using the "Feature to Polygon (Data Management)" tool by specifying the contour SHP file as the 'Input Features', to create the Coastal Clip SHP file. The clipped DEM TIFF was modified further using the "Clip (Data Management)" tool to remove coastal land areas by adding the DEM TIFF as the 'Input Raster' and the Coastal Clip SHP file as the 'Output Extent', creating the FloridaReefTract_ElevationSurface_OriginalClip TIFF.

Step 4: The clipped habitat SHP file (from Step 2) was modified by removing coastal land areas using the "Clip (Analysis)" tool with the habitat SHP file as the 'Input Features' and the Coastal Clip SHP file (from Step 3) as the 'Clip Features', creating the FloridaReefTract_Habitat_Clip SHP file.

Step 5: Mean future elevation for each habitat type in the FloridaReefTract_ElevationSurface_OriginalClip DEM was calculated by applying previously published mean erosion values in the Upper Florida Keys (UFK) that were projected 50 years into the future from the year 2002 (Yates and others, 2018, https://doi.org/10.5066/P9CI9LNH). Mean projected erosion rates in the UFK were compiled by habitat type into a comma separated values (CSV) file using Microsoft Excel 2016.

Step 6: Using the FloridaReefTract_Habitat_Clip SHP file in Esri ArcGIS, individual polygon SHP files were created for each habitat type using the "Select by Attribute" tool and exporting each habitat type as a separate SHP file. Individual DEMs were created from the FloridaReefTract_ElevationSurface_OriginalClip TIFF using the "Clip (Data Management)" tool to clip the full DEM to the extent of each habitat type SHP file by specifying the FloridaReefTract_ElevationSurface_OriginalClip as the 'Input Raster' and the habitat type SHP file as the 'Output Extent'. The "Raster Calculator (Spatial Analyst)" tool was used to project future erosion per habitat type by modifying individual habitat DEMs by adding or subtracting the corresponding 'Mean projected erosion 50 years from 2002 in the Upper Florida Keys (m/50 years)' from the 50_Year_FloridaReefTract_Seafloor_Projection_MeanErosion CSV file. Individual habitat DEMs were merged together using the "Mosaic to New Raster (Data Management)" tool with the 'Pixel Type' set to 32_BIT_FLOAT, 'Number of Bands' set to 1 and the 'Mosaic Operator' set to MEAN to create the final 50_Year_FloridaReefTract_Seafloor_Projection_DEM_MeanErosion TIFF file.

Added keywords section with USGS persistent identifier as theme keyword.