Ecological Archives C006-086-A1
Barbara L. Peckarsky, Angus R. McIntosh, Maruxa Ālvarez, and Jennifer M. Moslemi. 2015. Disturbance legacies and nutrient limitation influence interactions between grazers and algae in high elevation streams. Ecosphere 6:241. http://dx.doi.org/10.1890/es15-00236.1
Appendix A. Detailed description of streams where surveys and experiments were conducted.
Stream site selection
Twenty 50 m sampling sites (stream reaches) were selected in the upper East River drainage basin north of RMBL to provide a wide range of conditions for comparisons of stream size, flow-related disturbance, predation regimes, invertebrates and aquatic plant assemblages. Sampling sites were separated spatially to ensure they could be considered independent replicates (see map in Peckarsky et al. 2014).
Biological variables surveyed in 20 streams
During July and August 2006 we compared the composition of benthic invertebrates, aquatic plant assemblages, and predatory fish densities among streams with different disturbance regimes, and tested for associations among trophic levels. We estimated the relative proportions of diatoms, moss, and filamentous forms of algae at all 20 sites. We took five randomly located quadrat samples per site to estimate the percent of the stream bottom covered with diatoms, filamentous forms of algae, and moss when benthic algal biomass was at its peak (Table A1). We also disturbed the substrate inside a bottomless plastic bucket (diameter 30 cm) and collected organic matter in 100 mL cups to estimate the quantity of benthic organic matter (BOM) ash free dry mass (AFDM) and chlorophyll (chl a) to compare the standing crop of basal resources used by primary consumers in all streams (Table A3). Blanks in the table indicate streams with no data. The 10 streams that were used for the mesocosm experiment in 2007 are indicated in bold font, and the five streams used for in-stream channel experiment in 2010 are indicated in bold caps (Table A3).
We estimated macroinvertebrate density and species composition of grazers, shredders, and predators (the predominant functional feeding groups) using a collapsible modified box sampler with a stainless steel frame that covers an area of substrate = 0.104 m², and can be fitted with a standard Wilco ® drift net with mesh size 0.364 μm. This sampler was invented by David Hoffman at the University of Wisconsin. Plans are available from B. L. Peckarsky ([email protected]). We took five benthic samples totaling 0.52 m² of stream bed one time during July 2006 to estimate density and species composition of each functional group, the density of Chironomidae and Oligochaeta, which could not be identified beyond family and were thereby not classified into functional feeding groups, and total invertebrate density (Table A2).
We also estimated the density of predatory fish (mostly brook trout) using a backpacking electroshocker in all 20 streams in August 2006. Thirteen streams were completely fishless (Fig. 1A), and the other seven had a range of densities of brook trout. We have presented the mean density of brook trout in those seven streams based on annual fish surveys that have been conducted since 1996 (Table A3).
Multivariate axis of disturbance
We directly characterized variation among flow regimes by continuously monitoring all 20 sites using Trutrack stage height data loggers from 20042009 (Peckarsky et al. 2014). We used stage data from loggers, calibrated by manual estimates of discharge over as wide a range of stages as possible, to compare the hydrologic regimes among sites and combined metrics related to the hydrologic regime with empirical estimates of substrate instability and a qualitative assessment of stream bed stability to classify all 20 sites according to variation in hydrology and geomorphology (low numbers are stable and high numbers are most disturbed) using a principal components analysis (PC axis 1; Table A3) (Peckarsky et al. 2014).
Other abiotic variables
In July 2006 we also measured other abiotic variables in all 20 streams (Table A3) that could explain not only variation in disturbance (see Peckarsky et al. 2014 for additional variables measured), but also variation in the abundance of algae and invertebrate grazers. Drainage basin areas were estimated from digital elevation maps using Arc-GIS software. Maximum discharge (Qmax) in each stream over the 5-year period of the study was obtained from calibrated estimates from stage-height data loggers. We also collected water samples during two dates during July 2006 to estimate soluble reactive phosphorus-P (SRP), ammonium-N (NH4), and nitrate-N (NO3) following standard protocols described in the methods. Mean values of those two dates are displayed in Table A3.
Biological variables in experimental streams (2007 and 2010)
Before and at the conclusion of the mesocosm experiment (1213 July and 24 July 2007) we collected four rocks from each stream used for the experiment and estimated the biomass of chl-a using protocols described in the Methods section. Values represent the means of the two sampling dates (Table A4). We estimated the densities of macroinvertebrate grazers using the collapsible modified box sampler. We took five benthic samples totaling 0.52 m² of stream bed in July 2007 at the end of the mesocosm experiment, and five benthic samples at the end of the stream channel experiment in July 2010. Densities of grazing invertebrates collected in 2007 and 2010 are displayed in Figs. 1B, 2B, and 3B.
Abiotic variables in experimental streams
In each year when experiments were conducted (2007 and 2010) we measured abiotic variables that could be related to the accrual of algae and the strength of grazer-algal interactions, which were the response variables for the experiments. Disturbance regimes, drainage basin areas and maximum discharge in each stream over the 5-year period of the study are displayed in Table A3 (see stream names in bold font). Water temperatures were measured year-round using the Trutrack stage-height data loggers and Hobo data loggers during the summer months. We took one-time measurements of light levels (photosynthetically active radiation - PAR) during the mesocosm experiment in July 2007. We also collected water samples during the mesocosm experiment in 2007 to estimate soluble reactive phosphorus-P (SRP), ammonium-N (NH4), and nitrate-N (NO3) following standard protocols described in the methods. Results from 2006 experiments deploying nutrient diffusing substrates (NDS) were used to calculate an index of phosphorus limitation in the 10 streams used for the mesocosm experiment (Moslemi 2010). The P limitation response index represents the effect size of increases of algal biomass when P was added to the NDS (higher numbers indicate greater P limitation).
Conditions in the mesocosms
We used the densities and species composition of grazing mayflies and caddisflies sampled during 2006 (Fig. 1A) to estimate the relative proportions of each species of invertebrate grazer to stock in mesocosms to establish assemblages that reflected the legacy of disturbance in each stream. Analysis of the 2007 invertebrate samples (Fig. 1B) confirmed that those assemblages remained consistent from year to year. At the end of the experiment we counted and weighed the grazers recovered from each mesocosm to estimate final densities and biomasses of invertebrates grazing algae under ambient and enriched nutrient conditions (Tables A5 and A6), and used those values to calculate the grazer-algal interaction strengths.
Table A1. Proportions of aquatic plants in quadrats and densities of invertebrates (number /m²) in benthic samples from 20 streams surveyed in July 2006, sorted from least to most disturbed.
Stream |
Diatoms |
Moss |
Filaments |
Grazers |
Shredders |
Predators |
Chironomids |
Oligochaeta |
Total Invertebrates |
Mossy |
0.12 |
0.53 |
0 |
231 |
281 |
177 |
3085 |
444 |
2924 |
Talus |
0.49 |
0.13 |
0 |
210 |
204 |
1075 |
3110 |
52 |
4007 |
B3 |
0.50 |
0.45 |
0 |
444 |
154 |
681 |
1558 |
23 |
2001 |
Whiterock |
0.50 |
0.27 |
0.03 |
488 |
290 |
392 |
3975 |
21 |
2794 |
B9 |
0.86 |
0.06 |
0 |
190 |
171 |
131 |
348 |
75 |
598 |
Baldy |
0.63 |
0.32 |
0 |
456 |
323 |
202 |
498 |
33 |
902 |
Pika |
0.68 |
0.18 |
0 |
273 |
738 |
315 |
938 |
135 |
1435 |
B2 |
0.96 |
0 |
0.01 |
319 |
108 |
412 |
1102 |
423 |
1257 |
Mosca |
0.95 |
0 |
0.01 |
490 |
550 |
288 |
723 |
114 |
1404 |
Lower Benthette |
0.90 |
0.01 |
0.07 |
460 |
452 |
213 |
923 |
100 |
1432 |
Avery |
0.93 |
0 |
0 |
1633 |
277 |
262 |
492 |
98 |
1502 |
Marmot |
0.93 |
0 |
0.01 |
712 |
487 |
313 |
913 |
129 |
1517 |
Sylvanite |
0.92 |
0 |
0.02 |
1438 |
283 |
165 |
744 |
4 |
1491 |
Lower Bradley |
0.92 |
0.04 |
0 |
1058 |
412 |
338 |
448 |
185 |
1386 |
Upper East |
0.95 |
0 |
0.05 |
937 |
367 |
387 |
256 |
14 |
1071 |
Rustlers |
1.00 |
0 |
0 |
1537 |
429 |
306 |
363 |
0 |
1470 |
Gothic |
0.95 |
0 |
0 |
190 |
125 |
190 |
271 |
35 |
625 |
Quigley |
1.00 |
0 |
0 |
162 |
42 |
73 |
60 |
14 |
208 |
East River - Gothic |
0.96 |
0 |
0.04 |
1162 |
144 |
81 |
519 |
14 |
1085 |
Copper |
0.96 |
0 |
0.03 |
1519 |
415 |
246 |
885 |
89 |
1722 |
Table A2. Dominant taxa in each functional feeding group of invertebrates in 20 streams surveyed in July 2006, sorted from least to most disturbed.
Stream |
Grazers |
Shredders |
Predators |
Mossy |
Caddisflies (Neothrema) |
Zapada |
Sweltza, Rhyacophila, Flatworms |
Talus |
Caddisflies (Allomyia) |
Zapada |
Flatworms |
B3 |
Caddisflies (Allomyia, Neothrema) |
Zapada |
Sweltza, Kogotus, Rhyacophila, Flatworms |
Whiterock |
Caddisflies (Allomyia) |
Zapada |
Rhyacophila, Flatworms |
B9 |
Mayflies (Cinygmula), caddisflies (Neothrema) |
Zapada |
Sweltza, Kogotus |
Baldy |
Mayflies (Baetis) |
Zapada |
Sweltza, Kogotus |
Pika |
Caddisflies (Allomyia) |
Zapada |
Dichronota, Flatworms |
B2 |
Mayflies (Baetis, Cinygmula), caddisflies (Allomyia) |
Zapada |
Kogotus, Rhyacophila, Flatworms |
Mosca |
Mayflies (Cinygmula) |
Zapada |
Sweltza, Megarcys, Rhyacophila, Flatworms |
Lower Benthette |
Mayflies (Baetis), Caddisflies (Neothrema) |
Zapada |
Sweltza, Megarcys, Rhyacophila, Flatworms |
Avery |
Mayflies (Baetis) |
Zapada |
Hesperoperla, Kogotus, Flatworms |
Marmot |
Mayflies (Baetis, Cinygmula) |
Zapada |
Kogotus, Rhyacophila, Flatworms |
Sylvanite |
Mayflies (Baetis) |
Zapada |
Megarcys, Rhyacophila |
Lower Bradley |
Mayflies (Baetis) |
Zapada |
Sweltza, Megarcys, Rhyacophila, Flatworms |
Upper East |
Mayflies (Baetis, Cinygmula, Epeorus, Rhithrogena) |
Zapada |
Sweltza, Kogotus, Rhyacophila, Flatworms |
Rustlers |
Mayflies (Baetis, Cinygmula, Epeorus, Rhithrogena) |
Zapada |
Kogotus, Rhyacophila, Flatworms |
Gothic |
Mayflies (Cinygmula, Epeorus, Rhithrogena) |
Zapada |
Rhyacophila, Flatworms |
Quigley |
Mayflies (Baetis, Cinygmula, Epeorus, Rhithrogena) |
Zapada |
Megarcys, Rhyacophila |
East River -Gothic |
Mayflies (Baetis) |
Zapada |
Sweltza, Rhyacophila |
Copper |
Mayflies (Cinygmula, Epeorus, Rhithrogena) |
Zapada |
Sweltza, Megarcys, Rhyacophila |
Table A3. Conditions in the 20 streams surveyed in 2006, sorted from least to most disturbed (PC Axis 1; Peckarsky et al. 2014). Abbreviations are described in the text of Appendix A.
Stream |
Fish density (no/m²) |
BOM AFDM (g/m²) |
BOM chl a (mg/m²) |
NH4 (μg/L) |
NO3 (μg/L) |
SRP (μg/L) |
PC Axis 1 |
Drainage basin area (km²) |
Q max (m3/s) |
Mossy |
0 |
3.10 |
117.67 |
0.44 |
402 |
1.00 |
-2.54 |
0.12 |
0.02 |
Talus |
0 |
0.27 |
|
0.36 |
478 |
7.30 |
-2.33 |
0.04 |
0.03 |
B3 |
0 |
0.49 |
12.89 |
0.25 |
387 |
0.85 |
-1.75 |
0.001 |
0.04 |
Whiterock |
0 |
0.74 |
56.26 |
0.12 |
459 |
1.45 |
-1.68 |
1.08 |
0.07 |
B9 |
0 |
0.23 |
3.95 |
0.85 |
50 |
0.35 |
-1.21 |
0.25 |
0.04 |
Baldy |
0 |
0.11 |
|
0 |
288 |
3.30 |
-1.18 |
0.10 |
0.05 |
Pika |
0 |
0.49 |
5.41 |
0.44 |
495 |
2.45 |
-1.15 |
0.07 |
0.19 |
B2 |
0 |
0.27 |
4.46 |
0.93 |
377 |
1.05 |
-0.75 |
0.0003 |
0.53 |
Mosca |
0 |
0.09 |
0.55 |
0.97 |
153 |
2.60 |
-0.05 |
0.18 |
0.06 |
Lower Benthette |
0 |
0.16 |
7.31 |
1.82 |
270 |
1.90 |
0.18 |
0.27 |
0.50 |
Avery |
0.41 |
0.18 |
6.09 |
0.97 |
224 |
1.15 |
0.35 |
0.71 |
0.08 |
Marmot |
0 |
0.34 |
3.76 |
0.73 |
232 |
1.00 |
0.65 |
0.28 |
0.16 |
Sylvanite |
0 |
0.07 |
2.41 |
0.23 |
227 |
1.90 |
0.71 |
1.78 |
0.24 |
LOWER BRADLEY |
0.09 |
0.07 |
1.35 |
0.58 |
322 |
1.75 |
0.82 |
3.73 |
0.58 |
UPPER EAST |
0.13 |
0.05 |
2.33 |
1.16 |
459 |
2.00 |
1.12 |
7.64 |
1.60 |
RUSTLERS |
0.11 |
0.23 |
3.83 |
1.82 |
272 |
0.95 |
1.20 |
15.0 |
2.64 |
Gothic |
0 |
0.05 |
0.93 |
0.40 |
20 |
2.25 |
1.58 |
1.26 |
0.20 |
QUIGLEY |
0.01 |
0.09 |
0.97 |
1.27 |
474 |
3.50 |
1.79 |
1.80 |
1.16 |
East River -Gothic |
0.25 |
0.18 |
18.46 |
1.52 |
187 |
0.80 |
1.97 |
44.9 |
13.37 |
COPPER |
0.02 |
0.11 |
5.82 |
1.14 |
226 |
1.35 |
2.32 |
23.2 |
5.51 |
Table A4. Conditions of 10 streams used for the mesocosm experiment in 2007, sorted from least to most disturbed. Abbreviations are described in the text of Appendix A.
Stream |
Total algae |
Index of |
NH4 |
NO3 |
SRP |
Water |
PAR |
B3 |
13.80 |
0.19 |
0.45 |
157 |
1.00 |
5.5 |
414 |
B9 |
15.35 |
0.83 |
1.40 |
41.2 |
0.35 |
7.0 |
423 |
B2 |
16.41 |
0.61 |
1.45 |
160 |
0.25 |
6.1 |
215 |
Lower Benthette |
23.82 |
0.47 |
2.10 |
95.7 |
1.25 |
10.5 |
73 |
Marmot |
26.00 |
0.57 |
1.80 |
86.8 |
1.70 |
6.9 |
231 |
Lower Bradley |
42.73 |
0.07 |
1.60 |
118 |
2.05 |
7.2 |
827 |
Upper East |
21.29 |
0.39 |
3.25 |
60.7 |
1.70 |
10.4 |
1273 |
Quigley |
13.37 |
0.88 |
3.75 |
131 |
2.65 |
9.4 |
52 |
East River - Gothic |
31.18 |
0.39 |
2.90 |
66.2 |
2.20 |
11.6 |
564 |
Copper |
17.51 |
0.19 |
2.75 |
81.2 |
1.20 |
8.6 |
271 |
Table A5. Conditions in mesocosms with ambient (limited) nutrients adjacent to 10 streams sorted from least to most disturbed.
Stream |
Density Mayflies (no/m²) |
Density Caddisflies (no/m²) |
Density total grazers (no/m²) |
Weight Mayflies (g) |
Weight Caddisflies (g) |
Weight total grazers (g) |
Interaction strength (per grazer biomass) × 1000 |
B3 |
236 |
270 |
506 |
8.35 |
7.37 |
15.73 |
-6.24 |
B9 |
185 |
84 |
270 |
10.76 |
1.42 |
12.18 |
-5.70 |
B2 |
202 |
135 |
337 |
12.41 |
4.72 |
17.12 |
-1.29 |
Lower Benthette |
354 |
185 |
540 |
18.68 |
5.50 |
24.18 |
-1.18 |
Marmot |
573 |
0 |
573 |
39.86 |
0.00 |
39.86 |
-0.26 |
Lower Bradley |
118 |
0 |
118 |
8.84 |
0.00 |
8.84 |
-2.84 |
Upper East |
759 |
0 |
759 |
54.18 |
0.00 |
54.18 |
-0.44 |
Quigley |
84 |
0 |
84 |
6.70 |
0.00 |
6.70 |
0.68 |
East River - Gothic |
725 |
0 |
725 |
59.28 |
0.00 |
59.28 |
-0.65 |
Copper |
1231 |
0 |
1231 |
89.76 |
0.00 |
89.76 |
-0.57 |
Table A6. Conditions in mesocosms with enriched nutrients adjacent to 10 streams sorted from least to most disturbed.
Stream |
Density Mayflies (no/m²) |
Density Caddisflies (no/m²) |
Density total grazers (no/m²) |
Weight Mayflies (g) |
Weight Caddisflies (g) |
Weight total grazers (g) |
Interaction strength (per grazer biomass) × 1000 |
B3 |
202 |
270 |
472 |
8.62 |
7.28 |
15.90 |
-0.98 |
B9 |
185 |
84 |
270 |
11.96 |
1.66 |
13.62 |
0.53 |
B2 |
320 |
185 |
506 |
17.02 |
8.74 |
25.75 |
-0.67 |
Lower Benthette |
422 |
51 |
472 |
22.23 |
1.13 |
23.36 |
-0.01 |
Marmot |
590 |
0 |
590 |
46.21 |
0.00 |
46.21 |
-0.43 |
Lower Bradley |
691 |
0 |
691 |
49.66 |
0.00 |
49.66 |
-0.50 |
Upper East |
860 |
0 |
860 |
67.93 |
0.00 |
67.93 |
-0.79 |
Quigley |
101 |
0 |
101 |
8.28 |
0.00 |
8.28 |
0.42 |
East River - Gothic |
843 |
0 |
843 |
62.67 |
0.00 |
62.67 |
-0.21 |
Copper |
894 |
0 |
894 |
61.37 |
0.00 |
61.37 |
-0.93 |
Literature Cited
Peckarsky, B. L, A. R. McIntosh, S. C. Horn, A. Wilcox, K. McHugh, D. Booker, W. S. Brown, and M. Alvarez. 2014. Characterizing disturbance regimes of mountain streams. Freshwater Science 33:716730.