Elevation of marsh units in Cape Cod National Seashore salt marsh complex, Massachusetts

This process step and all subsequent process steps were performed by the same person, Zafer Defne, in ArcMap (ver. 10.5.1) using tools from ArcToolbox, unless otherwise stated. For complex operations, names of specific tools used are given in CAPITAL letters (any critical parameters used are given in parentheses, separated by a semicolon, immediately after the tool name). The input and output file names are provided in [square brackets] when necessary. Units for length and area calculations are meters (m) and square meters (m2) unless otherwise stated.
a) Set the data frame coordinate system and projection to NAD 1983 UTM Zone 19N. Prepare elevation dataset to be clipped to the boundaries polygon. MOSAIC TO NEW RASTER(Pixel type= 32 bit float; Mosaic operator=Mean; Cellsize=1 m) all of the CoNED raster to a new raster dataset with 1 m resolution [elev_mosaic.tif].
b) Calculate values for each marsh unit. DISSOLVE(Input features=[CMU_CACOp]; Dissolve field=FID_CMU; Statistics field=ATOT_M2 and FLG with Statistics type= First) to dissolve the vegetated and unvegetated classes in each marsh unit. FLG flags anomalous marsh units based on absence of vegetation, elevation and surface area, and is set to (-1) no vegetated area; (-10) marsh unit elevation higher than the 99.8 percentile; (-100) marsh unit elevation less than 0.2 percentile; (-1000) marsh unit surface area less than 900 m2; and (0) no flag. Combination of negative values indicate combination of flags. For example, -1001 indicates no vegetated area and area less than 900 m2. PROJECT(Input coordinate system=WGS 1984 Web Mercator Auxiliary Sphere; Output coordinate system=NAD 1983 UTM Zone 19N; Geographic transformation=WGS 1984 (ITRF00) to NAD 1983) the feature dataset to obtain dissolved marsh units [mu_diss.shp].
c) Clip the raster with the polygon features by using EXTRACT BY MASK to obtain [elev_idw_mu.tif].
d) ZONAL STATISTICS AS TABLE(Feature Zone Data=[mu_diss.shp]; zone field=FID_CMU; Value Raster=[elev_idw_mu.tif]; Output Table=mu; Statistics type=MEAN)
e) ADD JOIN(Layer name=[mu_diss.shp]; Input join field=FID_CMU; Join Table=mu; Output join field=FID_CMU) to marsh units layer based on the "FID_CMU" field
f) Export dataset as a shapefile [mu_elev.shp] and keep the 'mu_MEAN" and "mu_ATOT_M2" fields from the join operation.
g) Calculate elevation of the vegetated areas only. SELECT BY ATTRIBUTES FROM [CMU_CACOp] where "CLASS"= 'vegetated'. SELECT BY ATTRIBUTES FROM [CMU_CACOp] where "CLASS"= 'vegetated'. Follow steps (c) to (f) using the selected features to export [veg_elev.shp].
h) JOIN FIELD(Input table=[mu_elev.shp]; Input join field=mu_FID_CMU; Join table=[veg_elev.shp]; Output join field=vg_FID_CMU; Join fields=vg_AVEG_M2,vg_ATOT_M2,vg_MEAN). Change the field names to VG_ELEV and MU_ELEV respectively for the fields with the mean elevation for the vegetated area only and mean elevation for the entire marsh unit. Set VG_ELEV to -9999 if there is no vegetated area in the marsh unit. Name AVEG_M2 and ATOT_M2 for vegetated and total area, respectively, and export to a new shapefile. [mu_elev_CACO.shp].
i) SELECT BY ATTRIBUTES FROM [CMU_CACOp] where "CLASS"= 'vegetated', and join field "APGN_M2" from CMU_CACOp.shp on the feature dataset [mu_elev_CACO.shp].
j) Rearrange field names and change the projection for better performance of web services with online base maps. PROJECT(Input coordinate system=NAD 1983 UTM Zone 19N; Output coordinate system=WGS 1984 Web Mercator Auxiliary Sphere; Geographic transformation=WGS 1984 (ITRF00) to NAD 1983) the feature dataset to obtain the final conceptual marsh units [mu_elev_CACOp.shp].

Added keywords section with USGS persistent identifier as theme keyword.