Attribute_Accuracy_Report:
The accuracies of the four derived thematic maps for geomorphic setting, substrate type, vegetation type, and vegetation density were assessed for 2014–2015 conditions at 15 study areas, including Rockaway Peninsula. These maps were derived in the same manner as those distributed here (derived for Rockaway Peninsula based on 2012 conditions). We assume that accuracy of maps of 2014–2015 conditions reflect the accuracy of our methods more generally, and accuracy values are therefore relevant for the files distributed here (Rock12_SubType.tif, Rock12_VegType.tif, Rock12_VegDen.tif, Rock12_GeoSet.tif).
Accuracy of the 2014–2015 derived raster layers at 15 study sites was assessed by comparing them against 482-501 ground validation points. The number of validation points used depended on the dataset examined. Validation points were subset from the iPlover dataset (Sturdivant and others 2016). We followed the methodology recommended by the National Park Service (Lea and others 2010) to produce statistics for overall accuracy, chance agreement, and a kappa coefficient. Classifications were not altered when a difference between a ground validation point and the underlying raster cell was found.
For geomorphic setting (501 validation points), overall accuracy was 46% and chance agreement was 23%, resulting in a kappa coefficient of 0.30. For substrate type (490 validation points), overall accuracy was 96% and chance agreement was 51%, resulting in a kappa coefficient of 0.92. For vegetation type (482 validation points), overall accuracy was 93% and chance agreement was 46%, resulting in a kappa coefficient of 0.88. For vegetation density (490 validation points), overall accuracy was 85% and chance agreement was 39%, resulting in a kappa coefficient of 0.76. This is substantial to almost perfect accuracy based on the Landis and Koch (1977) accuracy thresholds for the kappa coefficient. Based on substrate type, vegetation density, and vegetation type, attribute accuracy for the supervised classification is estimated to be 85-96%. See Zeigler and others (2019) for more details on this validation exercise.
These raster layers (Rock12_SubType.tif, Rock12_VegType.tif, Rock12_VegDen.tif) show discrete landcover classes and were produced based on a supervised classification of aerial imagery (Rock12_SupClas.tif). Geomorphic settings contained in Rock12_GeoSet.tif were hand-digitized with reference to the orthoimagery and DEM (see Source Information). The data were reviewed using standard USGS review procedures. No checks for topological consistency in addition to those described in the Attribute Accuracy Report were performed on these data.
This dataset is clipped to a custom boundary and may not include the entire spatial extent of the original source dataset. However, the custom boundary spans the entire coverage of the study area relevant to the broader research program (see Zeigler and others, 2019 for more details). The data are therefore considered complete for the information presented as described in the abstract section. Users are advised to read the rest of the metadata record carefully for additional details.
Source_Information:
Source_Citation:
Citation_Information:
Originator:
Department of Commerce (DOC), National Oceanic and Atmospheric Administration (NOAA), National Ocean Service (NOS), National Geodetic Survey (NGS), Remote Sensing Division
Publication_Date: 20121115
Title:
Hurricane Sandy: Rapid Response Imagery of the Surrounding Regions
Geospatial_Data_Presentation_Form: Raster Digital Data (Aerial Imagery)
Publication_Information:
Publication_Place: Silver Spring, MD
Publisher: NOAA's Ocean Service, National Geodetic Survey (NGS)
Online_Linkage: https://geodesy.noaa.gov/storm_archive/storms/sandy/index.html
Online_Linkage:
Online_Linkage: http://ngs.woc.noaa.gov/storms/sandy/
Type_of_Source_Media: digital data
Source_Time_Period_of_Content:
Time_Period_Information:
Range_of_Dates/Times:
Beginning_Date: 20121031
Ending_Date: 20121106
Source_Currentness_Reference: ground condition
Source_Citation_Abbreviation: Orthoimage
Source_Contribution:
The orthoimage was used to classify the scene into landcover types. Source data were distributed at 0.35 m pixel resolution, in horizontal datum NAD83. Downloaded on 2/15/2016. Data were projected to UTM Zone 18N (EPSG:26918) using the ‘Project Raster’ tool in ArcToolbox (version 10.4.1).
Source_Information:
Source_Citation:
Type_of_Source_Media: digital data
Source_Time_Period_of_Content:
Time_Period_Information:
Range_of_Dates/Times:
Beginning_Date: 20121116
Ending_Date: 20121116
Source_Currentness_Reference: ground condition
Source_Citation_Abbreviation: DEM
Source_Contribution:
Elevation data used for visual reference while digitizing dunes and geomorphic feature types. Source data were downloaded from
https://coast.noaa.gov/dataviewer in horizontal datum NAD 1983, UTM Zone 18N (EPSG:26918), vertical datum NAVD88. Downloaded on 2/15/2016.
Process_Step:
Process_Description:
Rock12_SupClas.tif
Using the orthoimage, we conducted a supervised classification to delineate substrate and vegetation characteristics. Full methods are provided in the associated Methods OFR (Zeigler and others, 2019). For additional information and example figures, see Zeigler and others (2019). All steps were undertaken in ArcGIS and ArcToolbox version 10.4.1.
Because the signatures of some features closely resembled those of others (e.g., water and the reflection from buildings), we ran the classification in three stages and combined results into a single classification of landcover. We created sub-classifications for (1) bare sand and water, (2) marsh, and (3) vegetation and other features in the barrier island’s interior.
To divide the landscape into these three sub-classifications, we hand-digitized masking polygons in ArcGIS. All digitization was performed with the orthoimage for reference. For the marsh classification, we hand-digitized polygons of marsh platforms. For vegetation and inland features, we hand-digitized a polygon that encompassed vegetation and the barrier’s interior by delineating the boundary of open sand and vegetation. Finally, we hand-digitized polygons to encompass areas of human development. ‘Human development’ could include housing communities, commercial infrastructure, recreational facilities, major roads, and shoreline armoring (e.g., jetties). In sites without any human development, such as Cedar Island, we omitted the development component of the following steps.
We then conducted the classification for marsh. We used the Classification toolbar to hand-digitize training polygons for the following classes: (1) Water; (2) Marsh, vegetation or unknown cover; (3) Marsh, unvegetated sediment; and (4) Marsh, shrub or forest. We ran the interactive supervised classification routine in the Classification toolbar and clipped the resulting classification image along the boundaries of the marsh mask. In the Image Analysis window, we masked out areas that overlapped with development polygons, if present (Rock12_Development.shp in larger work). The resulting image matched the extent of the marsh polygons and contained cells at the original resolution of the orthoimage (0.35 m) with one of the four values represented by the training polygons.
Next, we created training polygons in the same manner that represented different types of vegetation in the inland portion of the island. Training samples represented the following classes: (1) Water; (2) Sand; (3) Herbaceous Vegetation or Shrub, assumed sparse to moderate cover (< 20% cover); (4) Herbaceous Vegetation or Shrub, assumed moderate to dense cover (>20% cover); and (5) Shrub or Forest, assumed moderate to dense cover. We ran the interactive supervised classification routine in the Classification toolbar and clipped the resulting classification image along the boundaries of the vegetation mask. In the Image Analysis window, we masked out areas that overlapped with the marsh or development mask. The resulting image matched the extent of the vegetation mask, minus areas overlapping with marsh and development, and contained cells at the original resolution (0.35 m) of the orthoimage. Cells had one of the five values represented by the training polygons.
We created training polygons in the same manner for the final sub-classification. These samples represented the following classes: (1) Water and (2) Sand. We ran the interactive supervised classification routine in the Classification toolbar. We masked out areas that overlapped with the marsh, development, or vegetation masks in the Image Analysis window. The resulting image covered areas not already covered by the marsh, vegetation, or development and contained cells at the aerial imagery’s original resolution (0.35 m) with a value for either sand or water.
Using the Merge tool in the Image Analysis window in ArcGIS, we combined the rasters for (1) marsh, (2) vegetation, (3) sand/water sub-classifications with the (4) rasterized development layer into a single raster (hereafter, the ‘supervised classification’). We changed the supervised classification’s resolution from the aerial imagery’s original cell size (0.35 m) to a 5 m cell size with the Resample tool in ArcToolbox. We used ‘majority’ as the resampling technique, where the new 5 m cell took on the most common value within the 5x5 m area.
In many instances, cells took on the value of ‘NoData’ around the boundaries of the marsh, vegetation, and development masks when we merged individual sub-classifications. We replaced these NoData values using statistical information from surrounding cells according to the following code in the Raster Calculator tool in ArcToolbox:
CON(isnull([raster]), FOCALMAJORITY([raster], rectangle, 5, 5), [raster]).
This code indicates that, if the value of a cell in the supervised classification is NoData, then replace that value with the value held by the majority of cells in the surrounding 5x5 cell window.
The final supervised classification was a raster layer with a 5 m cell size clipped to the extent of the study area with cells taking on one of nine possible values: (1) water, value = 1; (2) sand, value = 2; (3) herbaceous vegetation or low shrub, sparse to moderate cover, value = 9; (4) herbaceous vegetation or low shrub, moderate to dense cover, value = 3; (5) high shrub or forest, moderate to dense cover, value = 4; (6) marsh, sediment, value = 11; (7) marsh, vegetation or unknown cover, value = 7; (8) marsh, high shrub or forest, value = 14; and (9) development, value = 10. Note, this list contains all possible values in the supervised classification we considered within the larger research program (see Zeigler and others, 2019). Not all values may be represented in this supervised classification specifically.
Process_Date: 2018
Process_Step:
Process_Description:
Geomorphic settings raster (Rock12_GeoSet.tif)
We created the geomorphic settings raster dataset (Rock12_GeoSet.tif) by first delineating the boundaries of individual features (e.g., the edges of washover fans), combining these features into a single shapefile, and converting that shapefile to a raster. Each individual geomorphic feature was given a value for identification purposes (‘value’) as well as a priority number (‘priority’) on which to base the merging of all features into a single layer.
For additional information and example figures, see Zeigler and others (2019). All steps were conducted in ArcGIS and ArcToolbox version 10.4.1.
We began by hand-digitizing the boundaries for marsh (value = 6, priority = 6) based on the visual inspection of the orthoimage. The rear-edges of marsh boundaries were drawn such that they extended out past the study area boundary and included areas of open water that were directly across from marsh platforms.
We created a polygon for beach (value = 1, priority = 5) using the shoreline (rock12_shoreline.shp in the larger work) and custom study area boundary polygons. The custom study area boundary spans the entire coverage of the Rockaway Peninsula site relevant to the broader research program (see Zeigler and others, 2019 for more details). The boundary polygon’s spatial extent does not affect the quality or completeness of the dataset.
To create the beach polygon, we used the ‘Erase’ tool in ArcToolbox to mask out the shoreline polygons and the hand-digitized marsh polygons from a copy of the full study area polygon. This process essentially removed the features of marsh and anything interior of the shoreline from the study area boundary, leaving only the beach geomorphic setting in the final polygon. Together, the beach and marsh polygons covered the outermost edges of the study area.
We also hand-digitized the boundaries of dunes (value = 3, priority = 7) using the DEM (displayed in 1 m intervals), dune crest and dune toe points (rock12_DCpts.shp and rock12_DTpts.shp in larger work and Doran and others 2017), and the orthoimage as guides. We digitized the rear of dune complexes such that the boundary fell landward of dune crest points (if present) and along a clear decrease in elevation (according to visual inspection of the source DEM). The front of the dune was also drawn such that the boundary fell along a clear decrease in elevation and passed through dune toe points. We used the ‘Snapping’ feature in the Editing toolbar of ArcGIS to ensure that the front boundary of dune complexes coincided with the dune toe points. In some instances, particularly in washovers, dune crest points were present without obvious changes in elevation (in the DEM) or the presence of dune toe points. In these cases, we hand-digitized the dune boundary such that it tightly encompassed dune crest points.
Boundaries of the remaining geomorphic features were hand-digitized according to visual inspection of the orthoimage and in reference to boundaries already created for beach, marsh, and dunes. The rear, landward boundary of the backshore (value = 2, priority = 4) was digitized such that it overlapped with dune polygons (but did not extend past the landward-most boundary of the dune). When dunes were not present, the landward boundary of the backshore was drawn where non-vegetated sand ended and dense vegetation began. We extended the backshore polygon boundary seaward beyond the shoreline (or the boundary of marsh in the case of the back-barrier).
The characteristic fan-shaped boundaries of washovers (value = 4, priority = 3) were hand-digitized such that they overlapped with the boundaries of marsh along the back-barrier and with the boundaries of dune complexes on the seaward side of the island. In instances where dunes were not present seaward of a washover, we digitized the washover boundary along wrack if visible in the orthoimage.
We hand-digitized polygons for ridge/swale complexes (value = 7, priority = 2), when present, at the boundary of open sand and vegetation according to the orthoimage. Note: this geomorphic setting was not present at all sites and may not be represented.
For barrier interior (value = 5, priority = 1), the final geomorphic setting, we used the Erase tool in ArcToolbox to remove all other geomorphic settings (i.e., beach, backshore, dune, washover, marsh, and ridge/swale) from a copy of the study area polygon (not published). In this way, the barrier interior setting occurred only in the absence of any other geomorphic features.
Using the Merge tool in ArcToolbox, we combined the individual polygon shapefiles for beach, backshore, dune, washover, barrier interior, marsh, and ridge/swale into a single shapefile. We used the Integrate tool in ArcToolbox (distance = 1 m) to close small gaps between polygons.
We then converted the merged polygon shapefile to a raster with a 5 m cell size, using the priority value for each geomorphic setting type in the attribute table to determine which geomorphic setting type took precedence when two or more settings overlapped (where a raster cell takes the value of the feature with the highest priority value). We selected ‘maximum area cell assignment’ to determine cell values. The extent was set to the study area extent in the geoprocessing Environment Settings window so that the converted raster would be clipped to the extent of the study area.
The final raster provided a categorical map of geomorphic features with every 5x5 m cell categorized as one of the seven possible geomorphic settings.
Process_Date: 2018
Process_Step:
Process_Description:
Reclassification for the substrate type layer (Rock12_SubType.tif):
All steps were undertaken in ArcGIS and ArcToolbox version 10.4.1. For additional information and example figures, see Zeigler and others (2019).
Each of the classes contained in the supervised classification is associated with a substrate type, vegetation type, and vegetation density for later use in creating raster coverages for these characteristics. For substrate type, we reclassified the original supervised classification value to the following substrate types and associated values:
(1) water -> Water, 4444;
(2) sand -> Sand or ShellGravelCobble, 7777;
(3) herbaceous or shrub, sparse to moderate cover -> Sand or ShellGravelCobble, 7777;
(4) herbaceous or shrub, moderate to dense cover -> Sand or ShellGravelCobble, 7777;
(5) shrub or forest -> Unknown, 9999;
(6) marsh, sediment -> Sand or MudPeat, 1000;
(7) marsh, vegetation or unknown cover -> MudPeat, 3333;
(8) marsh, shrub or forest -> Unknown, 9999;
(9) development -> Development, 6666.
Note, this list contains all possible values for substrate types considered within the larger research program (see Zeigler and others, 2019). Not all values may be represented for this site.
We made one manual change to the reclassification of the supervised classification to create the substrate type layer. For raster cells that were classified as ‘beach’ in the geomorphic settings layer (Rock12_GeoSet.tif), we manually reclassified the substrate type for these cells to ‘Water’ with 'None' for both vegetation density and vegetation type. The beach geomorphic setting represented all raster cells outside of the shoreline polygons (rock12_shoreline.shp in larger work), and we assumed that these areas would be underwater at some point depending on the tide. We reclassified everything seaward of the full island shoreline as Water to maintain consistency in landcover classifications relative to tide fluctuations.
Process_Date: 2018
Process_Step:
Process_Description:
Reclassification for the vegetation density layer (Rock12_VegDen.tif):
All steps were undertaken in ArcGIS and ArcToolbox version 10.4.1. For additional information and example figures, see Zeigler and others (2019).
For vegetation density, we reclassified the original supervised classification value to the following vegetation densities and associated values:
(1) water -> None, 111;
(2) sand -> None or Sparse, 666;
(3) herbaceous or shrub, sparse to moderate cover -> Sparse or Moderate, 777;
(4) herbaceous or shrub, moderate to dense cover -> Moderate or Dense, 888;
(5) shrub or forest -> Moderate or Dense, 888;
(6) marsh, sediment -> None or Sparse, 666;
(7) marsh, vegetation or unknown cover -> Unknown, 9999;
(8) marsh, shrub or forest -> Moderate or Dense, 888;
(9) development -> Development, 555.
Note, this list contains all possible values for vegetation density classes considered within the larger research program (see Zeigler and others, 2019). Not all values may be represented in this dataset.
We made one manual change to the reclassification of the supervised classification to create the vegetation density layer. For raster cells that were classified as ‘beach’ in the geomorphic settings layer (Rock12_GeoSet.tif), we manually reclassified the substrate type for these cells as ‘Water’ with 'None' for both vegetation density and vegetation type. The beach geomorphic setting represented all raster cells outside of the shoreline polygons (rock12_shoreline.shp in larger work), and we assumed that these areas would be underwater at some point depending on the tide. This also reduced misclassifications where beach wrack (or dead materials washed up on to the beach) was incorrectly identified as vegetation, which was a common issue in the supervised classification process. Thus, we reclassified everything outside of the shoreline polygons as 'None', or lacking terrestrial vegetation.
Process_Date: 2018
Process_Step:
Process_Description:
Reclassification for the vegetation type layer (Rock12_VegType.tif):
All steps were undertaken in ArcGIS and ArcToolbox version 10.4.1. For additional information and example figures, see Zeigler and others (2019).
For vegetation type, we reclassified the original supervised classification value to the following vegetation types and associated values:
(1) water -> none, 11;
(2) sand -> none or herbaceous, 77;
(3) herbaceous or shrub, sparse to moderate cover -> herbaceous or shrub, 88;
(4) herbaceous or shrub, moderate to dense cover -> herbaceous or shrub, 88;
(5) shrub or forest -> shrub or forest, 99;
(6) marsh, sediment -> none or herbaceous, 77;
(7) marsh, vegetation or unknown cover -> unknown, 9999;
(8) marsh, shrub or forest -> shrub or forest, 99;
(9) development -> development, 55.
Note, this list contains all possible values for vegetation types considered within the larger research program (see Zeigler and others, 2019). Not all values may be represented for this site.
We made one manual change to the reclassification of the supervised classification to create the vegetation type layer. For raster cells that were classified as ‘beach’ in the geomorphic settings layer (Rock12_GeoSet.tif), we manually reclassified the substrate type for these cells as ‘Water’ with 'None' for both vegetation density and vegetation type. The beach geomorphic setting represented all raster cells outside of the shoreline polygons (rock12_shoreline.shp in larger work), and we assumed that these areas would be underwater at some point depending on the tide. We reclassified everything seaward of the full island shoreline as Water to maintain consistency in landcover classifications relative to tide fluctuations. This also reduced misclassifications where beach wrack (or dead materials washed up on to the beach) was incorrectly identified as vegetation, which was a common issue in the supervised classification process. For these reasons, we reclassified everything seaward of the full island shoreline as none, or lacking in terrestrial vegetation.
Process_Date: 2018
Process_Step:
Process_Description:
Added keywords section with USGS persistent identifier as theme keyword.
Process_Date: 20200810
Process_Contact:
Contact_Information:
Contact_Organization_Primary:
Contact_Organization: U.S. Geological Survey
Contact_Person: VeeAnn A. Cross
Contact_Position: Marine Geologist
Contact_Address:
Address_Type: Mailing and Physical
Address: 384 Woods Hole Road
City: Woods Hole
State_or_Province: MA
Postal_Code: 02543-1598
Contact_Voice_Telephone: 508-548-8700 x2251
Contact_Facsimile_Telephone: 508-457-2310
Contact_Electronic_Mail_Address: vatnipp@usgs.gov