A Model for Methylmercury Uptake and Trophic ... - Harvard University

Dec 11, 2017 - zooplankton reflected lengthening of the food web rather than size based effects on phytoplankton ... ide (SI Table S2). Dietary assimi...
0 downloads 0 Views 752KB Size
Subscriber access provided by READING UNIV

Article

A Model for Methylmercury Uptake and Trophic Transfer by Marine Plankton Amina Traore Schartup, Asif Qureshi, Clifton Dassuncao, Colin Pike Thackray, Gareth Harding, and Elsie M. Sunderland Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b03821 • Publication Date (Web): 11 Dec 2017 Downloaded from http://pubs.acs.org on December 14, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25

Environmental Science & Technology

A Model for Methylmercury Uptake and Trophic Transfer by Marine Plankton Amina T. Schartup1,2,*, Asif Qureshi2,3, Clifton Dassuncao1,2, Colin P. Thackray1, Gareth Harding4, and Elsie M. Sunderland1,2 1

Harvard John A. Paulson School of Engineering & Applied Sciences, Harvard University, Cambridge, MA, USA 02138

2

Department of Environmental Health, Harvard T. H. Chan School of Public Health, Harvard University, Boston, MA, USA 02214

3

Department of Civil Engineering, IIT Hyderabad, Kandi, Sangareddy, TS 502285, India

4

Bedford Institute of Oceanography, Dartmouth, NS, Canada

*

Corresponding author: 29 Oxford Street, Cambridge, MA 02138; [email protected]

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 26

26

Abstract

27

Methylmercury (MeHg) concentrations can increase by up to 100,000 times between seawater

28

and marine phytoplankton but levels vary across sites. To better understand how ecosystem

29

properties affect variability in planktonic MeHg concentrations, we develop a model for MeHg

30

uptake and trophic transfer at the base of marine food webs. The model successfully reproduces

31

measured concentrations in phytoplankton and zooplankton across diverse sites from the

32

Northwest Atlantic Ocean. Highest MeHg concentrations in phytoplankton are simulated under

33

low dissolved organic carbon (DOC) and ultra-oligotrophic conditions typical of open ocean

34

regions. This occurs because large organic complexes bound to MeHg inhibit cellular uptake and

35

cell surface area to volume ratios are greatest under low productivity conditions. Modeled

36

bioaccumulation factors for phytoplankton (102.4-105.9) are more variable than those for

37

zooplankton (104.6-106.2) across ranges in DOC (40-500 µM) and productivities (ultra-

38

oligotrophic to hyper-eutrophic) typically found in marine ecosystems. Zooplankton growth

39

dilutes their MeHg body burden but they also consume greater quantities of MeHg enriched prey

40

at larger sizes. These competing processes lead to lower variability in MeHg concentrations in

41

zooplankton compared to phytoplankton. Even under hyper-eutrophic conditions, modeled

42

growth dilution in marine zooplankton is insufficient to lower their MeHg concentrations,

43

contrasting findings from freshwater ecosystems.

44

ACS Paragon Plus Environment

2

Page 3 of 26

Environmental Science & Technology

45 46 47

Introduction

48

wildlife and is the only form of mercury (Hg) that biomagnifies in food webs.1,2 Marine fish and

49

shellfish account for more than 90% of the population-wide MeHg intake in the United States.3

50

Concentrations of MeHg in phytoplankton can be up to 100,000 times higher than seawater.

51

However, few studies have examined how ecosystem properties affect differences in MeHg

52

uptake and trophic transfer at the base of marine food webs.4,5 Here, we develop a new model for

53

MeHg bioaccumulation in marine phytoplankton and zooplankton.

54

Methylmercury (MeHg) is a potent neurotoxin and oxidative stressor in humans and

Marine phytoplankton are thought to accumulate MeHg mainly by passive uptake from

55

seawater (diffusion) across the cell membrane.6,7 Surface area to volume ratios of phytoplankton

56

are therefore critical for anticipating total uptake.8 When ecosystems are less productive, smaller

57

phytoplankton with greater surface area to volume ratios are more abundant to maximize nutrient

58

uptake efficiency, which effectively increases MeHg uptake.9,5 Conversely, higher nutrient status

59

can lead to a greater abundance of larger cells, reducing nutrient uptake at the base of the food

60

web10 and potentially also MeHg uptake.

61

Prior work suggests MeHg trophic transfer from phytoplankton to zooplankton is mainly

62

driven by dietary intake with smaller direct uptake from seawater (10-20%).11 Loss pathways for

63

MeHg in zooplankton include growth dilution and temperature dependent elimination in fecal

64

pellets.11,12 Herbivorous (small) zooplankton grow most rapidly at higher temperatures and when

65

sufficient nutrients are available.13 Omnivorous (large) zooplankton are opportunistic feeders and

66

preferentially consume the largest prey available (other smaller zooplankton) but will also

67

consume phytoplankton, depending on availability.14 Biomass dilution has been established as a

68

major factor controlling MeHg concentrations in zooplankton from eutrophic freshwater

ACS Paragon Plus Environment

3

Environmental Science & Technology

69

systems.6,15 Similar processes for marine ecosystems have been hypothesized but not directly

70

observed.4

71

Page 4 of 26

High concentrations of terrestrial dissolved organic matter may inhibit phytoplankton

72

MeHg uptake by forming large complexes with MeHg in seawater that reduce transport across

73

the cell membrane.16–19 By contrast, higher MeHg concentrations in marine zooplankton exposed

74

to elevated dissolved organic matter concentrations were reported in mesocosm experiments by

75

Jonsson et al.20 In these experiments, the algal community shifted from predominantly

76

autotrophic (phytoplankton) to heterotrophic (bacterioplankton) at high dissolved organic matter

77

concentrations. MeHg increases in zooplankton reflected lengthening of the food web rather than

78

size based effects on phytoplankton community composition. This example illustrates why

79

evaluating and diagnosing the role of simultaneously occurring factors that affect MeHg uptake

80

and trophic transfer is essential for understanding differences in MeHg bioaccumulation across

81

marine ecosystems.

82

Here we use experimental data and bioenergetic theory to develop a model that links

83

ecosystem properties to MeHg uptake in phytoplankton and trophic transfer to zooplankton.

84

Field measurements from coastal, shelf and pelagic regions of the Northwest Atlantic Ocean are

85

used to evaluate the model. We apply the model to quantitatively bound differences in MeHg

86

accumulation attributable to variability in primary productivity, trophic structure, and dissolved

87

organic carbon concentrations across ecosystems.

88 89

Methods

90

Model overview

ACS Paragon Plus Environment

4

Page 5 of 26

Environmental Science & Technology

91

We develop a non-steady state model to calculate changes in MeHg concentrations in

92

phytoplankton and zooplankton due to varying seawater MeHg concentrations and across

93

different life stages for zooplankton. The model solves coupled first-order differential equations

94

for zooplankton MeHg concentrations (wet weight) using one-step forward Euler-type numerical

95

integration. Dissolved MeHg concentrations are used to force the simulation and can be specified

96

deterministically or as distributions (probability density functions). Monitoring data are used to

97

specify probability density functions for temperature, DOC and Chl a distributions in each

98

ecosystem as input parameters that are simulated over 103 Monte Carlo iterations.

99

Phytoplankton

100

Our model includes three size classes of phytoplankton (picoplankton: 0.2-2 µm;

101

nanoplankton: 2 – 20 µm, and microplankton: 20-200 µm. The relative abundance of different

102

phytoplankton size classes is based on empirical relationships with surface (Chl a)

103

concentrations (Supporting Information (SI), Table S1).9

104

Phytoplankton MeHg concentrations are modeled using experimental net MeHg uptake

105

rates from Lee and Fisher5 for several marine phytoplankton species. Concentrations reflect

106

integrated uptake over 4 hours, after which the cells are assumed to achieve equilibrium.5 Lee

107

and Fisher5 find no relationship between MeHg uptake rates and temperature and nutrients but a

108

strong linear correlation between cell surface area (SA) to volume (V) ratios (Figure 1A; Table 1).

109

We combine this information with data showing an exponential decline in uptake (R2 = 0.78)

110

with increasing DOC concentrations from Luengen et al.17 (Figure 1B). The resulting expression

111

for phytoplankton MeHg uptake rate (U) varies as a function of both cell-size and DOC

112

concentrations in seawater (Table 1).

113

Zooplankton

ACS Paragon Plus Environment

5

Environmental Science & Technology

114

Page 6 of 26

We develop different model parameterizations for MeHg accumulation by herbivorous

115

(small) and omnivorous (large) zooplankton based on established theory of growth and diet

116

composition. For both, changes in MeHg concentrations (CZ, ng g-1) are based on first-order rates

117

for: (1) seawater uptake (kB, d-1), (2) dietary intake from consumed prey (kD, d-1), (3) fecal

118

elimination (kE, d-1), and (4) growth dilution (kG, d-1), as follows:

119

 

=  +  − ( +  ) ∙ 

[1]

120

Direct uptake from seawater (bioconcentration) is based on the dissolved MeHg

121

concentration in seawater and an octanol-water partition coefficient for methylmercuric chloride

122

(SI Table S2). Dietary assimilation efficiencies for MeHg have been shown to range between

123

50% and 70%11 and are simulated probabilistically using a uniform distribution. MeHg

124

elimination by zooplankton is parameterized as a function of body burden and seawater

125

temperature (SI Table S2). Ingestion rates for all zooplankton are based on energy needed for

126

growth and metabolic function.21

127

Growth is parameterized differently for herbivorous and omnivorous zooplankton.

128

Herbivorous zooplankton growth is driven by temperature and surface productivity (Chl a)

129

following the relationship developed for marine copepods by Hirst and Bunker (SI Table S2).13

130

The composition of food consumed and associated dietary MeHg uptake is based on seawater

131

filtration and the availability of different size classes of phytoplankton (SI Table S2).

132

Omnivorous (large) zooplankton growth is based on standard bioenergetics that

133

characterize energy requirements for growth (G, g d-1, Table S3). Energy obtained from prey

134

consumption (I, g g-1 d-1, SI Tables S3 and S4) is allocated toward metabolism (respiration (R)

135

and elimination (E), g g-1 d-1), and the difference is used to grow:

136

=

 

= ( −  − ) ∙ 

[3]

ACS Paragon Plus Environment

6

Page 7 of 26

Environmental Science & Technology

137

where M is zooplankton (g) mass and t is time (d). Species-specific von Bertalanffy growth

138

curves from the literature are used to constrain ingestion rates needed for growth.22 Omnivorous

139

zooplankton graze primarily on herbivorous zooplankton and occasionally on large

140

phytoplankton (microplankton).23 Diet composition is weighted toward the largest available prey

141

that falls below a predator-prey length ratio constraint derived from other studies.24–26 Prey size

142

consumed increases with length and gape size (SI Table S5).27 We assume omnivorous

143

zooplankton will preferentially consume microplankton of comparable size to copepods (