2013-14 Massachusetts Lidar-Derived Dune Crest Point Data

2013–14 lidar elevation data were downloaded from NOAA’s online data access viewer (https://coast.noaa.gov/dataviewer/#/) in LAS format (NAD83 2011, UTM Zone 19N (meters) and NAVD88 (meters)) for ten regions along the Massachusetts coast. These regions include: the North Shore, South Shore, Cape Cod Bay, Outer Cape, South Cape, part of the South Coast along Buzzard’s Bay, the southwestern, southern, and eastern coastlines of Martha’s Vineyard, the Vineyard Sound-facing coastline of Martha’s Vineyard, the Atlantic Ocean-facing coastlines of Muskeget, Tuckernuck, and Nantucket Islands, and the Nantucket Sound-facing coastlines of Muskeget, Tuckernuck, and Nantucket Islands. Each of the ten regions has its own coast-following reference line along which the lidar data are processed, and each reference line is divided into segments.
The lidar point data were processed segment by segment along each region’s reference line to extract the peak of the most seaward dune following methods developed by Stockdon and others, 2012. The elevation data for each segment were interpolated in MATLAB 2017a to create a shore-parallel grid with a resolution of 10 m in the longshore direction and 2.5 m in the cross-shore. Each 10 m wide grid row is a cross-shore profile along which dune features are extracted. Each segment grid begins at profile number one, and for each profile, the MATLAB dune crest extraction algorithm finds inflection points in the slope of the data, classifies them as peaks and troughs, and selects the highest peak as the crest of the dune. The vertical error, which is based on the scatter of the data, is also calculated for each point. See Stockdon and others, 2012 for details.
The MATLAB algorithm was run for all segments in each region. Within the program, the geographic coordinates (WGS84) were calculated for each point and the automatically-extracted dune crest data, which includes dune crest elevation, error, and location in both eastings/northings and latitude/longitude, were written to ArcGIS shapefiles using MATLAB’s shapewrite.m script and Keyhole Markup Language (kml) files.
The dune crest locations were manually examined to ensure that an appropriate and consistent alongshore feature was identified. This was done by simultaneously viewing the data in three ways: dune crest kml files were viewed on 2014 Google Earth imagery, the shapefiles were viewed on DigitalGlobe world imagery in ArcMap, and the dune crest and lidar data (interpolated and point data) were viewed for each grid profile in a MATLAB Graphical User Interface (GUI). The GUI allowed for the deleting or moving of dune crest points as necessary to produce a final shapefile of dune peaks for each region. In the absence of a dune, the beach berm or seaward edge of the headland, cliff, bluff or hard structure (e.g., road, parking lot, seawall) was extracted since those features would serve as the first line of defense during a storm. Much manual editing and deleting was required in developed regions as well as in regions with complex topography. The dune crest shapefiles for each of the ten regions were combined in ArcMap (ArcToolbox >> Data Management Tools >> General >> Merge) to create a final shapefile of dune crest position and elevation for the entire state.

A “Class” field was added to the dune crest attribute table to distinguish between actual dune crests and other features such as seawalls or headlands. The dune crest points were viewed over imagery (2014 Google Earth Imagery and DigitalGlobe world imagery in ArcMap), and each point was assigned one of four classes based on the material of the feature and its susceptibility to erosion: 1). The dune crest point is at the peak of a dune, the berm crest of a sandy beach, or the top seaward edge of a bluff. 1A). (1-Armored) The dune crest point is at the peak of an armored dune or the top seaward edge of an armored bluff. Armoring usually consisted of a rocky wall at the base of the feature. Note that sand fencing isn’t considered armoring so dunes with sand fencing were classified as (1), not (1A). 2). The dune crest point is at the crest of a rocky beach or at the top seaward edge of a cliff or headland. 3). The dune crest point is at the top seaward edge of a hard (non-natural) structure such as a road, parking lot, seawall, etc.).
A coordinate system transformation of the dune crest data was required for the points to correspond with the beach and marsh shoreline data that are part of this update to the Shoreline Change Project. Using ArcGIS version 10.5, the dune crest shapefile was transformed from the UTM zone 19 (NAD83 2011) coordinate system to the WGS84 geographic coordinate system (ArcToolbox >> Data Management Tools >> Projections and Transformations >> Project) using the “NAD_1983_HARN_To_NAD_1983_2011 + NAD_1983_HARN_To_WGS_1984” transformation. The DBF file of the shapefile was then saved as a CSV file using Microsoft Excel (version 2016).

Added keywords section with USGS persistent identifier as theme keyword.