Salt Marsh Monitoring in Jamaica Bay , New York from 2003 to 2013 : A Decade of Change from Restoration to Hurricane Sandy

This study used Quickbird-2 and Worldview-2, high resolution satellite imagery, in a multi-temporal salt marsh mapping and change analysis of Jamaica Bay, New York. An object-based image analysis methodology was employed. The study seeks to understand both natural and anthropogenic changes caused by Hurricane Sandy and salt marsh restoration, respectively. The objectives of this study were to: (1) document salt marsh change in Jamaica Bay from 2003 to 2013; (2) determine the impact of Hurricane Sandy on salt marshes within Jamaica Bay; (3) evaluate this long term monitoring methodology; and (4) evaluate the use of multiple sensor derived classifications to conduct change analysis. The study determined changes from 2003 to 2008, 2008 to 2012 and 2012 to 2013 to better understand the impact of restoration and natural disturbances. The study found that 21 ha of salt marsh vegetation was lost from 2003 to 2013. From 2012 to 2013, restoration efforts resulted in an increase of 10.6 ha of salt marsh. Hurricane Sandy breached West Pond, a freshwater environment, causing 3.1 ha of freshwater wetland loss. The natural salt marsh showed a decreasing trend in loss. Larger salt marshes in 2012 tended to add vegetation in 2012–2013 (F4,6 = 13.93, p = 0.0357 and R2 = 0.90). The study provides important information for the resource management of Jamaica Bay.


Introduction
Jamaica Bay, an estuary within the New York City (NYC) limits, is heavily influenced by urbanization.The salt marshes serve as an interface between the Bay and surrounding urban areas.Currently, over a dozen marsh islands span the Bay.Their landscapes are composed of mudflats, a variety of salt marsh plant species, sediment deposited to rebuild drowning salt marsh, transitional vegetation denoting the shift to upland, and human created upland areas.Salt marshes provide numerous ecological benefits such as high biodiversity, improved water quality, flood reduction, and carbon sequestration [1].The wetland ecosystems of New York State, including salt marshes, were reduced by 60% from 1780 to 1980 [2]. Nationally, salt marshes have been under particular stress with increasing rates of loss from 2004 to 2009 caused in part by coastal storms [3].In the past, these trends were exacerbated in the urban-impacted Jamaica Bay.Jamaica Bay's salt marsh loss is severe.Since 1951, approximately 60% of the Bay's salt marsh has converted into mudflats due to a combination of factors including a reduction in sediment supply, changes in tidal regime, nutrient enrichment and increased hydrogen sulfide concentrations [4].This estimate does not include areas of wetlands around the estuary lost to land filling and urbanization.From 1989 to 2003, Jamaica Bay's salt marshes were in rapid decline losing 13.4 ha/year [5].The nitrogen load of the Bay is one factor that may contribute to this high rate of loss [6].
Remote sensing is uniquely suited for monitoring coastal environments, due to the difficulty of in situ access and the high temporal resolution required to understand these dynamic landscapes [7].Remote sensing monitoring of the salt marsh landscape can be used to determine vegetation trends for the entire bay and individual islands, facilitating an assessment of restoration impacts.Remote sensing is an important tool for furthering our understanding of how Jamaica Bay's salt marshes are affected by anthropogenic and natural factors [8,9].This study used imagery data spanning a decade and two high resolution sensor systems.
In October 2012, Hurricane Sandy impacted the coast of New York and surrounding states with high winds and storm surge.It was a 500-year storm surge event at the Manhattan Battery [10].The boroughs of Brooklyn and Queens directly surrounding Jamaica Bay were inundated; the storm caused two million New Yorkers to lose power [11].This study seeks to understand the impact of Hurricane Sandy on salt marsh vegetation within Jamaica Bay.The salt marsh vegetation types of interest are smooth cordgrass (S. alterniflora), high marsh (a mixture of Distichlis spicata, Spartina patens, and Juncus gerardii) and the common reed (Phragmites australis).Successful management of Jamaica Bay is contingent on continuing to further our understanding of the change experienced by the Bay's salt marshes due to both natural disturbance and human impacts.
The objectives of this study were to: (1) document salt marsh changes that occurred in Jamaica Bay from 2003 to 2013; (2) determine the impact of Hurricane Sandy on salt marshes within Jamaica Bay; (3) evaluate this long-term monitoring methodology for the determination of change; and (4) evaluate the use of multiple sensor derived classifications to conduct change analysis.The combination of climate change, sea level rise and their impacts on natural disturbances are expected to have detrimental effects on coastal salt marshes [12]; thereby, enhancing the need for accurate remote sensing monitoring and assessment of coastal wetlands to inform decision-makers.

Study Area
Jamaica Bay is an urban estuary residing within the New York City boroughs of Brooklyn and Queens.Kings County, synonymous with Brooklyn, is the most populated county in New York State [13].Approximately 3704 ha of the Bay are managed by the National Park Service as Jamaica Bay National Wildlife Refuge, a subunit of Gateway National Recreation Area (Figure 1).The region has a humid continental climate with a mean temperature of approximately 10 • C.Over the last 150 years, anthropogenic impacts to Jamaica Bay have been extensive.The Bay's volume has increased 350% while surface area fell by approximately 4856 ha [14].In 2005, Waste Water Treatment Plants serving 1,610,990 people discharged into Jamaica Bay [15].Beginning in 2003, salt marsh islands including Big Egg, Yellow Bar, Rulers Bar, Black Bar, and Elders Point East and West (Figure 1) have undergone salt marsh restoration.After restoration, sites were monitored in situ for 5 years [4].These marsh restoration projects involved the deposition of dredge sediment from channels in the Bay onto the marsh surface then the transplanting and seeding of salt marsh vegetation [16].

Remote Sensing Data
High spatial resolution Quickbird-2 and Worldview-2 data were employed for salt marsh mapping and change analysis.The spatial resolutions of Worldview-2's multispectral and panchromatic sensors are 1.85 m and 0.42 m, respectively; Quickbird-2's resolutions are 2.6 m and 0.62 m, respectively.The Worldview-2 sensor collects eight spectral bands including the Coastal

Remote Sensing Data
High spatial resolution Quickbird-2 and Worldview-2 data were employed for salt marsh mapping and change analysis.The spatial resolutions of Worldview-2's multispectral and panchromatic sensors are 1.85 m and 0.42 m, respectively; Quickbird-2's resolutions are 2.6 m and 0.62 m, respectively.The Worldview-2 sensor collects eight spectral bands including the Coastal Blue, Blue, Green, Yellow, Red, Red Edge, Near Infrared 1 (NIR1), and Near Infrared 2 (NIR2).The Coastal Blue, Yellow, Red Edge, and NIR2 spectral bands of Worldview-2 have been shown to increase the accuracy of wetland vegetation classification [17].This study used Quickbird-2 imagery data acquired on 10 September 2003 and 9 September 2008, and Worldview-2 data acquired on 15 September 2012 and 19 September 2013.The imagery data were geo-rectified to the 2013 imagery.The data were also atmospherically corrected to top of atmosphere reflectance.
This study uses object-based image analysis (OBIA) which first divides an image into objects, using a segmentation algorithm, and then classifies those objects based on their spectral and spatial attributes [18].Object-based change detection (OBCD) utilizes image objects to conduct a change analysis between multiple time periods.The change analysis can be conducted with object attributes, classified objects, multi-temporal image objects, or a hybrid of these techniques [19].This study compared the classified 2003, 2008, 2012 and 2013 objects to understand restoration and Hurricane Sandy's impact on wetlands within Jamaica Bay.

Segmentation
An important component of OBIA classifications is the determination of segmentation scale, which determines the size and similarity of resulting image objects, and parameterization, i.e., the inclusion of texture [20].Texture is the use of a moving window to quantify measures that represent ideas such as coarseness and roughness [19].This study arrived at an appropriate segmentation scale with the comparison of multiple segmentation scales for each time period to maximize intra-segment homogeneity and intersegment heterogeneity [21,22].The parameterization of the resulting image objects included spectral values, texture, geospatial attributes, upland data, vegetation indices, and neighborhood and scene difference attributes (described in Section 2.4).Segmentation scale is the key to accurately mapping a landscape.Scale parameters can be arrived at through "trial-and-error".However, this method risks determining an inappropriate segmentation scale.Over or under segmenting an image can result in lower classification accuracy [23].In addition, segmentation scale can impact the land cover classes that can be accurately mapped [20].This study used the mean shift clustering approach to determine segmentation.Mean shift is a non-parametric segmentation algorithm which groups pixels based on their spectral mean in a feature space.The algorithm has improved accuracy when compared to other clustering techniques [24,25].Mean shift considers a spectral radius in the feature space as the scale parameter, which results in a hierarchical relationship between segmentation scales [26].These factors make the algorithm suitable for multiscale segmentation.
There are different methods for assessing the quality of segmentation.This study assessed segmentation scales with an index of intra-segment homogeneity and intersegment heterogeneity [21].Intersegment heterogeneity was assessed through computation of Global Moran's I that were normalized and then combined with the intra-segment homogeneity, as determined by normalized area controlled variance, to create a single parameter measuring segmentation quality [22].The mean shift segmentation parameters that were determined were minimum size and spectral radius.Minimum size refers to the fewest number of pixels that can compose a segment, and spectral radius is the distance in the feature which a pixel must be within to merge into the segment.Each image date was tested with the parameters from 5 to 50 spectral radii in increments of 1 and minimum size from 5 to 50 in increments of 5. Appropriate segmentation scale for the Worldview-2 2012 data was determined to be a spectral radius of 15 and a minimum size of 5.The 25% most over segmented objects were segmented again at a quantitatively determined appropriate scale of spectral radius 20 and minimum size 5.The same was done for 25% most under segmented objects, for which the appropriate scale was spectral radius 6 and minimum size 5.The appropriate scale for the Worldview-2 2013 data was determined to be spectral radius 22 and minimum size 20.The 25% most over segmented objects were segmented again at a quantitatively determined appropriate scale of spectral radius 27 and minimum size 5.The 25% most under segmented objects were re-segmented at a scale of spectral radius 7 and minimum size 5.The Quickbird-2 data were segmented at a spectral radius of 8 and a minimum size of 20 no additional levels of segmentation were done as this scale adequately captured the landscapes and spectral complexity of the Quickbird-2 data.
The classification was conducted with the Random Forest classifier.Random Forest is a non-parametric ensemble learning algorithm that has been demonstrated to achieve appropriate classification accuracy in a variety of landscapes [27][28][29].The 9 classes used in this study included 6 from a previous study of the Bay [8].These classes included water, mudflat, sand, high marsh, patchy S. alterniflora, and S. alterniflora (>50% vegetation cover).The two S. alterniflora classes were based on percent cover with patchy being between 10% and 49% vegetation cover and S. alterniflora (>50% vegetation cover) being above 50%.Salicornia species are present within the Bay as a small component of the salt marsh [30], and were not prevalent enough to classify on their own.Additional classes included in this study are wrack, upland vegetation, Phragmites, and shadow, however shadow was removed with a decision tree post-classification.The 2003 classification did not include wrack due to the limited separability of the class in those images.These additional classes were included to expand our understanding of the Bay and inform management decisions.

Object Attributes
Spectral attributes included the mean and standard deviation of all available spectral bands.The spatial variables computed were perimeter, area, and nodes.The panchromatic band was utilized to create Grey-Level Co-Occurrence Matrix (GLCM) textural measurements, including inverse difference moment, entropy, contrast, correlation, and uniformity.GLCM and other texture measures have been shown to improve classification accuracies in both Very High Resolution image classification [31] and object-based wetland classification [32].Red Edge-based vegetation indices, have been shown to more accurately discern differences between high density vegetation species [33].In this study, Worldview Vegetation Index (WVVI), Worldview Water Index (WVWI), Red Edge-based NDVI, NDVI, and Soil Adjusted Vegetation Index (SAVI) were calculated after pan-sharpening due to its benefits for detecting small vegetation patches (formulas in Table 1) [34].Ancillary data included an upland GIS layer created from a geomorphological map of Jamaica Bay [35] and Digital Elevation Model (DEM) derived from 2014 Topo-bathymetric Light Detection and Ranging [36].
Table 1.Vegetation Indices, including Worldview-2 Vegetation Index, Worldview-2 Water Index, Red Edge Vegetation Index, Normalized Difference Vegetation Index and Soil Adjusted Vegetation Index.

WVVI WVWI NDVI
Red Edge Vegetation Index SAVI Object neighborhoods, those objects that share a border with an object, and weights were calculated to determine the neighborhood difference of the mean spectral, textural and vegetation index attributes giving additional spatial context to the data [37].The final Worldview-2 image objects had 79 attributes including 3 spatial attributes, 18 texture attributes, 32 spectrally derived attributes, 7 elevation based, 18 vegetation index, and a binary upland variable (Table A3).The Quickbird-2 image objects had additional attributes including tasseled cap values but no Red Edge based NDVI.

Accuracy Assessment
The accuracy assessments were conducted for each classification by generating equalized random points.The number of points to generate was calculated with following equation [38].
where B is the Chi-squared distribution with 1 degree of freedom for the target error divided by the number of classes, ∏ i is the percent land cover of the most prevalent class and b is the desired confidence interval of that class.The calculation required over 750 test points to fulfill the accuracy assessment.The final test dataset was composed of 765 test points.The objects were classified by the user based on Worldview-2, Quickbird-2, and Google Earth historic imagery from each time period.Overall accuracy, the Kappa statistic, producer's accuracy, and user's accuracy were calculated to analyze the confusion matrix results [39,40].The study site was visited in 2014 and 2015 to verify the characteristics of the landscape and collect field reference data.The training samples and objects were extracted from Worldview-2 and Quickbird-2 imageries in combination with expert knowledge from the field visits.Land cover points were collected on each of the field visits.The point locations included areas in West Pond, Black Bank, Yellow Bar, JoCo, Elders Point, Canarsie Pol, and East High.The points were navigated to with a Trimble XH and the areas dominate vegetation community was verified.

Statistical Analysis
The finished classifications were utilized to determine change rates (ha/year) in the three time periods 2003-2008, 2008-2012 and 2012-2013.Jamaica Bay's unique salt marsh structure of individual islands led to their use for statistical analysis.The paired Wilcoxon signed rank test was utilized to test the differences between wrack extent throughout the Bay in 2008, 2012 and 2013 (Table A1).These extents were for each island for each year.The difference between percent change (∆%/year) of restored and natural salt marsh from 2012 to 2013 was tested with a student's t-test.Comparing all natural salt marshes to those where restoration was completed by the data acquisition in 2012.Before utilizing the t-test, normality was tested with the Shapiro-Wilkes statistic, which indicated normality could not be rejected with p values = 0.37 and 0.80 for restoration and natural, respectively.The natural salt marsh islands for all time periods were used to compare change rates (ha/year) with paired t tests.The time periods were tested with Shapiro-Wilkes for normality finding p values of 0.54, 0.29, 0.43 and 0.19 for 2003-2008, 2008-2012, 2012-2013 and 2003-2012 respectively.Linear regressions were used to understand the impact of salt marsh extent, latitude, and longitude on combined high marsh and both classes of S. alterniflora change rates (ha/year).Latitude and longitude were determined from the center point of each salt marsh island.

Results
The landscape was mapped accurately throughout the classification results (accuracies of 85.63, 85.2, 90.46 and 92.55 for 2003, 2008, 2012 and 2013, respectively).The overall accuracies were further analyzed by producer's accuracy and user's accuracy (Table 2).The 2003 data had an adequate overall accuracy of 85.6%, with vegetation classes exhibiting the lowest accuracies (Table 2).This led to a focus on comparing vegetated salt marsh and non-vegetated areas as most of the error was between the multiple classes of salt marsh.The three classes of vegetated (S. alterniflora classes, high marsh, and Phragmites), non-vegetated (water, mudflat, sand, and wrack) and upland were used for comparisons between periods unless stated otherwise.These three classes had overall accuracies of 96.09, 93.46, 93.46 and 96.73 for 2003, 2008, 2012 and 2013, respectively.

Wetland Change
The 2003 and 2013 classifications were compared to determine change between all classes (Table 3).From 2003 to 2013, 54.9 ha of sand, mudflat and water were converted into salt marsh.However, during that same period, 70.7 ha of high marsh and S. alterniflora were converted into sand, mudflat, or water.Salt marsh vegetation gains occurred in restoration sites, however, these were exceeded by losses in areas not subject to intervention (Figure 2).Elders Point East and West were restored during the study period, an example of restoration driven change in the Bay (Figure 3).West Pond was breached during Hurricane Sandy and areas of freshwater wetland and upland vegetation shifted to mudflat (Figure 4).From 2003 to 2013, 21 ha of salt marsh were lost, including both S. alterniflora classes, high marsh, and Phragmites.Smaller salt marshes such as Duck Point and Pumpkin Patch nearly disappeared (Figures 1 and 2).

Restored Islands: 2003-2013
Elders Point East and West were restored in 2006 and 2010, respectively [4].These islands were not being actively restored during the 2012-2013 period, however they did increase in salt marsh extent (Table A1).From 2012 to 2013, Yellow Bar, Rulers Bar, and Black Wall were the focus of significant restoration.In 2013, Yellow Bar added 8.0 ha of salt marsh vegetation, but had a negligible change in extent from 2003 to 2013.Yellow Bar's restoration process also added approximately 15 ha of mudflat, however, this does not account for the 32.5 cm higher tide in 2003 as determined from the Sandy Hook tidal gauge [41].From 2003 to 2008, restoration of Big Egg and Elders Point East were completed, resulting in increases in salt marsh extent of 4.0 and 9.5 ha, respectively.Big Egg subsequently lost 4.7 ha of salt marsh extent between 2008 and 2012.

Impact of Hurricane Sandy
West Pond (Figure 2) is a retention pond created during the construction of the Cross Bay Boulevard and an important resource for migratory birds [42] (Figure 4).Hurricane Sandy breached West Pond, resulting in saltwater intrusion into the freshwater environment [43].Prior to this breach, West Pond's wetlands were dominated by Phragmites australis.The area represents the most drastic change from Hurricane Sandy; alterations to the upland and freshwater wetlands are evident (Figure 4 and Table 4).Between 2003 and 2013, the JoCo site lost salt marsh vegetation going from 131.2 ha to 127.6 ha.However, from 2012 to 2013 vegetation increased (Table A1).This increase in vegetation was accompanied by a reduction in wrack across the Bay compared to both 2008 (W15 = 110, p < 0.003) and 2012 (W15 = 113, p = 0.0011).The area of wrack was reduced after Hurricane Sandy going from 2.2 ha to 0.5 ha.This in part accounts for the 3.6 ha increase of salt marsh vegetation observed in JoCo.The 2008 and 2012 classifications of JoCo had only 0.2 ha of overlapping wrack.

Restored Islands: 2003-2013
Elders Point East and West were restored in 2006 and 2010, respectively [4].These islands were not being actively restored during the 2012-2013 period, however they did increase in salt marsh extent (Table A1).From 2012 to 2013, Yellow Bar, Rulers Bar, and Black Wall were the focus of significant restoration.In 2013, Yellow Bar added 8.0 ha of salt marsh vegetation, but had a negligible change in extent from 2003 to 2013.Yellow Bar's restoration process also added approximately 15 ha of mudflat, however, this does not account for the 32.5 cm higher tide in 2003 as determined from the Sandy Hook tidal gauge [41].From 2003 to 2008, restoration of Big Egg and Elders Point East were completed, resulting in increases in salt marsh extent of 4.0 and 9.5 ha, respectively.Big Egg subsequently lost 4.7 ha of salt marsh extent between 2008 and 2012.

Impact of Hurricane Sandy
West Pond (Figure 2) is a retention pond created during the construction of the Cross Bay Boulevard and an important resource for migratory birds [42] (Figure 4).Hurricane Sandy breached West Pond, resulting in saltwater intrusion into the freshwater environment [43].Prior to this breach, West Pond's wetlands were dominated by Phragmites australis.The area represents the most drastic change from Hurricane Sandy; alterations to the upland and freshwater wetlands are evident (Figure 4 and Table 4).Between 2003 and 2013, the JoCo site lost salt marsh vegetation going from 131.2 ha to 127.6 ha.However, from 2012 to 2013 vegetation increased (Table A1).This increase in vegetation was accompanied by a reduction in wrack across the Bay compared to both 2008 (W 15 = 110, p < 0.003) and 2012 (W 15 = 113, p = 0.0011).The area of wrack was reduced after Hurricane Sandy going from 2.2 ha to 0.5 ha.This in part accounts for the 3.6 ha increase of salt marsh vegetation observed in JoCo.The 2008 and 2012 classifications of JoCo had only 0.2 ha of overlapping wrack.
JoCo salt marsh was the most stable during the time period analyzed (Table A1).The restoration salt marshes in 2012-2013 had a larger percentage increase of salt marsh vegetation than natural salt marshes (t 4 , p = 0.041).A2).However, 2012-2013 was the only period when a mean increase in natural salt marsh vegetation was observed.
We analyzed the natural salt marshes yearly change rates (ha/year) by linear regression for each time period.The only time period where salt marshes towards the eastern side of the Bay tended to gain vegetation was 2012-2013 (F 1,9 = 22.21, p < 0.002 and R 2 = 0.7116).Larger salt marshes in 2012 tended to gain vegetation in 2012-2013 (F 4,6 = 13.93,p < 0.036 and R 2 = 0.9028).From 2008 to 2012, larger salt marshes in 2008 tended to lose more vegetation (F 4,6 = 6.83, p < 0.011 and R 2 = 0.8199).From 2003 to 2008, no relationship was found between salt marsh extent and change (F 4,6 = 0.75, p = 0.41 and R 2 = 0.33).The switch in the direction of this relationship demonstrates different processes dominating the Bay between 2008-2012 and 2012-2013.

Accuracy Assessment
Confusion matrices were utilized to determine the performance of each of the classifications.The 2012 and 2013 classifications performed well in all vegetation classes of most interest including S. alterniflora and high marsh (Table 2).The lowest performing class was Phragmites, which is a difficult to classify land cover with overlap between many of the other classes spectrally and spatially.The 2003 and 2008 classifications had low salt marsh vegetation accuracy due to confusion between the salt marsh vegetation types.The 2003 and 2008 error was mitigated by focusing our analysis on change in vegetation not changes in particular types of vegetation.Overall, the Worldview-2 data were better suited for the specificity of this classification.

Discussion
Since the 1950s, salt marsh vegetation in Jamaica Bay has been in rapid decline and, in the early 2000s, restoration was deemed necessary to maintain the salt marsh.This study and past estimates of salt marsh change were compared to better understand vegetation trends.From 1989 to 2003, there was an estimated 13.4 ha of yearly loss [5].From 2003 to 2013, a yearly loss of 2.1 ha was observed.The long-term rate of salt marsh loss in the Bay slowed due in part to restoration, however, both tidal stage and nutrient inputs may have influenced this result.
The 2003 and 2013 data were collected at a tidal stage of −0.129 m and −0.454 m (North American Vertical Datum) [41].Between 2003 and 2013, the larger salt marsh islands appeared to gain vegetation in the interior and lose salt marsh on the edges (Figure 2).However, the difference in tidal stage of the data could be responsible for some vegetation increases between the two dates.The tidal stage of the 2012 data was −0.577 NAVD [41].The small tidal difference in 2012 and 2013 could result in less inundated vegetation in 2012.Tidal stage may have influenced the larger trends from 2003 to 2013, but was not a factor in the vegetation increase from 2012 to 2013.The impact of the tides on mapping salt marsh in Jamaica Bay should be further explored to account for this uncertainty.
Since the mid-2000s, the Bay has had a 30% reduction in nitrogen load [44].Nutrient enrichment in salt marsh systems can lead to creek bank collapse and conversion to mudflat [6].The Waste Water Treatment Plants in Jamaica Bay account for 89% of all nitrogen inputs into the Bay; due to the Bay's currents, the highest nitrogen concentrations were in the south and eastern sides of JoCo [15].The different responses of salt marshes in the Bay to nutrient enrichment was partly explained by lower elevation marshes having longer periods of inundation increasing decomposition and loss of organic matter [45].The nitrogen load reduction coincided with the slowing of salt marsh loss, however, the impact is unknown and in situ analysis would be necessary to explore this possible connection.

Restoration
In 2003, the first salt marsh restoration in Jamaica Bay began at Big Egg.The project utilized dredge sediment to increase marsh elevation and then S. alterniflora plugs were planted 50 cm apart [17].In 2006, Elders Point East's elevation was increased with dredge sediment and then vegetated with both plugs and hummock relocation, the removal followed by placement of the entire salt marsh platform on areas of restored elevation [4].The hummock relocation saves salt marsh that would be covered in dredge sediment, and provides vegetation to the restored area.In 2010, the elevation of Elders Point West was increased with dredge sediment and vegetated with a combination of hummock relocation, planting of high marsh species, seeding of S. alterniflora, and a test site with no planting [4].In early 2012, Yellow Bar was restored with dredge material and vegetated with a mix of hummock relocation and salt marsh seeding [46].In fall 2012, the elevation of Rulers Bar and Black Wall was increased with dredge sediment.In June 2013, a community effort added vegetation to these islands with plugs [46].This decade of restoration coincided with our study, and resulted in the evaluation of this methodology for understanding restoration.
Black Wall and Rulers Bar were restored between 2012 and 2013.These marsh islands showed no evidence of revegetation at the time of the 2013 mapping.The salt marsh vegetation of Black Wall and Rulers Bar was reduced from 2.7 to 1.2 ha while sand and mudflat increased from 11.2 to 18.2 ha.The loss of vegetation appeared to be connected with sediment deposition from restoration and lack of hummock relocation.However, the storm event could have exacerbated the loss.Rulers Bar lost nearly all salt marsh vegetation from 2012 to 2013 (Table A1).In June 2013, Plugs had been planted on Black Wall and Rulers Bar.However, the vegetation was sparse and classified as mudflat.
While restored salt marsh corresponded with a visual change, it may not represent a recovery of all the ecosystem services.Differences between natural and restored salt marshes include lower soil organic matter, insufficient nitrogen availability, stunted plant growth and increased susceptibility to herbivory [47].Field studies in Jamaica Bay have demonstrated some differences between restored and natural salt marshes, including a high percent of sand and less soil organic matter in the first 10 cm of soil [45].These differences and the unknown longevity of restored marshes are the reasons long-term monitoring is necessary.Big Egg and Elders Point East both demonstrated losses post restoration from 2008 to 2012, with a loss of 1.1 ha and 1.0 ha per year, respectively.Post-restoration losses demonstrate the need for further understanding of the underlying processes causing salt marsh loss in Jamaica Bay.The expected lifetime of a restored marsh could be estimated and used to inform management decisions.
Restoration planting occurred on Yellow Bar between the 2012 and 2013 data collections (Table 5).The restoration process added elevation and S. alterniflora to the northern area of the site.The restoration resulted in vegetation increasing from 18.2 to 26.3 ha.Elders Point's restoration was already complete in 2012, however, the combined vegetated extent of Elders Point East and West went from 13.0 to 14.5 ha.Restoration sites added vegetation in the post-storm growing season.Restoration sites did not appear to be negatively impacted by Hurricane Sandy.However, post-storm the Yellow Bar restoration required extensive repairs and replanting [46].The storm impacted Yellow Bar at an early stage of restoration, which led to a slowing of the process.However, considerable vegetation was gained in the post-storm growing season.

Hurricane Sandy
The response of salt marshes to storm events vary and include net elevation increases leading to vegetation growth [48] and accretion varying with distance to an inlet [49].The natural and restoration salt marshes responded differently to Hurricane Sandy.The analysis of natural salt marsh separated the restoration and storm impacts.In 2012-2013, larger salt marshes and those further from Rockaway inlet tended to gain vegetation.This is in agreement with past hurricane impacts which had a wide variation in sediment deposition and salt marsh response including edge erosion [50].The large salt marshes may have been less impacted by Hurricane Sandy, and captured more of the accompanying sediment pulse.
The response of vegetation in Jamaica Bay to Hurricane Sandy depended on the location and the ecosystem.Saltwater intrusion into freshwater ecosystems is a major source of storm event derived vegetation loss; evident in both coastal wetland environments [51] and forests [52].The survival and recovery of freshwater wetland vegetation depends on the species [53] and replanting of coastal forests can be limited by the increased soil salinity and herbivory [54].These long-term impacts emphasize the importance of monitoring the West Pond breach.Post-storm, both freshwater wetland and upland vegetation lost extent declining from 13.3 ha to 10.2 ha and 11.5 ha to 6.0 ha, respectively (Table 4).There were 2.9 ha of change from upland vegetation to freshwater wetland, which can be understood as a loss of vegetation biomass but not a complete loss of vegetation.Excluding those areas, 5.9 ha of freshwater wetland were lost.The majority of vegetation lost became mudflat.The loss of upland vegetation suggests the approximate extent of saltwater intrusion into the upland areas around West Pond.The environmental assessment of the site has resulted in the decision to close the breach and restore the freshwater wetland environment [55].This approach will create early successional habitat.Continued monitoring of West Pond is necessary to understand both the recovery of the freshwater ecosystem and unforeseen impacts of the management decision.

Wrack
Wrack is an important component of Jamaica Bay's landscape as persistent wrack deposits, for over four months of time, have a negative impact on the growth rate of all the principal marsh species [56].Storm events including hurricanes are understood as one of the causes of wrack accumulation [50].Mapping wrack accumulation pre-and post-storm enabled the evaluation of both the deposition and movement of wrack within Jamaica Bay.Post-storm there was less wrack on the salt marsh than in 2008 and 2012.When examining JoCo, it appears areas of wrack moved towards the center of the marsh island (Figure 5).If the same pattern occurred in islands with upland, wrack would have moved under the upland vegetation canopy.Throughout the Bay, most wrack became S. alterniflora, capturing the removal of wrack and regrowth of impacted salt marsh vegetation in the following growing season.These findings suggest recovery from wrack can be rapid, with storm events as a major driver in the deposition and distribution of the material throughout Jamaica Bay.

Long-Term Monitoring
The two most prevalent mapping protocols for wetland change analysis are the National Wetland Inventory (NWI) conducted by United States Fish and Wildlife Service (FWS) and Coastal Change Analysis Program (C-CAP) conducted by the National Oceanic and Atmospheric Administration (NOAA).These programs are each focused on mapping wetland change across the entire or the majority of the United States.The NWI is an estimate of the trends conducted by mapping a large number of randomly sampled plots which are interpreted based on aerial imagery [57].The methodology leads to trends in states or regions, however, these conclusions are not necessarily representative of rapidly changing sites such as Jamaica Bay.Between 2004 and 2008, the NWI estimated salt marsh increased in the Atlantic by 133 hectares, a negligible percent increase [3].Between 2003 and 2008, Jamaica Bay added 6.3 hectares of salt marsh vegetation, a 1.8% increase.The two estimates agree that an increase occurred, however, the NWI estimate lacks the precise location or magnitude of the restoration driven change.
The C-CAP utilizes Landsat, a 30 m spatial resolution sensor, to understand long-term change, however, accuracy reports showed confusion between water, consolidated shore, and emergent marsh [58].From 2001 to 2010, C-CAP's estuarine emergent wetland class maintained an extent of 674 ha in Jamaica Bay.During that time frame, Big Egg and Elders Point East were restored, which had no discernable change in the extent of estuarine emergent wetland class.Remote sensing with high resolution imagery has been successfully utilized for monitoring restoration [59].The coarse spatial and temporal resolution of C-CAP makes understanding storm events or restoration in Jamaica Bay difficult.Localized solutions are necessary for capturing a restoration baseline and then mapping at an appropriate temporal resolution to understand shifts between vegetation communities and long-term restoration trajectories.
Salt marsh losses are increasingly driven by sea level rise and high water events causing migration of S. alterniflora into areas previously composed of high marsh [60].In order to understand these shifts between vegetation communities, a specialized high resolution classification is necessary.When conducting analysis over large areas C-CAP and NWI programs are invaluable.However, a specialized protocol is preferable when presented with single study site and unique management issues.
The regular collection of satellite imagery is necessary for long-term monitoring.This can have a prohibitive cost, when using very high resolution satellite data.This study's five-year data collection interval and additional data collected following the storm event was adequate for understanding both the decadal trends and Hurricane Sandy's impact.Jamaica Bay is representative

Long-Term Monitoring
The two most prevalent mapping protocols for wetland change analysis are the National Wetland Inventory (NWI) conducted by United States Fish and Wildlife Service (FWS) and Coastal Change Analysis Program (C-CAP) conducted by the National Oceanic and Atmospheric Administration (NOAA).These programs are each focused on mapping wetland change across the entire or the majority of the United States.The NWI is an estimate of the trends conducted by mapping a large number of randomly sampled plots which are interpreted based on aerial imagery [57].The methodology leads to trends in states or regions, however, these conclusions are not necessarily representative of rapidly changing sites such as Jamaica Bay.Between 2004 and 2008, the NWI estimated salt marsh increased in the Atlantic by 133 hectares, a negligible percent increase [3].Between 2003 and 2008, Jamaica Bay added 6.3 hectares of salt marsh vegetation, a 1.8% increase.The two estimates agree that an increase occurred, however, the NWI estimate lacks the precise location or magnitude of the restoration driven change.
The C-CAP utilizes Landsat, a 30 m spatial resolution sensor, to understand long-term change, however, accuracy reports showed confusion between water, consolidated shore, and emergent marsh [58].From 2001 to 2010, C-CAP's estuarine emergent wetland class maintained an extent of 674 ha in Jamaica Bay.During that time frame, Big Egg and Elders Point East were restored, which had no discernable change in the extent of estuarine emergent wetland class.Remote sensing with high resolution imagery has been successfully utilized for monitoring restoration [59].The coarse spatial and temporal resolution of C-CAP makes understanding storm events or restoration in Jamaica Bay difficult.Localized solutions are necessary for capturing a restoration baseline and then mapping at an appropriate temporal resolution to understand shifts between vegetation communities and long-term restoration trajectories.
Salt marsh losses are increasingly driven by sea level rise and high water events causing migration of S. alterniflora into areas previously composed of high marsh [60].In order to understand these shifts between vegetation communities, a specialized high resolution classification is necessary.When conducting analysis over large areas C-CAP and NWI programs are invaluable.However, a specialized protocol is preferable when presented with single study site and unique management issues.
The regular collection of satellite imagery is necessary for long-term monitoring.This can have a prohibitive cost, when using very high resolution satellite data.This study's five-year data collection interval and additional data collected following the storm event was adequate for understanding both the decadal trends and Hurricane Sandy's impact.Jamaica Bay is representative of the future for increasingly populated coastal communities worldwide, necessitating continued remote sensing monitoring of the impact of urbanization on the Bay's salt marsh.Long-term monitoring requires additional exploration of the impact that multiple sensors have on change analyses.The switch from Quickbird-2 to Worldview-2 could be partly responsible for the change seen from 2008 to 2012.Quantifying this impact is a necessary step as we proceed into the third decade of commercially available very high resolution satellite imagery.

Conclusions
This study reiterates the importance of continuing salt marsh monitoring with high spatial resolution satellite data within Jamaica Bay.This long history of monitoring allows an understanding of salt marsh change, restoration, and natural disturbance.Despite 10 years of restoration, salt marshes in Jamaica Bay continue to decline, though the yearly rate of loss slowed from 13.4 ha during from 1989 to2003 to 2.1 ha during 2003-2013 [5].While Quickbird-2 data resulted in an adequate classification, a single scene of Worldview-2 was better suited to discern between salt marsh vegetation classes.The analysis of individual marsh islands elucidates the varied responses over the last 10 years such as the stabilization of JoCo and the near complete loss of Pumpkin Patch.
Hurricane Sandy influenced both the salt marsh and freshwater wetlands of Jamaica Bay.The 2013 growing season in the Bay appeared to be impacted by the hurricane.The greatest change in Jamaica Bay attributed to Hurricane Sandy was the breach of West Pond, which caused a die-off of both upland and freshwater wetland vegetation within this important bird habitat.In total 8.6 ha of vegetation was lost around West Pond.Continued monitoring of the site is necessary to understand the long-term recovery of this area.While outside of our study's target salt marsh protocol, the classification and change analysis was robust enough to interpret this landscape's change.
The vegetation loss in Jamaica Bay slowed over the study period.The salt marsh extent increased from 2012 to 2013, which can partly be accounted for by the restoration of Yellow Bar, movement of wrack off the salt marsh, and differences in phenology between the two dates.Significant vegetation loss occurred in smaller salt marsh islands and the West Pond area.
The dynamic nature and complexity of coastal wetlands makes monitoring with high temporal resolution important and necessary to understand change.This study demonstrates the feasibility of object-based classification and change detection using Worldview-2 data for mapping, monitoring and understanding salt marsh change in Jamaica Bay.The approach could be expanded to other coastal systems, with a focus on areas of restoration or periods of change.The decline of the salt marsh habitats in the Jamaica Bay is of concern from an ecological stand point and for the important role that coastal wetlands have in mitigating storm surge [61].Future research should explore the impact of tidal stage on vegetation extent within the salt marsh environments of Jamaica Bay.

Figure 1 .
Figure 1.The study area of Jamaica Bay, NYC, includes salt marsh islands as labeled on top of the pseudo color display of 2012 Worldview-2 imagery (NIR-1, G, B in RGB).Field photos illustrate: (a) the transition from Phragmites australis to salt marsh; (b) Isolated S. alterniflora patch; and (c) S. alterniflora 50%-100% cover.Salt marshes that have been restored at some point are indicated by a white border.

Figure 1 .
Figure 1.The study area of Jamaica Bay, NYC, includes salt marsh islands as labeled on top of the pseudo color display of 2012 Worldview-2 imagery (NIR-1, G, B in RGB).Field photos illustrate: (a) the transition from Phragmites australis to salt marsh; (b) Isolated S. alterniflora patch; and (c) S. alterniflora 50%-100% cover.Salt marshes that have been restored at some point are indicated by a white border.

Figure 3 .
Figure 3. Salt marsh of Elders Point East and West for 2003, 2008, 2012 and 2013.Figure 3. Salt marsh of Elders Point East and West for 2003, 2008, 2012 and 2013.

Figure 3 .
Figure 3. Salt marsh of Elders Point East and West for 2003, 2008, 2012 and 2013.Figure 3. Salt marsh of Elders Point East and West for 2003, 2008, 2012 and 2013.

Figure 4 .
Figure 4. Vegetation change from 2012 to 2013 of the West Pond area.

Figure 4 .
Figure 4. Vegetation change from 2012 to 2013 of the West Pond area.

Table 3 .
Change between 2003 and 2013 (ha).Areas that had no change between the two dates are in grey.

Table 4 .
Change of land cover classes between 2012 and 2013 for West Pond area.

Table 4 .
Change of land cover classes between 2012 and 2013 for West Pond area.