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.20.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.30.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 3041095 Bq kg-1 in 2011 and 1815 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 69115 Bq kg-1, 3952 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 103169 Bq kg-1, 78 128 Bq kg-1, 3960 Bq
316
kg-1, 2734 Bq kg-1, 1714 Bq kg-1, and 1510 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
2227 Bq kg-1, 1311 Bq kg-1, 114 Bq kg-1, and 71 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, 184295 Bq kg-1, 137188 Bq kg-1, and 111194 Bq kg-1 in 2011, 2012, 2013, 2014, and
332
2015, respectively. In sediments of offshore stations, the average 137Cs concentrations (SD)
333
were 8263 Bq kg-1, 7662 Bq kg-1, 6478 Bq kg-1, and 4539 Bq kg-1 in 2011, 2012, 2013, and
334
2014, respectively.
concentrations (SD) in sediments of coastal stations were 502611 Bq kg-1, 309362 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 2103 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