Articles | Volume 9, issue 4
Research article
03 Aug 2021
Research article |  | 03 Aug 2021

Assessing the effect of topography on Cs-137 concentrations within forested soils due to the Fukushima Daiichi Nuclear Power Plant accident, Japan

Misa Yasumiishi, Taku Nishimura, Jared Aldstadt, Sean J. Bennett, and Thomas Bittner

Topographic effects on Cs-137 concentrations in a forested area were quantitatively examined using 58 soil core samples collected in a village in Fukushima, Japan, which was directly impacted by the radioactive plume emitted during the 2011 Fukushima Daiichi Nuclear Power Plant (FDNPP) accident. In this study, five topographic parameters and two soil properties were evaluated as controls on the soil Cs-137 concentration using generalized additive models (GAMs), a flexible statistical method for evaluating the functional dependencies of multiple parameters. GAMs employing soil dry bulk density, mass water content, and elevation explained 54 % of the observed concentrations of Cs-137 within this landscape, whereas GAMs employing elevation, slope, and upslope distance explained 47 % of the observed concentrations, which provide strong evidence of topographic effects on Cs-137 concentrations in soils. The model fit analysis confirmed that the topographic effects are strongest when multiple topographic parameters and soil properties are included. The ability of each topographic feature to predict Cs-137 concentrations was influenced by the resolution of the digital elevation models. The movement of Cs-137 into the subsurface in this area near Fukushima was faster in comparison to regions affected by the Chernobyl Nuclear Power Plant accident. These results suggest that the effects of topographic parameters should be considered carefully in the use of anthropogenic radionuclides as environmental tracers and in the assessment of current and future environmental risks due to nuclear power plant accidents.

1 Introduction

On 11 March 2011, a 9.0 magnitude earthquake occurred near the northeast coast of Honshu, the largest island of Japan. A tsunami formed, arriving at the Fukushima Daiichi Nuclear Power Plant (FDNPP) approximately 1 h later. A subsequent power outage caused the water circulation pumps to fail. This then led to overheating of the water, meltdowns, and hydrogen explosions (IAEA, 2015; Mahaffey, 2014). Nearby towns, villages, and farmlands were contaminated by the atmospheric fallout from the radioactive plume. Because human and environmental exposure to levels of elevated radiation can cause adverse health effects (EPA, 2017; Wrixon et al., 2004), the Japanese government designated an area within a 20 km radius of the FDNPP to be evacuated the next day.

Table 1Soil sampling projects conducted in the Fukushima region following the 2011 FDNPP accident (through 2019). NA stands for not available.

Download XLSX

One of the radionuclides released to the atmosphere during the accident was cesium-137 (Cs-137). Cs-137 is an isotope of 55Cs, which emits gamma rays with a half-life of 30.17 years (US EPA, 2017). Because of its long half-life, Cs-137 concentration in the environment is still a concern in Japan. Many studies have been conducted in the Fukushima region since the accident to assess the radioactivity of the land surface and environs (Table 1). Those studies confirmed that the majority of FDNPP-derived Cs-137 was concentrated in the soils close to the ground surface and exponentially decreased with depth, and that Cs-137's subsurface migration speed varied according to land use, soil types, and soil chemistry. Few studies, however, examined how local topography would affect Cs-137 concentrations accumulated by atmospheric fallout or subsequent remobilization on the land surface.

While Cs-137 is a byproduct of the nuclear energy generation process and does not exist in the environment naturally (Amaral et al., 1998; IAEA, 2015; Tsoulfanidis, 2012), it has been used since the 1960s as an environmental tracer to understand surface soil and sediment movement (Walling and He, 1997). There are two primary pathways for Cs-137 once it is released in the environment. Cs-137 lacks a single valence electron and it is positively charged, which enables it to form electrovalent bonds with organic matter and mineral anions. Once Cs-137 is deposited on the ground, it becomes adsorbed by clay minerals (Claverie et al., 2019; Fan et al., 2014; Murota et al., 2016; Nagao, 2016; Nakao et al., 2008; Ohnuki and Kozai, 2013; Park et al., 2019; Ritchie and Ritchie, 1995). It is this adsorption process that makes Cs-137 an effective environmental tracer for modeling surface soil loss, reservoir sedimentation, and sediment yield (Bennett et al., 2005; Loughran et al., 1987; Lowrance et al., 1988; Mabit et al., 2007; Martz and De Jong, 1987; Pennock et al., 1995; Quine et al., 1997; Ritchie and Ritchie, 1995; Wallbrink et al., 2002; Walling et al., 2007; Xinbao et al., 1990). Water also provides an environmental pathway for Cs-137 transport; Cs-137 has a high water solubility and it can attach to sediment within surface waters (Iwagami et al., 2015; Osawa et al., 2018; Sakuma et al., 2018; Tsuji et al., 2016).

Topographic indices such as elevation and slope should play an important role in determining Cs-137 movement in the environment because these indices affect sediment transport and hydrologic processes (Catani et al., 2010; Chen et al., 1997; Griffiths et al., 2009; Hoover and Hursh, 1943; Roering et al., 1999, 2001; Roering, 2008; Rossi et al., 2014; Yang et al., 1998). For example, Zaslavsky and Sinai (1981) noted that surface concavity was the controlling factor in distributing soil water in a catchment. Heimsath et al. (1999) examined the relationships among cosmogenic nuclides, topographic curvature, and soil depth, and concluded that the variable thickness of soil is a function of topographic curvature. Studies using Cs-137 as a tracer found that topography affects runoff and infiltration and, hence, the concentration of Cs-137 in soils (Komissarov and Ogura, 2017; Martin-Garin et al., 2012; Ritchie and McHenry, 1990; Walling et al., 1995).

Figure 1Locations of the FDNPP and the study site (Basemaps: a ESRI, HERE, Garmin, © OpenStreetMap contributors, and the GIS community; b ESRI Japan).

Table 2A comparison of climates, topographies, and soil formations between the FDNPP accident-affected area and the CNPP accident-affected area.

Download Print Version | Download XLSX

Cs-137 accumulation patterns in the Fukushima region can be compared to the eastern European region affected by the Chernobyl Nuclear Power Plant (CNPP) accident in 1986. The CNPP and the FDNPP accidents are the only nuclear power plant accidents categorized as Level 7, which is the highest level on the International Nuclear and Radiological Event Scale (INES; IAEA, 2015; IAEA and INES, 2013; Mahaffey, 2014). The climate and topography of the Fukushima and Chernobyl regions differ (Table 2), and researchers suspect that Cs-137 concentrations on the landscape reflect these differences.

The overall goal of the research presented here is to understand the role of topography in determining the concentrations and distributions of Cs-137 on a landscape affected by the FDNPP accident. The objectives of the current study are as follows: (1) to define the concentrations and distributions of Cs-137 in surface and near-surface soil samples in a forested landscape directly impacted by the FDNPP accident and (2) to quantitatively assess the control of topographic indices and soil properties on Cs-137 concentrations in this forested landscape using predictive analytics. It is envisioned that the primary data collected will contribute toward advancing the knowledge and understanding of the environmental hazards associated with radioactive fallout, societal response, and remediation actions, as well as elucidate the movement and long-term persistence of radionuclides in terrestrial environments.

2 Methods

2.1 Study site topography and forests

Soil samples were collected in Iitate village, a forested region located 35 km northwest of the FDNPP (Fig. 1). The study site was under the plume released on 15 March 2011. The plume initially moved toward the southwest and then, following a change in wind direction, toward the northwest, resulting in wet deposition across northwestern Fukushima and other prefectures (IAEA, 2015; Fig. 2). The Fukushima region is underlain by Paleozoic metamorphic rocks and Paleozoic–Mesozoic igneous rocks (Forest Management Center, 2017). The topography is composed of mountains dissected by narrow streams and covered by deciduous and conifer trees with abundant litter on the forest floor.

Figure 2Timing and locations of the main Cs-137 deposition events following the explosions at the FDNPP (Map: IAEA, 2015).

Figure 3Study site topography and aerial images. (a) Numbered sample locations where multi-year samples were collected at the locations marked by * and contour lines are at 5 m intervals (Basemap: ESRI). (b) A large basin enclosing the longest slope (Basemap: ESRI).

The topography and sampling locations are shown in Fig. 3. The maximum change in elevation of the sampling points is 135 m, and the largest basin area in the study site is 0.56 km2. From late spring to autumn, the trees form a canopy over the hills, and the visibility of the sky from the ground is limited. On the south edge of the largest basin, the forest canopies are thin, and the nearby hills toward the southeast, in the direction of the FDNPP, are visible.

Some parts of the forests in the study site are not native. Overall, 31 % of Fukushima forests are planted (Fukushima Prefecture, 2021). Because no major forestry work had been conducted in the area after the accident, land use history was not considered in the analysis.

2.2 Soil sample collection

Soil samples were collected during the summer in 2016, 2017, and 2018. In 2016, samples were collected at locations 1 to 20, 44, and 45 (Fig. 3a) to cover accessible hillslope areas from the lowlands in a circular pattern. The 2017 sampling campaign sought to confirm anomalies observed in previous year and to check the Cs-137 concentration at the highest elevation by collecting samples at locations 11 to 14, 38, and 46. In 2018, the sampling activity focused on the southwest slope and the eastern side by collecting samples at 11 to 14, 21 to 37, 39 to 43, 45, and 46. Most sampling locations were on the southwest slope due to accessibility. In sum, there were 46 total sampling locations, multiple-year samples were collected at eight locations (averaged herein), and the total number of samples collected was 58. Coordinates were recorded at all sampling locations.

For sample collection, a sampler 5 cm in diameter and 30 cm long from Daiki Rika Kogyo Co., Ltd., Japan, was used. The circular tube was made of metal and contained a replaceable plastic liner. The sampler was tamped into the ground with a hammer. Once fully inserted, the sample was pulled out, and the plastic liner containing the soil was removed, sealed, and taken to the University of Tokyo.

Table 3Background levels of Cs-137 in soils in Japan before the FDNPP accident (the highest concentration at a measured soil depth). NA stands for not available.

Download Print Version | Download XLSX

2.3 Soil property and radioactivity measurements

In the laboratory, each plastic tube was opened, photographed, divided into 2 cm thick disks from the surface to a depth of 20 cm, and 2.5 cm thick disks from depths of 20–30 cm. All visible plant roots or rocks were removed. Each disk sample was placed on a plate, weighed, and dried in an oven for 24 h at 105 C. The drying time was extended as necessary to achieve consistent dryness across the samples. The dried sample was reweighed, placed in a mortar, disaggregated, and particles larger than 2 mm were removed. Mass water content (%) and soil dry bulk density (g cm−3) were calculated for each disk sample. Textural analysis (sand, silt, and clay) was performed on approximately one-half of all soil samples using the pipette method. The samples were selected so that various soil types, colors, locations, and elevations were represented. Organic matter in the soil samples was broken up by placing the sample into a beaker, adding hydrogen peroxide (H2O2), and placing the beaker into a 100 C water bath for 4–5 h. The next day, sodium hexametaphosphate ((NaPO3)6) was added to the samples, which were shaken with a sonic homogenizer. As soon as the sample was placed onto the desk, suspended sediment was collected using a pipette at 20 s for sand (≥0.05 mm) and 4 h for clay (<0.002 mm) to determine the texture class of the sample. Each fraction was dried in the oven and weighed, and its fractional percentage determined. Silt fraction percentage was calculated from the results at 20 s and 4 h. To verify the percentage of clay, selected samples were tested using SALD-7500nano (nano-particle size analyzer, Shimadzu Corporation, Kyoto, Japan).

Processed samples were stored in polyethylene vials and sent to the Isotope Facility for Agricultural Education and Research, University of Tokyo, for analysis. Radioactivity levels were measured with a NaI(Tl) scintillation automatic gamma counter (2480 WIZARD2 gamma counter, PerkinElmer Inc., Waltham, MA, USA), which was equipped with a well-type NaI(Tl) crystal 7.6 cm in diameter and 7.6 cm long, and covered with a 75 mm thick lead shield. Energy calibrations were performed using the 662 keV (kiloelectron volts) energy peak of gamma rays from Cs-137. For radiocesium, the detection limit was approximately 0.5 Bq. After each measurement, the radiation was separated into radiation emitted by Cs-137 and Cs-134 using the abundance ratio of Cs-137 to Cs-134 at the time of sampling. This ratio was obtained from the physical decay rates of the isotopes and the elapsed time from the accident to sampling, assuming that the ratio at the time of the FDNPP accident was 1:1 (Nobori et al., 2013; Tanoi et al., 2019). Although a gamma-ray spectrometer provides more precise Cs-137 measurement data, lower-resolution NaI with effective algorithms can be accurate, and it is preferred for its ruggedness, shorter time required for evaluation, and cost of operation (Burr and Hamada, 2009; Stinnett and Sullivan, 2013). The US Environmental Protection Agency allows the NaI method for gamma-ray measurement (EPA, 2012).

In this study, Cs-137 values are reported in two units: Bq kg−1 or Cs-137 activity measurement per volume, and Bq m−2 or mass depth (Bq kg−1× soil dry bulk density × sample thickness). Since soil bulk density varies across samples, mass depth indicates the functionality of soil bulk density as an explanatory variable, and it provides a means to compare concentration levels among samples with varied soil bulk densities (Kato et al., 2012; Miyahara et al., 1991; Rosén et al., 1999). Cs-137 values among samples were decay normalized to 28 June 2016 (see Appendix A), which enabled a comparison of measurements conducted on different dates, and the sum total of Cs-137 in the entire tube is called “core total”.

Background Cs-137 contamination data before the FDNPP accident were unavailable for the study site. According to previous work, Cs-137 activities in Japanese soils varied from ≤15 to 100 Bq kg−1 before the accident (Table 3). Although the 100 Bq kg−1 measurement was an outlier, this value was used as the conservative background contamination level in the following analysis.

2.4 Digital elevation models

This study employed two digital elevation models (DEMs). The 1 m resolution DEM was provided by the Forestry and Forest Products Research Institute, Japan. Its original datum was GCS JGD 2011 (Zone 9), and the data collection year was 2012. The 10 m resolution DEM was downloaded from the Geospatial Information Authority of Japan (, last access: 13 January 2018) website. Its file date was 1 October 2016, and the original datum was GCS JGD 2000. The coordinate projection of all GIS files, including DEMs, was set to UTM 54N (WGS 1984).

Analysis of geomorphic processes can be affected by spatial resolution (Claessens et al., 2005; Martinez et al., 2010; Miller et al., 2015; Sørensen and Seibert, 2007). Using multiple DEMs enables researchers to identify resolution dependency of specific topographic features (Gallant and Wilson, 1996; Kim and Lee, 2004; Moore et al., 1993).

The topographic parameters for analysis were selected based on the following assumptions.

  • Elevation: soil particles with adsorbed Cs-137 move down a sloped surface, suggesting that Cs-137 concentration in surface soils would be higher at lower elevations (Martin, 2000; Roering et al., 1999).

  • Slope: slope influences the downward mass movement on a surface (Roering et al., 1999).

  • Upslope distance to the basin edge: assuming Cs-137 continuously moves down a sloped surface, the longer the upslope distance, the higher the Cs-137 concentrations (Komissarov and Ogura, 2017; Roering et al., 1999, 2001).

  • Surface plan curvature and topographic wetness index (TWI): Cs-137 moves into the subsurface by infiltration (Schimmack et al., 1989, 1994; Teramage et al., 2014). Cs-137 concentrations would be higher where water ponds, and Cs-137 would migrate to a greater depth in the same temporal period, compared with other locations where water does not pond. Slope profile and surface plan curvature are influential factors in hydrology and soil transport on a sloped surface; hence, they would influence Cs-137 accumulation (Gessler et al., 1995; Heimsath et al., 1997; Momm et al., 2012; Moore et al., 1993; Tesfa et al., 2009). The topographic wetness index (TWI) is frequently used in surface hydrology assessments (Hengl and Reuter, 2008). This index reflects water's tendency to accumulate at any point in a catchment (Quinn et al., 1991), and it is calculated as the natural log of total upslope area divided by the tangent of hillslope gradient (Quinn et al., 1991, 1995).

The elevation of sampling locations was extracted from the DEMs using the raster extract function in the R package (R Core Team, 2015). SAGA GIS (Conrad et al., 2015) was used to compute slope and plan curvature because this function has multiple options for the curvature calculation setting. The method used for plan curvature calculation was the second-order polynomial based on the elevation values in nine surrounding cells (Zevenbergen and Thorne, 1987). The “D-Infinity Distance Up” tool of TauDEM was used to compute the upslope distance (Tarboton, 1997), which calculates the distance from each grid cell up to the ridge cells according to the reverse flow path directions. Here, the “minimum” distance and “surface” (total surface flow path) distance calculation options were used. TauDEM also was used to calculate TWI values. All topographic calculations in SAGA and TauDEM were saved as .tif files. The “raster::extract” function of the R package was used to extract each topographic value for sampling locations.

2.5 Generalized additive models

The sampling locations were influenced by terrain access, so they were not randomly selected or uniform, and the number of samples was small compared with the sampling area size. Thus, distance-based spatial analyses such as semivariograms would not be appropriate, and a regression approach was adopted.

Multiple regression methods (linear, polynomial, logarithm, and bi-splines) were run to check the relationships between Cs-137 measurements and the topographic parameters and the soil properties. The resulting regression plots showed fitting curves with multiple knots, and most of the topographic parameters had non-linear, complex relationships with Cs-137 accumulation patterns. Thus, generalized additive models (GAMs) were employed. GAM replaces the linear predictor of a linear additive model with a sum of smooth functions of predictor variables (non-parametric) (Hastie and Tibshirani, 1990; Wood, 2017; Wood et al., 2016). The smoothing parameter of GAM controls the trade-off between the smoothness of the fit and the closeness of the fit to the data (Wood, 2017; Eq. 1).

(1) g μ i = A i θ + f 1 x 1 i + f 2 x 2 i + f 3 x 3 i + f k x k i ,

where x1, x2, …, xk are explanatory variables (predictors), μiE(expected)(Yi) and YiEF(μi,ϕ)Yi are response variables, EF(μi,ϕ) denotes an exponential family distribution with mean μi and a scale parameter ϕ, Ai is a row of the model matrix (parametric), and θ is a corresponding parameter vector. Finally, fj are smooth functions (non-parametric) of the covariates, xk. In GAM algorithms, models are estimated by minimizing squared errors with an increasing penalty as the curves become less smooth.

GAMs allow flexible specification of the response variable dependence on predictors according to the smooth function rather than parametric relationships (Wood, 2017). Because a smooth function is estimated for each covariate, when multiple predictors are included in a model, it is a sum of these functions and a constant.

GAM users can examine the significance and the model fit of each predictor and the overall model performance with multiple indicators (model knots, effective degree of freedom, F statistic, p value, R2, and deviance-explained percentage). Deviance-explained percentage is the proportion of the total deviance explained by the model. This serves as a generalization of R2 in GAM. While overfitting can be a potential problem of GAM (Friedman and Stuetzle, 1981), GAMs work well when users aim to detect the significance of predictors, and the baseline relationship between the predictors and the response variable in a complex model environment (Hastie and Tibshirani, 1990; Linnik et al., 2020; Tesfa et al., 2009).

3 Results

All primary data collected in the field campaign, processed in the laboratory, and derived from the DEMs are summarized in Appendix B.

Figure 4Distribution of mass-based total Cs-137 in a core sample. (a) A unit of each bin is 200 k. (b) Distribution of total mass depth Cs-137 in a core sample in a log scale (red: Cs-137 kBq kg−1; gray: Cs-137 kBq m−2).


Figure 5(a) Depth distributions of Cs-137 concentrations in kBq kg−1 (red) and kBq m−2 (gray). Standard deviations (red) and COVs (gray) by depths in (b) kBq kg−1 and (c) kBq m−2.


3.1 Cs-137 concentration distributions among core samples

Figure 4a shows the binned distributions of core total Cs-137 values in kBq kg−1 (red) and kBq m−2 (gray), as well as their means (Fig. 4b in a log scale). When measured in kBq kg−1, the Cs-137 values clustered in six bins, and 41 % of the core samples were in the 100–300 kBq kg−1 interval (Fig. 4a: red). When measured in kBq m−2, Cs-137 values spread across 13 bins (Fig. 4a: gray). The average total core Cs-137 values among these samples were 267 kBq kg−1 and 1075 kBq m−2.

Figure 6Depth profiles of (a) mass water content (%) and (b) soil dry bulk density (g cm−3).


Cs-137 activities exponentially decreased with depth in both kBq kg−1 and kBq m−2 (Fig. 5a). These vertical profiles are similar to the results reported in previous studies (Fujii et al., 2014; Koarashi et al., 2012; Matsunaga et al., 2013; Takahashi et al., 2015; Tanaka et al., 2012; Teramage et al., 2014). The current samples were collected 5–7 years after the FDNPP accident, and the soils still held the largest amount of Cs-137 in the uppermost layer. The average depth of 90 % concentration (90 % of the entire Cs-137 in a core sample) was 6.1 cm in kBq kg−1 and 6.9 cm in kBq m−2. Standard deviations of Cs-137 activities were up to 115 kBq kg−1 and 335 kBq m−2 in the upper 2 cm depth (black dots in Fig. 5b and c). The relative variability among samples (coefficient of variations; COVs) was below 100 % in the top layer, and COVs increased toward the 14–16 cm depths and decreased again for both units (gray dots in Fig. 5b and c).

3.2 Soil properties

Figure 6 displays the averages of mass water content (%) and soil dry bulk density (g cm−3) for each depth layer. Gravel (below the mid-depth from two disks) and visible roots (at the ground surface from several samples) were removed, which decreased the dry bulk densities of select samples.

Soil texture affects the amount of adsorbed Cs-137 per unit mass of soil particles and accumulation patterns of Cs-137 in soils (Bennett et al., 2005; Giannakopoulou et al., 2007; Korobova et al., 2016; Walling and Quine, 1992). On average, the tested soils contained >50 % sand, and most of the samples were in the categories of sand, loamy sand, sandy loam, and loam (Appendix B).

The average mass water content percentage at the ground surface was above 100 % (Fig. 6a). The standard deviations of water content in the uppermost soil layers were >50 % and the deviation percentages decreased with depth. Some samples collected near the surface exceeded field capacity, likely due to poor drainage, soil texture, and/or high organic content.

Soil dry bulk density and its standard deviation increased with depth (Fig. 6b). The average dry bulk density of all disk samples was 0.45 g cm−3, and the densities ranged from 0.06 to 1.43 g cm−3. The Japanese agricultural soil profile physical properties database, Solphy-J, compiled data for 1800 categories of Japanese soils (Eguchi et al., 2011). This database shows that the soil dry bulk densities of Japanese orchard soils range from 0.80–1.44 g cm−3. Fujii et al. (2014) collected soil samples to 20 cm depth at five locations in Fukushima forests, and these soil dry bulk densities ranged from 0.3–0.7 g cm−3. The relatively low bulk density in this study agrees well with previous work.

3.3 Cs-137 concentration distribution on representative slope

Figures 7 and 8 organize the core total Cs-137 concentration of all samples by elevation and slope extracted with the 10 m DEM. Migration head depths have been added to Fig. 7, defined as the depth at which Cs-137 activity (Bq kg−1) decreased to the conservative background activity level (i.e., 100 Bq kg−1). This study considers the head depth as the deepest subsurface migration point of the FDNPP-derived Cs-137. These figures provide an overview of Cs-137 distribution patterns on the hillslopes. Numbers 1–4 are added to the graphs to indicate the highest concentration samples in Bq m−2. The elevation data were extracted from the 10 m resolution DEM.

Figure 7Core total Cs-137 in both units along the elevation (m). Numbers 1–4 indicate the samples with the four highest Cs-137 concentrations in kBq m−2. Migration head depths (cm) are the depths at which Cs-137 concentrations decreased to 100 Bq kg−1 in a core sample.


Figure 8Slope (degrees) of sampling points along the elevation (m). The numbers 1–4 correspond to the same numbers in Fig. 7.


Cs-137 activities in kBq m−2 units increase with elevation (p value is 0.03), but this trend was not observed for Cs-137 concentrations in kBq kg−1 (Fig. 7). The migration head depths of three of the samples with the highest Cs-137 concentrations had reached the end of the sampler (30 cm), and the average migration depth among all samples was 21.7 cm. Two of the four highest Cs-137 concentrations were at the highest elevation and near the bottom of a hill (Fig. 7), and they were locations with relatively low slopes (Fig. 8). The remaining two samples with the highest Cs-137 concentrations were at 24–26 slopes, which are two of the steepest sampling locations. Elevation and slope do not appear to explain the distribution of Cs-137 concentration, suggesting that additional topographic parameters may be required.

3.4 GAM assessment

Model assessment will focus on eight parameters: Cs-137, five topographic parameters, and two soil properties. This study used mgcv in the R package (Wood et al., 2016) for GAM calculations. The model setting was consistent for all model runs, except rotating predictors (see Appendix A).

First, GAMs were run with a single topographic parameter or one of the soil properties. Then models with parameter combinations were run to assess the interactions among the parameters. Combinations with more than three parameters were not tested to avoid overfitting and to identify parameters with distinct influences on Cs-137 concentrations. The total number of possible combinations was 63 – combinations of seven parameters, including five topographic parameters (elevation, slope, upslope distance, plan curvature, and TWI), and two soil properties (water content and soil dry bulk density), not considering the parameter orders (Fig. 9). Among the 63 combinations, three were only soil properties (combinations of two properties), 25 were only topographic parameter(s) (combinations of five predictors), and the remaining 35 were mixed soil properties and topographic parameters. Then, the same GAM with topographic parameter(s) was executed four times against Cs-137: topographic parameters extracted from the (1) 1 m and (2) 10 m DEMs, and with units of (3) Bq kg−1 and (4) Bq m−2. Soil parameters are not affected by DEM resolutions (Fig. 9).

Figure 9GAM parameter combinations.


GAMs were executed against each depth layer to check model performance as affected by changes by depth. Topographic parameters have one distinct value for an entire core sample, such as elevation and slope. Thus, when a model included a topographic parameter, the model check was a “one-to-many” relationship (topographic parameters to Cs-137 in multiple depth layers). Soil properties had a measurement for each depth layer. Thus, when a model included a soil property, the prediction was a “one-to-one” relationship (soil property in a depth layer to Cs-137 in the depth layer).

GAM accuracies were evaluated by three indicators: Akaike information criterion (AIC), generalized cross validation (GCV), and R2 (see Appendix A). These model parameters were extracted from mgcv outputs in the R package. The AIC is a weighted sum of the log likelihood of the model and the number of fitted coefficients. The lower the AIC indicates a better model fit (Bivand et al., 2008). GCV is a simplified version of cross validation that checks model fit by removing one data point at a time. GCV can be used in a similar way as AIC to measure the relative performance of multiple models.

Table 4The results of single topographic parameter GAMs for deviance-explained percentages and the number of depth layers where p values were equal to or less than 0.05. “( )” indicates standard deviation.

Download Print Version | Download XLSX

Table 5The results of single soil property parameter GAMs for deviance explained and the number of depth layers where p values were equal to or less than 0.05. “( )” indicates standard deviation.

Download Print Version | Download XLSX

Identifying outliers is vital for understanding model errors and capturing influential factors not included in the model. Cook's distance method was used for outlier identification (Cook, 1977; see Appendix A). Cook's distance identifies influential data points by checking how the fitted model parameters change when an outlier data point is removed. Cook's outliers were extracted from linear regression summary plots using the R package (Stevens, 1984).

There was no reference data from other sites to validate the model prediction accuracy. Therefore, Cs-137 concentrations were reverse-predicted by applying the two best-performing prediction models to the original data. Here, mgcv's “predict.gam( )” in the R package was used to generate the predictions.

3.5 GAM performance

In Tables 4 and 5, the percentage of deviance explained shows the average values for the 30 cm depth for single-parameter model results. The p value columns show the number of depth layers in which the average p values were equal to or less than the significance threshold (p≤0.05; e.g., “4/14” means that four depth layers had p≤0.05). Table 4 includes the model results with each topographic parameter that varied with DEM resolution. Table 5 displays the model results with each soil property parameter, whose values were unaffected by the DEM resolution. None of the single-parameter models returned average deviance-explained percentages greater than 36 %, which was the model with water content. The lowest average p value was 0.05, which was also the model with water content. Among topographic parameters, TWI was the most effective predictor with the 1 m DEM (27 %), and slope was the most effective predictor with the 10 m DEM (28 %). Elevation and slope explained Cs-137 deviance better with 10 m DEM than with 1 m DEM. Topographic parameters explained Cs-137 deviance at higher percentages in Bq kg−1 units than in Bq m−2 units, except plan curvature. Soil properties explained Cs-137 deviance at higher percentages than did the topographic parameters. Like the topographic parameters, the soil properties model explained Cs-137 deviance at higher percentages in Bq kg−1 units than in Bq m−2 units.

Figure 10Deviance-explained percentages using the 1 m DEM. (a) All parameter including combinations soil properties and topographic parameters. (b) Parameter combinations with only water content (%) and soil dry bulk density (g cm−3). (c) Parameter combinations with only topographic parameters.


Figure 11Deviance-explained percentages using the 10 m DEM. (a) All parameter combinations soil properties and topographic parameters. (b) Parameter combinations with only water content (%) and soil dry bulk density (g cm−3; the same graphic in Fig. 10 is repeated for the comparison purpose). (c) Parameter combinations with only topographic parameters.


Figures 10 and 11 display the deviance-explanation percentages of GAMs with combinations of parameters, where the red and gray dots represent the deviance-explanation percentages of Cs-137 in Bq kg−1 and Bq m−2 units, respectively. The average correlation index between the average deviation explanation percentages and p values throughout the depth was −0.64. That is, higher explanation percentages returned lower p values.

Table 6The most effective parameter combinations with the best deviance-explained percentages and the number of depth layers where p values were equal to or less than 0.05. “( )” indicates standard deviation. The two models given in bold, [A] and [B], returned the best model performance with the 10 m DEM and were used for reverse prediction accuracy checks.

* An abbreviation of “topography”.

Download Print Version | Download XLSX

The combination parameter models explained Cs-137 deviance in Bq kg−1 (red dots) at higher explanation percentages on average than deviance in Bq m−2 (gray dots) for both DEMs, and all six figures showed similar vertical S curves (Figs. 10 and 11). The deviance-explained percentages of some models were over 75 % for the top layer; the percentages explained decreased with depth, then increased again toward the 12–14 cm depth. These S-shaped curves resemble the vertical profiles of the COVs of Cs-137 measurements (Fig. 5b and c), and the curves were observable for models with only soil properties and only topographic parameters.

The most effective parameter combinations for explaining Cs-137 deviance are summarized in Table 6. The rows are separated into the model combinations with a mix of soil and topographic parameters, and the model combinations with only the topographic parameters. The most effective combinations for predicting Cs-137 deviance were “water content + bulk density + elevation” for models with both 1 m (53.56 %) and 10 m (54.03 %) DEMs. The most effective topographic parameter combinations were “slope + plan curvature + TWI” when using the 1 m DEM (40.62 %) and “elevation + slope + upslope distance” when using the 10 m DEM (46.70 %). The combination models returned higher explanation percentages using the 10 m DEM and Cs-137 in Bq kg−1 units.

Two models with the highest explanation percentages in Table 6 are identified as Model A and Model B. These will be used for reverse prediction accuracy checks.

Figure 12(a, b) Regression analysis of actual and predicted Cs-137 activities. (c, d) Model fit diagnosis results of the best-performing models (light gray: AIC; black: GCV score; dashed: deviance-explained percentage: loosely dash-dotted: R2). Higher deviance-explanation percentages and R2 values indicate how well the model explained, or approximated, the real data points. Lower AIC indicates a good model with less information loss (parsimony), and lower GCV suggests a good model with a smaller prediction error (Hastie et al., 2009; Wood, 2017).

The two models in Table 6 that had the best performance were applied to the original data to reverse predict Cs-137 concentrations with the 10 m DEM, and model accuracies were checked. Regression analyses between the actual and predicted Cs-137 values returned R2=0.54 for Model A with soil properties and elevation (Fig. 12a). Model B with only the topographic parameters returned lower R2=0.25 (Fig. 12b), overestimating two data points at 400–1000 Bq kg−1. With increased depth, model fit (AIC and GCV) continues to improve while the Cs-137 standard deviations decrease (Fig. 12c and d). The explanation power of the model (deviance explained and R2) was the highest at the mid-depth region, where the relative variance of Cs-137 was the largest (Fig. 5b and c).

Table 7Comparison of Cs-137 downward migration depth between the CNPP accident-affected area and Fukushima at about the same length of time after the nuclear plant accidents. Values in the current study are not decay normalized.

Download Print Version | Download XLSX

Four outliers were identified on the basis of Cook's criterion. Two of the outliers (samples 37 and 38 in Fig. 3) were at the highest elevation on the basin ridges with a flat surface. The third outlier (sample 23 in Fig. 3) was at a location with a lower TWI (in the first quantile) and a 26 slope. This sample did not show extreme values in Cs-137, water content, or soil dry bulk density, but the migration head depth of the sample was one of the deepest at 16 cm. The fourth outlier (sample 32 in Fig. 3) was near the basin edge and at a location with a 15 slope, and concave curvature. This sample had the highest water content in the uppermost layer (0–2 cm). Outlier samples 23 and 37 in Fig. 3 had the highest Cs-137 concentrations in Bq kg−1.

4 Discussion

4.1 Explanatory Power of GAM

The GAM results demonstrated about 50 % of Cs-137 deviance was explained by three parameters (Table 6). Model performance might be improved by including factors such as surface vegetation and subsurface organic matter (Claverie et al., 2019; Doerr and Münnich, 1989; Dumat et al., 1997, 2000; Fan et al., 2014; Korobova et al., 2016; Mabit et al., 2008; Staunton et al., 2002; Takenaka et al., 1998; Tatsuno et al., 2020). Topography itself did not explain the entire Cs-137 concentration patterns in the forest; however, the influence of topography on Cs-137 concentration patterns was measurable.

Elevation and slope are commonly used in geomorphic assessments of Cs-137 contamination (see references above). GAMs with only elevation or slope explained Cs-137 concentrations at less than 28.31 % in either unit (Bq kg−1 or Bq m−2). Despite its less significance for explaining Cs-137 deviance as a single parameter, elevation appeared seven times in the eight best-performing models (Table 6). The best models in all four categories (two by Cs-137 units; two by DEM resolutions) in Table 6 consisted of water content, soil dry bulk density, and elevation. The topographic parameters in the remaining four best-performing models varied. Elevation did explain Cs-137 deviance; however, its effects did not materialize by itself in the model result and needed to be combined with other supporting topographic features or soil properties. This indirect effect on Cs-137 deviance is true for other topographic parameters because the combination models returned higher explanation percentages for any single-parameter model. These results show that when the models account for both the characteristics of soils, as the Cs-137 retention medium, and topography, as the Cs-137 translocation driver, their prediction power increases.

At depths of 12–18 cm, the explanation percentages for the GAMs were above 75 % (Figs. 10 and 11), and Cs-137 displayed the largest COVs around the same depths (Fig. 5b and c). Because those depths were just above the average migration head depth, 21.7 cm, it is hypothesized that the downward movement of Cs-137 was affected by topography, and that these effects are retained as particles with adsorbed Cs-137 continued their interstitial movement downward.

The vertical profiles of the GAM explanation percentages were similar for soil properties and topographic parameters throughout the depth (Figs. 10 and 11). Pearson's correlation indices were calculated between the soil properties and topographic parameters. The only pairs that showed a correlation greater than 0.5 were water content at 12–16 cm depth vs. TWI, and soil dry bulk density at 2–4 cm vs. upslope distance. The explanation percentages by soil properties and topographic parameters were the result of a non-linear functional relationship between them.

The migration head depths in Fig. 7 did not show a distinct trend with elevation; however, the migration head depths for eight samples had already reached 30 cm deep. Table 7 compares the depths at which 90 % of the total Cs-137 concentration in soil samples was measured (not migration head depths) for the CNPP accident (Ivanov et al., 1997) and those reported here (values are not decay normalized). The present results show that Cs-137 concentrations had migrated deeper into the subsurface than in the region affected by the CNPP accident.

4.2 Model fit analysis

The two GAMs with the best overall performance, Model A with soil properties and elevation and Model B with only topographic parameters in Table 6, were explored to assess their predictive ability. While each GAM was run against all 14 depth layers, only the outputs for five depth layers are discussed here. These are 0–2, 4–6, and 8–10 cm to display model fits in the soils close to the surface where Cs-137 was concentrated, 14–16 cm to display the model fits in the mid-depth, and 20–22.5 cm to display the model fits at the average migration head depth. Figure 13 lists the fitted smoothing parameters of the model predictors for these two models.

Figure 13GAM smooth function outputs at five depth layers for the parameters in the two best-performing models. The x axis on each plot is the model predictor and its real value range. The y axis is the smoothing parameter, s(covariate, degree of freedom). Vertical dashed red lines indicate a point at which the smoothing parameter is the zero residual point.

For each plot in Fig. 13, the x axis is the model predictor and its value range, and the y axis is the smoothing parameter. The numbers on the y axis indicate the partial residuals from the fit. Thus, the y axis demonstrates an increase or decrease of the response values against its predicted mean. The gray area represents the 95 % confidence interval range, and the vertical dashed red line indicates the predictor value at which the smoothing parameter is zero mean.

The fitted curves for Model A (fitted smooth) show that water content and elevation were positively related to Cs-137 (Fig. 13a and c). Soil dry bulk density was negatively related to Cs-137 (Fig. 13b). Cs-137 increases as water content increases, although the increasing pattern is not linear in two depth intervals (Fig. 13a). The Cs-137 increases above the mean at around the water content 50 %–120 %. The samples in this study contained soils with mass water content above 100 %, and those soils had elevated Cs-137, when they were dried. Cs-137 decreased below the mean at a soil dry bulk density around 0.2–0.55 g cm−3 (Fig. 13b). GAMs returned higher explanation percentages when run against Cs-137 in Bq kg−1 unit than in Bq m−2 (Tables 4 to 6). These results indicated that soil dry bulk density, included in mass depth (Bq m−2) calculation, was not a significant factor in explaining Cs-137 concentrations. The model-fitted curves in Fig. 13b show the negative linear relationships between soil dry bulk density and Cs-137 predictions. Since soil dry bulk density increases with depth, Cs-137 is concentrated in low-density soils close to the surface. In this study site, soil dry bulk density had an influence on Cs-137 concentrations, which appeared when density was combined with other model parameters. Cs-137 and elevation showed positive linear correlations or U-shaped fitted curves depending on the depth (Fig. 13c). The bottom of the U shape was around the elevation of 610 m. One sample collected at the elevation of 613 m contained the highest core average soil dry bulk density, which also had the lowest core total Cs-137 concentration and the highest dry bulk density in the top 0–2 cm layer.

The fitted curves for Model B show that elevation was positively related to Cs-137 (Fig. 13d), and slope was negatively related to Cs-137 except in the 0–2 cm depth (Fig. 13e). Elevation was a common predictor in both models (Fig. 13c and d). The elevation's fitted curves in the two models showed that the same predictor's fitted smooths could vary depending on its relation to other predictors in the model. Slope showed a fitted curve near the surface (0–2 cm depth) with an increasing trend, but an opposite trend from other layers below (Fig. 13e). Three samples with the highest Cs-137 concentrations in the 0–2 cm depth were found at slopes steeper than 15. These samples pushed the trend upwards toward a steeper slope in the 0–2 cm depth range. Slope presented the fitted curves with multiple knots between 5–25 in the 14–16 cm depth. Yet, the fitted smooths did not move away from the mean, indicating that slope in this range contributed little to Cs-137 concentration. Upslope distance contributed to a Cs-137 increase toward the shortest upslope distance (at the ridge area) and the longest upslope distance (hill bottom area; Fig. 13f). The sample with the highest core total Cs-137 concentration was about 10 m below the ridge edge. Higher Cs-137 concentrations at the ridge area is counter-intuitive because Cs-137 was expected to move downslope with surface soil creep or runoff. The Cs-137 concentrations in five samples collected at the shortest upslope distances did not show any apparent relationships with water content, soil dry bulk density, or other topographic features.

On the basis of the model fits, predictions for the Cs-137 concentrations in this study site are as follows.

  1. Cs-137 concentrations generally increase with elevation, but the increase is not linear (Figs. 7 and 13c).

  2. Higher Cs-137 concentrations can be found at the highest elevation near the ridge or hill bottom (Fig. 13c and f).

  3. Higher Cs-137 concentrations will be observed at the locations with lower slopes, which could include flat areas at the ridge and near the hill bottom (Fig. 13e).

  4. Soils with higher water content might contain higher Cs-137 concentrations once they are dried (Fig. 13a).

  5. Higher Cs-137 concentrations can be found in lower soil dry bulk density close to the surface. Exceptions can be found at locations where water does not pond, even though soil dry bulk density is higher (Fig. 13b and c).

4.3 The effect of DEM resolution

The ratios between topographic parameter values extracted from the two DEMs are shown in Table 8. The last column of Table 8 indicates the DEM resolution where each topographic parameter appeared in the best-performing GAM models (Table 6). The elevation and slope values were consistent between the 1 and 10 m DEM, indicating that elevation did not change drastically over distances less than 10 m. Upslope distance had the largest difference between the DEMs. Plan curvature and TWI, which influence surface water ponding, appeared in the best-performing models with the 1 m DEM. Upslope distance performed better with the 10 m DEM, but the average upslope distance was 4 times greater with the 10 m DEM.

Table 8Topographic value difference ratio between the 1 and 10 m resolution DEMs (medians). The right column shows with which DEM, 1 or 10 m, each topographic parameter appeared in the best-performing model.

Download Print Version | Download XLSX

4.4 Study contribution

This study employed empirical approaches to derive predictive relationships among Cs-137 and soil property data and local topography. The statistical results show the power of detecting co-functionalities of topographic predictors on Cs-137 accumulations on the ground surface and with depth, which could enhance the understanding of the geomorphic mechanisms governing Cs-137 accumulations in a forested area. While model prediction errors might depend heavily on factors not included in this study, the current approach proposes a methodology to partition observed Cs-137 concentrations into explainable and unexplainable portions. The methodology and results presented here could contribute to ongoing research on understanding the distribution of radioactive contamination in complex terrains.

4.5 Study limitations

In this study, the Cs-137 measurement began in 2016, and the physical translocation of Cs-137 prior to 2016 was not considered as no soil monitoring was conducted. The only data available were the initial Cs-137 deposition estimates from the Japanese government and the US Department of Defense (MEXT, 2011), which was 1000–3000 kBq m−2 as of July 2011. The Cs-137 values are decay normalized and do not consider the physical movement of Cs-137 over the 3-year sampling period. While the effects of precipitation magnitude and frequency cannot be ignored, all soil sampling campaigns were conducted at approximately the same time of year (i.e., after the rainy season). All sampling campaigns avoided intense precipitation immediately prior to sample collections.

5 Conclusions

Following an earthquake near the northeast coast of Honshu, Japan, and a subsequent tsunami, the FDNPP suffered a catastrophic failure and a radioactive plume, including Cs-137, was released into the atmosphere. While much previous research focused on documenting Cs-137 concentrations in nearby environments, there has been no systematic and quantitative study that examined the effect of topographic indices and soil properties on the concentration and distribution of Cs-137 from this accident. To address this issue, 58 soil core samples were collected over a 3-year period from the ground surface to depths of 30 cm in Iitate village, a forested region located 35 km northwest of the FDNPP. The Cs-137 concentrations for these samples were assessed along with various topographic indices derived from digital elevation models and soil properties. Generalized additive models (GAMs) were then employed to quantitatively determine the singular and collective impact of topography, soil properties, and digital elevation model resolution on predicting the concentration of Cs-137 within this landscape. The primary conclusions of the study are as follows.

  1. In general, Cs-137 activities were highest near the ground surface and decreased exponentially with depth. The average depth of 90 % concentration was 6.1 cm in kBq kg−1 and 6.9 cm in kBq m−2, and the average migration depth for all samples was 21.7 cm.

  2. Higher concentrations of Cs-137 occurred at the highest elevations or near the bottom of hills.

  3. Elevation, slope, upslope distance, plan curvature, and TWI were derived for each sample location using 1 and 10 m DEMs. Resolution of the DEM produced marked differences in the upslope distance, plan curvature, and TWI indices determination, which then affected the Cs-137 predictions.

  4. Using two units for Cs-137 concentration, five topographic indices, and two soil properties, 63 GAMs were evaluated using one- to three-parameter combinations and two DEMs with different resolutions, which were assessed using several criteria. For the single-parameter models, none returned average deviance-explained percentages greater than 36 %. In contrast, the most effective model combinations for predicting Cs-137 deviance were “water content + bulk density + elevation” for the 1 m (53.56 %) and 10 m (54.03 %) DEMs, given the empirical data reported herein. The most effective topographic parameter combinations were “slope + plan curvature + TWI” when using the 1 m DEM (40.62 %) and “elevation + slope + upslope distance” when using the 10 m DEM (46.70 %). Concentrations of Cs-137 within this landscape were reverse-predicted using the best two GAMs (Model A: “water content + bulk density + elevation” and Model B: “elevation + slope + upslope distance”). For the data reported herein, Model A outperformed Model B.

While this study focused on a small forested region affected by the FDNPP accident, the results clearly show that Cs-137 concentrations in soils are strongly affected by landscape topography. This topographic effect should be given careful consideration in the use of anthropogenic radionuclides as environmental tracers and in the assessment of current and future environmental risks due to nuclear power plant accidents.

Appendix A: Equations

A1 R package mgcv GAM model settings

Model  gam(Cs-137 s(x1)+s(x2)+s(x3), data = data, method = “REML”, bs = “cr”, family = Gamma(link = log)), where “s” is a smooth function of covariates of the predictor, “method” specifies a fitting method, and “bs” specifies basis function (“cr”, cubic splines, in this model). “family” specifies a response variable distribution type. Here, “link” solves the problem of using linear models with non-normal data by linking the estimated fitted values to the linear predictor (Clark, 2013). This study used “Gamma(link = log)” to address complex, non-linear interactions among predictors.

A2 Cs-137 decay equations

The Cs-137 decay constant is calculated as follows (Eq. A1):

(A1) λ = ln 2 T 1 2 ,

where λ is the decay constant and T1/2 is the half-life. The half-life of Cs-137 is 30.17 years; thus, λ=0.023 yr−1. The radioactive decay formula is as follows (Eq. A2):

(A2) N 0 = N t e - λ t ,

where N0 is the original Cs-137 value, Nt is the Cs-137 value after time t, and λ is the decay constant (IAEA-TECDOC, 2003).

A3 Akaike information criterion (AIC)

(A3) AIC = - 2 l ( θ ^ ) + 2 p ( Wood , 2017 ) ,

where l(θ^) is log likelihood and p is the number of identifiable model parameters (usually, the dimension of θ).

A4 Generalized cross validation (GCV)

(A4) V g = n i = 1 n y i - f ^ i 2 [ n - tr ( A ) ] 2 ( Hastie et al . , 2009 ; Wood , 2017 ) ,

where Vg is the generalized cross-validation score, yi is the excluded data, f^ is the estimate from fitting to all the data, and tr(A) is the mean of the model matrix Aii. The matrix measures the degree that the ith datum influences the overall model fit.

A5 Cook's outliers

(A5) D i = β ^ ( - i ) - β ^ X X β ^ ( - i ) - β ^ ( p + 1 ) σ ^ 2 i = 1 , 2 , n ( Cook , 1977 ; Stevens , 1984 ) ,

where β^ is the vector of estimated regression coefficients with the ith data point deleted, p is the number of predictors, σ^2 is the residual variance of the full dataset.

Appendix B: Summary of all primary data collected

Table B1Summary of topographic indices for all core sample locations (58 core samples, 46 locations).

Download Print Version | Download XLSX

Table B2Summary of data collected for all core sample locations. Units include Cs-137: Bq kg−1, dry bulk density: g cm−3, mass water content: %, sand, silt, and clay content: %. Among 58 core samples, multi-year samples were decay normalized and averaged.

Download XLSX

Code availability

The programming applications used in this study were R and ArcMap and its add-on tools. Both are cited and credited in the body text. The authors used standard R codes and standard function tools on ArcMap. R code used for GAM regression is presented in the Appendix. No proprietary codes or applications that readers cannot access were used.

Data availability

The data used in this paper are presented in Tables B1 and B2.

Author contributions

MY conceived and designed the project. MY and TN conducted field data collection and laboratory analysis. MY conducted data analysis and graphing. TN, JA, and SJB provided critical feedback on data analysis. MY drafted the manuscript. TB advised on data concepts. All authors provided critical feedback to finalize the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Field sampling, soil sample lab tests, and radioactivity measurements were accomplished with the support of the Laboratory of Soil Physics and Soil Hydrology, Graduate School of Agricultural and Life Sciences, the University of Tokyo, Japan. Forestry and Forest Products Research Institute, Japan, provided the 1 m resolution DEM. Kinichi Okubo, a farmer in Iitate village, kindly provided his land for this research.

Financial support

This research has been supported by the NSF East Asia and Pacific Summer Institutes (EAPSI) 2016 (grant no. 1614049), the NSF Doctoral Dissertation Research Improvement Award – Geography and Spatial Sciences Program (DDRI-GSS) (grant no. 1819727), and the College of Arts and Sciences Dissertation Fellowship of the University at Buffalo.

Review statement

This paper was edited by Edward Tipper and reviewed by Elena Korobova and two anonymous referees.


Amaral, E., Amundsen, I., Barišić, D., Booth, P., Clark, D., Ditmars, J., Dlouhy, Z., Drury, N., Gehrche, K., and Gnugnoli, G.: Characterization of radioactively contaminated sites for remediation purposes, International Atomic Energy Agency, IAEA-Tecdoc, 1011-4289, International Atomic Energy Agency, Vienna, Austria, 1998. 

Ayabe, Y., Hijii, N., and Takenaka, C.: Effects of local-scale decontamination in a secondary forest contaminated after the Fukushima nuclear power plant accident, Environ. Pollut., 228, 344–353, 2017. 

Bennett, S. J., Rhoton, F. E., and Dunbar, J. A.: Texture, spatial distribution, and rate of reservoir sedimentation within a highly erosive, cultivated watershed: Grenada Lake, Mississippi, Water Resour. Res., 41, W01005,, 2005. 

Bivand, R. S., Pebesma, E. J., Gómez-Rubio, V., and Pebesma, E. J.: Applied spatial data analysis with R, Springer, New York, USA, 2008. 

Burr, T. and Hamada, M.: Radio-Isotope Identification Algorithms for NaI γ Spectra, Algorithms, 2, 339–360,, 2009. 

Catani, F., Segoni, S., and Falorni, G.: An empirical geomorphology-based approach to the spatial prediction of soil thickness at catchment scale, Water Resour. Res., 46, W05508,, 2010. 

Chen, Z.-S., Hsieh, C.-F., Jiang, F.-Y., Hsieh, T.-H., and Sun, I.-F.: Relations of soil properties to topography and vegetation in a subtropical rain forest in southern Taiwan, Plant Ecol., 132, 229–241,, 1997. 

Claessens, L., Heuvelink, G., Schoorl, J., and Veldkamp, A.: DEM resolution effects on shallow landslide hazard and soil redistribution modelling, Earth Surf. Proc. Land., 30, 461–477, 2005. 

Clark, M.: Generalized additive models: getting started with additive models in R, Center for Social Research, University of Notre Dame, Notre Dame, 35, 2013. 

Claverie, M., Garcia, J., Prevost, T., Brendlé, J., and Limousy, L.: Inorganic and Hybrid (Organic–Inorganic) Lamellar Materials for Heavy Metals and Radionuclides Capture in Energy Wastes Management – A Review, Materials, 12, 1399,, 2019. 

Climate-Data.Org: Prypiat Climate, available at:, last access: 25 November 2019. 

Conrad, O., Bechtel, B., Bock, M., Dietrich, H., Fischer, E., Gerlitz, L., Wehberg, J., Wichmann, V., and Bohner, J.: System for Automated Geoscientific Analyses (SAGA) v. 2.1.4, Geosci. Model Dev., 8, 1991–2007,, 2015. 

Cook, R. D.: Detection of influential observation in linear regression, Technometrics, 19, 15–18,, 1977. 

Doerr, H. and Münnich, K.: Downward movement of soil organic matter and its influence on trace-element transport (210Pb, 137Cs) in the soil, Radiocarbon, 31, 655–663,, 1989. 

Dumat, C., Cheshire, M., Fraser, A., Shand, C., and Staunton, S.: The effect of removal of soil organic matter and iron on the adsorption of radiocaesium, Eur. J. Soil Sci., 48, 675–683,, 1997. 

Dumat, C., Quiquampoix, H., and Staunton, S.: Adsorption of cesium by synthetic clay- organic matter complexes: effect of the nature of organic polymers, Environ. Sci. Technol., 34, 2985–2989,, 2000. 

Eguchi, S., Aoki, K., and Kohyama, K.: Development of agricultural soil-profile physical properties database, Japan: SolphyJ, in: Proceedings of the ASA-CSSA-SSSA International Annual Meetings, San Antonio, TX, USA, 16–19, 2011. 

Endo, S., Kimura, S., Takatsuji, T., Nanasawa, K., Imanaka, T., and Shizuma, K.: Measurement of soil contamination by radionuclides due to the Fukushima Dai-ichi Nuclear Power Plant accident and associated estimated cumulative external dose estimation, J. Environ. Radioact., 111, 18–27, 2012. 

Endo, S., Kajimoto, T. and Shizuma, K.: Paddy-field contamination with 134 Cs and 137 Cs due to Fukushima Dai-ichi Nuclear Power Plant accident and soil-to-rice transfer coefficients, J. Environ. Radioact., 116, 59–64, 2013. 

EPA: Radiological Laboratory Sample Analysis Guide for Incident Response – Radionuclides in Soil, US EPA, Montgomery, AL, USA, 2012. 

EPA: PAG Manual: Protective Action Guides and Planning Guidance for Radiological Incidents, US EPA, Washington, DC, USA, 2017. 

Fan, Q., Tanaka, M., Tanaka, K., Sakaguchi, A., and Takahashi, Y.: An EXAFS study on the effects of natural organic matter and the expandability of clay minerals on cesium adsorption and mobility, Geochim. Cosmochim. Ac., 135, 49–65,, 2014. 

Forest Management Center: Topography, geology, and climate of Fukushima Prefecture, available at:, last access: 1 July 2017. 

Friedman, J. H. and Stuetzle, W.: Projection pursuit regression, J. Am. Stat. Assoc., 76, 817–823, 1981. 

Fujii, K., Ikeda, S., Akama, A., Komatsu, M., Takahashi, M., and Kaneko, S.: Vertical migration of radiocesium and clay mineral composition in five forest soils contaminated by the Fukushima nuclear accident, Soil Sci. Plant Nutr., 60, 751–764,, 2014. 

Fujiwara, T., Saito, T., Muroya, Y., Sawahata, H., Yamashita, Y., Nagasaki, S., Okamoto, K., Takahashi, H., Uesaka, M., and Katsumura, Y.: Isotopic ratio and vertical distribution of radionuclides in soil affected by the accident of Fukushima Dai-ichi nuclear power plants, J. Environ. Radioact., 113, 37–44, 2012. 

Fukushima Prefecture: Fukushima Forest and Culture Exhibition Pictorial Record, available at:, last access: 3 March 2021. 

Gallant, J. C. and Wilson, J. P.: TAPES-G: a grid-based terrain analysis program for the environmental sciences, Comput. Geosci., 22, 713–722,, 1996. 

Gessler, P. E., Moore, I., McKenzie, N., and Ryan, P.: Soil-landscape modelling and spatial prediction of soil attributes, Int. J. Geogr. Inform. Syst., 9, 421–432,, 1995. 

Giannakopoulou, F., Haidouti, C., Chronopoulou, A., and Gasparatos, D.: Sorption behavior of cesium on various soils under different pH levels, J. Hazard. Mater., 149, 553–556, 2007. 

Griffiths, R. P., Madritch, M. D., and Swanson, A. K.: The effects of topography on forest soil characteristics in the Oregon Cascade Mountains (USA): Implications for the effects of climate change on soil properties, Forest Ecol. Manage., 257, 1–7,, 2009. 

Hastie, T., Tibshirani, R., and Friedman, J.: The elements of statistical learning: data mining, inference, and prediction, Springer Science & Business Media, New York, USA, 2009. 

Hastie, T. J. and Tibshirani, R. J.: Generalized additive models, Chapman & Hall/CRC, Washington, DC, USA, 1990. 

Heimsath, A., Dietrich, W., Nishiizumi, K., and Finkel, R.: The soil production function and landscape equilibrium, Nature, 388, 6640,, 1997. 

Heimsath, A. M., Dietrich, W. E., Nishiizumi, K., and Finkel, R. C.: Cosmogenic nuclides, topography, and the spatial variation of soil depth, Geomorphology, 27, 151–172,, 1999. 

Hengl, T. and Reuter, H. I.: Geomorphometry: concepts, software, applications, Newnes, Oxford, UK, 2008. 

Hoover, M. and Hursh, C.: Influence of topography and soil-depth on runoff from forest land, Eos Trans. Am. Geophys. Union, 24, 693–698,, 1943. 

IAEA: Fukushima Daiichi Accident Technical Volume 4/5 Radiological Consequences, Vienna, Austria, 2015. 

IAEA-TECDOC: Manual for reactor produced radioisotopes, International Atomic Energy Agency, Vienna, Austria, 2003. 

IAEA and INES: The International Nuclear and Radiological Event Scale User's Manual, International Atomic Energy Agency, Vienna, 2013. 

Ivanov, Y. A., Lewyckyj, N., Levchuk, S., Prister, B., Firsakova, S., Arkhipov, N., Arkhipov, A., Kruglov, S., Alexakhin, R., and Sandalls, J.: Migration of 137Cs and 90Sr from Chernobyl fallout in Ukrainian, Belarussian and Russian soils, J. Environ. Radioact., 35, 1–21,, 1997. 

Iwagami, S., Tsujimura, M., Onda, Y., Nishino, M., Konuma, R., Abe, Y., Hada, M., Pun, I., Sakaguchi, A., and Kondo, H.: Temporal changes in dissolved 137 Cs concentrations in groundwater and stream water in Fukushima after the Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 166, 458–465,, 2015. 

Japan Meteorological Agency: Data and References, available at:, last access: 2 September 2019. 

Kato, H., Onda, Y., and Teramage, M.: Depth distribution of 137 Cs, 134 Cs, and 131 I in soil profile after Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 111, 59–64,, 2012. 

Kim, S. and Lee, H.: A digital elevation analysis: a spatially distributed flow apportioning algorithm, Hydrol. Process., 18, 1777–1794,, 2004. 

Koarashi, J., Atarashi-Andoh, M., Matsunaga, T., Sato, T., Nagao, S., and Nagai, H.: Factors affecting vertical distribution of Fukushima accident-derived radiocesium in soil under different land-use conditions, Sci. Total Environ., 431, 392–401,, 2012. 

Komamura, M., Tsumura, A., Yamaguchi, N., Fujiwara, E., Konokata, N., and Kodaira, K.: Long-term monitoring and variability analysis of 90Sr and 137Cs in rice, wheat, and soils in Japan Agriculture, Forestry and Fisheries Research Council Secretariat, Ibaragi, Japan, 24, 1–21, 2006. 

Komissarov, M. and Ogura, S.: Distribution and migration of radiocesium in sloping landscapes three years after the Fukushima-1 nuclear accident, Euras. Soil Sci., 50, 861–871,, 2017. 

Korobova, E. M., Linnik, V. G., and Brown, J.: Distribution of artificial radioisotopes in granulometric and organic fractions of alluvial soils downstream from the Krasnoyarsk Mining and Chemical Combine (KMCC), Russia, J. Soils Sediments, 16, 1279–1287,, 2016. 

Lepage, H., Evrard, O., Onda, Y., Lefèvre, I., Laceby, J. P., and Ayrault, S.: Depth distribution of cesium-137 in paddy fields across the Fukushima pollution plume in 2013, J. Environ. Radioact., 147, 157–164, 2015. 

Linnik, V. G., Saveliev, A. A., and Sokolov, A. V.: Transformation of the Chernobyl 137 Cs Contamination Patterns at the Microlandscape Level as an Indicator of Stochastic Landscape Organization, Landscape Patterns in a Range of Spatio-Temporal Scales, Springer, Cham, Switzerland, 77–89, 2020. 

Loughran, R., Campbell, B., and Walling, D.: Soil erosion and sedimentation indicated by caesium 137: Jackmoor Brook catchment, Devon, England, Catena, 14, 201–212,, 1987. 

Lowrance, R., McIntyre, S., and Lance, C.: Erosion and deposition in a field/forest system estimated using cesium-137 activity, J. Soil Water Conserv., 43, 195–199, 1988. 

Mabit, L., Bernard, C., and Laverdière, M. R.: Assessment of erosion in the Boyer River watershed (Canada) using a GIS oriented sampling strategy and 137Cs measurements, Catena, 71, 242–249,, 2007. 

Mabit, L., Bernard, C., Makhlouf, M., and Laverdière, M.: Spatial variability of erosion and soil organic matter content estimated from 137Cs measurements and geostatistics, Geoderma, 145, 245–251,, 2008. 

Maekawa, A., Momoshima, N., Sugihara, S., Ohzawa, R., and Nakama, A.: Analysis of Cs and Cs distribution in soil of Fukushima prefecture and their specific adsorption on clay minerals, J. Radioanal. Nucl. Chem., 303, 1485–1489, 2015. 

Mahaffey, J.: Atomic Accidents: A History of Nuclear Meltdowns and Disasters: from the Ozark Mountains to Fukushima, Open Road Media, New York, USA, 2014. 

Martin, Y.: Modelling hillslope evolution: linear and nonlinear transport relations, Geomorphology, 34, 1–21,, 2000. 

Martinez, C., Hancock, G., Kalma, J., Wells, T., and Boland, L.: An assessment of digital elevation models and their ability to capture geomorphic and hydrologic properties at the catchment scale, Int. J. Remote Sens., 31, 6239–6257, 2010. 

Martin-Garin, A., Van Meir, N., Simonucci, C., Kashparov, V., and Bugai, D.: Quantitative assessment of radionuclide migration from near-surface radioactive waste burial sites: the waste dumps in the Chernobyl exclusion zone as an example, in: Radionuclide Behaviour in the Natural Environment, Elsevier, Cambridge, UK, 570–600, 2012. 

Martz, L. and De Jong, E.: Using Cesium-137 to assess the variability of net soil erosion and its association with topography in a Canadian prairie landscape, Catena, 14, 439–451,, 1987. 

Matsuda, N., Mikami, S., Shimoura, S., Takahashi, J., Nakano, M., Shimada, K., Uno, K., Hagiwara, S., and Saito, K.: Depth profiles of radioactive cesium in soil using a scraper plate over a wide area surrounding the Fukushima Dai-ichi Nuclear Power Plant, Japan, J. Environ. Radioact., 139, 427–434, 2015. 

Matsunaga, T., Koarashi, J., Atarashi-Andoh, M., Nagao, S., Sato, T., and Nagai, H.: Comparison of the vertical distributions of Fukushima nuclear accident radiocesium in soil before and after the first rainy season, with physicochemical and mineralogical interpretations, Sci. Total Environ., 447, 301–314,, 2013. 

MEXT: The 48th environmental radioactivity research result abstract papers, Ministry of Education, Culture, Sports, Science and Technology, Tokyo, Japan, 2006. 

MEXT: Regarding the production of radiological dosage distribution map by MEXT (Mapping radiocesium in soils), available at: (last access: 1 October 2017), 2011. 

Miller, B. A., Koszinski, S., Wehrhan, M., and Sommer, M.: Impact of multi-scale predictor selection for modeling soil properties, Geoderma, 239, 97–106, 2015. 

Miyahara, K. A., Kohara, T., Yusa, Y., and Sasaki, N.: Effect of bulk density on diffusion for cesium in compacted sodium bentonite, Radiochim. Acta, 52, 293–298,, 1991. 

Momm, H. G., Bingner, R. L., Wells, R. R., and Wilcox, D.: AGNPS GIS-based tool for watershed-scale identification and mapping of cropland potential ephemeral gullies, App. Eng. Agricult., 28, 17–29,, 2012. 

Moore, I. D., Gessler, P., Nielsen, G., and Peterson, G.: Soil attribute prediction using terrain analysis, Soil Sci. Soc. Am. J., 57, 443–452,, 1993. 

Murota, K., Saito, T., and Tanaka, S.: Desorption kinetics of cesium from Fukushima soils, J. Environ. Radioact., 153, 134–140,, 2016. 

Nagao, S.: Radionuclides released from nuclear accidents: Distribution and dynamics in soil, in: Environmental Remediation Technologies for Metal-Contaminated Soils, Springer, Tokyo, Japan, 43–65, 2016. 

Nakanishi, T., Kobayashi, N., and Tanoi, K.: Radioactive cesium deposition on rice, wheat, peach tree and soil after nuclear accident in Fukushima, J. Radioanal. Nucl. Chem., 296, 985–989, 2013. 

Nakanishi, T., Matsunaga, T., Koarashi, J., and Atarashi-Andoh, M.: 137Cs vertical migration in a deciduous forest soil following the Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 128, 9–14, 2014. 

Nakao, A., Thiry, Y., Funakawa, S., and Kosaki, T.: Characterization of the frayed edge site of micaceous minerals in soil clays influenced by different pedogenetic conditions in Japan and northern Thailand, Soil Sci. Plant Nutr., 54, 479–489,, 2008. 

Nobori, T., Tanoi, K., and Nakanishi, T.: Method of radiocesium determination in soil and crops using NaI (T1) scintillation counter attached with auto sampler, J. Soil Sci. Plant Nutr., 84, 182–186, 2013. 

Ohno, T., Muramatsu, Y., Miura, Y., Oda, K., Inagawa, N., Ogawa, H., Yamazaki, A., Toyama, C., and Sato, M.: Depth profiles of radioactive cesium and iodine released from the Fukushima Daiichi nuclear power plant in different agricultural fields and forests, Geochem. J., 46, 287–295, 2012. 

Ohnuki, T. and Kozai, N.: Adsorption behavior of radioactive cesium by non-mica minerals, J. Nucl. Sci. Technol., 50, 369–375,, 2013. 

Onda, Y., Kato, H., Hoshi, M., Takahashi, Y., and Nguyen, M.-L.: Soil sampling and analytical strategies for mapping fallout in nuclear emergencies based on the Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 139, 300–307,, 2015. 

Osawa, K., Nonaka, Y., Nishimura, T., Tanoi, K., Matsui, H., Mizogichi, M., and Tatsuno, T.: Quantification of dissolved and particulate radiocesium fluxes in two rivers draining the main radioactive pollution plume in Fukushima, Japan (2013–2016), Anthropocene, 22, 40–50,, 2018. 

Park, S.-M., Alessi, D. S., and Baek, K.: Selective adsorption and irreversible fixation behavior of cesium onto 2:1 layered clay mineral: A mini review, J. Hazard. Mater., 369, 569–576,, 2019. 

Pennock, D., Jong, E. D., and Lemmen, D.: Cesium-137-measured erosion rates for soils of five parent-material groups in southwestern Saskatchewan, Can. J. Soil Sci., 75, 205–210,, 1995. 

Quine, T. A., Govers, G., Walling, D. E., Zhang, X., Desmet, P. J., Zhang, Y., and Vandaele, K.: Erosion processes and landform evolution on agricultural land – new perspectives from caesium-137 measurements and topographic-based erosion modelling, Earth Surf. Proc. Land., 22, 799–816,<799::aid-esp765>;2-r, 1997. 

Quinn, P., Beven, K., Chevallier, P., and Planchon, O.: The prediction of hillslope flow paths for distributed hydrological modelling using digital terrain models, Hydrol. Process., 5, 59–79,, 1991. 

Quinn, P., Beven, K., and Lamb, R.: The in (a/tan/β) index: How to calculate it and how to use it within the topmodel framework, Hydrol. Process., 9, 161–182,, 1995. 

R Core Team: R Foundation for Statistical Computing, 2014, R: A language and environment for statistical computing, 2013, Vienna, Austria, 2015. 

Ritchie, J. and Ritchie, C.: 137 Cs use in erosion and sediment deposition studies: Promises and problems, International Atomic Energy Agency, Vienna, Austria, 1995. 

Ritchie, J. C. and McHenry, J. R.: Application of radioactive fallout cesium-137 for measuring soil erosion and sediment accumulation rates and patterns: a review, J. Environ. Qual., 19, 215–233,, 1990. 

Roering, J., Kirchner, J., and Dietrich, W.: Hillslope evolution by nonlinear, slope-dependent transport: Steady state morphology and equilibrium adjustment timescales, J. Geophys. Res.-Solid, 106, 16499–16513,, 2001. 

Roering, J. J.: How well can hillslope evolution models “explain” topography? Simulating soil transport and production with high-resolution topographic data, Geol. Soc. Am. Bull., 120, 1248–1262,, 2008. 

Roering, J. J., Kirchner, J. W., and Dietrich, W. E.: Evidence for nonlinear, diffusive sediment transport on hillslopes and implications for landscape morphology, Water Resour. Res., 35, 853–870,, 1999. 

Rosén, K., Öborn, I., and Lönsjö, H.: Migration of radiocaesium in Swedish soil profiles after the Chernobyl accident, 1987–1995, J. Environ. Radioact., 46, 45–66,, 1999. 

Rossi, G., Ferrarini, A., Dowgiallo, G., Carton, A., Gentili, R., and Tomaselli, M.: Detecting complex relations among vegetation, soil and geomorphology. An in-depth method applied to a case study in the Apennines (Italy), Ecol. Complex., 17, 87–98,, 2014. 

Saito, T., Makino, H., and Tanaka, S.: Geochemical and grain-size distribution of radioactive and stable cesium in Fukushima soils: implications for their long-term behavior, J. Environ. Radioac., 138, 11–18, 2014. 

Saito, K., Tanihata, I., Fujiwara, M., Saito, T., Shimoura, S., Otsuka, T., Onda, Y., Hoshi, M., Ikeuchi, Y., and Takahashi, F.: Detailed deposition density maps constructed by large-scale soil sampling for gamma-ray emitting radioactive nuclides from the Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 139, 308–319, 2015. 

Sakai, M., Gomi, T., Nunokawa, M., Wakahara, T., and Onda, Y.: Soil removal as a decontamination practice and radiocesium accumulation in tadpoles in rice paddies at Fukushima, Environ. Pollut., 187, 112–115, 2014. 

Sakuma, K., Tsuji, H., Hayashi, S., Funaki, H., Malins, A., Yoshimura, K., Kurikami, H., Kitamura, A., Iijima, K., and Hosomi, M.: Applicability of Kd for modelling dissolved 137Cs concentrations in Fukushima river water: Case study of the upstream Ota River, J. Environ. Radioact., 184, 53–62,, 2018. 

Schimmack, W., Bunzl, K., and Zelles, L.: Initial rates of migration of radionuclides from the Chernobyl fallout in undisturbed soils, Geoderma, 44, 211–218,, 1989. 

Schimmack, W., Bunzl, K., Dietl, F., and Klotz, D.: Infiltration of radionuclides with low mobility (137Cs and 60Co) into a forest soil. Effect of the irrigation intensity, J. Environ. Radioact., 24, 53–63,, 1994. 

Shiozawa, S., Tanoi, K., Nemoto, K., Yoshida, S., Nishida, K., Hashimoto, K., Sakurai, K., Nakanishi, T. M., Nihei, N., and Ono, Y.: Vertical concentration profiles of radioactive caesium and convective velocity in soil in paddy field in Fukushima, Radioisotopes, 60, 323–328,, 2011. 

Sørensen, R. and Seibert, J.: Effects of DEM resolution on the calculation of topographical indices: TWI and its components, J. Hydrol., 347, 79–89, 2007. 

Staunton, S., Dumat, C., and Zsolnay, A.: Possible role of organic matter in radiocaesium adsorption in soils, J. Environ. Radioact., 58, 163–173,, 2002. 

Stevens, J. P.: Outliers and influential data points in regression analysis, Psycholog. Bull., 95, 334–344,, 1984. 

Stinnett, J. and Sullivan, C. J.: An automated isotope identification algorithm using Bayesian statistics, in: 2013 IEEE Nuclear Science Symposium and Medical Imaging Conference (2013 NSS/MIC), 27 October–2 November 2013, Seoul, South Korea, 1–7,, 2013. 

Tagami, K., Uchida, S., Uchihori, Y., Ishii, N., Kitamura, H., and Shirakawa, Y.: Specific activity and activity ratios of radionuclides in soil collected about 20 km from the Fukushima Daiichi Nuclear Power Plant: Radionuclide release to the south and southwest, Sci. Total Environ. 409, 4885–4888, 2011. 

Taira, Y., Hayashida, N., Yamashita, S., Kudo, T., Matsuda, N., Takahashi, J., Gutevitc, A., Kazlovsky, A., and Takamura, N.: Environmental contamination and external radiation dose rates from radionuclides released from the Fukushima Nuclear Power Plant, Radiat. Protect. Dosimet., 151, 537–545, 2012. 

Takahashi, J., Tamura, K., Suda, T., Matsumura, R., and Onda, Y.: Vertical distribution and temporal changes of 137Cs in soil profiles under various land uses after the Fukushima Dai-ichi Nuclear Power Plant accident, J. Environ. Radioact., 139, 351–361,, 2015. 

Takata, Y., Kohyama, K., Obara, H., Maejima, Y., Ishitsuka, N., Saito, T., and Taniyama, I.: Spatial prediction of radioactive Cs concentration in agricultural soil in eastern Japan, Soil Sci. Plant Nutr., 60, 393–403, 2014. 

Takenaka, C., Onda, Y., and Hamajima, Y.: Distribution of cesium-137 in Japanese forest soils: Correlation with the contents of organic carbon, Sci. Total Environ., 222, 193–199,, 1998. 

Tanaka, K., Takahashi, Y., Sakaguchi, A., Umeo, M., Hayakawa, S., Tanida, H., Saito, T., and Kanai, Y.: Vertical profiles of iodine-131 and cesium-137 in soils in Fukushima prefecture related to the Fukushima Daiichi Nuclear Power Station accident, Geochem. J., 46, 73–76,, 2012. 

Tanaka, K., Iwatani, H., Takahashi, Y., Sakaguchi, A., Yoshimura, K., and Onda, Y.: Investigation of spatial distribution of radiocesium in a paddy field as a potential sink, PloS One, 8, e80794,, 2013. 

Tanoi, K., Nobori, T., Shiomi, S., Saito, T., Kobayashi, N. I., Leonhardt, N., and Nakanishi, T. M.: Chapter 3 Cesium Translocation in Rice, in: Agricultural Implications of the Fukushima Nuclear Accident (III), Springer, Singapore, 15–25, 2019. 

Tarboton, D. G.: A new method for the determination of flow directions and upslope areas in grid digital elevation models, Water Resour. Res., 33, 309–319,, 1997. 

Tatsuno, T., Hamamoto, S., Nihei, N., and Nishimura, T.: Effects of the dissolved organic matter on Cs transport in the weathered granite soil, J. Environ. Manage., 254, 109785,, 2020. 

Teramage, M. T., Onda, Y., Patin, J., Kato, H., Gomi, T., and Nam, S.: Vertical distribution of radiocesium in coniferous forest soil after the Fukushima nuclear power plant accident, J. Environ. Radioact., 137, 37–45,, 2014. 

Tesfa, T. K., Tarboton, D. G., Chandler, D. G., and McNamara, J. P.: Modeling soil depth from topographic and land cover attributes, Water Resour. Res., 45, W10438,, 2009. 

Tsoulfanidis, N.: Nuclear energy: selected entries from the encyclopedia of sustainability science and technology, Springer Science & Business Media, Goteborg, Sweden, 2012. 

Tsuji, H., Nishikiori, T., Yasutaka, T., Watanabe, M., Ito, S., and Hayashi, S.: Behavior of dissolved radiocesium in river water in a forested watershed in Fukushima Prefecture, J. Geophys. Res.-Biogeo., 121, 2588–2599,, 2016. 

Tsukada, H., Hasegawa, H., Hisamatsu, S. I., and Yamasaki, S. I.: Transfer of 137 Cs and stable Cs from paddy soil to polished rice in Aomori, Japan, J. Environ. Radioact., 59, 351–363,, 2002. 

US EPA: Radionuclide Basics: Cesium-137, Radiation Protection, available at:, last access: 22 August 2017. 

Wakai, N., Yamashita, J., Enomoto, T., Tadashi, H., Ono, T., and Maeda, M.: Factors affecting 137Cs concentration in wild plants and soils of different land use in Iitate village after the Fukushima nuclear power plant accident, Radiat. Safe. Manage., 18, 1–8,, 2019. 

Wallbrink, P., Roddy, B., and Olley, J.: A tracer budget quantifying soil redistribution on hillslopes after forest harvesting, Catena, 47, 179–201,, 2002. 

Walling, D. and He, Q.: Models for converting 137Cs measurements to estimates of soil redistribution rates on cultivated and uncultivated soils (including software for model implementation), Report to IAEA, University of Exeter, Exeter, UK, 315–341, 1997. 

Walling, D. and Quine, T.: The use of caesium-137 measurements in soil erosion surveys, in: Erosion and sediment transport monitoring programmes in river basins, 24–28 August 1992, Oslo, Norway, 210, 143–152, 1992. 

Walling, D., He, Q., and Quine, T.: Use of caesium-137 and lead-210 as tracers in soil erosion investigations, IAHS Publ., 229, 163–172, 1995. 

Walling, D., Zhang, Y., and He, Q.: Models for converting measurements of environmental radionuclide inventories (137Cs, Excess 210Pb, and 7Be) to estimates of soil erosion and deposition rates (including software for model implementation), Department of Geography, University of Exeter, Exeter, UK, 2007. 

Wood, S. N.: Generalized additive models: an introduction with R, Chapman and Hall/CRC, Boca Raton, Florida, 2017. 

Wood, S. N., Pya, N., and Säfken, B.: Smoothing parameter and model selection for general smooth models, J. Am. Stat. Assoc., 111, 1548–1563, 2016. 

Wrixon, A. D., Barraclough, I., and Clark, M. J.: Radiation, People And The Environment, International Atomic Energy Agency, Vienna, Austria, 2004. 

Xinbao, Z., Higgitt, D., and Walling, D.: A preliminary assessment of the potential for using caesium-137 to estimate rates of soil erosion in the Loess Plateau of China, Hydrolog. Sci. J., 35, 243–252,, 1990. 

Yamaguchi, N., Takada, Y., Hayashi, K., Ishikawa, S., Kuramata, M., Eguchi, S., Yoshikara, S., Sakaguchi, A., Asada, K., and Rota, W.: Behavior of radiocaesium in soil-plant systems and its controlling factor, Report 31, National Institute for Agro-Environmental Sciences, Ibaragi, Japan, 75–129, 2012. 

Yamamoto, M., Takada, T., Nagao, S., Koike, T., Shimada, K., Hoshi, M., Zhumadilov, K., Shima, T., Fukuoka, M., and Imanaka, T.: An early survey of the radioactive contamination of soil due to the Fukushima Dai-ichi Nuclear Power Plant accident, with emphasis on plutonium analysis, Geochem. J., 46, 4, 341–353, 2012. 

Yang, B., Onda, Y., Wakiyama, Y., Yoshimura, K., Sekimoto, H., and Ha, Y.: Temporal changes of radiocesium in irrigated paddy fields and its accumulation in rice plants in Fukushima, Environ. Pollut., 208, 562–570, 2016.  

Yang, C., Peterson, C., Shropshire, G., and Otawa, T.: Spatial variability of field topography and wheat yield in the palouse region of the Pacific Northwest, T. ASAE, 41, 17–27,, 1998. 

Yoshikawa, N., Obara, H., Ogasa, M., Miyazu, S., Harada, N., and Nonaka, M.: 137 Cs in irrigation water and its effect on paddy fields in Japan after the Fukushima nuclear accident, Sci. Total Environ., 481, 252–259, 2014. 

Zaslavsky, D. and Sinai, G.: Surface hydrology: I – explanation of phenomena, J. Hydraul. Div., 107, 1–16, 1981. 

Zevenbergen, L. W. and Thorne, C. R.: Quantitative analysis of land surface topography, Earth Surf. Proc. Land., 12, 47–56, 1987. 

Zheng, J., Tagami, K., Watanabe, Y., Uchida, S., Aono, T., Ishii, N., Yoshida, S., Kubota, Y., Fuma, S., and Ihara, S.: Isotopic evidence of plutonium release into the environment from the Fukushima DNPP accident, Scient. Rep., 2, 1–8, 2012. 

Short summary
Topographic effects on radioactive contamination in a forested area were quantitatively examined using soil core samples collected in a village in Fukushima, Japan. The results confirmed that local topography influences the contamination patterns in soils, and its effects vary depending on the combinations of the topographic parameters. This finding suggests that topographic characteristics should be considered carefully in future environmental radioactive risk assessments.