1
Development of a Multimetric Index for Integrated Assessment of Salt Marsh Ecosystem Condition
Jessica L. Nagel1, Hilary A. Neckles2*, Glenn R. Guntenspergen1, Erika N. Rocks3, Donald R. Schoolmaster Jr.4, James B. Grace4, Dennis Skidds5 and Sara Stevens5
1U.S. Geological Survey, Patuxent Wildlife Research Center, 12100 Beech Forest Rd, Laurel, MD 20708 USA
2U.S. Geological Survey, Patuxent Wildlife Research Center, 196 Whitten Road, Augusta, ME 04330 USA
3National Park Service, Northeast Coastal and Barrier Network, 54 Elm Street, Woodstock, VT 05091 USA
4U.S. Geological Survey, Wetland and Aquatic Research Center, 700 Cajundome Blvd., Lafayette, LA 70506 USA
5National Park Service, Northeast Coastal and Barrier Network, University of Rhode Island Coastal Institute in Kingston, 1 Greenhouse Rd., Kingston, RI 02881 USA
*corresponding author: e-mail: [email protected]; tel: (207) 626-6619; fax: (207) 622-8204
1
2
1
2
3456789
1011121314151617181920212223242526272829303132333435363738394041
2
Abstract
Tools for assessing and communicating salt marsh condition are essential to guide decisions
aimed at maintaining or restoring ecosystem integrity and services. Multimetric indices (MMIs)
are increasingly used to provide integrated assessments of ecosystem condition. We employed a
theory-based approach that considers the multivariate relationship of metrics with human
disturbance to construct a salt marsh MMI for five National Parks in the northeastern US. We
quantified the degree of human disturbance for each marsh using the first principal component
score from a principal components analysis of physical, chemical, and land-use stressors. We
then applied a metric-selection algorithm to different combinations of about 45 vegetation and
nekton metrics (e.g. species abundance, species richness, and ecological and functional
classifications) derived from multi-year monitoring data. While MMIs derived from nekton or
vegetation metrics alone were strongly correlated with human disturbance (r-values from -0.80 to
-0.93), an MMI derived from both vegetation and nekton metrics yielded an exceptionally strong
correlation with disturbance (r = -0.96). Individual MMIs included from one to five metrics. The
metric-assembly algorithm yielded parsimonious MMIs that exhibit the greatest possible
correlations with disturbance in a way that is objective, efficient, and reproducible.
Keywords: Multimetric index; Salt marsh; Ecosystem assessment; Ecological indicators;
Coastal wetlands
3
4
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
3
Introduction
Coastal wetlands occur at the interface between uplands and oceans within a narrow
range of elevation influenced by tidal characteristics. Salt marshes occupy as much as 45,000
km2 world-wide (Greenberg et al. 2006). These coastal wetlands provide essential livelihood
services to millions of people, as well as critical regulating services such as maintenance of water
quality, protection from storms and erosion, and carbon sequestration valued at US $10,000 per
hectare per year (Barbier et al. 2011; Hopkinson et al. 2012; Möller et al. 2014). However, as
much as half of global coastal wetlands have been lost due to human conversion for agriculture
and other anthropogenic changes in land use (Pendleton et al. 2012), and climate change,
declining water quality, and changes in sediment delivery rates associated with human activity
continue to affect the world’s remaining wetlands (Kirwan and Megonigal 2013).
In the US, direct human modification of coastal wetlands has led to extensive wetland
degradation (Kennish 2001; Bertness et al. 2002) through activities such as ditching and drainage
of the marsh platform; reduction of tidal exchange by dikes, roads, and water control structures;
disposal of dredge spoil; introduction of invasive species; and discharge of nitrogen, phosphorus,
oil and other contaminants (Roman et al. 2000). The impacts of such disturbances on salt marsh
integrity and their ecosystem services are far reaching (Daiber 1986; Wolfe 1996; Adamowicz
and Roman 2005). Salt marsh management and conservation decisions are increasingly focused
on maintaining and improving ecological integrity (MEA 2005). As a consequence, approaches
to monitor and assess ecosystem condition have gained widespread attention (e.g. Simenstad and
Thom 1996; Short et al. 2000; Wigand et al. 2001; Neckles et al. 2002, 2015).
5
6
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
4
Considerable research has been devoted to development of both quantitative and
functional approaches to assess ecosystem integrity and track ecosystem condition over time. In
recent decades, multimetric indices (MMIs) have been used to integrate numerous characteristics
of complex ecosystems into one-dimensional indices of ecosystem condition (Karr 1999; Hering
et al. 2006). By combining multiple indicators of human disturbance into a composite score, an
MMI may be more sensitive to disturbance than any of its individual component metrics (Lee et
al. 2004; Schoolmaster et al. 2012). MMIs for wetland ecosystem assessment have been based on
specific biota, such as fish (Simon 1998), zooplankton (Lougheed and Chow-Frazer 2002), and
plants (Lopez and Fennessy 2002; Mack 2004), and they have also been based on combinations
of taxa (plants, invertebrates, and fish; Wilcox et al. 2002; James-Pirri et al. 2014) or ecosystem
components (plants, various fauna, and physical environmental features; Wigand et al. 2011;
Langman et al. 2012; Miller et al. 2016). As development of wetland MMIs has proliferated,
they have been applied to increasingly larger scales. MMIs have been developed worldwide for
regional assessments of inland (e.g., Ortega et al. 2004; Boix et al. 2005; Miller et al. 2006;
Mereta et al. 2013) and coastal wetlands (e.g. Wigand et al. 2011; Staszak and Armitage 2013;
Nestlerode et al. 2014) and are integral to reporting results of the National Wetland Condition
Assessment (NWCA) in the conterminous US (USEPA 2016).
Traditionally, MMI development has involved exploratory analysis of ecosystem
attributes or metrics measured at sites representing a range of conditions, from a reference
standard to sites that are subject to high levels of anthropogenic disturbance (Karr 1999; Karr
and Chu 2000). Expert knowledge or various ordination techniques may be used to reduce the
number of candidate metrics, organize them into classes, and rank their importance in explaining
the variation among sites, and metrics then may be selected for inclusion in the MMI based on
7
8
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
5
the strength of bivariate relationships between disturbance and individual responses. However,
interactions among ecosystem stressors, environmental covariates, and ecological responses may
introduce unexplained variability into metric data sets, reducing the sensitivity of individual
metrics to perceived disturbance factors and obscuring determination of the “best” metrics to
include in the MMI (Schoolmaster et al. 2012; 2013). In such cases, Schoolmaster et al. (2012)
demonstrated that evaluating how metrics respond to disturbance as a group can improve the
potential to generate a maximally sensitive MMI. This may be particularly relevant for salt
marsh ecosystems, where relationships among agents of change, physical and chemical
environmental stressors, and ecological responses are inherently complex.
We applied a multivariate, algorithmic approach for MMI development to salt marsh
ecosystems in the northeastern US. As described originally by Schoolmaster et al. (2012; 2013),
this is an empirical method of MMI construction from a set of candidate metrics, in which the
unique contribution of each metric toward explaining the variation in human disturbance is
evaluated within the context of all candidates. Through an objective, transparent, and
reproducible process for metric selection, an MMI is produced that is statistically robust and
sensitive given the set of available data. This method of MMI construction was tested previously
on freshwater wetlands (Schoolmaster et al. 2013). We applied the algorithmic approach to salt
marshes within US National Park Service units distributed between Massachusetts and
Maryland. Our study shows the broad utility of this approach for developing robust assessment
tools for salt marsh ecosystems.
9
10
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
6
Methods
Salt Marsh Monitoring Data
Salt marsh monitoring data were collected from 33 marsh study units (MSUs) in five
National Parks within the NPS Northeast Coastal and Barrier Network (NCBN) along the
northeastern coast of the US (Figure 1, Table 1). MSUs were established at Cape Cod National
Seashore (NS), Sagamore Hill National Historic Site (NHS), Fire Island NS, Gateway National
Recreation Area (NRA), and Assateague Island NS. Individual MSUs were delineated to
encompass the contiguous marsh area bounded by geographic features, such as tidal creeks, salt
marsh-upland borders, grid ditches, shoreline, and other distinguishable features, using GIS. For
all parks except Cape Cod NS, all areas of salt marsh within the park were identified, all
potential MSUs were delineated, and the final sample of MSUs at each park was then randomly
selected from this population. At Cape Cod NS, we used existing MSUs that had been selected
as representative of the park’s salt marsh resources. MSU size at most parks ranged from 5-15
ha, but at Cape Cod NS the MSUs were generally larger (up to 168 ha).
MMI development was based on marsh vegetation and nekton (i.e., fish and free-
swimming crustaceans) monitoring data collected from 2008 – 2013 as part of the network-wide
vital signs monitoring program (Stevens et al. 2005). Monitoring occurred biennially in most
MSUs, yielding three years of data during this time period; however, for some MSUs only two
years of data were available (Table 1). Vegetation data were collected in all 33 MSUs and nekton
data were collected at only 24 of the MSUs (Table 1). To allow comparisons of MMIs developed
from different sized datasets (24 vs 33 MSUs) and metric types (vegetation only vs vegetation +
nekton), we generated four separate datasets: vegetation data alone from 33 MSUs (V33);
11
12
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
7
vegetation and nekton data combined from 24 MSUs (VN24); vegetation data alone from 24
MSUs (V24); and nekton data alone from 24 MSUs (N24). In addition, salt marsh monitoring
data were collected irregularly between 2003 and 2007 in five of the MSUs located in three parks
as part of a pilot program (Cape Cod NS, Gateway NRA, and Sagamore Hill NHS, Table 1). We
applied the MMIs to the data collected during the pilot program to evaluate changes in salt marsh
condition as revealed by different MMIs.
The percent cover of marsh vegetation and other cover types (e.g., water, bare ground,
wrack, and litter) was measured once during the peak of the growing season (early July through
early September) per sampling year. From 2008 – 2013, cover was measured using the Braun-
Blanquet method (Kent and Coker 1994) within approximately 50 randomized square-meter (1m
x 1m) quadrats along 10 randomized transects per MSU. Percent cover of each species was
estimated visually within cover classes (absent; <1 %; 1-5 %; 6-25 %; 26-50 %; 51-75 %; 76-
100 %), and we converted cover classes to class midpoints for analytical purposes. Before 2008,
cover was measured using either the Braun-Blanquet method or the point-intercept method, in
which cover was recorded as the proportion of points intercepted by each species in a grid of 50
or 100 points per square-meter quadrat (Roman et al. 2001). In addition to percent cover of
vegetation, the height of the introduced common reed, Phragmites australis, was measured when
present, as the distance from the marsh surface to the tallest point of the plant. If there were 20 or
fewer stems within a quadrat, all stems were measured; if there were more than 20 stems, all
stems within a randomly selected quarter of the quadrat were measured.
Nekton populations were monitored in shallow pools, creeks, and ditches within each
MSU twice during the growing season: once in early summer (mid-June through July) and once
in late summer-early fall (August through September). During each sampling event, between 5
13
14
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
8
and 20 stations were sampled per MSU, depending on the availability of pool, creek, or ditch
habitat. Wherever possible, a mix of habitats (e.g., pools and creeks) was sampled. Sampling
stations were established randomly around the perimeter of pools or along the length of creeks
and ditches. Shallow pools and creeks greater than 1 m wide were sampled using throw traps,
and ditches were sampled using ditch nets (James-Pirri et al. 2012). Each individual captured
was counted and identified to species and up to 15 individuals of each species were measured for
total length (fish and shrimp) or carapace width (crabs) at each station. Species counts were
converted to density (number m-2) by dividing by sampled area.
Candidate Response Metrics
We derived a wide variety of vegetation and nekton response metrics for potential
inclusion in each MMI. Our initial list included 50 vegetation metrics and 28 nekton metrics at
species and community levels that were suggested by other assessments of salt marsh ecosystem
condition (Konisky et al. 2006; Moore et al. 2010; Wigand et al. 2011; Neckles et al. 2013).
Species-based vegetation metrics were derived from estimates of percent cover or frequency of
seven species or species groups (Spartina alterniflora, Spartina patens, Distichlis spicata, Iva
frutescens, Phragmites australis, Salicornia spp., Scirpus and Schoenoplectus spp.) that are
broadly representative of different salt marsh environments. Species-based nekton metrics were
calculated using estimates of relative abundance (percent of total density) of five common native
species or species groups (Fundulus spp., Cyprinodon variegatus, Palaemonetes spp., Anguilla
rostrata, Menidia spp.) and one invasive species (Carcinus maenas). Species-level metrics also
included the height of the invasive plant Phragmites australis as an indicator of plant vigor (cf.
15
16
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
9
Howard and Mendelssohn 1999) and length of the common fish species Fundulus heteroclitus as
an indicator of life-history stage (Raposa et al. 2003).
Vegetation and nekton community metrics described the abundance and proportional
distribution of species within key ecological and taxonomic categories. Each plant species was
classified according to salt marsh zone (low marsh; high marsh; pannes, pools, and creeks; salt
marsh border; brackish border; or upland), salinity tolerance (high, medium, or low), and native
status (native or introduced) using information available from various sources (Duncan 1974;
Eleuterius 1990; Schafale and Weakley 1990; Haines and Vining 1998; Magee and Ahles 1999;
Silberhorn 1999; Shumway and Banks 2001; Swain and Kearsley 2001; Edinger et al. 2002;
Harrison 2004; Konisky et al. 2006; Tiner 2009; Sperduto 2011; Fleming and Patterson 2013;
NYNHP 2013; USDA 2014; Supplemental Table 1). Percent cover of vegetation community
classes in each quadrat was calculated by summing over the percent cover estimates of all
species within each community. Similarly, nekton species were classified according to broad
taxonomic group (e.g., fish, shrimp, crab, etc.), migrational status (resident or transient), and
native status (native or introduced) using various sources (NYNHP 2013; Chesapeake Bay
Program 2014; FishBase 2014; IRL Species Inventory 2014; Supplemental Table 2), and density
and relative abundance of nekton communities in each sample were derived by summing over all
species within each community. Total and relative (percent of total species) species richness
metrics for each vegetation and nekton community category were calculated at both the sample
and the MSU level. An additional nekton community metric, species dominance, was derived by
calculating the minimum number of species forming greater than 90% of abundance in the
sample (Hughes et al. 2002).
17
18
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
10
Temporal means of most species- and community-based metrics derived from the 2008-
2013 monitoring data were derived by averaging the annual MSU means. Frequency metrics
(proportion of occurrences relative to total sample number) were calculated both as the average
of annual means and as the total frequency of occurrence across all samples and years (2008-
2013). Within each dataset, metrics that contained a large proportion of zero values (> 66 % of
the total values within a metric) were omitted from further analysis (Stoddard et al. 2008;
Schoolmaster et al. 2013). In addition, metrics that duplicated another metric were removed, as
they would have contributed essentially identical information to the MMI. We identified
duplicate metrics by high values of inter-metric correlation (r > 0.95 or r < -0.95; Schoolmaster
et al. 2013). For each duplicate pair we retained the metric that was most interpretable in the
context of salt marsh ecosystem condition. Following these metric-reduction steps, the V33 data
set contained 27 vegetation metrics (Table 2), the V24 data set contained 28 vegetation metrics
(Table 2), the N24 data set contained 18 nekton metrics (Table 3), and the VN24 data set
contained 27 vegetation (Table 2) and 18 nekton metrics (Table 3).
Human Disturbance Index
We derived a Human Disturbance Index (HDI) to use for selecting response metrics for
inclusion in the MMI. We used existing aerial imagery and the 2006 National Land Cover
Database (NLCD; Fry et al. 2011) to quantify human impacts to each MSU in terms of physical,
chemical, and land-use stressors (Table 4). Categorical landscape metrics and classes were
adopted from the New England Rapid Assessment Method for assessing salt marsh condition
(Carullo et al. 2007) and aerial photographs of marshes representing each category were used as
a guide in classification. The metrics were scored by a trained technician using the most recently
19
20
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
11
available orthophotography in ArcMapTM 10.2 in consultation with park staff familiar with the
MSUs. Continuous metrics quantifying the amount of developed land cover within 150-m and 1-
km buffers of each MSU were calculated from NLCD data using ArcGISTM 10.2. The percent of
developed land area in each buffer zone was derived by aggregating developed land cover
classes in the NLCD (Developed Open Space, Developed Low Intensity, Developed Medium
Intensity, Developed High Intensity, and Barren Land) into a single category. Values were then
relativized to the size of the MSU by multiplying by (buffer area / MSU area) (Table 4). Only
one MSU (CACO_10) had measurable agricultural land in the buffer zone (relativized percent of
agricultural land in the 150-m buffer = 0.8 %), thus this land cover was not included in the HDI
formulations.
Categorical and continuous disturbance metrics were converted to an ordinal scale
following scoring criteria outlined in the New England Rapid Assessment Method (Carullo et al.
2007). Categorical metrics were assigned a rank score in which 0 represented undisturbed
conditions and 6 represented the most disturbed condition. Continuous variables were converted
to ranked values using three break-points in the distribution: values ≤ (mean - ½SD) were
assigned 0; values > (mean - ½SD) and ≤ mean were assigned 2; values > mean and ≤ (mean +
½SD) were assigned 4; and values > (mean + ½SD) were assigned 6.
A Principal Components Analysis (PCA) of the seven ranked disturbance metrics was
used to derive two HDIs: HDI 24 corresponded to the metric datasets using 24 MSUs (VN24,
V24, N24) and HDI 33 corresponded to the dataset using 33 MSUs (V33). A Spearman rank
correlation matrix of the disturbance metrics was used as input for the PCA and we evaluated
the Principal Components (PCs) by examining the eigenvectors, scree test, variance explained,
and comprehensibility (Kachigan 1982). The first Principal Component score (PC-1) for each
21
22
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
12
site was used as the basis of the HDI score (Falcone et al. 2010) as PC-1 accounted for the
largest part of the total variance inherent in the data. The PC-1 scores were transformed to HDI
by first reordering if necessary so that more positive values represented higher human
disturbance, and then rescaling to the HDI range of 0 (minimal human disturbance) to 100
(highest human disturbance). This type of disturbance index has been used successfully for
watershed and wetland assessment in the Great Lakes Basin (Brazner et al. 2007; Danz et al.
2007), in New England (Wigand et al. 2011), and throughout the western US (Falcone et al.
2010).
MMI Construction
We applied the MMI construction algorithm to the four sets of candidate metrics derived
from the 2008-2013 salt marsh monitoring data (V33, V24, VN24, N24; Table 2; Table 3).
Because the metrics have different units, all metrics were rescaled to unitless measures with
similar ranges using a continuous scoring method (Blocksom 2003):
m scaled=m-LU-L
where m = metric; L = 2.5 percentile values of m; U = 97.5 percentile values of m. Values of m
that were < L or > U were set to L or U, respectively. Metrics that were positively correlated
with HDI were reflected about their midpoints to ensure that the MMI would be negatively
correlated with disturbance. MMI assembly then followed a stepwise metric selection process to
maximize the likelihood of predicting the observed HDIs. Details of the MMI assembly
algorithm are provided by Schoolmaster et al. (2012; 2013). In general, each candidate metric
was used as a starting point, m1. The initial ml was added site-by-site to each of the remaining
metrics, mj, to find the combination m1 + mj yielding the strongest negative correlation with HDI.
23
24
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
13
This metric combination was then added to each of the remaining metrics, mj, to find the new
combination + mj yielding the strongest negative correlation with HDI. This metric selection
process was repeated until the addition of another mj failed to improve the ability to predict HDI
(likelihood ratio test, p > 0.05). This entire assembly process was repeated using each metric as
the starting point m1, thereby yielding as many candidate MMIs as there were candidate metrics.
Each candidate MMIj thus represented a model, HDI = β0 + β1 MMIj , with k parameters
(effectively, the number of metrics in MMIj + 2). The Akaike Information Criterion (AIC) of
each model was then calculated as the basis for comparing MMIs and selecting those with the
highest capacity to predict disturbance. For each candidate MMI, the difference between its AIC
and that of the MMI with the lowest AIC was determined (i.e., ΔAIC). All MMIs with ΔAIC < 2
were retained as having substantial support. If more than one MMI resulted for any dataset, a
single, model-averaged MMI was derived from the sum of the matrix products of the candidate
MMI and the Akaike weights for each candidate MMI (Schoolmaster et al. 2013). MMI scores
for each MSU were calculated by averaging the rescaled metric values for that MSU. The MMI
construction algorithm was written and run in the R statistical computing environment (R Core
Team 2016).
MMI Application
To explore how MMI formulation influenced perceived trends in salt marsh condition
over time, we applied the MMIs derived from the VN24 and V24 datasets to data collected
within five of the MSUs during a pilot study between 2003 and 2007 (Table 1). Values for
metrics included in each MMI were calculated for each year that data were available for a given
25
26
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
14
MSU. Metric values were then rescaled and combined as described above to generate MMI
scores for each year of pilot data.
Results
HDI Scores
The PCAs of the ranked disturbance metrics revealed a substantial drop in eigenvalue
magnitude, or the amount of variation explained by each PC, between the first and second PC for
each dataset (Table 5); PC-1 explained from 38 % (HDI 33) to 41 % (HDI 24) of the total
variation, and PC-2 explained 21 % (HDI 33) or 20 % (HDI 24). The metrics with the highest
loadings relative to the first PC were Fill / Fragmentation, Tidal Flushing and Tidal Restriction
(Table 5).
The HDI scores for MSUs were similar between the 24- and 33-MSU datasets (Figure 2).
There was substantial variation in HDI among MSUs. Approximately two thirds of the MSUs
exhibited relatively low HDI scores (i.e., HDI< 30 on a scale of 0 to 100). The MSUs with the
lowest HDI scores occurred on the barrier islands of Fire Island NS and Assateague Island NS.
Sites with highest HDI scores were located at Cape Cod NS (CACO_1, CACO_2, CACO_3,
CACO_6), the Long Island coast of Fire Island NS (FIIS_1, FIIS _2, FIIS _3) and Gateway NRA
(GSH_1, GSH_3).
MMI Development
Application of the MMI assembly algorithm to the vegetation-only metric sets yielded
individual MMIs consisting of three to five metrics. We refer to MMIs based on the metric set of
origin (e.g., MMI V33 was based on the V33 dataset). Single-letter subscripts refer to candidate
27
28
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
15
MMIs with ΔAIC < 2, and the subscript wt refers to the model-weighted MMI. All of the
vegetation-only MMIs were strongly correlated with HDI (r = -0.89 to -0.93; Table 6). A single
MMI (MMI V33) was constructed based on 33 MSUs, and two candidate MMIs with ΔAIC < 2
(MMI V24a, MMI V24b) were constructed based on 24 MSUs. Averaging the two MMI V24s
using model weights slightly improved the correlation with HDI (MMI V24wt; Table 6). There
was strong similarity among the metric sets included in the vegetation-only MMIs, with four
metrics identical between MMI V33 and MMI V24wt (Table 6). Metrics associated with the
prevalence of Phragmites australis and upland border communities were positively correlated
with HDI, and metrics describing the frequency of high-marsh species and the proportion of
species with high salinity tolerance were negatively correlated with HDI.
When nekton metrics alone were used as the basis for MMI construction, the assembly
algorithm identified five candidate MMIs composed of seven unique metrics, with from one to
three metrics each (MMI N24a-e; Table 7). Correlation of the nekton-only MMIs with HDI
ranged from -0.80 to -0.84. The model-weighted MMI composed of all seven metrics had a
correlation of -0.85 with HDI (MMI N24wt, Table 7). Metrics describing the relative abundance
of resident shrimp and the number of introduced species were most positively correlated with
HDI, and metrics associated with the relative abundance and number of species of resident fish
were most negatively correlated with HDI.
When both vegetation and nekton metrics were included in the MMI construction, a
single MMI (MMI VN24) resulted that was highly correlated with HDI (r = -0.96). MMI VN24
consisted of four metrics (two vegetation, two nekton), three of which had been selected during
the vegetation- and nekton-only constructions. The individual metrics were correlated with HDI
as follows: percent cover of native vegetation species, r = -0.44; percent cover of brackish border
29
30
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
16
vegetation species, r = 0.64; relative abundance of resident shrimp, r = 0.49; relative abundance
of resident fish, r = -0.79. Including both metric types in MMI construction resulted in the
strongest correlation with HDI in comparison to either vegetation or nekton metrics alone
(Figure 3).
In many cases, the relative differences in MMI scores across the MSUs were similar
between the individual MMI formulations (Figure 4). For example, all of the MSUs at
Assateague Island National Seashore and several at Fire Island National Seashore (FIIS_5, 6, 7,
9) generally had the highest MMI scores regardless of MMI formulation. However, some distinct
differences in the pattern of MMI scores across MSUs did emerge among the MMIs. In
particular, CACO_2 had among the lowest MMI scores based on MMI N24, but among the
highest scores based on the other MMI formulations. In addition, the nekton-only MMI
formulation (N 24) yielded the lowest MMI scores at SAHI_1 and FIIS 1, 2, and 3, but yielded
the highest MMI scores for most of the ASIS MSUs.
MMI Application to Earlier Years
For two of the MSUs with earlier-years’ monitoring data available, the MMI formulation
used to derive MMI scores did not affect the perceived trends in salt marsh condition over time:
within GSH_1, the MMI VN24 and the MMI V24 scores both decreased from 2003 to 2008-
2013, and within SAHI_1, they increased from 2004 to 2008-2013 (Figure 5). In contrast, for
three of the MSUs, the temporal pattern of MMI scores varied with the MMI formulation. Within
CACO_1 and CACO_2, MMI VN24 scores decreased from 2006 to 2008-13 while MMI V24
scores remained fairly stable over the same time period, and within CACO_3, MMI VN24
increased from 2005 to 2008-13 while MMI V24 decreased (Figure 5).
31
32
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
17
Discussion
The MMI construction algorithm (Schoolmaster et al. 2012; 2013) used in this study for
NPS northeastern salt marshes produced MMIs with very strong correlations to human
disturbance. This algorithmic approach identifies the metrics that, in combination, best reflect the
multivariate response to disturbance. The vegetation metrics that were positively correlated with
disturbance (frequency or percent cover of the introduced species Phragmites australis and
brackish border species, Table 6) are indicative of vegetation communities associated with lower
salinity (Supplemental Table 1) and likely reflected the degree to which tidal restriction reduces
integrity of northeastern salt marshes by reducing salt water influx (cf. Roman et al. 2000). In
addition, Phragmites australis tends to form monotypic stands that reduce wetland diversity and
compromise ecosystem services (Zedler and Kercher 2004). Vegetation metrics describing the
frequency, abundance, or percent composition of Distichlis spicata, Spartina patens, high
salinity tolerant species, and native species were all negatively correlated to disturbance (r = -
0.44 to -0.80, Table 6), and are indicative of natural salt marsh communities. The direction and
magnitude of the correlations of the selected nekton metrics to HDI (Table 7) were also
consistent with observations in northeastern salt marshes. For example, several studies have
shown that resident fish abundance is reduced and shrimp abundance is increased in salt marshes
in densely populated and developed watersheds (James-Pirri et al. 2014) or those with a high
degree of fragmentation or subsidence (Deegan 2002). Ultimately, the vegetation and nekton
metrics selected by the MMI assembly algorithm provide a basis for developing causal models of
these systems. It is important to note that the vegetation metrics that were positively correlated
with disturbance in our study of northeastern salt marshes would not indicate degraded
33
34
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
18
conditions in naturally-occurring lower salinity marshes, which are typically dominated by
brackish border species (Supplemental Table 1). These results emphasize the value of clearly
defining the system of interest and the candidate metrics to ensure MMI applicability and
interpretability.
The strong similarity in the vegetation-based MMIs regardless of sample size (MMI V24
and MMI V33, Table 6) suggested that the vegetation communities in the smaller, 24-MSU
dataset adequately represented those of the larger, 33-MSU group of northeastern salt marshes,
and lends confidence to the applicability of nekton metrics included in MMI N24 and MMI
VN24 beyond the 24 MSUs with nekton data. Which MMI provides the “best” integrative index
of northeastern salt marsh condition, however, depends on various considerations. Each MMI
assembly produced MMIs that were maximally sensitive to disturbance within the range of data
available. The exceptionally strong correlation between MMI scores and HDI that emerged when
the assembly algorithm included both vegetation and nekton metrics would argue for MMI
VN24 as the most robust indicator of human disturbance (Figure 3). However, the final MMIs
resulting from each of the other data sets also exhibited strong correlations to HDI (Figure 3) and
might be deemed preferable under certain conditions. In particular, time and resource constraints
often dictate the types of data that can be collected for salt marsh assessment, so an MMI based
on metrics that are the least costly to quantify could be warranted for long-term implementation.
Our analysis suggested that MMIs incorporating strictly vegetation metrics (MMI V33, MMI
V24wt) are highly sensitive to anthropogenic disturbance in these systems. MMI selection may
also be guided by specific management, monitoring, and assessment questions. For example,
estuarine restoration goals often focus partly on re-establishing the hydrologic connectivity
necessary for estuarine and marine fish to access salt marsh habitat (Simenstad and Thom 1996;
35
36
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
19
Dionne et al. 1999; Roman et al. 2002). In such instances, it would be desirable to include nekton
metrics in an MMI used to assess restoration responses. During MMI construction, all candidate
MMIs with Δ AIC < 2 were retained as viable alternatives, and two nekton-only MMIs relied on
one or two metrics alone (MMI N24a, MMI N24e, Table 7); these MMIs provide parsimonious,
highly focused assessment tools. Schoolmaster et al. (2013) described how certain metrics of
intrinsic management value could be weighted to increase the likelihood of being selected by the
MMI assembly algorithm. We did not employ such value-weighting in this study, but that
approach would offer further opportunity to develop MMIs that are highly relevant to defined
management goals.
MMIs are often developed as tools to detect changes in ecosystem condition over time or
to compare the existing condition of multiple systems in space. Our study revealed some
differences in salt marsh MMI trends (Figure 5) and spatial patterns (Figure 4) depending on
MMI formulation, illustrating the potential influence of the metrics used in MMI construction on
interpreting marsh condition. For example, at CACO_3, MMI VN24 scores increased after 2005,
whereas MMI V24 scores decreased during the same period (Figure 5). This MSU, located in
East Harbor Lagoon on outer Cape Cod, MA, underwent tidal restoration between 2001 and
2008, with a concomitant increase in the density of resident fish in marsh pools (Smith et al.
2010). Nekton often respond rapidly when hydrologic connections to salt marsh habitat are
reestablished (e.g. Simenstad and Thom 1996; Burdick et al. 1997), and by incorporating nekton
metrics, MMI VN24 reflected this improved marsh function. Vegetation monitoring from 2008 –
2013 in CACO_2, located in a tidally-unrestricted section of Hatches Harbor on outer Cape Cod,
revealed an abundance of salt marsh grasses with high salinity tolerance, contributing to high
MMI scores among formulations including vegetation metrics; nekton samples, however, were
37
38
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
20
dominated by introduced European green crabs (Carcinus maenas), contributing to a low MMI
N24 score (Figure 4). These differences reinforce the importance of considering monitoring and
assessment goals carefully in identifying, and possibly value-weighting, the candidate metrics
used during MMI construction.
The algorithm-based MMI is, by design, dependent also on the formulation of the HDI.
The PCA approach we used to generate the HDIs summarized the majority of the variation in the
disturbance metrics into a single value for each site. Consequently, the HDIs were related most
to metrics with the highest loadings relative to the first Principal Component score (e.g., Fill /
Fragmentation, Tidal Flushing, and Tidal Restriction). Using different approaches to generating
disturbance indices, such as an additive method (Carullo et al. 2007), might yield different HDIs
and, thus, different MMI metric sets and MMI scores across individual MSUs (Falcone et al.
2010). It is equally important to note that the MMIs characterized marsh condition in the context
of the anthropogenically-influenced disturbance metrics that entered into the HDI. We
incorporated disturbance metrics that quantified a combination of physical, chemical, and land-
use stressors that were easily calculated from existing aerial imagery in a cost-effective and
efficient manner. Other important aspects of marsh degradation, such as the rate of sea-level rise
relative to marsh surface elevation, are not represented in our HDI, and thus the MMI scores we
generated may omit potentially important components of marsh condition.
The MMIs generated here for northeastern salt marshes fit within the multilevel
framework for wetland assessment used widely in the US (Fennessy et al. 2007; Wardrop et al.
2007; Miller et al. 2016). Consisting of a landscape-scale analysis of remotely sensed data (level
1), a rapid field assessment of structural attributes (level 2), and intensive field measurements of
wetland biological, chemical, and physical characteristics (level 3), the multilevel framework
39
40
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
21
provides a comprehensive set of surveys that can be employed singly or in concert, depending on
the goals and purpose of the assessment and the resources available (Wardrop et al. 2007). The
majority of wetland rapid assessment methods include condition indicators focused on vegetation
communities (Fennessy et al. 2007), and many level-2 MMIs are vegetation based (e.g. Wigand
et al. 2011; Miller et al. 2016). The NWCA also relies heavily on a vegetation MMI (VMMI) for
reporting on wetland condition at the national scale (USEPA 2016). Our MMIs generated from
vegetation metrics only (MMI V24, MMI V33) are useful for reporting on salt marsh condition
at the level of resolution of rapid assessments and the NWCA VMMI. An important criterion of
wetland rapid assessments and NWCA monitoring is the amount of field effort demanded:
Fennessy et al. (2007) defined “rapid” assessments as requiring no more than two people a half
day to complete on site, and NWCA field sampling for each site is intended to be accomplished
in one day (USEPA 2011). The quadrat-based vegetation monitoring within each NPS MSU was
generally accomplished by two people in from 0.5 – 1 d per site, and sampling methods requiring
even less field time could be employed as desired to generate the plant species and plant
community metrics that were derived from these data (Table 2). Our nekton-based MMIs were
derived from intensive field measurements of biotic indicators of wetland trophic functions, and
would fall within level-3 wetland monitoring. At both levels, the MMI assembly algorithm
generated MMIs exhibiting the greatest possible correlation with disturbance.
The MMIs produced in this study can be used to assess condition and detect change in
northeastern National Park salt marshes, but we believe the algorithmic approach used here is
broadly applicable to salt marsh assessment. The causal networks linking human disturbance to
the measured biological metrics in salt marsh ecosystems are complex, and interactions among
agents of change, stressors, and responses may introduce unexplained variability and reduce
41
42
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
22
sensitivity of individual metrics to perceived disturbance. The assembly algorithm identified the
optimal metric set, reducing a long list of potential candidate metrics to the vital few using an
objective and transparent selection process. We expect this approach to MMI development
would be useful for salt marsh monitoring and assessment at a variety of scales, from
establishing local baseline conditions and management responses to evaluating regional status
and trends.
Acknowledgements
This work was supported, in part, by funding from USGS National Park Monitoring Program. In
addition, GRG and JGB acknowledge support from the USGS Climate and Land-Use Research
and Development Program. We are grateful to Holly Plaisted for assistance with data
preparation, to Jennifer Olker for analytical guidance, and to the many NPS professional staff
and summer research assistants who helped with vegetation and nekton sampling. We extend
sincere appreciation to Nick Danz and two anonymous reviewers for helpful comments that
improved the manuscript. All data used in this paper are publicly available through USGS at
https://dx.doi.org/10.5066/F7ZP449D. The use of trade, product, or firm names in this
publication is for descriptive purposes only and does not imply endorsement by the US
Government.
Literature Cited
Adamowicz, S.C., and C.T. Roman. 2005. New England salt marsh pools: A quantitative
43
44
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
23
analysis of geomorphic and geographic features. Wetlands 25: 279-288.
Barbier, E.B., S.D. Hacker, C. Kennedy, E.W. Koch, A.C. Stier, and B.R. Silliman. 2011. The
value of estuarine and coastal ecosystem services. Ecological Monographs 81: 169-193.
Bertness, M.D., P.J. Ewanchuk, and B.R. Silliman. 2002. Anthropogenic modification of New
England salt marsh landscapes. Proceedings of The National Academy of Sciences of The
United States of America 99: 1395-1398.
Blocksom, K.A. 2003. A performance comparison of metric scoring methods for a multimetric
index for mid-Atlantic highlands streams. Environmental Management 31: 670–682.
Boix, D., S. Gascón, J. Sala, M. Martinoy, J. Gifre, and X. D. Quintana. 2005. A new index of
water quality assessment in Mediterranean wetlands based on crustacean and insect
assemblages: The case of Catalunya (NE Iberian peninsula). Aquatic Conservation: Marine
and Freshwater Ecosystems 15: 635-651.
Brazner, J.C., N.P. Danz, G.J. Niemi, R.R. Regal, A.S. Trebitz, R.W. Howe, J.M. Hanowski,
L.B. Johnson, J.H. Ciborowski, C.A. Johnston, E.D. Reavie, V.J. Brady, and G.V. Sgro.
2007. Evaluation of geographic, geomorphic and human influences on Great Lakes wetland
indicators: A multi-assemblage approach. Ecological Indicators 7: 610-635.
Burdick, D.M., M. Dionne, R.M. Boumans, and F.T. Short. 1997. Ecological response to tidal
restorations of two northern New England salt marshes. Wetlands Ecology and Management
4: 129-144.
Carullo, M., B.K. Carlisle, and J.P. Smith. 2007. A New England rapid assessment method for
assessing condition of estuarine marshes: A Boston Harbor, Cape Cod and islands pilot
study. Boston: Massachusetts Office of Coastal Zone Management.
45
46
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
24
Chesapeake Bay Program. 2014. Bay Field Guide 2014. Online taxonomic database.
http://www.chesapeakebay.net/fieldguide. Accessed 15 January 2014.
Daiber, F.C. 1986. Conservation of tidal marshes. New York: Van Nostrand Reinhold.
Danz, N.P., G.J. Niemi, R.R. Regal, T. Hollenhorst, L.B. Johnson, J.M. Hanowski, R. Axler,
J.J.H. Ciborowski, T. Hrabik, V.J. Brady, J.R. Kelly, J.A. Morrice, J.C. Brazner, R.W.
Howe, C.A. Johnston, and G.E. Host. 2007. Integrated gradients of anthropogenic stress in
the U.S. Great Lakes basin. Environmental Management 39: 631–647.
Deegan, L.A. 2002. Lessons learned: The effects of nutrient enrichment on the support of
nekton by seagrass and salt marsh ecosystems. Estuaries 25: 727-742.
Dionne, M., F.T. Short, and D.M. Burdick. 1999. Fish utilization of restored, created, and
reference salt-marsh habitat in the Gulf of Maine. American Fisheries Society Symposium 22:
384–404.
Duncan, W.H. 1974. Vascular halophytes of the Atlantic and Gulf Coasts of North Carolina
north of Mexico. In Ecology of halophytes, eds. R.J. Reimold and W.H. Queen, 23-50. New
York: Academic Press.
Edinger, G.J., D.J. Evans, S. Gebauer, T.G. Howard, D.M. Hunt, and A.M. Olivero, eds. 2002.
Ecological communities of New York state, second edition. A revised and expanded edition of
Carol Reschke's Ecological communities of New York state. Albany: New York Natural
Heritage Program, New York State Department of Environmental Conservation.
Eleuterius, L. 1990. Tidal marsh plants. Gretna, LA: Pelican Publishing.
Falcone, J.A., D.M. Carlisle, and L.C. Weber. 2010. Quantifying human disturbance in
watersheds: Variable selection and performance of a GIS-based disturbance index for
predicting the biological condition of perennial streams. Ecological Indicators 10: 264-
47
48
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
25
273.
Fennessy, M.S., A.D. Jacobs, and M.E. Kentula. 2007. An evaluation of rapid methods for
assessing the ecological condition of wetlands. Wetlands 27: 543-560.
FishBase. 2014. A Global Information System on Fishes. Online taxonomic database.
http://www.fishbase.org/home.htm. Accessed January 2014.
Fleming, G.P., and K.D. Patterson. 2013. Natural communities of Virginia: Ecological groups
and community types. Richmond: Virginia Department of Conservation and Recreation
Natural Heritage Technical Report 13-16.
Fry, J.A., G. Xian, S. Jin, J.A. Dewitz,, C.G. Homer, L. Yang, C.A. Barnes, N.D. Herold, and
J.D. Wickham. 2011. Completion of the 2006 National Land Cover Database for the
conterminous United States. Photogrammetric Engineering and Remote Sensing 77: 858-
864.
Greenberg, R., J.E. Maldonado, S. Droege, and M.V. McDonald. 2006. Tidal marshes: A global
perspective on the evolution and conservation of their terrestrial vertebrates. BioScience
56: 675-685.
Haines, A., and T.F. Vining. 1998. Flora of Maine: A manual for identification of native and
naturalized vascular plants. Bar Harbor, ME: Thomas.
Harrison, J.W. 2004. Classification of vegetation communities of Maryland: First iteration.
Annapolis: NatureServe and Maryland Natural Heritage Program Maryland Department of
Natural Resources.
Hering, D., C.K. Feld, O. Moog, and T. Ofenböck. 2006. Cook book for the development of a
multimetric index for biological condition of aquatic ecosystems: Experiences from the
European AQEM and STAR projects and related initiatives. Hydrobiologia 566: 311-324.
49
50
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
26
Hopkinson, C.S., W. Cai, and X. Hu. 2012. Carbon sequestration in wetland dominated coastal
systems— a global sink of rapidly diminishing magnitude. Current Opinion in
Environmental Sustainability 4: 186–194.
Howard, R.J., and I.A. Mendelssohn. 1999. Salinity as a constraint on growth of oligohaline
marsh macrophytes. I. Species variation in stress tolerance. American Journal of Botany 86:
785-794.
Hughes, J.E., L.A. Deegan, J.C. Wyda, M.J. Weaver, and A. Wright. 2002. The effects of
eelgrass habitat loss on estuarine fish communities of southern New England. Estuaries 25:
235-249.
IRL (Indian River Lagoon) Species Inventory. 2014. Smithsonian Indian River Lagoon Species
Inventory. Online taxonomic database. http://www.sms.si.edu/irlspec/index.htm. Accessed
15 January 2014.
James-Pirri, M.J., C.T. Roman, and E.L. Nicosia. 2012. Monitoring nekton in salt marshes: A
protocol for the National Park Service’s long-term monitoring program, Northeast Coastal
and Barrier Network. Fort Collins, CO: National Park Service Natural Resource Report
NPS/NCBN/NRR—2012/579.
James-Pirri, M.J., J.L. Swanson, C.T. Roman, H.S. Ginsberg, and J.F. Heltshe. 2014. Ecological
thresholds for salt marsh nekton and vegetation communities. In Application of threshold
concepts in natural resource decision making, ed. G.R. Guntenspergen, 291-319. New York:
Springer.
Kachigan, S.K. 1982. Multivariate statistical analysis: A conceptual introduction. New York:
Radius Press.
Karr, J.R. 1999. Defining and measuring river health. Freshwater Biology 41: 221-234.
51
52
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
27
Karr, J.R., and E.W. Chu. 2000. Sustaining living rivers. Hydrobiologia 422/423: 1-14.
Kennish, M.J. 2001. Coastal salt marsh systems in the US: A review of anthropogenic impacts.
Journal of Coastal Research 17: 731-748.
Kent, M., and P. Coker. 1994. Vegetation description and analysis – a practical approach.
Chichester, UK: John Wiley and Sons.
Kirwan, M.L., and J.P. Megonigal. 2013. Tidal wetland stability in the face of human impacts
and sea-level rise. Nature 504: 53-60.
Konisky, R.A., D.M. Burdick, M. Dionne, and H.A. Neckles. 2006. A regional assessment of salt
marsh restoration and monitoring in the Gulf of Maine. Restoration Ecology 14: 516-525.
Langman, O.C., J.A. Hale, C.D. Cormack, M.J. Risk, and S.P. Madon. 2012. Developing
multimetric indices for monitoring ecological restoration progress in salt marshes. Marine
Pollution Bulletin 64: 820-835.
Lee, K-S., F.T. Short, and D.M. Burdick. 2004. Development of a nutrient pollution indicator
using the seagrass, Zostera marina, along nutrient gradients in three New England estuaries.
Aquatic Botany 78: 197-216.
Lopez, R.D., and M.S. Fennessy. 2002. Testing the floristic quality assessment index as an
indicator of wetland condition. Ecological Applications 12: 487-497.
Lougheed, V.L., and P. Chow-Fraser. 2002. Development and use of a zooplankton index of
wetland quality in the Laurentian Great Lakes basin. Ecological Applications 12: 474-486.
Mack, J.J. 2004. Integrated wetland assessment program, part 4: Vegetation index of biotic
integrity (VIBI) and tiered aquatic life uses (TALUs) for Ohio wetlands. Columbus: Ohio
Environmental Protection Agency, Wetland Ecology Group, Technical Report WET/2004-4.
53
54
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
28
Magee, D.W., and H.E. Ahles. 1999. Flora of the northeast: A manual of the vascular flora of
New England and adjacent New York. Amherst, MA: Univ. of Massachusetts Press.
MEA (Millennium Ecosystem Assessment). 2005. Ecosystems and human well-being: Current
state and trends - coastal systems. Washington: Island Press.
Mereta, S.T., P. Boets, L. De Meester, and P.L.M. Goethals. 2013. Development of a multimetric
index based on benthic macroinvertebrates for the assessment of natural wetlands in
southwest Ethiopia. Ecological Indicators 29: 510-521.
Miller, K.M., B.R. Mitchell, and B.J. McGill. 2016. Constructing multimetric indices and testing
ability of landscape metrics to assess condition of freshwater wetlands in the northeastern
US. Ecological Indicators 66: 143-152.
Miller, S.J., D.H. Wardrop, W.M. Mahaney, and R.P. Brooks. 2006. A plant-based index of
biological integrity (IBI) for headwater wetlands in central Pennsylvania. Ecological
Indicators 6: 290-312.
Möller, I., M. Kudella, F. Rupprecht, T. Spencer, M. Paul, B.K. van Wesenbeeck, G. Wolters, K.
Jensen, T.J. Bouma, M. Miranda-Lange, and S. Schimmels. 2014. Wave attenuation over
coastal salt marshes under storm surge conditions. Nature Geoscience 7: 727-731.
Moore, G.E., D.M. Burdick, C.R. Peter, and M. Dionne. 2010. Regional assessment of tidal
marsh restoration in New England using the Restoration Performance Index (RPI): Final
report. Gloucester, MA: NOAA Restoration Center.
Neckles, H.A., G.R. Guntenspergen, W.G. Shriver, N.P. Danz, W.A. Wiest, J.L. Nagel, and J.H.
Olker. 2013. Identification of metrics to monitor salt marsh integrity on National Wildlife
Refuges in relation to conservation and management objectives. Laurel, MD: USGS
Patuxent Wildlife Research Center.
55
56
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
29
Neckles, H.A., J.E. Lyons, G.R. Guntenspergen, W.G. Shriver, and S.C. Adamowicz. 2015. Use
of structured decision making to identify monitoring variables and management priorities for
salt marsh ecosystems. Estuaries and Coasts 38: 1215-1232.
Neckles, H.A., M. Dionne, D.M. Burdick, C.T. Roman, R. Buchsbaum, and E. Hutchins. 2002.
A monitoring protocol to assess tidal restoration of salt marshes on local and regional scales.
Restoration Ecology 10: 556-563.
Nestlerode, J.A., V.D. Hansen, A. Teague, and M.C. Harwell. 2014. Application of a three-tier
framework to assess ecological condition of Gulf of Mexico coastal wetlands. Environmental
Monitoring and Assessment 186: 3477–3493.
NYNHP (New York Natural Heritage Program). 2013. Online conservation guide. Available
from: http://acris.nynhp.org/. Accessed 8 February 2014.
Ortega, M., J. Velasco, A. Millán, and C. Guerrero. 2004. An ecological integrity index for
littoral wetlands in agricultural catchments of semiarid Mediterranean regions.
Environmental Management 33: 412-430.
Pendleton, L., D.C. Donato, B. C. Murray, S. Crooks, W.A. Jenkins, S. Sifleet, C. Craft, J.W.
Fourqurean, J. B. Kauffman, N. Marbà, and P. Megonigal. 2012. Estimating global “blue
carbon” emissions from conversion and degradation of vegetated coastal ecosystems. PLoS
One, 7(9), p.e43542.
Raposa, K.B., C.T. Roman, and J.F. Heltshe. 2003. Monitoring nekton as a bioindicator in
shallow estuarine habitats. Environmental Monitoring and Assessment 81: 239-255.
R Core Team. 2016. R: A language and environment for statistical computing. Vienna, Austria:
R Foundation for Statistical Computing, http://www.r-project.org/.
57
58
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
30
Roman, C.T., M.J. James-Pirri, and J.F. Heltshe. 2001. Monitoring salt marsh vegetation.
Wellfleet, MA: National Park Service Technical Report, Long-term Coastal Ecosystem
Monitoring Program at Cape Cod National Seashore.
Roman, C.T., N. Jaworski, F.T. Short, S. Findlay, and R.S. Warren. 2000. Estuaries of the
northeastern United States: Habitat and land use signatures. Estuaries 23: 743-764.
Roman, C.T., K.B. Raposa, S.C. Adamowicz, M.J. James-Pirri, and J.G. Catena. 2002.
Quantifying vegetation and nekton response to tidal restoration of a New England salt marsh.
Restoration Ecology 10: 450-460.
Schafale, M.P. and A.S. Weakley. 1990. Classification of the natural communities of North
Carolina, third approximation. Raleigh, NC: North Carolina Department of Environment,
Health, and Natural Resources, Natural Heritage Program.
Schoolmaster Jr., D.G., J.B. Grace, and E.W. Schweiger. 2012. A general theory of multimetric
indices and their properties. Methods in Ecology and Evolution 3: 773-781.
Schoolmaster Jr., D.R., J.B. Grace, E.W. Schweiger, G.R. Guntenspergen, B.R. Mitchell, K.M.
Miller, and A.M. Little. 2013. An algorithmic and information-theoretic approach to
multimetric index construction. Ecological Indicators 23: 14-23.
Short, F.T., D.M. Burdick, C.A. Short, R.C. Davis, and P.A. Morgan. 2000. Developing success
criteria for restored eelgrass, salt marsh and mud flat habitats. Ecological Engineering 15:
239-252.
Shumway, S.W., and C.R. Banks. 2001. Species distributions in interdunal swale communities:
the effects of soil waterlogging. American Midland Naturalist 145: 137-146.
Silberhorn, G.M. 1999. Common plants of the mid-Atlantic coast: A field guide. Baltimore: Johns
Hopkins University Press.
59
60
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
31
Simenstad, C.A., and R.M. Thom. 1996. Functional equivalency trajectories of the restored Gog-
Le-Hi-Te estuarine wetland. Ecological Applications 6: 38-56.
Simon, T.P. 1998. Modification of an index of biotic integrity and development of reference
condition expectation for dunal, palustrine wetland fish communities along the southern
shore of Lake Michigan. Aquatic Ecosystem Health and Management 1: 49-62.
Smith, S.M., K.C. Medeiros, K. Lee, S.E. Fox, and H. Bayley. 2010. Annual Report on
Estuarine Restoration at East Harbor (Truro, MA), Cape Cod National Seashore, 2009.
Philadelphia: National Park Service Natural Resources Report NPS/NER/NRR—2010.
Sperduto, D.D. 2011. Natural community systems of New Hampshire: 2nd Edition. Concord,
NH: New Hampshire National Heritage Bureau.
Staszak, L.A., and A.R. Armitage. 2013. Evaluating salt marsh restoration success with an index
of ecosystem integrity. Journal of Coastal Research 29: 410-418.
Stevens, S., B. Milstead, M. Albert, and G. Entsminger. 2005. Northeast Coastal and Barrier
Network vital signs monitoring plan. Boston: National Park Service Technical Report
NPS/NER/NRTR—2005/025.
Stoddard, J.L., A.T. Herlihy, D.V. Peck, R.M. Hughes, T.R. Whittier, and E. Tarquinio. 2008. A
process for creating multimetric indices for large-scale aquatic surveys. Journal of the North
American Benthological Society. 27: 878-891.
Swain, P.C., and J.B. Kearsley. 2001. Classification of the natural communities of
Massachusetts, version 1.3. Westborough, MA: Natural Heritage and Endangered Species
Program, Massachusetts Division of Fisheries and Wildlife.
61
62
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
32
Tiner, R.W. 2009. Field guide to tidal wetland plants of the northeastern United States and
neighboring Canada: Vegetation of beaches, tidal flats, rocky shores, marshes, swamps, and
coastal ponds. Amherst, MA: Univ. of Massachusetts Press.
USDA (U.S. Department of Agriculture). 2014. Natural Resources Conservation Service plants
database. http://plants.usda.gov/index.html. Accessed January 2014.
USEPA (U.S. Environmental Protection Agency). 2011. National wetland condition assessment:
Field operations manual. Washington: U.S. Environmental Protection Agency EPA-843-R-
10-001.
USEPA (U.S. Environmental Protection Agency). 2016. National wetland condition assessment:
2011 technical report. Washington: U.S. Environmental Protection Agency EPA-843-R-15-
006.
Wardrop, D.H., M.E. Kentula, S.F. Jensen, D.L. Stevens, Jr., K.C. Hychka, and R.P. Brooks.
2007. Assessment of wetlands in the upper Juniata watershed in Pennsylvania, USA using
the hydrogeomorphic approach. Wetlands 27: 432-445.
Wigand, C., R. Comeleo, R. McKinney, G. Thursby, M. Chintala, and M. Charpentier. 2001.
Outline of a new approach to evaluate ecological integrity of salt marshes. Human and
Ecological Risk Assessment 5: 1541-1554.
Wigand, C., B. Carlisle, J. Smith, M. Carullo, D. Fillis, M. Charpentier, R. McKinney, R.
Johnson, and J. Heltshe. 2011. Development and validation of rapid assessment indices of
condition for coastal tidal wetlands in southern New England, USA. Environmental
Monitoring and Assessment 182: 31-46.
63
64
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
33
Wilcox, D.A., J.E. Meeker, P.L. Hudson, B.J. Armitage, M.G. Black, and D.G. Uzarski. 2002.
Hydrologic variability and the application of index of biotic integrity metrics to wetlands: A
Great Lakes evaluation. Wetlands 22: 588-615.
Wolfe, R.J. 1996. Effects of open marsh water management on selected tidal marsh resources: A
review. Journal of the American Mosquito Control Association 12: 701-712.
Zedler, J.B., and S. Kercher. 2004. Causes and consequences of invasive plants in wetlands:
Opportunities, opportunists, and outcomes. Critical Reviews in Plant Sciences 23: 431–452.
65
66
726
727
728
729
730
731
732
34
Table 1 Location and size (area in hectares) of marsh study units monitored from 2008 – 2013 in NCBN parks, availability of vegetation (V) and nekton (N) data, and number of years sampled
NPS Park State Marsh Study Unit (MSU) Area (ha)
2008-13
Data # YearsAssateague Virginia ASIS_1 5.1 V,N 3Island NS ASIS_2 7.5 V,N 3
ASIS_3 6.7 V,N 3ASIS_4 8.3 V,N 3ASIS_5 8 V,N 3ASIS_6 5.1 V,N 3ASIS_7 7.4 V,N 3ASIS_8 5.4 V,N 3ASIS_9 6.6 V,N 3
Cape Cod NS Massachusetts CACO_1a 55.7 V,N 2CACO_2a 34.6 V,N 2CACO_3a 78.8 V,N 2CACO_4 72.2 V 2CACO_5 26.1 V 2CACO_6 47 V 2CACO_7 25.9 V 2CACO_8 5.6 V 2CACO_9 19.8 V 2CACO_10 112.6 V,N 2CACO_11 167.8 V 2
Fire Island NS New York FIIS_1 5.2 V,N 3FIIS_2 5.2 V,N 3FIIS_3 4.7 V,N 3FIIS_4 6.2 V 2FIIS_5 5 V,N 3FIIS_6 6 V,N 3FIIS_7 5.2 V,N 3FIIS_8 5.4 V,N 3FIIS_9 7.1 V,N 3
Gateway NRA, Sandy Hook Unit New Jersey GSH_1a 7 V,N 2
GSH_2 15.1 V 2GSH_3 6.5 V,N 2
Sagamore Hill NHS New York SAHI_1a 2.7 V,N 3a denotes units that were also sampled from 2003-7 as part of a pilot monitoring program
67
68
733
35
Table 2 Candidate vegetation metrics from NCBN parks for incorporation into salt marsh MMIs
Vegetation CommunitiesPercent cover of low marsh speciesPercent cover of high marsh speciesPercent cover of brackish border speciesPercent cover of salt marsh border speciesPercent cover of panne, pool and creek speciesPercent cover of native speciesPercent cover of high salinity-tolerant speciesSpecies richness (total no. spp. within each sample averaged over MSU)Native species richness (total no. native spp. within each sample averaged over MSU)Native species richness (total no. native spp. within MSU)Introduced species richness (total no. introduced spp. within each sample averaged over MSU)a
Introduced species richness (total no. introduced spp. within MSU)a
Percent native species richness (no. native spp. in MSU /total no. spp. in MSU)Percent introduced species richness (no. introduced spp. in MSU /total no. spp. in MSU)b
High salinity-tolerant species richness (no. high salinity-tolerant spp. within MSU)High salinity-tolerant species richness (no. high salinity-tolerant spp.in sample /total no. spp. in sample)Percent high salinity-tolerant species (no. high salinity-tolerant spp. in MSU /total no. spp. in MSU)Percent high salinity-tolerant species (no. high salinity-tolerant spp. in sample /total no. spp. in sample)
Vegetation speciesPercent frequency Spartina alternifloraPercent cover Spartina patensPercent frequency S. patensPercent cover Distichlis spicataPercent frequency D. spicataPercent cover Salicornia spp.Percent frequency Salicornia spp.Percent frequency Iva frutescensPercent cover Scirpus and Schoenoplectus sppb
Percent frequency Scirpus and Schoenoplectus sppb
Percent cover Phragmites australisPercent frequency P. australis
a Metric used only in the dataset consisting of vegetation data alone from 33 MSUs (V33)b Metric used only in the dataset from the 24-MSU subset including both vegetation and nekton data (V24, VN24)
69
70
734735
736737738739
36
Table 3 Candidate nekton metrics from NCBN parks for incorporation into salt marsh MMIs
Nekton communitiesRelative abundance of resident fish speciesRelative abundance of resident shrimp speciesRelative abundance transient fish speciesTotal density (no. m-2)Fish density (no. fish m-2)Decapod density (no. decapods m-2)a
Species richness (total no. spp. within each sample averaged over MSU)Native species richness (total no. native spp. within MSU) SiteIntroduced species richness (total no. introduced spp. within MSU)Resident fish species richness (total no. resident fish spp. within MSU)Resident fish species richness (total no. resident fish spp. within each sample averaged over MSU)Percent native species richness (no. native spp. in MSU /total no. spp. in MSU)Percent native species richness (no. native spp. in sample /total no. spp. in sample)Percent resident fish species richness (no. resident fish spp.in MSU /total no. spp. in MSU)
Nekton speciesRelative abundance of Fundulus spp.Relative abundance of Cyprinodon variegatusRelative abundance of Menidia spp.Length of Fundulus heteroclitus
aDecapod = Crabs and Shrimp
71
72
740741
742743
37
Table 4 Land use metrics used to calculate the Human Disturbance IndexMetric Classes/Rank Definition
Ditch density no (0), low (2), moderate (4), severe (6)
Extent of ditching and draining of the marsh unit
Tidal restriction no (0), yes (6) Potential restriction of normal flow and tidal range by features such as under-sized culverts or bridges, causeways, dikes, etc.
Tidal flushing well flushed (0), moderately flushed (3),poorly flushed (6)
Degree of hydrologic connection with marine waters
Fill/fragmentation no (0), low (2), moderate (4), severe (6)
Extent of salt marsh areas filled or fragmented (i.e. were once whole systems).
Relativized %Disturbed_150m
continuous converted to ranked variable
% disturbed landa in 150-m buffer * (area of buffer/area of MSU)
Relativized %Disturbed_1km
continuous converted to ranked variable
% disturbed landa in 1-km buffer * (area of buffer/area of MSU)
Point source discharges of pollutants
no (0), low (2), moderate (4), severe (6)
Outfalls, drains emptying into marsh
aAggregated NLCD classes: Developed Open Space, Developed Low Intensity, Developed Medium Intensity, Developed High Intensity, Barren Land and Agricultural
73
74
744
745746747748
38
Table 5 Principle Component loadings (PC1 – PC7) for each disturbance metric incorporated in the HDIs. Eigenvalue size for each Principle Component is also shown
HDI33
Disturbance Metrics PC1 PC2 PC3 PC4 PC5 PC6 PC7Ditch Density 0.25 -0.50 -0.29 0.67 0.34 0.19 -0.08Fill / Fragmentation -0.52 -0.18 -0.08 0.00 -0.39 0.60 -0.41Tidal Flushing -0.50 -0.07 -0.03 -0.20 0.69 -0.28 -0.39Tidal Restriction -0.43 -0.43 -0.32 -0.15 -0.07 -0.11 0.70Point Source Discharges of Pollutants -0.40 0.36 -0.12 0.64 -0.28 -0.46 -0.04Relativized % Disturbed_150m -0.26 0.08 0.73 0.30 0.26 0.31 0.36Relativized % Disturbed_1km 0.05 -0.62 0.50 0.00 -0.32 -0.45 -0.23
Eigenvalue Size 2.68 1.46 1.31 0.67 0.47 0.28 0.13
HDI24
Disturbance Metrics PC1 PC2 PC3 PC4 PC5 PC6 PC7Ditch Density -0.13 0.64 0.14 0.6 -0.43 0.13 0.03Fill / Fragmentation 0.51 0.22 0.12 -0.27 0.07 0.78 -0.01Tidal Flushing 0.49 -0.29 0.11 -0.03 -0.6 -0.21 0.51Tidal Restriction 0.45 0.35 0.34 -0.21 0.03 -0.52 -0.5Point Source Discharges of Pollutants 0.36 -0.23 0.27 0.65 0.55 -0.03 0.14Relativized % Disturbed_150m 0.32 -0.23 -0.64 0.31 -0.23 0.08 -0.53Relativized % Disturbed_1km 0.23 0.48 -0.6 -0.07 0.31 -0.25 0.44
Eigenvalue Size 2.90 1.40 1.23 0.80 0.38 0.18 0.10
75
76
749
39
Table 6 Metrics included in the MMIs and correlations of individual metrics with HDI scores for datasets including vegetation metrics only. Single-letter subscripts refer to candidate MMIs with ΔAIC < 2; wt refers to the model-weighted MMI. Dashes indicate metrics that were not included in the MMI for that dataset
Correlation with HDI
Metrics included in MMIMMI V33
MMI V24a
MMI V24b
MMI V24wt
Frequency of Distichlis spicata -0.52 -0.59 -0.59 -0.59Frequency of Phragmites australis 0.65 0.67 - 0.67Frequency of Spartina patens -0.56 - - -Percent cover of brackish border spp. - 0.64 - 0.64Percent cover of Phragmites australis - - 0.73 0.73Percent cover of salt marsh border spp. 0.12 - 0.12 0.12Percent of high salinity spp. at sample level -0.80 - -0.75 -0.75
Correlation of MMI with HDI -0.93 -0.89 -0.90 -0.91
77
78
750751752753
754
755756
40
Table 7 Metrics included in the MMIs and correlations of individual metrics with HDI scores for the dataset including nekton metrics only. Single-letter subscripts refer to candidate MMIs with ΔAIC < 2; wt refers to the model-weighted MMI. Dashes indicate metrics that were not included in the MMI
Correlation with HDI
Metrics included in MMIMMI N24a
MMI N24b
MMI N24c
MMI N24d
MMI N24e
MMI N24wt
Relative abundance of resident shrimp - 0.49 0.49 - - 0.49Relative abundance of resident fish -0.79 -0.79 - - -0.79 -0.79Total density - - - 0.30 - 0.30Introduced spp. at site level - - 0.59 - - 0.59Resident fish spp. at site level - - -0.58 - -0.58 -0.58Resident fish spp. at sample level - - - -0.64 - -0.64Percent of native spp. at sample level - -0.44 - -0.44 - -0.44
Correlation of MMI with HDI -0.80 -0.84 -0.84 -0.83 -0.80 -0.85
79
80
757758759760
761
762763
41
Figure Captions
Fig. 1 Location of parks sampled within the Northeast Coastal and Barrier Network. Black
triangles = MSUs with vegetation and nekton data; white circles = MSUs with vegetation data
only
Fig. 2 MSU scores for the 24- (HDI 24) and 33-MSU (HDI 33) Human Disturbance Indices
(HDI). MSUs are listed along the horizontal axis in latitudinal order from north to south. CACO
= Cape Cod National Seashore, SAHI = Sagamore Hill National Historic Site, FIIS = Fire Island
National Seashore, GSH = Gateway National Recreation Area – Sandy Hook Unit, and ASIS =
Assateague Island National Seashore
Fig. 3 MMI scores as a function of Human Disturbance Index for the (A) MMIVN(24), (B) model-
weighted MMIV(24), (C) model-weighted MMIN(24), and (D) MMIV(33) datasets. Correlations (r-
values) are also shown
Fig. 4 MMI scores by MSU for the V33, V24, N24, and VN24 datasets. MSUs are listed along
the horizontal axis in latitudinal order from north to south. CACO = Cape Cod National
Seashore, SAHI = Sagamore Hill National Historic Site, FIIS = Fire Island National Seashore,
GSH = Gateway National Recreation Area – Sandy Hook Unit, and ASIS = Assateague Island
National Seashore
Fig. 5 Change in MMI scores over time at 5 MSUs using the (A) MMI VN24 and (B) MMI V24
formulations. CACO = Cape Cod National Seashore, SAHI = Sagamore Hill National Historic
Site, and GSH = Gateway National Recreation Area – Sandy Hook Unit
81
82
764
765
766
767
768