points, transects, beach width: Barrier island geomorphology and shorebird habitat metrics at 50-m alongshore transects and 5-m cross-shore points: Edwin B. Forsythe NWR, NJ, 2010

Frequently anticipated questions:

• What does this data set describe?
  1. How might this data set be cited?
  2. What geographic area does the data set cover?
  3. What does it look like?
  4. Does the data set describe conditions during a particular time period?
  5. What is the general form of this data set?
  6. How does the data set represent geographic features?
  7. How does the data set describe geographic features?
• Who produced the data set?
  1. Who are the originators of the data set?
  2. Who also contributed to the data set?
  3. To whom should users address questions about the data?
• Why was the data set created?
• How was the data set created?
  1. From what previous works were the data drawn?
  2. How were the data generated, processed, and modified?
  3. What similar or related data should the user be aware of?
• How reliable are the data; what problems remain in the data set?
  1. How well have the observations been checked?
  2. How accurate are the geographic locations?
  3. How accurate are the heights or depths?
  4. Where are the gaps in the data? What is missing?
  5. How consistent are the relationships among the data, including topology?
• How can someone get a copy of the data set?
  1. Are there legal restrictions on access or use of the data?
  2. Who distributes the data?
  3. What's the catalog number I need to order this data set?
  4. What legal disclaimers am I supposed to read?
  5. How can I download or order the data?
• Who wrote the metadata?

What does this data set describe?

1. How might this data set be cited?
   Sturdivant, Emily J., Zeigler, Sara L., and Gutierrez, Benjamin T., 2019, points, transects, beach width: Barrier island geomorphology and shorebird habitat metrics at 50-m alongshore transects and 5-m cross-shore points: Edwin B. Forsythe NWR, NJ, 2010: data release DOI:10.5066/P944FPA4, U.S. Geological Survey, Coastal and Marine Geology Program, Woods Hole Coastal and Marine Science Center, Woods Hole, MA.

   Online Links:

   • https://doi.org/10.5066/P944FPA4
   • https://www.sciencebase.gov/catalog/item/5d0bc8fbe4b0941bde4fc5d1

   This is part of the following larger work.
   Sturdivant, Emily J., Zeigler, Sara L., Gutierrez, Benjamin T., and Weber, Kathryn M., 2019, Barrier island geomorphology and shorebird habitat metrics—Four sites in New York, New Jersey, and Virginia, 2010–2014: data release DOI:10.5066/P944FPA4, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://doi.org/10.5066/P944FPA4
   • https://www.sciencebase.gov/catalog/item/5be5c5bce4b0b3fc5cf8c7cb

   Other_Citation_Details:

   Suggested citation: Sturdivant, E.J., Zeigler, S.L., Gutierrez, B.T., and Weber, K.M., 2019, Barrier island geomorphology and shorebird habitat metrics—Four sites in New York, New Jersey, and Virginia, 2010–2014: U.S. Geological Survey data release, https://doi.org/10.5066/P944FPA4.

2. What geographic area does the data set cover?

   West_Bounding_Coordinate: -74.39145948
   East_Bounding_Coordinate: -74.09177767
   North_Bounding_Coordinate: 39.77959626
   South_Bounding_Coordinate: 39.43181496
   

3. What does it look like?

   https://www.sciencebase.gov/catalog/file/get/5d0bc8fbe4b0941bde4fc5d1/?name=ebf_pts_trans_ubw_browse.png (PNG)

   Example views of the 5-m points, which follow the cross-shore transects. The points were projected from the 'seg_x' and 'seg_y' and overlaid on the DEM. Two panels include the shoreline delineation (larger work). The points are colored by a different variable in each panel, as labeled.

4. Does the data set describe conditions during a particular time period?

   Calendar_Date: 2010

   Currentness_Reference:

   Ground condition measured by source data.

5. What is the general form of this data set?

   Geospatial_Data_Presentation_Form: tabular
   

6. How does the data set represent geographic features?
   1. How are geographic features stored in the data set?
      
   2. What coordinate system is used to represent geographic features?
      

      Grid_Coordinate_System_Name: Universal Transverse Mercator
      Universal_Transverse_Mercator:

      UTM_Zone_Number: 18
      Transverse_Mercator:

      Scale_Factor_at_Central_Meridian: 0.9996
      Longitude_of_Central_Meridian: -75.0
      Latitude_of_Projection_Origin: 0.0
      False_Easting: 500000.0
      False_Northing: 0.0
      

      Planar coordinates are encoded using coordinate pair
      Abscissae (x-coordinates) are specified to the nearest 0.6096
      Ordinates (y-coordinates) are specified to the nearest 0.6096
      Planar coordinates are specified in Meter
      The horizontal datum used is D_North_American_1983.
      The ellipsoid used is GRS_1980.
      The semi-major axis of the ellipsoid used is 6378137.0.
      The flattening of the ellipsoid used is 1/298.257222101.
      

      Vertical_Coordinate_System_Definition:

      Altitude_System_Definition:

      Altitude_Datum_Name: North American Vertical Datum of 1988 (NAVD88)
      Altitude_Resolution: 0.001
      Altitude_Distance_Units: meters
      Altitude_Encoding_Method: Attribute values
      

7. How does the data set describe geographic features?

   ebf_trans.shp

   Geometry of 878 shore-normal transects used as the base for sampling barrier island geomorphology. Various ID fields are included to permit later joining. NoData value for all: -99999. (Source: Producer defined)

   FID

   Internal feature number (Source: Esri) Sequential unique whole numbers that are automatically generated.

   Shape

   Feature geometry. (Source: Esri)

   Value	Definition
   Polyline	Polyline geometry

   sort_ID

   Unique identifier that orders transects sequentially along the shoreline. (Source: Producer defined)

   Range of values
   Minimum:	1
   Maximum:	878
   Units:	integer

   Shape_Length

   Length of the feature automatically calculated by ArcGIS. (Source: Producer defined)

   Range of values
   Minimum:	3829.4702575968727
   Maximum:	5086.298187009409
   Units:	meters

   TransectId

   Transect identifier corresponding to that used in NASC transects. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	1.0
   Maximum:	1646.0
   Units:	integer

   TransOrder

   Transect identifier corresponding to that used in NASC transects. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	1.0
   Maximum:	1676.0
   Units:	integer

   DD_ID

   Identifier that orders transects sequentially along the shoreline and is unique across all sites analyzed for Zeigler and others, 2018. (Source: Producer defined)

   Range of values
   Minimum:	30001
   Maximum:	30878
   Units:	integer

   ebf10_pts.csv

   Geomorphology metrics at 108974 discrete points along shore-normal transects. Transect-averaged values are also included. (Source: Producer defined)

   SplitSort

   Unique identifier that sorts the points by transect and distance from MHW. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0
   Maximum:	108973
   Units:	integer

   seg_x

   Point easting referenced to UTM Zone 18N, NAD 83. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557023.6634
   Maximum:	577632.6146999998
   Units:	meters

   seg_y

   Point northing referenced to UTM Zone 18N, NAD 83. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367190.025699999
   Maximum:	4402020.0723
   Units:	meters

   Dist_Seg

   Along-transect distance to the mean high water (MHW) shoreline position. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.0
   Maximum:	1965.5886022206423
   Units:	meters

   Dist_MHWbay

   Distance of the point from the bayside shoreline. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-1379.0903703261793
   Maximum:	1548.9698317746531
   Units:	meters

   DistSegDH

   Distance to the dune crest position along the same transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-312.81199001582485
   Maximum:	1887.0488643202004
   Units:	meters

   DistSegDL

   Distance to the dune toe position along the same transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-297.71378573696916
   Maximum:	1943.9976223051478
   Units:	meters

   DistSegArm

   Distance to the armoring position along the same transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-128.43318620452968
   Maximum:	1286.5545226098743
   Units:	meters

   ptZ

   Elevation referenced to NAVD88. Derived from the corresponding pixel in the 5-m DEM. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-0.25
   Maximum:	9.09999942779541
   Units:	meters

   ptSlp

   Mean slope of the corresponding 5-m pixel in the slope surface. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.0
   Maximum:	42.99135208129883
   Units:	percent

   ptZmhw

   Elevation referenced to MHW. Derived from attribute ptZ (ptZmhw = ptZ - MHW offset). NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-0.6799999999999999
   Maximum:	8.66999942779541
   Units:	meters

   GeoSet

   The geomorphic setting indicates the classification of the major topographic features within approximately 5 m. NoData value of -99999. (Source: Producer defined)

   Value	Definition
   1	Beach. The relatively thick and temporary accumulation of loose, water-borne material (usually well-sorted sand and pebbles, accompanied by mud, cobbles, boulders, and smoothed rock and shell fragments) that is in active transit along, or deposited on, the shore zone between the limits of low water and high water (Neuendorf and others, 2011). In this study, the beach geomorphic setting occurred between the ocean- or sound-side study area boundary and the mean high water line (approximated for by the upper-most wrack line in the field).
   2	Backshore. The upper, usually dry, zone of the shore or beach, lying between the high-water line of mean spring tides and the upper limit of shore-zone processes; it is acted upon by waves or covered by water only during exceptionally severe storms or unusually high tides (Neuendorf and others, 2011). In this study, the Backshore geomorphic setting occurred between the mean high water line and either (i) the dune toe, (ii) the edge of developed areas, or (iii) the edge of dense vegetation (or forest).
   3	Dunes. A low mound, ridge, bank, or hill of loose, windblown granular material (generally sand), either bare or covered by vegetation, capable of movement from place to place but retaining its characteristic shape (Neuendorf and others, 2011). In this study, “dune” also describes low-lying areas between dunes (or “interdune” regions) that are part of the larger dune complex.
   4	Washover. A fan of material deposited from the ocean landward on a mainland beach or barrier island, produced by storm waves breaking over low parts of the mainland beach or barrier and depositing sediment either landward (mainland beaches) or across a barrier island into the bay/sound (barrier islands). A washover typically displays a characteristic fan-like shape (Neuendorf and others, 2011).
   5	Barrier Interior. In this study, the Barrier Interior geomorphic setting described all areas spanning the interior boundary of the dunes (or backshore in the absence of dunes) on the ocean-side to the interior boundary of the marsh, dunes, or backshore on the back-barrier side. This setting was typically used to describe areas that did not fall into any other geomorphic setting (e.g., washovers, ridge/swale complexes).
   6	Marsh. A relatively flat, low-lying, intermittently water-covered area with generally halophytic grasses existing landward of a barrier island (Neuendorf and others, 2011).
   7	Ridge/Swale complex: Long subparallel ridges and swales aligned obliquely across the regional trend of the contours (Neuendorf et al. 2011).
   -99999	Missing value.

   SubType

   Substrate type: dominant classification of substrate within approximately 5 m. NoData value of -99999. Curly brackets ('{}') around two values are used in the modeling process. (Source: Zeigler and others, 2018)

   Value	Definition
   1000	Sand or Mud/Peat: In this study, wet, sandy substrates could not be differentiated from mud/peat in marshy areas. Therefore, we identified substrate as being either of these two types in the Substrate raster layer. Mud/Peat is a sticky, fine-grained, predominantly clay- or silt-sized marine detrital sediment (Neuendorf and others 2011). Sand included rock or mineral grains with diameters between 0.074 and 4.76 mm (Neuendorf et al. 2011).
   1111	Sand: predominantly sandy substrates that contain rock or mineral grains with diameters between 0.074 and 4.76 mm (Neuendorf et al. 2011) with no discernible shells fragments or large rock fragments.
   2222	Shell/Gravel/Cobble: substrates containing a mixture of sand, shell or rock fragments, or large rocks.
   3333	MudPeat: A sticky, fine-grained, predominantly clay- or silt-sized marine detrital sediment (Neuendorf and others 2011).
   4444	Water: Any location that (i) is always submerged (e.g., locations several meters into the ocean, bay, or inland water body), (ii) was submerged at the time aerial imagery was captured (i.e., intertidal regions of beaches), or (iii) was not submerged at the time aerial imagery was captured but was seaward of the shoreline polygon.
   6666	Development: Any location that fell within areas obviously influenced by anthropogenic activities (e.g., housing developments, paved roads or parking lots, recreational sports fields, etc.).
   {1111, 2222}	Sand or ShellGravelCobble: Substrate is either predominantly sand or a mix of sand with shell or rock fragments. Sand included rock or mineral grains with diameters between 0.074 and 4.76 mm (Neuendorf and others 2011) with no discernible shells fragments or large rock fragments. ShellGravelCobble described substrates containing a mixture of sand, shell or rock fragments, or large rocks.
   -99999	Unknown: Substrate type could not be determined based on aerial imagery

   VegDens

   The dominant type of vegetation within approximately 5 m. NoData value of -99999. Curly brackets ('{}') around two values are used in the modeling process. (Source: Zeigler and others, 2018)

   Value	Definition
   111	None: Areas lacking terrestrial vegetation of any type. Such areas were most frequently associated with the beach geomorphic setting (found seaward of the study area shoreline) assumed to be covered by water.
   222	Sparse: areas where vegetation was apparent and covered less than 20% of the 5x5-m raster cell.
   333	Moderate: areas where vegetation appeared to cover 20–90% of the 5x5-m raster cell.
   444	Dense: areas where vegetation appeared to cover greater than 90% of the 5x5-m raster cell.
   555	Development: In this study, we selected development as the vegetation density for any location that fell within areas obviously influenced by anthropogenic activities (e.g., housing developments, paved roads or parking lots, recreational sports fields, etc.).
   {111, 222}	None or Sparse: The given location either appeared to lack vegetation completely in the aerial imagery or, if vegetation was apparent, covered less than 20% of the approximately 5 x 5 m area designated by the transect point.
   {222, 333}	Sparse or Moderate: In this study, areas with sparse vegetation could not consistently be differentiated from areas containing moderate vegetation in the orthoimagery alone. Therefore, we identified vegetation as being either of these two types in the Vegetation Density raster layer. Areas described as having 'Sparse' vegetation appeared to have vegetation that covered less than 20% of the 5x5-m raster cell. In areas described as 'Moderate', vegetation appeared to cover 20-90% of the 5x5-m raster cell.
   {333, 444}	Moderate or Dense: The given location either appeared to have vegetation covering 20-90% or >90% of the approximately 5 x 5 m area designated by the transect point.
   -99999	Unknown: Vegetation density could not be determined based on the aerial imagery

   VegType

   The dominant type of vegetation within approximately 5 m. NoData value of -99999. Curly brackets ('{}') around two values are used in the modeling process. (Source: Zeigler and others, 2018)

   Value	Definition
   11	None: Areas lacking terrestrial vegetation of any type. Such areas were associated with the beach geomorphic setting (found seaward of the study area shoreline) assumed to be covered by water.
   22	Herbaceous: areas containing primarily herbaceous vegetation of the forb/herb growth habit (USDA 2015) and lacking shrubs, trees, or any other vegetation with woody stems (Neuendorf et al. 2011). In this study, the Herbaceous vegetation type typically described the vegetation cover found in Godfrey’s (1976) (i) ‘grassland’ ecological zone along the backshore and dunes, dominated by beach grasses (e.g., Ammophila breviligulata) or (ii) ‘intertidal marsh’ ecological zone dominated by cordgrass (e.g., Spartina patens).
   33	Shrub: Areas containing low (less than 5 m height), multi-stemmed woody plants of the subshrub or shrub growth habits (USDA 2015). In this study, the Shrub vegetation type typically described vegetation cover found in Godfrey’s (1976) heath-like ‘shrublands’ ecological zone in stable dune systems.
   44	Forest: Areas containing trees and tall (> 5 m) shrubs of the tree growth habit (USDA 2015b). In this study, the Forest vegetation type typically described vegetation cover found in Godfrey’s (1976) ‘woodlands–forests’ ecological zone found in barrier island interiors and dominated by deciduous (e.g., Quercus velutina), pine (e.g., Pinus rigida), and juniper (e.g., Juniperus virginiana) species.
   55	Development: In this study, we selected Development as the vegetation type for any location that fell within areas obviously influenced by anthropogenic activities (e.g., housing developments, paved roads or parking lots, recreational sports fields, etc.).
   {11, 22}	None or Herbaceous: The given location either lacked vegetation or contained primarily herbaceous vegetation of the forb/herb growth habit (USDA 2015). In this study, the herbaceous vegetation type typically described the vegetation cover found in Godfrey’s (1976) (i) ‘grassland’ ecological zone along the backshore and dunes, dominated by beach grasses (e.g., Ammophila breviligulata) or (ii) ‘intertidal marsh’ ecological zone dominated by cordgrass (e.g., Spartina patens).
   {22, 33}	Herbaceous or Shrub: The given location either contained primarily herbaceous vegetation of the forb/herb growth habit or low (less than 5 m height), multi-stemmed woody plants of the subshrub or shrub growth habits (USDA 2015). In this study, the herbaceous vegetation type typically described the vegetation cover found in Godfrey’s (1976) (i) ‘grassland’ ecological zone along the backshore and dunes, dominated by beach grasses (e.g., Ammophila breviligulata) or (ii) ‘intertidal marsh’ ecological zone dominated by cordgrass (e.g., Spartina patens). The shrub vegetation type typically described vegetation cover found in Godfrey’s (1976) heath-like ‘shrublands’ ecological zone in stable dune systems.
   {33, 44}	Shrub or Forest: The given location contained either tall (greater than 5 m height) multi-stemmed woody plants of the shrub growth habit or plants of the tree growth habit (USDA 2015). In this study, the shrub vegetation type typically described vegetation cover found in Godfrey’s (1976) heath-like ‘shrublands’ ecological zone in stable dune systems. The forest vegetation type typically described vegetation cover found in Godfrey’s (1976) ‘woodlands–forests’ ecological zone found in barrier island interiors and dominated by deciduous (e.g., Quercus velutina), pine (e.g., Pinus rigida), and juniper (e.g., Juniperus virginiana) species.
   -99999	Unknown: Vegetation type could not be determined based on aerial imagery

   sort_ID

   Unique identifier that orders transects sequentially along the shoreline. (Source: Producer defined)

   Range of values
   Minimum:	1.0
   Maximum:	878.0
   Units:	integer

   TransOrder

   Transect identifier corresponding to that used in NASC transects. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	1.0
   Maximum:	1676.0
   Units:	integer

   TransectId

   Transect identifier corresponding to that used in NASC transects. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	1.0
   Maximum:	1646.0
   Units:	integer

   DD_ID

   Identifier that orders transects sequentially along the shoreline and is unique across all sites analyzed for Zeigler and others, 2018. (Source: Producer defined)

   Range of values
   Minimum:	30001
   Maximum:	30878
   Units:	integer

   Azimuth

   Bearing of the NASC transect measured in degrees clockwise from North. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	260.58
   Maximum:	299.08
   Units:	degrees

   LRR

   Rate of shoreline change determined by linear regression. These are defined for each NASC transect. NoData value of -99999. (Source: Himmelstoss and others (2010))

   Range of values
   Minimum:	-8.61
   Maximum:	15.39
   Units:	meters/year

   SL_x

   Easting of position of MHW shoreline along transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557644.6747000003
   Maximum:	577574.6350999996
   Units:	meters

   SL_y

   Northing of position of MHW shoreline along transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367209.305199999
   Maximum:	4401290.4297
   Units:	meters

   Bslope

   Slope of foreshore at the shoreline point nearest to the transects within 25 m. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-0.151534
   Maximum:	-0.005069
   Units:	percent

   DL_x

   Easting of the nearest dune toe point within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557703.7824868169
   Maximum:	576860.6697092913
   Units:	meters

   DL_y

   Northing of the nearest dune toe point within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367258.169780839
   Maximum:	4400792.011568513
   Units:	meters

   DL_z

   Elevation of the nearest dune toe point within 25 meters of the transect, referenced to NAVD88. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.873698
   Maximum:	4.354614
   Units:	meters

   DL_zmhw

   Elevation of the nearest dune toe point within 25 meters of the transect, referenced to MHW. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.443698
   Maximum:	3.9246139999999996
   Units:	meters

   DL_snapX

   Easting of the nearest dune toe point within 25 meters of the transect, snapped to the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557694.8090529181
   Maximum:	576862.2902176064
   Units:	meters

   DL_snapY

   Northing of the nearest dune toe point within 25 meters of the transect, snapped to the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367243.500404998
   Maximum:	4400795.479067003
   Units:	meters

   DH_x

   Easting of the nearest dune crest point within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557594.1588149052
   Maximum:	577409.572759321
   Units:	meters

   DH_y

   Northing of the nearest dune crest point within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367240.233419549
   Maximum:	4401729.394236792
   Units:	meters

   DH_z

   Elevation of the nearest dune crest point within 25 meters of the transect, referenced to NAVD88. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	1.374724
   Maximum:	8.343055
   Units:	meters

   DH_zmhw

   Elevation of the nearest dune crest point within 25 meters of the transect, referenced to MHW. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.9447240000000001
   Maximum:	7.913055
   Units:	meters

   DH_snapX

   Easting of the nearest dune crest point within 25 meters of the transect, snapped to the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	557594.7290145863
   Maximum:	577411.4442311807
   Units:	meters

   DH_snapY

   Northing of the nearest dune crest point within 25 meters of the transect, snapped to the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4367241.128335746
   Maximum:	4401723.23637284
   Units:	meters

   Arm_x

   Easting of the nearest armoring structure within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	563303.7812571861
   Maximum:	575928.2168329637
   Units:	meters

   Arm_y

   Northing of the nearest armoring structure within 25 meters of the transect. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4375960.548755558
   Maximum:	4399355.507504804
   Units:	meters

   Arm_z

   Elevation of the nearest armoring structure within 25 meters of the transect, referenced to NAVD88. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	1.7999999523162842
   Maximum:	8.09999942779541
   Units:	meters

   Arm_zmhw

   Elevation of the nearest armoring structure within 25 meters of the transect, referenced to MHW. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	1.369999885559082
   Maximum:	7.669999599456787
   Units:	meters

   DistDH

   Horizontal distance along transect from the MHW shoreline to the dune crest (DH_snapX, DH_snapY). NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	17.727273090147442
   Maximum:	313.2066299614794
   Units:	meters

   DistDL

   Horizontal distance along transect from the MHW shoreline to the dune toe (DL_snapX, DL_snapY). NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4.213341862535684
   Maximum:	297.15870333767253
   Units:	meters

   DistArm

   Horizontal distance along transect from the MHW shoreline to the first armoring structure (Arm_snapX, Arm_snapY). NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	22.061101512599688
   Maximum:	128.55759914724618
   Units:	meters

   Dist2Inlet

   Alongshore distance to nearest inlet. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	11.689137858131241
   Maximum:	17272.980876352663
   Units:	meters

   WidthPart

   Horizontal distance between seaward shoreline and the first break in the land, such as intervening water bodies like back-barrier creeks. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	1.1192753067473384
   Maximum:	1549.0290738207536
   Units:	meters

   WidthLand

   Barrier Island Width (m). Width of the barrier island cross-section. Measurement only considers width of land and not intervening water bodies such as back barrier creeks. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	39.58940072207548
   Maximum:	1686.243290968214
   Units:	meters

   WidthFull

   Barrier Island Full Width (m). Horizontal distance along transect from seaward shoreline to farthest bayside shoreline. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	39.58940072207548
   Maximum:	1968.4164759325477
   Units:	meters

   uBW

   Upper beach width. The horizontal distance from MHW shoreline to the dune toe or equivalent. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	4.213341862535684
   Maximum:	297.15870333767253
   Units:	meters

   uBH

   Upper beach height. The vertical distance from MHW to the dune toe or equivalent. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	0.443698
   Maximum:	6.860000133514404
   Units:	meters

   ub_feat

   Feature (dune toe, dune crest, or armoring) that indicates the top of beach, used to calculate beach width and beach height. (Source: Producer defined)

   Value	Definition
   -99999	No feature available to calculate beach width.
   DL	Dune toe.
   DH	Dune crest.
   Arm	Shorefront armoring used in place of dune toe to calculated beach width.

   mean_Zmhw

   Mean elevation along transect, referenced to MHW tidal datum. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-0.24111111616094907
   Maximum:	2.520806391815985
   Units:	meters

   max_Zmhw

   Maximum elevation along transect, referenced to MHW tidal datum. NoData value of -99999. (Source: Producer defined)

   Range of values
   Minimum:	-0.13000001788139343
   Maximum:	8.66999942779541
   Units:	meters

   Construction

   Construction: denotes presence of construction-based human erosion management activities. Unknown is indicated with a NoData value of -99999. (Source: Producer defined)

   Value	Definition
   111	None-denotes no constructed features.
   222	Denotes presence of soft approaches (sand fencing, geotubes).
   333	Denotes presence of hard approaches (rip rap, seawall, jetty, groyne).
   444	Denotes presence of hard and soft approaches.

   Development

   Development intensity: denotes intensity of human development. Unknown is indicated with a NoData value of -99999. (Source: Producer defined)

   Value	Definition
   111	Denotes no human development.
   222	Denotes light development.
   333	Denotes moderate development.

   Nourishment

   Nourishment activity: denotes presence of nourishment-based human erosion management activities. Unknown is indicated with a NoData value of -99999. (Source: Producer defined)

   Value	Definition
   111	Denotes no nourishment.
   222	Denotes occasional nourishment.

   ebf10_ubw.tif values attribute table

   Values attribute table (ebf10_ubw.tif.vat.dbf), which indicates the beach width of every cell in the raster (ebf10_ubw.tif). The raster is composed of 5,031 x 7,764 grid cells. Raster cells outside the bounds of the data have a NoData value of 65535. (Source: Producer defined)

   Value

   Matches the sort_ID identifier values used to uniquely identify transects in the supplemented NASC transects and the 5-m points (Source: Producer defined)

   Range of values
   Minimum:	1
   Maximum:	878
   Units:	integer

   Count

   Number of pixels pertaining to the transect identified by the Value field. (Source: Producer defined)

   Range of values
   Minimum:	1205
   Maximum:	16126
   Units:	count

   uBW

   Upper beach width. The horizontal distance from the MHW shoreline to the dune toe or equivalent. The fill value for unavailable data is -99999. (Source: Producer defined)

   Range of values
   Minimum:	4.213341863
   Maximum:	297.1587033
   Units:	meters

   Entity_and_Attribute_Overview:

   This section provides a separate detailed entity and attribute information section for each dataset described in these metadata. The complete dataset is included in the 5-m points. The other files are included for convenience. Each point represents a 5 x 5 m sample area with the geomorphological characteristics specified. All calculations for length are in meter units and were based on the UTM zone 18N NAD83 projection. NoData is indicated with the value -99999. Please review the individual attribute descriptions for detailed information.

   Entity_and_Attribute_Detail_Citation: Methods Open-File Report by Zeigler and others, 2019
   

Who produced the data set?

1. Who are the originators of the data set? (may include formal authors, digital compilers, and editors)
   • Emily J. Sturdivant
   • Sara L. Zeigler
   • Benjamin T. Gutierrez
2. Who also contributed to the data set?
3. To whom should users address questions about the data?
   Emily J. Sturdivant

   U.S. Geological Survey

   Geographer

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   (508) 548-8700 x2230 (voice)

   (508) 457-2310 (FAX)

   esturdivant@usgs.gov

Why was the data set created?

These data provide samples of barrier island characteristics for use in Bayesian networks to model geomorphology and shorebird nesting habitat condition and change. Transects spaced every 50 m alongshore from the National Assessment of Shoreline Change (NASC; ebf_trans.shp) were extended and supplemented to use as the base sampling unit to compile and sample barrier island characteristics. Attributes are provided at 5-m sample points (ebf10_pts.csv) along each transect and include both values extracted for the entire transect and values extracted at each point. The attributes in the points file can be aggregated by transect and joined to the transect features by the field sort_ID. The beach width values calculated for each transect are also provided as a spatially continuous raster (ebf10_ubw.tif). See Zeigler and others (2019) for additional details.

How was the data set created?

1. From what previous works were the data drawn?

   NASC transects (source 1 of 3)

   Himmelstoss, Emily A., Kratzmann, Meredith G., Hapke, Cheryl J., Thieler, E. Robert, and List, Jeffrey, 2010, National Assessment of Shoreline Change: A GIS Compilation of Vector Shorelines and Associated Shoreline Change Data for the New England and Mid-Atlantic Coasts: Open-File Report 2010-1119, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://pubs.usgs.gov/of/2010/1119/
   • https://pubs.usgs.gov/of/2010/1119/data_catalog.html

   Type_of_Source_Media: digital data
   Source_Contribution:

   Shore-normal transects with long term shoreline change rates from the National Assessment of Shoreline Change (NASC) (NewJerseyN_LT.shp). The data are distributed as an Esri polyline shapefile referenced to World Geodetic System 1984 (WGS84). They were downloaded in 2017.

   Orthoimage (source 2 of 3)

   USDA/FSA - Aerial Photography Field Office, 2010, USDA-FSA-APFO NAIP MrSID Mosaic (NAIP10).

   Online Links:

   • https://www.fsa.usda.gov/programs-and-services/aerial-photography/imagery-programs/naip-imagery/ordering-naip-imagery/index
   • https://www.fsa.usda.gov/programs-and-services/aerial-photography/imagery-products/imagery-catalogs/index
   • mailto:apfo.sales@slc.usda.gov

   Other_Citation_Details:

   Imagery is available by request. Request NAIP10 1m imagery for Long Beach Island, Ocean, NJ and Pullen Island, Atlantic, NJ from apfo.sales@slc.usda.gov. Further instructions at the first link.

   Type_of_Source_Media: digital data
   Source_Contribution:

   Visual imagery used for digitizing shorefront development, tidal inlets, and for QA/QC. Source data were distributed at 1 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).

   DEM (source 3 of 3)

   Department of Commerce (DOC), National Oceanic and Atmospheric Administration (NOAA), National Ocean Service (NOS), Office for Coastal Management (OCM), and JALBTCX (Joint Airborne Lidar Bathymetry Technical Center of eXpertise), Unpublished material, 2010 US Army Corps of Engineers (USACE) Joint Airborne Lidar Bathymetry Technical Center of eXpertise (JALBTCX) Lidar: New Jersey (Topo): NOAA's Ocean Service, Office for Coastal Management (OCM), Charleston, SC.

   Online Links:

   • https://coast.noaa.gov/dataviewer/#/lidar/search/where:ID=1133
   • https://coast.noaa.gov/htdata/lidar1_z/geoid12a/data/1133
   • https://coast.noaa.gov/dataviewer
   • https://coast.noaa.gov
   • https://inport.nmfs.noaa.gov/inport/item/50090

   Type_of_Source_Media: digital data
   Source_Contribution:

   Elevation. 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.

2. How were the data generated, processed, and modified?

   Date: 2018 (process 1 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others 2019). The iPython notebook used for processing is distributed with this dataset (extractor_ebf10.ipynb).
   ebf_trans.shp
   NASC transects from Himmelstoss and others (2010) were modified for the purposes of this study. Transects were extended inland to cover the width of the island, and additional transects were added to fill alongshore gaps greater than 50 m.
   Transects were extended by using the last two vertices of each transect to programmatically place the end of the line 3,000 m beyond the end of the original line segment. Python (version 3) with the modules ArcPy and Collections in an ArcGIS Pro 2.0 environment were used for programming (see function core/functions_warcpy.ExtendLine in bi-transect-extractor (Sturdivant, 2019)). The process transformed the data from WGS84 to NAD83 UTM Zone 18N using WGS_1984_(ITRF00)_To_NAD_1983_2011 (WKID: 108354, accuracy: 0.1 m). It used the coordinates of the last two vertices of each feature and calculated a new second coordinate at a specified distance (3,000 m) along the prolongation of the line from the first coordinate.
   Next, transects were manually added to fill alongshore gaps greater than 50 m. To do so we copied the extended NASC transects, replaced their attribute values with fill values (-99999), moved groups of them to fill gaps, deleted any that were unchanged, and then merged the altered file back to the NASC transects. ID values were assigned that ordered the transects consecutively along the shoreline using the ArcGIS Spatial Sort tool in spatial increments.
   We eliminated transect overlap in certain locations by manually clipping the transects to the first intersection point with an overlapping transect. While doing so, we prioritized the original NASC transect geometries. Overlapping transects retained the azimuths of the original lines but in some cases were shortened.

   Date: 2018 (process 2 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others, 2019) and technical steps can be found in the source code (Sturdivant, 2019). The iPython notebook used to create the point dataset (extractor_ebf10.ipynb) is distributed in this data release. All steps described below were performed in the iPython notebook extractor_ebf10.ipynb, which used Python 3 and ArcPy distributed with ArcGIS Pro 2.0. Manual steps indicated in the notebook were executed in a session of ArcGIS Pro 2.0.
   ebf10_pts.csv, part 1
   First, we calculate values that apply to entire cross-island transects. Supplemented NASC transects are populated with the shoreline change rate values from the original NASC transects. Lines that are not present in the NASC transects are populated with fill values (-99999) for LRR, TRANSECTID, and TRANSORDER fields.
   MHW position and foreshore slope along transect:
   Each transect is assigned the foreshore slope (Bslope) from the nearest shoreline point within 25 m (see larger work). The MHW shoreline easting and northing locate the intersection of the transect with the oceanside shoreline. These values are populated for each transect as follows (using Python 3 distributed with ArcGIS Pro 2.0, especially the modules ArcPy, numpy, and pandas): (1) Create a line representing the oceanside shoreline by converting the shoreline polygons to lines, clipping them at the inlet lines, and selecting only those segments that intersect shoreline points; (2) get 'SL_x' and 'SL_y' at the intersection point between the transect and the oceanside shoreline (created from ebf10_shoreline.shp in larger work); (3) find the closest shoreline point to that intersection point; and (4) assign the slope value from the shoreline point to the transect ('Bslope').
   Dune positions along transects:
   'DL_x', 'DL_y', and 'DL_zMHW' are the easting, northing, and height above MHW, respectively, of the nearest dune toe point (ebf10_DLpts.shp in larger work) within 25 meters of the transect. Likewise, 'DH_x', 'DH_y', and 'DH_zMHW' are the easting, northing, and height above MHW, respectively, of the nearest dune crest point within 25 meters (ebf10_DCpts.shp in larger work). We find the XYZ positions of the nearest dune crest and dune toe within 25 meters for each transect using bi-transect-extractor (Sturdivant, 2019) to execute the following, repeated for both dune crest and dune toe: (1) evaluate the distance from the transect to every dune crest/toe point and find the shortest of these distances (distanceTo geometry method in ArcPy); (2) if the distance is less than 25 m, return the elevation from the point and the XY position of the point ‘snapped’ to the transect (snapToLine geometry method); (3) if there are no points within 25 m of the transect, populate the transect with null/fill values (pandas); and (4) convert the elevations to the MHW datum by applying the MHW offset (0.43 m based on Weber and others, 2005).
   'Arm_x', 'Arm_y', and 'Arm_zMHW' are the easting, northing, and height above MHW, respectively, where an artificial structure crosses the transect in the vicinity of the beach. These features are meant to supplement the dune toe dataset by providing an upper limit to the beach in areas where dune toe extraction was confounded by the presence of an artificial structure. Values are populated for each transect as follows: (1) prioritizing areas where dune toe features were not successfully extracted, use orthoimagery, supplemented with the DEM, to manually digitize line segments on the oceanside face of artificial impediments to sediment ("armoring"), such as sand-fencing, sandbags, seawalls, etc.; (2) get the positions of intersection between the digitized armoring lines and the transects (Intersect tool from the Overlay toolset); (3) extract the elevation value at each intersection point from the DEM (Extract Multi Values to Points tool from Spatial Analyst); and (4) convert the elevations to the MHW datum by applying the MHW offset (0.43 m based on Weber and others, 2005).
   Beach width and height:
   Upper beach width ('uBW') and upper beach height ('uBH') are calculated based on the difference in position between two points: the position of MHW along the transect ('SL_x', 'SL_y') and the dune toe position or equivalent (usually 'DL_x', 'DL_y'). All calculations are performed by the function calc_BeachWidth_fill in bi-transect-extractor (Sturdivant, 2019). In some cases, the dune toe is not appropriate to designate the "top of beach" so beach width and height are calculated from either the position of the dune toe, the dune crest, or the base of an armoring structure. The dune crest is only considered a possibility if the dune crest elevation (DH_zMHW) is less than or equal to 2.5 m. They are calculated as follows, relying primarily on numpy and the snapToLine geometry method in ArcPy, and using pandas for data storage/organization:
   (1) Find the position along the transect of an orthogonal line drawn to the dune point ('DL_x', 'DL_y' and 'DH_x', 'DH_y'). (2) Calculate distances from MHW to the position along the transect of the dune toe ('DistDL'), dune crest ('DistDH'), and armoring ('DistArm'). (3) Conditionally select the appropriate feature to represent "top of beach." Dune toe is prioritized. If it is not available and 'DH_zMHW' is less than or equal to maxDH, use dune crest. If neither of the dune feature positions satisfy the conditions and an armoring feature intersects with the transect, use the armoring position. If none of the three are possible, 'uBW' and 'uBH' are given a NoData value of -99999. (4) Copy the distance to shoreline and height above MHW ('Dist--', '---zMHW') to 'uBW' and 'uBH', respectively.
   Distance to inlet:
   Distance to nearest tidal inlet ('Dist2Inlet') is computed as alongshore distance of each sampling transect from the nearest tidal inlet. Inlets are manually delineated during the creation of the shoreline polygon file. This distance includes changes in the path of the shoreline rather than just a straight-line distance between each transect and the inlet and reflects sediment transport pathways. It is measured using the shoreline polygon(s) and the delineated tidal inlets (see shoreline data in larger work) as follows: (1) split the shoreline polygon(s) at the tidal inlets by converting the shoreline polygon(s) into a polyline feature class with the inlet lines included (Feature to Line in Data Management); (2) retain only the oceanside segments of the shoreline by deleting all segments that do not intersect any shoreline points (disjoint geometry method in ArcPy data access module); (3) if the shoreline is bounded on both sides by an inlet, measure the distance to both and assign the minimum distance of the two or if the shoreline meets only one inlet (meaning the study area ends before the island ends), use the distance to the only inlet (cut, disjoint, and length geometry methods and properties in ArcPy)
   Island widths:
   Barrier Island Width ('WidthLand') is calculated as the above-water distance between the back-barrier and seaward MHW shorelines. 'WidthLand' only included regions of the barrier within the shoreline polygon(s) (ebf10_shoreline.shp in larger work). We also measure the shore-to-shore extent of the island, which includes space occupied by waterways ('WidthFull') and the width of only the most seaward portion of land within the shoreline ('WidthPart'). These are calculated as follows: (1) clip the transect to the shoreline polygon(s) (Clip in the Analysis toolbox); (2) for 'WidthLand', get the length of the multipart line segment from "SHAPE@LENGTH" feature class attribute, which will include only the remaining portions of the transect; (3) for 'WidthPart', convert the clipped transect from multipart to singlepart and get the length of the first line segment, which should be the most seaward; (4) for 'WidthFull', calculate the distance between the first vertex and the last vertex of the clipped transect (Feature Class to NumPy Array with explode to points, pandas groupby, numpy hypot). Person who carried out this activity:
   

   Emily J. Sturdivant

   U.S. Geological Survey

   Geographer

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   508-548-8700 x2230 (voice)

   508-457-2310 (FAX)

   esturdivant@usgs.gov

   Date: 2018 (process 3 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others, 2019) and technical steps can be found in the source code (Sturdivant, 2019). The iPython notebook used to create the point dataset (extractor_ebf10.ipynb) is distributed in this data release. All steps described below were performed in the iPython notebook extractor_ebf10.ipynb, which used Python 3 and ArcPy distributed with ArcGIS Pro 2.0. Manual steps indicated in the notebook were executed in a session of ArcGIS Pro 2.0.
   ebf10_pts.csv, part 2
   Nourishment, Construction, Development:
   We manually assigned coded values for the attributes Nourishment, Construction, and Development by comparing the transect positions to ancillary datasets in ArcGIS 10.5. These datasets included the inventory of habitat modification (Rice 2015), available aerial imagery, and the development layer included in this data release. See the associated Methods OFR (Zeigler and others, 2019) for details and example figures. Person who carried out this activity:
   

   Benjamin T. Gutierrez

   U.S. Geological Survey

   Geologist

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   508-548-8700 x2289 (voice)

   508-457-2310 (FAX)

   bgutierrez@usgs.gov

   Date: 2018 (process 4 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others, 2019). All steps described below were performed in the iPython notebook extractor_ebf10.ipynb, which used Python 3 and ArcPy distributed with ArcGIS Pro 2.0. Manual steps indicated in the notebook were executed in a session of ArcGIS Pro 2.0.
   ebf10_pts.csv, part 3
   The point dataset represents 5 m sampling of the land along each shore-normal transect (ebf_trans.shp; see Zeigler and others, 2019). The 5-m point file is created from the supplemented NASC transects as follows: (1) Manually shorten overlapping transects to the first point of intersection. When a transect overlaps with a supplementary transect, shorten the supplementary transect to the intersection point. (2) Clip the transects to the shoreline polygon of the barrier island, retaining only those portions of the transects that represent land. (3) Create a point along each transect every 5 m starting from the ocean-side shoreline. Person who carried out this activity:
   

   Emily J. Sturdivant

   U.S. Geological Survey

   Geographer

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   508-548-8700 x2230 (voice)

   508-457-2310 (FAX)

   esturdivant@usgs.gov

   Date: 2018 (process 5 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others, 2019). All steps described below were performed in the iPython notebook extractor_ebf10.ipynb, which used Python 3 and ArcPy distributed with ArcGIS Pro 2.0 (Sturdivant, 2019). Manual steps indicated in the notebook were executed in a session of ArcGIS Pro 2.0.
   ebf10_pts.csv, part 4
   Point identifier:
   We populate the 5-m points with a numerical identifier ('SplitSort') of the 5-m data sampling points at a particular study site, sorted by order along shoreline and by distance from oceanside. 'SplitSort' values are populated by sorting the points by 'sort_ID' and 'Dist_Seg' (see below).
   Distances:
   'Dist_Seg' and 'Dist_MHWbay' measure the distance of the point from the seaward shoreline and bayside shoreline respectively. 'Dist_Seg' is calculated as the Euclidean distance between the point and the shoreline ('SL_x', 'SL_y'). 'Dist_MHWbay' is calculated by subtracting the 'Dist_Seg' value from the transect 'WidthPart'. This is performed by the function prep_points in bi-transect-extractor (Sturdivant, 2019).
   'DistSegDH', 'DistSegDL', and 'DistSegArm' measure the distance of each 5-m point from the dune crest, dune toe, and dune armoring position respectively along a particular transect. They are calculated as the Euclidean distance between the 5-m point and the given feature.
   Elevation and slope:
   'ptZ' and 'ptSlp' are the elevation (NAVD88) and slope at the 5-m cell corresponding to the point. We use the 5-m DEM to generate a slope surface (Slope tool in 3D Analyst). The elevation and slope values are assigned to the points using the Extract Multi Values to Points tool in Spatial Analyst, which is called by bi-transect-extractor (Sturdivant, 2019). 'ptZmhw' is calculated from 'ptZ' by subtracting the MHW offset (0.43 m based on Weber and others, 2005).
   Transect-averaged elevation:
   We calculate the per-transect mean and maximum barrier elevation (mean_Zmhw, max_Zmhw) from the 5-m elevations (ptZmhw). Mean barrier elevations are calculated for only those transects having less than 20 percent missing values within the 5-m points. Locations not satisfying this criterion are assigned a fill value. This is performed by the function aggregate_z in bi-transect-extractor (Sturdivant, 2019).
   Habitat variables:
   Variables for vegetation type (VegType), vegetation density (VegDens), substrate type (SubType), distance to foraging habitat (DisMOSH), and geomorphic setting (GeoSet) are populated from the corresponding raster layers provided with the larger work (files EF10_VegType.tif, EF10_VegDen.tif, EF10_SubType.tif, EF10_DisMOSH.tif, EF10_GeoSet.tif). They are assigned to the points using the Extract Multi Values to Points tool in Spatial Analyst, which is called by bi-transect-extractor (Sturdivant, 2019).
   The values are recoded as follows using bi-transect-extractor, which is documented in the file extractor_ebf10.ipynb.
   SubType: 7777:'{1111, 2222}', 1000:'{1111, 3333}’ VegType: 77:'{11, 22}', 88:'{22, 33}', 99:'{33, 44}' VegDens: 666: '{111, 222}', 777: '{222, 333}', 888: '{333, 444}', 999: '{222, 333, 444}'
   To calculate the values, the data are passed between pandas dataframe format and ArcGIS geodatabase feature class. Once all values have been calculated, the pandas dataframe is saved in comma-separated values (CSV) format. Curly brackets ('{}') around two values are used in the modeling process. Person who carried out this activity:
   

   Emily J. Sturdivant

   U.S. Geological Survey

   Geographer

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   508-548-8700 x2230 (voice)

   508-457-2310 (FAX)

   esturdivant@usgs.gov

   Date: 2018 (process 6 of 7)

   Full methods are provided in the associated Methods OFR (Zeigler and others, 2019). For the detailed documentation of the processing, see the files extractor_ebf10.ipynb and the source code bi-transect-extractor (Sturdivant, 2019).
   ebf10_ubw.tif
   The width of the upper beach (from mean high water (MHW) shoreline to either dune crest, dune toe, or coastal armoring/development) is calculated along the supplemented NASC transects.
   Calculate beach width along transect:
   To calculate the beach width for the transects, upper beach width ('uBW') is calculated based on the difference in position between two points: the position of MHW along the transect ('SL_x', 'SL_y') and the dune toe position or equivalent (usually 'DL_x', 'DL_y'). It is calculated using bi-transect-extractor (Sturdivant, 2019) and documented in extractor_ebf10.ipynb as follows:
   (1) Find the position along the transect of an orthogonal line drawn to the dune point ('DL_x', 'DL_y' and 'DH_x', 'DH_y'). (2) Calculate distances from MHW to the position along the transect of the dune toe ('DistDL'), dune crest ('DistDH'), and armoring ('DistArm'). (3) Conditionally select the appropriate feature to represent "top of beach." Dune toe is prioritized. If it is not available and 'DH_zMHW' is less than or equal to maxDH, use dune crest. If neither of the dune feature positions satisfy the conditions and an armoring feature intersects with the transect, use the armoring position. If none of the three are possible, 'uBW' is given a NoData value of -99999. (4) Copy the selected distance to shoreline and height above MHW ('Dist--') to 'uBW'.
   Assign values to raster
   A transect ID raster is created from the transect file. Transects are modified from the supplemented transect file (ebf_trans.shp) by manually shortening overlapping transects to the first point of intersection. When a transect overlaps with a supplementary transect, the supplementary transect is shortened to the intersection point. The beach width values from each transect are assigned to the cells that represent that transect using the JoinField tool in the data management toolbox, which is also an automated process performed by bi-transect-extractor.
   The raster file is exported to GeoTiff using the Export Raster tool. Raster cells outside the bounds of the data are assigned a NoData value of 65535. Person who carried out this activity:
   

   Emily J. Sturdivant

   U.S. Geological Survey

   Geographer

   384 Woods Hole Road

   Woods Hole, MA
   

   USA

   508-548-8700 x2230 (voice)

   508-457-2310 (FAX)

   esturdivant@usgs.gov

   Date: 10-Aug-2020 (process 7 of 7)

   Added keywords section with USGS persistent identifier as theme keyword. Person who carried out this activity:
   

   U.S. Geological Survey

   Attn: VeeAnn A. Cross

   Marine Geologist

   384 Woods Hole Road

   Woods Hole, MA
   

   508-548-8700 x2251 (voice)

   508-457-2310 (FAX)

   vatnipp@usgs.gov

3. What similar or related data should the user be aware of?
   Zeigler, Sara L., Sturdivant, Emily J., and Gutierrez, Benjamin T., 2019, Evaluating barrier island characteristics and piping plover (Charadrius melodus) habitat availability along the U.S. Atlantic coast—Geospatial approaches and methodology: Open-File Report 2019–1071, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://doi.org/10.3133/ofr20191071

   Other_Citation_Details:

   Details the methods used to process these data for use in barrier island and piping plover habitat modeling.

   Himmelstoss, Emily A., Kratzmann, Meredith, Hapke, Cheryl, Thieler, E. Robert, and List, Jeffrey, 2010, The National Assessment of Shoreline Change: A GIS Compilation of Vector Shorelines and Associated Shoreline Change Data for the New England and Mid-Atlantic Coasts: Open-File Report 2010-1119, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://pubs.usgs.gov/of/2010/1119/

   Weber, Kathryn M., List, Jeffrey H., and Karen L. M. Morgan, 2005, An Operational Mean High Water Datum for Determination of Shoreline Position from Topographic Lidar Data: Open-File Report 2005-1027, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://pubs.usgs.gov/of/2005/1027/

   Rice, Tracy M., 2015, Inventory of Habitat Modifications to Sandy Beaches in the U.S. Atlantic Coast Breeding Range of the Piping Plover (Charadrius melodus) prior to Hurricane Sandy: Maine to the North Shore and Peconic Estuary of New York: U.S. Fish and Wildlife Service, Hadley, Massachusetts.

   Online Links:

   • http://northatlanticlcc.org/teams/coastal-resiliency/projects/hurricane-sandy/increasing-resiliency-of-beach-habitats-and-species-1/InventoryofHabitatModificationstoSandyBeachesMENY_Rice2015.pdf
   • https://databasin.org/galleries/164daee0855c4228bb6fe8552e704558
   • http://northatlanticlcc.org/products/synthesis-of-tidal-inlet-and-beach-habitat-inventories
   • http://northatlanticlcc.org/teams/coastal-resiliency/projects/beach-and-tidal-inlet-habitat-inventories/beach-and-tidal-inlet-habitat-inventories

   Sturdivant, Emily J., 2019, bi-transect-extractor: software release DOI:10.5066/P915UYMY, U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://doi.org/10.5066/P915UYMY
   • https://code.usgs.gov/cmgp/bi-transect-extractor

   Other_Citation_Details:

   Custom-built package used for data processing. The iPython notebook file distributed with these data was part of the processing to create these data. The notebook relies on this software package.

   Doran, Kara J., Long, Joseph W., Birchler, Justin J, Brenner, Owen T., Hardy, Matthew W., Karen L. M. Morgan, Stockdon, Hilary F., and Torres, Miguel L., 2017, Lidar-derived Beach Morphology (Dune Crest, Dune Toe, and Shoreline) for U.S. Sandy Coastlines: U.S. Geological Survey, Reston, VA.

   Online Links:

   • https://doi.org/10.5066/F7GF0S0Z

   Other_Citation_Details:

   Affiliated datasets created using the same methods. A version of this geomorph points dataset will be available as 10CNT11_morphology.zip. These data were made available to us prior to publication. As a result some of the processing is different from what is described here.

   Neuendorf, Klaus K.E., James P. Mehl, Jr., and Jackson, Julia A., 2011, Glossary of geology. Fifth edition: The American Geosciences Institute, Alexandria, VA.

   Godfrey, Paul J., 1976, Comparative ecology of East Coast barrier islands—Hydrology, soil, vegetation, in Barrier islands and beaches: Technical proceedings of the 1976 Barrier Islands Workshop, Annapolis, Maryland, May 17–18, 1976 p. 5–31, The Conservation Foundation, Washington, D.C..

   U.S. Department of Agriculture (USDA), and USDA Natural Resources Conservation Service (NRCS), 2015, The PLANTS Database (http://plants.usda.gov, 13 January 2014): National Plant Data Team, Greensboro, NC.

How reliable are the data; what problems remain in the data set?

1. How well have the observations been checked?
   Attribute values at each point (ebf10_pts.csv) represent a 5 x 5 m square centered at the point. The values are compiled from multiple sources. Transect-averaged values, such as beach width, may represent 25 m on either side of the indicated transect. The following methods were used to validate attribute accuracy: symbolized display of transect and point attribute values overlaid on input datasets including elevation, island extent, geomorphic feature (ArcGIS Pro version 2.0); spot-checking of values at individual points and comparison with input datasets and topology (ArcGIS Pro version 2.0); and for random transects, plotting the cross-shore profile using elevation values from each point along the transect, distance values and the cross-shore position of the features they were calculated from (Matlab R2015b). These checks were performed by at least two operators.
   Metrics related to the geomorphic feature points (i.e. dune toe, dune crest, shoreline) do not necessarily provide accurate measurements for bay-side areas. For example, in rare situations a transect might intersect multiple ocean-facing shorelines, possibly due to the orientation of an inlet (e.g. transect 1024 at Fire Island, fiis_trans.shp in larger work), but the single dune crest point associated with the transect is simply the closest point, whether it be on the seaward portion of land or not. Thus, the distance to dune crest values may not be realistic for a portion of the points along the transect.
   ebf10_ubw.tif: Accuracy of the beach width value is dependent on the accuracy of the lineage data. Beach widths are calculated to the tenths of meters in precision. However, the values are only calculated every 50 m alongshore. Thus, a given cell value may represent the width of the beach at a transect as far as 25 m removed from the cell. Refer to the process steps for details. Beach width values were spot-checked while verifying the reliability of the transect values. Spot checks were performed by plotting the beach width on a cross-shore elevation profile with the available dune positions positioned on the plot as well. These values were always found to be in agreement unless reported elsewhere.
   More detail is provided elsewhere below (Process Steps and Entity and Attribute Information), in the associated Open-File Report (Zeigler and others, 2019), and in the source code (Sturdivant, 2019).
2. How accurate are the geographic locations?
   Depending on the attribute, the accuracy is assumed to be between 5 and 25 m. The horizontal positional accuracy is dependent on the accuracy of the source data and error incorporated during processing. Refer to the process steps for details.
3. How accurate are the heights or depths?
   The vertical accuracy of those attributes that incorporate vertical position is dependent on the digital elevation model and its source data as well as the MHW datum produced by Weber and others (2005).
4. Where are the gaps in the data? What is missing?
   ebf_trans.shp: Transects (ebf_trans.shp and base feature for ebf10_pts.csv and ebf10_ubw.tif) are spaced alongshore every 50 m within the study area. They cover the entire island width. Some transects may not intersect land in certain years. Curves in the generalized shoreline may cause transects to be spaced more or less densely than the standard 50 m alongshore, especially along the inland side of the study area. Transects may overlap. In such cases, the full extended transect length was preserved allowing overlap so that each transect could represent the full width of the island. Original NASC transects include ID values that match them to the original transect, but the supplemental transects have NoData for NASC fields.
   ebf10_pts.csv: Points were confined to the area within the shoreline polygon included in the larger work and located only along the transects. These exclude areas below mean high water (MHW) elevation on the seaward side and areas below mean tidal level (MTL) on the inland side. Those MHW and MTL contours were generalized to approximately 25 m so in some cases the points may represent submerged areas.
   The points sample every 5 m along shore-normal transects. Curves in the generalized shoreline may cause transects and the resulting points to be spaced more or less densely than the standard 50 m alongshore, especially along the inland side of the study site.
   Before creating the 5-m points, transect lines were manually edited to prevent overlapping transects. Where two transects overlapped, one was shortened to the first point of intersection with the other. When a National Assessment of Shoreline Change (NASC) transect overlapped with a supplementary transect, the supplementary transect was shortened to the intersection point.
   Where values could not be calculated due to lack of input data in the lineage dataset, a NoData value of -99999 was recorded for the attribute. Where a transect does not intersect land, it corresponds to a single point with NoData values for all fields. More detail is provided in the associated Open-File Report (Zeigler and others, 2019) and in the Entity and Attribute Information.
   ebf10_ubw.tif: Before converting the transect vectors to rasters, the transect lines were manually edited to prevent overlapping transects. Where two transects overlapped, one was shortened to the first point of intersection with the other. When a National Assessment of Shoreline Change (NASC) transect overlapped with a supplementary transect, the supplementary transect was shortened to the intersection point. Where values could not be calculated due to lack of input data in the lineage dataset, a fill value of -99999 was recorded for the attribute. More detail is provided in the associated Open-File Report (Zeigler and others, 2019) and in the Entity and Attribute Information.
5. How consistent are the relationships among the observations, including topology?
   These datasets consist of data compiled from multiple sources and aggregated spatially. 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. The primary vertical datum used is NAVD88, consistent with the source datasets lidar and geomorph points. Some fields (those with a 'mhw' suffix and explained in the Entity and Attribute section) include the elevation adjusted to the mean high water (MHW) datum calculated by Weber and others (2005) for the area.
   ebf10_pts.csv: With increasing distance of a point from the seaward face of the island, there is less likely to be a direct orthogonal relationship between the value and the shoreline. Prior to segmenting transects into 5-m points, some transects were shortened to eliminate overlap (see Completeness Report and first Process Step). Thus transect-based values at points near the bayside of the island may pertain to a more distant shoreline.
   ebf_trans.shp: This shapefile consists of line data compiled from NASC transects and manually added. No checks for topological consistency were performed on these data. Transect extent is based on the study area, rather than the shoreline, so that transect positions encompass the land in years with different shoreline positions. As a result, some transects may not intersect land in certain years.
   ebf10_ubw.tif: With increasing distance of a cell from the seaward face of the island, there is less likely to be a direct orthogonal relationship between the cell value and the beach because transects were shortened to eliminate overlap prior to assigning the line values to the raster grid. Cells may not have the same value as neighboring cells directly seaward because of the same overlap-elimination process applied to the transects.

How can someone get a copy of the data set?

Are there legal restrictions on access or use of the data?

Access_Constraints: none
Use_Constraints:

Public domain data from the U.S. Government are freely redistributable with proper metadata and source attribution. Please recognize the U.S. Geological Survey (USGS) as the source of this information.

1. Who distributes the data set? (Distributor 1 of 1)
   
   U.S. Geological Survey - ScienceBase

   Denver Federal Center, Building 810, Mail Stop 302

   Denver, CO
   

   USA

   1-888-275-8747 (voice)

   sciencebase@usgs.gov
2. What's the catalog number I need to order this data set? This dataset contains three individual datasets: point data as a generic comma-separated ASCII formatted data file (ebf10_pts.csv); an Esri shapefile with transect geometries and index values (ebf_trans.shp and other shapefile components); and a 16-bit GeoTIFF file containing the upper beach width values (ebf10_ubw.tif and associated value attribute table). The CSDGM FGDC metadata describing the datasets in XML format (ebf10_pts_trans_ubw_meta.xml) and the browse graphic (ebf_pts_trans_ubw_browse.png) are also included as well as the iPython notebook (extractor_ebf10.ipynb) used for processing. These datasets can be downloaded individually or packaged on-demand in a zip file (see the Digital Transfer Option section).
3. What legal disclaimers am I supposed to read?

   Neither the U.S. Government, the Department of the Interior, nor the USGS, nor any of their employees, contractors, or subcontractors, make any warranty, express or implied, nor assume any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, nor represent that its use would not infringe on privately owned rights. The act of distribution shall not constitute any such warranty, and no responsibility is assumed by the USGS in the use of these data or related materials. Although these data and associated metadata have been reviewed for accuracy and completeness and approved for release by the U.S. Geological Survey (USGS), and have been processed successfully on a computer system at the USGS, no warranty expressed or implied is made regarding the display or utility of the data for other purposes, nor on all computer systems, nor shall the act of distribution constitute any such warranty. The USGS or the U.S. Government shall not be held liable for improper or incorrect use of the data described and/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.

4. How can I download or order the data?
   • Availability in digital form:

     Data format:

     This dataset contains 5-m point data as a comma-separated values (CSV) table (75 MB), a polyline shapefile (0.158 MB), and a 16-bit GeoTIFF of beach width (80 MB). The GeoTIFF and its values attribute table must be stored together. Also included are the iPython notebook used for processing, the CSDGM metadata, and the browse graphic. in format CSV (version ArcGIS Pro 2.0) Esri polyline shapefile, comma-separated values table, and GeoTIFF Size: 155.158

     Network links:

     https://www.sciencebase.gov/catalog/item/5d0bc8fbe4b0941bde4fc5d1 https://www.sciencebase.gov/catalog/file/get/5d0bc8fbe4b0941bde4fc5d1 https://doi.org/10.5066/P944FPA4

   • Cost to order the data: None

5. What hardware or software do I need in order to use the data set?

   To utilize these data, the user must have software capable of reading a comma-delimited data file, shapefile format, and 16-bit GeoTIFF with associated values attribute table.

Who wrote the metadata?

Dates:

Last modified: 19-Nov-2021

Metadata author:

Emily J. Sturdivant

U.S. Geological Survey

Geographer

384 Woods Hole Road

Woods Hole, MA

USA

(508) 548-8700 x2230 (voice)

(508) 457-2310 (FAX)

whsc_data_contact@usgs.gov

Contact_Instructions:

The metadata contact email address is a generic address in the event the metadata contact is no longer with the USGS or the email is otherwise invalid.

Metadata standard:

FGDC Content Standards for Digital Geospatial Metadata (FGDC-STD-001-1998)