Time Dependence of the 137Cs ... - ACS Publications

Mar 21, 2016 - accident and deposited on the ground over the North Kanto .... border and shore-lines were obtained from Geospatial Information Authori...
0 downloads 0 Views 2MB Size
Subscriber access provided by Washington University | Libraries

Article 137

Time dependence of the Cs concentration in particles discharged from rice paddies to freshwater bodies after the Fukushima Daiichi NPP accident Kazuya Yoshimura, Yuichi Onda, and Taeko Wakahara Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b05513 • Publication Date (Web): 21 Mar 2016 Downloaded from http://pubs.acs.org on March 24, 2016

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

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

Page 1 of 29

Environmental Science & Technology

1

Time dependence of the 137Cs concentration in

2

particles discharged from rice paddies to freshwater

3

bodies after the Fukushima Daiichi NPP accident

4 5

Kazuya Yoshimura,*,†,‡ Yuichi Onda,‡ Taeko Wakahara,§

6



7

8F 1-29, Okitama-cho, Fukushima-shi, Fukushima 960-8034, Japan

8



9

Tennodai, Tsukuba, Ibaraki 305-8572, Japan

Sector of Fukushima Research and Development, Japan Atomic Energy Agency, Sahei Building

Center for Research in Isotopes and Environmental Dynamics, University of Tsukuba, 1-1-1

10

§

11

Technology, 3-5-8, Saiwai-cho, Fuchu, Tokyo 183-8509, Japan

Institute of Symbiotic Science and Technology, Tokyo University of Agriculture and

12 13

*Corresponding author (Kazuya Yoshimura)

14

Phone: +81-24-529-5560; e-mail address: [email protected]

15

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 29

16

ABSTRACT

17

The concentration of particulate

18

entering the water system, was studied following the Fukushima Daiichi Nuclear Power Plant

19

accident. To parameterize the concentration and to estimate the time dependence, paddy fields

20

covering various levels of

21

121). The particulate 137Cs concentration (kBq kg-SS−1) showed a significant correlation with the

22

initial surface deposition density (kBq m-2). This suggests that the entrainment coefficient (m2

23

kg-SS−1), defined as the ratio between the particulate

24

deposition density, is an important parameter when modeling

25

The entrainment coefficient decreased with time following a double exponential function. The

26

decrease rate constant of the entrainment coefficient was clearly higher than that reported for

27

other land uses and for river water. The difference in the decrease rates of the entrainment

28

coefficient suggests that paddy fields play a major role in radiocesium migration through the

29

water system. An understanding of the decrease rate of the entrainment coefficient of paddy

30

fields is therefore crucial to understand the migration of radiocesium in the water system.

137

137

Cs in paddy fields, which can be a major source of

137

Cs

Cs deposition were investigated over the period 2011–2013 (n =

137

Cs concentration and the initial surface 137

Cs wash-off from paddy fields.

31

ACS Paragon Plus Environment

2

Page 3 of 29

Environmental Science & Technology

32

1. INTRODUCTION

33

Large amounts of radiocesium, especially

34

Nuclear Power Plant (FDNPP) accident and deposited on the ground over the North Kanto and

35

the South Tohoku regions of Japan.1-3 Most of the

36

soil particles,4,5 and a part of this

37

physical processes such as soil erosion. It was then mainly transported in particulate form

38

through the river system.6-9 The migration of particulate 137Cs from terrestrial fields to the water

39

system is a serious issue when predicting changes in environmental contamination levels and for

40

evaluating the effect on water utilization in downstream region.

41

Paddy fields for growing rice are a major land use in South, East, and Southeast Asia. In

42

Fukushima prefecture, paddy fields cover an area of 1053 km2, accounting for 70% of

43

agricultural capacity.10 Large areas of the paddy in Fukushima prefecture have been affected by

44

137

45

rice have been confirmed to be within safety margins, a major proportion of the radiocesium

46

inventory in paddy fields might still be in place.11 Paddy fields have been shown to be one of the

47

most sensitive land uses to soil erosion in the region.9,11-13 Accordingly, they constitute the main

48

source of particle-bound radiocesium entering the rivers. The typical paddy cultivation cycle in

49

Japan and the corresponding field conditions are as follows. In spring, the dry field is flooded

50

with irrigation water, followed by soil puddling, which mixes the paddy soil to a depth of around

51

15 cm. This produces a muddy soil surface covered by highly turbid water with a 5–10 cm depth.

52

The field is then left fallow for 2–5 days as to allow the water level to adjust to 0–5 cm and for

53

soil compaction appropriate for rice planting to take place. The rice seedlings are then planted.

54

The paddy field is kept flooded throughout the irrigation period at a water depth varying

137

137

Cs, were released from the Fukushima Daiichi

137

Cs deposited over ground bound tightly to

Cs subsequently migrates to the water system through

Cs deposition following the FDNPP accident (Fig. 1). While the contamination levels of the

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 29

55

approximately 0–10 cm. The irrigation water is drained off in early autumn, followed by

56

harvesting. There is then a period from the end of autumn to early spring in which the fields are

57

not irrigated. The major causes of soil wash-off from paddy fields are runoff of turbid paddy

58

water during puddling, enforced drainage after puddling, runoff of paddy water during heavy

59

rainfall, and soil erosion during the non-irrigation period. Many studies have used monitoring

60

and modeling to evaluate the soil wash-off from paddy field.11,13,14-17 However, to estimate the

61

particulate

62

concentration in the discharged suspended solid (SS) is essential.

63

The

64

the catchment area to model the rate of

65

between the particulate

66

density in the catchment area have been reported.9,20 An entrainment coefficient (m2 kg-SS−1),

67

defined as the ratio of 137Cs concentration in the eroded soil (kBq kg-SS−1) and the initial surface

68

deposition density (kBq m-2), has been used to model the concentration of 137Cs in eroded soil.21-

69

23

70

surface deposition density can be used to estimate the particulate 137Cs concentration discharged

71

from paddy fields into the water system. Entrainment coefficients were investigated from

72

experimental runoff plots in land uses such as farmland, grassland, bare ground, and forests,

73

following the Chernobyl Nuclear Power Plant accident21-23 and the FDNPP accident,24 but the

74

coefficient has never been established for paddy fields. A decrease in

75

river and lake waters was reported in the period following the Chernobyl Nuclear Power Plant

76

accident,18,19,25,26 suggesting that the entrainment coefficient decreases with time. This decrease

77

is thought to reflect factors affecting radiocesium wash-off, such as susceptibility to soil erosion,

137

137

Cs wash-off from paddy fields, parameterization of the particulate

137

Cs

Cs concentration in river water has been normalized by the surface deposition density in 137

Cs transfer in the river.18,19 Significant correlations

137

Cs concentration in the SS in river water and the surface deposition

These studies suggest that the relationship between the

137

Cs concentration and the initial

137

Cs concentrations in

ACS Paragon Plus Environment

4

Page 5 of 29

Environmental Science & Technology

78

the chemical/physical properties of the soil (i.e., adsorption and desorption of radiocesium), and

79

progressive vertical migration of radiocesium. These factors in turn depended on the uses to

80

which the land is put. It is therefore important to evaluate the entrainment coefficient of paddy

81

fields.

82

The entrainment coefficient should vary with the particle size of the matrix, such as soil and SS,

83

because the specific surface area of the matrix directly affects the 137Cs adsorption capacity. The

84

fine fractions of soil and SS in river water and the sediment in river beds show higher

85

concentrations than those of coarser fractions.7-8, 27-30 The entrainment coefficient has also been

86

shown to vary depending on the particle size of SS in the river water.24 These findings suggest

87

that particle size needs to be considered to parameterize the entrainment coefficient.

88

This study monitored particulate

89

covering various deposition levels in the Fukushima prefecture. The aim of the study was to

90

provide a parameter that can be used to model the wash-off of

91

improve our understanding of the migration of

92

system.

137

137

Cs

Cs concentration in paddy fields extensively in the area

137

Cs from paddy fields and to

137

Cs from agricultural fields into the water

93

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 29

94

95 96

Figure 1. Distributions of 137Cs inventory and paddy fields in Fukushima prefecture. This map

97

was generated by ArcGIS 10 software. The prefectural border and shore-lines were obtained

98

from Geospatial Information Authority of Japan. The 137Cs distribution was derived from the

99

Third Airborne Monitoring Survey by Ministry of Education, Culture, Sports, Science and

100

Technology, Japan.37 The distribution of paddy field in 2009 was based on the National Land

101

Numerical Information service of the Ministry of Land, Infrastructure, Transport and Tourism.39

102 103

2. MATERIALS AND METHODS

ACS Paragon Plus Environment

6

Page 7 of 29

Environmental Science & Technology

104

2-1. Study sites and sampling

105

To estimate the entrainment coefficient in paddy fields at various deposition levels of

106

turbid paddy water was collected 0–2 days after puddling from 14 sites between May 6 and June

107

22, 2012, and from 27 sites between April 29 and May 23, 2013, including 12 of the sites

108

selected in 2012 (Fig. 2). At all sites, cultivation had been carried out continuously, including

109

2011. Sampling was carried out before seedling planting and while sedimentation of SS was in

110

progress. Coarse soil fractions, which are sedimented within 0–2 days of puddling, were not

111

sampled. Although puddling causes vertical and horizontal mixing of the paddy soil, the

112

concentration and inventory have been found to be heterogeneous.31,32 In this study, we therefore

113

collected turbid paddy water at four points in a paddy field and combined these samples.

114

In addition to the extensive sampling of turbid paddy water, particulate

115

the discharged SS due to puddling or rainfall were monitored at six experimental paddy fields to

116

evaluate the sequential decrease in the entrainment coefficient (Fig. 2). Monitoring at site 1 was

117

conducted between September 1 and October 24, 2011 and from May 14, 2012 to September 14,

118

2013. Monitoring at sites 2–6 was carried out from August 19, 2012 to September 23, 2013. The

119

monitoring at site 1 was a continuation of a previous study,11 which monitored two paddy fields,

120

one of which had been decontaminated, while the other one was not. In this study, we analyzed

121

only the data from the un-decontaminated paddy field. The SS discharged from paddy fields by

122

heavy rainfall and puddling was trapped using a time-integrated SS sampler33 installed at an

123

outlet of the paddy field, as shown in Figure 3. The trapped SS was collected in suspension at

124

intervals of 2–4 weeks. The monitoring system was based on the previous study of in site 1.11

125

Samples of turbid paddy water and the suspension were kept still for 2–3 days, and the

126

supernatant was removed. Next, the residue was dried at 105 °C for 24 h for analysis.

137

Cs,

137

Cs

137

Cs concentrations in

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 29

127

128 129

Figure 2. Sites where turbid paddy water was collected after puddling (open circles and double

130

circles) and SS discharge from paddy fields was monitored (open squares). A double circle

131

indicates that the site was used to collect turbid paddy water in both 2012 and 2013, while an

132

open circle denotes a site used only once in either 2012 or 2013. This map was generated by

133

ArcGIS 10 software. The prefectural border and shore-lines were obtained from Geospatial

134

Information Authority of Japan. The 137Cs distribution was derived from the Third Airborne

135

Monitoring Survey by Ministry of Education, Culture, Sports, Science and Technology, Japan.37

136

ACS Paragon Plus Environment

8

Page 9 of 29

Environmental Science & Technology

137 138

Figure 3. Cross-section of paddy field and location of the time-integrated SS sampler. The water

139

level ranged approximately 0–10 cm because of evaporation, weeping, and withdrawal of

140

irrigation water during the irrigation period.

141 142

2-2. Analysis of 137Cs

143

The particulate

144

germanium gamma-ray well detector (GCW2022S, Canberra–Eurisys, Meriden, U.S.A.)

145

equipped with an amplifier (PSC822, Canberra, Meriden, U.S.A.), and a multichannel analyzer

146

(DSA1000, Canberra, Meriden, U.S.A.). Gamma-ray emissions were measured at an energy of

147

662 keV. Analytical accuracy was certified by the World-Wide Proficiency Test using standard

148

soil samples from IAEA.34 Measurement was continued until the counting error in the

149

measurement of 137Cs activity was less than 10%. This required a measurement period of 1–24 h,

150

depending on the sample weight.

151

There were time lags in the date of sampling of turbid irrigation water after puddling in the 2012

152

and 2013 samples. Therefore, the particulate

153

first date of each sampling period. The

154

experimental paddy fields was also decay-corrected from the first date of sampling for each

155

period.

137

Cs concentration (kBq m2 kg-SS−1) was measured using a high-purity

137

137

Cs concentration was decay-corrected from the

Cs concentration in the SS samples from the six

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 29

156 157

2-3. Correction of particle size dependency

158

The particle size dependency of the

159

correction factor (P) based on specific surface area.9,35,36 We measured the particle size

160

distribution using a laser diffraction particle size analyzer (SALD-3100, Shimadzu Co., Ltd.,

161

Kyoto, Japan), and the specific surface area was estimated using a spherical approximation of

162

particles in each size class. To normalize the particle size dependency of the 137Cs concentration

163

in the SS, the concentration was divided by P, which was calculated as follows:

164

137

Cs concentration was corrected using a particle-size

P = (Sr/Ss)v, (1)

165

where Ss is the specific surface area of the criterial sample. In this study, Ss was 0.386 m2 g-dry−1,

166

determined from the SS discharge from the experimental paddy field during puddling at site 1 on

167

June 13, 2011. Sr denotes the specific surface area of each soil sample (tables S1-S2), and v is a

168

constant coefficient with the value 0.65.9,30

169

To ensure valid correction, a typical Ss should be applied. The average P obtained in this study

170

was 1.0 (standard deviation (SD) = 0.41, n = 121), indicating that the Ss value applied was

171

typical for correcting the 137Cs concentration in the SS discharged from the paddy fields.

172 173

2-4. Entrainment coefficient

174

The entrainment coefficient (Sc, m2 kg-SS−1) of 137Cs was calculated as follows:

175

Sc = Ct/A0,

(2)

ACS Paragon Plus Environment

10

Page 11 of 29

Environmental Science & Technology

137

176

where Ct is the particulate

Cs concentration with or without correction for particle size

177

dependency at time t (kBq m2 kg-SS−1), and A0 is the initial surface deposition density of

178

(kBq m−2). The initial surface deposition densities were based on the results of the Third

179

Airborne Monitoring Survey37 as the value on July 2, 2011. These were calibrated using a

180

proportional relationship between air dose rates estimated by the survey and the deposition

181

density of

182

The initial surface deposition density ranged 13–415 kBq m−2 in the paddy fields used to collect

183

turbid paddy water and 61–391 kBq m−2 in the six experimental paddy fields.

137

Cs

137

Cs on the ground surface measured from soil core samples taken at 2,200 points.

184 185

3. RESULTS AND DISCUSSION

186

3-1. Variation in particulate 137Cs concentrations

187

Figure 4 shows the variations in particulate

188

2012 and 2013 with the initial surface deposition density. The increase in uncorrected particulate

189

137

190

deposition density in both the 2012 and 2013 samples. This suggests that the initial surface

191

deposition density is an important factor in the particulate

192

and that the entrainment coefficient is an appropriate metric for estimating the particulate

137

Cs

193

concentration in the discharged SS. This is in line with reports of the particulate

137

Cs

194

concentration of SS in river water.9,20

137

Cs concentration in the turbid paddy waters in

Cs concentrations showed a statistically significant relationship with the initial surface

137

Cs concentration in paddy fields,

195

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 29

196 137

197

Figure 4. Relationship between initial surface deposition densities of

Cs obtained from the

198

Third Airborne Monitoring Survey37 and the particulate

137

199

turbid paddy water in 2012 and 2013. The particulate

137

200

particle size dependency is shown by the closed circles with a solid regression line, and the

201

corrected data is shown by the open circles with a dashed regression line. The particulate

202

concentrations collected from May 6 to June 22, 2012 and from April 29 to May 23, 2013 were

203

corrected for decay from the first date of each sampling period.

Cs concentrations found in the SS of

Cs concentration uncorrected for the

137

Cs

204

ACS Paragon Plus Environment

12

Page 13 of 29

Environmental Science & Technology

205

The particle size dependency of particulate 137Cs concentration of eroded soil, SS in river water,

206

and sediments in riverbeds have been well documented.7-8, 27-30 Particle size varies depending on

207

the location, processes and intensity of soil erosion. Difference in particulate size distribution

208

between the samples were therefore assumed to be reflected in variation in the entrainment

209

coefficient. The particulate size dependency of the particulate 137Cs concentration was corrected

210

and compared with the initial surface deposition density (shown by the open circles in Figure 4).

211

The corrected particulate 137Cs concentration also significantly correlated with the initial surface

212

deposition density. The determination coefficients were higher than those of un-corrected data.

213

The average entrainment coefficient and its standard deviation are given in Table 1. The standard

214

deviation decreased after correcting for the particle size dependency, which has also been

215

reported to improve the correlation between particulate

216

and the initial surface deposition density in the catchment area.9 This suggests that the particulate

217

size dependency of the particulate

218

considered in the parameterization of the entrainment coefficient in terrestrial environments.

219

The averaged entrainment coefficients calculated for the eroded soil samples from the six

220

experimental paddy fields are given in Table 2. The entrainment coefficient of the samples

221

collected at site 1 in the 2011 study11 was also calculated. The particulate

222

corrected for particle size dependency using the particle size distribution re-analyzed in this

223

study, and the initial surface deposition density (310 kBq m-2) estimated by the Third Airborne

224

Monitoring Survey37 were used for the calculation. The entrainment coefficients obtained in this

225

study (after September 23, 2011) varied by almost the same order of magnitude.

137

137

Cs concentration of SS in river water

Cs concentration is an important factor that should be

137

Cs concentration,11

226

ACS Paragon Plus Environment

13

Environmental Science & Technology

Page 14 of 29

227

3-2. Temporal decrease in entrainment coefficient

228

Figure 5 shows the distributions of the entrainment coefficient measured from turbid paddy

229

water in 2012 and 2013. The entrainment coefficient from the previous study, measured at site 1

230

on June 13, 2011 (0.184), is also plotted. The range, median, and averaged values decreased

231

between 2012 and 2013. The entrainment coefficient in 2011 was higher than those in 2012 and

232

2013. These results demonstrate a temporal decrease in the entrainment coefficient.

233

234 235

Figure 5. Boxplot of entrainment coefficients calculated from particulate 137Cs concentrations in

236

turbid paddy water corrected for particle size dependency. The box and the error bars show the

237

ranges from the first to third quartiles and from the minimum to the maximum values,

238

respectively. The line identified the median value, and the averaged entrainment coefficient is

239

represented by the closed circle. The entrainment coefficient calculated for the turbid paddy

240

water sample collected on June 13, 2013 in previous study11 is shown as an open circle.

ACS Paragon Plus Environment

14

Page 15 of 29

Environmental Science & Technology

241 242

Figure 6 shows the temporal decrease in entrainment coefficient taken from Tables 1 and 2. The

243

decrease occurred in two phases: a steep decrease in the first half year after the FDNPP accident,

244

followed by a slower decrease. Although the data 2011 data covered only site 1, the decrease in

245

the entrainment coefficient of

246

function as follows:

247

ܵ௖௧ = ‫ܣ‬଴ × ݁ ି௞భ ௧ + ‫ܤ‬଴ × ݁ ି௞మ ௧

248

where Sct is the entrainment coefficient at time t (yr), and A0 and B0 (m2 kg-SS−1) are constants

249

denoting the initial fractions of the entrainment coefficient that decreased steeply and slowly

250

over time with decrease rate constants of –k1 and –k2 (yr−1), respectively. The decrease rate

251

constants include the decay constant of 137Cs (0.023 yr−1). The data used to calculate A0 and –k1

252

came from only site 1, and the coefficients may be biased by the first one data of June 13, 2011.

253

This adds large uncertainties to the coefficient. The data were included in the analysis, however,

254

because the initial period after the accident is crucially important, and data on this period is

255

scarce. Assuming that the decrease half a year after the accident was represented by the −k2

256

obtained using the results of this study, the parameters of the paddy field were calculated as an

257

A0 of 1.7 m2 kg-SS−1 (σ range from 0.19 to 15 m2 kg-SS−1), a B0 of 0.051 m2 kg-SS−1 (σ range from

258

0.041 to 0.065 m2 kg-SS−1), a −k1 of −14 yr−1 (σ range from −20 to −7.5 yr−1), and a −k2 of −0.48

259

yr−1 (σ range from −0.60 to −0.35 yr−1).

260

The decreases in the

261

represented using a double or triple exponential function.19, 25−26 The results suggest that the

262

decrease in

137

137

Cs over time was derived by fitting a double exponential

(3)

137

Cs concentrations in European and Asian river waters were also

Cs concentration in a water system is characterized by both fast and slow phases.

ACS Paragon Plus Environment

15

Environmental Science & Technology

Page 16 of 29

263

Smith and coworkers used the triple exponential function to represent the decrease based on data

264

from the long-term monitoring over a period of 10 years following the Chernobyl Nuclear Power

265

Plant accident.19 The third decrease rate (i.e. -k3 in their report) was one order of magnitude

266

lower than the -k2 values in our study and in their own report, suggesting that the radiocesium

267

wash-off from paddy field will decrease with time. Therefore, long-term monitoring data is

268

needed to evaluate the rate of decrease in the entrainment coefficient in paddy fields and to

269

predict the migration of 137Cs through paddy fields.

270

271 137

272

Figure 6. Temporal decrease in entrainment coefficients obtained from the particulate

Cs

273

concentration corrected for the particle size dependency. Closed and open circles show the

274

averaged entrainment coefficients obtained from the turbid paddy water after puddling and from

275

the SS discharged from the six experimental paddy fields, respectively. Error bars represents

276

standard deviations. Gray circles and crosses represent the entrainment coefficient obtained at

277

site 1 from June to July 2012 and from the previous 2011 study.11 The double exponential

278

function fitted to the data is shown as a solid line.

ACS Paragon Plus Environment

16

Page 17 of 29

Environmental Science & Technology

279 280

While the temporal trends in the concentration of dissolved and total (dissolved and particulate)

281

137

282

accident,18,19,25,26 the decrease in particulate

283

Furthermore, the decreasing rate constant of the 137Cs concentration depends strongly on the time

284

span over which the constant is calculated.26 As a result, only limited information is available to

285

compare with our results. Temporal variation in the entrainment coefficient has been reported for

286

the eroded soil from experimental plots of farmland, grassland, and forest in the two years after

287

the FDNPP accident.24 The results showed no temporal decrease in the entrainment coefficient,

288

while the coefficient of the paddy fields decreased by approximately one order of magnitude

289

over the same period. The decrease in 137Cs concentration in river water has been observed after

290

the FDNPP accident.7 The entrainment coefficient of SS in water from the Abukuma River,

291

which is the largest river in Fukushima prefecture, has been evaluated from the ratio between the

292

particulate

293

rate of decrease in the entrainment coefficient of SS observed from one to three years after the

294

accident was 0.27 yr-1, or 56% of that obtained for paddy fields. These results suggest a rapid

295

decrease in particulate

296

those in land under other uses and in river water.

297

SS in river water is derived from a diversity of land uses. The entrainment coefficient of SS in

298

the river water decreased with time, while the decrease in the entrainment coefficient has only

299

been observed for paddy fields. This suggests that the paddy field is an important source of SS

300

and particulate radiocesium in river water, accounting for the decrease in the entrainment

301

coefficient of SS in river water. This large contribution by paddy fields to SS migration in the

Cs in river and lake waters were well documented after the Chernobyl Nuclear Power Plant

137

137

Cs concentrations were less precisely evaluated.

Cs concentration and the initial surface deposition density in the catchment.38 The

137

Cs concentrations in SS discharged from paddy fields, compared with

ACS Paragon Plus Environment

17

Environmental Science & Technology

Page 18 of 29

302

water system is consistent with the results of studies monitoring the SS flux from paddy fields.11-

303

13

304

therefore crucial for understanding the migration of radiocesium through the water system.

305

The paddy fields studied, with the exception of site 1, are located outside the evacuation zone,

306

and all the sites, including site 1, have been cultivated, every year including 2011. The

307

countermeasures used against radiocesium contamination of the paddy fields mainly comprised

308

fertilizer management and deep tillage to inhibit the migration of radiocesium into the rice plants.

309

The role of decontamination in the decrease in the entrainment coefficient has therefore been

310

limited. As the paddy fields were kept flooded for several days after puddling, the fine soil

311

fractions, which settle slowly, are assumed to form the surface layer and to be discharged

312

preferentially. Fine soil is well known to contain higher concentrations of 137Cs than coarse soil.7-

313

8, 27-30

314

paddy fields preferentially during the period immediately after the FDNPP accident. Annual

315

puddling would lead to mixing with deeper and less-contaminated soil, further decreasing the

316

particulate 137Cs concentration of the surface soil.31,32 Although additional research is necessary,

317

the preferential discharge of highly contaminated fine soil, and the physical mixing with deeper

318

and less-contaminated soil, can explain the faster rate of decrease in the particulate

319

concentration of paddy fields than that recorded for other land uses.

An understanding of the time dependence of the entrainment coefficient of paddy fields is

Highly contaminated fine soil was therefore assumed to have been removed from the

137

Cs

320

ACS Paragon Plus Environment

18

Page 19 of 29

Environmental Science & Technology

321

TABLES.

322

Table 1. Averaged entrainment coefficient obtained from turbid paddy water and their standard

323

deviation (SD). Both of the uncorrected and corrected data for the particle size dependency are

324

shown.

Sampling dates May 6, 2012 Apr. 29, 2013

Entrainment coefficient (uncorrected) m2 kg-SS-1 Ave SD 0.050 0.062 0.026 0.020

Entrainment coefficient (corrected) m2 kg-SS-1 Ave SD 0.037 0.023 0.023 0.013

325 326

Table 2. Entrainment coefficients obtained from SS discharged from paddy fields and their

327

standard deviation (SD). Both of the uncorrected and corrected data for the particle size

328

dependency are shown. The data from June 13, 2011 to July 27, 2012 was obtained from site 1

329

only. Italics indicate the data calculated for the

330

campaign at site 1.

Sampling dates Jun. 13, 2011 Jun. 20, 2011 Jul. 3, 2011 Jul. 4, 2011 Jul. 15, 2011 Jul. 22, 2011 Jul. 31, 2011 Aug. 6, 2011 Aug. 9, 2011 Aug. 18, 2011 Aug. 19, 2011 Aug. 28 2011 Sep. 23, 2011

Entrainment coefficient (uncorrected) m2 kg-SS-1 Ave SD 0.184 0.036 0.046 0.058 0.041 0.030 0.044 0.034 0.031 0.024 0.039 0.029 0.030

137

Cs concentration obtained in previous

Entrainment coefficient (corrected) m2 kg-SS-1 Ave SD 0.184 - 0.069 - 0.062 - 0.074 - 0.057 - 0.047 - 0.062 - 0.056 - 0.043 - 0.041 - 0.057 - 0.038 - 0.050 -

ACS Paragon Plus Environment

19

Environmental Science & Technology

Oct. 24, 2011 Jun. 1, 2012 Jun. 22, 2012 Jul. 6, 2012 Jul. 18, 2012 Jul. 27, 2012 Sep. 4, 2012 Sep. 17, 2012 Sep. 30, 2012 Oct. 15, 2012 Nov. 10, 2012 Dec. 9, 2012 Jan. 1, 2013 Feb. 5, 2013 May 9, 2013 Jun. 12, 2013 Jul. 14, 2013 Aug. 20, 2013 Sep. 14, 2013

0.021 0.011 0.023 0.017 0.017 0.004 0.027 0.029 0.037 0.029 0.026 0.025 0.034 0.033 0.029 0.022 0.022 0.015 0.015

- - - - - 0.017 0.019 0.027 0.021 0.018 0.016 0.022 0.016 0.025 0.012 0.013 0.012 0.012

0.038 0.017 0.028 0.021 0.021 0.008 0.028 0.031 0.036 0.029 0.027 0.022 0.027 0.025 0.029 0.021 0.022 0.016 0.014

Page 20 of 29

- - - - - - 0.016 0.019 0.023 0.022 0.021 0.014 0.014 0.008 0.018 0.010 0.013 0.010 0.011

331 332 333 334

AUTHOR INFORMATION

335

Corresponding Author

336

* Kazuya Yoshimura

337

Phone: +81-3-3592-2177

338

e-mail address: [email protected]

339

Present Addresses

340

§

341

8F 1-29, Okitama-cho, Fukushima-shi, Fukushima 960-8034, Japan

Sector of Fukushima Research and Development, Japan Atomic Energy Agency, Sahei Building

ACS Paragon Plus Environment

20

Page 21 of 29

Environmental Science & Technology

342 343

Author Contributions

344

K.Y. and J.T. carried out the collections of sample and data and measurement. K.Y. performed

345

analysis. K.Y. and Y.O. discussed the results and contributed to the manuscript preparation.

346 347

Notes

348

The authors declare no competing financial interest.

349 350

ACKNOWLEDGMENT

351

This work was conducted as a part of project entitled the Establishment of grasp method of long-

352

term effects caused by radioactive materials from the Fukushima Daiichi Nuclear Power Plant

353

accident, financially supported by the Ministry of Education, Culture, Sports, Science and

354

Technology and by the Nuclear Regulation Authority, Japan. We are grateful to reviewers for

355

their valuable comments to improve this article.

356 357

ASSOCIATED CONTENT

358

Supporting Information Available

359

Table S1. Specific surface area of SS in the turbid paddy water collected after puddling.

360

Table S2. Specific surface area of SS discharged from six monitoring paddy fields.

361

These materials are available free of charge via the Internet at http://pubs.acs.org.

362 363

ABBREVIATIONS

ACS Paragon Plus Environment

21

Environmental Science & Technology

364

Page 22 of 29

FDNPP, Fukushima Daiichi Nuclear Power Plant; SS, Suspended solid.

365 366

REFERENCES

367

(1) Morino, Y.; Ohara, T.; Watanabe, M.; Hayashi, S.; Nishizawa, M. Episode Analysis of

368

Deposition of Radiocesium from the Fukushima Daiichi Nuclear Power Plant Accident. Environ.

369

Sci. Technol. 2013, 47 (5), 2314-2322; DOI: 10.1021/es304620x.

370

(2) Saito, K.; Tanihata I.; Fujiwara, M.; Saito, T.; Shimoura, S.; Otsuka, T.; Onda, Y.; Hoshi,

371

M.; Ikeuchi, Y.; Takahashi, F.; and et al. Detailed deposition density maps constructed by large-

372

scale soil sampling for gamma-ray emitting radioactive nuclides from the Fukushima Daiichi

373

Nuclear

374

10.1016/j.jenvrad.2014.02.014.

Power

Plant

accident.

J.

Environ.

Radioact.

2015,

139,

308-19;

DOI:

375

(3) Mikami, S.; Maeyama, T.; Hoshide, Y.; Sakamoto, R.; Sato, S.; Okuda, N.; Demongeot, S.;

376

Gurriaran, R.; Uwamino, Y.; Kato, H.; and et al. Spatial distributions of radionuclides deposited

377

onto ground soil around the Fukushima Dai-ichi Nuclear Power Plant and their temporal change

378

until

379

10.1016/j.jenvrad.2014.09.010.

December

2012.

J.

Environ.

Radioact.

2015,

139,

320-343;

DOI:

380

(4) Qin, H.B.; Yokoyama, Y.; Fan, Q.H.; Iwatani, H.; Tanaka, K.; Sakaguchi, A.; Kanai, Y.;

381

Zhu, J.M.; Onda, Y.; Takahashi, Y. Investigation of cesium adsorption on soil and sediment

382

samples from Fukushima Prefecture by sequential extraction and EXAFS technique. Geochem. J.

383

2012, 46 (4), 297-302.

ACS Paragon Plus Environment

22

Page 23 of 29

Environmental Science & Technology

384

(5) Tanaka, K.; Sakaguchi, A.; Kanai, Y.; Tsuruta, H.; Shinohara, A.; Takahashi, Y.

385

Heterogeneous distribution of radiocesium in aerosols, soil and particulate matters emitted by the

386

Fukushima Daiichi Nuclear Power Plant accident: retention of micro-scale heterogeneity during

387

the migration of radiocesium from the air into ground and river systems. J. Radioanal. Nucl. Ch.

388

2013, 295 (3), 1927-1937; DOI: 10.1007/s10967-012-2160-9 2013b.

389

(6) Yamashiki, Y.; Onda, Y.; Smith, H.G.; Blake, W.H.; Wakahara, T.; Igarashi, Y.; Matsuura,

390

Y.; Yoshimura, K. Initial Flux of Sediment-associated Radiocesium to the Ocean from the

391

Largest River Impacted by Fukushima Daiichi Nuclear Power Plant. Sci. Rep. 2014, 4, 3714;

392

DOI: 10.1038/srep03714.

393

(7) Sakaguchi, A.; Tanaka, K.; Iwatani, H.; Chiga, H.; Fan, Q.; Onda, Y.; Takahashi, Y. Size 137

394

distribution studies of

Cs in river water in the Abukuma riverine system following the

395

Fukushima Dai-ichi Nuclear Power Plant accident. J. Environ. Radioact. 2015, 139, 379-389;

396

DOI: 10.1016/j.jenvrad.2014.05.011.

397

(8) Tanaka, K.; Iwatani, H.; Sakaguchi, A.; Takahashi, Y.; Fan, Q. Size-dependent distribution

398

of radiocesium in riverbed sediments and its relevance to the migration of radiocesium in river

399

systems after the Fukushima Daiichi Nuclear Power Plant accident. J. Environ. Radioact. 2015,

400

139, 390-397; DOI: 10.1016/j.jenvrad.2014.05.002.

401

(9) Yoshimura, K.; Onda, Y.; Sakaguchi, A.; Yamamoto, M.; Matsuura, Y. An extensive study

402

of the concentrations of particulate/dissolved radiocaesium derived from the Fukushima Dai-ichi

403

Nuclear Power Plant accident in various river systems and their relationship with catchment

404

inventory. J. Environ. Radioact. 2015, 139, 370-378; DOI: 10.1016/j.jenvrad.2014.08.021.

ACS Paragon Plus Environment

23

Environmental Science & Technology

Page 24 of 29

405

(10) Japan Statistical Yearbook 2015; Ch. 7; Statistics Bureau, Ministry of Internal Affairs and

406

Communications Japan, Tokyo, 2015; www.stat.go.jp/english/data/nenkan/back61/1431-07.htm

407

(accessed Oct. 28, 2015)

408

(11) Wakahara, T.; Onda, Y.; Kato, H.; Sakaguchi, A.; Yoshimura, K. Radiocesium discharge

409

from paddy fields with different initial scrapings for decontamination after the Fukushima Dai-

410

ichi Nuclear Power Plant accident. Environ. Sci. Process. Impacts 2014, 16 (11), 2580-2591;

411

DOI: 10.1039/c4em00262h.

412

(12) Harada, H.; Ota, T.; Shindo, H.; Kobayashi, H.; Ito, C. Estimations and Characteristics of

413

Water Pollutant Load from Large Scale Paddy Fields in Hachirogata Polder. Japanese Journal of

414

Soil Science and Plant Nutrition 2008, 79, 53-60 (in Japanese).

415

(13) Sudo, M.; Miki, T.; Masuda, Y. Research on Characteristics of Turbid Water Effluent

416

from Paddy Fields during the Paddling and the Transplanting Period: Based on a study of paddy

417

watershed in Uso river basin, Japan. Transactions of the Japanese Society of Irrigation,

418

Drainage and Rural Engineering 2009, 77 (2), 113-119 (in Japanese).

419

(14) Yamada, Y.; Igeta, A.; Nakashima, S.; Mito, Y.; Ogasahara, T.; Wada, S.; Ohno, T.;

420

Ueda, A.; Hyodo, F.; Imada, M.; and et al. The runoff of suspended, substances, nitrogen, and

421

phosphorus by enforced draining during the ploughing season; Experiments in paddy fields.

422

Japanese Journal of Limnology 2006, 67, 105- 112 (in Japanese).

423

(15) Matsui, H.; Fukunaga, R.; Shimizu, S.; Noda, K. Suspended sediment runoff from paddy

424

fields in rainfall-runoff events and applicability of USLE. Proceedings of hydraulic engineering,

425

JSCE. 2009, 53, 673-678 (in Japanese).

ACS Paragon Plus Environment

24

Page 25 of 29

Environmental Science & Technology

426

(16) Chen, S.K.; Liu, C.W.; Chen, Y.R. Assessing soil erosion in a terraced paddy field using

427

experimental measurements and universal soil loss equation. CATENA 2012, 95, 131-141; DOI:

428

10.1016/j.catena.2012.02.013.

429

(17) Chen, S.K.; Chen, Y.R.; Peng, Y.H. Experimental study on soil erosion characteristics in

430

flooded terraced paddy fields. Paddy Water Environ. 2013, 11 (1-4), 433-444; DOI:

431

10.1007/s10333-012-0334-2.

432

(18) Santschi, P.H.; Bollhader, S.; Zingg, S.; Luck, A.; Farbenkothen, K. The self-cleaning

433

capacity of surface waters after radioactive fallout. Evidence from European waters after

434

Chernobyl,

435

10.1021/es00074a009.

436

1986-1988.

Environ.

Sci.

Technol.

1990,

24

(4),

519-527;

DOI:

(19) Smith, J. T.; Wright, S. M.; Cross, M. A.; Monte, L.; Kudelsky, A. V.; Saxén, R.;

437

Vakulovsky, S.; Timms, D. N. Global Analysis of the Riverine Transport of

438

Environ. Sci. Technol. 2004, 38 (3), 850-857; DOI: 10.1021/es0300463.

90

Sr and

137

Cs.

439

(20) Tsuji, H.; Yasutaka, T.; Kawabe, Y.; Onishi, T.; Komai, T. Distribution of dissolved and

440

particulate radiocesium concentrations along rivers and the relations between radiocesium

441

concentration and deposition after the nuclear power plant accident in Fukushima. Water Res.

442

2014, 60, 15-27; DOI: 10.1016/j.watres.2014.04.024.

443 444

(21) Konoplev, A. V.; Bulgakov, A. A.; Popov, V. E.; Bobovnikova, Ts. I. Behaviour of longlived Chernobyl radionuclides in a soil-water system. Analyst 1992, 117 (6), 1041-1047.

ACS Paragon Plus Environment

25

Environmental Science & Technology

Page 26 of 29

445

(22) Garcia-Sanchez, L.; Konoplev, A.V.; Bulgakov, A. Radionuclide entrainment coefficients

446

by wash-off derived from plot experiments near Chernobyl. Radioprotection 2005, 40, S519-

447

S524; DOI: 10.1051/radiopro:2005s1-076.

448

(23) Garcia-Sanchez, L.; Konoplev, A.V. Watershed wash-off of atmospherically deposited

449

radionuclides: a review of normalized entrainment coefficients. J. Environ. Radioact. 2009, 100

450

(9), 774-778; DOI: 10.1016/j.jenvrad.2008.08.005.

451

(24) Yoshimura, K.; Onda, Y.; Kato H. Evaluation of radiocaesium wash-off by soil erosion

452

from various land uses using USLE plots. J. Environ. Radioact. 2015, 139, 362-369; DOI:

453

10.1016/j.jenvrad.2014.07.019.

454

(25) Smith, J. T.; Fesenko, S. V.; Howard, B. J.; Horrill, A. D.; Sanzharova, N. I.; Alexakhin, 137

455

R. M.; Elder, D. G.; Naylor, C. Temporal Change in Fallout

Cs in Terrestrial and Aquatic

456

Systems: A Whole Ecosystem Approach. Environ. Sci. Technol. 1999, 33 (1), 49-54; DOI:

457

10.1021/es980670t.

458

(26) Smith, J.T.; Clarke, R. T.; Saxén, R. Time-dependent behaviour of radiocaesium: A new

459

method to compare the mobility of weapons test and Chernobyl derived fallout. J. Environ.

460

Radioact. 2000, 49, 65-83; DOI: 10.1016/S0265-931X(99)00088-0.

461 462 463

(27) Livens, F.R.; Baxter, M.S. Particle size and radionuclide levels in some west Cumbrian soils. Sci. Total Environ. 1988, 70, 1-17; DOI: 10.1016/0048-9697(88)90248-3. (28) Walling, D. E.; Quine, T. A. Use of

137

Cs measurements to investigate soil erosion on

464

arable fields in the UK: potential applications and limitations. J. Soil Sci. 1991, 42, 147-165;

465

DOI: 10.1111/j.1365-2389.1991.tb00099.x.

ACS Paragon Plus Environment

26

Page 27 of 29

Environmental Science & Technology

466

(29) Walling, D. E.; Woodward, J. C. Use of radiometric fingerprints to derive information on

467

suspended sediment sources. In Erosion and Sediment Transport Monitoring Programmes in

468

River Basins; Bogen, J., Walling, D. E., Day, T., Eds.; International Association of Hydrological

469

Sciences Publications: Wallingford 1992; pp 153–164.

470

(30) He, Q.; Walling, D. E. Interpreting particle size effects in the adsorption of 210

471

unsupported

472

DOI: 10.1016/0265-931X(96)89275-7.

137

Cs and

Pb by mineral soils and sediments. J. Environ. Radioact. 1996, 30 (2), 117-137;

473

(31) Harada, N.; Nonaka, M. Soil radiocesium distribution in rice fields disturbed by farming

474

process after the Fukushima Dai-ichi Nuclear Power Plant accident. Sci. Total Environ. 2012,

475

438, 242-247; DOI: 10.1016/j.scitotenv.2012.05.032.

476

(32) Tanaka, K.; Iwatani, H.; Sakaguchi, A.; Takahashi, Y.; Onda, Y. Investigation of Spatial

477

Distribution of Radiocesium in a Paddy Field as a Potential Sink. PLoS ONE, 2013, 8 (11),

478

e80794; DOI: 10.1371/journal.pone.0080794.

479

(33) Phillips, J. M.; Russell, M. A.; Walling, D. E. Time-integrated sampling of fluvial

480

suspended sediment: a simple methodology for small catchments. Hydrol. Process. 2000, 14

481

(14), 2589–2602; DOI: 10.1002/1099-1085(20001015)14:143.0.CO;2-D.

482

(34) Report on the IAEA-CU-2006-03 World-wide Proficiency Test on the Determination of

483

Gamma Emitting Radionuclides (IAEA/AL/171); International Atomic Energy Agency: Vienna,

484

2007; https://nucleus.iaea.org/rpst/Documents/iaea_al_171.pdf. (accessed Oct. 28, 2015)

485 486

(35) Mizugaki, S.; Onda, Y.; Fukuyama, T.; Koga, S.; Asai, H.; Hiramatsu, S. Estimation of suspended sediment sources using

137

Cs and

210

Pbex in unmanaged Japanese cypress plantation

ACS Paragon Plus Environment

27

Environmental Science & Technology

Page 28 of 29

487

watersheds in southern Japan. Hydrol. Process. 2008, 22 (23), 4519-4531; DOI:

488

10.1002/hyp.7053.

489

(36) Yoshimura, K.; Onda, Y.; Fukushima, T. Sediment particle size and initial radiocesium

490

accumulation in ponds following the Fukushima DNPP accident. Sci. Rep. 2014, 4, 4514; DOI:

491

10.1038/srep04514.

492

(37) Corrections of the Results of Airborne Monitoring Surveys by MEXT and Ibaraki

493

Prefecture; Ministry of Education, Culture, Sports, Science and Technology: Tokyo, 2012;

494

http://radioactivity.nsr.go.jp/en/contents/4000/3171/24/1270_083114.pdf. (accessed Jan. 7, 2016)

495

(38) Report for the project of Nuclear Regulation Authority entitled “Establishment of grasp

496

method of long-term effects caused by radioactive materials from the Fukushima Daiichi

497

Nuclear Power Plant accident (in 2013); Japan Atomic Energy Agency: Ibaraki, 2015;

498

http://fukushima.jaea.go.jp/initiatives/cat03/entry06.html (accessed Oct. 28, 2015) (in Japanese)

499 500

(39) Ministry of Land, Infrastructure, National Land Numerical Information Website; http://nlftp.mlit.go.jp/ksj-e/index.html (accessed Oct. 28, 2015)

501 502

ACS Paragon Plus Environment

28

Page 29 of 29

503

Environmental Science & Technology

Table of Contents Graphic

504

ACS Paragon Plus Environment

29