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 2004–2009 (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 (12–13 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
pigments
(mg/m²)

Index of
P-limitation

NH4
(μg /L)

NO3
(μg/L)

SRP
(μg/L)

Water
temperature
(°C)

PAR
(μmol/s
·m²)

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:716–730.


[Back to C006-086]