Temporal changes in 137Cs concentrations in fish, sediments, and

Oct 15, 2018 - Retention patterns of 137Cs in pelagic fish were comparable to that in laboratory studies of fish in which there were no sustained 137C...
1 downloads 0 Views 679KB Size
Subscriber access provided by REGIS UNIV

Environmental Processes 137

Temporal changes in Cs concentrations in fish, sediments, and seawater off Fukushima Japan Cuiyu Wang, Robert M. Cerrato, and Nicholas S. Fisher Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b03294 • Publication Date (Web): 15 Oct 2018 Downloaded from http://pubs.acs.org on October 19, 2018

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 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 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.

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 29

Environmental Science & Technology

1

Main text of Temporal changes in 137Cs concentrations in fish,

2

sediments, and seawater off Fukushima Japan.

3

Cuiyu Wang, Robert M. Cerrato, and Nicholas S. Fisher

4

1 ACS Paragon Plus Environment

Environmental Science & Technology

5

Temporal changes in 137Cs concentrations in fish,

6

sediments, and seawater off Fukushima Japan

7

Cuiyu Wang*, Robert M. Cerrato, and Nicholas S. Fisher

8

School of Marine and Atmospheric Sciences, Stony Brook University

9

Stony Brook, New York 11794-5000

10 11

Abstract

12 13

We analyzed publicly-available data of Fukushima 137Cs concentrations in coastal fish, in surface

14

and bottom waters, and in surface marine sediments and found that within the first year of the

15

accident pelagic fish lost 137Cs at much faster rates (mean of ~1.3% d-1) than benthic fish (mean

16

of ~0.1% d-1), with benthopelagic fish having intermediate loss rates (mean of ~0.2% d-1). The

17

loss rates of 137Cs in benthic fish in the first year were more comparable to the decline of 137Cs

18

concentrations in sediments (0.03% d-1), and the declines in pelagic fish were more comparable

19

to the declines in seawater. Retention patterns of 137Cs in pelagic fish were comparable to that in

20

laboratory studies of fish in which there were no sustained 137Cs sources, whereas the

21

benthopelagic and benthic fish species retained 137Cs to a greater extent, consistent with the idea

22

that there is a sustained additional 137Cs source for these fish. These field data, based on 13,511

23

data points in which 137Cs was above the detection limit, are consistent with conclusions from

24

laboratory experiments that demonstrate that benthic fish can acquire 137Cs from sediments,

25

primarily through benthic invertebrates that contribute to the diet of these fish.

26 2 ACS Paragon Plus Environment

Page 2 of 29

Page 3 of 29

Environmental Science & Technology

27

Key words

28

cesium; radioactivity; Fukushima

29 30 31

Abstract art

First year after Fukushima accident March 2011- March 2012

Surface Water

Bottom Water

Sediments

Pelagic Fish

32

3 ACS Paragon Plus Environment

Benthopelagic Benthic Fish Fish

Environmental Science & Technology

33

Introduction

34 35

Damage to the Fukushima Dai-ichi nuclear power plant (FDNPP) in March 2011 resulted in a

36

massive release of radionuclides into Japanese coastal waters. Immediately after the accident,

37

various governmental agencies launched coordinated activities to comprehensively monitor the

38

fate of the released radionuclides from FDNPP to the marine environment and fishery products.

39

Until January 2017, over 3000 seawater samples, nearly 2000 sediment samples, and more than

40

10,000 marine fish collected off Fukushima prefecture were analyzed for Fukushima-derived

41

radionuclides. These data are generally available online to the public.

42 43

Surface seawater1, subsurface seawater2, marine sediments3-5 and biota in local5, 6 and distant

44

waters7 were found to be contaminated with cesium isotopes (134Cs and 137Cs). In marine fish, Cs

45

mainly concentrates in muscle8-10 which is the principal part consumed by seafood eaters.

46

Therefore, the relatively long half-life of 137Cs (t1/2 = 30.17 y) has led to continuing concern

47

about public health implications of these released nuclear fission products for seafood

48

consumers11. Fishery activities were restricted in the areas off the shore of Fukushima Prefecture,

49

many fish species were listed as restricted foodstuffs for shipment by the Japan Ministry of

50

Health, Labor and Welfare, and fish exports faced the restriction of foreign governments12. The

51

fishery industry consequently experienced a serious decline after the accident13.

52 53

Buesseler, citing the monitoring results within the first year after the accident, reported that

54

radiocesium concentrations were higher in benthic fish than in pelagic fish14. The Fisheries

55

Agency of Japan15 and several other independent studies6, 16 reported similar findings. A central

4 ACS Paragon Plus Environment

Page 4 of 29

Page 5 of 29

Environmental Science & Technology

56

hypothesis to help explain this pattern of fish contamination has been that marine sediments

57

serve as a continuous source of radiocesium for benthic animals, and benthic fish assimilate Cs

58

through a diet of invertebrates associated with the sediment17. This hypothesis is consistent with

59

the persistent radiocesium contamination in the sediments in the region18 and over 40 species of

60

benthic invertebrates caught in the region19.

61 62

The ocean discharge of radioisotopes from FDNPP peaked on April 6, 2011. In May 2011, the

63

levels of 137Cs in seawater immediately offshore of the power plant had declined by 3 orders of

64

magnitude20, likely due to dilution and dispersal by ocean currents and mixing processes21.

65

About six months after the accident, it was estimated that more than 90% of the radiocesium in

66

sediment had accumulated in regions shallower than 200m4. The concentration of 137Cs in marine

67

sediment was highly heterogeneous22. Most of these studies focused on biota or seawater or

68

sediment contamination separately. For example, spatial and temporal distributions of 137Cs were

69

reported only in sediment3, 4 or only in seawater1, 2, 20. Wada et al. reported the 134+137Cs decline

70

rates in pooled pelagic and benthic fish, and some representative species from both ecological

71

groups6. Iwata et al. compared the biological half-lives of 137Cs in 16 species of benthic

72

invertebrates and fish23. However, only a few studies evaluated the 137Cs contamination in

73

different matrices simultaneously and related 137Cs concentrations in these different components

74

of the marine ecosystem to one another. Sohtome et al. reported similar decline rates of 137Cs in

75

marine benthic invertebrates and sediments5. Tagami and Uchida estimated the contribution of

76

137Cs

in contaminated sediment and seawater to Japanese rockfish Sebastes cheni24.

77

5 ACS Paragon Plus Environment

Environmental Science & Technology

78

No studies have evaluated 137Cs patterns in marine fish of different ecological groups and related

79

them to contamination in ambient water and sediment. Thus, this study explores the patterns of

80

137Cs

81

contamination in surface, bottom seawater and in sediments.

contamination in different groups of fish and compares them with patterns of 137Cs

82 83

In general, marine animals acquire Cs from both aqueous and dietary sources17, 25-28. Depending

84

on prey selection, ingestion rate and contamination levels in water and diet, both aqueous25 and

85

dietary sources17, 25, 26 can contribute a significant portion of Cs to the overall body burden in an

86

animal28. Benthic invertebrates, many being the prey of benthic fish, live in or on marine

87

sediments, and have feeding behaviors tightly associated with sediment exposure. They may

88

serve as important conduits for radiocesium from sediment to benthic fish that feed on them17, 28.

89

Apart from acquiring 137Cs from their diet, bottom-dwelling fish may also acquire 137Cs from the

90

dissolved phase in bottom waters, which is derived from desorbed 137Cs from contaminated

91

sediment. Such routes of metal uptake in benthic animals have been described for a wide variety

92

of other contaminant metals in coastal ecosystems29, 30. However, the degree to which this occurs

93

with 137Cs in marine ecosystems remains poorly understood. Regardless of diet or aqueous

94

pathways, the sediments can potentially serve as the ultimate source of 137Cs for these bottom-

95

dwelling fish in waters off the east coast of Japan.

96 97

The data of 137Cs contamination in fish off Fukushima collected by Japanese authorities are the

98

most comprehensive datasets of radionuclide contamination in marine ecosystems. The goal of

99

this study was to use these valuable data to describe the patterns of 137Cs contamination in

100

individual fish species and in different ecological groups (pelagic, benthopelagic and benthic

6 ACS Paragon Plus Environment

Page 6 of 29

Page 7 of 29

Environmental Science & Technology

101

species) in eastern Japanese coastal waters. Further, it was envisioned that comparing the

102

temporal trends in 137Cs contamination in fish, seawater and sediment could provide important

103

inferences about key sources of 137Cs in fish of different ecological groups. As noted above, prior

104

studies reported declines of 137Cs levels over time in fish23. However, only for benthic

105

invertebrates5 and one species of benthic fish24 were 137Cs declines related to water and sediment

106

concentrations. Our study is the first to relate temporal declines of 137Cs concentrations in marine

107

fish with those in sediments and water in the field, and used a comprehensive field data set to test

108

the hypothesis that contaminated marine sediments have provided an ongoing source of

109

radiocesium for benthic biota. Our statistical analyses focused primarily on field data obtained

110

during the first year after the accident during which time all fish samples had radiocesium levels

111

above reported detection limits.

112 113

Materials and Methods

114 115

137Cs

116

from data available on the websites of the Fisheries Agency of Japan and Nuclear Regulation

117

Authority of Japan. Both 134Cs and 137Cs were decay-corrected to the time of the accident and

118

provided in these databases32,33 . As fish caught off Fukushima had higher 137Cs concentrations

119

than those from the prefectures to the north and the south14, we chose to analyze data only from

120

fish that were landed off Fukushima and from seawater and sediment stations from comparable

121

locations (Fig. 1). 137Cs concentrations were expressed in Bq kg-1 (on a wet weight basis) in both

122

fish and sediment, or in Bq L-1 in seawater32,33.

concentration data in fish, surface seawater, bottom seawater and sediment were collected

123

7 ACS Paragon Plus Environment

Environmental Science & Technology

124

While both Cs isotopes behave similarly, we only analyzed 137Cs because its radiological half-

125

life (30.1 years) is longer than that of 134Cs (2.1 years); consequently, radioactivity levels in biota,

126

seawater and sediment are less subject to radioactive decay during our study period (2011-2016).

127

As Cs mainly concentrates in the muscle of marine organisms8-10, we used muscle tissue

128

concentrations to represent the Cs concentration in fish. We quantified 137Cs loss patterns for the

129

first year after the accident to follow the initial 137Cs contamination in fish directly from the

130

power plant discharge, rather than 137Cs that may have undergone cycling within the marine

131

environment (including seawater and sediments). For analyses that reported lower limits of

132

detection, the 137Cs activities were greater than the detection limit for this radionuclide in all fish

133

captured within one year of the accident. That is, none of the fish samples had 137Cs

134

radioactivities below the detection limits, including for all 36 fish species (regardless of pelagic

135

or benthic habitats) considered in this study within the first year. The first 137Cs non-detects in

136

fish were observed at 381 days after the accident.

137 138

Over the 5-year period following the accident, 137Cs concentrations were analyzed in 102

139

samples of 3 selected species of pelagic fish, 963 samples of 7 selected species of benthopelagic

140

fish, and 7400 samples of 26 selected species of benthic fish. Seawater samples analyzed

141

included 2042 surface and 2104 bottom samples. Over this 5-year period, 91% of pelagic fish

142

samples, 70% of benthopelagic fish, and 60% of benthic fish had 137Cs activities below detection

143

limits; 28% of surface seawater samples and 25% of bottom seawater samples had 137Cs below

144

detection limits. It is recognized that the fish that had 137Cs concentrations below the detection

145

limit may have migrated into the sampling area from relatively uncontaminated regions and may

146

never have been contaminated with 137Cs; thus including them in regression analyses could

8 ACS Paragon Plus Environment

Page 8 of 29

Page 9 of 29

Environmental Science & Technology

147

complicate interpretation of the temporal trends in the 137Cs concentrations, so they were not

148

included in the regression analyses.

149 150

Fish were categorized as pelagic, benthopelagic and benthic fish based on their diet and

151

information provided by Fishbase34. Diet studies of the species of interest are listed in Table S1.

152

Sandlance is placed into the pelagic grouping since it feeds primarily in the water column,

153

although it is recognized that it burrows into marine sediments (Table S1). Species with more

154

than 10 observations were selected to ensure enough data for analysis. A total of 3 pelagic, 7

155

benthopelagic, and 26 benthic fish species were selected.

156 157

To analyze the temporal trend of 137Cs concentrations in fish, seawater and sediment, we used

158

log transformed linear regression analysis. We chose this regression function because

159

exploratory data analysis indicated that 137Cs concentrations in pelagic fish, surface seawater and

160

bottom seawater decreased exponentially over time and had highly skewed residual distributions.

161

To compare the loss coefficients in fish, seawater, and sediments, 137Cs concentrations were

162

natural log transformed and fit to a linear function of time:

163 164

ln(137Cs) = a + bt + 

(Equation 1)

165 166

where b is the loss coefficient (d-1) of the ln(137Cs) concentrations in samples, t represents the

167

time in days since the date of the FDNPP accident on 11 March 2011, and  is the residual,

168

which is assumed to be normally distributed with a mean of 0 and an unknown standard

169

deviation. In order to distinguish interspecific 137Cs depuration patterns within each ecological

9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 29

170

group, regression analysis was conducted for individual fish species. Additionally, regression

171

analyses of pooled species of the three ecological groups were performed to compare differences

172

among these groups.

173 174

To better understand the loss of 137Cs in fish off Fukushima over the first year after the accident,

175

we compared the estimated loss coefficients for different marine fish to depuration estimates

176

from various laboratory studies (Table S2). Published loss coefficients range from 0.008 to 0.05

177

d-1, therefore we considered a range of 0.01 and 0.05 d-1 in the present study.

178 179

Dunnett’s test35, 36 with adjusted degrees of freedom was used to compare the 137Cs loss

180

coefficients in benthic fish with sediments, in benthopelagic fish with surface seawater, bottom

181

seawater and sediments, and pelagic fish with surface seawater. The minimum significant

182

difference between the estimated loss coefficient 𝑏𝑖 of species 𝑖 (𝑖 = 1,…,𝑘) and the coefficient 𝑏𝑠

183

for sediment or seawater is given by

184

1/2

𝑀𝑆𝐷𝑖 = 𝑞𝛼,𝑘,𝑣(𝑠2𝑏𝑖 + 𝑠2𝑏𝑠)

(Equation 2)

185 186

where 𝑠𝑏𝑖 is the standard error of the loss coefficient of species i from a group of k species, 𝑠𝑏𝑠 is

187

the standard error of the loss coefficient in sediment or seawater, and 𝑞𝛼,𝑘,𝑣 is the critical value of

188

the Dunnett test with family error rate 𝛼, number of comparisons 𝑘, and degrees of freedom 𝑣37.

189

To account for heterogeneous residual error variances between fish, sediment, and seawater, the

190

degrees of freedom in the tests were adjusted to

191

10 ACS Paragon Plus Environment

Page 11 of 29

Environmental Science & Technology

192

𝑣=

(𝑠2𝑏𝑖 + 𝑠2𝑏𝑠)2 2

2

(𝑠2𝑏𝑖) (𝑛𝑖

(𝑠2𝑏𝑠)

― 2)

+ (𝑛 𝑠

(Equation 3)

― 2)

193 194

as suggested by Games and Howell38. In this equation, 𝑛𝑖 is the number of individuals of species

195

𝑖 and 𝑛𝑠 is the number of sediment or seawater samples. All analyses were carried out at a

196

family error rate of 0.05. Critical values for Dunnett’s test were obtained from the qNCDun()

197

function in the nCDunnett library of R39.

198 199

Results and Discussion

200

By analyzing the 137Cs concentrations in marine fish, seawater and sediment off Fukushima from

201

2011 to 2016, we gained an overview of 137Cs contamination and attenuation in the ecosystem

202

off Fukushima after the nuclear accident. Long-term contamination of 137Cs in many marine fish

203

species, seawater and marine sediment was observed, but the rate of decline of 137Cs

204

concentrations in fish varied among ecological groups and species during the first year after the

205

FDNPP accident. Generally, benthic fish were shown to have a lower rate of loss of 137Cs over

206

time than other fish.

207 208

Relating 137Cs declines in fish with water and sediment in the first year

209

When data for species in an ecological group were pooled, 137Cs concentrations in pelagic and

210

benthopelagic fish but not benthic fish declined more slowly over the first year than the seawater

211

they contacted (Fig. 2). The loss coefficient in 137Cs concentrations was -1.2% d-1 for the pelagic

212

group (p < 0.001) compared to -2.3% d-1 for surface seawater (p < 0.001) (Fig. 3). The loss

213

coefficient in 137Cs concentrations was -0.2% d-1 for the benthopelagic group (p < 0.05) and -0.9% 11 ACS Paragon Plus Environment

Environmental Science & Technology

214

d-1 in bottom seawater (p < 0.001). There was little difference in the loss coefficient for benthic

215

fish (-0.1% d-1, p < 0.01) and sediments (0.03% d-1), and notably the loss coefficient for

216

sediments was not significant (p = 0.78).

217 218

Loss coefficients were heterogeneous among species within an ecological group, probably

219

reflecting life history and habitat variability (Fig. 3), and while individual pelagic and benthic

220

species tended to reflect the pooled results, differences between individual and pooled species

221

were evident for benthopelagic fish. The loss coefficients in 137Cs concentrations ranged from -

222

0.9 to -1.5% d-1 in pelagic species and all were significantly different from surface seawater at a

223

family error rate of 0.05 (Table S3). Loss coefficients ranged from -0.8 to 0.4% d-1 in benthic

224

species and significantly differed from sediments in only 2 of 26 species (Table S4). In contrast

225

to these two ecological groups, the analysis of individual benthopelagic species did not reinforce

226

the pooled result. Only 1 of 7 species differed significantly from that of sediment, and only 3 of

227

7 species differed significantly from the loss coefficient for bottom seawater (Table S5). Loss

228

coefficients ranged from -0.07 to -0.63% d-1 in benthopelagic species.

229 230

The higher 137Cs concentrations typically observed in benthic fish might be caused by higher

231

initial 137Cs contamination following the accident, or slower depuration of the 137Cs. Most

232

studies focused on the depuration rates6, 24, 40 rather than the initial contamination levels. Our

233

analysis of the data (Figs. S1, S2) showed that, among the fish of the 3 ecological groups

234

considered here, the initial contamination of 137Cs was highest in the pelagic species Japanese

235

sandlance A. personatus (7200 Bq kg-1), followed by the benthic species olive flounder P.

236

olivaceus (2250 Bq kg-1) and the benthopelagic species brown hakeling P. maximowiczi (885 Bq

12 ACS Paragon Plus Environment

Page 12 of 29

Page 13 of 29

Environmental Science & Technology

237

kg-1). This comparison highlighted that the generally higher 137Cs concentrations found in

238

benthic fish was not due to the higher initial 137Cs contamination; therefore, slower net loss of

239

137Cs

is more likely to be the reason. This may be due to the proximity to a sustained source of

240

137Cs

(sediments) for the benthic fish, either through dermal contact, exposure to bottom waters

241

that remain contaminated because of 137Cs desorption from contaminated sediments, or

242

invertebrate prey contaminated from sediment exposure17.

243 244

The average 137Cs concentration in 2011 was highest in pooled pelagic fish, as the initial

245

contamination levels in the Japanese sandlance is very high and the sample size (n=102) is

246

smaller compared to that of the pooled benthopelagic (n=963) and benthic fish (n=7400) groups

247

(Fig. S1). However, in 2012 the average 137Cs concentration in pooled pelagic fish was lower

248

than that in both pooled benthopelagic and benthic species. There were no reports comparing the

249

contamination of 137Cs in pelagic or benthic fish to that in benthopelagic fish, which includes

250

many important commercial species, such as Pacific cod G. macrocephalus, nibe croaker N.

251

mitsukurii, John dory Z. faber, flounder P. japonicus, drum A. argentatus, dory Z. nebulosa and

252

brown hakeling P. maximowiczi. Many studies considered them as benthic species; however,

253

these species have significant portions of their diet from both the pelagic and benthic habitat,

254

thus they should be grouped independently. The average 137Cs concentration was lower in pooled

255

benthopelagic fish than that in benthic fish from 2011 to 2016 (Fig. S1).

256 257

Our analysis showed that during the first year after the accident, the loss coefficients of 137Cs

258

concentrations in marine fish were species-specific and ecological group-specific. While there

259

was considerable interspecific variability, regression analysis showed that within the first year of

13 ACS Paragon Plus Environment

Environmental Science & Technology

260

the accident, pelagic fish lost 137Cs at much faster rates (mean of ~1.3% d-1) than benthic fish

261

(mean of ~0.1% d-1), with benthopelagic fish having intermediate loss rates (mean of ~0.2% d-1).

262

Wada et al. reported radiocesium concentrations in field populations of several species of marine

263

fish in the region for which we calculated loss coefficients.6 These concentrations were similar to

264

those found in this study (Table S6), except for sandlance, possibly because their dataset did not

265

contain samples with extremely high Cs contamination levels. In this study, the decline rate (SE)

266

of 137Cs in Japanese rockfish S. cheni caught off Fukushima was 0.20.3% d-1 (Fig. 3). It was

267

around 0.2% d-1 in the same species caught off two stations that were within 10 km of the

268

FDNPP24. Matsumoto et al. reared 23 individuals of the same species contaminated by the

269

FDNPP accident and caught off Fukushima prefecture, and they reported40 a decline rate (SE)

270

of 0.30.02% d-1. Although the standard errors of the decline rates in our analyses were large

271

(0.3% d-1), the mean value (0.2% d-1) was comparable to that derived from field observations

272

reported by Tagami and Uchida24 and the laboratory study of Matsumoto et al.40

273 274

The rates of decline in 137Cs concentrations in benthic fish off Fukushima were lower than the

275

lower range of those found in laboratory studies (0.8 to 5% d-1, Table S2, Figs. 4, S3, S4),

276

possibly because of a sustained source of radiocesium in Japanese coastal waters and possibly

277

because most of the lab studies used juvenile rather than adult fish, and juveniles with higher

278

metabolic activities41 may be expected to have higher loss rates of 137Cs. Sustained sources could

279

be 137Cs runoff from the land, desorption from contaminated sediments, and consumption of

280

benthic prey contaminated from sediment. Pooling data for many species (especially true for

281

benthic fish species), greatly reduced the confidence intervals for the 137Cs loss rates due to the

282

increase in sample size. Therefore, it is important to acknowledge the wider ranges of the 137Cs 14 ACS Paragon Plus Environment

Page 14 of 29

Page 15 of 29

Environmental Science & Technology

283

loss rates in individual fish of the same ecological group while discussing the different patterns

284

of 137Cs concentrations in the three ecological groups. Thus, pooled data weigh each individual

285

fish of an ecological group equally, regardless of species, and this accounts for differences

286

between mean values among species within an ecological group and pooled values from that

287

group. For example, for benthopelagic species, the loss coefficient of 137Cs from pooled data is

288

0.2% d-1, whereas the average loss coefficient for individual species within this ecological group

289

is 0.37% d-1; similarly, the pooled benthic species had a loss coefficient of 0.1% d-1 and that of

290

individual species is 0.20% d-1 (Fig. 3).

291 292

All but one fish species had estimated loss coefficients that were less than the -1 to -5% d-1 range

293

observed in laboratory studies which primarily used juveniles (Figs. 3, 4); there is a relative

294

paucity of experimental data using adult fish. Of the 36 species analyzed only whitebait, a

295

pelagic fish, had an estimated loss coefficient that exceeded -1% d-1. The remaining species had

296

loss coefficients that ranged from -1 to 0.4% d-1. The loss rates of 137Cs were either larger than or

297

very close to 1% d-1 in pelagic fish (all p < 0.05); in contrast, 137Cs loss rates were all below 1%

298

d-1 in benthopelagic and benthic species (Figs. 3, 4, S3, S4).

299 300

137Cs

301

The 137Cs concentrations in fish, sediments, and seawater that were contaminated over the five

302

years following the March 2011 accident are shown in Fig. S1. Unlike the fish, most of the

303

sediment samples (99.9%) were above detection limits (n = 1978). Detection limits varied

304

considerably among datasets for fish, typically ranging between 5 and 10 Bq kg-1 wet wt, and for

305

water between 1 and 25 Bq L-1; the detection limit was 2 Bq kg-1 wet wt for sediments.

in fish and sediments after the first year

15 ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 29

306 307

The average concentration of 137Cs (SD) in pelagic fish whose 137Cs concentrations were above

308

the detection limit was 3041095 Bq kg-1 in 2011 and 1815 Bq kg-1 in 2012.

309

concentrations above 100 Bq kg-1 were found throughout the first 2.5 years following the

310

accident in benthopelagic fish. Average concentrations of 137Cs (SD) in those benthopelagic

311

fish whose 137Cs concentrations exceeded the detection limit were 69115 Bq kg-1, 3952 Bq kg-

312

1,

313

and 2016, respectively.

314

4 years following the accident in benthic fish. The average concentration of 137Cs (SD) in

315

pooled benthic fish (above the detection limit) was 103169 Bq kg-1, 78 128 Bq kg-1, 3960 Bq

316

kg-1, 2734 Bq kg-1, 1714 Bq kg-1, and 1510 Bq kg-1 in 2011, 2012, 2013, 2014, 2015, and

317

2016, respectively. 137Cs concentrations in representative individual species of the three

318

ecological groups are presented in Fig. S2.

137Cs

2227 Bq kg-1, 1311 Bq kg-1, 114 Bq kg-1, and 71 Bq kg-1 in 2011, 2012, 2013, 2014, 2015, 137Cs

concentrations above 100 Bq kg-1 were found throughout the first

319 320

In the seawater samples above the detection limit, average monthly concentrations decreased by

321

one to two orders of magnitude within the first year after the accident and then declined more

322

gradually thereafter. The average monthly 137Cs concentrations decreased from 12 to 0.08 Bq L-1

323

in surface seawater and from 9 to 1 Bq L-1 in bottom seawater during the first year (Fig. S1).

324

From 2012 to 2015, concentrations gradually decreased from 0.08 to 0.02 Bq L-1 in surface

325

seawater and from 0.1 to 0.01 Bq L-1 in bottom seawater (Fig. S1).

326 327

In contrast to fish and seawater, the decline in the average 137Cs in sediment samples above the

328

detection limit was gradual, but as expected, coastal sediments had much higher (up to 6-fold) 16 ACS Paragon Plus Environment

Page 17 of 29

Environmental Science & Technology

329

average concentrations and variability of 137Cs than offshore stations (Fig. S1). The average

330

137Cs

331

kg-1, 184295 Bq kg-1, 137188 Bq kg-1, and 111194 Bq kg-1 in 2011, 2012, 2013, 2014, and

332

2015, respectively. In sediments of offshore stations, the average 137Cs concentrations (SD)

333

were 8263 Bq kg-1, 7662 Bq kg-1, 6478 Bq kg-1, and 4539 Bq kg-1 in 2011, 2012, 2013, and

334

2014, respectively.

concentrations (SD) in sediments of coastal stations were 502611 Bq kg-1, 309362 Bq

335 336

Factors affecting 137Cs bioaccumulation

337

Benthic invertebrates living in marine sediments constitute the diet of many benthic fish species

338

and can serve as a conduit of radiocesium from sediment to fish17. The decline rates in

339

radiocesium concentrations over the first 2 years after the accident in invertebrates in the coastal

340

benthic food web off Fukushima were very similar to those in the sediments in the same region,

341

suggesting that sediments were a likely source of radiocesium for these benthic invertebrates5.

342

Our analysis showed that the change in 137Cs concentrations in most of the benthic and

343

benthopelagic species was comparable to that in the sediments, indicating that sediment could

344

have been the source of 137Cs for these benthic and benthopelagic fish. Sediments are major

345

repositories30 as well as possible sources42-44 of many metals and radionuclides in aquatic

346

ecosystems. When contaminants are present in a bioavailable form42-44, they can be assimilated

347

by deposit-feeding animals such as polychaetes, clams and fish ingesting sediment and

348

transferred to predators in demersal food chains17. The sediment-water distribution coefficient

349

(Kd) of Cs, which reflects Cs enrichment in sediment relative to water, is 2103 in the marine

350

environment45. Given this moderate Kd value of Cs, some of the 137Cs in the sediment may

351

desorb from sediment particles into overlying seawater and porewater, from which it can 17 ACS Paragon Plus Environment

Environmental Science & Technology

352

presumably be taken up by marine fish. Through these two routes (dietary and aqueous), marine

353

sediments can serve as the ultimate source of 137Cs for marine fish.

354 355

The concentration of 137Cs in fish can be affected by many ecological and biological factors,

356

such as the ambient K+ concentration46, temperature46, 47, trophic level48, animal size47, 48, and

357

diet48. Another hypothesis to explain the slower decline of 137Cs in benthopelagic and benthic

358

fish than in pelagic fish, in addition to sediments serving as a continuous source of 137Cs, is that

359

colder bottom waters (compared to surface waters) leads to slower depuration of 137Cs in

360

benthopelagic and benthic fish. The excretion of elements in animals is a basic metabolic

361

process49, on which temperature could have a profound effect. Many studies have shown that the

362

depuration rate of 137Cs in fish increases with temperature10, 27, 40, 46, 47, 50. Recent and previous

363

Q10 determinations, ranging from 1 to 4 for 137Cs loss from fish, indicate that although 137Cs loss

364

is inversely related to temperature, the influence of temperature is not sufficient to account for

365

the retention of this radionuclide in benthic fish for longer than 2 years (Wang et al. submitted).

366

This highlights the importance of an ongoing source of 137Cs contamination in the fish such as

367

acquisition from contaminated sediments.

368

18 ACS Paragon Plus Environment

Page 18 of 29

Page 19 of 29

369

Environmental Science & Technology

Acknowledgements

370 371

We thank J. Nishikawa and H. Kaeriyama for helping to compile the data on Cs concentrations

372

in fish and two anonymous reviewers for helpful comments on the manuscript. This study was

373

supported by Grants 3423 from the Gordon and Betty Moore Foundation and 269672 from the

374

European Commission to N. Fisher.

375

19 ACS Paragon Plus Environment

Environmental Science & Technology

376

Figure 1. Study area. The circle indicates the location of the Fukushima Daiichi Nuclear Power

377

Plant. The triangles denote the sample locations of sediment and seawater samples.

378

379

20 ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29

Environmental Science & Technology

380

Figure 2. 137Cs concentrations in pooled pelagic, pooled benthopelagic, pooled benthic fish,

381

surface seawater, bottom seawater and sediment during the first year after the accident, and their

382

fit to a linear function of time. Data below the detection limit were excluded. Equations are for

383

regressions of log-transformed data.

384 Pelagic species

Surface seawater

10000 Cs-137 concentration (Bq/L) (Bq L -1)

100.000

y = 305.51e-0.012x R! = 0.44842

1000 100 10 1

0

100

200

300

10.000 1.000 0.100 0.010 0.001

400

y = 53.442e-0.023x R! = 0.92123

!

"! !

Cs-137 concentration (Bq (Bq/L) L-1)

kg-1) (Bq/kg) Cs-137 concentration (Bq

100.000

y = 54.137e-0.002x R! = 0.01751 1000 100 10

!

"! !

#! !

$! !

kg-1) Cs-137 concentration (Bq (Bq/kg)

100 10 1

385

200

0.100 0.010

!

"! !

#! !

$! !

%! !

300

400

10,000

1000

100

1.000

Sediment

y = 66.017e-0.001x R! = 0.00505

0

%! !

y = 18.855e-0.009x R! = 0.87764

10.000

0.001

%! !

Benthic species 10000

$! !

Bottom seawater

Benthopelagic species 10000

1

#! !

300

400

1,000 100 10

y = 113.39e0.0003x R! = 0.0004

1 0

100

200

Time since accident (d)

Time since accident (d)

21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 29

386

Figure 3. Rates of change (% d-1) of 137Cs concentrations in individual species of pelagic,

387

benthopelagic, and pelagic fish, each fish ecological group, surface seawater, bottom seawater,

388

and sediments off Fukushima during the first year after the accident. Data points are means + 95%

389

confidence intervals.

390 22 ACS Paragon Plus Environment

Page 23 of 29

Environmental Science & Technology

391

Figure 4. 137Cs concentrations in representative pelagic, benthopelagic and benthic fish during

392

the first year after the accident. Solid lines indicate regressions of measured values for each

393

species. The dashed lines are the depuration curves of 137Cs in the fish with the y-intercept the

394

same as that derived from the observations off Fukushima, and the lower limit of the loss rates (1%

395

d-1) found in laboratory studies. The dotted lines represent the depuration curves of 137Cs with the

396

y-intercept the same as that derived from the observations off Fukushima, and the upper limit of

397

the loss rates (5% d-1) found in laboratory studies. Anchovy (Engraulis japonicus)

Whitebait

100

Japanese sandlance (Ammodytes personatus)

1000

10000 y = 269.89e-0.015x R! = 0.52987

84.097e -0.009x

y= R! = 0.46531

Pelagic

y = 552.47e-0.01x R! = 0.47517

1000

100 100

10 10

10 1

Benthopelagic

(Bq/kg) Cs-137 concentration (Bq kg-1)

1

!

"! !

#! !

$! !

%! !

Pacific cod (Gadus macrocephalus)

!

"! !

#! !

$! !

Nibe croaker (Nibea mitsukurii)

%! !

1000

1000

100

100

10

10

!

"! !

#! !

$! !

%! !

Brown hakeling (Physiculus maximowiczi)

10000

y = 351.05e-0.006x R! = 0.34513

y = 74.644e-0.003x R! = 0.0561

1

y = 143.73e-0.004x R! = 0.07038

1000 100

1

!

"! !

#! !

$! !

%! !

1

Flathead flounder (Hippoglossoides dubius)

100

10 1 !

"! !

#! !

$! !

%! !

Littlemouth flounder (Pleuronectes herzensteini)

1000

y = 42.242e-0.004x R! = 0.05236

!

"! !

y = 81.605e-0.004x R! = 0.18212

#! !

$! !

Fat greenling (Hexagrammos otakii)

10000

%! !

y = 129.85e-0.001x R! = 0.00482

1000

100

Benthic

100

10 10

1 0

100

200

300

400

10 1

1 0

100

200

300

400

0

100

200

300

Time since accident (d)

398

Exponential regression

1% loss rate

399

23 ACS Paragon Plus Environment

5% loss rate

400

Environmental Science & Technology

400

REFERENCES

401

1. Aoyama, M.; Tsumune, D.; Hamajima, Y. Distribution of Cs-137 and Cs-134 in the North

402

Pacific Ocean: impacts of the TEPCO Fukushima-Daiichi NPP accident. J. Radioanalyt. Nucl.

403

Chem. 2013, 296, 535-539.

404

2. Oikawa, S.; Takata, H.; Watabe, T.; Misonoo, J.; Kusakabe, M. Distribution of the

405

Fukushima-derived radionuclides in seawater in the Pacific off the coast of Miyagi, Fukushima,

406

and Ibaraki Prefectures, Japan. Biogeosciences 2013, 10, 5031-5047.

407

3. Kusakabe, M.; Oikawa, S.; Takata, H.; Misonoo, J. Spatiotemporal distributions of

408

Fukushima-derived radionuclides in nearby marine surface sediments. Biogeosciences 2013, 10,

409

5019-5030.

410

4. Otosaka, S.; Kato, Y. Radiocesium derived from the Fukushima Daiichi Nuclear Power Plant

411

accident in seabed sediments: initial deposition and inventories. Environ. Sci. Process. Impact

412

2014, 16, 978-990.

413

5. Sohtome, T.; Wada, T.; Mizuno, T.; Nemoto, Y.; Igarashi, S.; Nishimune, A.; Aono, T.; Ito,

414

Y.; Kanda, J.; Ishimaru, T. Radiological impact of TEPCO's Fukushima Dai-ichi Nuclear Power

415

Plant accident on invertebrates in the coastal benthic food web. J. Environ. Radioact. 2014, 138,

416

106-115.

417

6. Wada, T.; Nemoto, Y.; Shimamura, S.; Fujita, T.; Mizuno, T.; Sohtome, T.; Kamiyama, K.;

418

Morita, T.; Igarashi, S. Effects of the nuclear disaster on marine products in Fukushima. J.

419

Environ. Radioact. 2013, 124, 246-254.

420

7. Madigan, D. J.; Baumann, Z.; Fisher, N. S. Pacific bluefin tuna transport Fukushima-derived

421

radionuclides from Japan to California. PNAS 2012, 109, 9483-9486.

24 ACS Paragon Plus Environment

Page 24 of 29

Page 25 of 29

Environmental Science & Technology

422

8. Baptist, J. P.; Price, T. J. Accumulation and Retention of Cesium-137 by Marine Fishes. US

423

Government Printing Office: 1962.

424

9. Malek, M.; Nakahara, M.; Nakamura, R. Uptake, retention and organ/tissue distribution of

425

137Cs

426

10. Peters, E. L.; Newman, M. C. Cs-137 elimination by chronically-contaminated largemouth

427

bass (Micropterus salmoides). Health Phys. 1999, 76, 260-268.

428

11. Fisher, N. S.; Fowler, S. W.; Madigan, D. J. Perspectives and reflections on the public

429

reaction to recent Fukushima‐related radionuclide studies and a call for enhanced training in

430

environmental radioactivity. Environ. Toxicol Chem. 2015, 34, 707-709.

431

12. Ministry of Agriculture, Forestry and Fisheries. Lifting and Relaxation of the Import

432

Restrictions on Japanese Foods Following the Accident of Fukushima Daiichi Nuclear Power

433

Plant (54 Countries and Regions). 2017.

434

13. Fisheries Agency of Japan. Japan’s Fishery at a Glance. 2012.

435

14. Buesseler, K. O. Fishing for answers off Fukushima. Science 2012, 338, 480-482.

436

15. Fishery Agency. Report on the Monitoring of Radionuclides in Fishery Products (March

437

2011 - March 2014). http://www.mofa.go.jp/files/000040303.pdf

438

16. Nakata, K.; Sugisaki, H. Impacts of the Fukushima Nuclear Accident on Fish and Fishing

439

Grounds. Springer: 2015.

440

17. Wang, C.; Baumann, Z.; Madigan, D. J.; Fisher, N. S. Contaminated marine sediments as a

441

source of cesium radioisotopes for benthic fauna near Fukushima. Environ. Sci. Technol. 2016,

442

50, 13174-13174.

443

18. Yamamoto, K.; Tagami, K.; Uchida, S.; Ishii, N. Model estimation of Cs-137 concentration

444

change with time in seawater and sediment around the Fukushima Daiichi Nuclear Power Plant

by Japanese catfish (Silurus asotus Linnaeus). J. Environ. Radioact. 2004, 77, 191-204.

25 ACS Paragon Plus Environment

Environmental Science & Technology

445

site considering fast and slow reactions in the seawater-sediment systems. J. Radioanalyt. Nucl.

446

Chem. 2015, 304, 867-881.

447

19. Wada, T.; Fujita, T.; Nemoto, Y.; Shimamura, S.; Mizuno, T.; Sohtome, T.; Kamiyama, K.;

448

Narita, K.; Watanabe, M.; Hatta, N.; Ogata, Y.; Morita, T.; Igarashi, S. Effects of the nuclear

449

disaster on marine products in Fukushima: An update after five years. J. Environ. Radioact. 2016,

450

164, 312-324.

451

20. Buesseler, K.; Aoyama, M.; Fukasawa, M. Impacts of the Fukushima Nuclear Power Plants

452

on marine radioactivity. Environ. Sci. Technol. 2011, 45, 9931-9935.

453

21. Buesseler, K. O. Fukushima and ocean radioactivity. Oceanography 2014, 27, 92-105.

454

22. Higashi, H.; Morino, Y.; Furuichi, N.; Ohara, T. Ocean dynamic processes causing spatially

455

heterogeneous distribution of sedimentary caesium-137 massively released from the Fukushima

456

Daiichi Nuclear Power Plant. Biogeosciences 2015, 12, 7107-7128.

457

23. Iwata, K.; Tagami, K.; Uchida, S. Ecological half-lives of radiocesium in 16 species in

458

marine biota after the TEPCO’s Fukushima Daiichi Nuclear Power Plant accident. Environ. Sci.

459

Technol. 2013, 47, 7696-7703.

460

24. Tagami, K.; Uchida, S. Consideration on the long ecological half-life component of 137Cs in

461

demersal fish based on field observation results obtained after the Fukushima accident. Environ.

462

Sci. Technol. 2016, 50, 1804-1811.

463

25. Mathews, T.; Fisher, N. S. Dominance of dietary intake of metals in marine elasmobranch

464

and teleost fish. Sci. Total Environ. 2009, 407, 5156-5161.

465

26. Metian, M.; Pouil, S.; Hedouin, L.; Oberhansli, F.; Teyssie, J. L.; Bustamante, P.; Warnau,

466

M., Differential bioaccumulation of Cs-134 in tropical marine organisms and the relative

467

importance of exposure pathways. J. Environ. Radioact. 2016, 152, 127-135.

26 ACS Paragon Plus Environment

Page 26 of 29

Page 27 of 29

Environmental Science & Technology

468

27. Pan, K.; Wang, W.-X. Radiocesium uptake, trophic transfer, and exposure in three estuarine

469

fish with contrasting feeding habits. Chemosphere 2016, 163, 499-507.

470

28. Hansman, R.L.; Metian, M.; Pouil, S.; Oberhänsli, F.; Teyssié, J.L.; Swarzenski, P.W. A

471

double-tracer radioisotope approach to assess simultaneous bioaccumulation of caesium in the

472

olive flounder Paralichthys olivaceus. J. Environ. Radioact. 2018, 190, 141-148.

473

29. Bryan, G. W.; Langston, W. J. Bioavailability, accumulation and effects of heavy-metals in

474

sediments with special reference to united-kingdom estuaries - a review. Environ. Poll. 1992, 76,

475

89-131.

476

30. Luoma, S. N.; Rainbow, P. S. Metal Contamination in Aquatic Environments: Science and

477

Lateral Management. Cambridge University Press: 2008.

478

31. Smith, J.; Comans, R.; Beresford, N.; Wright, S.; Howard, B.; Camplin, W. Pollution:

479

Chernobyl's legacy in food and water. Nature 2000, 405, 141-141.

480

32. Fisheries Agency of Japan. 2017. Results of the monitoring on radioactivity level in fisheries

481

products. http://www.jfa.maff.go.jp/e/inspection/index.html.

482

33. Nuclear Regulation Authority of Japan. 2017. Readings of Sea Area Monitoring.

483

http://radioactivity.nsr.go.jp/en/list/205/list-1.html.

484

34. Froese, R.; Pauly, D. Eds. 2018. FishBase. www.fishbase.org

485

35. Dunnett, C. W. A multiple comparison procedure for comparing several treatments with a

486

control. J. Am. Stat. Assoc. 1955, 50, 1096-1121.

487

36. Zar, J. H. Biostatistical Analysis. Prentice Hall: 1999.

488

37. Sokal, R.; Rohlf, F. Biometry: The Principles and Practice of Statistics in Biological

489

Research 3rd ed. WH Freeman and Company. San Francisco, CA, USA 1995.

27 ACS Paragon Plus Environment

Environmental Science & Technology

490

38. Games, P. A.; Howell, J. F. Pairwise multiple comparison procedures with unequal n’s

491

and/or variances: a Monte Carlo study. J. Educat. Stat. 1976, 1, 113-125.

492

39. Broch, S. C.; Ferreira, D. F. Noncentral Dunnett's Test Distribution, 1.1.0; CRAN: 2015.

493

40. Matsumoto, A.; Shigeoka, Y.; Arakawa, H.; Hirakawa, N.; Morioka, Y.; Mizuno, T.

494

Biological half-life of radioactive cesium in Japanese rockfish Sebastes cheni contaminated by

495

the Fukushima Daiichi nuclear power plant accident. J. Environ. Radioact. 2015, 150, 68-74.

496

41. Gillooly, J.F.; Brown, J.H.; West, G.B.; Savage, V.M.; Charnov, E.L. Effects of size and

497

temperature on metabolic rate. Science, 2001, 293, 2248-2251.

498

42. Baumann, Z.; Fisher, N. S. Relating the sediment phase speciation of arsenic, cadmium, and

499

chromium with their bioavailability for the deposit‐feeding polychaete Nereis succinea. Environ.

500

Toxicol. Chem. 2011, 30, 747-756.

501

43. Luoma, S. N. Can we determine the biological availability of sediment-bound trace elements?

502

Hydrobiologia 1989, 176, 379-396.

503

44. Wang, W.-X.; Stupakoff, I.; Fisher, N. S. Bioavailability of dissolved and sediment-bound

504

metals to a marine deposit-feeding polychaete. Mar. Ecol. Prog. Ser. 1999, 178, 281-293.

505

45. IAEA. Technical Report Series 422; Sediment Distribution Coefficients and Concentration

506

Factors for Biota in the Marine Environment. International Atomic Energy Agency: Vienna,

507

Austria 2004.

508

46. Cocchio, L.; Beamish, F.; Rodgers, D. Effects of water chemistry and temperature on

509

radiocesium dynamics in rainbow trout, Oncorhynchus mykiss. Can. J. Fish. Aq. Sci. 1995, 52,

510

607-613.

511

47. Ugedal, O.; Jonsson, B.; Njåstad, O.; Næumann, R. Effects of temperature and body size on

512

radiocaesium retention in brown trout, Salmo trutta. Freshwater Biol. 1992, 28, 165-171.

28 ACS Paragon Plus Environment

Page 28 of 29

Page 29 of 29

Environmental Science & Technology

513

48. Kasamatsu, F.; Ishikawa, Y. Natural variation of radionuclide Cs-137 concentration in

514

marine organisms with special reference to the effect of food habits and trophic level. Mar. Ecol.

515

Prog. Ser. 1997, 160, 109-120.

516

49. Reichle, D. E.; Dunaway, P. B.; Nelson, D. J. Turnover and concentration of radionuclides in

517

food chains. Nucl. Safety 1970, 11, 43-55.

518

50. Doi, H.; Takahara, T.; Tanaka, K. Trophic position and metabolic rate predict the long-term

519

decay process of radioactive cesium in fish: a meta-analysis. PloS one 2012, 7, e29295.

520

51. Takahashi, K.; Masuda, R.; Yamashita, Y. Ontogenetic changes in the spatial learning

521

capability of jack mackerel Trachurus japonicus. J. Fish Biol. 2010, 77, 2315-2325.

522

52. Mathews, T.; Fisher N.S. Trophic transfer of seven trace metals in a four-step marine food

523

chain. Mar. Ecol. Prog. Ser. 2008, 367, 23-33.

524

53. Jeffree, R.A.; Warnau, M.; Teyssie, J.L. Markich, S.J. Comparison of the bioaccumulation

525

from seawater and depuration of heavy metals and radionuclides in the spotted dogfish

526

Scyliorhinus canicula (Chondichthys) and the turbot Psetta maxima (Actinopterygii: Teleostei).

527

Sci. Total Environ. 2006, 839-852.

528

54. Zhao, X.; Wang, W.X.; Yu, K.N.; Lam, P.K. Biomagnification of radiocesium in a marine

529

piscivorous fish. Mar. Ecol. Prog. Ser. 2001, 222, 227-237.

530

55. Bustamante, P.; Teyssié, J.L.; Fowler, S.W.; Warnau, M. Assessment of the exposure

531

pathway in the uptake and distribution of americium and cesium in cuttlefish (Sepia officinalis)

532

at different stages of its life cycle. J. Exp. Mar. Biol. Ecol. 2006, 331, 198-207.

29 ACS Paragon Plus Environment