Mercury Temporal Trends in Top Predator Fish of the Laurentian Great

Jun 5, 2017 - Mercury (Hg) concentration trends in top predator fish (lake trout and walleye) of the Great Lakes (GL) from 2004 to 2015 were determine...
0 downloads 0 Views 1MB Size
Subscriber access provided by Binghamton University | Libraries

Article

Mercury temporal trends in top predator fish of the Laurentian Great Lakes from 2004 to 2015: are concentrations still decreasing? Chuanlong Zhou, Mark D. Cohen, Bernard S. Crimmins, Hao Zhou, Timothy A. Johnson, Philip K. Hopke, and Thomas M. Holsen Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 05 Jun 2017 Downloaded from http://pubs.acs.org on June 7, 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 19

Environmental Science & Technology

1

Mercury temporal trends in top predator fish of the

2

Laurentian Great Lakes from 2004 to 2015: are

3

concentrations still decreasing?

4

Chuanlong Zhou1, Mark D. Cohen2, Bernard A. Crimmins1, Hao Zhou3, Timothy A. Johnson3, Philip K. Hopke1,

5

Thomas M. Holsen1, *

6

1

Department of Civil and Environmental Engineering, Clarkson University, Potsdam, NY 13676.

7

2

Air Resources Laboratory, U.S. National Oceanic and Atmospheric Administration, College Park, MD 20740.

8

3

Institute for a Sustainable Environment, Clarkson University, Potsdam, NY 13676.

9

*

Address: 204 Rowley Laboratories

10

Clarkson University

11

PO Box 5710

12

Potsdam, NY 13699-5710

13

Phone: 315-268-3851

14

Fax numbers: 315-268-7636

15

Email address: [email protected]

16

1 ACS Paragon Plus Environment

Environmental Science & Technology

17

Table of Contents Graphic (TOC) Art:

18 19

Abstract

20

Mercury (Hg) concentration trends in top predator fish (lake trout and walleye) of the Great Lakes (GL)

21

from 2004 to 2015 were determined by Kendall-Theil robust regression with a cluster-based age

22

normalization method in order to control for the effect of changes in lake trophic status. Combining data

23

from the GLs (except LE), a significant decreasing trend in the lake trout Hg concentrations was found

24

between 2004 and 2015 with an annual decrease of 4.1% per year, consistent with the decline in regional

25

atmospheric Hg emissions and water Hg concentrations. However, a breakpoint was detected with a

26

significant decreasing slope (-8.1% per year) before the breakpoint (2010) and no trend after the

27

breakpoint. When the lakes are examined individually Lakes Superior and Huron, which are dominated

28

by atmospheric Hg inputs and are more likely to respond to declining emissions from areas surrounding

29

the GL than the lower lakes, have significant decreasing trends with rates between 5.2% and 7.8% per

30

year from 2004 to 2015. These declining trends appear to be driven by decreasing regional atmospheric

31

Hg emissions although they may be partly counterbalanced by other factors including increasing local

32

emissions, food web changes, eutrophication, and responses to global climate change. Lakes Michigan,

33

Erie, and Ontario may have been more impacted by these other factors and have decreasing then non-

34

decreasing or increasing trends.

35

Keywords:

36

Fish Mercury, Great Lakes, Lake Trout, Trends, Age normalization

37 38 39 40 2 ACS Paragon Plus Environment

Page 2 of 19

Page 3 of 19

Environmental Science & Technology

41

Introduction

42

The Great Lakes watershed is a major water resource and fishery for North America.1 For this ecosystem,

43

mercury (Hg) pollution is recognized as a serious environmental and health concern.2-4 Hg can be

44

transformed among elemental, ionic, and organic species when it cycles within and among air, water, land

45

and the biosphere.5 Once transformed into methylmercury (MeHg) in aquatic ecosystems, it becomes

46

more toxic, highly bioaccumulative

47

through fish consumption.8, 9 Currently there are more than 3900 fish consumption advisories for Hg in

48

the U.S. including in the Great Lakes (GL),10, 12

6, 7

and a potential threat to aquatic ecosystems and to human health 11

and Hg is one of the major contaminants causing

49

restriction on consuming fish.

50

2014, which suggests that pregnant, and young children should eat fish with lower mercury concentration

51

in order to gain important developmental and health benefits.13 Hg is one of the analytes of Great Lakes

52

Fish Monitoring and Surveillance Program (GLFMSP), which monitors spatiotemporal trends of

53

bioaccumulative chemicals in the GL using top predator fish as biomonitors.2

54

Atmospheric Hg emissions and subsequent deposition can be strongly correlated with Hg levels in aquatic

55

biota and fish,

56

discharges of Hg containing effluents (e.g. chlor-alkali plant), were largely eliminated in the early 90s.25,

57

26

58

2005, largely due to controls of Hg emissions through state, regional, binational, and voluntary actions.14,

59

27, 28

60

from 2002-2010.29 However the benefits of these reductions could be offset by increasing global Hg

61

emissions which have been largely due to increases from Asian countries. 30

62

The trophic status of a lake, fish species, size, and age are important determinants of fish Hg

63

concentrations.31-33 For example, positive correlations have been reported between Hg concentration and

64

fish size or age of a given species, 6, 34 however, consumption of prey with a higher caloric content could

65

cause a lower Hg concentration in predator fish due to growth dilution.35 These variables need to be

66

analyzed and then controlled for to better determine concentration trends and allow informative lake-to-

67

lake comparisons.32, 36, 37 Fish growth rate has varied significantly in some GL ecosystems because of

68

changing lake-specific food webs caused primarily by invasive species and eutrophication.38, 39, 40 Invasive

69

species can also affect the methylation rate in a lake. For example, increased MeHg concentrations were

70

found in nearshore regions of Lake Michigan with prominent invasive mussel infestations and associated

71

filamentous benthic green algae. 41

72

Other factors related to eutrophication and climate change can also affect fish Hg levels. Since the

73

methylation process occurs primarily in organic-rich particles (in the water column) or surface sediment,

14-16

EPA and FDA issued updated draft advice on fish consumption in June

and are likely the main contemporary source of Hg to the GL17-24 as the direct

Atmospheric Hg emissions in the GL region declined by approximately 50 percent between 1990 and

These and subsequent reductions contributed to observed decreases in atmospheric Hg concentrations

3 ACS Paragon Plus Environment

Environmental Science & Technology

74

42, 43

75

by increasing the extent of anoxic conditions, favorable for methylating bacteria. Hg methylation rates

76

also increase with sediment and water column temperature leading to increasing bioaccumulation44, 45 ,

77

and so a warming climate can impact MeHg concentrations in fish. Climate change may also alter the Hg

78

wet deposition and tributary runoff inputs to the lakes due to changes in precipitation depth. 18 A warmer

79

climate can also increase lake sediment resuspension events due to increased and more intense storms,

80

and decreased ice cover, causing the release of Hg accumulated in sediment.8, 42, 46

81

In this study, Hg concentrations and fish age in top predator fish (lake trout for all of the lakes except for

82

Lake Erie (LE) where walleye are the top predator fish) were measured as part of the GLFMSP from 2004

83

to 2015.

84

Kendall’s tau test and Kendall-Theil robust regression, after age normalization using a cluster-based

85

method, to evaluate the potential impacts of changes in atmospheric Hg inputs, regional Hg emissions,

86

eutrophication, and climate change.

87

Methods

88

Sample collection and total Hg analysis

89

Lake trout (Salvelinus namaycush) samples (600-700 mm) were collected in Lakes Huron (LH),

90

Michigan (LM), Ontario (LO), and Superior (LS). Walleye (Sander vitreus) samples (400-500 mm) were

91

collected in LE. Fifty fish samples were collected annually and grouped into ten composites (based on

92

size similarity) from two sampling sites for each lake, one site for even years (shallow site), and the other

93

one for odd years (deep site). Note that all the sites are in offshore water ecosystems (both industrial area

94

and non-industrial areas in each lake).

95

DMA-80 Direct Mercury Analyzer (Milestone Srl, Bergamo, Italy) based on US EPA Method 7473. Other

96

sample collection, preparation and analysis details are described in the Supplemental Information (SI) and

97

by Zananski et al.2

98

Fish age analysis

99

Fish ages were determined by otoliths, fin clips, coded wire tags (CWT), and scales prior to 2012. A

100

maxillae estimation method, developed at the Michigan Department of Natural Resources, was added in

101

2012 and primarily used due to increased simplicity and accuracy.49 Ages used in this research were an

102

average number of these methods. Significant increasing annualized age trends (%/yr) were found in LH

103

(7.2% even years; 7.4% odd years), LM (4.8% even years) and LS sites (8.3% even years; 6.9% odd

104

years). Additional details regarding the fish age analysis are presented in the SI.

eutrophication can accelerate the production of MeHg by increasing the organic carbon loading, and

Hg concentration trends within and among lakes were determined using non-parametric

47, 48

Composites were kept frozen before being analyzed with a

105 4 ACS Paragon Plus Environment

Page 4 of 19

Page 5 of 19

Environmental Science & Technology

106

Age normalization method

107

The correlations between age and Hg concentration were found to be nonlinear and year-dependent

108

(Figures 1b and S4). These results are unlike previous reported results that could be described using linear

109

or logarithmic equations.

110

linear model or analysis of covariance (ANCOVA),

111

Hg concentration to a consistent age so that all the fish concentrations used for trends analysis likely

112

resulted from the same Hg exposure time, sampling years with similar bioaccumulation patterns were

113

grouped together (clustered) before normalization as follows: 1) an equation was used to describe the

114

annual bioaccumulation patterns; in this case, a linear model was used because it yielded the best fits to

115

the data with an average r2 in each year of 0.35 ranging from 0 ~ 0.85; 2) The distance between each

116

year’s bioaccumulation pattern was calculated based on the equations obtained in step 1; the distance was

117

determined by the average distance of each data point in one group to the other regression line (this

118

approach considers both the distance of data points between the groups and the similarity of the

119

regression lines); 3) hierarchical clustering was performed to define the clusters based on the distances

120

obtained in step 2; as hierarchical clustering will generate new group combinations, the distances among

121

new groups is recalculated in each step of hierarchical clustering; 4) linear regression equations were

122

calculated for each cluster obtained in step 3; and 5) Hg concentrations were normalized to an average

123

age (7.0 years) based on the equations obtained in step 4. A more detailed discussion and the algorithms

124

used are presented in the SI.

125

Trend analysis

126

The non-parametric Kendall’s tau test and Kendall-Theil robust line (Sen’s slope) were used for trends

127

and breakpoint analysis using SPSS 22 (2014 IBM Corporation) and KTRLine software developed by the

128

U.S. Geological Survey (USGS).

129

therefore, earlier fish Hg data (1999-2003) were not included in this analysis. An overall trend analysis

130

for all the lake trout data following age normalization (2004-2015) was performed. In this analysis, all the

131

lake trout data among sites and lakes (except Erie) were combined to develop a regional temporal trend.

132

LE was not included because a different fish species was collected during this period. The changes in fish

133

Hg concentrations are generally smooth, however, a breakpoint analysis can be used when times between

134

sampling periods are long. Trends were also analyzed lake-by-lake with the shallow and deep sites

135

analyzed separately if the age normalized concentrations had significant differences between the sites

136

based on a rank sum test. A more detailed discussion of the trend analysis methodology is presented in the

137

SI.

138

Hg emission inventory

6, 34

Therefore, traditional age normalization approaches, such as the general

54-62

6, 14, 33, 50-53

were not appropriate. To normalized fish

Age data were only available for fish from 2004 to present,

5 ACS Paragon Plus Environment

Environmental Science & Technology

139

U.S. Hg point-source air emissions data were obtained from the National Emissions Inventory (NEI)

140

(1999-2014)

141

emissions data from the National Pollutant Release Inventory (NPRI) were used (2000-2015).65

142

Additional information regarding these emissions datasets is provided in the SI.

143

Results and Discussion

144

Age Normalization. Age normalization was used to control for biases in concentration trends caused by

145

the differences in ages among the sampled fish. As shown in Figures 1 and S4, the clustering results for

146

each site fell into one or multiple clusters. One cluster was found for the LO and LM even years, and LE

147

and LO odd years, indicating that there was no significant change in the Hg bioaccumulation rates at

148

those sites over time.

149

bioaccumulation rates between 2004 and 2015 were significant. Detailed age normalization clustering

150

results including r2 values for each lake are presented in Figure S4 and Table S1. The LM odd years and

151

all of the sites in LH and LS have multiple clusters with clear patterns: 1) the clusters are separated by

152

breakpoint year, which means that Hg bioaccumulation rate in lake trout changed between the early years

153

and recent years; and 2) the slope and intercept of regression lines have decreased in recent years. These

154

results indicate that the recently sampled fish in the upper lakes (LH and LS) have both lower Hg

155

concentrations (smaller intercepts) and less Hg bioaccumulated over time (shallower slopes), indicating

156

that the recent fish have lower Hg concentrations at the same age when compared to the earlier samples.

157

Note that LE even years does not have a breakpoint year, but does have multiple clusters. The change in

158

the Hg bioaccumulation pattern can be caused by changes in Hg inputs,15 and/or changes in the food-web

159

structure and other ecosystem characteristics.35, 41

160

Overall Trend. A significant decreasing trend in lake trout Hg concentrations was found between 2004

161

and 2015 (Kendall’s tau = -0.32, p